In a previous paper, we studied the Higgs pair production in the standard model with the reaction $e^{+}e^{-}to t bar t HH$. Based on this, we study the Higgs pair production via $e^{+}e^{-}to b bar b HH$. We evaluate the total cross section of $bbar bHH$ and calculate the number total of events considering the complete set of Feynman diagrams at tree-level, and compare this process with the process $e^{+}e^{-}to t bar t HH$. The numerical computation is done for the energy which is expected to be available at a possible Next Linear $e^{+}e^{-}$ Collider with a center-of-mass energy $800, 1000, 1600$ $GeV$ and luminosity 1000 $fb^{-1}$.