New Model to Help Predict CO₂ Reservoir Stability
First large‐scale simulation of reactive transport at the pore scale
August 5, 2013
Contact: Margie Wylie, email@example.com, +1 510 486 7421
Geologic carbon sequestration—trapping carbon dioxide (CO2) from fossil-fuel-burning power plants in underground aquifers—is a key strategy for mitigating global warming. Developing tools and techniques to predict and control the long-term stability of these CO2 reservoirs is the mission of the DOE Energy Frontier Research Center (EFRC) for Nanoscale Control of Geologic CO2 (NCGC), a collaborative effort led by Lawrence Berkeley National Laboratory, and including three other national laboratories and four universities.
As CO2 is pumped underground, it can react chemically with underground minerals and brine in various ways, sometimes resulting in mineral dissolution and precipitation, which can change the porous structure of the aquifer. But predicting these changes is difficult because these processes take place at the pore scale and cannot be calculated using macroscopic models. The dissolution rates of many minerals have been found to be slower in the field than those measured in the laboratory; understanding this discrepancy requires modeling the pore-scale interactions between reaction and transport processes, then scaling them up to reservoir dimensions.
To meet this challenge, the NCGC team combined the SCIDAC‐developed Chombo platform for direct numerical simulation of flow and transport with the complex geochemistry code CrunchFlow. With the integrated Chombo-Crunch code, and using a NERSC NISE allocation on Hopper, they performed the first‐ever large‐scale simulation of reactive transport pore‐scale processes on real‐pore‐space geometry as obtained from experimental data.
Textures affecting flow, transport, and reaction—as might be found in porous rock—were explicitly introduced into the pore-scale model of a capillary tube packed with crushed calcite, with full resolution of all interfaces at the micron scale. This was the largest pore-scale reactive flow simulation ever attempted, involving 1.6 billion grid points, and requiring 8.7 million processor hours on Hopper and up to 138,000 cores. The grid resolution of the model is now very close to that of imaging experiments, enabling a key EFRC goal of combining simulation and experiment.
The new high-resolution model demonstrated that the mineral dissolution rate depends on the pore structure of the aquifer: the average reaction rate depends on the flow pattern within the porous rock and the concentration gradients that develop at the pore scale. Lack of uniformity in the flow field at the pore scale (due to varying rock structures) has the effect of decreasing the overall reactivity of the system, and this effect becomes more pronounced as the variety of rock structures increases and the fluid flow slows down. This result is the first step in explaining the discrepancy between laboratory experiments and field measurements of mineral dissolution rates.
The integration of microscopic imaging, pore-scale modeling, and carefully controlled experimentation will provide an unprecedented new understanding of how fundamental molecular-scale processes of flow, transport, and reaction play out at the pore scale. New insights into the pore-scale dynamics that merge under the far-from-equilibrium conditions associated with the injection of CO2 will greatly improve researchers’ ability to develop scientifically defensible predictive models of CO2.
Publication: S. Molins, D. Trebotich, C. I. Steefel, and C. Shen, “An Investigation of the Effect of Pore Scale Flow on Average Geochemical Reaction Rates using Direct Numerical Simulation,” Water Resources Research 48, W03527 (2012), doi:10.1029/2011WR011404.
About Computing Sciences at Berkeley Lab
The Lawrence Berkeley National Laboratory (Berkeley Lab) Computing Sciences organization provides the computing and networking resources and expertise critical to advancing the Department of Energy's research missions: developing new energy sources, improving energy efficiency, developing new materials and increasing our understanding of ourselves, our world and our universe.
ESnet, the Energy Sciences Network, provides the high-bandwidth, reliable connections that link scientists at 40 DOE research sites to each other and to experimental facilities and supercomputing centers around the country. The National Energy Research Scientific Computing Center (NERSC) powers the discoveries of 6,000 scientists at national laboratories and universities, including those at Berkeley Lab's Computational Research Division (CRD). CRD conducts research and development in mathematical modeling and simulation, algorithm design, data storage, management and analysis, computer system architecture and high-performance software implementation. NERSC and ESnet are DOE Office of Science User Facilities.
Lawrence Berkeley National Laboratory addresses the world's most urgent scientific challenges by advancing sustainable energy, protecting human health, creating new materials, and revealing the origin and fate of the universe. Founded in 1931, Berkeley Lab's scientific expertise has been recognized with 13 Nobel prizes. The University of California manages Berkeley Lab for the DOE’s Office of Science.
DOE’s Office of Science is the single largest supporter of basic research in the physical sciences in the United States, and is working to address some of the most pressing challenges of our time. For more information, please visit science.energy.gov.