南海深水多次波成像

南海海域水深可达5 km,是我国未来油气勘探最有前景的区域.南海海域地震资料含有各种类型多次波,在地震成像前需将多次波去除.然而与传统的一次反射波偏移比较,多次波携带更多的地下结构信息,多次波偏移可提供更好的地下成像.本文提出了一种修改传统逆时偏移的方法,利用多次波进行成像.该方法用包含一次波和多次波的记录代替震源子波,用预测多次波代替一次反射波作为输入数据.在逆时偏移处理流程中,将地表多次波记录沿时间延拓到各个深度层,通过互相关成像条件将包含一次波和多次波的地震记录沿时间延拓到相同深度层.将多次波偏移结果与一次波偏移结果的差应用目标函数进行匹配,弥补一次波偏移的不足.对Sigsbee2B数据做数值测试的结果表明与传统一次反射波偏移相比,该方法可以得到范围更大、均衡性更好的地下照明度.将该方法应用到南海深水实际资料成像,获得了南海深水地震资料多次波成像剖面.偏移结果表明,多次波成像可以弥补反射波照明的不足,由此说明多次波偏移是未来的一个重要研究方向.     

逆时偏移震源子波影响分析

逆时偏移作为一种先进的地震偏移成像方法,其成像结果的好坏取决于很多因素,其中成像奈件是关键的一个因素.成像条件一般采用互相关成像条件,它是将源波场沿时间进行正推,接收波场沿时间进行逆推,然后将源波场和接收波场进行互相关,从而得到偏移成像结果.源波场正推和接收波场逆推都是采用同一个波动方程,这里采用波动方程一阶应力-速度形式,数值计算方法采用交错网格有限差分方法.采用振幅补偿拉普拉斯滤波方法压制逆时偏移成像中的低频噪声.采用三次样条插值方法解决逆时偏移成像中的波形不光滑问题.在接收波场逆推过程中,是将地震记录作为边界条件.而在源波场正推过程中,需要给定震源子波作为初值条件.理论上来讲,震源子波为脉冲时是最佳的,但在实际当中难以实现.在源波场正推时,一般采用子波函数.子波函数形式非常多,最常用的是雷克子波.不同主频震源子波产生的源波场频率也不同,从而导致逆时偏移互相关成像效果也不同.那么震源子波主频对逆时偏移成像影响到底有多大呢?为了简便,本文基于不同主频雷克子波,初步分析了相应的逆时偏移成像效果,同时对速度模型存在不同误差时的情形也进行了试验分析.结果 表明:在不同速度误差情况下,震源子波主频越低,成像同相轴连续性越好,散射点聚焦能量越强;震源子波主频越高,虽然分辨率提高,但成像同相轴连续性变差,散射点聚焦能量变弱,同时背景噪声增强.本文研究结果对逆时偏移成像技术在实践中的应用具有一定的参考价值.     

地震多次波成像

地震资料含有各种类型多次波,而传统成像方法仅利用地震一次反射波成像,在地震成像前需将多次波去除.然而,多次波携带了丰富的地下结构信息,多次波偏移能够提供除反射波外的额外地下照明.修改传统逆时偏移方法,用包含一次反射波和多次波的原始记录代替震源子波,将SRME方法预测的表面多次波代替一次反射波作为输入数据,可将表面多次波成像.多次波成像的挑战和困难在于大量串扰噪声的产生,针对表面多次波成像中的成像噪声问题,将最小二乘逆时偏移方法与多次波分阶思想结合起来,发展可控阶数的表面多次波反演成像方法,有望初步实现高精度的表面多次波成像.在消除原始记录中的表面多次波后,通过逆散射级数方法预测得到层间多次波,将层间多次波作为逆时偏移方法的输入数据可将其准确归位到地下反射位置.数值实验表明,多次波成像能够有效地为地下提供额外照明,而可控阶表面多次波最小二乘逆时偏移成像方法几乎完全避免成像噪声.     

多次波分阶逆时偏移成像

本文深入分析了多次波逆时偏移的成像原理和串扰假象产生机制,并提出了多次波分阶逆时偏移策略,即首先应用自由界面多次波衰减方法对原始炮集记录中的多次波进行剔除得到一次波记录,然后应用一次波记录预测获得各阶多次波记录,最后将各阶多次波记录分别进行逆时偏移成像.模型实验结果表明,多次波分阶逆时偏移其各阶多次波成像剖面均可对真实界面正确成像且其能够有效压制串扰假象,其中一阶多次波逆时偏移剖面的成像精度最高,其深层成像质量明显优于常规多次波逆时偏移.     

基于解析时间波场外推与波场分解的逆时偏移方法研究

双程波方程逆时深度偏移是复杂介质高精度成像的有效技术,但其结果中通常包含成像方法引起的噪音和假象,一般的滤波方法会破坏成像剖面上的振幅,其中的假象也会给后续地质解释带来困扰.将波场进行方向分解然后实现入射波与反射波的相关成像能够有效地消除这类成像噪音,并提高逆时偏移成像质量.波传播方向的分解通常在频率波数域实现,它会占用大量的存储和计算资源,不便于在沿时间外推的逆时深度偏移中应用.本文提出解析时间波场外推方法,可以在时间外推的每个时间片上实现波传播方向的显式分解,逆时深度偏移中利用分解后的炮检波场进行对应的相关运算,实现成像噪音和成像信号的分离.在模型和实际数据上的测试表明,相比于常规互相关逆时偏移成像结果,本文方法能够有效地消除低频成像噪音和特殊地质构造导致的成像假象.     

基于波场分离的弹性波逆时偏移

尽管叠前逆时偏移成像精度高,但仅针对单一纵波的成像也可能形成地下介质成像盲区,由于基于弹性波方程的逆时偏移成像可形成多波模式的成像数据,因此弹性波逆时偏移成像可提供更为丰富的地下构造信息.本文依据各向同性介质的一阶速度-应力方程组构建震源和检波点矢量波场,再利用Helmholtz分解提取纯纵波和纯横波波场,使用震源归一化的互相关成像条件获得纯波成像,避免了直接使用坐标分量成像而引起的纵横波串扰问题.针对转换波成像的极性反转问题,文中提出一种共炮域极性校正方法.为有效节约存储成本,也提出一种适用于弹性波逆时偏移的震源波场逆时重建方法,在震源波场正传过程中,仅保存PML边界内若干层的速度分量波场,进而逆时重建出所有分量的震源波场.本文分别对地堑模型和Marmousi2模型进行了弹性波逆时偏移成像测试,结果表明:所提出的共炮域极性校正方法正确有效,基于波场分离的弹性波逆时偏移成像的纯波数据能够对复杂地下构造准确成像.     

基于边界存储的地震波场重构及逆时偏移成像应用研究

逆时偏移是基于双程波动理论,没有倾角和横向变速限制,能够解决复杂地质构造成像问题.在逆时偏移成像实现过程中,需要对源波场进行正推和对接收波场进行逆推,由于这两个波场外推方向不一致,需要对其中一个波场(源波场或接收波场)进行存储,这部分存储量很大,是制约逆时偏移成像应用的瓶颈.本文探讨了一种有效的边界存储方法,即先让源波场正推到最大记录时间,在研究区域边界存储所有时刻波场记录,同时记录最大时刻该区域全波场.然后将边界存储波场和最大时刻全波场分别作为边界条件和初始条件,进行波场逆推来重构源波场,这样就解决了源波场与接收波场外推方向不一致的问题,从而可以克服逆时偏移成像存储瓶颈,并将该方法应用到逆时偏移成像.通过简单层状模型和复杂逆掩断层模型试验,并和直接存储方法进行了对比,证明了该方法的有效性.     

基于自适应优化有限差分方法的全波VSP逆时偏移

与地面地震资料相比,VSP资料具有分辨率高、环境噪声小及能更好地反映井旁信息等优点.常规VSP偏移主要对上行反射波进行成像,存在照明度低、成像范围受限等问题.为了增加照明度、拓宽成像范围、提高成像精度,本文采用直达波除外的所有声波波场数据(全波),包括一次反射波、多次反射波等进行叠前逆时偏移成像.针对逆时偏移中的四个关键问题,即波场延拓、吸收边界条件、成像条件及低频噪声的压制,本文分别采用自适应变空间差分算子长度的优化有限差分方法(自适应优化有限差分方法)求解二维声波波动方程以实现高精度、高效率的波场延拓,采用混合吸收边界条件压制因计算区域有限所引起的人工边界反射,采用震源归一化零延迟互相关成像条件进行成像,采用拉普拉斯滤波方法压制逆时偏移中产生的低频噪声.本文对VSP模型数据的逆时偏移成像进行了分析,结果表明:自适应优化有限差分方法比传统有限差分方法具有更高的模拟精度与计算效率,适用于VSP逆时偏移成像;全波场VSP逆时偏移成像比上行波VSP逆时偏移的成像范围大、成像效果好;相对于反褶积成像条件,震源归一化零延迟互相关成像条件具有稳定性好、计算效率高等优点.将本文方法应用于某实际VSP资料的逆时偏移成像,进一步验证了本文方法的正确性和有效性.     

各向异性弹性波动方程多分量联合叠后逆时偏移

为了更好地描述地下介质的各向异性特性,实现准确的波场成像,开展了各向异性介质逆时偏移研究。采用高阶交错网格有限差分法推导出各向异性弹性波动方程叠后逆时深度偏移方程,证明各向异性介质向各向同性介质转化的条件,并研究Thom son参数对弹性波场的影响,合成了各向同性介质和各向异性介质多分量零偏移距剖面,实现了多分量模拟记录联合叠后逆时偏移。计算结果表明,能够准确地实现多分量波场叠后逆时成像,使地质层位中的断层、断点等目标成像清晰准确,且偏移成像精度较高,因此各向异性逆时偏移法是一种高精度的深度域偏移成像方法,可以指导实际天然地震资料和人工地震资料的数据处理。     

地表旋转控制照明叠前深度偏移

基于炮集数据的常规叠前深度偏移波场外推的计算量非常大,控制照明叠前深度偏移是基于快速的组合炮技术,通过在目的层上控制震源波场,高效实现对地质目标的高质量成像. 组合炮所需的合成算子,由在目的层定义的震源波场反向外推到地表建立. 本文提出一种更快速的地表旋转控制照明方法,直接在地表对震源波场进行面向目标的照明控制,实现目标区域高质量成像. 该方法只需对某一核合成算子进行简单的几何旋转建立合成算子函数,避免了复杂的波场外推. 这种新的控制照明方法与快速的一般组合炮偏移效率相同,但对复杂目标区域的成像能力显著增强. 应用本文方法对Marmousi模型进行计算,得到深层复杂构造非常理想的成像结果.     

基于干涉成像条件的探地雷达衰减补偿逆时偏移成像

探地雷达(GPR)衰减补偿逆时偏移方法采用互相关成像条件对实际地下结构成像时，由介质小尺度不均匀所产生的散射波和噪声能量在电磁波衰减补偿逆时外推中会急剧增大，严重影响成像分辨率.为此，本文提出一种基于干涉成像条件的GPR衰减补偿逆时偏移方法.该方法通过在逆时外推过程中，人为改变电磁波动方程中衰减项前的正负号，以保持逆时外推的时间对称性，补偿衰减的电磁波能量；干涉成像条件利用Wigner分布函数(WDF)对所有时刻的逆时外推电磁波场进行滤波，以有效压制逆时偏移剖面中的低频噪声及成像伪影.数值试验结果表明：相比于互相关成像条件，干涉成像条件可有效压制衰减补偿逆时偏移过程中由于背景小尺度不均匀性所产生的散射波和低频噪声，且几乎不影响目标体的成像能量；在计算效率相当的同时，成像分辨率更高.     

基于单程波偏移算子的地表相关多次波成像

在常规地震资料处理中,多次反射波被视为噪声并从地震数据中去除,以免在之后的地震资料解释中造成误解.而事实上,多次波也是地震信号,是照明波场的一部分,能够对地下构造成像的精度做出贡献.本文分析了多次波在传统单程波叠前深度偏移中产生构造假象的机制和表现,为实现基于单程波偏移算子的多次波成像,修改了单程波叠前深度偏移的边界条件,即将输入的震源波场用包含多次波的记录来替代,输入的记录波场用预测出的表层相关多次波来替代,实现了基于单程波偏移算子的地表相关多次波成像,并从理论上给出了其成像依据.通过基于二范式最小能量差原则求取的匹配因子,将多次波成像结果与一次波成像结果进行匹配叠加,应用多次波成像来弥补一次波成像的不足.简单模型验证了基于单程波偏移算子的多次波成像方法的有效性,最后对Sigsbee2B模型进行了一次波与多次波联合成像试算,盐边界高陡构造成像质量得到了明显改善.     

基于波场分离归一化成像条件的被动源数据多次波逆时偏移成像（英文）

随着地震勘探技术的进步,被动源数据越来越受到重视,噪声型被动源数据在自然界和工业生产中广泛存在,利用好被动源数据对随钻测井,大尺度构造勘探等领域有重要意义。本文针对噪声性被动源地震数据等特性,提出了基于波场分离归一化成像条件的被动源数据多次波逆时偏移成像方法。不同于地震干涉法,本方法不需要重构虚拟炮集,对噪声型被动源数据本身偏移成像。将逆时偏移结合多次波成像的边界条件应用于被动源地震数据,并采用波场分离的方法消除逆时偏移的低频噪声问题,同时采用能量归一化的成像条件,克服了被动源震源分布不均匀导致成像能量比例失调问题,同时增强了边缘和深部地层成像质量,且计算效率高。最后用Marmousi模型进行测试,证明了本方法的有效性。     

反射声波测井逆时偏移中的波场存储策略优化方法（英文）

本文从反射声波测井逆时偏移成像方法出发,考虑逆时偏移本身存在海量波场数据存储问题,探讨了波场存储优化策略,提出并实现了激发振幅成像条件逆时偏移和随机边界法逆时偏移两种避免波场数据存储的优化方法。首先,激发振幅成像条件方法利用计算激发时间作为成像时间判断,只需存储小部分波场,如几十个时间点的波场数据,通过实例证明甚至仅需两个时间点波场即可完成成像。传统逆时偏移通常需存储数千时间点的波场数据,与之相比数据存储和读写量可缩小几十倍甚至数千倍;其次,随机边界方法利用波场在随机介质中发生散射而非反射进行震源声场反传重构,从而避免正传波场数据的存储。数值计算表明两种方法均可有效减少反射声波逆时偏移的波场数据存储,与其它偏移算法和常规逆时偏移相比,在保证成像精度的同时提高了数据处理效率,为井旁缝洞体高效、高精度成像提供了重要的技术手段。     

TTI介质叠前逆时偏移成像研究综述

逆时偏移是基于双程波波动方程进行波场逆时外推,再利用成像条件提取成像值,进而实现偏移成像的一种方法.目前,RTM已经由各向同性介质推广到了TI介质,包括VTI介质和TTI介质.考虑了各向异性因素的逆时偏移技术能更好地满足日益增长的勘探精度的要求,进一步提高和改善成像效果.在回顾各向异性介质逆时偏移发展历史的基础之上,本文针对于国内外TTI介质逆时偏移当前的研究现状进行了概述,并对目前TTI介质逆时偏移方法中所存在的特有的数值稳定性和qSV波人为干扰压制问题的解决方法进行了归纳总结.最后对逆时偏移未来的发展方向和研究重点进行了展望.     

表面多次波最小二乘逆时偏移成像

使用相同的炮记录,多次波偏移能提供比反射波偏移更广的地下照明和更多的地下覆盖但是同时产生很多的串声噪声.相比传统逆时偏移,最小二乘逆时偏移反演的反射波成像结果具有更高的分辨率和更均衡的振幅.我们主要利用最小二乘逆时偏移压制多次波偏移产生的串声噪声.多次波最小二乘逆时偏移通常需要一定的迭代次数以较好地消除串声噪声.若提前将一阶多次波从所有阶数的多次波中过滤出来,使用相同的迭代次数,一阶多次波的最小二乘逆时偏移能够得到具有更高信噪比的成像剖面,而且能够提供与多次波最小二乘逆时偏移相似的有效地下结构成像.     

基于变网格的不同观测系统下的逆时偏移（英文）

逆时偏移被认为是处理陡倾角断层、回转波和逆掩推覆构造最有效的成像方法。考虑到当今勘探对地下复杂构造体、薄互层和隐蔽储层的高分辨率成像需求,井间、井中地震资料成为地面地震资料的一个很重要的补充。但由于地面地震、井中地震和井间地震存在很大的频率差异,本文将相应的变网格算法应用到偏移中。另外,也将Lanczos滤波引入到波场外推计算中,从而可以很好地消除变网格引起的人为反射。试算表明,对比地面、井中和井间的偏移,可以看到微构造、陡反射层都得到了更好地刻画和描述。基于变网格的逆时偏移在储层预测和岩性解释方面有广阔的应用前景。     

变阶数有限差分法及逆时偏移有效边界存储最优化策略

地震资料逆时偏移方法在计算过程中需要保存大量的波场信息,这种巨大的存储需求限制了该方法的实际应用.因此,研究能够显著降低波场存储量的存储策略对于推进逆时偏移方法的实用化具有重要的理论意义和实际价值.本文分别从时间和空间的角度提出了相应的波场存储策略:(1)提出变阶数有限差分算子,在计算波场左、右和下边界处时,将差分阶数由高阶递减至二阶,结合有效边界存储策略,可在不损失波场重建精度的同时大幅降低波场存储量;(2)借鉴checkpoint策略和缓冲区思想,提出了基于拉格朗日乘子法的参数最优化方案,在时间方向上压缩波场存储量.通过对二维盐丘模型进行震源波场重建和逆时偏移存储量测试,结果表明:本文方法可实现对震源波场的高精度重建,与传统方法相比,在不损失成像精度的前提下,新的存储策略可大幅降低逆时偏移方法的波场存储需求.     

一种基于逆时偏移的高效照明补偿成像方法

对地震波偏移成像而言,由于观测采集系统不规则、不完备采样以及地下复杂构造对波场传播的影响,在成像点处通常存在照明不均匀的现象,对偏移成像结果进行照明补偿是提高地质体成像精度的有效方法.本文给出了一种基于逆时偏移实现的高效照明补偿成像方法.在高频近似下,本文首先定义了衡量地下照明均匀程度的雅克比矩阵,并以此为基础给出照明补偿算子.通过波场的边界积分表达式,进一步将照明补偿算子利用外推波场进行表达,避免了照明补偿算子中高频格林函数的显式计算,以此构造出高效的照明补偿成像算法.合成数据和实际数据的数值实验结果验证了本文提出方法的正确性.     

多分量联合逆时偏移最佳匹配层吸收边界

有限空间内的波动方程逆时偏移需要利用有效的边界处理技术用以消除人工截断对偏移结果产生的影响。本文以横向各向同性介质弹性波速度-应力方程为基础,依据传统分裂式最佳匹配层（Perfect Matched Layer,PML）吸收边界技术的思想,推导了应用于逆时偏移的完全匹配层波动方程,并给出了其高阶交错网格有限差分格式。针对由边界处向计算区域内传播的＂反射波＂,以及地震记录排列两端地震同相轴突变对计算区域的影响这两方面问题,本文给出了逆时偏移中吸收层的布设方式。模型和实际资料的弹性波叠前多分量联合逆时深度偏移结果表明本文的边界处理方法取得了较好的吸收效果,获得了好的联合偏移成像结果。