油藏工程常用计算方法

1、. 油藏工程常用计算方法 ;. . 目 录 . 3.1、地层压降对气井绝对无阻流量的影响及预测 . 3 2、利用指数式和二项式确定气井无阻流量差异性研究 . 3.3、预测塔河油田油井产能的方法 . 4.、确定气井高速湍流系数相关经验公式 4 . 4、表皮系数分解5 . 5.、动态预测油藏地质储量方法简介 6 物质平衡法计算地质储量 . 56.1 水驱曲线法计算地质储量 . 76.2 产量递减法计算地质储量 . 86.3 旋回模型预测可采储量. 9 6.4Weng. 试井法计算地质储量 . 106.5 . 15 .7、油井二项式的推导及新型IPR方程的建立 . 15 .8、预测凝析气藏可采储量的

2、方法 . 169、水驱曲线 甲型水驱特征曲线 . 169.1 乙型水驱特征曲线 . 179.2 . 17.10、岩石压缩系数计算方法 . 18.11、地层压力及流压的确定 利用流压计算地层压力. 19. 11.1 利用井口油压计算井底流压. 19. 11.2 利用井口套压计算井底流压. 20. 11.3 利用复压计算平均地层压力的方法（压恢）. 22 11.4 地层压力计算方法的筛选. 22 11.5. . 23A递减分析12、RPS . 24、模型预测方法的原理13 . 25.、采收率计算的公式和方法 14 . 25 、天然水侵量的计算方法.15 稳定流法 . 2715.1 非稳定流法 .

3、2715.2 . 34.、注水替油井动态预测方法研究 16 . 38. 17、确定缝洞单元油水界面方法的探讨 ;. . 、地层压降对气井绝对无阻流量的影响及预测 1*，和其相应的绝对无阻流量就可以用下式计如果知道了气藏的原始地层压力pqAOFi?2* 下的气井绝对无阻流量：。算不同压力pppq?qRiAOFAOFR差异性研究、 2利用指数式和二项式确定气井无阻流量且随着无阻流量的增大两者差指数式确定的无阻流量大于二项式确定的无阻流量， 别越明显。当无阻流量小于50万时，两者相差不大。 、预测塔河油田油井产能的方法3FEp? b） （流压为0油井的绝对无阻流量：。?Jpq? oAOFR252.?

4、Kh.5430 ?JJ 采油指数，；平均地层压力（关井静压），MPa；p ooR?r?e?lnB? ooS?.5r?0?wmS870.?1FE?流动效率，； FE p?pwfRBq2.12oo?m 。?Kho?2pq?cdp? ：油嘴产量公式一（类达西定理推导）hot2 ：油嘴产量公式二（管流推导）ppq?ad?hot2dbptq? +经验）：油嘴产量公式三（试验 o0.5GOR33pp。参/m；mmGOR气油比，m油嘴，；回压，MPa油压，；MPadht数c，a和b可以通过拟合得到。 ;. . 、确定气井高速湍流系数相关经验公式 43878?1.11? 数据回归：K4706.?10?1-1?

5、 。气体高速湍流系数，m渗透率，KmD212B.559?10rh3w? 理论推导：?ZTgB二项式直线关系的斜率，T气藏温度，K。 、表皮系数分解 51?C?32 打开程度表皮系数S1：1h?0.339lnC2.498?3.884C1?.586C948?S2.D1CKhH?hhh。打开程度，C ；pDKrVwhr?perfperfb1?abrL?S10 射孔表皮系数，S2：，K?rKK?LK1DD2HDDVHVhL2perfperfbrlg?aa?ra?b?b ，2D21D1 ;. . ?4?B.71?102oo高速湍流表皮系数S3：。 （油）q?Dq?S o3o0.3878?hrKwo?p4

6、59.scgi。 （气）q?SDq go338780.?KhTrwpgisc钻完井污染表皮系数S4：S4StS1S2S3。 储量方法简介 、动态预测油藏地质6目前，国内外广泛用于油藏地质储量评价的动态法有物质平衡法、水驱曲线法、产量递减法、旋回模型法及试井法（压力恢复曲线拟合法、压降法和压力恢复法、试井综合法）等。探测半径法是试井法中运用的主要方法，物质平衡法是常用的地质储量动态计算方法，这些方法适用于不同的油藏地质和开发条件。下面分别简单介绍： 物质平衡法计算地质储量 6.1如果知道原始地层压力和累计采出量，试井中测到了目前地层压力，或者测试到了阶段压降和阶段采出量，就可以使用这种方法计算储

7、量。 这是物质平衡定律最直接体现。其实，在试井计算储量的其它方法中都遵循这个定律，只是表现的形式不同罢了。 油藏按驱动能量可划分为不同驱动类型。不管哪种驱动类型的油藏中的原始流体的总量必然遵守物质守恒的原则，其主要用途为：根据开发过程中的实际动态资料和流体物性资料预测各种类型油气藏的地质储量，预测油藏天然水侵量，开发过程中定产条件下的压力变化以及油藏最终采收率。 NWW表示累积注水），表示累积产水量（以下以104t表示累积产油量（104t）ppiBW分别为目前地层条件下水、原油及、104t）量（104t），表示水侵量（BBgeowBBS原油及天然气体积系数，、分别为原始地层条件下水、天然气体积

8、系数，、BwiowgiiiRRR分别表示生产油气比和溶解油气比及原始溶解油气比，、表示束缚水饱和度，、pssiCB?BB?B?(R?R)B，假定原始两相体积系数原油两相体积系数，和ftoiigotpsCS?Cfwwi?CC? ot1?SG表示气顶区天然，分别为岩石压缩系数和综合压缩系数，1/MPawi?P表示地层压降，MPa气地面体积，。 （1）未饱和油藏的物质平衡法计算储量 A封闭型弹性驱动油藏 ;. . NBoP?N PC?B 1） （ 地质储量为： （104t） toi B天然水驱和人工注水的弹性水压驱动油藏B)W?W?NBW?(wPeoip ） （2 地质储量为： （104t） ?N

9、PBC?toi ）饱和油藏物质平衡法计算储量（2 溶解气驱油藏AB?R)B?(RNgsiPtp4 ）（地质储量为：103t） （?N C?CSfwwiP?()?B)?B(B tititS?1wi 气顶气和溶解气驱动油藏B 地质储量为：)B?(RRNBgsiPtp4 4） （10（t?N C?CSmBfwwitiP?(B?B)?(1?mB(B?)?)B) tigitgtiS?B1wigiGBgi?m NB 为气顶区天然气气地下体积与含油区原油地下体积之比。oi C溶解气驱和人工注水驱动动油藏 地质储量为：B)W?W)B?(?NB?(RRWtisigPpp?N C?CSfwwiP?()?B)?B(

10、B tititS?1 5） （ （104t） wi 特别地，对于弹性水压驱动油藏，计算步骤如下： 公式法求解动态储量C-1 首先计算弹性产率。 对于封闭的未饱和油藏，即无边水时，B)W?WNB?(wipPo?N P?BC 地质储量为，或写为：toiPK?NB)?WWB?CP?(?BN （）61iwoptoipK，弹性产率可以衡量油田弹性能量的为弹性产率，单位地层压降下的产量，m31大小。采出液体的体积与注入水的体积之差（即地下亏空）与总压降成直线关系，称为 假想压降线，根据其斜率可以求出弹性产率。进而计算边水的水侵量。;. . 不封闭的未饱和油藏，如果有边水存在，并能弥补一定的地下亏空时，压降

11、与亏空曲线不是一条直线。并不是产率增大，而是边水入侵的影响。 边水水侵量为： （7） PW)B?K?WNB?(W?1ppoiweC-2图解法求解动态储量： NB?WpoP与作关系曲线图，利用成直线测点的斜率直线与纵轴截距为动态储Np P?BCtoi 量。进一步利用弹性产率求解动态储量，如果在开发初期，边水入侵速度小，甚至可以 忽略，则弹性水压驱动的方程式可简化为；PPNWNB?CB?K? （8）oioptp在图上如果没有边水入侵，随着亏空体积的不断增加，对于一个封闭的油藏其地下亏空体积与压降之间是直线关系。在座标原点引出的实际亏空切线，称这条直线为假想K。水侵量公式为：压降线。根据其斜率可求出

12、其弹性产率 W?NB?W?K?P （9） pope弹性水压驱动油藏的动态储量： NB?W?(W?W)BwoipeP （10） ?N P?BCtoi物质平衡法适用性条件：该方法不仅适用于均质油藏储量计算，也适用于非均质油藏储量计算，尤其对裂缝性油藏计算储量精度较其它动态方法高。但是前提是一是必须求准地层压力，二是地层压力变化要波及到整个油藏，三是保证累计产量是由同一油藏供给。 水驱曲线法计算地质储量 6.2水驱曲线法适合于高含水油田开发中后期计算原始地质储量。 普通水驱曲线有甲、乙、丙、丁型四种，开发人员比较了解，现在介绍新型水驱曲线，利用新型水驱曲线基本关系式： log?a?bN （11） p

13、p44N?t。为累积产油量，10式中，10为累积产液量， t，pp?和累左右），油田的累积产液量当水驱开发油田进入中期含水之后（含水率40%pN由下式，积产油量b经线性回归求得直线的斜率后在半对数坐标纸上呈直线关系。，p;. . mSoiow （12） 测算油田的地质储量：?N b.3032mS 原始含油饱和度，分数。为油水相对渗透率常数，owoi适用条件：水驱曲线应用的条件，从曲线来看是出现直线段，从油田生产看，则要无重大调整措施时，无重大调整措施，在油田生产是否相对稳定，求生产保持相对稳定，虽未出现直线段也可以校正后使用。水驱曲线不但可以以油田、油藏为单元使用，也可 以一单井或某些井组合使

14、用。 产量递减法计算地质储量6.3油田开发实践表明，无论何种储层类型、驱动方式，以及采用什么开发方式开发的油田，在其开发全过程中，产量一般要经历逐步上升、相对稳定和逐渐下降三个阶段，构成油田的开发模式图。油田何时进入产量递减阶段，主要取决于油藏的储层类型、驱动类型、稳产阶段的采出程度，以及开发调整和提高采油工艺技术和效果。统计表明，左右时，就开始了产量递减阶段。递减方式一般有60%水驱开发油藏当采出可采储量的 其中指数递减主要适用于弹性驱动和重力驱动的油藏。指数递减、双曲递减和调和递减。 ）指数递减（1t 某一开发阶段时间（年）产量与累积产量如下：Dt?eqq? ）（13 iqq?i 14）

15、（?NN? pipD4tqq, 。t分别为时间）（年）时产量和参考时间（年）原油产量（10iN,N t（年）时累积产量和截止参考时间（年）累积产量。分别为pipDt 。和分别为时间（可以用年或月等）和递减率（常数，小数）轴截Y14）式可看出，年产量与累积产量关系在直角坐标上存在一条直线，在由（以后，直线段发生弯曲，递减率变小。因此应用距即为可采储量。一般含水率高达90%以上公式预测有一定的偏差，略低于实际可采储量，预测的参数只能作为可采储量的下 限。 2（）双曲递减其递减类型称之当油田产量随时间的关系曲线在直角坐标上呈双曲线形态变化时， 为双曲线递减。主要适用于各种水驱油田，递减速度比指数递减

16、要缓慢一些。;. . qi?q （15） Dn(1?t) n1n?1qn1n? n i （16） )q)(q?N?N?(n n ipip1Dn?式中递减指数。 ?n?1?（3）调和递减 调和递减是递减指数时特定条件下的递减类型。产量递减速度低于双曲线递减1n?的速度，适用于递减阶段的后期。 ?1)?Dt?q(1q （17） iDNP?qlogq?log （18） i2.303qi适用性条件：与水驱法一样，适用于油田开发后期无重大调整措施，尤其是用调和法预测极限可采储量时，如何确定极限采油量是一个很难解决的问题。 6.4Weng旋回模型预测可采储量 对于资源有限体系，其初期、中期和后期开采的全过

17、程可以用翁氏旋回模型表述。y?yn?t00e?AtQ、 油田年产量： （19） ?t CQ为油气田年产量，104t/年；t为翁氏时间，年； y油气田某一生产年份；y0油气田某一生产参考年份 n、C为模型常数 1n?B?AC令、 （21） ?a C?(n?1)4)1(n?N?Ba?t ） （1022）则得到可采储量： R式中为Gamma函数。 )(n1?n?t?)log(Qt? （23） 经过简化则有关系：?logB （24） ?0.4343a?n?ntlog(Qt?)?值，直至获得较好的直线进行求解，选择合适的将产量数据按?a和B，最后得到段，对直线进行回归后，得到直线截距和斜率，由此求得We

18、ng旋;. . 回模型计算的油田可采储量。 适用性条件：该方法也是适用于油田开发所有时期，但是必须是中间过程中没有增加开发井数和开发层位等。 6.5试井法计算地质储量 -此部分为项目中摘出部分 动态储量是根据井生产取得的测试信息（如产量、累积采出量和压力数据）计算出来的。它与许多计算静态储量的物理参数（如面积、厚度、孔隙度等参数）没有直接关系。这就是说，一组合格的测试信息所对应的储量是确定的。当使用解释模型去拟合这组信息时，改变厚度、或者孔隙度时，计算出来的储量总是确定的。例如，当厚度变小时，模拟系统就会增大面积来弥补。使用科学的语言的来描述，它遵循物质平衡定律。从储层中索取了多少，就有一个对

19、应的信息即压力作出相应的反映，它是唯一的。就试井计算储量的方法来说，储层的孔隙空间是由压力和产量信息所唯一确定的。除了测试信息以外，影响试井计算储量的参数是含油饱和度和综合压缩系数，故需要认真仔细确定。对于可动油饱和度变为零的区域，即进入到边底水范围，这部分储量在计算时应该排除。所以对于动态储量，特别是试井计算的储量都要注意是否含有水储量，这个水储量是多大，必须认真对待，需要想办法排除。 6.5.1 压降曲线法计算储量 在试井计算储量的方法中，首推该种方法。因为，这种方法在理论上是精确的，我们通过油藏数值模拟的方法对该种方法进行过验证，计算误差小于1%。 当油气井以稳定产量开井生产，所测试的井

20、底流动压力随开井生产时间的关系曲线，称为压降曲线。其按压力随时间的动态变化，可以划分为非稳定阶段、过渡阶段和t分成三个阶段界限，非稳定阶段与过拟稳定阶段。对于封闭油气藏可以无量纲时间Dett=0.25；过渡阶段与拟稳定阶段的界限为渡阶段的界限为。无量纲时间关系式=0.1DeDe?3Kt6?103.?t 为： De2?rCet2?t?m为地层有效孔隙度，小数；为压降时间，式中：K为地层有效渗透率，hr；rC油藏半径，m。；为综合压缩系数，smpa为粘度，?；1/MPa et在下面过程中用到探测半径方法。这种方法是指试井过程中没有测到任何边界，可;. . Kt12?0.r 。以通过计算探测半径来计

21、算探测范围中的储量。探测半径为：（m）i?Ct探测半径是指井不稳定过程在没有遇到不渗透边界或定压边界的情形下所影响到 的最大范围。在探测半径处压力变化为零。 、非稳定阶段1Kt12?r0.计算探测范围地质储量。在非稳定非稳定阶段根据探测半径关系式i?Ctm可求得地层有效渗）（MPa/cycle阶段根据井底流动压力与时间半对数关系曲线的斜率?B.12q2 透率：?K mh 综合压缩系数为：C*t?CCCS?CCS? ， tofoitowiwSoitqp1920t0959q0.igo油气藏工程实，对于气体：得到控制地质储量为：（?N?G *CmmCtto 。页）用方法416 q 。式中为压降期间稳

22、定日产油量，t/do 、过渡阶段2过渡阶段，又称为非稳定阶段的晚期，它存在于非稳定阶段和拟稳定阶段之间，整个地层内的压力变化动态，尚未进入拟稳定条件。该阶段的长短，主要取决于封闭的油气藏大小和地层导压系数的数值。油气藏越小，而地层导压系数越大，则过渡段越短。bt?aplog(p?)b,ap 将测试数据绘成曲线，通过迭代计算可求得值。RwfR2230qpq0.11156igo，该阶段预测油井控制的地质储量为：对于气体：（油?NG? *a*a10bCbC10tto 418气藏工程实用方法页）。 、拟稳定阶段3随着生产时间的延续，该阶段又称为半稳定阶段。当探测半径达到油气藏边界之后，当油气井控制的油

23、气藏范围内任一点的地层其压力动态已偏离无限大作用地层的特征，压力降（包括井底流动压力）达到同步速率下降，油气井的压力动态已由过渡阶段转入 拟稳定阶段，此时井底流动压力随时间变化呈直线关系。mt?p ，可求得直线斜率）1曲线（如图将测试数据绘成* MPa/hr（）wf;. . q0416710.o 油气藏工程实践189页）。（该阶段预测油井控制的地质储量为：?N*Cm*tp.08334q0ig 189页）。对于气体：（油气藏工程实践?G*Cmt 图1 压力降落曲线示意图对于油气藏中多井同时生产情况下，当达到拟稳定流时，可用该方法进行整个油气 藏地质储量计算，但公式中所用产量为全油气藏整个平均日产

24、量。压力降落法适用性条件：压力降落法计算储量适用于定容有限封闭油气藏开发的 早、中、后期。 6.5.2 压力恢复曲线法计算储量q）开井生产，当其压力动态（对于有限封闭的油气藏，油井若以定原油产量t/do达到非稳定阶段，或过渡阶段，或拟稳定阶段之后，将井关闭测试压力恢复曲线，可分 以下三种情况。 关井前的压力动态处于非稳定阶段A.t?0.0539qso 该阶段未受到边界影响，预测油井控制的地质储量为：?N*Cmtoot?PMDH；再将）数值时，可先由Horner在确定时间法外推求得原始地层压力（hrsit?P ，得到相应的时间，即可作为的数值。法的直线外推到si 关井前的压力动态处于过渡阶段B.

25、0.11156qo 该阶段预测油井控制的地质储量为：?N*abC10to;. . log(p?p)?a?bt曲线，通过反复迭代计算方法可求得将测试数据绘成p,a,bwsRRp和p分别为地层压力和关井恢复压力，MPa。 值。wsRC.关井前的压力动态处在拟稳定阶段 0.00763Cq?tpAo 该阶段预测油井控制地质储量：?N *mCto其中为压力恢复半对数曲线斜率（MPa/cycle），为油藏形状因子，无量纲。CmA*?AS8277.Ctooi ?t pKCA 压力恢复法的适用条件与压力降落法相同。 压力曲线拟合法计算储量6.5.3 因为在封闭油藏中一口井关井测恢复没有拟稳这种方法特别适合于压

26、力恢复试井，态，只有稳定态即地层压力达到稳定。压力曲线拟合包括双对数曲线拟合、半对数曲线A，），通过与封闭油藏模型拟合可以确定油藏的面积2-4拟合和压力史曲线拟合（见图 。然后用容积法公式计算储量。?B/SAhN?ooo? 2A 式中。为原油比重，为油藏面积，mo但拟合的面积会随着输入的厚度从上式中看起来和静态计算储量的公式是一样的，或孔隙度变小而变大。即拟合的孔隙体积是恒定的，是由测试信息所确定的，已经包含 在测试信息中。压力曲线适用性条件：适用于封闭有限油气藏，即探测到整个边界或部分边界，对于部分边界情况，除了恢复数据外，还有压力历史数据，即通过整个压力史拟合确定其 它边界。;. . 图2

27、 双对数拟合曲线 图3 半对数拟合曲线图 图4 压力史拟合曲线图 ;. . 6.5.4 试井综合法测试过程中只测到了部分边界，许多情况下需要将试井分析结果与地质资料（包括构造图）结合起来确定含油气面积，再用容积法计算储量。许多试井资料只能用这种方法计算储量，特别是测到油水边界的情况下需要进行综合分析。例如，构造上倾部位测到了两条夹角相交的不渗透断层，低部位测到了油水边界，只能根据这些边界组合的几何形状，计算含油面积，然后使用上式计算储量。只要输入到解释系统中的参数与计算储量的参数是一致的，那么计算的储量是唯一的。变小解释系统的厚度，就会增大边界距离。总之，正如开始所说的，一组测试信息，对应的储

28、量总是唯一的。只有改变含 油饱和度或综合压缩系数，才会改变试井计算的储量。 适用性条件：适用于任何类型油气藏，但要有相应的试井解释模型。方程的建立、油井二项式的推导及新型IPR 7q6o?q 一点法确定的油井无阻流量：AOFpp?wfR1?1?48pR、预测凝析气藏可采储量的方法 8?M?24056GE，GE凝析油气体当量体积M：凝析油的分子量。 ooooo1?f。： 凝析气藏的干气摩尔分量fg gGEo1? GOR?830?GORog?凝析气井总井流相对密度： 。? wtGOR?GEo可用下式计算总井流的拟临界压力和临界温度： 2?07653.?086774.?0.3565?p ）MPa（w

29、twtpc2?72222.?39?T103.8889183.3333 ）（Kwtpcwt偏差系数（见油层物理-何更生106页）：求得拟临界压力和临界温度后就可以求得拟对比压力和拟对比温度，通过查图版求得Z值。或者通过下面的经验公式计算： ;. . 1?ad cp?a?Z? prexp(b)?50. 1010.36TT?0.92?a?1.390prpr6?p.3200660.?pr50.2 p?0.037T?b.0.62?023T? prprpr1T?exp20.T?0.8672?prprc?0.132?0.32lgT pr?20.7153?1.1285Td ?exp?0.4201Tprpr1?

30、。凝析油含量的计算方法： GOR?GEo、水驱曲线 9甲乙丙型水驱曲线的理论推导、地层含水饱和度和水油比关系式的推导、油田含水率和累计产油量的关系式推导、测算可采储量和采收率的关系式推导。 水驱特征曲线是指油田注水（或天然水驱）开发过程中，累积产油、累积产水和累积产液量之间的某种关系曲线。应用于天然水驱和人工注水开发油田的水驱曲线很多，经过多年的实践应用，普遍认为其中四种水驱特征曲线具有较好的实用意义，即甲型、乙型、丙型和丁型水驱特征曲线。本次计算用到了甲型、乙型水驱特征曲线，现列出这两种水驱特征曲线的关系式： 甲型水驱特征曲线 9.1lnW?a?bN 甲型水驱特征曲线：P11P434NWba

31、与曲线有关的系数。 累积产油量，10，累积产水量，10tm；pp11?f1?WN?lg?c 通过推导可得：? p1b1?f?W1fc?a?lnb。含水率。式中 W1111(3.8918?c)?N 当油田极限含水率取0.98可得油田可采储量的计算公式：时， 1Rb1利用甲型水驱特征曲线可以得到油田的地质储量的相关经验公式： ;. . 49690.?t。 ，地质储量，10b.905N?16N1乙型水驱特征曲线 9.2乙型水驱特征曲线： Nb?L?alnPP2243；，m与曲线有关的系数。累积产液量，10 Lbap22?11?lg?N?c 通过推导可得：? 2p1?bf?W2 。式中blnac?22

32、2当油田极限含水率取0.98时，可得油田可采储量的计算公式： 1(3.912?c)N? 2Rb2、岩石压缩系数计算方法10 C的主要方法有：Hall经验公式，目前计算岩石有效压缩系数Newman公式，f李传亮公式。 a、Hall经验公式 1953年，HNHall通过大量的实际测量数据，统计出来岩石压缩系数与孔隙度?4105846?2. 之间的关系曲线，得出经验公式如下：?C f43580.?为岩石孔隙度，百分数值； 式中，b.Newman公式 Newman（1991） 推导如下关系式： I）砂岩类的岩石孔隙度计算公式： ?1.?342859 )8721?55?C14.11510./(1fII）

33、石灰岩类的岩石孔隙度计算公式： ?0.69299 )4766421754429123?C./(?.?10f?为岩石孔隙度，百分数值；其中， ;. . c.李氏推导公式 ?C?C 2003年李传亮提出的岩石压缩系数计算公式如下： Sf?1 其中?岩石孔隙度，百分数值；固体骨架的压缩系数，为一常数。 CS3(1?2v)?C 在弹性变形条件下，固体物质的压缩系数可以用下式计算： sEs式中v固体物质的泊松比，无量纲；固体物质的弹性模量，MPa。 Es由于岩石的弹性模量（E）容易测得，可以通过下式计算岩石的平均弹性模量：EsE ?E s?14E一般在0.3左右，弹性模量岩石固体骨架物质的泊松比一般在M

34、Pa1010)?(1vs?(1?2v3)）。之间。（注：根据李传亮推导公式也可以为： ?C fEC 的试验测定d.f e岩石压缩系数计算方法比较图版经验公式法、李氏法和Hall目前确定岩石压缩系数的方法主要有实验测量法、 经验公式方法。Newman，21Cf三种计算方法比较：鉴于实测岩心资料中，的范围大概在：810-4MPa1时比较接近实510经验公式所得值范围变化较小，孔隙度在以上三种方法中，Hall所得岩石压缩系数明显偏小。而李氏公式中，Newman测值。公式所得值范围变化偏大，李传亮认为Hall图版显示了错误的关系曲线（李传亮编著的油藏工程原理中有详述）. Hall方法比其它方法可靠但实

35、际应用表明种处理方法，一是忽略岩石的压缩系数，由于岩石压缩系数的数值通常较小，有2二是严格计算出或测试出岩石压缩系数以免对动态储量的计算结果造成不必要的影响， 值。 、地层压力及流压的确定11获得地层压力可以通过实际静压来得要进行物质平衡计算必须有可靠的地层压力。到，但目前测试静压井的非常少，这时可以通过复压解释得到，实际油田开发中测试复压的资料也比较少。因此，需要研究如何利用其它相关数据（如流压、油压、产量等数;. . 据）确定目前生产状况下的地层压力，应用井口油套压力计算井底压力，然后估算地层压力。 利用 流压计算地层压力11.1应用现有的流压计算地层压力：利用已测得的流压计算地层压力，计

36、算时要用到产能方程：。 )p?p(q?Jwfr1 在实际应用时也可以变换写成：p?qp? rwfJp?q曲线图确定采油指数和地层压力） （已知流压和日产量也可以通过wfp）和井底流压方法思想：注意到上面给出的产能方程，有两个已知数（日产量qwfp）和地层压力，由于日产量和流压是变化的，要想求出地层两个未知数（采油指数Jr压力必须建立方程组，从而我们把采油指数和地层压力看成是暂时不变的值，为了保证这一点，我们在实际计算的过程中就要注意，选取日产量和对应流压的时候就要尽量找日产量有明显变化而时间间隔小的两个或多个生产段（保证每个生产段内日产量稳定，选取的生产段越多，计算的结果越精确）。 11.2利

37、用井口油压计算井底流压 实际数据中已测得的流压值是有限的，只有部分井有流压值，对于那些没有实测流压值的井来说，就要考虑利用地面上的数据来折算井底流压，考虑到每天都可以得到产量和井口油、套压，所以可以根据油压和套压来计算，下面这种方法首先就利用了井口油压。 方法思想：由于多相垂直管流中每相流体的参数及混合密度和流速等都随压力和温度而变，因而沿程压力梯度并不是常数。因此，多相管流需要分段计算，并预先求得相应段的流体物性参数。然而，这些参数又是压力和温度的函数，压力却又是计算中需要的未知数。所以，多相管流通常采用迭代法进行计算，即深度迭代和压力迭代。 下面是深度迭代的计算步骤： （1）任一点（井口或

38、井底）的压力p作为起点，任选一个合适的压力降作为计p?算压力的间隔。一般选=500-1000 kPa 。具体数值应根据流体流量（油井的气、液p?产量）、管长（井深）及流体性质确定。 （2）估计一个对应于的深度，以便根据温度梯度估算该段下端的温度T。 p?h?l;. . （1）计算出该管段的平均温度及平均压力，并确定在该温度和压力下的全部流体性质。 dp ）计算该段的压力梯度（2 dhdp （3）计算对应于的该段管长为； ?h?p/()p? dh进行比较，若两者之差超过允许范步中估计的步计算的（1）将第5与第2h?h?围，则以计算的作为估计值，重复进行2-5的步骤，直到计算的与估计的之差h?h?

39、内为止。 在允许的范围0n?1i?，3.n，2pi?p?p、Pi：及 （2）计算该段下端对应的深度Li h?L?0iii1i?以Li处的压力为起点，重复第2-7的步骤，计算下一段的深度和压力，直到各段的累加深度等于或大于管长时为止。 对于单相流：当油井的井口压力大于饱和压力的时候，油管内流体就是单相流，这p?p?p?p 时利用垂直管流计算井底流压公式：whHwffrppp井口油压MPa MPa 摩擦阻力式中：MPa 井底流压 frwhwfp和，三个数据，但在实际数据里从公式可以看出要计算流压，就要知道ppfrwhHpp也可以求得。 可以找到的只有。而对于单相管流而言，whH?H?p 公式为：H

40、3? 油管内单相流的密度kg/m 井深m Hp不容易确定，可以利用实测流压求出摩擦阻力。在利用摩擦阻力最后一个参数fr计算本井其他时间的流压。 11.3利用井口套压计算井底流压 井内各相流体分布如图： ;. . 井内各相流体分布 实际资料中流压数据非常少，气柱高度一般没有实测值，需要建立流压、套压及气 柱高度的关系，再利用这三者之间的关系得到更多的流压。气体造成的压力，从井口的套压折算到井底的流压实际分三部分（依次是井口套压， 混合液体造成的压力）。p?p?pp? 根据油气藏气体部分流压计算公式： 可用（2-12）式表示：cwf液气s16?22)e11.2893*10?pFQTZ(mscwga

41、as2?ep?p cgas5D得到油井套管内气体部分流压计算，得到我们只涉及油管生产的情形，上式中=0Qs2e?pp 公式：C气wfs2p?ep?p? 1-26（ ）式可写成：所以，cwf液?h*.003415gs2?p?epp?s)hp?0.00981(H?h?p? ；式中 ；wgcc液气气ZT 式中：?p ；MPa；气柱内气体相对密度，空气为1井口套压， gc；为平均气体偏差系数； 气体区平均温度， KZTh气柱高，m；H井深，m。 ?h)?(H?h?pp0.00981 gwflc?f)?(?f?1 其中owlww;. . 11.4利用复压计算平均地层压力的方法（压恢）Hassan-Kab

42、irArps-Smtih法，通过研究已有的各种地层压力计算方法，优化出了法又分为径向流和线性流两种。法能较好地应用于塔河油田，其中Kuchuk法和Kuchuk 有关方法的原理如下：)0?p(pwfws ：） 1）胡氏法（修正的Hassan-Kabir法）bpp? avgwst?）pws-pwf由上式可以看出，在普通直角坐标纸上，以关井恢复压力为纵坐标，以（的关系图，将得到一条直线，该直线外t）/t为横坐标，作出pws与（pws-pwf/ 时，纵坐标轴上的截距即为油井泄油面积内的平均地层压力。t=0推到（p-p）/wfwsc Hassan-Kabir法：2）?ap? wst?b?p与本公式不受边

43、界条件的影响。作上式表示一条水平渐进线为a的等轴双曲线，ws可得到正确值。用试错法假定不同的b，b和c1/（b+t）的关系图可线性回归得到aap? 的直线段。在无限大的关井时间里，井底压力趋近于平均地层压力。avg Kuchuk法（径向流）3）b 对于径向流有：?pp(?t) avgwst?pp 值。和1/t的关系图可线性回归得到-b作与avgws 法（线性流）Kuchuk4b?p(?t)p 对于线性流有：avgwst?pp和-b值。 ）1/（t0.5的关系图可线性回归得到作与avgws11.5地层压力计算方法的筛选 在以上几种压力计算可以分为两部分：一种方法是利用有关资料直接计算地层压力；另

44、一种方法是要根据地面上的数据先把井底流压计算出来进一步计算得到地层压力。要把地层压力计算的很准确，就需要实测的数据很精确，总结这上述几种方法，得出以下结论： （1）只要实测的流压数据符合是在产量平稳的时间段测得的，对应每个流压数据的平均日产量要和其他数据的有明显变化，保证每个时间段应该至少平稳1个月且间隔不能太长（一般为10天左右），那么利用已测得的流压计算得到的地层压力就很准确了。 （2）利用油压计算流压一般在计算单井单相流的时候比较准确。 （3）利用套压计算流压，需要实测的动液面数据，首先把流压确定的很准确再进一步推地层压力的时候也就很准确了。 ;. . （4）利用复压计算平均地层压力的方

45、法是将关井后测得的一些数据导入事先编好的软件里，得到一个压力恢复曲线，从而将地层压力回归得到，此方法外推地层压力比较准确。 、递减分析 Arps12根据递减指数取值的不同，在递减分析中包括三种递减类型，即：指数递减、双曲线递减、调和递减。这三种递减类型的主要关系式如下表所示： 递减类型 指双调 和 线 数 曲 递减指数n?00?n?1n?1 产量与时?1n?1/Dt?)?Dtt)?q(1q()t?nD(q(t)?q1eq(t)?q 间iiiiinEqqEq)qE(q?1产量与累iin?1lnN?n?1ii?qN?()q?N 积产量pippqDnD1?Diin递减指数；D瞬时递减率，%/mon或%/a； 44t/at/mon递减阶段的初始产量，q； 或1010i4Nt；自递减阶段初始时刻的累计产量，10或t时间的产量，mona；