def __init__(self, temporary, number): chp = CHPInternalCombustionEngine(temporary, number) heat_pump = HeatPump(temporary) elestorage = EleStorage(temporary) coldstorage = ColdStorage(temporary) absorption_chiller = DoubleEffectAbsorptionChiller(temporary) ice = InternalCombustionEngine(number, temporary) demand = DemandData() time = demand.E_sheetnrows - 1 cold_stor = 0 ele_stor = 0 fuel = [] ele_bought = [] ele_ice = [] cold_absorption_chiller = [] for t in range(0, time, Parameters.delttime): judge_ele = min(demand.E[t] - elestorage.get_E_out_max(ele_stor), ice.nominal) judge_cold = chp.get_cold_through_ele(judge_ele) if demand.C[t] > coldstorage.get_C_out_max(cold_stor) + judge_cold + heat_pump.nominal: self.judge = 0 break else: self.judge = 1 ele_drive = EleDrive(t, temporary, number, ele_stor) if ele_drive.ele == 0: cold_drive = ColdDrive(t, temporary, number, cold_stor) cold_stor = cold_drive.cold_stor self.cold = cold_drive.cold self.ele = ice.get_ele_out_through_cold(self.cold) pl = self.ele / ice.nominal ele_stor = elestorage.get_S(ele_stor, self.ele, ele_drive.elestorage_out) fuel.append(ice.get_fuel(pl)) ele_bought.append(heat_pump.get_E_in(cold_drive.heat_pump_out)) else: self.ele = ele_drive.ele pl = self.ele / ice.nominal exhaust_gas = ice.get_exhaust_gas_pl(pl) jacket_water = ice.get_jacket_water_pl(pl) self.cold = absorption_chiller.get_cold_out(exhaust_gas, jacket_water) cold_follow = ColdFollow(t, temporary, cold_stor, self.cold) cold_stor = cold_follow.cold_stor ele_stor = elestorage.get_S(ele_stor, 0, ele_drive.elestorage_out) fuel.append(ice.get_fuel(pl)) ele_bought.append(ele_drive.ele_bought + heat_pump.get_E_in(cold_follow.heat_pump_out)) ele_ice.append(self.ele) cold_absorption_chiller.append(self.cold) if cold_stor > coldstorage.nominal: cold_stor = coldstorage.nominal if ele_stor > elestorage.nominal: ele_stor = elestorage.nominal self.fuel = sum(fuel) self.ele_bought = sum(ele_bought) self.emission_calculate_ice = sum(ele_ice) self.emission_calculate_boiler = 0 self.emission_calculate_absorption_chiller = sum(cold_absorption_chiller) self.emission_calculate_grid = self.ele_bought
def __init__(self, temporary, number): chp = CHPInternalCombustionEngine(temporary, number) ice = InternalCombustionEngine(number, temporary) absorption_chiller = DoubleEffectAbsorptionChiller(temporary) gas_boiler = GasBoiler(temporary) elestorage = EleStorage(temporary) heatstorage = HeatStorage(temporary) demand = DemandData() time = demand.E_sheetnrows - 1 heat_stor = 0 ele_stor = 0 fuel = [] ele_bought = [] ele_ice = [] heat_gas_boiler = [] heat_absorption_chiller = [] for t in range(0, time, Parameters.delttime): ele = min(demand.E[t] - elestorage.get_E_out_max(ele_stor), ice.nominal) judge_heat = chp.get_heat_through_ele(ele) if demand.H[t] > heatstorage.get_H_out_max(heat_stor) + judge_heat + gas_boiler.nominal: self.judge = 0 break else: self.judge = 1 ele_drive = EleDrive(t, temporary, number, ele_stor) if ele_drive.ele == 0: heat_drive = HeatDriveChiller(t, temporary, number, heat_stor) self.heat = heat_drive.heat self.heat_gas_boiler = heat_drive.gas_boiler_out self.ele = ice.get_ele_out_through_heat_mode(self.heat) pl = self.ele / ice.nominal self.heat_absorption_chiller = ice.get_exhaust_gas_pl(pl) * absorption_chiller.COP_heat heat_stor = heat_drive.heat_stor ele_stor = elestorage.get_S(ele_stor, self.ele, ele_drive.elestorage_out) fuel.append(ice.get_fuel(pl) + gas_boiler.get_Fuel_in(heat_drive.gas_boiler_out)) ele_bought.append(0) else: self.ele = ele_drive.ele pl = self.ele / ice.nominal self.heat_absorption_chiller = ice.get_exhaust_gas_pl(pl) * absorption_chiller.COP_heat self.heat = self.heat_absorption_chiller + ice.get_jacket_water_pl(pl) heat_follow = HeatFollow(t, temporary, heat_stor, self.heat) heat_stor = heat_follow.heat_stor self.heat_gas_boiler = heat_follow.gas_boiler_out self.fuel = ice.get_fuel(pl) + gas_boiler.get_Fuel_in(heat_follow.gas_boiler_out) fuel.append(ice.get_fuel(pl) + gas_boiler.get_Fuel_in(heat_follow.gas_boiler_out)) ele_bought.append(ele_drive.ele_bought) ele_ice.append(self.ele) heat_gas_boiler.append(self.heat_gas_boiler) heat_absorption_chiller.append(self.heat_absorption_chiller) if heat_stor > heatstorage.nominal: heat_stor = heatstorage.nominal if ele_stor > elestorage.nominal: ele_stor = elestorage.nominal self.fuel = sum(fuel) self.ele_bought = sum(ele_bought) self.emission_calculate_ice = sum(ele_ice) self.emission_calculate_boiler = sum(heat_gas_boiler) self.emission_calculate_absorption_chiller = sum(self.heat_absorption_chiller) self.emission_calculate_grid = self.ele_bought
def __init__(self, t, temporary, number, ele_stor): ice = InternalCombustionEngine(number, temporary) demand = DemandData() elestorage = EleStorage(temporary) self.elestorage_in = 0 if demand.E[t] <= elestorage.get_E_out_max(ele_stor): self.status = 0 self.elestorage_out = demand.E[t] self.ele = 0 self.ele_bought = 0 elif ( (demand.E[t] > elestorage.get_E_out_max(ele_stor)) & (demand.E[t] <= elestorage.get_E_out_max(ele_stor) + ice.nominal)): self.status = 1 self.elestorage_out = elestorage.get_E_out_max(ele_stor) self.ele = demand.E[t] - self.elestorage_out self.ele_bought = 0 else: self.status = 1 self.elestorage_out = elestorage.get_E_out_max(ele_stor) self.ele = ice.nominal self.ele_bought = demand.E[t] - self.elestorage_out - self.ele self.ele_stor = elestorage.get_S(ele_stor, self.elestorage_in, self.elestorage_out)
def __init__(self, t, temporary, ele_stor, ele, season): # 写这里!!! demand = DemandData() elestorage = EleStorage(temporary) self.ele = ele if season == 0: demand_ele = demand.cold_E elif season == 1: demand_ele = demand.heat_E else: demand_ele = demand.transition_E if demand_ele[t] <= elestorage.get_E_out_max(ele_stor): self.elestorage_out = demand_ele[t] self.elestorage_in = self.ele self.ele_bought = 0 elif ( (demand_ele[t] > elestorage.get_E_out_max(ele_stor)) & (demand_ele[t] <= elestorage.get_E_out_max(ele_stor) + self.ele)): self.elestorage_out = elestorage.get_E_out_max(ele_stor) self.elestorage_in = self.ele + self.elestorage_out - demand_ele[t] self.ele_bought = 0 else: self.elestorage_out = elestorage.get_E_out_max(ele_stor) self.elestorage_in = 0 self.ele_bought = demand_ele[t] - self.ele - self.elestorage_out self.ele_stor = elestorage.get_S(ele_stor, self.elestorage_in, self.elestorage_out)
def __init__(self, temporary, t, ele_stor, season): elestorage = EleStorage(temporary) demand = DemandData() if season == 0: demand_ele = demand.cold_E elif season == 1: demand_ele = demand.heat_E else: demand_ele = demand.transition_E if elestorage.get_E_out_max(ele_stor) > 0: if demand_ele[t] <= elestorage.get_E_out_max(ele_stor): elestorage_ele_out = demand_ele[t] self.ele_bought_for_ele = 0 else: elestorage_ele_out = elestorage.get_E_out_max(ele_stor) self.ele_bought_for_ele = demand_ele[t] - elestorage_ele_out self.ele_stor = elestorage.get_S(ele_stor, 0, elestorage_ele_out) else: self.ele_stor = elestorage.get_S(ele_stor, 0, 0) self.ele_bought_for_ele = demand_ele[t]
def __init__(self, temporary, number): ele_stor = 0 elestorage = EleStorage(temporary) ice = InternalCombustionEngine(number, temporary) demand = DemandData() self.judge = 1 self.ele_bought = [] self.fuel = [] self.ele_ice = [] self.ele_waste = [] for t in range(6, 22, Parameters.delttime): ele = ice.nominal if demand.transition_E[t] <= ele: elestorage_ele_out = 0 elestorage_ele_in = ele - demand.transition_E[t] ele_bought = 0 elif (demand.transition_E[t] > ele) & (demand.transition_E[t] <= ele + elestorage.get_E_out_max(ele_stor)): elestorage_ele_out = demand.transition_E[t] - ele elestorage_ele_in = 0 ele_bought = 0 else: elestorage_ele_out = elestorage.get_E_out_max(ele_stor) elestorage_ele_in = 0 ele_bought = demand.transition_E[t] - ele - elestorage_ele_out ele_stor = elestorage.get_S(ele_stor, elestorage_ele_in, elestorage_ele_out) 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.fuel.append(ice.get_fuel(1)) self.ele_bought.append(ele_bought) for t in range(22, 24, Parameters.delttime): ele_during_shut_down = EleDuringShutDown(temporary, t, ele_stor, 2) ele_stor = ele_during_shut_down.ele_stor self.ele_bought.append(ele_during_shut_down.ele_bought_for_ele) for t in range(0, 6, Parameters.delttime): ele_during_shut_down = EleDuringShutDown(temporary, t, ele_stor, 2) ele_stor = ele_during_shut_down.ele_stor self.ele_bought.append(ele_during_shut_down.ele_bought_for_ele)
def __init__(self, t, temporary, ele_stor, ele): # 写这里!!! demand = DemandData() elestorage = EleStorage(temporary) self.ele = ele if demand.E[t] <= elestorage.get_E_out_max(ele_stor): self.elestorage_out = demand.E[t] self.elestorage_in = self.ele self.ele_bought = 0 elif ((demand.E[t] > elestorage.get_E_out_max(ele_stor)) & (demand.E[t] <= elestorage.get_E_out_max(ele_stor) + self.ele)): self.elestorage_out = elestorage.get_E_out_max(ele_stor) self.elestorage_in = self.ele + self.elestorage_out - demand.E[t] self.ele_bought = 0 else: self.elestorage_out = elestorage.get_E_out_max(ele_stor) self.elestorage_in = 0 self.ele_bought = demand.E[t] - self.ele - self.elestorage_out self.ele_stor = elestorage.get_S(ele_stor, self.elestorage_in, self.elestorage_out)
def __init__(self, temporary, number): ele_stor = 0 elestorage = EleStorage(temporary) ice = InternalCombustionEngine(number, temporary) demand = DemandData() time = demand.E_sheetnrows - 1 ele_bought = [] fuel = [] ele_ice = [] for t in range(0, time, Parameters.delttime): ele_drive = EleDrive(t, temporary, number, ele_stor) self.ele = ele_drive.ele ele_ice.append(self.ele) pl = self.ele / ice.nominal fuel.append(ice.get_fuel(pl)) ele_stor = elestorage.get_S(ele_stor, 0, ele_drive.elestorage_out) ele_bought.append(ele_drive.ele_bought) if ele_stor > elestorage.nominal: ele_stor = elestorage.nominal self.fuel = sum(fuel) self.ele_bought = sum(ele_bought)
def running_judge(temporary, number, season, mode, steam_or_not): heatstor = 0 coldstor = 0 elestor = 0 cold_stor_list = [] heat_stor_list = [] ele_stor_list = [] chp = CHP(temporary, number) gasboiler = GasBoiler(temporary) heatpump = HeatPump(temporary) absorption_chiller = AbsorptionChiller(temporary) boiler = Boiler(temporary) elestorage = EleStorage(temporary) heatstorage = HeatStorage(temporary) coldstorage = ColdStorage(temporary) totalfuel_gasturbine = 0 totalfuel_gasboiler = 0 totalele_powergrid = 0 total_ele_waste = 0 total_heat_waste = 0 total_cold_waste = 0 demand = DemandData() if season == 0: demand_ele = demand.cold_E demand_cold = demand.C if steam_or_not == 0: demand_heat = [] for hour in range(24): demand_heat.append(0) else: demand_heat = demand.cold_Steam elif season == 1: demand_ele = demand.heat_E if steam_or_not == 0: demand_heat = demand.H else: demand_heat = [] for hour in range(24): demand_heat.append(demand.H[hour] + demand.heat_Steam[hour]) demand_cold = [] for hour in range(24): demand_cold.append(0) else: demand_ele = demand.transition_E demand_cold = [] if steam_or_not == 0: demand_heat = [] for hour in range(24): demand_heat.append(0) demand_cold.append(0) else: demand_heat = demand.transition_Steam for hour in range(24): demand_cold.append(0) for t in range(0, demand.E_sheetnrows - 1, Parameters.delttime): cold_stor_list.append(coldstor) heat_stor_list.append(heatstor) ele_stor_list.append(elestor) if mode == 0: # 以冷定热再定电 chp_cold = min(demand_cold[t], chp.C_out_max) boiler_heat_out = absorption_chiller.get_H_in( chp_cold) / Parameters.k chp_heat = boiler_heat_out * (1 - Parameters.k) if demand_heat[t] > heatstorage.get_H_out_max( heatstor) + chp_heat + gasboiler.nominal: return 0 else: result = mode_cold_first(t, temporary, number, coldstor, heatstor, elestor, season, steam_or_not) elif mode == 1: # 以电定热再定冷 暂缺base load运行(运行模式已写好,就是还没加进来) chp_ele = min(demand_ele[t], chp.E_out_max) boiler_heat_in = chp_ele / chp.heat_ele_ratio chp_heat = boiler.get_H_out(boiler_heat_in) * (1 - Parameters.k) chp_cold = chp_heat / (1 - Parameters.k) * Parameters.k if ((demand_heat[t] > heatstorage.get_H_out_max(heatstor) + chp_heat + gasboiler.nominal) | (demand_cold[t] > coldstorage.get_C_out_max(coldstor) + chp_cold + heatpump.nominal)): return 0 else: result = mode_ele_first(t, temporary, number, coldstor, heatstor, elestor, season, steam_or_not) else: if (t >= 7) & (t <= 21): if season == 0: if demand_cold[t] > coldstorage.get_C_out_max( coldstor) + chp.C_out_max + heatpump.nominal: return 0 else: if demand_heat[t] > heatstorage.get_H_out_max( heatstor) + chp.H_out_max + gasboiler.nominal: return 0 else: if season == 0: if demand_cold[t] > coldstorage.get_C_out_max( coldstor) + heatpump.nominal: return 0 else: if demand_heat[t] > heatstorage.get_H_out_max( heatstor) + gasboiler.nominal: return 0 result = mode_base_load(t, temporary, number, coldstor, heatstor, elestor, season, steam_or_not) coldstorage_cold_in = result[0] coldstorage_cold_out = result[1] heatstorage_heat_in = result[2] heatstorage_heat_out = result[3] elestorage_ele_in = result[4] elestorage_ele_out = result[5] powergrid_ele_out = result[11] gasboiler_fuel = result[12] gasturbine_fuel = result[13] totalfuel_gasturbine = totalfuel_gasturbine + gasturbine_fuel totalfuel_gasboiler = totalfuel_gasboiler + gasboiler_fuel totalele_powergrid = totalele_powergrid + powergrid_ele_out elestor = elestorage.get_S(elestor, elestorage_ele_in, elestorage_ele_out) if elestor > elestorage.nominal: total_ele_waste = total_ele_waste + elestor - elestorage.nominal elestor = elestorage.nominal # 若冲进的电过多,则剩余部分浪费 heatstor = heatstorage.get_S(heatstor, heatstorage_heat_in, heatstorage_heat_out) if heatstor > heatstorage.nominal: total_heat_waste = total_heat_waste + heatstor - heatstorage.nominal heatstor = heatstorage.nominal coldstor = coldstorage.get_S(coldstor, coldstorage_cold_in, coldstorage_cold_out) if coldstor > coldstorage.nominal: total_cold_waste = total_cold_waste + coldstor - coldstorage.nominal coldstor = coldstorage.nominal output = (1, totalfuel_gasturbine, totalfuel_gasboiler, totalele_powergrid, total_ele_waste, total_heat_waste, total_cold_waste) return output
def __init__(self, temporary, number): demand = DemandData() heatstorage = HeatStorage(temporary) coldstorage = ColdStorage(temporary) elestorage = EleStorage(temporary) gas_boiler = GasBoiler(temporary) heat_pump = HeatPump(temporary) ice = InternalCombustionEngine(number, temporary) chp = CHPInternalCombustionEngine(temporary, number) heat_stor = 0 ele_stor = 0 cold_stor = 0 self.fuel = [] self.ele_bought = [] self.ele_ice = [] self.heat_jacket_water = [] self.heat_gas_boiler = [] self.cold_absorption_chiller = [] self.cold_heat_pump = [] self.cold_waste = [] self.heat_waste = [] self.ele_waste = [] for t in range(7, 23, Parameters.delttime): if ((demand.C[t] > chp.cold_out_max_exhaust_gas + coldstorage.get_C_out_max(cold_stor) + heat_pump.nominal) | (demand.H[t] > (ice.get_jacket_water_pl(1) + heatstorage.get_H_out_max(heat_stor) + gas_boiler.nominal))): self.judge = 0 break else: self.judge = 1 ele = ice.nominal heat = ice.get_jacket_water_pl(1) cold = chp.cold_out_max_exhaust_gas if demand.C[t] <= cold: coldstorage_cold_in = cold - demand.C[t] coldstorage_cold_out = 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] coldstorage_cold_in = 0 cold_heat_pump = 0 ele_needed_for_heat_pump = 0 else: coldstorage_cold_out = coldstorage.get_C_out_max(cold_stor) coldstorage_cold_in = 0 cold_heat_pump = demand.C[t] - coldstorage_cold_out - cold ele_needed_for_heat_pump = heat_pump.get_E_in( cold_heat_pump) cold_stor = coldstorage.get_S(cold_stor, coldstorage_cold_in, coldstorage_cold_out) if demand.transition_E <= ele: ele_bought_for_ele = 0 if ele_needed_for_heat_pump > 0: if ele - demand.transition_E[ t] >= ele_needed_for_heat_pump: elestorage_ele_in = ele - demand.transition_E[ t] - ele_needed_for_heat_pump elestorage_ele_out = 0 ele_bought_for_heat_pump = 0 else: elestorage_ele_in = 0 elestorage_ele_out = 0 ele_bought_for_heat_pump = ele_needed_for_heat_pump - ( ele - demand.transition_E[t]) else: ele_bought_for_heat_pump = 0 elestorage_ele_out = 0 elestorage_ele_in = ele - demand.transition_E[t] elif(demand.transition_E[t] > ele) \ & (demand.transition_E[t] < ele + elestorage.get_E_out_max(ele_stor)): ele_bought_for_ele = 0 elestorage_ele_out = demand.transition_E[t] - ele elestorage_ele_in = 0 ele_bought_for_heat_pump = ele_needed_for_heat_pump else: ele_bought_for_ele = demand.transition_E[ t] - ele - elestorage.get_E_out_max(ele_stor) elestorage_ele_out = elestorage.get_E_out_max(ele_stor) elestorage_ele_in = 0 ele_bought_for_heat_pump = ele_needed_for_heat_pump ele_stor = elestorage.get_S(ele_stor, elestorage_ele_in, elestorage_ele_out) ele_bought = ele_bought_for_heat_pump + ele_bought_for_ele if demand.H[t] <= heat: heatstorage_heat_out = 0 heatstorage_heat_in = heat - demand.H[t] heat_gas_boiler = 0 elif (demand.H[t] > heat) & (demand.H[t] <= heat + heatstorage.get_H_out_max(heat_stor)): heatstorage_heat_out = demand.H[t] - heat heatstorage_heat_in = 0 heat_gas_boiler = 0 else: heatstorage_heat_out = heatstorage.get_H_out_max(heat_stor) heatstorage_heat_in = 0 heat_gas_boiler = demand.H[ t] - heat - heatstorage.get_H_out_max(heat_stor) heat_stor = heatstorage.get_S(heat_stor, heatstorage_heat_in, heatstorage_heat_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 heat_stor > heatstorage.nominal: self.heat_waste.append(heat_stor - heatstorage.nominal) heat_stor = heatstorage.nominal else: self.heat_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.heat_jacket_water.append(heat) self.cold_absorption_chiller.append(cold) self.cold_heat_pump.append(cold_heat_pump) self.heat_gas_boiler.append(heat_gas_boiler) self.ele_bought.append(ele_bought) self.fuel.append( ice.get_fuel(1) + gas_boiler.get_Fuel_in(heat_gas_boiler)) if self.judge == 1: for t in range(23, 24, Parameters.delttime): cold_during_shut_down = ColdDuringShutDown( temporary, t, cold_stor) heat_during_sut_down = HeatDuringShutDown( temporary, t, heat_stor) ele_during_shut_down = EleDuringShutDown( temporary, t, ele_stor, 2) cold_stor = cold_during_shut_down.cold_stor heat_stor = heat_during_sut_down.heat_stor ele_stor = ele_during_shut_down.ele_stor self.cold_heat_pump.append( cold_during_shut_down.cold_heat_pump) self.heat_gas_boiler.append( heat_during_sut_down.heat_gas_boiler) self.fuel.append(heat_during_sut_down.fuel) self.ele_bought.append( cold_during_shut_down.ele_bought_for_heat_pump + ele_during_shut_down.ele_bought_for_ele) for t in range(0, 7, Parameters.delttime): cold_during_shut_down = ColdDuringShutDown( temporary, t, cold_stor) heat_during_sut_down = HeatDuringShutDown( temporary, t, heat_stor) ele_during_shut_down = EleDuringShutDown( temporary, t, ele_stor, 2) cold_stor = cold_during_shut_down.cold_stor heat_stor = heat_during_sut_down.heat_stor ele_stor = ele_during_shut_down.ele_stor self.cold_heat_pump.append( cold_during_shut_down.cold_heat_pump) self.heat_gas_boiler.append( heat_during_sut_down.heat_gas_boiler) self.fuel.append(heat_during_sut_down.fuel) self.ele_bought.append( cold_during_shut_down.ele_bought_for_heat_pump + ele_during_shut_down.ele_bought_for_ele)
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 = heatstorage.get_H_out_max(heat_stor) heat_gas_boiler = demand.H[ t] - heat - heatstorage.get_H_out_max(heat_stor) heat_stor = heatstorage.get_S(heat_stor, heatstorage_heat_in, heatstorage_heat_out) if demand.heat_E[t] <= ele: elestorage_ele_in = ele - demand.heat_E[t] elestorage_ele_out = 0 ele_bought = 0 elif (demand.heat_E[t] > ele) & (demand.heat_E[t] <= ele + elestorage.get_E_out_max(ele_stor)): elestorage_ele_in = 0 elestorage_ele_out = elestorage.get_E_out_max( ele_stor + ele - demand.heat_E[t]) ele_bought = 0 else: elestorage_ele_in = 0 elestorage_ele_out = elestorage.get_E_out_max(ele_stor) ele_bought = demand.heat_E[ t] - ele - elestorage.get_E_out_max(ele_stor) ele_stor = elestorage.get_S(ele_stor, elestorage_ele_in, elestorage_ele_out) if heat_stor > heatstorage.nominal: self.heat_waste.append(heat_stor - heatstorage.nominal) heat_stor = heatstorage.nominal else: self.heat_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.heat_gas_boiler.append(heat_gas_boiler) self.ele_bought.append(ele_bought) self.fuel.append( ice.get_fuel(1) + gas_boiler.get_Fuel_in(heat_gas_boiler)) if self.judge == 1: for t in range(23, 24): heat_during_shut_down = HeatDuringShutDown( temporary, t, heat_stor) ele_during_shut_down = EleDuringShutDown( temporary, t, ele_stor, 1) heat_stor = heat_during_shut_down.heat_stor ele_stor = ele_during_shut_down.ele_stor self.heat_gas_boiler.append( heat_during_shut_down.heat_gas_boiler) self.fuel.append(heat_during_shut_down.fuel) self.ele_bought.append(ele_during_shut_down.ele_bought_for_ele) for t in range(0, 7, Parameters.delttime): heat_during_shut_down = HeatDuringShutDown( temporary, t, heat_stor) ele_during_shut_down = EleDuringShutDown( temporary, t, ele_stor, 1) heat_stor = heat_during_shut_down.heat_stor ele_stor = ele_during_shut_down.ele_stor self.heat_gas_boiler.append( heat_during_shut_down.heat_gas_boiler) self.fuel.append(heat_during_shut_down.fuel) self.ele_bought.append(ele_during_shut_down.ele_bought_for_ele)
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)
def __init__(self, temporary, number): demand = DemandData() heatstorage = HeatStorage(temporary) coldstorage = ColdStorage(temporary) elestorage = EleStorage(temporary) gas_boiler = GasBoiler(temporary) heat_pump = HeatPump(temporary) ice = InternalCombustionEngine(number, temporary) absorption_chiller = DoubleEffectAbsorptionChiller(temporary) chp = CHPInternalCombustionEngine(temporary, number) time = demand.E_sheetnrows - 1 heat_stor = 0 ele_stor = 0 cold_stor = 0 self.fuel = [] self.ele_bought = [] self.ele_ice = [] self.heat = [] self.heat_gas_boiler = [] self.cold_absorption_chiller = [] self.cold_heat_pump = [] self.cold_waste = [] self.heat_waste = [] self.ele_waste = [] for t in range(0, time, Parameters.delttime): ele = min(demand.transition_E[t] - elestorage.get_E_out_max(ele_stor), ice.nominal) cold_heat = chp.get_cold_and_heat_through_ele(ele) judge_cold = cold_heat[0] judge_heat = cold_heat[1] if ((demand.C[t] > coldstorage.get_C_out_max(cold_stor) + judge_cold + heat_pump.nominal) | (demand.H[t] > heatstorage.get_H_out_max(heat_stor) + judge_heat + gas_boiler.nominal)): self.judge = 0 break else: self.judge = 1 ele_drive = EleDrive(t, temporary, number, ele_stor, 2) if ele_drive == 0: heat_drive = HeatDriveJW(t, temporary, number, heat_stor) gas_boiler_out = heat_drive.gas_boiler_out if heat_drive == 0: cold_drive = ColdDrive(t, temporary, number, cold_stor) cold = cold_drive.cold exhaust_gas = cold / absorption_chiller.COP_double pl = ice.get_pl_through_exhaust_gas(exhaust_gas) ele = pl * ice.nominal heat = ice.get_jacket_water_pl(pl) cold_stor = cold_drive.cold_stor heat_stor = heatstorage.get_S(heat_stor, heat, heat_drive.heatstorage_out) ele_stor = elestorage.get_S(ele_stor, ele, ele_drive.elestorage_out) self.cold_heat_pump.append(cold_drive.heat_pump_out) self.fuel.append(ice.get_fuel(pl)) self.ele_bought.append(heat_pump.get_E_in(cold_drive.heat_pump_out)) else: heat = heat_drive.heat pl = ice.get_pl_through_jacket_water(heat) ele = pl * ice.nominal exhaust_gas = ice.get_exhaust_gas_pl(pl) cold = exhaust_gas * absorption_chiller.COP_double cold_follow = ColdFollow(t, temporary, cold_stor, cold) cold_stor = cold_follow.cold_stor heat_stor = heat_drive.heat_stor ele_stor = elestorage.get_S(ele_stor, ele, ele_drive.elestorage_out) self.cold_heat_pump.append(cold_follow.heat_pump_out) self.fuel.append(ice.get_fuel(pl) + gas_boiler.get_Fuel_in(heat_drive.gas_boiler_out)) self.ele_bought.append(heat_pump.get_E_in(cold_follow.heat_pump_out)) else: ele = ele_drive.ele pl = ele / ice.nominal heat = ice.get_jacket_water_pl(pl) cold = ice.get_exhaust_gas_pl(pl) * absorption_chiller.COP_double heat_follow = HeatFollow(t, temporary, heat_stor, heat) cold_follow = ColdFollow(t, temporary, cold_stor, cold) gas_boiler_out = heat_follow.gas_boiler_out cold_stor = cold_follow.cold_stor heat_stor = heat_follow.heat_stor ele_stor = ele_drive.ele_stor self.cold_heat_pump.append(cold_follow.heat_pump_out) self.fuel.append(ice.get_fuel(pl) + gas_boiler.get_Fuel_in(heat_follow.gas_boiler_out)) self.ele_bought.append(ele_drive.ele_bought + heat_pump.get_E_in(cold_follow.heat_pump_out)) self.ele_ice.append(ele) self.heat.append(heat) self.heat_gas_boiler.append(gas_boiler_out) self.cold_absorption_chiller.append(cold) if cold_stor > coldstorage.nominal: self.cold_waste.append(cold_stor - coldstorage.nominal) cold_stor = coldstorage.nominal else: self.cold_waste.append(0) if heat_stor > heatstorage.nominal: self.heat_waste.append(heat_stor - heatstorage.nominal) heat_stor = heatstorage.nominal else: self.heat_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)
def running_judge(temporary, mode, season): heatstor = 0 coldstor = 0 elestor = 0 cold_stor_list = [] heat_stor_list = [] ele_stor_list = [] chp = CHP(temporary) gasboiler = GasBoiler(temporary) heatpump = HeatPump(temporary) absorption_chiller = AbsorptionChiller(temporary) boiler = Boiler(temporary) elestorage = EleStorage(temporary) heatstorage = HeatStorage(temporary) coldstorage = ColdStorage(temporary) totalfuel_gasturbine = 0 totalfuel_gasboiler = 0 totalele_powergrid = 0 demand = DemandData() if season == 0: demand_ele = demand.cold_E elif season == 1: demand_ele = demand.heat_E else: demand_ele = demand.transition_E for t in range(0, demand.E_sheetnrows - 1, Parameters.delttime): cold_stor_list.append(coldstor) heat_stor_list.append(heatstor) ele_stor_list.append(elestor) if mode == 0: # 以冷定热再定电 chp_cold = min(demand.C[t], chp.C_out_max) boiler_heat_out = absorption_chiller.get_H_in( chp_cold) / Parameters.k chp_heat = boiler_heat_out * (1 - Parameters.k) if demand.H[t] > heatstorage.get_H_out_max( heatstor) + chp_heat + gasboiler.nominal: return 0 else: result = mode_cold_first(t, temporary, coldstor, heatstor, elestor, season) else: # 以电定热再定冷 暂缺base load运行(运行模式已写好,就是还没加进来) chp_ele = min(demand_ele[t], chp.E_out_max) boiler_heat_in = chp_ele / chp.heat_ele_ratio chp_heat = boiler.get_H_out(boiler_heat_in) * (1 - Parameters.k) chp_cold = chp_heat / (1 - Parameters.k) * Parameters.k if ((demand.H[t] > heatstorage.get_H_out_max(heatstor) + chp_heat + gasboiler.nominal) | (demand.C[t] > coldstorage.get_C_out_max(coldstor) + chp_cold + heatpump.nominal)): return 0 else: result = mode_ele_first(t, temporary, coldstor, heatstor, elestor, season) coldstorage_cold_in = result[0] coldstorage_cold_out = result[1] heatstorage_heat_in = result[2] heatstorage_heat_out = result[3] elestorage_ele_in = result[4] elestorage_ele_out = result[5] powergrid_ele_out = result[11] gasboiler_fuel = result[12] gasturbine_fuel = result[13] totalfuel_gasturbine = totalfuel_gasturbine + gasturbine_fuel totalfuel_gasboiler = totalfuel_gasboiler + gasboiler_fuel totalele_powergrid = totalele_powergrid + powergrid_ele_out elestor = elestorage.get_S(elestor, elestorage_ele_in, elestorage_ele_out) if elestor > elestorage.nominal: elestor = elestorage.nominal # 若冲进的电过多,则剩余部分浪费 heatstor = heatstorage.get_S(heatstor, heatstorage_heat_in, heatstorage_heat_out) if heatstor > heatstorage.nominal: heatstor = heatstorage.nominal coldstor = coldstorage.get_S(coldstor, coldstorage_cold_in, coldstorage_cold_out) if coldstor > coldstorage.nominal: coldstor = coldstorage.nominal output = 1, totalfuel_gasturbine, totalfuel_gasboiler, totalele_powergrid return output