基于矢量有限元法的三维张量CSAMT正演模拟.PDF

2015 年03 月 APPLIED GEOPHYSICS （应用地球物理） 第12 卷，第1 期 基于矢量有限元法的三维张量CSAMT 正演模拟//Three-dimensional tensor controlled-source electromagnetic modeling based on the vector finite-element method, APPLIED GEOPHYSICS, 2015, 12(1), P. 35-46. DOI: 10.1007/s11770-014-00476-2 基于矢量有限元法的三维张量CSAMT 正演模拟 1 ◆1 2 3 4 胡英才 ，李桐林 ，范翠松 ，王大勇 ，李建平 （1. 吉林大学地球探测科学与技术学院，长春 130026 ；2. 中国地质调查局天津 地调中心，天津 300170 ；3. 中国地质科学院地球物理地球化学勘查研究所，廊 坊 065000 ；4 山东科技大学地质科学与工程学院，青岛 266590 ） 摘要：标量CSAMT 只适合一维及测量方向与构造方向垂直的二维情况，对于复 杂的三维地电结构，CSAMT 需采用张量测量。本文试图采用矢量有限元法实现三 维张量 CSAMT 的正演模拟。为了验证算法的正确性，本文在层状介质中计算了 三维 CSAMT 远区的电场，磁场及阻抗张量，并且与层状介质中的理论解进行了 比较，接着还模拟了均匀半空间中含有三维异常体的模型，并且分析了四个阻抗 张量、视电阻率及阻抗相位的响应特征。得出如下结论：采用矢量有限元法来模 拟三维张量CSAMT，其电磁场及阻抗张量的实虚部计算精度都比较高，并且该方 法本身满足电场法向不连续，不用进行散度校正。 关键词：可控源音频大地电磁法，正演，阻抗张量，矢量有限元法 引言 1997a, 1997b; Wang et al., 2014)。 在三维正演模拟中，相对于有限差分和积分 方程而言，有限元法因模拟精度高且可以不规则 可控源音频大地电磁法（简称CSAMT ）由 网格剖分而备受国内外众多学者的青睐。有限元 于采用人工源发射，具有抗干扰能力强和信噪比 法分为节点有限元和矢量有限元法等。 高等众多优点，在深部找矿，地下水勘查等领域 在节点有限元法的模拟中，Zhang (2009)等 得到了广泛的应用(Wu et al., 1996; Hu et al., 基于电场矢量波动方程实现了三维节点有限元 2012) 。以往在野外的数据采集中，CSAMT 大多 法的可控源电磁法模拟，其中利用伪 函数代替 采用标量测量，在资料处理及解释中，多数也只 源避免了源的奇异性，另外，为避免电场在界面 进行二维甚至是一维反演解释，随着地质勘探深 上的法向分量连续，通过增加一个罚项来做散度 度的加大，要求越来越高，对于三维的复杂地电 较正，减少了伪解的出现。Xu (2010)等采用电磁 结构或者测量方向