共查询到20条相似文献，搜索用时 52 毫秒
1.
The performance of a 3D prestack migration of the Kirchhoff type can be significantly enhanced if the computation of the required stacking surface is replaced by an efficient and accurate method for the interpolation of diffraction traveltimes. Thus, input traveltimes need only be computed and stored on coarse grids, leading to considerable savings in CPU time and computer storage. However, interpolation methods based on a local approximation of the traveltime functions fail in the presence of triplications of the wavefront or later arrivals. This paper suggests a strategy to overcome this problem by employing the coefficients of a hyperbolic traveltime expansion to locate triplications and correct for the resulting errors in the interpolated traveltime tables of first and later arrivals. 相似文献
2.
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. 相似文献
3.
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. 相似文献
4.
振幅随偏移距变化是描述储层特征的重要方法之一,保幅偏移方法就是使偏移剖面能够反映出振幅随偏移距的变化.本论文中的保幅偏移是以走时为基础,主要的方法是采用走时的双曲线展开法,通过走时的二阶空间导数来确定波前曲率.该方法通过建立在大网格上的走时表来确定插值系数,将大网格插值成为较为精细的网格,这样就节省了数据的存储空间.对于相同的网格密度,通过插值来计算走时表比采用程函方程有限差分法直接计算走时要节省5至6倍的时间.走时的插值系数还可以用来计算几何扩散因子、权函数,不仅提高了成像质量,还大大节省了计算时间. 相似文献
5.
Geometrical spreading plays an important role for amplitude preserving migration, which is a very time-consuming process. In order to achieve efficiency in terms of computational time and, particularly, storage space, we propose a method to determine geometrical spreading from coarsely gridded traveltime tables. The method is based on a hyperbolic traveltime expansion and provides also a fast and accurate algorithm for the interpolation of traveltimes, including the interpolation of complete shots. Examples demonstrate the applicability of the method to isotropic and anisotropic media. 相似文献
6.
7.
8.
A new ray-tracing method called linear traveltime interpolation (LTI) is proposed. This method computes traveltimes and raypaths in a 2D velocity structure more rapidly and accurately than other conventional methods. The LTI method is formulated for a 2D cell model, and calculations of traveltimes and raypaths are carried out only on cell boundaries. Therefore a raypath is considered to be always straight in a cell with uniform velocity. This approach is suitable to tomography analysis. The algorithm of LTI consists of two separate steps: step 1 calculates traveltimes on all cell boundaries; step 2 traces raypaths for all pairs of receivers and the shot. A traveltime at an arbitrary point on a cell boundary is assumed to be linearly interpolated between traveltimes at the adjacent discrete points at which we calculate traveltimes. Fermat's principle is used as the criterion for choosing the correct traveltimes and raypaths from several candidates routinely. The LTI method has been compared numerically with the shooting method and the finite-difference method (FDM) of the eikonal equation. The results show that the LTI method has great advantages of high speed and high accuracy in the calculation of both traveltimes and raypaths. The LTI method can be regarded as an advanced version of the conventional FDM of the eikonal equation because the formulae of FDM are independently derived from LTI. In the process of derivation, it is shown theoretically that LTI is more accurate than FDM. Moreover in the LTI method, we can avoid the numerical instability that occurs in Vidale's method where the velocity changes abruptly. 相似文献
9.
10.
11.
最短路径射线追踪算法,用预先设置的网络节点的连线表示地震波传播路径,当网络节点稀疏时,获得的射线路径呈之字形,计算的走时比实际走时系统偏大. 本文在波前扩展和反向确定射线路径的过程中,在每个矩形单元内,通过对某边界上的已知走时节点的走时进行线性插值,并利用Fermat原理即时求出从该边界到达其他边界节点的最小走时及其子震源位置和射线路径,发展了相应的动态网络算法. 从而克服了最短路径射线追踪算法的缺陷,大大提高了最小走时和射线路径的计算精度. 相似文献
12.
三维地震波走时计算技术是三维地震反演、层析成像、偏移成像等诸多地震数据处理技术中非常重要的正演计算工具.为了获得精度高且兼顾效率的三维走时计算方法:首先,在常规双线性插值公式推导过程中,充分利用平面波双线性假设的结论,获得了二元极小值超越方程的解析解,进而推导出了准确的局部走时计算公式,同时构造性地证明了该计算公式满足地震波的传播规律和Eikonal方程;其次,引入迎风差分的基本思想,提出迎风双线性插值的局部走时计算策略,该计算策略能简化算法、提高效率且保证无条件稳定性;然后,将上述计算公式和迎风双线性插值策略与常规快速推进法中的窄带技术结合,获得了一种新的基于快速推进迎风双线性插值法的三维地震波走时计算方法;最后,通过精度和效率分析检验了新算法的精度、效率和正确性,并通过计算实例验证了算法在面对复杂介质时的稳定性和有效性. 相似文献
13.
为更好地适应复杂构造的地震偏移成像,本文提出了一套快速射线追踪算法和一种高精度的走时外插计算方法.采用线性多步法的预测-校正公式求解射线追踪方程组,与传统的四阶Runge-Kutta法相比,提高了计算效率.在网格节点上的走时计算中,应用一种基于圆台的外插方法,该方法以射线的方向为轴确定圆台,将轴上的走时外插到圆台内的网格节点上.与传统的矩形体外插方法相比,圆台走时外插方法提高了计算精度,且具有更好的稳定性.另外,该方法利用稀疏分布的射线即可获得高精度的走时表,节省计算量,对复杂构造的偏移成像非常有利,尤其是三维偏移.最后通过逆散射偏移成像算例,验证了算法的有效性和适用性. 相似文献
14.
E. BRÜCKL 《Geophysical Prospecting》1987,35(9):973-992
We consider multiply covered traveltimes of first or later arrivals which are gathered along a refraction seismic profile. The two-dimensional distribution of these traveltimes above a coordinate frame generated by the shotpoint axis and the geophone axis or by the common midpoint axis and the offset axis is named a traveltime field. The application of the principle of reciprocity to the traveltime field implies that for each traveltime value with a negative offset there is a corresponding equal value with positive offset. In appendix A procedures are demonstrated which minimize the observational errors of traveltimes inherent in particular traveltime branches or complete common shotpoint sections. The application of the principle of parallelism to an area of the traveltime field associated with a particular refractor can be formulated as a partial differential equation corresponding to the type of the vibrating string. The solution of this equation signifies that the two-dimensional distribution of these traveltimes may be generated by the sum of two one-dimensional functions which depend on the shotpoint coordinate and the geophone coordinate. Physically, these two functions may be interpreted as the mean traveltime branches of the reverse and the normal shot. In appendix B procedures are described which compute these two functions from real traveltime observations by a least-squares fit. The application of these regressed traveltime field data to known time-to-depth conversion methods is straightforward and more accurate and flexible than the use of individual traveltime branches. The wavefront method, the plus-minus method, the generalized reciprocal method and a ray tracing method are considered in detail. A field example demonstrates the adjustment of regressed traveltime fields to observed traveltime data. A time-to-depth conversion is also demonstrated applying a ray tracing method. 相似文献
15.
16.
Tariq Alkhalifah 《Geophysical Prospecting》2002,50(4):373-382
A linearized eikonal equation is developed for transversely isotropic (TI) media with a vertical symmetry axis (VTI). It is linear with respect to perturbations in the horizontal velocity or the anisotropy parameter η. An iterative linearization of the eikonal equation is used as the basis for an algorithm of finite-difference traveltime computations. A practical implementation of this iterative technique is to start with a background model that consists of an elliptically anisotropic, inhomogeneous medium, since traveltimes for this type of medium can be calculated efficiently using eikonal solvers, such as the fast marching method. This constrains the perturbation to changes in the anisotropy parameter η (the parameter most responsible for imaging improvements in anisotropic media). The iterative implementation includes repetitive calculation of η from traveltimes, which is then used to evaluate the perturbation needed for the next round of traveltime calculations using the linearized eikonal equation. Unlike isotropic media, interpolation is needed to estimate η in areas where the traveltime field is independent of η, such as areas where the wave propagates vertically.
Typically, two to three iterations can give sufficient accuracy in traveltimes for imaging applications. The cost of each iteration is slightly less than the cost of a typical eikonal solver. However, this method will ultimately provide traveltime solutions for VTI media. The main limitation of the method is that some smoothness of the medium is required for the iterative implementation to work, especially since we evaluate derivatives of the traveltime field as part of the iterative approach. If a single perturbation is sufficient for the traveltime calculation, which may be the case for weak anisotropy, no smoothness of the medium is necessary. Numerical tests demonstrate the robustness and efficiency of this approach. 相似文献
Typically, two to three iterations can give sufficient accuracy in traveltimes for imaging applications. The cost of each iteration is slightly less than the cost of a typical eikonal solver. However, this method will ultimately provide traveltime solutions for VTI media. The main limitation of the method is that some smoothness of the medium is required for the iterative implementation to work, especially since we evaluate derivatives of the traveltime field as part of the iterative approach. If a single perturbation is sufficient for the traveltime calculation, which may be the case for weak anisotropy, no smoothness of the medium is necessary. Numerical tests demonstrate the robustness and efficiency of this approach. 相似文献
17.
Hui Sun Jian-Guo Sun Zhang-Qing Sun Fu-Xing Han Zhi-Qiang Liu Ming-Chen Liu Zheng-Hui Gao Xiu-Lin Shi 《应用地球物理》2017,14(1):56-63
3D traveltime calculation is widely used in seismic exploration technologies such as seismic migration and tomography. The fast marching method (FMM) is useful for calculating 3D traveltime and has proven to be efficient and stable. However, it has low calculation accuracy near the source, which thus gives it low overall accuracy. This paper proposes a joint traveltime calculation method to solve this problem. The method firstly employs the wavefront construction method (WFC), which has a higher calculation accuracy than FMM in calculating traveltime in the small area near the source, and secondly adopts FMM to calculate traveltime for the remaining grid nodes. Due to the increase in calculation precision of grid nodes near the source, this new algorithm is shown to have good calculation precision while maintaining the high calculation efficiency of FMM, which is employed in most of the computational area. Results are verified using various numerical models. 相似文献
18.
The first-order perturbation theory is used for fast 3D computation of quasi-compressional (qP)-wave traveltimes in arbitrarily anisotropic media. For efficiency we implement the perturbation approach using a finite-difference (FD) eikonal solver. Traveltimes in the unperturbed reference medium are computed with an FD eikonal solver, while perturbed traveltimes are obtained by adding a traveltime correction to the traveltimes of the reference medium. The traveltime correction must be computed along the raypath in the reference medium. Since the raypath is not determined in FD eikonal solvers, we approximate rays by linear segments corresponding to the direction of the phase normal of plane wavefronts in each cell. An isotropic medium as a reference medium works well for weak anisotropy. Using a medium with ellipsoidal anisotropy as a background medium in the perturbation approach allows us to consider stronger anisotropy without losing computational speed. The traveltime computation in media with ellipsoidal anisotropy using an FD eikonal solver is fast and accurate. The relative error is below 0.5% for the models investigated in this study. Numerical examples show that the reference model with ellipsoidal anisotropy allows us to compute the traveltime for models with strong anisotropy with an improved accuracy compared with the isotropic reference medium. 相似文献
19.
3D P‐wave traveltime computation in transversely isotropic media using layer‐by‐layer wavefront marching 
下载免费PDF全文
下载免费PDF全文 Jiangtao Hu Junxing Cao Huazhong Wang Xingjian Wang Renfei Tian 《Geophysical Prospecting》2018,66(7):1303-1314
Subsurface rocks (e.g. shale) may induce seismic anisotropy, such as transverse isotropy. Traveltime computation is an essential component of depth imaging and tomography in transversely isotropic media. It is natural to compute the traveltime using the wavefront marching method. However, tracking the 3D wavefront is expensive, especially in anisotropic media. Besides, the wavefront marching method usually computes the traveltime using the eikonal equation. However, the anisotropic eikonal equation is highly non‐linear and it is challenging to solve. To address these issues, we present a layer‐by‐layer wavefront marching method to compute the P‐wave traveltime in 3D transversely isotropic media. To simplify the wavefront tracking, it uses the traveltime of the previous depth as the boundary condition to compute that of the next depth based on the wavefront marching. A strategy of traveltime computation is designed to guarantee the causality of wave propagation. To avoid solving the non‐linear eikonal equation, it updates traveltime along the expanding wavefront by Fermat's principle. To compute the traveltime using Fermat's principle, an approximate group velocity with high accuracy in transversely isotropic media is adopted to describe the ray propagation. Numerical examples on 3D vertical transverse isotropy and tilted transverse isotropy models show that the proposed method computes the traveltime with high accuracy. It can find applications in modelling and depth migration. 相似文献
20.
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. 相似文献