We present observations of the ground state 1$_{0,1}$-0$_{0,0}$ rotational transition of HDO at 464.925 GHz and the 1$_{1,0}$-1$_{0,1}$ transition at 509.292 GHz toward the three high-mass star forming regions: G34.26+0.15, W49N, and W51e$_{1}$/e$_{2}$, carried out with the Caltech Submillimeter Observatory. The latter transition is observed for the first time from the ground. The spectra are modeled, together with observations of higher-energy HDO transitions, as well as submillimeter dust continuum fluxes from the literature, using a spherically symmetric radiative transfer model to derive the radial distribution of the HDO abundance in the target sources. The abundance profile is divided into an inner hot core region, with kinetic temperatures higher than 100~K, and a cold outer envelope with lower kinetic temperatures. The derived HDO abundance with respect to H$_2$ is (0.3-3.7)$\times 10^{-8}$ in the hot inner region ($T > 100 \mathrm{K}$) and (7.0-10.0)$\times 10^{-11}$ in the cold outer envelope. We also used two H$_{2}^{18}$O fundamental transitions to constrain the H$_{2}$O abundances in the outer envelopes. The HDO/H$_{2}$O ratios in these cold regions are found to be (1.8-3.1)$\times 10^{-3}$ and are consequently higher than in the hot inner regions of these sources.