Homogenization of Wall-Slip Gas Flow Through Porous Media

The permeability of reservoir rocks is most commonly measured with an atmospheric gas. Permeability is greater for a gas than for a liquid. The Klinkenberg equation gives a semi-empirical relation between the liquid and gas permeabilities. In this paper, the wall-slip gas flow problem is homogenized. This problem is described by the steady state, low velocity Navier–Stokes equations for a compressible gas with a small Knudsen number. Darcy's law with a permeability tensor equal to that of liquid flow is shown to be valid to the lowest order. The lowest order wall-slip correction is a local tensorial form of the Klinkenberg equation. The Klinkenberg permeability is a positive tensor. It is in general not symmetric, but may under some conditions, which we specify, be symmetric. Our result reduces to the Klinkenberg equation for constant viscosity gas flow in isotropic media.     

Does Klinkenberg's Law Survive Upscaling?

This work deals with the large-scale mathematical modelling of flow of gas at low pressure in porous media. At the pore scale, this type of flow is characterised by a wall-slip effect, which at the sample scale entails a dependence of permeability upon gas pressure. This latter property is described by Klinkenberg's law. The goal of the present work is to examine the robustness of this law, by determining whether it is still verified on a large-scale: upscaling is thus applied, starting with Klinkenberg's law at the local scale. A Klinkenberg's flow of gas in a two-constituent composite porous medium is considered, and the constituents are firstly assumed to be homogeneous. The cases of low and of high permeability contrast are successively examined. Upscaling is performed using the homogenisation method of multiple scale expansions. In both cases, the large-scale permeability tensor differs from its liquid counterpart. Except in the particular case of equal Klinkenberg factors, Klinkenberg's law is not verified at low permeability contrast. At high permeability contrast, the large-scale gas permeability verifies Klinkenberg's law. The case of heterogeneous constituents is then examined. It is shown that the large-scale permeability differs from its liquid counterpart, but it does not verify Klinkenberg's law.     

Heat transfer effects on particle motion under rarefied conditions

Gas–solid flows occurring on very small spatial scales (of the order of micro and nanometres) are of great relevance in a number of industrial applications. It is currently not well established how particle motion and filtration are affected by non-isothermal conditions at such scales. Furthermore, when the particle size is comparable to the mean free path of the gas, rarefaction effects become important. In the present work we investigate the effects of heat transfer and non-isothermal conditions on the motion of small particles in rarefied flow. For that purpose, a suitable framework is developed here as a generic multiphase DNS method for rarefied flows. The resulting model is valid for low particle Reynolds number flows, irrespective of the Biot number, and for particle Knudsen numbers up to unity in unbounded flow. Using this model, we show that there is different settling behaviour of particles with an internal heat source in rarefied and continuum cases of the carrier gas respectively. It is shown that the chances for thermal levitation and/or lifting up of a particle due to buoyancy effects are significantly reduced under rarefied conditions.     

Investigation of the rarefied gas flow around a circular cylinder in stationary and oscillatory regimes

Numerically, on the basis of the Krook kinetic equation, the rarefied gas flow around a circular cylinder is investigated in stationary and oscillatory regimes. The flows around a rotating cylinder and a cylinder with a nonuniformly heated surface are considered. The Knudsen numbers at which the lift force acting on the rotating cylinder changes sign are calculated. It is shown that at low Knudsen numbers a lift force acts on the nonuniformly heated cylinder.     

考虑多重运移机制耦合页岩气藏压裂水平井数值模拟

页岩作为典型的微纳尺度多孔介质,游离气与吸附气共存,传统的达西定律已无法准确描述气体在页岩微纳尺度的运移规律.基于双重介质模型和离散裂缝模型构建页岩气藏分段压裂水平井模型,其中基岩中考虑气体的黏性流、Knudsen 扩散以及气体在基岩孔隙表面的吸附解吸,吸附采用Langmuir等温吸附方程;裂缝中考虑黏性流和Knudsen扩散,在此基础上建立基岩-裂缝双重介质压裂水平井数学模型并采用有限元方法对模型进行求解.结果表明,基岩固有渗透率越小,表面扩散和Knudsen扩散的影响越大,反之则越小;人工裂缝的性质包括条数、开度、半长以及间距,主要影响压裂水平井生产早期,随着人工裂缝参数值的增加,压裂水平井产能增加,累产气量也越大.其次,页岩气藏压裂诱导缝和天然裂缝的发育程度对页岩气藏的产能有很大的影响,水平井周围只有人工裂缝,周围天然裂缝不开启或不发育时,页岩气藏的水平井的产能较低.      

Knudsen’s Permeability Correction for Tight Porous Media

Various flow regimes including Knudsen, transition, slip and viscous flows (Darcy’s law), as applied to flow of natural gas through porous conventional rocks, tight formations and shale systems, are investigated. Data from the Mesaverde formation in the United States are used to demonstrate that the permeability correction factors range generally between 1 and 10. However, there are instances where the corrections can be between 10 and 100 for gas flow with high Knudsen number in the transition flow regime, and especially in the Knudsen’s flow regime. The results are of practical interest as gas permeability in porous media can be more complex than that of liquid. The gas permeability is influenced by slippage of gas, which is a pressure-dependent parameter, commonly referred to as Klinkenberg’s effect. This phenomenon plays a substantial role in gas flow through porous media, especially in unconventional reservoirs with low permeability, such as tight sands, coal seams, and shale formations. A higher-order permeability correlation for gas flow called Knudsen’s permeability is studied. As opposed to Klinkenberg’s correlation, which is a first-order equation, Knudsen’s correlation is a second-order approximation. Even higher-order equations can be derived based on the concept used in developing this model. A plot of permeability correction factor versus Knudsen number gives a typecurve. This typecurve can be used to generalize the permeability correction in tight porous media. We conclude that Knudsen’s permeability correlation is more accurate than Klinkenberg’s model especially for extremely tight porous media with transition and free molecular flow regimes. The results from this study indicate that Klinkenberg’s model and various extensions developed throughout the past years underestimate the permeability correction especially for the case of fluid flow with the high Knudsen number.     

复杂微通道气体流动的格子Boltzmann模拟

构建了一个模拟复杂微通道内气体流动的多松弛格子Boltzmann模型。该模型采用动力学曲面滑移边界,考虑了微尺度效应和努森层影响。此外,为了更准确地描述微通道内气体的滑移速度,在模型中引入孔隙局部Kn数来代替平均Kn数。之后采用Poiseuille流对模型进行验证,模拟结果与用直接模拟蒙特卡洛方法和分子模拟结果吻合较好,证明了该模型模拟微通道内处于滑移区和过渡区气体流动的有效性。最后,采用该模型模拟多孔介质内气体渗流过程。结果表明,随着孔隙平均Kn数的增加,多孔介质内的高渗区域增加,且优先从小孔隙中开始增加,这是由于小孔隙中微尺度效应更加明显,相对大孔隙流动阻力更小所致。     

页岩气和致密砂岩气藏微裂缝气体传输特性

页岩气和致密砂岩气藏发育微裂缝,其开度多在纳米级和微米级尺度且变化大,因此微裂缝气体传输机理异常复杂.本文基于滑脱流动和努森扩散模型,分别以分子之间碰撞频率和分子与壁面碰撞频率占总碰撞频率的比值作为滑脱流动和努森扩散的权重系数,耦合这两种传输机理,建立了微裂缝气体传输模型. 该模型考虑微裂缝形状和尺度对气体传输的影响. 模型可靠性用分子模拟数据验证.结果表明:(1)模型能够合理描述微裂缝中所有气体传输机理,包括连续流动,滑脱流动和过渡流动;(2)模型能够描述不同开发阶段,微裂缝中各气体传输机理对传输贡献的逐渐变化过程;(3)微裂缝形状和尺度影响气体传输,相同开度且宽度越大的微裂缝,气体传输能力越强,且在高压和微裂缝大开度的情况下表现更明显.      

Relaxation processes in hypersonic rarefied binary gas mixture flow around a cylinder

Using the direct simulation Monte Carlo method, the hypersonic flow of a binary gas mixture around a cylinder is investigated over a wide rarefaction range: from an almost continuum regime (at the Knudsen number Kn = 0.01) to free-molecular flow. The effect of a small admixture of heavy diatomic particles in a light gas flow on the relaxation processes near the cylinder and the heat flux is studied.     

Bulk equations and Knudsen layers for the regularized 13 moment equations

The order of magnitude method offers an alternative to the Chapman-Enskog and Grad methods to derive macroscopic transport equations for rarefied gas flows. This method yields the regularized 13 moment equations (R13) and a generalization of Grad’s 13 moment equations for non-Maxwellian molecules. Both sets of equations are presented and discussed. Solutions of these systems of equations are considered for steady state Couette flow. The order of magnitude method is used to further reduce the generalized Grad equations to the non-linear bulk equations, which are of second order in the Knudsen number. Knudsen layers result from the linearized R13 equations, which are of the third order. Superpositions of bulk solutions and Knudsen layers show good agreement with DSMC calculations for Knudsen numbers up to 0.5.      

Numerical experiments on the gas flow between eccentric rotating cylinders

A numerical experiment was carried out on the gas flow field between two eccentric cylinders, one of which is rotating. Attention was paid to the presence of separated recirculating regions from the continuum to the rarefied regimes. The direct simulations were performed by means of a Monte Carlo (DSMC) method and bi‐polar co‐ordinates were adopted. The calculations were relative to isothermal walls at the same temperature. Streamlines and velocity profiles were evaluated as functions of the Knudsen number, of the Mach number and of the geometric parameters. The gas considered was argon. Copyright © 2000 John Wiley & Sons, Ltd.     

Intense evaporation of a gas from a two-dimensional periodic surface

Evaporation (or condensation) of a gas is said to be intense when the normal component of the velocity of the gas in the Knudsen layer has a value of the order of the thermal velocity of a molecule, c_T=(2kT/m)^1/2. In this case the distribution function of the molecules with respect to their velocities in the Knudsen layer differs from the equilibrium (Maxwellian) value by its own magnitude. As a result of this, over the thickness of the Knudsen layer the macroparameters also vary by their own magnitudes. So in order to obtain the correct boundary conditions for the Euler gas dynamic equations, it is necessary to solve the nonlinear Boltzmann equation in the Knudsen layer. The problem of obtaining such boundary conditions for the case of a plane surface was considered in [1–11]. In the present study this problem is solved for a two-dimensional periodic surface in the case when the dimensions of the inhomogeneities are of the order of the mean free path of the molecules and the inhomogeneities have a rectangular shape. The flow in the Knudsen layer becomes two-dimensional, and this leads to a considerable complication of the solution of the problem.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 2, pp. 132–139, March–April, 1985.In conclusion the author would like to express his gratitude to V. A. Zharov for his valuable advice, and also V. S. Galkin, M. N. Kogan, and N. K. Makashev for discussion of the results obtained.     

Comprehensive Seepage Simulation of Fluid Flow in Multi-scaled Shale Gas Reservoirs

Shale gas seepage behaviour is a multi-field/-scale problem and makes transient pressure analysis a very challenging task. Non-Darcy flow in nanopores is prominent due to the broken of continuity hypothesis. Slippage effect and Knudsen diffusion are two important seepage mechanisms in nanopores, while recent studies show surface diffusion is another important transporting mechanism on surface of nanopores. Porous kerogen system contains large amounts of dissolved gas, which should not be overlooked. In this study, a comprehensive mathematical model was established by pseudo-quadruple porosity medium conception, coupling the effects of slippage flow, Knudsen diffusion, surface diffusion, ad-/desorption and gas transferring from kerogen to nanopore system, while fluid flow in fractures/macropores is described by Darcy’s law. Transient pressure behaviours of a multiple fractured horizontal well in box-shaped shale gas reservoir were studied, with nine possible flow regimes divided and parameters sensitivity analysed. Adsorbed constant and dissolved constant were defined to reflect the amount of adsorbed gas and dissolved gas, respectively. Research shows that adsorbed gas and dissolved gas are two important gas storage forms, neither of which should be neglected. The study can not only help us understand fluid flow mechanisms in nanopores from microscopic perspective, but enable us to analyse production performance and determine key operational parameters from macroscopic perspective.     

Diffusiophoresis of an aerosol particle in a binary gas mixture

The diffusion force and rate are calculated for the diffusiophoresis of a spherical particle in a binary gas mixture by solving the gas–kinetic equations. Two schemes of diffusiophoresis are considered: constant–pressure diffusion and diffusion of one mixture component through the other fixed component. The problem is solved by the integral–momentum method at arbitrary Knudsen numbers. Diffuse scattering of the gas molecules on the particle surface is assumed. The Lorentzian and Rayleigh models of a binary gas mixture are considered. The dependences of the force and rate of diffusiophoresis on the Knudsen number and the other determining parameters are analyzed. The results obtained are compared with well–known experimental data.     

Contribution of the Dusty Gas Model to Permeability/Diffusion Tests on Partially Saturated Clay Rocks

Gas transfer experiments on claystone and numerical simulations have been conducted to enhance the knowledge of gas transport in nuclear waste repositories in the Callovo-Oxfordian clay formation in Bure, France. Laboratory Gas transfer experiments were performed with a specific device dedicated to very low permeability measurement (10^?23 to 10^?20 m²). Experiments were performed on both dry and close to saturation claystone. The Dusty Gas Model, based on multi-component gas transfer equations with Knudsen diffusion, was used to describe the experimental results. The parameters obtained are the effective permeability, the Knudsen diffusion (Klinkenberg effect) and molecular diffusion coefficients and the porosity accessible to gas. Numerical simulations were carried with various boundary conditions and for different gases (helium vs hydrogen) and were compared with experiments to test the reliability of the model parameters and to better understand the mechanisms involved in clays close to saturation. The numerical simulation fitted the experimental data well whereas simpler models cannot describe the complexity of the Knudsen/Klinkenberg effects. Permeabilities lie between 10^?22 and 10^?20 m². Claystones close to saturation have an accessible porosity to gas transfer that is lower than 0.1?C1% of the porosity. Analysis of the Klinkenberg effect suggests that this accessible pore network should be made of 50?C200?nm diameter pores. It represents pore networks accessible at capillary pressure lower than 4?MPa.     

Effect of rotational degrees of freedom of molecules on energy flux in attenuated gas

On the basis of model kinetic equations a solution is obtained by a numerical method for the flow of attenuated gas around a sphere. The effect of rotational degrees of freedom on the energy flux to the body is investigated. Values of the ratio between the energy flux Q and its free-molecular value Q* for monatomic and diatomic gases are compared; for the comparison, the dimensionless temperature of the body, the gas velocity at infinity, and the law of viscosity must be the same in the two cases. For sufficiently cold bodies (when the body temperature is below the equilibrium temperature for a diatomic gas) the difference between Q/Q* for monatomic and diatomic gases is insignificant. For a diatomic gas when the body temperature is close to equilibrium, the ratio Q/Q* is found to have a nonmonotonic dependence on the Knudsen force.Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 5, pp. 119–124, September–October, 1977.     

Weakly rarefied gas flows between nonuniformly heated parallel plates

The aim of this research is to establish the validity of the predictions of the theory of slow nonisothermal flows, to study the limits of applicability (with respect to the Knudsen number) of the conclusions reached and to determine the effect of the Knudsen layers on these flows on the basis of a numerical investigation of slow nonisothermal weakly rarefied gas flow in a plane infinite channel with weakly nonequilibrium heating of the walls and a finite wall temperature difference. The gas flow is described by a relaxation transport equation. The results obtained show how quickly, as the Knudsen number decreases, the solutions of the transport equation outside the Knudsen layers tend to the solution of the equations of gas dynamics of slow nonisothermal flows (and not to the solution of the Navier-Stokes equations).Translated from Izvestiya Akademii Nauk SSSR, Mekhanika Zhidkosti i Gaza, No. 1, pp. 115–121, January–February, 1988.     

基于稀薄效应的微气体径向轴承稳态性能

针对微气体轴承给出参考努森数的定义,根据空气不同温度时的黏度,得到参考努森数的分布范围;考虑气体稀薄效应,给出基于Burgdorfer一阶滑移速度边界的微气体径向轴承润滑Reynolds方程的修正形式; 采用有限差分法求解修正的Reynolds方程,得到不同参考努森数$Kn_0$, 轴承数以及轴颈偏心率情况下轴承的压力分布、无量纲承载能力及偏位角. 数值分析表明:随气体稀薄程度的增强,气体径向轴承的压力明显降低,无量纲承载力降低,而轴承偏位角增大. 当偏心率小于0.6时,轴承偏位角变化平缓,受$Kn_0$数的影响不明显. 当轴承数较小时,气体稀薄程度对轴承的无量纲承载力、偏位角影响较小.      

超薄膜磁头滑块气动力特性

采用有限差分法对广义润滑方程进行数值求解,计算出计算机磁头滑块压强场的分布情况。分析、研究了其稳态和动态气动力特性,并将计算结果分别与求解一阶、二阶修正雷诺方程所得到的结果进行了比较,得到如下三个结论:(1)当飞行高度很小,飞行速度较低时,必须采用广义润滑方程进行磁头滑块的气动力计算,(2)与广义润滑方程结果比较,求解一阶修正雷诺方程所得到的计算结果总是偏高,而求解二阶修正雷诺方程所得到的计算结果总是偏低。此外,还解决了大压缩数下数值失稳问题,使得压缩数可以计算到120万,足以适应任何实际工程的需要。