Pore-scale modeling is essential in understanding and predicting flow and transport properties of rocks. Generally, pore-scale modeling is dependent on imaging technologies such as Micro Computed Tomography (micro-CT), which provides visual confirmation into the pore microstructures of rocks at a representative scale. However, this technique is limited in the ability to provide high resolution images showing the pore-throats connecting pore bodies. Pore scale simulations of flow and transport properties of rocks are generally done on a single 3D pore microstructure image. As such, the simulated properties are only representative of the simulated pore-scale rock volume. These are the technological and computational limitations which we address here by using a stochastic pore-scale simulation approach. This approach consists of constructing hundreds of 3D pore microstructures of the same pore size distribution and overall porosity but different pore connectivity. The construction of the 3D pore microstructures incorporates the use of Mercury Injection Capillary Pressure (MICP) data to account for pore throat size distribution, and micro-CT images to account for pore body size distribution. The approach requires a small micro-CT image volume (7–19 mm3) to reveal key pore microstructural features that control flow and transport properties of highly heterogeneous rocks at the core-scale. Four carbonate rock samples were used to test the proposed approach. Permeability calculations from the introduced approach were validated by comparing laboratory measured permeability of rock cores and permeability estimated using five well-known core-scale empirical model equations. The results show that accounting for the stochastic connectivity of pores results in a probabilistic distribution of flow properties which can be used to upscale pore-scale simulated flow properties to the core-scale. The use of the introduced stochastic pore-scale simulation approach is more beneficial when there is a higher degree of heterogeneity in pore size distribution. This is shown to be the case with permeability and hydraulic tortuosity which are key controls of flow and transport processes in rocks.
more »
« less
Unexpected scaling of interstitial velocities with permeability due to polymer retention in porous media
Polymer retention from the flow of a polymer solution through porous media results in substantial decrease of the permeability; however, the underlying physics of this effect is unknown. While the polymer retention leads to a decrease in pore volume, here we show that this cannot cause the full reduction in permeability. Instead, to determine the origin of this anomalous decrease in permeability, we use confocal microscopy to measure the pore-level velocities in an index-matched model porous medium.We show that they exhibit an exponential distribution and, upon polymer retention, this distribution is broadened yet retains the same exponential form. Surprisingly, the velocity distributions are scaled by the inverse square root of the permeabilities. We combine experiment and simulation to show these changes result from diversion of flow in the random porous-medium network rather than reduction in pore volume upon polymer retention.
more »
« less
- Award ID(s):
- 2011754
- PAR ID:
- 10499992
- Publisher / Repository:
- American Physical Society
- Date Published:
- Journal Name:
- Physical Review Fluids
- Volume:
- 6
- Issue:
- 8
- ISSN:
- 2469-990X
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Despite their potential for the scalable production of mixed-matrix membranes (MMMs), the MMMs prepared by the polymer-modification-enabled in situ metal–organic framework formation (PMMOF) process showed a considerable reduction in gas permeability as the filler loading increased. It was hypothesized that a correlation existed between the decrease in permeability and the change in the properties of the polymer, such as free volume and chain flexibility, upon in situ MOF formation. Herein, we aim to address the permeability reduction by using a cross-linked polyimide (6FDA-DAM:DABA (3:2)). It was found the degree of cross-linking affected not only the properties of the polymer, but also the in situ formation of the ZIF-8 filler particles in the cross-linked polymer. The proper degree of cross-linking resulted in suppressing C3H6 permeability reduction, suggesting a possible strategy to overcome the issue of PMMOF. The swelling of the polymer followed by chain rearrangement during the PMMOF, as well as the structural rigidity of the polymer, were found to be critical in mitigating permeability reduction.more » « less
-
This study discusses turbulent suspension flows of non-Brownian, non-colloidal, neutrally buoyant and rigid spherical particles in a Newtonian fluid over porous media with particles too large to penetrate and move through the porous layer. We consider suspension flows with the solid volume fraction $${{\varPhi _b}}$$ ranging from 0 to 0.2, and different wall permeabilities, while porosity is constant at 0.6. Direct numerical simulations with an immersed boundary method are employed to resolve the particles and flow phase, with the volume-averaged Navier–Stokes equations modelling the flow within the porous layer. The results show that in the presence of particles in the free-flow region, the mean velocity and the concentration profiles are altered with increasing porous layer permeability because of the variations in the slip velocity and wall-normal fluctuations at the suspension-porous interface. Furthermore, we show that variations in the stress condition at the interface significantly affect the particle near-wall dynamics and migration toward the channel core, thereby inducing large modulations of the overall flow drag. At the highest volume fraction investigated here, $${{\varPhi _b}}= 0.2$$ , the velocity fluctuations and the Reynolds shear stress are found to decrease, and the overall drag increases due to the increase in the particle-induced stresses.more » « less
-
Abstract The coupled hydrology and mechanics of soft porous materials (such as clays, hydrogels, membranes, and biofilms) is an important research area in several fields, including water and energy technologies as well as biomedical engineering. Well‐established models based on poromechanics theory exist for describing these coupled properties, but these models are not adapted to describe systems with more than one characteristic length scale, that is, systems that contain both macropores and micropores. In this paper, we expand upon the well‐known Darcy‐Brinkman formulation of fluid flow in two‐scale porous media to develop a “Darcy‐Brinkman‐Biot” formulation: a general coupled system of equations that approximates the Navier‐Stokes equations in fluid‐filled macropores and resembles the equations for poroelasticity in microporous regions. We parameterized and validated our model for systems that contain either plastic (swelling clay) or elastic microporous regions. In particular, we used our model to predict the permeability of an idealized siliciclastic sedimentary rock as a function of pore water salinity and clay content. Predicted permeability values are well described by a single parametric relation between permeability and clay volume fraction that agrees with existing experimental data sets. Our novel formulation captures the coupled hydro‐chemo‐mechanical properties of sedimentary rocks and other deformable porous media in a manner that can be readily implemented within the framework of Digital Rock Physics.more » « less
-
Reactive transport modeling of subsurface environments plays an important role in addressing critical problems of geochemical processes, such as dissolution and precipitation of minerals. Current transport models for porous media span various scales, ranging from pore-scale to continuum-scale. In this study, we established an upscaling method connecting pore-scale and continuum-scale models by employing a deep learning methodology of Convolutional Neural Networks (CNNs). We applied Darcy-Brinkmann-Stokes (DBS) method to simulate the fluid flow and reactive transport in pore-scale models, which would act as constituents of a continuum-scale model. The datasets of spatial pore distribution of subvolume samples were used as the input for the deep learning model, and the continuum (Darcy)-scale parameters such as permeability, effective surface area, and effective diffusion coefficient were figured out as outputs (i.e., labels). By applying the trained models of the subvolumes in the entire sample volume, we generated the initial field of porosity, permeability, effective diffusion coefficient, and effective surface area for continuum-scale simulation of a mineral dissolution problem. We took an acid dissolution case as an example to utilize the outcomes of trained deep learning models as input data in the continuum-scale simulation. This work presents a comprehensive upscaling workflow, as bridging the findings of microscale simulations to the continuum-scale simulations of a reactive transport problem.more » « less
An official website of the United States government

