[發明專利]一種基于多元回歸的氣井早期采收率計算方法

有效

• 申請號: 202210019135.7
• 申請日: 2022-01-10
• 公開（公告）號: CN114547850B
• 公開（公告）日: 2023-03-28
• 發明（設計）人: 魏明強;杜怡鶴;段永剛;李舫;杜凌云;林興;李政瀾;任科屹;龍騰屹;方全堂
• 申請（專利權）人: 西南石油大學
• 主分類號: G06F30/20
• 分類號: G06F30/20;G06F119/08;G06F119/14
• 代理公司: 深圳峰誠志合知識產權代理有限公司 44525
• 代理人: 胡石開

1.一種基于多元回歸的氣井早期采收率計算方法，其特征在于，包括以下步驟：

S1、收集已開發區塊氣井的地層參數與生產數據；地層參數包括廢棄壓力、氣體溫度和井內氣體特性參數；井內氣體特性參數包括氣體粘度和偏差因子；

S2、基于地層參數建立氣井早期采收率計算模型，進而得到不同井型早期采收率計算模型的全微分形式及其邊界條件；

S3、采用多元回歸算法離散處理含氣區半徑的變化過程，得到不同井型井控半徑及其邊界條件；

S4、以氣井停產時的廢棄壓力對應的井口流量、溫度、粘度、偏差因子為經濟極限最低狀態參數，結合不同井型井控半徑、井控半徑邊界條件、氣體黏度邊界條件和廢棄壓力邊界條件，通過對步驟S2得到的不同井型早期采收率計算模型的全微分形式求解，得到氣井早期采收率解析解直接形式的計算方程；

S5、利用氣井早期采收率解析解直接形式的計算方程計算未開采氣井的早期采收率；

所述氣井早期采收率計算模型如下：

其中，

式中：R_wf為氣井早期采收率，％；q為氣井井口流量，m³/d；Δt為氣井生產時間段，d；N為氣藏地質儲量，m³；F()為氣井早期采收率計算函數；P_S為氣井停產時的地層壓力，MPa；T為氣體溫度，K；μ為地層氣體粘度，mPa·s；Z為氣體偏差因子，無量綱；K為含氣層巖石滲透率，D；h為含氣層厚度，m；C為標準常數；P為含氣區半徑r處的地層壓力，MPa；P_wf為井底流壓，MPa；r為含氣區半徑，m；r_w為井筒半徑，m；

所述不同井型早期采收率計算模型的全微分形式及其邊界條件如下：

式中，q_i為t_i時刻氣井井口流量，m³/d；Δt_i＝t_i-t_i-1,為生產時間段；r_i為時間點t_i所對應的含氣區半徑，m；ρ為氣體密度，kg/m³；R為氣體常數；M為氣體分子質量，g/mol；μ_s為氣井停產時所采出的氣體粘度，mPa·s；

步驟S3中所述不同井型井控半徑計算式如下：

井控半徑的邊界條件如下：

其中，

式中，F_i為生產時間段Δt_i內氣井累計采出氣量所占氣井總儲氣量的百分值，％；r_e為井控半徑，m；Ω為概率計算函數，％；Y代表不同變量；ω_i為自變量權重，指所測量的氣井生產時間段Δt_i占總生產時間的百分比，無量綱；為平均估計值；b為概率計算常數，無量綱；/表示概率對數輔助函數；n表示測量次數。

2.根據權利要求1所述的一種基于多元回歸的氣井早期采收率計算方法，其特征在于，步驟S4中所述氣井早期采收率計算模型的解析解直接形式如下：

式中：Z_S為氣井停產時所采出的氣體偏差因子，無量綱；P_P為原始地層壓力，MPa；Z_P為初始時刻t＝0時的氣體偏差因子，無量綱；α為計算系數；q_s為氣井停產時井口流量，m³/d；

T_S為氣井停產時所采出的氣體溫度，K；μ_s為氣井停產時所采出的氣體粘度，mPa·s；P_wfs為氣井停產時井底流壓，MPa。