[發明專利]一種用于預測頁巖氣資源的方法

有效

• 申請號: 201810289872.2
• 申請日: 2018-03-30
• 公開（公告）號: CN110318744B
• 公開（公告）日: 2022-01-21
• 發明（設計）人: 金之鈞;盛秀杰;李治昊;王鵬威;吳博
• 申請（專利權）人: 中國石油化工股份有限公司;中國石油化工股份有限公司石油勘探開發研究院
• 主分類號: E21B49/00
• 分類號: E21B49/00
• 代理公司: 北京聿宏知識產權代理有限公司 11372
• 代理人: 吳大建;何嬌

1.一種用于預測頁巖氣資源的方法，包括：

步驟S100：對待評價區域內的多個處進行測井以獲得每個井位處目的儲層的有機碳含量、烴源巖熱成熟度指標、頁巖厚度、頁巖密度、頁巖含氣量、埋深；

對所有頁巖密度求平均以獲得平均頁巖密度；

步驟S200：對待評價區域進行均勻的網格劃分，待評價區域內設置有井位的單元格為已知單元格，其他為未知單元格；

步驟S300：將每個已知單元格中的一個井位處的有機碳含量、烴源巖熱成熟度指標、頁巖厚度、頁巖含氣量、埋深的值依次分別作為目的儲層在該已知單元格范圍內有機碳含量、烴源巖熱成熟度指標、頁巖厚度、頁巖含氣量、埋深的均值；

步驟S400：隨機構建模擬所有未知單元格的地質屬性的順序，按順序依次對每個未知單元格的地質屬性進行模擬，

模擬第一個未知單元格的地質屬性則進行步驟S410，

模擬第二個未知單元格至最后一個未知單元格的地質屬性時，每對一個未知單元格的地質屬性進行模擬則進行步驟S420；

步驟S410：根據每個已知單元格的有機碳含量的均值、烴源巖熱成熟度指標的均值、頁巖厚度的均值、頁巖含氣量的均值、埋深的均值、各個單元格的位置和當前待計算的未知單元格的位置采用克里金插值法計算出該未知單元格中有機碳含量插值、烴源巖熱成熟度指標插值、頁巖厚度插值、頁巖含氣量插值、埋深插值；

步驟S420：根據每個已知單元格的有機碳含量的均值、烴源巖熱成熟度指標的均值、頁巖厚度的均值、頁巖含氣量的均值、埋深的均值，已經計算出的未知單元格中有機碳含量插值、烴源巖熱成熟度指標插值、頁巖厚度插值、頁巖含氣量插值、埋深插值，以及各個單元格的位置和當前待計算的未知單元格的位置采用克里金插值法計算出該未知單元格中有機碳含量插值、烴源巖熱成熟度指標插值、頁巖厚度插值、頁巖含氣量插值、埋深插值；

步驟S500：對每個已知單元格的有機碳含量的均值進行評價獲得有機碳評分值，對每個未知單元格的有機碳含量插值進行評價獲得有機碳評分值，對每個已知單元格的烴源巖熱成熟度指標的均值進行評價獲得熱成熟度評分值，對每個未知單元格的烴源巖熱成熟度指標插值進行評價獲得熱成熟度評分值，對每個已知單元格的頁巖厚度的均值進行評價獲得厚度評分值，對每個未知單元格的頁巖厚度插值進行評價獲得厚度評分值，對每個已知單元格的埋深的均值進行評價獲得經濟可采性評分，對每個未知單元格的埋深插值進行評價獲得經濟可采性評分；

步驟S600：將每個已知單元格的有機碳評分值、熱成熟度評分值、厚度評分值相乘獲得每個已知單元格的頁巖氣富集概率；

將每個未知單元格的有機碳評分值、熱成熟度評分值、厚度評分值相乘獲得每個未知單元格的頁巖氣富集概率；

步驟S700：將有機碳評分值、熱成熟度評分值、厚度評分值、經濟可采性評分均大于或等于0.5分的單元格的集合劃為有利區；

步驟S800：計算出有利區內每個已知單元格對應區域中的天然氣地質資源量和有利區內每個未知單元格對應區域中的天然氣地質資源量，

將有利區內每個已知單元格對應區域中的天然氣地質資源量與有利區內每個未知單元格對應區域中的天然氣地質資源量進行累加獲得累積天然氣地質資源量；

步驟S900：重復步驟S400至步驟S800預設次數，以獲得預設個數的累積天然氣地質資源量；

步驟S1000：以累積天然氣地質資源量數值從大到小的順序對所有累積天然氣地質資源量從序號1開始進行連續編號，按照算式(1)計算每個序號對應的超概率，

超概率＝(序號/N+(序號-1)/N)/2 (1)

其中，N是累積天然氣地質資源量的個數，

將每個累積天然氣地質資源量及其對應的超概率以坐標點的形式標記在以橫軸表示累積天然氣地質資源量數值、縱軸表示超概率的直角坐標系上，將直角坐標系上相鄰的兩個點用直線段連接起來得到超概率曲線，將超概率曲線上超概率等于百分之五十所對應的累積天然氣地質資源量作為待評價區域的天然氣地質資源量的最終估值。