Frequency regulation service in energy market is one of the profitable ancillary services. Its main objectives are to keep the tie-line power exchanges within the scheduled limits and to maintain zero change in frequency. Since frequency control in a complex and competitive environment is a tedious process, the control technique must be more intelligent and adaptive for a changing environment. More research work has been carried on Frequency control in restructured power system [1] - [7] . Possible issues in Frequency control in energy market were discussed in [1] . Vertically integrated structure and bilateral based scheme modeling were discussed in [2] [3] . Different approaches like pluralistic and hierarchical methods were discussed and quoted in [4] . A decentralized approach based on control theory was formulated and analyzed in [5] . Some researchers have proposed a frequency control logic based on fuzzy logic controller [6] , reinforcement learning [7] and artificial neural network [8] in interconnected power systems.

The Load Frequency Control (LFC) issue as an ancillary service represents an important role to maintain an acceptable level of efficiency, quality, and reliability in a deregulated power system environment. Many researchers have started to analyze possible new LFC schemes and regulation solutions, with paradigms suited for the energy market scenarios [4] - [9] . In this paper bilateral contract scheme under restructured power system is taken for analysis. The objective of the proposed controller in bilateral market is to reduce damping of frequency and tie-line power deviations oscillation for all admissible uncertainties and load changes. So, three different testing conditions are taken to test the robustness of the controllers.

This paper is organized as follows: Generalized bilateral scheme is presented in section 2 and Generation Participation Matrix is defined in this section. System taken for investigation is defined in section 3. Fuzzy rule based control scheme for a bilateral based restructured power system is presented in section 4. Automatic membership function tuning using PSO based controller is discussed in section 5. Simulation results for different cases are presented in Section 6. The objective of this work is to formulate the dynamic LFC and evaluate and analyze the LFC based on different performance measures under parametric uncertainties.

Energy market scenarios are adapting changes every day to achieve the goal of better performance and efficiency in restructured power systems. In this paper, bilateral based market is taken for analysis. A general configuration for the LFC [2] [3] in a deregulated environment is shown in Figure 1.

Bilateral based deregulated environment consist of an Independent System Operator (ISO), Distribution Companies (DISCOs), Generation Companies (GENCOs), and Transmission Companies (TRANSCOs). In the new environment, DISCOs may contract power from any GENCOs and ISO have to supervise these contracts

[9] . The “Generation Participation Matrix (GPM)” shows the participation factor of each Genco in the considered control areas, and each control area is determined by a Disco. The rows of a GPM correspond to Gencos and the columns to the control areas that contract power. For example, the GPM for a large scale power system with m control areas (Discos) and n Gencos, has the following structure, in which gpf_ij refers to ‘‘generation participation factor’’ and shows the participation factor of Genco i in total load following the requirement of DISCO j based on a specified bilateral contract [10] in Equation (1).

New information signals due to possible various contracts between Disco i and other Discos and Gencos are shown as dashed arrows in Figure 1.

where

where ΔP_di (Figure 1) is the area load disturbance, ΔP_Loc-i is the contracted load demand (contracted and uncontracted) in area i, and ΔP_tie-I,actual is the actual tie-line power in area i. Using Equation (2), the scheduled tie-line power (ΔP_{tie-I,scheduled}) can be calculated using Equation (7).

In the bilateral LFC structure, Control is highly decentralized. Each Load Matching Contract requires a separate control process, yet these control processes must cooperatively interact to maintain system frequency and minimize the area control error [4] , where

X_i Generator participation factor under contract.

Y_i Difference of import and export power.

Z_i Tie-line coefficient.

∆f_i frequency deviation.

B_i frequency bias.

∆P_{tie_i} net tie-line power flow.

∆P_di area load disturbance between areas i and j.

α area control error (ACE) participation factor.

∆P_Li contracted demand of area i.

∆P_mi power generation of Genco i.

∆P_{Loc_i} total local demand (contracted and uncontracted) in area i.

∆P_{tie_i,actual} actual ∆P_tie__i.

Consider the load frequency control problem for a two-area power system as shown in Figure 2. Each control area has the structure as shown in Figure 1. The control area considered here is simplified power systems consisting of two generating units of thermal reheat system [11] . According to [12] , the collective performance of all generators in the system is the interesting part in the analysis of LFC. Genco’s considered here are reheat thermal units. The inter-machine oscillations are not considered. It is assumed that the response of all generators to changes in system load are coherent and can be represented by an equivalent generator, which has an inertia constant and a damping constant equal to the sum of the inertia constants and damping constants of all the generating units.

The model is simple but captures the essential dynamics of a power system and has been widely used for LFC design purpose. The load frequency control problem for a multi-area power system requires that not only the frequency deviation of each area must return to its nominal value but also the tie-line power flows must return to their scheduled values. So a composite variable, the area control error (ACE), is used as the feedback variable to ensure the two objectives.

For Area #i, the area control error is defined as

and the feedback control for Area #i take the form

The transfer functions of the governor, the turbine, and the rotor inertia and load for Area #i are denoted by G_gi(s), G_ti(s), and G_pi(s), respectively. The transfer functions of the above mentioned are represented as

then the transfer function from u_i to ∆f_i can be easily found as

The transfer function for Area #i.

The system parameters taken for investigation [13] are given in Appendix A.

Due to the increasing complexity and changing structure of power systems, it needs an intelligent controller as an alternative to the conventional control methodology. Fuzzy Logic Control would serve the better role in load frequency control in restructured power system. Fuzzy controller consists of three main stages, namely the Fuzzification interface, the Inference rules engine and Defuzzification interface.

For LFC, ACE and its derivative are chosen as inputs and Proportional and Integral Gains are taken as control outputs of the fuzzy controller. Gaussian Membership Functions (MFs) is chosen for obtaining fast response from the system. The MFs were named NL (Negative Large), NM (Negative Medium), NS (Negative Small), Z (Zero), PS (Positive Small), PM (Positive Medium) and PL (Positive Large) followed by error(E) and Derivative error (DE). In [14] , Ziegler-Nichols proposed a heuristic method to quickly adjust the controller parameters P, PI, PID. The MFs ranges of the inputs (Ace & dAce) and control outputs (K_P & K_i) for area i are determined based on the experiments with classical (Ziegler-Nichols Tuned) PI controller. These two input signals are used as rule-antecedent in the formation of rule base, and the control outputs are used to represent the contents of the rule-consequent in performing the rule base [7] .

Each fuzzy rule is in the form,

where MF_nj is a membership function and y_n and z_n are singleton values.

where C_nj is mean and σ_nj is standard deviation of the membership function.

AND operation is used and is given by

The output K_P and K_I calculated by

“k” refers total number of rules.

Based on the knowledge gathered from classical (Z-N) controller, fuzzy rules are formed and the appropriate rules for K_i and K_p are given in Table 1.

PSO is a novel population based meta heuristic, which utilize the swarm intelligence generated by the cooperation and competition between the particles in a swarm and has emerged as a useful tool for engineering optimization. It has also been found to be robust in solving problems featuring nonlinearity, non-differentiability and high dimensionality.

In PSO, each particle is flown through the multidimensional search space, adjusting its position in search space according to its own experience or knowledge and that of neighboring particles. Therefore, a particle makes use of the best position encountered by itself and the best position of its neighbors to position itself toward an optimum solution. The effect is that particles fly toward the global minimum, while still searching a wide area around the best solution. The performance of each particle is measured according to a predefined fitness function which is related to the problem being solved [15] . The upper limit of the random value is a system parameter specified by the user. The larger the upper limit, the more the trajectory of the particles oscillates. Smaller values ensure smooth trajectories. This upper limit prevents particles from moving too rapidly from one region in search space to another. Vmax is usually initialized as a function of the range of the problem. The inertia weight controls the influence of previous velocities on the new velocity. Large inertia weights cause larger exploration of the search space, while smaller inertia weights focus the search on a smaller region [16] . In this proposed PSO, the inertia weight (w) is started with a large value, which is decreased over time.

PSO starts with a population of random solutions “particles” in a D-dimension space. The ith particle is represented by. Each particle keeps track of its coordinates in hyperspace, which are as-

sociated with the fittest solution it has achieved so far. The value of the fitness for particle i (pbest) is also stored

as. The global version of PSO keeps track of the overall best value (gbest), and its location, obtained thus far by any particle in the population. PSO consists of, at each step, changing the velocity of each particle toward its pbest and gbest according to Equation (15). The velocity of particle i is represented as. Acceleration is weighed by a random term, with separate random numbers being generated

for acceleration toward p best and g best. The position of the ith particle is then updated according to Equation (16) [17] . Figure 3 shows the flowchart of the proposed PSO algorithm.

It should be noted that choice properly fitness function is very important in synthesis procedure. Because different fitness functions promote different PSO behaviors, which generate fitness value providing a performance measure of the problem considered [1] . Integral of Squared Error (ISE), criterion integrates the square of the error over time. ISE will penalize large errors more than smaller ones (since the square of a large error will be much bigger). ISE is used as fitness function in the proposed PSO.

In order to illustrate the behavior of the proposed control strategy some simulations has been carried out. The system parameter of the test system is given in Appendix A.

Consider a system where all GENCOs in each area participate in LFC, i.e. ACE participation factors are,;,

Case I: In this case it is assumed that all Gencos are participating in the LFC task as per contract. It is assumed that a large step load 0.05 p.u MW is demanded by each DISCOs in all areas. Frequency response for

both areas, tie line power exchange are shown in Figures 4-6. The actual power generated by each Gencos reach the desired value is given in Table 2. The membership function parameters Mean (C_nj) and standard Deviation (σ_nj) before and after updation are given in Table 3.

Case II: It is assumed that the rotating mass and load pattern parameters D_i and M_i have uncertain values in each control area. In addition to case I, 25% decrease in parametric uncertainties is taken. The results are shown in Figures 7-9. Simulation results show that all the Hybrid Fuzzy controllers track the load fluctuations possibly good. Their overshoot & settling time is better than conventional Fuzzy controllers.

Case III: Consider Case II, in addition to the specified contracted load demand and 25% decrease in D_i and M_i, a bounded random step load change (∆P_di) as an uncontracted demand appears in each control area.

The main aim of this test is to check the robustness of the proposed controller against uncertainties and random load disturbance. The load pattern applied to areas is shown in Figure 10. The results are shown in Figure 11 and Figure 12. The performance measures like overshoot and settling time are compared in Table 4.

The optimum values of the Proportional and Integral Gains of the different controllers for the aforementioned cases are shown in Table 5.

In this paper, automatic tuning of membership function of fuzzy controller using PSO is applied to bilateral LFC scheme. PSO is a popular search algorithm; it utilizes the swarm intelligence of cooperation and fast convergence to find the optimum value of the parameters of membership function used in fuzzy logic controller. Simulation results show the effectiveness of the proposed controller in damping the frequency oscillations and tie-line power very fast with less undershoot and overshoot. It is also seen from simulation results that under normal circumstances, Fuzzy and Hybrid Fuzzy using PSO-based controllers produce the better results. During parameter uncertainties and large load disturbance, Hybrid Fuzzy penalizes excessively over damped and excessively under damped responses and reduces overshoots. It has been therefore shown that Hybrid Fuzzy controller produces better performance in all testing environments.

P. Anitha,P. Subburaj, (2016) Hybrid Fuzzy Controller Based Frequency Regulation in Restructured Power System. Circuits and Systems,07,759-770. doi: 10.4236/cs.2016.76065

;;;;;;;; D₁ = 0.015 pu×MW/Hz; D₂ = 0.016 pu×MW/Hz; M₁ = 0.16667; M₂ = 0.2012; B₁ = 0.3483; B₂ = 0.3827; T₁₂ = 0.20.