大地电磁倾子有限元数值模拟与稀疏反演研究

来源 :中国科学技术大学 | 被引量 : 0次 | 上传用户:Y644900
下载到本地 , 更方便阅读
声明 : 本文档内容版权归属内容提供方 , 如果您对本文有版权争议 , 可与客服联系进行内容授权或下架
论文部分内容阅读
大地电磁(MT)测深方法是一种将天然交变电磁场作为场源的地球物理勘测方法,一直以来在地球物理勘探以及地壳和上地幔的大尺度构造的基础研究中很重要的应用。传统的MT测深方法通过视电阻率和相位参数来定量研究地下介质的电性构造。然而,倾子资料对于异常体横向边界分辨率高的特点使得倾子数据也成为大地电磁勘探中的一个重要的解释参数。此外,在大地电磁的反演研究中,平滑正则化约束反演是最常见的反演方法,但是这类方法在处理深部区域的异常体时会出现异常假象的问题。然而,使用L1范数的正则化约束条件可以在一定程度上克服这种缺陷。鉴于上述两点,本文主要研究二维和三维大地电磁倾子的数值模拟,并分析讨论各种因素对于倾子响应的影响,同时研究二维大地电磁基于小波变换的稀疏反演,并比较讨论稀疏反演和平滑反演的各自特点。  在二维大地电磁倾子正演数值模拟中,从TE模式的Maxwell方程组出发,使用有限元加权余量法,开展了二维大地电磁倾子参数的正演数值模拟研究。相比较于传统的矩形剖分单元,采用的三角形剖分单元可以模拟各种复杂地形以及在重点研究区域的加密处理可以带来更高精度的求解。首先正演模拟了单个高、低阻两种异常体地电模型的倾子响应,对比和分析了倾子数据和视电阻率、相位数据在大地电磁测深中的各自特点,结果显示了倾子数据具有更高的横向分辨率。在此基础上,由于频率域Z轴倾子航空电磁法(ZTEM)结合了航空电磁法和倾子资料的特点,适用于探测一定深度的异常体的边界位置,于是模拟巷道模型的ZTEM倾子响应,并分析测线高度、异常体埋深、地形起伏、磁导率异常以及噪声误差等因素对于ZTEM倾子响应的影响,进一步加深了对于倾子响应特征的认识。  在三维大地电磁倾子正演数值模拟中,针对ZTEM方法,从Maxwell方程组出发,采用非结构四面体单元剖分,实现了三维ZTEM倾子数据的矢量有限元数值模拟。同时通过对于巷道模型的计算和试验定量分析了测线高度、异常体埋深、地形起伏、磁导率异常以及噪声误差等因素对于ZTEM倾子响应的影响。最后对于单个异常体模型和复合异常体模型,利用倾子参数具有较高的横向分辨率这一特点,实现了倾子数据对于简单和复杂地电模型目标体的水平位置的识别,为研究ZTEM方法和ZTEM反演奠定了基础。  在二维大地电磁稀疏反演中,由于由大地电磁观测数据得到地下电性结构的反演过程是一个非线性的病态问题,于是平滑约束被广泛用来估计正则化的解。在本文中,基于压缩感知算法的思想,提出一种以小波基函数表示稀疏模型的反演方法。反演的目标函数是通过最小化数据拟合差的L2范数和模型罚项的L1范数来实现。通过高低阻异常体地电模型和断层地电模型的反演模拟,验证了二维大地电磁稀疏约束反演的正确性和有效性。
其他文献
地震波速度参数贯穿于整个地震数据采集、处理和解释流程,在地震勘探的每个环节中,速度均起着至关重要的作用。因此精确求取地震波在地下介质中的传播速度,直是地震勘探中值
地球内部物质的成分以及其运移与地幔相变面有着密切的联系,然而了解地球内部的相变面又是一个非常复杂的问题,相变面在地震学上表现为速度间断面,地下速度间断面对地球科学的研
郯庐断裂带是贯穿我国东部NNE走向的一条深大断裂。郯庐断裂带安徽江苏段穿过了大别造山带、苏鲁造山带、长江中下游成矿带及合肥盆地。为了研究该区域地壳速度结构及变形特
碳酸盐岩地层是重要的油气勘探开发领域,世界油气探明储量的50%以上及产量的60%以上都存在于碳酸盐岩储层中。随着油气勘探的逐步深入和勘探水平的不断提升,碳酸盐岩储层研究的