The 1+2D nonlocal nonlinear Schrdinger equation can be transformed to the variational approach in a cylindrical coordinate system, and is applied to a model describing the propagation of optical beam in strongly nonlocal nonlinear media. By solving variational problems with expanding media response functions and assuming reasonable ansatz, the solution of the Laguerre-Gauss form is obtained. The Laguerre-Gaussian beams will form solitons or be reduced to Gaussian beams under certain conditions.