The proposed approach was applied to a realistic large scale power system to assess the impact of PST of the branches’ power flows.
4.2. Reduced Model
The hierarchical clustering, as defined in
Section 2.1, was applied to aggregate the network, using 300 different LMP scenarios.
The LMP scenarios were calculated using 300 operating points corresponding to different levels of net load for the overall system to capture the seasonal characteristics of the European load profile. Those 300 periods were uniformly picked from the annual European load curve, whose total energy values per country can be found in
Table 3. During this calculation, all PSTs of network were considered as an optimization variable that could vary between −30 and +30 degrees.
The goal of the network reduction methodology was to keep the information regarding network congestions, with a much smaller representation. For sake of simplicity, the aggregation was performed separately for each existing bidding zone and only for the buses at its interior, avoiding clusters that would share different bidding zones. Despite that, the methodology allowed aggregating buses regardless of the bidding zone definition.
In addition, since the goal of this approach was to study the impact of PST, a supplementary constraint was added r to avoid the aggregation of areas connected by a PST. Therefore, the algorithm would stop when all observations at the interior of a zone were grouped except for those connected by a PST.
This resulted in an equivalent model with 54 buses and 82 branches. The connectivity between clusters were defined to match those of the complete model. The same input data were used to define the reduced PTDF matrix, as detailed in
Section 2.2.
4.3. Impact of PST Modeling
Section 2.4 demonstrates the suitability of the PSDF matrix to represent PST in reduced network models. The goal was to assess the accuracy of the representation in a large-scale power system with realistic values. Therefore, the same approach was used: the full European transmission network was clustered into a set of reduced buses and represented by a static PTDF matrix, as stated above, and then the same tap optimization was applied in both the full network model (with an explicit representation of the PST) and in the reduced model (using the PSDF matrix). The results in terms of branches’ power flows were compared to determine the best PST representation.
Given the high complexity of analyzing the entire power system, for simplification purposes, we focused on analyzing a single PST located at the border between Germany and the Netherlands. As detailed in
Section 3, three different cases ertr used to assess the proposed methodology:
is a static PTDF matrix calculated having as input and does not explicitly include any variable to represent the PSTs.
is a static PTDF matrix calculated having as input and does not explicitly include any variable to represent the PSTs.
+ includes the previous PTDF matrix and a matrix.
The generation plan was firstly obtained by running an OPF with the full model and was then used to perform a load flow with the reduced model. The resulting flows of both the OPF and load flow were compared to assess the accuracy of the modeling.
Table 6 presents the error of the flows obtained with the reduced network model for
and
PTDF matrix modeling. The error was calculated using Equation (
12) with the evaluation scenarios
.
Table 6 demonstrates the impact of the optimization of the PST in the reduced model. The
matrix performed well when the input scenarios did not consider the optimization of the PST, but the error tended to increase when the
was used. On the other hand, the
matrix could reduce the error for the case where the PST was optimized
, with a RMSE of only 138.5 MW per branch per scenario, but the error rapidly increased when
was applied.
The results in
Table 6 show that the static PTDF matrix representation performed well under the scenarios from which it was built. When the “non PST optimized” injected powers (
) were applied to the “PST optimized” PTDF matrix (
), the results are less accurate than when applied to the “non PST optimized” PTDF matrix (
) representation and vice versa.
The same trend can be observed in
Table 7, which presents the VaR of 5% for the flows calculated using the reduced model. It can be remarked that the values of VaR were similar for both
and
when using the
, but more significant values arose when applying
to the
matrix.
Figure 2 and
Figure 3 show the error distribution for the
and
PTDF matrix, respectively, when applying
and
. Following the same trend observed in
Table 6 and
Table 8, the results demonstrate that, when the operational scenarios with no PSTs modeled (
) were applied to the non PST optimized PTDF matrix (
), the errors were lower than when the opposite occurred. With these two figures, the over fitting process that occurred in the optimization process defined in
Section 2.2 became clear, therefore, stressing the need for a more general representation of PSTs in the reduced models.
Finally, the reduced model including an explicit modeling of the PST (
+
) was tested. For the injected power set where PSTs were optimized (
), the PST of the reduced model mimicked its behavior; in other words, an injected power was multiplied by the
matrix. As presented in
Section 2.4, a relationship between degrees of the PST and injected power was established and the same power was injected by the PST in the reduced model.
Table 8 shows the error results for the case where an explicit modeling of the PST was done (
+
). As can be observed, for the case where the PST was not optimized (
), the error was the same as presented in
Table 6, as the injected power of the PST of the reduced model was set to zero.
When considering the case where the PST was optimized (), the error showed a slight reduction issue to the explicit PST modeling. In addition, when looking into the VaR, it was observed that, with the explicit modeling, the VaR tended to be similar independently of the considered case.
Figure 4 shows the error distribution for the representation using a PTDF + PSDF matrix. It was observed that the error obtained with the
was similar to the one presented in
Figure 2. For the
, the error was superior to the one obtained with the
PTDF matrix (as shown in
Figure 3) but inferior to the ones obtained with the
PTDF matrix.
Comparing both the results obtained with and without the explicit PST modeling, it was observed that relying only on the optimization of the PTDF matrix using a set of data where the PSTs were optimized tended to underperform the case where the PST were not optimized and vice versa. When adding the explicit modeling of the PSTs, the results for the case using were not as accurate as the ones obtained with the matrix, but were more accurate for the case where was applied.
Overall, the proposed model with an explicit modeling of the PSTs lost some accuracy for a specific case
, but compensated for the opposite case
, as highlighted by the VaR values presented in
Table 8.