Results of computational-analytic study of heat exchange processes in glass furnace inside are presented. Calculations for heating with four burners are performed using a three dimensional zonal model of heat exchange for three versions of heat load distribution by the burners as well as for the furnace heating with natural gas sprays.