A novel reactive smart tracer method, termed the kinetic interface-sensitive (KIS) tracer test, has been demonstrated in laboratory column experiments to enable measurement of the specific capillary-associated fluid–fluid interfacial area in dynamic two-phase flow displacement processes in porous media. Development of the tracer method towards effective application in real field conditions requires investigation of the influence of the porous media heterogeneity on the front size and the specific interfacial area, and, consequently, in how far a kinetic interface-sensitive tracer experiment, and the corresponding breakthrough curves, are affected. This study employs a two-dimensional Darcy-scale two-phase flow reactive transport model to investigate numerically the KIS tracer transport in heterogeneous porous media. Simulations were carried out for the primary drainage process in a domain formed of fine and coarse porous media. Various heterogeneity patterns, having different numbers of inclusions and different geometrical distributions, were studied. It is shown that the shape of the breakthrough curves can be used as an indicator for quantifying the displacement front roughness, the specific interfacial area in the domain, and the domain heterogeneity, e.g., the existence of preferential flow pathways inside the porous media. The results indicate that when the displacement front roughness is small, the concentration breakthrough curves exhibit a linear increase. The slopes of the breakthrough curves linearly depend on the fraction of the bulk volume occupied by the low-permeability sand inclusions. The volume-averaged specific interfacial area and the size of the transition zone can be determined from the slopes of the breakthrough curves.