Objectives
This notebook will cover some basic raw functionality at first, and then how they are encapsulated with easy to use functions.
The recommended way to use these is to through the Runner, and not by getting through the instanciation of class one by one.
import os
import sys
import grid2op
res = None
try:
from jyquickhelper import add_notebook_menu
res = add_notebook_menu()
except ModuleNotFoundError:
print("Impossible to automatically add a menu / table of content to this notebook.\nYou can download \"jyquickhelper\" package with: \n\"pip install jyquickhelper\"")
res
Though the Grid2Op
package can be used to perform many different tasks, these set of notebooks will be focused on the machine learning part, and its usage in a Reinforcement learning framework.
The reinforcement learning is a framework that allows to train "agent" to solve time dependant domain. We tried to cast the grid operation planning into this framework. The package Grid2Op
was inspired by it.
In a reinforcement learning (RL), there are 2 distinct entities:
These 2 entities exchange 3 main type of information:
A schematic representaiton of this is shown in the figure bellow (Credit: Sutton & Barto):
In this notebook, we will develop a simple Agent that takes some action (powerline disconnection) based on the observation of the environment.
In order to be initialized, an Agent need to know in which space he operates. For that, we load an Environement, based on the IEEE case14.
An example of this powergrid can be found in the package data. We import them here:
powergrid_path = grid2op.CASE_14_FILE
multi_episode_path = grid2op.CHRONICS_MLUTIEPISODE
names_chronics_to_backend = grid2op.NAMES_CHRONICS_TO_BACKEND
max_iter = 10
*NB* In order to work smoothly, the backend and the data in the files need to have the exact same name. As often the process to generate the data (we suppose that the equilibrium between productions and load is performed beforehand) is different and agnostic from the powergrid, it is not surpising to have the same physical object with different name in the data for the temporal series and for the powergrid file description. In order to simplify the matching, in this example we use the mapper names_chronics_to_backend
that is able to "convert" the names given in the data to the names of the objects in the powergrid description file. It's a mapper that makes the link between the two.
More detail about how the matching is performed can be found in the help of the ChronicsHandler.GridValue.initialize method here or in the file ChronicsHandler.py
In order to work, an Environement need to be fed with data. These data can be read from files for example. Some examples are also provided in this package. We import them:
from grid2op.Chronics import ChronicsHandler, Multifolder, GridStateFromFileWithForecasts
data_feeding = ChronicsHandler(chronicsClass=Multifolder,
path=multi_episode_path,
gridvalueClass=GridStateFromFileWithForecasts,
max_iter=max_iter)
data_fed
is now an instance that will automatically load the data and modify the powergrid accordingly. The process of reading is handled by this class, but the process of modifying the underlying powergrid is carried out by the Environment and performed by the Backend.
A Backend is a dedicated object that has the responsability to compute the resulting powerflow from given injections (poductions and loads). The possibility to implement your own Backend make the Grid2Op framework completely agnostic of the modeling of the powergrid you want to use, or the method to solve the powerflow.
A implementation of a Backend is provided with Pandapower.
from grid2op.Backend import PandaPowerBackend
backend = PandaPowerBackend()
backend
is now a variable that is able to compute powerflow and is able to emulate cascading failures. To be able to work properly and carry out the right computations, the worker need to be aware of some Parameters.
For this example, we will use the default parameters available. More information about the parameters that can be modified can be found in the help here, or in the file Parameters.py.
from grid2op.Parameters import Parameters
param = Parameters()
from grid2op.Environment import Environment
env = Environment(init_grid_path=powergrid_path,
chronics_handler=data_feeding,
backend=backend,
parameters=param,
names_chronics_to_backend=names_chronics_to_backend)
Creating an Environment will load the powergrid, load the data to feed it. Use the first row to initialize the powergrid and some other internal checking that the data are suited to the powergrid.
It can take a moment (usually a few seconds) to load it.
This environment can be greatly customize. We expose here only basic functionnalities. For more information it is advised to read the documentation here (if it has been built locally), to consult the official documentation on the internet or to consult the source code of Environment.py
All of the above can be done in calling one function that will handle the Environment create.
This environment is named case14_fromfile
.
To definie it we can call:
env = grid2op.make("case14_fromfile")
An Agent is the name given to the "operator" / "bot" / "algorithm" that will perform some modification of the powergrid when he faces some "observation".
Some example of Agents are provided in the file Agent.py.
A deeper look at the different Agent provided can be found in the 4_StudyYourAgent notebook (in progress). We suppose here we use the most simple Agent, the one that does nothing
from grid2op.Agent import DoNothingAgent
my_agent = DoNothingAgent(env.helper_action_player)
The performance of each Agent is assessed with the reward. For this example, the reward is a FlatReward that just computes how many times step the Agent has managed to performed before breaking any rules. For more control on this reward, it is recommended to use look at the document of the Environment class.
More example of rewards are also available at Rewards.py or on the official document or here (if the documentation has been built locally -- recommended).
done = False
time_step = int(0)
cum_reward = 0.
act = env.helper_action_player({}) # initialize the environement with a "do nothing action"
while not done:
obs, reward, done, info = env.step(act) # implement this action on the powergrid
act = my_agent.act(obs, reward, done) # chose an action to do, in this case "do nothing"
cum_reward += reward
time_step += 1
if time_step > max_iter:
break
We can now evaluate how well this agent is performing:
print("This agent managed to survive {} timesteps".format(time_step))
print("It's final cumulated reward is {}".format(cum_reward))
This agent managed to survive 11 timesteps It's final cumulated reward is 219.99783718276583
All the above steps have been detailed as a "quick start", to give an example of the main classes of the Grid2Op package. Having to code all the above is quite tedious, but offers a lot of flexibility.
Implementing all this before starting to evaluate an agent can be quite tedious. What we expose here is a much shorter way to perfom all of the above. In this section we will expose 2 ways:
For this section, we assume the same as before:
When only simple assessment need to be performed, the grid2op.main API is perfectly suited. This API can also be access with the command line:
python3 -m grid2op.main
We detail here its usage as an API, to assess the performance of a given Agent.
As opposed to building en environment from scratch (see the previous section) this requires much less effort: we don't need to initialize (instanciate) anything. All is carried out inside the Runner called by the main function.
We ask here 1 episode (eg. we play one scenario until: either the agent does a game over, or until the scenario ends). But this method would work as well if we more.
from grid2op.main import main
res = main(nb_episode=1,
agent_class=DoNothingAgent,
path_casefile=powergrid_path,
path_chronics=multi_episode_path,
names_chronics_to_backend=names_chronics_to_backend,
gridStateclass_kwargs={"gridvalueClass": GridStateFromFileWithForecasts, "max_iter": max_iter}
)
A call of the single 2 lines above will:
print("The results are:")
for chron_name, _, cum_reward, nb_time_step, max_ts in res:
msg_tmp = "\tFor chronics located at {}\n".format(chron_name)
msg_tmp += "\t\t - cumulative reward: {:.2f}\n".format(cum_reward)
msg_tmp += "\t\t - number of time steps completed: {:.0f} / {:.0f}".format(nb_time_step, max_ts)
print(msg_tmp)
The results are: For chronics located at /home/donnotben/Documents/Grid2Op_dev/getting_started/grid2op/data/test_multi_chronics/1 - cumulative reward: 10.00 - number of time steps completed: 10 / 10
This is particularly suited to evaluate different agents, for example we can quickly evaluate a second agent. This Agent will simulate the effect of disconnecting each powerline on the powergrid and take the best action found (its execution can take a long time, depending on the scenario and the amount of powerlines on the grid). The execution of the code below can take a few minutes
from grid2op.Agent import PowerLineSwitch
res = main(nb_episode=1,
agent_class=PowerLineSwitch,
path_casefile=powergrid_path,
path_chronics=multi_episode_path,
names_chronics_to_backend=names_chronics_to_backend,
gridStateclass_kwargs={"gridvalueClass": GridStateFromFileWithForecasts, "max_iter": max_iter}
)
print("The results are:")
for chron_name, _, cum_reward, nb_time_step, max_ts in res:
msg_tmp = "\tFor chronics located at {}\n".format(chron_name)
msg_tmp += "\t\t - cumulative reward: {:.2f}\n".format(cum_reward)
msg_tmp += "\t\t - number of time steps completed: {:.0f} / {:.0f}".format(nb_time_step, max_ts)
print(msg_tmp)
The results are: For chronics located at /home/donnotben/Documents/Grid2Op_dev/getting_started/grid2op/data/test_multi_chronics/1 - cumulative reward: 10.00 - number of time steps completed: 10 / 10
Using this API it's also possible to store the results for a detailed examination of the aciton taken by the Agent. Note that writing on the hard drive has an overhead on the computation time.
To do this, only a simple argument need to be added. An example can be found below (where the outcome of the experiment will be stored in the saved_experiment_donothing
directory):
res = main(nb_episode=1,
agent_class=DoNothingAgent,
path_casefile=powergrid_path,
path_chronics=multi_episode_path,
names_chronics_to_backend=names_chronics_to_backend,
gridStateclass_kwargs={"gridvalueClass": GridStateFromFileWithForecasts, "max_iter": max_iter},
path_save=os.path.abspath("saved_experiment_donothing")
)
print("The results are:")
for chron_name, _, cum_reward, nb_time_step, max_ts in res:
msg_tmp = "\tFor chronics located at {}\n".format(chron_name)
msg_tmp += "\t\t - cumulative reward: {:.2f}\n".format(cum_reward)
msg_tmp += "\t\t - number of time steps completed: {:.0f} / {:.0f}".format(nb_time_step, max_ts)
print(msg_tmp)
The results are: For chronics located at /home/donnotben/Documents/Grid2Op_dev/getting_started/grid2op/data/test_multi_chronics/1 - cumulative reward: 10.00 - number of time steps completed: 10 / 10
!ls saved_experiment_donothing/1
actions.npy episode_meta.json _parameters.json agent_exec_times.npy episode_times.json rewards.npy disc_lines_cascading_failure.npy observations.npy env_modifications.npy other_rewards.json
All the informations saved are showed above. For more information about them, please don't hesitate to read the documentation of Runner (if compiled locally) or consult the Runner.py file.
NB A lot more of informations about Action is provided in the 2_Action_GridManipulation notebook. In the 3_TrainingAnAgent (last section) there is an quick example on how to read / write action from a saved repository.
By default, the grid2op framework come with a few pre defined environment, each with its own properties. Some are suitable for discrete control, some more for continuous control etc.
Two environments can be used to get familiar with the platform. The first one is the "case5_example" that represents a tiny powergrid (with only 5 substations), the other one is the "case14_redisp" that is a transposition of the case14 ieee powergrid TODO describe it.
They can both be created easily:
env_case5 = grid2op.make("case5_example")
env_case14 = grid2op.make("case14_redisp")
# more data for this environment are available as a github release . You can download them easily with:
# python -m grid2op.download --name "case14_redisp" --path_save PATH\WHERE\YOU\WANT\TO\DOWNLOAD\DATA
# from a command line
/home/donnotben/Documents/Grid2Op_dev/getting_started/grid2op/MakeEnv.py:652: UserWarning: Your are using only 2 chronics for this environment. More can be download by running, from a command line: python -m grid2op.download --name "case14_redisp" --path_save PATH\WHERE\YOU\WANT\TO\DOWNLOAD\DATA
from grid2op.Runner import Runner
from grid2op.Agent import DoNothingAgent
runner = Runner(**env_case14.get_params_for_runner(), agentClass=DoNothingAgent)
res = runner.run(nb_episode=2, nb_process=1, max_iter=10)
for path_chron, episode_id, total_reward, nb_iter, max_iter in res:
print("Total reward for episode {} is {}".format(episode_id, total_reward))
Total reward for episode 0 is 10948.052798806984 Total reward for episode 1 is 11096.619599771308