计算曲面积分 $I=\iint_{\Sigma} x z \mathrm{~d} y \mathrm{~d} z+2 z y \mathrm{~d} z \mathrm{~d} x+3 x y \mathrm{~d} x \mathrm{~d} y$, 其中 $\Sigma$ 为曲面 $z=1-x^{2}-\frac{y^{2}}{4}(0 \leqslant z \leqslant 1)$ 的上侧.
$\text{A.}$
$\text{B.}$
$\text{C.}$
$\text{D.}$