• / 4
  • 下载费用:1 下载币  

基于探地雷达频谱反演法的薄层识别技术研究

关 键 词:
基于 雷达 频谱 反演 薄层 识别 技术研究
资源描述:
第 32 卷第 11 期 电 子 与 信 息 学 报 010 年 11 月 2010 基于探地雷达频谱反演法的薄层识别技术研究 秦 瑶①②陈 洁①方广有①阴和俊③①(中国科学院电子学研究所 北京 100190) ②(中国科学院研究生院 北京 100039) ③(中国科学院 北京 100864) 摘 要:衡量探地雷达对薄层的识别能力通常有两个指标:最小可识别层厚和反射系数的识别精度。该文提出的频谱反演法,通过对电磁波在多层介质中频谱传输函数的推导,得出回波频谱的正演计算模型,然后采用阻尼最小二乘法进行二参数反演。通过实验,证实了这种方法能够对层厚小于八分之一波长的薄层进行层厚和反射系数的有效识别。 关键词:探地雷达;薄层;层厚;反射系数;频谱反演 中图分类号: 文献标识码: A 文章编号:1009010)11on (00190, ②(00039, ③(00864, to of In is to of of is of is by a 引言层状介质是探地雷达 (测中常见的介质结构,如冰川、公路路基等[1 3]−。薄层是指层厚小于四分之一波长的层。 衡量 薄层的识别能力通常有两个指标:最小可识别层厚和反射系数的识别精度。 薄层的传统识别方法来源于地震勘探中广泛应用的振幅法[4,5],这种方法对于反射系数已知的介质层的探测十分有效。然而, 测中更普遍的情况是层厚及反射系数均未知,这时需要针对层状介质建立严谨的雷达波正演模型,进而采用多参数反演法进行薄层的识别。文献 [3]针对路面结构层介质模型,利用雷达时域回波实现了多参数2009到, 2010回 国家科技支撑计划 (2008助课题 通信作者:秦瑶 反演,但没有涉及最小层厚的分辨率问题。 本文提出的频谱反演法,首先通过对电磁波在多层介质中传输函数的推导, 得出严谨的正演模型,然后利用回波频谱进行反演,不仅能够准确得到薄层反射系数,而且提高了层厚的分辨率,从而增强了 薄层的识别能力。 2 时域回波分析 厚的影响 对于层状介质模型,如果目标层上下界面跟上下地层有较明显的电性差异, 则在该目标层的顶面、底面上都能形成反射波。 薄层厚度的识别能力就是区分上下界面反射波的能力,也就是垂直分辨率[6]。目前,对于 直分辨极限的定义,通常参照 即薄层的分辨率极限为 /4λ 。 设雷达发射信号为高斯脉冲,经收发天线,接收机采集到的信号为二阶高斯微分信号。图 1 为反第 11 期 秦 瑶等: 基于探地雷达频谱反演法的薄层识别技术研究 2761 射系数 时, 3 种不同层厚对应的时域回波。可以看出,当薄层厚度 l 较大时,介质上、下层反射的信号在回波中能够明显区分开,可以直接通过回波信号波峰与波谷的时间差值计算[5]。 随着薄层厚度的减少, 薄层上下界面的回波在时间轴上逐渐重叠,但波形的幅度逐渐增大,这时可以利用回波幅度求取目标层厚,这就是振幅法的探测原理。 射系数的影响 反射系数是表示层界面上反射强度的物理量,电磁波在不同介质的界面上的反射系数取决于两种介质的相对介电常数。两种介质的介电常数差别越大,反射系数的模值 R 也就越大, (全反射 ); 相反, 如果两介质的介电常数非常接近或一样,反射系数模值 R 就很小, (无反射 )。 反射系数对接收回波的影响,可由图 2 的仿真波形直观表述。可以看出,回波幅度随反射系数 由图 2 结合反射系数与介电常数的关系,我们能够得出这样的结论:界面两边介质的介电常数的差别越大,回波幅度越大;相反,如果两介质的介电常数非常接近或一样,回波幅度就很小,很难确定薄层的存在。因此,运用 别薄层要求薄层的介电特性必须与上下地层有较明显的差别。 3 正演模型的建立 状介质的传输函数 实现薄层参数的反演,首先需要推导层状介质的系统传输函数。设计简单的 3 层反射模型如图 3所示[7],源点和场点均在第 1 层中。以下推导过程中不考虑收发天线间的直达波,忽略层之间的多次散射回波。因此,图 3 只给出了电磁波在层状介质之间的主波路径示意图。 以 为例,设电磁波的入射方向与 z 方向夹角为1θ ,波在界面1d− ,2d− 上均会发生反射和折射,电场反射系数 R 和透射系数 T 分别为 1, 1,11, 1, 1,2ii z i z i +++−==+(1) 1,,11, 12, 1,2ii z i +++(2) 其中 i ,1, 1 1i += ;1 别为电磁波在介质 i 和介质 i+1 中的复波数,电磁波垂直入射时,10=;1别为介质i 和介质 i+1 的复介电常数。 区域 1,区域 2,区域 3 中的波可分别表示为 1 111211 1z d ik e+−⎡ ⎤=+⎢ ⎥⎣ ⎦(3) 2 222222 23z d ik e+−⎡ ⎤=+⎢ ⎥⎣ ⎦(4) 333 (5) 其中12R为广义反射系数,它表示第 1 界面1处所有上行波幅值和下行波幅值之比。由于 3 区延伸至无穷远,故式 (5)中只存在下行波。 未知数1A ,2A ,3A 和12R可由界面边界条件获得。区域 2 中的下行波为区域 1 中下行波的透射波和区域 2 中上行波的反射波的迭加;在上交界面1 处,该边界条件可写为 21 11 22 2122213zz d ik d ik d ik e T R AR e−=+ (6) 同时,区域 1 中的上行波是由区域 1 下行波的反射和区域 2上行波的透射形成的。 于是, 在界面1处,有如下约束条件: 11 11 22 212112 121 21223zz d ik d ik d ik e R AR e−=+(7) 由式 (6),2A 可借助于1A 表达,即 ( )()12122112 12221 231k d =−(8) 将式 (8)代入式 (7),得到 ( )()221221212 23 2112 12221 231d d =+−(9) 式 (9)即为 3 层介质时的广义反射系数,它包含了第1 层和第 2 层的反射影响。 由此递推 , 并利用关系式,1 =+ ,可以得到 N 层介质结构的广义反射系数,1为[8]( )()1, 11, 12,1 1,2,,21d i d i e++++−++++−++++=+(10) 对探地雷达测量环境而言,收发天线均置于空气中,因此最简单的测试模型应为 4 层,自上而下依次为:空气层、沙 (土 )层、异常层、沙 (土 )层。因 图 1 不同层厚 l 时的时域回波 ( ) 图 2 不同反射系数时的时域回波 ( /2l λ= ) 图 3 电磁波在 3 层介质中的主波路径示意图 2762 电 子 与 信 息 学 报 第 32 卷 此需要在 3层模型系统函数式 (9)的基础上, 将23R 用广义反射系数23R代替。 令区域 4 与区域 3 的交界面为3 ,由式 (10)得到23R的表达式为 ( )()332332223 3423223 341d d +=+(11) 假设第 4 层延伸至无穷远, 故该区内只存在下行波,满足: 34 34(12) 将式 (11),式 (12)代入式 (10),再利用关系式,1 1,ii i =− ,整理后得 () ()( )()12112212212 23 12 3412 1223 34 12 23 12 3411=+1iD −+−+−(13) 式中 ()12212; ()2332 322k=−= l⋅ ,32ld d=−,即为异常层的厚度。 由于探地雷达大多采用地面耦合的测量方式,因此,收发天线间的直达波与第 1 层地表的反射波12R 在时域上无法区分, 在滤波处理中会随直达波一起滤除,因此,对于地面耦合探地雷达的 4 层介质谱域传输函数应为 () ()( )()12112212212 23 12 341223 34 12 23 12 3411=1iD −+−+−(14) 为简化书写、明确物理意义,令式 (14)中的()12,。 (),l ,异常层反射系数为 R 时 4 层介质的谱域传输函数。同理,可按上述推导方法得到 N 层介质的谱域传输函数。 需要特别说明的是,式 (14)中的参数除了目标层厚度 l 和与薄层反射系数相关的目标层介电常数3ε 外,还有 f ,1ε ,1d ,2ε 和2d 。其中,当测量系统一定时,收发天线相对位置 发射频率 f , 空气介电常数1ε , 地面的位置坐标1d 均为已知量;对于一般的探测介质,2ε 也是可以采用点源反射, 方法预先测量得到的;2d 为异常体上表面的相对位置,可通过对时域波形的分析来确定。 波频谱 回波频谱的正演计算还需要用到雷达的入射波。由于雷达的型号、天线、测试环境都会对雷达入射波产生很大影响,因此入射波不能简单地用一个函数来定义。在每次试验前,需要采用金属板反射法提取雷达入射波[3]。具体实现方法是:将一块表面平整的金属板铺在介质表面,记录接收到的金属板全反射波形;再将天线对向天空,记录无目标时的回波波形;最后将两次记录的波形相减,即可得到雷达入射波。将其进行傅里叶变换,得到入射波频谱 A。因此, 4 层介质模型回波频谱 (),=A (15) 式 (15)即为雷达回波频谱的正演计算公式。 4 二参数反演 首先给出模型中待反演参数的初始假设值为 [ ],,,p=P " (16) 其中 M 为待反演参数个数。 设 Y 为测量回波的频谱值,用泰勒级数将 Y 展开,略去二阶以上项,得 001, 1, 2, , ; 1, 2, ,⎛⎞∂⎟⎜⎟=+ ⋅Δ⎜⎟⎜⎟⎜∂⎝⎠==∑P"" (17) 其中 N 为测量点数。将式 (17)写成矩阵形式: =−=Δ (18) 式中, , ]y=Y " 为测量数据矩阵;12[,, T,] 为假设模型的计算值;1122[, , −ε T,]" 为测量值与计算值之差列阵; Δ=P , , ]p" 为待反演参数增量阵列。 J 为 系统针对层厚和反射系数反演, 阵是二参数表达式如下: 111212⎤∂∂⎢ ⎥⎢ ⎥∂∂⎢ ⎥⎢ ⎥=⎢ ⎥⎢ ⎥⎢ ⎥∂∂⎢ ⎥⎢ ⎥∂∂⎣ ⎦J ## (19) 式 (18)是关于待反演参数修正步长 ΔP 的线性方程组。采用阻尼最小二乘法计算参数增量 ΔP (阻尼最小二乘法的实现过程可参见文献[9,10]),由此得到新的参数列阵Δ,再将这组新的参数代入正演模型,可求出新的预测数据 y ,新的差值列阵 ε ,以及新的 阵;再用这些结果得出新的参数增量由此构成迭代过程直到满足收敛条件为止。 反演迭代误差 e 的定义为:测量频谱值与反演频谱值之间相对误差绝对值的和,即 () ()(),, 1,2,⋅ =∑" (20) 其中 I 为迭代允许次数。当反演误差 ()于给定值,此时的为反演得到的薄层厚度和反射系数。 由于反演结果具有不唯一性,因此在反演过程中需要加上 0l > 的约束条件。在具体应用中,还可以利用待测介质的时域回波相位等先验信息来排除赝解。 第 11 期 秦 瑶等: 基于探地雷达频谱反演法的薄层识别技术研究 2763 5 实验验证 实验采用拉脱维亚 2c 雷达进行探测。在沙坑中 0.3 m 深处埋设厚度为 1.2 的木板。用 得沙4ε = , 。雷达天线中心频率为 磁波在沙中的波长 0.1 。因此,目标层厚可以用波长表示为 = 。图 4 为进行滤波后的 图可以较清楚判断目标层的位置。同时,由图 4 提取含层目标的第 150 道回波数据进行反演计算。首先判断出目标上表面距地表距离约为 0.3 m;再给出一组反演初始值,令两个待反演参数的步长均为 式 (19)构造 阵;之后代入式 (18)得到反演方程;接着用阻尼最小二乘法 (阻尼因子 ),迭代 10 次后输出反演结果,并计算层厚和反射系数的反演误差。 图 5 为金属板反射法得到的发射波形频谱。图6 , 时的反演频谱与测量频谱的对比图,由于实验噪声等影响,二者出现一定偏差。 表 1 为不同反演初始值迭代 10 次的反演结果。可以看出,薄层厚度和反射系数的反演值与真实值的误差均小于 6%,可以满足工程探测需求。 图 4 量图 图 5 发射波形频谱 图 6 计算频谱与测量频谱对比图 表1 不同初始值迭代10次的反演结果 初始值 反演值 误差 (%) 序号 层厚 (单位 :λ ) R 层厚 (单位 : λ ) R 层厚 R 1 3 结论 本文提出采用频谱反演法提高探地雷达对薄层的识别能力,通过对电磁波在分层媒质中频谱传输函数的推导,实现雷达回波频谱的正演计算,进而构造反演模型,采用阻尼最小二乘法实现了二参数反演。通过对实验数据的处理,证实了频谱反演法能够对层厚小于传统 准的薄层进行有效识别,使探地雷达对薄层的识别能力得到了提高。 参 考 文 献 [1] 黄饶 , 陈小宏 , 李景叶 . 层间多次波模拟与影响分析 [J]. 地球物理学进展 , 2009, 24(3): 981[2] of ]. [ 2008. [3] 张蓓 . 路面结构层材料介电特性及其厚度反演分析的系统识别方法 [D]. [博士论文 ], 重庆 : 重庆大学 , 2003. [4] . J]. 1983, 48(1): 12[5] 陈先红 . 基于叠前资料道集上的远近偏移距属性差异的薄层识别 [J]. 油气地球物理 , 2008, 6(4): 51[6] 袁明德 . 浅析探地雷达的分辨率 [J]. 物探与化探 , 2003, 27(1): 28[7] W, J, Y. PM ]. 182009: 1387[8] M]. 1990: 45[9] 张建华,刘振华,仵杰 . 电法测井原理与应用 [M]. 西安 : 西北大学出版社 , 2002: 189[10] 刘国林 . 非线性最小二乘与测量平差 [M]. 北京 : 测绘出版社 , 2002: 12秦 瑶: 女, 1981 年生,博士生,研究方向为雷达信号处理及成像技术研究 . 陈 洁: 男, 1978 年生,博士,助研,研究方向为雷达信号处理 及成像技术研究 . 方广有: 男, 1963 年生,研究员,博士生导师,中科院“百人计划”入选者,研究方向为超宽带电磁学及其工程应用、损耗介质中隐蔽目标的电磁探测、微波成像新技术和新方法和电磁场的数值计算方法等 .
展开阅读全文
  石油文库所有资源均是用户自行上传分享,仅供网友学习交流,未经上传用户书面授权,请勿作他用。
0条评论

还可以输入200字符

暂无评论,赶快抢占沙发吧。

关于本文
本文标题:基于探地雷达频谱反演法的薄层识别技术研究
链接地址:http://www.oilwenku.com/p-63741.html
关于我们 - 网站声明 - 网站地图 - 资源地图 - 友情链接 - 网站客服客服 - 联系我们
copyright@ 2016-2020 石油文库网站版权所有
经营许可证编号:川B2-20120048,ICP备案号:蜀ICP备11026253号-10号
收起
展开