Original article: Numerical simulation of two-phase flow in deformable porous media: Application to carbon dioxide storage in the subsurface

  • Authors:
  • O. Kolditz;S. Bauer;N. BöTtcher;D. Elsworth;U. -J. GöRke;C. -I. Mcdermott;C. -H. Park;A. K. Singh;J. Taron;W. Wang

  • Affiliations:
  • Department of Environmental Informatics, Helmholtz Centre for Environmental Research - UFZ Leipzig, Germany and Environmental System Analysis, TU Dresden, Germany;Geohydromodelling, University of Kiel, Germany;Groundwater Management, TU Dresden, Germany;Center for Geomechanics, Geofluids and Geohazards, PennState, USA;Department of Environmental Informatics, Helmholtz Centre for Environmental Research - UFZ Leipzig, Germany;School of Geosciences, Edinburgh University, UK;Department of Environmental Informatics, Helmholtz Centre for Environmental Research - UFZ Leipzig, Germany and Department of Geothermal Resources, Korea Institute of Geoscience and Mineral Resour ...;Department of Environmental Informatics, Helmholtz Centre for Environmental Research - UFZ Leipzig, Germany;Department of Environmental Informatics, Helmholtz Centre for Environmental Research - UFZ Leipzig, Germany and Center for Geomechanics, Geofluids and Geohazards, PennState, USA;Department of Environmental Informatics, Helmholtz Centre for Environmental Research - UFZ Leipzig, Germany

  • Venue:
  • Mathematics and Computers in Simulation
  • Year:
  • 2012

Quantified Score

Hi-index 0.00

Visualization

Abstract

In this paper, conceptual modeling as well as numerical simulation of two-phase flow in deep, deformable geological formations induced by CO"2 injection are presented. The conceptual approach is based on balance equations for mass, momentum and energy completed by appropriate constitutive relations for the fluid phases as well as the solid matrix. Within the context of the primary effects under consideration, the fluid motion will be expressed by the extended Darcy's law for two phase flow. Additionally, constraint conditions for the partial saturations and the pressure fractions of carbon dioxide and brine are defined. To characterize the stress state in the solid matrix, the effective stress principle is applied. Furthermore, the interaction of fluid and solid phases is illustrated by constitutive models for capillary pressure, porosity and permeability as functions of saturation. Based on this conceptual model, a coupled system of nonlinear differential equations for two-phase flow in a deformable porous matrix (H^2M model) is formulated. As the displacement vector acts as primary variable for the solid matrix, multiphase flow is simulated using both pressure/pressure or pressure/saturation formulations. An object-oriented finite element method is used to solve the multi-field problem numerically. The capabilities of the model and the numerical tools to treat complex processes during CO"2 sequestration are demonstrated on three benchmark examples: (1) a 1-D case to investigate the influence of variable fluid properties, (2) 2-D vertical axi-symmetric cross-section to study the interaction between hydraulic and deformation processes, and (3) 3-D to test the stability and computational costs of the H^2M model for real applications.