油藏數(shù)值模擬作業(yè)

1、中國石油大學(xué)(北京)研究生油藏數(shù)值模擬課程上機作業(yè)一油藏的儲層參數(shù)分布附圖所示，W1井網(wǎng)格位置為(4，8)，W2井位置為(8，5)。任務(wù)：獨立編制計算機程序完成第五章滲流問題的模擬求解。具體要求： (一) 使用顯式和隱式兩種方法。(二) 輸出油藏投產(chǎn)后如下時刻的壓力分布、w1井底流壓和w2井產(chǎn)油量：初始時刻、10天、1個月、2個月、1個季度、半年、1年、2年、5年；給出壓力剛好達到穩(wěn)定的時刻及其壓力分布、w1井底流壓和w2井產(chǎn)油量。(三) 對照(二)中內(nèi)容，比較分析顯式和隱式兩種方法的計算過程及結(jié)果有何不同。(四) 建議計算過程用國際單位制，輸出的結(jié)果中壓力用MPa，產(chǎn)量用m3/d。(五) 作

2、業(yè)完成形式。要求三個電子文檔：1. 綜合結(jié)果報告，word文檔；2顯式求解方法源程序；3隱式求解方法源程序。259222200190180185310240235228210195330290270250230205350300280259222200190180185340320290310240235228210195355335315310290270250230205325300280240210215340320290260235225355335315295275255附圖 油藏的儲層參數(shù)分布（1）厚度分布= 附圖中數(shù)據(jù)50，單位: m（2）滲透率分布= 附圖中數(shù)據(jù)，單位: （3）

3、孔隙度分布= (附圖中數(shù)據(jù)0.02+15)%，無單位。綜合結(jié)果報告一、基礎(chǔ)數(shù)據(jù)Pini=20MPa;=510-3Pas;C=210-4/MPa;W1: Q1=30m3/d ; K1=22810-3 m 2 ; H1 =4.56m; 1=0.1956W2: Pwf2=15MPa; K2=21010-3 m 2 ; H2 =4.2m; 2=0.192x=y=200m;re=0.208x=0.208200m=41.6m;rw=0.1m;顯示求解：t=1.7h;隱示求解：t=24h;穩(wěn)定條件：下一時間步與上一時間步壓力變化絕對值的最大值小于10-6MPa。二、顯式求解方法計算結(jié)果1、剛好達到穩(wěn)定時：程

4、序：P,d,Pwf1,Q2=explict(2000)剛好達到穩(wěn)定的時刻為：36.9天。此時壓力分布：單位：MPa。其中:W1井：P136.9d=18.4527MPa; W2井：P236.9d=17.7762MPa。19.263519.169418.969318.713218.753518.794119.331219.234918.988518.452718.754218.831619.451619.387619.238019.027718.988018.944619.610719.572319.478319.348319.203819.005018.657218.209318.236019.

5、798519.780519.720919.595719.422019.153618.706317.776218.260620.000020.000020.000019.812419.629919.403519.103218.757718.778320.000019.808319.645119.471919.307719.269620.000019.918719.843219.774019.715219.681520.000020.000020.000020.000020.000020.0000W1井底流壓：Pwf136.9d=16.8500(MPa)W2井產(chǎn)油量：Q236.9d=44.0837

6、(m3/d)2、初始時刻：程序：P,d,Pwf1,Q2=explict(0)壓力分布：單位：MPa。其中:W1井：P10=20MPa; W2井：P20=20MPa。202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020W1井底流壓：Pwf10=18.3973(MPa)W2井產(chǎn)油量：Q20=79.3954(m3/d)3、10天時：程序：P,d,Pwf1,Q2=explict(10)壓力分布：單位：MP

7、a。其中:W1井：P110d=18.5265MPa; W2井：P210d=17.7938MPa。19.337619.245419.048218.795518.839318.882219.397819.302919.059018.526518.831318.911119.505919.442719.295419.088119.051519.009419.649419.611319.519119.391419.249719.049018.689918.235218.260619.818519.800519.742419.621119.450419.182218.730917.793818.2809

8、20.000020.000020.000019.822619.645319.421619.121318.774318.795020.000019.815319.655319.483319.319219.281120.000019.921419.847619.779219.720519.687020.000020.000020.000020.000020.000020.0000W1井底流壓： Pwf110d =16.9238(MPa)W2井產(chǎn)油量：Q210d=44.3624(m3/d)4、1個月（約為第30天）：程序：P,d,Pwf1,Q2=explict(30)壓力分布：單位：MPa。其中:W

9、1井：P130d=18.4530MPa; W2井：P230d=17.7763MPa。19.263819.169618.969618.713518.753818.794419.331519.235218.988718.453018.754418.831919.451719.387819.238219.027918.988218.944819.610919.572519.478519.348419.204019.005118.657318.209418.236119.798619.780619.721019.595819.422119.153718.706417.776318.260720.000

10、020.000020.000019.812419.629919.403619.103218.757818.778320.000019.808419.645119.471919.307819.269620.000019.918719.843219.774119.715219.681620.000020.000020.000020.000020.000020.0000W1井底流壓：Pwf130d=16.8503(MPa)W2井產(chǎn)油量：Q230d=44.0846(m3/d)5、穩(wěn)定后由于顯式法求解達到穩(wěn)定的時間是36.9天，穩(wěn)定后壓力和產(chǎn)量不再隨時間變化，所以2個月、一個季度、半年、一年、2年、5年

11、的壓力分布、W1井底流壓和W2井產(chǎn)油量與36.9天時相同。6、結(jié)果小結(jié)表1顯式方法計算結(jié)果初始10天1個月穩(wěn)定時間（36.9天）2個月、1個季度、半年、時刻1年、2年、5年P(guān)1（MPa）2018.526518.45318.452718.4527P2（MPa）2017.793817.776317.776217.7762W1井底流壓Pwf1（MPa）18.397316.923816.850316.8516.85W2井產(chǎn)量Q2（m3/d）79.395444.362444.084644.083744.0837三、隱式求解方法計算結(jié)果1、剛好達到穩(wěn)定時：程序：P,d,Pwf1,Q2=implict(20

12、00)剛好達到穩(wěn)定的時刻為：53天。此時壓力分布：單位：MPa。其中:W1井：P153d=18.4527MPa; W2井：P253d=17.7762MPa。19.263519.169318.969318.713118.753518.794019.331219.234918.988418.452718.754118.831519.451519.387619.238019.027618.988018.944519.610719.572319.478319.348219.203819.005018.657218.209318.236019.798519.780519.720919.595719.42

13、2019.153518.706317.776218.260620.000020.000020.000019.812419.629919.403519.103218.757718.778320.000019.808319.645119.471819.307719.269620.000019.918719.843119.774019.715219.681520.000020.000020.000020.000020.000020.0000W1井底流壓：Pwf153d=16.8500(MPa)W2井產(chǎn)油量：Q253d=44.0835(m3/d)2、初始時刻：程序：P,d,Pwf1,Q2=implic

14、t(0)壓力分布：單位：MPa。其中:W1井：P10=20MPa; W2井：P20=20MPa。202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020W1井底流壓：Pwf10=18.3973(MPa)W2井產(chǎn)油量：Q20=79.3954(m3/d)3、10天時：程序：P,d,Pwf1,Q2=implict(10)壓力分布：單位：MPa。其中:W1井：P110d=18.5562MPa; W2井：P21

15、0d=17.8041MPa。19.365919.274619.079018.828318.874218.918519.423219.329119.086718.556218.863018.944119.526719.464119.318119.112619.078019.037019.664319.626519.535419.409219.269419.069118.706818.250218.276119.826219.808319.751119.632019.463319.196118.744017.804118.293820.000020.000020.000019.827119.652

16、719.430919.131318.784318.805720.000019.818819.660819.489919.326419.288620.000019.922819.850019.782219.723819.690720.000020.000020.000020.000020.000020.0000W1井底流壓： Pwf110d =16.9535(MPa)W2井產(chǎn)油量：Q210d=44.5265(m3/d)4、1個月（約為第30天）：程序：P,d,Pwf1,Q2=implict(30)壓力分布：單位：MPa。其中:W1井：P130d=18.4535MPa; W2井：P230d=17.

17、7764MPa。19.264319.170218.970218.714118.754418.795019.331919.235618.989218.453518.755018.832419.452119.388219.238619.028318.988718.945319.611119.572719.478819.348719.204319.005418.657518.209618.236219.798819.780819.721119.595919.422319.153918.706617.776418.260820.000020.000020.000019.812519.630019.40

18、3719.103418.757918.778420.000019.808419.645219.472019.307919.269720.000019.918719.843219.774119.715219.681620.000020.000020.000020.000020.000020.0000W1井底流壓：Pwf130d=16.8508(MPa)W2井產(chǎn)油量：Q230d=44.0865(m3/d)5、穩(wěn)定后由于隱式求解法達到穩(wěn)定的時間是53天，穩(wěn)定后壓力和產(chǎn)量不再隨時間變化，所以2個月、一個季度、半年、一年、2年、5年的壓力分布、W1井底流壓和W2井產(chǎn)油量與53天時相同。6、結(jié)果小結(jié)表2隱

19、式方法計算結(jié)果初始10天1個月穩(wěn)定時間（53天）2個月、1個季度、半年、1年、2年、5年時刻P1（MPa）2018.556218.453518.452718.4527P2（MPa）2017.804117.776417.776217.7762W1井底流壓Pwf1（MPa）18.397316.953516.850816.8516.85W2井產(chǎn)量Q2（m3/d）79.395444.526544.086544.083544.0835四、顯隱式方法對比1、結(jié)果對比表3顯式、隱式方法求解結(jié)果對比 結(jié)果時間、方法P1（MPa）P2（MPa）W1井底流壓Pwf1（MPa）W2井產(chǎn)量Q2（m3/d）初始時刻顯式202018.397379.3954隱式202018.397379.395410天顯式18.526517.793816.923844.3624隱式18.556217.804116.953544.52651個月顯式18.45317.776316.850344.0846隱式18.453517.776416.850844.0865穩(wěn)定時間顯式（3