基于非结构有限元法的三维大地电磁测深法各向异性正反演研究

论文摘要

大地电磁测深法（MT）是一种基于天然源的地球物理电磁测深法,它是地球深部构造研究的最重要的方法之一。当前三维大地电磁测深法在资源勘探、地质调查和浅地表至上地幔深度地质构造研究中获得了广泛应用。传统的三维大地电磁测深法正反演算法是基于各向同性理论,且采用规则的网格进行区域离散难以适应地形和复杂的构造。然而,越来越多的大地电磁测深法数据表明地下介质的电各向异性是广泛存在的,尤其是在深部构造中。采用各向同性理论进行反演和解释会造成解释的困难和错误。只有基于各向异性理论的反演和解释才能给出合理的结果。实际勘探研究中,地形和地下构造非常复杂,传统基于规则网格的正反演方法模拟和反演复杂地形和地下构造的能力有限,难以满足实际勘探需要,限制了其实用性。当前在三维复杂地形和地质条件下,大地电磁测深法响应正演模拟、大地电磁测深法各向异性响应特征识别及大地电磁测深法各向异性数据高效反演方法和处理策略都是亟待研究和解决的问题。为了解决以上问题,本文开展了基于非结构有限元的三维大地电磁法各向异性正反演研究。三维大地电磁测深法各向异性正演是三维各向异性反演的基础。高效的正演算法可以提高三维大地电磁测深法的计算速度和精度,是构造高效三维反演最重要的因素。在三维大地电磁测深法模拟中,为了得到高精度的数值解,需要对高度复杂的地质构造和界面进行精确建模。然而,传统的规则网格不能实现高精度地模拟任意复杂地质体和界面。最有效的解决办法就是使用非结构网格。非结构四面体网格可以灵活地构建任意的地形和地质构造,且可以灵活地实现局部优化。非结构四面体网格单元相比结构网格可以使用较少的单元拟合任意复杂的地形和地下构造。有限元法是一种适应非结构四面体网格的高效数值算法。本文采用非结构四面体单元有限元法系统地开发了各向异性介质条件下的三维大地电磁正演算法,该算法可以用于直接进行带地形的三维大地电磁测深法数据反演,克服了结构网格计算精度有限和难以拟合地形和地下复杂异常体的缺陷。在任意各向异性数值模拟部分,本文首先分析了在水平地形条件下不同各向异性参数对MT响应的影响特征,总结了根据视电阻率识别各向异性信息的方法。然后,针对在地形存在条件下根据常规的视电阻率响应难以识别地下介质各向异性信息的问题,本文根据各向异性响应的极性特征,提出使用极坐标下的视电阻率来识别地下介质的各向异性信息的方法,并证明了其在复杂地形存在的条件下仍然具有适用性。为了克服人工网格剖分的困难,尤其是地形和复杂构造和各向异性同时存在的地质情况,本文开发了基于加权后验误差的面向目标网格自适应算法,该方法能得到高质量的网格,从而大大提高有限元正演的计算精度。同时,本文分析了针对各向异性介质的网格自适应特征,进而研究了主轴电阻率对网格的自适应优化的影响规律,结果表明水平方向最小主轴电阻率决定了网格的自适应的疏密程度。高内存需求和低计算效率是当前三维电磁反演实际应用中面临的主要难题。针对三维大地电磁测深法反演问题,本文开发了基于有限内存拟牛顿法（L-BFGS）的反演算法。该算法在反演迭代过程中只需要求取目标函数的梯度,避免了求取海森矩阵,可以有效地降低内存需求。本文研究了基于非结构网格和各向异性条件下的目标函数构建及其梯度的求取方法,分别实现了各向同性和各向异性条件下的三维大地电测深法反演算法。通过带地形的和实际地质背景的三维大地电磁测深反演算例表明采用非结构网格的有限元法进行反演,可以避免复杂的地形校正过程,实现直接带地形反演,且可以建立有效的先验模型。利用本文算法对美国EarthScope项目的USArray的部分长周期数据进行各向同性反演,并和基于规则网格算法的WSINV3D程序的三维反演结果进行对比,验证了本文反演算法的有效性和适用性,同时表明本文基于非结构网格的算法可以恢复更多的细节信息。鉴于地下介质存在各向异性时,各向同性反演不能恢复真实模型的分布规律和电性特征,会造成反演解释的困难和错误。在三维大地电磁各向异性算法的基础上,本文提出了一种对三维大地电磁测深法各向异性数据进行反演解释的策略。首先,利用提出的极坐标下各向异性识别方法,定性地分析和识别地下介质各向异性特征。然后,若确定数据包含各向异性信息,采用三维各向异性算法进行反演,否则采用各向同性算法进行反演解释。其次,根据反演出的各向异性比值分布,确定各向异性构造的空间分布。最后,结合地质背景进行三维各向异性解释。通过带地形和不带地形的各向异性反演数值算例表明,本文提出的各向异性识别方法对复杂地形的各向异性数据是有效的,可以用来指导三维大地电磁测深法各向异性的识别和解释,避免了采用传统的各向同性理论容易造成的解释困难和错误;基于非结构网格和L-BFGS的反演算法,适合带地形条件下的各向异性反演;在三维大地电磁测深法各向异性反演中,ρx和ρy为可反演恢复参数,而ρz不可恢复参数。最后,本文利用三维各向异性反演算法和提出各向异性处理策略对Coprod2大地电测深法数据进行各向异性反演解释,进一步验证了各向异性反演算法和策略的有效性。数据反演结果建立了地下深部各向异性构造的模型,解释了NACP和TOBE的位置,并发现了新的各向异性体。对USArray实测数据三维各向异性研究,证明了Cascadia俯冲带上板块运动的应力方向为各向异性电阻率的良导方向,并建立了Cascadia俯冲带的动力学模型,对研究地壳深部构造应力及运动的研究及各向异性的成因具有一定的科学意义。本文基于非结构有限元法和有限内存拟牛顿法开展了三维大地电磁测深法各向异性正反演研究。本文实现了高效的任意各向异性介质中三维大地电磁测深法面向目标自适应有限元法正演模拟;提出了地形条件下三维大地电磁各向异性响应特征识别方法;建立了适用于实测数据的三维大地电磁测深法各向异性反演方法及处理策略。本文为三维大地电磁测深法各向异性介质正反演研究奠定了理论基础,为研究地下深部介质的电各向异性提供了技术支持。研究成果能够提高三维大地电磁测深法数据处理和解释的水平,对地球深部构造的研究具有一定的理论和实际意义。

论文目录

• 摘要
• Abstract
• 第1章 绪论
•   1.1 研究目的及意义
•   1.2 三维大地电磁测深法各向同性正反演技术发展现状
•   1.3 大地电磁法测深法各向异性正反演研究现状
•   1.4 自适应有限元法的发展现状
•   1.5 本文主要研究内容
•     1.5.1 各向异性介质下的大地电磁场的基本理论
•     1.5.2 三维大地电磁测深法各向异性正演模拟及识别
•     1.5.3 三维大地电测测深法面向目标网格自适应
•     1.5.4 三维大地电磁测深法各向同性反演算法
•     1.5.5 三维大地电磁测深法各向异性反演算法
•   1.6 本文创新点
• 第2章 各向异性介质中的大地电磁场基本理论
•   2.1 大地电磁场的场源
•   2.2 地下介质中电磁感应基本原理
•   2.3 电各向异性基本理论
•     2.3.1 电各向异性的起源
•     2.3.2 电导率张量
•     2.3.3 欧拉旋转
•   2.4 大地电磁测深法一维解析解
•     2.4.1 各向同性介质的解析解
•     2.4.2 任意各向异性介质的解析解
•   2.5 本章小结
• 第3章 三维大地电磁测深法各向异性正演算法
•   3.1 有限元法基本理论
•     3.1.1 边值问题及其求解
•     3.1.2 非结构网格离散
•     3.1.3 有限元法基函数
•   3.2 正演模拟算法
•     3.2.1 有限元控制方程的推导
•     3.2.2 单元矩阵计算
•     3.2.3 总体矩阵合成及压缩存储
•     3.2.4 强加边界条件
•   3.3 大地电磁测深法向异性响应函数
•     3.3.1 阻抗张量及视电阻率和相位
•     3.3.2 倾子响应
•   3.4 大地电磁测深法各向异性响应特征
•     3.4.1 精度验证
•     3.4.2 各向同性模型
•     3.4.3 主轴电阻率绕x轴旋转
•     3.4.4 主轴电阻率绕y轴旋转
•     3.4.5 主轴电阻率绕z轴旋转
•   3.5 大地电磁测深法各向异性识别
•     3.5.1 地形条件下的各向异性响应特征
•     3.5.2 极坐标下的视电阻率响应特征
•   3.6 本章小结
• 第4章 面向目标网格自适应技术
•   4.1 面向目标自适应算法
•     4.1.1 后验误差估计
•     4.1.2 加权系数
•     4.1.3 加权后验误差
•   4.2 网格自适应策略及流程
•   4.3 网格自适应结果分析
•     4.3.1 层状介质自适应结果
•     4.3.2 三维介质网格自适应结果
•     4.3.3 各向异性介质对自适应的影响
•   4.4 本章小结
• 第5章 三维大地电磁测深法各向同性反演研究
•   5.1 反演基本理论
•     5.1.1 反演问题的数学描述
•     5.1.2 反演的目标函数及其梯度
•     5.1.3 终止条件
•   5.2 三维最优化算法
•     5.2.1 牛顿法与高斯牛顿法
•     5.2.2 拟牛顿法
•     5.2.3 有限内存拟牛顿法
•   5.3 三维大地电磁测深法各向同性反演
•     5.3.1 目标函数及其梯度计算
•     5.3.2 水平地形模型
•     5.3.3 带地形模型
•     5.3.4 实际地质背景模型
•     5.3.5 实测数据反演
•   5.4 本章小结
• 第6章 三维大地电磁测深法各向异性反演研究
•   6.1 目标函数及其梯度计算
•   6.2 简单各向异性异常体的三维反演测试
•   6.3 多异常体各向异性的三维反演测深
•   6.4 带地形的各向异性体的三维反演测试
•   6.5 Coprod2 实测数据各向异性反演的应用研究
•     6.5.1 研究现状及其各向异性指示
•     6.5.2 极性图分析
•     6.5.3 三维各向异性反演
•   6.6 USArry实测数据三维各向异性反演的应用研究
•     6.6.1 各向异性分析
•     6.6.2 三维各向异性反演
•   6.7 本章小结
• 第7章 结论与展望
•   7.1 结论
•   7.2 展望
• 参考文献
• 作者简介及在学期间所取得的科研成果
• 致谢

文章来源

类型: 博士论文

作者: 曹晓月

导师: 殷长春

关键词: 地球物理电磁测深法,大地电磁测深法,三维正反演,非结构有限元法,面向目标自适应,有限内存拟牛顿法,各向异性介质

来源: 吉林大学

年度: 2019

分类: 基础科学,工程科技Ⅰ辑

专业: 地质学,地球物理学,矿业工程

单位: 吉林大学

分类号: P631.325

总页数: 181

文件大小: 17418K

下载量: 355