一种用于准单站配置的步进频探地雷达建模方法

张经纬 刘新 叶盛波 方广有

引用本文:
Citation:

一种用于准单站配置的步进频探地雷达建模方法

    作者简介: 张经纬(1992–),男,湖北荆州人,中国科学院电子学研究所博士研究生,主要研究方向为探地雷达系统建模及探地雷达信号处理,目前已发表多篇期刊论文及国际会议论文;刘 新(1992–),男,山东新泰人,中国科学院电子学研究所博士研究生,主要研究方向为伪随机编码雷达系统设计及其在穿墙成像中的应用;叶盛波(1985–),男,湖北黄陂人,现为中国科学院电子学研究所副研究员,主要研究方向为超宽带探地雷达系统设计,目前发表多篇论文;方广有(1963–),男,河南人,现为中国科学院电子学研究所研究员,主要研究方向为超宽带雷达成像技术、微波成像新方法和技术、 超宽带电磁场理论及工程应用,目前发表多篇高质量论文.
    通讯作者: 方广有, gyfang@mail.ie.ac.cn
  • 基金项目:

    中国科学院面上基金(CXJJ-17-M140),国家重点研发计划(2017YFF0107700)

  • 中图分类号: TN958

A Modeling Method for Quasi-monostatic Stepped Frequency Ground Penetrating Radar

    Corresponding author: FANG Guangyou, gyfang@mail.ie.ac.cn
  • Fund Project: Chinese Academy of Sciences Innovation Fund Grant (CXJJ-17-M140), The National Key R&D Program of China (2017YFF0107700)

    CLC number: TN958

  • 摘要: 探地雷达是一种被广泛使用的无损检测技术。利用探地雷达对分层媒质进行全波反演时,构建精确的探地雷达正演模型具有十分重要的意义。该文提出一种可用于准单站配置的步进频探地雷达的建模方法。在该模型中,探地雷达系统及其与分层媒质间的相互作用被表示成线性方程,天线对雷达信号的影响被表示为只与频率有关的传输函数。为验证模型准确性,该文在实验室条件下搭建了准单站配置的步进频探地雷达系统,并对已知厚度的石膏板和木板的雷达测量信号进行全波反演。反演结果表明:石膏板和木板的厚度估计误差均不超过0.3 mm,验证了所提出的正演模型具有高准确度。利用石膏板和木板搭建分层模型,该文进一步比较了准单站配置和单站配置步进频探地雷达系统对介电常数差异较小的分层媒质的反演性能。实验结果表明:准单站配置探地雷达能获得更精确的反演参数。通过对分界面反射信号的信噪比估计可知,准单站配置比单站配置探地雷达系统能获得高出约10 dB的信噪比,因此具有更好的反演性能。
  • 图 1  准单站探地雷达建模框图

    Figure 1.  Modeling diagram for the quasi-monostatic ground penetrating radar

    图 2  准单站配置和单站配置探地雷达系统

    Figure 2.  Configurations of quasi-monostatic and monostatic GPR

    图 3  简化模型和原始模型建模误差

    Figure 3.  Fitting errors of the simplified model and original model

    图 4  利用厚度已知的石膏板验证准单站探地雷达模型的精度

    Figure 4.  Validating the accuracy of the quasi-monostatic model with the thickness-known plasterboard

    图 5  石膏板与木板的实测波形与仿真波形的对比

    Figure 5.  Comparison of the measured and simulated signal for plasterboard and woodblock

    图 6  分界面反射信号的计算模型

    Figure 6.  Models used to calculate the reflections on the interface

    图 7  单站配置探地雷达和准单站配置探地雷达的信噪比对比

    Figure 7.  The comparison of the SNR between the monostatic and quasi-monostatic GPR

    表 1  准单站探地雷达石膏板反演结果

    Table 1.  Inversion results of the plasterboard with the quasi-monostatic GPR

    高度(mm)反演相对介电常数反演厚度(mm)厚度估计误差(mm)
    534.32.3228.70.3
    522.52.3328.30.1
    515.32.3228.60.2
    下载: 导出CSV

    表 2  准单站探地雷达木板反演结果

    Table 2.  Inversion results of the woodblock with quasi-monostatic GPR

    高度(mm)反演相对介电常数反演厚度(mm)厚度估计误差(mm)
    531.52.7219.20.2
    523.22.7318.80.2
    512.42.7418.70.3
    下载: 导出CSV

    表 3  简化模型对石膏板的反演结果

    Table 3.  Inversion results of the plasterboard with simplified model

    高度(mm)反演相对介电常数反演厚度(mm)厚度估计误差(mm)
    533.82.4328.00.4
    522.02.4527.60.8
    514.72.4427.90.5
    下载: 导出CSV

    表 4  准单站配置雷达3层结构反演结果

    Table 4.  Inversion results of the three-layer structure with quasi-monostatic GPR

    天线高度
    (mm)
    反演石膏板
    介电常数
    石膏板反演厚度
    (mm)
    石膏板厚度误差
    (mm)
    反演木板
    介电常数
    木板反演厚度
    (mm)
    木板厚度误差
    (mm)
    532.22.3428.10.32.7119.50.5
    524.22.3228.80.42.7219.40.4
    510.32.3128.90.52.7518.50.5
    下载: 导出CSV

    表 5  单站配置雷达3层结构反演结果

    Table 5.  Inversion results of the three-layer structure with monostatic GPR

    天线高度
    (mm)
    反演石膏板
    介电常数
    石膏板反演厚度
    (mm)
    石膏板厚度误差
    (mm)
    反演木板
    介电常数
    木板反演厚度
    (mm)
    木板厚度误差
    (mm)
    535.32.3127.80.62.6320.11.1
    525.22.3028.90.52.7818.10.9
    513.32.3327.70.72.6519.80.8
    下载: 导出CSV
  • [1] 刘澜波, 钱荣毅. 探地雷达: 浅表地球物理科学技术中的重要工具[J]. 地球物理学报, 2015, 58(8): 2606–2617. doi: 10.6038/cjg20150802LIU Lan-bo and QIAN Rong-yi. Ground penetrating radar: A critical tool in near-surface geophysics[J]. Chinese Journal of Geophysics, 2015, 58(8): 2606–2617. doi: 10.6038/cjg20150802
    [2] ZHANG J W, YE S B, YI L, et al. A Hybrid method applied to improve the efficiency of full-waveform inversion for pavement characterization[J]. Sensors, 2018, 18(9): 2916. doi: 10.3390/s18092916
    [3] KLĘSK P, KAPRUZIAK M, and OLECH B. Statistical moments calculated via integral images in application to landmine detection from Ground Penetrating Radar 3D scans[J]. Pattern Analysis and Applications, 2018, 21(3): 671–684. doi: 10.1007/s10044-016-0592-5
    [4] DEIANA R, LEUCCI G, and MARTORANA R. New perspectives on geophysics for archaeology: A special issue[J]. Surveys in Geophysics, 2018, 39(6): 1035–1038. doi: 10.1007/s10712-018-9500-4
    [5] LAHOUAR S and AL-QADI I L. Automatic detection of multiple pavement layers from GPR data[J]. NDT & E International, 2008, 41(2): 69–81. doi: 10.1016/j.ndteint.2007.09.001
    [6] LAMBOT S, SLOB E C, VAN DEN BOSCH I, et al. Modeling of ground-penetrating radar for accurate characterization of subsurface electric properties[J]. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42(11): 2555–2568. doi: 10.1109/TGRS.2004.834800
    [7] WARREN C and GIANNOPOULOS A. Creating finite-difference time-domain models of commercial ground-penetrating radar antennas using Taguchi’s optimization method[J]. Geophysics, 2011, 76(2): G37–G47. doi: 10.1190/1.3548506
    [8] HYUN S Y, KIM S Y, and KIM Y S. An equivalent feed model for the FDTD analysis of antennas driven through a ground plane by coaxial lines[J]. IEEE Transactions on Antennas and Propagation, 2009, 57(1): 161–167. doi: 10.1109/TAP.2008.2009650
    [9] VENKATARAYALU N V, GAN Y B, LEE R, et al. Application of hybrid FETD-FDTD method in the modeling and analysis of antennas[J]. IEEE Transactions on Antennas and Propagation, 2008, 56(9): 3068–3072. doi: 10.1109/TAP.2008.928809
    [10] ALKHALIFEH K, HISLOP G, OZDEMIR N A, et al. Efficient MoM simulation of 3-D antennas in the vicinity of the ground[J]. IEEE Transactions on Antennas and Propagation, 2016, 64(12): 5335–5344. doi: 10.1109/TAP.2016.2618482
    [11] GENTILI G G and SPAGNOLINI U. Electromagnetic inversion in monostatic ground penetrating radar: TEM horn calibration and application[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(4): 1936–1946. doi: 10.1109/36.851775
    [12] LAMBOT S, SLOB E, and VEREECKEN H. Fast evaluation of zero-offset green’s function for layered media with application to ground-penetrating radar[J]. Geophysical Research Letters, 2007, 34(21): L21405. doi: 10.1029/2007GL031459
    [13] The data sheet of the portable VNA[OL]. http://www.c-oppermountaintech.com/products/12/S5048/.
  • [1] 蒋彦雯邓彬王宏强秦玉亮庄钊文 . 基于双频联合处理的太赫兹InISAR成像方法. 雷达学报, 2018, 7(1): 139-146. doi: 10.12000/JR17109
    [2] 覃尧黄春琳陆珉徐伟 . 基于小波变换与主成分分析的探地雷达自适应杂波抑制方法研究. 雷达学报, 2015, 4(4): 445-451. doi: 10.12000/JR15013
    [3] 王杰沈明威吴迪朱岱寅 . 基于主瓣杂波高效配准的机载非正侧视阵雷达STAP 算法研究. 雷达学报, 2014, 3(2): 235-240. doi: 10.3724/SP.J.1300.2014.13122
    [4] 刘向阳杨君刚孟进张晓牛德智 . 低信噪比下基于Hough变换的前视阵列SAR稀疏三维成像. 雷达学报, 2017, 6(3): 316-323. doi: 10.12000/JR17011
    [5] 魏民李小波王理 . 减小距离模糊影响的机载双基地雷达配置方法. 雷达学报, 2017, 6(1): 106-113. doi: 10.12000/JR16082
    [6] 吴秉横纪奕才方广有 . 带反射腔的分布式电阻加载探地雷达天线设计与分析. 雷达学报, 2015, 4(5): 538-544. doi: 10.12000/JR15070
    [7] 洪永彬张勇鲁振兴黄巍 . 一种高效的基于对比度的步进频雷达运动补偿算法. 雷达学报, 2016, 5(4): 378-388. doi: 10.12000/JR16068
    [8] 万显荣易建新程丰饶云华龚子平柯亨玉 . 单频网分布式外辐射源雷达技术. 雷达学报, 2014, 3(6): 623-631. doi: 10.12000/JR14156
    [9] 雷文太梁琼谭倩颖 . 基于自动反相校正和峰度值比较的探地雷达回波信号去噪方法. 雷达学报, 2018, 7(3): 294-302. doi: 10.12000/JR17113
    [10] 尹德叶盛波刘晋伟纪奕才刘小军方广有 . 一种用于高速公路探地雷达的新型时域超宽带TEM喇叭天线. 雷达学报, 2017, 6(6): 611-618. doi: 10.12000/JR17004
    [11] 易建新万显荣赵志欣程丰柯亨玉 . 单频网CP-OFDM 信号外辐射源雷达的分载波杂波抑制方法 (英文). 雷达学报, 2013, 2(1): 1-13. doi: 10.3724/SP.J.1300.2013.13030
    [12] 路东伟阙肖峰綦鑫聂在平 . 米波频段飞机目标全极化散射仿真与特性分析. 雷达学报, 2016, 5(2): 182-189. doi: 10.12000/JR16030
    [13] 周子铂蒋李兵王壮 . 一种基于波程差补偿的InISAR图像配准方法. 雷达学报, 2018, 7(6): 758-769. doi: 10.12000/JR18070
    [14] 折小强仇晓兰雷斌张薇卢晓军 . 一种极化熵结合混合GEV模型的全极化SAR潮间带区域地物分类方法. 雷达学报, 2017, 6(5): 554-563. doi: 10.12000/JR16149
    [15] 赵勇胜赵拥军赵闯 . 联合角度和时差的单站无源相干定位加权最小二乘算法. 雷达学报, 2016, 5(3): 302-311. doi: 10.12000/JR15133
    [16] 赵博黄磊周汉飞张亮李强黄敏 . 基于单频时变阈值的1-bit SAR成像方法研究. 雷达学报, 2018, 7(4): 446-454. doi: 10.12000/JR18036
    [17] 向寅张凯胡程 . 基于NUFFT的调频步进频高分辨成像与目标识别算法. 雷达学报, 2015, 4(6): 639-647. doi: 10.12000/JR15083
    [18] 谭云华王李波李廉林 . 一种抑制探地/穿墙成像多径虚假目标的新型概率模型:数值研究(英文). 雷达学报, 2015, 4(5): 509-517. doi: 10.12000/JR15066
    [19] 王福来庞晨李永祯王雪松 . 一种同时全极化雷达正交多相编码波形设计方法. 雷达学报, 2017, 6(4): 340-348. doi: 10.12000/JR16150
    [20] 高敬坤邓彬秦玉亮王宏强黎湘 . THz全尺寸凸体粗糙目标雷达回波散射建模与成像仿真. 雷达学报, 2018, 7(1): 97-107. doi: 10.12000/JR17086
  • 加载中
图(7)表(5)
计量
  • 文章访问数:  405
  • HTML浏览量:  120
  • PDF下载量:  101
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-11-04
  • 录用日期:  2018-11-28
  • 网络出版日期:  2018-12-24
  • 刊出日期:  2019-08-28

一种用于准单站配置的步进频探地雷达建模方法

    通讯作者: 方广有, gyfang@mail.ie.ac.cn
    作者简介: 张经纬(1992–),男,湖北荆州人,中国科学院电子学研究所博士研究生,主要研究方向为探地雷达系统建模及探地雷达信号处理,目前已发表多篇期刊论文及国际会议论文;刘 新(1992–),男,山东新泰人,中国科学院电子学研究所博士研究生,主要研究方向为伪随机编码雷达系统设计及其在穿墙成像中的应用;叶盛波(1985–),男,湖北黄陂人,现为中国科学院电子学研究所副研究员,主要研究方向为超宽带探地雷达系统设计,目前发表多篇论文;方广有(1963–),男,河南人,现为中国科学院电子学研究所研究员,主要研究方向为超宽带雷达成像技术、微波成像新方法和技术、 超宽带电磁场理论及工程应用,目前发表多篇高质量论文
  • ①. 中国科学院大学   北京   100049
  • ②. 中国科学院电子学研究所   北京   100190
  • ③. 中国科学院电磁辐射与探测技术重点实验室   北京   100190
基金项目:  中国科学院面上基金(CXJJ-17-M140),国家重点研发计划(2017YFF0107700)

摘要: 探地雷达是一种被广泛使用的无损检测技术。利用探地雷达对分层媒质进行全波反演时,构建精确的探地雷达正演模型具有十分重要的意义。该文提出一种可用于准单站配置的步进频探地雷达的建模方法。在该模型中,探地雷达系统及其与分层媒质间的相互作用被表示成线性方程,天线对雷达信号的影响被表示为只与频率有关的传输函数。为验证模型准确性,该文在实验室条件下搭建了准单站配置的步进频探地雷达系统,并对已知厚度的石膏板和木板的雷达测量信号进行全波反演。反演结果表明:石膏板和木板的厚度估计误差均不超过0.3 mm,验证了所提出的正演模型具有高准确度。利用石膏板和木板搭建分层模型,该文进一步比较了准单站配置和单站配置步进频探地雷达系统对介电常数差异较小的分层媒质的反演性能。实验结果表明:准单站配置探地雷达能获得更精确的反演参数。通过对分界面反射信号的信噪比估计可知,准单站配置比单站配置探地雷达系统能获得高出约10 dB的信噪比,因此具有更好的反演性能。

English Abstract

    • 探地雷达作为一种高效的无损检测工具,在国防、考古、地质勘探、道路与桥梁检测等方面具有十分广泛的应用[14]。探地雷达的回波信号解释方法可以分为两类,一类是从信号中识别反射脉冲的幅值和时延,然后依据简化的传播模型对被检测媒质的厚度和介电常数进行反演[5];另一类是先对雷达信号在被检测媒质中的传播过程进行物理建模,然后利用最优化方法对回波信号进行全波反演,得到媒质的介电常数和厚度[6]。全波反演方法是从物理角度出发模拟雷达回波信号的产生机理,因此能获得更高的检测精度。传统的探地雷达建模方法对雷达信号的传播过程作了较强的简化假设,尤其是在建模过程中忽略了天线对雷达信号的影响,包括天线的辐射模式、增益、传播时间、色散以及天线与媒质间的相互作用等。

      为了准确模拟天线系统对雷达信号的影响,许多研究人员采用数值方法对探地雷达系统进行建模,包括时域有限差分方法(FDTD)[7,8]、有限元方法(FEM)[9]、矩量法(MOM)[10]等。这些数值方法的优势在于能够处理各种复杂的媒质结构,并且不需要对被检测媒质做过多的简化假设。但是数值建模方法的劣势在于计算效率低,这使得它们在需要对正演模型进行大量的重复计算的全波反演方法中不适用。探地雷达系统在特定的应用场景下也可以利用解析方法进行建模,如针对平面分层媒质的检测。在众多的解析建模方法中,基于电场积分方程的方法是最有前景的。例如,Gentili和Spagnolini采用与频率无关的电偶极子阵列对喇叭天线进行建模[11],但是在该模型中,没有考虑雷达信号在天线口面与媒质之间的多次反射。Lambot等人[6]针对单站配置(收、发天线重合)步进频探地雷达系统提出了一种建模的方法,该方法在平面分层媒质中的雷达回波模拟和介电参数反演等方面有着出色的表现。由于常用的探地雷达系统都采用准单站配置(收、发天线分离但离得较近),本文基于单站配置探地雷达的建模思想提出了一种用于准单站配置的步进频探地雷达正演模型。在该正演模型中,发射天线和接收天线被等效为水平电偶极子,天线对雷达信号的影响以及天线与分层媒质间的相互作用被表示为只与频域有关的传输函数。为验证提出的正演模型的建模精度,本文在实验室条件下利用矢量网络分析仪和喇叭天线搭建了步进频探地雷达系统,并依据正演模型对从厚度已知的石膏板和木板测得的雷达信号进行全波反演。为进一步比较准单站配置和单站配置的探地雷达系统的反演性能,本文利用石膏板、木板和金属板构建了3层结构模型,并对两种配置的探地雷达在3层结构中得到的回波信号进行了全波反演。通过对反演结果分析,证明了准单站配置探地雷达因其具有更高的信噪比而具有更好的反演性能。

    • 利用矢量网络分析仪搭建单站步进频探地雷达系统时,发射天线与接收天线重合。当分层媒质处于探地雷达的远场区域时,天线口面处接收到的回波近似为平面波,此时在探地雷达的正下视方向上天线的增益可以假设为只与频率有关。在该假设下,单站配置的步进频探地雷达信号可以用一个只与频率有关的线性方程进行建模[6]。在该建模方法中,探地雷达的信号${S_{11}}\left( \omega \right)$可以表示为

      $ {S_{11}}\left( \omega \right) = \frac{{b\left( \omega \right)}}{{a\left( \omega \right)}} = {H_{\rm {i}}}\left( \omega \right) + \frac{{H\left( \omega \right)G_{{\rm{xx}} }^{\,0}\left( \omega \right)}}{{1 - {H_{\rm {f}} }\left( \omega \right)G_{{\rm{xx}} }^{\,0}\left( \omega \right)}} $

      其中,$b\left( \omega \right)$$a\left( \omega \right)$分别表示在矢量网络分析仪校准面上的反射信号和发射信号;${H_{{\rm i}}}\left( \omega \right)$表示由矢量网络分析仪校准面与天线馈电点间的阻抗不匹配,天线内部阻抗渐变,以及天线口面处的阻抗突变产生的综合反射效应;$H\left( \omega \right) = {H_\rm t}\left( \omega \right){H_\rm r}\left( \omega \right)$为发射天线与接收天线的传输函数的乘积,用来表示天线的增益和传播时间;${H_{\rm {f}} }\left( \omega \right)$表示发生在天线口面与媒质间的多次反射;$G_{{\rm{xx}} }^{\, 0}\left( \omega \right)$为置于天线口面处的零偏置水平电偶极子的格林函数[12],用来表示雷达信号在分层媒质中的响应;$\omega $为角频率。

    • 本文提出了一种用于准单站配置的步进频探地雷达的正演模型。考虑到在准单站配置中,发射天线与接收天线间的偏置较小,因此在远场条件下接收天线处测得的电场可以作平面波近似,同时在正下视方向上天线的增益只与频率有关。在该假设条件下,发射天线与接收天线可以分别用一个位于天线口面中心处的水平电偶极子近似;其中天线的增益、相位延迟、收发天线间的互耦以及天线与分层媒质间的相互作用可以用一组只与频率有关的传输函数表示。依据这些假设,本文提出了按照图1所示的系统框图对准单站步进频探地雷达进行建模。

      图  1  准单站探地雷达建模框图

      Figure 1.  Modeling diagram for the quasi-monostatic ground penetrating radar

      根据图1所示,步进频探地雷达的信号可以表示为

      $ \begin{align} {S_{21}}\left( \omega \right) =& \frac{{b\left( \omega \right)}}{{a\left( \omega \right)}} = {T_0}\left( \omega \right) + \frac{{{T_{\rm{i}}}\left( \omega \right){T_{\rm{s}}}\left( \omega \right){G_{{\rm{xx}} }}\left( \omega \right)}}{{1 - {R_{\rm{s}}}\left( \omega \right)G_{{\rm{xx}} }^{\,0}\left( \omega \right)}} \\ {\rm{}}\\ {\rm{}}& + \frac{{{T_0}\left( \omega \right)r\left( \omega \right){T_{\rm{s}}}\left( \omega \right)T_{\rm{i}}{\, '}\left( \omega \right)G_{{\rm{xx}} }^{\,0}\left( \omega \right)}}{{1 - {R_{\rm{s}}}\left( \omega \right)G_{{\rm{xx}} }^{\,0}\left( \omega \right)}} \end{align} $

      其中,$b\left( \omega \right)$$a\left( \omega \right)$分别表示矢量网络分析仪校准面处的反射信号和发射信号;${T_0}\left( \omega \right)$表示发射天线与接收天线间的耦合系数;${T_{\rm{i}}}\left( \omega \right)$表示发射天线的发射传输函数;${T_{\rm{s}}}\left( \omega \right)$表示接收天线的接收传输函数;$G_{{\rm{xx}} }^0\left( \omega \right)$${G_{{\rm{xx}} }}\left( \omega \right)$分别表示水平分层媒质上方零偏置和单偏置水平电偶极子的格林函数。对于准单站配置的探地雷达系统,从发射天线直接传播到接收天线的信号有一部分会被反射出来,重新形成新的发射信号,系统框图中复系数$r\left( \omega \right)$表示这部分的反射效应;$T_{\rm{i}}{\, '}\left( \omega \right)$表示接收天线充当发射天线时的发射传输函数;正反馈分量${R_{\rm{s}}}\!\left( \omega \right)$表示雷达信号在接收天线与平面分层媒质间的多次反射。式(2)可以通过定义中间变量$T\left( \omega \right) = {T_{\rm{i}}}\left( \omega \right){T_{\rm{s}}}\left( \omega \right)$, $ {\rm TR}\left( \omega \right) =$$ r\left( \omega \right)T_{\rm{i}}{\,'}\left( \omega \right){T_{\rm{s}}}\left( \omega \right)$进行化简:

      $ \begin{align} {S_{21}}\left( \omega \right) =& {T_0}\left( \omega \right) \\ {\rm{}}\\ {\rm{}}&+ \frac{{T\left( \omega \right){G_{\rm xx}}\left( \omega \right) + {T_0}\left( \omega \right){\rm TR}\left( \omega \right)G_{{\rm{xx}} }^{\, 0}\left( \omega \right)}}{{1 - {R_{\rm{s}}}\left( \omega \right)G_{{\rm{xx}} }^{\, 0}\left( \omega \right)}} \end{align} $

      (3)

      传输函数${T_0}\left( \omega \right)$, $T\left( \omega \right)$, ${\rm TR}\left( \omega \right)$${R_{\rm{s}}}\left( \omega \right)$可以通过求解式(3)得到。由于${S_{21}}\left( \omega \right)$可以直接测得,当分层媒质的参数已知时,${G_{\rm xx}}\left( \omega \right)$$G_{{{\rm xx}} }^{\,0}\left( \omega \right)$可以精确计算,因此本文采用在足够大的金属板上方远场区域不同高度处进行测量构建方程组,这里足够大的金属板可以当作无限大平面理想电导体。由于有4个待求解传输函数,理论上构建4个方程即可。为了避免求得传输系数与雷达的测量高度有关,这里进行M次测量构成冗余方程组,其中M>>4。利用${S_{21,k}}\left( \omega \right)$表示第k(k=1, 2, $ ·\!·\!·$, M)次测得的雷达信号,$G_{{\rm{xx}} ,k}^{\, 0}\left( \omega \right)$${G_{{\rm{xx}} ,k}}\left( \omega \right)$为对应的格林函数,则得到的第k个方程为

      $ {S_{21,k}} \!=\! {T_0} \!+\! {S_{21,k}}G_{{\rm{xx}} ,k}^{\, 0}{R_{\rm{s}}} \!+\! {T_0}\left( {{\rm TR} \!-\! {R_{\rm{s}}}} \right)G_{{\rm{xx}} ,k}^{\, 0} \!+\! {G_{{\rm{xx}} ,k}}T $

      (4)

      因此线性方程组可以表示为以下矩阵形式:

      $ {{y}} = {{Bx}} $

      其中

      $ {{y}} = {\left[ {{S_{21,1}} ·\!·\!· {S_{21,k }}·\!·\!· {S_{21,M}}} \right]^{\rm T}}\hspace{60pt} $

      $ {{B}} = \left[ {\begin{array}{*{20}{c}} 1&{{S_{21,1}}G_{{\rm{xx}} ,1}^{\, 0}}&{G_{{\rm{xx}} ,1}^{\, 0}}&{{G_{{\rm{xx}} ,1}}} \\ \vdots & \vdots & \vdots & \vdots \\ 1&{{S_{21,k}}G_{{\rm{xx}} ,k}^{\, 0}}&{G_{{\rm{xx}} ,k}^{\, 0}}&{{G_{{\rm{xx}} ,k}}} \\ \vdots & \vdots & \vdots & \vdots \\ 1&{{S_{21,M}}G_{{\rm{xx}} ,M}^{\, 0}}&{G_{{\rm{xx}} ,M}^{\, 0}}&{{G_{{\rm{xx}} ,M}}} \end{array}} \right] $

      $ {{x}} = {\left[ {\begin{array}{*{20}{c}} {{T_0}}&{{R_{\rm{s}} }}&{{T_0}\left( {{\rm TR} - {R_{\rm{s}} }} \right)}&T \end{array}} \right]^{\rm T}}\hspace{20pt} $

      未知参数向量${{x}}$可以利用最小二乘方式进行求解

      $ {{x}} = {\left( {{{{B}}^{\rm H}}{{B}}} \right)^{ - 1}}{{{B}}^{\rm H}}{{y}} $

      其中,H表示共轭转置操作。

    • 分层媒质的参数反演问题本质上是一个寻找全局最优解${{z}} = \left[ {{\varepsilon _u},{\sigma _u},{h_u}} \right]$, u=1, 2, ···, N的优化问题,其中,${\varepsilon _u}$表示第u层媒质的介电常数,${\sigma _u}$表示电导率,${h_u}$表示第u层的厚度。通过定义目标函数${\rm{obj}} \left( {{z}} \right)$,则参数反演的过程就是使目标函数最小化的过程。由于探地雷达目标的反射信号限制在一定时间范围内,因此全波反演过程一般在时域进行。目标函数可以定义为

      $ {\rm{obj}} \left( {{z}} \right) = \frac{{{{\left| {{{{s}}^*}\!\left( t \right) - {{s}}\!\left( t \right)} \right|}^{\rm T}}\left| {{{{s}}^*}\!\left( t \right) - {{s}}\!\left( t \right)} \right|}}{{\left\| {{{{s}}^*}\!\left( t \right) - {{s}}\!\left( t \right)} \right\|}} $

      $ {{{s}}^*}\!\left( t \right) = {\rm{IFFT}}\left[ {{{{S}}^*}\!\left( \omega \right)} \right]\hspace{75pt} $

      $ {{s}}\!\left( t \right) = {\rm{IFFT}}\left[ {{{S}}\!\left( {\omega ,{{z}}} \right)} \right]\hspace{76pt} $

      其中,${{{s}}^*}\!\left( t \right)$是包含时域反演目标信号的向量;${{s}}\!\left( t \right)$是反演过程中利用未知参数向量${{z}}$仿真得到的时域波形向量;${\rm{IFFT}}\left[ \cdot \right]$表示逆傅里叶变换;${{{S}}^*}\!\left( \omega \right)$表示探地雷达系统测得的信号;${{S}}\!\left( {\omega ,{{z}}} \right)$表示仿真得到的探地雷达信号;$\left\| \cdot \right\|$表示向量的二范数;$\left| \cdot \right|$表示求模操作;上标T表示转置操作。在多数的电磁反演问题中,目标函数都是极度非线性的,其解空间中存在很多局部最优解,因此一般电磁反演问题都使用全局最优化算法,如遗传算法、差分进化算法等。

    • 本文利用便携式矢量网络分析仪模块S5048[13]和2个喇叭天线搭建步进频探地雷达系统。其中,矢网的工作频段为0.9~4.0 GHz,共256 个频率点,中频带宽为3 kHz,在测试前矢网经过校准件进行校准。喇叭天线长24 cm,口面尺寸为16 cm×30 cm,有效带宽为0.8~8.0 GHz。2个天线呈E面平行放置,天线口面中心点相距20 cm。为了得到正演模型的传输函数,利用该雷达在一个2.5 m×2.5 m大的金属板上方0.48~0.75 m高度范围内进行55 次不同高度的测量,如图2所示。为了量化建模误差,可以定义一个建模误差函数:

      图  2  准单站配置和单站配置探地雷达系统

      Figure 2.  Configurations of quasi-monostatic and monostatic GPR

      $ {j} \left( {{x}} \right) = \frac{{\left| {{{y}} - {{{y}}^*}} \right|{{\left| {{{y}} - {{{y}}^*}} \right|}^{\rm T}}}}{{{{\left\| {{y}} \right\|}^2}}} $

      其中,y是各个高度下金属板实测信号,${{{y}}^*}$是各个高度下金属板实测信号的最优逼近值,${{x}}$是传输函数的最小二乘解。为了验证模型中的参量$r\left( \omega \right)$的作用,可以通过将式(2)中的$r\left( \omega \right)$置零,得到一个简化的模型:

      ${S_{21}}\left( \omega \right) = {T_0}\left( \omega \right) + \frac{{{G_{{\rm{xx}} }}\left( \omega \right)}}{{1 - R_{\rm{s}}\! ' \!\left( \omega \right)mG_{{\rm{xx}} }^{\, 0}\left( \omega \right)}}{T{\; '}}\!\left( \omega \right) $

      利用同样的雷达测量数据可以对上述简化模型的传输参数进行求解。图3展示了参数简化模型和原始模型在各个频率上的建模误差。结果表明,去掉模型参量$r\left( \omega \right)$后,雷达系统的建模误差会上升,因此在模型中考虑参量$r\left( \omega \right)$的贡献对提升探地雷达系统建模的准确度是十分重要的。

      图  3  简化模型和原始模型建模误差

      Figure 3.  Fitting errors of the simplified model and original model

    • 为验证提出的正演模型的准确度,本文利用搭建的探地雷达系统分别在厚为28.4 mm面积为1.6 m×1.8 m的石膏板和厚度为19.0 mm面积为1.8 m×2 m的干木板上方不同高度处进行测量,并且用一个面积为2.5 m×2.5 m的覆铜板用来充当底部的边界条件,如图4所示。在对测得的雷达信号进行参数反演时,假定石膏板和干木板的介电常数与频率无关,且两者的电导率均可忽略。图5(a)图5(b)分别展示了利用反演参数得到的仿真波形和实测波形,从图中可以看出利用反演结果得到的仿真波形和实测波形具有很高的吻合度。

      图  4  利用厚度已知的石膏板验证准单站探地雷达模型的精度

      Figure 4.  Validating the accuracy of the quasi-monostatic model with the thickness-known plasterboard

      图  5  石膏板与木板的实测波形与仿真波形的对比

      Figure 5.  Comparison of the measured and simulated signal for plasterboard and woodblock

      表1表2分别展示了对石膏板和干木板的雷达测量波形进行全波反演得到的结果。

      高度(mm)反演相对介电常数反演厚度(mm)厚度估计误差(mm)
      534.32.3228.70.3
      522.52.3328.30.1
      515.32.3228.60.2

      表 1  准单站探地雷达石膏板反演结果

      Table 1.  Inversion results of the plasterboard with the quasi-monostatic GPR

      高度(mm)反演相对介电常数反演厚度(mm)厚度估计误差(mm)
      531.52.7219.20.2
      523.22.7318.80.2
      512.42.7418.70.3

      表 2  准单站探地雷达木板反演结果

      Table 2.  Inversion results of the woodblock with quasi-monostatic GPR

      表1中,反演得到的石膏板的平均相对介电常数为2.32,标准差为0.006。在表2中,反演得到的木板的平均相对介电常数为2.73,标准差为0.009。石膏板和木板的最大厚度反演误差均为0.3 mm。在表1表2中反演得到的介电常数的值都比较稳定,且对应的反演厚度误差都非常小,因此利用反演得到的相对介电常数的平均值代表其真实值是合理的。表1表2的反演结果说明本文提出的准单站步进频探地雷达的正演模型具有很高的精确度。为了进一步验证,在原始模型中考虑参数$r\left( \omega \right)$能够提升模型的精度,本文利用简化模型对同样的石膏板测量数据进行了全波反演。表3中为利用简化模型反演得到的石膏板参数。对比表1表3的结果可知,由简化模型反演得到的石膏板的介电常数略微大于原始模型得到的介电常数。这是因此,利用简化模型进行全波反演时,实际由原始模型中参数$r\left( \omega \right)$对反射波形产生的部分贡献也被错误地归因于简化模型中的其他参数,导致简化模型中得到的石膏板反射波形比实际石膏板的反射波形强度大,即等效反射系数大,而在进行全波反演时,利用的是最小化均方误差的优化准则,因此最终通过最优化方法反演得到的石膏板介电常数也会略大于真实值。

      高度(mm)反演相对介电常数反演厚度(mm)厚度估计误差(mm)
      533.82.4328.00.4
      522.02.4527.60.8
      514.72.4427.90.5

      表 3  简化模型对石膏板的反演结果

      Table 3.  Inversion results of the plasterboard with simplified model

    • 根据前人的研究成果和本文的实验结果表明,单站配置和准单站配置步进频探地雷达系统都能被精确地建模,且能获得良好的反演性能。为了进一步比较准单站配置和单站配置探地雷系统在介电常数差异较小的分层媒质中的反演性能,本文利用相同的矢量网络分析仪和天线在不改变频率测量点和中频带宽等配置参数的条件下搭建了单站步进频探地雷达系统。在实验室条件下利用相同的石膏板和木板搭建了3层结构模型,其中第1层为空气,第2层为石膏板,第3层为木板,并且用2.5 m×2.5 m的覆铜板充当分层媒质的底部边界条件。同样利用单站配置和准单站配置的探地雷达分别在3层结构上方不同高度处进行测量,并对测量波形进行反演。表4表5分别展示了准单站配置和单站配置的雷达测量波形的反演结果。从结果中可以看出,对于两个不同配置的探地雷达,反演得到的石膏板的介电常数都非常接近其真实值。对于木板的介电常数,由准单站配置的探地雷达得到的反演结果更接近真实值,且标准差更小。在石膏板和木板的厚度反演中,由准单站配置探地雷达系统得到的厚度反演误差整体更小。表4表5的结果表明,当分层媒质的介电常数差异较小时,准单站配置的探地雷达系统具有更好的反演性能。影响全波反演结果精度的一个主要因素就是信号的信噪比。为了理解准单站配置探地雷达系统的优势,这里分析了两个雷达系统获得的信号的信噪比。在单站配置中,矢量网络分析仪采用反射测量模式,其端口需要采用定向耦合器进行收发信号分离。由于定向耦合器的隔离度一般只有23 dB,因此在接收信号中有一部分是从发射端串扰过去的噪声。经校准过的矢量网络分析仪能获得有效的方向性增益约为46 dB[13],所以矢网在进行反射测量时其系统内部归一化噪声主要受系统的方向性增益影响。当矢网采用传输测量方式工作时,由于收发通道在物理上是隔离的,因此其归一化噪声主要受中频带宽影响。对于校准后的矢网,在其接收端口接上匹配负载可以实际测得其内部归一化的噪声水平。另外,当分层媒质的厚度和介电常数已知时,可以利用正演模型有效估计在不同的雷达配置下测量到的来自各分界面反射信号的归一化功率。这里利用图6(a)所示的3层结构模型来说明如何估计来自分界面的反射信号,由该3层媒质测得的反射信号可以看作是分别由分界面${z_2}$${z_3}$产生的反射信号的叠加。

      天线高度
      (mm)
      反演石膏板
      介电常数
      石膏板反演厚度
      (mm)
      石膏板厚度误差
      (mm)
      反演木板
      介电常数
      木板反演厚度
      (mm)
      木板厚度误差
      (mm)
      532.22.3428.10.32.7119.50.5
      524.22.3228.80.42.7219.40.4
      510.32.3128.90.52.7518.50.5

      表 4  准单站配置雷达3层结构反演结果

      Table 4.  Inversion results of the three-layer structure with quasi-monostatic GPR

      天线高度
      (mm)
      反演石膏板
      介电常数
      石膏板反演厚度
      (mm)
      石膏板厚度误差
      (mm)
      反演木板
      介电常数
      木板反演厚度
      (mm)
      木板厚度误差
      (mm)
      535.32.3127.80.62.6320.11.1
      525.22.3028.90.52.7818.10.9
      513.32.3327.70.72.6519.80.8

      表 5  单站配置雷达3层结构反演结果

      Table 5.  Inversion results of the three-layer structure with monostatic GPR

      图  6  分界面反射信号的计算模型

      Figure 6.  Models used to calculate the reflections on the interface

      如果介电常数${\varepsilon _2}$${\varepsilon _3}$相等,那么在分界面${z_3}$上将不会产生反射。此时3层分层模型退化成了2层分层模型,如图6(b)。这里用符号$S\left( \omega \right)$${S_1}\left( \omega \right)$分别表示由3层模型和2层模型得到雷达信号,因此单独来自于3层模型中的分界面${z_3}$的反射信号${S_2}\left( \omega \right)$可以表示为

      $ {S_2}\!\left( \omega \right) = S\!\left( \omega \right) - {S_1}\!\left( \omega \right) $

      假定天线高度为525 mm,利用表1表2中反演的石膏板(${\varepsilon _2}$)和木板(${\varepsilon _3}$)的介电常数,以及两者的真实厚度(${h_2},{h_3}$),依据式(15)可以准确计算来自石膏板和木板间的分界面的反射信号。图7(a)图7(b)分别展示了两种雷达配置下得到的归一化反射信号和归一化内部噪声,以及对应的信噪比。可以看出,准单站配置探地雷达能获得更高的信噪比,因此在反演中表现出更好的性能。

      图  7  单站配置探地雷达和准单站配置探地雷达的信噪比对比

      Figure 7.  The comparison of the SNR between the monostatic and quasi-monostatic GPR

    • 本文提出了一种用于准单站配置的步进频探地雷达的正演模型。通过在已知厚度的石膏板和木板上方不同高度处进行测量,并对测得的雷达回波进行全波反演,证明了该正演模型具有很高的建模精度。通过利用石膏板和木板构建参数已知的3层模型,并分别构建单站配置探地雷达和准单站配置探地雷达进行参数反演,证明了当分界面上下介质介电常数差异较小时,准单站配置探地雷达因为其底噪更低能获得更高的信噪比,因此反演性能更好。虽然本文中是利用VNA测试的结果验证准单站配置探地雷达比单站配置探地雷达系统信噪比高,反演性能好,但是该结论同样适用于其他形式构建的步进频探地雷达系统。

参考文献 (13)

目录

    /

    返回文章
    返回