We have developed a fine-scale model of the sandstone core acid flooding process by solving acid and mineral balance equations for a fully three-dimensional flow field that changed as acidizing proceeded. The initial porosity and mineralogy field could be generated in a correlated manner in three dimensions; thus, a laminated sandstone could be simulated. The model has been used to simulate sandstone acidizing coreflood conditions, with a 1 in. diam by 2 in. long core represented by 8000 grid blocks, each having different initial properties. Results from this model show that the presence of small-scale heterogeneities in a sandstone has a dramatic impact on the acidizing process. Flow field heterogeneities cause acid to penetrate much farther into the formation than would occur if the rock were homogeneous, as is assumed by standard models. When the porosity was randomly distributed (sampled from a normal distribution), the acid penetrated up to twice as fast as in the homogeneous case. When the porosity field is highly correlated in the axial direction, which represents a laminated structure, acid penetrates very rapidly into the matrix along the high-permeability streaks, reaching the end of the simulated core as much as 17 times faster than for a homogeneous case.