賴曉明

摘 要：針對國內(nèi)外高壓管匯件產(chǎn)品很難滿足我國頁巖氣商業(yè)化開采，面臨超高壓壓裂作業(yè)難題，在實際工程作業(yè)中陸續(xù)出現(xiàn)的一些質(zhì)量問題，本課題將采用CFD軟件仿真分析固-液兩相流動，運用數(shù)值模擬方法能獲得在各種工況下的完整數(shù)據(jù)，考慮固相顆粒之間的作用和固-液相之間的耦合作用力，分別對顆粒物含量較低時和固相顆粒含量較高的壓裂作業(yè)流體工況進行分析，優(yōu)化現(xiàn)有產(chǎn)品結(jié)構(gòu)，解決一些重要的國產(chǎn)管匯件使用壽命不穩(wěn)定的難題。

關(guān)鍵詞：頁巖氣；壓裂；CFD分析；數(shù)學(xué)模型

中圖分類號：TE345 文獻標(biāo)識碼：A 文章編號：1671-2064（2019）11-0097-02

0 引言

我國在油氣勘探開發(fā)的廣度和深度方面呈現(xiàn)出大幅的上漲態(tài)勢，特別是頁巖氣油氣資源的勘探開發(fā)，對壓裂設(shè)備提出了前所未有的挑戰(zhàn)，2015年左右，由于原有高壓管匯性能滿足不了作業(yè)要求，分別在新疆油田和川慶鉆探等作業(yè)區(qū)塊，連續(xù)發(fā)生了多次高壓管匯失效事故。另外，近20年來，計算機運用技術(shù)取得的飛速發(fā)展，其中流體力學(xué)計算（Computational Fluid Dynamics，CFD）方面的研究也取得了前所未有的進步，研究逐步成為了一種與實驗研究并存的有效手段。

1 研究目標(biāo)

以為內(nèi)徑為3寸，工作壓力約140Mpa的高壓彎頭為對象。這里為該樣式的高壓彎頭構(gòu)建了四種空間角度位置不同的兩彎頭相聯(lián)接方式，分別為兩彎頭所在平面夾角為0°，45°，90°和135°。針對這些不同的聯(lián)接方式，首先使用FLUENT所帶前置網(wǎng)格軟件GAMBIT建立相應(yīng)的計算區(qū)域，各聯(lián)接方式彎頭的計算區(qū)域均由4部分組成：進口直管段、彎管1、彎管2及出口直管段。

2 網(wǎng)格劃分

運用計算機對樣品進行三維建模，用GAMBIT軟件對模型進行本體材料的網(wǎng)格劃分，完成計算空間的離散創(chuàng)建。采用FLUENT軟件對流場進行有限體積法計算，實現(xiàn)計算的基礎(chǔ)和前提是網(wǎng)格劃分，其結(jié)果將直接影響到最終數(shù)值計算的穩(wěn)定性、高效性和結(jié)果的準確性。在GAMBIT軟件時，為使計算更加準確、快捷，可以采用六面體網(wǎng)格劃分，并在管壁處添加適當(dāng)?shù)倪吔鐚泳W(wǎng)格以便能更好的分析近管壁處的流動情況。所研究的四種不同聯(lián)接方式模型圖如圖1所示。

3 邊界條件及多相物性

把速度進口（velocity inlet）設(shè)定為管線進口端邊界條件，根據(jù)實測的流量值計算得到流場入口速度；管段出口邊界條件設(shè)置為壓力出口（pressure outlet）；管壁設(shè)置為無滑移壁面（no-slip wall）。計算中所用材料相關(guān)物性為：支撐劑為30/50目（視密度：3.22g/cm3，體積密度：1.8g/cm3）；100目（視密度：3.14g/cm3，體積密度：1.67g/cm3），壓裂液密度：1.0g/cm3，粘度：0.08pa·s。

4 多相流模型設(shè)置

當(dāng)使用離散相模型時，采用雙向耦合方式計算，液相湍流模型中選用標(biāo)準雙方程模型，顆粒運動阻力計算時選用球形阻力公式，顆粒粒徑采用均勻分布，管壁邊界條件設(shè)為反射邊界，選用隨機追蹤模型預(yù)測由液相湍流引起的顆粒彌散。

壓力速度耦合求解方法選用SIMPLEC方法，動量方程，湍動能及湍動能耗散率方程選用二階迎風(fēng)格式進行離散。各松弛因子選用默認設(shè)置。

當(dāng)使用雙流體模型時，選擇Eulerian雙流體模型，湍流選擇標(biāo)準多相混合模型（Mixture Multiphase Model）的雙方程模型。固相顆粒粘度選用syamlal-obrien模型，顆粒體積粘度選用lun-et-al模型，忽略顆粒摩擦粘度，固體壓力選用lun-et-al模型，徑向分布函數(shù)選用lun-et-al模型，相間曳力選用gidaspow模型。

壓力速度耦合求解方法選用PC-SIMPLE（Phase Coupled SIMPLE）方法，體積分數(shù)方程選用一階迎風(fēng)格式；動量方程，湍動能及湍動能耗散率方程選用二階迎風(fēng)格式進行離散。各松弛因子選用默認設(shè)置。

5 模擬結(jié)果及分析

對于目前確定的四種不同類型的彎管結(jié)構(gòu)，分別模擬了30、50、100目等三種不同粒度顆粒多相流的情況。計算所設(shè)液流速度為15m/s，顆粒在計算區(qū)域入口處的流速設(shè)為與液流相同，顆粒流量設(shè)為3.6Kg/s。

分析顆粒粒度取為30目，彎管夾角為0°時的流場分布。彎管部分外側(cè)的壓應(yīng)力高于內(nèi)側(cè)，在內(nèi)側(cè)壓應(yīng)力沿流動方向先降后升，在外側(cè)應(yīng)壓力沿流動方向先升后降。因此，液流速度會在靠近內(nèi)側(cè)處先增加后減少，在靠近外側(cè)處先減少后增加，液流相速度的急劇變化會加強顆粒相湍流程度，使得顆粒更易在管壁處發(fā)生沖蝕磨損。

由于固-液兩相存在密度差，在沿流動方向上顆粒濃度逐漸增大，在彎管處兩相出現(xiàn)了不同的運動形式，固-液兩相在離心力作用下分離使固相顆粒集中在偏向彎管外側(cè)，造成外側(cè)壁面的顆粒濃度大幅增加。另外，在該處彎管兩側(cè)各有一個近似對稱的顆粒濃度較高區(qū)域，這是彎管內(nèi)的二次流造成的。

圖2分別給出了當(dāng)固相顆粒粒度取30目的情況下，三種彎管夾角分別為45°、90°和135°時彎管內(nèi)的顆粒磨損率分布圖。從圖中可以看出磨損區(qū)域主要分布在彎管外側(cè)管壁附近，在內(nèi)側(cè)管壁區(qū)域幾乎沒有磨損發(fā)生，且隨著角度的增大，磨損嚴重區(qū)向彎管出口方向移動。

表1給出了三種不同粒度顆粒分別在四種不同彎管結(jié)構(gòu)中模擬計算的壓降及磨損率。從表中可以看出，顆粒粒度對壓降的影響不大，彎管4壓降彎管3彎管2彎管1；對于同種顆粒粒度來說，彎管2結(jié)構(gòu)的平均磨損率最大。

6 討論及小結(jié)

彎管的結(jié)構(gòu)形式，即在生產(chǎn)中表現(xiàn)為管匯件的連接方式，對管匯件的耐沖蝕性能有著顯著影響。不同連接方式，可能導(dǎo)致產(chǎn)品使用壽命出現(xiàn)倍數(shù)差距。

參考文獻

[1] 張繼信，樊建春，汪彤，姜心.壓裂液對高壓彎管沖蝕作用的數(shù)值分析[J].潤滑與密封，2013（4）：27-31.

[2] 陳春，聶松林，吳正江，李壯云.高壓水射流CFD仿真分析液壓[J].機床與液壓，2006（2）：103-105.

[3] 鄧世彪，邱杰，張宏，等.大型壓裂施工裝備優(yōu)化配套[J].石油機械，2011，39（9）：89-91.