油气田产量递减算法

1、、产量递减率的定义下面所要叙述和应用的递减率是瞬时递减率，和油田上常用的阶段递减率略有不同。另外还要叙述两种不 同的递减率定义方法。由于定义方法或表达方式不一样，经验公式的表达方法就不一样，用途也不一样。首先绘制产量与时间变化的关系曲线，如图6-2所示。从图中可以看出，产量是随时间而下降的，所谓递减率（declinerate）就是单位时间内的产量变化率，或单位时间内产量递减百分数，其方程式如下：1 9dQ Q?(6-1)式中a为产量递减率，用（mon-1）或（a-1）表 数表达和计算；dQ为 Qk Qk 1是从阶段初至 减值，其单位是（t/mon ）或（t/a）; dt为阶段初至 隔，月或年。

2、这就是递减率的定义。从公式（6-1）可以看出，递减率表示的是产量 一个小于1的数，单位是时间的倒数，这里的时间 所用的时间单位一致，例如，Q的单位是（t/mon）,递减率a也可以用图解法来确定，由图 6-2所 时刻的递减率时（此时瞬时产量为 Q1）,首先做（Q1,然后求切线的斜率为 Q/At,再除以此点的产率即得，这就是示，通常采用小 阶段末的产量递 阶段末的时间间下降的速度，是 单位应与产量中 则用月等等。示可见，当求t1 t1）点的切线，dQQdt1 QQ1 t(6-2)各油田的递减规律是不同的，同一油田在不同开发阶段的递减规律也不相同，因此首先需要对不同的递减 规律特性、表达方式和应用方

3、法有切实的了解。下面就对几种常见的和广泛应用的递减规律分别介绍。二、产量递减规律油田产量递减规律一般包括指数递减（exponentialdecline） 、双曲线递减（hyperbolicdecline） 、调和递减（harmonicdecline）和产量衰减曲线四种类型，产量与递减率的关系可用下式表示：aiQ tQin(6-3)式中Qi递减期初始产量，t/mon或t/a ;a初始时刻的递减率，1/mon或1/a;n=1时为调和递减，当 n=0时为指n一递减指数，是用于判断递减类型、确定递减规律的重要参数，当 数递减，当n=0.5时为产量衰减，当 0<n<1且nw0.5时为双曲线递

4、减。1 .指数递减规律当n=0时为指数递减规律，此时（6-3）式变为a=a,递减率a是一个常数，因此指数递减规律也叫常递减规律。这种形式的递减主要表现在某些封闭型弹性驱动油藏、重力驱动油藏和一些封闭气藏。它的优点是公式简单，使用方便，缺点是利用这一规律来预测油田今后的产量变化时一般不能外推很远，不然会引起较大的误差用它计算的可采储量通常比实际的偏低。由递减率的定义式（6-1）可得dQQ-a dt (6-4)如果当t=0时，即在开始递减的时刻产量为Qi,而在任一时刻的产量为Q （t）,则把公式（6-4）积分后,可以得到lnQ(t)Qiat (6-5)或者Q（t） Qie at（6-6）公式（6-

5、6）为指数递减的基本方程，对其两端取对数，得 ,二.：:-logQ tlog Qi %3026（6-7）logQ t A Bt （6-8）*I!方程式（6-8）表示，若油田产量服从指数递减规律，则在单对数坐标纸上，实际产量的对数与时间呈一直线关系，直线的截距 A即为log Qi,由此可求出公式（6-7）中的初始产量 Qi,而由直线的斜率 B可以求得产量递减率a为T0时，油田产量正好变为a=2.3026B （6-9） I由公式（6-7）还可以推导出一个重要概念，即递减周期的概念，设在某一时刻初始产量Qi的十分之一，则T0为一个递减周期。 I由定义将logQ t log Qi /10代入公式（6-

6、7）可得 aT0= 2.3026即2.3026To(6-10)丁 2.3026 -To 或 aa同样可定义半周期，即产量降为初始产量之半的时间Ti,同理可得出T1=0.69315/a； a0.69315(6-11)只要知道递减率则可计算出递减周期或半周期，同样根据递减周期或半周期的大小，都可以计算出递减率a。当所绘制的油田产量的对数与时间的关系为一直线时，则该油田产量服从指数递减规律。已知油田产量的 递减率时，可由公式（6-6）计算出今后任一时刻的产量Q （t）。相反地，也可预测产量递减到某一极限值时所经历的开发时间。根据累积产量的定义得：tNp 0Q（t）dt （6-12）代入Q （t）的表

7、达式并积分后得：Qiat ,、Np 1 e（6-13）a其中Np是从递减期开始起算的累积产量，以后的递减规律推导中相同。由公式（6-13）可知，当t趋近于无限大时，递减期最大累积产量，即递减期的剩余可采储量 Np应为Q / a。 p也就说，在求得递减期初始产量Qi和递减率以后，可采储量就可以计算出来。对公式（6-13）进行分析可以得出如下结果： j1）由于式（6-13）中的哥指数是负值，因而式中括号内第二项始终小于1,并随时间增大而逐步趋近于零，因此累积产量是逐步趋近于一恒定值的，按此公式得出的累积产量随时间变化曲线是一条减速递增曲线。2）有了递减期累积产量随时间变化公式，就可以计算今后任一时

8、刻的剩余累积产量Np,又知Qi/a为递减期可采储量，则公式（6-13）可写为：1Np ' a 1 e at （6-14） Qi它表示的是按递减期可采储量计算的采出程度随时间变化公式，其末值为1.0。3）公式（6-14）还可以写为另一形式，即 IIQ t Qi aND t （6-15） p它说明对于指数递减规律，产量与累积产量的关系是一条直线。2.双曲线递减规律当0<n<1且n w 0.5时为双曲线递减规律， 双曲线递减规律不但适用于各种天然驱动油藏，更主要的是适用于各种不同的水压驱动油藏。所以它有较广泛的适用性。双曲线递减规律的特点：递减率a随时间而变，而且越变越小，即愈接

9、近油田开发末期，递减愈慢 （如图6-3）。其表达式为：a t按照前面已经讲到的递减率的定义，则有:Qi(6-16)dQ _2_dt Q t(6-17)将上式代入公式(6-16)并经整理后可得aidt=-Qin , Q-(n+1)dQ若递减的初始条件是 t=0时Q=Qi,则任- 时间变化的表达式如下：aitQ叱巡(6-18)n n.& 4 2 O 8 6 4 .1 1-1常递减曲线；2双曲线递减曲线图6-3双曲线递减规律曲线1常递减曲线；2双曲线递减曲线时刻之产量随Q或者na1 1Qn最后可得双曲线型产量随时间变化的公式。QiQ t 万丁(6-19)1 nat不同油田产量的递减规律不同，

10、这表现在其初始递减率和递减指数不同，而又主要取决于n值的变化情况。影响n值的因素很多，如岩石性质、驱动方式、地质条件和开采方式等。所以n值变化范围很广，其选择和应用都比较麻烦。这样在实际应用中，更多地使用标准曲线对比法而不是解析方法来求解。根据累积产量的定义得：tN P t 0 Q t dt (6-20)代入Q(t)的表达式并积分后得：Qi1NP t i- 1 : (6-21)ai 1 n d i1 aint n根据递减率的定义可求得递减率随时间的变化规律dQ 1 一 口一一，a t代入Q(t)的表达式并求导，得到递减率的表达式为:dt Q t,ai一a t (6-22)1 nait从(6-2

11、2)式可以看出，双曲线递减的递减率与递减指数n和时间t有关，t越大，递减率a越小。3.调和递减规律在递减率的定义中，当n=1时为调和递减规律，调和递减规律可以看成是双曲线递减规律的一种特殊形式。调和递减规律的特点：递减率与递减产量成正比,其表达式为：a t Q t(6-23)aiQi调和递减产量随时间变化的公式为(令6-19中 n=1):Qi,、Q t (6-24)1 ait根据累积产量的定义得:tN P t 0 Q t dt (6-25)代入Q的表达式并积分后得:,Qi .,、NP t-ln 1ait (6-26)ai根据递减率的定义可得到递减率的表达式为:ai-一(6-27) ait从(6

12、-27)式可以看出，调和递减的递减率与时间t有关，t越大，递减率a越小。4.产量衰减规律将n=0.5代入(6-19)式得:当n=0.5时为产量衰减规律，产量衰减规律也是双曲线递减规律的一种特殊形式,QiQ t 2 (6-28)1 0.5ait上式可以进一步表达成:2式中：B 4Qiai c=2/ai递减类 型指数递减双曲线递减调和递减产量衰减递减指 数n=00<n<1,n 丰 0.5n=1n=0.5递减率a=a=常数a=ai(1+nait)-1a=ai(1+ait)-1a=ai(1+0.5ait)-1产量与 时间累积产量与时间当递减时间很长时，t>>c, t+c=t,上

13、式变为:表6-2四种递减类型对比表产量与 累积产 量Q(t)=(A-N p)2/B开发时 间t=2 (Qi/Q)1/2-1) /ai断其所属的递减类型，确定递减参数(a,ai和n),建立相关经验公式。目前经常采用的方法有，图解法、试凑法、曲线位移法和典型曲线拟合法。1 .图解法图解法就是将实际生产数据按照表所列的指数递减和调和递减的线性关系，画在相应的坐标纸上，若能得到一条直线，就表明它符合于那一种递减类型。反之, 若不成直线，它必然属于其它的递减类型。指数递减的半对数直线关系可写为:logQ=A Bt式中：A=logQi,或 Qi=10AB=a2303,或 a=2.303B调和递减的半对数直

14、线关系可写为:logQ=A BNp式中：A=logQi,或 Qi=10AB=a/2.303Qi,或 ai=2.303BQi当由图解法判定递减类型用线性回归法，确定直线的关系数，并由直线的截距和之后，需要利图6-5调和递减类型的产量与累积产量关系截距、斜率和相斜率确定Qi、a或ai的数值。此时，即可建立实用的相关经验公式。2 .试凑法试凑法又称为试差法(tryanderrormethod),它是处理矿场资料常用的一种方法。当用图解法已经确认不是指n值的大小。应数递减时，即可用此法，以判断到底是双曲线递减、调和递减还是产量衰减。主要的判断指标就是 用试凑法的主要关系式为：Qi/Q (t) n=1+

15、nait (6-43)若设：a=1,b=nai,或 ai=b/n则得：Qi/Q (t) n=a+bt假定不同的n值，求出Qi/Q (t)的不 与t的相应的值在直角坐标纸上作图，如图 的是一条直线，直线的斜率为nai,因而可求程中，假如值取得比正确的值大时，则成一条之，当取得的n值比正确的n值小时，则得到当由试凑法得到一条最佳直线，并确定 回归法求得该直线的截距和斜率，并求得ain信正通外值偏小n图6-6利用试凑法求解的关系图拜债偏大同数值，然后将此数据 6-6,正确的n值得到 得值ai。在试凑计算过 向上弯曲的曲线。反 一条向下弯曲的曲线。值之后，即可利用线性 的值。3 .曲线位移法由左向右位

16、移某一合适的距离，使其成为曲线位移法就是将画在双对数坐标纸上成曲线的产量与生产时间图,一条直线的方法。其原理为，将（6-19）式取常用对数后得:log Q1 .将（6-44）式改写为:log Q a blog t c (6-45)log Qi log 1 na(6-44) n式中：a log QiC1/n，或Qj 10aC 1/n图6-7曲线位移法求解的关系图b=1/n,或 n=1/b,c=1/ （ nai）由式（6-45）可以看出，某一正确的C值，可以使 Q与（t+c）的对应数值，在双对数坐标纸上得到一条直线，在给定的C值比正确的C值偏小时，所得到的仍是一条向右弯曲的曲线；反之，如果给定的C

17、值比正确的C值扁大时，则是一条向左弯曲的曲线（见图）当经过曲线的位移，得到一条直线之后，仍然按照（6-45）式进行线性回归，求得直线的截距，并由此确定Qi、n和ai的数值，以满足建立相关经验公式的需要。4 .典型曲线拟合（curve巾t）法将产量递减的公式改成为如下的无量纲形式:a：tQi /Q e i （指数递减）（6-46）" i 1 /n Qi /Q1 nait（双曲线递减）（6-47）Qi /Q 1 a（调和递减）（6-48） I/V利用（6-47）式和（6-48）式，当给定不同的 n值和ait值时，可以计算出不同的产量比（Qi/Q）。然后，将不同n值下的Qi/Q与ait的对

18、应值，画在双对数坐标纸上，即可得到理论的典型曲线图（见图6-8）。若将递减阶段的产量比 Qi/Q与相应的生产时间，画在与典型曲线比例尺相同，并放在典型曲线图上的透明纸 上。然后，在保持画有数据点的透明纸图的坐标，与其典型曲线图的坐标完全重合的条件下，水平向右滑动透明纸 图，使透明纸图上的数据点，能与某一条典型曲线达到最佳拟合为止。在达到最佳拟合之后，可在典型曲线图上直 接读得用以判断递减类型的n值，并可在取任一个Qi/Q值的条件下在典型曲线图的纵坐标上作一水平线，交于最佳拟合的那条典型曲线之后，再往下作垂线，交于典型曲线图的横坐标上得at的数值。最后，再将已得到的ait值除以与Qi/Q值相应的

19、t值，即得ai值。图6-8典型曲线图、应用举例1 .指数递减规律的应用油田产量递减规律分析和应用中的两种常见方法，一种是解析方法，另一种是标准曲线对比方法或拟合法。解析方法的步骤是首先在单对数坐标中绘出如图6-9所示的产量随时间变化的曲线。找出其直线段，求得其初始产量和递减期起算时间。为了计算递减率a可以先从图上计算其周期或半周期，然后按公式（6-9）或（6-10）计算其递减率。这样，产量变化公式就完全确定了。12臬生产黄（百万此）T5)图69产量与时间变化曲线解析方法的另一种作法是在直角坐标中绘制如图图610产量与累积产量变化曲线6-10所示的产量随累积产量变化的曲线。找出其斜率a和递减Np

20、时，其产量变化为初期的产量 Q和累积产量。然后就可以应用递减公式计算产量或累积产量变化。例如若取4 Q,则a=A Q/ Np实例6-1 :某油田一个开发区的产量变化数据由表6-3给出，从表6-3可知该油田产量是递减的，但是处于稳定递减期是1995年以后。试求该油藏的递减率，预测 2005、2006和2007年时的产量和累积产量，并预测其最终采 收率，已知其地质储量为2630 X 104t。表6-3某油田产量变化数据时间（a）19911992199319941995199619971998199920002001200520062007年产量104t/a85.0083.2081.3078.847

21、5.0061.0050.5041.3533.1027.3022.8010.158.316.80累积产量104t11 1卜11 11070.50126.70172.80210.60241.50267.10324.25333.45340.98总累积产量104t150.50234.40316.50396.60473.50508.00600.2646.30684.10715.00740.60797.75806.95814.48解：把日产量随时间变化曲线绘在单对数坐标中（见图6-11）。由图6-11可以看出，油田产量由 1995a开始,处于稳定递减阶段，从此以后，产量变化呈直线，因而可求出其半周期为Ti

22、=3.467a。代入公式（6-11）可得a=0.69315/3.467=0.20 (a-1)2.50数小的量油产年2.001.501.000.500.00时间（a）图6-11某油田产量变化曲线另由表6-3可以获得该油田的产量递减率a为75.00 33.10 八彳” / “a 0.199 (a-1)684.10 473.50用这两种方法求出的递减率在数值上有一定差别，但二者是接近的。代入产量变化公式可以得到 ;IQ(t)=Qie-at=75.00e-0.2(t-1995)Qi /0.2 1 19950.2 t 19954或 N po 1 e375 1 e(104t)a在公式中我们取递减的起始时间

23、t0=1995a,此时的产量 Qi=75X 104t/a,而累积产量 Npo=473.5 X M 由此得出不同时刻的产量及累积产量的预测值如下（表6-4)。表6-4某油田产量和累积产量预测表时间（a）200520062007产量，104t/a10.158.316.80累积产量，104t324.25333.45340.98总累积产量，104 t797.75806.95814.48由公式可以计算递减期最大累积产量为Qi/a=75/0.2=375 (104t)由此可得该开发区可采储量为Npmax=473.5+375=848.5 (104t)精心整理而最终采收率为上面叙述了解析方法的内容及应用，下面再

24、简要叙述运用标准曲线拟合确定递减规律的方法，为此首先要绘出 如图6-12所示的标准曲线图版，在其上分别对不同的递减率值绘出了无因次产量和无因次累积产量变化的曲线。对于一个实际油田， 先在与lgQt同样比例的坐标图上绘出其实际的产量变化曲线和累积产量（由递减初始时刻起算）变化曲线。把这张图与标准图版叠合在一起，然后上下平移，使实际产量变化的直线段与图版上的某一 条标准曲线重合，即可确定出该产量变化的递减率a。再根据直线的左端点，可获得初始产量Qi。再从这一初始时刻起，把累积产量变化曲线与图版上的某一标准曲线拟合，则在标准图版上纵坐标为1.0的水平线在实际曲线上所代表和指出的纵坐标（累积产量）就是

25、递减期可采储量，即Qi/a,而此标准曲线所指的递减率即a为该油田的递减率。实例6-2：仍以上面的实例为基础，说明这一方法的应用。 此时以1995年为起始点。将实例6-1的数据整理后，列于表6-5中。其初始产量 Qi=75X 104t/a,累积产量Npo=473.5 X 10试求该油藏的递减率及可采储量。表6-5某油田产量数据时间1996 年1997 年1998 年1999 年2000 年厂量，104t/d61.0050.5041.3533.1027.30累积产量，104t70.50126.70172.80210.60241.50解：将日产量和累积产量变化曲线在单对数坐标中绘出。把日产量变化曲线

26、和标准曲线对比，可见该曲线a=0.20。这样可得递减期可采储量为75X 104/0.2=375 （ 104t）。这一数值和前解析法得到的结果是相同的。,I2 .双曲线递减规律的应用实例6-3:某溶解气驱油田在生产后期的产量数据如下表6-6和图6-13所示。若油田地质储量为7230X 104t,求产量的初始递减率、递减指数和最终采收率。表6-6某溶解气驱油田的动态数据序号时间年产量，104t/a累积产量，104t11992 年96.33727.0021993 年76.61813.1031994 年60.93881.5941995 年48.46936.0351996 年39.50979.80619

27、97 年32.501015.1071998 年27.271044.8981999 年23.091069.9992000 年19.701091.33102001 年16.931109.61o o o O0 8 6 4N数对的量油产年O21.001990199219941996199820002002时间（a）图6-13油田产量变化曲线解：根据表6-6中数据作图6-13,将图6-13的产量变化曲线与标准曲线图版（其递减指数n=0.3）对比发现该生产曲线与图版上递减率aF0.18的标准曲线能很好地重合在一起，标准曲线原点在实测曲线上所指出的产量即为递减初始产量 Qi,而所指的时间即为递减起始时间，即图 6-13上的点6的坐标。由此得：Qi=32.50X 104t; t0=1997a而由标准曲线可得a=0.