共查询到20条相似文献,搜索用时 250 毫秒
1.
2.
本文从能量泛函的角度着手,基于速度梯度,采用最速下降法推导了基于偏微分方程的速度模型平滑公式,用于解决射线类偏移成像过程当中速度模型的平滑处理问题.同时针对偏微分方程速度模型平滑系数中阈值k的选取对速度模型空间结构的影响,在简单的高速体速度模型上分析了不同k值的选取对原始速度模型空间结构的改变,并通过射线路径和时间场的对比分析证明偏微分方程速度模型平滑处理相对于卷积算子平滑的优越性.最后通过在Marmousi、Sigsbee2A原始速度模型以及平滑处理后的速度模型上的成像结果分析验证偏微分方程模型平滑的有效性. 相似文献
3.
4.
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. 相似文献
5.
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. 相似文献
6.
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 相似文献
7.
本文以行星际大尺度扇形磁场为基础,建立了研究太阳宇宙线在行星际传播较合理的扇形坐标系,并导出了太阳宇宙线在这种坐标系中的传播方程。同时研究了太阳宇宙线的一种扩散模型,扩散方程的解依赖于耀斑源在太阳上的磁径度,并得到了粒子强度随时间变化的特性,与观测结果比较,理论与观测资料符合较好。 相似文献
8.
P. HUBRAL 《Geophysical Prospecting》1980,28(2):221-239
For a 3D velocity model of curved first order interfaces and layer velocities which are arbitrary smooth functions of the space coordinates, the normal moveout (NMO)-velocity can be computed by numerically integrating a system of first order ordinary differential equations for a hypothetical wavefront that originates at the normal incidence point of the normal ray and moves up along the ray to the common mid-point of the common datum point (CDP) profile. 相似文献
9.
10.
针对线性走时插值算法(LTI)不能正确追踪逆向传播射线的问题, 目前已提出多种改进算法, 如扩张收缩LTI算法、 循环计算LTI算法、 动态网络最短路径射线追踪算法等, 但这些算法的计算效率普遍偏低. 在分析各种改进LTI算法的优劣后, 本文提出了改进动态网络最短路径射线追踪算法. 该改进算法依据波的传播规律以及LTI算法的基本方程, 排除动态网络最短路径射线追踪算法中大量冗余节点计算, 并采用传统的二叉树堆排序算法对波前阵列节点进行管理. 数值算例表明, 本文提出的改进算法具有较高的计算效率, 其计算效率是动态网络最短路径射线追踪算法的4.5—30倍, 是原始LTI算法的2—6.5倍; 当动态网络最短路径射线追踪算法采用堆排序算法时, 改进算法的计算效率是其3.5—15倍. 相似文献
11.
最短路径射线追踪算法,用预先设置的网络节点的连线表示地震波传播路径,当网络节点稀疏时,获得的射线路径呈之字形,计算的走时比实际走时系统偏大. 本文在波前扩展和反向确定射线路径的过程中,在每个矩形单元内,通过对某边界上的已知走时节点的走时进行线性插值,并利用Fermat原理即时求出从该边界到达其他边界节点的最小走时及其子震源位置和射线路径,发展了相应的动态网络算法. 从而克服了最短路径射线追踪算法的缺陷,大大提高了最小走时和射线路径的计算精度. 相似文献
12.
走时计算已广泛用于地震建模、成像及速度分析等诸方面.基于地震波场的Hamilton力学性质,本文探索应用适于Hamilton力学的计算方法──辛几何算法对地震射线 进行走时及路径的计算.利用辛几何算法和属于耗散算法的四阶 Runge-Kutta算法进行射线 路径和走时的计算对比,结果表明这两类算法数值精度相当,但辛几何算法的速度却快了3 倍.本文还利用二阶Euler型辛差分格式对Marmousi模型进行了射线追踪计算,结果显示所 得射线具有良好的光滑性和较强的阴影区穿透能力. 相似文献
13.
A fast and robust method for two-point ray tracing in one-dimensional layered media is presented. This method is applicable to layered models with constant or linearly varying isotropic layer velocity. For given model properties and source and receiver positions, a ray path can be uniquely determined once its ray parameter (i.e. horizontal slowness) is known. The ray parameter can be obtained by numerically solving the nonlinear offset (i.e. source–receiver horizontal distance) equation using Newton's method, which generally works well at near and mid offsets. However, Newton's method becomes hard to converge at large offsets due to the oversensitivity of offset to ray parameter. Based on the analysis of the characteristic of the offset equation, a modified ray parameter is proposed and used to replace the generic ray parameter in numerical calculation. Numerical experiments show that the iteration process becomes stable and converges rapidly with the modified ray parameter. Moreover, a rational function that asymptotically approximates the shape of the offset equation is introduced for obtaining good initial estimates of the modified ray parameter. Numerical tests show that this method is robust in any situation, and an accurate ray parameter can be obtained within two or three iterations for a wide range of model velocity structure and source–receiver distance. Furthermore, the proposed two-point ray tracing method is easy to implement. 相似文献
14.
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. 相似文献
15.
本文给出了电离层中存在电离运动时的高频无线电波的多普勒频移的表达式,讨论了用高频多普勒阵研究电离层声重波扰动有关的若干电波传播问题。 相似文献
16.
TTI介质是石油地震勘探领域最常用的各向异性介质,快速计算TTI介质射线路径和走时信息有重要的研究意义.TTI介质传统运动学射线追踪方法一般基于任意弹性介质射线方程,利用Bond变换或者四阶张量变换来处理复杂的21个弹性参数,因而非常耗时.实际野外对称轴统一的TTI介质模型,一般可以看成VTI介质模型旋转一定角度获得.为此,本文推导了三维VTI介质射线追踪方程,提出先在本构坐标系中进行VTI介质射线追踪,再通过坐标旋转将射线路径旋转至观测坐标系中,获得TTI介质射线路径.数值模型计算表明该方法高效和精确,较传统方法效率提高了近4倍.在强各向异性等特殊情况下,体波波前面都与理论群速度面一致. 相似文献
17.
18.
19.
走时层析成像是根据首波走时重建声波、地震波或电磁波波速的层析技术,其所用到的射线路径依赖于波速分布。本文提出了一个依据费马原理寻找射线路径的新方法。费马原理指出,从发射点到接点的射线首波的走时极小,因此寻找射线路径的问题是一个最优化问题。 相似文献
20.
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. 相似文献