SEARCH · Engineering Papers

Results for “fluid injection”

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 253 records · Page 14

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Inorganic Alterations in Unconventional Shale Reservoirs: Importance of Additive and Base Fluid Chemistry

The effective development of unconventional petroleum systems requires the use of significant water resources. In an effort to reduce the consumption of freshwater resources for hydraulic fracturing, highly saline produced waters are increasingly recycled for use as a base fluid. However, there are significant knowledge gaps regarding potential water–rock interactions resulting from the introduction of produced waters and associated additives into shale reservoirs such as formation and deposition of mineral scale, which can negatively affect hydrocarbon production through wellbore restriction and damage to hydraulically generated fractures. To assess the impacts of field stimulation practices in the subsurface, a series of laboratory experiments were completed using (a) three distinct sedimentary rock formations of the Midland Basin (Texas, USA) and (b) additives with two different base fluids: municipal fresh water and clean brine. The experimental approach used relevant injection sequences and mixing ratios in specialized reactors for 3 weeks. Static pressurized experiments and nonpressurized time-resolved experiments were undertaken. The resulting solids and liquids were analyzed by using a variety of laboratory- and synchrotron-based techniques. The use of an acid spearhead (15% HCl) resulted in texturing of both clay-rich and calcareous shales, which can temporarily enhance porosity but subsequently result in mineral scale deposition. The primary matrix scale was Fe(III)-bearing phases, which occurred in all experiments regardless of base fluid chemistry. Additionally, strontium sulfate (SrSO 4 ) precipitated on shale surfaces when clean brines were used. It was concluded that clean brine was the main source of Sr 2+ species, while persulfate breaker degradation and oxidation of pyrite were the sources of SO 4 2– . Sulfate scaling was more pronounced in clay-rich shales, suggesting that Sr sorption is important for promoting celestite formation. This work demonstrates that mineral scale deposition is a complex phenomenon, whereby the type and proportions of various mineral phases are determined from reservoir alteration processes and coprecipitation of constituents from injection fluids. In conclusion, the experimental results shown here should be considered when evaluating different base fluids and additives in order to mitigate mineral precipitation in unconventional shale reservoirs, which could result in reservoir degradation.

Jew, Adam D. [SLAC National Accelerator Laboratory↗

NASA NTRS1993

Physiological water model development

The water of the human body can be categorized as existing in two main compartments: intracellular water and extracellular water. The intracellular water consists of all the water within the cells and constitutes over half of the total body water. Since red blood cells are surrounded by plasma, and all other cells are surrounded by interstitial fluid, the intracellular compartment has been subdivided to represent these two cell types. The extracellular water, which includes all of the fluid outside of the cells, can be further subdivided into compartments which represent the interstitial fluid, circulating blood plasma, lymph, and transcellular water. The interstitial fluid surrounds cells outside of the vascular system whereas plasma is contained within the blood vessels. Avascular tissues such as dense connective tissue and cartilage contain interstitial water which slowly equilibrates with tracers used to determine extracellular fluid volume. For this reason, additional compartments are sometimes used to represent these avascular tissues. The average size of each compartment, in terms of percent body weight, has been determined for adult males and females. These compartments and the forces which cause flow between them are presented. The kidneys, a main compartment, receive about 25 percent of the cardiac output and filters out a fluid similar to plasma. The composition of this filtered fluid changes as it flows through the kidney tubules since compounds are continually being secreted and reabsorbed. Through this mechanism, the kidneys eliminate wastes while conserving body water, electrolytes, and metabolites. Since sodium accounts for over 90 percent of the cations in the extracellular fluid, and the number of cations is balanced by the number of anions, considering the renal handling sodium and water only should sufficiently describe the relationship between the plasma compartment and kidneys. A kidney function model is presented which has been adapted from a previous model of normal renal function in man. To test the validity of the proposed kidney model, results predicted by the model will be compared to actual data involving injected or ingested fluids and subsequent urine flow rates. Comparison of the model simulation to actual data following the ingestion of 1 liter of water is shown. The model simulation is also shown with actual data following the intravenous infusion of hypertonic saline.

Doty, Susan↗

DOE OSTI2021

Numerical Modeling of CO2 Sequestration within a Five-Spot Well Pattern in the Morrow B Sandstone of the Farnsworth Hydrocarbon Field: Comparison of the TOUGHREACT, STOMP-EOR, and GEM Simulators

The objectives of this study were (1) to assess the fate and impact of CO2 injected into the Morrow B Sandstone in the Farnsworth Unit (FWU) through numerical non-isothermal reactive transport modeling, and (2) to compare the performance of three major reactive solute transport simulators, TOUGHREACT, STOMP-EOR, and GEM, under the same input conditions. The models were based on a quarter of a five-spot well pattern where CO2 was injected on a water-alternating-gas schedule for the first 25 years of the 1000 year simulation. The reservoir pore fluid consisted of water with or without petroleum. The results of the models have numerous broad similarities, such as the pattern of reservoir cooling caused by the injected fluids, a large initial pH drop followed by gradual pH neutralization, the long-term persistence of an immiscible CO2 gas phase, the continuous dissolution of calcite, very small decreases in porosity, and the increasing importance over time of carbonate mineral CO2 sequestration. The models differed in their predicted fluid pressure evolutions; amounts of mineral precipitation and dissolution; and distribution of CO2 among immiscible gas, petroleum, formation water, and carbonate minerals. The results of the study show the usefulness of numerical simulations in identifying broad patterns of behavior associated with CO2 injection, but also point to significant uncertainties in the numerical values of many model output parameters.

STOMP↗

DOE OSTI2025

Induced seismicity and geothermal energy production in the Salton Sea Geothermal Field, California

We analyze the relationship between geothermal energy production and seismic hazards in the Salton Sea Geothermal Field (SSGF) between 1972 and 2022. A clear increase in seismic activity accompanies geothermal energy production and is greatest to the east of the Brawley fault, where the amount of injection exceeds the amount of production. We estimate that, whereas there was a 2% chance of a M6.07 earthquake nucleating inside the SSGF within a fifty-year period prior to energy production (pre-production), there is an equal probability that a M6.70 earthquake nucleates there in a fifty-year period during which historically average energy production conditions persist (syn-production). Similarly, we estimate a 2% chance of a M7.15 and M6.92 occurring in the broader Brawley Seismic Zone (excluding the SSGF) within a fifty-year period during the pre- and syn-production eras, respectively. We find that linear regression models fail to reliably forecast the background seismicity rate as a function of fluid production and injection.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

Proxy-based Bayesian inversion of strain tensor data measured during well tests

Recent instrument developments have made it possible to measure the strain tensor caused by injecting or pumping fluid from aquifers or reservoirs, but the full value of these data is limited because the long runtimes of poroelastic forward models makes it impractical to use many inversion schemes. This limits the interpretation of strain data for managing the recovery of resources or storage of wastes in the subsurface. This paper describes a method of inverting deformation data using a poroelastic numerical simulator so the results can be used to manage reservoirs or aquifers. We developed a workflow designed to reduce the number of simulations sufficiently to make it feasible to use DREAMzs, an advanced Bayesian inversion method that translates the uncertainties from different sources into unbiased posterior parameter distributions and uncertainty envelopes around the field data. Using a KNN proxy model for the poroelastic simulator is key to reducing the overall computations, and the workflow includes a strategy for ensuring the proxy model results converge on the results from the simulator. The workflow is tested using an idealized example that verifies the ability to correctly identify parameters and characterize noise used to perturb the data. Field data from an injection test at an oil reservoir near Tulsa, Oklahoma, are also used to evaluate the efficacy of the workflow with a real dataset. The workflow identified 265 history matching solutions out of 1240 total simulation runs (21% acceptance ratio), where the results were used to characterize posterior parameter distribution and evaluate the prediction uncertainty. Furthermore, this workflow is significant because it enables strain tensor, or other geomechanical measurements to be interpreted to guide decision-making during energy and environmental processes in the subsurface.

42 ENGINEERING↗

DOE OSTI2024

A Novel Workflow for Coupled Simulation of Hydraulic Stimulation with Simultaneous Injection of Proppant

Enhanced Geothermal Systems (EGS) are located in impermeable host rock. Therefore, to achieve promising circulation, hydraulic stimulation is a key factor. Stimulation involves injecting a combination of fluids and proppant particles into existing fractures under high pressure. The effectiveness of stimulation depends on a variety of different factors, including the pressurization scheme and the optimal choice of carrying fluid and proppant particles to be used. For example, high-density particles can settle out of the injected slurry quickly, whereas lighter particles as well as a dense carrying fluid allow for a long transport into the reservoir. A numerical tool is therefore essential to take all these factors into account and to make robust predictions about the performance of the hydraulic stimulation. We present a novel computational application for coupled numerical simulation of hydraulic stimulation with concurrent injection of proppant in a fractured geothermal reservoir. This application is based on the open-source MOOSE framework that has been enriched with a set of equations to account for the circulation of proppant-laden slurries. It accounts for both the fluid and proppant components of the mixture, as well as particle-driven processes that include gravity settling, particle-particle interactions, and strong density and viscosity contrasts. A particular challenge is the mutual coupling with fracture mechanical processes in a propagating fracture using the FE code Defmod. Aperture and fluid pressure on one side and slurry density and viscosity changes on the other side interact and affect each other. Their basic understanding is therefore essential for the success of the stimulation. We compare our application to common numerical benchmark problems for both slurry flow and hydraulic stimulation.

15 GEOTHERMAL ENERGY↗

DOE OSTI2026

Stochastic fracture generation and thermo-hydro-mechanical modeling in an equivalent continuum framework for enhanced geothermal systems

Enhanced geothermal systems (EGS) involve fracturing low permeability material to establish well connectivity and then injecting and circulating fluid into the fractured subsurface for geothermal power production. Changes in fracture aperture from contraction of the cooling matrix rock may alter network connectivity and risk thermal short-circuiting. Thermo-hydro-mechanical (THM) models are a useful tool to study these processes. However, as fracture networks are complex, and data may be limited, fracture networks in THM models are often stochastically generated. Given reliance on stochastic fracture networks and THM modeling to represent the subsurface and assess productivity of EGS, increased understanding of the influence of such statistically derived fracture networks on flow and heat transport in THM models is needed. Here, a new fracture process model is developed in the reactive transport code PFLOTRAN to stochastically generate fracture families and simulate changes in fracture aperture over time due to temperature changes of the rock matrix. Sixty-four different fracture networks ranging from well to poorly-connected, are modeled in PFLOTRAN with and without mechanical processes (THM vs TH). Results indicate that for well-connected fracture networks, thermal short-circuiting is less of a concern due to the abundance of available alternative flowpaths. For poorly-connected fracture networks, inclusion of mechanical processes showed steep thermal drawdown coincident with increase in fracture aperture along developing colder flowpaths, demonstrating the risk of thermal short-circuiting. Simulations with additional, larger fractures engineered to establish connectivity in a poorly-fractured subsurface, indicate that while stochastic variation of fracture orientation of the background network had limited influence, such variation in the engineered fractures significantly affected flow and heat transport.

Discrete fracture networks (DFN)↗

DOE OSTI2021

Thermal loading analysis of the ring injection dump for the Spallation Neutron Source facility

The ring injection dump (RID) is the largest beam dump in the Spallation Neutron Source facility at Oak Ridge National Laboratory and accepts a fraction of the beam from the ring that is not captured in the ring during injection. Thermo-fluid modeling of the RID components was performed during the initial design by using the software ANSYS CFX to characterize the power rating for the RID. From the calculations, the power rating was lowered from 200 kW to 150 kW, at 1.3 GeV beam energy due to concerns regarding the heating of the concrete structure. The Proton Power Upgrade project scope included reevaluating the power limit to the RID in anticipation that increasing it would provide more operational flexibility. The focus was on the shielding because it was the limiting factor, and data collected from thermocouples installed on the structure were used to benchmark the analysis. A three-stage process was adopted in the latest study. First, the already existing steady-state calculations were validated and extended with the latest Monte Carlo N-Particle source term calculations. Second, transient calculations were conducted to capture the dynamic state of the system in response to the energy dumped to the RID during the following operational period of the beam from 2005–2019 (14 years), and sensitivity analyses of crucial parameters were performed to benchmark the model with reasonable accuracy. Finally, an idealized transient cycle analysis was conducted with realistic duty factors to predict the temperature distribution at higher beam powers. Overall, the results confirmed the accuracy of the original steady state model but allowed the development of a new validated transient model with higher accuracy for future analyses.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

DOE OSTI2022

Scaling Behavior of Thermally Driven Fractures in Deep Low‐Permeability Formations: A Plane Strain Model With 1‐D Heat Conduction

Abstract Injection of cold fluids through/into deep formations may cause significant cooling, thermal stress, and possible thermal fracturing. In this study, the thermal fracturing of low‐permeability formations under one‐dimensional heat conduction was investigated using a plane strain model. Dimensionless governing equations, with dimensionless fracture length , aperture , spacing , time , and effective confining stress , were derived. Solution of single thermal fracture was derived analytically, while solution of multiple fractures with constant (or dynamic) spacing were obtained using the displacement discontinuity method (and stability analysis). For single fracture, increases nonlinearly with and then transitions to scaling law , indicating that late‐time fracture length increases linearly with the square root of cooling time. For constantly spaced fractures, deviates from the single‐fracture solution at a later for a larger , showing slower propagation under inter‐fracture stress interaction . For dynamically spaced fractures, fracture arrest induced by stress interaction was determined by the stability analysis; the fully transient solution provides evolution of dimensionless fracture length, spacing, aperture, and pattern; a similar scaling law, with , obtained shows the effect of both stress interaction and fracture arrest. The solution and scaling law provide fast predictions for all reservoir and cooling conditions using (single) model parameter . Application to a geothermal site with demonstrates that thermal fractures reach 0.67, 6.25, and 78.00 m in length, 0.49, 2.30, and 13.00 m in spacing, and 0.43, 2.09, and 12.19 mm in aperture at 1, 100, and 10,000 days.

Chen, Bin↗

DOE OSTI2022

Shallow Aseismic Slip in the Delaware Basin Determined by Sentinel-1 InSAR

The Delaware Basin, Texas is currently a hot-spot of induced seismicity and ground deformation due to fluid extraction and injection associated with horizontal drilling techniques; however, the driving mechanism behind the seismicity and deformation remains under debate. Here, using vertical and east-west horizontal surface deformation measurements derived from Sentinel-1 interferometric synthetic aperture radar (InSAR), we show that the subsurface responds differently to oil and gas activity in the northern and southeastern portions of the basin. In the north, where there is little seismicity, deformation patterns display long-wavelengths and equidimensional patterns. In contrast, the southeast region hosts most of the seismicity and displays spatial deformation patterns with narrow linear features that strike parallel to the maximum principal horizontal stress and to trends in seismicity, suggesting movement along normal faults. We model a linear deformation feature using edge dislocations and show that the InSAR observations can be reproduced by slip on normal faults contained within the Delaware Mountain Group (DMG), the formation that hosts local wastewater injection and the majority of earthquakes. Our model consists of three parallel, high-angle normal faults, with two dipping toward one another in a graben structure. Slip magnitudes reach up to 25 cm and are spatially correlated with injection wells. Measured seismicity can only explain ~2% of the fault motion predicted by our fault model, suggesting that slip leading to the deformation is predominantly aseismic. We conclude that seismic and aseismic fault motion in the southeastern Delaware Basin is likely driven by wastewater injection near critically-stressed normal faults within the DMG.

58 GEOSCIENCES↗

DOE OSTI2021

Key Learnings from Hydraulic Fracturing Test Site - 2 (HFTS-2), Delaware Basin

Hydraulic fractures (HF) play a dominant role in oil and gas production from unconventional reservoirs, fracture characteristics, geometry, and extent are generally not well known and are typically inferred from various monitoring methods or subsurface models. The Hydraulic Fracturing Test Site-2 (HFTS-2) is a field-based research experiment performed in the Wolfcamp formation of the Permian (Delaware) Basin. HFTS-2, a unique joint industry program (JIP) that brings together government, academia, and industry, focused on field data acquisition and interpretation. The project has built a comprehensive HF research dataset that includes microseismic (MS), 1,500 ft of core, diagnostic formation injection tests (DFIT), fluid characterization, advanced logs, geochemistry, downhole gauges, varying completions, and permanent fiber optic (FO) cable in vertical and horizontal wells. These data were integrated and used to characterize HF geometry/dimensions and resulting depletion profiles, gain insights for completion design, and achieve a better understanding of parent well depletion effects on child well performance. The analysis and integration of the comprehensive HFTS-2 diagnostics dataset provide insights for optimal well placement and completions and improves understanding about the effects of parent well depletion. Furthermore, this dataset can be used to blind test and develop subsurface workflows for unconventional resources. Key learnings from HFTS-2 and the new technologies tested can easily be transferred to other unconventional plays.

58 GEOSCIENCES↗

DOE OSTI2023

NRAP-Open-IAM Multisegmented Wellbore Reduced-Order Model: Improvement and Quality Assurance

The multisegmented wellbore model (MSW) semi-analytically estimates the amount of CO 2 and brine leakage from a leaking legacy well by segmenting it into intervals to simulate site-specific stratigraphic and hydrogeologic properties. The model is a component of the National Risk Assessment Partnership Open-Source Integrated Assessment Model (NRAP-Open-IAM), which was developed to perform risk assessment for geologic CO 2 storage. The new wellbore leakage model, which uses deep learning networks for a caprock segment, was developed to enhance the analytical MSW. The model was trained and validated using a synthetic data set of Subsurface Transport Over Multiple Phases (STOMP) multiphase flow simulations from various geological, well attribute, and operational conditions to ensure its quality. The results demonstrate that the model is more accurate than the existing model in predicting the transport of two-phase fluids (brine and injected CO 2 ) through the well. This report provides a detailed explanation of the model development and quality assurance.

58 GEOSCIENCES↗

DOE OSTI2024

Pore architecture controls on mineral reactivity

Mineral dissolution rates measured in natural environments are much slower than those measured in laboratory settings. This project tested the hypothesis that the way fluid flows through rocks in natural systems creates areas where mineral dissolution is fast and areas where mineral dissolution is slow. This hypothesis was tested with a combination of laboratory experiments and numerical simulation. We demonstrated a separation of fluid flow pathways and rates of mineral dissolution in laboratory experiments for the first time using an experimental approach where we created synthetic rocks that have different ratios of connected and dead-end pathways for fluid flow. In laboratory experiments with higher proportions of dead-end pathways, the mineral dissolution rates were slower. We found that where fluid flows through connected pathways the continuous refreshing of fluid at the mineral surface creates conditions where dissolution is fast. In contrast, where fluid either flows slowly through poorly connected pathways or is stagnant in dead-end pathways, mineral dissolution is slow. The results from this project suggest that the overall slowing of rates of mineral dissolution is important when the proportion of dead-end pathways is greater than ~40%. This project informs our understanding of the way that fluids react with rocks in carbon dioxide sequestration and enhanced geothermal projects where fluids are purposefully injected into rocks for energy applications.

58 GEOSCIENCES↗

DOE OSTI2021

Findings from the NETL Well Integrity Workshop: Insights for multiple subsurface energy operations

Wells are a core aspect of subsurface energy infrastructure. Ensuring the resilience of that infrastructure relies on the integrity of wells as engineered pathways that safely extract, inject, and store fluids in the subsurface. To identify areas where focused research can yield new insights to improve integrity and extend the operational life of wells, the United States (US) Department of Energys (DOE) National Energy Technology Laboratory (NETL) organized and hosted a workshop on well integrity research needs in May and June of 2021. The three-day workshop consisted of six technical sessions and 18 presentations from international well integrity experts on topics of geologic storage of carbon dioxide, natural gas, and hydrogen, as well as hydrocarbon and geothermal energy production. There were over 350 registrants from more than 10 countries, 4 major oil & gas companies, 14 US states, 8 US DOE National Laboratories, US Departments of Interior, Transportation, and Energy, and regulators from eight of ten US Environmental Protection Agency (EPA) regions and US EPA headquarters. Through presentations, moderated discussions, and an attendee survey several important areas for future research were identified. In this presentation, we highlight the findings of the NETL well integrity workshop and discuss future well integrity research needs.

carbon capture utilization and storage (CCUS)↗

NASA NTRS1962

The Analysis of a Vortex Type Magnetohydrodynamic Induction Generator

Consideration it is given to the performance to the characteristics of an AC magnetohydrodynamic power generator, A rotating magnetic field is imposed on the vortex flow of an electrically conducting fluid, which is injected tangentially into an annulus formed by two nonconducting concentric cylinders and two nonconducting end plates. A perturbation technique is used to determine the two dimensional velocity and three dimensional electromagnetic field and current distributions. Finally, the generated power, the ohmic losses, the effective power and the electrical efficiency of the converter system are calculated.

MAGNETOHYDRODYNAMIC GENERATOR↗

NASA NTRS1987

Preparation of polystyrene microspheres for laser velocimetry in wind tunnels

Laser Velocimetry (L/V) had made great strides in replacing intrusive devices for wind tunnel flow measurements. The weakness of the L/V has not been the L/V itself, but proper size seeding particles having known drag characteristics. For many Langley Wind Tunnel applications commercial polystyrene latex microspheres suspended in ethanol, injected through a fluid nozzle provides excellent seeding but was not used due to the high cost. This paper provides the instructions, procedures, and formulations for producing polystyrene latex monodisperse microspheres of 0.6, 1.0, 1.7, 2.0, and 2.7 micron diameters. These are presently being used at Langley Research Center as L/V seeding particles.

Nichols, Cecil E., Jr.↗

NASA NTRS2002

Reassessing the Evidence for the Earliest Traces of Life

The isotopic composition of graphite is commonly used as a biomarker in the oldest (>3.5 Gyr ago) highly metamorphosed terrestrial rocks. Earlier studies on isotopic characteristics of graphite occurring in rocks of the approximately 3.8-Gyr-old Isua supracrustal belt (ISB) in southern West Greenland have suggested the presence of a vast microbial ecosystem in the early Archean. This interpretation, however, has to be approached with extreme care. Here we show that graphite occurs abundantly in secondary carbonate veins in the ISB that are formed at depth in the crust by injection of hot fluids reacting with older crustal rocks (metasomatism). During these reactions, graphite forms from the disproportionation of Fe(II)-bearing carbonates at high temperature. These metasomatic rocks, which clearly lack biological relevance, were earlier thought to be of sedimentary origin and their graphite association provided the basis for inferences about early life. The new observations thus call for a reassessment of previously presented evidence for ancient traces of life in the highly metamorphosed Early Archaean rock record.

vanZullen, Mark A.↗

NASA NTRS2017

Imaging a Free-Geometry Magma Reservoir Using FEM-Based Linear Inverse Modeling of InSAR Data: Application to Rabaul Caldera, Papua New Guinea

We test an innovative inversion scheme using Green’s functions from an array of pressure sources embedded in finite-element method (FEM) models to image, without assuming an a-priori geometry, the composite and complex shape of a volcano deformation source. We invert interferometric synthetic aperture radar (InSAR) data to estimate the pressurization and shape of the magma reservoir of Rabaul caldera, Papua New Guinea. The results image the extended shallow magmatic system responsible for a broad and long-term subsidence of the caldera between 2007 February and 2010 December. Elastic FEM solutions are integrated into the regularized linear inversion of InSAR data of volcano surface displacements in order to obtain a 3-D image of the source of deformation. The Green’s function matrix is constructed from a library of forward line-of-sight displacement solutions for a grid of cubic elementary deformation sources. Each source is sequentially generated by removing the corresponding cubic elements from a common meshed domain and simulating the injection of a fluid mass flux into the cavity, which results in a pressurization and volumetric change of the fluidfilled cavity. The use of a single mesh for the generation of all FEM models avoids the computationally expensive process of non-linear inversion and remeshing a variable geometry domain. Without assuming an a-priori source geometry other than the configuration of the 3-D grid that generates the library of Green’s functions, the geodetic data dictate the geometry of the magma reservoir as a 3-D distribution of pressure (or flux of magma) within the source array. The inversion of InSAR data of Rabaul caldera shows a distribution of interconnected sources forming an amorphous, shallow magmatic system elongated under two opposite sides of the caldera. The marginal areas at the sides of the imaged magmatic system are the possible feeding reservoirs of the ongoing Tavurvur volcano eruption of andesitic products on the east side and of the past Vulcan volcano eruptions of more evolved materials on the west side. The interconnection and spatial distributions of sources correspond to the petrography of the volcanic products described in the literature and to the dynamics of the single and twin eruptions that characterize the caldera. The ability to image the complex geometry of deformation sources in both space and time can improve our ability to monitor active volcanoes, widen our understanding of the dynamics of active volcanic systems and improve the predictions of eruptions.

Inverse theory↗