Controlling fracture propagation using deep reinforcement learning
Not Available
SEARCH · Engineering Papers
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.
Not Available
The presence of fractures in caprocks can pose increased risks in subsurface energy systems and processes like CO 2 sequestration by introducing high-permeability leakage paths. Fracture apertures and permeability can be altered through mineral dissolution and precipitation reactions, but the reactive evolution of fractures is not well understood. In fractures, minerals that are otherwise inaccessible to reactive fluids can become exposed, resulting in mineral reactions unpredicted by bulk formation data. This work seeks to understand the relationship between mineralogy and fracture formation to enhance our understanding of reactive fracture evolution and CO 2 leakage potential. Here, the mineral compositions of mechanically induced fracture surfaces in samples of the Mancos and Marcellus shales have been quantified and compared to those of the near-fracture matrices using imaging and bulk X-ray diffraction (XRD) data. In the Mancos shale, the concentrations of clay minerals are enhanced along fracture surfaces with respect to the bulk, and the fracture is most likely to form at kaolinite–kaolinite interfaces. Further, evaluation of the mineralogical spatial variability through cross-correlation analysis of the surrounding matrix in images of samples cut perpendicular to the fracture shows that clay is 16.7 times more likely to be present than carbonate minerals near the fracture surface. The high correlation persists roughly 200 μm into the surrounding matrix for the Mancos sample and implies that the fracture formed within a defined clay-rich lithofacies.
Explore the source record for details and available documents.
This work explores the application of the recently developed Koopman operator approach for model identification and feedback control of a hydraulic fracturing process. Controlling fracture propagation and proppant transport with precision is a challenge due in large part to the difficulty of constructing approximate models that accurately capture the characteristic moving boundary and highly-coupled dynamics exhibited by the process. Koopman operator theory is particularly attractive here as it offers a way to explicitly construct linear representations for even highly nonlinear dynamics. The method is data-driven and relies on lifting the states to an infinite-dimensional space of functions called observables where the dynamics are governed by a linear Koopman operator. Here this work considers two problems: (a) fracture geometry control, and (b) proppant concentration control. In both cases, an approximate linear model of the corresponding dynamics is constructed and used to design a model predictive controller (MPC). The manuscript shows that in the case of highly nonlinear dynamics, as observed in the proppant concentration, use of canonical functions in the observable basis fails. In such cases, a priori system knowledge can be leveraged to choose the required basis. The numerical experiments demonstrate that the Koopman linear model shows excellent agreement with the real system and successfully achieves the desired target values maximizing the oil and gas productivity. Additionally, due to its linear structure, the Koopman models allow convex MPC formulations that avoid any issues associated with nonlinear optimization.
Fluid flow through fractures is intimately linked to the fracture surfaces that define the void geometry through which fluids flow. Thus, an understanding of what controls fracture surface roughness is essential to the development of models for predicting fluid transport through fractured rock. The difficulty in predicting surface roughness arises from the complexity of rock which is inherently heterogeneous and nonuniform in composition, fabric, and structural components, even when samples are acquired from the same rock mass. Here, a benchmarked-simulation approach motivated from geo-architected 3D printed synthetic gypsum rocks is used to provide insight into the competing contributions from fabric and layering on fracture roughness formation. Simulation results from a discrete element model (Particle Flow Code, Itasca Consulting Group, Inc.) clearly indicate that the relative orientation between mineral layers and in-layer mineral fabric, and the variability in mineral bonding strengths determine whether anisotropic corrugated surfaces or isotropic surfaces are formed. Weak mineral layers oriented perpendicular to the applied load resulted in strong roughness anisotropy. Peak failure loads were found to vary up to 30% depending on the strength of the mineral fabric at the location of fracture initiation, which provides insight into the observed high variability in strength values of natural rock. The uniqueness of induced fracture roughness and peak failure load is intimately linked to layering, mineral fabric, and their distribution in the rock. These findings have important implications for any architected material fabricated through serial printing of layers with local compositional heterogeneity.
Fluid flow through fractured rock systems is governed in large part by the distribution, interconnectivity, and size of fracture apertures. In-situ stress is one of the primary factors controlling fracture aperture, and one that is altered significantly during high-pressure fluid injections or extractions. Interactions between stress, pore pressure, aperture, and fluid flow can result in complex and evolving poroelastic behavior with significant implications regarding the predictability and risk of developing and managing deep subsurface reservoirs (geothermal, fossil energy, and geologic carbon sequestration). In saturated crystalline rocks, bulk electrical conductivity is highly sensitive to porosity, and therefore to fracture aperture size and distribution. We demonstrate the use of time-lapse 3D electrical resistivity tomography for remotely monitoring stress induced changes in aperture distribution during high pressure injections into a fractured rock system at a scale of tens of meters. Results reveal a complex and continuously evolving stress field involving aperture dilations in the natural fracture system and aperture contractions in adjacent zones of shadow stress. Results provide information about the spatiotemporal changes in the system behavior and point to the potential of electrical imaging for autonomously and remotely monitoring evolving stress conditions by proxy through changes in bulk electrical conductivity.
Approximately half of global energy consumption is used for heating and cooling. Because fossil fuels are used to meet most of this demand, heating and cooling of buildings produces a large portion of global greenhouse gas (GHG) emissions. The proposed project seeks to demonstrate the feasibility of a seasonal thermal energy storage concept called fracture thermal energy storage (FTES), which has the potential to dramatically lower energy requirements for heating and cooling and improve the resilience of building energy systems. FTES is a technique for building a highly efficient heat exchanger by creating a carefully designed set of fractures in the ground below a building. This heat exchanger allows thermal energy to be stored over seasonal timeframes—for example, hot thermal energy that is easily captured in the summer or energy from waste heat sources—to meet heating needs during colder months. FTES offers a solution with a relatively small footprint and lower cost compared to currently operating aquifer thermal energy storage (ATES) systems and could be used in the many areas without aquifers suitable for ATES. If commercialized, FTES systems could exceed the 2.5 TWh of energy storage per year from the more than 2,800 ATES systems in operation worldwide, which range from 0.1 MW to 30 MW. The savings in CO2 emissions per year are also expected to match or exceed individual ATES projects, with the largest, a 30 MW system used to provide heating and cooling to the University of Technology in Eindhoven in the Netherlands, estimated to save 13,000 tons of CO2 emissions per year. FTES utilizes mature drilling and fracturing technology and therefore has the potential to be rapidly commercialized once demonstrated. The ability to construct and establish flow through an FTES heat exchanger has been demonstrated by a previous project and the potential for efficient, large-capacity energy storage has been shown using numerical models. However, no experimental validation of these numerical estimates of thermal energy storage has been made. The first crucial need to advance FTES technology is to identify the sensitivity of key metrics such as thermal energy storage and production rates, capacities, and efficiencies to design parameters such as the number of fractures, depth/temperature of fractures, size of fractures, and circulation rates. The second crucial need is experimental testing of achievable thermal performance with optimized system design parameters. The proposed scope of work seeks to systematically address these two critical needs through a highly complementary international collaboration spanning theory, laboratory, and mesoscale field evaluation. The proposed work plan calls for using dimensional analysis and existing state-of-the-art numerical simulators to design carefully scaled laboratory and 10-meter-scale field tests of the thermal efficiency of FTES. The existing advanced laboratory and intermediate-scale field testbeds that will be used for this project will allow for detailed monitoring of the system performance during the test and of how the performance changes across time and length scales. These results will determine the feasibility of full- scale FTES systems. If the thermal performance is consistent with model predictions, the results will provide a strong economic justification for rapid commercialization of FTES technology in a wide range of geographical areas.
The DEMO-FTES project investigated the feasibility of Fracture Thermal Energy Storage (FTES) as a seasonal energy storage solution in crystalline rock formations. FTES leverages hydraulically induced fractures to exchange heat between circulating fluids and the surrounding rock mass, enabling long-term thermal energy retention due to the high specific heat and low thermal conductivity of rock. This approach has the potential to reduce heating and cooling energy demands and enhance building energy resilience. The project combined dimensional analysis, numerical modeling, laboratory experiments, and meso-scale field tests to evaluate FTES performance and advance its technology readiness level from 3 to 5. Scaling analysis identified key dimensionless parameters governing heat transfer and fluid flow, ensuring laboratory and field tests were representative of larger-scale systems. Numerical simulations using TOUGH and iTOUGH2 frameworks supported experiment design and interpretation, modeling fracture geometry, thermal-hydraulic behavior, and thermo-mechanical coupling. Laboratory tests at EPFL involved creating single and multiple fractures in 25 cm cubic samples of Gabbro and Granite under true triaxial stress.
The DEMO-FTES project sought to demonstrate the thermal efficiency of fracture thermal energy storage (FTES) through numerical simulations, laboratory and meso-scale field tests. A detailed dimensional and scaling analysis was performed to identify key parameters and how they can be most effectively scaled to the laboratory and decameter scale. Numerical models were developed and used for three purposes: 1. Before field experiments, numerical modelling can be used to estimate fracture properties based on previous data from the EGS Collab experiment and then predict thermal hydrological behaviors of the fracture system with hot water injection/withdrawal, therefore, help to design the experiments (e.g., to decide the duration of the cycles based on the flow rate the pump can provide, and the estimated fracture properties); 2. After the field experiment, to estimate the system properties during the experiment (as the size and shape of a fracture could change over time), and understand why system performance is different than what has been predicted, i.e., to help understand the meso-scale test; and, 3. To model the lab experiments and estimate fracture properties and storage efficiency. Ultimately, the experiment and numerical models could shed light on the processes and uncertainty happening during fracture activation and help understand the scaling between lab and field test, and finally, the design and optimization of potential fracture thermal energy storage systems.
Abstract not provided.
Objectives/Scope: StimuFrac (US Patents 9,873,828 B2 and 9,447,315 B2), a CO2-reactive polymer aqueous solution [polyallylamine (PAA) 1wt% in water] combined with CO2, can be used as a less water-intensive fracturing fluid for enhanced geothermal systems (EGS). Our previous results show that in hot dry rock (HDR), PAA/CO2 fracturing fluids outperformed other fluids such as water, CO2, and CO2/water in generating large fractures with less fluid consumed. The objective of this work is to study the effect of initial water saturation of rock on the performance of StimuFrac fluid in ½ foot cubic rock samples and under representative EGS pressure/temperature conditions using cyclic and constant flow rate injection strategies. The fracturing results are compared with results using different fracturing fluids in terms of controlling fracture propagation rates, fracture hydraulic conductivity, breakdown pressures and volumes of fluids required. Methods/Procedures/Process: In all tests, water was initially injected into the rock to increase the water saturation before the fracturing processes to simulate actual geothermal reservoir conditions. For the cyclic injection, one complete cycle consisted of (1) a PAA slug (or water slug) injection followed by (2) CO2 injection to initiate the fracture. In the second step of the first cycle, the pressure of CO2 is increased until a maximum pressure is reached (fracture is initiated at this moment), and then the injection of CO2 is allowed for another 30 seconds to propagate the fracture. Then, the two-step cycle of PAA followed by CO2 injection (up to 2-4 mL/min) was continued. For the constant flow rate injection strategies, the initial water saturation was increased by injecting water at 1000 psi and 200°C for three days. After that, an initial slug of water, CO2, or PAA was injected and then fracturing was initiated and propagated by injecting CO2 at a constant flow rate. Applications/Significance/Novelty: The results of this study suggest that water saturation, especially near the wellbore region, will significantly affect the fracturing fluid transmission into the rock porous media by changing the relative permeability of CO2 or water, thus affecting the fracture initiation and propagation. In this study, fracturing with cyclic injection or constant flow rate injection methods were performed using three different kinds of fluids systems. These fluids are water, CO2, or CO2 with PAA. Splitting the rock samples in half after fracturing reveals that the fracture propagation is significantly limited under these high water saturated conditions compared to dry initial conditions: The fractures propagate less than 1/3 length of the distance from the wellbore to rock surface, and in some cases no fracture is generated. This may be caused by the fact that leak-off is dominating the fracturing process and the injected fluid flow rate is not high enough to overcome the leak-off even under high flow rate injection conditions. Additionally, CO2 could be leaking off into the wellbore annulus and this may be making it more difficult to generate pressure gradients away from the near-wellbore region.
Abstract Fractures control fluid flow, mass transport, and heat transfer in a geothermal reservoir. This makes accurate characterization of fracture networks a prerequisite for optimal design and control of a reservoir's exploitation. We develop a deep‐learning procedure to identify fracture locations via interpretation of temporally and spatially continuous downhole temperature measurements. A long short‐term memory fully convolutional network (LSTM‐FCN) is used both to capture long‐term dependencies in sequential temperature data and to distill local features around fractures. A wellbore and fractured‐reservoir thermal model is established to generate temperature data for network training. The trained LSTM‐FCN exhibits a unique ability to detect multiple fractures intersecting a borehole. We use the LSTM‐FCN algorithm to evaluate the effectiveness of different‐stage wellbore temperature measurements on fracture detection in a complex fractured system. Our experiments reveal that the use of various‐stage temperature information as an input feature set improves the robustness of fracture detection to noise interference. This study indicates the practical feasibility of obtaining accurate fracture‐network reconstructions from temperature signals, at reasonable computational cost.
Not Available
Brittle materials, such as oxide glasses, are usually very sensitive to flaws, giving rise to a macroscopic fracture strength that is much lower than that predicted by theory. The same applies to metallic glasses (MGs), with the important difference that these glasses can exhibit certain plastic strain prior to catastrophic failure. Here we consider the strongest metallic alloy known, a ternary Co 55 Ta10B 35 MG. In this work, we show that this macroscopically brittle glass is flaw-insensitive at the micrometer scale. This discovery emerges when testing pre-cracked specimens with self-similar geometries, where the fracture stress does not decrease with increasing pre-crack size. The fracture toughness of this ultra-strong glassy alloy is further shown to increase with increasing sample size. Both these findings deviate from our classical understanding of fracture mechanics, and are attributed to a transition from toughness-controlled to strength-controlled fracture below a critical sample size.
Fracture spatial patterns can strongly affect fluid flow in the subsurface. Proximity and distribution of fractures control reservoir flow behavior over various length scales. In many studies, however, simplified geometrical patterns are generated for fractures which may lead to unrealistic subsurface models. Here we introduce a new method for characterization and modeling of fracture spatial patterns based on outcrop observations. We use Ripley's K-function to characterize the arrangement of fracture barycenters and intersection points over various length scales. In addition, we use semivariograms to quantify spatial correlation in fracture intensity maps. Using this information, we develop a stochastic algorithm that generates two-dimensional fracture network realizations with spatial properties similar to those of a real fracture network measured in the field. Numerical simulation models indicate that the generated fracture realizations exhibit similar flow behaviors as that of the original fracture network. Furthermore, such modeling tools expand and improve our capability in building representative fracture models and in quantification of uncertainty in naturally and hydraulically fractured reservoirs.
In support of the U.S. DOE SubTER Crosscut initiative, we established a field test facility in a deep mine and designed and carried out in situ hydraulic fracturing experiments relevant to enhanced geothermal systems (EGS) in crystalline rock to characterize the stress field, understand the effects of rock fabric on fracturing, and gain experience in monitoring using geophysical methods. The project also included pre- and post-fracturing simulation and analysis, and laboratory measurements and experiments. The kISMET (permeability (k) and Induced Seismicity Management for Energy Technologies) site was established in the West Access Drift of the Sanford Underground Research Facility (SURF) 4757 ft (1450 m) below ground (on the 4850 ft level (4850L)) in phyllite of the Precambrian Poorman Formation. We drilled and continuously cored five near-vertical boreholes in a line on 3 m (10 ft) spacing, deviating the two outermost boreholes slightly to create a five-spot pattern around the test borehole centered in the test volume 40 m below the drift invert (floor) at a total depth of ~1490 m (4890 ft). Laboratory measurements of core from the center test borehole showed P-wave velocity heterogeneity along each core indicating strong, fine-scale (~1 cm or smaller) changes in the mechanical properties of the rock. Field measurements of the stress field by hydraulic fracturing showed that the minimum horizontal stress at the kISMET site averages 21.7 MPa (3146 psi) trending approximately N-S (356 degrees azimuth) and plunging slightly NNW at 12°. The vertical and horizontal maximum stresses are similar in magnitude at 42-44 MPa (6090-6380 psi) for the depths of testing, which averaged approximately 1530 m (5030 ft). Hydraulic fractures were remarkably uniform suggesting core-scale and larger rock fabric did not play a role in controlling fracture orientation. Analytical solutions suggest that the fracture radius of the large fracture (stimulation test) was more than 6 m (20 ft), depending on the unknown amount of leak-off.