首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 46 毫秒
1.
Summary. Use of paraxially approximated Gaussian beams continues to be actively pursued for construction of synthetic seismograms in complicated environments. How to select the beams in the stack remains a source of difficulty which has primarily been addressed by semi-heuristic considerations. In this paper, the classical example of line-source field reflection from a homogeneous half-space that can sustain a head wave is examined from a plane-wave spectral point of view. The individual beam fields are modelled exactly by the complex source point technique, which emphasizes the complex spectral content of these wave objects. The quality of the paraxial approximation of a typical reflected (Gaussian) beam characterized by different parameters is examined from this perspective, and is compared with uniform and non-uniform asymptotics generated from the exact beam field spectral integral. With this information as background, the reflected field for a real line-source is synthesized by beam superposition. Except for the immediate vicinity of the critical reflection angle, the well-known failure of narrow paraxial beams, no matter how densely stacked, to reproduce the head wave effects is shown to be due to the inadequate spectral content of these beams and not to the failure of beam stacking per se. When the rigorous solutions are used for the narrow-waist beams, even relatively few suffice to yield agreement with the exact solution. This circumstance emphasizes the importance of fully understanding the spectral implications of various beam stacking schemes.  相似文献   

2.
Summary. The Green's function, in a constant gradient medium, is derived for an explosive point source, in the frequency and the time domains. The analytical dynamic ray tracing (DRT) solution is rederived with conditions stated in Part I. The Gaussian beam (GB) solution is investigated. New beam parameters and conditions are defined. Comparisons between exact and approximate solutions are undertaken.
For both methods, DRT and GB, conditions of validity are explicit and quantitative. An accuracy criterion is defined in the time domain, and measures a global relative error. The range of validity is expressed in the form of two inequalities for the dynamic ray tracing method and of five inequalities for the Gaussian beam method. Results remain accurate at ray turning points. For the types of medium considered, the breakdown of the dynamic ray tracing method is smoother and better behaved than that of Gaussian beams. As examples, a vertical seismic profiling configuration, and a shallow earthquake are modelled, using Gaussian beams.  相似文献   

3.
Gaussian beams in elastic 2-D laterally varying layered structures   总被引：2，自引：0，他引：2  
Summary. In a paper by Červený & Pšenčik, high-frequency Gaussian beams in elastic 2-D, laterally inhomogeneous, smooth media were investigated as asymptotic high-frequency solutions of elastodynamic equations, concentrated close to rays of P - and S -waves. This paper generalizes the above results for 2-D, laterally inhomogeneous, layered structures. Gaussian beams concentrated close to any multiply-reflected, possibly converted, ray are investigated. Gaussian beams are regular everywhere, including caustic regions. The paraxial ray approximation, which allows the wavefield in the zero-order ray approximation to be evaluated not only directly on the ray, but also in its vicinity, is derived as a limiting case of the Gaussian beams.  相似文献   

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

5.
Summary. Some relations between Gaussian beams, complex rays and the analytic extension of the Green's function in smoothly inhomogeneous media are shown in this paper. It is found that: (1) a single Gaussian beam is a paraxial approximation of the analytical extension of the ray-approximated Green's function in smoothly inhomogeneous media by putting the source point into a complex space. The Gaussian beam approximation of the Green's function has an advantage in computational efficiency and stability and can avoid the singularity problems at caustics, but also introduces a parabolic approximation to the wavefront and an angle-dependent amplitude damping. Therefore the validity of the Gaussian beam approximation should be checked using other methods. (2) Complex-ray tracing, which does not involve the paraxial approximation, can also avoid the singularity problemsm though without the computational efficiency. Therefore, it should be used to verify the Gaussian beam approximation, whenever possible. (3) The decomposition of a plane wave into an ensemble of Gaussian beams is equivalent to approximating the Green's function (the kernel of the ray-Kirchhoff method) with a single Gaussian beam. This introduces a parabolic approximation to the wavefront and a Gaussian windowing for arrival angles which may cause some problems in modelling wave propagation and scattering and has no advantages over other methods. (4) The representation of a point source field by a superposition of Gaussian beams, on the other hand, is equivalent to approximating the Green's function with a bundle of overlapped Gaussian beams. This representation is similar to a Maslov uniform asymptotic representation. It has no caustics and has improved accuracies at the caustics for quasi-plane waves compared to the extended WKBJ method.  相似文献   

6.
Summary. Asymptotic ray theory (ART) fails in transition regions near critically reflected, bottom glancing or caustic-forming rays in a vertically inhomogeneous layered earth. These deficiencies are repaired here by replacing the transitional ray fields with guided modes plus truncation remainders. Exact ray-mode equivalences and their high-frequency asymptotic approximations are formulated, and their validity and efficiency are verified by numerical comparisons for SH motion in a two-layer earth model comprised of an inhomogeneous sediment above an homogeneous semi-infinite bedrock.  相似文献   

7.
Summary. An asymptotic procedure for the computation of wave fields in two-dimensional laterally inhomogeneous media is proposed. It is based on the simulation of the wave field by a system of Gaussian beams. Each beam is continued independently through an arbitrary inhomogeneous structure. The complete wave field at a receiver is then obtained as an integral superposition of all Gaussian beams arriving in some neighbourhood of the receiver. The corresponding integral formula is valid even in various singular regions where the ray method fails (the vicinity of caustic, critical point, etc.). Numerical examples are given.  相似文献   

8.
Some remarks on the Gaussian beam summation method   总被引：1，自引：0，他引：1  
Summary. Recently, a method using superposition of Gaussian beams has been proposed for the solution of high-frequency wave problems. The method is a potentially useful approach when the more usual techniques of ray theory fail: it gives answers which are finite at caustics, computes a nonzero field in shadow zones, and exhibits critical angle phenomena, including head waves. Subsequent tests by several authors have been encouraging, although some reported solutions show an unexplained dependence on the 'free' complex parameter ε which specifies the initial widths and phases of the Gaussian beams.
We use methods of uniform asymptotic expansions to explain the behaviour of the Gaussian beam method. We show how it computes correctly the entire caustic boundary layer of a caustic of arbitrary complexity, and computes correctly in a region of critical reflection. However, the beam solution for head waves and in edge-diffracted shadow zones are shown to have the correct asymptotic form, but with governing parameters that are explicitly ε-dependent. We also explain the mechanism by which the beam solution degrades when there are strong lateral inhomogeneities. We compare numerically our predictions for some representative, model problems, with exact solutions obtained by other means.  相似文献   

9.
Summary. Two methods of computing body wave synthetic seismograms in three-dimensional laterally varying media are discussed. Both these methods are based on the summation of Gaussian beams. In the first, the initial beam parameters are chosen at the source, in the second at the beam endpoints. Both these variants eliminate the ray method singularities. The expansion of the wavefield into plane waves may be considered as the limiting case of the first approach and the Chapman–Maslov method as the limiting case of the second approach. Computer algorithms are briefly described and numerical examples presented. In the first numerical example, the comparisons of the two approaches, based on summing Gaussian beams, with the reflectivity method indicate that the computed synthetic seismograms are satisfactorily accurate even in the caustic region. The next example suggests that the two methods discussed can be simply and effectively applied to 3-D laterally inhomogeneous structures.  相似文献   

10.
Summary. Several important applications of the paraxial ray approximation (PRA) to numerical modelling of high-frequency seismic body wavefields are discussed. The PRA can be used to evaluate the displacement vector not only directly on the ray, as in the standard ray method. but also approximately in the vicinity of this ray. The PRA also offers simple ways of approximate evaluation of paraxial rays, situated in the vicinity of the central ray, and of two-point ray tracing. A very important application of the PRA consists in a simple, fast and effective Computation of body-wave synthetic seismograms in general, 3-D, laterally inhomogeneous, layered structures. Examples of synthetic seismograms for 3-D structures, computed using the PRA, are presented.  相似文献   

11.
Summary. This paper discusses several aspects of the calculation of theoretical seismograms for two-dimensional inhomogeneous media with the method of Gaussian beams. The most important steps of this method, kinematic and dynamic ray tracing, can be performed very efficiently, if the model cross-section is subdivided into triangles with linear velocity laws. Each Gaussian beam is characterized by a complex beam constant ε which determines its width and phase-front curvature. Various possibilities to choose ε are discussed, including cases where beam properties at the beam endpoint (and not at the beginning) are prescribed; for instance, the beam width at the endpoint can be specified. In such cases the beam constant is a function of the radiation angle at the source, and the decomposition of a cylindrical wave into beams has to take this into account by weighting the beams differently, at least in principle. The exact weight function is derived and shown to be reasonably well approximated by the weight function, corresponding to angle-independent ε Theoretical seismograms are presented for a laterally heterogeneous model of the crust–mantle transition which is characterized by complications in the reflection from the transition and in the refraction from below. These complications are modelled by and large with success. The seismograms, however, depend to a certain extent on the choice of the beam constant. Moreover, according to the reciprocity principle calculations with source and receiver interchanged should have the same results as calculations for the original configuration. In practice this is not so, and the difference increases with the strength of lateral heterogeneities. Hence, for a successful application of Gaussian beams the model should not vary too strongly in lateral direction.  相似文献   

12.
  相似文献   

13.
Summary. An algorithm for the computation of travel times, ray amplitudes and ray synthetic seismograms in 3-D laterally inhomogeneous media composed of isotropic and anisotropic layers is described. All 21 independent elastic parameters may vary within the anisotropic layers. Rays and travel times are evaluated by numerical solution of the ray tracing equations. Ray amplitudes are determined by evaluating reflection/ transmission coefficients and the geometrical spreading along individual rays. The geometrical spreading is computed approximately by numerical measurement of the cross-sectional area of the ray tube formed by three neighbouring rays. A similar approximate procedure is used for the determination of the coefficients of the paraxial ray approximation. The ray paraxial approximation makes computation of synthetic seismograms on the surface of the model very efficient. Examples of ray synthetic seismograms computed with a program package based on the described algorithm are presented.  相似文献   

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

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

16.
Summary. Several approaches to computing body wave seismograms in 2–D and 3–D laterally inhomogeneous layered structures are suggested. They are based on the Gaussian beam method, which has been recently applied to the evaluation of time-harmonic high-frequency wavefields in inhomogeneous media. Three variants are discussed in some detail: the spectral method, the convolutory method and the wave-packet method. The most promising seems to be the wave-packet approach. In this approach, the wavefield, generated by a source, is expanded into a system of wave packets, which propagate along rays from the source in all directions. The wave packets change their properties due to diffusion, spreading, reflections/transmissions, etc. The resulting seismogram at any point of the medium is then obtained as a superposition of those packets which propagate close to the point. The final expressions in all the three methods are regular even in regions, in which the ray method fails, e.g. in the vicinity of caustics, in the critical region, at boundaries between shadow and illuminated regions, etc. Moreover, they are not as sensitive to the minor details of the medium as the ray method and, what is more, they remove the time-consuming two-point ray tracing from computations. Numerical examples of synthetic seismograms computed by the wave-packet approach are presented.  相似文献   

17.
Data recorded by the Italian Telemetered Seismic Network (ITSN) of the Istituto Nazionale di Geofisica (ING) have been widely used in recent years to image slab structures and to find evidence for active processes along the Italian Peninsula. However, the use of seismic data for geostructural purposes may be affected by the well-known trade-off between earthquake location and seismic-velocity parameters. Furthermore, the confidence ellipse predicted by standard procedures may be inadequate for the representation of the probable error of a computed localization. This paper evaluates the probable errors on the hypocentre determinations of the seismic events recorded by the ITSN, using a Monte Carlo method.
We compute synthetic arrival times using a 1-D velocity model appropriate as an average for the Italian area. The hypocentres used are all those recorded by the ITSN during the period January 1992 to March 1994 (1972 events). Station locations are those of the current ITSN configuration. The synthetic arrival times are perturbed with a Gaussian distribution of errors and input to ING's standard hypocentral location procedure, but using crustal velocities differing by 10 per cent from those used to generate them. Each simulation is repeated at least 30 times. Average absolute shifts of hypocentres are assessed in grid cells of linear dimension 33 km covering the whole Italian region.
For regions within the ITSN, shifts are typically 5–10 km in location and up to 20 km in depth. However, for offshore and coastal regions, they are much greater: 50 km or more in both location and depth (far exceeding the equivalent uncertainties quoted by ING bulletins). Possible consequences of this are highlighted by producing a cross-section of subcrustal hypocentres from the Adriatic to the Tyrrhenian Sea, where the large uncertainty in depth precludes any confident interpretation of dipping tectonic features.  相似文献   

18.
Summary. Asymptotic ray theory is applied to surface waves in a medium where the lateral variations of structure are very smooth. Using ray-centred coordinates, parabolic equations are obtained for lateral variations while vertical structural variations at a given point are specified by eigenfunctions of normal mode theory as for the laterally homogeneous case. Final results on wavefields close to a ray can be expressed by formulations similar to those for elastic body waves in 2-D laterally heterogeneous media, except that the vertical dependence is described by eigenfunctions of 'local' Love or Rayleigh waves. The transport equation is written in terms of geometrical-ray spreading, group velocity and an energy integral. For the horizontal components there are both principal and additional components to describe the curvature of rays along the surface, as in the case of elastic body waves. The vertical component is decoupled from the horizontal components. With complex parameters the solutions for the dynamic ray tracing system correspond to Gaussian beams: the amplitude distribution is bell-shaped along the direction perpendicular to the ray and the solution is regular everywhere, even at caustics. Most of the characteristics of Gaussian beams for 2-D elastic body waves are also applicable to the surface wave case. At each frequency the solution may be regarded as a set of eigenfunctions propagating over a 2-D surface according to the phase velocity mapping.  相似文献   

19.
Gaussian beams in two-dimensional elastic inhomogeneous media   总被引：1，自引：0，他引：1  
Summary. Asymptotic high-frequency solutions of elastodynamic equations in two-dimensional laterally inhomogeneous media concentrated close to rays of P - and S -waves are investigated. From a physical point of view, these vectorial solutions correspond to Gaussian beams; the amplitude distribution of their principal components is bell-shaped along the direction perpendicular to the ray. The principal component of the elastodynamic Gaussian beam is controlled by the parabolic equation, which has exactly the same form as the parabolic equation for scalar Gaussian beams. The elastodynamic Gaussian beams are regular everywhere, including caustics.  相似文献   

20.
We present new 3D seismic and well data from the Ebro Margin, NW Mediterranean Sea, to shed new light on the processes that formed the Messinian Erosion Surfaces (MES) of the Valencia Trough (Mediterranean Sea). We combine these data with backstripping techniques to provide a minimum estimate of the Messinian sea level fall in the EBRO Margin, as well as coupled isostasy and river incision and transport modeling to offer new constraints on the evolution of the adjacent subaerial Ebro Basin. Four major seismic units are identified on the Cenozoic Ebro Margin, based on the seismic data, including two major prograding megasequences that are separated by a major unconfirmity: the MES. The 3D seismic data provide an unprecedented view of the MES and display characteristic features of subaerial incision, including a drainage network with tributaries of at least five different orders, terraces and meandering rivers. The Messinian landscape presents a characteristic stepped‐like profile that allows the margin to be subdivided in three different regions roughly parallel to the coastline. No major tectonic control exists on the boundaries between these regions. The boundary between the two most distal regions marks the location of a relatively stable base level, and this is used in backstripping analysis to estimate the magnitude of sea level drop associated with the Messinian Salinity Crisis on the Ebro Margin. The MES on the Ebro Margin is dominated by a major fluvial system, that we identify here as the Messinian Ebro River. The 3D seismic data, onshore geology and modeling results indicate that the Ebro River drained the Ebro Basin well in advance of the Messinian.  相似文献   

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

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