石油与天然气化工  2021, Vol. 50 Issue (6): 91-95
大埋深油气管道防腐层破损检测技术研究
冀洁梅1 , 何川2 , 王文明2     
1. 中国石油物资沧州有限公司北京分公司;
2. 中国石油大学(北京)机械与储运工程学院
摘要:为了解决大埋深油气管道由于表面的腐蚀和破损对管道安全、高效、持续的运输造成的不利影响,针对大埋深管道进行防腐层破损检测研究。利用电磁波衰减的理论研究,对大埋深管道检测防腐层破损处的变化规律进行了分析;采用Comsol软件建立了相应的三维埋地管道仿真模型,模拟管道防腐层破损检测研究;通过对仿真模型的分析,得出了在不同埋深情况下、不同防腐层破损面积、不同频率的对应结果;研究确定了埋地管道防腐层更容易出现破损的位置。
关键词大埋深管道    防腐层破损检测    Comsol软件    仿真模型    
Study on damage detection of anticorrosive coating for deep buried oil and gas pipeline
Ji Jiemei1 , He Chuan2 , Wang Wenming2     
1. Beijing Branch, China Petroleum Materials Cangzhou Company Limited, Beijing, China;
2. College of Mechanical and Transportation Engineering, China University of Petroleum(Beijing), Beijing, China
Abstract: In order to solve the adverse effects of surface corrosion and damage of oil and gas pipelines with large buried depths on safe, efficient and continuous transportation of pipelines, this paper conducts researches of damage detection of anti-corrosion coating for deep buried pipelines. Through the theoretical research and simulation of electromagnetic wave attenuation, the law of changes in the detection of damage of the anti-corrosion coating with the deep buried pipeline is analyzed, and the corresponding three-dimensional buried pipeline simulation model is established using Comsol software to simulate the research of the pipeline anti-corrosion damage detection. Through the analysis of the simulation model, the corresponding results of different anti-corrosion layer damaged areas and different frequencies under different buried depths are obtained, and the location where the buried pipeline anticorrosive coating is more likely to be damaged is determined through pipeline corrosion test.
Key words: deep buried pipeline    damage detection of anticorrosive coating    Comsol software    simulation model    

随着社会的快速发展,城市建设规模不断扩大,大埋深管道在现代化城市建设中发挥着重要作用[1]。大埋深管道是指埋深较大的管道。由于大埋深管道长时间埋于地下,在温度、土壤、水分等环境的影响下,管道的防腐层容易出现老化、损坏、脱落,甚至导致管道破裂、穿孔,从而引发管道的泄漏,在给国家和社会带来巨大经济财产损失的同时,也会造成环境污染和灾难事故的发生[2-5]

为避免由于防腐层损坏导致重大事故的发生,应采取有效的方式检测大埋深油气管道防腐层的破损情况。在检测防腐层情况的实际工况下,目前运用的常见手段有直流电位梯度测试、PCM法、Pearson法、CIPS法等,通过检测结果可以掌握防腐涂层的具体信息[6-9]。然而受复杂管道环境以及不断升级的管道技术的影响,我国对管道外检测提出了更高的要求,单一化的检测模式已无法满足管道发展的需求[10-13]

针对当前大埋深油气管道防腐层破损检测技术存在的弊端,利用电磁波衰减理论,对大埋深管道检测防腐层破损处的变化规律进行了分析,同时,采用Comsol仿真软件建立了相应的三维埋地管道仿真模型来模拟管道防腐层破损检测研究。

1 大埋深管道检测防腐层破损处变化规律

为了研究大埋深管道检测防腐层破损处的变化规律,采用电磁波衰减理论进行分析,无界、无源的导电介质的麦克斯韦方程组表达见式(1)。

$ \left\{\begin{array}{l} \nabla \times H=\sigma E+j \omega \varepsilon E \\ \nabla \times E=-j \omega \mu H \\ \nabla \cdot H=0 \\ \nabla \cdot E=0 \end{array}\right. $ (1)

由(1)得:

$ \nabla \times H=j \omega\left(\varepsilon-j \frac{\sigma}{\omega}\right) E=j \omega \varepsilon_{c} E $ (2)

式中:▽为哈密顿算子(当▽后面接“×”时,▽为旋度; 当▽后面接“·”时,▽为散度; 当▽后面既不接“×”也不接“·”时,▽为梯度); H为磁场强度,A/m; E为电场强度, V/m; ε为相对介电常数,为F/m; μ为相对磁导率,H/m; ω为角速度, rad/s; j为虚数; σ为衰减系数,Np/m; εc为可导电媒质的复数介电常数,F/m,${\mathit{\varepsilon }_\mathit{c}} = \mathit{\varepsilon } - \mathit{j}\frac{\mathit{\sigma }}{\mathit{\omega }} = \mathit{\varepsilon }(1 - \mathit{j}\frac{\mathit{\sigma }}{{\mathit{\omega \varepsilon }}})$

通过对上述公式的推导和转化,得出的波动方程见式(3)和式(4)。

$ \left\{\begin{array}{l} \nabla^{2} E+\gamma^{2} E=0 \\ \nabla^{2} H+\gamma^{2} H=0 \end{array}\right. $ (3)
$ \gamma^{2}=\omega^{2} \mu \varepsilon_{c} $ (4)

假设在直角坐标系,一个沿着z轴正方向传播的均匀的电场强度大小只有x的一个分量Ex的平面电磁波,它的波动方程的一个解是:

$ E=e_{x} E_{0} e^{-j r z} $ (5)

式中: ex为沿x轴正向的单位矢量; E0为电场在地表时的初始振幅, V/m; z为笛卡尔坐标系的z轴坐标, m。

γ=β-jα,则:

$ E=e_{x} E_{0} e^{-j(\beta-j \alpha) z}=e_{x} E_{0} e^{-\alpha z} e^{-j \beta z} $ (6)

式中:α为电磁波衰减常数; β为相位常数; γ为传播常数。

由上式可得电场强度E随着因子z的增大而减小,所以电场强度的瞬时表达式为:

$ E(z, t)=e_{x} E_{\mathrm{m}} e^{-\alpha z} \cos \left(\omega t-\beta z+\varphi_{0}\right) $ (7)

式中: Em为电场强度的振幅值; φ0代表电场强度的初相角,t为时间,即:

$ E_{0}=E_{m} e^{-j \varphi_{0}} $ (8)

因为:

$ \gamma^{2}=\omega^{2} \mu \varepsilon_{c} $ (9)

所以:

$ (\beta-j \alpha)^{2}=\omega^{2} \mu\left(\varepsilon-j \frac{\sigma}{\omega}\right) $ (10)

故有:

$ \beta^{2}-\alpha^{2}-j 2 \alpha \beta=\omega^{2} \mu \varepsilon-j \omega \mu \sigma $ (11)

从而可以推出:

$ \beta^{2}-\alpha^{2}=\omega^{2} \mu \varepsilon $ (12)
$ 2 \alpha \beta=\omega \mu \sigma $ (13)

由以上的方程可以解得:

$ \left\{\begin{array}{l} \alpha=\omega \sqrt{\frac{\mu \varepsilon\left[\sqrt{1+\left(1+\frac{\sigma}{\omega \varepsilon}\right)^{2}}-1\right]}{2}} \\ \beta=\omega \sqrt{\frac{\mu \varepsilon\left[\sqrt{1+\left(1+\frac{\sigma}{\omega \varepsilon}\right)^{2}}+1\right]}{2}} \end{array}\right. $ (14)

时谐平面电磁波情况下,由麦克斯韦推出电场强度波动方程:

$ \nabla^{2} E+k_{c}^{2} E=0 $ (15)
$ k_{c}=\sqrt{\omega^{2} \mu \varepsilon-j \omega \mu \sigma}=\beta-j \alpha $ (16)

式中:kc为复波数,在垂直地面入射的平面电磁波条件下,可以由波动方程解出的电场强度的表达式如下:

$ E(z, t)=e_{x} E_{0} e^{-j(\beta-j \alpha)} $ (17)

电磁波在有损介质中传播会逐渐衰减,电场与磁场的最大幅度会随传播距离的增大按照e-αz的指数规律慢慢减小,衰减的快慢速度取决于衰减率α

$ \alpha=\omega \sqrt{\frac{\mu \varepsilon\left[\sqrt{1+\left(1+\frac{\sigma}{\omega \varepsilon}\right)^{2}}-1\right]}{2}} $ (18)

对衰减系数公式进行取极限:

$ \begin{array}{c} \lim \limits_{\omega \rightarrow \infty} \omega \sqrt{\frac{\mu \varepsilon\left[\sqrt{1+\left(1+\frac{\sigma}{\omega \varepsilon}\right)^{2}}-1\right]}{2}}=\\ \sqrt{\frac{\mu \varepsilon}{2}} \lim\limits _{\omega \rightarrow \infty} \sqrt{\frac{\left[\sqrt{1+\left(1+\frac{\sigma}{\omega \varepsilon}\right)^{2}}-1\right]}{\left(\frac{\sigma}{\omega \varepsilon}\right)^{2}}} \end{array} $ (19)

$\mathit{t} = {\left( {\frac{\mathit{\sigma }}{{\mathit{\omega \varepsilon }}}} \right)^2}$,则:

$ \text { 上式 }=\sigma \sqrt{\frac{\mu}{2 \varepsilon}} \lim \limits_{t \rightarrow 0} \sqrt{\frac{[\sqrt{1+t}-1]}{t}} $ (20)

继续通过等价无穷小替换,令$\mathop {{\rm{lim}}}\limits_{\mathit{t} \to 0} \sqrt {1 + \mathit{t}} - 1 = \frac{1}{2}\mathit{t}$,可得:

$ \sigma \sqrt{\frac{\mu}{2 \varepsilon}} \lim \limits_{t \rightarrow 0} \sqrt{\frac{\sqrt{1+t}-1}{t}}=\sigma \sqrt{\frac{\mu}{2 \varepsilon}} \lim \limits_{t \rightarrow 0} \sqrt{\frac{\frac{t}{2}}{t}}=\frac{\sigma}{2} \sqrt{\frac{\mu}{\varepsilon}} $ (21)

所以,得到电磁波衰减系数的极限为$\frac{\mathit{\sigma }}{2}\sqrt {\frac{\mathit{\mu }}{\mathit{\varepsilon }}} $

综上可知,电磁波衰减系数的极限是一个只与电导率、磁导率、介电常数有关的正数。其在传播过程中呈指数形式衰减。通常情况下,电磁波衰减系数会随频率的增大而增大,但当频率达到了一定数值以后,电磁波衰减系数基本上不会再随频率的增大而一直增加; 当频率增大到某个极值后,衰减系数会急剧下降。

2 仿真分析

通过Comsol软件对管道进行分析,埋地管道建模基本参数见表 1,模型材料属性见表 2。管道内的介质是空气,土壤模型采用半径为15 m、高10 m的圆柱体, 位于正中间的防腐层破损处俯视图为20 mm×20 mm正方形,用土壤材料填充破损处。

表 1    埋地管道建模基本参数

表 2    模型材料属性定义参数

为了分析土壤层表面的电流密度,在平行于管道轴向方向上设置一条路径,且该直线也正好处于管道防腐层破损处的正上方。通过检测该轴线上的电流密度大小来研究大埋深管道的具体变化规律。

当管道防腐层完好时,设置该路径与管道的垂直距离分别为0.2 m、3 m、6 m、10 m,对应的电流密度曲线见图 1。由图 1可以看出,当管道防腐层完好时,土壤层距离管道0.2 m、3 m、6 m平行管道轴线土壤层上直线的电流密度大小基本一致,但是随着土壤层上设置路径和管道轴线距离的增大,路径上的电流密度数值逐渐减小,而且深度越大,其数值越小,检测难度就会增加,对相应检测设备的要求也会越来越高。

图 1     防腐层完好时不同埋深的电流密度曲线图

当防护层破损时,不同埋深的电流密度曲线见图 2。从图 2可看出,当管道防腐层存在破损时(破损点在管道正中间,大小为20 mm×20 mm),土壤层距离管道0.2 m、3 m、6 m、10 m设置路径上的电流密度会先缓慢增大,在管道正中心达到最大值后再慢慢变小,可以根据峰值的大小和其突变判断防腐层破损点的位置; 在竖直方向上,随着距离的增大,峰值也在变小,即管道埋地越深,破损处产生的电流到地面的信号强度越微弱,检测难度越大。电流密度最大时(即在5 000 mm处)为防腐层破损点的位置,在管道埋地0.2 m处较明显。

图 2     防腐层破损时不同埋深的电流密度曲线图

在管道服役过程中,由于人为或环境因素,管道可能会在不同位置造成不同程度的损伤。为了研究上述情况下电流的传输方式,将防腐层破损点设置在管道中心5 m处,且破损区域边长分别为10 mm、15 mm、20 mm、25 mm、30 mm的正方形破损。埋深10 m时不同破损点的电流密度见图 3。由图 3可知,当管道防腐层存在破损点时,管道上的电流从管道防腐层的破损点处流进土壤,以破损点为球心往四面八方传播。随着管道防腐层破损面积的增大,流进土壤的电流也就越大,土壤层上路径的曲线峰值也就越高,破损点的正上方电流密度要比两边的电流密度大,曲线趋势呈一个反U字型; 边长为30 mm的正方形的破损区域处电流密度最大,由此方法可检测埋地管道防腐层的破损程度。

图 3     埋深10m时不同破损点的电流密度曲线图

除此之外,通过改变激励电流的频率也可以实现对大埋深管道防腐层破损时的检测研究,控制频率为8 Hz、32 Hz、64 Hz、128 Hz、256 Hz,检测电流密度情况,见图 4。由图 4可知,当管道埋深和防腐层上破损点相同时,施加相同大小的电流信号源,激励频率越小,土壤层上路径的电流密度越大。管道上泄漏的电流在低频情况下穿透土壤的能力越强,其峰值也随着频率的减少而增大。

图 4     埋深10 m不同频率的直线电流密度曲线图

综合以上研究,对大埋深油气管道防腐层破损规律进行了进一步分析。如图 5所示,通过钟表切分法将管道切分成多份,然后从中选取几个钟点方向的个点作为分析点,进而更清晰地分析管道防腐层破坏规律。以管道横切面圆心位置为原点,顺管道一周将管道切分为8块,从中选取5个点作为分析点,分别为A1(4点钟方向)、A2(3~9点钟方向)、A3(12点钟方向)、A4(10点钟方向)、A5(3点钟方向)。然后将带有破损的管道展开,如图 6所示。A2位置处的防腐层破损较为严重,A4处的防腐层破损较弱,管道A1~A5位置的防腐层破损程度如表 3所示。由此可得,在实际工况中,深埋管道的A1和A2位置较薄弱,更容易出现破损,应加强该处的防护。该方法能够较为便捷地发现油气管道防腐层相对薄弱的位置,对于进一步研究大埋深油气管道防腐层破损规律及解决措施提供了良好的参考。

图 5     管道位置切分示意图

图 6     展开管道破损位置分析图

表 3    埋地管道展开管道破损位置

3 结论

研究了大埋深管道的探测检测过程,通过利用电磁波衰减的理论研究,对大埋深管道检测防腐层破损处的变化规律进行了分析; 通过Comsol软件建立了相应的三维埋地管道仿真模型来模拟管道防腐层破损检测研究,得出以下结论:

(1) 当管道防腐层完好或者破损时,仿真结果所得曲线存在巨大的差异。当管道防腐层完好时,检测到土壤层轴线的电流密度和磁场强度大小基本一致,曲线图形是一条直线; 当管道防腐层破损以后,管道中的电流泄漏,在土壤层轴线上检测到的电流大小是一个倒U字型曲线,在破损位置的电流大小是最大的。

(2) 当大埋深管道防腐层存在不同破损点时,管道防腐层破损面积越大,流进土壤的电流也就越大,土壤层上路径的电流密度大小曲线峰值也就越高。

(3) 当其他基本条件都一样时,施加激励频率越小,对检测大埋深管道的破损情况越有利。

参考文献
[1]
赵洪涛. 城市地下管线探测技术方法及其应用[J]. 科技创新与应用, 2014(9): 286.
[2]
姚小静, 王威强, 张峰. 埋地长输管道防腐层综合检测技术[J]. 压力容器, 2006, 23(2): 49-51. DOI:10.3969/j.issn.1001-4837.2006.02.011
[3]
曹震峰, 张汉春, 葛如冰. 城市新型地下管线的探测方法及其应用[J]. 勘察科学技术, 2009(4): 58-60. DOI:10.3969/j.issn.1001-3946.2009.04.016
[4]
李学军, 洪立波. 城市地下管线探测与管理技术的发展及应用[J]. 城市勘测, 2010(4): 5-11. DOI:10.3969/j.issn.1672-8262.2010.04.001
[5]
刘吉东, 陈绍萍, 刘寒冰. 国内外油气管线腐蚀泄漏检测技术[J]. 油气田地面工程, 2003, 22(11): 48-49. DOI:10.3969/j.issn.1006-6896.2003.11.040
[6]
王春起. 大埋深管道位置探测技术探讨[J]. 城市勘测, 2009(增刊): 91-94.
[7]
章海生, 王春起, 黄智. 超常埋深燃气管道位置探测技术[J]. 煤气与热力, 2009, 29(10): B08-B10.
[8]
焦敬品, 何存富, 吴斌, 等. 管道超声导波检测技术研究进展[J]. 实验力学, 2002, 17(1): 1-9. DOI:10.3969/j.issn.1001-4888.2002.01.001
[9]
杨雄, 闫宏伟, 侯相荣, 等. 直轮式管道泄漏检测机器人弯管通过性能研究[J]. 石油机械, 2021, 49(4): 138-143.
[10]
杨琪. 埋地管道综合检测与评估技术研究及应用[D]. 北京: 北京化工大学, 2017.
[11]
赵书婷. 长输管道外检测技术现状及发展方向[J]. 工程技术, 2016(12): 305. DOI:10.3969/j.issn.2095-6630.2016.12.286
[12]
程浩. 长输天然气管道外检测综合技术研究与探讨[J]. 轻工科技, 2014, 30(11): 98-99.
[13]
王金龙, 何仁洋, 张海彬, 等. 海底管道检测最新技术及发展方向[J]. 石油机械, 2016, 44(10): 112-118.