一种稳定的波场延拓反Q滤波算法

反Q滤波算法中稳定性是关键。针对地下介质为层状Q模型结构,本文提出一种稳定的反Q滤波算法,基于波场延拓原理,对于每个常Q层反Q滤波算法分两步完成:(1)利用一种稳定的波场延拓算法,将地表波场记录直接延拓到当前层上;(2)在当前层内进行常Q反滤波。在步骤(1)中,利用稳定的波场延拓算法并结合稳定因子,将地表波场记录直接延拓到当前层上,避免了上覆层反Q滤波误差的累积。在步骤(2)中,对于当前层内的反Q滤波,通过对该层地震信号的Gabor谱分析,拾取时变的限幅频率,并确定相应于该频率值的限幅增益值来实现反Q滤波算法的稳定性。最后通过理论模型与实际数据验证了算法的优越性。     

两种抗噪反Q滤波方法对比及应用

本文就两种抗噪反Q滤波方法进行了讨论,其中对比了考虑时频域信噪比的反Q滤波方法和基于变稳定因子的反Q滤波方法.与稳定的反Q滤波方法相比,这两种方法可以抑制噪声,提高反Q滤波后地震数据的分辨率和信噪比.然而,它们基于不同的方法原理,并且在处理效果上存在一定的差异.为了探究反Q滤波方法抑制噪声的关键问题,回顾了方法的原理,对比和讨论这两种方法的差异,并通过理论模型数据和实际地震数据进行了比较.测试结果表明,反Q滤波方法抑制噪声的关键在于有效控制振幅补偿的频带范围.就算法而言,考虑时频域信噪比的反Q滤波方法可根据地震数据的信噪比水平选择不同的频带补偿范围,具有较强的灵活性,但需要先计算出地震数据的时频域信噪比;而基于变稳定因子的反Q滤波方法具有更加简单的算法,其振幅补偿规律本身即具有一定的噪声压制能力,避免了计算时频域信噪比.当选取合适的参数,二者的振幅补偿频带范围接近时,这两种抗噪反Q滤波方法可以得到相似的处理效果.     

基于等效加权模型的单步延拓反Q滤波方法研究及应用

通过建立等效加权模型并进行单步向上延拓反Q滤波,减弱目前常规Q值模型对反Q滤波效果的限制,达到进一步提高地震记录分辨率的目的.目前,常规Q值模型包括常Q模型和层状模型,常Q模型与复杂的地下情况相差较大,导致地震波能量衰减补偿效果较差;而基于层状模型进行反Q滤波则会导致噪声逐层放大,降低地震记录信噪比.基于等效加权模型进行单步向上延拓反Q滤波,可以有效补偿地震波振幅衰减并校正相位畸变,在保持信噪比的同时提高地震记录分辨率.通过模型测试,基于等效加权模型进行单步向上延拓反Q滤波,可以在有效压制噪声的同时增强能量衰减补偿效果;实际资料应用的结果表明:该方法不仅能有效补偿深层能量衰减,并且能较好的压制环境噪声,在提高地震资料分辨率的同时保持了信噪比.     

一种改进的广义稳定反Q滤波方法

由于地层的黏弹性和非均质性,地震波在传播过程中振幅衰减,相位畸变,导致地震资料分辨率和信噪比降低.反Q滤波是补偿地层衰减的一种常用方法.常规的反Q滤波方法的振幅补偿算子只与时间、频率有关,其振幅补偿频带范围、强度固定不变,不能够根据地震记录的自身特点灵活地进行振幅补偿.基于此,本文在广义稳定反Q滤波法的基础上,提出了一种改进的广义稳定反Q滤波方法.新方法提出了一种新的稳定因子,并通过参数p对振幅补偿算子进行了优化.其振幅补偿算子与地震资料信噪比、走时、频率、品质因子相适应,可以更加精确地对振幅进行补偿.理论数据测试和实际地震资料应用均表明,新方法在振幅补偿和噪声压制上优于广义稳定反Q滤波法,其处理后的地震资料振幅恢复保真度、分辨率和信噪比更高,对储层预测和精细解释工作具有重要意义.     

基于反演的衰减补偿方法(英文)

提高地震资料分辨率的一个有效途径就是衰减补偿,通过对地震波的衰减和频散效应进行校正,提高地震资料的分辨率。常规衰减补偿方法都是基于波场延拓的反Q滤波方法。本文利用Futterman衰减模型,导出了一种衰减介质中合成地震记录的计算方法,在此基础上将衰减补偿问题归结为一个Fredholm积分方程反问题,利用反演方法来实现衰减补偿。针对衰减补偿问题的不稳定性,利用Tikhonov正则化方法提高反演过程的稳定性,数值模拟资料和实际资料处理结果验证了方法的有效性。     

一种自适应增益限的反Q滤波

地层的Q吸收会造成地震波振幅衰减、相位畸变,分辨率和信噪比明显降低.反Q滤波可消除由于地层Q吸收造成的振幅衰减和相位畸变,从而提高地震资料的分辨率;但反Q滤波振幅补偿的数值不稳定性问题会严重降低地震资料的信噪比,并产生很多假象.截止频率法和稳定因子法反Q滤波振幅补偿方法虽可控制数值非稳定性问题,但振幅补偿函数的增益限为一个时不变的常数,且与地震数据动态范围无关,其经常会压制深层地震波的高频成分,反而降低地震资料的分辨率;因此,本文在研究截止频率法和稳定因子法的基础上,结合地震数据的动态范围对地震记录分辨率的影响,提出了一种自适应增益限的反Q滤波振幅补偿方法,其增益限和稳定因子都是时变的,且都自适应于地震数据有效频带的截止频率.合成数据和实际数据试算表明,本文的自适应增益限的反Q滤波方法可恢复地震信号有效频带范围内的能量,且能较好地控制数值非稳定性问题,最终获得高分辨率和高信噪比的地震数据.     

用稳定高效的反Q滤波技术提高地震资料分辨率

地震波在地下传播时受到衰减影响,衰减会导致地震波场高频能量的损失和相位畸变.反Q滤波可补偿大地衰减效应.已有的反Q滤波方法存在下列不足:频率域的算法由于算子长度较长,所以计算效率较低;时间域的算法,或者对地震记录上到时较晚的同相轴进行了过度的补偿,或者为防止过度补偿后来的振幅而在最大增益处进行限制,导致振幅的多解性,而且还会影响滤波器的相位效应.本文给出一种通过直接求解时间域的Q模型方程来进行反Q滤波的算法.由于采用带状矩阵解算器,所以具有较高的计算效率,理论数据和实际地震资料的试算结果证明,本方法对地震波的吸收衰减进行了出色的补偿,提高了地震资料的分辨率.     

地震波频散效应与反Q滤波相位补偿

本文以Futterman提出的地震波振幅衰减和相速度频散表达式为基础,从井震匹配的角度出发,推导了不进行反Q滤波、仅反Q滤波振幅补偿、仅反Q滤波相位补偿,以及反Q滤波相位与振幅同时补偿四种情况下,地震波速度频散与相速度及地震记录振幅谱间的表达式,并从理论上说明了反Q滤波相位补偿的必要性.通过相同观测系统,炸药震源和可控震源分别激发采集的零偏移距VSP资料实例,验证了本文所推导的速度频散表达式的合理性;通过井震标定实例,进一步说明了反Q滤波相位补偿,可有效消除地震子波速度频散,提升地面地震资料与零偏移距VSP走廊叠加剖面的匹配度,最终提高地震资料成果的可信度.     

一种起伏地表条件下的照明补偿方法

起伏的地表条件限制了采集孔径范围并造成深层地震照明不足,为改善该类地区的成像质量,本文提出了一种起伏地表条件下的照明补偿方法.首先,基于小波束波场延拓算子和逐步累加的外推方法在波场延拓过程中解决起伏地表面的影响,并引入空间滤波函数压制虚拟层内的偏移噪音;其次,利用局部指数标架对上、下行波场分解,得到局部角度域成像和照明补偿因子.再次,利用计算出的成像值和照明补偿因子,在局部倾角域完成照明补偿.SEG起伏地表模型测试证明了本方法的有效性,深层构造照明度明显加强,不同角度成像振幅更加均衡,该技术为提高起伏地表地区的成像品质提供了新的手段.     

基于褶积原理和改进广义S变换的震动波衰减补偿试验研究

为了恢复震动波能量在传播过程中产生的衰减损耗,提出基于褶积原理求取品质因子Q的方法与改进广义S变换相结合的反Q滤波法。通过震动波衰减补偿模型试验,对试验数据进行改进广义S变换的时频特性分析,得出了信号的能量分布情况以及时间频率对应关系;采用基于褶积原理求取品质因子的方法,得到时变Q值;对试验数据进行反Q滤波处理,使震动波能量得到了补偿。结果表明本文提出的反Q滤波法提高了对震动波能量衰减补偿的效果,拓宽了地震资料的频带,提高了地震资料分辨率,有利于进行高分辨率地震勘探、深部信号增强和油气藏预测工作的开展。     

叠前纵波和转换波地震资料Q值提取及反Q滤波(英文)

纵波和转换波联合的多波地震勘探技术是解决复杂油气勘探的有效技术,提高转换波的分辨率是其关键问题之一。影响转换波分辨率的主要原因是地层对地震波的吸收,有效地计算地层Q值、消除地层吸收对转换波传播的影响,是提高转换波分辨率的关键。本文提出了从叠前转换波道集中估算横波Q值的方法,并利用沿射线路径的波场延拓,将一种稳定有效的反Q滤波方法应用到叠前共炮点纵波和转换波道集的衰减补偿中。模型资料结果表明,本文提出的估算转换横波Q值的方法精度较高;模型资料和实际资料的吸收补偿结果表明,此稳定全反Q滤波能有效地提高叠前纵波和转换波资料的分辨率。     

基于Bregman迭代的复杂地震波场稀疏域插值方法

在地震勘探中,野外施工条件等因素使观测系统很难记录到完整的地震波场,因此,资料处理中的地震数据插值是一个重要的问题。尤其在复杂构造条件下,缺失的叠前地震数据给后续高精度处理带来严重的影响。压缩感知理论源于解决图像采集问题,主要包含信号的稀疏表征以及数学组合优化问题的求解,它为地震数据插值问题的求解提供了有效的解决方案。在应用压缩感知求解复杂地震波场的插值问题中,如何最佳化表征复杂地震波场以及快速准确的迭代算法是该理论应用的关键问题。Seislet变换是一个特殊针对地震波场表征的稀疏多尺度变换,该方法能有效地压缩地震波同相轴。同时,Bregman迭代算法在以稀疏表征为核心的压缩感知理论中,是一种有效的求解算法,通过选取适当的阈值参数,能够开发地震波动力学预测理论、图像处理变换方法和压缩感知反演算法相结合的地震数据插值方法。本文将地震数据插值问题纳入约束最优化问题,选取能够有效压缩复杂地震波场的OC-seislet稀疏变换,应用Bregman迭代方法求解压缩感知理论框架下的混合范数反问题,提出了Bregman迭代方法中固定阈值选取的H曲线方法,实现地震波场的快速、准确重建。理论模型和实际数据的处理结果验证了基于H曲线准则的Bregman迭代稀疏域插值方法可以有效地恢复复杂波场的缺失信息。     

波场延拓深度滤波方法

消除面波是地震数据处理中的一个重要内容．本文提出了基于15°波动方程的深度滤波方法．由于面波与有效反射波具有不同的传播深度,可利用波场向下延拓方法将二者进行波场分离．把向下延拓后的波场中集中在地表附近的面波能量切除后,再将波场重新延拓回原始的观测面,达到去除干扰的目的．实际资料处理显示：方法计算稳定,消除面波能力强,能更好地保持波场的有效成份和幅值,符合波场的实际传播状态,表明该方法正确可行．     

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

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

基于改进线性Radon变换的井孔地震上下行波场分离方法

反射波场分离是井孔地震资料处理中极其重要的一个环节,波场分离的质量直接影响成像结果的精度.不管是VSP还是井间地震资料,其反射波时距曲线都近似直线型,根据这一特征,本文提出一种改进的线性Radon变换方法来进行井孔资料的反射波上下行波场分离.该方法基于频率域线性Radon变换,通过引入一个新的变量λ来消除变换算子对频率的依赖性,避免了求取每一频率分量对应的不同变换算子,显著降低了计算成本;文中在求解该方法对应的最小二乘问题时,引入了发展较为成熟的高分辨率Radon变换技术来进一步提高波场分离的精度.采用本文方法进行井孔地震资料的上下行波场分离可以在保证分离精度的前提下有效地提高计算效率.根据上下行波在λ-f域内分布的特殊性,设计简单的滤波算子就可实现上下行波场的分离.最后通过合成数据试算以及实际资料处理（VSP数据和井间地震数据）验证了该方法的可行性和有效性.     

基于波动方程的上下缆地震数据鬼波压制方法研究

本文发展基于波动方程的上下缆鬼波压制方法,推导了上下缆地震波场频率波数域波动方程延拓合并公式.基于Fourier变换的波场解析延拓确保上下缆资料振幅相位的一致性,消除了长拖缆远偏移距信号的计算误差,同时具有较高的计算效率;上下缆地震波场的波动方程法合并有效解偶鬼波干涉,实现综合利用上下缆地震数据压制鬼波.理论模型数据和实际采集地震数据的测试表明了方法的有效性.     

基于两步SVD变换的零偏VSP资料上下行波场分离方法

垂直地震剖面(Vertical Seismic Profiling,VSP)资料处理中波场分离是关键问题之一.随着属性提取技术的发展,新的属性参数(例如Q值)提取技术对波场分离的保真性要求越来越高.本文改进了传统奇异值分解(Singular Value Decomposition,SVD)法,给出了一种对波场的动力学特征具有更好的保真性,可以作为Q值提取的预处理步骤的零偏VSP资料上下行波场分离方法.该方法通过两步奇异值分解变换实现:第一步,排齐下行波同相轴,利用SVD变换压制部分下行波能量;第二步,在剩余波场中排齐上行波同相轴,使用SVD变换提取上行波场.在该方法的实现过程中,压制部分下行波能量后的剩余波场中仍然存在较强的下行波干扰,使得上行波同相轴的排齐比较困难.本文给出了一种通过极大化多道数据线性相关程度(Maximize Coherence,MC)排齐同相轴的算法,在一定程度上解决了低信噪比下排齐同相轴的问题.将本文提出的方法用于合成数据和实际资料的处理,并与传统SVD法的处理结果进行对比,结果表明本文提出的波场分离方法具有良好的保真性,得到波场的质量明显优于传统SVD法.通过对本文方法和传统SVD法处理合成数据得到的下行波场提取Q值,然后进行对比可知,本文方法可以有效提高所提取Q值的准确性,适合作为Q值提取的预处理步骤.