计算曲面积分 $I=\iint_{\Sigma} \frac{x \mathrm{~d} y \mathrm{~d} z+y \mathrm{~d} z \mathrm{~d} x+z \mathrm{~d} x \mathrm{~d} y}{\left(x^2+y^2+z^2\right)^{\frac{3}{2}}}$, 其中 $\Sigma$ 是曲面 $1-\frac{z}{7}=\frac{(x-2)^2}{25}+\frac{(y-1)^2}{16}$ $(z \geqslant 0)$ 的上侧.