共查询到20条相似文献，搜索用时 31 毫秒
1.
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. 相似文献
2.
Vlastislav Červený 《Studia Geophysica et Geodaetica》2013,57(2):267-275
Recently, several expressions for the two-point paraxial travel time in laterally varying, isotropic or anisotropic layered media were derived. The two-point paraxial travel time gives the travel time from point S′ to point R′, both these points being situated close to a known reference ray Ω, along which the ray-propagator matrix was calculated by dynamic ray tracing. The reference ray and the position of points S′ and R′ are specified in Cartesian coordinates. Two such expressions for the two-point paraxial travel time play an important role. The first is based on the 4 × 4 ray propagator matrix, computed by dynamic ray tracing along the reference ray in ray-centred coordinates. The second requires the knowledge of the 6 × 6 ray propagator matrix computed by dynamic ray tracing along the reference ray in Cartesian coordinates. Both expressions were derived fully independently, using different methods, and are expressed in quite different forms. In this paper we prove that the two expressions are fully equivalent and can be transformed into each other. 相似文献
3.
Summary Among the various dynamic ray tracing systems described by ervený and Hron[1] is one particular linear system of second order that readily provides identical parameters (in the ray centred coordinate system) to those that fall out of the system by Popov and Peník[3, 4]. Hence the initial conditions of the latter system for sources and interfaces can easily be used to provide those for the linear system of second order. 相似文献
4.
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. 相似文献
5.
The coupling ray theory bridges the gap between the isotropic and anisotropic ray theories, and is considerably more accurate than the anisotropic ray theory. The coupling ray theory is often approximated by various quasi-isotropic approximations.Commonly used quasi-isotropic approximations of the coupling ray theory are discussed. The exact analytical solution for the plane S wave, propagating along the axis of spirality in the 1-D anisotropic oblique twisted crystal model, is then numerically compared with the coupling ray theory and its three quasi-isotropic approximations. The three quasi-isotropic approximations of the coupling ray theory are (a) the quasi-isotropic projection of the Green tensor, (b) the quasi-isotropic approximation of the Christoffel matrix, (c) the quasi-isotropic perturbation of travel times. The comparison is carried out numerically in the frequency domain, comparing the exact analytical solution with the results of the 3-D ray tracing and coupling ray theory software. In the oblique twisted crystal model, the three studied quasi-isotropic approximations considerably increase the error of the coupling ray theory. Since these three quasi-isotropic approximations do not noticeably simplify the numerical implementation of the coupling ray theory, they should deffinitely be avoided. The common ray approximations of the coupling ray theory do not affect the plane wave, propagating along the axis of spirality in the 1-D oblique twisted crystal model, and should be studied in more complex models. 相似文献
6.
I reformulate well-established systems for kinematic and dynamic ray tracing in 3D heterogeneous media with arbitrary anisotropy. Matrices of size 3 × 3, e.g. the Christoffel matrix, are substituted by six-component vectors, and the Christoffel matrix elements are expressed explicitly in terms of the elements of the 6 × 6 matrix of elastic coefficients, written in Voigt notation. Thereby, I find it easier to see the effects on the ray tracing systems of vanishing elastic coefficients and Christoffel matrix elements and of vanishing derivatives with respect to spatial coordinates and slowness vector components. The eigenvalue of the current wave and its derivatives with respect to the ray parameters are included explicitly, which may be favorable for ray tracing processes optimized with respect to speed. With the ANRAY program as a reference, I show that the new formulation requires less optimization than the conventional one. 相似文献
7.
Gjøystdal Håvar Iversen Einar Laurain Renaud Lecomte Isabelle Vinje Vetle Åstebøl Ketil 《Studia Geophysica et Geodaetica》2002,46(2):113-164
Throughout the last twenty years, 3D seismic ray modelling has developed from a research tool to a more operational tool that has gained growing interest in the petroleum industry. Various areas of application have been established and new ones are under development. Many of these applications require a modelling system with flexible, robust and efficient modelling algorithms in the core. The present paper reviews the basic elements of such a system, based on the open model concept and the wavefront construction technique. In the latter, ervený's dynamic ray tracing is an intrinsic part. The modelling system can be used for generating ray attributes and synthetic seismograms for realistic 3D surveys with tens of thousands of shots and receivers. Moreover, some other types of application areas are illustrated: Production of Green's functions for prestack depth migration and hybrid modelling (combined ray and finite-difference modelling), attribute mapping and illumination analysis, both for survey planning and interpretation. Finally, the concepts of isochron rays and velocity rays related to seismic isochrons have been introduced recently, with very interesting future applications. 相似文献
8.
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. 相似文献
9.
The insufficiently known phenomenon of MF-wave propagation from the Earths surface through the magnetosphere (guiding) to the conjugate hemisphere and back to the transmitter has been experimentally studied. Computer modelling fulfilled on the basis of ray tracing showed that guiding was possible only from area of the main ionospheric trough. The effect of MF guiding is most useful for the diagnostics of the plasmapause, poleward edge of the trough, the diffuse precipitation boundary and so on. 相似文献
10.
电流线追踪电位电阻率层析成像方法初探 总被引：4，自引：3，他引：4
电阻率层析成像技术尽管已有了一些比较好的结果，但从国内外发表的文章可以看出，基本上采用的都是有限元方法，而电阻率层析成像的核心问题也就是雅可比矩阵的求取问题。有限元方法能够很好地实现该问题的求解，但需要的计算机内存及计算时间相当的大，为此我们类比地震学中走时射线追踪技术，开展了电流线追踪电位电阻率层析成像方法研究。 相似文献
11.
Analytical expressions for the exact 2 × 2 one-way propagator matrix of a plane S wave, propagating along the axis of spirality in the simple 1-D anisotropic simplified twisted crystal model, are presented. The analytical equations are useful in testing the applicability and accuracy of various approximate wavefield modelling methods, especially of the coupling ray theory and of its various quasi-isotropic approximations and various numerical implementations.In addition to the exact analytical solution of the elastodynamic equation in the simplified twisted crystal model, the analytical solutions of the equations of the four ray methods are given. The ray methods are (a) the coupling ray theory, (b) the coupling ray theory with the quasi-isotropic perturbation of travel times, (c) the anisotropic ray theory, (d) the isotropic ray theory. These four approximate solutions of the elastodynamic equation are roughly compared with the exact solution. Both the exact analytical solution and the analytical ray-theory solutions in the simplified twisted crystal model are also helpful in debugging computer codes for various approximate wavefield modelling methods, especially for the coupling ray theory. 相似文献
12.
A review of the 6 × 6 anisotropic interface ray propagator matrix in Cartesian coordinates and within the framework of the Hamiltonian formalism shows that there is one unique propagator satisfying the symplectic property. This is essential, since the symplecticity furnishes an exact inverse, while an eigenvalue analysis indicates that the propagator may be arbitrarily ill-conditioned. As such, the symplectic interface propagator naturally connects to symplectic ray integration algorithms for smooth media, designed to maintain accuracy. Moreover, several ray invariants for smooth media remain invariant across interfaces. It is straightforward to derive expressions for the interface propagator, both explicit and implicit. Symplecticity is equivalent to the condition that the propagator preserves the eikonal constraint across the interface. The symplectic interface propagator complies with phase matching of the incident and reflected/transmitted ray field, and is therefore in accordance with the earlier derived 4 × 4 matrix in ray-centred coordinates. The symplectic property is related to the symmetry of the second derivative matrix of the reflected/transmitted traveltime field. Thanks to the analytic expression of the symplectic interface propagator, relating interface curvature directly to second derivatives of traveltimes observed at a datum level, numerous applications are available in the area of processing and inversion. 相似文献
13.
针对线性走时插值算法(LTI)不能正确追踪逆向传播射线的问题, 目前已提出多种改进算法, 如扩张收缩LTI算法、 循环计算LTI算法、 动态网络最短路径射线追踪算法等, 但这些算法的计算效率普遍偏低. 在分析各种改进LTI算法的优劣后, 本文提出了改进动态网络最短路径射线追踪算法. 该改进算法依据波的传播规律以及LTI算法的基本方程, 排除动态网络最短路径射线追踪算法中大量冗余节点计算, 并采用传统的二叉树堆排序算法对波前阵列节点进行管理. 数值算例表明, 本文提出的改进算法具有较高的计算效率, 其计算效率是动态网络最短路径射线追踪算法的4.5—30倍, 是原始LTI算法的2—6.5倍; 当动态网络最短路径射线追踪算法采用堆排序算法时, 改进算法的计算效率是其3.5—15倍. 相似文献
14.
The period of interplanetary, geomagnetic and solar disturbances of September 7–15, 2005, is characterized by two sharp increases of solar wind velocity to 1000 km/s and great Dst variation of the geomagnetic field (~140 nT). The time variations of theoretical and experimental geomagnetic thresholds observed during this strong geomagnetic storm, their connection with solar wind parameters and the Dst index, and the features of latitudinal behavior of geomagnetic thresholds at particular times of the storm were studied. The theoretical geomagnetic thresholds were calculated with cosmic ray particle tracing in the magnetic field of the disturbed magnetosphere described by Ts01 model. The experimental geomagnetic thresholds were specified by spectrographic global survey according to the data of cosmic ray registration by the global station network. 相似文献
15.
Vlastislav Červený Antoine B. Frangié Reviewer J. Vaněk 《Studia Geophysica et Geodaetica》1982,26(3):238-253
Summary The effects of the causal absorption on seismic wave fields are discussed. Several computer programs for the evaluation of seismic wave fields have been modified to include the causal absorption. The programs are based mostly on matrix and ray methods. The main attention is paid to ray methods. A simple plane wave algorithm is described, which can be used along the ray to compute elementary seismograms in dissipative media. For a certain class of signals, commonly used in seismology, its application does not require any convolution or frequency-domain computations, and is only slightly slower than the evaluation of elementary seismograms in non-dissipative media. The derived algorithm even describes the velocity dispersion related to absorption and the decrease of the prevailing frequency of the signal well. Numerical examples are presented, both exact and approximate. 相似文献
16.
A note on dynamic ray tracing in ray-centered coordinates in anisotropic inhomogeneous media 总被引：1，自引：0，他引：1
V. Červený 《Studia Geophysica et Geodaetica》2007,51(3):411-422
Dynamic ray tracing plays an important role in paraxial ray methods. In this paper, dynamic ray tracing systems for inhomogeneous
anisotropic media, consisting of four linear ordinary differential equations of the first order along the reference ray, are
studied. The main attention is devoted to systems expressed in a particularly simple choice of ray-centered coordinates, here
referred to as the standard ray-centered coordinates, and in wavefront orthonormal coordinates. These two systems, known from
the literature, were derived independently and were given in different forms. In this paper it is proved that both systems
are fully equivalent. Consequently, the dynamic ray tracing system, consisting of four equations in wavefront orthonormal
coordinates, can also be used if we work in ray-centered coordinates, and vice versa.
vcerveny@seis.karlov.mff.cuni.cz 相似文献
17.
本文采用随机介质描述地球内部大尺度背景场上存在的小尺度不均匀性和自组织结构;文中没有采用传统的层状网格结构,提出了块状结构来描述复杂的二维自组织结构的方法,分别以高斯型、指数型和von Karman型等自相关函数描述各向同性和各向异性非均匀分布自组织特征;采用射线追踪分析了不同分布特征自组织结构对地震波运动学和动力学特征参数的影响;结果表明,由于地球内部介质的自组织性存在,地震波射线轨迹可能发生明显的畸变; 不同偏移距处,地震反射振幅减弱或增强; 自组织结构从高斯型到von Karman型,在小尺度上表现更大的非均匀性,因此走时和振幅表征依次更强的平均效应. 相似文献
18.
Seismic ray path variations in a 3D global velocity model 总被引：2，自引：0，他引：2
A three-dimensional (3D) ray tracing technique is used to investigate ray path variations of P, PcP, pP and PP phases in a global tomographic model with P wave velocity changing in three dimensions and with lateral depth variations of the Moho, 410 and 660 km discontinuities. The results show that ray paths in the 3D velocity model deviate considerably from those in the average 1D model. For a PcP wave in Western Pacific to East Asia where the high-velocity (1-2%) Pacific slab is subducting beneath the Eurasian continent, the ray path change amounts to 27 km. For a PcP ray in South Pacific where very slow (−2%) velocity anomalies (the Pacific superplume) exist in the whole mantle, the maximum ray path deviation amounts to 77 km. Ray paths of other phases (P, pP, PP) are also displaced by tens of kilometers. Changes in travel time are as large as 3.9 s. These results suggest that although the maximal velocity anomalies of the global tomographic model are only 1-2%, rays passing through regions with strong lateral heterogeneity (in velocity and/or discontinuity topography) can have significant deviations from those in a 1D model because rays have very long trajectories in the global case. If the blocks or grid nodes adopted for inversion are relatively large (3-5°) and only a low-resolution 3D model is estimated, 1D ray tracing may be feasible. But if fine blocks or grid nodes are used to determine a high-resolution model, 3D ray tracing becomes necessary and important for the global tomography. 相似文献
19.
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. 相似文献