设 $\Sigma$ 为曲面 $4 x^2+y^2+z^2=1(z \geqslant 0)$ 的下侧, 计算曲面积分
$$
I=\iint_{\Sigma}(x+2 y) \mathrm{d} y \mathrm{~d} z+\frac{z}{\sqrt{x^2+y^2+z^2}} \mathrm{~d} z \mathrm{~d} x+\left(x^2-\frac{y}{\sqrt{x^2+y^2+z^2}}\right) \mathrm{d} x \mathrm{~d} y .
$$
$\text{A.}$
$\text{B.}$
$\text{C.}$
$\text{D.}$