Summary. Following a cursory review and synthesis of multipolar continua, the rudiments of graph theory, and granular mechanics, a graph-theoretic, energy-based homogenization is proposed for the systematic derivation of multipolar stress and kinematics in granular media. This provides a weakly non-local hierarchy of multipolar field equations for quasi-static mechanics based on polynomial representations of the kinematics of the type employed in past works. As an improvement on those works, a method is proposed for avoiding "overfitting" of fluctuations based on the so-called "Generalized Additive Method" of statistics. Among other results, it is shown that the standard formula for Cauchy stress in granular media may break down owing to multipolar effects, and that granular rotations in the typical granular medium should not lead to Cosserat effects, as the lowest-order departure from the simple-continuum model.