[发明专利]一种提高滩坝砂储层地震识别性能与描述精度的方法

权利要求书

1.一种提高滩坝砂储层地震识别性能与描述精度的方法，其特征在于包括以下步骤：

a、首先利用三参数小波对地震资料进行提高分辨率处理；三参数小波公式为

φ(t,Λ)=e-τ(t-β)2{p(Λ)[cos(σt)-k(Λ)]+iq(Λ)sin(σt)}]]>

其中，t表示时间，i表示虚数，σ为小波函数的调制频率，τ为能量衰减因子，β为能量延迟因子，（σ,τ,β∈R，σ,τ≥0），Λ=(σ,τ,β)记参数σ，τ，β集合；上式中的p(Λ)，q(Λ)和k(Λ)分别为

k(Λ)=e-σ24τ|cos(βσ)+iq(Λ)p(Λ)sin(βσ)|]]>

p(Λ)=|2τπ|14|4|e-σ22τ-e-3σ28τ|cos2(βσ)+1-e-σ22τ|-12]]>

q(Λ)=|2τπ|14|4|e-σ22τ-e-3σ28τ|sin2(βσ)+1-e-σ22τ|-12]]>

b、根据选取的三个小波参数Λ=(σ,τ,β),得到三参数小波函数φ(t,Λ)；设地震数据为X(t,i),i＝1,2，...,N，其中N为道数、t为地震道时间、i表示第i道，按道进行褶积处理得到

s(t,i)=X(t,i)*φ(t,Λ)，i=1,2，...,N

得到三参数小波处理后的数据；

c、对处理后的地震数据信号，首先估算滩坝砂储层位置的子波的频带范围[f_min,f_max]，则可求出中心频率f_c和频带宽度f_w，设子波在该频带内的能量占子波总能量的λ²，则

λ2=∫-fmaxfmax|W(f)|2df-∫-fminfmin|W(f)|2df∫-fN/2fN/2|W(f)|2df]]>

将该式转化为求解以下方程

Aw=λw

其中，A={a_m,n}_N×N，则得到

am,n=sin[2π(fc+fw/2)Δt(m-n)]π(m-n)-sin[2π(fc-fw/2)Δt(m-n)]π(m-n)]]>

其中，N为子波样点数，Δt为采样间隔；取λ接近于1所对应的特征向量w，记为{w}，{w}的频域能量限制在频带[f_min,f_max]以内，{w}为所要求取的多子波；

d、求取了标准子波以后，利用所求得的子波与地震道做相关处理

C(t,i)=Σt+T/2t+T/2s(τ,i)w(τ-t)dτ,i=1,2,...,N]]>

其中s为地震道，w为所求子波，C(t,i)便是相关信号，即为最后多子波变换数据；

e、对以上数据进行属性提取，对提取的属性进行多属性融合。