Predictive Reliability Assessment of Generation System

Due to increasing load and characteristic stagnation and fluctuations of existing generation systems capacity, the reliability assessment of generation systems is crucial to system adequacy. Furthermore, a rapid load increase could amount to a consequent sudden deficit in the generation supply before the next scheduled assessment. Hence, a reliability assessment is conducted at regular and close intervals to ensure adequacy. This study simulates and establishes the relationship between the load growth and generation capacity using the generation and load data of the IEEE reliability test system (IEEE RTS ‘96 standard). The generation capacity states and the risk model were obtained using the sequential Monte Carlo simulation (MCS) method. The load was gradually increased stepwise and is simulated against the constant generation capacity. In each case, the reliability index was recorded in terms of loss-of-load evaluation (LOLE). The recorded reliability index was thereafter fitted with the load-growth trend by the linear regression approach. A predictive assessment approach is thereafter proffered through the obtained fitting equation. In addition, a reliability threshold is effectively determined at a yield point for a reliability benchmark.


Introduction
The economic growth of nations is positively impacted by their electrical power supply adequacy. The menace of inadequacy could be disastrous and catastrophic both to the consumers' side and the utility company. This ranges from a huge loss of revenue to death in critical cases such as in the health and aviation sector, yet the insufficiency of generated power is still a natural factor during load growth. In addition, the tendency of generation plants to become derated during operation intervals is a natural factor. For instance, a hybrid generation system comprising of solar, wind, and hydro plants would have its solar capacity derated at night due to photon absence. In addition, the hydro and wind source would be most likely derated in the summer. Because the generation capacity is ideally either constant or derated while the load uniquely rises, electrical loads would experience deficit supply in the future. The reliability assessment of the generation system ensures that the available generation capacity adequately covers the consumers' load at all times [1]. It checks the susceptibility of the load growth to the generation capacity supply shortage. Thus, it is performed at frequent time intervals to assist power providers and planners to prevent loss of load [2,3]. This offers a timely generation The generation system consists of one or more generation source(s) known as the generation unit(s). Generation units could be "in service (up)" or "out of service (down)" states when needed due to planned or unplanned outages [17]. The generation system reliability is usually calculated and quantified using the reliability index which reveals the system reliability strength [18].

The Generation System Reliability Index
The impact of the continued load growth on the reliability of the generation system was examined using the loss-of-load evaluation (LOLE) reliability assessment metric. This is a widely accepted and most widely utilized probabilistic method in power generation reliability evaluation [19]. The LOLE presents the total time duration during which the load suffers power supply insufficiency within a period usually a year [11]. From Figure 2, the LOLE was elementarily calculated as shown in (1) and conventionally measured in hour-per-year (h/yr). 1 2 where 1 t means generation failure duration 1 and 2 t means generation failure duration 2.  The generation system consists of one or more generation source(s) known as the generation unit(s). Generation units could be "in service (up)" or "out of service (down)" states when needed due to planned or unplanned outages [17]. The generation system reliability is usually calculated and quantified using the reliability index which reveals the system reliability strength [18].

The Generation System Reliability Index
The impact of the continued load growth on the reliability of the generation system was examined using the loss-of-load evaluation (LOLE) reliability assessment metric. This is a widely accepted and most widely utilized probabilistic method in power generation reliability evaluation [19]. The LOLE presents the total time duration during which the load suffers power supply insufficiency within a period usually a year [11]. From Figure 2, the LOLE was elementarily calculated as shown in (1) and conventionally measured in hour-per-year (h/yr).
where t 1 means generation failure duration 1 and t 2 means generation failure duration 2.
Energies 2020, 13, x FOR PEER REVIEW 3 of 13 Figure 1. Power system reliability assessment path.
The generation system consists of one or more generation source(s) known as the generation unit(s). Generation units could be "in service (up)" or "out of service (down)" states when needed due to planned or unplanned outages [17]. The generation system reliability is usually calculated and quantified using the reliability index which reveals the system reliability strength [18].

The Generation System Reliability Index
The impact of the continued load growth on the reliability of the generation system was examined using the loss-of-load evaluation (LOLE) reliability assessment metric. This is a widely accepted and most widely utilized probabilistic method in power generation reliability evaluation [19]. The LOLE presents the total time duration during which the load suffers power supply insufficiency within a period usually a year [11]. From Figure 2, the LOLE was elementarily calculated as shown in (1) and conventionally measured in hour-per-year (h/yr).
where 1 t means generation failure duration 1 and 2 t means generation failure duration 2.

The Generation Model, Load Model and the Risk Model
To evaluate the reliability of generation systems, each generator is usually represented with a model. This reflects the performance of the generating unit sequentially in the mutually exclusive Energies 2020, 13, 4350 4 of 13 states, namely, up (operating) or down states [11]. This is as shown in Figure 3. A single generator's model is referred to as a 'unit model'. Individual unit models are combined by a superimposition method to obtain a resultant single generation model (resultant system capacity) as seen in Figure 4.

The Generation Model, Load Model and the Risk Model
To evaluate the reliability of generation systems, each generator is usually represented with a model. This reflects the performance of the generating unit sequentially in the mutually exclusive states, namely, up (operating) or down states [11]. This is as shown in Figure 3. A single generator's model is referred to as a 'unit model'. Individual unit models are combined by a superimposition method to obtain a resultant single generation model (resultant system capacity) as seen in Figure 4.  To model the generation capacity, the operation performance of the generating unit is characteristically expressed using the time-to-fail (TTF) and time-to-repair (TTR) parameters. TTF is the duration that a generating unit in the up-state lasts before entering the failed state. TTR is the duration that the generating unit remained in the down-state (failed state) before getting back to the up-state (operating state). TTF and TTR are illustrated in Figure 5 and are obtained as described in (2) and (3), respectively via the sequential Monte Carlo simulation (MCS) method [11].
where λ and µ are the failure rate and repair rate, respectively and U1 and U2 are sets of generated random numbers, 'x' ranging from zero (0) to one (1), where x = x1, x2, x3, ….., xn and 0 < x < 1. The random numbers (U1 and U2) are generated using a random number generator. The value of the failure rate λ is given by: λ = 1/MTTF (4) where MTTF is the mean time to fail. The value of the repair rate µ is described by:

The Generation Model, Load Model and the Risk Model
To evaluate the reliability of generation systems, each generator is usually represented with a model. This reflects the performance of the generating unit sequentially in the mutually exclusive states, namely, up (operating) or down states [11]. This is as shown in Figure 3. A single generator's model is referred to as a 'unit model'. Individual unit models are combined by a superimposition method to obtain a resultant single generation model (resultant system capacity) as seen in Figure 4.  To model the generation capacity, the operation performance of the generating unit is characteristically expressed using the time-to-fail (TTF) and time-to-repair (TTR) parameters. TTF is the duration that a generating unit in the up-state lasts before entering the failed state. TTR is the duration that the generating unit remained in the down-state (failed state) before getting back to the up-state (operating state). TTF and TTR are illustrated in Figure 5 and are obtained as described in (2) and (3), respectively via the sequential Monte Carlo simulation (MCS) method [11].
where λ and µ are the failure rate and repair rate, respectively and U1 and U2 are sets of generated random numbers, 'x' ranging from zero (0) to one (1), where x = x1, x2, x3, ….., xn and 0 < x < 1. The random numbers (U1 and U2) are generated using a random number generator. The value of the failure rate λ is given by: λ = 1/MTTF (4) where MTTF is the mean time to fail. The value of the repair rate µ is described by: To model the generation capacity, the operation performance of the generating unit is characteristically expressed using the time-to-fail (TTF) and time-to-repair (TTR) parameters. TTF is the duration that a generating unit in the up-state lasts before entering the failed state. TTR is the duration that the generating unit remained in the down-state (failed state) before getting back to the up-state (operating state). TTF and TTR are illustrated in Figure 5 and are obtained as described in (2) and (3), respectively via the sequential Monte Carlo simulation (MCS) method [11].
where λ and µ are the failure rate and repair rate, respectively and U 1 and U 2 are sets of generated random numbers, 'x' ranging from zero (0) to one (1), where x = x 1 , x 2 , x 3 , . . . .., x n and 0 < x < 1. The random numbers (U 1 and U 2 ) are generated using a random number generator. The value of the failure rate λ is given by: where MTTF is the mean time to fail. The value of the repair rate µ is described by: where MTTR is the mean time to repair. The Equations (2) and (3) are suitable for describing state intervals of stochastic systems with any number of discrete states [11]. In our paper, they were Energies 2020, 13, 4350 5 of 13 calculated using the generation system data from the IEEE RTS '96 standard as shown in Appendix A in Table A1. It consists of nine (9) generation sources, each consisting of several generating units. The values of λ and µ are also obtainable from the IEEE RTS standard for each of the generating units via (4) and (5). The states' durations of each generating unit, i.e., TTF and TTR, are sequentially calculated as follows: (1) For each of the random numbers generated in U 1 , the TTF is calculated using (2); (2) For each of the random numbers generated in U 2 , the TTR is calculated using (3); (3) The resultant series of TTF and TTR are thereafter graphically and sequentially joined (one after the other) (i.e., TTF, then TTR, then TTF, and so on) as shown in Figure 5. This is continued until the total value reaches 8736 representing the total number of hours of a year. This is performed for each of the entire generating units. In the end, they are combined by a superimposition method to obtain a single resultant system capacity as shown in Figure 4. This is the annual duration of the generation capacity states, otherwise known as the generation model. where MTTR is the mean time to repair. The Equations (2) and (3) are suitable for describing state intervals of stochastic systems with any number of discrete states [11]. In our paper, they were calculated using the generation system data from the IEEE RTS '96 standard as shown in Appendix A in Table A1. It consists of nine (9) generation sources, each consisting of several generating units. The values of λ and µ are also obtainable from the IEEE RTS standard for each of the generating units via (4) and (5). The states' durations of each generating unit, i.e., TTF and TTR, are sequentially calculated as follows: (1) For each of the random numbers generated in U1, the TTF is calculated using (2); (2) For each of the random numbers generated in U2, the TTR is calculated using (3); (3) The resultant series of TTF and TTR are thereafter graphically and sequentially joined (one after the other) (i.e., TTF, then TTR, then TTF, and so on) as shown in Figure 5. This is continued until the total value reaches 8736 representing the total number of hours of a year. This is performed for each of the entire generating units. In the end, they are combined by a superimposition method to obtain a single resultant system capacity as shown in Figure 4. This is the annual duration of the generation capacity states, otherwise known as the generation model. The load model is a graphic representation of the consumers' hourly chronological peak load for a period of one year. It, thus, consists of 8736 values which represent 8736 h of a year [11]. This is shown in Figure 6 through the red line in the graph of the superimposed generation and load model. The risk model represents the superimposition of the generation and the load model from which the LOLE reliability index is calculated [20]. The LOLE is the sum of the individual algebraic differences between the load level and the generation capacity at the points where the load level exceeds the generation capacity. This is as shown by the two red shaded regions in Figure 6. This is usually calculated in a yearly period. The generation system has a record of the total capacity for each hour, hence 8736 values corresponding to 8736 h of the year. In addition, the load system has a record of the consumers' hourly load usage level, hence, also 8736 values. Since both the generation capacity and the load level have the same graphic units (MW versus hours), it becomes easy to superimpose The load model is a graphic representation of the consumers' hourly chronological peak load for a period of one year. It, thus, consists of 8736 values which represent 8736 h of a year [11]. This is shown in Figure 6 through the red line in the graph of the superimposed generation and load model. where MTTR is the mean time to repair. The Equations (2) and (3) are suitable for describing state intervals of stochastic systems with any number of discrete states [11]. In our paper, they were calculated using the generation system data from the IEEE RTS '96 standard as shown in Appendix A in Table A1. It consists of nine (9) generation sources, each consisting of several generating units. The values of λ and µ are also obtainable from the IEEE RTS standard for each of the generating units via (4) and (5). The states' durations of each generating unit, i.e., TTF and TTR, are sequentially calculated as follows: (1) For each of the random numbers generated in U1, the TTF is calculated using (2); (2) For each of the random numbers generated in U2, the TTR is calculated using (3); (3) The resultant series of TTF and TTR are thereafter graphically and sequentially joined (one after the other) (i.e., TTF, then TTR, then TTF, and so on) as shown in Figure 5. This is continued until the total value reaches 8736 representing the total number of hours of a year. This is performed for each of the entire generating units. In the end, they are combined by a superimposition method to obtain a single resultant system capacity as shown in Figure 4. This is the annual duration of the generation capacity states, otherwise known as the generation model. The load model is a graphic representation of the consumers' hourly chronological peak load for a period of one year. It, thus, consists of 8736 values which represent 8736 h of a year [11]. This is shown in Figure 6 through the red line in the graph of the superimposed generation and load model. The risk model represents the superimposition of the generation and the load model from which the LOLE reliability index is calculated [20]. The LOLE is the sum of the individual algebraic differences between the load level and the generation capacity at the points where the load level exceeds the generation capacity. This is as shown by the two red shaded regions in Figure 6. This is usually calculated in a yearly period. The generation system has a record of the total capacity for each hour, hence 8736 values corresponding to 8736 h of the year. In addition, the load system has a record of the consumers' hourly load usage level, hence, also 8736 values. Since both the generation capacity and the load level have the same graphic units (MW versus hours), it becomes easy to superimpose The risk model represents the superimposition of the generation and the load model from which the LOLE reliability index is calculated [20]. The LOLE is the sum of the individual algebraic differences between the load level and the generation capacity at the points where the load level exceeds the generation capacity. This is as shown by the two red shaded regions in Figure 6. This is usually calculated in a yearly period. The generation system has a record of the total capacity for each hour, hence 8736 values corresponding to 8736 h of the year. In addition, the load system has a record of the consumers' hourly load usage level, hence, also 8736 values. Since both the generation capacity and the load level have the same graphic units (MW versus hours), it becomes easy to superimpose the load level on the generation capacity to calculate the resultant difference in their capacities (MW). The points Energies 2020, 13, 4350 6 of 13 (hours) where the generation capacity exceeds the load is regarded as the safe region, while, the points (hours) where the load level exceeds the generation capacity is recorded as the LOLE. The summation of these points (hours) of LOLE for a year gives the LOLE.

The Reliability Assessment and the Monte Carlo System
The LOLE reliability index was obtained from the reliability assessment via the MCS. The MCS method accommodates the natural stochastic factors in the generation system operating states considering the inherent fluctuations in the generating capacity [20][21][22]. The Monte Carlo model was developed to be dynamic in solving several kinds of problems instead of the problem having to match the model or the method [23,24]. This predicts the system's real-time behavioral characteristics in the simulated time and generates the probability distribution of each of the parameters [25,26]. The reliability assessment data were obtained from the IEEE reliability test system (IEEE RTS '96 standard). The generation and load data of the IEEE RTS are given in Appendices A and B, respectively. This was established in the IEEE power engineering society by the IEEE subcommittee on the application of probability methods [27,28]. The aim was to present a common test platform for the comparison of results obtained through various methods of power system reliability analysis. It is comprised of the following [11]: A. Total generation capacity [Appendix A- Table A1]  (1) The annual peak load of 2850 MW; (2) Weekly peak load presented in percent of the annual peak load. This comprises of 52 percentage values that represent the 52 weeks of a year (Table A2); Daily peak load presented in percent of the weekly peak load. This comprises of 7 percentage values that represent the 7 days of a week (Table A3); (4) Hourly peak load presented in percent of the daily peak load. This comprises of 24 percentage values that represent the 24 h of a day (Table A4).
The load composition of the IEEE-RTS is such that the load characteristics of every season of the year are considered. Thus, the inherent stochasticity in the load system is well accommodated in the standard.

Modeling the Generation Capacity and the Load
From the IEEE RTS standard data used in the simulation process, the derived generation model was superimposed with the load model and the risk model was consequently accessed. From the risk model, the LOLE reliability index was thereafter calculated. The simulation was performed in the MATLAB environment.
To further investigate the impact of the load growth on the available generation capacity adequacy, the load data were further increased stepwise and sequentially. They were increased from the original level (0%) to 30% fold by a tiny step size of 0.06%. Thus, the load multiplication factor ranges from 1.0 to 1.30 through a step size of 0.0006. In each step, the generation capacity model was simulated against the incrementing load and the LOLE was calculated as described in Figure 6. Hence, 501 LOLE values were recorded, i.e., (30%/0.06%) + 1 (the original load level).
In addition, in each step, the calculated LOLE value was translated into the "percentage of annum" metric, regarded as A%, described in (6), i.e., what percentage of annum is the calculated LOLE? The A% metric is introduced to quantify the realistic impact of each of the calculated LOLE values Energies 2020, 13, 4350 7 of 13 on the generation capacity adequacy. A yield point, defined by A y %, is thereafter established at a particular value of A% such that the values of A% that are above A y % are unacceptable. This is known as the reliability threshold of the generation capacity. It is the maximum accommodable limit of A% above which the generation system is pronounced unreliable. At this point, the simulation process is stopped, and the generation improvement is demanded. Hence, the values at A y % and below signify that the generation capacity is adequate and thus, reliable. The "percentage of annum" A% is described by: A% = (LOLE/8736)·100% (6) where LOLE means loss-of-load evaluation. Thus, if: A% < A y %, generation system is reliable, A% = A y %, generation system is reliable, and A% > A y %, generation system is unreliable.
A y % signifies the reliability threshold. The value of the A y % is determined by the power utility and it serves as a benchmark that defines her quality of service (QoS). A lower value of A y % implies higher QoS, thus, greater system reliability and vice versa. A flowchart for generation system reliability assessment and notification is shown in Figure 7. This depicts a smart approach to accommodate the load rise effect and timely generation system reinforcement through a notification alert system.
Energies 2020, 13, x FOR PEER REVIEW 7 of 13 known as the reliability threshold of the generation capacity. It is the maximum accommodable limit of A% above which the generation system is pronounced unreliable. At this point, the simulation process is stopped, and the generation improvement is demanded. Hence, the values at Ay% and below signify that the generation capacity is adequate and thus, reliable. The "percentage of annum" A% is described by: where LOLE means loss-of-load evaluation. Thus, if: A% < Ay%, generation system is reliable, A% = Ay%, generation system is reliable, and A% > Ay%, generation system is unreliable. Ay% signifies the reliability threshold. The value of the Ay% is determined by the power utility and it serves as a benchmark that defines her quality of service (QoS). A lower value of Ay% implies higher QoS, thus, greater system reliability and vice versa. A flowchart for generation system reliability assessment and notification is shown in Figure 7. This depicts a smart approach to accommodate the load rise effect and timely generation system reinforcement through a notification alert system.

Simulation Results and Discussions
The generation-load simulation process comprises of 501 LOLE values. An overview of the load increment at six major intervals (namely: 0%, 5%, 10%, 15%, 25% and 30%) against the available generation capacity is shown in Figure 8 in a generation-load graph. The graph of the load increment against the generation capacity could be seen ranging from (a) to (f) for 0% to 30% load increment, respectively. This is designated with distinct colors for clarity. The generation capacity is represented by the magenta color, while the various load increment is represented by other colors. The "load spikes" into the generation capacity region in (b) for a 5% increment is seen greater than the "load spikes" in (a) for a 0% load increment. Likewise, the "load spikes" in (c) for a 10% load increment is greater than that in (b) for a 5% increment. In addition, the spikes in (d) for a 15% increment is greater than that in (c) for 10% increment, and so on, until (f). These various "load spikes" into the generation

Simulation Results and Discussion
The generation-load simulation process comprises of 501 LOLE values. An overview of the load increment at six major intervals (namely: 0%, 5%, 10%, 15%, 25% and 30%) against the available generation capacity is shown in Figure 8 in a generation-load graph. The graph of the load increment against the generation capacity could be seen ranging from (a) to (f) for 0% to 30% load increment, respectively. This is designated with distinct colors for clarity. The generation capacity is represented by the magenta color, while the various load increment is represented by other colors. The "load spikes" into the generation capacity region in (b) for a 5% increment is seen greater than the "load spikes" in (a) for a 0% load increment. Likewise, the "load spikes" in (c) for a 10% load increment is greater than Energies 2020, 13, 4350 8 of 13 that in (b) for a 5% increment. In addition, the spikes in (d) for a 15% increment is greater than that in (c) for 10% increment, and so on, until (f). These various "load spikes" into the generation capacity are known as the LOLE for each load increment considered. They are the points where the consumer load surpasses the generation capacity. The yield point was established at A y % ≈ 15% (floor), thus, the simulation was stopped at this point. This was reached at approximately 30% load increment and the LOLE of 1380 h/yr, as seen in Figure 9. From the data sequence in the graphs of Figure 8, it is evident that the LOLE increases linearly with the load increment. Hence, a further simulation was performed to obtain the respective LOLEs from up to 500 steps of load increment. Figure 9 shows the graph of the respective LOLEs obtained from the generation-load model as a result of 500 points of load increment ranging from 0% to 30%.
Energies 2020, 13, x FOR PEER REVIEW 8 of 13 capacity are known as the LOLE for each load increment considered. They are the points where the consumer load surpasses the generation capacity. The yield point was established at Ay% ≈ 15% (floor), thus, the simulation was stopped at this point. This was reached at approximately 30% load increment and the LOLE of 1380 h/yr, as seen in Figure 9. From the data sequence in the graphs of Figure 8, it is evident that the LOLE increases linearly with the load increment. Hence, a further simulation was performed to obtain the respective LOLEs from up to 500 steps of load increment. Figure 9 shows the graph of the respective LOLEs obtained from the generation-load model as a result of 500 points of load increment ranging from 0% to 30%.  Since the relationship between LOLE and the load rise maintained a linear trend, a regression analysis is further introduced in the simulation process. Thus, the obtained LOLE values are simulated with the load rise through the least-squares method. The 501 simulation data points were split into 80% training and 20% testing datasets. The simulation was performed in the Python-Jupyter integrated development environment (IDE). The linear regression algorithm of Scikit-learn was used, and the random state of the train_test_split was seeded to 90 (random_state = 90). Through simulation, a relationship between the LOLE and the load rise was established. The fitting formula is given by: i LOLE = 45.843L + 5.296 where Li is load. Equation (7) depends on the uniqueness of the relationship, and values would vary based on the recorded generation and load data. The equation is wholly dependent on the load-growth history, Li of the particular consumers in the area under consideration. Also, the capacity and state of the available generation system is a contributing factor. Thus, different areas would record different simulation values. However, the relative variation of the independent parameter, Li, with the LOLE remains proportionally unique. In our simulation, the accuracy of the formula is tested with 101 data points of the test data, and the coefficient of determination (R-square) recorded is 0.99. Given the average deviation of the prediction in (8), the average deviation of the 101 data points is 0.63. The average deviation Dav is given by: where LOLEi is the actual LOLE. From the simulation result, the future generation system reliability assessments can be effectively performed using (7). From the generation system operation data and the recent loadgrowth history (data), a simulation is performed using the linear regression algorithm. Using the resulting fitting equation, the value of LOLE could be predicted at any instance given the current load information. The LOLE is then monitored as the load rises, by translating it to the A% equivalent as seen in (6). If A% exceeds the preset yield point (Ay%), known as the percentage of annum, the generation reinforcement is necessitated. At this point, the generation system reinforcement demand is declared critical. From our simulation, the tiny load increment step size of 0.06% adopted ensures high precision in the reliability assessment scheme. Thus, this guarantees high accuracy in the prediction through (7). This, therefore, would accommodate the smallest possible load growth through the robust feature. Since the relationship between LOLE and the load rise maintained a linear trend, a regression analysis is further introduced in the simulation process. Thus, the obtained LOLE values are simulated with the load rise through the least-squares method. The 501 simulation data points were split into 80% training and 20% testing datasets. The simulation was performed in the Python-Jupyter integrated development environment (IDE). The linear regression algorithm of Scikit-learn was used, and the random state of the train_test_split was seeded to 90 (random_state = 90). Through simulation, a relationship between the LOLE and the load rise was established. The fitting formula is given by:

Conclusions
where L i is load. Equation (7) depends on the uniqueness of the relationship, and values would vary based on the recorded generation and load data. The equation is wholly dependent on the load-growth history, L i of the particular consumers in the area under consideration. Also, the capacity and state of the available generation system is a contributing factor. Thus, different areas would record different simulation values. However, the relative variation of the independent parameter, L i , with the LOLE remains proportionally unique. In our simulation, the accuracy of the formula is tested with 101 data points of the test data, and the coefficient of determination (R-square) recorded is 0.99. Given the average deviation of the prediction in (8), the average deviation of the 101 data points is 0.63. The average deviation D av is given by: where LOLE i is the actual LOLE. From the simulation result, the future generation system reliability assessments can be effectively performed using (7). From the generation system operation data and the recent load-growth history (data), a simulation is performed using the linear regression algorithm. Using the resulting fitting equation, the value of LOLE could be predicted at any instance given the current load information. The LOLE is then monitored as the load rises, by translating it to the A% equivalent as seen in (6). If A% exceeds the preset yield point (A y %), known as the percentage of annum, the generation reinforcement is necessitated. At this point, the generation system reinforcement demand is declared critical. From our simulation, the tiny load increment step size of 0.06% adopted ensures high precision in the reliability assessment scheme. Thus, this guarantees high accuracy in the prediction through (7). This, therefore, would accommodate the smallest possible load growth through the robust feature.

Conclusions
Generation system reliability assessment is demanding and necessitates that it is performed at regular intervals. The system must be monitored to ascertain the reliability status of the existing generation system at any point in time to ensure continued adequacy in supply. A reliability predictive analysis is made using the generation system data and consumer load-growth history. The generation system reliability assessment propelled by such a predictive approach would facilitate timely response to the reliability demand of deficient generation capacity. Through the entire simulation approach, the LOLE is carefully bounded below the yield point. This is regarded as the safe region within which the generation system operation capacity is endorsed reliable. By this approach, the adequacy of the available generation system is predictively guaranteed. The yield point metric (A y ) would also serve as a basis to guide the consumers in the decision making process. The patronage-related decisions arising from the course of market competitions would rightly be made given the various available utility companies. The power utilities operating at lower A y % have a higher quality of service (QoS). They would certainly earn greater consumers' patronage in the electricity market. The entire system also promotes an optimal balance between system reliability and the cost of power supply.