skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Three-dimensional analysis of a faulted CO 2 reservoir using an Eshelby-Mori-Tanaka approach to rock elastic properties and fault permeability

Abstract

This work develops a three-dimensional multiscale model to analyze a complex CO 2 faulted reservoir that includes some key geological features of the San Andreas and nearby faults southwest of the Kimberlina site. The model uses the STOMP-CO 2 code for flow modeling that is coupled to the ABAQUS® finite element package for geomechanical analysis. A 3D ABAQUS® finite element model is developed that contains a large number of 3D solid elements with two nearly parallel faults whose damage zones and cores are discretized using the same continuum elements. Five zones with different mineral compositions are considered: shale, sandstone, fault damaged sandstone, fault damaged shale, and fault core. Rocks’ elastic properties that govern their poroelastic behavior are modeled by an Eshelby-Mori-Tanka approach (EMTA). EMTA can account for up to 15 mineral phases. The permeability of fault damage zones affected by crack density and orientations is also predicted by an EMTA formulation. A STOMP-CO 2 grid that exactly maps the ABAQUS® finite element model is built for coupled hydro-mechanical analyses. Simulations of the reservoir assuming three different crack pattern situations (including crack volume fraction and orientation) for the fault damage zones are performed to predict the potential leakage of CO 2more » due to cracks that enhance the permeability of the fault damage zones. Here, the results illustrate the important effect of the crack orientation on fault permeability that can lead to substantial leakage along the fault attained by the expansion of the CO 2 plume. Potential hydraulic fracture and the tendency for the faults to slip are also examined and discussed in terms of stress distributions and geomechanical properties.« less

Authors:
 [1];  [1];  [1];  [1]
  1. Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
Publication Date:
Research Org.:
Pacific Northwest National Lab. (PNNL), Richland, WA (United States)
Sponsoring Org.:
USDOE Office of Fossil Energy (FE)
OSTI Identifier:
1329338
Report Number(s):
PNNL-SA-117116
Journal ID: ISSN 1674-7755; AA7020000
Grant/Contract Number:
AC05-76RL01830
Resource Type:
Journal Article: Published Article
Journal Name:
Journal of Rock Mechanics and Geotechnical Engineering
Additional Journal Information:
Journal Volume: 8; Journal Issue: 6; Journal ID: ISSN 1674-7755
Publisher:
Chinese Society for Rock Mechanics and Engineering - Elsevier
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; 58 GEOSCIENCES; CO2 reservoir; geomechanical modeling; mineralogy; homogenization; fault; elastic properties; leakage; slip; permeability

Citation Formats

Nguyen, Ba Nghiep, Hou, Zhangshuan, Last, George V., and Bacon, Diana H. Three-dimensional analysis of a faulted CO2 reservoir using an Eshelby-Mori-Tanaka approach to rock elastic properties and fault permeability. United States: N. p., 2016. Web. doi:10.1016/j.jrmge.2016.06.007.
Nguyen, Ba Nghiep, Hou, Zhangshuan, Last, George V., & Bacon, Diana H. Three-dimensional analysis of a faulted CO2 reservoir using an Eshelby-Mori-Tanaka approach to rock elastic properties and fault permeability. United States. doi:10.1016/j.jrmge.2016.06.007.
Nguyen, Ba Nghiep, Hou, Zhangshuan, Last, George V., and Bacon, Diana H. 2016. "Three-dimensional analysis of a faulted CO2 reservoir using an Eshelby-Mori-Tanaka approach to rock elastic properties and fault permeability". United States. doi:10.1016/j.jrmge.2016.06.007.
@article{osti_1329338,
title = {Three-dimensional analysis of a faulted CO2 reservoir using an Eshelby-Mori-Tanaka approach to rock elastic properties and fault permeability},
author = {Nguyen, Ba Nghiep and Hou, Zhangshuan and Last, George V. and Bacon, Diana H.},
abstractNote = {This work develops a three-dimensional multiscale model to analyze a complex CO2 faulted reservoir that includes some key geological features of the San Andreas and nearby faults southwest of the Kimberlina site. The model uses the STOMP-CO2 code for flow modeling that is coupled to the ABAQUS® finite element package for geomechanical analysis. A 3D ABAQUS® finite element model is developed that contains a large number of 3D solid elements with two nearly parallel faults whose damage zones and cores are discretized using the same continuum elements. Five zones with different mineral compositions are considered: shale, sandstone, fault damaged sandstone, fault damaged shale, and fault core. Rocks’ elastic properties that govern their poroelastic behavior are modeled by an Eshelby-Mori-Tanka approach (EMTA). EMTA can account for up to 15 mineral phases. The permeability of fault damage zones affected by crack density and orientations is also predicted by an EMTA formulation. A STOMP-CO2 grid that exactly maps the ABAQUS® finite element model is built for coupled hydro-mechanical analyses. Simulations of the reservoir assuming three different crack pattern situations (including crack volume fraction and orientation) for the fault damage zones are performed to predict the potential leakage of CO2 due to cracks that enhance the permeability of the fault damage zones. Here, the results illustrate the important effect of the crack orientation on fault permeability that can lead to substantial leakage along the fault attained by the expansion of the CO2 plume. Potential hydraulic fracture and the tendency for the faults to slip are also examined and discussed in terms of stress distributions and geomechanical properties.},
doi = {10.1016/j.jrmge.2016.06.007},
journal = {Journal of Rock Mechanics and Geotechnical Engineering},
number = 6,
volume = 8,
place = {United States},
year = 2016,
month = 9
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record at 10.1016/j.jrmge.2016.06.007

Save / Share:
  • This work develops a three-dimensional multiscale model to analyze a complex CO 2 faulted reservoir that includes some key geological features of the San Andreas and nearby faults southwest of the Kimberlina site. The model uses the STOMP-CO 2 code for flow modeling that is coupled to the ABAQUS® finite element package for geomechanical analysis. A 3D ABAQUS® finite element model is developed that contains a large number of 3D solid elements with two nearly parallel faults whose damage zones and cores are discretized using the same continuum elements. Five zones with different mineral compositions are considered: shale, sandstone, faultmore » damaged sandstone, fault damaged shale, and fault core. Rocks’ elastic properties that govern their poroelastic behavior are modeled by an Eshelby-Mori-Tanka approach (EMTA). EMTA can account for up to 15 mineral phases. The permeability of fault damage zones affected by crack density and orientations is also predicted by an EMTA formulation. A STOMP-CO 2 grid that exactly maps the ABAQUS® finite element model is built for coupled hydro-mechanical analyses. Simulations of the reservoir assuming three different crack pattern situations (including crack volume fraction and orientation) for the fault damage zones are performed to predict the potential leakage of CO 2 due to cracks that enhance the permeability of the fault damage zones. Here, the results illustrate the important effect of the crack orientation on fault permeability that can lead to substantial leakage along the fault attained by the expansion of the CO 2 plume. Potential hydraulic fracture and the tendency for the faults to slip are also examined and discussed in terms of stress distributions and geomechanical properties.« less
  • This work applies a three-dimensional (3D) multiscale approach recently developed to analyze a complex CO 2 faulted reservoir that includes some key geological features of the San Andreas and nearby faults. The approach couples the STOMP-CO2-R code for flow and reactive transport modeling to the ABAQUS ® finite element package for geomechanical analysis. The objective is to examine the coupled hydro-geochemical-mechanical impact on the risk of hydraulic fracture and fault slip in a complex and representative CO 2 reservoir that contains two nearly parallel faults. STOMP-CO2-R/ABAQUS ® coupled analyses of this reservoir are performed assuming extensional and compressional stress regimesmore » to predict evolutions of fluid pressure, stress and strain distributions as well as potential fault failure and leakage of CO 2 along the fault damage zones. The tendency for the faults to slip and pressure margin to fracture are examined in terms of stress regime, mineral composition, crack distributions in the fault damage zones and geomechanical properties. Here, this model in combination with a detailed description of the faults helps assess the coupled hydro-geochemical-mechanical effect.« less
  • EMTA is a stand-alone computer program that has been developed for the computation of elastic properties and thermal expansion coefficients (thermoelastic properties) of discontinuous fiber composites. EMTA stands for the Eshelby-Mori-Tanaka approach. It implements the standard and modified Mori-Tanaka models that use the Eshelby's equivalent inclusion method. EMTA carries out the Eshelby-Mori-Tanaka homogenization procedure accounting for the constituents (fiber and matrix) properties such as the elastic properties and thermal expansion coefficients (CTEs) of the fibers and of the matrix. It also accounts for the constituents features such as fiber length and orientation distributions, fiber curvature, and imperfect fiber/matrix interfaces. Themore » outputs of an EMTA execution are the elastic properties (engineering constants) and CTEs of the as-formed composite in the defined material coordinate system. These results can readily be used in engineering applications and designs that require these properties.« less
  • Biomass thermochemical conversion, often done in fluidized beds, recently gained a lot of attention due to its potential to efficiently produce renewable liquid fuels. Optimization of reactor design and operating conditions, however, requires a fundamental understanding of bed dynamics. In this work, a numerical framework based on an Euler-Lagrange approach is developed and used to perform and analyze large-scale simulations of two- and three-dimensional periodic fluidized beds. Collisions are handled using a soft-sphere model. An efficient parallel implementation allows one to explicitly track over 30 million particles, which is representative of the number of particles found in lab-scale reactor, thereforemore » demonstrating the capability of Lagrangian approaches to simulate realistic systems at that scale. An on-the-fly bubble identification and tracking algorithm is used to characterize bubble dynamics for inlet velocities up to 9 times the minimum fluidization velocity. Statistics for gas volume fraction, gas and particle velocities, bed expansion, and bubble size and velocity, is compared across the two- and three-dimensional configurations, and comparison with literature data generally shows good agreement. The wide distribution of gas residence times observed in the simulations is linked to the different gas hold-up characteristics of the gas-solid system.« less