Limit Models of Pore Space Structure of Porous Materials for Determination of Limit Pore Size Distributions Based on Mercury Intrusion Data

This paper proposes the application of capillary and chain random models of pore space structure for determination of limit pore diameter distributions of porous materials, based on the mercury intrusion curves. Both distributions determine the range in which the pore diameter distribution of the investigated material occurs and defines the degree of inaccuracy of the method based on the mercury intrusion data caused by the indeterminacy of the sample shape and its pore space architecture. We derived equations describing the quasi-static process of mercury intrusion into the porous layer and porous ball with a random chain pore space structure and analysed the influence of the model parameters on the mercury intrusion curves. It was shown that the distribution of link length in the chain model of the pore space, random location of chain capillaries in the sample and the length distribution of the capillaries do not influence significantly the intrusion process. Therefore, a simple model of the mercury intrusion into the layer is proposed in which chain links of the pore space have random diameters and constant length. This model is used as a basic model of the intrusion process into a sample of any shape and size and with homogeneous and isotropic chain pore space architecture. The thickness of the layer then represents the mean length of chain capillaries in the sample. It was also proved that the capillary and chain models of pore space architecture are limit models of the network model identified in this paper with the pore architecture of the investigated sample. This justifies the application of both models for determination of limit cumulative distributions of pore diameters in porous materials based on the mercury intrusion data.

     

WAG Displacements of Oil-Condensates Accounting for Hydrocarbon Ganglia

During two-phase flow in porous media, non-wetting phase is present simultaneously in states of mobile connected continuum and of trapped isolated ganglia. Mass exchange between these two parts of non-wetting phase is going on by dissolution and diffusion of component in the wetting phase, so, compositions of non-wetting phase in both parts are different. Nevertheless, the traditional mathematical model for two-phase multicomponent transport in porous media assumes the homogeneous distribution of each component in the overall non-wetting phase. New governing equations honouring ganglia of non-wetting phase are derived. They are successfully verified by a number of laboratory tests. Analytical model is developed for miscible water-alternate-gas (WAG) displacement of oil-condensates. The modelling shows that the significant amount of oil-condensate is left in porous media after miscible WAG, while the traditional model predicts that the miscible displacement results in the total sweep.     

Pore-Network Modeling of Isothermal Drying in Porous Media

In this paper we present numerical results obtained with a pore-network model for the drying of porous media that accounts for various processes at the pore scale. These include mass transfer by advection and diffusion in the gas phase, viscous flow in the liquid and gas phases and capillary effects at the liquid--gas interface. We extend our work by studying the effect of capillarity-induced flow in macroscopic liquid films that form at the pore walls as the liquid--gas interface recedes. A mathematical model that accounts for the effect of films on the drying rates and phase distribution patterns is presented. It is shown that film flow is a major transport mechanism in the drying of porous materials, its effect being dominant when capillarity controls the process, which is the case in typical applications.     

Two-Phase Flow in Heterogeneous Porous Media with Non-Wetting Phase Trapping

This article presents a mathematical model describing flow of two fluid phases in a heterogeneous porous medium. The medium contains disconnected inclusions embedded in the background material. The background material is characterized by higher value of the non-wetting-phase entry pressure than the inclusions, which causes non-standard behavior of the medium at the macroscopic scale. During the displacement of the non-wetting fluid by the wetting one, some portions of the non-wetting fluid become trapped in the inclusions. On the other hand, if the medium is initially saturated with the wetting phase, it starts to drain only after the capillary pressure exceeds the entry pressure of the background material. These effects cannot be represented by standard upscaling approaches based on the assumption of local equilibrium of the capillary pressure. We propose a relevant modification of the upscaled model obtained by asymptotic homogenization. The modification concerns the form of flow equations and the calculation of the effective hydraulic functions. This approach is illustrated with two numerical examples concerning oil–water and CO₂–brine flow, respectively.     

Review of Steady-State Two-Phase Flow in Porous Media: Independent Variables,Universal Energy Efficiency Map,Critical Flow Conditions,Effective Characterization of Flow and Pore Network

In many applications of two-phase flow in porous media, a wetting phase is used to displace through a network of pore conduits as much as possible of a non-wetting phase, residing in situ. The energy efficiency of this physical process may be assessed by the ratio of the flow rate of the non-wetting phase over the total mechanical power externally provided and irreversibly dissipated within the process. Fractional flow analysis, extensive simulations implementing the DeProF mechanistic model, as well as a recent retrospective examination of laboratory studies have revealed universal systematic trends of the energy efficiency in terms of the actual independent variables of the process, namely the capillary number, Ca, and the flow rate ratio, r. These trends can be cast into an energy efficiency map over the (Ca, r) domain of independent variables. The map is universal for all types of non-wetting/wetting phase porous medium systems. It demarcates the efficiency of steady-state two-phase flow processes in terms of pertinent system parameters. The map can be used as a tool for designing more efficient processes, as well as for the normative characterization of two-phase flows, as to the predominance of capillary or viscous effects. This concept is based on the existence of a unique locus of critical flow conditions, for which the energy efficiency takes locally maximum values. The locus shape depends on the physicochemical characteristics of the non-wetting phase/wetting phase/porous medium system, and it shows a significant mutation as the externally imposed flow conditions change the type of flow, from capillary- to viscosity-dominated. The locus can be approached by an S-type functional form in terms of the capillary number and the system properties (viscosity ratio, wettability, pore network geometry, etc.), suggesting that formative criteria can be derived for flow characterization in any system. A new, extended definition of the capillary number is also proposed that effectively takes into account the critical properties of all the system constituents. When loci of critical flow conditions pertaining to processes with different viscosity ratio in the same pore network, are expressed in terms of this true-to-mechanism capillary number, they collapse into a unique locus. In this context, a new methodology for the effective characterization of pore networks is proposed.     

Dynamic Network Model of Mobility Control in Emulsion Flow Through Porous Media

Modeling the flow of emulsion in porous media is extremely challenging due to the complex nature of the associated flows and multiscale phenomena. At the pore scale, the dispersed phase size can be of the same order of magnitude of the pore length scale and therefore effective viscosity models do not apply. A physically meaningful macroscopic flow model must incorporate the transport of the dispersed phase through the porous material and the changes on flow resistance due to drop deformation as it flows through pore throats. In this work, we present a dynamic capillary network model that uses experimentally determined pore-level constitutive relationships between flow rate and pressure drop in constricted capillaries to obtain representative transient macroscopic flow behavior emerging from microscopic emulsion flow at the pore level. A parametric analysis is conducted to study the effect of dispersed phase droplet size and capillary number on the flow response to both emulsion and alternating water/emulsion flooding in porous media. The results clearly show that emulsion flooding changes the continuous-phase mobility and consequently flow paths through the porous media, and how the intensity of mobility control can be tuned by the emulsion characteristics.     

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.     

Constitutive equation of a gas-filled two-phase medium

A set of equations governing the consolidation of a two-phase medium consisting of a porous elastic skeleton saturated with a highly compressible liquid (gas), is described. The homogenization method was utilized to deduce the equations. For the equivalent macroscopic medium, mass and momentum conservation equations and the flow equation of pore liquid are presented. Sample material constants were calculated using laboratory test results which were carried out at the Institute of Geotechnics, Technical University of Wroclaw.     

Poromechanics of freezing materials

When subjected to a uniform cooling below the freezing point a water-infiltrated porous material undergoes a cryo-deformation resulting from various combined actions: (i) the difference of density between the liquid water and the ice crystal, which results in the initial build-up of an in-pore pressure at the onset of crystallization; (ii) the interfacial effects arising between the different constituents, which eventually govern the crystallization process in connection with the pore access radius distribution; (iii) the drainage of the liquid water expelled from the freezing sites towards the air voids; (iv) the cryo-suction process, which drives liquid water towards the already frozen pores as the temperature further decreases; (v) the thermomechanical coupling between the solid matrix, the liquid water and the ice crystal. We work out a comprehensive theory able to encompass this whole set of actions. A macroscopic approach first provides the constitutive equations of freezing poroelastic materials, including the interfacial energy effects. This approach reveals the existence of a thermodynamic state function—namely the liquid saturation degree as a function of the temperature only. The macroscopic ice-dependent poroelastic properties are then upscaled from the knowledge of the elastic properties of the solid matrix, of the pore access radius distribution, and of the capillary curve. The theory is finally illustrated by analysing quantitatively the effects of the cooling rate and of the pore radius distribution upon the cryo-deformation of water-infiltrated porous materials. The theory succeeds in accounting for the experimentally observed shrinkage of embedded air voids, while predicting the partial melting of the ice already formed when the cooling suddenly stops.     

非饱和颗粒材料的多孔连续体有效压力与有效广义Biot应力

多孔连续体理论框架下的非饱和多孔介质广义有效压力定义和Bishop参数的定量表达式长期以来存在争议,这也影响了对与其直接相关联的非饱和多孔介质广义Biot有效应力的正确预测.基于随时间演变的离散固体颗粒-双联液桥-液膜体系描述的Voronoi胞元模型,利用由模型获得的非饱和颗粒材料表征元中水力-力学介观结构和响应信息,文章定义了低饱和度多孔介质局部材料点的有效内状态变量:非饱和多孔连续体的广义Biot有效应力和有效压力,导出了其表达式.所导出的有效压力公式表明,非饱和多孔连续体的有效压力张量为各向异性,它不仅对非饱和多孔连续体广义Biot有效应力张量的静水应力分量的影响呈各向异性,同时也对其剪切应力分量有影响.文章表明,非饱和多孔连续体中提出的广义Biot理论和双变量理论的基本缺陷在于它们均假定反映非混和两相孔隙流体对固相骨架水力-力学效应的有效压力张量为各向同性.此外,为定义各向同性有效压力张量和作为加权系数而引入的Bishop参数并不包含对非饱和多孔连续体中局部材料点水力-力学响应具有十分重要效应的基质吸力.所导出的非饱和多孔介质广义Biot有效应力和有效压力公式(包括反映有效压力...     

Digitally Reconstructed Porous Media: Transport and Sorption Properties

The basic aim of this work is to present a combination of techniques for the reconstruction of the porous structure and the study of transport properties in porous media. The disordered structure of porous systems like random sphere packing, Vycor glass and North Sea chalk, is represented by three-dimensional binary images. The random sphere pack is generated by a standard ballistic deposition procedure, while the chalk and the Vycor matrices by a stochastic reconstruction technique. The transport properties (Knudsen diffusivity, molecular diffusivity and permeability) of the resulting 3-dimensional binary domains are investigated through computer simulations. Furthermore, physically sound spatial distributions of two phases filling the pore space are determined by the use of a simulated annealing algorithm. The wetting and the non-wetting phases are initially randomly distributed in the pore space and trial-and-error swaps are performed in order to attain the global minimum of the total interfacial energy. The effective diffusivities of the resulting domains are then computed and a parametric study with respect to the pore volume fraction occupied by each phase is performed. Reasonable agreement with available data is obtained in the single- and multi-phase transport cases.     

Theoretical and experimental aspects of vacuum impregnation of porous media using transparent etched networks

Vacuum impregnation is a process method in which air and native solution are removed from the porous space of a given porous material and replaced by an external solution. Vacuum impregnation is divided into two steps: Firstly, the porous material is immersed in a liquid solution and exposed to subatmospheric pressure for a given time to ensure that air trapped in the porous materials will be removed; secondly, atmospheric pressure is re-established and the external solution penetrates the pore structure of the porous material. The objective of this study was to describe the hydrodynamic mechanisms involved in vacuum impregnation of porous materials as a function of capillary number and viscosity ratio. To achieve the objectives proposed in the present study, a transparent glass micromodel 7.7 cm × 7.4 cm was first constructed using the photolithographic technique. In addition, a stainless steel vacuum tank was built. The tank top was covered with a transparent reinforced glass plate. The whole system was connected to a vacuum pump, and a conventional video camera was adapted to record the experiments. Liquid saturation was determined through the image analysis process. Capillary number and viscosity ratio were determined for the drainage and imbibition processes. For the systems studied, we conclude that transport mechanisms ranged between stable displacement and capillary fingering during the vacuum step (drainage) while transport mechanisms ranged between continuous capillary and discontinuous capillary domains during the atmospheric step (imbibition). Earlier work indicated that our proposed process should be even more efficient for realistically large systems.     

Ensemble phase averaged equations for multiphase flows in porous media. Part 2: A general theory

Most models for multiphase flows in a porous medium are based on a straightforward extension of Darcy’s law, in which each fluid phase is driven by its own pressure gradient. The pressure difference between the phases is thought to be an effect of surface tension and is called capillary pressure. Independent of Darcy’s law, for liquid imbibition processes in a porous material, diffusion models are sometime used. In this paper, an ensemble phase averaging technique for continuous multiphase flows is applied to derive averaged equations and to examine the validity of the commonly used models. Closure for the averaged equations is quite complicated for general multiphase flows in a porous material. For flows with a small ratio of the characteristic length of the phase interfaces to the macroscopic length, the closure relations can be simplified significantly by an approximation with a second order error in this length ratio. This approximation reveals the information of the length scale separation obscured during an averaging process and leads to an equation system similar to Darcy’s law, but with additional terms. Based on interactions on phase interfaces, relations among closure quantities are studied.     

Droplet Imbibition into Paper Coating Layer: Pore-Network Modeling Simulation

Liquid penetration into thin porous media such as paper is often simulated using continuum-scale single-phase Darcy’s law. The underlying assumption was that a sharp invasion front percolates through the layer. To explore this ambiguous assumption and to understand the controlling pore-scale mechanisms, we have developed a dynamic pore-network model to simulate imbibition of a wetting phase from a droplet into a paper coating layer. The realistic pore structures are obtained using the FIB-SEM imaging of the coating material with a minimum resolution of 3.5 nm. Pore network was extracted from FIB-SEM images using Avizo software. Data of extracted pore network are used for statistically generating pore network. Droplet sizes are chosen in the range of those applicable in inkjet printing. Our simulations show no sharp invasion front exists and there is the presence of residual non-wetting phase. In addition, penetration of different sizes of droplets of different material properties into the pore network with different pore body and pore throat sizes are performed. We have found an approximately linear decrease in droplet volume with time. This contradicts the expected \(\sqrt{t}\)-behavior in vertical imbibition that is obtained using macroscopic single-phase Darcy’s law. With increase in flow rate, transition of imbibition invasion front from percolation-like pattern to a more sharper one with less trapping of non-wetting phase is also reported. Our simulations suggest that the single-phase Darcy’s law does not adequately describe liquid penetration into materials such as paper coating layer. Instead Richards equation would be a better choice.     

Motion equations of the dynamic theory of consolidation from the point of view of the theory of transient pipe flows

An elastic fluid-saturated porous medium is modeled as a bundle of parallel cylindrical tubes aligned in a direction parallel to the fluid movement. The pore space is filled with viscous compressible liquid. A cell model and the theory of transient pipe flow are used to derive one-dimensional governing equations in such media. All macroscopic constants in these equations are defined by the individual material constants of the fluid and solid. The interaction force includes an additional term not found in Biot's theory.