均匀半空间视电阻率MATLAB代码,均匀半空间中心回线瞬变电磁法正演研究

第 12 卷 第 9 期 2012 年 3 月 1671—1815( 2012) 09-2128-04 科 学 技 术 与 工 程 Science Technology and Engineering Vol. 12 No. 9 Mar． 2012  2012 Sci. Tech. Engrg. 均匀半空间中心回线瞬变电磁法正演研究 晏冲为 李文尧 ( 昆明理工大学 国土资源与工程学院，昆明 650093) 摘 要 通过对均匀半空间中心回线装置瞬变电磁法表达式分析，可以为瞬变电磁法的反演提供理论正演数据、可以为瞬变电磁法的反演提供一定的方法，同时还可以为野外生产设计与室内数据筛选与剔除提供一定的参考。通过一个反演实例说明在晚期微小的二次电位误差将对视电阻率的计算产生很大的影响。编制了 matlab 程序计算各种工作装置和视电阻率理论正演二次电位衰减到噪声影响程度高的准确时间。对工程设计、评价视电阻计算的准确度及剔除受噪声影响严重的数据有一定的参考性。 关键词 中心回线瞬变电磁法 二次电位 衰减时间中图法分类号 P631. 34; 文献标志码 A 2012 年 1 月 4 日收到 地质调查项目( 1212011121273) 资助第一作者简介: 晏冲为( 1984—) ，男，云南富源人，硕士研究生，研究方向: 瞬变电磁法理论与仪器。E-mail:89534861@ qq. com。 瞬变电磁法( TEM) ［1，2］中用接收线圈接收的感应电压随时间呈早期信号幅值高且衰变速度快，晚期信号很弱衰减速度慢的特点，均匀半空间中心回线瞬变电磁法的感应电压表达式可以为理论研究提供无干扰的正演数据，同时还能为反应的研究提供一定的方法，均匀半空间的瞬变电磁场与地下电阻率之间存在着复杂的隐函数关系，无法获得全期视电阻率的解析表达式，目前通常使用各种近似定义和精确定义的方法来计算早期与晚期视电阻率或全期视电阻率［3］，全期视电阻率的计算很多学者提出了各种成熟的计算方法［3—7］，其中陈清礼［7］基于在 0. 08 ms ＜ t ＜1 s 的范围内，感应电动势随电阻率的增大而单调下降的特性提出了全期视电阻率的二分搜索算法，搜索出理论感应电动势与实际观测的电动势相符的视电阻率即是基于正演研究提供反演的一种方法的例子，现通过一个实例来讨论晚期干扰噪声对视电阻率计算的影响，并通过 mat-lab 编写了一个计算各种工作装置和视电阻率理论正演二次电位衰减到自主设定的噪声影响程度高的准确时间的程序，对于前期工程设计、后期数据 处理解释时剔除不合理晚期数据，对提高解释准确度有一定的参考意义。 1 均匀半空间下感应电动势的解析表达 中心回线或重叠回线方式是时间域瞬变电磁法大多采用的工作方法，1983 年美国地球物理学家P. Raab 和 F. Frisehknecht［8］推导出了中心回线装置的感应电压表达式为 V( t) = qIπ3/2ρ L3 ［3Φ( z) － ( 3z + 2z3) Φ · ( z) ］u( t) ( 1) 其中 q 为接收线圈的有效面积( m2) ，L 为发送回线边长( m) ，I 为发送电流强度( A) ，它为满足 I( t) = 0， t ＜ 0 I0， t ≥{ 0 ( 2) 的阶跃电流。Φ( z) 为概率积分， Φ( z) = ∫ z 0 2槡 π exp( － t2) dt ( 3) Φ · ( z) = 2 槡 π exp( － z2) ( 4) z 为瞬变场参数，τ 为扩散参数( m) z = 2 槡 πL/τ ( 5) τ = 2π 2ρ