We investigated the neutral Higgs boson pair production at the CERN Large Hadron Collider (LHC) in the SM with four families. We found that the gluon-gluon fusion mode is the most dominant one in producing neutral Higgs boson pair at the LHC, and it can be used to probe the trilinear Higgs coupling. If the heavy quarks of the fourth generation really exist within the SM, they can manifest their effect on the cross section of the Higgs pair production process at the LHC. Our numerical results show that there will be 2 × 104 neutral Higgs boson pair production events per year if the next generation heavy quarks really exist, while there will be only 2 × 103 events produced per year if there are only three families in the SM.