A Method for Load Classification and Energy Scheduling Optimization to Improve Load Reliability

With the large amount of distributed generation in use, the structure of the distribution system is increasingly complex. Therefore, it is necessary to establish a method to improve load reliability. Based on the reliability model of distributed generation, this paper investigates the time sequential simulation of a wind/solar/storage combined power supply system under off-grid operation. After classifying the load by power supply region, the load weight coefficient is established, which modifies the reliability index of the load point and system. The modified expected energy not supplied (EENS) is adopted as the objective function, and the particle swarm optimization algorithm is used to solving the optimal energy scheduling for improving the load reliability. Finally, the load reliability is calculated with a hybrid method. Using the IEEE-RBTS Bus 6 system as an example, the correctness and validity of the proposed method are verified as an effective way to improve load reliability.


Introduction
In recent years, with the large amount of distributed generations (DGs) access, the traditional single-source radial distribution network has been transformed into a multi-source distribution system.When the system fails, a certain amount of load can be separated from the power grid to form a power supply region off-grid.In the power supply region off-grid, DGs can be used to power the load to improve reliability for users [1][2][3].
With the distribution network structure and operation mode changing, the reliability calculation method has also been changed.Many scholars have done the relevant research for the reliability evaluation of distribution networks containing DGs.In [4], a method was proposed to assess the reliability of a power distribution network with DG and a battery swapping station (BSS).The IEEE three-feeder distribution system was used to verify the effectiveness of the method.The literature [5] proposed an analytical formulation for assessing the reliability impact of energy storage supporting DG in supply restoration of isolated network areas.
The access of DGs improves the load reliability in distribution network.Literature [6] proposed a hybrid approach that combines scenario selection and enumerative analysis.The test system simulation was used to illustrate the dependence of reliability on factors.Some recommendations on the parameter settings of a protection system were provided to enhance the system reliability.In [7], the sequential Monte Carlo method was applied to complete the reliability evaluation of the RBTS Bus 4 test system.The effects on the system reliability of battery storage (BS) capacity and vehicle-to-grid (V2G) technology, driving behavior, recharging mode, and penetration of electric vehicles (EVs) were of electric vehicles (EVs) were all investigated.Literature [8] proposed a constrained stochastic framework to maximize the expected profit of a microgrid operator.According to the proposed method, the impact of consumer participation demand response (DR) on system reliability is analyzed by using the conditional value-at-risk (CVaR) method.In [9], an intelligent particle swarm optimization (PSO) based search method is proposed.The reliability of the Institute of Electrical and Electronics Engineers reliability test system (IEEE RTS) was calculated by using the PSO search method.
Nevertheless, there are few researches on how to improve the load reliability through the energy scheduling optimization.According to the importance of the load, the objective function of the energy scheduling optimization model in the power supply region is established in this paper.With the PSO algorithm, it obtains the optimal energy scheduling scheme.According to the optimal energy scheduling scheme, the load reliability in the power supply region is improved.The structure flow chart of this paper is shown in Figure 1.As shown in Figure 1, the aim of this paper is to improve the load reliability in the power supply region by the energy scheduling optimization.The step-by-step description is as follows: Section 1: The main work and contributions of this paper are introduced.Section 2: This section researches the reliability models of DGs.After considering the random output characteristics of DGs, the timing sequence model of wind/solar/storage joint-power-supply is established.
Section 3: Considering the load weight coefficient, the reliability indexes of load point and system are modified.According to the importance of the load, the objective function of the energy As shown in Figure 1, the aim of this paper is to improve the load reliability in the power supply region by the energy scheduling optimization.The step-by-step description is as follows: Section 1: The main work and contributions of this paper are introduced.Section 2: This section researches the reliability models of DGs.After considering the random output characteristics of DGs, the timing sequence model of wind/solar/storage joint-power-supply is established.Section 3: Considering the load weight coefficient, the reliability indexes of load point and system are modified.According to the importance of the load, the objective function of the energy scheduling optimization model in the power supply region is established.With the PSO algorithm, the optimal energy scheduling scheme is obtained.
Section 4: The load reliability is calculated by the hybrid method.The reliability of load outside the power supply region and point of common coupling (PCC) are calculated by the Goal Oriented (GO) method, and the reliability of load in the power supply region is calculated by Monte Carlo (MC) simulation.During MC simulation in the power supply region, the load classification and energy scheduling optimization in Section 3 is adopted to improve the load reliability.
Section 5: With the example system, this section verifies that the energy scheduling strategy can effectively improve the load reliability in the power supply region under the different DGs capacity, energy scheduling schemes, and the load classification numbers.
Section 6: The results of the calculation are summarized.The primary contributions of the paper are as follows: (1) After considering the random output characteristics of DGs, the timing sequence model of wind/solar/storage joint-power-supply is established.(2) According to the load reduction and load weight coefficient, the reliability indexes of load point and system are modified.The objective function of energy scheduling optimization model in the power supply region is established.With the PSO algorithm, the optimal energy scheduling scheme is obtained.(3) The reliability of load is calculated by the hybrid method.With the load classification and energy scheduling optimization, the load reliability is improved.

Wind Turbine Output Power Model
Weibull distribution is used to describe the probability distribution characteristics of the wind speed in this paper [10].The probability density function is as follows: where v is the wind speed, c and k are the scale and shape parameters, respectively.The power output of the wind turbine (WT) depends on the wind speed.The relationship between output P w (v) and wind speed v of the wind turbine is represented as follows: where v ci is the cut-in wind speed, v co is the cut-out wind speed, v r is the rated wind speed, and P r is the WT rated output power.

Photovoltaic Output Power Model
Solar irradiance obeys Beta distribution over a certain time period [11,12].The probability density function is as follows: where α and β are the shape coefficients of Beta distribution, Γ is the gamma function, r and r max are the actual light intensity and the maximum light intensity in the time period, respectively.
For a photovoltaic (PV) power generation module, the total area and the total photoelectric conversion efficiency of the battery assembly are A and η, respectively.Then the battery assembly output power is P p = rAη, and the output power probability density function of the PV is as follows: where R p = Aηr max is the maximum output power of the PV generator.

Energy Storage System Model
An energy storage system can be used to smooth the random output of DGs.When the grid fails, the load in the power supply region is powered by the DGs and the energy storage system [13][14][15].During the load off-grid operation, the energy storage system is charged when the total power output of the WT generation and the PV unit is greater than the total energy consumed by the load.The energy storage system is discharged when the total power output of the WT generation and the PV unit is less than total power of the load.
Considering the constraints of charging and discharging power and capacity of energy storage system, the model of energy storage system is established as follows: where P in and P out are the charge and discharge power of energy storage, P ch-max and P dch-max are the maximum charging and discharging power of energy storage system, η c and η d are the charge and discharge efficiency of energy storage, respectively, E remain (t) is the remaining energy storage capacity at time t, P PV (t), P WT (t), and P L (t) are the power of the WT generation, PV unit and load at time t, E max and E min are the maximum capacity and the minimum capacity of energy storage, respectively.

Timing Sequence Simulation of Wind/Solar/Storage Joint-Power-Supply in Off-Grid Operation
A time sequential model of wind/solar/storage is established in this paper.With the time sequential simulation method, it takes one hour as the basic step length.When the load is in off-grid operation, load reduction and load outage represent the power supply reducing and the power supply stopping to the load, respectively.It defines the power supply degree k i (t), in which k i (t) represents the ratio of the load i power consumption after and before the load reduction.When the power supply to the load is too low, the load cannot maintain the basic operation.It defines the minimum of power supply degree k pt,i (t).When the k pt,i (t) is less than k i (t), the load i will be outage.
Figure 2 shows the time sequence simulation diagram of the wind/solar/storage joint-power-supply system.The abscissa represents the duration of off-grid operation; the ordinate represents the state of the load or energy storage; "1", "0" and "−1" of the load represent the operation, outage and excessive reduce respectively; and "1" and "0" of the energy storage represent operation and outage, respectively.
As shown in Figure 2, the load can continue to operate off-grid only when the load and energy storage are in state "1".The others cases will cause a load outage.After the fault components of the distribution network are repaired, the load can be switched to grid-connected operation.As shown in Figure 2, the load can continue to operate off-grid only when the load and energy storage are in state "1".The others cases will cause a load outage.After the fault components of the distribution network are repaired, the load can be switched to grid-connected operation.

Load Classification and Energy Scheduling Optimization
In this paper, the reliability indexes of load point and system are modified.According to the modified reliability index, the energy scheduling optimization model is established.With the PSO algorithm, the energy scheduling optimization problem is solved.The MC simulation method is used to simulate the power supply region off-grid operation.According to the simulation results, the energy scheduling scheme is modified again.Finally, the optimal energy scheduling scheme is obtained.

Modification of the Load Reliability Index Considering Load Curtailment
The power grid companies usually make compensation to users who implement load reduction or load outage by signing agreements with users [16].In the t-th hour, the user compensation amount ( ) C t depends on the power supply degree ( ) k t , and the relationship between the two is shown in Figure 3.In Figure 3, c ( ) C t is the user maximum compensation amount, and 0 ( ) k t is the load threshold at the t-th hour.The calculation method can effectively ensure that the user can obtain a reasonable amount of compensation after implementing load reduction or load outage, and an upper limit is set on the compensation amount at the same time.
According to the relationship between the load reduction of users and compensation amount, the reliability indexes of the load point are modified.The commonly used load point reliability

Load Classification and Energy Scheduling Optimization
In this paper, the reliability indexes of load point and system are modified.According to the modified reliability index, the energy scheduling optimization model is established.With the PSO algorithm, the energy scheduling optimization problem is solved.The MC simulation method is used to simulate the power supply region off-grid operation.According to the simulation results, the energy scheduling scheme is modified again.Finally, the optimal energy scheduling scheme is obtained.

Modification of the Load Reliability Index Considering Load Curtailment
The power grid companies usually make compensation to users who implement load reduction or load outage by signing agreements with users [16].In the t-th hour, the user compensation amount C(t) depends on the power supply degree k(t), and the relationship between the two is shown in Figure 3.
As shown in Figure 2, the load can continue to operate off-grid only when the load and energy storage are in state "1".The others cases will cause a load outage.After the fault components of the distribution network are repaired, the load can be switched to grid-connected operation.

Load Classification and Energy Scheduling Optimization
In this paper, the reliability indexes of load point and system are modified.According to the modified reliability index, the energy scheduling optimization model is established.With the PSO algorithm, the energy scheduling optimization problem is solved.The MC simulation method is used to simulate the power supply region off-grid operation.According to the simulation results, the energy scheduling scheme is modified again.Finally, the optimal energy scheduling scheme is obtained.

Modification of the Load Reliability Index Considering Load Curtailment
The power grid companies usually make compensation to users who implement load reduction or load outage by signing agreements with users [16].In the t-th hour, the user compensation amount ( ) C t depends on the power supply degree ( ) k t , and the relationship between the two is shown in Figure 3.In Figure 3, c ( ) C t is the user maximum compensation amount, and 0 ( ) k t is the load threshold at the t-th hour.The calculation method can effectively ensure that the user can obtain a reasonable amount of compensation after implementing load reduction or load outage, and an upper limit is set on the compensation amount at the same time.
According to the relationship between the load reduction of users and compensation amount, the reliability indexes of the load point are modified.The commonly used load point reliability In Figure 3, C c (t) is the user maximum compensation amount, and k 0 (t) is the load threshold at the t-th hour.The calculation method can effectively ensure that the user can obtain a reasonable amount of compensation after implementing load reduction or load outage, and an upper limit is set on the compensation amount at the same time.
According to the relationship between the load reduction of users and compensation amount, the reliability indexes of the load point are modified.The commonly used load point reliability indexes are failure frequency λ and outage time U.Only the outage time index of the load changes with the Energies 2018, 11, 1558 6 of 19 load power supply degree.When load reduction is considered, the power outage time of load i is defined as where ∆U i (t) represents the equivalent outage time of load i in the t-th hour, ∆U i is the equivalent outage time of load i when it is in off-grid operation, T max is the time from the start of off-grid operation to the return to grid-connected operation, and a i is a constant which is set according to the user maximum compensation.The bigger the compensation amount is, the larger the a i value is.

Modification of the System Reliability Index Considering Load Classification
In this paper, different loads are set to different weight coefficients according to the importance of the load.The more important the load is, the greater the weight coefficient is [17].The load is converted into the equivalent load according to its weight coefficient.The equivalent load P EL is defined as where β i is the weight coefficient of the i-th load, and P L is the power of the i-th load.A load with a large weight coefficient occupies a large proportion in the equivalent load, so the energy distribution prioritizes the important load.
The weight coefficient increases the influence of the important load on the system reliability index.According to the concept of equivalent load, the reliability indexes of the system are modified.Taking the system average interruption frequency index SAIFI as an example, the SAIFI is defined as where λ i is the average failure rate of load i, and N i is the number of users of load i.

Objective Function of Energy Scheduling
According to the load weight, the energy scheduling can improve load reliability when the power supply region is in off-grid operation.The expected energy not supplied (EENS) is chosen as the objective function to solve the energy scheduling optimization problem in this paper.Each time the power supply region commences off-grid operation, part of the load experiences outage and produces an incremental ∆EENS.The purpose of energy scheduling optimization is to make each ∆EENS a minimum so that the EENS is the minimum, to thus improve load reliability.
When the power supply region is operated off-grid, the load is powered by the DGs and the energy storage system.The power is distributed to the load according to the power supply degree.If the power supply degree is K = {k i (t)}, the assigned power of load i is k i (t)P L,i (t).To simplify the calculation, suppose that the power supply degree for load is the same during every hour on the off-grid operation, the power supply degree for all loads is set to K = {k i }, the power allocated by load i is k i P L,i (t).According to K, the ∆EENS for each off-grid operation is calculated as follows: where n is the number of loads in the power supply region.

Energy Scheduling Constraint Condition
(1) During the off-grid operation, the energy storage system is constrained by the maximum capacity and maximum charge and discharge conditions: (13) (2) Because of the lack of DGs and energy storage output power, load outage occurs in the power supply region.The load can't be restored to operation until grid-connected operation: where t cut,i is the moment when load i is powered off during off-grid operation.
(3) When the power supply degree of load is lower than the minimum power supply degree, the energy can't supply the power required for the basic operation of the load.Therefore, the load power supply degree can't be lower than a certain value: (4) During the off-grid operation, the total power consumption of the load should be less than the total power supply of the DGs and energy storage system: where P DG (t) is the total power generation of the distributed generation at time t.

Optimization of Energy Dispatching by the Particle Swarm Optimization Algorithm
The PSO algorithm is a stochastic algorithm that can be optimized in a complex nonlinear solution space.It is concise and easy to operate without too much parameter adjustment.It has been widely used in various optimization problems [18,19].The PSO algorithm is used to solve the energy scheduling optimization problem in this paper.Equations ( 13) and ( 14) are constraint conditions under the power supply region off-grid operation.During the calculation, K is corrected until constraints (13) and ( 14) are satisfied.The PSO algorithm is as follows: Step 1: The period from the start of a power supply region's off-grid operation to the recovery of grid-connected operation is T max , the weight coefficient of the load i is β i , and the minimum power supply degree of load i is k pt,i (t).During load off-grid operation, the power of the DGs and each load are P DG (1), P DG (2), . . ., P DG (T max ) and P L,i (1), P L,i (2), . . ., P L,i (T max ) respectively.
Step 2: Set all the load off-grid operation time is T max , the total power during the off-grid operation of each load P L,i can be calculated as follows: Step 3: The relationship between the objective function and K is determined according to (8), (9), and (12), and then considering the constraints of ( 15) and ( 16), using the PSO algorithm to solve the energy scheduling optimization and obtain the optimal energy scheduling scheme K = {k i }.
(a) Initialize the number of particles in the particle group.For particle velocity V and position K, the i-th particle is where c 1 and c 2 are learn factors, u 1 and u 2 are (0, 1) on the uniformly distributed random number, respectively.(d) Judge whether the number of cycles has been reached.If yes, the global extreme value gBest is the value of ∆EENS, and the corresponding solution is power supply degree K. Otherwise, let α = α + 1, and return to (b).
Step 4: According to the optimal scheduling scheme, the simulation method is used to simulate the off-grid operation status of the load, and judge whether there is a power shortage in the energy scheduling scheme: (a) Set the off-grid operation time to t = 1, and the energy storage remaining power to E remain (t) = E max .(b) Judge whether ∑ k i P L,i − P DG (t) > 0. If the condition is met, go to the (c).Otherwise, accumulated charge power of energy storage E remain (t) = E remain (t) + η c (P DG (t) − ∑ k i P L,i (t)), go to (e).(c) Judge whether η d (E remain (t) − E min ) > ∑ k i P L,i − P DG (t).If the condition is met, go to the (d).
Otherwise, cut off the load with the least weight coefficient in the off-grid operation load, use (19) to recalculate the load total demand power P L,i during the off-grid operation, and then return to Step 3.
(d) Judge whether η d P dch-max > ∑ k i P L,i − P DG (t).If yes, accumulated the discharge power of energy storage, make E remain (t) = E remain (t) − 1 η d (∑ k i P L,i (t) − P DG (t)), go to (f).Otherwise, cut off the load with the least weight coefficient in the off-grid operation load, use (19) to recalculate the load total demand power P L,i during the off-grid operation, and then return to Step 3. (e) Judge whether E remain (t) ≥ E max .If yes, make E remain (t) = E max .(f) Judge whether t = T max .If yes, the load is changed to grid-connected operation, and this off-grid operation ends.Otherwise, the accumulated off-grid operation time is t = t + 1, and return to (b).
Energies 2018, 11, 1558 The flow chart is shown in Figure 4.
Energies 2018, 11, x FOR PEER REVIEW 9 of 19 The flow chart is shown in Figure 4.

Load Reliability Evaluation Based on the Hybrid Method
The load reliability by hybrid method is calculated in this paper.The reliability of load outside the power supply region and PCC are calculated by the GO method, and the reliability of load in the power supply region is calculated by MC simulation.
The hybrid method combines the advantages of the analytical method and MC simulation method.Firstly, the complex distribution network is equated with a simple radial main feeder network by the analytical method.Then, the MC simulation method is used to evaluate the reliability of the simplified distribution network.Compared with the traditional MC simulation method, this method significantly improves the simulation efficiency and calculation accuracy [20].

Using GO Method to Calculate Load Reliability outside the Power Supply Region
The load outside the power supply region is only powered by the power grid without DGs.Therefore, calculating the reliability of the load grid-connected operation does not require considering the influence of the DGs supply.
There are problems in modeling complexity and the difficulty in balancing the accuracy and speed of calculation when traditional analytical methods are applied to the reliability assessment of complex distribution systems.The GO method is used to calculate the reliability of the PCC and the load outside the power supply region in this paper.The method is a success-oriented system probability analysis technology to improve reliability calculation accuracy and speed [21][22][23].Compared with the traditional analytic method, the GO method can fully consider the failure rate of switch action, the planned maintenance of the component, and the occurrence of multiple failures.According to the successful correlation between system devices, the reliability index of each load point is directly calculated.

Load Reliability Evaluation Based on the Hybrid Method
The load reliability by hybrid method is calculated in this paper.The reliability of load outside the power supply region and PCC are calculated by the GO method, and the reliability of load in the power supply region is calculated by MC simulation.
The hybrid method combines the advantages of the analytical method and MC simulation method.Firstly, the complex distribution network is equated with a simple radial main feeder network by the analytical method.Then, the MC simulation method is used to evaluate the reliability of the simplified distribution network.Compared with the traditional MC simulation method, this method significantly improves the simulation efficiency and calculation accuracy [20].

Using GO Method to Calculate Load Reliability outside the Power Supply Region
The load outside the power supply region is only powered by the power grid without DGs.Therefore, calculating the reliability of the load grid-connected operation does not require considering the influence of the DGs supply.
There are problems in modeling complexity and the difficulty in balancing the accuracy and speed of calculation when traditional analytical methods are applied to the reliability assessment of complex distribution systems.The GO method is used to calculate the reliability of the PCC and the load outside the power supply region in this paper.The method is a success-oriented system probability analysis technology to improve reliability calculation accuracy and speed [21][22][23].Compared with the traditional analytic method, the GO method can fully consider the failure rate of switch action, the planned maintenance of the component, and the occurrence of multiple failures.According to the successful correlation between system devices, the reliability index of each load point is directly calculated.
When components in the distribution system, such as transformers and transmission lines, fail, the power grid can't power to the load.The frequency probability P and time probability P U of each component's equivalent operation are calculated as follows: where λ is the failure rate, N is a time period (generally, 8760 h), r is the average failure duration of a component, and U is the component average power outage time.
For the load point i, the failure of the system components can be divided into three categories: (1) a component that requires repair can make the load restore the power supply, that is, the A class of components denoted as N 1 ; (2) a component whose fault can make the load restore the power supply by switching the switch, that is, the B class of components, denoted as N 2 ; and (3) a component that has no effect on the load, that is, the C class of components, denoted as N 3 .The success probability coefficients a A and a B and the time probability coefficients b A and b B of the class A and B components relative to load i are defined as follows: Considering that the reliable operation rate of the circuit breaker is not 100%, the success probability coefficient a C and the time probability coefficient b C of the C class components are defined as follows: where P i and P U,i are the successful operation of the frequency probability and time probability, and P b is the circuit breaker reliable action probability.
According to (21) to (23), the reliability indexes of load point n are defined as follows:

Using the Monte Carlo Simulation Method to Calculate Load Reliability in the Power Supply Region
The DGs output in the power supply region has a random nature.The MC simulation method can be used to calculate multiple uncertain events through multiple random sampling.The MC simulation method is used to calculate the reliability of load in the power supply region in this paper.
After sampling the external equivalent component and the components in the power supply region with the MC method, it judges whether the power supply region off-grid operates successfully.It uses the PSO algorithm to solve the energy scheduling optimization, and calculates the load reliability in the power supply region with the obtained power supply degree K.The specific steps of the algorithm are as follows: Step 1: Set the number of each load's failure time f [i] = 0, the outage time to U[i] = 0 and the simulation time to t = 0. Set the maximum simulation time t max .
Step 2: Generate a random number uniformly distributed on (0, 1) for the external equivalent component and each of the components within the power supply region.The time to failure (TTF) of each component is calculated as follows: where u is a random number uniformly distributed between (0, 1), and λ i is the annual average failure rate of the i-th component.
Step 3: Find the component with the smallest TTF, generate a random number uniformly distributed on (0, 1) for the component, and calculate the component's time to repair (TTR): where µ i is the repair rate of the i-th component.
Step 4: If the component is the external equivalent component, it can proceed as follows: (a) Generate a random number that is uniformly distributed on (0, 1), and calculate the success rate of off-grid operation to judge whether the power supply region successfully commences off-grid operation when the component faults.(b) If yes, generate a random number that is uniformly distributed on (0, 1), determine the time at which the fault occurred in a year, and calculate the power of the DGs and load in TTR hours from that moment.(c) According to the wind/solar/storage joint-power-supply and PSO algorithm, solve the energy scheduling optimization problem of the off-grid operation.According to the obtained K, accumulate the failure times outage time of the load in the power supply region.
If the component is not the external equivalent component, the outage time of the load is TTR, and the failure times and outage time of the load in the power supply region are calculated as follows: Step 5: Calculate the simulation time: time = time + TTF i + TTR i .
Step 6: Judge whether time > t max .If yes, calculate the load point reliability index, and end the cycle.Otherwise, proceed to Step 2.

Calculating the Reliability of the Entire Network by the Hybrid Method
In this paper, the GO method is used to calculate the reliability of PCC and the load in the power supply region when the load is in grid-connected operation.MC simulation is used to sample the component and simulate the off-grid operation of the power supply region.The PSO algorithm is used to solve the energy scheduling optimization problem.Then, according to the obtained K or TTR, the failure times and outage time of each load are calculated.Finally, the load reliability of the entire network is calculated.
The evaluation process of the hybrid method is shown in Figure 5.

The Example System Parameters
The Feeder 4 RBTS BUS6 system [24] is adopted as a test system in this paper.The example electricity distribution network is shown in Figure 6.
The feeder consists of 30 lines, 23 load points, 23 fuses, 23 transformers, 4 circuit breakers, and 1 isolating switch.The reliable action rate of the circuit breaker is 80%, the reliable operation rate of the fuse is 100%, the failure rate of the line is 0.065 f/km-year, the repair time of each line is 5 h, the transformer failure rate is 0.015 f/km-year, the repair time is 200 h, and the operation time of the isolating switch is 1 h.
In the power supply region 1 and the power supply region 2, there is a DG composed of a wind/solar/storage joint-power-supply system.The parameters of the two DGs are the same.The

The Example System Parameters
The Feeder 4 RBTS BUS6 system [24] is adopted as a test system in this paper.The example electricity distribution network is shown in Figure 6.
The feeder consists of 30 lines, 23 load points, 23 fuses, 23 transformers, 4 circuit breakers, and 1 isolating switch.The reliable action rate of the circuit breaker is 80%, the reliable operation rate of the fuse is 100%, the failure rate of the line is 0.065 f/km-year, the repair time of each line is 5 h, the transformer failure rate is 0.015 f/km-year, the repair time is 200 h, and the operation time of the isolating switch is 1 h.
In the power supply region 1 and the power supply region 2, there is a DG composed of a wind/solar/storage joint-power-supply system.The parameters of the two DGs are the same.The WT rated power is 1.5 MW.The cut-in, cut-out, and rated wind speed are 2.5, 22 and 10.5 m/s, respectively.The maximum output power of PV is 1.5 MW; the α and β are 0.230 and 0.915, respectively.The storage capacity is 2 MW, energy storage charge and discharge efficiency η c and η d are 80%.The weight coefficient of the load outside the power supply region is 1.0.The minimum power supply degree of the load in the power supply region is 0.8.

The Influence of DG Capacity on System Reliability Index
This paper uses reliability index EENS to study the influence of DGs capacity on system reliability index.
The capacity of the energy storage system is unchanged at 2 MW.The load classification, load reduction, and energy scheduling are not considered.The EENS changes with the capacity of the WT and the PV as shown in Figure 7.The three groups in Figure 7 are WT, PV, and WT + PV, where the capacity of the WT and PV is the value corresponding to the abscissa.In the group of WT + PV, the capacity of the WT is the same as the capacity of PV, whose capacity is the value corresponding to the abscissa.

The Influence of DG Capacity on System Reliability Index
This paper uses reliability index EENS to study the influence of DGs capacity on system reliability index.
The capacity of the energy storage system is unchanged at 2 MW.The load classification, load reduction, and energy scheduling are not considered.The EENS changes with the capacity of the WT and the PV as shown in Figure 7.The three groups in Figure 7 are WT, PV, and WT + PV, where the capacity of the WT and PV is the value corresponding to the abscissa.In the group of WT + PV, the capacity of the WT is the same as the capacity of PV, whose capacity is the value corresponding to the abscissa.From the Figure 7, it can be seen that the addition of DGs can improve system reliability.In the initial stage, DGs have more obvious effect on reliability improvement.When DGs capacity increases continuously, the effect of system reliability improvement gradually saturates.In this example, the effect of WT on reliability improvement is better than that of PV.
The EENS of the system, calculated by taking several sets of different WT + PV capacity combinations from Figure 7, can be shown in Figure 8. From Figure 8, it can be seen that the EENS decline in the initial stage of the energy storage joining system is obvious.As the capacity of energy storage systems continues to increase, the rate of decline in EENS tend to be flat.The impact of EENS on the WT+PV capacity is also decreasing.When the capacity of the energy storage increases, the EENS decreases stably.It can reflect the characteristics in that the energy storage can smooth the random output of the DG to improve the power supply quality.

Reliability Index of a Different Energy Scheduling Scheme
According to the hybrid method evaluation process, the reliability evaluation of the distribution network with DG is performed in three schemes: Scheme 1: no connecting DGs in the distribution network.Scheme 2: connecting DGs in the distribution network, no weight coefficient is set for the load, and all load weight coefficient are the same.From the Figure 7, it can be seen that the addition of DGs can improve system reliability.In the initial stage, DGs have more obvious effect on reliability improvement.When DGs capacity increases continuously, the effect of system reliability improvement gradually saturates.In this example, the effect of WT on reliability improvement is better than that of PV.
The EENS of the system, calculated by taking several sets of different WT + PV capacity combinations from Figure 7, can be shown in Figure 8. From the Figure 7, it can be seen that the addition of DGs can improve system reliability.In the initial stage, DGs have more obvious effect on reliability improvement.When DGs capacity increases continuously, the effect of system reliability improvement gradually saturates.In this example, the effect of WT on reliability improvement is better than that of PV.
The EENS of the system, calculated by taking several sets of different WT + PV capacity combinations from Figure 7, can be shown in Figure 8. From Figure 8, it can be seen that the EENS decline in the initial stage of the energy storage joining system is obvious.As the capacity of energy storage systems continues to increase, the rate of decline in EENS tend to be flat.The impact of EENS on the WT+PV capacity is also decreasing.When the capacity of the energy storage increases, the EENS decreases stably.It can reflect the characteristics in that the energy storage can smooth the random output of the DG to improve the power supply quality.

Reliability Index of a Different Energy Scheduling Scheme
According to the hybrid method evaluation process, the reliability evaluation of the distribution network with DG is performed in three schemes: Scheme 1: no connecting DGs in the distribution network.Scheme 2: connecting DGs in the distribution network, no weight coefficient is set for the load, and all load weight coefficient are the same.From Figure 8, it can be seen that the EENS decline in the initial stage of the energy storage joining system is obvious.As the capacity of energy storage systems continues to increase, the rate of decline in EENS tend to be flat.The impact of EENS on the WT+PV capacity is also decreasing.When the capacity of the energy storage increases, the EENS decreases stably.It can reflect the characteristics in that the energy storage can smooth the random output of the DG to improve the power supply quality.

Reliability Index of a Different Energy Scheduling Scheme
According to the hybrid method evaluation process, the reliability evaluation of the distribution network with DG is performed in three schemes: Scheme 1: no connecting DGs in the distribution network.Scheme 2: connecting DGs in the distribution network, no weight coefficient is set for the load, and all load weight coefficient are the same.Scheme 3: connecting DGs in the distribution network and the weights for the load are set.In the power supply region 1, the weights for loads 14, 15, 16, 17, and 18 are 1.0, 1.05, 1.1, 1.2, and 1.3 respectively.In the power supply region 2 the weights for loads 19, 20, 21, 22, and 23 are 1.0, 1.05, 1.1, 1.2, and 1.3 respectively.
Through the described model and algorithm, the reliability index of each load point and the entire system are calculated.Figure 9 and Table 1 list the reliability indexes of the load points and the reliability indexes of the system, respectively.
As shown in Figure 9, loads 1 to 13 are load points outside the power supply region, and their reliability is not affected by DGs.
Through the described model and algorithm, the reliability index of each load point and the entire system are calculated.Figure 9 and Table 1 list the reliability indexes of the load points and the reliability indexes of the system, respectively.
As shown in Figure 9, loads 1 to 13 are load points outside the power supply region, and their reliability is not affected by DGs.Based on the results of 1 and scheme 2 in Figure 9 and Table 1, it can be observed that after accessing the DGs, the failure rate and the annual average outage time of the load point in the power supply region are reduced.This outcome indicates that the access of the DGs can effectively improve load reliability.Compared with the reliability of the load points in scheme 2 and scheme 3, more power can be allocated for the load with a larger weight coefficient and their reliability is improved.Less power will be allocated and the reliability will be reduced for the load with a small weight coefficient.Compared with the traditional algorithm and modified algorithm in Scheme 3, the modified reliability indexes can improve the effect of the large weight coefficient load on the system reliability, which improves system reliability and verifies the effectiveness of the proposed method.

Comparison of Load Reliability Indexes with Different Classification Numbers
Taking the power supply region 2 as an example, the load reliability indexes with different classification numbers are studied.Under the condition of the same equivalent load, the weight coefficient of the load is changed, and energy scheduling optimization is adopted to compare load reliability with different classification numbers.The load reliability index in the power supply region is shown in Table 2.As shown in Table 2, the difference in load reliability with the same weight coefficient is small, and the difference in load reliability with a different weight coefficient is large.The larger the load weight coefficient is, the more power is allocated in the energy scheduling and the higher the reliability.Therefore, we can set a large weight coefficient for an important load to meet a demand for priority power supply for that load.

Comparison Load Reliability Indexes with Different Capacity Ratios
Taking the power supply region 2 as an example, the load reliability indexes with different capacity ratios are studied.In the case of energy scheduling optimization, the influence of different capacity ratios in different load classes on load reliability is discussed under the condition of the same load classification number as follows: Scheme 1: the load in the power supply region is divided into two types according to degree of importance: loads 19, 20, 21: weight coefficient 1.0, accounting for 58% of the total load; loads 22 and 23: weight coefficient 1.4, accounting for 42% of the total load; Scheme 2: the load in the power supply region is divided into three types according to degree of importance: loads 19 and 20: weight coefficient 1.0, accounting for 32% of the total load; loads 21, 22, and 23: weight coefficient 1.4, accounting for 68% of the total load; Scheme 3: the load in the power supply region is divided into four types according to degree of importance: load 19: weight coefficient 1.0, accounting for 14% of the total load; loads 20, 21, 22, and 23: weight coefficient 1.40, accounting for 86% of the total load.
The load reliability indexes in the power supply region are shown in Table 3.As shown in Table 3, in scheme 1, the load ratio of 1.2 is 42%, the failure rate of load 22 is 1.611 f/year, and the annual mean outage time is 8.692 h/year.In scheme 3, the load ratio of 1.2 is 86%, the failure rate of load 22 is 2.210 f/year, and the annual mean outage time is 11.979 h/year.It indicates that load 22 in scheme 1 has a higher reliability.The total power of the important load increases, and the average power allocation of each load is reduced.Then the reliability is reduced.For a load with a weight coefficient 1.0, such as load 19 in scheme 1, the failure rate is 2.831 f/year and the annual mean outage time is 16.357 h/year.In scheme 3, the load 19 failure rate and annual mean outage time are 2.950 f/year and 17.540 h/year, respectively.The output power of DGs and energy storage are mostly allocated to important loads, and the amount of loads assigned to a small weight coefficient is small.In order to improve the reliability of important loads, a reasonable weight coefficient should not only be set but also the proportion of various types of load should be considered.

Conclusions
This paper establishes a reliability model of DGs and an energy storage joint-power-supply based on modified load reliability indexes according to load classification, which have improved the reliability indexes of the load and system.The PSO algorithm is used to solve the optimal energy scheduling problem, and load reliability indexes under different load scheduling schemes, classification numbers, and capacity ratios are compared.The research shows the following:

Figure 1 .
Figure 1.The structure flow chart of this paper.

Figure 1 .
Figure 1.The structure flow chart of this paper.

Figure 3 .
Figure 3. Relationship between compensation amount ( ) C t and power supply degree ( ) k t .

Figure 3 .
Figure 3. Relationship between compensation amount ( ) C t and power supply degree ( ) k t .

Figure 3 .
Figure 3. Relationship between compensation amount C(t) and power supply degree k(t).

Figure 4 .
Figure 4. PSO algorithm used to solve the energy scheduling problem.

Figure 4 .
Figure 4. PSO algorithm used to solve the energy scheduling problem.

Energies 2018 ,
11, x FOR PEER REVIEW 13 of 19 WT rated power is 1.5 MW.The cut-in, cut-out, and rated wind speed are 2.5, 22 and 10.5 m/s, respectively.The maximum output power of PV is 1.5 MW; the α and β are 0.230 and 0.915, respectively.The storage capacity is 2 MW, energy storage charge and discharge efficiency c η and d η are 80%.The weight coefficient of the load outside the power supply region is 1.0.The minimum power supply degree of the load in the power supply region is 0.8.

Figure 7 .
Figure 7. EENS with different WT and PV capacities.

Figure 7 .
Figure 7. EENS with different WT and PV capacities.

Figure 9 .
Figure 9. Failure rate and annual average interruption duration of load points.(a) Failure rate of load points; (b) Annual mean outage time of load points.

Figure 9 .
Figure 9. Failure rate and annual average interruption duration of load points.(a) Failure rate of load points; (b) Annual mean outage time of load points.

Scheme 1 :
the load in the power supply region is divided into two types according to degree of importance as follows: loads 19, 20, 21: weight coefficient 1.0; loads 22 and 23: weight coefficient 1.488; Scheme 2: the load in the power supply region is divided into three types according to degree of importance as follows: loads 19 and 20: weight coefficient 1.0; load 21: weight coefficient 1.20; loads 22 and 23: weight coefficient 1.36; Scheme 3: the load in the power supply region is divided into four types according to degree of importance as follows: loads 19 and 20: weight coefficient 1.0; load 21: weight coefficient 1.20; load 22: weight coefficient 1.30; load 23: weight coefficient 1.40.
initialize the iteration number as α = 1, and set the maximum number of iterations α max .(b)According to the DGs and energy storage joint-power-supply model, calculate each particle's adaptive value Fit[i], that is, ∆EENS, using the previously described time sequence simulation.Then, correct the particle's position K i and compare it with the individual extremum pBest[i].If Fit[i] < pBest[i], use Fit[i] to replace pBest[i].Then, compare Fit[i] with the global extremum gBest.If Fit[i] < gBest, use Fit[i] to replace gBest.(c) Update the velocity and position of the particle: ). Initialize the individual extremum pBest[i] and the global