Wavelet upscaling based on piecewise bilinear approximation of the permeability field

A method for upscaling of permeability in heterogeneous porous media is presented. The upscaled field takes the form K = e ^Y , where Y, in two dimensions, is a piecewise bilinear function. The method is tested on a number of random permeability fields, with different integral scale/correlation length and variance. The numerical results show that this method conserves much more of the heterogeneous fingering than classical block-based upscaling methods, e.g., geometric mean.     

Fractal porous media II: Geometry of porous geological structures

Some geological structures are analysed and found to be fractal. An interesting feature is the very large range of scales involved; the spreading dimension is also measured for some of them. The consequences of these measurements on the analysis of transport processes in porous media are presented - the existence of fractal structures multiplies the variety of actual porous media.     

Effective permeability of cracked unsaturated porous materials

This study focuses on a theoretical estimation of the effective permeability of unsaturated cracked porous media. The closed-form flow solution around and in a superconductive crack, embedded in an infinite porous matrix under a far-field condition, is recalled first. Then the solution of flow around a completely unsaturated (empty) crack that is considered as an obstruction against the flow is determined. The flow solution for partially saturated crack in special configurations is obtained by superposition of the two basic solutions for superconductive and empty cracks. The contribution of an unsaturated crack, with a given saturation degree, to the effective permeability is estimated by using dilute upscaling scheme. Numerical results obtained by Finite Elements Method, are in good agreement with the theoretical results for weak crack densities but show the additional effect of cracks interaction for higher densities.     

The use of renormalization for calculating effective permeability

There is a need in the numerical simulation of reservoir performance to use average permeability values for the grid blocks. The permeability distributions to be averaged over are based on samples taken from cores and from logs using correlations between permeabilities and porosities and from other sources. It is necessary to use a suitable effective value determined from this sample. The effective value is a single value for an equivalent homogeneous block. Conventionally, this effective value has been determined from a simple estimate such as the geometric mean or a detailed numerical solution of the single phase flow equation.If the permeability fluctuations are small then perturbation theory or effective medium theory (EMT) give reliable estimates of the effective permeability. However, for systems with a more severe permeability variation or for those with a finite fraction of nonreservoir rock all the simple estimates are invalid as well as EMT and perturbation theory.This paper describes a real-space renormalization technique which leads to better estimates than the simpler methods and is able to resolve details on a much finer scale than conventional numerical solution. Conventional simulation here refers to finite difference (or element) techniques for solving the single phase pressure equation. This requires the pressure and permeability at every grid point to be stored. Hence, these methods are limited in their resolution by the amount of data that can be stored in core. Although virtual memory techniques may be used they increase computer time. The renormalization method involves averaging over small regions of the reservoir first to form a new averaged permeability distribution with a lower variance than the original. This pre-averaging may be repeated until a stable estimate is found. Examples are given to show that this is in excellent agreement with computationally more expensive numerical solution but significantly different from simple estimates such as the geometric mean.     

The effect of anisotropic permeability on free convective boundary layer flow in porous media

The effect of an anisotropic permeability on thermal boundary layer flow in porous media is studied. The convective flow is induced by a vertical, uniformly heated surface embedded in a fluid-saturated medium. A leading-order boundary layer theory is presented. It is shown that the thickness of the resulting boundary layer flow is different from that obtained in an isotropic porous medium. In general, an anisotropic permeability induces a fluid drift in the spanwise direction, the strength of which depends on the precise nature of the anisotropy. Conditions are found which determine whether or not the boundary layer flow is three-dimensional.     

Experimental study of two-phase flow in three sandstones. II. Capillary pressure curve measurement and relative permeability pore space capillary models

By means of the porous plate method and mercury porosimetry intrusion tests, capillary pressure curves of three different sandstones were measured. The testing results have been exploited jointly with three relative permeability models of the pore space capillary type (Burdine’s model type), these models are widely used and in rather distinct fields. To do so, capillary pressure has been correlated to saturation degree using six of the most popular relations encountered in the literature. Model predictions were systematically compared to the experimentally measured relative permeabilities presented in the first part of this work. Comparison indicated that the studied models underestimate the water relative permeability and over-estimate that of the non-wetting phase. Moreover, this modeling proves to be unable to locate the significant points that are the limits of fields of saturation where the variation of the relative permeabilities becomes consequent. We also showed that, if pore structure is modeled as a “bundle of capillary tubes”, model predications are independent of the capillary pressure curve measuring method.     

Higher-order correction of effective permeability of heterogeneous isotropic formations of lognormal conductivity distribution

Steady flow of an incompressible fluid takes place in a porous formation of spatially variable hydraulic conductivityK. The latter is regarded as a lognormal stationary random space function and Y=ln(K/K _G ), whereK _G is the geometric mean ofK, is characterized by its variance ² and correlation scale I. Exact results are known for the effective conductivityK _eff in one- and two-dimensional flows. In contrast, only a first-order term in a perturbation expansion in ² has been derived exactly for the three-dimensional flow. A conjecture has been made in the past onK _eff for any ², but it was not yet proved exactly. This study derived the exact nonlinear correction, i.e. the termO(⁴) ofK _eff, which is found to be the one resulting from the conjecture, strengthening the confidence in it. It is also shown that the self-consistent approximation leads to the exact results for one-dimensional and two-dimensional flows, but underestimates the nonlinear correction ofK _eff for in the three-dimensional case.     

The effective permeability of a heterogeneous porous medium

The effective (single-phase) permeability of an (infinite) heterogeneous porous medium is studied using a formalism of Green's functions. We give formal expressions for it in the form of a series expansion involving the microscopic random-permeability field many-body correlation functions of higher and higher order.The particular case of a log-normal medium of infinite extent is studied using field-theoretical methods. Using partial series resummation techniques, we derivea formula up to all orders in the local correlations which was first reckoned by many authors by means of a first-order calculation. The formula — which remains an approximation — works whatever the dimensionality of the space, and gives the following simple estimate for the effective permeability in 3 D:K _eff=k ^1/33. The method is general and the approximations can be systematically improved on when more complex situations are studied.Roman Letters D number of dimensions of the space in which the flow takes place - f(r) body force field,N - f(q) Fourier-transformed body-force field, Nm³ - G ₀(r, r) Green's function of the Laplace operator, m^–1 - g(k,r, r) velocity propagator before averaging, m^–1 - G(r, r) velocity propagator after averaging, m^–1 - j(r) a scalar dimensionless field - k(r) local value of the permeability at point r, m² - K _eff effective permeability - K _g geometric average of the local permeability, m² - l typical size of the averaging volume, m - L characteristic length of the porous medium or of the reservoir, m - L(r, r) projection operator, m^–2 - M(r, r) scattering operator, m^–3 - p(r) local value of the pressure, Nm^–2 - p(k,r, r) pressure propagator before averaging, m^–1 - P(r, r) pressure propagator after averaging, m^–1 - r position vector, m - r modulus of vectorr, m - unit vector pointing in the direction ofr - q Fourier wave vector, m^–1 - q modulus of the Fourier wave-vectorq, m^–1 - unit vector pointing in the direction ofq - projector over vector - 1 unit tensor - X(r) a local random variable - ¯X(r) volume averaged local random variable - X (r) ensemble averaged local random variable - V large-scale averaging volume, m³ - Z(j) generating functional of a random field - Z(r,j) modified generating functional of a random field - Z normalization factor Greek Letters ₀ average value of the logarithm of the permeability - (r) fluctuation of the logarithm of permeability at pointr - viscosity of the fluid, Nt/m² - (r–r) two-point correlation function of the fluctuations of the logarithm of the permeability - _k correlation length of the permeability correlation function, m - _u correlation length of the velocity correlation function, m     

A variational approach for the deformation of a saturated porous solid. A second-gradient theory extending Terzaghi's effective stress principle

Summary The principle of virtual power is used to derive the equilibrium field equations of a porous solid saturated with a fluid, including second density-gradient effects; the intention is the elucidation and extension of the effective stress principle of Terzaghi and Fillunger. In the context of a first density-gradient theory for a saturated solid we interpret the porewater pressure as a Lagrange multiplier in the expression for the deformation energy, assuring that the saturation constraint is verified. We prove that this saturation pressure is distributed among the constituents according to their respective volume fraction (Delesse law) only if they are both true density-preserving. We generalize the Delesse law to the case of compressible constituents. If a material-dependent effective stress contribution is to arise, it is, in general, nonvanishing simultaneously in both the solid and fluid constituents. Moreover, saturation pressure, effective stresses and compressibility constitutive equations determine the exchange volume forces. In a theoretical formulation without non-isotropic strain measures, second density-gradient effects must be incorporated, not only to accommodate for the equilibrium-solid-shear stress and the interaction among neighboring solid-matrix pores, but also to describe internal capillarity effects. The earlier are accounted for by a dependence of the thermodynamic energy upon the density-gradient of the solid, while the latter derives from a mixed density-gradient dependence. Examples illustrate the necessity of these higher gradient effects for properly posed boundary value problems describing the mechanical behaviour of the disturbed rock zone surrounding salt caverns. In particular, we show that solid second-gradient strains allow for the definition of the concept of static permeability, which is distinct from the dynamic Darcy permeability. Received 1 February 1999; accepted for publication 9 March 1999     

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

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

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

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

基于分形理论的高强混凝土动态损伤本构关系

基于高强混凝土（HSC）试块在SHPB冲击实验中的分形损伤演化规律,推导了HSC的分形损伤变量表达式,标定了HSC裂纹的分形维数范围。然后参考ZWT模型,并结合HSC实验过程中的应变率相关性、动态损伤特性及近似恒应变率,推导了分形损伤演化的HSC动态损伤本构方程。采用4组应变率工况下的C60、C80混凝土应力-应变曲线对本构方程进行验证,理论曲线和实验曲线吻合较好。

     

基于分维理论的爆炸焊接界面形貌的定量表征

以爆炸焊接界面形貌为研究对象,利用分形理论研究其界面特征。运用三维超景深显微镜获得界面形貌图像,利用Matlab图像分析技术对界面形貌进行二值化处理,基于分形理论计算图像分维值以及多重分维谱。由分维值及多重分形谱的物理意义可知,轮廓分维值是对界面不规则程度宏观度量,而多重分形谱能定量表征界面的起伏程度及高度分布。因而,通过对界面结构进行分形几何分析,可实现界面形貌定量表征,弥补当前定性分析的不足。     

基于分形的捷联惯组历次测试数据混合插值算法

针对捷联惯组历次测试数据小样本、非等间隔、非线性的问题,提出了一种基于分形插值的三次混合插值算法.通过第一次分形插值保证原始测试数据的变化趋势;通过第二次样条函数插值保证了插值的准确性,实现原始测试数据等间隔化;通过第三次分段线性插值,扩大样本容量,同时保证了原有测试序列的统计特性不变.实例分析表明,该算法很好的实现了对捷联惯组历次测试数据的等间隔处理和样本容量扩大,为捷联惯组历次测试数据小样本建模分析提供良好的基础.     

Effects of wall permeability on turbulence

In order to understand the effects of the wall permeability on turbulence near a porous wall, flow field measurements are carried out for turbulent flows in a channel with a porous bottom wall by a two-component particle image velocimetry (PIV) system. The porous media used are three kinds of foamed ceramics which have almost the same porosity (0.8) but different permeability. It is confirmed that the flow becomes more turbulent over the porous wall and tends to be turbulent even at the bulk Reynolds number of Re_b=1300 in the most permeable wall case tested. Corresponding to laminar to turbulent transition, the magnitude of the slip velocity on the porous wall is found to increase drastically in a narrow range of the Reynolds number. To discuss the effects of the wall roughness and the wall permeability, detailed discussions are made of zero-plane displacement and equivalent wall roughness for porous media. The results clearly indicate that the turbulence is induced by not only the wall roughness but the wall permeability. The measurements have also revealed that as Re_b or the wall permeability increases, the wall normal fluctuating velocity near the porous wall is enhanced due to the effects of the wall permeability. This leads to the increase of the turbulent shear stress resulting in higher friction factors of turbulence over porous walls.     

Relationship between space charge and effective pore size of nanoporous electrode in electric double-layer capacitor

The space-charge distributions of electric double-layer capacitors (EDLCs), which are energy storage devices, were examined with the pulsed electroacoustic (PEA) method. It was found that the experimental results could be influenced by the reflection and penetration of sound waves when the space-charge distributions of EDLCs were measured with the PEA method. This is because EDLCs have a five-layer structure consisting of three materials (aluminum, cellulose, and activated carbon). We calculated the reflect...     

On the use of conventional cocurrent and countercurrent effective permeabilities to estimate the four generalized permeability coefficients which arise in coupled,two-phase flow

In the case of coupled, two-phase flow of fluids in porous media, the governing equations show that there are four independent generalized permeability coefficients which have to be measured separately. In order to specify these four coefficients at a specific saturation, it is necessary to conduct two types of flow experiments. The two types of flow experiments used in this study are cocurrent and countercurrent, steady-state permeability experiments. It is shown that, by taking this approach, it is possible to define the four generalized permeability coefficients in terms of the conventional cocurrent and countercurrent effective permeabilities for each phase. It is demonstrated that a given generalized phase permeability falls about midway between the conventional, cocurrent effective permeability for that phase, and that for the countercurrent flow of the same phase. Moreover, it is suggested that the conventional effective permeability for a given phase can be interpreted as arising out of the effects of two types of viscous drag: that due to the flow of a given phase over the solid surfaces in the porous medium and that due to momentum transfer across the phase 1-phase 2 interfaces in the porous medium. The magnitude of the viscous coupling is significant, contributing at least 15% to the total conventional cocurrent effective permeability for both phases. Finally, it is shown that the nontraditional generalized permeabilities which arise out of viscous coupling effects cannot equal one another, even when the viscosity ratio is unity and the surface tension is zero.     

Influence of mixture sensitivity and pore size on detonation velocities in porous media

Experiments have been carried out to determine the dependence of the detonation velocity in porous media, on mixture sensitivity and pore size. A detonation is established at the top end of a vertical tube and allowed to propagate to the bottom section housing the porous bed, comprised of alumina spheres of equal diameter (1–32 mm). Several of the common detonable fuels were tested at atmospheric initial pressure. Results indicate the existence of a continuous range of velocities with change in Φ, spanning the lean and the rich propagation limits. For all fuels in a given porous bed, the velocity decreases from a maximum value at the most sensitive mixture near Φ≈1 (minimum induction length), toV/V _CJ≈0.3 at the limits. A decrease in pore size brings about a reduction inV/V _CJ and a narrowing of the detonability range for each fuel. For porous media comprised of spherical particles, it was possible to correlate the velocity data corresponding to a variety of different mixtures and for a broad range of particle sizes, using the following empirical expression:V/V _CJ=[1–0.35 log(d _c /d _p)]±0.1. The critical tube diameterd _c is used as a measure of mixture sensitivity andd _p denotes the pore diameter. An examination of the phenomenon at the composition limits, suggests that wave failure is controlled by a turbulent quenching mechanism.     

基于空间分布的重力场持续适配能力评估方法

重力匹配导航需要在重力场特征丰富的区域进行,现有的重力场适配性评估方法多针对单条测线或局部区域,未考虑到重力特征的空间分布。为了评估重力场持续适配能力,提出了将重力场适配半径作为评估重力场特征区域空间适配性能指标参数。从适配半径定义出发,阐述了其物理意义内涵,证明了其与匹配测线覆盖全区域的最小辐射半径的关系,并根据适配半径不同分布形态,推导了具体计算方法。适配性评估实例表明,中国南海重力场适配半径为95 n mile,说明重力匹配导航技术在南海范围具有较好适用性。适配半径直接反映了特征区域提供持续重力匹配导航的能力,可作为重力匹配导航系统工程应用重力场持续适配能力评估和航迹规划的依据。