Next Article in Journal
PRE-YOLO: A Lightweight Model for Detecting Helmet-Wearing of Electric Vehicle Riders on Complex Traffic Roads
Previous Article in Journal
Jaw Morphology and Factors Associated with Upper Impacted Canines: Case-Controlled Trial
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Distributed Cooperative Dispatch Method of Distribution Network with District Heat Network and Battery Energy Storage System Considering Flexible Regulation Capability

by
Xin Fu
1,
Shunjiang Yu
2,*,
Qibo He
1,
Long Wang
1,
Changming Chen
2,
Chengeng Niu
2 and
Zhenzhi Lin
2
1
State Grid Wuxi Power Supply Company of Jiangsu Electric Power Co., Ltd., Wuxi 214061, China
2
College of Electrical Engineering, Zhejiang University, Hangzhou 310027, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2024, 14(17), 7699; https://doi.org/10.3390/app14177699
Submission received: 12 August 2024 / Revised: 27 August 2024 / Accepted: 29 August 2024 / Published: 31 August 2024

Abstract

:
Flexible resources, including district heat networks (DHN) and battery energy storage systems (BESS), can provide flexible regulation capability for distribution networks (DN), thereby increasing the absorption capacity for renewable energy. In order to improve the operation economy of DN and ensure the information privacy of different operators, a distributed cooperative dispatch method of DN with DHN and BESS considering flexible regulation capability is proposed. First, a distributed cooperative dispatch framework of DN-DHN-BESS is constructed. Then, an optimal dispatch model of DHN under constant flow-variable temperature control strategy is established in order to utilize the heat storage capacity to provide flexible regulation capability for DN. Next, the optimal dispatch models of BESS and DN are established with the objective of minimizing the operation cost, respectively. Finally, a solution method based on the alternating direction multiplier method of distributed cooperative dispatch for DN-DHN-BESS is proposed. Case studies are performed on a system consisting of a 33-node DN and a 44-node DHN, and simulation results demonstrate that the proposed method differs from the centralized dispatch method by only 0.52% in the total system cost, and the proposed method reduces the total system cost by 34.5% compared to that of the independent dispatch method.

1. Introduction

With the increasing penetration of renewable energy sources, the need for flexibility in distribution networks (DN) over multiple time scales has increased [1]. In this context, flexible resources such as battery energy storage systems (BESS) and district heat networks (DHN) have become an important instrument for smoothing out power fluctuations of renewable energy sources and improving power quality due to their ability to transfer energy over time [2]. Therefore, it is important to study the cooperative operation method of DN with different types of flexible resources in order to improve the flexible regulation capability of the system and thus improve the operation economy of the system [3].
With the characteristics of fast power regulation and storage of electricity, battery energy storage systems play a large role in smoothing power fluctuations of intermittent energy sources, shaving peaks and filling valleys, improving voltage quality, and providing backup power sources, etc. It is the key to realizing flexible regulation of widely accessed distributed energy sources in DN. There have been more references on the cooperative dispatch of DN and BESS. In [4], an optimal dispatch method of BESS in unbalanced distribution networks is proposed so as to minimize power losses and simultaneously ensure that the network satisfies current and voltage constraints. In [5], an optimal dispatch method for BESS in DN considering peak-load transfer is proposed to improve the voltage distribution in DN. In order to mitigate the adverse effects of renewable energy sources in active distribution networks, energy management and the optimal dispatch method of BESS have been proposed in [6] to reduce operation costs, power losses, and voltage fluctuations in unbalanced distribution networks. In [7], a joint day-ahead and intraday optimal dispatch method for active distribution networks is proposed, which is capable of dispatching centralized and distributed energy storage based on their differences in capacity and responsiveness.
DHN is a potential flexible resource that can break the rigid constraints of real-time balance between heat supply and heat load by adjusting the water supply temperature [8,9,10]. This type of energy storage effect provides additional flexible regulation capability to DN to consume renewable energy. Several references have studied the cooperative dispatch of DN and DHN. In order to improve the absorption capacity of DN for wind power, an integrated electric-heat dispatch model considering the heat inertia of the DHN has been developed in [11], in which an integrated model for fully simulating the dynamic distribution of the district heat network has been proposed for the first time. In [12], a centralized optimal dispatch model for an integrated electric-heat energy system is developed by considering the use of electric boilers and heat storage tanks to improve the flexibility of combined heating and power (CHP) units for better consumption of wind power. In [13,14], pipelines in the DHN and buildings are further utilized as heat storage to provide flexibility in order to solve the problem of wind curtailment in the DN, and a dispatch model for CHP that takes into account the heat storage capacity of pipelines and buildings is developed. In [15], a coordinated operation strategy based on model predictive control for a multi-area integrated electric-heat energy system is proposed, which takes into account the heat storage capacity of DHN by taking into account the transmission delay of the pipelines and the heat inertia of the buildings.
All of the above references adopt the centralized cooperative dispatch methods, whereas the DN, DHN, and BESS are managed by the distribution network operator, the district heat network operator, and the energy storage aggregator, respectively, and the centralized cooperative dispatch of them is difficult to achieve in this context. On the one hand, centralized optimization requires getting system-wide data, which makes the size of the optimization model increase dramatically and leads to a large computational cost. On the other hand, the operators cannot share all the information in the dispatch process due to the requirement of information privacy. Therefore, how to realize distributed cooperative dispatch among DN, DHN, and BESS needs to be solved urgently.
The augmented Lagrangian relaxation (ALR) method is a more mature distributed algorithm with current applications, which decouples the multi-region optimization problem by introducing quadratic term relaxation coupling constraints in the objective function and iteratively updating the Lagrangian multipliers. Distributed algorithms based on ALR can be categorized into the fully distributed algorithms represented by the auxiliary problem principle [16], the alternating direction multiplier method (ADMM) [17], and the hierarchical distributed algorithms represented by analytical target cascading (ATC) [18]. Among the above distributed methods, since ADMM is a combination of the dyadic decomposition method and the ALR method, ADMM has the robustness of the multiplier method and the distributed computational ability of the dyadic decomposition and is also widely used in power system control and optimization for its excellent convergence, fast iteration speed, and protection of the subject’s privacy.
References [19,20,21] have investigated ADMM-based dispatch methods for distribution networks. In [19], a coordinated day-ahead dispatch method of multiple power distribution grids hosting stochastic resources based on ADMM is proposed to aggregate the power and energy flexibilities in an interconnected power distribution system. In [20], a fully decentralized hierarchical transactive energy framework for charging electric vehicles (EV), with local distributed energy resources (DER) in a power distribution system based on ADMM is used to ensure that the privacy of market participants is well preserved since the bid data of each participant are not exposed to others. In [21], a decentralized distributed convex optimal power flow model for power distribution systems based on ADMM is proposed, which is based on dividing the power grid network into subproblems representing individual areas by interchanging minimum network information. Currently, only a few references have investigated the distributed cooperative dispatch of distribution networks with district heat networks. In [22], a transactive energy-supported economic operation method for multi-energy complementary microgrids (MECM) is proposed to coordinate interconnected MECM in a regional integrated energy system (RIES). In [23], a decentralized demand management method based on the alternating direction method of multipliers algorithm for an industrial park with CHP units and thermal storage, which can protect private data of all participants while achieving solutions with high quality. However, the above references do not establish the accurate model of DHN that includes rate and temperature of mass flow and therefore do not allow the flexible regulation capability of DHN to be utilized. Moreover, the above references do not consider the distributed cooperative dispatch among the three main actors: DN, DHN, and BESS. In conclusion, there is still a lack of research on distributed cooperative dispatch for DN, DHN, and BESS.
In response to the above problems, a distributed cooperative dispatch method for DN with DHN and BESS that considers the flexible regulation capability is proposed. The main contributions are as follows:
(1)
A distributed cooperative dispatch framework for DN-DHN-BESS is constructed, which transforms the cooperative dispatch problem into the distributed optimization problem by decoupling the system into three subsystems at the boundary power.
(2)
An optimal dispatch model for DHN under constant flow-variable temperature (CF-VT) control strategy is established, which takes into account the heat storage capacity by inscribing the temperature quasi-dynamics in pipeline and node.
(3)
An optimal dispatch model for DN is established, which is transformed into a second-order cone programming (SOCP) model based on the second-order cone relaxation and linearization method for quadratic constraint to improve the solution efficiency.
(4)
A solution method based on ADMM of distributed cooperative dispatch for DN-DHN-BESS is proposed, to realize the independent solution and cooperative dispatch of three subsystems.
The rest of the paper is organized as follows: In Section 2, the distributed cooperative dispatch framework for DN-DHN-BESS is constructed. In Section 3, the optimal dispatch model for DHN under CF-VT control strategy is established. In Section 4, the optimal dispatch model for BESS is established. In Section 5, the optimal dispatch model for DN is established. In Section 6, the solution method based on ADMM of distributed cooperative dispatch for DN-DHN-BESS is proposed. In Section 7, case studies are conducted to demonstrate the effectiveness and efficiency of the proposed method. Conclusions are given in Section 8.

2. Distributed Cooperative Dispatch Framework for DN-DHN-BESS

The structure of the DN-DHN-BESS system studied in this paper is shown in Figure 1. The DN includes distributed photovoltaic (PV), distributed wind turbine (WT), and load and is dispatched by the distribution network operator. The DHN includes CHP units, gas boilers (GB), electric heaters (EH), and loads and is dispatched by the district heat network operator. The charging and discharging of the BESS are dispatched by the energy storage aggregator. The difference in boundary power is to distinguish the boundary power information exchanged by different subjects; for example, PD-E is used to represent the boundary power information exchanged by the distribution network, and PE-D represents the boundary power information exchanged by the energy storage.
In order to follow the autonomous operation characteristics and the information privacy of different operators, a distributed cooperative dispatch framework for DN-DHN-BESS is constructed in this paper, as shown in Figure 2. The coupling variable between the DN, DHN, and BESS is the boundary power, so under the premise of satisfying the boundary power consistency constraint in (1), the system can be decoupled into three subsystems, and thus the cooperative dispatch problem can be transformed into a distributed optimization problem. By independently optimizing the three subsystems and exchanging the information of boundary power for coordination, accurate results can be obtained.
{ P D - E = P E - D P D - H = P H - D

3. Optimal Dispatch Model for DHN under CF-VT Control Strategy

3.1. Objective of the Optimal Dispatch Model

The objective of the optimal dispatch model for DHN is to minimize the operation cost CDHN, including the cost of purchased electricity C e DHN , the cost of purchased gas C g DHN , and the cost of pollutant emission penalties C f DHN , as shown in (2).
min C DHN = C e DHN + C g DHN + C f DHN
C e DHN = t i N D P i , t H - D c t E , 2
C g DHN = t m N H ( G m , t CHP + G m , t GB ) c t G
C f DHN = t m N H k Ω P ( P m , t CHP λ k CHP + H m , t GB λ k GB ) c k env
where ND is the set of nodes in DN; NH is the set of nodes in DHN; ΩP is the set of pollutants; P i , t H - D is the power transmitted from the DN node i to the DHN at period t; P m , t CHP and G m , t CHP are the active power and gas consumption power of the CHP unit at DHN node m at period t, respectively; G m , t GB and H m , t GB are the gas consumption power and heat power of the GB at DHN node m at period t, respectively; c t E , 2 is the price of electricity sold from the DN at period t; c t G is the gas price at period t; λ k CHP and λ k GB are the emission factors of category k pollutant of CHP and GB, respectively; c k env is the treatment cost factor of category k pollutant.

3.2. Constraints of the Optimal Dispatch Model

In order to take into account, the heat storage capacity, the temperature quasi-dynamics in pipeline and node are considered in the power flow constraints of DHN, and pipeline temperature decay and delay constraints and node temperature mixing constraints, are established, respectively [24]. Further combined with the heat power balance constraint, operation constraints of the energy device, and boundary power constraints between DHN and DN, the complete constraints are formed.

3.2.1. Pipeline Temperature Decay and Delay Constraints

Considering the transmission delay and heat loss of mass flow in the pipeline, the pipeline temperature decay and delay constraints in the supply and return networks are established as shown in (6) and (7), respectively [25]:
T m n , t S , tail = T t am + ( T m n , t τ m n S , head T t am ) e λ m n L m n f m n S A m n ρ w c w ,   m n L H
T m n , t R , head = T t am + ( T m n , t τ m n R , tail T t am ) e λ m n L m n f m n R A m n ρ w c w ,   m n L H
τ m n = l m n A m n ρ w f m n S ,   m n L H
where LH is the set of pipelines in DHN; T m n , t S , head and T m n , t S , tail are the temperature at the head and end of pipeline mn in supply network at period t, respectively; T m n , t R , head and T m n , t R , tail are the temperature at the head and end of pipeline mn in return network at period t; T t am is the outdoor temperature at period t; f m n S and f m n R are the mass flow rate of pipeline mn in supply and return network, respectively; τmn is the transmission delay of pipeline mn; λmn is the heat conductivity factor of pipeline mn; lmn is the length of pipeline mn; Amn is the cross-sectional area of pipeline mn; cw is the specific heat capacity of water; ρw is the density of water.

3.2.2. Node Temperature Mixing Constraints

Considering that mass flow with different temperatures will mix at the nodes, the node temperature mixing constraints in the water supply and return networks are established as shown in (9) and (10) respectively.
T m , t S ( l m L H f l m S + f m ) = f m T m , t S , ex + l m L H f l m S T l m , t S , tail ,   m N H
T m , t R ( m n L H f m n R f m ) = f m T m , t R , ex + m n L H f m n R T m n , t R , head ,   m N H
where T m , t S and T m , t R are the temperature at node m in supply and return networks at period t, respectively; T m , t S , ex and T m , t R , ex are the supply and return temperatures of heat exchange station at node m at period t, respectively.
In addition, the temperature of the mass flow entering the pipeline from the nodes should be the same as the node temperature, thus the constraints shown in (11) and (12) are established. Moreover, there is a relation between the node temperature and the supply and return temperatures of the heat exchange station as (13) and (14) [26].
T m , t S = T m n , t S , head ,   m N H ,   m n L H
T m , t R = T l m , t R , tail ,   m N H ,   m n L H
T m , t S = T m , t S , ex ,   m N H
T m , t R = T m , t R , ex ,   m N S H
where N S H is the set of source nodes in DHN.

3.2.3. Heat Power Balance Constraints

H m , t N = c w f m ( T m , t S , ex T m , t R , ex ) ,   m N H
H m , t N = H m , t CHP + H m , t GB + H m , t EH H m , t load ,   m N H
where H m , t N is the heat power of node m at period t; H m , t CHP is the heat power of CHP unit at node m at period t; H m , t EH is the heat power of EH at node m at period t; H m , t load is the heat power of the load at node m at period t [27].

3.2.4. Operation Constraints of Energy Devices

The operational constraints of energy equipment include the operational constraints of CHP units, GB, and EH [28,29].
H m , t CHP = K e CHP P m , t CHP ,   m N H
H m , t CHP = K h CHP G m , t CHP ,   m N H
P min CHP S m CHP P m , t CHP P max CHP S m CHP ,   m N H
where K e CHP is the heat-to-electric ratio coefficient of the CHP unit; K h CHP is the gas-heat conversion efficiency of the CHP unit; S m CHP is the capacity of the CHP unit at node m; P max CHP and P min CHP are the upper and lower limits of the generation power of the CHP unit, respectively.
H m , t GB = K GB G m , t GB ,   m N H
H min GB S m GB H m , t GB H max GB S m GB ,   m N H
where KGB is the heat-gas conversion factor of the GB; S m GB is the capacity of GB at node m; H max GB and H min GB are the upper and lower limits of the heat power of the GB, respectively.
H m , t EH = K EH P m , t EH ,   m N H
H min EH S m EH H m , t EH H max EH S m EH ,   m N H
where P m , t EH is the active power of EH at node m at period t; KEH is the heat-electric conversion factor of the EH; S m EH is the capacity of the EH at node m; H max EH and H min EH are the upper and lower limits of the heat power of EH, respectively.

3.2.5. Boundary Power Constraints between DHN and DN

In order to facilitate distributed cooperative dispatch, the boundary power between the DHN and the DN is calculated as follows:
P i , t H - D = m N i H P m , t EH m N i H P m , t CHP ,   i N D
where N i H is the set of DHN nodes connected to DN node i.

4. Optimal Dispatch Model for BESS

4.1. Objective of the Optimal Dispatch Model

The objective of the optimal dispatch model for BESS is to minimize the operation cost of CBESS, as shown in the following:
min C BESS = t k N E ( P k , t ES , c P k , t ES , f ) c t E , 2
where NE is the set of BESSs; P k , t ES , c and P k , t ES , f are the charging power and discharging power of BESS k at period t, respectively.

4.2. Constraints of the Optimal Dispatch Model

4.2.1. Operation Constraints of BESS

Referring to [30], the operation constraints of BESS are established as follows. Equation (26) is the mutual exclusion constraint of charging and discharging states. Equations (27) and (28) are the charging and discharging power constraints. Equation (29) is the total power constraint. Equations (30) and (31) are the energy constraints.
x k , t ES , c + x k , t ES , f 1 ,   k N E
0 P k , t ES , c x k , t ES , c S k ES ,   k N E
0 P k , t ES , f x k , t ES , f S k ES ,   k N E
P k , t ES = P k , t ES , c P k , t ES , f ,   k N E
E k , t ES = E k , t 1 ES + η ES , c P k , t ES , c Δ t P k , t ES , f Δ t η ES , f ,   k N E
0 E k , t ES E k ES ,   k N E
where x k , t ES , c and x k , t ES , f are the charging and discharging states of the BESS k at period t, respectively; P k , t ES is the net power of the BESS k at period t; S k ES and E k ES are the power and capacity of the BESS k, respectively; E k , t ES is the capacity of the BESS k at period t; ηES,c and ηES,f are the charging and discharging efficiencies of the BESS, respectively.

4.2.2. Boundary Power Constraints between BESS and DN

In order to facilitate distributed cooperative dispatch, the boundary power between the BESS and the DN is calculated as follows:
P i , t E - D = k N i E P k , t ES ,   i N D
where N i E is the set of BESSs connected to DN node i.

5. Optimal Dispatch Model for DN

5.1. Objective of the Optimal Dispatch Model

The objective of the DN dispatch model is to minimize the operation cost CDN, including the cost of purchased electricity C b D N , the penalty cost for curtailment for PV and WT C f D N , and the revenue from the sale of electricity C s D N , as shown in the following:
min C DN = C b DN + C f DN C s DN
C b DN = t i N D P i , t SG c t E , 1
C f DN = t i N D [ y f PV ( P ˜ i , t PV P i , t PV ) + y f WT ( P ˜ i , t WT P i , t WT ) ]
C s DN = t i N D c t E , 2 ( P i , t D - H + P i , t D - E )
where P i , t SG is the active power transmitted from the superior grid to the DN node i at period t; c t E , 1 is the price of electricity sold by the superior grid at period t; P ˜ i , t PV and P i , t PV are the predicted and actual output of PV at DN node i at period t, respectively; P ˜ i , t WT and P i , t WT are the predicted and actual output of WT at DN node i at period t, respectively; y f PV and y f WT are the penalty cost coefficients of curtailed power for PV and WT, respectively; P i , t D - H is the active power transmitted from DN node i to DHN at period t; P i , t D - E is the active power transmitted from DN node i to BESS at period t.

5.2. Constraints of the Optimal Dispatch Model

The DN optimal dispatch model contains the power flow constraint and output constraint of PV and WT.

5.2.1. Power Flow Constraint of DN

Based on the Distflow model, the power flow constraints of DN are established as follows [31].
{ P i , t N = P i , t SG + P i , t PV + P i , t WT P i , t load P i , t D - E P i , t D - H Q i , t N = Q i , t SG + Q i , t PV + Q i , t WT Q i , t load ,   i N D
{ k i L D ( P k i , t L R k i I k i , t 2 ) + P i , t N = i j L D P i j , t L k i L D ( Q k i , t L X k i I k i , t 2 ) + Q i , t N = i j L D Q i j , t L ,   i N D
U i , t 2 U j , t 2 I i j , t 2 ( R i j 2 + X i j 2 ) 2 ( P i j , t L R i j + Q i j , t L X i j ) = 0 ,   i N D ,   i j L D
U i , t 2 = ( P i j , t L ) 2 + ( Q i j , t L ) 2 I i j , t 2 ,   i N D ,   i j L D
( P i j , t L ) 2 + ( Q i j , t L ) 2 ( S i j ) 2 ,   i j L D
( U min ) 2 U i , t 2 ( U max ) 2 ,   i N D
where LD is the set of lines in DN; P i , t N and Q i , t N are the active and reactive power at node i at period t, respectively; Q i , t SG is the reactive power purchased by node i from the superior grid at period t; P i , t PV and Q i , t PV are the active and reactive outputs of the PV at node i at period t, respectively; P i , t WT and Q i , t WT are the active and reactive outputs of the WT at node i at period t, respectively; P i , t load and Q i , t load are the active and reactive power of the load at node i at period t, respectively; P i j , t L and Q i j , t L are the active and reactive power transmitted by line ij at period t, respectively; Rij and Xij are the resistance and reactance of line ij, respectively; U i , t 2 is the squared value of the voltage at node i at period t; I i j , t 2 is the squared value of the current of line ij at period t; Sij is the capacity of line ij; Umin and Umax are the lower and upper limits of voltage, respectively.

5.2.2. Output Constraints of PV and WT

The output constraints of PV and WT are established as follows, where Equation (43) is the output constraint of PV and Equation (44) is the output constraint of WT.
0 P i , t PV P ˜ i , t PV ,   i N D
0 P i , t WT P ˜ i , t WT ,   i N D

5.2.3. Boundary Power Constraints between DN, DHN, and BESS

When DHN or BESS is connected to node i of DN, the boundary power is not 0, otherwise the boundary power is 0. Therefore, the constraints shown in (45) and (46) are established.
{ P i , t D - H = 0 if   N i HN = P i , t D - H 0 otherwise ,   i N D
{ P i , t D - E = 0 if   N i ES = P i , t D - E 0 otherwise ,   i N D

5.2.4. Second-Order Cone Relaxation and Linearization of Constraints

Equation (40) contains the product term of voltage and current. In order to improve the efficiency of the solution, (40) is first relaxed to the following form [32]:
( P i j , t L ) 2 + ( Q i j , t L ) 2 U i , t 2 I i j , t 2 ,   i N D ,   i j L D
Further, (47) is transformed into a second order cone constraint as shown in (48):
[ 2 P i j , t L 2 Q i j , t L I i j , t 2 U i , t 2 ] T I i j , t 2 + U i , t 2 ,   i N D ,   i j L D
In addition, (41) contains the squared terms of active and reactive power; in order to improve the solution efficiency, (41) is linearized based on the quadratic constraint linearization method, and two square constraints (49) are employed to substitute for (41) [32].
{ S i j P i j , t L S i j S i j Q i j , t L S i j 2 S i j P i j , t L + Q i j , t L 2 S i j 2 S i j P i j , t L Q i j , t L 2 S i j ,   i j L D

6. Solution Method of Distributed Cooperative Dispatch for DN-DHN-BESS

6.1. Modification of the Objective Function

Under the framework of ADMM, the boundary power consistency constraint needs to be relaxed and added as a penalty term to the objective functions of the optimal dispatch models of DN, DHN, and BESS so that the boundary power is as close as possible during the calculation process and finally achieves consistency [33]. The modified objective functions of the optimal dispatch models of DN, DHN, and BESS are shown in (50)–(52):
min F DN = C DN + t i { λ i , t D - H ( P i , t D - H P ¯ i , t D - H ) + [ ρ i , t D - H ( P i , t D - H P ¯ i , t D - H ) ] 2 } + t i { λ i , t D - E ( P i , t D - E P ¯ i , t D - E ) + [ ρ i , t D - E ( P i , t D - E P ¯ i , t D - E ) ] 2 }
min F HN = C HN + t i { λ i , t D - H ( P ¯ i , t D - H P i , t H - D ) + [ ρ i , t D - H ( P ¯ i , t D - H P i , t H - D ) ] 2 }
min F ES = C ES + t i { λ i , t D - E ( P ¯ i , t D - E P i , t E - D ) + [ ρ i , t D - E ( P ¯ i , t D - E P i , t E - D ) ] 2 }
where λ i , t D - H and λ i , t D - E are the Lagrange multipliers; ρ i , t D - H and ρ i , t D - E are the penalty parameters; P ¯ i , t D - H and P ¯ i , t D - E are the average value of boundary power, which can be calculated by (53).
{ P ¯ i , t D - H = 1 2 ( P i , t D - H + P i , t H - D ) P ¯ i , t D - E = 1 2 ( P i , t D - E + P i , t E - D )

6.2. Solution Process

The specific solution process includes the following five steps:
Step 1: Initialize the Lagrange multipliers, penalty parameters, and boundary power, and set the iteration flag r = 0.
Step 2: Let r = r +1 and calculate the average value of the boundary power according to (53).
Step 3: The optimal dispatch models for DN, DHN, and BESS are solved separately to obtain the dispatch results and update the boundary power.
Step 4: Calculate the original residuals ( M i , t D - H , ( r ) , M i , t D - E , ( r ) ) and pairwise residuals ( D i , t D - H , ( r ) , D i , t D - E , ( r ) ) according to (54) and (55), and check whether the residuals satisfy the convergence condition in (56). If satisfied, output the calculation results; otherwise, continue with step 5.
{ M i , t D - I , ( r ) = P i , t D - H , ( r ) P ¯ i , t D - H , ( r ) 2 M i , t D - E , ( r ) = P i , t D - E , ( r ) P ¯ i , t D - E , ( r ) 2
{ D i , t D - H , ( r ) = P ¯ i , t D - H , ( r ) P ¯ i , t D - H , ( r 1 ) 2 D i , t D - E , ( r ) = P ¯ i , t D - E , ( r ) P ¯ i , t D - E , ( r 1 ) 2
{ M max ( r ) = max { M i , t D - H , ( r ) , M i , t D - E , ( r ) } ε M D max ( r ) = max { D i , t D - H , ( r ) , D i , t D - E , ( r ) } ε D
where the superscript (r) is the value of the variable corresponding to the r-th iteration; εM and εD are the convergence margins of the original and pairwise residuals, respectively.
Step 5: Update the Lagrange multipliers and penalty parameters according to (57)–(60) and go to step 2.
ρ i , t D - H , ( r + 1 ) = { ξ 1 ρ i , t D - H , ( r ) max { M i , t D - H , ( r ) } > ν max { D i , t D - H , ( r ) } ρ i , t D - H , ( r ) ξ 2 max { D i , t D - H , ( r ) } > ν max { M i , t D - H , ( r ) } ρ i , t D - H , ( r ) other
λ i , t D - H , ( r + 1 ) = λ i , t D - H , ( r ) + ρ i , t D - H , ( r + 1 ) ( P i , t D - H , ( r ) P ¯ i , t D - H , ( r ) )
ρ i , t D - E , ( r + 1 ) = { ξ 1 ρ i , t D - E , ( k ) max { M i , t D - E , ( r ) } > ν max { D i , t D - E , ( r ) } ρ i , t D - E , ( r ) ξ 2 max { D i , t D - E , ( r ) } > ν max { M i , t D - E , ( r ) } ρ i , t D - E , ( r ) other
λ i , t D - E , ( r + 1 ) = λ i , t D - E , ( r ) + ρ i , t D - E , ( r + 1 ) ( P i , t D - E , ( r ) P ¯ i , t D - E , ( r ) )
where ξ1 and ξ2 are the update coefficients of penalty parameters; υ is the update judgment coefficient.

7. Case Studies

The case studies are completed on a computer with an Intel i7-14700K CPU and 32 GB of RAM. The programs are developed using Matlab R2023a and are solved by Gurobi 11.0.0. The relative gap of the MILP solver is 0.1%.

7.1. Test System

The proposed method is tested using the topology shown in Figure 3, containing a 33-node DN and a 44-node DHN [34]. In the DN, nodes 13 and 28 are connected to PV with parameters of 1600 kW and 600 kW, nodes 6 and 24 are connected to WT with parameters of 600 kW and 300 kW, and nodes 6, 24, and 28 are connected to BESS with parameters of 150 kW/300 kWh, 100 kW/200 kWh, and 200 kW/400 kWh, respectively. In the DHN, node 1 is a heat station equipped with 1000 kW CHP, 2000 kW GB, and 2000 kW EH, where CHP and EH are connected to the DN node 13.

7.2. Analysis of the Effectiveness of the Proposed Distributed Cooperative Dispatch Method

In order to verify the effectiveness of the proposed distributed cooperative dispatch method, the centralized dispatch method is used as a comparison [34]. Table 1 shows the comparison of the system costs of the two methods. Table 2 shows the comparison of the calculation times of the two methods.
The comparison in Table 1 shows that the operation costs of DN, DHN, and BESS of the two methods have a very small difference, and the proposed method differs from the centralized dispatch method by only 0.52% in the total system cost, which proves that the proposed distributed cooperative dispatch method does not fall into the local optimum during the solution process and has the same solution accuracy as the centralized dispatch method. The comparison in Table 2 shows that although the solution time of the proposed distributed cooperative dispatch method is much more than that of the centralized dispatch method, it is still within the acceptable range.
In summary, it can be concluded that the proposed distributed cooperative dispatch method can achieve the same solution accuracy as the centralized dispatch method.

7.3. Analysis of the Efficiency of the Proposed Distributed Cooperative Dispatch Method

In order to verify the efficiency of the proposed distributed cooperative dispatch method, the independent dispatch method is used as a comparison. The independent dispatch method indicates that the optimal dispatch models for DHN and BESS are first solved to obtain the boundary power, and then the optimal dispatch model for DN is solved to obtain the dispatch scheme.
Table 3 shows the comparison of the system costs of the two methods. Table 4 shows the comparison of curtailed energy for PV and WT of the two methods. Figure 4 shows the comparison of curtailed power for PV and WT of the two methods. Figure 5 shows the comparison of the power purchased from the superior grid by DN of the two methods.
As shown in Table 3, compared to the independent dispatch method, the operation cost of DN is reduced by 61.2%, the operation cost of DHN is increased by 22.5%, the profitability of BESS is reduced by 65.8%, and the total cost of the system is reduced by 34.5% of the proposed method. This is because the proposed method achieves efficient synergy among DN, DHN, and BESS, which cannot be achieved by an independent dispatch method. As shown in Table 4, compared with the independent dispatch method, the curtailed energy for PV of the proposed method is reduced by 7393.2 kWh. This is because in 7:00–18:00, the curtailed power for PV and WT of the proposed method is significantly smaller than that of the independent dispatch method, as shown in Figure 4. Further, the penalty cost for curtailment for PV and WT of the proposed method is reduced by 86.2% compared with the independent dispatch method. Due to the higher capacity of PV and WT power consumption of the proposed method, less power is purchased from the superior grid compared to the independent dispatch method, as shown in Figure 5. As a result, the cost of purchased electricity from the superior grid of the proposed method is reduced by 7.0% compared to the independent dispatch method. Ultimately, the total system cost of the proposed method is reduced by 34.5% compared to the independent dispatch method.
In summary, compared with the independent dispatch method, the proposed distributed cooperative dispatch method realizes the efficient synergy between DN and the two types of flexible resources, DHN and BESS, which can effectively improve the absorption capacity of DN for PV and WT and then reduce the operation cost.

7.4. Analysis of the Impact of Heat Storage Capacity on Cooperative Dispatch

In order to analyze the impact of heat storage capacity on cooperative dispatch, the following case is set up for comparison with the proposed method.
Case 1-1: adopt the exiting distributed cooperative dispatch method in [22], i.e., the heat storage capacity of DHN is not taken into account.
Table 5 shows the comparison of the system cost of the proposed method and Case 1-1. Table 6 shows the comparison of curtailed energy for PV and WT of the proposed method and Case 1-1. Figure 6 shows the comparison of the water supply temperature of DHN of the proposed method and Case 1-1. Figure 7 shows the comparison of the heat generation power of DHN for the proposed method and Case 1-1. Figure 8 shows the comparison of curtailed power for PV and WT of the proposed method and Case 1-1.
As seen in Table 5, compared to Case 1-1, the operation cost of DN of the proposed method is reduced by 18.1%, the operation cost of DHN is reduced by 2.0%, and the profitability of BESS is reduced by 52.5%, which ultimately results in a reduction of the total system cost by 10.2%. This is because the proposed method takes into account the heat storage capacity of DHN, which significantly increases the operational flexibility of the system. As shown in Figure 6 and Figure 7, under the proposed method, DHN can flexibly adjust the water supply temperature to change the heat generation power. In the period when the output of PV and WT is large, DHN raises the water supply temperature to increase the heat generation power, which converts the electric energy of PV and WT that cannot be absorbed by DN into the heat energy of water and increases the heat storage capacity. In the period when the output of PV and WT is small, DHN reduces the heat generation power and utilizes the heat storage to provide the auxiliary power support, which achieves the effect of increasing the absorption capacity of DN for PV and WT. As shown in Figure 8, the curtailed power for PV and WT of the proposed method is significantly smaller than that of Case 1-1 in several periods. So as shown in Table 6, the curtailed energy for PV of the proposed method is reduced by 1164.3 kWh compared to Case 1-1, and consequently the penalty cost for curtailment for PV and WT of the proposed method is reduced by 49.6% compared with that of Case 1-1.
From the above analysis, it can be concluded that the proposed method considering the heat storage capacity of DHN can provide additional regulation capacity for DN by flexibly adjusting the water supply temperature of DHN, which in turn improves the absorption capacity of DN for PV and WT and then reduces the system operation cost.

7.5. Analysis of the Impact of BESS Capacity on Cooperative Dispatch

In order to analyze the impact of BESS capacity on cooperative dispatch, the following four cases are set up for comparison.
Case 2-1: Reduce the capacity of the energy storage in Section 7.1 by 50%.
Case 2-2: The capacity of the energy storage in Section 7.1 is unchanged.
Case 2-3: Increase the capacity of the energy storage in Section 7.1 by 50%.
Case 2-4: Increase the capacity of the energy storage in Section 7.1 by 100%.
Table 7 shows the comparison of system costs for the four cases. Table 8 shows the comparison of curtailed energy for PV and WT of the four cases. Figure 9 shows the trend in total system cost and curtailed energy for PV and WT of the four cases.
As seen in Table 7, compared to Case 2-2, the operation cost of DN in Case 2-1 is increased by 3.9%, the operation cost of DHN is increased by 3.7%, and the profitability of BESS is reduced by 17.1%, which ultimately results in an increase of the total system cost by 3.9%. This is because when the BESS capacity is reduced, the regulation capacity it can provide to the DN is weakened, and therefore the absorption capacity of DN for PV and WT is reduced. As shown in Table 8, the curtailed energy for PV and WT of Case 2-1 is increased by 3.2% compared to that of Case 2-2. As seen in Table 7, compared to Case 2-2, the operation costs of DN of Case 2-3 and Case 2-4 are reduced by 0.32% and 0.42%, respectively; the operation costs of DHN of Case 2-3 and Case 2-4 are reduced by 0.73% and 0.80%, respectively; and the profitability of BESS of Case 2-3 and Case 2-4 is increased by 95.2% and 106.7%, respectively, which ultimately results in a reduction of the total system cost by 0.95% and 1.1%, respectively. This is because when the BESS capacity is increased, the regulation capacity it can provide to the DN is increased, and therefore the absorption capacity of DN for PV and WT is increased. As shown in Table 8, the curtailed energy for PV and WT of Cases 2-3 and 2-4 is reduced by 3.3% and 4.2% compared to that of Case 2-2.
In addition, as seen in Figure 9, with the increase of BESS capacity, the reduction of the total system cost and the curtailed energy for PV and WT gradually become slower. Therefore, increasing the BESS capacity within a certain range can reduce the total system cost and improve the dispatch economy.

8. Conclusions

In this paper, a distributed cooperative dispatch method for DN with DHN and BESS considering flexible regulation capability is proposed. According to the simulation results, the following conclusions are drawn:
(1) The proposed distributed cooperative dispatch method does not fall into the local optimum during the solution process and differs from the centralized dispatch method by only 0.52% in the total system cost, which indicates that it can achieve the same solution accuracy as the centralized dispatch method.
(2) Compared with the independent dispatch method, the proposed distributed cooperative dispatch method realizes the efficient synergy between DN and the two types of flexible resources, DHN and BESS, which can effectively improve the absorption capacity of DN for PV and WT and then reduce the operation cost of the system by 35.6%.
(3) The proposed method, considering the heat storage capacity of DHN, can provide additional regulation capacity for DN by flexibly adjusting the water supply temperature of DHN, which in turn improves the absorption capacity of DN for PV and WT and then reduces the system operation cost.
In the future, the authors will further consider flexible resources such as electric vehicles in distributed cooperative dispatch to provide more regulation capability for the distribution network and further investigate how to take into account the impact of uncertainty of renewable energy in the distribution network in cooperative dispatch. In addition, the authors will apply the proposed method to a specific case in Wuxi City in the future to further test the performance of the proposed methods in practical applications.

Author Contributions

Conceptualization, X.F., S.Y. and Z.L.; data curation, C.C. and C.N.; formal analysis, Q.H. and L.W.; funding acquisition, X.F., Q.H. and L.W.; investigation, X.F. and S.Y.; methodology, X.F., S.Y. and Q.H.; project administration, X.F. and Z.L.; resources, Q.H. and Z.L.; software, S.Y., L.W. and C.N.; supervision, Z.L.; validation, S.Y., C.C. and Z.L.; visualization, C.C. and C.N.; writing—original draft, X.F., S.Y. and Q.H.; writing—review and editing, S.Y., C.C. and Z.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Science and Technology Project of State Grid Jiangsu Electric Power Co., Ltd. (No. J2023083).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author(s).

Conflicts of Interest

Authors Xin Fu, Qibo He and Long Wang were employed by the company State Grid Wuxi Power Supply Company of Jiangsu Electric Power Co., Ltd. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. He, X.; Qiu, R.C.; Ai, Q.; Zhu, T. A Hybrid Framework for Topology Identification of Distribution Grid with Renewables Integration. IEEE Trans. Power Syst. 2021, 36, 1493–1503. [Google Scholar] [CrossRef]
  2. Showers, S.O.; Raji, A.K. Benefits and Challenges of Energy Storage Technologies in High Penetration Renewable Energy Power Systems. In Proceedings of the 2019 IEEE PES/IAS Power Africa, Abuja, Nigeria, 20 August 2019; pp. 209–214. [Google Scholar]
  3. Calero, F.; Cañizares, C.A.; Bhattacharya, K.; Anierobi, C.; Calero, I.; De Souza, M.F.Z.; Farrokhabadi, M.; Guzman, N.S.; Mendieta, W.; Peralta, D.; et al. A Review of Modeling and Applications of Energy Storage Systems in Power Grids. Proc. IEEE 2023, 111, 806–831. [Google Scholar] [CrossRef]
  4. Watson, J.D.; Watson, N.R.; Lestas, I. Optimized Dispatch of Energy Storage Systems in Unbalanced Distribution Networks. IEEE Trans. Sustain. Energy 2018, 9, 639–650. [Google Scholar] [CrossRef]
  5. Li, X.; Ma, R.; Gan, W.; Yan, S. Optimal Dispatch for Battery Energy Storage Station in Distribution Network Considering Voltage Distribution Improvement and Peak Load Shifting. J. Mod. Power Syst. Clean Energy 2022, 10, 131–139. [Google Scholar] [CrossRef]
  6. Gangwar, T.; Padhy, N.P.; Jena, P. Energy Management Approach to Battery Energy Storage in Unbalanced Distribution Networks. IEEE Trans. Ind. Appl. 2024, 60, 1345–1356. [Google Scholar] [CrossRef]
  7. Wang, G.; Wang, C.; Feng, T.; Wang, K.; Yao, W.; Zhang, Z. Day-Ahead and Intraday Joint Optimal Dispatch in Active Distribution Network Considering Centralized and Distributed Energy Storage Coordination. IEEE Trans. Ind. Appl. 2024, 60, 4832–4842. [Google Scholar] [CrossRef]
  8. Tian, X.; Lin, X.; Zhong, W.; Zhou, Y.; Cong, F. Optimal Dispatch of Integrated Electricity and Heating Systems Considering the Quality-Quantity Regulation of Heating Systems to Promote Renewable Energy Consumption. Energy 2024, 300, 131599. [Google Scholar] [CrossRef]
  9. Li, Z.; Wu, W.; Shahidehpour, M.; Wang, J.; Zhang, B. Combined Heat and Power Dispatch Considering Pipeline Energy Storage of District Heating Network. IEEE Trans. Sustain. Energy 2016, 7, 12–22. [Google Scholar] [CrossRef]
  10. Li, Z.; Wu, W.; Wang, J.; Zhang, B.; Zheng, T. Transmission-Constrained Unit Commitment Considering Combined Electricity and District Heating Networks. IEEE Trans. Sustain. Energy 2016, 7, 480–492. [Google Scholar] [CrossRef]
  11. Zheng, J.; Zhou, Z.; Zhao, J.; Wang, J. Integrated Heat and Power Dispatch Truly Utilizing Thermal Inertia of District Heating Network for Wind Power Integration. Appl. Energy 2018, 211, 865–874. [Google Scholar] [CrossRef]
  12. Chen, X.; Kang, C.; O’Malley, M.; Xia, Q.; Bai, J.; Liu, C.; Sun, R.; Wang, W.; Li, H. Increasing the Flexibility of Combined Heat and Power for Wind Power Integration in China: Modeling and Implications. IEEE Trans. Power Syst. 2015, 30, 1848–1857. [Google Scholar] [CrossRef]
  13. Dai, Y.; Chen, L.; Min, Y.; Chen, Q.; Hao, J.; Hu, K.; Xu, F. Dispatch Model for CHP With Pipeline and Building Thermal Energy Storage Considering Heat Transfer Process. IEEE Trans. Sustain. Energy 2019, 10, 192–203. [Google Scholar] [CrossRef]
  14. Gu, W.; Wang, J.; Lu, S.; Luo, Z.; Wu, C. Optimal Operation for Integrated Energy System Considering Thermal Inertia of District Heating Network and Buildings. Appl. Energy 2017, 199, 234–246. [Google Scholar] [CrossRef]
  15. Wu, C.; Gu, W.; Jiang, P.; Li, Z.; Cai, H.; Li, B. Combined Economic Dispatch Considering the Time-Delay of District Heating Network and Multi-Regional Indoor Temperature Control. IEEE Trans. Sustain. Energy 2018, 9, 118–127. [Google Scholar] [CrossRef]
  16. Khanabadi, M.; Fu, Y.; Gong, L. A Fully Parallel Stochastic Multiarea Power System Operation Considering Large-Scale Wind Power Integration. IEEE Trans. Sustain. Energy 2018, 9, 138–147. [Google Scholar] [CrossRef]
  17. Chen, R.; Sun, Y.; Xu, J.; Liao, S.; Wang, B.; Zhang, Z.; Ding, X.; Wu, Y. Fully Distributed Energy and Reserve Scheduling for Hybrid AC/DC Systems. IEEE Trans. Ind. Appl. 2023, 59, 2173–2185. [Google Scholar] [CrossRef]
  18. Huang, H.; Zhou, M.; Li, G. Coordinated Decentralized Dispatch of Wind-Power-Integrated Multi-Area Interconnected Power Systems Considering Multiple Uncertainties and Mutual Reserve Support. Power Syst. Technol. 2019, 43, 381–389. [Google Scholar]
  19. Gupta, R.; Fahmy, S.; Paolone, M. Coordinated Day-Ahead Dispatch of Multiple Power Distribution Grids Hosting Stochastic Resources: An ADMM-Based Framework. Electr. Power Syst. Res. 2022, 212, 108555. [Google Scholar] [CrossRef]
  20. Yang, J.; Wiedmann, T.; Luo, F.; Yan, G.; Wen, F.; Broadbent, G.H. A Fully Decentralized Hierarchical Transactive Energy Framework for Charging EVs With Local DERs in Power Distribution Systems. IEEE Trans. Transp. Electrif. 2022, 8, 3041–3055. [Google Scholar] [CrossRef]
  21. Biswas, B.D.; Kamalasadan, S. Decentralized Distributed Convex Optimal Power Flow Model for Power Distribution System Based on Alternating Direction Method of Multipliers. IEEE Trans. Ind. Appl. 2023, 59, 627–640. [Google Scholar] [CrossRef]
  22. Yang, Z.; Hu, J.; Ai, X.; Wu, J.; Yang, G. Transactive Energy Supported Economic Operation for Multi-Energy Complementary Microgrids. IEEE Trans. Smart Grid 2021, 12, 4–17. [Google Scholar] [CrossRef]
  23. Wei, J.; Zhang, Y.; Wang, J.; Wu, L.; Zhao, P.; Jiang, Z. Decentralized Demand Management Based on Alternating Direction Method of Multipliers Algorithm for Industrial Park with CHP Units and Thermal Storage. J. Mod. Power Syst. Clean Energy 2022, 10, 120–130. [Google Scholar] [CrossRef]
  24. Chen, C.; Zhang, T.; Wang, Y.; Sun, G.; Wu, C.; Lin, Z. Coordinated Siting and Sizing for Integrated Energy System Considering Generalized Energy Storage. Int. J. Electr. Power 2024, 155, 109619. [Google Scholar] [CrossRef]
  25. Zheng, W.; Zhu, J.; Luo, Q. Distributed Dispatch of Integrated Electricity-Heat Systems with Variable Mass Flow. IEEE Trans. Smart Grid 2023, 14, 1907–1919. [Google Scholar] [CrossRef]
  26. Xue, Y.; Shahidehpour, M.; Pan, Z.; Wang, B.; Zhou, Q.; Guo, Q.; Sun, H. Reconfiguration of District Heating Network for Operational Flexibility Enhancement in Power System Unit Commitment. IEEE Trans. Sustain. Energy 2021, 12, 1161–1173. [Google Scholar] [CrossRef]
  27. Saravi, V.S.; Kalantar, M.; Anvari-Moghaddam, A. Resilience-Constrained Expansion Planning of Integrated Power–Gas–Heat Distribution Networks. Appl. Energy 2022, 323, 119315. [Google Scholar] [CrossRef]
  28. Zhang, Y.; Li, J.; Ji, X.; Ye, P.; Yu, D.; Zhang, B. Optimal Dispatching of Electric-heat-hydrogen Integrated Energy System Based on Stackelberg Game. Energy Convers. Econ. 2023, 4, 267–275. [Google Scholar] [CrossRef]
  29. Hou, H.; Liu, P.; Xiao, Z.; Deng, X.; Huang, L.; Zhang, R.; Xie, C. Capacity Configuration Optimization of Standalone Multi-energy Hub Considering Electricity, Heat and Hydrogen Uncertainty. Energy Convers. Econ. 2021, 2, 122–132. [Google Scholar] [CrossRef]
  30. Jannati-Oskuee, M.-R.; Gholizadeh-Roshanagh, R.; Najafi-Ravadanegh, S. Robust Expansion Planning of Reliable and Resilient Smart Distribution Network Equipped by Storage Devices. Electr. Power Syst. Res. 2024, 229, 110081. [Google Scholar] [CrossRef]
  31. Pamshetti, V.B.; Singh, S.; Thakur, A.K.; Singh, S.P.; Babu, T.S.; Patnaik, N.; Krishna, G.H. Cooperative Operational Planning Model for Distributed Energy Resources with Soft Open Point in Active Distribution Network. IEEE Trans. Ind. Appl. 2023, 59, 2140–2151. [Google Scholar] [CrossRef]
  32. Zhang, L.; Yu, S.; Zhang, B.; Li, G.; Cai, Y.; Tang, W. Outage Management of Hybrid AC/DC Distribution Systems: Co-Optimize Service Restoration with Repair Crew and Mobile Energy Storage System Dispatch. Appl. Energy 2023, 335, 120422. [Google Scholar] [CrossRef]
  33. Chen, C.; Li, Y.; Qiu, W.; Liu, C.; Zhang, Q.; Li, Z.; Lin, Z.; Yang, L. Cooperative-Game-Based Day-Ahead Scheduling of Local Integrated Energy Systems With Shared Energy Storage. IEEE Trans. Sustain. Energy 2022, 13, 1994–2011. [Google Scholar] [CrossRef]
  34. Niu, C.; Chen, C.; Weng, G.; Ren, J.; Gao, Q.; Lin, Z. Optimal Scheduling Strategy for Integrated Electric-Thermal-Gas Energy Systems Considering Multiple Types of Energy Storage Resources. In Proceedings of the 2023 IEEE 7th Conference on Energy Internet and Energy System Integration (EI2), Hangzhou, China, 15 December 2023; pp. 5014–5020. [Google Scholar]
Figure 1. Structure of the DN-DHN-BESS system.
Figure 1. Structure of the DN-DHN-BESS system.
Applsci 14 07699 g001
Figure 2. Distributed cooperative dispatch framework for DN-DHN-BESS.
Figure 2. Distributed cooperative dispatch framework for DN-DHN-BESS.
Applsci 14 07699 g002
Figure 3. Test topology containing a 33-node DN and a 44-node DHN.
Figure 3. Test topology containing a 33-node DN and a 44-node DHN.
Applsci 14 07699 g003
Figure 4. Comparison of curtailed power for PV and WT of the proposed method and independent dispatch method.
Figure 4. Comparison of curtailed power for PV and WT of the proposed method and independent dispatch method.
Applsci 14 07699 g004
Figure 5. Comparison of the power purchased from the superior grid by the DN of the proposed method and independent dispatch method.
Figure 5. Comparison of the power purchased from the superior grid by the DN of the proposed method and independent dispatch method.
Applsci 14 07699 g005
Figure 6. Comparison of the water supply temperature of DHN of the proposed method and Case 1-1.
Figure 6. Comparison of the water supply temperature of DHN of the proposed method and Case 1-1.
Applsci 14 07699 g006
Figure 7. Comparison of the heat generation power of DHN of the proposed method and Case 1-1.
Figure 7. Comparison of the heat generation power of DHN of the proposed method and Case 1-1.
Applsci 14 07699 g007
Figure 8. Comparison of curtailed power for PV and WT of the proposed method and Case 1-1.
Figure 8. Comparison of curtailed power for PV and WT of the proposed method and Case 1-1.
Applsci 14 07699 g008
Figure 9. Trend in total system cost and curtailed energy for PV and WT of the four cases.
Figure 9. Trend in total system cost and curtailed energy for PV and WT of the four cases.
Applsci 14 07699 g009
Table 1. Comparison of system cost of the proposed method and centralized dispatch method.
Table 1. Comparison of system cost of the proposed method and centralized dispatch method.
Cost (103 CNY)Proposed MethodCentralized Dispatch Method
CDN9.459.38
CDHN13.7613.72
CBESS−0.105−0.117
Total23.1022.98
Table 2. Comparison of computation time of the proposed method and centralized dispatch method.
Table 2. Comparison of computation time of the proposed method and centralized dispatch method.
Proposed MethodCentralized Dispatch Method
Computation time75.46 s2.98 s
Table 3. Comparison of system cost of the proposed method and independent dispatch method.
Table 3. Comparison of system cost of the proposed method and independent dispatch method.
Cost (103 CNY)Proposed MethodIndependent Dispatch Method
CDN9.4524.33
CDHN13.7611.23
CBESS−0.105−0.307
Total23.1035.26
Table 4. Comparison of curtailed energy for PV and WT of the proposed method and independent dispatch method.
Table 4. Comparison of curtailed energy for PV and WT of the proposed method and independent dispatch method.
Curtailed Energy for PV (kWh)Curtailed Energy for WT (kWh)
Proposed method1185.30
Independent dispatch method8578.50
Table 5. Comparison of system cost of the proposed method and Case 1-1.
Table 5. Comparison of system cost of the proposed method and Case 1-1.
Cost (103 CNY)Proposed MethodCase 1-1
CDN9.4511.54
CHN13.7614.41
CES−0.105−0.221
Total23.1025.73
Table 6. Comparison of curtailed energy for PV and WT of the proposed method and Case 1-1.
Table 6. Comparison of curtailed energy for PV and WT of the proposed method and Case 1-1.
Curtailed Energy for PV (kWh)Curtailed Energy for WT (kWh)
Proposed method1185.30
Case 1-12349.60
Table 7. Comparison of system costs of the four cases.
Table 7. Comparison of system costs of the four cases.
Cost (103 CNY)Case 2-1Case 2-2Case 2-2Case 2-3
CDN9.829.459.429.41
CHN14.2713.7613.6613.65
CES−0.087−0.105−0.205−0.217
Total24.0023.1022.8822.84
Table 8. Comparison of curtailed energy for PV and WT of the four cases.
Table 8. Comparison of curtailed energy for PV and WT of the four cases.
Curtailed Energy for PV (kWh)Curtailed Energy for WT (kWh)
Case 2-11222.90
Case 2-21185.30
Case 2-31146.30
Case 2-41135.10
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Fu, X.; Yu, S.; He, Q.; Wang, L.; Chen, C.; Niu, C.; Lin, Z. Distributed Cooperative Dispatch Method of Distribution Network with District Heat Network and Battery Energy Storage System Considering Flexible Regulation Capability. Appl. Sci. 2024, 14, 7699. https://doi.org/10.3390/app14177699

AMA Style

Fu X, Yu S, He Q, Wang L, Chen C, Niu C, Lin Z. Distributed Cooperative Dispatch Method of Distribution Network with District Heat Network and Battery Energy Storage System Considering Flexible Regulation Capability. Applied Sciences. 2024; 14(17):7699. https://doi.org/10.3390/app14177699

Chicago/Turabian Style

Fu, Xin, Shunjiang Yu, Qibo He, Long Wang, Changming Chen, Chengeng Niu, and Zhenzhi Lin. 2024. "Distributed Cooperative Dispatch Method of Distribution Network with District Heat Network and Battery Energy Storage System Considering Flexible Regulation Capability" Applied Sciences 14, no. 17: 7699. https://doi.org/10.3390/app14177699

APA Style

Fu, X., Yu, S., He, Q., Wang, L., Chen, C., Niu, C., & Lin, Z. (2024). Distributed Cooperative Dispatch Method of Distribution Network with District Heat Network and Battery Energy Storage System Considering Flexible Regulation Capability. Applied Sciences, 14(17), 7699. https://doi.org/10.3390/app14177699

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