Search NASA⌕ Search

SEARCH · Search NASA

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.

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

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Stress Controls Rupture Extent and Maximum Magnitude of Induced Earthquakes

Abstract Seismic hazard forecasts of induced seismicity often require estimates of the maximum possible magnitude ( M max ). Empirical models suggest that maximum magnitudes, or expected number of earthquakes, are related to the volume of injected fluid. We perform a suite of 3D physics‐based earthquake simulations with rate‐ and state‐dependent friction, systematically varying the area of the pressurized region and the amplitude of the initial homogeneous or heterogeneous shear stress. Using the resulting catalog we explore the conditions that result in pressure‐controlled versus runaway ruptures that extend outside the pressurized zone. We find that proposed empirical scaling laws correctly predict M max when shear stresses are further from failure (≤90% of maximum shear stress) and for high amplitude stress fields. Runaway ruptures are observed for higher initial shear stresses and smoother stress fields. In these cases, runaway ruptures occur early after the onset of injection and rarely preceded by foreshock activity.

58 GEOSCIENCES↗
DOE OSTI2021

Seismic Diffusivity and the Influence of Heterogeneity on Injection‐Induced Seismicity

Abstract The spatiotemporal patterns of injection‐induced seismicity (IIS) are commonly interpreted with the concept of a triggering front, which propagates in a diffusion‐like manner with an associated diffusivity parameter. Here, we refer to this diffusivity as the “seismic diffusivity.” Several previous studies implicitly assume that seismic diffusivity is equivalent to the effective hydraulic diffusivity of the subsurface, which describes the behavior of the mean pressure field in heterogeneous porous media. Seismicity‐based approaches for hydraulic characterization or simulations of IIS using domains of homogeneous equivalent porous media are implicitly based on this assumed equivalence. However, seismicity is expected to propagate with the threshold triggering pressure, and thus not be controlled by the evolution of the mean pressure field. We present numerical simulations of fluid injection to compare the seismic and effective hydraulic diffusivities in heterogeneous formations (including fractured rock). The numerical model combines uncoupled, linear pressure diffusion with the Mohr‐Coulomb failure criterion to simulate IIS. We demonstrate that connected pathways of relatively high hydraulic diffusivity in heterogeneous media (particularly in fractured rock domains) allow the threshold triggering pressure to propagate more rapidly than predicted by the effective hydraulic diffusivity. As a result, the seismic diffusivity is greater than the effective hydraulic diffusivity in heterogeneous porous media, possibly by an order of magnitude or more. Additionally, we present a case study of IIS near Soultz‐sous‐Forêts where seismic diffusivity is found to be at least one order of magnitude larger than the effective hydraulic diffusivity.

Haagenson, Ryan↗
DOE OSTI2011

AFR-100 Sodium-CO 2 Reaction Analysis

Pressure-boundary failure within the heat-source heat exchanger (IHTS-fluid to secondary-fluid heat exchanger) in a sodium-cooled reactor will lead to a violent decompression on the secondary side and large pressure transients in the IHTS. If the secondary fluid is water or CO 2 , the pressure transients will be intensified by an energetic reaction between the sodium & water or sodium & CO 2 . The coupling of secondary-fluid injection, reaction dynamics, product gas bubble growth, and sodium system response results in a succession of short duration, high intensity pressure pulses travelling in the heat-source heat exchanger and IHTS. Each pulse has higher amplitude than the preceding one if no mitigating features such as rupture disks and gas spaces are provided in the system design. Pressure pulse absorption features and relief systems must be included in the IHTS design to prevent structural damage to the heat-source heat exchanger, and IHTS piping. In the case of the AFR-100 design, the boundary between the IHTS sodium and the secondary-side CO 2 within the sodium-CO 2 heat exchanger could potentially fail. The resulting fluid-hammer and sodium-CO 2 reaction must be mitigated. A variety of heat exchangers types have been proposed for the sodium-CO 2 side. Thus, the following analysis of the IHTS response to a sodium-CO 2 reaction includes a wide range of potential break sizes to cover this spectrum of heat exchanger types. The IHTS pressure response due to leak sizes ranging from 2.0 mm – 16.1 mm were evaluated. Double-ended guillotine tube breaks at the top, center, and bottom of the tube bundle were also considered. As expected, the peak reaction bubble pressure increased and rupture disk burst time decreased as the break size increased. As currently designed, all components, including the sodium-CO 2 heat exchanger tubes and shell and the IHX tubes, can withstand the maximum pressures calculated.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Characterization of induced seismicity patterns derived from internal structure in event clusters

Microseismicity induced by CO 2 injection at Decatur, Illinois, occurs in distinct clusters and shows no obvious correlation with the proceeding pressure front. Here, we analyze some of these clusters in more depth by using a waveform cross-correlation approach. With this approach we can associate about 1400 events from two clusters, with moment magnitudes between 1.1 and -1.7, with specific formations of much smaller vertical dimensions (tens of meters) than the depth resolution of traveltime-based event locations. The differentiation of reservoir and basement events, and the definition of subclusters by waveform correlation, rather than by location, helps to better analyze the spatiotemporal evolution of the events within a cluster. In the Decatur case, this is characterized by event migration from the reservoir into the adjacent basement. The spatial variation of Brune stress drop and Gutenberg b value exhibits signs of a fluid-driven triggering mechanism at the cluster level, revealing a punctual hydraulic connection between reservoir and basement, most likely associated with basement faults cutting into the reservoir. The observed clustering of microseismicity can thus be explained by the lateral heterogeneity of permeability and crustal strength and is overall consistent with a pressure-induced triggering mechanism. Hence, proper long-term risk mitigation for large-scale fluid injection close to the basement requires prior mapping of small subseismic basement-connected faults.

58 GEOSCIENCES↗
DOE OSTI2017

Evaluation of Stimulation by Shear Slip in Fractured Rock Using a 3D Coupled Thermo-Poro-Mechanical FEM

This paper presents a 3D finite element model for simulating shear slip on fractures in rock in response to coupled thermo-poro-mechanical processes. The modeling is part of a broader effort to investigate the role of pore pressure and cooling by cold water injection on shear slip and permeability increase in granitic rock using laboratory shearing experiments under triaxial conditions. In particular, 3D thermo-poro-mechanical finite element modeling and analysis of injection experiment in fractured rock has been carried out to analyze the role of pore pressure, temperature, fracture deformation and their interactions. 4-noded tetrahedral elements are employed for intact rock deformation and transport processes within the matrix. To represent the mechanical response of a fracture surface, zero thickness contact interface element is developed based on recently developed element partitioning algorithm and procedures for applying hydraulic pressure on the fracture surfaces. A Mohr-Coulomb type criterion is adapted to capture the slip on the fracture and to quantify its reversible and irreversible deformation. The numerical model has been calibrated using results from well controlled, advanced laboratory experiments. Excellent agreement between modeling and experimental observations is achieved. Simulation results illustrate that pore pressure increase and rock matrix/fracture surface cooling cause the fracture system to deform and slip. Fracture slip is promoted due to its normal stress reduction associated with cooling effect of cold fluid injection. The numerical model provides a physics-based understanding of the role of coupled processes on shear stimulation phenomenon and the resulting permeability enhancement.

58 GEOSCIENCES↗
DOE OSTI2023

Fracture-matrix fluid exchange in oil-bearing unconventional mudstones

The poromechanical properties of unconventional reservoir materials are in large part dictated by their mineralogy. Since these properties govern the response to stress experienced during hydraulic fracturing, fluid production, and fluid injection, they play a central role in the formation of microcracks or bedding delaminations which ultimately dominate mass transport. In this work we study access to the porosity of end member unconventional reservoir materials, where the end members are predominantly dictated by carbonate content. Access to the porosity is quantified using state of the art 3D x-ray computed tomography coupled with physics informed data analytics. Xenon gas, which attenuates x-rays, provides a spatiotemporal map of access to the porosity. The accessible porosity is quantified over a range of net confining stress relevant to the manmade disturbances listed above. These experiments demonstrate that heavily carbonated mudstones are nearly impermeable at the core (~ cm) scale, while carbonate free analogues afford better access to the microstructure. Consistent with previous qualitative 2D radiographs, access to the interior of the clastic mudstones is first observed along planar microcracks, followed by slow penetration into the surrounding matrix. Physics informed data analytics of the 3D tomography measurements presented here show that these microcracks do not permit uniform access to the adjacent rock matrix. In addition, variation of the effective pressure elucidates the mechanisms that govern fracture/matrix fluid exchange. Under conditions consistent with hydrocarbon production fluid accumulates in the immediate vicinity of the nearest microcrack. While there is clear evidence that, as intended, part of this accumulation is from the more distant matrix, fluid is also squeezed out of the microcrack. The fluid build-up at the microcrack indicates that migration out of the rock is hindered by the coupled poroelastic response of the microcrack and adjacent rock matrix. We show that these mechanisms ultimately account for the meager oil recovery factors realized in practice. These insights have implications for making reservoir scale predictions based on core scale observations, and provide a basis for devising new asset development techniques to access more porosity, and enhance fluid extraction. Finally, these findings shed light on key features and mechanisms that govern shale storage capacity, with relevance to other important industrial processes, such as geologic CO 2 storage.

58 GEOSCIENCES↗
DOE OSTI2017

Frictional stability-permeability relationships for fractures in shales

There is wide concern that fluid injection in the subsurface, such as for the stimulation of shale reservoirs or for geological CO 2 sequestration (GCS), has the potential to induce seismicity that may change reservoir permeability due to fault slip. However, the impact of induced seismicity on fracture permeability evolution remains unclear due to the spectrum of modes of fault reactivation (e.g., stable versus unstable). As seismicity is controlled by the frictional response of fractures, we explore friction-stability-permeability relationships through the concurrent measurement of frictional and hydraulic properties of artificial fractures in Green River shale (GRS) and Opalinus shale (OPS). We observe that carbonate-rich GRS shows higher frictional strength but weak neutral frictional stability. The GRS fracture permeability declines during shearing while an increased sliding velocity reduces the rate of permeability decline. By comparison, the phyllosilicate-rich OPS has lower friction and strong stability while the fracture permeability is reduced due to the swelling behavior that dominates over the shearing induced permeability reduction. Furthermore, we conclude that the friction-stability-permeability relationship of a fracture is largely controlled by mineral composition and that shale mineral compositions with strong frictional stability may be particularly subject to permanent permeability reduction during fluid infiltration.

58 GEOSCIENCES↗
DOE OSTI2019

Field experiments of surface water to groundwater recharge to characterize the mobility of uranium and vanadium at a former mill tailing site

Characterizing the mobility of uranium and vanadium in groundwater with a hydraulic connection to surface water is important to inform the best management practices of former mill tailing sites. In this study, the recharge of river water to the unsaturated and saturated zones of a uranium-contaminated alluvial aquifer was simulated in a series of forced-gradient single- and multi-well injection-extraction tests. The injection fluid (river water) was traced with natural and artificial tracers that included halides, fluorobenzoates, lithium, and naphthalene sulfonate to characterize the potential mass transport mechanisms of uranium and vanadium. The extraction fluid (river water/groundwater mixture) was analyzed for the tracers, uranium, and vanadium. The results from the tracers indicated that matrix diffusion was likely negligible over the spatiotemporal scales of the tests as evident by nearly identical breakthrough curves of the halides and fluorobenzoates. In contrast, the breakthrough curves of lithium and naphthalene sulfonate indicated that sorption by cation exchange and sorption to organic matter, respectively, were potential mass transport mechanisms of uranium and vanadium. Uranium was mobilized in the saturated zone containing gypsum (gypsum-rich zone), the vadose zone (vadose-rich zone), and the saturated zone containing organic carbon (organic-rich zone) whereas vanadium was mobilized only in the saturated gypsum-rich zone. The mechanisms responsible for the mobilization of uranium and vanadium were likely dissolution of uranium- and vanadium-bearing minerals and/or desorption from the gypsum-rich zone, flushing of uranium from the vadose-rich zone, and desorption of uranium from the organic-rich zone due to the natural contrast in the geochemistry between the river water and groundwater. The experimental design of this study was unique in that it employed the use of multiple natural and artificial tracers coupled with a direct injection of native river water to groundwater. Here, these results demonstrated that natural recharge and flooding events at former mill tailing sites can mobilize uranium, and possibly vanadium, and contribute to persistent levels of groundwater contamination.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Thermal Effects on Far-Field Distributed Acoustic Strain-Rate Sensors

Summary Fiber-optic cables cemented outside of the casing of an unconventional well measure crosswell strain changes during fracturing of neighboring wells with low-frequency distributed acoustic sensing (LF-DAS). As a hydraulic fracture intersects an observation well instrumented with fiber-optic cables, the fracture fluid injected at ambient temperatures can cool a section of the sensing fiber. Often, LF-DAS and distributed temperature sensing (DTS) cables are run in tandem, enabling the detection of such cooling events. The increasing use of LF-DAS for characterizing unconventional hydraulic fracture completions demands an investigation of the effects of temperature on the measured strain response by LF-DAS. Researchers have demonstrated that LF-DAS can be used to extract the temporal derivative of temperature for use as a differential-temperature-gradient sensor. However, differential-temperature-gradient sensing is predicated on the ability to filter strain components out of the optical signal. In this work, beginning with an equation for optical phase shift of LF-DAS signals, a model relating strain, temperature, and optical phase shift is explicitly developed. The formula provides insights into the relative strength of strain and temperature effects on the phase shift. The uncertainty in the strain-rate measurements due to thermal effects is estimated. The relationship can also be used to quantify uncertainties in differential-temperature-gradient sensors due to strain perturbations. Additionally, a workflow is presented to simulate the LF-DAS response accounting for both strain and temperature effects. Hydraulic fracture geometries are generated with a 3D fracture simulator for a multistage unconventional completion. The fracture width distributions are imported by a displacement discontinuity method (DDM) program to compute the strain rates along an observation well. An analytic model is used to approximate the temperature in the fracture. Using the derived formulae for optical phase shift, the model outputs are then used to compute the LF-DAS response at a fiber-optic cable, enabling the generation of waterfall plots including both strain and thermal effects. The model results suggest that before, during, and immediately following a fracture intersecting a well instrumented with fiber, the strain on the fiber drives the LF-DAS signal. However, at later times, as completion fluid cools the observation well, the temperature component of the LF-DAS signal can be equal to or exceed the strain component. The modeled results are compared to a published field case in an attempt to enhance the interpretation of LF-DAS waterfall plots. Finally, we propose a sensing configuration to identify the events when “wet fractures” (fractures with fluids) intersect the observation well.

Engineering↗
DOE OSTI2016

Increased stream discharge after the 3 September 2016 M w 5.8 Pawnee, Oklahoma earthquake

Earthquakes influence hydrogeological processes and properties in Earth's crust, some of which affect surface waters. Here we document increased discharge in a stream after the 3 September 2016 M w 5.8 earthquake near Pawnee, Oklahoma, an event likely induced by underground wastewater disposal. Discharge increased by an order of magnitude and remained elevated until the change was obscured by rain 1 week later. Given the earthquake magnitude and distance from the stream, by comparison with previous examples of responses to earthquakes, increased discharge after this earthquake is expected. While the mechanism increasing discharge cannot be confirmed, the observations require changes in physical properties of the subsurface. Fluid injection may thus influence hydrogeological properties of shallow groundwater systems and aquifers indirectly by inducing seismicity, if the induced seismic events are large enough.

58 GEOSCIENCES↗
DOE OSTI2016

Validating induced seismicity forecast models—Induced Seismicity Test Bench

Abstract Induced earthquakes often accompany fluid injection, and the seismic hazard they pose threatens various underground engineering projects. Models to monitor and control induced seismic hazard with traffic light systems should be probabilistic, forward‐looking, and updated as new data arrive. In this study, we propose an Induced Seismicity Test Bench to test and rank such models; this test bench can be used for model development, model selection, and ensemble model building. We apply the test bench to data from the Basel 2006 and Soultz‐sous‐Forêts 2004 geothermal stimulation projects, and we assess forecasts from two models: Shapiro and Smoothed Seismicity (SaSS) and Hydraulics and Seismics (HySei). These models incorporate a different mix of physics‐based elements and stochastic representation of the induced sequences. Our results show that neither model is fully superior to the other. Generally, HySei forecasts the seismicity rate better after shut‐in but is only mediocre at forecasting the spatial distribution. On the other hand, SaSS forecasts the spatial distribution better and gives better seismicity rate estimates before shut‐in. The shut‐in phase is a difficult moment for both models in both reservoirs: the models tend to underpredict the seismicity rate around, and shortly after, shut‐in.

Király‐Proag, Eszter↗
DOE OSTI2018

Effectiveness of supercritical-CO 2 and N 2 huff-and-puff methods of enhanced oil recovery in shale fracture networks using microfluidic experiments

Current oil recovery methods in hydraulically fractured shale reservoirs have a low recovery efficiency of about 10%. The objective of this work is to investigate the effectiveness of nitrogen and supercritical carbon dioxide in the huff-and-puff method for enhanced oil recovery as a means to re-energize reservoirs and improve recovery rates. Here, we conduct direct visualization experiments with a microfluidic system to reveal the mechanisms and to quantify the recovery rates of oil from fracture networks. We compared the effectiveness of water, nitrogen and supercritical carbon dioxide at reservoir conditions in a process mimicking the huff-and-puff method in both dead-end and connected fracture systems. The microfluidic chips were made of glass and placed in a confining pressure system pressurized to 10 MPa, 50 °C. The system was allowed to equilibrate, and then depressurized to simulate huff-and-puff oil recovery. Fluorescence microscope images were continuously taken to visualize and calculate residual oil saturation as a function of pressure drawdown. As the system was depressurized from 10 MPa, gas exsolution from the oil liquid phase, including bubble nucleation, growth, and coalescence, appeared to be the main driver for mobilizing oil from the fracture networks. Injection of supercritical CO 2 resulted in the highest recovery rate with an average end-point recovery of about 90% in the connected fracture network and 60% in the dead-end fracture network. N 2 has lower solubility in oil and hence showed a lower recovery rate of 40% in the connected fracture network and 25% in the dead-end fracture network. Injection of water had no effect on oil mobilization since water is insoluble, immiscible and incompressible. The main mechanism of enhanced recovery was gas exsolution from the liquid phase as pressure was decreased below the bubble point pressure. Because the gas was distributed throughout the oil phase, bubble nucleation, growth, coalescence, and elongation occurred throughout the fracture network. Expansion of the gas forced oil out of the network through piston displacement in continuous oil areas and film flow in the dispersed oil areas. Bubbles began as spheres and grew until they touched the fracture walls where they elongated along the fracture length. The bubble growth rate depended on local mass transfer from liquid to gas phase and gas volume expansion due to pressure drop. The efficiency of the huff-and-puff process is dependent on the solubility and miscibility of the injection fluid with oil. High gas solubility allows for more bubble nucleation, growth and expansion during the depressurization cycle.

58 GEOSCIENCES↗
DOE OSTI2024

Investigation of a high-temperature combination heat pump for lower-cost electrification in multifamily buildings

The development of space and water heating combination heat pumps capable of generating water temperatures high enough for convective heat emitters will enable more cost-effective and equitable decarbonization solutions for electrifying multifamily buildings. Here, in this paper, multifamily building models and a charge-sensitive mechanistic cycle model of a combination heat pump are developed, and the system performance is predicted based on the models. Unlike other state-of-the-art residential heat pumping equipment, the modeled combination heat pump using an economized, fluid-injected variable-speed compressor can achieve higher temperature lifts of 40° - 85°C, with lower installation costs and complexity. The model predicted heating coefficient of performance (COP h ) is 2.1 at an ambient temperature of -15°C with a high-temperature lift of nearly 85°C, and a seasonal coefficient of performance in heating mode (SCOP h ) ranges from 2 - 4 for different locations. The system shows 30% - 90% lower CO 2 eq emissions over a condensing gas boiler and 9% - 13% lower projected installation costs than two separate space and water heat pumping appliances.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Experimental evaluation of the role of redox during glauconite-CO 2 -brine interactions

Greensands formations are globally abundant sedimentary rocks rich in Fe clays (typically glauconite) that commonly contain natural hydrocarbon accumulations and may be important reservoirs for geologic storage of anthropogenic CO 2 . Diagenesis in greensands is commonly accompanied by the conversion of primary glauconite to siderite (FeCO 3 ), a process that could be exploited for the permanent trapping of CO 2 . Importantly, siderite formation after glauconite requires that the mostly oxidized Fe in the primary Fe clay minerals is reduced during diagenetic interactions. Here, we explore the effect of solution redox state on the stability of glauconite in sandstones with implications for the diagenetic and/or engineered formation of siderite. We performed two flow-through experiments on intact, glauconite-rich sandstone cores at 150 °C and 150 bar. Both experiments employed a 1 mol NaCl/kg, 0.1 mol NaHCO 3 /kg solution charged with ~0.58 mol CO 2 /kg solution, but the redox state of the injected fluid was manipulated between experiments in order to compare glauconite reactivity and siderite saturation state at oxidizing and reducing end-member conditions. After reaction with the oxidizing (O 2 (aq) ≈ 6 μmol/kg) fluid, chemical and Mӧssbauer spectroscopic analyses indicate the production of Fe(III)-oxy/hydroxide minerals from glauconite, whereas, in the reducing (H 2 (aq) ≈ 5–40 mmol/kg) experiment, thermodynamic calculations and coupled chemical, mineralogical, and Mӧssbauer analyses suggest glauconite dissolution and precipitation of an Fe(II) mineral, likely siderite, and minor magnetite formation. These experimental results, along with thermodynamic calculations, confirm that solution redox state is the master variable dictating siderite formation in greensands.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

TOUGH–RBSN simulator for hydraulic fracture propagation within fractured media: Model validations against laboratory experiments

This study presents coupled hydro-mechanical modeling of hydraulic fracturing processes in complex fractured media using a discrete fracture network (DFN) approach. The individual physical processes in the fracture propagation are represented by separate program modules: the TOUGH2 code for multiphase flow and mass transport based on the finite volume approach; and the rigid-body-spring network (RBSN) model for mechanical and fracture-damage behavior, which are coupled with each other. Fractures are modeled as discrete features, of which the hydrological properties are evaluated from the fracture deformation and aperture change. The verification of the TOUGH–RBSN code is performed against a 2D analytical model for single hydraulic fracture propagation. Subsequently, modeling capabilities for hydraulic fracturing are demonstrated through simulations of laboratory experiments conducted on rock-analogue (soda-lime glass) samples containing a designed network of pre-existing fractures. Sensitivity analyses are also conducted by changing the modeling parameters, such as viscosity of injected fluid, strength of pre-existing fractures, and confining stress conditions. The hydraulic fracturing characteristics attributed to the modeling parameters are investigated through comparisons of the simulation results.

58 GEOSCIENCES↗
DOE OSTI2017

CO 2 injection-induced fracturing in naturally fractured shale rocks

In this study, Niobrara shale cubes of 20 cm from Colorado were employed to investigate gas and supercritical CO 2 injection-induced fracturing in naturally fractured caprocks of deep aquifers/depleted reservoirs and fractured shale reservoirs. Under tri-axial stresses, gas or supercritical CO 2 was injected into the center of the cubes to induce fracturing. Real-time pressure and temperature, acoustic wave, pressure decay, fracture coloring, and gas fracturing were used to characterize the fracturing process and fracture morphology. Without pore pressure, CO 2 injection-induced fracturing occurred and completed instantly, accompanied by an evident temperature drop. Strongly bonded fractures barely affected transverse fracture propagation, whereas weakly bonded or open fractures arrested the injected fluid first and then allowed it to generate new fractures perpendicular to the minimum horizontal stress. Breakdown pressures for cubes with preexisting fractures using gas and supercritical CO 2 are much lower than both poroelastic predictions and slick-water fracturing pressure, and some are even lower than the minimum horizontal stress. This is attributed to unconformable preexisting fractures and the low viscosity of CO 2 . Moreover, decreasing tri-axial stress levels and increasing stress differences tend to lower the breakdown pressure. In conclusion, this study is instructive for understanding and tackling geomechanical issues related to CO 2 geological storage and fracturing of shale reservoirs.

58 GEOSCIENCES↗
DOE OSTI2022

Improved geophysical monitoring of carbon sequestration through parameter linkage to reservoir modeling

Predictive reservoir modeling, even if present in the form of only basic hydrogeological model assumptions, is expected to accompany the majority of carbon capture and sequestration monitoring activities. It thus represents a source of prior information about the migration of injected fluids that can benefit geophysical survey planning and ensuing monitoring. Constraining the imaging of geophysical monitoring data with reservoir modeling is preferable over standalone geophysical imaging because of additional complementary hydrogeological information. However, fully coupled hydrogeophysical data inversion for flow-modeling parameters that control saturation predictions is an involved process. Within the context of three-dimensional electromagnetic (EM) inversion of data from borehole-to-surface layouts, we employ a "poor people's" alternative. The approach constrains geophysical inversion parameters through saturation predictions. The coupling is realized through spatially variable lower and upper parameter bounds that scale with gas saturation magnitudes, the latter provided by reservoir modeling. In this work, enhancement of three-dimensional time-lapse plume EM imaging is demonstrated for simulated sequestration into a depleted gas reservoir.

58 GEOSCIENCES↗
DOE OSTI2025

A review of thermo-hydro-mechanical modeling of coupled processes in fractured rock: From continuum to discontinuum perspective

Coupled thermo-hydro-mechanical (THM) processes in fractured rock are playing a crucial role in geoscience and geoengineering applications. Diverse and conceptually distinct approaches have emerged over the past decades in both continuum and discontinuum perspectives leading to significant progress in their comprehending and modeling. This review paper offers an integrated perspective on existing modeling methodologies providing guidance for model selection based on the initial and boundary conditions. By comparing various models, one can better assess the uncertainties in predictions, particularly those related to the conceptual models. The review explores how these methodologies have significantly enhanced the fundamental understanding of how fractures respond to fluid injection and production, and improved predictive capabilities pertaining to coupled processes within fractured systems. It emphasizes the importance of utilizing advanced computational technologies and thoroughly considering fundamental theories and principles established through past experimental evidence and practical experience. The selection and calibration of model parameters should be based on typical ranges and applied to the specific conditions of applications. The challenges arising from inherent heterogeneity and uncertainties, nonlinear THM coupled processes, scale dependence, and computational limitations in representing field scale fractures are discussed. Realizing potential advances on computational capacity calls for methodical conceptualization, mathematical modeling, selection of numerical solution strategies, implementation, and calibration to foster simulation outcomes that intricately reflect the nuanced complexities of geological phenomena. Future research efforts should focus on innovative approaches to tackle the hurdles and advance the state-of-the-art in this critical field of study.

Coupling scheme↗