Using Richards-Wolf vectorial diffraction integral, three-dimensional light field complex amplitude function was derived when vectorial polarized Bessel-Gaussian beam passed through a high numerical-aperture system with primary coma aberration. In the case of radially and azimuthally polarized Bessel-Gaussian beam, focused optical field was simulated under different comatic aberration coefficients. The results show that the presence of primary coma produces a positional shift of the high intensity lobes and make the intensity distribution transformed in the high focused optical field, and the intensity distribution on focal plane and light maximum along the optical axis depend on primary coma and polarization. In the case of intensity distribution along the optical axis, polarization and primary coma aberration doesn’t change its symmetry.