Bi-level optimization model for integrated energy system considering the thermal comfort of heat customers. (15th December 2018)