An Optimal Source-Load Coordinated Restoration Method Considering Double Uncertainty

: The security of power system restoration is severely affected by uncertain factors, especially the start-up time of generating unit and the amount of load pick-up. Solving the optimization restoration problem is challenging since it needs to determine different priorities in which units and loads are restored with the consideration of double uncertainty. Therefore, an optimal source-load coordinated restoration method that is based on information gap decision theory (IGDT) is proposed. Firstly, the time-domain restoration characteristics of black-start unit (BSU), non-black-start unit (NBSU), and load are described with analysis of double uncertainty. On this basis, a coupled multi-objective optimization model is built with double uncertainty, in which source-load coordinated restoration is realized. Then, IGDT is adopted to convert the uncertainty optimization model to a certain one with robustness, which tolerates the most uncertainty and still meets the desired requirement. Finally, the optimization model is solved by non-dominated genetic algorithm II (NSGA-II). The effectiveness and robustness of the proposed method is further illustrated through a case study based on the IEEE 39-bus system.


Introduction
With the continuous development of power system, the operation of the system has been more complex than ever. Due to some occasional and inevitable factors, blackouts have become a serious threat to the power system [1,2]. Besides, the security of power system has been severely affected by uncertain factors. Therefore, the investigation on source-load coordinated restoration method considering double uncertainty plays supporting role in restoration decision, which is beneficial to ensure the stability and rapidity of power system in restoration process.
The procedure of power system restoration can be divided into three stages: black-start, reconfiguration of the network and load restoration [3,4]. There is no obvious sequence as three stages are interleaved with each other. It is necessary to restore part of the important loads to balance the output of the restarted units during the black-start stage, which ensures the stable operation of the restarted units and achieves the minimum steady output as soon as possible after connecting to the grid [5,6].
In recent years, theoretical investigations on unit and load restoration have attracted researchers' interests. To maximize the generating power output, the start-up of generating units need to be optimized [7]. Sun W. et al. applied mixed-integer linear programming (MILP) algorithm to identify the optimal start-up sequence of NBSUs [8]. An optimal restoration approach based on the Wide Area Measurement System (WAMS) was proposed by Nourizadeh S. et al. [9]. Gu X. et al. proposed an extended black-start restoration strategy, in view of the fact that some current black-start resources with large rated power can provide enough cranking power for more than one non-black-start unit (1) A coupled multi-objective optimization model is built, which takes account of the double uncertainty in restoration process and realizes source-load coordinated restoration. (2) IGDT is adopted to power system restoration in order to convert the uncertainty multi-objective optimization model to a certain one with robustness. The requirement of distribution characteristics and membership function of uncertainties are avoided. Moreover, the optimal solution tolerates the most uncertainty and still meets the desired requirement.
The rest of the paper is arranged as follows: Section 2 presents the double uncertainty analysis. Section 3 presents the proposed source-load coordinated restoration model considering double uncertainty. Section 4 presents the solution to the model based on IGDT and NSGA-II. A case study is shown in Section 5. The conclusions are drawn in the last section.

Analysis of Double Uncertainty
In the process of source-load coordinated restoration, the double uncertainty including the start-up time of generating unit and the amount of load pick-up severely affect restoration time, generating Energies 2018, 11, 558 3 of 18 power output, and stability. Therefore, a study of double uncertainty is conducted in this section. The time-domain restoration characteristics of BSU, NBSU, and load are described with analysis of double uncertainty, which lays the foundation for later modeling.

Uncertainty of Start-Up Time of Generating Unit
First of all, a brief analysis of BSU i is carried out. The time-domain restoration characteristic is shown in Figure 1.
Energies 2018, 11, x FOR PEER REVIEW 3 of 17 this section. The time-domain restoration characteristics of BSU, NBSU, and load are described with analysis of double uncertainty, which lays the foundation for later modeling.

Uncertainty of Start-Up Time of Generating Unit
First of all, a brief analysis of BSU i is carried out. The time-domain restoration characteristic is shown in Figure 1. where Bi K is the ramping rate of BSU i; e t is the entire restoration time for system; and, Bci t is the time when BSU i begins to ramp up. BSU is not involved in optimization, owing to only one BSU in a system after blackouts. Therefore, suppose where max Bi P is the maximum generating power output of BSU i.
In contrast, a brief analysis of NBSU j affected by uncertainty is carried out. The time-domain restoration characteristic is shown in Figure 2. NBSU j is being cranked, but not paralleled with system from Nsj t to Ncj t . After Ncj t , NBSU j begins to ramp up. Obviously, the length of time when NBSU j is being cranked but not paralleled with system ( )  where K Bi is the ramping rate of BSU i; t e is the entire restoration time for system; and, t Bci is the time when BSU i begins to ramp up. BSU is not involved in optimization, owing to only one BSU in a system after blackouts. Therefore, suppose t Bci = 0 can reduce complexity in later modeling. The generating power output function P BSUi (t) of BSU i can be expressed as: where P maxBi is the maximum generating power output of BSU i. In contrast, a brief analysis of NBSU j affected by uncertainty is carried out. The time-domain restoration characteristic is shown in Figure 2.
Energies 2018, 11, x FOR PEER REVIEW 3 of 17 this section. The time-domain restoration characteristics of BSU, NBSU, and load are described with analysis of double uncertainty, which lays the foundation for later modeling.

Uncertainty of Start-Up Time of Generating Unit
First of all, a brief analysis of BSU i is carried out. The time-domain restoration characteristic is shown in Figure 1.
where max Bi P is the maximum generating power output of BSU i.
In contrast, a brief analysis of NBSU j affected by uncertainty is carried out. The time-domain restoration characteristic is shown in Figure 2. NBSU j is being cranked, but not paralleled with system from Nsj t to Ncj t . After Ncj t , NBSU j begins to ramp up. Obviously, the length of time when NBSU j is being cranked but not paralleled with system ( )  NBSU j is being cranked, but not paralleled with system from t Nsj to t Ncj . After t Ncj , NBSU j begins to ramp up. Obviously, the length of time when NBSU j is being cranked but not paralleled with system (t Nsj ∼ t Ncj ) is affected by the temperature and pressure of the thermal system, which is positively related to the shutting down time of the NBSU j (0 ∼ t Nsj ). Thus, the start-up time of NBSU j is uncertain and t Ncj is used to describe the uncertainty time when NBSU j begins to ramp up. K Bi is ramping rate of NBSU j; P maxNj is the maximum generating power output of NBSU j; P Nstj is the cranking power requirement of NBSU j. The generating power output function P NBSUj (t) of NBSU j can be expressed as: Although the distribution characteristics and membership function of t Ncj are hard to obtain, the fluctuation range of the time from t Nsj to t Ncj is available by dispatcher. Therefore, t Ncj can be expressed as: where t suj is also uncertain, which represents the length of time when NBSU j is being cranked. t Esuj is predicted length of that time based on experience and historical data, and a 1 is the uncertainty radius of t suj . In other words, values of t suj whose deviation from t Esuj is nowhere greater than a 1 t Esuj .

Uncertainty of Amount of Load Pick-Up
The time-domain restoration characteristic about feeder k of load l is analyzed and shown in Figure 3.
Although the distribution characteristics and membership function of Ncj t are hard to obtain, the fluctuation range of the time from Nsj t to Ncj t is available by dispatcher. Therefore, Ncj t can be expressed as: where suj t is also uncertain, which represents the length of time when NBSU j is being cranked.
Esuj t is predicted length of that time based on experience and historical data, and 1 a is the uncertainty radius of suj t . In other words, values of suj t whose deviation from Esuj t is nowhere greater than 1 Esuj a t .

Uncertainty of Amount of Load Pick-Up
The time-domain restoration characteristic about feeder k of load l is analyzed and shown in Figure 3.  The time-domain restoration characteristic of load is simplified as a step function. Owing to the influence of precipitation, temperature, cold pick-up, and discrete access, lk P is uncertain, which represents the amount of load pick-up in feeder k of load l. Therefore, lk P is used to describe the uncertainty amount of load pick-up in feeder k of load l. Lclk t is the time when feeder k of load l is restored. The function of load pick-up ( ) Loadlk t P can be expressed as: Similarly, in order to accurately express the uncertainty amount of load pick-up for later modeling, lk P can be expressed as:   The time-domain restoration characteristic of load is simplified as a step function. Owing to the influence of precipitation, temperature, cold pick-up, and discrete access, P lk is uncertain, which represents the amount of load pick-up in feeder k of load l. Therefore, P lk is used to describe the uncertainty amount of load pick-up in feeder k of load l. t Lclk is the time when feeder k of load l is restored. The function of load pick-up P Loadlk (t) can be expressed as: Similarly, in order to accurately express the uncertainty amount of load pick-up for later modeling, P lk can be expressed as: where P Elk is predicted amount of load pick-up in feeder k of load l; a 2 is the uncertainty radius of P lk .

Source-Load Coordinated Restoration Model Considering Double Uncertainty
The main task of source-load coordinated restoration is to optimize the restoration sequence of NBSUs and loads. The interval between two units begin to restart is set as a time-step, in which important loads are coordinately restored to absorb excessive active power that is generated. The source-load coordinated restoration model is formulated as a multi-objective optimization problem; the optimization model, including objective function and constraints, are described in detail in this section.

Objective of Unit Restoration
The maximum generating power output of restarted units can be expressed as: where n is the number of BSUs in the system; m is the number of NBSUs in the system. It is not easy to solve the objective function containing the integral process. Therefore, the calculation of objective function is converted to that of the maximum area that is encircled by time-domain characteristics curve of units and X axis in this paper. The computational efficiency is obviously improved, which only needs simple algebraic operations. As mentioned above, the time when BSUs begin to ramp is the same as that when the system starts to restore. The first part in Equation (6) can be ignored, as it represents the generating output of BSUs, which is constant.
The second part in Equation (6) can be expressed as the area of white trapezoid in Figure 2. The following equation represents the generating power output of NBSU j.
where t Nsj is a decision variable to be solved. The maximum generating output of all the NBSUs in system can be expressed as:

Objective of Load Restoration
While the NBSUs are being restored, the important load should be taken into account to balance active power and accelerate restoration process. The maximum weighted amount of load pick-up can be expressed as: (9) where N k is the number of load nodes in the system; N p is the number of feeders belonging to load l. ω lk is the weight of load in feeder k of load l. The weight is determined by Analytic Hierarchy Process with the type of load by Wang D. et al. [18], or even set according to the need and preference of the dispatcher.
In summary, two objectives are coupled with each other through time variable t. The objective function of the model is as follow: Energies 2018, 11, 558 6 of 18

Cranking Power Constraint
The sum of cranking power for NBSUs to be restarted and the active power of loads pick-up must be less than the available generating power output provided by restarted units. Therefore, the cranking power constraint to NBSU j is as follows: There are four parts in Equation (11). The first and second parts represent generating power output provided by restored BSUs and NBSUs, respectively. The third part is the amount of load pick-up and the last part is the cranking power for NBSU j to be restarted. u j (t) indicates the state of NBSU j at time t. Only when u j (t) = 0 and u j (t + ∆t) = 1, NBSU j begins to ramp up. NBSU should keep operating state once restarted. Assuming that NBSU j no longer shuts down after restarted, it can be expressed as follow:

Maximum Power Constraint to Load Pick-Up at One Time
The maximum power of load to be restored at one time must meet the requirements of the system stability; otherwise, the excessive voltage drop and frequency deviation may lead to failure of system restoration. The maximum power constraint to load pick-up at one time is given by: where P BSUni and P NBSUnj are the active rated power of BSU i and NBSU j, respectively; ∆ f max is the permitted frequency drop; f di and f dj are the frequency response rate of BSU i and NBSU j, which can be obtained from reference [24].

Reactive Power Constraint
The reactive power generated by restoration paths must be less than the possible maximum reactive power that is absorbed by BSUs; otherwise, the excessive reactive power may lead to sustained overvoltage and transformer over excitation. The reactive power constraint is given, as follows: where N p is the number of transmission lines in the restoration paths; Q p is the charging reactive power of the transmission line p, Q Bimax is the maximum reactive power absorbed by BSU i. Besides, the charging reactive power of restoration paths may also lead to the self-excitation of BSU. In order to avoid this situation, the self-excitation constraint to BSU is expressed as: where K Qi and S Bi are the short circuit ratio and the rated power of BSU i, respectively.

System Operating Constraint
The system operating constraint includes voltage and reactive power output of units.
where Q Gi and Q Gj are the reactive power output of BSU i and NBSU j, respectively. Q Gimin and Q Gjmin are the minimum reactive power output of BSU i and NBSU j, respectively. Q Gimax and Q Gjmax are the maximum reactive power output of BSU i and NBSU j, respectively. V i , V j and V l are voltage amplitude of BSU i, NBSU j and load l, respectively. V imax , V jmax and V lmax are the permitted voltage upper limits of BSU i, NBSU j and load l, respectively. V imin , V jmin and V lmin are the permitted voltage lower limits of BSU i, NBSU j and load l, respectively.

Maximum Critical Hot-Start Time Constraint
To make the unit restarted within the maximum critical time, the maximum critical hot-start limitation of NBSU must be longer than the time when the NBSU begins to ramp up. Therefore, the maximum critical hot-start constraint to NBSU is given by: where t CH is the maximum critical hot-start limitation of NBSU.

The Solution to Source-Load Coordinated Restoration Model Considering Double Uncertainty
IGDT is adopted to convert uncertainty multi-objective optimization model to a certain one, which tolerates the most uncertainty and still meets the desired requirement. When compared with the fuzzy random method, IGDT neither needs distribution characteristics nor membership function of uncertainties. Finally, the specific restoration scheme can be solved by NSGA-II.

Information Gap Decision Theory
IGDT can build a robust model based on different requirements owing to the different preferences of decision-maker. The robustness is defined as the immunity of the minimum requirement satisfaction at presence of uncertain parameters [22].
For a given uncertainty optimization model: where X is the uncertainty set describing the uncertain input parameters. d is the decision variables set. B(X, d) represents the optimization objective. H(X, d) = 0 and G(X, d) ≤ 0 are quality and inequality constrains, respectively. The uncertainty is defined as the distance between what is known and what may happen in reality [22]. Therefore, the envelope bound model is used to represent the prior information about the uncertain parameters X, as follows [22]: where X E is the predicted value of the uncertain parameters X; a is the uncertainty radius of X; U(a, X E ) is the set of all values of X whose deviation from X E is nowhere greater than aX. The uncertainty optimization model of Equation (18) can be converted to the certainty model based on the predicted value of the uncertain parameters X as follows: By solving Equation (20), the basic value of objective function is B 0 . However, it is obvious that the optimal solution of Equation (20) is not reliable while the value of X may be different from the predicted value X E . In order to guarantee the optimization effect, decision maker sets a minimum requirement B as not surpassing a predefined limit.
where δ is a deviation factor, which is the degree of deviation between minimum requirement and the optimal solution of the certainty model. The more δ is, the bigger value of uncertainty radius of X can be tolerated. Based on these, the uncertainty optimization model of Equation (18) can be finally converted to solve the maximum value of uncertainty radius of X when the minimum requirement is always satisfied. The robust model can be expressed as: The decision scheme of Equation (22) can tolerate the maximum value of uncertainty radius of X and meet the minimum requirement. Owing to the uncertainty set X in the model has been replaced by predicted value X E , the original objective function and constraints need to be modified to correspond with the actual requirements.

Source-Load Coordinated Restoration Model Based on IGDT
The uncertainty optimization model and envelope bound model of double uncertainty have been established above. Therefore, the maximum generating power output B 1 and the weighted amount of load pick-up B 2 are calculated by Equations (8) and (9), based on the predicted value t Esuj and P Elk . Dispatcher set B c and B d as the minimum requirements of generating output and the amount of load pick-up, respectively.
According to the idea of IGDT, two objective functions are converted to constraints in a new robust model, as follows: Energies 2018, 11, 558 9 of 18 The generating power output is minimal when the output of each NBSU is the minimum. Therefore, Equation (24) should be modified as: The weighted amount of load pick-up is minimal when the amount of load restored in each feeder is the minimum. Therefore, Equation (25) should be modified as: Besides, the constraints also need to be modified. Similarly, Equation (11) should be modified as the cranking power constraint to NBSU j when the output of each NBSU is the minimum and the amount of load that is restored in each feeder is the maximum.
Equation (13) should be modified as maximum power constraint to load pick-up at one time when the amount of load restored in each feeder is the maximum.
Equations (12) and (14)- (17) do not need further modification, owing to no relationship with double uncertainty. In summary, source-load coordinated restoration model based on IGDT can be expressed as: (12) and (14) to Equations (17) and (26) to Equation (29) Source-load coordinated restoration scheme calculated by this model is robust to double uncertainty. More specifically, the restoration scheme can guarantee that the generating power output and the amount of load pick-up are no less than (1 − δ 1 )B 1 and (1 − δ 2 )B 2 when the length of cranking time fluctuates within (1 − a 1 , 1 + a 1 ) and the amount of load pick-up fluctuates within (1 − a 2 , 1 + a 2 ).

The Solution to Source-Load Coordinated Restoration Model Based on IGDT
NSGA-II which is one of the most effective and efficient algorithms for solving multi-objective problems is applied to solve the certain model with robustness [25]. Its main characteristics are the following: (1) Two fitness functions, given by a 1 and a 2 (2) Integer codification: each individual is described by the number of s(s = m + l × k) variables that are describing the time when NBSUs start to be cranked and loads are restored. Although time is not continuous and integer, Haffner S built conversion that can transform the codified variables back to the original variables [26], as follows: Based on NSGA-II algorithm, the specific solving progress shown in Figure 4 is as follows: (1) Input basic system data.
(2) Set up the deviation factors δ 1 and δ 2 , basic parameters of NSGA-II including number of individuals in the population n pop , maximum number of generation g max , crossover probability p c , mutation probability p m . (3) Generate 3n pop individuals randomly and among them select n pop with distinct characteristics and constraints satisfaction as parent population P t . (4) Create offspring population Q t from P t by using the tournament selection, crossover and mutation operators.  Based on NSGA-II algorithm, the specific solving progress shown in Figure 4 is as follows: (1) Input basic system data. (4) Create offspring population t Q from t P by using the tournament selection, crossover and mutation operators.

Case Study
The IEEE 39-bus shown in Figure 5 is used as a case study to verify the effectiveness of the proposed method. The optimization program was developed with the programming tool of MATLAB. Suppose that the unit at bus 33 is used as the BSU, with rated power 600 MW, cos ϕ = 0.9, Q Bmax = 0.3S N , K Q = 0.45, and it is successfully restarted at the instant 0 h as the beginning time of the system restoration. The rest buses are NBSUs and loads to be restored. All of the units are subjected to a maximum critical hot start time of 1.5 h. The numbers and parameters of the units are assumed as that in Table 1. The parameters of the loads are assumed as that in Table 2. The feeder of load is represented as X(y), in which X is the node number and y is the feeder number. The power factors for all the loads are set to 0.8.

Case Study
The IEEE 39-bus shown in Figure 5 is used as a case study to verify the effectiveness of the proposed method. The optimization program was developed with the programming tool of MATLAB. Suppose that the unit at bus 33 is used as the BSU, with rated power 600 MW, cos 0.9 , and it is successfully restarted at the instant 0 h as the beginning time of the system restoration. The rest buses are NBSUs and loads to be restored. All of the units are subjected to a maximum critical hot start time of 1.5 h. The numbers and parameters of the units are assumed as that in Table 1. The parameters of the loads are assumed as that in Table 2.
The feeder of load is represented as X(y), in which X is the node number and y is the feeder number. The power factors for all the loads are set to 0.8.    The parameters of NSGA-II are set as follows: the number of individuals in the population is 200, maximum number of generation is 20, crossover probability is 0.8 and the mutation probability is 0.1. Besides, ∆x i = 1 min. The entire restoration time for system is 2 h. The charging time for each transmission line is 2 min.

Result of Source-Load Coordinated Restoation
Firstly, the maximum generating power output B 1 and the weighted amount of load pick-up B 2 are calculated by Equations (8) and (9) based on the predicted value t Esuj and P Elk . The unique solution can be obtained, in which the generating power output is 797.25 MW, the weighted amount of load pick-up is 343.23. The units and loads coordinately restored from 0 min to 85 min. After 85 min, only loads are restored further. The source-load coordinated restoration scheme without considering double uncertainty is shown in Table 3. The interval between two units begin to restart is set as a time-step. Then, double uncertainty is considered by setting the deviation factors δ 1 and δ 2 to determine minimum requirements B c and B d . For example, δ 1 is assumed as 0.03 and δ 2 is assumed as 0.2. Through solving the multi-objective optimization model that is built in this paper, the Pareto-optimal front can be obtained, as shown in Figure 6. The bigger values of uncertainty radius 1 a and 2 a are, the more robust restoration process is.
In order to achieve the balance between 1 a and 2 a , the midpoint of Pareto-optimal (0.24, 0.13) is selected, the corresponding source-load coordinated restoration scheme is shown in Table 4. The amount of load pick-up is 257.27 MW to balance the active power in system, accounting for eight percent of the generating output when the system is in normal operation.

Analysis of Optimization Characteristics
In order to further investigate the optimization characteristics of this method, the load restoration and unit restoration are optimized with the consideration of uncertainty separately. The generating power output and the weighted amount of load pick-up are compared under different deviation factors.
Firstly, the optimization characteristic of unit restoration is investigated. 2 δ is assumed as 0, 1 δ is changed to determine the minimal requirement of generating power output. The midpoint of Pareto-optimal is selected to analyze. The values of uncertainty radius 1 a and corresponding unit restoration scheme are shown in Table 5. It can be seen that the generating power output is decreasing with the increase of the deviation factor 1 δ . The bigger values of uncertainty radius a 1 and a 2 are, the more robust restoration process is. In order to achieve the balance between a 1 and a 2 , the midpoint of Pareto-optimal (0.24, 0.13) is selected, the corresponding source-load coordinated restoration scheme is shown in Table 4. The amount of load pick-up is 257.27 MW to balance the active power in system, accounting for eight percent of the generating output when the system is in normal operation.

Analysis of Optimization Characteristics
In order to further investigate the optimization characteristics of this method, the load restoration and unit restoration are optimized with the consideration of uncertainty separately. The generating power output and the weighted amount of load pick-up are compared under different deviation factors.
Firstly, the optimization characteristic of unit restoration is investigated. δ 2 is assumed as 0, δ 1 is changed to determine the minimal requirement of generating power output. The midpoint of Pareto-optimal is selected to analyze. The values of uncertainty radius a 1 and corresponding unit restoration scheme are shown in Table 5. It can be seen that the generating power output is decreasing with the increase of the deviation factor δ 1 . Then, the optimization characteristic of load restoration is investigated similarly. δ 1 is assumed as 0, δ 2 is changed to determine the minimal requirement of weighted amount of load pick-up. The first time-step is taken as an example. The values of uncertainty radius a 2 and corresponding load restoration scheme are shown in Table 6. It can be seen that the load restoration schemes can be obtained through the optimization model that is proposed under different uncertainty radius a 2 . The dispatcher can look up the table to select the appropriate load restoration scheme according to the approximate fluctuation range of load pick-up based on historical data. A comparison between load restoration scheme considering double uncertainty and that considering ideal situation is shown in Figure 7.  Then, the optimization characteristic of load restoration is investigated similarly. 1 δ is assumed as 0, 2 δ is changed to determine the minimal requirement of weighted amount of load pick-up. The first time-step is taken as an example. The values of uncertainty radius 2 a and corresponding load restoration scheme are shown in Table 6. Table 6. The load restoration scheme in first time-step considering double uncertainty.  It can be seen from Figure 7 that except time-step 6 and 8, the weighted amount of load pick-up considering double uncertainty is lower than the ideal value not only in a single time-step, It can be seen from Figure 7 that except time-step 6 and 8, the weighted amount of load pick-up considering double uncertainty is lower than the ideal value not only in a single time-step, but also in the entire restoration process. The essence of IGDT is to sacrifice part of the load pick-up and enhance the robustness of the system restoration.
Finally, the relationship between uncertainty radius and deviation factors is analyzed. δ 1 − α 1 and δ 2 − α 2 are shown in Figure 8. but also in the entire restoration process. The essence of IGDT is to sacrifice part of the load pick-up and enhance the robustness of the system restoration. Finally, the relationship between uncertainty radius and deviation factors is analyzed. The deviation factor has a positive correlation with the maximum uncertainty radius. Specifically, the smaller requirement of generating output is, the longer the start-up time of a unit can be tolerated. The smaller requirement of weighted amount of load pick-up is, the greater fluctuation of load pick-up in each feeder can be tolerated.

Comparison of Optimization Characteristics
To further verify the advantages of the proposed source-load coordinated restoration method based on IGDT, the proposed method in this paper is compared to fuzzy chance constrained programming. Assuming that the start-up time of generating unit and the amount of load pick-up are fluctuated around the predicted value in (0.8, 1.2). Double uncertainty is expressed by triangular fuzzy parameters and confidence intervals are all assumed as 0.8. The weighted amount of load pick-up in each time-step based on different methods is recorded, as shown in Figure 9. If the constraints are not satisfied, then the weighted amount of load pick-up is 0.  The deviation factor has a positive correlation with the maximum uncertainty radius. Specifically, the smaller requirement of generating output is, the longer the start-up time of a unit can be tolerated. The smaller requirement of weighted amount of load pick-up is, the greater fluctuation of load pick-up in each feeder can be tolerated.

Comparison of Optimization Characteristics
To further verify the advantages of the proposed source-load coordinated restoration method based on IGDT, the proposed method in this paper is compared to fuzzy chance constrained programming. Assuming that the start-up time of generating unit and the amount of load pick-up are fluctuated around the predicted value in (0.8, 1.2). Double uncertainty is expressed by triangular fuzzy parameters and confidence intervals are all assumed as 0.8. The weighted amount of load pick-up in each time-step based on different methods is recorded, as shown in Figure 9. If the constraints are not satisfied, then the weighted amount of load pick-up is 0.
Energies 2018, 11, x FOR PEER REVIEW 15 of 17 but also in the entire restoration process. The essence of IGDT is to sacrifice part of the load pick-up and enhance the robustness of the system restoration. Finally, the relationship between uncertainty radius and deviation factors is analyzed.  The deviation factor has a positive correlation with the maximum uncertainty radius. Specifically, the smaller requirement of generating output is, the longer the start-up time of a unit can be tolerated. The smaller requirement of weighted amount of load pick-up is, the greater fluctuation of load pick-up in each feeder can be tolerated.

Comparison of Optimization Characteristics
To further verify the advantages of the proposed source-load coordinated restoration method based on IGDT, the proposed method in this paper is compared to fuzzy chance constrained programming. Assuming that the start-up time of generating unit and the amount of load pick-up are fluctuated around the predicted value in (0.8, 1.2). Double uncertainty is expressed by triangular fuzzy parameters and confidence intervals are all assumed as 0.8. The weighted amount of load pick-up in each time-step based on different methods is recorded, as shown in Figure 9. If the constraints are not satisfied, then the weighted amount of load pick-up is 0.  The source-load coordinated restoration scheme based on fuzzy chance constrained programming is shown in Table 7. It can be seen that the source-load coordinated restoration scheme based on fuzzy chance constrained programming can not meet the security constraints in time-step 4. The descending value of frequency is more than 0.5 Hz, which exceeds the limit.

Time-step
In summary, the source-load coordinated restoration method based on IGDT can tolerate the fluctuation of the start-up time of the generating unit and the amount of load pick-up. Although the weighted amount of load pick-up is smaller than fuzzy chance constrained programming, it is more robust in the overall restoration process.

Conclusions
An optimal source-load coordinated restoration model considering double uncertainty is proposed, in which unit and load are restored coordinately in the manner of different priorities. Firstly, double uncertainty including the start-up time of generating units and the amount of load pick-up is analyzed to build time-domain characteristics of unit and load. On this basis, a coupled multi-objective optimization model is established with uncertainties; the units to be restarted can be determined by maximizing the generating power output, while the loads to be restored can be determined by maximizing the weighted amount of load pick-up. Then, IGDT is adopted to convert uncertainty multi-objective optimization model to a certain one with robustness. NSGA-II is applied to solve the proposed model. Finally, The test results on the IEEE 39-bus system show that unit and load can be effectively and coordinately restored, which tolerates the most uncertainty and meet the minimal requirements. The proposed method improves the robustness of restoration and provides a useful way for decision making of the restoration scheme.