The structure of neutron-rich Ca isotopes is studied in the spherical Skyrme-Hartree-Fock-Bogoliubov (SHFB) approach with SLy5, SLy5+T, and 36 sets of TIJ parametrizations. The calculated results are compared with the available experimental data for the average binding energies, two-neutron separation energies and charge radii. It is found that the SLy5+T, T31, and T32 parametrizations reproduce best the experimental properties, especially the neutron shell effects at N = 20, 28 and 32, and the recently measured two-neutron separation energy of 56Ca. The calculations with the SLy5+T and T31 parametrizations are extended to isotopes near the neutron drip line. The neutron giant halo structure in the very neutron-rich Ca isotopes is not seen with these two interactions. However, depleted neutron central densities are found in these nuclei. By analyzing the neutron mean-potential, the reason for the bubble-like structure formation is given.