The theory of exponential dispersion models was applied to construct a stochastic model for heterogeneities in regional organ blood flow as inferred from the deposition of labeled microspheres. The requirements that the dispersion model be additive (or reproductive), scale invariant, and represent a compound Poisson distribution, implied that the relative dispersion (RD = standard deviation/mean) of blood flow should exhibit self-similar scaling in macroscopic tissue samples of masses m and m(ref) such that RD(m) = RD(m(ref)). (m/m(ref))(1-D), where D was a constant. Under these circumstances this empirical relationship was a consequence of a compound Poisson-gamma distribution that represented macroscopic blood flow. The model also predicted that blood flow, at the microcirculatory level, should also be heterogeneous but obey a gamma distribution-a prediction supported by observation.