[發明專利]儲層孔隙度預測方法及裝置

在審

• 申請號: 201911256155.0
• 申請日: 2019-12-10
• 公開（公告）號: CN112946754A
• 公開（公告）日: 2021-06-11
• 發明（設計）人: 夏亞良;魏小東;陳鑫;董建雄;安福利;楊建房
• 申請（專利權）人: 中國石油天然氣集團有限公司;中國石油集團東方地球物理勘探有限責任公司
• 主分類號: G01V1/50
• 分類號: G01V1/50
• 代理公司: 北京三友知識產權代理有限公司 11127
• 代理人: 周曉飛;許曼

1.一種儲層孔隙度預測方法，其特征在于，包括：

獲取多個單井在目的層段的儲層孔隙度和儲層厚度，以及目的層段地震反演縱波阻抗數據；

根據所述目的層段地震反演縱波阻抗數據，提取得到目的層段地震反演縱波阻抗分布平面數據，確定所述多個單井在目的層段的縱波阻抗值；

根據所述多個單井中每個單井在目的層段的儲層孔隙度和儲層厚度，確定每個單井在目的層段的儲層孔隙厚度積；

根據每個單井在目的層段的儲層孔隙厚度積和縱波阻抗值，擬合得到目的層段的儲層孔隙厚度積與目的層段的縱波阻抗之間的第一關系；

根據所述第一關系，確定目的層段的儲層孔隙度與目的層段的儲層厚度和目的層段的縱波阻抗之間的第二關系；

根據所述目的層段地震反演縱波阻抗數據，確定目的層段的儲層厚度分布平面數據；

根據所述目的層段地震反演縱波阻抗分布平面數據、所述目的層段的儲層厚度分布平面數據以及所述第二關系，確定目的層段的儲層孔隙度分布平面數據。

2.如權利要求1所述的方法，其特征在于，獲取目的層段地震反演縱波阻抗數據，包括：

獲取目的層段所在區域的地震資料以及測井資料；

根據所述測井資料，計算得到所述多個單井的縱波阻抗的測井曲線；

根據多個單井的縱波阻抗的測井曲線，得到多個單井中每個單井在目的層段的縱波阻抗的測井計算值；

根據所述地震資料以及測井資料，完成目的層段的地震地質標定，開展目的層段地震縱波阻抗反演或地質統計學反演，確定目的層段地震反演縱波阻抗數據；所述目的層段地震反演縱波阻抗數據中的每個單井在目的層段的縱波阻抗的地震反演值與每個單井在目的層段的縱波阻抗的測井計算值一致。

3.如權利要求2所述的方法，其特征在于，根據所述目的層段地震反演縱波阻抗數據，確定目的層段的儲層厚度分布平面數據，包括：

根據所述多個單井的縱波阻抗的測井曲線，確定儲層與非儲層的門檻值；所述門檻值為用于劃分儲層與非儲層的縱波阻抗值；

根據所述門檻值以及所述目的層段地震反演縱波阻抗數據，確定目的層段的儲層厚度分布平面數據。

4.如權利要求3所述的方法，其特征在于，還包括：

根據獲取的多個單井在目的層段的儲層厚度，對所述目的層段的儲層厚度分布平面數據進行校正。

5.如權利要求1所述的方法，其特征在于，按如下公式，根據多個單井中每個單井在目的層段的儲層孔隙度和儲層厚度，確定每個單井在目的層段的儲層孔隙厚度積：

其中，PH_i表示第i口單井在目的層段的儲層孔隙厚度積；POR_ij表示第i口單井在目的層段中第j段儲層的孔隙度；T_ij表示第i口單井在目的層段中第j段儲層的厚度；i表示單井編號；j表示儲層段序號；n表示目的層段中儲層段的總數。

6.如權利要求5所述的方法，其特征在于，按如下公式，擬合得到所述第一關系：

AI＝a×ln(PH)+b

其中，AI表示目的層段的縱波阻抗；PH表示目的層段的儲層孔隙厚度積；a和b表示經驗系數。

7.如權利要求6所述的方法，其特征在于，按如下公式，確定所述第二關系：

其中，POR'表示目的層段的儲層孔隙度；AI表示目的層段的縱波阻抗；T'表示目的層段的儲層厚度；c表示井震轉化系數，為一常數。