一种基于参数化模型的油气管道剩余强度预测方法
【专利摘要】本发明公开了一种基于参数化模型的油气管道剩余强度预测方法。根据已建立的受损管道三维模型，以管道材料抗拉极限强度作为判据对受损管道进行非线性有限元力学仿真分析，获得缺陷几何量对管道剩余强度的影响关系，根据缺陷几何量对管道剩余强度的影响关系构建用于定量化描述管道剩余强度的参数化模型，利用智能优化方法求解参数化模型的未知系数，将检测得到的缺陷几何量代入到已获得未知系数的参数化模型中求得管道剩余强度，完成对油气管道剩余强度的预测。本发明有效解决海底受损油气管道的安全性预测问题，可靠性高，可有效地用于海底油气受损管道的剩余强度分析，为管道的安全性预测提供理论依据。
【专利说明】
一种基于参数化模型的油气管道剩余强度预测方法
技术领域
[0001] 本发明涉及了一种管道检测预测方法，尤其是涉及了一种基于参数化模型的油气 管道剩余强度预测方法。
【背景技术】
[0002] 油气管道是油气开发的重要组成部分，是油气输送的主要手段，被认为是油气生 产系统中的生命线。由于受环境侵蚀、外部载荷作用等因素的影响势必导致油气管道系统 的损伤积累和抗力衰减，油气管道破裂泄漏的风险也悄然增加，防范油气管道的泄漏已成 为当务之急。
[0003] 从上世纪70年代起，国外就开展了受损管道极限承载力和失效机理的相关研究， 并制定了一系列的管道安全性预测规范。美国机械工程师协会颁布了具有代表性的管道安 全性预测规范ASME B31G《确定腐蚀管道剩余强度的手册》。随后，针对ASME B31G规范的保 守性等问题，O'Grady等对腐蚀面积和膨胀系数等做了相应的修正，被称为改进版的B31G。 挪威船级社通过对海底腐蚀管道进行了一系列的试验研究，与英国燃气公司联合颁布了腐 蚀管道剩余强度预测规范DNV-RP-F101。为了验证现有规范的可靠性，国内外学者开展了腐 蚀管道失效实验研究，并将试验结果与现有的规范进行比较，结果表明规范DNV-RP-F101和 ASME B31G均过于保守，而改进版的B31G则偏于危险。以上规范准则中所给出的计算模型均 未考虑缺陷宽度对于管道剩余强度的影响，这也是影响现有规范预测结果的影响因素之 〇
【发明内容】
[0004] 针对现有管道安全性预测规范中所存在不足，本发明提出了一种基于参数化模型 的油气管道剩余强度预测方法，能有效解决海底受损油气管道的安全性预测问题。
[0005] 本发明采用的技术方案是，如图1所示：
[0006] 根据已建立的受损管道三维模型，以管道材料抗拉极限强度作为判据对受损管道 进行非线性有限元力学仿真分析，获得缺陷几何量对管道剩余强度的影响关系，根据缺陷 几何量对管道剩余强度的影响关系构建用于定量化描述管道剩余强度的参数化模型，将非 线性有限元力学仿真分析的数据代入到参数化模型中利用智能优化方法求解参数化模型 的未知系数，获得未知系数的参数化模型给出了管道剩余强度与缺陷几何量之间的定量关 系，将检测得到的缺陷几何量代入到已获得未知系数的参数化模型中求得管道剩余强度， 完成对油气管道剩余强度的预测。
[0007] 所述的缺陷几何量对管道剩余强度的影响关系包括缺陷长度对管道剩余强度的 影响关系、缺陷深度对管道剩余强度的影响关系和缺陷宽度对管道剩余强度的影响关系。
[0008] 所述缺陷长度对管道剩余强度的影响关系为指数分布关系。
[0009] 所述缺陷深度对管道剩余强度的影响关系为抛物线分布关系。
[0010]所述缺陷宽度对管道剩余强度的影响关系为指数分布关系。
[0011] 所述参数化模型具体采用以下公式：
[0012]
[0013]其中，Pf表示管道剩余强度，D表示管道外径，t表示管道厚度，f表示缺陷几何量，d 表示缺陷深度，表示管道材料的抗拉极限强度。
[0014] 所述的缺陷几何量f采用以下公式表示：
[0015]
[0016] 其中，k、m和η分别为参数化模型的长度影哬因于、深度影响因子和宽度影响因子， L表示缺陷长度，d表示缺陷深度，w表示缺陷宽度。
[0017] 将所述缺陷几何量f的公式代入参数化模型后获得以下最终的参数化模型公式表 达：
[0018]
[0019] 其中，Pf表示管道剩余强度，D表示管道外径，t表示管道厚度，f表示缺陷几何量，〇b 表示管道材料的抗拉极限强度，k、m和η分别为参数化模型的长度影响因子、深度影响因子 和宽度影响因子，L表示缺陷长度，d表示缺陷深度，w表示缺陷宽度。
[0020] 所述的智能优化方法具体可采用遗传算法。
[0021] 非线性有限元力学仿真分析包括：
[0022] 利用非线性有限元方法研究缺陷长度对管道剩余强度的影响规律并进行统计分 析，分析结果如图2所示，得到的结论是，呈现指数分布；
[0023] 利用非线性有限元方法研究缺陷深度对管道剩余强度的影响规律并进行统计分 析，分析结果如图3所示，得到的结论是，呈现抛物线分布；
[0024] 利用非线性有限元方法研究缺陷宽度对管道剩余强度的影响规律并进行统计分 析，分析结果如图4所示，得到的结论是，呈现指数分布。
[0025] 本发明以管道材料的抗拉极限强度作为判据对受损管道进行非线性有限元力学 分析，目前的有限元数值计算中，通常采用会弹塑性本构模型，该模型未曾考虑材料的硬化 效应，而事实上材料的硬化效应对管道的失效会有很大的影响，为此在有限元计算模型中 引入了Ramberg-Osgood幂硬化应力-应变法则，使得分析结果更加接近于实际。
[0026]本发明经过实验验证，表明该模型比现有的管道预测标准具有更好的可靠性，可 有效地用于海底油气受损管道的剩余强度分析，为管道的安全性预测提供理论依据。
[0027]本发明具有的有益效果是：
[0028]本发明在有限元分析基础上结合相应的边界条件构建了所述参数化模型，综合考 虑了管道上缺陷的各个参数对于管道剩余强度的影响，克服了现有管道剩余强度预测规范 所存在的不足，其预测结果能够更加接近于实际情况，可提高受损管道预测结果的可靠性。 [0029]我国目前有将近6000公里长的海底油气管道，随着海洋油气资源的不断开发，海 底管网规模会日益扩大，国外公司对我国实施技术封锁，只提供检测服务，每公里收费动辄 数万美元。
[0030] 本发明能为我国海底管道的安全使用质量管理领域提供一套有效的、具有自主知 识产权的复杂载荷下海底受损油气管道剩余强度分析方法。不仅能在管道检测分析领域中 得到直接应用，而且相关共性技术进一步推广可方便延伸到其它对象应用，拓宽应用范围， 创造更大的社会和经济效益。
【附图说明】
[0031] 图1是本发明总体流程框图。
[0032] 图2是本发明失效载荷随缺陷长度变化示意图。
[0033] 图3是本发明失效载荷随缺陷深度的变化示意图。
[0034] 图4是本发明失效载荷随缺陷宽度变化示意图。
[0035] 图5是本发明模型参数k的最优解进化过程示意图。
[0036] 图6是本发明模型参数m的最优解进化过程示意图。
[0037] 图7是本发明模型参数η的最优解进化过程示意图。
[0038] 图8是本发明失效载荷预测结果比较示意图。
【具体实施方式】
[0039]以下结合附图和实施例对本发明作进一步说明。
[0040]本发明的实施例如下：
[0041] 本实施例采用的对象为受损的油气管道。
[0042] 假定在无缺陷状态下管道的失效载荷为Ρ，有缺陷状态下受损管道的失效载荷为 Pf，利用有限元计算方法分析了受损管道失效载荷Pf随着缺陷长度的变化情况，分析结果如 图2所示。从图中可以看出，随着缺陷长度的增加，失效载荷会随之变小，而且随着缺陷深度 的增加，缺陷长度对失效载荷的影响更为明显。
[0043] 失效载荷随管道缺陷深度的变化情况如图3所示，从图中可以看出，缺陷的深度对 失效载荷的影响较为明显，存在一定的非线性关系。为进一步分析缺陷宽度对于失效载荷 的影响程度，通过计算得到如图4所示的结果，从图中可以看出，缺陷宽度与失效载荷之间 也存在一定的非线性关系。
[0044] A)构建参数化模型：
[0045] 因此，将无缺陷状态下管道所能承受的载荷作为载荷上限：
[0046]
( 1 )
[0047] 其中，D表示管道外径，t表示管道厚度，〇b表示管道材料的抗拉极限强度；
[0048] 再将受损管道最小剩余壁厚作为全壁厚计算得到载荷下限：
[0049]
(2)
[0050] 其中，d表示缺陷深度。
[0051]根据极限载荷的上下限构造出本发明受损管道的失效载荷的参数化模型：
[0052]
(3)
[0053]式中，缺陷几何量函数f中包含了缺陷长度、缺陷深度等缺陷几何量信息。
[0054]结合失效载荷的实际情况，确定失效载荷模型的边界条件：当管道上无缺陷时，失 效载荷应等于载荷上限pmax;而当缺陷轴向长度Ld很大或者缺陷深度接近于壁厚时，缺陷几 何量函数等于〇,失效载荷接近于载荷下限p min。
[0055] 根据以上分析，将缺陷几何量函数f定义为：
[0056]
(4)
[0057] 其中，L表示缺陷长度，d表示缺陷深度，w表示缺陷宽度，D表示管道外径，t表示管 道厚度，k、m和η分别为参数化模型的长度影响因子、深度影响因子和宽度影响因子，可将有 限元计算的结果代入模型中，对模型中参数值进行辨识。
[0058] 将公式1、公式2和公式4代入公式3可得：
[0059]
：5)
[0060] Β)求解参数化模型：
[0061] 经典的非线性寻优方法如牛顿迭代法等采用了梯度下降的方式进行求解，此类方 法的局部搜索能力较强，全局搜索能力较弱，因此只有在初始值与真实值较为接近的情况 下，才能得到较好的迭代结果，然而在模型的参数辨识过程中，理想初始值的获取较为困 难。
[0062] 与传统的搜索方法不同，智能优化方法具有极强的全局搜索能力，比如遗传方法 是从随机产生的初始解开始搜索，通过一定的选择、交叉、变异等操作逐步迭代以产生新的 解，采用概率化的寻优方法，能够自适应地调整搜索方向，因此在多变量优化问题中优于传 统的迭代方法。群体中的每个个体代表问题的一个解，在方法中表示为染色体，染色体的好 坏通过使用度值来衡量，以适应度作为指标从上一代中选择优秀个体，通过交叉、变异等操 作形成下一代群体。经过若干代的进化之后，方法收敛于最优的染色体，即得到问题的最优 解。
[0063]以欧式距离最小作为目标函数，
[0064]
(6)
[0065] 式中，Pf为有限元方法计算得到的管道失效载荷，Pff是参数化模型计算得到的管 道失效载荷，N为所有限元方法计算得到的失效载荷数据个数，N表示参与计算的数据个数， i表示参与计算的数据序列，表示第i个有限元计算得到的管道失效载荷数据，表示第 i个参数模型计算得到的管道失效载荷。
[0066] 针对传统遗传方法容易早熟收敛于局部最优解或者不收敛的缺点，对其交叉操作 和变异操作过程进行一些改进。当种群中个体的适应度值比种群的平均适应度值低时，表 明该个体的性能较差，此时的交叉概率和变异概率应该取较大的值，反之亦然，因此需要根 据个体的特性选取合适的交叉概率和变异概率值与之对应，交叉概率P。和变异概率p m自适 应取值额表示为：
[0067] (7)
[0068] (8)
[0069] 式中，favg表示每一代群体的适应度平均值，fmax表示每一代群体中最大的适应度 值;f'表示参加交叉运算的两个个体中适应度较大的值;f表示变异个体的适应度值。p cl、 pc2分别表示交叉概率的上下极值，pml、pm2分别表示变异概率的上下极值，设定参数p cl = 0 · 9，Pc2 = 0 · 6，Pml = 0 · 1，Pm2 = 0 · 01 〇
[0070] 设定个体数目为180，最大遗传代数100，为提高收敛速度，采用十进制编码方式， 并以式(6)作为适应度函数，通过自适应遗传方法对模型参数进行优化求解，参数的最优解 进化过程如图5、图6、图7所示，参数k的最优解为-0.2976，参数m的最优解为0.3872，参数k 的最优解为-1.6852,将得到的参数值代入式(5)所示的失效载荷模型中，将模型计算所得 的失效载荷与实验实测值以及ASME B31G中所列的失效载荷公式的计算结果进行比较，计 算结果如图8所示，从图中可看出，模型所得到的计算结果相比较ASME B31G更接近于实测 值。
[0071] 上述【具体实施方式】用来解释说明本发明，而不是对本发明进行限制，在本发明的 精神和权利要求的保护范围内，对本发明作出的任何修改和改变，都落入本发明的保护范 围。
【主权项】
1. 一种基于参数化模型的油气管道剩余强度预测方法，其特征在于:根据已建立的受 损管道Ξ维模型，W管道材料抗拉极限强度作为判据对受损管道进行非线性有限元力学仿 真分析，获得缺陷几何量对管道剩余强度的影响关系，根据缺陷几何量对管道剩余强度的 影响关系构建用于定量化描述管道剩余强度的参数化模型，利用智能优化方法求解参数化 模型的未知系数，将检测得到的缺陷几何量代入到已获得未知系数的参数化模型中求得管 道剩余强度，完成对油气管道剩余强度的预测。2. 根据权利要求1所述的一种基于参数化模型的油气管道剩余强度预测方法，其特征 在于:所述的缺陷几何量对管道剩余强度的影响关系包括缺陷长度对管道剩余强度的影响 关系、缺陷深度对管道剩余强度的影响关系和缺陷宽度对管道剩余强度的影响关系。3. 根据权利要求2所述的一种基于参数化模型的油气管道剩余强度预测方法，其特征 在于:所述缺陷长度对管道剩余强度的影响关系为指数分布关系。4. 根据权利要求2所述的一种基于参数化模型的油气管道剩余强度预测方法，其特征 在于:所述缺陷深度对管道剩余强度的影响关系为抛物线分布关系。5. 根据权利要求2所述的一种基于参数化模型的油气管道剩余强度预测方法，其特征 在于:所述缺陷宽度对管道剩余强度的影响关系为指数分布关系。6. 根据权利要求1所述的一种基于参数化模型的油气管道剩余强度预测方法，其特征 在于:所述参数化模型具体采用W下公式：其中，Pf表示管道剩余强度，D表示管道外径，t表示管道厚度，f表示缺陷几何量，d表示 缺陷深度，Ob表示管道材料的抗拉极限强度。7. 根据权利要求6所述的一种基于参数化模型的油气管道剩余强度预测方法，其特征 在于:所述的缺陷几何量f采用W下公式表示：其中，k、m和η分别为参数化模型的长度影响因子、深度影响因子和宽度影响因子，L表 示缺陷长度，d表示缺陷深度，W表示缺陷宽度。8. 根据权利要求7所述的一种基于参数化模型的油气管道剩余强度预测方法，其特征 在于:将所述缺陷几何量f的公式代入参数化模型后获得W下最终的参数化模型公式表达：其中，Pf表示管道剩余强度，D表示管道外径，t表示管道厚度，f表示缺陷几何量，Ob表示 管道材料的抗拉极限强度，k、m和η分别为参数化模型的长度影响因子、深度影响因子和宽 度影响因子，L表示缺陷长度，d表示缺陷深度，W表示缺陷宽度。
【文档编号】G06F17/50GK105975678SQ201610287151
【公开日】2016年9月28日
【申请日】2016年5月4日
【发明人】周闻青, 周红明, 茅振华, 余松青, 叶欣, 陈挺, 卢歆
【申请人】浙江省计量科学研究院