基于渗透率连续变化的低渗透多孔介质非线性渗流模型

基于低渗透多孔介质渗透率的渐变理论,确定了能精确描述低渗透多孔介质渗流特征的非线性运动方程,并通过实验数据拟合．验证了非线性运动方程的有效性。非线性渗流速度关于压力梯度具有连续-阶导数,方便于工程计算;由此建立了低渗透多孔介质的单相非线性径向渗流数学模型,并巧妙采用高效的Douglas-Jones预估一校正有限差分方法求得了其数值解。数值结果分析表明：非线性渗流模型为介于拟线性渗流模型和达西渗流模型之间的一种中间模型或理想模型,非线性渗流模型和拟线性渗流模型均存在动边界;拟线性渗流高估了启动压力梯度的影响,使得动边界的移动速度比实际情况慢得多;非线性越强,地层压力下降的范围越小,地层压力梯度越陡峭,影响地层压力的敏感性减弱,而影响地层压力梯度的敏感性增强。     

低渗透煤层气藏中气-水两相不稳定渗流动态分析

针对低渗透煤层渗流问题,考虑启动压力梯度及其引起的动边界和动边界内吸附气解吸作用的渗流模型研究目前仅限于单相流,而更符合实际的气-水两相渗流动边界模型未见报道.本文综合考虑了煤层吸附气的解吸作用、气-水两相渗流、非达西渗流、地层应力敏感等影响因素,进行了低渗透煤层的气-水两相渗流模型研究.采用了试井技术中的"分相处理"方法,修正了两相渗流的综合压缩系数和流度,并基于含气饱和度呈线性递减分布的假设,建立了煤层气藏的气-水两相渗流耦合模型.该数学模型不仅可以描述由于低渗透煤层中渗流存在启动压力梯度而产生的可表征煤层有效动用范围随时间变化的移动边界,还可以描述煤层有效动用范围内吸附气的解吸现象以及吸附气解吸作用所引起的煤层含气饱和度的上升;为了提高模型精度,控制方程还保留了二次压力梯度项.采用了稳定的全隐式有限差分方法进行了模型的数值求解,并验证了数值计算方法的正确性,获得了模型关于瞬时井底压力与压力导数响应的双对数特征曲线,由此分析了各渗流参数的敏感性影响.本文研究结果可为低渗透煤层气藏开发的气-水两相流试井技术提供渗流力学的理论基础.     

考虑二次梯度项及动边界的双重介质低渗透油藏流动分析

在传统试井模型的非线性偏微分方程中根据弱可压缩流体的假设,忽略了二次梯度项,对于低渗透油藏这种方法是有疑问的.低渗透问题一个显著的特点就是流体的流动边界随着时间不断向外扩展.为了更好地研究双重介质低渗透油藏中流体的流动问题,考虑了二次梯度项及活动边界的影响,同时考虑了低渗透油藏的非达西渗流特征,建立了双重介质低渗透油藏流动模型.采用Douglas-Jones预估-校正差分方法获得了无限大地层定产量生产时模型的数值解,分别讨论了不同参数变化时压力的变化规律及活动边界随时间的传播规律,还分析了考虑和忽略二次梯度项影响时模型数值解之间的差异随时间的变化规律,做出了典型压力曲线图版,这些结果可用于实际试井分析.     

动边界双重介质油藏低速非达西渗流试井模型

裂缝性油藏中基质岩块的渗透率一般很低,大量岩心测试实验证实在基质岩块内的液体渗流和在一定含水饱和度下的气体渗流将偏离达西渗流,往往出现低速非达西渗流,表现出启动压力梯度以及流体流动边界不断向外扩展等特殊现象。本文充分考虑启动压力梯度与动边界的影响,建立了微可压缩双重介质油藏低速非达西渗流的试井数学模型,对时间和空间变量...     

低速非达西渗流动边界问题的积的分解

研究了低渗透油藏低速非达西径向流的动边界问题,给出了高精度的积分解,分析了启动压力梯度对压力分布的影响,发现启动压力梯度越大,井底附近压力下降越快,外边界传播越慢．     

低渗油层压裂水平井两相流研究

依据压裂水平井不同流动区域的流动规律, 将压裂水平井的渗流分为裂缝中的高速非达西流 动区、裂缝控制影响的椭圆渗流区、远离裂缝的基质非达西渗流区, 考虑启动压力梯度的影 响, 对压裂水平井两相渗流进行了分析, 得到了低渗透油层压裂水平井的产量公式. 研究结 果表明, 裂缝的导流能力越大, 压裂水平井的产量越高. 但随着开采时间的增加, 其产量递减幅度越 大; 压裂裂缝长度越小, 压裂水平井的初始产量越高. 但随着生产时间的推移, 压裂裂缝的 长度越大, 产量的递减幅度越小; 中间裂缝长, 两翼裂缝短的情况下, 压裂水平井的产量最 高.     

低渗透煤层气有限导流垂直裂缝井压力动态分析

作为典型的致密多孔介质,煤岩储层已被证实存在启动压力梯度.根据煤层气垂直裂缝井的双线性流动机制,综合考虑启动压力梯度和井筒储存效应的影响,建立了一个新的低渗透煤层气有限导流垂直裂缝井双线性流动数学模型,采用Laplace变换和Stehfest数值反演方法对数学模型进行了求解,并分析了无因次启动压力梯度等参数对无因次井底压力及其导数曲线的影响规律.分析结果表明:典型的低渗透煤层气垂直裂缝井双线性流动曲线可划分为早期续流段、双线性流段、煤层线性流段、过渡流段和煤层边界线性流段5个特征阶段,其中由于启动压力梯度存在的影响,无因次井底压力及其导数曲线自煤层线性流段开始出现明显上翘,且启动压力梯度值越大,曲线上翘趋势越明显;此外,煤层边界线性流段呈现为单位斜率的直线,而非1/2斜率的线性流段直线.这些结果表现出启动压力梯度对低渗透煤层气垂直裂缝井双线性流动的影响,可用于指导现场煤层气井试井分析.     

一种微管中启动压力梯度测量方法

为了更加准确地测量低渗透多孔介质中液体流动的启动压力梯度,本文设计了一种在微管中测量启动压力梯度的方法:将静态法和稳态流动实验相结合,研究微管内的去离子水从静态到流动状态整体的压力反应. 实验结果表明:由稳态流动实验得到微管启动压力梯度,是静态法所得启动压力梯度值的14.5倍,说明在以往通过动态实验值所推得启动压力梯度数值过大. 这可以解释为在油藏工程中,实验室测得的启动压力梯度往往过大,无法向现场推广的原因. 研究表明:可以将该方法用于低渗透岩心的启动压力梯度测量实验中,更加准确的获得低渗透油藏的启动压力梯度值. 本文的研究不仅证明了微米尺度下启动压力梯度的存在,而且给出了更加准确测量启动压力梯度的方法.     

多孔介质非线性渗流问题的摄动解

考虑变形多孔介质渗透参数(渗透率和孔隙度)与孔隙压力呈负指数变化的特点，建立了多孔介质渗流问题的数学模型，采用积分变换方法求出了一维非线性渗流问题的摄动解，并对常数渗透参数和指数渗透参数的渗流问题进行对比分析，计算结果表明：两者之间的差别较大，且渗透参数的变化对于流体渗流中后期过程有着重要的影响，但对渗流早期影响不大，这对于定量研究工程中非线性渗流问题模型参数的相对重要性提供了可靠的理论依据。     

缝洞型多孔介质中多相流的有限体积法数值模拟

本文研究的碳酸盐岩油藏储集体属于缝洞型多孔介质.这类缝洞型多孔介质由裂缝、溶蚀孔洞和低孔隙度低渗透率的基岩组成.裂缝是空隙流体流动的主要通道;溶蚀孔洞大小从几厘米到数米不等,渗透率和孔隙度都很高,是流体主要的储集空间.由于缝洞型多孔介质空隙空间的复杂性和强非均质性,数值计算中基本控制方程的空间离散应采用非结构化网格的计算模型.本文采用有限体积法模拟缝洞型多孔介质中多相流体的流动,并给出了相应的单元中心格式有限体积法的计算公式.裂缝介质和溶洞介质中单元间多相流体的流动考虑为高速非达西流,其质量通量采用Forchheimer定律计算.非线性方程的离散选取全隐式格式,并采用Newton-Raphson迭代进行求解.通过两个二维模型注水驱油的数值模拟,验证了本文方法的有效性.     

Effect of quadratic pressure gradient term on a one-dimensional moving boundary problem based on modified Darcy’s law

A relatively high formation pressure gradient can exist in seepage flow in low-permeable porous media with a threshold pressure gradient, and a significant error may then be caused in the model computation by neglecting the quadratic pressure gradient term in the governing equations. Based on these concerns, in consideration of the quadratic pressure gradient term, a basic moving boundary model is constructed for a one-dimensional seepage flow problem with a threshold pressure gradient. Owing to a strong nonlinearity and the existing moving boundary in the mathematical model, a corresponding numerical solution method is presented. First, a spatial coordinate transformation method is adopted in order to transform the system of partial differential equations with moving boundary conditions into a closed system with fixed boundary conditions; then the solution can be stably numerically obtained by a fully implicit finite-difference method. The validity of the numerical method is verified by a published exact analytical solution. Furthermore, to compare with Darcy’s flow problem, the exact analytical solution for the case of Darcy’s flow considering the quadratic pressure gradient term is also derived by an inverse Laplace transform. A comparison of these model solutions leads to the conclusion that such moving boundary problems must incorporate the quadratic pressure gradient term in their governing equations; the sensitive effects of the quadratic pressure gradient term tend to diminish, with the dimensionless threshold pressure gradient increasing for the one-dimensional problem.     

Analytical solution of a double moving boundary problem for nonlinear flows in one-dimensional semi-infinite long porous media with low permeability简

Based on Huang＇s accurate tri-sectional nonlin- ear kinematic equation （1997）, a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction： The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation （PDE） sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient （TPG） has a great effect on the external moving boundary but has little effect on the internal moving boundary.     

Analysis of Multiphase Non-Darcy Flow in Porous Media

Recent laboratory studies and analyses (Lai et al. Presented at the 2009 Rocky Mountain Petroleum Technology Conference, 14–16 April, Denver, CO, 2009) have shown that the Barree and Conway model is able to describe the entire range of relationships between flow rate and potential gradient from low- to high-flow rates through porous media. A Buckley and Leverett type analytical solution is derived for non-Darcy displacement of immiscible fluids in porous media, in which non-Darcy flow is described using the Barree and Conway model. The comparison between Forchheimer and Barree and Conway non-Darcy models is discussed. We also present a general mathematical and numerical model for incorporating the Barree and Conway model in a general reservoir simulator to simulate multiphase non-Darcy flow in porous media. As an application example, we use the analytical solution to verify the numerical solution for and to obtain some insight into one-dimensional non-Darcy displacement of two immiscible fluids with the Barree and Conway model. The results show how non-Darcy displacement is controlled not only by relative permeability, but also by non-Darcy coefficients, characteristic length, and injection rates. Overall, this study provides an analysis approach for modeling multiphase non-Darcy flow in reservoirs according to the Barree and Conway model.     

A New Non-Darcy Flow Model for Low-Velocity Multiphase Flow in Tight Reservoirs

The pore and pore-throat sizes of shale and tight rock formations are on the order of tens of nanometers. The fluid flow in such small pores is significantly affected by walls of pores and pore-throats. This boundary layer effect on fluid flow in tight rocks has been investigated through laboratory work on capillary tubes. It is observed that low permeability is associated with large boundary layer effect on fluid flow. The experimental results from a single capillary tube are extended to a bundle of tubes and finally to porous media of tight formations. A physics-based, non-Darcy low-velocity flow equation is derived to account for the boundary layer effect of tight reservoirs by adding a non-Darcy coefficient term. This non-Darcy equation describes the fluid flow more accurately for tight oil reservoir with low production rate and low pressure gradient. Both analytical and numerical solutions are obtained for the new non-Darcy flow model. First, a Buckley–Leverett-type analytical solution is derived with this non-Darcy flow equation. Then, a numerical model has been developed for implementing this non-Darcy flow model for accurate simulation of multidimensional porous and fractured tight oil reservoirs. Finally, the numerical studies on an actual field example in China demonstrate the non-negligible effect of boundary layer on fluid flow in tight formations.     

格子Boltzmann方法模拟多孔介质惯性流的边界条件改进

格子Boltzmann方法可以有效地模拟水动力学问题,边界处理方法的选择对于可靠的模拟计算至关重要.本文基于多松弛时间格子Boltzmann模型开展了不同边界条件下,周期对称性结构和不规则结构中流体流动模拟,阐述了不同边界条件的精度和适用范围. 此外,引入一种混合式边界处理方法来模拟多孔介质惯性流, 结果表明:对于周期性对称结构流动模拟,体力格式边界条件和压力边界处理方法是等效的,两者都能精确地捕捉流体流动特点; 而对于非周期性不规则结构,两种边界处理方法并不等价,体力格式边界条件只适用于周期性结构;由于广义化周期性边界条件忽略了垂直主流方向上流体与固体格点的碰撞作用,同样不适合处理不规则模型;体力-压力混合式边界格式能够用来模拟周期性或非周期性结构流体流动,在模拟多孔介质流体惯性流时,比压力边界条件有更大的应用优势,可以获得更大的雷诺数且能保证计算的准确性.      

A VARIATIONAL INEQUALITY APPROACH FOR NON-STEADY SEEPAGE FLOW THROUGH THREE-DIMENSIONAL FRACTURE NETWORK

为了求解裂隙岩体有自由面非稳定渗流问题,将Darcy定律延拓至整个研究区域,使得潜在溢出边界条件满足Signorini型边界条件,建立了三维裂隙网络非稳定渗流问题的抛物型变分不等式（parabolic variational inequality,PVI）提法,并证明其与偏微分方程（partial differential equation,PDE）提法的等价性,从而将自由面上的流量条件以及潜在溢出边界上的互补条件转化成自然边界条件,降低该问题求解难度。同时给出了基于PVI提法的有限元数值求解方法,通过与交叉裂隙模型理论解的对比分析,证明了该方法的正确性。最后将该方法对含复杂三维裂隙网络的边坡进行非稳定渗流分析,计算结果表明该方法对于复杂裂隙网络求解具有较强的可靠性和适应性。     

页岩气藏压裂水平井试井分析

页岩气藏资源丰富,开发潜力巨大,已成为目前研究的热点.与常规气藏相比,页岩气藏运移机制复杂,流动模式呈非线性,有必要考虑页岩气的吸附解吸,天然微裂缝的应力敏感性,人工裂缝内的非达西流等非线性因素对压裂水平井压力响应的影响. 基于双重介质和离散裂缝混合模型,分别采用Langmuir等温吸附方程描述吸附解吸,渗透率指数模型描述应力敏感,Forchheimer方程描述非达西效应,建立页岩气藏压裂水平井数值试井模型. 运用伽辽金有限元法对模型进行求解.根据试井特征曲线,划分流动阶段,着重分析非线性因素对压力响应的影响.结果表明:页岩气藏压裂水平井存在压裂裂缝线性流、压裂裂缝径向流、地层线性流、系统径向流及封闭边界影响5 种流动阶段.吸附解吸的影响发生窜流之后,Langmuir吸附体积增大,拟压力导数曲线凹槽更加明显,系统径向流出现时间与压力波传播到边界时间均延迟;天然裂缝系统的应力敏感性主要影响试井曲线的晚期段,拟压力和拟压力导数曲线均表现为上翘,应力敏感效应越强,上翘幅度越大;高速非达西效应对早期段影响较大,非达西效应越强,拟压力降幅度越大,试井曲线上翘.与解析解的对比以及矿场实例验证了模型的正确性与适用性.