(劳雷会议)地质雷达处理解释的一点体会杨峰
希尔伯特-黄变换在地质雷达资料处理中的应用
Hilbert Huang Transform and Applications of it in Ground Penetrating Radar Data Analysis and Processing
ABSTRACT
Hilbert-Huang Transform is an algorithm which apply huang transform and hilbert transform to original signal in proper order. It applys in Ground Penetrating Radar data analysis and processing due to its time frequency analysis advantages in order that (1)we can depict the Ground Penetrating Radar record wave characteristics much more finely in time frequency domain to obtain particular and benefit strata information;(2) we can get rid of noise in noise contained Ground Penetrating Radar data to obtain better quanlity,according to the decomposition advantages of huang transform. (3) we can extract natural Ground Penetrating Radar attribute profiles to verify or distinguish the wave impedance and tiny structures,according to the Ground Penetrating Radar data and the work flow for instantaneous attributes. The research method in this thesis is Hilbert-Huang Transform which is a combination of huang transform and hilbert transform.Huang transform can decompose an original signal to a series of intrinsic mode functions which its number is limited.We can get the meaningful instantaneous frequency to display exquisitely in time frequency domain after applying hilbert transform to IMFs. Many numerical algorithms correspond to Hilbert Huang Transform are stimulated and applied in actual Ground Penetrating Radar data processing due to its advantages and characteristics.Huang transform is a superior method because components decomposed by huang transform is much more pure than that of wavelet transform decomposed.The“end effect” phenomena are better suppressed after improving the
c0 GSSI软件RADAN地质雷达资料处理步骤
地质雷达软件RADAN用户手册美国地球物理测量系统公司美国劳雷工业公司2010年10月RADAN处理软件安装安装采集软件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地质雷达探测资料处理流程图数据备份,资料整理,资料处理,资料解释II GSSI处理软件功能模块介绍基本工具打开数据文件，显示雷达数据剖面。
保存数据文件，保存雷达剖面。
选择数据块，选择目标数据剖面。
剪切数据块，切除多余数据剖面。
保存数据块，单独保存雷达数据剖面。
复制剖面图像至剪贴板，地质雷达剖面制作图片。
编辑数据文件头，输入相关参数：标记间隔、扫描数、介电常数、信号位置；编辑标记信息、补充漏打的标记、删除多余标记信息。
线扫描显示方式、以灰阶图或者彩色图形式显示雷达剖面。
一般地质雷达数据处理步骤
一般地质雷达数据处理步骤分界面厚度变化时可用此法，一般不用2）有倾斜地层时可用此法3）使钢筋显示更清楚用此法⑹主要用此法的地方1）测工字钢个数，埋深，形态，间隔2）测空洞3）测钢筋网个数1.反褶积、一维频率滤波（取默认值。
垂直方向上出现一串时（等间隔的多次波）用此）。
Process→Deconvolution；Process→IIR Filter.2.偏移归位Process→Migration，选择偏移类型kirchhoff，调整曲线形态。
3.希尔伯特变化Process→Hilbert Xform，选phase显示瞬态相位信息。
4.添加地面高程信息，并利用高程归一化函数进行处理。
Process→SurfaceNorm。
5.静态校正Process→Static，mode选择manual手动调整方式。
6.文件拼接。
打开Radan软件，选择File→Append files。
7.通道合并，多通道资料对比分析。
打开Radan软件，选择File→Combinechannels。
8.交互式解释View→Interactive，生成*.lay文件。
步骤1）点2）如果从没解释时就选generate new pick file，如果是在原来的基础上对此文件进行解释就选pick file找到lay文件3）选目标体（如钢筋类的，解释后可以看出有多少根）：①在剖面上点右键---target options—new target—双击目标体名字----然后在target parameters里改各个要改的参数②在剖面上点右键---pick options---在pick options里填参数（若拾取工具选block时，在剖面上选一块然后点右键然后加点，）4）选分层①在剖面上点右键----layer options---改layer options里的参数然后确定②在剖面上点右键---pick options---在pick options里填参数（若拾取工具选block时，在剖面上选一块然后点右键然后加点，若当中有空的没有连起来则点右键，插值）5）在剖面上点右键----spreadsheet（表格）6）在剖面上点右键----save changes---current file---保存为lay文件7）用excel打开此lay文件（打开时分割符号选tab键和逗号），打开后去掉头文件然后画图。
c GSSI软件RADAN地质雷达处理步骤
地质雷达软件RADAN用户手册美国地球物理测量系统公司美国劳雷工业公司2010年10月RADAN处理软件安装安装采集软件RADAN66和RADAN5，并且激活采集软件输入软件序列号serialnumber输入处理软件产品ID代码:radan计算获取软件激活码Windows7系统安装radan5安装radan程序，找到setup.exe鼠标右键要求以系统管理员身份运行；RADAN软件第一次运行要以系统管理员身份打开。
Windows7系统调整显示效果选择控制面板->所有控制面板项->显示->更改配色方案->windows经典->高级，对话框如下：选择颜色项目->桌面->颜色->设置红绿蓝资料整理1打测量，布置网格和测线，数据采集2数据拷贝与备份：从地质雷达主机把数据复制在个人电脑上，并利用2种以上存储介质对原始数据进行备份。
3野外记录整理：整理野外记录本（包括各种参数，利用数码相机或者扫描仪对原始纪录扫描拍照，并制作成PDF格式文件便于日后随时查看野外现场原始资料），工作照片，收集的各种第三方资料（设计图纸、设计厚度、第三方检测资料），现场钻孔资料（里程桩号、芯样实物和照片、长度）。
利用钻孔资料反算电磁波传播速度或者材料介电常数。
4数据编辑与初步整理RADAN5资料处理RADAN6资料解释7图片制作8探测报告编写IGSSI地质雷达探测资料处理流程图数据备份,资料整理,资料处理,资料解释IIGSSI处理软件功能模块介绍基本工具打开数据文件，显示雷达数据剖面。
保存数据文件，保存雷达剖面。
选择数据块，选择目标数据剖面。
剪切数据块，切除多余数据剖面。
保存数据块，单独保存雷达数据剖面。
复制剖面图像至剪贴板，地质雷达剖面制作图片。
编辑数据文件头，输入相关参数：标记间隔、扫描数、介电常数、信号位置；编辑标记信息、补充漏打的标记、删除多余标记信息。
线扫描显示方式、以灰阶图或者彩色图形式显示雷达剖面。
地质雷达
探地雷达使用提纲1、适用范围及适用条件2、设计规范及收费标准3、不同地质情况的雷达波形特征1、适用范围及适用条件1.1适用范围：探地雷达法适用于基岩深度、水位深度、软土层厚度与深度，断裂构造等地质工程探查，城市路面塌陷、岩溶塌陷、土洞、滑坡面等地质灾害调查，地下水污染带监测，地基加固效果评价，路面、机场跑道、洞室衬砌检测，堤坝隐患，地下泄露，地下管线及其他埋设物探测，考古探查等。
1.2适用条件：（1）探测目的体与周边介质之间应存在明显介电常数差异，电性稳定，电磁波发射信号明显；（2）目的体在探测深度或距离范围内，其尺寸应满足探测分辨率的要求；（3）测线上天线经过的表面应相对平缓，无障碍，且易于天线移动；（4）测区内不应存在大范围金属构件、无线电发射频源等较强的电磁波干扰，或通过处理无法消除的干扰；（5）不应存在极低阻屏蔽层；（6）单孔或跨孔检测时不得有金属套管；2地质雷达测线测点设计规范及收费标准2.1测线测点设计规范2.1.1工程物探应根据任务要求、探测方法、目的物的规模与埋深等因素综合确定工作比例尺，测网布置应与工作比例尺一致，测网密度应能保证异常的连续、完整和便于追踪；2.1.2布置测线时，测线方向宜避开地形及其它干扰的影响，应垂直于或大角度相交于目的物或已知异常的走向，岩溶、采空区、防空洞等走向多变体的探测宜布设两组相互正交的测线；2.1.3测线长度应保证异常的完整和具有足够的异常背景；2.1.4探测范围内有已知点时，测线应通过或靠近该已知点的布设；2.1.5点测时，测点布设位置、测量应满足资料解释推断的需要；2.1.6工作比例尺确定后，宜参照表1选择测网密度。
表1 工作比例尺与测网密度2.2收费标准地质雷达探测收费参见《工程勘察设计收费标准》第7章——工程物探，收费标准见表2表2 地质雷达收费标准3、波形特征在工程勘察中，常见的不良地质现象有：断层破碎带、裂隙带、富水带、岩溶洞穴、岩性变化带等。
ReflexW读取GSSI劳雷探地雷达dzt文件及信号处理方法
ReflexW读取GSSI劳雷雷达数据和处理方法Dr.Zhang/ 2020.2.211雷达数据1.1雷达数据文件探地雷达的采集原理类似于地震，老一点的数据格式一般跟地震一样都是segy。
后来 GSSI 公司的 TerraSIRch SIR3000地质雷达系统(简称SIR-3000) 在国内普及起来后，格式dzt就成了主流格式（美国SIR系列探地雷达的数据格式）。
现在最新的情况是国内生产单位还是美国最多，其次是加拿大和意大利的。
国内其他常用的GPR数据还包括：DT格式（意大利RIS系列探地雷达的数据格式），DT1格式（加拿大Pulse-Ekko 系列探地雷达的数据格式）以及RD3格式(瑞典MALA系列探地雷达的数据格式)。
1.2 雷达文件dzt数据目前主要是劳雷GSSI SIR 3000和SIR 4000的数据格式。
2REFLEXW软件简介2.1 REFLEXW软件简介本文基于MATLAB开发的GUI界面，主要实现了REFLEXW软件的相应功能。
REFLEXW是地质雷达数据（类地震）数据处理及解释软件，应用于地质雷达的数据处理以及资料解释。
Reflexw软件兼容了世界上大多数雷达的数据格式，在欧美地区，Reflexw已经成为了地质雷达数据处理的标准软件。
随着地质雷达行业的发展，在国内也越来越多的人开始使用Reflexw软件。
软件特点：功能强大，可做多种滤波处理可对个雷达数据进行批量处理导入GPS数据，可绘制测线轨迹、修正地形可显示测线中的标记对不需要的雷达数据可进行删除可做2D剖面处理和3D时间切片处理。
2.2 REFLEXW软件功能分析通过REFLEXW软件，可以实现对探地雷达数据的读入，一维滤波，校正，二维滤波，波形图观察等功能，REFLEXW提供了较为全面的滤波手段，可以将探地雷达图谱处理的更加容易观察在REFLEXW软件中，在显示数据方面，也同样提供了大量的处理方式。
其中，最主要的处理方式为Plot Options。
地质雷达探测异常区域处理方案
地质雷达探测异常区域处理方案地质雷达探测异常区域处理方案一、工程简介本标段为周边小区及商铺密集。
标段工程概况见下表。
地质主要为素填土、粘土、砂层、砾石层、卯石层，粉质粘土、全风化变质砂岩、中风化变质砂岩、微风化变质砂岩，地下水丰富。
二、地质雷达探测1、雷达探测车站周边及隧道线路受影响范围进行了地质雷达探测，通过对地质雷达成果断面图的分析，结合现场实际情况，在地下地层中发现多处异常反射信号，推测存在疑似空洞或土层疏松等地质缺陷。
为进一步核实地层中是否存在空洞及其他地质缺陷，我司根据雷达探测结果对发现的异常反射信号进行现场调查核实，排除非地质缺陷区域、确定地质缺陷区域及除害。
及时对地质缺陷区域进行处理，降低地质灾害发生的机率和危害程度。
2、探测结果通过对地质雷达成果断面图的分析，结合现场实际情况，在地下地层中发现多处异常反射信号，综合已经揭露的地层岩性和地下管线分布特征等情况，推测存在以下地质缺陷。
异常信号里程位置线路开始里程（m）结束里程（m）深度区间（m）异常信号推测结果叠线段24+87424+8800.8m~2土层疏松25+12625+140 1.0m~2疑似空洞或土25+16125+168 1.0m~2疑似空洞或土右线25+20325+2090.8m~2疑似空洞或土左C125+10525+114 3.6m~5疑似空洞或土25+12325+128 1.0m~2土层疏松25+16025+1680.8m~2土层疏松25+19225+196 1.0m~2土层疏松25+24125+245 4.0m~5疑似空洞或土25+27925+285 4.0m~5疑似空洞或土线C225+24025+245 3.5m~4疑似空洞或土C225+28025+284 3.5m~4疑似空洞或土25+13225+135 1.0m~2疑似空洞或土C325+21225+217 1.0m~2土层疏松25+03325+036 1.0m~2土层疏松C425+04525+049 1.0m~2疑似空洞或土25+10425+111 1.0m~2疑似空洞或土25+13825+1440.4m~1疑似空洞或土25+19225+1970.6m~2疑似空洞或土25+23225+238 1.0m~2疑似空洞或土25+23725+242 3.0m~4疑似空洞或土左线C625+31325+3181.0m~2疑似空洞或土异常位置测线起止范围（m ）深度范围（m ）异常信号推断结果18~200.8~2.4疑似空洞41~430.7~2.5疑似空洞B178~810.5~1.0疑似脱空40~450.9~1.4疑似脱空162~166 5.0~7.0软弱土层175~2000.7~1.0疑似脱空217~2220.5~1.5疑似空洞B2240~2670.7~1.2疑似脱空B32~40.6~0.8疑似脱空9~10 3.2~7.0疑似空洞16~200.5~0.9疑似脱空30~350.5~1.3疑似空洞或土12~14 3.0~6.0土层疏松54~570.4~0.7疑似脱空B475~780.7~1.6疑似空洞99~100 1.0~2.0疑似空洞110~1160.7~1.2疑似脱空130~1330.6~1.7疑似空洞B4134~1380.8~1.2疑似脱空130~138 3.0~7.0土层疏松8~150.7~1.6疑似空洞34~380.7~1.6疑似空洞B551~550.9~2.0疑似空洞99~106 2.8~3.2土层疏松138~1480.6~1.7土层疏松严重40~46 1.9~3.5土层疏松严重58~61 2.8~10软弱带B674~810.7~1.2疑似脱空125~128 1.1~1.9疑似空洞13~14 3.0~6.0疑似空洞23~25 1.0~1.9疑似空洞64~70 1.1~2.1土层疏松B786~920.8~1.1疑似脱空103~1100.8~1.8土层疏松严重111~1170.8~1.8土层疏松严重119~1260.6~0.8疑似脱空4.5~6.7 3.1~5.0土层疏松B89.7~11 3.1~5.0土层疏松14.3~16 1.2~2.1疑似空洞根据统计表显示，疑似空洞或土层疏松区域埋深主要在0.5～3米之间，主要分布在南北两侧及区间左线正上方。
对一起原油储罐雷击爆炸火灾的反思
元
二 、事故责任 及处理
大 连 市 人 民政 府 对 此 次
油罐着火情况
水上消防２ １ 年第 ２ ｏ２ 期 ３ ９
消防 队员蹬罐灭 火
幸运 的是 ，大连 港公 安局 消 防支 队 、大 连市消防支队具有较为强大的灭火 作战力 量。 同时 由于经常性的辖区情况熟悉和演练 ，指挥 员对现场情况 比较 了解 ，在火灾发生时 ，立即 要求按照预案启动 固定消防设施 ，同时移动消
系统及 时启 动 ，很快将火扑灭 ，均没有造 成很
大影 响。
详细 的调查 ，对大连港集 团的应 急预案 以及 作 用进行 了全面评估 。
根 据调 查 、现 场勘 查 、技 术检 测 ，调 查
组推 断认 定在 Ｔ ３ 油罐 遭受 直击 雷 ，Ｔ ３ 油 ０ １ ０２ 罐遭受感应 雷后 ，油罐浮顶的一次密封钢板与 罐壁之 间 、二次密封导电靴与罐壁之间 的放 电
成立了 由大连市 安全生产监督 局牵 头 ，大连市 监察 局 、大连 市 总工会 、大连 市港 口与 口岸
消 防ห้องสมุดไป่ตู้ 队增援 。经奋力扑救 ，于 当 日２ 时许将 ０
局 、大连市气象局 、大连市公 安局 、大连市公 安消防局 、大连市金州新 区安监局 、大连市保
３８ 水 上消防 ２ ｌ 年 第 ２期 Ｏ２
加 ，为火灾扑救带来 了巨大 的困难 。
变 电所 ） 电。１ 月２ 日晚 ， 供 １ ２
大连地 区的雷 电灾 害在大连港 油 品码头公 司范围造成 了方 圆 ２ｋ 范 围内大量设 施损坏 。中 ｍ
心变 电所高 压配电的操作 电源
由设 在所 内 的ＭＫ Ｍ智 能高 —
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
地质雷达处理解释的一点体会中国矿业大学（北京校区）杨峰地质雷达应用得到广泛的发展和应用。
本文主要论述近几年来在地质雷达处理、解释等应用方面的其中一部分做一些探讨，希望得到广大专家的指导。
1 水平预测滤波从地质雷达采集的信号来看，存在普遍的水平同相轴信号的干扰。
这些水平干扰信号对不同深度信号影响效果不同，越深的信号影响越大。
这主要是由于高频电磁波在地层传播过程中存在指数形式的衰减和能量扩散等。
因而，在相对变化较小的水平干扰信号的作用下，深部反射信号信噪比明显低于浅部信号的信噪比。
因此对水平干扰信号的去除，就显得非常重要。
如何去除水平信号，其实方法也很多，常用的方法有：（1）背景道去除；（2）窗口滑动平均高通滤波；（3）二维滤波；（4）二维谱的反变化等。
不同的方法都有各自的有点，同时也有各自的缺点。
不同地区、不同仪器、不同勘探目的、不同采集方法可能都有不同的方法选取。
作者在研究去除水平信号过程中，也尝试了不同的数据处理方式，在大量数据试验的基础上，提出水平预测滤波，将信号预测和滤波结合在一起达到去除水平信号的目的。
一、常规不同水平信号去除方法的对比1．原始信号2、背景去除背景去除中背景噪声的计算是对背景道范围进行求均值运算得出的。
这种处理方式对由仪器本身或偶合差异引起的噪声具有较好的效果。
采用背景去噪。
背景去噪是一种常用的处理方法，尤其对数据量大的剖面运算，速度快，操作简单。
但是：在进行背景去噪之前，一定要将剖面上没有的数据（如：天线停滞采集的数据，隧道的蔽塞洞等）删除掉，避免这些信号对综合背景信号的干扰。
3、窗口滑动平均高通滤波滑动平均高通滤波，其实并不是真正意义上的滤波处理，其原理与背景去除相同，无非该方法的背景噪声是随道数窗口移动，对局部信号的突出有更显著的效果。
4、二维滤波二维滤波就是利用F－K域中视速度的不同来提取滤波因子，从而达到压制水平信号的目的。
5、二维谱反变化首先计算出信号的二维谱值，在对二维谱进行编辑，将去除干扰信号的谱值清楚，在通过反变化达到压制干扰信号的目的。
由于二维谱运算量较大，这种方法比较适合短剖面的处理。
二、水平预测滤波水平预测滤波只需要输入一个参数即可，即预测步长。
通过不同的预测步长就可以达到水平信号压制程度不同的目的。
预测步长越小，其水平信号压制能力越强，否则相反。
2 城市区域雷达勘探高压线缆信号的识别与去除在城市区域进行地质雷达勘探，无论屏蔽天线还是非屏蔽天线都会受到高压线缆的干扰，如果不对这些干扰信号识别和去除，可能将干扰信号错误地解释为地层信号，甚至整个剖面都不能解释。
这里讨论的识别和去除并不是针对所有的信号都有效，这里只是抛砖引玉，希望大家能提出更多的方法，扩大我们的思路。
一、高压电缆信号的识别当天线从高压电线杆经过时，在雷达剖面上会形成一个清晰的双曲线异常，如果不做分析，可能会将该异常解释为地下金属管线。
其实识别的方法很简单，就是利用反射双曲线弧度来求出异常双曲线的速度，根据速度大小，来确定异常的来源。
如果所求出速度大于0.2m/ns，则可以初步判断该异常来源高空，可以作为干扰波去除，相反可能来源于地下。
当采用反射双曲线弧度求速度的时候，一定要进行道间隔均匀化处理，即将所有设计的道间隔等距离化。
二、高压电缆干扰信号的去除如何有效地去除高压电缆信号，也是大家关心地问题。
这里提供地去除方法可以适用于任何突变干扰信号地去除。
其步骤如下：1．利用水平预测滤波去除水平信号这个步骤其目的就是把电缆线干扰信号提取出来，压制掉其它任何信号。
2．相位取反把经过水平预测滤波信号地相位取反，即正信号变成负信号，负信号变成正信号，其振幅不改变。
这个步骤处理后的电缆线干扰信号与原来剖面信号相反。
3．剖面叠加将相位取反的剖面与原始剖面进行叠加，由于电缆线干扰信号与原始剖面相位相反，叠加后被很好压制。
3 地质雷达剖面解释雷达剖面解释主要针对显示的剖面进行的，其主要内容有：标记管理及编辑、里程控制、层位追踪拾取、异常物的标定与统计等。
一、标记管理及编辑主要实现如下功能：1．在剖面上直接删除标记2．标间隔数据道的等间隔化3．标的规一化二、里程控制里程控制主要在实际应用中，需要将雷达剖面图与实际里程结合起来。
其主要功能如下：1．输入道间隔参数2．输入起始里程和里程方向三、层位追踪与拾取层位追踪拾取就是利用信号识别方法，从相邻道之间找到同一地层反射的波阻抗界面。
从层位追踪结果可以得到以下不同横坐标的层位数据：1．每一道对应的层位数据2．标记里程对应的层位数据3．道里程对应的层位数据四、异常物体的标记与统计该功能主要提供对剖面确定异常物在屏幕上的直接标定，例如裂隙、空洞等等。
通过标定的结果，可以根据标定的属性直接进行异常物的统计，统计包括如下内容：1．异常物体的范围和大小2．异常物体的数量4 地质雷达地下物体属性定性解释这里主要讨论公路路基基础地质雷达探测结果的分析，主要公路路基探测结果的数据比较充分。
采用地质雷达技术进行公路探测，对压实度差、孔隙度高的地段从雷达剖面来判断很难达到一致的标准，即对同一剖面，不同人得出的结果差异很大。
因此人为干扰因素很大，为了同一判断结果，我们在大量实测数据和钻孔资料的基础上进行公路路基属性解释的研究工作，虽然取得一定结果，但是与现实要求还是相差甚远，本次提出也只是提出一种思路，这种思路可能不是最佳的方法，敬请广大专家提出宝贵意见。
一、技术方法方法1：本次研究主要采用谱分析的方法，采用的谱信号主要有滑动平均谱和振幅谱。
方法2：采用K－L变化，提取信号的特征值及其对应的特征向量。
二、分析结果对比1．背景介质信号对应的谱值分布特征2．不同介质对应的谱值特征3．油面层和水稳层对应的谱值特征4．水稳层离析分析5．水稳层软夹层分析6．水稳层松散分析5 地质雷达超前预报定性解释研究本次主要采用三维显示技术来实现超前预报的定性解释研究，三维技术采用两种方式进行。
一、数据来源第一步：输入掌子面对应测线位置。
第二步：网格化测线数据。
采用三种数据网格化方式：距离反比、最小二乘法和克里格估值算法。
二、地质雷达三维平面解释在开发过程中，根据系统的功能需求把系统分为以下一些模块：总体控制系统：该模块含有三个子模块：图形放缩、参数设置、颜色管理模块，其中，图形放缩模块，是根据需要对图形的放大或缩小1.5倍；参数设置模块是对要显示的图形采取何种方式的参数进行设置；颜色处理：是对图形的颜色进行设置和管理。
另外，该模块还具有图形的存取，和对G3D文件的读取管理。
雷达文件管理：该模块是对雷达采集的原始数据进行管理，同时负责把原始的雷达文件转化为项目管理文件（.prj文件）并对各测线的坐标进行设定和管理。
G3D文件管理模块：负责把项目管理文件转化为G3D文件（其流程图见图），进行参数设置、网格划分和G3D文件的管理。
因为雷达采集的数据为单剖面文件，要想利用雷达文件进行数据拟合，则在拟合之前，把各个雷达文件转化为拟合所需要的G3D文件是至关重要，也是程序的重点内容之一。
把雷达文件转换为G3D文件的前提为：一、雷达文件的在采集时所有参数的设置应该相同，否则用其转换后的文件进行拟合的结果没有任何意义；二、文件的数目大于两个，因为如果文件数目太少其拟合的结果没有实际意义。
在上述两个条件的前提下从雷达文件转化成G3D文件的步骤为：(1) 把单个的雷达文件形成一个项目管理文件，并输入各文件的起始点坐标。
(2) 读取雷达项目管理文件并统计该项目中的测线文件数目。
(3) 对各测线进行检查，看其各测线的参数设置是否一致，如果一致则进行下一步，如果不一致返回出错信息。
(4) 读入测线的数据。
(5) 设置时间间隔，从各项目管理文件的各测线中读取该时刻的采样值，形成需要进行拟合的各时刻G3D文件。
(6) 设置拟合参数，选取拟合算法，对各时刻的文件分别进行拟合，给出拟合结果。
算法模块：模块包含三个子模块：距离反比、最小二乘法和克里格估值算法模块；其中根据试验知，运用克里格算法比其它两种算法效果要好的多。
绘图模块：根据需要绘制各种图形：反射强度的等直线图、颜色充填图和各种预览图、剖面图。
三、三维立体解释本软件主要含有三个显示功能：1、三维模型显示2、高程模型显示3、断面的平面显示本软件分为两个部分：高程数据模型和三维数据模型。
前者有两种显示方式——高程模型显示和断面的平面显示；后者仅有三维模型显示一种显示方式，但可对模型进行剖分。
本软件通过参数面板框调整显示参数，从而调整显示对象的显示方式等。
主要有四个参数面板框：1、3D：选择高程或平面显示、选择显示的组件、参数控制等；2、4D：选择显示的组件、剖分三维模型、剖分参数控制等；3、变换：缩放图形、旋转图形等；4、颜色：选择基本颜色、调节数据颜色等。
其中在高程数据模型下2不可用，在三维数据模型下1不可用。
四、三维显示的内容性质不同的处理内容、不同的分析内容均可实现其三维显示与解释。
6 公路评价系统本次解释主要采用公路评价标准的规范来编写的。
主要针对路面厚度进行评价。
7 感谢。