The survival analysis methods are usually conducted based on assumption that the population is homogeneity. However, generally, this assumption in most cases is unrealistic, because of unobserved risk factors or subject specific random effect. Disregarding the heterogeneity leads to unbiased results. So frailty model as a mixed model was used to adjust for uncertainty that cannot be explained by observed factors in survival analysis. In this paper, family of power variance function distributions that includes gamma and inverse Gaussian distribution were introduced and evaluated for frailty effects. Finally the proportional hazard frailty models with Weibull baseline hazard as a parametric model used for analyzing survival data of the colorectal cancer patients.