3.2. Future Energy Scenarios (FESs)
For a comprehensive analysis of the KSA’s future power system, three plausible FESs are developed based on the government’s renewables targets in the coming decade. In 2018, the KSA government issued a target to achieve 27.3 GW of local generation from RES by 2023 and 58.7 GW by 2030 [
21].
The 2030 target is considered to be the high progression scenario (HP). Another two FESs of the KSA’s power system with various penetrations of RESs have been developed considering the medium progression (MP) and low progression (LP) of the 2030 renewables target. These scenarios aim to detect when the overall system inertia will drop to critical levels that could endanger the grid’s stability. An estimation of the electricity generation mix was considered based on the available energy sources, renewables potential, and network capabilities. The energy mixes, and more importantly the synchronous and asynchronous generation capacity resulting from these FESs, are summarized in
Table 1 and
Table 2 for 2023 and 2030, respectively [
21,
24,
26].
The HP scenario represents the KSA government’s vision of achieving almost 25% of the generation mix from RESs by 2023 and 41% by 2030, as shown in
Table 1 and
Table 2, respectively. It is assumed the KSA government will make outstanding progress in achieving the projected target on time.
In the MP scenario, a medium progression in the projects of the 2030 renewables target is assumed. The MP scenario anticipates that 50% of the planned RES projects will be achieved by 2023 and 2030, respectively, due to financial and technical obstacles. It is assumed that non-synchronous generation will contribute to the generation mix by 13 GW and 28 GW in 2023 and 2030, respectively, as shown in
Table 1 and
Table 2.
The LP scenario assumes slow progression in the achievement of the 2030 renewables plan, due to some financial and technical obstacles to developing and integrating such large renewable energy projects into the existing power system. This scenario assumes that just 30% of the planned target will be achieved by 2023 and 2030, respectively. The RESs are expected to contribute to the energy mix with only 9.3 GW and 14 GW in 2023 and 2030, respectively, as shown in
Table 1 and
Table 2.
3.3. Simplified KSA Power System Model
The electrical frequency of a power system is the mean frequency of all of the synchronous generators that are directly coupled to the grid, which is directly proportional to their rotor speeds. Hence, the control of power system frequency can be expressed as the frequency control of one synchronous generator. The control of the system frequency is achieved by a turbine governor (speed governor) control. The speed governor adjusts the position of the main steam valve to maintain the speed of the rotor at the synchronous speed and ensure a balance between the mechanical and electrical torques [
30,
31].
For instance, a change in electrical load is reflected instantaneously as a change in the electrical torque output () of the generator. This change causes an imbalance between the turbine mechanical torque () and generator electrical torque (). The imbalance between and results in an imbalance in the net torque () and a change in the rotor speed (), according to the equation of motion. When is less than , decreases and the governor extends the valve inlet to increase . Similarly, when is greater than , increases and the governor shrinks the valve inlet to decrease .
The relationship between the net torque and the variation in the mechanical speed of the generator rotor can be described using the swing equation of the synchronous generator, as shown in Equation (1):
where
in [kg.m
2] is the combined moment of inertia of the generator and turbine.
The swing equation can be expressed in terms of inertia constant (H), which is defined as the ratio of the stored kinetic energy on the machine rotor (K.E) to the rating power of the generator (S
base) in [MVA], as follows:
Substituting for J in Equation (1) gives:
For a thorough study of the system frequency, it is more desirable to express this relationship in terms of mechanical and electrical power instead of torque. The relationship between
in [W] and torque
in [N.m] is shown in Equation (4), where
in [rad/sec] is the speed of the machine rotor [
30,
31]:
Considering a small deviation from the normal values, the above quantities can be rewritten as follows:
and equation can be rewritten considering their steady state and derivative terms, as shown in Equation (8).
The relationship between the derivative terms only can be written as follows:
In steady state,
and the synchronous speed (
) = 1 p.u.; hence, Equation (10) can be written as follows:
Equations (3) and (11) can now be expressed in terms of inertia constant as shown in Equation (12).
Loads in the power system are classified based on their sensitivity to the frequency deviation into two categories: frequency-independent loads, such as resistive loads (i.e., lighting and heating loads); and frequency-dependent loads, such as inductive loads (i.e., motor loads). Hence, the overall characteristics of the change in electrical power due to deviation in the system frequency can be written as follows:
where
is non-frequency-sensitive load change,
is a frequency-sensitive load change, and
is the load damping coefficient.
The load damping constant is defined as the percent of change in load for a 1% change in the system’s frequency. The value of D is typically between 1% and 2%. To illustrate, a value of D = 2 means that the load will change by 2% for a 1% change in the system frequency. Considering the effect of the load damping constant and taking Laplace transform to Equation (12), the transfer function of the relationship between power and variation in speed can be expressed as in Equation (14) and represented in a block diagram as shown in
Figure 5 [
30,
31]:
Analysis of frequency deviation required the use of complex numerical methods because of the nonlinear and time-varying nature of the power system. However, simple first-order transfer functions are used to analyze the power system frequency and control design. The load frequency control (LFC) model, which is widely used in academic power network dynamic studies, was modified to resemble the KSA’s power system [
30,
31].
Figure 6 gives a schematic diagram of the LFC model.
The adopted model is simulated using real data that were collected from various KSA electricity domains, including the SEC and Electricity & Cogeneration Regulatory Authority (ECRA) [
24,
26,
32]. Personal interviews with experts in the SEC have been conducted to ensure the ability of the adopted model to resemble the KSA’s power system. This research analyzes the initial period following power disturbance to evaluate the impact of a reduction in the system’s inertia on the frequency stability. Thus, this model does not include a representation of the secondary frequency control loop (AGC) or other control loops.
Given that all of the generating units in the KSA’s power system are thermal units, all of the generators are substituted with a single generation unit considering their total aggregated installed capacity. The KSA’s synchronous generation is represented by governing turbine transfer functions. The model includes a representation of a speed governor, turbine, rotating mass, and load for the provision of system frequency analysis.
Table 3 shows the values of the parameters used in the simplified LFC model of the KSA’s grid. Generator–turbine time constants (Tg, Ttr, Tt, and Tr) are set according to their typical values. According to the KSA grid code, all generation units must have a governor droop setting of 4% [
2]. Therefore, the droop setting is represented in the test model by the block gain 1/R and is set to 20 p.u. The value of the load damping coefficient (D) is set to 1, which means that a 1% change in system frequency would cause a 1% change in load.
The inertia constant of each generation unit is set according to the generator type to determine the rotational inertia contribution of each generation technology.
Table 4 gives the typical values of the inertia constant for different generation technologies at 3600 rpm [
33,
34]. As can be seen, OCGT and CCGT power plants would deliver most of the inertial response, while non-synchronous generators such as wind and solar PV would not offer any rotational inertial response because they are indirectly connected to the grid using power converters. The overall system inertia (H
eq) is estimated to be 4.15 s through aggregating the inertia constants of all generating units connected directly to the grid.
The accuracy of the adopted model was verified against the KSA’s power system through various real power disturbances. The parameters of the model have been adjusted to ensure that it resembles the dynamic behavior of the real KSA power system.
Figure 7 shows an example of the dynamic behavior of the adopted model when compared to the real KSA system.
The model was validated with a real frequency incident that occurred in the KSA grid when the system frequency dropped to almost 59.78 Hz following a sudden power outage in the eastern operation area, which resulted in a power loss of 1206 MW. The frequency response of the adopted model matched well with the response of the real system for the same power disturbance, which indicates the validity of the adopted model for resembling the KSA power system.