共查询到20条相似文献，搜索用时 593 毫秒
1.
Traveltime computation by wavefront-orientated ray tracing 总被引：1，自引：0，他引：1
For multivalued traveltime computation on dense grids, we propose a wavefront‐orientated ray‐tracing (WRT) technique. At the source, we start with a few rays which are propagated stepwise through a smooth two‐dimensional (2D) velocity model. The ray field is examined at wavefronts and a new ray might be inserted between two adjacent rays if one of the following criteria is satisfied: (1) the distance between the two rays is larger than a predefined threshold; (2) the difference in wavefront curvature between the rays is larger than a predefined threshold; (3) the adjacent rays intersect. The last two criteria may lead to oversampling by rays in caustic regions. To avoid this oversampling, we do not insert a ray if the distance between adjacent rays is smaller than a predefined threshold. We insert the new ray by tracing it from the source. This approach leads to an improved accuracy compared with the insertion of a new ray by interpolation, which is the method usually applied in wavefront construction. The traveltimes computed along the rays are used for the estimation of traveltimes on a rectangular grid. This estimation is carried out within a region bounded by adjacent wavefronts and rays. As for the insertion criterion, we consider the wavefront curvature and extrapolate the traveltimes, up to the second order, from the intersection points between rays and wavefronts to a gridpoint. The extrapolated values are weighted with respect to the distances to wavefronts and rays. Because dynamic ray tracing is not applied, we approximate the wavefront curvature at a given point using the slowness vector at this point and an adjacent point on the same wavefront. The efficiency of the WRT technique is strongly dependent on the input parameters which control the wavefront and ray densities. On the basis of traveltimes computed in a smoothed Marmousi model, we analyse these dependences and suggest some rules for a correct choice of input parameters. With suitable input parameters, the WRT technique allows an accurate traveltime computation using a small number of rays and wavefronts. 相似文献
2.
A new 3D wavefield modelling approach based on dynamic ray tracing is presented. This approach is called wavefront construction, and it can be used in 3D models with constant or smoothly varying material properties (S- and P-velocity and density) separated by smooth interfaces. Wavefronts consisting of rays arranged in a triangular network are propagated stepwise through the model. At each time step, the differences in a number of parameters are checked between each pair of rays on the wavefront. New rays are interpolated whenever this difference between pairs of rays exceeds some predefined maximum value. A controlled sampling of the wavefront at all time steps is thus obtained. Receivers are given multiple-event values by interpolation when the wavefronts pass them. The strength of the wavefront construction method is that it is robust and efficient. 相似文献
3.
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. 相似文献
4.
Seismic amplitude variations with offset contain information about the elastic parameters. Prestack amplitude analysis seeks to extract this information by using the variations of the reflection coefficients as functions of angle of incidence. Normally, an approximate formula is used for the reflection coefficients, and variations with offset of the geometrical spreading and the anelastic attenuation are often ignored. Using angle of incidence as the dependent variable is also computationally inefficient since the data are recorded as a function of offset. Improved approximations have been derived for the elastic reflection and transmission coefficients, the geometrical spreading and the complex travel-time (including anelastic attenuation). For a 1 D medium, these approximations are combined to produce seismic reflection amplitudes (P-wave, S-wave or converted wave) as a Taylor series in the offset coordinate. The coefficients of the Taylor series are computed directly from the parameters of the medium, without using the ray parameter. For primary reflected P-waves, dynamic ray tracing has been used to compute the offset variations of the transmission coefficients, the reflection coefficient, the geometrical spreading and the anelastic attenuation. The offset variation of the transmission factor is small, while the variations in the geometrical spreading, absorption and reflection coefficient are all significant. The new approximations have been used for seismic modelling without ray tracing. The amplitude was approximated by a fourth-order polynomial in offset, the traveltime by the normal square-root approximation and the absorption factor by a similar expression. This approximate modelling was compared to dynamic ray tracing, and the results are the same for zero offset and very close for offsets less than the reflector depth. 相似文献
5.
We describe two practicable approaches for an efficient computation of seismic traveltimes and amplitudes. The first approach is based on a combined finite‐difference solution of the eikonal equation and the transport equation (the ‘FD approach’). These equations are formulated as hyperbolic conservation laws; the eikonal equation is solved numerically by a third‐order ENO–Godunov scheme for the traveltimes whereas the transport equation is solved by a first‐order upwind scheme for the amplitudes. The schemes are implemented in 2D using polar coordinates. The results are first‐arrival traveltimes and the corresponding amplitudes. The second approach uses ray tracing (the ‘ray approach’) and employs a wavefront construction (WFC) method to calculate the traveltimes. Geometrical spreading factors are then computed from these traveltimes via the ray propagator without the need for dynamic ray tracing or numerical differentiation. With this procedure it is also possible to obtain multivalued traveltimes and the corresponding geometrical spreading factors. Both methods are compared using the Marmousi model. The results show that the FD eikonal traveltimes are highly accurate and perfectly match the WFC traveltimes. The resulting FD amplitudes are smooth and consistent with the geometrical spreading factors obtained from the ray approach. Hence, both approaches can be used for fast and reliable computation of seismic first‐arrival traveltimes and amplitudes in complex models. In addition, the capabilities of the ray approach for computing traveltimes and spreading factors of later arrivals are demonstrated with the help of the Shell benchmark model. 相似文献
6.
Generalized non‐hyperbolic approximation for qP‐wave relative geometrical spreading in a layered transversely isotropic medium with a vertical symmetry axis 
下载免费PDF全文
下载免费PDF全文 Compensation for geometrical spreading along the ray‐path is important in amplitude variation with offset analysis especially for not strongly attenuative media since it contributes to the seismic amplitude preservation. The P‐wave geometrical spreading factor is described by a non‐hyperbolic moveout approximation using the traveltime parameters that can be estimated from the velocity analysis. We extend the P‐wave relative geometrical spreading approximation from the rational form to the generalized non‐hyperbolic form in a transversely isotropic medium with a vertical symmetry axis. The acoustic approximation is used to reduce the number of parameters. The proposed generalized non‐hyperbolic approximation is developed with parameters defined by two rays: vertical and a reference rays. For numerical examples, we consider two choices for parameter selection by using two specific orientations for reference ray. We observe from the numerical tests that the proposed generalized non‐hyperbolic approximation gives more accurate results in both homogeneous and multi‐layered models than the rational counterpart. 相似文献
7.
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. 相似文献
8.
We present a new method of three-dimensional (3-D) seismic ray tracing, based on an improvement to the linear traveltime interpolation (LTI) ray tracing algorithm. This new technique involves two separate steps. The first involves a forward calculation based on the LTI method and the dynamic successive partitioning scheme, which is applied to calculate traveltimes on cell boundaries and assumes a wavefront that expands from the source to all grid nodes in the computational domain. We locate several dynamic successive partition points on a cell's surface, the traveltimes of which can be calculated by linear interpolation between the vertices of the cell's boundary. The second is a backward step that uses Fermat's principle and the fact that the ray path is always perpendicular to the wavefront and follows the negative traveltime gradient. In this process, the first-arriving ray path can be traced from the receiver to the source along the negative traveltime gradient, which can be calculated by reconstructing the continuous traveltime field with cubic B-spline interpolation. This new 3-D ray tracing method is compared with the LTI method and the shortest path method (SPM) through a number of numerical experiments. These comparisons show obvious improvements to computed traveltimes and ray paths, both in precision and computational efficiency. 相似文献
9.
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. 相似文献
10.
针对线性走时插值算法(LTI)不能正确追踪逆向传播射线的问题, 目前已提出多种改进算法, 如扩张收缩LTI算法、 循环计算LTI算法、 动态网络最短路径射线追踪算法等, 但这些算法的计算效率普遍偏低. 在分析各种改进LTI算法的优劣后, 本文提出了改进动态网络最短路径射线追踪算法. 该改进算法依据波的传播规律以及LTI算法的基本方程, 排除动态网络最短路径射线追踪算法中大量冗余节点计算, 并采用传统的二叉树堆排序算法对波前阵列节点进行管理. 数值算例表明, 本文提出的改进算法具有较高的计算效率, 其计算效率是动态网络最短路径射线追踪算法的4.5—30倍, 是原始LTI算法的2—6.5倍; 当动态网络最短路径射线追踪算法采用堆排序算法时, 改进算法的计算效率是其3.5—15倍. 相似文献
11.
Summary A new approximation of the velocity-depth distribution in a vertically inhomogeneous medium is suggested. This approximation guarantees the continuity of velocity and of its first and second derivatives and does not generate false low-velocity zones. It is very suitable for the computations of seismic wave fields in vertically inhomogeneous media by ray methods and its modifications, as it removes many false anomalies from the travel-time and amplitude-distance curves of seismic body waves. The ray integrals can be evaluated in a closed form; the resulting formulae for rays, travel times and geometrical spreading are very simple. They do not contain any transcendental functions (such asln (x) orsin
–1, (x)) like other approximations; only the evaluation of one square root and of certain simple arithmetic expressions for each layer is required. From a computational point of view, the evaluation of ray integrals and of geometrical spreading is only slightly slower than for a system of homogeneous parallel layers and even faster than for a piece-wise linear approximation. 相似文献
12.
在许多地震反演和偏移成像方法中,都要涉及到射线路径和旅行时的计算.本文将波前面三角形网格剖分和三维波前重建法射线追踪技术结合使用,实现了射线路径和旅行时的准确快速计算.三维波前重建法射线追踪过程中可以保证稳定合理的射线密度,克服了常规射线追踪方法存在阴影区的问题.波前面三角形网格剖分在描述和拆分波前面时更加准确有效,而且不需太多的网格数目,从而提高了射线追踪的精度和效率.该方法在三维复杂构造成像方面有独特的优势,目前在实际的Kirchhoff 偏移中的已经有相关应用. 相似文献
13.
Emil Blias 《Geophysical Prospecting》2013,61(3):574-581
I introduce a new explicit form of vertical seismic profile (VSP) traveltime approximation for a 2D model with non‐horizontal boundaries and anisotropic layers. The goal of the new approximation is to dramatically decrease the cost of time calculations by reducing the number of calculated rays in a complex multi‐layered anisotropic model for VSP walkaway data with many sources. This traveltime approximation extends the generalized moveout approximation proposed by Fomel and Stovas. The new equation is designed for borehole seismic geometry where the receivers are placed in a well while the sources are on the surface. For this, the time‐offset function is presented as a sum of odd and even functions. Coefficients in this approximation are determined by calculating the traveltime and its first‐ and second‐order derivatives at five specific rays. Once these coefficients are determined, the traveltimes at other rays are calculated by this approximation. Testing this new approximation on a 2D anisotropic model with dipping boundaries shows its very high accuracy for offsets three times the reflector depths. The new approximation can be used for 2D anisotropic models with tilted symmetry axes for practical VSP geometry calculations. The new explicit approximation eliminates the need of massive ray tracing in a complicated velocity model for multi‐source VSP surveys. This method is designed not for NMO correction but for replacing conventional ray tracing for time calculations. 相似文献
14.
本文采用随机介质描述地球内部大尺度背景场上存在的小尺度不均匀性和自组织结构;文中没有采用传统的层状网格结构,提出了块状结构来描述复杂的二维自组织结构的方法,分别以高斯型、指数型和von Karman型等自相关函数描述各向同性和各向异性非均匀分布自组织特征;采用射线追踪分析了不同分布特征自组织结构对地震波运动学和动力学特征参数的影响;结果表明,由于地球内部介质的自组织性存在,地震波射线轨迹可能发生明显的畸变; 不同偏移距处,地震反射振幅减弱或增强; 自组织结构从高斯型到von Karman型,在小尺度上表现更大的非均匀性,因此走时和振幅表征依次更强的平均效应. 相似文献
15.
16.
Two-dimensional seismic processing is successful in media with little structural and velocity variation in the direction perpendicular to the plane defined by the acquisition direction and the vertical axis. If the subsurface is anisotropic, an additional limitation is that this plane is a plane of symmetry. Kinematic ray propagation can be considered as a two-dimensional process in this type of medium. However, two-dimensional processing in a true-amplitude sense requires out-of-plane amplitude corrections in addition to compensation for in-plane amplitude variation. We provide formulae for the out-of-plane geometrical spreading for P- and S-waves in transversely isotropic and orthorhombic media. These are extensions of well-known isotropic formulae.
For isotropic and transversely isotropic media, the ray propagation is independent of the azimuthal angle. The azimuthal direction is defined with respect to a possibly tilted axis of symmetry. The out-of-plane spreading correction can then be calculated by integrating quantities which describe in-plane kinematics along in-plane rays. If, in addition, the medium varies only along the vertical direction and has a vertical axis of symmetry, no ray tracing need be carried out. All quantities affecting the out-of-plane geometrical spreading can be derived from traveltime information available at the observation surface.
Orthorhombic media possess no rotational symmetry and the out-of-plane geometrical spreading includes parameters which, even in principle, are not invertible from in-plane experiments. The exact and approximate formulae derived for P- and S-waves are nevertheless useful for modelling purposes. 相似文献
For isotropic and transversely isotropic media, the ray propagation is independent of the azimuthal angle. The azimuthal direction is defined with respect to a possibly tilted axis of symmetry. The out-of-plane spreading correction can then be calculated by integrating quantities which describe in-plane kinematics along in-plane rays. If, in addition, the medium varies only along the vertical direction and has a vertical axis of symmetry, no ray tracing need be carried out. All quantities affecting the out-of-plane geometrical spreading can be derived from traveltime information available at the observation surface.
Orthorhombic media possess no rotational symmetry and the out-of-plane geometrical spreading includes parameters which, even in principle, are not invertible from in-plane experiments. The exact and approximate formulae derived for P- and S-waves are nevertheless useful for modelling purposes. 相似文献
17.
The recursive nature of rays in blocky models can be exploited to solve some difficult problems in seismic modelling. Each segment of a ray travels from an initial point up to a reflecting interface, where it is split into reflected and transmitted ray segments, which each continue in a similar way. The tree structure that thus emanates is conveniently handled by a recursive scheme. Recursion allows an automatic generation of all phases on a seismogram, together with all information necessary to analyse or select them. By operating recursively with a ray cell, bounded by a pair of vicinal rays in 2D, or a triplet of vicinal rays in 3D, and two successive isochrons, the two-point ray-tracing problem is reduced to a simple interpolation. Also, the cellular approach allows for a stable and robust evaluation of dynamic ray quantities without any paraxial tracing, which is cumbersome in blocky models of realistic complexity. Geometric shadows are filled by recursively generated diffractions. The recursive ray tracer has found applications in the fast computation of Green's functions in target-oriented inversion and in phase identification in VSP. 相似文献
18.
Ｉｎｔｒｏｄｕｃｔｉｏｎ３－Ｄｓｅｉｓｍｉｃｔｏｍｏｇｒａｐｈｙｈａｓｂｅｅｎａｐｐｌｉｅｄｔｏｖａｒｉｏｕｓｇｅｏｐｈｙｓｉｃａｌｐｒｏｂｌｅｍｓ．ＡｋｉａｎｄＬｅｅ（１９７６）ａｎｄＨａｗｌｅｙｅｔａｌ．（１９８１）ｉｎｖｅｒｔｅｄ３－Ｄｍｏｄｅ... 相似文献
19.
KLAUS HELBIG 《Geophysical Prospecting》1990,38(2):189-220
Wavefront charts in anisotropic gradient media are a useful tool in ray geometric constructions, particular in shear-wave exploration. They can be constructed by: (i) a family of wavefronts that contains a vertical plane as member - it is convenient to choose constant time increments; (ii) tracing one ray that makes everywhere the angle with the normal to the wavefront that is required by the anisotropy of the medium; (iii) scaling this ray to obtain a set of rays with different ray parameters; (iv) shifting these rays (with wavefront elements attached) so that they pass through a common source point; (v) interpolating the wavefronts between the elements. The construction is particularly simple in linear-gradient media, since here all members of the family of wavefronts are planes. Since the ray makes everywhere the angle prescribed by the anisotropy with the normal of the (plane) wavefronts, the ray has the shape of the slowness curve rotated by ?π/2. For isotropic media the slowness curve is a circle, and thus rays are circular arcs. The circles themselves intersect in the source point and in a second point above the surface of the earth. This provides a simple proof that wavefronts emanating from a point source in an isotropic linear-gradient medium are spheres: inversion of the set of circular rays with the source as centre maps the pencil of circular rays into a pencil of straight lines passing through a point. A pencil of concentric spheres around this point is perpendicular to the pencil of straight lines. On inverting back the pencil of spheres is mapped into another pencil of spheres that is perpendicular to the circular rays. 相似文献
20.
Seismic ray tracing in layered media becomes complicated and demanding when modeling for multiple ray codes (reflection/transmission
sequences) and/or dense acquisition geometries. However, we observe some redundancies in current algorithms: (a) the same
layers are crossed repeatedly by similar ray segments, and (b) the effort of tracing through a layer is determined by variations
in the incoming wavefront rather than the medium. We deal with these redundancies by separating the modeling process in two
stages: (Stage 1) compute ray field maps representing all ray segments between each pair of adjacent interfaces, then (Stage
2) for each desired ray code assemble the complete ray field from ray segments by iterative lookup in the ray field maps. 相似文献