CN111177947B - 一种考虑多因素的co2腐蚀预测图版建立方法 - Google Patents

Abstract

本发明公开了一种考虑多因素的CO₂腐蚀预测图版建立方法，属于油气田腐蚀防护技术领域。所述方法包括：首先在不同温度、CO₂分压和氯离子浓度下，采用高温高压釜对钢材进行腐蚀失重测试得到腐蚀速率实测值；然后利用非线性拟合方法建立腐蚀速率与温度、CO₂分压和氯离子浓度的函数关系式和关系曲线图；最后综合上述得到的腐蚀速率实测值、函数关系式和关系曲线图绘制CO₂腐蚀预测图版。该方法考虑多种因素对CO₂腐蚀的影响，具有较高的预测精度和实用性，为预防油田金属设施的腐蚀失效提供了参考依据，为材质的安全使用提供了技术保障。

Description

一种考虑多因素的CO2腐蚀预测图版建立方法

技术领域

本发明属于油气田腐蚀防护技术领域，具体涉及一种考虑多因素的CO₂腐蚀预测图版建立方法。

背景技术

油气田开发过程中，CO₂气体对井筒内管材和地面设备会造成严重的腐蚀，导致金属设备发生腐蚀失效，甚至引发安全事故。室内模拟腐蚀是评估材料耐腐蚀性能的主要手段，但如何通过室内实验的有限测试数据得到材质的腐蚀规律是预判材质适用性和预防设备失效的关键，为动态评估设备及石油管的腐蚀状况，需绘制可靠的CO₂腐蚀预测图版对金属管材进行腐蚀寿命预测。

目前，针对油气田金属设施的腐蚀预测，中石油的赵永涛等人提出了“一种BP神经网络的碳钢腐蚀速率预测方法”，该预测方法结合试验数据建立BP神经网络模型，然后利用Matlab编程对BP神经网络进行训练和测试，最后再通过测试的BP神经网络模型对碳钢腐蚀速率进行预测。然而，该方法基于大量的实验测试数据，试验周期长，费用成本高昂。DW模型考虑了温度、压力、pH值和流速等腐蚀因素，当温度低于85℃时，预测结果与试验结果吻合得较好，此模型适用于预测低温环境下的腐蚀速率，但在高温和高pH值下与试验结果相差较大。NORSOK模型考虑了温度、CO₂分压、pH值和剪切力等因素，在不发生局部腐蚀的条件下，此模型适用于预测二氧化碳腐蚀环境下的最大腐蚀速率，其中CO₂分压和pH值在现场不易测得，并且温度有效输入变量的范围为20℃～150℃，CO₂分压有效输入变量的范围为0.01MPa～1MPa，然而现场腐蚀工况温度有时会高达200℃左右，CO₂分压有时会大于1MPa，使用范围相对较窄，这将无法满足现场腐蚀预测的要求。可见现有模型的应用范围具有一定的局限性，未综合考虑到温度，CO₂分压，氯离子浓度腐蚀因素带来的影响。

因此，绘制可靠的CO₂腐蚀速率预测图版，发明一种基于典型工况实验数据且能考虑多种腐蚀因素的CO₂腐蚀预测图版的建立方法显得十分必要。

发明内容

本发明的目的在于提供一种考虑多因素的CO₂腐蚀预测图版建立方法，以解决现有的CO₂腐蚀预测的问题，并在达到上述目的的同时，为油气田特定环境下的CO₂腐蚀预测提供依据。一种考虑多因素的CO₂腐蚀预测图版建立方法，其特征在于，包括如下步骤：

步骤1：根据CO₂腐蚀实验工况确定各项影响因素参数范围，其中温度参数为T₁、T₂…T_i…T_n,T_1<T_2<T_i<T_n，CO₂分压参数为

氯离子浓度参数为Ccl^- ₁、Ccl^- ₂…Ccl^- _i…Ccl^- _n，Ccl^- _1<Ccl^- _2<Ccl^- _i<Ccl^- _n)；

其中：n为实验选定的参数个数3～6个；T_i为参数范围内的某一个温度点，℃；

为参数范围内的某一个CO₂分压点，MPa；Ccl^- _i为参数范围内的某一个氯离子浓度点，ppm；

步骤2：在上述确定的参数范围内对同一材质进行腐蚀失重测试，建立腐蚀速率与腐蚀影响因素包括温度、CO₂分压和氯离子浓度的函数关系式和关系曲线图，具体步骤如下：

(a)对同一材质在相同CO₂分压

氯离子浓度Ccl^- _i和不同温度T₁、T₂…T_i…T_n下进行失重法腐蚀测试得到实测腐蚀速率值为v_T1、v_T2…v_Ti…v_Tn；

(b)对同一材质温度T_i、氯离子浓度Ccl^- _i和不同CO₂分压

下进行失重法腐蚀测试得到实测腐蚀速率值为v_P1、v_P2…v_pi…v_Pn；

(c)对同一材质在相同温度T_i、CO₂分压

和不同氯离子浓度Ccl^- ₁、Ccl^- ₂…Ccl^- _i…Ccl^- _n下进行失重法腐蚀测试得到实测腐蚀速率值为v_C1、v_C2…v_Ci…v_Cn；

其中：v_T1、v_T2…v_Ti…v_Tn为T₁、T₂…T_i…T_n温度下对应的腐蚀速率实测值，mm/a；v_P1、v_P2…v_pi…v_Pn为

分压下对应的腐蚀速率实测值，mm/a；v_C1、v_C2…v_ci…v_Cn为Ccl^- ₁、Ccl^- ₂…Ccl^- _i…Ccl^- _n氯离子浓度下对应的腐蚀速率实测值，mm/a；

(d)根据上述所得到的腐蚀速率值作出散点图，确定函数类型，利用非线性拟合方法得到材质的腐蚀速率与温度、CO₂分压和氯离子浓度的函数关系式(1)-(3)如下，并且相应拟合得到关系曲线图；

ν_Tcor＝f(T) (1)

ν_Ccor＝f(Ccl^-) (3)

f(T)式中：v_Tcor为不同温度下的腐蚀速率预测值，mm/a；f(T)为温度与腐蚀速率的关系函数；T为温度，℃；v_Pcor为不同腐蚀介质分压下的腐蚀速率预测值，mm/a；

为CO₂分压与腐蚀速率的关系函数；PCO₂为腐蚀介质分压，MPa；v_Ccor为不同氯离子浓度下的腐蚀速率预测值，mm/a；f(Ccl^-)为氯离子浓度与腐蚀速率的关系函数；Ccl^-为氯离子浓度，ppm；步骤3：根据腐蚀失重测试所得到的腐蚀速率值和步骤2所得到的函数关系式和关系曲线图，绘制CO₂腐蚀预测图版，具体过程如下：

(a)保证材质腐蚀影响因素氯离子浓度Ccl^- _i参数不变，在温度T₁、T₂…T_i…T_n和CO₂分压

参数变化范围内，按照一定的步长选取变量为t₁、t₂…t_w…t_m，t₁<t₂<t_w<t_m，带入步骤2拟合得到的腐蚀速率与温度的函数关系式(1)求得相对应的腐蚀速率值v_T1cor、v_T2cor…v_Twcor…v_Tmcor，同样地，按照一定的步长选取变量为

带入步骤2拟合得到的腐蚀速率与CO₂分压的函数关系式(2)求得相对应的腐蚀速率值v_P1cor、v_P2cor…v_pwcor…v_Pmcor，失重测试得到的腐蚀速率值v_T1、v_T2…v_Ti…v_Tn和v_P1、v_P2…v_pi…v_Pn，利用非线性拟合方法绘制腐蚀速率与温度和CO₂分压的CO₂腐蚀预测图版；

(b)保证材质腐蚀影响因素CO₂分压

参数不变，在温度T₁、T₂…T_i…T_n和氯离子浓度Ccl^- ₁、Ccl^- ₂…Ccl^- _i…Ccl^- _n参数变化范围内，按照一定的步长选取变量为c_cl-1、c_cl-2…c_cl-w…c_cl-m，c_cl-1<c_cl-2<c_cl-w<c_cl-m，带入步骤2拟合得到的腐蚀速率与氯离子浓度的函数关系(3)求得相对应的腐蚀速率值v_C1cor、v_C2cor…v_Cwcor…v_Cmcor，结合步骤(a)中t₁、t₂…t_w…t_m对应的腐蚀速率值v_T1cor、v_T2cor…v_Twcor…v_Tmcor、失重测试得到的腐蚀速率值v_T1、v_T2…v_Ti…v_Tn和v_C1、v_C2…v_ci…v_Cn，利用非线性拟合方法绘制腐蚀速率与温度和氯离子浓度的CO₂腐蚀预测图版；

(c)保证材质腐蚀影响因素温度T_i参数不变，在CO₂分压

和氯离子浓度Ccl^- ₁、Ccl^- ₂…Ccl^- _i…Ccl^- _n参数变化范围内，结合步骤(a)中

对应的腐蚀速率值v_P1cor、v_P2cor…v_pwcor…v_Pmcor、步骤(b)中c_cl-1、c_cl-2…c_cl-w…c_cl-m对应的腐蚀速率值v_C1cor、v_C2cor…v_Cwcor…v_Cmcor、失重测试得到的腐蚀速率值v_P1、v_P2…v_pi…v_Pn和v_C1、v_C2…v_ci…v_Cn，利用非线性拟合方法绘制中腐蚀速率与CO₂分压和氯离子浓度的CO₂腐蚀预测图版；

其中：m为划分的步长个数10～50个；t_w为选取的某一个温度点，℃；

为选取的某一个CO₂分压点，MPa；c_cl-w为选取的某一个氯离子浓度点，ppm；v_T1cor、v_T2cor…v_Twcor…v_Tmcor为t₁、t₂…t_w…t_m对应的腐蚀速率预测值，mm/a；v_P1cor、v_P2cor…v_pwcor…v_Pmcor为

对应的腐蚀速率预测值，mm/a；v_C1cor、v_C2cor…v_Cwcor…v_Cmcor为c_cl-1、c_cl-2…c_cl-w…c_cl-m对应的腐蚀速率预测值，mm/a。

本发明具有以下有益效果：

(1)该腐蚀速率预测图版基于部分典型工况的室内实验数据，所需数据易得，数据量不大，可节约测试成本；

(2)该腐蚀速率预测图版综合考虑温度、CO₂分压和氯离子浓度多种腐蚀因素的影响，具有较高的适用性，可克服经典模型考虑因素不全面的缺陷，可突破腐蚀影响因素有效输入范围较窄的局限性。

(3)该腐蚀速率预测图版基于模拟工况实验数据，能够较为准确预测腐蚀情况，可为油气田金属设施的腐蚀预测、防护和更换等提供技术依据。

附图说明

图1是一种考虑多因素的CO₂腐蚀预测图版绘制方法的流程图；

图2是D级杆腐蚀速率与温度的关系曲线；

图3是D级杆腐蚀速率与二氧化碳分压的关系曲线；

图4是D级杆腐蚀速率与氯离子浓度的关系曲线；

图5是氯离子浓度为1000ppm时，腐蚀速率与温度和二氧化碳分压的CO₂腐蚀预测图版；

图6是二氧化碳分压为0.1MPa时，腐蚀速率与温度和氯离子浓度的CO₂腐蚀预测图版；

图7是温度为90℃时，腐蚀速率与二氧化碳分压和氯离子浓度的CO₂腐蚀预测图版。

具体实施方式

下面结合附图和具体实施案例对本发明进行详细的描述。

实施案例：

模拟XX油田实际工况，对D级杆钢绘制腐蚀预测图版。该工况条件如下：井筒温度范围为60℃～120℃、CO₂分压范围为0.1MPa～0.3MPa、氯离子浓度范围为1000ppm～3000ppm，具体实施步骤如下所示：

步骤1：根据XX油田实际工况确定各项影响因素参数范围，其中温度参数为60℃、75℃、90℃、105℃、120℃，CO₂分压参数为0.1MPa、0.2MP、0.3MPa和氯离子浓度参数为1000ppm、2000ppm、3000ppm；

步骤2：在上述确定的参数范围内对D级杆进行腐蚀失重测试，建立腐蚀速率与腐蚀影响因素包括温度、CO₂分压和氯离子浓度的函数关系式和关系曲线图，具体步骤如下：

(a)对D级杆在相同CO₂分压0.1MPa、氯离子浓度1000ppm和不同温度60℃、75℃、90℃，105℃、120℃下进行失重法腐蚀测试得到实测腐蚀速率值0.22mm/a、0.60mm/a、0.94mm/a，0.68mm/a、0.38mm/a；

(b)对D级杆在相同温度90℃、氯离子浓度1000ppm和不同CO₂分压0.1MPa、0.2MP、0.3MPa下进行失重法腐蚀测试得到实测腐蚀速率值0.94mm/a、0.99mm/a、1.12mm/a；

(c)对D级杆在相同温度90℃、CO₂分压0.1MPa和不同氯离子浓度1000ppm、2000ppm、3000ppm下进行失重法腐蚀测试得到实测腐蚀速率值0.37mm/a、0.94mm/a、1.04mm/a；

(d)根据上述所得到的腐蚀速率值作出散点图，确定函数类型，利用非线性拟合方法得到材质的腐蚀速率与温度、CO₂分压和氯离子浓度的函数关系式如式(1)-(3)所示，并且相应拟合得到关系曲线图如图2-4所示；

式中：v_Tcor为不同温度下的腐蚀速率预测值，mm/a；T为温度，℃；v_Pcor为不同腐蚀介质分压下的腐蚀速率预测值，mm/a；P_CO2为腐蚀介质分压，MPa；v_Ccor为不同氯离子浓度下的腐蚀速率预测值，mm/a；Ccl^-为氯离子浓度，ppm；

步骤3：根据腐蚀失重测试所得到的腐蚀速率值和步骤2所得到的函数关系式和关系曲线图，绘制CO₂腐蚀预测图版，具体过程如下：

(a)保证D级杆腐蚀影响因素氯离子浓度1000ppm参数不变，在温度60℃～120℃和CO₂分压0.1MPa～0.3MPa参数变化范围内，按照一定的步长选取温度变量64℃、68℃、72℃、76℃、80℃、84℃、88℃、92℃、96℃、100℃、104℃、108℃、112℃、116℃,带入步骤2拟合得到的腐蚀速率与温度的函数关系式(1)求得相对应的腐蚀速率值0.32mm/a、0.42mm/a、0.53mm/a、0.64mm/a、0.74mm/a、0.83mm/a、0.88mm/a、0.91mm/a、0.89mm/a、0.84mm/a、0.76mm/a、0.66mm/a、0.55mm/a、0.44mm/a，同样地，按照一定的步长选取CO₂分压变量0.1135MPa、0.127MPa、0.1405MPa、0.154MPa、0.1675MPa、0.181MPa、0.1945MPa、0.208MPa、0.2215MPa、0.235MPa、0.2485MPa、0.262MPa、0.2755MPa、0.289MPa,带入步骤2拟合得到的腐蚀速率与CO₂分压的函数关系式(2)求得相对应的腐蚀速率值0.94mm/a、0.94mm/a、0.94mm/a、0.95mm/a、0.96mm/a、0.97mm/a、0.98mm/a、0.99mm/a、1.01mm/a、1.02mm/a、1.04mm/a、1.06mm/a、1.08mm/a、1.10mm/a、失重测试得到的温度对应的腐蚀速率值0.22mm/a、0.60mm/a、0.94mm/a、0.68mm/a、0.38mm/a和CO₂分压对应的腐蚀速率值0.94mm/a、0.99mm/a、1.12mm/a，利用非线性拟合方法绘制腐蚀速率与温度和CO₂分压的CO₂腐蚀预测图版如图5所示；

(b)保证材质腐蚀影响因素CO₂分压0.1MPa参数不变，在温度60℃～120℃和氯离子浓度1000ppm～3000ppm参数变化范围内，按照一定的步长选取氯离子浓度变量1142ppm、1282ppm、1426ppm、1568ppm、1710ppm、1852ppm、1994ppm、2136ppm、2278ppm、2420ppm、2562ppm、2704ppm、2846ppm、2988ppm,带入步骤2拟合得到的腐蚀速率与氯离子浓度的函数关系式(3)求得相对应的腐蚀速率值0.44mm/a、0.52mm/a、0.61mm/a、0.69mm/a、0.78mm/a、0.86mm/a、0.94mm/a、1.00mm/a、1.05mm/a、1.09mm/a、1.11mm/a、1.10mm/a、1.08mm/a、1.05mm/a，结合步骤(a)中64℃、68℃、72℃、76℃、80℃、84℃、88℃、92℃、96℃、100℃、104℃、108℃、112℃、116℃对应的腐蚀速率值0.32mm/a、0.42mm/a、0.53mm/a、0.64mm/a、0.74mm/a、0.83mm/a、0.88mm/a、0.91mm/a、0.89mm/a、0.84mm/a、0.76mm/a、0.66mm/a、0.55mm/a、0.44mm/a，失重测试得到的温度对应的腐蚀速率值0.22mm/a、0.60mm/a、0.94mm/a、0.68mm/a、0.38mm/a和氯离子浓度对应的腐蚀速率值0.37mm/a、0.94mm/a、1.04mm/a，利用非线性拟合方法绘制腐蚀速率与温度和氯离子浓度的CO₂腐蚀预测图版如图6所示；

(c)保证材质腐蚀影响因素温度90℃参数不变，在CO₂分压0.1MPa～0.3MPa和氯离子浓度1000ppm～3000ppm参数变化范围内，结合步骤(a)中0.1135MPa、0.127MPa、0.1405MPa、0.154MPa、0.1675MPa、0.181MPa、0.1945MPa、0.208MPa、0.2215MPa、0.235MPa、0.2485MPa、0.262MPa、0.2755MPa、0.289MPa对应的腐蚀速率值0.94mm/a、0.94mm/a、0.94mm/a、0.95mm/a、0.96mm/a、0.97mm/a、0.98mm/a、0.99mm/a、1.01mm/a、1.02mm/a、1.04mm/a、1.06mm/a、1.08mm/a、1.10mm/a，步骤(b)中1142ppm、1282ppm、1426ppm、1568ppm、1710ppm、1852ppm、1994ppm、2136ppm、2278ppm、2420ppm、2562ppm、2704ppm、2846ppm、2988ppm对应的腐蚀速率值0.44mm/a、0.52mm/a、0.61mm/a、0.69mm/a、0.78mm/a、0.86mm/a、0.94mm/a、1.00mm/a、1.05mm/a、1.09mm/a、1.11mm/a、1.10mm/a、1.08mm/a、1.05mm/a、失重测试得到CO₂分压对应的腐蚀速率值0.94mm/a、0.99mm/a、1.12mm/a和氯离子浓度对应的腐蚀速率值0.37mm/a、0.94mm/a、1.04mm/a，利用非线性拟合方法绘制中腐蚀速率与CO₂分压和氯离子浓度的CO₂腐蚀预测图版如图7所示；

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (1)

1.一种考虑多因素的CO₂腐蚀预测图版建立方法，其特征在于，包括如下步骤：

步骤1：根据CO₂腐蚀实验工况确定各项影响因素参数范围，其中温度参数为T₁、T₂…T_i…T_n，T_1<T_2<T_i<T_n，CO₂分压参数为

氯离子浓度参数为Ccl^- ₁、Ccl^- ₂…Ccl^- _i…Ccl^- _n，Ccl^- _1<Ccl^- _2<Ccl^- _i<Ccl^- _n；

其中：n为实验选定的参数个数3～6个；T_i为参数范围内的某一个温度点，℃；

为参数范围内的某一个CO₂分压点，MPa；Ccl^- _i为参数范围内的某一个氯离子浓度点，ppm；

步骤2：在上述确定的参数范围内对同一材质进行腐蚀失重测试，建立腐蚀速率与腐蚀影响因素包括温度、CO₂分压和氯离子浓度的函数关系式和关系曲线图，具体步骤如下：

(a)对同一材质在相同CO₂分压

氯离子浓度Ccl^- _i和不同温度T₁、T₂…T_i…T_n下进行失重法腐蚀测试得到实测腐蚀速率值为v_T1、v_T2…v_Ti…v_Tn；

(b)对同一材质温度T_i、氯离子浓度Ccl^- _i和不同CO₂分压

下进行失重法腐蚀测试得到实测腐蚀速率值为v_P1、v_P2…v_pi…v_Pn；

(c)对同一材质在相同温度T_i、CO₂分压

和不同氯离子浓度Ccl^- ₁、Ccl^- ₂…Ccl^- _i…Ccl^- _n下进行失重法腐蚀测试得到实测腐蚀速率值为v_C1、v_C2…v_Ci…v_Cn；

其中：v_T1、v_T2…v_Ti…v_Tn为T₁、T₂…T_i…T_n温度下对应的腐蚀速率实测值，mm/a；v_P1、v_P2…v_pi…v_Pn为

分压下对应的腐蚀速率实测值，mm/a；v_C1、v_C2…v_ci…v_Cn为Ccl^- ₁、Ccl^- ₂…Ccl^- _i…Ccl^- _n氯离子浓度下对应的腐蚀速率实测值，mm/a；

(d)根据上述所得到的腐蚀速率值作出散点图，确定函数类型，利用非线性拟合方法得到材质的腐蚀速率与温度、CO₂分压和氯离子浓度的函数关系式(1)-(3)如下，并且相应拟合得到关系曲线图；

ν_Tcor＝f(T) (1)

ν_Ccor＝f(Ccl^-) (3)

f(T)式中：v_Tcor为不同温度下的腐蚀速率预测值，mm/a；f(T)为温度与腐蚀速率的关系函数；T为温度，℃；v_Pcor为不同腐蚀介质分压下的腐蚀速率预测值，mm/a；

为CO₂分压与腐蚀速率的关系函数；

为腐蚀介质分压，MPa；v_Ccor为不同氯离子浓度下的腐蚀速率预测值，mm/a；f(Ccl^-)为氯离子浓度与腐蚀速率的关系函数；Ccl^-为氯离子浓度，ppm；步骤3：根据腐蚀失重测试所得到的腐蚀速率值和步骤2所得到的函数关系式和关系曲线图，绘制CO₂腐蚀预测图版，具体过程如下：

(a)保证材质腐蚀影响因素氯离子浓度Ccl^- _i参数不变，在温度T₁、T₂…T_i…T_n和CO₂分压

参数变化范围内，按照一定的步长选取变量为t₁、t₂…t_w…t_m，t₁<t₂<t_w<t_m，带入步骤2拟合得到的腐蚀速率与温度的函数关系式(1)求得相对应的腐蚀速率值v_T1cor、v_T2cor…v_Twcor…v_Tmcor，同样地，按照一定的步长选取变量

带入步骤2拟合得到的腐蚀速率与CO₂分压的函数关系式(2)求得相对应的腐蚀速率值v_P1cor、v_P2cor…v_pwcor…v_Pmcor，失重测试得到的腐蚀速率值v_T1、v_T2…v_Ti…v_Tn和v_P1、v_P2…v_pi…v_Pn，利用非线性拟合方法绘制腐蚀速率与温度和CO₂分压的CO₂腐蚀预测图版；

(b)保证材质腐蚀影响因素CO₂分压P_CO2i参数不变，在温度T₁、T₂…T_i…T_n和氯离子浓度Ccl^- ₁、Ccl^- ₂…Ccl^- _i…Ccl^- _n参数变化范围内，按照一定的步长选取变量为c_cl-1、c_cl-2…c_cl-w…c_cl-m，c_cl-1＜c_cl-2＜c_cl-w＜c_cl-m，带入步骤2拟合得到的腐蚀速率与氯离子浓度的函数关系式；(3)求得相对应的腐蚀速率值v_C1cor、v_C2cor…v_Cwcor…v_Cmcor，结合步骤(a)中t₁、t₂…t_w…t_m对应的腐蚀速率值v_T1cor、v_T2cor…v_Twcor…v_Tmcor、失重测试得到的腐蚀速率值v_T1、v_T2…v_Ti…v_Tn和v_C1、v_C2…v_ci…v_Cn，利用非线性拟合方法绘制腐蚀速率与温度和氯离子浓度的CO₂腐蚀预测图版；

(c)保证材质腐蚀影响因素温度T_i参数不变，在CO₂分压

和氯离子浓度Ccl^- ₁、Ccl^- ₂…Ccl^- _i…Ccl^- _n参数变化范围内，结合步骤(a)中

对应的腐蚀速率值v_P1cor、v_P2cor…v_pwcor…v_Pmcor、步骤(b)中c_cl-1、c_cl-2…c_cl-w…c_cl-m对应的腐蚀速率值v_C1cor、v_C2cor…v_Cwcor…v_Cmcor、失重测试得到的腐蚀速率值v_P1、v_P2…v_pi…v_Pn和v_C1、v_C2…v_ci…v_Cn，利用非线性拟合方法绘制中腐蚀速率与CO₂分压和氯离子浓度的CO₂腐蚀预测图版；

其中：m为划分的步长个数10～50个；t_w为选取的某一个温度点，℃；

为选取的某一个CO₂分压点，MPa；c_cl-w为选取的某一个氯离子浓度点，ppm；v_T1cor、v_T2cor…v_Twcor…v_Tmcor为t₁、t₂…t_w…t_m对应的腐蚀速率预测值，mm/a；v_P1cor、v_P2cor…v_pwcor…v_Pmcor为

对应的腐蚀速率预测值，mm/a；v_C1cor、v_C2cor…v_Cwcor…v_Cmcor为c_cl-1、c_cl-2…c_cl-w…c_cl-m对应的腐蚀速率预测值，mm/a。

Images