大地电磁二维自适应地形有限元正演模拟∫Ω-中国地球物理学会地球.PDFVIP

第 ３２卷　第 ３期 地　震　地　质 Ｖｏｌ．３２，Ｎｏ．３ ２０１０年９月 ＳＥＩＳＭＯＬＯＧＹＡＮＤＧＥＯＬＯＧＹ Ｓｅｐｔ．，２０１０ ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０２５３－４９６７．２０１０．０３．００４ 大地电磁二维自适应地形有限元正演模拟 刘　云　王绪本 （成都理工大学信息工程学院，成都　６１００５９） 摘　要　　复杂地形对大地电磁测深法（ＭＴ）的影响非常大。文中在吸收和修正前人工作基础上，提 出了一种自适应地形四边形网格剖分、单元内电磁场双二次插值数值模拟方法，能更好地模拟出复杂 地形和场值分布，并推导出复杂地形条件下有限元数值模拟算法、单元辅助场的计算和ＴＭ模式视电 阻率的局部定义法。模型计算表明，该数值模拟方法计算速度快，有较高的计算精度，能很好地模拟 ＭＴ复杂地形影响和异常体影响。 关键词　　大地电磁测深　自适应地形　四边形网格　有限元 中图分类号：Ｐ６３１ 文献标识码：Ａ 文章编号：０２５３－４９６７（２０１０）０３－０３８２－１０ ０　引言 大地电磁测深法（ＭＴ）已广泛应用于油气勘探、海洋及深部地球物理探测。在山区开展 ＭＴ 工作时，复杂地形条件对电磁场的影响非常大，给资料的处理和解释工作带来了很大困难，因此 有必要研究带地形的ＭＴ数值模拟方法。 在二维地电模型条件下，国内外学者采用过矩形单元、双线性插值和双二次插值有限单元 ＭＴ数值模拟（Ｐｈｉｌｉｐｅｔａｌ．，１９８６；徐世浙，１９９４；史明娟等，１９９７；阮百尧，１９９７；王绪本等， １９９９；陈小斌，２０００，２００５；马为等，２００８）。但是，在野外复杂地形环境中，场的变化往往随着起 伏地形和复杂地质体非线性变化，因此采用矩形网格、线性插值都很难精确模拟出野外实际地 形影响和电磁场的异常形态。本文基于史明娟等（１９９７）矩形单元、双二次插值有限单元法基 础上，对部分系数矩阵做了一些改动，并提出了自适应地形四边形网格剖分、电磁场双二次插值 二维 ＭＴ有限单元数值模拟，文中算例模型的网格采用陈小斌等（２０００）所示的网格。 １　二维 ＭＴ变分问题 如图１ａ所示，假设三维地质体电（磁）参数在地质构造方向（ｘ方向）上无变化，这时可将三 维模型问题转化为二维模型问题处理。此时，与二维 ＭＴ边值问题相应的变分问题（徐世浙， １９９４；史明娟等，１９９７）是  １ ２ １ ２ １ ２ Ｆ（ｕ）＝ ｄ ＋ ｋｕｄ （ｕ） － ｕ Ω τ Γ [ τ λ ]  ∫ ∫２  Ω ２ ２ ＣＤ  ｕ ＝１ （１）  ＡＢ   Ｆ（ｕ）＝０ δ 〔收稿日期〕　２０１０－０１－０６收稿，２０１０－０７－０６改回。 〔基金项目〕　国家８６３滚动课题 “起伏地形二三维电磁法正反演解释工作站”（２００６ＡＡ０６Ｚ１１０）资助。 ３期 刘　云等：大地电磁二维自适应地形有限元正演模拟 ３８３ 式（１）中，ｕ为Ｅ 或Ｈ， 为二维求解区域；对于ＴＥ模式（Ｅ 型）：上边界 ＡＢ处离地面足够 Ω ｘ ｘ