[發明專利]一種油氣儲層流體活動性預測方法

有效

• 申請號: 201810497306.0
• 申請日: 2018-05-22
• 公開（公告）號: CN108952695B
• 公開（公告）日: 2021-11-26
• 發明（設計）人: 李福來;盧莉
• 申請（專利權）人: 中國石油大學（華東）
• 主分類號: E21B49/00
• 分類號: E21B49/00;G06F30/20;G06F17/10
• 代理公司: 北京匯澤知識產權代理有限公司 11228
• 代理人: 程殿軍

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的值，得到流體活動性值。