1.一種強非均質性儲層中產水氣井的定產量遞減預測方法,其特征在于,該方法包括以下步驟:
S100、收集目標產水氣井的原始地層壓力pi、井口輸壓pt、點測靜壓數據pj與點測靜壓相對應的累產氣量Gpj、地層溫度Ti、井口溫度t、井筒產層中部深度h、井筒半徑rw、無阻流量qAOF、目前的累產氣量Gp、累產水量Wp、日產氣量qg、日產水量qw、氣樣的相對密度γg、氮氣摩爾分數yN2、二氧化碳摩爾分數yCO2、硫化氫摩爾分數yH2S,水樣的相對密度γw和氯化鈉摩爾分數yNaCl;
S200、基于每天的累產水量和每天的累產氣量,得出水驅常數a、水驅常數b,并得到目標產水氣井的甲型水驅公式;
S300、采用Blasingame圖版法擬合得出目標產水氣井的動態儲量G,由點測靜壓相對應的累產氣量除以目標產水氣井的動態儲量得出點測靜壓相對應的儲量采出程度Rj;
S400、收集目標產水氣井壓力恢復試井數據,進行壓力恢復試井分析,計算得出目標產水氣井所在儲層的非均質系數D,具體步驟為:第一、基于目標產水氣井壓力恢復試井得到的壓力隨試井的變化數據,采用雙重介質模型進行數據擬合,得到彈性儲容比ω與竄流系數λ;第二、根據擬合得到的彈性儲容比ω與竄流系數λ,代入計算儲層非均質系數D,其中,α為形狀因子,由目標產水氣井所在儲層中的巖心取心獲得,單位為m-2;rw為井筒半徑,單位為m;λ為竄流系數,無單位;ω為彈性儲容比,無單位;D為儲層非均質系數,無單位;
S500、根據收集到的氣體相對密度γg、原始地層壓力pi、點測靜壓數據p,采用D-A-K法求得原始地層壓力下的偏差因子zi、點測靜壓下的偏差因子z;
S600、結合水封氣物質平衡公式,采用牛頓迭代法計算得出水侵常數C,其中,p為點測靜壓數據,單位為MPa;z為點測靜壓下的偏差因子,無單位;pi為原始地層壓力,單位為MPa;zi為原始地層壓力下的偏差因子,無單位;D為儲層非均質系數,無單位;R為儲量的采出程度,無單位;C為水侵常數,無單位,具體步驟為:第一、由水封氣物質平衡公式得到水侵常數C為未知量的公式,其中,f(C)為表征水侵常數C的公式,無單位;第二、根據f(C)對水侵常數C進行求導,得到,其中f′(C)為f(C)對水侵常數C求導后的公式,無單位;第三、設定水侵常數C為1,帶入f(C)與f′(C),通過C減去f(C)與f′(C)的比值計算得出新的水侵常數C1;第四、計算C與C1的絕對差,若C與C1的絕對差小于0.00001,則C1為所求目標產水氣井的水侵常數;若C與C1的絕對差大于0.00001,則將C替換為C1,重新帶入f(C)與f′(C),計算得出全新的水侵常數C1,直至C與C1的絕對差小于0.00001,最終得出目標氣井的水侵常數C;
S700、對目標產水氣井進行定產量的遞減預測,得出目標產水氣井定產量的穩產年限,具體步驟為:第一、采用Hagedom-Brown法,將原始地層壓力pi、井口輸壓pt、地層溫度Ti、井口溫度t、井筒產層中部深度h、井筒半徑rw、日產氣量qg、日產水量qw、氣樣的相對密度γg、氮氣摩爾分數yN2、二氧化碳摩爾分數yCO2、硫化氫摩爾分數yH2S,水樣的相對密度γw和氯化鈉摩爾分數yNaCl帶入,求得井口輸壓下的井底流壓pwfmin,即穩產末期的井底流壓pwfmin;第二、根據一點法公式計算穩產末期的井底流壓pwfmin下的穩產末期地層壓力pmin;第三、由目標產水氣井目前的累產氣量除以目標產水氣井的動態儲量得出儲量的采出程度R,結合水封氣物質平衡公式與D-A-K法得出目前地層壓力p與目前地層壓力對應的壓縮因子z;第四、目標產水氣井以qg定量生產,以時間1天為迭代步長,疊加Gp得到新一天的累產氣量,代入目標產水氣井的甲型水驅公式計算新一天的累產水量,結合水封氣物質平衡公式與D-A-K法得到新一天的地層壓力,直到新一天的地層壓力小于等于穩產末期地層壓力pmin結束迭代,代入一點法公式反算井底流壓,繪制井底流壓隨時間的變化曲線,得到目標產水氣井定產量遞減預測曲線;第五、根據迭代結束時的時間除以365天得出目標產水氣井的穩產年限。