起伏地形下可控源音频大地电磁三维数值模拟-地球信息科学学报.PDF

第３２卷　第２期 现　代　地　质 Ｖｏｌ３２　Ｎｏ２ ２０１８年４月 ＧＥＯＳＣＩＥＮＣＥ Ａｐｒ２０１８  ＤＯＩ：１０１９６５７／ｊｇｅｏｓｃｉｅｎｃｅ１０００－８５２７２０１８０２１９ 起伏地形下可控源音频大地电磁 三维数值模拟 １，２ ２ ３ 钟苏美 ，林昌洪 ，谢裕春 （１上海勘测设计研究院有限公司，上海　２００４３４；２中国地质大学（北京）地球物理与信息技术学院，北京　１０００８３； ３北京城建勘测设计研究院有限责任公司，北京　１００１０１） 摘要：起伏地形对可控源音频大地电磁（ＣＳＡＭＴ）响应具有强烈的影响，因此在ＣＳＡＭＴ数据处理解释时需要考虑地形。 同时，实际的地下地质情况和地表的地形情况通常比较复杂，地质结构和地形大部分情况下都是三维的。在水平地表 三维有限差分ＣＳＡＭＴ数值模拟算法的基础上，推导了利用地下交错网格采样点处的总磁场计算起伏地形下空气－地下 介质分界面处的总电场和总磁场的表达式，从而实现了起伏地形下三维ＣＳＡＭＴ数值模拟算法。在算法实现过程中，采 用伪 函数代替麦克斯韦方程中的场源项和直接计算总场的策略，避免了原有的将总场分离成背景场和二次场的策略 δ 在复杂地质条件下难以选择合适背景电阻率的问题。为了直接模拟总场，起伏地形下三维ＣＳＡＭＴ数值模拟算法给出了 新的三维正演方程的边界条件。将模拟水平地表三维异常体和三维山峰地形两个理论模型得到的响应结果与前人算法 的计算结果进行对比，验证了所实现算法的正确性和有效性。 关键词：可控源音频大地电磁法；三维地形；数值模拟；伪 函数 δ 中图分类号：Ｐ６３１　　　文献标志码：Ａ　　　文章编号：１０００－８５２７（２０１８）０２－０３９８－０８ ＴｈｒｅｅｄｉｍｅｎｓｉｏｎａｌＭｏｄｅｌｉｎｇｏｆＣｏｎｔｒｏｌｌｅｄＳｏｕｒｃｅＡｕｄｉｏ ＭａｇｎｅｔｏｔｅｌｌｕｒｉｃＤａｔａｗｉｔｈＴｏｐｏｇｒａｐｈｙ １，２ ２ ３ ＺＨＯＮＧＳｕｍｅｉ ，ＬＩＮＣｈａｎｇｈｏｎｇ，ＸＩＥＹｕｃｈｕｎ （１Ｄｅｓｉｇｎ＆ＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅ，ＳｈａｎｇｈａｉＩｎｖｅｓｔｉｇａｔｉｏｎ，Ｓｈａｎｇｈａｉ　２００４３４，Ｃｈｉｎａ； ２ＳｃｈｏｏｌｏｆＧｅｏｐｈｙｓｉｃｓａｎｄＩｎｆｏｒｍａｔｉｏｎＴｅｃｈｎｏｌｏｇｙ，ＣｈｉｎａＵｎｉｖｅｒｓｉｔｙｏｆＧｅｏｓｃｉｅｎｃｅｓ，Ｂｅｉｊｉｎｇ　１０００８３，Ｃｈｉｎａ； ３ＢｅｉｊｉｎｇＵｒｂａｎＣｏｎｓｔｒｕｃｔｉｏｎ＆ＳｕｒｖｅｙｉｎｇＤｅｓｉｇｎＲｅｓｅａｒｃｈＩｎｓｔｉｔｕｔｅＣｏ，ＬＴＤ，Ｂｅｉｊｉｎｇ　１００１０１，Ｃｈｉｎａ） Ａｂｓｔｒａｃｔ：Ｔｈｅｔｏｐｏｇｒａｐｈｙｈａｓａｓｔｒｏｎｇｅｆｆｅｃｔｏｎｔｈｅｃｏｎｔｒｏｌｌｅｄｓｏｕｒｃｅａｕｄｉｏｍａｇｎｅｔｏｔｅｌｌｕｒｉｃ（ＣＳＡＭＴ）ｄａｔａ， ｔｈｅｒｅｆｏｒｅ，ｔｈｅｔｏｐｏｇｒａｐｈｙｓｈｏｕｌｄｂｅｃｏｎｓｉｄｅｒｅｄｗｈｅｎｉｎｔｅｒｐｒｅｔｉｎｇｔｈｅＣＳＡＭＴｄａｔａ．Ｔｈｅｒｅａｌｓｕｂｓｕｒｆａｃｅｇｅｏｌｏ ｇｉｃａｌｓｔｒｕｃｔｕｒｅａｎｄｔｈｅｔｏｐｏｇｒａｐｈｙｏｆｔｈｅｓｕｒｆａｃｅａｒｅｕｓｕａｌｌｙｍｏｒｅｃｏｍｐｌｅｘａｎｄａｒｅｔｈｒｅｅｄｉｍｅｎｓｉｏｎａｌｓｐａｃｅｓｉｎ ｍｏｓｔｃａｓｅｓ．Ｂａｓｅｄｏｎｔｈｅ