The modeling of heat and fluid flow inside two-phase geothermal wells is a vital task required for the study of the production performance. Gas void fraction is one of the crucial parameters required for a better prediction of pressure and temperature gradients in two-phase geothermal wells. This parameter affects the correct matching between simulated and measured data. Modeling of two-phase flow inside wells is complex because two phases exist concurrently (exhibiting various flow patterns that depend on their relative concentrations, the pipe geometry, and the mass flowrate). A reliable modeling requires the precise knowledge of the two-phase flow patterns (including their transitions and some flow parameters). In this work, ten empirical correlations were used to estimate the gas void fraction in vertical-inclined pipes, and to evaluate their effect on the prediction of two-phase flow characteristics of some Mexican geothermal wells. High quality downhole pressure/ temperature logs collected from four producing geothermal wells were studied [Los Azufres, Mich. (Az-18); Los Humeros, Pue. (H-1), and Cerro Prieto, B.C. (M-90 and M-201)]. The pressure/ temperature gradients were simulated using an improved version of the wellbore simulator GEOPOZO, and the gas void fraction correlations. The simulated results were statistically compared with measured field data.