Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “permeability models”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Numerical Simulations of Carbon Dioxide Storage Efficiency in Heterogeneous Reservoir Models

The U.S. Department of Energy’s National Energy Technology Laboratory (DOE-NETL) has been developing methods and tools (the online Carbon Dioxide Storage prospeCtive Resource Estimation Excel aNalysis (CO2-SCREEN) tool) to estimate carbon dioxide (CO2) storage potential in subsurface reservoirs. The CO2 storage efficiency terms are input in the tool to calculate storage potential in targeted reservoirs. In this effort, two CO2 storage efficiency terms were evaluated: volumetric displacement ( E V ) and microscopic displacement ( E d ). The first term deals with efficiency of CO2 propagation into an accessible reservoir volume, while the second term evaluates effectiveness of native fluid displacement with CO2. The interpreted well logs and core sample measurements were applied to create the heterogeneous reservoir models including geostatistical realizations of porosity and intrinsic permeability fields. Supercritical CO2 was injected over the course of 30 years into brine-saturated reservoir models for clastics, limestone, and dolomite lithologies and deltaic fluvial, aeolian, shallow marine, and reef depositional environments by means of varying reservoir parameters and injection scenarios. The reservoir models providing vertically heterogeneous petrophysical properties and designated as “layered reservoir models” (with homogeneous parameters along each layer of the model) were not determined to be a transition between the homogeneous and heterogeneous models in respect to storage efficiency. Another finding shows that high-efficiency factors do not necessarily mean increased CO2 storage; they rather indicate that the available volume and pore space are more fully utilized. The CO2 storage efficiency factors were evaluated dynamically at the select time points using P 10 ‐ P 50 ‐ P 90 percentiles. The results of this study show that the P 10 ‐ P 90 distribution for volumetric efficiency is wider when compared to the microscopic efficiency. It was found that where dominant buoyancy forces drive the plume to the top of a target formation, the volumetric efficiency is low. Tighter sandstone and carbonate formations show prevalence of capillary forces and better utilization of reservoir volume.

Myshakin, Evgeniy M.↗
DOE OSTI2023

Advective gas flow in bentonite: Development and comparison of enhanced multi-phase numerical approaches

Understanding the impact of repository gas, generated from degradation of waste and its interaction with the host rock, is essential when assessing the performance and safety function of long-term disposal systems for radioactive waste. Numerical models based on conventional multi-phase flow theory have historically been applied to predict the outcome and impact of gas flow on different repository components. However, they remain unable to describe the full complexity of the physical processes observed in water-saturated experiments (e.g., creation of dilatant pathways) and thus, the development of novel representations for their description is required when assessing fully saturated clay-based systems. This was the primary focus of Task A within the international cooperative project DECOVALEX-2019 (D-2019) and refinement of these approaches is the primary focus of this study (Task B in the current phase of DECOVALEX-2023). This paper summarises development of enhanced numerical representations of key processes and compares the performance of each model against high-quality laboratory test data. Experimental data reveals that gas percolation in water-saturated compacted bentonite is characterised by four key features: (i) a quiescence phase, followed by (ii) the gas breakthrough, which leads to a (iii) peak value, which is then followed by (iv) a negative decay. Three models based on the multiphase flow theory have been developed. These models can provide good initial values and reasonable responses for gas breakthrough (although some of them still predict a too-smooth response). Peak gas pressure values are in general reasonably well captured, although maximum radial stress differences are observed at 48 mm from the base of the sample. Here, numerical peak values of 12.8 MPa are predicted, whereas experimental values are about 11 MPa. These models are also capable of providing a reasonable representation of the negative pressure decay following peak pressure. However, other key specific features (such as the timing of gas breakthrough) still require a better representation. The model simulations and their comparison with experimental data show that these models need to be further improved with respect to model parameter calibration, the numerical representation of spatial heterogeneities in material properties and flow localisation, and the upscaling of the related physical processes and parameters. To further understand gas flow localisation, a new conceptual model has been developed, which shows that discrete channels can possibly be induced through the instability of gas-bentonite interface during gas injection, thus providing a new perspective for modeling gas percolation in low-permeability deformable media.

58 GEOSCIENCES↗
DOE OSTI2025

A comparative analysis of numerical approaches for the description of gas flow in clay-based repository systems: From a laboratory to a large-scale gas injection test

There is nowadays a consensus among many countries that geological disposal is a favourable solution for the long-term management. Although different host formations and different barrier systems are under consideration around the world, clay-based materials form an important component for waste isolation in most national programmes. Hence, a good comprehension of the effect of gas flow on the hydro-mechanical behaviour of clay-based soils is essential, both at laboratory and field scale. Task B under the international cooperative project DECOVALEX-2023 has recently shown that, after some enhancement, models can be employed to reproduce laboratory scale tests, even with different sample geometries 37 . However, further work is required to understand whether they can be applied to simulate a large-scale experiment. Up-scaling of models for the advective transport of gas through clay-based low permeable material presents a number of problems related to the difficulty in obtaining consistent hydrogeological parameters and constitutive relationships at both laboratory and field scale. Based on a unique dataset from a large-scale gas injection test (Lasgit) performed at the Äspö Hard Rock Laboratory (Sweden), Task B within DECOVALEX-2023 has explored the refinement of these numerical strategies applied to the simulation of gas flow. Work performed within the task reveals that codes do not need to be substantially modified from the laboratory models to reproduce full-scale tests: indeed, model parameters calibrated and validated at laboratory scale have been applied to predict field scale gas flow at Lasgit, including peak gas pressure and injected cumulative gas volume. By means of (1) the introduction of interfaces between blocks to reflect the experimental configuration and the (2) adjustment of some parameters (e.g., higher permeability), the updated models are able to represent most of the key features observed in the experimental data, even at a large scale.

Tamayo-Mas, E↗
DOE OSTI2020

Proppant embedment in coal and shale: Impacts of stress hardening and sorption

During methane production in CBM reservoirs, the influence of proppant embedment and permeability damage cannot be neglected – especially where the wall-rock is soft. Effective stresses are elevated during methane recovery, increasing both normal loading stress and confinement and simultaneously overprinting sorption-induced volumetric strains. Experiments and analytic modeling are conducted to define key mechanisms controlling these competitive effects. We independently measure overall sample compaction (external LVDT) and local strain (strain gauge) in the matrix to deconvolve proppant embedment in a propped fracture for different conditions of confining stress. The results show symptomatic behaviors of elastic (shale) and elastoplastic (coal) responses of embedment. Different from shale, the evolution of embedment is convex upwards with increased stress where indented depth increases more rapidly as loading stress increases under constant confinement. In addition, a stress-hardening effect is found to play a pivotal role in determining the characteristics of indentation, which are examined in terms of evolution profiles, deformation regimes, embedment slopes, curvatures, yield points and irreversible indentations. Based on the experimental observations a semianalytical model predicts indentation and the evolution of propped permeability under recreated in-situ stress conditions. A simplified case study is conducted to further illustrate the evolution of aperture and permeability of a propped fracture in CBM reservoirs. The modeling results suggest that proppant embedment is significantly overestimated if the variable stress-hardening (VSH) effect is neglected, especially when effective stress is large. Moreover, a decrease in indentation depth possibly occurs during late stage methane production, resulting in a reversal/recovery in fracture closure. This is because desorption-induced shrinkage becomes the predominant effect, causing an increase in aperture and a reduction in the indented volume of proppant. The resulting recovery in permeability implies that the propped coal fracture has the potential to optimally facilitate methane production as a pathway, even at high closure stresses generated by methane drainage.

01 COAL, LIGNITE, AND PEAT↗
DOE OSTI2020

Improved pore network models to simulate single-phase flow in porous media by coupling with lattice Boltzmann method

In this paper, different pore network models to simulate single-phase flow in porous media are built and their accuracy is evaluated. In addition to the conventional pore network model (CPNM) which consists of regular pore bodies and throat bonds, three improved pore network models (IPNMs) are developed allowing to better describing the real pore and throat geometry. The first improved pore network model (IPNM1) replaces the regular throat bond with a throat bond showing the real throat cross section. The second improvement (IPNM2) uses a series of sub-throat bonds with varying cross sections to better describe the real throat geometry, which is firstly proposed in this paper. The third model (IPNM3) extracts the real pore-throat-pore geometry without simplification. The conductance of fluid flow through these more realistic throat bonds is calculated by the lattice Boltzmann method (LBM). The accuracy and computational efficiency of the different pore network models are evaluated taking the LBM simulation over the whole porous medium as reference solution. The global permeability and detailed pressure distributions in the pores for the different pore network models are validated. The results show that the accuracy of the pore network model increases from CPNM to IPNM3, but at the expense of increasing computational cost. This study suggests that IPNM3 can replace a whole-domain LBM simulation with similar accuracy but much lower computational cost. As a first-order approximation the newly proposed IPNM2 is suggested as good compromise between accuracy and computational cost.

42 ENGINEERING↗
DOE OSTI2022

Geomechanical and Hydrogeological Evaluation of a Shallow Hydraulic Fracture at the Devine Fracture Pilot Site, Medina County, Texas

UT-Austin’s Devine Fracture Pilot Site, 50 miles southwest of San Antonio, Texas, has been targeted for a comprehensive, multidisciplinary development of fracture diagnostic techniques that are cross-validated by ground-truth data acquisition near a recently created, 175-ft-deep, horizontal hydraulic fracture (Ahmadian et al. 2018 Demonstration of proof of concept of electromagnetic geophysical methods for high resolution illumination of induced fracture networks. In Proceedings of the SPE Hydraulic Fracturing Technology Conference and Exhibition, The Woodlands, Texas, USA, 23–25 January 2018. SPE-189858-MS.). To evaluate the fracture diagnostic methods at this site, we conducted injection tests with a predefined volumetric flow-rate profile, resembling a diagnostic fracture injection test on September 2020. Subsequently, we developed hydrogeological and geomechanical models based on flow-rate and bottomhole-pressure measurements. History-matching efforts using a simplified layer-cake hydrogeological model resulted in the field-scale formation permeability of 9.87 × 10 –15 m 2 (10 mD) and Darcy-scale fracture permeability. The analysis of the bottomhole pressure and injection-rate history showed that (1) the newly created horizontal fracture was closed adjacent to the injection well pre-injection and (2) the initial pump-pressure increase at a nominal volumetric injection rate led to near-well fracture reopening, fluid conductivity increase, and abrupt injection-rate increase. To overcome hydrogeological-model limitations of predicting fracture reopening throughout injection, we extended the modeling to a finite-element, poroelastic analysis of horizontal-fracture growth using a cohesive-zone model. Using this fracture-reopening model, we improved the history match of the transient-pressure response during the experiment by adjusting the hydromechanical properties. Furthermore, post-injection pressure transient analyses helped reduce uncertainty in the overburden-stress gradient, and the initial hydraulic-fracturing simulation verified the plausibility of achieving the surveyed propped fracture area.

02 PETROLEUM↗
DOE OSTI2026

Determining the extent of potential fugitive fluid migration from geologic carbon storage in hydrocarbon-bearing reservoirs: Insights from one-dimensional numerical modeling

Numerical modeling of Geologic Carbon Sequestration in permeable reservoirs initially containing hydrocarbons is conducted using the multi-phase, multi-component thermohydrologic simulator TOGA (TOUGH Oil, Gas, Aqueous; TOUGH stands for Transport Of Unsaturated Groundwater and Heat), to determine how phase and composition of the original fluids influence the extent of the zone where upward fugitive fluid migration could potentially occur, denoted R f . The area within R f comprises regions of substantially elevated pressure and free-phase CO 2 saturation, where a breach in reservoir sealing capacity would lead to upward fugitive fluid migration. The model examines the conditions within the storage reservoir that could lead to fugitive flow, but does not model the fugitive flow itself. A one-dimensional radial model of the storage reservoir is used, and three initial phase conditions are considered: single-phase aqueous, two-phase gas-aqueous, and three-phase oil-gas-aqueous. Components that may be present are H 2 O, CO 2 , CH 4 , C 4 H 10 , and C 10 H 22 . The most important factors controlling Rf are (1) the initial gas-phase saturation within the reservoir, and (2) the lateral extent of multi-phase initial conditions, particularly CO 2 . The composition of liquid and gas phases has a secondary effect. The impact of reservoir depth, thickness, injection rate, and hydrologic properties are also briefly examined, with thickness (or equivalently injection rate) having the biggest effect. These results can help to understand important trends in potential response of CO 2 -EOR fields being considered for dedicated CO 2 storage.

CO₂ plume migration↗
DOE OSTI2023

Effects of anisotropic shale creep on the stress and permeability evolution of a geological nuclear waste repository

To ensure long-term safety and performance, geological nuclear waste repositories require low-permeability barriers such as bentonite buffers and/or shale host rock. Shale is not only known for its low permeability but also for its trend to undergo time-dependent deformation (i.e., creep), which could heal damage, but the effects of shale creep on the long-term performance of nuclear waste repositories have not been clearly understood. In particular, the anisotropic nature of shale (i.e., bedding) could have a significant effect on its creep behavior, and consequently, on the long-term performance of nuclear waste repositories. In this research, numerical simulations were carried out with the objective of showing the effects of anisotropic shale creep on the stress and permeability evolution of a generic geological nuclear waste repository in shale. The TOUGH-FLAC simulator, a thermo-hydromechanically (THM) coupled numerical code, was used for the simulations. To achieve the objective, comparisons were performed between the results of anisotropic shale creep simulations and those of different simulation cases, namely, no creep (i.e., elastic), isotropic creep, and long-term creep shale cases. Results of the comparisons show that the elastic and isotropic creep shale cases respectively led to the overestimation and underestimation of stress and permeability in the repository, whereas the long-term creep shale case, which accumulated greater creep in later periods than in earlier periods, helped to keep large shear and tensile stresses from developing while maintaining compressive spherical stress, resulting in consistently low permeability levels. These results indicate that performance assessments with elastic and isotropic creep formation models will provide the upper and lower bound estimates of stress and permeability, while more reasonable estimates will be given by an anisotropic creep formation model, and that shale with long-term creep characteristics will be beneficial in many aspects of the safety and performance of nuclear waste repositories.

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

Atomic view of photosynthetic metabolite permeability pathways and confinement in synthetic carboxysome shells

Carboxysomes are protein microcompartments found in cyanobacteria, whose shell encapsulates rubisco at the heart of carbon fixation in the Calvin cycle. Carboxysomes are thought to locally concentrate CO 2 in the shell interior to improve rubisco efficiency through selective metabolite permeability, creating a concentrated catalytic center. However, permeability coefficients have not previously been determined for these gases, or for Calvin-cycle intermediates such as bicarbonate (HCO$^–_3$), 3-phosphoglycerate, or ribulose-1,5-bisphosphate. Starting from a high-resolution cryogenic electron microscopy structure of a synthetic β-carboxysome shell, we perform unbiased all-atom molecular dynamics to track metabolite permeability across the shell. The synthetic carboxysome shell structure, lacking the bacterial microcompartment trimer proteins and encapsulation peptides, is found to have similar permeability coefficients for multiple metabolites, and is not selectively permeable to HCO$^–_3$ relative to CO 2 . To resolve how these comparable permeabilities can be reconciled with the clear role of the carboxysome in the CO 2 -concentrating mechanism in cyanobacteria, complementary atomic-resolution Brownian Dynamics simulations estimate the mean first passage time for CO 2 assimilation in a crowded model carboxysome. Despite a relatively high CO 2 permeability of approximately 10 -2 cm/s across the carboxysome shell, the shell proteins reflect enough CO 2 back toward rubisco that 2,650 CO 2 molecules can be fixed by rubisco for every 1 CO 2 molecule that escapes under typical conditions. The permeabilities determined from all-atom molecular simulation are key inputs into flux modeling, and the insight gained into carbon fixation can facilitate the engineering of carboxysomes and other bacterial microcompartments for multiple applications.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Characterizing Biomass Feedstock Transport Properties Using State of the Art Imaging and Computational Techniques

The microstructure of lignocellulosic biomass determines heat and mass transfer during conversion processes. We present a novel method for characterizing the transport properties of biomass using advanced imaging and computational techniques. The microstructure of two woody feedstocks, red oak and Douglas fir, before and after pyrolysis, is revealed using X-ray computed tomography (XCT). Transport properties are calculated from the XCT images, and principal permeability tensors are calculated using an immersed boundary-based finite volume solver to model gas flow through the geometries. We observe that the permeabilities of native biomass are distinctly anisotropic, however, this anisotropy is greatly reduced after pyrolysis.

adaptive mesh refinement↗
DOE OSTI2024

DECOVALEX-2023: Task F1 Final Report

DECOVALEX-2023 Task F is a comparison of models and methods for post-closure performance assessment (PA) of a deep geologic repository for radioactive waste. The general aims of Task F are to build confidence in the models, methods, and software used for PA and to stimulate additional research and development in PA methodologies. The task objectives are to motivate development of PA modelling skills and capabilities, to examine the influence of model choices on calculated repository performance, and to compare the uncertainties introduced by model choices to other sources of uncertainty. Task F involves no actual experiment or site. It is a PA modelling exercise that requires the conceptual development of hypothetical repository designs and geologic settings. Because three of the teams were interested in salt and the rest of the teams were interested in crystalline rock, Task F was split into two branches: Task F1 for crystalline rock and Task F2 for salt. This report is for Task F1, crystalline rock. Teams from seven countries (Canada, Czech Republic, Germany, Korea, Sweden, Taiwan, and United States) participated in Task F1. The teams worked together to define the features, events, and processes of the reference case repository and established a set of performance measures. In addition, they defined a set of benchmark problems designed to test and compare modelling capabilities for fracture flow and transport at different scales. The repository design and benchmark problems are documented in a Task Specification that evolved over time as the group honed the specifications. The benchmark problems verified that each team can aptly model flow and transport in fractured media in 1-, 2-, and 3-dimensions. Two general approaches were used for the 3-dimensional benchmarks: discrete fracture network (DFN) and equivalent continuous porous medium (ECPM). DFN modelling involves explicit meshing of each fracture while ECPM modelling aims to capture the effective porosity and directional permeability of each cell in a space-filling mesh as affected by intersecting fractures. In some models, a combination of the two is used, i.e., DFN for large known fractures and ECPM for the rest of the domain. Transport is solved by using either the advection-dispersion equation or particle tracking. Although some variation is observed among model breakthrough curves in the benchmark problems, there is strong agreement in breakthrough behaviour up to at least the 75 th percentile for all benchmarks. At the 90 th percentile, breakthrough results show larger differences, suggesting several models retain substantially higher fractions of tracer in regions of slower moving water. In addition to the flow and transport benchmarks, several teams completed the source term benchmark, verifying capabilities for modelling radionuclide decay and ingrowth, waste package breach, instant release fractions, fuel matrix degradation rates, and radionuclide solubility limitations. The reference case is conceptualized as a generic spent fuel repository at a depth of 450 m in fractured crystalline rock. The repository has 50 parallel backfilled drifts, each with 50 deposition holes 6 m apart. Each deposition hole contains a 4-PWR waste package and bentonite buffer. The rock domain is 5 km in length, 2 km in width, and 1 km in depth. It has 6 deterministic fractured deformation zones and a multitude of stochastic fractures. Teams generally used the ECPM approach for the entire rock or a hybrid approach in which the deterministic fracture zones are modelled with a DFN and the rest of the rock is modelled by ECPM. Of the reference case problems specified, only the results of the initial reference case problem are compared in this report. The initial problem focuses on transport from the deposition holes to the surface, i.e., it neglects waste package performance. Tracers are released at all waste package locations at time zero and tracked for their releases to the near field and ground surface. The water fluxes calculated at the ground surface entry and exit regions of the domain are similar for all models except for two that have considerably lower fluxes. For tracer transport, large differences are observed among models in the magnitude of tracer transported. Much of the difference appears to be due to how the repository is implemented and hence the different degrees of repository simplification. Models that exclude the drifts, buffer, and backfill from the domain tend to show greater release of tracers and radionuclides from the repository. The initial study presented here indicates that major differences in modelling important processes within the repository (e.g., diffusion through buffer and backfill) can produce broadly different release and transport results, especially when those processes are excluded. Even for the models that included all specified features, events, and processes, the results show significant differences and demonstrate the importance of examining multiple modelling approaches in performance assessment. The differences in results observed in this study are expected to motivate teams to either increase complexity in future versions of the reference case models or to improve methods to account for the effects of simplified features and processes. Either way, future improvements in these models are expected to produce results that more closely agree.

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

A Pro‐Angiogenic Immunoprotective Membrane for Cell Therapies

Abstract Immunoisolation strategies that rely on porous membranes play an important role in cell transplantation therapies to protect cells from the host's immune system. These membranes must possess immunoprotective properties while facilitating the transport of nutrients and cell products to maintain the functional integrity of encapsulated cells. An easy and scalable process is described to fabricate a dual function porous polymeric membrane that shields cells against immune cell attack and promotes vascularization to address the nutritional and oxygen requirements of transplanted cells. The fabrication process results in a membrane cross‐section with a gradient of nanopores to micropores that support cell immunoisolation and interfacial vascularization requirements, respectively. The membranes demonstrate excellent cell compatibility and effectively prevent T cell transmigration without compromising glucose diffusion and oxygen permeability. In a murine subcutaneous implantation model, membranes are stable for 60 days and exhibit significantly reduced fibrous capsules, with enhanced vascularization near the membrane. These porous polymeric membranes can potentially be used as pro‐angiogenic immunoprotective membranes for cell transplantation applications where maximizing cell viability and function is of critical importance.

Engineering↗
DOE OSTI2026

Computer simulation of reaction and transport of core-scale serpentinization of Fe-bearing olivine and evolution of geological hydrogen

We present reactive flow and transport simulations of core-scale serpentinization of forsteritic, Fe-bearing olivine under hydrothermal conditions. The model captures fluid flow, evolving porosity–permeability, and redox-controlled H 2 production in a cylindrical core over 55 days at 245 °C and 37 bar. Key processes include olivine dissolution, precipitation of serpentine and magnetite, and oxidation of Fe2+ to Fe 3+ . Here, results show strong coupling among flow velocity, alteration front propagation, and spatially heterogeneous H 2 generation. Elevated H2 concentrations align with inlet-localized magnetite precipitation, consistent with experiments. Secondary mineral formation reduces porosity and permeability, altering transport pathways. Reactive flow at the fluid–solid interface governs H 2 generation rates and distribution. Dissolved SiO2 promotes serpentine and talc formation by suppressing brucite, while bicarbonate extends reaction duration and moderates surface complexation effects.

Geologic hydrogen↗
DOE OSTI2024

Machine learning-based discovery of molecular descriptors that control polymer gas permeation

While machine learning has found increasing use in predicting the properties of polymeric materials with only a knowledge of chain architecture, determining the molecular factors underpinning properties (“interpretable AI”) has remained less well explored. We show that encoding chain chemistry in commonly employed formats, e.g., binary-valued fingerprints, leads to uniqueness issues during the hashing process to save storage space. This is because the hashing algorithm can map several chemical moieties into the same bit. These issues carry over into the ML algorithms, especially for “inverse” design and interpretable AI, and cannot be avoided by changing the length of the fingerprint. Using MACCS key featurizations of monomer repeats resolves some of these issues, and we show that a few substructures consistently appear in top features for maximizing permeability across several gases and ML models. These are carbon–carbon double bonds (as in polyacetylenes) especially when they are associated with methyl groups (found in branching architectures). Here these results, derived from the limited data set of ~ 500 polymers with experimental gas permeation data, are in agreement with physical insight and thus provide a robust foundation which could further enable study of these material classes through detailed experiments and simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Pore‐Scale Modeling of Reactive Transport with Coupled Mineral Dissolution and Precipitation

Abstract We present a new pore‐scale model for multicomponent advective‐diffusive transport with coupled mineral dissolution and precipitation. Both dissolution and precipitation are captured simultaneously by introducing a phase transformation vector field representing the direction and magnitude of the overall phase change. An effective viscosity model is adopted in simulating fluid flow during mineral dissolution‐precipitation that can accurately capture the velocity field without introducing any empirical parameters. The proposed approach is validated against analytical solutions and interface tracking simulations in simplified structures. After validation, the proposed approach is employed in modeling realistic rocks where mineral dissolution and precipitation are dominant at different locations. We have identified three regimes for mineral dissolution‐precipitation coupling: (a) compact dissolution‐precipitation where dissolution is dominant near the inlet and precipitation is dominant near the outlet, (b) wormhole dissolution with clustered precipitation where dissolution generates wormholes in the main flow paths and precipitation clogs the secondary flow paths, and (c) dissolution dominant where all solid grains are gradually dissolved. In the three regimes, the proposed approach provides reliable porosity‐permeability relationships that cannot be described well by traditional macroscale models. We find that the permeability can increase while the overall porosity decreases when the main flow paths are expanded by dissolution and adjacent pore spaces are clogged by precipitation.

58 GEOSCIENCES↗
DOE OSTI2023

Improving multiwell petrophysical interpretation from well logs via machine learning and statistical models

Well-log interpretation estimates in situ rock properties along well trajectory, such as porosity, water saturation, and permeability, to support reserve-volume estimation, production forecasts, and decision making in reservoir development. However, due to measurement errors, variability of well logs caused by multiple measurement vendors, different borehole tools, and nonuniform drilling/borehole conditions, estimations of rock properties with original well logs without proper preprocessing may not be accurate, especially in the context of multiwell estimation. Well-log normalization techniques such as two-point scaling and mean-variance normalization are commonly used to improve the robustness of multiwell rock-property estimation. However, these techniques do not consider the correlation between well logs and require subjective knowledge for their effective implementation. To reduce uncertainties and processing time associated with multiwell rock-property estimation from well logs, we develop discriminative adversarial (DA) and linear constraint models for well-log normalization and rock-property estimation. The DA neural network model developed for well-log normalization and interpretation can perform linear and nonlinear well-log normalization while considering the joint distribution of each well log and rock properties. However, the linear constraint model uses an ensemble of predictions from linear models to constrain well-log normalization and rock-property estimation. We also develop a divergence-based type well identification method to select type (training) wells for a test well based on the statistical similarity of associated well-log distributions instead of the interwell distance. We apply the DA model to perform well-log normalization and prediction of permeability for the Seminole San Andres Unit carbonate reservoir. Compared with the permeability predicted with the classical machine learning model without well-log normalization and models with two-point scaling normalization, the DA model yields the most accurate permeability prediction by decreasing the mean-squared error of permeability prediction by 20%–50%.

Geochemistry & Geophysics↗
DOE OSTI2021

Multimaterial Aerosol Jet Printed Magnetic Nanocomposites for Microwave Circuits

Here this article demonstrates the process and characterization of additively manufactured polymer-matrix magnetic nanocomposites. These nanocomposites are composed of nickel-zinc ferrite (Ni 0.5 Zn 0.5 Fe 2 O 4 ) nanoparticles in a polyimide matrix. These materials are printed using multimaterial aerosol jet printing where aerosols of nickel-zinc ferrite and polyimide inks are mixed in place to fabricate structures with nonunity relative permeability. To characterize these materials, we use energy-dispersive X-ray spectroscopy to estimate the content of ferrite in the composite, fitted simulation models, and ring resonator measurements to estimate the permeability. In fully aerosol jet printed 1-mm 2 inductors, we demonstrate an increase in inductance density of 40% in a 4.5-turn microstrip spiral inductor using this printed composite. These fully printed inductors demonstrate an inductance density of 4.2-nH/mm 2 . EDS analysis shows films with 18 wt.% Ni 0.5 Zn 0.5 Fe 2 O 4 . Finally, we show microstrip transmission lines on printed composite compared to the same on polyimide subjected to a dc magnetic bias, displaying a magnetic resonance behavior.

42 ENGINEERING↗