复杂地质层面的网格化模型及两点射线追踪方法

本文针对三维复杂介质的两点射线追踪的适用性和有效性,提出了网格模型和层状模型相结合的建模方式,文中利用规则网格化的层状模型对复杂地质单元进行数值描述,一方面避免了用纯粹复杂地质网格体描述时需要的大量内存,另一方面结合运用规则化的层状模型,必有利于提高计算效率,在计算单元自动生成的射线追踪方法基础上,本文结合波前法提出一种快速、有效的两点射线追踪方法,利用反射波波前在地表形成的射线方位角和倾角坐标体系,能快速、准确地收敛到所有从源点到接收点的射线路径,亦可有效地解决射线路径的多值解难题,从而可适用于复杂地质模型的特点.     

三维非均匀地质模型中的逐段迭代射线追踪

地震射线追踪是地震定位、层析成像、偏移等领域的重要正演环节.随着这些领域研究的深入,针对传统的网格结构和层状结构在描述复杂地质模型遇到的很大困难,我们采用大小不等、形状各异的地质块组成的集合体来描述三维复杂地质模型,并用三角形面片来描述地质块之间的物性间断面,理论上可以描述任意复杂的地质模型.为适应任意非均匀速度分布的地质模型,基于费马原理,本文发展了与之相适应的逐段迭代射线追踪方法.该方法属于弯曲法范畴,对路径点采用一阶显式增量修正,相对于传统的迭代法,高效省时.数值试验表明,联合逐段迭代法和伪弯曲法的射线追踪扰动修正方案在三维复杂非均匀块状模型中有适用性和高效性.     

基于波前面三角形网格剖分的波前重建法三维射线追踪

在许多地震反演和偏移成像方法中,都要涉及到射线路径和旅行时的计算.本文将波前面三角形网格剖分和三维波前重建法射线追踪技术结合使用,实现了射线路径和旅行时的准确快速计算.三维波前重建法射线追踪过程中可以保证稳定合理的射线密度,克服了常规射线追踪方法存在阴影区的问题.波前面三角形网格剖分在描述和拆分波前面时更加准确有效,而且不需太多的网格数目,从而提高了射线追踪的精度和效率.该方法在三维复杂构造成像方面有独特的优势,目前在实际的Kirchhoff 偏移中的已经有相关应用.     

Indirect-approach method for specified-end points seismic ray tracing in three dimensional inhomogeneous media

ＩｎｄｉｒｅｃｔａｐｐｒｏａｃｈｍｅｔｈｏｄｆｏｒｓｐｅｃｉｆｉｅｄｅｎｄｐｏｉｎｔｓｓｅｉｓｍｉｃｒａｙｔｒａｃｉｎｇｉｎｔｈｒｅｄｉｍｅｎｓｉｏｎａｌｉｎｈｏｍｏｇｅｎｅｏｕｓｍｅｄｉａＣＨＡＯＦＡＮＸＵ（徐朝繁）ＸＩＡＮＫＡＮＧＺＨ...     

Common-ray tracing and dynamic ray tracing for S waves in a smooth elastic anisotropic medium

Anisotropic common S-wave rays are traced using the averaged Hamiltonian of both S-wave polarizations. They represent very practical reference rays for calculating S waves by means of the coupling ray theory. They eliminate problems with anisotropic-ray-theory ray tracing through some S-wave slowness-surface singularities and also considerably simplify the numerical algorithm of the coupling ray theory for S waves. The equations required for anisotropic-common-ray tracing for S waves in a smooth elastic anisotropic medium, and for corresponding dynamic ray tracing in Cartesian or ray-centred coordinates, are presented. The equations, for the most part generally known, are summarized in a form which represents a complete algorithm suitable for coding and numerical applications.     

三维场地波动传播的快速射线追踪法

本文基于Snell定律和Fermat原理对三维任意界面情况下的两点间射线追踪问题进行了研究，从同一条射线满足相同的射线参数出发，推得一个适用于任意界面情况下计算反（折）射点的一阶近似公式。结合迭代技术，给出了三维场地条件下射线追踪迭代算法的计算格式，并进行了三维场地射线追踪模拟计算。计算表明：计算速度相当快，且其计算精度可以根据需要满足要求。     

Quasi-shear wave ray tracing by wavefront construction in 3-D, anisotropic media

Wavefront construction (WFC) methods provide robust tools for computing ray theoretical traveltimes and amplitudes for multivalued wavefields. They simulate a wavefront propagating through a model using a mesh that is refined adaptively to ensure accuracy as rays diverge during propagation. However, an implementation for quasi-shear (qS) waves in anisotropic media can be very difficult, since the two qS slowness surfaces and wavefronts often intersect at shear-wave singularities. This complicates the task of creating the initial wavefront meshes, as a particular wavefront will be the faster qS-wave in some directions, but slower in others. Analogous problems arise during interpolation as the wavefront propagates, when an existing mesh cell that crosses a singularity on the wavefront is subdivided. Particle motion vectors provide the key information for correctly generating and interpolating wavefront meshes, as they will normally change slowly along a wavefront. Our implementation tests particle motion vectors to ensure correct initialization and propagation of the mesh for the chosen wave type and to confirm that the vectors change gradually along the wavefront. With this approach, the method provides a robust and efficient algorithm for modeling shear-wave propagation in a 3-D, anisotropic medium. We have successfully tested the qS-wave WFC in transversely isotropic models that include line singularities and kiss singularities. Results from a VTI model with a strong vertical gradient in velocity also show the accuracy of the implementation. In addition, we demonstrate that the WFC method can model a wavefront with a triplication caused by intrinsic anisotropy and that its multivalued traveltimes are mapped accurately. Finally, qS-wave synthetic seismograms are validated against an independent, full-waveform solution.     

弯曲射线追踪中Dijkstra算法的改进与实现

文章针对图论中寻找最短路径的Dijkstra算法内存占用量大,效率低的缺点,对该算法进行了改进,修改后的算法计算效率是原来的四倍,内存使用量和图中节点数呈线性关系.在此基础上,用新算法求出了激发点和接收点的最短走时路径,并由激发接收点的旅行时结合联合迭代法对理论模型和实际场地进行了反演.结果表明：和直射线追踪相比,弯曲射线路径能更好地反演出地质体内部的速度场分布.理论模型和实际探测结果证实改进后的算法是有效的.     

Theory of anisotropic dynamic ray tracing in ray-centred coordinates

A 4×4-propagator matrix formalism is presented for anisotropic dynamic ray tracing, including the propagation across curved interfaces. The computations are organised in the same way as in ervený's well-known isotropic propagator matrix formalism. Attention is paid to cases where double eigenvalues of the Christoffel matrix result in unstable expressions in the dynamic ray tracing system, but where geometrical spreading is well-behaved.     

椭球坐标系下多震相地震射线追踪

地震射线追踪方法技术在地震学领域有着较为广泛的应用,然而大多数算法建立在直角坐标系或球坐标系下,实际地球并非完美的球体,而是两极略扁的椭球体,因此,球坐标系下计算结果与真实情况存在一定误差.传统的做法一般是在球坐标系下进行计算,而后进行椭球校正.本文提出了一种直接在椭球体模型中采用分区多步最短路径算法进行多震相地震射线追踪的方法技术,实现了椭球坐标系下多震相地震波射线路径追踪和走时计算.与解析解的对比表明：该算法具有较高的计算精度,适用于任意形状的椭球体,且不需要进行额外的走时校正.数值模拟结果表明,计算所得P波和PcP反射波的走时与AK135走时表的误差小于0.1 s.当震中距较大时,使用球对称模型和椭球体模型计算所得的走时差异显著,说明采用椭球坐标系的必要性.

     

Modeling,ray tracing,and block nonlinear travel-time inversion in 3D

We describe an integrated forward and inverse three-dimensional modeling system that can deal with complex geological structures. The system has been designed to handle large-scale problems by using a distributed approach. It uses seismic ray tracing for forward simulation, time-to-depth mapping, and nonlinear travel-time inversion.A novel decomposition method is our tool to attack large-scale problems in a parallel approach. The system is fully implemented and we demonstrate its performance with synthetic examples.This research was partially supported by the National Science Foundation under SBIR Grant III-9300992.     

基于牛顿迭代法求解群速度的地震波射线追踪模拟

地下介质中普遍存在着各向异性,当前基于各向异性的地震波射线追踪多是在弱各向异性介质中进行且采用群速度近似表示方法,这些近似方法在强各项异性介质中会导致很大误差而无法真正模拟地震波的传播规律.根据地下普遍存在各向异性的事实和地震波基本传播规律,提出利用牛顿迭代法高效求解群速度,基于Paraview平台 自动化构建三维地质...     

Effects of bathymetry on tsunami propagation: Application of ray tracing to tsunamis

Ray tracing of seismic surface waves is applied to tsunami propagation to examine bathymetric effect along its propagation path. Computations are made for trans-Pacific tsunamis and for near-field tsunamis in the Japan Sea. For tsunamis across the Pacific Ocean, the comparison to a uniform ocean shows that focusing and defocusing, due to bathymetry, are significant for some combinations of source and receiver. For example, the refraction of rays is predominant at the East Pacific Rise for the tsunami from Chile. The tsunamis in the Japan Sea are strongly affected by the shallow Yamato Rise. The predicted arrival time and amplitude distribution generally agree with the observations from an actual tsunami. Since the computation can be made very quickly, the method is useful for preliminary analysis of tsunami propagation, such as in an operational warning system or in the determination of computational area for finite-difference computation.     

A rapid and accurate two-point ray tracing method in horizontally layered velocity model

Introduction The determination of the ray path between a source and a receiver is fundamental to many seismic problems such as earthquake location, travel time inversion, computation of synthetic seismograms, depth migration, and seismic tomography. The calculation of travel times and their derivatives depend on the determination of ray paths, which is the most time-consuming work in seismic tomography. Therefore, a fast and accurate ray tracing method is particularly important. Many ray-traci…     

最短路径算法下三维层状介质中多次波追踪

本文使用改进后的最短路径算法（MSPM）结合分区多步计算技术实现了三维复杂层状起伏介质中的多次透射、反射及转换波波前传播的数值模拟,以及相应走时和射线路径的跟踪计算.其原理是将三维复杂层状模型按速度界面分成若干个独立的计算区域,采用分步计算技术进行多次波的跟踪计算.基于多次波是通过速度界面简单的入射、透射、反射及转换波按一定规律及原理的不同组合,因此可实施分区多步计算技术.数值模拟实例及误差分析表明分区多步计算技术具有单步最短路径算法中的诸多优点:算法简单、数值计算稳健、计算精度高、速度快及全球解等,因此是解决多次波跟踪计算行之有效的方法.     

一种最短路径射线追踪的快速算法

为提高最短路径射线追踪的精度,需要增加模型的剖分网格和离散节点,并增加子波传播方向,或者采用其他方法改善计算结果,这些处理会带来大量的额外计算.本文的快速算法改进了波前点的管理和子波传播的计算这两项耗时的工作,较大幅度地提高了传统算法的效率.在波前点的管理上,采用按时间步划分区间的方法,实现了波前点的桶排序管理,其效率高于传统方法中常用的堆排序算法. 在子波传播的计算上,利用斯奈尔定律,同时参考来自邻近节点的波的走时,来限定当前子波传播的有效区域,排除大量不需要计算的子波传播方向. 模型实算表明,本文快速算法的计算速度是传统方法的几倍至十多倍.     

Point-to-curve ray tracing

Point-to-curve ray tracing is an attempt at dealing with multiplicity of solutions to a generic boundary-value problem of ray tracing. In a point-to-curve tracing (P2C) the input parameters of the boundary-value problem (BVP), such as the ends of the ray, are allowed to vary along a curve. The solutions of the BVP automatically wander from one solution branch to another generating a nearly complete multi-valued solution of the BVPs.A procedure for transforming an arbitrary iterative algorithm, solving a ray tracing BVP to a corresponding P2C algorithm, is presented. Bifurcations of the solution curve of the P2C problem at caustics are studied and an algorithm for obtaining the bifurcating branches is developed. In particular, transition from real rays to complex rays in a caustic shadow offers an additional link between otherwise disconnected solution curves of the P2C problem. The topological structure of a generic solution curve and its implications for the algorithm are studied.     

A note on dynamic ray tracing in ray-centered coordinates in anisotropic inhomogeneous media

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     

多步法快速射线追踪与圆台型高精度走时外插计算

为更好地适应复杂构造的地震偏移成像,本文提出了一套快速射线追踪算法和一种高精度的走时外插计算方法.采用线性多步法的预测-校正公式求解射线追踪方程组,与传统的四阶Runge-Kutta法相比,提高了计算效率.在网格节点上的走时计算中,应用一种基于圆台的外插方法,该方法以射线的方向为轴确定圆台,将轴上的走时外插到圆台内的网格节点上.与传统的矩形体外插方法相比,圆台走时外插方法提高了计算精度,且具有更好的稳定性.另外,该方法利用稀疏分布的射线即可获得高精度的走时表,节省计算量,对复杂构造的偏移成像非常有利,尤其是三维偏移.最后通过逆散射偏移成像算例,验证了算法的有效性和适用性.

     

长白山天池火山研究区地壳三维界面与速度分布联合成像

迅速发展并得到广泛应用的空间深地震测深技术通过采用三维数据采集的观测系统,利用三维层析技术,以求获得区域地壳三维分层结构和三维速度分布图像。它与传统的二维宽角反射－折射剖面技术相结合,可以有效地研究区域性地壳结构,特别是壳内深断裂和低速层的空间展布特征。Kanasewich等(1985)给出了利用空间深地震测深资料重建地壳三维界面的方法,我们在其方法的基础之上,进一步提出地壳三维界面和速度分布联合反演的方法。正问题的计算是在Chander(1977)关于三维平界面的快速两点追踪算法上的改进,在获得界面三维反演结果的基础之上,利用剩余走时残差,采用模型不分块反演技术(Tarantola,Nercession 1986)重建地壳三维速度图像。 1998年国家地震局地球物理勘探中心在长白山火山区实施了三维深地震测深观测,目的是研究天池火山的岩浆系统。利用本次实验所获得的780余个PmP波走时数据,采用上述的方法重建了研究区莫霍界面和地壳三维速度分布图像。研究结果表明,本区莫霍界面由北东方向向南逐渐加深,在天池火山口下达最深,并且被一些可能存在的地壳厚度陡变带(或深断裂)所切割。在东西方向莫霍面由西向东缓缓加深,其变化较南北方向缓和。特别值得注意的是,存在着一条近北东方向的莫霍面深度陡变带(或深断裂带)从天池火山口西部穿过,相应位置与马鞍山—三道白河地堑型断裂相一致,该断裂带可能对天山火山喷发时岩浆的运移起到重要作用。深度为15km和25km的P波速度图像表明,在天池火山口下分布着近南北走向明显的低P波速度分布,其南北方向延伸的范围约为80－90公里左右。比较这两个深度上的低P波速度体的分布特点,可以看出这个低P波速度体尺度随深度逐渐变小,但在25km深度处仍清晰可见,这表明该区岩浆自上地幔侵入地壳的“痕迹”,这也意味着,长白山天池火山的岩浆系统极有可能延伸到上地幔或更深一些。