Statistics of Phase Reduction Models for Improving the Accuracy of the Finite Element Solution of Time-Harmonic Scattering Problems I: General Approach and Low-Order Models