By using similar notation in cylindrical coordinates, the three-dimensional distribution function of focused light field was presented when azimuthally and radially polarized Bessel-Gaussian beam propagated through a high numerical-aperture system with primary spherical aberration. Considering behavior of such beams, the transversal field distributions at focal plane and longitudinal field distributions through focal plane had been simulated at different spherical aberration coefficients. The results show that, in the particular case of a highly focused, azimuthally polarized Bessel-Gaussian beam, the inside radius of annular intensity distribution at the focal plane gradually decreases to a constant with an increase of spherical aberration coefficient, but the outside radius increases. The distance between diffraction focus and Gaussian focus increases more and more, transversal intensity do not take on symmetric distribution about the diffraction focus plane, and the spherical aberration cannot perfectly compensated by defocusing, while the change of light intensity of the radially polarized Bessel-Gaussian beam is similar to that of azimuthally polarized one under different primary spherical aberrations.