Introduction lectric power systems have long been perceived and exploited like national entities or regional areas. The interconnections between such zones being mainly used with the aim of help in case of failure in a nearby zone. In the new context of the electricity sector, this situation has changed since the interconnections are increasingly used in order to maximize the exchange of electricity through power exchanges recently introduced. This leads the transmission system operators (TSOs) to operate their systems increasingly close to their limits. In an individual way, each TSO is obliged to re-examine its way of doing the usual studies on its own network. In particular, a more refined modeling of the influence of neighboring systems on its own should be considered. Moreover, the process of extension of the synchronous area which continues today, led to a significant increase in the size of dynamical systems to study. Besides the difficulty posed by the size, new structural phenomena occur from this extension as the electromechanical oscillations of low frequency observed between generators remote network. Many dynamic equivalence techniques have been developed over the years [1]- [9]. More research focused on coherency equivalents and modal equivalents nowadays, however, they need detailed data of the external system which may be difficult to get in the power market environment. Some estimation methods are based on available information from the boundaries nodes and do not require any knowledge detail of the external network. The classical estimation methods are mainly based on the linearization of the system around an operating point with theoretical constraints validity [10] - [12]. Consequently, they are limited by their validity in nonlinear practical applications, offline and online applications. This paper proposes a new evolutionary approach for estimating parameters to determine a dynamic equivalent model of an external power system from synchronized measurements of disturbance occurred in internal system, which obtained by PMUs (Phasor Measurement Units). This approach was applied to 10 machines 39 buses New England test system [13]. Comparing the properties of the equivalent system with ones of the original system, the result shows that the reduced system can represent dynamic behaviors of the original system well, for any kind of disturbance. # II. # Dynamic Equivalents The equivalent is represented by a model with unknown parameters. In operation, the power system is often perturbed by small random disturbances. The estimation of unknown parameters of the equivalent model is operated by comparing the measurements obtained from the real network, and the same measurements made on the equivalent network by minimizing some objective function. # a) Dynamic equivalent model The basic requirement of the dynamic equivalents is that the response of the equivalent system can fit the original system dynamically and approximately when faults happen in the internal system. In our equivalent system, equivalent generator represents dynamic effect of the external system on the internal system located at the boundary bus (Fig. 1). When there are multiple boundary buses connecting the external system with the internal system, multiple generators are adopted to represent dynamic influence of the external system on the internal system, with an equivalent generator located at each boundary bus [14]. Parameters to be estimated are: ????,????,?????,??,??,???????, voltage regulator gain ?? ?? and its time constant ?? ?? . The parameters identification for nonlinear systems includes deterministic and evolutionary methods. Evolutionary methods offer greater ease of adjustment problems and other advantages over deterministic methods because they: -Require only the calculation of the objective function, without this last being forced to be continuous or differentiable. -Easily adaptable to multi-objective problems optimization, and complex systems with very important unknown number of parameter. -Have great theoretical probability to find the global optimum. -Offer great simplicity of implementation. An evolutionary approach based on genetic algorithms is presented in this document, to identify the dynamic equivalent model parameters, by minimizing the range of speed variation, electrical power variation and terminal voltage variation produced in all generators of the internal system between the original and the equivalent system. # i. Error function The original system equations are generally written: (5) Where: A= State matrix (4m× 4m). B = Control matrix (4m×2m). C = Output matrix (2m× 4m). ? Î?"?? Absolute angle rotor variation (rad. elec/sec). ? Î?"ð??"ð??" Angular velocity rotor variation (pu). ? Î?"????? Internal voltage variation (pu). ? Î?"??ð??"ð??"?? Internal excitation voltage variation (pu). ? ?????? Mechanical torque variation (pu). ? ?????? Order excitation systems variation (pu). ? Î?"???? Active power variation (pu). ? ?????? Change in terminal voltage generators (pu). The equations describing the dynamic equivalent system are expressed similar to those of the original system ( 5), but with a simplified structure and less parameters. (6) are all functions of ??, the equivalent system parameter vector to be identified. Therefore, an error function may be defined by: (7) Where Î?"?? represents the responses of the original system, which are directly measureable. And Î?" (??), must be calculated by simulating the equivalent system with the same disturbance. # ii. Objective function The idea is to minimize the error function (7) between the measure and the model output for all machines belonging to the internal network. The mathematical model of the optimization problem is a multi-objective function and can be formulated as follows: (8) Where: By weighting the measurements by a weighting factor ????, our objective function becomes: (9) Where: ?????????= number of machines belonging to the study system (internal system). iii. Encoding and initial population The encoding of individuals is an important parameter in population research methods. These are represented as a strings (chromosomes) containing characters or genes of a predetermined alphabet. There are different ways to code a solution. In our study, the individual is represented by eight parts of chromosome corresponding to the eight parameters to be estimated, each gene (parameter) is represented by its physical value, which means that, the real coding is adopted (fig. 2). The GAs requires an initial population to start search process. Applied methods generate randomly a set of solutions belonging to the following area: (10) Where: The individuals' number in the initial population is chosen such that: To form a group of (Local parents), in which each parent "i" satisfied the fitness function minimum: (13) Where: -A number of individual « ???????? » ; To form a group of (Global parents), in which each parent "k" satisfied the fitness function minimum: (14) Where: The local parents group includes the favorable solutions for each type of measure that corresponds to each machine belonging to the internal system. However the global parents group includes the favorable solutions for each type of measure that corresponds to all the machines belonging to the internal system. The population called, "Parents population" is made up of the two groups (Global parents & Local parents). # b. Multi-parent recombination Recombination is used mainly in evolutionary strategies. Contrary to k points crossover operators, which exchange information between two parents, the recombination creates the descendants, by weighting many parents components'. We define three weighting operators: Where: ??h?? ? Integer belonging to [1,8], which represents crossed chromosome order. ???? 1 ,???? 2 : Two individuals of local parents group. We generate ???????? individuals for construct the new population (generation) from the two groups' "local parents" and "global parents", with the following steps: -we choose a random pair of individuals from the local parents group. -a crossed chromosome is chosen randomly « chr ». -a random real number, "h", is generated between 0 and 1, to select the weighting operator (Fig. 4). (15) -we move to the next individual "i +1" when the number of weighting point "k" of the individual "i" is complete. 3. Accordance with the procedure for selection of parents, the best individuals evaluated by the fitness function (12) between parents population and this new population, construct the new parents population and clear the old parents population for the next generation (elitist strategy). III. # Algorithm The complete algorithm of the proposed method is given below: Step1 : Introduction of static and dynamic data of the complete system. Step2 : Power flow calculation, linearization, disturbance choice and network dynamic simulation. Step3 : Border node choice and equivalent system linear model construction from results of the power flow. Step4 : Generating the initial population according to the procedure specified in subsection (III.B.3). Step5 : Run equivalent system dynamic simulation and evaluate all individuals by fitness function (12). Step6 : Construct the population of parents following the procedure described in paragraph (III.B.5.a), with equations (13,14). Step7 : While the number of generation has not reached the maximum number ?? ð??"ð??"???? : Step7.1 : Generate the new population using described procedure in paragraph (III.B.5.b). Step7.2 : Run equivalent system dynamic simulation and evaluate all individuals by fitness function (12). # Application To validate the proposed approach, developed algorithm was used to build the dynamic equivalent of transmission network IEEE New-England 39 bus [13]. It represents a simplified New England transmission network (northeastern United States) that is part of real U.S. network. This network consists of 10 generators (PGtotal= 6.19 GW, QGtotal= 1.28 GVAr) and 39 bus with 19 load bus. The original system structure is as Fig. 5. new population (Fig. 3). In Fig. 5, the shaded part is the internal system and the rest is the external system to be reduced. There are a boundary bus 16 and two tie lines 16-15 and 16-17 between the internal and the external systems. Fig. 6 is the sketch of the equivalent system. II. These parameters were obtained after several tests with an appropriate adjustment. # ? > 0 # Table II : Control Parameter Values To allow a better parameters estimation, different disturbances were applied in different parts of the internal system. Each proposed disturbance has a period of 10s. An event is generated each time to the first second of the simulation. The algorithm verifies the equivalents behavior in period of 10s, by comparing its response with the response of the original system. If the equivalent response exceeds some level relatively to the original systems response in a new disturbance, the equivalent generator parameters are recalculated and updated. Each solution obtained by applying the previous disturbances must be part of the initial population of the next estimation procedure. The estimated parameters of the equivalent generator are presented in Table III. # b) Validation A small perturbation represented by disconnecting the line 3-18 is applied to bus18 of the equivalent system and the original system respectively at time, t=2s and lasted for about 12 cycles. Time domain dynamic simulations were performed both on the two systems. In order to verify the performance of the proposed dynamic equivalent approach in the case of the large disturbances, three-phase fault was applied to bus 30 of the equivalent system and the original system respectively at time t=2s, and the fault lasted for about 3 cycles. Time domain dynamic simulations were carried out on the original system and the equivalent system, and dynamic responses are shown in fig. 10 The simulation results of equivalent system with estimated parameters are encouraging and show the effectiveness of the developed algorithm and the ability of the estimated equivalent to reproduce the influence of the external system on internal system for small and severe disturbances. V. # Conclusion In this study, we proposed a nonlinear estimation method for parameters, based on an evolutionary algorithm. In order to obtain the external system dynamic equivalent, several dynamic simulation disturbances were applied to increase the accuracy of the model estimated by the developed algorithm. The proposed approach does not require data of external system. This approach requires only configuration information, settings and operating status of the internal system. Dynamic simulation was performed on both the original system and the equivalent system under different operating conditions, and the results show that the obtained equivalent system can represent the main dynamic characteristics of the original system well. Thus, the approach proposed is proved to be feasible and has potential for tackling the complex practical application in power system. # Global Journal of Researches in Engineering 1![Fig. 1.a : Original System Fig. 1.b : Equivalent System The equivalent generator model can be expressed by the following equations: (1) (2) (3) (4)](image-2.png "Fig. 1 .") ![Estimation of External Power System Dynamic Equivalent Parameters Global Journal of Researches in Engineering Volume ? ?? ? ? ? + ? ? ? ?? ? ? ? ? ? ? ?(?ð??" ? 1)) ? = ð??" 0 (?ð??" ? 1) ? ? ? ? ?? ? ? ? ? ?? ? ) ? ð??"? = 1 ? ? (? ? ? ??ð??" ? ? ? + ? ? ? ? ð??"? ) Î?"? = ??? + ??? Î?"? = ??? with : ?? = ?ð??" 1 ? ?ð??" ? ?? 1 ? ?? ? Î?"? ? ?1 ? Î?"? ? ?? Î?"? ð??"?1 ? Î?"? ð??"?? ? ?? = ?? ? 1 ? ?? ?? Î?"? ?1 ? Î?"? ?? ? ?? = ?ð??" 1 ? ?ð??" ? Î?"?? 1 ? Î?"?? ? Î?"? ?1 ? Î?"? ?? ? Î?"? ? = ? ? ?? ? + ??? Î?"? ? = ? ? ?? ? ? ? , ? ? , ?? ? ?? Î?"? ? , ? ? = ?? ? ?? ? ? ??? ? ? = ? 1 ? , ? 2 ? , ? , ? 3 ? ???? ? ? > 0 ? ? ? = ?? ? ? ? ?? ? ?, ? ? ð??"?? ?=?? ; ? = 1,2, ? ,3 ? ? ??? ??? ? ? = ? ? ?? ? ? ? ?? ? ?, ? ? ð??"?? ?=?? 3 ? ???? ?=1 ? > 0 ? = [? ? , ? ? , ? ? ? , ?, ?, ? ? ?? , ? ? , ? ? ] ?? = number of machines belonging to the complete system.](image-3.png "Nonlinear") ![,????,????)=3. ??= Total number of machines belonging to the complete network.??????= Number of machines belonging to the external system.A representation of Nind individuals (solutions) in an initial population is shown in fig.2.](image-4.png "") 2![Fig. 2 : Organization and coding of individuals in initial population. iv. Evaluation mechanism With a fitness function or adaptive function the evaluation is done in a closed interval [?? ?? ,??ð??"ð??" ???? ] using each time a new dynamic simulation of the equivalent system for each population individuals. (12)](image-5.png "Fig. 2 :") ![?ð??" ? ???? Î?"?? 1 ? Î?"?? ? ???? Î?"? ?1 ? Î?"? ?(? ???? ) ? ?? = ?ð??" 1 ? ?ð??" ? ???? Î?"?? 1 ? Î?"?? ? ???? Î?"? ?1 ? Î?"? ?(? ???? ) 2, ? , ? ??? * ? ? ??? ?????? ? = min ? =1..? ??? ð??" ? (??? ? ) ? = 1,2, ? , ? ??? * ? ? ??? ?????? ? = min ? =1..? ??? ð??" ? (??? ? ) ? ??? * ? ???? ?=1 ? = 1,2, ? , ? ??? a.](image-6.png "") 13![Fig. 3 : New population from the group Global Parents 2. For each individual "I " belonging to the new population: -A random integer number of weighted point "k" is generated between 1 and 8 * ?? ?????? . -In each weightings j=1,2?,k :-we choose a random pair of individuals from the local parents group.-a crossed chromosome is chosen randomly « chr ».-a random real number, "h", is generated between 0 and 1, to select the weighting operator (Fig.4).](image-7.png "1 .Fig. 3 :") 4![Fig. 4 : Recombination operator distribution's -the value obtained by the weighting operator is assigned to the individual's chromosome "i(chr)" of the new population, if this value satisfies the constraint:](image-8.png "Fig. 4 :") 3![Select the new individuals in parents population following described procedure in paragraph (III.B.5.b.3). Step8 : End while. Step9 : Select best individuals in parents population with the objective function (9), and display the results.](image-9.png "Step7. 3 :") 23![Estimation of External Power System Dynamic Equivalent Parameters Global Journal of Researches in Engineering Volume ?? = ??? ?? 1 ?? , ?? 2 ?? ? ?? 1 ?? ??? 2 ?? 2](image-10.png "Nonlinear 2 ?? 3") 5![Fig.5 : Original system IEEE-39 bus](image-11.png "Fig. 5 :") 6![Fig. 6 : Equivalents system](image-12.png "Fig. 6 :") ![Table III : Equivalent Generator Parameters Dynamic responses are shown in fig. 7 to fig. 9. The boundary node voltage is shown by fig. 7, The angle of the machine n°10 is shown in fig. 8 and angular velocity of the machine n°8 in fig. 9.](image-13.png "") . I F © 2012 Global Journals Inc. (US) 2012 July ð??" = F © 2012 Global Journals Inc. (US) 2012 July © 2012 Global Journals Inc. (US) 2012 July F © 2012 Global Journals Inc. (US) 2012 July Fig. 8 : Angle of the machine 10 * Dynamic reduction of large power systems for stability studies WangLei MKlein SYirga PKundurl IEEE Trans. Power Systems 12 2 1997 * Dynamic equivalent modeling of large power systems using structure preservation technique MLOurari LADessaint DVan-Que IEEE Trans on Power Systems 21 3 2006 * Coherency and aggregation techniques incorporating rotor and voltage dynamics Sung-Kwan LiuJoo LEChen-Ching Jong-WoongJones Choe IEEE Trans. Power Systems 19 2 2004 * Testing of the modal dynamic equivalents technique WWPrice EMGulachensk IEEE Trans on Power Apparatus and Systems 94 4 1975 IEEE Trans. measurements * A technique to reduce power systems electromechanical models JMRamfirez RGValle IEEE Transa on Energy Conversion 19 2 2004 * Application of Dynamic Equivalents Based on Identificaiton of Coherent Generator Group in Engineering YangJingping XuZheng Power System Techonology 29 17 2005 * Dynamic Equivalent Using Operating Data and Stochastic Modeling Ibrahim M A H Mostafa O M -Abiad A HEl IEEE PAS 95 5 1976 * Real Time External Equivalents for Static Security Analysis Alvarado F L IEEE TPAS 98 1979 * Power System Dynamics YuY 1983 Academic Press * Synchronous Machine Modeling Mechanization and System Performance Study JEFagan 1977 The University of Texas at Arlington Thèse de doctorat * Obtaining dynamic equivalents through the minimization of a line flows function JMRamirez Arredondo International Journal of Electrical Power & Energy Systems 21 1999 * Modal reduced dynamic equivalent model for analog type power system simulator KNojiri SSuzaki KTakenaka MGoto IEEE Trans. Power Systems 12 4 1997 * Estimation of external dynamic equivalents of a thirteen system MAYu Yaonan El-Sharkawi IEEE Trans. Power Apparatus and Systems 100 3 1981 * WWPrice DNEwart EMGulachenski RFSilva Dynamic equivalents from on-line