地质雷达正演的方法及实测岩溶技术的应用.docVIP

地质雷达正演的方法及实测岩溶技术的应用 　　 摘要：临沂地区地貌单元为河流冲洪积平原类型，下伏基岩主要为奥陶系石灰岩，该层岩溶裂隙主要在表层发育，存在串珠状裂隙的破碎带，常见充填粘土、碎石。地质雷达可以无损且快速的大范围探测岩溶，保障施工建设顺利进行。为提高探测质量和图像解译精度，正演模拟是必须的，作者以现场常见地质灾害为模型，正演结果和实测相比对，为复杂地质雷达图像的异常甄别提供了依据。 　　关键词：地质雷达 正演 岩溶 粘土 　　Abstract: the linyi area for rivers geomorphic units of flood plains rushed in type, fell under the bedrock mainly for the ordovician limestone, this layer karat fracture main in surface development, existence of fissures beads-strings fracture zone, common filling clay, gravel. Geological radar can nondestructive and fast detection range of karst, security construction smoothly. In order to improve the detection quality and image interpretation accuracy, the forward modeling is necessary, the author on the geological disasters common for model, the modeling results and the measured compared to, for the complex geological radar image abnormal screened to provide the basis. 　　Keywords: geological radar forward modeling karst clay 　　 　　 　　中图分类号：F407.1文献标识码：A 文章编号： 　　临沂地区地貌单元为河流冲洪积平原类型，下伏基岩主要为奥陶系石灰岩，该层岩溶裂隙主要在表层发育，存在串珠状裂隙的破碎带，常见充填粘土、碎石。岩溶发育造成的不良地质，对建筑物地基的影响是显而易见的。因此地基施工前，对基岩进行地质雷达探测，根据雷达图像，识别真异常的位置及范围，并采用合适的工艺对地基进行加固成了该区以后建筑工程必然面临的问题。 　　地质雷达是一种使用高频电磁波探测地下介质分布的非破坏性探测仪器。它通过剖面扫描的方式获得地下剖面的扫描图像。雷达通过在地面上移动的发射天线向地下发射高频电磁波，在地下旅行的电磁波遇到不同的电性界面时，就会发生反射、透射和折射。电介质间的电性差异越大，反射回波能量也越大。 　　反射到地面的电磁波被与发射天线同步移动的接收天线接收后，通过雷达主机精确地记录下反射回波到达的时间、相位、振幅、波长等特征，再通过信号叠加放大、滤波降噪、图像合成等数据加工处理手段，形成地下剖面的扫描图像。通过对雷达图像的判读，便可得到地下目标物的分布位置和状态。 　　一、地球物理特征 　　本次物探场区上部为第四系冲洪积的砂土、粘性土等覆盖层，厚度较小，下伏基岩主要为奥陶系石灰岩。现场上部覆盖层已经挖除，基坑底部入岩2～5m。 　　由现场勘察及开挖情况，可知，岩溶发育比较剧烈，溶孔、溶洞密集成串珠状，其内以充填粘土为主。地下水发育，局部裂隙贯通，成导水通道。 　　根据在本区地质雷达探测经验，灰岩的相对介电常数为6～8，电导率为10-9～2.5×10-2 s/m，电磁波速为0.11～0.12×10-9 m/s；粘土（湿）的相对介电常数为8～12，电导率为10-1～1 s/m，电磁波速为0.06×10-9 m/s；纯水的相对介电常数为81，电导率为10-4～3×10-2 s/m，电磁波速为0.033×10-9 m/s。 　　二、模型 　　本次正演程序为GPRMAX2D及MATLAB7.0，前者基于 FDTD 算法，该算法是K.S.Yee在1966年提出来的[1]。该方法直接求解依赖于时间变量的Maxwell旋度方程组。GPRMAX2D是爱丁堡大学的Dr Antonis Giannopoulos于1996年推出来的一种基于 FDTD算法和PML边界吸收条件的探地雷达正演数值模拟软件，用于探地雷达成像研究[2]。MATLAB7.0用于模型及正演