大地电磁测深一维正演-地电学实验报告讲诉

1、实验报告课程名称： 地电学课题名称： 大地电磁层状模型数值模拟实验 专业：地球物理学姓 名： xx班级：06xxxx完成日期： 2016 年 11 月 26 日I / 14地电学实验报告 Copyright060141-Weway目录一、实验名称 3二、实验目的 3三、实验要求 3四、实验原理 3五、实验题目 4六、实验步骤 4七、实验整体流程图 8八、程序及运行结果 9九、实验结果分析及体会 142/ 14地电学实验报告 Copyright060141-Weway、实验名称大地电磁层状模型数值模拟实验、 实验目的1）学习使用 Matlab编程，并设计 大地电磁层状模型一层， 二层，三层正演程

2、序2）在设计正演程序的基础上实现编程模拟3）MATLAB 软件基本操作和演示三、实验要求（1） 利用 MT 一维测深法及其相关公式，计算地面上的 pc 视电阻率和 ph 相位，绘 制视电阻率正演曲线和相位曲线并分析。（ 2） 利用 Matlab 软件作为来实现该实验。四、实验原理一）、正演的概念：正演是反演的前提。在实际地球物理勘探中，一些模型的参数是不容易确定 的，如埋藏在地下的地质体模型的岩性、厚度、产状等参数，我们把这些描述未 知模型的参数的集合定义为“模型空间” 。为了获得这些模型参数，可以利用那 些可以直接观测的量来推测， 而这些能够直接观测的量的集合则被称作 “数据空 间”。如果把

3、模型空间中的一个点定义为 m，把数据空间中的一个点定义为 d， 按照物理定律，可以把两者的关系写成式中，G 为模型空间到数据空间的一个映射。我们把给定模型m 求解数据 d的过程称为正演问题。二）、 MT 一维正演模型简介大地电磁法作为一种电磁类勘探方法， 它的模型参数为一组能够表征地球物 理勘探目标体的电性参数，即目标体电阻率和相应层的层厚度。所谓一维模型， 即介质在三维空间中沿两个方向上模型参数是不变的， 只在另一个方向上特征属3 / 14地电学实验报告 Copyright060141-Weway性会变化。 在此一维模型即指水平层状一维介质， 即介质只在沿垂直于地面上的 方向上电性（电阻率）

4、变化，在另外两个方向上保持不变的典型特征，所以就构 成一组电阻率不同的电性层， 抽象出来即是一组由电阻率及对应的层厚度构成的 电性层数。根据正演问题的概念，构成正演的元素有 3 个，即模型、测量数据和模型 到数据的映射。 对模型来说比较简单， 即为水平层状一维介质模型。 我们知道大 地电磁法属于一种天然的交变电磁场的地球物理勘探方法， 所以它的测量数据一 般为大地电磁场的电场和磁场分量。 而将以上两者联系起来的关系 映射则是二 者之间的物理规律， 由于大电磁场场源的性质， 可将大电磁场看作是垂直入射的 平面波， 通过地下介质传播到地面上。 在这个过程中， 大地电磁场遵循电磁场的 普遍规律，即

5、Maxwell 方程组。在大地电磁法中， 我们利用在地面上的视电阻率 和相位进行后续的解释工作， 所以正演的数据空间需转化为视电阻率和相位。 综 上所述， MT 一维正演即求解水平层状一维介质对垂直入射平面波在地面上的视 电阻率和相位响应。五、实验题目1、利用 MT 一维测深法及其相关公式，计算地面上的 pc 视电阻率和 ph 相位。2、绘制视电阻率正演曲线。3、绘制相位曲线。六、实验步骤大地电磁法一维正演具有以下的基本推导思路： 从大地电磁场满足的基本方程 麦克斯韦方程组出发， 结合大地电磁场的特点， 推导出单一方向的波动方程； 然后， 结合水平层状介质的边界条件，推导出能够表示地面波阻抗的

6、递推式；最后根据视 电阻率和相位的定义式，得出水平层状介质的大地电磁场响应函数（视电阻率和相 位）。我们知道麦克斯韦方程组有 4 个基本方程构成， 另外还有 3 个本构关系将 4个基 本方程联系起来，其具体的形式如式：4/ 14地电学实验报告 Copyright060141-Weway其中 E和 H 为电场强度和磁场强度，j 为电流密度， D 为电位移矢量， B 为磁感应强度， 、 、分别为电导率、磁导率和介电常数。由于大地电磁法应用的频率都很低，一般f10hz ，这时在导电介质的位移电流?D/?t与传导电流 j 相比可以忽略不计。所以麦克斯韦方程组可以简化为以下形式：并根据矢量分析公式可得出

7、考虑在谐变场的情况下，对上式前两式两边取旋度， 波动方程的形式由于是一维层状介质，所以在笛卡尔坐标系下电磁场在水平方向上是不变的，故只需研究沿 Z 轴向下方向上的电磁场分量。由波动方程上式知：其中 Km 为第 求解得：m层的复波数，所以，波阻抗Z 可求得为：其中 Z0m 为第 我们知道同一层内部积分常数m 层的特征阻抗：Am和 Bm 是相同的，因此层内不同深度处的波阻抗5 / 14地电学实验报告 Copyright060141-Weway可以通过积分常数联系起来。为此，将上作如下变换：? 则有：若取底面处波阻抗代入上式中求出 Bm/Am ，然后代入上上式求取顶面的波阻抗， 则可把同一层顶面和底

8、面的波阻抗联系起来，并消去积分常数Am 和 Bm 。记 Zm为第 m 层的顶面阻抗，底面的波阻抗等于第 M+1 层顶面的波阻抗，则结果如下：将上式代入阻抗的定义式：其中同样将上式写成如下形式：其中 Rm为第 m 层的反射系数。所以就得到了顶面波阻抗的递推公式： 最底层为 n 层。6 / 14地电学实验报告 Copyright060141-Weway而正演则是要求出在地面上的视电阻率和相位响应，对相位来说即是波阻抗相位，也就是波阻抗所对应的复数的幅角。对视电阻率来说，根据视电阻率的定义有：所以就有 n 层层状介质的视电阻率响应为：由特征阻抗公式及变换式可得：从以上的递推过程可以看出，根据反射系数

9、（波阻抗）的递推公式可以计算出地面上的视电阻率表达式以及阻抗相位的表达式，可用于进一步的程序实现。7/ 14地电学实验报告 Copyright060141-Weway七、实验整体流程图或算法为了测试该 MT 一维正演程序的应用效果，考虑选取几种典型的地电断面作为正演 程序的输入模型，即二层模型、三层模型。根据电性层各层电阻率的相互关系，二 层模型可以分为 G 型和 D 型，而三层模型则分为 A 型、H 型、K 型和 Q 型这四种 类型，至于多层层状（大于三层情况下）介质则可以分解为上述的几种简单类型。 我们已经知道，以上所述几种典型模型的视电阻率响应函数特征，如果将以上模型 输入到本次所写的程

10、序中，则可以作为测试本程序是否可行的依据。大地电磁一维测深模拟输入输出模块8 / 14地电学实验报告 Copyright060141-Weway大地电磁一维测深运算模块八、 程序及其运行结果MT 一维测深运算程序代码：9 / 14地电学实验报告 Copyright060141-WewayG型曲线为两层模型曲线，其各层电阻率的关系为 12，程序 正演时取 1=1000 m， 2=100 m， h1=1000m。正演理论结果如下 图 2-4 。从图上可以看出 D型曲线在短周期视电阻率较大，随着周期 T 变 长，视电阻率也相应的减小， 但 D 型曲线仍存在尾支渐进线， 渐近线与第 二层的真电阻率相等

11、。A型： A型曲线为三层模型曲线，其各层电阻率的关系为 123， 程序正演时取 1=10 m， 2=100 m， 3=1000 m， h1=h2=1000m。正演理论结果如下图 2-5 。从图上可以看出 A型曲线在短 周期视电阻率较小，随着周期 T变长，视电阻率也相应的增大，但 A 型 曲线仍存在尾支渐进线， 渐近线与第三层的真电阻率相等。 从以上特征来 看， A型与 G型曲线有相似的特点。11 / 14地电学实验报告 Copyright060141-WewayK型： K型曲线为三层模型曲线，其各层电阻率的关系为 13， 程序正演时取 1=10 m， 2=100 m，1=10 m ，h1=h2

12、=1000m。 正演理论结果如下图 2-6 。从图上可以看出 K 型曲线在短周期视电阻率较 小，随着周期 T 变长，视电阻率也相应的增大，达到一峰值后逐渐减小， 这一峰值与第二层的电阻率有关， 但并不等于该电阻率， 说明该峰值还受 相邻层的影响。 但 K型曲线仍存在尾支渐进线， 渐近线与第三层的真电阻 率相等。12 / 14地电学实验报告 Copyright060141-WewayH型： H型曲线为三层模型曲线，其各层电阻率的关系为 1223， 程序正演时取 1=1000 m， 2=100 m， 3=10 m ， h1=h2=1000m。正演理论结果如下图。从图上可以看出 Q型曲线在短周期 视电阻率较大， 随着周期 T变长，视电阻率也相应的减小， 但 Q型曲线仍 存在尾支渐进线，渐近线与第三层的真电阻率相等。13 / 14地电学实验报告 Copyright060141-Weway九、 实验体会实验中设计了几个模型参数， 经过正演计算以及 matla