幺　成 ， 李介普 ， 李　翔

(1.中國(guó)石油大學(xué)(北京) ， 北京　102249 ； 2.中國(guó)特種設(shè)備檢測(cè)研究院 ， 北京　100029)

?

石化管道沖蝕磨損的數(shù)值模擬分析

幺成1,2， 李介普1,2， 李翔2,1

(1.中國(guó)石油大學(xué)(北京) ， 北京102249 ； 2.中國(guó)特種設(shè)備檢測(cè)研究院 ， 北京100029)

以某煉廠柴油管線為研究對(duì)象，對(duì)柴油流經(jīng)的管線進(jìn)行流場(chǎng)分析。利用Fluent對(duì)該柴油管線進(jìn)行數(shù)值模擬計(jì)算。研究結(jié)果表明：該管線沖蝕最嚴(yán)重的部位在彎管與盲三通之間靠近彎管內(nèi)拱壁面中心位置，此處受到的剪切應(yīng)力最大。該研究可指導(dǎo)石化管道的檢測(cè)，可以有效提高輸氣管道系統(tǒng)運(yùn)行的安全性、可靠性，并能顯著節(jié)約管道系統(tǒng)的檢測(cè)成本。

管道 ； 沖蝕 ； 磨損 ； 數(shù)值模擬

管道有多種失效形式，其中由沖蝕磨損造成破壞引起的管道失效是最常見(jiàn)、最廣泛的破壞形式。管道的沖蝕磨損不僅嚴(yán)重影響了油氣集輸，同時(shí)也造成了巨大的經(jīng)濟(jì)損失。石化集輸管道系統(tǒng)中幾何結(jié)構(gòu)的復(fù)雜性導(dǎo)致其受力情況比較復(fù)雜。管道內(nèi)沖蝕磨損機(jī)理較為復(fù)雜，目前已有的沖蝕計(jì)算方法也是在實(shí)驗(yàn)的基礎(chǔ)上建立的，具有一定的局限性，根據(jù)經(jīng)驗(yàn)公式所得到的結(jié)果往往也存在爭(zhēng)議。因此對(duì)某石化管道沖蝕磨損進(jìn)行數(shù)值模擬分析具有重要的現(xiàn)實(shí)意義。

Blatt W等[1]綜合考慮流速、流型、質(zhì)量傳遞同沖蝕之間的相互作用，率先結(jié)合流體動(dòng)力學(xué)來(lái)研究管道沖蝕失效規(guī)律，提出實(shí)現(xiàn)沖蝕預(yù)測(cè)的可能性。Nesic S等[2]采用數(shù)值模擬方法，研究電極表面近壁面處流體作用力對(duì)沖蝕的影響，并取得了初步成果。季楚凌[3]以輸送稠油的彎管為研究對(duì)象，借助Fluent與Ansys軟件對(duì)彎管進(jìn)行熱流固耦合模擬計(jì)算。陳佳等[4]針對(duì)三通管建立了沖刷腐蝕的數(shù)學(xué)模型，運(yùn)用SIMPLE算法求解方程組。劉勇峰等[5]將計(jì)算流體的方法引入到彎管沖刷腐蝕的研究中，建立了彎管沖刷腐蝕的數(shù)學(xué)模型。并根據(jù)氣田集輸管道中的實(shí)際運(yùn)行參數(shù)，模擬三種工況下管道彎管中流體的運(yùn)動(dòng)。陳思維等[6]以濕氣集輸條件下管流流體動(dòng)力學(xué)數(shù)值模型為基礎(chǔ)，通過(guò)數(shù)值模擬并結(jié)合現(xiàn)場(chǎng)檢測(cè)結(jié)果，開(kāi)展了氣液兩相沖蝕規(guī)律、沖蝕失效風(fēng)險(xiǎn)預(yù)測(cè)研究。高萬(wàn)夫等[7]研究了輸送管道中氣固兩相流中微粒的直徑、濃度和材料種類等因素對(duì)管道彎頭磨損量以及磨損率的影響。趙燕輝等[8]借助CFD軟件對(duì)油氣混輸流體流經(jīng)T型管進(jìn)行流場(chǎng)和應(yīng)力分析，探討了流體組成、流體性質(zhì)和流動(dòng)參數(shù)等對(duì)沖刷腐蝕的影響。

國(guó)內(nèi)外學(xué)者對(duì)局部管件內(nèi)的流場(chǎng)分布進(jìn)行了大量研究，但是專門針對(duì)石化管道沖蝕磨損進(jìn)行數(shù)值分析的研究還比較少。通過(guò)研究可以得到該石化管道的沖蝕磨損最嚴(yán)重位置，對(duì)指導(dǎo)石化管道的檢測(cè)有重要幫助。

1　數(shù)值模擬分析

1.1理論模型

管道內(nèi)氣體中含液滴的流動(dòng)屬于典型的氣—液兩相流，其連續(xù)方程如下：

(1)

湍流模型采用適合流動(dòng)類型比較廣泛的RNG k-ε模型，k-ε方程如下：

Gb-ρε-YM

(2)

(3)

式中：GK是由于平均速度梯度引起的湍動(dòng)能產(chǎn)生；Gb是由于浮力影響引起的湍動(dòng)能產(chǎn)生；YM為可壓縮湍動(dòng)能膨脹對(duì)總的耗散率的影響；C1ε、C3ε為經(jīng)驗(yàn)常數(shù)，F(xiàn)LUENT中默認(rèn)C1ε=1.44、C3ε=0.09；σk、σε分別為湍動(dòng)能和湍動(dòng)能耗散率對(duì)應(yīng)的普朗特?cái)?shù)，F(xiàn)LUENT中默認(rèn)值為σk=1.0、σε=1.3。

1.2計(jì)算模型

根據(jù)某化工廠管道系統(tǒng)發(fā)生沖蝕的管線，建立如圖 1 所示的物理模型。包含一個(gè)入口、兩個(gè)盲三通、四個(gè)彎管、一個(gè)異徑管和一個(gè)出口，管道入口直徑為 100 mm、出口直徑80 mm、管道總長(zhǎng)6 700 mm，異徑管入口直徑100 mm、出口直徑80 mm、高300 mm。 輸氣管線的布局如圖1所示。

圖1　管線示意圖

1.3網(wǎng)格劃分

采用專用前處理軟件GAMBIT生成網(wǎng)格。劃分管線的網(wǎng)格時(shí)，忽視閥門、儀表等對(duì)管線的影響，采用古錢幣畫法對(duì)整條管線生成結(jié)構(gòu)化網(wǎng)格，網(wǎng)格數(shù)目為174 289。

1.4邊界條件

入口邊界條件：在圓形入口截面上給定法向時(shí)均速度v=5 m/s、湍流強(qiáng)度為5%、水力學(xué)直徑為100 mm。出口邊界條件是outflow，即認(rèn)為出口時(shí)的湍流是充分發(fā)展的。壁面邊界條件采用固體邊界(wall邊界)。

1.5介質(zhì)特性

采用單一柴油介質(zhì)對(duì)管線進(jìn)行沖蝕研究，柴油的密度為850 kg/m3，黏度為0.004 845 Pa·s。

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

通過(guò)fluent軟件對(duì)管線的穩(wěn)態(tài)流域進(jìn)行了數(shù)值模擬，模擬結(jié)果如下。

2.1模擬結(jié)果

圖2　管線的整體剪切力云圖

由圖2可以看出，對(duì)于管線磨損最嚴(yán)重的區(qū)域在所標(biāo)示的彎管與盲三通之間，其中靠近彎管內(nèi)側(cè)中心位置的壁面剪切應(yīng)力值最大，所以此處磨損最為嚴(yán)重。此外，彎管外拱下游及內(nèi)拱壁面所受切應(yīng)力相對(duì)較大，內(nèi)拱壁面處柴油速度和沖擊角度都相對(duì)較大，導(dǎo)致內(nèi)拱壁面處的切應(yīng)力值相對(duì)較大。而彎管外拱下游近壁面處柴油的速度較大，所以在該處彎管壁面所受剪切力較大。

圖3　最大剪切力發(fā)生區(qū)域放大顯示圖

由圖3可以看出，管線最大壁面剪切應(yīng)力為377 Pa，最大壁面剪切應(yīng)力出現(xiàn)的位置在盲三通入口處的壁面左側(cè)中間位置。

2.2磨損強(qiáng)度排序

通過(guò)對(duì)管線壁面剪切應(yīng)力的模擬能夠更清楚直觀的分析出管線發(fā)生磨損的位置和程度，可以對(duì)管線較長(zhǎng)且磨損狀況較為復(fù)雜的管線進(jìn)行預(yù)測(cè)。通過(guò)對(duì)管線的模擬可以得到管線的壁面剪切應(yīng)力云圖，即獲得管線的整體磨損情況。對(duì)該石化管線的磨損情況排序結(jié)果：4內(nèi)側(cè)>4外側(cè)>6>2>5>3=1，如圖4所示。

圖4　管件沖蝕結(jié)果大小排序

3　結(jié)論

某煉廠該柴油管線沖蝕最嚴(yán)重的位置在彎管與盲三通之間靠近彎管內(nèi)拱壁面中心位置(編號(hào)4)，此處壁面剪切應(yīng)力值最大，為377 Pa。對(duì)該輸氣管線的磨損情況排序結(jié)果為4內(nèi)側(cè)>4外側(cè)>6>2>5>3=1，該結(jié)果將指導(dǎo)該管線的測(cè)厚點(diǎn)布置和檢測(cè)。

[1]Blatt W,Heitz E.Hydromechanical measurements for erosion corrosion[J].National Association of Corrosion Engineers, 1993: 543-546.

[2]Nesic S,Postlethwaite J.The Canadian of Chem Eng[J].1991,69:698-702.

[3]季楚凌.稠油管道90°彎管流場(chǎng)及應(yīng)力分析[J].當(dāng)代化工,2015,44(2):401-404.

[4]陳佳,劉勇峰. 三通管沖刷腐蝕數(shù)值計(jì)算[J].當(dāng)代化工,2013,42(1):76-78.

[5]劉勇峰,吳明,趙玲,等.凝析氣田集輸管道彎管沖刷腐蝕數(shù)值計(jì)算[J].腐蝕與防護(hù),2012,33(2):132-135.

[6]陳思維,覃明友,劉德緒,等.高含硫氣田濕氣集輸管道沖蝕風(fēng)險(xiǎn)預(yù)測(cè)研究[J].天然氣與石油,2015,33(1):80-83，12.

[7]高萬(wàn)夫,鄭雁軍,崔立山,等.管道彎頭磨損特性的研究[J].石油化工高等學(xué)校學(xué)報(bào),2003,16(4):56-60.

[8]趙燕輝,張濤,張義貴,等.集輸管道T型管內(nèi)沖刷腐蝕數(shù)值模擬[J].當(dāng)代化工,2014,43(11):2457-2459.

Numerical Simulation Analysis on Erosion of Petrochemical Pipeline
YAO Cheng1,2， LI Jiepu1,2， LI Xiang2,1
(1.China University of Petroleum-Beijing ， Beijing102249 , China ； 2.China Special Equipment Inspection and Research Institute ， Beijing100029 , China)

The flow field of diesel oil flowing throughout the pipeline is analyzed.The pipeline is simulated using Fluent software.The results show that the most serious erosion of the pipeline is close to the center of the wall between the elbow and the tee,which is the biggest shear stress.The results can guide the operation of the petrochemical pipeline.The reliability and security for the operation of the petrochemical pipeline system are effectively improved.Inspection cost is saved for petrochemical pipeline system due to the investigation.

pipeline ; erosion ; wear ; numerical simulation

1003-3467(2016)07-0019-03

2016-05-10

幺成(1992-)，男，碩士，從事化工機(jī)械研究工作，電話：18510864998。

TE832，TQ050.2

A