[發明專利]軸對稱介質各向異性參數確定方法及裝置

有效

• 申請號: 202010111495.0
• 申請日: 2020-02-24
• 公開（公告）號: CN113296153B
• 公開（公告）日: 2023-04-25
• 發明（設計）人: 王兆磊;公亭;王飛;顧小弟;孟曉夢
• 申請（專利權）人: 中國石油天然氣集團有限公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/30
• 分類號: G01V1/30
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 周曉飛;許曼

1.一種軸對稱介質各向異性參數確定方法，其特征在于，包括：

獲取軸對稱介質的地震資料，根據所述地震資料，反演迭代求取軸對稱介質的各向同性縱波速度，進行疊前深度偏移得到初始地震地層厚度；

根據所述初始地震地層厚度、實際地層厚度和巖石物理允許范圍，得到縱波速度調節參數；

根據所述縱波速度調節參數和所述各向同性縱波速度，得到初始各向異性縱波速度；

循環執行以下步驟，直至各向異性疊前深度偏移的共成像點偏移道集中同相軸拉平，確定新的縱波速度調節參數為軸對稱介質各向異性參數：

根據所述初始各向異性縱波速度和所述縱波速度調節參數，反演迭代得到新的各向異性縱波速度，進行各向異性疊前深度偏移得到新的地震地層厚度，根據所述新的地震地層厚度和所述實際地層厚度得到新的縱波速度調節參數；

將新的縱波速度調節參數和新的各向異性縱波速度作為下一次循環的初始各向異性縱波速度和縱波速度調節參數；

根據所述新的地震地層厚度和所述實際地層厚度得到新的縱波速度調節參數，包括：

根據所述新的地震地層厚度和所述實際地層厚度，得到井震誤差；

根據所述井震誤差和所述實際地層厚度，計算得到縱波速度調節參數的校正值；

根據所述縱波速度調節參數和所述縱波速度調節參數的校正值，累加求和得到新的縱波速度調節參數；

按照如下公式，根據所述井震誤差和所述實際地層厚度，計算得到縱波速度調節參數的校正值：

其中，δ_i表示縱波速度調節參數的校正值；i表示反演迭代的重復次數；Δh表示井震誤差；h₀表示實際地層厚度；

按照如下公式，根據所述縱波速度調節參數和所述縱波速度調節參數的校正值，累加求和得到新的縱波速度調節參數：

其中，δ表示新的縱波速度調節參數；n表示反演迭代的重復總次數；δ₀表示初始縱波速度調節參數。

2.如權利要求1所述的方法，其特征在于，根據所述初始地震地層厚度、實際地層厚度和巖石物理允許范圍，得到縱波速度調節參數，包括：

根據所述初始地震地層厚度和實際地層厚度，求取初始縱波速度調節參數；

根據所述巖石物理允許范圍，調整初始縱波速度調節參數，得到縱波速度調節參數。

3.如權利要求2所述的方法，其特征在于，根據所述初始地震地層厚度和實際地層厚度，求取初始縱波速度調節參數，包括：

根據獲取的測井資料，計算得到實際地層厚度；

根據所述實際地層厚度和所述初始地震地層厚度，求取初始縱波速度調節參數。

4.如權利要求3所述的方法，其特征在于，按照如下公式，根據所述實際地層厚度和所述初始地震地層厚度，求取初始縱波速度調節參數：

其中，δ₀表示初始縱波速度調節參數；h表示初始地震地層厚度；h₀表示實際地層厚度。

5.如權利要求2所述的方法，其特征在于，根據巖石物理允許范圍，調整初始縱波速度調節參數，得到縱波速度調節參數，包括：

判斷所述初始縱波速度調節參數是否超出巖石物理允許范圍，若超出巖石物理允許范圍，調整初始縱波速度調節參數，得到縱波速度調節參數；

若未超出巖石物理允許范圍，將所述初始縱波速度調節參數作為縱波速度調節參數。

6.如權利要求1所述的方法，其特征在于，按照如下公式，根據所述縱波速度調節參數和所述各向同性縱波速度，得到初始各向異性縱波速度：

其中，v_p表示初始各向異性縱波速度；v_p0表示各向同性縱波速度；δ表示縱波速度調節參數。