Open access peer-reviewed chapter

Distributed Sources Optimal Sites and Sizes Search in Large Power Systems

By Mustafa Mosbah, Redha Djamel Mohammedi and Salem Arif

Submitted: March 1st 2020Reviewed: November 27th 2020Published: January 19th 2021

DOI: 10.5772/intechopen.95266

Downloaded: 35

Abstract

The integration of renewable sources into the power system has now become an unavoidable necessity for these technical and economic advantages and for the protection of the environment. In this chapter, a study is given for the integration of the Distributed Source (DS) in an optimal way and this by looking for the best location (sites) and the best power to be injected (size). The optimization technique used is based on genetic algorithms under technical and safety constraints, with the aim of minimizing active network losses and maximizing voltage stability. These objective functions are handled as a single and multi-objective problem. This study is applied on the standard IEEE 30 bus network under the MATLAB code.

Keywords

  • distributed source
  • optimal site and size
  • active power losses
  • voltage stability improvement
  • genetic algorithm

1. Introduction

Non-conventional or clean production is increasing rapidly around the world, thanks to it’s advantages such as reduced environmental impact, small size and it is a renewable energy [1]. The traditional power system is characterized by the unidirectional of the power flow because the energy comes from centralized sources [2]. These sources are generally based on fossil resources that are exhaustible and polluting for the environment.

The transmission power system spreads on long distances, which leads to losses in the lines by Joule effect on the one hand, and on the other hand, it requires huge investments and waste lands with long achievement times, which has led researchers to think about other solutions to face these problems [3].

Among these solutions is to have sources of electrical energy close to the consumers (local production). With the liberalization of the electricity market and the evolution of decentralized source (DS) technology in recent years, an increased trend towards their use has emerged. DS is defined as small producers based on renewable or conventional sources installed at different points in the power system, either at the transmission or distribution level [4]. The rate of DS integration tends to increase progressively in several countries [5].

In spite of the various advantages of DS, its sources present disadvantages when they are inserted into the power system, such as frequency instability caused by the intermittency of some renewable sources (e.g. photovoltaic and wind power) [6]. They also present some constraints, such as exceeding the thermal limits of power lines, increased Joule effect losses and the dysfunctioning of electrical protection devices and the exceeding of voltages at connection points. These constraints are due to the wrong choice of size (maximum power) and site of the DS, which requires the search for the best sites and adequate sizes [7].

A number of researchers have studied the problem of optimal site and size of DS to distribution networks [4, 8]. Example, in researchers [9] presented a technique a review of optimal DS placement in distribution network. DS site and size search techniques can be based on artificial intelligence techniques, metaheuristic techniques or deterministic techniques [10]. While most papers deal with this problem as a single-objective problem and some papers have focused on the problem of multi-objective optimizations [11]; these studies have considered a variety of objective functions including voltage profile improvement, losses minimization, reliability index minimization (SAIFI, SAIDI and END), fuel cost minimization, greenhouse gas emission minimization, maximization of DS penetration rate and maximization of voltage stability [12, 13, 14, 15, 16, 17, 18]. Some researchers have examined the problem of optimal sitting and size of DS to transmission networks [19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29].

Based on the literature search carried out, it was noted that the optimal integration of DS into distribution networks is largely discussed, however, there is less literature on the integration of DS into transmission networks.

The objective of this chapter is to study the problem of determining the Optimal Sites and Sizes of DS (OSSDS) in the transmission power system while taking into account the various constraints of the system (technical and security constraints). This is achieved by using a metaheuristic optimization technique such as the Genetic Algorithm (GA) technique. The objective functions considered in this work are the minimization of active power losses and the voltage stability improvement. These objective functions are treated as mono-objective and multi-objective (the objective functions are combined to a single objective function via weighting factors). This study has been applied on the IEEE 30 bus network under MATLAB code. For this reason, this chapter is organized according to the following plan: Section 2 will present the mathematical formulas of the OSSDS problem, i.e. the objective functions, the different constraints. Section 3 will present definitions on the method used (GA) and their application on OSSDS. Section 4 will give the description of the IEEE 30 node network and the limitations of the study framework. Section 5 will present the simulation, interpretation and analysis of the results obtained. Finally, the conclusion of the chapter and some perspectives.

2. Formulation of the OSSDS problem

The objective of this work is to research what is the optimal power to be injected by DS and the bus of their insertions that gives us the best performance of the power system considering the imposed constraints. To reach this objective we must define the fitness function and the constraints of equality and inequality which will be detailed in the following sections [30, 31].

2.1 Objective function (OF)

The OSSDS is formulated as a single and multi-objective problem using two objective functions.

2.1.1 First objective: minimization of active power losses

The first concern of the power system operator is to minimize active power losses in the line, this is expressed by the following equation:

TPLIk=k=1NBrkIk2E1

where rkand Ik, are the resistance and the current of the transmission line k respectively. TPLIrepresent the Total Power Losses Index of the network.

2.1.2 Second objective: improvement of the line voltage stability index

To assure that the operating point of the power system is far from the voltage collapse point, the voltage stability index must be improved. Among the effective indices that have been proposed in the literature is the Line Voltage Stability Index (LVSI) [32]. The LVSI of the line between bus iand jis given by the following equation:

LVSIk=4XijQjVisinθijδij21E2

where Xij, Qj, Vi, θij, and δijare the line reactance, is the reactive power at receiving end, the voltage at receiving end, is the line impedance angle between bus iand jand voltage phase angle difference between bus iand busj, respectively.

The line which represents an LVSI value close to 1 will be considered, the most critical line and may lead to the network collapse.

For multiobjective optimization, the two objectives functions are combined in a single linear function using weighting factors. Mathematically, this is given in the following form:

FMO=α1TPLIk+α2LVSIkE3

where α1and α2represent the weighting factors of active power losses index and the line voltage stability index, respectively.

2.2 Equality constraints

Equality constraints represent power balance equations between generation and demand. For a transmission power system in presence of DS, the active and reactive power equality constraints can be expressed as follows [33]:

i=1NGPGi+i=1NDSPDSii=1NDPDii=1NBPli=0i=1NGQGi+i=1NDSQDSii=1NDQDii=1NBQli=0E4

where PGQGis the active and reactive power of the conventional generator, PDSQDSis the active and reactive power produced by DS. This source is capable of delivering the active and reactive power. PDQDrepresent the active and reactive power demand at the load bus, PlQlrepresent the active and reactive power losses.

2.3 Inequality constraints

These constraints represent the physical limits of the lines, conventional generators and DS as also the security limits of the voltages of the network busses. They are expressed by the following equations:

SbiSbimaxfori=1..NBE5
PGiminPGiPGimaxfori=1.NGE6
QGiminQGiQGimaxfori=1..NGE7
ViminViVimaxfori=1NE8
SiminDSSiDSSimaxDSfori=1.NDSE9
i=1NDSPDSiτi=1NDPDifori=1NDSandi=1..NDE10

where Sbiis the apparent power that transits via the line between the nodes iand j. Sbimax, is the maximum limit of the line (thermal limit). PGimin, PGimax, PGiminand PGimaxare the minimum and maximum active and reactive powers of the ith conventional generator.SminDSand SmaxDSare the minimum and maximum production powers of the DS. N,ND,NG,NBand NDSare the number of network busses, number of load busses, the number of generators, the number of branches and the number of DS, respectively. Viminand Vimaxare minimum and maximum voltage in bus i.

For reasons of power system security, the network operator has limited the penetration rate (τ%) of DS. This limit is set depending on the robustness of the network (each network has its specific). This rate is calculated by the ratio of the total active power of the DS to the total active power demanded or the total power generated by conventional sources multiplied by one hundred; this is given by the following Equation [34]:

τ%=i=1NDSPDSii=1NDPDi100E11

2.4 Constraint processing

To optimize the site and size of the DS, it is important to mention that the control variables are generated within their allowable limits using a random strategy using a metaheuristic technique. During the optimization process, it is possible to come across solutions that are unfeasible due to exceeding the voltage limit, the thermal limit of the lines or the power limit of the reference generator, in this case the OF is penalized. This OF is reformulated as follows [33]:

Fp=i=1nOFfi+kPPslackPslacklim2+kVi=1NLVLiVLilim2+ksi=1NBSbiSbilim2E12

where kp, kVand ksare the penalty factors of the active power produced by the reference generator (slack bus), the voltage at the load nodes and the penalty constrained of the thermal limit of the lines, respectively.

The limits of the various variables are determined by the following equation:

Tilim=TimaxifTiminifTiifTi>TimaxTi<TiminTiminTiTimaxE13

T, represents the variables Pslack, VLiand Sbi.

The values of the penalty factors are determined by empirical means. After several tests it is decided that, kV=ks=10,000 and kp=1000.

From the above mathematical formulation, it can be noted that the use of such a model to solve real size problems is practically not possible with a classical approach. Consequently, in order to solve practical problems, it is necessary to use metaheuristic methods, that is, the genetic algorithms method.

3. Genetic algorithm method

GA is a metaheuristic optimization technique inspired by natural selection, and genetics developed by Holland-John, who conceived and realized an idea on how to transform the characteristics of natural evolution into a computer program [35]. The algorithm is based on a set of possible solutions randomly initialized in the search space. Individuals are represented by their design variables or by chromosome coding. Some solutions from the first population are used to form a new population based on genetic operators (crossover, mutation and selection). The goal is for the new population to be better than the previous one. The solutions that will be used to form new solutions are randomly selected according to their merit represented by an objective function specific to the problem posed, which should be minimized or maximized, so the better the individual, the greater his chances of surviving and reproducing, until the stop criterion is satisfied.

3.1 Application of the GA method

The goal of the GA method is to determine the optimal site and size of the DS to be integrated into the power system while minimizing the objective function under imposed constraints. Initially the vector of state variables and the vector of control variables are expressed as below:

The vector χTof state variables is composed of the active generation power at the reference bus PG1the load bus voltages VLi, the reactive power produced by conventional generators QGi, the apparent power transiting through transmission lines Sli. This vector χTcan be expressed by:

χT=PG1VL1VLNLQG1QGNGSl1SlNlE14

The vector vTof the control variables shows the active power outputs of the DS (size) SDSand the site or location of the DS LDS. Hence, this vector can be represented by:

vT=SDS1SDS2.SDSNDSLDS1LDS2..LDSNDSE15

Figure 1 indicates the chromosome structure employed in this study.

Figure 1.

Control variable vector structure.

Step 1: In determining the site and size of DS, the genetic algorithm method was suggested. The main steps for researching the site and size of DS are as follows:

Step 2: Run a power flow and determine the various network parameters in the absence and presence of DS, using the Newton Raphson method.

Step 3: Select the GA parameters (number of generations, population size, crossover and mutation probability) and randomly generate the values of the sites and sizes between their limits using Eqs. E16, E17 and E18 (creation of the initial population P0).

P0=X1X2Xi.XnpopulationE16
Lwi=round2+randNbus2E17
Swi=PminDSi+randPmaxDSiPminDSiE18

where roundrepresents the value of a number that is rounded to the next integer, and randrepresents a random number uniformly distributed between 0 and 1 [24].

Start the iteration meter t=0.

Step 4: Run the power flow in the presence of the DS and evaluate the objective function for each individual.

Step 5: Examine all network constraints using Eqs. E4 to E8. If the latter are satisfied, go to the next step. If not, penalize the OF using Eq. E12 and go to the next step.

Step 6: Generate new populations following the laws of GA (crossover, mutation and selection), increase the generation meter k=k+1and repeat steps 4 to step 6up to the stop criteria (maximum number of generations).

Step 7: Extract the best individual and show the results (site, size and various network parameters).

Figure 2 shows the flowchart of the search for the site and size of DS to be incorporated into the power system.

Figure 2.

DS site and size search flowchart.

4. Description of the studied network

The IEEE 30 bus network is used by power system researchers to test the effectiveness of their programs simulation. This system is composed of 30 buss, 41 transmission lines, 04 transformers, 06 conventional generators with a total generating capacity of 435 MW with a reactive capacity of −95 MVAr to 520 MVAr and 21 loads with a total power demand of 283.4 MW and 126.2 MVAr [35]. Bus No. 1 is taken as the reference (slack bus). Figure 3 shows the single line diagram of the network studied.

Figure 3.

Schematic diagram of the IEEE 30 bus network.

4.1 Study framework

Before beginning the presentation of the simulation results, it is necessary to cite the limitations of the study framework under consideration. The Newton Raphson method is used via the MATPOWER software to calculate the power flow. For the GA method the existing Toolbox in the MATLAB library (GA function) will be used. The voltage limits of the network busses are limited between 0.95 pu- 1.1 pu. DS are considered as sources capable of delivering active and reactive power, using a power factor of 0.8. DS integration are modeled as PQ busses (negative loads). It is important to note that all load nodes are considered candidate busses for DS sites. It is important to note that the number of DS to be integrated into the network is chosen in advance. After integration of a DS up to 10 DS, it is observed that, the best number from the viewpoint of improving the network parameters is 05 DS. For this reason, this study considers the insertion of five DS. The limits of the active power delivered by each DS are between 0 and 100 MW. And the penetration level limit of the total power delivered by the DS is 0.5i=1NDPDiMW that is i=1NDSPDSi0.5i=1NDPDi. After several tests on the studied network, it is noted that, the optimal parameters of the GA method are: the number of chromosomes (population size) is fixed at 100, the maximum number of generations is 50, the probability of mutation is 0.01 and the probability of crossover is 0.9. The multiobjective optimization considered in this study is based on the method of aggregate objectives through the weighting factors. It is important to note that, the choice of the values of these factors are dependent on the network operator and on the importance of the index to be improved. For this reason, the values of the weighting factors are set as follows: α1=0.8 and α2=0.2. In this study four cases are considered: the first case shows the simulation results for the base case (Execute the optimal power flow without DS). Second case, presents the results of the minimization of active losses only. Third case, presents the results of the voltage stability index improvement. Fourth case, presents the simulation results of the multi-objective case, it is the minimization of the bi-objectives at the simultaneous (minimization of losses and voltage stability index improvement).

5. Simulation results and interpretations

After the execution of the elaborated program, the simulation results achieved are shown bellow. Table 1, represents all the simulation results obtained for the five cases studied, as also the physical power limits of conventional generators. Figure 4, shows the situation of apparent power transmitted via transmission lines in all the cases studied. Figure 5 shows the voltage profile of several different cases.

ParametersLimitsCase 01Case 02Case 03Case 04
MaxMin
PG1(MW)20050182.6452.4760.6150.46
PG2(MW)802050.5850.5850.5850.58
PG3(MW)501518.5218.5218.5218.52
PG4(MW)351018.0918.0918.0918.09
PG5(MW)301010.4010.4010.4010.40
PG6(MW)401013.2613.2613.2613.26
QG1(MVAr)200−20−8.5616.946.1617.97
QG2(MVAr)100−2023.47−8.88−15.30−8.98
QG3(MVAr)50−1532.043.2223.75−5.10
QG4(MVAr)60−1549.29−5.34−2.80−4.82
QG5(MVAr)50−105.39−10.00−6.98−7.92
QG6(MVAr)60−152.74−14.26−15.00−12.29
Total Conventional Active Production435115293.49163.33171.46161.31
Total Conventional Reactive Production520−15104.37−18.32−10.16−21.14
Total Load (MW, MVAr)(283.40, 126.20)
DS Site(Best Location)******747
91217
122020
182323
302930
DS Size
(MW, MVAr)
**(46.96, 35.22)(36.40, 27.30)(67.52, 50.64)
(38.58, 28.93)(52.97, 39.73)(27.48, 20.61)
(5.63, 4.22)(6.07, 4.55)(10.36, 7.77)
(18.09, 13.56)(11.37, 8.52)(4.60, 3.45)
(13.60, 10.20)(9.73, 7.30)(15.05, 11.28)
TPLI
(MW)
10.892.7984.6202.945
LVSI0.1810.1580.0670.149
Loading Parameter (pu)2.14.45.75.2

Table 1.

Summary of network simulation results for the various cases.

** without DS integration.

Figure 4.

Apparent power of the transmission lines for several cases.

Figure 5.

Voltage values of the network bus for several cases.

From the results obtained by running the optimal power flow without DS (case 01), it is found that the total active losses are 10.89 MW, the LVSI is 0.181 and loading parameter (LP) is 2.1 pu.

In the case of the minimization the TPLI only (case 2), it can be seen that the total active losses are reduced to 2.798 MW, which represents a reduction rate of 74.3%, the LVSI has been improved with 12.7% and loading parameter represents a rate of 52.27%, that is after the integration of 05 DS with a total active power of 122.86 MW and reactive power of 92.13 MVAr distributed to the different sites (busses 7, 9, 12, 18 and 30) respecting all the constraints of the network.

When improving the LVSI only (case 3), the simulation results show that busses 4, 12, 20, 23 and 29 are selected as the best sites for DS power installations producing a total power of 116.54 MW and 87.4 MVAr. This implementation has enhanced the LVSI parameter from the value 0.181 to 0.067, which represents a rate of 62.9%, the losses are reduced to 4.62 MW, which is 57.5% and a loading parameter improvement rate of 63.15%.

The multiobjective optimization (case 04) proved that the results give a compromise between the different values of the minimized objective functions. In this case, the TPLI is improved by 72.9%, LVSI has been minimized by 17.67% and loading parameter has been increased by 59.61%. The advantage of Multi-objective optimization consists in improving several parameters of the electrical network and having the best compromise.

Figure 4 shows that all line powers are below the thermal limit power, which demonstrates the consideration of line stresses.

Figure 5 shows that all the bus voltages are included between the limits 0.95 pu and 1.1 pu. Figure 6 illustrates the PV curves before and after the integration of the DS for the three cases.

Figure 6.

(a) LP in case 2, (b) LP in case 3 (c) LP in case 4.

Figure 6 shows the positive influence of DS integration on the voltage stability margin of the IEEE 30 bus network. This figure also shows that case 3 represents the best value of LM, because the objective function only maximizes the voltage stability.

6. Conclusion

This chapter solves the problem of searching for the optimal location and size of DS in the transmission power system considering the required constraints, using the GA technique. In this work, various objective functions are targeted, minimizing active losses and voltage stability improvement. The last ones are treated as a mono and multi objective problem. The simulations carried out have shown that the results are dependent on the minimized objective function. The results achieved show the efficiency of the GA method. The perspective of this study is to include the constraint related to the transient stability of the rotor angles of conventional generators, in order to maintain the stability of the system during DS disconnection.

Acknowledgments

I thank very sincerely the members of the laboratory of LACoSERE, Electrical Department of the Amar Telidji University, Laghouat, Algeria and laboratory of LAADI, Electrical Department of Djelfa University, Djelfa, Alegria for their participations in the realization of this work.

The data (lines, consumptions, generators) of the IEEE 30 bus network used in this study are shown below [36]. Table 2 shows the bus data of the IEEE network studied. Table 3 represents the line data. Table 4 shows the data of the conventional generators.

Bus iTypePd(MW)Qd(MVAr)GsBsAreaVmVaBaseKVZoneVmaxVmin
1*3000011.06013211.10.95
2221.712.70011.043−5.4813211.10.95
312.41.20011.021−7.9613211.10.95
417.61.60011.012−9.6213211.10.95
5294.2190011.01−14.3713211.10.95
61000011.01−11.3413211.10.95
7122.810.90011.002−13.1213211.10.95
8230300011.01−12.113211.10.95
91000011.051−14.38111.10.95
1015.8201911.045−15.973311.10.95
112000011.082−14.391111.10.95
12111.27.50011.057−15.243311.10.95
132000011.071−15.241111.10.95
1416.21.60011.042−16.133311.10.95
1518.22.50011.038−16.223311.10.95
1613.51.80011.045−15.833311.10.95
17195.80011.04−16.143311.10.95
1813.20.90011.028−16.823311.10.95
1919.53.40011.026−173311.10.95
2012.20.70011.03−16.83311.10.95
21117.511.20111.033−16.423311.10.95
221000011.033−16.413311.10.95
2313.21.60011.027−16.613311.10.95
2418.76.704.311.021−16.783311.10.95
251000011.017−16.353311.10.95
2613.52.30011−16.773311.10.95
271000011.023−15.823311.10.95
281000011.007−11.9713211.10.95
2912.40.90011.003−17.063311.10.95
30110.61.90010.992−17.943311.10.95

Table 2.

Bus data of IEEE 30 bus power system.

Base MVA = 100

BusiBusjR(pu)x(pu)b(pu)RateA
(Slim MVA)
RateBRateCRatioAngleStatus
120.01920.05750.052813000001
130.04520.16520.040813000001
240.0570.17370.03686500001
340.01320.03790.008413000001
250.04720.19830.041813000001
260.05810.17630.03746500001
460.01190.04140.0099000001
570.0460.1160.02047000001
670.02670.0820.01713000001
680.0120.0420.0093200001
6900.208065000.97801
61000.556032000.96901
91100.20806500001
91000.1106500001
41200.256065000.93201
121300.1406500001
12140.12310.255903200001
12150.06620.130403200001
12160.09450.198703200001
14150.2210.199701600001
16170.05240.192301600001
15180.10730.218501600001
18190.06390.129201600001
19200.0340.06803200001
10200.09360.20903200001
10170.03240.084503200001
10210.03480.074903200001
10220.07270.149903200001
21220.01160.023603200001
15230.10.20201600001
22240.1150.17901600001
23240.1320.2701600001
24250.18850.329201600001
25260.25440.3801600001
25270.10930.208701600001
282700.396065000.96801
27290.21980.415301600001
27300.32020.602701600001
29300.23990.453301600001
8280.06360.20.04283200001
6280.01690.05990.0133200001

Table 3.

Line data of IEEE 30 bus power system.

BusPg
(MW)
Qg
(MVAr)
Qmax
(MVAr)
Qmin (MVAr)VgS BaseStatusPmax (MW)Pmin (MW)
100200−201.06100120050
250.58460100−201.04510018020
518.5227050−151.0210015015
818.0865060−151.02910013510
1110.4038050−101.0610013010
1313.2553060−151.0610014012

Table 4.

Generator data of IEEE 30 bus power system.

© 2021 The Author(s). Licensee IntechOpen. This chapter is distributed under the terms of the Creative Commons Attribution 3.0 License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

How to cite and reference

Link to this chapter Copy to clipboard

Cite this chapter Copy to clipboard

Mustafa Mosbah, Redha Djamel Mohammedi and Salem Arif (January 19th 2021). Distributed Sources Optimal Sites and Sizes Search in Large Power Systems, Renewable Energy - Technologies and Applications, Tolga Taner, Archana Tiwari and Taha Selim Ustun, IntechOpen, DOI: 10.5772/intechopen.95266. Available from:

chapter statistics

35total chapter downloads

More statistics for editors and authors

Login to your personal dashboard for more detailed statistics on your publications.

Access personal reporting

Related Content

This Book

Next chapter

New Generation Aero Combustor

By Jushan Chin and Jin Dang

Related Book

First chapter

Introductory Chapter: An Overview of PEM Fuel Cell Technology

By Tolga Taner

We are IntechOpen, the world's leading publisher of Open Access books. Built by scientists, for scientists. Our readership spans scientists, professors, researchers, librarians, and students, as well as business professionals. We share our knowledge and peer-reveiwed research papers with libraries, scientific and engineering societies, and also work with corporate R&D departments and government entities.

More About Us