低速非达西渗流动边界移动规律 @@.pdfVIP

东 北 石 油 大 学 学 报 第40卷 第2期 2016年4月 UNIVERSITY V01．40No．2 2016 JOURNALOFNORTHEASTPETROLEUM Apr． DOI 10．3969／j．issn．2095—4107．2016．02．009 低速非达西渗流动边界移动规律 王建俊1，鞠斌山1，罗二辉2 100083；2．中国石油勘探开发研究院，北京 100083) (1．中国地质大学(北京)能源学院．北京 摘要：低渗透油藏非达西渗流动边界的产生归因于启动压力梯度。在考虑动边界的基础上，建立包含表皮效应和 井筒存储的单一、以Warren—Root模型为基础的双重和三重介质低速非达西渗流数学模型；并通过定义无量纲压力得到 模型的Laplace解，利用Stehfest数值反演获得动边界的传播规律，对比动边界与无限大、封闭及定压等三类外边界条件 下的压力及压力导数分布特征。结果表明：对于低速非达西渗流，启动压力梯度越大，动边界向外传播越幔；在双重介质 和三重介质中，动边界移动出现短暂的停滞，曲线上表现为水平段或平缓段；动边界主要影响压力及压力导数曲线的中 后期，而无限大、封闭及定压等三类外边界主要影响曲线的后期形态。 关键词：非达西渗流；启动压力梯度；动边界；单一介质；双重介质；三重介质 中图分类号：TE312 文献标识码：A 文章编号：2095—4107(2016)02—0071一07 0 引言 达西定律已广泛应用于多孔介质的渗流力学问题，当渗流速度过快或过慢时，达西定律描述的渗流速 度与压力梯度之间的线性关系不再成立，偏离达西定律线性关系的渗流问题可以归结为非达西渗流。低 速非达西渗流是指流速较慢而导致偏离达西定律的流动，如低渗透储层或深部弱透水层[1-2]。李凡华 等[33研究表明，当启动压力梯度值大于0时，压力的传播并非瞬时达到无穷远，可以采用数值离散化方法 计算动边界随时间的变化；贾振岐等[4]通过实验证明，低渗透油藏的孔隙越小、喉道越窄、固液表面的分子 作用力越强，启动压力越高；程时清等[5]对低速非达西动边界问题进行积分求解；闫栋栋等[61通过气泡法 岩心测试表明流度越低，启动压力梯度越高；计秉玉等[7]运用稳态逐次逼近方法，求解定流量和定流压条 件下的非达西渗流模型；王晓冬等[8]采用数理方法，求解包含启动压力梯度和动边界的非线性不定常渗流 数学模型；王军磊等[9]利用叠加原理，求解考虑启动压力梯度的七点井网稳态压力方程，表明启动压力梯 度越大，有效动用程度越低。 当储层中存在溶洞、裂缝时，渗流通道不再是单一的孔隙介质，而是形成双重介质和三重介质。多重 介质中低速非达西渗流问题更加复杂。姜瑞忠等[1们利用油藏工程方法，研究裂缝存在时形成的优势渗流 通道。针对双重介质，罗二辉、程林松、XuJianchun等考虑启动压力梯度，得到非达西渗流的特征[1卜13]。 罗二辉等[1们利用解析方法，得到三重介质中压力分布，并分析敏感因素。由于低速非达西渗流具有高度 非线性，难以获得解析解，目前还没有针对不同介质中动边界的移动规律的研究。在考虑动边界的基础 上，基于启动压力梯度的非达西渗流，笔者建立包含表皮效应和井筒存储的单一、以Warren—Root模型为 基础的双重[1 获得压力分布的解析解表达式，分析不同介质类型中启动压力梯度对动边界传播规律的影响，并对比动边 界与无限大、封闭及定压等三类外边界条件下的压力及压力导数的变化动态特征。 1 低速非达西渗流数学模型 1．1假定条件 (1)平面径向流动，储层厚度为h；(2)介质分别为单一、双重和三重孔隙介质且各向同性，外边界为无 收稿日期：2015—11—13；编辑：刘丽丽 基金项目：国家科技重大专项(2011ZX05016-006)；中央高校基本科研业务费专项资金项目(35832015034) 作者简介：王建俊(1986一)，女，博士研究生，主要从事油气藏工程与渗流机理方面的研究。