注水开发储采比研究报告——东北石油大学最终版本.doc

东 北 石 油 大 学课 程 设 计课 程 海洋油气工程课程设计 题 目 注水开发油田储采比变化特征的研究 院 系 石油工程学院 专业班级 海油12-1班 学生姓名 于洪楠 学生学号 00000000000 指导教师 00000000000 2015年7月24日东北石油大学课程设计任务书课程 海洋油气工程课程设计题目 注水开发油田储采比变化特征的研究专业 海洋油气工程 姓名 于洪楠 学号 00000000000主要内容、基本要求、主要参考资料等主要内容：(1) 假设产量上升阶段产量变化符合指数关系，推导产量上升阶段储采比与产量之间的关系；(2) 推导稳产阶段储采比与开发时间的关系；(3) 假设油田产量递减符合双曲递减，推导递减阶段储采比与产量之间的关系；(4) 假设油田产量可由威布尔预测模型给出，推导全程储采比与开发时间的关系；(5) 根据给定的可采储量和产量历史数据，计算不同时刻累积产量和储采比，并绘制产量、累积产量和储采比随时间的变化曲线；(6) 根据(5)中计算结果数据，编制计算机程序，进行线性回归，得到产量递减阶段的递减指数、威布尔预测模型中各参数及(1)(4)中各关系式的系数；(7) 绘制产量上升阶段储采比与时间的关系曲线；(8) 绘制稳产阶段储采比与开发时间的关系曲线；(9) 绘制递减阶段储采比与产量之间的关系曲线；(10) 绘制全程储采比与开发时间的关系曲线。基本要求：1)基础数据某油田或区块的可采储量和产量历史数据。2)设计要求该专题设计最终要求是，学生通过自选上述基础数据，利用所学知识完成规定设计内容，编制相应软件，并提交规范设计报告。主要参考资料1 陈元千.油藏工程实践M.石油工业出版社. 完成期限 2015年7月6日2015年7月24日 指导教师 李占东 崔晓娜 专业负责人 杨二龙 2015年7月 6日目 录1 前 言11.1 设计的目的意义11.2 设计的主要内容12 基础数据23 基础理论33.1 产量上升阶段储采比变化特征33.2 稳产阶段储采比变化特征43.3 递减阶段储采比变化关系43.4 全程储采比变化关系43.5 确定a，b和c常数方法44 程序设计框图54.1 计算机程序设计框图54.2 计算机程序55 设计结果及分析6认识与结论7参考文献8附录：计算机程序与结果91 前 言1.1 设计的目的意义储采比是油田开发的重要指标，它不仅直接反映了油田的开发程度，而且与油田的稳产形式密切相关。因此本次课程设计运用指数模型、双曲线模型和威布尔模型，对于注水开发油田的产量上升阶段提出了储采比与产量的关系、稳产阶段储采比与时间的关系、递减阶段储采比与产量关系和全程储采比与时间的关系。并通过实例运用此方法和验证此方法的有效性。1.2 设计的主要内容(1) 假设产量上升阶段产量变化符合指数关系，推导产量上升阶段储采比与产量之间的关系；(2) 推导稳产阶段储采比与开发时间的关系；(3) 假设油田产量递减符合双曲递减，推导递减阶段储采比与产量之间的关系；(4) 假设油田产量可由威布尔预测模型给出，推导全程储采比与开发时间的关系；(5) 根据给定的可采储量和产量历史数据，计算不同时刻累积产量和储采比，并绘制产量、累积产量和储采比随时间的变化曲线；(6) 根据(5)中计算结果数据，编制计算机程序，进行线性回归，得到产量递减阶段的递减指数、威布尔预测模型中各参数及(1)(4)中各关系式的系数；(7) 绘制产量上升阶段储采比与时间的关系曲线；(8) 绘制稳产阶段储采比与开发时间的关系曲线；(9) 绘制递减阶段储采比与产量之间的关系曲线；（10）绘制全程储采比与开发时间的关系曲线。2 基础数据某油田的开发数据如下：2-1 某油田产量历史数据可采储量（万吨）180419.1序号年年产量（万吨）2519776894.311953113.62619786564.221954343.8271979619131955655.62819805784.1419561033.72919815353.1519571466.63019824907.6619581944.33119834456.5719592456.63219844008.5819602993.73319853571.19196135453419863150.91019624099.93519872753.21119634647.33619882382.41219645176.53719892041.21319655676.73819901731.71419666137.63919911454.41519676549.74019921209.31619686904.5411993995.31719697195421994810.81819707631.7431995653.71919717631.7441996521.62019727631.7451997411.82119737631.7461998321.72219747631.7471999248.62319757391.54820001902419767172.7492001143.73 基础理论3.1 产量上升阶段储采比变化特征Weibull(威布尔)于1939年提出的统计分布模型已成为生命试验和可靠性理论研究的基础。该模型的分布密度表示为： (3-1)式中 f(x)威布尔分布的分布密度函数；x分布变量，根据实际问题，分布区间为0；控制分布形态的形状参数；控制分布峰位和峰值的尺度参数。若对(3-l)式进行积分，在x为0区间内，可以得到Weibu11的分布函数值等于1，推证如下： （3-2）为将Weibull分布模型用于油气田开发指标的预测，将(3-l)式改写为 (3-3)式中 Q油气田的年产量，104t/a（油）或108m3/a（气）； t油气田的开发时间，a； C由Weibull分布模型转换为油气田开发实用模型的模型转换常数。 油气田的累积产量表达式为： （3-4）式中 NP油气田的累积产量，104t或108t(油)；108m(气)。若产量上升阶段的产量随时间的关系符合如下指数递减关系： （3-5）则产量上升阶段的累积产量可表示为： （3-6）将（3-6）式代入（3-7）式得： (3-7)将（3-8）式代入储采比的定义式中得到下式： (3-8)对于产量上升阶段，因此有： (3-9)对（3-9）式两边取常用对数得： (3-10)则（3-10）式可改写为： (3-11)其中：。3.2 稳产阶段储采比变化特征在稳产阶段某一开发时间的储采比表示为： (3-12)在稳产阶段某一开发时间的累积产量可表示为： (3-13)将（3-13）式代入（3-12）式可得，在稳产阶段储采比与开发时间的关系式为： (3-14)则（3-14）式可写成： (3-15)其中：。3.3 递减阶段储采比变化关系进入递减阶段之后的储采比可表示为： (3-16)进入递减阶段累积产量与产量之间的关系如下： (3-17)则（3-17）式可以简化为： (3-18)其中：。将（3-22）式代入（3-20）式可得： (3-19)由于（3-22）式可知： (3-20)注水开发油田的可采储量可表示为： (3-21)由（3-20）式和（3-21）式可得： (3-22)将（3-22）式代入（3-23）式可得： (3-23)将（3-23）式两边取常用对数得： (3-24)则（3-24）式可改写为： (3-25)其中：。 3.4 全程储采比变化关系假设全程产量、累积产量、和可采储量随时间的关系可由威布尔预测模型表示： (3-26) (3-27) (3-28)将（3-26）、（3-27）、（3-28）式代入储采比的定义式中得到下式： (3-29)将（3-29）式两边取常用对数得： (3-30)则（3-30）式可简化为： (3-31)其中：。3.5 确定和常数方法由产量随时间的关系可由威布尔预测模型表示： (3-32)将（3-32）式改写成： (3-33)将（3-33）是等号两端取常用对数得： (3-34)则（3-34）式可简化为： (3-35)其中：那么： (3-36) (3-37)根据油田开发的实际产量和相应的时间的数据，给出不同的值，利用（3-35）式进行线性迭代试差，能得到相关系数最高的值，即为所求的值。由线性回归得到直线的斜率和截距值后，再由（3-36）式和（3-37）式求出和的数值。目录4 程序设计框图4.1 计算机程序设计框图1、最小二乘法程序设计框图，见图4-1。输入x和y值分别计算计算计算直线的斜率计算直线的截距计算线性相关系数输出拟合公式开始结束图4-1 最小二乘法程序设计框图2、产量上升阶段程序设计框图，如图4-2所示开始输入产量上升阶段起始点计算调用最先二乘法进行直线拟合绘出散点图和拟合直线图利用直线截距求出可采储量结束图4-2 产量上升阶段程序设计框图2、稳产阶段程序设计框图，如图4-3所示开始输入稳产阶段起始点调用最先二乘法进行直线拟合绘出散点图和拟合直线图结束图4-3 稳产阶段程序设计框图3、递减阶段程序设计框图，如图4-4所示。开始输入产量递减阶段起始点计算调用最先二乘法进行直线拟合绘出散点图和拟合直线图利用直线斜率求出递减指数n结束图4-4 递减阶段程序设计框图5、全程程序设计框图，如图4-5所示。开始输入全程起始点计算调用最先二乘法进行直线拟合绘出散点图和拟合直线图利用直线斜率求出b值输入b值调用最先二乘法进行直线拟合绘出散点图和拟合直线图利用直线截距与斜率求出值结束相关度好不好否是·图4-5 全程程序设计框图4.2 计算机程序计算机程序代码，见附录。5 设计结果及分析1、利用已知数据求得累积产量和储采比，见表5-1。表5-1利用已知数据求得累积产量和储采比可采储量（万吨）180419.1序号年年产量（万吨）Np（万吨）w11953113.6113.61587.221954343.8457.4523.431955655.61113273.5419561033.72146.7172.5519571466.63613.3120.6619581944.35557.689.9719592456.68014.270.2819602993.711007.956.691961354514552.946.81019624099.918652.839.51119634647.323300.133.81219645176.528476.629.41319655676.734153.325.81419666137.640290.922.81519676549.746840.620.41619686904.553745.118.3171969719560940.116.61819707631.768571.814.71919717631.776203.513.72019727631.783835.212.72119737631.791466.911.72219747631.799098.610.72319757391.5106490.110.02419767172.7113662.89.32519776894.3120557.18.72619786564.2127121.38.12719796191133312.37.62819805784.1139096.47.12919815353.1144449.56.73019824907.6149357.16.33119834456.5153813.66.03219844008.5157822.15.63319853571.1161393.25.33419863150.9164544.15.03519872753.2167297.34.83619882382.4169679.74.53719892041.2171720.94.33819901731.7173452.64.03919911454.41749073.84019921209.3176116.33.6411993995.3177111.63.3421994810.8177922.43.1431995653.7178576.12.8441996521.6179097.72.5451997411.8179509.52.2461998321.7179831.21.8471999248.6180079.81.4482000190180269.80.8492001143.7180413.50.02、绘制产量、累积产量和储采比随时间的变化曲线：图5-1产量与时间关系曲线图5-2 累积产量和时间的关系曲线 图5-3 储采比和时间的关系曲线3.产量上升阶段储采比与产量的关系曲线，如图5-4所示。图5-4 产量上升阶段储采比与产量的关系曲线4、稳产阶段储采比与开发时间的关系曲线，如图5-5所示。图5-5 稳产阶段储采比与开发时间的关系曲线5、递减阶段储采比与产量的关系曲线，如图5-6所示。 图5-6 递减阶段储采比与产量的关系曲线6、全程储采比与开发时间的关系曲线，如图5-7所示。图5-7 全程储采比与开发时间的关系曲线7、求威布尔模型中的和的值关系曲线，如图5-8所示。图5-8 求威布尔模型中的和的值关系曲线8、产量上升阶段储采比与产量的变化曲线，如图5-9所示。图5-9 产量上升阶段储采比与产量的变化曲线9、稳产阶段储采比与开发时间的关系曲线，如图5-10所示。图5-10 稳产阶段储采比与开发时间的关系曲线10、递减阶段储采比与产量的关系曲线，如图5-11所示。图5-11 递减阶段储采比与产量的关系曲线11、全程储采比与开发时间的关系曲线，如图5-12所示。图5-12 全程储采比与开发时间的关系曲线5.2 设计结果及分析1、产量上升阶段：产量上升阶段储采比和产量呈现双对数的直线关系，直线关系式为。符合其理论直线关系为。2、稳产阶段：稳产阶段储采比与时间呈直线递减关系，直线关系式为。其理论直线关系为，经对比可知理论描述的稳产阶段储采比的变化特征较为合理。3、递减阶段：递减阶段储采比和产量呈现双对数的直线关系，直线关系式为，利用直线斜率可求得递减指数。其理论直线关系为，经对比可知理论描述的递减阶段储采比的变化关系较为合理。并由此推断得的递减指数应为正确值。4、全程：全程储采比和时间呈双对数的直线关系，直线关系式可表示为。其理论直线关系为经对比可知理论描述的全程储采比的变化规律合理有效。5、求威布尔模型中的和的值，及利用和的值求得可采储量值。直线关系为：。根据拟合的直线的截距和斜率可求出：由于，将和的值代入求可采储量公式得可采储量为：由威布尔模型预测可采储量值：，实际给的已知，相对误差为2.38%，这说明威布尔模型预测可采储量有效。5.3 计算结果的认识根据各阶段用实际数据的拟合的直线关系与相应的理论关系式进行对比分析，可知各阶段描述的储采比与时间的变化关系合理有效。威布尔预测模型预测注水开发油田可采储量与实际值相当误差较小，说明此法预测的可采储量较为准确。威布尔预测模型能很好的预测注水开发油田产量、累积产量、可采储量、储采比随时间的变化关系，但是不能预测油田的含水率、产水量、产液量及累积产水量和累积产液量，而这些开发指标却是水驱开发油田所需要预测的。乙型水驱曲线法是油藏工程中重要的预测方法，但它只能预测累积产水量和累积产油量之间的关系，却不能预测开发指标与时间的关系，而油田开发指标的预测，都离不开与时间的关系。因此如果将两者结合起来进行预测，不仅可以完善各种的理论，而且更有利于预测出更多的油田开发的相关参数。因此我们应该用威布尔模型和乙型水驱曲线联解来预测油田开发指标。5.4结论1、产量上升阶段：产量上升阶段储采比和产量呈现双对数的直线关系。2、稳产阶段：稳产阶段储采比与时间呈直线递减关系。3、递减阶段：递减阶段储采比和产量呈现双对数的直线关系。4、全程：全程储采比和时间呈双对数的直线关系。5、威布尔模型预测可采储量值：由威布尔模型预测可采储量值：，实际给的已知，相对误差为2.38%，这说明威布尔模型预测可采储量有效。参考文献1. 陈元千，双曲线递减的简化及确定可采储量的截距法，天然气工业J，1994.14（4），32-37;2. 陈元千，赵庆飞，油气田储采比变化关系的研究，断块油气田J，1999.6（6），23-26；3. 陈元千，胡建国，预测油气田产量和可采储量weibull模型，新疆石油地质，1995.16（3），250-255；4. 陈元千，油田可采储量计算方法，新疆石油地质J，2000.12（2），130-137；5. 陈元千，油气藏工程实践，石油工业出版社M,2005.11;附录（计算机程序）一、程序运行界面1、数据导入及计算结果，如附图1所示。附图1 数据导入及计算结果2、绘制产量、累积产量和储采比随时间的变化曲线程序界面，如附图2所示.。附图2 绘制产量、累积产量和储采比随时间的变化曲线程序界面3、产量上升阶段的计算结果和图形输出，如附图3所示。附图3 产量上升阶段的计算结果和图形输出4、稳产阶段的计算结果和图形输出，如附图4所示。附图4 稳产阶段的计算结果和图形输出5、递减阶段的计算结果和图形输出，如附图5所示。附图5 递减阶段的计算结果和图形输出6、全程的计算结果和图形输出，如附图6所示附图6 全程的计算结果和图形输出7、利用威布尔模型求可采储量的计算结果和图形输出，如附图7所示附图7 利用威布尔模型求可采储量的计算结果和图形输出二、计算机程序代码Form1代码：Dim t() As Single, nian() As Single, Q() As Single, NP() As Single, w() As SingleDim A As Single, B As Single, QS As Single, QAa As SingleDim xx As String, yy As StringDim XMAX As Single, YMAX As SingleDim XMIN As Single, YMIN As SingleDim x() As Single, y() As SingleSub ercheng(x() As Single, y() As Single, m As Integer, n As Integer) '定义最小二乘法过程进行线性拟合 Dim i As Integer, R As Single Dim ER() As Single ReDim ER(5, m - 1 To n) XMAX = x(m): YMAX = y(m) '求出横坐标和纵坐标的最大，最小值。画图使用XMIN = x(m): YMIN = y(m) 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), "0.00"), yy; "(" i; ") = " Format(ER(2, i), "0.00") If i = 49 Then CurrentX = Form1.Width / 2 CurrentY = 0 ElseIf i > 49 Then CurrentX = Form1.Width / 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) / _ 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 Picture2.Print "拟合公式为：" '输出拟合公式和相关性系数 If A = 1 Then Picture2.Print yy; "=" xx Else Picture2.Print yy; "=" Format(A, "0.0000"); xx; End If If B > 0 Then Picture2.Print "+" Format(B, "0.0000") ElseIf B < 0 Then Picture2.Print Format(B, "0.0000") End If Picture2.Print Picture2.Print "线性相关系数R=" Format(Abs(R), "0.0000") Else Picture2.Print "拟合公式为：" If A = 1 Then Picture2.Print yy; "=" xx Else Picture2.Print yy; "=" Format(A, "0.0000"); xx; End If If B > 0 Then Picture2.Print "+" Format(B, "0.0000") ElseIf B < 0 Then Picture2.Print Format(B, "0.0000") End If Picture2.Print "线性相关系数R=" Format(Abs(R), "0.0000") End IfEnd SubPrivate Sub Command1_Click() '导入数据，求解累积产量和储采比 Dim m As Integer Dim i As Integer Picture1.Visible = False Form1.Cls Picture2.Cls YZNR = InputBox("输入已知NR值", "数据", 180419.1) Print "时间", "年份", "产量", "累积产量", "储采比"m = InputBox("输入已知数据的个数", "数据", 49) If m > 1 ThenOn Error Resume Next ReDim t(m), nian(m), Q(m), NP(m), w(m) Open "E:已知数据.txt" For Input As #1 Do Until EOF(1) For j = 1 To 49 If j < 50 Then ReDim Preserve t(j), Q(j) Input #1, t(j), nian(j), Q(j) End If Next j Loop Close #1 For i = 1 To m NP(i) = NP(i - 1) + Q(i) w(i) = (YZNR - NP(i) / Q(i) Print t(i), nian(i), Q(i), Format(NP(i), "0.0"), Format(w(i), "0.0") If i = 49 Then '换列输出 CurrentX = Form1.Width / 2 CurrentY = 0 ElseIf i > 49 Then CurrentX = Form1.Width / 2 End IfNext i Else Print "输入已知数据的个数必须大于一" End If Command2.Enabled = True '控制按钮可用性 Command3.Enabled = True Command4.Enabled = True Command5.Enabled = TrueEnd SubPrivate Sub Command2_Click() '产量上升阶段 Dim NR As Single Dim m As Integer, n As Integer Form1.Cls Picture1.Visible = True Print xx = "lgQ"yy = "lgw" m = InputBox("输入数据的起点", "上升阶段", 1) '输入数据起点和终点 n = InputBox("输入数据的终点", "上升阶段", 17) ReDim x(m To n), y(m To n) For i = m To n x(i) = Log(Q(i) / Log(10) y(i) = Log(w(i) / Log(10) Next i Call ercheng(x(), y(), m, n) '调用最小二乘法过程 Picture1.Cls '绘图 Picture1.Scale (-2, Int(YMAX) + 2)-(Int(XMAX) + 2, -2) '定义坐标系 Picture1.DrawWidth = 2