首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
Summary. An analytical solution is obtained for the E-polarization problem of electromagnetic induction in two adjacent half-sheets underlain by a uniform conducting half-space. In this mode the inducing magnetic field is assumed horizontal, uniform and perpendicular to the discontinuity. The same model was previously solved under B-polarization by Dawson & Weaver. The present solution then completes the study of two-dimensional induction in the described model. Further, it extends both the analytic E-polarization solution of Weidelt by the inclusion of an underlying conductor and that of Raval, Weaver & Dawson by the inclusion of arbitrary conductance values for the two surface sheets. The solution may be used as an idealized model of the coast effect and allows detailed study of the field behaviour near the discontinuity. The horizontal magnetic field on each side of the surface layer has a finite jump discontinuity at the interface and the vertical magnetic field exhibits a logarithmic singularity there. If the right-hand conductance (say) becomes infinite, the horizontal magnetic field exhibits an algebraic singularity as the coastline is approached from the right, while the vertical magnetic field does likewise from the left. Calculations are presented for the same two models as discussed in B-polarization by Dawson & Weaver and the results are compared to values obtained from a more general numerical scheme. The electric current distribution inside the conducting half-space is depicted for the second model.  相似文献   

2.
Summary. An existing 2-D integral equation method for modelling electromagnetic induction in a thin sheet at the surface of a uniform half-space can be generalized to deal with a layered half-space by the inclusion of an extra term in the integral equation. The results obtained are shown to be in excellent agreement with finite difference solutions to the same modelling problem.  相似文献   

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

4.
The relation between the seafloor electric field and the surface magnetic field is studied. It is assumed that the fields are created by a 2-D ionospheric current distribution resulting in the E-polarization. The layered earth below the sea water is characterized by a surface impedance. The electric field at the seafloor can be expressed either as an inverse Fourier transform integral over the wavenumber or as a spatial convolution integral. In both integrals the surface magnetic field is multiplied by a function that depends on the depth and conductivity of the sea water and on the properties of the basement. The fact that surface magnetic data are usually available on land, not at the sea surface, is also considered. Test computations demonstrate that the numerical inaccuracies involved in the convolution method are negligible. The theoretical equations are applied to calculate the seafloor electric fields due to an ionospheric line current or associated with real magnetic data collected by the IMAGE magnetometer array in northern Europe. Two different sea depths are considered: 100 m (the continental shelf) and 5 km (the deep ocean). It is seen that the dependence of the electric field on the oscillation period is weaker in the 5 km case than for 100 m.  相似文献   

5.
We present a spectral-finite-element approach to the 2-D forward problem for electromagnetic induction in a spherical earth. It represents an alternative to a variety of numerical methods for 2-D global electromagnetic modelling introduced recently (e.g. the perturbation expansion approach, the finite difference scheme). It may be used to estimate the effect of a possible axisymmetric structure of electrical conductivity of the mantle on surface observations, or it may serve as a tool for testing methods and codes for 3-D global electromagnetic modelling. The ultimate goal of these electromagnetic studies is to learn about the Earth's 3-D electrical structure.
Since the spectral-finite-element approach comes from the variational formulation, we formulate the 2-D electromagnetic induction problem in a variational sense. The boundary data used in this formulation consist of the horizontal components of the total magnetic intensity measured on the Earth's surface. In this the variational approach differs from other methods, which usually use spherical harmonic coefficients of external magnetic sources as input data. We verify the assumptions of the Lax-Milgram theorem and show that the variational solution exists and is unique. The spectral-finite-element approach then means that the problem is parametrized by spherical harmonics in the angular direction, whereas finite elements span the radial direction. The solution is searched for by the Galerkin method, which leads to the solving of a system of linear algebraic equations. The method and code have been tested for Everett & Schultz's (1995) model of two eccentrically nested spheres, and good agreement has been obtained.  相似文献   

6.
The construction of effective methods for electromagnetic modelling   总被引：1，自引：0，他引：1  
Summary. This paper deals with the further development of finite-difference methods for electromagnetic field modelling in two-and three-dimensional cases. The main feature of the approach suggested here is the application of generalized asymptotic boundary conditions valid with the accuracy (1/ρN), where ρ is the distance from the heterogeneities. The finite-difference approximation of problems under solution is made using the balance method, which results in 5-point difference schemes in the 2-D case and 7-point difference schemes in the 3-D case. To solve the linear system of difference equations the successive over-relaxation (SOR) method is used, the relaxation factor being chosen during the iteration procedure. In view of the vectorial character of the problem for the 3-D case, a successive blocked over-relaxation method (SBOR) is applied.
The model's validity is based on the comparison of the fields accounted at the ground surface with those computed by the integral transformation of excessive currents, determined in the heterogeneity region using the finite-difference scheme.  相似文献   

7.
Summary. A conducting slab of finite thickness divided into three segments of different conductivities and overlying a perfect conductor is proposed as a suitable two-dimensional 'control' model for testing the accuracy of the various numerical modelling programs that are available for calculating the fields induced in the Earth by an external, time-varying magnetic source. An analytic solution is obtained for this control model for the case of the magnetic field everywhere parallel to the conductivity boundaries ( B -polarization). Values of the field given by this solution for a particular set of model parameters are calculated at selected points on the surface and on a horizontal plane inside the conductor, and are tabulated to three figure accuracy for reference. They are used to check the accuracy of the results given by the finite difference program of Brewitt-Taylor & Weaver and the finite element program of Kisak & Silvester for the same model. Improved formulae for calculating the derived electric field components in B -polarization are first developed for incorporation in the finite difference program, and these give surface electric fields within 1 per cent of the analytic values, while all three field components inside the conductor are calculated to better than 96 per cent accuracy by the finite difference program. The results given by the finite element program are not quite so satisfactory. Errors somewhat greater than 10 per cent are present and although the program requires much less disk space it takes rather more CPU time to complete the calculations.  相似文献   

8.
The 2004 M = 9.2 Sumatra–Andaman earthquake profoundly altered the state of stress in a large volume surrounding the ∼1400 km long rupture. Induced mantle flow fields and coupled surface deformation are sensitive to the 3-D rheology structure. To predict the post-seismic motions from this earthquake, relaxation of a 3-D spherical viscoelastic earth model is simulated using the theory of coupled normal modes. The quasi-static deformation basis set and solution on the 3-D model is constructed using: a spherically stratified viscoelastic earth model with a linear stress–strain relation; an aspherical perturbation in viscoelastic structure; a 'static' mode basis set consisting of Earth's spheroidal and toroidal free oscillations; a "viscoelastic" mode basis set; and interaction kernels that describe the coupling among viscoelastic and static modes. Application to the 2004 Sumatra–Andaman earthquake illustrates the profound modification of the post-seismic flow field at depth by a slab structure and similarly large effects on the near-field post-seismic deformation field at Earth's surface. Comparison with post-seismic GPS observations illustrates the extent to which viscoelastic relaxation contributes to the regional post-seismic deformation.  相似文献   

9.
Magnetotelluric data from the backarc of the Central Andes in NW Argentinawere re-examined by employing impedance tensor decomposition and 2-D inversion and modelling techniques. The data in the period range of 50–15 000 s were collected on a profile of 220 km length reaching from the Eastern Cordillera across the Santa Barbara System to the Andean foreland of the Argentinean Chaco.
After a dimensionality analysis, data from most sites were treated as regional 2-D. The exception was the eastern section of the profile, where the magnetotelluric transfer functions for periods ≤ 1000 s reflect a 3-D earth. Application of two tensor decomposition schemes yielded a regional strike direction of N–S, which is the azimuth of the Central Andean mountain chains. Several 2-D models were obtained by pseudo- and full 2-D Occam inversion schemes. Special emphasis was placed on the inversion of phase data to reduce the influence of static shifts in the apparent resistivity data. The smooth inversion models all show a good conductor at depth. A final model was then calculated using a finite element forward algorithm.
The most prominent feature of the resulting model is a conductor which rises from depths of 180 km below the Chaco region to 80 km beneath the Santa Barbara System and the Eastern Cordillera. Its interpretation as a rise of the electrical asthenosphere is supported by seismic attenuation studies. Magnetotelluric results, surface heat-flow distribution in the area, and the electrical properties of crustal and mantle rocks suggest that the upper mantle is predominantly ductile beneath the Eastern Cordillera and the western Santa Barbara System. This generally agrees with anelastic seismic attenuation models of the area and is useful in discriminating between models of Q quality factor distribution.  相似文献   

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

11.
Rapid relaxation inversion of CSAMT data   总被引：7，自引：0，他引：7  
In this paper an inversion algorithm for controlled-source audio frequency magnetotelluric data is presented. This algorithm combines 2.5-D finite element forward modelling with the concepts of rapid relaxation inversion of magnetotelluric data. The inversion uses the same technique to compute sensitivities as the rapid relaxation inversion, and these approximate sensitivities are validated by comparison with exact 2.5-D sensitivities. The comparison shows that the approximate sensitivities are similar in shape to the exact sensitivities when transmitter–receiver offsets are greater than one skin depth in the Earth. The magnitudes of the two sensitivities differ but the variations with depth are similar. Tests of the algorithm on synthetic data and field data provide promising results.  相似文献   

12.
This paper presents a geometrically based algorithm for computing synthetic seismograms for energy transmitted through a 3-D velocity distribution. 3-D ray tracing is performed to compute the traveltimes and geometrical spreading (amplitude). The formulations of both kinematic and dynamic ray-tracing systems are presented. The two-point ray-tracing problem is solved by systematically updating the initial conditions and adjusting the ray direction until the ray intersects the specified endpoint. The amount of adjustment required depends on the derivatives of the position with respect to the given starting angles between consecutive rays. The algorithm uses derivatives to define the steepest-descent direction and to update the initial directions. The convergence rate depends on the complexity of the model.
Test seismograms compare favourably with those from a 2-D asymptotic ray theory algorithm and a 3-D Gaussian-beam algorithm. The algorithm is flexible in modelling arbitrary source and recorder geometries for various smoothly varying 3-D velocity distributions. The algorithm is further tested by simulating surface-to-tunnel vibroseis field data. Shear waves as well as compressional waves may be approximately included. Application of the algorithm to a data set from the Rainier Mesa of the Nevada Test Site produced a good fit to the transmitted (first arrival) traveltimes and amplitudes, with approximately 15 per cent variation in the local 3-D velocity.  相似文献   

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

14.
Summary. The forward solution of the general two-dimensional problem of induction in a model earth comprising a uniformly conducting half-space covered by a thin sheet of variable integrated conductivity is obtained. Unlike some previous treatments of similar problems, the method presented here does not require the field to be separated into its normal and anomalous parts. Both the E - and B -polarization modes of induction are considered and in each case the solution is expressed in terms of the horizontal component of the electric field satisfying, on the surface of the conductor, a singular integral equation whose kernel is a well-known analytic function. A recently published solution of the coast effect is included as a special case. The numerical procedure for solving the integral equations is described and some illustrative calculations are presented.  相似文献   

15.
New techniques for improving both the computational and imaging performance of the three-dimensional (3-D) electromagnetic inverse problem are presented. A non-linear conjugate gradient algorithm is the framework of the inversion scheme. Full wave equation modelling for controlled sources is utilized for data simulation along with an efficient gradient computation approach for the model update. Improving the modelling efficiency of the 3-D finite difference (FD) method involves the separation of the potentially large modelling mesh, defining the set of model parameters, from the computational FD meshes used for field simulation. Grid spacings and thus overall grid sizes can be reduced and optimized according to source frequencies and source–receiver offsets of a given input data set. Further computational efficiency is obtained by combining different levels of parallelization. While the parallel scheme allows for an arbitrarily large number of parallel tasks, the relative amount of message passing is kept constant. Image enhancement is achieved by model parameter transformation functions, which enforce bounded conductivity parameters and thus prevent parameter overshoots. Further, a remedy for treating distorted data within the inversion process is presented. Data distortions simulated here include positioning errors and a highly conductive overburden, hiding the desired target signal. The methods are demonstrated using both synthetic and field data.  相似文献   

16.
Summary. In an earlier work, mathematical formulation on computing the electromagnetic response of an arbitrarily shaped three-dimensional inhomogeneity in a layered earth had been worked out using an integral equation technique. The method has been used to show its efficacy by computing numerical results. Introducing suitable changes of variables the secondary contributions to Green's dyadic are put in the form of convolution integrals and are computed using a digital linear filtering scheme. The matrix equation is solved for the unknown electric fields in the inhomogeneity. The scattered fields are then calculated at the surface of the Earth using the appropriate Green's dyadic. The performance of the computations has been shown by comparing the numerical results with those obtained by analogue modelling as well as by other numerical schemes. The use of digital linear filtering saves an enormous amount of computer time.
The effects of varying excitation-frequency, conductivity of the host medium and that of the overburden have been studied in detail for a horizontal loop system traversing over a two-layered earth with a prismatic inhomogeneity situated in the lower conducting half space.  相似文献   

17.
A seismic re fraction/wide-angle reflection profile is analysed for the presence of correlated events ('phases'). The correlation problem is formulated in terms of temporally, spatially and frequency-local complex covariances. For robustness, the method concentrates on phase rather than amplitude information. This allows a computationally efficient algorithm that can make allowance for signal correlation length and can model curved wavefronts. A statistical test based on residual phase misfit across the analysed subarray is used to assess the probability that a detected event represents a real correlated signal.
With our chosen analysis parameters and confidence level (over 99.9 per cent). 1222 events were detected in the data. Using simple techniques based on 1-D earth models, detected events are associated with a small number of particular wave types. In this way, we have succeeded in classifying almost 95 per cent of the detected events. Those that remain describe those components of the data that are inconsistent with our simple ray paths in the 1-D assumption and with our prescribed tolerance. These include reverberations, near-surface guided waves and reflected waves from strongly laterally inhomogeneous structures. According to our modelling, about 25 per cent of the detected events are consistent with simple P -wave reflected energy, and these are to a very large extent (over 85 per cent) distinct from all the other wave-type models we have used. A direct mapping of the detected events into the offset-depth domain reveals dear internal and external consistencies among the detections for the various wave types. Estimated earth structure is consistent with models from previous analyses based on much larger data sets.
We have thus succeeded in extracting correlated events from the data and decomposing these, approximately but meaningfully, into distinct classes (ray paths)  相似文献   

18.
We design a numerical algorithm for wave simulation in a borehole due to multipole sources. The stress–strain relation of the formation is based on the Kelvin–Voigt mechanical model to describe the attenuation. The modelling, which requires two anelastic parameters and twice the spatial derivatives of the lossless case, simulates 3-D waves in an axisymmetric medium by using the Fourier and Chebyshev methods to compute the spatial derivatives along the vertical and horizontal directions, respectively. Instabilities of the Chebyshev differential operator due to the implementation of the fluid–solid boundary conditions are solved with a characteristic approach, where the characteristic variables are evaluated at the source central frequency. The algorithm uses two meshes to model the fluid and the solid. The presence of the logging tool is modelled by imposing rigid boundary conditions at the inner surface of the fluid mesh. Examples illustrating the propagation of waves are presented, namely, by using monopoles, dipoles and a quadrupoles as sources in hard and soft formations. Moreover, the presence of casing and layers is considered. The modelling correctly simulates the features—traveltime and attenuation—of the wave modes observed in sonic logs, namely, the P and S body waves, the Stoneley wave, and the dispersive S waves in the case of multipole sources.  相似文献   

19.
The perfectly matched layer (PML) absorbing boundary condition is incorporated into an irregular-grid elastic-wave modelling scheme, thus resulting in an irregular-grid PML method. We develop the irregular-grid PML method using the local coordinate system based PML splitting equations and integral formulation of the PML equations. The irregular-grid PML method is implemented under a discretization of triangular grid cells, which has the ability to absorb incident waves in arbitrary directions. This allows the PML absorbing layer to be imposed along arbitrary geometrical boundaries. As a result, the computational domain can be constructed with smaller nodes, for instance, to represent the 2-D half-space by a semi-circle rather than a rectangle. By using a smooth artificial boundary, the irregular-grid PML method can also avoid the special treatments to the corners, which lead to complex computer implementations in the conventional PML method. We implement the irregular-grid PML method in both 2-D elastic isotropic and anisotropic media. The numerical simulations of a VTI lamb's problem, wave propagation in an isotropic elastic medium with curved surface and in a TTI medium demonstrate the good behaviour of the irregular-grid PML method.  相似文献   

20.
Calculation of electromagnetic sensitivities in the time domain   总被引：1，自引：0，他引：1  
The speed of calculating sensitivities for 3-D conductivity structures for time- domain electromagnetic methods is significantly improved by applying the reciprocity theorem directly in the time domain. The sensitivities are obtained by convolving the electric field in the subsurface due to a transmitter at the surface with the electric field impulse response due to another transmitter, which replaces the original receiver. The acceleration compared to the classical perturbation method is approximately P/R , where P is the number of model parameters and R is the number of receiver positions. If the sensitivity has to be calculated very close to the receiver, approximate sensitivities can be obtained using an integral condition. Comparisons with the classical perturbation approach show that the method gives accurate results. Examples using transmitter–receiver configurations from a long-offset transient electromagnetics survey demonstrate the usefulness of sensitivities for the evaluation of resolution properties.  相似文献   

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

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