[發明專利]一種計算干熱巖資源潛力區生熱率的方法和系統

在審

• 申請號: 202210360711.4
• 申請日: 2022-04-07
• 公開（公告）號: CN114690270A
• 公開（公告）日: 2022-07-01
• 發明（設計）人: 郭淑君;連晟;程正璞;孫晟;任政委;謝興隆;田蒲源;王春輝
• 申請（專利權）人: 中國地質調查局水文地質環境地質調查中心
• 主分類號: G01V11/00
• 分類號: G01V11/00
• 代理公司: 北京高沃律師事務所 11569
• 代理人: 趙興華

1.一種計算干熱巖資源潛力區生熱率的方法，所述干熱巖資源潛力區中包括已完成鉆探的若干個鉆井，其特征在于，所述方法包括以下步驟：

采集干熱巖資源潛力區的多組巖心，并對多組巖心分別進行巖樣制備，得到多組測試巖樣；所述多組巖心為分別從干熱巖資源潛力區每個鉆井中取出的指定地層時代的巖心；所有鉆井的指定地層時代均相同；所述指定地層時代為所述干熱資源潛力區具有代表性的地層；

獲取所述測試巖樣在實驗室環境下的縱波波速實驗室測量值和生熱率實驗室測量值；所述實驗室環境包括溫度和壓力；

確定第一導數關系式、第二導數關系式以及縱波波速-生熱率對數經驗關系式；所述第一導數關系式為溫度對所述縱波波速實驗室測量值的導數關系式；所述第二導數關系式為壓力對所述縱波波速實驗室測量值的導數關系式；所述縱波波速-生熱率對數經驗關系式為所述縱波波速實驗室測量值和所述生熱率實驗室測量值之間的對數關系式；

對干熱巖資源潛力區的每兩個鉆井之間進行層析成像反演，得到每兩個鉆井之間的縱波波速原位測量剖面；所述縱波波速原位測量剖面包括若干個縱波波速原位測量值；

根據所述第一導數關系式以及所述第二導數關系式，確定校正因子；

根據所述校正因子對所述縱波波速原位測量值進行校正，得到實驗室環境下的縱波波速值；

根據所述實驗室環境下的縱波波速值以及所述縱波波速-生熱率對數經驗關系式確定所述縱波波速原位測量值對應的生熱率測量值，得到每兩個鉆井之間的生熱率剖面；

根據各個生熱率剖面進行插值細化，得到所述干熱巖資源潛力區各測量點的生熱率。

2.根據權利要求1所述的方法，其特征在于，所述對干熱巖資源潛力區的每兩個鉆井之間進行層析成像反演，得到每兩個鉆井之間的縱波波速原位測量剖面，具體包括：

利用地震層析成像單元按照層級間距逐層采集每兩個鉆井之間的地震數據；每兩個鉆井中任意一個鉆井為激發井，另一個鉆井為接收井；

對逐層采集得到的所述地震數據進行井間地震層析成像反演，依次得到兩個鉆井之間每層的縱波波速原位測量值，兩個鉆井之間的所有縱波波速原位測量值構成兩個鉆井之間的縱波波速原位測量剖面。

3.根據權利要求2所述的方法，其特征在于，利用射線層析法和波動方程層析法對逐層采集得到的所述地震數據進行井間地震層析成像反演。

4.根據權利要求2所述的方法，其特征在于，在對逐層采集得到的所述地震數據進行井間地震層析成像反演，依次得到兩個鉆井之間每層的縱波波速原位測量值之后，所述方法還包括：

利用鉆井地層資料對得到的每層的縱波波速原位測量值進行校驗；所述鉆井地層資料為所述鉆井在進行鉆探時所記錄的地層數據，包括所述鉆井各地層的縱波波速鉆孔測量值；

若所述縱波波速原位測量值與所述縱波波速鉆孔測量值之差大于設定閾值，則修改所述層級間距，并重新利用地震層析成像單元按照層級間距逐層采集每兩個鉆井之間的地震數據。

5.根據權利要求1所述的方法，其特征在于，所述根據所述第一導數關系式以及所述第二導數關系式，確定校正因子，具體包括：

確定縱波波速原位測量值處的溫度變化量和壓力變化量；所述溫度變化量為溫度隨深度的變化量，所述壓力變化量為壓力隨深度的變化量；

根據所述溫度變化量、所述壓力變化量、所述第一導數關系式以及所述第二導數關系式，計算得到校正因子。

6.根據權利要求5所述的方法，其特征在于，根據下式計算得到校正因子：

其中，ΔT表示所述溫度變化量，ΔP表示所述壓力變化量，v_p室表示所述縱波波速實驗室測量值，T_室和P_室分別表示實驗室環境下的溫度和壓力。

7.根據權利要求6所述的方法，其特征在于，根據下式將所述縱波波速原位測量值校正為實驗室環境下的縱波波速值：

v_p室＝v_p原·f

其中，v_p原表示所述縱波波速原位測量值。