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

Title: Approximate solutions for diffusive fracture-matrix transfer: Application to storage of dissolved CO 2 in fractured rocks

Abstract

Analytical solutions with infinite exponential series are available to calculate the rate of diffusive transfer between low-permeability blocks and high-permeability zones in the subsurface. Truncation of these series is often employed by neglecting the early-time regime. Here in this paper, we present unified-form approximate solutions in which the early-time and the late-time solutions are continuous at a switchover time. The early-time solutions are based on three-term polynomial functions in terms of square root of dimensionless time, with the first coefficient dependent only on the dimensionless area-to-volume ratio. The last two coefficients are either determined analytically for isotropic blocks (e.g., spheres and slabs) or obtained by fitting the exact solutions, and they solely depend on the aspect ratios for rectangular columns and parallelepipeds. For the late-time solutions, only the leading exponential term is needed for isotropic blocks, while a few additional exponential terms are needed for highly anisotropic rectangular blocks. The optimal switchover time is between 0.157 and 0.229, with highest relative approximation error less than 0.2%. The solutions are used to demonstrate the storage of dissolved CO 2 in fractured reservoirs with low-permeability matrix blocks of single and multiple shapes and sizes. These approximate solutions are building blocks for developmentmore » of analytical and numerical tools for hydraulic, solute, and thermal diffusion processes in low-permeability matrix blocks.« less

Authors:
ORCiD logo [1];  [1];  [2];  [1]
  1. Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States). Energy Geosciences Division
  2. Montana State Univ., Bozeman, MT (United States). Big Sky Carbon Sequestration Partnership
Publication Date:
Research Org.:
Lawrence Berkeley National Lab. (LBNL), Berkeley, CA (United States); Montana State University
Sponsoring Org.:
USDOE Office of Science (SC); USDOE Office of Fossil Energy (FE)
OSTI Identifier:
1393610
Alternate Identifier(s):
OSTI ID: 1430906
Grant/Contract Number:  
AC02-05CH11231; FC26-05NT42587; FE0023323
Resource Type:
Accepted Manuscript
Journal Name:
Water Resources Research
Additional Journal Information:
Journal Volume: 53; Journal Issue: 2; Journal ID: ISSN 0043-1397
Publisher:
American Geophysical Union (AGU)
Country of Publication:
United States
Language:
English
Subject:
54 ENVIRONMENTAL SCIENCES; diffusion; mass transfer; analytical solutions; approximate solutions; fractured rocks; solubility trapping

Citation Formats

Zhou, Quanlin, Oldenburg, Curtis M., Spangler, Lee H., and Birkholzer, Jens T. Approximate solutions for diffusive fracture-matrix transfer: Application to storage of dissolved CO2 in fractured rocks. United States: N. p., 2017. Web. doi:10.1002/2016WR019868.
Zhou, Quanlin, Oldenburg, Curtis M., Spangler, Lee H., & Birkholzer, Jens T. Approximate solutions for diffusive fracture-matrix transfer: Application to storage of dissolved CO2 in fractured rocks. United States. doi:10.1002/2016WR019868.
Zhou, Quanlin, Oldenburg, Curtis M., Spangler, Lee H., and Birkholzer, Jens T. Thu . "Approximate solutions for diffusive fracture-matrix transfer: Application to storage of dissolved CO2 in fractured rocks". United States. doi:10.1002/2016WR019868. https://www.osti.gov/servlets/purl/1393610.
@article{osti_1393610,
title = {Approximate solutions for diffusive fracture-matrix transfer: Application to storage of dissolved CO2 in fractured rocks},
author = {Zhou, Quanlin and Oldenburg, Curtis M. and Spangler, Lee H. and Birkholzer, Jens T.},
abstractNote = {Analytical solutions with infinite exponential series are available to calculate the rate of diffusive transfer between low-permeability blocks and high-permeability zones in the subsurface. Truncation of these series is often employed by neglecting the early-time regime. Here in this paper, we present unified-form approximate solutions in which the early-time and the late-time solutions are continuous at a switchover time. The early-time solutions are based on three-term polynomial functions in terms of square root of dimensionless time, with the first coefficient dependent only on the dimensionless area-to-volume ratio. The last two coefficients are either determined analytically for isotropic blocks (e.g., spheres and slabs) or obtained by fitting the exact solutions, and they solely depend on the aspect ratios for rectangular columns and parallelepipeds. For the late-time solutions, only the leading exponential term is needed for isotropic blocks, while a few additional exponential terms are needed for highly anisotropic rectangular blocks. The optimal switchover time is between 0.157 and 0.229, with highest relative approximation error less than 0.2%. The solutions are used to demonstrate the storage of dissolved CO2 in fractured reservoirs with low-permeability matrix blocks of single and multiple shapes and sizes. These approximate solutions are building blocks for development of analytical and numerical tools for hydraulic, solute, and thermal diffusion processes in low-permeability matrix blocks.},
doi = {10.1002/2016WR019868},
journal = {Water Resources Research},
number = 2,
volume = 53,
place = {United States},
year = {2017},
month = {1}
}

Journal Article:
Free Publicly Available Full Text
Publisher's Version of Record

Citation Metrics:
Cited by: 5 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Mathematical modeling of tracer behavior in short-term experiments in fissured rocks
journal, July 1990


A non-local two-phase flow model for immiscible displacement in highly heterogeneous porous media and its parametrization
journal, December 2013


First-passage-time transfer functions for groundwater tracer tests conducted in radially convergent flow
journal, January 2000


Field-scale effective matrix diffusion coefficient for fractured rock: Results from literature survey
journal, August 2007


A numerical dual-porosity model with semianalytical treatment of fracture/matrix flow
journal, July 1993

  • Zimmerman, Robert W.; Chen, Gang; Hadgu, Teklu
  • Water Resources Research, Vol. 29, Issue 7
  • DOI: 10.1029/93WR00749

Behavior of the mass transfer coefficient during the MADE-2 experiment: New insights: MASS TRANSFER DURING THE MADE-2 EXPERIMENT
journal, February 2008

  • Guan, Jianyong; Molz, Fred J.; Zhou, Quanlin
  • Water Resources Research, Vol. 44, Issue 2
  • DOI: 10.1029/2007WR006120

Tracer experiments in fractured rocks: Matrix diffusion and the validity of models
journal, August 1993

  • Maloszewski, Piotr; Zuber, Andrzej
  • Water Resources Research, Vol. 29, Issue 8
  • DOI: 10.1029/93WR00608

Numerical Simulation of Water-Oil Flow in Naturally Fractured Reservoirs
journal, December 1976

  • Kazemi, H.; Merrill, L. S.; Porterfield, K. L.
  • Society of Petroleum Engineers Journal, Vol. 16, Issue 06
  • DOI: 10.2118/5719-PA

Dead-End Pore Volume and Dispersion in Porous Media
journal, March 1964

  • Coats, K. H.; Smith, B. D.
  • Society of Petroleum Engineers Journal, Vol. 4, Issue 01
  • DOI: 10.2118/647-PA

Contaminant transport in fractured porous media: Analytical solutions for a system of parallel fractures
journal, December 1982


CO2-H2O mixtures in the geological sequestration of CO2. I. Assessment and calculation of mutual solubilities from 12 to 100°C and up to 600 bar
journal, August 2003

  • Spycher, Nicolas; Pruess, Karsten; Ennis-King, Jonathan
  • Geochimica et Cosmochimica Acta, Vol. 67, Issue 16, p. 3015-3031
  • DOI: 10.1016/S0016-7037(03)00273-4

Interpretation of injection–withdrawal tracer experiments conducted between two wells in a large single fracture
journal, September 2004


CO2-H2O mixtures in the geological sequestration of CO2. II. Partitioning in chloride brines at 12–100°C and up to 600 bar
journal, July 2005


On the theory of tracer experiments in fissured rocks with a porous matrix
journal, July 1985


Matrix-fracture transfer shape factors for dual-porosity simulators
journal, November 1995


Multi-rate mass transfer modeling of two-phase flow in highly heterogeneous fractured and porous media
journal, May 2016


The Grimsel (Switzerland) migration experiment: integrating field experiments, laboratory investigations and modelling
journal, February 1996


Accurate early-time and late-time modeling of countercurrent spontaneous imbibition: ACCURATE MODELING OF SPONTANEOUS IMBIBITION
journal, August 2016

  • March, Rafael; Doster, Florian; Geiger, Sebastian
  • Water Resources Research, Vol. 52, Issue 8
  • DOI: 10.1002/2015WR018456

A general real-time formulation for multi-rate mass transfer problems
journal, January 2009


Multiple-Rate Mass Transfer for Modeling Diffusion and Surface Reactions in Media with Pore-Scale Heterogeneity
journal, June 1995

  • Haggerty, Roy; Gorelick, Steven M.
  • Water Resources Research, Vol. 31, Issue 10
  • DOI: 10.1029/95WR10583

A Semi-Analytical Solution for Large-Scale Injection-Induced Pressure Perturbation and Leakage in a Laterally Bounded Aquifer–Aquitard System
journal, October 2008

  • Zhou, Quanlin; Birkholzer, Jens T.; Tsang, Chin-Fu
  • Transport in Porous Media, Vol. 78, Issue 1
  • DOI: 10.1007/s11242-008-9290-0

CO2 migration and pressure evolution in deep saline aquifers
journal, September 2015

  • Birkholzer, Jens T.; Oldenburg, Curtis M.; Zhou, Quanlin
  • International Journal of Greenhouse Gas Control, Vol. 40
  • DOI: 10.1016/j.ijggc.2015.03.022

Contaminant transport in fractured porous media: Analytical solution for a single fracture
journal, June 1981

  • Tang, D. H.; Frind, E. O.; Sudicky, E. A.
  • Water Resources Research, Vol. 17, Issue 3
  • DOI: 10.1029/WR017i003p00555

Numerical simulations on the influence of matrix diffusion to carbon sequestration in double porosity fissured aquifers
journal, July 2009


Mass Transfer Studies in Sorbing Porous Media I. Analytical Solutions1
journal, January 1976


Transport upscaling in heterogeneous aquifers: What physical parameters control memory functions?: HETEROGENEITY EFFECTS ON MEMORY FUNCTION
journal, December 2008

  • Willmann, M.; Carrera, J.; Sánchez-Vila, X.
  • Water Resources Research, Vol. 44, Issue 12
  • DOI: 10.1029/2007WR006531

Exact integral solutions for two-phase flow
journal, March 1990


Evaluation of the applicability of the dual-domain mass transfer model in porous media containing connected high-conductivity channels: MASS TRANSFER IN POROUS MEDIA
journal, December 2007

  • Liu, Gaisheng; Zheng, Chunmiao; Gorelick, Steven M.
  • Water Resources Research, Vol. 43, Issue 12
  • DOI: 10.1029/2007WR005965

Modeling the transport of solutes influenced by multiprocess nonequilibrium
journal, September 1989

  • Brusseau, M. L.; Jessup, R. E.; Rao, P. S. C.
  • Water Resources Research, Vol. 25, Issue 9
  • DOI: 10.1029/WR025i009p01971

Upscaling solute transport in naturally fractured porous media with the continuous time random walk method: UPSCALING TRANSPORT IN FRACTURED ROCK WITH CTRW
journal, December 2010

  • Geiger, S.; Cortis, A.; Birkholzer, J. T.
  • Water Resources Research, Vol. 46, Issue 12
  • DOI: 10.1029/2010WR009133

Modeling Basin- and Plume-Scale Processes of CO2 Storage for Full-Scale Deployment
journal, December 2009


Basic concepts in the theory of seepage of homogeneous liquids in fissured rocks [strata]
journal, January 1960

  • Barenblatt, G. I.; Zheltov, Iu. P.; Kochina, I. N.
  • Journal of Applied Mathematics and Mechanics, Vol. 24, Issue 5
  • DOI: 10.1016/0021-8928(60)90107-6

A Practical Method for Modeling Fluid and Heat Flow in Fractured Porous Media
journal, February 1985

  • Pruess, Karsten
  • Society of Petroleum Engineers Journal, Vol. 25, Issue 01
  • DOI: 10.2118/10509-PA

Evidence of Multi-Process Matrix Diffusion in a Single Fracture from a Field Tracer Test
journal, June 2006

  • Zhou, Quanlin; Liu, Hui-Hai; Bodvarsson, Gudmundur S.
  • Transport in Porous Media, Vol. 63, Issue 3
  • DOI: 10.1007/s11242-005-1123-9

On the late-time behavior of tracer test breakthrough curves
journal, December 2000

  • Haggerty, Roy; McKenna, Sean A.; Meigs, Lucy C.
  • Water Resources Research, Vol. 36, Issue 12
  • DOI: 10.1029/2000WR900214

Tracer tests in a fractured dolomite: 1. Experimental design and observed tracer recoveries
journal, May 2001

  • Meigs, Lucy C.; Beauheim, Richard L.
  • Water Resources Research, Vol. 37, Issue 5
  • DOI: 10.1029/2000WR900335

A method for quick assessment of CO2 storage capacity in closed and semi-closed saline formations
journal, October 2008

  • Zhou, Quanlin; Birkholzer, Jens T.; Tsang, Chin-Fu
  • International Journal of Greenhouse Gas Control, Vol. 2, Issue 4
  • DOI: 10.1016/j.ijggc.2008.02.004

On matrix diffusion: formulations, solution methods and qualitative effects
journal, June 1998

  • Carrera, Jesús; Sánchez-Vila, Xavier; Benet, Inmaculada
  • Hydrogeology Journal, Vol. 6, Issue 1
  • DOI: 10.1007/s100400050143

Tracer tests in a fractured dolomite: 2. Analysis of mass transfer in single-well injection-withdrawal tests
journal, May 2001

  • Haggerty, Roy; Fleming, Sean W.; Meigs, Lucy C.
  • Water Resources Research, Vol. 37, Issue 5
  • DOI: 10.1029/2000WR900334

What controls the apparent timescale of solute mass transfer in aquifers and soils? A comparison of experimental results: CONTROLLING THE TIMESCALE OF SOLUTE MASS TRANSFER
journal, January 2004

  • Haggerty, Roy; Harvey, Charles F.; Freiherr von Schwerin, Claudius
  • Water Resources Research, Vol. 40, Issue 1
  • DOI: 10.1029/2002WR001716

A time-convolution approach for modeling heat exchange between a wellbore and surrounding formation
journal, December 2011