共查询到20条相似文献，搜索用时 46 毫秒
1.
The behaviour of the actual polarization of an electromagnetic wave or elastic S–wave is described by the coupling ray theory, which represents the generalization of both the zero–order isotropic and anisotropic ray theories and provides continuous transition between them. 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. In a generally anisotropic or bianisotropic medium, the actual wave paths may be approximated by the anisotropic–ray–theory rays if these rays behave reasonably. In an approximately uniaxial (approximately transversely isotropic) anisotropic medium, we can define and trace the SH (ordinary) and SV (extraordinary) reference rays, and use them as reference rays for the prevailing–frequency approximation of the coupling ray theory. In both cases, i.e. for the anisotropic–ray–theory rays or the SH and SV reference rays, we have two sets of reference rays. We thus obtain two arrivals along each reference ray of the first set and have to select the correct one. Analogously, we obtain two arrivals along each reference ray of the second set and have to select the correct one. In this paper, we suggest the way of selecting the correct arrivals. We then demonstrate the accuracy of the resulting prevailing–frequency approximation of the coupling ray theory using elastic S waves along the SH and SV reference rays in four different approximately uniaxial (approximately transversely isotropic) velocity models. 相似文献
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.
Luděk Klimeš 《Studia Geophysica et Geodaetica》2016,60(3):391-402
For a given stiffness tensor (tensor of elastic moduli) of a generally anisotropic medium, we estimate to what extent the medium is transversely isotropic (uniaxial) and determine the direction of its reference symmetry axis expressed in terms of the unit reference symmetry vector. If the medium is exactly transversely isotropic (exactly uniaxial), we obtain the direction of its symmetry axis. We can also calculate the first–order and second–order spatial derivatives of the reference symmetry vector which may be useful in tracing the reference rays for the coupling ray theory. The proposed method is tested using various transversely isotropic (uniaxial) and approximately transversely isotropic (approximately uniaxial) media. 相似文献
4.
The common ray approximation considerably simplifies the numerical algorithm of the coupling ray theory for S waves, but may introduce errors in travel times due to the perturbation from the common reference ray. These travel-time errors can deteriorate the coupling-ray-theory solution at high frequencies. It is thus of principal importance for numerical applications to estimate the errors due to the common ray approximation.We derive the equations for estimating the travel-time errors due to the isotropic and anisotropic common ray approximations of the coupling ray theory. These equations represent the main result of the paper. The derivation is based on the general equations for the second-order perturbations of travel time. The accuracy of the anisotropic common ray approximation can be studied along the isotropic common rays, without tracing the anisotropic common rays.The derived equations are numerically tested in three 1-D models of differing degree of anisotropy. The first-order and second-order perturbation expansions of travel time from the isotropic common rays to anisotropic-ray-theory rays are compared with the anisotropic-ray-theory travel times. The errors due to the isotropic common ray approximation and due to the anisotropic common ray approximation are estimated. In the numerical example, the errors of the anisotropic common ray approximation are considerably smaller than the errors of the isotropic common ray approximation.The effect of the isotropic common ray approximation on the coupling-ray-theory synthetic seismograms is demonstrated graphically. For comparison, the effects of the quasi-isotropic projection of the Green tensor, of the quasi-isotropic approximation of the Christoffel matrix, and of the quasi-isotropic perturbation of travel times on the coupling-ray-theory synthetic seismograms are also shown. The projection of the travel-time errors on the relative errors of the time-harmonic Green tensor is briefly presented. 相似文献
5.
Common-ray tracing and dynamic ray tracing for S waves in a smooth elastic anisotropic medium 总被引：1，自引：0，他引：1
L. Klimeš 《Studia Geophysica et Geodaetica》2006,50(3):449-461
Anisotropic common S-wave rays are traced using the averaged Hamiltonian of both S-wave polarizations. They represent very
practical reference rays for calculating S waves by means of the coupling ray theory. They eliminate problems with anisotropic-ray-theory
ray tracing through some S-wave slowness-surface singularities and also considerably simplify the numerical algorithm of the
coupling ray theory for S waves.
The equations required for anisotropic-common-ray tracing for S waves in a smooth elastic anisotropic medium, and for corresponding
dynamic ray tracing in Cartesian or ray-centred coordinates, are presented. The equations, for the most part generally known,
are summarized in a form which represents a complete algorithm suitable for coding and numerical applications. 相似文献
6.
7.
Determination of the ray vector (the unit vector specifying the direction of the group velocity vector) corresponding to a given wave normal (the unit vector parallel to the phase velocity vector or slowness vector) in an arbitrary anisotropic medium can be performed using the exact formula following from the ray tracing equations. The determination of the wave normal from the ray vector is, generally, a more complicated task, which is usually solved iteratively. We present a first-order perturbation formula for the approximate determination of the ray vector from a given wave normal and vice versa. The formula is applicable to qP as well as qS waves in directions, in which the waves can be dealt with separately (i.e. outside singular directions of qS waves). Performance of the approximate formulae is illustrated on models of transversely isotropic and orthorhombic symmetry. We show that the formula for the determination of the ray vector from the wave normal yields rather accurate results even for strong anisotropy. The formula for the determination of the wave normal from the ray vector works reasonably well in directions, in which the considered waves have convex slowness surfaces. Otherwise, it can yield, especially for stronger anisotropy, rather distorted results. 相似文献
8.
Luděk Klimeš 《Studia Geophysica et Geodaetica》2018,62(2):255-260
For a given stiffness tensor (tensor of elastic moduli) of a generally anisotropic medium, we can estimate the extent to which the medium is transversely isotropic, and determine the direction of its reference symmetry axis. In this paper, we rotate the given stiffness tensor about this reference symmetry axis, and determine the reference transversely isotropic (uniaxial) stiffness tensor as the average of the rotated stiffness tensor over all angles of rotation. The obtained reference transversely isotropic (uniaxial) stiffness tensor represents an analytically differentiable approximation of the given generally anisotropic stiffness tensor. The proposed analytic method is compared with a previous numerical method in two numerical examples. 相似文献
9.
TTI介质是石油地震勘探领域最常用的各向异性介质,快速计算TTI介质射线路径和走时信息有重要的研究意义.TTI介质传统运动学射线追踪方法一般基于任意弹性介质射线方程,利用Bond变换或者四阶张量变换来处理复杂的21个弹性参数,因而非常耗时.实际野外对称轴统一的TTI介质模型,一般可以看成VTI介质模型旋转一定角度获得.为此,本文推导了三维VTI介质射线追踪方程,提出先在本构坐标系中进行VTI介质射线追踪,再通过坐标旋转将射线路径旋转至观测坐标系中,获得TTI介质射线路径.数值模型计算表明该方法高效和精确,较传统方法效率提高了近4倍.在强各向异性等特殊情况下,体波波前面都与理论群速度面一致. 相似文献
10.
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. 相似文献
11.
The common-ray approximation eliminates problems with ray tracing through S-wave singularities and also considerably simplifies
the numerical algorithm of the coupling ray theory for S waves, but may introduce errors in travel times due to the perturbation
from the common reference ray. These travel-time errors can deteriorate the coupling-ray-theory solution at high frequencies.
It is thus of principal importance for numerical applications to estimate the errors due to the common-ray approximation applied.
The anisotropic-common-ray approximation of the coupling ray theory is more accurate than the isotropic-common-ray approximation.
We derive the equations for estimating the travel-time errors due to the anisotropic-common-ray (and also isotropic-common-ray)
approximation of the coupling ray theory. The errors of the common-ray approximations are calculated along the anisotropic
common rays in smooth velocity models without interfaces. The derivation is based on the general equations for the second-order
perturbations of travel time. 相似文献
12.
横向各向同性介质是地球内部广泛分布的一种各向异性介质.针对这种介质,我们对各向同性介质的最小走时树走时模拟方法进行了推广,推广后的方法可适用于非均匀、对称轴任意倾斜的横向各向同性介质模型.为保证计算效率,最小走时树的构建采用了一种子波传播区域随地震波传播动态变化的改进算法.对于弱各向异性介质,我们使用了一种新的地震波群速度近似表示方法,该方法基于用射线角近似表示相角的思想,对3种地震波（qP, qSV和qSH）均有较好的精度.应用本文地震波走时模拟方法对均匀介质、横向非均匀介质模型进行了计算,并将后者结果与弹性波方程有限元方法的模拟结果进行了对比,结果表明两者符合得很好.本文方法可用于横向各向同性介质的深度偏移及地震层析成像的深入研究. 相似文献
13.
In-seam seismic surveys with channel waves have been widely used in the United Kingdom and elsewhere to map coal-seams and to detect anomalous features such as dirt bands, seam thinning and thickening, and particularly in-seam faulting. Although the presence of cleat-induced anisotropy has been recognized in the past, almost all previous analyses have assumed homogeneous isotropic or transversely isotropic coal-seams. Channel waves, however, exhibit properties which cannot be fully explained without introducing anisotropy into the coal-seam. In particular, Love-type channel waves are observed for recording geometries where, in a homogeneous isotropic or transversely isotropic structure, the source would not be expected to excite transverse motion. Similarly, modes of channel-wave propagation display the coupled three-component motion of generalized modes in anisotropic substrates, which would not be expected for Rayleigh and Love wave motion in isotropy or in transversely isotropic media with azimuthal isotropy. We model the observed in-seam seismic channel waves with synthetic seismograms to gain an understanding of the effects of cleat-induced anisotropy on the behaviour of channel waves. The results show a reasonable good match with the observations in traveltime, relative amplitudes, dispersion characteristics and particle motions. We demonstrate that anisotropy in the surrounding country rocks contributes significantly to the coupling of channel wave particle motion, although its effect is not as strong as the anisotropy in the coal-seam. We conclude that the effects of cleat- and stress-induced anisotropy are observed and can be modelled with synthetic seismograms, and that anisotropy must be taken into account for the detailed interpretation of channel waves. 相似文献
14.
The objective is to provide, in one single paper, a complete collection of equations governing kinematic and dynamic ray tracing related to a symmetry plane of an anisotropic medium. Well known systems for kinematic ray tracing and in-plane dynamic ray tracing are reformulated for the purpose of clarity, by taking advantage of a vector representation of the Christoffel matrix elements and related quantities. A generalized formula is derived for the integrand in out-of-plane dynamic ray tracing, pertaining to a monoclinic medium. Integrands corresponding to non-tilted orthorhombic and transversely isotropic media are obtained as special cases. 相似文献
15.
本文拓展了一种模拟地震波在地球核幔边界D″区各向异性介质中传播的数值方法:谱元-简正振型耦合方法(CSEM).该方法通过在球对称各向同性介质空间采用简正振型方法,在各向异性的D″区采用谱元方法,并在两种介质的边界采用"DtN"算子耦合的策略计算一维模型PREM(见文献[1])或修改后的D″区横向各向同性VTI-PREM模型的理论地震图.模拟所得数值解结果与采用简正振型方法得到的解析解结果进行对比以验证方法的精度.在中国科学院地球深部结构重点实验室高性能计算机上使用128个CPU计算得到的结果显示,在10-5～0.125 Hz的频率范围内谱元简正振型法得到的波形与简正振型方法能很好拟合.此外,对于VTI介质结构模型,谱元简正振型法能够准确模拟S波分裂现象,从而验证了谱元简正振型耦合方法对各向异性介质中地震波传播数值模拟是一种有效的方法. 相似文献
16.
研究了三维弱各向异性近似下,利用伪P波(伪纵波)模拟弹性波场P分量在倾斜对称轴的横向各向同性(TTI)介质中的传播过程,并对比了分别基于弹性Hooke定律、弹性波投影和运动学色散方程所建立的三种二阶差分伪P波方程的正演特点.目前这些伪P波方程数值计算主要采用规则网格差分,但是规则网格在TTI模拟中有低效率、低精度以及不稳定的缺点.为了提高计算的精度,本文构建出相应方程的交错网格有限差分格式.通过对比伪P波方程在三维TTI介质中不同的数值模拟的表达形式,本文认为基于色散方程所建立的伪P波方程在模拟弹性波中P波传播的过程中具有最小的噪声.本文分析不同的各向同性对称轴空间角度的频散特征,并引入适当的横波速度维持计算的稳定.二维模型算例表明,本文提出的交错网格正演算法可以得到稳定光滑的伪P波正演波场.使用本文交错网格算法对二维BP TTI模型的逆时偏移也具有较稳定的偏移结果. 相似文献
17.
研究与实践表明,对于长偏移距、宽方位地震数据,忽略各向异性会明显降低成像质量,影响储层预测与描述的精度.针对典型的横向各向同性(TI)介质,本文面向深度域构造成像与偏移速度分析的需要,研究基于射线理论的局部角度域叠前深度偏移成像方法.它除了像传统Kirchhoff叠前深度偏移那样输出成像剖面和炮检距域的共成像点道集,还遵循地震波在成像点处的局部方向特征、基于扩展的脉冲响应叠加原理获得入射角度域和照明角度域的成像结果.为了方便快捷地实现TI介质射线走时与局部角度信息的计算,文中讨论和对比了两种改进的射线追踪方法:一种采用从经典各向异性介质射线方程演变而来的由相速度表征的简便形式;另一种采用由对称轴垂直的TI(即VTI)介质声学近似qP波波动方程推导出来的射线方程.文中通过坐标旋转将其扩展到了对称轴倾斜的TI(即TTI)介质.国际上通用的理论模型合成数据偏移试验表明,本文方法既适用于复杂构造成像,又可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎. 相似文献
18.
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. 相似文献
19.
Anisotropic reverse-time migration for tilted TI media 总被引：1，自引：0，他引：1
Seismic anisotropy in dipping shales results in imaging and positioning problems for underlying structures. We develop an anisotropic reverse‐time depth migration approach for P‐wave and SV‐wave seismic data in transversely isotropic (TI) media with a tilted axis of symmetry normal to bedding. Based on an accurate phase velocity formula and dispersion relationships for weak anisotropy, we derive the wave equation for P‐wave and SV‐wave propagation in tilted transversely isotropic (TTI) media. The accuracy of the P‐wave equation and the SV‐wave equation is analyzed and compared with other acoustic wave equations for TTI media. Using this analysis and the pseudo‐spectral method, we apply reverse‐time migration to numerical and physical‐model data. According to the comparison between the isotropic and anisotropic migration results, the anisotropic reverse‐time depth migration offers significant improvements in positioning and reflector continuity over those obtained using isotropic algorithms. 相似文献
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. 相似文献