FDTD算法实现及典型地下结构的GPR剖面特征分析

来源 :同济大学海洋与地球科学学院 同济大学 | 被引量 : 0次 | 上传用户:adzqx2009
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
随着现代电子技术和数据处理技术的发展及应用,探地雷达适用条件已从最初的冰层、盐矿等弱耗介质逐渐扩展到土层、煤层、岩层等有耗介质,其应用范围已经涵盖了考古、矿产资源勘探、岩土勘查、无损检测及工程建筑物结构调查等诸多领域,为经济建设做出了巨大贡献。但是,随着应用范围的拓广,探地雷达勘查的对象也呈现了复杂化,由于缺少理论模型正演结果的对比,限制了解释人员对复杂结构在GPR剖面上表现特征的认识,传统的看图识字式的解释方法对这些复杂对象的GPR剖面解释有时显得无能为力。遵循刘光鼎院士提出的综合地质地球物理研究原则中的正演与反演相结合的思路,对这些复杂勘查对象的GPR反射剖面记录进行数值模拟,分析它们在雷达剖面上显示的特征,为实测GPR剖面的解释提供依据就显得十分必要,这也是提高雷达解决复杂问题能力的探索方向之一。 本文结合工程实践的需要,对当前实际工程勘查中迫切需要探明的复杂地下结构进行归纳总结,建立了包括层状、圆状、阶梯状、断裂破碎带以及组合模型等几大类具有代表性的典型地电模型,使用FDTD方法对建立的几大类典型地电模型的GPR反射剖面记录进行正演合成,并分析、总结了它们在GPR剖面上显示出的特征规律,在雷达剖面图像上显示的波形特征与模型的形态、电性特征之间建立了一定的联系,为实测GPR剖面的解释提供了依据。 为了模拟雷达波反射数据剖面时,前人从依赖时间的Maxwell旋度方程出发,利用二阶精度的中心差分近似把旋度方程中的微分算符直接转化为差分形式,推导出正演模拟的时域有限差分方程。为了保证算法的精确性,这里对FDTD方法解的稳定性与数值频散问题进行了分析。在数值解稳定条件和理想频散关系的基础上设定时间采样间隔和空间剖分网格的尺寸,采用Yee氏网格进行网格的剖分,设置正演模型的空间电磁参数。为了在有限的研究区域获得无限空间的波场模拟,防止雷达波在截断边界引起的超强反射,在截断边界处引入廖氏吸收边界条件,并取得了良好的效果。 利用自行开发的二维GPR正演软件对一系列典型地下结构的雷达剖面进行了正演合成,得到了较好的模拟效果,验证了模拟方法的可行性和软件的正确性。通过对正演合成剖面的分析,指出了不同类型的地下电性结构在GPR反射剖面上表现出的特征和规律。与实际资料的对比,说明正演模拟揭示的GPR剖面上的特征在实际资料中也是存在的,这就为以后工程实践中类似结构的GPR探测剖面的解释提供了依据。 通过模型试验发现,在利用FDTD方法进行GPR正演模拟的过程中,为了得到质量较好的GPR正演剖面,重点是合适模型的抽象和合适边界条件的加入。文中通过设置断裂破碎带中介质的介电常数在一定范围内随机变化,使得正演模拟得到的GPR剖面较好地揭示了断裂破碎带在实测雷达剖面上显示的特征。
其他文献
在二维地震偏移成像中,染色算法被证实可以有效地提高弱照明区域的信噪比。染色算法应用于逆时偏移过程中,与目标成像区域相关的染色波场和常规的震源波场同步正向传播。分别
本文主要分两个部分。第一部分主要讲述了两个内容:一是大气边界层研究的意义,介绍了目前用于边界层研究的各种方法以及用激光雷达研究边界层的状况;二是研制机载激光雷达的意义以及目前国内外机载激光雷达的研究状况。第二部分为本人所做的研究工作,包含四块内容。(一)机载大气探测激光雷达的研制通过大气回波信号和信噪比的模拟计算,确定了机载大气探测激光雷达的主要技术参数;在详细叙述和探讨了机载大气探测激光雷达总体
大地电磁法是一种以天然电磁场为场源的电磁勘探方法,成本低,工作简便,不受高阻层屏蔽,尤其对低阻层较敏感,勘探频带宽。因此,大地电磁法在研究深部构造中有着独特的优势。大
反褶积作为地震资料处理的三大关键技术之一,在整个处理过程中起着十分重要的作用。反褶积结果与所采用的方法及其假设条件息息相关,通常我们都是假设地震子波为最小相位的,而在