瞬变电磁法拟地震成像研究进展

瞬变电磁法是工程地球物理勘察中应用较多的一种勘探方法之一,瞬变电磁法拟地震成像方法研究是当前电磁探测理论与应用研究的热门.文中详细叙述了瞬变电磁法拟地震成像的两种不同思路,一种是建立时间—频率转换经验公式,把TEM数据转换成平面波数据,借鉴大地电磁法拟地震解释思路进行成像,另一种是通过积分变换,把扩散场变成虚拟波场,把电磁衰减信号变成拟地震子波,然后按照地震解释法进行解释,并分析了各自存在的的问题.指出下一步研究的方向.     

电性源瞬变电磁地空逆合成孔径成像

电性源地空瞬变电磁法具有工作效率高、勘探深度大、采集信号信噪比高、适用于地形地质条件复杂地区等优点.但是,到目前为止,由于尚未建立起该方法的解释系统,大大制约了该方法的发展.本文旨在建立起完整的地空电磁探测系统,丰富整个探测系统的理论.本文围绕地空瞬变电磁法全域视电阻率定义、瞬变电磁虚拟波场的克希霍夫偏移成像、逆合成孔径成像方法三个科学问题进行了系统研究.提出了用磁场强度定义全域视电阻率的迭代算法,理论模型试验结果表明计算出的视电阻率曲线首支趋于第一层电阻率,尾支趋于最后一层电阻率,实现了全空域、全时域视电阻率的计算;在先前研究的基础上,实现了适合电性源地空装置的瞬变电磁虚拟波场的克希霍夫偏移成像;采用相关迭加技术,实现孔径内多测点数据合成,将传统的单点处理方式发展成为逐点推移多次覆盖的逆合成孔径处理方法.层状模型试验表明:(1)全域视电阻率能够光滑、完整、渐变地反映出模型的电性信息变化;(2)当改变三层模型中间层电阻率时,全域视电阻率曲线随着参数的改变分异明显,对电性层的识别容易且直观;(3)由于在电阻率计算中同时考虑了接收机高度、偏移距、时间等各参数的影响,全域视电阻率可实现全空域、全时域的视电阻率计算.含水采空区的复杂模型算例表明:(1)根据不同测线的全域视电阻率结果可以看出,在靠近采空区的位置,全域视电阻率断面可以清晰地反映出采空区的空间位置,随着测线离采空区越来越远,采空区异常越来越弱直至消失;(2)波场变换和偏移成像的结果显示存在两个电性差异较大的界面,上界面指示地表,由于空气和大地之间的电性差异较大,故该界面波场信号反映强烈,遍布整个区域,下界面异常信号则主要集中在中部,向外逐渐减弱,指示采空区;(3)逆合成孔径成像结果表明地表界面在合成前后没有变化,而采空区异常合成后范围明显变小,且异常边界清晰,指示的采空区位置与模型吻合很好.本文借助于逆合成孔径雷达成像的基本思想,建立了一套电性源瞬变电磁地空逆合成孔径成像方法.基于反函数思想结合迭代算法提出的电性源地空瞬变电磁法的全域视电阻率定义方法,实现了全空域、全时域的视电阻率计算;借鉴瞬变电磁拟地震偏移成像算法,实现了瞬变场的三维成像;借鉴逆合成孔径雷达的思想,提出电性源瞬变电磁地空逆合成孔径算法,进一步提高了成像的分辨率.采空区模型算例表明相关叠加合成确实具有增强有用信号、提高信噪比、提高分辨率的诸多优点,证实了瞬变电磁地空逆合成孔径成像方法的有效性.     

瞬变电磁拟地震子波宽度压缩研究

导电媒质中电磁场所满足的扩散方程与无耗媒质中电磁场所满足的无阻尼波动方程之间存在着数学转换关系式,可以实现瞬变电磁法拟地震资料解释.但是转换出的虚拟波往往波形较宽,使分辨能力下降,影响了瞬变电磁拟地震方法的发展和推广.本文针对虚拟波形展宽这一现象进行了研究,认为虚拟波形展宽的原因,不是由于波在介质中传播时能量损耗所致,...     

电性源瞬变电磁有限差分偏移成像方法

目前,瞬变电磁数据的主要处理结果是给出地下介质的（视）电阻率信息,较难给出清晰的地质界面信息.本文基于有限差分偏移成像原理,实现了电性源瞬变电磁虚拟波场的偏移成像技术.根据瞬变电磁扩散场与虚拟波场之间的数学关系,本文使用精细积分法对第一类Fredholm积分方程求解,然后消除直达波,得到只含有地层反射波的虚拟波场.在此基础上,本文从时间域波动方程出发,对上行波的一级近似方程反向外推,得到时间-空间域二维与三维的波场递推公式.通过求解地下空间的波场,并使用反褶积方法削弱虚拟波场的波形展宽效应,从而提高对地质界面的纵向分辨率,获得较为准确的地质界面信息.本文分别采用均匀半空间模型、D型模型、HK型模型以及两层含水采空区模型对本文方法进行验证,并对实测数据进行处理,偏移成像的结果与地质资料相吻合,说明瞬变电磁有限差分偏移成像方法在识别电性界面以及地层特征方面具有可行性.

     

瞬变电磁虚拟波场高分辨成像技术研究

本文提出了从瞬变电磁场到波场的优化算法.目的是通过对现有渡场变换过程进行插值的办法提高波场变换的分辨率,从而为后期数据处理的质量打好坚实的基础.在进行波场反变换的过程中,将正则化算法用于逆变换过程的计算中,通过采用偏差原理和Newton迭代格式选出最优的正则化参数,使得反变换所得到的波场稳定、可靠.解决了积分系数过多而...     

瞬变电磁虚拟波场的三维曲面延拓成像

根据瞬变电磁场扩散方程与波动方程间存在的数学对应形式,将已知瞬变电磁场数据转换为虚拟波场数据,在此基础上,提出了用克希霍夫积分法将所转换的波场从地面向地下反向延拓的方法.在延拓的数值计算中,为提高计算效率,采用三维边界元技术,把边界积分分解为诸三角单元积分的积分之和, 建立了曲面延拓方程式;通过地面测点的波场值,求出地下某一点的波场值,实现了地表为曲面的向下延拓成像计算.通过对理论模型计算和实际资料处理,证明了该方法可以增强瞬变电磁法识别地下电性分界面的能力,使瞬变电磁法对地下目标体的三维精细探测成为可能.     

矿井瞬变电磁法超前探测合成孔径成像研究

在矿井瞬变电磁超前探测数据处理和解释中,为了更好地识别巷道前方含水体电性分界面,利用扩散场与波动场间的函数关系,将矿井瞬变电磁法超前探测中瞬变电磁波场转换为拟地震波场,实现波场转换.同时,为了解决波场转换后信号幅值偏弱,以增加信噪比利于解释,利用雷达的合成孔径方法对波场转换后的拟地震信号进行相关叠加处理.将矿井瞬变电磁法巷道超前探测的同一测点不同探测方向的数据作为孔径数据进行合成孔径成像,以突出弱异常而提高信噪比,结合应用实例验证了方法的效果.研究表明,矿井瞬变电磁法合成孔径成像方法能够提高电性界面信息的信噪比和界面几何分辨率,对矿井瞬变电磁法精细超前探测含水构造具有重要的理论意义和应用价值.     

基于高性能瞬变电磁辐射源的城市地下空间多分辨成像方法研究

城市地下空间的高精度、多分辨探测是城市地下合理开发的前提.由于城市探测环境的特殊性,需要探测方法具备抗干扰能力强、分辨能力好、分异多尺度地下目标等特性,现有的物探方法难以兼顾城市地下空间的探测需要.本文利用三维矢量有限元方法为正演手段,在瞬变电磁高性能辐射源进行微分脉冲扫描的基础上,对探测数据进行多时窗的扫时波场变换,将微分脉冲扫描后的多分辨响应信息,进行多分辨信息提取.同时通过地震探测中的多次覆盖处理,提高探测场对地下目标的分辨能力,最后对多次叠加后的虚拟波场进行拟地震偏移成像,最终实现城市地下空间的高精度探测.

     

基于高性能瞬变电磁辐射源的城市地下空间多分辨成像方法研究

城市地下空间的高精度、多分辨探测是城市地下合理开发的前提.由于城市探测环境的特殊性,需要探测方法具备抗干扰能力强、分辨能力好、分异多尺度地下目标等特性,现有的物探方法难以兼顾城市地下空间的探测需要.本文利用三维矢量有限元方法为正演手段,在瞬变电磁高性能辐射源进行微分脉冲扫描的基础上,对探测数据进行多时窗的扫时波场变换,将微分脉冲扫描后的多分辨响应信息,进行多分辨信息提取.同时通过地震探测中的多次覆盖处理,提高探测场对地下目标的分辨能力,最后对多次叠加后的虚拟波场进行拟地震偏移成像,最终实现城市地下空间的高精度探测.     

弹性波克希霍夫积分偏移法

P、S波同时偏移法是为了适应复杂构造地区勘探情况而开始研究的一种偏移方法。在前人工作的基础上,应用弹性波克希霍夫积分偏移的基本原理,改进了纵波（PP）及转换波（PS）的同时克希霍夫积分偏移方法,给出了一种近于实用的计算机流程和程序;并用其分别对射线理论合成反射波资料、有限元模拟地震波传播资料进行了检验,证明此种方法是有效的。在解决复杂构造地区的地震波成象方面,该方法具有很大的潜力。     

等值反磁通瞬变电磁成像技术及工程应用

为了实现对可溶岩区公路路线范围内的隐伏岩溶发育规模及形态进行高精度勘察,为后期处治设计提供翔实可靠的地质依据,在公路塌陷区域岩溶勘察中引入了一种新技术——等值反磁通瞬变电磁成像技术.文中首先对等值反磁通瞬变电磁成像技术的基本原理进行简述;然后论述了该方法的4种主要技术特点;紧接着对该方法的质量控制原则及数据处理流程等关键技术进行分析;最后利用该方法在惠清高速公路岩溶塌陷区探测中进行实际应用,成像结果准确地反映出了岩土界面的位置及变化形态,精确探测出了岩溶发育区的形态及规模,为后期塌陷区域的处治设计提供了可靠的基础资料.通过理论分析及实际应用结果表明:等值反磁通瞬变电磁法可对公路穿过区域隐伏岩溶发育情况进行精确成像,可快速准确的为岩溶塌陷区的后期处治设计提供基础资料.     

瞬变电磁BEDS-FDTD三维正演新算法及稳定性验证

常规FDTD的时间步长需满足Courant-Friedrich-Lewy（CFL）稳定性条件,导致时间迭代步过多,正演非常耗时.针对上述问题,本文采用Backward Euler（BE）差分方法近似Maxwell方程组中场对时间的一阶导数,使时间步长突破CFL限制;引入Direct Splitting（DS）策略将电磁场分量解耦,并将大型稀疏矩阵降阶和重构为一系列低阶且主对角占优的三对角矩阵,加快方程求解效率;而且,为了减少模型尺寸,针对上述改进的全新方程,本文采用双线性变换方法（bilinear transform,BT）推导了复频率偏移完全匹配层（complex frequency shifted perfectly matched layer,CFS-PML）吸收边界,形成瞬变电磁三维正演全新算法：BEDS-FDTD.首先采用von Neumann方法测试了新算法在有耗介质中、非均匀时间步长下的稳定性.之后,将BEDS-FDTD算法用于模型实验,并将模拟结果与层状模型的半解析解进行了对比,实验结果证明了新算法的计算精度能够满足需要;同时对新算法的计算效率进行了分析,当采用Tesla A100 GPU计算时,50×50×50网格数目的模型仅用10 s,即使网格数目增加到200×200×200,也仅需224 s.最后,将BEDS-FDTD算法用于复杂三维模型计算.

     

3-D approximate inversion for transient electromagnetic soundings

IntroductionThe issues of forward and inverse interpretation for homogeneous half space models and layered models of the Earth have well been solved in trAnsient electromagnetic sounding (TEM), butinterpreting data from more complicated 2-D or 3-D models are difficult because their responsefunctions are complicated and expensive to calculate. For the sake of simplification, based on thetheories of smoke ring (Nabighian, 1979), we regard current lines in the Earth induced by transmitting sou…     

On the exploration of a marine aquifer offshore Israel by long-offset transient electromagnetics

The existence of aquifers extending from land beneath the sea floor up to a distance of several kilometres has been observed and examined all over the world. The coastal aquifer of Israel is a heavily used groundwater reservoir which has to be constantly monitored to ensure the drinking water supply. Former land-based electromagnetic measurements show that it is, in several places, blocked to seawater intrusion and is consequently a candidate for submarine extension. Multicomponent long-offset transient electromagnetic measurements were carried out offshore on the coast of Israel. We deployed a 400-m-long grounded dipole as transmitter and several electric and magnetic receivers on the sea floor up to a distance of 4.8 km from the coast. Altogether, we deployed 8 transmitter positions and received data sets at 14 receiver stations onshore and offshore, with offsets of mostly 400 and 800 m. In this paper, we present the survey and 1D Occam and Marquardt inversions of the offshore horizontal electric components in the broadside and inline configuration. In addition, the vertical magnetic component in the broadside position is also considered. Only single inversions, both single offset and single component, were used to detect the aquifer under sea sediments. We prove the submarine existence of the Israeli coastal aquifer up to a distance to the coast of approximately 3.2 to 3.6 km using all measured long-offset transient electromagnetic components. In addition, we present modelling studies with synthetic data derived from a subsurface model adjusted to our survey area with very shallow water from 10 to 50 m. As part of the planning before the survey, a parameter study of the expected subsurface, the examination of the airwave phenomenon and the justification for our 1D inversion strategy are shown. More detailed eigen parameter analyses are added to explain the measured data.     

Integrating long-offset transient electromagnetics (LOTEM) with seismics in an exploration environment1

The applications of electromagnetics have increased in the past two decades because of an improved understanding of the methods, improved service availability, and the increased focus of exploration in the more complex reservoir characterization issues. For electromagnetic methods surface applications for hydrocarbon Exploration and Production are still a special case, while applications in borehole and airborne research and for engineering and environmental objectives are routine. In the past, electromagnetic techniques, in particular deep transient electromagnetics, made up a completely different discipline in geophysics, although many of the principles are similar to the seismic one. With an understanding of the specific problems related to data processing initially and then acquisition, the inclusion of principles learned from seismics happened almost naturally. Initially, the data processing was very similar to seismic full-waveform processing. The hardware was also changed to include multichannel acquisition systems, and the field procedures became very similar to seismic surveying. As a consequence, the integration and synergism of the interpretation process is becoming almost automatic. The long-offset transient electromagnetic (LOTEM) technique will be summarized from the viewpoint of its similarity to seismics. The complete concept of the method will also be reviewed. An interpretation case history that integrates seismic and LOTEM from a hydrocarbon area in China clearly demonstrates the limitations and benefits of the method.     

利用三维高斯射线束成像进行地震定位

常规的地震定位方法通常需要拾取地震记录的初至,当初至不明显或被较高水平的噪声淹没时精度较低.本文采用基于三维高斯射线束的偏移成像方法对震源进行定位,较好地解决了该问题.通过三维高斯射线束对台站记录进行偏移归位,并将各台站成像结果的交点作为地震能量释放的中心位置;当各台站成像结果不能交于一点时,采用三维空间高斯滤波方法可实现震源位置的自动获取.提出的变网格计算方案极大地减少了计算量,显著地提高了成像精度和计算效率.利用首都圈地震台网数据,对涿鹿、滦县以及房山三个地震事件进行试算,结果表明:基于变网格三维高斯束偏移成像的地震定位方法自动化程度很高,而且具有较好的抗噪能力,特别适合处理低信噪比资料的地震定位问题.     

Deep exploration using long-offset transient electromagnetics: interpretation of field data in time and frequency domain

In the framework of the Deep Electromagnetic Soundings for Mineral Exploration project, we conducted ground-based long-offset transient-electromagnetic measurements in a former mining area in eastern Thuringia, Germany. The large-scale survey resulted in an extensive dataset acquired with multiple high-power transmitters and a high number of electric and magnetic field receivers. The recorded data exhibit a high data quality over several decades of time and orders of magnitude. Although the obtained subsurface models indicate a strong multi-dimensional subsurface with variations in resistivity over three orders of magnitude, the electrical field step-on transients are well fitted using a conventional one-dimensional inversion. Due to superimposed induced polarization effects, the transient step-off data are not interpretable with conventional electromagnetic inversion. For further interpretation in one and two dimensions, a new approach to evaluate the long-offset transient-electromagnetic data in frequency domain is realized. We present a detailed workflow for data processing in both domains and give an overview of technical obstructions that can occur in one domain or the other. The derived one-dimensional inversion models of frequency-domain data show strong multi-dimensional effects and are well comparable with the conventional time domain inversion results. To adequately interpret the data, a 2.5D frequency-domain inversion using the open source algorithm MARE2DEM (Modeling with Adaptively Refined Elements for 2-D EM) is carried out. The inversion leads to a consistent subsurface model with shallow and deep conductive structures, which are confirmed by geology and additional geophysical surveys.     

重力和重力梯度数据三维相关成像

本文提出了重力异常三维相关成像方法和重力梯度数据三维相关成像方法,并提出了基于异常分离的三维相关成像方法来提高成像分辨率.通过合成Y型岩脉模型和合成多个直立长方体组合模型的重力异常和重力梯度数据试验分析,验证了本文三维相关成像方法可显示出异常地质体的空间赋存状态和等效剩余质量分布,具有良好的纵向和横向分辨率.