Cluster variation method (CVM) was applied to calculate phase equilibria of metal-hydrogen systems. Two subjects are introduced in the present report. One is the summary of previous studies on the Pd-H system, and it is demonstrated that a single CVM free energy formula can systematically derive information of phase equilibria, intrinsic stability, and short range order diffuse intensities. The second subject is the theoretical calculations of superabundant vacancy (SAV) formation. Within the square approximation of the CVM, it is shown that abundant vacancies are introduced with the absorption of hydrogen when the interaction between vacancy and hydrogen is considered.