The Effect of Carbonate Cement Type on Sandstone Matrix Acidizing

Document Type : Research Paper


Department of Petroleum Engineering, Faculty of Chemical Engineering, Tarbiat Modares University, Tehran, Iran


This study aims to numerically determine the roles of the geochemical reactions during the injection of a strong acid into a sandstone sample. As a case study, we used laboratory results of hydrochloric acid (HCl) injection into a sandstone core plug sample from the literature. As the exact cement composition of the implemented sandstone was not available, two probable cement compositions were considered (i.e., calcite and dolomite cement). A fully-implicit model, coded in Python, was used to simulate the underlying geochemical reactions during the HCl injection (i.e., equilibrium and kinetical reactions). In addition, the reactive surface area and porosity-permeability changes of the rock sample were included in the model. The modelling results show that dolomite cement matched better than calcite cement with the experimental acidizing data. A perfect effluent pH prediction was therefore achieved when the reactive surface area was considered as a function of mineral volume fraction. Moreover, a detailed analysis of the dissolution/precipitation rate of different minerals involved in simulations was provided. The presented model improves our understanding of sandstone acidizing by determining dominant reactions.


  1. Sadeghnejad S, Enzmann F, Kersten M. Digital rock physics, chemistry, and biology: challenges and prospects of pore-scale modelling approach. Applied Geochemistry. 2021;131:105028.
  2. Khojastehmehr M, Bazargan M. A new geochemical reactive transport model for sandstone acidizing. Computers & Geosciences. 2022;166:105178.
  3. Hou B, Qui K, Chen M, Jin Y, Chen K. The wellbore collapse on sandstone formation during well test with matrix acidizing treatment. Petroleum science and technology. 2013;31(3):237-49.
  4. Li X, Gomaa A, Nino-Penaloza A, Chaudhary S, editors. Integrated Carbonate Matrix Acidizing Model for Optimal Treatment Design and Distribution in Long Horizontal Wells. SPE Production and Operations Symposium; 2015: Society of Petroleum Engineers.
  5. Dong R, Wang Q, Wheeler MF, editors. Prediction of mechanical stability of acidizing-induced wormholes through coupled hydro-chemo-mechanical simulation. 53rd US Rock Mechanics/Geomechanics Symposium; 2019: OnePetro.
  6. Shafiq MU, Mahmud HB. Sandstone matrix acidizing knowledge and future development. Journal of Petroleum Exploration and Production Technology. 2017;7(4):1205-16.
  7. Elakneswaran Y, Takeya M, Ubaidah A, Shimokawara M, Okano H, Nawa T, editors. Integrated geochemical modeling of low salinity waterflooding for enhanced oil recovery in carbonate reservoir. International Petroleum Technology Conference; 2020: OnePetro.
  8. Khurshid I, Al-Shalabi EW, Afgan I, Al-Attar H. A numerical approach to investigate the impact of acid-asphaltene sludge formation on wormholing during carbonate acidizing. Journal of Energy Resources Technology. 2022;144(6).
  9. Choi SK, Ermel YM, Bryant SL, Huh C, Sharma MM, editors. Transport of a pH-sensitive polymer in porous media for novel mobility-control applications. SPE/DOE Symposium on Improved Oil Recovery; 2006: Society of Petroleum Engineers.
  10. Benson IP, Nghiem LX, Bryant SL, Huh C, editors. Development and use of a simulation model for mobility/conformance control using a pH-sensitive polymer. SPE Annual Technical Conference and Exhibition; 2007: Society of Petroleum Engineers.
  11. Kazempour M, Alvarado V. Geochemically based modeling of ph-sensitive polymer injection in berea sandstone. Energy & Fuels. 2011;25(9):4024-35.
  12. Blanc P, Lassin A, Piantone P, Azaroual M, Jacquemet N, Fabbri A, et al. Thermoddem: A geochemical database focused on low temperature water/rock interactions and waste materials. Applied Geochemistry. 2012;27(10):2107-16.
  13. Parkhurst DL, Appelo C. Description of input and examples for PHREEQC version 3: a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations. US Geological Survey; 2013. Report No.: 2328-7055.
  14. Balashov VN, Guthrie GD, Hakala JA, Lopano CL, Rimstidt JD, Brantley SL. Predictive modeling of CO2 sequestration in deep saline sandstone reservoirs: Impacts of geochemical kinetics. Applied Geochemistry. 2013;30:41-56.
  15. Steefel C, Appelo C, Arora B, Jacques D, Kalbacher T, Kolditz O, et al. Reactive transport codes for subsurface environmental simulation. Computational Geosciences. 2015;19(3):445-78.
  16. Younesian-Farid H, Sadeghnejad S. Geochemical performance evaluation of pre-flushing of weak and strong acids during pH-triggered polymer flooding. Journal of Petroleum Science and Engineering. 2019;174:1022-33.
  17. Koochakzadeh A, Younesian-Farid H, Sadeghnejad S. Acid pre-flushing evaluation before pH-sensitive microgel treatment in carbonate reservoirs: Experimental and numerical approach. Fuel. 2021;297:120670.
  18. Altree-Williams A, Brugger J, Pring A, Bedrikovetsky P. Coupled reactive flow and dissolution with changing reactive surface and porosity. Chemical Engineering Science. 2019;206:289-304.
  19. Steefel CI. CrunchFlow. Softw. Model. Multicomponent React. Flow Transp. User's Man. . Lawrence Berkeley Natl Lab, Berkeley USA. 2009.