Modeling solute transport in heterogeneous porous media faces two challenges: scale dependence of dispersion and reproducing mixing separately from spreading. Both are crucial since real applications may require km scales whereas reactions, often controlled by mixing, may occur at the pore scale. Methods have been developed in response to these challenges, but none has satisfactorily characterized both processes. In this paper, we propose a formulation based on the Water Mixing Approach extended to account for velocity variability. Velocity is taken as an independent variable, so that concentration depends on time, space and velocity. Therefore, we term the formulation the Multi-Advective Water Mixing Approach. A new mixing term between velocity classes emerges in this formulation. We test it on Poiseuille’s stratified flow using the Water Parcel method. Results show high accuracy of the formulation in both dispersion and mixing. Moreover, the mixing process exhibits Markovianity in space even though it is modeled in time