Probabilistic Assessment of Power System Performance Quality

In a recently published work by the authors, a novel framework was developed and applied for assessment of reliability and quality performance levels in real-life power systems with practical large-scale sizes. The new assessment methodology is based on three metaphors (dimensions) representing the relationship between available generation capacities and required demand levels. The developed reliability and performance quality indices were deterministic in nature. That is, they represent one operating state (a snapshot of the system conditions) in which the required demand as well as the generation and transmission capacities are known with 100% certainty. In real life, however, load variations occur randomly so as the contingencies which cause some generation and/or transmission capacities to be lost (become unavailable). In other words, neither the load levels nor the generation or transmission capacities are known with absolute certainty. They are rather subject to random variations and, consequently, the calculated reliability and performance quality indices are all subject to random variations where only expected values of these indices can be evaluated. This paper presents a major extension to the previously published work by developing a theory and formulas for computing the expected values of different system reliability and performance quality indices. In this context, a “contingency scenario” or a system “demand level” are regarded, in a more general sense, as a “state”, which occurs with certain probability and represents a given demand value and availability pattern of various capacities in the system. The work of this paper provides a practical and meaningful methodology for real-life assessment of power system reliability and performance quality levels. Practical applications are also presented, for demonstration purposes, to the Saudi electricity power grid.


Introduction
Maintaining a continuous and sufficient power supply to the customers at a reasonable cost is the prime objective of electric power companies around the world.In this regard, power system cost-effectiveness, security, adequacy and reliability analyses have become a major concern in today's highly-competitive business environment of power utility planning and operations [1][2][3].In a recent paper by the authors [4], a novel framework was developed and applied for assessment of reliability and quality performance levels in real-life power systems with practical large-scale sizes.The new assessment methodology is based on three metaphors (dimensions) representing the relationship between available generation capacities and required demand levels.The first metaphor defines whether or not the capacity exists, the second metaphor defines whether or not the capacity is needed, and the last metaphor defines whether or not the capacity can reach (delivered to) the demand.The eight possible combinations associated with the 0/1 (Yes/No) values of the three metaphors would, in turn, define a set of powerful systemwide performance quality measures relating to generation deficiency, redundancy, bottling, etc.The developed reliability and performance quality indices were deterministic in nature.That is, they represent one operating state (a snapshot of the system conditions) in which the required demand as well as the generation and transmission capacities are known with 100% certainty.
In real life, however, load variations occur randomly so as the contingencies which cause some generation and/or transmission capacities to be lost (become unavailable).In other words, neither the load levels nor the generation or transmission capacities are known with absolute certainty.They are rather subject to random variations and, consequently, the calculated reliability and performance quality indices are all subject to random variations where only expected values of these indices can be evaluated.
Methods for computing probabilistic contingency-based reliability and performance quality indices have previously been published in the literature [5][6][7].These methods are based on a combined contingency analysis and reliability evaluation scheme which integrates both the contingency effect and its probability of occurrence into one routine of analysis.In the present research work, similar analysis will be used to compute the expected values of different system reliability and performance quality indices.In this context, a "contingency scenario" or a system "demand level" are regarded, in a more general sense, as a "state", which occurs with certain probability and represents a given demand value and availability pattern of various capacities in the system.
The work of this paper presents a major extension to the previously published work [4] by developing a theory and formulas for computing the expected values of different system reliability and performance quality indices.In this context, a "contingency scenario" or a system "demand level" are regarded, in a more general sense, as a "state", which occurs with certain probability and represents a given demand value and availability pattern of various capacities in the system.
The work of this paper provides a practical and meaningful methodology for real-life assessment of power system reliability and performance quality levels.Practical applications are also resented in the Saudi electricity power grid.

Performance Quality Framework
In the framework presented in [4], three metaphors (dimensions) were introduced to represent the relationship between certain system generation capacity and the demand.These metaphors relate to the following demand fulfillment issues: 1) Need of capacity for demand fulfillment; 2) Existence of capacity (availability for demand fulfillment); 3) Ability of capacity to reach the demand.The first metaphor defines whether or not the capacity is needed, the second metaphor defines whether or not the capacity exists, and the last metaphor defines whether or not the capacity can reach (delivered to) the demand.The eight possible combinations associated with the 0/1 (Yes/No) values of the three metaphors would, in turn, define a set of powerful system-wide performance quality measures, namely: 1) Utilized: A given capacity is said to be utilized if it is needed (for demand fulfillment), exists, and can reach the demand; 2) Bottled: A given capacity is said to be bottled if it is needed (for demand fulfillment) and exists, but cannot reach the demand; 3) Shortfall: A given capacity is said to be shortfall if it is needed (for demand fulfillment) and, anyhow, does not exist and can reach the demand; 4) Deficit: A given capacity is said to be deficit if it is needed (for demand fulfillment) but, however, does not exist and cannot reach the demand; 5) Surplus: A given capacity is said to be surplus if it is not needed (for demand fulfillment) although exists and can reach the demand; 6) Redundant: A given capacity is said to be redundant if it is not needed (for Demand fulfillment) although exists but, anyhow, cannot reach the demand; 7) Spared: A given capacity is said to be spared if it is not needed (for demand fulfillment) and, anyhow, does not exist although can reach the demand; 8) Saved: A given capacity is said to be saved if it is no needed (for demand fulfillment) and, anyhow, does not exist and cannot reach the demand.
We note here that the above performance quality measures are associated with different combinations (topples) of the three quality metaphors, namely, "existence", "need" and "ability to reach the demand".The corresponding quality state of a given capacity can be represented by a threevalue expression of either a "Yes/No" or "1/0" type indicating the true/false value associated with each quality metaphor.
The evaluation of the above quality indices requires the knowledge of the following data types for the demand and various system facilities: 1) The value of demand required to be supplied; 2) The value of generation capacity as well as the maximum site capacity (the limit of potential increase in existing generation capacity); 3) The value of transmission capacity.

Linear Program Formulation
In the computational scheme of [4], the integrated system quality assessment is performed via solving a master linear programming problem [8] in which a feasible power flow is established which minimizes the total system non-served load subject to capacity limits and flow equations.The master linear program, which utilizes the network bus incidence matrix A, is formulated as   In the above master linear program, T P = vector of n T elements representing transmission branch capacities; L P = vector of n L elements of peak bus loads; G P = vector of n G elements representing generator capacities.
Also, in the above master linear program (1), P L , P G , and P T are n L , n G and n T column vectors representing the actual load bus powers (measured outward), generator bus powers (measured inwards) and transmission line powers (measured as per the network bus incidence matrix A), respectively.The solution of the above linear program provides a more realistic (less conservative) flow pattern in view of the fact that when load curtailments are anticipated, all system generation resources would be re-dispatched in such a way which minimizes such load cuts.The feasible flow pattern established from the Master Linear Program is then used to evaluate various integrated system quality indices through a set of closely related sub-problems.

Implementation Mechanisms
For real life power systems with practical sizes, the quality indices cannot be evaluated by inspection.An appropriate computerized scheme is needed in order to properly evaluate various quality indices according to their stated definitions.The master linear program presented before forms the bases for analyzing and evaluating the quality indices.For example, the Load Supply Reliability can be evaluated as follows: where the bus loads at the solution of the master linear program are termed as , and P l denotes the solution load value at bus (l).

 
1 l P On the other hand, generation quality indices are defined in terms of the previously defined "1/0" states indicating the (Needed, Exists, Can-reach) true/false values associated with each quality metaphor.We shall use the symbol Q g ijk to indicate the generation quality index state.Also, in the following expressions, we shall use Min{x, y, •••, z} to indicate the minimum of x, y, •••, z.The notation <x> will be used to denote Max{0, x}, that is the maximum of x and zero (=x if x > 0, or 0 otherwise).For example, the Utilized Generation Capacity index is given by 111 Utilized Capacity needed, exists, can reach Similarly, the Bottled Generation Capacity index is given by where the generation output values P g are calculated at the solution of the linear program with open limits on the loads.

Probabilistic Reliability Indices
The power system can be described, for the purpose of composite reliability and performance quality assessment, by the three-component model, as shown in Figure 1, in which generation, transmission and load are considered as multi-state elements of the power system.For a given operating state m, the values of the network variables will be the solution of the maximum load-supply optimization problem described in the previous section.Also, let f m be the probability of operating state m (the sum of f m for all m, including base-case scenario is 1).Then, the following three system-wide reliability indices may be defined: 1) Loss of Load Probability where represents the system loss of load probability for any operating state m (load level, loss of generation and/or transmission capacities) in the power grid, Copyright © 2012 SciRes.EPE represents the loss of load probability at bus  for operating state m, variations of system demand level as well as forced out-ted fo hand, the randomness in the generation an

Applications to SEC Power System
where n L is the number of load buses in the system, represents the expected value of load not-served at bus , represents the expected value of Load Not-Served at bus  for the operating state m, and Load not served at buslfor operating state m LNS m   which is obtained from the solution of the Master Linear Program (1).
3) Expected Energy Not Served represents the expected value of energy not served at a bus , and represents the expected value of energy not served at bus  for operating state m, represents the energy not served at bus  for operating state m, and T (m) denotes the time duration of operating state m.

Probabilistic Performance Quality Indices
Probabilistic performance quality indices can be calculated using the previously derived formulas based on the solution of the master linear program (1) subject to random ages in various generation and transmission facilities.
For example, the load variations, which are accoun r using the so-called "load-duration curves" can be used to calculate the expected value of the Load Not-Served (LNS), which is widely known as the Expected Load Not-Served (eLNS).
On the other d transmission capacity availability are accounted for using the so-called forced-outage rates (or availability rates) associated with various facilities.Consequently, the expected values of the performance quality indices Q g 111, Q g 110, Q g 101, etc., denoted by eQ g 111, eQ g 110, eQ g 101, etc., can be evaluated using the modeled randomness of the system load as well as the generation and transmission capacity availabilities.

SEC Quality Performance Indices
In a recently completed industry supported stu tions were conducted on a practical power system comprising a portion of the interconnected Saudi power grid.The power system consists of two main regions, namely the Central region and the Eastern region.The two systems are interconnected through two 380 kV and one 230 kV double-circuit lines.Four zones are identified in the present analysis, three in the Central region (Riyadh, Qassim and Hail zones) and one in the Eastern region.In this application, three reliability and quality performance indices are considered, namely the system Load Not-Served (LNS), Utilized Generation Capacity (Q g 111) and the Bottled Generation Capacity (Q g 110).In the present work, a particular focus will be on Qassim and Hail zones for demonstration purposes.The system models used for these two zones are shown in Figures 2(a) and (b), respectively.Table 1 outlines the network data in terms of generation and transmission facilities as well as system loads.The results obtain ns.For example, the results obtained for the isolated network scenario of Qassim zone (Figure 3(a)) show that the Load Not-Served is non-zero even for relatively low    at a demand level of 1840 MW to reach 2400 MW when the demand level is 4410 MW.This problem is clearly mitigated in the interconnected network scenario of Qassim zone (Figure 3(b)), where generation support from Riyadh zone becomes available.In this case, the Load Not-Served stays at zero value for all demand levels up to 2620 MW where it starts to increase slowly to reach 70 MW at a demand level of 3370 MW before it starts to increase sharply afterwards to reach about 2000 MW at a demand level of 5610 MW.
The Utilized Gen olated network scenario of Qassim zone (Figure 3(a)) increases continuously with the required demand level until it saturates at about 2000 MW when the required demand reaches 2943 MW when no more available generation can be utilized.This situation is avoided-as expectedin the interconnected network scenario of Qassim zone (Figure 3(b)) where the Utilized Generation Capacity increases continuously to reach, for example, 3600 MW at demand level of 5610 MW as more generation support becomes available.The Bottled Generation Capacity index Q g 110, for the isolated network scenario of Qassim zone (Figure 3(a)), decreases continuously with the required demand until it disappears at a demand level of 2575 MW.In the case of the interconnected network scenario of Qassim zone (Figure 3(b)), however, the Bottled Generation Capacity coincides with the Load Not-Severed for all required demand levels up to 4115 MW.After this level, the Bottled Generation Capacity starts to decrease continuously.
As was stated before, the Hail network under in underground cables and power transformers) as well as 46 loads.In order to simplify the probabilistic assessment, the total available generation is represented by three equivalent units, two of which represent the internal available generation within Hail and the third represents the interconnection support from Qassim.The availability rate of the generating unit is 0.9825.On the other hand, only the five major transmission lines in Hail are considered in the probabilistic assessment with equal availability rate of 0.985.The other transmission lines are assumed to be available all the time.
Using the results of the probabilistic analysis of Hail system, the di s reliability and performance quality indices can be evaluated and displayed.These discrete density functions show the overall probabilities of occurrence associated with certain values of the system performance indices.The probability density function of the Load Not Served (LNS) for Hail network at required load level 150 MW is depicted in Figure 6 (from Figure 6(a)) while there is a probability of 0.0310 MW.On the other hand, there is a probability of 0.92 that the Surplus Generation Capacity (Q g 011) in Hail (at load level 200 MW) is equal to, or greater than 1175 MW.
When all required load levels are considered with their respective probabilities of occurrence (as per the load du n curve), an overall estimate of the expected value of reliability and performance quality indices can be obtained as shown in   1 that the Hail generation capacity is 594 MW (base value) while the required base load is 655 MW, one may conclude that the 64% expected value for the utilized generation is a reflection of the relatively large variations of the required load during different periods of the year, as is evident from the load duration curve of Figure 5.

Conclu
The work of this p the previously publi formulas for computing the expected values of different system reliability and performance quality indices.The reliability and performance quality indices, when evaluated at a given load level and a certain scenario of available generation and transmission capacities, would provide indication on system performance for only such a particular system condition (snapshot).However, the novel formulation presented in this paper can accommodate the randomness associated with the load level as well as the availability of generation and transmission capacities.In this case, expected values of reliability indices, such as the Expected Load Not-Served (eLNS), as well as the expected values of performance quality indices, such as Utilized Generation Capacity (eQ g 111), Bottled Generation Capacity (eQ g 110), Shortfall Generation Capacity (eQ g 101), Deficit Generation Capacity (eQ g 100), Surplus Generation Capacity (eQ g 011), Redundant Generation Capacity (eQ g 010), Spared Generation Capacity (eQ g 001) and Saved Generation Capacity (eQ g 000) can be calculated using the load duration curve as well as the availability rates of generation and transmission facilities in the system.
The results of the probabilistic analysis of Hail system have NS) in Hail being zero (when the required load level is 150 MW) is 0.9387, while there is a probability of 0.0310 that the unsupplied load is equal to, or greater than 25 MW.owledgements the Saudi Electricity Com-
dy, applica-(a) and (b) depict the variation of nd, show 3-dimensi ed reveal several important observatio (2d) and denotes the scheduled (required) load at load bus .Also, in Equation (3.8), M s denotes the number of all possible states. o

Fig- ures 3
(a) and (b), on the other hand, summarize the results of the performance quality measures applied to the SEC power system for various system status (isolated or connected) of each zone.In particular, Figures 3 the variation of quality indices (LNS, Q g 111 and Q g 110) with the required load level of the Qassim isolated and interconnected network, respectively.Figures4(a) and (b), on the other ha onal graphs depicting the variation of Utilized Generation Capacity index Q g 111 with both load and generation capacity levels of the Hail isolated and interconnected network, respectively.

Figure 3 .
Figure 3. (a) Variation of qu of the Qassim isolated network; ality indices (LNS, Q g 111 and Q g 110) with the required load level (b) Variation of quality Indices (LNS, Q g 111 and Q g 110) with the required load level of the Qassim interconnected network.

Table 2 .
The overall expected value of the Load Not Served (LNS) for Hail is 4.01 MW, w ad in Hail.On the other hand, overall expected value of the Bottled Generation Capacity (Q g 110) is 3.86 (MW),

Table 2 . Expedited values of reliability and performance quality indices for Hail network.
Utilized Generation Capacity (Q g 111) is 378.1 (MW), which represents almost 64% of the total generation capacity in Hail.In other words, about 36% of Hail available generation capacity is expected to be unutilized.Recalling from Table