Recovering magnetic susceptibility from electromagnetic data over a one-dimensional earth

While the inversion of electromagnetic data to recover electrical conductivity has received much attention, the inversion of those data to recover magnetic susceptibility has not been fully studied. In this paper we invert frequency-domain electromagnetic (EM) data from a horizontal coplanar system to recover a 1-D distribution of magnetic susceptibility under the assumption that the electrical conductivity is known. The inversion is carried out by dividing the earth into layers of constant susceptibility and minimizing an objective function of the susceptibility subject to fitting the data. An adjoint Green's function solution is used in the calculation of sensitivities, and it is apparent that the sensitivity problem is driven by three sources. One of the sources is the scaled electric field in the layer of interest, and the other two, related to effective magnetic charges, are located at the upper and lower boundaries of the layer. These charges give rise to a frequency-independent term in the sensitivities. Because different frequencies penetrate to different depths in the earth, the EM data contain inherent information about the depth distribution of susceptibility. This contrasts with static field measurements, which can be reproduced by a surface layer of magnetization. We illustrate the effectiveness of the inversion algorithm on synthetic and field data and show also the importance of knowing the background conductivity. In practical circumstances, where there is no a priori information about conductivity distribution, a simultaneous inversion of EM data to recover both electrical conductivity and susceptibility will be required.     

Time-domain electromagnetic migration in the solution of inverse problems

Time-domain electromagnetic (TDEM) migration is based on downward extrapolation of the observed field in reverse time. In fact, the migrated EM field is the solution of the boundary-value problem for the adjoint Maxwell's equations. The important question is how this imaging technique can be related to the solution of the geoelectrical inverse problem. In this paper we introduce a new formulation of the inverse problem, based on the minimization of the residual-field energy flow through the surface or profile of observations. We demonstrate that TDEM migration can be interpreted as the first step in the solution of this specially formulated TDEM inverse problem. However, in many practical situations this first step produces a very efficient approximation to the geoelectrical model, which makes electromagnetic migration so attractive for practical applications. We demonstrate the effectiveness of this approach in inverting synthetic and practical TDEM data.     

Candidate precursors: pulse-like geoelectric signals possibly related to recent seismic activity in Japan

Telemetric network observations of pulse-like geoelectric charge signals using a vertical dipole buried under the ground were undertaken at various observation sites over a wide area in Japan from April 1996. From continuous records of the signals during the six months following that, we attempted to select anomalous signals, possibly related to seismic electric activity. Special attention was paid to the elimination of noise resulting from industrial and meteorological electric activity, comparison with other electromagnetic signals in the VLF band and the relation between the precursor period and the distance from the eventual main shock that occurred in Japan. Four candidate precursor electric signals, which were not contaminated by industrial and meteorological electric activity, were then selected, of which the second appeared before the Akita-ken Nairiku-nanbu earthquake swarm, for which the maximum M = 5.9 occurred on 1996 August 11, and the third and fourth before the Chiba-ken Toho-oki earthquake, M = 6.6, on 1996 September 11. A tentative qualitative model explaining why the candidate precursory signal is related to stress building up before an earthquake is presented in terms of the electrification of gases released from fracturing rocks immediately prior to the main shock.     

A general correspondence principle for time-domain electromagnetic wave and diffusion fields

A general correspondence principle is presented that relates any time-domain electromagnetic diffusion field to an electromagnetic wavefield in a 'corresponding' configuration. The principle applies to arbitrarily inhomogeneous and anisotropic media and arbitrary transmitters and receivers. For the correspondence between the two types of electromagnetic fields to hold, the electric conductivity in the diffusive case and the permittivity in the wavefield case should have the same spatial variation, while the permeability distributions in space in the two cases are to be identical. Essential steps in the derivation of the correspondence principle are the use of the time Laplace transformation of causal signals, taken at real, positive values of the transform parameter, the Schouten-Van der Pol theorem in the theory of the Laplace transformation, and the reliance upon Lerch's theorem of the uniqueness of the interrelation between causal field quantities and their time-Laplace-transform representations at real, positive values of the transform parameter. Correspondence is then established between the tensorial Green's functions in the two cases, where the Green's functions are the point-receiver responses (either electric or magnetic field) to point-transmitter excitations (either electric- or magnetic-current source).
Through the correspondence principle, all transient electromagnetic wavefields (where losses are neglected) have as a counterpart a transient diffusive electromagnetic field (where the electric displacement current is neglected). The interrelation yields the tool to compare quantitatively the potentialities of the two types of fields in transient electromagnetic geophysical prospecting.
Finally, a general medium-parameter scaling law for time-domain electromagnetic wavefields is presented.     

用于多频电磁波测井的圆锥螺旋式有源接收天线

对用负阻抗变换电路同圆锥螺旋天线相结合，构成的有源接收天线进行了理论分析。并给出了该有源天线稳定工作的条件和设计原理。结果表明，用该方法设计的有源接收天线具有很宽的频带、较高的灵敏度和很小的尺寸，是多频电磁波测井较为理想的接收天线。     

覆盖层下三维异常体对垂直磁偶极源的电磁响应

针对地下半空间存在三维异常体及覆盖层情况，计算了垂直磁偶极子的场，计算结果表明σ１与σ２接近时，覆盖层影响减小；σ１增加时覆盖层影响增大；对于不同的σ２／σ１值，场幅值对高频（ｆ≥２ｋＨｚ）和低频（ｆ≤１０Ｈｚ）都明显减小，只有在２ｋＨｚ≥ｆ≥１０Ｈｚ之间，场幅值可达最大。还表明垂直磁偶极子的场具有较强的岩性识别和分层能力     

碳酸盐岩地区电磁波场的特性分析

形成岩溶的高阻碳酸盐岩对电磁波的吸收较弱，而且成层沉积的特定环境使碳酸盐岩地区岩石的各向异性现象较明显，由于地空界面的存在及岩溶洞壁两侧介质电性差异明显，界面清晰，又使电磁波出现强烈的反射及散射现象，从理论上对碳酸盐岩地区电磁波传播的特有现象进行了研究和分析，并给出结论。     

瞬变电磁场的奇点展开法

在分析电磁系统自然谐振现象的基础上,简述了瞬变电磁场奇点展开法的物理基础、基本理论、极点提取方法,并考查了奇点展开法在地球物理中的应用可能性。     

双频电磁波电导率层析成象

本文论述双频透射电磁波电导率层析成象方法。设｜Ｅ１｜，｜Ｅ２｜为在远场区观测的同一电偶极天线的两个辐射电场强度。如地质介质的电导率连续变化，观测场强的辐射频率ω１，ω２之差较小，则有如下的电导率层析成象方程｜Ｅ１｜－｜Ｅ２｜＝２０ｌｏｇ（ω１／ω２）＋２０ｌｏｇ（ｅ）μ／２（ω２－ω１）∫σ～（ｒ）ｄｒ式中σ＾＝σ，σ为介质的电导率。双频透射电磁波电导率层析成象方法克服了长期困扰电磁波吸收层析成象的初始辐射场强计算问题，应用实例表明它可获得更精确的层析图。