1887

Abstract

Summary

Underground geological storage (UGS) of CO2 in saline aquifer or oil reservoirs is an effective means to reduce CO2 emission at scale. To evaluate these UGS processes and understand the long-term fate of the injected CO2, we need a simulator that can account for multiphase equilibrium involving CO2, speciation reactions in brine, and the reactions with minerals. The calculation algorithms for multiphase geochemical reactions are essential to the robustness and efficiency of such a simulator. We applied the modified RAND method (SPE 182706-PA) to electrolyte systems to calculate phase equilibrium together with speciation reactions and mineral dissolution/precipitation. Modified RAND is a non-stoichiometric approach for simultaneous chemical and phase equilibrium calculation. The method linearizes the species chemical potential and eventually uses the elemental chemical potentials as the main independent variables. This greatly reduces the size of the equations for geochemical systems with many species and reactions. Modified RAND is more structured than the classical methods, for which we need to reselect the independent variables during the calculation to reduce round-off errors, and thus more suitable for UGS in oil reservoirs, where both hydrocarbon phase equilibrium and brine-mineral reactions are important. It is 2nd-order and its solution can be guided by minimizing the Gibbs energy. Modified RAND can be applied directly to geochemical systems at a fixed overall composition. Some geochemical applications, however, require analysis at constant chemical potential of a neutral species (e.g., CO2) or a charged species (e.g., H+), the latter case expressed usually as constant pH. We also extended modified RAND to those open systems. For the former, a new state-function can be constructed through the Legendre transform and the obtained algorithm is an energy minimization. For the latter, the problem is no longer minimization but we can still formulate a 2nd-order convergent algorithm. We tested the modified RAND algorithms with phase equilibrium cases relevant to UGS in closed systems, open systems with specified CO2 fugacity, and open systems with specified pH. Modified RAND provides a more efficient solution than the classical equation solving approach used in PHREEQC. The algorithms for closed and open systems exhibit 2nd-order convergence in all the tested cases. We then integrated modified RAND into a 1-D simulator and included the kinetic reactions, and compared the simulator with PHREEQC for 1-D geochemical simulations. The study provides the foundation for a future reactive transport simulator using modified RAND for the core multiphase reaction calculation.

Loading

Article metrics loading...

/content/papers/10.3997/2214-4609.202035103
2020-09-14
2024-04-23
Loading full text...

Full text loading...

References

  1. Bethke, C. M.
    [2007] Geochemical and Biogeochemical Reaction Modeling.Cambridge University Press, 2 ed.
    [Google Scholar]
  2. Brinkley, S. R.
    [1947] “Calculation of the equilibrium composition of systems of many constituents,” The Journal of Chemical Physics, 15(2), 107–110.
    [Google Scholar]
  3. Castier, M.
    [2009], “Solution of the isochoric-isoenergetic flash problem by direct entropy maximization,” Fluid Phase Equilibria, 276(1), 7–17.
    [Google Scholar]
  4. Greiner, H.
    [1991] “An efficient implementation of Newton’s method for complex non ideal chemical equilibria,” Computers and Chemical Engineering, 15(2), 115–123.
    [Google Scholar]
  5. Jiang, X.
    [2011] “A review of physical modelling and numerical simulation of long-term geological storage of 2,”Applied Energy, 88(11), 3557–3566.
    [Google Scholar]
  6. Kulik, D. A., Wagner, T., Dmytrieva, S. V., Kosakowski, G., Hingerl, F. F., Chudnenko, K. V., and Berner, U. R.
    [2013] “GEM-Selektor geochemical modeling package: Revised algorithm and GEMS3K numerical kernel for coupled simulation codes,” Computational Geosciences, 17(1), 1–24.
    [Google Scholar]
  7. Leal, A. M., Blunt, M. J., and LaForce, T. C.
    [2014] “Efficient chemical equilibrium calculations for geochemical speciation and reactive transport modelling,” Geochimica et Cosmochimica Acta, 131, 301–322.
    [Google Scholar]
  8. [2015] “A chemical kinetics algorithm for geochemical modelling,” Applied Geochemistry, 55, 46–61.
    [Google Scholar]
  9. Liu, P., Zhang, T., and Sun, S.
    [2019], “A tutorial review of reactive transport modeling and risk assessment for geologic CO 2 sequestration,” Computers and Geosciences, 127, 1–11.
    [Google Scholar]
  10. Marini, L.
    [2007] Geological Sequestration of Carbon Dioxide: Thermodynamics, Kinetics, and Reaction Path Modeling.Elsevier, 1s. ed.
    [Google Scholar]
  11. Mathieu, P.
    [2006], The IPCC special report on carbon dioxide capture and storage.
    [Google Scholar]
  12. Medeiros, F. A., Stenby, E. H., and Yan, W.
    [2020] “State Function-Based Flash Specifications for Open Systems in the Absence or Presence of Chemical Reactions,” submitted to the AIChE Journal.
    [Google Scholar]
  13. Michelsen, M. L.
    [1999], “State function based flash specifications,” Fluid Phase Equilibria, 158–160, 617–626.
    [Google Scholar]
  14. Michelsen, M. L., and Mollerup, J. M.
    [2008], Thermodynamic Models: Fundamentals & Computational Aspects.Holte, Denmark: Tie-Line Publications, 2n. ed.
    [Google Scholar]
  15. Parkhurst, D. L., and Appelo, C. A. J.
    [2013] “Description of input and examples for PHREEQC version 3: a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations,” tech. rep., Reston, VA.
    [Google Scholar]
  16. Paterson, D., Michelsen, M. L., Yan, W., and Stenby, E. H.
    [2018], “Extension of modified RAND to multiphase flash specifications based on state functions other than (T,P),”Fluid Phase Equilibria, 458. 288–299.
    [Google Scholar]
  17. Paterson, D.
    [2019] Flash Computation and EoS Modelling for Compositional Thermal Simulation of Flow in Porous Media. Springer Theses, Springer International Publishing.
    [Google Scholar]
  18. Peng, C., Crawshaw, J. P., Maitland, G. C., and Trusler, J. P.
    [2015], “Kinetics of calcite dissolution in CO2-saturated water at temperatures between (323 and 373) K and pressures up to 13.8MPa,”Chemical Geology, 403, 74–85.
    [Google Scholar]
  19. Tsanas, C., Stenby, E. H., and Yan, W.
    [2017] “Calculation of Multiphase Chemical Equilibrium by the Modified RAND Method,” Industrial and Engineering Chemistry Research, 56(41), 11983–11995.
    [Google Scholar]
  20. [2019] “Calculation of multiphase chemical equilibrium in electrolyte solutions with non-stoichiometric methods,” Fluid Phase Equilibria, 482, 81–98.
    [Google Scholar]
http://instance.metastore.ingenta.com/content/papers/10.3997/2214-4609.202035103
Loading
/content/papers/10.3997/2214-4609.202035103
Loading

Data & Media loading...

This is a required field
Please enter a valid email address
Approval was a Success
Invalid data
An Error Occurred
Approval was partially successful, following selected items could not be processed due to error