We present an extensive analysis of the 850 μm polarization maps of the SCUBA Polarimeter Legacy (SCUPOL) Catalogue produced by Matthews et al., focusing exclusively on the molecular clouds and star-forming regions. For the sufficiently sampled regions, we characterize the depolarization properties and the turbulent-to-mean magnetic field ratio of each region. Similar sets of parameters are calculated from two-dimensional synthetic maps of dust-emission polarization produced with three-dimensional magnetohydrodynamics (MHD) numerical simulations scaled to the S106, OMC-2/3, W49, and DR21 molecular cloud polarization maps. For these specific regions, the turbulent MHD regimes retrieved from the simulations, as described by the turbulent Alfvén and Sonic Mach numbers, are consistent within a factor one to two with the values of the same turbulent regimes estimated from the analysis of Zeeman measurements data provided by Crutcher. Constraints on the values of the inclination angle α of the mean magnetic field with respect to the line of sight are also given. The values obtained from the comparison of the simulations with the SCUPOL data are consistent with the estimates made by using two observational methods provided by other authors. Our main conclusion is that simple, ideal, isothermal, and non-self-gravitating MHD simulations are sufficient in order to describe the large-scale observed physical properties of the envelopes of this set of regions.