文章快速检索     高级检索
  浙江大学学报(工学版)  2018, Vol. 52 Issue (9): 1828-1836  DOI:10.3785/j.issn.1008-973X.2018.09.025
0

引用本文 [复制中英文]

刘钰, 石战结, 王帮兵, 王征宇, 余天祥. 探地雷达子波确定性稀疏脉冲反褶积技术[J]. 浙江大学学报(工学版), 2018, 52(9): 1828-1836.
dx.doi.org/10.3785/j.issn.1008-973X.2018.09.025
[复制中文]
LIU Yu, SHI Zhan-jie, WANG Bang-bing, Wang Zheng-Yu, YU Tian-xiang. Deterministic-wavelet sparse spike deconvolution technique for ground penetrating radar data[J]. Journal of Zhejiang University(Engineering Science), 2018, 52(9): 1828-1836.
dx.doi.org/10.3785/j.issn.1008-973X.2018.09.025
[复制英文]

基金项目

国家社会科学基金重大招标资助项目(13&ZD192);国家自然科学基金资助项目(41676181);中央高校基本科研业务费专项资助项目(ZJUSK2016007,2018QNA3013);浙江省文物保护科技资助项目(2013011,2015019,2017013);中石化石油工程地球物理有限公司胜利分公司科技资助项目(10200274-17-ZC0613-0018)

作者简介

刘钰(1992—),男,硕士,从事探地雷达数据处理及反演方法研究.
orcid.org/0000-0003-0391-7212.
E-mail:liu_yu@zju.edu.cn.

通信联系人

石战结,男,副教授.
orcid.org/0000-0001-9283-9070.
E-mail:shizhanjie@zju.edu.cn
.

文章历史

收稿日期:2017-07-08
探地雷达子波确定性稀疏脉冲反褶积技术
刘钰1, 石战结2,1, 王帮兵1, 王征宇3, 余天祥1     
1. 浙江大学 地球科学学院,浙江 杭州 310027;
2. 浙江大学 文化遗产研究院,浙江 杭州 310058;
3. 杭州市文物考古研究所,浙江 杭州 310008
摘要: 为了提升探地雷达数据的纵向分辨率,并解决传统反褶积方法对噪声敏感的问题,提出一种针对探地雷达信号的稀疏脉冲反褶积处理方法. 通过试验获取子波信息,基于确定性的子波,采取正则化的反演方法,从信号中获得主要的反射系数. 通过2个数值试验以及1个实际案例,对比分析处理前后的单道信号、剖面、频谱等情况,并与传统反褶积方法的结果做比较,结果表明:该方法相对于传统方法能较大地改善分辨率,并且受噪声影响较小,得到的结果成像效果更好. 该方法经验证是一项有效、可行的探地雷达信号处理技术.
关键词: 探地雷达(GPR)    子波    反褶积    反演    高分辨率    
Deterministic-wavelet sparse spike deconvolution technique for ground penetrating radar data
LIU Yu1 , SHI Zhan-jie2,1 , WANG Bang-bing1 , Wang Zheng-Yu3 , YU Tian-xiang1     
1. School of Earth Sciences, Zhejiang University, Hangzhou 310027, China;
2. Institute of Cultural Heritage, Zhejiang University, Hangzhou 310058, China;
3. Hangzhou Institute of cultural relics and Archaeology, Hangzhou 310008, China
Abstract: A sparse spike deconvolution method for GPR signals was put forward in order to improve the vertical resolution of ground penetrating radar (GPR) data and solve the noise-sensitive problem of traditional deconvolution technique. Firstly the wavelet was acquired by experiment. Then a constraint inversion method using the deterministic wavelet was carried out to obtain the dominating reflectivity series from signals. The signals of single channel, profiles, spectrums, etc. before and after the process were analyzed through two numerical simulation experiments, as well as a practical case. Results were compared with that of the traditional deconvolution. The analyzing and comparing results show that the method can improve the resolution very well and is less affected by noise, which can provide a better image of the target. Our method is proved to be an effective and feasible GPR signal processing technique.
Key words: ground penetrating radar (GPR)    wavelet    deconvolution    inversion    high resolution    

探地雷达(ground penetrating radar, GPR)通常利用高频电磁波在媒质电磁特性不连续处的反射达到近地表地球物理探测的目的[1],是一项快速无损且分辨率较高的方法. 与其他调查方法相比,GPR不仅在异常体平面展布信息上有较好反映,还可以提供更多深度上的信息,这是GPR在考古、环境调查和工程等应用中广受青睐的重要原因[2]. 然而,GPR发射的脉冲信号在传播过程中会因高频成分的衰减而成为一个带宽有限的子波,这导致原始资料的分辨率受到限制. 为了解决这个问题,需要对信号作反褶积处理,以压缩子波、拓宽频带,完成纵向分辨率的提高.

在地震勘探领域,已经有很多关于反褶积的研究和应用[3-4]. 随着GPR的广泛应用,越来越多的反褶积技术被尝试用于GPR数据处理. 首先是一些基于统计学子波的方法[5-8],在一些领域取得了较好的效果[9-11]. 后来Xia等[12-13]在做沉积研究时,认为这些方法依据的假设很难适用,提出了利用确定性子波作反褶积,并成功将其应用于碳酸盐地层研究上. 不过,Annan[14]认为这些经典的地震反褶积方法在大量实际工作中还是鲜有好的应用效果,并指出部分原因在于子波高频成分在传播过程中的衰减,需要用反Q滤波等技术对此进行补偿[15],后来国内苏茂鑫等[16-17]也针对补偿方法作了有益的尝试. 近年来,一些新的方法被陆续提出,包括基于贝叶斯估计子波的频率域反褶积[18]、时变子波的反褶积[19]等,国内也不乏相关的研究[20-24]. 不过随着GPR应用的逐渐推广,反褶积技术的发展仍然赶不上分辨率的需求,例如薄层识别的问题变得越来越多[25]. 值得注意的是,一些反褶积方法基于反射系数可由稀疏脉冲表示的假设,对分辨率有着明显的提升,不过对子波的处理则仍然采用估计的策略[18, 26-27].

GPR剖面常常包含各种成因复杂的噪声,而常规的反褶积方法对噪声比较敏感,这也是现在反褶积技术没有大规模推广应用到实际工作中的主要原因之一. 另外,传统的反褶积需要假设子波是最小相位以及反射系数序列是白噪声序列,这在GPR工作中经常是难以满足的. 如果不考虑这2种假设,为了提取地下复杂介质信息的主要特征,将地下地层的反射系数看作可由一系列稀疏分布的脉冲表示,则GPR信号可以理解为子波与一个综合反射系数序列褶积的结果,这个综合反射系数序列包括了稀疏分布的、代表地下地层的大反射系数和随机分布、由噪声引起的小反射系数. 稀疏脉冲反褶积就是基于该条件约束的一种反演策略[28-32],其对噪声的敏感度较小,通过迭代优化能够从信号中恢复主要反射系数的准确位置和大小,并且突破了分辨率的限制,极大提升了频带的宽度. 在地震勘探领域,稀疏脉冲反褶积已被广泛用于石油工业的波阻抗反演当中[33],对于子波的获取,主要是靠井震标定工作来完成. 关于GPR波阻抗反演的研究近年来也有所发展[34-38],但对于子波的处理方法尚无细致的研究. 不同于地震勘探中震源激发条件的差异化,GPR特定型号的主机和天线所激发的子波应该是比较稳定的. 虽然GPR的子波相比地震子波面临着更大的频率损失,但通过前文提到的补偿方法,仍然能够假设子波是稳定的、时不变的. 因此本文提出,通过类似Xia等[12-13]的实验获取确定性子波,避免估计子波因噪声干扰等问题引起的误差;利用地震勘探中经典的约束和优化策略做稀疏脉冲反褶积,从含噪信号中提取高分辨率的有效信息,甚至突破GPR信号理论上最小分辨率的限制.

本文介绍GPR确定性子波的获取手段,以及稀疏脉冲反褶积的方法原理;通过对2个理论模型的分析,验证该方法在抗噪能力和薄层识别效果上与传统方法相比的优势;利用该方法对实际数据进行处理,得到高分辨率的结果,证明该方法的可行性.

1 确定性子波的获取

不同的GPR产品有着不同的仪器系统,因此需要对特定型号的主机以及所选天线作专门的试验和研究. 本文以PulseEKKO pro系统搭配50 M板状天线为例. 为了获取干净的GPR子波,并且减少能量衰减和频率损失,需要选择一个空旷的场地,避免来自周围环境或者人文噪声的干扰,将收发天线垂直于地面以固定间距互相正对底面(如图1所示),以便采集到空气中的直达波. 这种情况下,来自草地的反射相对于空气直达波来说非常弱,可以忽略[12]. 天线之间的最佳间隔则需要经过试验来确定,对于PulseEKKO pro系统,则可以通过搜索模式来完成. 如果间隔较近,可能信号中会包含电磁波在天线之间周期性回荡而产生的交混回响,因此需要调试不同间距的效果. 选择好最优间距后,在实际采集时需要选择较多的叠加次数,从而使随机噪声的影响减少到最小. 具体的例子将在第4章进行介绍.

图 1 探地雷达50 MHz天线确定性子波的数据采集工作照片 Fig. 1 Work photo of acquiring GPR 50 MHz antennas' deterministic wavelet
2 稀疏脉冲反褶积

根据褶积模型,GPR信号 ${{d}_t}$ 可以看作是一个有限带宽的静态子波 ${{{w}}_t}$ 与地下反射系数序列 ${{r}_t}$ 的褶积,又因为野外采集不可避免会遇到噪声 ${n_t}$ nt,即有

${{{d}}_t} = {{{w}}_t}{\rm{*}}{{{r}}_t} + {{{n}}_t},\, t = 1,2, \cdots ,{N} .$ (1)

从式:可以看出,信号的频带受到了子波的影响. 为了从 ${{d}_t}$ 中恢复 ${{r}_t}$ ,可以通过已知子波或者统计学方法,设计一个相当于褶积作用反过程的滤波器并与GPR信号褶积. 最经典的反褶积方法是基于维纳滤波器的最小平方反褶积[3],滤波器的设计准则是使其实际输出与期望输出的误差平方和最小. 不过,该方法假设地震反射系数序列是白噪声序列,并且子波是最小相位的,从而在未知子波的条件下用记录的自相关和与期望输出的互相关,代替子波的自相关和与期望输出的互相关. 当期望输出为一个零延时或有一定延时脉冲时,该方法被称为脉冲反褶积. 实际情况中前面2个假设通常难以成立,另外GPR信号不可避免会包含未知的噪声,因此设计出滤波器的效果经常难以令人满意.

根据式(1)求取 ${{r}_t}$ 可以被认为是一个高度病态的反问题,需要采取正则化的方法来求解[29],下面介绍具体的过程. 首先假设反射系数序列是稀疏分布的,计算值和观测值的残差为

${{{e}}_t} = \mathop \sum \limits_i \,{{{w}}_{t - i}}{{{r}}_i} - {{{d}}_t},\,i = 1,2, \cdots ,{{N}}.$ (2)

反问题的目的就是找到符合条件的 ${{r}_t}$ ,使下面的目标函数最小:

${{J}} = {{{J}}_r} + {\rm{\mu }}{{{J}}_x}.$ (3)

该目标函数的第一项:

${{{J}}_r} = \mathop \sum \limits_t\, \frac{1}{2}{{{e}}_t}^2,$

对残差进行最小平方约束,保证所求结果与GPR信号符合;第二项为正则化项,对反射系数进行稀疏约束,能够保证求解结果稀疏,提升解的稳定性, ${{\mu }}$ 称作阻尼系数,需要选取适当的值进行稀疏控制. ${{{J}}_x}$ 的选取可以有不同的策略,本文采取L1范数正则化方法,即

${{{J}}_x} = \mathop \sum \limits_t \,\left| \,{{{{r}}_t}} \right|.$ (4)

为J达到最小,可由解决下式来完成:

$\frac{{\partial {{J}}}}{{\partial {{{r}}_i}}} = \mathop \sum \limits_t \,\left({{{e}}_t}{{{w}}_{t - i}}\right) + {\rm{\mu }}\mathop \sum \limits_i \,\frac{{\left| {{{{r}}_i}} \right|}}{{{{{r}}_i}}} = 0.$ (5)

根据上式可以推导出:

$\mathop \sum \limits_t \mathop \sum \limits_i\, \Big[\Big({{{w}}_{t - i}}{{{w}}_{t - i}} + {\rm{\mu }}\frac{1}{{\left| {{{{r}}_i}} \right|}}\Big){{{r}}_i}\Big] = \mathop \sum \limits_t\, \Big({{{w}}_{t - i}}{{{d}}_t}\Big).$ (6)

将上式转化成矩阵形式,有:

$\left( {{{R + Q}}} \right){{x = g}},\,{{Q}} = {\rm{diag}}\,\left( {{\rm{\mu }}\frac{1}{{\left| {{{{r}}_i}} \right|}}} \right).$ (7)

其中,diag表示将括号中的列向量转换成相应的对角矩阵; ${{x}}$ 即为反射系数序列 ${{r}_t}$ 的列向量. 根据式(7)求解 ${{x}}$ 的方法有很多,由于GPR数据的规模通常并不大,矩阵的求逆运算是可以接受的,本文采用迭代重复加权最小二乘法算法,以求取较高精度的解. 该方法的步骤具体如下:

1)根据已知子波和GPR信号计算矩阵 ${{R}}$ ${{g}}$ ,选取适当的阻尼系数 ${\rm{\mu }}$ ,利用脉冲反褶积结果 ${{{x}}_s}$ 计算初始的矩阵:

${{{Q}}^{\left( 0 \right)}} = {\rm{diag}}\,\,\left( {{\rm{\mu }}\frac{1}{{\left| {{{{x}}_s}} \right|}}} \right).$

2)根据 ${{{Q}}^{\left( {{\rm{k}} - 1} \right)}}$ 计算 ${{{x}}^{\left( {\rm{k}} \right)}} = {\left( {{{R}} + {{{Q}}^{\left( {{\rm{k}} - 1} \right)}}} \right)^{ - 1}}{{g}}$ ,并迭代更新:

${{{Q}}^{\left( {\rm{k}} \right)}} = {\rm{diag}}\,\,\left( {{\rm{\mu }}\frac{1}{{\left| {{{{x}}^{\left( {\rm{k}} \right)}}} \right|}}} \right).$

其中,k为迭代的次数;

3)给出一个固定的次数,步骤2)的迭代次数达到该值后停止并输出结果 ${{x}}$ ; 也可以设置一个阈值,当步骤2)结果的目标函数小于该值时停止迭代.

3 数值试验

为了验证子波确定条件下的稀疏脉冲反褶积效果,设计2个简单的层状模型,如图2所示. 模型中共有A和B两种介质:背景介质A的相对介电常数为16,电导率为0.001 S/m,相对磁导率为1;目标水平层B的相对介电常数为9,电导率为0.001 S/m,相对磁导率为1. 考虑到GPR信号理论上垂直分辨率极限为1/4波长(这里即2 m),设置图2(a)中的介质B层厚为1个波长,图2(b)中的介质B层厚为1/8个波长. 正演模拟采用时间域有限差分的方法[1],设置激励源为50 MHz中心频率的雷克子波[39],采用2 m收发距以共偏移距方式采集.

图 2 数值模拟采用的层状模型 Fig. 2 Layer models used for numerical simulation

图2(a)模型产生的数据经过增益和去直达波处理后(如图3(a)所示,其中A为信号振幅),对各道数据加上20%主频50 MHz的随机噪声(如图3(b)所示). 考虑到各道信号的相似性,并且为了放大细节特征,下面主要分析剖面中抽取的第5道信号(如图3(c)所示). 图3(d)显示了对该道作常规的脉冲反褶积处理后的结果,可以看出,由于噪声的影响,产生了很多没有意义的尖脉冲. 数据分辨率提升的同时,信噪比大大降低,数据解译工作变得更加困难. 即使对数据再作频域滤波处理(如图3(e)所示),也很难挽回信噪比的损失,并且真实反射系数位置附近的信号比较杂乱,已经很难推知真实的信息. 再对图3(c)数据基于已知的子波作稀疏脉冲反褶积(如图3(f)所示),得到的结果只有2个大小相近、方向相反的大脉冲和1个较小的脉冲. 2个大脉冲即真实的反射序列,而小脉冲则是由于噪声引起的误差,不过这仍在可以接受的范围内. 不难看出,在稀疏控制项的约束下,稀疏脉冲反褶积的结果对噪声不那么敏感,可以在噪声较为严重的情况下较为准确地找出最主要的反射系数.

图 3 稀疏脉冲反褶积与常规脉冲反褶积的抗噪能力对比 Fig. 3 Comparison on anti-noise capability between sparse spike deconvolution and conventional spike deconvolution

图2(b)模型产生的数据经过增益和去直达波处理后(如图4(a)所示),剖面仍然难以反映薄层的信息. 从剖面中抽取第5道信号(如图4(b)所示),可以发现由于层厚小于理论上的纵向分辨率的极限,来自水平层上 Anal 下界面的反射波发生了比较严重的干涉,传统的处理方法已经很难发挥作用. 对该道信号基于已知子波作稀疏脉冲反褶积处理,结果得到了分离得比较好的2个尖脉冲,可以很好地对薄层的上 Anal 下界面进行识别(如图4(c)所示). 稀疏脉冲反褶积前后的谱振幅表明,处理后的信号在频带被拓宽的基础上,有效带宽上的能量也得到增强(如图5所示,Amp为信号的振幅比例). 上面的分析说明,稀疏脉冲反褶积通过已知子波反演真实的反射系数,能够大大地提升数据的分辨率,甚至突破GPR信号理论上最小分辨率的限制.

图 4 稀疏脉冲反褶积的薄层识别能力分析 Fig. 4 Analysis on sparse spike deconvolution’s capability of thin-layer recognition
图 5 稀疏脉冲反褶积前、后振幅谱对比 Fig. 5 Comparison of spectrums before and after sparse spike deconvolution
4 实例分析

本研究的实例工作区域位于杭州市吴越国捍海塘某发掘地点(如图6所示). 发掘工作表明,捍海塘主体为土塘,塘体之后有粗加工过的木材框架,以土、石为填充材料,结合竹篱、竹编等辅助材料形成整体,从而加固捍海塘[40].

图 6 捍海塘某发掘地点照片 Fig. 6 Photograph of excavated seawall

在工作区域布置一条测线,横切发掘后回填的捍海塘所在位置. 根据已知的捍海塘大致深度,选取50 MHz天线进行共偏移距测量. 在探测工作开始前,利用第1章所讲的方法试验1 、1.5 和2 m的天线间距,最终选择2 m间距采集了子波的信号. 应该注意的是,对原始子波的处理是有必要的,因为作稀疏脉冲反褶积之前信号需要作基本处理,子波处理要参考GPR信号的处理过程来进行. 这里对子波信号作零时校正、Dewow、预测反褶积和频域滤波处理. 信号中子波后面的噪声被压制但仍存在,为避免子波拖尾带来的误差,切除60 ns之后的非有效信息,得到最终的确定性子波(如图7所示),供后面稀疏脉冲反褶积所用.

图 7 获取到的确定性子波 Fig. 7 Obtained deterministic wavelet

采集到的原始剖面经过基本处理后(如图8(a)所示)已能较好地反映捍海塘石质结构的顶底界面,即21~26 m处位于约110 ns和170 ns的反射波组. 由于采用的是低频天线,信号的波长比较长,石质结构顶界面的细致位置还无法查明. 另外一些噪声,如剖面右部分的线性相干噪声[41],虽然被压制但仍旧明显,影响后面的解译工作. 通过上面获取到的确定性子波以及第2章所介绍的稀疏脉冲反褶积方法,对该剖面作子波确定性稀疏脉冲反褶积处理.

图 8 子波确定性稀疏脉冲反褶积前、后剖面对比 Fig. 8 Comparison of profiles before and after deterministic-wavelet sparse spike deconvolution

处理结果更为清晰地反映了捍海塘石质结构顶部界面的位置(如图8(b)所示),并且没有受到线性相干噪声的影响. 虽然一些微弱的有效信号受到了压制,但总体的成像效果更好. 观察处理前、后的振谱幅(如图9所示)可知,经过稀疏脉冲反褶积后的数据频带被明显拓宽,包含的信息较之前更为丰富. 综合上述分析,认为子波确定性稀疏脉冲反褶积针对该剖面的处理是十分有意义的.

图 9 子波确定性稀疏脉冲反褶积前、后振幅谱对比 Fig. 9 Comparison on spectrums before and after deterministic-wavelet sparse spike deconvolution
5 结 论

(1)理论和实践表明,本文提出的探地雷达子波确定性稀疏脉冲反褶积技术在分辨率和信噪比上都能够有较好的效果. 采集到的确定性子波较为可信,以此对信号作稀疏脉冲反褶积处理,能够大大地提升数据的纵向分辨率;在稀疏控制项的约束下,反演得到的解比较稳定,对噪声敏感度低,可以在噪声较为严重的情况下较为准确地找出最主要的反射系数.

(2)本文作稀疏脉冲反褶积采取的策略是L1范数约束,迭代重复加权最小二乘算法求解优化问题. 反演方法的创新可以推动本研究的发展,因此可以尝试使用其他正则化方法约束,以及其他的最优化算法,或者与其他方法的数据联合反演[42],这里不再赘述.

(3)所提出的方法基于确定性子波,即假设子波是时不变的,这与补偿方法的研究密切相关. 另外基于时变子波对探地雷达信号作反褶积也是今后一个有意义的方向.

(4)稀疏脉冲反褶积虽然有一定的抗噪能力,但对于噪声问题过于严重的剖面还是很难奏效. 因此需要将稀疏脉冲反褶积和一些滤波方法有机地组合起来,方能发挥最佳效果.

(5)稀疏脉冲反褶积得到的反射系数如果质量较好,经过递推或道积分后可以得到电磁波波阻抗,进而推出介电常数、含水量等参数,这对于地球物理探测来说具有更大的意义.

参考文献
[1]
曾昭发, 刘四新, 冯晅. 探地雷达原理与应用[M]. 北京: 电子工业出版社, 2010: 1-2. http://www.bookask.com/book/197917.html
[2]
SHI Z, TIAN G, HOBBS R W, et al. Magnetic gradient and ground penetrating radar prospecting of buried earthen archaeological remains at the Qocho City site in Turpan, China[J]. Near Surface Geophysics, 2015, 13(5): 1-8.
[3]
YILMAZ Ö. Seismic data analysis: Processing, inversion, and interpretation of seismic data [M]. Tulsa: Society of Exploration Geophysicists, 2001. https://www.researchgate.net/publication/49181873_Seismic_data_analysis
[4]
石战结, 田钢, 赵文轲, 等. 超浅层三维地震勘探技术应用[J]. 浙江大学学报: 工学版, 2013, 47(5): 912-917.
SHI Zhan-jie, TIAN Gang, ZHAO Wen-ke, et al. Application on ultra-shallow 3D seismic exploration technology[J]. Journal of Zhejiang University: Engineering Science, 2013, 47(5): 912-917.
[5]
JUDGE A S, TODOESCHUCK J P, LAFLECHE P T, et al. Analysis of ground-probing radar data: predictive deconvolution[J]. Canadian Geotechnical Journal, 1991, 28(1): 134-139. DOI:10.1139/t91-014
[6]
TURNER G. Subsurface radar propagation deconvolution[J]. Geophysics, 1994, 59(2): 215. DOI:10.1190/1.1443583
[7]
FISHER S C, STEWART R R, JOL H M. Ground penetrating radar (GPR) data enhancement using seismic techniques[J]. Journal of Environmental and Engineering Geophysics, 1996, 1(2): 89-96. DOI:10.4133/JEEG1.2.89
[8]
PORSANI M J, URSIN B. Mixed-phase deconvolution of seismic and ground penetrating radar data [C] // SEG Technical Program Expanded Abstracts 1996. Denver: Society of Exploration Geophysicists, 1996: 1603–1606. https://www.researchgate.net/publication/259497221_Mixed-phase_deconvolution_of_seismic_and_ground_penetrating_radar_data
[9]
ARCONE S A, LAWSON D E, DELANEY A J, et al. Ground-penetrating radar reflection profiling of groundwater and bedrock in an area of discontinuous permafrost[J]. Geophysics, 1998, 63(5): 1573-1584. DOI:10.1190/1.1444454
[10]
MORAN M L, GREENFIELD R J, ARCONE S A, et al. Delineation of a complexly dipping temperate glacier bed using short-pulse radar arrays[J]. Journal of Glaciology, 2000, 46(153): 274-286. DOI:10.3189/172756500781832882
[11]
PIRRI M, LOMBARDO P, BUCCIARELLI T. Performance assessment of wide-band pulsed GPR for archaeological applications [C] // Geoscience and Remote Sensing Symposium. Honolulu: IEEE, 2000: 2385–2387. https://www.researchgate.net/publication/251791078_Performance_assessment_of_wide-band_pulsed_GPR_for_archaeological_applications
[12]
XIA J, FRANSEEN E K, MILLER R D, et al. Improving ground-penetrating radar data in sedimentary rocks using deterministic deconvolution[J]. Journal of Applied Geophysics, 2003, 54(1/2): 15-33.
[13]
XIA J, FRANSEEN E K, MILLER R D, et al. Application of deterministic deconvolution of ground-penetrating radar data in a study of carbonate strata[J]. Journal of Applied Geophysics, 2004, 56(3): 213-229. DOI:10.1016/S0926-9851(04)00054-0
[14]
ANNAN A P. GPR Methods for Hydrogeological Studies [M]. Dordrecht: Springer, 2005: 185–213. https://www.researchgate.net/publication/226466848_GPR_Methods_for_Hydrogeological_Studies
[15]
IRVING J D, KNIGHT R J. Removal of wavelet dispersion from ground-penetrating radar data[J]. Geophysics, 2003, 68(3): 960-970. DOI:10.1190/1.1581068
[16]
苏茂鑫, 李术才, 薛翊国, 等. 基于反褶积的探地雷达高分辨率处理方法[J]. 浙江大学学报: 工学版, 2010, 44(6): 1201-1206.
SU Mao-xin, LI Shu-cai, XUE Yi-guo, et al. Ground penetrating radar's high resolution processing method based on deconvolution[J]. Journal of Zhejiang University: Engineering Science, 2010, 44(6): 1201-1206.
[17]
张先武, 高云泽, 方广有. 消除探地雷达数据的子波衰减和频散的反滤波方法[J]. 地球物理学报, 2014, 57(3): 932-938.
ZHANG Xian-wu, GAO Yun-ze, FANG Guang-you. An inverse filtering method for removing the wavelet attenuation and dispersion of ground]enetrating radar data[J]. Chinese Journal of Geophysics, 2014, 57(3): 932-938.
[18]
SCHMELZBACH C, SCHERBAUM F, TRONICKE J, et al. Bayesian frequency-domain blind deconvolution of ground-penetrating radar data[J]. Journal of Applied Geophysics, 2011, 75(4): 615-630. DOI:10.1016/j.jappgeo.2011.08.010
[19]
ECONOMOU N, VAFIDIS A. GPR data time varying deconvolution by kurtosis maximization[J]. Journal of Applied Geophysics, 2012, 81: 117-121. DOI:10.1016/j.jappgeo.2011.09.004
[20]
张志禹, 刘亚丽, 汪文秉. 探地雷达信号的反褶积研究[J]. 西安理工大学学报, 2007, 23(2): 123-126.
ZHANG Zhi-yu, LIU Ya-li, WANG Wen-bing. Research on the deconvolution of ground penetrating radar signal[J]. Journal of Xi'an University of Technology, 2007, 23(2): 123-126. DOI:10.3969/j.issn.1006-4710.2007.02.003
[21]
石刚, 田养军, 杨晓华. 基于波形切除反褶积的探地雷达信号处理效果[J]. 物探与化探, 2012, 36(6): 981-984.
SHI Gang, TIAN Yang-jun, YANG Xiao-hua. The effet of ground penetrating radar signal processing based on waveform removal deconvolution[J]. Geophysical & Geochemical Exploration, 2012, 36(6): 981-984.
[22]
石刚, 田养军, 王亚琼, 等. 基于模糊分形脉冲反褶积的探地雷达信号处理[J]. 长安大学学报: 自然科学版, 2012, 32(5): 79-84.
SHI Gang, TIAN Yang-jun, WANG Ya-Qiong, et al. Signal processing of ground-penetrating radar on the principle of deconvolution based on fuzzy fractal pulse[J]. Journal of Chang'an University: Natural Science Edition, 2012, 32(5): 79-84.
[23]
张建林, 石刚. 路用探地雷达反褶积法应用研究[J]. 中外公路, 2013, 33(4): 85-88.
ZHANG Jian-lin, SHI Gang. Research on application of road ground penetrating radar data's deconvolution[J]. Journal of China & Foreign Highway, 2013, 33(4): 85-88. DOI:10.3969/j.issn.1671-2579.2013.04.022
[24]
石刚, 谢永利, 杨晓华, 等. 探地雷达信号的最小平方反褶积处理效果[J]. 长安大学学报: 自然科学版, 2014, 34(4): 104-108.
SHI Gang, XIE Yong-li, YANG Xiao-hua, et al. Effect of GPR signal processing based on least square deconvolution[J]. Journal of Chang'an University: Natural Science Edition, 2014, 34(4): 104-108.
[25]
ZHAO S, SHANGGUAN P, AL-QADI I L. Application of regularized deconvolution technique for predicting pavement thin layer thicknesses from ground penetrating radar data[J]. NDT & E International, 2015, 73: 1-7.
[26]
CHAHINE K, BALTAZART V, WANG Y, et al. Blind deconvolution via sparsity maximization applied to GPR data[J]. European Journal of Environmental And Civil Engineering, 2011, 15(4SI): 575-586.
[27]
LI L. Sparsity-promoted blind deconvolution of ground-penetrating radar (GPR) data[J]. IEEE Geoscience and Remote Sensing Letters, 2014, 11(8): 1330-1334. DOI:10.1109/LGRS.2013.2292955
[28]
OBRIEN M S, SINCLAIR A N, KRAMER S M. Recovery of a sparse spike time-series by L1 norm deconvolution[J]. IEEE Transactions on Signal Processing, 1994, 42(12): 3353-3365. DOI:10.1109/78.340772
[29]
SACCHI M D. Reweighting strategies in seismic deconvolution[J]. Geophysical Journal International, 1997, 129(3): 651-656. DOI:10.1111/gji.1997.129.issue-3
[30]
VELIS D R. Stochastic sparse-spike deconvolution[J]. Geophysics, 2008, 73(1): R1-R9. DOI:10.1190/1.2790584
[31]
王万里, 杨午阳, 魏新建, 等. 随机稀疏脉冲非线性反褶积[J]. 地球物理学进展, 2014(4): 1780-1784.
WANG Wan-li, YANG Wu-yang, WEI Xin-jian, et al. Stochastic sparse spike nonlinear deconvolution[J]. Progress in Geophysics, 2014(4): 1780-1784.
[32]
梁东辉, 陈生昌. 基于L0范数稀疏约束的地震数据反褶积[J]. 石油物探, 2014, 53(4): 397-403.
LIANG Dong-hui, CHEN Sheng-chang. Deconvolution of seismic data based on L0 norm sparse constraint[J]. Geophysical Prospecting for Petroleum, 2014, 53(4): 397-403. DOI:10.3969/j.issn.1000-1441.2014.04.004
[33]
刘金连, 石战结, 李晨曦. 含泥质碳酸盐岩储层反演技术[J]. 浙江大学学报: 工学版, 2012, 46(2): 372-376.
LIU Jin-lian, SHI Zhan-jie, LI Chen-xi. Inverse technology of clay-bearing carbonate reservoir[J]. Journal of Zhejiang University: Engineering Science, 2012, 46(2): 372-376.
[34]
ZENG Z, XIONG C, JING L, et al. Recursive impedance inversion of ground-penetrating radar data in stochastic media[J]. Applied Geophysics, 2015, 12(4): 615-625. DOI:10.1007/s11770-015-0514-0
[35]
刘钰, 石战结. 基于探地雷达有限带宽阻抗反演的薄层识别技术研究[C]//2016年中国地球科学联合学术年会. 北京: 中国地球物理学会, 2016: 962-963.
LIU Yu, SHI Zhan-jie. Research on thin layer identification based on band-limited impedance inversion [C] // Annual Meeting of Chinese Geoscience Union. Beijing: Chinese Geophysical Society, 2016: 962-963. http://www.researchgate.net/profile/Yu_Liu332/publication/321713685_Thin_layer_identification_based_on_impedance_inversion_of_ground_penetrating_radar_data/links/5a3e70910f7e9ba868934aec/Thin-layer-identification-based-on-impedance-inversion-of-ground-penetrating-radar-data.pdf
[36]
LIU Y, SHI Z, WANG B, et al. GPR impedance inversion for imaging the paleochannel: a case study at Mudu city site in Suzhou, China [C] // SEG Technical Program Expanded Abstracts 2017. Houston: Society of Exploration Geophysicists, 2017: 5115–5118. https://www.onepetro.org/conference-paper/SEG-2017-17659889
[37]
LIU Y, SHI Z, WANG B, et al. GPR impedance inversion for imaging and characterization of buried archaeological remains: a case study at Mudu city cite in Suzhou, China[J]. Journal of Applied Geophysics, 2018, 148: 226-233. DOI:10.1016/j.jappgeo.2017.12.002
[38]
LI J, ZENG Z, LIU C, et al. A study on lunar regolith quantitative random model and lunar penetrating radar parameter inversion[J]. IEEE Geoscience and Remote Sensing Letters, 2017, 14(11): 1953-1957. DOI:10.1109/LGRS.2017.2743618
[39]
张海燕, 李庆忠. 几种常用解析子波的特性分析[J]. 石油地球物理勘探, 2007, 42(6): 651-657.
ZHANG Hai-yan, LI Qing-zhong. Analysis on feature of common analytic wavelets[J]. Oil Geophysical Prospecting, 2007, 42(6): 651-657. DOI:10.3321/j.issn:1000-7210.2007.06.008
[40]
唐俊杰, 郎旭峰. 浙江杭州五代吴越捍海塘遗址[J]. 大众考古, 2015, 2: 14-15.
TANG Jun-jie, LANG Xu-feng. Remains of Wuyue seawall at Five Dynasties in Hangzhou, Zhejiang[J]. Popular Archaeology, 2015, 2: 14-15.
[41]
田钢, 林金鑫, 王帮兵, 等. 探地雷达地面以上物体反射干扰特征模拟和分析[J]. 地球物理学报, 2011, 54(10): 2639-2651.
TIAN Gang, LIN Jin-xin, WANG Bang-bing, et al. Simulation and analysis reflections interference from above surface objects of ground penetrating radar[J]. Chinese Journal of Geophysics, 2011, 54(10): 2639-2651. DOI:10.3969/j.issn.0001-5733.2011.10.021
[42]
SHI Z, HOBBS R W, MOORKAMP M, et al. 3-D cross-gradient joint inversion of seismic refraction and DC resistivity data[J]. Journal of Applied Geophysics, 2017, 141: 54-67. DOI:10.1016/j.jappgeo.2017.04.008