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.     

A Geometrical Pore Model for Estimating the Microscopical Pore Geometry of Soil with Infiltration Measurements

This paper presents a simple geometrical pore model designed to relate characteristic pore radii of the porous network of soils with macroscopic infiltration parameters. The model composed of a stack of spherical hollow elements is described with two radii values: the pore access radius and the actual pore radius. The model was compared to cylindrical pore models and its mathematical consistency was assessed. Soil sorptivity S and the second parameter A of the Philip infiltration equation (1957), have been determined by numerically simulated infiltration. A diagram and an empirical relation have been set in order to relate the pore access and pore radii to the infiltration parameters S and A. The consistency of the model was validated by comparing the predicted sorptivity and hydraulic conductivity values, with the widely used unsaturated soil hydraulic functions (van Genuchten, 1980). The model showed good agreement with experimental infiltration data, and it is therefore concluded that the use of a model with two radii improves the relation between microscopic pore size and macroscopic infiltration parameters.     

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.     

Pore water pressures in secondary consolidation

This paper describes the work carried out to investigate the pore pressures occurring in secondary consolidation. A theoretical approach and an experimental technique was developed in order to conduct the study.By considering compression to occur only due to water leaving the soil it was possible to derive an expression for the dissipation of pore pressure in the secondary phase. By further simplified assumptions which are based on experimental observations, the above general solution was reduced to a simple formula which predicted the observed behaviour of pore water pressures during secondary consolidation.     

Correction to: Pore Scale Numerical Modelling of Geological Carbon Storage Through Mineral Trapping Using True Pore Geometries

Transport in Porous Media -     

Pore structures and transport properties of sandstone

To quantitatively analyze the macroscopic properties of the flow in porous media by means of the continuum approach, detailed information (velocity and pressure fields) on the microscopic scale is necessary. In this paper, the numerical solution for incompressible, Newtonian flow in a diverging-converging representative unit cell (RUC) is presented. A new solution procedure for the problem is introduced. A review of the accuracy of the computational method is given.Nomenclature A _ff ^* area of entrance and exit of RUC - A _fs ^* interfacial area between the fluid and solid phases - d throat diameter of RUC (m) - D pore diameter of RUC (m) - i, j unit vector for RUC - L ^* wave length of a unit cell - L _p pore length of RUC (m) - L _t throat length of RUC (m) - n unit outwardly directed vector for the fluid phase - p ^* fluid pressure - ^* cross-sectional mean pressure - _en ^* entrance cross-sectional mean pressure - Re _d Reynolds number - x ^*, r^* cylindrical coordinates - u ^*, v^* velocity - u _cl ^* centerline velocity - _d mean velocity at the throat of RUC (m/s) - _D mean velocity at the large segment of RUC (m/s) Greek viscosity coefficient (Ns/m²) - _p excess momentum loss factor defined in (4.1) - fluid density (kg/m³) - ^* stream function - ^* vorticity - dimensionless circulation defined in (2.7) Symbols - the mean value - * dimensionless quantities     

Pore space network characterization with sub-voxel definition

Morphological measurements in 3D for pore space characterization (connectivity pore-body/throat classification, shape factors, virtual fluid intrusion) are based on computed intensive digital-thinning operations for skeletonization and medial axis extraction from 3D digital images. We present an alternative method that is measurably faster and allows sub-voxel definition of the pore space network. The method allows extracting—based on morphological considerations only—the centered and shortest stream-lines—i.e., the paths—to follow in order to go through the pore space from one given point to another and to exit. In addition the method penalizes long and narrow pore-throats in favor of short stubby/ones—i.e., it has a built-in exemplification capacity. It exploits well-established mathematical methods successfully applied in medical endoscopy.     

Multi-scale Image-Based Pore Space Characterisation and Pore Network Generation: Case Study of a North Sea Sandstone Reservoir

Transport in Porous Media - In this paper, we examine the pore space geometry and topology of a North Sea sandstone reservoir rock based on multi-scale scanning electron microscopy. The reservoir...     

Coalescence of Pore Columns by Domain Switching

The present paper studies the coalescence of pore columns in ferroelectric ceramics driven by back and forth domain switching under cyclic electric field. A finite element method that incorporates mass transfer capacity is formulated to simulate the evolution of point defects subjected to the kinetics of pore surface diffusion and domain wall migration. The merge of point defects provides a mechanism for the vacancy agglomeration that leads to the formation of large pores or microcracks.The project supported by the National Natural Science Foundation of China (10332020 and 10121202), and the National 973 Project (2004CB619304).     

Three-Phase Effective Contact Angle in a Model Pore

We consider a model pore (2D) in which a sharp interface between two fluids contact a third fluid which wets the solid boundary. If the configuration is capillary dominated, the geometry can be determined analytically in terms of the effective contact angle. This angle depends not only on the interfacial tensions, but also on the capillary pressures. However, if the height of the cusp formed by the wetting fluid is much smaller than the pore width, the effective contact angle is a simple function of the interfacial tensions. It turns out to be the same function as in the case of an undeformed wetting layer of molecular thickness. The analytical expression for the effective contact angle has been confirmed by a numerical technique, known as the lattice-Boltzmann method. This method, in turn, has been validated with Neumann's law for the three-phase contact angles.     

Capillary Imbibition and Pore Characterisation in Cement Pastes

The kinetics of capillary imbibition in ordinary Portland cement pastes has been studied experimentally and theoretically. Nuclear magnetic resonance stray field imaging (STRAFI) has been used to record water concentration profiles for various ingress times. The profiles follow a t law and thus a master curve can be formed using the Boltzmann transformation. The distribution of pore sizes within the sample as measured by NMR cryoporometry shows a prominent peak at 100Å. A computer model of the pore structure was developed consisting of a lattice of interconnecting pores with a size distribution consistent with the cryoporometry results. The Hagen–Poiseuille law was used to describe the kinetics of the water in this pore structure. The best agreement between the computer simulations and the experimental master curve was obtained by using a narrower range of pore sizes than indicated by the cryoporometry results.     

Pore Scale Modeling of Rate Effects in Imbibition

We use pore scale network modeling to study the effects of flow rate and contact angle on imbibition relative permeabilities. The model accounts for flow in wetting layers that occupy roughness or crevices in the pore space. Viscous forces are accounted for by solving for the wetting phase pressure and assuming a fixed conductance in wetting layers. Three-dimensional simulations model granular media, whereas two-dimensional runs represent fracture flow.We identify five generic types of displacement pattern as we vary capillary number, contact angle, and initial wetting phase saturation: flat frontal advance, dendritic frontal advance, bond percolation, compact cluster growth, and ramified cluster growth. Using phase diagrams we quantify the range of physical properties under which each regime is observed. The work explains apparently inconsistent experimental measurements of relative permeability in granular media and fractures.     

Bubble Migration Velocity in a Uniform Pore Network

Transport in Porous Media - Gas bubbles can be generated naturally or introduced artificially in porous media. Gas bubble migration through porous media governs the rate of gas emission to the...     

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.

     

Pore Volume and Porosity Changes under Uniaxial Strain Conditions

Expressions for the changes that occur in the pore volume and the porosity of a porous rock, due to changes in the pore pressure, overburden stress, and temperature, are derived within the context of the linearised theory of poroelasticity. The resulting expressions are compared to the commonly used equations proposed by Palmer and Mansoori, and it is shown that their expressions are consistent with the exact expressions if their factor f is identified with \((1+\nu )/3(1-\nu )\), where \(\nu \) is the Poisson’s ratio of the porous rock. Finally, the first derivation is given, within the context of the fully coupled linearised theory of poroelasticity, that under uniaxial strain, the partial differential equation that governs the evolution of the pore pressure is a pure diffusion equation, with a total compressibility term that (exactly) equals the sum of the fluid compressibility and the uniaxial pore volume compressibility.     

Stochastic Pore Network Generation from 3D Rock Images

Pore networks can be extracted from 3D rock images to accurately predict multi-phase flow properties of rocks by network flow simulation. However, the predicted flow properties may be sensitive to the extracted pore network if it is small, even though its underlying characteristics are representative. Therefore, it is a challenge to investigate the effects on flow properties of microscopic rock features individually and collectively based on small samples. In this article, a new approach is introduced to generate from an initial network a stochastic network of arbitrary size that has the same flow properties as the parent network. Firstly, we characterise the realistic parent network in terms of distributions of the geometrical pore properties and correlations between these properties, as well as the connectivity function describing the detailed network topology. Secondly, to create a stochastic network of arbitrary size, we generate the required number of nodes and bonds with the correlated properties of the original network. The nodes are randomly located in the given network domain and connected by bonds according to the strongest correlation between node and bond properties, while honouring the connectivity function. Thirdly, using a state-of-the-art two-phase flow network model, we demonstrate for two samples that the rock flow properties (capillary pressure, absolute and relative permeability) are preserved in the stochastic networks, in particular, if the latter are larger than the original, or the method reveals that the size of the original sample is not representative. We also show the information that is necessary to reproduce the realistic networks correctly, in particular the connectivity function. This approach forms the basis for the stochastic generation of networks from multiple rock images at different resolutions by combining the relevant statistics from the corresponding networks, which will be presented in a future publication.     

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.     

Determination of Nanoparticle Macrotransport Coefficients from Pore Scale Processes

Nanoparticle transport in porous media is modeled using a hierarchical set of differential equations corresponding to pore scale and macroscale. At the pore scale, movement and interaction of a single particle with a solid matrix is modeled using the advection–dispersion–sorption equation. A single nanoparticle entering the space encounters viscous, diffusion and surface forces. Surface forces (electrostatic and van der Waals forces) between nanoparticles and mineral grains appear as sorption propensity on solid matrix boundary condition. These local events are then transformed into a macroscale continuum by imposing periodic boundary conditions for contiguous unit cells representing porous media and using a scheme of moment analysis. At the macroscale, propagation and retention of particles are characterized by three position-independent coefficients: mean nanoparticle velocity vector \({\bar{\mathbf{U}}}^*\), macroscopic dispersion coefficient \({\bar{\mathbf{D}}}^*\), and mean nanoparticle retention rate constant \({\bar{K}}^*\). The modeling results are validated with a set of nanoparticle transport tests in porous microchips. We also present simulations of realistic porous media, where an actual image of sandstone samples is processed into binary tones. The representative unit cells are constructed from the resulting binary image by searching for areas within the sample with maximum similarities to the whole sample in terms of porosity and specific surface area, which are found to show strong correlations with the resulting \({\bar{\mathbf{U}}}^*\) and \({\bar{K}}^*\), respectively.