Next Article in Journal
Effect of Crack Orientation on Fatigue Life of Reinforced Concrete Bridge Decks
Next Article in Special Issue
An Intelligent Hybrid Energy Management System for a Smart House Considering Bidirectional Power Flow and Various EV Charging Techniques
Previous Article in Journal
Cyclic p-y Curves of Monopiles in Dense Dry Sand Using Centrifuge Model Tests
Previous Article in Special Issue
Virtual Inertia and Mechanical Power-Based Control Strategy to Provide Stable Grid Operation under High Renewables Penetration
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Parameterized Modeling and Planning of Distributed Energy Storage in Active Distribution Networks

1
Key Laboratory of Smart Grid of Ministry of Education, Tianjin University, Tianjin 300072, China
2
State Grid Economic and Technological Research Institute Co. Ltd., Beijing 102209, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2019, 9(8), 1643; https://doi.org/10.3390/app9081643
Submission received: 14 March 2019 / Revised: 6 April 2019 / Accepted: 15 April 2019 / Published: 20 April 2019
(This article belongs to the Special Issue Smart Home and Energy Management Systems 2019)

Abstract

:
In recent years, distributed energy storage (DES) has experienced rapid growth and has been widely applied in active distribution networks (ADNs). Owing to the close correlation between the characteristics and the application scenarios, DES modeling needs to be parameterized separately for various application demands. In this paper, a parameterized model for optimal DES planning in ADNs is proposed. The typical scenarios for DES planning are generated by the clustering technique, containing the patterns of load demand, wind turbine output and photovoltaic output. Secondly, an optimal planning model of DES considering parameterized characteristics is established, which is essentially a mixed integer non-linear optimization problem. Then, the model is converted to a mixed-integer second-order cone programming model, which can be solved efficiently by available commercial software. Finally, case studies on the modified IEEE 33-node system and IEEE 123-node system verify the efficiency of the proposed method, and the effects of DES planning are validated by two evaluation indexes.

1. Introduction

With the integration of distributed generators (DGs) and flexible loads, traditional distribution systems are evolving into active distribution networks (ADNs) [1]. The integration of renewable DGs contributes to power loss reduction, power supply reliability enhancement, reduction of pollution gases emission, etc. [2,3], while bringing new challenges to the planning and operation of distribution system like voltage violation [4], bi-direction power flow [5], power flow fluctuations [6]. In recent years, distributed energy storage (DES) as an important measure to alleviate the above problems has developed rapidly in both economic and technical aspects [7,8]. Compared with large-scale centralized energy storage plants, DES has fewer restrictions on the geographical conditions of installation location and performs an increasingly important part in various application scenarios of ADNs [9].
By the criterion of discharge duration, the DES technologies can be divided into two categories, energy-type storage and power-type storage. The energy-type storage technologies, including pumped hydraulic storage (PHS), compressed air energy storage (CAES), and large-scale battery storage, are applied to provide long-term electricity support such as arbitrage [10], loss reduction [11], congestion alleviation [12], and long time-scale voltage control [13]. The latter type of energy storage, known as power-type storage, is suitable for distributed energy resources smoothing [14], power quality management [15], frequency regulation [16] and short time-scale voltage control [17]. Supercapacitor or ultracapacitor energy storage (SCES), flywheel energy storage (FES), superconducting magnetic energy storage (SMES) and small-scale battery storage belong to this category.
However, the high costs of DES, especially high investment costs, are the major obstacles limiting further development. Although some new technologies such as lead-carbon battery and LiFePO4 (LFP) battery have largely mitigated the economic problem [18], DES is not as widely used as other traditional electricity equipment. As a result, it is of great significance to optimize the sizing and allocation of DES to maximize the economic benefits in different scenarios.
Many works in the literature have studied on the sizing and allocation of DES in ADNs. A combined genetic algorithm and sequential quadratic programming method is proposed in [19], and further testify the potential economic benefits of energy storage systems. Considering the voltage deviation, congestion, network losses and load shifting, a multi-objective planning model is proposed in [20]. A bi-level programming method is proposed in [21], and optimal planning and operation of energy storage system are solved by upper and lower programming respectively. Considering the influence of soft open point (SOP) and network configuration, a second-order cone programming (SOCP) method to determine the optimal siting and sizing DES is established in [22]. In [23], a relaxed convex optimal power flow (OPF) model is proposed, and solved by benders decomposition method. The optimal methods of resources allocation in the fields of information communication and industrial manufactory also provide ideas for DES planning. For instance, a resource service sharing model of cloud manufactory (CMfg) based on the Gale–Shapley algorithm is proposed in [24], indicating great advantages in promoting the resource-sharing utilization rate. In [25], to support the capacity sharing issue among independent firms, an advanced framework based on game theory and fuzzy logic is proposed.
The type selection of energy storage is a critical issue in DES planning [26]. The characteristics of different DES types vary in economics and technology. Owing to the close correlation between the characteristics and the application scenarios, the DES modeling needs to be parameterized separately for various application demands. For instance, low cycle efficiency will increase the cost of effective output power; and a low cycle life adds to the total cost owing to high-frequency equipment updates. A single type of energy storage device has its unique advantages, but it is difficult to incorporate high efficiency, long life, and many other requirements in one device. Therefore, coordinated planning of multiple types of energy storage systems can adequately cover the potential technical advantages of various DES, and meet the demand for energy storage performance in distribution networks. On the other hand, with the development of DES, the cost of energy storage systems will be significantly reduced, and the economic benefits in distribution networks will become more promising. For the reason that the cost of different DES types varies enormously, it is necessary to compare the costs of different types of energy storage in different scenarios and choose the optimal DES selection to maximize economic benefits. However, the quantization parameters of DES are complicated and challenging to parameterize, and it will increase the difficulty of the calculation due to the parameter multiplication of various DES types. In summary, there is a need to propose a parameterized model, which can transform the efficiency, life cycle, and other characteristics of multiple DES types into parameterized representation. The parameterized model will provide a quantitative description for different typed of DES, making it possible to scientifically and comprehensively consider the specific characteristics of different DES types in optimal planning.
This paper proposes a parameterized model for optimal DES planning in ADNs. The main contributions are summarized as follows:
(1)
The economic and technical characteristics of DES are comprehensively categorized and analyzed. Then, the influences of DES characteristics on the operation and planning are also considered to provide a parameterized modeling for DES planning in ADNs.
(2)
An optimal planning model of DES considering parameterized characteristics of DES is proposed in this paper. The proposed model considers multiple types of DES, which can solve the problem with coordinated planning. By applying the convex relaxation technique and introducing the relevant auxiliary variables, the proposed model is transformed into an effectively solved mixed-integer second-order cone programming (MISOCP) model.
The remainder of this paper is organized as follows. In Section 2, a set of parameterized data is defined to apply in DES planning in ADNs. In Section 3, considering two different application scenarios for DES, the parameterized modeling of DES planning in ADNs is proposed. Section 4 proposes the method that transforms the original model into an effectively solved mixed-integer second-order cone programming (MISOCP) model. In Section 5, numerical results on the modified IEEE 33-node system and modified IEEE 123-node system demonstrate the effectiveness of the proposed method, and the effects of DES are evaluated by two evaluation indexes. Finally, conclusions are drawn in Section 6.

2. Parameterized Representation of Distributed Energy Storage (DES) Characteristics

In general, DES normally exhibits different economic and technical characteristics [27]. The common characteristics of energy storage are the capital cost of DES, lifetime in cycles and years, power/energy density, and maximum depth of discharge (DoD). The investment cost of DES includes the cost of energy storage and power inverters, which is the most important factor affecting the economic benefits of DES. The cycle efficiency indicates the round-trip efficiency in one cycle operation. Low cycle efficiency will increase the cost of effective output power, further increasing the total energy purchase cost from the external grid. Maximum depth of discharge is the ability to discharge the total energy of DES in a discharge duration. High DoD means a higher total energy output at the same energy capacity. Meanwhile, DES with high DoD is more economical on the same capacity of total energy output. The lifetime demonstrates the durability of DES. Low cycle life will result in quickly life decreasing during the frequent switching of discharge-charge state and add to the total cost owing to high-frequency equipment updates. Overall, the characteristics of DES largely affect the results of DES planning. Therefore, all these characteristics of DES need to be taken into account and parameterized while modeling.
In this paper, three representative types of DES are considered, and the parameters are shown in Table 1 [7,10]. The three types of DES have their own features: a lead-acid battery is the most widely used rechargeable battery owing to its low investment cost; a Li-ion battery has the highest efficiency; and a vanadium redox flow battery (VRB) has an exceptionally long lifetime. As a consequence, parametrized modeling is needed to search their applicability in different demand scenarios.

3. Parameterized Modeling of Optimal DES Planning in Active Distribution Networks (ADNs)

3.1. Typical Scenario Generation

By clustering techniques, the method is developed to capture the daily patterns of load demand, wind turbine output and photovoltaic output. With the clustering method, typical scenarios and their occurrence probabilities can be obtained [28]. The purpose of clustering is to group the data objects into multiple clusters in regards to the similarity. The objects within a cluster are similar, whereas the objects of different clusters are dissimilar.
The k-means algorithm is one of the most popular clustering techniques in the fields of mathematical statistics, pattern recognition, machine learning and data mining. It can well reflect the geometric and statistical significance of clustering [29].
The calculation steps of the k-means clustering algorithm are as follows:
Step 1: Given the number of clusters x, randomly select x objects from the cluster data as cluster centers;
Step 2: According to the principle that the distance from the cluster center is the smallest, the remaining objects are assigned to the corresponding categories;
Step 3: After all the objects are divided by category, the cluster center is recalculated; the object with the smallest sum of distances to other objects in the same category is determined as the center of the current class;
Step 4: Repeat Steps (2) and (3) until the cluster center no longer changes, the cluster ends. Then, the clustering result is obtained.
By the above steps, the typical scenarios for DES planning are generated by the clustering technique, containing the patterns of load demand, wind turbine output and photovoltaic output.

3.2. Objective Function

In order to clearly show the feasibility and applicability of the proposed method in different application scenarios, two objective functions are considered and calculated respectively. Moreover, to meet the comprehensive planning requirements in real scenarios, the individual objects can be easily combined and solved by the multi-objective approaches.
In reality, the DGs and DES may have different ownerships, such as belonging to the distribution company, the third-party, and the consumers. The differences in ownerships may affect their observability, controllability, and dispatch cost. For simplicity, all the DGs and DES in this paper are assumed to be owned by the distribution company or the distribution system operator (DSO), which means that they all can be dispatched without extra costs besides the operational cost. The dispatch cost of DGs that belong to the third-party can be easily considered in the proposed method with an additional cost item.

3.2.1. Economic Benefits of ADN

The minimum annual comprehensive cost of ADN is set as the objective function, which is expressed as follows:
min C 1 = C OPE + C INV
where C OPE is the annual operation cost of ADN and C INV is the investment cost of DES, which are formulated as follows:
C OPE = s Ω S ( t N T λ t P s , t SUB Δ t ) p s
C INV = d ( 1 + d ) y m ( 1 + d ) y m 1 m Ω type i N N ( C m POW y i , m S m unit + C m ENE z i , m E m unit )
where Ω S is the set of scenarios and Ω type is the set of DES types. N T is the total periods of the time horizon. N N is the total number of the nodes. C m POW and C m ENE are the capital cost for 1 kW power capacity and 1 kWh energy capacity of the m th DES. S m unit and E m unit are the unit power capacity and energy capacity of the m th DES. y i , m and z i , m are the total number of power units and energy units of the m th DES at node i . P s , t SUB is the active power at the substation in the s th scenario. Δ t is the time interval, λ t is the time-of-use electricity price, p s is the probability of the s th scenario, d is the discount rate of DES, and y m is the lifetime of the m th DES.
The annual operation cost of an ADN is represented by the cost of total energy purchasing from the external grid under the electricity time-of-use tariff. The benefits brought by DES in power loss reduction can be reflected in the cost of energy purchasing from the external grid. The investment cost includes the cost of energy storage and power inverter.
As the power source of the distribution network, the substation node can be conveniently equipped with DES to optimize the power flow between the distribution network and the external bulk power grid. However, in this paper we mainly focus on the optimization within the distribution network by the optimal planning of different DESs. For this application, the effectivity of DES that installed at the source node is limited. Therefore, all the nodes except the substation node are considered as candidate locations for DES planning in this paper.

3.2.2. Power Fluctuation Smoothing

Utilizing the fast charge/discharge characteristics, DES can effectively smooth the rapid power fluctuation brought by renewable DGs comprising wind turbines (WTs) and photovoltaic generators (PVs). Considering the investment cost of DES and the ability to smooth power fluctuation of DGs, the objective function is expressed as follows:
min C 2 = λ f C FLU + C INV
where λ f is the cost of DG power fluctuation. To demonstrate the performance of DES in the application scenario of power fluctuation smoothing, λ f is the transformed parameter indicating the mathematical relationship between the power fluctuation and the investment cost of DES. The purpose of the transformation process is to make the objective function more clear. λ f is designed according to [30]. C FLU is the power fluctuation of DG, which is expressed as:
C FLU = s Ω S { m Ω type i N N t N T [ ( P s , t , i DG + P s , t , i , m DES ) ( P s , t 1 , i DG + P s , t 1 , i , m DES ) ] 2 } p s
where P s , t , i DG is the active power injection by DG at node i in the s th scenario, and P s , t , i , m DES is the active power injection by the m th DES at node i in the s th scenario.
The power fluctuation of DG is determined by the accumulation of power difference between adjacent moments. Considering that the DES integrated with DGs can be best utilized to smooth the power fluctuation, only the nodes with DG installed are considered as candidate locations for DES planning in this application scenario.

3.3. System Power Flow Constraint

i k Ω b P s , t , i k = j i Ω b ( P s , t , j i r j i I s , t , j i 2 ) + P s , t , i
i k Ω b Q s , t , i k = j i Ω b ( Q s , t , j i x j i I s , t , j i 2 ) + Q s , t , i
P s , t , i = P s , t , i DG + m Ω type P s , t , i , m DES P s , t , i LOAD
Q s , t , i = Q s , t , i DG + m Ω type Q s , t , i , m DES Q s , t , i LOAD
P s , t SUB = i j Ω b P s , t , i j ,     i Ω sub
U s , t , i 2 U s , t , j 2 2 ( r i j P s , t , i j + x i j Q s , t , i j ) + ( r i j 2 + x i j 2 ) I s , t , i j 2 = 0
I s , t , i j 2 = P s , t , i j 2 + Q s , t , i j 2 U s , t , i 2
where Ω b is the set of all branches and Ω sub is the set of substation nodes. P s , t , i j and Q s , t , i j are the active and reactive power flow of branch i j in the s th scenario. P s , t , i and Q s , t , i are the total active and reactive power injection at node i in the s th scenario. I t , i j is the branch current magnitude and U t , i is the node voltage magnitude. r i j and x i j are the resistance and reactance of branch i j . P s , t , i LOAD and Q s , t , i LOAD is the active and reactive power load at node i in the s th scenario.
Constraints (6) and (7) represent the power balance of node i in period t at the s th scenario. The power injection of node i in period t in the s th scenario can be described as (8) and (9). Constraint (10) represents the active power balance at the substation node. The voltage drop equation over each branch is expressed as (11). The current magnitude of each branch can be determined by using (12).

3.4. Secure Operation Constraint

( U i min ) 2 U s , t , i 2 ( U i max ) 2
0 I s , t , i j 2 ( I i j max ) 2
where U i min and U i max are the lower and upper limit of statutory voltage at node i . I i j max is the upper limit of statutory current at branch i j . The limits on nodal voltage and current are described as constraints (13) and (14).

3.5. DES Operation Constraint

( P s , t , i , m DES ) 2 + ( Q s , t , i , m DES ) 2 y i , m S m unit
P s , t , i , m DES , L = A m DES ( P s , t , i , m DES ) 2 + ( Q s , t , i , m DES ) 2
E s , t + Δ t , i , m DES = E s , t , i , m DES ( P s , t , i , m DES + P s , t , i , m DES , L ) Δ t
z i , m E m unit S O C m min E s , t , i , m DES z i , m E m unit S O C m max
t N T | P s , t , i , m DES | Δ t 2 z i , m E m unit 1 365 N m M y m
E s , i , m , t = T DES = z i , m E m unit S O C s , m , t = 0
where P s , t , i , m DES and Q s , t , i , m DES are the active and reactive power injection by the m th DES at node i in the s th scenario. E s , t , i , m DES is the energy stored of the m th DES at node i in the s th scenario. A m DES is the loss coefficient of the m th DES. S O C m min and S O C m max are the minimum and maximum state of charge limit of the m th DES. N m M is the cycle life of the m th DES.
Constraint (15) represents the converter capacity limit of DES. The power loss of DES is considered in constraint (16). Constraint (17) determines that the energy stored in DES in period ( t + 1 ) depends on the previous energy stored in DES and the charge/discharge power of the time interval. Constraint (18) represents the maximum and minimum state of charge (SOC) limits. Weighted energy throughput method is adopted to describe the cycle life of DES in constraint (19) [31]. The constraint relies on the principle that an electrochemical cell can exchange a finite amount of charge during its lifespan. This value can be assessed through cycling tests and calculated as the charge that pass through a cell during a complete discharge-charge cycle multiplied by the total number of cycles that a battery can perform before depletion. Constraint (20) implies that the energy stored in DES at the end of one day has to equal its initial value. It is noteworthy that the power of DES is determined by the minimum power of the energy storage equipment and the power converter. To make the DES modeling more generalized and simplified, an ideal power matching between the energy storage device and its power converter is assumed in this paper. Therefore, constraint (15) can be considered as the power limit of the entire energy storage system.

3.6. DES Planning Constraint

m Ω type i N N y i , m S m unit S BGT
m Ω type i N N z i , m E m unit E BGT
y i , m S m unit S BGT δ i
i N N δ i n DES
where S BGT and E BGT are the maximum power capacity and energy capacity of DES in planning. δ i is a binary variable indicating the location of DES at node i . n DES is the maximum number of DES installing nodes.
Constraints (21) and (22) represent the maximum limits on power and energy capacity of DES in planning. The constraint that represents whether DES is installed at node i can be described as (23). When δ i = 1 , DES is installed at node i ; when δ i = 0 , no DES is installed at node i . The limit on the maximum number of nodes that DES can be installed is expressed as constraint (24).
As a consequence, parameterized model for optimal DES planning in ADN is finally expressed as (25), which is a mixed-integer non-linear programming (MINLP) model.
MINLP   Model :   min C 1   or   min C 2 s . t .   ( 6 ) ( 24 )

4. Mixed-Integer Second-Order Cone Programming (MISOCP) Model Conversion

In this section, to solve the MINLP model more effectively, the original model is converted into a MISOCP model by convex relaxation. First, let u s , t , i and i s , t , i j denote the quadratic terms U s , t , i 2 and I s , t , j i 2 Linearized functions are expressed as follows:
i k Ω b P s , t , i k = j i Ω b ( P s , t , j i r j i i s , t , i j ) + P s , t , i
i k Ω b Q s , t , i k = j i Ω b ( Q s , t , j i x j i i s , t , i j ) + Q s , t , i
u s , t , i u s , t , j 2 ( r i j P s , t , i j + x i j Q s , t , i j ) + ( r i j 2 + x i j 2 ) i s , t , i j = 0
i s , t , i j = P s , t , i j 2 + Q s , t , i j 2 u s , t , i
( U i min ) 2 u s , t , i ( U i max ) 2
0 i s , t , i j ( I i j max ) 2
The operation constraints of DES in (15) and (16) can be transformed into the rotated quadratic cone constraints:
( P s , t , i , m DES ) 2 + ( Q s , t , i , m DES ) 2 2 y i , m S m unit 2 y i , m S m unit 2
( P s , t , i , m DES ) 2 + ( Q s , t , i , m DES ) 2 2 P s , t , i , m DES , L 2 A m DES P s , t , i , m DES , L 2 A m DES
Although the absolute value term | P s , t , i , m DES | in constraint (19) is nonlinear, it can be replaced by a linear program according to Chebyshev approximation problem [32]. An auxiliary variable P s , t , i , m AUX is introduced to replace | P s , t , i , m DES | , and constraint (19) can be transformed to constraint (34). Constraints (35)–(37) are added to bound the value of P s , t , i , m AUX , making constraint (34) linearized. Then, the constraint (19) is transformed into constraints (34)–(37), which are expressed as:
t N T P s , t , i , m AUX Δ t 2 z i , m E m unit 1 365 N m M y m
P s , t , i , m AUX P s , t , i , m DES
P s , t , i , m AUX P s , t , i , m DES
P s , t , i , m AUX 0
Constraint (29) is relaxed to the inequality constraint (38), then transformed into second-order cone constraint (39) [33]:
P s , t , i j 2 + Q s , t , i j 2 i s , t , i j u s , t , i
| | [ 2 P s , t , i j 2 Q s , t , i j i s , t , i j u s , t , i ] T | | 2 i s , t , i j + u s , t , i
After the conic relaxation, the original MINLP model is converted into a MISOCP model to realize an efficient calculation. The MISOCP model is expressed as:
MISOCP   Model :   min C 1   or   min C 2 s . t .   ( 8 ) ( 10 ) ,   ( 17 ) ,   ( 18 ) ,   ( 20 ) ( 24 ) ,   ( 26 ) ( 28 ) ,   ( 30 ) ( 37 ) ,   ( 39 )

5. Case Study

In this section, the effectiveness of the proposed method is analyzed and verified on a modified IEEE 33-node system [34]. The voltage level is 12.66 kV, the total active load is 3715 kW, and the total reactive load is 2300 kvar. The test case is shown in Figure 1.
Three types of DES are available to be selected in the planning, and the parameters are shown in Table 1. Five wind turbines and three photovoltaic generators are installed, as shown in Table 2. The investment costs of DES are annualized based on a discount rate of 8%. Time-of-use electricity price is shown in Table 3.
Assuming that the annual load curve, annual WT output curve and PV output curve in the distribution network are shown in Figure 2, Figure 3 and Figure 4, respectively. The typical scenarios for DES planning are generated by the clustering technique, containing the patterns of load demand, wind turbine output and photovoltaic output, which are shown in Figure A1 in the Appendix A. It is assumed that the maximum power and energy capacity of DES in planning are 1 MVA/4 MWh. The total number of nodes that DES can be installed with is up to 4.
The proposed model is implemented with the General Algebraic Modeling System (GAMS) optimization software and is solved by the commercial solver GUROBI [35]. GUROBI is an optimization package which has been widely applied in solving the MISOCP problem. The calculation is carried out on a PC with Intel Core i5 3.20 GHz CPU and 4 GB RAM.

5.1. Planning Results

5.1.1. Economic Benefits Improvement of ADN

The planning results are shown in Table 4, Table 5 and Table 6. It can be seen that different types of DES have different performances. A lead-acid battery and Li-ion battery are both helpful but with different economic benefits in this case; a VRB is not selected to install because of the negative profit. In coordinated planning, the lead-acid battery and Li-ion battery are selected to allocate at the nodes 31 and 32 simultaneously. From Table 6, when only the Li-ion battery is allocated, the operation cost of ADN is lowest. However, the annual comprehensive cost is lowest in the coordinated planning, reducing by $10,378; the annual operation cost of the ADN is reduced by $67,881, a decrease of 5.21%. The price of lead-acid batteries is relatively low, but the cycle efficiency and the life cycle of the Li-ion battery are high. The two DESs have complementary advantages, and the economics of system operation is improved by coordinated planning. Although VRB has the highest cycle life compared to the other two DESs, VRB energy storage is expensive and has not been selected to be allocated in ADNs.

5.1.2. Power Fluctuation Smoothing

The planning results are shown in Table 7 and Table 8. The coordinated planning results are the same as the planning of only VRB in Table 7. It can be seen that the power fluctuation cost is reduced by $49,902 when a VRB is allocated, a decrease of 68.18%. Considering the investment cost of DES, C 2 is reduced by $18,702, a decrease of 25.43%.
In this case, because of the rapid power fluctuation of renewable DGs, DES needs to charge and discharge frequently to meet the demand for power smoothing. VRB has an exceptionally long lifetime compared with the other two DES types, making it is quite suitable for charging/discharging with high frequency. Although a lead-acid battery has an advantage in investment cost, it is not selected to install owing to its low life cycle.
With the development of DES technology, capital cost will be further reduced, and the economic benefits of DES will become more promising. In addition, DES can also perform the auxiliary service functions such as peak shaving, frequency support, and voltage regulation. Considering the above economic and environmental benefits, the practical value of DES will be further improved.

5.2. Evaluation of DES Planning

5.2.1. Economic Evaluation

The rate of return is set as the economic evaluation index, considering the benefit from reducing annual operation cost of ADN and the investment cost of DES, which is expressed as follows:
R R = B C INV C INV × 100 %
where B is the benefit from reducing the operation cost of ADN, which are formulated as follows:
B = C OPE , BEF C OPE , AFT
where C OPE , BEF is the operation cost of ADN before DES installed, and C OPE , AFT is the operation cost after DES installed.
The results of the economic evaluation are shown in Table 9. The results show that a lead-acid battery has better economic benefits than a Li-ion battery. As VRB is not allocated owing to its high investment cost, it has no economic improvement. Furthermore, compared to the planning of single-type DES, the coordinated planning of DES significantly improves the economic benefits of ADNs.

5.2.2. Power Smoothing Evaluation

In the case of power fluctuation smooth, the power fluctuation coefficient is chosen as the evaluation index, which is the ratio of the power after DES smooth to original power output. The power fluctuation coefficient is expressed as:
f FLU = s Ω S { m Ω type i N N t N T [ ( P s , t , i DG + P s , t , i , m DES ) ( P s , t 1 , i DG + P s , t 1 , i , m DES ) ] } p s s Ω S { m Ω type i N N t N T [ P s , t , i DG + P s , t 1 , i DG ] } p s
The evaluation results of power fluctuation smooth in Section 5.1.2 are shown in Table 10. The results show that VRB is more capable in power smooth of DG than Li-ion battery. As lead-acid battery is not allocated owing to its low cycle life, it has no ability for DG power smoothing.

5.3. Scalability Verification

The modified IEEE 123-node system is tested to verify the scalability of the proposed method on large-scale ADNs. The detailed parameters can refer to [36]. Three WTs and six PVs are integrated into the networks, of which the basic installation parameters are shown in Table 11. The test case is shown in Figure 5.
The results of the economic benefits case study are shown in Table 12 and Table 13. It can be seen from the results that a lead-acid battery is selected to allocate in ADN and the economic benefits are improved remarkably.
The results of the power fluctuation smooth by each DES type are shown in Table 14 and Table 15. The coordinated planning results are the same as the planning of VRB in Table 14. It can be seen from the results that the power fluctuation cost is reduced $58,495 by the integration of VRB, a decrease of 49.73%. The Li-ion battery is also applied to power smoothing, but the effect is inferior to VRB.

6. Conclusions

In this paper, a parameterized model for optimal DES planning in ADNs is presented. The typical scenarios for DES planning are generated by a clustering technique, containing the patterns of load demand, wind turbine output and photovoltaic output. To effectively solve the problem, the original MINLP model is transformed into a MISOCP model by convex conversion. The calculation and evaluation results show that through the presented method, the economic benefits of the distribution network have been significantly improved and the power fluctuation of DG has been smoothed through the optimal planning of DES.
With the rapid development of energy storage technology, it is foreseeable that the energy storage cost will be greatly reduced, and the applications of DES will become more and more extensive. The proposed method considers the economic and technical characteristics of DES and can achieve optimal planning with multiple types of DES in ADNs. It can provide reasonable schemes based on the parameterized model and support the distribution network planners to make decisions on the optimal sizing, allocation, and type selection of DES. It also shows great potential in real applications with different temporal and spatial scales, providing more effective support for renewable resource accommodation, voltage and frequency regulation, and self-healing control by detailed type selection and reasonable planning of DES.

Author Contributions

T.Z., H.Y. and G.S. conceived and designed the study; T.Z. and H.Y. performed the study; C.S. and P.L. reviewed and edited the manuscript; T.Z. wrote the paper; H.Y. and G.S. made the critical and final reviewing.

Funding

This research was funded by the National Natural Science Foundation of China, grant number 51807132 and U1866207.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

Sets
Ω b set of all branches
Ω sub set of substation nodes
Ω s set of scenarios
Ω type set of DES types
Indices
i , j indices of nodes
t indices of time periods
s indices of scenarios
m indices of DES types
Variables
P s , t , i j active power flow of branch i j in the s th scenario
Q s , t , i j reactive power flow of branch i j in the s th scenario
P s , t , i total active power injection at node i in the s th scenario
Q s , t , i total reactive power injection at node i in the s th scenario
I t , i j , i t , i j branch current magnitude and its square
U t , i , u t , i node voltage magnitude and its square
P s , t , i DG active power injection by DG at node i in the s th scenario
Q s , t , i DG reactive power injection by DG at node i in the s th scenario
P s , t SUB active power at the substation in the s th scenario
P s , t , i , m DES active power injection by the m th DES at node i in the s th scenario
Q s , t , i , m DES reactive power injection by the m th DES at node i in the s th scenario
P s , t , i , m DES , L active power losses of the m th DES at node i in the s th scenario
E s , t , i , m DES energy stored of the m th DES at node i in the s th scenario
y i , m total number of power unit of the m th DES at node i
z i , m total number of energy unit of the m th DES at node i
Parameters
N T total periods of the time horizon
N N total number of the nodes
P s , t , i LOAD active power load at node i in the s th scenario
Q s , t , i LOAD reactive power load at node i in the s th scenario
r i j resistance of branch i j
x i j reactance of branch i j
U i max upper limit of statutory voltage at node i
U i min lower limit of statutory voltage at node i
I i j max upper limit of statutory current at branch i j
S m unit unit power capacity of the m th DES
E m unit unit energy capacity of the m th DES
C m POW capital cost for 1 kW power capacity of the m th DES
C m ENE capital cost for 1 kWh energy capacity of the m th DES
A m DES loss coefficient of the m th DES
S O C m max maximum state of charge limit of the m th DES
S O C m min minimum state of charge limit of the m th DES
λ t time-of-use electricity price
λ f cost of DG power fluctuation
N m M life cycle (cycles) of the m th DES
y m lifetime (years) of the m th DES
d discount rate of DES
D o D m max maximum depth of discharge of the m th DES
Δ t time interval
p s probability of the s th scenario
S BGT maximum power capacity of DES in planning
E BGT maximum energy capacity of DES in planning
δ i binary variable of DES installing at node i
n DES maximum number of DES installing nodes

Appendix A

Figure A1. Typical scenarios generation for DES planning.
Figure A1. Typical scenarios generation for DES planning.
Applsci 09 01643 g0a1

References

  1. Wang, C.S.; Song, G.Y.; Li, P.; Ji, H.R.; Zhao, J.L.; Wu, J.Z. Optimal siting and sizing of soft open points in active electrical distribution networks. Appl. Energy 2017, 189, 301–309. [Google Scholar] [CrossRef]
  2. Rueda-Medina, A.C.; Padilha-Feltrin, A. Distributed Generators as Providers of Reactive Power Support—A Market Approach. IEEE Trans. Power Syst. 2013, 28, 490–502. [Google Scholar] [CrossRef]
  3. Sedghi, M.; Ahmadian, A.; Aliakbar-Golkar, M. Optimal storage planning in active distribution network considering uncertainty of wind power distributed generation. IEEE Trans. Power Syst. 2016, 31, 304–316. [Google Scholar] [CrossRef]
  4. Li, P.; Ji, H.R.; Wang, C.S.; Zhao, J.L.; Song, G.Y.; Ding, F.; Wu, J.Z. Coordinated Control Method of Voltage and Reactive Power for Active Distribution Networks Based on Soft Open Point. IEEE Trans. Sustain. Energy 2017, 8, 1430–1442. [Google Scholar] [CrossRef]
  5. Li, P.; Ji, H.R.; Wang, C.S.; Zhao, J.L.; Song, G.Y.; Ding, F.; Wu, J.Z. Optimal operation of soft open points in active distribution networks under three-phase unbalanced conditions. IEEE Trans. Smart. Grid 2019, 10, 380–391. [Google Scholar] [CrossRef]
  6. Li, C.; Liu, X.B.; Cao, Y.J.; Zhang, P.; Shi, H.Q.; Ren, L.Y.; Kuang, Y.H. A time-scale adaptive dispatch method for renewable energy power supply systems on islands. IEEE Trans. Smart. Grid 2016, 7, 1069–1078. [Google Scholar] [CrossRef]
  7. Saboori, H.; Hemmati, R.; Ghiasi, S.M.S.; Dehghan, S. Energy storage planning in electric power distribution networks–A state-of-the-art review. Renew. Sustain. Energy Rev. 2017, 79, 1108–1121. [Google Scholar] [CrossRef]
  8. Ji, H.R.; Wang, C.S.; Li, P.; Ding, F.; Wu, J.Z. Robust operation of soft open points in active distribution networks with high penetration of photovoltaic integration. IEEE Trans. Sustain. Energy 2019, 10, 280–289. [Google Scholar] [CrossRef]
  9. Shen, X.W.; Shahidehpour, M.; Han, Y.D.; Zhu, S.Z.; Zheng, J.H. Expansion Planning of Active Distribution Networks With Centralized and Distributed Energy Storage Systems. IEEE Trans. Sustain. Energy 2016, 8, 126–134. [Google Scholar] [CrossRef]
  10. Huang, S.J.; Wu, Q.W.; Oren, S.S.; Li, R.Y.; Liu, Z.X. Distribution locational marginal pricing through quadratic programming for congestion management in distribution networks. IEEE Trans. Power Syst. 2015, 30, 2170–2178. [Google Scholar] [CrossRef]
  11. Ji, H.R.; Wang, C.S.; Li, P.; Song, G.Y.; Yu, H.; Wu, J.Z. Quantified analysis method for operational flexibility of active distribution networks with high penetration of distributed generators. Appl. Energy 2019, 239, 706–714. [Google Scholar] [CrossRef]
  12. Poullikkas, A. A comparative overview of large-scale battery systems for electricity storage. Renew. Sustain. Energy Rev. 2013, 27, 778–788. [Google Scholar] [CrossRef]
  13. Chen, Q.X.; Zhao, X.Y.; Gan, D.H. Active-reactive scheduling of active distribution system considering interactive load and battery storage. Proc. CSEE 2017, 2, 29. [Google Scholar] [CrossRef]
  14. Lin, W.; Jin, X.D.; Mu, Y.F.; Jia, H.J.; Xu, X.D.; Yu, X.D.; Zhang, B. A two-stage multi-objective scheduling method for integrated community energy system. Appl. Energy 2018, 216, 428–441. [Google Scholar] [CrossRef]
  15. Hung, D.Q.; Mithulananthan, N.; Bansal, R.C. Integration of PV and BES units in commercial distribution systems considering energy loss and voltage stability. Appl. Energy 2014, 113, 1162–1170. [Google Scholar] [CrossRef]
  16. Wan, C.; Zhao, J.; Song, Y.H.; Xu, Z.; Lin, J.; Hu, Z.C. Photovoltaic and solar power forecasting for smart grid energy management. CSEE J. Power Energy Syst. 2016, 1, 38–46. [Google Scholar] [CrossRef]
  17. Dong, X.H.; Mu, Y.F.; Xu, X.D.; Jia, H.J.; Wu, J.Z.; Yu, X.D.; Yan, Q. A charging pricing strategy of electric vehicle fast charging stations for the voltage control of electricity distribution networks. Appl. Energy 2018, 225, 857–868. [Google Scholar] [CrossRef]
  18. Aneke, M.; Wang, M.H. Energy storage technologies and real life applications – A state of the art review. Appl. Energy 2016, 179, 350–377. [Google Scholar] [CrossRef]
  19. Carpinelli, G.; Celli, G.; Mocci, S.; Mottola, M.; Pilo, F.; Proto, D. Optimal integration of distributed energy storage devices in smart grids. IEEE Trans. Smart. Grid 2013, 4, 985–995. [Google Scholar] [CrossRef]
  20. Nick, M.; Cherkaoui, R.; Paolone, M. Optimal allocation of dispersed energy storage systems in active distribution networks for energy balance and grid support. IEEE Trans. Power Syst. 2014, 29, 2300–2310. [Google Scholar] [CrossRef]
  21. Xiao, J.; Zhang, Z.Q.; Bai, L.Q. Determination of the optimal installation site and capacity of battery energy storage system in distribution network integrated with distributed generation. IET Gener. Trans. Distrib. 2016, 10, 601–607. [Google Scholar] [CrossRef]
  22. Bai, L.Q.; Jiang, T.; Li, F.X.; Chen, H.H.; Li, X. Distributed energy storage planning in soft open point based active distribution networks incorporating network reconfiguration and DG reactive power capability. Appl. Energy 2018, 210, 1082–1091. [Google Scholar] [CrossRef]
  23. Nick, M.; Cherkaoui, R.; Paolone, M. Optimal planning of distributed energy storage systems in active distribution networks embedding grid reconfiguration. IEEE Trans. Power Syst. 2018, 33, 1577–1590. [Google Scholar] [CrossRef]
  24. Liu, Y.K.; Zhang, L.; Tao, F.; Wang, L. Resource service sharing in cloud manufacturing based on the Gale–Shapley algorithm: Advantages and challenge. Int. J. Comput. Integr. Manuf. 2017, 30, 420–432. [Google Scholar] [CrossRef]
  25. Argoneto, P.; Renna, P. Supporting capacity sharing in the cloud manufacturing environment based on game theory and fuzzy logic. Enterp. Inf. Syst. 2016, 10, 193–210. [Google Scholar] [CrossRef]
  26. Zidar, M.; Georgilakis, P.S.; Hatziargyriou, N.D.; Capuder, T.; Škrlec, D. Review of energy storage allocation in power distribution networks: applications, methods and future research. IET Gener. Trans. Distrib. 2015, 10, 645–652. [Google Scholar] [CrossRef]
  27. Luo, X.; Wang, J.H.; Dooner, M.; Clarke, J. Overview of current development in electrical energy storage technologies and the application potential in power system operation. Appl. Energy 2015, 137, 511–536. [Google Scholar] [CrossRef] [Green Version]
  28. Wang, C.S.; Jiao, B.Q.; Guo, L.; Yuan, K.; Sun, B. Optimal planning of stand-alone microgrids incorporating reliability. J. Modern Power Syst. Clean Energy 2014, 2, 195–205. [Google Scholar] [CrossRef] [Green Version]
  29. Kanungo, T.; Mount, D.M.; Netanyahu, N.S.; Piatko, C.D.; Silverman, R.; Wu, A.Y. An efficient k-means clustering algorithm: analysis and implementation. IEEE Trans. Pattern Anal. Mach. Intell. 2002, 24, 881–892. [Google Scholar] [CrossRef]
  30. Wang, S.X.; Wang, K.; Teng, F.; Strbac, G.; Wu, L. An affine arithmetic-based multi-objective optimization method for energy storage systems operating in active distribution networks with uncertainties. Appl. Energy 2018, 223, 215–228. [Google Scholar] [CrossRef]
  31. Sauer, D.U.; Wenzl, H. Comparison of different approaches for lifetime prediction of electrochemical systems—Using lead-acid batteries as example. J. Power Sour. 2008, 176, 534–546. [Google Scholar] [CrossRef]
  32. Boyd, S.; Vandenberghe, L. Convex Optimization; Cambridge University Press: Cambridge, UK, 2004. [Google Scholar]
  33. Farivar, M.; Low, S.H. Branch flow model: Relaxations and convexification (Part I). IEEE Trans. Power Syst. 2013, 28, 2554–2564. [Google Scholar] [CrossRef]
  34. Baran, M.E.; Wu, F.F. Network reconfiguration in distribution systems for loss reduction and load balancing. IEEE Trans. Power Del. 1989, 4, 1401–1407. [Google Scholar] [CrossRef]
  35. Brook, A.; Kendrick, D.; Meeraus, A. GAMS, a user’s guide. ACM Signum Newsl. 1988, 23, 10–11. [Google Scholar] [CrossRef]
  36. Chen, X.; Wu, W.C.; Zhang, B.M. Robust restoration method for active distribution networks. IEEE Trans. Power Syst. 2016, 31, 4005–4015. [Google Scholar] [CrossRef]
Figure 1. Structure of the modified IEEE 33-node system.
Figure 1. Structure of the modified IEEE 33-node system.
Applsci 09 01643 g001
Figure 2. Annual load curve.
Figure 2. Annual load curve.
Applsci 09 01643 g002
Figure 3. Annual power output of wind turbines (WTs).
Figure 3. Annual power output of wind turbines (WTs).
Applsci 09 01643 g003
Figure 4. Annual power output of photovoltaic generators (PVs).
Figure 4. Annual power output of photovoltaic generators (PVs).
Applsci 09 01643 g004
Figure 5. Structure of the modified IEEE 123-node system.
Figure 5. Structure of the modified IEEE 123-node system.
Applsci 09 01643 g005
Table 1. Typical parameters of distributed energy storage (DES) [7,10].
Table 1. Typical parameters of distributed energy storage (DES) [7,10].
ParameterLead-Acid BatteryLi-Ion BatteryVRB
Capital cost of power converter ($/kW)505050
Capital cost of energy storage ($/kWh)125200250
Cycle efficiency (%)909575
Cycle life (cycles)3000500010,000
Lifetime (years)101215
Maximum depth of discharge (%)709070
Table 2. Parameters of distributed generators (DGs).
Table 2. Parameters of distributed generators (DGs).
ParameterWind TurbinesPhotovoltaic Generators
Location101617303371327
Capacity (kVA)500300200200300500300400
Table 3. Parameters of time-of-use electricity price.
Table 3. Parameters of time-of-use electricity price.
ParameterOn-PeakMid-PeakOff-Peak
Time span16:00–22:008:00–15:001:00–7:00, 23:00–24:00
Electricity price ($/kWh)0.1730.1040.050
Table 4. Planning results of economic benefits improvement in IEEE 33-node system.
Table 4. Planning results of economic benefits improvement in IEEE 33-node system.
TypeLocationPower Capacity (kVA)Energy Capacity (kWh)
Only lead-acid battery14120610
17100450
303901750
32210930
Only Li-ion battery17120990
301601420
3180690
33100900
Vanadium redox flow battery (VRB)
Table 5. Coordinated planning results of economic benefits improvement in IEEE 33-node system.
Table 5. Coordinated planning results of economic benefits improvement in IEEE 33-node system.
TypeLocationPower Capacity (kVA)Energy Capacity (kWh)
Lead-acid battery1450550
1770810
3120260
3260690
Li-ion battery311301040
3280650
VRB
Total14, 17, 31, 324104000
Table 6. Annual costs of DES planning in IEEE 33-node system.
Table 6. Annual costs of DES planning in IEEE 33-node system.
ParameterOperation Cost of ADN ($)Investment Cost of DES ($) C 1   ( $ )
Without DES1,303,6631,303,663
Only lead-acid battery allocated1,253,52145,8461,299,367
Only Li-ion battery allocated1,226,92471,1721,298,096
Coordinated planning1,235,78257,5031,293,285
Table 7. Planning results of power fluctuation smooth in IEEE 33-node system.
Table 7. Planning results of power fluctuation smooth in IEEE 33-node system.
TypeLocationPower Capacity (kVA)Energy Capacity (kWh)
Only Li-ion battery102601380
161601240
17180540
33400840
Only VRB102901450
162001240
17160490
33350820
Only lead-acid battery
Table 8. Power smooth results of DES planning in IEEE 33-node system.
Table 8. Power smooth results of DES planning in IEEE 33-node system.
Power Fluctuation Cost ($)Investment Cost of DES ($) C 2   ( $ )
Without DES73,19173,191
Only Li-ion battery allocated24,69244,12568,817
Only VRB allocated23,28931,20054,489
Table 9. Results of the economic evaluation in IEEE 33-node system.
Table 9. Results of the economic evaluation in IEEE 33-node system.
ParameterLead-Acid BatteryLi-Ion BatteryVRBCoordinated Planning
R R 9.378.70018.05
Table 10. Power smoothing evaluation results in IEEE 33-node system.
Table 10. Power smoothing evaluation results in IEEE 33-node system.
ParameterLead-Acid BatteryLi-Ion BatteryVRB
f FLU 1.000.130.11
Table 11. Parameters of DGs.
Table 11. Parameters of DGs.
ParameterWind TurbinesPhotovoltaic Generators
Location289210833428697111116
Capacity (kVA)80012001000300300200200400600
Table 12. Planning results of economic benefits improvement in IEEE 123-node system.
Table 12. Planning results of economic benefits improvement in IEEE 123-node system.
TypeLocationPower Capacity (kVA)Energy Capacity (kWh)
Only lead-acid battery3270300
471101220
5880760
1191501720
Only Li-ion battery
Only VRB
Table 13. Annual costs of DES planning in IEEE 123-node system.
Table 13. Annual costs of DES planning in IEEE 123-node system.
Operation Cost of ADN ($)Investment Cost of DES ($) C 1   ( $ )
Without DES2,087,8632,087,863
Only lead-acid battery allocated2,032,36945,8462,078,215
Table 14. Planning results of power fluctuation smoothing in IEEE 123-node system.
Table 14. Planning results of power fluctuation smoothing in IEEE 123-node system.
TypeLocationPower Capacity (kVA)Energy Capacity (kWh)
Only Li-ion battery923501350
1083101280
111150630
116190740
Only VRB923801390
1083301320
111120610
116170680
Only lead-acid battery
Table 15. Power smoothing results of DES planning in IEEE 33-node system.
Table 15. Power smoothing results of DES planning in IEEE 33-node system.
Power Fluctuation Cost ($)Investment Cost of DES ($) C 2   ( $ )
Without DES117,620117,620
Only Li-ion battery allocated69,44344,125113,568
Only VRB allocated59,12531,20090,325

Share and Cite

MDPI and ACS Style

Zhao, T.; Yu, H.; Song, G.; Sun, C.; Li, P. Parameterized Modeling and Planning of Distributed Energy Storage in Active Distribution Networks. Appl. Sci. 2019, 9, 1643. https://doi.org/10.3390/app9081643

AMA Style

Zhao T, Yu H, Song G, Sun C, Li P. Parameterized Modeling and Planning of Distributed Energy Storage in Active Distribution Networks. Applied Sciences. 2019; 9(8):1643. https://doi.org/10.3390/app9081643

Chicago/Turabian Style

Zhao, Tianyu, Hao Yu, Guanyu Song, Chongbo Sun, and Peng Li. 2019. "Parameterized Modeling and Planning of Distributed Energy Storage in Active Distribution Networks" Applied Sciences 9, no. 8: 1643. https://doi.org/10.3390/app9081643

APA Style

Zhao, T., Yu, H., Song, G., Sun, C., & Li, P. (2019). Parameterized Modeling and Planning of Distributed Energy Storage in Active Distribution Networks. Applied Sciences, 9(8), 1643. https://doi.org/10.3390/app9081643

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