共查询到20条相似文献，搜索用时 343 毫秒
1.
In the computation of paraxial travel times and Gaussian beams, the basic role is played by the second-order derivatives of
the travel-time field at the reference ray. These derivatives can be determined by dynamic ray tracing (DRT) along the ray.
Two basic DRT systems have been broadly used in applications: the DRT system in Cartesian coordinates and the DRT system in
ray-centred coordinates. In this paper, the transformation relations between the second-order derivatives of the travel-time
field in Cartesian and ray-centred coordinates are derived. These transformation relations can be used both in isotropic and
anisotropic media, including computations of complex-valued travel times necessary for the evaluation of Gaussian beams. 相似文献
2.
The paper outlines the most important results of the paraxial complex geometrical optics (CGO) in respect to Gaussian beams
diffraction in the smooth inhomogeneous media and discusses interrelations between CGO and other asymptotic methods, which
reduce the problem of Gaussian beam diffraction to the solution of ordinary differential equations, namely: (i) Babich’s method,
which deals with the abridged parabolic equation and describes diffraction of the Gaussian beams; (ii) complex form of the
dynamic ray tracing method, which generalizes paraxial ray approximation on Gaussian beams and (iii) paraxial WKB approximation
by Pereverzev, which gives the results, quite close to those of Babich’s method. For Gaussian beams all the methods under
consideration lead to the similar ordinary differential equations, which are complex-valued nonlinear Riccati equation and
related system of complex-valued linear equations of paraxial ray approximation. It is pointed out that Babich’s method provides
diffraction substantiation both for the paraxial CGO and for complex-valued dynamic ray tracing method. It is emphasized also
that the latter two methods are conceptually equivalent to each other, operate with the equivalent equations and in fact are
twins, though they differ by names.
The paper illustrates abilities of the paraxial CGO method by two available analytical solutions: Gaussian beam diffraction
in the homogeneous and in the lens-like media, and by the numerical example: Gaussian beam reflection from a plane-layered
medium. 相似文献
3.
R. L. Nowack 《Pure and Applied Geophysics》2003,160(3-4):487-507
— In this paper, an overview of the calculation of synthetic seismograms using the Gaussian beam method is presented accompanied by some representative applications and new extensions of the method. Since caustics are a frequent occurrence in seismic wave propagation, modifications to ray theory are often necessary. In the Gaussian beam method, a summation of paraxial Gaussian beams is used to describe the propagation of high-frequency wave fields in smoothly varying inhomogeneous media. Since the beam components are always nonsingular, the method provides stable results over a range of beam parameters. The method has been shown, however, to perform better for some problems when different combinations of beam parameters are used. Nonetheless, with a better understanding of the method as well as new extensions, the summation of Gaussian beams will continue to be a useful tool for the modeling of high-frequency seismic waves in heterogeneous media. 相似文献
4.
5.
The 4 × 4 T -propagator matrix of a 3D central ray determines, among other important seismic quantities, second-order (parabolic or hyperbolic) two-point traveltime approximations of certain paraxial rays in the vicinity of the known central ray through a 3D medium consisting of inhomogeneous isotropic velocity layers. These rays result from perturbing the start and endpoints of the central ray on smoothly curved anterior and posterior surfaces. The perturbation of each ray endpoint is described only by a two-component vector. Here, we provide parabolic and hyperbolic paraxial two-point traveltime approximations using the T -propagator to feature a number of useful 3D seismic models, putting particular emphasis on expressing the traveltimes for paraxial primary reflected rays in terms of hyperbolic approximations. These are of use in solving several forward and inverse seismic problems. Our results simplify those in which the perturbation of the ray endpoints upon a curved interface is described by a three-component vector. In order to emphasize the importance of the hyperbolic expression, we show that the hyperbolic paraxial-ray traveltime (in terms of four independent variables) is exact for the case of a primary ray reflected from a planar dipping interface below a homogeneous velocity medium. 相似文献
6.
This study shows that the use of the first-order additional components of the ray method in the seismic wave field modeling is easy and that it can bring a substantial improvement of the standard ray results obtained with the zero-order ray approximation only. For the calculation of a first-order additional component, spatial derivatives of the parameters of the medium and spatial derivatives of the zero-order ray amplitude term are necessary. The evaluation of the former derivatives is straightforward; the latter derivatives can be calculated approximately from neighboring rays by substituting the derivatives by finite differences. This allows an effective calculation of the first-order additional terms in arbitrary laterally varying layered media.The importance of the first-order additional terms is demonstrated by the study of individual higher-order terms of the ray series representing elementaryP andS elastodynamic Green functions for a homogeneous isotropic medium. The study shows clearly that the consideration of the first-order additional terms leads to a more substantial decrease of the difference between approximate and exact elementary Green functions than any other higher-order term. With this in mind, effects of the first-order additional terms on the ray synthetic seismograms for aVSP configuration are studied. It is shown that the use of the additional terms leads to such phenomena, unknown in the zero-order approximation of the ray method, like quasi-elliptical and transverse polarization of a singleP wave or longitudinal polarization of a singleS wave. 相似文献
7.
A. B. Druzhinin 《Pure and Applied Geophysics》1996,148(3-4):637-683
A comprehensive approach, based on the general nonlinear ray perturbation theory (Druzhinin, 1991), is proposed for both a fast and accurate uniform asymptotic solution of forward and inverse kinematic problems in anisotropic media. It has been developed to modify the standard ray linearization procedures when they become inconsistent, by providing a predictable truncation error of ray perturbation series. The theoretical background consists in a set of recurrent expressions for the perturbations of all orders for calculating approximately the body wave phase and group velocities, polarization, travel times, ray trajectories, paraxial rays and also the slowness vectors or reflected/transmitted waves in terms of elastic tensor perturbations. We assume that any elastic medium can be used as an unperturbed medium. A total 2-D numerical testing of these expressions has been established within the transverse isotropy to verify the accuracy and convergence of perturbation series when the elastic constants are perturbed. Seismological applications to determine crack-induced anisotropy parameters on VSP travel times for the different wave types in homogeneous and horizontally layered, transversally isotropic and orthorhombic structures are also presented. A number of numerical tests shows that this method is in general stable with respect to the choice of the reference model and the errors in the input data. A proof of uniqueness is provided by an interactive analysis of the sensitivity functions, which are also used for choosing optimum source/receiver locations. Finally, software has been developed for a desktop computer and applied to interpreting specific real VSP observations as well as explaining the results of physical modelling for a 3-D crack model with the estimation of crack parameters. 相似文献
8.
Einar Iversen 《Pure and Applied Geophysics》1996,148(1-2):287-317
The motivation for this paper is to provide expressions for first-order partial derivatives of reflection point coordinates, taken with respect to model parameters. Such derivatives are expected to be useful for processes dealing with the problem of estimating velocities for depth migration of seismic data.The subject of the paper is a particular aspect of ray perturbation theory, where observed parameters—two-way reflection time and horizontal components of slowness, are constraining the ray path when parameters of the reference velocity model are perturbed. The methodology described here is applicable to general rays in a 3D isotropic, heterogeneous medium. Each ray is divided into a shot ray and a receiver ray, i.e., the ray portions between the shot/receiver and the reflection point, respectively. Furthermore, by freezing the initial horizontal slowness of these subrays as the model is perturbed,elementary perturbation quantities may be obtained, comprising derivatives of ray hit positions within theisochrone tangent plane, as well as corresponding time derivatives. The elementary quantities may be estimated numerically, by use of ray perturbation theory, or in some cases, analytically. In particular, when the layer above the reflection point is homogeneous, explicit formulas can be derived. When the elementary quantities are known,reflection point derivatives can be obtained efficiently from a set of linear expressions.The method is applicable for a common shot, receiver or offset data sorting. For these gather types, reflection point perturbationlaterally with respect to the isochrone is essentially different. However, in theperpendicular direction, a first-order perturbation is shown to beindependent of gather type.To evaluate the theory, reflection point derivatives were estimated analytically and numerically. I also compared first-order approximations to true reflection point curves, obtained by retracing rays for a number of model perturbations. The results are promising, especially with respect to applications in sensitivity analysis for prestack depth migration and in velocity model updating. 相似文献
9.
Paraxial ray methods for anisotropic inhomogeneous media 总被引：1，自引：0，他引：1
A new formalism of surface-to-surface paraxial matrices allows a very general and flexible formulation of the paraxial ray theory, equally valid in anisotropic and isotropic inhomogeneous layered media. The formalism is based on conventional dynamic ray tracing in Cartesian coordinates along a reference ray. At any user-selected pair of points of the reference ray, a pair of surfaces may be defined. These surfaces may be arbitrarily curved and oriented, and may represent structural interfaces, data recording surfaces, or merely formal surfaces. A newly obtained factorization of the interface propagator matrix allows to transform the conventional 6 × 6 propagator matrix in Cartesian coordinates into a 6 × 6 surface-to-surface paraxial matrix. This matrix defines the transformation of paraxial ray quantities from one surface to another. The redundant non-eikonal and ray-tangent solutions of the dynamic ray-tracing system in Cartesian coordinates can be easily eliminated from the 6 × 6 surface-to-surface paraxial matrix, and it can be reduced to 4 × 4 form. Both the 6 × 6 and 4 × 4 surface-to-surface paraxial matrices satisfy useful properties, particularly the symplecticity. In their 4 × 4 reduced form, they can be used to solve important boundary-value problems of a four-parametric system of paraxial rays, connecting the two surfaces, similarly as the well-known surface-to-surface matrices in isotropic media in ray-centred coordinates. Applications of such boundary-value problems include the two-point eikonal, relative geometrical spreading, Fresnel zones, the design of migration operators, and more. 相似文献
10.
Prestack depth migration is a key technology for imaging complex reservoirs in media with strong lateral velocity variations. Prestack migrations are broadly separated into ray-based and wave-equation-based methods. Because of its efficiency and flexibility, ray-based Kirchhoff migration is popular in the industry. However, it has difficulties in dealing with the multi-arrivals, caustics and shadow zones. On the other hand, wave-equation-based methods produce images superior to that of the ray-based methods, but they are expensive numerically, especially methods based on two-way propagators in imaging large regions. Therefore, reverse time migration algorithms with Gaussian beams have recently been proposed to reduce the cost, as they combine the high computational efficiency of Gaussian beam migration and the high accuracy of reverse time migration. However, this method was based on the assumption that the subsurface is isotropic. As the acquired azimuth and maximum offsets increase, taking into account the influence of anisotropy on seismic migration is becoming more and more crucial. Using anisotropic ray tracing systems in terms of phase velocity, we proposed an anisotropic reverse time migration using the Gaussian beams method. We consider the influence of anisotropy on the propagation direction and calculate the amplitude of Gaussian beams with optimized correlation coefficients in dynamic ray tracing, which simplifies the calculations and improves the applicability of the proposed method. Numerical tests on anisotropic models demonstrate the efficiency and accuracy of the proposed method, which can be used to image complex structures in the presence of anisotropy in the overburden. 相似文献
11.
在给出真实模型和相应光滑背景模型的情况下,如何计算扰动模型(散射体)产生的散射波场是一个有实际意义的正演问题.在Gabor变换域描述散射体,且入射波场为短时宽带信号时,散射波场可以在频率域用高斯束或时间域用高斯波包描述.相对于波动方程方法,高斯束和高斯波包的计算效率更高;背景模型光滑时,高斯束和高斯波包方法的精度也接近波动方程方法.文中导出了声波假设下应用高斯束和高斯波包计算散射波的方法.测试分析了高斯波包的计算精度.给出了一般散射体的散射波模拟策略.同时针对一个理论模型完成了本文方法计算散射波的实验,实验结果表明高斯波包散射波计算方法是有效可行的. 相似文献
12.
ADI plus interpolation: accurate finite-difference solution to 3D paraxial wave equation 总被引：2，自引：0，他引：2
Yanghua Wang 《Geophysical Prospecting》2001,49(5):547-556
An accurate finite-difference solution is developed for the paraxial wave equation in 3D seismic migration. The conventional alternating-direction-implicit (ADI) scheme used in migration causes errors, because the variables in the migration problem are complex-valued, not real-valued, and the imaginary part of the higher-order spatial derivatives cannot be ignored. The accuracy of the 3D paraxial extrapolator is preserved by (i) retaining these higher-order terms so that it does not produce the apparent azimuthal anisotropy in conventional migration, and (ii) filtering the non-physical evanescent waves during the downward extrapolation. The implementation of the accurate solution consists of two steps: firstly, the application of ADI to solve two tridiagonal systems sequentially, and secondly, an interpolation between the extrapolated wavefields of successive extrapolation levels. The method is computationally efficient as it uses the ADI scheme and, in addition, couples the correction for azimuthal anisotropy and the suppression of evanescent waves into a single operation, the interpolation step. 相似文献
13.
Luděk Klimeš 《Studia Geophysica et Geodaetica》2016,60(3):451-470
We consider the partial derivatives of travel time with respect to both spatial coordinates and perturbation parameters. These derivatives are very important in studying wave propagation and have already found various applications in smooth media without interfaces. In order to extend the applications to media composed of layers and blocks, we derive the explicit equations for transforming these travel–time derivatives of arbitrary orders at a general smooth curved interface between two arbitrary media. The equations are applicable to both real–valued and complex–valued travel time. The equations are expressed in terms of a general Hamiltonian function and are applicable to the transformation of travel–time derivatives in both isotropic and anisotropic media. The interface is specified by an implicit equation. No local coordinates are needed for the transformation. 相似文献
14.
高斯束逆时偏移是一种兼具计算效率和成像精度的深度域成像方法,能够面向目标成像.地下介质中黏滞性普遍存在,利用传统各向同性或完全弹性的成像方法处理黏滞性探区的数据会降低分辨率,并导致成像位置不准确和振幅欠估计等问题.本文在高斯束逆时偏移的基础上,通过对震源点和检波点处的波场进行衰减补偿,并结合高斯束求解时的角度信息,实现了黏声介质角度域高斯束逆时偏移方法.最后通过模型和实际资料试算对本文方法的正确性和适用性进行了验证.试算结果表明:相比于声波高斯束逆时偏移,本文方法能够对黏滞性引起的吸收衰减进行有效补偿,同时提取的角度域共成像点道集(ADCIGs)不仅可以用于分角度叠加成像压制成像噪声,而且能够为后续的偏移速度分析提供支撑. 相似文献
15.
16.
The partial derivatives of travel time with respect to model parameters are referred to as perturbations. Explicit equations for the second-order and higher-order perturbations of travel time in both isotropic and anisotropic media are derived. The perturbations of travel time and its spatial derivatives can be calculated by simple numerical quadratures along rays. 相似文献
17.
In this paper, ray theoretical amplitudes and travel times are calculated in slightly perturbed velocity models using perturbation analysis. Also, test inversions using travel time and amplitude are computed. The pertubation method is tested using a 3-D velocity model for NORSAR having velocity variations up to 8.0 percent. The perturbed amplitudes are found to be in excellent agreement with the calculated ray amplitudes. Velocity inversions based on travel time and amplitude are next investigated. Perturbation analysis using linearized ray equations is efficiently used to compute amplitude derivatives with respect to model parameters. The trial linearized inversions use smaller velocity variations of 1.7 percent to avoid possible effects due to ray shift, even though the perturbation analysis is valid for larger variations. The trial 2-D inversion results show that linearized amplitude inversions are complementary and not redundant to travel time inversions, even in smoothly varying models. 相似文献
18.
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. 相似文献
19.
各向异性射线理论基础上的局部角度域叠前深度偏移方法能够为深度域构造成像与基于角道集的层析反演提供有力支撑,但是对于复杂地质构造而言,高斯度叠前深度偏移在不失高效、灵活等特点的情况下,具有明显的精度优势.为此,本文研究局部角度域理论框架下的高斯束叠前深度偏移方法.为提高算法效率与实用性,文中讨论了一种从经典弹性参数表征的各向异性介质运动学和动力学射线方程演变而来的由相速度表征的简便形式,并提出了一种比较经济的各向异性高斯束近似合成方案.结合地震波局部角度域成像原理,讨论一种适合高斯束偏移的角度参数计算方法.国际上通用的理论模型合成数据试验表明：相比局部角度域Kirchhoff叠前深度偏移成像方法,本文方法具有更高的成像精度与抗噪能力,既适用于复杂构造成像,也可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎. 相似文献
20.
The moveout approximations play an important role in seismic data processing. The standard hyperbolic moveout approximation is based on an elliptical background model with two velocities: vertical and normal moveout. We propose a new set of moveout approximations based on a perturbation series in terms of anellipticity parameters using the alternative elliptical background model defined by vertical and horizontal velocities. We start with a transversely isotropic medium with a vertical symmetry axis. Then, we extend this approach to a homogeneous orthorhombic medium. To define the perturbation coefficients for a new background, we solve the eikonal equation with horizontal velocities in transversely isotropic medium with a vertical symmetry axis and orthorhombic media. To stabilise the perturbation series and improve the accuracy, the Shanks transform is applied for all the cases. We select different parameterisations for both velocities and anellipticity parameters for an orthorhombic model. From the comparison in traveltime error, the new moveout approximations result in better accuracy comparing with the standard perturbation‐based methods and other approximations. 相似文献