1

Load Forecasting Based Distribution System Network Reconfiguration
–A Distributed Data-Driven Approach
Yi Gu1 , Student Member, IEEE, Huaiguang Jiang2, Member, IEEE, Jun Jason Zhang1 , Senior Member,
IEEE, Yingchen Zhang2 , Senior Member, IEEE, Eduard Muljadi2, Fellow, IEEE, and Francisco J. Solis3 ,
1

arXiv:1711.10690v1 [cs.SY] 29 Nov 2017

Dept. of Electrical and Computer Engineering, University of Denver, Denver, CO, 80210
2
National Renewable Energy Laboratory, Golden, CO 80401
3
School of Mathematical and Natural Sciences, Arizona State University, Glendale, AZ, 85306

Abstract—In this paper, a short-term load forecasting approach
based network reconfiguration is proposed in a parallel manner.
Specifically, a support vector regression (SVR) based short-term
load forecasting approach is designed to provide an accurate load
prediction and benefit the network reconfiguration. Because of
the nonconvexity of the three-phase balanced optimal power flow,
a second-order cone program (SOCP) based approach is used
to relax the optimal power flow problem. Then, the alternating
direction method of multipliers (ADMM) is used to compute
the optimal power flow in distributed manner. Considering the
limited number of the switches and the increasing computation
capability, the proposed network reconfiguration is solved in a
parallel way. The numerical results demonstrate the feasible and
effectiveness of the proposed approach.

Index terms— Electrical distribution system, network reconfiguration, alternating direction method of multipliers, support vector regression, semidefinite relaxation programming,
convex optimization, optimal power flow, short-term load
forecasting
N OMENCLATURE
G
ω

b
C
γ
ε
Rrisk
ϕ1 , ϕ2
θ1 , θ2
ηgΩ
αΩ
i3
νΩ
i3
ηΩ
i3
Vi

The graph for PDS with a node set and a
link set: G = [V, E].
ω is a linear combination coefficient of
SVR, 21 ω T ω indicates the flatness of the
regression coefficients.
b is an offset coefficient of SVR.
C is a trade-off parameter of SVR.
γ is a parameter of (Gaussian) radial basis
function.
ε is an adjustable precision parameter,
which indicates the training error threshold.
The objective function of designed for the
SVR.
Two positive acceleration coefficients for
velocity updates of each particle in the PSO.
Two independently random variables with
uniformly distributed range (0, 1) for velocity updates of each particle in the PSO.
The best position vector among all particles
in PSO.
The vector of position of the i3 th particle
in PSO.
The vector of velocity of the i3 th particle in
PSO.
The vector of the best historical position of
the i3 th particle in PSO.
The complex voltage on bus i.

Ii
zij
Sij
Pij
si

The complex current to bus i.
The complex impedance between bus i and
bus j, zij = rij + ixij .
The branch flow Sij = Pij + iQij .
The real power Pij = |Iij |2 rij .
Power injection on bus i.
I. I NTRODUCTION

The network reconfiguration is a convenient way to reduce
the system loss, maintain the voltage profiles, and balance the
system load by opening and closing the limited number of
switches [1]–[4]. The developing load forecasting technology
can provide the load deviations in the periods between two
scheduled time points, which is usually ignored by traditional
network reconfiguration approaches. Considering the stochastic deviation of load profiles in distribution system, it is imperative to build a dynamic and efficient network reconfiguration
approach [5]–[10]. In this paper, a support vector regression
(SVR) based short-term load forecasting approach is employed
to cooperate with the network reconfiguration to minimize the
system loss.
Although most parameters in the SVR can be solved in a
convex manner, several parameters (can be defined as hyperparameters) cannot be determined by a similar way. The
optimization of these hyper-parameters is an indispensable
issue in SVR, which dramatically impacts the efficiency and
performance of the forecaster [11]–[14]. Specifically, a twostep parameters optimization approach is proposed with grid
traverse algorithm (GTA) and particle swarm optimization
(PSO). Considering the distributed computation frameworks
in [15]–[18], the hyper-parameter optimization of SVR is designed based on ”Mapreduce” idea to optimize the parameters
in parallel and reduce the computation time in this paper.
In this paper, we are focusing on the three-phase balanced
distribution system, and many heuristic methods are used to
solve this problem. However, they are easily trapped into
local minimums and heavily impacted by the initial points [7],
[19], [20]. To overcome this shortage, a distributed network
reconfiguration approach is proposed for the distribution systems with convex relaxation. Specifically, a second-order cone
program (SOCP) based approach is used to relax the threephase balanced power flow problem [21], [22]. Then, the
alternating direction methods of multipliers (ADMM) provides
a distributed method to solving the relaxed optimal power flow
problem [23]. With the decomposition and dual ascent alternating, the ADMM can solve the relaxed network reconfiguration

problem with a quickly convergence speed. Considering the
increasing computation capability, a parallel searching manner
is designed to traverse all the permutations of the switches, and
determine the best configuration to minimize the system loss.
The paper is organized as follows. In Section II, the
flowchart of the proposed approach is introduced. In Section III, the SVR based approach is used to short-term forecast
the distribution loads. In Section IV, based on the branch flow
model, the OPF problem is relaxed with SOCP and solved with
ADMM. In Section V, the numerical results are presented to
validate the proposed approach. The conclusion is presented
in Section VI.
II. T HE A RCHITECTURE

OF THE

can be built to minimize the forecast error with the soft margin
as follows:
n
X
1 T
(ξi1 + ξi∗1 )} (1)
{ ω ω+C
Rrisk =
min
ε,ω,ξi1 ,ξi∗ ,C,b,γ 2
1
i =1
1

Subject to
 ′

 Li1 − f (xi1 ) 6 ε + ξi1 ,
− L′i1 + f (xi1 ) 6 ε + ξi∗1 ,


ξi1 , ξi∗1 > 0.

(2)

where in (1),f is a Kernel based regression function, i1 is
a time index, 12 ω T ω indicates the flatness of the regression
coefficients, the second item is the tube violation, C is a tradeoff coefficient between the first two items, ξi1 and ξi∗1 indicates
the two training errors. In [11], the risk function (1) with the
constraint (2) can be derived to a dual problem with KarushKuhn-Tucker (KKT) condition. However, the parameters γ, C,
and ε are still need to be determined, which are the critical
factors to the performance of the forecaster [11], [24]. The
detail derivative of the SVM or SVR and its dual forms can
be found in [11], [24]. Then, as shown in Fig. 1, a two-step
based parameter optimization approach is designed to compute
the optimal parameters.

P ROPOSED A PPROACH

As shown in Fig. 1, the proposed approach consists of
two parts, the short-term load forecasting and distribution
system network reconfiguration. Both of them are designed
in distributed manner to effectively organize computation
resources and reduce the time consumption.
As shown in the left part of Fig. 1, in the part of shortterm load forecasting, the optimization of the hyper-parameters
contains GTA and PSO to avoid be trapped into local minimum. In the first step, the global solution space is spitted by
the GTA into the local solution spaces. In the Map phase,
because of the independency among the local solution spaces,
they are traversed by GTA in parallel. In the Reduce phase,
one or several local solution spaces are selected with the
minimum training errors. In the second step, the selected
local solution spaces are optimized by the PSO in the similar
manner. Then, the optimal parameters can be generated after
comparison in the Reduce phase. Finally, the short-term load
can be forecasted with the optimized hyper-parameters.
As shown in the right part of Fig. 1, in the part of distribution system network reconfiguration, a three-phase balanced
distribution system model is built with the forecasted load
profiles. The SOCP is used to relax the three-phase balanced
optimal power flow problem into a convex problem. In the
Map phase, the ADMM is used to compute the three-phase
balanced optimal power flow in parallel. Considering the very
high speed of ADMM, in the Reduce phase, all the statuses
of the switches are traversed and compared to get the global
optimization results. Considering the limited number of the
switches in the distribution system, all the permutations can
be traversed in a parallel manner to determine the optimal
configuration. Finally, the system loss can be reduced with
the optimized distribution reconfiguration results.

B. Two-step parameter optimization
Because the parameters γ, C, and ε cannot be solved
with the convex optimization, they are be defined as hyperparameters in [25]. There are several approaches are proposed
for the hyper-parameters such as random search, sequential
search, and Gaussian process [25], [26]. Considering the
complexity and feasibility, a grid traverse search based twostep hyper-parameter optimization is proposed for the SVR
based short-term load forecasting [12].
1) First Step: the GTA Procedure: As shown in Fig. 1,
the GTA procedure is the first step for the hyper-parameters
optimization, which aims to traverse the global solution space
into one or several local solution spaces. In the second step,
the local solution spaces can be searched with the PSO based
approach. The proposed approach is based on the increasing
computation capability and new computer cluster cooperation
soughs, for example Mapreduce. In the first step, the three
hyper-parameters are initialized with their upper bounds, lower
bounds, and grid searching steps. Then, a traversing vector H
can be generated as a finite multi-Cartesian product, which is
critical for the Mapreduce process. Hj2 is an element in H.
For each Hj2 , the loss function of SVR Rrisk can be computed
independently. As shown in Fig. 1, they can be computed in
parallel to reduce the computation time. In the last step, the
minimum Rrisk is selected. In addition, if several Hj2 are
selected, all of them are transmitted to the second step for PSO
optimization. In real-applications, if there is high requirement
for the time consumption, the GTA can provide good enough
parameter optimization.
2) Second Step: the PSO Procedure: In this paper, the PSO
procedure is designed as a “fine” optimization for the hyperparameters, which can be implemented as the scenarios with
less time consumption requirements.

III. SVR-BASED S HORT- TERM L OAD F ORECASTER

A. SVR Formulation
In this part, the SVR-based short-term load forecaster is
trained by the collected historical data to get the optimal hyperparameters. The objective function of the Kernel based SVR
2

Convex Relaxation

Data

Three-phase Unbalanced
Distribution System Network
Reconfiguration Model

Global Solution Space of SVR

Reduce

Local
Solution
Space
1

Local
Solution
Space
2

Local
Solution
Space

Work
Node 1

Work
Node 2

Work
Node q1

Second-order Cone Programming
Based Relaxation

p1

Objective Function of ADMM

GTA Results Comparison and Sublocal Solution Spaces Generation
Map

Particle Swarm Optimization

Map

Local
Solution
Space
1

Local
Solution
Space

Work
Node 1

Work
Node q2

Local
Solution
Space
1

Local
Solution
Space
2

Local
Solution
Space

Work
Node 1

Work
Node 2

Work
Node q3

p3

p2

Optimal AC Power Flow Results
Comparison

Reduce

Optimal Parameters
Determination

Reduce

Alternating Direction Method of Multipliers

Grid Traverse Algorithm

Map

Final Results

Forecasting with the
Optimal Parameters

Short-term Load Forecasting

Distribution System Network Reconfiguration

Fig. 1. The flowchart of the proposed approach.

Algorithm 1 GTA for Hyper-parameter Optimization
Objective: Shrink the global solution spaces into one or
several local solution spaces.

IV. D ISTRIBUTION S YSTEM N ETWORK
R ECONFIGURATION

Initialization: Hyper-parameters initialization and multiCartesian product generation for the GTA.

The topology of a distribution system can be represented in a
graph with buses and branches: G = [V, E]. Then the branch
flow model can be built as follows [21], [23]:
X
X
Sij −
(Ski − lki zki ),
(4a)
si =

Grid Traverse Searching: For each core or each process, assign the Hj2 to compute the Rrisk , which can be
computed in parallel with the Mapreduce model.

j

k

2
vj = vi − 2(rij Pij + xij Qij ) + (rij
+ x2ij )lij ,

Determine Local Solution Space: Collected all the
results, and select the local solution spaces with minimum
Rrisk .

lij =

αΩ
i3 (t)

Ω
+ ϕ2 θ2 (η Ω
g − αi3 (t
Ω
=αΩ
i3 (t − 1) + ν i3 (t),

+
2

Q2ij )/vj ,

(4b)
(4c)

2

where lij := |Iij | , vi := |Vi | , Sij , Pij , Qij and zij indicate
the complex power flow, active power, reactive power, and
impedance on branch ij ∈ E, Sij = Pij + iQij , Pij = |Iij |2 rij
and zij = rij + ixij .
During the operation of network reconfigurations, the topology of the distribution system is keeping radial and avoid
any loops. Considering the characteristics of the three-phase
balanced distribution system, the SOCP relaxation inequalities
can be represented as follows [21], [23]:

After initialization of the particles, for particle i3 , its velocity and position can be computed as following:
Ω
Ω
Ω
νΩ
i3 (t) =ν i3 (t − 1) + ϕ1 θ1 (η i3 − αi3 (t − 1))

(Pij2

(3a)

− 1)),
(3b)

|Sij |2
≤ lij ,
vi

where the acceleration coefficients are defined as ϕ1 and ϕ2 ,
θ1 and θ2 can be seemed as two independently weightiness
coefficients, the best historical position and the best position
Ω
are defined as η Ω
i3 and η g , respectively, and hyper-parameter
Ω
Ω = [γ C ε]. αi3 (t) and ν Ω
i3 (t) are the position and velocity
vectors, respectively.

(5)

where (5) can be used to instead of (4c) as the inequalities
constraints. In this paper, the objective function is defined as
total line loss as follows:
X
F =
Pij ,
(6)
E

3

25

Density(%)

20

81.21%

15

10

5

0

-20

-15

-10

-5

-3

0

3

5

10

15

20

Forecast Error (%)

(a)

(b)

Fig. 2. (a) The IEEE 123-bus based distribution system, (b) the error distribution of short-term load forecasting.

No.
1.
2.
3.
4.
5.
6.
7.
8.

Scenario
Load Increasing
Load Increasing
Load Increasing
Load Increasing
Load Decreasing
Load Decreasing
Load Decreasing
Load Decreasing

TABLE I
R ESULTS OF N ETWORK R ECONFIGURATION
Bus No. Opened Switches
Original PLoss
New PLoss
83
TS-3, TS-4
54.2 Kw
36.7 Kw
300
TS-2, TS-3
42.4 Kw
31.5 Kw
95
TS-3, TS-4
78.5 Kw
67.0 Kw
49
TS-1, TS-4
17.6 Kw
13.4 Kw
47
TS-2, TS-4
39.4 Kw
22.3 Kw
108
TS-1, TS-4
24.1 Kw
20.1 Kw
250
TS-3, TS-4
29.5 Kw
26.8 Kw
56
TS-1, TS-2
35.2 Kw
33.7 Kw

Then, considering the independency of each configuration
(permutation), all the permutation can be implemented into
different cores or processes and computed independently,
which dramatically reduces the computation time and keeps
the convexity to get the guaranteed optimization results.

where the constraints contain (4a), (4b), (5), and the basic
physical constraints:
Vi,min ≤ Vi ≤ Vi,max ,

(7a)

Iij ≤ Iij,max .

(7b)

Considering the ADMM, the objective function (6) with the
constraints (4a), (4b), (5), (7a), and (7b) can be decomposed
into a dual problem. The detail derivatives of the ADMM can
be found in [21], [23].
During the parallel traverse of all statuses of the switches,
the topology of the distribution system is keeping as a radial
network without any loops, which can be formulated as [1]:
rank(A) = N − d,
X
aij = N − d,

V. N UMERICAL R ESULTS
As shown in Fig. 2(a), the test bench is based on the
IEEE 123-bus distribution system, and four initially opened tie
switches TS-1, TS-2, TS-3 and TS-4 are added to make the
system topology changeable, and the detail information can
be found in [1], [27]. A test load data set is from a partner
utility’s distribution feeder.

(8a)
(8b)

A. Short-term Load Forecasting

E

X

aij = Mk − 1.

Loss Reduction
32.28 %
25.71 %
14.65 %
23.86 %
43.40 %
16.59 %
9.15 %
4.27 %

The test load data contains four seasons of one year, and 30
days are selected for each season. The proposed SVR based
load forecasting approach is used for 1-hour-ahead sliding
window forecasting with 1 second resolution. The training data
is 5 times as the test data. The distribution of the forecaster
errors is shown in Fig. 2(b). The mean absolute precentage
error (MAPE) is 2.23%, normalized root-mean-square error
(NRMSE) is 4.03%, and more than 80% of the errors are
accumulated between (-3.1%, 3.1%).

(8c)

Ek

where A is the adjacency matrix of the graph G, d is the
number of slack bus, N is the number of buses V, Mk is the
number of branches in path Ek , and aij is an element of A:
(
1, if bus i and bus j are connected,
aij =
(9)
0, else.
Considering the limited number of switches, the proposed
approach is designed to traverse all the permutations and
determine the optimal configuration of the distribution system.
For example, the modified IEEE 123-bus system with 4
switches indicates 16 scenarios with all the permutations of the
switches [1]. With the topology constraints discussed above,
the number of scenarios can be reduced in different scenarios.

B. Network Reconfiguration
As shown in Table I, considering the load increasing, for
example, in scenario 1, with the forecasting results, there
is a load increasing 20.31% in bus 83. The system loss
reduces 32.28% with the proposed approach. For the load
decreasing, for example, in scenario 5, with the forecasting
4

[7] Y. Gu, H. Jiang, Y. Zhang, and D. W. Gao, “Statistical scheduling of
economic dispatch and energy reserves of hybrid power systems with
high renewable energy penetration,” in 2014 48th Asilomar Conference
on Signals, Systems and Computers, 2014, pp. 530–534.
[8] M. Cui, J. Zhang, Q. Wang, V. Krishnan, and B.-M. Hodge, “A datadriven methodology for probabilistic wind power ramp forecasting,”
IEEE Transactions on Smart Grid, 2017.
[9] X. Zhu, J. Yan, and N. Lu, “A probabilistic-based pv and energy storage
sizing tool for residential loads,” in Transmission and Distribution
Conference and Exposition (T&D), 2016 IEEE/PES. IEEE, 2016, pp.
1–5.
[10] Y. Chen and J. Borken-Kleefeld, “No x emissions from diesel passenger
cars worsen with age,” Environmental science & technology, vol. 50,
no. 7, pp. 3327–3332, 2016.
[11] C. J. Burges, “A tutorial on support vector machines for pattern
recognition,” Data Mining and Knowledge Discovery, vol. 2, no. 2, pp.
121–167, 1998.
[12] H. Jiang, Y. Zhang, E. Muljadi, J. Zhang, and W. Gao, “A short-term
and high-resolution distribution system load forecasting approach using
support vector regression with hybrid parameters optimization,” IEEE
Transactions on Smart Grid, 2016.
[13] Y. Chen and J. Borken-Kleefeld, “Real-driving emissions from cars
and light commercial vehicles–results from 13 years remote sensing at
zurich/ch,” Atmospheric Environment, vol. 88, pp. 157–164, 2014.
[14] Y. Jiang, J. Xu, Y. Sun, C. Wei, J. Wang, D. Ke, X. Li, J. Yang,
X. Peng, and B. Tang, “Day-ahead stochastic economic dispatch of wind
integrated power system considering demand response of residential
hybrid energy system,” Applied Energy, vol. 190, pp. 1126–1137, 2017.
[15] M. Zaharia, M. Chowdhury, M. J. Franklin, S. Shenker, and I. Stoica,
“Spark: Cluster computing with working sets.” HotCloud, vol. 10, no.
10-10, p. 95, 2010.
[16] Y. Gu, H. Jiang, Y. Zhang, J. J. Zhang, T. Gao, and E. Muljadi,
“Knowledge discovery for smart grid operation, control, and situation
awarenessła big data visualization platform,” in North American Power
Symposium (NAPS), 2016. IEEE, 2016, pp. 1–6.
[17] W. Yan, W. Gao, T. Gao, D. W. Gao, S. Yan, and J. Wang, “Distributed
cooperative control of virtual synchronous generator based microgrid,”
2017.
[18] H. Jiang, L. Huang, Y. Zhang, J. J. Zhang, and D. W. Gao, “Spatialtemporal characterization of synchrophasor measurement systems - A
big data approach for smart grid system situational awareness,” in 2014
Conference Record of the Forty Eighth Asilomar Conference on Signals,
Systems and Computers (ASILOMAR), 2014, pp. 750–754.
[19] Z. Tang and K. K. Bagchi, “Globally convergent particle swarm optimization via branch-and-bound,” Computer and Information Science,
vol. 3, no. 4, pp. 60–71, 2010.
[20] F. He, Y. Gu, J. Hao, J. J. Zhang, J. Wei, and Y. Zhang, “Joint real-time
energy and demand-response management using a hybrid coalitionalnoncooperative game,” in 2015 49th Asilomar Conference on Signals,
Systems and Computers. IEEE, 2015, pp. 895–899.
[21] Q. Peng and S. H. Low, “Distributed algorithm for optimal power flow
on a radial network,” in 2014 IEEE 53rd Annual Conference on Decision
and Control (CDC). IEEE, 2014, pp. 167–172.
[22] H. Jiang, J. J. Zhang, D. W. Gao, Y. Zhang, and E. Muljadi, “Synchrophasor based auxiliary controller to enhance power system transient
voltage stability in a high penetration renewable energy scenario,” in
2014 IEEE Symposium Power Electronics and Machines for Wind and
Water Applications (PEMWA), 2014, pp. 1–7.
[23] S. H. Low, “Convex relaxation of optimal power flowłpart i: Formulations and equivalence,” IEEE Transactions on Control of Network
Systems, vol. 1, no. 1, pp. 15–27, 2014.
[24] C.-C. Chang and C.-J. Lin, “Libsvm: a library for support vector
machines,” ACM Transactions on Intelligent Systems and Technology
(TIST), vol. 2, no. 3, p. 27, 2011.
[25] J. Bergstra and Y. Bengio, “Random search for hyper-parameter optimization,” Journal of Machine Learning Research, vol. 13, no. Feb, pp.
281–305, 2012.
[26] J. S. Bergstra, R. Bardenet, Y. Bengio, and B. Kégl, “Algorithms
for hyper-parameter optimization,” in Advances in Neural Information
Processing Systems, 2011, pp. 2546–2554.
[27] F. Ding and K. A. Loparo, “Hierarchical decentralized network reconfiguration for smart distribution systemsłpart i: Problem formulation and
algorithm development,” IEEE Transactions on Power Systems, vol. 30,
no. 2, pp. 734–743, 2015.

TABLE II
P ERFORMANCE COMPARISON
Methods
Loss Reduction
Computation Time (s)
GA based traditional
17.87%
107
Proposed Approach
21.24%
30

results, there is a load decreasing 55.4% in bus 47. The system
loss reduces 43.40% with the proposed approach. The average
loss reduction for the load increasing scenarios is 24.13%,
the average loss reduction for the load decreasing scenarios is
18.35%, and the total average for all scenarios is 21.24%.
C. Comparison
As shown in Table II, compared with the traditional network
reconfiguration approach with the genetic algorithm (GA), the
proposed approach has less computation time and more loss reduction. Furthermore, the proposed approach is more intuitive,
and convenient for implementation in different programming
language such as python and Matlab.
VI. C ONCLUSION
In this paper, a short-term load forecasting based network
reconfiguration is proposed to reduce the distribution system
loss dynamically. Instead of the static load measurements
at the scheduled time spots, the short-term load forecasting
approach can provide the accurate future load profiles, which
contains more information for the network reconfiguration.
With the SOCP, the OPF of the three-phase balanced distribution system can be relaxed into a convex problem and solved
with ADMM. The optimal reconfiguration is generated by the
parallel computation with traversing all the permutations of
the switche status. The whole proposed approach is designed
as a distributed computation approach.
In the next step, with the increasing penetration of renewable energies, the wind and solar power will be considered,
and other system such as traffic system with electrical vehicles
will be also considered in the future research.
R EFERENCES
[1] F. Ding, H. Jiang, and J. Tan, “Automatic distribution network reconfiguration: An event-driven approach,” in Power and Energy Society General
Meeting (PESGM), 2016. IEEE, 2016, pp. 1–5.
[2] Y. Luo, M. Panwar, M. Mohanpurkar, and R. Hovsapian, “Real time
optimal control of supercapacitor operation for frequency response,” in
Power and Energy Society General Meeting (PESGM), 2016. IEEE,
2016, pp. 1–5.
[3] L. Bai, F. Li, H. Cui, T. Jiang, H. Sun, and J. Zhu, “Interval optimization
based operating strategy for gas-electricity integrated energy systems
considering demand response and wind uncertainty,” Applied Energy,
vol. 167, pp. 270–279, 2016.
[4] X. Fang, Q. Hu, F. Li, B. Wang, and Y. Li, “Coupon-based demand
response considering wind power uncertainty: A strategic bidding model
for load serving entities,” IEEE Transactions on Power Systems, vol. 31,
no. 2, pp. 1025–1037, 2016.
[5] Y. Zhang, R. Yang, K. Zhang, H. Jiang, and J. J. Zhang, “Consumption
behavior analytics-aided energy forecasting and dispatch,” IEEE Intelligent Systems, vol. 32, no. 4, pp. 59–63, 2017.
[6] M. Cui, D. Ke, Y. Sun, D. Gan, J. Zhang, and B.-M. Hodge, “Wind
power ramp event forecasting using a stochastic scenario generation
method,” IEEE Transactions on Sustainable Energy, vol. 6, no. 2, pp.
422–433, 2015.

5

