A kinetic Monte Carlo (kMC) simulation with a cluster-multiple labeling technique was carried out to study the dynamic behavior of Si-C clusters during the early stage of epitaxial growth. In the model, a crystal lattice was established to fix the physical location of atoms and interatomic bonding. Events, such as adatoms adsorption, diffusion on the planar surface, and attachment to and detachment from clusters of adatoms were considered in the model. The Si and C atoms were treated individually to achieve more elaborate information for dynamic behavior of the clusters in the crystal surface. Moreover, a cluster multiple labeling technique, also known as the Hoshen - Kopelman algorithm was used to identify Si-C clusters. Then the properties of the clusters were calculated and analyzed by the algorithm. The simulation results showed that the number of Si-C clusters gradually decreases with increasing substrate temperature and decreasing deposition rate, and that the average size of Si-C clusters at high temperatures is much larger than that at low temperatures. Furthermore, the Si:C ratio on the growing surface has a significant effect on the surface morphology. At Si-rich conditions, the cluster density is relatively high, and the cluster shape is not homogeneously distributed. While with the increase of C mole fraction, the cluster density gradually decreases and the cluster shape becomes more uniform. Therefore, a higher concentration of C in deposition flux is beneficial to the early stage of nucleation. These findings are consistent with predictions for the case of the epitaxial growth of SiC.