劉 波,王樹立,楊 燕,李靖博,汪敏慧,馮旻祎

(常州大學(xué) 江蘇省油氣儲運技術(shù)重點實驗室，常州 213016)

涂層孔隙率對鋼質(zhì)套管穿越段管道陰保電位影響規(guī)律的數(shù)值模擬

劉 波,王樹立,楊 燕,李靖博,汪敏慧,馮旻祎

(常州大學(xué) 江蘇省油氣儲運技術(shù)重點實驗室，常州 213016)

采用COMSOL Multiphysics軟件，對管道極化電位和管地電位進行模擬計算，從套管有涂層和無涂層的兩個方面，研究了涂層孔隙率對套管處埋地管道的陰保電位影響規(guī)律。結(jié)果表明：套管內(nèi)管道涂層質(zhì)量越差，管道所需的極化電流就越多，管道的陰保效果就越差；加涂層套管在環(huán)形空間內(nèi)的電位比無涂層套管的更正，但在地表處卻稍微偏負，所以僅僅依靠管地電位來判斷管道的極化電位是有誤差的；當(dāng)套管無涂層時，只要環(huán)形空間內(nèi)有電解質(zhì)，其對管道的極化電流無屏蔽作用，所以利用密間隔電位測量(CIPS)、直流電壓梯度法(DCVG)和交流電壓梯度法(ACVG)等地面檢測方法可以判斷出套管內(nèi)管道防腐層的質(zhì)量和陰極保護的等級狀況。

涂層孔隙率；埋地管道；陰極保護；模擬

在我國城市化急劇擴張和經(jīng)濟快速發(fā)展的背景下，管道與越來越多的人類公共走廊交叉穿越，按照現(xiàn)行的油氣管道穿越工程設(shè)計規(guī)范，穿越鐵路或二級及二級以上公路時，應(yīng)在套管或涵洞內(nèi)敷設(shè)穿越管段[1]。由于套管與管道之間無電解質(zhì)，所以陰極保護電流和檢測設(shè)備的發(fā)射信號無法抵達和反饋，使得該處的管道成為腐蝕的易發(fā)點和外腐蝕直接評價技術(shù)(ECDA)檢測的盲點[2]。因此，管道安全評估行業(yè)急需利用數(shù)值模擬技術(shù)預(yù)測管道的腐蝕電位，從而決定該處的陰極保護等級。COMSOL Multiphysics是一款大型的高級數(shù)值仿真軟件，在各個科學(xué)領(lǐng)域和工程計算中都有廣泛的應(yīng)用，其中在金屬電化學(xué)方面的研究頗多[3-5]。

目前，油氣管道的防腐蝕系統(tǒng)絕大多數(shù)采用陰極保護和防腐涂層聯(lián)合保護，陰保電流通過涂層材料的空隙進入管道從而建立起陰極保護。對于套管穿越段，管道涂層采用加強級防護，而套管外表面是否加涂層一般視現(xiàn)場情況而定。管道保護電位的大小及其分布狀況與多種因素有關(guān)，其中與管道涂層質(zhì)量有著密切的關(guān)系[6-7]。涂層質(zhì)量越好，管道所需的陰保電流就越小，陽極使用壽命就越長，但是隨著管道運行時間的推移，涂層的質(zhì)量會有所下降，此時會影響到管道陰極保護等級狀況[8-10]。本工作結(jié)合一條500 m長的管道，運用COMSOL Multiphysics軟件對管道和套管的涂層孔隙率與管道陰保電位(極化電位)以及管地電位之間的影響規(guī)律進行計算，為套管穿越段管道保護方法與檢測技術(shù)提供參考。

1 控制方程與邊界條件

對于均一、各向同性的體系而言，所建立的幾何模型在溶液中的靜電控制方程滿足拉普拉斯方程，如式(1)所示。

式中：φi為幾何中ith域的電解質(zhì)靜電勢(幾何體中有兩個域：套管內(nèi)、外的土壤)。

管道表面的總電流(Ip)與流經(jīng)管道和電解質(zhì)界面的電流密度有關(guān)，如式(2)所示。

式中：Dp為主管道的外徑；inetpi為ith域內(nèi)流過管道與電解質(zhì)界面的凈電流密度。

管道金屬壁面沿軸向的金屬電位Ψp，遵循歐姆定律，如式(3)所示。

式中：ρFe為管道的電阻率；t為管道壁厚。

把式(3)代入式(2)中得到管道中金屬電位的控制方程，如式(4)所示。

以式(1)和式(4)為模型獲取的腐蝕電位即為所需要解決的控制方程。

根據(jù)控制方程所需的邊界條件，對模型進行物理定義，包括管道和套管表面的電流密度、各電解質(zhì)及土壤的電阻率、陰陽極極化斜率、與管道平行方向的地表面設(shè)置無限域等，土壤上表面和四周均采用絕緣邊界條件，二者的電流設(shè)為零，邊界條件可定義為φi=0。

2 模型建立及驗證

2.1 模型建立與參數(shù)設(shè)置

根據(jù)現(xiàn)場套管處管道的埋設(shè)特點，利用COMSOL Multiphysics軟件對管道表面的電位進行模擬。第一步，選擇三維、一次電流分布、穩(wěn)態(tài)的計算環(huán)境。第二步，利用Geometry建立分析模型，先建立套管、管道和陽極模型，而后建立土壤環(huán)境，具體幾何參數(shù)見圖1。第三步，對模型進行物理定義，包括管道和套管涂層孔隙率、涂層和土壤的電阻率、陽極電位以及邊界無限域。第四步，對已生成的模型采用自由剖分四面體網(wǎng)格劃分，劃分后的模型見圖2。最后對計算結(jié)果進行后處理，得出設(shè)定位置的結(jié)果值。

涂層的孔隙率是指整個有涂層的管道上，裸露金屬的表面積占總管段表面積的百分比，例如孔隙率為10-4相當(dāng)于1 m2的涂層表面有1 cm2的裸露面積。套管和管道涂層的孔隙率分別設(shè)為10-5，10-4，10-3，土壤和涂層的電阻分別為150Ω和1.5×10-6Ω，陽極電位為-1.5 V，交換電流密度為0.1 mA/m2。

圖1 幾何模型示意圖(單位：m)Fig. 1 Sketch of the model geometry

圖2 管道及土壤網(wǎng)格劃分Fig. 2 Grid generation of pipeline and soil

2.2 模型驗證

通過試驗測套管內(nèi)管道極化電位，對此模型進行驗證，試驗裝置如圖3所示。管道長8 m，距管道左端2 m處設(shè)有套管，套管全長4 m，管道下方設(shè)有犧牲陽極，每0.5 m有一電位測試點，模擬結(jié)果與測試結(jié)果如圖4。

由圖4可知，該模型模擬電位的最大誤差為0.246%，平均誤差為0.087%，在允許的范圍之內(nèi)。由于模型是對實際情況的簡化，而且本身計算過程是迭代過程，所以會產(chǎn)生一定的誤差，但誤差在允許范圍內(nèi)。驗證試驗結(jié)果表明采用COMSOL Multiphysics所建立的模型是可用于研究套管內(nèi)管道的電位分布。

圖3 試驗裝置示意圖Fig. 3 Schematic diagram of experimental equipment

圖4 模擬電位與實測電位對比Fig. 4 Comparison between simulated potential and measured potential

3 結(jié)果與討論

3.1 管道涂層孔隙率的影響

對埋地管道的現(xiàn)場檢測，通常利用參比電極在管道上方沿管道鋪設(shè)方向獲取管道極化電位，但是由于管道與土壤間存在IR降，使得檢測數(shù)據(jù)并不等于真實電位。圖5為套管無涂層時，管道涂層孔隙率對管地電位的影響。結(jié)果表明，管地電位和管道極化電位有很明顯的差別，而且隨著涂層孔隙率的減小，二者差別增大。涂層孔隙率越小，管道與地表之間的電阻

圖5 管道涂層孔隙率對管地電位的影響(套管無涂層)Fig. 5 Effect of coating porosity of pipeline on pipe-to-soil potential (no coating on casing)

越大，IR降則變化越大，所以二者差別越大。雖然涂層和土壤都有一定的電阻，但可以看出，涂層孔隙率的變化是造成電位降的主要原因。

從圖5中還可以看出，管道極化電位和管地電位在套管段沒有任何變化，表明無涂層的鋼質(zhì)套管對環(huán)形空間內(nèi)管道的陰極保護沒有影響。相比土壤和涂層的電阻，金屬套管的極化電阻很小，而且被分散到整個表面，所以外部陰保電流可以通過套管流入管道，同樣對該特殊段起到保護作用。

圖6為套管無涂層時，管道涂層孔隙率對管道極化電位的影響。從圖6中可以看出，隨著管道涂層孔隙率的增大，管道極化電位越來越正，這主要是因為隨著涂層孔隙率的增大，管道需要更多的陰保電流，但是陽極所能提供的保護電位以及其表面積無法改變，所以管道的極化電位就朝著正方向移動。該模擬結(jié)果與現(xiàn)場實際相吻合，通常隨著運行時間的推移，其管道涂層的電阻越來越小，即孔隙率越來越大，所需的陰保電流就會增加，致使管道的極化電位達不到保護的標(biāo)準(zhǔn)。當(dāng)管道涂層孔隙率由1×10-5增大到5×10-5時，管道的保護電位就正移了約100 mV，低于標(biāo)準(zhǔn)規(guī)定的-850 mV，達不到陰極保護的效果。

圖6 管道涂層孔隙率對管道極化電位的影響(套管無涂層)Fig. 6 Effect of coating porosity of pipeline on polarized potential of pipeline (no coating on casing)

由圖6還可見，管道在穿越套管段時，其電位并沒有改變，同樣表明無涂層的鋼質(zhì)套管對環(huán)形空間內(nèi)管道的陰極保護電位無影響，而管道整體涂層孔隙率的變化對管道極化電位有很明顯的影響。

圖7為套管涂層孔隙率10-5時，管道涂層孔隙率對管道極化電位的影響。從圖7中可以看出，雖然套管涂層孔隙率保持不變，但隨著環(huán)形空間內(nèi)管道涂層孔隙率的不斷增大，管道的整體電位朝著正方向偏移，表明環(huán)形空間內(nèi)管道涂層的質(zhì)量不僅對其內(nèi)部電位有影響，同時也影響整個管道的極化電位。當(dāng)環(huán)形空間內(nèi)的管道無涂層或涂層孔隙率相對套管涂層更大時，外部陰保電流更容易通過此處流回管道，同時需要更大陰保電流來達到陰極保護的目的。由于外部陽極所提供的陰保電流是一定的，所以造成該特殊段管道的陰保電位下降，大部分電流通過涂層孔隙率更大的地方流回，也致使套管外管道整體的陰保電位下降。

圖7 管道涂層孔隙率對管道極化電位的影響(套管涂層孔隙率10-5)Fig. 7 Effect of coating porosity of pipeline on polarized potential of pipeline (casing coating porosity of 10-5)

管道涂層孔隙率的改變對套管內(nèi)管道極化電位有很明顯的影響，在現(xiàn)場，如果環(huán)形空間內(nèi)進入水或土壤等電解質(zhì)，此時可以利用直流電壓梯度法(DCVG)、密間隔電位測量(CIPS)和交流電壓梯度法(ACVG)等地面檢測方法判斷管道涂層是否存在漏點和陰極保護的等級。

3.2 套管涂層孔隙率的影響

套管涂層孔隙率對管道極化電位的影響如圖8所示。由圖8可見，與套管無涂層時相比，套管有涂層時，管道在套管外的電位輕微偏負，而且隨著套管涂層孔隙率的減小，管道極化電位偏移越大，而套管內(nèi)管道極化電位更偏向正方向。外部陰保電流是通過套管流入管道的，當(dāng)套管有涂層時，限制了部分陰保電流流入管道，所以在某種程度上，可以認為套管屏蔽了陰保電流，當(dāng)套管涂層孔隙率越小時，流入環(huán)形空間內(nèi)管道的電流就越少，該處管道的極化效果就越差，所以套管內(nèi)管道極化電位更正。

圖8 套管涂層孔隙率對管道極化電位的影響(管道涂層孔隙率為10-5)Fig. 8 Effect of coating property of casing on polarized potential of pipeline (pipeline coating property of 10-5)

圖9為套管涂層孔隙率對管地電位的影響。從圖9中可以看出，套管涂層的孔隙率越小，管地電位就越負。當(dāng)套管涂層孔隙率在10-5時，管地電位在套管處稍微偏負，但數(shù)值很小，所以在地表檢測套管內(nèi)管道極化電位時，該誤差可以忽略不計。當(dāng)套管有涂層時，管地電位在套管段基本無明顯變化，這是由于套管屏蔽了大部分從套管內(nèi)部反饋到外部的電流信號，所以才使得地表處的電位沒有明顯變化。這表明在工程檢測現(xiàn)場，當(dāng)套管有涂層時，僅僅通過管地電位的大小變化是很難判斷出套管內(nèi)部管道的陰保狀況和腐蝕環(huán)境，需要做進一步的調(diào)查研究。

圖9 套管涂層孔隙率對管地電位的影響(管道涂層孔隙率為10-5)Fig. 9 Effect of coating property of casing on pipe-to-soil potential (pipeline coating property of 10-5)

4 結(jié)論

(1) 利用COMSOL Multiphysics軟件建立的陰極保護模型，可以模擬研究套管內(nèi)管道的電位變化，該模型可以預(yù)測一些現(xiàn)場無法檢測的管道極化電位，能為管道完整性評估提供參考依據(jù)。

(2) 管道涂層的質(zhì)量對管道的腐蝕防護影響很大，環(huán)形空間內(nèi)涂層質(zhì)量越差，管道得到的陰極保護就越少，陽極損耗的陰保電流就越多，這與現(xiàn)場檢測情況相符合，而且該種情況在套管內(nèi)尤為明顯。

(3) 加涂層套管在環(huán)形空間內(nèi)的電位比無涂層的更正，但在地表處卻稍微偏負。此情況可以說明，在不了解套管環(huán)形空間內(nèi)管道的真正極化電位時，利用管地電位來判斷管道的陰保電位是有誤差的。

(4) 計算結(jié)果表明，當(dāng)套管無涂層時，只要環(huán)形空間內(nèi)有電解質(zhì)，其對管道的陰保電流無屏蔽作用，所以利用CIPS、DCVG和ACVG等地面檢測方法可以判斷出套管內(nèi)管道防腐層的質(zhì)量和陰極保護的等級狀況。

[1] GB 50423-2007 油氣輸送管道穿越工程設(shè)計規(guī)范[S]. 北京:中國計劃出版社,2008.

[2] 張珂,史國福,寧尚鋒,等. 鋼質(zhì)套管對埋地管道陰極保護的影響[J]. 腐蝕與防護,2007,28(11):580-582.

[3] 趙琳. 金屬套管腐蝕檢測方法與技術(shù)研究[D]. 西安:西安石油大學(xué),2013.

[4] 王雷,董麗娜. 基于COMSOL Multiphysics的雜散電流腐蝕仿真分析[J]. 新技術(shù)新工藝,2014(1):22-24.

[5] DESHPANDE K B. Numerical modeling of micro-galvanic corrosion[J]. Electrochimica Acta,2011,56(4):1737-1745.

[6] 李彬,袁宗明,秦燕,等. 在役埋地管線涂層質(zhì)量評價方法[J]. 石油化工腐蝕與防護,2004,21(5):52-54.

[7] 段行知,胡劍,杜艷霞,等. 淺析不同涂層管線聯(lián)合陰極保護時電位分布的影響因素[C]//2013年中國油氣田腐蝕與防護技術(shù)科技創(chuàng)新大會論文集.[出版地不詳]：[出版者不詳]，2013:416-420.

[8] 張盈盈,韓文禮,李愛貴，等. 外防腐蝕涂層對管道陰極保護電流屏蔽的研究進展[J]. 腐蝕與防護,2010,31(10):787-789.

[9] 郭娟麗,竇宏強,廖煜炤，等. 鋼制管道外防腐蝕涂料的應(yīng)用進展[J]. 腐蝕與防護,2011,32(12):980-983.

[10] 唐誼平,李建新,黃子陽，等. 長輸管線外防護涂層技術(shù)的現(xiàn)狀與發(fā)展探討[J]. 腐蝕與防護,2009,30(12):860-864.

Numerical Simulation of The lnfluence Rule of Coating Porosity on Cathodic Potential of Pipeline in Steel Casing

LIU Bo, WANG Shu-li, YANG Yan, LI Jing-bo, WANG Min-hui, FENG Min-yi

(Jiangsu Key Laboratory of Oil-Gas Storage and Transportation Technology, Changzhou Universtiy, Changzhou 213016, China)

Using COMSOL Multiphysics software, the influence rule of coating porosity on polarized potential of buried pipelines and pipe-to-soil potential was investigated from two aspects of casing with coating and without coating. Results show that the worse the pipeline coating quality is, the more polarization current is required, so the results of cathodic protection is more undesirable. The potential in coating casing annular space is more positive than that in uncoated casing, but more negative on ground surface, so it is inaccurate to determine the polarization potential of pipeline just rely on the pipe-to-soil potential. When the casing has no coating, there is no shield influence to polarization current of pipe as long as electrolyte is in the annular space. So CIPS, DCVG, ACVG can be used to determine the pipe coating quality and cathodic protection level.

coating porosity; buried pipe; cathodic protection; simulation

2015-00-00

江蘇省科技廳資助項目(BY2014037-33)

王樹立(1957-)，教授，博士，從事油氣管道腐蝕與防護研究，0519-83293890，wsl@cczu.edu.cn

10.11973/fsyfh-201612003

TE988

A

1005-748X(2016)12-0961-05