14-17 May 2018
New Orleans
US/Central timezone

Modelling and simulation of reactive dissolution during acidization of fractured carbonate rocks

17 May 2018, 11:59
New Orleans

New Orleans

Poster + 3 Minute Pitch MS 2.02: Modeling and simulation of subsurface flow at various scales Parallel 10-H


Dr Liu Piyang


Acidizing is a commonly used stimulation treatment for carbonate reservoir by injecting acid into near-wellbore formation to remedy the damage and to create operational-optimal deep-penetrating and narrow high permeability channels, termed as wormholes, by mineral dissolution. Wormhole formation in un-fractured matrix has been widely studied by numerical modelling, but little work has been done for fractured cases which may cause wormholes to propagate quite differently. In this study, a continuum model with explicit fractures is developed, where the governing equations are discretized by the finite-volume method. This model produces the correct dissolution patterns on a 2-D un-fractured domain discretized by Delaunay triangulation. Then it is used to examine wormhole formation in cases with single, multiple, and complex fractures, including characteristics such as fracture orientation, aperture, connection, and distribution. It is found that neither the volume of acid required nor the dissolution structure is affected by fractures with small aperture. If the fracture aperture is large enough, the dissolution structure is dominated by the fracture distribution at intermediate injection rates but is independent with the fracture distribution at very low or very high injection rates. The presence of fractures and the increase of fracture aperture have no influence on the type of dissolution patterns. The volume of acid required to break through the core decreases with the increase of the fracture aperture, which is particularly noticeable when the dominate wormhole is created. For the fractured porous media of not well-connected fractures, the final dissolution structure of wormhole is determined by the fracture and matrix porosity distributions together. As for the well-connected fractures, the matrix heterogeneity has no contribution to the creation of the wormhole, and the dissolution front propagates along the fracture trace. Finally, for the field-scale application, it is found that the optimum injection rate is almost unchanged from un-fractured medium to fractured medium if the acid-mineral reaction system is same.


[1] Panga M K, Ziauddin M, Balakotaiah V. Two‐scale continuum model for simulation of wormholes in carbonate acidization[J]. AlChE J, 2005, 51(12): 3231-3248.
[2] Kalia N, Balakotaiah V. Effect of medium heterogeneities on reactive dissolution of carbonates[J]. Chem Eng Sci, 2009, 64(2): 376-390.
[3] Maheshwari P, Ratnakar R R, Kalia N, et al. 3-D simulation and analysis of reactive dissolution and wormhole formation in carbonate rocks[J]. Chem Eng Sci, 2013, 90(0): 258-274.
[4] Liu P, Yao J, Couples G D, et al. Numerical modelling and analysis of reactive flow and wormhole formation in fractured carbonate rocks[J]. Chem Eng Sci, 2017, 172: 143-157.
[5] Liu P, Yao J, Couples G D, et al. Modelling and simulation of wormhole formation during acidization of fractured carbonate rocks[J]. Journal of Petroleum Science and Engineering, 2017, 154: 284-301.

Acceptance of Terms and Conditions Click here to agree

Primary authors

Dr Liu Piyang Prof. Yao Jun (China University of Petroleum (East China))

Presentation Materials

There are no materials yet.