首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 62 毫秒
1.
A method for calculating synthetic seismograms in laterally varying media   总被引：2，自引：0，他引：2  
Summary An effective algorithm for computing synthetic seismograms in laterally inhomogeneous media has been developed. The method, based on zero-order asymptotic ray theory, is primarily intended for use in refraction and reflection studies and provides an economical means of seismic modelling.
A given smoothed velocity-depth-distance model is divided into small squares with constant seismic parameters and first-order interfaces are represented by an arbitrary number of dipping linear segments. The computation of ray propagation and amplitudes through such a model does not involve complicated analytic expressions and therefore minimizes computer time.
Amplitudes are determined by geometrical spreading of spherical wave-fronts and energy partitioning at interfaces. Synthetic seismograms calculated for laterally homogeneous models are in good agreement with those obtained by the Reflectivity Method.  相似文献   

2.
Reflection coefficients for weak anisotropic media   总被引：1，自引：0，他引：1  
The interaction of plane elastic waves with a plane boundary between two anisotropic elastic half-spaces is investigated. The anisotropy dealt with in this study is of a general type. Explicit expressions for energy-related reflection and transmission coefficients are derived. They represent an approximation which is valid for a small deviation of the elastic parameters from isotropy.
Classical perturbation theory is applied on a 6times6 non-symmetric real eigenvalue problem to calculate first-order corrections for the polarization and stress of the plane waves. The explicit solution of the isotropic problem is used as a reference case. Degenerate perturbation theory is used to consider the splitting of the isotropic S -wave into two anisotropic qS-waves. The boundary conditions for two half-spaces in welded contact lead to a 6times6 system of linear equations. A correction to the isotropic solution is calculated by linearization. The resultant coefficients are functions of horizontal slowness, Lamé parameters and densities of the reference media, and of the perturbation of the elasticity tensors from isotropy.  相似文献   

3.
We construct a catalogue of all the possible elementary point sources of seismic waves. There are three general classes of sources, two spheroidal and one toroidal. We consider excitation functions for these point-like sources as well as for sources of finite size in far-, intermediate- and near-field for an infinite homogeneous isotropic medium. The sources corresponding to seismic-moment tensors for the second-, third- and fourth-ranks are considered in more detail; we identify 10 different seismic sources in this range: one monopole, two or three dipoles, three quadrupoles, etc. For the step-function of the scalar seismic-moment release, the amplitude spectrum for the third-rank sources is proportional to the angular frequency ω in the region below the corner frequency ω cr. The fourth-rank sources have an ω 2 spectrum in the same range. The possibility of separate and simultaneous inversion of seismic body-wave data and static deformation data for sources of different order is discussed. Some equivalent-force moment higher-rank sources are 'shielded' by lower-rank sources of the same order; the former sources cannot be inverted from seismic data without additional assumptions. Because of their simple radiation pattern, the lower order multipoles, i.e. the monopole and dipoles, are the first sources other than the double-couple which should be considered for inversion.  相似文献   

4.
In this article the interaction of plane waves with a weak-contrast interface between two weakly anisotropic half-spaces is investigated. The anisotropy dealt with is of a general type. The stress–displacement vectors of the plane waves are calculated by perturbation theory. By assuming that the jump in elastic parameters and density across the interface is small, one can derive a simple expression for the R qPqP coefficient. In cases in which the wave motion is restricted to a symmetry plane of an anisotropic medium, simple expressions for the R qSVqSV and R SHSH coefficients are also derived.  相似文献   

5.
We derive asymptotic formulae for the toroidal and spheroidal eigenfrequencies of a SNREI earth model with two discontinuities, by considering the constructive interference of propagating SH and P-SV body waves. For a model with a smooth solid inner core, fluid outer core and mantle, there are four SH and 10 P-SV ray parameters regimes, each of which must be examined separately. The asymptotic eigenfrequency equations in each of these regimes depend only on the intercept times of the propagating wave types and the reflection and transmission coefficients of the waves at the free surface and the two discontinuities. If the classical geometrical plane-wave reflection and transmission coefficients are used, the final eigenfrequency equations are all real. In general, the asymptotic eigenfrequencies agree extremely well with the exact numerical eigenfrequencies; to illustrate this, we present comparisons for a crustless version of earth model 1066A.  相似文献   

6.
Summary. The limitations of asymptotic wave theory and its geometrical manifestations are newly formalized and scrutinized. Necessary and sufficient conditions for the existence of acoustic and seismic rays and beams in general inhomogeneous media are expressed in terms of new physical parameters: the threshold frequency ω0 associated with the P/S decoupling condition, the cut-off frequency ωc associated with the radiation-zone condition, the total curvature of the wavefront and the Fresnel-zone radius.
The analysis is facilitated with the introduction of a new ancillary functional – the hypereikonal which is capable of representing ordinary as well as evanescent waves. The hypereikonal is the natural extension of the eikonal theory.
With the aid of the above new parameters, simple conditions are obtained for the decoupled far field, the decoupled near field, two point dynamic ray tracing, paraxial wavefields and Gaussian beams.  相似文献   

7.
We derive both 3-D and 2-D Fréchet sensitivity kernels for surface-wave group-delay and anelastic attenuation measurements. A finite-frequency group-delay exhibits 2-D off-ray sensitivity either to the local phase-velocity perturbation  δ c / c   or to its dispersion  ω(∂/∂ω)(δ c / c )  as well as to the local group-velocity perturbation  δ C / C   . This dual dependence makes the ray-theoretical inversion of measured group delays for 2-D maps of  δ C / C   a dubious procedure, unless the lateral variations in group velocity are extremely smooth.  相似文献   

8.
Summary. The symmetry relations between the reflection and transmission coefficients for plane elastic waves incident upon an arbitrary horizontally stratified medium are derived by a novel approach. Previous results, particularly for a single interface, are obtained as special cases of this treatment.
In addition, for perfectly elastic media, projection operators for travelling and evanescent waves are introduced and used to derive a number of new relationships between the reflection and transmission coefficients.  相似文献   

9.
Summary Reflectivity and ray theories are united to produce a hybrid technique of computing synthetic seismograms for a plane layered medium in subcritical regions. Numerical experiments have indicated that this technique is useful when the depth structure is one composed of thick layers separated by finely layered zones. As the theory for wave propagation in a plane layered medium is well known, the simple SH case is investigated so that the basic idea of the method may be conveyed without an excess of mathematics that would be necessitated if the P-SV problem were considered.
In computing the ray-reflectivity seismogram, the thick layers are treated using asymptotic ray theory while the thin-layered zones are treated as quasiinterfaces where analogues of reflection and transmission coefficients called reflectivities and transmittivities are calculated utilizing a Thomson-Haskell formulation. A stationary phase approximation is employed when evaluating the integral which gives the displacement due to an arbitrary ray propagating in the thick layers of the above-mentioned medium, and the validity of this approximation is discussed.
A comparison of ray, numerical integration (reflectivity) and ray-reflectivity synthetic sections indicates that this method yields quite acceptable results for subcritical reflection work and is suitable for application in seismic interpretation as individual arrivals associated with ray-paths in the thick layers may be identified. Furthermore, the method is quite cost efficient and may be extended to a medium where the thick layers are non-planar using asymptotic ray theory in these layers.  相似文献   

10.
We present approximate displacement and energy PP and PS reflection/transmission coefficients for weak-contrast interfaces in general weakly anisotropic elastic media. The coefficients were obtained by applying first-order perturbation theory and then expressed in a compact and relatively simple form. The formulae can be used for arbitrary orientations of the incidence plane and interface, without the need to transform the elasticity parameters to a local Cartesian coordinate system. The accuracy of the approximate formulae is illustrated for the PS reflection coefficient for two synthetic models. For these models, we also study the possibility of using the approximate PP reflection coefficient in the inverse problem.  相似文献   

11.
Okada (1992) provided expressions for the displacement and strain fields due to a finite rectangular source in an elastic, homogeneous and isotropic half-space. Starting with these results, we applied the correspondence principle of linear viscoelasticity to derive the quasi-static displacement, strain and stress fields in a viscoelastic, homogeneous and isotropic half-space. We assume that the medium deforms viscoelastically with respect to both the shear and the normal stresses but keeps a constant bulk modulus; in particular, the shear modulus relaxes as Maxwell fluid. We presented the viscoelastic effect on displacement, displacement gradient and stress fields, for a choice of parameter values. The viscoelastic effect due to the sudden dislocation reaches a limit value after about 10 times the Maxwell time. The expressions obtained here provide tools for the study of viscoelastic relaxation of lithosphere associated with seismic and volcanic phenomena.  相似文献   

12.
Seismic waves in a stratified half space   总被引：22，自引：0，他引：22  
Summary. The response of a stratified elastic half space to a general source may be represented in terms of the reflection and transmission properties of the regions above and below the source. For P-SV and SH waves and both buried sources and receivers, convenient forms of the response may be found in which no loss of precision problems arise from growing exponential terms in the evanescent regime. These expressions have a ready physical interpretation and enable useful approximations to the response to be developed. The reflection representation leads to efficient computational procedures for models composed of uniform layers, which may be extended in an asymptotic development to piecewise smooth models.  相似文献   

13.
Non-linear Bayesian joint inversion of seismic reflection coefficients   总被引：2，自引：0，他引：2  
Inversion of seismic reflection coefficients is formulated in a Bayesian framework. Measured reflection coefficients and model parameters are assigned statistical distributions based on information known prior to the inversion, and together with the forward model uncertainties are propagated into the final result. This enables a quantification of the reliability of the inversion. Quadratic approximations to the Zoeppritz equations are used as the forward model. Compared with the linear approximations the bias is reduced and the uncertainty estimate is more reliable. The differences when using the quadratic approximations and the exact expressions are minor. The solution algorithm is sampling based, and because of the non-linear forward model, the Metropolis–Hastings algorithm is used. To achieve convergence it is important to keep strict control of the acceptance probability in the algorithm. Joint inversion using information from both reflected PP waves and converted PS waves yields smaller bias and reduced uncertainty compared to using only reflected PP waves.  相似文献   

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

15.
The phase velocity and the attenuation coefficient of compressional seismic waves, propagating in poroelastic, fluid-saturated, laminated sediments, are computed analytically from first principles. The wavefield is found to be strongly affected by the medium heterogeneity. Impedance fluctuations lead to poroelastic scattering; variations of the layer compressibilities cause inter-layer flow (a 1-D macroscopic local flow). These effects result in significant attenuation and dispersion of the seismic wavefield, even in the surface seismic frequency range, 10–100 Hz. The various attenuation mechanisms are found to be approximately additive, dominated by inter-layer flow at very low frequencies. Elastic scattering is important over a broad frequency range from seismic to sonic frequencies. Biot's global flow (the relative displacement of solid frame and fluid) contributes mainly in the range of ultrasonic frequencies. From the seismic frequency range up to ultrasonic frequencies, attenuation due to heterogeneity is strongly enhanced compared to homogeneous Biot models. Simple analytical expressions for the P -wave phase velocity and attenuation coefficient are presented as functions of frequency and of statistical medium parameters (correlation lengths, variances). These results automatically include different asymptotic approximations, such as poroelastic Backus averaging in the quasi-static and the no-flow limits, geometrical optics, and intermediate frequency ranges.  相似文献   

16.
We use Monte Carlo Markov chains to solve the Bayesian MT inverse problem in layered situations. The domain under study is divided into homogeneous layers, and the model parameters are the conductivity of each layer. We use an a priori distribution of the parameters which favours smooth models. For each layer, the a priori and a posteriori distributions are digitized over a limited set of conductivity values.
  The Markov chain relies on updating the model parameters during successive scanning of the domain under study. For each step of the scanning, the conductivity is updated in one layer given the actual value of the conductivity in the other layers. Thus we designed an ergodic Markov chain, the invariant distribution of which is the a posteriori distribution of the parameters, provided the forward problem is completely solved at each step.
  We have estimated the a posteriori marginal probability distributions from the simulated successive values of the Markov chain. In addition, we give examples of complex magnetotelluric impedance inversion in tabular situations, for both synthetic models and field situations, and discuss the influence of the smoothing parameter.  相似文献   

17.
The dynamic coalescence of two mode II cracks on a planar fault is simulated here using the elastodynamic boundary integral equation method. We focus on the complexity of the resultant slip rate and seismic radiation in the crack coalescence model (CCM) and on the reconstruction of a single crack model (SCM) that can reproduce the CCM waveforms from heterogeneous source parameters rather than coalescence. Simulation results reveal that localized higher slip rates are generated by coalescence as a result of stress interaction between the approaching crack tips. The synthesized seismic radiation exhibits a distinct coalescence phase that has striking similarities to stopping phases in the radiation and propagation properties. The corresponding SCM yields a singular increase in the stress drop distribution, which is accompanied by a sudden decrease in it across the point of coalescence in the CCM. This implies that the generation of high-frequency radiation is more efficient from coalescence than from stopping, although both phenomena exhibit the same strong  ω−2  -type displacement spectra.  相似文献   

18.
Summary. The propagation of surface waves in a laterally varying medium can be described by representing the wavetrain as a superposition of modal contributions for a reference structure. As the guided waves propagate through a heterogeneous zone the modal coefficients needed to describe the wavetrain vary with position, leading to interconversions between modes and reflection into backward travelling modes. The evolution of the modal terms may be described by a set of first-order differential equations which allow for coupling to both forward and backward travelling waves; the coefficients in these equations depend on the differences between the actual structure and the reference structure. This system is established using the orthogonality properties of the modal eigenfunctions and is valid for SH -waves, P - SV -waves and full anisotropy.
The reflected and transmitted wavefields for a region of heterogeneity can be related to the incident wave by introducing reflection and transmission matrices which connect the modal coefficients in these fields to those in the incident wavetrain. By considering a sequence of models with increasing width of heterogeneity we are able to derive a set of Ricatti equations for the reflection and transmission matrices which may be solved by initial value techniques. This avoids an awkward two-point boundary value problem for a large number of coupled equations. The method is demonstrated for 1 Hz Lg - and Sn -waves in a multilayered model for which there are 19 coupled modes.
The method is applicable to three-dimensional heterogeneity, and we are able to show that the interconversion between Love and Rayleigh waves, in the presence of gradients in seismic properties transverse to the propagation path, leads to a net rate of increase of the transverse components of the seismogram at the expense of the other components.  相似文献   

19.
This study describes an examination of surface gravity changes caused by dislocations within a 3-D heterogeneous earth. This new theory is described using six independent dislocations: a vertical strike-slip, two vertical dip-slips perpendicular to each other, and three tensile openings on three perpendicular planes. A combination of the six independent dislocations is useful to compute coseismic gravity changes resulting from an arbitrary seismic source at an arbitrary position. Based on the 3-D lateral inhomogeneous P -wave velocity model, we deduce the 3-D density and S -wave velocity models using the relation of Karato. Finally, numerical computations are performed for a location south of Japan (30°N, 135°E). We calculate the coseismic gravity changes resulting from the six independent dislocations for source depths of 100, 300 and 637 km, respectively. Numerical results show that the maximum 3-D effect varies concomitantly with the dislocation type and the source depth. For seismic problems, the effect of elastic parameter  μ  is dominant.  相似文献   

20.
Summary. The equations describing seismic waves in a stratified earth have a number of symmetry properties, one of which has recently been used by Garmany to derive a simple expression for the inverse of the matrix of eigenvectors appearing in the solution of the equations. We review these symmetries of the wave equation in several notations to demonstrate that the property used by Garmany is distinct from the usual symmetries found in the seismological literature. Like the others, the new symmetry has implications for the reflection and transmission properties of a medium. These implications are briefly reviewed in order to show how the new symmetry is placed relative to the others. A limited discussion of the physical origins of the symmetries is given and, though the overall picture is incomplete, it is noted that the new symmetry yields conservation of energy for reflection/transmission at a single interface in all types of media (anisotropic, lossy, etc.).  相似文献   

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

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