James E. Bean

Computational thermal, chemical, fluid, and solid mechanics for geosystems management.

This document summarizes research performed under the SNL LDRD entitled - Computational Mechanics for Geosystems Management to Support the Energy and Natural Resources Mission. The main accomplishment was development of a foundational SNL capability for computational thermal, chemical, fluid, and solid mechanics analysis of geosystems. The code was developed within the SNL Sierra software system. This report summarizes the capabilities of the simulation code and the supporting research and development conducted under this LDRD. The main goal of this project was the development of a foundational capability for coupled thermal, hydrological, mechanical, chemical (THMC) simulation of heterogeneous geosystems utilizing massively parallel processing. To solve these complex issues, this project integrated research in numerical mathematics and algorithms for chemically reactive multiphase systems with computer science research in adaptive coupled solution control and framework architecture. This report summarizes and demonstrates the capabilities that were developed together with the supporting research underlying the models. Key accomplishments are: (1) General capability for modeling nonisothermal, multiphase, multicomponent flow in heterogeneous porous geologic materials; (2) General capability to model multiphase reactive transport of species in heterogeneous porous media; (3) Constitutive models for describing real, general geomaterials under multiphase conditions utilizing laboratory data; (4) General capability to couple nonisothermal reactive flow with geomechanics (THMC); (5) Phase behavior thermodynamics for the CO2-H2O-NaCl system. General implementation enables modeling of other fluid mixtures. Adaptive look-up tables enable thermodynamic capability to other simulators; (6) Capability for statistical modeling of heterogeneity in geologic materials; and (7) Simulator utilizes unstructured grids on parallel processing computers.

Geomechanical analyses to investigate wellbore/mine interactions in the Potash Enclave of Southeastern New Mexico.

Geomechanical analyses have been performed to investigate potential mine interactions with wellbores that could occur in the Potash Enclave of Southeastern New Mexico. Two basic models were used in the study; (1) a global model that simulates the mechanics associated with mining and subsidence and (2) a wellbore model that examines the resulting interaction impacts on the wellbore casing. The first model is a 2D approximation of a potash mine using a plane strain idealization for mine depths of 304.8 m (1000 ft) and 609.6 m (2000 ft). A 3D wellbore model then considers the impact of bedding plane slippage across single and double cased wells cemented through the Salado formation. The wellbore model establishes allowable slippage to prevent casing yield.