Keywords
|
deregulation, contingency, appraisal, line outage, generator outage. |
INTRODUCTION
|
The power system consists of generation, transmission, distribution bundled together. In deregulation unbundling of power system network for efficiency, reliability and least price of power to the customer [1]. Under deregulation minimum price of power transfer to the utility, that time more number of buyers to buy the power from the generation. But all the transmission lines and generators have some limits [2]. Whenever the demand of power is maximum than compared with the transfer limits, the line will be damaged. The power demand is reduced then the generator is reliving from the power system network [3-6]. The system security will be collapsed. So the secured dispatch scheduling of power market is important due to open access and competition. Before security assessment the contingency analysis is needed [7]. The contingency limits are based on system operator experience. But this methodology not predicts the security limits. Before secured dispatch scheduling we have to analysis the contingency under single outage of line, generator and multiple outage of combination of both line and generator [8-10]. Whenever the maximum violation in the element some line and generator get damaged. So the contingency analysis is essential. Various methodologiesare used for contingency analysis such as fast decoupled method, DC power flow, neural network and fuzzy theory [10]. But full NewtonâÃâ¬ÃŸs method is high accurate, reliable and better solution. |
In this paper presents the contingency analysis of test system by using full NewtonâÃâ¬ÃŸs method in power world simulator and two important factor of line outage distribution factor and generation shift factor. In state estimation based contingency analysis more complex. So the student friendly software of power world simulator is used for analysis. The proposed algorithm of NewtonâÃâ¬ÃŸs method based contingency selection is presented. The mathematical formulation is in section II. The contingency analysis algorithm is in power world simulator in section III. The description of test system is in section IV. In section V includes simulation results with description and conclusion from the result in section VI. |
PROBLEM FORMULATION
|
The power world simulator can be set to use a full Newton solution or use a DC load flow method to analyse each contingency. The full Newton approach is not as fast as a DC load flow, but the results tend to be significantly more accurate and allow for gauging voltage/VAR effects. The Newton solution method (also called Newton-Rapson method) is more efficient for large power systems. The number of iteration required to obtain a solution is independent of a system size but more functional evaluation are required at each iteration. |
Equation for bus Admittance matrix |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e001.gif) |
In above equation j includes bus i expressing this equation in polar form, we have |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e002.gif) |
The complex power at bus |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e003.gif) |
Substituting from (2) for Ii in (3) |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e004.gif) |
Separating the real and imaginary parts |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e004.gif) |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e006.gif) |
Equation (5) and (6) constitute of nonlinear algebraic equation in terms of the independent variables, voltage magnitude in per unit and phase angle in radians. |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e007.gif) |
In above equation, bus 1 is assumed to be slack bus. The jacobian matrix gives the linearized relationship between small changes in voltage angle Δδi(k) and voltage magnitude Δ|Vi(k)| with small changes in real and reactive power ΔPi(k) and ΔQi(k) elements of jacobian matrix are the partial derivatives of (5) and (6) evaluated at Δδi(k) and Δ|Vi(k)| |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e008.gif) |
Accordingly there are (n-1) real power constraints and (n-1-m) reactive power constraints and the jacobian matrix is the order of (2n-2-m) (2n-2-m). |
J1 is the order of (n-1) x (n-1) |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e009.gif) |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e010.gif) |
J2 is the order of (n-1) x (n-1-m) |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e011.gif) |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e012.gif) |
J3 is the order of (n-1-m) x (n-1) |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e013.gif) |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e014.gif) |
J4 is the order of (n-1-m) x (n-1-m) |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e015.gif) |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e016.gif) |
The terms ΔPi(k) and ΔQi(k) are difference between the schedule and calculated values, known as the power residuals, given by |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e017.gif) |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e018.gif) |
The new estimated for bus voltage is |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e019.gif) |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e020.gif) |
Another solution method in simulator is DC load flow. These factors can be derived in a variety of ways and basically come down to two types. |
a) Generation shift factors |
b) Line outage distribution factor |
The generation shift factors are designated Ali and have the following definition |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e021.gif) |
The line outage distribution factors are used in a similar manner, only they apply to the testing for overloads when transmission circuits are lost. By defining the line outage distribution factor has the following meanings. |
![image](https://www.rroij.com/articles-images/IJAREEIE-675-e022.gif) |
PROPOSED APPROACH
|
Before the secured dispatch scheduling the contingency analysis is the important one for selection of contingency element in the maximum violation of the system network. The following steps to involve in the contingency analysis. |
A. CONTINGENCY ANALYSIS ALGORITHM
|
Step 1: Draw the Simulink one line diagram in new case window of power world simulator for the given power System in edit mode. |
Step 2:Save the case with apt name. |
Step 3:Select run mode. |
Step 4:Play or Run the one line diagram in tool menu. |
Step 5:Select CONTINGENCY ANALYSIS in tool menu, then the contingency analysis dialogue box is open. |
Step 6:Right click on label and select auto insert contingencies through insert special option. |
Step 7:Verify that single transmission line or transformer is selected. |
Step 8:If can limit the contingencies inserted to only those meeting a definedfilter. |
Step 9:We want to insert contingencies for all branches and generators so nofiltering is desired. |
Step 10: To check the following conditions |
a) Remove the checkmark in use area/ zone filters. |
b) Verify no other options are selected. |
Step 11: Click do insert contingencies button to accept the all contingencies. |
Step 12: Click YES to get the contingencies. |
Step 13: Now the contingency analysis dialog shows contingencies. |
a) Right click on the list display on the contingency tap and select insert special and click auto insert to the local menu |
b) Select single generating unit then click the do insert contingencies button. Click YES to complete. |
Step 14:The auto insert tool did not insert a contingency for the generator connected to the slack bus. |
Step 15:Click „start runâÃâ¬ÃŸ on the contingencies tab click start on summary tab or Run contingency. |
Step 16:Select the maximum violation of contingency analysis taken for account in the secured dispatch of deregulation of power market. |
TEST SYSTEM
|
The contingency analysis of 14-bus test system is shown below when the power flow is running on the power world simulator. The percentage of power flow is mentioned in power flow diagram. It consists of five generators for dispatch of power. |
SIMULATION RESULTS
|
The contingency analysis is classified as single contingency, multiple contingency. In IEEE test system having 110 contingency. But in maximum violation of single contingency of line and generator is taken in account in table1.which depends on maximum branch percentage and minimum voltage the contingency element is taken in account for security dispatch of power. |
The combination of both generator and line contingency is also in security analysis. The maximum violation only taken in account in table 2. The L_ indicates line and G_ indicates generator. |
The generation shift factor describe a generator power sensitivity under the contingency condition is shown in table 3. It consists of 5 generators and their sensitivity. |
Line outage distribution factor (LODF) is used to approximate the change in the flow on one line caused by the outage of a second line. It is a real time operation. The percentage of LODF is tabulated in each line in table 4. |
CONCLUSIONS
|
The Full NewtonâÃâ¬ÃŸs method based contingency analysis algorithm is high accuracy and better efficiency. The contingency analysis in power world simulator is easy to run the power system and more reliable than compared with state estimation based contingency analysis. The security limits described from maximum violation of the element of test system and sensitivity analysis of both line outage distribution factor and generation shift factor. In feature power world simulator based contingency analysis is widely used for secured dispatch scheduling and the demand response improvement of deregulated power market. Voltage profile is varied which depends on violated and unviolated line is shown in Figure 2 and 3. |
Tables at a glance
|
![Click Here Table icon](https://www.rroij.com/articles-images/table-icon.gif) |
![Click Here Table icon](https://www.rroij.com/articles-images/table-icon.gif) |
![Click Here Table icon](https://www.rroij.com/articles-images/table-icon.gif) |
![Click Here Table icon](https://www.rroij.com/articles-images/table-icon.gif) |
Table 1 |
Table 2 |
Table 3 |
Table 4 |
|
|
Figures at a glance
|
![Click Here Figure 1](https://www.rroij.com/articles-images/IJAREEIE-675-g001.gif) |
![Click Here Figure 2](https://www.rroij.com/articles-images/IJAREEIE-675-g002.gif) |
Figure 1 |
Figure 2 |
|
|
References
|
- Amit Kumar Roy “contingency analysis in power system” Thapuar University Patiala.2009.
- Stott B and Alsac O, “Fast decoupled load flow”, IEEE Transactions on Power Apparatus and Systems, Vol. PAS-91, No. 5, pp 859- 869, May 1974.
- Peterson N.M, Tinney W.F and Bree D.W, “Iterative linear AC power flow solution for fast approximate outage studies”, IEEE Transactions on Power apparatus and system, vol.PAS-91, No.5pp 2048-2058, October 1972.
- Richard D.Christie, AnjanBose “load frequency control issues in power sytem operations after deregulation” University of Washington.1995.
- Saavedra, O.R., "Solving the security constrained optimal power flow problem in a distributed computing environment," Generation, transmission and distribution,IEEE proceedings,vol.143,No.6 pp 593-598,1996.
- Singh S.N and Srivastava S.C, “Improved voltage and reactive distribution factor for outage studies”, IEEE Transactions on Power systems, Vol. 12, No.3, pp 1085-1093 Auguest 1997.
- Rajasekaran S. and Vijayalakshmi G.A., “Neural Networks, Fuzzy logic and GeneticAlgorithm Synthesis and applications”, learning Private limited 2010.
- Salah EldeenGasim Mohamed1, AbdelazizYousif Mohamed1 and Yousif Hassan Abdelrahim “Power System Contingency Analysis to detect Network weakness” ZaytoonehUniversity International Engineering Conference on design and Innovation in Infrastructure 2012.
- SivaramaKarthikVijapurapu “contingency analysis of power system in presence of geomagneticallyinduced currents” university of Kentucky 2013
- Sidhu T.S and Cui L., “Contingency Screening for Steady State Security Analysis Byusing FFT and Artificial Neural Networks”,IEEE Transactions on Power system,Vol. 15,No. 1, pp.421-426, February 2000.
|