Thermal load prediction in district heating systems