In light of the recent LHC Higgs search data, we investigate the pair production of a SM-like Higgs boson around 125GeV in the MSSM and NMSSM. We first scan the parameter space of each model by considering various experimental constraints, and then calculate the Higgs pair production rate in the allowed parameter space. We find that in most cases the dominant contribution to the Higgs pair production comes from the gluon fusion process and the production rate can be greatly enhanced, maximally 10 times larger than the SM prediction (even for a TeV-scale stop the production rate can still be enhanced by a factor of 1.3). We also calculate the χ 2 value with the current Higgs data and find that in the most favored parameter region the production rate is enhanced by a factor of 1.45 in the MSSM, while in the NMSSM the production rate can be enhanced or suppressed (σ SUSY/σ SM varies from 0.7 to 2.4).