大地电磁法及其探测中的应用-北京大学学术报告.ppt

某地区实测的MT视电阻率和相位曲线 与地下资源 （石油、媒、 金属矿、地热等） 密切相关 地壳深部结构密切相关 深层原因 与地下资源 （石油、媒、 金属矿、地热等） 密切相关 与地震现象 密切相关 岩石的电导率及其与地震和地下资源的关系 影响岩石导电性的因素 岩石结构 物质组成 含水量 温 度 压 强 无物理理－化学变化 物理－化学变化 孔隙度越大，连通性越好，导电性越好 含水量越大，导电性越好 温度越高，导电性越好 流体静压力越高，导电性越差 矿物成分不同，导电性不同 热脱水反应 高温熔融 使岩石的导电性迅速增加 部分熔融变可极大地增加岩石的导电性 地震 地震引起的视电阻率幅度的变化 五大连池火山区电性结构 的三维成象图 长白山火山区NE测线的 二维反演结果 长白山火山区电性结构 的三维成象图 青藏高原东北缘的大地电磁探测 吉林地热田探测的二维反演结果 Lg(Res/Ωm) 存在的一些问题及初步的研究 高速度高精度的三维正反演问题 畸变问题和静位移问题 二维反演中极化模式选取的问题 地形影响问题 反演结果的有效深度问题 反演结果的非唯一性问题 可视化集成系统的开发 应用问题：怎么与大陆动力学直接结合起来？ 野外实际地电结构一般是三维的！ 实测数据二维反演结果 二维反演中数据极化模式的选择 (蔡军涛，陈小斌，赵国泽，未发表） 异常体沿走向X的延伸是变化的，分别取4、6、8、10、20、40km 二维反演结果的比较 2.地形的影响 （未发表） 纯地形模型 反演网格对地形的拟合 TE和TM的互换 反演结果的有效深度 ——复杂模型趋肤深度的计算 （未发表） 反演结果的有效深度：剖面长度的影响 （杨静、陈小斌、赵国泽，中国地球物理年刊，2009） 反演结果的非唯一性分析 可视化集成系统的开发 可视化集成系统的开发 汶川地震初步的动力学模型 前山断裂 中央断裂 后山断裂 岷江断裂 龙日坝断裂 正演指的是什么？ 正演指的是对于一个给定的模型，在一定激发源的作用下，根据一定的物理原理求其响应的过程； 在地球物理中，正演一般对应于求解一数理模型的边值或初值问题； 正演指的是什么？ 正演指的是对于一个给定的模型，在一定激发源的作用下，根据一定的物理原理求其响应的过程； 在地球物理中，正演一般对应于求解一数理模型的边值或初值问题； 正演指的是什么？ 正演指的是对于一个给定的模型，在一定激发源的作用下，根据一定的物理原理求其响应的过程； 在地球物理中，正演一般对应于求解一数理模型的边值或初值问题； 在大地电磁测深中，正演问题是，假设大地的电性结构已知，当均匀平面电磁波垂直入射到大地表面时，求在地表产生的电磁响应。 二维模型：场可解耦为两组模式 二维正演：边值问题 TM模式： TE模式： 二维正演方法 1）积分方程法 2）有限差分法 3）有限元法 有限元直接迭代算法（陈小斌，1999，2000） 二维正演计算的简单实例 TE模式 TM模式 视电阻率分布 二维正演的简单算例 TE模式 TM模式 阻抗相位分布 复杂模型的正演结果 三维正演：边值问题 在三维情况下，电磁场不能解耦成两组独立的场，这时必须直接求解矢量波方程。 以上方程包含了一个隐含条件：求解域的电阻率是分块均匀的。 三维正演方法 1）积分方程法 2）交错网格有限差分法 3）矢量有限元法（棱边元） 如何理解反演？ 反演是指根据实测的数据来反推产生这些数据的系统内在信息的一种数学物理过程。 反演的两个基本条件：实测的数据和一个先验模型系统。 通常的最小二乘多项式拟合就可以看成是一个反演过程。参与拟合的数据就是反演中实测的数据，“多项式”这种函数形式就是“先验模型系统”。 如何理解反演？ 对于大地电磁测深而言，“实测的数据”就是在地表实测的视电阻率、相位等数据；“先验模型系统”是对地球电导率模型的假设（一维、二维还是三维？），以及在此假设基础上的正演实现过程。更明确的说，这里的“先验模型系统”就是指的是“一维正演”过程、“二维正演”过程或“三维正演”过程。 对于大地电磁测深而言，所谓待反演的“系统内在信息”指的就是电导率结构。 大地电磁测深反演就是根据地表实测的视电阻率、相位等数据来求取大地深部电导率结构的过程，该电导率结构的正演响应能极好地拟合视电阻率、相位等实测数据。 ？ ？ 正演 反演 二维正演计算的简单实例 TE模式 TM模式 视电阻率分布 二维正演的简单算例 TE模式 TM模式 阻抗相位分布 复杂模型的正演结果 手工量板法 反演问题和反演方法的分类 反演问题主要分两类：线性问题和非线性问题。大地电磁测深反演属于非线性反演问题。 反演方法也有线性反演和非线性反演之分。 线性反演方法是针对线性反演问题发展起来的，但也被广泛应用于解决非线性问题，这时称为非线