Search NASA⌕ Search

Engineering topics

Celia, Michael A.

Publications and source records attributed to Celia, Michael A..

13 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

A guideline for appropriate application of vertically-integrated modeling approaches for geologic carbon storage modeling

Mathematical modeling is an essential tool for answering questions related to geologic carbon storage (GCS). The choice of modeling approach depends on the type of questions being asked. In this paper we discuss a series of approaches with a hierarchical complexity including vertically-integrated single-phase flow approaches, vertically-integrated multi-phase flow approaches (with and without vertical equilibrium assumption), three-dimensional multi-phase flow approaches, and fully-coupled multi-phase flow approaches that couple flow with geochemistry and/or geomechanics. Three spatial scales are used to categorize the questions to be addressed by modeling: regional scale (encompasses CO 2 plume extent and majority of area of pressure impact of one or more injection operations), site scale (includes the CO 2 plume extent and some of the area impacted by the pressure increase of a single injection site), and well scale (the immediate vicinity of an injection well). A set of guidelines is developed to help modelers choose the most appropriate modeling approach, and show when simpler modeling approaches may be the better choice. Vertically-integrated single-phase flow models are the most appropriate choice at both the site and regional scales, if the pressure impact outside of the CO 2 plume is of interest. Vertically-integrated multi-phase flow models should be chosen at the regional scale, if the locations of CO 2 plumes are of interest, and at the site scale if vertical segregation of CO 2 and brine is fast or vertical heterogeneity in properties can be presented by distinct, continuous layers. Three-dimensional multi-phase flow models are the appropriate choice at the well and site scales for cases with significant vertical flow components of CO 2 and brine. Fully-coupled multi-phase flow models should only be chosen if pore-space alteration through geochemistry or geomechanics feeds back to fluid flow.

42 ENGINEERING↗
DOE OSTI2019

Multiscale Modeling of CO 2 Migration and Trapping in Fractured Reservoirs with Validation by Model Comparison and Real-Site Applications (Final Technical Report)

This report documents the accomplishments achieved during the project titled “Multiscale Modeling of CO 2 Migration and Trapping in Fractured Reservoirs with Validation by Model Comparison and Real-Site Applications” funded by the US Department of Energy, Office of Fossil Energy. The objectives of the project were to develop modeling capabilities for predicting CO 2 and brine migration in fractured reservoirs during geologic carbon storage and to investigate the feasibility of carbon storage in fractured reservoirs using the newly developed modeling capabilities. To achieve these objectives new mass transfer functions were developed to adapt existing dual-continuum approaches – a commonly used approach in hydrocarbon reservoir modeling – to the CO 2 -brine system. In the dual-continuum approach the undisturbed rock matrix and the fractures are modeled as two separate continua, which are coupled through mass exchange of CO 2 and brine between the two continua. Mass transfer functions represent the mass exchange between fractures and the rock matrix, without requiring highly resolved numerical grids of individual fractures, thus allowing large modeling domains relevant to geologic carbon storage questions (tens of meters to hundreds of kilometers). Mass transfer functions for the exchange of both free-phase CO 2 and brine between fractures and rock matrix were developed for both the injection phase and the post-injection phase. During the injection phase the fractures are filled with CO 2 , while the rock matrix is initially filled with brine. Therefore, a mass transfer function was developed that is based on vertical displacement of brine by CO 2 in the rock matrix due to gravitational forces. Once CO 2 injection ceases and CO 2 migrates away from the direct vicinity of the injection site, the fractures will begin to re-fill with brine, while the rock matrix will have high CO 2 saturations. For these conditions a mass transfer functions was developed that is based on spontaneous imbibition, as the pore spaces in the rock matrix are much smaller than in the fractures. For formations with low permeability or high capillary entry pressure of the rock matrix, free-phase CO 2 may not be able to enter the rock matrix. However, free-phase CO 2 at the interface between fractures and rock matrix can dissolve into brine resident in the rock matrix. Therefore, additional mass transfer functions were developed based on aqueous-phase diffusion of CO 2 through the rock matrix. These mass transfer functions are based on existing solutions for diffusion and were expanded to be applicable to a variety of block shapes relevant to dual-continuum modeling. Implementations of the newly developed mass transfer functions showed that they give accurate results when compared to highly resolved models. In addition to developing new transfer functions, a new vertically-integrated dual-continuum modeling framework was developed. This development takes advantage of fast vertical segregation of CO 2 and brine in the fractures due to the density difference between CO 2 and brine and the high permeability of the fractures. Therefore, a vertically-integrated vertical equilibrium model for the fractures is coupled with a vertically-distributed model for the rock matrix. The newly developed approach significantly reduces the computational effort, but is only accurate for fracture permeabilities of 100 mD or higher. Finally, the newly developed models were applied to real and hypothetical geologic carbon storage sites to investigate the feasibility of using fractured reservoirs for carbon storage. The general results were, that storage of CO 2 – both as free-phase and dissolved – is feasible, but that storage capacities are likely lower than in unfractured systems. While the high permeability of the fractures reduces the injection pressure (or permits higher injection rates), fast migration of CO 2 may limit the overall storage volume, because CO 2 may escape laterally with not enough of the rock matrix being swept for effective storage.

20 FOSSIL-FUELED POWER PLANTS↗
DOE OSTI2016

Model complexity and choice of model approaches for practical simulations of CO 2 injection, migration, leakage and long-term fate

This report documents the accomplishments achieved during the project titled “Model complexity and choice of model approaches for practical simulations of CO 2 injection,migration, leakage and long-term fate” funded by the US Department of Energy, Office of Fossil Energy. The objective of the project was to investigate modeling approaches of various levels of complexity relevant to geologic carbon storage (GCS) modeling with the goal to establish guidelines on choice of modeling approach.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Axisymmetric flows from fluid injection into a confined porous medium

In this work, we study the axisymmetric flows generated from fluid injection into a horizontal confined porous medium that is originally saturated with another fluid of different density and viscosity. Neglecting the effects of surface tension and fluid mixing, we use the lubrication approximation to obtain a nonlinear advection-diffusion equation that describes the time evolution of the sharp fluid-fluid interface. The flow behaviors are controlled by two dimensionless groups: M , the viscosity ratio of displaced fluid relative to injected fluid, and Γ , which measures the relative importance of buoyancy and fluid injection. For this axisymmetric geometry, the similarity solution involving R 2 /T (where R is the dimensionless radial coordinate and T is the dimensionless time) is an exact solution to the nonlinear governing equation for all times. Four analytical expressions are identified as asymptotic approximations (two of which are new solutions): (i) injection-driven flow with the injected fluid being more viscous than the displaced fluid ( Γ $\ll$ 1 and M < 1) where we identify a self-similar solution that indicates a parabolic interface shape; (ii) injection-driven flow with injected and displaced fluids of equal viscosity ( Γ $\ll$ 1 and M = 1), where we find a self-similar solution that predicts a distinct parabolic interface shape; (iii) injection-driven flow with a less viscous injected fluid ( Γ $\ll$ 1 and M > 1) for which there is a rarefaction wave solution, assuming that the Saffman-Taylor instability does not occur at the reservoir scale; and (iv) buoyancy-driven flow ( Γ $\gg$ 1) for which there is a well-known self-similar solution corresponding to gravity currents in an unconfined porous medium [S. Lyle et al. “Axisymmetric gravity currents in a porous medium,” J. Fluid Mech. 543, 293–302 (2005)]. The various axisymmetric flows are summarized in a Γ - M regime diagram with five distinct dynamic behaviors including the four asymptotic regimes and an intermediate regime. The implications of the regime diagram are discussed using practical engineering projects of geological CO 2 sequestration, enhanced oil recovery, and underground waste disposal.

42 ENGINEERING↗
DOE OSTI2015

Flow regimes for fluid injection into a confined porous medium

We report theoretical and numerical studies of the flow behaviour when a fluid is injected into a confined porous medium saturated with another fluid of different density and viscosity. For a two-dimensional configuration with point source injection, a nonlinear convection–diffusion equation is derived to describe the time evolution of the fluid–fluid interface. In the early time period, the fluid motion is mainly driven by the buoyancy force and the governing equation is reduced to a nonlinear diffusion equation with a well-known self-similar solution. In the late time period, the fluid flow is mainly driven by the injection, and the governing equation is approximated by a nonlinear hyperbolic equation that determines the global spreading rate; a shock solution is obtained when the injected fluid is more viscous than the displaced fluid, whereas a rarefaction wave solution is found when the injected fluid is less viscous. In the late time period, we also obtain analytical solutions including the diffusive term associated with the buoyancy effects (for an injected fluid with a viscosity higher than or equal to that of the displaced fluid), which provide the structure of the moving front. Numerical simulations of the convection–diffusion equation are performed; the various analytical solutions are verified as appropriate asymptotic limits, and the transition processes between the individual limits are demonstrated.

58 GEOSCIENCES↗
DOE OSTI2015

Multiphase Modeling of Geologic Carbon Sequestration in Saline Aquifers

Abstract Geologic carbon sequestration ( GCS ) is being considered as a climate change mitigation option in many future energy scenarios. Mathematical modeling is routinely used to predict subsurface CO 2 and resident brine migration for the design of injection operations, to demonstrate the permanence of CO 2 storage, and to show that other subsurface resources will not be degraded. Many processes impact the migration of CO 2 and brine, including multiphase flow dynamics, geochemistry, and geomechanics, along with the spatial distribution of parameters such as porosity and permeability. In this article, we review a set of multiphase modeling approaches with different levels of conceptual complexity that have been used to model GCS . Model complexity ranges from coupled multiprocess models to simplified vertical equilibrium (VE) models and macroscopic invasion percolation models. The goal of this article is to give a framework of conceptual model complexity, and to show the types of modeling approaches that have been used to address specific GCS questions. Application of the modeling approaches is shown using five ongoing or proposed CO 2 injection sites. For the selected sites, the majority of GCS models follow a simplified multiphase approach, especially for questions related to injection and local‐scale heterogeneity. Coupled multiprocess models are only applied in one case where geomechanics have a strong impact on the flow. Owing to their computational efficiency, VE models tend to be applied at large scales. A macroscopic invasion percolation approach was used to predict the CO 2 migration at one site to examine details of CO 2 migration under the caprock.

Bandilla, Karl W.↗
DOE OSTI2013

Effects of Pore Structure Change and Multi-Scale Heterogeneity on Contaminant Transport and Reaction Rate Upscaling

This project addressed the scaling of geochemical reactions to core and field scales, and the interrelationship between reaction rates and flow in porous media. We targeted reactive transport problems relevant to the Hanford site – specifically the reaction of highly caustic, radioactive waste solutions with subsurface sediments, and the immobilization of 90 Sr and 129 I through mineral incorporation and passive flow blockage, respectively. We addressed the correlation of results for pore-scale fluid-soil interaction with field-scale fluid flow, with the specific goals of (i) predicting attenuation of radionuclide concentration; (ii) estimating changes in flow rates through changes of soil permeabilities; and (iii) estimating effective reaction rates. In supplemental work, we also simulated reactive transport systems relevant to geologic carbon sequestration. As a whole, this research generated a better understanding of reactive transport in porous media, and resulted in more accurate methods for reaction rate upscaling and improved prediction of permeability evolution. These scientific advancements will ultimately lead to better tools for management and remediation of DOE legacy waste problems.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Effects of Pore Structure Change and Multi-Scale Heterogeneity on Contaminant Transport and Reaction Rate Upscaling

This project addressed the scaling of geochemical reactions to core and field scales, and the interrelationship between reaction rates and flow in porous media. We targeted reactive transport problems relevant to the Hanford site – specifically the reaction of highly caustic, radioactive waste solutions with subsurface sediments, and the immobilization of 90 Sr and 129 I through mineral incorporation and passive flow blockage, respectively. We addressed the correlation of results for pore-scale fluid-soil interaction with field-scale fluid flow, with the specific goals of (i) predicting attenuation of radionuclide concentration; (ii) estimating changes in flow rates through changes of soil permeabilities; and (iii) estimating effective reaction rates. In supplemental work, we also simulated reactive transport systems relevant to geologic carbon sequestration. As a whole, this research generated a better understanding of reactive transport in porous media, and resulted in more accurate methods for reaction rate upscaling and improved prediction of permeability evolution. These scientific advancements will ultimately lead to better tools for management and remediation of DOE’s legacy waste problems.

58 GEOSCIENCES↗
DOE OSTI2011

EMSL Pore Scale Modeling Challenge/Workshop

The design and implementation of cost-effective remediation strategies for subsurface contamination not only relies on descriptive models of contaminant distribution but predictive models of contaminant fate and transport. The development of such predictive models is of great importance to many government agencies including the Department of Energy (DOE), Department of Defense, and the Environmental Protection Agency. Accurate prediction of the fate and transport of environmental contaminants in the subsurface has driven development of increasingly sophisticated subsurface flow and transport models. Even without considering spatial and temporal scalability, contaminant transport at the field scale is complex, and this almost inevitably results in an especially sophisticated model being developed for a specific field site. Unfortunately, these models are not easily generalizable to other field sites or and can only be used by the model authors to generate meaningful results. Over time, this results in the existence of numerous specialized models that are not readily comparable to each other and difficult to cross-validate. There is a lack of fundamental experimental observations at the pore scale that can be used to support model development and validation at the field scale.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗