注水开发储采比研究报告——东北石油大学最终版(编辑修改稿)

注水开发储采比研究报告——东北石油大学最终版(编辑修改稿)内容摘要：

8 0 0 . 0 0 01 0 0 0 . 0 0 01 2 0 0 . 0 0 01 4 0 0 . 0 0 01 6 0 0 . 0 0 01 8 0 0 . 0 0 00 1 0 0 0 2 0 0 0 3 0 0 0 4 0 0 0 5 0 0 0 6 0 0 0 7 0 0 0 8 0 0 0储采比时间产量上升阶段储采比与产量的变化曲线系列 1 图 59 产量上升阶段储采比与产量的变化曲线 稳产阶段储采比与开发时间的关系曲线，如图 510 所示。 图 510 稳产阶段储采比与开发时间的关系曲线 东北石油大学课程设计（报告） 18 递减阶段储采比与产量的关系曲线，如图 511 所示。 图 511 递减阶段储采比与产量的关系曲线 1全程储采比与开发时间的关系曲线，如图 512 所示。 图 512 全程储采比与开发时间的关系曲线 设计结果及分析 产量上升阶段： 产量上升阶段储采比和产量呈现双对数 的直线关系，直线关系式为lg 1. 06 85 lg 5. 44 02Q   。 符合其理论直线关系为 lg lgAQ 。 稳产阶段： 东北石油大学课程设计（报告） 19 稳 产 阶 段 储 采 比 与 时 间 呈 直 线 递 减 关 系 ， 直 线 关 系 式 为 59 8 34 .30 33t   。 其理论直线关系为 Et，经对比可知理论描述的稳产阶段储采比的变化特征较为合理。 递减阶段： 递减阶段 储采比和 产量呈 现双对 数的直线 关系， 直线关 系式为lg 0. 48 73 1 lg 1. 02 33Q ，利用直线斜率可求得递减指数 。 其理论直线关系为 lg lg Q   ，经对比可知理论描述的递减阶段储采比的变化关系较为合理。 并由此推断得的递减指数应为正确值。 全程： 全程储采 比和时间 呈双对 数的直 线关系， 直线关 系式可 表示为lg lg    。 其理论直线关系为 lg lgA B t  经对比可知理论描述的全程储采比的变化规律合理有效。 求威布尔模型中的 ,ab和 c 的值，及利用 ,ab和 c 的值求得可采储量 RN 值。 直线关系为： 5 2 . 7 2 2 91 . 7 2 2 9l g 8 . 0 4 2 1 0 1 . 9 7 0 3Q tt    。 根据拟合的直线的截距和斜率可求出： α 1. 970310 10    511 5 3 9 9 . 3 12 . 3 0 3 2 . 3 0 3 8 . 0 4 2 1 0c    由于  ，将 ,ab和 c 的值代入求可采储量 RN 公式得可采储量为： 4R 9 3 . 3 9 5 3 9 9 . 3 1 1 8 5 1 7 7 . 5 ( 1 0 t )1 2 . 7 2 2 9acN b    由威布尔模型预测可采储量 RN 值： 4R (10 t )N  ，实际给的已知4R (10 t )N  ，相对误差为 %，这说明威布尔模型预测可采储量有效。 计算结果的 认识 根据各阶段用实际数据的拟合的直线关系与相应的理论关系式进行对比分析，可知各阶段描述的储采比与时间的变化关系合理有效。 威布尔预测模型预测注水开发油田可采储量与实际值相当误差较小，说明此法预测的可采储量较为准东北石油大学课程设计（报告） 20 确。 威布尔预测模型能很好的预测注水开发油田产量、累积产量、可采储量、储采比随时间的变化关系，但是不能预测油田的含水率、产水量、产液量及累积产水量和累积产液量，而这些开发指标却 是水驱开发油田所需要预测的。 乙型水驱曲线法是油藏工程中重要的预测方法，但它只能预测累积产水量和累积产油量之间的关系，却不能预测开发指标与时间的关系，而油田开发指标的预测，都离不开与时间的关系。 因此如果将两者结合起来进行预测，不仅可以完善各种的理论，而且更有利于预测出更多的油田开发的相关参数。 因此我们应该用威布尔模型和乙型水驱曲线联解来预测油田开发指标。 结论 产量上升阶段： 产量上升阶段储采比和产量呈现双对数的直线关系。 稳产阶段： 稳产阶段储采比与时间呈直线递减关系。 递减阶段： 递 减阶段储采比和产量呈现双对数的直线关系。 全程： 全程储采比和时间呈双对数的直线关系。 威布尔模型预测可采储量 RN 值： 由威布尔模型预测可采储量 RN 值： 4R (10 t )N  ，实际给的已知4R (10 t )N  ，相对误差为 %，这说明威布尔模型预测可采储量有效。 东北石油大学课程设计（报告） 21 参考文献 1. 陈元千，双曲线递减的简化及确定 可采储量的截距法，天然气工业 [J， （ 4）， 3237。 2. 陈元千，赵庆飞，油气田储采比变化关系的研究，断块油气田 [J]， （ 6），2326； 3. 陈元千，胡建国，预测油气田产量和可采储量 weibull 模型，新疆石油地质，（ 3）， 250255； 4. 陈元千，油田可采储量计算方法，新疆石油地质 [J]， （ 2）， 130137； 5. 陈元千，油气藏工程实践，石油工业出版社 [M],。 东北石油大学课程设计（报告） 22 附录（计算机程序） 一、 程序运行界面 数据导入及计算结果，如附图 1 所示。 附图 1 数据导入及计算结果 绘制产量、累积产量和储采比随时间的变化曲线程序界面，如附图 2 所示 .。 东北石油大学课程设计（报告） 23 附图 2 绘制产量、累积产量和储采比随时间的变化曲线程序界面 产量上升阶段的计算结果和图形输出，如附图 3 所示。 附图 3 产量上升阶段的计算结果和图形输出 东北石油大学课程设计（报告） 24 稳产阶段的计算结果和图形输出，如附图 4 所示。 附图 4 稳产阶段的计算结果和图形输出 递减阶段的计算结果和 图形输出，如附图 5 所示。 东北石油大学课程设计（报告） 25 附图 5 递减阶段的计算结果和图形输出 全程的计算结果和图形输出，如附图 6 所示 附图 6 全程的计算结果和图形输出 利用威布尔模型求可采储量 RN 的计算结果和图形输出，如附图 7 所示 东北石油大学课程设计（报告） 26 附图 7 利用威布尔模型求可采储量 RN 的计算结果和图形输出 二、计算机程序代码 Form1 代码： Dim t() As Single, nian() As Single, Q() As Single, NP() As Single, w() As Single Dim A As Single, B As Single, QS As Single, QAa As Single Dim xx As String, yy As String Dim XMAX As Single, YMAX As Single Dim XMIN As Single, YMIN As Single Dim x() As Single, y() As Single Sub ercheng(x() As Single, y() As Single, m As Integer, n As Integer) 39。 定义最小二乘法过程进行线性拟合 Dim i As Integer, R As Single Dim ER() As Single ReDim ER(5, m 1 To n) XMAX = x(m): YMAX = y(m) 39。 求出横坐标和纵坐标的最大，最小值。 画图使用 XMIN = x(m): YMIN = y(m) 东北石油大学课程设计（报告） 27 For i = m To n If XMAX x(i) Then XMAX = x(i) If YMAX y(i) Then YMAX = y(i) If XMIN x(i) Then XMIN = x(i) If YMIN y(i) Then YMIN = y(i) Next i For i = m To n ER(1, i) = x(i) ER(2, i) = y(i) ER(3, i) = ER(1, i) * ER(2, i) ER(4, i) = ER(1, i) ^ 2 ER(5, i) = ER(2, i) ^ 2 Print xx。 (。 i。 ) =。 Format(ER(1, i), ), yy。 (。 i。 ) =。 Format(ER(2, i), ) If i = 49 Then CurrentX = / 2 CurrentY = 0 ElseIf i 49 Then CurrentX = / 2 End If ER(1, m 1) = ER(1, m 1) + ER(1, i) ER(2, m 1) = ER(2, m 1) + ER(2, i) ER(3, m 1) = ER(3, m 1) + ER(3, i) ER(4, m 1) = ER(4, m 1) + ER(4, i) ER(5, m 1) = ER(5, m 1) + ER(5, i) Next i ER(1, m 1) = ER(1, m 1) / (n m + 1) ER(2, m 1) = ER(2, m 1) / (n m + 1) R = (ER(3, m 1) (n m + 1) * ER(1, m 1) * ER(2, m 1)) / 东北石油大学课程设计（报告） 28 _ Sqr((ER(4, m 1) (n m + 1) * ER(1, m 1) ^ 2) * (ER(5, m 1) (n m + 1) * ER(2, m 1) ^ 2)) A = ((n m + 1) * ER(1, m 1) * ER(2, m 1) ER(3, m 1)) / ((n m + 1) * ER(1, m 1) ^ 2 ER(4, m 1)) B = ER(2, m 1) A * ER(1, m 1) If i 49 Then 拟合公式为： 39。 输出拟合公式和相关性系数 If A = 1 Then yy。 =。 xx Else yy。 =。 Format(A, )。 xx。 End If If B 0 Then +。 Format(B, ) ElseIf B 0 Then Format(B, ) End If 线性相关系数 R=。 Format(Abs(R), ) Else 拟合公式为： If A = 1 Then yy。 =。 xx Else yy。 =。 Format(A, )。 xx。 End If If B 0 Then 东北石油大学课程设计（报告） 29 +。 Format(B, ) ElseIf B 0 Then Format(B, ) End If 线性相关系数 R=。 F。