The Transient Flow Analysis of Fluid in a Fractal, Double-Porosity Reservoir

The mathematical model for transient fluid flow in porous media is based in general on mass conservation principle. Because of the small compressibility of formation fluid, the quadratic term of pressure gradient is always ignored to linearize the non-linear diffusion equation. This may result in significant errors in model prediction, especially at large time scale. In order to solve this problem, it may be necessary to keep the quadratic term in the non-linear equations. In our study, the quadratic term is reserved to fully describe the transient fluid flow. Based on this rigorous treatment, the mathematical models are established to analyze the transient flow behavior in a double porosity, fractal reservoir with spherical and cylindrical matrix. In addition, Laplace transformation method is employed to solve these mathematical models and the type curves are provided to analyze the pressure transient characteristics. This study indicates that the relative errors in calculated pressure caused by ignoring the quadratic term may amount to 10?% in a fractal reservoir with double porosity, which can??t be neglected in general for fractal reservoirs with double porosity at large time scale.     

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

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

分形油藏不稳定渗流问题的精确解

研究了分形油藏无限大地层和有界地层渗流模型，引入了一类有限广义Ｈａｎｋｅｌ变换，利用这种变换和Ｗｅｂｅｒ变换，在井底定流量和定压生产时，对无限大地层及有界地层（包括封闭和定压地层）六种情况，求得了实空间解析解用双参数（ｄｆ，ｄｓ）来刻画分形油藏的分形特性，分析了分形油藏压力动态特征以及分形参数和边界对压力动态的影响     

Generalized flow analysis of non-Newtonian visco-elastic fluid flow through fractal reservoir

Ｆｒａｃｔａｌｇｅｏｍｅｔｒｙｉｓａｐｏｗｅｒｆｕｌｔｏｏｌｔｏｄｅｓｃｒｉｂｅｃｏｍｐｌｅｘｐｈｅｎｏｍｅｎｏｎ．Ｅｓｐｅｃｉａｌｌｙｉｔｉｓａｐｐｒｏｐｒｉａｔｅｔｏｓｃａｌｅｔｈｅｎｏｎｕｎｉｆｏｒｍｉｔｙａｎｄｎｏｎｓｅｑｕｅｎｃｅｏｆｐｏｒｏｕｓｍｅｄｉａ．Ｉｆｔｈｅｍｅｃｈａｎｉｃｓｏｆｆｌｕｉｄｆｌｏｗｔｈｒｏｕｇｈｐｏｒｏｕｓｍｅｄｉａｉｓｓｔｕｄｉｅｄｂｙｕｓｉｎｇｆｒａｃｔａｌ,ｔｈｅｄｉｓｃｅｒｎｉｂｌｅａｎｄｃｏｇｎｉｔｉｖｅａｂｉｌｉｔｙｆｏｒｐｏｒｏｕｓｍｅｄｉａａｎ…     

幂律流体在多孔介质中径向渗流的分形模型

基于分形理论及毛细管模型,本文研究了非牛顿幂律流体在各向同性多孔介质中径向流动问题,推导了幂律流体径向有效渗透率的分形解析表达式．研究结果表明,幂律流体径向有效无量纲渗透率模型和Chang and Yortsos’s模型吻合很好;同时还得出幂律流体径向有效渗透率随孔隙率、幂指数的增加而增加,随迂曲度分形维数的增加而减少．     

分形理论在煤层气非平衡吸附模型中的应用

传统的煤层气动力学模型均是建立在欧几里得几何基础上的,难以描述煤层孔隙结构的复杂性及形状的不规则性。本文以分形理论为基础,通过引入分形维数来刻画煤层孔隙结构的复杂性并考虑煤层的吸附特性、双重介质特征及介质的变形,建立基于Fick第二定律的分形介质煤层气非稳态渗流数学模型。由于流动方程的强非线性,结合各类边界条件用正则摄动法和Laplace变换得到模型在拉氏空间上的近似解析解,再利用Laplace数值反演求得实空间上的数值解。对参数进行敏感性分析并绘制了典型压力曲线,这些结果为煤层气开采提供了理论依据和试井方法。     

Flow Modeling of Well Test Analysis for Porous–Vuggy Carbonate Reservoirs

Porous–vuggy carbonate reservoirs consist of both matrix and vug systems. This paper represents the first study of flow issues within a porous–vuggy carbonate reservoir that does not introduce a fracture system. The physical properties of matrix and vug systems are quite different in that vugs are dispersed throughout a reservoir. Assuming spherical vugs, symmetrically distributed pressure, centrifugal flow of fluids and considering media that is directly connected with wellbore as the matrix system, we established and solved a model of well testing and rate decline analysis for porous–vuggy carbonate reservoirs, which consists of a dual porosity flow behavior. Standard log–log type curves are drawn up by numerical simulation and the characteristics of type curves are analyzed thoroughly. Numerical simulations showed that concave type curves are dominated by interporosity flow factor, external boundary conditions, and are the typical response of porous–vuggy carbonate reservoirs. Field data interpretation from Tahe oilfield of China were successfully made and some useful reservoir parameters (e.g., permeability and interporosity flow factor) are obtained from well test interpretation.     

变形双重介质广义流动分析

对于碳酸盐油藏和低渗油藏的渗流问题，传统的研究方法都是假设地层渗透率是常数，这假设，对于地层渗透率是压力敏感的情况，对压力的空间变化和瞬时变化将导致较大的误差。本文研究了应力敏感地层中双重介质渗流问题的压力不稳定响应，不仅考虑了储层的双重介质特征，而且考虑了应力敏感地层中介质的变形，建立了应力敏感地层双重介质的数学模型，渗透率依赖于孔隙压力变化的流动方程是强非线性的，采用Douglas-Jones预估－校正法获得了只有裂缝发生形变定产量生产时无限大地层的数值解及定产量生产岩块与裂隙同时发生形变时无限大地层的数值解，并探讨了变形参数和双重介质参数变化时压力的变化规律，给出几种情况下典型压力曲线图版，这些结果可用于实际试井分析。     

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

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

Stable Propagation of Saturation Overshoots for Two-Phase Flow in Porous Media

Propagation of saturation overshoots for two-phase flow of immiscible and incompressible fluids in porous media is analyzed using different computational methods. In particular, it is investigated under which conditions a given saturation overshoot remains stable while moving through a porous medium. Two standard formulations are employed in this investigation, a fractional flow formulation and a pressure–saturation formulation. Neumann boundary conditions for pressure are shown to emulate flux boundary conditions in homogeneous media. Gravity driven flows with Dirichlet boundary conditions for pressure that model infiltration into heterogeneous media with position-dependent permeability are found to exhibit pronounced saturation overshoots very similar to those seen in experiment.     

Reservoir fluid flow through vertical fracture of thick stratum with underlying fluid interface

An exact solution is obtained for the problem of steady-state filtration of a heavy dense incompressible fluid in a thin, infinitely deep, inclined reservoir having a crack of given depth along the reservoir rise. The region of filtration of the lighter liquid (oil) has an impermeable upper boundary in the form of a horizontal fault line. Below the filtration region there is a free boundary, below which lies the region of stationary fluid (bottom water). The interface of the fluids, the fissure profile, and the reservoir fluid flow rate are determined from the solution of the problem on the basis of the given parameters (permeability of the reservoir and of the material filling the fissure, viscosity of the filtering fluid, specific weight of the upper and lower fluids, depth of the fissure, pressure differential between a point at the fissure and a point at the interface of the fluids). In the case when the thin reservoir is a vertical filtering layer, the considered flow is interpreted as the motion of the reservoir fluid through a vertical fissure of a thick reservoir (half-space) in the presence of an underlying fluid interface. The problem is solved in finite form with the aid of known analytic functions using integrals of the Cauchy type. The fundamental solution is first found of the special problem of flow with a point singularity. The fundamental solution is also of independent importance as an extension of the solution of certain known problems [1–4].     

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.     

Microscopic mechanism of periodical electroosmosis in reservoir rocks

Based on the electric double layer (EDL) theory and the momentum equation governing the electroosmosis flow, this paper presents an analytical solution to the periodical electroosmosis with a parallel straight capillary bundle model of reservoir rocks to reveal the microscopic mechanism of the electroosmotic flows in rocks. The theory shows that both the frequency dispersion characteristics of the macroscopic electroosmotic Darcy velocity in unsealed rocks and the electroosmotic pressure coefficient in sealed rocks depend on the porosity and electrochemical properties of reservoir rocks. The mathematical simulation indicates that the distribution of the periodical electroosmotic velocity is wavelike in the rock pore. The greater the porosity is, the greater electroosmotic the Darcy velocity and the smaller electroosmotic pressure coefficient are generated. The module values of the electroosmotic Darcy velocity and the electroosmotic pressure coefficient increase with the decreasing solution concentration or the increasing cation exchange capacity without affecting the phase of the electroosmotic Darcy velocity.     

考虑诱导缝多段压裂水平井非均质改造压力动态模型研究

为了准确模拟致密油藏水平井大规模压裂形成复杂裂缝网络系统和非均质储层井底压力变化,建立考虑诱导缝矩形非均质储层多段压裂水平井不稳定渗流数学模型,耦合裂缝模型与储层模型得到有限导流裂缝拉普拉斯空间井底压力解,对两种非均质储层模型分别利用数值解、边界元和已有模型验证其准确性.基于压力导数曲线特征进行流动阶段划分和参数敏感性分析,得到以下结果:和常规压裂水平井井底压力导数曲线相比较,理想模式下,考虑诱导缝影响时特有的流动阶段是综合线性流阶段、诱导缝向压裂裂缝“补充”阶段、储层线性流动阶段和拟边界控制流阶段.诱导缝条数的增加加剧了综合线性流阶段的持续时间,降低了流体渗流阻力,早期阶段压力曲线越低;当诱导缝与压裂裂缝导流能力一定时,裂缝导流能力越大,线性流持续时间越长;当所有压裂裂缝不在一个区域时,沿井筒方向两端区域低渗透率弱化了低渗区域诱导缝流体向压裂裂缝“补充”阶段,因此,沿井筒方向两端区域渗透率越低,早期阶段压力曲线越高;当所有压裂裂缝在一个区域时,渗透率变化只影响径向流阶段之后压力曲线形态,外区渗透率越低,早期径向流阶段之后压力曲线越高.通过实例验证,表明该模型和方法的实用性和准确性.     

Exact solutions for nonlinear transient flow model including a quadratic gradient term

The models of the nonlinear radial flow for the infinite and finite reservoirs including a quadratic gradient term were presented. The exact solution was given in real space for flow equation including quadratic gradiet term for both constant-rate and constant pressure production cases in an infinite system by using generalized Weber transform.Analytical solutions for flow equation including quadratic gradient term were also obtained by using the Hankel transform for a finite circular reservoir case. Both closed and constant pressure outer boundary conditions are considered. Moreover, both constant rate and constant pressure inner boundary conditions are considered. The difference between the nonlinear pressure solution and linear pressure solution is analyzed. The difference may be reached about 8% in the long time. The effect of the quadratic gradient term in the large time well test is considered.     

基于分形理论的煤层气非达西流动分析

在考虑了煤层的分形特征和启动压力梯度影响的基础上,建立了无限大煤层中气体低速非达西流动的数学模型,并求得了量纲为一的井底压力的Laplace空间解析解,并根据数值求解结果绘制了典型的井底压力动态曲线。研究结果表明:在定产量生产的情况下,分形维数和量纲为一的启动压力梯度对早期井底压力动态无显著影响;在生产的中后期,由于受二者的影响,压力导数曲线上的径向流水平直线段消失;量纲为一的井储系数的影响主要表现在续流阶段,而吸附系数则主要影响煤基质向裂隙扩散的过渡阶段。     

Flow processes with change of regime in fractured reservoirs

Experimental and industrial observations indicate a strong nonlinear dependence of the parameters of the flow processes in a fractured reservoir on its state of stress. Two problems with change of boundary condition at the well — pressure recovery and transition from constant flow to fixed bottom pressure — are analyzed for such a reservoir. The latter problem may be formulated, for example, so as not to permit closure of the fractures in the bottom zone. For comparison, the cases of linear [1] and nonlinear [2] fractured porous media and a fractured medium [3] are considered, and solutions are obtained in a unified manner using the integral method described in [1]. Nonlinear elastic flow regimes were previously considered in [3–6], where the pressure recovery process was investigated in the linearized formulation. Problems involving a change of well operating regime were examined for a porous reservoir in [7].Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 3, pp. 67–73, May–June, 1991.     

Effect of Hydrogen Sulfide on Condensation and Flow of Gas Condensate Fluids in Porous Media

Condensation and flow experiments were conducted at subsurface conditions in a glass micromodel using reservoir fluids with and without the hydrogen sulfide component. It has been noted that the formation of the condensing phase as well as modes of condensate flow are similar for both fluids. Furthermore, an additional condensate transport mechanism, termed lamella flow, was observed with the sour fluid. It has been concluded that core flow experiments conducted with sweet reservoir fluid should reproduce the flow of sour fluid to a large extent.     

Flow of Particulate-Fluid Suspension in a Channel with Porous Walls

The problem of the dispersed particulate-fluid two-phase flow in a channel with permeable walls under the effect of the Beavers and Joseph slip boundary condition is concerned in this paper. The analytical solution has been derived for the longitude pressure difference, stream functions, and the velocity distribution with the perturbation method based on a small width to length ratio of the channel. The graphical results for pressure, velocity, and stream function are presented and the effects of geometrical coefficients, the slip parameter and the volume fraction density on the pressure variation, the streamline structure and the velocity distribution are evaluated numerically and discussed. It is shown that the sinusoidal channel, accompanied by a higher friction factor, has higher pressure drop than that of the parallel-plate channel under fully developed flow conditions due to the wall-induced curvature effect. The increment of the channel’s width to the length ratio will remarkably increase the flow rate because of the enlargement of the flow area in the channel. At low Reynolds number ranging from 0 to 65, the fluids move forward smoothly following the shape of the channel. Moreover, the slip boundary condition will notably increase the fluid velocity and the decrease of the slip parameter leads to the increment of the velocity magnitude across the channel. The fluid-phase axial velocity decreases with the increment of the volume fraction density.     

Nonlinear Pressure Diffusion in Flow of Compressible Liquids Through Porous Media

In the flow of liquids through porous media, nonlinear effects arise from the dependence of the fluid density, porosity, and permeability on pore pressure, which are commonly approximated by simple exponential functions. The resulting flow equation contains a squared gradient term and an exponential dependence of the hydraulic diffusivity on pressure. In the limiting case where the porosity and permeability moduli are comparable, the diffusivity is constant, and the squared gradient term can be removed by introducing a new variable y, depending exponentially on pressure. The published transformations that have been used for this purpose are shown to be special cases of the Cole–Hopf transformation, differing in the choice of integration constants. Application of Laplace transformation to the linear diffusion equation satisfied by y is considered, with particular reference to the effects of the transformation on the boundary conditions. The minimum fluid compressibilities at which nonlinear effects become significant are determined for steady flow between parallel planes and cylinders at constant pressure. Calculations show that the liquid densities obtained from the simple compressibility equation of state agree to within 1% with those obtained from the highly accurate Wagner-Pru?  equation of state at pressures to 20 MPa and temperatures approaching 600 K, suggesting possible applications to some geothermal systems.