基于Born波路径的高斯束初至波波形反演

为了提高表层速度反演精度,本文提出了一种新的波形反演方法.该方法只利用初至波波形信息以减少波形反演对初始模型的依赖性,降低反演多解性与稳定性.由于只利用初至波波形信息,所以该方法利用高斯束计算格林函数和正演波场,以减少正演计算量.为了避免庞大核函数的存储,该方法基于Born波路径,利用矩阵分解算法实现方向与步长的累加计算.将此基于Born波路径的初至波波形反演方法应用于理论模型实验,并与声波方程全波形反演和初至波射线走时层析方法相对比,发现该方法的反演效果略低于全波形反演方法,但明显优于传统初至波射线走时层析方法,而计算效率却与射线走时层析相当.同时,相对于全波形反演,本文方法对初始模型的依赖性也有所降低.     

基于声波测井资料的波动方程走时层析反演速度建模

地震勘探目标逐渐由浅层(<2000 m)转向中深层(2000～3500 m),高精度速度建模是实现精确地震成像的关键手段之一.目前,初至波或早至波速度层析建模方法已较为成熟,然而对深层速度建模更为重要的反射波层析建模方法依然面临走时拾取困难及计算量大等难题.针对当前中深层走时层析速度建模方法面临的难题,本文提出了一种基于声波测井资料的波动方程全波走时层析速度建模方法及策略.该方法首先通过早至波层析反演获得浅层及大尺度速度信息;其次,通过提取初次成像剖面的构造倾角信息,利用声波测井资料进行稀疏反演,进一步获得更高精度的速度模型;最后,以此为初始模型,进行全波层析反演获得最终偏移速度模型.数值结果表明,本文所提出的速度建模方法及策略不仅可以避免传统反射波层析建模方法因采用Born近似反偏移所带来的大计算量问题,且大幅提高了中深层速度建模的精度,进而能有效地提高中深层复杂构造区域的地震成像质量.     

高斯波包反射走时速度反演方法

扰动高斯波包理论指出,在Gabor域描述模型的扰动成分,且入射波场为短时宽带信号时,扰动波场可在时间域通过高斯波包算子描述.在此基础上通过拟合反射波的走时,提出一种速度反演方法.反射波走时残差利用地震道局部波形的互相关函数表示,以走时残差的二范数作为目标函数,优化目标函数实现对速度场的反演.基于一阶Born近似,利用扰动高斯波包理论推导出目标函数对速度场的梯度是本文理论部分的核心内容.梯度包括两部分:正传的背景波场与反传的扰动高斯波包之间的互相关,反传的背景波场和正传的扰动高斯波包之间的互相关.梯度表达式中背景波场和扰动波场均利用高斯波包算子模拟.计算梯度的具体算法中,如何模拟扰动波场,以及如何计算反射波的走时残差是两个要点,文中对此做了详细的讨论.数值实验进一步阐述了反演的实现策略,实验结果表明高斯波包反射走时速度反演方法和实现策略有效可行,并得到了理想的反演结果.     

弹性矢量波成像域联合层析反演

在地震弹性矢量波场框架下,推导了多波联合层析速度反演方程以及走时残差与角道集剩余曲率的转换关系式,提出了一种利用成像域角道集更新P波、S波速度的走时层析反演方法.其实现过程可以概括为:将弹性波多分量数据作为输入,基于高斯束实现矢量波场成像并提取角道集,利用层析反演方程求解慢度更新量,最终获得多波联合反演结果.模型试算和实际资料处理验证了该方法的反演效果,能够为弹性矢量波联合深度偏移提供高质量的叠前速度场.     

2D共炮时间域高斯波束偏移

针对传统射线方法在奇异区成像精度不高,而2D频率域高斯波束叠前深度偏移需要计算成像点处每个频率的格林函数,影响计算效率的问题,本文通过使用复走时代替实走时,改变频率域下成像公式的积分顺序,给出了在时间域下进行高斯波束偏移的方法和计算公式.本文使用复杂数值模型验证了2D时间域高斯波束叠前偏移方法的正确性,并同传统射线偏移成像结果做了对比.对比结果表明时间域高斯波束偏移在成像精度上优于传统射线偏移.     

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

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

双平方根波动方程偏移速度分析

传统的剩余校正（RMO）偏移速度分析方法基于走时原理,在陡倾角和欠照明地区,因为不能得到充分的角度域信息而失效.本文将展示一种基于波场延拓理论的偏移速度分析方法,即波动方程偏移速度分析（WEMVA）.这种方法先利用成像优化方法获得剩余成像,再利用剩余成像反演剩余速度.此类方法继承了波动方程偏移方法的优点和缺点.波动方程偏移速度分析是一种线性反演方法,它要求对Born近似的展开序列作一阶截断.高阶部分的丢失必然带来巨大的截断误差,因此剩余成像必须也进行线性化,以适应大速度扰动和大延拓步长.因此,在此类算法中,剩余成像的获取和线性化是偏移速度分析的关键.在叠前偏移算子中,因为双平方根算子的数学表达式更为简洁,所以本文基于对波动方程偏移速度分析初步讨论,并通过模型验证其原理.     

高斯束逆时偏移方法及应用

随着油气勘探的不断深入,叠前深度偏移技术逐渐成为解决复杂构造区域成像问题的重要手段.考虑到逆时偏移精度高但效率低,而高斯束偏移具有高效、灵活等优势,本文将逆时偏移的理论思想应用于高斯束偏移中,发展了高斯束逆时偏移,即以Kirchhoff偏移为基础,利用高斯束叠加积分计算Green函数.其不仅保留了Kirchhoff偏移的高效性和灵活性,克服了其阴影区、焦散区、单次波至等缺陷,而且还具有波动方程偏移对陡倾角的成像优势,解决了其计算量大的问题.模型和实际资料试算验证了高斯束逆时偏移方法的正确性、有效性和实用性,可以作为一项实用化技术应用于偏移成像和层析速度建模中.     

初至波菲涅尔体地震层析成像

根据地震波传播的有限频理论,对于某个特定震相的观测信息,不仅射线路径上的点对该信息具有影响,射线领域上的其他点对接收信息也具有影响,这种影响可以用核函数来表达.本文基于波动方程的Born近似与Rytov近似,给出了非均匀介质情况下初至地震波振幅与走时菲涅尔体层析成像单频、带限层析核函数的计算方法.通过对均匀介质情况下初至地震波菲涅尔体层析成像核函数解析表达式的理论模型实验与分析,给出了不同维度振幅、走时单频菲涅尔体的空间分布范围,进而给出了带限菲涅尔体边界的确定方法.将本文的走时菲涅尔体层析成像理论应用于表层速度结构反演中,理论模型试验与实际资料处理结果表明,初至波菲涅尔体地震层析成像方法比传统的初至波射线层析成像理论具有更高的反演精度.     

双平方根单程波动方程叠前τ偏移方法

本文将常规双平方根（DSR）单程波动方程从深度域变换到双程垂直走时（τ）域，由此推导出可从数学上实现“沉降观测”的单程波DSR传播算子. 其递归波场延拓算法包含波数域针对常速背景的相移处理和空间域针对横向速度扰动的相位校正，可以应对上覆地层速度横向变化对构造成像的影响. 结合零炮检距、零时间成像条件，提出了在τ域进行波场延拓与成像的DSR方程叠前偏移新方法. 为了克服其全三维偏移算法在实际应用中可能面临的困难，本文采用稳相近似，在crossline常炮检距偏移理论基础上推导了实用的共方位角叠前τ偏移方法. 数值试验表明，DSR方程叠前τ偏移在强横向非均匀介质中的成像精度与分辨率优于传统的时间域成像技术.     

初至波相位走时层析

近地表速度结构通常是利用射线走时层析或菲涅尔体走时层析等反演方法得到的,但它们的目标函数仍利用射线走时残差构建,导致反演精度不高.为此,本文提出了基于散射积分算法的初至波相位走时层析成像方法.该方法的核心是:(1)提出了依赖于频率的相位走时概念;(2)利用依赖于频率的相位走时信息,而非单一的无限频率射线走时;(3)发展了一种改进的相位展开方法,即通过监测相位不连续性和2π周期判定来消除相位折叠现象;(4)考虑了地震波传播的有限频特征,即基于波动理论而非传统的射线路径或有限空间的菲涅尔体构建核函数.通过利用Overthrust模型的数值实验及与传统射线走时层析和菲涅尔体走时层析的对比表明:本文提出的方法是一种有效的初至波走时反演方法.同时,基于Overthrust模型的数值试验还证明了下列结论,即通过挖掘更多的走时信息的确可以获得更高的反演精度和分辨率.     

Wave-equation migration velocity analysis. II. Subsalt imaging examples

Subsalt imaging is strongly dependent on the quality of the velocity model. However, rugose salt bodies complicate wavefield propagation and lead to subsalt multipathing, illumination gaps and shadow zones, which cannot be handled correctly by conventional traveltime‐based migration velocity analysis (MVA). We overcome these limitations by the wave‐equation MVA technique, introduced in a companion paper, and demonstrate the methodology on a realistic synthetic data set simulating a salt‐dome environment and a Gulf of Mexico data set. We model subsalt propagation using wave paths created by one‐way wavefield extrapolation. Those wave paths are much more accurate and robust than broadband rays, since they inherit the frequency dependence and multipathing of the underlying wavefield. We formulate an objective function for optimization in the image space by relating an image perturbation to a perturbation of the velocity model. The image perturbations are defined using linearized prestack residual migration, thus ensuring stability, relative to the first‐order Born approximation assumptions. Synthetic and real data examples demonstrate that wave‐equation MVA is an effective tool for subsalt velocity analysis, even when shadows and illumination gaps are present.     

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.     

Wave equation tomographic velocity inversion method based on the Born/Rytov approximation

This paper discusses Born/Rytov approximation tomographic velocity inversion methods constrained by the Fresnel zone. Calculations of the sensitivity kernel function and traveltime residuals are critical in tomographic velocity inversion. Based on the Born/Rytov approximation of the frequency-domain wave equation, we derive the traveltime sensitivity kernels of the wave equation on the band-limited wave field and simultaneously obtain the traveltime residuals based on the Rytov approximation. In contrast to single-ray tomography, the modified velocity inversion method improves the inversion stability. Tests of the near-surface velocity model and field data prove that the proposed method has higher accuracy and Computational efficiency than ray theory tomography and full waveform inversion methods.     

Numerical Analysis of the Effect of Mesoscale Eddies on Seismic Imaging

The influence of the mesoscale eddy on seismic wave propagation and seismic imaging in deep sea is investigated. Based on fundamental fluid equations, an appropriate partial differential equation is derived for the acoustic pressure field in water with eddies, including current effects. Seismic wavefields and synthetic seismograms in the center of the eddy are simulated. Numerical experiments demonstrate that velocity variations caused by the eddy can lead to traveltime perturbations. Further, in seismic images, the reflectors below the water layer are positioned incorrectly due to the perturbation of the eddy and this image perturbation depends linearly on the migration velocity of the layer below the corresponding reflector. The zero-offset seismic profiling throughout the affected area of the eddy shows that the maximum traveltime perturbation appears at the center of the eddy and the structure of horizontal reflectors below the water layer are distorted.     

Numerical modeling of Gaussian beam propagation and diffraction in inhomogeneous media based on the complex eikonal equation

Gaussian beam is an important complex geometrical optical technology for modeling seismic wave propagation and diffraction in the subsurface with complex geological structure. Current methods for Gaussian beam modeling rely on the dynamic ray tracing and the evanescent wave tracking. However, the dynamic ray tracing method is based on the paraxial ray approximation and the evanescent wave tracking method cannot describe strongly evanescent fields. This leads to inaccuracy of the computed wave fields in the region with a strong inhomogeneous medium. To address this problem, we compute Gaussian beam wave fields using the complex phase by directly solving the complex eikonal equation. In this method, the fast marching method, which is widely used for phase calculation, is combined with Gauss–Newton optimization algorithm to obtain the complex phase at the regular grid points. The main theoretical challenge in combination of this method with Gaussian beam modeling is to address the irregular boundary near the curved central ray. To cope with this challenge, we present the non-uniform finite difference operator and a modified fast marching method. The numerical results confirm the proposed approach.     

有限频率地震层析成像方法及研究进展

本文从基于射线理论的走时地震层析成像的发展引出有限频率地震层析成像, 说明两者最为本质的差别是有限频率灵敏度核函数。 概述了有限频率层析成像方法和基本原理, 详细阐述了国内外有关有限频率层析成像方法的研究进展和应用成果, 强调有限频率层析成像对全球地幔柱研究的重要意义, 最后指出有限频率层析成像方法在理论上的诸多优点目前还没有在实际应用中得到充分体现, 还存在一些问题需要解决。