Engineering Papers⌕ Search

Engineering topics

Cusini, Matteo

Publications and source records attributed to Cusini, Matteo.

Closing the Loop Between In Situ Stress Complexity and EGS Fracture Complexity (Quarterly Research Performance Progress Report: Q3Y3)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. To start, we have been working towards adjusting numerical parameters to reproduce DFIT tests with different fracture geometries. Then, we have focused modeling the DFIT conducted at well 58-32 (cycle 4 in zone 2). Specifically, we have carefully analyzed the influence of the aperture-contact stress relationship and investigated how two different laws affect model results. We have also introduced thermal effects in the model and rerun the simulations to estimate the effect of cooling on the DFIT response. Due to the small amount of fluid injected during a DFIT, thermal effects seem to be moderate.

15 GEOTHERMAL ENERGY↗

GEOS: A performance portable multi-physics simulation framework for subsurface applications

GEOS is a simulation framework focused on solving tightly coupled multi-physics problems with an initial emphasis on subsurface reservoir applications. Currently, GEOS supports capabilities for studying carbon sequestration, geothermal energy, hydrogen storage, and related subsurface applications. The unique aspect of GEOS that differentiates it from existing reservoir simulators is the ability to simulate tightly coupled compositional flow, poromechanics, fault slip, fracture propagation, and thermal effects, etc. Extensive documentation is available on the GEOS documentation pages (GEOS Documentation, 2024). Note that GEOS, as presented here, is a complete rewrite of the previous incarnation of the GEOS referred to in (Settgast et al., 2017).

58 GEOSCIENCES↗

Utah FORGE 2-2446: Closing the Loop Between In-Situ Stress Complexity and EGS Fracture Complexity - 2024 Annual Workshop Presentation

This is a presentation on Closing the Loop Between In-Situ Stress Complexity and EGS Fracture Complexity by Lawrence Livermore National Laboratory, presented by Matteo Cusini. The video discusses the combination of high-fidelity simulations and true-triaxial block fracturing tests to explore the intricate relationship between in-situ stress and hydraulic fracture patterns. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 14, 2024.

15 GEOTHERMAL ENERGY↗

FORGE STRESS annual report

The project's goal is to combine high-fidelity numerical models and true-triaxial block fracturing tests at high temperatures to understand the relationship between in situ stress, thermal effects, wellbore orientations and hydraulic fracture patterns. The numerical models are calibrated against field data, such as well pressures and microseismic data, and employed to estimate the in situ stress at the FORGE site. Laboratory experiments investigate the complex physics driving hydraulic fracture nucleation in EGS, enhancing understanding of the role of parameters like temperature, well orientation, and stress. Additionally, they are employed to validate some of the numerical tools used in the project. The project will have a significant impact by: (1) improving the characterization of the in-situ stress field at FORGE; (2) demonstrating the use of high-fidelity modeling tools for EGS; (3) providing a unique set of high-temperature hydraulic fracturing results to identify key components for in-situ stress estimation, validate current theories, and propose new ones; (4) offering a validated set of numerical tools within an open-source simulation framework, GEOS, that will be available to any future user.

15 GEOTHERMAL ENERGY↗

Simulation of Multiphase Flow and Poromechanical Effects Around Injection Wells in CO 2 Storage Sites

In geological CO 2 storage operations, wellbore deformations and leakage pathways formations can occur around injection and abandoned wells subjected to high rates and long-term CO 2 injection. To guide engineering design and prevent CO 2 leakage risks, a full understanding of the underlying physics and robust numerical models is necessary to evaluate the response of underground formations in the near wellbore region and in the reservoir. In this study, a multi-scale and multi-physics open-source simulator (GEOS) is used to simulate multiphase flow and poromechanical deformations over time in three dimensions. The governing equations for mechanical deformations of the rock body and multiphase compositional fluid flow within the rock matrix are solved with a fully coupled finite element and finite volume approach. The Drucker–Prager model with friction hardening is applied to simulate elastoplastic deformation and a multiphase fluid model with power-law correlations for relative permeability is used to model the migration of CO 2 plume, which are coupled with numerical implicit scheme. Simulation results are verified against multiple analytical solutions for multiphase flow and wellbore problems, thus demonstrating the accuracy of this advanced simulator. In two engineering applications, here we highlight the impact of elastoplastic deformation and coupled modeling for assessing induced displacements and stress perturbations, which are more pronounced in the near wellbore regions. This work focuses on short-term processes in the vicinity of injection wells where stress evolutions, rock deformations and multiphase compositional flow and transport are simulated jointly to ensure wellbore stability and prevent damage. This fully coupled geomechanical model can simulate multiphase flow and any associated poromechanical effects within the CO 2 storage site and in the surrounding formations. Such a large-scale, long-term, multi-physics simulation model is useful in many ways: it can guide operational decisions for CO 2 injection, assess the containment potential and risks of a site, and analyze the wellbore stability and integrity during and after CO 2 injection.

58 GEOSCIENCES↗

Phase Field Modeling of Near-Wellbore Hydraulic Fracture Nucleation and Propagation

The goal of Task 2 of the project titled "Closing the loop between in situ stress complexity and EGS fracture complexity" is to model fracture nucleation and propagation in the near wellbore region and, in combination with the experimental study conducted as part of Task 4, investigate the relationship between the in situ stress and the near-wellbore fracture patterns. To this end, We have identified phase field as a promising modeling approach to model the nucleation and propagation process in the near wellbore region. In phase field, fractures are represented by a diffuse variable, i.e the damage. Compared to approaches that model fractures as sharp interfaces, phase field does not require any re-meshing or element insertion algorithm, which makes the treatment of complex fracture geometries straightforward, even with simple meshes.

15 GEOTHERMAL ENERGY↗

Quarterly Research Performance Progress Report (Q8)

As part of Task 1, we have started by testing our modeling capabilities by reproducing isothermal DFIT simulations presented in the literature. Once satisfied with the results we have started by targeting the modeling of the DFITs at conducted at well 58-32. We have a identified a specific test (cycle 4 in zone 2) as the most interesting to be model with GEOS hydraulic fracturing module. Thus, we have first produced results with an isothermal model and adjusted model parameters to get a satisfying match with field pressure data. The, we have added thermal effects and compared the modeling results with and without thermal effects to estimate how thermal effects may influence test interpretation. Models seem to suggest that, for small volumes of fluid, thermal effects are moderate. In Task 2, we have adapted GEOS phase-field formulation to be able to simulate near-wellbore hydraulic fracture nucleation and propagation. We have devised a novel formulation that, compared to other existing ones, incorporates rock strengths. We have submitted a journal publication about our work. We are currently employing this phase-field formulation to model the experiments taking place at U Pitt and help us understand the effect of various parameters. In Task 3, we have built a model of the region surround well 16A and have started modeling stage 3 stimulation because of its simpler planar geometry. After calibrating simulation parameters using known analytical solutions, we have simulated the stage 3 stimulation using our isothermal hydraulic fracturing module, varying the permeability field, the stress conditions including different physics to get a better understanding of the numerical challenges and of the effects of varying these parameters on the simulation results. In Task 4 laboratory experimentation, a set of specialized drilling and injection tools has been customized and constructed to accommodate an inclined well with an orientation of up to 30 degrees relative to material anisotropy or principal stress axes. These inclined samples have also undergone thermal stress and hydraulic fracturing at a temperature of 190 degrees Celsius. Furthermore, both vertical and deviated sampling testing setups enable an extended analysis of post-peak pressure behaviors, facilitating post-test pressure analyses such as the G-function, step rate, and fracture reopening measurements. Thus, the key components of in-situ stress estimation can be extracted and validated through our experiment, providing a solid foundation for validating existing in-situ stress estimation theories or proposing new ones. Simultaneously, we are integrating computer vision techniques with traditional experimental fracture observation methods such as multi-overcore/slicing and water-penetration fracture observation. This combination will prove beneficial in populating the hydraulic fracture patterns database, generated under challenging EGS conditions. This approach aims to deepen our understanding of the complexities in EGS reservoirs and pave the way for future data-driven investigations. Additionally, PITT has also equipped the ELE International compression machine, which is now prepared for conducting indirect tensile and fracture toughness tests. These tests will aid in characterizing how rock fabrics influence the resulting fracture patterns. Additionally, we have completed the required personnel training and gained access to Scanning Electron Microscopy (SEM) and Energy Dispersive Spectroscopy (EDS) for conducting more detailed characterization and analysis of rock fabrics, as well as the examination of thermal and hydraulically induced fracture patterns. Thus, the PITT team has effectively demonstrated the capabilities of our experimental apparatuses in exploring the thermal effects, well deviation angles, material anisotropy, and operational choices (such as circulation rate, injection fluid viscosity, and injection rate) and their impact on pressure responses and fracture trajectories under the Utah FORGE conditions.

15 GEOTHERMAL ENERGY↗

Go/No-Go Decision Point 1 Report

To satisfy the G/NG point 1 we have employed GEOS hydraulic fracturing module to simulate one of the Diagnostic Fracture Injection Tests conducted at well 58-32. The hydraulic fracturing solver relies on a finite element discretization of the mechanics fully coupled to a finite volume discretization of the fluid flow. Fractures are explicitly modeled as lower dimensional manifolds and discretized by 2D elements sitting at the boundary of cell elements. We have calibrated our model to reproduce the DFIT corresponding to Zone 2 Cycle 4 of the injection tests that took place at well 58-32.

15 GEOTHERMAL ENERGY↗

A phase-field model for hydraulic fracture nucleation and propagation in porous media

Many geo-engineering applications, for example, enhanced geothermal systems, rely on hydraulic fracturing to enhance the permeability of natural formations and allow for sufficient fluid circulation. Over the past few decades, the phase-field method has grown in popularity as a valid approach to modeling hydraulic fracturing because of the ease of handling complex fracture propagation geometries. However, existing phase-field methods cannot appropriately capture nucleation of hydraulic fractures because their formulations are solely energy-based and do not explicitly take into account the strength of the material. Thus, in this work, we propose a novel phase-field formulation for hydraulic fracturing with the main goal of modeling fracture nucleation in porous media, for example, rocks. Built on the variational formulation of previous phase-field methods, the proposed model incorporates the material strength envelope for hydraulic fracture nucleation through two important steps: (i) an external driving force term, included in the damage evolution equation, that accounts for the material strength; (ii) a properly designed damage function that defines the fluid pressure contribution on the crack driving force. In conclusion, the comparison of numerical results for two-dimensional test cases with existing analytical solutions demonstrates that the proposed phase-field model can accurately model both nucleation and propagation of hydraulic fractures. Additionally, we present the simulation of hydraulic fracturing in a three-dimensional domain with various stress conditions to demonstrate the applicability of the method to realistic scenarios.

15 GEOTHERMAL ENERGY↗

Utah FORGE 2-2446: Closing the Loop Between In-situ Stress Complexity and EGS Fracture Complexity - Workshop Presentation

This is a presentation on the Closing the Loop Between In-situ Stress Complexity and EGS Fracture Complexity project by Lawrence Livermore National Laboratory, presented by Dr. Matteo Cusini. The project's objective was to employ a combination of high-fidelity simulations and true-triaxial block fracturing tests at high temperature to explore the intricate relationship between in-situ stress and hydraulic fracture patterns and better characterize the in-situ stress at Utah FORGE. This presentation was featured in the Utah FORGE R&D Annual Workshop on September 7, 2023. The workshop provided a valuable opportunity to explore the progress made in each of the 17 Research and Development projects funded under Solicitation 2020-1 which aim to enhance our understanding of the crucial factors influencing the development of Enhanced Geothermal Systems (EGS) reservoirs and resources.

15 GEOTHERMAL ENERGY↗

Closing the loop between in-situ stress complexity and EGS fracture complexity

The goal of the project is to employ a combination of high-fidelity simulations and true-triaxial block fracturing tests at high temperature to explore the intricate relationship between in-situ stress and hydraulic fracture patterns and better characterize the in-situ stress at Utah FORGE. Laboratory experiments are used to validate numerical models and better understand the impact of various parameters (e.g., temperature, well orientation, stress, etc) on the hydraulic fracture nucleation process in EGS. The project will have a significant impact as it will (1) improve the characterization of the in-situ stress field at FORGE; (2) demonstrate how high-fidelity modeling tools can be employed to analyze data for EGS; (3) provide a unique set of high-temperature hydraulic fracturing experimental results for EGS to help identifying the most important components for in-situ stress estimation and to validate the current in-situ stress estimation theory and establish a strong foundation for proposing a new one; (4) provide a validated set of numerical tools in an open source simulation framework (GEOS) that can be used to model EGS.

15 GEOTHERMAL ENERGY↗

A multi-resolution approach to hydraulic fracture simulation

Abstract We present a multi-resolution approach for constructing model-based simulations of hydraulic fracturing, wherein flow through porous media is coupled with fluid-driven fracture. The approach consists of a hybrid scheme that couples a discrete crack representation in a global domain to a phase-field representation in a local subdomain near the crack tip. The multi-resolution approach addresses issues such as the computational expense of accurate hydraulic fracture simulations and the difficulties associated with reconstructing crack apertures from diffuse fracture representations. In the global domain, a coupled system of equations for displacements and pressures is considered. The crack geometry is assumed to be fixed and the displacement field is enriched with discontinuous functions. Around the crack tips in the local subdomains, phase-field sub-problems are instantiated on the fly to propagate fractures in arbitrary, mesh independent directions. The governing equations and fields in the global and local domains are approximated using a combination of finite-volume and finite element discretizations. The efficacy of the method is illustrated through various benchmark problems in hydraulic fracturing, as well as a new study of fluid-driven crack growth around a stiff inclusion.

58 GEOSCIENCES↗

Simulation of coupled multiphase flow and geomechanics in porous media with embedded discrete fractures

In fractured natural formations, the equations governing fluid flow and geomechanics are strongly coupled. Hydrodynamical properties depend on the mechanical configuration, and they are therefore difficult to accurately resolve using uncoupled methods. In recent years, significant research has focused on discretization strategies for these coupled systems, particularly in the presence of complicated fracture network geometries. In this work, we explore a finite-volume discretization for the multiphase flow equations coupled with a finite-element scheme for the mechanical equations. Fractures are treated as lower dimensional surfaces embedded in a background grid. Interactions are captured using the embedded discrete fracture model (EDFM) and the embedded finite element method (EFEM) for the flow and the mechanics, respectively. This nonconforming approach significantly alleviates meshing challenges. EDFM considers fractures as lower dimension finite volumes that exchange fluxes with the rock matrix cells. The EFEM method provides, instead, a local enrichment of the finite-element space inside each matrix cell cut by a fracture element. Both the use of piecewise constant and piecewise linear enrichments are investigated. They are also compared to an extended finite element approach. One key advantage of EFEM is the element-based nature of the enrichment, which reduces the geometric complexity of the implementation and leads to linear systems with advantageous properties. Synthetic numerical tests are presented to study the convergence and accuracy of the proposed method. It is also applied to a realistic scenario, involving a heterogeneous reservoir with a complex fracture distribution, to demonstrate its relevance for field applications.

58 GEOSCIENCES↗