Pore Network Analysis of Resistivity Index for Water-Wet Porous Media

Pore network analysis is used to investigate the effects of microscopic parameters of the pore structure such as pore geometry, pore-size distribution, pore space topology and fractal roughness porosity on resistivity index curves of strongly water-wet porous media. The pore structure is represented by a three-dimensional network of lamellar capillary tubes with fractal roughness features along their pore-walls. Oil-water drainage (conventional porous plate method) is simulated with a bond percolation-and-fractal roughness model without trapping of wetting fluid. The resistivity index, saturation exponent and capillary pressure are expressed as approximate functions of the pore network parameters by adopting some simplifying assumptions and using effective medium approximation, universal scaling laws of percolation theory and fractal geometry. Some new phenomenological models of resistivity index curves of porous media are derived. Finally, the eventual changes of resistivity index caused by the permanent entrapment of wetting fluid in the pore network are also studied.Resistivity index and saturation exponent are decreasing functions of the degree of correlation between pore volume and pore size as well as the width of the pore size distribution, whereas they are independent on the mean pore size. At low water saturations, the saturation exponent decreases or increases for pore systems of low or high fractal roughness porosity respectively, and obtains finite values only when the wetting fluid is not trapped in the pore network. The dependence of saturation exponent on water saturation weakens for strong correlation between pore volume and pore size, high network connectivity, medium pore-wall roughness porosity and medium width of the pore size distribution. The resistivity index can be described succesfully by generalized 3-parameter power functions of water saturation where the parameter values are related closely with the geometrical, topological and fractal properties of the pore structure.     

Effect of Vibrations on Pore Fluid Distribution in Porous Media

Understanding the role of shuttle vibrations in pore fluid distribution is an essential task in the exploration of plant growth in root modules aboard space flights. Results from experimental investigations are reported in this paper on the distribution of immiscible fluid phases in glass beads under vibrations. Hexadecane, a petroleum compound immiscible with and lighter than water, was used in the experiments. The higher freezing point of Hexadecane (18 °C) allowed the solidification of the entrapped blobs in the presence of water in porous media, so that their size distribution can be obtained. van Genuchten function, commonly used to express moisture retention curves, is found to be an adequate fit for blob size distribution at residual saturation. The effect of vibrations on the fate (mobilization, stranding, or breakup) of a solitary ganglion in porous media was studied using a network model. A mobility criterion considering viscous, gravity, and capillary forces was developed to determine the fate of a solitary ganglion in a porous medium. It is concluded that the effect of vibrations is to increase the likelihood of breakup and mobilization of blobs entrapped in porous media at residual saturation. The pore fluid distributions after vibrations are less uniform than those before vibrations.     

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.     

Pore Size of Shale Based on Acyclic Pore Model

The transport properties inside a nanosize conduit, such as viscosity and density, deviate from those reported at unconfined conditions with an identical pressure and temperature. The deviation from the nominal value is well studied for the simple topology in nanofluidics. A straight circular tube is an example of a simple topology. The pore space inside the matrix of a shale formation is also a nanofluidic system, because its characteristic size is smaller than or equal to 100 nm, but it includes a complex structure. With this in mind, we determine the pore-throat and pore-body size distributions of different shales whose data are available in the literature. The main objective is to quantitatively distinguish the two sizes, whose importance is overlooked in the study of shale formations. The pore-throat size distribution is determined from mercury injection capillary pressure measurements, and the pore-body size distribution is derived from nitrogen adsorption–desorption. The acyclic pore model, which is physically representative of the pore space at the core scale, allows us to interpret the petrophysical measurements. Our study of different shales shows that the average pore-body size is usually larger than 20 nm; thus, there is no need to account for the pore proximity (confinement). The pore-throat size distributions of different shales usually fall below 20 nm, which entails the modification of a transport property that is relevant to the formation resistance against the flow.     

Condensation and isothermal water transfer in cement mortar Part I — Pore size distribution,equilibrium water condensation and imbibition

The pore size distribution of cement mortar is studied in relation to water sorption experiments with the help of mercury intrusion and nitrogen sorption. The importance of adsorbed water is pointed out. Isothermal imbibition experiments at four temperatures are presented. The temperature-dependence of the mass transfer coefficients is compared to the one predicted by the classical model. Significant discrepancies are noticed. On the basis of the knowledge of the pore structure, a modelisation of the transfer process at moderate water content is proposed. It particularly takes into account Knudsen's vapor diffusion and effects of the presence of a discontinuous capillary phase interacting with vapor diffusion.     

Pore Characterization Methodology by Means of Capillary Sorption Tests

In this study, a methodology is proposed for obtaining information about the porous structure of materials by analysing data supplied by a hydric absorption test. A model of vertical cylindrical pores is used for studying the variation of absorbed mass versus time. The proposed methodology requires an estimation of the numerical interval in which the radius distribution must be calculated and a minimum amount of experimental data. The experimental tests are developed with a system, designed and built in our laboratory, that allows capture the amount of data necessary to use the calculus methodology proposed. This methodology is based on the derivative of normalized mass in relation to the square root of time of the capillary absorption tests. The simplicity and low cost, in many cases, of the capillary absorption tests compared to other experimental techniques such as mercury porosimetry or SEM, make it feasible to design a simple methodology to obtain valuable information about the structure of a porous material, pore size distribution, and tortuosity factor.     

Capillary Hysteresis in Neutrally Wettable Fibrous Media: A Pore Network Study of a Fuel Cell Electrode

Hysteresis in the saturation versus capillary pressure curves of neutrally wettable fibrous media was simulated with a random pore network model using a Voronoi diagram approach. The network was calibrated to fit experimental air-water capillary pressure data collected for carbon fibre paper commonly used as a gas diffusion layer in fuel cells. These materials exhibit unusually strong capillary hysteresis, to the extent that water injection and withdrawal occur at positive and negative capillary pressures, respectively. Without the need to invoke contact angle hysteresis, this capillary behaviour is re-produced when using a pore-scale model based on the curvature of a meniscus passing through the centre of a toroid. The classic Washburn relation was shown to produce erroneous results, and its use is not recommended when modelling fibrous media. The important effect of saturation distribution on the effective diffusivity of the medium was also investigated for both water injection and withdrawal cases. The findings have bearing on the understanding of both capillarity in fibrous media and fuel cell design.     

Computing the Longtime Behaviour of NMR Propagators in Porous Media Using a Pore Network Random Walk Model

We present a pore network model combined with a random walk algorithm allowing the simulation of molecular displacement distributions in porous media as measured by NMR. A particular feature of this technique is the ability to probe the time evolution of these distributions. The objective is to predict the displacement behaviour for time intervals larger than the experimental observation time and explore the asymptotic dispersion regime at long times. Starting from 3D micro-CT images, we computed the variance of displacement distributions of water molecules in a Fontainebleau sand and found very good agreement of the time evolution of the variance with experimental data, without fitting parameter. The model confirms a weak superdispersion in the asymptotic regime. In addition, we conclude that, since pore network models do not take into account small scale features of the porous medium (e.g., surface roughness and grain shape), the origin of the observed superdispersion is mainly due to the topology and geometry of the porous medium.     

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.     

Numerical Simulation of Porosity in Cements

The pores in cementitious materials, their sizes and connectivity have an important influence on the durability of concrete. Several microstructural models have been developed to simulate the three-dimensional pore network in cement pastes. In this article, microstructures with the $\upmu $ μ ic model are compared with experimental results. It is seen that despite having a resolution for the capillary pores very close to reality, the experimentally observed breakthrough diameter from mercury intrusion is much lower than the values obtained by applying an algorithm of mercury intrusion to the simulated microstructure. The effect of some of the most important input parameters on the pore sizes in the simulated microstructure explored. The phenomenon which seems best able to explain this discrepancy is that C–S–H is not in fact a phase with a smooth surface as represented in microstructural models, but a phase which grows as needles into the pore space, leading to very small water-filled capillary pores from quite young ages. The results demonstrate it will be extremely challenging to represent the porosity of real microstructures in microstructural models on the scale of hundreds of microns necessary to study macroscopic transport.     

Pore-Level Modeling of Gas and Condensate Flow in Two- and Three-Dimensional Pore Networks: Pore Size Distribution Effects on the Relative Permeability of Gas and Condensate

We present a mechanistic model of retrograde condensation processes in two- and three-dimensional capillary tube networks under gravitational forces. Condensate filling-emptying cycles in pore segments and gas connection–isolation cycles are included. With the pore-level distribution of gas and condensate in hand, we determine their corresponding relative permeabilities. Details of pore space and displacement are subsumed in pore conductances. Solving for the pressure field in each phase, we find a single effective conductance for each phase as a function of condensate saturation. Along with the effective conductance for the saturated network, the relative permeability for each phase is calculated. Our model porous media are two- and three-dimensional regular networks of pore segments with distributed size and square cross-section. With a Monte Carlo sampling we find the optimum network size to avoid size effects and then we investigate the effect of network dimensionality and pore size distribution on the relative permeabilities of gas and condensate.     

Simulation of capillary pressure curves using bond correlated site percolation on a simple cubic network

A stochastic approach to network modelling has been used to simulate quasi-static immiscible displacement in porous media. Both number-based and volume-based network saturation results were obtained. Number-based results include: number-based saturation curves for primary drainage, secondary imbibition and secondary drainage, fluid distribution data, and cluster trapping history. Using pore structure data of porous media, it is possible to convert the number-based curves to capillary pressure — saturation relationships. Pore size distribution functions and pore shapes which are thought to closely represent Berea sandstone samples were used to predict the capillary curves. The physical basis of these calculations is a one-to-one correspondence between the cumulative node and bond index fractions in the network analysis, and the cumulative number-based distributions of pore body and pore throat diameters, respectively. The oil-water capillary pressure curve simulated for primary drainage closely resembles those measured experimentally. The agreement between the simulated and the measured secondary imbition and secondary drainage curves is less satisfactory.     

Pore Network Modelling of Electrical Resistivity and Capillary Pressure Characteristics

In order to model petrophysical properties of hydrocarbon reservoir rocks, the underlying physics occurring in realistic rock pore structures must be captured. Experimental evidence showing variations of wetting occurring within a pore, and existence of the so-called 'non-Archie' behaviour, has led to numerical models using pore shapes with crevices (for example, square, elliptic, star-like shapes, etc.). This paper presents theoretical derivations and simulation results of a new pore space network model for the prediction of petrophysical properties of reservoir rocks. The effects of key pore geometrical factors such as pore shape, pore size distribution and pore co-ordination number (pore connectivity) have been incorporated into the theoretical model. In particular, the model is used to investigate the effects of wettability and saturation history on electrical resistivity and capillary pressure characteristics. The petrophysical characteristics were simulated for reservoir rock samples. The use of the more realistic grain boundary pore (GBP) shape allows simulation of the generic behaviour of sandstone rocks, with various wetting scenarios. The predictions are in close agreement with electrical resistivity and capillary pressure characteristics observed in experiments.     

Permeability of Pore Networks Under Compaction

The characteristic pore length fixes the scale of permeability of a porous medium. For pore networks undergoing strong random compaction, this length becomes singular at transition porosities, revealing a change in the microstructure of the porespace controlling the transport. Nodal balances and lattice Boltzmann simulations of flow in pore networks under compaction show that the scaling between permeability and porosity changes near the transition porosities. Simulation results are compared with experimental permeability data from transparent two-dimensional micromodels of networks decorated with the same pore size distribution. Permeability?Cporosity data of media undergoing smooth compaction is well described by a single power law. Under strong compaction, however, the scaling between permeability and porosity is possible by traits only, the scaling exponent changes notably at given transition porosities. These behaviors are common to a wealth of permeability?Cporosity data thus far unexplained.     

Spontaneous Inertial Imbibition in Porous Media Using a Fractal Representation of Pore Wall Rugosity

Considering the separable phenomena of imbibition in complex fine porous media as a function of timescale, it is noted that there are two discrete imbibition rate regimes when expressed in the Lucas–Washburn (L–W) equation. Commonly, to account for this deviation from the single equivalent hydraulic capillary, experimentalists propose an effective contact angle change. In this work, we consider rather the general term of the Wilhelmy wetting force regarding the wetting line length, and apply a proposed increase in the liquid–solid contact line and wetting force provided by the introduction of surface meso/nanoscale structure to the pore wall roughness. An experimental surface pore wall feature size regarding the rugosity area is determined by means of capillary condensation during nitrogen gas sorption in a ground calcium carbonate tablet compact. On this nano size scale, a fractal structure of pore wall is proposed to characterize for the internal rugosity of the porous medium. Comparative models based on the Lucas–Washburn and Bosanquet inertial absorption equations, respectively, for the short timescale imbibition are constructed by applying the extended wetting line length and wetting force to the equivalent hydraulic capillary observed at the long timescale imbibition. The results comparing the models adopting the fractal structure with experimental imbibition rate suggest that the L–W equation at the short timescale cannot match experiment, but that the inertial plug flow in the Bosanquet equation matches the experimental results very well. If the fractal structure can be supported in nature, then this stresses the role of the inertial term in the initial stage of imbibition. Relaxation to a smooth-walled capillary then takes place over the longer timescale as the surface rugosity wetting is overwhelmed by the pore condensation and film flow of the liquid ahead of the bulk wetting front, and thus to a smooth walled capillary undergoing permeation viscosity-controlled flow.