1255--回线源瞬变电磁测深快速反演方法.docxVIP

第 27卷　第 4期 物探化探计算技术 　　 2005年 11月 　 文章编号: 1001—1749 (2005) 04—0289—03 　 　 　 回线源瞬变电磁测深快速反演方法 陆桂福, 张　强, 黄力军 (中国地质科学院 地球物理地球化学勘查研究所, 河北 廊坊 065000) 摘　要: 回线源瞬变电磁测深法 (简称瞬变电磁法或 T EM 法)是以不接地回线为场源,在回线 内、外测量电磁场的一种电磁测深方法,已经广泛用于固体矿产、工程地质和水文地质调查工作 之中 [ 1, 2 ] 。目前,回线源瞬变电磁法主要采用中心回线和重叠回线两种测量装置,测量回线内垂 直磁场产生的感应电动势。“烟圈 ”反演方法是瞬变电磁反演的主要计算方法,“烟圈 ”反演利用 地表瞬变电磁响应与某时刻电流源镜像等效原理,计算勘探深度和对应电阻率。这里给出了回 线源瞬变电磁测深一维“烟圈 ”反演解释方法的基本原理、计算方法和应用实例。结果表明,一 维“烟圈 ”反演解释结果纵向分辨率较高,能反映出地下电性变化情况,可以用于瞬变电磁测深 资料解释工作。 关键词: 回线源; 瞬变电磁测深; 烟圈反演理论; 快速反演方法 中图分类号: P63113+ 25　　文献标识码: A 0　前言 目前,回线源瞬变电磁法主要测量回线内垂直 磁场产生的感应电动势,资料解释大多采用多道电 位 (感应电动势)剖面。由于多道电位剖面解释还 处于定性解释阶段,因而影响瞬变电磁法的实际应 用效果。针对这一问题,我们在实际解释工作中, 图 1　均匀大地磁场与回线源 TEM装置示意图 F ig. 1　Sketch m ap of the magnetic fie ld and 将实测感应电动势转换成磁场 [ 3, 4 ] ,利用磁场计算 [ 3, 4 ] ,利用磁场计算 sou rce w ire on a homogeneous ha lf sp ace 全区视电阻率,采用“烟圈 ”反演方法进行一维反 1. 1　计算视电阻率 演。 在电阻率为 ρ,均匀导电半无限空间表面上, 1　反演原理 [ 5 ] 载有阶跃电流水平电偶极子产生的时间域磁场,其 解析表达式为 烟圈反演理论认为,感应涡流场在地表引起的 Idly 　hz ( x, y) = 3 · 4πr 瞬变电磁响应为地下各个涡流层的总效应,这种效 应可近似地用向下传播地电流环等效。这些电流 　　　 (1 - 3 2 ) e rf ( u) + 2u 3 - u2 ·e u π (1) 环好像是发射源吹出的“烟圈 ”,其形状与发射源 当 u很大时 (早期) ,式 (1)可以写成 相同,随着采样时间延长而向下延伸,地表瞬变电 2 ) / (4π·r3 ) 　hz ( t) = Idly (1 - 3 /2u 当 u很小时 (晚期) ,式 (1)可以写成 磁响应可以与某时刻电流源镜像等效,利用这个等 效原理计算勘探深度和对应电阻率,见图 1。 收稿日期: 2004 - 10 - 02 ·290· 物探化探计算技术 27卷 　hz ( t) = Idly (8u 3 / π) · 3 / π) · 　　　(1 /15 - u 2 /35 + u4 /126) / (4π·r3 ) 2 /35 + u4 /126) / (4π·r3 ) 式中 　I是阶跃电流的振幅; dl是偶极长度; t是 采样时间; u = [μ 2 / (4ρt) ]1 /2 ; r = ( x2 + y2 ) 1 /2 。 0 r采样时间; u = [μ 2 / (4ρt) ]1 /2 ; r = ( x2 + y2 ) 1 /2 。 N 　Ψ j = 6 j = 6 m kI ( dj ) /N (3) k ( tj ) - H H k =1 式中 　Hk ( tj )是第 k个接收点、第 j次采样的实测 m 的阶跃响应; Hk ( dj )是像在源下方深度 dj 时,在第 I k个接收点处引起的磁场; N 为总采样点数,这里 根据实际发射源的形状对式 (1)积分,可算出 相应装置的磁场函数,利用实测的磁场与磁场函数 dj dj - 1 0。 1. 3　确定电阻率与深度之间的函数关系 反演拟合,可以计算出对应时间的视电阻率。 在得到离散镜像源位置 dj = d ( tj )以后,需要 1. 2　利用实测数据计算镜像源深度 确定电阻率和勘探深度 。在非均匀半空间情况下, 在 TEM测量中, 接收线圈可以接收一个或几 视电阻率不能直接计算勘探深度,需将视电阻率带 个不同位置发射源产生的感应电动势。计算图 2 入深度函数, 利用深度函数对时间的导数计算速 所示接地源测量装置的反演,这里假设发射源沿着 度。地下电流“烟圈 ”扩散速度直接取决于地下电 剖面