Motivated by the central phenomenon of energy cascades in wave turbulence theory, we construct non-equilibrium statistical steady states (NESS), or invariant measures, for a simplified model derived from the nonlinear Schrödinger (NLS) equation with external forcing and dissipation. This new perspective to studying energy cascades, distinct from traditional analyses based on kinetic equations and their cascade spectra, focuses on the underlying statistical steady state that is expected to hold when the cascade spectra of wave turbulence manifest. In the full generality of the (infinite dimensional) nonlinear Schrödinger equation, constructing such invariant measures is more involved than the rigorous justification of the Kolmogorov-Zakharov (KZ) spectra, which itself remains an outstanding open question despite the recent progress on mathematical wave turbulence. Since such complexity remains far beyond the current knowledge (even for much simpler chain models), we confine our analysis to a three-mode reduced system that captures the resonant dynamics of the NLS equation, offering a tractable framework for constructing the NESS. For this, we introduce a novel approach based on solving an elliptic Feynman-Kac equation to construct the needed Lyapunov function.