首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 15 毫秒
1.
Joint inversion of receiver function and surface wave dispersion observations   总被引：16，自引：0，他引：16  
We implement a method to invert jointly teleseismic P wave receiver functions and surface wave group and phase velocities for a mutually consistent estimate of earth structure. Receiver functions are primarily sensitive to shear wave velocity contrasts and vertical traveltimes, and surface wave dispersion measurements are sensitive to vertical shear wave velocity averages. Their combination may bridge resolution gaps associated with each individual data set. We formulate a linearized shear velocity inversion that is solved using a damped leastsquares scheme that incorporates a priori smoothness constraints for velocities in adjacent layers. The data sets are equalized for the number of data points and physical units in the inversion process. The combination of information produces a relatively simple model with a minimal number of sharp velocity contrasts. We illustrate the approach using noisefree and realistic noise simulations and conclude with an inversion of observations from the Saudi Arabian Shield. Inversion results for station SODA, located in the Arabian Shield, include a crust with a sharp gradient near the surface (shear velocity changing from 1.8 to 3.5 km s1 in 3 km) underlain by a 5kmthick layer with a shear velocity of 3.5 km s1 and a 27kmthick layer with a shear velocity of 3.8 km s1, and an upper mantle with an average shear velocity of 4.7 km s1. The crustmantle transition has a significant gradient, with velocity values varying from 3.8 to 4.7 km s1 between 35 and 40 km depth. Our results are compatible with independent inversions for crustal structure using refraction data.  相似文献   

2.
A tomographic inversion technique that inverts traveltimes to obtain a model of the subsurface in terms of velocities and interfaces is presented. It uses a combination of refraction, wide-angle reflection and normal-incidence data, it simultaneously inverts for velocities and interface depths, and it is able to quantify the errors and trade-offs in the final model. The technique uses an iterative linearized approach to the non-linear traveltime inversion problem. The subsurface is represented as a set of layers separated by interfaces, across which the velocity may be discontinuous. Within each layer the velocity varies in two dimensions and has a continuous first derivative. Rays are traced in this medium using a technique based on ray perturbation theory, and two-point ray tracing is avoided by interpolating the traveltimes to the receivers from a roughly equidistant fan of rays. The calculated traveltimes are inverted by simultaneously minimizing the misfit between the data and calculated traveltimes, and the roughness of the model. This 'smoothing regularization' stabilizes the solution of the inverse problem. In practice, the first iterations are performed with a high level of smoothing. As the inversion proceeds, the level of smoothing is gradually reduced until the traveltime residual is at the estimated level of noise in the data. At this point, a minimum-feature solution is obtained, which should contain only those features discernible over the noise.
The technique is tested on a synthetic data set, demonstrating its accuracy and stability and also illustrating the desirability of including a large number of different ray types in an inversion.  相似文献   

3.
Summary. The computational effectiveness of travel-time inversion methods depends on the parameterization of a 3-D velocity structure. We divide a region of interest into a few layers and represent the perturbation of wave slowness in each layer by a series of Chebyshev polynomials. Then a relatively complex velocity structure can be dcscribed by a small set of parameters that can be accurately evaluated by a linearized inversion of travel-time residuals. This method has been applied to artificial and real data at small epicentral distances and in the teleseismic distance range. The corresponding matrix equations were solved using singular value decomposition. The results suggest that the method combines resolution with computational convenience.  相似文献   

4.
A self-adaptive automated parametrization approach is suggested for the sequential inversion of controlled-source seismic tomography and gravity data. The velocities and interfaces are parametrized by their Haar wavelet expansion coefficients. Only those coefficients that are well constrained by the data, as measured by the number of rays that cross the corresponding wavelet function support area and their angular coverage, are inverted for, others are set to zero. This approach results in a reasonable distribution of resolution throughout the model even in cases of irregular ray coverage and does overcome the trade-off between different types of model parameters. A modified sequential inversion approach is suggested to join the traveltimes and gravity anomalies inversion. An algorithm is developed that inverts for smooth velocity and density variations inside the seismic layer, the position of its bottom interface as well as for optimal values of the velocity-to-density regression coefficients. The algorithm makes use of direct (diving), reflected and head (critically refracted) wave traveltimes. The algorithm workflow is demonstrated on a synthetic data example.  相似文献   

5.
The presence of anisotropy requires that tomographic methods be generalized to account for anisotropy. This generalization allows geological structure to be correctly imaged and allows the anisotropic parameters to be estimated. Use of isotropic inversion for imaging anisotropic structures gives systematic trends in the traveltime and polarization residuals. However, due to the limited directional coverage, the traveltimes along may not be sufficient to study the anisotropic properties of the structure. Polarizations can provide independent information on the structure. Traveltime and polarization inversion are applied to synthetic examples simulating VSP experiments. Transverse isotropy and 1-D structure are assumed. Plots of traveltime and polarization residuals are an important tool to detect the anomalies due to the presence of anisotropy. For receivers located in anisotropic layers, polarization residuals display consistent anomalies of several degrees. The synthetic examples show that even the simple 1-D problem is difficult, when using direct arrivals only. Large a posteriori errors in anisotropic parameters are obtained by traveltime inversion in layers where available incidence angles are less than 45°. Resolution of the tomographic image of VSP data is greatly improved by a combination of traveltime and polarization information. In order to obtain accurate inversion results, the measurement error of polarization data should be kept to within a few degrees.  相似文献   

6.
Seismic traveltimes and amplitudes in reflection-seismic data show different dependences on the geometry of reflection interfaces, and on the variation of interval velocities. These dependences are revealed by eigenanalysis of the Hessian matrix, defined in terms of the Fréchet matrix and its adjoint associated with different norms chosen in the model space. The eigenvectors and eigenvalues of the Hessian clearly show that for reflection tomographic inversion, traveltime and amplitude data contain complementary information. Both for reflector-geometry and for interval-velocity variations, the traveltimes are sensitive to the model components with small wavenumbers, whereas the amplitudes are more sensitive to the components with high wavenumbers. The model resolution matrices, after the rejection of eigenvectors corresponding to small eigenvalues, give us some insight into how the addition of amplitude information could potentially contribute to the recovery of physical parameters.
In order to cooperatively invert seismic traveltimes and amplitudes simultaneously, we propose an empirical definition of the data covariance matrix which balances the relative sensitivities of different types of data. We investigate the cooperative use of both data types for, separately, interface-geometry and 2-D interval-velocity variations. In both cases we find that cooperative inversions can provide better solutions than those using traveltimes alone. The potential benefit of including amplitude-data constraints in seismic-reflection traveltime tomography is therefore that it may be possible to resolve the known ambiguity between the reflector-depth uncertainty and the interval-velocity uncertainty better.  相似文献   

7.
P-SH conversion is commonly observed in teleseismic P waves, and is often attributed to dipping interfaces beneath the receiver. Our modelling suggests an alternative explanation in terms of flat-layered anisotropy. We use reflectivity techniques to compute three-component synthetic seismograms in a 1-D anisotropic layered medium. For each layer of the medium, we prescribe values of seismic velocities and hexagonally symmetric anisotropy about a common symmetry axis of arbitrary orientation. A compressional wave in an anisotropic velocity structure suffers conversion to both SV -and SH -polarized shear waves, unless the axis of symmetry is everywhere vertical or the wave travels parallel to all symmetry axes. The P-SV conversion forms the basis of the widely used 'receiver function' technique. The P-SH conversion occurs at interfaces where one or both layers are anisotropic. A tilted axis of symmetry and a dipping interface in isotropic media produce similar amplitudes of both direct ( P ) and converted ( Ps ) phases, leaving the backazimuth variation of the P-Ps delay as the main discriminant. Seismic anisotropy with a tilted symmetry axis leads to complex synthetic seismograms in velocity models composed of just a few flat homogeneous layers. It is possible therefore to model observations of P coda with prominent transverse components with relatively simple 1-D velocity structures. Successful retrieval of salient model characteristics appears possible using multiple realizations of a genetic-algorithm (GA) inversion of P coda from several backazimuths. Using GA inversion, we determine that six P coda recorded at station ARU in central Russia are consistent with models that possess strong (> 10 per cent) anisotropy in the top 5 km and between 30 and 43 km depth. The symmetry axes are tilted, and appear aligned with the seismic anisotropy orientation in the mantle under ARU suggested by SKS splitting.  相似文献   

8.
We use teleseismic three-component digital data from the Trabzon, Turkey broadband seismic station TBZ to model the crustal structure by the receiver function method. The station is located at a structural transition from continental northeastern Anatolia to the oceanic Black Sea basin. Rocks in the region are of volcanic origin covered by young sediments. By forward modelling the radial receiver functions, we construct 1-D crustal shear velocity models that include a lower crustal low-velocity zone, indicating a partial melt mechanism which may be the source of surfacing magmatic rocks and regional volcanism. Within the top 5 km, velocities increase sharply from about 1.5 to 3.5 km s−1. Such near-surface low velocities are caused by sedimentation, extending from the Black Sea basin. Velocities at around 20 km depth have mantle-like values (about 4.25 km s−1 ), which easily correlate to magmatic rocks cropping out on the surface. At 25 km depth there is a thin low-velocity layer of about 4.0 km s−1. The average Moho velocity is about 4.6 km s−1, and its depth changes from 32 to 40 km. Arrivals on the tangential components indicate that the Moho discontinuity dips approximately southwards, in agreement with the crustal thickening to the south. We searched for the solution of receiver functions around the regional surface wave group velocity inversion results, which helped alleviate the multiple solution problem frequently encountered in receiver function modelling.
Station TBZ is a recently deployed broadband seismic station, and the aim of this study is to report on the analysis of new receiver function data. The analysis of new data in such a structurally complex region provides constraining starting models for future structural studies in the region.  相似文献   

9.
Velocity estimation remains one of the main problems when imaging the subsurface with seismic reflection data. Traveltime inversion enables us to obtain large-scale structures of the velocity field and the position of seismic reflectors. However, as the media currently under study are becoming more and more complex, we need to know the finer-scale structures. The problem is that below a certain range of velocity heterogeneities, deterministic methods become difficult to use, so we turn to a probabilistic approach. With this in view, we characterize the velocity field as a random field defined by its first and second statistical moments. Usually, a seismic random medium is defined as a homogeneous velocity background perturbed by a small random field that is assumed to be stationary. Thus, we make a link between such a random velocity medium (together with a simple reflector) and seismic reflection traveltimes. Assuming that the traveltimes are ergodic, we use 2-D seismic reflection geometry to study the decrease in the statistical traveltime fluctuations as a function of the offset (the source–receiver distance). Our formulae are based on the Rytov approximation and the parabolic approximation for acoustic waves. The validity and the limits are established for both of these approximations in statistically anisotropic random media. Finally, theoretical inversion procedures are developed for the horizontal correlation structure of the velocity heterogeneities for the simplest case of a horizontal reflector. Synthetic seismograms are then computed (on particular realizations of random media) by simulating scalar wave propagation via finite difference algorithms. There is good agreement between the theoretical and experimental results.  相似文献   

10.
Inversion of seismic attributes for velocity and attenuation structure   总被引：1，自引：0，他引：1  
We have developed an inversion formuialion for velocity and attenuation structure using seismic attributes, including envelope amplitude, instantaneous frequency and arrival times of selected seismic phases. We refer to this approach as AFT inversion for amplitude, (instantaneous) frequency and time. Complex trace analysis is used to extract the different seismic attributes. The instantaneous frequency data are converted to t * using a matching procedure that approximately removes the effects of the source spectra. To invert for structure, ray-perturbation methods are used to compute the sensitivity of the seismic attributes to variations in the model. An iterative inversion procedure is then performed from smooth to less smooth models that progressively incorporates the shorter-wavelength components of the model. To illustrate the method, seismic attributes are extracted from seismic-refraction data of the Ouachita PASSCAL experiment and used to invert for shallow crustal velocity and attenuation structure. Although amplitude data are sensitive to model roughness, the inverted velocity and attenuation models were required by the data to maintain a relatively smooth character. The amplitude and t * data were needed, along with the traveltimes, at each step of the inversion in order to fit all the seismic attributes at the final iteration.  相似文献   

11.
Summary. Six aftershock sequences in Iran and Turkey are relocated using existing teleseismic data. Two of these are in the Zagros mountains where local fieldwork has failed to detect subcrustal seismicity but published teleseismic locations show depths greater than 100 km. All apparently deep events are shown to be small and badly recorded with poor depth resolution. There is thus no evidence for active lithospheric subduction in the Zagros.
Relocations of other sequences in Iran and Turkey are used with fault plane solutions, satellite photographs and surface faulting to provide new insight on the geometry of faulting and crustal deformation of those regions. Linear seismic trends from these sequences are shown to cut older geological structures and do not always bear a simple relation to surface faulting. In such cases aftershock activity may be on primary buried faults whose behaviour is not simply revealed in surface structure and deformation.
A linearized inversion scheme is used to investigate the trade-off between resolution and uncertainty in the hypocentral parameters. The ultimate resolution of teleseismic locations is shown to be limited by the quality of arrival time data.  相似文献   

12.
Two-dimensional magnetotelluric inversion   总被引：2，自引：0，他引：2  
Summary. When complex structure is encountered in magnetotelluric surveys, interpretation by locally fitted layered models is of questionable validity. However, when the processed data show two-dimensional structure, numerical inversion schemes for two-dimensional models may be constructed as an aid to regional data interpretation.
The two-dimensional magnetotellurics inversion problem is here formulated in a way that may be applied to many problems. A resulting computer program is analysed carefully in terms of its cost relative to that of simpler layered modelling.
As an example, the method is applied to some field data where the interpretive advantages of the program become evident.  相似文献   

13.
The BABEL marine seismic experiment has been carried out to investigate the lithospheric structure and antecedent tectonic signatures of the Baltic Shield, including the Archaean-Proterozoic collisional structure in the northern part of the Gulf of Bothnia.
Lithospheric seismic-reflection streamer data and simultaneously recorded wide-angle reflection and refraction data collected in the Gulf of Bothnia as part of the BABEL project have been used for 3-D modelling. The distribution of land stations around the Gulf provides a good 3-D ray coverage of the PMP reflection data recorded at the eight stations in the area and allows an estimation of strikes and dips of the Moho boundary in the area. The traveltimes of reflected phases are calculated using a method that utilizes the finite-difference solution of the eikonal equation. The Moho wide-angle-reflection (PMP) traveltimes are modelled using an inversion method. A 2-D model from the Gulf of Bothnia extended into the third dimension is used as an initial model. During the inversion the velocity is kept constant and only the Moho boundary is allowed to vary. To estimate the strike of the Moho boundary and the stability of the inversion, two initial models with different strikes are examined.
The results indicate that the Moho depth in the Gulf of Bothnia undulates and has a maximum depth of 55 km in the south, rising to 42 km in the north. The Moho depth variations seem to be step-like. This change in the Moho depth coincides with the location of the presumed fossil subduction zone in the area. The crustal-thickness variations seem to be well approximated by a nearly 2-D structure striking parallel to a postulated subduction zone immediately to the south of the Skellefte area. The presence of the step at the crust/mantle boundary can be interpreted as a result of a plate-collision event at about 2 Ga.  相似文献   

14.
This paper presents a method to invert underside-reflection ( P d P or S d S arrivals) data for lateral depth variations of upper-mantle discontinuities, combining traveltime and amplitude data. The method greatly improves the resolution of small-scale undulations obtained by existing imaging methods and does not suffer from the long-wavelength biases that are likely to be present in currently available models. Existing inversion methods account for the large size of the Fresnel zone of underside reflections, but not for its complexity, arising from the mini-max traveltime nature of PP- and SS -related waves. This neglect results in long-wavelength artefacts from small-scale undulations of the discontinuities, obscuring true long-wavelength depth variations. The inversion method presented in this paper uses a complex-valued sensitivity kernel, derived from the representation of underside reflections through a Kirchhoff integral formulation. The sensitivity kernel accounts for the varying sensitivity of the waveforms to discontinuity structure over the Fresnel zone. The method is applied to a large, synthetic data set. The data set consists of P d P amplitudes and traveltimes. The results show that the new inversion method resolves depth variations on a lateral scale that is smaller than the size of the Fresnel zone of individual underside reflections (but larger than the dominant wavelength), retaining the resolution of large-scale variations. The results presented here suggest that the discontinuity depth variations induced by slab penetration of the 670 discontinuity could be resolved by current broad-band P 670 P data sets.  相似文献   

15.
A new algorithm is presented for the integrated 2-D inversion of seismic traveltime and gravity data. The algorithm adopts the 'maximum likelihood' regularization scheme. We construct a 'probability density function' which includes three kinds of information: information derived from gravity measurements; information derived from the seismic traveltime inversion procedure applied to the model; and information on the physical correlation among the density and the velocity parameters. We assume a linear relation between density and velocity, which can be node-dependent; that is, we can choose different relationships for different parts of the velocity–density grid. In addition, our procedure allows us to consider a covariance matrix related to the error propagation in linking density to velocity. We use seismic data to estimate starting velocity values and the position of boundary nodes. Subsequently, the sequential integrated inversion (SII) optimizes the layer velocities and densities for our models. The procedure is applicable, as an additional step, to any type of seismic tomographic inversion.
We illustrate the method by comparing the velocity models recovered from a standard seismic traveltime inversion with those retrieved using our algorithm. The inversion of synthetic data calculated for a 2-D isotropic, laterally inhomogeneous model shows the stability and accuracy of this procedure, demonstrates the improvements to the recovery of true velocity anomalies, and proves that this technique can efficiently overcome some of the limitations of both gravity and seismic traveltime inversions, when they are used independently.
An interpretation of field data from the 1994 Vesuvius test experiment is also presented. At depths down to 4.5 km, the model retrieved after a SII shows a more detailed structure than the model obtained from an interpretation of seismic traveltime only, and yields additional information for a further study of the area.  相似文献   

16.
Summary. A method of synthetic seismogram computation for teleseismic SV -waves is developed in order to treat quantitatively SV -waves in problems of body wave source inversion and source—receiver structure studies. The method employs WKBJ theory for a generalized ray in a vertically inhomogeneous half-space and the propagator matrix technique for waves in near-surface homogeneous layers. Wavenumber integration is done along the real axis of the wavenumber plane and anelasticity is included by using complex velocity in all regions of the earth model. The near-surface source structure is taken into account in the computation for the case of the shallow source by allowing a point source to be located in the homogeneous layers. Source and receiver area structures are also allowed to differ. A general moment tensor point source is considered.  相似文献   

17.
Summary. Moment tensor inversion methods can be applied with success in the determination of source properties of simple earthquakes. However, these methods utilize the assumption of a point source, which is inadequate for modelling many complicated, shallow earthquakes. For complex earthquakes, an inversion using finite faulting models is desirable but the number of parameters involved requires that a good starting model be found or that independent constraints be placed on some of the parameters. A method is presented for low-pass filtering both the data and Green's functions, passing only signals with wavelengths greater than the dimension of the entire fault. The filter tends to smooth complications in the waveforms and allows application of the point source moment tensor inversion. This method is applied to body waves from the 1978 Thessaloniki, Greece, earthquake, the 1971 San Fernando earthquake and to a multiple-point source synthetic model of the San Fernando event. For the Thessaloniki event, although a multiple-source mechanism has been suggested, inversion results before and after filtering were essentially identical, indicating that a point source mechanism is sufficient in modelling the long-period, teleseismic body waves. In the case of the San Fernando earthquake, the point source Green's functions were incapable of simultaneously modelling the P - and SH -waves. Inversion of P -waves alone resulted in extreme parameter resolution problems, but allowed constraint in one axis of the moment tensor and suggested an overall source time function. Inversion of a synthetic San Fernando data set yielded similar results, but allowed an investigation of the shortcomings of the method under controlled circumstances. Although the results may require substantial interpretation, the method presented represents a simple first step in the analysis of complex earthquakes.  相似文献   

18.
Summary. Various factors can make it difficult to explain observations of Love- and Rayleigh-wave dispersion with the same relatively simple isotropic model. These factors include systematic errors which might occur in determinations of observed group and phase velocities, lateral variations in structure along the path of travel, and the attempt to explain observations with a model comprised of only a small number of thick layers. The last of these factors is illustrated by an inversion of dispersion data in the central United States where shear-wave anisotropy had previously been invoked as one way to explain incompatible Love- and Rayleigh-wave velocities. It is shown that the data can be satisfied equally well by an isotropic model consisting of several thin layers.
In cases where the incompatibility of Love- and Rayleigh-wave data might be produced by intrinsic anisotropy, it is necessary to invert those data using an anisotropic theory rather than by separate isotropic inversions of Love and Rayleigh waves. Inversions of fundamental-mode data for a region of the Pacific, assuming anisotropic media in which the layers are transversely isotropic with a vertical axis of symmetry, lead to models which are highly non-unique. Even if the inversions solve only for shear velocities in the litho-sphere and asthenosphere it is not possible, without supplementary information, to ascertain the depth interval over which anisotropy occurs or to determine the thickness of the lithosphere or asthenosphere with much precision.  相似文献   

19.
选取新疆奇台县的134个土壤样本,利用土壤反射率对数的一阶导数光谱分别对4 种小波函数进行多层离散分解,采用PLSR方法分别建立了土壤速效钾含量的反演模型,并对其精度值进行检验。结果表明：小波分解获得的各层低频系数以1～3层较高,而其余各层则较低。所有函数分解的6层中,均以第2层低频系数建模的精度最高,随着分解层数（>2层）的增加,其精度值和显著性明显降低。相同尺度下,采用4种小波函数的低频系数构建的反演模型的精度差异较小,而Bior1.3为最优函数;基于Bior 1.3分解的ca2低频系数建模的R2达0.964,RMSE仅为8.19 mg·kg-1,且为极显著水平,为最佳反演模型,经样本检验后发现,此模型可用以快速、准确估算土壤高光谱速效钾含量。  相似文献   

20.
The Massif Central, the most significant geomorphological unit of the Hercynian belt in France, is characterized by graben structures which are part of the European Cenozoic Rift System (ECRIS) and also by distinct volcanic episodes, the most recent dated at 20 Ma to 4000 years BP. In order to study the lithosphere-asthenosphere system beneath this volcanic area, we performed a teleseismic field experiment.
During a six-month period, a joint French-German team operated a network of 79 mobile short-period seismic stations in addition to the 14 permanent stations. Inversion of P -wave traveltime residuals of teleseismic events recorded by this dense array yielded a detailed image of the 3-D velocity structure beneath the Massif Central down to 180 km depth. The upper 60 km of the lithosphere displays strong lateral heterogeneities and shows a remarkable correlation between the volcanic provinces and the negative velocity perturbations. The 3-D model reveals two channels of low velocities, interpreted as the remaining thermal signature of magma ascent following large lithospheric fractures inherited from Hercynian time and reactivated during Oligocene times. The teleseismic inversion model yields no indication of a low-velocity zone in the mantle associated with the graben structures proper. The observation of smaller velocity perturbations and a change in the shape of the velocity pattern in the 60–100 km depth range indicates a smooth transition from the lithosphere to the asthenosphere, thus giving an idea of the lithosphere thickness. A broad volume of low velocities having a diameter of about 200 km from 100 km depth to the bottom of the model is present beneath the Massif Central. This body is likely to be the source responsible for the volcanism. It could be interpreted as the top of a plume-type structure which is now in its cooling phase.  相似文献   

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

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