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

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

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

3、0-/znr（3）孔隙度分布=（附圖中數(shù)據(jù)X002+15）%,無單位。任課教師：劉刀田綜合結(jié)果報告Pint = 20MPa;p = 5 X 1()7 Pas； C = 2X 10-/MPa；1： Qi = 30m3/d；Ki = 228 X 103 “m2； 弘W2： Pwf2 = ISMPa；K2 =210x10-3 “m2.=456m; 01= 0,1956H2=4.2m; 02 = 0.192Ax = Ay = 200m;Vg = 0208Ax = 0.208 X 200rn = 41.6m； 心=O.lzn;顯示求解：At = 1.7h；隱示求解：At = 24h;穩(wěn)定條件：下一時間步

4、與上一時間步壓力變化絕對值的最大值小J lO-MPao二、顯式求解方法計算結(jié)果1、剛好達到穩(wěn)定時:程序：P,aPwfl,Q2=explict(2000) 剛好達到穩(wěn)定的時刻為：36.9天。 此時壓力分布；單位：MPa。其中:W1 井：P嚴9d = i84527MPa; W?井：P茫刃=177762MPa。19219118918.718.718.763569469313253594119319.218.9118418.7312349885w42316W1井底流壓:P詳d = 168500(MPa)W2井產(chǎn)油最:02% 加=44,0837 (mVci)2. 初始時刻:程岸:P,d,Pwfl,Q2=

5、explict(0) 壓力分布：單位；MPao其中:W1 井：Pf = 20MPa; 井：Ff = 20MPa.202020202020202020：02020任課教師：劉/田202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020202020 W1井底流壓:= 18.3973 (MPa) W2井產(chǎn)汕最:Q2 = 793954(m3/d)3、10天時:程岸：P,d,Pwfl,Q2j=explict(10) 壓力分布：單位；MPJU其中:Wi井J Pfd = i

6、85265MPa; W?井：P/% = 177938MPa。19319.219J)18.7181837645448295539382219319319D1818997802959031311119占19419.21918J18.2109725785484040051573094361 W1井底流壓:Pi = 168503 (MPa) W2井產(chǎn)汕量:QzOd =440846(m3/d)5、穩(wěn)定后由于顯式法求解達到穩(wěn)定的時間是36.9天，穩(wěn)定后床力和產(chǎn)最不再隨時間變化，所以2個月.一個季度、半年.一年、2年、5年的壓力分 布.W1井底流壓和W2井產(chǎn)油最與36.9天時相同。6. 結(jié)果小結(jié)1顯式方法計

7、算結(jié)果初始 時刻10天1個月穩(wěn)定時間 （36.9 天）2個月、1個季度、半年、1年、2年、5年P(guān)l (MPa)201& 526518. 45318. 452718. 4527?2 (MPa)2017. 793817. 776317. 776217. 7762井底流壓Pwfl (MPa)1& 397316. 92381& 85031& 8516. 85W,井產(chǎn)量Q2 (m7d)79. 395444. 362444. 084644. 083744.0837三.隱式求解方法計算結(jié)果K剛好達到穩(wěn)定時: 程序:P,d,Pwfl ,Q2=iiiiplict(2000) 剛好達到穩(wěn)定的時刻為：53天。 此時

8、壓力分布：單位：MPa。其中:Wi井：P嚴d = 184527MPa; W?井：= r7.7762MPa.1919418918718-718-763569369313153594019319F1891184ET18.718$31：34988454131519419319.219659746790283742185任課教師：劉刀田19419319J)18518$189232291867063044119519419319419J0191126764118112780J701919A19.519419.21918.71&218.26432653540926946910685027611931919

9、.71919419J18.718.226208351132063396144093820112OD20 X)19$19j619-419418-718$00000000027152730931384305720Q19319爲19.419319-200018860889926488620Q1951919-719-719600022850082223890720Q201120 J)20B20 J)2OJ0000000000000000000 W1井底流壓:=16.9535 (MPa) W2井產(chǎn)油最:Q2IM = 44.5265 (m3/ci)4、1個月(約為第30天)：程岸：P,d,Pwfl,Q2=

10、iinplict(30) ffi力分布：單位：MPju其中:W1 井：PfM = 18.4535MPa； W?井：= 17.7764MPao19J19418918718.718.764370270214154495019-319.218911841851831935689255032419419319.21901891895218823862838874531919519.419319.219D11172778847043OM18.218.209636218必57519J98819780819721119595919422319.153918.7066|l75制18.260820D20X)20

11、4)19319A19.419.118.7187000000000125300or03457978420X)19$1919.419319.200008445272007969720X)1991919/719-719300018743274115281620Q20 JI201120Q20D20A000000000000000000 W1井底流壓:Pi = 168508(MPa) W2井產(chǎn)汕最:=44.0865 (m3/d)由于隱式求解法達到穩(wěn)定的時間是53天，穩(wěn)定后床力和產(chǎn)量不再隨時間變化，所以2個刀、一個季度、半年、一年、2年.5年的壓力分布、W1井底流壓和W2井產(chǎn)油最與53天時相同。6、結(jié)果小

12、結(jié)表2Rft氏方法il算結(jié)果初始 時刻10天1個月穩(wěn)定時間 （53 天）2個月、1個季度、半年、1年、2年、5年P(guān)l (MPa)2018. 556218. 45351& 452718. 4527P2 (MPa)2017. 804117. 776417. 776217. 7762W1井厳流用Pwfl (MPa)18. 397316. 953516. 850816. 8516. 85W2井產(chǎn)Q2 (m7d)79. 395444. 526544. 08654-1. 083544. 0835表3顯式、隱式方法求解結(jié)果對比結(jié)果時間、方法P1 (MPa)P： (MPa)乳井底流壓Pwfl （MPa）W?井

13、產(chǎn)量Q2 (ra7d)初始 時刻顯式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)定 時間顯式（36. 9 天）18.452717. 776216. 8544. 0837隱式 （53 天）18. 452717. 776216. 8544. 0835穩(wěn)定 后顯式18.452717

14、. 776216. 8544. 0837隱式18.452717. 776216. 8544. 0835從表3中可以看出: 兩種方法穩(wěn)定前，同一時刻，1）隱式方法計算的地層壓力P1 P2比顯式方法計算的地層壓力P1、P2髙；2）相應(yīng)的隱式方法計算的W1 井底流壓和W2井產(chǎn)量也比顯式方法計算的高一些; 顯式方法計算的床力達到穩(wěn)定的時間為36.9天，隱式為53天: 兩種方法都達到穩(wěn)定后-1）隱式一部分地層壓力分布比顯式低一些（由于只保留四位有效數(shù)字，看起來不是很明顯）；2）相應(yīng)的穩(wěn)定產(chǎn)最 隱式較顯式稍低； 以上結(jié)論可以說明，在達到同樣的精度時，隱式方法計算結(jié)果更為準確。2、計算過程對比 顯式差分方法:假設(shè)第n時間步的值己知，則方程中只有一個未知最P/j+i,可以直接計算得到。在t=0的初值給定后，便可利用顯式差分方程逐點逐層對整個 問題進行計算。求解步驟比較簡單，但是穩(wěn)定性較差，本問題中