Correlator expansion approach to stationary states of weakly coupled cavity arrays