Copper and zinc removal from aqueous solution by chelating resin was investigated theoretically and experimentally in the present study. A multistage process was proposed as an alternative for enhancement of the heavy removal of the single-stage process. Heavy metal mass balance equations with empirical Freundlich adsorption isotherm were developed to represent the multistage process and the theoretical model permits determination of the inter-stage heavy metal concentrations and the total amount of chelating resin required for achieving a desired level of heavy metal removal. Optimization of the linearized theoretical model shows that equal division of the total amount of chelating resin among all stages of the multistage process yields the best results in terms of saving of chelating resin for a given heavy metal removal or enhanced heavy metal removal fur a given total amount of chelating resin. Experimental tests were also conducted to establish the equilibrium adsorption of heavy metal by the chelating resin and to empirically verify the advantages of the multistage adsorption process.