螺旋坐标下的因子分解合成炮叠前深度偏移

在合成炮叠前深度偏移的基础上，提出了一种高效的合成炮叠前深度偏移方法，即在螺旋坐标下用因子分解进行波场外推的混合法. 用因子分解进行波场外推分为因果过程和反因果过程两个显式求解过程. 这种螺旋坐标下的显式求解过程,提高了波场外推的效率. 根据相位编码原理,对多个射线参数的合成波场进行了编码叠加,基于射线参数实现了MPI并行计算,进一步提高了计算效率. 在推导了有关公式并进行定量分析之后,对Marmousi复杂模型进行了计算,并作了比较,结果表明本文方法具有精度高和速度快的特点,可用于实际资料的计算.     

高精度混合法叠前深度偏移及其并行实现

叠前深度偏移是复杂构造成像的重要手段.本文基于波场分裂理论,首先给出了波场延拓的一般方程,即一个上下行波的耦合方程组.通常所用的波场延拓方程就是该耦合方程组的特例.根据平方根算子的近似,推导了一种新的高精度混合法偏移方法,用分裂法即可进行计算.通过对一个横向强烈变速模型的叠后偏移及Marmousi复杂构造模型的叠前偏移计算,说明了方法的有效性,精度较高.采用MPI并行编程实现并行计算,提高了计算效率.该方法可用于对横向强烈变速的复杂构造的精确成像.     

基于CLG光学流和波场分解的逆时偏移角度道集提取方法研究

角度域共成像点道集是衔接叠前地震数据与储层特征的重要桥梁,对地震偏移成像与储层描述具有重要意义.与克希霍夫偏移和单向波动方程偏移相比,逆时偏移是复杂地区最精确的成像方法.高效稳健地生成逆时偏移角度道集目前仍然是一个挑战.本文主要讨论如何采用光学流方法高效、高质量地提取角度道集.在逆时偏移波场外推过程中,光学流方法可以估计波场传播方向.其中Lucas-Kanade(LK)和Horn-Schunck(HS)方法是光学流方法中两种典型的方法.LK光学流方法是一种局部方法,该方法依赖于局部点的梯度值,但是容易出现奇异现象,HS光学流方法属于全局方法,波场方向估计依赖于整个波场,易受噪声影响,对异常值比较敏感,导致整体波场方向计算精度不高.本文提出采用局部和整体结合(Combining Local and Global,CLG)的光学流方法估计波场传播方向.该方法可以有效地提高波场方向的精度,并且简单高效,便于并行处理.对比HS光学流方法,CLG光学流方法几乎不增加额外的计算量.另外,为了弱化光学流方法无法处理波前重叠问题,本文利用解析波场和方向滤波对波场进行方向分解,仅需波场的空间傅里叶变换即可实现任意波场方向分解,将分解后的波场分别估计波场反向,提取成像结果.进一步地,在估计反射张角和方位角时,本文提出有效的归一化方法和改进的最小二乘除法,提高角度估计的精度和稳定性.最后,理论和实际资料例证了本文提出方法的有效性.     

基于格林函数理论的波场预测和鬼波压制方法

作为一种特殊的噪声,鬼波对一次波的波形及频带宽度产生极大的影响,鬼波压制是提高海上地震资料分辨率及保真度的重要因素.以格林公式为基础,详细论述了基于格林函数理论的鬼波压制方法,在不需要地下介质信息的条件下,进行地震数据驱动鬼波压制,并根据"Double Dirichlet"(双狄利克雷)边界条件,预测压力波场和垂直速度波场.建立了基于格林函数理论鬼波压制的处理流程,数值模拟和实际资料处理结果表明,基于格林函数理论鬼波压制方法在很好地去除鬼波的同时极大地拓宽了地震资料的频带,尤其提升了低频端能量,有利于后续资料的处理解释.     

裂缝诱导双相HTI介质地震波场错格伪谱法模拟与波场特征分析

裂缝诱导的双相具有水平对称轴的横向各向同性(HTI)介质模型是由一组平行排列的垂直裂缝嵌入到统计各向同性的流体饱和多孔隙岩石中而组成的,它综合考虑了裂缝型储层岩石的各向异性和孔隙性.高精度的地震波场数值模拟技术是研究该介质中地震波传播规律的主要方法.本文结合错格伪谱法和时间分裂法,求解描述该介质中地震波传播的一阶速度-应力方程.模拟了单层和双层模型中的地震波场,并对其进行了特征分析.研究结果表明:错格伪谱法能有效消除标准网格伪谱法波场模拟结果中出现的数值伪影现象,与时间分裂法结合能够获得稳定的、高精度的模拟结果;裂缝诱导双相HTI介质中的地震波场兼具裂缝各向异性介质和双相介质中传播的地震波的波场特征.     

瞬变电磁虚拟波场的三维曲面延拓成像

根据瞬变电磁场扩散方程与波动方程间存在的数学对应形式,将已知瞬变电磁场数据转换为虚拟波场数据,在此基础上,提出了用克希霍夫积分法将所转换的波场从地面向地下反向延拓的方法.在延拓的数值计算中,为提高计算效率,采用三维边界元技术,把边界积分分解为诸三角单元积分的积分之和, 建立了曲面延拓方程式;通过地面测点的波场值,求出地下某一点的波场值,实现了地表为曲面的向下延拓成像计算.通过对理论模型计算和实际资料处理,证明了该方法可以增强瞬变电磁法识别地下电性分界面的能力,使瞬变电磁法对地下目标体的三维精细探测成为可能.     

共炮检距二维弹性波地震剖面的偏移方法

在均匀各向同性介质的分界面上入射纵波时，除产生纵波反射Ｐ－Ｐ波外，还会产生转换横波的反射Ｐ－ＳＶ波。同样，当入射波为ＳＶ横波时，在分界面上会同时产生ＳＶ－ＳＶ反射波和ＳＶ－Ｐ反射转换波，对这种二维弹性波地震记录如何进行偏移是正在探讨的问题。为了能够处理实际地震数据，本文提出在共炮检距剖面上进行二维弹性波记录的偏移方法。具体做法是，首先用弹性波方程将两分量的地震记录分解为Ｐ－Ｐ（或ＳＶ－ＳＶ）波和Ｐ－ＳＶ（或ＳＶ－Ｐ）波，然后用我们已导出的适合于处理共炮检距剖面的叠前偏移方法对它们进行处理。用模拟数据进行的试验说明方法正确，效果良好。     

The design of visco‐acoustic weighted L1‐error frequency–space wavefield extrapolators

Seismic imaging is an important step for imaging the subsurface structures of the Earth. One of the attractive domains for seismic imaging is explicit frequency–space (f – x) prestack depth migration. So far, this domain focused on migrating seismic data in acoustic media, but very little work assumed visco‐acoustic media. In reality, seismic exploration data amplitudes suffer from attenuation. To tackle the problem of attenuation, new operators are required, which compensates for it. We propose the weighted L ₁ ‐error minimisation technique to design visco‐acoustic f – x wavefield extrapolators. The L ₁ ‐error wavenumber responses provide superior extrapolator designs as compared with the previously designed equiripple L ₄ ‐norm and L_∞‐norm extrapolation wavenumber responses. To verify the new compensating designs, prestack depth migration is performed on the challenging Marmousi model dataset. A reference migrated section is obtained using non‐compensating f – x extrapolators on an acoustic dataset. Then, both compensating and non‐compensating extrapolators are applied to a visco‐acoustic dataset, and both migrated sections are then compared. The final images show that the proposed weighted L ₁ ‐error method enhances the resolution and results in practically stable images.     

基于一步法波场延拓的正演模拟和逆时偏移成像

通常工业界实现逆时偏移算法时采用有限差分数值方法模拟地震波场,波场模拟常常受稳定性条件限制,且易产生数值频散,成像精度降低.本文引入了一步法波场延拓方法,首先构建声波传播算子,借助Chebyshev多项式和Jacobi-Anger展开式近似传播算子中的e指数项,进而实现波场递推,该方法时间步长的选取不受稳定性条件限制而且不存在空间频散现象.本文将一步法波场延拓方法用于逆时偏移成像的波场模拟,并提出双缓冲区存储策略,在不增加计算量的前提下,大幅降低了逆时偏移方法的波场存储量.波场模拟和逆时偏移成像测试表明,本文提出的一步法波场延拓方法模拟地震波场精度高,消除了频散影响,可在较大时间步长的情况下实现高精度波场模拟;提出的基于一步法波场延拓的逆时偏移方法成像质量好;基于双缓冲区存储策略的逆时偏移成像方法存储成本低.     

基于新的惩罚因子算法的波场重构反演

全波形反演是一种高精度的反演方法,其目标函数是一个强非线性函数,易受局部极值影响,而且反演过程计算量较大.波场重构反演是近几年提出的一种改进的全波形反演理论.该反演方法通过将波动方程作为惩罚项引入到目标函数中,通过拓宽解的寻找空间减弱了局部极小值的影响,而且反演过程不需要计算伴随波场,提高了计算效率.但该反演方法一直缺少准确的惩罚因子算法,直接影响到该方法的准确度.本文将波场重构反演拓展到时间域并利用梯度法进行波场重构.频率域的惩罚因子用来加强波动方程的约束,而时间域惩罚因子表现为调节模拟波场和实际波场的权重因子.为此,我们根据约束优化理论,在波动方程准确以及重构波场与反演参数解耦的假设下,提出以波动方程为目标函数的新的惩罚因子算法.根据波形反演在应用时普遍存在的噪音干扰、子波错误和低频信息缺失的情况下,应用部分Sigsbee2A模型合成数据对本文提出的算法进行实验.数值实验结果表明：基于新的惩罚因子算法,在其他信息不准确的情况下,波场重构反演可以给出高精度的反演结果.