低渗透多孔介质渗流动边界模型的解析与数值解

考虑启动压力梯度的低渗透多孔介质非达西渗流模型属于强非线性动边界问题, 分别利用相似变量变换方法和基于空间坐标变换的有限差分方法, 对内边界变压力情况下、考虑启动压力梯度的一维低渗透多孔介质非达西渗流动边界模型进行了精确解析与数值求解研究. 研究结果表明:该动边界模型存在唯一的精确解析解, 且所求得的精确解析解可严格验证数值解的正确性;且当启动压力梯度值趋于零时, 非达西渗流动边界模型的精确解析解将退化为达西渗流情况下的精确解析解. 由求解结果作出的非零无因次启动压力梯度下的地层压力分布曲线表现出紧支性特点, 其与达西渗流模型的有显著不同. 因此, 研究低渗透多孔介质中非稳态渗流问题时, 应该考虑动边界的影响. 研究内容完善了低渗透多孔介质的非达西渗流力学理论, 为低渗透油气藏开发的试井解释与油藏数值模拟技术提供了理论基础.      

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

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

变形三重介质低渗透油藏垂直裂缝井的非线性流动分析

不仅考虑低渗透油藏具有启动压力梯度的渗流特征,还考虑应力敏感地层中介质的变形;发展了Cinco-Ley H.提出的有限导流垂直裂缝井双线性流理论,将流体在垂直裂缝与地层中形成的流动划分为两个区域—垂直裂缝中的线性流区域和变形三重介质低渗透油藏中的非线性流区域;由此建立了变形三重介质低渗透油藏有限导流垂直裂缝井的非线性流...     

分形煤层气藏数值模拟研究——Douglas-Jones预估-校正法的应用

传统的煤层气动力学模型多建立在欧几里得基础上,难以描述煤层气孔隙结构的复杂性和形状的不规则性。为此,以分形理论为基础,通过引入煤层基质和裂缝的分形维数来刻画煤层气孔隙结构的复杂性和吸附特性,建立了双重分形介质渗流模型,采用Douglas-Jones预估-校正法对非线性方程组进行离散,获得了无限大地层和有限地层定产量生产时拟稳态吸附模型的差分方程,求得数值解。结果表明,Douglas-Jones预估-校正法可以有效解决这类非线性模型的求解问题,获得无限大地层定产量生产时变形双重分形介质模型的数值解;分析各种分形参数下的煤层压力动态,做出了典型压力曲线图。对无限大地层,初期分形维数对压力影响不大,后期分形维数越小,压力越高。对有限地层,初期分形维数的影响明显,且分形维数越大,压力越低。压力随分形指数的减小量呈现先增大后减小的趋势,在末期压力平稳趋向同一值。     

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

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

??????????????????????о?

根据岩心实验提出了两参数的非线性渗流方程,并在此基础上结合Buckly-Leverett 方程建立了非线性非混相驱替方程. 通过对该方程求解分析可知,与线性达西渗流相比,非 线性渗流的驱油相率较差;同一含水饱和度下的含水率更高;而储层中各点的压力梯度更大, 导致所需注采压差更大. 说明有效的开发低渗透油藏需要比中高渗油藏更小的井距和更高的 注采压力体系.     

具有井筒储集的变形介质双渗模型的压力分析

为更好地研究碳酸盐油藏和低渗油的渗流问题,引入渗透率模数,考虑应力敏感地层中介质的变形,介质的双孔隙度、双渗秀率特征,同时考虑井筒储集的影响,建立新的数学模型。渗透率依赖于孔隙压力变化的流动方程是强非线性的,模型采用Douglas-Jones预估-校正法获得了无限大地层及有界封闭地层的数值解,形成了新的理论图版,并利用这些图版对模型中的有关参数进行了敏感性分析。     

具有井简储集的变形介质双渗模型的压力分析

为更好地研究碳酸盐油藏和低渗油的渗流问题，引入渗透率模数，考虑应力敏感地层中介质的变形，介质的双孔隙度、双渗秀率特征，同时考虑井筒储集的影响，建立新的数学模型。渗透率依赖于孔隙压力变化的流动方程是强非线性的，模型采用Douglas—Jones预估-校正法获得了无限大地层及有界封闭地层的数值解，形成了新的理论图版，并利用这些图版对模型中的有关参数进行了敏感性分析。     

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

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

饱和多孔介质中的混合有限元法和有限应变下应变局部化分析

对基于Biot理论的饱和多孔介质中动力-渗流耦合分析提出了一个耦合场混合元．固相位移、应变和有效应力以及流相压力、压力梯度和Darcy速度在单元内均处理为独立变量分别插值．基于胡海昌-Washizu三变量广义变分原理给出的饱和多孔介质动力-渗流耦合问题控制方程的单元弱形式，导出了单元公式．进一步导出了考虑压力相关非关联塑性的非线性单元公式和发展了相应的一致性算法．对几何非线性分析，采用了共旋公式途径．数值结果例题显示所发展耦合场混合元模拟大应变下由应变软化引起以应变局部化为特征的渐进破坏现象的性能．     

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.     

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.     

Mathematical model of two-phase fluid nonlinear flow in low-permeability porous media with applications

IntroductionItisasuccessfulexampleinadevelopmentstoryofscienceandtechnologyformechanicsoffluidsinporousmediatocombinewithengineeringtechnology .Fieldsinfluencedbythemechanicsinvolveddevelopmentofoil_gasandgroundwaterresources,controlonseawaterintrusionandsubsidenceandgeologichazards,geotechnicalengineeringandbioengineering ,andairlineindustry[1～ 7].Aproblemonnonlinearflowinlow_permeabilityporousmediaisbutonlyabasiconeindifferentkindsofengineeringfields,butalsooneoffrontlineresearchfieldsofmod…     

A Fractal Model for Oil Transport in Tight Porous Media

Tight porous media are mainly composed of micro/nano-pores and throats, which leads to obvious microscale effect and nonlinear seepage characteristics. Based on the capillary bundle model and the fractal theory, a new nonlinear seepage equation was deduced, and a further fractal permeability model was obtained for oil transport in tight porous media by considering the effect of the boundary layer. The predictions of the model were then compared with experimental data to demonstrate that the model is valid. This model clarifies the oil transport mechanisms in tight porous media: the effective permeability is no longer a constant value and is governed by properties of tight porous media and oil. Furthermore, parameters influencing effective permeability were analyzed. The model can accurately present the seepage characteristics of the oil in tight porous media and provide a reliable basis for the development of unconventional reservoirs.     

电磁场?渗流场耦合作用下离子液体多孔介质流动模型

离子液体是一类可调控、多功能的绿色环保材料, 具有良好的电磁场响应, 有望应用于调控水驱油路径. 在分析离子液体在毛细管中电磁场响应机理的基础上, 建立了电磁场?渗流场耦合作用下离子液体多孔介质流动模型. 通过理论推导与数值分析发现: 电磁场?渗流场耦合作用下毛细管流量大小主要由离子液体电导率与黏度的比值(内因)、电磁场强度与压力梯度(外因)两方面决定; 电磁场产生的洛伦兹力对离子液体施加一个电磁驱动压强, 形成一个类似压力梯度的电磁驱动等效压力梯度, 从而改变离子液体的流量, 当电磁场强度为2.0 × 104 V/m·T时, 电磁场在电导率为0.5 S/m的离子液体上可形成10 kPa/m电磁驱动等效压力梯度. 通过调整电磁场方向即可控制离子液体在多孔介质中的流动方向, 解决常规注水利用压力差难以控制流动路径的难题, 为离子液体智能驱油提供理论依据, 且电磁场产生的热效应会影响离子液体的流动能力及潜在驱油效率.      

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

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

Numerical Modelling of Free-Surface Flow-Porous Media Interaction Using Divergence-Free Moving Particle Semi-Implicit Method

A divergence-free moving particle semi-implicit method is introduced for free-surface flow through porous media. Numerical incompressibility is conserved by solving additional pressure Poisson equation (PPE). Depending on current particle coordinates, a porosity-based factor is introduced to incorporate the effect of solid volume inside the porous domain. A hybrid formulation containing specified boundary condition and PPE is utilized on free-surface particles. The current framework is tested for four different problems. The first problem shows the effect of the proposed factor in vertical flow through a rectangular porous block and its representative volume change for different phases. Second and third problems validate the numerical model for dam break through a rectangular block of homogeneous porous media. In the fourth problem, flow through a trapezoidal porous block consisting of different porous media with variable effective porosity and permeability is simulated.