一种多尺度电磁测深方法

电磁测深方法是地球物理勘探的重要方法.按照处理方式分为时间域和频率域电磁测深方法两类.但实质上时间域方法和频率域方法是等效的.而由于观测方式的不同,有时效果不能完全等同.本文介绍一种广义的时间域电磁测深方法.这种方法的最大优点就是在不改变激励和观测装置的前提下,通过改变激励波形,从而实现分段多尺度的电磁测深处理.     

编码电磁测深

利用逆重复M序列伪随机信号良好的自相关特性,提出相关检测地电传输特性的编码电磁勘探法.分析了编码电磁测深原理及参数提取方法.同时记录发送电流信号和多收发距电磁场响应,对源信号和场信号进行相关运算,解卷积分离接收系统响应后可得到大地的频率特性或时间特性,以此实现地电断面的精细探测.在时间域,大地冲激响应和阶跃响应含有丰富的地电信息.通过冲激响应的峰值时间或阶跃响应的晚期渐近值可估计地电阻率分布.基于层状模型的大地冲激响应和阶跃响应正演计算结果表明,编码电磁测深法对大埋深薄层目标体有精细的分辨能力.可应用于能源、矿产资源、水资源、环境地质及工程地质勘察,有良好的应用前景.     

Use of marine controlled-source electromagnetic sounding for sub-basalt exploration

The problem of imaging sedimentary structure obscured by high-velocity layers, such as carbonate, basalt or salt, using conventional seismic techniques is well known. When this problem is encountered in offshore areas, marine electromagnetic data can provide valuable, alternative and complementary constraints on the structure. We concentrate on the use of marine controlled-source electromagnetic (CSEM) sounding in the frequency domain. The feasibility of mapping sub-basalt sedimentary structure using this technique is explored by means of modelling studies which mimic a type of survey which could readily be employed in practice. As a starting point the problem is addressed in terms of 1D resistivity structure. We show that sub-basalt sediments can be detected and their depth of burial quantified to within 200 m in the examples shown. The presence of small-scale inhomogeneities in the basalt (which cause much of the scattering in seismic data) is shown to have no appreciable effect on the ability of the CSEM data to detect the sediments. The modelling is then extended to 2.5 dimensions. Again the presence of sub-basalt sediments can be readily detected and their properties and geometry quantified, especially if the electromagnetic data are combined with constraints from complementary geophysical techniques.     

从瞬变电磁测深数据到平面电磁波场数据的等效转换

中心回线源装置的瞬变电磁场是一种涡流场,在地下主要以扩散形式传播.平面电磁波在地下介质中传播时具有反射、折射等电磁学特性,为了借鉴平面电磁波场测深成熟的解释方法,文中开展了从瞬变电磁涡流场测深数据到平面波场数据的快速等效转换的研究.通过大量的理论模型的正演计算、曲线对比,误差统计,对瞬变电磁场、平面波场在地下传播机制特性的分析,建立了一种由瞬变电磁测深视电阻率数据向平面波场测深视电阻率数据转换的时间-频率等效对应关系.利用这一关系式,把扩散场数据快速转换成等效平面波场测深数据. 为进一步进行瞬变电磁拟平面波成像解释打下理论基础.     

天然电磁辐射测深技术的应用

天然电磁辐射测深技术是一项利用天然电磁场探测地球介质相对电阻率与介质埋深的地球物理探测技术．此项技术主要探测在天然电磁场作用下在近地表形成的综合电场的垂直分量,根据其频谱的局部变化构制岩层相对电阻率与埋藏深度的直方图,进而对地下岩层的性质和分布特征作出合理的分析和判断．本文主要以探测实例说明本项技术的应用效果及其适用的范围．     

Limiting apparent-resistivity values for dipping-bed earth models

Consideration is given to some special features of normalized apparent-resistivity (NAR) curves, resulting from the deployment of Wenner electrode arrays on the surface above dipping earth structures. Limiting values of the potential are derived when a direct-current source is located at points on the surface above a simple two-region dipping-bed earth model and their influence on corresponding NAR curve characteristics is investigated. Particular attention is given to the exploitation of such features to provide a new and direct approach to model parameter estimation, either as an alternative to traditional curve-matching techniques or as a source of supporting information when other earth model characterization methods have been employed. Throughout, the emphasis is on the single-dipping-bed model, but application to more complex structures is discussed, including examples of two dipping beds, dipping dikes and more general tilted unconformities.     

一种实用的时间域电磁测深解释方法

电磁测深方法,特别是可控主动源方式的电磁测深方法成了研究地球深部构造的一种非常有效的方法.由于电磁波在导电介质中传播时,其幅度是按照指数规律衰减的,在电磁测深的资料处理与解释中,往往反映地下深部信息的数值非常弱小,无法直观的反映深部构造信息.我们尝试一种全新的电磁感应模型来对时间域的电磁测深进行理论模拟,并对实际测深资料进行解释,取得了非常理想的效果.说明我们的模型假设是成功的.     

Results of electromagnetic sounding with JASC submarine cable

We present the results of long-term deep geoelectric studies using the JASC (Japan Sea Cable) submarine communication cable in the region of the Sea of Japan. In the 2D inversion of the amplitude and phase’s apparent resistivity curves and the frequency dependences of the tipper, we invoked the geological and geophysical information about the region and on-shore electromagnetic observations to fit the model to the observations. The resulting geoelectrical cross section of the region of the Sea of Japan along the JASC cable obtained in this way agrees well with the experimental data. The upper part of the section contains a conductive block beneath the bottom of the Central basin of the Sea of Japan at a depth of 10–40 km, a fault submerging below the continent in the marginal part and a deep fault in the continental region. In the lower portions of the cross section, the high-resistivity block interrupts the continuity of the horizontal conductive layers beneath the Yamato Uplift, and the conductive bottom part of the geoelectric cross section submerges under the continent. In the continental segment of the cross section, there is a large block with reduced electric resistivity, which is located between the conductive layers at a depth interval of 200–560 km. We analyze the characteristic features of the geoelectric cross section and the deep section imaged by seismic tomography in the region of the Sea of Japan.     

频率域2维电磁测深资料的RRI反演

本文将大地电磁反演中的RRI方法应用于线源频率测深模型资料的反演中.当背景电导率变化很小时,可用变化前的电场来近似代替变化后的电场,从而在反演方程的推导过程中将有源电磁场中的源项消掉,得到和大地电磁场相同的反演方程,使有源电磁波的反演也可以应用RRI方法.反演过程中所需要的模型资料通过有限元方法得到,该资料不需做近场校正,直接用适合于有源电磁场的RRI方法反演,避免了近场校正带来的误差.数值模型结果证明该方法是可行的.最后,用RRI方法讨论了当源和目标区间存在低阻异常体时只对目标区反演的可行性,对实际工作的解释有一定的指导意义.     

Nonlocal electromagnetic response: Implications in the regional sounding

Due regard was given to the nonlocal electromagnetic responses in the processing of the data of the International Baltic Electromagnetic Array Research (BEAR) 1998 Project. Maps of lateral distributions of the apparent electrical resistivity for different time periods were compiled. In these maps, the main tectonic boundaries between the Lapland-Kola, Karelian, and Svecofennian tectonic provinces within the Baltic Shield appear as zones with reduced apparent electrical resistivity. There is supposedly some correlation between the positions of the maxima in electrical resistivity and the maximum Moho depths. In the same way, the data from the network of geomagnetic observatories in the European region (the INTERMAGNET Project) were processed for 2004–2005, and maps of the laterally inhomogeneous distribution of the apparent resistivity in Europe have been compiled for different time periods. The correlation is observed between the positions of the minima in the electrical resistivity and the minima in the thickness of the lithosphere.     

Inversion of electromagnetic data: An overview of new techniques

This paper explores some of the newer techniques for acquiring and inverting electromagnetic data. Attention is confined primarily to the 2d magnetotelluric (MT) problem but the inverse methods are applicable to all areas of EM induction. The basis of the EMAP technique of Bostick is presented along with examples to illustrate the efficacy of that method in structural imaging and in overcoming the deleterious effects of near-surface distortions of the electric field. Reflectivity imaging methods and the application of seismic migration techniques to EM problems are also explored as imaging tools. Two new approaches to the solution of the inverse problem are presented. The AIM (Approximate Inverse Mapping) inversion of Oldenburg and Ellis uses a new way to estimate a perturbation in an iterative solution which does not involve linearization of the equations. The RRI (Rapid Relaxation Inverse) of Smith and Booker shows how approximate Fréchet derivatives and sequences of 1d inversions can be used to develop a practical inversion algorithm. The overview is structured to provide insight about the latest inversion techniques and also to touch upon most areas of the inverse problem that must be considered to carry out a practical inversion. These include model parameterization, methods of calculating first order sensitivities, and methods for setting up a linearized inversion.     

The inversion of the data of magnetotelluric sounding and electromagnetic frequency sounding with horizontal electric dipole

On the basis of the dispersion relation of magnetotelluric response functions (MTRF), a filter coefficient algorithm has been made, with which the corresponding impedance phase data can be estimated using a set of apparent resistivity data. The tests of theoretical models and observed magnetotelluric (MT) data show that this algorithm is effective. Comparing the impedance phase estimated using dispersion relation with the observed phase, it can be checked whether the dispersion relation between the observed apparent resistivities and phase data was satisfied. The use of phase data corrected using the dispersion relation in the joint inversion for MT impedance is advantageous to obtain more reliable inversion results. The problems on the one-dimensional joint inversion for the (MT) apparent resistivity and the apparent resistivity of the frequency electromagnetic sounding (FEMS) with horizontal electric dipole, whose observed frequency bands are linked up each other, are studied. The observed data of two kinds of electromagnetic (EM) methods at two sites are used to inverse, the comparison with the drilling data show the results are more reliable. To supply the phase data of FEMS using the dispersion relation, for the apparent resistivity-phase data and impedance real part-imaginary part apparent resistivities of two kinds of EM methods the imitated MT joint inversions are made, and more similar results also are obtained. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 91–96, 1993. The projects sponsored by the Chinese Joint Seismological Science Foundation.     

瞬变电磁接收装置对浅层探测的畸变分析与数值剔除

从理论上对瞬变电磁系统的接收装置频率特性进行分析,得出了瞬变信号在早期发生畸变的根本原因,揭示出接收装置在发射电流关断期间和电流关断后对一次场和二次场的影响关系,给出基于导电球体模型下接收装置的谐振频率与最小取样时间的关系图,进而讨论了接收装置对浅层目标体探测的影响.当瞬变电磁测量系统在电流关断开始时刻记录全程瞬变响应和发射电流波形时,如果接收装置的位置和谐振频率已知,就可以通过数值计算方法从根本上剔除接收装置对早期瞬变信号的影响,即使接收装置存在过渡过程,也同样可以实现近地表浅层的探测,从而缩短瞬变电磁法的浅层勘探盲区,提高近地表的探测分辨率和精度.     

瞬变电磁低频近似Maxwell方程的CPML吸收边界及施加方法

为保证瞬变电磁场晚期信号计算正确, 在时域有限差分正演中通常采用低频近似, 即将磁场散度方程显式的包含在电磁场迭代方程中, 因而H_z的求解方式与H_x、H_y不同, 导致常规CPML吸收边界无法施加.本文重新推导了CPML介质中的迭代方程, 并将其与常规介质中的迭代方程在形式上进行了统一, 提出了一种适用于瞬变电磁低频近似下三维FDTD的CPML边界条件及施加方法.首先采用均匀半空间模型来验证本方法的有效性, 发现模拟结果与解析解吻合较好, 且模拟用时为采用Dirichlet边界FDTD用时的一半; 而且, 反射误差实验验证了本文的CPML边界吸收效率较高.之后, 模拟了典型的三层模型, 发现本方法的模拟结果与线性数字滤波解吻合程度较好.最后采用本方法计算了航空瞬变电磁带地形的复杂模型, 实验证明了本文提出的CPML具有广泛的适用性.

     

电偶源频率电磁测深激发极化效应研究

本文研究含激电效应(IP)电偶源频率电磁测深问题.首先,通过对电偶源频率电磁测深电磁场分量分辨率的分析,提出利用电场Ex分量提取IP信息的合理性;然后,给出了同时存在激电和电磁效应时的理论计算公式和算法,对一些典型地电断面进行理论计算并对计算结果进行分析;最后,提出提取IP信息的几种可能方案.这些工作不仅为野外实测资料IP信息的提取提供理论依据,同时为进一步实现电磁响应中感应效应与激电效应分离提供基础.     

Direct and inverse problems in electromagnetic sounding of three-dimensional heterogeneous medium

The methods for solving three-dimensional (3D) direct and inverse problems of electromagnetic sounding are considered. It is shown that the method of integral equations is an efficient instrument for mathematical modeling of electromagnetic fields in 3D heterogeneous media. Adaptation of the integral equation technique to the solution of inverse problems is demonstrated.     

On reducing ambiguity in the interpretation of transient electromagnetic sounding data1

An attempt to resolve non-uniqueness in the interpretation of transient electromagnetic (TEM) sounding data using measured data alone is made. It is shown in the various examples studied that sufficiently early time measurements can be the determining factor in reducing the ambiguity caused by model equivalence. The early delay times thus play a dual role in transient soundings: they are responsible for resolving shallow structures and they may eliminate the ambiguity in the interpretation of geoelectric parameters of deeper targets. This is illustrated by the results of a follow-up TEM survey at the Dead Sea coast of Israel where the use of supplementary early time measurements allowed non-uniqueness in the determination of the depth to fresh/saline groundwater interface to be resolved.     

Usage of regularization techniques for effective electron collision frequency determination from oblique ionospheric sounding

Computer modelling is used to investigate the possibility of determining ionospheric parameters from slightly oblique ionospheric soundings, using absorption data for decametric radio waves of different polarization. It is shown that with mean square measurement errors of 0.5 dB, and using regularization algorithms to solve the inverse problems, electron collision frequency profiles can be obtained for the night F-region with errors of less than 30%. Both temperatures of electrons and neutrals are also determined to within 10%.     

High-voltage rectifier of the Energy-2 generator for electromagnetic sounding and monitoring of earthquake source zones

The high-voltage rectifier is described developed for the Energy-2 generator, with a capacity of 200 kW in which a step-up power transformer is used as a converter. The Energy-2 generator is intended for solving problems of precision deep electromagnetic monitoring of seismoactive regions of the Earth’s crust to find earthquake precursors. Theoretical investigation and numerical simulation of the high-voltage rectifier are carried out and parameters of its elemental base are optimized. All the high-voltage rectifier components, heat sink, and forced cooling system were manufactured on the basis of the developed documentation. The high-voltage rectifier was successfully tested as part of the“Energy-2 generator in August 2009 during the “FENIKS-2009” experiment.