謝明1,2,3  吳貴陽(yáng)1,2,3  陳偉4

1. 中國石油西南油氣田分公司天然氣研究院； 2. 國家能源高含硫氣藏開(kāi)采研發(fā)中心；

3. 中國石油天然氣集團公司高含硫氣藏開(kāi)采先導試驗基地； 4. 中國石化西南油氣分公司川東北采氣廠(chǎng)

摘  要：為了定量描述液體在攪拌釜內產(chǎn)生的剪切力與現場(chǎng)管線(xiàn)受到的液體剪切力的對應關(guān)系，通過(guò)fluent軟件模擬計算了液體在直管段、彎管及攪拌釜內的剪切力分布情況。結果表明，液體在直管段內的剪切力呈階梯狀分布，而彎管內的液體剪切力在彎管內側和流體流向12點(diǎn)鐘方向最大，攪拌釜內的最大液體剪切力則位于槳葉處。攪拌釜內最大液體剪切力與彎管內最大剪切力相當，而釜壁所受剪切力情況與流速較低時(shí)的直管段內所受剪切力情況相當。據此，可開(kāi)展不同條件下的沖刷腐蝕模擬評價(jià)試驗。

關(guān)鍵詞：fluent； 數值模擬； 剪切力； 流速

頁(yè)巖氣開(kāi)發(fā)過(guò)程中，大量石英砂被注入地下裂縫作為支撐劑使用。隨著(zhù)頁(yè)巖氣井投產(chǎn)，井底部分石英砂隨氣液介質(zhì)從井底進(jìn)入地面集輸系統。生產(chǎn)過(guò)程中，管線(xiàn)受氣、液、固三相流體沖擊與介質(zhì)腐蝕的協(xié)同作用，造成腐蝕加速，即沖刷腐蝕。對于靜態(tài)流體，金屬在腐蝕介質(zhì)作用下產(chǎn)生腐蝕產(chǎn)物膜覆蓋在金屬表面，具有一定的緩蝕效果，而流體的沖擊力造成的腐蝕產(chǎn)物膜脫落，則進(jìn)一步加劇了腐蝕。部分學(xué)者基于數值模擬計算方法，研究了一定工況條件下，導致腐蝕概率變大的原因[1,2]。目前，研究沖刷腐蝕的設備主要有動(dòng)態(tài)高溫高壓釜試驗裝置、旋轉電極試驗裝置、射流沖刷腐蝕試驗裝置、管流式?jīng)_刷腐蝕試驗裝置、科里奧利沖蝕試驗裝置[3]，最常用的為動(dòng)態(tài)高溫高壓試驗裝置，但利用該裝置模擬生產(chǎn)管線(xiàn)內介質(zhì)流型、流態(tài)的適應性并沒(méi)有深入研究案例。本文通過(guò)模擬計算液體在攪拌釜內產(chǎn)生的剪切力與現場(chǎng)管線(xiàn)受到的液體剪切力的對應關(guān)系，有利于針對現場(chǎng)生產(chǎn)情況開(kāi)展不同條件下的沖刷腐蝕模擬評價(jià)試驗。

1 數學(xué)描述

流體控制采用標準k-ε雙方程模型，標準k-ε模型自從被Launder和Spalding提出后，就變成工程流場(chǎng)計算中的主要工具，它是個(gè)半經(jīng)驗公式。標準k-ε雙方程模型假定湍動(dòng)粘度μi是各向同性的，其基本輸運方程為：

2 數值模擬

2.1 物理模型

直管段管內徑為200 mm，粗糙度0.05 mm，長(cháng)度2 m，水為流動(dòng)介質(zhì)，網(wǎng)格采用三維六面體網(wǎng)格劃分方式，整個(gè)模型從入口到出口采用掃略方法劃分，網(wǎng)格尺寸為2.0 mm，網(wǎng)格節點(diǎn)數為199 980，網(wǎng)格單元為192 357，如圖1所示。

圖1  直管段模型網(wǎng)格劃分

彎管內徑為200 mm，粗糙度0.05mm，中心距305 mm，水為流動(dòng)介質(zhì)，網(wǎng)格采用三維六面體網(wǎng)格劃分方式，整個(gè)模型從入口到出口采用掃略方法劃分，網(wǎng)格尺寸為2.0mm，網(wǎng)格節點(diǎn)數為435 305，網(wǎng)格單元為421 360，如圖2所示。

圖2  彎管模型網(wǎng)格劃分

攪拌器槳葉的直徑為60 mm，槳葉厚度10 mm，槳葉寬度3 mm，將攪拌槽整體分為轉子區和定子區兩部分，如圖3所示。采用非結構網(wǎng)格四面體網(wǎng)格對轉子區和定子區分別進(jìn)行劃分，網(wǎng)格尺寸為0.5 mm，轉子區網(wǎng)格節點(diǎn)數為15 074，網(wǎng)格單元為77 158。定子區網(wǎng)格節點(diǎn)數為125 515，網(wǎng)格單元為684 133。

圖3  轉子區與定子區

2.2 直管段模擬計算

利用fluent軟件模擬液體流速為1 m/s、3 m/s、5 m/s、8 m/s、10 m/s條件下直管段管壁剪切力受力情況，

結果可知，直管段管壁剪切力沿流動(dòng)方向逐漸減小，隨流速增加管壁剪切力逐漸增大。當液體流速達到10 m/s時(shí)，管壁最大剪切力可達到236 Pa，如圖4所示。

圖4  直管段管壁最大剪切力隨液體流速變化趨勢

2.3 彎管模擬計算

以DN200彎管為研究對象，模擬液體流速為1 m/s、3 m/s、5 m/s、8 m/s、10 m/s條件下管壁剪切力受力情況（液體流向從左至右，不考慮重力影響），

結果可知彎管管壁剪切力受力分布不均，其受力最大點(diǎn)為內側管壁和流出段外側管壁，流速增加管壁剪切力逐漸增大。當液體流速達到10 m/s時(shí)，管壁最大剪切力可達到500 Pa，如圖5所示。

圖5  彎管管壁最大剪切力隨液體流速變化趨勢

2.4 攪拌釜模擬計算

實(shí)驗室攪拌釜內徑100 mm，攪拌器直徑60 mm，寬度3 mm，在100 r/min、500 r/min、1000 r/min、1500 r/min、2000 r/min、2500 r/min轉速條件下，計算分析釜壁剪切力與中心截面速度變化情況，

通過(guò)攪拌釜模擬計算可知，隨著(zhù)攪拌槳轉速提高，釜壁最大剪切力在轉速為2500 r/min時(shí)，可達140 Pa，而攪拌槳外沿的剪切力最大可達464 Pa。當釜內液體流速在4 ~5m/s，攪拌釜中心流速可達7 ~8 m/s，但流體分布不均勻，如圖6、7所示。

圖6  不同轉速下釜壁最大剪切力

圖7  不同轉速下攪拌槳外沿剪切力變化

3  結論

基于數值模擬計算結果，結合頁(yè)巖氣生產(chǎn)工況，可利用動(dòng)態(tài)高溫高壓試驗裝置開(kāi)展不同條件下的沖刷腐蝕模擬評價(jià)試驗。

（1）攪拌槳周?chē)�液體流速可達4 ~5 m/s，中心流速可達7 ~8 m/s，可用以模擬該流速條件下生產(chǎn)管線(xiàn)內的腐蝕狀況。

（2）攪拌釜釜壁受到的剪切力與直管段在6 ~7m/s和彎管段在4 ~5m/s液體流速條件下受到的剪切力相當，而攪拌槳中心受到的剪切力完全可以模擬液體流速在10 m/s條件下，直管段和彎管所受到的剪切力。

基于頁(yè)巖氣生產(chǎn)井底出砂的情況，結合上述數值模擬計算方法，可對頁(yè)巖氣集輸系統腐蝕情況開(kāi)展不同溫度、壓力、流速、介質(zhì)、砂含量等工況條件下的腐蝕模擬試驗。

參考文獻：

[1] 劉勇峰，吳明，趙玲等.凝析氣田集輸管道彎管沖刷腐蝕數值計算[J].腐蝕與防護，2012，33（2）：132-135.

[2] 茅俊杰.氣液兩相流管道沖刷腐蝕的研究[D].山東大學(xué)，2012.

[3] 朱娟，張喬斌，陳宇等.沖刷腐蝕的研究現狀[J].中國腐蝕與防護學(xué)報，2014，34（3）：199-207.

作者：謝明，男， 1989年生，四川南充人，中國石油西南油氣田分公司天然氣研究院腐蝕與防護研究所助理工程師，主要從事頁(yè)巖氣集輸系統腐蝕狀況研究工作。

基金項目：中國石油西南油氣田分公司2017年科研項目“長(cháng)寧頁(yè)巖氣集輸系統含砂條件下腐蝕分析與控制措施研究”（20170304-07）。

《管道保護》2018年第1期（總第38期）