[发明专利]一种管道腐蚀人工智能预警方法有效
申请号: | 202111178474.1 | 申请日: | 2021-10-10 |
公开(公告)号: | CN113806902B | 公开(公告)日: | 2022-08-19 |
发明(设计)人: | 曾德智;韩雪;金龙;张新;于晓雨;赵春兰;仝春玥;汪宙峰;董宝军;喻智明 | 申请(专利权)人: | 西南石油大学 |
主分类号: | G06F30/18 | 分类号: | G06F30/18;G06F30/27;G06N3/04;G06N3/08;G06F113/14;G06F119/04 |
代理公司: | 暂无信息 | 代理人: | 暂无信息 |
地址: | 610500 四*** | 国省代码: | 四川;51 |
权利要求书: | 查看更多 | 说明书: | 查看更多 |
摘要: | |||
搜索关键词: | 一种 管道 腐蚀 人工智能 预警 方法 | ||
1.一种管道腐蚀人工智能预警方法,其特征在于,包括如下步骤:
步骤1:收集管道的基础数据:
管道外径Dw,mm;(2)管道内径Dn,mm;(3)管道材料的屈服强度σs,MPa;(4)管道壁厚d,mm;(5)管道设计压力P,MPa;(6)管道腐蚀裕量C,mm;(7)管道生产运行时间Ts,a;(8)管道设计使用年限Tu;
步骤2:划分管道腐蚀缺陷区域:
检测管道腐蚀缺陷按照轴向和环向对区域进行网格划分:轴向划分m份分别为C1、C2…Ci…Cm,环向划分n份分别为L1、L2…Lj…Ln,从而将腐蚀缺陷离散划分为m×n个壁厚测量点Aij(i=1、2、3…m;j=1、2、3…n);
其中:m为轴向划定的区域个数,C1、C2…Cm为轴向划定区域的每一份对应一个缺陷的轴向测量点;n为环向划定的区域个数,L1、L2…Ln为轴向划定区域的每一份对应一个缺陷的环向测量点;Aij为腐蚀缺陷离散划分的m×n个壁厚测量点;
步骤3:管道继续服役条件的判定,具体步骤:
(a)通过公式(1)求解轴向要求最小壁厚通过公式(2)求解环向要求最小壁厚将计算结果和带入公式(3)确定管道最小要求壁厚tmin:
式中:为轴向要求最小壁厚,mm;为环向要求最小壁厚,mm;tmin为管道最小要求壁厚,mm;
(b)统计管道腐蚀缺陷区域的m×n个壁厚测量点的壁厚值aij,其中测量得到壁厚值最小的是amin,通过公式(4)求解所有被测点壁厚的平均值tam:
式中:aij为腐蚀缺陷离散划分的m×n个测量点的壁厚值,mm;tam为所有被测点壁厚的平均值,mm;
(c)通过公式(5)求解管道的剩余壁厚比Rt:
式中:Rt为管道的剩余壁厚比;
(d)求解管道轴向最大允许腐蚀缺陷的长度L:
若Rt≥0.793,则
若Rt<0.793,则
其中:L为轴向最大允许腐蚀缺陷的长度值,mm;
(e)确定管道安全服役条件:
若La≤L,则管道可以继续服役;
若LaL且tam-C≥0.9tmin时,则管道可以继续服役;
若LaL且tam-C<0.9tmin时,则管道不可以继续服役,判定为一级预警等级,预警标识显示为红色;
其中:La为管道壁厚断面腐蚀缺陷的轴向长度,mm;
步骤4:预测管道剩余寿命:
(a)基于在线监测预测管道剩余寿命:
1)预测管道腐蚀速率,具体步骤:
①整理现场腐蚀监测数据集:某一时间段的监测时间数据集α(t1、t2、t3…tn);监测时间对应的腐蚀速率数据集β(V1、V2、V3…Vn);
其中:t1、t2、t3…tn为按照时间排序间隔一个步长的监测时间;V1、V2、V3…Vn为监测时间对应的腐蚀速率监测值,mm/a;
②观察腐蚀监测数据集补充缺失数据:
任取腐蚀监测时间数据集α中的3个时间ta、tb、tc,取腐蚀速率数据集β中对应的腐蚀速率监测值Va、Vb、Vc,带入公式(6)得到缺失时间tx对应的腐蚀速率值Vx:
式中:ta、tb、tc为监测时间数据集α中的任意三个时间;Va、Vb、Vc为腐蚀速率数据集β中ta、tb、tc对应的腐蚀速率监测值,mm/a;tx为缺失的监测时间;Vx为tx对应的腐蚀速率,mm/a;
③构建长短时记忆神经网络:
I完整的腐蚀监测数据集:包含腐蚀数据缺失值的监测时间数据集α’(t1、t2、t3…tx…tn)和对应的腐蚀速率数据集β’(V1、V2、V3…Vx…Vn),其中α’为输入值,β’为输出值;
II长短时记忆神经网络模型内部结构包含遗忘门、输入门、输出门:
i上一时间ti-1的腐蚀速率监测Vi-1经过ti时间长短时神经网络的遗忘门f(ti)式(7)进行更新遗忘:
式中:f(ti)表示为遗忘门函数;σ为神经网络激活函数,此处选用sigmoid函数;wf与uf为遗忘门权重系数矩阵;bf为网络偏置值;ti为目标预测时间;表示ti时间对应的时间数据;ti-1为监测时间数据集α’的某一时间;Vi-1表示对应的腐蚀速率数据集β’中的腐蚀速率,mm/a;
iiti对应的时间数据ti-1时间对应的腐蚀监测值Vi-1进入ti时间长短时神经网络的输入门i(ti)式(8)与备选内容式(9):
式中:i(ti)表示输入门的决定系数;σ为神经网络激活函数,此处选用sigmoid函数;wi与ui为输入门决定系数权重系数矩阵;bi表示输入门决定系数矩阵的偏置值;表示输入门备选内容;wc表示输入门备选内容权重矩阵;bc表示输入门备选内容的偏置值;tanh表示双曲正切激励函数;
iii利用ti时刻的遗忘门f(ti)与输入门决定系数i(ti)、备选内容进行神经元状态当前的更新,得到ti时刻的神经元更新函数Cti式(10):
式中:表示ti-1时间的神经更新函数;表示ti时间的神经更新函数;
ivti时间的神经元更新函数ti对应的时间数据和上一时间ti-1时刻对应的腐蚀速率监测值Vi-1经过ti时间的输出门式(11),输出腐蚀速率预测值Vi式(12):
式中:表示输出门决定系数;σ表示神经网络激活函数,此处选用sigmoid函数;wo表示输出门权重矩阵;bo表示输出门的偏置值;Vi表示目标时间ti的腐蚀速率预测值,mm/a;
2)将上述预测得到的腐蚀速率预测值Vi、腐蚀管道最小要求壁厚tmin和所有被测点壁厚的平均值tam带入公式(13)求解管道剩余寿命T'L:
式中:T'L为基于在线监测管道剩余寿命,a;
K为安全系数,当La≤L时,K取值为1;当La>L时,K取值为0.9;
(b)基于超声波测厚预测管道剩余寿命:
1)估算管道腐蚀速率Vμ式(14):
式中:Vμ为腐蚀速率估值,mm/a;△d为同点测厚前后壁厚之差,mm;△T为测厚前后时间差,a;
2)求解管道计算壁厚ε式(15):
式中:ε为钢管计算壁厚,mm;P为设计压力,MPa;σs为钢管屈服强度,MPa;
3)将腐蚀速率估值Vμ和管道计算壁厚ε带入计算公式(16)求解腐蚀剩余寿命T″L:
式中:T″L为基于超声波测厚管道剩余寿命,a;
(b)确定管道剩余寿命TL:
TL=min(T′L,T″L) (17)
式中:TL为管道剩余寿命,a;
步骤5:管道安全服役条件的判定:
(a)若TL≥Tu-Ts,则管道可以继续安全服役,判定为五级预警等级,预警标识显示为绿色;
(b)若TL<Tu-Ts,则需进一步判定腐蚀预警等级;
其中:Ts为管道生产运行时间,a;Tu为设计年限,a;
步骤6:腐蚀预警等级设定条件:
(a)若TsTLTu-Ts时,则判定为四级预警等级,预警标识显示为蓝色;
(b)若10TLTu时,则判定为三级预警等级,预警标识显示为黄色;
(c)若3TL10时,则判定为二级预警等级,预警标识显示为橙色;
(d)若TL3时,则判定为一级预警等级,预警标识显示为红色;
步骤7:预警处理:
(a)对于一级预警,腐蚀程度非常严重,需更换腐蚀管道;
(b)对于二级预警,腐蚀程度严重,管道需停止运行并进行修理;
(c)对于三级预警,腐蚀程度比较重,管道降压运行并进行修理;
(d)对于四级预警,注意观察管道腐蚀情况;
(e)对于五级预警,表示管道腐蚀控制情况较好,管道可以安全生产运行。
该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于西南石油大学,未经西南石油大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服】
本文链接:http://www.vipzhuanli.com/pat/books/202111178474.1/1.html,转载请声明来源钻瓜专利网。