CN113609640B - 一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法 - Google Patents

一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法 Download PDF

Info

Publication number
CN113609640B
CN113609640B CN202110726421.2A CN202110726421A CN113609640B CN 113609640 B CN113609640 B CN 113609640B CN 202110726421 A CN202110726421 A CN 202110726421A CN 113609640 B CN113609640 B CN 113609640B
Authority
CN
China
Prior art keywords
pipeline
corrosion
current
subway
buried pipeline
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
Application number
CN202110726421.2A
Other languages
English (en)
Other versions
CN113609640A (zh
Inventor
王承涛
李威
王禹桥
杨雪峰
许少毅
范孟豹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology CUMT
Original Assignee
China University of Mining and Technology CUMT
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology CUMT filed Critical China University of Mining and Technology CUMT
Priority to CN202110726421.2A priority Critical patent/CN113609640B/zh
Publication of CN113609640A publication Critical patent/CN113609640A/zh
Application granted granted Critical
Publication of CN113609640B publication Critical patent/CN113609640B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/04—Ageing analysis or optimisation against ageing

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)
  • Pipeline Systems (AREA)

Abstract

本发明公开了一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法，适用于地铁维护使用。包其括杂散电流干扰分布计算模块、坐标转换模块、缺陷极限尺寸计算模块和累计腐蚀深度实验及分布计算模块；通过充分模拟杂散电流腐蚀作用变化特点并结合电化学实验结果，能够通过计算沿管线纵向方向不同位置的腐蚀缺陷深度预测在当前腐蚀环境下的管线剩余寿命，对于准确评估和预测管线使用状态、正确确定管线维护策略、保证地铁线路旁长输管线的安全服役具有重要的工程意义。

Description

一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法
技术领域
本发明涉及一种埋地管线剩余寿命预测方法，尤其适用于一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法。
背景技术
在天然气、石油和城市燃气等高危可燃介质输送管线靠近或跨越地铁线路时，会受到泄漏杂散电流的电化学腐蚀影响。这一类电化学腐蚀会导致管壁减薄，甚至穿孔，极大地威胁此类埋地金属管线的运行寿命。由于天然气、石油和城市燃气埋地管线在泄漏时极易引起爆炸、爆燃事故，尤其是对于人口密集的城市区域而言，此类事故所造成的因此有效的针对地铁杂散电流干扰范围内的石油管线进行剩余寿命预测，对于提高天然气、石油和城市燃气管线的安全运行具有极大的工程实际意义。
受腐蚀管线的剩余寿命预测具有重要的工程意义，为及时维修替换受损管线提供了重要参考信息，能够有效降低管线泄漏事故风险，提高油气管网系统的运输可靠性。现有内检测技术只能获取当前状态下的管线壁厚信息，难以对壁厚减薄的未来发展趋势进行精准预测，因而亟需一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法。
发明内容
针对现有技术的不足之处，提供一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法，能够结合具体环境下的腐蚀过程，有效确定性地预测埋地管线在受杂散电流腐蚀作用下的剩余服役时间，
为实现上述目的，本发明的一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法，利用预测系统进行埋地管线剩余寿命预测，所述的预测系统包括杂散电流干扰分布计算模块、坐标转换模块、缺陷极限尺寸计算模块和累计腐蚀深度实验及分布计算模块，其中：
杂散电流干扰分布计算模块：用以计算在地铁机车牵引电流和位置不断变化下的动态干扰，对牵引电流曲线进行离散化处理以近似计算单位腐蚀时间Δt下对埋地管的腐蚀深度；用于获得埋地管线在单位腐蚀时间Δt内恒定杂散电流流入所造成的电位偏移分布和因杂散电流所造成土壤电位的偏移；
坐标转换模块：当地铁隧道与埋地管线互相不为平行关系时，用以将隧道沿长度方向上的离散坐标转换至埋地管线沿长度方向上的离散坐标；
缺陷极限尺寸计算模块：通过管线内压、管线直径、管线壁厚、体积扩张系数、管线流动压力、管线的屈服极限和强度极限计算埋地管线极限壁厚Hl；
累计腐蚀深度实验及分布计算模块：用以计算杂散电流作用下累计腐蚀缺陷深度在埋地管线上的分布，在△t离散时间内考虑自然腐蚀作用与杂散电流腐蚀作用的累积腐蚀深度Hd为累积腐蚀缺陷在管线上分布的最小值通过比对管线极限壁厚Hl即可预测的管线剩余寿命。
首先利用杂散电流干扰分布计算模块计算在地铁机车牵引电流和位置不断变化下的动态干扰，对牵引电流曲线进行离散化处理以近似计算单位腐蚀时间Δt下对埋地管的腐蚀深度，从而获得埋地管线在单位腐蚀时间Δt内恒定杂散电流流入所造成的电位偏移分布和因杂散电流所造成土壤电位的偏移；然后利用坐标转换模块将隧道沿长度方向上的离散坐标转换至埋地管线沿长度方向上的离散坐标；通过缺陷极限尺寸计算模块计算埋地管线的极限壁厚；最后利用累计腐蚀深度实验及分布计算模块计算杂散电流动态干扰下累计腐蚀缺陷深度在埋地管线上的分布，在离散时间△t内考虑自然腐蚀作用与杂散电流腐蚀作用的累积腐蚀深度Hd为累积腐蚀缺陷在管线上分布的最小值通过比对当前检测到的腐蚀深度即可预测的管线剩余寿命。
离散化处理具体为对地铁机车牵引曲线f(t)和地铁机车位置曲线G(t)进行离散化处理：首先通过地铁监测系统数据获得地铁机车牵引电流曲线为f(t)，地铁机车位置曲线为G(t)，获得地铁隧道中的固定参数为：钢轨纵向电阻、排流网纵向电阻、埋地管线纵向电阻、钢轨对排流网过渡电阻、排流网对管线过渡电阻，其中通过电阻网络分布模型计算得到埋地管线内电流计算模型I3(·)，埋地管线内电流计算模型I3(·)为f(t)和G(t)的函数，获取一个区间周期内的地铁机车运行时间t，利用电阻网络模型：{I3(f(t),G(t)，S1’),I3(f(t),G(t)，S2’),…,I3(f(t),G(t)，Si’),…,I3(f(t),G(t)，Sn’)}＝I3(f(t),G(t)，S)，获得沿埋地管线长度方向上的杂散电流分布，S’表示埋地管线长度方向上的绝对任意位置，S1’,S2’,…,Si’,…,Sn’分别表示在隧道坐标系上的绝对坐标；单位腐蚀时间Δt与一个区间周期内地铁机车运行时间t的关系为：nΔt＝t，n为将地铁机车运行时间t离散为Δt的个数。
离散化腐蚀时间Δt的个数n的具体计算方法为：
初始设定参数：将收敛精度p定为0.01；
步骤1：首先根据当前迭代次数ni(ni＝i)，进行如下计算：
步骤2：将ni+1进行累加：ni+1＝ni+1进行迭代计算
步骤3：利用计算L(ni+1):
步骤4：重复步骤1至步骤3，将计算结果L(ni+1)与上一次的计算结果L(ni)做差，若满足L(ni+1)-L(ni)≤p，则判断迭代结束，ni+1则为所求离散化位腐蚀时间Δt的个数n；若L(ni+1)-L(ni)>p，则重复步骤2至步骤3，直至求出满足L(ni+1)-L(ni)≤p，t为地铁机车运行时间，f(t)为地铁机车牵引电流曲线
电阻网络模型为：
组建杂散电流泄漏过程的电阻网络模型，分别包含钢轨纵向导体层、排流网纵向导体层、埋地管线纵向导体层和接地层，四层导体组成了电阻网络，包括钢轨电流I1(S)、排流网电流I2(S)、埋地管线电流I3(S)、钢轨电位U1(S)、排流网电位U2(S)和埋地管线电位U3(S)，利用以下微分方程组求解：
其中，r1为钢轨纵向电阻，r2为排流网纵向电阻，r3为埋地管线纵向电阻，R12为钢轨对排流网过渡电阻，R23为排流网对埋地管线过渡电阻，R34为埋地管线对大地过渡电阻。
地铁隧道泄漏至管线中的杂散电流I3(f(t),G(t),S)通过电阻网络模型计算获得。
坐标转换方法如下：
以地铁隧道轴向方向为原始坐标轴x轴，坐标轴水平面为x-y面，坐标轴垂直面为x-z面，地铁隧道轴向上离散坐标为C1,C2,…,Ci,…,Cm，C1为隧道任意牵引区间范围内起始坐标，C2位隧道长度方向上的第二个坐标，Ci为隧道范围内任意一点离散坐标，Cm为隧道内任意牵引区间范围内的终点坐标，m为离散坐标数量；
设埋地管线与地铁隧道x轴在x-y面上投影夹角为α，埋地管线与地铁隧道x轴在x-z面上投影夹角为β，那么埋地管线在以自身为局部坐标系的相对坐标{S1,S2,…,Si,…,Sm}，管线在隧道坐标系上的绝对坐标{S1’,S2’,…,Si’,…,Sm’}表示为：
{S1’,S2’,…,Si’,…,Sm’}＝{S1cosαcosβ,S2cosαcosβ,…,Sicosαcosβ,…,Smcosαcosβ}；
S1,S2,…,Si,…,Sm即为坐标转换后的埋地管线沿长度方向上离散坐标，S1为区间内所覆盖埋地管线的相对坐标起点，Si为区间内所覆盖埋地管线范围内任意一相对坐标点，Sn为区间内所覆盖埋地管线的相对坐标终点。S1’,S2’,…,Si’,…,Sm’为管线上沿长度方向上的绝对离散坐标，S1为区间内所覆盖埋地管线的绝对坐标起点，Si为区间内所覆盖埋地管线范围内任意一绝对坐标点，Sm为区间内所覆盖埋地管线的绝对坐标终点。
累计腐蚀深度实验及分布计算具体为：
首先通过埋地管线周围土壤的杂散电流腐蚀实验和自然腐蚀实验，两个实验中所用工作电极均与预测目标管线使用同种钢材，杂散电流腐蚀实验中将管线钢工作电极和石墨电极埋置于土壤中，埋置深度均为150mm，通过可控直流电源通电实现杂散电流腐蚀过程，其中管线钢工作电极与石墨电极相距200mm且相对布置，工作电极的尺寸为10mm×10mm×5mm，石墨电极的尺寸为10mm×10mm×2mm，通电电流为is；自然腐蚀实验中将管线钢工作电极埋置于土壤中，埋置深度为150mm，管线钢工作电极的尺寸为10mm×10mm×5mm，模拟管线的自然腐蚀过程；根据杂散电流腐蚀实验和自然腐蚀实验，通过三维数字显微镜在通电电流大小为is作用的情况下确定离散时间Δt内的管线钢工作电极表面平均腐蚀深度ΔH1和自然腐蚀作用下的平均腐蚀深度ΔH2，根据Kn＝ΔH2/ΔH1确定自然腐蚀因子Kn，通过拟合获得通电电流is与Kn的关系Kn(is)，通过拟合通电电流is与ΔH1/Δt获得通电电流与腐蚀速率的关系C(is)。其次，计算在地铁机车运行时间t至t+△t离散时间杂散电流作用下累计缺陷深度在管线上的分布{Ha(S1’),Ha(S2’),…,Ha(Si’),…,Ha(Sm’)}，计算方法为：
Ha(S'i)＝C(Is(f(t+Δt),G(t+Δt),Si’))·Δt
最终地铁机车运行时间t至t+△t离散时间内考虑自然腐蚀作用与杂散电流腐蚀作用的累积腐蚀深度Hd为累积腐蚀缺陷在管线上分布的最小值，即为：
其中，min{·}为取最小值的函数。
重复上述过程，直至其中j为满足 时累积腐蚀深度累加的次数，Hi为预测时刻的管线当前腐蚀深度，则m×△t为所预测的管线剩余寿命。
有益效果：本方法能够有效结合模拟结果与系统参数预测埋地管线在杂散电流腐蚀作用下的剩余服役时间。通过多个模块综合考虑埋地管线当前腐蚀信息和外界干扰信号特征，通过高度模拟腐蚀环境以提高预测精度，能够有效且准确的计算获得埋地管线的剩余寿命，其步骤简单，实施方便，且预测精度高，对于地铁线路附近埋地管线的腐蚀防护具有重要工程意义。
附图说明
图1为本发明受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测系统框图；
图2为杂散电流电阻网络模型示意图；
图3为地铁机车牵引曲线离散化流程图；
图4为地铁隧道与埋地管线坐标在x-y平面上的转换示意图；
图5为地铁隧道与埋地管线坐标在x-z平面上的转换示意图；
图6为埋地管线钢杂散电流腐蚀挂片实验系统图；
图7为埋地管线钢自然腐蚀挂片实验系统图。
具体实施方式
下面结合附图对本发明作进一步详细说明，
本发明的受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法，其特征在于使用的预测系统包括杂散电流干扰分布计算模块、坐标转换模块、缺陷极限尺寸计算模块和累计腐蚀深度分布计算模块，各个模块共同构成了一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法，如图1所示。
杂散电流干扰分布计算模块：用以计算在机车牵引电流和位置不断变化下的动态干扰，对牵引电流曲线进行离散化处理以近似计算单位腐蚀时间Δt下的腐蚀深度；用于获得埋地管线在单位腐蚀时间Δt内恒定杂散电流流入所造成的电位偏移分布和因杂散电流所造成土壤电位的偏移；
坐标转换模块：当地铁隧道与埋地管线互相不为平行关系时，用以将隧道沿长度方向上的离散坐标转换至埋地管线沿长度方向上的离散坐标；
缺陷极限尺寸计算模块：通过管线内压、管线直径、管线壁厚、体积扩张系数、管线流动压力、管线的屈服极限和强度极限计算管线极限壁厚Hl
累计腐蚀深度实验及分布计算模块：用以计算杂散电流动态干扰下累计缺陷深度在管线上的分布，在△t离散时间内考虑自然腐蚀作用与杂散电流腐蚀作用的累积腐蚀深度Hd为累积腐蚀缺陷在管线上分布的最小值通过比对管线极限壁厚Hl即可预测的管线剩余寿命。
组建杂散电流泄漏过程的电阻网络模型，分别包含钢轨纵向导体层、排流网纵向导体层、埋地管线纵向导体层和接地层，四层导体组成了电阻网络，如图2所示，包括钢轨电流I1(S)、排流网电流I2(S)、埋地管线电流I3(S)、钢轨电位U1(S)、排流网电位U2(S)和埋地管线电位U3(S)，利用以下微分方程组求解：
其中，r1为钢轨纵向电阻，r2为排流网纵向电阻，r3为埋地管线纵向电阻，R12为钢轨对排流网过渡电阻，R23为排流网对埋地管线过渡电阻，R34为埋地管线对大地过渡电阻。
进一步需对地铁机车牵引曲线和地铁机车位置曲线进行离散化处理，假设地铁机车牵引电流曲线为f(t)，地铁机车位置曲线为G(t)，f(t)和G(t)均可通过地铁监测系统数据获得，地铁隧道中的固定参数为：钢轨纵向电阻、排流网纵向电阻、埋地管线纵向电阻、钢轨对排流网过渡电阻、排流网对管线过渡电阻，埋地管线内电流计算模型I3(·)为f(t)和G(t)的函数，I3(·)可通过电阻网络分布模型进行计算，t为一个区间周期内的地铁机车运行时间，沿管线长度方向上的杂散电流分布可以通过电阻网络模型获得：{I3(f(t),G(t)，S1’),I3(f(t),G(t)，S2’),…,I3(f(t),G(t)，Si’),…,I3(f(t),G(t)，Sn’)}＝I3(f(t),G(t)，S)，S’表示管线长度方向上的绝对任意位置，S1’,S2’,…,Si’,…,Sn’分别表示在隧道坐标系上的绝对坐标；单位腐蚀时间Δt与一个区间周期内机车运行时间t的关系为：nΔt＝t，n为离散化Δt的个数。
离散化位腐蚀时间Δt的个数n计算方法流程如图3所示，具体实施方式如下：
初始设定参数：p(收敛精度,设定为0.01)
步骤1：首先根据当前迭代次数ni(ni＝i)，进行如下计算：
步骤2：将ni+1进行累加：ni+1＝ni+1进行迭代计算
步骤3：利用计算L(ni+1):
步骤4：重复步骤1至步骤3，将计算结果L(ni+1)与上一次的计算结果L(ni)做差，若满足L(ni+1)-L(ni)≤p，则判断迭代结束，ni+1则为所求离散化位腐蚀时间Δt的个数n；若L(ni+1)-L(ni)>p，则重复步骤2至步骤3，直至求出满足L(ni+1)-L(ni)≤p。
进一步，地铁隧道泄漏至埋地管线中的杂散电流I3(f(t),G(t),S)通过电阻网络模型计算获得。
在地铁隧道与埋地管线互相不为平行关系时，需要将隧道沿长度方向上的离散坐标转换至埋地管线沿长度方向上的离散坐标，坐标转换方法如下：
以地铁隧道轴向方向为原始坐标轴x轴，坐标轴水平面为x-y面，如图4所示，坐标轴垂直面为x-z面，如图5所示，地铁隧道轴向上离散坐标为C1,C2,…,Ci,…,Cm，C1为隧道任意牵引区间范围内起始坐标，C2位隧道长度方向上的第二个坐标，Ci为隧道范围内任意一点离散坐标，Cn为隧道内任意牵引区间范围内的终点坐标，m为离散坐标数量；
假设埋地管线与地铁隧道x轴在x-y面上投影夹角为α，埋地管线与地铁隧道x轴在x-z面上投影夹角为β，那么，管线在以自身为局部坐标系的相对坐标{S1,S2,…,Si,…,Sm}，管线在隧道坐标系上的绝对坐标{S1’,S2’,…,Si’,…,Sm’}表示为：
{S1’,S2’,…,Si’,…,Sm’}＝{S1cosαcosβ,S2cosαcosβ,…,Sicosαcosβ,…,Smcosαcosβ}；
S1,S2,…,Si,…,Sm即为坐标转换后的埋地管线沿长度方向上离散坐标，S1为区间内所覆盖埋地管线的相对坐标起点，Si为区间内所覆盖埋地管线范围内任意一相对坐标点，Sn为区间内所覆盖埋地管线的相对坐标终点。S1’,S2’,…,Si’,…,Sm’为管线上沿长度方向上的绝对离散坐标，S1为区间内所覆盖埋地管线的绝对坐标起点，Si为区间内所覆盖埋地管线范围内任意一绝对坐标点，Sm为区间内所覆盖埋地管线的绝对坐标终点。
最后通过埋地管线周围土壤的杂散电流腐蚀实验和自然腐蚀实验，实验系统如图6和图7所示，两个实验中所用工作电极均来自于预测目标管线同种钢材，试验箱均为亚克力板材料，杂散电流腐蚀实验中将管线钢工作电极和石墨电极埋置于土壤中，埋置深度均为150mm。管线钢工作电极通过加持装置只保留一个平面进行杂散电流腐蚀，通过可控直流电源通电实现杂散电流腐蚀过程，管线钢工作电极与石墨电极相距200mm且相对布置，工作电极的尺寸为10mm×10mm×5mm，石墨电极的尺寸为10mm×10mm×2mm，通电电流为is；自然腐蚀实验中将管线钢工作电极埋置于土壤中，埋置深度为150mm，管线钢工作电极的尺寸为10mm×10mm×5mm，模拟管线的自然腐蚀过程。实验中所用土壤需使用细筛除去大尺寸的石子。实验过程中需在实验箱上方加盖以保证土壤湿度稳定。
根据杂散电流腐蚀实验和自然腐蚀实验，能够通过三维数字显微镜在通电电流大小为is作用的情况下确定时间Δt内的管线钢工作电极表面平均腐蚀深度ΔH1和自然腐蚀作用下的平均腐蚀深度ΔH2，根据Kn＝ΔH2/ΔH1确定自然腐蚀因子Kn，通过拟合获得通电电流is与Kn的关系Kn(is)，通过拟合通电电流is与ΔH1/Δt的关系获得通电电流与腐蚀速率的关系C(is)。其次，计算地铁机车运行时间t至t+△t离散时间杂散电流作用下累计缺陷深度在管线上的分布{Ha(S1’),Ha(S2’),…,Ha(Si’),…,Ha(Sm’)}，计算方法为：
Ha(S'i)＝C(Is(f(t+Δt),G(t+Δt),Si’))·Δt
最终地铁机车运行时间t至t+△t离散时间内考虑自然腐蚀作用与杂散电流腐蚀作用的累积腐蚀深度Hd为累积腐蚀缺陷在管线上分布的最小值，即为：
其中，min{·}为取最小值的函数。
从t＝0时开始计算，重复上述过程，直至其中j为满足/>时累积腐蚀深度累加的次数，Hi为预测时刻的管线当前腐蚀深度，则m×△t为所预测的管线剩余寿命。

Claims (4)

1.一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法，其特征在于：利用预测系统进行埋地管线剩余寿命预测，所述的预测系统包括杂散电流干扰分布计算模块、坐标转换模块、缺陷极限尺寸计算模块和累计腐蚀深度实验及分布计算模块，其中：
杂散电流干扰分布计算模块：用以计算在地铁机车牵引电流和位置不断变化下的动态干扰，对牵引电流曲线进行离散化处理以近似计算单位腐蚀时间Δt下对埋地管的腐蚀深度；用于获得埋地管线在单位腐蚀时间Δt内恒定杂散电流流入所造成的电位偏移分布和因杂散电流所造成土壤电位的偏移；
坐标转换模块：当地铁隧道与埋地管线互相不为平行关系时，用以将隧道沿长度方向上的离散坐标转换至埋地管线沿长度方向上的离散坐标；
缺陷极限尺寸计算模块：通过管线内压、管线直径、管线壁厚、体积扩张系数、管线流动压力、管线的屈服极限和强度极限计算埋地管线极限壁厚Hl；
累计腐蚀深度实验及分布计算模块：用以计算杂散电流作用下累计腐蚀缺陷深度在埋地管线上的分布，在△t离散时间内考虑自然腐蚀作用与杂散电流腐蚀作用的累积腐蚀深度Hd为累积腐蚀缺陷在管线上分布的最小值通过比对管线极限壁厚Hl即可预测管线剩余寿命；
首先利用杂散电流干扰分布计算模块计算在地铁机车牵引电流和位置不断变化下的动态干扰，对牵引电流曲线进行离散化处理以近似计算单位腐蚀时间Δt下对埋地管的腐蚀深度，从而获得埋地管线在单位腐蚀时间Δt内恒定杂散电流流入所造成的电位偏移分布和因杂散电流所造成土壤电位的偏移；然后利用坐标转换模块将隧道沿长度方向上的离散坐标转换至埋地管线沿长度方向上的离散坐标；通过缺陷极限尺寸计算模块计算埋地管线的极限壁厚；最后利用累计腐蚀深度实验及分布计算模块计算杂散电流动态干扰下累计腐蚀缺陷深度在埋地管线上的分布，在离散时间△t内考虑自然腐蚀作用与杂散电流腐蚀作用的累积腐蚀深度Hd为累积腐蚀缺陷在管线上分布的最小值通过比对当前检测到的腐蚀深度即可预测的管线剩余寿命；
离散化处理具体为对地铁机车牵引曲线f(t)和地铁机车位置曲线G(t)进行离散化处理：首先通过地铁监测系统数据获得地铁机车牵引电流曲线为f(t)，地铁机车位置曲线为G(t)，获得地铁隧道中的固定参数为：钢轨纵向电阻、排流网纵向电阻、埋地管线纵向电阻、钢轨对排流网过渡电阻、排流网对管线过渡电阻，其中通过电阻网络分布模型计算得到埋地管线内电流计算模型I3(·)，埋地管线内电流计算模型I3(·)为f(t)和G(t)的函数，获取一个区间周期内的地铁机车运行时间t，利用电阻网络模型：{I3(f(t),G(t)，S1’),I3(f(t),G(t)，S2’),…,I3(f(t),G(t)，Si’),…,I3(f(t),G(t)，Sn’)}＝I3(f(t),G(t)，S)，获得沿埋地管线长度方向上的杂散电流分布，S’表示埋地管线长度方向上的绝对任意位置，S1’,S2’,…,Si’,…,Sn’分别表示在隧道坐标系上的绝对坐标；单位腐蚀时间Δt与一个区间周期内地铁机车运行时间t的关系为：nΔt＝t，n为将地铁机车运行时间t离散为Δt的个数；
离散化腐蚀时间Δt的个数n的具体计算方法为：
初始设定参数：将收敛精度p定为0.01；
步骤1：首先根据当前迭代次数ni(ni＝i)，进行如下计算：
步骤2：将ni+1进行累加：ni+1＝ni+1进行迭代计算
步骤3：利用计算L(ni+1):
步骤4：重复步骤1至步骤3，将计算结果L(ni+1)与上一次的计算结果L(ni)做差，若满足L(ni+1)-L(ni)≤p，则判断迭代结束，ni+1则为所求离散化位腐蚀时间Δt的个数n；若L(ni+1)-L(ni)>p，则重复步骤2至步骤3，直至求出满足L(ni+1)-L(ni)≤p，t为地铁机车运行时间,f(t)为地铁机车牵引电流曲线；
累计腐蚀深度实验及分布计算具体为：
首先通过埋地管线周围土壤的杂散电流腐蚀实验和自然腐蚀实验，两个实验中所用工作电极均与预测目标管线使用同种钢材，杂散电流腐蚀实验中将管线钢工作电极和石墨电极埋置于土壤中，埋置深度均为150mm，通过可控直流电源通电实现杂散电流腐蚀过程，其中管线钢工作电极与石墨电极相距200mm且相对布置，工作电极的尺寸为10mm×10mm×5mm，石墨电极的尺寸为10mm×10mm×2mm，通电电流为is；自然腐蚀实验中将管线钢工作电极埋置于土壤中，埋置深度为150mm，管线钢工作电极的尺寸为10mm×10mm×5mm，模拟管线的自然腐蚀过程；根据杂散电流腐蚀实验和自然腐蚀实验，通过三维数字显微镜在通电电流大小为is作用的情况下确定离散时间Δt内的管线钢工作电极表面平均腐蚀深度ΔH1和自然腐蚀作用下的平均腐蚀深度ΔH2，根据Kn＝ΔH2/ΔH1确定自然腐蚀因子Kn，通过拟合获得通电电流is与Kn的关系Kn(is)，通过拟合通电电流is与ΔH1/Δt获得通电电流与腐蚀速率的关系C(is)；其次，计算在地铁机车运行时间t至t+△t离散时间杂散电流作用下累计缺陷深度在管线上的分布{Ha(S1’),Ha(S2’),…,Ha(Si’),…,Ha(Sm’)}，计算方法为：
Ha(S′i)＝C(Is(f(t+△t),G(t+△t),Si’))·△t
最终地铁机车运行时间t至t+△t离散时间内考虑自然腐蚀作用与杂散电流腐蚀作用的累积腐蚀深度Hd为累积腐蚀缺陷在管线上分布的最小值，即为：
其中，min{·}为取最小值的函数；
重复上述过程，直至其中j为满足/> 时累积腐蚀深度累加的次数，Hi为预测时刻的管线当前腐蚀深度，则m×△t为所预测的管线剩余寿命。
2.根据权利要求1所述的受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法，其特征在于：地铁隧道泄漏至管线中的杂散电流I3(f(t),G(t),S)通过电阻网络模型计算获得。
3.根据权利要求1所述的受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法，其特征在于电阻网络模型为：
组建杂散电流泄漏过程的电阻网络模型，分别包含钢轨纵向导体层、排流网纵向导体层、埋地管线纵向导体层和接地层，四层导体组成了电阻网络，包括钢轨电流I1(S)、排流网电流I2(S)、埋地管线电流I3(S)、钢轨电位U1(S)、排流网电位U2(S)和埋地管线电位U3(S)，利用以下微分方程组求解：
其中，r1为钢轨纵向电阻，r2为排流网纵向电阻，r3为埋地管线纵向电阻，R12为钢轨对排流网过渡电阻，R23为排流网对埋地管线过渡电阻，R34为埋地管线对大地过渡电阻。
4.根据权利要求1所述的受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法，其特征在于，坐标转换方法如下：
以地铁隧道轴向方向为原始坐标轴x轴，坐标轴水平面为x-y面，坐标轴垂直面为x-z面，地铁隧道轴向上离散坐标为C1,C2,…,Ci,…,Cm，C1为隧道任意牵引区间范围内起始坐标，C2位隧道长度方向上的第二个坐标，Ci为隧道范围内任意一点离散坐标，Cm为隧道内任意牵引区间范围内的终点坐标，m为离散坐标数量；
设埋地管线与地铁隧道x轴在x-y面上投影夹角为α，埋地管线与地铁隧道x轴在x-z面上投影夹角为β，那么埋地管线在以自身为局部坐标系的相对坐标{S1,S2,…,Si,…,Sm}，管线在隧道坐标系上的绝对坐标{S1’,S2’,…,Si’,…,Sm’}表示为：{S1’,S2’,…,Si’,…,Sm’}＝{S1cosαcosβ,S2cosαcosβ,…,Sicosαcosβ,…,Smcosαcosβ}；
S1,S2,…,Si,…,Sm即为坐标转换后的埋地管线沿长度方向上离散坐标，S1为区间内所覆盖埋地管线的相对坐标起点，Si为区间内所覆盖埋地管线范围内任意一相对坐标点，Sn为区间内所覆盖埋地管线的相对坐标终点；S1’,S2’,…,Si’,…,Sm’为管线上沿长度方向上的绝对离散坐标，S1为区间内所覆盖埋地管线的绝对坐标起点，Si为区间内所覆盖埋地管线范围内任意一绝对坐标点，Sm为区间内所覆盖埋地管线的绝对坐标终点。
CN202110726421.2A 2021-06-29 2021-06-29 一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法 Active CN113609640B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110726421.2A CN113609640B (zh) 2021-06-29 2021-06-29 一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110726421.2A CN113609640B (zh) 2021-06-29 2021-06-29 一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法

Publications (2)

Publication Number Publication Date
CN113609640A CN113609640A (zh) 2021-11-05
CN113609640B true CN113609640B (zh) 2023-09-29

Family

ID=78303843

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110726421.2A Active CN113609640B (zh) 2021-06-29 2021-06-29 一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法

Country Status (1)

Country Link
CN (1) CN113609640B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115096806B (zh) * 2022-06-21 2024-04-12 徐州中矿传动轨道科技有限公司 隧道钢混结构杂散电流腐蚀高风险区域识别方法及系统

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101071098A (zh) * 2007-06-19 2007-11-14 广州市煤气公司 地下钢质燃气管网管道腐蚀预测系统
KR20090003903A (ko) * 2007-07-05 2009-01-12 한국가스공사 매설 배관의 외부부식 직접평가 방법
CN102590459A (zh) * 2012-03-09 2012-07-18 中国石油天然气股份有限公司 一种埋地管道腐蚀检测评价系统和方法
CN104122196A (zh) * 2014-07-15 2014-10-29 北京工业大学 应力条件下油气管道用高强钢直流杂散电流腐蚀评价方法
CN107991873A (zh) * 2017-12-11 2018-05-04 中国矿业大学 一种基于卷积神经网络及bp神经网络的地铁杂散电流泄漏等级预测方法
CN108827867A (zh) * 2018-04-23 2018-11-16 中国矿业大学 一种基于数据挖掘技术的地铁杂散电流腐蚀特征信息的提取方法
CN109374513A (zh) * 2018-09-29 2019-02-22 卢俊文 管道外防腐层检测分析系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101071098A (zh) * 2007-06-19 2007-11-14 广州市煤气公司 地下钢质燃气管网管道腐蚀预测系统
KR20090003903A (ko) * 2007-07-05 2009-01-12 한국가스공사 매설 배관의 외부부식 직접평가 방법
CN102590459A (zh) * 2012-03-09 2012-07-18 中国石油天然气股份有限公司 一种埋地管道腐蚀检测评价系统和方法
CN104122196A (zh) * 2014-07-15 2014-10-29 北京工业大学 应力条件下油气管道用高强钢直流杂散电流腐蚀评价方法
CN107991873A (zh) * 2017-12-11 2018-05-04 中国矿业大学 一种基于卷积神经网络及bp神经网络的地铁杂散电流泄漏等级预测方法
CN108827867A (zh) * 2018-04-23 2018-11-16 中国矿业大学 一种基于数据挖掘技术的地铁杂散电流腐蚀特征信息的提取方法
CN109374513A (zh) * 2018-09-29 2019-02-22 卢俊文 管道外防腐层检测分析系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
长输管道腐蚀缺陷检测技术与应用;周方舟等;《仪表电信》;第35卷(第3期);第72-74页 *

Also Published As

Publication number Publication date
CN113609640A (zh) 2021-11-05

Similar Documents

Publication Publication Date Title
CN104122196B (zh) 应力条件下油气管道用高强钢直流杂散电流腐蚀评价方法
Wang et al. Stray current distributing model in the subway system: A review and outlook
CN104674228B (zh) 一种定向钻穿越管道阴极保护的检测方法
CN108562782A (zh) 一种杂散电流的获取方法及杂散电流的检测系统
Shou et al. Numerical analysis of the mechanical behaviors of pressurized underground pipelines rehabilitated by cured-in-place-pipe method
CN113609640B (zh) 一种受地铁杂散电流腐蚀影响的埋地管线剩余寿命预测方法
CN113433057A (zh) 一种埋地金属管道杂散电流腐蚀速率预测装置及方法
CN106355320B (zh) 一种灰色系统海底油气管道坠物腐蚀可靠性评估方法
CN101696758B (zh) 非开挖施工大口径管道腐蚀控制方法
Charalambous et al. Influence of soil structures on corrosion performance of floating-DC transit systems
Wang et al. Remaining lifetime assessment of gas pipelines subjected to stray current interference using an integrated electric-electrochemical method
Zhichao et al. Evaluation of metro stray current corrosion based on finite element model
CN106777760A (zh) 一种基于数值分析的金属管道杂散电流预测方法
Li et al. Effect of soil salt content on stray current distribution in urban rail transit
Wang et al. A probabilistic-based model for dynamic predicting pitting corrosion rate of pipeline under stray current interference
CN112287589B (zh) 一种埋地钢质管道防腐层破损点定位有限元方法
Bortels et al. Three-dimensional boundary element method and finite element method simulations applied to stray current interference problems. A unique coupling mechanism that takes the best of both methods
Xi et al. Reliability analysis of buried polyethylene pipeline subject to traffic loads.
Feng et al. Influence of urban rail transit on corrosion of buried steel GAS pipeline
Bongiorno et al. Variability of measured railway track conductance due to test setup
CN201237085Y (zh) 一种用于燃气管道风险管理的装置
Wang et al. Assessment of interference scope induced by stray current in the metro depot for the corrosion risk of buried metal pipeline
Krissa et al. CP Management of Multiple Pipeline Right-Of-Ways
CN112213564B (zh) 一种铁路土壤电阻率测量方法
Desjardins Corrosion rate and severity results from in-line inspection data

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