遗传算法在上地幔速度结构研究中的应用

遗传算法是近年来发展较快的一种求解非线性优化问题的有效方法。本文通过对遗传算法基本原理的介绍，对该方法的特点进行了分析。我们采用ＷＫＢＪ理论地震图作正演，遗传算法作反演，对用体波波形反演上地幔速度结构的方法进行了研究。在ＷＫＢＪ理论地震图的计算中，通过计算主要射线的平均吸收特征时间，考虑了衰减随距离的变化，通过对不同震源引入不同的虚拟界面，同时对多个地震的波形记录进行反演。探讨了观测误差对反演结果的影响。对理论记录的反演表明，用遗传算法研究上地幔速度结构具有较好的效果。     

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

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

用于多频电磁波测井的负阻抗变换式有源接收天线

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

The glacilacustrine sedimentary environment of Bowser Lake in the northern Coast Mountains of British Columbia, Canada

Bowser Lake, a fiord lake in the northern Coast Mountains of British Columbia, contains a thick Holocene fill consisting mainly of silt and clay varves. These sediments were carried into the lake by proglacial Bowser River which drains a high-energy, heavily glacierized basin. Sedimentation in the lake is controlled by seasonal snow and ice melt, by autumn rainstorms, and by rare, but very large jökulhlaups from glacier-dammed lakes in the upper Bowser River basin which complicate environmental inferences from the sedimentary record. Sediment is dispersed through the deep western part of the lake by energetic turbidity currents. The turbidity currents apparently do not overtop a sill that separates the western basin from much shallower areas to the east. Large amounts of silt and clay are deposited from suspension in the eastern part of the lake, but sediment accumulation rates there are much lower than to the west. Several strong acoustic reflectors punctuate the varved fill in the western basin; these may be thick or relatively coarse beds deposited during jökulhlaups or exceptionally large storms. The contemporary sediment yield to Bowser Lake, estimated from sediments in the lake, is about 360 t km^-2a^-1. This is a relatively high value, but it is less than yields insome other, similar montane basins with extensive snow and ice cover.The most likely explanation for the difference is that large amounts of sediment have been, and continue to be, stored on the Bowser delta andin small proglacial lakes.     

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.