基于特征能量梯度算子的时间域全波形反演

全波形反演中构建常规梯度算子过程中需要三步骤:震源子波的正向传波场传播,波场残差的反传波场和波场互相关构建梯度算子,其过程存在数据量大、效率低等缺点,为提高反演的效率,本文针对常规时间域梯度算子进行优化,提出了基于特征能量的梯度算法.在正传过程中计算每一个网格点上的正传波场的最大激发能量及其对应的时间步,保存一个子波时间长度的利用特征能量以构建梯度算子.在构建梯度算法中利用保存的子波长度的特征能量进行构建梯度算子.该算法无需保存震源的正传波场,可以减少运算过程中的磁盘读写,提高全波形反演的计算效率.在Mamousi模型梯度测试和实际资料的反演测试中表明:该算法在可以保证梯度算子的精度,具有数据读写量小的优点,效率高的优点.     

起伏地形条件下二维声波频率域全波形反演（英文）

频率域全波形反演充分利用全波场的振幅、相位以及频率信息,采用较少的频率便能反演得到精度很高的速度模型。本文以有限单元法为基础,对起伏地形条件下二维声波频率域全波形反演进行了研究。在正演算法中,针对截断边界问题,并考虑多频率联合反演中计算区域采用同一套剖分网格的需求,提出了一种适用于起伏地形的衰减边界条件算法。该算法的核心思想是在控制方程波数项中引入衰减因子,通过一定方式调节衰减因子使得声波在衰减层中充分衰减,达到压制截断边界影响的目的。根据指数衰减规律,文中推导出了一种新的衰减因子计算公式,并给出了不同频率条件下衰减层厚度计算公式;在反演算法中,采用共轭梯度法求解高斯牛顿反演迭代方程组,避免直接求解雅克比矩阵和HeSSian矩阵带来的巨额计算量,并采用相同的反演模型,对比分析了不同初始模型和频率组合对全波形反演结果的影响。起伏地形模型数值模拟和全波形反演数值试验表明,本文提出的指数衰减边界条件算法和基于该算法的全波形反演算法具有很好的应用效果。     

频率域全波形反演中关于复频率的研究

Laplace-Fourier域全波形反演可以利用简单的初始模型,从缺失低频信息的地震数据中得到长波长速度模型.Laplace-Fourier域全波形反演等价于本文的复频率全波形反演,但二者的实现方式不同,因此研究复频率全波形反演,可以为二者的对比研究并发展更有效的方法奠定重要基础.本文首先比较用线性增加模型作为初始模型时几个包含不同高低频成分的频率组的反演效果,再比较结合复频率之后各个频率组的反演效果,从简单模型和复杂模型的测试中都可以看出这种复频率+频率反演的方式对反演效果有明显改善.     

基于混合自适应遗传算法的稳健全波形反演

地震全波形反演理论与技术虽已得到了广泛研究,但周波跳跃等问题的存在严重制约了常规全波形反演方法的实用化进程.基于遗传算法的全波形反演方法能够在一定程度上较好地缓解常规全波形反演面临的初始模型依赖性问题,但是当前方法仍存在收敛性和巨大的计算量问题.本文提出一种混合自适应遗传算法(HAGA),并提出基于HAGA的稳健全波形反演方法,该方法将基于HAGA的反演与基于共轭梯度法的常规全波形反演交替迭代进行,其可兼顾反演计算效率与精度.数值测试结果表明,局部与全局优化交替迭代的全波形反演方法集合了局部优化反演的高效与全局优化反演的稳定的优点,大大降低了全波形反演对初始模型的依赖性,能够有效的缓解常规全波形反演的周波跳跃问题.     

基于常分数阶拉普拉斯算子的黏声方程重建速度与衰减参数的全波形反演方法

地震波在地下介质传播过程中由于非弹性衰减的存在将导致能量损失和相位变化,精确的速度与衰减参数建模对油气识别、提高强衰减介质中地震波成像的质量都起着至关重要的作用.常分数阶拉普拉斯算子黏声方程由于完全分离的速度频散项与振幅衰减项的优势,以及在强非均质衰减介质中可以高精度求解的特点,已被应用于速度与衰减参数的建模中.本文将二阶常分数阶拉普拉斯算子黏声方程拆分为等价的一阶方程组,并在此一阶方程组的基础上推导出新的梯度公式与伴随方程,建立了一种新的速度与衰减参数同时重建的全波形反演方法.相较于原二阶常分数阶拉普拉斯算子黏声方程建立的全波形反演流程,数值实验表明,新建立的反演流程可以有效避免原梯度数值计算中的噪声,尤其是可以有效提高衰减参数梯度的反演精度,从而显著提高反演的收敛速度与反演精度.     

基于源波场重构的子波自适应频移算子

地震信号中的多尺度信息对于分辨率、成像精度和反演结果有非常重要的意义,本文提出一种基于直达波模拟波场重建震源子波的地震数据频移算子,以期能应用于基于波动的地震信号多尺度分频.频移过程主要分为两步,第一步利用直达波反传构建震源子波,第二步借助震源子波和频移算子进行多尺度分频.与常规数字滤波器相比,频移算子突破了滤波造成信号波形特征改变和震源子波形态畸变的限制,频移地震数据与数值模拟地震数据完美匹配,同时频移算子具有理论子波自适应功能,更易于后续波动类应用展开.最后通过层状模型和实际资料进行测试,并与常规滤波器结果对比证明本方法的准确性和稳定性.     

基于反射地震数据的时频域包络反演

包络信号含有丰富的低频分量,即使在地震数据缺失低频条件下,包络目标函数也能有效缓解全波形反演的周期跳跃现象.但是,当初始速度模型较为平滑时,观测数据中的反射地震事件在模拟数据中没有与之相对应的波形,导致包络反演初期无法很好地利用反射波信号进行速度建模.本文提出基于反射地震数据的时频域包络反演方法,通过结合反射波全波形反演理论,构建反射波时频域包络目标函数,来提高包络反演的速度建模精度.本文首先利用Gabor变换获取时频域地震数据,并提取振幅信息,即为时频域包络信号.然后,推导反射波时频域包络反演的伴随震源和梯度算子.Marmousi模型数据测试结果表明,基于反射地震数据的时频域包络反演方法可以为全波形反演提供一个较好的初始速度模型.     

基于反射地震数据的时频域包络反演

包络信号含有丰富的低频分量,即使在地震数据缺失低频条件下,包络目标函数也能有效缓解全波形反演的周期跳跃现象.但是,当初始速度模型较为平滑时,观测数据中的反射地震事件在模拟数据中没有与之相对应的波形,导致包络反演初期无法很好地利用反射波信号进行速度建模.本文提出基于反射地震数据的时频域包络反演方法,通过结合反射波全波形反演理论,构建反射波时频域包络目标函数,来提高包络反演的速度建模精度.本文首先利用Gabor变换获取时频域地震数据,并提取振幅信息,即为时频域包络信号.然后,推导反射波时频域包络反演的伴随震源和梯度算子.Marmousi模型数据测试结果表明,基于反射地震数据的时频域包络反演方法可以为全波形反演提供一个较好的初始速度模型.     

伪二维弹性波联合反演近地表的速度和衰减

利用弹性波的初至波和面波,应用交叉梯度算子,联合反演了近地表的二维纵横波速度和衰减参数,并提出了采用一维弹性波正演模拟,应用二维Tikhonov正则化,同时反演出二维速度模型和衰减模型的方法。理论模型测试和实际数据应用结果均表明本文算法极大地提高了计算效率,同时能够反演出可靠的速度模型和衰减模型。      

基于Hinge损失函数的垂向全变差约束全波形反演

全波形反演可以为叠前深度偏移成像提供更高精度的速度模型,但该方法具有较强的非线性,对初始速度模型的依赖性较强,尤其是在实际应用中,地质条件复杂多变,速度变化不连续,增加了反演非线性程度,常常使反演陷入局部极小值,影响反演的精度.全变差约束在图像去噪领域应用广泛,属于非光滑约束,在去噪过程中能有效的保留图像的不连续界面和边缘信息.本文提出基于Hinge损失函数的垂向全变差约束全波形反演方法,在全变差约束的基础上,利用Hinge损失函数控制模型的更新方向,并使用原-对偶混合梯度算法进行求解,给出这一优化问题的迭代格式,有效提高了对地下不连续界面的重构精度,同时也降低反演对初始速度模型的依赖程度.数值算例证明：与常规全波形反演方法相比,基于全变差约束的全波形反演方法可以有效的重构速度模型中的不连续界面,尤其对高速体边缘的重构效果更明显,但该方法对初始速度模型的依赖性仍然较强;基于Hinge损失函数的垂向全变差约束全波形反演方法降低了对初始速度模型的依赖程度,可以从一个较差的初始模型通过循环迭代的方式最终得到同样精确的速度模型,较好的重构了高速体边缘和不连续界面.     

Suppressing non‐Gaussian noises with scaled receiver wavefield for reverse‐time migration: comparison of different approaches

Numerical implementation of the gradient of the cost function in a gradient‐based full‐ waveform inversion (FWI) is essentially a migration operator used in wave equation migration. In FWI, minimizing different data residual norms results in different weighting strategies of data residuals at receiver locations prior to back‐propagation into the medium. In this paper, we propose different scaling methods to the receiver wavefield and compare their performances. Using time‐domain reverse‐time migration (RTM), we show that compared to conventional algorithms, this type of scaling is able to significantly suppress non‐Gaussian noise, i.e., outliers. Our tests also show that scaling by its absolute norm produces better results than other approaches.     

基于各向异性全变分约束的多震源弹性波全波形反演

多震源编码技术可以提高全波形反演的计算效率,但同时会引入串扰噪声使反演结果质量降低. 全变分约束可以有效地压制层内噪声并突出模型界面,其与多震源技术的结合,能在大大提高弹性波全波形反演效率的同时提高反演质量. 本文提出了一种高效的动态多震源全波形反演策略,可以在离散串扰噪声的同时保证照明的均匀性. 根据残留串扰噪声的分布特征,构建了与之匹配的基于各向异性全变分约束的弹性波全波形反演方法. 为了减少周期跳跃效应,将基于稀疏约束的低频重构算法应用于弹性波地震记录,给出了利用快速梯度投影算法求解各向异性全变分约束的全波形反演流程. 模型数据测试结果表明本文的方法不仅能有效地抑制多震源方法导致的串扰噪声,同时也能降低观测数据中的噪声对反演结果的影响.     

Seismic envelope inversion: reduction of local minima and noise resistance

Waveform inversion met severe challenge in retrieving long‐wavelength background structure. We have proposed to use envelope inversion to recover the large‐scale component of the model. Using the large‐scale background recovered by envelope inversion as new starting model, we can get much better result than the conventional full waveform inversion. By comparing the configurations of the misfit functional between the envelope inversion and the conventional waveform inversion, we show that envelope inversion can greatly reduce the local minimum problem. The combination of envelope inversion and waveform inversion can deliver more faithful and accurate final result with almost no extra computation cost compared to the conventional full waveform inversion. We also tested the noise resistance ability of envelope inversion to Gaussian noise and seismic interference noise. The results showed that envelope inversion is insensitive to Gaussian noise and, to a certain extent, insensitive to seismic interference noise. This indicates the robustness of this method and its potential use for noisy data.     

Full waveform inversion and the truncated Newton method: quantitative imaging of complex subsurface structures

Full waveform inversion is a powerful tool for quantitative seismic imaging from wide‐azimuth seismic data. The method is based on the minimization of the misfit between observed and simulated data. This amounts to the solution of a large‐scale nonlinear minimization problem. The inverse Hessian operator plays a crucial role in this reconstruction process. Accounting accurately for the effect of this operator within the minimization scheme should correct for illumination deficits, restore the amplitude of the subsurface parameters, and help to remove artefacts generated by energetic multiple reflections. Conventional minimization methods (nonlinear conjugate gradient, quasi‐Newton methods) only roughly approximate the effect of this operator. In this study, we are interested in the truncated Newton minimization method. These methods are based on the computation of the model update through a matrix‐free conjugate gradient solution of the Newton linear system. We present a feasible implementation of this method for the full waveform inversion problem, based on a second‐order adjoint state formulation for the computation of Hessian‐vector products. We compare this method with conventional methods within the context of 2D acoustic frequency full waveform inversion for the reconstruction of P‐wave velocity models. Two test cases are investigated. The first is the synthetic BP 2004 model, representative of the Gulf of Mexico geology with high velocity contrasts associated with the presence of salt structures. The second is a 2D real data‐set from the Valhall oil field in North sea. Although, from a computational cost point of view, the truncated Newton method appears to be more expensive than conventional optimization algorithms, the results emphasize its increased robustness. A better reconstruction of the P‐wave velocity model is provided when energetic multiple reflections make it difficult to interpret the seismic data. A better trade‐off between regularization and resolution is obtained when noise contamination of the data requires one to regularize the solution of the inverse problem.     

Correlation‐based reflection waveform inversion by one‐way wave equations

Reflection full waveform inversion can update subsurface velocity structure of the deeper part, but tends to get stuck in the local minima associated with the waveform misfit function. These local minima cause cycle skipping if the initial background velocity model is far from the true model. Since conventional reflection full waveform inversion using two‐way wave equation in time domain is computationally expensive and consumes a large amount of memory, we implement a correlation‐based reflection waveform inversion using one‐way wave equations to retrieve the background velocity. In this method, one‐way wave equations are used for the seismic wave forward modelling, migration/de‐migration and the gradient computation of objective function in frequency domain. Compared with the method using two‐way wave equation, the proposed method benefits from the lower computational cost of one‐way wave equations without significant accuracy reduction in the cases without steep dips. It also largely reduces the memory requirement by an order of magnitude than implementation using two‐way wave equation both for two‐ and three‐dimensional situations. Through numerical analysis, we also find that one‐way wave equations can better construct the low wavenumber reflection wavepath without producing high‐amplitude short‐wavelength components near the image points in the reflection full waveform inversion gradient. Synthetic test and real data application show that the proposed method efficiently updates the background velocity model.     

一种基于高效迭代解法的频率域全波形反演

巨大的计算量是制约全波形反演（FWI）生产实用化的难题之一.为此,本文提出了一种高效的波场迭代解法,将其应用于频率域常密度声波方程FWI,并给出了详细的反演流程.通过建立用于波场迭代的目标函数,推导相应梯度、步长公式,新方法将反演中波场正传和残差波场反传过程转化为无约束优化问题,从理论上分析了新方法的计算效率显著高于常规FWI.在数值试验中,本文方法通过几次迭代便能获得高精度的正传、残差反传波场,收敛速度明显高于未经预处理的GMRES方法.进一步引入高效编码策略,新方法的计算时间约为常规编码FWI的1/8,与理论分析结果吻合（波场迭代次数为8,模型未知量个数约为7万）,且波场迭代次数为6时,反演效果已与常规编码FWI相近.     

时间二阶积分波场的全波形反演

通过对波场的时间二阶积分运算以增强地震数据中的低频成分,提出了一种可有效减小对初始速度模型依赖性的地震数据全波形反演方法—时间二阶积分波场的全波形反演方法.根据散射理论中的散射波场传播方程,推导出时间二阶积分散射波场的传播方程,再利用一阶Born近似对时间二阶积分散射波场传播方程进行线性化.在时间二阶积分散射波场传播方程的基础上,利用散射波场反演地下散射源分布,再利用波场模拟的方法构建地下入射波场,然后根据时间二阶积分散射波场线性传播方程中散射波场与入射波场、速度扰动间的线性关系,应用类似偏移成像的公式得到速度扰动的估计,以此建立时间二阶积分波场的全波形迭代反演方法.最后把时间二阶积分波场的全波形反演结果作为常规全波形反演的初始模型可有效地减小地震波场全波形反演对初始模型的依赖性.应用于Marmousi模型的全频带合成数据和缺失4Hz以下频谱成分的缺低频合成数据验证所提出的全波形反演方法的正确性和有效性,数值试验显示缺失4Hz以下频谱成分数据的反演结果与全频带数据的反演结果没有明显差异.     

基于混叠震源和频率组编码的频率域自适应全波形反演(英文)

全波形反演是一种高精度的地震成像方法,可以对地下介质物性参数模型进行准确的重构。然而在实际应用中,尤其是在三维复杂介质反演中,计算成本太大是该方法的一个重要缺陷。将混叠震源技术引入到频率域全波形反演中可以大幅度地降低计算成本,提高反演效率。但是使用震源编码技术也带来了两个问题:一方面,参与编码的各个震源之间会产生"串扰噪声",导致反演结果中出现假象;另一方面,基于震源编码的频率域全波形反演方法周围噪声较为敏感,使该方法对含噪数据反演质量较差。本文引入一种频率组编码方法来压制"串扰噪声",并基于震源编码技术提出一种频率域自适应全波形反演方法,通过一个与频率相关的自适应选择机制,将常规频率域全波形反演方法和基于震源编码的全波形反演方法联合起来,在保证反演质量的同时也最大程度地提高了反演效率。     

2D acoustic‐elastic coupled waveform inversion in the Laplace domain

Although waveform inversion has been intensively studied in an effort to properly delineate the Earth's structures since the early 1980s, most of the time‐ and frequency‐domain waveform inversion algorithms still have critical limitations in their applications to field data. This may be attributed to the highly non‐linear objective function and the unreliable low‐frequency components. To overcome the weaknesses of conventional waveform inversion algorithms, the acoustic Laplace‐domain waveform inversion has been proposed. The Laplace‐domain waveform inversion has been known to provide a long‐wavelength velocity model even for field data, which may be because it employs the zero‐frequency component of the damped wavefield and a well‐behaved logarithmic objective function. However, its applications have been confined to 2D acoustic media. We extend the Laplace‐domain waveform inversion algorithm to a 2D acoustic‐elastic coupled medium, which is encountered in marine exploration environments. In 2D acoustic‐elastic coupled media, the Laplace‐domain pressures behave differently from those of 2D acoustic media, although the overall features are similar to each other. The main differences are that the pressure wavefields for acoustic‐elastic coupled media show negative values even for simple geological structures unlike in acoustic media, when the Laplace damping constant is small and the water depth is shallow. The negative values may result from more complicated wave propagation in elastic media and at fluid‐solid interfaces. Our Laplace‐domain waveform inversion algorithm is also based on the finite‐element method and logarithmic wavefields. To compute gradient direction, we apply the back‐propagation technique. Under the assumption that density is fixed, P‐ and S‐wave velocity models are inverted from the pressure data. We applied our inversion algorithm to the SEG/EAGE salt model and the numerical results showed that the Laplace‐domain waveform inversion successfully recovers the long‐wavelength structures of the P‐ and S‐wave velocity models from the noise‐free data. The models inverted by the Laplace‐domain waveform inversion were able to be successfully used as initial models in the subsequent frequency‐domain waveform inversion, which is performed to describe the short‐wavelength structures of the true models.