一种基于PCA-ABC-SVM模型的海底原油管道腐蚀速率预测方法

技术领域

本发明涉及海底油气管道输送技术领域，具体涉及一种基于PCA-ABC-SVM模型的海底原油管道腐蚀速率预测方法。

背景技术

随着我国海洋油气技术的不断发展，海底原油管道的数量不断增加。海底原油管道承担着海洋原油运输的重要任务，然而，由于恶劣的海洋环境，加之输送介质，服役年限等特点，管道腐蚀损伤长度不断增加，因腐蚀引起的管道穿孔、泄漏爆炸事故时有发生，造成大量人员伤亡、经济损失和严重的环境污染。统计显示，在海底管道事故类型占比中，由腐蚀造成的事故占比高达37％。腐蚀速率预测能够从统计角度提供可靠的管道腐蚀缺陷数据，是降低海底原油管道腐蚀事故发生概率，减缓腐蚀事故后果的有效手段。建立海底原油管道腐蚀速率预测模型，对保障管道长期平稳服役具有重要意义。

现有对于海底原油管道腐蚀速率的预测模型包括单一预测模型的预测方法和多种模型结合的预测方法。由于影响海底管道腐蚀的因素众多且关系错综复杂，且很多用于模型中的数据本身的获取难度较大，使得依赖固定部分数据的单一预测模型的其预测可靠性欠佳。

其次，多种模型结合的预测方法由于不同模型自身的优缺点，对于原始数据的利用是存在差异的，难以将其直接结合；且由于部分数据监测困难，在实际应用中常常难以获取，因此，简单的对多个预测模型进行整合，存在需要数据量较大，对数据采集要求较高，在实际应用中，常常难以满足的问题，因而也难以达到真正有效评估管道腐蚀速率的目的。

发明内容

本发明针对现有技术中的上述不足，提供了一种能够解决现有技术中海底原油管道腐蚀速率预测模型对基础数据要求量大且预测可靠性不足的问题的基于PCA-ABC-SVM模型的海底原油管道腐蚀速率预测方法。

为解决上述技术问题，本发明采用了下列技术方案：

提供了一种基于PCA-ABC-SVM模型的海底原油管道腐蚀速率预测方法，其包括如下步骤：包括如下步骤：

S1、获取待评估的海底原油管道包含实际腐蚀速率值在内的检测数据；

S2、根据获取的检测数据，通过PCA算法建立海底原油管道腐蚀指标体系；

S3、根据海底原油管道腐蚀指标体系，通过SVM算法和ABC算法得到训练完毕的腐蚀速率预测模型；

S4、将检测数据代入训练完毕的腐蚀速率预测模型中，得到腐蚀速率预测结果。

本发明提供的上述基于PCA-ABC-SVM模型的海底原油管道腐蚀速率预测方法的主要有益效果在于：

本发明针对海底原油管道，提出基于PCA-ABC-SVM模型的海底原油管道腐蚀速率预测方法，能够在较小样本数据下更加真实地反映海底原油管道腐蚀速率；相对于现有基于常规检测数据的预测方法所需要的至少数千组数据，本方法所需要的数据量显著降低。

实际应用中，能够提供大量数据的海底管道，其实际工作时间较长，在数据检测和获取方面有各种缺陷，因而实际上难以提供现有方法中需要的大量数据。而本方案对数据量要求较小，对数据种类的要求也更低，即使年份较早的海底管道，也能满足本方案的数据量要求，因而本方案的适用范围远大于现有技术。

且通过PCA算法对检测数据进行预先分析，将实际对腐蚀速率贡献较低的参数剔除，一方面，进一步减小了对数据量的要求，贴合实际应用中的需求；另一方面，由于PCA算法筛选后的数据都是对腐蚀速率有较高实际贡献的，其预测准确性也显著提高。

相对于现有技术中直接通过腐蚀速率值对未来腐蚀速率进行预测的方法，本方法获取的数据及其变化能更好的反映影响腐蚀速率的各因素各自的变化，由于腐蚀速率本身就是受大量因素影响的一个计算结果，因而直接依靠腐蚀速率值进行预测的方法，其缺乏对影响腐蚀速率的实际各因素的评估，其在中长期预测的实际可靠性将显著低于本方案。

进一步地，本方案通过将人工蜂群算法，即ABC算法，与SVM算法结合，与传统的粒子群算法、遗传算法等群智能优化算法相比，全局寻优能力强，收敛速度较快，从而显著提高SVM算法的准确性和可靠性。

本发明可准确、快速的预测海底原油管道腐蚀速率，且可靠性较高，能为海底原油管道腐蚀失效风险预警提供科学依据与技术支持。

附图说明

图1是本方案提供的基于PCA-ABC-SVM模型的海底原油管道腐蚀速率预测方法的流程图。

图2是海底原油管道腐蚀速率预测方法与各算法模型的关系图。

图3是用于预测海底原油管道腐蚀速率的SVM算法的结构示意图。

图4是ABC算法优化海底原油管道腐蚀速率预测模型参数的流程图。

图5是SVM算法训练的模型对海底原油管道腐蚀速率的预测结果。

图6是PCA-ABC-SVM算法训练的模型对海底原油管道腐蚀速率的预测结果。

具体实施方式

下面结合附图对本发明作进一步说明：

如图1所示，其为本方案提供的基于PCA-ABC-SVM模型的海底原油管道腐蚀速率预测方法的结构示意图。

本发明的基于PCA-ABC-SVM模型的海底原油管道腐蚀速率预测方法包括如下步骤：

S1、获取待评估的海底原油管道包含实际腐蚀速率值在内的检测数据。

进一步地，检测数据至少包括管道内含水率、二氧化碳含量、硫化氢含量、氯离子含量、钙镁离子含量、溶解氧含量、pH值、温度、压力、流速和实际腐蚀速率值。

在本方案的一些实施例中，检测数据至少有100组。

由于实际应用中不同管道的数据采集难度不同，因此，能够获得的检测数据量也各不相同。理想情况下，检测数据越多，腐蚀速率的预测准确性越高。但通过反复验证，本方案要求的最低数据量为每一种数据至少包括100项。此时的预测精度即能满足一般预测要求。

而现有技术中的方法，通常数据需求量在300组以上，否则起可靠性和准确性就无法通过一般预测要求。

S2、根据获取的检测数据，通过PCA算法建立海底原油管道腐蚀指标体系。

如图2所示，海底原油管道腐蚀指标体系中的指标为基于所述检测数据进行主成分分析后，提取主成分分析得到的累计贡献率为设定值的至少两个指标。

具体的，通过主成分分析(principal component analysis，简称PCA)算法，对检测数据进行预先分析，将实际对腐蚀速率贡献较低的参数剔除，减小后续分析中的数据需求量，同时降低计算复杂度，提高分析效率。

进一步地，建立海底原油管道腐蚀指标体系的方法包括如下步骤：

S2-1、设检测数据的数据集为X＝{x

其中，x

S2-2、对原始矩阵进行标准化处理，建立标准化矩阵Y：

其中，

S2-3、根据标准化矩阵Y，计算其对应的相关阵R：

其中，y

S2-4、计算相关阵R的特征值和特征向量，并对将特征值降序排列处理，满足：

其中，λ

S2-5、根据特征值，计算其对应的贡献率与累计贡献率，其计算方法为：

式中，l

S2-6、提取累计贡献率前m个指标为主成分指标，根据主成分指标对应数据的样本数据集，建立海底原油管道腐蚀指标体系。

贡献率越高，其对腐蚀速率的影响越大。由于腐蚀速率值本身是受多个因素复合影响下的结果，一般不考虑其对腐蚀速率的影响。

优选的，当指标的累计贡献率达到85％～98％时，将该指标作为主成分指标，以保证分析的全面性。

S3、根据海底原油管道腐蚀指标体系，通过SVM算法和ABC算法得到训练完毕的腐蚀速率预测模型。

其中，支持向量机(support vector machines，简称SVM)算法和人工蜂群(Artificial Bee Colony，简称ABC)算法均为现有算法。

如图3所示，在本方案中，海底原油管道腐蚀速率模型由输入层、支持向量层、输出层组成。

其中，输入层为前述海底原油管道腐蚀指标体系中的检测数据；支持向量层为通过训练得到的内核函数k(x

一般的，海底原油管道腐蚀速率模型的原始模型为可通过SVM算法得到。

进一步地，得到训练完毕的腐蚀速率预测模型的方法包括如下步骤：

S3-1、通过SVM算法初步构建海底原油管道腐蚀速率预测模型。

S3-2、将海底原油管道腐蚀指标体系中的样本数据集划分为训练样本集和测试样本集。

进一步地，划分为训练样本集和测试样本集的方法为：

以海底原油管道腐蚀指标体系为数据基础进行随机排序处理；

将随机排序处理后的海底原油管道腐蚀指标体系中的数据的一部分作为训练样本集，剩下的部分作为测试样本集。

在本方案的一些实施例中，将随机排序处理后的海底原油管道腐蚀指标体系中的数据的80％作为训练样本集，20％作为测试样本集。以保证最大限度的处理

S3-3、通过训练样本集中的数据对海底原油管道腐蚀速率预测模型进行训练，并通过ABC算法对SVM算法的关键参数进行优化，得到优化训练后的海底原油管道腐蚀速率预测模型。

其中，SVM算法的参数选择对SVM算法训练后模型的预测精度有很大的影响。现有技术中，SVM算法使用网格搜索法或梯度下降法对其参数进行选择，然而，这些方法不仅耗时且导致非最优参数。本方案中，在通过SVM算法对海底原油管道腐蚀速率预测模型进行训练的同时，通过ABC算法对SVM算法的参数进行优化选择，从而能够比现有其他算法能更有效地得到SVM算法的参数的全局最优解，因而，将通过这一过程得到海底原油管道腐蚀速率预测模型，称为优化训练后的海底原油管道腐蚀速率预测模型。

进一步地，ABC算法即人工蜂群算法，其常规步骤如图4所示，其包括三个基本要素，分别是蜜源、被雇佣的蜜蜂及未被雇佣的蜜蜂。蜜源的位置代表一种优化方案，蜜源的含蜜量表示优化方案的适应度值；被雇佣的蜜蜂即雇佣蜂，任务是搜索新蜜源，并以一定概率与其他蜜蜂分享蜜源信息，数量与蜜源数量相等；未被雇佣的蜜蜂包括跟随蜂与侦察蜂两种，跟随蜂是在蜂巢周围等待雇佣蜂分享信息并寻找新蜜源的蜜蜂，在蜂群中数量与雇佣蜂相等；侦察蜂是在蜂巢附近随机搜索新蜜源的蜜蜂。

在本方案的一些实施例中，ABC算法分别对海底原油管道腐蚀速率预测模型中的惩罚因子C与核参数γ值进行优化。

这两个参数在常规SVM算法中，需要更多的样本数据，才能保证取值的可靠性。通过ABC算法，在显著减小样本数据量的同时，还能提高取值结果的可靠性。

进一步地，优化训练后的海底原油管道腐蚀速率预测模型中的参数的方法包括如下步骤：

A1、设置ABC算法优化框架，确定需要优化的模型参数，进一步设置迭代代数，计算极限值以及优化参数的搜索上界与搜索下界，随机初始化NP个方案，其计算方法为：

其中，x

A2、ABC算法中的雇佣蜂更新初始方案，其计算方法为：

其中，l≠i，rand[-1，1]为在[-1，1]范围内的随机数。

计算方案的适应度值，其计算方法为：

其中，fitness

随后，比较初始方案与更新方案的适应度值，若更新方案的适应度值优于初始方案，则基于贪婪法储存更新方案，否则，保留初始方案。

A3、在所有雇佣蜂完成搜索后，根据雇佣蜂分享的蜜源信息，ABC算法中的跟随蜂计算所得的概率P

根据概率P

A4、若存在某个方案在设置的极限值循环内多次迭代适应度值没有改进，则放弃该方案，雇佣蜂转化为侦察蜂，重复步骤A1至步骤A3，直至得到最优方案，完成ABC算法。

在重复步骤A1至步骤A3的过程中得到的最高适应度值所对应的方案，即最优方案。

A5、ABC算法结束之后，将获取的惩罚因子C与核参数γ值耦合入SVM模型，完成对训练后的海底原油管道腐蚀速率预测模型中的参数的优化。

在实际操作中，也可以根据需要增加需要ABC算法优化的参数，以保证预测结果的可靠性和准确性。

优选的，在训练模型之前，需要对训练样本集和测试样本集中的样本数据进行归一化处理，从而避免奇异样本数据引起的不良误差影响。

S3-4、将测试样本集中的数据代入到优化后的腐蚀速率预测模型中，得到包含腐蚀速率值的预测结果。

S3-5、计算预测结果的误差。

具体的，计算预测结果的误差的方法为：

分别计算预测结果与测试样本集中的实际腐蚀速率值的平均绝对误差、均方根误差和决定系数，其中，平均绝对误差、均方根误差和决定系数的计算方法分别为：

其中，y

S3-6、判断所述误差是否满足对应的设定条件，当误差不满足对应的设定条件时，重复对海底原油管道腐蚀速率预测模型的优化训练，直至误差满足对应的设定条件。

S4、将检测数据代入训练完毕的腐蚀速率预测模型中，得到腐蚀速率预测结果。

通过将检测数据代入训练完毕的腐蚀速率预测模型中，预测结果与实际腐蚀速率测量值符合性较好，说明训练完毕的腐蚀速率预测模型中预测能力较强，准确性较好。本方案分析速度快，能为海底原油管道腐蚀失效风险预警提供科学依据与技术支持。

下面是基于本发明提供的海底管道原油腐蚀速率预测方法的实施例：

本实施例基于某区域海底管道的实际采集数据，全部数据共107组，剔除其中数据显著异常和欠缺的7组，实际用于预测的数据共100组，下面基于前述海底管道原油腐蚀速率预测方法的步骤，进行描述：

S1、获取待评估的海底原油管道包含实际腐蚀速率值在内的检测数据。

本实施例中监测到的海底原油管道腐蚀速率相关的影响因素包括：含水率(X

表1海底原油管道包含实际腐蚀速率值在内的检测数据

S2、根据获取的检测数据，通过PCA算法建立海底原油管道腐蚀指标体系。

由于腐蚀速率值本身是受多个因素复合影响下的结果，一般不考虑其对腐蚀速率的影响。表2所示为PCA算法计算所得的累计贡献率。

表2主成分累计贡献率

一般累计贡献率的取值为85％～98％，以保证分析的全面性。在本实施例中，取累计贡献率达到95％的前四项主成分作为海底原油管道腐蚀指标体系，以在分析全面性的同时保证计算效率。

S3、根据海底原油管道腐蚀指标体系，通过SVM算法和ABC算法得到训练完毕的腐蚀速率预测模型。

其中，腐蚀速率预测模型可以通过SVM算法建立。

下面是通过本方案训练完毕的腐蚀速率预测模型的结果误差，其中MAE和RMSE的设定要求均为小于10％，R

表3训练完毕的腐蚀速率预测模型的预测误差

由表3可知，本方案通过PCA算法、ABC算法与SVM算法的结合，得到腐蚀速率预测模型训练能满组全部设定要求。

S4、将检测数据代入训练完毕的腐蚀速率预测模型中，得到腐蚀速率预测结果。

如图5和图6所示，通过本方案PCA-ABC-SVM算法训练得到的腐蚀速率预测模型，相对于现有仅通过SVM算法预测的模型，其在结果准确性上有非常显著的优势。

上面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。