[發明專利]非線性AVO反演的方法和裝置

在審

• 申請號: 201610140300.9
• 申請日: 2016-03-11
• 公開（公告）號: CN107179545A
• 公開（公告）日: 2017-09-19
• 發明（設計）人: 陳天勝;劉春園;劉韜;劉炯;季玉新
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發研究院
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京思創畢升專利事務所11218
• 代理人: 孫向民

1.一種非線性縱波AVO反演的方法，該方法包括：

得到疊前角道集數據，所述疊前角道集數據包括N個入射角θ_i對應的縱波反射系數R_pp(θ_i)，i＝1,2,...,N；

將分別展開為縱波速度變化率r_α、橫波速度變化率r_β、密度變化率r_ρ的泰勒級數，將其代入Zoeppritz方程組得到公式5：

R_pp(θ)＝R_pp1(θ)+R_pp2(θ)+...+R_ppK(θ)， 公式5

其中，θ表示任意入射角，R_ppk(θ)僅包含縱波速度變化率r_α、橫波速度變化率r_β、密度變化率r_ρ的k階項，k＝1,2,...,K，K>1，α₁、β₁、ρ₁分別表示界面上層的縱波速度、橫波速度和密度，α₂、β₂、ρ₂分別表示界面下層的縱波速度、橫波速度和密度；

基于疊前角道集數據和公式5反演r_α、r_β、r_ρ。

2.根據權利要求1所述的非線性縱波AVO反演的方法，其中，反演r_α、r_β、r_ρ包括：

對r_α、r_β、r_ρ進行泰勒級數展開，得到：

r_α＝r_α1+r_α2+...+r_αL， 公式9

r_β＝r_β1+r_β2+...+r_βL， 公式10

r_ρ＝r_ρ1+r_ρ2+...+r_ρL， 公式11

其中，r_αl、r_βl、r_ρl分別表示r_α、r_β、r_ρ的泰勒級數展開中的l階項；

將r_α、r_β、r_ρ的泰勒級數展開代入公式5，得到公式12：

R_pp(θ)＝f(r_α1、r_α2、……、r_αL、r_β1、r_β2、……、r_βL、r_ρ1、r_ρ2、……、r_ρL)； 公式12

令公式12左側和右側的同階項相等，并使用疊前角道集數據，得到r_αl、r_βl、r_ρl，l＝1,2,...,L，L>1；

基于公式9至公式11得到r_α、r_β、r_ρ。