致密砂岩气藏水锁损害的孔隙分形预测模型.docxVIP

致密砂岩气藏水锁损害的孔隙分形预测模型 在钻井、采矿、生产和天然气生产的过程中，水基工作液入侵造成的水锁破坏是致密砂岩气藏最重要的破坏类型。水锁损害一旦发生则难以完全解除,会严重影响气藏的发现、储层的评价和后期的开发,因此对储层水锁损害程度进行准确的预测尤为重要。迄今为止,国内外学者建立了多种水锁损害的预测方法,包括实验数据拟合预测、灰色预测法、神经网络预测、相圈闭系数法,APTi指数法等,这些方法主要是根据实验数据的回归,或者基于孔隙度和渗透率等参数建立计算模型,而水锁损害受到储层岩性、物性、孔隙结构、侵入流体等多种因素的控制,因此目前的预测方法并不能较好的反映水锁损害程度。采用泡压法测试得到储层岩石的孔隙结构,能够反应不同含水饱和度下的渗流特征,在此基础上进行分形特征研究,建立了水锁损害程度预测新方法。 1 泡点压力法测试储层岩石的连通喉道对渗流能力的影响 常规的压汞法是测试储层岩石孔隙结构的最常用方法,但是压汞法测试了包括连通孔隙、盲孔在内的几乎所有岩石孔隙的体积分布特征,并以此计算不同孔喉对渗流的贡献,而决定储层渗流能力的关键在于连通喉道的大小和分布,可见压汞法测试的孔径分布并不能直接反应储层的渗透特征。利用孔隙内的两相平衡和渗透性质为基础的泡点压力法以气体和现场工作流体作为测试介质,属于无损检测技术,可测定储层岩石的连通喉道对渗流能力的分布。 泡点压力法测试时首先将样品饱和润湿流体,以气体在样品一侧施加一个从零开始逐渐增大的压力,孔喉由大到小逐次被打开,再测试相同压力范围内的非润湿条件下的流量变化特征,假定一定压力下气体流量与被打开连通喉道的面积成正比,流量在润湿条件下与非润湿条件下之比(R)反映被打开孔隙面积的分率: 式中Fw为样品在润湿条件下被一定压力打开连通喉道时的流量;FG为样品在非润湿条件下与Fw相同压力时的流量,可见该比率为一定压力下被打开连通喉道对渗流能力的贡献率。将R对喉道直径(r)微分,得到连通喉道对渗流能力贡献分布函数: 测得不同压力润湿样品和非润湿样品的流量数据后,根据以上两式可以得出样品不同有效喉道对渗流贡献率的分布。 2 流量q的密度函数 储层岩石的孔隙结构具有强烈的不规则性,受到多种因素的控制,使得经典几何学为基础的孔隙结构模型与真实的孔隙相差太远。分形理论是研究不规则形体的自相似性及其复杂程度的理论,分形维数作为描述分形特征的重要参数,为定量描述储层孔隙结构提供了有效一条有效的途径。 根据分形的几何原理,若储层孔径分布具有分形特征,则储层中直径大于r的孔隙数量(N)与r有如下幂函数关系: 式中rmax为储层中最大孔隙直径;f(r)为孔径分布密度函数;a为比例常数;D为孔隙分形维数。将式(3)对r求导得: 根据泡点压力法测试原理,实验过程中当压差升高到p时,喉道从小到大逐次被打开,对应的流量Q可以由Hagen-poiseuille方程的积分形式表示: 其中μ为液体黏度;l为岩心厚度,rmax为最大连通喉道直径;r为压差为Δp时打开的最小连通喉道直径;n为被打开的直径为r的连通喉道数。 假设多孔介质孔的总数为N,公式(5)中n可以表示为: 式中f(r)为连通喉道分布函数,将公式(6)代入公式(5)得: 将连通喉道分布函数(4)代入上式可以得到: 代入Laplace方程得: 积分得: 设,为常数,则有: 该式即为在压差Δp下通过直径大于r的连通喉道流量大小。同理,所有连通喉道被打开后(即非润湿条件下),相同压力下通过整个岩心孔喉的流量为: 则公式(11)和公式(12)之比为压差Δp下通过直径大于r的孔喉流量占过流面积流量比率W为: 相对应,在未饱和润湿相时,压差Δp下通过直径小于r的孔喉流量所占过流面积比例为: 当rmin趋于0,且远小于rmax和r时,上式简化为: 两边取对数: 将Laplace方程代入公式(15),可以得出驱替压力与渗流比率关系式: 利用泡压法分别测试润湿条件和非润湿条件下不同压差时的流量数据,根据上述公式计算不同直径连通喉道对渗流能力的比率W′,将lgW′与lgr或者lgp作图,结合公式(16)和公式(17),通过线性回归得出斜率,即可得出分形维数(D)。 3 储层水锁损害特征 水锁的本质是由于水相占据了不同的渗流空间,导致气相渗流能力降低。 运用泡压法测试得到的孔隙结构计算获取的分形维数进行水锁预测就有以下几方面的优点:(1)计算方法基于连通喉道的渗流分布,本身就反映了储层岩石在不同含水饱和度下的渗流特征;(2)气体泡压法测试连通喉道渗流分布时,首先饱和润湿介质(可以是地层水或者实际工作流体),再由大到小逐次打开孔喉,这一过程与实际储层遭到水相侵入后的气液流动过程相符合;(3)由于使用了储层岩心和实际工作流体进行实验测试,因而该分形维数包含了岩石物性、孔隙结构、黏土矿