In acidizing of carbonate reservoir, because of heterogeneous permeability in the reservoir, injected acid flows preferentially to the high_perm regions and it reacts with the surface of the reservoir which leads to weak stimulation of low_perm regions. Therefore, use of acid diverting technique for uniform stimulation is practical. Injection of viscoelastic surfactant_based self diverting (SDVA) acid is an effective method to divert acid to low_perm regions. Micelles in SDVA acid has the ability that with increasing pH, converts from spherical form to rod form which results in a sudden increase of fluid viscosity. Viscosity increasment performs temporally barrier inside high_perm region and diverts acid to low_perm region. In this study, wormholes distribution is simulated using two_scale continuum model with CFD. The effect of injection rate on dissolution pattern and required acid volume for single core breakthrough using conventional HCl acid and self diverting acid is analyzed and optimum injection rate is aquired. Simulation results of single core shows that when HCl acid is injected, injection pressure reduces uniformly with increasment of permeability due to dissolution while in the case of using self diverting acid, injection pressure is non uniform because of variation of viscosity. For showing diversion of acid, double core is used so that flow diversion-time plot for the two systems is compared together. Simulation results of dissolution pattern of double core shows when SDVA acid used, wormholes penetrates in higher depth in the low_perm regions compared to HCl acid and when breakthrough happens in high_perm regions, diffusion of wormholes in low_perm regions stops for both acid systems. Finally, wormholes distribution is simulated in 3D coordinates for cylinder core, simulation results is compared with experimental resulte, and their compatibility is proven. Keyworde : wells damage, acidizing, wormhole, acid diversion