A three-phase-lag (TPL) model is proposed to describe heat transfer in a finite domain skin tissue with temperature dependent metabolic heat generation. The Laplace transform method is applied to solve the problem. Three special types of heat flux are applied to the boundary of skin tissue for thermal therapeutic applications. The depth of tissue is influenced by the different oscillation heat flux. The comparison between the TPL and dual-phase-lag (DPL) models is analyzed and the effects of phase lag parameters (τq, τt and τv) and material (k*) on the tissue temperature distribution are presented graphically.