The charge, proton and neutron density distributions along with nuclear properties were calculated by Hartree-Fock approach with Skyrme force interaction for isotopic Pb chain (88 <= N <= 184). The effects of correlation on neutron skin thicknesses by obtaining bulk and surface contributions were analyzed by three different approaches. The occurrence of the nuclei with bubble structure due to central depletion in nucleonic was investigated for Pb isotopes as a function of relative neutron excess I = (N - Z)/A. The important role of the bubble effect in heavy region was explained by the relation between Coulomb and nn-interaction. The single particle energy levels were determined by each j = l -/+ 1/2 state of stable Pb isotopes, as well as charge form factors F-ch(k). Besides, the average neutron-proton delta V-np and residual neutron-proton delta(np) interactions of the Pb isotopes were calculated by the theoretical binding energies. The fluctuations in the obtained results were considered in detail because its variation may give a good criterion for the mass model approaches.