In this paper, we study the $\frac{1}{H}$-variation of stochastic divergence integrals $X_t = \int_0^t u_s δB_s$ with respect to a fractional Brownian motion $B$ with Hurst parameter $H < \frac{1}{2}$. Under suitable assumptions on the process u, we prove that the $\frac{1}{H}$-variation of $X$ exists in $L^1(Ω)$ and is equal to $e_H \int_0^T|u_s|^H ds$, where $e_H = \mathbb{E}|B_1|^H$. In the second part of the paper, we establish an integral representation for the fractional Bessel Process $\|B_t\|$, where $B_t$ is a $d$-dimensional fractional Brownian motion with Hurst parameter $H < \frac{1}{2}$. Using a multidimensional version of the result on the $\frac{1}{H}$-variation of divergence integrals, we prove that if $2dH^2 > 1$, then the divergence integral in the integral representation of the fractional Bessel process has a $\frac{1}{H}$-variation equals to a multiple of the Lebesgue measure.