[發明專利]AVO地震正演計算方法

有效

• 申請號: 201310454638.8
• 申請日: 2013-09-29
• 公開（公告）號: CN103487831A
• 公開（公告）日: 2014-01-01
• 發明（設計）人: 鄒文;張洞君;陶正喜;何光明;周晶晶;吳秋波
• 申請（專利權）人: 中國石油集團川慶鉆探工程有限公司地球物理勘探公司
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 北京銘碩知識產權代理有限公司 11286
• 代理人: 譚昌馳;張川緒

1.一種AVO地震正演計算方法，所述AVO地震正演計算方法包括：獲取目的層段的地質數據，其中，所述目的層段包括三個層，三個層中的第一層與三個層中的第二層之間形成第一界面，第二層與三個層中的第三層之間形成第二界面，

（1）利用獲取的地質數據，通過動態等效介質理論，計算在不同頻率下的三個層的縱波速度和橫波速度；

（2）利用第一層的縱波速度和橫波速度和第二層的縱波速度和橫波速度，計算第一界面在不同頻率不同入射角下的反射系數，作為第一界面的AVO地震反射系數；

（3）利用第二層的縱波速度和橫波速度和第三層的縱波速度和橫波速度，計算第二界面在不同頻率不同入射角下的反射系數，作為第二界面的AVO地震反射系數；

（4）在所述目的層段內沿測井深度方向選擇多個采樣點，對選擇的多個采樣點的每個采樣點執行以下步驟：利用三個層的縱波速度、第一界面的AVO地震反射系數、第二界面的AVO地震反射系數，通過波場延拓得到不同入射角的AVO地震疊前角道集。

2.如權利要求1所述的AVO地震正演計算方法，其中，計算在不同頻率下的三個層的縱波速度和橫波速度的步驟包括：構建彈性張量矩陣，利用彈性張量矩陣計算在不同頻率下的三個層的縱波速度和橫波速度。

3.如權利要求1所述的AVO地震正演計算方法，其中，利用如下公式計算第一界面和第二界面中的任意界面在不同頻率不同入射角下的反射系數：

R_P(ω,θ)=A(ω)+B(ω)sin²θ+C(ω)tan²θsin²θ

其中，

A(ω)=12[ΔVP(ω)VP(ω)+Δρρ],]]>

B(ω)=12ΔVP(ω)VP(ω)-4[VS(ω)VP(ω)]2ΔVS(ω)VS(ω)-2[VS(ω)VP(ω)]2Δρρ,]]>

C(ω)=12ΔVP(ω)VP(ω),]]>

式中，R_p(ω,θ)為所述任意界面在不同頻率不同入射角下的反射系數，ω為角頻率，ω=2πf，f為頻率，θ為入射角，V_P(ω)為所述任意界面之上的層與所述任意界面之下的層的縱波速度的平均，V_S(ω)為所述任意界面之上的層與所述任意界面之下的層的橫波速度的平均，ΔV_P(ω)為所述任意界面之上的層與所述任意界面之下的層的縱波速度之差，ΔV_S(ω)為所述任意界面之上的層與所述任意界面之下的層的橫波速度之差，ρ為所述任意界面之上的層與所述任意界面之下的層的密度的平均，Δρ為所述任意界面之上的層與所述任意界面之下的層的密度之差。