三维大地电磁地形影响的有限元数值模拟.pdfVIP

三维大地电磁地形影响的有限元数值模拟 师学明+㈣歌田久司@王家映①吴文鹂②王晓明①张云舻 ①中国地质大学地球物理与空问信息学院武汉430074． (孚日本东京大学地震研究所海半球观测研究中心 东京113—0032 ③中国地质科学院地球物理与地球化学勘查研究所廊坊065000) 摘要 我国南方地区，矿体丰富．但在这些地区开展电磁方法，必然要受到由于地形起伏的影 响，和地质体本身的三维复杂构造，严重困扰了大地电磁数据的资料解释．二维地形影响与 张翔等，1999)，但是三维地形MT影响的数值模拟，难度极大，目前只有4篇文献见报道．徐 等(2003)等用有限差分和矩形剖分，在微机上实现了MT三维地形影响的数值模拟．用有 限元进行三维MT地形影响的数值模拟的文献未见报道．因此，开展三维地形影响有限元模 拟理论方面的研究，不仅具有重要的理论意义，而且对带地形的三维MT反演具有十分重要 的实用价值． 本文从理论上综合矢量有限元方法和散度校正技术，基于电场提出了带电流散度校正的 finiteelementmethodcombinedwith correctionbasedonthe 矢量有限元算法(Vector divergence electric 际标准三维地电模型(COMMEMI3D．1和3D一2)模型进行了试算．对COMMEMI3D．1和3D．2 模型的数值模拟结果与国际同行已发表的结果基本一致，说明了该方法计算精度高，计算结 果准确可靠．vFEE++方法采用了矢量有限元消除了标量有限元方法中的“矢量寄生虫”问题， 同时也采用了散度校正技术来消除电流散度较大的伪解，从而保证电场解满足物理条件．在 带散度校正的矢量有限元方法过程中，矢量有限元方程、磁场散度方程和散度校正技术交错 式相互进行．散度校正既消除了电流散度较大的伪解，还极大地提高了三维MT矢量有限元 的求解速度． 在此基础上，对二维和三维大地电磁(MT)山峰地形影响进行了数值模拟．数值模拟结果 表明，二维山峰地形对TM极化模式影响较大，对TE极化模式较小，而三维山峰地形对XY 和YX模式的视电阻率都产生都比较大，受三维山峰地形影响，在山顶处视XY和YX模式的 视电阻率都呈现明显的低阻畸变特征． 关键词大地电磁(MT)三维地形影响有限单元法数值模拟散度校正 致谢 本工作为国家自然科学基金(批准号、国土资源部重点科技攻关项目(批准 号、中国石油天然气集团公司中青年创新基金项目(批准号：2002t70112)和湖 奖学金资助． 200 主要参考文献 PE，StodtJ Aand L Two—dimensional in modeled 【I]Wannamaker Rijo topographicresponsemagnetotelludcs 131—2144． usingfiniteelement．Geophysics，1986，51(II)：2 [2】徐世浙，王庆乙，王军用边界单元法模拟二维地形对大地电磁场的影响．地球物理学 报，1992．35(1)：380—388 [3】晋光文’剥、洁，王继军地形对大地电磁资料的影响地震地质，1997，19(4)：363—369． f41王绪本，李永年，高永才．大地电磁测深二维地形影响及其校正方法研究物探化探计算技 术，1999，21(4)：3271． 【5】 张翔．带地形的大地电磁测深联合反演．石油地球物理勘探，1999。34(2)：190—196． 【6】 徐世浙历百饶，周辉等．大地电磁场三维地形影响的数值模拟．中国科学(D辑)，1997，27(1)：15～20 [71 陈伯舫，侯作中，范目华．有限差分法计算三维地形影响的电磁感应．地震学报，1998．20(5)：541—544． SeamaN．Anew forthe ofseafloor in Kand [8】Baba