基于响应面法的CO2驱替中N80钢腐蚀影响因素交互作用研究

《Materials Today Communications》：A Study on the Interaction of Corrosion Influencing Factors of N80 Steel in CO2 Flooding Based on Response Surface Methodology

　　

CO₂驱替技术的应用导致CO₂在水中溶解形成酸性溶液，从而引发严重的腐蚀挑战。本研究采用Design Expert软件，以CO₂压力、温度和矿化度为自变量，腐蚀速率为因变量，通过高温高压釜实验进行正交试验。利用响应面法（Response Surface Methodology, RSM）建立了腐蚀速率模型。通过方差分析（ANOVA）、响应曲面图和等高线图分析了影响腐蚀的关键因素和腐蚀模式。研究结果表明，该模型具有很高的可靠性，主要因素的影响顺序为：矿化度 > 压力 > 温度，并且这些因素之间存在交互作用。此外，研究人员选择了腐蚀速率最大、最小和中间的三个样品进行微观形貌、表面粗糙度和三维形貌分析。结果表明，氯离子容易引发点蚀。

**研究背景与问题**
CO₂驱替作为一种提高原油采收率的技术，在油田中广泛应用。然而，注入的CO₂在地层水中溶解会形成碳酸，导致井下管柱和集输管线面临严重的腐蚀问题，特别是对于广泛使用的N80钢。CO₂腐蚀受多种因素影响，包括温度、CO₂分压和矿化度（盐度）等。尽管已有大量研究探讨了单一因素对腐蚀的影响，但在CO₂驱替的实际工况下，这些因素往往同时存在并发生复杂的交互作用。目前的研究存在一些局限性：部分研究侧重于缓蚀剂优化而忽略了工艺参数的影响；一些研究主要针对非CO₂环境，缺乏对CO₂驱替条件下腐蚀的系统研究；另一些研究则忽视了高压和高矿化度的协同效应，未能阐明CO₂驱替条件下温度、压力与矿化度之间的交互作用机制。因此，系统研究这些关键操作参数之间的交互作用及其对N80钢腐蚀行为的影响，对于准确预测腐蚀风险、制定有效的防护策略至关重要。本研究旨在填补这一空白，论文发表在《Materials Today Communications》上。

**关键技术方法**
为开展研究，研究人员主要采用了以下关键技术方法：首先，利用响应面法中的中心复合设计（Central Composite Design, CCD）在Design Expert软件中设计了包含温度（70-90°C）、CO₂压力（3-5 MPa）和矿化度（35,000-55,000 mg/L）三个因素的五水平正交实验方案。其次，使用高温高压釜模拟实际工况，在设定的流速（1 m/s）下对N80钢试样进行为期72小时的动态失重实验。最后，综合运用多种表征技术对腐蚀产物和形貌进行分析，包括扫描电子显微镜（SEM）结合能谱仪（EDS）分析微观形貌和元素分布、X射线衍射（XRD）分析腐蚀产物物相、激光共聚焦显微镜（LSCM）进行三维形貌观察和表面粗糙度测量。

**研究结果**
**3.1. 基于RSM的统计分析与建模**
基于失重实验数据，通过方差分析建立了腐蚀速率（CR）的二次多项式预测模型。模型回归显著（P < 0.0001），失拟项不显著（P = 0.1030），表明模型拟合良好，可靠性高。方差分析显示，矿化度和CO₂压力对均匀腐蚀速率有极显著影响（P < 0.01），温度有显著影响（P < 0.05）。根据F值，各因素主效应的影响顺序为：矿化度 > 压力 > 温度。此外，矿化度-温度、压力-矿化度之间的交互作用对腐蚀速率有极显著影响，温度-压力交互作用有显著影响。贡献率分析表明，矿化度是腐蚀速率的主要控制因素，占总效应的52.20%。模型的确定系数R²为0.9944，调整后R²为0.9880，预测值与实际值之间的相对误差低于7%，验证了模型的准确性和稳健性。

**3.2. 各腐蚀影响因素的交互效应分析**
通过响应曲面图和等高线图分析了因素间的交互效应。分析表明，在恒定的矿化度和CO₂压力下，腐蚀速率随温度升高呈现先增加后降低的趋势。在70-90°C范围内，温度的升高一方面加速了离子传输和H⁺的传递，促进了腐蚀；另一方面也提高了腐蚀产物FeCO₃的沉积速率，有助于形成致密的保护膜，从而抑制腐蚀。在恒定的温度和矿化度下，腐蚀速率随CO₂压力升高而增加，这是因为更高的压力增加了CO₂的溶解度和溶液中的H⁺浓度，降低了pH值。在恒定的温度和CO₂压力下，腐蚀速率与矿化度（35,000-55,000 mg/L范围内）呈正相关，因为矿化度增加提高了溶液电导率，加速了电荷转移。等高线图的椭圆形表明各因素之间存在明显的交互作用。

**3.3. 结果验证**
利用回归模型确定了腐蚀速率最高的优化条件（温度83.6°C，CO₂压力5 MPa，矿化度55,000 mg/L），并进行了平行验证实验。实测腐蚀率为25.627 ± 0.706 mm/a，与模型预测值25.401 mm/a吻合良好（误差在5%以内），证实了响应面法优化结果的可靠性。

**3.4. 腐蚀产物膜分析**
**3.4.1. 宏观形貌**
不同腐蚀环境下N80钢的宏观形貌差异显著。部分样品（如2#和10#）表现为相对平坦的表面，以均匀腐蚀为主；而另一些样品（如8#）则显示出明显的局部腐蚀，存在深而宽的坑蚀（台地腐蚀），保护膜破坏严重。
**3.4.2. 微观形貌**
选取腐蚀速率最低（A，条件：100°C，4 MPa，45,000 mg/L）、中等（B，条件：80°C，4 MPa，55,000 mg/L）和最高（C，条件：90°C，5 MPa，55,000 mg/L）的三个样品进行SEM观察。低腐蚀速率样品A表面腐蚀产物附着松散，存在小孔。中、高腐蚀速率样品B和C表面则生成了立方晶形的FeCO₃产物，颗粒较大，堆积更致密。这表明在80°C左右，腐蚀产物膜的生长状态发生了转变，影响了其保护性能。
**3.4.3. 腐蚀产物膜厚度**
截面形貌显示，低腐蚀速率样品A的腐蚀产物膜存在孔隙和裂纹，而中、高腐蚀速率样品B和C的腐蚀膜更为致密，外表面有明显的腐蚀产物颗粒。这与高矿化度溶液（含Cl^-和HCO₃^-）能渗透膜层并与基体作用，沉积产物从而增强膜致密性的观点一致。
**3.4.4. 腐蚀产物分析**
EDS面扫描和线扫描分析表明，所有样品的腐蚀产物层均含有明显的O和Fe元素，以及钢中固有的Cr、Mn等元素。XRD分析进一步确认，中、高腐蚀速率样品的主要腐蚀产物为FeCO₃，而低腐蚀速率样品的FeCO₃衍射峰较弱，表明其结晶度较低、晶粒较小。
**3.4.5. 三维形貌与表面粗糙度**
去除腐蚀产物膜后，利用LSCM分析基体表面的点蚀情况。结果表明，腐蚀速率越高，点蚀坑越深，表面粗糙度范围越大。在腐蚀最严重的条件下（C），点蚀深度接近200 μm，计算得到的点蚀速率高达约29.46 mm/a。这揭示了在高矿化度（高Cl^-含量）和高CO₂压力环境下，局部腐蚀的严重性。

**3.5. 腐蚀机理**
基于以上分析，研究人员提出了氯离子诱导的局部腐蚀机理。首先，CO₂溶解生成H₂CO₃并电离出H⁺，与Fe反应生成FeCO₃沉淀，形成腐蚀产物膜。在适宜条件下，致密的FeCO₃膜能有效阻挡腐蚀性离子，降低均匀腐蚀速率。然而，在高矿化度或极端温度下，FeCO₃膜生长缓慢、结构疏松多孔，无法有效阻挡Cl^-等离子的侵入。Cl^-通过膜层缺陷在基体表面局部富集，破坏钝化膜，使该区域转变为活化态阳极，与周围大面积的阴极区构成“大阴极-小阳极”电偶腐蚀电池，引发点蚀萌生。点蚀一旦形成，坑内Fe²⁺水解产生H⁺，Cl^-为维持电中性向内迁移，形成自催化的酸性环境，加速金属溶解，导致蚀坑不断加深和扩大。

**讨论与结论**
**讨论部分**综合了实验结果，深入探讨了温度、CO₂压力和矿化度三个因素对N80钢CO₂腐蚀行为的单独及交互影响。研究证实了响应面法在量化多因素交互作用方面的有效性，并揭示了腐蚀产物膜的结构与保护性能随环境条件变化的规律。高矿化度环境中的Cl^-是诱发和加速局部腐蚀的关键因素，其与CO₂腐蚀的协同作用显著加剧了材料失效风险。研究结果对理解CO₂驱替恶劣工况下的腐蚀机制具有重要价值。
**结论部分**翻译如下：
(1) 利用响应面法建立了CO₂驱替腐蚀速率的预测模型。通过残差分析和方差分析验证了模型的准确性和稳定性。预测值与实际值之间的相对误差始终低于7%，为实际的管道腐蚀防护策略提供了有价值的见解。
(2) 腐蚀速率预测模型的方差分析结果表明，矿化度是N80钢腐蚀的主要控制因素。此外，CO₂压力与矿化度之间，以及矿化度与温度之间存在显著的交互作用。温度与CO₂压力之间的交互作用也很显著。
(3) N80试样腐蚀形貌的微观检查结果与基于响应面法得出的腐蚀模型相符，证实了该方法在研究腐蚀机理方面的实用性和有效性。