解:设球面上一点$P(x,y,z)$,法向量为$\overrightarrow{n}(\frac{x}{2},\frac{y}{2},z)$
而切线方程为
$\frac{x}{2}X+\frac{y}{2}Y+zZ=1,$
$(0,0,0)$到$P(x,y,z)$的距离为
$d(x,y,x)=\frac{1}{\sqrt{(\frac{x}{2})^2+(\frac{y}{2})^2+z^2}}=\frac{1}{\rho }.$
则
$I=\iint_S\frac{z}{d(x,y,z)}dS=\iint_S z\rho dS.$
而
$\rho =\sqrt{(\frac{x}{2})^2+(\frac{y}{2})^2+z^2}=\frac{x}{2}\frac{\frac{x}{2}}{\rho }+\frac{y}{2}\frac{\frac{y}{2}}{\rho }+z\frac{z}{\rho }=\frac{x}{2}\cos \alpha +\frac{y}{2}\cos \beta +z\cos \gamma .$
故
$I=\iint_S z\rho dS=\iint_Sz(\frac{x}{2}\cos \alpha +\frac{y}{2}\cos \beta +z\cos \gamma)dS.$
补曲面
$S':z=0,\frac{x^2}{2}+\frac{y^2}{2}\leq 1.$取下侧。
而
$\iint_{S'}z(\frac{x}{2}\cos \alpha +\frac{y}{2}\cos \beta +z\cos \gamma)dS=0.$
由高斯公式,
$\begin{align*}\therefore I&=\iint_S\frac{z}{d(x,y,z)}dS\\\\&=\iiint_\Omega (\frac{1}{2}+\frac{1}{2}+2)zdV\\\\&=3\int_{0}^{1}zdx\iint_{D(z)}dxdy\\\\&=3\pi\int_{0}^{1}2z(1-z^2)dz\\\\&=\frac{3}{2}\pi.\end{align*}$