瞬变电磁法拟地震偏移成像研究(李貅,薛国强著)思维导图
瞬变电磁法正反演问题研究进展
第23卷 第4期地 球 物 理 学 进 展V ol.23 N o.42008年8月(页码:1165~1172)P ROG RESS IN G EOP HY SICSA ug. 2008瞬变电磁法正反演问题研究进展薛国强1, 李 貅2, 底青云1(1.中国科学院地质与地球物理研究所,北京100029; 2.长安大学地质工程与测绘工程学院,西安710054)摘 要 对瞬变电磁法的方法发展概况和仪器研制状况做出了综述性评价.对瞬变电磁法正反演问题的研究成果进行了系统总结.目前的数值模拟正演方法主要有一维滤波系数法,三维积分方程法,二维,三维有限差分法,2.5维有限元法等,主要的反演方法有:一维浮动薄板解释法,人机对话自动反演法,烟圈理论解释法,神经网络反演法,成像类反演等,论述了瞬变电磁法各种计算方法的特点.瞬变电磁法的正反演发展趋势主要是研究三维正反演的计算方法和目标体成像系统.关键词 瞬变电磁法,正演问题,反演问题中图分类号 P631 文献标识码 A 文章编号 1004-2903(2008)04-1165-08Research progress in TEM forward modelingand inversion calculationXU E Guo -Qiang 1, LI Xiu 2, DI Qing -Yun 1(1.I nstitute of Geolog y and Geop hysics ,Chinese A cad emy of science ,Beij ing ,100029,China;2.S chool of Ge ology and S ur v ey Eng ineering ,Chang 'a n Univ ersity ,X i 'a n 710054,China)Abstract We g ive a g ener al rev iew o f recent r esear ch prog ress abo ut the T r ansient Electro magnetic M ethod (T EM )and instrument .T he achievement o f T EM for wa rd modeling and inv ersion hav e been summar ized .A t pr esent,the for war d numer ical calculation metho ds include 1D dig ital filter method,3D integ ral equat ion metho d,2D and 3D differ -ence met ho d,2.5D infinite element met ho d.T he main 1D inv erse met ho ds include float plate method,auto -inv erse metho d,smoking theor y,artificial neural netw or k,imag ing inver sion.We present the featur e o f ev ery numerical calcu -latio n metho d.T he developing directio n of fo rw ard and inver se st udies mainly is the 3D calculat ion method and the targ et imag ing sysy tem.Keywords tr ansient electr omag netic metho d,forw ard pro blem,inver sion problem收稿日期 2008-03-10; 修回日期 2008-06-20.基金项目 国家自然科学基金重点项目(50539080)和国家自然科学基金项目(40774066)联合资助.作者简介 薛国强,男,1966年生,1989年于西安地质学院获学士学位,2002年于长安大学获硕士学位,于2005年于西安交通大学获博士学位.研究方向为电磁探测理论与应用(E -mail:ppxueguoqiang@ )0 引 言电磁场理论的应用已经遍及地学、生命科学、医学、空间科学、信息科学等几乎所有的技术科学领域,同时这些工程技术领域对电磁理论研究也不断地提出各种新的要求.电磁法勘探是基于研究电磁波在导电介质中传播特性,从而达到研究地下地质体赋存特性的目的.通过天然或人工场源在大地中激励的交变电磁场,研究电磁场的空间和时间分布,分析观测到的电磁场信号,得到地下目标体的电性分布特征的一种地球物理方法.瞬变电磁测深法(Transient electr omagneticm ethod,简称T EM)是电磁法勘探中应用较广的一种,是近年来在工程地质勘察中普遍应用的时间域电磁探测方法.它是利用阶跃波或其它脉冲电流场源激励,在大地产生过渡过程场,断电瞬间在大地中形成涡旋交变电磁场,测量这种由地下介质产生的二次感应电磁场随时间变化的衰减特性,从测量得到的异常信号中分析出地下不均匀体的导电性能和位置,从而推断矿体、工程基础、地下水、地质灾害、工程病态等地下目标体的分布性态.该技术具有灵敏度高、分辨率强、探测深度大、灵活多变适应性强地球物理学进展23卷以及轻便、快速、廉价诸多优点,近年来发展十分迅猛,应用前景十分广阔.目前,瞬变电磁法已经成为地球物理探测领域内的重要方法之一.已广泛应用于水利、交通、城建、环保、考古等部门.成功地解决了大量实际问题[1,2].但是应用领域的问题越来越多并且越来越难,发展瞬变电磁法精细探测技术是一次机遇和挑战[3~5].本文详细列举了瞬变电磁法各种正反演方法,并对方法进行了回顾和展望.1方法概述20世纪30年代,最先提出利用电流脉冲激发供电偶极形成时域电磁场的是美国科学家,当时利用不同电导率地层界面电磁波的反射与地震反射波信号的相似性,进行了大量的实验和比较.最早提瞬变电磁法工作方法的前苏联科学家,当时采用的是远区工作模式[7,8].到了20世纪50~60年代,前苏联科学家成功地完成了瞬变电磁法的一维正、反演,建立了瞬变电磁法的解释理论和野外工作方法之后,瞬变电磁法才开始进入实用阶段.20世纪60年代以后,当意识到时间域电磁测深法可以利用远远小于期望探测深度的收发距时,该方法有了一个快速发展.随之,/短偏移0、/晚期0、/近区0、等技术研究迅速发展起来.美国等西方国家在20世纪70~ 80年代之间,短偏移法一直处于研究和试验阶段,未被广泛运用,而长偏移法已得到了应用,特别是在地热调查和地壳结构的深部调查中.随后一些专家对瞬变电磁法的一维正反演及方法技术进行了大量研究[9~11].20世纪80年代后随着计算机技术的发展,欧美学者在二三维正演模拟技术方面,发表了大量的论文[12~17].与此同时,前苏联学者提出电磁波拟地震波的偏移方法,吸取了/偏移成像0的广义概念,在电磁法中确定了正则偏移和解析法偏移两种方法.80年代末,从激发极化现象理论出发,研究了时间域瞬变电磁法的激电效应特征及影响,解释了瞬变电磁法晚期电磁响应的变号现象,并对三维极化体的瞬变电磁响应特征进行了数值计算.[18~20].在我国,对瞬变电磁法的研究始于20世纪70年代,朴化荣、曾孝箴等人,他们将脉冲式航电仪用于地质填图和找矿中;方文藻、李貅等将大回线源瞬变电磁测深法广泛用于地热和地下水调查、工程调查和地质灾害调查,又将瞬变电磁法用于大地电磁测深曲线的静校正,取得了良好的效果.蒋邦远等将瞬变电磁法用于普查勘探良导金属矿,随后又研制出了DCM-1型电磁脉冲瞬变系统;牛之琏等将T EM法用于金属矿勘探上,取得了明鲜的效果,并与智通研究所合作研制并生产了SD-1型智能化瞬变电磁仪,白云仪器厂在此基础上研制了M SD-1, BYF5M SD1瞬变电磁系统;中国地质科学院地球物理地球化学勘察研究所研制了WDC-1,WDC-2瞬变电磁仪器,后又研制了IGGETEM-20瞬变电磁系统;西安强源物探研究所研制了LC,EMRS-1, EM RS-2瞬变电磁仪;北京矿产地质研究所王庆乙教授研制了T EMS-3S瞬变电磁仪器;吉林大学林君教授研制了ATEM-2瞬变电磁仪器,重庆奔腾数控技术研究所研制了WT EM系统;但是目前国内仪器的稳定性,重复性都不如国外仪器好.国内主要的进口仪器有加拿大GEONICS公司生产的PRO-TEM系统(PROTEM-37,47,57,67);PH OEN IX 公司生产的V6,V8系统;美国ZONGE公司生产的GDP-32系统等.以前,瞬变电磁法只局限于金属矿勘探,1992年以后随着仪器的智能化与数字化,瞬变电磁法开始步入工程、环境、灾害地质调查中,如探测地下采空区,陷落柱等煤田灾害,划分地下断层、寻找地下水,金属矿产勘探、石油、煤炭等非金属矿产调查、工程场地地质勘察、隧道超前地质预报等领域.取得了良好的效果.瞬变电磁法以其独特的优点广泛用于资源勘探和工程勘察中.但针对一些具体的精细探测问题如:高速公路和铁路建设中的隧道超前地质预报精细探测,大型重要古墓的墓室结构精细探测,大型水坝隐患精细探测等,常用的瞬变电磁法的分辨率受到限制,探测效果受到影响.从目前看,在理论研究、仪器研制方面处于初级阶段.虽已解决了一维正反演问题,但在二、三维的研究成果还未达到应用程度.传统的瞬变电磁探测方法对地下目标体的评价精度低,一般采用二次衰减曲线和由此算得和的视电阻率值及视纵向电导值作为解释参数,用一维的计算公式得到深度)视电阻率,或者深度)视电导率二维断面图来进行解释.多数仍处于定性和/看图识字0的阶段.对瞬变电磁法测深资料定量解释还局限于单点一维反演,很多情况下是靠解释人员的工作经验及地质先验知识来对测深结果做出判断,人为性较大.随着探测分辨率和精度要求的提高,在原有解释方法理论基础上对正反演问题做更深入地研究,探索新的反演理论并将其进行系统化,建立系统的T EM解释正反演理论,同时使瞬变电磁法的解释11664期薛国强,等:瞬变电磁法正反演问题研究进展向三维方向迈进,使这一方法能更好地解决一些高难的精细探测问题.当然,加强仪器研制,观测方式的改进,微信号拾取方法研究等方面也要加强.2 正演问题研究正问题是根据给定的地球模型求解地球物理电磁场理解值,反问题是根据实际测量的地球物理电磁场数据定性或者定量解释出地球内部的结构的过程.2.1 一维正演问题对瞬变电磁法一维情况下的正演计算大多采用先在频率域进行讨论,得到层状介质下的电磁响应表达式,然后把讨论结果转换到时间域来.从频率域到时间域响应的转换,可以采用的方法有:GAVER -STEH FEST 逆拉氏变换方法、延迟谱方法、线性数字滤波方法、折线化正余弦变换法.线性数字滤波技术是将汉克尔变换转换成卷积形式,离散化后形成数字滤波器.滤波器的系数由已知变换对求出[21].这一变换主要涉及H ankel 变换,H ankel 变换方法有两种,一种采用线性数字滤波技术,利用H ankel 系数反积分方程离散化,并把积分形式变成求和形式.(流程图如图1)计算精度与系数的个数多少有关系.另外一种办法是把足够长的积分区间分成两部分,求出贝塞尔函数的两个零点,在各区间内采用高斯积分求积,然后求和,这一算法精度较高,但是效率较低,图1 一维正演数值计算流程图Fig.1 F lo w chart o f 1D fo rw ard calcalat ion折线化正余弦变换法是利用正余弦函数的导数性质和分部积分法则将积分转换成为对核函数导数的的正余弦变换.对导数进行差商近似,将积分区间分段,并在每个段内用折线来逼近积分核函数,使核函数的二次导数变成一系列D 函数的和.该算法对缓变函数效果较好.2.2 二维正演问题二维数值计算多采用有限差分法进行.用两个无限长直导线近似作为发送回线源,可用均匀半空间的解析解在t >0时刻将源转化为初始条件加入.从反映电磁场基本规律的麦克斯韦方程组出发,导出时域电场的齐次扩散方程,对所研究的空间区域作差分离散,利用准静态近似处理空中边界,在地-空界面向上延拓一个网格,采用五点差分显式格式,在每一时间步计算网格空间各点的场量,然后进行时间的逐步递推,就能直接模拟电磁波的传播及其与地质体的相互作用过程,使电磁场的时域特性被直接反映出来[22].(流程图如图2)图2 二维正演有限差分法程序流程图F ig.2 Flow char t of 2D farw ard ca lculatio nusing the difference method由于源的处理和边界条件的确定是建立在在水平均匀半空间基础上的,目前的时域微分方程算法尚存在着不足,无法分析TEM 场对浅异常体的响1167地球物理学进展23卷应;无法有效地分析地形对T EM的影响,地形与异常体的相互作用还不甚清楚;由于对源的特殊处理,使得计算结果对浅部地质体的反映不佳.有限元法对频率域计算比较有利.因为它可以很方便地处理急剧变化的和倾斜的电导率分界面和地形等问题.由于频域电磁勘探中的二维定解问题相对简单,且用有限元法求解十分方便,因此,先在频域中用有限元法求解,然后变换到时间域,是瞬变电磁二维正演的一种有效途径[21].2.32.5维正演问题2.5维时间域电磁场数值模拟问题是目前国际上尚未妥善解决的计算地球物理疑难问题之一.我国从20世纪80年代开始着手研究2.5维电磁场的数值模拟,采用有限元法对时间域[23]和频率域电磁场的2.5维响应进行正演数值模拟,采用快速松驰算法实现三维源二维地质结构的CSAMT数值计算[24],尽管占用计算机资源较大,但是确实给多维反演研究开辟了方向.为避免过于庞杂的三维正演计算量,可以首先在频率域计算电磁场响应值,然后再把计算结果变到时间域采用三角形网格代替传统的矩形网格算法,导出了中心回线瞬变电磁2.5维二次场(纯异常)的有限单元计算公式.2.4三维正演问题3D反演是提高TEM资料解释效果的必由之路,3D模型正演又是反演的基础,因此,改进与完善TEM法3D模型正演,对进一步提高TEM资料解释水平和方法的应用效果具有理论和现实意义.目前国内还没有应用较好的3D正演软件.三维瞬变电磁场正演计算只能用数值方法,主要方法有三类:有限差分法,有限元法和积分方程法.前两种方法要求对所计算的全部区域进行离散化,所占用的计算机容量较大.后者只要对异常区域进行离散化,存贮量大大减小.积分方程是最早实现3D场模拟的数值算法[19,25,26],因为积分方程法只需要计算小体积异常区的场,不必计算整个区域的场.这一点使积分方程算法在数值计算的早期,比微分方程法具有更大的优越性.但是积分方程法求解要遇到某些更困难的数学问题,仅适合模拟简单模型.在计算比较复杂的模型,如层状大地中的3D异常体时,往往不能用直接时域积分方程求解,而要由频域积分方程的计算结果经傅立叶变换到时域.尽管几乎在所有情况下,电磁场的频率域特性和时间域特性之间可以通过傅氏变换一一对应,然而也有一些情况,除了转换中的计算精度问题以外,两者之间有着很微妙的差别.要模拟复杂的地质构造,还要依赖于时域微分方程的方法.和积分方程算法不同,在开放域的地球物理问题中,微分方程算法要解决边界条件和源的处理问题.瞬变电磁三维直接时域有限差分法数模拟已经有报道[27].而且还进一步地通过设定等效位移电流,使原适合波动场的时域有限差分算法(FDT D)能够适用于扩散场问题.直接时域算法的引入,展现了电磁场在地下随时间传播的全过程,直接模拟电磁波与地下异常体的相互作用,使TEM 场的时间特性被直接反映出来,从而给复杂的物理过程描绘出清晰的物理图象.直接在时间域中求解和从频率域中转换,两者的计算量基本相当,前者精度较高,尤其是晚期,但是计算复杂.后者由于频率域的研究已有很好的基础,所以做起来简便.也可以按照二维时间域有限差分的思路进行三维有限差分计算[27],即把三维源问题做特殊处理后,作为初始条件加入迭代方程,然后进行计算.3反演问题研究3.1浮动薄板解释法[28~34]它是一种根据视纵向电导曲线的特征值直观地划分地层的近似解释方法.因此称为/视纵向电导解释法0,也有人把该法形象地称为/浮动薄板解释法0.水平薄板模型是瞬变电磁场正演计算中唯一能用初等函数解析表示的地电模型.随着时间的推移,瞬变电磁场向地层深处传播.因此,可将每个瞬间观测到的电磁信号等效为某一/浮动0导电薄层产生的场,从而直接把观测的垂直磁场分量时间导数矩阵转换为电导率)深度值.根据电磁理论,我们可用一导电平面来代替地下均匀介质,然后用镜像法可以方便地求出空间任一点的感应磁场.可形象地理解为:随时间的增减,等效导电平面以速度1L0R上下/浮动0(R为电导率,L0为磁导率).当时间增大时,它逐渐下沉,当时间减小时,它又逐渐/上浮0.这样便可以用一块随时间的变化而/沉浮0的/载流0导电平面近似代替回线源中的均匀大地,从而方便的求出地表任一点的异常场.这就是/浮动薄板法0的基本原理.这一方法应用较早.3.2烟圈理论解释法[35]在均匀大地上,敷设输入阶跃电流的回线,当发11684期薛国强,等:瞬变电磁法正反演问题研究进展送回线中电流突然断开时,在下半空间中就要被激励起感应涡流场以维持在断开电流以前存在的磁场,此瞬间的电流集中于地表附近,并按r-4规律衰减(r 为中心至观测点的距离).随后,面电流开始扩散到下半空间中,在切断电流后的任一晚期时间里,感应涡流呈多个层壳的0环带0形,并形成一系列与发送回线同形状并且向下及向外扩散的/电流环0,通常称之为/烟圈0.大地感应涡流在地表面产生的电磁场可近似地用圆形电流环表示.这些电流环就像由发射回线吹出的/烟圈0,其半径随着时间增大而扩大,其深度随时间延长而加深.这就提示我们:当计算均匀半空间的地面瞬变电磁响应时,可以用某一时刻的镜像电流环来代替.随着时间的延长,涡流场向下及向外扩展.依据计算的结果,涡流场极大值将沿47b 倾斜锥面扩展,计算均匀半空间的瞬变电磁响应时,可以把/烟圈0看作一系列的二次发送线圈,很容易地计算出在某时刻沿地面测线的响应值,以及在某个测点的响应值随时间变化的规律.在层状介质中,仍然保持同样的/烟圈0效应,只是/烟圈0的传播将逐渐局限于导电地层中.3.3 人机对话自动反演法.先根据地质资料及定性分析结果,给出初始模型进行正演计算,将正演计算结果用改进的阻尼最小二乘法与实测数据对比拟合.如果拟合结果不满足要求,就修改模型参数,再进行正演计算,然后再对计算结果对比,如此重复直到满意为至.为了减少多解性的影响,采用了可行方向法,控制计算参数的变化范围.以上过程全部在人的控制下由计算机自动完成.瞬变电磁测深数据拟合问题可以表示为如下的约束条件:min U (x ),x -[x i [ x ,(1)其中U x 为目标函数,x =x 1,x 2,,x mT,x --为x i 的下界, x 为x i 的上界.令:x 1i=x 0i +$x i ,t \1x 0i+t $x i ,t <1(2)其中$x i 由下式决定:A TA +A I $x =A TB式中A 为雅可比矩阵,I 为单位阵,A 为阻尼因子,其中t 由下式决定:t =x i -x 0i/$x i ,$x i >E0,-E [$x i [E x ---x 0i/$x i ,$x i <-E(3)其中E 为给定小正数.把满足上述三式的x 1i作为下次迭代新的初值,如此迭代运算,直到求出x 1,x 2,,x m 的最佳值为至.3.4 人工神经网络反演法目前传统电磁资料反演方法理论相对来说比较复杂,不易学习掌握;后期数据处理量大,计算复杂,很难进行实时反演.人工神经网络特别适合处理不确定性和非结构化信息.地质勘探和开发中的大量信息就是这种信息,同时,它也避开了具体复杂的电磁场计算,只要经过适当的学习训练就能够解决那些复杂的实际问题;而且它还具有学习记忆功能,能够一边工作一边学习,使得瞬变电磁法的反演工作具备了延续性和继承性,便于推广.人工神经网络是模拟人脑机理和功能的一种新型计算机和人工智能技术,在数据处理中避免了数据分析和建模中的困难,采用拟人化的方法进行处理,特别适合于不确定性和非结构化信息处理,因此在地质勘探中具有重要意义.人工神经网络反演法,它不要求工作人员有丰富的工作经验,它避开了具体复杂的电磁场计算,只要经过适当的学习训练就能够解决复杂的实际问题,而且还具有学习记忆功能,它一边工作一边学习,使得瞬变电磁法的反演工作具有延续性和继承性.随着专家系统的不断完善,该方法将有广阔的发展前景.3.5 成像类反演由于地震勘探研究相对比较成熟,在20世纪80年代后期,人们提出了在电磁勘探数据中采用拟地震解释法,进行成像处理.大地电磁法和瞬变电磁法都取得了一定的进展.成像类反演主要有两种,一种是时频等效转换方法,即:通过一个经验公式,把TEM 数据等效转换成平面波场数据[36~38],借用M T 数据的拟地震思路,求取反射系数序列进行成像[39].另外一种就是波场转换方法[40~44],并进一步进行拟地震偏移成像[8,45].在场源激励下,地下介质中产生涡流,在阶跃波断开后,涡流不会立即消失,而是有一个过渡过程,在这个过程中,由地下介质所产生的二次感应场经过了一个由无到有,由小到大,到极大,到衰减,再到1169地球物理学进展23卷无的过程.对于地下某一深度z,磁场微分量d B/d t 初始值为零,经过建场,到达极大值,最后衰减到零.某一频率或者某一时间的电磁场分布在地下的任何深度内.从地球物理勘探的角度看,电磁探测深度与仪器的检测灵敏度、地电情况、围岩情况、噪声电平等因素有关.在比较理想的情况下,也可能探测到埋深相当于几倍的趋肤深度(或者扩散深度)的地下地质体,在复杂地质情况下,也可能探测不到埋藏在趋肤深度(或者扩散深度)以内的地质体.但是,总的来说,对于同一介质的同一深度,扩散场的视电阻率与平面波场的视电阻率对此深度地电性结构应该有相同的反映.在一维近似的情况下,把趋肤深度和扩散深度等效认为是电磁场的探测深度一种从瞬变电磁测深数据向平面波场转换的时间t频率f对应关系.210f=t.(4)忽略位移电流以后,低频电磁场满足扩散方程,不能直接对电磁探测数据进行成像,扩散方程中的波数含有一次项,波动方程中的波数含有二次项,两个方程的形式不同.但是,电磁扩散方程与地震波动方程间存在有趣的数学对应形式,从波场到时域场的波场正变换式H m(t)=12P t3Q]0S e-S24t U(S)d S.(5)变换得到的虚拟波场,不仅满足波动方程,而且还类似于地震子波一样,具有传播、反射、透射特征.瞬变电磁场偏移成像问题与地震勘探中的弹性波偏移成像问题有相似之处,但也不完全相同.如果将扩散的瞬变电磁场变换为波场处理,将该波场从地面向地下反向外推进行偏移成像,形成瞬变电磁偏移方法,将对提高瞬变电磁场的分辨率具有重要意义.虽然变换出来的波场在形式上与地震波场一样都满足波动方程,但是由于两种波场的物理背景不同,它们之间存在一些重要区别,前者是与感应的瞬变电磁衰减曲线相对应的/反射0子波,是虚拟的,而后者是在弹性介质中传播的客观存在的地震子波.另外,虚拟波场在每一介质中的传播速度不仅与本介质的电导率有关,而且还受相邻介质电导率的影响,这与地震波场是不同的.瞬变电磁法的三维反演问题将是研究的热门.由于计算的复杂性,国内研究较少,国外学者研究相对较多,不断探索更多的反演方法[46~51].国内对时频联合反演进行了研究,并取得初步成果[52].其它不同方法的三维正反演问题也得到了研究[53~57],其研究成果会对瞬变电磁法有所借鉴.但随时着计算技术的不断发展,三维反演问题得到解决,瞬变电磁法的勘探精度会得到很大的提高.4研究展望在理论与应用中,瞬变电磁法已以取得了很大的进展,但是一些焦点问题还没有解决,一些研究还没有真正达到实用化,还需要研究者做更多的突破.在注意发展研究复杂地电条件下二、三维问题正反演的同时,更应注意实际应用效果,发展瞬变电磁测深资料与其它勘探资料的联合反演;投入力量研究瞬变电磁拟地震的偏移成像技术;在回线源非共中心点情况下,对水平分量的正演模拟计算和资料反演都研究的不够.需要加强水平分量资料的观测,及多分量数据的正反演问题研究.瞬变电磁法目前大多停留在一维反演阶段,三维反演则处于方法研究或者特定方式的应用阶段,所以,要实现三维反演的实用化,成果推广化,还需要走一段很长的道路.在三维正演计算中,大多研究场源为电偶极子或者接地长导线形式,对于应用较广的矩形回线源情况,由于其复杂性,正演计算研究相对较少.另外,对航空瞬变电磁法,海洋瞬变电磁法,井地瞬变电磁法,井下瞬变电磁法的正反演理论与方法的研究还远远不够.只有这些研究有足够多的突破,才有可以推进瞬变电磁法的全面发展,展示该方法最美好的一面.参考文献(References):[1]胡祥云,杨迪坤,刘少华,胡正旺.环境与工程地球物理的发展趋势[J].地球物理学进展,200621(2):598~604.H u X Y,Yang D K,L iu S H,H u Z W.Th e developing trendsof environmental and engineering g eophysics[J].Progres s in geoph ysics,200621(2):(598~604).[2]曾昭发,刘四新,刘少华.环境与工程地球物理的新进展[J].地球物理学进展,2004,19(3):486~491.Zeng Z F,Liu S X,Liu S H.T he new progress of en viron men-tal and en gineering geophysics[J].Progress in geop hysics, 2004,19(3):486~491.[3]赵永贵.中国工程地球物理研究的进展与未来[J].地球物理学进展,2002,17(2):305~309.Zhao Y G.T he Pres ent S itu ation and Future of Engineerin g Geophys ics in China[J].Progres s in geophysics,2002,17(2): 305~309.[4]滕吉文.21世纪地球物理学的机遇与挑战[J].地球物理学进展,2003,19(2):208~215.T eng J W.Opportunity challenge and developing frontiers:1170。
瞬变电磁_
《地球探测与信息技术基础》课程作业题目：瞬变电磁法在地球勘探上的应用姓名：周桥立班级：064101学号：20101003648授课教师：胡祥云2013 年 04 月 20日摘要瞬变电磁法是近年来电法勘探领域一种重要方法，是根据地壳中岩石或者矿体的导电性及介电性等电学性质的差异，研究电磁场的空间或时间分布规律，从而解决各种地质问题。
目前已经发展为探测油气、金属和非金属矿产的一种重要方法，并且在深部地质构造研究、工程勘察、油气、矿产、水、地热勘探等领域得到广泛应用。
目录1．概述 (4)2. 瞬变电磁法 (5)3．正演问题的研究 (8)4．瞬变电磁勘探的应用 (12)5．总结 (14)参考文献 (15)一、概述电磁场理论的应用已经遍及地学、生命科学、医学、空间科学、信息科学等几乎所有的技术科学领域,同时这些工程技术领域对电磁理论研究也不断地提出各种新的要求. 电磁法勘探是基于研究电磁波在导电介质中传播特性,从而达到研究地下地质体赋存特性的目的. 通过天然或人工场源在大地中激励的交变电磁场,研究电磁场的空间和时间分布,分析观测到的磁场信号,得到地下目标体的电性分布特征的一种地球物理方法。
瞬变电磁测深法( Transient elect romagneticmet hod ,简称TEM) 是电磁法勘探中应用较广的一种,是近年来在工程地质勘察中普遍应用的时间域电磁探测方法. 它是利用阶跃波或其它脉冲电流场源激励,在大地产生过渡过程场,断电瞬间在大地中形成涡旋交变电磁场,测量这种由地下介质产生的二次感应电磁场随时间变化的衰减特性,从测量得到的异常信号中分析出地下不均匀体的导电性能和位置,从而推断矿体、工程基础、地下水、地质灾害、工程病态等地下目标体的分布性态. 该技术具有灵敏度高、分辨率强、探测深度大、灵活多变适应性强以及轻便、快速、廉价诸多优点,近年来发展十分迅猛,应用前景十分广阔. 目前,瞬变电磁法已经成为地球物理探测领域内的重要方法之一. 已广泛应用于水利、交通、城建、环保、考古等部门. 成功地解决了大量实际问题。
电性源瞬变电磁地空逆合成孔径成像
电性源瞬变电磁地空逆合成孔径成像李貅;张莹莹;卢绪山;姚伟华【摘要】电性源地空瞬变电磁法具有工作效率高、勘探深度大、采集信号信噪比高、适用于地形地质条件复杂地区等优点.但是,到目前为止,由于尚未建立起该方法的解释系统,大大制约了该方法的发展.本文旨在建立起完整的地空电磁探测系统,丰富整个探测系统的理论.本文围绕地空瞬变电磁法全域视电阻率定义、瞬变电磁虚拟波场的克希霍夫偏移成像、逆合成孔径成像方法三个科学问题进行了系统研究.提出了用磁场强度定义全域视电阻率的迭代算法,理论模型试验结果表明计算出的视电阻率曲线首支趋于第一层电阻率,尾支趋于最后一层电阻率,实现了全空域、全时域视电阻率的计算;在先前研究的基础上,实现了适合电性源地空装置的瞬变电磁虚拟波场的克希霍夫偏移成像;采用相关迭加技术,实现孔径内多测点数据合成,将传统的单点处理方式发展成为逐点推移多次覆盖的逆合成孔径处理方法.层状模型试验表明:(1)全域视电阻率能够光滑、完整、渐变地反映出模型的电性信息变化;(2)当改变三层模型中间层电阻率时,全域视电阻率曲线随着参数的改变分异明显,对电性层的识别容易且直观;(3)由于在电阻率计算中同时考虑了接收机高度、偏移距、时间等各参数的影响,全域视电阻率可实现全空域、全时域的视电阻率计算.含水采空区的复杂模型算例表明:(1)根据不同测线的全域视电阻率结果可以看出,在靠近采空区的位置,全域视电阻率断面可以清晰地反映出采空区的空间位置,随着测线离采空区越来越远,采空区异常越来越弱直至消失;(2)波场变换和偏移成像的结果显示存在两个电性差异较大的界面,上界面指示地表,由于空气和大地之间的电性差异较大,故该界面波场信号反映强烈,遍布整个区域,下界面异常信号则主要集中在中部,向外逐渐减弱,指示采空区;(3)逆合成孔径成像结果表明地表界面在合成前后没有变化,而采空区异常合成后范围明显变小,且异常边界清晰,指示的采空区位置与模型吻合很好.本文借助于逆合成孔径雷达成像的基本思想,建立了一套电性源瞬变电磁地空逆合成孔径成像方法.基于反函数思想结合迭代算法提出的电性源地空瞬变电磁法的全域视电阻率定义方法,实现了全空域、全时域的视电阻率计算;借鉴瞬变电磁拟地震偏移成像算法,实现了瞬变场的三维成像;借鉴逆合成孔径雷达的思想,提出电性源瞬变电磁地空逆合成孔径算法,进一步提高了成像的分辨率.采空区模型算例表明相关叠加合成确实具有增强有用信号、提高信噪比、提高分辨率的诸多优点,证实了瞬变电磁地空逆合成孔径成像方法的有效性.【期刊名称】《地球物理学报》【年(卷),期】2015(058)001【总页数】12页(P277-288)【关键词】瞬变电磁;逆合成孔径成像;全域视电阻率;克希霍夫偏移成像【作者】李貅;张莹莹;卢绪山;姚伟华【作者单位】长安大学地质工程与测绘学院,西安 710054;长安大学地质工程与测绘学院,西安 710054;长安大学地质工程与测绘学院,西安 710054;长安大学地质工程与测绘学院,西安 710054【正文语种】中文【中图分类】P6311 引言随着国民经济的飞速发展，为了缓解我国矿产资源的供求矛盾，解决地形起伏剧烈的山区和湖泊、沼泽等地面人员难以到达地区的详细地质调查工作，迫切需要提出一种适合这些地区的工作方式及其配套的解释方法（杨兵，2013；葛振华等，2009）.地空瞬变电磁法利用地表铺设的电性源或回线源发射大功率瞬变电磁场，在空中采用无人机或飞艇携带探头对磁场响应进行采集.相比较航空瞬变电磁法，由于发射源被置于地表，观测信号的信噪比更高，此外由于采用无人机或飞艇在空中对电磁响应信号进行采集，大大提高了该工作方式的安全系数；而相比较地面瞬变电磁法，这种方法大大提高了工作效率，特别是在地形条件复杂的山区、湖泊、沼泽等不易开展地面瞬变电磁法工作的区域，该法更能凸显出其独特的优势，是一种极具发展潜力的电法勘探新方法.地空瞬变电磁法的产生、发展和应用，必然会给常规的地面瞬变电磁法带来革命性的变化，测点密度将大大增加，信息采集将实现全区三维覆盖，辐射场源可以是单一的，更可以是多源的，解释方法可采用三维拟地震偏移成像解释法.如果采用多源大功率发射，有望把勘探深度提高到2～3km深，这样就可以将瞬变电磁法引入到深部探测，用于三维地质填图.实际上地空电磁探测的工作方式，早在1988年Nabighian就在水平电偶极源工作方式的基础上提出了（Nabighian，1998）；Mogi等（1998，2009）研制了GREATEM系统并成功将其应用于火山结构及海侵调查；Smith等（2001）对航空、地空和地面瞬变电磁法采集的数据进行了比较；Verma等（2010）分析了均匀半空间模型接地电性源地空瞬变电磁法的磁场垂直分量响应特征；阳贵红等（2012）对时域电性源地空电磁探测的数据预处理进行了研究；李肃义等（2013）研究了电性源时域地空电磁数据小波去噪方法.到目前为止，由于三维解释方法的研究跟不上，一直没有建立起完整的地空电磁探测体系.地空瞬变电磁探测解释理论体系中所涉及的科学问题有：地空瞬变电磁法全域视电阻率定义；瞬变电磁虚拟波场的克希霍夫偏移成像；逆合成孔径成像方法.本文在已有的研究基础上，整合了TEM三维拟地震成像技术，建立了完整的地空瞬变电磁探测解释理论体系，为整个探测系统的建立奠定了理论基础.2 全域视电阻率定义在地空探测系统中，全域视电阻率定义的问题十分重要，在实际应用上，为大大提高测量效率我们希望得到距离上不分远近、时间上不分早晚、不分采集高度的全域视电阻率，解决全域视电阻率定义的问题，地空探测系统就可以在全空域、全时域进行三维数据采集，同时配以三维的数据处理和三维的解释技术，使瞬变电磁法的观测技术和解释水平全面提升到一个更高的水平.为了给出全域视电阻率定义，在均匀大地表面放置一个水平电偶极子AB（如图1所示），发射源是位于地表沿X方向放置的长度为ds的电偶极子AB，其中通有谐变电流I＝I0e－iωt，坐标原点O位于电偶极子中心，Z轴朝下，空中任意一点M，在地表的投影P距原点距离为r，r与X轴的夹角为φ.图1 均匀大地表面电偶极子示意图Fig.1 Diagram of electric dipole on the homogeneous earth surface从电磁场满足的亥姆霍兹方程出发，结合电磁场在地表处所满足的边界条件，并根据磁场和电矢量位之间的旋度关系，可得均匀大地表面水平电偶源在空中任意一点产生的频率域磁场垂直分量为其中PE＝Ids，I是电流强度，单位是A，ds是电偶极子长度，，在忽略位移电流的情况下＝－iωμσi称为介质的波数.已经证实，当采用阶跃波作为时间域电磁法的发射波形时，频率域磁场和时间域磁场存在如下对应关系：可见（1）式是关于电阻率ρ的函数，显然（2）式也是电阻率的函数，下面分析一下（2）式中磁场强度随电阻率的变化关系.图2所示为均匀半空间模型中，空中100m处，长度为1000m的电性源在不同偏移距下对应不同时刻时间域磁场垂直分量随电阻率变化的曲线形态，从图中可以看出，在偏移距为50～3000m时，电阻率在0.1～10000Ωm范围内，时间域磁场垂直分量仍然是关于电阻率参数的一致单调函数，这就是说一个电阻率值唯一地对应着一个时间域磁场垂直分量值，因此，借助于反函数存在定理的思想，可以利用时间域磁场垂直分量进行全域视电阻率定义.图2 均匀半空间模型中不同偏移距下对应不同时刻时间域磁场垂直分量随电阻率变化（a）偏移距50m；（b）偏移距500m；（c）偏移距1500m；（d）偏移距3000m.Fig.2 Vertical component of the magnetic field at different times with different offsets of homogeneous half－space（a）Offset 50m；（b）Offset 500m；（c）Offset 1500m；（d）Offset 3000m.图3 不同高度不同偏移距下改变第二层电阻率时两层模型全域视电阻率曲线（a）高度50m；偏移距500m；（b）高度50m；偏移距3000m；（c）高度100m；偏移距500m；（d）高度100m，偏移距3000m.Fig.3 Full field apparent resistivity of two－layered models with various resistivities of the second layer at both short and long offsets at different receiving heights（a）Receiving height 50m，offset 500m；（b）Receiving height 50m，offset 3000m；（c）Receiving height 100m，offset 500m；（d）Receiving height 100m，offset 3000m.给定初值ρ0，在ρ0的邻域内对Hz（t）进行泰勒展开，即在ρ0的邻域内采用“以直代曲”的方法对曲线进行近似，即保留上述泰勒展开的前两项，亦即线性项部分：经整理我们可以得到电阻率的近似表达式写成迭代格式，并辅以下角标τ以表示为视电阻率：其中反复迭代下去，直至满足其中，Hz（ρ，r，t）是已知的某一位置某一时刻对应的时间域磁场垂直分量，是电阻率为的均匀半空间在该位置该时刻产生的时间域磁场垂直分量，取ε＝0.000001～0.0001.（6）式即为电性源瞬变电磁地空探测系统的视电阻率迭代定义式.下面给出几个典型的层状模型验证上述全域视电阻率定义式的有效性.图3所示为空中50m和100m处，长度为1000m的电性源在不同偏移距下，改变第二层电阻率时得到的二层模型（模型参数：ρ1＝100Ωm，h1＝20m，ρ2＝2，5，10，30，80Ωm，200，500，800Ωm）全域视电阻率曲线，从图中可以看出，无论是在空中50m还是在100m处，以及靠近源的位置还是在远离源的位置，所得的全域视电阻率曲线不仅在早期和晚期趋于模型的第一层和第二层电阻率，还能在过渡期反映出模型的电性信息变化.图4所示为不同偏移距（500m和3000m）下，改变第二层电阻率时得到的三层模型（模型参数：H：电阻率ρ1＝100Ωm，ρ2＝5，10，30，50Ωm，ρ3＝100Ωm，层厚h1＝20m，h2＝10m；K：电阻率ρ1＝10Ωm，ρ2＝30，60，100，200Ωm，ρ3＝10Ωm，层厚h1＝10m，h2＝20m；A：电阻率ρ1＝30Ωm，ρ2＝50Ωm，60，70，80Ωm，ρ3＝100Ωm，层厚h1＝10m，h2＝20m；Q：电阻率ρ1＝100Ωm，ρ2＝20，30，50，80Ωm，ρ3＝10Ωm，层厚：h1＝20m，h2＝10m）全域视电阻率曲线.由图可见，所有的全域视电阻率测深曲线首支均趋于ρ1，尾支均趋于ρ3，曲线中部随着参数的改变分异明显，与用电磁响应定义的视电阻率曲线比较，对电性层的识别更容易、更简单.3 克希霍夫曲面延拓成像根据Lee等（1989），Hoop等（1996），陈本池等（1999）及李貅等（2001，2005，2011，2013），克希霍夫积分偏移处理的是波场，所以我们首先要将时间域磁场通过波场变换转换为虚拟波场，利用瞬变电磁场与虚拟波场间存在的数学关系式进行波场变换：式中，Hm（r，t）为时域瞬变响应扩散场，U（r，τ）为以波速传播的虚拟波场，自变量τ是时间平方根的量纲.（9）式是典型的第一类Fredholm型算子方程，已知瞬变场求波场显然具有“不适定性”.采用最优化技术和预条件正则化共轭梯度算法求解其反变换问题，可以实现从全域瞬变场到虚拟波场的变换.在获得了观测平面上瞬变电磁虚拟波场值和地下空间的连续速度分布后，依据其场满足的波动方程，可给出Kirchhoff积分解（李貅等，2010）：其中为某一时刻地下瞬变电磁虚拟波场延拓点的瞬时速度，可通过等效导电平面法和近点线性三维空间插值算法，获得三维速度数据体（陈本池等1999；李貅等2001，2010）.由于偏移是一个求取记录的逆过程，所以可设波场的上行波为，则（10）式可写为这即是波场的向下延拓.图4 第二层电阻率变化时不同偏移距下三层模型全域视电阻率曲线（a）H型，偏移距500m；（b）K型，偏移距500m；（c）A型，偏移距500m；（d）Q型，偏移距500m；（e）H型，偏移距3000m；（f）K型，偏移距3000m；（g）A型，偏移距3000m；（h）Q型，偏移距3000m.Fig.4 Full field apparent resistivity of three－layered models with different resistivities of the second layer at two typical offsets（a）H model，offset 500m；（b）K model，offset 500m；（c）A model，offset 500m；（d）Q model，offset 500m；（e）H model，offset 3000m；（f）K model，offset 3000m；（g）A model，offset 3000m；（h）Q model，offset 3000m.对上式应用三维边界元技术就可以实现对虚拟波场进行偏移，这一过程实际上就是获取记录的逆过程，这样就可以确定地下反射面上作为二次虚拟波源的空间位置，进而可以获得地下地质体的空间分布图像.4 逆合成孔径成像根据李貅等（2012）和戚志鹏等（2012），航空瞬变电磁中的合成孔径方法借助合成孔径雷达的思想，把尺寸较小的真实天线孔径用数据处理的方法合成一个较大的等效孔径的发射线圈，达到提高分辨率、增加勘探深度的目的.逆合成孔径雷达通常是指固定的雷达对运动的目标进行纵向和横向高分辨率的观测成像，形式上与合成孔径雷达正好相反，但实质相同（Berrizzi等，1996；保铮等，2001；李宁，2011）.瞬变电磁地空系统通常是发射源在地面固定不动，也可以是多个源同时发射，接收探头在空中移动观测，显然与逆合成孔径雷达的观测装置相似，完全可以借鉴逆合成孔径雷达的思想.通过对空中观测点周边一定范围内的观测信号进行相关叠加处理，相当于合成了一个较大孔径的接收线圈，进而提高对地下探测的纵向和横向分辨率.图5是地空瞬变电磁逆合成孔径成像装置的示意图，A、B是位于地表的电性源供电电极，以此向地下发送瞬变电磁场，到达地质体后会感生出二次涡流，断电后，在空中用磁探头可以接收到由二次涡流产生的二次磁场随时间衰减的电磁响应信号.由于同一地质体产生的电磁响应信号在L孔径范围内是相关的，将L孔径的中点作为参考点i，此点的波场值可表示为U（ri，τ），其中ri为i点到L内（编号为：－N，…，N）某点的距离，τ为相对移动量.把参考点与L孔径内各点相关计算得到的最大相关系数作为权系数，将该权系数分别乘上各点的波场值，叠加到中心点，最终可得到中心点的合成值为图5 逆合成孔径示意图Fig.5 Diagram of inverse synthetic aperture imaging其中为参考点与孔径L内某点两列波场的最大相关系数，（k＝－N，…，－1，0，1，…，N）为最大相关系数ρmax（ri，τm）对应的移动量值.在一条测线上依次移动，则可得到第i＋1，i＋2，i＋3，… 点为中心点的合成值.5 模型计算为了验证本文方法的有效性，设计了如图6的充水采空区模型，采空区赋存于一层煤层之中，煤层上下分别为细砂岩和砂质页岩，底层为石英砂岩，详细的地电参数与地质体尺寸列于表1.图6 充水采空区复杂模型示意图Fig.6 Diagram of a complex model with one mined out water body表1 充水采空区复杂模型参数Table 1 Parameters of a complex model with one mined out water body地层尺寸（长×宽×高）电性参数／Ωm1610m×1600m×100m 100煤层410m×410m×10m 500充水采空区100m×100m×10m 10砂质页岩1610m×1600m×40m 150石英砂岩细砂岩1610m×1600m×800m 800图7给出了该复杂模型的俯视图，结合三维示意图6可以明确定位煤层与采空区的位置.模型设计采用301×301×100的10m立方体均匀网格，采用电性源长度AB＝100m，发射电偶源中点设在模型正中，在地面上空100m处取三条剖面LineX81、LineX86、LineX91（如图7所示），偏移中心点的距离分别为0m、50m、100m，对选定的三条剖面绘制全域视电阻率断面等值线图（图8），由图8可见，LineX81测线位于采空区中部，低值视电阻率等值线明显地反映出了采空区的空间位置，LineX86测线位于采空区边上，对充水采空区的低阻异常仍有显示，LineX91测线已远离采空区，采空区的异常基本消失.这说明用本文方法定义的全域视电阻率图像准确地反映了采空区的空间分布位置.用经三维正演计算获得的地面上空100m处的磁场强度数据进行波场变换，其结果如图9所示，存在两个界面，上层界面是地表面的显示，由于空气和大地间的电性差较大，故波场信号反映强烈，且遍布整个区域；下层界面高强度异常信号集中在中部，向外围逐渐减弱，显然是采空区的显示.图7 采空区复杂模型俯视图与发射源位置及测线布置Fig.7 Top view of the complex model with one mined out water bodies and position of the transmitting source and the survey line arrangement图8 选定测线视电阻率断面图（a）LineX81测线；（b）LineX86测线；（c）LineX91测线.Fig.8 Apparent resistivity contours of the selected lines图9 三维波场变换成像结果（a）三维立体图；（b）切片图.Fig.9 Results of 3－D wave field transform（a）Three dimensional map；（b）Slice map.图10 三维偏移成像结果（a）三维立体图；（b）切片图.Fig.10 Results of 3－D migration imaging（a）Three dimensional map；（b）Slice map.对于逆合成孔径的计算可以采用先合成后偏移的方法，即对波场变换后的数据进行逆合成孔径计算后，再进行偏移成像；也可以采用先偏移后合成的方法，即先对波场变换后的数据进行偏移成像（如图10所示），然后再进行逆合成孔径计算.图11为先偏移后合成的成像结果，由图11可见，地表界面合成前后没有变化，但采空区异常合成后范围明显变小，而且异常边界清楚，可见逆合成孔径效果明显.6 结论（1）本文基于反函数思想结合迭代算法提出的电性源地空瞬变电磁法的全域视电阻率定义方法，实现了在全空域、全时域的视电阻率计算，模型计算结果表明：视电阻率曲线对地层的电性特征分异明显.图11 逆合成孔径成像结果（a）三维立体图；（b）切片图；（c）去掉高能反射结果图.Fig.11 Results of inverse synthetic aperture imaging（a）Three dimensional map；（b）Slice map；（c）Map without high energy reflection.（2）理论模型试验表明，用磁场强度定义的视电阻率曲线没有假异常，首支渐近线趋近于第一层电阻率，尾支渐近线趋近于最后一层电阻率.（3）瞬变电磁拟地震偏移成像算法的提出，实现了瞬变场的三维成像，使得瞬变电磁法的解释技术提高到一个新水平.（4）为进一步提高成像分辨率，借鉴逆合成孔径雷达的思想，提出了电性源瞬变电磁地空逆合成孔径算法.经采空区数值模型实验验证，相关叠加合成确实具有增强有用信号、提高信噪比、提高分辨率的诸多优点，证明了瞬变电磁地空逆合成孔径成像方法的有效性.ReferencesYang B.2013.Analysis and discussion on dependency of overseas market and supply of China non－ferrous metals resources.Mineral Exploration（in Chinese），4（1）：8－11.Ge Z H，Zheng Y J，Xiao R G.2009.Features and suggestions of import ＆export trade of mineral products in China.Resources＆Industries（in Chinese），11（3）：103－109.Yang G H.2012.Data processing research on electrical－source of time domain Ground－airborne electromagnetic （in Chinese）［M aster′s thesis］.Changchun：Jilin University.Nabighian N M.1988.Electromagnetic Methods in Applied Geophysics－Theory（Volume 1）.Tulsa OK Society of Exploration.Mogi T，Tanaka Y，Kusunoki K，et al.1998.Development of grounded electrical－source airborne transient EM （GREATEM ）.Exploration Geophysics，29：61－64，doi：10.1071／EG998061.Mogi T，Kusunoki K，Kaieda H，et al.2009.Grounded electricalsource airborne transient electromagnetic（GREATEM）survey of Mount Bandai，north－eastern Japan.Exploration Geophysics，40：1－7，doi：10.1071／EG08115.Smith R S，Annan A P，Mcgowan P D.2001.A comparison of data from airborne，semi－airborne，and ground electromagneticsystems.Geophysics，66（5）：1379－1385，doi：10.1190／1.1487084. Verma S K，Mogi T，Allah S A.2010.Response characteristics of GREATEM system considering a half－space model.20thIAGA WG 1.2on Electromagnetic Induction in the Earth，Giza Egypt，September 18－24. Li S Y，Lin J，Yang G H，et al.2013.Ground－Airborne electromagnetic signals de－noising using a combined wavelet transformalgorithm.Chinese J.Geophys.（in Chinese），56（9）：3145－3152，doi：10.6038／cjg20130927.Qiang J K，Luo Y Z，Tang J T，et al.2010.The algorithm of alltime apparent resistivity for Airborne Transient Electromagnetic（ATEM）survey.Process in Geophysics （in Chinese），25（5）：1657－1661，doi：10.3969／j.issn.1004－2903.2010.05.018.Spies B R，Eggers D E.1986.The use and misuse of apparent resistivity in electromagnetic methods.Geophysics.，7（51）：1462－1471，doi：10.1190／1.1442194，doi：10.1190／1.1442753.Piao H R.1990.Theory of Electromagnetic Sounding（in Chinese）.Beijing：Geological Publishing House.Wu Q.2012.Forward modeling of large loop transient electromagnetic field and wave－field transform theory′s study（in Chinese）［Master′s thesis］.Xi′an：Chang′an University.Lee K H，Liu G，Morrison H F.1989.A new approach to modeling the electromagnetic response of conductive media.Geophysics.，54（9）：1180－1192.Li X，Xue G Q，Song J P，et al.2005.An optimize method for transient electromagnetic field－wave field conversion.Chinese J.Geophys.，（in Chinese），48（5）：1185－1190.Qi Z P，Li X，Wu Q，et al.2013.A new algorithm for full－timedomain wave－field transformation based on transient electromagneticmethod.Chinese J.Geophys.，（in Chinese），56（10）：3581－3595，doi：10.6038／cjg20131033.Chen B C，Li J M，Zhou F T.1999.Wave－field conversion method for transient electromagnetic field.Oil Geophysical Prospecting （in Chinese），34（5）：539－545.Xue G Q，Li X，Qi Z P，et al.2011.Study of sharpening the TEM pseudo－seismic wave－form.Chinese J.Geophys.，（in Chinese），54（5）：1384－1390，doi：10.3969／j.issn.0001－5733.2011.05.02.Li X，Guo W B，Hu J P.2001.The method and application effects of pseudo －seismic interpretation of TEM.Journal of Xi′an Engineering University （in Chinese），23（3）：42－45.Chen B C，Li J M，Zhou F T.1999.Quasi wave equation migration of transient electromagnetic field.Oil Geophysical Prospecting （in Chinese），34（5）：546－554.Li X，Qi Z P，Xue G Q，et al.2010.Three dimensional curved surface continuation image based on TEM pseudo wave－field.Chinese J.Geophys.（in Chinese），53（12）：3005－3011，doi：10.3969／j.issn.0001－5733.2010.12.025.Qi Z P.2012.3Dcontinuation imaging technology with a synthetic aperture based on the transient electromagnetic method［Ph.D.thesis］（in Chinese）.Xi′an：Chang′an University.Li X，Xue G Q，Liu Y A，et al.2012.A research on TEM imaging method based on synthetic－aperture technology.Chinese J.Geophys.，（in Chinese），55（1），333－340，doi：10.6038／j.issn.0001－5733.2012.01.034.de Hoop A T.1996.Transient electromagnetic vs.seismic prospecting—acorrespondence principle.Geophysical Prospecting，44（6）：987－995，doi：10.1111／j.1365－2478.1996.tb00187.x.Berizzi F，CorsiniG.1996.Autofocusing of inverse synthetic aperture radar images using contrast optimization.IEEE Transactions On Aerospace And Electronic Systems.32（3）：1185－1191，doi：10.1109／7.532282.Li N.2011.Study on inverse synthetic aperture radar imaging methods ［Master′s thesis］（in Chinese）.Nanjing：Nanjing University of Aeronautics and Astronautics.Bao Z，Zhang S H，Xing M D.2001.Inverse synthetic aperture radar imaging of maneuvering target.Acta Aeronautica et Astronautica Sinica （in Chinese），22（Supplement）：S6－S15.附中文参考文献杨兵.2013.中国有色金属矿产对外依存度与资源可供性之辨析.矿产勘查，1：8－11.葛振华，郑勇军，肖荣阁.2009.中国矿产品进出口贸易特点及建议.资源与产业，11（6）：103－109.阳贵红.2012.时域电性源地＿空电磁探测数据预处理研究［硕士论文］.长春：吉林大学.李肃义，林君，阳贵红等.2013.电性源时域地空电磁数据小波去噪方法研究.地球物理学报，56（9）：3145－3152，doi：10.6038／cjg20130927.强建科，罗延钟，汤井田等.2010.航空瞬变电磁法的全时域视电阻率计算方法.地球物理学进展，25（5）：1657－1661，doi：10.3969／j.issn.1004－2903.2010.05.018.朴化荣.1990.电磁测深法原理.北京：地质出版社.吴琼.2012.大回线源电磁场正演与波场变换理论研究［硕士论文］.西安：长安大学.李貅，薛国强，宋建平等.2005.从瞬变电磁场到波场的优化算法.地球物理学报，48（5）：1185－1190.戚志鹏，李貅，吴琼等.2013.从瞬变电磁扩散场到拟地震波场的全时域反变换算法.地球物理学报，56（10）：3581－3595，doi：10.6038／cjg20131033.陈本池，李金铭，周凤桐，1999.瞬变电磁场的波场变换算法.石油地球物理勘探，34（5）：539－545.薛国强，李貅，戚志鹏等.2011.瞬变电磁拟地震子波宽度压缩研究.地球物理学报，54（5）：1384－1390，doi：10.3969／j.issn.0001－5733.2011.05.028.李貅，郭文波，胡建平.2001.瞬变电磁测深快速拟地震解释方法及应用效果.西安工程学院学报，23（3）：42－45.陈本池，李金铭.1999.瞬变电磁场拟波动方程偏移成像.石油地球物理勘探，34（5）：546－554.李貅，戚志鹏，薛国强等.2010.瞬变电磁虚拟波场的三维曲面延拓成像.地球物理学报，53（12）：3005－3011，doi：10.3969／j.issn.0001－5733.2010.12.025.戚志鹏.2012.瞬变电磁三维合成孔径延拓成像方法研究［博士论文］.西安：长安大学.李貅，薛国强，刘银爱等.2012.瞬变电磁合成孔径成像方法研究.地球物理学报，55（1）：333－340，doi：10.6038／j.issn.0001－5733.2012.01.034.李宁.2011.逆合成孔径雷达成像方法研究［博士论文］.南京：南京航空航天大学. 保铮，张守宏，邢孟道.2001.机动目标的逆合成孔径雷达成像.航空学报，22（增刊）：S6－S15.。
瞬变电磁原理与应用PPT课件
0.2622 0.3171 0.382
0E+0
0.4594
440 460 480 500 520 540 560 580 600 620 640 660 680 700 720 740 760 780 800 0.5542
1E+3 8E+2
0.6691 0.8064
4E+2
0.9712
0E+0
440 460 480 500 520 540 560 580 600 620 640 660 680 700 720 740 760 780 800
2. 对于近地表浅层探测时， “烟圈”理论计算的视深度不在 适用，在浅层探测时计算结果严重偏离实际的深度，计算 的深度从距地表20米至40米之间，探测明显存在着盲区；
3. 对浅层1号坑的1米深的低阻异常反映不出来，但对于2号 坑的2米深的低阻异常，反映清晰，但异常体的深度位置 在视深度－25米至40米位置，与实际的目标体埋深不符。
4.0E+4
0.0974
3.0E+4 2.0E+4 1.0E+4
0.1198 0.1473 0.1798
0.2172
440 460 480 500 520 540 560 580 600 620 640 660 680 700 720 740 760 780 800 8E+3
6E+3 4E+3 2E+3
-20
22 20
-30
19 18
17
-40
16 15
14
-50
13 12
-60
11 10
8
-70
7 6
5
瞬变电磁法理论与应用研究进展_薛国强
第22卷 第4期地 球 物 理 学 进 展V ol .22 N o .42007年8月(页码:1195～1200)P ROG RESS IN G EOP H YSICSA ug . 2007瞬变电磁法理论与应用研究进展薛国强1, 李 貅2, 底青云1(1.中国科学院地质与地球物理研究所,北京100029;2.长安大学地质工程与测绘工程学院,西安710054)摘 要 对国内外瞬变电磁法的方法研究概况、理论研究、仪器状况以及在应用领域的研究进展情况做出了综述性评价.瞬变电磁法的理论研究主要涉及正演方法、反演方法、资料处理方法.瞬变电磁法的应用领域包括地面、海洋、航空以及地下等不同工作场地.同时指出瞬变电磁法的发展趋势为研究三维正反演、多分量观测、成像技术、微信号分析等,以及注重发展海洋、井下、航空瞬变电磁法等.关键词 瞬变电磁法,理论与应用,研究进展中图分类号 P631 文献标识码 A 文章编号 1004-2903(2007)04-1195-06The progress of TEM in theory and applicationXUE Guo -Qiang 1, LI Xiu 2, DI Qing -Yun1(1.I nstitute o f Geo log y and Geoph ysics ,Chinese A cademy o f Sciences ,Beijing 100029,Ch ina ;2.S chool o f Geo log y and S ur vey En gineer ing ,Chang ,an Un iver sity ,X i ,an 710054,China )A bstract We give the to tal condition of T EM study include techno lo gy ,theo ry ,instrument and application ,the the -o ry study mainly include fo rw ard calculation metho d ,inver se problem and da ta pro ce ssing scheme .T he applicatio n o f it include ea rth sur face ,sea and aerial survey ,subsurface explora tion .We also give advance prog ress the o f T EM ,include study ing 3D fo rwa rd and inv erse calcula tion ,multi -co mpo nent survey techno lo gy ,imaging study ,and w eak signal a nalyses ,as well as ,sea ,tunnel ,aerial T EM develo ping .Keywords tr ansient elect romag netic me tho d ,theo ry and applicatio n ,study pro g ress收稿日期 2007-04-10; 修回日期 2007-06-20.基金项目 中国博士后基金(2005038388),中国科学院王宽诚博士后奖励基金,中国科学院知识创新项目((KZCZ -yw -113).)国家自然科学基金重点项目(50539080)联合资助.作者简介 薛国强,男,1966年生,1989年于西安地质学院获学士学位,2002年于长安大学获硕士学位,于2005年于西安交通大学获博士学位.研究方向为电磁探测理论与应用.(E -mail :ppxueguoqiang @ )0 引 言瞬变电磁法(TEM )是一种时间域的电磁探测方法.介质在一次电流脉冲场激励下会产生涡流,在脉冲间断期间涡流不会立即消失,在其周围空间形成随时间衰减的二次磁场.二次磁场随时间衰减的规律主要取决于异常体的导电性、体积规模和埋深,以及发射电流的形态和频率[1～6].因此,我们可以通过接收线圈测量的二次场空间分布形态,了解异常体的空间分布[1～6].瞬变电磁法在上世纪30年代最早由前苏联科学家提出,当时采用的是远区工作模式.但是,利用电流脉冲激发供电偶极形成时域电磁场的是1933年由美国科学家L .W .Blau 最先提出,当时利用不同电导率地层界面电磁波的反射与地震反射波信号的相似性,进行了大量的实验和比较.到了上世纪50～60年代,前苏联科学家成功地完成了瞬变电磁法的一维正、反演,建立了瞬变电磁法的解释理论和野外工作方法之后,瞬变电磁法才开始进入实用阶段.上世纪60年代以后,当意识到时间域电磁测深法可以利用远远小于期望探测深度的收发距时,该方法有了一个快速发展.随之,“短偏移”、“晚期”、“近区”、等技术研究迅速发展起来.美国等西方国家在上世纪70～80年代之间,短偏移法一直处于研究和试验阶段,未被广泛运用,而长偏移法已得到了应用,特别是在地热调查和地壳结构的调查中.随后一些专家对瞬变电磁法的一维正反演及方法技术进行了大量研究[7].我国于70年代初开始研究TEM ,长地 球 物 理 学 进 展22卷春地质学院、中南工业大学、西安地质学院、中国有色金属工业总公司矿产地质研究院和西安物化探研究院等单位先后研制了电磁系统,并进行了理论和方法技术研究,主要采用近区方式的中心回线法和重迭回线法进行工作,取得了一批有价值的研究成果及大量成功的实际应用实例[8].这一时期代表作有朴化荣著《电磁测深法原理》,牛之链著《脉冲瞬变电磁法原理》及方文藻著《瞬变电磁测深法原理》.在最近几年的瞬变电磁法应用实践中,还出现了蒋邦远编著的《实用近区磁源瞬变电磁法勘探》,李貅编著的《瞬变电磁测深的理论与应用》.近几年,瞬变电磁法的理论与应用研究越来越活跃.特别是在工程勘察领域的应用.1 理论研究进展1.1 正演计算对瞬变电磁法一维情况下的正演计算大多采用先在频率域进行讨论,得到层状介质下的电磁响应表达式,然后把讨论结果转换到时间域来.从频率域到时间域响应的转换,可以采用的方法有:GAV ER-S TEH FES T逆拉氏变换方法,延迟谱方法,线性数字滤波方法,折线化正余弦变换法[1,2].二维数值计算多采用有限差分法进行.用两个无限长直导线近似作为发送回线源,可用均匀半空间的解析解在t>0时刻将源转化为初始条件加入.对所研究的空间区域作差分离散,利用准静态近似处理空中边界,在地-空界面向上延拓一个网格,采用五点差分显式格式[9],在每一时间步计算网格空间各点的场量,然后进行时间的逐步递推,就能直接模拟电磁波的传播及其与地质体的相互作用过程.有限元法对频率域计算比较有利.因为它可以很方便地处理急剧变化的和倾斜的电导率分界面和地形等问题.由于频域电磁勘探中的二维定解问题相对简单,且用有限元法求解十分方便,因此,先在频域中用有限元法求解,然后变换到时间域,是瞬变电磁二维正演的一种有效途径.2.5维人工源时间域电磁场数值模拟问题是目前国际上尚未妥善解决的计算地球物理疑难问题之一.我国从上世纪80年代开始着手研究2.5维电磁场的数值模拟,采用有限元法对时间域和频率域电源瞬变电磁场的2.5维响应进行正演数值模拟.在频率域已经实现三维源二维地质结构的数值计算[10,11],可以首先在频率域计算电磁场响应值,然后再把计算结果变到时间域[12,13,14].三维瞬变电磁场正演计算只能用数值方法,主要方法有三类:有限差分法[15,16]、有限元法[14]和积分方程法[17].前两种方法要求对所计算的全部区域进行离散化,所占用的计算机容量较大.第三者只要对异常区域进行离散化,存贮量大大减小.因为积分方程法只需要计算小体积异常区的场,不必计算整个区域的场.这一点使积分方程算法比微分方程法具有更大的优越性.但是积分方程法求解要遇到某些更困难的数学问题,仅适合模拟简单模型.在开放域的地球物理问题中,微分方程算法要解决边界条件和源的处理问题.瞬变电磁三维直接时域有限差分法数模拟已经有报道[15,16].而且还进一步地通过设定等效位移电流,使原适合波动场的时域有限差分算法(FDTD)能够适用于扩散场问题.直接时域算法的引入,展现了电磁场在地下随时间传播的全过程,使T EM场的时间特性被直接反映出来,从而给复杂的物理过程描绘出清晰的物理图象.目前国内还没有应用较好的3D正演软件. 1.2 反演问题研究浮动薄板解释法是一种根据视纵向电导曲线的特征值直观地划分地层的近似解释方法.因此称为“视纵向电导解释法”,也有人把该法形象地称为“浮动薄板解释法”.烟圈理论解释法在均匀大地上,当发送回线中电流突然断开时,在下半空间中就要被激励起感应涡流场以维持在断开电流以前存在的磁场,大地感应涡流在地表面产生的电磁场可近似地用圆形电流环表示.这些电流环就像由发射回线吹出的“烟圈”,其半径随着时间增大而扩大,其深度随时间延长而加深.这就提示我们:当计算均匀半空间的地面瞬变电磁响应时,可以用某一时刻的镜像电流环来代替.人机对话自动反演法先根据地质资料及定性分析结果,给出初始模型进行正演计算,将正演计算结果用改进的阻尼最小二乘法与实测数据对比拟合.如果拟合结果不满足要求,就修改模型参数,再进行正演计算,然后再对计算结果对比,如此重复直到满意为至.人工神经网络反演是模拟人脑机理和功能的一种新型计算机和人工智能技术,在数据处理中采用拟人化的方法进行处理,特别适合于不确定性和非结构化信息处理.它不要求工作人员有丰富的工作经验,它避开了具体复杂的电磁场计算,只要经过适11964期薛国强,等:瞬变电磁法理论与应用研究进展当的学习训练就能够解决复杂的实际问题,而且还具有学习记忆功能,它一边工作一边学习,使得瞬变电磁法的反演工作具有延续性和继承性.随着专家系统的不断完善,该方法将有广阔的发展前景.联合时-频分析反演是借助于信号分析领域内的新成果发展而成的技术,它同时分析TEM信号的时间和频率域响应曲线特性,从而进行定性和定量分析解释[18].TEM接收到的电压衰减曲线形态比较相近,一般情况下计算出全程视电阻率以及其它电性参数,再进行分析解释,这种常规方法对TEM信号的利用率不高.联合时间频率分析可以表示涡流场在不同时刻,不同频率分量在地下激发和衰减过程.设观测的时间曲线和由此进行快速付氏变换所得的频率曲线为已知,给定初始反演参数,采用广义逆矩阵迭代超定方程,结果表明,用两种曲线进行联合反演所得的参数精度比常规解释精度高.成像类反演主要有两种,一种是时频等效转换方法,即:通过一个经验公式,把TEM数据等效转换成平面波场数据,借用M T数据的拟地震思路,求取反射系数序列进行成像[19].另外一种就是波场转换方法,并进一步进行拟地震偏移成像.如果将扩散的瞬变电磁场变换为波场处理,将该波场从地面向地下反向外推进行偏移成像,形成瞬变电磁偏移方法[20],将对提高瞬变电磁场的分辨率具有重要意义.瞬变电磁法的三维反演问题将是研究的热门.由于计算的复杂性,国内研究较少,国外学者研究相对较多.但随时着计算技术的不断发展,三维反演问题得到解决,瞬变电磁法的勘探精度会得到很大的提高.1.3 资料处理方法瞬变电磁法的资料工作,主要是根据瞬变电磁响应的时间特性和剖面曲线特性,以及工区的地质地球物理特征,通过分析研究,划分出背景场及异常场.在有条件的情况下,应该进行数据全区转换[21],小波去噪等处理,很有必要利用一维、二维、三维正演数值计算得到工作地区的先验知识,便于资料分析解释,尽可能地利用钻井资料,区域地质资料.目前的瞬变电磁探测方法对地下目标体的评价精度低,对于所绘制的电性参数图件,一般采用二次电压衰减曲线和由此算得的视电阻率值及视纵向电导参数来进行解释.有时候也利用多测道剖面图,时间常数,二次感应电压与背景场感应电压的比值等参数等进来定性分析解释,但多数情况下仍处于定性和“看图识字”的阶段.对瞬变电磁法测深资料定量解释还局限于单点一维反演,且反演效果不佳,很多情况下是靠解释人员的工作经验及地质先验知识来对测深结果做出判断,人为性较大.其分辨率远不能满足实际工程地质的需要,这一状况也限制了方法的推广应用.2 应用状况按照工作空间来分,TEM可以分为:地面TEM.地下TEM、航空TEM、海洋TEM.在我国,前者应用较多,后三者的应用较少.2.1 地面勘察瞬变电磁测深法以其直接测二次场;可以进行近区观测;受静态位移的影响小;测地工作简单,工效高等特点广泛应用于金属矿勘探、油气田勘探和煤田地质勘探、工程地质勘察中.瞬变电磁法的应用还有:海洋调查;水文地质、环境工程灾害地质调查;高速公路、铁路的路线选择;水库、水坝、机场、桥梁以及高层建筑的选址、地下洞体的工程勘察;地下水、地热资源的探查;对长江、黄河等大江大河堤坝的质量探查;对古代地下建筑群位置和规模的无损探测;隧道掌子面超前预报、煤矿超前探测.国外的地球物理工作者把瞬变电磁测深法用于冰冻层、海水入浸等地质方面调查.以前,瞬变电磁法只局限于金属矿勘探,1992年以后随着仪器的智能化与数字化,瞬变电磁法开始步入工程、环境、灾害地质调查中,如探测地下采空区、陷落柱等煤田灾害,划分地下断层、寻找地下水,金属矿产勘探、石油、煤炭等非金属矿产调查、工程场地地质勘察、隧道超前地质预报等领域.取得了良好的效果.地面主要应用的装置形式有磁源重迭回线法、大回线源法、中心回线法、电偶源、LOTEM 应用较少.解决了大量工程地质问题.取得了很多实际应用和理论研究成果.2.2 地下瞬变电磁法2.2.1 矿井瞬变电磁法由于井下施工环境的限制,不可能采用大线圈装置形式,只能采用边长小于3m的多匝小线框,工作效率高.井下瞬变电磁法在施工过程中,具有一定的方向特性.在探测过程中,根据地质任务的不同,可调整发射线框与煤层底板之间的角度.可以将线圈直立于巷道内,当线圈面平行巷道掘进前方,可1197地 球 物 理 学 进 展22卷进行超前探测,或探测煤层内的含水构造异常体;也可以将线圈置于巷道底板,并设置一定的夹角,通过探测了解煤层底板下一定深度范围内含水异常体垂向和横向的发育规律.发射线框和接收线框分别为匝数不等且完全分离的2个独立线框,以便与地下(前方)异常体产生最佳耦合响应.在工作面回采前,利用现有巷道对相邻工作面未掘巷道实施侧向探测,确定隐伏含水异常范围及位置,为掘进中防、探水工作提供技术依据.减少、取代和部分取代现行超前探测工作[22,23].2.2.2 隧道掌子面瞬变电磁法随着我国铁路、公路隧道及南水北调引水隧洞建设的加快,隧道(洞)的勘查设计时间比较短,在隧道(洞)工程建设开发之前,很难提供足够的时间来进行地面地球物理勘探和岩土工程地质勘察,况且由于受到复杂地质和水文地质条件的影响,不良地质情况预报无论在理论方法和技术手段上都没有得到有效解决,难以应对隧道施工中可能出现的水体灾害问题.长安大学,西安煤炭研究所率先用瞬变电磁法进行隧道掌子面前方不良地质体超前预报.取得了较好的效果.类半平面问题在掌子面高和宽大于5倍的回线尺寸的情况下,近似认为是半平面.引进了微分电导成像技术,解决掌子面前方含水构造探测问题.采用3米长的回线框作为发送,直立地贴近于掌子面,用特制的接收磁探头(固定于回线正中心)接收来自掌子面前方的二次感应电压信号,通过一系列各种隧道效应校正,得到较好的资料.通过近年来几十个隧道探测,取得了较多的经验和效果. 2.3 海洋瞬变电磁法目前用于海洋油气勘探的可控源电磁方法主要有频率域和时间域电磁法两类.海洋电磁勘探时一般将接收装置布设于测量区域的海底,发射源为电偶极源,用船拖着在距观测区域一定的距离范围(2～10km)内移动,实现多偏移距测量.频率域方法与时间域方法的差别仅在于发射波形和电磁响应函数不同,而在解释参数上时间域电磁法可以获得地层的充电率参数,为研究地层的含油气性提供比电阻率更为直接的信息.国内在这一方面的理论研究[24,25]多于实际应用.3 仪器状况过去勘探队伍和设备主要集中于地质队,科研院所等,仪器主要为国产仪器为主.当时的国产商品化的仪器主要产自西安、廊坊、长沙等地.现在一些非物探专业部门,独立法人的企业,还有个人购买了进口仪器和国产仪器.受市场化作用的影响,国产仪器受冲击力较大.市场上大的项目还是依靠进口仪器来完成.进口的国外的仪器主要有加拿大GEON-ICS公司生产的PROT EM系统(PROTEM-37,47, 57,67);PH OENIX公司生产的V6,V8系统;美国ZONGE公司生产的GDP-32系统等.澳大利亚SIROT EM仪器和前苏联的M PPO仪器.前面三款仪器在国内应用较多,后两者仪器应用较少.国外仪器最大特点是接收机全智能化,形成多功能的工作站.目前,国际上商业用的固定翼时间域航空电磁系统主要有澳大利亚产TEM PEST、加拿大产GEOTEM R、GEOTEM DEEPTM、M EGA TEM、澳大利亚产QU ES TEM R450、南非产SPEC TREM 等系统.国内应用很少.在我国,70年代初期开始着手研究瞬变电磁仪器系统,最早研制并投入生产的是地矿部物化探研究所,其生产的的仪器型号为WDC-1、WDC-2,后又研制了IGGE TEM-20瞬变电磁系统;1988年西安物化探研究所采用脉冲压缩技术研制成功大功率的LC-1系统并投入生产,后来又研制EM RS-1、EM RS-2瞬变电磁仪.1992年长沙高新技术产业区智通新技术研究所与中南工业大学合作生产出SD-1型、SD-2型仪器,白云仪器厂在此基础上研制了M SD-1,BYF5MSD1瞬变电磁系统.1996年石油天然总公司与西安石油仪器厂开始研制用于深部探测的大功率高精度瞬变电磁仪器.北京矿产地质研究所王庆乙教授研制了TEM S-3S瞬变电磁仪器;吉林大学林君教授研制了ATEM-2瞬变电磁仪器[26,27],重庆奔腾数控技术研究所研制了WTEM 系统;限于我国电子技术工艺水平原因,虽然一些厂家在生产用于深部勘探的大功率仪器和用于浅部探测的小功率仪器,但在生产工艺上、原器件焊接技术上和性能稳定上都与国外仪器有很大的差距.目前国内没有生产航空电磁法仪器的厂家.接收磁探头方面:GDP32配有接收磁芯磁探头.V6,V8系列仪器配有专用的空心接收结圈,不太方便.国内已经有专门生产此探头的单位,可以生产浅、中、深三种不同型号的磁探头.受市场和国外仪器的影响,国内仪器需要向高性能,高分辨方向发展,乃是当务之急.将地震勘探仪器及方法中已采用的先进技术有针对性性地引入,发展多通道电磁法、阵列式观测、实现多次空间迭加和多次覆盖.将电子工业中业中新型技术有针11984期薛国强,等:瞬变电磁法理论与应用研究进展对性地引入到仪器设备中如:卫星定位、同步、大功率管控制技术、微波通讯技术等.设备器件耐高温高寒技术;仪器的系列化、多功能化;发展成熟的先进的数据处理技术、可视化软件、高温超导磁探头[28],提高磁场探测的灵敏度.4 研究展望(1)理论研究方面:研究高分辨率的数据采集系统;研究三维精细解释方法,发展对地下目标体的3D反演;以可视化技术为主体的计算机技术的应用;研究三分量技术;研究多孔径瞬变电磁探测理论.发展成像技术及瞬变电磁测深资料与其它勘探资料的联合反演.瞬变电磁法的激电效应特征研究.(2)研究弱信号技术:包括瞬变电磁早期信号的研究和晚期信号研究.瞬变电磁发射线圈和接收线圈之间存在互感,以及它们自身都有自感,其等效回路相当于RLC电路,自感信号与电磁响应信号迭加在一起,如何对这种信号进行分离,如何确定有效观测起始时间,从而确定瞬变电磁最小探测深度.当然瞬变电磁的探测深度与仪器的分辨能力、地下介质导电情况、围岩特性、目标体与围岩差异等因素有关.晚期信号能量较弱,最后被背景噪声淹没,研究有效的信噪分离方法,以增大瞬变电磁法的探测精度和深度.(3)方法技术方面:接收装置阵列化.瞬变电磁法拟地震的时间域和空间域多次覆盖技术研究是十分重要的.另外,多分量的观测技术研究也不应忽视.研究多源发射聚焦探测方法.(4)仪器方面:采集系统图像化.发展多通道大功率、多功能、智能化电测系统,高灵敏度高性能探头的研制,是现代仪器研制开发的难点.仪器向数字化、多功能化、智能化、自动化、轻便化方向发展.特别是研制多分量瞬变电磁法系列(航空、地面、井下).考虑到岩层各向异性对测量结果的影响,以及大地电磁干扰场是一种水平极化场,因此通常选择测量磁场的垂直分量.观测水平分量可以研究断面的各向异性,这一方面的理论研究有待发展.(5)应用方面:除了通常应用于金属矿产及石油资源的勘查外,还应用在地下水、地热环境及工程勘查中.另外,井下瞬变电磁法的研究和应用还处于起步阶段,瞬变电磁法的隧道超前地质预报及深部构造精细探测方面应该加强研究.还要大力研究海洋瞬变电磁法和深部剩余油瞬变电磁法精细探测.井—地瞬变电磁勘探系统,已经取得初步成果,希望在油田,矿产资源领域得到较好应用.石油天然气勘查.深部油藏监测,剩余油勘察等仍然是瞬变电磁法的一个难题.所以,要实现三维反演的实用化,成果推广化,还需要走一段很长的道路.5 结 论我国对瞬变电磁法的理论研究与生产实践已有30年的历史,但目前在理论研究、仪器研制方面尚处于初级阶段,而仪器大多是在国外仪器基础上的开发和改进,加之应用推广并不普遍,因此,在我国瞬变电磁法研究和普及还要有很长的路走.总的来说,瞬变电磁法取得了很大进展,在勘察领域得到相对广泛的认可,但是,还要研究非共中心点情况下的瞬变电磁正反演;还需要加强水平分量资料的观测,及多分量数据的正反演问题研究;复杂介质情况下,比如连续介质、地形起伏介质;瞬变电磁正演计算还研究的不够深入.拟地震成像解释的研究工作,虚拟波场速度分析,拟地震定量反演解释的研究才刚刚开始.特别是二维地电断面情况下的拟地震解释还处于探索阶段.另外,对于拟地震解释的精度评价研究做的还不够深入.在三维正演计算中,大多研究场源为电偶极子或者接地长导线形式,对于应用较广的矩形回线源情况,由于其复杂性,正演计算研究相对较少.航空瞬变电磁法,国外有开展,我国开展较少.感 谢 课题组王妙月、石昆法研究员,王光杰副研究员,王若博士对一些问题的探讨与帮助.中国地质大学王家映教授,江苏大学闫述研究员对一些疑难问题的研究和指导.参 考 文 献(References):[1] 朴化荣.电磁测深法原理[M].北京:地质出版社,1990,[2] 方文藻,李予国,李貅.瞬变电磁测深法原理[M].西安:西北工业大学出版社,1993[3] 牛之琏.时间域电磁法原理[M].中南工业大学出版社,1993,[4] 何继善.电法勘探的发展和展望[J].地球物理学报,1997,40(增刊):308～316.[5] 蒋邦远.实用近区磁源瞬变电磁法勘探[M].地质出版社,1998.[6] 李貅.瞬变电磁测深的理论与应用[M].西安:陕西科学技术出版社,200291-4,102～105.[7] 薛国强.回线源瞬变电磁法对地成像方法研究[D].博士学位论文,西安:西安交通大学,2005.[8] 李貅.瞬变电磁虚拟波场的三维曲面延拓成像研究[D].博士学位论文,西安:西安交通大学生,20051199。
第三章瞬变电磁法
图3.1.1 瞬变电磁法原理示意图
117
由于瞬变电磁法测量的是导体内涡流的过渡过程，观测是在脉冲间歇期间进行 的，不存在一次场源的干扰，观测参数为纯二次场，是电磁法中唯一可采用同点装置 的方法，探测目标耦合最紧，获得的响应最强。磁性源激发，不受接地条件限制。 在金属矿勘查中， 主要用于寻找良导性的致密块状、 团块状、 网脉状硫化物矿体。 但对于浸染状硫化物矿体的探测效果不佳。 瞬变电磁法具有以下特点： （1）由于 TEM 法接收的是纯二次场，因而不受一次场的影响； （2）可以采用高密度时序采样，纵向分辨率较高； （3）穿透低阻覆盖能力强，勘探深度大； （4）发射用不接地回线，不受地表接地条件限制； （5）一般矿山主要干扰是电场，相对 TEM 干扰较小。
119
甚至重要的作用。应用此图时，回线长可看着任意比例尺，并以它来归一测点距、深 度和异常体大小。
图 3.1.4 通过主剖面的垂直剖面下方的一次场磁力线分布图
2、正常场 正常场是剖面测量中的一个术语，它是指局部异常响应的背景。典型的正常场就 是均匀非磁性导电半空间表面的瞬变响应。层状大地也可视为相当某一电阻率的半空 间。 （1）重叠回线 设半径为 a 的单匝圆回线铺设于均匀非磁性导电半空间表面上，在 t=0 的瞬间， 回线中的电流 I 阶跃地下降为零： I t ≤ 0 I (t ) = 0 t > 0 则均匀半空间的电动势响应为： (3.1.4)
§3.2 野外工作方法
3.2.1 常用装置类型及功能 常用装置类型及功能
瞬变电磁法用于找矿勘查能够较准确地确定地质体的倾向、埋深、走向等。野外 工作装置形式繁多，并是电磁法中唯一能进行同点发射—接收的方法。根据勘查任务 的不同可非常灵活地选用装置，常用的装置组合有以下几种（图 3.2.1） 。
瞬变电磁法拟地震成像研究进展
第２２卷 第３期地 球 物 理 学 进 展Ｖｏｌ．２２ Ｎｏ．３２００７年６月（页码：８１１～８１６）ＰＲＯＧＲＥＳＳ ＩＮ ＧＥＯＰＨＹＳＩＣＳＪｕｎｅ ２００７瞬变电磁法拟地震成像研究进展李 貅１， 薛国强２， 郭文波３（１．长安大学地质工程与测绘工程学院，西安７１００５４； ２．中国科学院地质与地球物理研究所，北京１０００２９；３．西安交通大学电子与信息工程学院，西安７１００４９）摘 要 瞬变电磁法是工程地球物理勘察中应用较多的一种勘探方法之一，瞬变电磁法拟地震成像方法研究是当前电磁探测理论与应用研究的热门．文中详细叙述了瞬变电磁法拟地震成像的两种不同思路，一种是建立时间—频率转换经验公式，把ＴＥＭ数据转换成平面波数据，借鉴大地电磁法拟地震解释思路进行成像，另一种是通过积分变换，把扩散场变成虚拟波场，把电磁衰减信号变成拟地震子波，然后按照地震解释法进行解释，并分析了各自存在的的问题．指出下一步研究的方向．关键词 瞬变电磁法，拟地震，成像，波场转换中图分类号 Ｐ６３１ 文献标识码 Ａ 文章编号 １００４ ２９０３（２００７）０３ ０８１１ ０６犚犲狊犲犪狉犮犺狆狉狅犵狉犲狊狊犻狀犜犈犕狆狊犲狌犱狅 狊犲犻狊犿犻犮犻犿犪犵犻狀犵ＬＩＸｉｕ１， ＸＵＥＧｕｏ ｑｉａｎｇ２， ＧＵＯＷｅｎ ｂｏ３（１．犛犮犺狅狅犾狅犳犌犲狅犾狅犵狔犪狀犱犛狌狉狏犲狔犈狀犵犻狀犲犲狉犻狀犵，犆犺犪狀犵，犃狀犝狀犻狏犲狉狊犻狋狔，犡犻，犪狀７１００５４，犆犺犻狀犪；２．犐狀狊狋犻狋狌狋犲狅犳犌犲狅犾狅犵狔犪狀犱犌犲狅狆犺狔狊犻犮狊，犆犺犻狀犲狊犲犃犮犪犱犲犿狔狅犳犛犮犻犲狀犮犲狊，犅犲犻犼犻狀犵１０００２９，犆犺犻狀犪；３．犛犮犺狅狅犾狅犳犈犾犲犮狋狉狅狀犻犮犪狀犱犐狀犳狅狉犿犪狋犻狅狀犈狀犵犻狀犲犲狉犻狀犵，犡犻，犪狀犑犻犪狅狋狅狀犵犝狀犻狏犲狉狊犻狋狔，犡犻，犪狀７１００４９，犆犺犻狀犪）犃犫狊狋狉犪犮狋 Ｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｍｅｔｈｏｄ（ＴＥＭ）ｉｓｏｎｅｏｆｐｏｐｕｌａｒｕｓｅｄｍｅｔｈｏｄｉｎｅｎｇｉｎｅｅｒｉｎｇａｎｄｇｅｏｐｈｙｓｉｃａｌｅｘｐｌｏｒａｔｉｏｎ，ｔｈｅｓｔｕｄｙｏｆＴＥＭｐｓｅｕｄｏ ｓｅｉｓｍｉｃｉｍａｇｉｎｇｉｓｔｈｅｐｏｐｕｌａｒｓｕｂｊｅｃｔｉｎｔｈｅｔｈｅｏｒｙａｎｄａｐｐｌｉｃａｔｉｏｎｏｆｅｌｅｃｔｒｏ ｍａｇｎｅｔｉｃｄｅｔｅｃｔｉｎｇ．Ｉｎｔｈｉｓｐａｐｅｒ，ｗｅｉｎｔｒｏｄｕｃｅｔｗｏｗａｙｓｏｆＴＥＭｐｓｅｕｄｏ ｓｅｉｓｍｉｃｉｍａｇｉｎｇ，ｏｎｅｗａｙｉｓｔｒａｎｓｆｏｒｍｉｎｇＴＥＭｄａｔａｔｏｐｌａｎｅｗａｖｅｄａｔａ，ａｎｄｆｉｎｉｓｈｉｎｇＴＥＭｐｓｅｕｄｏ ｓｅｉｓｍｉｃｉｍａｇｉｎｇｂａｓｉｎｇｏｎＭＴｐｓｅｕｄｏ ｓｅｉｓｍｉｃｍｅｔｈｏｄ．Ａｎ ｏｔｈｅｒｗａｙｉｓｔｒａｎｓｆｏｒｍｉｎｇｄｉｆｆｕｓｉｏｎｆｉｅｌｄｉｎｔｏｐｓｅｕｄｏ ｗａｖｅｆｉｅｌｄ，ａｎｄｔｒａｎｓｆｏｒｍｉｎｇｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｄｅｃａｙｓｉｇｎａｌｉｎｔｏｐｓｅｕｄｏ ｓｅｉｓｍｉｃｗａｖｅｌｅｔ，ａｎｄｆｉｎｉｓｈＴＥＭｐｓｅｕｄｏ ｓｅｉｓｍｉｃｉｍａｇｉｎｇｂａｓｉｎｇｏｎｓｅｉｓｍｉｃｍｅｔｈｏｄ．Ｗｅａｎａｌｙｓｅｓｔｈｅｅｘｉｓｔｉｎｇｐｒｏｂｌｅｍｉｎｔｈｅｔｗｏｉｍａｇｉｎｇｗａｙａｎｄｇｕｉｄｅｔｈｅｄｉｒｅｃｔｉｏｎｏｆｎｅｘｔｓｔｕｄｙ．犓犲狔狑狅狉犱狊 ｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｍｅｔｈｏｄ，ｐｓｅｕｄｏ ｓｅｉｓｍｉｃ，ｉｍａｇｉｎｇ，ｗａｖｅ ｆｉｅｌｄｃｏｎｖｅｒｓｉｏｎ收稿日期 ２００６ ０６ ０７ 修回日期 ２００７ ０３ ２０．基金项目 国家自然科学基金重点项目（５０５３９０８０），中国科学院知识创新工程重要方向资助项目（ＫＺＣ２ ｙｗ １１３），中国博士后基金（２００５０３８３８８）和中国科学院王宽诚博士后奖励基金联合资助．作者简介 李貅，男，１９５８年生，博士，长安大学教授．研究方向为电磁探测理论与应用．（Ｅ ｍａｉｌ：ｌｉｘｉｕ＠ｃｈｄ．ｅｄｕ．ｃｎ）０ 瞬变电磁方法及特点瞬变电磁法的测量原理是利用不接地回线（或电偶源）向地下发送一次脉冲磁场（或电场），即在发射回线上供一个电流脉冲方波，方波后沿下降的瞬间，将产生一个向地下传播的一次瞬变磁场，在该磁场的激励下，在地质体内产生涡流，其大小取决于该地质体的导电能力，导电能力强则感应涡流强．在一次场消失后，涡流不能立即消失，它将有一个过渡过程（衰减过程），该过渡过程又产生一个衰减的二次场向地下传播．在地表用接收线圈接收二次磁场，该二次磁场的变化，将反映地下介质的电性情况，在接收机中按不同的延迟时间测量二次感应电动势，得到二次场随时间衰减的信号［１，２］．由于瞬变电磁测深法是在一次场断电后测量纯二次场，不存在一次场的干扰，此外，从富里叶变换可知，一个阶跃形脉冲实际上是由各种高频和低频谐波叠加而成的，产生的场是一种宽频带电磁波场，因此与频率域电磁法相比，瞬变电磁测深法具有以下优点［３，４，５］：地 球 物 理 学 进 展２２卷（１）断电后观测二次场，可以进行近区观测，减少旁侧影响，增强电性分辨能力；（２）可用加大功率的方法增强二次场信号，提高信噪比，从而增加勘探深度；（３）穿透高阻层能力强；（４）由于采用人工源方法，随机干扰影响小；（５）采用重叠回线装置工作，可以减小地形影响；（６）线圈形状、方位要求相对不严格，测地工作简单，工效高；（７）由于测量磁场，受静态位移的影响小；（８）通过多次脉冲频号场的重复观测叠加和空间域多次覆盖技术的应用，可以提高信噪比和观测精度；（９）可以通过选择不同的时窗窗口进行观测，有效地压制各种噪声，可以获得不同勘探深度，使剖面与测深工作与一体．１ 瞬变电磁法的发展方向在我国瞬变电磁法的研究、开发与实践已经有三十年的历史，在方法和仪器等方面取得了很大的进步［６～１２］，并且在实际中得到较广泛的应用［１３］．但从目前看，在理论研究、仪器研制方面处于初级阶段［１４］．一方面理论研究虽已解决了一维正反演问题，但在二三维的研究还未达到应用程度［１１］．国内仪器生产的规模性和持续性较差，加之应用推广并不普遍，受市场经济的影响，国内越来越多的生产单位和科研院所都购置了国外生产的多功能、多方法仪器，使国内自成体系的方法研究和发展受到制约．因此，可以认为瞬变电磁法的发展方向可以概括为以下几个方面［１５，１６］．１．１ 理论研究方面近年来，瞬变电磁法的研究和应用比较活跃［１７～２１］．在注意发展研究复杂地电条件下二三维问题正反演的同时［２２］，更应注意实际应用效果，对瞬变电磁拟地震的偏移成像技术，瞬变电磁法各种干扰信号的分离技术研究和解释方法是一个值得注意的研究方向．当前国内、外学者寻找一个数学上的处理办法，将瞬变电磁场变换成波场，即将电磁响应中与传播有关的特征提取出来，而把电磁波传播过程中与频散、衰减有关的特征进行压制或去除．如果能够实现的话就能将瞬变电磁场的求解问题转化为大家所熟悉的波动方程的求解问题，就能将地震偏移成像技术、Ｂｏｒｎ近似反演技术、脉冲谱反演技术、层析成像反演技术等用于瞬变电磁场的反演解释中［２３］．基于电磁波在导电介质中的传播和地震波在弹性介质中传播的运动相似性基础上建立起来的反演解释方法这些年来有了很大发展，取得了卓有成效的研究成果，在一定程度上缩小了解的非唯一问题，扩大了电磁法的应用范围，并可希望从观测资料中提取常规反演方法得不到的更多、更准确的地下信息．１．２ 方法技术方面瞬变电磁法拟地震的多次覆盖技术研究是十分重要的，另外，多分量、阵列式观测技术研究也不应忽视．１．３ 仪器方面主要是发展多通道，大功率、多功能、智能化电测系统，高灵敏度，高性能探头的研制，是现代电磁仪研制开发的难点．受市场和国外仪器的冲击，国内仪器需要向高性能，高分辨方向发展，乃是当务之急．２ 瞬变电磁法拟地震成像的焦点和前沿课题２．１ 基于时 频等效转换的瞬变电磁成像技术近年来，利用瞬变电磁场勘探石油、地热源和各种矿产资源的理论和应用研究工作在不断发展，对于探测埋在地下的低阻异常体，电磁法已证明是一种有效的方法［２４，２５］．但是，由于在分层、有耗媒质中电磁现象的复杂性，目前对实测数据的解释水平仍很低．近年来人们正深入研究二维、三维的复杂模型，力图更准确地描述大地中的瞬变电磁现象，并设法从电磁响应中获得地下结构的局部形状和尺寸等高分辨信息［４］．在远区情况下的瞬变电磁法拟地震成像方法相对容易［２６］．回线源瞬变电磁法是一种近区观测的电磁探测方法．由于场源的特殊性，在介质中传播的电磁场是扩散场．在研究回线源瞬变电磁法对地成像时，这一问题是不容回避的问题，这样，由扩散场向平面波场数据转换就成了问题研究的关键［２７］．文献［２７］从大量模型计算入手，通过对两种场源测深正演数据的分析、对比，以及对两种场在地下介质中传播的特性分析，建立一种从瞬变电磁测深数据向平面波场转换的时间—频率对应关系．通过大量的理论模型正演计算、曲线对比、误差分析，结合场的特性分析，从两种场的穿透深度及反映地下电性结构一致性角度，经过详细推导，得出瞬变电磁测深视电阻率数据可以转换成平面波场视电２１８３期李 貅，等：瞬变电磁法拟地震成像研究进展阻率数据的结论，构造了由时间到频率的转换关系式：２１０犳＝狋．（１）这一关系可以从两种场在地下介质中传播的时间、速度、深度等方面加以说明．电磁场在地下介质中的扩散速度与地下介质的电阻率及时间有关系．在阶跃波断开后某一时间延迟，扩散深度与速度的关系如下［２７］：犇（狋）＝∫狋０犞犱（ρ，狋′）ｄ狋′，（２）式中，犇（狋）是扩散深度，犞犱（ρ，狋）是扩散速度，ρ是介质的电阻率，狋′是断电后的时间延迟．为了推导瞬变电磁场在地下的传播速度及深度，我们分析瞬变电磁场在地下的建场机制及过程．在场源激励下，地下介质中产生涡流，在阶跃波断开后，涡流不会立即消失，而是有一个过渡过程，在这个过程中，由地下介质所产生的二次感应场经过了一个由无到有，由小到大，到极大，到衰减，再到无的过程．对于地下某一深度ｚ，磁场微分量ｄＢ／ｄ狋初始值为零，经过建场，到达极大值，最后衰减到零．某一频率或者某一时间的电磁场分布在地下的任何深度内．从地球物理勘探的角度看，电磁探测深度与仪器的检测灵敏度、地电情况、围岩情况、噪声电平等因素有关．在比较理想的情况下，也可能探测到埋深相当于几倍的趋肤深度（或者扩散深度）的地下地质体，在复杂地质情况下，也可能探测不到埋藏在趋肤深度（或者扩散深度）以内的地质体．但是，总的来说，对于同一介质的同一深度，扩散场的视电阻率与平面波场的视电阻率对此深度地电性结构应该有相同的反映．在一维近似的情况下，把趋肤深度和扩散深度等效认为是电磁场的探测深度．这样，我们令两个深度相等，从而可以找到频率与时间的类似公式（１）的关系式．文献［１６］提出了瞬变电磁法对地成像理论和成像数值计算方法：根据所构造的扩散场向平面波场转换时的时间频率关系式，把时间转换成频率，这样就完成了扩散场数据向平面波场的转换．由转换后的视电阻率值求出平面波场的波阻抗，再由计算出来的波阻抗，建立关于电性层反射系数序列的方程组，再通过线性规划法，求出反射系数序列，最终以反射系数为参数绘制成像剖面．给出了瞬变电磁法对地成像数值计算步骤：（１）给出各种不同地电模型，分别进行正演计算．计算出视电阻率值，根据不同的地表电阻率，对早期道数据进行校正．（２）通过转换关系，把时间延迟变成频率，对理论模型正演数据或者实测视电阻率值进行域的变换，把时间域扩散场视电阻率值变成平面波场视电阻率值．（３）由平面波场视电阻率值在频率域求出波阻抗．（４）以波阻抗为参数，构建方程组．（５）用线性规划法求出反射系数序列．（６）最终以反射系数为参数进行成像．２．２ 基于波场转换的瞬变电磁成像技术由于瞬变电磁场满足的微分方程事实上是一个扩散方程，因而不能采用目前大家熟悉的波动方程求解方法．所谓瞬变电磁场的波场变换是指：通过数学积分变换，将满足扩散方程的时域瞬变电磁场转换为满足波动方程的波场，然后借助于地震中发展起来的一些比较成熟的成像方法技术，求解被探目标体的物性和几何参数［２８］［２９］［３０］．Ｗｅｉｄｅｌｔ［３１］于１９７２年，Ｋｕｎｅｔｚ［３２］于１９７２年，以及Ｌｅｖｙ［３３］等于１９８８年，Ｌｅｅ［３４］等于１９８９年的一些研究成果，都揭示了在层状大地介质中，电磁扩散方程与地震波动方程间存在有趣的数学对应形式，但他们研究问题的着眼点都是将对应地电模型的波场模拟结果变换成时域电磁响应．但是更能激起学者们研究兴趣的与此恰相反，即采用波场逆变换，将已知时域场转换为波场，这将有利于偏移以及更加复杂的成像技术的应用．从波场到时域场的波场正变换式［３４］犎犿（狋）＝１∫∞０τｅ－τ２４狋犝（τ）ｄ狋，（３）这一变换过程称为正问题（ｄｉｒｅｃｔｐｒｏｂｌｅｍｓ）．如果反过来，已知时域场求波场，则称为反问题（ｉｎ ｖｅｒｓｅｐｒｏｂｌｅｍｓ）．众所周知，反问题往往和不适定性（ｉｌｌ ｐｏｓｅｄｎ ｅｓｓ）紧密相联，各种各样的反问题，不仅出现于地球物理问题中，而且出现于数学本身．波场正变换式是典型的第一类Ｆｒｅｄｈｏｌｍ型算子方程，而它的反问题有一个重要的特征，就是所谓的“不适定性”．关于“适定”（ｗｅｌｌ ｐｏｓｅｄｎｅｓｓ）和“不适定”（ｉｌｌ ｐｏｓｅｄｎｅｓｓ）的概念是Ｈａｄａｍａｒｄ为了描述数学物理问题与定解条件的合理搭配，于２０世纪初引入的［３５］．瞬变电磁场所满足的扩散方程主要刻画电磁涡流场的感应扩散特征，而基于扩散方程的偏移成像３１８地 球 物 理 学 进 展２２卷及反演方法，一般对电性界面的分辨能力较差，如何寻找到一个数学上的处理办法，将瞬变电磁场变换成“波场”，即提取出电磁响应中与传播有关的特征，压制或去除电磁波传播过程中与频散、衰减有关的特征．将瞬变电磁场的求解问题转化为波动方程的求解问题，就能将地震偏移成像技术；Ｂｏｒｎ近似反演技术；层析成像反演技术等用于瞬变电磁场的反演解释中［３６］．文献［３７］提出了从瞬变电磁场到波场的优化算法．在波场的正变换中，在保证瞬变电磁场的计算精度的同时，应用两步最优化算法，成功地控制和减少了积分系数个数和离散数字积分的采样点个数，解决了在波场反变换式中，由于积分系数过多而产生的欠定方程组的问题，同时改善了第一类算子方程的不适定性．在波场反变换中，将正则化算法用于逆变换过程的计算中，通过采用偏差原理和Ｎｅｗｔｏｎ迭代格式选出最优的正则化参数，使得反变换所得到的波场稳定、可靠．通过对数值计算结果与已知波场函数对比，证明了该方法的有效性和实用性．大量的理论计算表明：这种变换得到的虚拟波场，不仅满足波动方程，而且还类似于地震子波一样，具有传播、反射、透射特征．通过理论研究已经解决了瞬变电磁场的波场变换问题，这就为瞬变电磁场的偏移成像创造了条件．但瞬变电磁场偏移成像问题与地震勘探中的弹性波偏移成像问题有相似之处，但也不完全相同，需要指出：ＴＥＭ虚拟波场是经过数学变换得到的，不是客观存在，因此，还有其特殊性．２．３ 逆时偏移成像方法２．３．１ ｓｔｒａｔｔｏｎ ｃｈｕ积分偏移成像方法Ｚｈｄａｎｏｖ等［３８～４０］人，借鉴了地震勘探中的逆时偏移概念，经过十余年来的系统研究，对时间域的瞬变电磁场进行了逆时偏移成像的深入研究，提出了偏移电磁场的概念，并且建立在逆时偏移电磁场基础上对二，三维反演问题也展开了最新的研究．瞬变电磁场的偏移应用的是ｓｔｒａｔｔｏｎ ｃｈｕ积分：－１４π∫狋′０Γ｛（狀·犈）Δ犌 ＋（狀×犈）×Δ犌 ＋（狀×犎）σ犌 ｝·ｄ狊ｄ狋＝犎（狉′，狋），（４）其中犌＝（μ０σ）１／２２π１／２（狋′－狋）３／２ｅｘｐ－μ０σ狘狉′－狉狘４狘狋′－狋狘狌（狋′－狋）．（５）这里狌（狋′－狋）＝１ 狋＜狋′０ 狋＞｛狋′，（６）Ｇ 是引入的格林函数，犈和犎是地表Γ上的电磁场，狀是Γ的外法线方向的单位矢量．（４）式表明，可以通过将地表观测到的电磁场值向下延拓，得到地下电磁场的分布，延拓后的电磁场值在界面处发生突变，由此可以确定出地下介质的地电构造．Ｚｈ ｄａｎｏｖ等人主要对孤立和多个电性异常体的模型做了实验研究，成像效果较好．２．３．２ ｋｉｒｃｈｈｏｆｆ积分偏移成像方法将瞬变电磁场利用２．２中所述的算法进行波场变换得到“波场”，将波场分析的原理即地震波场从地面向地下反向外推进行偏移成像的方法用于对瞬变电磁场的解释，形成瞬变电磁场的偏移方法．也就是对观测电磁数据向下延拓，构建产生响应的地下电阻率图像．文献［１５］中给出了克希霍夫积分法进行电磁波场偏移成像处理的数值计算方法．波场在地下传播可以用波动方程描述，其方程为：Δ２狌－１狏２ ２狌狋２＝犉上述方程的克希霍夫积分解为：狌（狓，狔，狕，狋）＝－１４π犙［狌］ 狀１（）狉－１狉 狌 ［］｛狀－１狏狉 狉 狀 狌 ［］｝狋ｄ犙＋犉狉０，（７）对于狌（狓，狔，狕，狋），当把狋改变－狋时，既是狑（狓，狔，狕，狋）＝狌（狓，狔，狕，－狋）时，狑（狓，狔，狕，狋）仍可满足同样的波动方程．对于狑是时间向前的问题，对于狌就是时间“倒退”的问题．可以把反射界面的各点看为同时激发上行波的源点．这样，我们可以把地面上的接收点作为二次发射源，将这些信息值时间“倒退”到原来状态，寻找反射界面的波场函数，以确定反射界面．设自激自收的上行波为犌（狓，狔，狕０，狋），它是地下反射界面作为源点发射的波场犵（狓，狔，狕，狋）在地面上狕＝狕０上的值，由（７）式可得：犵（狓，狔，狕，狋）＝－１４π 犙０狀１（）狉－１狉 狌 ｛狀 －１狏狉 狉 狀 ｝狋犌ξ，η，ζ０，狋＋狉（）狏ｄ犙＋犉狉０，（８）上式中：犌ξ，η，ζ０，狋＋狉（）狏取狋＋狉狏，是由于现在考虑的是波动的逆过程；事实上是由犵（狓，狔，狕，狋）在地表４１８３期李 貅，等：瞬变电磁法拟地震成像研究进展面上的值犌ξ，η，ζ０，（）狋求出地下的波场值，从而确定反射界面，这即是波场的向下延拓．将三维边界元技术引入到克希霍夫积分计算中来，其最大的好处是能够实现三维曲面延拓计算，这在地震勘探中是很难实现的，这是由于公式（８）中存在 犌狀项，地震勘探中无法进行地表垂向梯度测量，相反瞬变电磁法利用梯度探头可以很容易地实现地表垂向梯度测量．另外，在速度分析中，可将瞬变电磁法中的反演技术融入其中，一方面，实现连续速度分析；另一方面，也可将反演与偏移成像相结合，提高解释的可靠程度．通过对时间域的瞬变电磁场进行了逆时偏移成像的深入研究，将瞬变电磁虚拟波场从地面向地下反向外推进行偏移成像，建立了有瞬变电磁法特点的偏移成像方法．３ 结 论提高探测精度是研究瞬变电磁法对地成像的目的．文中分析了瞬变电磁法的发展方向及拟地震解释研究的焦点问题．详细叙述了瞬变电磁法拟地震成像的几种不同思路．当前研究中存在的问题是：拟地震解释中的速度分析研究不完善，两种转换方法的成像结果是一种近似解释方法，对复杂界面的成像有待于进一步研究，在实用性方面及应用推广方面还需要做更多的工作．低频电磁数据经数值计算转换成虚拟波动场数据时，虚拟波场波形随着偏移距和地层电导率的增大，存在严重的波形展宽效应．使得计算得到的虚拟波场数据分辨能力不足，对多个界面的多个反射分辨不理想．所以瞬变电磁法拟地震解释技术还要做更多的研究．参 考 文 献（Ｒｅｆｅｒｅｎｃｅｓ）：［１］ 李金铭．电法勘探新进展［Ｍ］．北京：地质出版社，１９９５．［２］ 李貅．瞬变电磁测深的理论与应用［Ｍ］．西安：陕西科学技术出版社，２００２．［３］ 何继善．电法勘探的发展和展望［Ｊ］．地球物理学报．１９９７，４０（增刊）：３０８～３１６．［４］ 王家映．大地电磁拟地震解释［Ｍ］．石油工业出版社．１９９５．［５］ 方文藻，李予国，李貅．瞬变电磁测深法原理［Ｍ］．西安：西北工业大学出版社，１９９３．，［６］ 嵇艳鞠，林君，王忠．瞬变电磁接收装置对浅层探测的畸变分析与数值剔除［Ｊ］．地球物理学进展，２００７，２２（１）：２６２～２６７．［７］ 周逢道，林君，刘长胜，等．浅海底瞬变电磁法接收天线频率特性研究———带宽对浅层分辨的影响［Ｊ］．地球物理学进展，２００６，２１（４）：１３４２～１３４５．［８］ 嵇艳鞠，林君，朱凯光，等．利用瞬变电磁技术进行地下水资源勘察［Ｊ］．地球物理学进展，２００５，２０（３）：８２８～８３３．［９］ 嵇艳鞠，林君，于生宝，王忠，王静．ＡＴＴＥＭ系统中电流关断期间瞬变电磁场响应求解的研究［Ｊ］．地球物理学报，２００６，４９（６）：１８８４～１８９０．［１０］ 刘长胜，林君．海底表面磁源瞬变响应建模及海水影响分析［Ｊ］．地球物理学报，２００６，４９（６）：１８９１～１８９８．［１１］ 闫述．瞬变电磁场的直接时域数值分析［Ｊ］．地球物理学报，２００２，４５（２）：２７５～２８４．［１２］ 闫述，陈明生．瞬变电磁场资料的联合时 频分析解释［Ｊ］．地球物理学报，２００５，４８（１）：２０３～２０８．［１３］ 王华军，罗延钟．中心回线瞬变电磁法２．５维有限单元算法［Ｊ］．地球物理学报，２００３，４６（６）：８５５～８６２．［１４］ ＬｉｎＪ，ＣｈｅｎｇＤＦ，ＬｉｕＧＸ，犲狋犪犾．Ｄｅｖｅｌｏｐｍｅｎｔｏｆｔｈｅａｒｒａｙｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｓｙｓｔｅｍａｎｄｉｔｓｅｘｐｅｒｉｍｅｎｔｓ［Ｊ］．Ｊｏｕｒｎａｌｏｆｇｅｏｓｃｉｅｎｔｉｆｉｃｒｅｓｅａｒｃｈｉｎｎｏｒｔｈｅａｓｔａｓｉｎ．２０００，３（１）：９９～１０６．［１５］ 李貅．瞬变电磁虚拟波场的三维曲面延拓成像研究［Ｄ］．博士学位论文，西安交通大学，２００５．［１６］ 薛国强．回线源瞬变电磁法对地成像方法研究［Ｄ］．博士学位论文，西安交通大学，２００５．［１７］ 罗延钟，等．时间域航空电磁法一维正演研究［Ｊ］．地球物理学报，２００３，４６（５）：７１９～７２４．［１８］ 熊彬，等．瞬变电磁２．５维反演中灵敏度矩阵计算方法［Ｊ］．地球物理学进展，２００４，１９（３）：６１６～６２０．［１９］ 吴钦宝，等．中国东部煤矿深部开采中的地质勘探技术［Ｊ］．地球物理学进展，２００５，２０（２）：３７０～３７３．［２０］ 白登海，等．时间域瞬变电磁法中心方式全程视电阻率的数值计算［Ｊ］．地球物理学报，２００３，４６（５）：６９７～７０４．［２１］ 王赤军，等．７７ＫＳＱＵＩＤ磁强计在ＴＥＭ测量上的应用研究［Ｊ］．地球物理学报，１９９９，４２（９）：１４３～１４８．［２２］ 熊彬，罗延钟．电导率分块均匀的瞬变电磁２．５维有限元数值模拟［Ｊ］．地球物理学报，２００６，４９（２）：５９０～５９７［２３］ 于鹏，王家林，吴健生．大地电磁场成像方法综述与新进展［Ｊ］．地球物理学进展．２００３，１８（１）：５３～５８．［２４］ 李茂荣，等．西部地下水资源勘查的综合地球物理问题—以鄂尔多斯盆地白垩系地下水勘查为例［Ｊ］．地球物理学进展，２００１，１６（３）：１１９～１２６．［２５］ 张保祥，刘春华．瞬变电磁法在地下水勘查中的应用综述［Ｊ］．地球物理学进展．２００４，１９（３）５３７～５４２．［２６］ 郭文波，李貅，薛国强．瞬变电磁快速成像解释系统研究［Ｊ］．地球物理学报．２００５，４８（６）：１４００～１４０５．［２７］ 薛国强，李貅，宋建平．从瞬变电磁测深数据向平面波场数据的等效转换［Ｊ］．地球物理学报，２００６，４９（５），１５３９～１５４５．［２８］ 陈本池，李金铭，周凤桐．瞬变电磁场拟波动方程偏移成像［Ｊ］．石油地球物理勘探，１９９９，３４（５）：５４６～５５４．［２９］ 陈本池，周凤桐，李金铭．瞬变电磁场的波场变换研究［Ｊ］．物探与化探．１９９９，２３（３）：１９５～２０１．［３０］ ＴａｅＪｌ，ＪｕｎｇＨＳ，ＨｅｅＪＫ，ＹｏｏｎｈＳ，ＫｉＨａｌｉｅ．Ｅｌｅｃｔｒｏｍａｇ ｎｅｔｉｃｒａｖｅｒｔｉｎｅＴｏｍｏｇｒａｐｈｙＵｓｉｎｇＰｒｏｂａｂｉｌｉｔｙＷａｖｅ ｆｉｅｌｄ５１８地 球 物 理 学 进 展２２卷Ｔｒａｎｓｆｏｒｍ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ．２００２，６７（１）：６７～６９．［３１］ ＷｅｉｄｅｌｔＰ．Ｔｈｅｉｎｖｅｒｓｅｐｒｏｂｌｅｍｏｆｇｅｏｍａｇｎｅｔｉｃｉｎｄｕｃｔｉｏｎ：Ｚｅｉｔ．ｆｕｒ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ．，１９７２，３８，２５７～２９８．［３２］ ＫｕｎｅｔｚＧ．Ｐｒｏｃｅｓｓｉｎｇａｎｄｉｎｔｅｒｐｒｅｔａｔｉｏｎｏｆｍａｇｎｅｔｏｔｅｌｌｕｒｉｃｓｏｕｎｄｉｎｇｓ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，１９７２，３７，１００５～１０２１．［３３］ ＬｅｖｙＳ，ＯｌｄｅｎｂｕｒｇＤ，ＷａｎｇＪ．Ｓｕｂｓｕｒｆａｃｅｉｍａｇｉｎｇｕｓｉｎｇｍａｇ ｎｅｔｏｔｅｌｌｕｒｉｃｄａｔａ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，１９８８，５３：１０４～１１７．［３４］ ＬｅｅＫＨ．ＬｉｕＧ，ＭｏｒｒｉｓｏｎＨＦ．Ａｎｅｗａｐｐｒｏａｃｈｔｏｍｏｄｅｌ ｉｎｇｔｈｅｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｒｅｓｐｏｎｓｅｏｆｃｏｎｄｕｃｔｉｖｅｍｅｄｉａ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，１９８９，５４（６）：１１８０～１１９２．［３５］ 肖庭延，于慎根，王彦飞，反问题的数值解法［Ｍ］．北京：科学出版社，２００３．１ ３５，１４３～１５３．［３６］ ＬｅｅＳ，ＭｅｍｅｃｈａｎＧＡ．Ｐｈａｓｅ ｆｉｅｌｄｉｍａｇｉｎｇ：Ｔｈｅｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｅｑｕｉｖａｌｅｎｔｏｆｓｅｉｓｍｉｃｍｉｇｒａｔｉｏｎ［Ｊ］．Ｇｅｏｐｈｙｓｉｃｓ，１９８７，５２（５）：６７８～６９３．［３７］李貅，薛国强，宋建平，郭文波．从瞬变电磁场到波场的优化算法［Ｊ］．地球物理学报，２００５，４８（５）：１１８５～１１９０．［３８］ ＺｈｄａｎｏｖＭＳ，ＢｏｏｋｅｒＪＲ．Ｕｎｄｅｒｇｒｏｕｎｄｉｍａｇｉｎｇｂｙｅｌｅｃｔｒｏ ｍａｇｎｅｔｉｃｍｉｇｒａｔｉｏｎ［Ｊ］．Ｉｎ：６３ｒｄＡｎｎ．Ｉｎｔｅｒｎａｔ．Ｍｔｇ．Ｅｘｐｌ．Ｇｅｏｐｈｙｓ．，ＥｘｐａｎｄｅｄＡｂｓｔｒａｃｔｓ．１９９３，３５５～３５７．［３９］ ＺｈｄａｎｏｖＭＳ，ＬｉＷＤ．２Ｄｆｉｎｉｔｅ ｄｉｆｆｅｒｅｎｃｅｔｉｍｅｄｏｍａｉｎｅｌｅｃ ｔｒｏｍａｇｎｅｔｉｃｍｉｇｒａｔｉｏｎ［Ｊ］．Ｉｎ：６７ｔｈＳＥＧＥＭ２．１．１９９７，３７０～３７３．［４０］ ＤｅＨｏｏｐ，Ａｄｒｉａｎ’ｓＴ．Ｔｒａｎｓｉｅｎｔｅｌｅｃｔｒｏｍａｇｎｅｔｉｃｓｅｉｓｍｉｃｐｒｏｓｐｅｃｔｉｎｇ ａｃｏｒｒｅｓｐｏｎｄｅｎｃｅｐｒｉｎｃｉｐｌｅ［Ｊ］．ＧｅｏｐｈｙｓｉｃａｌＰｒｏｓｐｅｃｔｉｎｇ．１９９６，４４（６）：９８７～９９５．６１８。
用瞬变电磁法识别地下电性分界面
文章编号:100121986(2003)0520049202用瞬变电磁法识别地下电性分界面薛国强,宋建平,李 貅(西安交通大学电子与信息工程学院,陕西西安 710049)摘要:提高瞬变电磁法资料的解释精度一直是该方法推广应用的关键。
文中提出了基于反射系数的瞬变电磁地下介质电性分界面的识别方法,并给出了数学模型。
由所测磁场数据在频率域内数字滤波技术计算波阻抗,再通过域的变换在时间域内建立含有反射系数的线性方程组,并解出各时间道所对应的反射系数,最终以此值为参数绘制剖面图。
两个模型和一个实测数据的数值模拟计算结果表明,本文提出的方法正确、有效,并大大提高了瞬变电磁对地下电性层的分辨能力。
关 键 词:反射系数;瞬变电磁;电性界面;识别中图分类号:P63113+25 文献标识码:A1 引言瞬变电磁法(TE M)具有便携带,易操作,采集快,探测深度变化范围大,生产成本低,对地形条件要求少等特点。
在某些特殊情况下,如弱地震反射区、地形复杂地区等,可弥补地震勘探的不足。
提高瞬变电磁法的解释精度,扩大其应用范围,一直是地球物理工作者不断追求的目标[6],也是地质、工程等给该方法提出的急待解决的课题。
传统的瞬变电磁法资料解释手段单一,主要是通过视电阻率值这一参数来对地下电性分布情况作出判断,这往往使人为因素增大。
所以,有必要寻求一种更直观地反映地下介质分层情况的多元参数。
本文就是在分析瞬变电磁在水平层状介质下电磁响应的特性的基础上,提出了以反射系数为特征量,进行瞬变电磁电性界面的识别。
即从野外实测的电压衰减数据进行一系列计算,求出反射系数序列,并以此为参数绘制断面图,从而判定地下电性界面。
从对两个三层模型和实测资料的处理结果看,用文中提出的方法绘制的剖面清晰直观地显示出了地下电性层分界面。
2 基本理论[1,4,5,6]瞬变电磁场的求解是通过在频率域中将麦克斯韦方程组转化为亥姆霍兹方程求解,然后利用域的转换技术将频率域的结果转换到时间域来。
TEM瞬变电磁法简述
TEM瞬变电磁法简述瞬变电磁法或称时间域电磁法（Transient Electromagnetic Method，简称TEM），是以地壳中岩（矿）石的导电性与导磁性差异为主要物质基础，根据电磁感应原理，以不接地回线（磁偶源）向被测地质体发射脉冲式电场作为场源（一次场）。
以此来激励地下介质的二次涡流场，并对二次场进行观测。
在发射脉冲的间隙利用接收回线（线圈）接收二次场，通过分析二次场随时间的变化特征，来获取地下介质的电性特征（电阻率），推断目标体的空间赋存位置、产状、埋深等信息。
图1瞬变电磁法原理图如图1所示，在地面布设发送回线，并给发送回线上供一个电流脉冲方波，在一次磁场的激励下，地质体将产生涡流，在一次场消失后，该涡流不会立即消失，它将有一个过渡(衰减)过程。
该过渡过程又产生一个衰减的二次磁场向地表传播，在回线一定范围内接收回线接收二次磁场。
1.2 TEM如何实现测深在瞬变过程早期阶段，高频谐波占主导地位。
由于高频的趋肤效应，涡旋电流主要集中在导电介质的表层附近且阻碍电磁场向地质体深处传播。
所以早期阶段主要反映地质体断面上部地质信息。
随着时间的推移，高频成分被导电介质吸收，从而低频成分占主导地位。
它在导电地质体中激发出很强的涡旋电流。
然而由于热损耗，这些涡旋电流场很快就消失了。
在瞬变过程的晚期，局部地质体中的涡流实际上全部消失，而在各个地层中的涡流磁场之间连续的相互作用使场均匀化和使电流均匀分布，晚期场将依赖于断面的总纵向电导。
1.3 TEM如何探测地质体信息在发送一次脉冲磁场的间歇期间，观测由地质体受激励引起的涡流产生的随时间变化的感应二次场的强度。
地质体介质被激励所感应的二次涡流场的强弱决定于地质体介质所耦合的一次脉冲磁场磁力线的多少，即二次场的大小与地下介质的电性有关：（1）低阻地质体感应二次场衰减速度缓慢，二次场电压较大；（2）高阻地质体感应二次场衰减速度较快，二次场电压较小。
根据二次场衰减曲线的特征，就可以判断被测地质体的电性、性质、规模和产状等，由于瞬变电磁仪接收的信号是二次涡流场的电动势（即二次电位），因此，瞬变电磁作为一种时间域的人工源地球物理电磁感应探测方法，是根据地质构造本身存在的物性差异来间接判断相关地质现象的一种有效的地质勘探手段。
瞬变电磁
瞬变电磁法1、概述顺便电磁法（TEM）属于时间域电磁法，它是的原理是根据地壳中岩石或者矿体的导电性及介电性等电学性质的差异，以不接地的回线或者是连接地线通上脉冲电流为场源，地下发射一次脉冲磁场，在一次脉冲磁场间歇期间利用线圈或接地电极观测地下介质中引起的二次感应涡流场，从而探测介质电阻率的一种方法。
其基本工作方法是：于地面或空中设置通以一定波形电流的发射线圈，从而在其周围空间产生一次电磁场，并在地下导电岩矿体中产生感应电流：断电后，感应电流由于热损耗而随时间衰减，有一个瞬变的过程。
可以通过判断和分析二次的时空变化特征，来判断地下地质体的电性特征，找出其位置，产状和埋深等特征。
具有可以同时的具有时间和空间的可分性、探测深度达、分辨率高、信息丰富等优点。
近几十年来，我国科学技术快速进步，经济迅猛发展，各项基础建设稳步展开，对于各种矿产资源、能源、地下水资源等的需求快速增加。
同时，各项建设中遇到了许多工程问题，如公路建设中的地下空洞、煤田开采中的陷落柱、隧道开挖中的突水问题等等。
这些因素在一定程度上制约着我国经济的发展，而顺便电磁法的出现，利用其测量方面的优势，已经发展成为探测油气、金属和非金属矿产的一种重要方法，并且在深部地质构造研究，工程勘察、油气、矿产、水、地热勘探等领域得到了广泛的应用。
可以很好地保证资源供给，减少经济损失，加快建设进度。
2、研究现状2.1、研究历史对勘测工程工作的种种困难，把瞬变电磁法应用到地质勘探中的想法在上世纪30年代就有人提出来。
最初的时域电磁法是利用到了L.W.Blan在1993年获得专利，用电磁脉冲激发提供电偶极形成电场。
随后在前苏联有人提出了瞬变电磁测深法。
在50年代，前苏联、加拿大、美国等国已经开始就瞬变电磁法的理论与应用技术进行了深入的研究，同时期由J.R.Wait 提出了使用瞬变电磁场法寻找导电矿体的理念。
前苏联也基本已经建立了瞬变电磁法与野外施工的技术方法，更在70、80年代开展了大量的测量工作，特别是在二维和三维测量的方面就有了很大的进步，这使的瞬变电磁法在地质勘探上运用有了很大的发展。
瞬变电磁法理论的应用现状及研究进展初探
瞬变电磁法理论的应用现状及研究进展初探摘要：随着改革开放的不断深化，科技水平得到了极大地提高，作为时间域电磁法的一种表现形式，瞬变电磁法在我国多个领域得到了广范的应用。
文章主要对瞬变电磁法理论应用的现状和研究进展进行了简要分析。
关键词：瞬变电磁法；含义；现状；进展瞬变电磁法诞生于20世纪50年代的西方，经过多个国家地球物理学家对其理论、技术等多方面的深入分析探究，尤其是20世纪的70、80年代，前苏联在这方面开展了大量的试验、测量工作。
80年代后，瞬变电磁法得到了广泛的应用，目前主要应用的领域有：油气勘测、矿产勘查、工程勘察、环境调研、考古探测、军事探测等。
随着科技的不断进步，设备种类也越来越多，越来越细化，推出了适应不同领域的不同设备类型。
1 瞬变电磁法的含义瞬变电磁法是时间域电磁法的重要组成部分，它的生成原理主要是通过不接地线或接地线源向地下传送一次脉冲电磁场，在一次电磁场的作用下，地下导体内部受感应产生涡旋电流，然后在一次脉冲磁场的间隙期间，涡流电流产生的二次磁场不会随着一次磁场的消失而马上消失，继而产生了一个瞬变过程。
通过线圈或接地电极来观测二次磁场，进而观察它与时间之间的关系，得出结论判断地下导体的电性分布结构和空间形态。
2 瞬变电磁法理论的应用现状从工作时间上进行划分，可以将瞬变电磁法分为四种：地面瞬变电磁法、地下瞬变电磁法、航空瞬变电磁法、海洋瞬变电磁法。
在我国瞬变电磁法的应用过程中，主要使用的是地面瞬变电磁法。
2.1 瞬变电磁法在地面勘测中的应用现状瞬变电磁测深法具有能够近区观测、直接测二次场、工作效率高及受静态位移影响小的特性。
随着科技的发展，瞬变电磁仪器智能化、数字化技术水平得到了极大的提升。
瞬变电磁法也不在仅仅局限于金属矿勘测，还广泛地应用于工程建筑、环境等方面的地质研究调查，如对地下采空区、陷落柱等煤田灾害的探测，油气田的勘测等都取得了很好的成效，在地面勘测中除了瞬变电磁法之外，还应用了磁源重迭回线法、大回线源法等。
简述瞬变电磁法(TEM)的研究现状、进展与问题
摘 要 : 近年来 ,由于电子技术的快速发展 、勘查市场的旺盛需求以及瞬变电磁法 ( TEM )的自 身优点 ,使得 TEM 方法从仪器到方法技术均获长足进展 。在此 ,简述了 TEM 仪器及方法技术的 研究现状与进展 ,总结了 TEM 的应用领域 ,分析了 TEM 方法目前存在的主要问题 ,探讨了 TEM 急需解决的问题及发展方向 。
·110·
物探化探计算技术
29卷
[ 11 ] 熊彬 ,罗延钟 ,强建科 . 瞬变 电磁 2. 5 维反演 中灵敏 度矩阵计算方法 ( Ⅰ) [ J ]. 地 球物理 学进展 . 2004, 19 ( 3) : 6161
[ 12 ] 雷达 . 瞬变电磁法数据提取纯异常 方法技术 [ J ]. 物 探与化探 , 2004, 28 ( 4) : 3201
[ 8 ] 李貅 , 郭文 波 , 胡建平 . 瞬 变电磁 测深 快速 拟地 震解 释方法及应用效果 [ J ]. 西安工程学院学报 . 2001, 23 (03) : 42.
[ 9 ] 罗四维 . 人 工神经 网络 建造 [ M ]. 北京 :中 国铁 道出 版社 . 2000: 103.
[ 10 ] 陈 易玖. 小框瞬变电磁法 ( TEM )的勘探 深度 [ J ]. 湖 北 地矿. 2003, 17 ( 4) : 38.
[ 13 ] 薛云峰 ,孙晓暾 ,孙雅 芳. 瞬变 电磁电 阻率成 像法探 测岩体渗漏 通道的 研究 及应用 [ J ]. 水 文 地质 工程 地质 . 2004, 31 ( 3) : 961
[ 14 ] 郭 文波 ,宋建平 ,李正宇 . 小波 分析在 T EM 资料处 理中的应用 [ J ]. 矿产与地质 , 2006 , 20 (4～5).
勘测 ,风化层厚度勘查 ,水文地质分层等。 (2)断裂调查 : 基岩断裂及活动断裂勘查 ; 环
瞬变电磁法基本理论
瞬变电磁法基本理论发表时间：2018-12-12T16:56:08.060Z 来源：《基层建设》2018年第29期作者：刘晨阳[导读] 瞬变电磁法[1][2]（Transient Electromagnetic Methods）的基础是电磁感应原理，场源为人工源，因为研究的是响应场与时间的关系，又被称为时间域电磁法（Time Domain Electromagnetic Methods）。
吉林大学建设工程学院工程地质系长春 1300261瞬变电磁法的原理瞬变电磁法[1][2]（Transient Electromagnetic Methods）的基础是电磁感应原理，场源为人工源，因为研究的是响应场与时间的关系，又被称为时间域电磁法（Time Domain Electromagnetic Methods）。
其人工场源分为2类：电偶源（即接地回线）和电磁源（即不接地回线）。
利用人工场源向地下发射一次脉冲场，在其激发下的瞬间，产生一个向回线法线方向传播的一次磁场（即一次场），在一次场激励下，地质体将产生涡流。
在一次场消失以后，涡流不会马上消失，它会有一个衰减的过程，此过程会产生一个衰减的二次磁场，并继续传播，再由接收回线接收二次场。
这样，通过分析二次磁场的信息变化，就可以得到地质体的电性分布情况。
图1 瞬变电磁法的原理图[3]对于层状大地的瞬变电磁场，在垂直磁偶源的情况下，其均匀半空间的计算公式如下[4]：1.1烟圈效应针对瞬变电磁法激发的磁场，美国科学家M.N.Nahighian 提出[4][5]：感应涡流场在地表的磁场值为地下各个“环带”涡流层的总效应，该效应等效于一个电流环。
人们形象的将其称为“烟圈效应”。
涡流的形状为向下、向外扩散的同心环状，涡流与地面成47°扩散，其密度最大值的扩散方向与地面夹角大概30°。
其扩散公式如下[7]:其原因可以做如下解释：从上面的原理部分我们可以知道，在一次场激励下，地质体将产生涡流，涡流会引发一个磁场。