基于压缩感知的高斯频率域提高薄互层分辨率的方法技术

本发明专利技术提供一种基于压缩感知的高斯频率域提高薄互层分辨率的方法，包括：进行薄互层正演模拟，分析其地震响应特性；利用傅里叶变换把地震资料从时间域变换到频率域；引入对信号具有去噪和平滑作用的高斯函数，把频率域地震资料变换到高斯频率域；对原始地震资料利用多道复赛谱技术提取地震子波；利用提取的地震子波和高斯频率域地震资料进行基于压缩感知的稀疏反演，得到稀疏反射系数；基于提取的地震子波进行拓频处理得到宽频地震子波，把宽频地震子波与稀疏反射系数进行褶积，获得高分辨率地震资料。该基于压缩感知的高斯频率域提高薄互层分辨率的方法可以提高地震资料分辨率，解决薄互层识别的难题。

A method to improve the resolution of thin interbed based on compressed sensing in Gaussian frequency domain

【技术实现步骤摘要】
基于压缩感知的高斯频率域提高薄互层分辨率的方法
本专利技术涉及地震资料处理
，特别是涉及到一种基于压缩感知的高斯频率域提高薄互层分辨率的方法。
技术介绍
随着勘探技术的突破前勘探目标从常规油气藏转到复杂油气藏，尤其以薄互层为主。常规地震资料处理技术对薄互层的识别能力有限，无法对地震资料进行精细描述。常规的拓频技术方法存在噪声问题，而高斯函数利用“加窗”思想处理信号，去除信号高频随机噪声和实现信号平滑。将这一思想借鉴到基于压缩感知的稀疏反演处理中，可以有效保障信号的信噪比和光滑作用，利于储层的精细描述。再利用子波频带补偿拓展方法，提高储层的分辨率，对含薄互层储层进行精细评价具有至关重要的作用。为此我们专利技术了一种新的基于压缩感知的高斯频率域提高薄互层分辨率的方法，解决了以上技术问题。
技术实现思路
本专利技术的目的是提供一种能够更好地压制地震资料的高频随机噪声，提高地震资料的信噪比的基于压缩感知的高斯频率域提高薄互层分辨率的方法。本专利技术的目的可通过如下技术措施来实现：基于压缩感知的高斯频率域提高薄互层分辨率的方法，该基于压缩感知的高斯频率域提高薄互层分辨率的方法包括：步骤1，进行薄互层正演模拟，分析其地震响应特性；步骤2，利用傅里叶变换把地震资料从时间域变换到频率域；步骤3，引入对信号具有去噪和平滑作用的高斯函数，把频率域地震资料变换到高斯频率域；步骤4，对原始地震资料利用多道复赛谱技术提取地震子波；步骤5，利用提取的地震子波和高斯频率域地震资料进行基于压缩感知的稀疏反演，得到稀疏反射系数；步骤6，基于提取的地震子波进行拓频处理得到宽频地震子波，步骤7，把宽频地震子波与稀疏反射系数进行褶积，获得高分辨率地震资料。本专利技术的目的还可通过如下技术措施来实现：在步骤1，根据地震资料进行正演模型分析，研究以薄互层为储层地震资料的基本特性，给实际薄互层地震资料提高分辨率提供理论指导；特别是通过对地震资料的频谱信息研究，得到地震资料的频带信息范围，为子波拓频处理进行指导。在步骤2，利用傅里叶变换，把地震资料y从时域变换到频率进行处理，其傅里叶变换结果为：Y(w)＝F[y(t)]＝∫y(t)e-iwtdt式中，F表示傅里叶变换，Y(w)是频谱，t表示时间。在步骤3，加窗高斯傅里叶变换为：G(w)＝∫y(t)e-iwtdt×∫g(t)e-iwtdt式中，表示高斯窗函数，σ表示曲线宽度，控制有效信号频段。在步骤3，对于含噪地震记录来说，当信噪比低的时候，σ应选大一些，这样基于高斯函数对高频噪声具有很好的作用，对信号具有更明显的平滑效果。在步骤4，根据原始地震资料利于多道复赛谱技术，对地震资料进行子波提取，多道复赛谱利于加权平均的思想把地震资料的频谱取对数得到复赛谱的结果叠加在一起，再利用矩形窗提取合适的子波复赛谱结果，最后反变换得到时域子波；其复赛谱为：式中，Ln表示对数函数，i表示第i道地震记录，N表示总采样点个数，Y表示地震频谱，P表示地震复赛谱。在步骤4，对于模型数据，不需要提取子波，直接利用正演子波进行后续处理；对于实际资料，由于子波未知，需要基于多道复赛谱提取子波，窗函数不易过大或过小，只要包含主要信息即可。在步骤5，利用基于压缩感知的稀疏反演，根据地震子波和高斯频率域地震结果，获得稀疏反射系数，其关键是其L1范数的约束可以促进反射系数的稀疏表示；其公式为：频率域：高斯频率域：式中，x表示稀疏反射系数，W表示子波对角矩阵和部分傅里叶矩阵的组合，这里采用基追踪去噪算法求解，λ表示正则项因子，控制结果的稀疏效果，而G表示高斯频率域的地震谱。在步骤5，基于基追踪去噪算法求解目标函数的λ取值为最大振幅值的10％，这样既可以保证弱信号也能有效去除噪声干扰。在步骤6，利用提取的地震子波w(t)进行拓频处理，根据频带信息，拓宽子波的频带信息，得到宽频地震子波w′(t)，这里拓频处理，根据子波频带信息比较，选取合适的宽频带即可。在步骤7，利用宽频地震子波w′(t)与稀疏反射系数x(t)褶积，获得新的高分辨率地震记录S′(t)；其褶积公式为：S′(t)＝w′(t)*x(t)＝F-1[F(w(t))×F(x(t))]式中，F表示傅里叶变换，F-1表示傅里叶反变换，w(t)表示宽带子波。本专利技术中的基于压缩感知的高斯频率域提高薄互层分辨率的方法，主要用于地震资料目标处理以及提高地震资料分辨率，在噪声存在的情况下无法准确反演反射系数，并不能真实薄互层等的有效信息。为了在提高分辨率的同时保证信噪比，本专利技术利用高斯函数的去噪和平滑特征，应用到频率域地震资料进行信号去噪平滑处理，能够更好地压制地震资料的高频随机噪声，提高地震资料的信噪比，使地震资料特征更明显，特别是利于薄互层中弱信号和噪声的区分等。附图说明图1为本专利技术的基于压缩感知的高斯频率域提高薄互层分辨率的方法的一具体实施例的流程图；图2为本专利技术的一具体实施例中一维反射系数模型的示意图；图3为本专利技术的一具体实施例中一维合成地震记录的示意图；图4为本专利技术的一具体实施例中25Hz零相位雷克子波的示意图；图5为本专利技术的一具体实施例中25Hz零相位雷克子波对应的频谱的示意图；图6为本专利技术的一具体实施例中10％随机噪声的地震记录的示意图；图7为本专利技术的一具体实施例中含噪地震记录对应的频谱的示意图；图8为本专利技术的一具体实施例中加高斯窗处理后的地震记录(δ＝1)的示意图；图9为本专利技术的一具体实施例中地震记录处理对应的高斯域频谱(δ＝1)的示意图；图10为本专利技术的一具体实施例中加高斯窗处理后的地震记录(δ＝3)的示意图；图11为本专利技术的一具体实施例中地震记录处理对应的高斯域频谱(δ＝3)的示意图；图12为本专利技术的一具体实施例中基于压缩感知的频率域稀疏反演的示意图；图13为本专利技术的一具体实施例中基于压缩感知的高斯频率域稀疏反演(δ＝1.2)的示意图；图14为本专利技术的一具体实施例中宽带50Hz雷克子波的示意图；图15为本专利技术的一具体实施例中宽带50Hz雷克子波的频谱的示意图；图16为本专利技术的一具体实施例中基于压缩感知的频域稀疏反演拓频处理的示意图；图17为本专利技术的一具体实施例中频域处理对应的频谱的示意图；图18为本专利技术的一具体实施例中基于压缩感知的高斯频域稀疏反演拓频处理的示意图；图19为本专利技术的一具体实施例中高斯频域处理对应的频谱的示意图；图20为本专利技术的一具体实施例中二维反射系数模型的示意图；图21为本专利技术的一具体实施例中二维合成地震记录的示意图；图22为本专利技术的一具体实施例中25Hz的雷克子波的示意图；图23为本专利技术的一具体实施例中25Hz的雷克子波的频谱的示意图；图24为本专利技术的一具体实施例中含10％随机噪声二维地震模型的示意图；图25为本本文档来自技高网...

【技术保护点】
1.基于压缩感知的高斯频率域提高薄互层分辨率的方法，其特征在于，该基于压缩感知的高斯频率域提高薄互层分辨率的方法包括：/n步骤1，进行薄互层正演模拟，分析其地震响应特性；/n步骤2，利用傅里叶变换把地震资料从时间域变换到频率域；/n步骤3，引入对信号具有去噪和平滑作用的高斯函数，把频率域地震资料变换到高斯频率域；/n步骤4，对原始地震资料利用多道复赛谱技术提取地震子波；/n步骤5，利用提取的地震子波和高斯频率域地震资料进行基于压缩感知的稀疏反演，得到稀疏反射系数；/n步骤6，基于提取的地震子波进行拓频处理得到宽频地震子波，/n步骤7，把宽频地震子波与稀疏反射系数进行褶积，获得高分辨率地震资料。/n

【技术特征摘要】
1.基于压缩感知的高斯频率域提高薄互层分辨率的方法，其特征在于，该基于压缩感知的高斯频率域提高薄互层分辨率的方法包括：
步骤1，进行薄互层正演模拟，分析其地震响应特性；
步骤2，利用傅里叶变换把地震资料从时间域变换到频率域；
步骤3，引入对信号具有去噪和平滑作用的高斯函数，把频率域地震资料变换到高斯频率域；
步骤4，对原始地震资料利用多道复赛谱技术提取地震子波；
步骤5，利用提取的地震子波和高斯频率域地震资料进行基于压缩感知的稀疏反演，得到稀疏反射系数；
步骤6，基于提取的地震子波进行拓频处理得到宽频地震子波，
步骤7，把宽频地震子波与稀疏反射系数进行褶积，获得高分辨率地震资料。


2.根据权利要求1所述的基于压缩感知的高斯频率域提高薄互层分辨率的方法，其特征在于，在步骤1，根据地震资料进行正演模型分析，研究以薄互层为储层地震资料的基本特性，给实际薄互层地震资料提高分辨率提供理论指导；特别是通过对地震资料的频谱信息研究，得到地震资料的频带信息范围，为子波拓频处理进行指导。


3.根据权利要求1所述的基于压缩感知的高斯频率域提高薄互层分辨率的方法，其特征在于，在步骤2，利用傅里叶变换，把地震资料y从时域变换到频率进行处理，其傅里叶变换结果为：
Y(w)＝F[y(t)]＝∫y(t)e-iwtdt
式中，F表示傅里叶变换，Y(w)是频谱，t表示时间。


4.根据权利要求3所述的基于压缩感知的高斯频率域提高薄互层分辨率的方法，其特征在于，在步骤3，加窗高斯傅里叶变换为：
G(w)＝∫y(t)e-iwtdt×∫g(t)e-iwtdt
式中，表示高斯窗函数，σ表示曲线宽度，控制有效信号频段。


5.根据权利要求4所述的基于压缩感知的高斯频率域提高薄互层分辨率的方法，其特征在于，在步骤3，对于含噪地震记录来说，当信噪比低的时候，σ应选大一些，这样基于高斯函数对高频噪声具有很好的作用，对信号具有更明显的平滑效果。


6.根据权利要求1所述的基于压缩感知的高斯频率域提高薄互层分辨率的方法，其特征在于，在步骤4，根据原始地震资料利于多道复赛谱技术，对地震...

【专利技术属性】
技术研发人员：王延光，汪浩，刘立彬，梁鸿贤，毕丽飞，程远峰，揭景荣，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司胜利油田分公司物探研究院，
类型：发明
国别省市：山东;37