多孔介质渗透率的一种新分形模型

多孔介质的渗流特性是油气藏工程、地下水资源利用、高放废物深地质处置等实际工程领域的热门研究问题．基于分形理论及多孔介质由一束面积大小不等的椭圆形毛细管组成的假设,本文建立了流体在分形多孔介质中渗流时的绝对渗透率及相对渗透率的分形渗透率模型．结果表明,绝对渗透率是最大和最小孔隙面积、分形维数、形状因子ε的函数,且当ε =1时,本文模型可以简化成Yu与Cheng模型;而非饱和多孔介质的相对渗透率与饱和度和多孔介质微结构参数有关．将本文提出的渗透率分形模型预测与实验测量数据及其他模型结果进行对比,显示它们整体吻合很好．     

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

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

基于渗透率连续变化的低渗透多孔介质非线性渗流模型

基于低渗透多孔介质渗透率的渐变理论,确定了能精确描述低渗透多孔介质渗流特征的非线性运动方程,并通过实验数据拟合．验证了非线性运动方程的有效性。非线性渗流速度关于压力梯度具有连续-阶导数,方便于工程计算;由此建立了低渗透多孔介质的单相非线性径向渗流数学模型,并巧妙采用高效的Douglas-Jones预估一校正有限差分方法求得了其数值解。数值结果分析表明：非线性渗流模型为介于拟线性渗流模型和达西渗流模型之间的一种中间模型或理想模型,非线性渗流模型和拟线性渗流模型均存在动边界;拟线性渗流高估了启动压力梯度的影响,使得动边界的移动速度比实际情况慢得多;非线性越强,地层压力下降的范围越小,地层压力梯度越陡峭,影响地层压力的敏感性减弱,而影响地层压力梯度的敏感性增强。     

多孔介质输运性质的分形分析研究进展

首先对多孔介质输运性质的传统实验测量、解析分析和数值模拟计算研究进展作了扼要的评述.然后,着重综述采用分形理论和方法研究多孔介质输运性质分析解的理论、方法和所取得的进展.最后,指出采用分形理论和方法有可能解决其它尚未解决的有关多孔介质输运性质的若干课题和方向.      

An overview on nonlinear porous flow in low permeability porous media

This paper gives an overview on nonlinear porous flow in low permeability porous media, reveals the microscopic mechanisms of flows, and clarifies properties of porous flow fluids. It shows that, deviating from Darcy's linear law, the porous flow characteristics obey a nonlinear law in a low-permeability porous medium, and the viscosity of the porous flow fluid and the permeability values of water and oil are not constants. Based on these characters, a new porous flow model, which can better describe low permeability reservoir, is established. This model can describe various patterns of porous flow, as Darcy's linear law does. All the parameters involved in the model, having definite physical meanings, can be obtained directly from the experiments.     

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.     

Study of the effect of capillary pressure on the permeability of porous media embedded with a fractal-like tree network

In this paper, the analytical expressions for permeability of (both saturated and unsaturated) porous media embedded with a fractal-like tree network are presented based on fractal theory and technique when the capillary pressure is taken into account. Both the dimensionless effective permeability and the relative permeability of the composites, which are defined as porous media embedded with a fractal-like tree network in this work, are derived and found to be a function of saturation, the capillary pressure and microstructural parameters of the networks. The relative permeabilities predicted by the present fractal model are compared with the available experimental data and a fair agreement between them is found.     

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

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

孤立波与多孔介质结构物的非线性相互作用

基于精确至Ｏ（εμ＾２,μ＾４）的多孔介质无压渗流模型方程和均匀流体质波动的Ｂｏｕｓｓｉｎｅｓｑ方程,本文对孤立波与多孔介质结构物的相互作用了较系统的数值实验。控制方程采用基于有限差分方程离散,在时域上采用了预估－校正方法进行了时间积分。在求解演化方程的过程中,引入“内迭代”过程实现流体域和多孔介质交界面的连接条件。结果表明孤立波在多孔介质上的反射与在不可渗透的界面上的反射类似,形成反向的孤立波但     

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

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

Low Pressure Gas Percolation Characteristic in Ultra-low Permeability Porous Media

Low pressure gas percolation characteristic in ultra-low permeability porous media is investigated in this article through core flow experiments. The results show that the wall-slip layer covers more than 10% of the average porous channel radius on account of minimum pore size when the permeability is below 0.1 × 10^?3μ m ² order, and seepage behavior is contrasted to that in mid-high permeability pore media. When the gas pressure is not high enough, the flow regime turns into transitional flow instead of slip flow, and nonlinear relationship between the measured gas permeability and the reciprocal of average pressure exists. The gas measuring permeability experiment would be influenced by the non-linear relationship. If Klinkenberg-corrected method is applied to speculate the equivalent liquid permeability, the extrapolated value will become less or minus. Simultaneously, actual gas flow velocity at the outlet is beyond the calculated value with Klinkenberg formula. A new gas seepage model based on the general slip boundary condition is derived from the homogenization technique in this article. At last the flow model is examined to be suitable for representing the gas flow behavior in ultra-low permeability media and estimating the absolute permeability from single-point, steady-states measurements.     

基于多变量小样本的渗流代理模型及产量预测方法

多孔介质渗流过程中存在的多尺度、多变量、多物理场耦合的非线性渗流问题给复杂渗流机理的表征及数学模型求解提出了巨大的挑战, 综合考虑地下多孔介质耦合渗流过程中关键力学问题的渗流模型往往需要在计算效率和计算精度之间权衡. 近年来, 基于油田多数据的渗流代理模型为高效求解多变量非线性渗流问题提供了思路, 而渗流代理模型在实际油田中的应用往往由于记录不全, 操作不当等因素受到小样本数据的限制. 针对这一问题, 本文提出了一种基于地质?油藏?工艺的多数据小样本渗流代理模型的产量预测方法. 通过填补缺失值, 独热编码分类数据, 数据对数化及标准化等一系列数据预处理方法, 形成了油田产量预测数据库; 经过随机劈分数据集、十折交叉验证, 测试了三种渗流代理模型的预测效果. 结果表明, 三种代理模型的决定系数均超过0.8, 模型预测结果与实际数据较为吻合; 对于小样本多变量的油田数据, 合适的数据预处理方法对模型预测效果影响显著; 经过数据标准化后, 随机森林算法表现最好, 能快速准确预测石油产量(均方误差0.12, 决定系数0.87).      

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.     

Investigation of Generalized Relative Permeability Coefficients for Electrically Assisted Oil Recovery in Oil Formations

Evaluation of relative permeability coefficients is one of the key steps in reliable simulation of two-phase flow in porous media. An extensive body of work exists on evaluation of these coefficients for two-phase flow under pressure gradient. Oil transport under an applied electrical gradient in porous media is also governed by the principles of two-phase flow, but is less understood. In this paper, relative permeability coefficients under applied electric field are evaluated for a specific case of two- phase fluid flow in water-wet porous media, where the second fluid phase is oil. It is postulated that the viscous drag on the oil phase, exerted by the electro-osmotic flow of the water phase, is responsible for the transport of oil in the absence of a pressure gradient. Reliable prediction of the flow patterns necessitates accurate representation and determination of the relative permeability coefficients under the electrical gradient. The contribution of each phase to the flow is represented mathematically, and the relative permeability coefficients are evaluated through electro-osmotic flow measurements conducted on oil bearing rock cores.     

Effective Solution through the Streamline Technique and HT-Splitting for the 3D Dynamic Analysis of the Compositional Flows in Oil Reservoirs

In this article we suggest a new phenomenological mathematical model for the groundwater transport of colloid particles through porous media which is able to describe some significant effects experimentally observed but not captured within the framework of the classic approach. Our basic idea is to consider both the pure water and the colloid suspension as two thermodynamic phases. Using the network models of porous media, we simulated numerically the transport process at the pore-scale. By averaging the result derived, we have obtained the relative permeabilities for both phases, the percolation threshold for suspension flow, and the effective suspension viscosity. Due to specific laws of colloid particles repartition between various classes of pores, the relative permeability of suspension happens to be a highly nonlinear function of saturation, very far from the diagonal straight line. This determines a difference between the macroscale phase velocities. The suspension velocity is shown to be higher than that of water in major cases, only if the colloid particles are not too large. The suggested model predicts and describes in a closed form the effect of colloid transport facilitation observed experimentally.     

孔隙渗流对环面复层含油轴承润滑性能的影响

含油轴承基体中油液的渗流特性对轴承油膜润滑性能影响显著. 以不同孔隙率分布的环面复层含油轴承为研究对象,利用达西定律建立复层含油轴承基体中流体渗流的控制方程,在极坐标下建立环面复层含油轴承系统渗流润滑模型,研究复层环面副系统中油膜压力分布规律,分析轴承结构参数及孔隙渗流行为对油膜润滑性能的影响. 结果表明:复层含油轴承的流体动压力主要发生在环形摩擦面间,从摩擦界面到轴承底面,流体压力逐渐由外缘向圆心部位传导,流体动压力作用面积逐渐增大,压力峰值逐渐降低;随着倾角增大,摩擦界面间的油膜动压效应增强,油膜润滑性能变好;随着表层渗透率或厚度减小,摩擦界面间的油膜的渗流效应减弱,油膜润滑性能提高;与普通单层含油轴承相比,含致密表层的复层含油轴承能降低整体孔隙率,防止过多轴承间隙油液渗入多孔介质,提高轴承润滑性能. 研究工作为明晰环面复层含油轴承渗流行为及润滑机理提供一定理论依据.      

On the Permeability of Fractal Tube Bundles

The permeability of a porous medium is strongly affected by its local geometry and connectivity, the size distribution of the solid inclusions, and the pores available for flow. Since direct measurements of the permeability are time consuming and require experiments that are not always possible, the reliable theoretical assessment of the permeability based on the medium structural characteristics alone is of importance. When the porosity approaches unity, the permeability?Cporosity relationships represented by the Kozeny?CCarman equations and Archie??s law predict that permeability tends to infinity and thus they yield unrealistic results if specific area of the porous media does not tend to zero. The aim of this article is the evaluation of the relationships between porosity and permeability for a set of fractal models with porosity approaching unity and a finite permeability. It is shown that the tube bundles generated by finite iterations of the corresponding geometric fractals can be used to model porous media where the permeability?Cporosity relationships are derived analytically. Several examples of the tube bundles are constructed, and the relevance of the derived permeability?Cporosity relationships is discussed in connection with the permeability measurements of highly porous metal foams reported in the literature.     

Three-Dimensional Lattice Boltzmann Simulations of Single-Phase Permeability in Random Fractal Porous Media with Rough Pore–Solid Interface

Single-phase permeability k has intensively been investigated over the past several decades by means of experiments, theories and simulations. Although the effect of surface roughness on fluid flow and permeability in single pores and fractures as well as networks of fractures was studied previously, its influence on permeability in a random mass fractal porous medium constructed of pores of different sizes remained as an open question. In this study, we, therefore, address the effect of pore–solid interface roughness on single-phase flow in random fractal porous media. For this purpose, we apply a mass fractal model to construct porous media with a priori known mass fractal dimensions \(2.579 \le D_{\mathrm{m}} \le 2.893\) characterizing both solid matrix and pore space. The pore–solid interface of the media is accordingly roughened using the Weierstrass–Mandelbrot approach and two parameters, i.e., surface fractal dimension \(D_{\mathrm{s}}\) and root-mean-square (rms) roughness height. A single-relaxation-time lattice Boltzmann method is applied to simulate single-phase permeability in the corresponding porous media. Results indicate that permeability decreases sharply with increasing \(D_{\mathrm{s}}\) from 1 to 1.1 regardless of \(D_{\mathrm{m}}\) value, while k may slightly increase or decrease, depending on \(D_{\mathrm{m}}\), as \(D_{\mathrm{s}}\) increases from 1.1 to 1.6.