A numerical method of moments for solute transport in physically and chemically nonstationary formations: linear equilibrium sorption with random Kd