Next Article in Journal
Detection of Sensitive Data to Counter Global Terrorism
Previous Article in Journal
Shock Reduction through Opposing Jets—Aerodynamic Performance and Flight Stability Perspectives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Two-Stage Stochastic Programming Scheduling Model for Hybrid AC/DC Distribution Network Considering Converters and Energy Storage System

1
State Grid Energy Saving Service Co., Ltd., Beijing 100052, China
2
Nanrui group company/State Grid Electric Power Research Institute, Nanjing 210000, China
3
State Grid Electric Power Research Institute Wuhan Efficiency Evaluation Co., Ltd., Wuhan 430074, China
4
School of Electrical Engineering and Automation, Fuzhou University, Fuzhou 350116, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2020, 10(1), 181; https://doi.org/10.3390/app10010181
Submission received: 14 November 2019 / Revised: 12 December 2019 / Accepted: 20 December 2019 / Published: 25 December 2019
(This article belongs to the Section Energy Science and Technology)

Abstract

:
The development of DC distribution network technology has provided a more efficient way for renewable energy accommodation and flexible power supply. A two-stage stochastic scheduling model for the hybrid AC/DC distribution network is proposed to study the active-reactive power coordinated optimal dispatch. In this framework, the wind power scenario set is utilized to deal with its uncertainty in real time, which is integrated into the decision-making process at the first stage. The charging/discharging power of ESSs and the transferred active/reactive power by VSCs can be adjusted when wind power uncertainty is observed at the second stage. Moreover, the proposed model is transformed into a mixed integer second-order cone programming optimization problem by linearization and second-order cone relaxation techniques to solve. Finally, case studies are implemented on the modified IEEE 33-node AC/DC distribution system and the simulation results demonstrate the effectiveness of the proposed stochastic scheduling model and solving method.

1. Introduction

With the rapid growth of renewable energy and energy storage systems integrated in the distribution network, the development of a DC distribution network has drawn more and more attention [1]. Compared to the AC distribution network, the DC distribution network has the advantages of reducing active power losses, and improving the power quality and utilization rate of electric energy [2,3]. As a result, the hybrid AC/DC distribution network will become a new tendency of intelligent distribution network development in the future [4,5]. By combining the advantages of two kinds of distribution networks, the hybrid distribution network can achieve flexible power supply for different types of electric load and renewable energy accommodation [6,7].
To handle the variations of renewable energy and load demand, a multi-timescale coordinated stochastic voltage/var control approach was developed and a mixed-integer quadratic programming model was built for the IEEE 33-bus radial AC distribution [8]. The authors of [9] proposed a day-ahead reactive power dispatch method of AC distribution networks by considering power forecast errors of renewable energy, and a dynamic preliminary-coarse-fine adjustment strategy was designed to achieve the optimal scheduling of distributed generators (DGs) and capacitor banks (CBs).
At present, research of the hybrid AC/DC distribution network is in the elementary stage and needs to expand and deepen. Voltage source converters (VSCs) have been utilized to convert AC lines into DC lines to improve DG accommodations [10]. To realize the coordinated voltage regulation in a hybrid AC/DC distribution system, a priority-based real-time control strategy was proposed according to the voltage control effect of active-reactive power adjustment in [11]. Multiple battery energy storage systems (ESSs) have been considered in distribution networks for real-time voltage regulation in [12]. Furthermore, [13] proposed an optimal control management of ESSs to mitigate the fluctuation and intermittence of renewable energy. Taking into consideration prediction errors, the optimal operation of ESSs in a distribution system was developed and solved in a two-level framework to alleviate the net load uncertainties [14,15]. For the planning problems of ESSs, [16] introduced a new index to quantify wind power fluctuations, and the impact of different ESSs configurations on the solutions were analyzed. A scenario-based chance-constrained planning approach was developed to handle the joint planning of multiple technologies of ESSs, and an easy-to-implement variant of Benders decomposition algorithm was proposed to solve the chance-constrained optimization problem [17]. In order to improve the reliability of electricity supplies with wind power integration, a population meta-heuristics algorithm was adopted to determine the minimum capacity of kinetic energy storage in [18].
To consider the randomness and uncertainty of wind power output, this paper proposes a two-stage stochastic scheduling model to realize active-reactive power coordinated economic dispatch for the hybrid AC/DC distribution network with the integration of capacitor banks and an energy storage system. The first stage is to determine the power setpoints of generator units, switching number of CBs, and charging/discharging state of ESSs, and the second stage is to adjust the transferred active/reactive power of VSCs and charging/discharging power of ESSs based on the realized wind power scenarios. To sum up, the schematic framework of this study is shown in Figure 1, for which the corresponding research contents are described in the following sections.
The contributions of this paper can be summarized as follows:
(1) The capacitor bank, energy storage system, and voltage source converter are coordinated as the active adjustment measures simultaneously in the hybrid AC/DC distribution system. By this means, a comprehensive and practical framework considering active and reactive powers coordination dispatch is established.
(2) A two-stage stochastic programming model is established to achieve active-reactive power coordinated economic dispatch for the hybrid distribution system, in which wind power uncertainty is described as several wind power scenarios obtained by scenario generation and reduction techniques. As a result, the decision-making at the first stage not only depends on the day-ahead wind power predictions but is also affected by wind power uncertainty in the real-time scheduling stage. In other words, the realized wind power scenario set in real time is integrated into the decision process for the two-stage scheduling problem.
(3) The proposed stochastic programming model is a nonlinear optimization problem, which is transformed into a mixed integer second-order cone programming problem by the linearization and second-order cone relaxation techniques for solving. Moreover, simulation results about the modified IEEE 33-node distribution system demonstrate the effectiveness of the proposed stochastic scheduling model.
The remainder of this paper is organized as follows. The power flow models for the AC distribution network, DC distribution network, and voltage source converter are established in Section 2. The two-stage stochastic scheduling model for the hybrid AC/DC distribution system is proposed in Section 3. Case studies and simulation results about the modified IEEE 33-node distribution system are shown in Section 4. Finally, conclusions are drawn in Section 5.

2. Power Flow Model for Hybrid AC/DC Distribution Network

2.1. AC Distribution Network Model

The DistFlow equations are used to describe the power flows in AC distribution system [19]:
P i , t AC = k ( i , : ) AC P k , t AC l ( : , i ) AC [ P l , t AC ( I l , t AC ) 2 r l ]   i N AC ,
Q i , t AC = k ( i , : ) AC Q k , t AC l ( : , i ) AC [ Q l , t AC ( I l , t AC ) 2 x l ]   i N AC ,
( U i , t AC ) 2 = ( U j , t AC ) 2 2 ( r k P k , t AC + x k Q k , t AC ) + ( I k , t AC ) 2 ( r k 2 + x k 2 ) k ( j , i ) AC ,
( I k , t AC ) 2 ( U j , t AC ) 2 = ( P k , t AC ) 2 + ( Q k , t AC ) 2 k ( j , i ) AC ,
( P k , t AC ) 2 + ( Q k , t AC ) 2 ( S k , max AC ) 2 k ( j , i ) AC ,
U i , min AC U i , t AC U i , max AC   i N AC ,
where AC represents the set of AC branches; 𝒩AC represents the set of AC nodes; t is the scheduling time; k(i,;) represents branch k starting with node i and l(:,i) represents the branch l ending with node i; rl and xl are the resistance and reactance of branch l; P i , t AC and Q i , t AC are the injection active and reactive power at node i; I k , t AC , P k , t AC , and Q k , t AC are the phase current, active, and reactive power at branch k; U j , t AC and U i , t AC are the phase voltage at node j and i; U i , min AC and U i , max AC are the minimum and maximum voltage at node i; and S k , max AC is the maximum load flow of AC branch k.

2.2. DC Distribution Network Model

The DC distribution system can be described as:
P i , t DC = k ( i , : ) DC P k , t DC l ( : , i ) DC [ P l , t DC ( I l , t DC ) 2 r l ]   i N DC ,
( U i , t DC ) 2 = ( U j , t DC ) 2 2 r k P k , t DC + r k 2 ( I k , t DC ) 2 k ( j , i ) DC ,
( I k , t DC ) 2 ( U j , t DC ) 2 = ( P k , t DC ) 2 k ( j , i ) DC ,
U i , min DC U i , t DC U i , max DC i N DC ,
S k , max DC P k , t DC S k , max DC k ( j , i ) DC ,
where DC represents the set of DC branches; 𝒩DC represents the set of DC nodes; P i , t DC is the injection active power at node i; I k , t DC and P k , t DC are the phase current and active power at branch k; U i , min DC and U i , max DC are the minimum and maximum voltage at DC node i; and S k , max DC is the maximum load flow of DC branch k.

2.3. Voltage Source Converter Model

The schematic diagram of the voltage source converter (VSC) is shown in Figure 2.
As shown in Figure 1, r i VSC and x i VSC represent the equivalent resistance and reactance of the VSC connected with node j in the AC network. The branch from node j to AC side node i of the VSC can be regarded as the AC branch in the AC network. The VSC model can be described as [20,21]:
P i , t VSC ( AC ) = P k , t AC ( I k , t AC ) 2 r i VSC i N SVC ,
Q i , t SVC ( AC ) = Q k , t AC ( I k , t AC ) 2 x i VSC i N SVC ,
( U i , t VSC ( AC ) ) 2 = ( U j , t AC ) 2 2 ( r i VSC P k , t AC + x i VSC Q k , t AC ) + ( I k , t AC ) 2 [ ( r i VSC ) 2 + ( x i VSC ) 2 ] ,
3 P i , t VSC ( AC ) = P i , t VSC ( DC ) i N SVC ,
U i , t VSC ( AC ) = 3 3 μ M i , t U i , t VSC ( DC ) i N SVC ,
where P i , t VSC ( AC ) and Q i , t SVC ( AC ) are the input active and reactive power at the AC side of the VSC while P i , t VSC ( DC ) is the output active power at the DC side of the VSC; U i , t VSC ( AC ) and U i , t VSC ( DC ) are the voltage at the AC and DC sides of VSC i; μ is the utilization rate of the DC voltage; and Mi,t is the modulation degree of SVC i within the range of 0 to 1.

2.4. Second-Order Cone Relaxation for the Distribution System Model

For the above model of the hybrid AC/DC distribution network, and the nonlinear terms existing in Equations (1)–(5), (7)–(9), and (12)–(14), the following variables are defined to eliminate the quadratic terms of voltage and current for the AC/DC power flow.
( I k , t AC ) 2 = I ˜ k , t AC ,   ( I k , t DC ) 2 = I ˜ k , t DC ,
( U i , t AC ) 2 = U ˜ i , t AC ,   ( U i , t DC ) 2 = U ˜ i , t DC .
Then, Equations (4), (5), and (9) can be converted to second-order cone constraints as follows [22,23]:
[ 2 P k , t AC 2 Q k , t AC I ˜ k , t AC U ˜ j , t AC ] T 2 I ˜ k , t AC + U ˜ j , t AC   ,
[ P k , t AC Q k , t AC ] T 2 S k , max AC ,
[ 2 P k , t DC I ˜ k , t DC U ˜ j , t DC ] T 2 I ˜ k , t DC + U ˜ j , t DC .
By this means, the hybrid AC/DC distribution system model is transformed as a second-order cone programming problem.

3. Two-Stage Stochastic Scheduling Model for the Hybrid AC/DC Distribution Network

The two-stage stochastic scheduling model is established to cope with wind power uncertainty. The decision at the first stage is applicable for all the wind power scenarios, and the respective adjustment decision for each wind power scenario is made at the second stage.

3.1. Scenario Generation and Reduction

The wind power prediction error is regarded as a random variable, which obeys a normal distribution, then the wind power output in real time can be expressed by [24]:
P ˜ t W = P ^ t W + Δ P ˜ t W ,   Δ P ˜ t W ~ G ( 0 , δ t 2 ) ,
where Δ P ˜ t W is the wind power prediction error at time t; P ^ t W is the day-ahead wind power prediction value; δt is the standard deviation of the prediction error; and G represents the normal distribution function.
The scenario generation and reduction techniques are implemented to handle wind power uncertainty. Since the Monte Carlo simulation method based on simple random sampling suffers from a higher computational burden, the Latin hypercube sampling (LHS) is adopted to obtain the original wind power scenario set, which can achieve good coverage from the entire distribution of the random variable [25].
The procedures of LHS can be described as follows:
Step 1: Assuming the number of samples for the random variable, Δ P ˜ t W , is N, the interval [0, 1] can be equally divided into N intervals, and a uniformly distributed random number within the range of [(k−1)/N, k/N] is generated for k = 1 to N;
Step 2: Disorder the order of the above n random numbers obtained by Step 1;
Step 3: The n random numbers by Step 2 are regarded as the probability values of samples, and the value of the random variable, Δ P ˜ t W , can be obtained according to the inverse function of the cumulative distribution function (CDF) for Δ P ˜ t W :
Δ P t ( k ) = F t 1 ( k 0.5 N ) ,   k = 1 , , N ,
where Ft(·) represents the CDF of random variable Δ P ˜ t W , and Δ P t ( k ) is the ith sample value of Δ P ˜ t W .
Cholesky decomposition (CD) is introduced to alleviate the undesired correlations between random variables at different time points [26].
The sample matrix for Δ P ˜ t W by LHS can be expressed as:
Δ P = [ Δ P 1 ( 1 ) Δ P 1 ( 2 ) Δ P t ( N ) Δ P 2 ( 1 ) Δ P 2 ( 2 ) Δ P t ( N ) Δ P T ( 1 ) Δ P T ( 2 ) Δ P t ( N ) ] T × N .
The procedures of CD are as follows:
Step 1: A T × N ordering matrix, L, is generated, and the element at the tth row and ith column represents the position of the sample, Δ P t ( k ) , in the matrix, Δ P , to be arranged;
Step 2: The elements of each row in Δ P are arranged according to the ordering matrix, L;
Step 3: The corresponding correlation matrix of L noted as ρL can be calculated and decomposed by the Cholesky decomposition:
ρ L = D D T .
Step 4: A T × N matrix, G, can be obtained by:
G = D 1 L .
Step 5: The elements of each row in the matrix, L, are arranged according to the order of elements in the corresponding row in G;
Step 6: The elements of each row in the matrix, Δ P , are arranged according to the updated ordering matrix, L.
Then, the simultaneous backward reduction method is utilized to decrease the number of wind power scenarios [27]. For the matrix, Δ P , obtained by the LHS and CD techniques, the original scenario k noted as Δ P ( k ) is a vector composed of the elements of the kth row in Δ P , for which the occurrence probability is π(k). The Kantorovich distance of scenarios Δ P ( k )  and   Δ P ( j ) can be defined as:
d k , j = Δ P ( k ) Δ P ( j ) 2 .
The scenario number in the reduced scenario set is assigned as Ns. The procedures of the simultaneous backward reduction method can be described as follows:
Step 1: Delete the scenario, k, which satisfies the following condition:
π ( k ) π ( j ) min k j d k , j = min m { 1 , , N } π ( m ) { min n { 1 , , N } , n m π ( n ) d n , m } .
Step 2: The scenario number is set by N = N − 1, and the scenario, k*, closest to the deleted scenario, k, is chosen according to:
d k , k = min s k π ( k ) π ( s ) d s , k .
Step 3: The occurrence probability of scenario k* is set as:
π ( k ) = π ( k ) + π ( k ) .
Step 4: If N > Ns, return to Step 1; otherwise, the scenario reduction process is finished.
As a result, the reduced scenario set for wind power can be obtained as:
P = { P t , s W , s = 1 , 2 , , N s } .
The parameters, N and Ns, are set to 1000 and 10, respectively. The standard deviation of the wind power prediction error at time point t is set to 20% of P ^ t W . The original and reduced wind power scenario sets are shown in Figure 3.

3.2. Objective Function

The day-ahead scheduling cost includes the generating cost of gas-fired units in the first stage and the power purchasing cost, penalty cost of power fluctuation from the substation, network loss cost, load shedding, and wind curtailment cost:
min t i c G P i , t G + s π s t c L ( i N AC 3 L i , t , s Cur + i N DC L i , t , s Cur ) + s π s t c W ( i N AC 3 W i , t , s Cur + i N DC W i , t , s Cur ) + s π s t c Loss ( k AC 3 ( I k , t , s AC ) 2 r k + k DC ( I k , t , s DC ) 2 r k ) + s π s t i Ω Sub ( c t Sub P i , t , s Sub + c Ramp | P i , t , s Sub P i , t 1 , s Sub | )
where cG is the unit generating cost of gas-fired units; cL is the load shedding price; cW is the wind curtailment price; cLoss is the network loss price; c t Sub is the electricity price of the main grid at time t; cRamp is the penalty price for power fluctuation; ΩSub is the node set of substation; and πs is the occurrence probability of the wind power scenario, s.

3.3. Constraints for the Two-Stage Stochastic Programming

3.3.1. Operation Constraints of the Capacitor Bank

Q i , t CB = y i , t CB Q step CB ,
0 y i , t CB N max CB ,
t | y i , t + 1 CB y i , t CB | y max CB ,
where y i , t CB is the operation number of CB; N max CB is the maximum operation number of CB; Q step CB is the capacity of single group CB; and y max CB is the maximum adjustable number of CB in one day.
It can be seen that the operation number of CB is determined in the first stage, which is applicable for all wind power scenarios in the second stage. In addition, the nonlinear constraint (Equation (35)) can be linearized as follows [28]:
{ y i , t + 1 CB y i , t CB N max CB b t CB y i , t CB y i , t + 1 CB N max CB b t CB t b t CB y max CB ,
where b t CB is the binary variable indicating whether the operation number of CB has changed.

3.3.2. Operation Constraint of the Energy Storage System

y i , t ch + y i , t dis 1 ,
0 P i , t , s dis y i , t dis P max ESS ,
0 P i , t , s ch y i , t ch P max ESS ,
E i , t + 1 , s ESS = E i , t , s ESS + η ch P i , t , s ch η dis P i , t , s dis ,
0.2 E max ESS E i , t , s ESS 0.9 E max ESS ,
t | y i , t + 1 ch y i , t ch | y max ESS ,
t | y i , t + 1 dis y i , t dis | y max ESS ,
where y i , t ch and y i , t dis represent the charging and discharging status of the energy storage system; ηch and ηdis are the charging and discharging efficiency; P max ESS is the maximum charging/discharging power; E max ESS is the maximum storage capacity; E i , t , s ESS is the storage capacity of ESS at time t under scenario s; P i , t , s dis and P i , t , s ch represent the discharging and charging power output of ESS at time t under scenario s; and y max ESS is the maximum charging/discharging number in one day.
It can be seen that the charging/discharging state of ESS noted as y i , t ch and y i , t dis is not affected by wind power fluctuations under different scenarios, which remains unchanged during the whole scheduling period. On the other hand, the charging and discharging power output under wind power scenario s are defined as P i , t , s ch and P i , t , s dis , respectively, which indicates that the charging/discharging power of ESS depends on the possible wind power scenarios. As a result, the charging/discharging state of ESS is regarded as the decision variables in the first stage, which is determined before the possible wind power scenarios happened. Moreover, the charging/discharging power output of ESS is considered as the decision variable in the second stage, which is used to alleviate the wind power fluctuation in real time.

3.3.3. Relevant Constraints of Two Stages

For the two-stage scheduling model, the relevant constraints of the first and second stages are described as follows:
P i , t , s AC = P i , t , s Sub + P i , t , s W + P i , t G + P i , t , s dis P i , t , s ch P i , t L i N AC ,
Q i , t , s AC = Q i , t , s Sub + Q i , t , s W + Q i , t G + Q i , t CB Q i , t L i N AC ,
P i , t DC = P i , t , s W + P i , t , s VSC ( DC ) P i , t L i N DC ,
where P i , t , s Sub and Q i , t , s Sub represent the active and reactive power from the substation connected with node i at time t under scenario s; P i , t , s W and Q i , t , s W represent the active and reactive power from the wind generator under scenario s; P i , t G and Q i , t G represent the active and reactive power output by the gas-fired units; and P i , t L and Q i , t L represent the active and reactive power load at node i.
In this study, the gas consumption plan by gas-fired units is determined before wind power uncertainty is observed and will not adjust in real time. Consequently, the power output of gas-fired units cannot be changed under different wind power scenarios at the second stage.
Furthermore, the hybrid AC/DC distribution network model in Section 2 is established for the certain wind power scenario. In other words, all the variables in Equations (1)–(21) should be changed to the corresponding variables for different wind power scenarios. Taking the DC distribution network as an example, U i , t DC , I k , t DC , and P k , t DC should be replaced by U i , t , s DC , I k , t , s DC , and P k , t , s DC , respectively. Additionally, the power flow model for the DC distribution network under scenario s can be expressed as follows:
{ P i , t , s DC = k ( i , : ) DC P k , t , s DC l ( : , i ) DC ( P l , t , s DC I ˜ l , t , s DC r l )   i N DC U ˜ i , t , s DC = U ˜ j , t , s DC 2 r k P k , t , s DC + I ˜ k , t , s DC r k 2 k ( j , i ) DC [ 2 P k , t , s DC I ˜ k , t , s DC U ˜ j , t , s DC ] T 2 I ˜ k , t , s DC + U ˜ j , t , s DC ( U i , min DC ) 2 U ˜ i , t , s DC ( U i , max DC ) 2   i N DC S k , max DC P k , t , s DC S k , max DC k ( j , i ) DC .
The power flow model for the AC distribution network and VSC under the wind power scenario, s, can be expressed in a similar way, but the corresponding formulas are not listed here to save page.

4. Case Studies

The case study was conducted based on the modified IEEE 33-node hybrid AC/DC distribution system [19]. As shown in Figure 4, the solid red lines represent AC branches and the dashed blue lines represent DC branches. The rated voltages of the AC and DC distribution network are 12.66 and 15 kV, respectively. The voltage limitation range was set to [0.95, 1.05] p.u. There are two gas-fired units with the capacity of 300 kW, and three VSCs and four wind generators with the capacity of 500 kW in the hybrid system. The electricity price of main grid is 0.12 $/kWh at hour 7–22 and 0.06 $/kWh at other times. The operation parameters of different facilities and the cost coefficients are listed in Table 1, Table 2, Table 3, Table 4 and Table 5. The unit generating cost of gas-fired units and power fluctuation cost in Table 5 were obtained from [29], the unit network loss cost was obtained from [21], and the unit penalty costs for wind curtailment and load shedding can be found in [30]. The load demand and wind power curves are shown in Figure 5. The power output predictions of four wind generators are the same, and the wind generators in the AC distribution network are running at a constant power factor of 0.85. In this study, the solving tool Yalmip was adopted to develop the above application programs and then solved by Cplex in the MATLAB platform [31].
The following cases were considered to demonstrate the effectiveness of proposed model:
Case 1:
Day-ahead scheduling for the hybrid distribution system with CB considering the wind power prediction curve;
Case 2:
Day-ahead scheduling for the hybrid distribution system with CB and ESS considering the wind power prediction curve;
Case 3:
Two-stage stochastic scheduling for the hybrid distribution system with CB considering th wind power scenario set; and
Case 4:
Two-stage stochastic scheduling for the hybrid distribution system with CB and ESS considering the wind power scenario set.
The simulation results for the above cases are listed in Table 6. When the active and reactive power load at all the nodes were increased by 12.5%, the simulation results for the above four cases are listed in Table 7.
From Table 6, it can be seen that: (1) The power loss cost, the power purchasing cost, and power fluctuation penalty costs of the main grid in Case 2 were reduced prominently compared to Case 1, and a similar situation occurs between Case 3 and Case 4, which indicates that considering ESSs can improve the economical performance of the dispatch solution; (2) the various costs in Case 1 are less than the corresponding values in Case 3, and a similar situation exists between Cases 2 and 4. Yet, the deterministic scheduling model is adopted in Cases 1 and 2 while the two-stage stochastic scheduling model is implemented in Cases 3 and 4. Apparently, wind power uncertainty is not taken into consideration in the deterministic scheduling model, which is inconsistent with the actual situation. By contrast, the wind power scenario set has been introduced by the proposed model in Cases 3 and 4 to describe wind power fluctuation in real time. In other words, the influence of wind power uncertainty on the dispatch solution is integrated into the decision process. As a result, the two-stage stochastic scheduling model has increased the operation cost compared with the deterministic scheduling model, but the reliability of the dispatch solution has been improved.
From Table 7, it can be found that load shedding happened in Cases 1 and 3 with the increase of the active and reactive power load, in which ESSs are not integrated in the hybrid AC/DC distribution network. Considering ESSs in Cases 2 and 4 can prominently enhance the economy and reliability of system operation compared with Cases 1 and 3.
Furthermore, the operation status of regulating facilities, such as VSC and ESS, is investigated as follows. For the deterministic scheduling model in Case 2, the transferred active/reactive power by the VSC facilities and the charging/discharging power by ESS are shown in Figure 6 and Figure 7, respectively. It can be seen that: (1) VSC3 operates under the working mode of inverting while VSC1 and VSC2 operate under the working mode of rectifying; (2) as shown in Figure 7, the positive values of power denote that ESS works in the charging state and the negative values mean it works in the discharging state. Combined with the load curves in Figure 5, during the load valley periods of 1 to 8, ESS stored the surplus energy generated by wind generators and the main grid. During the load increasing period of 10 to 20, ESS discharged the previous stored energy to satisfy the balance of supply and demand. By this means, the power fluctuation of the main grid can be alleviated. Meanwhile, ESS can store energy when the electricity price is lower and discharge energy when the electricity price is higher to reduce the power purchasing cost from the main grid.
For the two-stage stochastic scheduling model in Case 4, the transferred active and reactive power by the VSC facilities under different scenarios are shown in Figure 8, and the voltage profiles under different wind power scenarios at hour 17 are shown in Figure 9. From Figure 8, it can be seen that the deviation of transferred power by VSC1 among different wind power scenarios is less than that by VSC2 and VSC3. The reason is that the location nodes of VSC2 and VSC3 are very close to wind generators, the power uncertainty of which has a direct impact on the transmission power of VSC. Combining Figure 8 and Figure 9, it can be found that the node 12 has the highest voltage value for the DC nodes from 8 to 18. The reason is that during the whole scheduling period, SVC2 transferred a large amount of active power to the DC distribution network through node 12, which can be regarded as the DC power supply.
Furthermore, the storage capacity curves by ESS in Cases 2 and 4 when the power load increased by 12.5% are shown in Figure 10. It can be seen that: (1) For the two-stage stochastic scheduling model in Case 4, the charging and discharging periods of ESS are consistent during the whole scheduling period while the charging and discharging power changed with different wind power scenarios; and (2) the charging/discharging period of ESS1 and ESS2 for the deterministic scheduling model in Case 2 is shorter compared to that in Case 4. For instance, ESS1 is not in use at hours 19 to 22 and ESS2 is not in use at hours 16 to 20 in Case 2 while the above ESSs are in the discharging/charging mode at these hours.
In addition, to verify the performance of the proposed model, the charging/discharging state of ESSs obtained by the deterministic scheduling model in Case 2 is implemented for the two-stage stochastic scheduling model in Case 4. The simulated results indicate that the power purchasing cost and power loss cost are similar to the corresponding values while the power fluctuation penalty cost has increased by $36.56 compared to the corresponding value shown in Table 7. It can be drawn that for the same charging/discharging state of ESSs, the two-stage stochastic scheduling model can make more economic decisions to reduce the power fluctuation penalty cost. On the other hand, the scheduling solution by the two-stage stochastic scheduling model has greater decision reliability under the consideration of different wind power scenarios in real time. As a result, the scheduling strategy achieved by the proposed model shows superior performance to ensure the economic and stable operation of the hybrid AC/DC distribution system considering wind power uncertainty.

5. Conclusions

This paper proposed a two-stage stochastic scheduling model for the hybrid AC/DC distribution system considering wind power uncertainty. With the purpose of achieving the active-reactive power coordinated optimal dispatch, the capacitor bank, energy storage system, and voltage source converter were introduced as the active adjustment measures. By using the linearization and second-order cone relaxation techniques for the above distribution network model, the original nonlinear optimization problem could be transformed into a mixed integer second-order cone programming problem to solve. Then, the modified IEEE 33-node hybrid AC/DC distribution system was employed to verify the performance of the proposed model. Th simulation results demonstrated that the two-stage stochastic scheduling model can achieve effective decision-making to ensure the economic and stable operation of the hybrid AC/DC distribution system with wind power integration.

Author Contributions

Conceptualization and supervision, Y.Z.; methodology, W.G.; resources, W.H.; software, Z.Q.; validation, M.Y.; writing—original draft, P.K.; writing—review and editing, F.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by Technology Project of State Grid Energy Saving Company (No. 52789918000F).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Zhao, B.; Zhao, Y.; Wang, Y.; Liu, G.; Song, Q.; Yuan, Z.; Yao, S. Energy internet system based on flexible medium-voltage DC distribution. Proc. CSEE 2015, 35, 4834–4851. [Google Scholar]
  2. Salomonsson, D.; Sannino, A. Low-voltage DC distribution system for commercial power systems with sensitive electronic loads. IEEE Trans. Power Deliv. 2007, 22, 1620–1627. [Google Scholar] [CrossRef]
  3. Starke, M.; Li, F.; Tolbert, L.M.; Ozpineci, B. AC vs. DC distribution: Maximum transfer capability. In Proceedings of the 2008 IEEE Power and Energy Society General Meeting—Conversion and Delivery of Electrical Energy in the 21st Century, 20–24 July 2008; IEEE: Pittsburgh, PA, USA, 2008. [Google Scholar]
  4. Wang, P.; Goel, L.; Liu, X.; Choo, F.H. Harmonizing AC and DC: A hybrid AC/DC future grid solution. IEEE Power Energy Mag. 2013, 11, 76–83. [Google Scholar] [CrossRef]
  5. Eajal, A.A.; Shaaban, M.F.; Ponnambalam, K.; El-Saadany, E.F. Stochastic centralized dispatch scheme for AC/DC hybrid smart distribution systems. IEEE Trans. Sustain. Energy 2017, 7, 1046–1059. [Google Scholar] [CrossRef]
  6. Ahmed, H.M.A.; Eltantawy, A.B.; Salama, M.M.A. A generalized approach to the load flow analysis of AC–DC hybrid distribution systems. IEEE Trans. Power Syst. 2018, 33, 2117–2127. [Google Scholar] [CrossRef]
  7. Chaudhary, S.K.; Guerrero, J.M.; Teodorescu, R. Enhancing the capacity of the AC distribution system using DC interlinks-a step toward future DC grid. IEEE Trans. Smart Grid 2015, 6, 1722–1729. [Google Scholar] [CrossRef]
  8. Xu, Y.; Dong, Z.; Zhang, R.; Hill, D.J. Multi-timescale coordinated voltage/var control of high renewable-penetrated distribution systems. IEEE Trans. Power Syst. 2017, 32, 4398–4408. [Google Scholar] [CrossRef]
  9. Zhang, L.; Tang, W.; Liang, J.; Cong, P.; Cai, Y. Coordinated day-ahead reactive power dispatch in distribution network based on real power forecast errors. IEEE Trans. Power Syst. 2015, 31, 2472–2480. [Google Scholar] [CrossRef]
  10. Zhang, L.; Tang, W.; Liang, J.; Li, G.; Cai, Y.; Yan, T. A medium voltage hybrid AC/DC distribution network and its economic evaluation. In Proceedings of the 12nd IET International Conference on AC and DC Power Transmission, Beijing, China, 28–29 May 2016. [Google Scholar]
  11. Zhang, L.; Chen, Y.; Shen, C.; Tang, W.; Jun, L. Coordinated voltage regulation of hybrid AC/DC medium voltage distribution networks. J. Mod. Power Syst. Clean Energy 2018, 6, 463–472. [Google Scholar] [CrossRef] [Green Version]
  12. Wang, L.; Liang, D.H.; Crossland, A.F.; Taylor, P.C.; Jones, D.; Wade, N.S. Coordination of multiple energy storage units in a low-voltage distribution network. IEEE Trans. Smart Grid 2015, 6, 2906–2918. [Google Scholar] [CrossRef] [Green Version]
  13. Li, X.; Yao, L.; Hui, D. Optimal control and management of a large-scale battery energy storage system to mitigate fluctuation and intermittence of renewable generations. J. Mod. Power Syst. Clean Energy 2016, 4, 593–603. [Google Scholar] [CrossRef] [Green Version]
  14. Mahdavi, S.; Hemmati, R.; Jirdehi, M.A.; Department of Electrical Engineering, Kermanshah University of Technology. Two-level planning for coordination of energy storage systems and wind-solar-diesel units in active distribution networks. Energy 2018, 151, 954–965. [Google Scholar] [CrossRef]
  15. Zheng, Y.; Zhao, J.; Song, Y.; Luo, F.; Meng, K.; Qiu, J.; Hill, D.J. Optimal operation of battery energy storage system considering distribution system uncertainty. IEEE Trans. Sustain. Energy 2018, 9, 1051–1060. [Google Scholar] [CrossRef]
  16. Kazari, H.; Oraee, H.; Pal, B.C. Assessing the Effect of Wind Farm Layout on Energy Storage Requirement for Power Fluctuation Mitigation. IEEE Trans. Sustain. Energy 2019, 10, 558–568. [Google Scholar] [CrossRef]
  17. Li, Y.; Wang, J.; Gu, C.; Liu, J.; Li, Z. Investment optimization of grid-scale energy storage for supporting different wind power utilization levels. J. Mod. Power Syst. Clean Energy 2019, 7, 1721–1734. [Google Scholar] [CrossRef] [Green Version]
  18. Tomczewski, A.; Kasprzyk, L. Optimisation of the structure of a wind farm—Kinetic energy storage for improving the reliability of electricity supplies. Appl. Sci. 2018, 8, 1439. [Google Scholar] [CrossRef] [Green Version]
  19. Baran, M.E.; Wu, F.F. Network reconfiguration in distribution systems for loss reduction and load balancing. IEEE Trans. Power Deliv. 1989, 4, 1401–1407. [Google Scholar] [CrossRef]
  20. Wang, S.; Chen, S.; Xie, S. Security-constrained coordinated economic dispatch of energy storage systems and converter stations for AC/DC distribution networks. Autom. Electr. Power Syst. 2017, 41, 85–91. [Google Scholar]
  21. Ma, X.; Guo, R.; Wang, L.; Huang, J.; Liu, J.; Fang, C. Day-ahead scheduling model for AC/DC active distribution network based on second-order cone programming. Autom. Electr. Power Syst. 2018, 42, 144–150. [Google Scholar]
  22. Farivar, M.; Low, S.H. Branch flow model: Relaxations and convexification—Part I. IEEE Trans. Power Syst. 2013, 28, 2554–2564. [Google Scholar] [CrossRef]
  23. Low, S.H. Convex relaxation of optimal power flow, Part II: Exactness. IEEE Trans. Control. Netw. 2014, 1, 177–189. [Google Scholar] [CrossRef] [Green Version]
  24. Zhang, Y.; Liu, K.; Liao, X.; Qin, L.; An, X. Stochastic dynamic economic emission dispatch with unit commitment problem considering wind power integration. Int. Trans. Electr. Energy 2017, 28, e2472. [Google Scholar] [CrossRef]
  25. Mckay, M.D.; Beckman, R.J.; Conover, W.J. A comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics 1979, 21, 239–245. [Google Scholar]
  26. Harbrecht, H.; Peters, M.; Schneider, R. On the low-rank approximation by the pivoted Cholesky decomposition. Appl. Numer. Math. 2012, 62, 428–440. [Google Scholar] [CrossRef] [Green Version]
  27. Wu, L.; Shahidehpour, M.; Li, T. Stochastic security-constrained unit commitment. IEEE Trans. Power Syst. 2007, 22, 800–811. [Google Scholar] [CrossRef]
  28. Gao, H.; Wang, L.; Liu, J.; Wei, Z. Integrated day-ahead scheduling considering active management in future smart distribution system. IEEE Trans. Power Syst. 2018, 33, 6049–6061. [Google Scholar] [CrossRef]
  29. Wang, J.; Xu, J.; Liao, S.; Sima, L.; Sun, Y.; Wei, C. Coordinated optimization of integrated electricity-gas energy system considering uncertainty of renewable energy output. Autom. Electr. Power Syst. 2019, 43, 2–9. [Google Scholar]
  30. Zhang, Y.; Le, J.; Zheng, F.; Liu, K. Two-stage distributionally robust coordinated scheduling for gas-electricity integrated energy system considering wind power uncertainty and reserve capacity configuration. Renew. Energy 2019, 135, 122–135. [Google Scholar] [CrossRef]
  31. Lofberg, J. YALMIP: A toolbox for modeling and optimization in MATLAB. In Proceedings of the 2004 IEEE International Symposium on Computer Aided Control Systems Design, Taipei, China, 2–4 September 2004; pp. 284–289. [Google Scholar]
Figure 1. The schematic framework of this study.
Figure 1. The schematic framework of this study.
Applsci 10 00181 g001
Figure 2. Schematic diagram of the voltage source converter.
Figure 2. Schematic diagram of the voltage source converter.
Applsci 10 00181 g002
Figure 3. The original and reduced wind power scenario sets.
Figure 3. The original and reduced wind power scenario sets.
Applsci 10 00181 g003
Figure 4. The modified IEEE 33-node hybrid AC/DC distribution network.
Figure 4. The modified IEEE 33-node hybrid AC/DC distribution network.
Applsci 10 00181 g004
Figure 5. Wind power output and load curves.
Figure 5. Wind power output and load curves.
Applsci 10 00181 g005
Figure 6. The transferred active and reactive power by VSC in Case 2.
Figure 6. The transferred active and reactive power by VSC in Case 2.
Applsci 10 00181 g006
Figure 7. The charging/discharging power of energy storage system in Case 2.
Figure 7. The charging/discharging power of energy storage system in Case 2.
Applsci 10 00181 g007
Figure 8. The transferred active and reactive power by voltage source converter in Case 4.
Figure 8. The transferred active and reactive power by voltage source converter in Case 4.
Applsci 10 00181 g008
Figure 9. Voltage profiles at hour 17 in Case 4.
Figure 9. Voltage profiles at hour 17 in Case 4.
Applsci 10 00181 g009
Figure 10. Storage capacity of ESS in Cases 2 and 4.
Figure 10. Storage capacity of ESS in Cases 2 and 4.
Applsci 10 00181 g010
Table 1. Parameters of capacitor bank.
Table 1. Parameters of capacitor bank.
Bus Node Q s t e p C B /kVar N m a x C B y m a x C B
26053
Table 2. Parameters ofenergy storage system.
Table 2. Parameters ofenergy storage system.
Bus Node P m a x E S S /kW E m a x E S S /(kWh) E i n t E S S /%ηchηdis y m a x E S S
14, 283001800300.93811.10663
Table 3. Parameters ofvoltage source converter.
Table 3. Parameters ofvoltage source converter.
Bus Node Q m a x S V C ( A C ) /kVar S m a x S V C ( A C ) /kVAμ
7, 21, 3060015000.866
Table 4. Parameters of gas-fired units.
Table 4. Parameters of gas-fired units.
Bus Node P m i n G /kW P m a x G /kW Q m a x G /kW
6, 2330300150
Table 5. Parameters for cost.
Table 5. Parameters for cost.
ParameterscGcLcWcLosscRamp
values/($/kWh)0.04210.20.0320.04
Table 6. Simulation results for different cases.
Table 6. Simulation results for different cases.
CasePower Purchasing Cost ($)Power Fluctuation Penalty Cost ($)Power Loss (kW)Load Shedding(kW)
Case 12463.86178.81186.470
Case 22343.8357.95181.620
Case 32487.71240.71225.780
Case 42427.2575.05203.120
Table 7. Simulation results for different cases.
Table 7. Simulation results for different cases.
CasePower Purchasing Cost ($)Power Fluctuation Penalty Cost ($)Power Loss (kW)Load Shedding(kW)
Case 13230.61178.56273.3862.23
Case 23118.2752.73268.270
Case 33249.85235.76319.21110.34
Case 43134.9680.01276.650

Share and Cite

MDPI and ACS Style

Kang, P.; Guo, W.; Huang, W.; Qiu, Z.; Yu, M.; Zheng, F.; Zhang, Y. Two-Stage Stochastic Programming Scheduling Model for Hybrid AC/DC Distribution Network Considering Converters and Energy Storage System. Appl. Sci. 2020, 10, 181. https://doi.org/10.3390/app10010181

AMA Style

Kang P, Guo W, Huang W, Qiu Z, Yu M, Zheng F, Zhang Y. Two-Stage Stochastic Programming Scheduling Model for Hybrid AC/DC Distribution Network Considering Converters and Energy Storage System. Applied Sciences. 2020; 10(1):181. https://doi.org/10.3390/app10010181

Chicago/Turabian Style

Kang, Peng, Wei Guo, Weigang Huang, Zejing Qiu, Meng Yu, Feng Zheng, and Yachao Zhang. 2020. "Two-Stage Stochastic Programming Scheduling Model for Hybrid AC/DC Distribution Network Considering Converters and Energy Storage System" Applied Sciences 10, no. 1: 181. https://doi.org/10.3390/app10010181

APA Style

Kang, P., Guo, W., Huang, W., Qiu, Z., Yu, M., Zheng, F., & Zhang, Y. (2020). Two-Stage Stochastic Programming Scheduling Model for Hybrid AC/DC Distribution Network Considering Converters and Energy Storage System. Applied Sciences, 10(1), 181. https://doi.org/10.3390/app10010181

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop