被动源电磁测深三维有限元数值模拟.pdfVIP

第37卷第6期 成都理工大学学报(自然科学版) V01．37No．6 2010年12月 JOURNALOFCHENGDUUNIVERSITYOFTECHN()LOGY(Science＆．Technology [文章编号]1671—9727(2010)06—0654—06 被动源电磁测深三维有限元数值模拟 张继锋1 汤井田2 王 烨3 肖 晓2 (1．长安大学地质工程与测绘学院，西安710054；2．中南大学信息物理工程学院，长沙410083； 700) 3．河南省洛宁县人民政府。洛阳471 [摘要]从基于电场的双旋度方程出发，针对被动源电磁测深法中的三维正演问题进行了有限 元数值模拟。首先推导了被动源电磁法中电场所满足的矢量波动方程，然后从变分原理出发， 得到其交分形武，通过有限元分析，最后形成大型复系数线性方程组。采用LDLl‘预处理共轭 梯度法，降低了矩阵的条件数，提高了解方程组的速度。建立了一个典型的3层地电模型，得 到了地表电磁场分布及其主测线视电阻率和阻抗相位，视电阻率平均相对误差为2．06％，阻 抗相位平均相对误差为0．734％，满足实际要求，同时也验证了程序的正确性。 [关键词]电磁测深；有限元，数值模拟I大地电磁 [分类号]P631．325 [文献标识码]A 被动源电磁测深主要指依靠天然场源的各种 有模拟复杂地电模型的能力，且通用性强，因而在 频段进行大地电磁测深，由于不需要人工场源。野 地球物理电磁法正反演中具有更为广泛的应用前 外数据采集方便，探测深度大，因而得到了广泛的 景。Bader【l纠(2001)采用库伦规范下的矢量位和 应用。目前，大地电磁测深的数据处理和正反演 标量位公式，结合散度条件，采用四面体网格剖分 算法以一维、二维算法为主，发表的论文也数以百 和局部加密技术，对井中可控源电磁法响应进行 Mit- 计[1~4]；而三维电磁法正反演由于数值计算过于了三维数值模拟，取得了良好的效果。Yuji 复杂并受计算机内存、计算时间和计算量的限 制【51，实际应用比较少，还处于探索阶段。 起来，采用T—Q公式，对2个三维COMMEMI模 近年来，随着计算机技术的快速发展，许多学 型进行了正演计算和误差分析，验证了程序的正 者对三维电磁法数值模拟进行了大量的研究，目 确性。MyungJinNam【14J(2007)提出基于六面 et 体的矢量有限元方法模拟三维大地电磁响应，并 前以积分方程法(Wannamaker，1991L61；Xiong and a1．，1997[7])，有限差分法(Mackie，Smith采用等参单元方法对三维地形进行了模拟，为进 and Madden，1994[8]；Sasaki，1999[9]．Newman 一步实现带地形的三维电磁模拟奠定了基础。在 et J(2000)采用矢量有限元方法对水 Alumbaugh，2002【10])和有限单元法(Zunoubi国内，闫述等【15 eta1．，2001 nZ．。1999c11]；BadeaL12J；Mitsuhata平电偶源下的电磁响应进行了三维正演，给出了 and Uchida，2004[¨])为主流。以变分原理为基相应的结果。黄临平和戴世坤【53(2003)推导了地 础的有