基于小波包变换的频谱成像方法技术

本发明专利技术公开了一种基于小波包变换的频谱成像方法，该方法包括以下步骤：在偏移后的三维地震数据体上进行加密小波包变换，将地震记录分解为按频带从小到大顺序排列的分频地震记录；根据小波包变换结果计算分频瞬时振幅谱、分频瞬时相位谱、瞬时平均频率、瞬时带宽、瞬时吸收系数五种三维属性数据体；在地震解释层位的控制下，进行属性切片。本发明专利技术通过复赛谱技术消除反射系数谱对地震子波谱估算的影响，在此基础上估算的瞬时吸收系数更为合理；三种不同频率属性的融合显示利用颜色混合特有的效果来突出目标体，能够有效改善单一频率属性对目标的分辨效果。

Spectral imaging method based on wavelet packet transform

The invention discloses a spectrum imaging method based on wavelet packet transform, the method comprises the following steps: encryption and wavelet packet transform in seismic data volume after offset, the seismic record is decomposed into frequency band according to seismic records from small to large order; according to the calculation results of wavelet packet transform frequency instantaneous amplitude spectrum, instantaneous frequency phase spectrum, instantaneous frequency, instantaneous bandwidth, instantaneous average absorption coefficient of five 3D attribute data; in the control of seismic interpretation horizon, attribute slice. The invention eliminates the reflectance spectrum of seismic wavelet spectrum estimation by cepstrum technique, instantaneous estimation based on absorption coefficient is more reasonable; three different frequency attribute fusion display to highlight the target utilization effect of color mixing characteristic, can effectively improve the resolution of single frequency attribute on the target.

【技术实现步骤摘要】
基于小波包变换的频谱成像方法
本专利技术涉及石油勘探开发技术，尤其涉及一种基于小波包变换的频谱成像方法。
技术介绍
根据薄层调谐原理，不同频率振幅信息反映了不同厚度地层的调谐特征，因此地震记录实质上是不同地层厚度信息的叠合，也就是不同地质体的综合反映。对于薄砂层、河道和砂坝等隐蔽性油气藏，由于受其它地质体信息的干涉，在地震数据体上难以清晰识别，而频谱成像通过时频分析技术对地震数据进行频谱分解，对储层响应较好的频率成分进行成像，从而有利于目标地质体的识别。频谱成像技术的技术核心是时频分析技术。目前，商业软件采用的时频分析技术基本上是以短时窗傅立叶变换和小波变换为主。众所周知，虽然短时窗傅立叶变换能对信号的局部特征进行分析，但它存在时窗固定，适应性较差的弊病。小波变换在短时窗傅立叶变换基础上进行了发展，它克服了短时窗傅立叶变换对时窗的限制，但它对频率的划分比较粗糙，在成像方面精度不够，从而影响了该项技术的推广应用。吸收系数往往通过地震子波在介质中传播时的高频衰减进行估算，而已知的只有地震记录的振幅谱，它是地震子波和反射系数的综合结果。在估算地震子波谱时应消除反射系数谱的影响。单一频率属性通常用于识别特定厚度的储层，而储层厚度在空间上是变化的，因此单一频率属性图在储层预测时存在一定局限。
技术实现思路
本专利技术要解决的技术问题在于针对现有技术中的缺陷，提供一种基于小波包变换的频谱成像方法，用于解决针对短时窗傅立叶变换时窗固定和小波变换对频率的划分比较粗糙的问题。本专利技术解决其技术问题所采用的技术方案是：一种基于小波包变换的频谱成像方法，包括以下步骤：1)在偏移后的三维地震数据体上进行加密小波包变换，将地震记录分解为按频带从小到大顺序排列的分频地震记录；2)根据小波包变换结果计算分频瞬时振幅谱、分频瞬时相位谱、瞬时平均频率、瞬时带宽、瞬时吸收系数五中三维属性数据体；3)在地震解释层位的控制下，进行属性切片；包括两种可选切片方式：以某一时间间隔δt沿层切片或将两个层位间数据剖分N等份切片；4)每种属性可以单独分析，也可以对三种不同频率属性进行融合分析。按上述方案，所述对偏移归位后地震数据进行按频带排序的加密小波包变换；定义共轭滤波器组hn和gn，满足：gk＝(-1)kh1-k(3)其中，n,k,l∈Z，离散地震信号sk在小波包基下的按频带排序的加密塔式分解方法为：如果n为偶数，则如果n为奇数，则式中，n＝0,1,…2l-1；l为小波包分解层数；为小波包分解系数，其中：(sk为地震信号)。通过以上变换，将地震记录分解为按频带从小到大顺序排列的分频地震记录。按上述方案，步骤2)中分频瞬时振幅谱和分频瞬时相位谱参数确定方法如下：设S(t,a)是信号h(t)的小波包变换结果，其中a是频率通道号。则其解析信号为X(t,a)＝S(t,a)+iH(t,a)(6)其中，H(t,a)是S(t,a)的希尔伯特变换，i为虚数单位，即这里，*表示褶积运算；将X(t,a)写成指数形式为X(t,a)＝A(t,a)eiφ(t,a)(8)其中，A(t,a)为瞬时振幅谱，φ(t,a)为瞬时相位谱，则φ(t,a)＝arg(X(t,a))+2kπ，k＝0,±1,±2,…(10)arg(X(t,a))为解析信号X(t,a)的辐角主值。其范围为-π＜arg(X(t,a))＜π设k＝0，根据解析信号所在的象限位置按以下变换来计算瞬时相位谱：在第一象限：S(t,a)＞0，H(t,a)＞0，在第二象限：S(t,a)＜0，H(t,a)＞0，在第三象限：S(t,a)＜0，H(t,a)＜0，在第四象限：S(t,a)＞0，H(t,a)＜0，按上述方案，步骤2)中从瞬时振幅谱A(t,a)可以得到功率谱密度函数P(t,a)：P(t,a)＝A2(t,a)(11)因为a对应着一定频带范围的中心频率，在这里不妨将P(t,a)记为P(t,ω)。按上述方案，步骤2)中瞬时平均频率和瞬时带宽的确定方法如下:瞬时平均频率定义为在某一时刻功率谱密度函数的期望值。即：该值能够反映地层厚度变化和高频能量的衰减情况。瞬时带宽定义为：瞬时带宽能够指示总体吸收效果。按上述方案，步骤2)中瞬时吸收系数的确定方法如下:瞬时吸收系数的提取包括两个过程：1、地震子波频谱的提取；2、吸收系数的提取。下面分述这两个过程：1)地震子波频谱的提取设地震记录为子波与反射系数的褶积，即s(t)＝w(t)*f(t)(14)其中，对上式进行傅里叶变换得:S(w)＝W(w)F(w)(15)取对数,得:ln(S(w))＝lnW(w)+lnF(w)(16)对上式进行反傅里叶变换得复赛谱:地震子波相对于反射系数序列而言是比较平滑的、变化较慢的低频部分,所以它的复赛谱主要集中在时间原点附近，而反射系数序列是按时间均匀分布的，通过低通滤波获得地震子波的复赛谱，通过上述过程的逆过程可以由地震子波的复赛谱得到其振幅谱A(w)；2)吸收系数的提取在某一时刻，地震子波的振幅谱通常可以用图3表示，地震子波的振幅谱在谱的高频段(ωd～ωN，ωd-主频，ωN-Nyquist频率)谱近似满足指数函数(图3粗黑线条所示)A(ω)＝Cexp(-Q(t)ω)(18)其中：C为常数，Q(t)为瞬时吸收系数，ω为角频率。因此，可以对振幅谱的高频段进行指数拟合来模拟吸收系数Q(t)值。式(18)经过对数变换可以化为线性拟合。ln(A(ω))＝lnC-Q(t)ω(19)令离散化后的地震子波高频段的(频率，对数振幅谱)对为(f1,A1)，(f2,A2)，(f3,A3)，…，(fn,An)，则其中：fi为频率，Ai为对应的对数振幅谱。这样就得到了瞬时吸收系数。本专利技术产生的有益效果是：本专利技术通过采用按频带排序的加密小波包变换技术，提高了频带划分精度；研究了分频瞬时吸收系数反演技术，在估算地震子波谱方面方法更为合理；通过将三种不同频率属性叠合显示在一张图上，提升了通过频谱成像技术对储层的识别能力。附图说明下面将结合附图及实施例对本专利技术作进一步说明，附图中：图1是本专利技术实施例的方法流程图；图2是本专利技术的小波包分解算法示意图；图3是本专利技术的地震子波的振幅谱；图4是本专利技术实施例的沿层切片图；图5是本专利技术实施例三频率属性RGB融合显示图。具体实施方式为了使本专利技术的目的、技术方案及优点更加清楚明白，以下结合实施例，对本专利技术进行进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本专利技术，并不用于限定本专利技术。如图1所示，本专利技术实施例基于小波包变换的频谱成像方法主要包括以下步骤：S1：在偏移后三维地震数据体上进行按频带排序的加密小波包变换；定义共轭滤波器组hn和gn，满足：gk＝(-1)kh1-k(3)其中，n,k,l∈Z，离散地震信号sk在小波包基下的按频带排序的加密塔式分解算法为：如果n为偶数，则如果n为奇数，则式中，n＝0,1,…2l-1；l为小波包分解层数；为小波包分解系数，其中：(sk为地震信号)。通过以上变换，将地震记录分解为按频带从小到大顺序排列的分频地震记录。S2：计算分频瞬时振幅谱、分频瞬时相位谱、瞬时平均频率、瞬时带宽、瞬时吸收系数；分频瞬时振幅谱和分频瞬时相位谱:设S(t,a)是信号h(t)的小波包变换结果，其中a是频率通道号。则其解析信号为X(t,a)本文档来自技高网...

【技术保护点】
一种基于小波包变换的频谱成像方法，其特征在于，包括以下步骤：1)在偏移后的三维地震数据体上进行加密小波包变换，将地震记录分解为按频带从小到大顺序排列的分频地震记录；2)根据小波包变换结果计算分频瞬时振幅谱、分频瞬时相位谱、瞬时平均频率、瞬时带宽、瞬时吸收系数五种三维属性数据体；3)在地震解释层位的控制下，进行属性切片；包括两种可选切片方式：以某一时间间隔δt沿层切片或将两个层位间数据剖分N等份切片；4)每种属性可以单独分析，也可以对三种不同频率属性进行融合分析，获得成像结果。

【技术特征摘要】
1.一种基于小波包变换的频谱成像方法，其特征在于，包括以下步骤：1)在偏移后的三维地震数据体上进行加密小波包变换，将地震记录分解为按频带从小到大顺序排列的分频地震记录；2)根据小波包变换结果计算分频瞬时振幅谱、分频瞬时相位谱、瞬时平均频率、瞬时带宽、瞬时吸收系数五种三维属性数据体；3)在地震解释层位的控制下，进行属性切片；包括两种可选切片方式：以某一时间间隔δt沿层切片或将两个层位间数据剖分N等份切片；4)每种属性可以单独分析，也可以对三种不同频率属性进行融合分析，获得成像结果。2.根据权利要求1所述的基于小波包变换的频谱成像方法，其特征在于，步骤1)中所述按频带排序的加密小波包变换实现如下；定义共轭滤波器组hn和gn，满足：gk＝(-1)kh1-k(3)其中，n,k,l∈Z，离散地震信号sk在小波包基下的按频带排序的加密塔式分解方法为：如果n为偶数，则如果n为奇数，则式中，n＝0,1,…2l-1；l为小波包分解层数；为小波包分解系数，其中：sk为地震信号。3.根据权利要求1所述的基于小波包变换的频谱成像方法，其特征在于，步骤2)中所述瞬时吸收系数的确定方法如下:瞬时吸收系数的提取包括两个过程：1)地震子波频谱的提取；2)吸收系数的提取；1)地震子波频谱的提取设地震记录为子波与反射系数的褶积，即s(t)＝w(t)*f(t)(6)其中，对上式进行傅里叶变换得:S(w)＝W(w)F(w)(7)取对数,得:ln(S(w))＝lnW(w)+lnF(w)(8)对上式进行反傅里叶变换得复赛谱:地震子波相对于反射系数序列而言是比较平滑的、变化较慢的低频部分,所以它的复赛谱主要集中在时间原点附近，而反射系数序列是按时间均匀分布的，通过低通滤波获得地震...

【专利技术属性】
技术研发人员：周开明，刘贤红，陈鹏，
申请(专利权)人：中国石油化工股份有限公司江汉油田分公司物探研究院，
类型：发明
国别省市：湖北,42