Engineering topics

Oldenburg, Curtis M.

Publications and source records attributed to Oldenburg, Curtis M..

8 records

DOE OSTI2022

Computational Tools and Workflows for Quantitative Risk Assessment and Decision Support for Geologic Carbon Storage Sites: Progress and Insights from the U.S. DOE’s National Risk Assessment Partnership

The 2005 Intergovernmental Panel on Climate Change (IPCC) Special Report on CCS raised the profile of CO2 capture and storage (CCS) as an important technology for reducing greenhouse gas (GHG) emissions. CCS is now recognized as a key component of most climate change mitigation scenarios. Since publication of that report the international research, development, and deployment (RD&D) community has advanced key technical aspects, clarified regulatory requirements, explored value chain and infrastructure solutions, and developed incentive paradigms to enable and promote large-scale deployment of CCS. These efforts have included research to better characterize geologic storage resources, to improve injection performance and storage efficiency, to assess and manage subsurface environmental risks, and to advance monitoring technologies to assure system conformance. These efforts have helped to build confidence in the viability of geologic carbon storage (GCS), but stakeholder concerns about long-term risks and liability associated with GCS remain a hurdle to broad acceptance and large-scale deployment of CCS. Since 2010, the U.S. DOE’s National Risk Assessment Partnership (NRAP) – a research collaboration between five contributing national laboratories – has worked to establish and demonstrate methods and tools to quantify and manage the subsurface environmental risks associated with GCS, amidst uncertainty. This work supports the Office of Fossil Energy and Carbon Management Carbon Transport and Storage Program’s goal of advancing safe and secure commercial-scale GCS deployment. To address the technical challenge of simulating the physical response of the GCS site to large-scale CO2 injection, NRAP has adopted an approach that relies on coupling computationally efficient reduced-order and/or data-driven proxy models of important system components (i.e., storage reservoir, sealing caprock, leakage pathways, intermediate formations, overlying groundwater aquifers, and the atmosphere) in integrated assessment framework. That integrated model of the physical system is complemented with fit-for purpose functionality to support site characterization and risk-related decisions. The recently released NRAP Phase II toolset includes the Open-Source Integrated Assessment Model (NRAP-Open-IAM) for evaluation of trends in leakage risk and potential impact, tools to support monitoring design optimization (Designs for Risk Evaluation and Management – DREAM v3.0 and Passive Seismic Monitoring Tool - PSMT), and tools for state of stress evaluation (State-of-Stress Analysis Tool - SOSAT) and forecasting induced seismicity risk. The NRAP team has also released a pair of reports describing conceptual workflows to incorporate physics-based, quantitative risk assessment into many of the design, planning, operation, and closure decisions for GCS projects. An online catalogue highlights published studies where these tools and methods are demonstrated. In this presentation, the utility of these products to assess risks and address key stakeholder questions will be highlighted through examples, and related insights about the safety and security of geologic carbon storage in qualified storage sites will be discussed. The prospect of rapid, large-scale deployment of GCS technology to aggressively reduce anthropogenic CO2 emissions requires careful consideration of interference between multiple commercial-scale storage projects within a basin. Going forward, NRAP is expanding and adapting site-scale risk quantification tools and methods to enable assessment of risks and inform management decisions for basin-scale deployment. Increasingly, this work will leverage next-generation approaches for surrogate modelling, fast prediction, and advanced visualization enabled by machine learning and artificial intelligence to promote virtual learning, scenario evaluation, and augment risk-based decision making.

quantitative risk assessment, geologic carbon stor↗

DOE OSTI2022

Dynamic risk assessment for geologic CO 2 sequestration

At a geologic CO 2 sequestration (GCS) site, geologic uncertainty usually leads to large uncertainty in the predictions of properties that influence metrics for leakage risk assessment, such as CO 2 saturations and pressures in potentially leaky wellbores, CO 2 /brine leakage rates, and leakage consequences such as changes in drinking water quality in groundwater aquifers. The large uncertainty in these risk-related system properties and risk metrics can lead to over-conservative risk management decisions to ensure safe operations of GCS sites. The objective of this work is to develop a novel approach based on dynamic risk assessment to effectively reduce the uncertainty in the predicted risk-related system properties and risk metrics. We demonstrate our framework for dynamic risk assessment on two case studies: a 3D synthetic example and a synthetic field example based on the Rock Springs Uplift (RSU) storage site in Wyoming, USA. Results show that the U.S. National Risk Assessment Partnership’s Open Source Integrated Assessment Model (NRAP-Open-IAM) coupled with a conformance evaluation can be used to effectively quantify and reduce the uncertainty in the predictions of risk-related system properties and risk metrics in GCS.

58 GEOSCIENCES↗

DOE OSTI2022

Coupled Hydromechanical Modeling of Induced Seismicity From CO 2 Injection in the Illinois Basin

Injection of CO 2 for geologic carbon sequestration into deep sedimentary formations involves fluid pressure increases that engage hydromechanical processes that can cause seismicity by activation of existing faults. In this work, we use a coupled multiphase fluid flow and geomechanical simulator to model spatiotemporal fluid pressure and stress changes in order to study the poroelastic effect of CO 2 injection on faults in crystalline basement rock below the injection zone. The seismicity rate along features interpreted to be basement faults is modeled using Dieterich's rate-and-state earthquake nucleation model. The methodology is applied to microseismicity detected during CO 2 injection into the Mount Simon formation during the Illinois Basin—Decatur Project. The modeling accurately captures an observed reduction in seismicity rate when the injection in the second well was into a slightly shallower zone above the base of the Mount Simon formation. Moreover, the modeling shows that it is important to consider poroelastic stress changes, in addition to fluid pressure changes for accurately modeling of the observed seismicity rate.

58 GEOSCIENCES↗

DOE OSTI2022

Downwind dispersion of CO 2 from a major subsea blowout in shallow offshore waters

Growing interest in offshore geologic carbon sequestration (GCS) motivates risk assessment of large-scale subsea CO 2 well blowouts or pipeline ruptures. For major leaks of CO 2 from wells or pipelines, significant fluxes of CO 2 may occur from the sea surface depending on water depth. In the context of risk assessment of human health and safety, we have used previously simulated coupled well-reservoir and water column model results as a source term for dense gas dispersion of CO 2 above the sea surface. The models are linked together by one-way coupling, that is, output of one model is used as input to the next model. These first-of-their-kind coupled flow results are applicable to assessing the hazard of CO 2 to people at and downwind of the sea surface location of emission. Hazard is quantified by plotting the downwind dispersion length (DDL), which we define in the study as the distances from the emission source to the point at which the emitted CO 2 has been diluted to 5% and 1.5% in air by volume. Here results suggest that large-scale blowouts in shallow water (10 m) may cause hazardous CO 2 plumes extending on the order of several hundred meters downwind. Details of the modeling show DDL has a maximum for windspeed (at an elevation of 10 m) of approximately 5 m/s, with smaller DDL for both weaker and stronger winds. This is explained by the fact that wind favors transport but also causes dispersion; therefore there is a certain wind speed that maximizes DDL.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1992

Hybrid model for solidification and convection

A hybrid model for continuum phase-change problems is presented. The hybrid model accounts for flow in regions of concentrated mush, dilute mush, and single-phase liquid. Scale analysis shows that, in dilute mush, the Blake-Kozeny-Carman relation may lead to inaccuracy in the continuum formulation for certain values of the Rayleigh and Darcy numbers. The hybrid model uses arctangent switching functions to switch on the Darcy flow and variable viscosity terms depending on the local value of the fraction solid. Two-dimensional example calculations suggest that the hybrid model more accurately accounts for transport in the dilute mush region.

Oldenburg, Curtis M.↗

NASA NTRS1991

Numerical modeling of solidification and convection in a viscous pure binary eutectic system

The solidification and convection of the pure binary eutectic silicate system diopside-anorthite (Di-An) is numerically modeled. A mass-weighted enthalpy of fusion is used to account for the second solid phase (An) which crystallizes at the solidus temperature. Variable under-relaxation is used to aid convergence of the momentum equations in the implementation of the SIMPLER algorithm used to solve the two-dimensional continuum conservation equations. Numerical experiments of the solidification of Di80 melt show that a large temperature drop occurs across the solid and mush regions which decreases convective vigor in the liquid. Interesting compositional segregation patterns are produced during the solidification of Di80.

Oldenburg, Curtis M.↗

NASA NTRS1989

Magma zonation - Effects of chemical buoyancy and diffusion

Numerical simulations and scale analysis are used to assess the viability of the marginal box-filling mechanism for producing compositional zonation in magma bodies. Scale analysis and two-dimensional numerical experiments both show that box-filling occurs provided a critical ratio of compositional-to-thermal buoyancy is exceeded. This critical ratio depends on the ratio of thermal-to-chemical diffusivity; application of this result to magma bodies suggests that box-filling may occur for components with relatively high-chemical diffusivities such as water. However, box-filling will not produce significant zonation for components with small chemical diffusivities, such as silica, unless diffusive coupling increases silica diffusivity.

Spera, Frank J.↗

NASA NTRS1989

Dynamic mixing in magma bodies - Theory, simulations, and implications

The magma-mixing process is different from the mantle mixing process in that the mixing components of magma are dynamically active, with the melt density depending strongly on composition. This paper describes simulations of time-dependent variable-viscosity double-diffusive convection which were carried out to investigate quantitatively the mixing dynamics of magma in melt-dominated magma bodies. Results show that the dynamics of double-diffusive convection can impart complex patterns of composition, through time and space. The mixing time depends nonlinearly on many factors, including heat flux driving convection, the rate of diffusion of chemical species, the relative importance of thermal and chemical buoyancy, the viscosities of the mixing components, and the shape of the magma body.

Oldenburg, Curtis M.↗