史晶瑩 陳燕才 蔡曉明 高翔 徐恒元 梁昌晶

1大慶油田設(shè)計(jì)院有限公司

2河北華北石油工程建設(shè)有限公司

3渤海鉆探第四鉆井工程分公司

4河北華北石油港華勘察規(guī)劃設(shè)計(jì)有限公司

隨著油氣田的深入開發(fā)，氣藏的出砂量越來(lái)越大，高壓高流速的天然氣攜帶固體顆粒以固定的沖擊角度不斷磨損管道內(nèi)壁、閥門和彎頭等部件，長(zhǎng)時(shí)間的沖刷腐蝕會(huì)對(duì)管道造成侵蝕或刻蝕破壞，其中彎頭因涉及到二次流問(wèn)題，沖蝕磨損最為嚴(yán)重[1]。因此，研究不同流動(dòng)參數(shù)對(duì)彎管的沖蝕規(guī)律及固體顆粒的運(yùn)動(dòng)軌跡對(duì)于預(yù)防輸氣管道攜砂沖蝕具有重要意義[2-3]。

國(guó)內(nèi)外諸多學(xué)者采用旋轉(zhuǎn)圓柱電極、管流試驗(yàn)裝置以及噴射試驗(yàn)裝置進(jìn)行了沖蝕磨損試驗(yàn)，并根據(jù)試驗(yàn)結(jié)果總結(jié)了大量的沖蝕模型[4-8]，但試驗(yàn)的可重復(fù)性不高，且試驗(yàn)設(shè)計(jì)具有盲目性，因此目前多采用Fluent、Comsol 等流體力學(xué)軟件進(jìn)行沖蝕數(shù)值模擬。張孟昀等[9]通過(guò)加入DPM離散相模型，研究了相同條件下盲管與彎管的流場(chǎng)分布與沖蝕情況，發(fā)現(xiàn)彎管的耐蝕性遠(yuǎn)小于盲管；李滄等[10]采用Fluent 軟件對(duì)輸氣管道氣-固兩相流進(jìn)行了模擬計(jì)算，發(fā)現(xiàn)最大沖蝕速率隨入口流速和固體顆粒質(zhì)量流量呈正比關(guān)系；曹學(xué)文等[11]分別研究了砂粒參數(shù)、流動(dòng)參數(shù)對(duì)沖蝕速率的影響，發(fā)現(xiàn)流速與沖蝕速率呈正相關(guān)，出口壓力變化對(duì)沖蝕范圍及沖蝕速率影響不大；YAP Y J[12]采用計(jì)算流體力學(xué)軟件研究了顆粒直徑、管徑比對(duì)彎管的沖蝕影響，得到管徑比的增加會(huì)減小沖蝕區(qū)域的結(jié)論。以上研究對(duì)彎管沖蝕規(guī)律的研究具有重要意義，但均基于特定工況且采用單一變量控制法，對(duì)不同流動(dòng)參數(shù)之間的變化規(guī)律認(rèn)識(shí)不足，且未見對(duì)最大沖蝕位置的預(yù)測(cè)和不同重力方向?qū)_蝕的影響。因此，采用Fluent軟件進(jìn)行氣-固雙相耦合，選用E/CRC 沖蝕模型對(duì)固體顆粒的流場(chǎng)分布、沖蝕速率及沖蝕位置進(jìn)行模擬，以期為管道的沖蝕預(yù)測(cè)提供理論依據(jù)和實(shí)際參考。

1 理論模型

1.1 連續(xù)相方程

甲烷作為連續(xù)相，在Eulerian 坐標(biāo)系下求解N-S方程得到流場(chǎng)分布，其中連續(xù)性方程和動(dòng)量方程為

式中：u為流體流速，為x、y、z三個(gè)方向的分量，m/s；ρ為甲烷密度，kg/m3；t為運(yùn)行時(shí)間，s；p為壓力，Pa；τ為應(yīng)力張量；g為重力加速度，m/s2；μ為甲烷黏度，Pa·s；S為離散相相對(duì)于連續(xù)相的附加源相；I為單位張量。

1.2 離散相方程

固體顆粒作為離散相，在Lagrangian 坐標(biāo)系下通過(guò)牛頓第二定律確定顆粒的運(yùn)動(dòng)軌跡，公式為

式中：up為固體顆粒速度，m/s；ReS為相對(duì)雷諾數(shù)；ρp為顆粒密度，kg/m3；dp為顆粒直徑，m；CD為曳力系數(shù)；Fy為y軸方向其他的作用力；a1、a2、a3為離子光滑常數(shù)，可參考文獻(xiàn)[13]。

1.3 沖蝕模型

Fluent 中內(nèi)置了多種經(jīng)驗(yàn)及半經(jīng)驗(yàn)的沖蝕模型，其中E/CRC模型計(jì)算二次流彎管壁厚減薄的精度較高，且該模型充分考慮了沖擊角度、材料粗糙化對(duì)壁面的影響，因此采用該模型計(jì)算沖蝕速率，公式為

式中：ER為沖蝕速率，kg/(m2·s)；mp為顆粒質(zhì)量，kg；C(dp)為粒徑函數(shù)，取1.8×10-9；b為速度指數(shù)，取2.6；Aface為壁面面積，m2；f()θ為沖擊角函數(shù)，采用分段線性的方式進(jìn)行定義，在沖擊角0°、20°、30°、45°、90°時(shí)的函數(shù)值分別為0、0.8、1、0.5、0.4。

1.4 壁面回彈模型

固體顆粒在與壁面撞擊過(guò)程中會(huì)有能量損失和傳遞，反彈速度會(huì)小于沖擊速度，因此通過(guò)建立回彈模型預(yù)測(cè)動(dòng)量變化趨勢(shì)，采用恢復(fù)系數(shù)表征固體顆粒在碰撞前后動(dòng)量的變化。輸氣管道大多為碳鋼管線，F(xiàn)order等[14]利用鐵介質(zhì)建立的回彈模型公式為

式中：eN和eT分別為法向恢復(fù)系數(shù)和切向恢復(fù)系數(shù)。

2 數(shù)值模型

2.1 彎管模型

彎管模型直徑D=60 mm，為了避免回流對(duì)精度的影響，水平和豎直管段長(zhǎng)度均為10D=600 mm，曲率半徑為1.5D=90 mm（圖1）。

2.2 網(wǎng)格無(wú)關(guān)性分析

網(wǎng)格質(zhì)量和數(shù)量直接影響模擬計(jì)算的精度和速度，在流速6 m/s、顆粒直徑100 μm、顆粒質(zhì)量流量0.2 kg/s 的條件下，對(duì)不同網(wǎng)格數(shù)量下的沖蝕速率進(jìn)行適用性分析（圖2）。當(dāng)網(wǎng)格數(shù)量小于1 002 567 時(shí)，沖蝕速率呈無(wú)規(guī)律上下波動(dòng)趨勢(shì)，當(dāng)網(wǎng)格數(shù)量大于1 662 882 時(shí)，沖蝕速率基本保持穩(wěn)定，因此在保證計(jì)算精度的前提下，選用網(wǎng)格數(shù)量為1 662 882。對(duì)入口和出口壁面進(jìn)行膨脹層設(shè)計(jì)，劃分5 層，生長(zhǎng)因子1.2，同時(shí)對(duì)彎曲部位進(jìn)行網(wǎng)格加密（圖3）。

圖1 彎管示意圖（流動(dòng)方式1）Fig.1 Schematic diagram of elbow(flow mode 1)

圖2 網(wǎng)格無(wú)關(guān)性分析Fig.2 Grid independence analysis

2.3 邊界條件和求解方式

對(duì)固體顆粒屬性進(jìn)行簡(jiǎn)化處理，假設(shè)顆粒為密度均勻的球形，且在射入后無(wú)旋轉(zhuǎn)運(yùn)動(dòng)，同時(shí)忽略顆粒受到的虛擬質(zhì)量力、壓力梯度力和升力等。模型及邊界條件設(shè)定見表1。入口及出口的湍流強(qiáng)度和水力半徑分別為5%和60 mm，常溫常壓下甲烷密度0.668 kg/m3，固體顆粒密度2 650 kg/m3。

圖3 彎頭處的網(wǎng)格劃分Fig.3 Mesh division at the elbow

表1 模型及邊界條件設(shè)定Tab.1 Model and boundary condition setting

選用壓力-速度穩(wěn)態(tài)求解器對(duì)邊界條件進(jìn)行求解，利用SIMPLE 算法，為了得到較好的收斂效果，壓力、動(dòng)量、湍動(dòng)能采用二階迎風(fēng)算法，設(shè)置監(jiān)視器中殘差范圍為10-5。

3 結(jié)果與討論

3.1 流場(chǎng)分布

在流速4 m/s，顆粒直徑100 μm，顆粒質(zhì)量流量0.6 kg/s 的條件下，彎管的軸向截面分布云圖見圖4。水平管道的壓力和速度分布基本穩(wěn)定，受流體黏滯力的影響，速度在近壁面處的數(shù)值較小，流體經(jīng)過(guò)彎頭時(shí)固體顆粒發(fā)生直接碰撞并反彈，二次流、速度場(chǎng)和壓力場(chǎng)發(fā)生變化。其中壓力較大處出現(xiàn)在彎管外側(cè)，此時(shí)流體在彎管處產(chǎn)生離心力，對(duì)彎管內(nèi)側(cè)的牽引作用使流動(dòng)方向發(fā)生改變，內(nèi)側(cè)出現(xiàn)負(fù)壓，從內(nèi)側(cè)到外側(cè)壓力值呈遞增趨勢(shì)，內(nèi)側(cè)的壓力轉(zhuǎn)化為動(dòng)能；速度較大處出現(xiàn)在彎管內(nèi)側(cè)，因此顆粒在彎管外側(cè)和豎直管道內(nèi)發(fā)生碰撞和滑動(dòng)，故只有彎管外側(cè)和沿豎直管道方向發(fā)生沖蝕，內(nèi)側(cè)不發(fā)生沖蝕。

圖4 軸向截面分布云圖Fig.4 Distribution cloud map of axial section

3.2 流速對(duì)沖蝕的影響

在顆粒直徑100 μm，顆粒質(zhì)量流量分別為0.2、0.4、0.6、0.8、1.0 kg/s 的條件下，對(duì)不同流速下的沖蝕速率進(jìn)行模擬，結(jié)果見圖5。在顆粒質(zhì)量流量一定的條件下，沖蝕速率隨流速的增大而增大，且在小于4 m/s 時(shí)，質(zhì)量流量的變化對(duì)沖蝕速率影響不大，而大于4 m/s 時(shí)，沖蝕速率提升較快，故4 m/s 為臨界流速，在實(shí)際生產(chǎn)中尤其是顆粒質(zhì)量流量較大時(shí)應(yīng)避免超過(guò)此流速。在顆粒質(zhì)量流量較小時(shí)，沖蝕速率隨流速呈線性變化，此時(shí)流體所攜帶的顆粒數(shù)量有限，即使初始動(dòng)能較大，對(duì)管壁的損傷依然有限；在顆粒質(zhì)量流量較大時(shí)，顆粒與壁面之間碰撞的概率變大，沖蝕速率隨流速呈指數(shù)變化。

圖5 流速對(duì)沖蝕的影響Fig.5 Effect of velocity on erosion

顆粒質(zhì)量流量0.6 kg/s 的條件下，流速2、6 和10 m/s 的沖蝕云圖見圖6。在低速情況下，顆粒動(dòng)能較小，形成的沖蝕區(qū)域較小且集中在水平管段底部，隨著流速的增大，連續(xù)相對(duì)離散相的作用力增加導(dǎo)致顆粒的動(dòng)能增大，沖蝕區(qū)域變大且逐漸向豎直管線兩側(cè)擴(kuò)展。其中流速為2、4、6、8、10 m/s 時(shí)的最大沖蝕角度分別為18°、32°、43°、51°、55°，顆粒與管道壁從直接碰撞轉(zhuǎn)為滑動(dòng)碰撞，二次碰撞的區(qū)域變大，與圖6云圖中的高亮沖蝕區(qū)域相符。

圖6 不同流速下的沖蝕云圖Fig.6 Erosion cloud map at different flow rates

3.3 顆粒直徑對(duì)沖蝕的影響

在流速6 m/s，顆粒質(zhì)量流量分別為0.2、0.4、0.6、0.8、1.0 kg/s的條件下，對(duì)不同顆粒直徑的沖蝕速率進(jìn)行模擬，結(jié)果見圖7。在顆粒質(zhì)量流量一定的條件下，隨顆粒直徑增大，沖蝕速率先增大，在顆粒直徑超過(guò)100 μm 后，逐漸呈平穩(wěn)趨勢(shì)。這是由于氣-固兩相耦合作用時(shí)，小顆粒與壁面之間的碰撞更充分，而大顆粒轉(zhuǎn)移到氣體中的動(dòng)能更多，流體的攜帶能力下降，顆粒之間的有效碰撞增多，大顆粒與壁面之間的碰撞機(jī)會(huì)和強(qiáng)度均減少，同時(shí)大顆粒因重力造成的停滯作用對(duì)后續(xù)顆粒的沖蝕影響較大，因此沖蝕速率保持平穩(wěn)。

圖7 顆粒直徑對(duì)沖蝕的影響Fig.7 Effect of particle diameter on erosion

在顆粒質(zhì)量流量0.6 kg/s 的條件下，隨著顆粒直徑的增大，最大沖蝕速率和沖蝕面積均有所增加，但最大沖蝕位置基本不變。其中顆粒直徑為50、100、150、200、250 μm 時(shí)的最大沖蝕角度分別為47°、46°、46°、45°、45°，均在45°～47°之間，由于曳力作用使小顆粒的碰撞角增大，導(dǎo)致大顆粒的最大沖蝕角度略小于小顆粒，但總體上顆粒直徑對(duì)最大沖蝕位置的影響不大。圖8為顆粒直徑50、150、250 μm的沖蝕云圖。

圖8 不同顆粒直徑下的沖蝕云圖Fig.8 Erosion cloud map of different particle diameters

圖9 不同流動(dòng)方式示意圖Fig.9 Diagram of different flow modes

3.4 重力方向?qū)_蝕的影響

在實(shí)際工況中，管道的布置形式及重力場(chǎng)的方向直接影響沖蝕區(qū)域。分別對(duì)4種重力場(chǎng)的的彎管沖蝕進(jìn)行模擬，其中一種攜砂氣體流動(dòng)方式如圖1所示，定義為流動(dòng)方式1，其余3 種攜砂氣體流動(dòng)方式見圖9。

圖10 不同位置的沖蝕云圖Fig.10 Erosion cloud maps at different locations

在流動(dòng)方式1 和流動(dòng)方式4 中，重力對(duì)顆粒運(yùn)動(dòng)起減速作用；流動(dòng)方式2 和流動(dòng)方式3 中，重力對(duì)顆粒運(yùn)動(dòng)起加速作用，沖蝕云圖見圖10。流動(dòng)方式2 和流動(dòng)方式3 對(duì)彎管的沖蝕程度遠(yuǎn)大于流動(dòng)方式1 和4，其中流動(dòng)方式4 受到的沖蝕磨損最輕微。位置3流體攜帶顆粒撞擊彎管后，由于重力作用，勢(shì)能增大，二次流能量較大；流動(dòng)方式2流體先沿重力方向豎直向下運(yùn)動(dòng)，勢(shì)能較大，撞擊彎管后將損失一部分能量，故二次流能量減小；流動(dòng)方式1流體在撞擊彎管后為克服重力作用，二次流能量更??；流動(dòng)方式4流體需先克服重力作用，然后再撞擊彎管，二次流能量最小。因此，當(dāng)入口或出口方向與重力方向起正向疊加作用時(shí)，沖蝕較嚴(yán)重，沖蝕嚴(yán)重程度由大到小依次為流動(dòng)方式3、流動(dòng)方式2、流動(dòng)方式1、流動(dòng)方式4。

重力場(chǎng)的方向不僅會(huì)影響最大沖蝕速率，也會(huì)影響沖蝕區(qū)域。對(duì)比流動(dòng)方式1 和流動(dòng)方式3，在入口水平流動(dòng)的過(guò)程中，由于重力作用，在靠近彎頭入口區(qū)域的管道下方存在一定的輕微磨損，但迎角較小，主要以切削作用為主；流動(dòng)方式1在進(jìn)入彎頭區(qū)域后立刻造成了較大的沖蝕，而流動(dòng)方式3在進(jìn)入彎頭區(qū)域的后方才造成較大沖蝕，這與重力方向帶來(lái)的顆粒偏流有關(guān)，且顆粒的動(dòng)能越大，差異越明顯。對(duì)比流動(dòng)方式2 和流動(dòng)方式4，在入口豎直流動(dòng)的過(guò)程中，無(wú)論重力與彎管布置方向的關(guān)系以及攜砂天然氣的流動(dòng)方向如何，因無(wú)偏流，入口直管段幾乎沒有沖蝕損傷，但流動(dòng)方式2與重力方向相同，因此在進(jìn)入彎頭區(qū)域后立刻造成了較大的沖蝕，沖蝕區(qū)域集中在彎頭及出口管道上。

4 結(jié)論

通過(guò)對(duì)彎管內(nèi)攜砂氣體流速、顆粒質(zhì)量流量、顆粒直徑等進(jìn)行模擬分析，發(fā)現(xiàn)顆粒質(zhì)量流量較小時(shí)，沖蝕速率隨流速呈線性變化且變化幅度不大，顆粒質(zhì)量流量較大時(shí)，沖蝕速率隨流速呈指數(shù)變化；對(duì)比了不同流速、不同顆粒直徑的最大沖蝕位置，其中最大沖蝕位置隨流速的增加而增大，而顆粒直徑對(duì)最大沖蝕位置的影響不大；分析了4種不同重力場(chǎng)條件下的沖蝕情況，當(dāng)入口或出口方向與重力方向起正向疊加作用時(shí)，沖蝕較嚴(yán)重。