基于CFD軟件的管道彎頭沖蝕數(shù)值模擬及彎頭測厚布點(diǎn)優(yōu)化

廖飛龍,劉　炯,周詠琳,趙　瀟,萬　夫,周兆明

(1.中國石油集團(tuán)川慶鉆探工程有限公司 安全環(huán)保質(zhì)量監(jiān)督檢測研究院,廣漢 618300；2.中石油西南管道公司 蘭州輸油氣分公司,蘭州 730000)

?

基于CFD軟件的管道彎頭沖蝕數(shù)值模擬及彎頭測厚布點(diǎn)優(yōu)化

廖飛龍1,劉炯1,周詠琳1,趙瀟2,萬夫1,周兆明1

(1.中國石油集團(tuán)川慶鉆探工程有限公司 安全環(huán)保質(zhì)量監(jiān)督檢測研究院,廣漢 618300；2.中石油西南管道公司 蘭州輸油氣分公司,蘭州 730000)

彎頭作為改變介質(zhì)流向的管件，沖蝕是其失效的主要原因。應(yīng)用計(jì)算流體動力學(xué)CFD軟件對彎頭的沖蝕過程進(jìn)行了模擬仿真，通過分析仿真結(jié)果提出了彎頭測厚布點(diǎn)的優(yōu)化方案?，F(xiàn)場的檢測結(jié)果表明：優(yōu)化后的布點(diǎn)方案能夠全面掌握彎頭外弧面的沖蝕情況，避免了漏檢的發(fā)生。

CFD;彎頭;沖蝕;檢測;優(yōu)化

壓力管道是具有爆炸危險的特種承壓設(shè)備，其內(nèi)部承載著高壓、高溫、低溫、易燃易爆、有毒或者腐蝕性的介質(zhì)。隨著服役時間的增加,因沖蝕磨損減薄導(dǎo)致管道失效而引發(fā)的事故越來越多[1-3]。因此如何發(fā)現(xiàn)管道存在的缺陷，確保管道的正常運(yùn)行成為近幾年來亟待解決的問題。

在管路系統(tǒng)中，彎頭是改變介質(zhì)方向的部件。當(dāng)介質(zhì)經(jīng)過彎頭時，其速度、壓力會發(fā)生較大的變化，從而導(dǎo)致流體的流動不穩(wěn)定[4]。目前針對彎頭檢測的主要手段是定點(diǎn)測厚，對于彎頭現(xiàn)場檢測的布點(diǎn)主要是根據(jù)標(biāo)準(zhǔn)或者經(jīng)驗(yàn)來確定的，缺乏一定的針對性。因此，筆者以新疆某含砂氣井井口的彎頭為例，基于ANSYS-CFD軟件對含砂氣沖蝕彎頭的過程進(jìn)行數(shù)值模擬分析，根據(jù)彎頭的沖蝕情況提出彎頭測厚布點(diǎn)的優(yōu)化方案。

1　流體理論與數(shù)學(xué)模型

1.1能量守恒方程

能量守恒定律是包含有熱交換的流動系統(tǒng)必須滿足的基本定律。該定律可表述為：微元體中能量的增加率等于進(jìn)入微元體的凈熱流量加上體積力與表面力對微元體所做的功。該定律實(shí)際上是熱力學(xué)第一定律,即：流體的能量通常是內(nèi)能、動能和勢能三項(xiàng)之和，內(nèi)能和溫度之間存在一定關(guān)系，因此可得到以溫度為變量的能量守恒方程：

(1)

式中：cp為比熱容；ρ為液體的密度;u為液體的速度;T為溫度；k為流體的傳熱系數(shù)；ST為流體的內(nèi)熱源及由于粘性作用流體機(jī)械能轉(zhuǎn)換為熱能的部分,有時稱ST為粘性耗散能。

1.2多相流分析

多相流分析中，各相之間在宏觀尺度上混合，該混合尺度遠(yuǎn)小于網(wǎng)格尺度，但是遠(yuǎn)大于分子尺度[5]。所有的相占有同一空間體積，在控制體內(nèi)假設(shè)每一相占有的體積大小用變量體積分?jǐn)?shù)來表示；每一相有自己的流場參數(shù)，各相通過相間的能量傳輸、動量傳輸、質(zhì)量傳輸模型耦合。每一相具有各自的質(zhì)量、動量、能量傳輸方程。

(1) 質(zhì)量連續(xù)方程為：

(2)

(2) 動量傳輸方程為：

(3)

式中：:Uβ為β相的速度;μα為α相的剪切粘性系數(shù)。

連續(xù)相和離散相之間的相互作用力是通過相間速度差、連續(xù)相的屬性及相間界面積進(jìn)行計(jì)算的，作用在顆粒/分散相上的力包括曳力和非曳力。其中,曳力通過無量綱曳力系數(shù)與顆粒雷諾數(shù)相關(guān)。

(4)

(5)

式中：:CD為曳力系數(shù);Re為雷諾數(shù);dp為顆粒直徑。

在CFX(CFD中的一種)軟件中：

(6)

1.3Finnie磨損腐蝕模型

壁面的磨損腐蝕主要是由于流體介質(zhì)中的顆粒碰撞侵蝕作用引起的。關(guān)于顆粒、顆粒的碰撞和壁面的屬性存在著比較復(fù)雜的函數(shù)關(guān)系[6]。對于幾乎所有金屬來說，磨損腐蝕程度隨著顆粒的入侵角和速度變化而變化，其函數(shù)關(guān)系是：

(7)

入侵角度是入侵的顆粒軌跡與壁面的夾角。指數(shù)n通常取2.3～2.5間的值。

Finnie磨損腐蝕模型[6]將磨損速率和入侵顆粒的動能聯(lián)系起來，取n的值為2，即：

(8)

(9)

(10)

2　含砂氣沖蝕彎頭的數(shù)值模擬分析及結(jié)果

2.1建立CFD模型

選用最常用的90°彎頭作為研究對象。以新疆某含砂氣井井口的彎頭為例，彎頭公稱直徑為114 mm，彎頭壁厚為8 mm。設(shè)置速度為入口邊界條件，壓力為出口邊界條件。由于湍流運(yùn)動中流體質(zhì)點(diǎn)的運(yùn)動情況十分復(fù)雜，因此借助試驗(yàn)數(shù)據(jù)用經(jīng)驗(yàn)公式進(jìn)行近似表達(dá):使用One Seventh Power Law(1/7次冪方律),其表達(dá)式為U=Wmax(1-r/Rmax)1/7(U為入口速度;Wmax為最大速度;Rmax為最大管徑;r為距離管中心的距離)，通過插入Expression函數(shù)的方法來設(shè)置。圖1為彎頭的結(jié)構(gòu)模型和流體模型。

圖1　彎頭的結(jié)構(gòu)模型和流體模型

2.2仿真結(jié)果與分析

(1) 壓力場

圖2為彎頭壓力的分布云圖。由圖可見:壓力在彎頭的外弧面處達(dá)到最大值，沿著彎頭徑向逐漸降低，至彎頭內(nèi)弧面時達(dá)到最小值。從數(shù)值上看,壓力在彎頭內(nèi)弧面處出現(xiàn)了負(fù)壓。當(dāng)管道內(nèi)介質(zhì)通過彎頭到達(dá)外弧面的位置時，介質(zhì)的流向發(fā)生改變，因此在外弧面位置會產(chǎn)生最大的壓力值。

圖2　彎頭壓力的分布云圖

(2) 砂粒的運(yùn)動軌跡

圖3為砂粒相在管道中的運(yùn)動軌跡。顆粒進(jìn)入彎頭之前做直線運(yùn)動，進(jìn)入彎頭后，砂粒以一定的動量沖擊彎頭的內(nèi)壁面，由圖可看到砂粒主要與彎頭的外弧面發(fā)生碰撞。相對于外弧面,內(nèi)弧面區(qū)域基本上沒有砂粒。在這一過程中，砂粒以不同的角度與外弧內(nèi)壁面發(fā)生碰撞，并且在該區(qū)域聚集了大量的砂粒，因此彎頭的外弧面區(qū)域是極易發(fā)生沖蝕的地方。

圖3　砂粒相的粒子運(yùn)動軌跡

(3) 沖蝕特性

通過在CFX-Post中創(chuàng)建Iso-Clip面(截取斷面)來提取彎頭上受到?jīng)_蝕的區(qū)域,如圖4所示。通過彎頭的沖蝕區(qū)域云圖來看，彎頭的沖蝕區(qū)域是以外弧面為中心向彎頭的兩頰面擴(kuò)散開的。

圖4　彎頭的沖蝕區(qū)域云圖

3　彎頭測厚布點(diǎn)方案優(yōu)化及現(xiàn)場應(yīng)用

3.1彎頭測厚布點(diǎn)方案優(yōu)化

目前，檢測彎頭的主要手段是使用超聲波測厚儀對其進(jìn)行測厚。如圖5所示，在現(xiàn)場彎頭檢測過程中，一般在每一測厚截面上,選擇上下左右4個測厚方位，即順介質(zhì)流向在同一截面上看，測厚點(diǎn)順時針方向分布在0點(diǎn)、3點(diǎn)、6點(diǎn)、9點(diǎn)四個方位，即圖中A,B,C,D 4個測厚點(diǎn)。它的確定原則遵循右手法則:右手大拇指指向介質(zhì)流向,4指彎曲的方向即是測厚點(diǎn)A-B-C-D的順序方向。而對于彎頭，彎頭的外弧面為A面。在同一個時鐘點(diǎn)(A、B、C、D方向)上，測厚點(diǎn)順介質(zhì)流向以秩序連續(xù)編號來布置。

圖5　現(xiàn)用的測厚布點(diǎn)方案

通過分析仿真結(jié)果可知,彎頭的沖蝕區(qū)域主要在外弧面。因此，現(xiàn)有的測厚布點(diǎn)方案不能完全涵蓋彎頭的沖蝕區(qū)域,從而易導(dǎo)致漏檢，存在一定的隱患和風(fēng)險。因此,筆者基于CFX的仿真結(jié)果對彎頭測厚布點(diǎn)的方案進(jìn)行優(yōu)化。優(yōu)化方案如圖6所示，即在彎頭外弧面的中心線(圖中的A方向)和彎頭頰面(圖中B方向)之間增加一個布點(diǎn)方向，定義為AB向。同理，在A向和D向之間增加AD向，從而可針對彎頭外弧面的沖蝕區(qū)域進(jìn)行較為全面的檢測。

圖6　測厚布點(diǎn)優(yōu)化方案

3.2現(xiàn)場檢測結(jié)果

采用測厚布點(diǎn)優(yōu)化方案,對新疆某油田部分單井的二級節(jié)流后生產(chǎn)管線進(jìn)行超聲波測厚檢測。該生產(chǎn)管線的材料是316L鋼，規(guī)格是φ114 mm×11 mm，其中彎頭的初始壁厚為14 mm。在現(xiàn)場檢測工作中發(fā)現(xiàn)有6個彎頭的最小壁厚出現(xiàn)在外弧面的中心線和彎頭頰面之間的區(qū)域，如圖7所示，圖中圓點(diǎn)標(biāo)示處為該彎頭最小壁厚的位置。其中有一個彎頭的厚值減薄為8.55 mm，不能通過強(qiáng)度校核而需立即更換。因此現(xiàn)場的檢測結(jié)果表明，使用該方案對彎頭進(jìn)行布點(diǎn)測厚能較為全面地掌握彎頭的沖蝕情況，有效地避免彎頭的漏檢情況。

圖7　彎頭最小壁厚標(biāo)示圖

4　結(jié)論

(1) 彎頭的外弧面處存在著最大壓力值，含砂氣攜帶著砂粒在該區(qū)域大量聚集，并且在彎頭外弧面區(qū)域沖蝕磨損最嚴(yán)重。

(2) 通過在彎頭外弧面的中心線和彎頭頰面之間增加布點(diǎn)方向的方法,對測厚布點(diǎn)方案進(jìn)行優(yōu)化?，F(xiàn)場的檢測結(jié)果表明,使用優(yōu)化后的方案能夠全面掌握彎頭外弧面的沖蝕情況，從而有效避免彎頭的漏檢。

[1]魏秀芝. 氣力輸送顆粒在輸送管彎頭中的運(yùn)動及磨損[J].黑龍江石油化工,1996(2):36-39.

[2]高福慶. 管道檢測的必要性[J].管道技術(shù)與設(shè)備，1998(1):40-42.

[3]郎需慶，趙志勇，宮宏，等.油氣管道事故統(tǒng)計(jì)分析與安全運(yùn)行對策[J].安全、健康和環(huán)境，2006,6(10)：15-17.

[4]曾涌捷.天然氣管道彎頭沖蝕失效機(jī)理研究[J].石油與化工設(shè)備,2011,14(2):44-46.

[5]李明高，李明.ANSYS 13.0流場分析技術(shù)及應(yīng)用實(shí)例[M]. 北京: 機(jī)械工業(yè)出版社，2012:349.

[6]謝龍漢，趙新宇，張炯明.ANSYS CFX流體分析及仿真[M].北京：電子工業(yè)出版社，2012:152-157.

Numerical Simulation of Erosion in Elbow Based on CFD and Optimization of Thickness Measurement

LIAO Fei-long1,LIU Jiong1,ZHOU Yong-lin1,ZHAO Xiao2,WAN Fu1,ZHOU Zhao-ming1

(1.Safety,Environment ,Quality Supervision and Testing Research Institute, CNPC Chuanqing Drilling Engineering Co.,Ltd.,Guanghan 618300,China;2.Lanzhou Oil and Gas Transmission Sub-company,PetroChina Southwest Pipeline Company,Lanzhou 730000,China)

The role of elbow is to change the flow of medium. The failure of elbow is mainly due to erosion. The numerical simulation of the erosion process about four direction of medium in elbow and horizontal elbow is established by using CFX. The optimization of thickness measurement is based on the simulation analysis outcomes. The on-site testing results show that the optimized method can get full grasp of the erosion status of outer surface of elbow. The undetected errors are avoided effectively.

CFD;Elbow;Erosion;Detection;Optimization

2015-12-17

廖飛龍(1988-)，男，助理工程師，主要從事特種設(shè)備及鉆機(jī)的檢驗(yàn)工作。

廖飛龍,E-mail：liaofeilong311@163.com。

10.11973/wsjc201608008

TG115.28

A

1000-6656(2016)08-0034-04