储量计算方法的原理.docx

弹性二相法是气藏开发早期或试采阶段的一种重要储量评价方法。对于其推导过程，目前相关书籍和文献均直接从拟稳态阶段Pwf2与t的直线关系式出发，然后根据该直线斜率的表达式推导出储量计算式。本章对弹性二相法进行完整推导时发现，所谓的Pwf2与t直线关系是基于两个假设前提而建立的；此外，本章还研究了气藏采出程度对该法精度的影响规律。2.1弹性二相法的原理及误差分析2.1.1 弹性二相法的原理分析当气井以恒定产量生产，并进入拟稳态后，任一时刻t的产能方程为：(2.1)式中，为t时刻的平均地层压力，MPa；Pwf为t时刻的井底流压，MPa；q为转换成标准状态下的日产气量，m3/d；μ为t时刻井底流压和平均地层压力的平均值所对应的气体粘度，mPa·s；Z为t时刻井底流压和平均地层压力的平均值所对应的偏差因子[20]；T为t时刻的储层温度，K；k为储层的径向渗透率，mD；h为储层的有效厚度，m；re为泄气区域的外边界半径，m；rw为井筒半径，m；S为表皮；D为紊流系数，d/m3。当气藏在较短时间内达到拟稳态，假设气体、岩石和束缚水的压缩性在短期内可忽略，则有： (2.2)式中，Pi为原始地层压力，MPa；Gi为原始地质储量转换为地面标准条件下的体积，m3；t为从投产到目前的累计生产时间，d。结合式（2.1）和（2.2）可消去平均地层压力得出： (2.3)整理可得： (2.4)假设：(2.5)则式（2.4）可简化为：(2.6)令 (2.7)(2.8)则式（2.6）可简化为：(2.9)这表明当气藏进入拟稳态时，井底流压平方与时间呈直线关系，如下图2.1所示。图2.1 气井定产生产时Pwf2与t关系曲线示意图[14]Fig.2.1Schematic plot ofPwf2vs.t for a gas well producing at a constant rate根据生产数据拟合得出β后，利用下式可求储量。 (2.10)1.2.2 流动物质平衡法传统的物质平衡法又被称为压降法，是利用气藏平均压力（P/Z）和累计产气量（GP）之间的关系直线在GP坐标上的截距确定原始地质储量的[20]。该方法要求多次关井测定平均储层压力，这一要求在低渗气藏中通常难以满足。因此，依赖气藏平均压力数据的压降法在低渗气藏开发中不甚适用。1995年，Mattar首次提出气藏流动物质平衡法，该方法要求气藏进入边界控制流，并以恒定产量生产一定的时间，此时井底流压Pwf/Zwf与累计产气量GP之间呈线性关系，并且该直线平行于平均地层压力与累计产气量之间的直线（P/Z~GP）。因此过原始地层压力点（Pi/Zi）作Pwf/Zwf~GP直线的平行线，利用该平行线在GP坐标上的截距即可确定原始地质储量。该方法解决了低渗气藏不关井测压给动态储量评价带来的困难，在国内外得到广泛应用[21-24]。但是该法要求气井达到边界控制流阶段，并且产量保持恒定，这使该法在变产量生产或者产量波动的情况下并不适用；此外，本文在运用该方法时发现Pwf/Zwf~GP直线与P/Z~GP直线并非总是平行的，即流动物质平衡法本身存在误差，其原因有待深入研究。1999年，Ram G. Agarwal和 David C. Gardner等采用物质平衡拟时间，建立了产量~累产量的递减典型曲线（物质平衡拟时间是由Palacio和Blasingame提出的，可有效处理气体PVT参数随压力的变化以及气井流量变化两大问题）。当生产数据与产量~累产量递减典型曲线达到最佳拟合效果时，可准确评价气井控制储量[15]。该法解决了气井变流量生产时的储量评价问题，但是由于计算典型曲线横坐标量时，用到了等效井筒半径参数。使用这个参数恰恰制约了该方法的适用范围和精度，因为在近井污染、气体高速非达西效应或复杂结构井等的影响下，难以准确评价该参数。该参数的误差将引起典型曲线横坐标量的计算误差，进而影响典型曲线评价储量的精度。2003年，L.Mattar等将Agarwal-Gardner的产量-累计产量递减典型曲线法进行了改进，建立了归一化产量与归一化累产之间的线性方程，该直线在横坐标上的截距即为原始地质储量。由于在计算该线性方程的横、纵坐标量时，没有用到等效半径这个参数，因此避免了该参数的误差所引起的储量评价误差，此外还增强了直观显示的效果[3]。L.Mattar改进了该方法，并运用到目前广为使用的Fast·RTA软件中，但仍将其命名为AG FMB（即Agarwal-Gardner流动物质平衡法）。由Mattar经过改进所建立的Agarwal-Gardner变流量流动物质平衡法在发表之后被广泛采用。但是，该法仅在单井开采一个气藏单元时较为可靠，当多井开采一个气藏单元时，往往存在投产时间不一致、部分井关井、修井、产量调整等现象，而这些现象均会引起单