CN113626970B - 公共管廊管道腐蚀剩余寿命评估方法及*** - Google Patents
公共管廊管道腐蚀剩余寿命评估方法及*** Download PDFInfo
- Publication number
- CN113626970B CN113626970B CN202010383246.7A CN202010383246A CN113626970B CN 113626970 B CN113626970 B CN 113626970B CN 202010383246 A CN202010383246 A CN 202010383246A CN 113626970 B CN113626970 B CN 113626970B
- Authority
- CN
- China
- Prior art keywords
- corrosion
- pipeline
- parameter
- distribution model
- corrosion depth
- 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.)
- Active
Links
- 230000007797 corrosion Effects 0.000 title claims abstract description 209
- 238000005260 corrosion Methods 0.000 title claims abstract description 209
- 238000000034 method Methods 0.000 title claims abstract description 43
- 230000001186 cumulative effect Effects 0.000 claims abstract description 23
- 238000012163 sequencing technique Methods 0.000 claims abstract description 10
- 238000010586 diagram Methods 0.000 claims description 9
- 238000012360 testing method Methods 0.000 claims description 9
- 230000014509 gene expression Effects 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 7
- 229910000831 Steel Inorganic materials 0.000 claims description 5
- 239000010959 steel Substances 0.000 claims description 5
- 238000013507 mapping Methods 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 abstract description 14
- 238000011156 evaluation Methods 0.000 description 9
- 238000013461 design Methods 0.000 description 6
- 238000001514 detection method Methods 0.000 description 4
- 230000011218 segmentation Effects 0.000 description 4
- 239000000463 material Substances 0.000 description 3
- 238000001276 Kolmogorov–Smirnov test Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 238000010998 test method Methods 0.000 description 2
- PPBRXRYQALVLMV-UHFFFAOYSA-N Styrene Natural products C=CC1=CC=CC=C1 PPBRXRYQALVLMV-UHFFFAOYSA-N 0.000 description 1
- 238000003556 assay Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004880 explosion Methods 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 231100000572 poisoning Toxicity 0.000 description 1
- 230000000607 poisoning effect Effects 0.000 description 1
- 239000010865 sewage Substances 0.000 description 1
- 125000003011 styrenyl group Chemical group [H]\C(*)=C(/[H])C1=C([H])C([H])=C([H])C([H])=C1[H] 0.000 description 1
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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
- G06Q50/265—Personal security, identity or safety
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- 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]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/30—Computing systems specially adapted for manufacturing
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Computer Security & Cryptography (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Educational Administration (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Economics (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
本发明公开了一种公共管廊管道腐蚀剩余寿命评估方法,包括:将被检测管道分成多个管段,对所有管段的腐蚀深度值的最大值进行排序,获取被检测管道的腐蚀深度值的样本序列,以及获取任一腐蚀深度值的累积概率;采用广义极值分布模型并进行参数估计,根据确定的参数选择被检测管道的腐蚀数据分布模型;应用所确定的腐蚀数据分布模型预测不同可靠度下被检测管道的极限腐蚀深度;根据被检测管道在不同可靠度下的极限腐蚀深度确定被检测管道在不同可靠度下的腐蚀剩余寿命。本发明还公开了公共管廊管道腐蚀剩余寿命评估***,包括:数据获取模块、状态分析模块和状态预测模块。本发明能够有效预测管廊管道腐蚀深度及评估管廊管道剩余寿命。
Description
技术领域
本发明涉及管道安全评价技术领域,尤其涉及公共管廊管道腐蚀剩余寿命评估方法及***。
背景技术
管廊管道分布密集,情况复杂,且大部分管道负责运输高危或中危介质,一旦发生管道失效,不仅会影响园区正常运行,还易导致火灾、***或中毒事故,甚至造成多米诺效应,对周围群众的生命财产安全构成严重威胁。其中,造成公共管廊管道发生失效的一个重要原因在于所发生的腐蚀,管廊管道发生腐蚀后将会降低管道的承压能力,可能会引起管道泄漏,最终造成严重的生产事故。因此,如何准确评估管廊管道的腐蚀速率,预测管道的腐蚀深度,评估管道的腐蚀剩余寿命已成为目前管廊管道安全评价中所需要重点解决的问题。
对于管道腐蚀深度的预测主要是通过对检测数据进行统计分析,采用合理的模型对管道材料的腐蚀深度进行预测。已有的研究多为采用某一数据分布模型对不同条件下的材料的腐蚀情况进行研究。统一的数据分布模型无法同时契合对应不同管道的不同检测数据的各类情况,而且考虑到管廊管道的特殊性、复杂性,现有技术的管道腐蚀深度预测方法显然已经不合适。而且目前对于管廊管道腐蚀情况的研究也较为不足,缺乏一种较为完整且可靠的腐蚀深度预测方法,以及管廊管道剩余寿命的评估方法。
发明内容
本发明的目的在于提供一种公共管廊管道腐蚀剩余寿命评估方法及***,有效预测管廊管道腐蚀深度及评估管廊管道剩余寿命。
实现上述目的的技术方案是:
一种公共管廊管道腐蚀剩余寿命评估方法,包括:
腐蚀数据处理及选取:将被检测管道分成多个管段,对所有管段的腐蚀深度值的最大值进行排序,获取被检测管道的腐蚀深度值的样本序列,以及获取任一腐蚀深度值的累积概率;
选择腐蚀数据分布模型:采用广义极值分布模型并进行参数估计,根据确定的参数选择被检测管道的腐蚀数据分布模型为Gumble分布模型、Frechet分布模型或者Weibull分布模型;
预测腐蚀深度:应用所确定的腐蚀数据分布模型预测不同可靠度下被检测管道的极限腐蚀深度;
预测腐蚀剩余寿命:根据被检测管道在不同可靠度下的极限腐蚀深度确定被检测管道在不同可靠度下的腐蚀剩余寿命。
优选的,所述的腐蚀数据处理及选取中,
每个管段的腐蚀深度值按由小到大的顺序进行排序,根据排序结果得到该管段的腐蚀深度值的最大值;
对所有管段的腐蚀深度值的最大值按从小到大的顺序进行排序,形成被检测管道的腐蚀深度值的样本序列:
D={x1,x2,x3,…xi,…,xm};
其中,xi表示被检测管道的第i个腐蚀深度值;m表示m个管段;
对于第i个腐蚀深度值,其累积概率为:
对被检测管道的所得到的所有(xi,F(xi))数据点进行做图,并画出曲线,得到被检测管道腐蚀深度的经验分布图。
优选的,所述的选择腐蚀数据分布模型中,
广义极值分布模型的表达式为:
其中,α为位置参数,β为尺度参数,k为广义极值模型分布参数,x表示随机变量;
k参数的值在预设的0附近的范围值内时,则被检测管道的腐蚀数据分布符合Gumble分布,其累积概率函数为:
其中,α为位置参数,β为尺度参数;
k参数的值大于0且在预设的0附近的范围值外时,则被检测管道的腐蚀数据分布符合Frechet分布,其累积概率函数为:
其中,α为位置参数,β为尺度参数,γ为形状参数;X表示随机变量的最大值;
k参数的值小于0且在预设的0附近的范围值外时,则被检测管道的腐蚀数据分布符合Weibull分布,其累积概率函数为:
其中,α为位置参数,β为尺度参数,γ为形状参数;X表示随机变量的最大值;
参数k通过将广义极值分布模型的表达式所述的曲线进行拟合得到。
优选的,在选择腐蚀数据分布模型之后,对选定的腐蚀数据分布模型适用性通过拟合优度检验进行判断。
优选的,所述的拟合优度检验的方法包括K-S检验方法、作图法、回归法、X2检验以及A-D检验。
优选的,所述的选择腐蚀数据分布模型中,
参数估计采用概率权重矩法:
其中,ωr表示概率权重矩;F表示累积概率;r表示实数;dF表示积分时的变量;
得到位置参数α及尺度参数β的确定方法,两个参数的表达式如下:
其中,Γ表示Gamma函数;b0,b1表示无偏估计量。
优选的,所述的预测腐蚀剩余寿命中,
计算得到管道最小允许剩余厚度,如下所示:
ΔA=A-Amin
其中,Amin为最小允许剩余厚度;p为管道运行压力;D为管道外径;SYMS为管线钢屈服强度;A为管道壁厚;ΔA为管道临界腐蚀深度;
结合不同可靠度下的腐蚀深度值,确定不同可靠度下的腐蚀进展曲线,如下:
x=c·tn
其中,x表示管道腐蚀深度值;c表示系数;t表示管道结构开始破坏的使用时间;n为时间常数,取0.3-0.5;
计算得到管道在不同可靠度下的剩余寿命。
本发明的公共管廊管道腐蚀剩余寿命评估***,包括:数据获取模块、状态分析模块和状态预测模块,
所述数据获取模块获取目标管段不同位置的管壁腐蚀减薄数据,统计整理腐蚀数据,并传输给所述状态分析模块;
所述状态分析模块采用广义极值分布模型针对管道腐蚀状况分析结果,按被检测目标管道的不同管段进行分类整理;
所述状态预测模块根据状态分析模块所得结果对管廊管道的状态进行预测,预测管道腐蚀的发展状况,获得管道在不同可靠度下剩余使用寿命。
本发明的有益效果是:本发明针对管廊管道的特殊性、复杂性,通过腐蚀数据处理和选择腐蚀数据分布模型,从而利用不同管道的不同检测数据的符合模型预测不同可靠度下被检测管道的极限腐蚀深度,完整并可靠地确定被检测管道在不同可靠度下的腐蚀剩余寿命。其中,选取腐蚀深度值的最大值排序得到样本序列,为后续结果的准确性提供保障。利用广义极值分布模型确定被检测管道最合适的腐蚀数据分布模型,提高最终结果的准确率和效率。通过拟合优度检验确定所选择的腐蚀数据分布模型对于所检测管道腐蚀数据的拟合适用性,保证选择最优的腐蚀数据分布模型。从而为处于复杂条件下的管廊管道提供一种能准确评价管道运行状态,预测管道剩余寿命的方法,提高管廊管道的安全运行能力,为化工园区的安全运营提供技术保障。
附图说明
图1是本发明的公共管廊管道腐蚀剩余寿命评估方法的流程图;
图2本发明的公共管廊管道腐蚀剩余寿命评价***;
图3是本发明中管廊管道1的走势图;
图4是本发明中管廊管道2的走势图;
图5是本发明中管廊管道3的走势图;
图6是本发明中管廊管道1腐蚀数据拟合图;
图7是本发明中管廊管道2腐蚀数据拟合图;
图8是本发明中管廊管道3腐蚀数据拟合图。
具体实施方式
下面将结合附图对本发明作进一步说明。
请参阅图1,本发明的公共管廊管道腐蚀剩余寿命评估方法,包括下列步骤:
步骤S01,腐蚀数据处理及选取。
首先对被检测管道的整体长度进行确定,并根据检测管道线路的具体状况,对被检测管道进行分段处理。在分段过程中需具体结合实际情况,合理确定管道的分段数,分段数不宜过多或过少。
之后,被检测管道进行平均分段后,对被检测管道的每一管段检测得到的腐蚀深度进行统计处理。具体地,对于每个管段的腐蚀深度值按由小到大的顺序进行排序,根据排序结果得到该管段的腐蚀深度最大值。通过对数据统计后,得到被检测管道每个管段的腐蚀深度值的最大值。
之后,对被检测管道的所有管段的腐蚀深度值的最大值进行排序。具体地,所采用的排序方法同样为将数据按从小到大的顺序进行排序。通过进行排序,形成得到被检测管道的腐蚀深度值的样本序列:
D={x1,x2,x3,…xi,…,xm} (1)
在样本序列中:
x1≤x2≤x3≤…≤xi≤…≤xm (2)
其中,xi表示被检测管道的第i个腐蚀深度值,即第i个大的管段的腐蚀深度值的最大值;m表示m个管段。
之后,根据排序得到的被检测管道腐蚀深度值的样本序列,可以得到任意测量值的累积概率。具体地,对于第i个腐蚀深度值,其累积概率为:
对被检测管道的每个管段的腐蚀深度值的最大值,应用式(3),计算得到任意腐蚀深度值xi所对应的累积概率F(xi)。
最后,对被检测管道的所得到的所有(xi,F(xi))数据点进行做图,并画出曲线,得到被检测管道腐蚀深度的经验分布图。此处得到的曲线用于下述步骤中的拟合。
步骤S02,腐蚀数据分布模型选择。
对腐蚀数据分布进行分析所采用的为广义极值分布模型,具体地,该模型的表达式如下:
在该模型中,涉及到三个参数,α为位置参数,β为尺度参数,k为广义极值模型分布参数。x表示随机变量。广义极值分布模型综合了Gumble分布模型、Frechet分布模型及Weibull分布模型三种分布模型,并通过参数k确定具体选择的模型。如下:
1)当广义极值分布模型中的k参数的值接近于0(在预设的0附近的范围值内)时,则被检测管道的腐蚀数据分布符合Gumble分布。Gumble分布,即第一渐近分布的最大分布或双指数分布,分析局部腐蚀的最大值。其累积概率函数为:
在模型中,α为位置参数,β为尺度参数。
2)当广义极值分布模型中的k参数的值大于0且较大(在预设的0附近的范围值外)时,则被检测管道的腐蚀数据分布符合Frechet分布。Frechet分布,其累积概率函数为:
在模型中,α为位置参数,β为尺度参数,γ为形状参数。X表示随机变量的最大值。
3)当广义极值分布模型中的k参数的值小于0且较小(在预设的0附近的范围值外)时,则被检测管道的腐蚀数据分布符合Weibull分布。Weibull分布,其累积概率函数为:
在模型中,具体地,α为位置参数,β为尺度参数,γ为形状参数。
在本实施例中参数k通过将广义极值模型与被检测管道所得腐蚀数据进行拟合得到。即:应用式(4)对步骤S01最后得到的曲线进行拟合,即可得到式(4)中的参数k,从而确定选择哪个分布模型是最合适的。k参数确定后,即可确定被检测管道腐蚀数据所符合的具体模型分布。
步骤S03,分布模型参数估计。
对腐蚀数据分布模型中涉及的参数进行估计时,采用概率权重矩法,具体地,如下:
其中,ωr表示概率权重矩;F表示累积概率;r表示实数;dF表示积分时的变量;
基于概率权重矩法得到位置参数α及尺度参数β的确定方法,两个参数的表达式如下:
其中,Γ表示Gamma函数;b0,b1表示无偏估计量。
步骤S04,拟合优度检验。对被检测管道的腐蚀数据进行处理及选取,腐蚀数据分布模型选择后,需要确定所选择的分布模型对于所检测管道腐蚀数据的拟合适用性。具体地,
对选定的腐蚀数据分布模型适用性的判断可通过拟合优度检验进行。通过拟合优度检验的结果确定所选择的分布模型是否满足要求,若不满足要求,则需再选择其他分布模型。对于拟合优度检验需选择合适的方法,本申请中涉及的拟合优度检验方法包括K-S检验方法、作图法、回归法、X2检验及A-D检验等。对于具体被检测管道腐蚀数据,采用一种或几种拟合优度检验方法进行检验。其中:
K-S检验法是指:将需要做统计分析的数据和另一组某一分布的标准数据进行对比,求得它和标准数据之间的偏差的方法。一般在K-S检验中,先计算两组观察数据的累积分布函数,然后计算这两个累积分布函数对应值的差的绝对值中的最大值。最后通过查表以确定该最大值是否落在所要求对应的置信区间内。若该值在对应的置信区间内,说明被检测的数据满足条件,反之亦然。
作图法是指:可以在图中比较直观的比较出经验分布与广义极值分布、Gumbel分布和两参数Weibull分布的拟合程度,通过累积概率分布图可以比较出和经验分布的拟合优度最好的分布形式,通过P-P散点图可以进一步来验证该分布与经验分布的拟合程度最好,P-P散点图判断的标准是,数据点分布在直线y=x附近即可认为该分布与经验分布拟合程度比较好。
步骤S05,腐蚀深度预测。
对被检测管道腐蚀数据实施步骤S01及S02后,确定了腐蚀数据分布模型。对被检测管道腐蚀数据实施步骤S03后,确定了腐蚀数据分布模型涉及的模型参数。对被检测管道腐蚀数据实施S04后,确定了所选择分布模型的适用性。基于步骤S01-S04的结果,应用所确定的腐蚀数据分布模型预测不同可靠度下管道的极限腐蚀深度。
步骤S06,腐蚀剩余寿命预测。对被检测管道腐蚀数据实施步骤S01至S05后,可根据所确定的管道在不同可靠度下的极限腐蚀深度确定管道在不同可靠度下的腐蚀剩余寿命。
具体地,首先计算得到管道最小允许剩余厚度,如下所示:
ΔA=A-Amin (12)
其中,Amin为最小允许剩余厚度,mm;p为管道运行压力,MPa;D为管道外径,mm;SYMS为管线钢屈服强度,MPa;A为管道壁厚,mm;ΔA为管道临界腐蚀深度,mm。
计算得到管道最小允许剩余厚度后结合不同可靠度下的腐蚀深度值,确定不同可靠度下的腐蚀进展曲线,如下:
x=c·tn (13)
其中,x—管道腐蚀深度值,mm;c—系数;t—管道结构开始破坏的使用时间,年。n为时间常数,由腐蚀体系等因素决定,一般取0.3-0.5。
计算得到不同可靠度下的腐蚀进展曲线中的参数后,进一步确定管道的腐蚀剩余寿命,即管道结构开始破坏的适用时间。对被检测管道腐蚀数据实施步骤S01-S06后,可计算得到管道在不同可靠度下的剩余寿命。
请参阅图2,本发明的公共管廊管道腐蚀剩余寿命评估***,包括:数据获取模块1、状态分析模块2和状态预测模块3。
采用公共管廊管道无损检测技术对目标管段进行检测,确定管道的腐蚀状况,得到各类数据。数据获取模块1获取目标管段不同位置的管壁腐蚀减薄数据,统计整理腐蚀数据,并传输给状态分析模块2。
状态分析模块2采用广义极值分布模型针对管道腐蚀状况分析结果,按被检测目标管道的不同管段进行分类整理。
状态预测模块3根据状态分析模块2所得结果对管廊管道的状态进行预测,预测管道腐蚀的发展状况,获得管道在不同可靠度下剩余使用寿命。
为提高评估***的应用效率,对数据获取模块1、状态分析模块2和状态预测模块3采用程序化方式进行。
下面结合一个具体实施例,对本申请的实施具体步骤及实际应用效果进行详细说明。
管廊管道腐蚀数据处理及选取:所进行分析的管廊管道编号为管道1,2,3。
管廊管道1的基本参数为:设计压力2.1MPa,设计温度65℃,工作压力0.8MPa,工作温度12℃,输送介质苯乙烯,管材20#钢,管长1447.5m,DN150x7.1,腐蚀裕量2.6mm,管廊管道1走势图如图3所示。
管廊管道2的基本参数为:长为1100m,DN130x6.9,20#钢,腐蚀裕量2.5mm。设计压力2.2MPa,工作压力0.9-2.0MPa,工作温度15℃,设计温度为70℃,输送介质为氢气,实际使用时间为4.5年,管廊管道2走势图如图4所示。
管廊管道3的基本参数为:输送化工园区污水,设计压力6.0MPa,设计温度80℃,管材20#钢,管长5496.6m,DN100x7.1,实际使用时间为10.3年,管廊管道3走势图如图5所示。
管廊管道1,2,3的腐蚀深度统计汇总如表1,2,3所示。
局部腐蚀深度区间范围 | 频数 |
1-2.5 | 16 |
2.6-2 | 8 |
4.1-4.5 | 3 |
4.6-5.0 | 1 |
5.1-5.5 | 4 |
>5.5 | 1 |
表1管廊管道1腐蚀深度统计汇总表
局部腐蚀深度区间范围 | 频数 |
0.5-1 | 6 |
2.1-2.5 | 17 |
2.6-2 | 5 |
>2 | 5 |
表2管廊管道2腐蚀深度统计汇总表
局部腐蚀深度区间范围 | 频数 |
0.5-1 | 7 |
2.1-2.5 | 26 |
表3管廊管道3腐蚀深度统计汇总表
管廊管道腐蚀数据分布模型确定:
对管廊管道1,2,3腐蚀深度统计结果进行分布模型及参数确定,首先判断管道1局部腐蚀深度数据是否服从Gumbel分布。
对式(5)两端同时进行2次取对数,得到:
以局部腐蚀深度为横坐标,为纵坐标,将各数据点画在直角坐标系中,如图6,7,8所示。
对于管廊管道1,如图6所示,线性拟合度R2为0.989,约为1,则管廊管道1腐蚀深度概率服从Gumbel极值I型分布。
具对于管廊管道2,如图7所示,线性拟合度R2为0.987,约为1,则管廊管道2腐蚀深度概率服从Gumbel极值I型分布。
对于管廊管道3,如图8所示,线性拟合度为R2为0.974,约为1,则管廊管道3腐蚀深度概率服从Gumbel极值I型分布。
管廊管道在不同可靠度下的腐蚀深度:
具体地,对于管廊管道1,2,3在不同可靠度下的极限腐蚀深度,如表4所示。
表4不同可靠度下管道腐蚀深度
管廊管道在不同可靠度下的腐蚀剩余寿命评估:
根据式(11)及式(12)计算得到管道最小允许剩余厚度。
下面对管廊管道3的腐蚀剩余寿命评估流程进行详细介绍,管廊管道1和2的腐蚀剩余寿命评估流程相同。
具体地,对于管廊管道3,可靠度为99%下管道的可使用寿命评估:
首先根据管道投入使用的时间和99%可靠度下的预测腐蚀深度,计算腐蚀进展曲线公式(13)中的系数c值。
6.231=c·10.30.3 (14)
由此可计算得出c=3.095,故腐蚀进展曲线公式为:
x=3.095·t0.3 (15)
管道厚度为7.1mm,最小允许剩余厚度为0.54mm,故:
7.1-0.54=3.095·t0.3 (16)
计算得到t=12.23,由此依照管道的运行条件,通过计算可知在可靠度为99%下管道的可使用寿命为12.23年。
具体地,对于管廊管道3,可靠度为99.9%下管道的可使用寿命评估:
首先根据管道投入使用的时间和99.9%可靠度下的预测腐蚀深度,计算腐蚀进展曲线公式中的系数c值。
6.534=c·10.30.3 (17)
由此可计算得出c=3.246,故腐蚀进展曲线公式为
x=3.246·t0.3 (18)
管道厚度为7.1mm,最小允许剩余厚度为0.54mm,故
7.1-0.54=3.246·t0.3 (19)
计算得到t=10.44,由此依照管道的运行条件,通过计算可知在可靠度为99.9%下管道的可使用寿命为10.44年。
具体地,对于管廊管道3,可靠度为99.99%下管道的可使用寿命评估:
首先根据管道投入使用的时间和99.99%可靠度下的预测腐蚀深度,计算腐蚀进展曲线公式中的系数c值。
6.837=c·10.30.3 (20)
由此可计算得出c=3.396,故腐蚀进展曲线公式为
x=3.396·t0.3 (21)
管道厚度为7.1mm,最小允许剩余厚度为0.54mm,故
7.1-0.54=3.396·t0.3 (22)
计算得到t=8.972,由此依照管道的运行条件,通过计算可知在可靠度为99.99%下管道的可使用寿命为8.972年。
上述管廊管道1,2,3,应用本申请提出的腐蚀剩余寿命评估方法,计算得到的腐蚀剩余寿命评估结果,如表5所示。
表5管廊管道1,2,3在不同可靠度下的腐蚀剩余寿命评估结果
以上实施例仅供说明本发明之用,而非对本发明的限制,有关技术领域的技术人员,在不脱离本发明的精神和范围的情况下,还可以作出各种变换或变型,因此所有等同的技术方案也应该属于本发明的范畴,应由各权利要求所限定。
Claims (3)
1.一种公共管廊管道腐蚀剩余寿命评估方法,其特征在于,包括:
腐蚀数据处理及选取:将被检测管道分成多个管段,对所有管段的腐蚀深度值的最大值进行排序,获取被检测管道的腐蚀深度值的样本序列,以及获取任一腐蚀深度值的累积概率;
选择腐蚀数据分布模型:采用广义极值分布模型并进行参数估计,根据确定的参数选择被检测管道的腐蚀数据分布模型为Gumble分布模型、Frechet分布模型或者Weibull分布模型;
预测腐蚀深度:应用所确定的腐蚀数据分布模型预测不同可靠度下被检测管道的极限腐蚀深度;
预测腐蚀剩余寿命:根据被检测管道在不同可靠度下的极限腐蚀深度确定被检测管道在不同可靠度下的腐蚀剩余寿命;
所述的腐蚀数据处理及选取中,
每个管段的腐蚀深度值按由小到大的顺序进行排序,根据排序结果得到该管段的腐蚀深度值的最大值;
对所有管段的腐蚀深度值的最大值按从小到大的顺序进行排序,形成被检测管道的腐蚀深度值的样本序列:
D={x1,x2,x3,…xi,…,xm};
其中,xi表示被检测管道的第i个腐蚀深度值;m表示m个管段;
对于第i个腐蚀深度值,其累积概率为:
对被检测管道的所得到的所有(xi,F(xi))数据点进行做图,并画出曲线,得到被检测管道腐蚀深度的经验分布图;
所述的选择腐蚀数据分布模型中,
广义极值分布模型的表达式为:
其中,α为位置参数,β为尺度参数,k为广义极值模型分布参数,x表示随机变量;
k参数的值在预设的0附近的范围值内时,则被检测管道的腐蚀数据分布符合Gumble分布,其累积概率函数为:
其中,α为位置参数,β为尺度参数;
k参数的值大于0且在预设的0附近的范围值外时,则被检测管道的腐蚀数据分布符合Frechet分布,其累积概率函数为:
其中,α为位置参数,β为尺度参数,γ为形状参数;X表示随机变量的最大值;
k参数的值小于0且在预设的0附近的范围值外时,则被检测管道的腐蚀数据分布符合Weibull分布,其累积概率函数为:
其中,α为位置参数,β为尺度参数,γ为形状参数;X表示随机变量的最大值;
参数k通过将广义极值分布模型的表达式所述的曲线进行拟合得到;
所述的选择腐蚀数据分布模型中,
参数估计采用概率权重矩法:
其中,ωr表示概率权重矩;F表示累积概率;r表示实数;dF表示积分时的变量;
得到位置参数α及尺度参数β的确定方法,两个参数的表达式如下:
其中,Γ表示Gamma函数;b0,b1表示无偏估计量;
所述的预测腐蚀剩余寿命中,
计算得到管道最小允许剩余厚度,如下所示:
△A=A-Amin
其中,Amin为最小允许剩余厚度;p为管道运行压力;D为管道外径;SYMS为管线钢屈服强度;A为管道壁厚;△A为管道临界腐蚀深度;
结合不同可靠度下的腐蚀深度值,确定不同可靠度下的腐蚀进展曲线,如下:
x=c·tn
其中,x表示管道腐蚀深度值;c表示系数;t表示管道结构开始破坏的使用时间;n为时间常数,取0.3-0.5;
计算得到管道在不同可靠度下的剩余寿命。
2.根据权利要求1所述的公共管廊管道腐蚀剩余寿命评估方法,其特征在于,在选择腐蚀数据分布模型之后,对选定的腐蚀数据分布模型适用性通过拟合优度检验进行判断。
3.根据权利要求2所述的公共管廊管道腐蚀剩余寿命评估方法,其特征在于,所述的拟合优度检验的方法包括K-S检验方法、作图法、回归法、X2检验以及A-D检验。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010383246.7A CN113626970B (zh) | 2020-05-08 | 2020-05-08 | 公共管廊管道腐蚀剩余寿命评估方法及*** |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010383246.7A CN113626970B (zh) | 2020-05-08 | 2020-05-08 | 公共管廊管道腐蚀剩余寿命评估方法及*** |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113626970A CN113626970A (zh) | 2021-11-09 |
CN113626970B true CN113626970B (zh) | 2024-01-09 |
Family
ID=78377289
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010383246.7A Active CN113626970B (zh) | 2020-05-08 | 2020-05-08 | 公共管廊管道腐蚀剩余寿命评估方法及*** |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113626970B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114492980B (zh) * | 2022-01-21 | 2022-09-02 | 中特检深燃安全技术服务(深圳)有限公司 | 一种城镇燃气埋地管道腐蚀风险的智能预测方法 |
CN115345879B (zh) * | 2022-10-18 | 2023-03-24 | 济宁康盛彩虹生物科技有限公司 | 基于图像分析高压釜内壁腐蚀度并预测其使用寿命的方法 |
CN116467857B (zh) * | 2023-03-30 | 2023-09-22 | 淄博市特种设备检验研究院 | 一种基于数据处理的压力管道参数检验管理*** |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102590459A (zh) * | 2012-03-09 | 2012-07-18 | 中国石油天然气股份有限公司 | 一种埋地管道腐蚀检测评价***和方法 |
CN104807966A (zh) * | 2015-04-30 | 2015-07-29 | 上海化学工业区公共管廊有限公司 | 一种管廊管道剩余强度与剩余寿命计算方法 |
CN106372299A (zh) * | 2016-08-29 | 2017-02-01 | 西安建筑科技大学 | 一种基于逆高斯过程和em‑pf的管道剩余寿命预测方法 |
CN108536926A (zh) * | 2018-03-21 | 2018-09-14 | 西安建筑科技大学 | 基于改进自适应gev分布的腐蚀油气管道剩余寿命分析方法 |
CN108875214A (zh) * | 2018-06-19 | 2018-11-23 | 中国计量大学 | 一种基于Gumbel极值分布的储罐剩余寿命预测方法 |
CN109374513A (zh) * | 2018-09-29 | 2019-02-22 | 卢俊文 | 管道外防腐层检测分析*** |
CN110309577A (zh) * | 2019-06-26 | 2019-10-08 | 西安建筑科技大学 | 一种基于im和lmle-bu算法的海底管道剩余寿命预测方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8589087B2 (en) * | 2010-07-28 | 2013-11-19 | General Electric Company | Systems, methods, and apparatus for monitoring corrosion or corrosive contaminants associated with liquid fuel |
-
2020
- 2020-05-08 CN CN202010383246.7A patent/CN113626970B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102590459A (zh) * | 2012-03-09 | 2012-07-18 | 中国石油天然气股份有限公司 | 一种埋地管道腐蚀检测评价***和方法 |
CN104807966A (zh) * | 2015-04-30 | 2015-07-29 | 上海化学工业区公共管廊有限公司 | 一种管廊管道剩余强度与剩余寿命计算方法 |
CN106372299A (zh) * | 2016-08-29 | 2017-02-01 | 西安建筑科技大学 | 一种基于逆高斯过程和em‑pf的管道剩余寿命预测方法 |
CN108536926A (zh) * | 2018-03-21 | 2018-09-14 | 西安建筑科技大学 | 基于改进自适应gev分布的腐蚀油气管道剩余寿命分析方法 |
CN108875214A (zh) * | 2018-06-19 | 2018-11-23 | 中国计量大学 | 一种基于Gumbel极值分布的储罐剩余寿命预测方法 |
CN109374513A (zh) * | 2018-09-29 | 2019-02-22 | 卢俊文 | 管道外防腐层检测分析*** |
CN110309577A (zh) * | 2019-06-26 | 2019-10-08 | 西安建筑科技大学 | 一种基于im和lmle-bu算法的海底管道剩余寿命预测方法 |
Non-Patent Citations (1)
Title |
---|
基于GRA-FOA-ELM的腐蚀管道剩余寿命预测;骆正山 等,;《消防科学与技术》;第39卷(第4期);529-532 * |
Also Published As
Publication number | Publication date |
---|---|
CN113626970A (zh) | 2021-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113626970B (zh) | 公共管廊管道腐蚀剩余寿命评估方法及*** | |
CN110020399B (zh) | 一种管道内腐蚀评价位置的确定方法 | |
CN112883538A (zh) | 一种埋地原油管道的腐蚀预测***及方法 | |
US10942083B2 (en) | System and method for the prediction of leakage in a pipeline | |
CN115077618A (zh) | 一种核级合金钢弯头的质量检测方法及*** | |
CN111553316A (zh) | 一种检测核级管道汽蚀故障的方法 | |
CN116617616A (zh) | 一种消防管道监测方法、***、终端设备及存储介质 | |
CN114528700B (zh) | 一种含腐蚀坑油管剩余强度的确定方法 | |
CN114201840B (zh) | 管段应力腐蚀开裂风险识别方法、装置、设备及介质 | |
CN110490493B (zh) | 一种基于分安全系数的含缺陷奥氏体不锈钢压力容器安全评定方法 | |
Wright et al. | Evaluation of Corrosion Growth Prediction Methodologies Using Burst Pressure Comparisons From Repeated In-Line Inspections | |
KR20190083131A (ko) | 위험도를 기반으로 한 배관감육 관리방법 | |
Anghel et al. | Risk assessment for vessels affected by corrosion | |
CN118243558A (zh) | 一种采用包络线评价高分子材料降解性能的方法及*** | |
KR102627849B1 (ko) | 감육 배관 우선 검사 라인 선정 시스템 및 방법 | |
CN112748005B (zh) | 一种铝合金型材薄弱位置的测算方法 | |
CN118225999B (zh) | 一种大直径压力钢管焊缝无损检测方法及*** | |
RU2796240C1 (ru) | Способ определения степени износа оборудования под воздействием коррозии | |
Amaya-Gómez et al. | Spatial and temporal corrosion degradation modelling with a Lévy process based on ILI (in-line) inspections | |
CN113268827B (zh) | 一种基于等效损伤路径含裂纹压力容器剩余寿命预测方法 | |
CN115166024B (zh) | 一种磁记忆检测异质金属涂层结合面损伤程度的方法 | |
CN111144002B (zh) | 一种供水管网灰口铸铁管爆管预测的方法 | |
Tomar et al. | Qualification of ILI performance in accordance with API 1163 and the potential impact for management of pipeline integrity | |
AMAYA-GÓMEZ | Multiscale framework for corrosion and reliability assessment of inspected pipelines | |
CN114910380A (zh) | 一种用于确定井下油管剩余使用寿命的方法 |
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 |