SEARCH · Search NASA
Results for “reservoir geomechanics”
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.
Geomechanical response due to nonisothermal fluid injection into a reservoir
The geomechanical response of a porous reservoir due to injection of fluid can result from a complex interplay between the changes in porepressure and temperature near the wellbore. As a result, predictions are usually made using either simplified analytical models, which may apply unrealistic assumptions in order to produce a tractable model, or detailed numerical simulations that can be computationally expensive.
Understanding poromechanical response of a biogenic coalbed methane reservoir
Abstract Biogenic coalbed methane (BCBM) reservoirs aim to produce methane from in situ coal deposits following microbial conversion of coal. Success of BCBM reservoirs requires economic methane production within an acceptable timeframe. The work reported here quantifies the findings of previously published qualitative work, where it was found that bioconversion induces strains in the pore, matrix and bulk scales. Using imaging and dynamic strain monitoring techniques, the bioconversion induced strain is quantified here. To understand the effect of these strains from a reservoir geomechanics perspective, a corresponding poromechanical model is developed. Furthermore, findings of imaging experiments are validated using core-flooding flow experiments. Finally, expected field-scale behavior of the permeability response of a BCBM operation is modeled and analyzed. The results of the study indicated that, for Illinois coals, bioconversion induced strains result in a decrease in fracture porosity, resulting in a detrimental permeability drop in excess of 60% during bioconversion, which festers itself exponentially throughout its producing life. Results indicate that reservoirs with high initial permeability that will support higher Darcian flowrates, would be better suited for coal bioconversion, thereby providing a site-selection criteria for BCBM operations.
Convergence analysis of single rate and multirate fixed stress split iterative coupling schemes in heterogeneous poroelastic media
Recently, the accurate modeling of flow–structure interactions has gained more attention and importance for both petroleum and environmental engineering applications. Of particular interest is the coupling between subsurface flow and reservoir geomechanics. Different single rate and multirate iterative and explicit coupling schemes have been proposed and analyzed in the past. In addition, Banach fixed point contraction results were obtained for iterative coupling schemes, and conditionally stable results were obtained for explicit coupling schemes. In this work, we will consider the mathematical analysis of the single rate and multirate fixed stress split iterative coupling schemes for spatially heterogeneous poroelastic media. We will re–establish the contractivity for both schemes in the localized case, and we will show that heterogeneities come at the expense of imposing more restricted conditions on the number of fine flow time steps that can be taken within one coarse mechanics time step in the multirate case. Our mathematical analysis is supplemented by numerical simulations validating our derived upper bounds. Finally, to the best of our knowledge, this is the first rigorous mathematical analysis of the multirate fixed–stress split iterative coupling scheme in heterogeneous poroelastic media.
Joint state-parameter estimation for the reduced fracture model via the united filter
Here, in this paper, we introduce an effective United Filter method for jointly estimating the solution state and physical parameters in flow and transport problems within fractured porous media. Fluid flow and transport in fractured porous media are critical in subsurface hydrology, geophysics, and reservoir geomechanics. Reduced fracture models, which represent fractures as lower-dimensional interfaces, enable efficient multi-scale simulations. However, reduced fracture models also face accuracy challenges due to modeling errors and uncertainties in physical parameters such as permeability and fracture geometry. To address these challenges, we propose a United Filter method, which integrates the Ensemble Score Filter (EnSF) for state estimation with the Direct Filter for parameter estimation. EnSF, based on a score-based diffusion model framework, produces ensemble representations of the state distribution without deep learning. Meanwhile, the Direct Filter, a recursive Bayesian inference method, estimates parameters directly from state observations. The United Filter combines these methods iteratively: EnSF estimates are used to refine parameter values, which are then fed back to improve state estimation. Numerical experiments demonstrate that the United Filter method surpasses the state-of-the-art Augmented Ensemble Kalman Filter, delivering more accurate state and parameter estimation for reduced fracture models. This framework also provides a robust and efficient solution for PDE-constrained inverse problems with uncertainties and sparse observations.
Analysis of Completion Design Impact on Cluster Efficiency and Pressure-Based Well Communication in HFTS-2 Delaware Basin
The Hydraulic Fracturing Test Site 2 (HFTS-2) is a joint industry project in the Delaware basin to advance hydraulic fracturing understanding and improve productivity in shale reservoirs. The project integrates multi-disciplinary approaches to evaluate different completion designs, well spacing, inter-well communication, and stimulated rock volume, among other factors. This paper focused on two major areas related to hydraulic fracture performance. First, an analysis of different completion designs on cluster efficiency based on near well Distributed Acoustic Sensing (DAS). Second, an evaluation of well and completion designs on inter-well fracture driven interactions (FDIs) based on downhole pressure monitoring. Fluid/sand distribution and cluster efficiency analyses were based on near wellbore DAS data collected from two adjacent horizontal wells completed in two different landing zones in the Wolfcamp formation. These wells had different completion designs aiming to evaluate the effect of normal vs. extended stage lengths, perforation hole tapering and limited entry. Standard deviation from ideal fluid/sand distribution and waterfall plots were used to evaluate cluster efficiency for each design and stage. Inter-well FDIs analysis was conducted among the horizontal wells and a vertical monitor well. One horizontal well served as the monitor well while the other horizontal well was being treated. The vertical well was instrumented with downhole pressure and temperature gauges to aid monitoring fracture height growth. The pressure response during and after fracturing was characterized based on maximum pressure increase value and slope. Pressure response vs. FDIs trigger factors such distance, cluster efficiency and stage fluid volume were also analyzed. Based on the different completion and perforation designs tested, DAS analysis suggests that limited entry design worked best. Extreme limited entry showed the potential of high perforation erosion and reduced cluster efficiency. Here, the limited entry and tapered perforation design demonstrated potential to improve the cluster efficiency for extended stage lengths. Pressure monitoring across formation units proved to be critical to understand fracture interactions and fracture vertical growth. Pressure communication across different formation units during hydraulic fracturing operation indicate fractures grew upwards during Wolfcamp wells fracturing. However, this pressure communication dissipated over time. High intensity FDIs were recorded when the frac stages Downloaded from http://onepetro.org/URTECONF/proceedings-pdf/21URTC/1-21URTC/D011S018R001/2477607/urtec-2021-5289-ms.pdf/1 by Carol Worster on 28 February 2022 URTeC 5289 were closer to the pressure gauge location in the monitor wells. Some of these stages that produced high intensity FDIs also had high fluid volume per cluster and low cluster efficiency. The multi-disciplinary and high-quality data collected from HFTS-2 helped to further understand why completion approaches such as limited entry and tapered perforation design are successful in improving cluster efficiency. The DAS data combined with downhole high-resolution pressure measurements also helped to quantify the effect of lower cluster efficiency data on the incidence and intensity of FDIs.
Developing Upscaling Approach for Swarming Hydraulic Fractures Observed at Hydraulic Fracturing Test Site through Multiscale Simulations
We show that this work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.
Field Experiments and Model Validation: The EGS Collab Project
Implementing enhanced geothermal systems (EGS) will require improvements in understanding stimulation of crystalline rock to create appropriate flow pathways, and the ability to effectively simulate both the stimulation and the flow and transport processes in the resulting fracture network. The US Department of Energy (DOE) is addressing these and other challenges at multiple scales. The EGS Collab project, addressed here, is performing tests and modeling at the 10 m scale. The FORGE project is performing tests at the full reservoir scale. The EGS Collab team created an underground testbed at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, at a depth of approximately 1.5 km to examine hydraulic fracturing (Experiment 1). We are currently building a second testbed at SURF at a depth of about 1.25 km aimed at investigating shear stimulation (Experiment 2). In Experiment 1, we drilled eight boreholes in a well-characterized phyllite and installed geophysical sensors in six of them to create a well-instrumented testbed to allow careful monitoring of stimulation events and flow tests. Numerical simulation was used to answer key experimental design questions, to forecast fracture propagation trajectories and extents, and to analyze and evaluate results both in near-real-time and in detailed process studies. Stimulations performed in this testbed allowed quantification of processes occurring during stimulation and the examination of dynamic flow occurrences. Long-term ambient temperature and chilled water flow tests were performed in addition to many tracer tests to examine system behavior. Our second testbed, targeted at shear stimulation, is currently being built at the SURF Facility at a depth of about 1.25 km in amphibolite under a different set of stress and fracture conditions than Experiment 1.
Coupled geochemical-geomechanical processes in CO2 sequestration reservoirs in Southeast, US
Coupled geochemical-geomechanical processes in CO2 sequestration reservoirs in Southeast, US
Coupled geochemical-geomechanical processes in CO2 sequestration reservoirs in Southeast, US
Coupled geochemical-geomechanical processes in CO2 sequestration reservoirs in Southeast, US
Geomechanical response of CO2 sequestration reservoirs due to different carbonate dissolution patterns
Geomechanical response of CO2 sequestration reservoirs due to different carbonate dissolution patterns
Coupled geochemical-geomechanical alteration of sandstone in CO2 sequestration reservoirs.
Coupled geochemical-geomechanical alteration of sandstone in CO2 sequestration reservoirs.
Validation of strongly coupled geomechanics and gas hydrate reservoir simulation with multiscale laboratory tests
In this work, we validate a coupled flow-geomechanics simulator for gas hydrate deposits, named T+M AM , performing two meter-scale laboratory experiments of gas hydrates for production by depressurization, replicating the gas hydrate deposit in the Ulleung Basin, East Sea, South Korea. The first experiment with a sand-only specimen is a 1D 1 m-scale depressurization test based on the excess gas method, which represents the grain coating hydrate growth. On the other hand, the second is a 3D 1.5 m-scale test with the excess water method for a sand-mud alternating layer system, representing the pore filling hydrate growth. We measure production and displacement at the top with different depressurization levels. In particular, the 3D test exhibits high coupling strength of substantial deformation induced by incompressibility of water and high deformability of the specimen. For validation, we match pressure, flow rate, and displacement between the experimental data and numerical results. Thus, we identify that T+M AM is a reliable simulator, which can be applied to fields in both permafrost and deep oceanic hydrate deposits of strongly coupled flow and geomechanics systems. This validation also implies that other coupled simulators based on the same coupling formulation as T+M AM can be validated when individual flow and geomechanics simulators are stable and reliable.
Observations and Modeling of Fiber-Optics Strain on Hydraulic Fracture Height Growth in HFTS-2
Understanding fracture height growth can be of great significance to optimizing field development and improving recovery. The Hydraulic Fracturing Test Site 2 (HFTS-2) has provided a unique opportunity and an advanced dataset to allow us to observe and understand fracture geometries rigorously. Low frequency distributed acoustic sensing (LF-DAS) data from a vertical well in HFTS-2 showed three key observations: (i) excessive upward height growth (>1000 ft) and limited downward growth of the hydraulic fractures during pumping, (ii) considerable additional upward fracture height growth (~300 ft) after well shut in, and (iii) very complex LF-DAS strain rate patterns for a small fiber-to-stage offset. Advanced geomechanical modeling was performed to simulate the hydraulic fracture propagation and the resulting strain responses in the vertical direction. The modeling results demonstrated asymmetric upward and downward fracture height growths as observed in HFTS-2 with a similar upward height growth rate. Simulated waterfall plots of vertical strain rate showed distinct patterns for different fiber-to-fracture distances. The upward-growing fracture tip can be clearly identified by the interfaces between compressing and extending zones. It was also found that the complex strain rate patterns observed in HFTS-2 for small fiber-stage offsets were not caused by the mechanical layering but possibly result from the simultaneous propagation of multiple hydraulic fractures at different rates. Furthermore, the simulation results improved the understating of the HFTS-2 LF-DAS data and the simulated strain rate patterns could also serve as templates for fracture height interpretation from LF-DAS data in future.
Method to account for natural fracture induced elastic anisotropy in geomechanical characterization of shale gas reservoirs
Shale has been usually recognized as a transverse isotropic (TI) medium in conventional geomechanical log interpretation due to its laminated nature. However, when natural fractures exist in the shale rock, additional elastic anisotropy is introduced, converting laminated Shale to an orthorhombic (OB) medium. Previous studies illustrate that neglecting the natural fracture induced anisotropy in shale geomechanical log interpretation could lead to inaccurate evaluations of elastic moduli and in-situ stresses. In this paper, a new method is developed to account for the natural fracture induced anisotropy in geomechanical log interpretation based upon the TI acoustic model developed by the author and a characterization technique of elastic wave anisotropy (Sayers, 1991). The new OB model incorporates the four acoustic log data inputs and five modeling constraints in a nonlinear optimization algorithm to solve for the nine independent stiffness coefficients of an OB rock, and further to solve for the geomechanical properties and in-situ stress profiles in an OB formation. The new method was validated with a Marcellus Gas Shale field case. Both the new OB model and the conventional TI model were applied to interpret the minimum horizontal stress profile for the same formation. By comparing the results, the OB model is more robust than the TI from two aspects. First, the average stress magnitude predicted by the OB model is closer to the one measured by the Diagnostic Fracture Injection Test (DFIT). Second, the OB model predicts a more obvious stress barrier between the lower Marcellus and upper Onondaga Limestone than the TI model does. Finally, the predicted stress barrier is consistent with the observation of the microseismic events of a horizontal well drilled and completed nearby, which reveals that no hydraulic fracture propagates downward through the bottom boundary of Marcellus Shale into the underlying Onondaga Limestone.
Surrogate models for development of unconventional shale reservoirs by an integrated numerical approach of hydraulic fracturing, flow and geomechanics, and machine learning
We develop well-completion surrogate models by taking an integrated workflow of hydraulic fracturing, flow, geomechanics, and machine learning simulation. There are three steps in the proposed workflow. First, history-matching processes are conducted with the field data including pumping and production data for characterization. Second, full-physics simulation is performed with various parameters of the field development (e.g., cluster spacing, clusters per stage, pumping rates and times, amount of proppant, and well spacing) to generate multiple simulation results by changing the parameters of the completion design with well-known hydraulic fracturing, reservoir, geomechanics simulators to calculate fracture geometry, reservoir depressurization, induced stress changes. The workflow is demonstrated over a field in the Southern Midland Basin. Here, we take two completion scenarios: a single well case followed by a multi-well case. Finally, a Long Short-Term Memory (LSTM) machine learning algorithm is employed to create surrogate models that can replicate the full-physics simulation results. Furthermore, results show that the trained models applied in the single well and multi-well cases for a particular geological system can provide good accuracy close to those provided by full-physics simulations. Specifically, the site-specific surrogate models can predict fracture parameters (length, height, and surface area) and cumulative production accurately with computational efficiency, suggesting our proposed workflow can be used as a pragmatic tool for expediting the well completion optimization process.
Mechanisms for Microseismicity Occurrence Due to CO 2 Injection at Decatur, Illinois: A Coupled Multiphase Flow and Geomechanics Perspective
Here, we numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO 2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologically consistent model that honors key stratigraphic horizons and 3D fault surfaces interpreted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO 2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO 2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection interval to the basement is essential to explain the destabilization of the regions where microseismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO 2 sequestration operations.
Geomechanical analysis of the geothermal reservoir at San Emidio, Nevada
Not Available