岩矿石电磁频散规律及其在电磁波法勘探中的影响和校正

本文以大量的实测资料为依据,总结了岩矿石电磁频散的规律性。结合电磁波勘探方法中的实际问题,讨论了频散对相位常数α和吸收系数β所产生的影响及误差。提出了在实际工作中对频散进行校正的必要性,并且给出了校正公式。     

Deconvolution of teleseismic recordings for mantle structure

Delineation of detailed mantle structure frequently requires the separation of source signature and structural response from seismograms recorded at teleseismic distances. This deconvolution problem can be posed in a log-spectral domain where the operation of time-domain convolution is reduced to an additive form. The introduction of multiple events recorded at many stations leads to a system of consistency equations that must be honoured by both the source time functions and the impulse responses associated with propagation paths between sources and receivers. The system is inherently singular, and stabilization is accomplished through the supply of an initial estimate of the source time function. Although alternative choices exist, an effective estimate is derived from the eigenimage associated with the largest eigenvalue in a singular-value decomposition of the suite of aligned seismograms corresponding to a given event. The relation of the deconvolution scheme to simultaneous least-squares deconvolution is examined. Application of the methodology to broadband teleseismic P waveforms recorded on the Canadian National Seismograph Network demonstrates the retrieval of effective Green's functions including secondary phases associated with upper-mantle structure.     

Centroid moment tensor inversion of low-frequency seismic spectra using Green's functions for aspherical earth models

We present a new method for centroid moment tensor (CMT) inversion, in which we employ the Green's function computed for aspherical earth models using the Direct Solution Method. We apply this method to CMT inversion of low-frequency seismic spectra for the 1994 Bolivia and 1996 Flores Sea deep earthquakes. The estimated centroid locations agree well with those obtained by multiple-shock analyses using body-wave data. This shows that it is possible to obtain reliable CMT solutions by analyses of low-frequency seismic spectra using accurate Green's functions computed for present 3-D earth models.     

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.     

Source investigation of a small event using empirical Green's functions and simulated annealing

We propose a two-step inversion of three-component seismograms that (1) recovers the far-field source time function at each station and (2) estimates the distribution of co-seismic slip on the fault plane for small earthquakes (magnitude 3 to 4). The empirical Green's function (EGF) method consists of finding a small earthquake located near the one we wish to study and then performing a deconvolution to remove the path, site, and instrumental effects from the main-event signal.
The deconvolution between the two earthquakes is an unstable procedure: we have therefore developed a simulated annealing technique to recover a stable and positive source time function (STF) in the time domain at each station with an estimation of uncertainties. Given a good azimuthal coverage, we can obtain information on the directivity effect as well as on the rupture process. We propose an inversion method by simulated annealing using the STF to recover the distribution of slip on the fault plane with a constant rupture-velocity model. This method permits estimation of physical quantities on the fault plane, as well as possible identification of the real fault plane.
We apply this two-step procedure for an event of magnitude 3 recorded in the Gulf of Corinth in August 1991. A nearby event of magnitude 2 provides us with empirical Green's functions for each station. We estimate an active fault area of 0.02 to 0.15 km² and deduce a stress-drop value of 1 to 30 bar and an average slip of 0.1 to 1.6 cm. The selected fault of the main event is in good agreement with the existence of a detachment surface inferred from the tectonics of this half-graben.