Developing an Effective Model for Shale Gas Flow in Nano-scale Pore Clusters based on FIB-SEM Images

NASA Astrophysics Data System (ADS)

Jiang, W. B.; Lin, M.; Yi, Z. X.; Li, H. S.

2016-12-01

Nano-scale pores existed in the form of clusters are the controlling void space in shale gas reservoir. Gas transport in nanopores which has a significant influence on shale gas' recoverability displays multiple transport regimes, including viscous, slippage flow and Knudsen diffusion. In addition, it is also influenced by pore space characteristics. For convenience and efficiency consideration, it is necessary to develop an upscaling model from nano pore to pore cluster scale. Existing models are more like framework functions that provide a format, because the parameters that represent pore space characteristics are underdetermined and may have multiple possibilities. Therefore, it is urgent to make them clear and obtained a model that is closer to reality. FIB-SEM imaging technology is able to acquire three dimensional images with nanometer resolution that nano pores can be visible. Based on the images of two shale samples, we used a high-precision pore network extraction algorithm to generate equivalent pore networks and simulate multiple regime (non-Darcy) flow in it. Several structural parameters can be obtained through pore network modelling. It is found that although the throat-radius distributions are very close, throat flux-radius distributions of different samples can be divided into two categories. The variation of tortuosity with pressure and the overall trend of throat-flux distribution changes with pressure are disclosed. A deeper understanding of shale gas flow in nano-scale pore clusters is obtained. After all, an upscaling model that connects absolute permeability, apparent permeability and other characteristic parameters is proposed, and the best parameter scheme considering throat number-radius distribution and flowing porosity for this model is selected out of three schemes based on pore scale results, and it can avoid multiple-solution problem and is useful in reservoir modelling and experiment result analysis, etc. This work is supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB10020302), the National Natural Science Foundation of China (41574129), the Key Instrument Developing Project of the CAS (ZDYZ2012-1-08-02), the 973 Program (2014CB239004).

Investigation of Coupled model of Pore network and Continuum in shale gas

NASA Astrophysics Data System (ADS)

Cao, G.; Lin, M.

2016-12-01

Flow in shale spanning over many scales, makes the majority of conventional treatment methods disabled. For effectively simulating, a coupled model of pore-scale and continuum-scale was proposed in this paper. Based on the SEM image, we decompose organic-rich-shale into two subdomains: kerogen and inorganic matrix. In kerogen, the nanoscale pore-network is the main storage space and migration pathway so that the molecular phenomena (slip and diffusive transport) is significant. Whereas, inorganic matrix, with relatively large pores and micro fractures, the flow is approximate to Darcy. We use pore-scale network models (PNM) to represent kerogen and continuum-scale models (FVM or FEM) to represent matrix. Finite element mortars are employed to couple pore- and continuum-scale models by enforcing continuity of pressures and fluxes at shared boundary interfaces. In our method, the process in the coupled model is described by pressure square equation, and uses Dirichlet boundary conditions. We discuss several problems: the optimal element number of mortar faces, two categories boundary faces of pore network, the difference between 2D and 3D models, and the difference between continuum models FVM and FEM in mortars. We conclude that: (1) too coarse mesh in mortars will decrease the accuracy, while too fine mesh will lead to an ill-condition even singular system, the optimal element number is depended on boundary pores and nodes number. (2) pore network models are adjacent to two different mortar faces (PNM to PNM, PNM to continuum model), incidental repeated mortar nodes must be deleted. (3) 3D models can be replaced by 2D models under certain condition. (4) FVM is more convenient than FEM, for its simplicity in assigning interface nodes pressure and calculating interface fluxes. This work is supported by the Strategic Priority Research Program of the Chinese Academy of Sciences (XDB10020302), the 973 Program (2014CB239004), the Key Instrument Developing Project of the CAS (ZDYZ2012-1-08-02), the National Natural Science Foundation of China (41574129).

Reply to 'Comments on upscaling geochemical reaction rates usingpore-scale network modeling' by Peter C. Lichtner and Qinjun Kang

DOE Office of Scientific and Technical Information (OSTI.GOV)

Li, Li; Peters, Catherine A.; Celia, Michael A.

2006-05-03

Our paper "Upscaling geochemical reaction rates usingpore-scale network modeling" presents a novel application of pore-scalenetwork modeling to upscale mineral dissolution and precipitationreaction rates from the pore scale to the continuum scale, anddemonstrates the methodology by analyzing the scaling behavior ofanorthite and kaolinite reaction kinetics under conditions related to CO2sequestration. We conclude that under highly acidic conditions relevantto CO2 sequestration, the traditional continuum-based methodology may notcapture the spatial variation in concentrations from pore to pore, andscaling tools may be important in correctly modeling reactive transportprocesses in such systems. This work addresses the important butdifficult question of scaling mineral dissolution and precipitationreactionmore » kinetics, which is often ignored in fields such as geochemistry,water resources, and contaminant hydrology. Although scaling of physicalprocesses has been studied for almost three decades, very few studieshave examined the scaling issues related to chemical processes, despitetheir importance in governing the transport and fate of contaminants insubsurface systems.« less

An Image-based Micro-continuum Pore-scale Model for Gas Transport in Organic-rich Shale

NASA Astrophysics Data System (ADS)

Guo, B.; Tchelepi, H.

2017-12-01

Gas production from unconventional source rocks, such as ultra-tight shales, has increased significantly over the past decade. However, due to the extremely small pores ( 1-100 nm) and the strong material heterogeneity, gas flow in shale is still not well understood and poses challenges for predictive field-scale simulations. In recent years, digital rock analysis has been applied to understand shale gas transport at the pore-scale. An issue with rock images (e.g. FIB-SEM, nano-/micro-CT images) is the so-called "cutoff length", i.e., pores and heterogeneities below the resolution cannot be resolved, which leads to two length scales (resolved features and unresolved sub-resolution features) that are challenging for flow simulations. Here we develop a micro-continuum model, modified from the classic Darcy-Brinkman-Stokes framework, that can naturally couple the resolved pores and the unresolved nano-porous regions. In the resolved pores, gas flow is modeled with Stokes equation. In the unresolved regions where the pore sizes are below the image resolution, we develop an apparent permeability model considering non-Darcy flow at the nanoscale including slip flow, Knudsen diffusion, adsorption/desorption, surface diffusion, and real gas effect. The end result is a micro-continuum pore-scale model that can simulate gas transport in 3D reconstructed shale images. The model has been implemented in the open-source simulation platform OpenFOAM. In this paper, we present case studies to demonstrate the applicability of the model, where we use 3D segmented FIB-SEM and nano-CT shale images that include four material constituents: organic matter, clay, granular mineral, and pore. In addition to the pore structure and the distribution of the material constituents, we populate the model with experimental measurements (e.g. size distribution of the sub-resolution pores from nitrogen adsorption) and parameters from the literature and identify the relative importance of different physics on gas production. Overall, the micro-continuum model provides a novel tool for digital rock analysis of organic-rich shale.

Pore and Continuum Scale Study of the Effect of Subgrid Transport Heterogeneity on Redox Reaction Rates

DOE Office of Scientific and Technical Information (OSTI.GOV)

Liu, Yuanyuan; Liu, Chongxuan; Zhang, Changyong

2015-08-01

A micromodel system with a pore structure for heterogeneous flow and transport was used to investigate the effect of subgrid transport heterogeneity on redox reaction rates. Hematite reductive dissolution by injecting a reduced form of flavin mononucleotide (FMNH2) at variable flow rates was used as an example to probe the variations of redox reaction rates in different subgrid transport domains. Experiments, pore-scale simulations, and macroscopic modeling were performed to measure and simulate in-situ hematite reduction and to evaluate the scaling behavior of the redox reaction rates from the pore to macroscopic scales. The results indicated that the measured pore-scale ratesmore » of hematite reduction were consistent with the predictions from a pore scale reactive transport model. A general trend is that hematite reduction followed reductant transport pathways, starting from the advection-dominated pores toward the interior of diffusion-dominated domains. Two types of diffusion domains were considered in the micromodel: a micropore diffusion domain, which locates inside solid grains or aggregates where reactant transport is limited by diffusion; and a macropore diffusion domain, which locates at wedged, dead-end pore spaces created by the grain-grain contacts. The rate of hematite reduction in the advection-dominated domain was faster than those in the diffusion-controlled domains, and the rate in the macropore diffusion domain was faster than that in the micropore domain. The reduction rates in the advection and macropore diffusion domains increased with increasing flow rate, but were affected by different mechanisms. The rate increase in the advection domain was controlled by the mass action effect as a faster flow supplied more reactants, and the rate increase in the macropore domain was more affected by the rate of mass exchange with the advection domain, which increased with increasing flow rate. The hematite reduction rate in the micropore domain was, however, not affected by the flow rate because molecular diffusion limits reductant supply to the micropore domain interior. Domain-based macroscopic models were evaluated to scale redox reaction rates from the pore to macroscopic scales. A single domain model, which ignores subgrid transport heterogeneity deviated significantly from the pore-scale results. Further analysis revealed that the rate expression for hematite reduction was not scalable from the pore to porous media using the single domain model. A three-domain model, which effectively considers subgrid reactive diffusion in the micropore and macropore domains, significantly improved model description. Overall this study revealed the importance of subgrid transport heterogeneity in the manifestation of redox reaction rates in porous media and in scaling reactions from the pore to porous media. The research also supported that the domain-based scaling approach can be used to directly scale redox reactions in porous media with subgrid transport heterogeneity.« less

Dynamic Pore-Scale Imaging of Reactive Transport in Heterogeneous Carbonates at Reservior Conditions

NASA Astrophysics Data System (ADS)

Menke, Hannah; Bijeljic, Branko; Andrew, Matthew; Blunt, Martin

2014-05-01

Sequestering carbon in deep geologic formations is one way of reducing anthropogenic CO2 emissions. Carbon capture, Utilization, and Storage (CCUS) in carbonate reservoirs has the added benefit of mobilizing more oil for extraction, increasing oil reservoir yield, and generating revenue while also mitigating climate change. The magnitude, speed, and type of dissolution are dependent the intrinsic properties of the rock. Understanding how small changes in the pore structure affect dissolution is paramount for successful predictive modelling both on the pore-scale and for up-scaled reservoir simulations. We propose an experimental method whereby both 'Pink Beam' synchrotron radiation and a Micro-CT lab source are used in dynamic X-ray microtomography to investigate the pore structure changes in carbonate rocks of varying heterogeneity at high temperatures and pressures. Four carbonate rock types were studied, two relatively homogeneous carbonates, Ketton and Mt. Gambier, and two very heterogeneous carbonates, Estalliades and Portland Basebed. Each rock type was imaged under the same reservoir and flow conditions to gain insight into the impact of heterogeneity. A 4-mm carbonate core was injected with CO2-saturated brine at 10 MPa and 50oC for 2 hours. Depending on sample heterogeneity and X-ray source, tomographic images were taken at between 30-second and 20-minute time-resolutions and a 4-micron spatial resolution during injection. Changes in porosity, permeability, and structure were obtained by first binning and filtering the images, then binarizing them with watershed segmentation, and finally extracting a pore/throat network. Furthermore, pore-scale flow modelling was performed directly on the binarized image and used to track velocity distributions as the pore network evolved. Significant differences in dissolution type and magnitude were found for each rock type. The most homogeneous carbonate, Ketton, was seen to have predominately uniform dissolution with minor dissolution rate differences between the pores and pore throats. This was not true for the heterogeneous carbonates, Estalliades and Portland Basebed, which formed wormholes. Pore-scale modelling of flow directly on the voxels showed the differences in the evolution of complex flow fields with changes in dissolution regime. The PDFs of normalized velocity for uniform dissolution showed that the maximum pore velocity within the system decreased as dissolution occurred. This is due to dissolution enlarging pores and pore throats. However, in the wormholing regime, there was a large increase in maximum velocity once the wormhole broke through the length of the core and a preferential flow path was created. Additionally, this study serves as a unique benchmark for pore-scale reactive transport modelling directly on the binarized Micro-CT images. This dynamic pore-scale imaging method offers advantages in helping fully explain the dominant physical and chemical processes at the pore scale so that they may be up-scaled to the reservoir scale for increased accuracy in model prediction.

Multiscale Simulation of Porous Ceramics Based on Movable Cellular Automaton Method

NASA Astrophysics Data System (ADS)

Smolin, A.; Smolin, I.; Eremina, G.; Smolina, I.

2017-10-01

The paper presents a model for simulating mechanical behaviour of multiscale porous ceramics based on movable cellular automaton method, which is a novel particle method in computational mechanics of solid. The initial scale of the proposed approach corresponds to the characteristic size of the smallest pores in the ceramics. At this scale, we model uniaxial compression of several representative samples with an explicit account of pores of the same size but with the random unique position in space. As a result, we get the average values of Young’s modulus and strength, as well as the parameters of the Weibull distribution of these properties at the current scale level. These data allow us to describe the material behaviour at the next scale level were only the larger pores are considered explicitly, while the influence of small pores is included via the effective properties determined at the previous scale level. If the pore size distribution function of the material has N maxima we need to perform computations for N - 1 levels in order to get the properties from the lowest scale up to the macroscale step by step. The proposed approach was applied to modelling zirconia ceramics with bimodal pore size distribution. The obtained results show correct behaviour of the model sample at the macroscale.

Multi-Scale Multi-Physics Modeling of Matrix Transport Properties in Fractured Shale Reservoirs

NASA Astrophysics Data System (ADS)

Mehmani, A.; Prodanovic, M.

2014-12-01

Understanding the shale matrix flow behavior is imperative in successful reservoir development for hydrocarbon production and carbon storage. Without a predictive model, significant uncertainties in flowback from the formation, the communication between the fracture and matrix as well as proper fracturing practice will ensue. Informed by SEM images, we develop deterministic network models that couple pores from multiple scales and their respective fluid physics. The models are used to investigate sorption hysteresis as an affordable way of inferring the nanoscale pore structure in core scale. In addition, restricted diffusion as a function of pore shape, pore-throat size ratios and network connectivity is computed to make correct interpretation of the 2D NMR maps possible. Our novel pore network models have the ability to match sorption hysteresis measurements without any tuning parameters. The results clarify a common misconception of linking type 3 nitrogen hysteresis curves to only the shale pore shape and show promising sensitivty for nanopore structre inference in core scale. The results on restricted diffusion shed light on the importance of including shape factors in 2D NMR interpretations. A priori "weighting factors" as a function of pore-throat and throat-length ratio are presented and the effect of network connectivity on diffusion is quantitatively assessed. We are currently working on verifying our models with experimental data gathered from the Eagleford formation.

Dynamic Pore-Scale Imaging of Reactive Transport in Heterogeneous Carbonates at Reservoir Conditions Across Multiple Dissolution Regimes

NASA Astrophysics Data System (ADS)

Menke, H. P.; Bijeljic, B.; Andrew, M. G.; Blunt, M. J.

2014-12-01

Sequestering carbon in deep geologic formations is one way of reducing anthropogenic CO2 emissions. When supercritical CO2 mixes with brine in a reservoir, the acid generated has the potential to dissolve the surrounding pore structure. However, the magnitude and type of dissolution are condition dependent. Understanding how small changes in the pore structure, chemistry, and flow properties affect dissolution is paramount for successful predictive modelling. Both 'Pink Beam' synchrotron radiation and a Micro-CT lab source are used in dynamic X-ray microtomography to investigate the pore structure changes during supercritical CO2 injection in carbonate rocks of varying heterogeneity at high temperatures and pressures and various flow-rates. Three carbonate rock types were studied, one with a homogeneous pore structure and two heterogeneous carbonates. All samples are practically pure calcium carbonate, but have widely varying rock structures. Flow-rate was varied in three successive experiments by over an order of magnitude whlie keeping all other experimental conditions constant. A 4-mm carbonate core was injected with CO2-saturated brine at 10 MPa and 50oC. Tomographic images were taken at 30-second to 20-minute time-resolutions during a 2 to 4-hour injection period. A pore network was extracted using a topological analysis of the pore space and pore-scale flow modelling was performed directly on the binarized images with connected pathways and used to track the altering velocity distributions. Significant differences in dissolution type and magnitude were found for each rock type and flowrate. At the highest flow-rates, the homogeneous carbonate was seen to have predominately uniform dissolution with minor dissolution rate differences between the pores and pore throats. Alternatively, the heterogeneous carbonates which formed wormholes at high flow rates. At low flow rates the homogeneous rock developed wormholes, while the heterogeneous samples showed evidence of compact dissolution. This study serves as a unique benchmark for pore-scale reactive transport modelling directly on the binarized Micro-CT images. Dynamic pore-scale imaging methods offer advantages in helping explain the dominant processes at the pore scale so that they may be up-scaled for accurate model prediction.

Pore-scale simulation of calcium carbonate precipitation and dissolution under highly supersaturated conditions in a microfludic pore network

NASA Astrophysics Data System (ADS)

Yoon, H.; Dewers, T. A.; Valocchi, A. J.; Werth, C. J.

2011-12-01

Dissolved CO2 during geological CO2 storage may react with minerals in fractured rocks or confined aquifers and cause mineral precipitation. The overall rate of reaction can be affected by coupled processes among hydrodynamics, transport, and reactions at pore-scale. Pore-scale models of coupled fluid flow, reactive transport, and CaCO3 precipitation and dissolution are applied to account for transient experimental results of CaCO3 precipitation and dissolution under highly supersaturated conditions in a microfluidic pore network (i.e., micromodel). Pore-scale experiments in the micromodel are used as a basis for understanding coupled physics of systems perturbed by geological CO2 injection. In the micromodel, precipitation is induced by transverse mixing along the centerline in pore bodies. Overall, the pore-scale model qualitatively captured the governing physics of reactions such as precipitate morphology, precipitation rate, and maximum precipitation area in first few pore spaces. In particular, we found that proper estimation of the effective diffusion coefficient and the reactive surface area is necessary to adequately simulate precipitation and dissolution rates. As the model domain increases, the effect of flow patterns affected by precipitation on the overall reaction rate also increases. The model is also applied to account for the effect of different reaction rate laws on mineral precipitation and dissolution at pore-scale. Reaction rate laws tested include the linear rate law, nonlinear power law, and newly-developed rate law based on in-situ measurements at nano scale in the literature. Progress on novel methods for upscaling pore-scale models for reactive transport are discussed, and are being applied to mineral precipitation patterns observed in natural analogues. H.Y. and T. D. were supported as part of the Center for Frontiers of Subsurface Energy Security, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Award Number DE-SC0001114. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Multi-scale modeling of multi-component reactive transport in geothermal aquifers

NASA Astrophysics Data System (ADS)

Nick, Hamidreza M.; Raoof, Amir; Wolf, Karl-Heinz; Bruhn, David

2014-05-01

In deep geothermal systems heat and chemical stresses can cause physical alterations, which may have a significant effect on flow and reaction rates. As a consequence it will lead to changes in permeability and porosity of the formations due to mineral precipitation and dissolution. Large-scale modeling of reactive transport in such systems is still challenging. A large area of uncertainty is the way in which the pore-scale information controlling the flow and reaction will behave at a larger scale. A possible choice is to use constitutive relationships relating, for example the permeability and porosity evolutions to the change in the pore geometry. While determining such relationships through laboratory experiments may be limited, pore-network modeling provides an alternative solution. In this work, we introduce a new workflow in which a hybrid Finite-Element Finite-Volume method [1,2] and a pore network modeling approach [3] are employed. Using the pore-scale model, relevant constitutive relations are developed. These relations are then embedded in the continuum-scale model. This approach enables us to study non-isothermal reactive transport in porous media while accounting for micro-scale features under realistic conditions. The performance and applicability of the proposed model is explored for different flow and reaction regimes. References: 1. Matthäi, S.K., et al.: Simulation of solute transport through fractured rock: a higher-order accurate finite-element finite-volume method permitting large time steps. Transport in porous media 83.2 (2010): 289-318. 2. Nick, H.M., et al.: Reactive dispersive contaminant transport in coastal aquifers: Numerical simulation of a reactive Henry problem. Journal of contaminant hydrology 145 (2012), 90-104. 3. Raoof A., et al.: PoreFlow: A Complex pore-network model for simulation of reactive transport in variably saturated porous media, Computers & Geosciences, 61, (2013), 160-174.

Water permeability in hydrate-bearing sediments: A pore-scale study

NASA Astrophysics Data System (ADS)

Dai, Sheng; Seol, Yongkoo

2014-06-01

Permeability is a critical parameter governing methane flux and fluid flow in hydrate-bearing sediments; however, limited valid data are available due to experimental challenges. Here we investigate the relationship between apparent water permeability (k') and hydrate saturation (Sh), accounting for hydrate pore-scale growth habit and meso-scale heterogeneity. Results from capillary tube models rely on cross-sectional tube shapes and hydrate pore habits, thus are appropriate only for sediments with uniform hydrate distribution and known hydrate pore character. Given our pore network modeling results showing that accumulating hydrate in sediments decreases sediment porosity and increases hydraulic tortuosity, we propose a modified Kozeny-Carman model to characterize water permeability in hydrate-bearing sediments. This model agrees well with experimental results and can be easily implemented in reservoir simulators with no empirical variables other than Sh. Results are also relevant to flow through other natural sediments that undergo diagenesis, salt precipitation, or bio-clogging.

Investigating the relative permeability behavior of microporosity-rich carbonates and tight sandstones with multiscale pore network models

NASA Astrophysics Data System (ADS)

Bultreys, Tom; Stappen, Jeroen Van; Kock, Tim De; Boever, Wesley De; Boone, Marijn A.; Hoorebeke, Luc Van; Cnudde, Veerle

2016-11-01

The relative permeability behavior of rocks with wide ranges of pore sizes is in many cases still poorly understood and is difficult to model at the pore scale. In this work, we investigate the capillary pressure and relative permeability behavior of three outcrop carbonates and two tight reservoir sandstones with wide, multimodal pore size distributions. To examine how the drainage and imbibition properties of these complex rock types are influenced by the connectivity of macropores to each other and to zones with unresolved small-scale porosity, we apply a previously presented microcomputed-tomography-based multiscale pore network model to these samples. The sensitivity to the properties of the small-scale porosity is studied by performing simulations with different artificial sphere-packing-based networks as a proxy for these pores. Finally, the mixed-wet water-flooding behavior of the samples is investigated, assuming different wettability distributions for the microporosity and macroporosity. While this work is not an attempt to perform predictive modeling, it seeks to qualitatively explain the behavior of the investigated samples and illustrates some of the most recent developments in multiscale pore network modeling.

Pore-Scale Determination of Gas Relative Permeability in Hydrate-Bearing Sediments Using X-Ray Computed Micro-Tomography and Lattice Boltzmann Method

NASA Astrophysics Data System (ADS)

Chen, Xiongyu; Verma, Rahul; Espinoza, D. Nicolas; Prodanović, Maša.

2018-01-01

This work uses X-ray computed micro-tomography (μCT) to monitor xenon hydrate growth in a sandpack under the excess gas condition. The μCT images give pore-scale hydrate distribution and pore habit in space and time. We use the lattice Boltzmann method to calculate gas relative permeability (krg) as a function of hydrate saturation (Shyd) in the pore structure of the experimental hydrate-bearing sand retrieved from μCT data. The results suggest the krg - Shyd data fit well a new model krg = (1-Shyd)·exp(-4.95·Shyd) rather than the simple Corey model. In addition, we calculate krg-Shyd curves using digital models of hydrate-bearing sand based on idealized grain-attaching, coarse pore-filling, and dispersed pore-filling hydrate habits. Our pore-scale measurements and modeling show that the krg-Shyd curves are similar regardless of whether hydrate crystals develop grain-attaching or coarse pore-filling habits. The dispersed pore filling habit exhibits much lower gas relative permeability than the other two, but it is not observed in the experiment and not compatible with Ostwald ripening mechanisms. We find that a single grain-shape factor can be used in the Carman-Kozeny equation to calculate krg-Shyd data with known porosity and average grain diameter, suggesting it is a useful model for hydrate-bearing sand.

Pore-Scale Simulation and Sensitivity Analysis of Apparent Gas Permeability in Shale Matrix

PubMed Central

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N.

2017-01-01

Extremely low permeability due to nano-scale pores is a distinctive feature of gas transport in a shale matrix. The permeability of shale depends on pore pressure, porosity, pore throat size and gas type. The pore network model is a practical way to explain the macro flow behavior of porous media from a microscopic point of view. In this research, gas flow in a shale matrix is simulated using a previously developed three-dimensional pore network model that includes typical bimodal pore size distribution, anisotropy and low connectivity of the pore structure in shale. The apparent gas permeability of shale matrix was calculated under different reservoir pressures corresponding to different gas exploitation stages. Results indicate that gas permeability is strongly related to reservoir gas pressure, and hence the apparent permeability is not a unique value during the shale gas exploitation, and simulations suggested that a constant permeability for continuum-scale simulation is not accurate. Hence, the reservoir pressures of different shale gas exploitations should be considered. In addition, a sensitivity analysis was also performed to determine the contributions to apparent permeability of a shale matrix from petro-physical properties of shale such as pore throat size and porosity. Finally, the impact of connectivity of nano-scale pores on shale gas flux was analyzed. These results would provide an insight into understanding nano/micro scale flows of shale gas in the shale matrix. PMID:28772465

Pore-Scale Simulation and Sensitivity Analysis of Apparent Gas Permeability in Shale Matrix.

PubMed

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N

2017-01-25

Extremely low permeability due to nano-scale pores is a distinctive feature of gas transport in a shale matrix. The permeability of shale depends on pore pressure, porosity, pore throat size and gas type. The pore network model is a practical way to explain the macro flow behavior of porous media from a microscopic point of view. In this research, gas flow in a shale matrix is simulated using a previously developed three-dimensional pore network model that includes typical bimodal pore size distribution, anisotropy and low connectivity of the pore structure in shale. The apparent gas permeability of shale matrix was calculated under different reservoir pressures corresponding to different gas exploitation stages. Results indicate that gas permeability is strongly related to reservoir gas pressure, and hence the apparent permeability is not a unique value during the shale gas exploitation, and simulations suggested that a constant permeability for continuum-scale simulation is not accurate. Hence, the reservoir pressures of different shale gas exploitations should be considered. In addition, a sensitivity analysis was also performed to determine the contributions to apparent permeability of a shale matrix from petro-physical properties of shale such as pore throat size and porosity. Finally, the impact of connectivity of nano-scale pores on shale gas flux was analyzed. These results would provide an insight into understanding nano/micro scale flows of shale gas in the shale matrix.

Upscaling of reaction rates in reactive transport using pore-scale reactive transport model

NASA Astrophysics Data System (ADS)

Yoon, H.; Dewers, T. A.; Arnold, B. W.; Major, J. R.; Eichhubl, P.; Srinivasan, S.

2013-12-01

Dissolved CO2 during geological CO2 storage may react with minerals in fractured rocks, confined aquifers, or faults, resulting in mineral precipitation and dissolution. The overall rate of reaction can be affected by coupled processes among hydrodynamics, transport, and reactions at the (sub) pore-scale. In this research pore-scale modeling of coupled fluid flow, reactive transport, and heterogeneous reaction at the mineral surface is applied to account for permeability alterations caused by precipitation-induced pore-blocking. This work is motivated by the observed CO2 seeps from a natural analog to geologic CO2 sequestration at Crystal Geyser, Utah. A key observation is the lateral migration of CO2 seep sites at a scale of ~ 100 meters over time. A pore-scale model provides fundamental mechanistic explanations of how calcite precipitation alters flow paths by pore plugging under different geochemical compositions and pore configurations. In addition, response function of reaction rates will be constructed from pore-scale simulations which account for a range of reaction regimes characterized by the Damkohler and Peclet numbers. Newly developed response functions will be used in a continuum scale model that may account for large-scale phenomena mimicking lateral migration of surface CO2 seeps. Comparison of field observations and simulations results will provide mechanistic explanations of the lateral migration and enhance our understanding of subsurface processes associated with the CO2 injection. This work is supported as part of the Center for Frontiers of Subsurface Energy Security, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Award Number DE-SC0001114. Sandia National Laboratories is a multi-program laboratory managed and operated by Sandia Corporation, a wholly owned subsidiary of Lockheed Martin Corporation, for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-AC04-94AL85000.

Global sensitivity analysis of multiscale properties of porous materials

NASA Astrophysics Data System (ADS)

Um, Kimoon; Zhang, Xuan; Katsoulakis, Markos; Plechac, Petr; Tartakovsky, Daniel M.

2018-02-01

Ubiquitous uncertainty about pore geometry inevitably undermines the veracity of pore- and multi-scale simulations of transport phenomena in porous media. It raises two fundamental issues: sensitivity of effective material properties to pore-scale parameters and statistical parameterization of Darcy-scale models that accounts for pore-scale uncertainty. Homogenization-based maps of pore-scale parameters onto their Darcy-scale counterparts facilitate both sensitivity analysis (SA) and uncertainty quantification. We treat uncertain geometric characteristics of a hierarchical porous medium as random variables to conduct global SA and to derive probabilistic descriptors of effective diffusion coefficients and effective sorption rate. Our analysis is formulated in terms of solute transport diffusing through a fluid-filled pore space, while sorbing to the solid matrix. Yet it is sufficiently general to be applied to other multiscale porous media phenomena that are amenable to homogenization.

A New Dual-Pore Formation Factor Model: A Percolation Network Study and Comparison to Experimental Data

NASA Astrophysics Data System (ADS)

Tang, Y. B.; Li, M.; Bernabe, Y.

2014-12-01

We modeled the electrical transport behavior of dual-pore carbonate rocks in this paper. Based on experimental data of a carbonate reservoir in China, we simply considered the low porosity samples equivalent to the matrix (micro-pore system) of the high porosity samples. For modeling the bimodal porous media, we considered that the matrix is homogeneous and interconnected. The connectivity and the pore size distribution of macro-pore system are varied randomly. Both pore systems are supposed to act electrically in parallel, connected at the nodes, where the fluid exchange takes place, an approach previously used by Bauer et al. (2012). Then, the effect of the properties of matrix, the pore size distribution and connectivity of macro-pore system on petrophysical properties of carbonates can be investigated. We simulated electrical current through networks in three-dimensional simple cubic (SC) and body-center cubic (BCC) with different coordination numbers and different pipe radius distributions of macro-pore system. Based on the simulation results, we found that the formation factor obeys a "universal" scaling relationship (i.e. independent of lattice type), 1/F∝eγz, where γ is a function of the normalized standard deviation of the pore radius distribution of macro-pore system and z is the coordination number of macro-pore system. This relationship is different from the classic "universal power law" in percolation theory. A formation factor model was inferred on the basis of the scaling relationship mentioned above and several scale-invariant quantities (such as hydraulic radius rH and throat length l of macro-pore). Several methods were developed to estimate corresponding parameters of the new model with conventional core analyses. It was satisfactorily tested against experimental data, including some published experimental data. Furthermore, the relationship between water saturation and resistivity in dual-pore carbonates was discussed based on the new model.

An investigation into preserving spatially-distinct pore systems in multi-component rocks using a fossiliferous limestone example

NASA Astrophysics Data System (ADS)

Jiang, Zeyun; Couples, Gary D.; Lewis, Helen; Mangione, Alessandro

2018-07-01

Limestones containing abundant disc-shaped fossil Nummulites can form significant hydrocarbon reservoirs but they have a distinctly heterogeneous distribution of pore shapes, sizes and connectivities, which make it particularly difficult to calculate petrophysical properties and consequent flow outcomes. The severity of the problem rests on the wide length-scale range from the millimetre scale of the fossil's pore space to the micron scale of rock matrix pores. This work develops a technique to incorporate multi-scale void systems into a pore network, which is used to calculate the petrophysical properties for subsequent flow simulations at different stages in the limestone's petrophysical evolution. While rock pore size, shape and connectivity can be determined, with varying levels of fidelity, using techniques such as X-ray computed tomography (CT) or scanning electron microscopy (SEM), this work represents a more challenging class where the rock of interest is insufficiently sampled or, as here, has been overprinted by extensive chemical diagenesis. The main challenge is integrating multi-scale void structures derived from both SEM and CT images, into a single model or a pore-scale network while still honouring the nature of the connections across these length scales. Pore network flow simulations are used to illustrate the technique but of equal importance, to demonstrate how supportable earlier-stage petrophysical property distributions can be used to assess the viability of several potential geological event sequences. The results of our flow simulations on generated models highlight the requirement for correct determination of the dominant pore scales (one plus of nm, μm, mm, cm), the spatial correlation and the cross-scale connections.

Micro-computed tomography pore-scale study of flow in porous media: Effect of voxel resolution

NASA Astrophysics Data System (ADS)

Shah, S. M.; Gray, F.; Crawshaw, J. P.; Boek, E. S.

2016-09-01

A fundamental understanding of flow in porous media at the pore-scale is necessary to be able to upscale average displacement processes from core to reservoir scale. The study of fluid flow in porous media at the pore-scale consists of two key procedures: Imaging - reconstruction of three-dimensional (3D) pore space images; and modelling such as with single and two-phase flow simulations with Lattice-Boltzmann (LB) or Pore-Network (PN) Modelling. Here we analyse pore-scale results to predict petrophysical properties such as porosity, single-phase permeability and multi-phase properties at different length scales. The fundamental issue is to understand the image resolution dependency of transport properties, in order to up-scale the flow physics from pore to core scale. In this work, we use a high resolution micro-computed tomography (micro-CT) scanner to image and reconstruct three dimensional pore-scale images of five sandstones (Bentheimer, Berea, Clashach, Doddington and Stainton) and five complex carbonates (Ketton, Estaillades, Middle Eastern sample 3, Middle Eastern sample 5 and Indiana Limestone 1) at four different voxel resolutions (4.4 μm, 6.2 μm, 8.3 μm and 10.2 μm), scanning the same physical field of view. Implementing three phase segmentation (macro-pore phase, intermediate phase and grain phase) on pore-scale images helps to understand the importance of connected macro-porosity in the fluid flow for the samples studied. We then compute the petrophysical properties for all the samples using PN and LB simulations in order to study the influence of voxel resolution on petrophysical properties. We then introduce a numerical coarsening scheme which is used to coarsen a high voxel resolution image (4.4 μm) to lower resolutions (6.2 μm, 8.3 μm and 10.2 μm) and study the impact of coarsening data on macroscopic and multi-phase properties. Numerical coarsening of high resolution data is found to be superior to using a lower resolution scan because it avoids the problem of partial volume effects and reduces the scaling effect by preserving the pore-space properties influencing the transport properties. This is evidently compared in this study by predicting several pore network properties such as number of pores and throats, average pore and throat radius and coordination number for both scan based analysis and numerical coarsened data.

Microfluidic Experiments Studying Pore Scale Interactions of Microbes and Geochemistry

NASA Astrophysics Data System (ADS)

Chen, M.; Kocar, B. D.

2016-12-01

Understanding how physical phenomena, chemical reactions, and microbial behavior interact at the pore-scale is crucial to understanding larger scale trends in groundwater chemistry. Recent studies illustrate the utility of microfluidic devices for illuminating pore-scale physical-biogeochemical processes and their control(s) on the cycling of iron, uranium, and other important elements 1-3. These experimental systems are ideal for examining geochemical reactions mediated by microbes, which include processes governed by complex biological phenomenon (e.g. biofilm formation, etc.)4. We present results of microfluidic experiments using a model metal reducing bacteria and varying pore geometries, exploring the limitations of the microorganisms' ability to access tight pore spaces, and examining coupled biogeochemical-physical controls on the cycling of redox sensitive metals. Experimental results will provide an enhanced understanding of coupled physical-biogeochemical processes transpiring at the pore-scale, and will constrain and compliment continuum models used to predict and describe the subsurface cycling of redox-sensitive elements5. 1. Vrionis, H. A. et al. Microbiological and geochemical heterogeneity in an in situ uranium bioremediation field site. Appl. Environ. Microbiol. 71, 6308-6318 (2005). 2. Pearce, C. I. et al. Pore-scale characterization of biogeochemical controls on iron and uranium speciation under flow conditions. Environ. Sci. Technol. 46, 7992-8000 (2012). 3. Zhang, C., Liu, C. & Shi, Z. Micromodel investigation of transport effect on the kinetics of reductive dissolution of hematite. Environ. Sci. Technol. 47, 4131-4139 (2013). 4. Ginn, T. R. et al. Processes in microbial transport in the natural subsurface. Adv. Water Resour. 25, 1017-1042 (2002). 5. Scheibe, T. D. et al. Coupling a genome-scale metabolic model with a reactive transport model to describe in situ uranium bioremediation. Microb. Biotechnol. 2, 274-286 (2009).

Hydraulic conductivity of variably saturated porous media: Film and corner flow in angular pore space

NASA Astrophysics Data System (ADS)

Tuller, Markus; Or, Dani

2001-05-01

Many models for hydraulic conductivity of partially saturated porous media rely on oversimplified representation of the pore space as a bundle of cylindrical capillaries and disregard flow in liquid films. Recent progress in modeling liquid behavior in angular pores of partially saturated porous media offers an alternative framework. We assume that equilibrium liquid-vapor interfaces provide well-defined and stable boundaries for slow laminar film and corner flow regimes in pore space comprised of angular pores connected to slit-shaped spaces. Knowledge of liquid configuration in the assumed geometry facilitates calculation of average liquid velocities in films and corners and enables derivation of pore-scale hydraulic conductivity as a function of matric potential. The pore-scale model is statistically upscaled to represent hydraulic conductivity for a sample of porous medium. Model parameters for the analytical sample-scale expressions are estimated from measured liquid retention data and other measurable medium properties. Model calculations illustrate the important role of film flow, whose contribution dominates capillary flow (in full pores and corners) at relatively high matric potentials (approximately -100 to -300 J kg-1, or -1 to 3 bars). The crossover region between film and capillary flow is marked by a significant change in the slope of the hydraulic conductivity function as often observed in measurements. Model predictions are compared with the widely applied van Genuchten-Mualem model and yield reasonable agreement with measured retention and hydraulic conductivity data over a wide range of soil textural classes.

Pore-scale spectral induced polarization (SIP) signaturesassociated with FeS biomineral transformations

DOE Office of Scientific and Technical Information (OSTI.GOV)

Slater, Lee; Ntarlagiannis, Dimitrios; Personna, Yves R.

2007-10-01

The authors measured Spectral Induced Polarization (SIP) signatures in sand columns during (1) FeS biomineralization produced by sulfate reducing bacteria (D. vulgaris) under anaerboci conditions, and (2) subsequent biomineral dissolution upon return to an aerobic state. The low-frequency (0.1-10 Hz peak) relaxations produced during biomineralization can be modeled with a Cole-Cole formulation, from which the evolution of the polarization magnitude and relaxation length scale can be estimated. They find that the modeled time constant is consistent with the polarizable elements being biomineral encrused pores. Evolution of the model parameters is consistent with FeS surface area increases and pore-size reduction duringmore » biomineral growth, and subsequent biomineral dissolution (FeS surface area decreases and pore expansion) upon return to the aerobic state. They conclude that SIP signatures are diagnostic of pore-scale geometrical changes associated with FeS biomineralization by sulfate reducing bacteria.« less

Evaluating the hydraulic and transport properties of peat soil using pore network modeling and X-ray micro computed tomography

NASA Astrophysics Data System (ADS)

Gharedaghloo, Behrad; Price, Jonathan S.; Rezanezhad, Fereidoun; Quinton, William L.

2018-06-01

Micro-scale properties of peat pore space and their influence on hydraulic and transport properties of peat soils have been given little attention so far. Characterizing the variation of these properties in a peat profile can increase our knowledge on the processes controlling contaminant transport through peatlands. As opposed to the common macro-scale (or bulk) representation of groundwater flow and transport processes, a pore network model (PNM) simulates flow and transport processes within individual pores. Here, a pore network modeling code capable of simulating advective and diffusive transport processes through a 3D unstructured pore network was developed; its predictive performance was evaluated by comparing its results to empirical values and to the results of computational fluid dynamics (CFD) simulations. This is the first time that peat pore networks have been extracted from X-ray micro-computed tomography (μCT) images of peat deposits and peat pore characteristics evaluated in a 3D approach. Water flow and solute transport were modeled in the unstructured pore networks mapped directly from μCT images. The modeling results were processed to determine the bulk properties of peat deposits. Results portray the commonly observed decrease in hydraulic conductivity with depth, which was attributed to the reduction of pore radius and increase in pore tortuosity. The increase in pore tortuosity with depth was associated with more decomposed peat soil and decreasing pore coordination number with depth, which extended the flow path of fluid particles. Results also revealed that hydraulic conductivity is isotropic locally, but becomes anisotropic after upscaling to core-scale; this suggests the anisotropy of peat hydraulic conductivity observed in core-scale and field-scale is due to the strong heterogeneity in the vertical dimension that is imposed by the layered structure of peat soils. Transport simulations revealed that for a given solute, the effective diffusion coefficient decreases with depth due to the corresponding increase of diffusional tortuosity. Longitudinal dispersivity of peat also was computed by analyzing advective-dominant transport simulations that showed peat dispersivity is similar to the empirical values reported in the same peat soil; it is not sensitive to soil depth and does not vary much along the soil profile.

Macroscopic modeling of heat and water vapor transfer with phase change in dry snow based on an upscaling method: Influence of air convection

NASA Astrophysics Data System (ADS)

Calonne, N.; Geindreau, C.; Flin, F.

2015-12-01

At the microscopic scale, i.e., pore scale, dry snow metamorphism is mainly driven by the heat and water vapor transfer and the sublimation-deposition process at the ice-air interface. Up to now, the description of these phenomena at the macroscopic scale, i.e., snow layer scale, in the snowpack models has been proposed in a phenomenological way. Here we used an upscaling method, namely, the homogenization of multiple-scale expansions, to derive theoretically the macroscopic equivalent modeling of heat and vapor transfer through a snow layer from the physics at the pore scale. The physical phenomena under consideration are steady state air flow, heat transfer by conduction and convection, water vapor transfer by diffusion and convection, and phase change (sublimation and deposition). We derived three different macroscopic models depending on the intensity of the air flow considered at the pore scale, i.e., on the order of magnitude of the pore Reynolds number and the Péclet numbers: (A) pure diffusion, (B) diffusion and moderate convection (Darcy's law), and (C) strong convection (nonlinear flow). The formulation of the models includes the exact expression of the macroscopic properties (effective thermal conductivity, effective vapor diffusion coefficient, and intrinsic permeability) and of the macroscopic source terms of heat and vapor arising from the phase change at the pore scale. Such definitions can be used to compute macroscopic snow properties from 3-D descriptions of snow microstructures. Finally, we illustrated the precision and the robustness of the proposed macroscopic models through 2-D numerical simulations.

Pore-Scale Modeling of Pore Structure Effects on P-Wave Scattering Attenuation in Dry Rocks

PubMed Central

Li, Tianyang; Qiu, Hao; Wang, Feifei

2015-01-01

Underground rocks usually have complex pore system with a variety of pore types and a wide range of pore size. The effects of pore structure on elastic wave attenuation cannot be neglected. We investigated the pore structure effects on P-wave scattering attenuation in dry rocks by pore-scale modeling based on the wave theory and the similarity principle. Our modeling results indicate that pore size, pore shape (such as aspect ratio), and pore density are important factors influencing P-wave scattering attenuation in porous rocks, and can explain the variation of scattering attenuation at the same porosity. From the perspective of scattering attenuation, porous rocks can safely suit to the long wavelength assumption when the ratio of wavelength to pore size is larger than 15. Under the long wavelength condition, the scattering attenuation coefficient increases as a power function as the pore density increases, and it increases exponentially with the increase in aspect ratio. For a certain porosity, rocks with smaller aspect ratio and/or larger pore size have stronger scattering attenuation. When the pore aspect ratio is larger than 0.5, the variation of scattering attenuation at the same porosity is dominantly caused by pore size and almost independent of the pore aspect ratio. These results lay a foundation for pore structure inversion from elastic wave responses in porous rocks. PMID:25961729

Application of a pore-scale reactive transport model to a natural analog for reaction-induced pore alterations

DOE PAGES

Yoon, Hongkyu; Major, Jonathan; Dewers, Thomas; ...

2017-01-05

Dissolved CO 2 in the subsurface resulting from geological CO 2 storage may react with minerals in fractured rocks, confined aquifers, or faults, resulting in mineral precipitation and dissolution. The overall rate of reaction can be affected by coupled processes including hydrodynamics, transport, and reactions at the (sub) pore-scale. In this work pore-scale modeling of coupled fluid flow, reactive transport, and heterogeneous reactions at the mineral surface is applied to account for permeability alterations caused by precipitation-induced pore-blocking. This paper is motivated by observations of CO 2 seeps from a natural CO 2 sequestration analog, Crystal Geyser, Utah. Observations alongmore » the surface exposure of the Little Grand Wash fault indicate the lateral migration of CO 2 seep sites (i.e., alteration zones) of 10–50 m width with spacing on the order of ~100 m over time. Sandstone permeability in alteration zones is reduced by 3–4 orders of magnitude by carbonate cementation compared to unaltered zones. One granular porous medium and one fracture network systems are used to conceptually represent permeable porous media and locations of conduits controlled by fault-segment intersections and/or topography, respectively. Simulation cases accounted for a range of reaction regimes characterized by the Damköhler (Da) and Peclet (Pe) numbers. Pore-scale simulation results demonstrate that combinations of transport (Pe), geochemical conditions (Da), solution chemistry, and pore and fracture configurations contributed to match key patterns observed in the field of how calcite precipitation alters flow paths by pore plugging. This comparison of simulation results with field observations reveals mechanistic explanations of the lateral migration and enhances our understanding of subsurface processes associated with the CO 2 injection. In addition, permeability and porosity relations are constructed from pore-scale simulations which account for a range of reaction regimes characterized by the Da and Pe numbers. Finally, the functional relationships obtained from pore-scale simulations can be used in a continuum scale model that may account for large-scale phenomena mimicking lateral migration of surface CO 2 seeps.« less

A multi-scale Lattice Boltzmann model for simulating solute transport in 3D X-ray micro-tomography images of aggregated porous materials

NASA Astrophysics Data System (ADS)

Zhang, Xiaoxian; Crawford, John W.; Flavel, Richard J.; Young, Iain M.

2016-10-01

The Lattice Boltzmann (LB) model and X-ray computed tomography (CT) have been increasingly used in combination over the past decade to simulate water flow and chemical transport at pore scale in porous materials. Because of its limitation in resolution and the hierarchical structure of most natural soils, the X-ray CT tomography can only identify pores that are greater than its resolution and treats other pores as solid. As a result, the so-called solid phase in X-ray images may in reality be a grey phase, containing substantial connected pores capable of conducing fluids and solute. Although modified LB models have been developed to simulate fluid flow in such media, models for solute transport are relatively limited. In this paper, we propose a LB model for simulating solute transport in binary soil images containing permeable solid phase. The model is based on the single-relaxation time approach and uses a modified partial bounce-back method to describe the resistance caused by the permeable solid phase to chemical transport. We derive the relationship between the diffusion coefficient and the parameter introduced in the partial bounce-back method, and test the model against analytical solution for movement of a pulse of tracer. We also validate it against classical finite volume method for solute diffusion in a simple 2D image, and then apply the model to a soil image acquired using X-ray tomography at resolution of 30 μm in attempts to analyse how the ability of the solid phase to diffuse solute at micron-scale affects the behaviour of the solute at macro-scale after a volumetric average. Based on the simulated results, we discuss briefly the danger in interpreting experimental results using the continuum model without fully understanding the pore-scale processes, as well as the potential of using pore-scale modelling and tomography to help improve the continuum models.

Adaptive Multiscale Modeling of Geochemical Impacts on Fracture Evolution

NASA Astrophysics Data System (ADS)

Molins, S.; Trebotich, D.; Steefel, C. I.; Deng, H.

2016-12-01

Understanding fracture evolution is essential for many subsurface energy applications, including subsurface storage, shale gas production, fracking, CO2 sequestration, and geothermal energy extraction. Geochemical processes in particular play a significant role in the evolution of fractures through dissolution-driven widening, fines migration, and/or fracture sealing due to precipitation. One obstacle to understanding and exploiting geochemical fracture evolution is that it is a multiscale process. However, current geochemical modeling of fractures cannot capture this multi-scale nature of geochemical and mechanical impacts on fracture evolution, and is limited to either a continuum or pore-scale representation. Conventional continuum-scale models treat fractures as preferential flow paths, with their permeability evolving as a function (often, a cubic law) of the fracture aperture. This approach has the limitation that it oversimplifies flow within the fracture in its omission of pore scale effects while also assuming well-mixed conditions. More recently, pore-scale models along with advanced characterization techniques have allowed for accurate simulations of flow and reactive transport within the pore space (Molins et al., 2014, 2015). However, these models, even with high performance computing, are currently limited in their ability to treat tractable domain sizes (Steefel et al., 2013). Thus, there is a critical need to develop an adaptive modeling capability that can account for separate properties and processes, emergent and otherwise, in the fracture and the rock matrix at different spatial scales. Here we present an adaptive modeling capability that treats geochemical impacts on fracture evolution within a single multiscale framework. Model development makes use of the high performance simulation capability, Chombo-Crunch, leveraged by high resolution characterization and experiments. The modeling framework is based on the adaptive capability in Chombo which not only enables mesh refinement, but also refinement of the model-pore scale or continuum Darcy scale-in a dynamic way such that the appropriate model is used only when and where it is needed. Explicit flux matching provides coupling betwen the scales.

Capillary filling rules and displacement mechanisms for spontaneous imbibition of CO2 for carbon storage and EOR using micro-model experiments and pore scale simulation

NASA Astrophysics Data System (ADS)

Chapman, E.; Yang, J.; Crawshaw, J.; Boek, E. S.

2012-04-01

In the 1980s, Lenormand et al. carried out their pioneering work on displacement mechanisms of fluids in etched networks [1]. Here we further examine displacement mechanisms in relation to capillary filling rules for spontaneous imbibition. Understanding the role of spontaneous imbibition in fluid displacement is essential for refining pore network models. Generally, pore network models use simple capillary filling rules and here we examine the validity of these rules for spontaneous imbibition. Improvement of pore network models is vital for the process of 'up-scaling' to the field scale for both enhanced oil recovery (EOR) and carbon sequestration. In this work, we present our experimental microfluidic research into the displacement of both supercritical CO2/deionised water (DI) systems and analogous n-decane/air - where supercritical CO2 and n-decane are the respective wetting fluids - controlled by imbibition at the pore scale. We conducted our experiments in etched PMMA and silicon/glass micro-fluidic hydrophobic chips. We first investigate displacement in single etched pore junctions, followed by displacement in complex network designs representing actual rock thin sections, i.e. Berea sandstone and Sucrosic dolomite. The n-decane/air experiments were conducted under ambient conditions, whereas the supercritical CO2/DI water experiments were conducted under high temperature and pressure in order to replicate reservoir conditions. Fluid displacement in all experiments was captured via a high speed video microscope. The direction and type of displacement the imbibing fluid takes when it enters a junction is dependent on the number of possible channels in which the wetting fluid can imbibe, i.e. I1, I2 and I3 [1]. Depending on the experiment conducted, the micro-models were initially filled with either DI water or air before the wetting fluid was injected. We found that the imbibition of the wetting fluid through a single pore is primarily controlled by the geometry of the pore body rather than the downstream pore throat sizes, contrary to the established capillary filling rules as used in current pore network models. Our experimental observations are confirmed by detailed lattice-Boltzmann pore scale computer simulations of fluid displacement in the same geometries. This suggests that capillary filling rules for imbibition as used in pore network models may need to be revised. [1] G. Lenormand, C. Zarcone and A. Sarr, J. Fluid Mech. 135 , 337-353 (1983).

Pore space analysis of NAPL distribution in sand-clay media

USGS Publications Warehouse

Matmon, D.; Hayden, N.J.

2003-01-01

This paper introduces a conceptual model of clays and non-aqueous phase liquids (NAPLs) at the pore scale that has been developed from a mathematical unit cell model, and direct micromodel observation and measurement of clay-containing porous media. The mathematical model uses a unit cell concept with uniform spherical grains for simulating the sand in the sand-clay matrix (???10% clay). Micromodels made with glass slides and including different clay-containing porous media were used to investigate the two clays (kaolinite and montmorillonite) and NAPL distribution within the pore space. The results were used to understand the distribution of NAPL advancing into initially saturated sand and sand-clay media, and provided a detailed analysis of the pore-scale geometry, pore size distribution, NAPL entry pressures, and the effect of clay on this geometry. Interesting NAPL saturation profiles were observed as a result of the complexity of the pore space geometry with the different packing angles and the presence of clays. The unit cell approach has applications for enhancing the mechanistic understanding and conceptualization, both visually and mathematically, of pore-scale processes such as NAPL and clay distribution. ?? 2003 Elsevier Science Ltd. All rights reserved.

Comparison of Pore-Network and Lattice Boltzmann Models for Pore-Scale Modeling of Geological Storage of CO2 in Natural Reservoir Rocks

NASA Astrophysics Data System (ADS)

Kohanpur, A. H.; Chen, Y.; Valocchi, A. J.; Tudek, J.; Crandall, D.

2016-12-01

CO2-brine flow in deep natural rocks is the focus of attention in geological storage of CO2. Understanding rock/flow properties at pore-scale is a vital component in field-scale modeling and prediction of fate of injected CO2. There are many challenges in working at the pore scale, such as size and selection of representative elementary volume (REV), particularly for material with complex geometry and heterogeneity, and the high computational costs. These issues factor into trade-offs that need to be made in choosing and applying pore-scale models. On one hand, pore-network modeling (PNM) simplifies the geometry and flow equations but can provide characteristic curves on fairly large samples. On the other hand, the lattice Boltzmann method (LBM) solves Navier-Stokes equations on the real geometry but is limited to small samples due to its high computational costs. Thus, both methods have some advantages but also face some challenges, which warrants a more detailed comparison and evaluation. In this study, we used industrial and micro-CT scans of actual reservoir rock samples to characterize pore structure at different resolutions. We ran LBM models directly on the characterized geometry and PNM on the equivalent 3D extracted network to determine single/two-phase flow properties during drainage and imbibition processes. Specifically, connectivity, absolute permeability, relative permeability curve, capillary pressure curve, and interface location are compared between models. We also did simulations on several subsamples from different locations including different domain sizes and orientations to encompass analysis of heterogeneity and isotropy. This work is primarily supported as part of the Center for Geologic Storage of CO2, an Energy Frontier Research Center funded by the U.S. Department of Energy, Office of Science and partially supported by the International Institute for Carbon-Neutral Energy Research (WPI-I2CNER) based at Kyushu University, Japan.

Importance of filter’s microstructure in dynamic filtration modeling of gasoline particulate filters (GPFs): Inhomogeneous porosity and pore size distribution

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gong, Jian; Stewart, Mark L.; Zelenyuk, Alla

The state-of-the-art multiscale modeling of GPFs including channel scale, wall scale, and pore scale is described. The microstructures of two GPFs were experimentally characterized. The pore size distributions of the GPFs were determined by mercury porosimetry. The porosity was measured by X-ray computed tomography (CT) and found to be inhomogeneous across the substrate wall. The significance of pore size distribution with respect to filtration performance was analyzed. The predictions of filtration efficiency were improved by including the pore size distribution in the filtration model. A dynamic heterogeneous multiscale filtration (HMF) model was utilized to simulate particulate filtration on a singlemore » channel particulate filter with realistic particulate emissions from a spark-ignition direct-injection (SIDI) gasoline engine. The dynamic evolution of filter’s microstructure and macroscopic filtration characteristics including mass- and number-based filtration efficiencies and pressure drop were predicted and discussed. The microstructure of the GPF substrate including inhomogeneous porosity and pore size distribution is found to significantly influence local particulate deposition inside the substrate and macroscopic filtration performance and is recommended to be resolved in the filtration model to simulate and evaluate the filtration performance of GPFs.« less

Importance of filter’s microstructure in dynamic filtration modeling of gasoline particulate filters (GPFs): Inhomogeneous porosity and pore size distribution

DOE PAGES

Gong, Jian; Stewart, Mark L.; Zelenyuk, Alla; ...

2018-01-03

The state-of-the-art multiscale modeling of gasoline particulate filter (GPF) including channel scale, wall scale, and pore scale is described. The microstructures of two GPFs were experimentally characterized. The pore size distributions of the GPFs were determined by mercury porosimetry. The porosity was measured by X-ray computed tomography (CT) and found to be inhomogeneous across the substrate wall. The significance of pore size distribution with respect to filtration performance was analyzed. The predictions of filtration efficiency were improved by including the pore size distribution in the filtration model. A dynamic heterogeneous multiscale filtration (HMF) model was utilized to simulate particulate filtrationmore » on a single channel particulate filter with realistic particulate emissions from a spark-ignition direct-injection (SIDI) gasoline engine. The dynamic evolution of filter’s microstructure and macroscopic filtration characteristics including mass- and number-based filtration efficiencies and pressure drop were predicted and discussed. In conclusion, the microstructure of the GPF substrate including inhomogeneous porosity and pore size distribution is found to significantly influence local particulate deposition inside the substrate and macroscopic filtration performance and is recommended to be resolved in the filtration model to simulate and evaluate the filtration performance of GPFs.« less

Importance of filter’s microstructure in dynamic filtration modeling of gasoline particulate filters (GPFs): Inhomogeneous porosity and pore size distribution

DOE Office of Scientific and Technical Information (OSTI.GOV)

Gong, Jian; Stewart, Mark L.; Zelenyuk, Alla

The state-of-the-art multiscale modeling of gasoline particulate filter (GPF) including channel scale, wall scale, and pore scale is described. The microstructures of two GPFs were experimentally characterized. The pore size distributions of the GPFs were determined by mercury porosimetry. The porosity was measured by X-ray computed tomography (CT) and found to be inhomogeneous across the substrate wall. The significance of pore size distribution with respect to filtration performance was analyzed. The predictions of filtration efficiency were improved by including the pore size distribution in the filtration model. A dynamic heterogeneous multiscale filtration (HMF) model was utilized to simulate particulate filtrationmore » on a single channel particulate filter with realistic particulate emissions from a spark-ignition direct-injection (SIDI) gasoline engine. The dynamic evolution of filter’s microstructure and macroscopic filtration characteristics including mass- and number-based filtration efficiencies and pressure drop were predicted and discussed. In conclusion, the microstructure of the GPF substrate including inhomogeneous porosity and pore size distribution is found to significantly influence local particulate deposition inside the substrate and macroscopic filtration performance and is recommended to be resolved in the filtration model to simulate and evaluate the filtration performance of GPFs.« less

Visualizing and Quantifying Pore Scale Fluid Flow Processes With X-ray Microtomography

NASA Astrophysics Data System (ADS)

Wildenschild, D.; Hopmans, J. W.; Vaz, C. M.; Rivers, M. L.

2001-05-01

When using mathematical models based on Darcy's law it is often necessary to simplify geometry, physics or both and the capillary bundle-of-tubes approach neglects a fundamentally important characteristic of porous solids, namely interconnectedness of the pore space. New approaches to pore-scale modeling that arrange capillary tubes in two- or three-dimensional pore space have been and are still under development: Network models generally represent the pore space by spheres while the pore throats are usually represented by cylinders or conical shapes. Lattice Boltzmann approaches numerically solve the Navier-Stokes equations in a realistic microscopically disordered geometry, which offers the ability to study the microphysical basis of macroscopic flow without the need for a simplified geometry or physics. In addition to these developments in numerical modeling techniques, new theories have proposed that interfacial area should be considered as a primary variable in modeling of a multi-phase flow system. In the wake of this progress emerges an increasing need for new ways of evaluating pore-scale models, and for techniques that can resolve and quantify phase interfaces in porous media. The mechanisms operating at the pore-scale cannot be measured with traditional experimental techniques, however x-ray computerized microtomography (CMT) provides non-invasive observation of, for instance, changing fluid phase content and distribution on the pore scale. Interfacial areas have thus far been measured indirectly, but with the advances in high-resolution imaging using CMT it is possible to track interfacial area and curvature as a function of phase saturation or capillary pressure. We present results obtained at the synchrotron-based microtomography facility (GSECARS, sector 13) at the Advanced Photon Source at Argonne National Laboratory. Cylindrical sand samples of either 6 or 1.5 mm diameter were scanned at different stages of drainage and for varying boundary conditions. A significant difference in fluid saturation and phase distribution was observed for different drainage conditions, clearly showing preferential flow and a dependence on the applied flow rate. For the 1.5 mm sample individual pores and water/air interfaces could be resolved and quantified using image analysis techniques. Use of the Advanced Photon Source was supported by the U.S. Department of Energy, Basic Energy Sciences, Office of Science, under Contract No. W-31-109-Eng-38.

Accounting for microbial habitats in modeling soil organic matter dynamics

NASA Astrophysics Data System (ADS)

Chenu, Claire; Garnier, Patricia; Nunan, Naoise; Pot, Valérie; Raynaud, Xavier; Vieublé, Laure; Otten, Wilfred; Falconer, Ruth; Monga, Olivier

2017-04-01

The extreme heterogeneity of soils constituents, architecture and inhabitants at the microscopic scale is increasingly recognized. Microbial communities exist and are active in a complex 3-D physical framework of mineral and organic particles defining pores of various sizes, more or less inter-connected. This results in a frequent spatial disconnection between soil carbon, energy sources and the decomposer organisms and a variety of microhabitats that are more or less suitable for microbial growth and activity. However, current biogeochemical models account for C dynamics at the macroscale (cm, m) and consider time- and spatially averaged relationships between microbial activity and soil characteristics. Different modelling approaches have intended to account for this microscale heterogeneity, based either on considering aggregates as surrogates for microbial habitats, or pores. Innovative modelling approaches are based on an explicit representation of soil structure at the fine scale, i.e. at µm to mm scales: pore architecture and their saturation with water, localization of organic resources and of microorganisms. Three recent models are presented here, that describe the heterotrophic activity of either bacteria or fungi and are based upon different strategies to represent the complex soil pore system (Mosaic, LBios and µFun). These models allow to hierarchize factors of microbial activity in soil's heterogeneous architecture. Present limits of these approaches and challenges are presented, regarding the extensive information required on soils at the microscale and to up-scale microbial functioning from the pore to the core scale.

FINAL REPORT: Mechanistically-Base Field Scale Models of Uranium Biogeochemistry from Upscaling Pore-Scale Experiments and Models

DOE Office of Scientific and Technical Information (OSTI.GOV)

Wood, Brian D.

2013-11-04

Biogeochemical reactive transport processes in the subsurface environment are important to many contemporary environmental issues of significance to DOE. Quantification of risks and impacts associated with environmental management options, and design of remediation systems where needed, require that we have at our disposal reliable predictive tools (usually in the form of numerical simulation models). However, it is well known that even the most sophisticated reactive transport models available today have poor predictive power, particularly when applied at the field scale. Although the lack of predictive ability is associated in part with our inability to characterize the subsurface and limitations inmore » computational power, significant advances have been made in both of these areas in recent decades and can be expected to continue. In this research, we examined the upscaling (pore to Darcy and Darcy to field) the problem of bioremediation via biofilms in porous media. The principle idea was to start with a conceptual description of the bioremediation process at the pore scale, and apply upscaling methods to formally develop the appropriate upscaled model at the so-called Darcy scale. The purpose was to determine (1) what forms the upscaled models would take, and (2) how one might parameterize such upscaled models for applications to bioremediation in the field. We were able to effectively upscale the bioremediation process to explain how the pore-scale phenomena were linked to the field scale. The end product of this research was to produce a set of upscaled models that could be used to help predict field-scale bioremediation. These models were mechanistic, in the sense that they directly incorporated pore-scale information, but upscaled so that only the essential features of the process were needed to predict the effective parameters that appear in the model. In this way, a direct link between the microscale and the field scale was made, but the upscaling process helped inform potential users of the model what kinds of information would be needed to accurately characterize the system.« less

A Generalized Hybrid Multiscale Modeling Approach for Flow and Reactive Transport in Porous Media

NASA Astrophysics Data System (ADS)

Yang, X.; Meng, X.; Tang, Y. H.; Guo, Z.; Karniadakis, G. E.

2017-12-01

Using emerging understanding of biological and environmental processes at fundamental scales to advance predictions of the larger system behavior requires the development of multiscale approaches, and there is strong interest in coupling models at different scales together in a hybrid multiscale simulation framework. A limited number of hybrid multiscale simulation methods have been developed for subsurface applications, mostly using application-specific approaches for model coupling. The proposed generalized hybrid multiscale approach is designed with minimal intrusiveness to the at-scale simulators (pre-selected) and provides a set of lightweight C++ scripts to manage a complex multiscale workflow utilizing a concurrent coupling approach. The workflow includes at-scale simulators (using the lattice-Boltzmann method, LBM, at the pore and Darcy scale, respectively), scripts for boundary treatment (coupling and kriging), and a multiscale universal interface (MUI) for data exchange. The current study aims to apply the generalized hybrid multiscale modeling approach to couple pore- and Darcy-scale models for flow and mixing-controlled reaction with precipitation/dissolution in heterogeneous porous media. The model domain is packed heterogeneously that the mixing front geometry is more complex and not known a priori. To address those challenges, the generalized hybrid multiscale modeling approach is further developed to 1) adaptively define the locations of pore-scale subdomains, 2) provide a suite of physical boundary coupling schemes and 3) consider the dynamic change of the pore structures due to mineral precipitation/dissolution. The results are validated and evaluated by comparing with single-scale simulations in terms of velocities, reactive concentrations and computing cost.

Assessing the utility of FIB-SEM images for shale digital rock physics

NASA Astrophysics Data System (ADS)

Kelly, Shaina; El-Sobky, Hesham; Torres-Verdín, Carlos; Balhoff, Matthew T.

2016-09-01

Shales and other unconventional or low permeability (tight) reservoirs house vast quantities of hydrocarbons, often demonstrate considerable water uptake, and are potential repositories for fluid sequestration. The pore-scale topology and fluid transport mechanisms within these nanoporous sedimentary rocks remain to be fully understood. Image-informed pore-scale models are useful tools for studying porous media: a debated question in shale pore-scale petrophysics is whether there is a representative elementary volume (REV) for shale models? Furthermore, if an REV exists, how does it differ among petrophysical properties? We obtain three dimensional (3D) models of the topology of microscale shale volumes from image analysis of focused ion beam-scanning electron microscope (FIB-SEM) image stacks and investigate the utility of these models as a potential REV for shale. The scope of data used in this work includes multiple local groups of neighboring FIB-SEM images of different microscale sizes, corresponding core-scale (milli- and centimeters) laboratory data, and, for comparison, series of two-dimensional (2D) cross sections from broad ion beam SEM images (BIB-SEM), which capture a larger microscale field of view than the FIB-SEM images; this array of data is larger than the majority of investigations with FIB-SEM-derived microscale models of shale. Properties such as porosity, organic matter content, and pore connectivity are extracted from each model. Assessments of permeability with single phase, pressure-driven flow simulations are performed in the connected pore space of the models using the lattice-Boltzmann method. Calculated petrophysical properties are compared to those of neighboring FIB-SEM images and to core-scale measurements of the sample associated with the FIB-SEM sites. Results indicate that FIB-SEM images below ∼5000 μm3 volume (the largest volume analyzed) are not a suitable REV for shale permeability and pore-scale networks; i.e. field of view is compromised at the expense of detailed, but often unconnected, nanopore morphology. Further, we find that it is necessary to acquire several local FIB-SEM or BIB-SEM images and correlate their extracted geometric properties to improve the likelihood of achieving representative values of porosity and organic matter volume. Our work indicates that FIB-SEM images of microscale volumes of shale are a qualitative tool for petrophysical and transport analysis. Finally, we offer alternatives for quantitative pore-scale assessments of shale.

Direct pore-scale reactive transport modelling of dynamic wettability changes induced by surface complexation

NASA Astrophysics Data System (ADS)

Maes, Julien; Geiger, Sebastian

2018-01-01

Laboratory experiments have shown that oil production from sandstone and carbonate reservoirs by waterflooding could be significantly increased by manipulating the composition of the injected water (e.g. by lowering the ionic strength). Recent studies suggest that a change of wettability induced by a change in surface charge is likely to be one of the driving mechanism of the so-called low-salinity effect. In this case, the potential increase of oil recovery during waterflooding at low ionic strength would be strongly impacted by the inter-relations between flow, transport and chemical reaction at the pore-scale. Hence, a new numerical model that includes two-phase flow, solute reactive transport and wettability alteration is implemented based on the Direct Numerical Simulation of the Navier-Stokes equations and surface complexation modelling. Our model is first used to match experimental results of oil droplet detachment from clay patches. We then study the effect of wettability change on the pore-scale displacement for simple 2D calcite micro-models and evaluate the impact of several parameters such as water composition and injected velocity. Finally, we repeat the simulation experiments on a larger and more complex pore geometry representing a carbonate rock. Our simulations highlight two different effects of low-salinity on oil production from carbonate rocks: a smaller number of oil clusters left in the pores after invasion, and a greater number of pores invaded.

Effects of pore-scale physics on uranium geochemistry in Hanford sediments

DOE Office of Scientific and Technical Information (OSTI.GOV)

Hu, Qinhong; Ewing, Robert P.

Overall, this work examines a key scientific issue, mass transfer limitations at the pore-scale, using both new instruments with high spatial resolution, and new conceptual and modeling paradigms. The complementary laboratory and numerical approaches connect pore-scale physics to macroscopic measurements, providing a previously elusive scale integration. This Exploratory research project produced five peer-reviewed journal publications and eleven scientific presentations. This work provides new scientific understanding, allowing the DOE to better incorporate coupled physical and chemical processes into decision making for environmental remediation and long-term stewardship.

Micro/Nano-pore Network Analysis of Gas Flow in Shale Matrix

PubMed Central

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N.; Gao, Shengyan

2015-01-01

The gas flow in shale matrix is of great research interests for optimized shale gas extraction. The gas flow in the nano-scale pore may fall in flow regimes such as viscous flow, slip flow and Knudsen diffusion. A 3-dimensional nano-scale pore network model was developed to simulate dynamic gas flow, and to describe the transient properties of flow regimes. The proposed pore network model accounts for the various size distributions and low connectivity of shale pores. The pore size, pore throat size and coordination number obey normal distribution, and the average values can be obtained from shale reservoir data. The gas flow regimes were simulated using an extracted pore network backbone. The numerical results show that apparent permeability is strongly dependent on pore pressure in the reservoir and pore throat size, which is overestimated by low-pressure laboratory tests. With the decrease of reservoir pressure, viscous flow is weakening, then slip flow and Knudsen diffusion are gradually becoming dominant flow regimes. The fingering phenomenon can be predicted by micro/nano-pore network for gas flow, which provides an effective way to capture heterogeneity of shale gas reservoir. PMID:26310236

Micro/Nano-pore Network Analysis of Gas Flow in Shale Matrix.

PubMed

Zhang, Pengwei; Hu, Liming; Meegoda, Jay N; Gao, Shengyan

2015-08-27

The gas flow in shale matrix is of great research interests for optimized shale gas extraction. The gas flow in the nano-scale pore may fall in flow regimes such as viscous flow, slip flow and Knudsen diffusion. A 3-dimensional nano-scale pore network model was developed to simulate dynamic gas flow, and to describe the transient properties of flow regimes. The proposed pore network model accounts for the various size distributions and low connectivity of shale pores. The pore size, pore throat size and coordination number obey normal distribution, and the average values can be obtained from shale reservoir data. The gas flow regimes were simulated using an extracted pore network backbone. The numerical results show that apparent permeability is strongly dependent on pore pressure in the reservoir and pore throat size, which is overestimated by low-pressure laboratory tests. With the decrease of reservoir pressure, viscous flow is weakening, then slip flow and Knudsen diffusion are gradually becoming dominant flow regimes. The fingering phenomenon can be predicted by micro/nano-pore network for gas flow, which provides an effective way to capture heterogeneity of shale gas reservoir.

A streamline splitting pore-network approach for computationally inexpensive and accurate simulation of transport in porous media

DOE Office of Scientific and Technical Information (OSTI.GOV)

Mehmani, Yashar; Oostrom, Martinus; Balhoff, Matthew

2014-03-20

Several approaches have been developed in the literature for solving flow and transport at the pore-scale. Some authors use a direct modeling approach where the fundamental flow and transport equations are solved on the actual pore-space geometry. Such direct modeling, while very accurate, comes at a great computational cost. Network models are computationally more efficient because the pore-space morphology is approximated. Typically, a mixed cell method (MCM) is employed for solving the flow and transport system which assumes pore-level perfect mixing. This assumption is invalid at moderate to high Peclet regimes. In this work, a novel Eulerian perspective on modelingmore » flow and transport at the pore-scale is developed. The new streamline splitting method (SSM) allows for circumventing the pore-level perfect mixing assumption, while maintaining the computational efficiency of pore-network models. SSM was verified with direct simulations and excellent matches were obtained against micromodel experiments across a wide range of pore-structure and fluid-flow parameters. The increase in the computational cost from MCM to SSM is shown to be minimal, while the accuracy of SSM is much higher than that of MCM and comparable to direct modeling approaches. Therefore, SSM can be regarded as an appropriate balance between incorporating detailed physics and controlling computational cost. The truly predictive capability of the model allows for the study of pore-level interactions of fluid flow and transport in different porous materials. In this paper, we apply SSM and MCM to study the effects of pore-level mixing on transverse dispersion in 3D disordered granular media.« less

Multiscale modeling of fluid flow and mass transport

NASA Astrophysics Data System (ADS)

Masuoka, K.; Yamamoto, H.; Bijeljic, B.; Lin, Q.; Blunt, M. J.

2017-12-01

In recent years, there are some reports on a simulation of fluid flow in pore spaces of rocks using Navier-Stokes equations. These studies mostly adopt a X-ray CT to create 3-D numerical grids of the pores in micro-scale. However, results may be of low accuracy when the rock has a large pore size distribution, because pores, whose size is smaller than resolution of the X-ray CT may be neglected. We recently found out by tracer tests in a laboratory using a brine saturated Ryukyu limestone and inject fresh water that a decrease of chloride concentration took longer time. This phenomenon can be explained due to weak connectivity of the porous networks. Therefore, it is important to simulate entire pore spaces even those of very small sizes in which diffusion is dominant. We have developed a new methodology for multi-level modeling for pore scale fluid flow in porous media. The approach is to combine pore-scale analysis with Darcy-flow analysis using two types of X-ray CT images in different resolutions. Results of the numerical simulations showed a close match with the experimental results. The proposed methodology is an enhancement for analyzing mass transport and flow phenomena in rocks with complicated pore structure.

Multiscale modelling of dual-porosity porous media; a computational pore-scale study for flow and solute transport

NASA Astrophysics Data System (ADS)

de Vries, Enno T.; Raoof, Amir; van Genuchten, Martinus Th.

2017-07-01

Many environmental and agricultural applications involve the transport of water and dissolved constituents through aggregated soil profiles, or porous media that are structured, fractured or macroporous in other ways. During the past several decades, various process-based macroscopic models have been used to simulate contaminant transport in such media. Many of these models consider advective-dispersive transport through relatively large inter-aggregate pore domains, while exchange with the smaller intra-aggregate pores is assumed to be controlled by diffusion. Exchange of solute between the two domains is often represented using a first-order mass transfer coefficient, which is commonly obtained by fitting to observed data. This study aims to understand and quantify the solute exchange term by applying a dual-porosity pore-scale network model to relatively large domains, and analysing the pore-scale results in terms of the classical dual-porosity (mobile-immobile) transport formulation. We examined the effects of key parameters (notably aggregate porosity and aggregate permeability) on the main dual-porosity model parameters, i.e., the mobile water fraction (ϕm) and the mass transfer coefficient (α). Results were obtained for a wide range of aggregate porosities (between 0.082 and 0.700). The effect of aggregate permeability was explored by varying pore throat sizes within the aggregates. Solute breakthrough curves (BTCs) obtained with the pore-scale network model at several locations along the domain were analysed using analytical solutions of the dual-porosity model to obtain estimates of ϕm and α. An increase in aggregate porosity was found to decrease ϕm and increase α, leading to considerable tailing in the BTCs. Changes in the aggregate pore throat size affected the relative flow velocity between the intra- and inter-aggregate domains. Higher flow velocities within the aggregates caused a change in the transport regime from diffusion dominated to more advection dominated. This change increased the exchange rate of solutes between the mobile and immobile domains, with a related increase in the value of the mass transfer coefficient and less tailing in the BTCs.

Dispersion upscaling from a pore scale characterization of Lagrangian velocities

NASA Astrophysics Data System (ADS)

Turuban, Régis; de Anna, Pietro; Jiménez-Martínez, Joaquín; Tabuteau, Hervé; Méheust, Yves; Le Borgne, Tanguy

2013-04-01

Mixing and reactive transport are primarily controlled by the interplay between diffusion, advection and reaction at pore scale. Yet, how the distribution and spatial correlation of the velocity field at pore scale impact these processes is still an open question. Here we present an experimental investigation of the distribution and correlation of pore scale velocities and its relation with upscaled dispersion. We use a quasi two-dimensional (2D) horizontal set up, consisting of two glass plates filled with cylinders representing the grains of the porous medium : the cell is built by soft lithography technique, wich allows for full control of the system geometry. The local velocity field is quantified from particle tracking velocimetry using microspheres that are advected with the pore scale flow. Their displacement is purely advective, as the particle size is chosen large enough to avoid diffusion. We thus obtain particle trajectories as well as lagrangian velocities in the entire system. The measured velocity field shows the existence of a network of preferential flow paths in channels with high velocities, as well as very low velocity in stagnation zones, with a non Gaussian distribution. Lagrangian velocities are long range correlated in time, which implies a non-fickian scaling of the longitudinal variance of particle positions. To upscale this process we develop an effective transport model, based on correlated continous time random walk, which is entirely parametrized by the pore scale velocity distribution and correlation. The model predictions are compared with conservative tracer test data for different Peclet numbers. Furthermore, we investigate the impact of different pore geometries on the distribution and correlation of Lagrangian velocities and we discuss the link between these properties and the effective dispersion behavior.

A level set method for determining critical curvatures for drainage and imbibition.

PubMed

Prodanović, Masa; Bryant, Steven L

2006-12-15

An accurate description of the mechanics of pore level displacement of immiscible fluids could significantly improve the predictions from pore network models of capillary pressure-saturation curves, interfacial areas and relative permeability in real porous media. If we assume quasi-static displacement, at constant pressure and surface tension, pore scale interfaces are modeled as constant mean curvature surfaces, which are not easy to calculate. Moreover, the extremely irregular geometry of natural porous media makes it difficult to evaluate surface curvature values and corresponding geometric configurations of two fluids. Finally, accounting for the topological changes of the interface, such as splitting or merging, is nontrivial. We apply the level set method for tracking and propagating interfaces in order to robustly handle topological changes and to obtain geometrically correct interfaces. We describe a simple but robust model for determining critical curvatures for throat drainage and pore imbibition. The model is set up for quasi-static displacements but it nevertheless captures both reversible and irreversible behavior (Haines jump, pore body imbibition). The pore scale grain boundary conditions are extracted from model porous media and from imaged geometries in real rocks. The method gives quantitative agreement with measurements and with other theories and computational approaches.

Hybrid discrete-continuum modeling for transport, biofilm development and solid restructuring including electrostatic effects

NASA Astrophysics Data System (ADS)

Prechtel, Alexander; Ray, Nadja; Rupp, Andreas

2017-04-01

We want to present an approach for the mathematical, mechanistic modeling and numerical treatment of processes leading to the formation, stability, and turnover of soil micro-aggregates. This aims at deterministic aggregation models including detailed mechanistic pore-scale descriptions to account for the interplay of geochemistry and microbiology, and the link to soil functions as, e.g., the porosity. We therefore consider processes at the pore scale and the mesoscale (laboratory scale). At the pore scale transport by diffusion, advection, and drift emerging from electric forces can be taken into account, in addition to homogeneous and heterogeneous reactions of species. In the context of soil micro-aggregates the growth of biofilms or other glueing substances as EPS (extracellular polymeric substances) is important and affects the structure of the pore space in space and time. This model is upscaled mathematically in the framework of (periodic) homogenization to transfer it to the mesoscale resulting in effective coefficients/parameters there. This micro-macro model thus couples macroscopic equations that describe the transport and fluid flow at the scale of the porous medium (mesoscale) with averaged time- and space-dependent coefficient functions. These functions may be explicitly computed by means of auxiliary cell problems (microscale). Finally, the pore space in which the cell problems are defined is time and space dependent and its geometry inherits information from the transport equation's solutions. The microscale problems rely on versatile combinations of cellular automata and discontiuous Galerkin methods while on the mesoscale mixed finite elements are used. The numerical simulations allow to study the interplay between these processes.

Pores-scale hydrodynamics in a progressively bio-clogged three-dimensional porous medium: 3D particle tracking experiments and stochastic transport modelling

NASA Astrophysics Data System (ADS)

Morales, V. L.; Carrel, M.; Dentz, M.; Derlon, N.; Morgenroth, E.; Holzner, M.

2017-12-01

Biofilms are ubiquitous bacterial communities growing in various porous media including soils, trickling and sand filters and are relevant for applications such as the degradation of pollutants for bioremediation, waste water or drinking water production purposes. By their development, biofilms dynamically change the structure of porous media, increasing the heterogeneity of the pore network and the non-Fickian or anomalous dispersion. In this work, we use an experimental approach to investigate the influence of biofilm growth on pore scale hydrodynamics and transport processes and propose a correlated continuous time random walk model capturing these observations. We perform three-dimensional particle tracking velocimetry at four different time points from 0 to 48 hours of biofilm growth. The biofilm growth notably impacts pore-scale hydrodynamics, as shown by strong increase of the average velocity and in tailing of Lagrangian velocity probability density functions. Additionally, the spatial correlation length of the flow increases substantially. This points at the formation of preferential flow pathways and stagnation zones, which ultimately leads to an increase of anomalous transport in the porous media considered, characterized by non-Fickian scaling of mean-squared displacements and non-Gaussian distributions of the displacement probability density functions. A gamma distribution provides a remarkable approximation of the bulk and the high tail of the Lagrangian pore-scale velocity magnitude, indicating a transition from a parallel pore arrangement towards a more serial one. Finally, a correlated continuous time random walk based on a stochastic relation velocity model accurately reproduces the observations and could be used to predict transport beyond the time scales accessible to the experiment.

Pore-scale modeling of Capillary Penetration of Wetting Liquid into 3D Fibrous Media: A Critical Examination of Equivalent Capillary Concept

NASA Astrophysics Data System (ADS)

Palakurthi, Nikhil Kumar; Ghia, Urmila; Comer, Ken

2013-11-01

Capillary penetration of liquid through fibrous porous media is important in many applications such as printing, drug delivery patches, sanitary wipes, and performance fabrics. Historically, capillary transport (with a distinct liquid propagating front) in porous media is modeled using capillary-bundle theory. However, it is not clear if the capillary model (Washburn equation) describes the fluid transport in porous media accurately, as it assumes uniformity of pore sizes in the porous medium. The present work investigates the limitations of the applicability of the capillary model by studying liquid penetration through virtual fibrous media with uniform and non-uniform pore-sizes. For the non-uniform-pore fibrous medium, the effective capillary radius of the fibrous medium was estimated from the pore-size distribution curve. Liquid penetration into the 3D virtual fibrous medium at micro-scale was simulated using OpenFOAM, and the numerical results were compared with the Washburn-equation capillary-model predictions. Preliminary results show that the Washburn equation over-predicts the height rise in the early stages (purely inertial and visco-inertial stages) of capillary transport.

Impact of spatially correlated pore-scale heterogeneity on drying porous media

NASA Astrophysics Data System (ADS)

Borgman, Oshri; Fantinel, Paolo; Lühder, Wieland; Goehring, Lucas; Holtzman, Ran

2017-07-01

We study the effect of spatially-correlated heterogeneity on isothermal drying of porous media. We combine a minimal pore-scale model with microfluidic experiments with the same pore geometry. Our simulated drying behavior compares favorably with experiments, considering the large sensitivity of the emergent behavior to the uncertainty associated with even small manufacturing errors. We show that increasing the correlation length in particle sizes promotes preferential drying of clusters of large pores, prolonging liquid connectivity and surface wetness and thus higher drying rates for longer periods. Our findings improve our quantitative understanding of how pore-scale heterogeneity impacts drying, which plays a role in a wide range of processes ranging from fuel cells to curing of paints and cements to global budgets of energy, water and solutes in soils.

Measurements of pore-scale flow through apertures

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chojnicki, Kirsten

Pore-scale aperture effects on flow in pore networks was studied in the laboratory to provide a parameterization for use in transport models. Four cases were considered: regular and irregular pillar/pore alignment with and without an aperture. The velocity field of each case was measured and simulated, providing quantitatively comparable results. Two aperture effect parameterizations were considered: permeability and transmission. Permeability values varied by an order of magnitude between the cases with and without apertures. However, transmission did not correlate with permeability. Despite having much greater permeability the regular aperture case permitted less transmission than the regular case. Moreover, both irregularmore » cases had greater transmission than the regular cases, a difference not supported by the permeabilities. Overall, these findings suggest that pore-scale aperture effects on flow though a pore-network may not be adequately captured by properties such as permeability for applications that are interested in determining particle transport volume and timing.« less

Topological characteristics underpin intermittency and anomalous transport behavior in soil-like porous media

NASA Astrophysics Data System (ADS)

Holzner, M.; Morales, V.; Willmann, M.; Jerjen, I.; Kaufmann, R.; Dentz, M.

2016-12-01

Continuum models of porous media are based on the validity of the Darcy equation for fluid and Fick's law for scalar fluxes on a representative elementary volume. Fluctuations of pore-scale flow and scalar transport are averaged out and represented in terms of effective parameters such as hydrodynamic dispersion. However, the intermittent behavior of pore-scale flow impacts on the nature of particle and scalar transport, and it determines the way dissolved substances mix and react. The understanding of the origin of these processes is of both fundamental and practical importance in applications ranging from reactive transport in groundwater flow to diffusion in fuel cells or biological systems. A central issue in porous medium flow is therefore to relate intermittent behavior of Lagrangian velocity at pore scale imposed by the complex pore network geometry to transport properties at larger scales. Lagrangian measurements in porous systems are nonetheless scarce and most experimental techniques do not provide access to all three velocity components. In this contribution we report 3D measurements of Lagrangian velocity in soil-like porous media. We complement these measurements with detailed X-ray scans of the pore network. We find sharp velocity transitions close to pore throats, and low flow variability in the pore bodies, which gives rise to stretched exponential Lagrangian velocity and acceleration distributions characterized by a sharp peak at low velocity and a superlinear evolution of particle dispersion. We demonstrate that porosity and pore size distribution alone cannot explain the observed features of the flow. Rather, anomalous transport is better interpreted in terms of how pores of various geometries are interconnected. We reproduce the main observations using a continuous-time random walk (CTRW) model revealing the main features that control the system and showing the potential of this simple model to capture transport in complex geometries.

Persistent Homology fingerprinting of microstructural controls on larger-scale fluid flow in porous media

NASA Astrophysics Data System (ADS)

Moon, C.; Mitchell, S. A.; Callor, N.; Dewers, T. A.; Heath, J. E.; Yoon, H.; Conner, G. R.

2017-12-01

Traditional subsurface continuum multiphysics models include useful yet limiting geometrical assumptions: penny- or disc-shaped cracks, spherical or elliptical pores, bundles of capillary tubes, cubic law fracture permeability, etc. Each physics (flow, transport, mechanics) uses constitutive models with an increasing number of fit parameters that pertain to the microporous structure of the rock, but bear no inter-physics relationships or self-consistency. Recent advances in digital rock physics and pore-scale modeling link complex physics to detailed pore-level geometries, but measures for upscaling are somewhat unsatisfactory and come at a high computational cost. Continuum mechanics rely on a separation between small scale pore fluctuations and larger scale heterogeneity (and perhaps anisotropy), but this can break down (particularly for shales). Algebraic topology offers powerful mathematical tools for describing a local-to-global structure of shapes. Persistent homology, in particular, analyzes the dynamics of topological features and summarizes into numeric values. It offers a roadmap to both "fingerprint" topologies of pore structure and multiscale connectedness as well as links pore structure to physical behavior, thus potentially providing a means to relate the dependence of constitutive behaviors of pore structures in a self-consistent way. We present a persistence homology (PH) analysis framework of 3D image sets including a focused ion beam-scanning electron microscopy data set of the Selma Chalk. We extract structural characteristics of sampling volumes via persistence homology and fit a statistical model using the summarized values to estimate porosity, permeability, and connectivity—Lattice Boltzmann methods for single phase flow modeling are used to obtain the relationships. These PH methods allow for prediction of geophysical properties based on the geometry and connectivity in a computationally efficient way. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia LLC, a wholly owned subsidiary of Honeywell International Inc. for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

EVALUATION AND ANALYSIS OF MICROSCALE FLOW AND TRANSPORT DURING REMEDIATION

EPA Science Inventory

The design of in-situ remediation is currently based on a description at the macroscopic scale. Phenomena at the pore and pore-network scales are typically lumped in terms of averaged quantities, using empirical or ad hoc expressions. These models cannot address fundamental rem...

Pore scale study of multiphase multicomponent reactive transport during CO 2 dissolution trapping

DOE Office of Scientific and Technical Information (OSTI.GOV)

Chen, Li; Wang, Mengyi; Kang, Qinjun

Solubility trapping is crucial for permanent CO 2 sequestration in deep saline aquifers. For the first time, a pore-scale numerical method is developed to investigate coupled scCO 2-water two-phase flow, multicomponent (CO 2(aq), H +, HCO 3 –, CO 3 2 – and OH –) mass transport, heterogeneous interfacial dissolution reaction, and homogeneous dissociation reactions. Pore-scale details of evolutions of multiphase distributions and concentration fields are presented and discussed. Time evolutions of several variables including averaged CO 2(aq) concentration, scCO 2 saturation, and pH value are analyzed. Specific interfacial length, an important variable which cannot be determined but is requiredmore » by continuum models, is investigated in detail. Mass transport coefficient or efficient dissolution rate is also evaluated. The pore-scale results show strong non-equilibrium characteristics during solubility trapping due to non-uniform distributions of multiphase as well as slow mass transport process. Complicated coupling mechanisms between multiphase flow, mass transport and chemical reactions are also revealed. Lastly, effects of wettability are also studied. The pore-scale studies provide deep understanding of non-linear non-equilibrium multiple physicochemical processes during CO 2 solubility trapping processes, and also allow to quantitatively predict some important empirical relationships, such as saturation-interfacial surface area, for continuum models.« less

Pore scale study of multiphase multicomponent reactive transport during CO 2 dissolution trapping

DOE PAGES

Chen, Li; Wang, Mengyi; Kang, Qinjun; ...

2018-04-26

Solubility trapping is crucial for permanent CO 2 sequestration in deep saline aquifers. For the first time, a pore-scale numerical method is developed to investigate coupled scCO 2-water two-phase flow, multicomponent (CO 2(aq), H +, HCO 3 –, CO 3 2 – and OH –) mass transport, heterogeneous interfacial dissolution reaction, and homogeneous dissociation reactions. Pore-scale details of evolutions of multiphase distributions and concentration fields are presented and discussed. Time evolutions of several variables including averaged CO 2(aq) concentration, scCO 2 saturation, and pH value are analyzed. Specific interfacial length, an important variable which cannot be determined but is requiredmore » by continuum models, is investigated in detail. Mass transport coefficient or efficient dissolution rate is also evaluated. The pore-scale results show strong non-equilibrium characteristics during solubility trapping due to non-uniform distributions of multiphase as well as slow mass transport process. Complicated coupling mechanisms between multiphase flow, mass transport and chemical reactions are also revealed. Lastly, effects of wettability are also studied. The pore-scale studies provide deep understanding of non-linear non-equilibrium multiple physicochemical processes during CO 2 solubility trapping processes, and also allow to quantitatively predict some important empirical relationships, such as saturation-interfacial surface area, for continuum models.« less

Pore scale study of multiphase multicomponent reactive transport during CO2 dissolution trapping

NASA Astrophysics Data System (ADS)

Chen, Li; Wang, Mengyi; Kang, Qinjun; Tao, Wenquan

2018-06-01

Solubility trapping is crucial for permanent CO2 sequestration in deep saline aquifers. For the first time, a pore-scale numerical method is developed to investigate coupled scCO2-water two-phase flow, multicomponent (CO2(aq), H+, HCO3-, CO32- and OH-) mass transport, heterogeneous interfacial dissolution reaction, and homogeneous dissociation reactions. Pore-scale details of evolutions of multiphase distributions and concentration fields are presented and discussed. Time evolutions of several variables including averaged CO2(aq) concentration, scCO2 saturation, and pH value are analyzed. Specific interfacial length, an important variable which cannot be determined but is required by continuum models, is investigated in detail. Mass transport coefficient or efficient dissolution rate is also evaluated. The pore-scale results show strong non-equilibrium characteristics during solubility trapping due to non-uniform distributions of multiphase as well as slow mass transport process. Complicated coupling mechanisms between multiphase flow, mass transport and chemical reactions are also revealed. Finally, effects of wettability are also studied. The pore-scale studies provide deep understanding of non-linear non-equilibrium multiple physicochemical processes during CO2 solubility trapping processes, and also allow to quantitatively predict some important empirical relationships, such as saturation-interfacial surface area, for continuum models.

Pore-scale modeling of capillary trapping in water-wet porous media: A new cooperative pore-body filling model

NASA Astrophysics Data System (ADS)

Ruspini, L. C.; Farokhpoor, R.; Øren, P. E.

2017-10-01

We present a pore-network model study of capillary trapping in water-wet porous media. The amount and distribution of trapped non-wetting phase is determined by the competition between two trapping mechanisms - snap-off and cooperative pore-body filling. We develop a new model to describe the pore-body filling mechanism in geologically realistic pore-networks. The model accounts for the geometrical characteristics of the pore, the spatial location of the connecting throats and the local fluid topology at the time of the displacement. We validate the model by comparing computed capillary trapping curves with published data for four different water-wet rocks. Computations are performed on pore-networks extracted from micro-CT images and process-based reconstructions of the actual rocks used in the experiments. Compared with commonly used stochastic models, the new model describes more accurately the experimental measurements, especially for well connected porous systems where trapping is controlled by subtleties of the pore structure. The new model successfully predicts relative permeabilities and residual saturation for Bentheimer sandstone using in-situ measured contact angles as input to the simulations. The simulated trapped cluster size distributions are compared with predictions from percolation theory.

Quantifying uncertainty and computational complexity for pore-scale simulations

NASA Astrophysics Data System (ADS)

Chen, C.; Yuan, Z.; Wang, P.; Yang, X.; Zhenyan, L.

2016-12-01

Pore-scale simulation is an essential tool to understand the complex physical process in many environmental problems, from multi-phase flow in the subsurface to fuel cells. However, in practice, factors such as sample heterogeneity, data sparsity and in general, our insufficient knowledge of the underlying process, render many simulation parameters and hence the prediction results uncertain. Meanwhile, most pore-scale simulations (in particular, direct numerical simulation) incur high computational cost due to finely-resolved spatio-temporal scales, which further limits our data/samples collection. To address those challenges, we propose a novel framework based on the general polynomial chaos (gPC) and build a surrogate model representing the essential features of the underlying system. To be specific, we apply the novel framework to analyze the uncertainties of the system behavior based on a series of pore-scale numerical experiments, such as flow and reactive transport in 2D heterogeneous porous media and 3D packed beds. Comparing with recent pore-scale uncertainty quantification studies using Monte Carlo techniques, our new framework requires fewer number of realizations and hence considerably reduce the overall computational cost, while maintaining the desired accuracy.

Pore-scale Modeling of CO2 Local Trapping in Heterogeneous Porous Media with Inter-granular Cements

NASA Astrophysics Data System (ADS)

Wang, D.; Li, Y.

2017-12-01

Based on pore-scale modeling of CO2/brine multiphase flow in heterogeneous porous media with inter-granular cements, we numerically analyze the effects of cement-modified pore structure on CO2 local trapping. Results indicate: 1) small pore throat is the main reason for causing CO2 local trapping in front of low-porosity layers (namely dense layers) formed by inter-granular cements; 2) in the case of the same pore throat size, the smaller particle size can increase the number of flow paths for CO2 plume and equivalently enhances local permeability, which may counteract the impediment of high capillary pressure on CO2 migration to some extent and consequently disables CO2 local capillary trapping; 3) the isolated pores by inter-granular cements can lead to dramatic reduction of CO2 saturation inside the dense layers, whereas the change of connectivity of some pores due to the cements can increase CO2 accumulation in front of the dense layers by lowering the displacement area of CO2 plume.

Unifying Pore Network Modeling, Continuous Time Random Walk Theory and Experiment - Accomplishments and Future Directions

NASA Astrophysics Data System (ADS)

Bijeljic, B.

2008-05-01

This talk will describe and highlight the advantages offered by a methodology that unifies pore network modeling, CTRW theory and experiment in description of solute dispersion in porous media. Solute transport in a porous medium is characterized by the interplay of advection and diffusion (described by Peclet number, Pe) that cause spreading of solute particles. This spreading is traditionally described by dispersion coefficients, D, defined by σ 2 = 2Dt, where σ 2 is the variance of the solute position and t is the time. Using a pore-scale network model based on particle tracking, the rich Peclet- number dependence of dispersion coefficient is predicted from first principles and is shown to compare well with experimental data for restricted diffusion, transition, power-law and mechanical dispersion regimes in the asymptotic limit. In the asymptotic limit D is constant and can be used in an averaged advection-dispersion equation. However, it is highly important to recognize that, until the velocity field is fully sampled, the particle transport is non-Gaussian and D possesses temporal or spatial variation. Furthermore, temporal probability density functions (PDF) of tracer particles are studied in pore networks and an excellent agreement for the spectrum of transition times for particles from pore to pore is obtained between network model results and CTRW theory. Based on the truncated power-law interpretation of PDF-s, the physical origin of the power-law scaling of dispersion coefficient vs. Peclet number has been explained for unconsolidated porous media, sands and a number of sandstones, arriving at the same conclusion from numerical network modelling, analytic CTRW theory and experiment. Future directions for further applications of the methodology presented are discussed in relation to the scale- dependent solute dispersion and reactive transport. Significance of pre-asymptotic dispersion in porous media is addressed from pore-scale upwards and the impact of heterogeneity is discussed. The length traveled by solute plumes before Gaussian behaviour is reached increases with an increase in heterogeneity and/or Pe. This opens up the question on the nature of dispersion in natural systems where the heterogeneities at the larger scales will profoundly increase the range of velocities in the aquifer, thus considerably delaying the asymptotic approach to Gaussian behaviour. As a consequence, the asymptotic behaviour might not be reached at the field scale.

Unifying Pore Network Modeling, Continuous Time Random Walk (CTRW) Theory and Experiment to Describe Impact of Spatial Heterogeneities on Solute Dispersion at Multiple Length-scales

NASA Astrophysics Data System (ADS)

Bijeljic, B.; Blunt, M. J.; Rhodes, M. E.

2009-04-01

This talk will describe and highlight the advantages offered by a novel methodology that unifies pore network modeling, CTRW theory and experiment in description of solute dispersion in porous media. Solute transport in a porous medium is characterized by the interplay of advection and diffusion (described by Peclet number, Pe) that cause dispersion of solute particles. Dispersion is traditionally described by dispersion coefficients, D, that are commonly calculated from the spatial moments of the plume. Using a pore-scale network model based on particle tracking, the rich Peclet-number dependence of dispersion coefficient is predicted from first principles and is shown to compare well with experimental data for restricted diffusion, transition, power-law and mechanical dispersion regimes in the asymptotic limit. In the asymptotic limit D is constant and can be used in an averaged advection-dispersion equation. However, it is highly important to recognize that, until the velocity field is fully sampled, the particle transport is non-Gaussian and D possesses temporal or spatial variation. Furthermore, temporal probability density functions (PDF) of tracer particles are studied in pore networks and an excellent agreement for the spectrum of transition times for particles from pore to pore is obtained between network model results and CTRW theory. Based on the truncated power-law interpretation of PDF-s, the physical origin of the power-law scaling of dispersion coefficient vs. Peclet number has been explained for unconsolidated porous media, sands and a number of sandstones, arriving at the same conclusion from numerical network modelling, analytic CTRW theory and experiment. The length traveled by solute plumes before Gaussian behaviour is reached increases with an increase in heterogeneity and/or Pe. This opens up the question on the nature of dispersion in natural systems where the heterogeneities at the larger scales will significantly increase the range of velocities in the reservoir, thus significantly delaying the asymptotic approach to Gaussian behaviour. As a consequence, the asymptotic behaviour might not be reached at the field scale. This is illustrated by the multi-scale approach in which transport at core, gridblock and field scale is viewed as a series of particle transitions between discrete nodes governed by probability distributions. At each scale of interest a distribution that represents transport physics (and the heterogeneity) is used as an input to model a subsequent reservoir scale. The extensions to reactive transport are discussed.

Pore-scale investigation on the response of heterotrophic respiration to moisture conditions in heterogeneous soils

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yan, Zhifeng; Liu, Chongxuan; Todd-Brown, Katherine E.

The relationship between microbial respiration rate and soil moisture content is an important property for understanding and predicting soil organic carbon degradation, CO 2 production and emission, and their subsequent effects on climate change. This paper reports a pore-scale modeling study to investigate the response of heterotrophic respiration to moisture conditions in soils and to evaluate various factors that affect this response. X-ray computed tomography was used to derive soil pore structures, which were then used for pore-scale model investigation. The pore-scale results were then averaged to calculate the effective respiration rates as a function of water content in soils.more » The calculated effective respiration rate first increases and then decreases with increasing soil water content, showing a maximum respiration rate at water saturation degree of 0.75 that is consistent with field and laboratory observations. The relationship between the respiration rate and moisture content is affected by various factors, including pore-scale organic carbon bioavailability, the rate of oxygen delivery, soil pore structure and physical heterogeneity, soil clay content, and microbial drought resistivity. Simulations also illustrates that a larger fraction of CO 2 produced from microbial respiration can be accumulated inside soil cores under higher saturation conditions, implying that CO 2 flux measured on the top of soil cores may underestimate or overestimate true soil respiration rates under dynamic moisture conditions. Overall, this study provides mechanistic insights into the soil respiration response to the change in moisture conditions, and reveals a complex relationship between heterotrophic microbial respiration rate and moisture content in soils that is affected by various hydrological, geochemical, and biophysical factors.« less

Tracer counterpermeation analysis of diffusivity in finite-length nanopores with and without single-file dynamics

DOE PAGES

Ackerman, David M.; Evans, James W.

2017-01-19

Here, we perform a tracer counterpermeation (TCP) analysis for a stochastic model of diffusive transport through a narrow linear pore where passing of species within the pore is inhibited or even excluded (single-file diffusion). TCP involves differently labeled but otherwise identical particles from two decoupled infinite reservoirs adsorbing into opposite ends of the pore, and desorbing from either end. In addition to transient behavior, we assess steady-state concentration profiles, spatial correlations, particle number fluctuations, and diffusion fluxes through the pore. From the profiles and fluxes, we determine a generalized tracer diffusion coefficient D tr(x), at various positions x within themore » pore. D tr(x) has a plateau value in the pore center scaling inversely with the pore length, but it is enhanced near the pore openings. The latter feature reflects the effect of fluctuations in adsorption and desorption, and it is also associated with a nontrivial scaling of the concentration profiles near the pore openings.« less

Tracer counterpermeation analysis of diffusivity in finite-length nanopores with and without single-file dynamics

NASA Astrophysics Data System (ADS)

Ackerman, David M.; Evans, James W.

2017-01-01

We perform a tracer counterpermeation (TCP) analysis for a stochastic model of diffusive transport through a narrow linear pore where passing of species within the pore is inhibited or even excluded (single-file diffusion). TCP involves differently labeled but otherwise identical particles from two decoupled infinite reservoirs adsorbing into opposite ends of the pore, and desorbing from either end. In addition to transient behavior, we assess steady-state concentration profiles, spatial correlations, particle number fluctuations, and diffusion fluxes through the pore. From the profiles and fluxes, we determine a generalized tracer diffusion coefficient Dtr(x ) , at various positions x within the pore. Dtr(x ) has a plateau value in the pore center scaling inversely with the pore length, but it is enhanced near the pore openings. The latter feature reflects the effect of fluctuations in adsorption and desorption, and it is also associated with a nontrivial scaling of the concentration profiles near the pore openings.

Multiscale pore structure and constitutive models of fine-grained rocks

NASA Astrophysics Data System (ADS)

Heath, J. E.; Dewers, T. A.; Shields, E. A.; Yoon, H.; Milliken, K. L.

2017-12-01

A foundational concept of continuum poromechanics is the representative elementary volume or REV: an amount of material large enough that pore- or grain-scale fluctuations in relevant properties are dissipated to a definable mean, but smaller than length scales of heterogeneity. We determine 2D-equivalent representative elementary areas (REAs) of pore areal fraction of three major types of mudrocks by applying multi-beam scanning electron microscopy (mSEM) to obtain terapixel image mosaics. Image analysis obtains pore areal fraction and pore size and shape as a function of progressively larger measurement areas. Using backscattering imaging and mSEM data, pores are identified by the components within which they occur, such as in organics or the clastic matrix. We correlate pore areal fraction with nano-indentation, micropillar compression, and axysimmetic testing at multiple length scales on a terrigenous-argillaceous mudrock sample. The combined data set is used to: investigate representative elementary volumes (and areas for the 2D images); determine if scale separation occurs; and determine if transport and mechanical properties at a given length scale can be statistically defined. Clear scale separation occurs between REAs and observable heterogeneity in two of the samples. A highly-laminated sample exhibits fine-scale heterogeneity and an overlapping in scales, in which case typical continuum assumptions on statistical variability may break down. Sandia National Laboratories is a multimission laboratory managed and operated by National Technology and Engineering Solutions of Sandia LLC, a wholly owned subsidiary of Honeywell International Inc. for the U.S. Department of Energy's National Nuclear Security Administration under contract DE-NA0003525.

DOE Office of Scientific and Technical Information (OSTI.GOV)

Yoon, Hongkyu; Major, Jonathan; Dewers, Thomas

Dissolved CO 2 in the subsurface resulting from geological CO 2 storage may react with minerals in fractured rocks, confined aquifers, or faults, resulting in mineral precipitation and dissolution. The overall rate of reaction can be affected by coupled processes including hydrodynamics, transport, and reactions at the (sub) pore-scale. In this work pore-scale modeling of coupled fluid flow, reactive transport, and heterogeneous reactions at the mineral surface is applied to account for permeability alterations caused by precipitation-induced pore-blocking. This paper is motivated by observations of CO 2 seeps from a natural CO 2 sequestration analog, Crystal Geyser, Utah. Observations alongmore » the surface exposure of the Little Grand Wash fault indicate the lateral migration of CO 2 seep sites (i.e., alteration zones) of 10–50 m width with spacing on the order of ~100 m over time. Sandstone permeability in alteration zones is reduced by 3–4 orders of magnitude by carbonate cementation compared to unaltered zones. One granular porous medium and one fracture network systems are used to conceptually represent permeable porous media and locations of conduits controlled by fault-segment intersections and/or topography, respectively. Simulation cases accounted for a range of reaction regimes characterized by the Damköhler (Da) and Peclet (Pe) numbers. Pore-scale simulation results demonstrate that combinations of transport (Pe), geochemical conditions (Da), solution chemistry, and pore and fracture configurations contributed to match key patterns observed in the field of how calcite precipitation alters flow paths by pore plugging. This comparison of simulation results with field observations reveals mechanistic explanations of the lateral migration and enhances our understanding of subsurface processes associated with the CO 2 injection. In addition, permeability and porosity relations are constructed from pore-scale simulations which account for a range of reaction regimes characterized by the Da and Pe numbers. Finally, the functional relationships obtained from pore-scale simulations can be used in a continuum scale model that may account for large-scale phenomena mimicking lateral migration of surface CO 2 seeps.« less

Effect of sub-pore scale morphology of biological deposits on porous media flow properties

NASA Astrophysics Data System (ADS)

Ghezzehei, T. A.

2012-12-01

Biological deposits often influence fluid flow by altering the pore space morphology and related hydrologic properties such as porosity, water retention characteristics, and permeability. In most coupled-processes models changes in porosity are inferred from biological process models using mass-balance. The corresponding evolution of permeability is estimated using (semi-) empirical porosity-permeability functions such as the Kozeny-Carman equation or power-law functions. These equations typically do not account for the heterogeneous spatial distribution and morphological irregularities of the deposits. As a result, predictions of permeability evolution are generally unsatisfactory. In this presentation, we demonstrate the significance of pore-scale deposit distribution on porosity-permeability relations using high resolution simulations of fluid flow through a single pore interspersed with deposits of varying morphologies. Based on these simulations, we present a modification to the Kozeny-Carman model that accounts for the shape of the deposits. Limited comparison with published experimental data suggests the plausibility of the proposed conceptual model.

Pore-scale modeling of saturated permeabilities in random sphere packings.

PubMed

Pan, C; Hilpert, M; Miller, C T

2001-12-01

We use two pore-scale approaches, lattice-Boltzmann (LB) and pore-network modeling, to simulate single-phase flow in simulated sphere packings that vary in porosity and sphere-size distribution. For both modeling approaches, we determine the size of the representative elementary volume with respect to the permeability. Permeabilities obtained by LB modeling agree well with Rumpf and Gupte's experiments in sphere packings for small Reynolds numbers. The LB simulations agree well with the empirical Ergun equation for intermediate but not for small Reynolds numbers. We suggest a modified form of Ergun's equation to describe both low and intermediate Reynolds number flows. The pore-network simulations agree well with predictions from the effective-medium approximation but underestimate the permeability due to the simplified representation of the porous media. Based on LB simulations in packings with log-normal sphere-size distributions, we suggest a permeability relation with respect to the porosity, as well as the mean and standard deviation of the sphere diameter.

Ultra-Long Time Dynamics of Contaminant Plume Mixing Induced by Transient Forcing Factors in Geologic Formations

NASA Astrophysics Data System (ADS)

Rajabi, F.; Battiato, I.

2016-12-01

Long term predictions of the impact of anthropogenic stressors on the environment is essential to reduce the risks associated with processes such as CO2 sequestration and nuclear waste storage in the subsurface. On the other hand, transient forcing factors (e.g. time-varying injection or pumping rate) with evolving heterogeneity of time scales spanning from days to years can influence transport phenomena at the pore scale. A comprehensive spatio-temporal prediction of reactive transport in porous media under time-dependent forcing factors for thousands of years requires the formulation of continuum scale models for time-averages. Yet, as every macroscopic model, time-averaged models can loose predictivity and accuracy when certain conditions are violated. This is true whenever lack of temporal and spatial scale separation occurs and it makes the continuum scale equation a poor assumption for the processes at the pore scale. In this work, we consider mass transport of a dissolved species undergoing a heterogeneous reaction and subject to time-varying boundary conditions in a periodic porous medium. By means of homogenization method and asymptotic expansion technique, we derive a macro-time continuum-scale equation as well as expressions for its effective properties. Our analysis demonstrates that the dynamics at the macro-scale is strongly influenced by the interplay between signal frequency at the boundary and transport processes at the pore level. In addition, we provide the conditions under which the space-time averaged equations accurately describe pore-scale processes. To validate our theoretical predictions, we consider a thin fracture with reacting walls and transient boundary conditions at the inlet. Our analysis shows a good agreement between numerical simulations and theoretical predictions. Furthermore, our numerical experiments show that mixing patterns of the contaminant plumes at the pore level strongly depend on the signal frequency.

Fluid displacement fronts in porous media: pore scale interfacial jumps, pressure bursts and acoustic emissions

NASA Astrophysics Data System (ADS)

Moebius, Franziska; Or, Dani

2014-05-01

The macroscopically smooth and regular motion of fluid fronts in porous media is composed of numerous rapid pore-scale interfacial jumps and pressure bursts that involve intense interfacial energy release in the form of acoustic emissions. The characteristics of these pore scale events affect residual phase entrapment and transport properties behind the front. We present experimental studies using acoustic emission technique (AE), rapid imaging, and liquid pressure measurements to characterize these processes during drainage and imbibition in simple porous media. Imbibition and drainage produce different AE signatures (AE amplitudes obey a power law). For rapid drainage, AE signals persist long after cessation of front motion reflecting fluid redistribution and interfacial relaxation. Imaging revealed that the velocity of interfacial jumps often exceeds front velocity by more than 50 fold and is highly inertial component (Re>1000). Pore invasion volumes reduced deduced from pressure fluctuations waiting times (for constant withdrawal rates) show remarkable agreement with geometrically-deduced pore volumes. Discrepancies between invaded volumes and geometrical pores increase with increasing capillary numbers due to constraints on evacuation opportunity times and simultaneous invasion events. A mechanistic model for interfacial motions in a pore-throat network was developed to investigate interfacial dynamics focusing on the role of inertia. Results suggest that while pore scale dynamics were sensitive to variations in pore geometry and boundary conditions, inertia exerted only a minor effect on phase entrapment. The study on pore scale invasion events paints a complex picture of rapid and inertial motions and provides new insights on mechanisms at displacement fronts that are essential for improved macroscopic description of multiphase flows in porous media.

Modeling of viscoelastic properties of nonpermeable porous rocks saturated with highly viscous fluid at seismic frequencies at the core scale

NASA Astrophysics Data System (ADS)

Wang, Zizhen; Schmitt, Douglas R.; Wang, Ruihe

2017-08-01

A core scale modeling method for viscoelastic properties of rocks saturated with viscous fluid at low frequencies is developed based on the stress-strain method. The elastic moduli dispersion of viscous fluid is described by the Maxwell's spring-dash pot model. Based on this modeling method, we numerically test the effects of frequency, fluid viscosity, porosity, pore size, and pore aspect ratio on the storage moduli and the stress-strain phase lag of saturated rocks. And we also compared the modeling results to the Hashin-Shtrikman bounds and the coherent potential approximation (CPA). The dynamic moduli calculated from the modeling are lower than the predictions of CPA, and both of these fall between the Hashin-Shtrikman bounds. The modeling results indicate that the frequency and the fluid viscosity have similar effects on the dynamic moduli dispersion of fully saturated rocks. We observed the Debye peak in the phase lag variation with the change of frequency and viscosity. The pore structure parameters, such as porosity, pore size, and aspect ratio affect the rock frame stiffness and result in different viscoelastic behaviors of the saturated rocks. The stress-strain phase lags are larger with smaller stiffness contrasts between the rock frame and the pore fluid. The viscoelastic properties of saturated rocks are more sensitive to aspect ratio compared to other pore structure parameters. The results suggest that significant seismic dispersion (at about 50-200 Hz) might be expected for both compressional and shear waves passing through rocks saturated with highly viscous fluids.