TEM中心回线法计算考虑关断时间的全区视电阻率.pdfVIP

第32卷第6期 物 探 与 化 探 Vo1．32．No．6 2008年 12月 GEOPHYSICAL GEOCHEMICALEXPLORATION Dee．．2008 TEM 中心回线法计算考虑关断时间的全区视 电阻率 杨 生 (有色金属矿产地质调查中心 北京中色物探有限公司，北京 100045) 摘 要：详细推导了TEM法中心回线装置发射电流为斜阶跃波形条件下全区视电阻率迭代反演计算方法。对正演 过程中对关断时间的积分采用了解析表示。对反演过程 中求导运算，通过变量变换，将对电阻率的求导变换为对 时问的求导，从而也给出简单的解析关系。避开了数值积分运算和差分替代导数的近似，提高了计算精度和速度。 为克服迭代反演过程中要多次计算概率积分的耗时问题 ，采用标准数据库的方法一劳永逸地加以解决，也大大提 高了计算精度和速度。对实测资料进行了计算，并通过算例对一些问题进行了讨论。 关键词：瞬变电磁法 ；中t9回线装置；全区视电阻率计算方法；关断时间 中图分类号：P631 文献标识码：A 文章编号：1000—8918(2008)06—0647—05 瞬变电磁法 (简称 TEM)因测激励场产生的二 延迟时问；err(z)为概率积分 (也称误差函数)。 次涡流场在地表的响应，对地下低阻体十分敏感，所 式(1)中，观测量B和P是隐函数关系，所以通 以是金属矿勘探的首选方法之一。与所有其他 电阻 常的作法是用 和 一。。时的渐近表达式来计算 率方法类似，瞬变电磁法的视 电阻率导入也是根据 视电阻率，这就出现了远区(一0时)和近区(一o。 均匀半空间的电阻率 (P)和测量参数(m)间的函数 时)的公式： 关系P=-厂(m)，将其推广到非均匀半空间(一维、二 P。l = (1／31)·B ·( 竹) ， (2) 维和三维地电条件)时，计算出的P称为视 电阻率， P l = (1／．rr)(i_t／t) (，·S／20B) 。(3) 它定性地反应了地下介质的电性分布特征。 式 (2)称为远区公式，式 (3)称为近区公式。要注意 关于全区视 电阻率的计算问题，已有较深入的 的是在瞬变电磁法 中，“区”和 “期”的称谓经常混 研究 ¨ ，而在本文 中，笔者在考虑关断时间的同 用，“区”是从场的特征描述的，“期”是从时问的角 时，对计算过程中的积分运算和求导运算尽量采用 度描述的，两者一致，远区就等于早期，近区就等于 了解析方法，不仅避开 了数值积分和差分对计算精 晚期 ，但在讨论瞬变延时曲线特征时要用到 “早期 度的影响，同时也大大提高了运算速度。笔者对计 取样道”和 “晚期取样道”这样的述语，而 “早期取样 算过程中的各个环节作了详尽的交待，力求实用和 道”和 “早期”，“晚期取样道”和 “晚期”是不同的2 可操作 。 个概念 ，为了避免叙述上的混乱，文 中采用 “区”的 称谓 。 1 问题的提出 数值计算表明，在Z5时用式 (2)近似式 (1) 目前的瞬变电磁法，观测的是瞬变磁场对时间 和当Z0．29时用式(3)近似式(1)，误差在容许的 的导数 (dB／dt)，以B表示，中心回线装置在均匀半 范围内 。在 中心回线法 中，Z5的情况很少碰 空间下有E4]： 到，而主要用的是式 (3)的近区公式，但在实际工作