vav.py 11 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253
  1. # -*- coding: utf-8 -*-
  2. from operator import attrgetter
  3. from typing import Dict, List, Tuple
  4. import numpy as np
  5. from fastapi import HTTPException
  6. from httpx import AsyncClient
  7. from loguru import logger
  8. from sqlalchemy.orm import Session
  9. from app.controllers.equipment.controller import EquipmentController
  10. from app.crud.space.weight import get_weights_by_vav
  11. from app.schemas.equipment import VAVBox, FCU
  12. from app.schemas.sapce_weight import SpaceWeight
  13. from app.schemas.space import Space
  14. from app.services.platform import DataPlatformService
  15. from app.services.transfer import Duoduo, SpaceInfoService, Season
  16. from app.utils.date import get_time_str
  17. class VAVController(EquipmentController):
  18. def __init__(self, equipment: VAVBox):
  19. super(VAVController, self).__init__()
  20. self.equipment = equipment
  21. async def get_strategy(self):
  22. strategy = 'Plan A'
  23. for space in self.equipment.spaces:
  24. for eq in space.equipment:
  25. if isinstance(eq, FCU):
  26. strategy = 'Plan B'
  27. break
  28. return strategy
  29. async def build_virtual_temperature(self) -> Tuple[float, float]:
  30. target_list, realtime_list = [], []
  31. buffer_list = []
  32. strategy = await self.get_strategy()
  33. for space in self.equipment.spaces:
  34. if not np.isnan(space.temperature_target):
  35. target_list.append(space.temperature_target)
  36. realtime_list.append(space.realtime_temperature)
  37. if strategy == 'Plan B':
  38. for eq in space.equipment:
  39. if isinstance(eq, FCU):
  40. buffer = (4 - eq.air_valve_speed) / 4
  41. buffer_list.append(buffer)
  42. break
  43. total_target = buffer_list + target_list
  44. total_realtime = buffer_list + realtime_list
  45. if total_target and total_realtime:
  46. target_result = np.array(total_target).sum() / len(target_list)
  47. realtime_result = np.array(total_realtime).sum() / len(realtime_list)
  48. self.equipment.setting_temperature = target_result
  49. else:
  50. target_result, realtime_result = np.NAN, np.NAN
  51. return target_result, realtime_result
  52. async def get_supply_air_flow_set(self, temperature_set: float, temperature_realtime: float) -> float:
  53. if np.isnan(temperature_set) or np.isnan(temperature_realtime):
  54. supply_air_flow_set = 0.0
  55. else:
  56. temperature_supply = self.equipment.supply_air_temperature
  57. if np.isnan(temperature_supply):
  58. temperature_supply = 19.0
  59. ratio = abs(1 + (temperature_realtime - temperature_set) / (temperature_set - temperature_supply))
  60. supply_air_flow_set = self.equipment.supply_air_flow * ratio
  61. supply_air_flow_set = max(self.equipment.supply_air_flow_lower_limit, supply_air_flow_set)
  62. supply_air_flow_set = min(self.equipment.supply_air_flow_upper_limit, supply_air_flow_set)
  63. self.equipment.supply_air_flow_set = supply_air_flow_set
  64. self.equipment.virtual_target_temperature = temperature_set
  65. self.equipment.virtual_realtime_temperature = temperature_realtime
  66. return supply_air_flow_set
  67. async def run(self):
  68. temperature_set, temperature_realtime = await self.build_virtual_temperature()
  69. await self.get_supply_air_flow_set(temperature_set, temperature_realtime)
  70. self.equipment.running_status = True
  71. def get_results(self):
  72. return self.equipment
  73. class VAVControllerV2(VAVController):
  74. def __init__(self, equipment: VAVBox, weights: List[SpaceWeight], season: Season):
  75. super(VAVControllerV2, self).__init__(equipment)
  76. self.weights = weights
  77. self.season = season
  78. async def build_virtual_temperature(self) -> None:
  79. valid_spaces = []
  80. weights = []
  81. for sp in self.equipment.spaces:
  82. if sp.realtime_temperature > 0.0 and sp.temperature_target > 0.0:
  83. valid_spaces.append(sp)
  84. for weight in self.weights:
  85. if weight.space_id == sp.id:
  86. weights.append(weight)
  87. if valid_spaces:
  88. weights = sorted(weights, key=lambda x: x.temporary_weight_update_time)
  89. if weights[-1].temporary_weight_update_time > get_time_str(60 * 60 * 2, flag='ago'):
  90. # If someone has submitted a feedback in past two hours, meet the need.
  91. weight_dic = {weight.space_id: 0.0 for weight in weights}
  92. weight_dic.update({weights[-1].space_id: weights[-1].temporary_weight})
  93. else:
  94. weight_dic = {weight.space_id: weight.default_weight for weight in weights}
  95. total_weight_value = 0.0
  96. for v in weight_dic.values():
  97. total_weight_value += v
  98. if total_weight_value > 0:
  99. weight_dic = {k: v / total_weight_value for k, v in weight_dic.items()}
  100. else:
  101. weight_dic.update({list(weight_dic.keys())[0]: 1.0})
  102. try:
  103. virtual_target, virtual_realtime = 0.0, 0.0
  104. for sp in valid_spaces:
  105. virtual_target += sp.temperature_target * weight_dic.get(sp.id)
  106. virtual_realtime += sp.realtime_temperature * weight_dic.get(sp.id)
  107. except KeyError:
  108. logger.error(f'{self.equipment.id} has wrong vav-space relation')
  109. raise HTTPException(status_code=404, detail='This VAV box has wrong eq-sp relation')
  110. self.equipment.virtual_target_temperature = virtual_target
  111. self.equipment.virtual_realtime_temperature = virtual_realtime
  112. else:
  113. self.equipment.virtual_target_temperature = np.NAN
  114. self.equipment.virtual_realtime_temperature = np.NAN
  115. async def rectify(self) -> Tuple[float, float]:
  116. bad_spaces = list()
  117. for sp in self.equipment.spaces:
  118. if (sp.realtime_temperature > max(27.0, sp.temperature_target) or
  119. sp.realtime_temperature < min(21.0, sp.temperature_target)):
  120. bad_spaces.append(sp)
  121. if bad_spaces:
  122. virtual_diff = self.equipment.virtual_target_temperature - self.equipment.virtual_realtime_temperature
  123. if self.season == Season.cooling:
  124. bad_spaces = sorted(bad_spaces, key=attrgetter('diff'))
  125. worst = bad_spaces[0]
  126. if worst.diff <= 0:
  127. if worst.diff < virtual_diff:
  128. self.equipment.virtual_target_temperature = worst.temperature_target
  129. self.equipment.virtual_realtime_temperature = worst.realtime_temperature
  130. else:
  131. self.equipment.virtual_target_temperature = min(21.0, worst.temperature_target) + 0.5
  132. self.equipment.virtual_realtime_temperature = worst.realtime_temperature
  133. elif self.season == Season.heating:
  134. bad_spaces = sorted(bad_spaces, key=attrgetter('diff'), reverse=True)
  135. worst = bad_spaces[0]
  136. if worst.diff >= 0:
  137. if worst.diff > virtual_diff:
  138. self.equipment.virtual_target_temperature = worst.temperature_target
  139. self.equipment.virtual_realtime_temperature = worst.realtime_temperature
  140. else:
  141. self.equipment.virtual_target_temperature = max(27.0, worst.temperature_target) - 0.5
  142. self.equipment.virtual_realtime_temperature = worst.realtime_temperature
  143. return self.equipment.virtual_target_temperature, self.equipment.virtual_realtime_temperature
  144. async def run(self) -> None:
  145. await self.build_virtual_temperature()
  146. temperature_set, temperature_realtime = await self.rectify()
  147. await self.get_supply_air_flow_set(temperature_set, temperature_realtime)
  148. self.equipment.running_status = True
  149. async def fetch_vav_control_params(project_id: str, equipment_id: str) -> Dict:
  150. async with AsyncClient() as client:
  151. duo_duo = Duoduo(client, project_id)
  152. platform = DataPlatformService(client, project_id)
  153. season = await duo_duo.get_season()
  154. served_spaces = await duo_duo.get_space_by_equipment(equipment_id)
  155. space_objects = []
  156. realtime_supply_air_temperature_list = []
  157. for sp in served_spaces:
  158. sp_id = sp.get('id')
  159. transfer = SpaceInfoService(client, project_id, sp_id)
  160. current_target = await transfer.get_current_temperature_target()
  161. realtime_temperature = await platform.get_realtime_temperature(sp_id)
  162. related_equipment = await transfer.get_equipment()
  163. equipment_objects = []
  164. for eq in related_equipment:
  165. if eq.get('category') == 'ACATFC':
  166. speed = await platform.get_fan_speed(eq.get('id'))
  167. temp_fcu_params = {'id': eq.get('id'), 'air_valve_speed': speed}
  168. fcu = FCU(**temp_fcu_params)
  169. equipment_objects.append(fcu)
  170. elif eq.get('category') == 'ACATAH':
  171. realtime_supply_air_temperature_list.append(
  172. await platform.get_realtime_supply_air_temperature(eq.get('id'))
  173. )
  174. temp_space_params = {
  175. 'id': sp_id,
  176. 'realtime_temperature': realtime_temperature,
  177. 'equipment': equipment_objects,
  178. 'temperature_target': current_target,
  179. 'diff': current_target - realtime_temperature
  180. }
  181. space = Space(**temp_space_params)
  182. space_objects.append(space)
  183. realtime_supply_air_temperature = np.array(realtime_supply_air_temperature_list).mean()
  184. realtime_supply_air_flow = await platform.get_realtime_supply_air_flow(equipment_id)
  185. lower_limit, upper_limit = await platform.get_air_flow_limit(equipment_id)
  186. vav_params = {
  187. 'id': equipment_id,
  188. 'spaces': space_objects,
  189. 'supply_air_temperature': realtime_supply_air_temperature,
  190. 'supply_air_flow': realtime_supply_air_flow,
  191. 'supply_air_flow_lower_limit': lower_limit,
  192. 'supply_air_flow_upper_limit': upper_limit,
  193. 'season': season
  194. }
  195. return vav_params
  196. @logger.catch()
  197. async def get_vav_control_v1(project: str, equipment_id: str) -> VAVBox:
  198. vav_params = await fetch_vav_control_params(project, equipment_id)
  199. vav = VAVBox(**vav_params)
  200. vav_controller = VAVController(vav)
  201. await vav_controller.run()
  202. regulated_vav = vav_controller.get_results()
  203. return regulated_vav
  204. @logger.catch()
  205. async def get_vav_control_v2(db: Session, project_id: str, equipment_id: str) -> VAVBox:
  206. vav_params = await fetch_vav_control_params(project_id, equipment_id)
  207. vav = VAVBox(**vav_params)
  208. weights = get_weights_by_vav(db, equipment_id)
  209. vav_controller = VAVControllerV2(vav, [SpaceWeight.from_orm(weight) for weight in weights], vav_params['season'])
  210. await vav_controller.run()
  211. regulated_vav = vav_controller.get_results()
  212. return regulated_vav