本发明属于埋地油气管道风险评估技术领域，涉及一种基于rs-pso-grnn的埋地油气管道土壤腐蚀预测方法。
背景技术：
近年来，随着我国地下管线铺设量增大，老管线服役时间过长，导致管道事故频发，埋地管道的安全性越来越受到人们广泛的关注。土壤腐蚀作为埋地管道面临失效风险的主要原因，由于其土壤成分繁多、结构复杂及腐蚀形式多样，使得准确评价管道的腐蚀状况并对管道进行安全性检测变得十分困难，同时管道所经某些地区，地广人稀，一旦发生泄漏，很难定位，可能引发连锁反应，造成极大的财产损失。因此，研究并掌握土壤的腐蚀规律和理化性质，建立具有普适性的土壤腐蚀预测新方法，以获取各监测点的管道年均腐蚀壁厚。
技术实现要素：
本发明的目的在于克服上述现有技术的缺点，提供了一种基于rs-pso-grnn的埋地油气管道土壤腐蚀预测方法，该方法能够准确获取各监测点的管道年平均腐蚀壁厚。
为达到上述目的，本发明所述的基于rs-pso-grnn的埋地油气管道土壤腐蚀预测方法包括以下步骤：
1)识别管道土壤腐蚀因素，再根据管道土壤腐蚀因素构建埋地管道土壤腐蚀指标体系，然后获取埋地管道土壤腐蚀的原始样本指标，并构建原始样本集{xij|i＝1,2,l,m,j＝1,2,l,n}，其中，xij为第i个管道样本的第j个指标值，m及n分别为监测点管道序号总数及管道土壤腐蚀指标的维数；
2)对步骤1)得到的原始样本指标值进行离散化处理，再根据离散化后的条件属性集及决策属性集建立决策表；
3)采用rs理论中的属性约简删除决策表中的冗余指标，得管道核心指标样本集；
4)对步骤3)得到的管道核心指标样本集进行归一化处理；
5)将步骤4)得到的归一化处理后的管道核心指标样本集分为训练样本集及测试样本集，再将训练样本作为grnn输入的训练集；
6)构建rs-pso-grnn土壤腐蚀预测模型，再将步骤5)得到的训练集输入至rs-pso-grnn土壤腐蚀预测模型中，得各监测点的管道年均腐蚀壁厚。
步骤4)中归一化处理的结果为：
其中，max(xj)及min(xj)分别为第j个指标的最大值及最小值，
步骤6)中的rs-pso-grnn土壤腐蚀预测模型以grnn输出的均方根误差最小为优化目标函数，以改进粒子群算法为求解该优化目标函数的方法。
优化目标函数为：
其中，n为样本数，
改进粒子群算法具体为：
1a)分别引入非线性函数及余弦函数对粒子群算法的惯性权重及学习因子进行改进；
2a)设初始化种群个数为n，最大迭代次数为n，随机生成n个初始粒子x＝(x1,x2,lxn)，其中，第i个粒子的位置向量及速度向量分别为xi＝(xi1,xi2,lxid)t及vi＝(vi1,vi2,lvid)t，其历史迭代的最优位置pi＝(pi1,pi2,lpid)t，全体粒子的迭代最优位置pg＝(pg1,pg2,lpgd)t，d为单个样本的维度；
3a)以grnn的均方根误差为优化目标，计算各粒子的初始适应度值，再将计算得到的所有粒子的初始适应度进行对比，得到并记录初始适应度值最优的粒子的位置向量及速度向量；
4a)通过更新公式更新粒子的速度及位置；
5a)计算当前各粒子的适应度值，再将当前粒子的适应度值与上一次迭代得到的粒子的适应度值进行比较，并按照择优保留公式进行择优保留，然后更新各粒子的位置及速度；
6a)判断迭代次数是否大于等于预设值、迭代精度是否达到预设目标，当迭代次数大于等于预设值时或者当迭代精度达到预设目标时，则输出最优结果，否则，则转至步骤4a)。
步骤5a)中的择优保留公式为：
其中，
步骤4a)中的更新公式为：
其中，d＝1,2,l,d，i＝1,2ln，k为当前迭代次数，r1,r2在[0，1]区间随机取值。
步骤1a)中，分别引入非线性函数及余弦函数对粒子群算法的惯性权重及学习因子进行改进，其中，
其中，
本发明具有以下有益效果：
本发明所述的基于rs-pso-grnn的埋地油气管道土壤腐蚀预测方法在具体操作时，先识别管道土壤腐蚀因素，再构建埋地管道土壤腐蚀指标体系，以提高本发明的适应性，另外，本通过采用rs理论中的属性约简删除决策表中的冗余指标，以减少预测的误差，最后构建rs-pso-grnn土壤腐蚀预测模型，再将步骤5)得到的训练集输入至rs-pso-grnn土壤腐蚀预测模型中，得各监测点的管道年均腐蚀壁厚，操作简单，预测的精度较高，经验证，与常规的bp土壤腐蚀预测模型及rs-grnn土壤腐蚀预测模型相比，本发明的准确率及预测精度均优于其它两种常规模型，表明本发明具有较强的工程适用性及良好的预测性能。
进一步，本发明分别引入非线性函数及余弦函数对粒子群算法的惯性权重及学习因子进行改进，以提高粒子的迭代寻优能力，有效避免粒子后期容易早熟而陷入局部最优等问题，迭代性能较优。
附图说明
图1是本发明的架构图；
图2为本发明的grnn结构图；
图3为粒子位置更新图；
图4为管壁腐蚀图；
图5为管壁表层脱落图；
图6为本法明的埋地管道土壤腐蚀指标体系图；
图7为仿真试验中3种模型相对误差的对比图；
图8为仿真试验中3种模型预测效果的对比图。
具体实施方式
下面结合附图对本发明做进一步详细描述：
参考图1，本发明所述的基于rs-pso-grnn的埋地油气管道土壤腐蚀预测方法包括以下步骤：
1)识别管道土壤腐蚀因素，再根据管道土壤腐蚀因素构建埋地管道土壤腐蚀指标体系，然后获取埋地管道土壤腐蚀的原始样本指标，并构建原始样本集{xij|i＝1,2,l,m,j＝1,2,l,n}，其中，xij为第i个管道样本的第j个指标值，m及n分别为监测点管道序号总数及管道土壤腐蚀指标的维数；
2)对步骤1)得到的原始样本指标值进行离散化处理，再根据离散化后的条件属性集及决策属性集建立决策表；
3)采用rs理论中的属性约简删除决策表中的冗余指标，得管道核心指标样本集；
4)对步骤3)得到的管道核心指标样本集进行归一化处理；
5)将步骤4)得到的归一化处理后的管道核心指标样本集分为训练样本集及测试样本集，再将训练样本作为grnn输入的训练集；
6)构建rs-pso-grnn土壤腐蚀预测模型，再将步骤5)得到的训练集输入至rs-pso-grnn土壤腐蚀预测模型中，得各监测点的管道年均腐蚀壁厚。
步骤4)中归一化处理的结果为：
其中，max(xj)及min(xj)分别为第j个指标的最大值及最小值，
步骤6)中的rs-pso-grnn土壤腐蚀预测模型以grnn输出的均方根误差最小为优化目标函数，以改进粒子群算法为求解该优化目标函数的方法。
优化目标函数为：
其中，n为样本数，
改进粒子群算法具体为：
1a)分别引入非线性函数及余弦函数对粒子群算法的惯性权重及学习因子进行改进；
2a)设初始化种群个数为n，最大迭代次数为n，随机生成n个初始粒子x＝(x1,x2,lxn)，其中，第i个粒子的位置向量及速度向量分别为xi＝(xi1,xi2,lxid)t及vi＝(vi1,vi2,lvid)t，其历史迭代的最优位置pi＝(pi1,pi2,lpid)t，全体粒子的迭代最优位置pg＝(pg1,pg2,lpgd)t，d为单个样本的维度；
3a)以grnn的均方根误差为优化目标，计算各粒子的初始适应度值，再将计算得到的所有粒子的初始适应度进行对比，得到并记录初始适应度值最优的粒子的位置向量及速度向量；
4a)通过更新公式更新粒子的速度及位置；
5a)计算当前各粒子的适应度值，再将当前粒子的适应度值与上一次迭代得到的粒子的适应度值进行比较，并按照择优保留公式进行择优保留，然后更新各粒子的位置及速度；
6a)判断迭代次数是否大于等于预设值、迭代精度是否达到预设目标，当迭代次数大于等于预设值时或者当迭代精度达到预设目标时，则输出最优结果，否则，则转至步骤4a)。
步骤5a)中的择优保留公式为：
其中，
步骤4a)中的更新公式为：
其中，d＝1,2,l,d，i＝1,2ln，k为当前迭代次数，r1,r2在[0，1]区间随机取值。
步骤1a)中，分别引入非线性函数及余弦函数对粒子群算法的惯性权重及学习因子进行改进，其中，
其中，
仿真试验
所选管道为中俄原油管道，中国境内全长933.11km，2011年投入使用，大部分管段埋入地下，且沿途气候恶劣，地质复杂，多处已发生较为严重的腐蚀，如附图4及图5所示，本发明选取该管道典型监测点中的35组腐蚀数据进行实证分析及研究，将管道年均腐蚀壁厚作为预测结果，其部分数据如下表1所示。
表1
设grnn非线性回归公式主方程为：
其中，估计值
由附图7、图8、表2及表3可知，本发明与bp模型相比，经rs理论处理的样本数据摒弃了土壤腐蚀因素间的冗杂性，rs-pso-grnn模型的均方根误差、平均相对误差及希尔系数分别低于bp模型的6.49％、9.87％及2.25％，说明rs理论能充分挖掘高维数据的耦合信息及冗余规律，提高了模型的预测精度。本发明与rs-svm模型相比，其三种指标值分别降低了3.31％、4.49％及0.12％，表明grnn比svm更加适合埋地管道土壤腐蚀的建模研究。
表2
表3
本发明说明书中未作详细描述的内容属于本领域专业技术人员周知的现有公开技术。
以上实施方式仅用于说明本发明，而并非对本发明的限制。尽管为说明目的公开了本发明的相关实施例和附图，但是本领域的技术人员可以理解；在不脱离本发明及所附的权利要求的精神和范围内，各种替换、变化、修改都是可能的。因此，所有等同的技术方案也属于本发明的范畴，本发明的专利保护范围应由权利要求限定，而不应局限于最佳实施例和附图所公开的内容。