CN113468757B - 一种基于模糊随机理论评估腐蚀天然气管道可靠性的方法 - Google Patents
一种基于模糊随机理论评估腐蚀天然气管道可靠性的方法 Download PDFInfo
- Publication number
- CN113468757B CN113468757B CN202110805563.8A CN202110805563A CN113468757B CN 113468757 B CN113468757 B CN 113468757B CN 202110805563 A CN202110805563 A CN 202110805563A CN 113468757 B CN113468757 B CN 113468757B
- Authority
- CN
- China
- Prior art keywords
- fuzzy
- variable
- pipeline
- establishing
- random
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 title claims abstract description 66
- 238000000034 method Methods 0.000 title claims abstract description 56
- 239000003345 natural gas Substances 0.000 title claims abstract description 33
- 238000005260 corrosion Methods 0.000 claims abstract description 48
- 230000007797 corrosion Effects 0.000 claims abstract description 47
- 230000008569 process Effects 0.000 claims abstract description 20
- 238000006731 degradation reaction Methods 0.000 claims abstract description 18
- 230000015556 catabolic process Effects 0.000 claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims abstract description 8
- 238000006243 chemical reaction Methods 0.000 claims abstract description 8
- 238000005315 distribution function Methods 0.000 claims abstract description 8
- 238000007476 Maximum Likelihood Methods 0.000 claims abstract description 4
- 238000012545 processing Methods 0.000 claims description 7
- 238000010606 normalization Methods 0.000 claims description 6
- 238000000342 Monte Carlo simulation Methods 0.000 description 7
- 238000004458 analytical method Methods 0.000 description 5
- 230000006378 damage Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 230000003628 erosive effect Effects 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000005653 Brownian motion process Effects 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005530 etching Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 231100000572 poisoning Toxicity 0.000 description 1
- 230000000607 poisoning effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 230000036962 time dependent Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
一种基于模糊随机理论评估腐蚀天然气管道可靠性的方法,包括以下步骤:收集影响天然气管道腐蚀的参数,并将这些参数分为退化变量、模糊变量和随机变量;根据Modified B31G准则建立管道腐蚀的模糊随机极限状态方程,建立基于Gamma过程的退化变量的退化轨道模型,利用极大似然函数估计Gamma过程的参数,利用当量正态转化法将退化轨道模型进行正态化处理;建立模糊变量服从模糊正态分布并建立相应的隶属函数,利用当量概率密度法将模糊变量转化为随机变量;建立随机变量的正态分布函数;将处理后的变量带入模糊随机极限状态方程,计算得出腐蚀天然气管道的可靠性。本发明的方法精度高、计算速度快,能够较好的应用于实际工程中。
Description
技术领域
本发明涉及天然气输送技术领域,具体涉及一种基于模糊随机理论评估腐蚀天然气管道可靠性的方法。
背景技术
天然气作为最重要的清洁能源之一,确保其安全可靠的运输是一项重要的任务。作为五大运输方式之一的管道运输,已经成为天然气运输的首选方案,对国民经济的发展、社会安全的保障以及国防建设的增强都有着极其重要的作用。
管道失效的原因不外乎第三方破坏、施工或材料缺陷、介质腐蚀、土体移动等原因,在这众多的失效原因中,介质腐蚀是其中不容忽视的一个重要问题。根据工程经验,在天然气管道投产运营的初期阶段,管道可以完好无损的正常工作,但伴随着时间的流逝,管道内壁的腐蚀缺陷开始逐步加深,当腐蚀深度超过安全界限时,管道就会有泄漏失效的危险。在天然气的输送过程中,管道穿越的地形非常复杂,不仅管道外部的土壤环境多变、地质活动频繁,而且管道内部运输的介质也具有腐蚀性,这些因素都会危害管道的安全运行,而天然气管道一旦发生介质泄露,就容易污染环境、导致沿线居民中毒,甚至引发火灾、爆炸等危险,从而造成严重的财产损失甚至危及生命安全。
为了确保天然气管道安全可靠的运行,工作人员需要对管道其进行定期的检测,分析管道内部产生缺陷的原因,并由研究人员对管道进行可靠性分析,根据分析结果制定对应的维修维护措施,从而使得天然气管道能在服役期间安全运行。
天然气管道一旦失效会产生严重的后果,而由腐蚀引起的失效是天然气管道失效的重要原因。由于管道的腐蚀问题具有非常强的复杂性和不确定性,所以工作人员有必要细致分析管道腐蚀的产生与发展的内部机理,加强对天然气管道的腐蚀速率和可靠性的研究与探讨,并根据分析结果采取对应的防护措施,才能有效地避免管道失效。
对于高寿命产品而言,很多学者都采用了退化轨道对其进行可靠性分析。该理论的优点是可以充分利用产品在使用过程中的参数变化信息,提高了信息利用率。但为了简化模型,很多学者都将产品的退化轨道选定为维纳过程,但并非所有的退化过程在有限时间上的变化都服从正态分布,所以该模型有一定的局限性。因此,有必要针对影响天然气管道腐蚀的各种因素进行分类处理,建立不同的影响模型,使得最终的结果更加准确,为管道的安全性提供更加准确的数据支持。
发明内容
为解决至少一个上述问题,本发明的目的在于提供一种基于模糊随机理论评估腐蚀天然气管道可靠性的方法,其精度高、计算速度快。
本发明提供的技术方案是,一种基于模糊随机理论评估腐蚀天然气管道可靠性的方法,包括以下步骤:
步骤1、收集影响天然气管道腐蚀的参数,并将这些参数分为退化变量、模糊变量和随机变量,所述退化变量为内壁腐蚀深度,所述模糊变量为模糊安全系数,所述随机变量为管道屈服强度、管道外径、管道壁厚、腐蚀长度和管道运行压力;
步骤2、根据Modified B31G准则建立管道腐蚀的模糊随机极限状态方程,
步骤3、建立基于Gamma过程的退化变量的退化轨道模型,利用极大似然函数估计Gamma过程的参数,利用当量正态转化法将退化轨道模型进行正态化处理;
步骤4、建立模糊变量服从模糊正态分布并建立相应的隶属函数,利用当量概率密度法将模糊变量转化为随机变量;
步骤5、建立随机变量的正态分布函数;
步骤6、将步骤3-5中的变量带入步骤2中的模糊随机极限状态方程,计算得出腐蚀天然气管道的可靠性。
所述极限状态方程为
当J≥50时,M=0.032J+0.33。
本发明的技术效果是:
本发明将影响腐蚀天然气管道可靠性的各变量,按照其各自特性分为退化变量、随机变量和模糊变量,后将退化变量和模糊变量正态化,最终建立模糊随机极限状态方程并计算得出管道的可靠性,采用本发明的方法,其精度高、计算效率高。
附图说明
图1为本发明的流程框图。
具体实施方式
下面结合实施例及附图,对本发明作进一步地的详细说明。
为使本发明实施方式的目的、技术方案和优点更加清楚,下面将结合本发明实施方式中的附图,对本发明实施方式中的技术方案进行清楚、完整地描述,显然,所描述的实施方式是本发明一部分实施方式,而不是全部的实施方式。基于本发明中的实施方式,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施方式,都属于本发明保护的范围。因此,以下对在附图中提供的本发明的实施方式的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施方式。
实施例:
步骤1、收集影响天然气管道腐蚀的参数,并将这些参数分为退化变量、模糊变量和随机变量,所述退化变量为内壁腐蚀深度,所述模糊变量为模糊安全系数,所述随机变量为管道屈服强度、管道外径、管道壁厚、腐蚀长度和管道运行压力;
对于上述变量,内壁腐蚀深度需要收集持续一段时间内的数据及其变化情况。
步骤2、根据Modified B31G准则建立管道腐蚀的极限状态方程,并将退化变量、模糊变量以及随机变量分别处理后的结果带入极限状态方程中,计算得出腐蚀天然气管道的可靠性。
具体的,模糊极限状态方程为如式13所示:
当J≥50时,M=0.032J+0.33。
步骤3、建立基于Gamma过程的退化变量的退化轨道模型,利用极大似然函数估计Gamma过程的参数,利用当量正态转化法将退化轨道模型进行正态化处理;
其具体过程为:
1、首先,建立内壁腐蚀深度d基于Gamma过程的退化轨道模型,如式2所示,
式中,Γ(·)为Gamma函数;
T为时间;
η(T)为Gamma过程中的形状参数,与时间有关;
λ为Gamma过程中的尺度参数,与时间无关;
d(T)为腐蚀深度,与时间有关。
2、在上述的退化轨道模型中,首先需要根据前期收集到的内壁腐蚀深度d的数据确定Gamma过程的待定参数,待定参数通常包括η(t)和λ,其中,一般假定η(t)=ct(c>0),其中,c为常数,此时,对Gamma过程中的参数η(t)和λ的估计变成了对c、λ的估计。
设Ti时刻的腐蚀深度为di,且i=(1,…,n),0=T0<T1<T2<…<Tn,di(i=1,2,…,n)。
则退化增量δi=di-di-1(i=1,2,…,n),参数c和λ的似然函数如式3所示
式中,Ti为测量时间,s;d为腐蚀深度,mm;δi为退化增量,mm;
对上述的似然函数中的参数c、λ求导并化简后,可得到如式4所示的下述方程组
采用当量正态转化法对内壁腐蚀深度d进行当量正态化处理,得到内壁腐蚀深度d的当量均值和当量方差,具体步骤如下所示:
式中,Φ-1(·)为标准正态分布函数的反函数;
对于Gamma分布d~Ga(η(T),λ)而言,其均值的计算公式如式7所示,方差的计算公式如式8所示
式中,μ(d(T))为腐蚀深度d的均值;
将式7和式8分别带入式5和式6,同时采用腐蚀深度d替换变量Xi,即可将非正态变量腐蚀深度d转化为等效正态变量。
步骤4、建立模糊变量服从模糊正态分布并建立相应的隶属函数,利用当量概率密度法将模糊变量转化为随机变量;
式中x0表示该模糊变量隶属函数的中间值;
a表示安全系数,是该隶属函数的自变量;
α,β表示模糊正态分布中的两个参数;
G1为左参照函数,表示在中间值左侧时模糊正态隶属函数的表达式;
G2为右参照函数,表示在中间值右侧时模糊正态隶属函数的表达式。
Φ(·)为标准正态分布的分布函数。
则当量模糊安全系数aT的当量均值如式11所示:
当量标准差的计算公式如式12所示:
步骤5、建立随机变量的正态分布函数;
具体的,随机变量包括管道屈服强度σy、管道外径D、管道壁厚t、腐蚀长度L和管道运行压力Pop,通过求解其正态分布函数,求得这些随机变量的均值和标准差。
步骤6、将步骤3中的内壁腐蚀深度d当量正态化值、步骤4中模糊安全系数的当量模糊安全系数以及步骤5中其余随机变量的值带入式1中,利用基于验算点法的一次二阶矩理论,并结合matlab进行编程处理,即可完成对天然气管道的模糊随机可行性分析。
下面为了进一步说明本发明的方法,采用现场生产的例子进行说明。
某天然气输送管线,2005年建成投产,其管段基本信息如表1所示
表1管线基本信息
项目 | 参数 | 项目 | 参数 |
管道级别 | GC1 | 传输介质 | 天然气 |
防腐层材料 | 三层PE | 管道材质 | L360钢 |
内壁厚度 | 8mm | 设计压力 | 4MPa |
同时,连续测量其腐蚀深度的数据如表2所示
表2不同测量点腐蚀深度测量数据(mm)
在管道的腐蚀失效问题中,导致管道可靠性急剧下降的主要原因是管道内壁腐蚀最严重的点产生了腐蚀穿孔现象。通过统计观察这十个检测点的腐蚀深度随时间的变化规律可以发现,腐蚀程度最严重的点为第九个检测点,所以将该点的腐蚀数据作为本文中天然气管道退化轨道的退化数据。
表3极限状态方程中的参数分布
按照上述方法进行计算,得到其可靠度计算结果,如表4所示,
表4基于模糊随机极限状态方程的可靠度计算结果
为了验证本方法的准确性,采用本领域用于验证准确性的蒙特卡洛法进行进一步的验证,最终结果如表4所示,
表4基于蒙特卡洛法的可靠度验证结果
以2025年12月的可靠度值为例,给出本发明的方法、蒙特卡洛法以及常规一次二阶矩法之间的对比,具体如表5所示。
表5模糊随机理论和蒙特卡洛法比较示例
方法 | 可靠度 | 运算次数 | 相对误差 |
模糊随机理论 | 0.997102 | 52 | 0.19% |
蒙特卡洛数值模拟 | 0.995191 | 100000 | — |
常规一次二阶矩法 | 0.627718 | 33 | 36.9% |
从表3-表5可以看出:蒙特卡洛法的虽然能计算出更为精确的可靠度,但是其运算次数远超其它方法,所以处理复杂的工程问题时耗时太长,只能作为检验其它方法精确度的理论依据;常规一次二阶矩法中,将所有变量设为符合正态分布的随机变量,可看出其误差较大,精度较低;基于模糊随机理论计算出的可靠度,与蒙特卡洛数值模拟的结果吻合效果较好,同时其运算次数远远小于蒙特卡洛法,具有良好的工程实际意义。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明实施例揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (4)
1.一种基于模糊随机理论评估腐蚀天然气管道可靠性的方法,其特征在于,包括以下步骤:
步骤1、收集影响天然气管道腐蚀的参数,并将这些参数分为退化变量、模糊变量和随机变量,所述退化变量为内壁腐蚀深度,所述模糊变量为模糊安全系数,所述随机变量为管道屈服强度、管道外径、管道壁厚、腐蚀长度和管道运行压力;
步骤2、根据Modified B31G准则建立管道腐蚀的模糊随机极限状态方程;
步骤3、建立基于Gamma过程的退化变量的退化轨道模型,利用极大似然函数估计Gamma过程的参数,利用当量正态转化法将退化轨道模型进行正态化处理;
步骤4、建立模糊变量服从模糊正态分布并建立相应的隶属函数,利用当量概率密度法将模糊变量转化为随机变量;
步骤5、建立随机变量的正态分布函数;
步骤6、将步骤3-4计算得出的变量以及步骤5中的变量带入步骤2中的模糊随机极限状态方程,通过一次二阶矩法,计算得出腐蚀天然气管道的可靠性。
3.根据权利要求1所述的方法,其特征在于,所述步骤3的具体操作为:
S31、建立内壁腐蚀深度基于Gamma过程的退化轨道模型;
S32、建立Gamma过程的似然函数,代入内壁腐蚀深度随时间的测量值,估算得出Gamma过程的形状参数和尺度参数;
S33、利用当量正态转化法将退化轨道模型进行当量正态化处理,求出其当量均值和当量方差,同时设定内壁腐蚀深度阈值,可计算。
4.根据权利要求1所述的方法,其特征在于,步骤4中,其具体步骤为:
S41、建立模糊安全系数的模糊正态隶属函数;
S42、采用当量概率密度法对模糊安全系数进行处理,得到模糊安全系数的当量概率密度函数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110805563.8A CN113468757B (zh) | 2021-07-16 | 2021-07-16 | 一种基于模糊随机理论评估腐蚀天然气管道可靠性的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110805563.8A CN113468757B (zh) | 2021-07-16 | 2021-07-16 | 一种基于模糊随机理论评估腐蚀天然气管道可靠性的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113468757A CN113468757A (zh) | 2021-10-01 |
CN113468757B true CN113468757B (zh) | 2022-08-30 |
Family
ID=77880725
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110805563.8A Expired - Fee Related CN113468757B (zh) | 2021-07-16 | 2021-07-16 | 一种基于模糊随机理论评估腐蚀天然气管道可靠性的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113468757B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114046456B (zh) * | 2021-11-23 | 2024-02-02 | 重庆大学 | 融合模糊推理和神经网络的腐蚀评估方法及系统 |
CN115308113B (zh) * | 2022-07-01 | 2024-07-19 | 西南石油大学 | 基于多元非线性拟合的页岩气集输管道腐蚀速率预测方法 |
CN116384817B (zh) * | 2023-03-30 | 2024-01-09 | 沧州市农林科学院 | 一种用于评价植物抗逆性的方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108460230B (zh) * | 2018-03-21 | 2022-07-08 | 西安建筑科技大学 | 基于数据融合的腐蚀油气管道剩余寿命及可靠性预测方法 |
CN109299544A (zh) * | 2018-09-26 | 2019-02-01 | 东北大学 | 一种基于gamma过程的机械零件二元相关退化可靠性评估方法 |
CN110929453A (zh) * | 2019-11-18 | 2020-03-27 | 西安电子科技大学 | 基于Copula函数失效相关系统动态模糊可靠性分析方法 |
CN112949190A (zh) * | 2021-03-08 | 2021-06-11 | 西南石油大学 | 基于r6-fad和bp-mcs的含裂纹缺陷管道的模糊随机可靠度评估方法 |
CN112949209B (zh) * | 2021-03-26 | 2022-05-17 | 北京航空航天大学 | 一种退化速率-波动联合更新的弹用密封橡胶贮存寿命评估方法 |
-
2021
- 2021-07-16 CN CN202110805563.8A patent/CN113468757B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN113468757A (zh) | 2021-10-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113468757B (zh) | 一种基于模糊随机理论评估腐蚀天然气管道可靠性的方法 | |
CN103455682B (zh) | 一种预测高温高压井腐蚀套管剩余寿命的方法 | |
CN104061445B (zh) | 一种基于神经网络的管道泄漏检测方法 | |
CN116629712B (zh) | 基于pso-bp神经网络的海底盾构隧道施工质量风险评估方法 | |
CN106355320B (zh) | 一种灰色系统海底油气管道坠物腐蚀可靠性评估方法 | |
Shojai et al. | Probabilistic modelling of pitting corrosion and its impact on stress concentrations in steel structures in the offshore wind energy | |
CN104063588A (zh) | 基于多源数据融合的管道腐蚀缺陷尺寸的预测系统及方法 | |
Okoro et al. | Multi-criteria risk assessment approach for components risk ranking–The case study of an offshore Wave Energy Converter | |
CN108920421A (zh) | 一种点腐蚀作用下油气管道系统的随机失效概率评估方法 | |
CN115358102B (zh) | 引入内检测数据的地震作用下腐蚀管道失效概率分析方法 | |
CN115186590A (zh) | 一种管道腐蚀剩余寿命预测方法 | |
Thodi et al. | The development of posterior probability models in risk‐based integrity modeling | |
Xu et al. | An improved regularized particle filter for remaining useful life prediction in nuclear plant electric gate valves | |
CN113705609A (zh) | 一种油气管道风险判定模型的构建方法和装置 | |
CN109977563B (zh) | 基于tv-hsmm的海底油气管道剩余寿命预测方法 | |
CN114021078A (zh) | 一种大坝监测量最佳统计模型优选方法 | |
CN111341396A (zh) | 一种大气环境对材料腐蚀安全评估方法及系统 | |
CN112241843A (zh) | 一种带腐蚀缺陷油套管剩余强度评价方法 | |
CN105224707A (zh) | 一种盐穴地下储气库溶腔体积收缩率的预测方法 | |
Wang et al. | Multi-parameter maximum corrosion depth prediction model for buried pipelines based on GSCV-XGBoost | |
AU2021101888A4 (en) | System and method for corrosion prediction in oil and gas pipeline | |
CN114201840B (zh) | 管段应力腐蚀开裂风险识别方法、装置、设备及介质 | |
Syzrantsev et al. | Improvement of calculating methods for trunk pipelines reliability | |
Yusof et al. | Markov chain model for predicting pitting corrosion damage in offshore pipeline | |
CN109784590B (zh) | 一种基于cagm(1,1)-bpnn的在役油气管道腐蚀预测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220830 |
|
CF01 | Termination of patent right due to non-payment of annual fee |