Open access peer-reviewed chapter

Particle Swarm Optimization Solution for Power System Operation Problems

By Mostafa Kheshti and Lei Ding

Submitted: June 2nd 2017Reviewed: November 13th 2017Published: December 20th 2017

DOI: 10.5772/intechopen.72409

Downloaded: 290

Abstract

Application of particle swarm optimization (PSO) algorithm on power system operation is studied in this chapter. Relay protection coordination in distribution networks and economic dispatch of generators in the grid are defined as two of power system-related optimization problems where they are solved using PSO. Two case study systems are conducted. The first case study system investigates applicability of PSO on providing proper overcurrent relay settings in the grid, while in the second case study system, the economic dispatch of a 15-unit system is solved where PSO successfully provides the optimum power output of generators with minimum fuel costs to satisfy the load demands and operation constraints. The simulation results in comparison with other methods show the effectiveness of PSO against other algorithms with higher quality of solution and less fuel costs on the same test system.

Keywords

  • power system
  • economic dispatch
  • relay coordination
  • particle swarm optimization

1. Introduction

Electric power system is the most complex man-made system, and the modern society depends heavily on continuous and reliable operation of this system to supply electricity to commercial, residential and industrial consumers. Operation of the power grid involves a balanced platform in generation, transmission and distribution, which costs billions of dollars to run. The reliable and continuous availability of electricity with minimum costs is the major objective of utility grids and energy providers. Two of the important complex problems in power system are economic dispatch and power system protection.

The power plants and utility grids need to allocate the available generation units in an efficient and economical way to respond to the load demand in order to provide continuous power supply in stable conditions and with minimum power production costs. This is addressed as economic dispatch (ED) [1, 2]. With the practical constraints on the generators, finding optimum power outputs with minimum fuel costs is challenging.

In addition, as the occurrence of failures and faults in the power grid is inevitable, the entire power system must be protected. The relay protection scheme is designed to detect faults and isolate the faulty parts of the grid from the healthy sections in order to mitigate the consequences of the faults and maintain continuity of service. If a fault occurs, the nearest corresponding relays must operate as fast as possible to clear the fault. If due to any reason these primary relays fail to react, their backup relays must operate and accomplish the task. Directional overcurrent relays (DOCRs) are a suitable and economical protection scheme for distribution systems [3]. The protection design of DOCRs is based on two parameters, time multiplier setting (TMS) and plug setting (PS). Proper settings of TMS and PS allow a primary relay to clear the faults in its protection zone as fast as possible and in case of failure, its backup relay operates immediately after a time interval to clear the fault. TMS and PS values of each relay must be coordinated with other backup relays, where again relays act with different current settings, which make the coordination a complex task. Each pairs of relays include four variables (TMS, PS) and the complexity of coordination will be intense in bigger systems with more relays and constraints. Due to the complex interconnection of the distribution systems and also nonlinear characteristics of operation time of relays, finding best relay settings could be very difficult.

Considering the non-convex and nonlinear nature of these problems, traditional methods fail to feasibly or optimally solve them. Therefore, evolutionary algorithms have gained more attentions as solutions to such optimization problems. Some of the recent related works on ED problem have been studied with the metaheuristic methods such as Genetic Algorithm (GA) [4], Particle Swarm Optimization (PSO) [4], Imperialist Competitive Algorithm (ICA) [5], Artificial Bee Colony (ABC) [6], Bacterial Foraging Optimization (BF) [7], Hybrid Harmony Search with Arithmetic Crossover (ACHS) [8], GA with a special class of ant colony optimization (GAAPI) [9] and so on. The modified and hybrid models of PSO such as Modified PSO (MPSO) [10], guaranteed convergence PSO (GCPSO) [10], Species-based Quantum Particle Swarm Optimization (SQPSO) [11], Iteration PSO (IPSO) [12], Parallel PSO with Modified Stochastic Acceleration Factors (PSO-MSAF) [13], Distributed Sobol PSO and Tabu Search Algorithm (DSPSO-TSA) [14], Self-Organizing Hierarchical PSO (SOH-PSO) [15], Passive Congregation-based PSO (PC-PSO) [15] and Simple PSO (SPSO) [15] have also been employed to address the ED problem.

Application of metaheuristic algorithms on power system protection and particularly, DOCR coordination in distribution networks has been introduced in literature such as PSO [3, 16], Harmony Search Algorithm (HSA) [17], Cuckoo Algorithm [18], chaotic firefly algorithm [19], differential evolution [20] and so on.

In this chapter, PSO is applied as a solution to the introduced power system operation problems, namely ED and DOCR coordination. The rest of the chapter is organized as follows: in Section 2, these power system problems are defined and formulated as optimization problems. PSO algorithm is explained in Section 3. In Section 4, PSO is applied in two case study systems to conduct the performance and feasibility of this method. Finally, Section 5 concludes the chapter with the results in pervious sections.

2. Problem formulation

In this section, overcurrent relay coordination and economic dispatch problems are formulated separately as optimization problems.

2.1. Relay coordination problem

In a protection scheme, each primary relay should operate as fast as possible to clear the fault in a system. If the operation time of the relay takes longer than an acceptable time, the damage on the faulty equipment would be severe with serious consequences. In other words, minimizing the total operation time of relays decreases the risk and stress on the protected apparatus, which can be depicted as an optimization objective function:

OF=mini=1nwitiE1

where tiis the operation time of relay Ri, wiis the probability of the occurrence of fault on transmission line in the zone of protection, and it is normally set to 1; nis the total number of relays in the system.

Generally, the operation time of DOCRs is defined in (2):

ti=λ×TMSiIFiPSiη1+LE2

where IFiis the fault current seen by the appropriate relay Riafter being transformed through the secondary winding of corresponding current transformer (CT). Depending on the type of relays, the characteristic constants λ, L and η are selected [16]. In this chapter, continuous form of TMS and PS is considered with relay type of standard inverse definite minimum time (IDMT). Based on that, all the relays in the system are assumed identical with a common characteristic function approximated by:

ti=0.14×TMSiIiPSi0.021E3

To ensure that the operation time of an individual relay is proper enough to mitigate the damage impact of faults on the apparatus, the time must be within an acceptable range:

timintitimax;i=1,,nE4

where, respectively, timinand timaxare the minimum and maximum operating time of the relay Ri. Each overcurrent relay has a manufactured TMS range to provide controllability of response to faults with different speeds. As shown in (3), tiis proportional to TMS values. Also, the PS has nonlinear effect on the operating time. Within a security margin and to avoid maloperation of an individual relay with normal load or slight overload current, the minimum pickup current setting is selected bigger than the maximum load current. The maximum plug setting is chosen not greater than the minimum fault current [19]. Therefore, there are constraints on TMS and PS as follows:

TMSiminTMSiTMSimax;i=1,,nE5
PSiminPSiPSimax;i=1,,nE6

where TMSimax, PSimax, TMSiminand PSiminare the maximum and minimum values of TMS and PS of relay Ri. Although the constraints in Eqs. (4)(6) seem to provide satisfactory limits on performance of each relay, they are not enough to guarantee correct performance of the protection scheme as the coordination between primary-backup relays has not been considered.

The constraints can only ensure the operation of an individual relay not a primary-backup pair. To coordinate adjacent relays as primary and backup relays, the primary relay should operate as fast as possible within its acceptable boundaries. If it fails to act, its backup relay needs to take over the tripping action with a minimum time. The minimum operating time of the backup relay must be small but yet bigger than the operating time of primary relay. Therefore, a coordination time interval (CTI) is added to the constraints to satisfy the proper coordination scheme.

tjbackuptiprimaryCTIE7

where tiand tjare the operation time of primary and backup relays, respectively. CTI depends upon the relay type, circuit breaker speed, relay over-travel time and the safety factor time for CT saturation, setting errors, contact gaps and so on. According to the IEEE standard [21], CTI is set to 0.2 s for the digital relays. Figure 1 shows the backup-primary pair of relays in a radial network.

Figure 1.

Primary-backup relation between relays in a distribution system.

2.2. Economic dispatch

How to allocate the available generators in the grid to respond to the load demand with minimum fuel costs is an economical aspect of power dispatch in the electric power system that annually costs millions of dollars to operate.

In a practical ED optimization, the generator constraints and network limits such as ramp rate limit, the prohibited zones of operation, generation capacity constraints and valve point effects are considered. Single quadratic equation is used to formulate the ED optimization problem:

minFt=i=1mFPi=i=1mαi+βiPi+γiPi2E8

where αi, βiand γiare the cost equation coefficients of unit i, Piis the output power in MW and FPiis the cost function of that unit in $/h. The index mdenotes the number of generators in a system.

The fuel cost function in (8) as the objective function of ED problem here is associated with practical constraints. Considering the ramp rate limits, the momentary output power of a generator cannot exceed its previously generated power more than a certain amount of URi, the up-ramp rate limit and neither can it be less than that of the previously generated power by more than a certain amount of DRi, the down-ramp rate limit of the generator. Therefore:

PiPi0URiE9
Pi0PiDRiE10

where Piis the current power output and Pi0refers to the previous power output of generator i. URiand DRirepresent the up-ramp limit and down-ramp limit of the generator i, respectively, in which the current generated power has the following constraint in MW/t:

Pi0DRiPiPi0URiE11

The input-output curves of the generation units have separate operation zones. The prohibited zones of operations are due to the operation of steam valve or the shaft bearing vibration of the generators. Therefore, the generated power is within the feasible zones of operation and outside the prohibited zones. For a generator i:

PiPiminPiPi,1lPi,j1uPiPi,j1lPi,nuPiPimax,j=2,3,,n;i=l,,mE12

where jis the number of prohibited zones of operation for unit i, Pi,jland Pi,juare the lower and upper boundaries, respectively, of prohibited zone jof generator i, Piminand Pimaxare minimum and maximum power capacity, respectively, of a generator that produces. From (11) and (12), it can be deduced that the operational power of a generator must be within the constraint in (13) and also be out of the prohibited zones in (12):

maxPiminPi0DRiPiminPimaxPi0+URiE13

Total delivered power from the units needs to meet the power demand and transmission loss in the grid.

i=1mPi=PD+PL;imE14

where PDand PLrepresent the demand power and the power loss in the gird, respectively. The overall power loss of the committed units is based on the output power, which is formulated by Bmatrix coefficients known as Kron’s formula:

PL=i=1mj=1mPiBijPj+i=1mB0iPi+B00E15

The power loss itself cannot be more than some permissible values:

PLf,kPLf,kmax;k=j,,LE16

where the real power flow of line jis represented with PLf,kand k is the number of transmission lines in a system. The power loss cannot be more than a maximum value of PLf,kmax.

3. Particle swarm optimization algorithm

PSO algorithm is a nature inspired method from social behavior of bird flocking and fish schooling, which is first introduced by Erberhart and Kennedy in 1995 [22]. As a population-based stochastic optimization technique using swarm intellects in the search space, this technique is based on interaction of swarm of particles. Every particle includes two values of position and velocity that are updated during the iteration runs by considering each particle’s best experience (best position) and the best achieved experience (global position) of all particles.

The update of position and velocity of the particles must be processed, and it has to follow Eqs. (17)(19) for satisfying the constraint of an optimization problem. Each particle movement is based on the changes in its position and its velocity:

pik+1=pik+vik+1E17

where pik+1and pikare the position of particle iin the iteration k+1and k, respectively, vik+1is the velocity of the particle in k+1iteration. A particle’s velocity is defined as follows:

vik+1=w×vik+c1×rand1×pbestipik+c2×rand2×gbestpikE18

where pbestiis the best so far position of the particle ias the best experience, while gbestis the best position among the whole swarm with all the particles in movement as global experience. c1and c2are the weighting factors, while rand1and rand2are two random numbers between zero and one. The parameter wis the inertia factor varying between wminwmax, as shown in (19), which is a linear decreasing inertia weight in this chapter.

w=wmaxwmaxwminitermax×kE19

where kand itermaxare the current iteration and the maximum number of iteration during simulations, respectively.

The general steps of PSO on solving an optimization problem are as follows:

  1. Set initial parameters wmin, wmax, c1and c2.

  2. Generate initial populations having initial positions pand velocities v.

  3. Set iteration k=1.

  4. Calculate fitness of particles Fik=fpik,iand find the index of the best particle b.

  5. Select pbestik=pik,iand set gbestk=pbk.

  6. Update inertia factor: w=wmaxwmaxwminitermax×k.

  7. Update the velocity and position of particles.

    vik+1=w×vik+c1×rand1×pbestipik+c2×rand2×gbestpik;ipik+1=pik+vik+1;i

  • Calculate fitness Fik+1=fpik+1,iand obtain the index of the best particle b1.

  • Update pbestfor all particles.

    ifFik+1<Fikthenpbestik+1=pik+1elsepbestik+1=pbestik

  • Update gbestof the population

    ifFb1k+1<Fbkthengbestk+1=pbestb1k+1andb=b1elsegbestk+1=gbestk

  • If k<Itermaxthen k=k+1and go to step 6 else go to step 12.

  • Print gbestkas optimum solution.

  • 4. Simulation results

    Application of PSO on solving the defined problems in previous sections is validated here. In the first case study system, the overcurrent relay coordination problem in a distribution network is solved. The size of population and iteration numbers of PSO are 60 and 100, respectively. In the second case study system, the ED problem is addressed with population size of 100 and 500 iterations. c1, c2, wminand wmaxin both case study systems are the same and set to 2, 2, 0.4 and 0.9, respectively.

    4.1. Case study 1: relay protection coordination

    A 15-node radial network system including a total of 13 loads is considered, and PSO is employed to determine the optimal settings of all 28 digital over-current relays shown in Figure 2. Identical digital relays are used with same current transformer (CT) ratio 500:1. The constraint values of the relays are as follows: timinand timaxfor each relay is set to 0.1 and 4 s, respectively; TMSminand TMSmaxconstraints are 0.1 and 1.1, respectively; PSminand PSmaxare 0.5 and 2.5, respectively [16].

    Figure 2.

    A 15-node distribution system as case study one.

    The primary-backup relationships of the relays in the system are shown in Table 1. According to the objective function in (3), the maximum fault currents sensed by the relays are also required. Therefore, the case study system in Figure 2 has been modeled in DigSILENT PowerFactory software with simulating three-phase faults occurring in front of each relay. The collected data have been shown in Table 2. The load parameters of the radial network are shown in Table 3.

    PrimaryBackupPrimaryBackupPrimaryBackupPrimaryBackup
    18151322
    291162321
    3110171324
    411918255
    531219326
    613120275
    7514211928

    Table 1.

    Primary-backup pair relationship of the relays.

    RelayFault current (kA)RelayFault current (kA)RelayFault current (kA)RelayFault current (kA)
    122.77884.636157.7224.636
    211.532911.509165.79234.632
    311.509107.71177.7243.866
    47.71117.7185.79255.785
    57.7125.79197.7264.636
    65.791311.509205.79275.785
    75.785147.71215.785284.636

    Table 2.

    Maximum fault currents.

    LoadActive power (MW)Reactive power (Mvar)LoadActive power (MW)Reactive power (Mvar)
    141.580.50.2
    210.2920.8
    32.511011
    432112.50.9
    510.51210.5
    630.71332
    710.5

    Table 3.

    Load parameters in case study one.

    After data collection in DigSILENT PowerFactory software, the optimal settings (TMS and PS) of all 28 overcurrent relays are obtained by solving (3) subject to the constraints in (4)(7) using the PSO in MATLAB software. The algorithm was executed 100 times to achieve accurate results. Table 4 shows the best TMS and PS settings of the relays to provide a reliable protection scheme in the distribution system. The total operation time of the relays in this system is 26.189 s. Figure 3 illustrates the convergence curve of PSO in solving the objective function.

    RelayTMSPSRelayTMSPS
    10.743940.54999150.426320.79222
    20.322161.2551160.543960.67839
    30.34161.5153170.217940.61987
    40.552452.4141180.120371.1896
    50.363420.60859190.319970.66671
    60.43091.0805200.565661.1323
    70.111881.2611210.263760.60754
    80.27261.822220.143131.8741
    90.453951.1148230.124960.7978
    100.153591.3933240.553021.3848
    110.321330.50089250.222540.65529
    120.264182.4727260.29041.6316
    130.320881.9075270.134710.77885
    140.738992.1497280.102772.4925
    OF (s)26.189

    Table 4.

    Obtained TMS and PS values of relays by PSO.

    Figure 3.

    Convergence of PSO solution in case study one.

    4.2. Case study 2: economic dispatch

    A 15-unit test system is used to investigate the feasibility of PSO in solving the nonsmooth economic dispatch considering transmission losses, ramp rate limits and the prohibited operation zones of the thermal generators. The cost curves data and operation limits of the 15-unit system are shown in Table 5. The B loss coefficients to calculate the power loss can be found in [4].

    UnitPiminPimaxPi0αiβiγiURiDRiProhibited zones (MW)
    115045540067110.10.00029980120
    215045530057410.20.00018380120[185 225] [305 335] [420 450]
    3201301053748.800.001126130130
    4201301003748.800.001126130130
    51504709046110.400.00020580120[180 200] [305 335] [390 420]
    613546040063010.100.00030180120[230 255] [365 395] [430 455]
    71354653505489.80.00036480120
    8603009522711.20.00033865100
    92516210517311.20.00080760100
    102516011017510.70.00120360100
    1120806018610.20.0035868080
    122080402309.900.0055138080[30 40] [55 65]
    1325853022513.10.0003718080
    1415552030912.10.0019295555
    1515552032312.40.0044475555

    Table 5.

    Generation unit characteristics of a 15-unit system in case study two.

    This system has many local minima with high dimensionality that draws realistic analysis for practical applications. PSO is applied on the 15-unit test system, and the results are compared with best results in literature on the same system: ACHS [8], SQPSO [11], ABC [6], IPSO [12], PSO-MSAF [13], DSPSO-TSA [14], ICA [5], GAAPI [9], MPSO [10], SOH-PSO [15], GCPSO [10], PC-PSO [15], BF [7], SPSO [15], PSO [4] and GA [4].

    The simulation is tested for 100 times to ensure reliable analysis. Table 6 shows the optimum results of each method for the 15-unit system. Figure 4 illustrates the convergence of PSO method in 100 different trials, while the final fuel costs in 100 trials are shown in Figure 5. Table 7 shows the best economic dispatch of power using PSO for ED optimization in the 15-unit system.

    MethodBest cost ($/h)Mean ($/h)
    PSO32701.28232704.10578
    ACHS [8]32706.650032706.65
    SQPSO [11]32706.674032708.4457
    ABC [6]32707.8532707.95
    IPSO [12]32709.0032784.5
    PSO-MSAF [13]32713.0932759.64
    DSPSO-TSA [14]32715.0632724.63
    ICA [5]32715.4305NA
    GAAPI [9]32732.95NA
    MPSO [10]32738.4177NA
    SOH-PSO [15]32751.3932,878
    GCPSO [10]32764.4616NA
    PC-PSO [15]32775.36NA
    BF [7]32784.502432796.81
    SPSO [15]32798.69NA
    PSO [4]32,85833,039
    GA [4]33,11333,228

    Table 6.

    Comparison of PSO results with other methods in case study two.

    Figure 4.

    Convergence curve of PSO over 100 different trials.

    Figure 5.

    Cost distribution with PSO over 100 different trials.

    Unit (MW)PSOICA [5]GCPSO [10]MPSO [10]GA [4]PSO [4]ABC [6]SOH-PSO [15]
    P1454.9963455449.89252455415.3108439.1162455455
    P2379.9998380366.99066380359.7206407.9727380380
    P3130130130130104.4250119.6324130130
    P4129.995413013013074.9853129.9925130130
    P5169.9999167.4174170170380.2844151.0681169.9997170
    P6459.9999460460460426.7902459.9978460459.96
    P7430430430430341.3164425.5601430430.00
    P866.1794113.473775.8846092.7278124.786798.569971.9698117.53
    P964.948525.155550.2268943.0282133.1445113.493659.179877.90
    P10159.2255155.3478160140.193889.2567101.1142159.8004119.54
    P1179.999680808060.057233.91168054.50
    P1279.99018077.870638049.999879.95838080.00
    P1325.0001252527.640338.771325.004225.002425.00
    P1415.00051515.831220.761041.942541.414015.005617.86
    P1515.00291539.6614622.272422.644535.614015.001415
    Total power2660.3382661.3942661.358062661.62352668.42262.42735.9592662.29
    Power loss30.33831.29130.8659329.97838.278232.430630.959132.28
    Cost ($)32701.28232715.430532764.461632738.41778331133285832707.8532751.39

    Table 7.

    Optimum solution of PSO compared with other methods in case study two.

    The optimum cost of 15-unit system with the proposed PSO solution is 32701.282 ($), which has better result than other methods. GA has the most deviating results compared with other hybrid and improved methods in this test system. Also, the mean value of final fuel cost of generators using PSO over 100 trials is less than minimum values obtained by other method, which indicates higher quality of solution and better performance of PSO compared with other algorithms on the same test system. The power loss in the grid obtained from PSO is less than other algorithms, which shows better dispatching scheme using PSO. The best convergence of PSO to the minimum fuel cost of generators in the grid while satisfying the constraints is shown in Figure 6.

    Figure 6.

    Convergence of PSO with best fuel cost result for ED problem in case study two.

    In case study two, PSO achieves better results when compared with other hybrid or improved methods. It is worth mentioning that the maximum iteration number is 500 and the population size is 100. In most of the quoted methods, the iteration and population sizes vary that can affect the final results. For example, the PSO in [4] has population size and iteration number of 100 and 500, respectively. The population size, iteration, crossover rate, mute rate and crossover parameter of GA [4] are 100, 200, 0.8, 0.01 and 0.5, respectively.

    Different system configuration and programming language frameworks can also influence the results in which MATLAB 2015Ra was used for programming.

    5. Conclusions

    Distribution network relay coordination and the economic dispatch of generators in the electric power system were modeled as optimization problems. Particle swarm optimization (PSO) was successfully employed to solve the defined problems where two case study systems were conducted to validate the results. In the first case study system, PSO provided proper relay settings that allow all the relays in a system to perform with high reliability and accuracy. In the second case, the optimal power outputs of thermal generators in the grid were scheduled to satisfy the load demands and other practical constraints on the generators and the grid with minimum fuel costs. The compared results with other methods demonstrated higher quality of solution, and less fuel costs obtained by PSO. The general performance of PSO in this chapter indicates applicability of this method on practical power system-related problems that are difficult to be handled by conventional methods.

    How to cite and reference

    Link to this chapter Copy to clipboard

    Cite this chapter Copy to clipboard

    Mostafa Kheshti and Lei Ding (December 20th 2017). Particle Swarm Optimization Solution for Power System Operation Problems, Particle Swarm Optimization with Applications, Pakize Erdoğmuş, IntechOpen, DOI: 10.5772/intechopen.72409. Available from:

    chapter statistics

    290total 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

    Stochastic Greedy-Based Particle Swarm Optimization for Workflow Application in Grid

    By Ruey-Maw Chen and Yin-Mou Shen

    Related Book

    First chapter

    Genetic Algorithm Optimization of an Energy Storage System Design and Fuzzy Logic Supervision for Battery Electric Vehicles

    By Stefan Breban

    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