Numerical modeling and analysis of the matrix acidizing process in fractured sandstone rocks with the Extended–FEM

Significant natural fractures develop in deep sandstone reservoirs. However, the previous numerical simulation for matrix acidizing in sandstone rocks always focused on the reactive flow in porous media. A series of responses between multiple types of acids and minerals in fracture and matrix pores have been a significant setback for seeking the mechanism of acid flow in fractured sandstone rocks. This paper, established a multistage reactive-transport model for fractured sandstone rocks based on the two-scale continuum model to explore the effects of fractures on the reactive flow. The weak forms of fluid flow and solute transport equations are derived using the Galerkin method to couple the matrix and fracture domains, and the extended finite element method was used to solve the discretization model. Additionally, we presented numerical simulations under 2D linear flow conditions with specific and sensitive analyses about fracture and matrix properties. Numerical cases showed that the wormhole-shaped structure generated by acid dissolution is hard to develop even in highly heterogeneous fractured sandstone rocks due to the low reaction rate between mud acid and minerals. However, fractures that are not parallel to the flow direction can accelerated acid flow in the formation and reduced Si(OH)₄ precipitation. Increasing the injection rate can not change the dissolution patterns of sandstone, but expanded the acid diffusion range and decreased Si(OH)₄ precipitation in fractured sandstone.