$f(z)=\frac{1}{(z^2+1)^2}$ behaves like $\frac{1}{z^4}$ for large $z$, hence if $\gamma_R$ is a circle centered at the origin with radius $R$, counter-clockwise oriented,
$$ \lim_{R\to +\infty}\oint_{\gamma_R}f(z)\,dz = 0 \tag{ML} $$
On the other hand $f(z)$ is a meromorphic function with double poles at $z=\pm i$. By Cauchy's integral formula, if $\Gamma_R$ is the closed contour, counter-clockwise oriented, going from $-R$ to $R$ with a straight line, then from $R$ to $-R$ with a half-circle in the upper half-plane, for any $R>1$:
$$ \oint_{\Gamma_R}f(z)\,dz = 2\pi i\cdot\text{Res}\left(f(z),z=i\right)\tag{RES}$$
and since $z=i$ is a double pole for $f(z)=\frac{1}{(z-i)^2 (z+i)^2}$,
$$ \text{Res}\left(f(z),z=i\right) = \lim_{z\to i}\frac{d}{dz}\left((z-i)^2 f(z)\right) = \lim_{z\to i}-\frac{2}{(z+i)^3} = -\frac{i}{4}\tag{DZ} $$
so that $\oint_{\Gamma_R}f(z)\,dz $ equals $\frac{\pi}{2}$ for any $R$ large enough. By $(ML)$, $\frac{\pi}{2}$ is the value of the original integral.
Now, a much shorter, real-analytic proof. For any $a>0$ we have
$$ I(a) = \int_{0}^{+\infty}\frac{1}{x^2+a^2}\,dx = \frac{\pi}{2a}\tag{CauchyPDF}$$
hence by differentiating both sides with respect to $a$ (and exploiting differentiation under the integral sign) we get:
$$ I'(a) = -\int_{0}^{+\infty}\frac{2a}{(x^2+a^2)^2}\,dx = -\frac{\pi}{2a^2}\tag{DUIS} $$
and it is enough to evaluate the last identity at $a=1$.