共查询到20条相似文献，搜索用时 296 毫秒
1.
2.
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. 相似文献
3.
4.
5.
T. W. Dawson 《Geophysical Journal International》1996,124(2):577-590
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. 相似文献
6.
A conservative staggered-grid finite difference method is presented for computing the electromagnetic induction response of an arbitrary heterogeneous conducting sphere by external current excitation. This method is appropriate as the forward solution for the problem of determining the electrical conductivity of the Earth's deep interior. This solution in spherical geometry is derived from that originally presented by Mackie et al. (1994 ) for Cartesian geometry. The difference equations that we solve are second order in the magnetic field H , and are derived from the integral form of Maxwell's equations on a staggered grid in spherical coordinates. The resulting matrix system of equations is sparse, symmetric, real everywhere except along the diagonal and ill-conditioned. The system is solved using the minimum residual conjugate gradient method with preconditioning by incomplete Cholesky decomposition of the diagonal sub-blocks of the coefficient matrix. In order to ensure there is zero H divergence in the solution, corrections are made to the H field every few iterations. In order to validate the code, we compare our results against an integral equation solution for an azimuthally symmetric, buried thin spherical shell model ( Kuvshinov & Pankratov 1994 ), and against a quasi-analytic solution for an azimuthally asymmetric configuration of eccentrically nested spheres ( Martinec 1998 ). 相似文献
7.
Pierre Langlois Michel C. Chouteau & Léonard Bolduc 《Geophysical Journal International》2000,140(3):509-520
Continuous magnetotelluric measurements were made over a period of 600 days, with 100-m-, 30-km- and 100-km-long dipoles and a period range of 40–4000 s. Data analysis for different dipole lengths indicates the presence of static shift at various scales. It is shown that the longer the telluric dipole, the less statically shifted the resistivity curves; nevertheless, static shifts can still be present due to geological structures causing anomalies exhibiting wavelengths comparable to the dipole length. Also, a relationship is observed between the coherence and the main impedance components. This relation is explained in terms of signal-to-noise ratio. A way to reduce the bias on the impedance estimates is suggested. The apparent resistivities and phases computed from three different impedance estimates using 100-km-long dipoles are then compared to those observed in similar studies made near our observation region. 相似文献
8.
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 相似文献
9.
D. Livelybrooks 《Geophysical Journal International》1993,114(3):443-458
10.
E. Ngakosso A. Straub M. Saillard P. Vincent 《Geophysical Journal International》1998,135(3):1028-1044
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. 相似文献
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. 相似文献
11.
12.
13.
14.
Recovering magnetic susceptibility from electromagnetic data over a one-dimensional earth 总被引：1，自引：0，他引：1
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. 相似文献
15.
16.
17.
18.
The conductivity structure of the Earth's mantle was estimated using the induction method down to 2100 km depth for the Europe–Asia region. For this purpose, the responses obtained at seven geomagnetic observatories (IRT, KIV, MOS, NVS, HLP, WIT and NGK) were analysed, together with reliable published results for 11 yr variations. 1-D spherical modelling has shown that, beneath the mid-mantle conductive layer (600–800 km), the conductivity increases slowly from about 1 S m−1 at 1000 km depth to 10 S m−1 at 1900 km, while further down (1900–2100 km) this increase is faster. Published models of the lower mantle conductivity obtained using the secular, 30–60 yr variations were also considered, in order to estimate the conductivity at depths down to the core. The new regional model of the lower mantle conductivity does not contradict most modern geoelectrical sounding results. This model supports the idea that the mantle base, situated below 2100 km depth, has a very high conductivity. 相似文献
19.