SEARCH · Search NASA

Results for “change point”

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.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Catalytic Water Electrolysis by Co–Cu–W Mixed Metal Oxides: Insights from X-ray Absorption Spectroelectrochemistry

Mixed metal oxides (MMOs) are a promising class of electrocatalysts for the oxygen evolution reaction (OER) and hydrogen evolution reaction (HER). Despite their importance for sustainable energy schemes, our understanding of relevant reaction pathways, catalytically active sites, and synergistic effects is rather limited. Here, we applied synchrotron-based X-ray absorption spectroscopy (XAS) to explore the evolution of the amorphous Co–Cu–W MMO electrocatalyst, shown previously to be an efficient bifunctional OER and HER catalyst for water splitting. Ex situ XAS measurements provided structural environments and the oxidation state of the metals involved, revealing Co 2+ (octahedral), Cu + / 2+ (tetrahedral/square-planar), and W 6+ (octahedral) centers. Operando XAS investigations, including X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS), elucidated the dynamic structural transformations of Co, Cu, and W metal centers during the OER and HER. The experimental results indicate that Co 3+ and Cu 0 are the active catalytic sites involved in the OER and HER, respectively, while Cu 2+ and W 6+ play crucial roles as structure stabilizers, suggesting strong synergistic interactions within the Co–Cu–W MMO system. In conclusion, these results, combined with the Tafel slope analysis, revealed that the bottleneck intermediate during the OER is Co 3+ hydroperoxide, whose formation is accompanied by changes in the Cu–O bond lengths, pointing to a possible synergistic effect between Co and Cu ions. Our study reveals important structural effects taking place during MMO-driven OER/HER electrocatalysis and provides essential experimental insights into the complex catalytic mechanism of emerging noble-metal-free MMO electrocatalysts for full water splitting.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Measured Long-Term Solar Irradiance for Climate Studies

Solar radiation is primarily measured using high-quality radiometers (e.g., pyranometers and pyrheliometers). These instruments need to be calibrated regularly (every two years at a minimum). They are also susceptible to various sources of uncertainties. Therefore, rigorous data quality assessment is required to obtain high-confidence data from these radiometers. This is especially true when the data is used to understand climatic trends and/or extreme weather events. In this study, we attempted to create a continuous and reliable dataset by correcting the underlying data which we believe contains bias due to reference pyranometers swap procedures. These biases can be significant which can be up to two percentage points thereby influencing the interpretation of climate changes and/or extreme events. This study elaborates on the biases and methods to correct those biases.

data integrity↗

DOE OSTI2024

Statistical Behavior of Low-Amplitude Power System Point-on-Wave Measurements

The power grid is undergoing massive changes to ensure resiliency and reliability in a more decentralized world. Distributed energy resources are becoming a prominent source of generation, potentially leading to a lack of centralized generation sources. Due to these new behaviors and system topologies, it is important to install measurement devices that are 1) accurate and 2) self-aware of their measurement quality. In this paper, a residential-scale microgrid is used to generate voltage and current waveforms, captured by Verivolt and National Instruments measurement equipment. A least-squares approach is used to separate the “clean” signals from the noise. Finally, Gaussian mixture modeling is used to approximate noise distributions, and it is shown these higher-order distribution estimates are a better fit to voltage and current noise profiles than single-mode Gaussian estimates.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

Property-Graded Stainless Steel 316L by Selective Laser Melting: Characterization & Design

Abstract The purpose of this research work is to characterize and inform the design of (mechanical) property-graded bulk structures made from a single metallic alloy via a laser powder bed fusion (LPBF) process, with an end goal of creating repeatable/reproducible functionally-graded additively manufactured (FGAM) parts. This paper specifically investigates the manufacture of stainless steel (SS) 316L structures via a pulsed selective laser melting (SLM) process, and the underlying causes of property variations (within a functionally-acceptable range) through various material characterization techniques. For this, a design of experiments spanning the volumetric energy density (VED) based process parameter design space was utilized to investigate the range of functionally-acceptable physical/mechanical properties achievable in SS 316L. Five sample conditions (made via different process parameter combinations) were down-selected for in-depth microstructure analysis and mechanical/physical property characterization; these were suitably selected to impart a wide and controllable property range (209–318 HV hardness, 90–99.9% relative density, and 154–211 GPa modulus). It was observed that property variations resulted from combinations of porosity types/amounts, martensitic phase fractions, and grain sizes. Based on these findings, property-graded standard test specimens were designed and manufactured for further investigation—tensile specimens having a monotonic hardness change along its gauge length, four-point bending specimens with varying elastic moduli as a function of the distance from the neutral axis, and Moore’s rotating beam fatigue specimens with moduli variations based on the distance from the center. Altogether, this work lays the foundation for understanding and designing the local and global mechanical performance of FGAM bulk structures.

Engineering↗

DOE OSTI2025

Energetic Electron Transport in Magnetized Plasma with Magnetic Islands (Final Technical Report)

The main scientific goal of this project was to investigate the interactions between energetic electrons and magnetic islands in magnetized plasmas from lab to space. The proposed plan of addressed the following specific objectives: (1) Determine if energetic electrons are trapped by magnetic islands. Specifically, test the hypothesis that energetic electrons are trapped near island O-points due to stable orbits and de-confined near island X-points due to X-point tangles. (2) Determine how electron transport changes with island width and location. Specifically, test the hypothesis that electrons are accelerated by contracting magnetic islands through a Fermi acceleration process. (3) Assess how electron transport is affected by island dynamics, including island rotation, bifurcation, and/or island overlap. Specifically, test the hypothesis that electrons are deconfined and possibly accelerated during island rotations, bifurcations, and overlap due to stochastization of the magnetic field lines.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

Hot-Carrier Photoluminescence in Low-Temperature CsSnI3 Nanocrystals

In the halide perovskites, the B-site metal lone-pairs electrons are strongly suspected to underpin many of the interesting properties of the inorganic halide perovskites. Recently, a stable low-temperature monoclinic polar phase was predicted for CsSnI3, opening the possibility of direct investigation of both these lone pair electrons and a ferroelectric distorted structure. To date there are no other reports of such a structure in CsSnI3, and a known low-temperature monoclinic structure in CsSnBr3 remains unexplored. We have found optical evidence of a transformation occurring around 240 K in CsSnI3 nanocrystals, with several changes in optical behavior below this transition point, including novel high-energy photoluminescene and new states in the transient absorption spectrum. We have successfully characterized the optical properties of this low-temperature phase and found evidence for a polar, monoclinic structure. Discovery of a stable monoclinic polar structure in the halide perovskites opens many new potential directions for further research and electronics applications.

14 SOLAR ENERGY↗

DOE OSTI2025

Mesoscale Modeling of Dislocation Cell Structure Evolution and Radiation-Induced Segregation in Additively Manufactured Austenitic Stainless Steel

Structural alloys under irradiation develop radiation-induced segregation (RIS) at point defect sinks, which lead to undesired changes in the alloy's properties. Additively manufactured austenitic stainless steels are expected to yield a distinct response to irradiation damage owing to their unique as-printed dislocation cell structure. In this talk, we present the development of a mesoscale model for dislocation cell structure evolution coupled with RIS using the MOOSE framework. The evolution of network dislocations within spatially distinct cells and cell walls are modeled using climb-mediated edge annihilation and generation processes. Starting with initial microstructures of as-printed non-equilibrium segregations at dislocation cell walls, our 2D and 3D simulations explore the effects of irradiation temperature, dose rate, dislocation sink bias, and specimen thickness on microstructure evolution. The results will be compared against experimental characterizations of in-situ and ex-situ ion-irradiated samples, and the implications of evolving cell structure on irradiation damage response will be discussed.

additive manufacturing↗

DOE OSTI2025

Impact of lead on an axisymmetric, single bead blown powder DED overhung geometry

Metal additive manufacturing can be utilized for the near net-shape manufacture of components in a layer wise technique. Traditionally, the buildup process is conducted in the direction of gravity or bottom to top in the vertical orientation; however, with the availability of commercial systems with additive heads or fixturing that can index from vertical, and the need to manufacture larger parts without modifying available systems, the limitation of manufacturing from bottom to top is removed enabling the buildup of larger components printed at an angle from vertical. The effect of gravity on the melt pool and as-deposited component quality when printing off vertical is unknown in literature, especially for axisymmetric components. This understanding is critical for the advancement of manufacturing components of increased size in 4 and 5-axis additive systems. This investigation utilizes blown powder directed energy deposition to evaluate the change in as-printed geometry when the start point is altered in relation to gravity while the part rotates to manufacture an axisymmetric component. The objectives of this work are to determine the impact of the deposition location on the geometric variability on axisymmetric components. This investigation tests the hypothesis that differences in layer height exist due to a change in catchment efficiency when the deposition location is moved to a tangent of the round geometry due to a change in the melt pool dynamics. It was found the ideal deposition location when printing at 26.6-degrees from vertical was at −90-degrees from the top center point of the component, along the tangent, where gravity was pushing the melt pool down, but the rotation of the part was pulling the deposited material towards the top center of the component. This work provides an understanding of layer height stability and catchment efficiency to guide print orientation strategy for high-aspect ratio components. It was found the −90-degree lead deposition had the best layer height stability at 0.6 % as compared to the top-center at 6.3 % and +90-degree deposition location at 9.2 % for the nominal programmed layer height. The change in layer height also effected the final diameter the greatest for the top center deposition location where the diameter diverged by 2.5 % during the overbuild condition and converged by 0.6 %. In conclusion, this finding will increase the manufacturing efficiency of axisymmetric components by increasing the passive stability of the printing process for the successful manufacture of parts without the need for perfectly calibrated manufacturing parameters.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Combined ambient vibration and surface displacement measurements for improved progressive failure monitoring at a toppling rock slab in Utah, USA

Seismic resonance and surface displacement measurements can be implemented in tandem to improve landslide characterization and progressive failure monitoring. Crack aperture data are frequently used in rock slope stability monitoring and often exhibit recognizable trends prior to failure, such as accelerated crack opening. Alternatively, ambient resonance data offer multiple parameters including modal frequencies, damping, and polarization that can be monitored alongside crack aperture and may respond differently to environmental forcings and complex failure evolution. We analyzed data from continuous ambient vibration monitoring and concomitant crack aperture measurements at the Courthouse Mesa instability, a large toppling sandstone slab in Utah, USA. Three years of data revealed crack aperture increases of 2–4 mm/year with no clearly detectable irreversible changes in modal parameters, including frequency. Annually, frequency and displacement varied by 29% and 19% of the mean, respectively, with average and maximum daily frequency fluctuations of 6.5% and 16%, respectively. These reversible cyclic changes were primarily temperature-driven, but annually, frequency was in-phase with temperature whereas crack aperture lagged temperature changes by ~37 days. Polarization and damping also varied seasonally but were less strongly correlated with temperature. Conceptual 3D finite element modeling demonstrated consistent frequency decreases associated with crack propagation but variable changes in crack aperture measured at a single point; i.e., crack propagation did not always result in increased crack opening but always generated a resonance frequency decrease. Taken together, our data suggest a possible thermal wedging-ratcheting mechanism at the Courthouse Mesa instability, where annual thermoelastic crack closure is impeded by debris infill but the absence of downward crack propagation during the monitoring period is evidenced by no permanent resonance frequency changes. Our study demonstrates that combined seismic resonance and crack aperture data provide an improved description of rock slope instability behavior, supporting refined characterization and monitoring of changes accompanying progressive failure.

58 GEOSCIENCES↗

DOE OSTI2025

Development of an ERT‐Based Framework for Bentonite Buffers Monitoring From Laboratory Tests: 1. Characterizing Thermal–Hydrological–Mechanical Processes

Abstract Bentonite clay is widely used in engineered barrier systems for the permanent disposal of high‐level radioactive waste due to its low permeability, high swelling capacity, and thermal stability. However, the complex thermal‐hydrological‐mechanical (THM) processes induced by heating from decaying radioactive waste and hydration from surrounding rock can lead to heterogeneous changes that are difficult to measure and predict. This study develops an Electrical Resistivity Tomography (ERT)‐based framework for monitoring THM processes, progressing from sample‐scale to bench‐scale tests, to inform field‐scale applications. Sample‐scale tests analyzed small bentonite samples under controlled variations in water content, temperature, and porosity to establish fundamental resistivity relationships. Bench‐scale tests involved larger bentonite columns subjected to heating (up to 200°C) and hydration under controlled pressure, simulating repository conditions. ERT measurements, complemented by X‐ray CT imaging, temperature monitoring, and tracing sensors, revealed coupled THM processes, such as hydration‐induced compression, swelling, and thermal gradients, leading to complex resistivity patterns. The results demonstrate the potential of ERT for capturing THM‐induced resistivity changes, though challenges remain in upscaling and quantitative analysis. This study evaluates laboratory test capabilities and proposes future improvements for understanding THM‐induced resistivity responses. A conceptual framework for ERT implementation in field‐scale monitoring is presented, synthesizing findings from both scales and exploring how ERT data can inform long‐term modeling and reduce prediction uncertainties. Overall, this ERT‐based framework offers a robust method for monitoring bentonite buffers, aiding in early issue detection and supporting the safe long‐term disposal of radioactive waste in geological repositories, while highlighting the need for future development. Plain Language Summary Bentonite clay is crucial in engineered barrier systems (EBS) for containing high‐level radioactive waste due to its ability to absorb water, swell, seal and remain stable under high temperatures. When bentonite absorbs water and heats up from radioactive decay, it experiences complex changes in its physical and mechanical properties. Understanding these changes is important for ensuring the long‐term safety and effectiveness of EBS. This study used Electrical Resistivity Tomography (ERT), a non‐invasive method that measures electrical conductivity to monitor these changes during laboratory experiments. The ERT data revealed significant variations in resistivity corresponding to changes in water content, temperature, and density, providing detailed spatial and temporal insights into the behavior of bentonite. These findings enhance our ability to predict the long‐term performance of bentonite barriers, ensuring the safe containment of radioactive waste. By improving our understanding of bentonite's behavior, this research supports the development of more reliable and effective barrier systems for radioactive waste disposal, protecting the environment and public health. Key Points ERT monitoring was employed to capture resistivity changes in bentonite during controlled heating and hydration experiments, providing insights into THM processes ERT data reveal significant resistivity changes correlated with water content, temperature, and mechanical effects, enhancing the understanding of THM dynamics in bentonite This study explores the potential of the framework for application in field‐scale EBS monitoring, emphasizing the need for integrating additional geophysical methods for comprehensive subsurface imaging

Chen, Hang↗

DOE OSTI2023

Microorganisms, climate change, and the Sustainable Development Goals: progress and challenges

The forthcoming UN summit marks the halfway point to 2030 and presents an important milestone in global efforts to address various challenges, including those related to climate change and environmental preservation. The UN SDGs include several related to microorganisms and climate change. Importantly, microbiology research is key to understanding and mitigating climate change, and in maintaining the health of terrestrial and aquatic ecosystems (SDGs 13, 14 and 15).

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2025

Accounting for electron-beam-induced warping of molecular nanocrystals in MicroED structure determination

High-energy electrons induce sample damage and motion at the nanoscale to fundamentally limit the determination of molecular structures by electron diffraction. Using a fast event-based electron counting (EBEC) detector, we characterize beam-induced, dynamic, molecular crystal lattice reorientations (BIRs). These changes are sufficiently large to bring reciprocal lattice points entirely in or out of intersection with the sphere of reflection, occur as early events in the decay of diffracted signal due to radiolytic damage, and coincide with beam-induced migrations of crystal bend contours within the same fluence regime and at the same illuminated location on a crystal. These effects are observed in crystals of biotin, a series of amino acid metal chelates, and a six-residue peptide, suggesting that incident electrons inevitably warp molecular lattices. The precise orientation changes experienced by a given microcrystal are unpredictable but are measurable by indexing individual diffraction patterns during beam-induced decay. Reorientations can often tilt a crystal lattice several degrees away from its initial position before irradiation, and for an especially beam-sensitive Zn(II)-methionine chelate, are associated with dramatic crystal quakes prior to 1 e − Å −2 electron beam fluence accumulates. Since BIR coincides with the early stages of beam-induced damage, it echoes the beam-induced motion observed in single-particle cryoEM. As with motion correction for cryoEM imaging experiments, accounting for BIR-induced errors during data processing could improve the accuracy of MicroED data.

Vlahakis, Niko (ORCID:0000000250920265)↗

DOE OSTI2021

Hole-like Fermi surface in the overdoped non-superconducting Bi 1.8 Pb 0.4 Sr 2 CuO 6+δ

In high-temperature cuprate superconductors, the anti-ferromagnetic spin fluctuations are thought to have a very important role in naturally producing an attractive interaction between the electrons in the d-wave channel. The connection between superconductivity and spin fluctuations is expected to be especially consequential at the overdoped end point of the superconducting dome. In some materials, that point seems to coincide with a Lifshitz transition, where the Fermi surface changes from the hole-like centered at $(\pi, \pi)$ to the electron-like, centered at the Γ point causing a loss of large momentum anti-ferromagnetic fluctuations. In this work, we study the doping dependence of the electronic structure of Bi 1.8 Pb 0.4 Sr 2 CuO$_{6+\delta}$ in angle-resolved photoemission and find that the superconductivity vanishes at lower doping than that at which the Lifshitz transition occurs. This requires a more detailed re-examination of a spin fluctuation scenario.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

Chemical changes from N-doped graphene and Metal-Organic Frameworks to N-G/MOF composites for improved electrocatalytic activity

Integrating N doped graphene (N G) with Metal Organic Frameworks (MOFs) enhances catalytic activity for oxygen reduction reaction (ORR), often exceeding both the performances of their precursors and, in some cases, even PGM based catalysts. However, the factors driving this improved catalytic activity in N G/MOF composites remain unexplored, particularly from the perspective of the chemical changes. To investigate the chemical changes in N G and MOF upon their integration and the implications of these changes on ORR catalytic activity, an N G/MOF was synthesized from N G with a MOF (ZIF 8) following a mechanochemical wet ball milling process. The N G, ZIF 8, and N G/MOF samples were examined for changes in elemental composition, chemical state of carbon, different nitrogen and carbon bonds, and other chemical interactions. In the N G/MOF catalyst, compared to its N G and ZIF 8 precursors, the relative oxygen content increased, indicating the formation of additional oxygen containing groups. The C 1s peak shifted to a lower binding energy in N G/MOF, suggesting changes in the overall chemical or oxidation state of the carbon atoms. Besides, the increase in pyridinic N functional groups in N G/MOF points to the formation of additional active sites. Furthermore, the formation of C Zn bonds in N G/MOF suggests the probable emergence of single atom Zn sites, while the increase in C=O bonds points to the formation of carboxyl or carbonyl groups. These chemical changes could be linked to the enhanced electrocatalytic activity of the N G/MOF composite for ORR. In conclusion, this study may also be beneficial for other research focused on developing composite catalysts involving various N G and MOF materials.

25 ENERGY STORAGE↗

DOE OSTI2025

Dual Photoluminescence in Low-Temperature Phase of CsSnI 3 Nanocrystals

The expression of metal lone-pair electrons is hypothesized to underpin many of the interesting properties of inorganic halide perovskite semiconductors. Recently, a stable low-temperature monoclinic polar phase was predicted for CsSnBr 3 and CsSnI 3 , opening the possibility of direct investigation of a ferroelectric distorted structure compared to the undistorted structure. To date, there have been no experimental reports of such a structure in CsSnI 3 , and the low-temperature optical properties of CsSnI 3 nanocrystals have remained unexplored. Here we report optical and structural evidence of a phase transition around 240 K in 8.9 nm CsSnI 3 nanocrystals. Several changes in optical behavior occur below this transition point, including high-energy photoluminescence (PL) that emits concurrently with the exciton PL. The emergence of this high-energy PL is correlated with X-ray diffraction (XRD) and differential scanning calorimetry (DSC) supporting a phase transition from the orthorhombic structure between 240-200 K. Transient absorption measurements show an increase in the excited state lifetimes, i.e., slowed carrier cooling, at 200 K when photoexciting with photon energies above the high-energy state, consistent with slowed carrier cooling and emergence of high-energy PL. We hypothesize that the slowed carrier cooling is distinctive to this phase transition that modifies both the electronic and phonon structures that dictate excited-state carrier dynamics, and we discuss these changes.

14 SOLAR ENERGY↗

DOE OSTI2020

Topological Regularization via Persistence-Sensitive Optimization

Optimization, a key tool in machine learning and statistics, relies on regularization to reduce overfitting. Traditional regularization methods control a norm of the solution to ensure its smoothness. Recently, topological methods have emerged as a way to provide a more precise and expressive control over the solution, relying on persistent homology to quantify and reduce its roughness. All such existing techniques back-propagate gradients through the persistence diagram, which is a summary of the topological features of a function. Their downside is that they provide information only at the critical points of the function. We propose a method that instead builds on persistence-sensitive simplification and translates the required changes to the persistence diagram into changes on large subsets of the domain, including both critical and regular points. This approach enables a faster and more precise topological regularization, the benefits of which we illustrate with experimental evidence.

Nigmetov, Arnur↗

DOE OSTI2022

A Mountaintop View Requires Minimal Sorting: A Faster Contour Tree Algorithm

Consider a scalar field f : M → R, where M is a triangulated simplicial mesh in R d . A level set, or contour, at value v is a connected component of f –1 (v). As v is changed, these contours change topology, merge into each other, or split. Contour trees are concise representations of f that track this contour behavior. The vertices of these trees are the critical points of f, where the gradient is zero. The edges represent changes in the topology of contours. It is a fundamental data structure in data analysis and visualization, and there is significant previous work (both theoretical and practical) on algorithms for constructing contour trees. Suppose M has n vertices, N facets, and t critical points. A classic result of Carr, Snoeyink, and Axen (2000) gives an algorithm that takes O(n log n+Nα(N)) time (where α(·) is the inverse Ackermann function). A further improvement to O(t log t + N) time was given by Chiang et al. All these algorithms involve a global sort of the critical points, a significant computational bottleneck. Unfortunately, lower bounds of Ω(t log t) also exist. We present the first algorithm that can avoid the global sort and has a refined time complexity that depends on the contour tree structure. Intuitively, if the tree is short and fat, we get significant improvements in running time. For a partition of the contour tree into a set of descending paths, P, our algorithm runs in O($\Sigma$ pϵP |p| log |p| + tα(t) + N). This is at most O(t log D + N), where D is the diameter of the contour tree. Moreover, it is O(tα(t) + N) for balanced trees, a significant improvement over the previous complexity. Our algorithm requires numerous ideas: partitioning the contour tree into join and split trees, a local growing procedure to iteratively build contour trees, and the use of heavy path decompositions for the time complexity analysis. There is a crucial use of a family of binomial heaps to maintain priorities, ensuring that any comparison made is between comparable nodes of the contour tree. We also prove lower bounds showing that the $\Sigma$ pϵP |p| log |p| complexity is inherent to computing contour trees.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗