A three factor, three-level Box-Behnken experimental design combined with response surface modeling (RSM) and quadratic programming (QP) was employed for investigation of lead and copper removal from oil field brine by potential sorption method. Three independent variables (initial pH of solution ranging from 7 to 13, temperature ranging from 20 to 24 degrees C, and reaction time ranging from 20 to 40min) were selected and a second-order polynomial regression equation was then derived to predict responses. The output functions of RSM method were used as an objective function for maximizing average particle size (APS) and minimizing the lead and copper concentration in solution. The particle morphology and different phases of calcium carbonate particles precipitated from pure solution, synthetic solution, and oil field brine were observed by scanning electron microscopy (SEM) and x-ray diffraction pattern (XRD) respectively. The pH=12.68, reaction time=27min, and temperature=24 degrees C were obtained as an optimum condition. Lead and copper concentration in industrial oil field brine were reduced from 2.911 to 0.127ppm (95.63% removal) and 0.476 to 0.025ppm (94.74% removal), respectively, by experiment with optimal condition.