Laguerre-Gaussian (LG) beams are characterized by an azimuthal index or topological charge (m), associated with the orbital angular momentum, and by a radial index (p), which represents the number of the rings in the intensity distribution. We present a detailed, systematic study of the first-order phase statistics of the speckle fields created when LG beams of different order interact with random phase screens with different optical roughness. The phase properties of the LG speckle fields are studied in both the Fresnel and the Fraunhofer regimes using the equiprobability density ellipse formalism such that analytical expressions can be derived for the phase statistics.