The statistical modeling was introduced to satisfy various experimental conditions on the sorption of heavy metals (Pb, Cu, Cd, and Zn) by clay minerals, i.e. kaolinite, illite and chlorite. The Box-Behnken model designed statistically was applied to determine a relative impact among three variables such as pH, HCO3(or K) concentration and initial concentration of heavy metals. The SAS program was used to obtain the statistical solution by surface response analysis. The results of a statistical sorption modelling indicated that pH is a strong impact of the variables influencing the sorption of heavy metals. A relative effect between an initial concentration of heavy metals and bicarbonate(or K) concentration is dependent on solution condition. The sorption edge of heavy metals as function of pH shows sigmoidal curve, and a great increase in the range of pH 6～8. The sorption sequence among heavy metals is Cu>Pb>>Zn>Cd. The solution chemistry exerts greater influence on the sorption of heavy metals rather than the crystal chemistry of clay minerals. The potassium exerts some effect into a sorption competition with heavy metals. The research suggests that the statistical modeling is an effective method to demonstrate sorption results in three dimension and to reduce the effort of batch sorption experiment.