導管架平臺外加電流陰極保護數(shù)值模擬計算研究

劉福國,武素茹

(1.中國海洋石油總公司海洋石油工程股份有限公司,天津 300451;2.天津出入境檢驗檢疫局化礦金屬材料檢測中心,天津 300456)

導管架平臺外加電流陰極保護數(shù)值模擬計算研究

劉福國1,武素茹2

(1.中國海洋石油總公司海洋石油工程股份有限公司,天津 300451;2.天津出入境檢驗檢疫局化礦金屬材料檢測中心,天津 300456)

數(shù)值模擬技術是求解陰極保護電位模型的有效方法,數(shù)值模擬計算可以有效地預測陰極保護被保護體表面電位分布,并以此驗證陽極布置效果。該文采用大型防腐蝕數(shù)值模擬分析軟件 Beasy對小型導管架式平臺外加電流陰極保護進行數(shù)值模擬分析研究,并進行物理模型試驗。結(jié)果表明數(shù)值模擬計算能夠模擬分析陽極布置以及被保護體表面電位分布,輔助陽極的位置對導管架表面的電位分布有一定的影響,輔助陽極位置距離導管架越近,導管架表面電位分布越不均勻。采用物理模型試驗對數(shù)值模擬計算進行的驗證表明,物理模型試驗與數(shù)值模擬計算結(jié)果吻合較好。

外加電流 陰極保護 數(shù)值模擬 導管架

外加電流陰極保護是基于電化學腐蝕原理的一種防腐蝕手段,也是國內(nèi)外公認的防止鋼質(zhì)構(gòu)筑物腐蝕的經(jīng)濟有效的措施[1]。在陰極保護工程中,電位是監(jiān)視和控制陰極保護效果的一項重要指標,因此掌握被保護體表面的電位分布是非常重要的。在傳統(tǒng)的陰極保護工程中,大多采用實際測量或經(jīng)驗估計的方法來獲得被保護體表面的電位分布[2]。然而,對于某些被保護結(jié)構(gòu),如海底管道、海洋平臺、深埋的鋼樁或鋼管、大罐的罐底等,實地測量技術難度很大或費用昂貴,對于新項目更不可能事先在現(xiàn)場測試,而經(jīng)驗公式又無法對復雜結(jié)構(gòu)表面的電位分布進行準確的預測[3],可見傳統(tǒng)的方法已難以滿足越來越高的安全可靠性和經(jīng)濟性的要求。隨著電化學和計算機技術的發(fā)展,人們嘗試采用數(shù)值計算方法來獲取被保護體表面的電位和電流分布狀況,這已成為陰極保護領域中十分活躍的一個方面[4]。

1 數(shù)值模擬計算數(shù)學模型[5]

在陰極保護中通常研究穩(wěn)態(tài)電位場。對于已經(jīng)達到穩(wěn)態(tài)的陰極保護體系,其中的電位分布可歸結(jié)為橢圓型方程的混合邊值問題。均勻?qū)щ娊橘|(zhì)中穩(wěn)態(tài)的陰極保護體系的電位分布為Laplace方程混合邊值問題,電場 (V)內(nèi)每一點的電位 (E)應符合下述邊界條件:

S1:E=E0(第一類邊界條件,電位為常數(shù)的電極表面 S1)

式中:E——電位,V;

σ——腐蝕介質(zhì)電導率,S/m;

i(E)——金屬的極化行為作為邊界條件;n——為邊界內(nèi)法線方向。

運用邊界元法將Laplace方程轉(zhuǎn)化成如下形式的邊界積分方程:

對整個邊界來說,將各單元疊加,可以用矩陣的形式來表示:

式中:H——電位矩陣;

G——電流密度矩陣。

將邊界條件帶入,用迭代法求解上述方程,可以得出各點電位和電流密度。

2 外加電流陰極保護數(shù)值模擬

2.1 試驗方法

該文采用Beasy軟件進行導管架外加電流陰極保護數(shù)值模擬計算。Beasy軟件是英國 CM BEASY公司開發(fā)的大型防腐蝕軟件,在石油化工、海洋工程等領域應用較多,此軟件是基于邊界元程序開發(fā)的,邊界元程序軟件能夠在不影響計算精度的前提下,降低計算維數(shù),減少計算量,特別適用于海洋工程大型結(jié)構(gòu)物的防腐蝕數(shù)值模擬計算。首先采用邊界元軟件 Beasy對導管架 (圖 1)建模 (圖 2)、網(wǎng)格剖分及外加電流陰極保護數(shù)值模擬。外加電流輔助陽極平均分布到 4個斜面,并遠離導管架固定于海水底面,Ag/AgCl參比電極放置到導管架上端中間位置,并始終浸沒于海水中。根據(jù)數(shù)值模擬輔助陽極大小、距離導管架遠近及電位值進行物理模型試驗。試驗水池尺寸 2.4 m×2.4 m×1.6 m,恒電位儀選用T D3691型(天津市中環(huán)電子儀器公司),電位、電流測量儀選用福泰克 9205A數(shù)字萬用表,輔助陽極選用管狀鈦基金屬氧化物陽極,試驗溫度為 (29±1)℃。

圖 2 采用 Beasy軟件對導管架建模Fig.2 Platfor m model by usingBeasy software

輔助陽極距離導管架遠近也影響陰極保護效果,輔助陽極距離導管架選擇 0.2,0.4和 0.6 m分別進行導管架外加電流陰極保護數(shù)值模擬計算,見圖 3～圖 5。輔助陽極電位設置為 -1.60 V。

圖 3 陽極距模型 0.2 m數(shù)值模擬電位分布Fig.3 Numerical simulation potential distribution of 0.2m be tween anode and model

圖 4 陽極距模型 0.4 m數(shù)值模擬電位分布Fig.4 Numerical simulation potential distribution of 0.4m be tween anode and model

圖 5 陽極距模型 0.6 m數(shù)值模擬電位分布Fig.5 Numerical s imulation potential distribution of 0.6m be tween anode and model

2.2 陰極保護數(shù)值模擬結(jié)果與討論

由導管架外加電流陰極保護數(shù)值模擬結(jié)果可知,輔助陽極的位置對導管架表面的電位分布有一定的影響,輔助陽極位置距離導管架越近,導管架表面電位分布越不均勻。但是都在挪威船級社(DNV)所規(guī)定的范圍內(nèi)[6]。采用數(shù)值模擬計算可以有效預測被保護體表面電位的分布情況,能給腐蝕工程師提供參考。采用物理模型試驗對數(shù)值模擬計算進行驗證。

2.3 陰極保護物理模型試驗結(jié)果與討論

根據(jù)數(shù)值模擬計算結(jié)果,選擇輔助陽極距離導管架 0.4 m時進行物理模型試驗 (圖 6),測量點見圖 7,待導管架表面電位穩(wěn)定后進行測量,表面電位分布見表 1。導管架外側(cè)表面測量點位 A,對應的內(nèi)側(cè)測量點位為A?！?/p>

圖 6 導管架陰極保護物理模型試驗Fig.6 Model experiment of platfor m cathodic protection

由物理模型試驗測量數(shù)據(jù)可知,導管架表面電位分布 -973～-1 002 mV,由上到下電位變負,數(shù)值模擬計算結(jié)果表面電位分布在 -971～-998 mV,電位分布也是由上到下變負,兩者吻合較好。

表 1 陰極保護物理模型試驗電位測量值Table 1 Potential values of cathodic protection model experiment mV

3 結(jié) 論

外加電流陰極保護數(shù)值模擬計算可以有效地預測被保護體表面電位分布。數(shù)值模擬計算結(jié)果表明陽極的分布對被保護體表面電位的分布有一定的影響,陽極距離導管架 0.6 m時保護效果最好,物理模型試驗與數(shù)值模擬計算結(jié)果吻合較好。

[1] 于楠,梁成浩,吳建華,等 .基于縮比模型理論的船舶外加電流陰極保護系統(tǒng)水線下表面電位的研究[J].腐蝕與防護,2006,5(28):369-372.

[2] Degiorgi V G,W immer S A.Geometric details and modeling accuracy requirements for shipboard impressed current cathodic protection system modeling[J].Engineering Analysis with Boundary Elements,2005(29):15-28.

[3] Douglas P R,Mark E O.A mathematical model for the cathodic protection of tank bottoms[J].Corrosion science,2005,47(3):849-868.

[4] 杜艷霞,張國忠,李健 .陰極保護電位分布的數(shù)值計算[J].中國腐蝕與防護學報,2008,28(1):53-58.

[5] 孫吉星.海洋結(jié)構(gòu)物陰極保護優(yōu)化模型及數(shù)值模擬計算[D].青島:中國海洋大學,2006.

[6] DETNORSKE VER ITAS(DNV).RP-B401-2005.Catho dic Protection Design[S].Nor way:DETNORSKE VER ITAS,2005.

Study on Numerical S imulation Calculation for Applied Current Cathodic protection of Platform

L iu Fuguo1,W u Suru2
(1.China Offshore O il Engineering Co.,L td.,Tianjin 300451;2.Chem icalsM inerals&M etallicM aterials Inspection Center of Entry-Exit Inspection&Quarantine Bureau,Tianjin300456)

Numerical s imulation calculation is an effective method to study the sacrificial anode cathodic protection,which can effectively predict the potential distribution on the surfaces of protected materials and verify the effectiveness of anode arrangement.The numerical s imulation calculation software-Beasy is used to study the applied current cathodic protection of platform,and the physicalmodel experiment is performed.The results show that numerical simulation calculation can analyze the anode distribution and potential distribution on the surface.The potential distribution is affected by anode position.The closer the anode is to the platform,the greater the ununiformity of potential distribution of platform.The data obtained from physicalmodel experiment fitwell with those from numerical simulation calculation.

applied current,cathodic protection,numerical simulation,platform

TG174.3

A

1007-015X(2011)06-0009-04

2011-02- 23;修改稿收到日期:2011-09-27。

劉福國,男,2008年畢業(yè)于中國海洋大學,博士學位,2010年中海研究總院博士后,現(xiàn)在中國海洋石油總公司海洋石油工程股份有限公司設計公司主要從事海洋結(jié)構(gòu)物防腐設計與科研工作。E-mail:liufg@mail.cooec.com.cn

(編輯 寇岱清)