共查询到20条相似文献，搜索用时 31 毫秒
1.
We have formulated a 3-D inverse solution for the magnetotelluric (MT) problem using the non-linear conjugate gradient method. Finite difference methods are used to compute predicted data efficiently and objective functional gradients. Only six forward modelling applications per frequency are typically required to produce the model update at each iteration. This efficiency is achieved by incorporating a simple line search procedure that calls for a sufficient reduction in the objective functional, instead of an exact determination of its minimum along a given descent direction. Additional efficiencies in the scheme are sought by incorporating preconditioning to accelerate solution convergence. Even with these efficiencies, the solution's realism and complexity are still limited by the speed and memory of serial processors. To overcome this barrier, the scheme has been implemented on a parallel computing platform where tens to thousands of processors operate on the problem simultaneously. The inversion scheme is tested by inverting data produced with a forward modelling code algorithmically different from that employed in the inversion algorithm. This check provides independent verification of the scheme since the two forward modelling algorithms are prone to different types of numerical error. 相似文献
2.
3.
The electrical conductivity of the oceanic upper mantle 总被引：2，自引：0，他引：2
4.
Finite-difference modelling of magnetotelluric fields in two-dimensional anisotropic media 总被引：4，自引：0，他引：4
An algorithm for the numerical modelling of magnetotelluric fields in 2-D generally anisotropic block structures is presented. Electrical properties of the individual homogeneous blocks are described by an arbitrary symmetric and positive-definite conductivity tensor. The problem leads to a coupled system of partial differential equations for the strike-parallel components of the electromagnetic field. E x , and H x These equations are numerically approximated by the finite-difference (FD) method, making use of the integro-interpolation approach. As the magnetic component H x , is constant in the non-conductive air, only equations for the electric mode are approximated within the air layer. The system of linear difference equations, resulting from the FD approximation, can be arranged in such a way that its matrix is symmetric and band-limited, and can be solved, for not too large models, by Gaussian elimination. The algorithm is applied to model situations which demonstrate some non-trivial phenomena caused by electrical anisotropy. In particular, the effect of 2-D anisotropy on the relation between magnetotelluric impedances and induction arrows is studied in detail. 相似文献
5.
b
The effect of anisotropy on the distribution of Earth's conductivity is evaluated by calculating the electromagnetic response of multilayered 2-D structures. The electric and magnetic fields are expanded in terms of Fourier series, the coefficients being obtained by applying the corresponding boundary conditions on each interface, given by arbitrary analytical functions. Then the results are used to analyse some particular structures. 相似文献
The effect of anisotropy on the distribution of Earth's conductivity is evaluated by calculating the electromagnetic response of multilayered 2-D structures. The electric and magnetic fields are expanded in terms of Fourier series, the coefficients being obtained by applying the corresponding boundary conditions on each interface, given by arbitrary analytical functions. Then the results are used to analyse some particular structures. 相似文献
6.
7.
8.
9.
An ellipticity criterion in magnetotelluric tensor analysis 总被引：2，自引：1，他引：2
10.
Summary. In recent years telluric sounding has been replaced by MT (magnetotellurics). However, several new purely telluric parameters, besides the traditional Jacobian, have been shown to be efficient geophysical indicators of lateral conductivity variations. A set of typical two-dimensional structures is analysed to demonstrate the resolving power of the new indicators. For such telluric studies, a wide frequency band is a great asset, especially because the indicators are best displayed as pseudo-sections in the frequency domain. However, a wide frequency range is easily achieved when only the telluric field needs to be measured. In MT the magnetic sensing coils often severely reduce the available bandwidth. 相似文献
11.
12.
13.
J. P. Narayan 《Geophysical Journal International》2001,146(2):269-281
14.
15.
16.
17.
18.
19.
Transfer functions: smooth robust estimates by least-squares and remote reference methods 总被引：2，自引：0，他引：2
J. C. Larsen 《Geophysical Journal International》1989,99(3):645-663
20.
Generalized Born scattering of elastic waves in 3-D media 总被引：1，自引：0，他引：1
It is well known that when a seismic wave propagates through an elastic medium with gradients in the parameters which describe it (e.g. slowness and density), energy is scattered from the incident wave generating low-frequency partial reflections. Many approximate solutions to the wave equation, e.g. geometrical ray theory (GRT), Maslov theory and Gaussian beams, do not model these signals. The problem of describing partial reflections in 1-D media has been extensively studied in the seismic literature and considerable progress has been made using iterative techniques based on WKBJ, Airy or Langer type ansätze. In this paper we derive a first-order scattering formalism to describe partial reflections in 3-D media. The correction term describing the scattered energy is developed as a volume integral over terms dependent upon the first spatial derivatives (gradients) of the parameters describing the medium and the solution. The relationship we derive could, in principle, be used as the basis for an iterative scheme but the computational expense, particularly for elastic media, will usually prohibit this approach. The result we obtain is closely related to the usual Born approximation, but differs in that the scattering term is not derived from a perturbation to a background model, but rather from the error in an approximate Green's function. We examine analytically the relationship between the results produced by the new formalism and the usual Born approximation for a medium which has no long-wavelength heterogeneities. We show that in such a case the two methods agree approximately as expected, but that in a media with heterogeneities of all wavelengths the new gradient scattering formalism is superior. We establish analytically the connection between the formalism developed here and the iterative approach based on the WKBJ solution which has been used previously in 1-D media. Numerical examples are shown to illustrate the examples discussed. 相似文献