双程声波方程逆时深度偏移

从双程声波方程出发,在交错网格空间中推导了地震波逆时延拓的高阶有限差分算子,依据最佳匹配层(PML)的方程分裂思路,得到了一阶声波方程的PML边界条件及其高阶差分格式,采用零时间成像条件和上行、下行波场互相关成像条件,实现了声波方程的叠后与叠前逆时深度偏移。逆时偏移对sigsbee_2b模型理论数据的偏移成像得到了满意效果。     

各向异性介质拟声波叠前逆时偏移及成像条件的应用

考虑到基于各向同性假设的传统逆时偏移方法的局限,研究了各向异性介质叠前逆时偏移方法。由于各向异性弹性波逆时偏移中,存在着计算效率低和波场分离困难等不足,各向异性介质拟声波逆时偏移更具优势。从声学近似的胡克定律和运动方程出发,研究了VTI介质拟声波波动方程,并引入旋转坐标系将其推广为TTI介质情形。然后利用交错网格高阶有限差分方法,构建逆时偏移所需的正向和逆时波场外推算子,并引入完全匹配层(PML)边界条件来吸收边界反射能量。最后应用优化的归一化互相关成像条件,实现各向异性拟声波逆时偏移成像。模型试算结果表明,考虑各向异性因素的逆时偏移方法能够有效地改善成像质量,应用优化的成像条件得到的成像效果,明显优于常规的归一化互相关成像条件。     

叠前弹性波逆时深度偏移及波场分离技术探讨

这里探讨了Sun和McMechan提出的叠前弹性波标量逆时深度偏移方法,即在地表附近,对地表接收到的弹性波波场分量进行波场逆时延拓,然后分别对在波场延拓过程中通过波场分离得到的纯纵、横波反射波场进行偏移。在模型中,每个网格点的纵、横波成像时间,为震源到该网格点的纵波初至时间。利用纵波速度模型和声波传播方程的有限差分解,对反射纵波进行逆时延拓和成像,利用横波速度模型和声波传播方程的有限差分解,对转换横波进行逆时延拓和成像,并在偏移前进行极性校正。然后,在Sun和McMechan偏移成像方法的基础上,提出了直接偏移成像方法,即对地表接收到的弹性波波场分量在整个模型范围内进行波场逆时延拓。在波场延拓过程中,对符合成像条件的网格点进行波场分离、偏移成像和转换波极性校正。最后,通过数值模拟对二种偏移成像方法的结果进行了比较。     

任意广角方程逆时偏移的脉冲响应及模型试算

采用高阶差分法求解任意广角声波方程,推导了各向同性介质中任意广角声波方程叠前逆时深度偏移的差分算子,给出了适用于任意广角声波方程逆时延拓的吸收边界条件和其计算方法,在此基础上实现了任意广角声波方程的高精度逆时延拓。分析了不同模型条件下任意广角声波方程叠前逆时深度偏移的脉冲响应,分析结果表明,任意广角方程逆时偏移算子可压制层间反射波,减少偏移噪声,理论模型的成像结果表明该方法能够对复杂构造模型精确成像。     

逆时偏移双向照明波场分离成像条件

通过对波场分离互相关成像条件与炮-检点双向照明互相关成像条件进行研究,将两种方法在噪声消除和补偿成像振幅畸变方面的优势进行结合,提高成像分辨率,改进现有成像条件。其主要步骤是将震源和检波点波场分离成它们的单程波传播分量,采用双向照明补偿互相关成像条件对需要的波场进行成像。通过对模型试算,验证了笔者提出的成像条件能有效的压制噪音,补偿成像振幅。     

基于屏方法的弹性波叠前深度偏移

弹性波成像是多分量地震勘探需要解决的关键技术之一。借鉴声波方程中傅立叶有限差分算子中的有限差分补偿方法,对弹性波相屏传播算子进行了大角度校正。在求取成像值时,采用了矢量成像条件,即通过极化矢量,将各分量投影极化方向,然后再成像,从而使波场中的信息得到充分利用,同时解决了P-SV波成像过程中的极性不一致问题。这样,便可以叠加多炮的成像结果,从而提高成像质量。利用模型数据,对大角度相位校正后的弹性屏延拓算子和矢量成像条件进行了试验,成像效果理想。      

复杂地表地震资料叠前深度偏移成像

分析了复杂地表地区复杂地震波场的形成原因,探讨了基于高角度时空域单程波动方程适用于起伏地表之下构造成像的逆时叠前深度偏移方法。从声波方程有限差分法合成起伏地表上的炮集记录知,当激发点位于低速层中,共炮点道集中低速层表面和高速层表面接收到的信号强度相差较大:低速层表面的信号很强,而高速层表面的信号很弱,低速层中产生了较强的槽波。将该逆时叠前深度偏移方法应用于起伏地表地震波场的偏移处理。虽然在偏移成像前既没有压制与低速层有关的槽波,又没有压制随机噪声,但偏移剖面上界面清晰、位置正确,断层面也得到了很好的聚焦。     

地震波逆时偏移中两种成像条件应用效果对比

成像条件是实现地震波场准确成像的关键环节,而优化逆时成像条件在改善逆时成像应用效果中的作用尤为明显。为此,从零井源距和非零井源距脉冲响应出发,系统对比了常规相关法逆时成像条件和基于行波分离法逆时成像条件的逆时偏移脉冲响应波场差异,分析了逆时偏移低波数噪声的成因机理及波场规律,同时以复杂断陷理论模型和SZ工区实际三维地震资料为例,系统对比了两种逆时成像条件的成像应用效果。研究结果表明,采用基于行波分离法逆时成像条件的成像结果能够有效衰减逆时偏移低波数背景噪声,恢复掩盖在这种背景噪声之下的有效地震成像信息,地层细节刻画更加清晰,这可为当前复杂波场和复杂构造高精度地震成像提供方法指导。     

基于双相介质弹性波波场分离的逆时偏移方法

随着地震勘探的发展,以单相介质模型为基础的声波逆时偏移和弹性波逆时偏移很难比较精确地描述地下介质和满足复杂的工程要求,尤其是在处理地下含流体多孔介质时,所以需要引入双相介质模型。首先基于Biot理论双相介质一阶速度-应力弹性波方程,利用散度场和旋度场的波场分离理论将固相和流相中的纵横波场解耦分离,并由波场快照图得出固相与流相中的纵波波场与横波波场均完成分离,而快纵波与慢纵波相互耦合均存在于纵波波场中。采用高阶交错网格有限差分方法,结合PML吸收边界条件,推导了具有波场分离的双相介质弹性波正传与反传的偏移延拓算子,结合互相关成像条件,可以得到双相介质弹性波固相与流相中的各分量偏移成像剖面。最后对模型使用该双相介质弹性波波场分离的逆时偏移方法进行偏移成像试算并取得成功,效果较常规双相介质弹性波逆时偏移成像效果有了明显提升。     

减小概率成像分辨宽度的一种方法

概率成像方法是利用测量的地面电场与电场扫描函数进行归一化的互相关运算,从而对地下电荷进行定位。为了提高这种方法的分辨率,从而提高对电荷的定位精度,根据正点电荷模型的电位剖面分布,提出了概率成像分辨宽度的概念,根据这个概念,提出了一种减小分辨宽度的方法-功率谱法,即由原先的场函数的归一化互相关运算改为场强平方后的归一化互相关运算。通过对算例的正反演研究,证明了确实是一种较好的改进方式。     

弹性波逆时偏移及其成像条件

从弹性波动方程出发,推导了二维各向同性介质情况下弹性波逆时传播的高阶差分格式,实现了弹性波在数值空间中的逆时延拓。从程函方程出发,采用逆时差分格式求取网格空间中各点的直达波旅行时,以此作为弹性波逆时偏移的成像条件(激发时间成像条件),实现多波多分量资料的逆时偏移。文中给出了几个典型模型的理论记录的偏移结果,数值试验表明,逆时偏移能够使地表接收到的波场准确归位,提高资料处理的精度。     

Stretching correction for amplitude-preserving vector wavefield reverse-time migration

The migration of multi-wave seismic data is aimed at obtaining the P- and S-wave imaging results of the amplitude preserving. But the P- and S-wave stretching effect produced by the reverse time migration of the elastic wave equation will not only reduce the vertical resolution of the migration results and the amplitude preserving of the large reflection angle. In this paper, the reverse time migration technique of amplitude preserving vector wave-field separating is used. Based on the analysis of the stretch mechanism and the influencing factors of stretch magnitude, the paper gave the stretch correcting factors. Then, realize the stretch correction method at the time that after the reverse extrapolation and before the imaging by solving the problem which is how to calculate the P-wave and Ps-wave propagation directions of imaging points at different times. The stretch correction method can improve the vertical resolution and amplitude fidelity of the imaging results and provide high fidelity input data for seismic data interpretation and inversion.© 2019 China Geology Editorial Office.     

偏移速度分析与偏移成像

偏移速度分析和偏移成像是地震资料处理的两个重要组成部分.目前时间偏移技术比较成熟,而深度偏移技术也在逐步完善,时间域偏移成像主要推崇叠前时间偏移法.采用沿层叠加速度分析技术可获得层面上准确的叠加速度,并通过倾角校正、叠前时间偏移和CRP反偏移速度分析逐步优化速度,得到一个符合地质规律、准确的均方根速度场;通过深度偏移方法的研究,总结了建立精确偏移速度场的方法,并提出了一种地震资料处理的思路,即基于射线追踪的Kirchhoff偏移和基于波场延拓的波动方程偏移的结合,使偏移速度分析和偏移成像在应用效果和效率上得到了很大的提高.     

随机介质中的探地雷达叠前逆时偏移成像

为了更好地分析探地雷达(GPR)叠前逆时偏移算法对实际地下随机分布介质的成像效果,采用随机过程的谱分解和混合型自相关函数理论构建了不同自相关长度的GPR随机介质模型。基于时域有限差分法(FDTD)构建了GPR叠前逆时偏移成像算法,其中FDTD用于计算正传和反传电磁波场,归一化互相关成像条件用于获取逆时偏移成像剖面。在此基础上,利用该算法对两个随机介质模型的多偏移距正演数据进行计算,并与相应背景介质为均匀介质的逆时偏移结果进行对比。结果表明:电磁波在随机介质中散射强烈,反射波扭曲变形、不连续,形成了明显的随机扰动,致使逆时偏移成像剖面的空间分辨率更低,低频噪声更强;自相关长度是影响随机介质中异常体成像效果的主要因素,自相关长度越小,异常体的成像越清晰,自相关长度越大,异常体的成像效果越差,且不易被识别。     

基于保幅波动方程的广义高阶屏地震偏移成像方法

原始地震数据中饱含丰富的反映相位的走时信息和反映反射率的振幅信息。基于波场延拓理论的波动方程保幅地震偏移成像,是在给出正确构造成像位置的同时也给出真实反射振幅的有效完善。基于全声波方程,利用严格的解耦理论进行单程波动保幅分解,得到一个由波场传播项与振幅补偿项构成的,在走时与振幅上满足全声波方程对应的程函方程与输运方程的保幅单程波动方程;利用摄动理论进行单平方根算子渐进展开,推导出基于保幅波动方程的广义高阶屏地震偏移算子方程。模型测试和实际资料处理表明,该方法不但可以凭借更准确的相位归位和散射能量聚焦提高构造成像精度,而且输出了能更正确反映地下反射属性的能量信息,从而可以为更深层次的勘探开发,提供地球物理技术支撑。     

二维声波方程稳健迭代速度反演

从二维声波方程初、边值问题出发,通过引进一个广义幂指数误差分布函数,以及对波动方程的Ｌｉｐｐｍａｎ－Ｓｃｈｗｉｎｇｅｒ方程的解进行Ｂｏｎ近似,导出一个残差加权迭代最小二乘（稳健迭代）算法,进而实现对地下介质剖面的速度结构以演计算。模拟计算表明,本方法具有较高的精度和较强的抗干扰能力,是求解声波方程反演问题的一种有效方法,同时也为地下介质的速度分布成像提供一种新的技术。     

地震勘探基于波场分离的逆时偏移去噪方法

基于双程波动方程的逆时偏移（RTM）具有原理简单、无倾角限制、可以适用于速度任意变化的模型,并具有对各类波正确成像等优点,是当今最好的偏移方法之一.本文阐述了逆时偏移存在的缺陷并分析了低频噪声的产生机制,重点探讨了使用波场分离的逆时偏移去噪方法,并通过模型的试算证明了该方法的有效性.