[發明專利]基于入射角的AVO近似公式進行屬性提取的方法

在審

• 申請號: 201510348231.6
• 申請日: 2015-06-23
• 公開（公告）號: CN104977610A
• 公開（公告）日: 2015-10-14
• 發明（設計）人: 張廣智;鄭靜靜;印興耀
• 申請（專利權）人: 中國石油大學（華東）
• 主分類號: G01V1/28
• 分類號: G01V1/28
• 代理公司: 山東濟南齊魯科技專利事務所有限公司 37108
• 代理人: 鄭向群

1.基于入射角的AVO近似公式進行屬性提取的方法，其特征在于：實現過程包括：1）基于入射角的AVO近似公式的提出；2）基于入射角的AVO近似公式的理論推導；3）應用常規方法提取P、G，Rs的屬性；4）將屬性歸一化；5）由井數據或者研究區域的經驗公式確定系數b值，然后應用公式（14）????????????????????????????????????????????????和公式（15）對常規方法提取的屬性進行修正；6）將修正結果與井數據進行對比分析，根據分析結果對修正結果進行校正，得到最終的屬性；其中，P是由零炮檢距構成的地震道，即縱波疊加道，代表對反射界面兩側的波阻抗變化的響應；斜率G為梯度疊加道，代表對橫波速度、縱波速度和體密度變化的響應，是由入射角近似公式計算得到的，G0代表對橫波速度、縱波速度和體密度變化的響應，由常規近似公式計算的得到；Rs表示由入射角近似公式計算的橫波阻抗反射系數，為由常規近似公式計算的橫波阻抗反射系數；

具體的求取方法為：?????????????

1）由Aki和Shuey近似可得：

??????（1）

其中，

，，

，，，

R（θ）表示縱波反射系數，和分別為反射界面上、下地層介質的密度，和分別是反射界面上下上、下地層的縱波速度，和分別是上、下地層的橫波速度，和分別是地層的入射角和透射角

2）根據Snell定律有：

????????（2）

對公式（1）左邊項進行泰勒展開，及

???????（3）

將公式（3）代入公式（2）可得：

??????（4）

又因為：

??????（5）

所以可得：

??????（6）

當采用一階近似時，上式可以寫成：

????（7）

令

則在此基礎上有

而當和相差不大時，

，

3）這樣可以得到基于入射角的AVO近似公式：

????（8）

4）Gardner等人于1974提出的密度與縱波速度之間的經驗公式的形式為[6]：

其中，a和b為系數(一般取a=0.31,b=0.25)

由上式可得

則

????（9）

又因為

代入上式得

所以

即：

將上式代入公式（8）得：

（10）

即:

???????(11)

其中，

，

5）一般由Shuey近似計算截距和梯度屬性P和G（進一步可計算橫波阻抗反射系數Rs等屬性）的公式為：

????（12）

其中，是根據平均角計算出的梯度

6）在實際應用公式（11）計算屬性時，由于求取角度不方便，常常用入射角代替進行估算，即應用下式計算P和Rs：

????（13）

比較式(11)和式(13)可得：

??????（14）

所以，計算橫波阻抗反射系數Rs的公式為：

??????（15）

這里的P和R_S0由公式（13）計算得到。