计算曲面积分
$$
I=\iint_{\Sigma} 2(1-x y) \mathrm{d} y \mathrm{~d} z+(x+1) y \mathrm{~d} z \mathrm{~d} x-4 y z^2 \mathrm{~d} x \mathrm{~d} y,
$$

其中 $\Sigma$ 是弧段 $\left\{\begin{array}{l}z=\sqrt{x-1}, \\ y=0\end{array},(1 \leqslant x \leqslant 3)\right.$ 绕 $x$ 轴旋转一周所得的旋转曲面, $\Sigma$ 上任一点的法向量与 $x$ 轴正向夹角大于 $\frac{\pi}{2}$.