**World Journal of Engineering and Technology**

Vol.07 No.03(2019), Article ID:94440,26 pages

10.4236/wjet.2019.73032

Availability, Performance and Reliability Evaluation for PV Distributed Generation

Abdulrahman K. Al-Sefri, Abdullah M. Al-Shaalan^{ }

King Saud University, Riyadh, KSA

Copyright © 2019 by author(s) and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: June 21, 2019; Accepted: August 17, 2019; Published: August 20, 2019

ABSTRACT

Nowadays, renewable energy resources play an important role in replacing conventional energy resources such as fossil fuel by integrating solar, wind and geothermal energy. Photovoltaic energy is one of the very promising renewable energy resources which grew rapidly in the past few years, it can be used to produce electric energy through photovoltaic process. The primary objective of the research proposed in this paper is to facilitate the increasing penetration levels of PV systems in the electric distribution networks. In this work, the PV module electrical model is presented based on the mathematical equations and was implemented on MATLAB to simulate the non-linear characteristics I-V and P-V curves with variable input parameters which are irradiance and temperature based on Riyadh region. In addition, the reliability evaluation of distribution networks, including distributed generators of solar photovoltaic (PV) with varying output power capacity is presented also. The Monte Carlo simulation algorithm is applied to test the distribution network which is RBTS Bus 2 and the same has been conducted on the original case of distribution network substation 7029 which is located at KSA Riyadh. The two distribution networks have been modified to include the PV’s distributed generators. The distributed generators contribute to supply a part of the load during normal mode and supply the entire load during component failure or failure of grid operation supply. The PV stochastic models have been used to simulate the randomness of these resources. Moreover, the study shows that the implementation and integration of renewable resources as distributed generations have improved the reliability of the distribution networks.

**Keywords:**

Renewable Resources, Photovoltaic (PV) Energy, Distributed Generation (DG), Distribution Networks

1. Introduction

The electric power system has been used to generate, transmit and distribute power since the advent of electricity evolution and the foremost important objective of an electrical utility is to deliver economical, reliable, and quality power to its consumer. The electrical power network is extremely complex, a failure may result in loss of power to a large number of customers or sometimes catastrophic events such as blackouts and it is difficult to analyze the entire system at once. For this reason, the power system is divided into three functional zones such as generation facility, transmission facility and distribution facility to evaluate the reliability of the system. These functional zones in series can be considered as the hierarchical levels of the power system reliability studies as shown in Figure 1 [1]. Hierarchical level I deals with the reliability evaluation of only generation system and the main indices analyzed in this study are loss of load expectation (LOLE), loss of energy expectation (LOEE), failure frequency and failure duration. Hierarchical level II is composed of generation and transmission facilities and it deals with the reliability evaluation of both generation and transmission systems. Hierarchical level III includes all of the three functional zones and refers to the complete electric power system. It deals with the reliability evaluation of the generation, transmission and distribution systems.

The conventional structure of electrical power systems has been developed mainly to become as following arrangement shown below in Figure 2. The electric power is generated in large generating plants at a relatively small number of locations which will here be called central stations.

The incorporation and integration of non-conventional or renewable energy sources in the network result in a new term called Distributed Generation (DG). DG can be defined as small-scale power generation units of electricity (a few kilowatts kW) connected directly to the grid, distribution network, on consumer side of the meter to serve a customer on site and at the same time to provide the support to distribution networks or work standalone with different rating levels. Table 1 presented the generation of DG at different levels [2].

The integration of DGs becomes the most economical solution to meet the increased demand due to load growth in the conventional system. While providing environmentally friendly energy and while helping to meet the increasing load economically. It can play an important role in electrical power network

Figure 1. Functional zones and hierarchical levels of power system.

Figure 2. Conventional large electric power system.

Table 1. Generation of DG at different level [2].

and by employing this technology in power systems has big advantages as following:

· Improving the reliability of the electrical network by providing an alternative source during power disturbance events.

· Reducing the electrical transmission power losses caused by the power traveling through long transmission lines and high voltage transformers.

· Providing better voltage support, used to supporting peak load and supply electricity during peak periods.

· Improving the power quality and enabling consumers to select the source of energy based on the cost and awareness of the environmental issues.

· Relieving the congested of electrical transmission networks and reducing the need to expand the electrical transmission networks.

· DGs considered as a viable alternative of energy storage alternative when interruptions are frequent. It works as Backup supply to ensure the uninterrupted electricity.

Currently, industrial countries generate most of their electricity in large centralized facilities, such as fossil fuel powered by coal or gas, nuclear, large solar power plants or hydropower plants. However, modern embedded systems can provide these traits with automated operation and renewables, such as sunlight, wind, and geothermal. Moreover, the last few years, a number of factors have led to an increased interest in distributed generation schemes such as availability of modular generating plant, ease of finding sites for smaller generators, short construction time and lower capital costs for smaller plant, reduction in gaseous emissions (mainly CO_{2}). Due to the DG sources have a potential solution for some issues, like the deregulation in power system and improving the performance of distribution system as well as reliability evaluation, that is making it more popular. The presence of the DGs, especially when the DG share is significantly high, will obviously impact on way of power system operation. The distribution networks can be designed for (radial) unidirectional power flow, and the electric power systems with DGs spread across the distribution network is shown in Figure 3.

In recent years there has been an increase in the global trend in general and particular in Kingdom of Saudi Arabia (KSA) towards finding and exploiting sources of renewable energy to generate electricity. This may be due to the urgent need for new and renewable sources of energy that are not fading and depleting. One of the reasons for this is that these sources limit the emission of toxic gases as well as saving energy costs and conserving their stocks from fossil sources (oil and natural gas) and not relying on it total. Due to the high prices of fossil fuels with the rapid rises in oil prices, capital costs of conventional power plants, escalation in power demand due to rapid growth in population and industrialization in the last several years especially in Saudi Arabia.

The conventional generation, low quality fuels and the methods of generation typical in KSA is a chief cause of environmental pollution and impacts human health through emissions of harmful gases that are remains a threat to public health around the world. Therefore, it is essential to find green power source that reduce the gasses emission and preserves oil in Saudi Arabia. So, the renewable energy

Figure 3. Electric power system with the presence of DGs.

resources have been globally accepted for power generation, utilizing and integrated them with conventional generation around the globe.

The solar radiation is variable in different parts of the world and the highest value of solar radiation located at Sunbelt. Saudi Arabia is most suitable country over the world for using solar power and its solar irradiances are among the highest in the world because it located in the “Global Sunbelt”, a geographic region situated between 35˚N and 35˚S. The generally characterized by high solar irradiation which can be witnessed in Figure 4.

KSA has widespread desert land, year-round clear skies and abundant solar resource which it to become one of the largest solar photovoltaic (PV) energy producers in future and a world leader in renewable electricity generation. Studies show that the use of solar equipment in KSA is very suitable and can easily provide part of the energy that nation needs.

The Several sources, including Snapshot of Global PV Markets 2016 [4], state that the total installed PV capacity at the end of 2016 was at least 303 GW. In addition, there is an estimated capacity of between 3000 GW and 10,000 GW expected to be installed by the year 2030 [5]. As part of the kingdom’s 2030 vision, to build the largest PV generation plant over the world, and KSA aims to cut down overreliance on oil and gas as the principal sources of energy in the generation of electricity and KSA government is backing the process of energy source diversification to keep the oil for next generations. The KSA targeting to generate 3.45 GW of renewable energy by 2020 and 9.5 GW by 2022 [6]. In recent years the KSA was is ranked as the 12^{th} largest consumer of total primary energy in year 2013, of which about 60% is petroleum based and the rest was accounted to natural gas. Energy consumption per capita is twice that of Europe and three times the world average [6]. By the Year 2032, the electricity peak demand is anticipated to surpass 120 GW [7]. The King Abdullah City for Atomic and Renewable Energy (KACARE) program seeks to ensure that half of the electricity generated in Saudi Arabia comes from renewable sources by 2032 [8].

Figure 4. Best countries for using solar power source [3].

2. Solar Radiation

The solar resource and weather conditions of a place are often intermittent and prone to vary from one year to another year, so the datasets need to be compiled over a number of years to create historical site data. The statistical model of solar radiation is based on data for the solar irradiance at KSA-Riyadh, during period Jan 2013-Jan 2016 which is available in King Abdullah City for Atomic and Renewable Energy’s (K.A.CARE’s) Renewable Resource Monitoring and Mapping (RRMM). Figure 5 shows the calculated distribution of the hourly average global, direct, and diffuse Solar Irradiance. Figure 6 shows the calculated monthly distribution average GHI of solar irradiance throughout the past time frame in Riyadh city.

Figure 5. Hourly average DNI, DHI and GHI of solar irradiance (Wh/m^{2}).

Figure 6. Monthly average GHI of solar irradiance (Wh/m^{2}).

The solar radiation in Riyadh is characterized by long radiation period with high irradiance during summer and short radiation period with low irradiance during winter due to the differences in the sun’s height at the summer and winter solstice. The solar irradiance is fluctuated from month to month based on the sun orbit. The highest value of monthly solar irradiance was for May, Jun, and July with maximum value 8134 (Wh/m^{2}) for Jun.

3. DG Model

It is well known that due to the fact that solar insolation is intermittent and the output power of PV systems is not deterministic. For this reason need for a stochastic model to simulate PV outputs. A stochastic model is a simulation-based technique to describe a non-deterministic behavior and the randomness of the system. The statistical data of solar insolation and temperature have been collected from (K.A.CARE’s) and implemented in our study.

a) Photovoltaic Power System Generation

The Solar photovoltaic (PV) is widely used as renewable energy source. The solar of PV has a high reliability in modules (>20 years), and high public acceptance. PV is the technology that generates direct current (DC) electrical power measured in Watts (W) or kilo Watts (kW) from semiconductors when they are illuminated by photons. As long as light is shining on the solar cell, it produces both a current and a voltage to generate electric power which is known as Photoelectric Effect [9]. A general representation of a grid-connected solar PV system is shown in Figure 7. The most essential components in a grid-connected solar PV system are solar PV generator which consists of a large number of solar arrays, DC-DC converter, inverters, filter and step up transformers. A solar cell module produces only a small amount of current and voltage, and in order to produce a large amount of electric power, the solar cell modules are connected into arrays. The output voltage from a PV array changes with solar radiation and ambient temperature. Moreover, to connect the PV system to the transmission

Figure 7. The schematic drawings of PV power system component.

grid, the output DC voltage from PV system should be first regulated by a DC-DC converter. Then DC-AC inverter converts the DC power produced by the PV modules into AC power. A filter is used to eliminate harmonics. The power electronic components have the tasks to guarantee safe and efficient operation, to track the maximum power point (MPP) of the PV system, and to maintain power quality of the PV system output. The transformer is an essential when the generated electricity is to be injected into the utility grid.

b) The PV Module

The ideal model of a solar cell consists of a photocurrent source and diode with an internal resistance connected to a parallel resistor. The current representing the photo generated current (I_{gc}) in parallel with a diode which is directly proportional to the light falling (solar irradiation) on the cell. A PV cell is physically modelled using the equivalent circuit model as represented in Figure 8.

Many researchers have been examined the general mathematical model description of a solar cell current-voltage output characteristic over the past three decades [10]. According to well know Kirchhoff’s current law, the current flowing in the circuit can be modelled by the following Equation (1) [11].

${I}_{PV}={I}_{sc}-{I}_{o}={I}_{sc}-{I}_{D}\left[{\text{e}}^{\left(\frac{q{V}_{D}}{nK{T}_{cell}}\right)}-1\right]-\frac{V+I{R}_{S}}{{R}_{P}}$ (1)

where: I_{PV} is the PV output current. I_{sc} is the cell’s short circuit current at 25˚C and 1 KW/m^{2}. q is the electron charge equal to 1.602 × 10^{−19} C. k is the Boltzmann constant (1.381 × 10^{−23} J/K). V_{D} is the voltage across the diode terminals. V is the output voltage. R_{S} and R_{P} are the series and parallel resistances. n is the ideality factor, also known as the quality factor or sometimes the emission coefficient and usually takes values in the range 1 to 2, which depends on the construction and semiconductor material listed in Table 2.

The cell’s saturation current (I_{o}) varies with cell’s temperature and depends on the cell temperature. As described by following Equations (2)-(5).

${I}_{gc}=\left[{\mu}_{sc}\left({T}_{cell}-{T}_{r}\right)+{I}_{SC}\right]G$ (2)

Figure 8. The PV equivalent circuit model.

Table 2. Ideality factor n dependence on PV technology [12].

${I}_{D}={I}_{o}\left({\text{e}}^{\frac{q{V}_{D}}{nK{T}_{cell}}}-1\right)$ (3)

${I}_{o}={I}_{o\alpha}{\left(\frac{{T}_{cell}}{{T}_{r}}\right)}^{3}{\text{e}}^{\left[\frac{e{v}_{d}}{KF}\left(\frac{1}{{T}_{r}}-\frac{1}{{T}_{c}}\right)\right]}$ (4)

${I}_{o\alpha}=\frac{{I}_{SC}}{{\text{e}}^{\left(\frac{e{v}_{o\alpha}}{KF{T}_{c}}\right)}}$ (5)

where: I_{D} is the diode current. T_{cell} is the cell’s absolute temperature. T_{r} is the cell’s reference temperature. I_{oα} is the cell’s reverse saturation current at solar radiation and reference temperature. G is the solar radiation in kW/m^{2}.

A solar cell alone can produce only output power from 1 to 2 Watt [13]. Furthermore, to increase the output power, the cell connected in series-parallel configuration on a module then array. A PV array is a group of several PV module and these arrays can be of sizes ranging from a few hundred watts to hundreds of kilowatts. The interconnection of PV modules in series and parallel or various other connections forms a PV array. The PV array composed of N_{S} series and N_{P} parallel connected modules. And in order to calculate the generated photovoltaic current and the output power of a PV module consisted of (N_{s} × N_{p}) cells as following Equation (6).

${I}_{PV}={N}_{p}{I}_{gc}-{N}_{p}{I}_{o}\left[{\text{e}}^{\frac{e\left(\frac{{v}_{d}}{{R}_{s}}+s\frac{I{R}_{s}}{{R}_{p}}\right)}{{N}_{s}KF{T}_{c}}}-1\right]-\frac{\frac{{N}_{p}{V}_{d}}{{N}_{s}}}{{R}_{p}}$ (6)

c) Non-Linear Characteristics of PV’s

The Photovoltaics have nonlinear characteristics, where the performance and output power are directly affected with the change of the operating conditions solar irradiance, temperature and the angle of the sun. Usually the PV manufacturer supply their products with a data sheet that contains values of current and voltage for three conditions namely, the short circuit, the open circuit, and the maximum power for a given set of reference condition. The major inputs for the proposed PV model were solar irradiation, PV panel temperature and PV manufactures data sheet information. In this study, the Astronergy CHSM6610P Polycrystalline PV Module is taken as example to simulate it. The data sheet of the simulated module is shown in Table 3.

An insolation-oriented model of PV module is built by using MATLAB to illustrate and verify the nonlinear I-V and P-V output characteristics of PV module. Both I-V and P-V output characteristics of PV module at various insolation and temperatures based on Riyadh region data are carried out. Figure 9 and Figure 10 show the effect of changing the temperature and solar irradiance on PV’s output current, voltage, and power as well.

It is clear from the previous figures that the output power of PV’s is directly proportional with the amount of solar irradiance falling on PV’s panel, and inversely proportional with its temperature. We can see the effect of change in solar irradiation on PV characteristics from Figure 9. It is illustrated that as we increase the solar irradiation short circuit current increases, and current increases from

Table 3. Astronergy CHSM6610P polycrystalline specifications.

Figure 9. The irradiance changing effect on I-V and P-V output characteristics.

Figure 10. The temperature changing effect on I-V and P-V output characteristics.

1.8 A to 8.3 A approximately. The effect of variation of solar irradiation on P-V characteristics result that as solar irradiation increases, power generated increases also. Increase in power is mainly due to increment in current. Decreasing the irradiance will reduce the overall performance of the PV module. On other hand, we observe from Figure 10 that, with an increase in the temperature the maximum power output of P-V curves are decreases whilst the short circuit current decrease also. As the temperature increases, the efficiency of the panel decreases. We can see the variation of temperature on I-V curves mostly effects voltage, as we increase the temperature, voltage decreases but current remains almost unaltered.

d) The PV Output Power Model

The PV output power system was proposed to assess the reliability of distribution system containing renewable energy resources (solar energy generation plants with different capacity). The PV system output power can be calculated by the following Equation (7) [14] :

$P\left(out\right)=\{\begin{array}{l}\frac{{\eta}_{C}}{K}\ast S\ast I{\left(t\right)}^{2}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}0<I\left(t\right)\le K\\ {\eta}_{C}\ast S\ast I\left(t\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}I\left(t\right)>K\end{array}$ (7)

where ( ${\eta}_{C}$ ) is the efficiency of the PV system, K is a threshold and $I\left(t\right)$ is the hourly solar insolation. The value of ( ${\eta}_{C}$ ) is not constant when $I\left(t\right)$ is less than or equal, it is proportional with the solar insolation. That means when the solar insolation increase the PV system efficiency will increase and vice versa. The following equation can express the hourly solar insolation (8) [12].

$I\left(t\right)=\{\begin{array}{l}{I}_{\mathrm{max}}\left(-\frac{1}{36}{t}^{2}+\frac{2}{3}t-3\right)\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}6\le t\le 18\\ 0\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}\text{\hspace{0.17em}}0\le t<6\text{\hspace{0.05em}}\text{\hspace{0.17em}}\text{and}\text{\hspace{0.17em}}\text{\hspace{0.05em}}18\le t\le 24\ge 0\end{array}$ (8)

4. Load Model

The load modelling is an important part of power system modelling process and has a significant effect on power system simulation results. The load model is created by using daily or monthly or yearly peak loads with respect to time in seconds or minutes or hours. The energy consumption can be displayed with load curve as shows the amount of demand used per period of time but Load Duration Curve (LDC) shows the amount of maximum ranking until the minimum demand is used. LDC is the basic tool used in the analysis of electric power systems such as estimating the operating cost of resource plans, and as tool to integrate the demand side management in the planning of electricity generation and enables to evaluate the operation of the power system more accurately than the traditional approach. Figure 11 shows LDC displaying various load-related variables. This load duration curve has a shape of a monotonously decreasing function and in this load duration curve, we can know that the slope varies according to time.

Over the last several years, the peak of electricity demand in Kingdom of Saudi Arabia (KSA) is continuously increasing. The Saudi Electricity Company (SEC) operates an interconnected transmission system for all main areas in KSA. SEC system is divided into four operating areas, named as Central, Eastern, Western and Southern Operating areas. Since Riyadh city is the major load center in Central operating area it faces unique issues in every day operation of the system [15]. Figure 12 and Figure 13 show the measured average values of daily and monthly electrical energy consumption in Riyadh city in the year 2018. It can be observed that the summer months from June to October have consistently higher loads than the non-summer months.

The behavior of power system loads is a frequent pattern during normal conditions, and a time varying load model can be developed by using historical data. The model used in this study is the hourly load curve which can be converted to load duration curve by arranging the data in descending order. The following equation used to find the predicated the load for any load point at any desired time (9).

Figure 11. Load duration curve displaying various load-related variables.

Figure 12. Variation of daily average load demand of Riyadh city in the year 2018.

Figure 13. Variation of monthly average load demand of Riyadh city in the year 2018.

${P}_{i}={W}_{h}\left(h\right)\times {W}_{m}\left(m\right)\times {P}_{Li}$ (9)

where, ${W}_{h}\left(h\right)$ = hourly weight factor, ${W}_{m}\left(m\right)$ = monthly weight factor, and ${P}_{Li}$ = peak load for load point.

5. System Reliability Indices

The basic function of an electric power system is to supply customers with reasonably economical and reliable electricity. The reliability level of one distribution system is often evaluated from two aspects which are consumer level and system level. These two frequently used reliability indices for users of load points are as follow. The first one is average failure rate (λ_{i}) which is number of outages/year and defined as the probability of failure occurrence during a specific period for load point. The second one is average annual interruption time (U_{i}) which is sum of the outages time/year and defined as the average interruption time of load point in a specific period. These two indices can be expressed as following Equations (10)-(11) [16].

$\lambda ={\displaystyle \sum _{i=1}^{n}{\lambda}_{i}}$ (10)

$U={\displaystyle \sum _{i}^{n}{\lambda}_{i}{r}_{i}}$ (11)

where λ_{i} is the failure rate of the series components from the source point to load point, n is the total number of components which affect load point, r_{i} is average restoration time of network component to restore load point i due to the failure of component j.

The commonly reliability induces of the system defined as functions of average failure rate over the total number of customers, and average interruption time [16]. The reliability system indices measure by some factors which can be aggregated to provide an appreciation of the system performance using a series of system indices such as System Average Interruption Frequency Index (SAIFI) which refers to the ratio of the total number of customer interruptions to total number of customers served and is the number of permanent interruption that customers would experience in one year. System Average Interruption Duration Index (SAIDI) which refers to the ratio of the sum of customer interruption durations to total number of customers and usually measures the duration of permanent interruption that customers would experience in one year. Customer Average Interruption Duration Index (CAIDI) which refers to the ratio of the sum of customer interruption durations to total number of customer interruptions and is the average time needed to restore service to the average customer per sustained interruption. Expected Energy Not Supply (ϵENS) measures the total of energy interruption that customers would experience in one year and its dimension (MWh/year). These indices can be calculated by following Equations (12)-(14) [16] [17].

$\text{SAIFI}=\frac{{{\displaystyle \sum}}_{i=1}^{n}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\lambda}_{i}{N}_{i}}{{{\displaystyle \sum}}_{i=1}^{n}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{N}_{i}}$ (12)

$\text{SAIDI}=\frac{{{\displaystyle \sum}}_{i=1}^{n}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{U}_{i}{N}_{i}}{{{\displaystyle \sum}}_{i=1}^{n}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{N}_{i}}$ (13)

$\text{CAIDI}=\frac{{{\displaystyle \sum}}_{i=1}^{n}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{U}_{i}{N}_{i}}{{{\displaystyle \sum}}_{i=1}^{n}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\lambda}_{i}{N}_{i}}=\frac{\text{SAIDI}}{\text{SAIFI}}$ (14)

$\u03f5\text{ENS}={\displaystyle \sum}_{i=1}^{n}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{E}_{i}{N}_{i}$ (15)

where N_{i} is number of customers at load points i, U_{i} is the annual outage time, N_{i} is the number of customers at load point i, E_{i} is the average of average interruption energy per load point.

6. Monte Carlo Simulation

Monte Carlo Simulation (MCS) has been used in reliability evaluation of distribution systems to simulate the failures due to that failures in power system networks are random in nature. It is a powerful method for solving a complex system and often used in complex mathematical calculations, stochastic process simulations, engineering system analysis and reliability calculation. MCS is a probabilistic method that can be used to predict the behavior of the system components which in real time will be all different in varying degrees, including the number of failures, times to failure, restoration times. The Simulation techniques used to estimate the reliability indices by simulating the actual process and random behavior of the system. The time sequential simulation is one of the MCS types used when the system behavior depends on past events to examine and predict behavior patterns in order to obtain the probability distributions of the various reliability parameters. The system reliability indices can be obtained from the artificial history that needed in time sequential simulation, and this can be obtained by generating the up (where expressed the element is in the operating state) and down (where expressed the element is inoperable due to failure) times randomly for the system elements. These artificial histories depend on the reliability parameters the elements and the system operating/restoration modes.

Time to failure (TTF) or failure time (FT) is the duration that it would take the component to fail or the time during the element remains in the up state. This time is predicted randomly and calculated by the following Equation (16) [16] [17].

$\text{TTF}=-\frac{1}{\lambda}\times \mathrm{ln}\left(n\right)$ (16)

where λ is failure rate of system component and n is a random number (range from 0 to1).

Time to repair (TTR) or time to replace (TTR) is the time required to repair a failed component or the time during the element remains in the down state and it is used to indicate the cycle time of failure. Also, this time is predicted randomly by the following (17) [16] [17].

$\text{TTR}=-\frac{1}{\mu}\times \mathrm{ln}\left(n\right)$ (17)

where μ is repair rate of system component.

The process of transiting from the up state to the down state is the failure process and can be caused by the failure of an element or by the removal of elements for maintenance time. The parameters TTF and TTR are random variables and it is obvious from Equations (14) and (15) that TTF and TTR follow exponential distributions. So, to predict the artificial history of system components, TTF and TTR can be generated to cover simulation times (e.g. one year) in chronological order. MCS have to be performed for a large number of scenarios, and the simulation time can be expanded to be a very long time (e.g. a thousand years or more) depending on the case study and also the desired accuracy. After that, the average can be calculated. Figure 14 shows an example of simulated component operating/restoration history of a component.

The element failures might have influence on one or additional load points and the most difficult in the simulation is to search out the load points and their failure period that affected by the failure of an element which are dependent on the network configuration. And in order to get a clear vision of how adding DGs impact on the reliability evaluation results of networks. We have used two reliability

Figure 14. Element operating/repair history [16].

evaluation algorithms which are MV network that does not contain DGs of PV and a MV network that containing DGs of PV.

As the main purpose of this study is to evaluate the reliability of MV network system with renewable energy resource which is the PV distributed generation with vary output power capacity to see their impact on the reliability, the following assumptions were created that shouldn't have a major impact on the results:

1) The permanent faults only are included in the study and all protection devices operate with success to isolate the faults.

2) The primary main feeder failures only are included in the analysis and every section protected by a breaker to isolate the faults.

3) It takes one hour to transfer the loads from the failed feeder to a neighboring feeder through a normal operating point.

4) Each circuit breaker is controlled by a bi-directional protection device.

The simulation procedures of reliability evaluation of distribution networks do not contain and contain renewable energy of PV distributed generators are shown in Figure 15 and Figure 16. It shows the algorithms of distribution system reliability evaluation indices without and with PVs distributed generators respectively.

7. Case Study Application

Many researchers have been studied the distribution network reliability and used the RBTS Bus 2 or Bus 4. Moreover, substation number 7029 MV distribution network in Riyadh region has selected similarly to conduct the reliability evaluation study and compare the results with RBTS Bus 2. These networks offer the information needed to conduct a reliability study. The RBTS bus 2 is supplied by two 33/11 KV, 16 MVA transformers whereas the substation 7029 is supplied by two 33/13.8 KV, 20 MVA transformers. The 0.415 KV and 0.400 KV low voltage customers are supplied via 11/0.415 KV and 13.8/0.400 KV transformers, whereas the 11 kV and 13.8 KV customers are supplied directly. For the reliability analysis the 33 kV supply has been considered 100% reliable and the feeders are operated as radial feeders but connected as a mesh through normally open sectionalizing points. Following a fault on a feeder, the ring main units permit the sectionalizing

Figure 15. Flowchart of distribution system reliability evaluation without DGs of PV.

Figure 16. Flowchart of distribution system reliability evaluation with DGs of PV.

point to be moved and customers to be supplied from alternative supply points. Figure 17 and Figure 18 show the modified single line diagram for the distribution networks system of RBTS bus 2 and substation 7029 with DGs of PV.

The data of RBTS Bus 2 distribution network may be found in reference [18]. The modified RBTS and substation 7029 have one type of distributed generator, is photovoltaic solar energy with vary of output power capacity. The maximum output power of PV1, PV2 and PV3 is 200 KW, 400 KW and 800 KW. The failure rate of PV1, PV2 and PV3 generators are taken as 0.11, 0.15 and 0.21 f/yr respectively. The average repair time for DGs of PV is taken as 32, 31 and 30 hr respectively. The all distributed generators of RBTS and substation 7029 contain all feeders (1, 2, 3 and 4) and are supplied by PV1, PV2 and PV3, except feeder 2 of RBTS and substation 7029 due to this feeder has less numbers of load points. The distributed generators of PV can operate in islanded mode, if there is enough

Figure 17. Distribution system for RBTS Bus 2 with PV DG.

Figure 18. Distribution network for MV substation 7029 with PV of DG.

power generation to supply the feeders loads during a permanent fault. We have used MCS in this work to evaluate the reliability indices of RBTS bus 2 and substation 7029 without DGs of PV and with DGs of PV. In addition, the impact of implementing these technologies in the reliability assessment of the distribution network is investigated. The annual average rate, average interruption time, and annual interruption energy of load points for RBTS and substation 7029 in two cases without and with DGs of PV are shown in Figures 19-21 respectively.

Figure 19. Annual average failure rate for all load points of RBTS and S/S7029.

Figure 20. Average interruption time for all load points of RBTS and S/S7029.

The RBTS and Substation 7029 failure rate of the feeder section is a function of the length. However, each section has a different length which resulted in the failure rates of the short feeder (F2) are low compared to the long feeders (F1, F2 and F4) of RBTS network. In addition, the failure of loads on the short feeder (F2) is also low compared to loads on the long feeders (F1, F3, F4) in the substation 7029. In the RBTS bus 2 and substation 7029 without distributed generators, it is obvious that the load points located at the end of the main feeder have high failure rates compared to the loads in the beginning of main feeder, because permanent faults result in isolating these load points from the main source. On

Figure 21. Average interruption energy for all load points of RBTS and S/S7029.

the other hand, the load points with distributed generators of PV located at the end of the main feeder have low failure rates. This reduction in failure rate is due to the excess generation capacity provided by the DGs during the outage of the main sources. From Figure 19, it is obvious that there is no reduction and improvement in the average failure rates for load points 1, 2, 8, 10, 11 and 16 in the RBTS and load points 1, 9, 13, 25 and 26 in the substation 7029 since a probabilistic simulation was used in this study.

The average interruption time is a function of failure rates and average restoration time of RBTS and substation 7029 network components. Therefore, the all previous observations are applicable to Figure 20. As a result, it is observed from Figure 20 that RBTS and substation 7029 with DGs of PV can reduce the average interruption time for most of load points slightly. Comparing the result between RBTS Bus 2 and substation 7029 it demonstrated that the average interruption time of substation 7029 is higher than RBTS Bus 2 due to has a greater numbers of load points and long lines. The calculated average interruption time for all load points of substation 7029 without and with DGs of PV is 32.305 per hour and 30.556 per hour. Whilst, the RBTS Bus 2 has total of average interruption time without and with DGs of PV are 17.847 per hour and 17.309 per hour.

The average interruption energy of all load points is described as a function of the average interruption time and the load of each load point. Therefore, all of the previous observations and clarifications are applicable to Figure 21. It is obvious that most of the load points have a little improvement in interruption energy after adding DGs of PV.

Table 4 shows the results of the overall reliability indices of the entire system for RBTS Bus 2 and substation 7029 for the two cases (without and with DGs of PV). All cases are using the method proposed in this study.

It is observed that from the previous study and results, the implementation of DGs of PV in distribution networks can improve the reliability of distribution networks by offering a backup source when the main source is not available. In addition, the overall reliability indices have improved after adding DGs of PV which is shown in Table 4. Based on the previous Figures 19-21 we observe that the reliability improvement associated in the presence of network with distributed generators not only depend on the size of DGs but the location of DGs is play an important role also. The DGs of PV integration have a positive impact on whole system reliability irrespective of the penetration level or the existing conditions.

The summary of these system reliability indices (SAIFI, SAIDI, CAIDI and ϵENS) are shown in Figure 22. SAIFI has least impact and lowest interruptions per customer on RBTS Bus 2 distribution network. The reliability indices of distribution networks (SAIDI) affected very slightly after adding DGs of PV. RBTS Bus 2 has the highest impact on system reliability indices (CAIDI) than S.S 7029.

Table 4. Reliability indices of RBTS Bus 2 and substation 7029.

Figure 22. Comparison of system reliability indices for two distribution system networks.

As the faulted feeder and all load points suffer power outage when a fault occurs either at any location between feeders and main bus or at any feeder circuit breaker or on the main bus. Thus, it has highest impact on ϵENS. We observed that, there is a slight improvement after adding DGs of PV for two distribution networks. RBTS Bus 2 has better system reliability (ϵENS) than S.S 7029. These reliability indices (SAIFI, SAIDI, CAIDI and ϵENS) can be used to determine the system performance of the composite distribution system. The system reliability indices can also be used to make assessment for the severity of system failures on future reliability analysis.

8. Conclusions

The implementation of renewable resources in distribution networks is promising in many environment and economic aspects, such as reducing green gas emissions, reducing the power losses on distribution networks and improving reliability of power services. The distributed generators of PV could pave the way to integrate solar energy in distribution systems, which can deal with different modes of operation such as interconnected mode and islanded mode. The excess of PV generation does not lead to any negative impact on the system reliability, but they surely have significant economic impacts. Due to high correlation of PV output power during the peak hours of the load, the PV integration results in reduced net loads which lead to higher effective load carrying capability.

In this study, the reliability evaluation was conducted on a test distribution network which is RBTS Bus 2. The data of RBTS Bus 2 have been used in this study because it offers very detailed information for each distribution components (transformers, breakers, lines and busbars) such as the failure rate and average restoration time. The same study has been conducted on the original case of distribution network substation 7029 which is located at KSA Riyadh. Since the original RBTS Bus 2 and substation 7029 did not contain of any distributed generators, the two systems have been modified to include the PV’s with different power output capacity. The probabilistic techniques used to evaluate the reliability of distribution networks containing PV as a distributed generator which is random in its output power production. MCS is a very useful tool, it requires only basic information to generate the artificial history of the distribution system components by generating large and random numbers of scenarios of the system element. The failure rates and average restoration times can be easily found by calculating the frequency and the duration of the down time of each load point in the system. The failures rates and average restoration times were assigned to all distributed generator resources, these indices were considered in reliability assessment studies. The common indices such as SAIFI, SAIDI, CAIDI and ϵENS have been calculated to evaluate the reliability of the distribution network in two cases without and with DGs of PV. The impact of the DGs of PV has been investigated and the simulation studies have shown that the implementation of DGs in the distribution system can improve the reliability of the system. In the RBTS containing DGs of PV SAIFI, SAIDI, CAIDI and EENS were improved by 8.57% 3.23% 5.84% 4.29% respectively. In the Substation 7029 containing DGs of PV SAIFI, SAIDI, CAIDI and EENS were improved by 19.5% 2.25% 21.5% 3.9%respectively.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

Cite this paper

Al-Sefri, A.K. and Al-Shaalan, A.M. (2019) Availability, Performance and Reliability Evaluation for PV Distributed Generation. World Journal of Engineering and Technology, 7, 429-454. https://doi.org/10.4236/wjet.2019.73032

References

- 1. Roy, B. and Jonnavithula, S. (1996) A Test System for Teaching Overall Power System Reliability Assessment. IEEE Transactions on Power Systems, 11, 1670-1676. https://doi.org/10.1109/59.544626
- 2. Viral, R. and Khatod, D.K. (2012) Optimal Planning of Distributed Generation Systems in Distribution System: A Review. Renewable and Sustainable Energy Reviews, 16, 5146-5165. https://doi.org/10.1016/j.rser.2012.05.020
- 3. Energy Trend PV Price Qut. http://pv.energytrend.com/pricequotes
- 4. International Energy Agency (2016) Snapshot of Global Photovoltaic Markets. 1-16.
- 5. Haegel, N.M., Margolis, R., Buonassisi, T., Feldman, D., Froitzheim, A., Garabedian, R., Green, M., Glunz, S., Henning, H.-M., Holder, B., Kaizuka, I., Kroposki, B., Matsubara, K., Niki, S., Sakurai, K., Schindler, R.A., Tumas, W., Weber, E.R., Wilson, G., Woodhouse, M. and Kurtz, S. (2017) Terawatt-Scale photovoltaics: Trajectories and Challenges. Science, 356, 141-143. https://doi.org/10.1126/science.aal1288
- 6. International Renewable Energy Agency (2016) Renewable Energy Market Analysis: The GCC Region. IRENA, Abu Dhabi.
- 7. Vision 2030 (2016) Kingdom of Saudi Arabia, April 2016.
- 8. U.S. Energy Information Administration (2017) Country Analysis Brief: Saudi Arabia. http://www.eia.gov/beta/international/analysis_includes/countries_long/ Saudi_Arabia/saudi_arabia.pdf
- 9. Wu, T.F., Shen, C.L., Nein, H.S. and Li, G.F. (2005) A Single-Phase 3 W Inverter with Grid Connection and Active Power Filtering Based on Nonlinear Programming and Fast-Zero-Phase Detection Algorithms. IEEE Transactions on Power Electronics, 20, 218-226. https://doi.org/10.1109/TPEL.2004.839786
- 10. Firdaus, M. and Daniyal, H. (2013) Modeling and Simulation of Photovoltaic Module with MPPT. 2nd Annual International Conference on Sustainable Energy and Environmental Sciences, Singapore, 25-26 February 2013.
- 11. Natsheh, E.M., Blackhurs, E.J. and Albarbar, A. (2011) PV System Monitoring and Performance of a Grid Connected PV Power Station Located in Manchester-UK. IET Conference on Renewable Power Generation, Edinburgh, 6-8 September 2011, 1-6.
- 12. Liang, H., Su, J. and Liu, S. (2010) Reliability Evaluation of Distribution System Containing Microgrid. China International Conference Electricity Distribution (CICED) 2010 Proceedings, Nanjing, 3-6 September 2010, 1-7.
- 13. Masters, G.M. (2004) Renewable and Efficient Electric Power Systems. 2nd Edition, Wiley-IEEE Press, Hoboken, NJ. https://doi.org/10.1002/0471668826
- 14. Cha, S., Jeon, D., Lee, I. and Kim, J. (2004) Reliability Evaluation of Distribution System Connected Photovoltaic Generation Considering Weather Effects. 2004 International Conference on Probabilistic Methods Applied to Power Systems, Ames, IA, 12-16 September 2004, 451-456.
- 15. SEC (2015) Electrical Data 2000-2013. Saudi Electricity Company. https://www.se.com.sa/en-us/Pages/ElectricalData.aspx
- 16. Faulin, J., Juan, A.A.,·Martorell, S. and Ramírez-Márquez, J.-E. (2010) Simulation Methods for Reliability and Availability of Complex Systems. Springer, Dordrecht, Heidelberg. https://doi.org/10.1007/978-1-84882-213-9
- 17. Billinton, R. and Allan, R.N. (1996) Reliability Evaluation of Power Systems. 2nd Edition, Plenum Press, New York. https://doi.org/10.1007/978-1-4899-1860-4
- 18. Allan, R.N., Billinton, R., Sjarief, I., Goel, L. and So, K.S. (1991) A Reliability Test System for Educational Purposes Basic Distribution System Data and Results. IEEE Transactions on Power Systems, 6, 813-820. https://doi.org/10.1109/59.76730