低渗透油藏中毛管力在水油两相渗液中的作用

研究在水湿低渗透油藏中，毛管力对水油两相渗流的影响，首先建立一维渗流模型和相应的渗流微分控制方程，用数值差分的方法进行求解，并绘制产量图和采出程度图，分析启动压力梯度和毛管力对水驱油效果的影响。分析表明：在低渗透油藏中，水湿油藏的毛管力和启动压力梯度使油藏开发难度增大，生产成本增加；毛管力使油藏的无水采收率和采收期减小。     

低渗透非均质油藏径向渗流特征研究

利用实验手段回归出启动压力梯度与渗透率的相关式,并用此式得出非均质油藏的启动压力梯度分布,建立了考虑启动压力梯度的三维三相径向坐标下的渗流数学模型和有限差分方程。有限差分方程中每个网格块上的启动压力梯度可以考虑为不同值,网格块上的启动压力梯度值可由初始化参数场给定,编制了相应的软件,并对比了低渗非均质储层不考虑启动压力梯度、启动压力梯度为常数和启动压力梯度非均匀分布时的压力分布特征。     

非达西渗流研究进展

越来越多的资料表明非达西渗流是普遍存在的现象,而且这种现象对某些实际工程的影响不可忽略,必须对其进行总结和研究.针对非达西渗流的两个研究领域--低速非达西渗流和高.速非达西渗流的运动方程等研究成果分别总结;对低速非达西渗流试验方法以及启动压力梯度、拟启动压力梯度和临界压力梯度的产生原因和影响因素进行了论述,且分类介绍了其运动方程;简要介绍了高速非达西渗流的试验方法,并讨论了通过实验曲线建立其运动方程的一般过程,并分类介绍了运动方程;最后概述目前非达西渗流研究的各种问题.     

非线性流下井筒压降影响的水平井产能分析

 大量研究表明,气体渗流受到启动压力梯度、滑脱效应、应力敏感等因素影响,针对气藏的特征,基于稳定渗流理论,在考虑气体高速非达西渗流的基础上,通过气藏渗流与井筒流动耦合,推导建立同时考虑滑脱效应、启动压力梯度、应力敏感等因素共同影响的气藏-水平井井筒耦合模型,并通过实例分析研究了各因素对井筒流量分布的影响.研究结果表明：非线性流特征对井筒流量在井筒跟端分布影响较大;井筒管径越小,绝对粗糙度越大,在靠近水平井的趾端部分井筒流量减小幅度越大,而在靠近水平井的跟端部分井筒流量上升幅度越大;非均质性越强,水平井长度越长,井筒中流量减小幅度越大.     

确定低渗岩芯气体启动压力梯度的一种新方法*

大量实验表明气体在低渗透气藏中渗流存在启动压力梯度。由于气体渗流的启动压力梯度难以直接测定，目前大多通过低压下低速渗流实验，测量气体渗流速度，建立气体渗流速度与压力平方差的关系，根据其非线性特征求解气体渗流的启动压力梯度。但低压下气体滑脱效应的影响在地层条件下几乎不存在，因此，这种方法求得的启动压力梯度值与实际值存在偏差。为了更准确地确定气体启动压力梯度值，提出了一种消除滑脱效应确定气体启动压力梯度的新方法，该方法确定的启动压力梯度为低渗透气田的合理有效开发提供了更为准确的依据。     

安塞油田长10储层渗流特征研究

为了全面研究陕北低渗透油田—安塞油田长10油藏渗流机理,解决开发过程中渗流阻力大、注水效果差等问题,进一步改善油田开发效果,应用室内高精度微量泵和高精度压差传感器对长10储层岩心开展了油水渗流特征、水驱油效率、启动压力梯度等多项实验研究,并推导了低渗透油田IPR新方程,结果表明该储层原油的物性非均质性较强,流体渗流不遵循达西定律,存在启动压力梯度且启动压力梯度与储层的渗透率成反比;油水两相渗流特征反映束缚水饱和度较高,油水两相流动范围窄,具有典型的低渗亲水油藏特征,因此可通过降低启动压力梯度、增大生产压差来改善开发效果。     

叠加原理不能求解含启动压力梯度渗流方程

由于含启动压力梯度的渗流方程不是齐次方程,使用叠加原理对其进行求解时会引起误差。首先建立含启动 压力梯度的渗流方程,然后基于PEBI 网格进行数值求解,最后将数值计算的试井曲线与基于叠加原理的试井曲线进 行对比。数值求解的计算结果表明：启动压力梯度使得关井后的井底压力很快稳定,并且井底压力不能恢复到原始压 力,同时,关井段的压力导数曲线下掉。基于叠加原理的计算结果表明：启动压力梯度使关井后的井底压力持续上升、 数百天乃至数千天都不能使恢复压力稳定,并且压力导数上翘。这不仅与数值解有很大的差别,而且也与封闭边界压 力导数下掉的趋势相悖。这表明叠加原理不适用来求解含启动压力梯度的渗流方程。     

低渗油藏油气两相渗流的理论模型

依据已有的实验结果 ,建立了低渗油藏油气两相渗流的理论模型 .分析表明 ,当启动压力梯度为零时 ,就得到中高渗流油气两相渗流表达式 ,因此统一了油气两相达西流和非达西流的渗流方程 .启动压力梯度的存在是影响低渗油藏开采率低的主要因素 ,当其达到一定值时产量会特别低     

特低渗透油藏启动压力梯度研究

针对单相和油水两相流体启动压力梯度的产生机理和对渗流规律的影响问题,通过稳态“压差—流量”测定启动压力梯度法,分别用模拟油、地层水、注入水和蒸馏水对我国某油田119块典型特低渗岩芯进行驱替实验,得到了单相和油水两相渗流的启动压力梯度值,分析了启动压力梯度的变化规律。物理模拟和理论研究都表明,不同渗流介质测得的启动压力梯度与渗透率均成幂函数关系,并且幂指数近似为-1。同时,由于毛管力和贾敏效应作用,油水两相启动压力梯度远远大于单相渗流的启动压力梯度。结合油田实际,说明了启动压力梯度在特低渗透油气藏开发中的应用。     

低渗透储层流体非线性渗流机理及特征分析

基于微尺度渗流效应，分析了低渗透油藏非线性渗流机理，根据力学平衡方程和压汞实验数据，建立了低渗透储层流体非线性渗流数学模型，分析了非线性渗流特征，并引入动态阻力梯度概念，实现对低渗储层非线性渗流的动态表征。结果表明，边界层、流体屈服应力、边界流体体相流体间的表面力对低渗透渗流形成的附加阻力，是非线性渗流的主要原因；非线性渗流特征在不同的驱替压力梯度下一直存在，在低驱替压力梯度下明显，高驱替压力梯度下弱化；采用动态阻力梯度代替启动压力梯度对非线性渗流进行表征更加合理；随着驱替压力梯度的增加，动态驱替梯度曲线特征为先从一个比较高初始值瞬间降落，之后逐渐增大；动态阻力梯度与驱替压力梯度的比值先增大后减小。     

低渗透油藏聚合物驱启动压力梯度研究

为了描述低渗透油藏聚合物驱的渗流规律及启动压力梯度特征,借鉴水驱启动压力梯度的试验测定方法,将压差-流量法和毛细管平衡法相结合,在不同储层渗透率和体系黏度下开展低渗透储层聚合物驱渗流试验;在此基础上,基于非线性渗流理论量化表征聚合物驱启动压力梯度,建立聚合物驱非线性渗流系数与体系黏度、储层渗透率和流度的量化关系。结果表明,随着储层渗透率的降低或体系黏度的增大,相同渗流速度下聚合物驱的渗流阻力增加、压力梯度增大;聚合物渗流曲线用非线性方法表征吻合较好,且随着储层渗透率降低、体系黏度增大或流度降低,聚合物驱非线性渗流系数均增大;流态图版包含对不流动区、非线性渗流区和拟线性渗流区的定量描述。     

Experimental study on moving boundaries of fluid flow in porous media

Researches on the boundary shape of fluid flow in porous media play an important role in engineering practices, such as petroleum exploitation, nuclear waste disposal and groundwater contamination. In this paper, six types of artificial porous samples （emery jade） with different porosities are manufactured. With the background of slow flow in porous media, laboratory experiments are carried out by observing the movement of five types of fluids with different dynamic viscosities in various types of porous media. A digital video recorder is employed to record the complete process of the fluid flow in the porous media. Based on the digital photos of the moving boundaries of fluid flow in porous media, the average displacement and fractal dimension of the moving boundary are estimated for different combinations of porosity and dynamic viscosity. Moreover, the evolution behavior of the average velocity and fractal dimension of the moving boundary with time is known. The statistical relations of the average velocity, the fractal dimension of the moving boundary and the porosity of porous media and the dynamic viscosity of fluids are proposed in this paper. It is shown that the front shape of the moving boundary of fluid flow in porous media is an integrated result of the porosity of porous media and the dynamic viscosity of fluids.     

水驱粘塑性非牛顿原油渗流特性的数值模拟

针对水驱粘塑性非牛顿原油在地层下呈强非线性渗流特征,建立了描述其流动过程的数学方程; 采用全隐式离散、牛顿－ 拉弗松方法迭代线化,稳定双共轭梯度方法求解,对特定初边条件下二维油藏油水两相流动问题作了数值模拟．结果表明,水驱非牛顿流体驱油效率明显低于水驱牛顿流体     

两种新的不稳定渗流点源函数

对两类有直线型补给边界的一维不稳定渗流问题进行了研究 ,给出了直线型补给边界的定义和表达式。应用积分变换方法解析求解两类一维有源和补给边界的不稳定渗流问题 ,得到两种新的不稳定渗流点源函数 ,该函数可直接用于构造试井分析的数学模型。以改进的矩形油藏试井分析数学模型为例 ,说明了新点源函数的应用方法和表现特征。分析和计算结果表明 ,渗透系数θ反映了补给边界相对渗透能力 ,而传统的定压、封闭边界是直线型补给边界的极限形式 ,渗透系数取给定值时的压力诊断曲线位于定压和封闭边界两种情形之间。     

大孔道内流体流动规律的物理模拟实验

 油藏大孔道内部形态结构既不同于多孔介质也不同于常见管道, 因此, 大孔道内流体流动规律也应不同于常规渗流及粗糙管流, 为方便表述, 本文定义流体通过大孔道的流动为“洞流”。为深入研究大孔道洞流规律, 以因次分析法为基础, 建立大孔道洞流规律物理模型的相似准数, 通过采用不同比例目数的砂粒, 制备了9 个大孔道洞流规律物理模型, 研究不同壁面突出度条件下大孔道的压降及流量关系。研究结果表明:随着流量的增加, 9 个不同突出度大孔道洞流物理模型的流体流动阻力压力梯度均随流量的增加而逐渐增大, 且突出度越大, 压力梯度的增幅越大;当流速增加到一定程度时, 大孔道洞流规律由线性流转变为非线性流;在相同横截面积及相同流动截面积条件下, 洞流的压降梯度小于渗流, 但大于粗糙管流。     

牛顿流体在充满多孔介质圆管内的结构流特征

引入研究多孔介质流体流动的体积平均方法和体均原理,采用Darcy阻力经验公式和Forch-heimer修正,得到了多孔介质内流体流动的Brinkman-Forchheimer广义Darcy定律.提出了牛顿流体在充满多孔介质管道内的结构流概念,并通过边界层动量积分计算证实了流核的存在,获得了流核流速和流核半径的计算结果;分析了Darcy数和惯性参数对速度剖面的影响程度.     

电场作用下多孔介质中单相流体的渗流速度

研究了外电场对多孔介质中单相流体渗流速度的影响。根据渗流力学和电渗理论,推导了外电场作用下多孔介质中单相流体渗流方程,建立了有限元数值模型,采用有限元方法进行了数值模拟。计算结果与实验结果拟合很好,二者之间的相对误差在5%左右。在此基础上探讨了外电场对多孔介质中单相流体渗流速度的影响,结果表明:在压力保持不变情况下电场可以将流体渗流速度增大1~7.5倍;相同电位梯度情况下,压力梯度越大,渗流速度之比越小。     

渗流方程的渗透率自适应权重网格粗化算法

为了提高渗流方程的计算速度和精度,将渗透率自适应网格技术应用于三维非均匀非稳态渗流方程的网格粗化算法中。对于渗透率或孔隙度变化异常的区域,采用精细网格直接求解其压力分布;而在其他区域,采用不均匀网格粗化的方法计算其压力分布。用自适应权重网格粗化算法计算了三维非均匀非稳态渗流场的压力分布。结果表明,三维非均匀非稳态渗流方程的三维不均匀自适应网格粗化算法的解在渗透率或孔隙度异常区域的压力分布规律非常逼近精细网格算法的解,在其他区域的压力分布规律非常逼近粗化算法的解。与采用精细网格算法相比,其计算速度大大提高。     

论多孔介质中流体流动问题的数值模拟方法

简要介绍了多孔介质中流动问题的数值模拟方法及其发展现状 ,对各种方法的优劣及适用范围进行了评价。对多孔介质纯流体耦合流动区域的数值模拟方法及发展现状进行了综合对比。多孔介质模型方法可以用于多孔介质中流体流动问题的数值模拟 ,并且非常适于进行多孔介质纯流体耦合流动区域的数值模拟。     

两相多孔介质拟静态问题的一种有限元解法

对于流体饱和两相多孔介质的拟静态问题。在忽略流体粘性的情况下,从场方程中消去流体相速度变量,得到以固体相位移和孔隙压力为基本变量的控制场方程,并给出相应力值和切值问题的提法。进而采用加权残值有了元法导出基于ｕ＾ｓ－ｐ变量的混合有限元公式,最后了系统方程的解法。