[发明专利]一种油气储层流体活动性预测方法

权利要求书

1.一种油气储层流体活动性预测方法，其特征在于，其包括的步骤是：

第一步，得到反射系数的线性近似表示式：

设有式(1-1)，

式(1-1)中，R为反射系数，F为流体活动性；ρ_f为流体密度；ω为地震波角频率；R₀、R₁分别是反射系数渐近表示的零阶项和一阶项，是储层岩石和流体性质相关的函数；

对式(1-1)做关于流体活动性F的一阶Taylor展开，得到反射系数R的线性近似表示式：

重新整理式(2)，将式(2)写作：

R(F,ω)＝A₀(F₀,ω)+B₀(F₀,ω)·F (3)，

式(3)中，F₀是储层中所含流体的流体活动性特征值，且A₀、B₀具体表示式如下：

第二步，做分频处理：将任意频率处的反射系数值表示为R_ωi＝[R_ωi(t₁)R_ωi(t₂)..R_ωi(t_n)]^T,其中，t₁,t₂,...,t_n表示不同时间深度，从1到n，利用式(3)计算任意频率处的反射系数值：

记A_ωi＝[A_0ωi A_0ωi L A_0ωi]^T，F＝[F(t₁) F(t₂)L F(t_n)]^T，式(5)写为：

R_ωi＝A_ωi+B_ωi·F (6)；

第三步，得到重组后的正演算子：将式(6)代入非平稳的地震褶积模型S(ω_i)＝W_i·R_ωi，移项整理得到：

S(ω_i)-S'(ω_i)＝W_i·B_ωi·F (7)，

式(7)中S(ω_i)＝W_i·R_ωi，S'(ω_i)＝W_i·A_ωi，最终得到重组后的正演算子如下：

第四步，在贝叶斯框架下构建目标函数：加入先验概率分布引入的稀疏正则化项和低频模型约束项，最终得到目标泛函为：

式(9)中，σ_n²和σ_F²分别为噪声分布和反演模型参数的协方差；ξ和D分别是低频先验和正则化矩阵；λ_l为模型约束系数；

将式(9)求取模型参数F的梯度，并令其为0，得到模型参数最终的非线性反演方程，使用最小二乘法求取F的值，得到流体活动性值。