Kinetic equation approach for beam thermalization in multimoded nonlinear photonic networks