共查询到20条相似文献，搜索用时 15 毫秒
1.
Due to the lateral heterogeneity of the upper layers of the Earth, paths of surface waves deviate from arcs of great circles. Because of the sphericity of the Earth, the paths intersect on a hemisphere opposite to the epicenter and form caustics consisting of two branches, with their tangent point being a cusp. For this reason, the field of surface waves cannot be analyzed in terms of the ray theory at distances larger than 90°. The asymptotic approach to the analysis of the field in the vicinity of such caustics is very ill-suited for numerical implementation. The difficulties of such an approach to the field calculation are aggravated by the fact that such caustics are superimposed in some regions. Therefore, it is suggested to use the theorem of representation, according to which the field within a certain contour is expressed as an integral whose integrand contains values of the function itself, its derivative along the normal to the contour, and Green’s function. The field on the contour (the circle bounding a hemisphere centered at the epicenter) is calculated by the ray method because rays do not intersect on this hemisphere. These data are used for the construction of the field on the opposite hemisphere assumed to be homogeneous, which enables the construction of Green’s function for this hemisphere. This limitation is not very stringent because the configuration of rays and caustics on this hemisphere is mainly determined by the field on the circle. The integral in the representation theorem is calculated numerically. Numerical examples are presented for models in which one caustic or two superimposed caustics form. These calculations yield constraints on variations in the amplitude and phase of the wave. Rayleigh wave fields are also calculated for a model of the real Earth. It is shown that, at some points, the Rayleigh wave spectrum can be strongly distorted because caustics corresponding to different periods differ in shape. 相似文献
2.
We describe the behaviour of the anisotropic–ray–theory S–wave rays in a velocity model with a split intersection singularity. The anisotropic–ray–theory S–wave rays crossing the split intersection singularity are smoothly but very sharply bent. While the initial–value rays can be safely traced by solving Hamilton’s equations of rays, it is often impossible to determine the coefficients of the equations of geodesic deviation (paraxial ray equations, dynamic ray tracing equations) and to solve them numerically. As a result, we often know neither the matrix of geometrical spreading, nor the phase shift due to caustics. We demonstrate the abrupt changes of the geometrical spreading and wavefront curvature of the fast anisotropic–ray–theory S wave. We also demonstrate the formation of caustics and wavefront triplication of the slow anisotropic–ray–theory S wave.Since the actual S waves propagate approximately along the SH and SV reference rays in this velocity model, we compare the anisotropic–ray–theory S–wave rays with the SH and SV reference rays. Since the coupling ray theory is usually calculated along the anisotropic common S–wave rays, we also compare the anisotropic common S–wave rays with the SH and SV reference rays. 相似文献
3.
Asymptotic methods provide an efficient way to compute seismograms in heterogeneous media. However, zeroth-order ray theory, the simplest of the asymptotic methods, often fails because of the presence of caustics. Maslov theory is an extension of zeroth-order ray theory, which gives a uniformly valid expression of the wavefield everywhere, including the caustics. This result is given in terms of an integral of ray data over one or two ray parameters. It is shown in this paper how geometrical arrivals are constructed in the one and two-parameter Maslov integrals.In practice Maslov seismograms have been computed using only one ray parameter. However, in three-dimensional media two parameters are needed to uniquely define a ray. In this paper we present an efficient algorithm to compute two-parameter Maslov integrals. The Maslov integral is evaluated by computing the frequency-to-time Fourier transform prior to integration over the ray parameters. The wavefield is then discretized by smoothing with a boxcar function. The resulting expression, which only requires the results of ordinary kinematic and dynamic ray tracing, cen be computed efficiently and robustly. A numerical example is given that illustrates the use of this algorithm. 相似文献
4.
This paper is the second in a sequel of two papers and dedicated to the computation of paraxial rays and dynamic characteristics along the stationary rays obtained in the first paper. We start by formulating the linear, second‐order, Jacobi dynamic ray tracing equation. We then apply a similar finite‐element solver, as used for the kinematic ray tracing, to compute the dynamic characteristics between the source and any point along the ray. The dynamic characteristics in our study include the relative geometric spreading and the phase correction due to caustics (i.e. the amplitude and the phase of the asymptotic form of the Green's function for waves propagating in 3D heterogeneous general anisotropic elastic media). The basic solution of the Jacobi equation is a shift vector of a paraxial ray in the plane normal to the ray direction at each point along the central ray. A general paraxial ray is defined by a linear combination of up to four basic vector solutions, each corresponds to specific initial conditions related to the ray coordinates at the source. We define the four basic solutions with two pairs of initial condition sets: point–source and plane‐wave. For the proposed point–source ray coordinates and initial conditions, we derive the ray Jacobian and relate it to the relative geometric spreading for general anisotropy. Finally, we introduce a new dynamic parameter, similar to the endpoint complexity factor, presented in the first paper, used to define the measure of complexity of the propagated wave/ray phenomena. The new weighted propagation complexity accounts for the normalized relative geometric spreading not only at the receiver point, but along the whole stationary ray path. We propose a criterion based on this parameter as a qualifying factor associated with the given ray solution. To demonstrate the implementation of the proposed method, we use several isotropic and anisotropic benchmark models. For all the examples, we first compute the stationary ray paths, and then compute the geometric spreading and analyse these trajectories for possible caustics. Our primary aim is to emphasize the advantages, transparency and simplicity of the proposed approach. 相似文献
5.
《中国科学:地球科学(英文版)》2021,(6)
The perfectly matched layer(PML) boundary condition has been proven to be effective for attenuating reflections from model boundaries during wavefield simulation. As such, it has been widely used in time-domain finite-difference wavefield simulations. The conventional PML has poor performance for near grazing incident waves and low-frequency reflections. To overcome these limitations, a more complex frequency-shifted stretch(CSF) function is introduced, which is known as the CFSPML boundary condition and can be implemented in the time domain by a recursive convolution technique(CPML). When implementing the PML technique to second-order wave equations, all the existing methods involve adding auxiliary terms and rewriting the wave equations into new second-order partial differential equations that can be simulated by the finite-difference scheme, which may affect the efficiency of numerical simulation. In this paper, we propose a relatively simple and efficient approach to implement CPML for the second-order equation system, which solves the original wave equations numerically in the stretched coordinate. The spatial derivatives in the stretched coordinate are computed by adding a correction term to the regular derivatives. Once the first-order spatial derivatives are computed, we computed the second-order spatial derivatives in a similar way; therefore, we refer to the method as two-step CPML(TS-CPML). We apply the method to the second-order acoustic wave equation and a coupled second-order pseudo-acoustic TTI wave equation. Our simulations indicate that amplitudes of reflected waves are only about half of those computed with the traditional CPML method, suggesting that the proposed approach has computational advantages and therefore can be widely used for forwarding modeling and seismic imaging. 相似文献
6.
A new method of numerical computation of elastic wavefields in regions containing caustics is tested. The method is an extension of the asymptotic ray theory (ART). The essential features of the method consist of the application of expressions which are well defined at caustics and expressed in terms of ray tracing combined with complex ray tracing in caustic shadows. The method and an outline of the underlying theory are briefly presented, followed by a comparison with finite differences on a test model involving a caustic cusp. The comparison reveals the unexpectedly high degree of accuracy of the new method. 相似文献
7.
L. Klimeš 《Studia Geophysica et Geodaetica》2006,50(3):417-430
Explicit equations for the spatial derivatives and perturbation derivatives of amplitude in both isotropic and anisotropic
media are derived. The spatial and perturbation derivatives of the logarithm of amplitude can be calculated by numerical quadratures
along the rays.
The spatial derivatives of amplitude may be useful in calculating the higher-order terms in the ray series, in calculating
the higher-order amplitude coefficients of Gaussian beams, in estimating the accuracy of zero-order approximations of both
the ray method and Gaussian beams, in estimating the accuracy of the paraxial approximation of individual Gaussian beams,
or in estimating the accuracy of the asymptotic summation of paraxial Gaussian beams. The perturbation derivatives of amplitude
may be useful in perturbation expansions from elastic to viscoelastic media and in estimating the accuracy of the common-ray
approximations of the amplitude in the coupling ray theory. 相似文献
8.
Andrei M. Ionov 《Geophysical Prospecting》2013,61(1):104-119
The dynamic response of a semi‐infinite fluid‐filled borehole embedded in an elastic half‐space under a concentrated normal surface load is analysed in the long‐wavelength limit. The solution of the problem is obtained with integral transforms in the form of a double integral with respect to the slowness and frequency. The partial P‐ and SV‐wave responses are further transformed to path integrals along Cagniard paths in the complex slowness plane. Unlike the traditional Cagniard‐de Hoop technique based on the Laplace transform of time dependence, this paper is based on the Fourier transform. The tube‐wave response is presented as a causal integral over a slowness range. The resultant representation in the time‐domain is suitable for the numerical evaluation of the complete response in the fluid‐filled borehole, especially at large distances. Asymptotic analysis of seismic phases arising in the borehole is performed on the basis of the obtained solution. The complete asymptotic wavefield consists in P‐ and SV‐waves, the Rayleigh wave and the low‐frequency Stoneley (tube) wave. Pressure synthetics obtained by the use of the asymptotic formulas are shown to be in good agreement with straightforward calculations. 相似文献
9.
Paraxial ray methods have found broad applications in the seismic ray method and in numerical modelling and interpretation
of high-frequency seismic wave fields propagating in inhomogeneous, isotropic or anisotropic structures. The basic procedure
in paraxial ray methods consists in dynamic ray tracing. We derive the initial conditions for dynamic ray equations in Cartesian
coordinates, for rays initiated at three types of initial manifolds given in a three-dimensional medium: 1) curved surfaces
(surface source), 2) isolated points (point source), and 3) curved, planar and non-planar lines (line source). These initial
conditions are very general, valid for homogeneous or inhomogeneous, isotropic or anisotropic media, and for both a constant
and a variable initial travel time along the initial manifold. The results presented in the paper considerably extend the
possible applications of the paraxial ray method. 相似文献
10.
11.
-- A time-domain pure-state polarization analysis method is used to characterize surface waves traversing California parallel to the plate boundary. The method is applied to data recorded at four broadband stations in California from twenty-six large, shallow earthquakes which occurred since 1988, yielding polarization parameters such as the ellipticity, Euler angles, instantaneous periods, and wave incident azimuths. The earthquakes are located along the circum-Pacific margin and the ray paths cluster into two groups, with great-circle paths connecting stations MHC and PAS or CMB and GSC. The first path (MHC-PAS) is in the vicinity of the San Andreas Fault System (SAFS), and the second (CMB-GSC) traverses the Sierra Nevada Batholith parallel to and east of the SAFS. Both Rayleigh and Love wave data show refractions due to lateral velocity heterogeneities under the path, indicating that accurate phase velocity and attenuation analysis requires array measurements. T he Rayleigh waves are strongly affected by low velocity anomalies beneath Central California, with ray paths bending eastward as waves travel toward the south, while Love waves are less affected, providing observables to constrain the depth extent of anomalies. Strong lateral gradients in the lithospheric structure between the continent and the ocean are the likely cause of the path deflections. 相似文献
12.
矢量波场弹性波Kirchhoff偏移 总被引：2，自引：0，他引：2
Based on Kuo and Dai＇s vectorial wave-field extrapolation equations, we derive new Kirchhoff migration equations by introducing unit vectors which represent the ray directions at the imaging points of the reflected P- and PS converted-waves. Furthermore, using the slope of the events on shot records and a ray racing procedure, mirror-image reflection points are found and the reflection data are smeared along the Fresnel zone. The migration method proposed in this paper solves two troublesome imaging problems caused by limited receiving aperture and migration artifacts resulting from wave propagation at the velocities of non original wave type. The migration method is applied successfully with model data, demonstrating that the new method is effective and correct. 相似文献
13.
In downhole microseismic monitoring, accurate event location relies on the accuracy of the velocity model. The model can be estimated along with event locations. Anisotropic models are important to get accurate event locations. Taking anisotropy into account makes it possible to use additional data – two S-wave arrivals generated due to shear-wave splitting. However, anisotropic ray tracing requires iterative procedures for computing group velocities, which may become unstable around caustics. As a result, anisotropic kinematic inversion may become time consuming. In this paper, we explore the idea of using simplified ray tracing to locate events and estimate medium parameters. In the simplified ray-tracing algorithm, the group velocity is assumed to be equal to phase velocity in both magnitude and direction. This assumption makes the ray-tracing algorithm five times faster compared to ray tracing based on exact equations. We present a set of tests showing that given perforation-shot data, one can use inversion based on simplified ray-tracing even for moderate-to-strong anisotropic models. When there are no perforation shots, event-location errors may become too large for moderately anisotropic media. 相似文献
14.
The coupling ray theory is usually applied to anisotropic common reference rays, but it is more accurate if it is applied to reference rays which are closer to the actual wave paths. If we know that a medium is close to uniaxial (transversely isotropic), it may be advantageous to trace reference rays which resemble the SH–wave and SV–wave rays. This paper is devoted to defining and tracing these SH and SV reference rays of elastic S waves in a heterogeneous generally anisotropic medium which is approximately uniaxial (approximately transversely isotropic), and to the corresponding equations of geodesic deviation (dynamic ray tracing). All presented equations are simultaneously applicable to ordinary and extraordinary reference rays of electromagnetic waves in a generally bianisotropic medium which is approximately uniaxially anisotropic. The improvement of the coupling–ray–theory seismograms calculated along the proposed SH and SV reference rays, compared to the coupling–ray–theory seismograms calculated along the anisotropic common reference rays, has already been numerically demonstrated by the authors in four approximately uniaxial velocity models. 相似文献
15.
The standard ray theory (RT) for inhomogeneous anisotropic media does not work properly or even fails when applied to S-wave
propagation in inhomogeneous weakly anisotropic media or in the vicinity of shear-wave singularities. In both cases, the two
shear waves propagate with similar phase velocities. The coupling ray theory was proposed to avoid this problem. In it, amplitudes
of the two S waves are computed by solving two coupled, frequency-dependent differential equations along a common S-wave ray.
In this paper, we test the recently developed approximation of coupling ray theory (CRT) based on the common S-wave rays obtained
by first-order ray tracing (FORT). As a reference, we use the Fourier pseudospectral method (FM), which does not suffer from
the limitations of the ray method and yields very accurate results. We study the behaviour of shear waves in weakly anisotropic
media as well as in the vicinity of intersection, kiss or conical singularities. By comparing CRT and RT results with results
of the FM, we demonstrate the clear superiority of CRT over RT in the mentioned regions as well as the dangers of using RT
there. 相似文献
16.
N. T. Afanasiev A. A. Zheonykh M. K. Ivelskaya V. I. Sazhin M. V. Tinin V. E. Unuchkov 《Journal of Atmospheric and Solar》2001,63(18)
This paper is concerned with the oblique propagation of decametric radio waves in the ionosphere with random electron density irregularities. Effective parameters are introduced for calculating the influence of irregularities on the wave field structure. A technique is proposed for determining these parameters from measurements of statistical characteristics of the signal in the vicinity of a regular caustic. The technique uses asymptotic expressions obtained using the interference integral method and perturbation theory, as well as matching them to the numerical solution on the basis of the method of characteristics. A global semi-empirical model that is updated for current ionospheric conditions is used to specify the background medium. The proposed technique has been tested using data from a number of mid-latitude paths. Results obtained in this study testify that the technique deserves a practical implementation. 相似文献
17.
George D. Manolis Tsviatko V. Rangelov 《Soil Dynamics and Earthquake Engineering》2006,26(10):952-959
Geological media are invariably non-homogeneous, which complicates considerably the analysis of seismically induced wave propagation phenomena. Thus, closed-form solutions in the form of Green's functions are difficult to construct, but are quite valuable in their own right and often play the role of kernels in boundary integral equation formulations that are used for the solution of complex boundary-value problems of engineering importance. In this work, we examine in some detail the types of wave-like equations that result from vector decomposition of the equations of motion for the infinitely extending non-homogeneous continuum, which would be a first step for evaluating Green's functions. Specifically, an eigenvalue analysis is first performed, followed by computations using the finite difference method for a specific example involving a soil layer with quadratically varying material parameters. The aforementioned wave-like equations, defined in terms of dilatational and rotational strains, are originally coupled. Their uncoupling involves use of algebraic transformations, which are in turn valid for certain restricted categories of non-homogeneous materials. Numerical solution of these equations clearly shows attenuation patterns and phase changes that are manifested as the incoming wave disturbance is continuously scattered by non-constant material stiffness values encountered along the propagation path. 相似文献
18.
We present a new ray bending approach, referred to as the Eigenray method, for solving two‐point boundary‐value kinematic and dynamic ray tracing problems in 3D smooth heterogeneous general anisotropic elastic media. The proposed Eigenray method is aimed to provide reliable stationary ray path solutions and their dynamic characteristics, in cases where conventional initial‐value ray shooting methods, followed by numerical convergence techniques, become challenging. The kinematic ray bending solution corresponds to the vanishing first traveltime variation, leading to a stationary path between two fixed endpoints (Fermat's principle), and is governed by the nonlinear second‐order Euler–Lagrange equation. The solution is based on a finite‐element approach, applying the weak formulation that reduces the Euler–Lagrange second‐order ordinary differential equation to the first‐order weighted‐residual nonlinear algebraic equation set. For the kinematic finite‐element problem, the degrees of freedom are discretized nodal locations and directions along the ray trajectory, where the values between the nodes are accurately and naturally defined with the Hermite polynomial interpolation. The target function to be minimized includes two essential penalty (constraint) terms, related to the distribution of the nodes along the path and to the normalization of the ray direction. We distinguish between two target functions triggered by the two possible types of stationary rays: a minimum traveltime and a saddle‐point solution (due to caustics). The minimization process involves the computation of the global (all‐node) traveltime gradient vector and the traveltime Hessian matrix. The traveltime Hessian is used for the minimization process, analysing the type of the stationary ray, and for computing the geometric spreading of the entire resolved stationary ray path. The latter, however, is not a replacement for the dynamic ray tracing solution, since it does not deliver the geometric spreading for intermediate points along the ray, nor the analysis of caustics. Finally, we demonstrate the efficiency and accuracy of the proposed method along three canonical examples. 相似文献
19.
地震诱导电磁现象是国内外地学领域十分关注的前沿问题,前人对地震波和电磁场耦合波场的认识主要是基于规则模型获得的.为研究含起伏地表和地下界面的地层中震电波场激发、传播特性,本文采用有限元软件COMSOL Multiphysics模拟点震源激发的电磁场.首先给出频率域二维SHTE模式震电耦合方程组,然后利用COMSOL软件建立计算模型,并求解出点力源激发震电波场的频率域响应,最后利用FFT变换得到地震波场和电磁场的时间域波形.模拟结果表明,震电波场中存在三种类型的电磁信号,第一种是震源直接激发的电磁波;第二种是地震波在分界面处激发的电磁波(包括自由表面、地下不同介质分界面);第三种是伴随地震波的同震信号,前两种电磁波比地震波更早到达远处观测台站,对地震预警有重要意义.此外,研究还发现:当地震波传播至地表并沿着地表传播时,在地表附近空气层中同样记录到了伴随地震波传播的电磁扰动信号,该信号与相同水平源距条件下、地下观测点接收到的电磁信号相同,这与前人的一些观测结果相符.本文研究结果为今后地震电磁信号的解释提供了理论证据. 相似文献
20.
Diffraction and anelasticity problems involving decaying, “evanescent” or “inhomogeneous” waves can be studied and modelled using the notion of “complex rays”. The wavefront or “eikonal” equation for such waves is in general complex and leads to rays in complex position-slowness space. Initial conditions must be specified in that domain: for example, even for a wave originating in a perfectly elastic region, the ray to a real receiver in a neighbouring anelastic region generally departs from a complex point on the initial-values surface. Complex ray theory is the formal extension of the usual Hamilton equations to complex domains. Liouville's phase-space-incompressibility theorem and Fermat's stationary-time principle are formally unchanged. However, an infinity of paths exists between two fixed points in complex space all of which give the same final slowness, travel time, amplitude, etc. This does not contradict the fact that for a given receiver position there is a unique point on the initial-values surface from which this infinite complex ray family emanates.In perfectly elastic media complex rays are associated with, for example, evanescent waves in the shadow of a caustic. More generally, caustics in anelastic media may lie just outside the real coordinate subspace and one must trace complex rays around the complex caustic in order to obtain accurate waveforms nearby or the turning waves at greater distances into the lit region. The complex extension of the Maslov method for computing such waveforms is described. It uses the complex extension of the Legendre transformation and the extra freedom of complex rays makes pseudocaustics avoidable. There is no need to introduce a Maslov/KMAH index to account for caustics in the geometrical ray approximation, the complex amplitude being generally continuous. Other singular ray problems, such as the strong coupling around acoustic axes in anisotropic media, may also be addressed using complex rays.Complex rays are insightful and practical for simple models (e.g. homogeneous layers). For more complicated numerical work, though, it would be desirable to confine attention to real position coordinates. Furthermore, anelasticity implies dispersion so that complex rays are generally frequency dependent. The concept of group velocity as the velocity of a spatial or temporal maximum of a narrow-band wave packet does lead to real ray/Hamilton equations. However, envelope-maximum tracking does not itself yield enough information to compute synthetic seismogramsFor anelasticity which is weak in certain precise senses, one can set up a theory of real, dispersive wave-packet tracking suitable for synthetic seismogram calculations in linearly visco-elastic media. The seismologically-accepiable constant-Q rheology of Liu et al. (1976), for example, satisfies the requirements of this wave-packet theory, which is adapted from electromagnetics and presented as a reasonable physical and mathematical basis for ray modelling in inhomogeneous, anisotropic, anelastic media. Dispersion means that one may need to do more work than for elastic media. However, one can envisage perturbation analyses based on the ray theory presented here, as well as extensions like Maslov's which are based on the Hamiltonian properties. 相似文献