In this paper, vacancy clustering formation and diffusion of germanium substrate was studied. The analysis method was adopted Monte Carlo method. At temperatures higher than melting point, fewer clusters formed, but there was less variation in the number of clusters than at lower temperatures, as the time increased. Equilibrium diffusivities in the clustering region were $10^2$ lower than those of free vacancies in the initial stage of kinetic lattice Monte Carlo simulations. They were expressed according to three temperature regimes: at temperatures above 1,100 K, at temperatures of 1,100-900 K, and at temperatures below 900 K. The effective mean migration energy, 1.1 eV, closely coincided with that of the 1.0-1.2 eV in experiments.