The statistical model was introduced to satisfy various experimental condition on the sorption of heavy metals (Pb, Cu, Cd, and Zn) by bentonite. The Box-Behnken model designed statistically was applied to determine relative impact among three variables such as pH, HCO$_3$ contents and heavy metal concentrations on the sorption. The SAS program was used to obtain the statistical solution. The statistical surface response analysis indicates that initial concentration of heavy metals and pH have an important effect on the sorption, and bicarbonate is not a serious variable. The sorption capability about heavy metals of bentonite is in the order of Pb>Cu>Zn>Cd. The precipitation as hydroxyl and carbonate complexes of heavy metals was thermodynamically analyzed as major mechanism of sorption under alkaline pHs and high bicarbonate solution. It was found that there is a little difference between the model prediction on the precipitation of heavy metals and the results of batch sorption experiment. The thermodynamic data of the programs have to revise to obtain the best fit condition between the model prediction and the experimental results.