王彥然，肖 杰，范永昭，張寅暉，田 源，李 珊

(1. 中國石油西南油氣田公司天然氣研究院，四川 成都 610213；2.國家能源高含硫氣藏開采研發(fā)中心，四川 成都 610213；3. 中石油集團公司高含硫氣藏開采先導試驗基地，四川 成都 610213；4. 中國石油西南油氣田公司氣田開發(fā)管理部，四川 成都 610051)

0 前 言

酸性氣田的管線在CO2和H2S的腐蝕作用下，常常具有較高的腐蝕速率，導致管道的壁厚快速減薄，甚至在短時間內(nèi)腐蝕穿孔，給油氣生產(chǎn)帶來了嚴重的損失。因此，準確預測CO2和H2S腐蝕下管線的腐蝕速率對于開展腐蝕防控工作具有重要意義。

近幾十年來，國內(nèi)外不少學者致力于CO2和H2S腐蝕速率預測模型的研究。隨著研究思路不斷拓寬和對已有成果的不斷完善，預測模型考慮的影響因素越來越全面。由于不同模型的建模方法和考慮的因素存在差異，所以它們各自的適用條件也不相同。本文綜述了CO2腐蝕和CO2/H2S共存體系的腐蝕預測模型研究進展，對選擇模型解決實際問題和預測模型的研究具有指導意義。

1 CO2腐蝕速率預測模型

CO2腐蝕速率預測模型可分為3大類型：一是基于實驗室數(shù)據(jù)或現(xiàn)場數(shù)據(jù)建立的經(jīng)驗模型；二是結合熱力學、動力學原理和實驗數(shù)據(jù)建立的半經(jīng)驗模型；三是通過分析微觀機理建立的機理模型。

1.1 經(jīng)驗模型

根據(jù)低溫實驗室數(shù)據(jù)和高溫現(xiàn)場數(shù)據(jù)總結得到的Norsok M506模型是最基礎的經(jīng)驗模型，該模型考慮了溫度、總壓、CO2含量、pH值、管壁剪切應力、流速等因素，提出了pH值和管壁切應力的計算模塊[1]。該模型合理地評估了產(chǎn)物膜的保護性，高溫條件下的預測精度比低溫條件更高，適用條件為溫度100～150 ℃，pH值3.5～6.5[2]。該模型的局限性包括：第一，只能預測均勻腐蝕；第二，管壁切應力計算模塊是基于單相流理論設計的；第三，未涉及產(chǎn)物膜形成的動力學過程及其厚度、孔隙率、機械強度等的影響[3]。這些局限性使得該模型對于局部腐蝕和多相流的情況不適用。Norsok M506模型預測的腐蝕速率見式(1)：

(1)

式中：vcorr是腐蝕速率，mm/a；Kt是與溫度和腐蝕產(chǎn)物膜相關的常數(shù)；S是管壁切應力，Pa；fCO2是CO2的逸度，bar；f(pH)t是溶液pH值對腐蝕速率的影響因子。

但在油氣輸送中，多相流更普遍，所以有必要借助多相流的理論成果，建立多相流條件下對管道腐蝕的預測模型[4]。Jepson通過研究多相流腐蝕介質(zhì)中的化學過程，提出了一種適用于水平段塞流條件的CO2腐蝕速率預測模型。該模型考慮了與段塞流和原油相關的參數(shù)的影響，還特別對有無產(chǎn)物膜分2種情況進行選擇[4]。同時，該模型對pH值敏感，pH值大于5時，預測值可能比Norsok M506模型低。Jepson模型預測的腐蝕速率見式(2)：

(2)

Corpos模型也是研究多相流腐蝕的模型，它借助Norsok模型，根據(jù)管道沿線的pH值、潤濕狀態(tài)等信息，計算管道沿線不同位置的腐蝕速率，在低含水率條件下，該模型的預測結果低于Norsok模型。還有針對氣液兩相流特點提出的Lafayette 模型，它根據(jù)流態(tài)差異，分分層流、間歇流和環(huán)狀流3種情況分別建立腐蝕速率的計算公式[5]。

表1分析了經(jīng)典的CO2腐蝕預測經(jīng)驗模型，供讀者選擇模型時參考。

表1 經(jīng)典的CO2腐蝕預測經(jīng)驗模型的分析情況

1.2 半經(jīng)驗模型

最早建立CO2腐蝕預測模型的是De Waard和Millams[6]，他們設置恒定pH值的實驗條件，研究CO2分壓和溫度的影響，于1975年建立了De Waard 75的半經(jīng)驗模型，該模型忽略了很多影響因素，造成預測值偏大。此后每多考慮一個因素，就引入一個與之對應的校正因子降低預測值，基于這種思路，De Waard 91模型將pH值、總壓力、烴、腐蝕產(chǎn)物膜等因素引入模型[7]，在此基礎上，De Waard 93模型又考慮了流速的影響[8]。之后再提出的De Waard 95模型是最全面、應用范圍最廣的版本，它關注了腐蝕反應中的質(zhì)量傳遞過程，結合腐蝕反應動力學和流動CO2傳質(zhì)過程2方面，而且還特別考慮了材料的差異，提出了與材料性質(zhì)有關的影響因子；但是，該模型保守估計了產(chǎn)物膜的保護作用，而高溫和高pH值條件下，產(chǎn)物膜可能會較大程度地抑制腐蝕，導致該模型得到的預測值遠高于真實值[9，10]，所以，100～150 ℃的條件下，De Waard 95模型的預測精度低于Norsok M506模型，經(jīng)過不斷改進，De Waard 95模型形成以下公式：

(3)

式中：vcorr是腐蝕速率，mm/a；T是溫度，K；PCO2為CO2分壓，bar；pHactual是溶液的pH值；pHCO2是在實際的溫度和壓力下，飽和水溶液的pH值。

De Waard模型的提出為建立半經(jīng)驗模型提供了方法和依據(jù)，國內(nèi)外學者通過借鑒De Waard 模型的思想建立了其它的CO2腐蝕半經(jīng)驗模型。

Cassandra模型[11]以De Waard 模型為基礎，引入成膜溫度的概念表征產(chǎn)物膜的影響，認為溫度高于成膜溫度時，產(chǎn)物膜會抑制腐蝕，溫度對預測值的影響很小。相對于De Waard 模型，該模型的不足是忽略了油膜的保護作用。李頔[12]總結得到B. Mishra模型借鑒了De Waard的實驗方法，考慮了溫度、CO2分壓和H+濃度的影響，其它因素的影響通過常數(shù)C修正。但是，該模型存在忽略了產(chǎn)物膜的影響的缺陷。

ECE模型[13]是基于De Waard95 模型的理論和Norway 的IFE環(huán)流測試數(shù)據(jù)建立的。該模型特別考慮了原油的緩蝕作用，以乳化臨界點判斷油和水相態(tài)，高于乳化臨界點說明含水率不高，形成油包水相，預測的腐蝕速率值較低，反之，說明含水率較高，得到的腐蝕速率也較大。Predict模型也充分考慮了產(chǎn)物膜和油膜的抑制作用，認為產(chǎn)物膜的保護作用隨pH值升高而提高，所以pH值大于5的條件下，Predict模型的預測精度會受影響，預測值偏低[14]。

田洋陽等[15]對溫度和pH值等參數(shù)的影響進行深入分析后，提出了De Waard 95模型的改進形式。崔鉞等[16]認為De Waard模型未考慮流場的誘導，故綜合De Waard模型、沖蝕模型和流場模擬結果，建立了流場誘導下的CO2腐蝕預測模型，應用時可精確至流場誘導位置預測腐蝕速率。

表2總結了經(jīng)典的CO2腐蝕預測半經(jīng)驗模型的主要考慮因素，供讀者選擇模型時參考，表2中打√表示建立該模型所考慮的相關因素，空白處表示不考慮該因素。

表2 經(jīng)典的CO2腐蝕預測半經(jīng)驗模型的主要考慮因素

1.3 機理模型

Nesic等[17，18]、Gray等[19]等對機理模型的研究較多。最初通過分析CO2腐蝕的陰陽極反應，結合流速、pH值等因素，建立了與De Waard模型預測結果相近的機理模型。此后，再進一步將表面電化學和化學過程、微觀傳質(zhì)過程與模型聯(lián)系起來，又考慮了產(chǎn)物膜的保護性，修正了模型，拓寬了適用范圍。2003年，Nesic等[18]在之前的模型基礎上，又深入分析了產(chǎn)物膜厚度和孔隙度的影響規(guī)律，考察了物質(zhì)沉積和溶解的平衡過程，建立了FeCO3沉積增長的動力學模型，因為pH值和溫度對產(chǎn)物膜特征影響較大，所以該模型對pH 值和溫度較敏感。盡管該模型對產(chǎn)物膜性質(zhì)的認識還不夠到位，但提出的成膜傾向因子對認識腐蝕機理具有指導意義，且該模型在低流速條件下的預測精度高于De Waard 95模型。

Tulsa模型屬于單相流機理模型，未考慮油膜的影響。該模型的理論基礎是傳質(zhì)動力學和電化學反應，同時適用于沒有產(chǎn)物膜生成、有產(chǎn)物膜覆蓋和產(chǎn)物膜破損的情況。模型對流速和pH值較敏感，pH值高于5的條件下產(chǎn)物膜保護性較好，所以模型預測的腐蝕速率很小。Gray等[19]研究了恒定 pH 值的CO2飽和溶液中H2CO3還原的電荷傳遞規(guī)律，并考慮了CO2緩慢水解的影響，建立了相應的機理模型。后來經(jīng)修正，將適用的pH值范圍拓寬為6～10。但是，該模型沒有對產(chǎn)物膜方面的機理進行研究。Dayalan等[20]提出了計算電極表面反應物濃度的方法，為研究FeCO3沉積模型奠定了理論基礎。王獻昉等[21]以Nesic模型的理論為基礎，考慮溫度和流速對產(chǎn)物膜的影響，提出了新的CO2腐蝕產(chǎn)物膜影響因子計算公式。針對酸性氣體溶于管道頂部的冷凝液中引起腐蝕的現(xiàn)象，Zhang等[22]建立了Top of Line模型。郭少強等[23]研究了管道頂部液滴冷凝過程、液滴內(nèi)液相化學變化和CO2腐蝕的規(guī)律，利用鐵元素守恒將3個過程統(tǒng)一起來，建立了預測油氣輸送管線頂部CO2腐蝕的模型。

以上機理模型的建立都從腐蝕的本質(zhì)出發(fā)，預測精度不受氣田腐蝕環(huán)境影響，擁有更廣泛的適用性，但機理模型尚未與氣田生產(chǎn)相結合，在實際應用中，較難從現(xiàn)場獲取到機理模型需要的數(shù)據(jù)，所以目前生產(chǎn)中仍普遍采用經(jīng)驗模型或半經(jīng)驗模型預測腐蝕[24-25]。

綜上所述，CO2腐蝕在經(jīng)驗模型、半經(jīng)驗模型和機理模型3方面都已建立了一些經(jīng)典的腐蝕速率預測模型。早期建立的模型為后續(xù)的研究提供了方法，通過不斷的改進完善，各種新模型考慮的因素越來越全面，也更加貼近酸性氣田的實際腐蝕環(huán)境。

2 H2S/CO2共存體系的腐蝕速率預測模型

在酸性氣田的實際腐蝕環(huán)境中，存在H2S的環(huán)境往往也伴隨著CO2，所以有關單一H2S腐蝕的預測模型很少，往往直接建立H2S/CO2共存體系的腐蝕預測模型。

目前，對于H2S/CO2共存體系，預測腐蝕速率的思路通常是：先根據(jù)分壓比判斷哪類氣體在腐蝕反應中占主導，將它作為主要因素，另一種氣體則作為影響因素校正計算結果。

張清等[26]分別推導了CO2主導腐蝕過程和H2S主導腐蝕過程時的預測模型，再疊加2個分模型得到總模型[式(4)]，該模型將2類氣體的分壓參數(shù)直接引入公式，打破了利用影響因子修正預測值的模式，為研究H2S / CO2共存體系的腐蝕速率預測模型奠定了基礎。

lnvcorr=k(lnPCO2+lnPH2S)+a(lnPH2S)2+blnPH2S+c

(4)

式中：vcorr是腐蝕速率，mm/a；PH2S是H2S分壓，MPa；PCO2為CO2分壓，MPa；k、a、b、c為待定常數(shù)，與氣體分壓、溫度、流速等因素有關。

艾志久等[27]借助這類方法，更新了H2S為主導，CO2作為影響因素時的預測模型，見式(5)：

lnvcorr=KlnPCO2+(K+A)lnPH2S+B(lnPH2S)2+C

(5)

式中：vcorr是腐蝕速率，mm/a；PH2S是H2S分壓，MPa；PCO2為CO2分壓，MPa；K、A、B、C是待定常數(shù)，與氣體分壓、溫度、流速等因素有關。

但是這類模型基本只對H2S分壓和CO2分壓2個因素進行了較深入的考慮，而沒有直接將H2S/CO2共存體系下的腐蝕速率與溫度、流速、腐蝕時間等變量關聯(lián)起來。因此，近年來不少學者基于以前的研究成果，針對某些高含硫氣田的腐蝕環(huán)境，嘗試建立同時考慮氣體分壓、溫度等多種因素的預測模型或提出了某些因素的校正方法?？紤]更多因素后，在特定的腐蝕環(huán)境下，修正模型具有更高的預測精度。

基于腐蝕量與時間在雙對數(shù)坐標中的線性關系，張智等[28]設計實驗確定相關參數(shù)，得到了揭示腐蝕速率與時間之間關系的預測模型，但該模型只能計算均勻腐蝕速率。也有學者嘗試將其它變量與腐蝕速率的關系引入經(jīng)典模型。廖柯熹等[29]以H2S主導腐蝕反應的模型為基礎，得到了組合腐蝕速率與溫度、流速的關系式。陳迪等[30]將溫度和流速2個因素引入總預測模型，以正交實驗結果確定待定系數(shù)，建立了不同氣體分壓、流速和溫度等工況條件下的含硫濕天然氣管道腐蝕速率預測模型，已通過正交實驗得到2種預測模型：

(6)

(7)

式中，vcorr是腐蝕速率，mm/a；PH2S是H2S分壓，MPa；PCO2為CO2分壓，MPa；R是氣體常數(shù)；T是溫度，K；v是液體流速，m/s。

H2S/CO2共存體系的腐蝕規(guī)律較復雜，為綜合考慮氣體分壓、溫度等因素，一些新方法也逐漸用于腐蝕預測，BP神經(jīng)網(wǎng)絡即是一種常用于預測H2S/CO2共存體系腐蝕速率的算法，在特定的腐蝕環(huán)境下，該模型的預測精度能夠與常規(guī)方法建立的模型相當，如果進一步與改進粒子群優(yōu)化過程結合，可能得到更高的預測精度[21-32]。還有粒子群算法、支持向量機等算法也逐漸用于腐蝕速率預測模型研究[33-35]。

H2S/CO2共存體系下的腐蝕速率預測模型相對較少，早期的模型中基本只明確考慮了分壓的影響，近期的研究逐漸將更多的影響因素考慮到模型中，并借助一些算法建立部分模型，取得了初步的成果。但是，H2S/CO2共存體系的腐蝕預測模型都只停留在經(jīng)驗和半經(jīng)驗階段，未能將腐蝕機理的研究成果充分運用到建立腐蝕預測模型中，而且大多預測模型與油氣田生產(chǎn)現(xiàn)場的結合度不夠。

3 結束語

CO2腐蝕預測相對較成熟，從最初僅考慮了溫度、CO2分壓2個因素，到之后通過引入流速、產(chǎn)物膜、油膜等關鍵因素校正，再到充分結合生產(chǎn)實際，在模型中體現(xiàn)流態(tài)、材料、產(chǎn)物膜的沉積規(guī)律甚至電化學反應機理等的影響，預測模型越來越完善，實用性不斷提高。不同模型的建模思路往往存在差異，適用條件也不同，應用時應選擇最符合生產(chǎn)現(xiàn)場實際情況的模型。H2S/CO2共存體系的預測模型研究起步較晚，模型中考慮的因素很少，比較有代表性的模型是根據(jù)2類氣體的分壓比建立的。部分學者將其它影響因素引入分壓模型進行補充或采用BP神經(jīng)網(wǎng)絡等新算法建立在某些腐蝕條件下更適用的新模型。但是，目前H2S/CO2共存體系的預測模型都是經(jīng)驗模型或半經(jīng)驗模型，未來還需更多地結合酸性氣田的實際情況，深入地分析各種腐蝕環(huán)境中的主要影響因素以及其作用規(guī)律，以獲得形式更簡單、考慮因素更全面、對特定腐蝕環(huán)境預測精度更高的腐蝕速率預測模型。