A model is established to calculate the steam quality distribution along the well bore during steam injection of thermal recovery wells, based on related thermodynamics and hydromechanics principles.