自动拾取的成像空间域走时层析速度反演

本文研究了一种成像空间域走时层析速度反演方法,利用波动方程双平方根算子叠前深度偏移提取的角度域共成像点道集(ADCIGs)作为速度分析道集,并基于剩余曲率自动拟合策略获取高精度的走时差。文中不仅给出ADCIGs道集剩余曲率自动拟合拾取的主要思路,而且给出成像空间域层析反演的实现。模型和实际资料试算验证了该方法具有较高的反演精度和计算效率。     

基于改进的频率域拉冬变换衰减多次波方法

在叠前深度偏移提取的角度域共成像道集中可以对多次波进行衰减,对于给定一个准确偏移速度模型,多次波和一次波具有不同的角度域时差,根据这一特征,实现一次波和多次波的分离。本文在提取的角度域共成像道集中,应用改进的频率域拉冬变换,分离一次波和多次波,由于每个角度域共成像道集皆包含复杂的三维地震波场传播效应,所以该方法具有处理三维地震资料和复杂地下构造的能力。相比于波动理论衰减多次波和传统拉冬变换衰减多次波方法,该方法在提高效率的同时,保证了衰减多次波的准确性。     

逆散射级数和抛物线Radon变换联合的层间多次波压制策略

层间多次波给地震成像和解释带来极大挑战,单独使用基于波动方程预测或基于滤波的方法压制层间多次波效果都不理想,前者对数据的品质要求高且运算量大,而后者则很难处理规律性较差的层间多次波.基于此,提出了一种将逆散射级数法与抛物线Radon变换相结合的层间多次波压制新策略.首先在抛物线Radon域采用自适应非线性滤波对CRP数...     

խ��λ����˫ƽ��������ƫ��Ӧ���о�

双平方根（DSR）方程为波动方程偏移提供了一种新的理论框架。基于“沉降观测”概念的DSR方程叠前深度偏移已经成为了一类重要的地震波成像方法。与单程波方程偏移方法相比，其计算效率要高很多，便于输出角度域的共成像道集，有利于开展波动方程偏移速度分析。DSR方程全偏移方法对计算机条件要求太高，且对地震勘探中大量存在的窄方位三维地震数据不太适用。文章结合DSR方程叠前深度偏移的波场传播算子与成像条件，讨论了它在窄方位三维地震数据成像中的应用问题，对比分析了几种稳相近似DSR方程偏移方法的特点与适用条件。数值试例表明，基于DSR方程的窄（共）方位角叠前深度偏移方法对盐丘（或礁体）、古潜山等复杂构造具有很强的成像能力，值得在石油天然气勘探中推广应用。     

混合域三维VSP波动方程深度偏移方法

 常规VSP成像方法利用VSP-CDP转换或基尔霍夫偏移，存在保幅性差及成像精度低等问题。混合域三维VSP波动方程深度成像方法是将VSP反射P波在深度域直接成像，利用分步傅里叶法在混合域进行波场外推，即在频率—波数域实现波场的相移，在频率—空间域实现时移校正。此法的优点在于能保持横向变速地质构造波动方程波场外推精度的同时提高计算效率。理论模型证明此方法可行，实际资料处理效果明显。     

海洋多次波特征和压制方法初探

多次波问题是海洋地震勘探中最突出的问题之一。多次波的存在,影响地震成像的真实性和可靠性,干扰地震资料的解释。如何有效地压制多次波是解决问题的关键。多次波压制技术分为基于有效波和多次波之间差异的滤波方法和基于波动方程的多次波预测减去法两大类。总结了海洋多次波的特征和识别方法以及各种压制多次波的常用方法,用模拟合成地震记录验证并讨论了局部相干滤波方法压制多次波的效果。     

基于奇异值分解的角度域去噪方法

受地层调谐效应和地震数据品质等因素的影响,角度域共成像点道集(angle domain common image gathers,ADCIGs)存在不同程度的波形拉伸和随机噪声干扰。为了提高ADCIGs及其叠加剖面的成像效果,提出了基于奇异值分解的角度域去噪方法。首先对叠前偏移输出的ADCIGs进行奇异值分解,然后对奇异值进行归一化修正,采用累计贡献率的方法确定降噪阶次,从而实现角度域内的信噪分离和噪声压制。在确定降噪阶次时,采用累计贡献率的方法可以直观地判断各奇异值分量对数据的贡献,便于快速选择降噪阶次。理论模型和实际数据的测试处理结果表明,基于奇异值分解的角度域去噪方法适用于具有水平同相轴的ADCIGs,它能有效分离角度域内的随机干扰,并且能压制高角度处的频率畸变,改善大角度数据的品质。对ADCIGs进行基于奇异值分解的角度域去噪,可进一步提高该叠前道集的精度,从而有效改善角度域叠加剖面的信噪比和分辨率,也为基于叠前道集的速度分析和叠前反演提供了更为准确的数据基础。     

基于角度域共成像点道集的层析速度反演方法研究

碳酸盐岩探区的复杂构造成像与速度反演是目前地震数据处理难点问题。碳酸盐岩构造复杂、目标层极深、断点发育,常规处理方法很难得到比较好的成像与速度反演结果。为此提出基于角度域共成像点道集(ADCIGs)、面向碳酸盐岩储层的层析速度分析方法,利用波动方程双平方根算子叠前深度偏移提取的角度域共成像点道集作为速度分析道集,能够更为准确地反映速度和深度的耦合关系,并且减少了假象等的干扰,使得速度分析结果精度更高。模型和实际结果表明,该方法能够较好地解决碳酸盐岩探区的速度反演问题。     

多次波分阶模拟及最小二乘偏移成像

地震正演模拟是地震数据反演及成像中极其重要的部分。基于波动方程的有限差分法是目前应用最广泛的正演模拟方法,该方法虽然精度高,但不能模拟某一特定成分的波场,缺乏灵活性。基于频率空间域单程传播算符的波场模拟方法（波场分阶模拟算符）不同于有限差分法,该方法可以分别模拟一次波与多次波（包含表面多次波和层间多次波）。在常规地震处理中,多次波往往被视为干扰信息,但多次波中包含着丰富且重要的地下结构信息,利用好这些信息将极大地提高地下结构的成像质量。通过分析多次波在常规偏移成像中的影响,利用波场分阶模拟得到一阶表面多次波以及一次波,并进行表面多次波最小二乘偏移成像研究,提高一次波的成像效果,压制表面多次波成像中的串扰,提高地下结构的成像质量。     

用多次波进行深度成像

多次波深度成像是一种叠前深度偏移方法,它是将多次波作为更准确的边界成图和振幅恢复的有效信号。这一思路与基于模型的多次波压制技术和逆时深度偏移有着部分的联系。常规的逆时偏移是使用计算所记录地震道和射线追踪的反射波传播的双程波动方程,或者是使用计算直达波旅行时的程函方程(激励时间成像原理),因而逆时偏移与大多数其他的只需一次反射波的单程波动方程方法或射线追踪偏移方法大致相同。因为如果不减弱一次波就几乎不可能衰减多次波,所以需要寻求一种在资料处理中建设性地使用多次波而不是想着去破坏它们的方法。此外,如果把一个双程全波动方程用于偏移中,那么多次波和其他类型的非反射波就能加强边界成像和振幅恢复。新方法通过有限差分技术解出了波在正向和反向方向上传播的双程波动方程。因此,它可考虑所有类型的声波,例如反射(一次波和多次波)、折射、绕射、透射波和这些波的任意组合。 在成像过程中,所有这些不同类型波场的衰减都发生在生成或改变它们的边界上。这一成像过程经过了四个主要步骤。首先,利用双程标量全波动方程使源函数(子波)从震源位置向所有方向正向推进;其次,用双程标量全波动方程使炮点道集中的记录道从全部检波器点向所有方向反向进展;然后,对正向和反向传播波场进     

聚束滤波法消除海上地震资料的多次波

本文简单介绍了用聚束滤波方法消除地震资料中的多次波的基本原理及合成资料处理结果；详细描述了用聚束滤波方法消除两套南海深海地区地震资料中多次波的实施步骤与效果。对于珊瑚礁产生的大范围强烈多次波，采用多次使用聚束滤波方法分别消除不同范围的多次波，取得了较高分辨率和较小计算量的优化效果。处理结果表明，聚束滤波方法可以消除相对较强的多次波，并增强较弱的一次波。对莺东南盆地一块海底多次波相当严重的深水地震资料的处理结果表明，聚束滤波方法对于消除远道和近道多次波的效果都相当好，并且优于常规的二维滤波方法对多次波的压制效果。     

基于3D匹配滤波器和伪地震数据算法的多次波自适应相减方法

多次波自适应相减是预测减去法压制多次波的关键步骤。为进一步去除残余多次波,基于常规2D匹配滤波方法,文中引入3D匹配滤波器,同时利用多个预测多次波道集以匹配原始数据。针对3D匹配滤波器可能造成的一次波损伤现象,利用相同的3D匹配滤波器同时拟合多个原始数据道集;同时,引入伪地震数据算法求解对一次波施加Huber范数最小化约束的优化问题,不需满足一次波与多次波正交的假设,能有效分离一次波与多次波。另外,在整个迭代过程中,伪地震数据算法只需利用Cholesky分解算法进行一次矩阵分解,计算效率较高。模型和实际数据的处理结果表明,与基于一次波能量最小化的3D匹配滤波器方法和基于伪地震数据算法的2D匹配滤波器方法相比,所提方法能更好地均衡一次波保护与多次波分离。     

基于-次波逆时理论的层间多次波衰减方法

针对层间多次波与-次波的动校正速度差异较小,应用基于动校速度差异的常规方法无法将层间多次波彻底去除的难题,本文对基于波动理论的层间多次波预测相减方法进行了研究。通过对地震数据进行-定的前期处理,选定某-地层,并沿该层面进行褶积运算,预测与其相关的所有层间多次波;然后应用传统的最小二乘自适应相减算法,压制层间多次波。模拟数据的测试结果验证了方法的有效性。     

基于共散射点道集的多次波压制方法研究

多次波压制是地震数据处理中的重要环节。一种常用的多次波压制方法是根据一次波和多次波在成像空间的曲率差异并借助Radon变换进行分离和压制。成像空间既可以是叠前时间域也可以是叠前深度域。考虑到对复杂介质的适应性以及叠前处理对计算效率的要求,提出基于共散射点(CSP)道集的高分辨率Radon变换多次波压制方法。该方法只需要一个简单的初始速度场,就可将常规共中心点道集(CMP)映射到覆盖次数更高、炮检距覆盖范围更广的共散射点道集,然后在该道集上应用高分辨率双曲Radon变换,可较好地分离一次波和多次波。相对于以水平层状介质为假设条件的CMP道集,CSP道集更适应复杂地质构造,且时距关系满足双曲规律。模型和实际资料测试结果表明,该方法可以较好地实现速度谱能量团的聚焦,有利于较复杂地质条件下的多次波压制。     

波动方程叠前深度偏移角度道集求取方法

针对常规直接转换法和直射线叠前时间偏移(PSTM)法提取的角度道集无法满足地质构造复杂地区的AVA分析与反演精度的问题,研究基于共炮检距道集波动方程保幅叠前深度偏移方法,提取角度域共成像点道集。首先,对地震数据进行叠前深度偏移得到炮检距域共成像点道集;其次,在频率-波数域对炮检距域共成像点道集利用快速插值映射法抽取角度域共成像点道集。通过与直接转换法和直射线PSTM法提取的角度道集比较,分析不同方法求得的角度道集信息的准确性和AVA特征。数值试算和实际资料处理结果表明:共炮检距道集波动方程保幅叠前深度偏移方法求取的角度道集的振幅相对保真,角度范围更广,更适用于AVA分析。     

数据增广的编解码卷积网络地震层间多次波压制方法

层间多次波压制是地震资料去噪领域的一项前沿技术挑战，对获取高质量数据、了解地下真实构造具有十分重要的意义。现有的层间多次波压制方法耗时长，对人工参数调整要求高，处理低信噪比数据时可能导致层间多次波泄露。为此，提出一种基于数据增广的编解码卷积神经网络层间多次波压制方法。首先，利用基于虚同相轴的层间多次波压制方法从原始数据中估计出一次波和层间多次波，生成一次波标签数据。然后，构建两种增广训练集：一方面，通过改变训练样本中层间多次波的振幅、极性及旅行时，进行层间多次波波场数据的增广，提高层间多次波压制网络的泛化能力；另一方面，通过对原始数据添加不同信噪比的高斯噪声进行噪声注入的数据增广，提高网络的抗噪性能。最后，结合去噪卷积神经网络(DnCNN)和U形全卷积神经网络(U-Net)的优势搭建了适合层间多次波压制的深层编、解码网络，进行神经网络训练和预测。合成数据和实际数据的处理结果表明，该方法能够有效压制地震层间多次波并保护一次波，具有较强的泛化能力和抗噪性能，可显著提高计算效率。     

琼东南深海盆地三维地震多次波组合压制方法

在琼东南深海盆地,浅层为海相沉积特征,中、深层为陆相沉积特征,原始地震资料中多次波不仅异常发育而且复杂.压制多次波的方法有很多(褶积类、速度滤波类、波动方程类),但都具有一定的局限性,不能适应复杂情况下的多次波压制.为此,对各种压制多次波的方法进行了分析,针对琼东南深海盆地三维地震资料多次波的特点,采用调频滤波、抛物线Radon变换、加权中值滤波等方法,在不同的处理阶段分别对海底多次波、层间多次波和畸变多次波进行了压制,获得了好的效果.应用这种组合方式,不仅多次波被完全压制,一次有效反射得到了清晰的显示,而且没有损害有效信号.     

弹性波逆时偏移的角度域共成像点道集提取方法

对于多分量地震勘探,由于介质中纵、横波的耦合作用,多波角道集存在能量串扰和低波数噪声干扰等问题。针对横波Poynting矢量构建过程中存在的应力串扰问题,基于解耦延拓方程及伪横波应力构建方法,对质点振动速度和应力进行解耦,实现了无能量串扰的纵波和横波Poynting矢量构建;在基于Poynting矢量提取PP波反射角和PS波反射角的基础上,提出了基于弹性逆时偏移的PP波和PS波角道集提取方法;针对PP波角度域共成像点道集存在的低波数噪声干扰问题,提出了角度域共成像点道集的角度衰减叠加策略,提高了偏移成像精度。应用模型数据验证了方法的正确性。     

迭代抛物Radon变换法分离一次波与多次波

Radon变换法是进行一次波与多次波分离的常用手段,最小平方约束下的频率域抛物Radon变换将t-x域数据转换到Radon域后,因存在剪刀状发散的截断效应,用传统方法难以彻底分离一次波和多次波。针对这一缺陷,提出了迭代抛物Radon变换法,即在Radon域截取一次波聚焦点附近很小区域内的数据为初始数据,经过Radon反变换和正变换后得到新的Radon域数据,然后用初始数据覆盖对应的小区域,经过迭代,最终得到保幅效果较好的一次波,且几乎不含多次波。利用相同的方法,也可以得到几乎不含一次波的多次波。最后通过理论模型和实际资料的处理,验证了本文方法的正确性和有效性。     

基于双平方根方程的双域保幅地震偏移方法

从全声波方程出发,基于严格的解耦理论进行单程波保幅分解,得到保幅的单程波动方程;然后按照叠前深度偏移中的波场双重向下延拓理论实现共炮域公式到CMP域公式的转换,推导出由双平方根方程定义的保幅单程波动方程;基于反问题求解中常用的摄动理论,把速度场分裂为层内常速背景和变速扰动,分别在频率波数域和频率空间域求得波场深度延拓的偏移时移量及振幅校正系数,得到最终的双域保幅波场延拓算子。理论模型试算和实际资料处理表明,该方法不但具有较高成像精度,而且可以为后续的地震属性分析提供更真实的振幅信息。