model-master.zip_python_综合能源系统

# _*_ coding: utf-8 _*_ from cn.modelICE.Parameters import Parameters from cn.modelICE.model.InternalCombustionEngine import InternalCombustionEngine from cn.modelICE.model.CHP import CHPInternalCombustionEngine from cn.modelICE.model.ColdStorage import ColdStorage from cn.modelICE.model.EleStorage import EleStorage from cn.modelICE.model.GasBoiler import GasBoiler from cn.modelICE.model.HeatPump import HeatPump from cn.modelICE.model.HeatStorage import HeatStorage from cn.modelICE.util.DemandData import DemandData class SeasonColdBL: def __init__(self, temporary, number): chp = CHPInternalCombustionEngine(temporary, number) heat_pump = HeatPump(temporary) elestorage = EleStorage(temporary) coldstorage = ColdStorage(temporary) ice = InternalCombustionEngine(number, temporary) demand = DemandData() cold_stor = 0 ele_stor = 0 self.fuel = [] self.ele_bought = [] self.ele_ice = [] self.cold_absorption_chiller = [] self.cold_heat_pump = [] self.cold_waste = [] self.ele_waste = [] for t in range(7, 23, Parameters.delttime): ele = ice.nominal cold = chp.cold_out_max if demand.C[t] > cold + coldstorage.get_C_out_max(cold_stor) + heat_pump.nominal: self.judge = 0 break else: self.judge = 1 if demand.C[t] <= cold: cold_stor = coldstorage.get_S(cold_stor, cold - demand.C[t], 0) cold_heat_pump = 0 ele_needed_for_heat_pump = 0 elif (demand.C[t] > cold) & (demand.C[t] <= (cold + coldstorage.get_C_out_max(cold_stor))): coldstorage_cold_out = coldstorage.get_C_out_max(cold_stor) + cold - demand.C[t] cold_stor = coldstorage.get_S(cold_stor, 0, coldstorage_cold_out) cold_heat_pump = 0 ele_needed_for_heat_pump = 0 else: coldstorage_cold_out = coldstorage.get_C_out_max(cold_stor) cold_heat_pump = demand.C[t] - coldstorage.get_C_out_max(cold_stor) - cold ele_needed_for_heat_pump = heat_pump.get_E_in(cold_heat_pump) cold_stor = coldstorage.get_S(cold_stor, 0, coldstorage_cold_out) if demand.cold_E[t] <= ele: ele_bought_for_ele = 0 if ele_needed_for_heat_pump > 0: if ele - demand.cold_E[t] >= ele_needed_for_heat_pump: elestorage_ele_in = ele - demand.cold_E[t] - ele_needed_for_heat_pump ele_bought_for_heat_pump = 0 else: elestorage_ele_in = 0 ele_bought_for_heat_pump = ele_needed_for_heat_pump - (ele - demand.cold_E[t]) else: elestorage_ele_in = ele - demand.cold_E[t] ele_bought_for_heat_pump = 0 ele_stor = elestorage.get_S(ele_stor, elestorage_ele_in, 0) elif (demand.cold_E[t] > ele) & (demand.cold_E[t] <= (ele + elestorage.get_E_out_max(ele_stor))): ele_bought_for_ele = 0 ele_bought_for_heat_pump = ele_needed_for_heat_pump elestorage_ele_out = ele + elestorage.get_E_out_max(ele_stor) - demand.cold_E[t] ele_stor = elestorage.get_S(ele_stor, 0, elestorage_ele_out) else: ele_bought_for_heat_pump = ele_needed_for_heat_pump ele_bought_for_ele = demand.cold_E[t] - ele - elestorage.get_E_out_max(ele_stor) elestorage_ele_out = elestorage.get_E_out_max(ele_stor) ele_stor = elestorage.get_S(ele_stor, 0, elestorage_ele_out) if cold_stor > coldstorage.nominal: self.cold_waste.append(cold_stor - coldstorage.nominal) cold_stor = coldstorage.nominal else: self.cold_waste.append(0) if ele_stor > elestorage.nominal: self.ele_waste.append(ele_stor - elestorage.nominal) ele_stor = elestorage.nominal else: self.ele_waste.append(0) self.ele_ice.append(ele) self.cold_absorption_chiller.append(cold) self.cold_heat_pump.append(cold_heat_pump) self.ele_bought.append(ele_bought_for_heat_pump + ele_bought_for_ele) self.fuel.append(ice.get_fuel(1)) if self.judge == 1: for t in range(23, 24, Parameters.delttime): # 把以下cold，ele，heat部分写成class cold_during_shut_down = ColdDuringShutDown(temporary, t, cold_stor) ele_during_shut_down = EleDuringShutDown(temporary, t, ele_stor, 0) cold_stor = cold_during_shut_down.cold_stor ele_stor = ele_during_shut_down.ele_stor self.ele_ice.append(0) self.cold_absorption_chiller.append(0) self.cold_heat_pump.append(cold_during_shut_down.cold_heat_pump) self.ele_bought.append(ele_during_shut_down.ele_bought_for_ele + cold_during_shut_down.ele_bought_for_heat_pump) self.fuel.append(0) for t in range(0, 7, Parameters.delttime): cold_during_shut_down = ColdDuringShutDown(temporary, t, cold_stor) ele_during_shut_down = EleDuringShutDown(temporary, t, ele_stor, 0) cold_stor = cold_during_shut_down.cold_stor ele_stor = ele_during_shut_down.ele_stor self.ele_ice.append(0) self.cold_absorption_chiller.append(0) self.cold_heat_pump.append(cold_during_shut_down.cold_heat_pump) self.ele_bought.append(ele_during_shut_down.ele_bought_for_ele + cold_during_shut_down.ele_bought_for_heat_pump) self.fuel.append(0) class SeasonHeatBL: # 注意不同！！！无热蒸汽需求，只有空间热和热水,exhaust gas进入锅炉制热，jacket water直接供热 def __init__(self, temporary, number): chp = CHPInternalCombustionEngine(temporary, number) ice = InternalCombustionEngine(number, temporary) gas_boiler = GasBoiler(temporary) elestorage = EleStorage(temporary) heatstorage = HeatStorage(temporary) demand = DemandData() heat_stor = 0 ele_stor = 0 self.fuel = [] self.ele_bought = [] self.ele_ice = [] self.heat_gas_boiler = [] self.heat_waste = [] self.ele_waste = [] for t in range(7, 23, Parameters.delttime): if demand.H[t] > chp.heat_out_max + heatstorage.get_H_out_max(heat_stor) + gas_boiler.nominal: self.judge = 0 break else: self.judge = 1 ele = ice.nominal heat = chp.heat_out_max if demand.H[t] <= heat: heatstorage_heat_in = heat - demand.H[t] heatstorage_heat_out = 0 heat_gas_boiler = 0 elif(demand.H[t] > heat) & (demand.H[t] <= heat + heatstorage.get_H_out_max(heat_stor)): heatstorage_heat_in = 0 heatstorage_heat_out = heatstorage.get_H_out_max(heat_stor) + heat - demand.H[t] heat_gas_boiler = 0 else: heatstorage_heat_in = 0 heatstorage_heat_out