大地电磁层状模型数值模实验实验报告.doc

实验报告 课程名称： 地电学 课题名称： 大地电磁层状模型数值模拟实验 专 业： 地球物理学 姓 名： xx 班 级： 06xxxx 完成日期： 2016 年 11月26日 目录 实验名称……………………………………………………3 实验目的……………………………………………………3 实验要求……………………………………………………3 实验原理……………………………………………………3 实验题目……………………………………………………4 实验步骤……………………………………………………4 实验整体流程图……………………………………………8 程序及运行结果 一、实验名称 大地电磁层状模型数值模拟实验 二、 实验目的 学习使用Matlab编程，并设计大地电磁层状模型一层，二层，三层正演程序 在设计正演程序的基础上实现编程模拟 MATLAB软件基本操作和演示 . 三、 实验要求 利用MT一维测深法及其相关公式，计算地面上的pc视电阻率和ph相位，绘制视电阻率正演曲线和相位曲线并分析。 利用Matlab软件作为来实现该实验。 四、 实验原理 （一）、正演的概念： 正演是反演的前提。在实际地球物理勘探中，一些模型的参数是不容易确定的，如埋藏在地下的地质体模型的岩性、厚度、产状等参数，我们把这些描述未知模型的参数的集合定义为“模型空间”。为了获得这些模型参数，可以利用那些可以直接观测的量来推测，而这些能够直接观测的量的集合则被称作“数据空间”。如果把模型空间中的一个点定义为m，把数据空间中的一个点定义为d，按照物理定律，可以把两者的关系写成 式中，G为模型空间到数据空间的一个映射。我们把给定模型m求解数据d的过程称为正演问题。 （二）、MT一维正演模型简介 大地电磁法作为一种电磁类勘探方法，它的模型参数为一组能够表征地球物理勘探目标体的电性参数，即目标体电阻率和相应层的层厚度。所谓一维模型，即介质在三维空间中沿两个方向上模型参数是不变的，只在另一个方向上特征属性会变化。在此一维模型即指水平层状一维介质，即介质只在沿垂直于地面上的方向上电性（电阻率）变化，在另外两个方向上保持不变的典型特征，所以就构成一组电阻率不同的电性层，抽象出来即是一组由电阻率及对应的层厚度构成的电性层数。 根据正演问题的概念，构成正演的元素有3个，即模型、测量数据和模型到数据的映射。对模型来说比较简单，即为水平层状一维介质模型。我们知道大地电磁法属于一种天然的交变电磁场的地球物理勘探方法，所以它的测量数据一般为大地电磁场的电场和磁场分量。而将以上两者联系起来的关系—映射则是二者之间的物理规律，由于大电磁场场源的性质，可将大电磁场看作是垂直入射的平面波，通过地下介质传播到地面上。在这个过程中，大地电磁场遵循电磁场的普遍规律，即Maxwell方程组。在大地电磁法中，我们利用在地面上的视电阻率和相位进行后续的解释工作，所以正演的数据空间需转化为视电阻率和相位。综上所述，MT一维正演即求解水平层状一维介质对垂直入射平面波在地面上的视电阻率和相位响应。 五、 实验题目 1、利用MT一维测深法及其相关公式，计算地面上的pc视电阻率和ph相位。 2、绘制视电阻率正演曲线。 3、绘制相位曲线。 六、 实验步骤 大地电磁法一维正演具有以下的基本推导思路：从大地电磁场满足的基本方程—麦克斯韦方程组出发，结合大地电磁场的特点，推导出单一方向的波动方程；然后，结合水平层状介质的边界条件，推导出能够表示地面波阻抗的递推式；最后根据视电阻率和相位的定义式，得出水平层状介质的大地电磁场响应函数（视电阻率和相位）。 我们知道麦克斯韦方程组有4个基本方程构成，另外还有3个本构关系将4个基本方程联系起来，其具体的形式如式： 其中E和H为电场强度和磁场强度，j为电流密度，D为电位移矢量，B为磁感应强度，σ、μ、ε分别为电导率、磁导率和介电常数。 由于大地电磁法应用的频率都很低，一般f10hz，这时在导电介质的位移电流?D/?t与传导电流j相比可以忽略不计所以麦克斯韦方程组可以简化为以下形式 考虑在谐变场的情况下，对上式前两式两边取旋度，并根据矢量分析公式可得出波动方程的形式 由于是一维层状介质，所以在笛卡尔坐标系下电磁场在水平方向上是不变的，故只需研究沿Z轴向下方向上的电磁场分量。由波动方程上式知： 其中Km为第m层的复波数， 求解得： 所以，波阻抗Z可求得为： 其中Z0m为第m层的特征阻抗： 我们知道同一层内部积分常数Am和Bm是相同的，因