复杂地电模型的探地雷达时域有限差分正演-硬件和射频工程师.PDF

维普资讯 第27卷 第4期 物探化探计算技术 2005年11月 文章编号：1001—1749(2005)04_-0279—05 复杂地 电模型的探地雷达时域有限差分正演 喻振华，冯德山，戴前伟，何继善 (中南大学信息物理工程学院，长沙 410083) 摘 要：目前探地雷达正演模拟，都是针对简单的层状模型、圆状空洞、正方形空洞等单一的规 则模型，而对于地下弯曲的界面或 “V”字形等复杂模型的正演，实现起来较为困难。然而，地下 构造是复杂的，因此，对复杂地 电模型的正演研究显得更为重要。这里以麦克斯韦两个旋度方程 为基本出发点，运用K．S．Yee的空间网格模型理论和时域有限差分法的基本原理，推导出二维 空间的探地雷达正演方程组，通过讨论数值频散关系及其产生原因，推导出符合探地雷达实际传 播规律的理想频散关系。为解决正演模拟时截断边界处的超强反射，采用了Mur超吸收边界务 件，并以自制的探地雷达正演模拟程序为依托，对两个复杂的探地雷达模型进行了正演模拟，得 到了正演剖面图。从其后模拟实例中的两正演剖面图可以看出，边界截断处的干扰波被大大地 减少，消除边界反射后的雷达剖面能更好地指导工作人员进行地质解释 ，从而达到了把雷达的正 演研究深入到复杂地电模型中去的目的，使正演研究更符合实际的地质情况。 关键词：探地雷达；时域有限差分法；正演模拟；理想频散关系 中图分类号：P631．4 文献标识码：A 为反演与解释基础的复杂地 电模型的探地雷达正 0 引 言 演模拟技术，就成了探地雷达理论研究的主要内容 之一 。 ’探地雷达是根据高频脉冲电磁波在地下介质 探地雷达的正演模拟方法主要有两类：即波动 中的传播特性，特别是在不同介电特性介质之间的 方程法和几何射线法。波动方程数值模拟方法实 反射及绕射等波动规律，来探测地下结构和特性的 质上是求解雷达波波动方程，因此，模拟的雷达波 地球物理方法。近年来，随着数字处理技术和电子 场包含了雷达波传播的所有信息，但其计算速度相 技术的飞速发展，探地雷达以其高的探测分辨率、 对于几何射线法要慢。几何射线法也就是射线追 快速无损的特点，已广泛应用于地质找矿、环境监 踪法-3．4J，它将雷达波波动理论简化为射线理论， 测、水文地质调查、工程勘查 及考古 等众多领 主要考虑的是雷达波传播的运动学特征，其主要优 域，成为浅层勘探的有力工具 J。但是雷达的理 点有概念明确、显示直观、运算方便、适应性强；其 论研究与实际应用相比，具有明显的滞后性。如探 缺点是缺少雷达波的动力学信息。对于复杂构造 地雷达的正演模拟还停留在只对常见的简单模型 的正演，波动方程模拟一直占有重要地位，但也存 的合成，资料解释也多以人工肉眼、手动识别为主， 在计算量非常大 的缺点。而 时域有 限差分法 故对于复杂的雷达地电模型的正演及异常体 自动 (FDTD)是近年来发展起来的，基于Maxwell方程 识别的研究成为重中之重。在研究中，解释人员要 的波动方程法，具有直接时域计算、节约存储空间 达到自动而精确地对复杂地电模型进行解析，就必 和计算时间、适合并行计算、计算程序通用性强等 须事先了解地质体的雷达反射剖面特征。所以，作 优点。目前 ，它已在探地雷达正演模拟中取得了良 基金项 目：国家自然科学基金资助03—02) 收稿 日期 ：20o5一o3—28 维普资讯 · 280 · 物探化探计算技术