Cluster variational method for critical concentration calculations of site-diluted Ising model