def testRandomSawToothStorage(self):
lp = self.lp
solar = self.solar
lp.add_nondispatchable_sources(solar)
storage = GridStorage(STORAGE, 0)
lp.add_storage(storage)
for storage_efficiency in np.linspace(0.1, 1.0, 4):
storage.storage_efficiency = storage_efficiency
self.assertTrue(lp.solve())
# build up storage profile
demand_profile = self.demand_profile
golden_storage_profile = []
last_storage_value = 0.0
for d in reversed(demand_profile):
next_storage_value = (last_storage_value + d) / storage_efficiency
golden_storage_profile.append(next_storage_value)
last_storage_value = next_storage_value
golden_storage_profile.reverse()
# First generation comes from solar, so set golden_storage[0] = 0.
golden_storage_profile[0] = 0
golden_solar_profile = self.impulse * golden_storage_profile[1]
npt.assert_allclose(solar.get_solution_values(),
golden_solar_profile)
npt.assert_allclose(storage.get_solution_values(),
golden_storage_profile,
atol=1e-7)
def testFreeStorageMeansCheapestSource(self):
"""Verify that free storage selects the cheapest energy supply."""
solar = GridSource(SOLAR, 2.0e6, 0)
wind = GridSource(WIND, 2.2e6, 0)
storage = GridStorage(STORAGE, 0)
lp = self.lp
lp.add_nondispatchable_sources(solar, wind)
lp.add_storage(storage)
self.assertTrue(lp.solve())
npt.assert_almost_equal(wind.get_solution_values(),
np.zeros(2))
npt.assert_almost_equal(solar.get_solution_values(),
np.array([0.0, 2.0]))
npt.assert_almost_equal(storage.get_solution_values(),
np.array([1.0, 0.0]))
npt.assert_almost_equal(storage.source.get_solution_values(),
np.array([1.0, 0.0]))
npt.assert_almost_equal(storage.sink.get_solution_values(),
np.array([0.0, 1.0]))
def testSimpleStorageWithStorageEfficiency(self):
"""Free Storage should backup Solar."""
solar = GridSource(SOLAR, 2.0e6, 0)
wind = GridSource(WIND, 5.0e6, 0)
ng = GridSource(NG, 1.0e10, 0)
se = math.pow(0.5, 1.0 / 3)
storage = GridStorage(STORAGE, 0, storage_efficiency=se)
lp = self.lp
lp.add_nondispatchable_sources(solar, wind)
lp.add_dispatchable_sources(ng)
lp.add_storage(storage)
self.assertTrue(lp.solve())
npt.assert_almost_equal(solar.get_solution_values(),
np.array([4.5, 0, 0, 2.25]))
npt.assert_almost_equal(wind.get_solution_values(),
np.zeros(4))
npt.assert_almost_equal(ng.get_solution_values(),
np.zeros(4))
npt.assert_almost_equal(storage.get_solution_values(),
np.array([0, 1.5, 1.5 * se, 1.5 * se * se]))
npt.assert_almost_equal(storage.source.get_solution_values(),
np.array([0, 0, 0, 0.75]))
npt.assert_almost_equal(storage.sink.get_solution_values(),
np.array([1.5, 0, 0, 0]))
def testSimpleStorageWithSourceLimit(self):
"""Tests Solar, Wind, Storage. Solar limited so wind makes up the rest."""
solar = GridSource(SOLAR, 2.0e6, 0, max_energy=3)
wind = GridSource(WIND, 5.0e6, 0)
ng = GridSource(NG, 1.0e10, 0)
storage = GridStorage(STORAGE, 0)
lp = self.lp
lp.add_nondispatchable_sources(solar, wind)
lp.add_dispatchable_sources(ng)
lp.add_storage(storage)
self.assertTrue(lp.solve())
npt.assert_almost_equal(solar.get_solution_values(),
np.array([2.0, 0, 0, 1.0]))
npt.assert_almost_equal(wind.get_solution_values(),
np.array([3.0, 0, 0, 0.0]))
npt.assert_almost_equal(ng.get_solution_values(),
np.zeros(4))
npt.assert_almost_equal(storage.get_solution_values(),
np.array([0, 2.0, 2.0, 2.0]))
npt.assert_almost_equal(storage.source.get_solution_values(),
np.array([0, 0, 0, 2.0]))
npt.assert_almost_equal(storage.sink.get_solution_values(),
np.array([2.0, 0, 0, 0]))
def testSimpleStorageWithStorageLimit(self):
"""Free Storage should backup Solar."""
solar = GridSource(SOLAR, 2.0e6, 0)
wind = GridSource(WIND, 5.0e6, 0)
ng = GridSource(NG, 1.0e10, 0)
storage = GridStorage(STORAGE, 0, max_storage=0.5)
lp = self.lp
lp.add_nondispatchable_sources(solar, wind)
lp.add_dispatchable_sources(ng)
lp.add_storage(storage)
self.assertTrue(lp.solve())
npt.assert_almost_equal(solar.get_solution_values(),
np.array([5.0, 0, 0, 2.5]))
npt.assert_almost_equal(wind.get_solution_values(),
np.zeros(4))
npt.assert_almost_equal(ng.get_solution_values(),
np.zeros(4))
npt.assert_almost_equal(storage.get_solution_values(),
np.array([0, 0.5, 0.5, 0.5]))
npt.assert_almost_equal(storage.source.get_solution_values(),
np.array([0, 0, 0, 0.5]))
npt.assert_almost_equal(storage.sink.get_solution_values(),
np.array([0.5, 0, 0, 0]))
def testStorageNameplateCost(self):
"""Keep increasing storage costs until ng finally wins out."""
wind = GridSource(WIND, 1.0e6, 0)
storage = GridStorage(STORAGE, 0)
ng = GridSource(NG, 4.6e6, 0)
lp = self.lp
lp.add_nondispatchable_sources(wind)
lp.add_dispatchable_sources(ng)
lp.add_storage(storage)
self.assertTrue(lp.solve())
npt.assert_almost_equal(wind.get_solution_values(),
np.array([2.0, 0.0]))
npt.assert_almost_equal(ng.get_solution_values(),
np.zeros(2))
# Change costs. Total wind + storage cost is now (2 + 1)E6.
# Still less than ng costs.
storage.charge_nameplate_cost = 1.0e6
self.assertTrue(lp.solve())
npt.assert_almost_equal(wind.get_solution_values(),
np.array([2.0, 0.0]))
npt.assert_almost_equal(ng.get_solution_values(),
np.zeros(2))
# Change Costs. Total wind + storage cost is now (2 + 1 + 1)E6.
# Still less than ng costs.
storage.storage_nameplate_cost = 1.0e6
self.assertTrue(lp.solve())
npt.assert_almost_equal(wind.get_solution_values(),
np.array([2.0, 0.0]))
npt.assert_almost_equal(ng.get_solution_values(),
np.zeros(2))
# Change costs. Total wind + storage cost is now
# (2 + 1 + 1 + 1)E6. Now more than ng costs.
storage.discharge_nameplate_cost = 1.0e6
self.assertTrue(lp.solve())
npt.assert_almost_equal(wind.get_solution_values(),
np.zeros(2))
npt.assert_almost_equal(ng.get_solution_values(),
np.array([1.0, 1.0]))
def testCircularStorageFirstHourSource(self):
"""Verify that storage from first hour affects last hour."""
wind = GridSource(WIND, 2.0e6, 0)
storage = GridStorage(STORAGE, 0)
lp = self.lp
lp.add_nondispatchable_sources(wind)
lp.add_storage(storage)
self.assertTrue(lp.solve())
npt.assert_almost_equal(wind.get_solution_values(),
np.array([2.0, 0.0]))
npt.assert_almost_equal(storage.get_solution_values(),
np.array([0.0, 1.0]))
npt.assert_almost_equal(storage.source.get_solution_values(),
np.array([0.0, 1.0]))
npt.assert_almost_equal(storage.sink.get_solution_values(),
np.array([1.0, 0.0]))
def testPostProcessing(self):
lp = self.lp
wind = GridSource(WIND, 2.0e6, 0, is_rps_source=True)
ng = GridSource(NG, 1.0e6, 1.0e6)
storage = GridStorage(STORAGE, 1, 1, 1)
lp.add_nondispatchable_sources(wind)
lp.add_dispatchable_sources(ng)
lp.add_storage(storage)
lp.rps_percent = 20
self.assertTrue(lp.solve())
wind_solution = MockPostProcessingGridSourceSolver(wind)
storage_sink_solution = MockPostProcessingGridSourceSolver(storage.sink)
wind_solution.solution_values -= 0.1
with self.assertRaises(DemandNotSatisfiedError):
lp._post_process()
# After reset it shouldn't error out.
wind_solution.reset()
try:
lp._post_process()
except RuntimeError:
self.fail()
storage_sink_solution.solution_values += 20
with self.assertRaises(DemandNotSatisfiedError):
lp._post_process()
# After reset it shouldn't error out.
storage_sink_solution.reset()
try:
lp._post_process()
except RuntimeError:
self.fail()
lp.rps_demand *= 100
with self.assertRaises(RpsPercentNotMetError):
lp._post_process()