油藏数值模拟作业….docx

任课教师：刘月田 中国石油大学(北京) 研究生《油藏数值模拟》课程上机作业一 油藏的储层参数分布附图所示，W1井网格位置为(4，8)，W2井位置为(8，5)。 任务：独立编制计算机程序完成第五章渗流问题的模拟求解。具体要求： (一) 使用显式和隐式两种方法。 (二) 输出油藏投产后如下时刻的压力分布、w1井底流压和w2井产油量：初始时刻、10天、1个月、2个月、1个季度、半年、1年、2年、5年；给出压力刚好达到稳定的时刻及其压力分布、w1井底流压和w2井产油量。 (三) 对照(二)中内容，比较分析显式和隐式两种方法的计算过程及结果有何不同。 (四) 建议计算过程用国际单位制，输出的结果中压力用MPa，产量用m3/d。 (五) 作业完成形式。要求三个电子文档： 1. 综合结果报告，word文档； 2．显式求解方法源程序； 3．隐式求解方法源程序。 259222200190180185310240235228210195330290270250230205350300280259222200190180185340320290310240235228210195355335315310290270250230205325300280240210215340320290260235225355335315295275255 附图 油藏的储层参数分布 （1）厚度分布= 附图中数据÷50，单位: m （2）渗透率分布= 附图中数据，单位: （3）孔隙度分布= (附图中数据×0.02+15)%，无单位。 综合结果报告 一、基础数据 Pini=20MPa; μ=5×10-3Pa?s; C=2×10-4/MPa; W1: Q1=30m3/d ; K1=228×10-3 μm 2 ; H1 =4.56m; ? 1=0.1956 W2: Pwf2=15MPa; K2=210×10-3 μm 2 ; H2 =4.2m; ? 2=0.192 ?x=?y=200m; re=0.208?x=0.208×200m=41.6m; rw=0.1m; 显示求解：?t=1.7h; 隐示求解：?t=24h; 稳定条件：下一时间步与上一时间步压力变化绝对值的最大值小于10-6MPa。 二、显式求解方法计算结果 1、刚好达到稳定时： 程序：[P,d,Pwf1,Q2]=explict(2000)  = 1 \* GB3 ①刚好达到稳定的时刻为：36.9天。  = 2 \* GB3 ②此时压力分布：单位：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.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.0000③W1井底流压： Pwf136.9d=16.8500(MPa)  = 4 \* GB3 ④W2井产油量： Q236.9d=44.0837(m3/d) 2、初始时刻： 程序：[P,d,Pwf1,Q2]=explict(0)  = 1 \* GB3 ①压力分布：单位：MPa。 其中:W1井：P10=20MPa; W2井：P20=20MPa。 202020202020202020202020202020