二维电导率各向异性介质电磁感应效应的有限元算法案例.docxVIP

二维电导率各向异性介质电磁感应效应的有限元算法 李予国 摘要 二维电性各向异性大地电磁场通过有限元方法计算获得。模型由包含各向异性块的层状介质构造组成。每个块或者层可能是由3×3的电导率张量限定。正演问题可以归结为两个耦合的关于垂直-平行于场分量的和。它们由有限元法数值地解出。线性有限元的解系是用预条件共轭梯度法求取的。之后，地表垂直-平行的场分量和，利用样条插值通过对和进行数值微分得到。 二维有限元算法通过对比二维有限差分的解证明是正确的。三个类型的模型的大地电磁响应用来证明各向异性的影响：水平，垂直，倾斜各向异性。第四个模型用来模拟在剪切和俯冲带各向异性的影响。这些模型的响应模拟了电阻率曲线在长周期和具有明显的主对角线元素的张量主抗存在室的分离，正如以前观测到的。 1、简介 近些年来，越来越多的关注投入到了对电性各向异性的电磁感应方面上的研究，尤其是为了试图完全理解对更长周期的大地电磁的观测数据。加拿大地区的大地电磁测量显示了电性的各向异性在下地壳和上地幔（Kellett等，1992；Mareschal等1995）。大的来自围绕德国深钻点（KTB）各向异性大地电磁曲线解释为上地壳高度的电各向异性（Eisel 和 Haak 1999）。Rasmussen (1988)使用一个深地壳层内的各向异性模型来解释瑞典南部的大地电磁横断面数据。层状构造的各向异性效应首先被O’Brien 和Morrison (1967)研究。Reddy 和 Rankin (1975)首先开始研究二维各向异性模型，他们只考虑水平各向异性的工作。最近，Osella 和Martinelli (1993)计算了一些圆滑不规则边界有一定的主轴旋转角的大地电磁响应。Schmucker (1994)提出一个计算不均匀薄的覆盖体在层状半空间的电磁感应效应，其中可以含有一个或者两个各向异性的电导层。 用有限元方法，Pek 和 Verner (1997) 还有 Weidelt (1996) 模拟了广义二维、三维各向异性构造，分别具有任意的主轴旋转角。 本文中，二维模型被重新研究，但是使用有限元（FE）的方法。首先，我们详细描述有限元程序的数值近似。之后，我们展示有限元程序模拟不同的简单测试模型的响应。我们的结果和Pek Verner (1997)的有限差分程序进行了对比。最后，我们计算三种类型的各向异性介质的大地电磁响应：水平，垂直，以及倾斜各向异性。我们最后以模拟一个地质剪切带或者板块俯冲带的地质背景下的模型结束。本文这样的构造是为了展示一个有限元策略去实现二维各向异性介质的电磁感应。尽管有限差分对于这类问题是有效的，有限元解仍然是需要的，因为有限元和有限差分都有自己的特殊的优势还可以彼此互相核对。另外，有限元方法可以模拟非矩形线性构造的真实地球构造。 2、边界值问题 考虑到图1中的二维模型。阴影异常区嵌入到一个简单的层状构造中，这个构造包括n层而且底部的层延伸到无穷远。为了简单，异常区域没有显示内边界作为一个定向独立的电导率，它至少和周边的第j层的电导率张量不同。然后我们的程序允许更一般的模型，异常区域可以细分成不同的不均匀块而且可以接触不同的层。 异常区域延伸到y正向到无穷远，例如模型可以出现到一个新的正常构造。在这种状态下算法排除源的感应。模型在走向方向x方向是不变的。感应电磁场在x方向也是不变的，尽管如此场的矢量也有三个分量对于极化的磁源平行或者垂直走向，由于各向异性。因此，通常区分TE和TM模式的各向同行构造就变得不合理了。由此就有结果，当主磁场矢量垂直轴向，电场将有主场方向的分量，导致电荷产生于界面（除了倾斜各向异性的例子）。 假设一个时谐变量，在准静态近似的电磁场的情况下 (1) 其中是自由空间的磁导率，且 为电性电导率张量。它是对称的，当旋转到主轴方向（x,y,z）上去，可以表示为 图1.本文中考虑的二维各向异性模型 在特别的二维例子中，方程（1）消退为常数电导率的均匀块 很明显如果垂直-平行分量和得到了，剩下的分量和，和可由和的空间导数得到。方程（2）-（7）可以得到两个关于分量和的二阶偏微分方程： 从方程（8）和（9）很明显看出各向异性和独立的垂直-平行分量和通过一阶偏微分方程耦合在一起。因此，在这种情况下不存在独立的TE、TM模式对于广义各向异性情况。因此，这些方程必须同时利用和求解。 之后叙述了不同的空间形式的各向异性的模拟结果。如图1一个均匀的异常区域具有各向同性的