class Simulation:
def __init__(self):
self.b = Battery(constants.battery_capacity,constants.battery_max_charge,constants.battery_max_discharge,constants.nbatt,constants.nbatt_c,constants.nbatt_d,\
constants.battery_cost, constants.life_time, constants.round_trip)
self.g = GasTurbine(constants.gas_turbine_max, constants.microgas_turbine_om_cost, constants.fual_cost, constants.co2_coe, constants.co2_cost, \
constants.so2_coe, constants.so2_cost, constants.no_coe, constants.no_cost)
self.l = Load(constants.shortage_cost)
self.p = PV(constants.pv_max, constants.pv_om_cost)
self.w = WindTurbine(constants.wind_turbine_max,
constants.wind_turbine_om_cost)
self.m = MEMS(self.b, self.g, self.l, self.p, self.w)
self.bm = BaseMEMS(self.b, self.g, self.l, self.p, self.w)
self.l.set_forecast(constants.load_important_forecast,
constants.load_transferable_forecast)
self.p.set_forecast([ir / 0.2 * 1000 for ir in constants.pv_forecast])
self.w.set_forecast(
[0.2 * wind_speed**3 for wind_speed in constants.wind_forecast])
def simulate(self, battery_powers, gas_turbine_powers):
self.b.set_powers(battery_powers)
self.g.set_powers(gas_turbine_powers)
for t in range(24): # 24 hours
self.m.control(t)
def simulate_base(self):
for t in range(24): # 24 hours
self.bm.control(t)
# cost
def get_objective_1(self):
total_cost = 0
for t in range(24):
# maintenance cost
om_cost = self.w.maintenance_factor*self.w.get_power(t)+self.p.maintenance_factor*self.p.get_power(t)\
+self.g.maintenance_factor*self.g.get_power(t)
fual = self.g.fual_cost * self.g.get_power(t)
total_cost = total_cost + om_cost + fual
return total_cost / 1000
# pollution
def get_objective_2(self):
pollution_cost = 0
for t in range(24):
hourly_cost = self.g.get_power(t)*self.g.co2_coe*self.g.co2_cost \
+self.g.get_power(t)*self.g.so2_coe*self.g.so2_cost \
+self.g.get_power(t)*self.g.no_coe*self.g.no_cost
pollution_cost += hourly_cost
return pollution_cost
def get_objective_3(self):
total_discharge = 0
battery_powers = self.b.powers
for t in range(24):
if battery_powers[t] < 0: # charging
total_discharge -= battery_powers[t]
return total_discharge / 10
def plot_power(self):
plt.plot(self.p.powers, label='Solar')
plt.plot(self.w.powers, label='Wind Turbine')
plt.plot(self.b.powers, label='Battery')
plt.plot(self.g.powers, label='Microgas Turbine')
plt.plot([
load1 + load2
for load1, load2 in zip(self.l.important, self.l.transferable)
],
label='Load')
plt.legend()
plt.xlabel("Time (h)")
plt.ylabel("Power (kW)")
plt.show()
