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

3.
An inversion method is presented for the reconstruction of interface geometry between two or more crustal layers from teleseismic traveltime residuals. The method is applied to 2-D models consisting of continuous interfaces separating constant-velocity layers. The forward problem of determining ray paths and traveltimes between incident wave fronts below the structure and receivers located on the Earth's surface is solved by an efficient and robust shooting method. A conjugate gradient method is employed to solve the inverse problem of minimizing a least-squares type objective function based on the difference between observed and calculated traveltimes. Teleseismic data do not accurately constrain average vertical structure, so a priori information in the form of layer velocities and average layer thicknesses is required. Synthetic tests show that the method can be used to reconstruct interface geometry accurately, even in the presence of data noise. Tests also show that, if layer velocities and initial interface positions are poorly chosen, lateral structure is still recoverable. The inversion method was applied to previously published teleseismic data recorded by an in-line array of portable seismographs that traversed the northern margin of the Musgrave Block, central Australia. The solution based on interface parametrization is consistent with models given by other studies that used the same data but different methods, most notably the standard tomographic approach that inverts for velocity rather than interface structure.  相似文献   

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

5.
Collocated magnetotelluric (MT) and seismic profiling is emerging as a necessary combined approach for deep and near-surface imaging but the resulting experimental data are typically interpreted separately since no production programs exist for multidimensional joint inversion of MT and seismic data. We present a joint 2-D inversion approach for imaging collocated MT and seismic refraction data with cross-gradient structural constraints. We describe the main features of the algorithm and first apply it to synthetic data generated for a hypothetical complex geological model. For the synthetic data, we find that the scheme leads to models with remarkable structural resemblance and improved estimates of electrical resistivity and seismic velocity. We apply the scheme to near-surface field data to test the consistency of a previously suggested resistivity–velocity interrelationship and its potential use for subsurface lithofacies discrimination or structural classification. The MT-seismic relationship is found to be in excellent accord with that derived previously for DC resistivity and seismic data set at the test site. Our results suggest that joint MT-seismic cross-gradient imaging leads to improved characterization of heterogeneous geological targets at near-surface to mantle depths.  相似文献   

6.
20 magnetotelluric (MT) soundings were collected on the Isle of Skye, Scotland to provide a high-resolution three-dimensional (3-D) electrical resistivity model of a volcanic province within the framework of a project jointly interpreting gravity, seismic, geological and MT data. The full 3-D inversion of the MT data jointly interpreted with gravity data reveals upper crustal structure. The main features of the model are interpreted in conjunction with previous geological mapping and borehole data. Our model extends to 13 km depth, several kilometres below the top of the Lewisian basement. The top of the Lewisian basement is at approximately 7–8 km depth and the topography of its surface was controlled by Precambrian rifting, during which a 4.5 km thick sequence of Torridonian sediments was deposited. The Mesozoic sediments above, which can reach up to 2.2 km thick, have small-scale depocentres and are covered by up to 600 m of Tertiary lava flows. The interpretation of the resistivity model shows that 3-D MT inversion is an appropriate tool to image sedimentary structures beneath extrusive basalt units, where conventional seismic reflection methods may fail.  相似文献   

7.
Inversion of time domain three-dimensional electromagnetic data   总被引：7，自引：0，他引：7  
We present a general formulation for inverting time domain electromagnetic data to recover a 3-D distribution of electrical conductivity. The forward problem is solved using finite volume methods in the spatial domain and an implicit method (Backward Euler) in the time domain. A modified Gauss–Newton strategy is employed to solve the inverse problem. The modifications include the use of a quasi-Newton method to generate a pre-conditioner for the perturbed system, and implementing an iterative Tikhonov approach in the solution to the inverse problem. In addition, we show how the size of the inverse problem can be reduced through a corrective source procedure. The same procedure can correct for discretization errors that inevidably arise. We also show how the inverse problem can be efficiently carried out even when the decay time for the conductor is significantly larger than the repetition time of the transmitter wave form. This requires a second processor to carry an additional forward modelling. Our inversion algorithm is general and is applicable for any electromagnetic field  ( E , H , d B / dt )  measured in the air, on the ground, or in boreholes, and from an arbitrary grounded or ungrounded source. Three synthetic examples illustrate the basic functionality of the algorithm, and a result from a field example shows applicability in a larger-scale field example.  相似文献   

8.
Topographic effects due to irregular surface terrain may prevent accurate interpretation of magnetotelluric (MT) data. Three-dimensional (3-D) topographic effects have been investigated for a trapezoidal hill model using an edge finite-element method. The 3-D topography generates significant MT anomalies, and has both galvanic and inductive effects in any polarization. This paper presents two different correction algorithms, which are applied to the impedance tensor and to both electric and magnetic fields, respectively, to reduce topographic effects on MT data. The correction procedures using a homogeneous background resistivity derived from a simple averaging method effectively decrease distortions caused by surface topography, and improve the quality of subsurface interpretation. Nonlinear least-squares inversion of topography-corrected data successfully recovers most of structures including a conductive or resistive dyke.  相似文献   

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

10.
The inversion of high-resolution geoid anomaly maps derived from satellite altimetry should allow one to retrieve the lithospheric elastic thickness, T e , and crustal density, c . Indeed, the bending of a lithospheric plate under the load of a seamount depends on both parameters, and the associated geoid anomaly is correspondingly dependent on the two parameters. The difference between the observed and modelled geoid signatures is estimated by a cost function, J , of the two variables, T e and c . We show that this cost function forms a valley structure along which many local minima appear, the global minimum of J corresponding to the true values of the lithospheric parameters. Classical gradient methods fail to find this global minimum because they converge to the first local minimum of J encountered, so that the final parameter estimate strongly depends on the starting pair of values ( T e ,   c ). We here implement a non-linear optimization algorithm to recover these two parameters from altimetry data. We demonstrate from the inversion of synthetic data that this approach ensures robust estimates of T e and c by activating two search phases alternately: a gradient phase to find a local minimum of J , and a tunnelling phase through high values of the cost function. The accuracy of the solution can be improved by a search in an iteratively restricted parameter subspace. Applying our non-linear inversion to the Great Meteor Seamount geoid data, we further show that the inverse problem is intrinsically ill-posed. As a consequence, minute geoid (or gravity) data errors can induce large changes in any recovery of lithospheric elastic thickness and crustal density.  相似文献   

11.
Summary. The Backus-Gilbert method has been extended to the estimation of the seismic wave velocity distribution in 2-D or 3-D inhomogeneous media from a finite set of travel-time data. The method may be applied to the inversion of body wave as well as surface wave data. The problem of determining a local average of the unknown velocity corrections may be reduced to a choice of a suitable δ-ness criterion for the averaging kernel. For 2-D and 3-D inhomogeneous media the simplest criterion is to minimize a sum of 'spreads' over all the coordinates. The use of this criterion requires the solution (the averaged velocity corrections) to be represented as a sum of functions, each of which depends only on one coordinate. This is a basic restriction of the method. In practice it is possible to achieve good agreement between the solution and a real velocity distribution by a reasonable choice of the coordinate system.
Numerical tests demonstrate the efficiency of the method. Some examples of the application of the method to the inversion of real seismological data for body and surface waves are given.  相似文献   

12.
Signature of remnant slabs in the North Pacific from P-wave tomography   总被引：1，自引：0，他引：1  
A 3-D ray-tracing technique was used in a global tomographic inversion in order to obtain tomographic images of the North Pacific. The data reported by the Geophysical Survey of Russia (1955–1997) were used together with the catalogues of the International Seismological Center (1964–1991) and the US Geological Survey National Earthquake Information Center (1991–1998), and the recompiled catalogue was reprocessed. The final data set, used for following the inversion, contained 523 430 summary ray paths. The whole of the Earth's mantle was parametrized by cells of 2° × 2° and 19 layers. The large and sparse system of observation equations was solved using an iterative LSQR algorithm.
A subhorizontal high-velocity anomaly is revealed just above the 660 km discontinuity beneath the Aleutian subduction zone. This high-velocity feature is observed at latitudes of up to ~70°N and is interpreted as a remnant of the subducted Kula plate, which disappeared through ridge subduction at about 48 Ma. A further positive velocity perturbation feature can be identified beneath the Chukotka peninsula and Okhotsk Sea, extending from ~300 to ~660 km depth and then either extending further down to ~800 km (Chukotka) or deflecting along the 660 km discontinuity (Okhotsk Sea). This high-velocity anomaly is interpreted as a remnant slab of the Okhotsk plate accreted to Siberia at ~55 Ma.  相似文献   

13.
A Magnetotelluric (MT) sounding was carried out at a site in south-east Queensland, in the Clarence-Moreton Basin. The synoptic recordings were taken over a period of four months at sampling frequencies from 500 Hz to 5 × 10-5 Hz. The resulting data was analysed by the stationary cross-frequency and the Cone kernel time-frequency distribution (TFD) methods of MT analysis. The results were compared as apparent resistivities on a daily basis for frequencies above 1 Hz, as well as over all the available data. The TFD MT apparent-resistivity results were more stable and less noisy on an daily basis than the cross-frequency results. Similarly the TFD analysis gave less noisy results than the cross-frequency analysis when all available data was processed. Application of these new non-stationary analysis techniques to MT processing should decrease the bias error problem of the MT methods and so increase reliability and repeatability of MT soundings.  相似文献   

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

15.
Electromagnetic investigations are usually intended to examine regional structures where induction takes place at a given period range. However, the regional information is often distorted by galvanic effects at local conductivity boundaries. Bahr (1985) and Groom & Bailey (1989) developed a physical distortion model for decomposing the MT impedance tensor, based upon local galvanic distortion of a regional 2-D electromagnetic field. We have extended their method to predict the magnetic variation fields created at an array of sites. The magnetic response functions at periods around 1000 s may be distorted by large-scale inhomogeneities in the upper or middle crust. In this period range, the data measured by a magnetometer array contain common information that can be extracted if the data set is treated as a unit, for example by using hypothetical event analysis. With this technique it is always possible to recover the regional strike direction from distorted data, even if a strong, spatially varying regional vertical field component is present in the data set. The determination of the regional impedance phases, on the other hand, is far more sensitive to deviations from the physical distortion model.
The approach has been used to investigate the Iapetus data set. For the array, which covers an area of 200  km × 300  km in northern England/southern Scotland, the technique revealed a common regional strike azimuth of ca . N125° E in the period range 500–2000  s. This direction differs from the strike indicated by the induction arrows, which seem influenced mainly by local current concentrations along the east–west-striking Northumberland Trough and a NE–SW-striking mid-crustal conductor. Both impedance phases are positive and differ by ca . 10°, which supports the assumptions of distortion fields in the data set and that the regional structure is 2-D.  相似文献   

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

17.
We develop a systematic approach to the phase identification of late-arriving groups in 2-D seismic data. Waveforms in the same traveltime branch are grouped, and synthetic traveltimes for all phases are calculated using an initial approximation to the 2-D structure. For each group, we identify the two synthetic phases providing the smallest RMS residuals. If their ratio is less than some predetermined threshold, then the group's phase is ambiguous and both assignments must be tested by traveltime inversion. If there are n unidentified groups, we construct 2 n phase tables and perform a traveltime inversion on every plausible phase assignment. The phase table that provides the highest value of the posterior probability density is taken as correct, and a 2-D velocity model is constructed from the data. This approach is shown to be effective and efficient on both simulated and real data. In addition, the residuals associated with late-arriving groups provide a means of identifying deficiencies in the initial model.  相似文献   

18.
We investigate the use of general, non- l 2 measures of data misfit and model structure in the solution of the non-linear inverse problem. Of particular interest are robust measures of data misfit, and measures of model structure which enable piecewise-constant models to be constructed. General measures can be incorporated into traditional linearized, iterative solutions to the non-linear problem through the use of an iteratively reweighted least-squares (IRLS) algorithm. We show how such an algorithm can be used to solve the linear inverse problem when general measures of misfit and structure are considered. The magnetic stripe example of Parker (1994 ) is used as an illustration. This example also emphasizes the benefits of using a robust measure of misfit when outliers are present in the data. We then show how the IRLS algorithm can be used within a linearized, iterative solution to the non-linear problem. The relevant procedure contains two iterative loops which can be combined in a number of ways. We present two possibilities. The first involves a line search to determine the most appropriate value of the trade-off parameter and the complete solution, via the IRLS algorithm, of the linearized inverse problem for each value of the trade-off parameter. In the second approach, a schedule of prescribed values for the trade-off parameter is used and the iterations required by the IRLS algorithm are combined with those for the linearized, iterative inversion procedure. These two variations are then applied to the 1-D inversion of both synthetic and field time-domain electromagnetic data.  相似文献   

19.
Surface wave tomography of the Barents Sea and surrounding regions   总被引：1，自引：0，他引：1  
The goal of this study is to refine knowledge of the structure and tectonic history of the European Arctic using the combination of all available seismological surface wave data, including historical data that were not used before for this purpose. We demonstrate how the improved data coverage leads to better depth and spatial resolution of the seismological model and discovery of intriguing features of upper-mantle structure. To improve the surface wave data set in the European Arctic, we extensively searched for broad-band data from stations in the area from the beginning of the 1970s until 2005. We were able to retrieve surface wave observations from regional data archives in Norway, Finland, Denmark and Russia in addition to data from the data centres of IRIS and GEOFON. Rayleigh and Love wave group velocity measurements between 10 and 150 s period were combined with existing data provided by the University of Colorado at Boulder. This new data set was inverted for maps showing the 2-D group-velocity distribution of Love and Rayleigh waves for specific periods. Using Monte Carlo inversion, we constructed a new 3-D shear velocity model of the crust and upper mantle beneath the European Arctic which provides higher resolution and accuracy than previous models. A new crustal model of the Barents Sea and surrounding areas, published recently by a collaboration between the University of Oslo, NORSAR and the USGS, constrains the 3-D inversion of the surface wave data in the shallow lithosphere. The new 3-D model, BARMOD, reveals substantial variations in shear wave speeds in the upper mantle across the region with a nominal resolution of 1°× 1°. Of particular note are clarified images of the mantle expression of the continent-ocean transition in the Norwegian Sea and a deep, high wave speed lithospheric root beneath the Eastern Barents Sea, which presumably is the remnant of several Palaeozoic collisions.  相似文献   

20.
The inversion of recent borehole temperatures has proved to be a successful tool to determine ancient ground surface temperature histories. To take into account heterogeneity of thermal properties and their non-linear dependence on temperature itself, a versatile 1-D inversion technique based on a finite-difference approach has been developed. Regularization of the generally ill-posed problem is obtained by an appropriate version of Tikhonov regularization of variable order. In this approach, a regularization parameter has to be determined, representing a trade-off between data fit and model smoothness. We propose to select this parameter by generalized cross-validation. The resulting technique is employed in case studies from the Kola ultradeep drilling site, and another borehole from northeastern Poland. Comparing the results from both sites corroborates the hypothesis that subglacial ground surface temperatures as met in Kola often are much higher than the ones in areas exposed to atmospheric conditions (Poland).  相似文献   

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

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