廖飛龍,劉　炯,周詠琳,趙　瀟,萬(wàn)　夫,周兆明

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

?

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

廖飛龍1,劉炯1,周詠琳1,趙瀟2,萬(wàn)夫1,周兆明1

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

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

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

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

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

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

1.1能量守恒方程

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

(1)

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

1.2多相流分析

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

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

(2)

(2) 動(dòng)量傳輸方程為：

(3)

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

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

(4)

(5)

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

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

(6)

1.3Finnie磨損腐蝕模型

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

(7)

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

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

(8)

(9)

(10)

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

2.1建立CFD模型

選用最常用的90°彎頭作為研究對(duì)象。以新疆某含砂氣井井口的彎頭為例，彎頭公稱直徑為114 mm，彎頭壁厚為8 mm。設(shè)置速度為入口邊界條件，壓力為出口邊界條件。由于湍流運(yùn)動(dòng)中流體質(zhì)點(diǎn)的運(yùn)動(dòng)情況十分復(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為距離管中心的距離)，通過(guò)插入Expression函數(shù)的方法來(lái)設(shè)置。圖1為彎頭的結(jié)構(gòu)模型和流體模型。

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

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

(1) 壓力場(chǎng)

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

圖2　彎頭壓力的分布云圖

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

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

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

(3) 沖蝕特性

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

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

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

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

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

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

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

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

3.2現(xiàn)場(chǎng)檢測(cè)結(jié)果

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

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

4　結(jié)論

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

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

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

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

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

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

[5]李明高，李明.ANSYS 13.0流場(chǎng)分析技術(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