During injection treatments, bottomhole pressure measurements may significantly mismatch modeling results. We devise a computationally effective technique for interpretation of fluid injection in a wellbore interval with multiple geological layers based on the bottomhole pressure measurements. The permeability, porosity and compressibility in each layer are initially setup, while the skin factor and partitioning of injected fluids among the zones during the injection are found as a solution of the problem. The problem takes into account Darcy flow and chemical interactions between the injected acids, diverter fluids and reservoir rock typical in modern matrix acidizing treatments. Using the synchronously recorded injection rate and bottomhole pressure, we evaluate skin factor changes in each layer and actual fluid placement into the reservoir during different pumping jobs: matrix acidizing, water control, sand control, scale squeezes and water flooding. The model is validated by comparison with a simulator used in industry. It gives opportunity to estimate efficiency of a matrix treatment job, role of every injection stage, and control fluid delivery to each layer in real time. The presented interpretation technique significantly improves accuracy of matrix treatments analysis by coupling the hydrodynamic model with records of pressure and injection rate during the treatment.