# A Parallel Implementation of Unscheduled Flow Control in Interconnected Power Systems

Authors: G. Ozdemir Dag; Mustafa Bagriyanik
Journal: Mathematical Problems in Engineering
 (2012)

---

## Abstract

The unscheduled power flow problem needs to be minimized or controlled as soon as possible in a deregulated power system since the transmission systems are mostly operated at their power-carrying limits or very close to it. The time spent for simulations to determine the current states of all the system and control variables of the interconnected power system is important. Taking necessary action in case of any failure of equipment or any other occurrence of an undesired situation could be critical. Using supercomputing facilities and parallel computing techniques together decreases the computation time greatly. In this study, a parallel implementation of a multiobjective optimization approach based on both genetic algorithms and fuzzy decision making to manage unscheduled flows is presented. Parallel computation techniques are applied using supercomputers (high-performance computers). The proposed method is applied to the IEEE 300 bus test system. Two different cases for some parameters of GA are considered to see the power of parallel computation technique. Then the simulation results are presented.


---

## Body

## 1. Introduction

After the liberalization of the electricity markets, the operation of the power system has fundamentally changed. How much power flows on which transmission lines, which company uses the other’s transmission lines, and/or the amount and duration of the transmission line usage have all become very important issues.The path electrical energy takes depends on physical laws. That is, Kirchhoff’s current and resistance laws determine the path for the energy to flow. Energy takes the shortest path (in terms of resistivity) instead of a contracted path. Several power flow cases could exist in a transmission network depending on the system topology and operating conditions. When the actual power flow exceeds system-operating limit for a transmission path, the transmission system operators must immediately mitigate the transmission overloads to reduce the actual power flow across the path. The deviation of actual electric power flows in transmission circuits from the scheduled (expected) power flows is called unscheduled flows, and these flows may result in the congestion problems and increase the active and reactive power transmission losses in the lines depending on the power to be transmitted. This can be best explained using Figure1. System B buys power from system A; however, power travels through system C occupying different transmission lines. One needs to control (i.e., redirect the power to the transmission line directly connecting A and B) or prevent power flows between A and C. The high power losses, low efficiencies, or the path power travels through (occupying transmission lines) before arriving to the loads has not seen as an important problem in a government-controlled power system [1–3].Figure 1
The loop flow.Under the competitive market conditions, the power flow control in an interconnected power system is now seen as a serious problem, since these flows affect the operation of the electric power system by introducing transmission bottlenecks, congesting flow paths, increasing the transmission losses, decreasing reliability, causing blackouts, and so forth. When an interconnected power transmission system is operated at or close to its power-carrying limits, the transmission congestion restricts additional power to be transferred on that line. In the deregulated power systems, Transmission System Operators (TSOs) have to provide an adequate set of paths from each buyer (seller) to multiple sellers (buyers). A transmission-pricing scheme can be used in order to achieve an efficiently operating market. Pricing unscheduled flows of electricity may also contribute to the reliability of the network. This action may also increase both competition and transmission capacity [4]. The prevention and/or control of unscheduled flows become important due to the need of prevention or reduction of fines associated with the unscheduled flows. The management of these flows requires remedial actions such as switching the lines for new topology and setting transformers and generators for new operating conditions.There have been a number of studies in the literature that attempt to solve the problem. The operating problems caused by unscheduled flows and the complexities of these problems were highlighted in [2]. The paper showed that a stable power flow might circulate in a ring AC power system. Suryanarayanan et al. [3] proposed an approach based on Lp-norms to estimate the unscheduled flows occurring in a wide area-interconnected system. They used contribution factors to assess the participation of generating utilities among causes of unscheduled flow. In [5], a modification to the formula for determining the contribution of generation companies to unscheduled flows in electric power networks was presented. Applications of contribution factor and game theoretic approach for managing unscheduled flows are described in [6]. The effects of parallel paths in system network topology and a survey to explore the present state of practices used to determine the transfer capability issues are well studied in [7]. The authors in [8] proposed an algorithm based on the security constrained optimal power flow for the coordinated control of several phase shifters by one operator with the objective of reducing the unscheduled flow through its system. A method based on power system state estimation was proposed in [9]. The authors presented a monitoring tool to track unscheduled flows due to power transactions not reported to the control area operator. In [10], a genetic algorithm-based procedure was designed for the topological optimization of the network against unscheduled flows. The method considered the status of substation breakers and the location (and angle) of phase-shifting transformers as the control variables.There have been quite a few applications of fuzzy set theory to various power system problems [10, 11]. With the advent of computational power and tools, researchers tried new methods such as genetic algorithms and evolutionary algorithms on power system problems as well [12]. We have proposed a fuzzy-set-theory-based genetic algorithm for the solution of the unscheduled power flow [13]. Since the power systems include thousands of buses resulting in large dimensional matrix operations, the use of parallel environments and techniques will shorten the solutions time. This has enabled us to solve the problem in a multicore environment. Parallel and distributed computing has been applied to speed up certain calculations in power system operation and control. The state of the art on parallel processing applications to power systems was presented in [14]. The authors in [15] presented an efficient parallel GA (EPGA) for the solution of large-scale OPF with practical generators constraints and with consideration shunt FACTS devices.We proposed a fuzzy-set-theory-based parallel genetic algorithm procedure to control and/or prevent unscheduled flows. The problem is formulated as a multiobjective problem subject to operational and electrical constraints and it is handled in fuzzy environment since in practice small variations of power systems variables (bus voltages, line currents, etc.) from their limit values can be tolerated. When the classical methods are used to control loop flow in a rigid manner, a feasible solution may not be found. However, using the proposed method, one can find a feasible solution that maximizes the satisfaction degree. Once the problem is formulated, we solve it using parallel genetic algorithms. The methods using serial GA have restricted applicability in a practical environment because they have a number of problems such as high time complexity. Since the power systems include thousands of buses resulting in large dimensional matrix operations needing the use of iterative solutions methods, the proposed method generates high-quality solutions in a short time. The solution of the problem in a multicore environment also enables us to consider more paths and loops. In the next section, the proposed approach is described. In Sections3–5, the basic principles of genetic algorithms, parallel computation, and parallel genetic algorithms are introduced. In Section 6, the simulation results are given and discussed. Finally conclusions are provided.


## 2. Proposed Approach

In general, a classical multiobjective optimization problem can be formularized combining multiple objectives into a single-objective function by using some user-defined parameters as follow:Minimize(2.1)F=w1⋅f1(x,u)+w2⋅f2(x,u)+⋯+wn⋅fn(x,u)Such that(2.2)g(x,u)=0,h(x,u)≤0,
where x represents system variables, and u represents control variables. The equality constraints are g(x,u)=0, and the inequality constraints are h(x,u)≤0.In the proposed multi-objective optimization model, the minimization of unscheduled line flows is the main objective. The method also tries to keep both active power losses and reactive power losses as small as possible. The constraints on transmission line flows (MVA limits), voltage magnitudes of load buses, and the settings of control variables such as transformer taps, VAR compensators, and generator voltages are inequality constraints. The equality constraints are the real and reactive power balances of the power system. The problem is formulated as follows.


### 2.1. Objectives

#### 2.1.1. The Minimization of Unscheduled Line Flows

The power flow from busi to bus j is(2.3)Pij+jQij=V̇iİij*=Vi(Vi*-Vj*)yij*+Vi2yi0*,
where Pij is the active power fed into the line from bus i and Qij is the reactive power fed into the line from bus j; V̇i is the phasor voltage of ith bus; yij is the admittance of the series elements; yi0 is the half of the line charging admittance.The proposed method minimizes unscheduled real power flows through the contracted paths. This objective is handled in fuzzy environment by using designated fuzzy membership functions.


#### 2.1.2. The Minimization of Power Losses

The active power loss in transmission lines can be calculated as (Figure2)(2.4)PL=∑k=1nlGk(Vi2+Vj2-2ViVjcos⁡⁡(δi-δj)).Figure 2
Transmission line.The reactive power loss in transmission lines can be calculated as(2.5)QL=∑k=1nlBk(Vi2+Vj2-2ViVjcos⁡⁡(δi-δj)),
where Gk is the conductance of the kth line that connects bus i to bus j; Bk is the susceptance of the kth line that connects bus i to bus j;Vi is the voltage magnitude of ith bus; δi is the phase angle of ith bus; nl is the total number of lines.Both objectives, the minimization of real power losses and the minimization of reactive power losses, are handled in fuzzy environment by using fuzzy membership functions.


### 2.2. Constraints

#### 2.2.1. Equality Constraints

The total generated power should be equal to total load demand plus the total losses in the system. The equality constraints are active and reactive power-balanced equations for each generator and load bus. These equations are also known as power flow equations. The equality constraints are nonlinear equations that can be solved by the Newton-Raphson method called the power flow program.The active power injected toith bus in a power system is(2.6)Pi=∑j=1nViVjYijcos⁡⁡(-θij-δj+δi),i=1,2,…,n,
and the reactive power injected to ith bus in a power system is(2.7)Qi=∑j=1nViVjYijsin⁡⁡(-θij-δj+δi),i=1,2,…,n,
where n is the total number of buses, and Ẏij=Yij∠θij is the transfer admittance between bus i and bus j.


#### 2.2.2. Inequality Constraints

The inequality constraints are limits on physical devices in the power system as well as the limits created to ensure system security. These constraints except for control variables limits are fuzzified by using membership functions. The inequality constraints on control variables were considered as crisp constraint during the solution by GA. These crisp constraints are(2.8)Vg,imin⁡≤Vg,i≤Vg,imax⁡,Qg,imin⁡≤Qg,i≤Qg,imax⁡,QC,imin⁡≤QC,i≤QC,imax⁡,Pg,imin⁡≤Pg,i≤Pg,imax⁡,timin⁡⁡≤ti≤timax⁡⁡,
where Vg,i is the generator bus voltage of ith bus, Qg,i is the generator reactive power generation of ith bus, QC,i is the shunt capacitor reactive power generation of ith bus, Pg,i is the generator active power generation of ith bus, and ti is the transformer tap setting of ith bus.The power transfer capacity of a transmission line is limited by thermal rating of the line (MVA rating). To maintain system security, bus voltages have maximum and minimum magnitudes. These limits are handled in fuzzy environment by using fuzzy membership functions.A fuzzy-set-theory-based genetic algorithm procedure was proposed to control and/or prevent unscheduled flows [13]. In fuzzification step, the objectives and constraints were fuzzified by using designated membership functions, where these functions should be properly constructed to guarantee the overall performance. The fuzzy membership functions define the degrees of membership that a crisp value has in a fuzzy set. The membership functions scale both the objectives and the constraints to the unit interval [0,1], so that the satisfactory of the objectives and the constraints can be compared without being influenced by their original values. The closer the membership is to one, the better the solution is for that objective or constraint. In fuzzy decision making, the optimal solution is defined to be the one with the highest degree of membership, and thus the optimization problem becomes that of maximizing the satisfaction with the solution, subject to the crisp and fuzzy constraints [16, 17]. In our approach, the objective functions were minimization of the unscheduled line flows on the contracted paths and minimization of both total active losses and total reactive losses. The constraints (voltages remaining within the limits, line flows remaining within the limits, etc.) were also handled in fuzzy environment by using designated fuzzy membership functions. Since piecewise linear membership functions such as trapezoidal or triangular functions were used in many literatures [11, 18] because of their simplicity and efficiency with respect to computability, they were considered for the objectives and constraints except for the case of the unscheduled flows.The exponential membership function described by four parameters(a,b,c,d) with four breakpoints of the shape, given in Figure 3, has been selected for the fuzzification of the unscheduled flows since an earlier study showed that more satisfying results can be obtained using this kind of membership function [19]. The membership function μL,ij(Pij) which belongs to the MW flow (line flow) through the line between buses i and bus j is defined as(2.9)μL,ij(Pij)={1e(b-a)(ePij-a-1)ifaij<Pij<bij,1ifbij<Pij<cij,e(c-Pij)(d-Pij)ifcij<Pij<dij,0ifotherwise,
where aij<bij<cij<dij must hold.Figure 3
The fuzzy memberships function for the loop MW flows (the exponential form).A membership functionμPL(PL) which is considered for the active power transmission losses is shown in Figure 4, wherePLundesired represents the unsatisfaction limit for the total active power transmission losses. The degree of satisfaction for the objective is zero for any value of PL greater thanPLundesired. We want to keep the objective below  unsatisfaction limit. The system operators taking into account both their experiences and system operation costs can determine these satisfaction limits on the membership function. The membership functions considered for reactive power transmission losses μQL(QL) are similar to μPL(PL).Figure 4
The fuzzy memberships function for the active power losses.The constraints on transmission line flows (MVA limits) are fuzzified using memberships function, which is shown in Figure5. Here, the degree of satisfaction for the MVA flow through the line between bus i and bus j is zero for any value of the MVA flow greater than the line’s current carrying capacity limit (thermal limit) (MVAijunsatisfied). The system operators, who can take into account the feasible operating conditions, can determine the desired MVA flow limit of a specific line MVAijdesired.Figure 5
The fuzzy memberships function for MVA flows.The trapezoidal membership functions shown in Figure6 were used to fuzzify the constraints on bus voltage magnitudes. These membership functions are defined according to the voltage violation of buses to respond to the practical situations as much as possible. In power systems operation, the unsatisfaction tolerance, d, is set to be the allowable maximum bus voltage violation [18].Figure 6
Fuzzy memberships function for bus voltage magnitudes.


## 2.1. Objectives

### 2.1.1. The Minimization of Unscheduled Line Flows

The power flow from busi to bus j is(2.3)Pij+jQij=V̇iİij*=Vi(Vi*-Vj*)yij*+Vi2yi0*,
where Pij is the active power fed into the line from bus i and Qij is the reactive power fed into the line from bus j; V̇i is the phasor voltage of ith bus; yij is the admittance of the series elements; yi0 is the half of the line charging admittance.The proposed method minimizes unscheduled real power flows through the contracted paths. This objective is handled in fuzzy environment by using designated fuzzy membership functions.


### 2.1.2. The Minimization of Power Losses

The active power loss in transmission lines can be calculated as (Figure2)(2.4)PL=∑k=1nlGk(Vi2+Vj2-2ViVjcos⁡⁡(δi-δj)).Figure 2
Transmission line.The reactive power loss in transmission lines can be calculated as(2.5)QL=∑k=1nlBk(Vi2+Vj2-2ViVjcos⁡⁡(δi-δj)),
where Gk is the conductance of the kth line that connects bus i to bus j; Bk is the susceptance of the kth line that connects bus i to bus j;Vi is the voltage magnitude of ith bus; δi is the phase angle of ith bus; nl is the total number of lines.Both objectives, the minimization of real power losses and the minimization of reactive power losses, are handled in fuzzy environment by using fuzzy membership functions.


## 2.1.1. The Minimization of Unscheduled Line Flows

The power flow from busi to bus j is(2.3)Pij+jQij=V̇iİij*=Vi(Vi*-Vj*)yij*+Vi2yi0*,
where Pij is the active power fed into the line from bus i and Qij is the reactive power fed into the line from bus j; V̇i is the phasor voltage of ith bus; yij is the admittance of the series elements; yi0 is the half of the line charging admittance.The proposed method minimizes unscheduled real power flows through the contracted paths. This objective is handled in fuzzy environment by using designated fuzzy membership functions.


## 2.1.2. The Minimization of Power Losses

The active power loss in transmission lines can be calculated as (Figure2)(2.4)PL=∑k=1nlGk(Vi2+Vj2-2ViVjcos⁡⁡(δi-δj)).Figure 2
Transmission line.The reactive power loss in transmission lines can be calculated as(2.5)QL=∑k=1nlBk(Vi2+Vj2-2ViVjcos⁡⁡(δi-δj)),
where Gk is the conductance of the kth line that connects bus i to bus j; Bk is the susceptance of the kth line that connects bus i to bus j;Vi is the voltage magnitude of ith bus; δi is the phase angle of ith bus; nl is the total number of lines.Both objectives, the minimization of real power losses and the minimization of reactive power losses, are handled in fuzzy environment by using fuzzy membership functions.


## 2.2. Constraints

### 2.2.1. Equality Constraints

The total generated power should be equal to total load demand plus the total losses in the system. The equality constraints are active and reactive power-balanced equations for each generator and load bus. These equations are also known as power flow equations. The equality constraints are nonlinear equations that can be solved by the Newton-Raphson method called the power flow program.The active power injected toith bus in a power system is(2.6)Pi=∑j=1nViVjYijcos⁡⁡(-θij-δj+δi),i=1,2,…,n,
and the reactive power injected to ith bus in a power system is(2.7)Qi=∑j=1nViVjYijsin⁡⁡(-θij-δj+δi),i=1,2,…,n,
where n is the total number of buses, and Ẏij=Yij∠θij is the transfer admittance between bus i and bus j.


### 2.2.2. Inequality Constraints

The inequality constraints are limits on physical devices in the power system as well as the limits created to ensure system security. These constraints except for control variables limits are fuzzified by using membership functions. The inequality constraints on control variables were considered as crisp constraint during the solution by GA. These crisp constraints are(2.8)Vg,imin⁡≤Vg,i≤Vg,imax⁡,Qg,imin⁡≤Qg,i≤Qg,imax⁡,QC,imin⁡≤QC,i≤QC,imax⁡,Pg,imin⁡≤Pg,i≤Pg,imax⁡,timin⁡⁡≤ti≤timax⁡⁡,
where Vg,i is the generator bus voltage of ith bus, Qg,i is the generator reactive power generation of ith bus, QC,i is the shunt capacitor reactive power generation of ith bus, Pg,i is the generator active power generation of ith bus, and ti is the transformer tap setting of ith bus.The power transfer capacity of a transmission line is limited by thermal rating of the line (MVA rating). To maintain system security, bus voltages have maximum and minimum magnitudes. These limits are handled in fuzzy environment by using fuzzy membership functions.A fuzzy-set-theory-based genetic algorithm procedure was proposed to control and/or prevent unscheduled flows [13]. In fuzzification step, the objectives and constraints were fuzzified by using designated membership functions, where these functions should be properly constructed to guarantee the overall performance. The fuzzy membership functions define the degrees of membership that a crisp value has in a fuzzy set. The membership functions scale both the objectives and the constraints to the unit interval [0,1], so that the satisfactory of the objectives and the constraints can be compared without being influenced by their original values. The closer the membership is to one, the better the solution is for that objective or constraint. In fuzzy decision making, the optimal solution is defined to be the one with the highest degree of membership, and thus the optimization problem becomes that of maximizing the satisfaction with the solution, subject to the crisp and fuzzy constraints [16, 17]. In our approach, the objective functions were minimization of the unscheduled line flows on the contracted paths and minimization of both total active losses and total reactive losses. The constraints (voltages remaining within the limits, line flows remaining within the limits, etc.) were also handled in fuzzy environment by using designated fuzzy membership functions. Since piecewise linear membership functions such as trapezoidal or triangular functions were used in many literatures [11, 18] because of their simplicity and efficiency with respect to computability, they were considered for the objectives and constraints except for the case of the unscheduled flows.The exponential membership function described by four parameters(a,b,c,d) with four breakpoints of the shape, given in Figure 3, has been selected for the fuzzification of the unscheduled flows since an earlier study showed that more satisfying results can be obtained using this kind of membership function [19]. The membership function μL,ij(Pij) which belongs to the MW flow (line flow) through the line between buses i and bus j is defined as(2.9)μL,ij(Pij)={1e(b-a)(ePij-a-1)ifaij<Pij<bij,1ifbij<Pij<cij,e(c-Pij)(d-Pij)ifcij<Pij<dij,0ifotherwise,
where aij<bij<cij<dij must hold.Figure 3
The fuzzy memberships function for the loop MW flows (the exponential form).A membership functionμPL(PL) which is considered for the active power transmission losses is shown in Figure 4, wherePLundesired represents the unsatisfaction limit for the total active power transmission losses. The degree of satisfaction for the objective is zero for any value of PL greater thanPLundesired. We want to keep the objective below  unsatisfaction limit. The system operators taking into account both their experiences and system operation costs can determine these satisfaction limits on the membership function. The membership functions considered for reactive power transmission losses μQL(QL) are similar to μPL(PL).Figure 4
The fuzzy memberships function for the active power losses.The constraints on transmission line flows (MVA limits) are fuzzified using memberships function, which is shown in Figure5. Here, the degree of satisfaction for the MVA flow through the line between bus i and bus j is zero for any value of the MVA flow greater than the line’s current carrying capacity limit (thermal limit) (MVAijunsatisfied). The system operators, who can take into account the feasible operating conditions, can determine the desired MVA flow limit of a specific line MVAijdesired.Figure 5
The fuzzy memberships function for MVA flows.The trapezoidal membership functions shown in Figure6 were used to fuzzify the constraints on bus voltage magnitudes. These membership functions are defined according to the voltage violation of buses to respond to the practical situations as much as possible. In power systems operation, the unsatisfaction tolerance, d, is set to be the allowable maximum bus voltage violation [18].Figure 6
Fuzzy memberships function for bus voltage magnitudes.


## 2.2.1. Equality Constraints

The total generated power should be equal to total load demand plus the total losses in the system. The equality constraints are active and reactive power-balanced equations for each generator and load bus. These equations are also known as power flow equations. The equality constraints are nonlinear equations that can be solved by the Newton-Raphson method called the power flow program.The active power injected toith bus in a power system is(2.6)Pi=∑j=1nViVjYijcos⁡⁡(-θij-δj+δi),i=1,2,…,n,
and the reactive power injected to ith bus in a power system is(2.7)Qi=∑j=1nViVjYijsin⁡⁡(-θij-δj+δi),i=1,2,…,n,
where n is the total number of buses, and Ẏij=Yij∠θij is the transfer admittance between bus i and bus j.


## 2.2.2. Inequality Constraints

The inequality constraints are limits on physical devices in the power system as well as the limits created to ensure system security. These constraints except for control variables limits are fuzzified by using membership functions. The inequality constraints on control variables were considered as crisp constraint during the solution by GA. These crisp constraints are(2.8)Vg,imin⁡≤Vg,i≤Vg,imax⁡,Qg,imin⁡≤Qg,i≤Qg,imax⁡,QC,imin⁡≤QC,i≤QC,imax⁡,Pg,imin⁡≤Pg,i≤Pg,imax⁡,timin⁡⁡≤ti≤timax⁡⁡,
where Vg,i is the generator bus voltage of ith bus, Qg,i is the generator reactive power generation of ith bus, QC,i is the shunt capacitor reactive power generation of ith bus, Pg,i is the generator active power generation of ith bus, and ti is the transformer tap setting of ith bus.The power transfer capacity of a transmission line is limited by thermal rating of the line (MVA rating). To maintain system security, bus voltages have maximum and minimum magnitudes. These limits are handled in fuzzy environment by using fuzzy membership functions.A fuzzy-set-theory-based genetic algorithm procedure was proposed to control and/or prevent unscheduled flows [13]. In fuzzification step, the objectives and constraints were fuzzified by using designated membership functions, where these functions should be properly constructed to guarantee the overall performance. The fuzzy membership functions define the degrees of membership that a crisp value has in a fuzzy set. The membership functions scale both the objectives and the constraints to the unit interval [0,1], so that the satisfactory of the objectives and the constraints can be compared without being influenced by their original values. The closer the membership is to one, the better the solution is for that objective or constraint. In fuzzy decision making, the optimal solution is defined to be the one with the highest degree of membership, and thus the optimization problem becomes that of maximizing the satisfaction with the solution, subject to the crisp and fuzzy constraints [16, 17]. In our approach, the objective functions were minimization of the unscheduled line flows on the contracted paths and minimization of both total active losses and total reactive losses. The constraints (voltages remaining within the limits, line flows remaining within the limits, etc.) were also handled in fuzzy environment by using designated fuzzy membership functions. Since piecewise linear membership functions such as trapezoidal or triangular functions were used in many literatures [11, 18] because of their simplicity and efficiency with respect to computability, they were considered for the objectives and constraints except for the case of the unscheduled flows.The exponential membership function described by four parameters(a,b,c,d) with four breakpoints of the shape, given in Figure 3, has been selected for the fuzzification of the unscheduled flows since an earlier study showed that more satisfying results can be obtained using this kind of membership function [19]. The membership function μL,ij(Pij) which belongs to the MW flow (line flow) through the line between buses i and bus j is defined as(2.9)μL,ij(Pij)={1e(b-a)(ePij-a-1)ifaij<Pij<bij,1ifbij<Pij<cij,e(c-Pij)(d-Pij)ifcij<Pij<dij,0ifotherwise,
where aij<bij<cij<dij must hold.Figure 3
The fuzzy memberships function for the loop MW flows (the exponential form).A membership functionμPL(PL) which is considered for the active power transmission losses is shown in Figure 4, wherePLundesired represents the unsatisfaction limit for the total active power transmission losses. The degree of satisfaction for the objective is zero for any value of PL greater thanPLundesired. We want to keep the objective below  unsatisfaction limit. The system operators taking into account both their experiences and system operation costs can determine these satisfaction limits on the membership function. The membership functions considered for reactive power transmission losses μQL(QL) are similar to μPL(PL).Figure 4
The fuzzy memberships function for the active power losses.The constraints on transmission line flows (MVA limits) are fuzzified using memberships function, which is shown in Figure5. Here, the degree of satisfaction for the MVA flow through the line between bus i and bus j is zero for any value of the MVA flow greater than the line’s current carrying capacity limit (thermal limit) (MVAijunsatisfied). The system operators, who can take into account the feasible operating conditions, can determine the desired MVA flow limit of a specific line MVAijdesired.Figure 5
The fuzzy memberships function for MVA flows.The trapezoidal membership functions shown in Figure6 were used to fuzzify the constraints on bus voltage magnitudes. These membership functions are defined according to the voltage violation of buses to respond to the practical situations as much as possible. In power systems operation, the unsatisfaction tolerance, d, is set to be the allowable maximum bus voltage violation [18].Figure 6
Fuzzy memberships function for bus voltage magnitudes.


## 3. Genetic Algorithms

Genetic algorithms (GAs) are nature inspired stochastic and population-based search and optimization methods. Their mechanism is based on natural selection and natural genetics. The evaluation function plays the role of the environment that rates solutions according to their fitness, that is, how good they are in solving the search problem for a population [20].GAs work very well on both continuous and discrete combinatorial problems. There are many ways to modify the basic algorithm and many parameters according to a specific problem [21]. They are less likely to get stuck at local optima than gradient search methods.There are three important stages to apply genetic algorithms. They are the following:(1)
definition of the objective function,(2)
definition and implementation of the genetic representation,(3)
implementation of the genetic operators (selection, crossover, mutation).Once these three have been defined, the basic genetic algorithm should work well. To improve performance, many different variations of the operators can be tested [22].Outline and the main operators of the basic GA [23] are as follows.(1)
[Start] Generate random populationP of n chromosomes (population, suitable solutions for the problem).(2)
[Fitness] Evaluate the fitnessf(x) of each chromosome x in the population.(3)
[New population] Create a new population by repeating following steps until the new population is complete.
(i)
[Selection] Select two parent chromosomes from a population according to their fitness (the better fitness, the bigger chance to be selected).(ii)
[Crossover] With a crossover probability, cross over the parents to form a new offspring (children). If no crossover was performed, offspring is an exact copy of parents.(iii)
[Mutation] With a mutation probability, mutate new offspring at each locus (position in chromosome).(iv)
[Accepting] Place new offspring in a new population.(4)
[Replace] Use new generated population for a further run of algorithm.(5)
[Test] If the end condition is satisfied, stop, and return the best solution in current population.(6)
[Loop] Go to step 2.In fuzzy decision making, the optimal solution, which is the fuzzy decisionμD, is given as an intersection of the fuzzy sets describing the constraints (μc(y)) and the objectives (μg(y)):(3.1)μD=min⁡⁡⌊μg(y),μc(y)⌋.
The fuzzy decision μD represents the worst satisfaction value or worst degree of fuzzy membership belongs to either system state variable, constraints, or objectives of the problem. Overall purpose is to maximize the worst degree of membership:(3.2)max⁡μD.
Evaluation function which is called fitness needs to be defined. The fitness function for our problem is defined as follows:(3.3)fitness=11+μD.
In GA formulation, the objective function is the minimization of fitness function:(3.4)min⁡⁡fitness.The control variables such as transformer tap settings, generation bus voltage magnitudes, and generator reactive power generations with their upper and lower limits are used to create chromosomes in GA. An example is given as follows:(3.5)x=⌊ti,Vj,Pk⌋{i=1,2,3,…,l,j=1,2,3,…,ng,k=1,2,3,…,m},
where l is the number of tap changing transformers, ng is the number of generator buses, m is the number of the generators in the system, and xis a possible solution set to the problem.Population consists of many sets of solution. The number of solution sets is defined as population size. Within GA, a load flow analysis is performed for each population. The results of load flow analysis (bus voltages, line flows, losses, etc.) are passed to fuzzy decision making process, where a membership value is assigned for all constraints and objectives. The minimum of those membership values is then tried to be maximized by GA. The process continues up to the end condition defined. One way of defining an end condition is setting a maximum iteration number.


## 4. Parallel Computation

The main disadvantage of GAs is the high CPU time execution, and this can cause heavy computation demands in a large-scale power system case, resulting in large dimensional matrix operations. One way to overcome these computation demands is to use parallel algorithms in multiprocessor computers or in distributed environments.Parallel computation is evolution of serial computation [24]. In general sense problems are partitioned and solved concurrently. This can be in many ways, some of which are instruction level partitioning, data partitioning, task partitioning, and so forth. To solve a computational problem, parallel computing simply involves the simultaneous use of many computing resources. Computing resources come in the form of either a single computer with multiple processors or multiple computers (computing nodes) connected by a fast network.One of the main purposes of using parallel computing is to solve a problem faster or solve a big problem that is otherwise unsolvable in an acceptable time. Parallel computing historically intended to model difficult scientific and engineering problems that exist in the real world [25].In most cases, there are both parallelizable and unparallelizable parts in the problem formulations. The speedup of a program using multiple processors in parallel computing is limited by the time needed for the serial part of the program. There are also several factors that limit the speedup such as communication time for sending messages. The scheduling of tasks on processors and communication between them consume extra CPU resources. Amdahl’s Law says that speedup of a program is limited by the serial part of the computation time of the program [26]:(4.1)speed=1(P/N)+S,
where N represents the number of processor, P represents the parallelization rate of the program, and S represents the serial computation rate of the program. Theoretically, if N goes to an infinitely large number of processor, computational time spent for the parallel part of the program goes to zero. Then the maximum speed will be(4.2)speedN→∞=1S.


## 5. Parallel Genetic Algorithms

GAs are very powerful search and optimization techniques that can be applied to problems in many different areas such as sciences, business, and engineering [27]. GAs simulate the evolution of a population of individuals, which are possible solutions of a specific problem, through the evolution principle. They can produce an acceptable solution to a problem in a reasonable time. However when they are applied to very hard and large problems, it takes quit a long time. In that case, there might be some effort to reduce the time into a reasonable range. One way of reducing execution time is to implement GAs in parallel environment. These topics are well discussed in literature [25, 27–31].There are some important questions to be answered when GAs are parallelized. These are the following.(i)
How fitness computation and mutation would be done?(ii)
How would selection operator be applied, local or global?(iii)
Would there be a single or multiple populations?(iv)
If multiple populations were used, how would individuals be exchanged?There are four major types of parallel GAs, namely, master-slave parallelization, coarse-grained parallel GA, fine-grained parallel GA, and hybrid parallel GA [30, 31]. The first type, which is one of the more popular types in practice, was used in this study. If function evaluation is expensive, master-slave implementation (global parallel GA) becomes more efficient. This type uses a single population. This is usually implemented using master-slave programs, where the master stores the population and applies the GA operators which are selection, crossover, and mutation. The slaves evaluate the fitness for the assigned individuals to them. Exploration of the search space is in the same way of serial GA. That is why it is called as global parallel GA. It is also known as master-slave model or distributed fitness evaluation. Master and slave processes are presented in Figure 7.Figure 7
Master-slave parallel GA.The execution time consists of two parts. One is computation time for evaluation of the fitness of a part of the population assigned to each individual slave processor. The other is the communication time in between master and slaves. The master-slave parallel GA is particularly easy to implement. The flow chart of the proposed procedure is shown in Figure8. Master and slaves are also presented in this figure.Figure 8
Flow chart of the proposed procedure.


## 6. Simulation Results

The simulations are performed using the utilities provided by the high-performance computing lab in Computational Science and Engineering Department in Istanbul Technical University. The computing system specifications are given in Table1.Table 1
The computing system specifications.
Computing systemDistributed cluster, with 37 nodes, each consisting of 2 CPU’s having3.40GHz CPU, and 2GB memoryNetwork typeMyrinet-2000Operating systemRHEL 4.4Resource management and scheduling softwareLSF (load sharing facility)LibraryMPICH-GM Myrinet 1.2.7The algorithm was implemented in MATLAB environment using toolboxes such as Genetic Algorithm Toolbox, Parallel Computing Toolbox, and Distributed Computing Server, and MATPOWER which is an open source power flow program was also used [32, 33]. Master-slave parallelization method was chosen among the types of parallel GAs. When a problem demands higher computation time, it fits better [34].The method is applied to IEEE 300-bus test system having 69 generators and 411 lines. The detailed data of the system can be found in literature. Three lines given in Table2 are considered as the contracted paths, and the scheduled flows for the lines are also given in the table. In the simulations, d=0.05 (see Figure 6) is taken in fuzzification of the voltage magnitudes. PLunsatified and QLunsatified are taken as 1.8 times of active and reactive power losses of base case, respectively. The ranges are taken large so that they would not be the dominative ones among the other objectives and constraints. Generator reactive limits are also taken into account. The parameter values used in GA solution are given in Table 3. The initial population is taken from the base case power flow solution of the system. Simulations are performed for two different cases. In the first case, population size and iteration number are taken as 100. In the second case, both population size and iteration number are taken as 150.Table 2
The contracted paths and their scheduled flows.
Path 1Path 2Path 3Line number4448117From bus3762To bus150131144Scheduled flow [MW]10060100(a,b,c,d)(85,95,105,115)(40,55,65,80)(90,95,105,110)Table 3
The parameter values used in GA.
Case 1Case 2The population size100150Iteration number100150The population typeReal valued (double vector),The selection typetournament (size is 4),chromosome length231The elite number2The crossover typeTwo point crossoversThe crossover rate0.8The mutation functionAdapt-feasibleWhen evaluating a parallel system, it is important to know how much performance gain is achieved by parallelizing a given application over a serial implementation [34]. In the first case, the rate of parallelization of our method was figured out as 78% after a series of simulations. We applied Amdahl’s law to obtain the theoretical maximum speedup. In this case, P is 0.78, S is 0.22, and from (2.7), the theoretical maximum speedup can be calculated as 4.5455. In other words, the maximum acceleration from the distributed computing for our program will be 4.5455 times that of 100% serial simulation.When a distributed computing system is used, the speed of program can be calculated by using the following equation:(6.1)speed=T1TN,
where N is the number of processors, T1is the execution time for running the program completely in serial manner, and TN is the execution time for running the program by using the N processors.Efficiency is a measure of the fraction of time for which a processing element is usefully employed; it is defined as the ratio of speedup to the number of processing elements. In an ideal parallel system, speedup is equal toN and efficiency is equal to one. In practice, speedup is less than N and efficiency is between zero and one, depending on the effectiveness with which the processing elements are utilized [34]. The efficiency of the system used for our program can be calculated from (6.2):(6.2)efficiency=(speedN)100%.As seen in Figure9, there is no improvement in terms of the execution time after the usage of fourteen processors even though we continue increasing the number of processor.Figure 9
Execution time with different numbers of processors for Case 1.In Figure10, we see the change of speed versus processor number. From this figure, it is easy to find that the speed can be increased greatly by using parallel implementation. But the implementation achieved low efficiency when the numbers of processor were increased.Figure 10
Speedup with different numbers of processors for Case 1.Figure12 shows the execution time for Case 2. The measured times were wall-clock time (elapsed time). As seen from Figures 9 and 12, computation times fluctuate. This was due to changes in the usage of the network at the moment of the program execution. As shown in Figure 13, reasonable speedup was obtained when the number of processors was increased. In Figure 14, the change of efficiency versus processor number can be seen. Even though we continue to use more processors, the efficiency of the program is not increasing due to not utilization of the full capacity of each processor. In the second case, the rate of parallelization of our program went up to 90%. In this case, the theoretical maximum speed obtained using Amdahl’s law will be 10 times that of 100% of serial simulation. For this case the figures for time, speed, and efficiency are given in Figures 12, 13, and 14, respectively.The system solutions for the objectives and the constraints are given in Tables4 and 5 for the Case 1 and Tables 6 and 7 for the Case 2. When we compare the base case results to the simulation results, we see that our method is successful. As seen from Tables 4 and 6, we reached the targeted values in the given ranges. This is better understood from the averaged fuzzy membership values given in Tables 5 and 7.Table 4
Averaged values for objectives and constraints (population size is 100, iteration number is 100).
Path 1 [MW]Path 2 [MW]Path 3 [MW]Ploss [MW]Qloss [MVAR]Base case92.9046.04−101.16408.3165504.18Target10060−100Test results101.153354.0599−97.4213412.25465526.4Best fitness0.7741Table 5
Averaged fuzzy membership values for objectives and constraints (population size is 100, iteration number is 100).
min(μ)μPath 1μPath 2μPath3μPlossμQlossBase case0.00010.11630.00011.0000.44510.4451Test results0.37500.82550.43820.83500.43910.4422Table 6
Averaged values for objectives and constraints (population size is 150, iteration number is 150).
Path 1 [MW]Path 2 [MW]Path 3 [MW]Ploss [MW]Qloss [MVAR]Base case92.9046.04−101.16408.3165504.18Target10060−100Test results103.379054.8188−95.7833405.30455492.9Best fitness0.7741Table 7
Averaged fuzzy membership values for objectives and constraints (population size is 150, iteration number is 150).
min(μ)μPath 1μPath 2μPath 3μPathμQlossBase case0.00010.11630.00011.0000.44510.4451Test results0.32020.97920.65440.56320.44850.4456As seen from Table5, the smallest membership value for the base case is improved from the level of 0.0001 to 0.3750. The membership values for the lines where  we try to control power flow are improved from 0.1163 to 0.8255 for the first path and from 0.0001 to 0.4382 for the second path. The membership value for the third path is decreased from 1.000 to 0.8350. Only decrease is occurred for the third path but this decrease is not significant when we consider other improvements.There is no bad effect on the membership values of both active and reactive power losses. This situation was among the objectives of our problem. During the control of power flows, active and reactive power losses should not increase significantly.Similar improvements are seen for the second case in which population size and iteration number are increased. These are seen from Tables6 and 7.The increment on iteration number and population size has an effect only on the total computation time. This is seen on Figures9 and 12. The efficiency for the population size 100 is getting below 40% after utilizing eight processors as seen from Figure 11. The efficiency for the population size 150 is getting below 40% after utilizing fourteen processors as seen from Figure 14. The efficiency value below 40% is not recommended.Figure 11
Efficiency with different numbers of processors for Case 1.Figure 12
Execution time with different numbers of processors for Case 2.Figure 13
Speedup with different numbers of processors for Case 2.Figure 14
Efficiency with different numbers of processors for Case 2.


## 7. Conclusion

This paper presents a multiobjective optimization approach based on both genetic algorithms and fuzzy decision making to manage unscheduled flows. The remedial actions for prevention and/or control of unscheduled flows can be determined by the proposed method. Solving the optimization problem with serial genetic algorithm may be time-consuming when the number of system states is large, since power flow analysis is performed for each system state (each population). Using supercomputing facilities and parallel computing techniques together decreases the computation time greatly.To speed up the computing process, parallel implementation of genetic algorithm is proposed in this study. As seen from the test results, when the population size is increased, the amount of computation job per processor is increased as well. In that case, each processor is utilized better. Some improvement on efficiency is also seen. The large portion of the computation time is spent for the function evolution in the problem. This is due to the nature of our problem. Master-slave PGA is perfectly matched with our needs.We conclude that by using a larger system, we would gain more benefit from using a high-performance computing facility that allows parallel implementation of our program to get the system solution faster than the serial corresponding.



---

*Source: 376291-2012-04-09.xml*