25维海洋电磁正演模拟的基本原理.docxVIP

2.5维海洋电磁正演模拟的基本原理所谓的2.5维问题指的是二维模型三维源问题。【赵光茂.带地形的复电阻率2_5维电磁场正反演研究，博士论文，2009:第二章】二维模型的电性结构如图1所示，将原点取在地面，方向与走向平行，且该方向上的电导率()、介电常数()和磁导率()不变，在平面内这些参数可以变化。电偶极源如图中箭头线所示。异常体 图1 二维介质模型1 2.5 维电磁场耦合微分方程推导由电磁场理论，麦克斯韦方程在时谐场的复矢量形式（假设电磁场是以变化的谐变场）可知：将关系式：，带入并整理得：式中，是电场矢量；为磁场矢量；为电流密度；为外加的电源矢量；为角频率；为磁导率；为电导率；为介电常数；虚数单位。按构造方向（y）按以下形式作傅立叶变换（为波数）所以:由由、式得：由、式得：其中：,,,。将、、带入、可得关于走向的、的对偶微分方程： 可得：如果从上式解出和，则利用辅助场、、的表达式就可计算出其它电场和磁场分量。于是 2.5 维电磁场的计算问题归结为求 域电场和磁场分量定解问题的数值解。场源等效偶极子的电磁场在源点附近一般是奇异的，也就是说源点附近的电磁场变化是非常迅速的。即使有精细的网格剖分、采用高阶插值函数，也不能很好的模拟场源附近的变化情况。通常的解决办法是通过对背景场和二次场给予分别求解。但是这种办法也有不足，它对于地表起伏或复杂地质体存在时不能找到合适的背景构造。一般可采用一个近似函数来等效电偶极子场源的作用，将源分布于一定范围内，避免在源点处产生奇异性，从而不用区分背景场和异常场，直接求解总场值。这个近似函数称为伪函数。于是一个电偶极源的三维函数表示为。式中、分别表示源点的坐标，参数控制源分布的宽度和幅值。网格剖分用有限元法求解微分方程的定解问题首先要对研究区域进行剖分，节点与单元遍及整个研究区域，将连续函数的求解化为节点或面上的离散函数值来进行求解。因此，网格的大小、网格的类型和网格和疏密，对函数值的计算会产生重要的影响。采用任意四边形单元剖分，通过调整水平分布各节点的位置来实现起伏地表的模拟，断面剖分网格如图所示。空气空气海水海水海水围岩层围岩层二维地电断面剖分网格示意图但是它在公共边界上不满足位移连续性条件，而且对形状各不相同的单元进行计算，计算机自动实现起来也较为麻烦。为解决这一问题，必须采用等参数坐标变换，对任意形状的四边形网格进行规则化处理。为了更接近实际情况，网格剖分中地形变化或者曲线边界采用合适形状的多边形单元进行剖分。这种剖分形式可以模拟复杂的地形和边界，调节网格的疏密结构，但对于任意四边形的单元直接进行数值计算，工作量过大，且不利于计算机自动化实现。为了将局部（自然）坐标中几何形状规则的单元转换成总体（笛卡儿）坐标中几何形状扭曲的单元，以满足对一般形状求解域进行离散化的需要，必须建立一个坐标变换。利用等参单元来处理，在规则单元域内进行计算，不管各个积分形式的矩阵表示的被积函数如何复杂，仍然可以方便地采用标准化的数值积分方法计算。对于本文中剖分网格中的任意四边形单元进行等参变换，本文将直角坐标系内不规则的网格单元转换成自然坐标体系下的正方形网格，便于计算，同时又容易实现起伏地形的模拟。下图就是等参变换的图示。图中右边为剖分形成的不规则四边形单元，左边为单元等参变换后在新坐标系中的映象。为了方便新坐标系内的数值计算，通常将正方形单元设置为关于ξ，η坐标轴的坐标原点中心对称，并取正方形四边的ξ和η值等于±1。坐标下不规则四边形单元的八个节点分别与ξ，η坐标下正方形单元的八个节点对应。确定对应规则后，选择适当的的插值函数，将求解域中的积分表达式完成离散化。单元等参变化图(a)的局部自然坐标系中的正方形单元称为基本单元。图 (b)中的在整体坐标系中的任意四边形单元称为实际单元。两单元内的参变量函数采用相同数目的节点参数和相同的形函数进行变换。可以在局部自然坐标系中的规则单元上进行单元分析，然后再映射到实际直角坐标系中的单元上进行计算，等参单元同时具有计算精度高和适用性好的特点。系数矩阵4.1加权余量法(地球物理中的有限单元法-徐世浙 p300,伽辽金法比泛函法简单，且应用更广，因为有时有些方程找不到合适的泛函)加权余值法要求在整个求解域中加权的残数(余值)和为零，假定方程为：则余值为，乘以加权系数，则加权余值为，从微分方程出发，得到线性方程组。从有限元基本单元出发，获得相应的形函数。将形函数作为加权系数代入上式子中(加权系数为时称为伽辽金法)，将各单元的近似解和形函数代入上式，可获得个方程的线性方程组。[陈小斌，大地电磁正反演新算法研究及资料处理与解释的可视化集成系统开发,博士论文，2003：p41]4.2方程离散通过式可得到余量：由伽辽金法，可得目标函数：式中是第个单元中第个节点的型函数。通过 G