考虑关断时间效应的瞬变电磁一维反演-物探与化探.PDFVIP

　第２９卷第３期 物　探　与　化　探 Ｖｏｌ．２９，Ｎｏ．３　 　２００５年６月 ＧＥＯＰＨＹＳＩＣＡＬ＆ＧＥＯＣＨＥＭＩＣＡＬＥＸＰＬＯＲＡＴＩＯＮ Ｊｕｎ．，２００５　 考虑关断时间效应的瞬变电磁一维反演 １ １ ２ 杨云见 ，王绪本 ，何展翔 （１．成都理工大学，四川 成都　６１００５９；２．东方地球物理公司，河北 涿州　０７２７５１） 摘 要：为克服因关断时间校正而带来的误差，尝试直接对包含关断时间效应的瞬变电磁资料进行反演。正演采用 将斜阶跃波离散为多个阶跃波之和的方法，反演采用经典的马奎特法。对模型试算的结果表明，采用该方法直接 对斜阶跃响应进行反演，结果比较理想，计算速度也可以接受。 关键词：瞬变电磁法；关断时间；一维反演 中图分类号：Ｐ６３１．３　　　文献标识码：Ａ　　　文章编号：１０００－８９１８（２００５）０３－０２３４－０３ 　　瞬变电磁法（ＴＥＭ）是利用接地导线或不接地 时，斜阶跃波分割为５０个小阶跃波的计算精度是足 回线向地下发送一次脉冲电磁场，在一次脉冲电磁 够的。在计算量上，由于仍是只需计算 １组汉克 场的间歇期间（断电），通过观测与研究二次涡流场 变换，再加上ｎ次正弦变换，因此与阶跃响应的计算 随时间的变化规律来探测介质电性的物探方法。由 相比，计算量增加不大。 于该方法观测的是纯二次场，可在近区观测，具有分 层能力强、精度高等特点，近年来在矿产勘查、水文 环境等领域得到了广泛的应用。但由于正演计算较 复杂、计算量大等原因，目前实际工作中瞬变电磁法 的处理解释还主要是采用以基于阶跃响应的一维反 演或是“烟圈”法等近似反演。由于实际观测的资 料均是含有斜阶跃效应，因此瞬变电磁资料的处理 解释需要作关断时间校正，关断时间校正的误差就 图１　斜阶跃波离散为阶跃波示意 直接影响着反演解释的精度。笔者针对当前使用最 广的中心回线法，尝试对观测资料不作关断时间校 正，直接对斜阶跃响应进行一维反演，以克服因关断 时间校正而带来的误差。 １　斜阶跃响应的一维正演 当前斜阶跃响应的计算方法较多，笔者采取将 斜阶跃波分割为ｎ个小电流元，每一小电流元用阶 跃波来代替斜阶跃波的方法，随着ｎ取值的增大，离 散的阶跃波就逼近斜阶跃波，如图１所示。这种方 法计算简单，且可计算任意关断波形的响应。图２ 图２　斜阶跃响应与阶跃响应 是在供电电流１０Ａ、圆回线半径２００ｍ、半空间电阻 对于阶跃瞬变响应，笔者采用数字滤波计算汉 率为１００ ·ｍ条件下，阶跃响应与 ｎ分别取５０及 Ω 克尔变换、将反富式变换转化为正弦变换的计算方 ５００，关断时间为５０ ｓ的斜阶跃响应的正演计算结 μ 法。计算分２步：先由汉克尔变换计算频率域的响 果。图２中ｎ取５０及５００的斜阶跃响应曲线几乎 应，再由正弦变换计算时间域的响应。（１）式为频 重合，说明将斜阶跃波分割为ｎ个小阶跃波的计算 率域中求解Ｈ 的汉克尔变换式，滤波系数采用参考 ｚ 方法是可行的；另一方面也说明当关断时间为５０ ｓ μ 文献［２］中的汉克 变换滤波系数；（２）式为从频率 收稿日期：２００４－１２－０５ 　３期 杨云见等：考虑关断时间效应的瞬变电磁一维反演 ·２３５· 域到时间