儲罐底板陰極保護技術現狀及其發展

王 東

(中國石化集團勝利油田，山東 德州 251508)

?

儲罐底板陰極保護技術現狀及其發展

王 東

(中國石化集團勝利油田，山東 德州 251508)

對于儲罐底板內外側的腐蝕，陰極保護是最為經濟有效的防腐措施。根據測量的電位來判斷儲罐底板是否得到有效保護，把保護電位作為調整陰極保護系統運行的重要參數。遇到罐底板中心部位電位不容易檢測的時候，為了保證儲罐底板中心得到充分保護，需要利用數學方法，根據陰極保護電位數值計算理論，建立陰極保護的數學模型，求解得到整個儲罐底板的電位分布。此舉措對于提高整個儲罐底板的陰極保護效果，延長站內儲罐的使用壽命有非常重要的意義。

儲罐底板 陰極保護 電位分布 邊界元法

儲罐是聯合站、集油站內接收、儲存油品的重要設備，儲罐的正常運行對整個站場的安全、穩定、長期生產運行起著十分重要的作用。隨著站內儲罐運行時間延長，儲罐底板安全度下降，一旦發生腐蝕穿孔，會導致油品泄漏直接造成經濟損失，也會對周圍的環境造成嚴重的破壞，甚至引發火災與爆炸等重大安全事故。

陰極保護投資成本僅為儲罐建設成本的1%，可以成倍地延長更換儲罐底板的周期。根據石油天然氣行業標準SY/T 0088-2006 《鋼質儲罐罐底外壁陰極保護技術標準》與NACE RP0193-2001 《在級碳鋼儲罐底的外部陰極保護》的規定，在施加陰極保護的情況下，使用銅/飽和硫酸銅參比電極(CSE)測得罐/地極化電位應為-0.85～-1.2 V。為精確測量保護電位，要求盡可能地降低土壤IR降的影響，參比電極應安置在盡量靠近被測試構件的位置。

1 儲罐底板陰極保護技術現狀

從陰極保護技術誕生到現在，已經發展了190多年。目前許多發達國家已經將在大型油庫、輸油和輸氣站場等強制實施區域性陰極保護寫入法律，并且做到了陰極保護系統和主體工程同時設計、同時施工和同時投產。

國外從20世紀60年代就開始了儲罐底板外側陰極保護技術的應用與研究，國內開始應用儲罐底板的陰極保護只有20多年的歷史。

20世紀70年代，區域性陰極保護技術首次在中國石化和勝利油田東辛采油廠得到利用，站內的儲罐和埋地管線都得到了有效保護[1]。1989年，福建省東山油罐區的5座儲油罐把犧牲陽極材料鋅基埋設在儲罐周邊1m外地下土層內，采用犧牲陽極方法對儲罐底板進行陰極保護。由于犧牲陽極材料的驅動電位較低，所以只適用于小型的儲罐。目前對于儲罐底板的陰極保護工程，國內外多選擇強制電流陰極保護技術。

外加強制電流陰極保護技術是防止儲罐底板外側腐蝕的有效手段，網狀陽極是目前最為理想的儲罐底板外加強制電流陰極保護系統的陽極形式，給儲罐底板提供均勻的保護電流，使儲罐底板中心和周圍的電位分布均勻。對于已建儲罐底板施加陰極保護，常采用深井陽極的形式，此形式下儲罐底板保護電位沿直徑方向分布不均勻。

目前可以通過儲罐底板下部埋設了參比電極的方法，測試儲罐底板外側的保護效果，而長效參比電極埋設在儲罐底板周圍時，則只能測得儲罐底板邊緣的保護電位。因此，為了實現對深井陽極形式下儲罐底板外側的保護效果的準確預測，有必要對儲罐底板外側的電位分布進行研究。

2 儲罐底板陰極保護模擬技術

2.1 電位解析方法研究

解析計算公式比較直觀，在工程應用中比較方便。長距離埋地管道陰極保護電位和電流密度分布的解析式已經在工藝計算和設計當中得到廣泛應用，儲罐底板陰極保護電位分布不能直接套用管道的陰極保護電位分布公式。

在整個儲罐底板上的電流密度非常均勻的假設前提下，Morgan推導得到了儲罐底板陰極保護的電位分布公式，電位與距罐中心距離呈線性關系，但是與實際電位分布存在很大差異，只適用于罐底板直徑非常小的情況。后來Morgan考慮了儲罐底板電流密度的不均勻性問題，從帶電圓盤電荷分布公式出發，推導得到距罐底中心距離為處罐底的保護電流密度為：

(1)

式中:x為距罐底中心的距離，m；r為罐底板半徑，m；is為罐底平均電流密度，A/m2。

李湘怡[2]等將被保護儲罐底板和其相鄰土壤介質當作兩個體系，一個為陽極電場，另一個為陰極電場，陽極電場是由陽極流向大地的無窮遠處產生的，而陰極電場是由大地的無窮遠處流入儲罐底板時產生的。利用穩流電場理論推導出罐底陰極保護電位分布模型解析公式，并驗證了該數學模型，利用簡化的解析公式可以方便計算儲罐底板電位分布，為儲罐底板陰極保護的評價和檢測提供理論支撐。

φ=Um-(Ua+Uc)

(2)

式中:Um為罐底板上流動的電流在罐底板產生的電位，V；Ua為陽極流向大地無窮遠處的電流在罐底板產生的電位，V；Uc為大地無窮遠處流入罐底板的電流在罐底板產生的電位，V。

翁永基[3]等將式(1)級數展開，省略高次項，引入描述電流密度不均勻性的經驗常數k，得到距離罐底中心x處電流密度為：

(3)

式中:y為距罐底中心的比例距離，m；i0為罐底中心的電流密度，A/m2；k為經驗常數，取30。

李湘怡以式(3)為基礎，根據三種恒流電場理論，分別計算三種穩流電場的電位分布，疊加各個電場電位分布，罐底板上流動的電流在罐底板產生的電位太小，可以忽略，當不考慮陽極電位影響時，可以得到陰極保護條件下儲罐電位分布解析公式：

(4)

式中：ρ土壤電阻率，Ω·m。

以三種恒流電場的疊加理論為基礎，梁宏進一步討論了單陽極、雙陽極、三陽極和四陽極等幾種陽極形式的陽極電位分布，得出單陽極形式下罐底板電位分布的不均勻程度比其他三種形式下大。耿曉梅也利用了三種穩流電場理論建立了站庫追加陰極保護系統的設計計算模型，并利用序列二次規劃方法優化了陽極位置，使被保護區域在最小的陽極電流情況下得到保護。

杜艷霞[4]改進了式(1)，引入電流密度修正系數，采用灰色關聯法確定了陽極埋深、數量、罐底板直徑是影響的主要因素，得到距離罐底中心x處電流密度為：

(5)

式中：δ為電流密度修正系數，無量綱。

根據電場疊加理論，應用電磁學基本原理和儲罐底板陰極保護系統的幾何特征，推導了罐旁單側布置深井陽極和罐周圍對稱布置兩口深井陽極情況下儲罐底板外側陰極保護電位分布的解析公式。根據數值模擬結果，數據擬合得到了與陽極埋深、數量、罐底板直徑的關系式。

儲罐底板陰極保護電位分布解析公式基本都是采用三種恒流電場理論，未考慮儲罐底板的非線性極化特性。同時解析公式受到陽極的布置方式和陽極數量的限制，對罐旁單側深井陽極陰極保護電位分布的研究較多，其它復雜陽極形式下的儲罐底板陰極保護電位分布研究較少。

2.2 電位數值模擬技術研究

隨著電化學技術和計算機技術的快速發展，為人們獲取儲罐底板陰極保護電位的分布提供了有效的途徑，人們嘗試采用數值計算的方法來獲取被保護體表面的電位分布。20世紀80年代，國外已有學者開始利用計算機技術進行陰極保護設計工作。如今人們可以利用現成的計算程序軟件包如FLUENT，ANSYS和MATLAB進行數值計算，節省編制復雜程序的時間。

2.2.1 數學模型

(1)控制方程

基于恒電流理論，對深井陽極儲罐底板電位分布公式進行了推導，推導得到Poisson方程為電位的控制方程，根據實際狀況，以研究區域絕緣面電流為零，無窮遠處電位為零，實驗室實測的極化曲線為儲罐底板外側陰極保護電位分布問題的邊界條件，建立了儲罐底板陰極保護電位分布數學模型。

(6)

當所研究的恒定電場中存在恒定的點源，電位分布滿足Poisson方程，在給定邊界條件時，可通過求解該方程可以得到電位的分布。

(2)邊界條件

根據儲罐底板陰極保護電位分布的物理模型，所研究的區域為Ω，被表面Γ所包圍，總邊界Γ包括：無限大地面表面邊界Γd、儲罐底板外表面邊界Γc、虛設半無限域土壤球冠邊界Γ∞。

地表面邊界Γd與大氣相接觸，外加陰極保護電流只能在土壤介質中流動，不能流入大氣，因此可以把地表面當作絕緣面，其邊界條件滿足：

(7)

可以忽略陽極對無限大土壤表面的影響，所以在無窮遠處表面的電流密度也為零，即滿足以下邊界條件：

(8)

陰極的邊界條件采用描述電流密度和電位之間函數關系的極化函數來定義：

J=f(φ-φeq)

(10)

式中：φ為儲罐底板施加陰極保護后的極化電位值；φeq為儲罐底板在土壤環境中的自腐蝕電位，f(φ-φeq)為陰極極化函數。

隨著電化學技術的發展，人們開始認識到極化電位和電流密度有關，人們可以通過電化學實驗測得描述這種關系的陰極極化曲線。Ashworth提出建立正確的極化模型將直接影響到整個陰極保護模型的準確性。葛艾天根據Butler-Volmer極化動力學原理，考慮陽極極化過程由活化極化控制，陰極極化過程由活化極化和濃差極化混合控制，推導出了對于管道和鋼質管件采用的陰極邊界條件。蔣卡克采用BEASY邊界元數值模擬軟件計算了實際工程中，外加電流陰極保護系統中，輔助陽極采用恒電流密度與陽極極化曲線作為邊界條件模擬誤差對比分析。

一般陰極電位和電流密度的關系是非線性的，其模擬處理直接影響計算結果的精度和計算的時間，常采用的是三種方法：一是公式法，Danson認為電位和電流密度之間存在類似塔菲爾方程的關系；二是Newton-Raphson迭代法，該方法是將非線性擬合得到的函數用泰勒級數展開，只取展開項的前兩項作為函數的近似值；三是分選線性擬合方法，將實測極化曲線分為若干段，每一段都進行線性擬合，分段線性擬合對擬合實驗測得的極化曲線有非常好的靈活性。

2.2.2 數值模擬技術應用現狀

可用包含邊界條件的數學物理方程來描述陰極保護場域的數學模型，通過數值計算方法求解偏微分方程，用網格離散連續場域，包含邊界的微分方程被含有限未知數的代數方程組所代替，通過求解代數方程組得到未知節點上的函數值大小。目前應用于腐蝕場的數值計算方法主要包括有限差分法、有限元法和邊界元法。

邱楓[5]等運用有限元法，利用計算機編寫了程序，計算了儲罐底板外側陰極保護時的電位和電流分布，討論了土壤電阻率、儲罐底板涂層、陽極距儲罐近陽極點的距離、陽極極化值、陽極的數量和布置方式等多種因素對電位分布的影響。

馬偉平[6]等從陰極保護數學模型控制方程和熱流體能量方程相似性入手，把恒定電場問題推廣到研究熱流場。討論如何設定儲罐底板、計算區域邊緣和輔助陽極的邊界條件，利用兩種電流密度分布假設，對深井陽極儲罐底板的陰極保護進行了數值模擬。

杜艷霞[7]考慮陰極條件由活化極化和濃差極化混合控制，改進了儲罐底板外側陰極保護電位分布數學模型，利用國際上流行的軟件包FLUENT模擬計算了儲罐底板外側的陰極保護電位分布，與室內建立的實驗模擬裝置測量結果對比分析，模擬數據和實驗數據比較吻合，驗證了求解方法的準確性。進一步考察了實際生產中常用的六種陽極埋設方式儲罐底板陰極保護電位分布規律以及保護效果進行了對比分析。并考察了土壤電阻率、儲罐底板直徑、接觸介質性質、陽極埋深和數量、與儲罐水平距離等因素對電位分布的影響。

齊建濤[8]等采用了邊界元方法對儲罐罐底水相區的陰極保護進行了數值模擬和仿真，通過調整外層陽極的數量和陽極布置的位置來對陰極保護系統進行了優化，提高了側壁的保護效果，改善了儲罐底板上的過保護問題，使水相區電位和電流密度分布也更加均勻。

張豐[9]等采用陰極保護領域應用比較廣泛的BEASYCP軟件對深井陽極和近網狀陽極對儲罐底板電位分布進行分析，通過電化學實驗測得A3鋼在土壤溶液、土壤、干燥瀝青砂和含水瀝青砂中的極化曲線。

夏宗春[10]采用邊界元法計算了罐底板外側陰極保護電位分布，討論了多種影響電位分布的因素，包括土壤電阻率大小，儲罐罐底板外側的表面狀態，陽極的數量和分布情況等。

國內外對埋地管道電位分布數值模擬的文獻較多，由于受陽極布置方式、儲罐底板非線性極化特性等因素影響，對深井陽極罐底板電位分布進行數值模擬的文獻相對較少。隨著電化學技術和計算機技術的發展，采用數值計算方法具有顯著的優勢，能夠更準確、高效地計算儲罐底板的電位分布，研究陽極設置參數、罐底板直徑和土壤電阻率對電位分布的影響，也成為國內外研究的熱點。

3 結論和建議

(1)強制電流陰極保護技術是防止儲罐底板外側腐蝕的有效手段，目前大型儲罐的陰極保護技術雖然得到了大規模的推廣應用，但是不夠完善。

(2)由于受介質分布不均勻、陽極布置方式、儲罐底板非線性極化特性等因素影響，儲罐底板陰極保護電位分布解析方法方面的研究十分困難。

(3)采用數值計算的方法來獲取被保護體表面的電位分布。相對于傳統設計方法，計算機技術可以高效、準確地判斷儲罐底板陰極保護的效果。

[1] 馬偉平, 張國忠. 區域性陰極保護技術研究進展[J]. 油氣儲運, 2005, 24( 9): 38-42.

[2] 李相怡, 翁永基. 陰極保護電位分布數學模型的研究及其應用[J]. 腐蝕科學與防護技術, 1998, 10(2): 64-69.

[3] 翁永基. 陰極保護設計中的模型研究及其應用[J]. 腐蝕科學與防護技術, 1999, 11(2): 99-110.

[4] 杜艷霞. 金屬儲罐底板外側陰極保護電位分布規律研究[D]. 青島：中國石油大學(華東), 2007.

[5] 邱楓, 徐乃欣. 鋼質貯罐底板外側陰極保護時的電位分布[J]. 中國腐蝕與防護學報, 1996, 16(1): 29-36.

[6] 馬偉平, 張國忠, 梁昌華，等. 兩點法電位檢測儲罐底板陰極保護效果評價[J]. 油氣儲運, 2005, 24(4): 40-42.

[7] 杜艷霞, 張國忠. 儲罐底板外側陰極保護電位分布的數值模擬[J]. 中國腐蝕與防護學報, 2006, 26(6): 346-349.

[8] 齊建濤, 李焰. 海運石腦油儲罐水相區陰極保護電場的數值模擬[J]. 海洋科學, 2012, 36(11): 79-80.

[9] 張豐, 陳洪源, 李國棟，等. 數值模擬在管道和站場陰極保護中的應用[J]. 油氣儲運, 2011, 30(3): 208-211.

[10] 夏宗春. 地面鋼質貯罐底板外側陰極保護電位分布計算[D]. 大連：大連理工大學, 2000.

(編輯 張向陽)

Overview and Development of Tank Bottom Cathodic Protection Technology

WangDong

(LinpanOilProductionPlantofShengliOilField,Dezhou251508,China)

At present, cathodic protection is the most effective and economic anti-corrosion measure for the internal and external sides of bottom plates. Whether the tank bottom plate has been effectively protected or not can be determined according to the measured electric potential parameters. Protection potential is an important parameter to adjust the operation of cathodic protection system. When it is difficult to get the protective electric potential of the tank center bottom and it is important to protect the tank center bottom plate, it is necessary to establish mathematical model of cathodic protection to get the electric potential distribution of the whole bottom plate by using numerical simulation based upon numerical calculation theory for cathodic protection electric potential. It is very significant for improving effectiveness of cathodic protection and extending the service life of tanks in the station.

tank bottom plate, cathodic protection, potential distribution, boundary element method

2016-01-06；修改稿收到日期：2016-02-25。

王東(1988-)，碩士，勝利油田臨盤采油廠壓氣一站主要從事腐蝕與安全方面研究。E-mail:shengli_2006@163.com