Search NASA⌕ Search

SEARCH · Search NASA

Results for “Reservoir flow model”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Leveraging operational information from wastewater injection wells to evaluate CO 2 injection performance for carbon storage applications in the Appalachian Basin

Abstract Geologic parameters, geophysical logging, injection testing, and operational metrics from wastewater injection wells were integrated to develop a preliminary design of a carbon storage facility in the Appalachian Basin. A scattered group of 10–20 commercial wastewater injection wells dispose off produced water from oil and gas wells in the region, utilizing a sequence of stacked deep saline formations for injection zones. These wastewater injection wells provide practical benchmarks for understanding the feasibility of carbon dioxide (CO 2 ) storage. Geologic models were developed based on characterization data from the wastewater injection wells. Reservoir simulations were calibrated according to injection testing and operational data from the wastewater injection wells. Long‐term operational data on injection flow rates and pressures measured in the wastewater injection wells were especially useful to evaluate the performance of carbon storage applications. The simulations were used to estimate injection pressures, radius of CO 2 saturation, and pressure response for industrial scale CO 2 storage applications. Results were also used to provide a design basis in terms of number of injection wells, well spacing, area of review, injection system components, monitoring plan, and CO 2 pipeline distribution system. The analysis demonstrates that there is sufficient injectivity in the deep saline formations in the west‐central Appalachian Basin to store commercial volumes of anthropogenic CO 2 . The geologic system appears suitable for supporting CO 2 injection rates of 0.5–1.0 million metric tons per year at injection pressures below formation fracture pressure in a single well. The long‐term operational data of wastewater injection wells within the study area suggested a lower permeability‐thickness values than indicated by initial reservoir tests. A workflow for developing realistic permeability values for input into reservoir simulations is presented. © 2020 Society of Chemical Industry and John Wiley & Sons, Ltd.

Valluri, Manoj↗
DOE OSTI2021

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.

artificial intelligence↗
DOE OSTI2021

Subsurface Characterization and Machine Learning Predictions at Brady Hot Springs: Preprint

Subsurface data analysis, reservoir modeling, and machine learning (ML) techniques have been applied to the Brady Hot Springs (BHS) geothermal field in Nevada, USA to further characterize the subsurface and assist with optimizing reservoir management. Hundreds of reservoir simulations have been conducted in TETRAD-G and CMG STARS to explore different injection and production fluid flow rates and allocations and to develop a training data set for ML. This process included simulating the historical injection and production since 1979 and prediction of future performance through 2040. ML networks were created and trained using TensorFlow based on multilayer perceptron (MLP), long short-term memory (LSTM), and convolutional neural network (CNN) architectures. These networks took as input selected flow rates, injection temperatures, and historical field operation data and produced estimates of future production temperatures. This approach was first successfully tested on a simplified single fracture doublet system, followed by the application to the BHS reservoir. Using an initial BHS dataset with 37 simulated scenarios, the trained and validated network predicted the production temperature for 6 production wells with the mean absolute percentage error of less than 8%. In a complementary analysis effort, the principal component analysis applied to 13 BHS geological parameters revealed that vertical fracture permeability shows the strongest correlation with fault density and fault intersection density. A new BHS reservoir model was developed considering the fault intersection density as proxy for permeability. This new reservoir model helps to explore under-exploited zones in the reservoir. A data gathering plan to obtain additional subsurface data was developed; it includes temperature surveying for three idle injection wells, at which the reservoir simulations indicate high bottom-hole temperatures. The collected data assist with calibrating the reservoir model and may lead to converting these wells to producers to access under-exploited zones in the reservoir. Data gathering activities are planned for the first quarter of 2021.

40 EE - Geothermal Technologies Office (EE-4G)↗
DOE OSTI2022

CO 2 zonal injection rate allocation and plume extent evaluation through wellbore temperature analysis

Temperature analysis during a pause in injection operations, known as warmback analysis, has been used in the petroleum industry for evaluating the injection conformance and estimating the location of the flooded front in applications, such as waterflooding oil reservoirs. Here in this work, methods are introduced to extend the application of temperature warmback analysis to estimate the zonal CO 2 injection rate and zonal CO 2 plume extent during geologic CO 2 storage in a saline aquifer. First, novel analytical solutions are developed to model transient temperature in the aquifer during the injection and subsequent shut-in periods considering two-phase flow (gaseous CO 2 and aqueous brine) conditions in the aquifer. The solution involves a discretization of the aquifer into regions; the energy and mass conservation equations for the regions are solved simultaneously considering appropriate boundary conditions at the interfaces. Two solutions techniques are presented: multi-region and three-region solutions. Inverse models are developed accordingly to evaluate the injection profile and estimate the extent of the plume front in the reservoir during the injection period. The multi-region solution results in an inversion approach that requires regression analysis. However, the three-region formulation results in a simple graphical technique for inverse modeling. The analytical solutions are validated against a thermally coupled reservoir simulation tool using different synthetic cases for CO 2 injection in deep saline aquifers. The results of the developed solutions provide a good match with numerical results during forward and inverse modeling.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Advances and challenges in CO2 foam technologies for enhanced oil recovery in carbonate reservoirs

Utilization of surfactants for generating foam has the potential for fluid-mobility control during CO2 flooding processes, leading to improved oil recovery. This manuscript is a review of the most important aspects for the design of CO2-foams with mobility control in carbonate reservoirs, including recent advances in novel surfactants and analytical techniques for analysis of surfactant adsorption on carbonate minerals and their thermal stability. Several key parameters and properties regarding foam transport in porous media are reviewed such as, the minimum pressure gradient required for foam generation, the effects of the partition coefficient of a CO2-soluble surfactant on foam transport in rock cores, visualization of foam flow by either microfluidic or PET/CT imaging processes in heterogeneous or fractured systems and the effect of oil on foam transport under reservoir conditions. Also are discussed the most recent advances of two major foam modeling methods, including the semi-empirical STARS foam model and the population balance model, as well as CO2 foam pilot tests and the factors governing their successes.

Jian, Guoqing↗
DOE OSTI2021

Predicting Thermal Performance of an Enhanced Geothermal System From Tracer Tests in a Data Assimilation Framework

Abstract Predicting the thermal performance of an enhanced geothermal system (EGS) requires a comprehensive characterization of the underlying fracture flow patterns from practically available data such as tracer data. However, due to the inherent complexities of subsurface fractures and the generally insufficient geological/geophysical data, interpreting tracer data for fracture flow characterization and thermal prediction remains a challenging task. The present study aims to tackle the challenge by leveraging a data assimilation method to maximize the utilization of information inherently contained in tracer data, and meanwhile maintain the flexibility to handle various uncertainties. A tracer data interpretation framework was proposed with the following three components integrated: (a) We use principal component analysis (PCA) to reduce the dimensionality of model parameter space. (b) We use ES‐MDA (ensemble smoother with multiple data assimilation) to invert for fracture aperture/flow fields and obtain posterior model ensembles for uncertainty quantification. Various data types are assimilated jointly to improve the predictive ability of the posterior ensemble. (c) The inverted fracture aperture fields are then incorporated into reservoir models to predict thermal performance. We developed a field‐scale EGS model to verify the ability of the framework to characterize highly heterogeneous fracture aperture/flow fields and predicting thermal performance. We also applied the framework to a mesoscale field experiment to demonstrate its potential application in real‐world geothermal reservoirs. The results indicate that the proposed framework can effectively retrieve fracture flow information from tracer data for thermal prediction and uncertainty quantification, and thus provide informative guidance for EGS optimization and risk management.

15 GEOTHERMAL ENERGY↗
DOE OSTI2020

Inferred inflow forecast horizons guiding reservoir release decisions across the United States

Abstract. Medium- to long-range forecasts often guide reservoir release decisions to support water management objectives, including mitigating flood and drought risks. While there is a burgeoning field of science targeted at improving forecast products and associated decision support models, data describing how and when forecasts are applied in practice remain undeveloped. This lack of knowledge may prevent hydrological modelers from developing accurate reservoir release schemes for large-scale, distributed hydrology models that are increasingly used to assess the vulnerabilities of large regions to hydrological stress. We address this issue by estimating seasonally varying, regulated inflow forecast horizons used in the operations of more than 300 dams throughout the conterminous United States (CONUS). For each dam, we take actual forward observed inflows (perfect foresight) as a proxy for forecasted flows available to the operator and then identify for each week of the year the forward horizon that best explains the release decisions taken. Resulting “horizon curves” specify for each dam the inferred inflow forecast horizon as a function of the week of the water year. These curves are analyzed for strength of evidence for contribution of medium- to long-range forecasts in decision making. We use random forest classification to estimate that approximately 80 % of large dams and reservoirs in the US (1553±50 out of 1927 dams with at least 10 Mm3 storage capacity) adopt medium- to long-range inflow forecasts to inform release decisions during at least part of the water year. Long-range forecast horizons (more than 6 weeks ahead) are detected in the operations of reservoirs located in high-elevation regions of the western US, where snowpack information likely guides the release. A simulation exercise conducted on four key western US reservoirs indicates that forecast-informed models of reservoir operations may outperform models that neglect the horizon curve – including during flood and drought conditions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Vanishing adsorption limit of Riemann problem solutions for the polymer model

Here we examine the vanishing adsorption limit of solutions of Riemann problems for the Glimm–Isaacson model of chemical flooding of a petroleum reservoir. A contact discontinuity is deemed admissible if it is the limit of traveling waves or rarefaction waves for an augmented system that accounts for weak chemical adsorption onto the rock. We prove that this criterion justifies the admissibility criteria adopted previously by Keyfitz–Kranzer, Isaacson–Temple and de Souza–Marchesin, provided that the fractional flow function depends monotonically on chemical concentration. We also demonstrate that the adsorption criterion selects the undercompressive contact discontinuities required to solve the general Riemann problem in an example model with non-monotone dependence.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

An extended trajectory-mechanics approach for calculating two-phase flow paths

A technique originating in quantum dynamics is used to derive a trajectory-based, semi-analytical solution for two-phase flow. The partial differential equation governing the evolution of the aqueous phase is equivalent to a family of ordinary differential equations defined along a path through the porous medium. The trajectories may be found by solving the differential equations directly or by post-processing the output of a numerical solution to the full set of governing equations. The trajectories, which differ from conventional streamlines, are found to bend downward in response to gravitational forces. The curvature is more pronounced as the dip of the porous layer containing the flow increases. Subtle changes in the relative permeability curve can lead to significant variations in the trajectories. The ordinary differential equation for the trajectory provides an expression for the travel time along the path. The expression produces a semi-analytical approximation to the model parameter sensitivities, the partial derivatives of the travel times with respect to changes in the permeability model. The semi-analytical trajectory-based sensitivities generally agree with those computed using a numerical reservoir simulator and a perturbation approach. The sensitivities are useful in tomographic imaging algorithms designed to estimate the spatial variation in permeability within a porous medium using multiphase observations.

58 GEOSCIENCES↗
DOE OSTI2022

Voronoi Meshing to Accurately Capture Geological Structure in Subsurface Simulations

Mesh generation lies at the interface of geological modeling and reservoir simulation. Highly skewed or very small grid cells may be necessary to accurately capture the geometry of geological features, but the resulting poorly scaled or small grid cells can have a substantial negative impact on simulator accuracy and speed. One way to minimize numerical errors caused by gridding complex structures is to simulate on high-quality Voronoi meshes, which reduce grid orientation effects in fluid flow. This work presents a complete methodology to create Voronoi simulation grids, model fluid flow in complex geological systems, and visualize the results. A recently developed Voronoi meshing method that can automatically generate provably good unstructured meshes that conform to input surfaces creating closed volumes is used. Initially an analytical benchmark simulation is presented to validate the quality of the meshes and simulation results and demonstrate the superiority of simulation results using Voronoi meshes over flexed-hexahedral meshes on a domain with internal features. Next, meshes are created for test structures representing four of the most common geological features in the subsurface: layering, pinch-out, an interior lens that tapers to zero thickness on all sides and a fault with offset. Two benchmark flow simulations are run for each test structure. Finally, a realistic geological example for CO 2 injection into an anticline is simulated. Three realizations of the Voronoi mesh at the same resolution are generated for the simulations. Each mesh is highly refined near the injection wells and coarse in areas of less interest. These three meshes are used to model the CO 2 plume in the subsurface as it migrates to the top of the structure and then fills downward. Simulations on the meshes with randomly generated elements inside the input volumes each give slightly different fingering patterns for the viscous-unstable buoyant gas flow. The results presented in this work show a promising step towards utilizing fully automated Voronoi meshing for subsurface flow simulations in complex geology.

58 GEOSCIENCES↗
DOE OSTI2024

River–aquifer interactions enhancing evapotranspiration in a semiarid riparian zone: A modelling study

The hydrologic flows across the river–aquifer interface play an important role in groundwater dynamics and biogeochemical reactions within the subsurface; however, little is known about the effects of river–aquifer interactions on land surface processes. In this study, we developed a fully coupled three-dimensional (3D) land surface and subsurface model at a high resolution (~1 km) that accounts for high-frequency hydrologic exchange flow conditions to investigate how river–aquifer interactions modulate surface water budgets in the Upper Columbia-Priest Rapids watershed, a typical semiarid watershed located in the northwestern United States where river stage fluctuates in response to reservoir releases changing. Our results show that the spatiotemporal dynamics of river–aquifer interactions are highly heterogeneous, driven mainly by river-stage fluctuations. Adding 6.64 × 10 6 m 3 year –1 of water over the watershed from the river to groundwater owing to the lateral flow, river–aquifer interactions led to an increase in soil evaporation and transpiration supplied by higher soil moisture content, particularly in deeper subsurface. In a hypothetic future scenarios where a 5-m rise in river stage was assumed, the hydrologic flow exchange rates were intensified, resulting in higher surface water over the entire watershed. Overall, lateral flow induced by river–aquifer exchanges leads to an increase in evapotranspiration of ~75% in the historical period and of ~83% in the hypothetical future scenario. Finally, our study demonstrates the potential of coupled model as an effective tool for understanding river–aquifer–land surface interactions, and indicates that river–aquifer interactions fundamentally alter the water balance of the riparian zone for the semiarid watershed and will likely become more frequent and intense in the future under the effects of climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

On the Feasibility of Deep Geothermal Wells Using Numerical Reservoir Simulation

This study examines the geothermal energy extraction potential from the basement rock within the Denver–Julesburg Basin, focusing on the flow performance and heat extraction efficiency of different geothermal well configurations. It specifically compares U-shaped, V-shaped, inclined V-shaped, and pipe-in-pipe configurations against enhanced geothermal system setups. Through numerical modeling, we evaluated the thermal behavior of these systems under various operational scenarios and fracture conditions. The results suggest that while closed-loop systems offer moderate temperature increases, Enhanced geothermal system configurations show substantial potential for high-temperature extraction. This underscores the importance of evaluating well configurations in complex geological settings. The insights from this study aid in strategic geothermal energy planning and development, marking significant advancements in geothermal technology and setting a foundation for future explorations and optimizations.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Experiments and Modeling of Proppant Embedment and Fracture Conductivity for the Caney Shale, Oklahoma, USA

ABSTRACT: The ultimate aim of hydraulic fracturing is to have a long and conductive flow path that extends from the wellbore into the formation. The effective fracture length is part of a hydraulically propped fracture which contributes to production. The difficulty in achieving economical production targets from shale reservoirs is at the forefront in many exploration companies. Fracture conductivity loss is related to; proppant embedment under depletion, proppant crushing, damage as a result of fracturing fluid, fines migration and proppant-pack permeability-damage are some of the factors that contribute to production decline after hydraulic fracturing in shale reservoirs. The Caney Shale is a calcareous organic-rich mudrock. Various studies have investigated the effect that clay on shale well productivity, however, there is currently no literature on the Caney shale in relation to horizontal wells; all available literature exists in vertical wells as well as on formations of the Caney that are shallow in comparison to an emerging play which is twice the depth. In this paper we investigate stress-dependent fracture conductivity of proppant-filled fractures and proppant embedment in Caney shale through laboratory and modeling studies. API fracture conductivity tests were conducted using 2% KCl on five locations within the Caney shale that consisted of selecting three brittle(reservoir) zones and two ductile zones. Confining pressures range from 1,000 psi to 12,000 psi at 210°F. Conductivity, permeability as well as embedment were measured during the test. Our experimental results have confirmed that improved fracture conductivity is attributed to; proppant size, the increase in porosity of the proppant pack, closure pressure changes and the reduction in fracture conductivity are a function of many factors such as fracture closure stress. The findings from this study could help the stimulation design by providing new insights into the critical factors that are to be determined to facilitate the choice of proppants as well as fracturing fluids for long term production and recovery from shale reservoirs. 1 INTRODUCTION The development of low permeability formations, like shales, has been aided by hydraulic fracturing of horizontal wells (Radonjic et al., 2020). Hydraulic fracturing fluid is injected at a high pressure to induce tensile fractures that can link to and stimulate natural fractures (Katende et al., 2021a,b). Preserving adequate conductivity in hydraulic fractures over the life of the wells is required for economic production; nevertheless, conserving such conductivity can be difficult in some circumstances, particularly in soft, clay-rich formations (Wang et al., 2021). Proppant particles help to keep the fractures open when the pumping stops and the fracturing fluid returns to the wellbore, producing one or more propped hydraulic fractures of varying length, breadth, and height (Katende et al., 2021a). The proppant pack within the hydraulic fracture boosts well output by providing a greater permeability flowpath for hydrocarbons (Duenckel et al., 2016). Proppant in the fracture is under complicated stress conditions, and the interplay between the rock formation and the proppant pack has a significant impact on proppant-pack permeability (Karazincir et al., 2019). Proppant may be embedded (Katende et al., 2021a) in the rock or crushed into small pieces if the proppant size and strength characteristics are not specified appropriately, resulting in a loss in proppant-pack permeability and fracture aperture, and consequently a fall in well output.

Katende, A.↗
DOE OSTI2020

Three-dimensional fracture continuum characterization aided by surface time-domain electromagnetics and hydrogeophysical joint inversion—proof-of-concept

Efficient and safe production of hydraulically fractured reservoirs benefits from the prediction of their geometrical attributes. Geophysical methods have the potential to provide data that are sensitive to fracture geometries, alleviating the typically sparse nature of in situ reservoir observations. Moreover, surface-based methods can be logistically and economically attractive since they avoid operational interference with the injection well infrastructure. Here, we investigate the potential of the surface-based time-domain electromagnetic (EM) method. EM methods can play an important role owing to their sensitivity to injection-induced fluid property changes. Two other advantageous factors are the EM signal-enhancing effect of vertical steel-cased wells and the fact that injected proppants can be enhanced to produce a stronger electrical conductivity contrast with the reservoir’s connate fluid. Nevertheless, an optimal fracture characterization will no doubt require the integration of EM and reservoir injection and production data. We hence carry out our investigations within a hydrogeophysical parameter estimation framework where EM data and injection flow rates are combined in a fully coupled way. Given the interdisciplinary nature of coupled hydrogeophysical inverse modeling, we dedicate one section to laying out key aspects in a didactic manner.

58 GEOSCIENCES↗
DOE OSTI2021

Coordination and control – limits in standard representations of multi-reservoir operations in hydrological modeling

Abstract. Major multi-reservoir cascades represent a primary mechanism for dealing with hydrologic variability and extremes within institutionally complex river basins worldwide. These coordinated management processes fundamentally reshape water balance dynamics. Yet, multi-reservoir coordination processes have been largely ignored in the increasingly sophisticated representations of reservoir operations within large-scale hydrological models. The aim of this paper is twofold, namely (i) to provide evidence that the common modeling practice of parameterizing each reservoir in a cascade independently from the others is a significant approximation and (ii) to demonstrate potential unintended consequences of this independence approximation when simulating the dynamics of hydrological extremes in complex reservoir cascades. We explore these questions using the Water Balance Model, which features detailed representations of the human infrastructure coupled to the natural processes that shape water balance dynamics. It is applied to the Upper Snake River basin in the western US and its heavily regulated multi-reservoir cascade. We employ a time-varying sensitivity analysis that utilizes the method of Morris factor screening to explicitly track how the dominant release rule parameters evolve both along the cascade and in time according to seasonal high- and low-flow events. This enables us to address aim (i) by demonstrating how the progressive and cumulative dominance of upstream releases significantly dampens the ability of downstream reservoir rules' parameters to influence flow conditions. We address aim (ii) by comparing simulation results with observed reservoir operations during critical low-flow and high-flow events in the basin. Our time-varying parameter sensitivity analysis with the method of Morris clarifies how independent single-reservoir parameterizations and their tacit assumption of independence leads to reservoir release behaviors that generate artificial water shortages and flooding, whereas the observed coordinated cascade operations avoided these outcomes for the same events. To further explore the role of (non-)coordination in the large deviations from the observed operations, we use an offline multi-reservoir water balance model in which adding basic coordination mechanisms drawn from the observed emergency operations is sufficient to correct the deficiencies of the independently parameterized reservoir rules from the hydrological model. These results demonstrate the importance of understanding the state–space context in which reservoir releases occur and where operational coordination plays a crucial role in avoiding or mitigating water-related extremes. Understanding how major infrastructure is coordinated and controlled in major river basins is essential for properly assessing future flood and drought hazards in a changing world.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Hopping mediated transport between finite pools of redox proteins

Transport reactions in biology involve the flow of particles—electrons, ions, or molecules—between reservoirs. Here, we explore how electron transport between finite reservoirs depends on the nature of the reservoirs, including their size, occupancy, and interactions. We compare the transport kinetics produced by narrowband and wideband infinite reservoir models (described earlier) with a finite narrowband reservoir model. The transport between finite reservoirs is found to depend on both the initial charge distribution and the number of carriers present. Whether or not a steady-state transport regime is accessed prior to reaching the equilibrium charge distribution depends on these initial conditions.

10 SYNTHETIC FUELS↗
DOE OSTI2020

Big Sky Regional Carbon Sequestration Partnership (Phase III Final Scientific/Technical Report)

The Big Sky Carbon Sequestration Partnership (BSCSP) pursued a Phase III demonstration project at Kevin Dome in north central Montana. Kevin Dome covers approximately 700 square miles and is a naturally occurring CO 2 reservoir that is flanked by oil and gas fields. The carbon dioxide (CO 2 ) is in the upper Devonian Duperow (carbonate) formation and does not reach the spill point of the dome; therefore, the dome has potential as a CO 2 sequestration reservoir, a CO 2 supply, or as both if anthropogenic sources and enhanced oil recovery (EOR) operations are associated with the dome. Kevin Dome could potentially act as a buffer to continue accepting anthropogenic CO 2 when EOR flooding operations are interrupted or completed. The project objective was to produce one million tonnes of CO 2 from the gas cap of Kevin Dome, pipe it laterally, inject, and re-store it in the brine leg of the same formation to test the hub / buffer storage concept. This was to be accomplished by drilling up to five production wells, building a short pipeline and compression facilities, and drilling an injection well and several monitoring wells. BSCSP commenced outreach and site characterization activities including acquiring baseline data for near-surface insurance monitoring, acquiring 3-dimensional, 9-component surface seismic over the project area, drilling two test wells (one in the production area and one in the injection area), coring key intervals, and performing comprehensive logging. Well tests of those wells revealed two barriers to the project. The production (Danielson 33-17) well was perforated in multiple zones but failed to produce any significant CO 2 . This was despite being drilled in the near vicinity of a historic well that had produced 3700 MCF per day in a drill stem test. Modeling indicated that this was likely due to a phase change during production causing a temperature drop resulting in hydrate and/or water ice formation that clogged the formation. Tests of the injection zone (Wallewein 22-1) well indicated total dissolved solids (TDS) slightly below the EPA required 10,000 parts per million (ppm) for a Class VI underground injection control permit. While the project was initiated before Class VI rules were promulgated, and this was an experimental project (seemingly qualified for a Class V permit), the Environmental Protection Agency (EPA) indicated that injection would require a Class VI permit. The low salinity result was unexpected as contours plotted based on regional formation water quality data indicated an expected TDS above 20,000 ppm, and wells between the recharge zone and the Wallewein well tested above 10,000 ppm. Faced with the inability to obtain an injection permit, the demonstration project could not proceed. However, the project had generated valuable samples and data on a large natural analog including 32 sq. mi. of 3-D, 9-C seismic, 430 ft. of carbonate core covering seven different depositional environments taken from areas with, and without the presence of CO 2 , 30 ft. of core of two caprocks, a tight carbonate and an anhydrite, a full set of modern logs on both wells, and well tests. DOE decided to re-scope the project around completing studies utilizing this data. This report covers both the initial scope and the re-scope (Task / Section numbers preceded with an R). While the report covers a wide range of project activities, highlights of this work include: Development of a geostatic model using neural nets to match well logs to facies and using multi-waveform seismic to inform reservoir heterogeneity; Unique mechanical testing of permeability – stress relationship in two caprock materials; Development of full waveform inversion to generate a high resolution velocity model; Model development for dual permeability (fracture and matrix) systems to better account for matrix-matrix interactions; Joint seismic wave inversion (including the first quadr-joint inversion) exhibiting better imaging of a challenging reservoir zone in stiff rock; Core flow and core flood results on a reactive carbonate; and Innovative laboratory measurements of seismic response of fractured core as a function of fluid fill.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

A screening model for predicting injection well pressure buildup and plume extent in CO 2 geologic storage projects

We present a new screening model for predicting injection well pressure buildup and CO 2 plume migration for CO 2 geologic sequestration projects. The model requires only limited information and is quite accurate when compared to detailed simulation results. Such models can assist project developers during the early days of project planning (e.g., for 45Q related projects), and also help regulators perform simple checks against detailed numerical models. The screening model consists of two correlations: one for injectivity index in terms of the slope of the CO 2 fractional flow curve, and the second for total storage efficiency within the footprint of plume as a function of gravity number and slope of the CO 2 fractional flow curve. Using these two correlations, as well as a knowledge of some basic reservoir characteristics and estimates of fluid properties from standard correlations, the injection-well pressure buildup and CO 2 plume extent in the formation can be readily estimated. Finally, the new correlations show a good match with the results of the underlying simulations, and also provide good agreement with independent calculations for two example problems.

42 ENGINEERING↗