The batch population balance model (PBM) was solved using QMOM-MOCH combinational method with specified assumptions for reconstruction of whole crystal size distribution (CSD) of SAPO-34 molecular sieve. Numerical solution of this model was done using programming in MATLAB software environment. Also during this modeling kinetics parameters of secondary nucleation and size-dependent growth were identified by definition of an optimization problem and with help of experimental data about the SAOP-34 CSD and supersaturation solution concentration. The result which that the PBM associated with secondary nucleation and size-dependent growth rates presents an acceptable prediction from dynamic evolution of SAPO-34 CSD at the crystallization time.