Fix P/Q targets of converter stations in optimization

docs/optimizer/inputs.md

@@ -3,7 +3,12 @@
 ## Network data
 
 Files with the prefix `ampl_` contain the data and the parameters of the network on which the reactive OPF is executed.
-These files are obtained by using the [V2 of the extended version of PowSyBl AMPL export](https://github.com/powsybl/powsybl-core/blob/main/ampl-converter/src/main/java/com/powsybl/ampl/converter/version/ExtendedAmplExporterV2.java), which is the default version.
+
+These files are obtained by using the [V2 of the extended version of PowSyBl AMPL export](https://github.com/powsybl/powsybl-core/blob/main/ampl-converter/src/main/java/com/powsybl/ampl/converter/version/ExtendedAmplExporterV2.java), which is the default version.  
+
+Note that from release `0.10.0` of OpenReac, the active target of VSC and LCC converter stations is calculated using HVDC line active set point 
+and the converter mode, both specified in `ampl_network_hvdc.txt`. The losses related to rectifier/inverter conversion, and HVDC line **are ignored**.
+
 
 ## Configuration of the run
 

open-reac/src/main/resources/openreac/acopf.mod

@@ -180,9 +180,9 @@ sum{(qq,k,n) in BRANCHCC} base100MVA * V[k] * Red_Tran_Act_Dir[qq,k,n]
 # Loads
 + sum{(c,k) in LOADCC} load_PFix[1,c,k]     # Fixed value
 # VSC converters
-+ sum{(v,k) in VSCCONVON} vscconv_P0[1,v,k] # Fixed value
++ sum{(v,k) in VSCCONVON} vscconv_targetP[v] # Fixed value
 # LCC converters
-+ sum{(l,k) in LCCCONVON} lccconv_P0[1,l,k] # Fixed value
++ sum{(l,k) in LCCCONVON} lccconv_targetP[l] # Fixed value
 = 0; # No slack variables for active power. If data are really too bad, may not converge.
 
 
@@ -220,7 +220,7 @@ sum{(qq,k,n) in BRANCHCC} base100MVA * V[k] * Red_Tran_Rea_Dir[qq,k,n]
 # VSC converters
 - sum{(v,k) in VSCCONVON} vscconv_qvar[v,k]
 # LCC converters
-+ sum{(l,k) in LCCCONVON} lccconv_Q0[1,l,k] # Fixed value
++ sum{(l,k) in LCCCONVON} lccconv_q0[1,l,k] # Fixed value
 # Slack variables
 + if k in BUSCC_SLACK then
 (- base100MVA * V[k]^2 * slack1_shunt_B[k]  # Homogeneous to a generation of reactive power (condensator)

open-reac/src/main/resources/openreac/commons.mod

@@ -72,25 +72,56 @@ set BATTERYCC := setof {(1,b,n) in BATTERY : n in BUSCC} (b,n);
 # Warning: units with Ptarget=0 are considered as out of order
 set UNITON := {(g,n) in UNITCC : abs(unit_Pc[1,g,n]) >= Pnull};
 
+# Active and reactive targets of converter stations
+# Warning: the losses are ignored
+set LCCCONV_NUM := setof{(t,lcc,bus) in LCCCONV}lcc;
+set VSCCONV_NUM := setof{(t,vsc,bus) in VSCCONV}vsc;
+param lccconv_targetP {LCCCONV_NUM};
+param vscconv_targetP {VSCCONV_NUM};
+for {(1,h) in HVDC} {
+  # case of VSC converter stations
+  if (hvdc_type[1,h] == 1) then {
+    if (hvdc_convertersMode[1,h] == "SIDE_1_RECTIFIER_SIDE_2_INVERTER") then {
+      let vscconv_targetP[hvdc_conv1[1,h]] := hvdc_targetP[1,h];
+      let vscconv_targetP[hvdc_conv2[1,h]] := -hvdc_targetP[1,h];
+    } else {
+      let vscconv_targetP[hvdc_conv1[1,h]] := -hvdc_targetP[1,h];
+      let vscconv_targetP[hvdc_conv2[1,h]] := hvdc_targetP[1,h];
+    }
+  }
+  # case of LCC converter stations
+  if (hvdc_type[1,h] == 2) then {
+    if (hvdc_convertersMode[1,h] == "SIDE_1_RECTIFIER_SIDE_2_INVERTER") then {
+      let lccconv_targetP[hvdc_conv1[1,h]] := hvdc_targetP[1,h];
+      let lccconv_targetP[hvdc_conv2[1,h]] := -hvdc_targetP[1,h];
+    } else {
+      let lccconv_targetP[hvdc_conv1[1,h]] := -hvdc_targetP[1,h];
+      let lccconv_targetP[hvdc_conv2[1,h]] := hvdc_targetP[1,h];
+    }
+  }
+}
+check {lcc in LCCCONV_NUM}: lccconv_targetP[lcc] != NaN;
+check {vsc in VSCCONV_NUM}: vscconv_targetP[vsc] != NaN;
+
 #
 # VSC converter stations
 #
 set VSCCONVON := setof{(t,v,n) in VSCCONV:
   n in BUSCC
-  and abs(vscconv_P0[t,v,n]  ) <= PQmax
+  and abs(vscconv_targetP[v])  <= PQmax
   and abs(vscconv_Pmin[t,v,n]) <= PQmax
   and abs(vscconv_Pmax[t,v,n]) <= PQmax
-  and vscconv_P0[t,v,n] >= vscconv_Pmin[t,v,n]
-  and vscconv_P0[t,v,n] <= vscconv_Pmax[t,v,n]
+  and vscconv_targetP[v] >= vscconv_Pmin[t,v,n]
+  and vscconv_targetP[v] <= vscconv_Pmax[t,v,n]
   } (v,n);
 
 #
 # LCC converter stations
 #
 set LCCCONVON := setof{(t,l,n) in LCCCONV:
   n in BUSCC
-  and abs(lccconv_P0[1,l,n]) <= PQmax
-  and abs(lccconv_Q0[1,l,n]) <= PQmax
+  and abs(lccconv_targetP[l]) <= PQmax
+  and abs(lccconv_q0[1,l,n])  <= PQmax
   } (l,n);
 
 

open-reac/src/main/resources/openreac/connected_component.run

@@ -122,16 +122,16 @@ if 1 in LOG_INFO then {
 let temp1 := sum{(c,n) in LOADCC} load_PFix[1,c,n];
 let temp2 := sum{(g,n) in UNITON} unit_Pc[1,g,n];
 let temp2 := temp2 + sum{(b,n) in BATTERYCC} battery_p0[1,b,n];
-let temp3 :=  (sum{(vscconv,n) in VSCCONVON} vscconv_P0[1,vscconv,n])+(sum{(l,k) in LCCCONVON} lccconv_P0[1,l,k]);
+let temp3 :=  (sum{(vscconv,n) in VSCCONVON} vscconv_targetP[vscconv])+(sum{(l,k) in LCCCONVON} lccconv_targetP[l]);
 let global_initial_losses_ratio := (temp2-temp1-temp3)/(temp1+temp3);
 
 printf{LOG_INFO} "HVDC injections (homogeneous to loads):\n";
 for {(v,n) in VSCCONVON}
   printf{LOG_INFO} "VSC converter %Q in %Q: P0=%.1fMW is fixed, Q is variable\n",
-  vscconv_id[1,v,n],substation_id[1,bus_substation[1,n]],vscconv_P0[1,v,n];
+  vscconv_id[1,v,n],substation_id[1,bus_substation[1,n]],vscconv_targetP[v];
 for {(l,n) in LCCCONVON}
   printf{LOG_INFO} "LCC converter %Q in %Q: P0=%.1fMW is fixed, Q0=%.1fMvar is fixed\n",
-  lccconv_id[1,l,n],substation_id[1,bus_substation[1,n]],lccconv_P0[1,l,n],lccconv_Q0[1,l,n];
+  lccconv_id[1,l,n],substation_id[1,bus_substation[1,n]],lccconv_targetP[l],lccconv_q0[1,l,n];
 printf{LOG_INFO} "Sum of HVDC conv.  H: %.0f MW\n", temp3;
 printf{LOG_INFO} "Sum of loads       C: %.0f MW\n", temp1;
 printf{LOG_INFO} "Sum of generations P: %.0f MW\n", temp2;

open-reac/src/main/resources/openreac/dcopf.mod

@@ -53,8 +53,8 @@ subject to ctr_balance{PROBLEM_DCOPF, n in BUSCC}:
   + sum{(c,n) in LOADCC} load_PFix[1,c,n]
   + sum{(qq,n,m) in BRANCHCC} activeflow[qq,n,m] # active power flow outgoing on branch qq at bus n
   - sum{(qq,m,n) in BRANCHCC} activeflow[qq,m,n] # active power flow entering in bus n on branch qq
-  + sum{(vscconv,n) in VSCCONVON} vscconv_P0[1,vscconv,n]
-  + sum{(l,n) in LCCCONVON} lccconv_P0[1,l,n]
+  + sum{(vscconv,n) in VSCCONVON} vscconv_targetP[vscconv]
+  + sum{(l,n) in LCCCONVON} lccconv_targetP[l]
   =
   balance_pos[n] - balance_neg[n];
 

open-reac/src/main/resources/openreac/dcopf.run

@@ -107,7 +107,7 @@ for {(g,n) in UNITON : abs(P_dcopf[g,n]-unit_Pc[1,g,n]) >= temp2*0.99}
 let temp1 := sum{(c,n) in LOADCC} load_PFix[1,c,n];
 let temp2 := sum{(g,n) in UNITON} P_dcopf[g,n];
 let temp2 := temp2 + sum{(b,n) in BATTERYCC} battery_p0[1,b,n];
-let temp3 :=  (sum{(vscconv,n) in VSCCONVON} vscconv_P0[1,vscconv,n])+(sum{(l,k) in LCCCONVON} lccconv_P0[1,l,k]);
+let temp3 :=  (sum{(vscconv,n) in VSCCONVON} vscconv_targetP[vscconv])+(sum{(l,k) in LCCCONVON} lccconv_targetP[l]);
 printf{LOG_INFO} "Sum of HVDC conv.  H: %.0f MW\n", temp3;
 printf{LOG_INFO} "Sum of loads       C: %.0f MW\n", temp1;
 printf{LOG_INFO} "Sum of generations P: %.0f MW\n", temp2;

open-reac/src/main/resources/openreac/reactiveopfoutput.run

@@ -115,7 +115,7 @@ printf{(v,n) in VSCCONVON} "%i;%i;%i;%Q;%.3f;%.1f;%.1f;%.1f;\n",
   vscconv_vregul[1,v,n],
   V[n],
   vscconv_qvar[v,n],
-  vscconv_P0[1,v,n],
+  vscconv_targetP[v],
   vscconv_qvar[v,n]
   > (fileOut);
 close (fileOut);