(本题满分 10 分) 计算曲面积分 $I = \iint_{\Sigma} xz \, \mathrm{d}y \, \mathrm{d}z + 2zy \, \mathrm{d}z \, \mathrm{d}x + 3xy \, \mathrm{d}x \, \mathrm{d}y$ ，其中 $\Sigma$ 为曲面 $z = 1 - x^2 - \frac{y^2}{4} (0 \leqslant z \leqslant 1)$ 的上侧.