Characteristics Analysis of Ground Penetrating Radar Signals for Groundwater Pipe Leakage Environment
-
摘要: 探地雷达(GPR)在地下水管渗漏的检测中具有良好的应用前景。前期研究表明:地下水管的渗漏会在雷达剖面中形成震荡信号,但其形成机理尚不明晰。为揭示不同材质地下管道渗漏后的探地雷达信号特征的成因,该文结合物理模型试验与渗流场-电磁场数值模拟分析干砂中PVC管和金属管渗漏前后雷达信号特征、渗漏后震荡信号的形成机理和电磁波传播路径。结果表明:地下水管发生渗漏后,管道周围区域出现一定分层状态,电磁波在传播过程中存在更多界面反射和界面间的多次波。PVC管渗漏后,管道顶部、底部反射信号和爬行波信号在渗漏区中多次反射形成复杂的震荡双曲线信号,而金属管渗漏后管壁与渗漏区间存在多次反射。研究成果可为探地雷达在地下水管渗漏探测实际应用提供技术支持。Abstract: Ground Penetrating Radar (GPR) has a good prospect in the detection of underground water pipeline leakage. Previous studies have shown that the underground water pipes leakage can form oscillating hyperbolic signals in GPR profile, but its formation mechanism is not clear. To reveal the formation mechanism of GPR signal characteristics after underground pipeline leakage of different materials, physical model test and multi-physics numerical simulation are carried out to analyzes the characteristics of GPR signals before and after leakage, the formation mechanism of oscillating hyperbolic signal of PVC/metal pipe after leakage and the propagation path of electromagnetic wave. The results show that a stratified leakage zone appears around the underground water pipe after leakage, and more interface reflections and multiple waves bouncing between interfaces occur during the propagation of electromagnetic waves. After the leakage of PVC pipe, creeping wave signals and the reflected signals at the top and bottom of the PVC pipe are reflected many times in the leakage zone to form complex oscillatory hyperbolic signals, while there are multiple reflections between the pipe outer wall and the leakage zone after the metal pipe leakage. The results are helpful to promote the practical applications of GPR to underground pipeline leakage detection and related data interpretation.
-
Key words:
- Ground Penetrating Radar (GPR) /
- Pipe detection /
- Leakage localization
-
1. 引言
紫外(UltraViolet, UV)波长在200~280 nm的波段,由于臭氧分子的强烈吸收作用,使得该波段的紫外光在近地面几乎衰减为零,称该波段为“日盲区”[1]。无线紫外光通信主要利用该波段的紫外光作为载体,通过大气分子、气溶胶等微粒的散射作用进行信息的传递[2],具有非直视(Non-Line-Of-Sight, NLOS)、高保密性、低窃听率、抗干扰能力强等优点,特别适用于无人机编队飞行、海军舰艇编队、陆军装甲编队、导弹车车队在无线电静默和复杂环境下的军事隐秘通信[3,4]。但由于紫外光的散射通信以及大气湍流的偏折作用,大气信道衰减严重、接收信号微弱、具有高路径损耗等,使得紫外光通信传输距离有限[5]。因此,通过多跳中继方式来延长无线紫外光通信距离成为研究的热点。
目前已有学者对紫外光多跳中继通信进行了相关研究,文献[6]在紫外光通信网中通过多跳中继方式增加了通信范围,并节省了功率消耗。文献[7]研究了紫外光多跳中继通信中节点连通问题,适当调整节点密度、发射功率和数据速率,使得紫外光多跳中继通信网络中中继孤立节点存在的概率趋近于零。文献[8]针对紫外光通信中光源器件发射功率受限和大气信道衰减严重造成传输距离近的问题,提出采用多跳中继方式构建紫外光长距离通信链路,提高了系统功率利用率。文献[9]研究结果表明,选择合适的无线紫外光通信系统配置对提高多跳中继通信系统性能是至关重要的。
上述研究都是建立在节点等距离分布下的性能研究,而对于紫外光多跳中继通信系统中节点随机分布的最优跳数问题研究较少。文献[10]为减少无线传感器网络节点的能量消耗,使用紫外光作为信息载体,在等距离分布下得出了使能量消耗最低的最优跳数表达式。但研究未考虑紫外光通信系统配置问题,而紫外光通信的整体性能高度依赖于系统配置,特别是收发仰角[9]。因此,本文基于无线紫外光非直视单次散射模型,研究了车队中车辆在随机分布状态下无线紫外光多跳中继通信的收发仰角与最优跳数的关系,根据信道容量和路径损耗得出两者之间的近似表达式,并分析了不同收发仰角下的系统性能。
2. 无线紫外光多跳中继通信系统模型
2.1 无线紫外光非直视单次散射模型分析
由于大气对紫外光强散射作用,使紫外光通信信道能够构成非直视通信链路。实际中紫外光通信是多次散射的,但在短距离通信中,单次散射传输为主[11],因此本文以单次散射模型作为研究基础。紫外光单次散射模型如图1所示[12],βT是发射仰角,βR是接收仰角,θT是发散角,θR是接收视场角,V是有效散射体,θS是散射角,r是收发端基线距离。
发射端TX以βT和θT的角度向空间发射光信号,光信号在有效散射体V内散射后,接收端RX以βR和θR的角度进行光信号接收,则紫外光单次散射通信的接收光功率为[12]
Pr=PtArKsPsθRθ2Tsin(βT+βR)32π3rsinβT(1−cosθT2)⋅exp[−Ker(sinβT+sinβR)sin(βT+βR)] (1) 其中,Pt是发射光功率,Ar是接收孔径面积,Ks是大气散射系数,Ps是散射角θs的相函数,Ke是大气信道衰减系数,且Ke=Ka+Ks, Ka是大气吸收系数。
非直视紫外光散射通信的路径损耗可表示为发射光功率与接收光功率的比值,如式(2)所示
L=PtPr=32π3rsinβT(1−cosθT2)ArKsPsθRθ2Tsin(βT+βR)⋅exp[Ker(sinβT+sinβR)sin(βT+βR)] (2) 式(2)为路径损耗计算表达式,只适用于收发仰角较小的状态,当收发仰角较大时,不再适用于分析非直视紫外光散射通信。而在实际应用中当通信距离小于1 km时,通常采用路径损耗简化公式[13]
L=ξrα (3) 其中,
ξ 是路径损耗因子,α 是路径损耗指数,其都是与收发仰角有关的参数。2.2 无线紫外光多跳中继节点随机分布结构模型
本文研究的紫外光多跳通信结构模型,考虑1维N跳网络结构且每跳的距离是随机性的,其无线紫外光多跳中继节点随机分布模型如图2所示,当
r1=r2=···=rN=d/N 时,即为等距离分布。紫外光多跳中继通信系统有一个源节点S,一个目的节点D, N–1个中继节点Ri (i=1, 2, ···, N–1),且Ri随机分布在S到D的距离范围内。假设所有节点采用半双工通信方式,考虑经典N跳时分解码转发协议,每个中继节点Ri接收Ri–1发送的信息,并解码转发到Ri+1节点。假设中继节点位置是独立随机变量,并设Xi是第i个中继节点在S到D间的随机位置,则Xi的概率密度函数为[14]
f(Xi)={1rδ, r(i−δ2)≤Xi≤r(i+δ2)0,其他 (4) 其中,r=d/N是等分距离,d是S到D的实际距离,
δ∈[0,1] 是随机偏移范围,其表征中继节点位置的随机性或者不确定性。定义
ri=Xi−Xi−1(i=2,3,···,N−1) ,特别地,r1=X1 ,rN=d−XN−1 ,则ri的累积分布函数表示为对于i=1或N时,
F(ri)={0,ri≤r(1−δ2)1rδ(ri−r(1−δ2)), r(1−δ2)≤ri≤r(1+δ2)1,ri≥r(1+δ2) (5) 对于i=2, 3, ···, N–1时,
F(ri)={0,ri≤r(1−δ)12(1rδ)2[ri−r(1−δ)]2, r(1−δ)≤ri≤r12(1rδ)2[r(1+δ)−ri]2, r≤ri≤r(1+δ)1,ri≥r(1+δ) (6) 3. 无线紫外光多跳中继节点随机分布最优跳数分析
根据紫外光NLOS单次散射模型,其紫外光通信在量子极限条件下的信噪比为[12]
γSNR=ηfηrPtλ2hcBL=μξ−1r−α (7) 其中,
μ=ηfηrPtλ2hcB ,ηf 和ηr 分别表示滤光片透过率和光电倍增管(Photo Multiplier Tube, PMT)探测效率,λ 是紫外波长,h是普朗克常数,c是光速,B=Kec2π 是紫外光通信信道带宽。利用香农公式可以计算出紫外光NLOS通信的信道容量[12],则紫外光单跳通信的频谱效率为
ηsh=log2(1+μξ−1d−α) (8) 在源节点到目的节点等距离分布N跳通信传输时,假设在任何时间点上只有一个节点进行传输,即在接收端无干扰,且每个节点在1/N时间传输相同的信息量,则每跳频谱效率是单跳通信的N倍,其等距离分布的频谱效率为
ηeq=1Nlog2(1+μξ−1(dN)−α) (9) 从式(9)可以看出频谱效率与跳数和收发仰角有关,因此研究使用近似理想路径路由计算方法[15],通过最大化频谱效率得到紫外光多跳中继通信系统节点等距离分布下的最优跳数近似表达式为
Nop=argmaxηeq≈[(2ε−1γ)1/α]+ (10) 其中,
ε=α+w(−αe−α)ln2 是仅取决于路径损耗指数的常数,w(⋅) 是郎伯W函数的主分支[16],γ=μξ−1d−α 是单跳情况下的接收信噪比。从式(10)可以看出最优跳数受到收发仰角的影响,在基于等距离分布的分析下,对中继节点随机分布时最优跳数进行了分析。在节点随机分布时,系统性能将取决于N跳中最长的一跳[17],则在此情况下频谱效率度量将是最长一跳距离上的平均值,如式(11)所示
¯η=Ermax[1Nlog2(1+μξ−1r−αmax)] (11) 其中,
rmax=maxi=1,⋯,Nri ,E[⋅] 是期望算子。为了求解节点随机分布下最优跳数值,需要求得
¯η 的闭合表达式,由式(5)和式(6)的累积分布函数可得r≤rmax≤r(1+δ) (12) 对式(11)应用Jensen不等式
E[f(x)]≥f(E[x]) ,则式(11)变换如式(13)所示E[log2(1+μξ−1r−αmax)]≥log2(1+μξ−1(E[rmax])−α) (13) 根据式(13),联合式(11)和式(12)求解得到
¯η 的上下界1Nlog2(1+μξ−1(dN(1+δ))−α)≤¯η≤1Nlog2(1+μξ−1(dN)−α) (14) 根据等距离分布最优跳数近似值的分析过程,求解式(14)的
¯η 下界最大化,即节点随机分布的最优跳数近似求解式如式(15)所示N′op≈[(2ε−1γ′)1/α]+ (15) 其中,
γ′=μξ−1d−α(1+δ)−α 是最大距离单跳情况下节点随机分布的接收信噪比。在节点随机分布情况下,通过式(15)求解得到的最优跳数值是源节点到目的节点多跳中继通信中的最大值。当
δ=0 时,式(15)等于式(10),最优跳数值达到最小值,即紫外光通信链路的频谱效率有最大值。4. 仿真结果及分析
4.1 仿真实验参数设置
通常在车队中车辆与车辆之间的距离根据车速应保持在30~100 m,而整个车队的最大距离不超过1~3 km[18]。基于此参数规定和最优跳数的分析,对节点等距离分布和随机分布的两种情况进行了仿真分析。其仿真实验的主要参数见表1。
表 1 系统主要仿真参数参数 数值 紫外波长 260 nm PMT探测效率 0.3 滤光片透过率 0.6 吸收系数 0.802×10–3 m–1 米氏散射系数 0.284×10–3 m–1 瑞利散射系数 0.266×10–3 m–1 普朗克常数h 6.6×10–34 4.2 节点等距离分布最优跳数仿真分析
在车队中车辆等间隔匀速运动时,通常满足节点等距离分布状态,即
δ=0 。根据不同收发仰角下ξ 和α 的取值分析[14],在Pt=30 mW, d=500 m和等距离分布情况下,本文对βT<βR, βT=βR和βT>βR 3种不同状态的收发仰角的单跳通信、多跳通信和最优跳通信进行了仿真,在d=500 m的多跳通信中,设定最大通信跳数为9跳,即车队有10辆车辆,其他仿真参数如表1所示,分析了3种通信方式的性能。从图3、图4和图5可以看出,当发射功率大于18 mW时,针对不同的收发仰角求得的最优跳数值都要比单跳通信和9跳通信获得较好的通信传输能力。另外,对比分析了相同最优跳数值下的不同收发仰角的频谱效率。从图6(a)和图6(b)可以看出,在相同最优跳数值下,紫外光通信的收发仰角对通信传输能力有着较大的影响。并且在紫外光多跳中继通信系统中,当有相同跳数时,使用小发射仰角和大接收仰角,能够使紫外光多跳通信系统获得更好的通信性能。
基于等距离分布模型的分析结果,将本文等距离分布下的最大频谱效率计算方法与最优能量计算方法进行对比分析。在βT<βR的紫外光几何结构配置下,不同方法下最优跳数和通信传输能力的对比情况如图7和图8所示。从图7可以看出两种方法的效果能够基本达到一致。当发射仰角固定时,增大接收仰角角度,其最优跳数随着通信距离的增加,最大频谱效率计算方法比最优能量计算方法的跳数值多一个数量级。
图8是根据图7中d=600 m时的最优跳数的分界点值的性能对比。此时,最大频谱效率计算方法最优跳数为6,最优能量计算方法最优跳数为5。从图8可以看出,当采用小功率传输时,最大频谱效率计算方法的通信传输能力优于最优能量计算方法,同时也达到节约功率的需求。在等距离分布模型下,当源节点到目的节点的通信距离确定时,根据无线紫外光通信的收发仰角得到适当的中继数,能够使整体通信传输能力达到最大。
4.3 节点随机分布最优跳数仿真分析
当车队在行驶中根据路况和调配信息,随时调整车速,此过程导致车队处于非匀速状态,不再满足等距离分布。本节对节点随机分布下的最优跳数进行了仿真分析。主要仿真参数如表1所示,在发射功率Pt=30 mW, S到D的距离d=500 m下,对比分析了不同随机偏移范围下的不同收发仰角的频谱效率与跳数的关系。从图9(a)、图9(b)、图9(c)和图9(d)可以看出,在紫外光多跳通信系统中,随着
δ 的增加系统总体性能呈下降趋势。并且不同收发仰角在不同随机偏移范围下,都存在一个使总体系统的频谱效率下界达到最大值的最优跳数值。通过图9进一步可以看出在紫外光中继通信系统中,使用小发射仰角和大接收仰角系统配置,能够使紫外光多跳通信系统获得更好的传输能力。为进一步说明随机分布模型的实际应用,本文在Pt=30 mW和
δ=0.2 下,对比分析了βT<βR和βT>βR下的收发端距离变化的频谱效率与跳数的关系,如图10所示。从图10可以看出,频谱效率随着跳数的增加存在最大值,并且随着S到D距离的增加,当达到最优跳数时,频谱效率不再随着跳数的增加有较大的变化。进一步说明当紫外光长距离通信时,并不是跳数越多其通信性能就越好。5. 结束语
本文针对车队中车辆在等距离和随机分布状态下的最优跳数进行了研究。基于非直视紫外光单次散射模型、信道容量和路径损耗,依据使频谱效率最大化原则,得出计算收发仰角与频谱效率的近似关系表达式。通过分析等距离分布最优跳数,求得随机分布的最优跳数近似表达式。仿真结果表明,收发仰角影响着紫外光多跳中继通信系统的最优跳数。不同随机偏移范围和不同收发仰角都对应特定的最优跳数,与最优能量计算方法相比,最大频谱效率计算方法在发射功率小于25 mW时有更好的信息传输能力,并达到节约功率的需求。在紫外光长距离通信时,选取合适的中继数及小发射仰角和大接收仰角的几何结构配置,不仅能够提高车队间无线紫外光多跳中继通信系统的传输能力,也满足车队间保持稳定可靠的隐秘通信需求。
-
表 1 孔压和饱和度之间的关系
参数 数值 饱和度 0.2 0.3 0.45 0.7 0.9 1 孔压(kPa) –14 –7.4 –4.2 –2.5 –1.3 0 -
[1] 城乡建设统计年鉴[Z]. 住房与城乡建设部, 2020.Statistical Yearbook of Urban and Rural Construction[Z]. Ministry of Housing and Urban-Rural Development, 2020. [2] 王帅超. 城市地下管道渗漏引起的路面塌陷机理分析与研究[D]. [硕士论文], 郑州大学, 2017.WANG Shuaichao. Analysis and study on subsidence mechanisms of road caused by leakage of urban underground pipeline[D]. [Master dissertation], Zhengzhou University, 2017. [3] 城市供水统计年鉴[Z]. 统计局, 2018.Statistical Yearbook of Urban Water Supply[Z]. Statistics, 2018. [4] YU Yicheng, SAFARI A, NIU Xudong, et al. Acoustic and ultrasonic techniques for defect detection and condition monitoring in water and sewerage pipes: A review[J]. Applied Acoustics, 2021, 183: 108282. doi: 10.1016/j.apacoust.2021.108282 [5] ALMEIDA F, BRENNAN M, JOSEPH P, et al. On the acoustic filtering of the pipe and sensor in a buried plastic water pipe and its effect on leak detection: An experimental investigation[J]. Sensors, 2014, 14(3): 5595–5610. doi: 10.3390/s140305595 [6] LIU Zheng and KLEINER Y. State of the art review of inspection technologies for condition assessment of water pipes[J]. Measurement, 2013, 46(1): 1–15. doi: 10.1016/j.measurement.2012.05.032 [7] GAO Yan, BRENNAN M J, JOSEPH P F, et al. On the selection of acoustic/vibration sensors for leak detection in plastic water pipes[J]. Journal of Sound and Vibration, 2005, 283(3/5): 927–941. doi: 10.1016/j.jsv.2004.05.004 [8] DATTA S and SARKAR S. A review on different pipeline fault detection methods[J]. Journal of Loss Prevention in the Process Industries, 2016, 41: 97–106. doi: 10.1016/j.jlp.2016.03.010 [9] 廖光伟, 张春萍, 武治国, 等. 基于双参数监测的供水管网漏损实验研究[J]. 供水技术, 2018, 12(6): 17–21. doi: 10.3969/j.issn.1673-9353.2018.06.005LIAO Guangwei, ZHANG Chunping, WU Zhiguo, et al. Experimental study on leakage of water supply network based on the two-parameter monitoring[J]. Water Technology, 2018, 12(6): 17–21. doi: 10.3969/j.issn.1673-9353.2018.06.005 [10] BACH P M and KODIKARA J K. Reliability of infrared thermography in detecting leaks in buried water reticulation pipes[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2017, 10(9): 4210–4224. doi: 10.1109/jstars.2017.2708817 [11] LIU Hai, SHI Zhenshi, LI Jianhui, et al. Detection of road cavities in urban cities by 3D ground-penetrating radar[J]. Geophysics, 2021, 86(3): WA25–WA33. doi: 10.1190/geo2020-0384.1 [12] CATALDO A, DE BENEDETTO E, CANNAZZA G, et al. Enhancement of leak detection in pipelines through time-domain reflectometry/ground penetrating radar measurements[J]. IET Science, Measurement & Technology, 2017, 11(6): 696–702. doi: 10.1049/iet-smt.2016.0310 [13] CATALDO A, PERSICO R, LEUCCI G, et al. Time domain reflectometry, ground penetrating radar and electrical resistivity tomography: A comparative analysis of alternative approaches for leak detection in underground pipes[J]. NDT & E International, 2014, 62: 14–28. doi: 10.1016/j.ndteint.2013.10.007 [14] GAO Lei, SONG Hantao, LIU Hanlong, et al. Model test study on oil leakage and underground pipelines using ground penetrating radar[J]. Russian Journal of Nondestructive Testing, 2020, 56(5): 435–444. doi: 10.1134/s1061830920050058 [15] LAI W W L, CHANG R K W, SHAM J F C, et al. Perturbation mapping of water leak in buried water pipes via laboratory validation experiments with high-frequency ground penetrating radar (GPR)[J]. Tunnelling and Underground Space Technology, 2016, 52: 157–167. doi: 10.1016/j.tust.2015.10.017 [16] OCAÑA-LEVARIO S J, CARREÑO-ALVARADO E P, AYALA-CABRERA D, et al. GPR image analysis to locate water leaks from buried pipes by applying variance filters[J]. Journal of Applied Geophysics, 2018, 152: 236–247. doi: 10.1016/j.jappgeo.2018.03.025 [17] 胡群芳, 郑泽昊, 刘海, 等. 三维探地雷达在城市市政管线渗漏探测中的应用[J]. 同济大学学报:自然科学版, 2020, 48(7): 972–981. doi: 10.11908/j.issn.0253-374x.19395HU Qunfang, ZHENG Zehao, LIU Hai, et al. Application of 3D ground penetrating radar to leakage detection of urban underground pipes[J]. Journal of Tongji University:Natural Science, 2020, 48(7): 972–981. doi: 10.11908/j.issn.0253-374x.19395 [18] LAU P K W, CHEUNG B W Y, LAI W W L, et al. Characterizing pipe leakage with a combination of GPR wave velocity algorithms[J]. Tunnelling and Underground Space Technology, 2021, 109: 103740. doi: 10.1016/j.tust.2020.103740 [19] 沈宇鹏, 董淑海, 王卿, 等. 城市供水管道渗漏程度的渗流模型分析与探地雷达信号正演[J]. 工程地质学报, 2016, 24(S1): 422–429. doi: 10.13544/j.cnki.jeg.2016.s1.062SHEN Yupeng, DONG Shuhai, WANG Qing, et al. Different degrees of water pipelines leakage model ananysis and GPR signal forward simulation[J]. Journal of Engineering Geology, 2016, 24(S1): 422–429. doi: 10.13544/j.cnki.jeg.2016.s1.062 [20] AYALA-CABRERA D, HERRERA M, IZQUIERDO J, et al. GPR-based water leak models in water distribution systems[J]. Sensors, 2013, 13(12): 15912–15936. doi: 10.3390/s131215912 [21] 费康, 彭劼. ABAQUS岩土工程实例详解[M]. 北京: 人民邮电出版社, 2017: 177FEI Kang and PENG Jie. ABAQUS geotechnical Examples are Explained in Detail[M]. Beijing: People's Post and Telecommunications Publishing House, 2017: 177. [22] 董泽君, 鹿琪, 冯晅, 等. 探地雷达测量土壤含水量的应用研究[J]. 地球物理学进展, 2017, 32(5): 2207–2213. doi: 10.6038/pg20170549DONG Zejun, LU Qi, FENG Xuan, et al. Estimation of soil water content using ground penetrating radar[J]. Progress in Geophysics, 2017, 32(5): 2207–2213. doi: 10.6038/pg20170549 [23] MCCUTCHEON M C, FARAHANI H J, STEDNICK J D, et al. Effect of soil water on apparent soil electrical conductivity and texture relationships in a dryland field[J]. Biosystems Engineering, 2006, 94(1): 19–32. doi: 10.1016/j.biosystemseng.2006.01.002 [24] WARREN C, GIANNOPOULOS A, and GIANNAKIS I. gprMax: Open source software to simulate electromagnetic wave propagation for Ground Penetrating Radar[J]. Computer Physics Communications, 2016, 209: 163–170. doi: 10.1016/j.cpc.2016.08.020 [25] 刘素贞, 董硕, 张闯, 等. 基于电磁超声爬波法对空腔结构部位裂纹缺陷的检测[J]. 高电压技术, 2019, 45(7): 2119–2125. doi: 10.13336/j.1003-6520.hve.20190628008LIU Suzhen, DONG Shuo, ZHANG Chuang, et al. Detection of crack defects in cavity structures based on electromagnetic ultrasonic creeping wave method[J]. High Voltage Engineering, 2019, 45(7): 2119–2125. doi: 10.13336/j.1003-6520.hve.20190628008 [26] DEMIRCI S, YIGIT E, ESKIDEMIR I H, et al. Ground penetrating radar imaging of water leaks from buried pipes based on back-projection method[J]. NDT & E International, 2012, 47: 35–42. doi: 10.1016/j.ndteint.2011.12.008 [27] 王楠. 基于任意曲面建模技术的一致性几何绕射理论方法[D]. [博士论文], 西安电子科技大学, 2007.WANG Nan. Study on the uniform geometrical theory of diffraction method with targets constructed by arbitrary surfaces[D]. [Ph. D. dissertation], Xidian University, 2007. -