,,2

(1. 中國石油大學(xué)(華東) 儲(chǔ)運(yùn)與建筑工程學(xué)院，青島 266580; 2. 中國石化青島液化天然氣有限責(zé)任公司,青島 266001)

專論

管道陰極保護(hù)數(shù)值計(jì)算方法的應(yīng)用進(jìn)展

寇杰1,尹雪明1,2

(1. 中國石油大學(xué)(華東) 儲(chǔ)運(yùn)與建筑工程學(xué)院，青島 266580; 2. 中國石化青島液化天然氣有限責(zé)任公司,青島 266001)

針對(duì)管道陰極保護(hù)數(shù)值模擬的應(yīng)用現(xiàn)狀，對(duì)陰極保護(hù)體系的控制方程以及三種常用邊界條件做了簡單的介紹。同時(shí)介紹了三種常用的數(shù)值計(jì)算方法以及國內(nèi)外研究學(xué)者對(duì)這三種數(shù)值計(jì)算方法的應(yīng)用情況，并對(duì)比分析了三種數(shù)值解法的優(yōu)缺點(diǎn)。最后指出了目前管道陰極保護(hù)數(shù)值模擬技術(shù)存在的不足,并提出了相應(yīng)的解決措施。

陰極保護(hù)；數(shù)值計(jì)算；有限差分法；有限元法；邊界元法

KOU Jie1, YIN Xueming1,2
(1. College of Pipeline and Civil Engineering, China University of Petroleum, Qingdao 266580, China;
2. Sinopec Qingdao LNG, Qingdao 266001, China)

目前，國內(nèi)外主要采用涂層與陰極保護(hù)技術(shù)相結(jié)合的方法來防止管道的腐蝕[1]。但管道所處環(huán)境復(fù)雜、管網(wǎng)中管道數(shù)量多、管道連接方式多樣，很難通過傳統(tǒng)的實(shí)地測量、室內(nèi)試驗(yàn)等方法來預(yù)測管道的保護(hù)效果，同時(shí)經(jīng)驗(yàn)公式在使用上的局限性也容易造成部分管道出現(xiàn)過保護(hù)或欠保護(hù)。因此，數(shù)值模擬技術(shù)逐漸發(fā)展起來并得到了廣泛的應(yīng)用，相比傳統(tǒng)的陰極保護(hù)設(shè)計(jì)方法，數(shù)值模擬技術(shù)可以更加準(zhǔn)確、高效地預(yù)測被保護(hù)管道的電位分布情況以及評(píng)價(jià)陰極保護(hù)的效果。

管道陰極保護(hù)系統(tǒng)的數(shù)值模擬即對(duì)被保護(hù)管道及其所處環(huán)境進(jìn)行合理假設(shè)，建立數(shù)學(xué)模型并進(jìn)行求解，通過計(jì)算機(jī)的數(shù)值計(jì)算和圖像輸出，定量描述管道電位分布情況，從而達(dá)到解決問題的目的[2]。常用的數(shù)值計(jì)算方法[3]有：有限差分法(FDM)、有限元法(FEM)以及邊界元法(BEM),這三種數(shù)值計(jì)算方法各有優(yōu)勢和不足。本工作基于各位專家學(xué)者對(duì)這三種常用的管道陰極保護(hù)數(shù)值模擬技術(shù)的研究，對(duì)數(shù)值求解方法的優(yōu)缺點(diǎn)以及管道陰極保護(hù)數(shù)值模擬技術(shù)存在的不足進(jìn)行了分析，并指出了相應(yīng)的解決措施,以期為指導(dǎo)今后的管道陰極保護(hù)數(shù)值計(jì)算研究提供借鑒和參考。

1 陰極保護(hù)體系的數(shù)學(xué)模型及邊界條件

陰極保護(hù)體系的數(shù)學(xué)模型有分布型模型和時(shí)變型模型兩種[4]，時(shí)變型模型主要用于研究陰極垢層的形成以及陰極極化隨時(shí)間的變化關(guān)系，模型建立時(shí)需要考慮溫度、土壤含水量等隨時(shí)間變化的環(huán)境參數(shù)。因?yàn)槿狈?duì)其規(guī)律的理解，所以時(shí)變型模型的研究難度較大，目前尚沒有合適的數(shù)值計(jì)算方法，因此對(duì)陰極保護(hù)體系的研究主要集中在分布型模型上。

分布型模型主要研究陰極保護(hù)體系達(dá)到穩(wěn)定狀態(tài)時(shí)管道上電位、電流分布及他們之間的關(guān)系[5]，可用靜態(tài)場理論來分析處理。在模型建立時(shí)需做如下假設(shè)[6]：環(huán)境介質(zhì)均勻、單一。當(dāng)電場為有源場時(shí)，穩(wěn)態(tài)陰極保護(hù)體系的電位分布可用泊松方程來描述[3],見式(1)。

式中:qs為場源；σ為區(qū)域內(nèi)介質(zhì)電導(dǎo)率。

電場為無源場時(shí)，方程變?yōu)槔绽狗匠蘙3]

微分方程必須具備相應(yīng)的邊界條件才能得到唯一解，作為數(shù)學(xué)模型的重要組成部分，邊界條件的選取對(duì)數(shù)值計(jì)算結(jié)果影響很大。因此，根據(jù)實(shí)際情況選擇合適的邊界條件是陰極保護(hù)數(shù)值計(jì)算的關(guān)鍵。陰極保護(hù)體系常用的邊界條件有以下三種[7]：

(1) 第一類邊界條件(Dirichlet條件),給定邊界上的電位值，即Φ(x,y)=ΦA(chǔ)。

(2) 第二類邊界條件(Neumann條件),給定邊界上的電流密度值，即?Φ(x,y)/?n(x,y)=qA。

(3) 第三類邊界條件(混合邊界條件),給定邊界上電位與電流密度之間的關(guān)系，即?Φ(x,y)/?n(x,y)=f(Φ)。

其中，Φ(x,y)為點(diǎn)(x,y)處的電位值，n(x,y)為在該點(diǎn)處表面的法線方向，qA為恒定的電流值，ΦA(chǔ)為恒定的電位值，f(Φ)為極化電流密度函數(shù)。

2 常用的數(shù)值計(jì)算方法

2.1 有限差分法(FDM)

20世紀(jì)60年代以來，有限差分法逐步應(yīng)用到電化學(xué)計(jì)算領(lǐng)域中，直到80年代，STROMMEN等[8]首次用有限差分法計(jì)算了陰極保護(hù)系統(tǒng)中被保護(hù)構(gòu)件表面的電位分布，并將其引入了陰極保護(hù)計(jì)算領(lǐng)域中。有限差分法可以近似求解偏微分方程的邊值問題，其基本原理[3]是用差分項(xiàng)[Φ(x+Δx)-Φ(x)]/Δx近似代替微分方程中的微分項(xiàng)dΦ/dx，從而求解微分方程。

有限差分法利用規(guī)則的網(wǎng)格對(duì)不規(guī)則的計(jì)算區(qū)域進(jìn)行劃分，其網(wǎng)格劃分模型如圖1所示，復(fù)雜區(qū)域的網(wǎng)格生成可能占總計(jì)算時(shí)間的大部分，網(wǎng)格的質(zhì)量對(duì)計(jì)算精度的影響很大，一般情況下，網(wǎng)格數(shù)越多，其得到的近似解精度越高，但當(dāng)網(wǎng)格足夠細(xì)密時(shí)，再進(jìn)一步加密網(wǎng)格對(duì)數(shù)值計(jì)算結(jié)果基本上沒有影響。

圖1 有限差分法網(wǎng)格劃分模型Fig. 1 Meshing model of finite difference method

張鳴鏑等[9]用有限差分法計(jì)算了海泥介質(zhì)中海底管道表面電位的分布情況及其隨時(shí)間的變化，同時(shí)對(duì)管道進(jìn)行了試驗(yàn)分析，試驗(yàn)結(jié)果與數(shù)值計(jì)算結(jié)果吻合良好，驗(yàn)證了有限差分法用于海底管道陰極保護(hù)電位分布計(jì)算的可行性，同時(shí)證明了有限差分法能夠用于海底管道陰極保護(hù)系統(tǒng)的監(jiān)測。

錢海軍等[10]采用有限差分法，使用“漏電”等效電路模型模擬管內(nèi)保護(hù)電位分布，并且還通過編寫程序?qū)Υ罂趶捷斔軆?nèi)陰極保護(hù)電位分布進(jìn)行了計(jì)算，結(jié)果與實(shí)測值十分吻合。

對(duì)于一維、二維計(jì)算場域,采用有限差分法所得計(jì)算結(jié)果可靠，但對(duì)于三維復(fù)雜結(jié)構(gòu)的陰極保護(hù)體系來說，因有限差分法采用折線來處理不規(guī)則的邊界，容易導(dǎo)致計(jì)算結(jié)果不收斂、計(jì)算精度降低等問題。因此，隨著被保護(hù)構(gòu)件復(fù)雜程度的增大以及人們對(duì)計(jì)算精度要求的提高，有限差分法在陰極保護(hù)體系上的應(yīng)用逐漸減少。

2.2 有限元法(FEM)

自20世紀(jì)70年代，有限元法逐漸被應(yīng)用到管道腐蝕防護(hù)的陰極保護(hù)設(shè)計(jì)中。有限元法是改進(jìn)后的有限差分法，是變分原理在差分方法中的應(yīng)用，其可用任意形狀的單元來劃分計(jì)算區(qū)域，從而便于復(fù)雜和彎曲邊界的處理，同時(shí)也提高了計(jì)算精度。有限元法不僅能夠得到管道表面上的電位、電流分布，也可以得到某個(gè)特定范圍內(nèi)的電位、電流分布情況。有限元法的網(wǎng)格劃分模型如圖2所示，有限元法可以針對(duì)每個(gè)不同的單元來設(shè)置其環(huán)境參數(shù)，從而更好地貼合實(shí)際。

圖2 有限元法網(wǎng)格劃分模型Fig. 2 Meshing model of finite element method

KENNELLEY等[11]在建模時(shí)考慮了涂層缺陷問題，得到了涂層缺陷對(duì)陰極保護(hù)電位分布的影響。PARSA等[12]指出土壤電阻率對(duì)電位分布有較大的影響，并得到了不同土壤電阻率下的最優(yōu)外加電流密度。KASPER等[13]首次指出穩(wěn)態(tài)溫度場和恒定電磁場有相似的自變量、因變量和邊界條件，并對(duì)笛卡爾坐標(biāo)系下熱傳導(dǎo)數(shù)學(xué)模型的控制方程進(jìn)行了類比離散，得到了海水中鋼棒陰極保護(hù)的電位分布。MARCASSOLI等[14]建立了海底管道陰極保護(hù)系統(tǒng)的二維模型，陽極采用恒電位控制，陰極假設(shè)有電絕緣涂層，并結(jié)合Butler-Volmer方程設(shè)置了鋼表面的涂層缺陷，最后根據(jù)有限元法計(jì)算了海底管道的電位分布情況以及涂層缺陷對(duì)管道電位分布的影響等。CHIN等[15]計(jì)算了二維穩(wěn)態(tài)涂層缺陷縫隙內(nèi)電化學(xué)環(huán)境改變時(shí)的電流分布。其將模型的控制方程在軸對(duì)稱柱坐標(biāo)系下離散為二階非線性差分方程，僅在徑向方向設(shè)定 ，從而減少了計(jì)算量。

在國內(nèi)，一些商業(yè)化的可用于有限元分析的軟件如COMSOL、ANSYS、FEMLAB、FEPG、ABAQUS等已廣泛用于管道陰極保護(hù)的計(jì)算中。

李佳奇[2]通過COMSOL Multiphysics軟件對(duì)埋地管道陰極保護(hù)電位分布進(jìn)行了系統(tǒng)計(jì)算，同時(shí)把傳統(tǒng)陰極保護(hù)設(shè)計(jì)方法與數(shù)值模擬方法進(jìn)行了對(duì)比分析，并提出將兩種方法結(jié)合對(duì)提高陰極保護(hù)的質(zhì)量、減少設(shè)計(jì)失誤率有重要意義。李成杰[16]通過COMSOL Multiphysics軟件建立了深海溫躍層環(huán)境中海管陰極保護(hù)的數(shù)學(xué)模型，得到了海管陰極保護(hù)電位分布情況，其建模時(shí)采用了動(dòng)態(tài)陰極邊界條件，考慮了溫度、與陽極之間的距離以及保護(hù)時(shí)間對(duì)管道電位分布的影響，解決了不同時(shí)間下鈣沉積層對(duì)陰極極化性能影響的問題。

宋成立等[17]將薄壁管等效成實(shí)心圓柱體，通過ANSYS軟件計(jì)算了采氣廠管道陰極保護(hù)體系的電位分布，系統(tǒng)評(píng)價(jià)了管線的運(yùn)行狀況并優(yōu)化了陰極保護(hù)體系。潘柳依等[18]通過ANSYS軟件研究了采氣廠中管道敷設(shè)垂直間距、涂層電阻率等對(duì)管道陰極保護(hù)電位分布的影響，成功優(yōu)化了欠保護(hù)管道的陰極保護(hù)系統(tǒng)，為管網(wǎng)的優(yōu)化設(shè)計(jì)提供了依據(jù)。徐勛[19]采用ANSYS軟件研究了外界參數(shù)對(duì)管道陰極保護(hù)系統(tǒng)電位分布的影響，提出計(jì)算域深度和寬度的增大會(huì)使整個(gè)管面的電位分布趨于均勻，而相對(duì)長度對(duì)管道表面電位影響不大等觀點(diǎn)。李丹丹等[20]利用ANSYS軟件分析了管道間陰極保護(hù)的干擾問題，建模時(shí)將兩條同溝敷設(shè)的管道等效成一條大直徑的管道，降低了計(jì)算量，劃分網(wǎng)格時(shí)在電場變化劇烈的區(qū)域?qū)W(wǎng)格模型進(jìn)行了局部加密處理，提高了計(jì)算的準(zhǔn)確度。

紀(jì)俊剛[21]利用FEMLAB軟件編制了計(jì)算程序，實(shí)現(xiàn)了海底管道陰極保護(hù)電位分布的計(jì)算。其采用了四面體單元進(jìn)行網(wǎng)格劃分，并設(shè)置了最大單元尺寸以及網(wǎng)格單元大小的增長因子來調(diào)整網(wǎng)格數(shù)量，有效控制了計(jì)算量、保證了求解精度。

孫吉星[22]在求解超長海管陰極保護(hù)電位分布的過程中,通過設(shè)置X、Y、Z軸比例因子、用正多面體來代替圓柱體進(jìn)行網(wǎng)格劃分、充分利用待解決問題的對(duì)稱性以及可共用數(shù)據(jù)等方法實(shí)現(xiàn)了管道的準(zhǔn)確剖分并有效控制了計(jì)算量。

目前，有限元法在管道陰極保護(hù)的設(shè)計(jì)上得到了普遍的應(yīng)用，有限元法適用性強(qiáng)，特別適合解決幾何和物理?xiàng)l件比較復(fù)雜的問題，便于編制標(biāo)準(zhǔn)化程序和工程應(yīng)用。但其仍存在模型計(jì)算域只能為有限域、計(jì)算量大、計(jì)算時(shí)間長等問題，在計(jì)算超長管道、結(jié)構(gòu)復(fù)雜管道等的陰極保護(hù)電位時(shí)與其他方法(如邊界元法)相結(jié)合會(huì)得到更精確的結(jié)果。

2.3 邊界元法(BEM)

20世紀(jì)80年代，邊界元法逐步在我國得到應(yīng)用，目前已經(jīng)成功用于埋地管道、海底管道以及近海石油平臺(tái)等領(lǐng)域，其在計(jì)算金屬表面電位分布上實(shí)現(xiàn)了優(yōu)化設(shè)計(jì)。邊界元法是在經(jīng)典積分方程的基礎(chǔ)上，吸收了有限元法的離散技術(shù)而發(fā)展起來的計(jì)算方法。其基本思想是用積分方程來求解微分方程。邊界元法的網(wǎng)格劃分模型如圖3所示，由于邊界元法只需要對(duì)邊界進(jìn)行離散化處理，因此可將計(jì)算域的維數(shù)降低一維，從而使得輸入數(shù)據(jù)量和代數(shù)方程組的未知量大大減少，有利于計(jì)算速率和計(jì)算精度的提高。

圖3 邊界元法網(wǎng)格劃分模型Fig. 3 Meshing model of boundary element method

TELLES等[23-25]采用邊界元法解決無限域中的陰極保護(hù)問題，成功得到了艦船陰極保護(hù)的電位分布。SANTIAGO等[26]等成功將與時(shí)間有關(guān)的極化曲線作為邊界條件融入到陰極保護(hù)數(shù)學(xué)模型中，為后人的研究提供了重要的依據(jù)。DIAZ等[27]優(yōu)化了不同形狀模型的輔助陽極的位置，從而準(zhǔn)確地設(shè)計(jì)了電流均勻分布的各模型的陰極保護(hù)系統(tǒng)。

ORAZEM等[28]采用線性單元對(duì)模型進(jìn)行網(wǎng)格劃分，得到涂層缺陷對(duì)管道陰極保護(hù)電位分布的影響。BRAEIL等[29]在建模時(shí)通過試樣連接管道來模擬涂層缺陷，并采用混合網(wǎng)格對(duì)陽極和管道進(jìn)行網(wǎng)格劃分，得到了帶有涂層缺陷埋地管道的電位分布，同時(shí)使用電流切斷法對(duì)試樣的電位降進(jìn)行了瞬時(shí)測量，并研究了試樣電位降與土壤電阻率的關(guān)系。

BRICHAU等[30]建立了埋地管網(wǎng)陰極保護(hù)的OKAPPI模型，在管道軸向應(yīng)用有限元法，在土壤無限大區(qū)域應(yīng)用邊界元法，將兩種計(jì)算方法結(jié)合使用得到了管道軸向土壤電位分布和軸向、徑向電流分布情況。MUHAREMOVIC等[31]采用有限元法和邊界元法研究了犧牲陽極系統(tǒng)的計(jì)算參數(shù)并對(duì)電位分布的影響因素進(jìn)行了探究。LACERDA等[32]通過雙向邊界元法，研究了細(xì)長管道中陰極保護(hù)電位和電流密度的分布情況。

伍欣[33]針對(duì)川氣東送管道防腐蝕層的破壞以及陰極保護(hù)系統(tǒng)不起作用的問題，采用BEASY CP軟件開展現(xiàn)場檢測分析，得到了管道相應(yīng)的保護(hù)電位準(zhǔn)則，為陰極保護(hù)的日常管理維護(hù)提供了依據(jù)。趙雷亮[34]采用BEASY軟件研究了陽極埋深、距管道距離以及埋設(shè)方式等對(duì)管道電位分布的影響。同時(shí)研究了站外干線陰極保護(hù)對(duì)站場陰極保護(hù)的干擾問題，并證明了當(dāng)輔助陽極移動(dòng)到一定距離時(shí)(200 m),其對(duì)站場管網(wǎng)的影響可以忽略。

劉立祺等[35]提出了一種三維管道邊界元模型，用軸線代替管道并將其離散成線單元，積分在每個(gè)線單元對(duì)應(yīng)的管道圓柱面上進(jìn)行，有效減少了離散難度和單元數(shù)量。同時(shí)采用一種新型的快速多級(jí)邊界元算法求解方程組，提高了計(jì)算速率。張東東等[36]也將快速多級(jí)算法應(yīng)用于邊界元法中，通過對(duì)二維管道表面電位分布的計(jì)算證明了此算法應(yīng)用于陰極保護(hù)領(lǐng)域的可行性，為應(yīng)用邊界元法解決大規(guī)模復(fù)雜構(gòu)件的數(shù)值計(jì)算問題提供了可行的思路。

胡舸[37]基于邊界元法利用MATLAB軟件開發(fā)了腐蝕電場陰極保護(hù)的仿真求解軟件，該軟件可以準(zhǔn)確地計(jì)算均勻電解質(zhì)中的電位分布，且對(duì)安裝有犧牲陽極的海底管線腐蝕電場電位分布的測量有普適性。梁成浩等[38]采用極化曲線作為陰陽極的邊界條件，并對(duì)其進(jìn)行三階段線性擬合處理。利用MATLAB軟件編制程序計(jì)算了凍土層碳鋼熱管陰極保護(hù)系統(tǒng)中的電位分布，為穩(wěn)定凍土層路基和管道的陰極保護(hù)設(shè)計(jì)提供了依據(jù)。

侯靜等[39]通過邊界元法對(duì)海底管道的陰極保護(hù)設(shè)計(jì)進(jìn)行了評(píng)估，并計(jì)算了管道涂層破損率對(duì)管道電位分布的影響。陳靜[40]采用線性擬合的陰極極化曲線作為邊界條件，利用邊界元法、鏡像法以及軸對(duì)稱法求得了管道的電位分布并優(yōu)化了深井陽極的位置，從而得到了達(dá)到保護(hù)要求且電位分布均勻的區(qū)域性陰極保護(hù)優(yōu)化方案。

近幾年，邊界元法因具有可以將計(jì)算域進(jìn)行降維處理、所需數(shù)據(jù)量小、計(jì)算時(shí)間短、計(jì)算精度高等優(yōu)點(diǎn)得到了廣泛的關(guān)注。如由英國Computational Mechanics BEASY集團(tuán)研發(fā)的BEASY CP軟件就是其中的典型代表，從應(yīng)用角度來說，BEASY CP軟件具有通用性強(qiáng)、邊界條件設(shè)定簡單等優(yōu)點(diǎn)。

3 三種數(shù)值解法的比較

以上三種方法在研究陰極保護(hù)體系，尤其是在海底管道陰極保護(hù)體系上都取得過成功，但有限差分法和有限元法的共同點(diǎn)是計(jì)算時(shí)必須對(duì)全部計(jì)算域進(jìn)行網(wǎng)格劃分，致使所需數(shù)據(jù)量大、計(jì)算量大、計(jì)算精度低，而邊界元法只需對(duì)被保護(hù)管道的邊界進(jìn)行網(wǎng)格劃分，克服了其他兩種方法的缺點(diǎn)，成為陰極保護(hù)設(shè)計(jì)最具前景的數(shù)值計(jì)算方法，但其無法對(duì)非均勻介質(zhì)系統(tǒng)進(jìn)行計(jì)算。邊界元法與有限差分法以及有限元法的比較如表1所示。

4 管道陰極保護(hù)數(shù)值模擬方法存在的不足及發(fā)展趨勢

利用數(shù)值計(jì)算對(duì)陰極保護(hù)體系的電位分布進(jìn)行研究已經(jīng)逐漸得到了國內(nèi)外專家學(xué)者的重視，從某種意義上說，陰極保護(hù)已經(jīng)進(jìn)入了數(shù)值仿真時(shí)代。但是應(yīng)用數(shù)值模擬計(jì)算陰極保護(hù)電位分布時(shí)仍存在以下問題：

(1) 數(shù)學(xué)模型的建立問題

目前用于數(shù)值計(jì)算的模型大都是分布型模型，但實(shí)際陰極保護(hù)體系的狀態(tài)是隨時(shí)間發(fā)生變化的，腐蝕介質(zhì)也并非均勻，因此在數(shù)學(xué)模型的建立方面有待進(jìn)一步研究，應(yīng)在計(jì)算模型中考慮介質(zhì)不均勻性以及環(huán)境參數(shù)如土壤含水率、溫度等因素周期性變化對(duì)管道陰極保護(hù)電位分布的影響。

表1 有限差分法、有限元法以及邊界元法的對(duì)比Tab. 1 The comparison of finite difference method,finite element method and boundary element method

(2) 陰陽極邊界條件的處理

在陰極保護(hù)電位分布的數(shù)值計(jì)算中陰陽極極化曲線的確定是至關(guān)重要的邊界條件，但是由于缺乏對(duì)輔助陽極導(dǎo)電機(jī)理、陰極極化以及垢層隨時(shí)間變化規(guī)律的深入研究，迄今尚未建立完善的陽極和陰極表面極化模型，因此陰陽極的極化特性也是今后研究的重點(diǎn)。同時(shí)，通過測量金屬材料在各種介質(zhì)中的極化數(shù)據(jù)，建立極化曲線數(shù)據(jù)庫，也是陰極保護(hù)仿真技術(shù)走向工程應(yīng)用的一個(gè)發(fā)展方向。

(3) 求解方法單一

目前常用的數(shù)值方法中，邊界元法的計(jì)算域?yàn)闊o限域、無法處理非均質(zhì)問題，有限元法正好相反。而實(shí)際問題的求解域往往是半無限域且求解域中介質(zhì)也并非均勻分布，若采用單一的方法來解決問題，只能將某一邊界簡化處理、假設(shè)介質(zhì)處于均勻或分區(qū)域均勻狀態(tài)，這就使其與實(shí)際的陰極保護(hù)體系差別較大，從而產(chǎn)生誤差。一種新的研究動(dòng)向是將幾種方法組合使用，從而提高解決實(shí)際問題的能力。

(4) 數(shù)值計(jì)算技術(shù)存在不足

包含多條管道且管道連接方式復(fù)雜的區(qū)域性陰極保護(hù)技術(shù)已逐漸受到關(guān)注，而如今的數(shù)值計(jì)算技術(shù)無法精確處理復(fù)雜的問題，如邊界元法是目前最具前景的計(jì)算方法，但其形成的系數(shù)矩陣是非對(duì)稱滿陣，求解效率隨求解規(guī)模的增大快速下降，其存儲(chǔ)量和計(jì)算量嚴(yán)重制約著邊界元法的發(fā)展。近幾年，劉立祺等[35-36]提出了一種隨著快速多級(jí)算法發(fā)展起來的快速多級(jí)邊界元解法，其可以加快求解速率，為解決復(fù)雜的陰極保護(hù)體系提供了可行的思路和新的發(fā)展方向。

(5) 國內(nèi)缺乏實(shí)用的陰極保護(hù)設(shè)計(jì)商業(yè)軟件

目前一些專家學(xué)者如孫吉星[22]、陳靜[40]等都對(duì)陰極保護(hù)系統(tǒng)的電位分布展開了研究，并編寫了一些針對(duì)性的程序，同時(shí)部分商業(yè)軟件如ANSYS、MATLAB等都可以用于陰極保護(hù)電位分布的計(jì)算，但是其針對(duì)性不強(qiáng)，這就帶來了一定的系統(tǒng)誤差，國內(nèi)目前還沒有針對(duì)陰極保護(hù)設(shè)計(jì)的綜合專業(yè)軟件，所以完善商業(yè)軟件的研究也是今后的重點(diǎn)。

5 結(jié)束語

目前，在油氣管道事業(yè)飛速發(fā)展的同時(shí)，管道面臨的腐蝕風(fēng)險(xiǎn)也不斷增加，作為控制腐蝕的主要措施，陰極保護(hù)在防止管道腐蝕方面的作用也顯得越發(fā)重要。數(shù)值模擬計(jì)算的出現(xiàn)為陰極保護(hù)系統(tǒng)的設(shè)計(jì)提供了新思路，但同時(shí)也存在著一定的缺點(diǎn)與不足。本文介紹了三種常用的數(shù)值計(jì)算方法以及國內(nèi)外研究學(xué)者對(duì)這三種數(shù)值計(jì)算方法的應(yīng)用情況，并對(duì)比分析了其各自的優(yōu)缺點(diǎn)，同時(shí)還指出了目前管道陰極保護(hù)數(shù)值模擬技術(shù)存在的不足,并提出了相應(yīng)的解決措施，為以后的研究提供了參考資料，對(duì)陰極保護(hù)數(shù)值計(jì)算的長遠(yuǎn)發(fā)展具有積極的影響。

[1] 劉凱,馬麗敏,陳志東,等. 埋地管道的腐蝕與防護(hù)綜述[J]. 管道技術(shù)與設(shè)備,2007(4):36-38.

[2] 李佳奇. 埋地金屬管道的管外陰極保護(hù)數(shù)值模擬研究[D]. 西安:西安石油大學(xué),2015.

[3] 馬偉平,張國忠,白宗成,等. 區(qū)域性陰極保護(hù)技術(shù)研究進(jìn)展[J]. 油氣儲(chǔ)運(yùn),2005,24(9):38-42,48.

[4] 馬含悅,杜磊,楊陽,等. 區(qū)域陰極保護(hù)技術(shù)概況及其發(fā)展[J]. 腐蝕與防護(hù),2014,35(5):425-429.

[5] 翁永基. 陰極保護(hù)設(shè)計(jì)中的模擬研究及其應(yīng)用[J]. 腐蝕科學(xué)與防護(hù)技術(shù),1999,11(2):99-110.

[6] 郝宏娜,李自力,王太源,等. 陰極保護(hù)數(shù)值模擬計(jì)算邊界條件的確定[J]. 油氣儲(chǔ)運(yùn),2011,30(7):504-507,472.

[7] 郝宏娜,李自力,王太源,等. 非均勻土壤中陰極保護(hù)電位數(shù)值模擬計(jì)算[J]. 化工機(jī)械,2012,39(4):494-496.

[8] STROMMEN R,RODLAND A. Computerized techniques applied in design of offshore cathodic protection systems[J]. Materials Performance,1981,20(4):15-20.

[9] 張鳴鏑,杜元龍,殷正安,等. 有限差分法計(jì)算海底管道陰極保護(hù)時(shí)的電位分布[J]. 中國腐蝕與防護(hù)學(xué)報(bào),1994,14(1):77-81.

[10] 錢海軍,劉小光,張樹霞,等. 管內(nèi)陰極保護(hù)的數(shù)值模擬(Ⅱ)——有限差分法計(jì)算大口徑管內(nèi)的電位分布[J]. 化工機(jī)械,1997,24(5):33-35,62.

[11] KENNELLEY K J,BONE L,ORAZEM M E. Current and potential distribution on acoated pipeline with holidays Part 1-Model and experimental verification[J]. Corrosion,1993,49(3):199-210.

[12] PARSA M H,ALLAHKARAM S R,GHOBADI A H. Simulation of cathodic protection potential distributions on oil well casings[J]. Journal of Petroleum Science and Engineering,2010,72(3):215-219.

[13] KASPER R G,APRIL M G. Electrogalvanic finite element analysis of partially protected marine structures[J]. Corrosion,1983,39(5):181-188.

[14] MARCASSOLI P,BONETTI A,LAZZARI L,et al. Modeling of potential distribution of subsea pipeline under cathodic protection by finite element method[J]. Materials & Corrosion,2015,66(7):619-626.

[15] CHIN D T,SABDE G M. Modeling transport process and current distribution in a cathodically protected crevice[J]. Corrosion,2000,56(8):783-793.

[16] 李成杰. 深海陰極保護(hù)過程研究及其在溫躍層環(huán)境下數(shù)學(xué)模型的建立[D]. 青島:中國海洋大學(xué),2014.

[17] 宋成立,淡勇. 基于有限元法的埋地管道陰極保護(hù)研究[J]. 全面腐蝕控制,2013,27(10):42-45.

[18] 潘柳依,雷寶剛,范錚,等. ANSYS有限元法在管道陰極保護(hù)中的應(yīng)用[J]. 材料保護(hù),2014,47(3):45-47,58,70.

[19] 徐勛. 基于有限元法的陰極保護(hù)系統(tǒng)數(shù)值模擬研究[D]. 大連:大連理工大學(xué),2014.

[20] 李丹丹,畢武喜,祁惠爽,等. 交叉并行管道陰極保護(hù)干擾數(shù)值模擬[J]. 油氣儲(chǔ)運(yùn),2014,33(3):287-291.

[21] 紀(jì)俊剛. 在海底管道的陰極保護(hù)中數(shù)學(xué)模型的建立與應(yīng)用研究[D]. 青島:中國海洋大學(xué),2008.

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

[23] TELLES J C F,WROBEL L C,MANSUR W J,et al. Boundary elements for cathodic protection problems[M]. Berlin:Springer-Verlag,1985:73-83.

[24] CHUANG J M. Numerical solution of nonlinear boundary value problems arising in corrosion and electroplating modeling with applications to 3D ships and marine structures[D]. Nova Scotia:Technical University of Nova Scotia,1986.

[25] ZAMANI N G,CHUANG J M,PORTER J F. BEM simulation of cathodic protection system employed in infinite electrolytes[J]. International Journal for Numerical Methods in Engineering,1987,24(3):605-620.

[26] SANTIAGO J A F,TELLES J C F. On boundary elements for simulation of cathodic protection systems with dynamic polarization curves[J]. International Journal of Numerical Methods in Engineering,1997,40(14):2611-2627.

[27] DIAZ E S,ADEY R. Optimising the location of anodes in cathodic protection systems to smooth potential distribution[J]. Advances in Engineering Software,2005,36(9):591-598.

[28] ORAZEM M E,ESTEBAN J M,DEGERSTEDT R M,et al. Mathematical models for cathodic protection of an underground pipeline with coating holidays:Part 1-Theoretical development[J]. Corrosion,1997,53(4):264-272.

[29] BRAEIL S L D C,TELLS J C F,MIRANDA L R M. Simulation of coating failure on a cathodically protected pipeline-Experimental and numerical results[J]. Corrosion,2000,56(56):1180-1188.

[30] BRICHAU F,DECONINCK J. A numerical-model for cathodic protection of buried pipes[J]. Corrosion,1994,50(1):39-49.

[31] MUHAREMOVIC A,ZILDZO H,BEHLILOVIC N,et al. Numerical model for calculation of parameters of cathodic protection system with galvanic anodes[J]. International Symposium on Information,Communication and Automation Technologies,2009,33(8):985-1000.

[32] LACERDA L A D,SILVA J M D,LZARIS J. Dual boundary element formulation for half-space cathodic protection analysis[J]. Engineering Analysis with Boundary Elements,2007,31(31):559-567.

[33] 伍欣. 基于檢查片檢測的川氣東送管道陰極保護(hù)有效性評(píng)價(jià)[D]. 四川:西南石油大學(xué),2015.

[34] 趙雷亮. 陽極布置方式對(duì)陰極保護(hù)效果的影響[D]. 青島:中國石油大學(xué)(華東),2013.

[35] 劉立祺,王海濤. 基于三維管道模型的快速邊界元法在陰極保護(hù)分析中的應(yīng)用[J]. 清華大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,55(9):1003-1009.

[36] 張東東,杜翠薇,程學(xué)群,等. 快速多極邊界元法在陰極保護(hù)求解電位分布中的應(yīng)用[C]//中國科學(xué)技術(shù)協(xié)會(huì)年會(huì).[出版地不詳]：[出版者不詳],2010.

[37] 胡舸. 海底管線腐蝕檢測與腐蝕預(yù)測的研究[D]. 重慶:重慶大學(xué),2007.

[38] 梁成浩,袁傳軍,黃乃寶. 邊界單元法計(jì)算凍土層管道陰極保護(hù)時(shí)的電位分布[J]. 大連海事大學(xué)學(xué)報(bào),2011,37(4):109-112,116.

[39] 侯靜,常世君,鄭國良. 數(shù)值模擬計(jì)算在海底管道陰極保護(hù)設(shè)計(jì)中的應(yīng)用[J]. 新技術(shù)新工藝,2015(5):15-18.

[40] 陳靜. 輸氣站場區(qū)域性陰極保護(hù)數(shù)值模擬及陽極位置優(yōu)化研究[D]. 四川:西南石油大學(xué),2009.

AppliedDevelopmentofNumericalCalculationMethodsinCathodicProtectionofPipelines

As to the situation about numerical simulation for cathodic protection of pipelines, the equations of cathodic protection system and their three common boundary conditions are presented. Moreover, three numerical calculation methods are described and their latest applications are summarized. Also, the advantages and disadvantages about the numerical calculation methods are analyzed. Finally, the problems existing in this field and their solutions are proposed.

cathodic protection; numerical calculation; finite difference method; finite element method; boundary element method

10.11973/fsyfh-201711001

2016-05-03

國家自然科學(xué)基金(51301201)； 山東省自然科學(xué)基金(ZR2013EMQ014)

寇 杰(1969-),教授,博士,主要從事多相管流及油氣田集輸技術(shù)、油氣儲(chǔ)運(yùn)系統(tǒng)安全工程、油氣長距離管輸技術(shù)的研究工作,13969877559,chuyunk@126.com

TG174.41

A

1005-748X(2017)11-0823-06