Evaluation of the Simultaneous Operation of the Mechanisms for Cross-Border Interchange and Activation of the Regulating Reserves

: This article examines the mechanisms for cross-border interchange of the regulating reserves (RRs), i.e., the imbalance-netting process (INP) and the cross-border activation of the RRs (CBRR). Both mechanisms are an additional service of frequency restoration reserves in the power system and connect different control areas (CAs) via virtual tie-lines to release RRs and reduce balancing energy. The primary objective of the INP is to net the demand for RRs between the cooperating CAs with different signs of interchange power variation. In contrast, the primary objective of the CBRR is to activate the RRs in the cooperating CAs with matching signs of interchange power variation. In this way, the ancillary services market and the European balancing system should be improved. However, both the INP and CBRR include a frequency term and thus impact the frequency response of the cooperating CAs. Therefore, the impact of the simultaneous operation of the INP and CBRR on the load-frequency control (LFC) and performance is comprehensively evaluated with dynamic simulations of a three-CA testing system, which no previous studies investigated before. In addition, a function for correction power adjustment is proposed to prevent the undesirable simultaneous activation of the INP and CBRR. In this way, area control error (ACE) and scheduled control power are decreased since undesired correction is prevented. The dynamic simulations conﬁrmed that the simultaneous operation of the INP and CBRR reduced the balancing energy and decreased the unintended exchange of energy. Consequently, the LFC and performance were improved in this way. However, the impact of the INP and CBRR on the frequency quality has no unambiguous conclusions.


Motivation and Literature Review
The mechanisms for cross-border interchange and activation of the regulating reserves (RRs) are evolving due to the expensive balancing energy, and are included in the European Union's current regulations [1,2]. They are put in operation in continental Europe by the members of the European network of transmission system operators for electricity (ENTSO-E) [3]. Since the first implementation of the cross-border interchange of the RRs, i.e., the imbalance-netting process (INP), in 2008, the cumulative value of all the netted imbalances amounted to more than €600 million by the third quarter of 2020 [4]. The total monthly volume of netted imbalances for September 2020 was 698.69 GWh, which amounts to €13.38 million. Moreover, the monthly avoided positive and negative RRs activations amount to a minimum of 10% to as much as 85%. The further development of the INP with the functionality of the cross-border activation of the RRs (CBRR), will additionally reduce the activation of the RRs and increase the associated savings. There has been a surge in the application of machine learning and statistical framework to solve similar problems focused in this paper. The authors in [24] explore the influencing factors of consumer purchase intention of cross-border e-commerce based on a wireless network and machine learning in order to provide decision support for the operation of e-commerce and to promote the better development of cross-border e-commerce. Several model-based experimental design techniques have been developed for design in domains with partial available data about the underlying process. The authors in [25] focus on a powerful class of model-based experimental design called the mean objective cost of uncertainty. To achieve a scalable objective-based experimental design, a graph-based mean objective cost of uncertainty with Bayesian optimization framework is proposed. A thorough review of the issues of data localization and data residency is given in [26], in addition to clarifying cross-border data flow restrictions and the impact of cross-border data flows in Asia.

Contribution and Structure of the Paper
Generally, the INP and the CBRR should have a positive impact on the LFC and performance. However, the quality of the frequency is continually decreasing [27]. Therefore, the impact of the INP and the CBRR on the frequency quality, the LFC, and performance in a three-CA test system was examined separately in [20,22] with dynamic simulations. In [21], the impact of INP on power-system dynamics is shown and an eigenvalue analysis of a two-CA system is conducted. The impact of CBRR on the power-system dynamics is shown in [23], and a modified implementation of the CBRR is proposed that has no impact on the system's eigendynamics. This article extends these earlier results with an in-depth evaluation of the simultaneous operation of the mechanisms for cross-border interchange and activation of the RRs, which was not studied before. Dynamic simulations are performed for all the cases where the simultaneous operation of the INP and the CBRR is possible. In addition, a function for correction-power adjustment is proposed as one of the contributions of this paper, since small delays in demand power sign change could cause undesired simultaneous activation of the INP and the CBRR. In this way, ACE and scheduled control power are decreased, since udesired correction is prevented. A basic schematic diagram of the LFC, INP, and CBRR and the function for correction power adjustment is given in Figure 1 (right), where the order of operation is clearly seen. Note here that the correction power is the output of the adjustment block. As far as we know, no researchers have examined the impact of the simultaneous operation of the INP and CBRR on the LFC and performance.
This article consists of the following parts: In Section 2, the elemental concepts of the LFC, the INP and the CBRR are described. Simultaneous operation of the INP and the CBRR is also described. Additionally, a function for correction power adjustment is proposed as one of the main contributions of this article, which prevents undesired correction. Section 3 describes indicators for evaluation of the frequency quality, LFC and performance, rate of change of frequency (RoCoF), balancing energy, unintended exchange of energy and energy exchange. In Section 4, a three-CA test system with the INP and the CBRR is described. Two types of test cases were performed with the dynamic simulations, i.e., step change of the load and the random load variation. The primary contribution of this article is given in Section 5, where the impact is given of the simultaneous operation of the INP and the CBRR on the frequency quality, the LFC, and performance. Lastly, Section 6 outlines the main conclusions and outlines future work.

LFC
An interconnected power system consists of a large number of CAs that are connected via tie-lines. Each individual TSO maintains the frequency of each CA within predefined standard limits and the tie-line power flows with neighboring CAs within prespecified tolerances, which are the main objectives of the LFC [28]. This is generally accomplished by reducing the ACE, which is, for the CA i, defined as where ∆P i = (P ai −P si ) and ∆ f i = ( f ai − f si ) denote interchange power variation and frequency deviation, respectively. Note that P ai and f ai denote actual values, while P si and f si denote scheduled values. Furthermore, B i is the frequency-bias coefficient [1]. Note that ACE i > 0 denotes that the generation is higher than the load; hence, the CA is denoted as "long". Similarly, a CA is denoted as "short" when ACE i < 0. The basic LFC framework of the CA i is shown in Figure 2, where LPF denotes a low pass filter and SH a sample and hold, with the value of a sampling time T s = 2 s. A negative control-feedback is characterized as −1 gain and PI denotes a proportionalintegral controller. Scheduled control power ∆P sci denotes the output of LFC, which is appointed to the participating control units that change the electrical control power ∆P ei appropriately. Neglecting the losses, then ∆P ei is, for the CA i, defined as where ∆P Li denotes the load-power variation. Note that ∆P ei is well-known as the balancing energy, whereas, instead of LFC reserve, the term RR is generally used. Figure 2. Schematic block diagram representation of the LFC in the CA i.

INP
The basic principle of the INP operation and a steady-state correction value calculation with the INP for three CAs is given in [20]. A control-demand concept is used for the INP and a schematic block diagram representation is shown in Figure 3   N can be different. The cooperating CAs are connected to the "red" summator, forming virtual tie-lines, as shown in Figure 3. The input variables to the "red" summator are the demand powers of the cooperating CAs, i.e., P di , P dj , . . . , P dN . The demand power of the CA i characterizes the maximum interchange power for the CA i among the cooperating CAs and is defined as according to [3,29]. Introducing (1) and (2) in (3) gives the following relation: The power imbalance between generation and load in addition to The output variables of the "red" summator are the correction powers of the cooperating CAs, i.e., P INP cori , P INP corj , . . . , P INP corN , determined with a delay of T s due to the SH. The correction power of the CA i characterizes the maximum interchange power for the CA i among the cooperating CAs with a different sign of ACE i , and is included as where the terms in brackets denote ACE i . Moreover, only CAs with different signs of demand power, i.e., sign(P di ) = sign(P dj ), can net the demand for balancing energy. If two or more cooperating CAs are "short", then CBRR is used instead of the INP and vice versa [23]. Therefore, the cooperating CAs must be "short" and "long" in order to net the demand power through the INP. Hence, the balancing energy in CAs that net the balancing energy from the cooperating CAs can be reduced, and simultaneously the RR is released. The P INP cori , P INP corj , . . . , P INP corN is determined by the INP optimization module, considering numerous target functions, as given in [20].
Considering N CAs, then the P INP cori is, for the CA i, expressed as Considering (4), then the P INP cori between N CAs is, for the CA i, expressed as In this way, the correction power of the CA i compensates the load variation that is varied by the frequency variation of the cooperating CAs. From a system point of view, this corresponds to additional frequency-based feedback and cross-couplings with cooperating CAs, which inseparably changes the eigendynamics of the CA i [21].

CBRR
The basic principle of the CBRR operation and a steady-state correction-value calculation with the CBRR for three CAs is given in [22]. The same control-demand concept is used for the CBRR as is currently used for the INP, and a schematic block diagram representation is shown in Figure 4 Figure 4.
The input variables to the "red" summator are the demand powers of the cooperating CAs, i.e., P di , P dj , . . . , P dN . The demand power of the CA i characterizes the maximum activation power for the CA i among the cooperating CAs, defined as (3), according to [3,29].  (1) and (2) in (3), (4) is obtained. The power imbalance between generation and load in addition to K CBRR

Similar to INP, by introducing
The output variables of the "red" summator are the correction powers of the cooperating CAs, i.e., P CBRR cori , P CBRR corj , . . . , P CBRR corN , determined with a delay of T s due to the SH. The correction power of the CA i characterizes the maximum activation power for the CA i among the cooperating CAs with matching sign of ACE i , and is included as where the terms in brackets denote ACE i . Moreover, only CAs with matching sign of demand power, i.e., sign(P di ) = sign(P dj ), can activate the demand for balancing energy. If any of the cooperating CAs are "long" and the others are "short", then INP is used instead of the CBRR and vice versa [21]. Therefore, the cooperating CAs must be either "short" or "long", depending on whether a positive or negative CBRR is activated. Hence, the balancing energy in CAs that activates the balancing energy in the cooperating CAs can be reduced, and simultaneously the RR is released. The P CBRR cori , P CBRR corj , . . . , P CBRR corN is determined by the CBRR optimization module, considering numerous target functions, as given in [22].
Considering N CAs, then the P CBRR cori is, for the CA i, expressed as Considering (4), then the P CBRR cori between N CAs is, for the CA i, expressed as Similar to the INP, the correction power of the CA i compensates the load variation that is varied by the frequency variation of the cooperating CAs. From a system point of view, this corresponds to an additional frequency-based feedback and cross-couplings with cooperating CAs, which inseparably changes the eigendynamics of the CA i [23].

Simultaneous Operation of the INP and the CBRR
In the cooperating CAs, ∆P Li changes randomly and continuously. Consequently, cases of ∆P di sign changes can occur, resulting in a continuous switching between the INP and the CBRR, which causes undesirable ∆P cori sign changes. Such a situation might occur when the signs of ∆P di and ∆P dj are changed with a short time delay. Therefore, a function for P cori adjustment is proposed as one of the contributions of this article.
A schematic block diagram representation for P cori adjustment in relation to P di of the CA i is shown in Figure 5. The signs of the two successive samples, i.e., P di,k−1 and P di,k , are compared with the relational operator, whose output is connected to a switch, marked with "c". Two states are possible, i.e., • State 1: sign(P di,k−1 ) = sign(P di,k ) and • State 2: sign(P di,k−1 ) = sign(P di,k ).
For State 1, the switch position is "1" and the output variable is P ADJ cori = P cori . For State 2, the switch position is "2" and the output variable is P ADJ cori = 0.
Relational operator Switch ADJUST ADJ cori P P cori Figure 5. Schematic block diagram representation for P cori adjustment in the CA i.
The implementation of the P cori adjustment in the INP and the CBRR framework is shown in Figure 6. The simultaneous operation of the INP and the CBRR, considering the P cori adjustment, can be described using the pseudo-code, as shown in Alogorithm 1.

Algorithm 1:
The code of simultaneous operation of the INP and the CBRR. Figure 6. Schematic block diagram representation for P cori adjustment with the INP and the CBRR in the CA i.
An example of the simultaneous operation of the INP and the CBRR, with and without the P cori adjustment, is shown in Figure 7. The values of the loads were set in such a way that the simultaneous operation of the INP and the CBRR was possible. In the time interval of 80-100 s, the INP operated between CA 1 -CA 3 and CA 2 -CA 3 , whereas the CBRR operated between CA 1 -CA 2 . During this time interval, P cori due to the INP was possible in all three CAs, whereas P cori due to CBRR was possible only in CA 1 and CA 2 . At t =100 s, a simultaneous step change of the load in CA 1 and CA 3 was applied, whereas, in CA 2 , it was applied with a delay, i.e., at t = 100.05 s, as seen at t = 102 s due to SH with T s = 2 s. Time responses of P di and P cori at t =100 s with one step-size activation of the INP in CA 1 and CA 3 , and the CBRR in CA 3 is clearly seen in Figure 7a (without P cori adjustment). Note that P cori due to CBRR should be zero in CA 3 . Additionally, at t =102 s, one step-size activation of the CBRR in CA 2 is also seen. However, in Figure 7b (with P cori adjustment), at t = 100 s and t = 102 s, the value of P cori was zero in all CAs. In this way, delayed P di sign changes have no impact on the switching between the INP and the CBRR. Moreover, the P cori variation is significantly reduced using the proposed adjustment.  Figure 7. Time responses of P di and P cori for a three-CA testing system without P cori adjustment (a) and with P cori adjustment (b), where "w INP" is with the INP and "w CBRR" is with the CBRR.

Indicators for Evaluation of LFC, INP, and CBRR Provision
The impact of the simultaneous operation of the INP and CBRR on the frequency quality, the LFC, and performance is evaluated using 15-min averages [1].

Performance Indicators
Frequency quality is evaluated with the standard deviation of ∆ f i , denoted as σ ∆ f i . In addition, the LFC and performance is evaluated with the standard deviation of ACE i , denoted as σ ACEi [1,30].

RoCoF
RoCoF is the time derivative of the power system's frequency, i.e., d f i dt [31]. The mean value of RoCoF i , denoted as µ RoCoF i , is evaluated individually for positive and negative values, denoted as µ RoCoF i+ and µ RoCoF i− .

Standard Deviation and Mean Value of RRs
RRs assist in active power balance to correct the imbalance in the transmission grid and lead the power system frequency to the normal frequency range [27]. The standard deviation and mean value of ∆P sci are calculated individually for positive and negative values, denoted as σ ∆P sci+ , σ ∆P sci− and µ ∆P sci+ , µ ∆P sci− .

Balancing Energy
The balancing energy enables TSOs to cost-effectively compensate for power and voltage variation in the transmission grid [2]. It is the actual electrical control power that is, for a particular period of time, calculated as ∆W ei = t 0 ∆P ei dt. Individual positive and negative values, denoted as ∆W ei+ and ∆W ei− , are calculated.

Unintended Exchange of Energy
The unintended exchange of energy is determined by the difference between interchange power variation and correction power [2] that is, for a particular period of time, calculated as ∆W uni = t 0 (∆P i − P cori )dt. Individual positive and negative values, denoted as ∆W uni+ and ∆W uni− , are calculated.

Energy Exchange
The energy exchange through the INP and CBRR is defined as the actual interchanged or activated power between cooperating CAs [32], that is, for a particular period of time, calculated as W cori = t 0 P cori dt. Additionally, positive and negative values are calculated, individually for the INP, denoted as W INP cori+ , W INP cori− , and individually for the CBRR, denoted as W CBRR cori+ , W CBRR cori− .

Dynamic Simulations
A three-CA test system was used for the dynamic simulations, where CA 1 -CA 2 and CA 2 -CA 3 were connected by physical tie-lines, whereas CA 1 -CA 3 were not connected with a tie-line. Moreover, all three CAs were connected with virtual tie-lines due to the INP and the CBRR. A Matlab/SIMULINK model was used, where the dynamic simulations were performed with a 50 ms step size.

Dynamic Model
The basic schematic block diagram representation of a single CA, characterized with a linearized, low-order, time-invariant model, is shown in Figures 8 and 9 [29,33]. Note that the INP and CBRR implementation is not shown. The generator-load dynamic is described by the rotor inertia H i and the damping D i . Moreover, three different types of the turbine-governor systems were considered, i.e., hydraulic, steam reheat, and steam non-reheat. A constant droop characteristic R ni was assumed. In addition, the ramping rate and the participation factors α ni of the control units were also taken into account. The tie-line between the connected CAs is described by the synchronizing coefficient T ij [34]. Furthermore, a 1st-order LPF is modeled by a time constant T LPFi , while the PI controller is modeled by a gain K ri and a time constant T ri .

Parameters
Common parameter values were set for a three-CA test system as given in Table 1 [28,29]. Note that the frequency-bias coefficient was determined as a constant, i.e., B i = 1/R 1i + 1/R 2i + 1/R 3i + D i . The model parameters were set equally for all three CAs. One cycle of the LFC, INP, and CBRR was incorporated with T s =2 s.

Test Cases
The maximum possible compensation with the INP and CBRR was considered. Dynamic simulations were performed so that the loads of individual CAs were altered during the simulation. In addition, two types of test cases were performed, i.e., step change of the load and the random load variation.
Moreover, the inertia time constant H i , the tie-line parameter T ij , and the droop characteristic R i have a considerable impact on frequency quality according to [35,36]. Therefore, different values of H i , T ij , and R i were used to show the impact of the simultaneous operation of the INP and CBRR on the indicators for evaluation of LFC, INP, and CBRR provision.

Step Change of Load
The values of the loads were set in such a way that the simultaneous operation of the INP and the CBRR was possible. At t = 0 s, a simultaneous step change of the loads was applied and the magnitudes were set as ∆P L1 = 0.06 pu, ∆P L2 = 0.07 pu and ∆P L3 = −0.08 pu. In addition, at t =100 s, the magnitudes were set as ∆P L1 = −0.06 pu and ∆P L3 = 0.08 pu, whereas, at t = 100.05 s, the magnitude was set as ∆P L2 = −0.07 pu. The resulting load is shown in Figure 10a. Consequently, the INP operated between CA 1 -CA 3 and CA 2 -CA 3 , while the CBRR operated only between CA 1 -CA 2 . Note that this case is used in Section 2.4. Step changes of ∆P Li (a) and random ∆P Li variations (b1-b3) for a three CA test system.

Random Load Variation
The random load was modeled as a linear, stochastic, time-invariant, first-order system with two components [37]. A low-frequency component captures the trend changes with a quasi-period of 10-30 min, whereas the residual component captures fluctuations with a quasi-period of several minutes. Measurements of an open-loop ACE in an undisclosed CA were used to determine the model parameters. The resulting normalized load for a three-CA test system with different signs is shown in Figure 10(b1-b3), where the random load was changed every 60 s for 24 h. The statistical parameters of the random loads for all three CAs are given in Table 2, where µ Li and σ Li denote the meanvalue and the standard deviation of the i-th load, while the correlation between the i-th and j-th loads is denoted as ρ Lij . The correlation ρ Lij is extremely small, which allows both the INP and CBRR to operate simultaneously.

Results
Dynamic simulations with and without the INP and CBRR were performed for a three-CA test system. The impact of the simultaneous operation of the INP and CBRR on the frequency quality, the LFC, and performance was evaluated with the obtained results. Note that the results shown in this section refer to a three-CA test system, whereas the basic principle is applicable to N CAs as shown in Section 2. In addition, the results cannot be generalized to the dynamics of the INP and CBRR.

Step Change of Load
The time responses to the step change of ∆P Li are shown in Figures 11-13. In Figure 11 (left), it is clear that the frequency deviations ∆ f i in all three CAs appeared following a step change of ∆P Li that was applied. After the first step change, ∆ f 1 and ∆ f 2 were negative due to the positive step change of ∆P Li , whereas ∆ f 3 was positive due to the negative step change of ∆P Li . Note that after the second step change, the signs were opposite. The primary frequency control decreased |∆ f i | in about 15-25 s after the step change of ∆P Li ; then, LFC decreased |∆ f i | slowly. The results show that the impact of the INP and CBRR on ∆ f i is not significant.
The impact of the INP and CBRR is shown more obviously in Figures 11 (right) and 12 (left). In all three CAs, the values of ACE i , ∆P sci and ∆P ei were decreased with the INP and CBRR. Furthermore, the INP and CBRR clearly increased ∆P i , due to the increased tie-line power flow between the CAs.
The signs of P di and P cori are opposite, as shown in Figure 12 (right). A 2 s time delay is seen, due to SH with T s =2 s. Note that the time responses of P di and P cori were already described in Section 2.4 and Figure 7b, where the same case was performed.
The time responses ACE i and ∆P sci with and without the P cori adjustment function are shown in Figure 13. Without the P cori adjustment, ACE i was increased and ∆P sci was undesirably increased, due to switching between the INP and the CBRR, which caused undesirable ∆P cori sign changes. Clearly, in all three CAs, the values of ACE i and ∆P sci were decreased with the P cori adjustment. Note that, in Figure 13, the same example is performed as described in Section 4.2.1 and shown in Figures 11 and 12. The difference can be seen because, in Figures 11 and 12, the comparison without and with the INP and CBRR is shown, whereas, in Figure 13, the comparison without and with the P cori adjustment is shown.  Figure 13. Time responses of ACE i (left) and time responses of ∆P sci (right) for a three-CA test system, where "wo ADJ" is without and "w ADJ" is with the P cori adjustment function.

Random Load Variation
Simulations of the simultaneous operation of the INP and CBRR were performed to show the impact of an individual mechanism. The simulations were also performed separately, i.e., operation of only the INP or only the CBRR. The results are given in Tables 3-5. There are no unambiguous conclusions about the impact of the mechanisms on σ ∆ f i , and the differences are extremely small, as shown in Table 3. However, both mechanisms reduce σ ACEi , INP slightly more than CBRR. The reduction is most pronounced when both mechanisms operate simultaneously. Both mechanisms reduce the |µ RoCoFi+ | and |µ RoCoFi− |, with the reduction being most pronounced when both mechanisms operate simultaneously according to Table 3. However, there are no unambiguous conclusions as to which mechanism reduces |µ RoCoFi+ | and |µ RoCoFi− | more, and, in most cases, it is the INP.
When both mechanisms operate simultaneously, |µ Psci+ | and |µ Psci− | are greatly reduced according to Table 3. However, the results of the separate operation of the mechanisms show that the impact of INP is greater than the impact of CBRR. This is expected, as the CBRR only activates the RRs in the cooperating CAs. Moreover, the reduction of |σ Psci+ | and |σ Psci− | is most noticeable when both mechanisms operate simultaneously. However, the results of the separate operation of the mechanisms show that the impact of the INP is greater than the impact of the CBRR.
The conclusions for |∆W ei+ | and |∆W ei− | are similar to |µ Psci+ | and |µ Psci− | according to Table 4, which is expected, as this indicator describes the response of the control units.
When both mechanisms operate simultaneously, |∆W uni+ | and |∆W uni− | are reduced according to Table 4, except in one case where only |∆W uni+ | was increased, while |∆W uni− | was reduced considerably. Furthermore, the results of the separate operation of the mechanisms show that the INP almost completely eliminates unintentional deviations, while the impact of the CBRR is not very pronounced.
When both mechanisms operate simultaneously, |W INP cori+ |, |W INP cori− |, |W CBRR cori+ | and |W CBRR cori− | are slightly reduced compared to separate operation of the INP and CBRR according to Table 5. This is due to the P cori adjustment mechanism, which is only required in the case of simultaneous operation of the INP and CBRR.
Moreover, simulations of the simultaneous operation of the INP and CBRR were performed for different values of H i , T ij , and R i . The results are given in  There are no unambiguous conclusions about the impact of the mechanisms on σ ∆ f i , and the differences are extremely small, as shown in Figure 14 (left). In addition, the impact of H i and T ij is not clear, whereas the increase of R i results in an increase of σ ∆ f i . However, both mechanisms significantly reduce σ ACEi , as shown in Figure 14 (right), whereas H i , T ij and R i have no impact on σ ACEi .  Generally, both mechanisms reduce the |µ RoCoFi+ | and |µ RoCoFi− |, as shown in Figure 15. In addition, an increase of H i and R i results in a decrease of |µ RoCoFi+ | and |µ RoCoFi− |, whereas an increase of T ij results in an increase of |µ RoCoFi+ | and |µ RoCoFi− |. Note that, when T ij = 1/15 pu/Hz, |µ RoCoFi+ | and |µ RoCoFi− | is increased with the mechanisms.
When both mechanisms operate simultaneously, |µ Psci+ | and |µ Psci− | are greatly reduced, as shown in Figure 16. However, H i , T ij and R i have no impact on |µ Psci+ | and |µ Psci− |.  Figure 16. Average values of µ ∆P sci+ (left) and average values of µ ∆P sci− (right) for different values of H i , T ij and R i , where "wo INP/CBRR" is without and "w INP/CBRR" is with the INP and CBRR.

Conclusions
This article discusses the simultaneous operation of the INP and CBRR, which no previous studies have investigated. Extensive dynamic simulations of a three-CA test system with the simultaneous operation of the INP and CBRR were performed to evaluate their impact on the frequency quality, the LFC, and performance.
The results confirmed the conclusions in [20][21][22][23], where the INP and CBRR were analyzed separately. The results of the step change of load and the random load variation confirmed that the impact of the INP and CBRR on frequency deviations has no unambiguous conclusions. In addition, the impact of inertia and synchronizing coefficient is not clear, whereas the increase of droop characteristic results in an increase of frequency deviations. However, both mechanisms reduce the ACE deviations-the INP slightly more than CBRR. The reduction is most pronounced when both mechanisms operate simultaneously. Moreover, the function for correction power adjustment additionally decreased ACE and scheduled control power, which prevents undesirable switching between the INP and CBRR. Both mechanisms also reduce the RoCoF, and the reduction is most pronounced when both mechanisms operate simultaneously. However, there are no unambiguous conclusions as to which mechanism reduces the RoCoF more, and, in most cases, it is the INP. Increase of inertia and droop characteristic results in a decrease of RoCoF, whereas increase of synchronizing coefficients results in an increase of RoCoF. When both mechanisms operate simultaneously, the scheduled control power is greatly reduced and the impact of the INP is greater than the impact of the CBRR. Similarly, the balancing energy as well as the unintended exchange of energy are greatly reduced when both mechanisms operate simultaneously. However, there is no impact of inertia, synchronizing coefficient, and droop characteristic on scheduled control power. Due to the function for correction power adjustment, which prevents undesirable activation of the INP and CBRR, energy exchange was slightly reduced, as expected. Because of the reduced unintended exchange of energy, beneficial economic consequences can be anticipated when the INP and CBRR operate simultaneously.
One of the tough challenges for all researchers in this domain is the dynamic dimensioning of RRs, considering the INP and CBRR. This article clearly shows that the simultaneous operation of the INP and CBRR reduces the activation of the RRs, which is currently not considered in the reserve dimensioning process.  Data Availability Statement: This study did not report any data.

Conflicts of Interest:
The authors declare no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript: