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

二维电导率各向异性介质电磁感应效应的有限元算法 李予国 摘要 二维电性各向异性大地电磁场通过有限元方法计算获得。模型由包含各向异性块的层状介质构造组成。每个块或者层可能是由3×3的电导率张量限定。正演问题可以归结为两个耦合的关于垂直-平行于场分量的和。它们由有限元法数值地解出。线性有限元的解系是用预条件共轭梯度法求取的。之后，地表垂直-平行的场分量和，利用样条插值通过对和进行数值微分得到。 二维有限元算法通过对比二维有限差分的解证明是正确的。三个类型的模型的大地电磁响应用来证明各向异性的影响：水平，垂直，倾斜各向异性。第四个模型用来模拟在剪切和俯冲带各向异性的影响。这些模型的响应模拟了电阻率曲线在长周期和具有明显的主对角线元素的张量主抗存在室的分离，正如以前观测到的。 1、简介 近些年来，越来越多的关注投入??了对电性各向异性的电磁感应方面上的研究，尤其是为了试图完全理解对更长周期的大地电磁的观测数据。加拿大地区的大地电磁测量显示了电性的各向异性在下地壳和上地幔（Kellett等，1992；Mareschal等1995）。大的来自围绕德国深钻点（KTB）各向异性大地电磁曲线解释为上地壳高度的电各向异性（Eisel 和 Haak 1999）。Rasmussen (1988)使用一个深地壳层内的各向异性模型来解释瑞典南部的大地电磁横断面数据。层状构造的各向