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

摘要

本发明公开了一种提高滩坝砂储层地震识别性能与描述精度的方法，属于地震勘探领域。其是先采用三参数小波来提高地震分辨薄互层的能力，再通过多子波技术对地震数据进行变换，得到成果数据，后期利用沿层切片或者属性分析等技术，以提高滩坝砂储层的地震识别性能与描述精度。本发明能有效提高薄层识别能力，提高主频，展宽频带范围。经多子波变化后，能更加突出滩坝砂储层信息，经属性融合显示后，能大大提高后期储层描述的精度。

主权项

一种提高滩坝砂储层地震识别性能与描述精度的方法，其特征在于包括以下步骤：a、首先利用三参数小波对地震资料进行提高分辨率处理；三参数小波公式为 φ ( t , Λ ) = e - τ ( t - β ) 2 { p ( Λ ) [ cos ( σt ) - k ( Λ ) ] + iq ( Λ ) sin ( σt ) } 其中，t表示时间，i表示虚数，σ为小波函数的调制频率，τ为能量衰减因子，β为能量延迟因子，（σ,τ,β∈R，σ,τ≥0），Λ=(σ,τ,β)记参数σ，τ，β集合；上式中的p(Λ)，q(Λ)和k(Λ)分别为 k ( Λ ) = e - σ 2 4 τ | cos ( βσ ) + iq ( Λ ) p ( Λ ) sin ( βσ ) | p ( Λ ) = | 2 τ π | 1 4 | 4 | e - σ 2 2 τ - e - 3 σ 2 8 τ | cos 2 ( βσ ) + 1 - e - σ 2 2 τ | - 1 2 q ( Λ ) = | 2 τ π | 1 4 | 4 | e - σ 2 2 τ - e - 3 σ 2 8 τ | sin 2 ( βσ ) + 1 - e - σ 2 2 τ | - 1 2 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、对处理后的地震数据信号，首先估算滩坝砂储层位置的子波的频带范围[fmin,fmax]，则可求出中心频率fc和频带宽度fw，设子波在该频带内的能量占子波总能量的λ2，则 λ 2 = ∫ - f max f max | W ( f ) | 2 df - ∫ - f min f min | W ( f ) | 2 df ∫ - f N / 2 f N / 2 | W ( f ) | 2 df 将该式转化为求解以下方程Aw=λw其中，A={am,n}N×N，则得到 a m , n = sin [ 2 π ( f c + f w / 2 ) Δt ( m - n ) ] π ( m - n ) - sin [ 2 π ( f c - f w / 2 ) Δt ( m - n ) ] π ( m - n ) 其中，N为子波样点数，Δt为采样间隔；取λ接近于1所对应的特征向量w，记为{w}，{w}的频域能量限制在频带[fmin,fmax]以内，{w}为所要求取的多子波；d、求取了标准子波以后，利用所求得的子波与地震道做相关处理 C ( t , i ) = Σ t + T / 2 t + T / 2 s ( τ , i ) w ( τ - t ) dτ , i = 1,2 , . . . , N 其中s为地震道，w为所求子波，C(t,i)便是相关信号，即为最后多子波变换数据；e、对以上数据进行属性提取，对提取的属性进行多属性融合。