威布尔模型与乙型水驱曲线的联解法——开发课程设计报告(编辑修改稿)

威布尔模型与乙型水驱曲线的联解法——开发课程设计报告(编辑修改稿)内容摘要：

年产油量、计算油田累积产油量、计算含水率如下表： 表 41 数据表 时间（ a） 实际产量 预测产量 实际累计产量 预测累计产量 实际含水率 预测含水率 1 0 2 0 3 0 4 0 5 0 6 0 7 0 8 0 9 0 10 11 12 13 14 15 16 17 18 19 20 87 21 2747 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 10 41 42 43 44 45 46 47 48 49 50 51 计算可采储量： NR =（ 104t） ，计算最高年产量发生的时间： ， 计算最高年产量： Qmax=(104t)。 预测的理论产油量、累积产油量和含水率，与实际产油量、累积产油量和含水率进行对比。 按照（ 318）式的关系在同一坐标系中绘制实际产量与预测产量及其对应时间 t 之间的关系曲线， 绘于图 43。 图 43 产量与时间关系曲线 11 按照（ 319）式的关系在同一坐标系中绘制实际累计产量与预测累 计产量及其对应时间 t 之间的关系曲线，绘于图 44。 图 44 累计产量与时间关系曲线 按照（ 328）式的关系在同一坐标系中绘制实际含水率与预测含水率及其对应时间 t之间的关系曲线，绘于图 45。 图 45 含水率与时间关系曲线 12 认识与结论 通过本次课程设计，得出如下认识与结论： （ 1）通过此次 油田开发的 课程设计， 我 基本掌握了威布尔（ Weibull）模 型与乙型水驱曲线联解法，并能初步对油田进行预测分析。 （ 2） 通过这次的课程设计，培养了我独立思考和动手操作的能力，也使我在日常的学习与讨论当中增长了许多课外知识，认识到数学模型对油气田的产量预测的重要性。 虽然在课程设计的过程中遇到了一些问题，但经过反复的思考与检查找出原因所在并进行改正。 （ 3）通过计算机的编程学习与操作，使我对 VB 语言有了更进一步的认识和了解。 同时，提高了使用计算机技应用的能力以及解决实际问题的工程应用能力，使理论与实践相结合。 （ 4） 现在能熟练掌握、运用威布尔（ Weibull）模型与乙型 水驱曲线联解法的原理、推导及公式，可以对油气田的年产量、累计产量、含水率、可采储量、最高年产量、最高年产量发生的时间以及与最高年产量对应的累积产量进行预测。 （ 5） 这次课程设计，使我认识到学好计算机的重要性与实践操作的必要性， 在老师的指导下和与 同学的讨论中，不断地发现问题、不断地解决问题、不断地领悟、不断地获取，使自己的程序不断完善。 （ 6）课程设计诚然是一门专业课，它给了我很多的专业知识以及专业上的提升。 课程设计同时又是一门讲道课，一门思辨课，给了我许多道，给了我许多思，得到很大的提升空间。 13 参考文献 [1]陈元千 .预测油气田产量的 Weibull 模型 [J].新疆石油地质， 1995， 16（ 3）： 250255. [2]童宪章 .天然水驱和人工注水油藏的统计规律探讨 [J].石油勘探与开发， 1978， 4（ 6）：3864. [3] 陈元千 .一种新型水驱曲线关系式的推导及应用 [J].石油学报， 1993， 14（ 2）： 6573. [4] 陈元千 .水驱曲线关系式的推导 [J]. 石油学报， 1985， 6（ 2）： 6978. 14 附录：计算机程序与结果 计算 机 程序 Option Base 0 Const e = Dim a, t(200), qt(200), t0(200), q0(200), qw(200), ye(200), wp(200) Dim b, c, d Dim afw, qmax, tm, NR, x, y, a1, a2, r, rmax, bmax, aa1, bb1, cc1, ub, aa, ab, ar, br, yt, an, bn, mn Dim s(), px(), py() Dim sx, sy, sxy, sx2, sy2 Dim np(200), lp(200) Dim yq0(200), ynp(200), yfw(200), ynr, fw(200) Private Sub c1_Click() yt = 50 Rem 源数据输入 Open + \ For Input As 1 ub = 1 Do Until EOF(1) ub = ub + 1 Input 1, ye(ub), t(ub), q0(ub), qw(ub), np(ub), wp(ub), lp(ub) lp(ub) = Log(lp(ub)) / Log(10) Loop Close 1 Private Sub c5_Click() Rem 乙型参数计算 For an = ub To 2 Step 1 39。 求相关系数 r,趋势线 sx = 0: sy = 0: sxy = 0: sx2 = 0: sy2 = 0 l = ub an For i = l To an + l sx = sx + np(i) Next i 15 For i = l To an + l sy = sy + lp(i) Next i For i = l To an + l sxy = sxy + np(i) * lp(i) Next i For i = l To an + l sx2 = sx2 + np(i) ^ 2 Next i For i = l To an + l sy2 = sy2 + lp(i) ^ 2 Next i a2 = ((an + 1) * sxy sx * sy) / ((an + 1) * sx2 sx ^ 2) a1 = (sx2 * sy sx * sxy) / ((an + 1) * sx2 sx ^ 2) r = ((an + 1) * sxy sx * sy) / (((an + 1) * sx2 sx ^ 2) ^ * ((an + 1) * sy2 sy ^ 2) ^ ) 39。 取最大值 r If r ar Then ar = r Next an aa = a1 ab = a2 = aa 39。 截距，参数 A = ab 39。 斜率，参数 B = ar 39。 最大相关系数 Rem weibull 参数计算 bn = 12 l = ub bn bmax = 0 rmax = 0 16 For b = 1 To 2 Step sx = 0: sy = 0: sxy = 0: sx2 = 0: sy2 = 0 For i = l To l + bn t(i) = i + 1 Next i For i = l To l + bn qt(i) = Log(q0(i) / (t(i) ^ b)) / Log(10) t0(i) = t(i) ^ (b + 1) Next i For i = l To l + bn sx = sx + t0(i) Next i For i = l To l + bn sy = sy + qt(i) Next i For i = l To l + bn sxy = sxy + qt(i) * t0(i) Next i For i = l To l + bn sx2 = sx2 + t0(i) ^ 2 Next i For i = l To l + bn sy2 = sy2 + qt(i) ^ 2 Next i r = Abs(((bn + 1) * sxy sx * sy) / (((bn + 1) * sx2 sx ^ 2) ^ * ((bn + 1) * sy2 sy ^ 2) ^ )) 39。 最大相关系数 If r rmax Then 39。 r 取最大值 rmax = r bmax = b 17 a2 = ((bn + 1)。