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

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

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

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

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

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

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

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

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

10.
Summary. The Gaussian beam method has recently been introduced into synthetic seismology to overcome shortcomings of the ray method, especially in transition regions due to focusing or diffraction where ray theory fails. One proceeds by discretizing the initial data as a superposition of paraxial Gaussian beams, each of which is then traced through the seismic environment. Since Gaussian beam fields do not diverge in ray transition regions, they are 'uniformly regular' although the quality of this regularity depends on the beam parameters and on the 'numerical distance' which defines the extent of the transitional domain. However, when Gaussian beam patches are used to simulate non-Gaussian initial data, there arise ambiguities due to choice of patch size and location, beam width, etc., which are at the user's disposal. The effects of this arbitrariness have customarily been explored by trial and error numerical experiment but no quantitative recommendations have emerged as yet. As a step toward a priori predictive capability, it is proposed here to perform a systematic study on analytically tractable prototype models of how the parameters and location of a single beam affect the quality of the observed seismic field, especially in ray transition regions. The conversion of ordinary ray fields into beam fields in canonical configurations can be accomplished conveniently by displacing a real source point into a complex coordinate space. Thus, the desired beam solutions can be obtained directly from available ray, and even paraxial ray, fields. Complex ray theory and its implications are reviewed here, with an emphasis on improvements of beam tracking schemes employed at present.  相似文献   

11.
Summary. A high-frequency asymptotic integral expansion of a time-harmonic wavefield into Gaussian beams was derived in a previous paper by Klimeš. The discretization error caused by replacing this integral superposition by a discrete summation of Gaussian beams is estimated in this paper.  相似文献   

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

13.
Finite-frequency sensitivity kernels for head waves   总被引：2，自引：0，他引：2  
Head waves are extremely important in determining the structure of the predominantly layered Earth. While several recent studies have shown the diffractive nature and the 3-D Fréchet kernels of finite-frequency turning waves, analogues of head waves in a continuous velocity structure, the finite-frequency effects and sensitivity kernels of head waves are yet to be carefully examined. We present the results of a numerical study focusing on the finite-frequency effects of head waves. Our model has a low-velocity layer over a high-velocity half-space and a cylindrical-shaped velocity perturbation placed beneath the interface at different locations. A 3-D finite-difference method is used to calculate synthetic waveforms. Traveltime and amplitude anomalies are measured by the cross-correlation of synthetic seismograms from models with and without the velocity perturbation and are compared to the 3-D sensitivity kernels constructed from full waveform simulations. The results show that the head wave arrival-time and amplitude are influenced by the velocity structure surrounding the ray path in a pattern that is consistent with the Fresnel zones. Unlike the 'banana–doughnut' traveltime sensitivity kernels of turning waves, the traveltime sensitivity of the head wave along the ray path below the interface is weak, but non-zero. Below the ray path, the traveltime sensitivity reaches the maximum (absolute value) at a depth that depends on the wavelength and propagation distance. The sensitivity kernels vary with the vertical velocity gradient in the lower layer, but the variation is relatively small at short propagation distances when the vertical velocity gradient is within the range of the commonly accepted values. Finally, the depression or shoaling of the interface results in increased or decreased sensitivities, respectively, beneath the interface topography.  相似文献   

14.
Summary. A class of elastic transition zones are modelled by considering a homogeneous half space overlying an inhomogeneous half space with a bounded and monotonically increasing profile for the rigidity modulus and constant Poisson's ratio and density. Reflected P waves due to a compressional point source in the upper half space are studied in the frequency and time domains by means of numerical contour integration in the complex k plane and the Fast Fourier Transform (FFT). Results from the exact fourth-order elasticity theory are compared with those from the approximate decoupled equations for P and SV waves. Agreement is observed between the two theories at high frequencies beyond the caustic range.  相似文献   

15.
An introduction to Maslov's asymptotic method   总被引：3，自引：0，他引：3  
Summary. Familiar concepts such as asymptotic ray theory and geometrical spreading are now recognized as an asymptotic form of a more general asymptotic solution to the non-separable wave equation. In seismology, the name Maslov asymptotic theory has been attached to this solution. In its simplest form, it may be thought of as a justification of disc-ray theory and it can be reduced to the WKBJ seismogram. It is a uniformly valid asymptotic solution, though. The method involves properties of the wavefronts and ray paths of the wave equation which have been established for over a century. The integral operators which build on these properties have been investigated only comparatively recently. These operators are introduced very simply by appealing to the asymptotic Fourier transform of Ziolkowski & Deschamps. This leads quite naturally to the result that phase functions in different domains of the spatial Fourier transform are related by a Legendre transformation. The amplitude transformation can also be inferred by this method. Liouville's theorem (the incompressibility of a phase space of position and slowness) ensures that it is always possible to obtain a uniformly asymptotic solution. This theorem can be derived by methods familiar to seismologists and which do not rely on the traditional formalism of classical mechanics. It can also be derived from the sympletic property of the equations of geometrical spreading and canonical transformations in general. The symplectic property plays a central role in the theory of high-frequency beams in inhomogeneous media.  相似文献   

16.
Summary. The usual asymptotic methods used to correct the high-frequency solutions of the wave equation are unsatisfactory as they do not give the low-frequency, partial reflections expected from a region of high velocity gradient. A new iterative solution is obtained which uses the first term of the Langer asymptotic expansion as the zeroth iterate. This satisfactorily gives the partial reflections from a region of high velocity gradient, even when they are generated near the turning point of the ray. Although the results are somewhat complicated in the frequency domain, in the time domain all types of wave interaction are described by six universal time functions. For any problem, these functions are scaled in time according to the depth of the interaction, and in strength according to the magnitude of the coupling parameter. Numerical results and approximations are given for these functions. Coupling parameters are investigated for acoustic and elastic waves in a plane model, and acoustic and elastic-gravitational waves in a spherical model. The same universal time functions allow the excitation of elastic waves to be studied when the source is in a region of high velocity gradient or is near the wave's turning point. Results are given for a moment tensor, point source in plane and spherical models.  相似文献   

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

18.
Summary. As high-frequency elastic waves propagate through real media, it is common for caustics and focusing to occur. Typically, rays may envelop a caustic surface in space, or exceptionally they may all coalesce at a focal point. In strong motion seismology, the observed large fluctuations in peak acceleration and intensity of ground shaking may just be a consequence of focusing and caustics created by waves propagating through irregularly shaped sedimentary basins. These basins, acting as deformed optical lenses, are capable of producing a complex network of patches and seemingly isolated pockets of intensified damage or high intensity shaking where the caustic intersects the ground surface. We adapt methods from optics and catastrophe theory to study the properties of caustics induced by typical sedimentary basins. Several hypothetical examples are shown that reflect the fact that these properties are useful to assess quantitatively the degree of wavefield amplification to be expected. A good correlation is found between actual damage patterns and caustic locations computed for the Caracas, Venezuela earthquake of 1967.  相似文献   

19.
Summary. We propose a simplified method for the calculation of near field accelerograms. It is based upon the hypothesis that, in the course of dynamic faulting, the dominating part of the seismic radiation is emitted by the rupture front. As the rupture moves smoothly it radiates continuously, generating the low-frequency part of the field. High-frequency waves are produced by jumps in the rupture velocity and abrupt changes in the stress intensity factor. The wave-front discontinuities created in this fashion are evaluated by asymptotic methods and may be propagated away from the source by ray theoretical methods. We apply our technique to the evaluation of asymptotic near field accelerograms for a circular fault buried in a half-space. The agreement with numerical accelerograms calculated by full-wave theory is very satisfactory. Two problems are given particular emphasis: (1) the phase shifts introduced by focusing and (2) a simpler method, based on dislocation theory, is proposed for the calculation of the radiation coefficients from a discontinuously moving rupture front.  相似文献   

20.
We study properties of the energy-flux vector and other related energy quantities of homogeneous and inhomogeneous time-harmonic P and S plane waves, propagating in unbounded viscoelastic anisotropic media, both analytically and numerically. We propose an algorithm for the computation of the energy-flux vector, which can be used for media of unrestricted anisotropy and viscoelasticity, and for arbitrary homogeneous or inhomogeneous plane waves. Basic part of the algorithm is determination of the slowness vector of a homogeneous or inhomogeneous wave, which satisfies certain constraints following from the equation of motion. Approaches for determination of a slowness vector commonly used in viscoelastic isotropic media are usually difficult to use in viscoelastic anisotropic media. Sometimes they may even lead to non-physical solutions. To avoid these problems, we use the so-called mixed specification of the slowness vector, which requires, in a general case, solution of a complex-valued algebraic equation of the sixth degree. For simpler cases, as for SH waves propagating in symmetry planes, the algorithm yields simple analytic solutions. Once the slowness vector is known, determination of energy flux and of other energy quantities is easy. We present numerical examples illustrating the behaviour of the energy-flux vector and other energy quantities, for homogeneous and inhomogeneous plane P , SV and SH waves.  相似文献   

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

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