一种地震资料自适应高分辨处理方法技术

本发明专利技术涉及一种地震资料自适应高分辨处理方法；资料采集处理得到叠后或者偏移后道集；选取地震资料有效带宽范围；在有效带宽范围内除以子波的频谱得到反射系数的部分谱序列；设计自适应滤波器，对部分谱序列进行滤波；估计对应所选时间位置的反射幅度估计；估计最终得到全部的反射系数序列；结合测井资料，井旁道处理，对剖面的应用效果进行评价，得到对整个地震数据的高保真和高分辨处理结果；在得到的更宽频带的道集的基础上，进行进一步的构造，岩性，微裂缝发育解释和评价，寻找页岩气勘探开发有利区域；结合已有探井，地质各方面资料，提供新的井位设计依据；方法可以进行更精细的页岩气储层解释与评价。

【技术实现步骤摘要】
【专利摘要】本专利技术涉及；资料采集处理得到叠后或者偏移后道集；选取地震资料有效带宽范围；在有效带宽范围内除以子波的频谱得到反射系数的部分谱序列；设计自适应滤波器，对部分谱序列进行滤波；估计对应所选时间位置的反射幅度估计；估计最终得到全部的反射系数序列；结合测井资料，井旁道处理，对剖面的应用效果进行评价，得到对整个地震数据的高保真和高分辨处理结果；在得到的更宽频带的道集的基础上，进行进一步的构造，岩性，微裂缝发育解释和评价，寻找页岩气勘探开发有利区域；结合已有探井，地质各方面资料，提供新的井位设计依据；方法可以进行更精细的页岩气储层解释与评价。【专利说明】
本专利技术属于地震资料数据处理方法，是一种基于非参数谱估计理论，通过模拟相对干扰的统计特征，设计自适应滤波器组对其进行压制，从而对不同时间位置的反射波幅度进行稳定准确地估计，从而提高地震数据剖面分辨率，高保真地拓宽地震数据频带的高分辨率处理技术。
技术介绍
我国海相页岩探区储层具有多层系、强改造、高成熟度等明显特征，主要分布在南方山地地区，工区表层和地质条件复杂，相对于常规油气，页岩气勘探开发难度更大。页岩储层分布相对复杂，泥页岩地层虽然很厚，可以达到几百上千米，但是达到有效的页岩气储层标准的只有几十甚至几米。有机碳含量的变化会引起地层密度、速度的变化，但是幅度非常小。这样一来，对页岩气储层的探测和识别时对地震资料的质量，特别是分辨率和保真度方面提出了更高的要求。针对页岩气储层的特征，需要开展具有高保真度的提高分辨率的地震数据处理技术研究。 常规的提高地震数据分辨率方法主要考虑对薄层的识别，以反褶积为代表(黄绪德，1992)，大多基于地震记录的单道的褶积模型(Robinson和Treitel，1980)，在频率域则表现为地震道频谱为子波与反射系数频谱的乘积(Castagna，2004)，主要考虑通过反滤波提高地震记录的视觉分辨率(Porsani和Ursin, 1998)。由于噪声的影响仅能恢复有效通带范围内反射系数的部分频谱，地震数据分辨率仍然由频带范围所决定，在分辨能力和信噪比之间存在矛盾(Levy和Fullagar, 1981)。随着反演理论研究的进展,稀疏反裙积(朱振宇和刘洪，2005)和反演方法(张繁昌，刘杰，印兴耀等，2008)也成为了识别薄层厚度以及顶底界面位置的重要手段。稀疏性的约束反演技术在寻找强反射，以及抗噪声方面表现出了较大的优势(刘喜武，宁俊瑞，张改兰，2009)。 提高地震数据分辨率方法是提取地质体反射特征的某种途径。基于不同的数学模型，提高地震数据分辨率的方法也表现为不同的特征(Wiggins，1985 ；Velis, 2008)，反问题理论越来越多地被应用于提高分辨率方法中(王家映，2002)，反演结果依赖于所选择的正演模型中相应的匹配观测数据中近似误差的分布以及对应评价准则的选择(张宏兵，尚作萍，杨长春等，2005)，忽略系统误差的具体形式而单纯去“匹配”的做法是存在“风险”的(Parker, 1977)。另外，约束条件的建立往往其“指向化”特征太明显，即不同的约束条件下反演结果会表现不同的特征，导致反演得到的模型参数估计结果是无法预测，并且是非常不稳定的(Sacchi, Velis和Cominguez, 1992)。从理论上来讲,先验条件是人们对所反演模型参数特征的一个初始的猜测，这个初始的猜测应该近似于真解，即一般情况下应该是符合实际物理规律的。如果这些条件过于“指向化”，这样得到的结果与真实情况往往是有很大的偏差的(Ulrych，2001)。 综上所述，现有地震数据高分辨率处理方法往往在获得高分辨处理结果的同时牺牲了很多有用的地震信息，实际上这些信息存在于地震道波形的微小变化之中，现有的地震数据高分辨方法特别是反演类算法对岩性的微小变化不敏感，高分辨处理后的地震资料在获得了地震数据视觉上的高分辨后丢失了很多有用的储层信息。特别是对于页岩储层的多层系特征以及岩性变化特点，这一大类高分辨率处理方法并不适用。 2002年，Vetterli提出了有限更新率信号的采样和重建理论，这一理论解释了有限带宽情况下反射系数重建的可能性。Vetterli (2002)在重建脉冲序列信号时所用的方法正是参数化谱分析中的一系列方法，比如Prony方法(Kay和Marple，1981)，子空间方法(Kung等，1983 ；Maravi, 2004)等。基于模型的参数类方法的估计性能都十分依赖于假设的数学模型和实际过程的符合程度，在模型不符时常常会得到失败的估计结果。实际上，针对于地震道反射系数的频谱模型特征，非参数化方法表现出更加抗噪的特性(Stoica和Moses，1997)。首先提出的最大似然(Capon)估计方法(Capon，1969)初始主要用于地震数据阵列信号处理中地震波的传播特征的研究，先被解释为一种基于最大似然的谱估计方法。后来，这一方法被用于对有限采样长度的不同频率的正弦信号分别来估计其振幅和相位，估计时首先模拟数据中噪声和相对干扰的协方差，然后设计相应的自适应滤波器使其能够无损地滤过特定的正弦信号的同时又尽可能地压制噪声和其它频率正弦的干扰，从而得到准最大似然估计结果。实验表明Capon估计方法在具有较高的估计稳定性的同时，其对频率接近的正弦信号分量有很高的分辨能力，且在某种意义上该方法可以在统计稳定性和频率分辨率之间进行调节，但问题在于其估计结果有偏差(有偏差原因在于对干扰协方差统计特征的粗略模拟)。1996年，Li和Stoica提出振幅与相位估计方法，其在目标探测和合成孔径雷达(SAR)成像中有广泛的应用，通过对干扰协方差的模拟设计自适应滤波器组来估计目标。相对于Capon方法的有偏差估计特性，振幅与相位估计方法是一种近无偏差估计，对干扰协方差的模拟更加准确。 因此，我们将这种思想应用于反射系数的频谱模型，对地下的反射幅度进行合理估计，在没有加入人为的先验信息的情况下，通过设计自适应的滤波器组压制相对干扰，可以得到高分辨率的反射系数幅度的估计结果。处理后的剖面或者三维道集资料在表现为较高分辨率较宽频带宽度的同时，同时表现为更准确的反射系数的幅值估计，提供了更丰富的幅度和相位信息，为进一步的针对页岩气储层特征的解释和评价工作提供资料保证。
技术实现思路
针对现有地震数据提高分辨率方法在频带拓宽效果以及资料保真中存在的问题，本专利技术以获得高保真度和高分辨率的地震道集为目标，借鉴和利用了国际上先进的统计信号处理方法，基于自适应滤波理论对地震资料进行高分辨处理。在高保真度的要求下，处理得到高分辨的数据体，数据频带得到拓宽，并且包含更高频更细节的地震信息。在此基础上，进行高分辨的属性提取和解释，进行页岩气储层分布的识别和储层地质条件的预测评价，以形成一整套的，针对页岩气藏的，涵盖页岩储层资料的地球物理方法的工作流程和技术体系。 针对页岩储层密度、速度的变化幅度非常小带来的弱反射特征，自适应高分辨率处理方法基于自适应统计滤波理论，设计相应的自适应滤波器组，在最大限度地保持原有地震资料信息的前提下，对地震道集进行处理，估计地下反射幅度。 在地震记录的有效频带范围内，用地震记录的频谱本文档来自技高网...

【技术保护点】
一种地震资料自适应高分辨处理方法，其特征在于：方法包括如下步骤：1)对页岩探区采集得到的地震资料进行常规处理和偏移，得到叠后或者偏移后道集；2)针对二维或者三维道集，选取目标页岩层层位，提取子波，并计算子波的频谱，并根据的资料特征，信噪比选取地震资料有效带宽范围；3)抽取单道地震记录，计算其频谱，在有效带宽范围内除以子波的频谱从而得到反射系数的部分谱序列；4)基于反射系数的部分谱序列与时间域反射脉冲的对应关系，选定某一时间位置，模拟相对干扰的统计特征，相对干扰由部分谱中其他时间位置反射脉冲对应的正弦和噪声构成，从而设计自适应滤波器，设计准则为使其滤波输出尽量接近所选定时间对应的频率域正弦，同时压制相对干扰，然后应用所设计的滤波器对部分谱序列进行滤波；5)对滤波后得到的近似正弦序列即对应于所选时间位置估计其幅度，根据其与时间域反射脉冲的对应关系，其估计即为对应所选时间位置的反射幅度估计；6)选择不同的时间位置，即针对不同的时间位置并设计相对应的不同的自适应滤波器来估计最终得到全部的反射系数序列；7)结合测井资料，井旁道处理，对剖面的应用效果进行评价，从而对滤波器组的相关参数进行调整优化效果，最后逐道进行步骤3‑6即可得到对整个地震数据的高保真和高分辨处理结果；8)在得到的更宽频带的道集的基础上，进行进一步的构造，岩性，微裂缝发育解释和评价工作，寻找页岩气勘探开发有利区域；9)结合已有探井，地质各方面资料，提供新的井位设计依据。...

【技术特征摘要】

【专利技术属性】
技术研发人员：郭锐，余刚，张宇生，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：河北;13