A physical and mathematical model of the reduced kinetics is presented describing heterogeneous detonation in suspensions non-uniform in particle concentration. The model is based on the heterogeneous media approaches, semi-empirical laws of ignition and combustion, and data on the dependence of the detonation velocity on particle concentration. Formation of suboxides and incomplete combustion of aluminum are taken into account integrally. The dependence of the heat release of chemical reactions and the fraction of unburnt particles on the initial composition is determined from the solution of the stationary problem of the structure of the detonation wave. In the calculations of unsteady detonation flows, it is supposed to solve an additional equation for the spatial distribution of initial concentrations. The problems of initiation and development of cellular detonation in flat channels in suspensions of micron-sized aluminum particles are studied. Dependences of the cell size on particle concentration in uniform suspensions are determined. The flow patterns of cellular structures, the forms of the leading front, and the propagation velocities in channels with longitudinal or transversal gradients of particle concentration are analyzed.