基于走时的保幅偏移方法

振幅随偏移距变化是描述储层特征的重要方法之一,保幅偏移方法就是使偏移剖面能够反映出振幅随偏移距的变化.本论文中的保幅偏移是以走时为基础,主要的方法是采用走时的双曲线展开法,通过走时的二阶空间导数来确定波前曲率.该方法通过建立在大网格上的走时表来确定插值系数,将大网格插值成为较为精细的网格,这样就节省了数据的存储空间.对于相同的网格密度,通过插值来计算走时表比采用程函方程有限差分法直接计算走时要节省5至6倍的时间.走时的插值系数还可以用来计算几何扩散因子、权函数,不仅提高了成像质量,还大大节省了计算时间.     

Determination of geometrical spreading from traveltimes

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.     

表驱三维角度域Kirchhoff叠前时间偏移成像方法

针对传统Kirchhoff叠前时间偏移方法的一些不足,以及振幅随入射角、方位角变化(AVA/AVAZ)分析的需要,本文提出一种基于射线理论的三维叠前时间偏移角度域成像方法.它通过横向均匀介质中稳健的射线追踪建立单程波走时和传播角度的数值表,然后在此基础上估算反射波双程走时以及在界面处传播的方位角和入射角,最后基于脉冲响应叠加原理获得三维构造图像和方位\|角度域共成像点道集.与传统方法不同之处在于,上述过程均考虑了地震波在垂向变速介质中的射线弯曲效应和三维传播特征,有利于准确提取随入射角和方位角变化的振幅和时差信息.理论模型合成数据和实际地震资料测试结果展示了方法的优越性与实用性.     

Secondorder interpolation of traveltimes

To carry out a 3D prestack migration of the Kirchhoff type is still a task of enormous computational effort. Its efficiency can be significantly enhanced by employing a fast traveltime interpolation algorithm. High accuracy can be achieved if secondorder spatial derivatives of traveltimes are included in order to account for the curvature of the wavefront. We suggest a hyperbolic traveltime interpolation scheme that permits the determination of the hyperbolic coefficients directly from traveltimes sampled on a coarse grid, thus reducing the requirements in data storage. This approach is closely related to the paraxial ray approximation and corresponds to an extension of the wellknown     method to arbitrary heterogeneous and complex media in 3D. Application to various velocity models, including a 3D version of the Marmousi model, confirms the superiority of our method over the popular trilinear interpolation. This is especially true for regions with strong curvature of the local wavefront. In contrast to trilinear interpolation, our method also provides the possibility of interpolating source positions, and it is 56 times faster than the calculation of traveltime tables using a fast finitedifference eikonal solver.     

基于快速推进迎风双线性插值法的三维地震波走时计算

三维地震波走时计算技术是三维地震反演、层析成像、偏移成像等诸多地震数据处理技术中非常重要的正演计算工具.为了获得精度高且兼顾效率的三维走时计算方法:首先,在常规双线性插值公式推导过程中,充分利用平面波双线性假设的结论,获得了二元极小值超越方程的解析解,进而推导出了准确的局部走时计算公式,同时构造性地证明了该计算公式满足地震波的传播规律和Eikonal方程;其次,引入迎风差分的基本思想,提出迎风双线性插值的局部走时计算策略,该计算策略能简化算法、提高效率且保证无条件稳定性;然后,将上述计算公式和迎风双线性插值策略与常规快速推进法中的窄带技术结合,获得了一种新的基于快速推进迎风双线性插值法的三维地震波走时计算方法;最后,通过精度和效率分析检验了新算法的精度、效率和正确性,并通过计算实例验证了算法在面对复杂介质时的稳定性和有效性.     

Second-order interpolation of later-arrival traveltimes

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.     

井间地震射线层析成像

本文采用射线走时对井间介质的地震波速进行层析成像，反演出井间介质的异常分布和构造．作者用快速弯曲射线迫踪法求解初值问题，用两步法和插值逼近法结合求解两点问题，避免了复杂运算，提高了运算效率和精度．反问题的数值计算用ＬＳＱＲ法．为了求取速度任意分布、对比度大的较为复杂介质模型的槽确初至时，本文采用动态程序法。该法灵活不受地质界面条件约束，且对介质速度分布未加任何限制。本文给出了几个算例．文中层析成像试验结果表明，此种层析成像算法是成功和有效的。     

三维复杂山地多级次群推进迎风混合法多波型走时计算

三维复杂山地条件下的各种地震波型的走时计算技术,可以直接用于复杂山地区域地震波运动学特性的分析、地震数据采集观测系统的设计以及直接基于三维复杂地表的地震数据处理技术的研发.为了在三维复杂地表条件下准确、灵活且稳定地计算各种地震波型的走时,提出一种多级次群推进迎风混合法.该算法利用不等距迎风差分法简洁稳定地处理三维复杂地表及附近的局部走时计算问题,利用计算精度不错的迎风双线性插值法处理绝大部分均匀正方体网格中的局部走时计算问题,利用群推进法模拟三维复杂地表条件下地震波前的扩展问题,利用多级次算法处理各种类型的地震波的走时计算问题.算法分析和计算实例表明:新方法具有很好的计算精度与效率,且能灵活稳定地处理三维复杂地表复杂介质条件下的多波型走时计算问题.     

EFFICIENT 3-D MIGRATION IN TWO STEPS*

The addition of a third dimension to seismic data acquisition offers significant improvement of our reflection images of geologic structures, but potentially at a substantial increase in data processing cost. Here we study an approach to migration of three-dimensional (3-D) stacked data that can considerably reduce the burden of computation and data manipulation. Cascading two-dimensional (2-D) migrations in orthogonal lateral directions, we obtain an efficient 3-D scheme that is exact for homogeneous media. As applied to examples of synthetic data from inhomogeneous media, the scheme introduces errors well below those attributable, in practice, to uncertainties in migration velocity. Application of this two-step approach, using both Kirchhoff-summation and finite-difference algorithms, to a 3-D seismic land survey over an area of known hydrocarbon traps additionally proves the importance of 3-D migration in stratigraphic imaging.     

基于Chebyshev多项式的弯曲射线Kirchhoff叠前时间偏移

弯曲射线Kirchhoff叠前时间偏移被认为是一种综合了叠前时间偏移效率和叠前深度偏移精度的方法.本文以保精度地减少高阶Kirchhoff叠前时间偏移走时计算量为目标,在分析了Chebyshev正交多项式性质的基础上,建立了Chebyshev多项式约简系数表,进而用模拟退火法对转换系数进行分段优化,从而实现了在大炮检距...     

密点速度分析技术在塔河油田西南部低幅构造中的应用

随着地震勘探精度的不断提高,常规速度分析方法已不能满足精细构造分析、储层描述及各向异性研究等对速度精度的要求.研究应用高密度叠前速度求取技术,自动拾取地震记录规则波的走时、振幅、波场相干性及炮点射线参数、检波点射线参数等信息,求取不同偏移距道上波峰时间处的反射波速度,并针对地质条件,采用空变时变速度滤波、频率滤波、倾角参数滤波等手段提高速度精度.将该技术应用于塔河油田西南部TFT地区1200km2叠前三维地震资料的精细速度分析,落实了主要目的层低幅度构造,确定了构造形态和要素;通过差异速度分析、AVO反射能量包络差及速度变化规律研究,对泥盆系东河砂岩进行了储层异常区带划分.     

模拟退火方法在三维速度模型地震波走时反演中的应用

采用块状建模以及三角形拼接的界面描述方式,并通过立方体速度网格线性插值获得块体内部的速度分布。正演过程中采用逐段迭代射线追踪方法计算三维复杂地质模型中的射线走时,并采用模拟退火方法进行了三维模型中的地震波走时反演研究。模型测试结果表明,使用的射线追踪和走时反演算法有效。     

数据域特征波反射走时反演

地震波的运动学信息(走时、斜率等)通常用于宏观速度建模.针对走时反演方法,一个基本问题是走时拾取或反射时差的估计.对于成像域反演方法,可以通过成像道集的剩余深度差近似计算反射波时差.在数据域中,反射地震观测数据是有限频带信号,如果不能准确地确定子波的起跳时间,难以精确地确定反射波的到达时间.另一方面,如果缺乏关于模型的先验信息,则很难精确测量自地下同一个反射界面的观测数据同相轴和模拟数据同相轴之间的时差.针对走时定义及时差测量问题,首先从叠前地震数据的稀疏表达出发,利用特征波场分解方法,提取反射子波并估计局部平面波的入射和出射射线参数.进一步,为了实现自动和稳定的走时拾取,用震相的包络极值对应的时间定义反射波的到达时,实现了立体数据中间的自动生成.理论上讲,利用包络极值定义的走时大于真实的反射波走时,除非观测信号具有无限带宽(即delta脉冲).然而,走时反演的目的是估计中-大尺度的背景速度结构,因此走时误差导致的速度误差仍然在可以接受的误差范围内.利用局部化传播算子及特征波聚焦成像条件将特征波数据直接投影到地下虚拟反射点,提出了一种新的反射时差估计方法.既避免了周期跳跃现象以及串层等可能性,又消除了振幅因素对时差测量的影响.最后,在上述工作基础之上,提出了一种基于特征波场分解的新型全自动反射走时反演方法(CWRTI).通过对泛函梯度的线性化近似,并用全变差正则化方法提取梯度的低波数部分,实现了背景速度迭代反演.在理论上,无需长偏移距观测数据或低频信息、对初始模型依赖性低且计算效率高,可以为后续的全波形反演提供可靠的初始速度模型.理论和实际资料的测试结果证明了本文方法的有效性.     

PS转换波界面二次源法射线追踪

界面二次源法是最近提出的一种最小走时射线追踪方法,尤其适合层状介质中走时和射线路径的计算.该方法相对于传统的最小走时树方法（如Moser法）,仅在物性界面上设置二次源,射线路径的方向只在层界面处发生改变,该方法最大程度地消除了射线路径的锯齿状现象,同时也避免了低变速区的射线路径多值现象,因此,它具有更高的追踪精度和效率.本文采用界面二次源法在各向同性介质中实现了PS转换波射线追踪,理论模型的计算证实了界面二次源法追踪PS转换波的准确性和高效性,同时该方法在各向异性介质中也很好地追踪出分离的PSV波和PSH波,因此该方法有利于横波分裂在地震勘探中的研究和应用     

TI介质局部角度域高斯束叠前深度偏移成像

各向异性射线理论基础上的局部角度域叠前深度偏移方法能够为深度域构造成像与基于角道集的层析反演提供有力支撑,但是对于复杂地质构造而言,高斯度叠前深度偏移在不失高效、灵活等特点的情况下,具有明显的精度优势.为此,本文研究局部角度域理论框架下的高斯束叠前深度偏移方法.为提高算法效率与实用性,文中讨论了一种从经典弹性参数表征的各向异性介质运动学和动力学射线方程演变而来的由相速度表征的简便形式,并提出了一种比较经济的各向异性高斯束近似合成方案.结合地震波局部角度域成像原理,讨论一种适合高斯束偏移的角度参数计算方法.国际上通用的理论模型合成数据试验表明：相比局部角度域Kirchhoff叠前深度偏移成像方法,本文方法具有更高的成像精度与抗噪能力,既适用于复杂构造成像,也可为TI介质深度域偏移速度分析与模型建立提供高效的偏移引擎.     

复杂山地随机介质GMM-ULTI法射线追踪

对复杂山地介质的非均质性以及介质中地震波运动学特征进行深入研究,对于提高复杂山地区域地震勘探的效果有着重要的理论意义和实际价值.为了研究复杂山地非均质性和该介质中地震波的一些运动特性,提出了一种复杂山地随机介质的建模方法和一种新的射线追踪算法.与常规算法相比,复杂山地随机介质的生成方法采用更贴近实际介质特点的梯度介质作为背景介质,并在模型生成过程中加入地形修正步骤;新提出的GMM-ULTI射线追踪算法,充分融合群推进法、迎风思想、走时插值法的优势,采用先计算走时后追踪射线路径的两步策略完成射线追踪.算法分析与计算实例表明:复杂山地随机介质的生成方法能灵活、精细且更贴近实际地刻画复杂山地介质的非均质特点;新射线追踪算法兼顾精度和效率、能无条件稳定且灵活地适应复杂山地随机介质的特点;同时基于对几个模型试算结果的分析也得出了复杂山地随机介质中的地震波的一些传播规律.     

2D inversion of refraction traveltime curves using homogeneous functions

A method using simple inversion of refraction traveltimes for the determination of 2D velocity and interface structure is presented. The method is applicable to data obtained from engineering seismics and from deep seismic investigations. The advantage of simple inversion, as opposed to ray‐tracing methods, is that it enables direct calculation of a 2D velocity distribution, including information about interfaces, thus eliminating the calculation of seismic rays at every step of the iteration process. The inversion method is based on a local approximation of the real velocity cross‐section by homogeneous functions of two coordinates. Homogeneous functions are very useful for the approximation of real geological media. Homogeneous velocity functions can include straight‐line seismic boundaries. The contour lines of homogeneous functions are arbitrary curves that are similar to one another. The traveltime curves recorded at the surface of media with homogeneous velocity functions are also similar to one another. This is true for both refraction and reflection traveltime curves. For two reverse traveltime curves, non‐linear transformations exist which continuously convert the direct traveltime curve to the reverse one and vice versa. This fact has enabled us to develop an automatic procedure for the identification of waves refracted at different seismic boundaries using reverse traveltime curves. Homogeneous functions of two coordinates can describe media where the velocity depends significantly on two coordinates. However, the rays and the traveltime fields corresponding to these velocity functions can be transformed to those for media where the velocity depends on one coordinate. The 2D inverse kinematic problem, i.e. the computation of an approximate homogeneous velocity function using the data from two reverse traveltime curves of the refracted first arrival, is thus resolved. Since the solution algorithm is stable, in the case of complex shooting geometry, the common‐velocity cross‐section can be constructed by applying a local approximation. This method enables the reconstruction of practically any arbitrary velocity function of two coordinates. The computer program, known as godograf , which is based on this theory, is a universal program for the interpretation of any system of refraction traveltime curves for any refraction method for both shallow and deep seismic studies of crust and mantle. Examples using synthetic data demonstrate the accuracy of the algorithm and its sensitivity to realistic noise levels. Inversions of the refraction traveltimes from the Salair ore deposit, the Moscow region and the Kamchatka volcano seismic profiles illustrate the methodology, practical considerations and capability of seismic imaging with the inversion method.     

3D P‐wave traveltime computation in transversely isotropic media using layer‐by‐layer wavefront marching

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.     

起伏地形下的高精度反射波走时层析成像方法

全球造山带及中国大陆中西部普遍具有强烈起伏的地形条件.复杂地形条件下的地壳结构成像问题像一面旗帜引领了当前矿产资源勘探和地球动力学研究的一个重要方向.深地震测深记录中反射波的有效探测深度可达全地壳乃至上地幔顶部,而初至波通常仅能探测上地壳浅部.为克服和弥补初至波探测深度的不足,本文基于前人对复杂地形条件下初至波成像的已有研究成果,采用数学变换手段将笛卡尔坐标系的不规则模型映射到曲线坐标系的规则模型,并将快速扫描方法与分区多步技术相结合,发展了反射波走时计算和射线追踪的方法.进而利用反射波走时反演,实现起伏地形下高精度的速度结构成像,从而为起伏地形下利用反射波数据高精度重建全地壳速度结构提供了一种全新方案.数值算例从正演计算精度、反演中初始模型依赖性、反演精度、纵横向分辨率以及抗噪性等方面验证了算法的正确性和可靠性.     

复杂地表条件下基于线性插值和窄带技术的地震波走时计算

为了在复杂地表条件下实现地震波走时计算,提出了一种基于线性插值和窄带技术的走时计算新方法.其中,线性插值用于局部走时计算,窄带技术用于局部波前捕获和追踪.为了逼近起伏地表,采用三角网和矩形网相结合的方法对速度模型进行剖分.为了得到局部走时计算公式,利用费马（Fermat）原理和关于入射点位置的限定条件.有关编程实践和数值试验表明：新方法不仅可以有效、灵活地处理地表高程的剧烈变化,而且还具有很好的适应性和稳定性,得到的计算结果满足波前传播规律.