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

摘要

本发明涉及勘探开发技术领域，具体公开了一种油气储层流体活动性预测方法，其包括的步骤是：第一步，得到反射系数的线性近似表示式；第二步，做分频处理；第三步，得到重组后的正演算子；第四步，在贝叶斯框架下构建目标函数：加入先验概率分布引入的稀疏正则化项和低频模型约束项，最终得到目标泛函，将目标泛函求取模型参数F的梯度，并令其为0，得到模型参数最终的非线性反演方程，使用最小二乘法求取F的值，得到流体活动性值。本发明是一种贝叶斯理论框架下流体活动性的频变反演预测方法，本发明方法能提高储层预测与流体识别的可靠性，使得计算流体活动性物理意义更加明确，结果更为可靠。

主权项

1.一种油气储层流体活动性预测方法，其特征在于，其包括的步骤是：第一步，得到反射系数的线性近似表示式：设有式(1‑1)，式(1‑1)中，R为反射系数，F为流体活动性；ρ_f为流体密度；ω为地震波角频率；R₀、R₁分别是反射系数渐近表示的零阶项和一阶项，是储层岩石和流体性质相关的函数；对式(1‑1)做关于流体活动性F的一阶Taylor展开，得到反射系数R的线性近似表示式：重新整理式(2)，将式(2)写作：R(F,ω)＝A0(F0,ω)+B0(F0,ω)·F  (3)，式(3)中，F0是储层中所含流体的流体活动性特征值，且A0、B0具体表示式如下：第二步，做分频处理：将反射系数表示为Rω＝[Rω1 Rω2 … RωM]T，M表示频率个数，利用式(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)＝Wi·Rωi，移项整理得到：S(ωi)‑S'(ωi)＝Wi·Bωi·F  (7)，式(7)中S(ωi)＝Wi·Rωi，S'(ωi)＝Wi·Aωi，最终得到重组后的正演算子如下：第四步，在贝叶斯框架下构建目标函数：加入先验概率分布引入的稀疏正则化项和低频模型约束项，最终得到目标泛函为：式(9)中，σn2和σF2分别为噪声分布和反演模型参数的协方差；ξ和D分别是低频先验和正则化矩阵；λl为模型约束系数；将式(9)求取模型参数F的梯度，并令其为0，得到模型参数最终的非线性反演方程，使用最小二乘法求取F的值，得到流体活动性值。