Numerical Simulation of the Wormhole Propagation in Fractured Carbonate Rocks during Acidization Using a Thermal-Hydrologic-Mechanics-Chemical Coupled Model

Acidizing is a widely adopted approach for stimulating carbonate reservoirs. The two-scale continuum (TSC) model is the most widely used model for simulating the reactive process in a carbonate reservoir during acidizing. In realistic cases, there are overburden pressure and pore pressure at present. When the injected acid reacts with the rock, the dissolution of the rock and the consumption of the acid in the pore will break the mechanical balance of the rock. Many experimental studies show that cores after acidizing have lower strength. However, it is still not clear how the deformation of rocks by the change of ground stress influences the acidizing dynamics. For fractured carbonate reservoirs, fractures play a leading role in the flow of injected acid, which preferentially flows into the fractures and dissolves the fracture walls. The effect of the combined action of rock mechanical balance broken and fracture wall dissolution on the formation of wormholes in fractured carbonate reservoirs remains to be studied. To address the above-mentioned issues, a thermal-hydrologic-mechanical-chemical coupled model is presented based on the TSC model for studying the wormhole propagation in fractured carbonate reservoirs under practical conditions. Linear and radial flow cases are simulated to investigate the influences of fracture distribution, reaction temperature, and effective stress on acidizing dynamics. The simulation results show that more wormhole branches are formed by acidizing if the fractures are perpendicular to the flow direction of acid. Temperature is a key parameter affecting the acidification dissolution patterns, so the influence of temperature cannot be ignored during the acidification design. As the effective stress of the formation increases, the diameter of the wormhole gradually decreases, and the branching decreases. More acid is needed for the same stimulation result under higher effective stress.