Search NASA⌕ Search

SEARCH · Search NASA

Results for “Heterogeneous media”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A gradient-based deep neural network model for simulating multiphase flow in porous media

We report simulation of multiphase flow in porous media is crucial for the effective management of subsurface energy and environment-related activities. The numerical simulators used for modeling such processes rely on spatial and temporal discretization of the governing mass and energy balance partial-differential equations (PDEs) into algebraic systems via finite-difference/volume/element methods. These simulators usually require dedicated software development and maintenance, and suffer low efficiency from a runtime and memory standpoint for problems with multi-scale heterogeneity, coupled-physics processes or fluids with complex phase behavior. Therefore, developing cost-effective, data-driven models can become a practical choice, and in this work, we choose deep learning approaches as they can handle high dimensional data and accurately predict state variables with strong nonlinearity. In this paper, we describe a gradient-based deep neural network (GDNN) constrained by the physics related to multiphase flow in porous media. We tackle the nonlinearity of flow in porous media induced by rock heterogeneity, fluid properties, and fluid-rock interactions by decomposing the nonlinear PDEs into a dictionary of elementary differential operators. We use a combination of operators to handle rock spatial heterogeneity and fluid flow by advection. Since the augmented differential operators are inherently related to the physics of fluid flow, we treat them as first principles prior knowledge to regularize the GDNN training. We use the example of pressure management at geologic CO 2 storage sites, where CO 2 is injected in saline aquifers and brine is produced, and apply GDNN to construct a predictive model that is trained with physics-based simulation data and emulates the physics process. We demonstrate that GDNN can effectively predict the nonlinear patterns of subsurface responses, including the temporal and spatial evolution of the pressure and saturation plumes. We also successfully extend the GDNN to convolutional neural network (CNN), namely gradient-based CNN (GCNN), and validate its capability to improve the prediction accuracy. GDNN has great potential to tackle challenging problems that are governed by highly nonlinear physics and enable the development of data-driven models with higher fidelity.

42 ENGINEERING↗
DOE OSTI2025

Knowledge Oriented Graph Unified Transformer (KOGUT) v0.1

KOGUT — Knowledge Oriented Graph Unified Transformer KOGUT implements the Relational Graph Transformer (RelGT) architecture for knowledge graph link prediction in biological domains, with a primary focus on microbial growth media prediction. While the original RelGT (arXiv:2505.10960) targets relational tables, time series, and multi-table databases, KOGUT adapts this architecture for heterogeneous biological knowledge graphs, providing first-in-class AI predictive models for microbial cultivation. Key Adaptations Beyond Original RelGT: - Knowledge Graph Focus: Applied to biological KGs with semantic node types (taxa, chemicals, media, phenotypes, environments) versus generic relational database tables, trained on the KG-Microbe knowledge graph (1.3M entities, 2.9M edges, 24 relation types). - Multimodal Node Encoding: Integrates node labels, categories, descriptions, and synonyms from KG metadata through learned embedding layers—adapting relational column features to graph node attributes with textual semantics. - Extended K-Hop Subgraph Strategy: Optimized neighborhood sampling (3-hop default, configurable up to 200 nodes) tuned for sparse biological networks, building on the original local-global attention framework with biological relation preservation. - Biolink Predicate Preservation: Type-specific transformations for 24 biological edge semantics (occurs_in, consumes, produces, has_phenotype, subclass_of) beyond standard relational foreign keys, enabling multi-relation link prediction. - Inductive Learning Support: Enables zero-shot predictions for novel taxa through feature-based embeddings (temperature, oxygen requirements, gram stain, cell shape), extending the original transductive relational benchmark scope to uncultured microorganisms. CheapSOTA Performance Optimizations (This Distribution): - VQ-EMA Centroid Attention: Vector quantization with exponential moving average for improved global context modeling (+5-10% MRR improvement). - HDF5 Precomputed Data Loading: One-time preprocessing of k-hop subgraphs to eliminate redundant graph traversals (2-5× training speedup). - Distributed Data Parallel Training: Multi-GPU support for scaling to larger knowledge graphs (tested on 4× NVIDIA A100 GPUs at NERSC Perlmutter). - Mixed Precision Training: Automatic mixed precision (AMP) for memory efficiency and faster training. Advantages Over Standard Knowledge Graph Embedding Models: Combines RelGT's proven multi-element tokenization (features, type, hop, structure) with graph-native biological representations, enabling interpretable link prediction across heterogeneous entities that standard embedding models (TransE, RotatE, ComplEx) and table-based transformers cannot directly model. Achieves near-perfect performance on microbial growth media prediction (MRR: 0.9966, Precision@1: 0.9932, Hit@10: 1.0000) while maintaining explainability through attention-based reasoning over biological pathways. Training Data: - KG-Microbe merged knowledge graph: 1,379,337 nodes, 2,960,472 edges - 24 biological relation types including taxonomic hierarchies, metabolic interactions, phenotype associations, and environmental relationships - Primary prediction task: Growth media suitability for microbial taxa (biolink:occurs_in, 50K edges) - Multi-relation capability: Predicts links for any of the 24 relation types, including chemical consumption/production, phenotype associations, and taxonomic classification Citation: Original RelGT Architecture: Dwivedi et al., "Relational Graph Transformer", arXiv:2505.10960, 2025 KOGUT Implementation: Knowledge Oriented Graph Unified Transformer for Microbial Growth Media Prediction Developed at Lawrence Berkeley National Laboratory (LBNL) Trained on NERSC Perlmutter supercomputer

Joachimiak, Marcin [Lawrence Berkeley National Lab↗
DOE OSTI2021

The effects of in-situ emulsion formation and superficial velocity on foam performance in high-permeability porous media

Foam has been considered as a mobility control agent for Enhanced Oil Recovery (EOR) applications in heterogeneous and fractured reservoirs. Foam enhances sweep efficiency by enabling the injected gas to flow from high-permeability to low-permeability regions. The foam injection velocity is an essential factor affecting foam mobility control and flow behavior. Furthermore, the presence of remaining oil adversely impacts foam generation leading to a failure in governing mobility control. In this study, foam performance in establishing mobility control effect and favorable apparent viscosity was experimentally investigated under a range of total superficial velocities in the presence and absence of crude oil during foam propagation in high permeability sand packs representing a propped fracture network. In the absence of oil, a pronounced effect of total superficial velocity on the mobility control performance was observed at the transition foam quality regime. Co-injection of oil and foam led to noticeable foam destabilization, as expected. However, a further increase in oil fractional flow resulted in larger overall flow resistance. We found that the formation of in-situ surfactant/oil emulsion contributed to the high apparent viscosity exhibited at the steady-state condition. These findings suggest that the steady-state flow resistance established during foam flooding at high oil saturations should be interpreted as the mobility control effect resulting from emulsion formation rather than the development of strong foam.

02 PETROLEUM↗
DOE OSTI2022

A Stochastic Transport Model for the Cumulative Number of Fissions and Deposited Fission Energy

The stochastic theory of neutron transport is extended to describe the cumulative distribution of fission numbers and deposited fission energy in a multiplying assembly. Solutions for the probability distributions are obtained using analytical approximations and Monte Carlo simulation in lumped geometry and in symmetric homogeneous and heterogeneous spheres. The results show the development of a power-law tail in the steady state fission number and deposited energy distributions when the medium is critical, independent of the fission neutron multiplicity distribution and domain heterogeneity. In contrast, the asymptotic decay is faster than exponential in subcritical media due to rapid chain extinction and in supercritical media due to the increasing probability of chain divergence. Here, a formal asymptotic analysis of the problem in lumped geometry with an arbitrary fission neutron multiplicity confirms the existence of power-law tails at critical.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Nanoparticle Tracking to Probe Transport in Porous Media

From the granular and fractured subsurface environment to highly engineered polymer membranes used in pharmaceutical purification, porous materials are ubiquitous in nature and industrial applications. In particular, porous media are used extensively in processes including water treatment, pharmaceutical sterilization, food/beverage processing and heterogeneous catalysis, where hindered mass transport is either essential to the process or a necessary, but undesirable limitation. Unfortunately, there are currently no universal models capable of predicting mass transport based on a description of the porous material because real porous materials are complex, and because many coupled dynamic mechanisms (e.g. adsorption, steric effects, hydrodynamic effects, electrostatic interactions, etc.) give rise to the observed macroscopic transport phenomena. While classical techniques, like nuclear magnetic resonance and dynamic light scattering, provide useful information about mass transport in porous media at the ensemble level, they provide limited insight into the microscopic mechanisms that give rise to complex phenomena such as anomalous diffusion, hindered pore-space accessibility, and unexpected retention under flow, among many others. To address this issue, we have developed refractive index matching imaging systems, combined with single-particle tracking methods, allowing the direct visualization of single particle motion within a variety of porous materials. Here, we summarize our recent efforts to advance the understanding of nanoparticle transport in porous media using single-particle tracking methods in both fundamental and applied scenarios. First, we describe the basic principles for two-dimensional and three-dimensional single-particle tracking in porous materials. Then, we describe concrete examples of nanoparticle transport in porous materials from two perspectives: (1) understanding fundamental elementary particle transport processes in porous media, including pore accessibility and cavity escape, which limit transport in porous media; (2) facilitating applications in industrial processes, e.g. by understanding the mechanisms of particle fouling and remobilization in filtration membranes. Finally, we provide an outlook of opportunities associated with investigating other types of mass transport in confined environments using single-particle tracking methods, including electrophoretic and self-propelled motion.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1986

Finite difference simulations of seismic scattering Implications for the propagation of short-period seismic waves in the crust and models of crustal heterogeneity

Synthetic seismographs that were obtained by the finite difference method are presently applied to the study of elastic and acoustic wave scattering in two-dimensional media with random spatial variations in seismic velocity. The seismograms are analyzed to determine the variation in travel times and waveforms across arrays of receivers. The random media with Gaussian and exponential correlation functions considered differ in the spectral falloff of their velocity fluctuations at wavelengths smaller than 2pi times the correlation distance. It is found that alternative models of crustal heterogeneity can be tested by improved measurements of the frequency dependence of the crustal Q at frequencies greater than about 1 Hz, assuming that scattering is responsible for most of the attenuation at such frequencies.

Frankel, A.↗
DOE OSTI2020

Use of an Alternative Conceptual Model of Vadose Zone Heterogeneity to Evaluate Past Tank Leaks and Other Unplanned Releases within a Tank Farm at the Hanford Site - 20061

Washington State Department of Ecology (Ecology) requested that U.S. Department of Energy (DOE)-Office of River Protection (ORP) consider an evaluation of effects of fine-grained thin sediment layers on transport with a separate alternative conceptual model in its evaluation of the potential impact of tank leaks and other unplanned releases within the tank farm at the Hanford Site. For this alternative model, Ecology recommended that the model be developed based on the general framework of fine-grained units identified by a stakeholder group in their interpretations of variability in moisture content data collected in the vicinity of Waste Management Area (WMA) C. In discussions with DOE-ORP, Ecology acknowledged that the underlying data and interpretations of the occurrence and lateral continuity of the fine-grained thin layers identified by a stakeholder group are uncertain. However, DOE-ORP agreed to Ecology's recommendation and has provided support for the requested evaluation that involved development of an alternative model based on the general framework of the unpublished report by the stakeholder. DOE-ORP considered this evaluation to be a hypothetical evaluation of vadose zone heterogeneities at WMA C. General observations from the range of simulation cases examined in the evaluation of the effects of hypothetical vadose zone heterogeneities at WMA C are as follows. - The movement of the center of mass of the simulated plumes was generally vertically downward below the source for all simulations, including those that incorporated the hypothetical fine-grained units. - All simulations that incorporated hypothetical heterogeneity produced additional plume spreading over what was produced in simulations using Equivalent Homogeneous Media (EHM) model(s)a. The spreading resulted in a broadening of the fringes of the plume, resulting in a wider region of low concentration, but lower peak concentrations associated with the center of mass of the plume. - Simulations that used the silty-sand hydraulic properties, suggested by Ecology for the hypothetical fine-grained units, generally produced similar spreading and slightly lower peak mass flux at the water table when compared to the EHM modeling results. a The EHM-based models used at WMA C do not explicitly include small-scale fine-grained heterogeneities used in the model advocated by Ecology. - Simulations that used the silty-sand hydraulic properties, suggested by Ecology for the hypothetical fine-grained units, generally produced less spreading and an earlier arrival of mass flux at the water table when compared to the use of another set of hydraulic properties from a silty-sand sample collected at a nearby disposal facility. - The EHM representation of the vadose zone generally produced higher peak mass flux and an earlier occurrence of peak fluxes at the water table compared to all analyses incorporating additional hypothetical heterogeneity. Results of this alternative model evaluation of past leaks provided some insight into the transport effects of vadose zone heterogeneities at WMA C that helped resolve Ecology's comments and issues related to the effects of vadose zone heterogeneities on past leaks and losses from the WMA C tank farms area. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Peridynamic Models for Random Media Found by Coarse Graining

Using coarse graining, the upscaled mechanical properties of a solid with small scale heterogeneities are derived. The method maps internal forces at the small scale onto peridynamic bond forces in the coarse grained mesh. These upscaled bond forces are used to calibrate a peridynamic material model with position-dependent parameters. These parameters incorporate mesoscale variations in the statistics of the small scale system. The upscaled peridynamic model can have a much coarser discretization than the original small scale model, allowing larger scale simulations to be performed efficiently. The convergence properties of the method are investigated for representative random microstructures. In conclusion, a bond breakage criterion for the upscaled peridynamic material model is also demonstrated.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stochastic Inversion of Gaussian Random Media Using Transverse Coherence Functions for Reflected Waves: Theory and Method

The transverse coherence functions (TCFs) of phase and amplitude fluctuations of a seismic wave are powerful to estimate the spatial distribution, length scales, and strength of random heterogeneities. However, TCFs have been formulated for transmitted waves only, not for reflected waves. In this paper, we derive reflection TCFs for Gaussian random media. Furthermore, we propose to invert for Gaussian random media using the reflection TCFs based on the grid search. We validate the new reflection TCF formulas using 2D finite-difference numerical experiments. The numerical example also illustrates the feasibility and efficiency of the inversion. The stochastic inversion using reflected waves can be used in both exploration and global seismology.

58 GEOSCIENCES↗
NASA NTRS1998

Spiral waves are stable in discrete element models of two-dimensional homogeneous excitable media

The spontaneous breakup of a single spiral wave of excitation into a turbulent wave pattern has been observed in both discrete element models and continuous reaction-diffusion models of spatially homogeneous 2D excitable media. These results have attracted considerable interest, since spiral breakup is thought to be an important mechanism of transition from the heart rhythm disturbance ventricular tachycardia to the fatal arrhythmia ventricular fibrillation. It is not known whether this process can occur in the absence of disease-induced spatial heterogeneity of the electrical properties of the ventricular tissue. Candidate mechanisms for spiral breakup in uniform 2D media have emerged, but the physical validity of the mechanisms and their applicability to myocardium require further scrutiny. In this letter, we examine the computer simulation results obtained in two discrete element models and show that the instability of each spiral is an artifact resulting from an unphysical dependence of wave speed on wave front curvature in the medium. We conclude that spiral breakup does not occur in these two models at the specified parameter values and that great care must be exercised in the representation of a continuous excitable medium via discrete elements.

NASA Discipline Cardiopulmonary↗
DOE OSTI2024

Effect of Nucleation Heterogeneity on Mineral Precipitation in Confined Environments

Abstract The formation of new mineral phases in confined environments, especially in porous media, is crucial for various geological processes like mineralization and diagenesis. The nucleation and precipitation of minerals are initiated at the microscale through fluid‐rock interaction, where dissolution of primary phases leads to supersaturated conditions and nucleation and growth of secondary ones. Previous research has focused primarily on either precipitation or nucleation, without fully exploring their combined impact. Our study introduces a computational framework that integrates classical nucleation theory with the micro‐continuum method. We validated our model by comparing with experiments, and discovered that different surface nucleation rate changes the mode of precipitation from a preferential to uniform precipitate textures. Furthermore, our study uncovered that the conventional deterministic precipitation method tends to underestimate the permeability of the porous matrix. In contrast, the new framework significantly improves model accuracy by incorporating preferential precipitation and heterogeneous nucleation.

58 GEOSCIENCES↗
DOE OSTI2021

Striving to translate shale physics across ten orders of magnitude: What have we learned?

Shales will play an important role in the successful transition of energy from fossil-based resources to renewables in the coming decades. Aside from being a significant source of low-carbon intensity fuels, like natural gas, they also serve as geologic seals of subsurface formations that may be used to isolate nuclear waste, sequester CO 2 , or store intermittent energy (e.g., solar hydrogen). Despite their importance, shales pose significant engineering and environmental challenges due to their nanoporous structure and extreme heterogeneity that spans at least ~10 orders of magnitude in spatial scale. Two challenges inhibit a system-level understanding: (1) the physics of fluid flow and phase behavior in shales are poorly understood due to the dominant molecular interactions between minerals and fluids under confinement, and (2) the apparent lack of scale separation that prevents a reliable (closed) description of the physics at any single scale of observation. In this review, we focus on the latter issue and discuss scale translation, which in its broadest sense is transforming data or simulations from one spatiotemporal scale to another. While effective scale translation is not exclusive to shales, but all geologic porous media, the need for it is especially acute in shales given their high degree of heterogeneity. Classical theories like homogenization, while indispensable, fail when scales are not separated. Other methods, like numerical upscaling, scale-translate in only one direction: small to large, but not the reverse, called downscaling. However, the confluence of advances in three areas are bringing challenging problems such as shales within reach: increased computational power and scalable algorithms; high-resolution imaging and multi-modal data acquisition; and machine learning to process massive amounts of data. While these advances equip geoscientists with a wide array of experimental and computational tools, no individual tool can probe the entire gamut of heterogeneity in shales. Their effective use, therefore, requires an ability to bridge between various data types obtained at different scales. The aim of this review is to present a coherent account of computational and experimental methods that may be used to achieve just that, i.e., to perform scale translation. We provide a broader definition of scale translation, one that transcends classical homogenization and upscaling methods, but is consistent with them and accommodates notions like downscaling and data translation. After a brief introduction to homogenization, we review hybrid methods, numerical upscaling and its recent extensions, multiscale computing, high-resolution imaging, and machine learning. We place particular emphasis on multiscale computing and propose an algorithmic framework to bridge between the pore (micro) and Darcy (macro) scales. Throughout the paper, we draw comparisons between the various methods and highlight their (often hidden) similarities, differences, benefits, and pitfalls. We finally conclude with two case studies on shales that exemplify some of the methods presented.

58 GEOSCIENCES↗
DOE OSTI2024

Effects of temperature inversion on densification in chemical vapor infiltration

In a classical chemical vapor infiltration (CVI) process, the competing effects of chemical kinetics and reagent gas transport lead to non–uniform depositions such that outer layers of a preform densify faster leaving the core highly porous. Currently, CVI must be performed at a sufficiently low temperature to achieve good densification quality which leads to high processing time and cost. Volumetric heating of the preform, especially through microwaves, can create temperature inversion such that the core is hotter than the outer surface and potentially, overcome the challenges associated with isothermal CVI. Direct numerical simulations (DNS) of densification under various such temperature distributions indicate that microwave heating in CVI processing can lead to better (uniform) densification of porous preforms. Here the role of key parameters describing the temperature distributions on the densification behavior is investigated. Strategic temporal control of the temperature distribution shows that processing times can be reduced by almost half while maintaining a good densification quality similar to that of low–temperature isothermal processing. Inside–out densification due to the inverted temperature profile is a key distinguishing characteristic of microwave assisted CVI.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Subsurface Spectroscopy in Heterogeneous Materials Using Self-Healing Laser Beams

Self-healing optical beams are a class of propagation modes that can recover their beam shapes after distortion or partial blockage. This self-healing property makes them attractive for use in applications involving turbid media as they can—in theory—penetrate further into these materials than standard Gaussian beams. In this paper, we characterize the propagation of two different self-healing beams (Bessel and Airy) through a solid scattering material with different scatterer concentrations and find that both beams do recover after scattering for samples below a threshold scatterer concentration. Additionally, we test the applicability of both beam shapes for improved sub-surface spectroscopy in heterogeneous materials using fluorescent particles and find that there is an average fluorescence intensity enhancement of 1.3× using self-healing beams versus a standard Gaussian beam.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

Scaling Analysis of Two–Phase Flow in Fractal Permeability Fields

Fluid mixing in permeable media is essential in many practical applications. The mixing process is a consequence of velocity fluctuations owing to geological heterogeneities and mobility contrast of fluids. Heterogeneities in natural rocks are often spatially correlated, and their properties, such as permeability, may be described using fractal distributions. This work models the fractal characteristics of such permeability fields in which the covariance function is expressed as a power-law function. A generalized scaling relation is derived relating various fractal permeability fields using the magnitude of their fluctuations. Here, this relation reveals the self-similar behavior of two-phase flow in such permeable media. To that end, a recently developed, high-resolution numerical simulator is employed to validate the analytically derived scaling relations. Two flow problems are considered in which flow is governed by 1) a linear, and 2) a nonlinear transport equation. Due to the probabilistic representation of the fractal permeability fields, a sensitivity study is conducted for each flow scenario to determine the number of realizations required for statistical convergence. Scaling analysis is performed using ensemble averages of simulated saturation profiles and their mixing lengths. Results support the validity of the developed scaling relation across the range of investigated flow conditions at intermediate times. The dynamics of linear flow in the asymptotic regime is affected by the correlation structure of heterogeneity. In nonlinear flow, scaling behavior appears to be dominated by the degree of nonlinearity.

58 GEOSCIENCES↗
DOE OSTI2020

Markovian transport processes in a heterogeneous, variably saturated watershed: A multi-domain spatial Markov model

Transport processes in watersheds remain challenging to describe accurately with upscaled models because of their intricate complexities at multiple scales that can lead to multi-modal, non-Fickian breakthrough curves. One of the recent advances in solute transport modeling has been the Lagrangian spatial Markov model (SMM), which describes transport using velocity increments and models the correlated transitions between them. The method has been applied previously to saturated porous media and fractured media, but these successful applications suggest it may be useful for describing watershed scale processes as well. An existing, 3-D, heterogeneous variably saturated model of a headwaters catchment and classical random walk particle tracking were used to generate a forward simulation and the resulting transport was analyzed with a spatial Markov framework. Existing SMM methods were used to analyze the transport paths through the complex, unconfined domain and we found that the vadose zone plays an important role in the correlation structure. As such, a new kind of SMM is proposed, termed a multi-domain SMM (MD-SMM), that uses separate correlation models for auto-transitions in the vadose zone and the saturated zone, and for the cross-transitions between them. The MD-SMM improves representation of transport relative to the single domain model and can naturally delineate times when a particle is in the vadose zone or the saturated zone, adding an additional degree of realism to the upscaled model. Lastly, we show that the transition probability matrix for transport in the saturated portion of this unconfined system is approximately Markovian, which is an important validation of the proposed framework.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Flow Channeling in Fracture Networks: Characterizing the Effect of Density on Preferential Flow Path Formation

Flow channelization is a commonly observed phenomenon in fractured subsurface media where the flow of fluids is restricted primarily to highly transmissive fracture networks surrounded by a low-permeability rock matrix. The multiscale structural heterogeneity of these networks results in multiscale flow channelization where preferential flow paths form at length scales ranging from the entire system down to the subfracture size. We present an analysis of how one of the largest scales in fractured media, the network density, influences the degree of flow channeling that occurs using an ensemble of semigeneric three-dimensional discrete fracture network (DFN) simulations. We construct 10 DFNs, whose fracture lengths follow a power law distribution, at four densities for a total of 40 networks. We characterize their structure in terms of the network topology and geometry. Eulerian and Lagrangian observations of the steady-state flow fields obtained within the networks are used to quantify the degree of flow channelization at the network scale. We introduce a measure for the importance-ranking/hierarchy of different flow paths in the network using graph-based analysis of Lagrangian transport by which the degree of flow channeling between networks is compared. These flow observations are then linked to the structural properties of the networks. In general, network-scale flow channeling decreases as the network density increases. However, at low densities, there is more uniform flow within the entire connected network than in high-density networks. We also demonstrate how standard transport observables can be used to infer the degree of flow channelization occurring within a fracture network.

58 GEOSCIENCES↗
DOE OSTI2021

Impacts of Pore Network-Scale Wettability Heterogeneity on Immiscible Fluid Displacement: A Micromodel Study

The mixed-wet nature of reservoir formations imposes a wide range of rock wettability from strong resident-fluid wetting to strong invading-fluid wetting. The characteristics of two-phase flow in porous media composed of mixed-wetting surfaces remain poorly understood. In this study, we investigated the displacement of resident ethylene glycol (EG) by hexane in two mixed-wet micromodels of identical 2.5-D geometry heterogeneity, with uniformly or heterogeneously distributed patches strongly wetting to hexane. These patches are mixed among pores with unaltered EG-wetting surfaces. Along with control tests in the originally EG-wet micromodel, here we show the classic fingering and transitions in flow regimes at $logCa$ (capillary number) from -7.2 to -3.9. Moreover, pore-scale distributions of wettability and their spatial correlation influence displacement efficiency. In the two mixed-wet micromodels, we found (a) an increase of steady-state hexane saturation at the end of experiments by up to 0.12 in the capillary fingering regime and a decrease of at most by 0.06 in the viscous fingering regime, compared to the EG-wet micromodel, and (b) dispersed and fragmented hexane distribution after displacement. Brine drainage during supercritical CO 2 (scCO 2 ) injections in these micromodels occurs with lower wettability contrasts, and under similar viscosity ratios and interfacial tensions resulted in higher displacement efficiency relative to displacement of EG by hexane. While mixed-wettability can enhance displacement efficiency compared to uniform wettability, the dynamics of immiscible fluids in strong mixed-wet reservoirs are expected to be less pronounced in contributing to the efficiency of geological CO 2 sequestration, oil recovery, and remediation of hydrocarbon-contaminated aquifers.

2.5-D micromodel↗