确定反射系数频散、反射能量及识别气层的方法及装置

摘要

本发明公开了一种确定反射系数频散、反射能量及识别气层的方法及装置，包括：对地震数据进行保幅处理，并提取角道集；根据地震资料的品质选择合适角度的角道集；通过叠前反演方法确定纵波速度、横波速度和密度；确定地层纵波品质因子和横波品质因子；确定反射系数频散；确定反射能量系数；确定反射能量；根据钻井和试井情况，确定门槛值；将反射能量低于门槛值的地层识别为气层。采用本发明可以将地层中含气后引起的速度频散及能量变化通过反射系数准确地反映出来，准确确定反射能量，据此判别地下储层的含气性，有效减低气层识别的多解性。

主权项

一种反射系数频散确定方法，其特征在于，包括：对地震数据进行保幅处理，并提取角道集；根据地震资料的品质选择合适角度的角道集；通过叠前反演方法确定纵波速度、横波速度和密度；确定地层纵波品质因子和横波品质因子；按如下公式确定反射系数频散Rf：Rf=R(θ,Q,f)=R(θ)+R(Qp,f)·sec2θ+R(Qp,Qs,f)·sin2θ其中：R(θ)=A+Bsin2θ+Csin2θtg2θ A = 1 2 ( ΔV p V p ‾ + Δρ ρ ‾ ) , B = Δ V p 2 V p ‾ - 4 V s 2 V p 2 Δ V s V s ‾ - 2 V s 2 V p 2 Δρ ρ ‾ , C = 1 2 Δ V p V p ‾ Δ V p V p ‾ = 2 V p ′ - V p V p ′ + V p , Δρ ρ ‾ = 2 ρ ′ - ρ ρ ′ + ρ , V s 2 V p 2 = 1 2 ( V s ′ 2 V p ′ 2 + V s 2 V p 2 ) Vp和Vp'分别为地下任意两层地层的纵波速度，Vs和Vs'分别为地下任意两层地层的横波速度，ρ和ρ'分别为地下任意两层地层的密度，θ为入射角； R ( Q p , f ) = 1 2 π ( 1 Q p - 1 Q p ′ ) · ln ( f f 0 ) + i · 1 4 ( 1 Q p - 1 Q p ′ ) R ( Q p , Q s , f ) = ( 1 2 B - C ) ( f f 0 ) 2 π ( 1 Q s - 1 Q p ) · cos ( 1 Q s - 1 Q p ) + ( 1 2 B - C ) ( f f 0 ) 2 π ( 1 Q s ′ - 1 Q p ′ ) · cos ( 1 Q s ′ - 1 Q p ′ ) + ( C - B ) - 4 D π ( 1 Q s - 1 Q s ′ ) · ln ( f f 0 ) + i · [ ( 1 2 B - C ) ( f f 0 ) 2 π ( 1 Q s - 1 Q p ) · sin ( 1 Q s - 1 Q p ) + ( 1 2 B - C ) ( f f 0 ) 2 π ( 1 Q s ′ - 1 Q p ′ ) · sin ( 1 Q s ′ - 1 Q p ′ ) - 2 D ] Qp和Qp'分别为地下任意两层地层的纵波品质因子，Qs和Qs'分别为地下任意两层地层的横波品质因子，f为分析频率，f0为参考频率，D=Vs2/Vp2。