SEARCH · Search NASA

Results for “Reservoir Modeling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

At least 235 records · Page 13

All sources NASA NTRS DOE OSTI

DOE OSTI2025

Modeling supercritical CO2 injection induced rupture of a minor fault embedded in a poroelastic layered reservoir-caprock system

CO2 injection for geologic carbon sequestration involves hydromechanical processes that lead to changes in fluid pressure and stresses that can activate existing faults. This paper presents a new method and workflow of modeling fault activation considering more complex three-dimensional geometry of natural faults using the TOUGH-FLAC multiphase fluid flow and geomechanical simulator. In this method and workflow, FLAC3D mechanical interfaces and TOUGH3 finite volume elements are discretized using computer aided design and gridding software along with a tailored mesh translation routine. The method and workflow are demonstrated with a model of a curved minor fault embedded in a poro-elastic layered reservoir-caprock system. The model is used for a comprehensive sensitivity analysis of fault responses to fault length, injection mass rate, injection schedule, well-fault distance, and well locations versus fault location. Four metrics (CO2 plume, shear state of fault, pressure and stress path at fault monitoring points) are selected to assess CO2 migration, pressure change, and the reactivation of faults. The results reveal that CO2 can bypass around the tip of the minor impermeable fault, building up pressure and poro-elastic stress on both sides that tends to impede fault rupture. Our study shows the benefit of carefully designing the injection to achieve the targeted final storage volume, starting at a relatively low rate for considerable time, and then ramping up the injection rate to the full rate of injection. The initial low injection has two distinct benefits: (1) it allows for the formation of an extensive CO2 plume with a much higher mobility through a low viscosity that will result in a lower pressure for a given injection rate, and (2) it allows for gradual build-up of horizontal poro-elastic stress within the reservoir that will tend to impede activation of steeply dipping faults. The injection scenario starting at a low injection rate, denoted here as conservative injection, can significantly reduce the risk of fault activation as high fluid mobility and reservoir strengthening poro-elastic stress has been established long before reaching the peak injection rates. Moreover, simultaneous injection in two injection wells on both sides of fault can provide further reservoir strengthening through poro-elastic stress buildup acting on a fault under normal faulting stress regime. The findings presented in the paper can provide practical and effective guidance on long-term, safe, and reliable geological CO2 storage.

Cao, Meng↗

DOE OSTI2022

Optimization of Water-Alternating-CO2 Injection Field Operations Using a Machine-Learning-Assisted Workflow

Summary This paper will present a robust workflow to address multiobjective optimization (MOO) of carbon dioxide (CO2)-enhanced oil recovery (EOR)-sequestration projects with a large number of operational control parameters. Farnsworth unit (FWU) field, a mature oil reservoir undergoing CO2 alternating water injection (CO2-WAG) EOR, will be used as a field case to validate the proposed optimization protocol. The expected outcome of this work would be a repository of Pareto-optimal solutions of multiple objective functions, including oil recovery, carbon storage volume, and project economics. FWU’s numerical model is used to demonstrate the proposed optimization workflow. Because using MOO requires computationally intensive procedures, machine-learning-based proxies are introduced to substitute for the high-fidelity model, thus reducing the total computation overhead. The vector machine regression combined with the Gaussian kernel (Gaussian-SVR) is used to construct proxies. An iterative self-adjusting process prepares the training knowledge base to develop robust proxies and minimizes computational time. The proxies’ hyperparameters will be optimally designed using Bayesian optimization to achieve better generalization performance. Trained proxies will be coupled with multiobjective particle swarm Optimization (MOPSO) protocol to construct the Pareto-front solution repository. The outcomes of this workflow will be a repository containing Pareto-optimal solutions of multiple objectives considered in the CO2-WAG project. The proposed optimization workflow will be compared with another established methodology using a multilayer neural network (MLNN) to validate its feasibility in handling MOO with a large number of parameters to control. Optimization parameters used include operational variables that might be used to control the CO2-WAG process, such as the duration of the water/gas injection period, producer bottomhole pressure (BHP) control, and water injection rate of each well included in the numerical model. It is proved that the workflow coupling Gaussian-SVR proxies and the iterative self-adjusting protocol is more computationally efficient. The MOO process is made more rapid by squeezing the size of the required training knowledge base while maintaining the high accuracy of the optimized results. The outcomes of the optimization study show promising results in successfully establishing the solution repository considering multiple objective functions. Results are also verified by validating the Pareto fronts with simulation results using obtained optimized control parameters. The outcome from this work could provide field operators an opportunity to design a CO2-WAG project using as many inputs as possible from the reservoir models. The proposed work introduces a novel concept that couples Gaussian-SVR proxies with a self-adjusting protocol to increase the computational efficiency of the proposed workflow and to guarantee the high accuracy of the obtained optimized results. More importantly, the workflow can optimize a large number of control parameters used in a complex CO2-WAG process, which greatly extends its utility in solving large-scale MOO problems in various projects with similar desired outcomes.

Energy & Fuels↗

DOE OSTI2022

Deep learning-accelerated 3D carbon storage reservoir pressure forecasting based on data assimilation using surface displacement from InSAR

Fast forecasting of the reservoir pressure distribution during geologic carbon storage (GCS) by assimilating monitoring data is a challenging problem. Due to high drilling cost, GCS projects usually have spatially sparse measurements from few wells, leading to high uncertainties in reservoir pressure prediction. To address this challenge, we use low-cost Interferometric Synthetic-Aperture Radar (InSAR) data as monitoring data to infer reservoir pressure build up. We develop a deep learning-accelerated workflow to assimilate surface displacement maps interpreted from InSAR and to forecast dynamic reservoir pressure. Employing an Ensemble Smoother Multiple Data Assimilation (ES-MDA) framework, the workflow updates three-dimensional (3D) geologic properties and predicts reservoir pressure with quantified uncertainties. We use a synthetic commercial-scale GCS model with bimodally distributed permeability and porosity to demonstrate the efficacy of the workflow. A two-step CNN-PCA approach is employed to parameterize the bimodal fields. The computational efficiency of the workflow is boosted by two residual U-Net based surrogate models for surface displacement and reservoir pressure predictions, respectively. The workflow can complete data assimilation and reservoir pressure forecasting in half an hour on a personal computer.

25 ENERGY STORAGE↗

DOE OSTI2021

Selective penetration behavior of microgels in superpermeable channels and reservoir matrices

Gel treatment is an effective way to attack excessive water production in many mature oilfields around the world. Selective penetration is desired for successful gel treatments. That is, gel materials should easily penetrate the target zones (i.e., channeling features such as superpermeable channels) without entering/damaging the nontarget zones (i.e., reservoir matrices or oil zones). This study revealed that presence of threshold penetration pressure (ΔP th ) was responsible for selective penetration behavior of tested microgels. The concept of ΔP th was utilized to figure out favorable working conditions for effective gel treatments. Microgel dispersions were injected into superpermeable (super-k) sandpacks (mimicking super-k channels in reservoirs, 60–221 darcies), heterogeneous models with super-k channels (79–230 darcies), and sandstone cores (mimicking reservoir matrices, 50–5000 md). The results demonstrated that a minimum differential driving pressure (i.e., threshold penetration pressure, ΔP th ) was required to push microgel particles to penetrate channels or matrices. The critical penetration behavior was closely related to the particle/pore size ratio. Low ΔP th at smaller particle/pore ratios was beneficial to allow easy penetration of gel materials into the channeling zones. On the contrary, high ΔPth at larger particle/pore ratios was desirable to prevent gel materials from massively invading and damaging the matrices. Instead, the gel particles accumulated at the inlet surface, and a gel cake was gradually formed. The cake further prevented the invasion of the gels. The cake could be removed by chemical breakers to resume the injectivity/productivity of the matrices. Correlations were developed to describe the relationship between ΔP th and particle/pore ratio. A distinct transition was identified at the particle/pore ratio of about 3. This research could help identify the favorable conditions to achieve successful gel treatments. In an effective conformance treatment, the particle/pore ratio in the channel should be sufficiently low to allow easy penetration of gel materials into the channel (e.g., particle/pore ratio<2 in this study). Meanwhile, the particle/pore ratio in the matrix should be large enough to support a high ΔP th and thus prevent massive gel invasion into the matrix. This study advances the current pore scale studies (a single particle passing through a single channel) to Darcy-scale characterization.

02 PETROLEUM↗

DOE OSTI2023

Hydropower Energy Storage Capacity (HESC) Dataset

The Hydropower Energy Storage Capacity (HESC) Dataset catalogs characteristics that are relevant to evaluating reservoir storage and estimates of energy storage capacity based on varying levels of detail. Hydropower dams and reservoirs were included based on information from the National Inventory of Dams (NID; USACE, 2021) and Global Reservoir and Dam (GRanD v1.3) and Existing Hydropower Assets datasets. These data provide a foundation for understanding available resources at existing hydropower facilities and their potential to provide storage of energy and more flexible generation. Estimates of energy storage capacity include: • Level 1 – nominal energy storage capacity based on maximum storage capacities and hydraulic head • Level 2 – nominal energy storage capacity based on historical models or observations of reservoir volume and hydraulic head. These estimates are provided based on capacity from the entire historical period as well as monthly values. • Level 3 – modeled energy generation based on volume-elevation relationships, historical storage, observed/modeled inflows, and hydraulic capacity of turbines and calculated both as overall and on a monthly basis. • Level 4 – modeled energy generation incorporating information from Level 3 and operational constraints. For facilities where installed capacity is known, there are also estimates for discharge duration (the length of time when a facility could provide generation at a given capacity).

13 HYDRO ENERGY↗

DOE OSTI2021

Kimberlina 1.2 Velocity Models and Seismic Data

Kimberlina 1.2 Velocity model and synthetic seismic data, produced in collaboration of teams at the National Energy Technology Laboratory, Los Alamos National Laboratory, and Lawrence Livermore National Laboratory through the National Risk Assessment Partnership. Data is associated with the following publication: Zheng Zhou, Youzuo Lin, Zhongping Zhang, Yue Wu, Zan Wang, Robert Dilmore, and George Guthrie, "A Data-Driven CO2 Leakage Detection Using Seismic Data and Spatial-Temporal Densely Connected Convolutional Neural Networks," International Journal of Greenhouse Gas Control, Vol 90, 2019. The Kimberlina 1.2 Velocity models were produced by Zan Wang, Robert Dilmore, William Harbert, and Lianjie Huang at NETL. The following citations are directly related to the creation of the velocity models: Wang, Z. Harbert, W., Dilmore, R., Huang, L. Modeling of time-lapse seismic monitoring using CO2 leakage simulations for a model CO2 storage site with realistic geology: Application in assessment of early leak-detection capabilities. International Journal of Greenhouse Gas Control. V. 76, September 2018, Pages 39-52. https://doi.org/10.1016/j.ijggc.2018.06.011 Wang, Z., Dilmore, R., Harbert, W. Inferring CO2 saturation from synthetic surface seismic and downhole monitoring data using machine learning for leakage detection at CO2 sequestration sites. International Journal of Greenhouse Gas Control, V. 100, September 2020. https://doi.org/10.1016/j.ijggc.2020.103115 The velocity models were built based on the Kimberlina 1.2 aquifer impact data which is associated with the following publications: Buscheck, T.A., Mansoor, K., Yang, X., Wainwright, H., and Carroll, S. (2019). Downhole pressure and chemical monitoring for CO2 and brine leak detection in aquifers above a CO2 storage reservoir. International Journal of Greenhouse Gas Control. 91. 102812. 10.1016/j.ijggc.2019.102812. Xianjin Yang, Thomas A. Buscheck, Kayyum Mansoor, Zan Wang, Kai Gao, Lianjie Huang, Delphine Appriou, Susan A. Carroll, Assessment of geophysical monitoring methods for detection of brine and CO2 leakage in drinking water aquifers, International Journal of Greenhouse Gas Control, Volume 90, 2019, 102803, ISSN 1750-5836, https://doi.org/10.1016/j.ijggc.2019.102803 The synthetic seismic data was produced by Youzuo Lin and team at LANL, and are associated with the following citations: Jordan, P. D., and J. L. Wagoner. Characterizing Construction of Existing Wells to a CO2 Storage Target: The Kimberlina Site, California. Zheng Zhou, Youzuo Lin, Zhongping Zhang, Yue Wu, Zan Wang, Robert Dilmore, and George Guthrie, "A Data-Driven CO2 Leakage Detection Using Seismic Data and Spatial-Temporal Densely Connected Convolutional Neural Networks," International Journal of Greenhouse Gas Control, Vol 90, 2019.

CO2 Leakage↗

DOE OSTI2020

Utah FORGE: Discrete Fracture Network (DFN) Data

The FORGE team is making these fracture models available to researchers wanting a set of natural fractures in the FORGE reservoir for use in their own modeling work. They have been used to predict stimulation distances during hydraulic stimulation at the open toe section of well 16A(78)-32. These fracture sets are fully stochastic and do not contain the deterministic set that matches the pilot well 58-32 FMI data. Well 58-32 has been completed and 16A(78)-32 is to be drilled as part of Phase 3. The original .fab files are not included due to redundancy. The *.fabgz data for the 800m and 1200m depth areas are in the native FracMan format and have been compressed using Gzip. Filtered data for the 800m depth area includes .csv spreadsheets, native FracMan (.fab), and GOCAD (.ts) files that are in a compressed zip format. The file titled "SGW 2020 Finnila and Podgorney DFN fracture files on GDR.pdf" is a description of the data and should be reviewed prior to data use.

15 GEOTHERMAL ENERGY↗

DOE OSTI2021

Utah FORGE Well 16A(78)-32 Simplified Discrete Fracture Network Data

15 GEOTHERMAL ENERGY↗

DOE OSTI2024

Utah FORGE 5-2615: Determination and Analysis of Thermo-poromechanical Response of Fractured Rock - 2024 Annual Workshop Presentation

This is a presentation on the Determination and Modeling-Informed Analysis of Thermo-poromechanical Response of Fractured Rock for Application to FORGE by the University of Oklahoma, presented by Ahmad Ghassemi. This video presentation discusses how to improve understanding and control of coupled thermoporomechanical (or thermo-hydro-mechanical-THM) processes in reservoir development, and to demonstrate its role in interpretations of the fracture closure pressure. This presentation was featured in the Utah FORGE R&D Annual Workshop on August 13-15, 2024.

15 GEOTHERMAL ENERGY↗

DOE OSTI2023

Simple Statistical Models for Predicting Overpressure Due to CO2 and Low-Salinity Waste-Fluid Injection into Deep Saline Formations

Deep saline aquifers have been used for waste-fluid disposal for decades and are the proposed targets for large-scale CO2 storage to mitigate CO2 concentration in the atmosphere. Due to relatively limited experience with CO2 injection in deep saline formations and given that the injection targets for CO2 sometimes are the same as waste-fluid disposal formations, it could be beneficial to model and compare both practices and learn from the waste-fluid disposal industry. In this paper, we model CO2 injection in the Patterson Field, which has been proposed as a site for storage of 50 Mt of industrial CO2 over 25 years. We propose general models that quickly screen the reservoir properties and calculate pressure changes near and far from the injection wellbore, accounting for variable reservoir properties. The reservoir properties we investigated were rock compressibility, injection rate, vertical-to-horizontal permeability ratio, average reservoir permeability and porosity, reservoir temperature and pressure, and the injectant total dissolved solids (TDS) in cases of waste-fluid injection. We used experimental design to select and perform simulation runs, performed a sensitivity analysis to identify the important variables on pressure build-up, and then fit a regression model to the simulation runs to obtain simple proxy models for changes in average reservoir pressure and bottomhole pressure. The CO2 injection created more pressure compared to saline waste-fluids, when similar mass was injected. However, we found a more significant pressure buildup at the caprock-reservoir interface and lower pressure buildup at the bottom of the reservoir when injecting CO2 compared with waste-fluid injection.

Ansari, Esmail↗

DOE OSTI2021

Hydropower Energy Storage Capacity Dataset

The Hydropower Energy Storage Capacity (HESC) Dataset catalogues estimates of nominal energy storage capacity based on varying levels of detail. Dams and reservoirs selected were selected based on those reported in the National Inventory of Dams (NID 2019) and/or the Global Reservoir and Dam (GRanD v1.3) datasets. These data provide a foundation for understanding available resources at existing hydropower facilities and their potential to provide storage of energy and more flexible generation. Current estimates include Level 1 (based on maximum storage capacities and hydraulic head) and Level 2 (based on historical models or observations of reservoir volume and hydraulic head). For facilities where installed capacity is known, there are also estimates for discharge duration or the length of time when a facility could provide generation at a given capacity. Essential information used to calculate the energy storage capacity and discharge duration (volume, hydraulic head, and details about the sources or records used to obtain those parameters) and summaries of historical generation (for context) are also included.

13 HYDRO ENERGY↗

DOE OSTI2020

Assessment of Enhanced Oil Recovery Using Carbon Dioxide in Michigan Basin Silurian Pinnacle Reefs

The Michigan team at Western Michigan University has provided technical analysis and support of the enhanced oil recovery (EOR) using carbon dioxide (CO2) in the Northern Michigan Silurian Pinnacle reefs. This report discussed the relevant data about reservoir properties, including lithologic and depositional characteristics of the reservoir formations and porosity and fluid flow properties through different compartments in the reservoir. The team has also developed models that represent the geological and physical characteristics of the reservoir and seal system.

EOR↗

DOE OSTI2020

Assessment of Enhanced Oil Recovery Using Carbon Dioxide in Michigan Basin Silurian Pinnacle Reefs

The Midwest Regional Carbon Sequestration Partnership (MRCSP) has incorporated the work of geologic research teams (Geoteams) in its regional characterization, project planning and carbon dioxide (CO 2 ) injection implementation work since the partnership was established by the U.S. Department of Energy (DOE) in 2003. Over this 16-year period, the cohort of Geoteams has grown from five to ten states and has contributed to the characterization of geologic sequestration opportunities, refinement of reservoir and seal data, and supported injection efforts through both predictive and post-injection assessments. The Michigan team at Western Michigan University has provided technical analysis and support of the enhanced oil recovery (EOR) using carbon dioxide (CO 2 ) in the Northern Michigan Silurian Pinnacle reefs. With Core Energy, LLC as the industry partner, this project has assessed the potential for EOR in known fields in the Michigan Pinnacle reefs and the capacity for CO2 storage in the reef reservoirs as they are depleted and after reservoir abandonment. This project managed by Core Energy, LLC, involves the only commercial CO 2 Enhanced Oil Recovery floods within the MRCSP region. The Silurian pinnacle reef trend shows significant potential for CO 2 -EOR: this play is relatively new, with approximately 60% of the wells in the largest Niagaran fields still accessible. Most wells have modern logs and there are many cores, core analyses, pressure tests and produced water analyses available for study from these fields. These carbonate reservoirs in Michigan have produced sufficient oil to warrant evaluation for additional recovery efforts. The role of the Michigan team was to collect and analyze all relevant data about reservoir properties, including lithologic and depositional characteristics of the reservoir formations and porosity and fluid flow properties through different compartments in the reservoir. The team has also developed models that represent the geological and physical characteristics of the reservoir and seal system.

01 COAL, LIGNITE, AND PEAT↗

DOE OSTI2020

Reservoirs Modify River Thermal Regime Sensitivity to Climate Change: A Case Study in the Southeastern United States

Seasonal thermal stratification in reservoirs changes the thermal regime of regulated river systems as well as stream temperature responses to climate change. Cold releases from reservoir hypolimnion can depress downstream river temperature during warm seasons. Recent large-scale climate change studies on stream temperature have largely ignored reservoir thermal stratification. In this study, we used established models to develop a framework which considers water demand and reservoir regulation with thermal stratification, and applied this model framework to the southeastern US. About half of all 271 reservoirs in our study area retain strong thermal stratification by the 2080s (2070-2099) under RCP8.5 even as median residence times decrease to 60 days from 69 days in the historic period (1979-2010). Reservoir impacts on downstream temperatures become slightly weaker in the future because of higher air temperature and stronger solar radiation. We defined a “cooling potential” to quantify the thermal energy that a water body can absorb before exceeding a water temperature threshold. In the future, higher river temperatures will reduce the cooling potential for all river segments, but more so for river segments minimally impacted by thermal stratification. Reservoir impacts on cooling potential remain strong for river segments downstream of reservoirs with strong thermal stratification. We conducted a sensitivity analysis to evaluate the robustness of our findings to errors in the hydrological simulations. While river segments subject to reservoir regulation are more sensitive to errors in hydrology than those without regulation impacts, our overall findings do not materially change due to these errors.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2024

Predicting Critical Transitions in Multiscale Data

Predicting the dynamics of complex nonlinear systems remains a challenging problem both in dynamical systems theory as well as real world science and engineering applications. Data-driven methods utilizing the latest advances in machine learning (ML) provide a promising new paradigm for this task. Our work centered on Reservoir Computing (RC), which has shown itself to be capable of skillfully predicting chaotic dynamics in multiscale systems. In the first part of the work, the focus is on how to improve predictions of critical transitions in a class of slow-fast metastable systems in which the equations are known. An additional goal was to determine whether a relationship exists between RC and Koopman operator theory, to improve the efficiency and broaden the applicability of the approach. In the second part of this work, a variation on the RC model known as Reconstructive Reservoir Computing (RRC) is applied to real-world data to identify anomalies.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2022

Supplementary Data for "Evaluation of the Economic Implications of Varied Pressure Drawdown Strategies Generated Using a Real-time, Rapid Predictive, Multi-fidelity Model for Unconventional Oil and Gas Wells" by Bello, K., Vikara, D., Sheriff, A., Viswanathan, H., Carr, T., Sweeney, M., O'Malley, D., Marquis, M., Vactor, R.T., and Cunha, L.

The Bello et al. study evaluates the impact of contrasting pressure drawdown on gas productivity and the resulting economics of a well in the Marcellus Shale of the Appalachian Basin. This research applies a techno-economic analysis approach to help identify potential ways pressure management strategies can be used to improve cumulative recovery of hydraulically fractured horizontal wells while maintaining project profitability. Gas production forecast outlook scenarios of the Marcellus Shale Energy and Environment Laboratory Laboratory's MIP-3H well were generated under varying pressure drawdowns using two approaches: 1) a novel physics-informed machine learning (PIML) workflow and 2) via traditional reservoir simulation in Computer Modeling Group’s (CMG) GEM Compositional & Unconventional Simulator. Cash flow and other economic metrics of interest were compiled on the production outlook using the U.S. Department of Energy's (DOE) National Energy Technology Laboratory (NETL) Unconventional Shale Well Economic Model (UShWEM).The sheets within this Microsoft ExcelTM workbook provide the economic metric outputs for the baseline condition and the one-at-a-time (OAT) sensitivity analysis of UShWEM's input parameters for each of the production scenarios evaluated.

Fracture Network Model↗

DOE OSTI2020

Hydraulic fracturing to enhance injectivity and storage capacity of CO 2 storage reservoirs: Benefits and risks

Several potential CO 2 storage reservoirs have been found to have insufficient porosity and permeability to support cost effective commercial-scale injection. As a result, the use of hydraulic fracturing to enhance injectivity and storage capacity of CO 2 storage reservoirs was explored. Previous modeling studies indicate that fracturing can increase storage capacity by modest to significant amounts (10%–35%), depending on model assumptions. Simulations completed as part of this study confirm that for a range of horizontal well lengths, number of fractures, fracture geometries, and fracture properties, injectivity is improved and capacity increases by 13%–71% over the base case unfractured vertical well. Intuitively, increasing the well length and number of fractures had a corresponding impact on increased capacity. Here, fracture area (i.e., the fracture height multiplied by its width) was an important parameter for increasing capacity but the specific geometry (e.g., the ratio of height-to-width) was unimportant. The most important aspect that affected capacity was the ability of a fracture to connect high permeability horizontal zones in the reservoir. This would allow a single well to access both layers, thereby maximizing total storage capacity of the reservoir while likely leading to an overall increase of the CO 2 footprint, which is an important consideration for Class VI UIC permitting. The results of this work demonstrate that hydraulic fracturing is an attractive option to consider when faced with an underperforming geologic carbon storage site that is at risk of causing a project to fail.

03 NATURAL GAS↗

DOE OSTI2021

NRAP-Open-IAM: Open Wellbore Component (V.2.0)

The Open Wellbore component of NRAP-Open-IAM is applicable to the calculation of Area of Review at a carbon storage site and fast estimation of dense gas dispersion from multiple continuous CO 2 surface leakage sources for risk assessment. The Open Wellbore Component of NRAP-Open-IAM is a lookup table model based on the drift-flux approach. The National Research Assessment Partnership’s second-generation integrated assessment model, NRAP-Open-IAM, is open-source software written in Python for use in performance and quantitative risk assessment of geologic carbon storage (GCS) systems. The look-up table for open well leakage is an updated version the previous model which provides finer resolution of reservoir depth and more accurate fluid properties than the previous open well look-up table, especially in calculating the gas (CO 2 -rich) phase properties. The component model input parameters include the reservoir transmissivity, aquifer transmissivity, brine salinity, well radius, and well top and bottom depths. The reservoir transmissivity and brine salinity were set to values consistent with the reservoir component, and the well top and bottom depths varied with location. The possible outputs from the Open Wellbore component are leakage rates of CO 2 and brine to either an aquifer or atmosphere, depending on the depth of the well top. Several python test scripts are presented that demonstrate that the Open Wellbore component works correctly. Use of the open wellbore model in NRAP-Open-IAM may result in a large CO 2 leakage rates, comparable to the leakage rates of CO 2 blow out.

54 ENVIRONMENTAL SCIENCES↗