The aim of this paper is to describe a methodology to simulate scale squeeze treatments in acid stimulated carbonate reservoir. Acid stimulation treatments refer to matrix acidizing, which aims to improve productivity in tight carbonate reservoirs. They generally show a high porosity, but very low permeability. Matrix acidizing treatments have been applied successfully in numerous occasions, they consist of the injection of an acid, generally, hydrochloric acid, mainly because of the high reactivity with carbonate formations. When the acid is injected into the formation, a number of wormholes or highly conductive channels are formed, which are the result of two very distinctive, but interconnected processes. On the one hand, there are chemical reactions between the carbonate minerals and the acid and on the other hand, there are the fluid dynamics of the injected acid, i.e. the fluid loss from wormhole to formation and the fluid distribution in possibly highly complex wormhole geometries. There have been numerous studies investigating the geometry of the wormhole growth pattern in radial and linear laboratory experiments, as well as stochastic simulation. They concluded that wormholes grow in a certain pattern in both axial and angular directions. Although the exact wormhole pattern will strongly depend on the permeability anisotropy and heterogeneity, it is reasonable to assume that the dominant wormholes are expected to grow symmetrically, and that the region dominated by each wormhole is approximately 90° around the wellbore.
To simulate squeeze treatments in an acid stimulated well with a corresponding wormhole pattern, as described above, a reservoir simulator is used. The reservoir simulator describes the pressure field and consequently the propagation of scale inhibitor (SI) along the wormholes, but also into the matrix. The final step is to determine how deep the SI propagates into the matrix, which was used to determine fully a specialized near wellbore model for scale treatment design.