对数与积分计算Hz定义中心回线瞬变电磁法全区视电阻率.pdf

对数与积分计算Hz定义中心回线 瞬变电磁法全区视电阻率 杨云见“王绪本’何展翔’ (1成都理工大学6100592冻方地球物理公司072751) 摘要：针对里羔衰减迅速、时间及幅值变化范围大的特点，探讨在对数域内插值、积分计算地，再由Hz df 来定义中心回线瞬变电磁法全区视电阻率。模型试算结果表明：只要选取合适的外推时问长度，该方法可得 到精度较高、连续的视电阻率定义，在延时较短的情况下也可得到较好的结果。 关键词：对数域插值、积分中心回线瞬变电磁法垂直磁场全区视电阻率 引言 瞬变电磁法(TEM)是利用接地导线或不接地回线向地下发送一次脉冲电磁场，在一次场断 电后，通过观测二次涡流场随时间的变化规律来探测介质电性的一种物探方法。由于观测的 是纯二次场，可在近区观测，具有旁侧影响小、异常响应强等特点。在各种观测装置中，中 心回线法是与探测目标耦合最好的方法之一，也是最常用的装置之一…，已广泛应用于矿产、 水、工、环等勘探领域。目前中心回线瞬变电磁法观测的是感应电动势(也即磁场随时问的 变化率(主当)，以下旦当用符号膏，表示)，由于技术上的限制而难以直接观测磁场H：…。 。 df df 由日，定义的视电阻率通常分早期和晚期，定义全区视电阻率存在多解甚至在无解的情况，难 以在全时间段给出连续的电性变化信息。而由磁场H：定义视电阻率不存在多解的问题，在全 区段上解唯一，故由日，计算出H：来定义全区视电阻率将是一条有效的途径。但由于日：的 观测时间长度有限，由日．计算H：总存在误差。针对全区视电阻率定义问题，电磁勘探工作 者们提出不少解决方法，如：苏朱刘针对H，提出了虚拟视电阻率的定义方法陀1，解决了全区 定义视电阻率的连续性问题，但所定义的虚拟视电阻率与早、晚期定义的电阻率有差异；黄 力军采用多项式拟和电动势、用解析积分的方法来求取Hzbl，从而唯一的定义全区视电阻率， 该方法计算的H：与理论计算得到的H；基本吻合。本文针对H，衰减迅速、时间及幅值变化范 围大的特点，探讨在对数域插值、积分计算啦来定义视电阻率，以期得到精度较高、全时间 段连续的全区视电阻率。 一、地的计算 中心回线法瞬变电磁观测的是垂直线圈接收的感应电动势V(t)，将v(o对接收线圈面积作 归一化并除于磁导率‰即可得到日，。依据电磁学原理，相应磁场强度Hz(t)为： 鼍’ H：=JH：dt (1) 结合中心回线瞬变电磁法工作特点，可知(I)式的积分具有如下几个特点 287 ①厅，的时间、幅值变化范围大，可以达到几个数量级，时间的变化范围经常在n·us。 n·lOOms、幅值变化范围经常在n-100A／(m·s)一noe．6A／(m·s)。 ②日，变化快，在均匀半空间中，Ⅳ，以时间t的．5／3次幂衰减。 ③观测时间的变化范围虽然很大，但确是有限长的，而(1)式的积分要求却是无限长。 这些特点决定了(1)式积分的关键是解决精度问题。针对该式积分的特点本文从两个方 面来考虑：①调整日，时间、幅值的变化范围，使之处于同一个量级上；②减缓日，的变化趋 势，提高插值、外推的精度。通过这两个方面的考虑可知：对观测时间、幅值取对数，在对 数域中进行插值、积分将是一个有效的方法。 对(1)式中H，的观测时间、幅值取对数，设： U=lg(H一)，．V=lg(t) 则： J10“10”dv (2) H：=JH：dr=rio“dt=lnl0 ，t Ig(t) (2)式即将(1)式算术域内的积分转到了对数域。对该式进行数值积分，·可取为最后一 个时间点的对数值加上系数a，即相当于沿时间往后外推了104倍。本文插值采用三次样条函 数插值法。 二、H：的全区视电阻率定义 将非均匀介质中观测的结果按均匀半空间的公式计算电阻率即为视电阻率，它是