一种山地管道裂纹缺陷剩余强度的评价方法与流程

1.本发明涉及管道评价领域，特别地，涉及一种山地管道裂纹缺陷剩余强度的评价方法。


背景技术：

2.管道途经人口密集区、高后果区的管道发生事故，可能会导致重大的人员伤亡和巨大的经济损失。因此油气管道的可靠性对社会稳定和经济发展具有极大影响。
3.西南山区管道途径区域山高谷深、河流纵横、林木密布、地质灾害频发、地震发育密集，导致穿跨越管段、沉管管段较多，同时管道易受外部环境的影响产生裂纹等管道损伤，管道失效发生风险较大，安全生产压力大。因此，对含裂纹管道的精确评价对管道安全运行具有重要意义。


技术实现要素：

4.本发明的目的在于解决现有技术存在的上述不足中的至少一项。例如，本发明的目的之一在于对易发生地质灾害区域的管段进行更加精确的评价。
5.为了实现上述目的，本发明提供了一种山地管道裂纹缺陷剩余强度的评价方法。
6.所述方法可包括以下步骤：确定山地管道裂纹缺陷评价模型；确定山地管道安全系数的修正因子；结合所述修正因子，确定山地管道安全系数的范围；根据所述评价模型和所述安全系数的范围，进行山地管道裂纹缺陷剩余强度的评价。
7.进一步地，所述评价模型可包括：
[0008][0009]
其中，kr＝k
ι
/k
mat
为韧性比，ki为应力强度因子，k
mat
为材料的断裂韧性；lr＝σ
ref
/σy为载荷比，σ
ref
为参比应力，σy为材料的屈服强度；l
rmax
为评估曲线的截止线，σu为材料的抗拉强度。
[0010]
进一步地，所述确定修正因子的步骤可包括：选择风险因素组成数据集；利用mipca模型进行归一化处理；采用waspas方法对裂纹缺陷的风险因素进行综合分析；根据综合分析结果，确定所述修正因子。
[0011]
进一步地，风险因素是与裂纹缺陷相关的风险因素。
[0012]
进一步地，利用下式进行所述综合分析：
[0013]
其中，qi为第i个观测点的综合评价值，λ为的贡献度，λ＝0,
…
,1， wj为第j
个属性集的权重，为标准化后的第j个属性集cj中的第i个观测点的分值。
[0014]
进一步地，根据下式来确定所述修正因子：
[0015]
其中，a为修正因子，min是所有观测点的综合评价值中的最小值，max 是所有观测点的综合评价值中的最大值。
[0016]
进一步地，所述根据下式来确定：
[0017]
其中，c
ij
表示第j个属性集cj中的第i个观测点的分值。
[0018]
进一步地，所述进行归一化处理的步骤可包括：计算风险因素的互信息矩阵；计算所述互信息矩阵的特征值，并对其进行排列，找到对应的特征向量；计算互信息的主成分；计算主成分的贡献率，进而确定特征的维度。
[0019]
进一步地，根据下式确定所述安全系数的范围：
[0020][0021]
其中，sf为安全系数，p为设计压力，ph为最小水压试验压力，maop 为最大允许操作压力，p0为操作压力，f为设计系数，a为修正因子。
[0022]
进一步地，所述进行山地管道裂纹缺陷剩余强度评价的步骤可包括：根据所述安全系数修正所述评价模型，利用修正后的评价模型来进行所述评价。
[0023]
进一步地，根据所述安全系数修正所述评价模型中的应力强度因子，其中，k
isf
＝ki×
sf，k
isf
为修正后的应力强度因子，sf为安全系数，ki为修正前的应力强度因子。
[0024]
进一步地，修正后评价模型中的kr＝k
ιsf
/k
mat
。
[0025]
进一步地，所述山地管道所处的山地条件包括：高丘陵地貌敷设占75～ 80％，平原谷地地貌敷设占20％～25％。例如，高丘陵地貌敷设占78％，平原谷地地貌敷设占22％。
[0026]
与现有技术相比，本发明的有益效果包括以下内容中的至少一项：
[0027]
(1)本发明通过增加山地管道安全系数修正因子来进行更严苛的约束，能够较为准确的反映山地管道在裂纹缺陷影响下的管道状况。
[0028]
(2)本发明通过增加山地安全系数修正因子，能够有效提高易发生地质灾害区域管段的评价精确性。
[0029]
(3)本发明针对含缺陷裂纹管道剩余强度所进行的精确评价，对延长管道服役寿命具有指导意义。
附图说明
[0030]
通过下面结合附图进行的描述，本发明的上述和其他目的和/或特点将会变得更加清楚，其中：
[0031]
图1示出了本发明的山地管道裂纹缺陷剩余强度的评价方法的一个流程示意图。
[0032]
图2示出了本发明的失效评估图的一个示意图。
具体实施方式
[0033]
在下文中，将结合示例性实施例来详细说明本发明的一种山地管道裂纹缺陷剩余强度的评价方法。
[0034]
山地管道的裂纹缺陷问题非常突出，山地管道的裂纹型缺陷主要是指平面缺陷，根部半径较为尖锐；而裂纹状缺陷包括平面缺陷、焊缝上的未熔合和未焊透、尖锐的沟槽状局部腐蚀以及环境开裂相关的分枝状裂纹。裂纹缺陷特征变化较大，这与裂纹的成因、材料以及环境有关。裂纹可以在管道外表面萌生，在长度和表面方向发生扩展。沿表面扩展的方向垂直于环向应力，导致裂纹在管道轴向连接在一起。考虑到西南管道所受周围复杂载荷的特性，在确定裂纹型缺陷评价评价模型的同时，选择安全系数计算方式，增加山地安全系数修正因子，用于对易发生地质灾害区域管段进行更加精确的评价。
[0035]
示例性实施例1
[0036]
图1示出了本发明的山地管道裂纹缺陷剩余强度的评价方法的一个流程示意图。如图1所示，所述山地管道裂纹缺陷剩余强度的评价方法可包括以下步骤：
[0037]
步骤s10：确定裂纹缺陷失效应力评价模型。
[0038]
对于裂纹型缺陷的评价，本发明可以选用api 579-2007评价方法，而且其分类保持不变。评价级别均为1级、2级和3级。在1级评价中，依据管道的规格、缺陷的位置和管道的工作温度来选择相应的评价曲线对缺陷的可接受性进行评估。1级、2级和3级三种评价级别互为独立，由于1级评价考虑因素不够全面，3级评价使用过程较为繁琐，为了保证评价的效率及准确性，以及评价对象的适用性及后续安全系数的修正工作，本发明仅针对2级评价级别进行研究。
[0039]
在2级评价中，采用了失效评估图技术。失效评估图如图2所示，评估曲线方程如式为：
[0040][0041]
其中，kr＝k
ι
/k
mat
为韧性比，ki为应力强度因子，与管道承受的压力和缺陷的尺寸大小有关，k
mat
为材料的断裂韧性；
[0042]
lr＝σ
ref
/σy为载荷比，σ
ref
为参比应力，σy为材料的屈服强度；
[0043]
为评估曲线的截止线，σu为材料的抗拉强度。
[0044]
结合图2，当以(lr,kr)为坐标的评估点o落在评估曲线左下方时，则在管道当前工作压力下，缺陷可以接受。否则，缺陷不能接受。
[0045]
步骤s20：确定山地安全系数的修正因子。
[0046]
在该步骤中，主要是利用mipca模型和waspas方法来进行，其中， mipca用于筛选出与管道划痕凹陷缺陷相关的风险因素；waspas根据筛选出的风险因素对不同观测点的综合评价值进行确定。
[0047]
如图1所示，该步骤可包括：选择风险因素组成数据集；利用mipca模型进行归一化处理；采用waspas方法对裂纹缺陷的风险因素进行综合分析；根据综合分析结果，确定所述修正因子。
[0048]
在本实施例中，主成分分析(principal component analysis,pca)是一种将多个变量通过线性变换以选出较少个数重要变量的一种多元统计分析方法。但是由于在实际数据集中，变量之间的关系不仅是线性关系，还存在非线性关系，因此，互信息方法就被引入作为特征处理的一种新的方式。由于互信息是以信息理论为基础，具有反映变量之间全部信息的优点，因此互信息与主成分分析的结合具有更好的变量选择优势，也就形成了mipca模型。 mipca的计算过程包括：
[0049]
(1)假设共选择出p个风险因素(即自变量)组成数据集x，x＝[x1,x2,
…
,x
p
]，则他们之间的互信息矩阵为：
[0050][0051]
(2)计算互信息矩阵的特征值，并按照降序对其进行排列，找到对应的特征向量，公式如下：
[0052]
b'∑i
x
b＝λ
ꢀꢀꢀꢀꢀ
(2)
[0053]
其中，b(b1,b2,
…
,b
p
,)是一个关于特征向量bβ的矩阵，b'是b的转置，λ(μ1,μ2,
…
,μ
p
,)是一个包含特征值μ的对角阵。
[0054]
(3)计算互信息的主成分，公式如下：
[0055]
z＝b'x
ꢀꢀꢀꢀꢀ
(3)
[0056]
其中，z(z1,z2,
…
,z
p
)是一个关于主成分的矩阵，zk＝b'kxk(k＝1,2,
…
,p)。
[0057]
(4)计算特征的维度m，公式如下：
[0058][0059]
其中，σk是第k个主成分的贡献率；
[0060][0061]
其中，δk是前k个主成分的贡献率之和，一般来说，当δk达到85％～95％时，m＝k。
[0062]
在本实施例中，waspas主要包括三个步骤，分别代表三个优化：
[0063]
(1)可以实现指标的精确评估，计算公式如下：
[0064][0065]
其中，wj为第j个属性集的权重；wj为可认为是上述主成分的贡献度σk，或者，wj可由ahp求得；
[0066]
n表示属性集的总数；q
i(1)
表示第i个观测点的第一个评价值；表示标准化后的第j个属性集cj中的第i个观测点的分值，计算公式如下：
[0067]
[0068]
其中，c
ij
表示第j个属性集cj中的第i个观测点的分值。
[0069]
其中，属性集是数据集经过mipca模型处理后与决策属性相关性较高的条件属性集。
[0070]
(2)能够突出当前数据对模型准确性的贡献度，计算公式如下：
[0071][0072]
其中，q
i(2)
表示第i个观测点的第二个评价值。
[0073]
(3)将(1)与(2)结果相加，从而实现指标评价与数据贡献度的结合，提高评价结果的准确性，计算公式如下：
[0074][0075]
其中，qi为第i个观测点的综合评价值，λ、1-λ分别为q
i(1)
和q
i(2)
的贡献度，λ＝0,
…
,1。
[0076]
在本实施例中，安全系数的修正因子计算公式如下：
[0077][0078]
其中，a是安全系数的修正因子，qi为第i个观测点的综合评价值，min 是所有观测点的综合评价值中的最小值，max是所有观测点的综合评价值中的最大值。
[0079]
步骤s30：确定山地管道安全系数范围。
[0080]
考虑到山地管道所受载荷复杂，山地地理环境复杂，油气管道所受威胁变数更大等影响，本发明根据更加保守的asme b31g-2012来确定安全系数，同时增加山地管道安全系数修正因子进行更严苛的约束，较为准确的反映山地管道在裂纹缺陷影响下管道状况，即：
[0081][0082]
其中，p为设计压力，ph为最小水压试验压力，maop为最大允许操作压力，sf为安全系数，p0为操作压力，f为设计系数，a为修正因子。pf为预测失效压力，为设计压力p与设计系数f的比值，即
[0083]
步骤s40：考虑山地安全系数修正因子的裂纹缺陷管道剩余强度评价。
[0084]
具体地，该步骤可包括：根据步骤s10中的评价模型和步骤s30中的安全系数的范围，进行山地管道裂纹缺陷剩余强度的评价。
[0085]
在本实施例中，在评价模型中，可考虑山地安全系数后的山地管道应力强度因子，即k
isf
＝ki×
sf，其中，k
isf
为考虑了安全系数的应力强度因子， sf为安全系数，ki为考虑安全系数的应力强度因子。
[0086]
尽管上面已经结合示例性实施例及附图描述了本发明，但是本领域普通技术人员应该清楚，在不脱离权利要求的精神和范围的情况下，可以对上述实施例进行各种修改。