Some experiments on a 5 substations test case

To demonstrate the use of the grid2op framework, we propose here to show how to make some quick studies on a purely fictious test case, a 5 bus system.

This system should not be used for deep research purpose. It is provided here as an example and a tool to get used to how some analyses can be performed using grid2op.

First, we can create a environment:

In [1]:
%matplotlib inline
import grid2op
env = grid2op.make("case5_example")

Inspect the powergrid

In grid2op, to make more concrete the powergrid, we also added some visual capabilities to represent, for example, the names of the objects on the grid, or the thermal limit, or any other data you want. All of that is for now only available in the matplotlib "plot helper" and can be accessed as followed.

Plotting the layout of the graph

This utility function allows you to visualize where are the object located on the map. Each object has its own color, and is represented with its name and its id (id starts at 0, following python convention).

In [2]:
from grid2op.Plot import PlotMatplotlib
plot_helper = PlotMatplotlib(env.observation_space)
_ = plot_helper.plot_layout()

Visualizing data on the grid

It is, with the same method, also possible to visualize any data on the powergrid, for example, we can inspect the thermal limit (line property) as follow:

In [3]:
_ = plot_helper.plot_info(line_info=env._thermal_limit_a,
                         colormap="line")

The above plot can be improved. It has lots of useless informations for the purpose (all the names of the loads, generetors and substations). They can be remove simply with:

In [4]:
_ = plot_helper.plot_info(line_info=env._thermal_limit_a,
                          load_info=[None for el in range(env.n_load)],
                          gen_info=[None for el in range(env.n_gen)],
                          sub_info=[None for el in range(env.n_sub)],
                         colormap="line")

Similarly, it is possible to display some data about the generators, for example the maximum value of power they can produce.

In [5]:
_ = plot_helper.plot_info(line_info=[None for el in range(env.n_line)],
                          load_info=[None for el in range(env.n_load)],
                          gen_info=env.gen_pmax,
                          sub_info=[None for el in range(env.n_sub)],
                         colormap="gen")

Create an agent

Now we can see how well the "do nothing" agent (the most basic imaginable) can performe, using some "gym like" methods

In [6]:
from grid2op.Agent import DoNothingAgent
my_agent = DoNothingAgent(env.action_space)

Now do the standard gym loop, and we save all the observations:

In [7]:
all_obs = []
obs = env.reset()
all_obs.append(obs)
reward = env.reward_range[0]
done = False
nb_step = 0
while True:
    action = my_agent.act(obs, reward, done)
    obs, reward, done, _ = env.step(action)
    if done:
        break
    all_obs.append(obs)
    nb_step += 1

We can check if the episode has been completed, or if there has been a game over:

In [8]:
print("Number of timesteps computed: {}".format(nb_step))
print("Total maximum number of timesteps possible: {}".format(env.chronics_handler.max_timestep()))
Number of timesteps computed: 97
Total maximum number of timesteps possible: 2016

As we see here, there is a game over: the agent successfully managed 97 timesteps, while the episode could have lasted 2016.

Let's try to investigate this, for example by plotting the last observation. First we need to create a utilities that will make the plot:

In [9]:
from grid2op.Plot import PlotPlotly
plot_helper = PlotPlotly(observation_space=env.observation_space)
---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
<ipython-input-9-5ddf1aab6591> in <module>
      1 from grid2op.Plot import PlotPlotly
----> 2 plot_helper = PlotPlotly(substation_layout=env.graph_layout,
      3                       observation_space=env.observation_space)

AttributeError: 'Environment' object has no attribute 'graph_layout'
In [10]:
last_obs = all_obs[-1]
fig = plot_helper.plot_observation(last_obs)
fig.show()
---------------------------------------------------------------------------
AttributeError                            Traceback (most recent call last)
<ipython-input-10-d86f2716abe0> in <module>
      1 last_obs = all_obs[-1]
----> 2 fig = plot_helper.plot_observation(last_obs)
      3 fig.show()

AttributeError: 'PlotMatplotlib' object has no attribute 'plot_observation'

As we can see, the last obsevation is pretty clear: 4 powerlines have been disconnected, thus isolating the load on the bottom right. This lead to a game over.

It's also possible, of course, to inspect the previous state, just before this one:

In [ ]:
previous_obs = all_obs[-2]
fig = plot_helper.plot_observation(previous_obs)
fig.show()

And now we can know the cause for this problem: all the powerlines that could provide power on the bottom right load are overloaded in this situation, so the protection worked and disconnect them to prevent hurting the surroundings.