含缺陷油氣管道在土壤環境下應力腐蝕數值模擬

摘 要 基于有限元數值模擬的方法，分析研究高強度管線鋼在應力作用下，不同土壤離子（Cl-、HCO3-、SO42-）濃度和不同溫度時的電化學腐蝕特征。結果表明：隨著離子濃度和溫度的增加，缺陷處的腐蝕電流密度均呈現增加的趨勢；在塑性應變階段，由于缺陷中心應力集中，機械-電化學效應顯著，缺陷中心位置相較于兩端腐蝕更加嚴重，中心位置的腐蝕電流密度顯著增加。

關鍵詞 油氣管道 管線鋼 應力腐蝕 機械-電化學效應 離子濃度 溫度

中圖分類號 TQ055.8+1" "文獻標志碼 A" "文章編號 0254?6094（2024）05?0723?05

基金項目：陜西延長石油（集團）有限責任公司項目（批準號：ycsy2023ky?B?65）資助的課題；山東省油氣儲運安全重點實驗室開放基金資助的課題；中央高校基本科研業務費專項（批準號：19CX05007A）資助的課題。

作者簡介：李騰飛（1986-），助理工程師，從事城市燃氣安全管理、油氣儲運系統安全管理工作，434888798@qq.com。

引用本文：李騰飛，張成斌，馬春迅，等.含缺陷油氣管道在土壤環境下應力腐蝕數值模擬［J］.化工機械，2024，51（5）：723-727.

管道是五大運輸方式與九大基礎設施網絡的重要組成部分，管道運輸具有良好的安全性和經濟性，現已成為油氣能源主要的輸送方式。截至2022年底，我國管道總里程達18.5×104 km，初步形成橫跨東西、縱貫南北、覆蓋全國、聯通海外的油氣“全國一張網”［1］。隨著人們對能源需求量的日益攀升，長輸油氣管道朝著大口徑、高強度和高壓力的方向發展［2］。埋地管道主要采用防腐層和外加電流陰極保護聯合防護［3～5］。然而，在復雜的土壤環境中，由于防腐層破損老化、土壤環境物化性質和電化學腐蝕綜合因素的作用，管材在遠低于材料屈服強度下發生應力腐蝕開裂，造成管道失效［6，7］。張秀云等主要通過電化學實驗研究X100管線鋼在鹽漬土壤模擬溶液中的腐蝕行為，發現不同質量分數的SO42-下，陽極均表現出活性溶解，說明SO42-對腐蝕的發生起到促進作用［8，9］。王彬彬等利用電化學工作站測得電化學極化曲線和阻抗譜，研究X70鋼在不同Cl-濃度的酸性紅壤中的腐蝕行為，結果表明，低Cl-濃度下發生局部腐蝕，且隨著濃度的增加點蝕坑數量和尺寸增加［10］。楊霜等利用電化學阻抗技術，研究了X80鋼在不同溫度的酸性模擬紅壤溶液中的腐蝕性能，結果表明土壤電阻和電荷轉移電阻隨溫度的升高而降低［11］。梁錕等利用電化學實驗測試技術，研究不同溫度下X90鋼母材及其焊縫在近中性土壤中的腐蝕行為，結果表明，隨溫度的增加，腐蝕電位負移、腐蝕傾向增加，腐蝕電流密度增大，腐蝕速率加快，相同情況下焊縫腐蝕比母材更為嚴重［12］。WANG Y X等研究不同應變下X80鋼在0.62 mol/L的NaCl溶液中的腐蝕行為，隨著應變增加，腐蝕活性的增加并不是線性變化的，而是呈現先增加后減小的趨勢［13］。宮克等研究X90管線鋼在沈陽土壤環境中不同HCO3-濃度下的應力腐蝕特性，研究表明，HCO3-主要通過改變鈍化膜的保護特性影響管道腐蝕，腐蝕速率隨著HCO3-濃度的增加呈現先增大后減小的趨勢，HCO3-濃度為7%時腐蝕速率達到峰值［14］。然而，當前高強度管線鋼土壤應力腐蝕大多數局限在電化學腐蝕的實驗研究中［15～17］，缺少對機械-電化學協同效應的綜合考慮，有關管線鋼應力腐蝕的數值模擬研究鮮有報道。筆者應用COMSOL多物理場耦合軟件，考慮機械-電化學協同效應建立含缺陷管道有限元模型，分析高強度管線鋼在應力作用下，不同土壤離子（Cl-、HCO3-、SO42-）濃度和不同溫度時的電化學腐蝕特性。

1 含缺陷管道有限元模型

含缺陷管道外表面與土壤環境的有限元二維幾何模型如圖1所示。取X90管線鋼長200 mm，壁厚19.6 mm，管道缺陷長20 mm，深5.88 mm（管道壁厚的30%），腐蝕電解質溶液區域選擇高100 mm作為模擬區域。

同時，基于COMSOL有限元仿真軟件，采用三角形網格類型，實現數值模擬計算域的網格劃分，如圖2所示。

2 結果與分析

2.1 Cl-濃度對應力腐蝕的影響

在0.1%和0.3%的應變下，不同Cl-濃度（0.024 5、0.049 0、0.073 5、0.098 0 mol/L）時缺陷處的陽極腐蝕電流密度如圖3所示。從圖3a可知，在0.1%應變條件下，缺陷中心和缺陷兩端的腐蝕電流密度幾乎重合，差異很小，陽極腐蝕電流密度變化為4.765、5.170、5.744、6.070 μA/cm2，機械電化學效應不顯著。從圖3b的結果來看，在0.3%應變條件下，缺陷中心陽極腐蝕電流密度變化為7.408、8.063、8.686、10.312 μA/cm2，呈現不斷增加的趨勢。

Cl-是一種腐蝕性離子，隨著Cl-濃度的增加，土壤的侵蝕作用不斷加強，破壞金屬基體導致金屬發生陽極溶解，加快腐蝕進程［18］，腐蝕電流密度逐漸增大。0.1%應變下從缺陷中心到缺陷兩端比較平滑，說明此階段管道處于彈性變形階段，機械-電化學效應對陽極反應的影響效果不顯著；應變為0.3%時，缺陷中心腐蝕電流密度明顯大于缺陷兩端的，管道處于塑性變形階段，缺陷中心的應力較大，超過材料的屈服強度，故在缺陷中心處的機械-電化學效應更加顯著。

2.2 HCO3-濃度對腐蝕的影響

在0.1%和0.3%的應變下，不同HCO3-濃度（0.029、0.058、0.087、0.116 mol/L）時缺陷處陽極腐蝕電流密度如圖4所示。從圖4a可以看出，在0.1%應變條件下，缺陷中心和缺陷兩端的腐蝕電流密度相差很小，陽極腐蝕電流密度變化為9.851、14.000、18.466、19.548 μA/cm2，但隨著HCO3-濃度的增加，侵蝕作用不斷加強，加速腐蝕速率進程，腐蝕電流密度逐漸增大。從圖4b可以看出，在0.3%應變條件下，缺陷中心處陽極腐蝕電流密度變化為14.014、19.354、23.973、25.461 μA/cm2，呈現不斷增加的趨勢。在0.1%應變下，從缺陷中心到缺陷兩端比較平滑，說明在此階段管道處于彈性變形階段，機械-電化學效應對陽極反應的影響效果不顯著。在0.3%應變下，缺陷中心腐蝕電流密度明顯大于兩端的腐蝕電流密度，說明此階段材料處于塑性變形，缺陷中心的應力較大，超過屈服強度，故在中心處的機械-電化學效應更加顯著。

2.3 SO42-濃度對應力腐蝕的影響

在0.1%和0.3%的應變下，不同SO42-濃度（0.002 5、0.005 0、0.007 5、0.010 0 mol/L）時缺陷處陽極腐蝕電流密度如圖5所示。從圖5a可看出，在0.1%應變條件下，缺陷中心和缺陷兩端的腐蝕電流密度變化不明顯，不同SO42-濃度下，陽極腐蝕電流密度為3.329、3.547、4.467、5.005 μA/cm2。隨著SO42-濃度的增加，陽極腐蝕電流密度呈現不斷增加的趨勢，侵蝕作用不斷加強，金屬基體穩定性變差，加快腐蝕速率［19］。從缺陷中心到缺陷兩端比較平滑，說明在此階段管道處于彈性變形階段，機械-電化學效應對陽極反應的效果不顯著。從圖5b可以看出，在0.3%應變條件下，缺陷中心處陽極腐蝕電流密度變化為5.244、5.926、7.566、8.312 μA/cm2，陽極腐蝕電流密度呈現不斷增加的趨勢。缺陷中心腐蝕電流密度明顯大于兩端的腐蝕電流密度，且隨SO42-濃度增加差距逐漸增大。說明當應變為0.3%時，管道處于塑性變形階段，缺陷中心的應力較大，超過材料的屈服強度，故在缺陷中心處的機械-電化學效應更加顯著。

2.4 溫度對腐蝕的影響

在0.1%和0.3%的應變下，不同溫度（25、40、50、60 ℃）時缺陷處陽極腐蝕電流密度如圖6所示。從圖6a的結果來看，在0.1%應變條件下，缺陷中心和缺陷兩端的腐蝕電流密度幾乎重合，陽極腐蝕電流密度變化為2.509、14.260、18.319、28.261 μA/cm2。從圖6b的結果來看，在0.3%應變條件下，缺陷中心處陽極腐蝕電流密度變化為4.228、19.549、27.287、37.967 μA/cm2，呈現不斷增加的趨勢。在0.1%應變下，從缺陷中心到缺陷兩端比較平滑，說明在此階段管道處于彈性變形階段，機械-電化學效應對陽極反應的效果并不顯著［20］。當應變為0.3%時，缺陷中心腐蝕電流密度明顯大于兩端的腐蝕電流密度，說明管道處于塑性變形階段，缺陷中心的應力較大，超過材料的屈服強度，故在缺陷中心處的機械-電化學效應更加顯著。隨著溫度的增加，陽極腐蝕電流密度呈現不斷增加的趨勢，溫度升高物質傳輸和擴散速度加快，溶液中離子擴散到電極表面的阻力減小［21］，擴散速率增加，電荷轉移電阻減小，破壞金屬基體導致金屬發生陽極溶解，加速腐蝕速率進程。

3 結束語

通過研究土壤離子（Cl-、HCO3-、SO42-）和溫度對管道缺陷處腐蝕的影響，發現隨著離子濃度和溫度增加，缺陷處腐蝕電流密度均呈現增加的趨勢，其中溫度比離子濃度對于腐蝕影響更加敏感，腐蝕速率增加更顯著。

管道在發生彈性應變和塑性應變階段對應力腐蝕影響程度不同。0.1%應變（彈性應變），在缺陷中心和兩端受到應力差異不顯著，腐蝕電流密度變化差異不顯著，此應變下缺陷處的機械-電化學效應不顯著。0.3%應變（塑性應變），在缺陷中心位置應力較大，機械-電化學效應更加顯著，中心處的腐蝕電流密度顯著增加。

參 考 文 獻

［1］" "陳朋超.油氣管網安全狀態監測傳感系統構建與創新發展［J］.油氣儲運，2023，42（9）：998-1008.

［2］" "劉宇，張立忠，高維新.管線鋼的歷史沿革及未來展望［J］.油氣儲運，2022，41（12）：1355-1362.

［3］" "李秋揚，趙明華，張斌，等.2020年全球油氣管道建設現狀及發展趨勢［J］.油氣儲運，2021，40（12）：1330-1337；1348.

［4］" "丁銳，姚寶慧，方孝斌.長輸地埋油氣管道腐蝕因素分析與防護對策探討［J］.應用化工，2019，48（12）：2972-2977.

［5］" "FENG Y C，CHENG Y F.An intelligent coating doped with inhibitor?encapsulated nanocontainers for corrosion protection of pipeline steel［J］.Chemical Engineering Journal，2017，315：537-551.

［6］" "LI X G，ZHANG D W，LIU Z Y，et al.Materials science：Share corrosion data［J］.Nature，2015，527：441-442.

［7］" "馬曉鳳，李媛媛，趙凱峰，等.埋地管道防腐保溫層下的電化學腐蝕行為［J］.西安石油大學學報（自然科學版），2021，36（3）：105-112.

［8］" "張秀云，石志強，王彥芳，等.X100管線鋼在鹽漬土壤模擬溶液中的腐蝕行為［J］.中國腐蝕與防護學報，2015，35（1）：33-37.

［9］" "MOAYED M H，NEWMAN R C.Aggressive effects of pitting “inhibitors” on highly alloyed stainless steels［J］.Corrosion Science，1998，40（2/3）：519-522.

［10］" "王彬彬，王軍，盧志偉，等.紅壤模擬溶液中Cl-對X80管線鋼腐蝕行為的影響［J］.全面腐蝕控制，2016，30（11）：54-58.

［11］" "楊霜，唐囡，閆茂成，等.溫度對X80管線鋼酸性紅壤腐蝕行為的影響［J］.中國腐蝕與防護學報，2015，35（3）：227-232.

［12］" "梁錕，雒設計，徐飛.溫度對X90管線鋼及其焊縫電化學腐蝕行為的影響［J］.腐蝕與防護，2020，41（8）：23-27；69.

［13］" "WANG Y X，ZHAO W M，AI H，et al.Effects of strain on the corrosion behaviour of X80 steel［J］.Corrosion Science，2011，53（9）：2761-2766.

［14］" "宮克，吳明，張勝.HCO3-對X90管線鋼應力腐蝕行為的影響［J］.中國腐蝕與防護學報，2021，41（5）：727-731.

［15］" LYNCH S P. Environmentally assisted cracking：Overview of evidence for an adsorption?induced localised?slip process［J］.Acta Metallurgica，1988，36（10）：2639-2661.

［16］" "FAYEZIOGHANI A，DEKKER R，SLUYS L J.Verification，validation，and parameter study of a computational model for corrosion pit growth adopting the level?set method.Part Ⅱ：Stress corrosion［J］.Materials Today Communications，2022，33：104210.

［17］" "JONES R H.Analysis of hydrogen?induced subcritical intergranular crack growth of iron and nickel［J］.Acta Metallurgica et Materialia，1990，38（9）：1703-1718.

［18］" "QIN G J，HUANG Y F，WANG Y H，et al.Pipeline condition assessment and finite element modeling of mechano?electrochemical interaction between corrosion defects with varied orientations on pipelines［J］.Tunnelling and Underground Space Technology，2023，136：105101.

［19］" "QIN G J， CHENG Y F， ZHANG P.Finite element modeling of corrosion defect growth and failure pressure prediction of pipelines［J］.International Journal of Pressure Vessels and Piping，2021，194：104509.

［20］" "CAO Y G，ZU Y Z，ZHEN Y，et al.Determination of the true stress?strain relations of high?grade pipeline steels based on small punch test correlation method［J］.International Journal of Pressure Vessels and Piping，2022，199：104739.

［21］" "XU L Y，CHENG Y F.Development of a finite element model for simulation and prediction of mechanoelectrochemical effect of pipeline corrosion［J］.Corrosion Science，2013，73：150-160.

（收稿日期：2023-11-17，修回日期：2024-09-11）

Numerical Simulation of Stress Corrosion of Oil and Gas

Pipelines with Defects in the Soil Environment

LI Teng?fei1， ZHANG Cheng?bin2， MA Chun?xun3， ZHANG Chen3，

ZHOU Yang?yang3， ZHANG Yu?hong3， BI Hai?sheng3

（1. PetroChina Jovo Natural Gas Co.， Ltd.； 2. Research Institute， Shaanxi Yanchang Petroleum （Group） Company Ltd .；

3. College of Electromechanical Engineering， Qingdao University of Science and Technology）

Abstract" "Having the FEM simulation method based to analyze electrochemical corrosion characteristics of high?strength steel pipeline stressed with different soil ion（Cl-， HCO3-， SO42-） concentrations and at different temperatures was implemented. The results show that， with the increase of both ion concentration and temperature， the corrosion current density at the defect tends to increase； in the plastic strain stage， the stress becomes concentrated in the defect center and the mechanical?electrochemical effect gets to be significant and the corrosion at the defect center becomes more seriously than that at both ends and the corrosion current density at the defect center becomes significantly increased.

Key words" "oil amp; gas pipeline， pipeline steel，stress corrosion， mechanical?electrochemical effect， ion concentration， temperature