TEM正演响应计算的几种频时域转换方法对比-物探与化探.PDFVIP

　 第４０卷第４期 物　 探　 与　 化　 探 Ｖｏｌ．４０，Ｎｏ．４　 　 ２０１６年８月 ＧＥＯＰＨＹＳＩＣＡＬ ＆ ＧＥＯＣＨＥＭＩＣＡＬ ＥＸＰＬＯＲＡＴＩＯＮ Ａｕｇ．，２０１６　 ｄｏｉ：１０．１１７２０／ ｗｔｙｈｔ．２０１６．４．１７ － 李锋平，杨海燕，邓居智，等．ＴＥＭ正演响应计算的几种频 时域转换方法对比［Ｊ］．物探与化探，２０１６，４０（４）：７４３ ７４９．ｈｔｔｐ：／ ／ ｄｏｉ．ｏｒｇ／ １０．１１７２０／ ｗｔｙｈｔ．２０１６．４．１７ Ｌｉ ＦＰ，ＹａｎｇＨＹ，ＤｅｎｇＪＺ，ｅｔａｌ．Ｃｏｍｐａｒｉｓｏｎｏｆ ｓｅｖｅｒａｌｆｒｅｑｕｅｎｃｙ⁃ｔｉｍｅｔｒａｎｓｆｏｒｍａｔｉｏｎｍｅｔｈｏｄｓｆｏｒＴＥＭｆｏｒｗａｒｄｍｏｄｅｌｉｎｇ［Ｊ］．ＧｅｏｐｈｙｓｉｃａｌａｎｄＧｅｏｃｈｅｍｉ⁃ － ｃａｌ Ｅｘｐｌｏｒａｔｉｏｎ，２０１６，４０（４）：７４３ ７４９．ｈｔｔｐ：／ ／ ｄｏｉ．ｏｒｇ／ １０．１１７２０／ ｗｔｙｈｔ．２０１６．４．１７ ＴＥＭ正演响应计算的几种频 时域 转换方法对比 李锋平，杨海燕，邓居智，汤洪志，刘旭华，赵海娇 （东华理工大学放射性地质与勘探技术国防重点学科实验室，江西 南昌　 ３３００１３） 摘 要：在瞬变电磁的一维正演响应模拟中，常用的方法是先在频率域中求解，之后将结果转换到时间域。 但该方 － 法在晚期的计算精度通常不高，因此，使用５种频 时域转换方法（正、余弦变换的数值滤波算法，Ｇ⁃Ｓ逆拉普拉斯 算法、正、余弦变换的折线逼近法）进行了计算，与解析解对比，得出余弦变换的数值滤波算法在晚期计算中精度最 高的结论，并对这五种转换方法产生误差问题的原因进行了讨论和分析。 本研究有利于瞬变电磁一维正演响应的 高精度计算，使其在多维计算中得到更好的应用。 关键词：瞬变电磁；频 时域转换；数值滤波；Ｇ⁃Ｓ逆拉普拉斯算法；折线逼近 － － － 中图分类号：Ｐ６３１　 　 　 文献标识码：Ａ　 　 　 文章编号：１０００ ８９１８（２０１６）０４ ０７４３ ０７ 　 　 瞬变电磁法正演响应的计算一般都是先求得频 法，比如费龙算法、折线逼近算法和数值滤波算 率域的解，然后由频率域转换到时间域得到瞬态 法［１３ １４］。 国内长期致力于研究瞬变电磁法的专家 ［１ ５］ 朴化荣教授对 Ｇａｖｅｒ⁃Ｓｔｅｈｆｅｓｔ 进行了改进，提出Ｇａ⁃ 解 。 目前，可用于频 时域转换的方法较多，且 ［１５］ 各具特色。 傅立叶变换是最早应用于其中的方法， ｖｅｒ⁃Ｓｔｅｈｆｅｓｔ逆拉普拉斯算法 ；陈乐寿、王光锷于 但傅立叶变换过程中要求频率的采样点数必须满足 １９９１年用数字滤波方法计算了瞬变电磁的响 ［１６］ 采样定理，这样频率采样点会增多，从而造成计算量 应 ；李吉松、朴化荣对考夫曼 Ａ．Ａ．和凯勒 Ｇ．Ｖ．提 ［６］ 加大，计算速度减慢 。 Ａｎｄｅｒｓｏｎ于１９７４年应用改