Feasibility Study of Wind Farm Grid-Connected Project in Algeria under Grid Fault Conditions Using D-Facts Devices

The use of renewable energy such as wind power is one of the most affordable solutions to meet the basic demand for electricity because it is the cleanest and most efficient resource. In Algeria, the highland region has considerable wind potential. However, the electrical power system located is this region is generally not powerful enough to solve the problems of voltage instability during grid fault conditions. These problems can make the connection with the eventual installation of a wind farm very difficult and inefficient. Therefore, a wind farm project in this region may require dynamic compensation devices, such as a distributed-flexible AC transmission system (D-FACTS) to improve its fault ride through (FRT) capability. This paper investigates the implementation of shunt D-FACTS, under grid fault conditions, considering the grid requirements over FRT performance and the voltage stability issue for a wind farm connected to the distribution network in the Algerian highland region. Two types of D-FACTSs considered in this paper are the distribution static VAr compensator (D-SVC) and the distribution static synchronous compensator (D-STATCOM). Some simulation results show a comparative study between the D-SVC and D-STATCOM devices connected at the point of common coupling (PCC) to support a wind farm based on a doubly fed induction generator (DFIG) under grid fault conditions. Finally, an appropriate solution to this problem is presented by sizing and giving the suitable choice of D-FACTS, while offering a feasibility study of this wind farm project by economic analysis.


Introduction
In Algeria, the first attempt to connect the wind energy conversion system (WECS) to the electricity distribution network dates back to 1957, with the installation of a 100-kW wind turbine at the Grands Vents site (Algiers) by the French designer Andreau.
Nowadays, the depletion of fossil fuels reserves in Algeria, fluctuations in oil price and the location of energy resources are causing instability in energy policy.
In addition, the use of fossil fuels for conventional power plants triggers alarms of an environmental disaster.Currently, to reduce the harmful impact of conventional resources and improve Algerian energy efficiency, the energy policy program announced by the Ministry of Mines and Energy aims, by 2030, to produce 40% electrical energy from renewable resources [1].For WECS, the power to be produced over the period 2012-2022 is estimated at approximately 516 MW, of which 10 MW are installed at Kabertene (70 km from Adrar) in the Algerian desert [1][2][3][4].This pilot wind farm consists of 12 wind turbines with a unit capacity of 0.85 MW, the energy produced will be injected to the 30/220 kV step up transformer situated in the same locality [1], as shown in Figure 1.Currently, the Algerian electrical grid code does not consider WECS.In the region of Adrar, the electrical grid is not interconnected with the north; it is a local grid (or micro-grid).Therefore, this program of the energy policy must be accompanied by continual development of wind energy technology and optimization techniques, looking for better options concerning reduced costs, improvement regarding wind turbine performances, reliability of electrical groups and electrical grid integration.
Appl.Sci.2018, 8, x 2 of 22 In addition, the use of fossil fuels for conventional power plants triggers alarms of an environmental disaster.Currently, to reduce the harmful impact of conventional resources and improve Algerian energy efficiency, the energy policy program announced by the Ministry of Mines and Energy aims, by 2030, to produce 40% electrical energy from renewable resources [1].For WECS, the power to be produced over the period 2012-2022 is estimated at approximately 516 MW, of which 10 MW are installed at Kabertene (70 km from Adrar) in the Algerian desert [1][2][3][4].This pilot wind farm consists of 12 wind turbines with a unit capacity of 0.85 MW, the energy produced will be injected to the 30/220 kV step up transformer situated in the same locality [1], as shown in Figure 1.Currently, the Algerian electrical grid code does not consider WECS.In the region of Adrar, the electrical grid is not interconnected with the north; it is a local grid (or micro-grid).Therefore, this program of the energy policy must be accompanied by continual development of wind energy technology and optimization techniques, looking for better options concerning reduced costs, improvement regarding wind turbine performances, reliability of electrical groups and electrical grid integration.In the period of 2009-2010, Sebaa Ben Miloud F et al. [6] and Himri et al. [7] undertook the first study to identify a suitable site in Adrar region for the wind farm installation.In addition, Himri et al. [7] used data of wind speed over a period of nearly 10 years to assess the potential of wind power stations in two southern Algerian regions, namely, Timimoun and Tindouf.In [8], a study of the wind potential in seven southern Algerian sites was undertaken, from west to east, Tindouf, Bechar, Adrar, and Ghardaia, In Amenas and In Salah (Tamanrasset).In [9] wind speed data was collected over a period of almost 5 years, from three selected stations in northern Algeria.Within this context, some studies in the Algerian high plateau region were performed in [10][11][12].However, the authors did not consider the integration issue of the wind farm into the electrical grid and it is well known that the electrical grid influences greatly the performance of wind farm installation and production.
At the present time, doubly fed induction generators (DFIG) are the most used in WECS [13,14] and especially in the Algerian wind farm at Kabertene in Adrar.Simple induction generators have some weaknesses such as reactive power absorption and uncontrolled voltage during variable rotor speed.These complications are avoided by the installation of DFIG and power converters or power drives [15,16].The particular feature of the DFIG is that the injected power by the rotor converter is only a small part from the total provided power with its stator directly connected to the electrical grid [17][18][19].Hence, the size, the cost and losses of the power converter are optimized compared to a full-size power converter of the other generators.
One of the most important considerations in a wind farm grid-connected project is fault ridethrough (FRT) capability, where the energy grid is often weak and the DFIG is frequently working In the period of 2009-2010, Sebaa Ben Miloud F et al. [6] and Himri et al. [7] undertook the first study to identify a suitable site in Adrar region for the wind farm installation.In addition, Himri et al. [7] used data of wind speed over a period of nearly 10 years to assess the potential of wind power stations in two southern Algerian regions, namely, Timimoun and Tindouf.In [8], a study of the wind potential in seven southern Algerian sites was undertaken, from west to east, Tindouf, Bechar, Adrar, and Ghardaia, In Amenas and In Salah (Tamanrasset).In [9] wind speed data was collected over a period of almost 5 years, from three selected stations in northern Algeria.Within this context, some studies in the Algerian high plateau region were performed in [10][11][12].However, the authors did not consider the integration issue of the wind farm into the electrical grid and it is well known that the electrical grid influences greatly the performance of wind farm installation and production.
At the present time, doubly fed induction generators (DFIG) are the most used in WECS [13,14] and especially in the Algerian wind farm at Kabertene in Adrar.Simple induction generators have some weaknesses such as reactive power absorption and uncontrolled voltage during variable rotor speed.These complications are avoided by the installation of DFIG and power converters or power drives [15,16].The particular feature of the DFIG is that the injected power by the rotor converter is only a small part from the total provided power with its stator directly connected to the electrical grid [17][18][19].Hence, the size, the cost and losses of the power converter are optimized compared to a full-size power converter of the other generators.
One of the most important considerations in a wind farm grid-connected project is fault ride-through (FRT) capability, where the energy grid is often weak and the DFIG is frequently working under grid faults when the wind farm is located relatively far from this electrical grid [20].Therefore, many research works focus on studying the dynamic behaviors of wind farms during and after the clearance of the grid fault conditions without disconnection from the electrical grid.In [21], several methods employed to improve the FRT capability of the fixed-speed wind turbines are based on induction generators.In [20], an enhanced application to overcome grid fault conditions is studied for a wind farm based on DFIG.FRT control of wind turbines with DFIGs under symmetrical voltage dips is presented in [21].In [22] a flexible AC transmission system (FACTS) system for DFIG to reduce the effects of grid faults is proposed.
The present paper can extend the aforementioned research works.This paper, shows the feasibility of installing a wind farm in an Algerian highland region, confirmed by some data of wind potential in the selected geographical location in the first part, which is an input for the wind power system.On the other hand, another important aspect is the weakness of the electrical grid, which is often an obstacle in many countries that have established wind energy projects.Consequently, during the feasibility study, some techniques to identify the cost-effectiveness of areas for the wind farms installations, the possible electrical path of distribution lines, and their corresponding estimated cost are used, and the incorporation of electrical devices such as distributed FACTS (D-FACTS) technology is considered.Furthermore, the investors may be confident to fund this possible project when these technical difficulties are taken into consideration.Then, in order to ensure the economic success of the future wind farm project in the highland region, an accurate study by some simulation results, showing the interaction between wind turbine generators and the electrical grid in this region with the impact of the D-FACTS systems, is undertaken, which has not been previously done.In this study, the Algerian electrical grid code could be considered in simulations similar to that of the Spanish grid code as shown in Figure 2.
Appl.Sci.2018, 8, x 3 of 22 under grid faults when the wind farm is located relatively far from this electrical grid [20].Therefore, many research works focus on studying the dynamic behaviors of wind farms during and after the clearance of the grid fault conditions without disconnection from the electrical grid.In [21], several methods employed to improve the FRT capability of the fixed-speed wind turbines are based on induction generators.In [20], an enhanced application to overcome grid fault conditions is studied for a wind farm based on DFIG.FRT control of wind turbines with DFIGs under symmetrical voltage dips is presented in [21].In [22] a flexible AC transmission system (FACTS) system for DFIG to reduce the effects of grid faults is proposed.The present paper can extend the aforementioned research works.This paper, shows the feasibility of installing a wind farm in an Algerian highland region, confirmed by some data of wind potential in the selected geographical location in the first part, which is an input for the wind power system.On the other hand, another important aspect is the weakness of the electrical grid, which is often an obstacle in many countries that have established wind energy projects.Consequently, during the feasibility study, some techniques to identify the cost-effectiveness of areas for the wind farms installations, the possible electrical path of distribution lines, and their corresponding estimated cost are used, and the incorporation of electrical devices such as distributed FACTS (D-FACTS) technology is considered.Furthermore, the investors may be confident to fund this possible project when these technical difficulties are taken into consideration.Then, in order to ensure the economic success of the future wind farm project in the highland region, an accurate study by some simulation results, showing the interaction between wind turbine generators and the electrical grid in this region with the impact of the D-FACTS systems, is undertaken, which has not been previously done.In this study, the Algerian electrical grid code could be considered in simulations similar to that of the Spanish grid code as shown in Figure 2.

Algerian Wind Potential
This section discusses a method for determining the production of wind energy at different sites in Algeria in order to choose a suitable site for a cost-effective energy installation.Thus, we have both the average wind speed and the power produced by wind turbines; we can combine them to calculate the energy produced by these wind turbines.Furthermore, in this paper, five selected geographical locations (altitude, latitude and longitude) shown in Table 1 and Figure 3 were obtained from the National Meteorological Office (NMO) [24].

Algerian Wind Potential
This section discusses a method for determining the production of wind energy at different sites in Algeria in order to choose a suitable site for a cost-effective energy installation.Thus, we have both the average wind speed and the power produced by wind turbines; we can combine them to calculate the energy produced by these wind turbines.Furthermore, in this paper, five selected geographical locations (altitude, latitude and longitude) shown in Table 1 and Figure 3 were obtained from the National Meteorological Office (NMO) [24].These wind speed data are collected only at 10 m of altitude, measured using a type of anemometer cup and vane.However, the action of the wind speed at the turbine (tower height over 70 m) is very complex, and includes both deterministic effects (wind and shadow average round), and stochastic fast varying wind speed is turbulent.In fact, wind speed describing these variations is usually measured in the lower atmosphere using either instrumented towers or tethered balloons, which have not been available in previous stations [25,26].
The wind speed is the most important aspect of wind potential; in fact, the annual variation of the long-term average wind speed provides a good understanding of the long-term trend of wind speed and gives confidence to investors on the availability of wind energy in the years ahead [27].Figure 4 provides the average wind speed during five years of data collection at 5 stations in Algeria, which are considered in this study [7,34].These wind speed data are collected only at 10 m of altitude, measured using a type of anemometer cup and vane.However, the action of the wind speed at the turbine (tower height over 70 m) is very complex, and includes both deterministic effects (wind and shadow average round), and stochastic fast varying wind speed is turbulent.In fact, wind speed describing these variations is usually measured in the lower atmosphere using either instrumented towers or tethered balloons, which have not been available in previous stations [25,26].
The wind speed is the most important aspect of wind potential; in fact, the annual variation of the long-term average wind speed provides a good understanding of the long-term trend of wind speed and gives confidence to investors on the availability of wind energy in the years ahead [27].Figure 4 provides the average wind speed during five years of data collection at 5 stations in Algeria, which are considered in this study [7,.These wind speed data are collected only at 10 m of altitude, measured using a type of anemometer cup and vane.However, the action of the wind speed at the turbine (tower height over 70 m) is very complex, and includes both deterministic effects (wind and shadow average round), and stochastic fast varying wind speed is turbulent.In fact, wind speed describing these variations is usually measured in the lower atmosphere using either instrumented towers or tethered balloons, which have not been available in previous stations [25,26].
The wind speed is the most important aspect of wind potential; in fact, the annual variation of the long-term average wind speed provides a good understanding of the long-term trend of wind speed and gives confidence to investors on the availability of wind energy in the years ahead [27].Figure 4 provides the average wind speed during five years of data collection at 5 stations in Algeria, which are considered in this study [7,34].

Connection Issue of a Wind Farm to the Electrical Grid
In this paper, the proposed wind farm in a highland region with average power is considered as a decentralized generator unit, which is most often connected to the distribution network and that differs from centralized generator units.In Algeria, the electrical distribution grids are the most important infrastructure of the whole power system, which is considered as the final interface that leads to most industrial and domestic customers.These distribution grids are operated in ranges of voltages below 50 kV, which is the voltage level of the Medium Voltage (MV) and Low Voltage (LV) ranges.Moreover, in the Algerian distribution grid, the nominal voltage of the MV is 10 kV and 30 kV.These voltage levels allow a good compromise to limit the voltage drops, minimizing the number of source positions (connecting to High Voltage (HV)/MV power station) and reduce the inherent constraints to high voltages (investment costs, protection of property and persons).Moreover, Algerian distribution grids are, in most cases, radially networked.The map of the western Algerian electrical grid is shown in Figure 5.This figure shows the structure of the High Voltage B 220 kV transmission lines, while the substations and power plants are also shown in this figure [1].The structure of High Voltage A 60 kV distribution lines is shown in Figure 6.

Connection Issue of a Wind Farm to the Electrical Grid
In this paper, the proposed wind farm in a highland region with average power is considered as a decentralized generator unit, which is most often connected to the distribution network and that differs from centralized generator units.In Algeria, the electrical distribution grids are the most important infrastructure of the whole power system, which is considered as the final interface that leads to most industrial and domestic customers.These distribution grids are operated in ranges of voltages below 50 kV, which is the voltage level of the Medium Voltage (MV) and Low Voltage (LV) ranges.Moreover, in the Algerian distribution grid, the nominal voltage of the MV is 10 kV and 30 kV.These voltage levels allow a good compromise to limit the voltage drops, minimizing the number of source positions (connecting to High Voltage (HV)/MV power station) and reduce the inherent constraints to high voltages (investment costs, protection of property and persons).Moreover, Algerian distribution grids are, in most cases, radially networked.The map of the western Algerian electrical grid is shown in Figure 5.This figure shows the structure of the High Voltage B 220 kV transmission lines, while the substations and power plants are also shown in this figure [1].The structure of High Voltage A 60 kV distribution lines is shown in Figure 6.Based on the structure of the electrical grid on the west side near the proposed wind farm (see Figure 6), this radial electrical grid of 60 kV can be reconfigured; it is then simulated using the Power System Analysis Toolbox (PSAT) software with actual electrical grid parameters and consumer profiles at the peak load of each bus, with a centralized generation source Tiaret City (TIARC) power plant.Simulation of the latter gives the results shown in Figures 7 and 8.More details on the overall simulation of the west Algerian grid can be found in [34].Based on the structure of the electrical grid on the west side near the proposed wind farm (see Figure 6), this radial electrical grid of 60 kV can be reconfigured; it is then simulated using the Power System Analysis Toolbox (PSAT) software with actual electrical grid parameters and consumer profiles at the peak load of each bus, with a centralized generation source Tiaret City (TIARC) power plant.Simulation of the latter gives the results shown in Figures 7 and 8.More details on the overall simulation of the west Algerian grid can be found in [34].Based on the structure of the electrical grid on the west side near the proposed wind farm (see Figure 6), this radial electrical grid of 60 kV can be reconfigured; it is then simulated using the Power System Analysis Toolbox (PSAT) software with actual electrical grid parameters and consumer profiles at the peak load of each bus, with a centralized generation source Tiaret City (TIARC) power plant.Simulation of the latter gives the results shown in Figures 7 and 8.More details on the overall simulation of the west Algerian grid can be found in [34].Based on the structure of the electrical grid on the west side near the proposed wind farm (see Figure 6), this radial electrical grid of 60 kV can be reconfigured; it is then simulated using the Power System Analysis Toolbox (PSAT) software with actual electrical grid parameters and consumer profiles at the peak load of each bus, with a centralized generation source Tiaret City (TIARC) power plant.Simulation of the latter gives the results shown in Figures 7 and 8.More details on the overall simulation of the west Algerian grid can be found in [34].The existence of an SNVI industrial site in the bus of N 2 between the city of Tiaret and the city of Tissemsilt can justify the presence of the voltage drop across this line as shown in Figure 7a.This resulted in demand of an excessive reactive power and active line losses due to the high-fluctuated demand of energy to the industrial site, as shown in Figures 7 and 8.

Distribution Flexible Alternative Current Transmission System (D-FACTS)
Shunt D-FACTS devices can be classified into two main categories, namely the variable impedance type such as the distribution static var compensator (D-SVC) and the switching converter type such as the distribution static synchronous compensator (D-STATCOM).
The configuration of the D-SVC connected to the distribution grid is shown in Figure 9.
Appl.Sci.2018, 8, x 7 of 22 The existence of an SNVI industrial site in the bus of N 2 between the city of Tiaret and the city of Tissemsilt can justify the presence of the voltage drop across this line as shown in Figure 7a.This resulted in demand of an excessive reactive power and active line losses due to the high-fluctuated demand of energy to the industrial site, as shown in Figures 7b and 8.

Distribution Flexible Alternative Current Transmission System (D-FACTS)
Shunt D-FACTS devices can be classified into two main categories, namely the variable impedance type such as the distribution static var compensator (D-SVC) and the switching converter type such as the distribution static synchronous compensator (D-STATCOM).
The configuration of the D-SVC connected to the distribution grid is shown in Figure 9.

Constitution of the Wind Farm and the Location of the Shunt D-FACTS
The wind farm connected to the distribution system proposed in this section is shown in Figure 11, which consists of eight DFIGs with 1.5 MW of power for each wind turbine.These generators are connected between them to a voltage level of 30 kV by a step-up transformer 690 V/30 kV with 4 MVA of power for each generator.Then a 45 km line that is connected to the source substation 60 kV through another step-up transformer 30 kV/60 kV with 47 MVA of power.For this study, these lines are modelled with the π model.
Based on the work done in [21,35] the simulation results obtained with a D-FACTS provides an effective support to the bus voltage to which it is connected.Therefore, in this study, D-FACTS is placed at the PCC for two reasons:

•
The location for the reactive power support should be as close as possible to the point at which the carrier is necessary because of the variation in the voltage and, therefore, power loss (Joule loss) in the distribution line associated with reactive power flow, • In the studied system, the effect of the change in voltage is most common in this bus.

Constitution of the Wind Farm and the Location of the Shunt D-FACTS
The wind farm connected to the distribution system proposed in this section is shown in Figure 11, which consists of eight DFIGs with 1.5 MW of power for each wind turbine.These generators are connected between them to a voltage level of 30 kV by a step-up transformer 690 V/30 kV with 4 MVA of power for each generator.Then a 45 km line that is connected to the source substation 60 kV through another step-up transformer 30 kV/60 kV with 47 MVA of power.For this study, these lines are modelled with the π model.
Based on the work done in [21,35] the simulation results obtained with a D-FACTS provides an effective support to the bus voltage to which it is connected.Therefore, in this study, D-FACTS is placed at the PCC for two reasons:

•
The location for the reactive power support should be as close as possible to the point at which the carrier is necessary because of the variation in the voltage and, therefore, power loss (Joule loss) in the distribution line associated with reactive power flow, • In the studied system, the effect of the change in voltage is most common in this bus.

Simulation Results
In this section, simulations are performed on Matlab/Simulink, to show the impact of D-FACTS on the ability to control the voltage at the PCC between the electric distribution grid and the wind farm, which is described in the previous section (Figure 11).According to the previous section of the wind potential in the Algerian highland region, the considered wind speed in the simulation starts at 8 m/s and then reaches 9 m/s.The parameters of generators and D-FACTS are presented in the Appendix A.
For the more detailed study, the proposed system devices, the D-SVC and D-STATCOM structures used in the context of this paper are the same as the SVC and STATCOM structures, which are presented in [36,37], giving their associated models with their appropriate control schemes.In addition, the detailed model with some simulations for WECS based on DFIG in power system dynamics are described in [38][39][40][41].Generally, a short-circuit fault has a significant effect on the wind farm; a voltage drop is caused even if the fault is located near or far from the PCC or the wind farm.This voltage drop at the PCC leads to an over-current in the rotor circuit of DFIG, and fluctuations in the DC bus voltage.Therefore, the rotor side converter (RSC) of the DFIG should be blocked to avoid being damaged by overcurrent in the rotor circuit.
The block diagram of the simulated wind farm connected to the electrical distribution grid is shown in Figure 12.

Simulation Results
In this section, simulations are performed on Matlab/Simulink, to show the impact of D-FACTS on the ability to control the voltage at the PCC between the electric distribution grid and the wind farm, which is described in the previous section (Figure 11).According to the previous section of the wind potential in the Algerian highland region, the considered wind speed in the simulation starts at 8 m/s and then reaches 9 m/s.The parameters of generators and D-FACTS are presented in the Appendix A.
For the more detailed study, the proposed system devices, the D-SVC and D-STATCOM structures used in the context of this paper are the same as the SVC and STATCOM structures, which are presented in [36,37], giving their associated models with their appropriate control schemes.In addition, the detailed model with some simulations for WECS based on DFIG in power system dynamics are described in [38][39][40][41].Generally, a short-circuit fault has a significant effect on the wind farm; a voltage drop is caused even if the fault is located near or far from the PCC or the wind farm.This voltage drop at the PCC leads to an over-current in the rotor circuit of DFIG, and fluctuations in the DC bus voltage.Therefore, the rotor side converter (RSC) of the DFIG should be blocked to avoid being damaged by overcurrent in the rotor circuit.
The block diagram of the simulated wind farm connected to the electrical distribution grid is shown in Figure 12.In order to study the behavior and the impact of electrical faults in the distribution grid on the wind farm, worst-case scenarios of grid faults were assumed and included in the simulation.Therefore, the entire system was tested under two types of grid faults:

•
Line to line electrical grid fault.

•
Voltage drop at the 60 kV bus.
The wind speed is considered as constant during the grid fault period, except for this disturbance; and generators and the electrical distribution grid are considered to be working in ideal conditions (no disturbances and no parameter variations in the studied system).

Simulation Results of the Line-to-Line Electrical Grid Fault
In this section, we consider that the phases "b" and "c" at the PCC at P1 come into accidental contact.Then, the same grid fault is considered at a distance of 45 km from the PCC to the point P2 (Figure 11).Simulation results for this grid fault are shown in Figures 13-15.
The active powers at the PCC with a short temporary grid fault of two-phase to ground are shown in Figure 14.The reactive powers at the PCC with a short temporary grid fault of two-phase to ground are shown in Figure 15.
According to Figure 13, which reveals that without the use of D-FACTS systems the voltage at the PCC exceeds the acceptable voltage level 1 pu due to the voltage swell.However, using D-FACTS devices such as D-SVC and the D-STATCOM these undesirable effects are corrected.In addition, voltages at the PCC during a temporary grid fault presented in this figure show that, without the use of D-FACTS and when the grid fault is at the point P1, the voltage at the PCC drops to the value of 0.48 pu, which is less than the acceptable value.Thus, when the grid fault is at the point P2 without D-FACTS, the voltage drops to 0.52 pu.However, with the presence of D-FACTS devices, when this type of grid fault is at point P1, the voltage at the PCC drops to 0.63 with a slight fluctuation.In addition, when the fault is at the point P2 with D-FACTS, the voltage is maintained at 0.71 pu.Moreover, it is noticed that the voltage at the PCC, when using D-STATCOM many oscillations are mitigated compared to using the D-SVC.
According to Figure 14, which shows that during the occurrence of the same grid fault type in both points P1, P2 and without the presence of D-FACTS systems, no active power is supplied.Then, when the grid fault of 1000 ms duration exceeds the limit (see Figure 2), the wind farm is disconnected from the grid.Moreover, the installation of D-FACTS systems at the PCC guarantees the wind farm commissioning during and after this type of grid fault at these points (P1, P2) without disconnecting from the electrical distribution grid, providing the active power of 10.6 MW.Therefore, in the presence of D-FACTS systems, production of active power by the wind farm is uninterrupted and in the absence of these systems the wind farm is disconnected from the grid by triggering the protection system.
From Figure 15, it is noticed that in the absence of D-FACTS systems and when the grid fault is located at the points P1 and P2, no exchange of reactive power is provided to the electrical distribution grid.However, in the presence of D-FACTS systems, it provides almost the same amount of reactive power, 7.09 MVAr, when the grid fault is located at point P1 and 8.17 MVAr when the grid fault is located at point P2.Indeed, these injected reactive powers are required for compensation to maintain the stability of the wind farm with the voltage at the PCC around the acceptable value.Therefore, the wind farm is kept in service during and after this type of grid fault without disconnecting from the electrical distribution grid.Thus, the use of D-STATCOM has a capacity to compensate faster than using the D-SVC and the peaks of the injected reactive powers are eliminated.

Simulation Results of the Voltage Drop at the 60 kV Bus
The main purpose of this test is to study how a remote grid fault from the PCC may affect the operation of a wind farm; this fault affects the source of 60 kV, which is far from the PCC where the wind farm based on the DFIGs is connected to the grid.Hence, a temporary voltage drop of 50% is applied to the source for the duration of 500 ms at t = 10 s, and then the same grid fault for 1000 ms of duration.

Simulation Results of the Voltage Drop at the 60 kV Bus
The main purpose of this test is to study how a remote grid fault from the PCC may affect the operation of a wind farm; this fault affects the source of 60 kV, which is far from the PCC where the wind farm based on the DFIGs is connected to the grid.Hence, a temporary voltage drop of 50% is applied to the source for the duration of 500 ms at t = 10 s, and then the same grid fault for 1000 ms of duration.
Figure 16 shows that during this type of grid fault and without the presence of D-FACTS, the voltage at the PCC drops to 0.44.Therefore, the protection system will be triggered and the wind farm will be disconnected if the fault duration exceeds the electrical interconnection grid code for the wind turbine systems (see Figure 2).However, in the same figure, it is shown that during this grid fault and with the presence of D-FACTS systems, the voltage at the PCC is maintained around 0.88 pu with a transient peak without triggering the protection system.Thus, by using the D-STATCOM, transient peaks are reduced and the time response is faster than by using the D-SVC.
Appl.Sci.2018, 8, x 15 of 22 Figure 16 shows that during this type of grid fault and without the presence of D-FACTS, the voltage at the PCC drops to 0.44.Therefore, the protection system will be triggered and the wind farm will be disconnected if the fault duration exceeds the electrical interconnection grid code for the wind turbine systems (see Figure 2).However, in the same figure, it is shown that during this grid fault and with the presence of D-FACTS systems, the voltage at the PCC is maintained around 0.88 pu with a transient peak without triggering the protection system.Thus, by using the D-STATCOM, transient peaks are reduced and the time response is faster than by using the D-SVC.The active powers at the PCC are presented in Figure 17. Figure 17 shows that, when the D-FACTS devices are not installed at the PCC, the wind farm cannot maintain its connection to the grid during the grid fault that lasts 1000 ms because the protection systems are triggered and the wind farm is disconnected.However, after the installation of D-FACTS devices and with the same grid fault type and the same duration of grid fault, the wind farm can return to the steady state and inject the active power of 10.6 MW to the grid.
From the results shown in Figure 18, it is noticed that without the presence of a compensation system, the wind farm is operating in a weak electrical grid due to its normal behavior and there is no reactive power exchange with the electrical grid.However, in the same situation with the D-FACTS devices, the necessary reactive power is provided to the grid of 12.6 MVAr with D-SVC and 9.7 MVAr with the D-STATCOM.Thus, a very fast and significant fluctuation is observed with the use of the D-SVC compared to the use of the D-STATCOM.The active powers at the PCC are presented in Figure 17. Figure 17 shows that, when the D-FACTS devices are not installed at the PCC, the wind farm cannot maintain its connection to the grid during the grid fault that lasts 1000 ms because the protection systems are triggered and the wind farm is disconnected.However, after the installation of D-FACTS devices and with the same grid fault type and the same duration of grid fault, the wind farm can return to the steady state and inject the active power of 10.6 MW to the grid.
From the results shown in Figure 18, it is noticed that without the presence of a compensation system, the wind farm is operating in a weak electrical grid due to its normal behavior and there is no reactive power exchange with the electrical grid.However, in the same situation with the D-FACTS devices, the necessary reactive power is provided to the grid of 12.6 MVAr with D-SVC and 9.7 MVAr with the D-STATCOM.Thus, a very fast and significant fluctuation is observed with the use of the D-SVC compared to the use of the D-STATCOM.

Economic Analysis of D-FACTS Systems
The global market for FACTS and D-FACTS systems is expected to reach $1,386,010,000 in 2018, it had already reached $912,850,000 in 2012 [12].The D-SVC is the most widely used solution in the world market, followed by fixed capacitor banks.However, devices such as D-STATCOM are one customized solution for specific requirements of the distribution network.Obviously, some D-FACTSs are relatively expensive because they consist of many components such as advanced power electronics components, thyristors, reactors, capacitor banks, switches, protection systems and

Economic Analysis of D-FACTS Systems
The global market for FACTS and D-FACTS systems is expected to reach $1,386,010,000 in 2018, it had already reached $912,850,000 in 2012 [12].The D-SVC is the most widely used solution in the world market, followed by fixed capacitor banks.However, devices such as D-STATCOM are one customized solution for specific requirements of the distribution network.Obviously, some D-FACTSs are relatively expensive because they consist of many components such as advanced power electronics components, thyristors, reactors, capacitor banks, switches, protection systems and

Economic Analysis of D-FACTS Systems
The global market for FACTS and D-FACTS systems is expected to reach $1,386,010,000 in 2018, it had already reached $912,850,000 in 2012 [12].The D-SVC is the most widely used solution in the world market, followed by fixed capacitor banks.However, devices such as D-STATCOM are one customized solution for specific requirements of the distribution network.Obviously, some D-FACTSs are relatively expensive because they consist of many components such as advanced power electronics components, thyristors, reactors, capacitor banks, switches, protection systems and control systems.In this section, the range of the cost of the key features is often taken from the company Siemens and the Electric Power Research Institute (EPRI) with the database specified in [13], as shown in Figure 19.
control systems.In this section, the range of the cost of the key features is often taken from the company Siemens and the Electric Power Research Institute (EPRI) with the database specified in [13], as shown in Figure 19.Generally, the cost of a D-FACTS system has two components: the installation costs and operating expenses.The total cost of the entire installed systems comprises the equipment price and the delivery and installation of these systems.The operating cost includes the cost of maintenance and service.Specifically, the operating cost of these devices is approximately 5% to 10% of the total installation cost.Therefore, the cost functions for the D-SVC and D-STATCOM are developed as follows [15]:   From the above table, we see that the cost of D-FACTS devices (D-SVC and D-STATCOM) is much more expensive compared to capacitors due to the cost of the control devices and the complexity of the design and application of D-FACTS systems.The D-STATCOM is the source of reactive power compensation that is more expensive because of the used power electronics components like the Insolated Gate Bipolar Transistor (IGBT).Generally, the cost of a D-FACTS system has two components: the installation costs and operating expenses.The total cost of the entire installed systems comprises the equipment price and the delivery and installation of these systems.The operating cost includes the cost of maintenance and service.Specifically, the operating cost of these devices is approximately 5% to 10% of the total installation cost.Therefore, the cost functions for the D-SVC and D-STATCOM are developed as follows [15]: where s is the operating range of D-FACTS devices kVAr.The marginal cost per kVAr of the installed D-FACTS devices decreases as the operating rate of capacity increases.An overall cost for reactive power 100 MVAr, D-SVC ranges from $60 to $100 per kVAr.Although the D-SVC has sophisticated components such as thyristors, inductors and capacitors, it has a control structure that is relatively simple.Similarly, based on Figure 19, the overall cost of a D-STATCOM varies from $100 to $130 per kVAr and 100 MVAr of operating range.The costs of the installed parallel D-FACTS devices are shown in the Table 2 [16]: From the above table, we see that the cost of D-FACTS devices (D-SVC and D-STATCOM) is much more expensive compared to capacitors due to the cost of the control devices and the complexity of the design and application of D-FACTS systems.The D-STATCOM is the source of reactive power compensation that is more expensive because of the used power electronics components like the Insolated Gate Bipolar Transistor (IGBT).
In this study, the Figure 20 summarizes the performance of both D-FACTS types (D-STATCOM and D-SVC), also by comparing the amount of injected reactive power (MVAr) and the installation cost of these devices ($).This comparison provides a basis for system integration D-FACTS in a wind farm consisting of the DFIG type of generator helping to achieve a better balance between performance and cost in the condition of specific defects.The D-STATCOM provides a very effective reactive power compensation with respect to the D-SVC for all fault conditions.However, the D-SVC has a capacity option for a large amount of reactive power to be injected during a severe fault condition of the source.In the case of a two-phase ground fault (worst event of the grid fault applied in this study), the cost of installation is an important factor to consider.Therefore, in this paper, the economic analysis aims to compare the total cost of two types of parallel D-FACTS connected to a wind farm based on DFIGs, as presented in Figure 21.
According to Figure 21, it is clear that the use of D-STATCOM to maintain the wind farm in service during grid fault conditions is an expensive application compared with D-SVC due to the use of the transformer and the cost of power electronics [17,18].Consequently, one can conclude that the D-STATCOM is more cost-effective compared to D-SVC for voltage support at the PCC and the wind farm connection in the event of the most severe grid fault conditions.The D-STATCOM provides a very effective reactive power compensation with respect to the D-SVC for all fault conditions.However, the D-SVC has a capacity option for a large amount of reactive power to be injected during a severe fault condition of the source.In the case of a two-phase ground fault (worst event of the grid fault applied in this study), the cost of installation is an important factor to consider.Therefore, in this paper, the economic analysis aims to compare the total cost of two types of parallel D-FACTS connected to a wind farm based on DFIGs, as presented in Figure 21.In this study, the Figure 20 summarizes the performance of both D-FACTS types (D-STATCOM and D-SVC), also by comparing the amount of injected reactive power (MVAr) and the installation cost of these devices ($).This comparison provides a basis for system integration D-FACTS in a wind farm consisting of the DFIG type of generator helping to achieve a better balance between performance and cost in the condition of specific defects.
Reactive power (MVar) The D-STATCOM provides a very effective reactive power compensation with respect to the D-SVC for all fault conditions.However, the D-SVC has a capacity option for a large amount of reactive power to be injected during a severe fault condition of the source.In the case of a two-phase ground fault (worst event of the grid fault applied in this study), the cost of installation is an important factor to consider.Therefore, in this paper, the economic analysis aims to compare the total cost of two types of parallel D-FACTS connected to a wind farm based on DFIGs, as presented in Figure 21.
According to Figure 21, it is clear that the use of D-STATCOM to maintain the wind farm in service during grid fault conditions is an expensive application compared with D-SVC due to the use of the transformer and the cost of power electronics [17,18].Consequently, one can conclude that the D-STATCOM is more cost-effective compared to D-SVC for voltage support at the PCC and the wind farm connection in the event of the most severe grid fault conditions.According to Figure 21, it is clear that the use of D-STATCOM to maintain the wind farm in service during grid fault conditions is an expensive application compared with D-SVC due to the use of the transformer and the cost of power electronics [17,18].Consequently, one can conclude that the D-STATCOM is more cost-effective compared to D-SVC for voltage support at the PCC and the wind farm connection in the event of the most severe grid fault conditions.
Figure 22 shows the cost breakdown of a proposed 12 MW wind farm installation.
Figure 22 shows the cost breakdown of a proposed 12 MW wind farm installation.According to this figure, it can be seen that the proposed solution based on a D-FACTS system represents only 4% of the overall cost of the wind turbine installation.On the other hand, this solution offers good performance at the wind farm, ensuring its connection with the electrical grid and the reliability of the wind energy conversion system.

Conclusions
The aim of this paper is to investigate the feasibility study of the wind farm project in an Algerian highland region.Nevertheless, the study in this paper shows that Tiaret's electrical grid is susceptible to host the proposed wind farm in order to benefit from the wind potential.Therefore, in this paper, the application of the D-FACTS systems as the appropriate solution for the electrical grid connection issue and to accomplish the uninterrupted operation of a wind farm based on DFIG during the lineto-line fault and the voltage drop at 60 kV has been investigated.The D-FACTS is connected at the PCC where the wind farm is connected to the grid, to provide necessary reactive power for voltage support of the wind farm.Based on the simulation results, the stability improvement of the wind farm through the incorporation of the D-SVC or the D-STATCOM has been illustrated.In addition, it can be concluded that using D-STATCOM reduces the complexity of controlling the wind turbinegenerators, improves the time response of reactive power compensation, and corrects for the lack of a wind turbine.Future work will evaluate the impact of different scenarios for wind farms integrated into the Algerian electrical grid with other renewable energy sources and their electricity price.

Figure 2 .
Figure 2. Fault ride-through (FRT) profile according to the Spanish grid code.Reproduced from reference [23].

Figure 4 .
Figure 4. Annual wind speed in Algeria.

Figure 3 .
Figure 3. Wind speed variations at different locations in Algeria.

Figure 4 .
Figure 4. Annual wind speed in Algeria.

Figure 4 .
Figure 4. Annual wind speed in Algeria.

Figure 5 .
Figure 5. Map of the western Algerian electrical grid.Reproduced from reference [1].Figure 5. Map of the western Algerian electrical grid.Reproduced from reference [1].

Figure 5 .
Figure 5. Map of the western Algerian electrical grid.Reproduced from reference [1].Figure 5. Map of the western Algerian electrical grid.Reproduced from reference [1].

22 Figure 6 .
Figure 6.Western electrical grid near to the proposed wind farm.

Figure 7 .
Figure 7.In each bus of the 60-kV electrical grid in Tiaret region: (a) voltage amplitudes; (b) active powers.

Figure 8 .
Figure 8. Reactive power in each bus of the 60-kV electrical grid in the Tiaret region.

Figure 6 .
Figure 6.Western electrical grid near to the proposed wind farm.

Figure 6 .
Figure 6.Western electrical grid near to the proposed wind farm.

Figure 7 .
Figure 7.In each bus of the 60-kV electrical grid in Tiaret region: (a) voltage amplitudes; (b) active powers.

Figure 8 .
Figure 8. Reactive power in each bus of the 60-kV electrical grid in the Tiaret region.

Figure 7 .
Figure 7.In each bus of the 60-kV electrical grid in Tiaret region: (a) voltage amplitudes; (b) active powers.

Figure 6 .
Figure 6.Western electrical grid near to the proposed wind farm.

Figure 7 .
Figure 7.In each bus of the 60-kV electrical grid in Tiaret region: (a) voltage amplitudes; (b) active powers.

Figure 8 .
Figure 8. Reactive power in each bus of the 60-kV electrical grid in the Tiaret region.

Figure 8 .
Figure 8. Reactive power in each bus of the 60-kV electrical grid in the Tiaret region.

Figure 9 .
Figure 9.Control configuration of the distribution static var compensator (D-SVC) connected to the point of common coupling (PCC) with a wind farm.

Figure 9 . 22 Figure 10 .
Figure 9.Control configuration of the distribution static var compensator (D-SVC) connected to the point of common coupling (PCC) with a wind farm.This figure shows a D-SVC consisting of a thyristor switched capacitor (TSC) part composed by two switching thyristors connected with capacitive reactance X TSC ; the other part is the thyristor-controlled reactor (TCR) composed by two thyristors connected with an impedance of an inductive reactance branch X TCR .By controlling the angle thyristors (the angle with respect to the zero crossing of the phase voltage), the device is able to control the amplitude of the voltage at the point of common coupling (PCC) due to the changes in the angle resulting mainly in changes of the current.Therefore, the amount of the reactive power consumed by the inductor L, for an angle of α = 90 • , the inductive circuit is activated, whereas for α = 180 • , this inductive circuit is off.The configuration of the D-STATCOM connected to the distribution grid is shown in Figure10.This figure shows the different blocks constituting this configuration of the control strategy, which consists of: A phase lock loop (PLL) for synchronization of the component of the positive sequence voltage with the primary voltage of the power distribution grid.An external control loop consists of controlling the DC bus voltage and the grid voltage.The outputs of the voltage controllers are the current references for the current controllers.The internal current control loop consists of the current

Figure 10 .
Figure 10.Control configuration of the distribution static synchronous compensator (D-STATCOM) connected to the PCC with a wind farm.

Figure 11 .
Figure 11.Structure of the studied system based on a wind farm and a distributed-flexible AC transmission system (D-FACTS) connected to the Algerian distribution grid.

1 Figure 11 .
Figure 11.Structure of the studied system based on a wind farm and a distributed-flexible AC transmission system (D-FACTS) connected to the Algerian distribution grid.

Figure 12 .
Figure 12.Block diagram of the proposed system on Matlab/Simulink.Figure 12. Block diagram of the proposed system on Matlab/Simulink.

Figure 12 .
Figure 12.Block diagram of the proposed system on Matlab/Simulink.Figure 12. Block diagram of the proposed system on Matlab/Simulink.

Figure 14 .
Figure 14.(a) Active power at the PCC during line-to-line electrical grid fault; (b) zoom of active power at the PCC during line-to-line electrical grid fault.

Figure 14 .
Figure 14.(a) Active power at the PCC during line-to-line electrical grid fault; (b) zoom of active power at the PCC during line-to-line electrical grid fault.

Figure 15 .
Figure 15.(a) Reactive power at the PCC during line-to-line electrical grid fault; (b) zoom of reactive power at the PCC during line-to-line electrical grid fault.

Figure 15 .
Figure 15.(a) Reactive power at the PCC during line-to-line electrical grid fault; (b) zoom of reactive power at the PCC during line-to-line electrical grid fault.

Figure 16 .
Figure 16.Voltages at the PCC during the voltage drop at the 60 kV bus.

Figure 16 .
Figure 16.Voltages at the PCC during the voltage drop at the 60 kV bus.

Figure 17 .Figure 18 .
Figure 17.Active power at the PCC during the voltage drop at the 60 kV bus.

Figure 17 .
Figure 17.Active power at the PCC during the voltage drop at the 60 kV bus.

Figure 17 .Figure 18 .
Figure 17.Active power at the PCC during the voltage drop at the 60 kV bus.

Figure 18 .
Figure 18.Reactive power at the PCC during the voltage drop at the 60 kV bus.

Figure 19 .
Figure 19.Cost of the operation range for different D-FACTS devices.Reproduced from [14].
s is the operating range of D-FACTS devices kVAr.The marginal cost per kVAr of the installed D-FACTS devices decreases as the operating rate of capacity increases.An overall cost for reactive power 100 MVAr, D-SVC ranges from $60 to $100 per kVAr.Although the D-SVC has sophisticated components such as thyristors, inductors and capacitors, it has a control structure that is relatively simple.Similarly, based on Figure19, the overall cost of a D-STATCOM varies from $100 to $130 per kVAr and 100 MVAr of operating range.The costs of the installed parallel D-FACTS devices are shown in the

Figure 19 .
Figure 19.Cost of the operation range for different D-FACTS devices.Reproduced from [14].

Figure 20 .
Figure 20.Reactive power injected into the PCC by using D-SVC and D-STATCOM with different types of grid faults.

Figure 21 .
Figure 21.Cost of D-FACTS systems from different types of grid faults.

Figure 20 .
Figure 20.Reactive power injected into the PCC by using D-SVC and D-STATCOM with different types of grid faults.

Figure 20 .
Figure 20.Reactive power injected into the PCC by using D-SVC and D-STATCOM with different types of grid faults.

Figure 21 .
Figure 21.Cost of D-FACTS systems from different types of grid faults.

Figure 21 .
Figure 21.Cost of D-FACTS systems from different types of grid faults.

Figure 22 .
Figure 22.Breakdown of the costs of the wind project.

Table 1 .
Coordinates of stations at different Algerian sites.

Table 1 .
Coordinates of stations at different Algerian sites.

Table 1 .
Coordinates of stations at different Algerian sites.
MonthsFigure 3. Wind speed variations at different locations in Algeria.

Table 2 .
Costs of different reactive power compensators.

Table 2 .
Costs of different reactive power compensators.