油藏的储层参数分布附图所示独立编制计算机程序完成第.docx

油藏的储层参数分布附图所示独立编制计算机程序完成第.docx

《油藏的储层参数分布附图所示独立编制计算机程序完成第.docx》由会员分享，可在线阅读，更多相关《油藏的储层参数分布附图所示独立编制计算机程序完成第.docx（16页珍藏版）》请在冰豆网上搜索。

油藏的储层参数分布附图所示独立编制计算机程序完成第

油藏的储层参数分布附图所示。

独立编制计算机程序完成第五章渗流问题的模拟求解。

具体要求：

（一）使用显式和隐式两种方法。

（二）输出油藏投产后如下时刻的压力分布、w1井底流压和w2井产油量：

初始时刻、10天、1个月、2个月、1个季度、半年、1年、2年、5年；给出压力刚好达到稳定的时刻及其压力分布、w1井底流压和w2井产油量。

（三）对照

（二）中内容，比较分析显式和隐式两种方法的计算过程及结果有何不同。

（四）建议计算过程用国际单位制，输出的结果中压力用MPa，产量用m3/d。

（五）2010年5月15日之前交作业。

（六）作业完成形式。

要求三个电子文档：

1.综合结果报告，word文档；

2．显式求解方法源程序；

3．隐式求解方法源程序。

259

222

200

190

180

185

310

240

235

228

210

195

330

290

270

250

230

205

350

300

280

259

222

200

190

180

185

340

320

290

310

240

235

228

210

195

355

335

315

310

290

270

250

230

205

325

300

280

240

210

215

340

320

290

260

235

225

355

335

315

295

275

255

附图油藏的储层参数分布

（1）厚度分布=附图中数据÷50，单位:

m

（2）渗透率分布=附图中数据，单位:

（3）孔隙度分布=（附图中数据×0.02+15）%，无单位。

（4）其它参数见课上笔记。

一、概述

本次作业在MATLAB上完成，输出结果显式在11.txt文件中，隐式在22.txt文件中。

使用显式和隐式两种方法。

显式方法包括三个函数文件：

clear，outputs，outputa。

其中clear为主函数，实现压力和产量计算；outputs和outputa实现输出。

隐式方法包括四个函数文件：

implicitmain，AG，outputs，outputa。

其中implicit为主函数，实现压力和产量计算；AG函数计算系数矩阵A和右端项G，outputs和outputa实现输出

二、显式运算结果

x方向步长（m）：

200.0000

y方向步长（m）：

200.0000

时间步长（h）：

2.0000

============================================================================================

时刻tt（d）：

0

压力分布（MPa）：

20.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.0000

定产井w1井底流压（MPa）：

18.490414

定压井w2产量（m^3/d）：

112.454746

============================================================================================

时刻tt（d）：

10

压力分布（MPa）：

19.164119.059518.882518.958819.029119.0587

19.232919.086618.643218.961319.056919.0748

19.383919.313219.172219.175219.149719.0980

19.571319.535919.458219.382619.248419.053118.813518.774218.7945

19.780619.764519.707619.576019.356919.018718.627218.715918.8087

20.000020.000020.000019.774219.486519.020718.040418.676618.9127

20.000019.708519.412319.075219.162719.2556

20.000019.872819.745819.641219.641619.6529

20.000020.000020.000020.000020.000020.0000

定产井w1井底流压（MPa）：

17.133637

定压井w2产量（m^3/d）：

68.382573

============================================================================================

时刻tt（d）：

30

压力分布（MPa）：

19.087318.980718.800618.873218.939718.9667

19.164019.016218.570018.884418.976418.9917

19.327619.256019.112519.112319.083119.0297

19.531319.495519.415719.337619.200119.005718.775118.738918.7598

19.759919.743819.685319.549419.326918.988018.598818.687118.7791

20.000020.000020.000019.763619.470519.002618.024918.656318.8907

20.000019.701419.402019.064019.149419.2412

20.000019.870019.741319.636019.635619.6463

20.000020.000020.000020.000020.000020.0000

定产井w1井底流压（MPa）：

17.060370

定压井w2产量（m^3/d）：

68.032260

============================================================================================

时刻tt（d）：

60

压力分布（MPa）：

19.087018.980318.800318.872818.939318.9664

19.163719.015918.569718.884118.976118.9913

19.327319.255819.112319.112019.082919.0294

19.531119.495319.415619.337419.199919.005518.775018.738718.7597

19.759819.743719.685219.549319.326818.987918.598718.687018.7790

20.000020.000020.000019.763619.470519.002618.024818.656218.8906

20.000019.701419.402019.063919.149419.2411

20.000019.870019.741319.635919.635619.6462

20.000020.000020.000020.000020.000020.0000

定产井w1井底流压（MPa）：

17.060072

定压井w2产量（m^3/d）：

68.030902

============================================================================================

时刻tt（d）：

90

压力分布（MPa）：

19.087018.980318.800318.872818.939318.9664

19.163719.015918.569718.884118.976118.9913

19.327319.255819.112319.112019.082919.0294

19.531119.495319.415619.337419.199919.005518.775018.738718.7597

19.759819.743719.685219.549319.326818.987918.598718.687018.7790

20.000020.000020.000019.763619.470519.002618.024818.656218.8906

20.000019.701419.402019.063919.149419.2411

20.000019.870019.741319.635919.635619.6462

20.000020.000020.000020.000020.000020.0000

定产井w1井底流压（MPa）：

17.060072

定压井w2产量（m^3/d）：

68.030901

============================================================================================

时刻tt（d）：

180

压力分布（MPa）：

19.087018.980318.800318.872818.939318.9664

19.163719.015918.569718.884118.976118.9913

19.327319.255819.112319.112019.082919.0294

19.531119.495319.415619.337419.199919.005518.775018.738718.7597

19.759819.743719.685219.549319.326818.987918.598718.687018.7790

20.000020.000020.000019.763619.470519.002618.024818.656218.8906

20.000019.701419.402019.063919.149419.2411

20.000019.870019.741319.635919.635619.6462

20.000020.000020.000020.000020.000020.0000

定产井w1井底流压（MPa）：

17.060072

定压井w2产量（m^3/d）：

68.030901

============================================================================================

时刻tt（d）：

360

压力分布（MPa）：

19.087018.980318.800318.872818.939318.9664

19.163719.015918.569718.884118.976118.9913

19.327319.255819.112319.112019.082919.0294

19.531119.495319.415619.337419.199919.005518.775018.738718.7597

19.759819.743719.685219.549319.326818.987918.598718.687018.7790

20.000020.000020.000019.763619.470519.002618.024818.656218.8906

20.000019.701419.402019.063919.149419.2411

20.000019.870019.741319.635919.635619.6462

20.000020.000020.000020.000020.000020.0000

定产井w1井底流压（MPa）：

17.060072

定压井w2产量（m^3/d）：

68.030901

============================================================================================

时刻tt（d）：

720

压力分布（MPa）：

19.087018.980318.800318.872818.939318.9664

19.163719.015918.569718.884118.976118.9913

19.327319.255819.112319.112019.082919.0294

19.531119.495319.415619.337419.199919.005518.775018.738718.7597

19.759819.743719.685219.549319.326818.987918.598718.687018.7790

20.000020.000020.000019.763619.470519.002618.024818.656218.8906

20.000019.701419.402019.063919.149419.2411

20.000019.870019.741319.635919.635619.6462

20.000020.000020.000020.000020.000020.0000

定产井w1井底流压（MPa）：

17.060072

定压井w2产量（m^3/d）：

68.030901

============================================================================================

时刻tt（d）：

1800

压力分布（MPa）：

19.087018.980318.800318.872818.939318.9664

19.163719.015918.569718.884118.976118.9913

19.327319.255819.112319.112019.082919.0294

19.531119.495319.415619.337419.199919.005518.775018.738718.7597

19.759819.743719.685219.549319.326818.987918.598718.687018.7790

20.000020.000020.000019.763619.470519.002618.024818.656218.8906

20.000019.701419.402019.063919.149419.2411

20.000019.870019.741319.635919.635619.6462

20.000020.000020.000020.000020.000020.0000

定产井w1井底流压（MPa）：

17.060072

定压井w2产量（m^3/d）：

68.030901

============================================================================================

压力达到稳定时刻（d）：

70

压力分布（MPa）：

19.087018.980318.800318.872818.939318.9664

19.163719.015918.569718.884118.976118.9913

19.327319.255819.112319.112019.082919.0294

19.531119.495319.415619.337419.199919.005518.775018.738718.7597

19.759819.743719.685219.549319.326818.987918.598718.687018.7790

20.000020.000020.000019.763619.470519.002618.024818.656218.8906

20.000019.701419.402019.063919.149419.2411

20.000019.870019.741319.635919.635619.6462

20.000020.000020.000020.000020.000020.0000

定产井w1井底流压（MPa）：

17.060072

定压井w2产量（m^3/d）：

68.030901

三、隐式运算结果

x方向步长（m）：

200.0000

y方向步长（m）：

200.0000

时间步长（h）：

2.0000

============================================================================================

时刻tt（d）：

0

压力分布（MPa）：

20.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.0000

20.000020.000020.000020.000020.000020.0000

定产井w1井底流压（MPa）：

18.490414

定压井w2产量（m^3/d）：

112.454746

============================================================================================

时刻tt（d）：

10

压力分布（MPa）：

19.138919.031618.849318.918618.982419.0080

19.208819.059418.609518.919119.006319.0189

19.362019.287719.138619.129519.090919.0307

19.553919.514519.426319.332019.174018.963818.716518.677418.6993

19.770619.751119.682919.520719.253818.901318.522618.617118.7141

20.000020.000020.000019.722019.322218.858117.936818.577918.8220

20.000019.378619.161218.921219.061419.1766

20.000018.988919.414219.499819.570619.6051

20.000020.000020.000020.000020.000020.0000

定产井w1井底流压（MPa）：

17.099934

定压井w2产量（m^3/d）：

66.051928

============================================================================================

时刻tt（d）：

30

压力分布（MPa）：

19.051918.942318.756418.821518.880818.9034

19.130618.979518.526418.831818.914718.9243

19.298219.222919.070919.058019.015218.9529

19.508519.468619.378119.280819.118918.909618.672318.636418.6588

19.747119.727619.657519.490519.21