首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 443 毫秒
1.
2.
3.
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.  相似文献   

4.
An effective and accurate technique for the numerical solution of 2-D electromagnetic scattering problems with 3-D sources is presented. This solution introduces a set of the usual boundary integral equations and uses a scalar Green's function. In this scalar version, the unknowns of the problem are the boundary values of the longitudinal fields and their normal derivatives in the Fourier domain. A generalization of the usual boundary integral formulation enables us to handle a large class of models composed of piecewise homogeneous domains, including contiguous domains, multiply-connected domains and unbounded domains. This formulation involves the solution of a system of linear equations, and results in a significant saving in computation time in comparison with other rigorous methods.
  The requirements for the numerical implementation of this solution are described in detail. Numerical tests were carried out using the important example of electromagnetic tomography. The specific symmetry properties of the response function in this case are illustrated. Numerical accuracy is verified over a large frequency range, up to 1  MHz.  相似文献   

5.
6.
7.
8.
9.
10.
11.
A 3-D electromagnetic inversion scheme has been successfully applied to crosswell electromagnetic data collected at the Richmond Field Station near Berkeley, California. By comparing images of data collected before and after the injection of 50000 gallons of salt water, a 3-D image of the plume has been developed, which shows the location of zones of maximum permeability surrounding the injection well through which the salt water has migrated. A resolution analysis has determined that the location of the plume is fairly accurate. However, the image of the geology will be distorted due to incomplete data coverage. This latter problem is further complicated by the fact that data residuals indicate that the wells deviate from the vertical. These deviations have been shown to cause artefacts within the images and thus further reduce the accuracy of the images with respect to the geology.  相似文献   

12.
The standard 1-D inversion approach for the interpretation of transient electromagnetic (TEM) data usually fails in the presence of near-surface conductivity anomalies. Since multidimensional inversion codes are not routinely available, the only alternative to discarding the data may be trial-and-error forward modelling. We interpret data from a long-offset transient electromagnetic (LOTEM) survey which was carried out in 1995 in the Odenwald area, using 2-D finite-difference modelling. We focus on a subsegment of the LOTEM profile, which was shot with two different electric dipole transmitters. A model is found which consistently explains the electric and magnetic field data at eight locations for both transmitters. First, we introduce a conductive dyke under the receiver spread to explain sign reversals in the magnetic field transients. A conductive slab under one of the transmitters is required to obtain a reasonable quantitative fit for that transmitter. Consideration of the electric field data then requires a modification of the layered earth background. Finally, we study the response of a crustal conductor, which was the original target of the survey. The data are sensitive to the conductor, and for the investigated subset of the data the fits are slightly better without the conductive layer.  相似文献   

13.
A theoretical approach to the propagation of interacting cracks   总被引:2,自引:0,他引:2  
We propose a scheme to compute interaction effects between two randomly oriented cracks under compressive stresses and we discuss the role crack interactions play in the crack coalescence process. Stress intensity factors are computed by using an iterative technique based on the method of successive approximations. Once crack propagation occurs, curved wing cracks grow from the initial crack tips. The stress intensity factors at the wing crack tips are calculated as the sum of two terms: a component for a single wing crack subjected to both the applied stresses and the interaction effect, and a component due to the sliding of the initial crack. We have applied our procedure to various crack geometries. Our results show that interaction effects act on the crack propagation path. For cracks under tension, our approach correctly predicts the curving, hook-shaped paths of interacting cracks that have been observed in various materials. For en echelon compressive cracks, interaction effects depend on the geometry of stepping. For right-stepping cracks, no mode I crack coalescence occurs. A mixedmode propagation criterion may be introduced to check whether coalescing secondary shear fractures initiate. For left-stepping cracks, depending on whether or not there is overlapping, crack coalescence is achieved by tension wing cracks at the inner crack tips. Without overlapping, the growing wing cracks delimit a region where a tensile secondary fracture may develop and lead to coalescence. These results are consistent with previous work and show that our procedure may be now extended to a population of cracks.  相似文献   

14.
An iterative solution to the non-linear 3-D electromagnetic inverse problem is obtained by successive linearized model updates using the method of conjugate gradients. Full wave equation modelling for controlled sources is employed to compute model sensitivities and predicted data in the frequency domain with an efficient 3-D finite-difference algorithm. Necessity dictates that the inverse be underdetermined, since realistic reconstructions require the solution for tens of thousands of parameters. In addition, large-scale 3-D forward modelling is required and this can easily involve the solution of over several million electric field unknowns per solve. A massively parallel computing platform has therefore been utilized to obtain reasonable execution times, and results are given for the 1840-node Intel Paragon. The solution is demonstrated with a synthetic example with added Gaussian noise, where the data were produced from an integral equation forward-modelling code, and is different from the finite difference code embedded in the inversion algorithm  相似文献   

15.
Telluric distortion occurs when electric charges accumulate along near-surface inhomogeneities. At low frequencies, the electric currents associated with these charges can be neglected compared to currents induced deeper in the Earth. At higher frequencies, the magnetic fields associated with these currents may be significant. Some parameters describing the distortion magnetic fields can be estimated from measured magneto-telluric impedance matrices. For regional magnetic fields aligned with regional strike directions, parameters associated with the distortion magnetic field component parallel to the regional magnetic field are undeterminable, whereas parameters associated with the distortion magnetic field component perpendicular to the regional magnetic field can be estimated. Optimal estimates are straightforward even for the realistic case of measurement errors that are correlated between elements of a measured impedance matrix. In a simple example of a 1-D anisotropic model with anisotropy direction varying with depth, the modelling of distortion magnetic fields results in regional impedance estimates corresponding more closely to the responses of uncoupled isotropic models, allowing sensible interpretation of an additional one and a half decades of data.  相似文献   

16.
We investigate the reconstruction of a conductive target using crosswell time-domain electromagnetic tomography in the diffusive limit. The work is a natural extension of our ongoing research in the modification of time-domain methods for the rugged marine mid-ocean-ridge environment, an environment characterized by extreme topography and pronounced variations in crustal conductivity on all scales. We have proved both in theory and in practice that 'traveltime', the time taken for an electromagnetic signal to be identified at a receiver following a change of current in the transmitter, is an excellent, robust estimator of average conductivity on a path between transmitter and receiver. A simple estimate of the traveltime for a parallel electric dipole-dipole system is the time at which the derivative of the electric field with respect to logarithmic time at the receiver reaches its maximum. We have derived the fundamental relationship between the traveltime and the conductivity of the medium for a uniform whole-space. We have applied the concept of the traveltime inversion to the related crosswell problem and demonstrated reconstructions of finite targets based on tomographic analyses. Results show that the crosswell time-domain electromagnetic tomography can supply useful information, such as the location and shape of a conductive target.  相似文献   

17.
18.
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.  相似文献   

19.
The mode-matching method is used to obtain an exact analytical solution to the problem of B -polarization induction in two adjacent thin half-sheets, lying on a conducting layer that is terminated by a perfect conductor at finite depth. These components of the model represent, respectively, the Earth's conducting surface layers, crust, and mantle. In dimensionless variables, the model has three independent parameters, these being the two thin-sheet conductances and the layer thickness. The mode-matching solution obtained in this paper is shown to be identical lo that derived via the Wiener-Hopf method in a companion paper (Dawson 1996), and so provides additional verification of that solution. As was shown in the companion paper, the solution for the present model contains, as special limiting cases, those for three models considered earlier by various authors. The second part of the present paper addresses the solutions for the electric fields in the non-conducting half-space above the conductors, which represents the atmosphere. In the final part, sample numerical calculations are presented to illustrate the solution.  相似文献   

20.
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号