Search NASA⌕ Search

SEARCH · Search NASA

Results for “permeability variation”

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 19 records

Assessing reservoir performance for geologic carbon sequestration in offshore saline reservoirs

Carbon capture and storage (CCS) is a technological strategy to reduce CO 2 emissions from hydrocarbon-based energy production and non-energy industrial emitters, like fertilizer and cement production. To expand the geographic opportunities for CCS, offshore geologic settings are increasingly considered for saline CO 2 storage, particularly where legacy oil and gas infrastructure may be repurposed for CO 2 transport and injection. In this context, the Gulf of Mexico (GoM), United States (U.S.), may offer tremendous opportunity for CCS because the surrounding states are among the largest CO 2 producers in the U.S. and there has been extensive oil and gas development over the last half century. Nevertheless, there remains significant uncertainty in how GoM reservoir/seal systems will respond to industrial-scale CO 2 injections, particularly in the context of meso-scale (1–10s of meters) permeability variations that are difficult to identify prior to injection. Such permeability variations can have significant impacts on CO 2 migration and fluid pressure propagation, both of which govern the overall effectiveness of the CCS project. This research uses ensemble simulation methods to quantify the influence of spatially variable and a priori unknown permeability fields on CO 2 plume development and fluid pressure propagation when CO 2 is injected at 1 MMt per year for ten years. Results show a number of characteristic patterns that reveal the substantial influence of both near- and far-field permeability on reservoir performance. By focusing on basin characteristics that are common among offshore basins worldwide, results from this study suggest that offshore CCS may be a feasible carbon management strategy in regions where offshore basins are proximal to industrial CO 2 sources.

42 ENGINEERING↗

Distributed Acoustic Sensing to Estimate the Permeability

Optical fiber in a borehole can be interrogated with distributed acoustic sensors (DAS) to capture fracture displacements with the potential to map surrounding fracture networks. We designed a laboratory experiment to test the capability of DAS to determine borehole flow characteristics, and we show that for the first time DAS can be used to remotely estimate permeability. Optical fiber was wrapped around a bead filled pipe and the pressure drop and flow velocity were measured to directly calculate permeability. A machine learning model using statistical features from continuous DAS estimated the bulk permeability. Fluid interactions with the permeable material demonstrate insufficient resolution using DAS amplitude-based measurements for estimating pressure drop to infer permeability. Variations in the spectral domain relate DAS measurements to the pressure drop and provide consistent permeability estimates. Resolution with DAS is sufficient to estimate permeability and provides a reliable method to monitor at depth in borehole conditions.

58 GEOSCIENCES↗

Three‐Dimensional Permeability Inversion Using Convolutional Neural Networks and Positron Emission Tomography

Abstract Quantification of heterogeneous multiscale permeability in geologic porous media is key for understanding and predicting flow and transport processes in the subsurface. Recent utilization of in situ imaging, specifically positron emission tomography (PET), enables the measurement of three‐dimensional (3‐D) time‐lapse radiotracer solute transport in geologic media. However, accurate and computationally efficient characterization of the permeability distribution that controls the solute transport process remains challenging. Leveraging the relationship between local permeability variation and solute advection behavior, an encoder‐decoder based convolutional neural network (CNN) is implemented as a permeability inversion scheme using a single PET scan of a radiotracer pulse injection experiment as input. The CNN can accurately capture the 3‐D spatial correlation between the permeability and the radiotracer solute arrival time difference maps in geologic cores. We first test the inversion accuracy using synthetic test datasets and then test the accuracy on a suite of experimental PET imaging datasets acquired on four different geologic cores. The network‐predicted permeability maps from the geologic cores are used to parameterize forward numerical models that are directly compared with the experimental PET imaging data. The results indicate that a single trained network can generate robust 3‐D permeability inversion maps in seconds. Numerical models parameterized with these permeability maps closely capture the experimentally observed solute arrival time behavior. This work provides an unprecedented approach for efficiently characterizing multiscale permeability heterogeneity in complex geologic samples.

54 ENVIRONMENTAL SCIENCES↗

Bayesian Poroelastic Aquifer Characterization From InSAR Surface Deformation Data. Part I: Maximum A Posteriori Estimate

Characterizing the properties of groundwater aquifers is essential for predicting aquifer response and managing groundwater resources. In this work, we develop a high-dimensional scalable Bayesian inversion framework governed by a three-dimensional quasi-static linear poroelastic model to characterize lateral permeability variations in groundwater aquifers. In this work, we determine the maximum a posteriori (MAP) point of the posterior permeability distribution from centimeter-level surface deformation measurements obtained from Interferometric Synthetic Aperture Radar (InSAR). The scalability of our method to high parameter dimension is achieved through the use of adjoint-based derivatives, inexact Newton methods to determine the MAP point, and a Mat´ern class sparse prior precision operator. Together, these guarantee that the MAP point is found at a cost, measured in number of forward/adjoint poroelasticity solves, that is independent of the parameter dimension. We apply our methodology to a test case for a municipal well in Mesquite, Nevada, in which InSAR and GPS surface deformation data are available. We solve problems with up to 320,824 state variable degrees of freedom (DOFs) and 16,896 parameter DOFs. A consistent treatment of noise level is employed so that the aquifer characterization result does not depend on the pixel spacing of surface deformation data. Our results show that the use of InSAR data significantly improves characterization of lateral aquifer heterogeneity, and the InSAR-based aquifer characterization recovers complex lateral displacement trends observed by independent daily GPS measurements.

54 ENVIRONMENTAL SCIENCES↗

Predictive Inverse Model for Advective Heat Transfer in a Short–Circuited Fracture: Dimensional Analysis, Machine Learning, and Field Demonstration

Identifying fluid flow maldistribution in planar geometries is a well–established problem in subsurface science/engineering. Of particular importance to the thermal performance of enhanced (or “engineered”) geothermal systems is identifying the existence of nonuniform (i.e., heterogeneous) permeability and subsequently predicting advective heat transfer. Here, machine learning via a genetic algorithm (GA) identifies the spatial distribution of an unknown permeability field in a two–dimensional Hele–Shaw geometry (i.e., parallel plates). The inverse problem is solved by minimizing the L2 norm between simulated residence time distribution (RTD) and measurements of an inert tracer breakthrough curve (BTC) (C–Dot nanoparticle). Principal component analysis (PCA) of spatially correlated permeability fields enabled reduction of the parameter space by more than a factor of 10 and restricted the inverse search to reservoir–scale permeability variations. Thermal experiments and tracer tests conducted at the mesoscale Altona Field Laboratory (AFL) demonstrate that the method accurately predicts the effects of extreme flow channeling on heat transfer in a single bedding–plane rock fracture. However, this is only true when the permeability distributions provide adequate matches to both tracer RTD and frictional pressure loss. Without good agreement to frictional pressure loss, it is still possible to match a simulated RTD to measurements, but subsequent predictions of heat transfer are grossly inaccurate. Here, the results of this study suggest that it is possible to anticipate the thermal effects of flow maldistribution, but only if both simulated RTDs and frictional pressure loss between fluid inlets and outlets are in good agreement with measurements.

42 ENGINEERING↗

AnisONet: A deep neural operator-based anisotropic permeability upscaler from pore to Darcy scale

Directional permeability variations, which govern directional fluid flow in porous media with anisotropy, are important to accurately predict flow behavior, reactive transport, and fluid–solid interactions for various processes such as enhanced geothermal systems, energy storage devices, and biological systems. However, the intricate architecture of porous media makes it difficult to predict directional permeabilities. In this work, we present a novel machine learning (ML) framework, AnisONet, built upon an integration of a convolutional neural network, Swin transformer, and the deep operator network architecture, designed to predict anisotropic permeability and upscale predictions to larger spatial domains. First, AnisONet was evaluated with three classes of two-dimensional (2D) porous media, including synthetic circular and elliptical grains and natural sandstone grains from micro-computed tomography images. A lattice Boltzmann model (LBM) was used to calculate directional permeabilities at every 10° angle, producing 19 data points per image of porous media. AnisONet is then trained to predict permeability as a function of rotation angle. AnisONet showed strong predictive capability of directional permeability. Second, we tested our model for five upscaling cases with a large image size in the finite-element method (FEM) for 2D Darcy flow with various permeability tensor construction methods. Overall, upscaled permeability tensors in FEM simulations produce a reasonably good match with LBM results, highlighting the importance of selecting appropriate tensor formation strategies for accurate permeability upscaling. AnisONet, as a directional permeability estimator, could be further developed for more complex geometries, with the potential to develop a foundational ML model for various applications in porous media.

42 ENGINEERING↗

Red-green-bleached redox interfaces in the proximal Permian Cutler red beds: implications for regional fluid alteration

Siliciclastic strata of the Colorado Plateau attract attention for their striking red, green, bleached, and variegated colors that potentially record both early depositional and later diagenetic events. We investigated the proximal-most strata of the Paradox Basin, from their onlap contact with the Precambrian basement of the Uncompahgre Plateau to the younger Cutler strata exposed within 10 km of the Uncompahgre Plateau to attempt to understand the significance of the striking colors that occur here. These strata preserve a complex geology associated with buried paleorelief and sediment-related permeability variations at a major basin-uplift interface. Strata exposed within ~1.5 km of the onlap contact exhibit a pervasive drab color in contrast to the generally red colors that predominate farther from this front. In-between, strata commonly host variegated red/green/bleached intercalations. Thin-section petrography, SEM, XRD, Raman spectroscopy, Mössbauer spectroscopy, and whole-rock geochemistry of samples representing different color variations from demonstrate that water–rock interactions charged the rocks with Fe(II) that persists primarily in the phyllosilicate fraction. Color variations reflect grain-size differences that allowed the reduction of fluids from regional fault and basement/fill contacts to permeate coarser-grained Cutler sediments. Hematite and chlorite occur in both red and green sediments but are absent in the bleached sediments. Pervasive hematite in both red and green layers suggests that sediments were hematite-rich before later alteration. Chlorite and smectite are elevated in green samples and inversely correlated with biotite content. Green coloration is generally associated with 1) coarser grain sizes, 2) spatial association with basement contacts, 3) elevated smectite and/or chlorite, 4) less total Fe but greater Fe(II)/Fe(III) primarily in the phyllosilicate fraction, and 5) uranium enrichment. The bleached coloration reflects the removal of pigmentary Fe(III) oxide, while the green coloration is due to the removal of pigmentary hematite and the abundance of Fe(II)-bearing phyllosilicates. Abundant mixed-layer and swelling clays such as smectite, illite/smectite, and chlorite/smectite (including tosudite) dominate the mineralogy of the clay fraction. These results are consistent with other studies demonstrating fault-associated fluid alteration in the Paradox Basin region. However, the pervasive greening was not observed in many of these studies and appears to reflect the unique aspects of the paleovalley system and the importance of biotite alteration to Fe(II)-bearing phyllosilicates.

58 GEOSCIENCES↗

Simulation of Reverse Osmosis Membrane Compaction Using Material Point Method

Reverse Osmosis (RO) is a promising technology to address the impending water-crisis in the upcoming decades. RO at high salinities and high pressures is challenging due to membrane compaction that changes its porosity and permeability. In this study, we present a simulation methodology for membrane structural mechanics to understand pore size distribution and permeability variations under high pressure. We use the material point method (MPM), that solves the solid mechanics equations in a Lagrangian framework. MPM provides many features that make it well-suited for simulating mesoporous membranes. The Lagrangian framework allows for large deformations, easy integration of constitutive models and direct import of complex geometries as particles. The spatial discretization in our MPM solver is achieved using linear or cubic-spline shape functions while the time integration is carried out using the explicit Euler method. A series of images containing detailed pore structures obtained from X-ray tomography is first converted to a collection of material points. Compressive loads are applied to the top layer of the membrane to simulate the application of pressure. The membrane deformation and pore size distribution before and after load application are reported and compared with the experimental measurements. The presentation discusses the numerical methods used, the performance of the solver on high-performance computing machines, and the results of membrane compact in detail.

compaction↗

3D Simulations of gas injection on callovo-oxfordian claystone assuming spatial heterogeneity and anisotropy

A series of gas injection tests on Callovo-Oxfordian (COx) claystone from the Bure underground research laboratory (URL) in France were carried out at the British Geological Survey (BGS). The tests were performed using a triaxial apparatus specifically designed to capture small volumetric strains induced by the injected gas flow and consequent material dilatancy. The long-duration experiments were monitored throughout. Measurements also included pressure, stresses (axial and radial stresses prescribed for each test stage), rate of gas inflow, gas outflow volume as well as pore-pressures observed at various points of the sample. A coupled hydro-gas-mechanical 3D numerical model has been developed to simulate the tests. Initial permeability is assumed heterogeneous throughout the specimen and embedded fractures are incorporated in the formulation. Gas pressure-induced deformations during the test lead to variations of permeability due to changes in matrix porosity and, especially, fracture aperture as well as fracture orientation due to material anisotropy. A programme of sensitivity analyses involving the variation of different aspects and parameters of the model contributes to a better understanding of the phenomena and highlights its complexity. The model is able to reproduce the observed behaviour of the tests.

58 GEOSCIENCES↗

Grain detachment and transport clogging during mineral dissolution in carbonate rocks with permeable grain boundaries

In this paper, we explore the effect of permeable grain boundaries on carbonate rock dissolution in order to improve our understanding of grain detachment and migration during reactive flow. To do so we investigated the effects of grain size and enhanced reactivity along grain boundaries on global dissolution kinetics. Variations of permeability and porosity were calculated, and exponential relationships were observed in fractured rocks. Our model employed a reactive transport framework based on the Darcy-Brinkman equation to simulate calcite dissolution in carbonate rocks composed of microporous grains and, included fluid transport along grain boundaries. The model includes fluid flow, solute transport by advection–diffusion, heterogeneous reaction between fluid and minerals and grain detachment with subsequent grain transport in macropores. The migration of solid particles due to dissolution was based on cluster analysis and local movement, and the results show that grain detachment can lead to significant decreases in permeability due to the clogging of transport pathways. Microporous media with smaller grain sizes (fine grains in this study) showed a higher average reaction rate than those with coarser grains. In addition, as the overall rates of geochemical processes are commonly affected by the presence of texture heterogeneities such as fractures, a single fracture (macropore) introduced into a microporous matrix composed of permeable grains and grain boundaries was modelled as a large channel connecting the inlet and outlet of the simulation domain. It was found that macropore clogging by grain detachment temporary decreased permeability, lowering the long-term global reaction rates. Lastly, we observed that increases in flow rate can reduce detachment and fracture clogging by reducing local dissolution along grain boundaries.

58 GEOSCIENCES↗

Vertical Characterization of Flow and Contaminant Concentration in New Hanford Extraction Wells

Effective groundwater remediation should focus on high-concentration, low-permeability zones within the aquifer, as these areas often present significant challenges to contaminant removal. Characterizing the vertical distribution of contamination in an aquifer is therefore critical for optimizing and monitoring pump-and-treat (P&T) systems to meet target cleanup timelines. Flow patterns driven by hydraulic gradients and variations in permeability are described to provide insight into contaminant mass removal from long-screened wells. The objective of this multi-phase work is to provide actionable recommendations for improving remedy performance monitoring and P&T operations, which ongoing evaluations will address. This report summarizes the preliminary results of vertical flow and contaminant distribution monitoring in new extraction wells at the Hanford Site.

54 ENVIRONMENTAL SCIENCES↗

A micromechanical hydro-mechanical-damage coupled model for layered rocks considering multi-scale structures

An anisotropic Hydro-Mechanical-Damage coupled model was established in the framework of thermodynamics for saturated layered rocks, which contain arbitrarily-distributed cracks of smaller scales and orientated bedding planes of much larger scales. Here, the proposed model considers the anisotropic damage growth and frictional sliding of cracks and shear sliding and dilatancy of bedding planes, as well as the interaction between cracks and bedding planes. Numerical implementation of the nonlinear coupled model was developed based on TOUGHREACT, which is a well-established simulator for multiphase fluid flow and reactive transport analysis. The proposed model was first validated with laboratory water injection test results on rock blocks, and then used to study the effects of injection rates, in situ stresses and bedding planes on fluid injection processes at the field scale, demonstrating the importance to consider the variation of multi-scale structures for better understanding the coupled Hydro-Mechanical responses in layered rocks.

58 GEOSCIENCES↗

Unsteady- and Steady-State Relative Permeability Study with X-ray and Acoustic Monitoring for CO 2 Storage in Deep Saline Aquifers

In this study, we link the multiphase flow measurements with controls of sedimentary structures (e.g., heterogeneity and anisotropy) on relative permeability to variations in ultrasonic velocities for two deep saline aquifer formations (Entrada and Bluff Sandstones) in the San Juan Basin of the Southwestern USA. The rock specimens were extracted from outcrop sites near Durango, CO, USA. They have distinct differences in grain size, cementation composition, and individual chemical amounts, despite both formations being eolian sandstones. We performed a series of unsteady- and steady-state CO 2 -brine relative permeability experiments under capillary-controlled displacement rates. Unsteady-state experiments were conducted at 71 °C and 9.65 MPa; steady-state experiments were conducted at 85 °C and 22.8 MPa and 89 °C and 24.1 MPa for the Bluff and Entrada Sandstones, respectively. During the unsteady experiments, X-ray computed tomography was used to visualize multiphase flow in porous media and quantify saturations during brine drainage under various flow rates. Scan images and saturation profiles indicate that the CO 2 distribution in the pore volume was strongly impacted by the presence of high-angle cross-laminations, heterogeneous rock structure, and direction of bedding orientation. Those factors contribute to dramatic and quick initial breakthroughs and affect the overall saturation dynamics. Steady-state relative permeability tests were conducted at net flow rates of 1 mL/min for both brine drainage and imbibition. During the experimental steps, the CO 2 fractional flow was increased and decreased for both drainage and imbibition scenarios to mimic the front when CO 2 contacts brine and behind the front when brine enters space previously occupied by CO 2 . It was found that compressional velocity decreased, while shear waves slightly increased as brine saturation decreased. The hysteresis effects for the relative permeability and acoustic velocities were distinct. Furthermore, the CO 2 /brine front stability is quantified by applying a mobility ratio approach to spot saturations at which the boundary line between fluids becomes uneven. The results presented in this work can potentially boost the quality and precision of forecasts for the CO 2 storage projects in which the vertical and horizontal core-scale heterogeneity and anisotropy impact the plume migration within host reservoirs.

carbon dioxide (CO2)↗

Nanorod Pair Complexes Manipulated via Magnetic Casimir Forces

Here, we demonstrate a route to modulate Casmir–Lifshitz interactions between anisotropic nanoparticles and magnetic fluids. By semiclassical quantum electrodynamics, we study ground state dispersion forces for cylindrical dielectric nanorods made of polystyrene (PS) and zinc oxide (ZnO) embedded in toluene-based host media with gold-coated magnetite nanoparticles and also predict magnetic contributions to the fully retarded excited state interaction. The variation in magnetic permeability enables tuning between repulsive and attractive interactions, and measurable magnetic Casimir traps are predicted between a pair of ZnO–PS nanoparticles whose equilibrium position can be modulated over an order of magnitude with a small variation in the size of the magnetite nanoparticle. This provides an alternative magnetic Casimir-effect pathway to reversibly tune quantum electromagnetic forces at the nanoscale for the assembly and enhancement of colloidal stability.

Excited-state interaction↗

Permeability-controlled migration of induced seismicity to deeper depths near Venus in North Texas

Abstract Migration of seismic events to deeper depths along basement faults over time has been observed in the wastewater injection sites, which can be correlated spatially and temporally to the propagation or retardation of pressure fronts and corresponding poroelastic response to given operation history. The seismicity rate model has been suggested as a physical indicator for the potential of earthquake nucleation along faults by quantifying poroelastic response to multiple well operations. Our field-scale model indicates that migrating patterns of 2015–2018 seismicity observed near Venus, TX are likely attributed to spatio-temporal evolution of Coulomb stressing rate constrained by the fault permeability. Even after reducing injection volumes since 2015, pore pressure continues to diffuse and steady transfer of elastic energy to the deep fault zone increases stressing rate consistently that can induce more frequent earthquakes at large distance scales. Sensitivity tests with variation in fault permeability show that (1) slow diffusion along a low-permeability fault limits earthquake nucleation near the injection interval or (2) rapid relaxation of pressure buildup within a high-permeability fault, caused by reducing injection volumes, may mitigate the seismic potential promptly.

58 GEOSCIENCES↗

A comparative investigation of the effect of gas type on foam strength and flow behavior in tight carbonates

Foam generation technique has been practically applied to overcome gas mobility issue during enhanced oil recovery (EOR) processes. Incorporated with the utilized gas, the capability of foaming agent to generate persistent foam determines the flow performance under reservoir conditions. Here, in this study, we first identified the surfactants with excellent phase behavior and foaming properties and, second, evaluated their bulk foam performance with different gases, including N 2 , CO 2 , and hydrocarbon gas (CH 4 ) under high-pressure and high-temperature conditions. In bulk, the foam generated with the high-density gas, i.e., CO 2 , exhibited the greatest foamability but the lowest foam stability. With amidoamine oxide-based surfactant, CH 4 -foam was found to be more stable than N 2 -foam, while the opposite was observed for sulfobetaine-based surfactant. The most effective surfactant was later used to perform flow experiments on two tight carbonate rocks, including Minnesota Norther Cream Buff and Edwards limestone, with permeabilities of ~0.8 and 16 mD, respectively. The core-flooding experiments showed that the optimum quality of CO 2 -foam was lower than that of N 2 - and CH 4 -foams. For a given flow rate and fixed foam quality, the steady-state strength followed the descending order of CH 4 -, N 2 -, and CO 2 -foam. However, in general, the difference in foam strength between N 2 - and CH 4 -foams was insignificant when the total flow rate was varied. We also found that the variation in rock permeability had a greater effect on the strength of the CO 2 -foam compared to those generated with N 2 and CH 4 . In the presence of oil, the shear-thinning behavior of all foams was found more prominent. The increase in oil saturation had a more adverse effect on the strength of N 2 - and CH 4 -foams, particularly in the higher permeability rock. The results indicated that the foam strength, which is key in the effectiveness of foam EOR applications, is greatly influenced by several factors such as gas density, gas solubility in water, the presence of oil, and rock permeability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An extended trajectory-mechanics approach for calculating two-phase flow paths

A technique originating in quantum dynamics is used to derive a trajectory-based, semi-analytical solution for two-phase flow. The partial differential equation governing the evolution of the aqueous phase is equivalent to a family of ordinary differential equations defined along a path through the porous medium. The trajectories may be found by solving the differential equations directly or by post-processing the output of a numerical solution to the full set of governing equations. The trajectories, which differ from conventional streamlines, are found to bend downward in response to gravitational forces. The curvature is more pronounced as the dip of the porous layer containing the flow increases. Subtle changes in the relative permeability curve can lead to significant variations in the trajectories. The ordinary differential equation for the trajectory provides an expression for the travel time along the path. The expression produces a semi-analytical approximation to the model parameter sensitivities, the partial derivatives of the travel times with respect to changes in the permeability model. The semi-analytical trajectory-based sensitivities generally agree with those computed using a numerical reservoir simulator and a perturbation approach. The sensitivities are useful in tomographic imaging algorithms designed to estimate the spatial variation in permeability within a porous medium using multiphase observations.

58 GEOSCIENCES↗

Caprock Remains Water Wet Under Geologic CO 2 Storage Conditions

Abstract Carbon storage technology is primarily targeted in saline formations, which is a porous rock matrix filled with brine, sealed with a low permeability caprock. There are significant variations of CO 2 wetting properties, typically reported in the literature as contact angle of CO 2 and brine interacting with a rock material, suggesting that CO 2 could become wetting under geostorage conditions and negatively impact containment effectiveness. Here, we performed the first controlled laboratory measurements of CO 2 ‐brine contact angles on shale rocks from low permeability sealing formations with distinctive mineralogic properties—calcite‐rich, quartz‐rich, and dolomite‐rich. We targeted temperatures at 40° and 100°C, pressures at 8.3, 34.5, and 62.1 MPa, and salinity at 35,000 and 260,000 ppm. Results show no significant change in contact angle with mineralogy, temperature, pressure, salinity, and CO 2 bubble size. We conclude that caprocks will remain water‐wet at geologic CO 2 storage conditions and will not lose their capillary sealing capacity.

Tapriyal, Deepak↗