建立二维非定常Stokes方程的时间二阶精度的稳定化Crank-Nicolson (CN)有限体积元格式, 并给出其稳定化CN有限体积元解的误差估计。数值实验说明时间二阶精度的稳定化CN有限体积元格式比时间一阶精度格式更优越, 从而表明稳定化CN有限体积元格式对于求解非定常Stokes方程的数值解是有效可行的。