The Impact of Replacing Synchronous Generation with Renewable Energy Technologies on the Transient Stability in the Mangystau Power System

1. Introduction

The problem of climate change causes more and more concern among people every year. The main cause of climate change is greenhouse gas emissions, which have been increasing since the industrial revolution. According to Intergovernmental Panel on Climate Change [1], the bulk of the contribution to global warming has been on CO₂ emissions from fossil fuel combustion and methane. Furthermore, the largest amount greenhouse gas emissions come from electricity and heat sector, accounting for approximately 16 GtCO₂ in 2019 [2].

The global trend in energy development is towards the rapid decommissioning of coal-fired power plants and greater integration of renewable energy sources into electrical networks. Kazakhstan has also been actively involved in the Paris Agreement since 2016, reaching impressive recent achievements and launching new initiatives in the sphere of tacking the climate change [3]. According to Ministry of Energy of the Republic of Kazakhstan [4], at the end of 2022, 130 renewable energy stations operated in the Republic with a cumulative installed capacity of 2400 MW. The volume of electricity generated accounted for 5.11 billion kWh which is 4.53 % of the total electricity production. In the last year only, 12 renewable energy power plants were commissioned with a total installed capacity of 385 MW and in 2023, it is planned to implement an additional 15 facilities with a capacity of 276 MW.

While renewable energy penetration is associated with positive effects on environmental issues, it is fraught with difficulties that need to be addressed. One of those significant technical disadvantages is lack of rotating inertia since they are mainly connected to the grid via power electronic inverters [5]. When, for instance, wind turbines supersede massive rotation synchronous generators operating prevalently in coal plants, it will cause huge danger in terms of grid frequency and rotor angle stability. Makolo et al [6] investigated the impact of inertia on the grid frequency and claimed that increased penetration of low-inertia technologies will lead to a considerable frequency drop and increased rate of change of frequency (RoCoF). One of the examples is a blackout in South Australia in 2016 [7], when the frequency dropped under 47 Hz. The similar emergency occurred in 2018, when Queensland and South Australia systems were eventually separated [8].

The focus of this paper is on the impact of reduced inertia, entailed by high wind energy penetration, on the grid rotor angle transient stability in Mangystau region, Kazakhstan. The loss of the inertial capacity of the system means changes in the parameters in the rotor swing equation [9], and therefore changes in the conditions for the transient stability of parallel operation of power plants and electrically connected systems. The loss of synchronism (the transient stability violation) is identified as the most dangerous operation mode since it is accompanied by deep voltage fluctuations at the nodal points of the electrical network [10,11] which ultimately leads to the shutdown of auxiliary units of thermal and/or nuclear power plants [12].

Significant number of research has been performed in the rotor angle stability problem with penetration of inverter-connected renewable energy. Naik et al [13] investigated the role of inertia in critical clearing time of a transmission line 3-phase short circuit using the potential energy boundary surface method. Liu [14] applied the center of inertia method and described the transient behavior of doubly fed induction generator (DFIG) and its contribution in terms of active and reactive power to the short circuit. He [15] explored the inverter-based generation (IBG) and proposed a sufficient transient stability criterion using a “delta-power-frequency” model which allowed them to discover new mechanisms of dynamic behavior different from those of conventional synchronous generators. The notion that the aggregated inertia constant in high renewable energy system becomes time-variant because of high volatility of the power dispatch was researched and concluded by Ulbig [16]. This, in turn, appears to inevitably amplify the grid instability of power systems.

Collados-Rodriguez [17,18] analyzed transient stability and steady-state stability limits with implementation of voltage source converters (VSC) while reducing the installed capacity of synchronous generators and evaluated the minimum amount of synchronous generation that provides the grid stability. Vittal [19] investigated an interesting relationship between wind generation reactive power control and the rotor angle stability of the synchronous generators in the system. It turns out that voltage and reactive power control is supposed to ease the reactive power burden and voltage issues of synchronous generators. This fact was also proved by [20]. They compared the impact of different DFIGs control techniques on the power systems’ small signal and transient stability and concluded that DFIGs with voltage control mode can enhance the synchronous generators’ rotor angle stability situation compared to those with “power-factor” control mode. Tina [21] investigated potential solutions, such as static var compensator (SVC) and static synchronous compensator (STATCOM), to improve the rotor angle and frequency stability with high DFIG penetration.

All the above studies considered a three-phase short circuit and its subsequent elimination without involving the operation of automatic reclosing. However, in the Unified Power System of the Soviet Union, of which the Power System of Kazakhstan was a part until 1991, the standardization of transient stability was carried out by specifying a list of disturbances under which transient stability should not be violated.

The methodology of this paper applies a 220 kV transmission line’s three-phase short circuit with following operation of protection and tripping the line. An automatic reclose scheme is applied with following line trip again. This scenario is called a “failed automatic reclose” since the line is closed on still existing three-phase short circuit.

The purpose of this work is to:

a)

Study the impact of replacing existing synchronous generators in the Mangystau power system with DFIG generation on the rotor angle transient stability.

b)

Determine the maximum allowable amount of replaced generating capacity at stations with renewable energy technologies.

c)

Explore two methods that can be used to improve the transient stability of the system with a high share of DFIG generation.

2. Methodology for Transient Stability Analysis

It is needed for analysis:

-

replace the dynamic model of generators from a synchronous machine model to the DFIG model shown in Figure 1 (Power World Corporation, [22]). In this case, the inertial constant of the replaced stations becomes equal to zero.

-

the sequence for replacing synchronous generators is Station 2 (250 MW Generation, H=11.5 sec), Station 1 (75 MW Generation, H=1.998 sec) and G-3 Station 3 (20 MW Generation, H=2.604 sec);

-

to analyze the transient stability of the synchronous generators remaining in operation during a three-phase short circuit on the line Station 3 220 kV - Substation 4 220 kV and its unsuccessful automatic reclosure (automatic reclose);

-

transient stability analysis was carried out using PowerWorld software.

In addition, 2 main assumptions are made:

-

the intermittency of power generation from DFIGs is neglected. In other words, all DFIG technologies are equipped with storage systems.

-

the location of DFIG stations is neglected due to the fact that there are no real projects for the implementation of wind power plants of such relatively high power. Replacement of stations on synchronous machines using DFIG technologies is carried out exactly at the location of the replaced stations;

-

the work of automatic excitation control (AEC) and primary frequency control via governors is neglected.

Figure 2 shows the output of active power from DFIG, synchronous generators and equivalent inertia of the Mangystau part of the power system in relation to the share of DFIG technology integration. Thus, with 0 % wind integration, 475 MW is released from stations on synchronous machines and the Unified Power System. Station 2 supplies 250 MW and when it is replaced by DFIG, it turns out that 250 MW is generated from DFIG, and 225 MW from synchronous machines (53 % of generation from renewable energy sources). This reduces the equivalent inertial constant since Station 2 had an inertia of 11.5 seconds. At 73 % of generation from renewable energy sources, 130 MW is supplied from synchronous machines and the Unified Power System (20 MW each from G-1 and G-2 thermal power plants and 90 MW from the Unified Power System).

The general philosophy of transient stability during a short circuit is to ensure that the decelerating area in the angular characteristic of the generator is greater than the acceleration area. Let us consider the behavior of the angular characteristics of one generator operating on a large system with single machine infinite bus (SMIB). Figure 3 shows the circuits of a synchronous generator connected to the SMIB in pre-emergency, emergency, and post-emergency mode. Figure 4 shows the angular characteristic of this generator.

The angular characteristic is determined by the formula:

$$P=\frac{E^{\prime }\cdot V}{X}\sin \delta$$

(1)

where $E^{\prime }$ – is the value of the electromotive force (EMF) of the generator in question; V – is the voltage value on the receiving bus of infinite power; X – is the equivalent resistance between generator and SMIB; δ – is the difference between the angle of the EMF of the generator and the voltage on the SMIB.

Thus, the angular characteristic (I) corresponds to the pre-emergency mode, when all power lines are turned on and no accident occurs. Angle δ from formula (1) corresponds to point 1 in Figure 4. Consequently, the value of the difference in the angle between the generator and the SMIB in pre-emergency mode is equal to δ₁. Further, when a three-phase short circuit occurs, due to a considerable voltage drop and power output decrease, the circuit and angular characteristic take the form corresponding to (II). At the same time, the electrical power supplied by the generator is significantly reduced, and the turbine power has not yet had time to decrease during the transient process. As a result, the turbine begins to accelerate the generator rotor, which leads to an increase in its relative angle. The increase in the angle value continues until the moment when the short circuit is eliminated by the action of the protection. At the moment when the damaged line is disconnected, the circuit and angular characteristic take the form corresponding to (III), since one line was disconnected and the equivalent resistance increased, which led to a decrease in the transmitted power. At point 4, as can be seen from the figure, the electrical power exceeds the turbine power. This, in turn, leads to sharp decelerating of the generator rotor. However, the rotor angle still increases up to point 5. At this moment, automatic re-closing occurs for an unresolved short circuit. Process 1–5 is repeated at points 5–9. Point 9 is the moment when all the stored kinetic energy has dissipated, and the generator is pulled into synchronous operation with the power system. The rotor angle goes to the value indicated at point 1’. However, if point 9 crossed the critical value of the rotor deflection angle (δ_crit), then the synchronous machine would fall out of synchronization with the power system and an asynchronous mode would arise. The condition for maintaining transient stability is that the total decelerating areas (S3 and S4) be greater than the total acceleration areas (S1 and S2). The formula for meeting this condition is given below:

$${\displaystyle {\int }_{{\delta }_1}^{{\delta }_3}\left(P_{mech}-P_{2\max }\sin \delta \right)d\delta +}{\displaystyle {\int }_{{\delta }_5}^{{\delta }_7}\left(\left(P_{mech}-P_{2\max }\sin \delta \right)d\delta \right)}\le {\displaystyle {\int }_{{\delta }_3}^{{\delta }_5}\left(P_{3\max }\sin \delta -P_{mech}\right)d\delta +}{\displaystyle {\int }_{{\delta }_7}^{{\delta }_9}\left(\left(P_{3\max }\sin \delta -P_{mech}\right)d\delta \right)}$$

(2)

In the case when the generator is very close to the short circuit, its supplied electrical power is equal to zero. From the equation of rotor motion [23], it is possible to derive the value of the maximum deviation of the rotor angle with a known and constant tripping time of the damaged section for the scenario when automatic reclose is not used:

$${\delta }_{\max }=\frac{\pi f_0}{2H}{P}_{mech}{t}^2+{\delta }_0$$

(3)

where f₀ – is the nominal network frequency; $H$ is the generator inertia constant; $P_{mech}$ is the turbine power (and generator electrical power in pre-emergency mode); t – is the time to turn off the damaged area.

From Eq. (3) it is clear that with a constant relative angle and power output in the pre-emergency mode, as well as a constant short circuit tripping time, the value of the maximum deviation of the relative angle depends inversely on the inertial constant. In other words, the lower the inertial constant of the generator/system, the higher the maximum angle deviation during the transient process will be. However, in this work, it is not possible to calculate relative angles using Eq. (3) due to the fact that the system has several synchronous machines, and the emergency operation process is very complex. To calculate the stability of multi-machine systems with complex transient processes, the method of numerical integration of the rotor motion equation is used (Glover et al., 2015; Saadat, 1999). Modern software such as DigSilent PowerFactory, ETAP, PSSE, PowerWorld Simulator etc., allow to carry out stability calculations in systems of any complexity.

To assess transient stability, the Transient Stability Index is used, determined by the formula (Shi et al., 2009):

$$TSI=\frac{{360}^{\circ }-{\rho }_{\max }}{{360}^{\circ }+{\rho }_{\max }}\cdot 100$$

(4)

where ${\rho }_{\max }$ – is the maximum deviation of the rotor angle relative to the base machine (equivalent to the Kazakhstan system on Substation 4). A system is considered dynamically stable if the transient stability index is greater than zero.

3. Materials and Methods

The goal is achieved by solving the following tasks:

1)

Model a real diagram of the Mangystau power system.

2)

Linking the inertial constant to the nominal parameters of existing synchronous generators of the Mangistau energy system.

3)

Set fault: three-phase short circuit on the 220 kV intersystem transmission line with subsequent tripping action by relay protection and automatic reclosure action. In this case, the line is switched on for an unresolved short circuit, and the line is switched off again by the action of the protection.

4)

Identify the maximum possible share of DFIG generation that can replace the synchronous generators.

5)

Identify the influence of the automatic reclosure response time delay on a faulted line, as a means of enhancing transient stability of the power system with high DFIG integration.

Identify the influence of the magnitude of the load shedding as a means of improving transient stability under conditions of high DFIG integration.

This section may be divided by subheadings. It should provide a concise and precise description of the experimental results, their interpretation, as well as the experimental conclusion that can be drawn.

3.1. Configuration of the Mangystau Energy System

Main components of the system (see Figure 5):

-

the electrical power system under study consists of three main power plants: Station 1, 2 and 3.

-

there are Open Switchgears 110 kV at Station 1 and 2 while Station 3 has both Open Switchgear 110 and 220 kV.

-

the given transmission network is connected to the rest of the Unified Power System of Kazakhstan through Lines 1 and 2 with voltage of 220 kV. All the generators in Stations 1, 2 and 3 are rotating in synchronism with the aggregated Unified Power System. Substation 4 is the slack bus.

Figure 5 shows a single-line diagram of the Mangystau power system. This diagram was built based on the map of the unified power system of Kazakhstan, which is publicly available on the website of the system operator (Kazakhstan Electricity Grid Operating Company, [25]). In addition, all rating data, resistance and dynamic characteristics of generators, transformers and power lines (power lines) were taken from the reference book for designing electrical power systems.

Composition of turbogenerators:

Table 1 shows the composition of turbogenerators for Station 1, 2 and 3 with their nominal active capacity. All technical characteristics such as power factor, transient and synchronous reactance, leakage reactance, inertia, and rpm (rotation per minute) have been taken from the power systems design handbooks.

3.2. Aggregating and Calculation of Network Parameters of the Mangystau System

Equivalent inertia is the sum of all inertial constants reduced to one common basis. The resistances of all generators are collected in parallel into one equivalent. The details of dynamic aggregation have been described in detail in (Germond and Podmore, [26]). Thus, in this work, the Station 1 and Station 2 generators are presented as one equivalent in order to simplify the circuit. At Station 3, each generator block is modeled separately, since they are the most powerful and largest, have the largest inertial constants, and are located in close proximity to the short circuit point.

A. 

Calculation of inertia of turbogenerators

The moment of inertia of synchronous generators is given in units - T∙m² and is denoted by the letter J. However, for this analysis it is necessary to convert the moment of inertia into an inertial constant (measuring unit – seconds):

$$H=\frac{W_k}{S_{nom}}\left(c\right)$$

(5)

where $W_k$ – is the stored kinetic energy of rotating motion (Joule); $S_{nom}$ – is the rated apparent power of the generator (Volt∙Аmps).

The kinetic energy of rotating motion is determined by the generally accepted formula:

$$W_k=\frac{1}{2}J{\omega }^2$$

(6)

where $\omega =2\pi f$ – is the angular velocity; $f$ – is the nominal network frequency, 50 Hz.

In order to determine the dependence of maintaining transient stability in relation to the rated power of the system, it is necessary to bring all the inertial constants of all turbogenerators under one basis. The basis is taken equal to Sb= 100 МVА, which allows all calculations to be carried out in relative units:

$$H^{Sb=100MVA}=H^{Snom}\cdot \frac{S_{nom}}{S_b}\left(c\right)$$

(7)

The equivalent inertial constant is obtained by adding all the inertias of the turbogenerators connected to the power system:

$${\displaystyle \sum H={\displaystyle \sum H_i}}\left(c\right)$$

(8)

where i – is the turbogenerator index.

B. 

Calculation of system element resistances

Generator reactance

The values of transition and synchronous resistance in relative units are taken from the reference data of the generator. Transient reactance is taken to analyze the transient processes. Synchronous resistance is used to calculate normal conditions and steady-state short-circuit current after saturation of the generator magnetic system. For each of these parameters, resistances are calculated using the formulas below:

$$x_g={x^{\prime }}_d\left(x_d\right)\cdot \frac{S_b}{S_{nom}}$$

(9)

Transformer reactance:

$$x_t=\frac{x_t\%}{100}\cdot \frac{S_b}{S_{nom}}$$

(10)

The $x_t\%$ value is calculated from the short circuit voltage (%), which is indicated in the transformer data sheet. For a two-winding transformer:

$$x_t\%=u_K\%$$

(11)

For a three-winding transformer:

$$\begin{array}{c}{x}_{tB}\%=0.5\left(u_{KB-H}+u_{KB-C}-u_{KC-H}\right)\\ x_{tC}\%=0.5\left(u_{KB-C}+u_{KC-H}-u_{KB-H}\right)\\ x_{tH}\%=0.5\left(u_{KB-H}+u_{KC-H}-u_{KB-C}\right)\end{array}$$

(12)

Power line inductive resistance:

$$x_l=x_{sp}l\frac{S_b}{U_b^2}$$

(13)

where $x_{sp}$ – is the wire specific reactance (Ohm/km); l – is the line length (km); S_b – is the base power (MVA); U_b – is the base voltage, which is selected for a specific section of the network (kV).

The equivalent resistance of the system is calculated by connecting all generators connected to the network in parallel, connecting transformers and power lines in series and parallel, and converting the equivalent star to delta and vice versa. Equivalent circuit conversion formulas are presented in Glover et al (2015).

Table 2 shows the results of calculations of equivalent generator resistances and inertial constants. Inertial constants are given under a common basis. When decommissioning a specific generator unit(s), their inertial constants will be subtracted from the total. Table 3 and Table 4 present the electrical parameters of connecting power lines and transformers.

3.3. Setting the Pre-Emergency Mode

The task of the pre-emergency mode is to solve the problem of power flows throughout the entire circuit. To do this, it is necessary to set load on the nodes, active power of the generators and terminal voltage, and determine the slack bus. Substation 4 220 kV is considered as slack bus with voltage of 1 per unit and angle equal to 0 degrees. Table 5 and Table 6 show the operating data of loads and generators.

4. Results

Simulations of changes in the electrical and electromechanical parameters of the system during the transient process were performed in accordance with the methodology described in Chapter 4. The main system parameters for analysis in this chapter are the rotor angle of the generators in the Mangystau system relative to the Unified Energy System, voltage levels at each circuit node and electrical frequency on the nodes. Below are descriptions of each scenario:

1)

0 % generation from DFIG - the base scenario in which the generators of each station are synchronous. Equivalent inertia 22.26 sec;

2)

53 % of generation from DFIG is the replacement of synchronous generators with Station 2. In this scenario, Station 2 is removed from the analysis due to the fact that it is no longer a synchronous generator, but a DFIG connected through a power inverter. Equivalent inertia 10.76 sec;

3)

68 % of generation from DFIG - replacement of synchronous generators at Station 1. Remaining synchronous machines - G-1, G-2 and G-3 at Station 3. Equivalent inertia 8, 76 sec;

4)

73 % of generation from DFIG - replacement of one synchronous generator with Station 3 (G-3). In this scenario, the system is supported only by two Station 3 synchronous generators (G-1 and G-2). Equivalent inertia 6.16 sec;

5)

73 % generation from DFIG and an increase in the automatic reclose time by 0.5 sec - a scenario similar to 4);

6)

73 % generation from DFIG and disconnection of 20 MW load on Substation 2 - a scenario similar to 4), but with part of the load disconnection from special load shedding automation (SLSA) as a means of preventing instability.

Scenarios 1-3

From Figure 6, Figure 7 and Figure 8 the transient stability is maintained up to a generation fraction of DFIG of 68 %. The deviation of the relative angles of the rotors of synchronous generators is maintained within the permissible values so that the generators do not fall out of synchronism. Nevertheless, it is seen that the steady relative rotor angles changed due to the 220 kV line eventually tripped, resulting in a change in the circuit as shown in Figure 3. The relative angle of the synchronous generators took the value 1' as described in Figure 4.

The voltages at each node experience a transient mode during a short circuit, dropping significantly twice due to the occurrence of a short circuit and turning the line back on for an unresolved short circuit. However, after the second disconnection of the line from the protections, the voltage levels return to their original state.

The preservation of synchronous operation is also visible from the graphs of frequency changes in each node of the power system. In all three scenarios, the frequency undergoes a transient process, changing in value, but eventually returns to synchronous operation, where the frequency is the same at all nodes. However, an increase in frequency is observed, since this work neglects the work of Primary Frequency Control and Automatic Excitation Regulation, which was mentioned in the methodology chapter.

Scenario 4

At 73 % generation from DFIG, the Mangystau system becomes dynamically unstable to a three-phase short circuit with a failed automatic reclose. From Figure 9 it is clear that when such a severe emergency occurs, the synchronous generators G-1 and G-2 remaining in operation are physically unable to maintain stability in the system and fall out of synchronism. Since these generators operate at the same power plant, they enter an asynchronous mode with one electrical swing center (ESC), which is located at the intersystem transmission line 220 kV. The relative angle of the rotors exceeds a critical value, and then it is no longer possible to return to the original state.

Voltage levels begin to fluctuate and eventually a deep drop in voltage occurs, making it unsuitable for consumers and power plants to operate. It should be noted that the ESC is located closer to the 220 kV Substation No. 4 buses, since a deeper voltage reduction occurs there at the initial stage of the asynchronous mode.

From the frequency change graph, the frequency at substation No. 4 increases, while at other nodes of the Mangystau system it decreases. This indicates that the system is indeed asynchronous. In addition, we can conclude that the Mangystau system in this mode is deficient, while the rest of the Unified Energy System is excess. From the technical literature, in the deficient part of the power system, special load shedding automation (SLSA) must be provided as a means of increasing stability (Berkovich M., 1991), while in the excess part, generation must be limited by influencing the speed controllers of turbines at power plants.

It is worth noting that even if the transient stability is maintained at 53 % and 68 % wind generation, the transient stability index gradually weakens as the equivalent inertial constant decreases. This can be seen from Figure 10, where at 0 % DFIG share, all 5 synchronous machines shown have a high level of transient stability index (above 70 %). However, with the gradual replacement of generators with DFIG technologies, the remaining synchronous machines in operation lose a certain percentage of the stability index. This is because with reduced equivalent inertia, the maximum rotor deflection angle increases (see Equation 3). So, for example, in scenario 4), when the two synchronous machines remaining in operation (G-1 and G-2) fall out of synchronism, the relative rotor angle exceeds the critical value and goes to infinity until the emergency automation operates. In this regard, the stability index of these generators reaches a negative value of almost 100 %.

Scenarios 5-6

If a short circuit occurs and the line is successfully disconnected from the relay protection, the system circuit changes. However, the operating mode of the power system returns to the steady state, since the damaged element is taken out of operation. From the graph of voltage changes in asynchronous mode (see Figure 9b) it can be seen that before the automatic reclose operation, the voltage levels tend to return to their original value. However, the line is then switched on with unresolved short circuit. This leads to the fact that the power system cannot withstand such a strong disturbance and falls out of synchronization.

If the delay time of automatic reclose operation is increased at the same fraction of DFIG, the power system will have more time to return to normal steady state before restarting. Thus, the increased reclosing operating time has a positive effect on transient stability.

Thus, in Figure 11b and Figure 12b we can see a comparison of the transient process of the relative rotor angles and voltage levels at the nodes at 73 % generation from DFIG and taking measures to increase the operating time of automatic reclose. An increase in the automatic reclose time delay by 0.5s had a positive effect on the nature of the transient process and allowed the two remaining synchronous generators in the system to maintain transient stability and synchronism with the Unified Energy System.

Figure 13 describes the dependence of the transient stability index and the operating time of the automatic restart. Initially, the operating time was 3 sec after the line trip. By increasing the time delay by 0.5 sec, the stability index increases to 40 %.

As for load shedding, this measure reduces the power supplied by synchronous generators by an amount equivalent to the load being switched off. From Eq. (3), the lower the released active power of the generator, the smaller the maximum rotor deflection angle relative to a large power system. Thus, special automatic load shedding can improve transient stability in the power system.

From Figure 11c and 12c it is clear that disconnecting the load at substation No. 2 in the amount of 20 MW had a positive effect on the nature of the transient process and made it possible to maintain synchronism.

Figure 14 illustrates the dependence of the transient stability index in relation to the magnitude of the disconnected load. Thus, starting from the switched load value of 20 MW, the dynamic stability index sharply increases to 50 %, which means that the condition for maintaining dynamic stability in the power system is met.

5. Conclusions

• In this steady state regime, the share of energy generation from wind units should not exceed 68 % or 325 MW out of a total generation of 475 MW (of which 90 MW is the flow from the unified power system) in the Mangystau power system.
• The minimum output from synchronous generators should not be less than 60 MW out of 475 MW of total generating capacity, taking into account the flow through 220 kV intersystem lines.
• The minimum number of synchronous generators to ensure transient stability during the disturbance under study is 3 synchronous generators at “Station No. 3” with a rated power of 200 MW and a total inertial constant of 8.76 sec.
• When the third synchronous generator at “Station No. 3” (G-3) is taken out of operation and replaced with a wind station, transient stability is violated due to a critical lack of inertia.
• The means of increasing dynamic stability are increasing the response time of automatic reclose and/or disconnecting part of the load (special load shedding automation). Figure 11 and Figure 12 show that at 73 % wind generation (or 345 MW), both methods lead to the prevention of transient stability violation.

This section is not mandatory but may be added if there are patents resulting from the work reported in this manuscript.