Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fracture network”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

From Points to Planes: A Workflow for Converting Three‐Dimensional Point Cloud Data Into Discrete Fracture Network Flow and Transport Models

We present the Point cLoud Algorithm for NEtwork Extraction of Discrete Fracture Networks (PLANE-DFN), a point cloud–based algorithm for automatic fracture network extraction designed to support discrete fracture network (DFN) modeling workflows. PLANE-DFN segments three-dimensional fracture planes from raw point cloud data using RANdom SAmple Consensus coupled with statistical outlier removal and density-based clustering to isolate individual fracture features. Each candidate plane is constrained against site-specific structural constraints based on strike and dip. After segmentation, each fracture is converted into a 2-D convex polygon suitable for meshing and simulation. The PLANE-DFN algorithm is validated by comparing geometric and flow and transport data against data from dfnWorks simulations with ensembles of plane-fit networks. We find that the flow and transport in plane-fit networks are comparable to dfnWorks-generated networks when realistic network geometry is maintained. The PLANE-DFN algorithm provides an automated and streamlined workflow to transform point clouds of data into DFN network geometry.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Characterizing the impact of finite matrix block size on conservative particle transport through three-dimensional fracture networks

Mass transfer of solutes between fractures and the surrounding rock matrix exerts a noticeable signature on the tail of travel time distributions. When the width of the matrix is assumed to be infinite and advective transport through the fracture is sufficiently fast, the tails of the travel time distributions exhibit a classically expected slope of ψ(t) ∝ t -3/2 . However, studies have yet to characterize how solute transfer between fractures via diffusion through finite matrix blocks influences the tail’s slope in three-dimensional fractured media. Here, in this study, we assess the impact of finite matrix block size on breakthrough curve shape at different spatio-temporal scales by con ducting particle tracking simulations in three-dimensional discrete fracture networks. We consider a variety of hydrodynamic and geostructural proper ties to determine their relative impact on the resulting travel time distributions. We observe that the impact of matrix diffusion through a finite block on travel time distributions is similar to that of an infinite matrix block when the fracture spacing is sufficiently large, matrix diffusion is relatively weak, or transport is considered at an early control plane distance. We observe that the converse of these conditions, results in deviations from the classical ψ(t) ∝ t -3/2 scaling. These results provide a first step toward developing a metric to assess when finite block size effects are expected to significantly influence transport.

58 GEOSCIENCES↗
DOE OSTI2025

FracML: A Machine Learning Based Tool to Quantify Reservoir Scale Fracture Network for CO2 Storage

Poster on “FRACML: A Machine Learning Based Tool to Quantify Reservoir Scale Fracture Network for CO2 Storage” for the CCUS 2025 conference held in Houston, Texas March 3-5, 2025. The accurate characterization of subsurface fracture networks is essential for the secure operation of carbon capture, utilization, and storage (CCUS) projects. A thorough understanding of the spatial distribution of subsurface faults and fractures is crucial for predicting CO2 plume evolution and minimizing risks such as potential leakage into overlying formations or induced seismicity. In this context, robust fracture network quantification plays a pivotal role in reservoir management, providing the data necessary to fine-tune operational parameters, and ensure the environmental and economic viability of CCUS projects. As part of the U.S. Department of Energy’s SMART (Science-informed Machine Learning for Accelerating Real-time Decisions in Subsurface Applications) initiative, we focused on the development and application of a machine learning-based tool (FRACML) designed to quantify and map fracture networks using real-world (non-synthetic) data from an active CO2 injection site. Our objective is to demonstrate the utility of this tool in improving operational efficiency and safety across CCUS sites.

artifical intelligence / machine learning (AI/ML)↗
DOE OSTI2025

Learning the factors controlling mineral dissolution in three-dimensional fracture networks: applications in geologic carbon sequestration

We perform a set of high-fidelity simulations of geochemical reactions within three-dimensional discrete fracture networks (DFN) and use various machine learning techniques to determine the primary factors controlling mineral dissolution. The DFN are partially filled with quartz that gradually dissolves until quasi-steady state conditions are reached. At this point, we measure the quartz remaining in each fracture within the domain as our primary quantity of interest. We observe that a primary sub-network of fractures exists, where the quartz has been fully dissolved out. This reduction in resistance to flow leads to increased flow channelization and reduced solute travel times. However, depending on the DFN topology and the rate of dissolution, we observe substantial variability in the volume of quartz remaining within fractures outside of the primary subnetwork. This variability indicates an interplay between the fracture network structure and geochemical reactions. We characterize the features controlling these processes by developing a machine learning framework to extract their relevant impact. Specifically, we use a combination of high-fidelity simulations with a graph-based approach to study geochemical reactive transport in a complex fracture network to determine the key features that control dissolution. We consider topological, geometric and hydrological features of the fracture network to predict the remaining quartz in quasi-steady state. We found that the dissolution reaction rate constant of quartz and the distance to the primary sub-network in the fracture network are the two most important features controlling the amount of quartz remaining. This study is a first step towards characterizing the parameters that control carbon mineralization using an approach with integrates computational physics and machine learning.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Fracture Network Quantification for Enhanced Reservoir Characterization in CO2 Storage Sites

Presentation slides on “Fracture Network Quantification for Enhanced Reservoir Characterization in CO2 Storage Sites” for the CCUS 2025 conference held in Houston, Texas March 3-5, 2025. Accurate reservoir characterization is fundamental for the safe and efficient operation of carbon capture, utilization, and storage (CCUS) projects. This process is not only critical during initial site selection but also plays a key role in continuously updating reservoir geomodels throughout injection operations to effectively manage subsurface storage sites. Our primary objective is to develop a comprehensive understanding of fracture networks within CO2 injection reservoirs, which are critical for predicting the behavior of injected fluids, ensuring storage integrity post-injection, and mitigating risks such as leakage or induced seismicity.

artifical intelligence / machine learning (AI/ML)↗
DOE OSTI2025

Fracture Network Quantification during CO2 Injection

This is the presentation prepared for the ARMA 2025 (59th US Rock Mechanics/Geomechanics Symposium) Conference held in Santa Fe, New Mexico, June 8-11, 2025. Accurate mapping and quantification of these networks are essential to ensure the integrity of CO2 storage reservoirs, understand and reduce potential leakage, and maintain long-term environmental safety. This study presents a novel machine learning-driven approach, integrated with geomechanical analysis, to quantify fracture networks and assess their spatial distribution during CO2 injection. This paper combines microseismic monitoring data with principles of hydraulic diffusivity and geomechanical analysis to characterize reservoir scale fracture network. The novelty of our approach lies in its capacity to assimilate time-dependent pressure data and microseismicity into a cohesive framework, which not only identifies microseismic triggering fronts but also tracks fracture distribution during active injection. Besides, leveraging image log data and analysis our approach also provides another angle of the insights to solidate the fracture networks understanding and geomechanical impacts. Key results from our study include the detection of over 100 distinct fracture clusters across the injection site, with fracture orientations strongly correlated with the prevailing in-situ stress field.

CO2 storage and sequestration↗
DOE OSTI2025

Fracture Network Quantification during CO2 Injection

This is the conference paper accompanying an oral presentation at the ARMA 2025 (59th US Rock Mechanics/Geomechanics Symposium) Conference held in Santa Fe, New Mexico, June 8-11, 2025. Accurate mapping and quantification of these networks are essential to ensure the integrity of CO2 storage reservoirs, understand and reduce potential leakage, and maintain long-term environmental safety. This study presents a novel machine learning-driven approach, integrated with geomechanical analysis, to quantify fracture networks and assess their spatial distribution during CO2 injection. This paper combines microseismic monitoring data with principles of hydraulic diffusivity and geomechanical analysis to characterize reservoir scale fracture network. The novelty of our approach lies in its capacity to assimilate time-dependent pressure data and microseismicity into a cohesive framework, which not only identifies microseismic triggering fronts but also tracks fracture distribution during active injection. Besides, leveraging image log data and analysis our approach also provides another angle of the insights to solidate the fracture networks understanding and geomechanical impacts. Key results from our study include the detection of over 100 distinct fracture clusters across the injection site, with fracture orientations strongly correlated with the prevailing in-situ stress field.

CO2 storage and sequestration↗
DOE OSTI2025

Discrete fracture network model benchmarks developed and applied in a DECOVALEX-2023 repository performance assessment study

This study presents newly developed benchmarks for modeling flow and transport within discrete fracture networks (DFNs) and useful methods for analyzing the results. The new benchmarks are designed to test modeling approaches for use in probabilistic performance assessment models of deep geologic repositories in fractured rock. The benchmarks simulate flow and transport through a 1 km 3 block of fractured rock. The first simulates migration of a short pulse of tracer through a simple network of four intersecting fractures. The second adds 1089 stochastically generated fractures. The third changes the pulse to a continuous point source. Evaluation of model performance relies on moment analysis and comparison of the results of different models. The expected nondimensional first moment of the conservative tracer for each benchmark is 1. The benchmarks were simulated by teams from Canada, Czechia, Germany, Korea, Sweden, Taiwan, and the United States as part of a DECOVALEX-2023 study (decovalex.org). The teams used various approaches, including explicit DFN modeling, DFN upscaling to an equivalent continuous porous medium (ECPM), and a combination of both methods. Transport mechanisms are modeled using either the advection-dispersion equation or particle tracking. Results demonstrate strong agreement among the models in breakthrough behavior up to the 75th percentile. Significant deviations in first moments and well-clustered outputs led to the identification of inaccuracies in several models. Such findings exemplify the benefit of exercising these benchmarks and using the presented methods to test DFN flow and transport models.

Benchmark↗
DOE OSTI2026

Effects of Dissolution Regimes on Flow Channelization and Solute Transport in 3D Fracture Networks: Insights From Graph‐Based Reactive Transport Modeling

We investigate how mineral dissolution reshapes flow pathways and solute transport in three‐dimensional discrete fracture networks using a computationally efficient graph‐based reactive transport model. The DFNs are inspired by field‐site observations of fractured carbonate and represent realistic connectivity and structural heterogeneity. Flow is simulated with the Reynolds equation, and dissolution follows first‐order kinetics with diffusive limitations captured through an effective mass‐transfer coefficient. By systematically varying two key dimensionless parameters, the effective Damköhler number (Da), governing reaction versus advection rates, and a transport parameter (Da), analogous to the Thiele modulus, distinct flow channelization regimes emerge: mildly channelized at low G, highly channelized at intermediate Da, and extreme wormhole formation at high Da and low G. Eulerian and Lagrangian analyses, including breakthrough curves, particle tortuosity, dispersivity, and flow channeling indicators quantitatively characterize the progression of dissolution‐driven network restructuring. Across all regimes, initial fracture heterogeneity persists. The results underscore how the interplay between this initial structure, advection, reaction, and diffusion critically shapes subsurface flow pathways, with implications for applications ranging from groundwater remediation to enhanced geothermal systems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Predicting Flow in Fracture Networks With Quantum Algorithms

Uncertainty quantification plays a crucial role in the modeling of subsurface flow. For instance, uncertainties in the properties of geologic fracture networks significantly impact flow, requiring numerous simulations to accurately estimate quantities of interest. However, each simulation is computationally expensive because it requires solving a large linear system to capture features that involve both small and large fractures. An example is in percolation, where the interaction of many small fractures (which cumulatively can have a large surface area) with the rock matrix must be modeled precisely. Quantum computing is an emerging tool with the potential to address this issue. Quantum algorithms offer a significant speedup in solving linear systems, achieving efficiencies that are challenging to match with classical approaches. These classical approaches include direct solvers, such as LU decomposition, and iterative methods, notably preconditioned conjugate gradient, commonly used in subsurface modeling to solve large sparse systems. However, applying quantum algorithms to geologic fracture flow requires careful attention to algorithmic and problem-specific constraints to fully realize this quantum advantage. In this work we describe a quantum algorithm for generalized Monte Carlo applications with a quadratic speedup over the classical approaches which can be combined with the quantum speedup, currently under investigation, for solving quantum linear systems for subsurface flow. We show that for quantum algorithms the computational cost of estimating a quantity of interest for a statistical ensemble of networks is roughly the same as that of a single realization, essentially implying that one can get uncertainty quantification for free.

58 GEOSCIENCES↗
DOE OSTI2026

Geologic stress modulates fluid mixing at fracture intersections

Fracture intersections are critical links that enable flow and transport in subsurface fracture networks, and their behavior strongly influences fluid mixing in a network. Although all subsurface fractures are subjected to geological stress, we lack a fundamental understanding of how fracture intersection geometry evolves under stress and how these changes influence fluid mixing. Here, we combine 3D printing, 3D X-ray tomographic imaging, and 3D pore-scale numerical simulations to reveal stress-induced changes in intersection geometry and their impact on mixing. Mixing is found to be strongly affected by partial closure of an intersection under stress. As an intersection closes, the void area for fluid flow and diffusion decreases leading to substantial deviations between conventional mixing models and full pore-scale modeling. To address this, we propose a modified mixing model that accounts for intersection deformation, which is essential for accurate modeling of solute transport and mixing through fracture networks.

15 GEOTHERMAL ENERGY↗
DOE OSTI2025

Mechanistic understanding of carbon mineralization in fracture systems using microfluidics

Carbon mineralization in mafic and ultramafic rocks presents an opportunity for permanent carbon storage in the Earth's subsurface. However, due to their lower permeability, pre-existing fracture networks are key for mineralization to occur. Therefore, to fully develop this technology, a mechanistic understanding of the mineralization behavior in fractures with the consideration of hydrodynamic components is required. We use high-pressure microfluidics to investigate key mechanisms influencing dissolution–precipitation in a fracture network. The experiments were conducted in micromodels made of natural rocks with a comb-shaped flow channel to mimic a fracture network. This enabled studying the effect of injection rate on coupled dissolution–precipitation in advection and diffusion-dominated flow paths. We used gypsum carbonation as an analog reaction to allow for realistic experimental time frames due to its rapid reaction kinetics. The experimental work is coupled with high-fidelity numerical simulations to enhance our understanding of the parameters affecting the mineralization reaction. Our results demonstrate the importance of flow rate on the rate and nature of the gypsum carbonation reaction revealing that higher flow rates enable deeper penetration of the mineral precipitation front into the dead-end channels. This is an important finding since for sustained mineralization in a fracture network, precipitation in dead-ends while still allowing for flowing fractures is critical. Detailed characterization of the precipitates showed that lower flow rates led to porous and loose precipitates in the form of aragonite while higher flow rates mimicked supersaturation behavior leading to the formation of calcite. The reactive transport simulations further demonstrated the significance of flow velocity in advection-dominated channels to influence the efficiency of carbon mineralization in diffusion-dominated channels, potentially clogging of dead-end channels. These findings highlight the need for coupling chemical, mechanical, and hydrodynamic processes to evaluate the nature and extent of carbon mineralization in fractured media critical for permanent storage in mafic and ultramafic formations. This research further highlights the need for more investigation in potential subsurface fracture generation techniques to aid carbon mineralization.

25 ENERGY STORAGE↗
DOE OSTI2026

Sequential Fracture Activation and Stress Evolution During EGS Stimulation at Utah FORGE Revealed by Waveform Cross‐Correlation

Mapping fracture networks in Enhanced Geothermal Systems (EGS) is essential for optimizing reservoir performance, yet complex fracture evolution during stimulation remains difficult to resolve. This study examines the evolution of microseismicity and fracture networks during stage 3 of the 2022 EGS stimulation at the Utah Frontier Observatory for Research in Geothermal Energy site. We map the fracture network represented by 20 clusters of seismic events identified by waveform similarities with cross-correlation. We characterize their geometric properties such as strike, dip, length, and width, and analyze the time evolution of activated fractures. The results reveal a systematic fracture evolution: early activation of pre-existing natural fractures, complex network development during peak injection, and continued activation of less favorably oriented fractures post-injection. Magnitude calibration using the Principal Component Analysis of cross-correlated waveforms improves relative amplitude measurements, refining estimations of the Gutenberg-Richter b-values with spatial variations in b-values suggesting stress re-distribution across the stimulated area. Analysis of the stress state of selected fractures further shows that fractures requiring higher excess pore pressure primarily activate at the end of injection and post-injection, highlighting stress transfer due to pore pressure as a dominant triggering mechanism. These findings provide insights into fracture propagation, stress evolution, and seismic hazard assessment in EGS reservoirs.

Asirifi, Richard [Texas A & M Univ., College Stati↗
DOE OSTI2026

Fracture Intersections under Stress: Laboratory Data and Code [Data set]

The connectivity of natural and induced fractures governs the injection and withdrawal of fluids from subsurface reservoirs. Connectivity depends on intersections that control how fluids mix and move through the entire system. Here, we present data sets from 3D X-ray microscopy measurements of simple fracture networks under stress. 3D printing was used to create prismatic blocks that formed fracture networks composed of 2 orthogonal fractures. The network orientation was either "x" or "+" relative to an applied vertical stress. 3D data sets were collected for normal loads of 25, 100 and 200 Newtons for samples with fracture surfaces with either correlated or uncorrelated asperity distributions. The file contains data from the 12 samples analyzed along with an example code used to extract the intersection geometry. Additional experimental details can be found in the manuscript "Geologic Stress Modulates Fluid Mixing at Fracture Intersections" (10.1038/s43247-026-03525-9)and supplemental information to appear in Communications Earth & Environment in 2026.

02 PETROLEUM↗
DOE OSTI2024

Complex Fluid‐Driven Fractures Caused by Crack‐Parallel Stress

Abstract Managing fluid‐driven fracture networks is crucial for subsurface resource utilization, yet the current understanding of the key controlling factors remains insufficient. While geologic discontinuities have been shown to significantly influence fracture network complexity, this study identifies another major contributor. We conducted a new set of experiments using a transparent true triaxial cell, which enabled video recording of the temporal evolution of fluid‐driven fracture paths. Using pseudo‐2D samples without macroscale structural discontinuities, we observed multiple occurrences of hydraulic fracture curving and branching under anisotropic boundary stresses. We proposed a theoretical model demonstrating that the stress parallel to the crack line in the solid matrix near the crack tip (i.e., the T ‐stress) accounts for the observed fracture curving behavior. This finding suggests that T ‐stress is an additional mechanism contributing to the complexity of fluid‐driven fracture networks in the subsurface, besides the geologic discontinuities.

58 GEOSCIENCES↗