c0-GSSI软件RADAN地质雷达资料处理步骤20120725

. . . 地质雷达软件RADAN用户手册美国地球物理测量系统公司美国劳雷工业公司2010年10月69 / 69RADAN处理软件安装安装采集软件RADAN66和RADAN5，并且激活采集软件输入软件序列号serial number输入处理软件产品ID代码:radan计算获取软件激活码Windows 7 系统安装radan 5安装radan程序，找到setup.exe鼠标右键要求以系统管理员身份运行；RADAN软件第一次运行要以系统管理员身份打开。Windows 7 系统调整显示效果选择控制面板-所有控制面板项-显示-更改配色方案-windows经典-高级，对话框如下：选择颜色项目-桌面-颜色-设置红绿蓝资料整理1打测量，布置网格和测线，数据采集2数据拷贝与备份：从地质雷达主机把数据复制在个人电脑上，并利用2种以上存储介质对原始数据进行备份。3野外记录整理：整理野外记录本（包括各种参数，利用数码相机或者扫描仪对原始纪录扫描拍照，并制作成PDF格式文件便于日后随时查看野外现场原始资料），工作照片，收集的各种第三方资料（设计图纸、设计厚度、第三方检测资料），现场钻孔资料（里程桩号、芯样实物和照片、长度）。利用钻孔资料反算电磁波传播速度或者材料介电常数。4数据编辑与初步整理RADAN5资料处理RADAN6资料解释7图片制作8探测报告编写I GSSI地质雷达探测资料处理流程图数据备份,资料整理,资料处理,资料解释3文件编辑-剪切4剖面方向调整5距离归一化6添加起始里程桩号7剖面水平拉伸、压缩8调整地面时间零点10时间深度转换2读入数据文件(*.dzt)1打开RADAN软件11增益调整12叠加去噪13背景去除14一维垂直滤波15反褶积16偏移归位17希尔伯特变换静态校正18高程修正频谱分析速度分析19剖面追加21交互式解释II GSSI处理软件功能模块介绍基本工具打开数据文件，显示雷达数据剖面。保存数据文件，保存雷达剖面。选择数据块，选择目标数据剖面。剪切数据块，切除多余数据剖面。保存数据块，单独保存雷达数据剖面。复制剖面图像至剪贴板，地质雷达剖面制作图片。编辑数据文件头，输入相关参数：标记间隔、扫描数、介电常数、信号位置；编辑标记信息、补充漏打的标记、删除多余标记信息。线扫描显示方式、以灰阶图或者彩色图形式显示雷达剖面。波形加变面积显示方式、以波列图形式显示雷达剖面。线扫描+波形显示方式波形显示方式，以波形方式显示雷达剖面，做频谱分析，显示雷达工作频段。三维数据体显示时间切片显示方式多通道显示方式、多个雷达剖面做比对处理和信息显示。交互式解释、进行地质解释，绘制地质剖面图、给出钢筋位置。调整显示参数，设置各种显示参数，保存、调用显示参数文件，设置长标记。打印雷达剖面关于RADAN程序版本信息命令参数按钮剖面向左滑动显示剖面向右滑动显示执行命令暂停停止并保存数据宏文件命令按钮新建宏文件打开宏文件编辑宏文件参数关闭宏文件软件功能按钮剖面水平拉伸与压缩距离归一化，给出雷达剖面距离信息。高程修正，消除地形剧烈起伏对地质雷达剖面的影响。静态校正信号位置调整，时间置0.算术运算，调整雷达数据大小。增益调整（自动增益、指数增益）、信号放大或缩小增益恢复，去除现场数据采集时的增益影响，恢复雷达数据到原始信息。IIR滤波，水平滤波、垂直滤波FIR滤波，水平滤波，消除随机噪音、水平干扰，和垂直滤波提高频率信息。反褶积运算，作用有：1压缩子波提高垂直分辨率；2消除多次波。偏移归位处理，提高地质雷达剖面水平分辨率希尔伯特变换二维空间滤波速度分析，计算电磁波在地下传播速度自动提取局部峰值信号多通道雷达资料运算GSSI软件RADAN地质雷达资料处理步骤1. 打开软件RADAN，选择文件夹ViewCustomizeDirectories. 2. 编辑文件属性，去除只读属性。打开文件FileOpen(*.dzt)。3. 扫描信息预编辑：利用图标EditSelect, 选择一段扫描剖面，切除多余扫描信息Cut，或者保存特定扫描剖面 Save。4. 文件测量方向掉转。打开文件，选择FileSave As ，打勾。5. 添加距离信息。测量轮测量直接获取距离概念。连续测量方式加距离需要三步A) 编辑文件头的距离信息EditFile Header, 扫描/米scans/m, 米/标记m/mark，B)编辑用户标记，C)并利用距离归一化函数进行处理,ProcessDistance Norm。6. 添加里程信息.A)EditFile Header 3D optionX start输入里程起点坐标. 7. 水平刻度调整。 ProcessHorizontal scale.叠加stacking、抽道skipping、加密stretching。8. 确定地面反射波信号位置EditFile Headerposition(ns)。9. 调整信号延时信息，找到地面ProcessCorrect Positiondelta pos (ns).10. 设置和修改介电常数，计算深度信息EditFile HeaderDielConstant。11. 信号振幅自动增益调整 ProcessRange Gain。选择参数：增益类型Gain Type为自动增益automatic，增益点数number of points为5。12. 水平相关分析，消除雪花噪音干扰。 ProcessFIR Filter，选择水平叠加Horizontal Stacking(scans)为5。13. 背景去除，显示构造特征。 ProcessFIR Filter，选择背景去除Background Removal(scans)为1023。14. 一维频率滤波 ProcessIIR Filter.15. 反褶积、一维频率滤波。 ProcessDeconvolution； ProcessIIR Filter.16. 偏移归位 ProcessMigration，选择偏移类型kirchhoff，调整曲线形态。17. 希尔伯特变化 ProcessHilbert Xform，选phase显示瞬态相位信息。18. 添加地面高程信息，并利用高程归一化函数进行处理。ProcessSurface Norm。19. 静态校正 ProcessStatic，mode选择manual手动调整方式。20. 文件拼接。打开Radan软件，选择FileAppend files。21. 通道合并，多通道资料对比分析。打开Radan软件，选择FileCombine channels。22. 交互式解释ViewInteractive，生成*.lay文件。23. 绘制地质剖面图.利用电子表格Excel或者Surfer 8、CAD软件绘制地质图件。III GSSI软件RADAN地质雷达资料处理 4调整测线方向 5 加距离 6起点桩号 8地面对应的时间 10深度1. 打开软件RADAN，选择文件夹.视图自定义文件目录. 2. 文件显示。文件打开(*.dzt)。文件显示，换颜色。波形显示选择第二个图标，采用波形方式显示数据剖面 选择显示器图标，设置波形参数，比例为间隔2倍，如32-16。标准为正，大小为0图像局部显示与放大选择显示器图标，弹出显示参数设置窗口。设置开始采样和完毕采样参数，仅仅显示地质雷达剖面在竖直方向的局部信息，相当于局部放大。垂向局部显示剖面原始图局部放大图参数设置（start sample 700, end sample 1500）3. 扫描信息预编辑利用图标编辑选择, 选择一段扫描剖面，切除多余扫描信息删除，或者保存特定扫描剖面保存。剪切处理须知：有标记的位置不能切掉，原因是标记代表了距离概念。切除的是同一个位置多采集的扫描信息。4. 文件测量方向掉转打开文件，文件另存为-方向反转，打勾。目的是调整测线方向，重新排列测线方向，按照习惯从小桩号到大桩号。测线方向方向反转操作命令上图原始图 以下图是反转后的图像5. 添加距离信息。测量轮测量直接获取距离概念。连续测量方式加距离需要三步A) 编辑文件头的距离信息编辑文件头, 扫描/米scans/m, 米/标记m/mark，B)编辑用户标记，C)处理距离归一。视图-自定义-数据库-临时数据库;开启标记编辑功能图5-0图5-A图5-b图5-c原始标记保存标记标记类型转换距离归一化6. 添加里程信息.编辑文件头三维选项X起点，输入里程起点坐标。原始数据添加里程数据RADAN5006添加里程桩号Y起点坐标，表示测线坐标RADAN6.0添加里程桩号调整起始里程信息7. 重采样：垂直方向和水平方向。7-1垂直方向重采样。菜单:编辑-文件头-采样/扫描,512修改为1024,2048.7-2水平刻度调整处理水平缩放.叠加、抽道、加密。打开软件（不打开数据文件），视图-自定义-数据库关闭标记信息数据库对话框水平刻度调整地质雷达测量原始剖面8. 确定地面反射波信号位置编辑文件头信号位置(纳秒)，如 -2.5。建议：数据采集时，主机中所用的参数保持固定：比如position/offset信号位置和延时。9. 调整信号延时信息找地面处理信号位置调整信号移动(ns)。10. 深度计算设置和修改介电常数，计算深度信息编辑文件头介电常数。介电常数 时间 深度C=0.3m/nst=t目标-t地面=t目标- 0=t目标11. 信号振幅调整野外现场采集数据前对信号进行放大调整。Gain。SIR-3000增益设置SIR-20增益设置地质雷达剖面记录的文件经过了增益放大，现场调整先自动调整，切换至手动模式，人工调整增益（信号放大倍数），手动锁定增益值；进行数据采集，即记录的不是最原始的数据，而是经过增益放大处理以后的数据。若现场采集时的信号调整足够清晰，则利用后处理软件直接显示数据。调整显示增益即可。自动增益调整处理增益调整自动增益，增益点数为5。放大倍数，值一般选择2-5。自动增益调整自动增益处理比对指数增益参数设置：手动设置增益点数，调整增益值大小原始数据局部指数增益原始数据压制局部干扰-近电磁波场干扰利用指数增益压制电磁波近场的强干扰。12. 水平相关分析，消除雪花噪音干扰。现场数据采集。SIR-3000主机叠加3-5次SIR-20数据采集设置对话框处理FIR滤波水平叠加(扫描)5。原始数据5次平滑数据13. 背景去除，显示构造特征。针对的目标是局部不规则物体。对于水平层状结构地下物体探测（比如地下分界面），不适用背景去除法。处理FIR滤波背景去除(扫描)为1023。滤波参数选择原始数据处理结果14. 一维频率滤波，消除低频电磁干扰 处理IIR滤波。先做频谱分析,再做频率高通滤波。仪器工作过程：利用地质雷达探测地下结果或者地质构造，在测量中，发射命令由主机发出，经由电缆传给发射机，发射天线接收到命令以后开始发射电磁波至地下空间；电磁波在地下空间中传播遇到地下界面发射发射，反射回来的电磁波被接收天线捕捉到；接收命令由主机发出，经由电缆传输给接收机，接收机接到命令以后就把接收到的信号反馈给主机。电磁波传播过程，5次用到滤波：a)发射天线发射一定频率带宽的电磁波至地下空间；b)电磁波在地下传播过程中，有一部分能量会耗散掉，高频能量容易损失，从而留下中频和低频的电磁波继续传播，这样反射回来的电磁波以中低频为主；c)接收天线（一定尺寸的）也只能接收一定带宽的电磁波;d)在采集软件中，我们对接收到的信号进行滤波处理；e)在后处理软件中，我们也可以再利用处理软件进行频谱分析和带通滤波。原始数据高通滤波(参数160兆赫兹)频谱图与地质雷达原始记录曲线频谱图-波形图-线扫描图低频信号垂直滤波(MHz) 高通40原始数据与高通滤波数据高通滤波一般针对100兆天线做处理。15. 反褶积、频率滤波尖脉冲反褶积原理：频率越高，周期越短，电磁波垂直分辨率越高。处理反褶积算子长度； 处理IIR滤波。尖脉冲反褶积参数选择示意图反褶积参数设置水泥板测试原始剖面水泥板测试反褶积剖面900兆天线原始数据900兆天线反褶积处理数据多次波反射信号图谱16. 偏移归位天线结构图：双天线收发距示意Tx-Rx 900MHz 0.1524m 400MHz 0.16m电磁波反射示意图红点:目标体原理：双曲线形成机理浅层反射目标体收发距 深层目标体目标体反射形成双曲线形态。天线之间存在收发距，电磁波传播过程中需要考虑收发距的影响。 处理偏移归位偏移类型：克希霍夫，调整曲线形态。双曲线拟合 Er =(0.3/v).2根据电磁波传播速度大小，反算出材料的介电常数信息。双曲线拟合界面偏移对话框速度9.5117厘米/纳秒,代表周围材料中电磁波传播速度信息双曲线拟合介电常数计算电磁波干扰信号排除方法：偏移归为电磁波传播速度为0.3米/纳秒，表示双曲线波形为空气中物体的反射信号投影到地质雷达剖面中。干扰排除地下管线与高压线/空中建筑物比对图地下管线与高压线/空中建筑物比对图管线与高压线模拟-速度拟合图建议在三维测量数据处理中采用偏移归位处理，以提高横向分辨率。偏移前数据剖面偏移后数据剖面左：水平剖面1，右：剖面2上：剖面1，下：剖面217. 希尔伯特变换处理希尔伯特变换。18. 添加地面高程信息。处理高程修正。利用版本5和6添加高程信息的步骤一样，操作略有不同。对于版本5操作如下2步。添加高程数值与高程修正。第一步，点图标打开”标记编辑”对话框，为每个扫描点输入对应的高度信息(Z)，按完成。第二步，点图标，弹出高程修正图标，如以下图，按确定并保存结果18-1。如果地形相对比较平缓，选择垂直比列尺1:1。在野外现场采集数据时，建议把记录长度range略微放大，在数据处理时，利用RADAN软件把数据剖面下移一部分（详见第9步），再做地形校正18-2。18-1高程修正对比图18-2修正图（改善）地形校正剖面图19. 文件拼接。打开Radan软件，文件文件添加。文件拼接所需要的工作。要求在垂直方向是同一个频率的天线测试的雷达剖面才能拼接，所有的测量参数相同，如采样率samples，信号起始位置position，以与记录长度range，甚至是增益。而在水平方向上，2个文件相连接或者测线有重合部分，水平方向上的采样信息要基本一致，比如都是时间模式测量或者都是距离测量的剖面。打开软件RADAN，而不打开任何数据文件，直接选择菜单”文件-添加文件”，弹出对话框。依次选择对应的文件。按照测量的先后顺序依次选择对应的雷达剖面，并按完成按钮，生成雷达剖面并且保存一个文件名，就是拼接后的雷达剖面。对拼接后的文件做进一步处理。20. 通道合并通道合并的前提：2个文件的参数基本一样，如采样点数，剖面长度与扫描数相同。多通道资料对比分析。打开Radan软件，选择菜单”文件合并通道”。选择文件选择添加按钮或者双击文件名，依次选择并添加文件，形成一个文件列表。按”完成”按钮，软件提示保存多通道文件，如以下图。保存多通道文件双通道雷达图片选择菜单按钮多通道计算按钮多通道计算窗口通道1减去通道2的数据，选择命令sub，按OK保存文件如下。通道分离：选择radan-文件-另存为通道分离，并且保存多个文件，如file012a.dzt, file012b.dzt,file012c.dzt,file012d.dzt等等。对于相减运算后的数据，如file012a.dzt，需要做相位取反运算才能和原来文件的相位保持一致。选择菜单按钮，弹出下面对话框，选择算术函数”取反”命令，并且选择确定保存一个新文件即可。21. 交互式解释生成*.lay文本文件视图交互式解释。原始数据交互式解释参数选择界面1界面2界面1和2电子表格隧道探测图像解释波形分析22图片制作，截取并输出地质雷达剖面。