探地雷达数据采集与解释.doc
- 文档编号:75713
- 上传时间:2022-10-02
- 格式:DOC
- 页数:71
- 大小:25.98MB
探地雷达数据采集与解释.doc
《探地雷达数据采集与解释.doc》由会员分享,可在线阅读,更多相关《探地雷达数据采集与解释.doc(71页珍藏版)》请在冰豆网上搜索。
探地雷达数据采集
以及解释
山东大学岩土中心
第1章.探地雷达简介
1.1工作基本原理
探地雷达(GroundPenetratingRadar,简称GPR)是利用频率介于106~109Hz的无线电波来确定地下介质的一种地球物理探测仪器。
随着微电子技术和信号处理技术的不断发展,探地雷达技术被广泛应用于工程地质勘察、建筑结构调查、公路工程质量检测、地下管线探测等众多领域。
探地雷达的基本原理如图1所示。
发射天线将高频短脉冲电磁波定向送入地下,电磁波在传播过程中遇到存在电性差异的地层或目标体就会发生反射和透射,接收天线收到反射波信号并将其数字化,然后由电脑以反射波波形的形式记录下来。
对所采集的数据进行相应的处理后,可根据反射波的旅行时间、幅度和波形,判断地下目标体的空间位置、结构及其分布。
探地雷达是在对反射波形特性分析的基础上来判断地下目标体的,所以其探测效果主要取决于地下目标体与周围介质的电性差异、电磁波的衰减程度、目标体的埋深以及外部干扰的强弱等。
其中,目标体与介质间的电性差异越大,二者的界面就越清晰,表现在雷达剖面图上就是同相轴不连续。
可以说,目标体与周围介质之间的电性差异是探地雷达探测的基本条件。
图1探地雷达基本原理
1.2电磁波传播特征
探地雷达的电磁脉冲在介质中的传播速度为:
其中c为电磁波在空气中的传播速度,为介质的介电常数,常见介质的介电常数如表1所示。
材质
相对介电常数
材质
相对介电常数
粉质粘土
6
水
81
干砂
3~5
灰岩
4~8
湿砂
20~30
花岗岩
4~7
金属
300
砂岩
6
PVC塑料
3.3
页岩
5~15
混凝土
6.4
淤泥
5~30
空气
1
海水
80
粘土
5~40
表1各种常见介质的介电常数
电磁波脉冲在地质界面上的反射系数为:
根据电磁脉冲的传播规律,在地质界面上如果反射系数为负,则相位与发射脉冲相反,若反射系数为正,则相位与反射脉冲一致。
如图2和图3,可以清除看到反射波相位的变化规律。
图2
图3
1.3雷达的分辨率
对于地质雷达的探测方式.它的分辨率也是一个必须了解的内容.地质雷达的分辨率包括垂直分辨率和水平分辨率。
地质雷达的垂直分辨率主要由地质雷达的波长的二分之一决定。
从波的传播规律可知,识别目标体的尺度一般需大于1/2波长,假设垂育最小可分辨的层的厚度为Dm。
则它的计算式为:
其中,c为电磁波在真空中的传播速度。
可见频率越高,介质的介电常数越大,Dm越小.即垂直可分辨层的厚度越薄,垂直分辨率越高。
地质雷达的水平分辨率是指地质雷达在水平方向上所能分辨的最小异常体的尺寸的能力。
通常用Fresenel来表示,当反射面的埋深为H,发射,接受天线的间距远小于H时,
RF就是水平分辨率的最小尺度。
从计算公式可以看出.当目标体埋深越大,雷达波频率越低,波长越长,则RF越大,水平分辨率越低.反之,水平分辨率越高。
第2章.雷达数据的采集
雷达数据采集的步骤如下:
1.仪器安装调试;
2.现场地质以及其它情况记录;
3.布置测线;
4.参数设置;
5.参数现场校核,如不合格,重新调试参数;
6.数据采集。
下面将各步骤的注意事项说明。
2.1仪器安装调试
①必须在断电的状态下进行安装和拆卸;
②如果现场空气比较潮湿或者有水,注意防潮防水,以免短路导致损坏电路;
2.2现场地质记录
记录现场的地质条件,以免现场的干扰在雷达记录中造成假象,扰乱最终的解释。
2.3布置测线
如果已知探测对象的大体走向,尽量使测线与其走向正交;
2.4参数设置
参数设置是关系到采集数据质量的关键工作,各项参数如下:
Collect菜单(如图4):
注:
①增益:
地质雷达发射的电磁波在介质中传播过程中,在电性(介电常数)分界面上会发生发射,有一部分电磁波继续向下传播,传播过程中电磁波能量会被介质吸收。
随着深度的增加,电磁波能量减弱,信号幅度相应地减小,不利于信号识别和辨认。
为了能更好地识别信号特征,采用增益(gain)函数来提高信号的幅度,使得信号的细微变化更容易显示和识别。
增益菜单有两个选项:
手动/自动Manual/Auto,增益点数Points,GP1、GP2、GP3…GP5。
在采集主机的屏幕的右半部分有一个示波器窗口,除了显示波形之外,还有一条红色的曲线,该曲线就是增益曲线,曲线的转折点点就是增益控制点。
适当的增益函数会提高信号的可视性,但是增益过大会出现削波现象,应该避免削波现象。
②位置确定:
雷达剖面的最顶端并不是第一个反射面(比如在地面探测时,剖面的顶端并不是地面),这主要是因为:
第一,系统延时,即主机给出发射指令到天线开始发射的延迟时间,第二,直达波,即有发射天线直接到接受天线的电磁波。
为了更加精确的定位,应该去除这两方面的干扰。
但是由于确定零线较为困难。
可以在探测时在地面放置一根电缆与测线正交,天线经过电缆时在剖面上会记录下电缆的位置,通过识别电缆就可以确定零线的位置了。
图4collect菜单的注解
Playbacck菜单(如图5)
该菜单用于浏览已经采集的文件,在该菜单下有一个SCAN和PROCESS子菜单,通过子菜单可以可以进行一些参数修改和数据处理,以便得到更好的显示效果。
注:
这些修改和处理只是用于显示,并不会改变原数据。
图5Playbacck菜单的注解
Output菜单(见图6)
该菜单的主要功能就是设置数据的显示参数和数据的导出。
图6Output菜单的注解
System菜单参数设置(见图7)
图7System菜单的注解
2.5参数现场调试
调试技巧:
①如果出现削波(振幅过大)或者振幅过小的现象,需要调节增益设置;
②如果出现雪花现象(特别是深部),需要增加叠加次数;
③如果存在已知深度的目标体并且该目标体可以在雷达剖面上识别的话,可以通过该目标体的深度来反算波速,进而求出介电常数。
这样求出来的介电常数比较接近真实值。
如果没有已知深度的目标体,可以打钻确定一个目标体,然后量测其深度,用同样的方法来反算介电常数。
2.6数据采集
注意事项:
①天线拖动过程中要匀速,并且不要发生跳动现象,跳动现象会造成非常大的干扰。
②打标记要及时,尽量不要重复。
③隧道拱顶或者拱腰部位探测时,可以使用台车(由装载机拖动台车),汽车架子和装载机架子,如图8。
其中台车最安全,汽车架子最不平稳。
台车(由装载机拖动)
汽车架子
装载机架子
图8各种架子
第3章数据分析
好
不
好
成果解释
报告编写
确定处理流程和参数
数据处理
增益处理
IIR滤波
FIR滤波
F-K滤波
算术运算
反卷积
偏移处理
静态校正
希尔伯特变换
……
处理效果
绘制测线布置图
尺寸调整
零线确定
数据分析
频谱分析
希尔伯特变换
数据分析和数据处理的目的就是压制干扰,突出有效信号,提高信噪比,这是进行成果解释的前提,只有进行仔细的处理,才能获得良好的效果。
图9为数据分析和处理的流程图,其中数据分析和数据处理需要进行多次调试对比,才能达到较好的处理效果。
本章将对各个步骤进行详细讲解。
图9雷达数据分析和处理流程
3.1绘制测线布置图和尺寸调整
3.1.1绘制测线布置图
内业工作的第一步就是将现场的测线布置草图绘制成正式文档。
由于现场雷达数据文件是以编号命名的,所以每条测线所对应的文件编号需要记录在文档中,以防混乱。
如下图10,就是雷达探测桥洞的测线布置图,可见在图中标注了测线的方向以及测线所对应的文件编号,应该尽量多的将现场信息反映在布置图上,有助于在解释过程中识别干扰,必要时可以做文字说明。
图10雷达测线布置图
3.1.2雷达数据剖面的尺寸调整
由于在数据采集过程中存在以下的情况:
①采用时间连续采集模式的时候天线走速不均匀,导致标记(等距离标记)之间的道数不一样,甚至差别非常大;
②地表起伏比较大,容易在雷达剖上造成假象;
③由于直达波的存在,使得剖面的最顶部并不是地面的反映,使得深度产生误差。
针对以上的情况,需要进行尺寸调整,各功能模块如下:
剖
面
尺寸
调整
水平尺寸调整
垂直尺寸调整
距离正常化
水平长度的缩放
文件拼接
深度正常化
零线调整
距离正常化
距离正常化,该项功能允许你在标记(必须是等距标记)之间建立等长尺寸,也就是要求在标记之间每单位距离上的尺寸等同,或道数一样。
在没有测量轮的情况下以连续模式采集数时,天线移动速度难以保持恒速,这样就需要利用距离正常化的功能。
联合标记或者距离标记(注:
关于Mark数据库的性质与使用见软件说明书的的第二章)设定后,该项功能就会通过增或删的方式来修正每个标记之间的道数,当然同时也修正了采集速度(天线运行速度)。
参数输入框
注意:
①在运行该功能之前,必须确保标记信息是正确的(无重复标记,无丢失标记,首尾标记都存在,所有的用户标记都已经转化为距离标记或联合标记)。
②Scans/unit,unit/mark都必须在头文件中设置好,以便运行该功能。
③unit/mark是根据测量时的设定来设置的。
④Scans/unit需要用鼠标来清点每个标记之间的道数(注意:
因为在雷达剖面中显示的是已经做了叠加处理的数据,若每个标志之间的道数为3,而叠加次数为16,则实际每个标志之间的道数为二者的乘积48,切记!
!
)
⑤该功能运行之后,原数据中在第一个标志之前的部分已经被cutoff了。
水平尺寸缩放
可以通过水平尺寸缩放中的叠加,去除,添加功能来修改雷达数据。
参数设置对话框
注意每次只能这三项功能中的一个。
叠加(stacking):
使用该功能可以对数据进行简单的滑动平均处理。
该功能就是将所指定的几道平均叠加之后输出一道数据。
去除(skipping):
选定此项,你可以将指定的道去除,比如,你输入参数1,就会每隔一道去除一道数据(因此,数据被压缩到原来的二分之一)
添加(stretching):
选择此项,将拓展水平尺寸。
该功能将计算出每相邻两道的平均值(或指定道数),然后将平均道添加到已有数据中。
拼接文件
有时因为场地条件所限,我们不得不分部采集数据,在后处理中为了将各部分数据连接在一起,就要使用该功能了。
①选择File>AppendFile.
②选择所要连接的各个文件。
③点击Done,完成连接,重新命名加以保存。
表面位置调整(非常重要)
由于系统延时和直达波存在,使得整个剖面的最顶部并不是地面的位置。
确定地面的位置,对于精确的深度定位来说非常关键,但是如何确定,目前还没有定论,下面列举三个常用的方法。
①根据RADAN介绍,百分之九十的情况下把直达波的第一个正峰位置作为地面。
②根据华东院资料,将直达波的第二个波瓣作为地面,如下图
③根据经验,可以在探测时在起始部位放置一根电缆,在后处理时在剖面上识别出该电缆,这样就可以确定地面位置了。
下面以一个探测剖面为例说明一下这三种方法的用法与区别(如图11):
图11零线的确定
如
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 雷达 数据 采集 解释