[發明專利]一種油氣藏型地下儲氣庫多周期注采動態庫容計算方法

有效

• 申請號: 202210580896.X
• 申請日: 2022-05-25
• 公開（公告）號: CN114970153B
• 公開（公告）日: 2023-07-25
• 發明（設計）人: 楊棽垚;戚志林;李繼強;嚴文德;黃小亮;袁迎中
• 申請（專利權）人: 重慶科技學院
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06F119/14
• 代理公司: 重慶蘊博君晟知識產權代理事務所(普通合伙) 50223
• 代理人: 王玉芝

1.一種油氣藏型地下儲氣庫多周期注采動態庫容計算方法，其特征在于，包括以下步驟：

獲取不同物性實驗巖樣的巖石常壓孔隙度、不同實驗凈應力和所述不同實驗凈應力對應的覆壓孔隙度；

基于所述不同物性實驗巖樣的巖石常壓孔隙度、所述不同實驗凈應力和所述不同實驗凈應力對應的覆壓孔隙度建立所述實驗巖樣的不可逆孔隙度損失率模型，包括，基于所述巖石不同凈應力及其對應的覆壓孔隙度擬合覆壓孔隙度比值與無因次凈應力關系，得到第一擬合參數；

基于所述巖石覆壓孔隙度擬合孔隙度應力敏感指數與常壓孔隙度關系，得到第二擬合參數；

基于所述第一擬合參數和所述第二擬合參數建立所述實驗巖樣的不可逆孔隙度損失率模型；

基于所述不同物性實驗巖樣的巖石常壓孔隙度和不同實驗巖樣的初始覆壓孔隙度建立巖石初始覆壓孔隙度計算模型，包括，建立巖石初始覆壓孔隙度與巖石常壓孔隙度之間的轉換關系式，根據巖石初始覆壓孔隙度與巖石常壓孔隙度轉換關系式，通過常壓孔隙度計算得到巖石的初始覆壓孔隙度模型；

基于所述不可逆孔隙度損失率模型和所述巖石初始覆壓孔隙度計算模型建立初始凈應力下的巖石覆壓孔隙度計算模型，計算模型為：

式中：為i號實驗巖樣在第j個儲氣庫運行周期內初始凈應力下的覆壓孔隙度；

基于所述初始凈應力下的巖石覆壓孔隙度計算模型建立儲氣庫注采動態庫容計算模型，包括，整理不可逆孔隙度損失率計算模型、巖石初始覆壓孔隙度與巖石常壓孔隙度轉換關系以及初始凈應力下的覆壓孔隙度計算模型，根據巖石室內實驗常壓孔隙度，計算得到任意儲氣庫運行周期初始凈應力下的覆壓孔隙度，考慮孔隙度應力敏感性和孔隙度交變應力敏感性的儲氣庫多周期注采動態庫容評價模型為：

式中：G為儲氣庫庫容，m3；A_g為儲氣庫儲層面積，m2；h為儲氣庫儲層厚度，m；S_wc為束縛水飽和度；B_gi為天然氣體積系數；

將儲氣單元平均參數輸入所述儲氣庫注采動態庫容計算模型得到研究目標儲氣庫的儲層孔隙度和儲氣庫運行周期的儲氣庫庫容。

2.如權利要求1所述的油氣藏型地下儲氣庫多周期注采動態庫容計算方法，其特征在于，

所述獲取不同物性實驗巖樣的巖石常壓孔隙度、所述不同實驗凈應力和所述不同實驗凈應力對應的覆壓孔隙度的具體方式為：

對實驗巖樣開展常壓孔隙度測試實驗，得到常壓孔隙度；

對所述實驗巖樣開展多周期循環載荷覆壓孔隙體積實驗，得到覆壓孔隙體積；

基于所述覆壓孔隙體積計算所述實驗巖樣在不同凈應力狀態下對應的巖石常壓覆壓孔隙度。

3.如權利要求1所述的油氣藏型地下儲氣庫多周期注采動態庫容計算方法，其特征在于，

所述基于所述巖石覆壓孔隙度擬合覆壓孔隙度比值與無因次凈應力關系，得到第一擬合參數的具體方式為：

將所述巖石覆壓孔隙度除以初始覆壓孔隙度，得到覆壓孔隙度比值；

將凈應力除以初始凈應力得到無因次凈應力；

基于所述覆壓孔隙度比值和所述無因次凈應力，得到第一擬合參數。

4.如權利要求3所述的油氣藏型地下儲氣庫多周期注采動態庫容計算方法，其特征在于，

所述基于所述巖石常壓覆壓孔隙度擬合孔隙度應力敏感指數與常壓孔隙度關系，得到第二擬合參數的具體方式為：

采用冪函數回歸分析建立所述凈應力增大過程應力敏感指數與常壓孔隙度的關系模型，得到凈應力增大過程應力敏感指數關系模型；

將所述凈應力減小過程的系數進行算術平均處理，建立凈應力減小過程應力敏感指數的關系模型，得到凈應力減小過程應力敏感指數關系模型；

基于所述凈應力增大過程應力敏感指數關系模型和所述凈應力減小過程應力敏感指數關系模型得到第二擬合參數。