半径-斜率域加权反假频地震数据重建

为减小地震数据缺失给地震后续处理工作带来的影响,需要对地震数据进行插值重建.针对反假频插值重建这个难点问题,进行了相关研究,并由此提出了一种改进的R-P（半径-斜率）域加权反假频地震数据插值重建方法.该方法将F-K（频率-波数）谱变换到R-P域,在R-P域设计一个权函数并将其作用于每次的迭代插值过程.通过模型数据和实际数据的测试,证明了该方法具有较好的反假频插值重建能力.     

深度域叠前弹性波阻抗反演（英文）

叠前深度偏移的广泛应用产生了大量的深度域地震数据,提供了深度域地震处理的数据基础,并为探索直接基于深度域地震数据的储层弹性参数反演提供了可能性。由于深度域地震子波存在与速度密切相关的特征,直接提取深度域地震子波存在困难,本文将深度域速度变换与脊回归方法结合,实现了深度域地震子波的确定性提取,将提取的深度域常速地震子波构建成依赖速度的空变地震子波,基于加权叠加思想,提出了深度域地震记录合成方法。再结合非线性共轭梯度算法,建立了基于深度域地震数据和模型驱动的深度域叠前波阻抗反演流程,不再进行额外的时间域和深度域之间的数据转换。深度域合成地震数据与实际地震资料的叠前波阻抗反演取得了较好的结果,表明该方法是可行的。该方法可作为深度域正演分析和叠前弹性参数反演的基础方法,为深度域反演方法的稳定性、精度和效率的进一步研究创造了条件。     

地震数据反假频规则化方法研究

由于诸多因素的影响,地震数据沿空间方向通常是稀疏采样的,因此引起较为严重的空间假频.本文提出一种反假频地震数据规则化的方法,采用Fourier变换域加权范数带限重建方法完成低频数据重建,利用自适应频谱加权范数的正则化项约束方程的解,将地震数据的带宽和谱形状作为先验信息,具有较好的低频重建特性.文中采用共轭梯度算法求解方程,而后利用重建的低频数据信息,应用频带延拓的方法重建高频数据,未知的高频带信息由重建的低频带信息构建.本方法在完成地震数据规则化的同时,可有效去除地震数据中的空间假频干扰.理论模型和实际资料处理均表明文中所提出的反假频地震数据规则化方法是有效可行的.     

基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法

针对利用地震道进行相对波阻抗反演中遇到的横向连续性难以保持、初始子波容错度差以及随机噪声干扰影响反演结果等问题,提出了一种基于矩阵Toeplitz稀疏分解的相对波阻抗反演方法.该方法将地震数据剖面的Toeplitz稀疏分解问题分解为两个子反演问题,其一以Toeplitz子波矩阵元素为待反演的参数,用Fused Lasso方法求解,可保证子波具有紧支集且是光滑的;其二以稀疏反射系数矩阵元素为待反演参数,用基于回溯的快速萎缩阈值迭代算法求解,大大降低了目标函数中参数选择的难度.通过交替迭代求解上述两个子反演问题可将地震数据剖面因式分解为一个Toeplitz子波矩阵和一个稀疏反射系数矩阵;然后由反射系数矩阵递推反演可以得到高分辨率的相对波阻抗剖面;利用测井资料加入低频分量后,也可得到高分辨率的绝对波阻抗剖面.Marmousi2模型生成的合成记录算例和实际地震资料算例均表明:本文方法可以从带限地震数据中有效地反演相对波阻抗,反演结果分辨率高并且能够很好地保持地震数据的横向连续性;即使在初始估计子波存在误差和地震数据被随机噪声污染的情况下也能取得较好的效果.     

地震低频信息缺失特征分析及拓频方法研究（英文）

地震低频信息能够提高分辨率与成像精度,改善反演质量,甚至直接进行油气检测,需要对其进行有效保护与拓展。对于子波而言,缺失低频信息会导致主瓣幅度降低、第一旁瓣幅度增加,并出现次级旁瓣呈周期震荡衰减的现象;从合成地震记录和典型地质模型来看,低频缺失会产生假同相轴,造成分辨率提高的假象,且模型不同位置的特征存在一定差异;对缺失低频的模型数据进行波阻抗反演,会造成构造失真、岩性改变的假象,特别是高陡构造和薄互层。针对缺失低频的地震资料,本文还研究了基于压缩感知与稀疏约束的拓频方法,开发了相应的模块,并对实际CIP道集进行处理,取得了较好的应用效果。     

稀疏反演用于消除气泡和增强频谱

由无气泡气枪震源形成的简单波形可以大大简化子波相位函数的确定和控制,可以提高地震数据的地层学分层的可靠性。本文中,我们提出一种创新方法,即通过利用频率-波数域稀疏反演法用于消除气泡和增强频谱。与著名的气枪震源不同,我们使用目标源的概念。目标源是单一波瓣且无气泡的气枪震源。为了反演使用类似目标源获得的地震数据,我们计算一个估计问题。该方法的基本思想是利用卷积和反卷积,由于稳定性因素在时空域存在随机噪声。我们提出当通过频率-波数(f-k)域值约束做反卷积时,利用迭代消除随机噪声。与传统的维纳滤波法相比,此方法可以获取更为接近完美的结果,并消除其他噪声和人工干扰。我们利用一个线性事件合成数据和一个更为真实的Marmousi模型实例来演示此方法的性能。结果表明,此方法能有效地消除气泡影响并填补频谱凹槽。     

基于LO范数最小化的地球物理数据稀疏重构(英文)

在实际的地球物理数据采集工作中,会因为多方面的客观原因导致数据缺失,对缺失数据进行插值重构是地球物理数据处理和解释的基础问题。基于地球物理数据自身或在变换域内的稀疏性,将地球物理数据的重构转化为稀疏优化问题可提高数据重构的精确度与稳定性。本文建立了LO范数最小化的地球物理数据稀疏重构模型,针对不同规模、不同特征的地球物理数据引入了两种不同类型的LO范数最小优化问题的近似求解算法,即基于LO范数最小化的迭代再加权最小二乘算法与具有快速收敛性的快速迭代硬阈值法。理论分析与数值试验表明,将迭代再加权最小二乘算法应用到位场数据重构中可发挥其收敛速度快,计算时间短,精度高的优势,而快速迭代硬阈值法更适合处理地震数据,相对于传统的迭代硬阈值法计算效率有了很大的提高。     

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

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

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

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

宽频数据在致密砂岩储层预测中的应用

西湖凹陷深层致密砂岩储层具有良好的勘探开发前景,受埋深影响,目的层地震资料品质较差.A构造通过斜缆宽频采集和处理获取宽频地震数据,提升了资料品质,然而应用常规子波提取方法对宽频数据进行子波提取并反演计算纵波阻抗,结果与井上实测数值差异较大,影响储层的定量解释.针对这一问题,提出统计性子波和确定性子波相结合的长短子波合并宽频子波提取方法,提取的宽频子波比常规子波低频丰富、旁瓣小,能更真实地反映地震信息,约束稀疏脉冲反演的纵波阻抗结果与测井曲线吻合度更高.基于宽频数据和常规数据分别进行约束稀疏脉冲弹性波阻抗反演,预测A构造优质储层分布,经已钻井证实,宽频数据比常规数据储层预测精度高,预测的储层展布特征与研究区地质沉积认识一致.结果表明:这种基于宽频子波提取的宽频资料应用方法有效降低了致密砂岩储层预测的多解性.