Renormalization group and variational treatments for a d-dimensional Sine-Gordon model with correlated random fields