.idea/DRL-for-Energy-Systems.iml

@@ -2,7 +2,7 @@
 <module type="PYTHON_MODULE" version="4">
   <component name="NewModuleRootManager">
     <content url="file://$MODULE_DIR$" />
-    <orderEntry type="jdk" jdkName="Remote Python 3.9.18 (sftp://chenxd@124.16.151.196:22121/home/chenxd/miniconda3/envs/grid/bin/python3.9)" jdkType="Python SDK" />
+    <orderEntry type="jdk" jdkName="rl-microgrid" jdkType="Python SDK" />
     <orderEntry type="sourceFolder" forTests="false" />
   </component>
   <component name="PyDocumentationSettings">

.idea/deployment.xml

@@ -1,22 +0,0 @@
-<?xml version="1.0" encoding="UTF-8"?>
-<project version="4">
-  <component name="PublishConfigData" autoUpload="Always" serverName="chenxd@124.16.151.196:22121 password" remoteFilesAllowedToDisappearOnAutoupload="false">
-    <serverData>
-      <paths name="chenxd@124.16.151.196:22121 password">
-        <serverdata>
-          <mappings>
-            <mapping deploy="/tmp/pycharm_project_78" local="$PROJECT_DIR$" />
-          </mappings>
-        </serverdata>
-      </paths>
-      <paths name="root@124.16.151.196:10531 password">
-        <serverdata>
-          <mappings>
-            <mapping deploy="/mnt/chenxd/DRL-for-Energy-Systems" local="$PROJECT_DIR$" />
-          </mappings>
-        </serverdata>
-      </paths>
-    </serverData>
-    <option name="myAutoUpload" value="ALWAYS" />
-  </component>
-</project>

.idea/misc.xml

@@ -3,5 +3,5 @@
   <component name="Black">
     <option name="sdkName" value="rl-microgrid" />
   </component>
-  <component name="ProjectRootManager" version="2" project-jdk-name="Remote Python 3.9.18 (sftp://chenxd@124.16.151.196:22121/home/chenxd/miniconda3/envs/grid/bin/python3.9)" project-jdk-type="Python SDK" />
+  <component name="ProjectRootManager" version="2" project-jdk-name="rl-microgrid" project-jdk-type="Python SDK" />
 </project>

DDPG.py

@@ -105,7 +105,7 @@ if __name__ == '__main__':
         record = test_one_episode(env, agent.act, agent.device)
         eval_data = pd.DataFrame(record['system_info'])
         eval_data.columns = ['time_step', 'price', 'netload', 'action', 'real_action', 'soc', 'battery', 'gen1', 'gen2',
-                             'gen3', 'temperature', 'irradiance', 'unbalance', 'operation_cost']
+                             'gen3', 'unbalance', 'operation_cost']
     if args.save_test_data:
         test_data_save_path = f'{args.cwd}/test_data.pkl'
         with open(test_data_save_path, 'wb') as tf:

PPO.py

@@ -12,34 +12,33 @@ os.environ['KMP_DUPLICATE_LIB_OK'] = 'TRUE'
 script_name = os.path.basename(__file__)
 
 
+# after adding layer normalization, it doesn't work
 class ActorPPO(nn.Module):
     def __init__(self, mid_dim, state_dim, action_dim, layer_norm=False):
         super().__init__()
-        self.layer_norm = layer_norm
         self.net = nn.Sequential(nn.Linear(state_dim, mid_dim), nn.ReLU(),
                                  nn.Linear(mid_dim, mid_dim), nn.ReLU(),
                                  nn.Linear(mid_dim, mid_dim), nn.Hardswish(),
-                                 nn.Linear(mid_dim, action_dim))
+                                 nn.Linear(mid_dim, action_dim), )
+
         # the logarithm (log) of standard deviation (std) of action, it is a trainable parameter
         self.a_logstd = nn.Parameter(torch.zeros((1, action_dim)) - 0.5, requires_grad=True)
         self.sqrt_2pi_log = np.log(np.sqrt(2 * np.pi))
+        if layer_norm:
+            self.layer_norm(self.net)
 
-        if self.layer_norm:
-            self.apply_layer_norm()
-
-    def apply_layer_norm(self):
-        def init_weights(layer):
-            if isinstance(layer, nn.Linear):
-                nn.init.orthogonal_(layer.weight, 1.0)
-                nn.init.constant_(layer.bias, 0.0)
-
-        self.net.apply(init_weights)
+    @staticmethod
+    def layer_norm(layer, std=1.0, bias_const=0.0):
+        for i in layer:
+            if hasattr(i, 'weight'):
+                torch.nn.init.orthogonal_(i.weight, std)
+                torch.nn.init.constant_(i.bias, bias_const)
 
     def forward(self, state):
-        return self.net(state).tanh()  # action.tanh() limit the data output of action
+        return self.net(state).tanh()  # action.tanh() # in this way limit the data output of action
 
     def get_action(self, state):
-        a_avg = self.forward(state)  # too big for the action
+        a_avg = self.net(state)  # too big for the action
         a_std = self.a_logstd.exp()
 
         noise = torch.randn_like(a_avg)
@@ -47,7 +46,7 @@ class ActorPPO(nn.Module):
         return action, noise
 
     def get_logprob_entropy(self, state, action):
-        a_avg = self.forward(state)
+        a_avg = self.net(state)
         a_std = self.a_logstd.exp()
 
         delta = ((a_avg - action) / a_std).pow(2) * 0.5
@@ -64,21 +63,19 @@ class ActorPPO(nn.Module):
 class CriticAdv(nn.Module):
     def __init__(self, mid_dim, state_dim, _action_dim, layer_norm=False):
         super().__init__()
-        self.layer_norm = layer_norm
         self.net = nn.Sequential(nn.Linear(state_dim, mid_dim), nn.ReLU(),
                                  nn.Linear(mid_dim, mid_dim), nn.ReLU(),
                                  nn.Linear(mid_dim, mid_dim), nn.Hardswish(),
                                  nn.Linear(mid_dim, 1))
-        if self.layer_norm:
-            self.apply_layer_norm()
+        if layer_norm:
+            self.layer_norm(self.net, std=1.0)
 
-    def apply_layer_norm(self):
-        def init_weights(layer):
-            if isinstance(layer, nn.Linear):
-                nn.init.orthogonal_(layer.weight, 1.0)
-                nn.init.constant_(layer.bias, 0.0)
-
-        self.net.apply(init_weights)
+    @staticmethod
+    def layer_norm(layer, std=1.0, bias_const=0.0):
+        for i in layer:
+            if hasattr(i, 'weight'):
+                torch.nn.init.orthogonal_(i.weight, std)
+                torch.nn.init.constant_(i.bias, bias_const)
 
     def forward(self, state):
         return self.net(state)  # Advantage value
@@ -119,6 +116,7 @@ class AgentPPO:
 
         self.cri_optim = torch.optim.Adam(self.cri.parameters(), learning_rate)
         self.act_optim = torch.optim.Adam(self.act.parameters(), learning_rate) if self.ClassAct else self.cri
+        del self.ClassCri, self.ClassAct  # why del self.ClassCri and self.ClassAct here, to save memory?
 
     def select_action(self, state):
         states = torch.as_tensor((state,), dtype=torch.float32, device=self.device)
@@ -131,8 +129,8 @@ class AgentPPO:
         last_done = 0
         for i in range(target_step):
             action, noise = self.select_action(state)
-            # the step of cut action is finally organized into the environment
-            state, next_state, reward, done, = env.step(np.tanh(action))
+            state, next_state, reward, done, = env.step(
+                np.tanh(action))  # here the step of cut action is finally organized into the environment.
             trajectory_temp.append((state, reward, done, action, noise))
             if done:
                 state = env.reset()
@@ -142,8 +140,8 @@ class AgentPPO:
         self.state = state
 
         '''splice list'''
-        # store 0 trajectory information to list
-        trajectory_list = self.trajectory_list + trajectory_temp[:last_done + 1]
+        trajectory_list = self.trajectory_list + trajectory_temp[
+                                                 :last_done + 1]  # store 0 trajectory information to the list
         self.trajectory_list = trajectory_temp[last_done:]
         return trajectory_list
 
@@ -151,12 +149,12 @@ class AgentPPO:
         """put data extract and update network together"""
         with torch.no_grad():
             buf_len = buffer[0].shape[0]
-            # decompose buffer data
-            buf_state, buf_action, buf_noise, buf_reward, buf_mask = [ten.to(self.device) for ten in buffer]
+            buf_state, buf_action, buf_noise, buf_reward, buf_mask = [ten.to(self.device) for ten in
+                                                                      buffer]  # decompose buffer data
 
             '''get buf_r_sum, buf_logprob'''
             bs = 4096  # set a smaller 'BatchSize' when out of GPU memory: 1024, could change to 4096
-            buf_value = [self.cri_target(buf_state[i:i + bs]) for i in range(0, buf_len, bs)]
+            buf_value = [self.cri_target(buf_state[i:i + bs]) for i in range(0, buf_len, bs)]  #
             buf_value = torch.cat(buf_value, dim=0)
             buf_logprob = self.act.get_old_logprob(buf_action, buf_noise)
 
@@ -319,7 +317,7 @@ if __name__ == '__main__':
         agent = args.agent
         env = args.env
         agent.init(args.net_dim, env.state_space.shape[0], env.action_space.shape[0], args.learning_rate,
-                   args.if_per_or_gae, layer_norm=True)
+                   args.if_per_or_gae)
 
         cwd = args.cwd
         gamma = args.gamma
@@ -373,7 +371,7 @@ if __name__ == '__main__':
         record = test_one_episode(env, agent.act, agent.device)
         eval_data = pd.DataFrame(record['system_info'])
         eval_data.columns = ['time_step', 'price', 'netload', 'action', 'real_action', 'soc', 'battery', 'gen1', 'gen2',
-                             'gen3', 'temperature', 'irradiance', 'unbalance', 'operation_cost']
+                             'gen3', 'unbalance', 'operation_cost']
     if args.save_test_data:
         test_data_save_path = f'{args.cwd}/test_data.pkl'
         with open(test_data_save_path, 'wb') as tf:

SAC.py

@@ -109,7 +109,7 @@ if __name__ == '__main__':
         record = test_one_episode(env, agent.act, agent.device)
         eval_data = pd.DataFrame(record['system_info'])
         eval_data.columns = ['time_step', 'price', 'netload', 'action', 'real_action', 'soc', 'battery', 'gen1', 'gen2',
-                             'gen3', 'temperature', 'irradiance', 'unbalance', 'operation_cost']
+                             'gen3', 'unbalance', 'operation_cost']
     if args.save_test_data:
         test_data_save_path = f'{args.cwd}/test_data.pkl'
         with open(test_data_save_path, 'wb') as tf:

TD3.py

@@ -105,7 +105,7 @@ if __name__ == '__main__':
         record = test_one_episode(env, agent.act, agent.device)
         eval_data = pd.DataFrame(record['system_info'])
         eval_data.columns = ['time_step', 'price', 'netload', 'action', 'real_action', 'soc', 'battery', 'gen1', 'gen2',
-                             'gen3', 'temperature', 'irradiance', 'unbalance', 'operation_cost']
+                             'gen3', 'unbalance', 'operation_cost']
     if args.save_test_data:
         test_data_save_path = f'{args.cwd}/test_data.pkl'
         with open(test_data_save_path, 'wb') as tf:

environment.py

@@ -175,7 +175,7 @@ class ESSEnv(gym.Env):
 
         # process_elements(pv, lambda x: x, self.data_manager.add_pv_element)
         process_elements(price, lambda x: max(x / 10, 0.5), self.data_manager.add_price_element)
-        process_elements(load, lambda x: x * 3, self.data_manager.add_load_element)
+        process_elements(load, lambda x: x * 5, self.data_manager.add_load_element)
         process_elements(irradiance, lambda x: x, self.data_manager.add_irradiance_element)
         process_elements(temperature, lambda x: x - 273.15, self.data_manager.add_temperature_element)
         process_elements(wind, lambda x: x, self.data_manager.add_wind_element)

tools.py

@@ -46,37 +46,39 @@ def optimization_base_result(env, month, day, initial_soc):
 
     m = gp.Model("UC")
 
-    # 设置系统变量
+    # set system variables
     on_off = m.addVars(NUM_GEN, period, vtype=GRB.BINARY, name='on_off')
     gen_output = m.addVars(NUM_GEN, period, vtype=GRB.CONTINUOUS, name='output')
-    # 设置充放电约束
+    # set constrains for charge/discharge
     battery_energy_change = m.addVars(period, vtype=GRB.CONTINUOUS, lb=env.battery.max_discharge,
                                       ub=env.battery.max_charge, name='battery_action')
-    # 设置外部电网与能源系统交换约束
+    # set constrains for exchange between external grid and distributed energy system
     grid_energy_import = m.addVars(period, vtype=GRB.CONTINUOUS, lb=0, ub=env.grid.exchange_ability, name='import')
     grid_energy_export = m.addVars(period, vtype=GRB.CONTINUOUS, lb=0, ub=env.grid.exchange_ability, name='export')
     soc = m.addVars(period, vtype=GRB.CONTINUOUS, lb=0.2, ub=0.8, name='SOC')
 
-    # 1. 添加平衡约束
+    # 1. add balance constrain
     m.addConstrs(((sum(gen_output[g, t] for g in range(NUM_GEN)) + pv[t] + wind[t] + grid_energy_import[t] >= load[t] +
                    battery_energy_change[t] + grid_energy_export[t]) for t in range(period)), name='powerbalance')
-    # 2. 添加发电机最大/最小功率约束
+    # 2. add constrain for pmax pmin
     m.addConstrs((gen_output[g, t] <= on_off[g, t] * p_max[g] for g in range(NUM_GEN) for t in range(period)),
                  'gen_output_max')
     m.addConstrs((gen_output[g, t] >= on_off[g, t] * p_min[g] for g in range(NUM_GEN) for t in range(period)),
                  'gen_output_min')
-    # 3. 添加上升和下降约束
+    # 3. add constrain for ramping up ramping down
     m.addConstrs((gen_output[g, t + 1] - gen_output[g, t] <= ramping_up[g] for g in range(NUM_GEN)
                   for t in range(period - 1)), 'ramping_up')
     m.addConstrs((gen_output[g, t] - gen_output[g, t + 1] <= ramping_down[g] for g in range(NUM_GEN)
                   for t in range(period - 1)), 'ramping_down')
-    # 4. 添加电池容量约束
+    # 4. add constrains for SOC
     m.addConstr(battery_capacity * soc[0] == battery_capacity * initial_soc +
                 (battery_energy_change[0] * battery_efficiency), name='soc0')
     m.addConstrs((battery_capacity * soc[t] == battery_capacity * soc[t - 1] +
                   (battery_energy_change[t] * battery_efficiency) for t in range(1, period)), name='soc update')
-    # 设置成本函数
-    # 发电机成本
+    # 5. add constrain for pv output
+    m.addConstrs((pv[t] >= 0 for t in range(period)), name='pv_output_min')
+    # set cost function
+    # 1 cost of generator
     cost_gen = gp.quicksum(
         (a_para[g] * gen_output[g, t] * gen_output[g, t] + b_para[g] * gen_output[g, t] + c_para[g] * on_off[g, t]) for
         t in range(period) for g in range(NUM_GEN))