Approximate Solutions for Non-stationary Gas Permeability Tests

Full and approximate models describing changes in gas pressure for non-stationary single- and double-reservoir methods of testing of gas permeability of porous materials were formulated. The models include effects of compressibility and slip of gas. For derivation of the approximate models it was assumed that a transient mass flux of gas is constant along the tested sample (spatially homogeneous). The assumption allowed finding analytical solutions for evolution of reservoir pressures as well as the transient spatial distributions of pore pressure in tested samples. The full model was solved with help of the finite element method and the results were compared with the derived approximations. The comparison shows very good agreement of predictions of both models. Parametric studies have given the evaluation of the role of selected material parameters (permeability, porosity and slip factor) and parameters of the experimental system (initial pressure, volume of reservoirs), especially important from the point of view of identification procedures of the material parameters.     

End effects in rarefied gas outflow into vacuum through a long tube

The problem of steady outflow of rarefied gas from a reservoir into vacuum through a long cylindrical tube of circular cross-section at a constant temperature is considered on the basis of the kinetic S-model. The kinetic equation is solved numerically using a conservative second-order method. The basic calculated quantity is the gas flow rate through the channel; the flowfields are also studied. The solutions obtained are compared with the known results.     

Hydraulic model for oil and gas reservoirs

This paper considers a two-dimensional hydraulic model which describes gas or oil flow to a horizontal well with hydraulic fractures and takes into account the reservoir geometry and fluid flow between the reservoir and the well. A computational algorithm is proposed, and calculations for gas and oil reservoirs are performed. A comparison of the calculation results and the solutions of the corresponding problems in a three-dimensional formulation show that the calculations using the approximate hydraulic model yield reasonably accurate results.     

Friedrichs systems for systems of wave equations and shear waves in a three-dimensional elastic medium

This paper considers a two-dimensional hydraulic model which describes gas or oil flow to a horizontal well with hydraulic fractures and takes into account the reservoir geometry and fluid flow between the reservoir and the well. A computational algorithm is proposed, and calculations for gas and oil reservoirs are performed. A comparison of the calculation results and the solutions of the corresponding problems in a three-dimensional formulation show that the calculations using the approximate hydraulic model yield reasonably accurate results.     

Gas Flow in Porous Media With Klinkenberg Effects

Gas flow in porous media differs from liquid flow because of the large gas compressibility and pressure-dependent effective permeability. The latter effect, named after Klinkenberg, may have significant impact on gas flow behavior, especially in low permeability media, but it has been ignored in most of the previous studies because of the mathematical difficulty in handling the additional nonlinear term in the gas flow governing equation. This paper presents a set of new analytical solutions developed for analyzing steady-state and transient gas flow through porous media including Klinkenberg effects. The analytical solutions are obtained using a new form of gas flow governing equation that incorporates the Klinkenberg effect. Additional analytical solutions for one-, two- and three-dimensional gas flow in porous media could be readily derived by the following solution procedures in this paper. Furthermore, the validity of the conventional assumption used for linearizing the gas flow equation has been examined. A generally applicable procedure has been developed for accurate evaluation of the analytical solutions which use a linearized diffusivity for transient gas flow. As application examples, the new analytical solutions have been used to verify numerical solutions, and to design new laboratory and field testing techniques to determine the Klinkenberg parameters. The proposed laboratory analysis method is also used to analyze data from steady-state flow tests of three core plugs from The Geysers geothermal field. We show that this new approach and the traditional method of Klinkenberg yield similar results of Klinkenberg constants for the laboratory tests; however, the new method allows one to analyze data from both transient and steady-state tests in various flow geometries.     

Analytical Solutions to Gas Flow Problems in Low Permeability Porous Media

In most of conventional porous media the flow of gas is basically controlled by the permeability and the contribution of gas flow due to gas diffusion is ignored. The diffusion effect may have significant impact on gas flow behavior, especially in low permeability porous media. In this study, a dual mechanism based on Darcy flow as well as diffusion is presented for the gas flow in homogeneous porous media. Then, a novel form of pseudo pressure function was defined. This study presents a set of novel analytical solutions developed for analyzing steady-state and transient gas flow through porous media including effective diffusion. The analytical solutions are obtained using the real gas pseudo pressure function that incorporates the effective diffusion. Furthermore, the conventional assumption was used for linearizing the gas flow equation. As application examples, the new analytical solutions have been used to design new laboratory and field testing method to determine the porous media parameters. The proposed laboratory analysis method is also used to analyze data from steady-state flow tests of three core plugs. Then, permeability (k) and effective diffusion coefficient (D _e) was determined; however, the new method allows one to analyze data from both transient and steady-state tests in various flow geometries.     

On Vertical Diffusion of Gases in a Horizontal Reservoir

Exact and approximate solutions to vertical diffusion in gravity-stable, ideal gas mixtures in gas reservoirs, depleted oil reservoirs, or drained aquifers are presented, and characteristic times of diffusion are obtained. Our solutions also can be used to test numerical simulators that model diffusion after gas injection. First, we consider isothermal, countercurrent vertical diffusion of carbon dioxide and methane in a horizontally homogeneous reservoir. Initially, the bottom part of the reservoir, with no flow boundaries at the top and bottom, is filled with CO₂ and the upper part with CH₄. At time equal zero, the two gases begin to diffuse. We obtain the exact solution to the initial and boundary-value problem using Fourier series method. For the same problem, we also obtain an approximate solution using the integrated mass balance method. The latter solution has a particularly simple structure, provides a good approximation and retains the important features of the exact solution. Its simplicity allows one to perform calculations that are difficult and non-transparent with the Fourier series method. It also can be used to test numerical algorithms. Furthermore, we consider diffusion of CO₂ with partitioning into connate water. We show that at reservoir pressures the CO₂ retardation by water cannot be neglected. The diffusion-retardation problem is modelled by a non-linear diffusion equation whose self-similar solution is obtained. Finally, we obtain a self-similar solution to a nonlinear diffusion problem. This solution is a good approximation at early times, before the diffusing gases reach considerable concentrations at the top and bottom boundaries of the reservoir.     

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

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

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.     

Pressure transients during injection: Constant rate and convolution solutions

Pressure transient testing during water injection is undertaken to evaluate the injection potential of a well. If properly interpreted, it can yield information such as effective mobilities of fluids, wellbore damage, and residual oil saturation. This is best done by the simultaneous use of downhole flow-rate and pressure measurements. Analytical solutions obtained under various assumptions for pressure response of an injection well are investigated. For a constant downhole flow rate, it is demonstrated that exact solutions may be obtained for an infinite reservoir during both the injection and the falloff periods. Due to the inherent nonlinearity of the problem, the constant rate solutions are not readily extended for the general case of varying flow rates. Therefore, we have employed an approximate technique. This technique can be used with an altered form of convolution and permits calculation of the pressure response for arbitrary rate data. The range of parameters under which this method may be used are also identified. Computational methods related to convolution are presented. The numerical techniques developed in this paper are sufficiently general that they may be applied to similar well-testing problems involving single-phase flow.     

Mass flow rate measurements in gas micro flows

The main objective of this experimental investigation on the gas flow slip regime is to measure the mass flow rate in isothermal steady flows through cylindrical micro tubes. Two technical procedures devoted to mass flow rate measurements are compared, and the measured values are also compared with the results yielded by different approximated analytical solutions of the gas dynamics continuum equations. Satisfactory results are obtained and the way is clearly opened to measuring mass flow rates for higher Knudsen numbers, over all the micro flow transitional regime.     

Three-Phase Flow Modelling Using Pore-Scale Capillary Pressures and Relative Permeabilities for Mixed-Wet Media at the Continuum-Scale

When regions of three-phase flow arise in an oil reservoir, each of the flow parameters, i.e. capillary pressures and relative permeabilities, are generally functions of two phase saturations and depend on the wettability state. The idea of this work is to generate consistent pore-scale based three-phase capillary pressures and relative permeabilities. These are then used as input to a 1-D continuum core- or reservoir-scale simulator. The pore-scale model comprises a bundle of cylindrical capillary tubes, which has a distribution of radii and a prescribed wettability state. Contrary to a full pore-network model, the bundle model allows us to obtain the flow functions for the saturations produced at the continuum-scale iteratively. Hence, the complex dependencies of relative permeability and capillary pressure on saturation are directly taken care of. Simulations of gas injection are performed for different initial water and oil saturations, with and without capillary pressures, to demonstrate how the wettability state, incorporated in the pore-scale based flow functions, affects the continuum-scale displacement patterns and saturation profiles. In general, wettability has a major impact on the displacements, even when capillary pressure is suppressed. Moreover, displacement paths produced at the pore-scale and at the continuum-scale models are similar, but they never completely coincide.     

Similarity solutions for the flow behind an exponential shock in a non-ideal gas

Similarity solutions for the flow of a non-ideal gas behind a strong exponential shock driven out by a piston (cylindrical or spherical) moving with time according to an exponential law are obtained. Similarity solutions exist only when the surrounding medium is of constant density. Solutions are obtained, in both the cases, when the flow between the shock and the piston is isothermal or adiabatic. It is found that the assumption of zero temperature gradient brings a profound change in the density distribution as compare to that of the adiabatic case. Effects of the non-idealness of the gas on the flow-field between the shock and the piston are investigated. The variations of density-ratio across the shock and the location of the piston with the parameter of non-idealness of the gas are also obtained.     

A note on sinusoidal motion of a viscoelastic non-Newtonian fluid

In this note, the exact solutions of velocity field and associated shear stress corresponding to the flow of second-grade fluid in a cylindrical pipe, subject to a sinusoidal shear stress, are determined by means of Laplace and finite Hankel transform. These solutions are written as sum of steady-state and transient solutions, and they satisfy governing equations and all imposed initial and boundary conditions. The corresponding solutions for the Newtonian fluid, performing the same motion, can be obtained from our general solutions. At the end of this note, the effects of different parameters are presented and discussed by showing flow profiles graphically.     

Semi-Analytical Solutions for a Partially Penetrated Well with Wellbore Storage and Skin Effects in a Double-Porosity System with a Gas Cap

We have studied the effect of a constant top pressure on the pressure transient analysis of a partially penetrated well in an infinite-acting fractured reservoir with wellbore storage and skin factor effects. Semi-analytical solutions of a two-dimensional diffusivity equation have been obtained by using successive applications of the Laplace and modified finite Fourier sine transforms. Both pseudo-steady-state and transient exchanges between the matrix and the fractures have been considered. Solutions are presented that can be used to generate type curves for pressure transient analysis or can be used as a forward model in parameter estimation. The presented analysis has applications in well testing of fractured aquifers and naturally fractured oil reservoirs with a gas cap.     

基于遗传算法的低压液压管路模型参数识别

为了合理预测伴随气泡和气穴的低压液压管路压力瞬态脉动,提出了用改进遗传算法对低压液压管路压力瞬态脉动模型进行参数辨识的新方法.给出了用来描述管路流动特性的瞬态脉动数学模型,建立了用来计算伴随气泡和气穴的液压管路瞬态下气泡体积和气穴体积的数学模型.构造了基于最小二乘法的适应度模型,探讨了遗传操作方式及算法终止准则,采用了算术交叉同线性逼近相结合的改进算术交叉算子进行交叉操作,给出了模型参数寻优的算法流程.实现了对低压液压管路压力瞬态脉动数学模型的参数识别,得到了参数优化后的低压液压管路压力瞬态脉动模型.仿真结果与实验数据的比较表明在低压液压管路瞬态模型中,用改进遗传算法来识别模型中的未知参数的方法是可行的、有效的.     

Peristaltic flow of Walter’s B fluid in endoscope

The peristaltic flow of a Walter’s B fluid in an endoscope is studied.The problem is modeled in a cylindrical coordinate system.The main theme of the present analysis is to study the endoscopic effects on the peristaltic flow of the Walter’s B fluid.To the best of the authors’ knowledge,no investigation has been made so far in the literatures to study the Walter’s B fluid in an endoscope.Analytical solutions are obtained using the regular perturbation method by taking δ as a perturbation parameter.The appro...     

Solution of the problem of a flow of rarefied gas of constant density around a semiinfinite plate by the integral diffusion method

Several papers [1–4] have proposed approximate diffusion models which can be used to examine the transport process in a rarefied gas where the mean free path is large and transport is not determined by the local gradient of the particular quantity.In this paper the integral diffusion model [2] is used to solve the problem of determination of the friction stress and velocity of a flow of an incompressible gas around a plane semi-infinite plate in the whole range of Knudsen numbers. The obtained solution is compared with published solutions and experimental data [9].     

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

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

非均质流固耦合介质轴对称动力问题时域解

将地基视为流固两相介质并考虑其非均质成层特性，推导了多层地基动力问题时域解．文中首先建立了一组解耦的两相介质动力控制方程；而后利用Ｌａｐｌａｃｅ－Ｈａｎｋｅｌ变换推导了单层地基象空间初参数解答；再利用初参数法及传递矩阵技术导出了任意多层地基瞬态解的一般解析算式．本文获得的解答可方便地退化为现有理想弹性介质的解答