Next Article in Journal
Analysis of Natural Vibration Characteristics of Metal Roof Panel of Large-Span Standing-Seam Metal Roof System
Next Article in Special Issue
Peer-to-Peer Transactive Computation–Electricity Trading for Interconnected Virtual Power Plant Buildings
Previous Article in Journal
Research of Heat Tolerance and Moisture Conditions of New Worked-Out Face Structures with Complete Gap Spacings
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Cluster Scheduling of Active–Reactive Power for Distribution Network Considering Aggregated Flexibility of Heterogeneous Building-Integrated DERs

1
Guizhou Power Grid Co., Ltd., Guiyang 550002, China
2
State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources, North China Electric Power University, Beijing 102206, China
*
Author to whom correspondence should be addressed.
Buildings 2023, 13(11), 2854; https://doi.org/10.3390/buildings13112854
Submission received: 12 September 2023 / Revised: 3 November 2023 / Accepted: 8 November 2023 / Published: 14 November 2023
(This article belongs to the Special Issue Strategies for Building Energy Efficiency)

Abstract

:
This paper proposes an active–reactive power collaborative scheduling model with cluster division for the flexible distributed energy resources (DERs) of smart-building systems to resolve the high complexity of the centralized optimal scheduling of massive dispersed DERs in the distribution network. Specifically, the optimization objective of each cluster is to minimize the operational cost, the power-loss cost, and the penalty cost for flexibility deficiency, and the second-order cone-based branch flow method is utilized to convert the power-flow equations into linearized cone constraints, reducing the nonlinearity and heavy computation burden of the scheduling model. Customized virtual battery models for building-integrated flexible DERs are developed to aggregate the power characteristics of flexible resources while quantifying their regulation capacities with time-shifting power and energy boundaries. Moreover, a cluster division algorithm considering the module degree index based on the electrical distance and the flexible balance contribution index is formulated for cluster division to achieve information exchange and energy interaction in the distribution network with a high proportion of building-integrated flexible DERs. Comparative studies have demonstrated the superior performance of the proposed methodology in economic merits and voltage regulation.

1. Introduction

1.1. Relevant Background

The construction sector is the third-largest energy consumer, accounting for 33% of global energy demand and 16% of global CO2 emissions in 2022 [1]. Considering the growing tendency of urbanization progress and the continuous development of the power industry, the proportion of construction energy consumption is expected to increase by 40% in 2030 and even 50% in 2050 [2]. The aggregation of flexible DERs in buildings is recognized as a promising alternative method to reduce CO2 emissions and energy consumption, promoting the profound integration of distributed energy resources (DERs) and building energy supply [3]. Smart-building systems equipped with power-access points could achieve flexible resource aggregation and exploit the regulation capacity of these resources sufficiently, intensifying the interdependency between buildings and electricity networks [3,4]. Nevertheless, the location of building-integrated photovoltaics (BIPVs), building-integrated electric vehicles (BIEVs), building-integrated energy storage (BIES), building-integrated temperature-controlled load (BITC), and other types of comprehensive energy production and storage equipment in smart-building systems are generally scattered, resulting in the high complexity of centralized optimization scheduling due to the lack of information exchange and energy interaction [4,5,6]. Consequently, the optimization scheduling of a distribution network with cluster division for flexible resources integrated into smart-building systems would be ideal for the high-penetrated deployment of DERs to meet carbon-neutral targets.
One typical method studied in [7] for cluster optimization scheduling is to use the second-order neighborhood similarity matrix and modular matrix to extract relatively more comprehensive complex network feature information for network division. Another example is a chromosome-encoding method considering the modular index based on the electrical design to satisfy the demand of distribution network programmers [8]. In further research, a K-means clustering algorithm based on the variation sensitivity of voltage to power between nodes to solve the voltage overlimit and voltage fluctuations is used [9]. Moreover, to tackle difficulties in uncertainty and centralized control, a multi-time rolling dynamic cluster division method with the equilibrium degree of reactive power for active distribution networks was developed [10]. However, the cluster division in these studies only involves a single index without considering comprehensive indexes, including flexibility regulation capacity and electrical connection. In addition, the active–reactive power collaborative optimization scheduling of distribution networks is acknowledged as an effective strategy for handling the uncertainties from dispersive DERs. A method for the coordinated optimal operation scheduling of active distribution networks was developed to meet the electricity demand of the building itself with interconnected electric vehicles and integrated RESs [11]. To implement optimal power dispatch and significantly reduce the overall operation cost of the microgrid, a method for optimal energy and power management of microgrids consisting of mega buildings, plug-in electric vehicles (PEVs), and renewable energy sources (RESs) was designed [12]. Furthermore, an energy management system (EMS) for microgrids of building prosumers based on a hierarchical multi-agent system (MAS) was studied in [13] to minimize the operation cost of a microgrid, considering many operation and technical constraints.
Integration and quantification for the regulation capacity of flexible resources are regarded as effective means to take advantage of the flexibility of smart-building systems [4,14,15,16,17,18]. Several important investigations have been reported on the virtual battery model [14], virtual synchronous machine model [15], flexible supply-and-demand balance model [16], and node power model [17] for the aggregation of flexible resources. A generic and scalable approach considering zonotopic sets for flexible energy systems was studied in [18] to describe their flexibility quantitatively. The authors use a constraint space-superposition method based on Minkowski Sum to seek the maximal inner cube constraints for the virtual synchronous machine of a heat pump cluster [19]. Moreover, the authors in [20] developed a day-ahead optimal scheduling model based on the generalized virtual battery model to determine the energy schedule and reserve capacity of electric vehicles and air conditioning. Additionally, a flexibility assessment can represent the abundance of flexible resources in a distribution network and has been comprehensively investigated with the insufficient ramping resource expectation index [21], the network flexibility index [22], the technical uncertain scenario flexibility index [23], and the technical–economic uncertainty scenario flexibility index [24]. A quantitative assessment method for power system flexibility based on probabilistic optimal power flow was developed to improve system flexibility and operating economy [21,22]. The authors in [23] studied a new distribution system flexibility evaluation method based on cloud models to address the difficulty in fully characterizing the flexibility status of a system. Furthermore, to quantify the transmission capacity of power grids and effectively improve the flexibility of power systems, an evaluation method of supply-and-demand balance considering flexible carrying capacity was studied in [24]. However, the comprehensive effect of the indexes on the flexibility state of the system is usually neglected in these studies.
Aggregation techniques, such as the k-order approximate model [25], the multi-time-scale approximate model [26], the real-time aggregation flexibility feedback method [27], and the Minkowski Sum [19] have been exploited to simplify the constraints of the aggregation feasible domains and realize the flexibility of the system. The authors studied a virtual power plant aggregation and operation mechanism to effectively guarantee the stability of the aggregation results [14]. An adaptive and robust optimization algorithm was adopted to accurately quantify the maximum power support range of the heat grid to the power grid [17]. After the aggregation of the flexible resources, the coordination and optimization of the distribution network with cluster division based on the flexible resources in the smart-building system can be achieved. Nevertheless, the potential active and reactive power support capacitates from the clusters of building-integrated flexible DERs for distribution networks are not involved.

1.2. Contribution

In this paper, an active–reactive power collaborative scheduling model with cluster division is proposed for the flexible DERs of smart-building systems to resolve the high complexity of centralized optimal scheduling for massively dispersed DERs in the distribution network. The detailed contributions of this paper are summarized as follows:
  • The customized virtual battery models are developed for building-integrated flexible DERs to aggregate the power characteristics while quantifying the operation and regulation capacity with the time-shifting energy state, the energy boundary, and the power boundary. A constrained space-superposition method with Minkowski Sum is exploited to derive the summation of the regulation capacity of building-integrated flexible DERs.
  • A cluster division algorithm considering the structural and functional cluster division indexes for the building-integrated flexible DERs is proposed to automatically divide the distribution network with a high proportion of renewable energy and flexible DERs into the optimal clusters. The structural cluster division index refers to the module degree index based on the electrical distance, which can ensure the close electrical connection of the nodes within the clusters, and the functional cluster division index represents the flexible balance contribution index, measuring the matching degree between the distributed power consumption demands and the supply of building-integrated flexible DERs.
  • An optimal active–reactive power collaborative scheduling model for the clusters of building-integrated flexible DERs is formulated to minimize the operational cost, the power-loss cost, and the penalty cost for flexibility deficiency, in which the second-order cone-based branch flow method is exploited to decompose the nonlinearity and nonconvexity of the power-flow equation and the nodal voltage magnitude constraint into the linearized cone model by introducing the intermediate variable.

2. Flexible Regulation Capacity of Flexible DERs in a Smart-Building System

2.1. Smart-Building System

Figure 1 depicts the typical modality of a smart-building system among a highly renewable-penetrated distribution network. The smart-building system is generally equipped with cogeneration systems, distributed photovoltaics, electric vehicles, heat pumps, and energy storage to satisfy the demand for electricity and the cold and heat load in the building [5]. To fully utilize the space and time flexibility of building-integrated flexible DERs, the establishment of a virtual battery model for flexible DERs is necessary, which can quantify and integrate the flexibility and regulation capacity of these resources.

2.2. Battery Model of ES in a Smart-Building System

BIES can transfer the electricity generated by BIPVs to balance the energy of the power system, improving the reliability of system operation and power supply quality. A battery model of ES in a smart-building system (1), consisting of charge power constraints, discharge power constraints, and energy state constraints, is developed to describe the operating characteristic mathematically:
P ch , i , t , min ES P ch , i , t ES P ch , i , t , max ES
P dis , i , t , min ES P dis , i , t ES P dis , i , t , max ES
E i , t , min ES E i , t ES E i , t , max ES
E i , t + 1 ES = E i , t ES + P i , t ES Δ t = E i , t ES + ( η ch ES P ch , i , t ES P dis , i , t ES η dis ES ) Δ t
where P c h , i , t E S   and P d i s , i , t E S are the charging and discharging power of BIES connected at node i and time t ; P c h , i , t , m i n E S , P c h , i , t , m a x E S , P d i s , i , t , m i n E S , and P d i s , i , t , m a x E S denote the lower and upper bounds of the charging and discharging power of BIES; E i , t E S   and P i , t E S are the operation capacity and the actual power of BIES; E i , t , m a x E S and E i , t , m i m E S are the threshold of the operation capacity of BIES; η c h E S and η d i s E S are the charge and discharge efficiency coefficients of BIES, and t is the time interval.
The flexible regulation capacity of ES in smart-building systems with energy time-shift capacity and power regulation capacity can be expressed as
F i , t ES , sup , up = P i , t , max ES P i , t ES
F i , t ES , sup , dn = P i , t ES P i , t , min ES
where F i , t E S , s u p , u p and F i , t E S , s u p , d n denote the upgraded flexible energy supply and the downgraded flexible energy supply of BIES connected at node i and time t , deriving from the lower and upper bounds of the actual power P i , t , m i n E S and P i , t , m a x E S of BIES.

2.3. Virtual Battery (VB) Model in a Smart-Building System

The VB model can accurately describe the operation and regulation characteristics of flexible DERs by constraining energy state, energy boundary, and power boundary, which can be expressed as
P i , t , min G P i , t G P i , t , max G
E i , t , min G E i , t G E i , t , max G
E i , t + 1 G = E i , t G + P i , t G Δ t + Δ E i , t G
where P i , t G and E i , t G are the regulation power and the flexibility reserve energy of flexible DERs in a smart-building system connected at node i and time t ; P i , t , m i n G , P i , t , m a x G , E i , t , m i n G , and E i , t , m a x G denote the lower and upper bounds of the regulation power and the flexibility reserve energy; (3c) represents the state of the flexibility reserve energy for the flexible DERs, and E i , t G is the effect of other factors on the electric energy of the VB model.

2.4. VB Model of EV in a Smart-Building System

In general, the power battery of BIEV can achieve power interaction with the grid through flexibly switching between charging, placing, and discharging states. The grid-connected operation status of BIEV is illustrated in Figure 2, where BIEV with an initial operating capacity of E i , t 0 E V are connected to the grid with the maximum charge and discharge power of P c h , i , t , m a x E V   and   P d i s , i , t , m a x E V at time t 0 , reaching the expected operating capacity of E i , t E V when off the grid at time t 1 .
Thus, the virtual battery model of BIEV can be formulated as follows:
P ch , i , t , min EV P ch , i , t EV P ch , i , t , max EV
P di s , i , t , min EV P di s , i , t EV P di s , i , t , max EV
E i , t 0 EV P dis , i , t , max EV ( t t 0 ) E i , t EV E i , t 0 EV + P ch , i , t , max EV ( t t 0 )
E i , t + 1 EV = E i , t EV + P i , t EV Δ t = E i , t EV + ( η ch EV P ch , i , t EV P dis , i , t EV η dis EV ) Δ t
where P c h , i , t , m i n E V and P d i s , i , t , m i n E V are the minimum thresholds of the charging and discharging power of BIEV connected at node i and time t ; P c h , i , t E V , P d i s , i , t E V P i , t E V , and E i , t E V denote the charging power, discharging power, actual power, and operation capacity; η c h E V and η d i s E V are the charge and discharge efficiency coefficients of BIEV. (4c) and (4d) represent the boundary of operation capacity and the operation state.
Then, considering that the power boundary P i , t , m a x E V and P i , t , m i n E V   of BIEV are affected by the operation capacity, the upgraded flexible energy supply F i , t E V , s u p , u p and the downgraded flexible energy supply F i , t E V , s u p , d n of EV in a smart-building system can be expressed as
F i , t EV , sup , up = min { P i , t , max EV P i , t EV , E i , t , max EV E i , t EV Δ t }
F i , t EV , sup , dn = min { P i , t EV P i , t , min EV P , E i , t EV E i , t , min EV Δ t }

2.5. VB Model of TCL in a Smart-Building System

The equivalent thermal parameter model is adopted to describe the thermal dynamics characteristics of TCL, which equates the internal environment, external environment, and heating (cooling) capacity to circuit devices to analyze the dynamic change relationship of the indoor temperature and power [28,29]. As illustrated in Figure 3, the heating capacity Q i , t 1 of the temperature control load can be defined as the sum of the heat absorbed by the room Q i , t 2   and the air convection heat Q i , t 3 .
(1d) shows that the energy of the BIES is equal to the superposition of the energy in the past state, the storage energy, and the release energy in the current state. Similarly, the indoor temperature T i n , i , t T C L of the TCL is equal to the summation of the indoor temperature in the past state, warming, and cooling in the current state. When ambient temperature T o u t , i , t T C L is higher than indoor temperature, the temperature of the heating equipment T h e a t , i , t T C L and the ambient temperature together transfers the heat quantity from high temperature to low temperature, increasing indoor temperature and the energy of the VB model. On the contrary, when the indoor temperature is higher in the summer, the cooling capacity of TCL works to reduce the indoor temperature, which is equivalent to the decrease in the energy of the VB model, as follows:
{ Q i , t 1 = Q i , t 2 + Q i , t 3 Q i , t 1 = A 1 η TCL ( T heat , i , t TCL T in , i , t TCL ) Q i , t 2 = ρ c V ( T in , i , t + 1 TCL T in , i , t TCL ) Q i , t 3 = A 2 η Wall ( T in , i , t TCL T out , i , t TCL ) η COP = P i , t TCL Δ t Q i , t 1
P i , t , min TCL P i , t TCL P i , t , max TCL
T base σ T in , i , t TCL T base + σ
T in , i , t + 1 TCL = ( 1 A 2 η Wall ρ c V ) T in , i , t TCL + A 2 η Wall ρ c V T out , i , t TCL + P i , t TCL ρ c V η COP Δ t
where A 1 and A 2 are the heat-dissipation area of BITC connected at node i and time t and the wall; η T C L , η W a l l , and η C O P denote the heating efficiency, wall surface thermal radiation rate, and energy conversion efficiency; ρ , c , and V are the indoor air density, air heat capacity, and room volume. (6b) and (6c) represent the boundary of actual power and indoor temperature, where T b a s e and σ are the indoor temperature baseline and the dead zone value, respectively.
The upgraded flexible energy supply F i , t T C L , s u p , u p and the downgraded flexible energy supply F i , t T C L , s u p , d n of TCL in smart-building systems can be expressed as
F i , t TCL , sup , up = P i , t , max TCL P i , t TCL
F i , t TCL , sup , dn = P i , t TCL P i , t , min TCL

3. Active–Reactive Power Collaborative Optimization Scheduling of a Distribution Network Cluster with the Aggregation of Flexible DERs in a Smart-Building System

3.1. Cluster Division Index of Building-Integrated Flexible DERs

Structural and functional cluster division indexes of building-integrated flexible DERs are designed for a distribution network with a high penetration of renewable energy to achieve independent autonomy within clusters and coordinated interaction among clusters. To ensure the close electrical connection of nodes within the clusters, the module degree index based on electrical distance is adopted to describe the connection strength of nodes in the structure [8]. The functional index refers to the flexible balance contribution index, which can evaluate the regulation capacity of building-integrated flexible DERs to ensure the power balance under a high penetration rate of distributed energy access.

3.1.1. Module Degree Index

Generally, electrical distance L i j is exploited to measure the tightness of electrical coupling between nodes in the distribution network based on the sensitivity relationship of node power variation and node voltage [8], and thus the module degree index ξ for the cluster division structure strength of complex networks is formulated as follows:
ξ = 1 2 m i j ( μ i j r i r j 2 m ) δ ( i , j )
{ μ i j = 1 L i j L i j , max m = i j μ i j r i = i μ i j
L i j = ( d i 1 d j 1 ) 2 + ( d i 2 d j 2 ) 2 + ... + ( d i n d j n ) 2
d i j = ( S i i P S i j P ) + ( S i i Q S i j Q )
Δ U = S P Δ P + S Q Δ Q
where μ i j , m , and r i denote the edge weight connecting node i and j , the sum of the edge weight in the network, and the edge weight connecting node i ; δ i , j , a 0–1 variable, equals 1 if Node 1 and Node 2 are in the same cluster. Since there is a coupling relationship between the two nodes and the surrounding nodes, L i j can be formulated with the combined effect d i j of the power change in node i and j as shown in (8c) and (8d); L i j , m a x is the maximum threshold of electrical distance; (8e) indicates the sensitivity relationship between node voltage and node-injected power; U , P , and Q are the variation of voltage amplitude, active power, and reactive power; S P and S Q are the sensitivity matrix of voltage-active power and voltage-reactive power.

3.1.2. Flexible Balance Contribution Index

The regulation capacity of building-integrated flexible DERs refers to the flexible energy supply capacity required to satisfy the flexible demands of net load in smart-building systems, which can be divided into upgraded flexible energy supply capacity and downgraded flexible energy supply capacity. The flexible balance contribution index φ can measure the matching degree between distributed power consumption demands and the supply of the building-integrated flexible DERs by calculating the equilibrium proportion of flexible DERs bearing node flexibility demands at a certain moment as follows:
φ = 1 n ( i = 1 n α i , t G , up + i = 1 n α i , t G , dn )
{ α i , t G , up = F i , t G , sup , up F i , t dem α i , t G , dn = F i , t G , sup , dn F i , t dem       G { ES , EV , TCL }
F i , t dem = P i , t Load P i , t PV
where α i , t G , u p and α i , t G , d n denote the upgraded and downgraded flexible balance contribution degree. (9c) represents the flexible demands of net load in smart-building systems, deriving from the difference between load demands P i , t L o a d   and BIPV generation.
Then, the cluster division index γ of building-integrated flexible DERs considering the module degree index and the flexible balance contribution index with different weights a 1 and a 2 can be calculated from Formulas (8a) and (9a),
γ = a 1 ξ + a 2 φ
The critical procedures of the cluster division process of building-integrated flexible DERs are summarized as follows:
  • Each node in the distribution network is regarded as a separate cluster initially to calculate the cluster division index γ of building-integrated flexible DERs based on Formula (10);
  • Node j is randomly selected from the remaining nodes to merge with node i , forming a new cluster k ( i , j ) , and then the variation of the cluster division index γ = γ γ can be derived from the division index of cluster k ( i , j ) . When γ reaches the maximum positive value, the two nodes can be divided into the same cluster;
  • The new cluster k ( i , j ) will be regarded as a new node to repeat the second procedure. The division procedures will continue until the nodes in the network cannot merge and the cluster division index γ of building-integrated flexible DERs reaches the maximum.

3.2. Active–Reactive Power Collaborative Optimization Scheduling within Clusters

3.2.1. Optimization Objective

The active–reactive power collaborative optimization scheduling within clusters can be formulated as a multi-object optimization model based on the flexibility of building-integrated flexible DERs, while the uncertainties in PV outputs, ES performance, EV parking behaviors, TCL heating/cooling behaviors, and SVC regulation capacity are considered and represented by various stochastic scenarios based on historical data in [3,30,31]. Specifically, each cluster k aims to minimize operational cost f k 1 , power-loss cost f k 2 , and penalty cost for flexibility deficiency f k 3 , as shown in (11a), (12), and (13a). f k 1 is composed of the maintenance cost C k , i , t , s P V of BIPVs due to the renewable generation curtailment, the lifetime degradation cost C k , i , t , s E S of BIESs, the compensation cost C k , i , t , s C L of BIEVs and BITCs stemmed from power regulation and the operation cost C k , i , t , s S V C of static var compensator (SVC) under scenario s . f k 2 can be derived from the power-flow calculation during the scheduling circulation, where U k , i , t , s , G k , i , t , s , B k , i , t , s , and δ k , i j , s represent the nodal voltage magnitudes, the conductance, the susceptance, and the phase angle of line i j at time t within cluster k under scenario s ; λ k , s L o s s is the compensation price of the network loss. The penalty cost f k 3 for the flexibility deficiency is introduced when the flexibility supply cannot satisfy the load demands within the clusters, including upgraded flexibility deficiency F k , t , s d e f i , u p and downgraded flexibility deficiency F k , t , s d e f i , d n under scenario s , as follows:
f k 1 = t Γ s = 1 N s π s ( i Ω PV N k C k , i , t , s PV + i Ω ES N k C k , i , t , s ES + i Ω CL N k C k , i , t , s CL + i Ω SVC N k C k , i , t , s SVC ) i , k , s
C k , i , t , s PV = λ k , s PV P k , i , t , s PV = λ k , s PV ( P k , i , t , s PV 0 Δ P k , i , t , s PV )
C k , i , t , s ES = λ k , s ES | P k , i , t , s ES |
C k , i , t , s CL = λ k , s EV | P k , i , t , s EV | + λ k , s TCL | P k , i , t , s TCL |
C k , i , t , s SVC = λ k , s SVC | Q k , i , t , s SVC |
f k 2 = t Γ s = 1 N s π s ( λ k , s Loss ( i , j ) N k G k , i j , s ( ( U k , i , t , s ) 2 + ( U k , j , t , s ) 2 2 U k , i , t , s U k , j , t , s cos δ k , i j , s ) ) i , k , s
f k 3 = t Γ s = 1 N s π s ( λ k , s defi ( F k , t , s defi , up + F k , t , s defi , dn ) ) i , k , s
F k , t , s defi , up = i N k Δ F k , i , t , s up = { i N k | F k , i , t , s G , sup , up F k , i , t , s dem | Δ F k , i , t , s up 0 0 Δ F k , i , t , s up 0
F k , t , s defi , dn = i N k Δ F k , i , t , s dn = { i N k | F k , i , t , s G , sup , dn F k , i , t , s dem | Δ F k , i , t , s dn 0 0 Δ F k , i , t , s dn 0
where Ω P V , Ω E S , Ω C L , and Ω S V G are the node sets of BIPV, BIES, BIEV, BITC, and SVC; N k and K indicate the k -th distributed resource cluster and the total number of divided clusters; N s is the number of scenarios for stochastic optimization; π s is the probability of scenario s , and the sum of probabilities for all scenarios is equal to 1. Due to the uncertainties in the balancing market prices, the unit operation prices λ k P V , λ k E S , λ k E V , λ k T C L , and λ k S V C of BIPV, BIES, BIEV, BITC, and SVC connected at node i and time t within cluster k are also uncertain and can be represented by stochastic scenarios with associated occurrence probabilities; P k , i , t , s P V , P k , i , t , s P V 0 , and Δ P k , i , t , s P V denote the actual generation, predicted generation, and generation curtailment of BIPV. Because of the frequent charging–discharging behaviors of building-integrated flexible DERs, the aging of these resources is inevitable where P k , i , t , s E S , P k , i , t , s E V , and P k , i , t , s T C L represent the actual power of BIES, BIEV, and BITC; Q k , i , t , s S V C is the reactive power output of SVC. (13b) and (13c) indicate the upgraded and downgraded flexibility deficiency deriving from the flexibility margin Δ F k , i , t , s u p and Δ F k , i , t , s d n of the negative value, and λ k , s d e f i is the penalty factor for flexibility deficiency.
Considering operational cost, power-loss cost, and penalty cost for flexibility deficiency comprehensively, the multi-objective function with different weights λ 1 , λ 2 , and λ 3 can be represented as
min f k = λ 1 f k 1 + λ 2 f k 2 + λ 3 f k 3

3.2.2. Operation Constraints

The Newton–Raphson method is adopted to describe power flows within the clusters, as shown in (15a). Constraint (15b) represents the active and reactive power balance at each bus, and constraint (15c) imposes the upper and lower limits on the nodal voltage magnitudes, U k , i , t , m a x and U k , i , t , m i n .
{ Δ P k , i , t , s U k , i , t , s ( i , j ) N k U k , j , t , s ( G k , i j , s cos δ k , i j , s + B k , i j , s sin δ k , i j , s ) = 0 Δ Q k , i , t , s U k , i , t , s ( i , j ) N k U k , j , t , s ( G k , i j , s sin δ k , i j , s B k , i j , s cos δ k , i j , s ) = 0
{ Δ P k , i , t , s = P k , , i , t , s PV + P k , i , t , s ES P k , i , t , s EV P k , i , t , s TCL P k , i , t , s Load Δ Q k , i , t , s = Q k , i , t , s PV + Q k , i , t , s ES + Q k , i , t , s SVC
U k , i , t , min U k , i , t , s U k , i , t , max
The power output curtailment of BIPV is supposed to be lower than the specified maximum thresholds of active and reactive power, Δ P k , i , t , m a x P V and Δ Q k , i , t , m a x P V . Moreover, the sufficient power factor c o s θ s under scenario s for the operation of a photovoltaic inverter is needed to limit the reactive power circulation, where the amount of reactive power injected or absorbed is constrained by the bounds in (16a) and (16b),
( P k , i , t , s PV ) 2 + ( Q k , i , t , s PV ) 2 ( S k , i , t , s PV ) 2
cos θ s ( P k , i , t , s PV ) 2 + ( Q k , i , t , s PV ) 2 P k , i , t , s PV
0 Δ P k , i , t , s PV Δ P k , i , t , max PV
0 Δ Q k , i , t , s PV Δ Q k , i , t , max PV
Also, the states of reactive power of SVC should be constrained within their lower and upper bounds Q k , i , t , m i n S V C and Q k , i , t , m a x S V C , as follows:
Q k , i , t , min SVC Q k , i , t , s SVC Q k , i , t , max SVC
In addition, the sum of the regulation capacity of building-integrated flexible DERs based on the constrained space-superposition method with Minkowski Sum can be expressed as Formula (18a), and the state of the total sum of the regulation capacity P k , t , s F l e x under scenario s should always be limited within the allowable lower and upper bounds P k , t , m i n F l e x and P k , t , m a x F l e x ,
P k , t , s Flex = P k , t , s ES P k , t , s EV P k , t , s TCL = i N k P k , i , t , s ES + i N k P k , i , t , s EV + i N k P k , i , t , s TCL
P k , t , min Flex P k , t , s Flex P k , t , max Flex
{ P k , t , max Flex = i N k P i , t , max ES + i N k P i , t , max EV + i N k P i , t , max TCL P k , t , min Flex = i N k P i , t , min ES + i N k P i , t , min EV + i N k P i , t , min TCL

3.2.3. Solution Methodology

Due to the nonlinearity and nonconvexity of the power-flow equations increasing the heavy computation burden of the optimization model, the second-order cone-based branch flow method is exploited to decompose this challenging problem into the linearized cone model by introducing intermediate variable X k , i , t , s , Y k , i , t , s , and Z k , i , t , s under scenario s and rotary cone constraint, as shown in (19).
X k , i , t , s = ( U k , i , t , s ) 2
Y k , i , t , s = U k , i , t , s U k , j , t , s cos δ k , i j , s
Z k , i , t , s = U k , i , t , s U k , j , t , s sin δ k , i j , s
X k , i , t , s X k , j , t , s = ( Y k , i , t , s ) 2 + ( Z k , i , t , s ) 2
Then, the power-loss cost function (12a), the power-flow Equation (15a), and the nodal voltage magnitudes constraint (15c) can be reformulated as follows:
f k 2 = t Γ s = 1 N s π s ( λ k , s Loss ( i , j ) N k G k , i j , s ( X k , i , t , s + X k , j , t , s 2 Y k , i , t , s ) ) i , k , s
{ Δ P k , i , t , s ( G k . i i , s X k , i , t , s + ( i , j ) N k j 1 n ( G k , i j , s Y k , i , t , s + B k , i j , s Z k , i , t , s ) ) = 0 Δ Q k , i , t , s ( B k , i i , s X k , i , t , s ( i , j ) N k j 1 n ( B k , i j , s Y k , i , t , s G k , i j , s Z k , i , t , s ) ) = 0
( U k , i , t , min ) 2 X k , i , t , s ( U k , i , t , max ) 2
When solving the mixed-integer second-order cone programming, a polyhedral approximate description of the second-order cone relaxation is exploited to describe formula (19d), and the original problem with a non-convex feasible domain is relaxed into a convex second-order cone-feasible domain using second-order cone relaxation, which can be solved quickly and effectively using a common solver, as follows:
X k , i , t , s X k , j , t , s ( Y k , i , t , s ) 2 + ( Z k , i , t , s ) 2
( Y k , i , t , s ) 2 + ( Z k , i , t , s ) 2 + ( X k , i , t , s X k , j , t , s 2 ) 2 X k , i , t , s + X k , j , t , s 4
( Y k , i , t , s ) 2 + ( Z k , i , t , s ) 2 ω k , i , t , s
( ω k , i , t , s ) 2 + ( X k , i , t , s X k , j , t , s 2 ) 2 X k , i , t , s + X k , j , t , s 4
As depicted in Figure 4, assuming that all feasible solutions of the non-convex feasible domain are on the second-order cone plane, the feasible domain is enlarged to convex cones after relaxation. Therefore, the solution within both the convex second-order cone and the original feasible domain is the optimal solution to the original problem.
In conclusion, the flowchart of the active–reactive power collaborative scheduling model with cluster division is shown in Figure 5.

4. Discussion

4.1. System Data

The proposed active–reactive power collaborative optimization scheduling model with cluster division is tested on a modified IEEE 33-bus distribution system, as depicted in Figure 6.
The installed capacities of BIPVs in nodes 1, 6, 13, 17, 18, 24, and 29 are set as 0.5 MW, 1.0 MW, 0.2 MW, 1.5 MW, 0.8 MW, 1.2 MW, and 0.4 MW to simulate the diverse PV endowments. The predicted energy outputs of BIPVs derived from the actual values in [6,30,31], as well as the estimated 24 h loads from [32,33], are utilized to generate possible scenarios. The random operation scenario is exploited to represent the uncertainty of PV output, ES performance, EV behavior, and TCL heating/cooling behaviors [34]. The initial random scenario tree has 5000 scenarios, and a scenario-reduction method is applied to decrease the number of scenarios to 50 to maintain the uncertainty approximation of the system while accelerating the calculation speed of the stochastic model. The prediction errors of BIPV generation and node load follow Gaussian distributions with zero means and a 10% standard deviation of the prediction values [31,32,35]. The detailed technical specifications of the building-integrated equipment and the distribution network are summarized in Table 1 [3,4,35].
The appropriate weights for the cluster division index of building-integrated DERs and the multi-object function are set as 0.41, 0.29, 0.63, 0.48, and 0.25 after many simulation tests. The optimization scheduling of the distribution network with cluster division for flexible DERs integrated into a smart-building system is implemented on an hourly basis over a 24 h horizon and solved by the CPLEX Optimizer on the MATLAB platform.

4.2. Simulation Results of the Random Scenarios

For the uncertainty problem in the active–reactive power collaborative optimization scheduling within clusters, multiple possible scenarios are derived from the basis values of forecasting results. Based on the sampling of forecasting errors on PV output, ES performance, EV behaviors, TCL heating/cooling behaviors, and electric load, a set of scenarios can be formed through Monte Carlo simulations. Taking the PV output at node 24 as an example, it can be found from Figure 7a that the initial random scenario is densely distributed around the predicted values. After the removal of similar scenarios by the scenario-reduction method, the scenarios retaining essential properties of the initial scenarios can be depicted in Figure 7b.

4.3. Comparative Results and Analysis

Three comparative schemes are performed for in-depth analysis of the effectiveness of the proposed methodology.
  • Scheme 1 performs the proposed optimization scheduling of the distribution network with cluster division for flexible DERs integrated into smart-building systems in Section 3.
  • Scheme 2 adopts the centralized optimization scheduling of the distribution network without considering the cluster division.
  • Scheme 3 is the initial distribution network before the optimization scheduling.

4.3.1. Cluster Division of Building-Integrated Flexible DERs

Considering the differences in node-voltage amplitude, power injection, and load demands at different times, the detailed data of each node at noon is exploited for cluster division. Figure 8 depicts the function curves for the cluster division index of building-integrated DERs corresponding to the number of different clusters. It can be found that the index function achieves the maximum value γ = 0.92 when the number of clusters is divided into 6. Therefore, the result of the optimal cluster division at K = 6 is shown in Figure 9.

4.3.2. Flexible Regulation Capacity Results of Building-Integrated Flexible DERs

Figure 10 illustrates the flexible regulation capacity of building-integrated flexible DERs based on the flexibility constraint space superimposed by Minkowski Sum. It can be observed from Figure 10 that the self-consistent rate is allowed to reflect the contribution of different flexibility resources to the flexibility balance of the power grid. For instance, the flexibility power region of building-integrated flexible DERs in Cluster 2 is larger than that of Cluster 1. In addition, the sum of the total amount of BIPV generation and the regulation capacity of building-integrated flexible DERs can satisfy the total load demands within Cluster 2, indicating that Cluster 2 achieves complete self-consistency of active power. On the contrary, the total load demands within Cluster 1 are greater than the sum of the total amount of BIPV generation and the regulation capacity of building-integrated flexible DERs at time t = 4 and t = 20 , thus giving rise to flexibility power deficiency.

4.3.3. Optimization Scheduling Results with Schemes 1–3

The power-loss and the node-voltage amplitude with Schemes 1–3 are illustrated in Figure 11, where the node-voltage at noon is selected for analysis since the generation of BIPVs reaches the peak at that time.
It can be found that the optimized system power loss with Scheme 1 is lower than that with Scheme 2 due to the comprehensive scheduling of the BIPV generation and the regulation capacity of building-integrated flexible DERs within clusters. In addition, the voltage-regulation effect of the centralized active power–reactive power cooperative optimization of the distribution network is unsatisfactory, maintaining the node voltage over the limit. However, since building-integrated flexible DERs are involved in the cluster power regulation, the node-voltage amplitude with Scheme 1 is within the specified range. Table 2 summarizes the statistical data of comparative economic performance results with Schemes 1 and 2, including the total cost, the operational cost, the power-loss cost, and the flexibility deficiency cost. Compared to Scheme 2, Scheme 1, considering the flexibility of building-integrated flexible DERs, can better estimate the flexibility regulation capacity, and thus, the operational cost and flexibility deficiency cost can be reduced by 11.69% and 9.21%, respectively. Consequently, Scheme 1 outperforms Schemes 2 and 3 on the economic merits and voltage regulation.

5. Conclusions

In this paper, an optimal active–reactive power collaborative scheduling model for building flexible DER clusters is proposed to resolve the uncertainties from the renewable energy in the distribution network. The key findings of this study are as follows: (1) The developed cluster division algorithm considering structural and functional indexes can achieve independent autonomy within clusters and coordinated interaction among clusters, contributing to decreasing the whole network loss; (2) Since the customized virtual battery models integrate and quantify the regulation capacity of flexible resources, the flexibility deficiency cost can be reduced by 9.21%; (3) The developed optimal active–reactive power collaborative scheduling model can achieve a better performance in economic merits and voltage regulation than the centralized optimization scheduling of a distribution network without considering the cluster division, whose operational cost can be decreased by 11.69%. In this article, the building-integrated flexible DERs, including BIES, BIEV, and BITC, were acknowledged as the resources on the demand side. Further research will focus on the collaborative flexibility regulation capacity of building-integrated flexible DERs on the supply side and demand side.

Author Contributions

Conceptualization, Y.F. and K.Z.; Methodology, S.H.; Software, J.Z.; Validation, L.Y., Y.F. and K.Z.; Formal analysis, S.H.; Investigation, Y.L.; Resources, J.Z.; Data curation, L.Y.; Writing—original draft preparation, Y.L.; Writing—review and editing, Y.L. and K.Z.; Visualization, Y.L.; Supervision, K.Z.; Project administration, Y.F.; Funding acquisition, Y.F., S.H., J.Z. and L.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was jointly supported by the Science and Technology Project of Guizhou Power Grid under Grant GZKJXM20222149 and National Key R&D Programs of China under Grant 2022YFB2403400.

Data Availability Statement

The data presented in this study are available in article.

Conflicts of Interest

Authors Yu Fu, Shuqing Hao, Junhao Zhang and Liwen Yu were employed by Guizhou Power Grid 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.

Nomenclature

Sets, Indices, and Function P k , t , m i n F l e x ,   P k , t , m a x F l e x Allowable lower and upper bounds of the regulation capacity
GSet of flexible resources Q k , i , t , m i n S V C ,   Q k , i , t , m a x S V C Lower and upper bounds of the reactive power of SVC
i , j Index of nodes Δ P k , i , t , m a x P V ,   Δ Q k , i , t , m a x P V Specified maximum thresholds of the active and reactive power connected at node i and time t within cluster k
k Index (set) of clustersVariables
t Index (set) of time slots P c h , i , t E S , P d i s , i , t E S Charging and discharging power of BIES connected at node i and time t
s Index (set) of scenarios
Ω P V ,   Ω E S ,   Ω C L ,   Ω S V C The node sets of BIPV, BIES, BIEV, BITC, and SVC
f k 1 Operational cost of the kth cluster E i , t E S Operation capacity of BIES connected at node i and timee t
f k 2 Power-loss cost of the kth cluster E i , t + 1 E S Operation state of BIES connected at node i and time t
f k 3 Penalty cost for flexibility deficiency of the kth cluster
Q i , t 1 Heating capacity of the temperature control load P i , t E S Actual power of BIES
Q i , t 2 Heat absorbed by the room P i , t G Regulation power of the flexible resources connected at node i and time t
Q i , t 3 Air convection heat E i , t G Flexibility reserve energy of the flexible resources connected at node i and time t
ξ The module degree index E i , t + 1 G State of the flexibility reserve energy of the flexible resources connected at node i and time t
φ Flexible balance contribution index
γ Building-integrated flexible resource cluster division index E i , t G Effect of other factors on the electric energy of the VB model
F i , t G , s u p , u p ,   F i , t G , s u p , d n Upgraded flexible energy supply and downgraded flexible energy supply of the flexible resources connected at node i and time t
Parameters F i , t E S , s u p , u p ,   F i , t E S , s u p , d n Upgraded flexible energy supply and downgraded flexible energy supply of BIES connected at node i and time t
t 0 BIEV grid connection time P c h , i , t E V ,   P d i s , i , t E V Charging and discharging power of BIEV connected at node i and time t
t 1 Off-grid time of BIEV P i , t E V Actual power of BIEV connected at node i and time t
Δ t Interval length of time slot E i , t E V Operation capacity of BIEV connected at node i and time t
η c h E V ,   η d i s E V Charge and discharge efficiency of BIEV E i , t + 1 E V Operation state of BIEV connected at node i and time t
A 1 ,   A 2 Heat-dissipation area of BITC and the wall F i , t E V , s u p , u p ,   F i , t E V , s u p , d n Upgraded flexible energy supply and downgraded flexible energy supply of BIEV connected at node i and time t
η T C L ,   η W a l l η C O P Heating efficiency, wall surface thermal radiation rate, and energy conversion efficiency T i n , i , t T C L Indoor temperature at time t
ρ Indoor air density T i n , i , t + 1 T C L Indoor temperature at time   t + 1
c Air heat capacity
V Room volume P i , t T C l Actual power of BITC connected at node i and time t
T b a s e Temperature baseline U The variation of the voltage amplitude
λ k , s L o s s Compensation price of network loss of the kth cluster
λ k , s d e f i Penalty factor for the flexibility deficiency of the kth cluster
λ k , s P V ,   λ k , s E S ,   λ k , s E V ,   λ k , s T C L ,   λ k , s S V C Unit operation prices of BIPV, BIES, BIEV, BITC, and SVC of the kth cluster P The variation of active power
P c h , i , t , m i n E S , P c h , i , t , m a x E S Lower and upper bounds of the charging power of BIES Q The variation of reactive power
P d i s , i , t , m i n E S , P d i s , i , t , m a x E S Lower and upper bounds of the discharging power of BIES d i j Combined effect of the power change at node j on node i
E i , t , m a x E S ,   E i , t , m i m E S Threshold of the operation capacity of BIES L i j Electrical distance between node i and node j
P i , t , m i n E S ,   P i , t , m a x E S Lower and upper bounds of the actual power of BIES μ i j The edge weight connecting node i and j
P i , t , m i n G ,   P i , t , m a x G Lower and upper bounds of the regulation power r i Edge weight connecting node i
E i , t , m i n G ,   E i , t , m a x G Lower and upper bounds of flexibility reserve energy α i , t G , u p ,   α i , t G , d n Upgraded and downgraded flexible balance contribution degree connected at node i and time t
P c h , i , t , m i n E V ,   P c h , i , t , m a x E V Minimum and maximum charging power of BIEV F i , t d e m Flexible demands of net load in the smart-building system connected at node i and time t
P i , t L o a d Load demands connected at node i and time t
P d i s , i , t , m i n E V ,   P d i s , i , t , m a x E V Minimum and maximum discharging power of BIEV P i , t P V BIPV generation connected at node i and time t
E i , t 0 E V Initial operating capacity of BIEV at time  t 0 C k , i , t , s P V Maintenance cost of BIPVs
P i , t , m i n E V ,   P i , t , m a x E V Lower and upper bounds of the actual power of BIEV
T o u t , i , t T C L Ambient temperature C k , i , t , s E S Maintenance cost of BIESs
T h e a t , i , t T C L Temperature of the heating equipment
P i , t , m i n T C L ,   P i , t , m a x T C L Lower and upper bounds of the actual power of BITC
σ Dead zone value C k , i , t , s C L Compensation cost of BIEVs and BITCs
S P ,   S Q The sensitivity matrix of voltage-active power and voltage-reactive power C k , i , t , s S V C Operation cost of static var compensator (SVC)
L i j , m a x Maximum threshold of the electrical distance U k , i , t , s The nodal voltage magnitudes connected at node i and time t within cluster k under scenario s
δ i , j A 0–1 variable equals 1 if Node 1 and Node 2 are in the same cluster
a 1 ,   a 2 The module degree index and the flexible balance contribution index, respectively G k , i j , s ,   B k , i j , s ,   δ k , i j , s Conductance, susceptance, and phase angle of line ij within cluster k under scenario s
N k The k th distributed resource cluster F k , t , s d e f i , u p ,   F k , t , s d e f i , d n Upgraded flexibility deficiency and downgraded flexibility deficiency at time t with cluster k under scenario s
N s Total number of scenarios P k , i , t , s E S ,   P k , i , t , s E V ,   P k , i , t , s T C L ,   P k , i , t , s L o a d Actual power of BIES, BIEV, BITC, load demands connected at node i and time t within cluster k under scenario s
λ 1 ,   λ 2 ,   λ 3 Weighting factor Q k , i , t , s S V C ,   Q k , i , t , s E S ,
Q k , i , t , s E V
Reactive power output of SVC, BIES, and BIEV connected at node i and time t within cluster k under scenario s
U k , i , t , m a x ,   U k , i , t , m i n Upper and lower limits on the nodal voltage magnitudes Δ F k , i , t , s u p ,   Δ F k , i , t , s d n Flexibility margin of the negative value connected at node i and time t within cluster k under scenario s

References

  1. IEA. World Energy Outlook 2022:C. Available online: https://iea.blob.core.windows.net/assets/75cd37b8-e50a-4680-bfd7-0424e04a1968/WorldEnergyOutlook2022.pdf (accessed on 10 November 2022).
  2. Xue, Q.W.; Wang, Z.J.; Chen, Q.Y. Multi-objective optimization of building design for life cycle cost and CO2 emissions: A case study of a low-energy residential building in a severe cold climate. Build. Simul. 2022, 15, 83–98. [Google Scholar] [CrossRef]
  3. Bhattarai, B.P.; Cerio Mendaza, I.D.; Myers, K.S.; Bak-Jensen, B.; Paudyal, S. Optimum Aggregation and control of spatially distributed flexible resources in smart grid. IEEE Trans. Smart Grid 2018, 9, 5311–5322. [Google Scholar] [CrossRef]
  4. Tsaousoglou, G.; Sartzetakis, I.; Makris, P.; Efthymiopoulos, N.; Varvarigos, E.; Paterakis, N.G. Flexibility aggregation of temporally coupled resources in real-time balancing markets using machine learning. IEEE Trans. Ind. Informat. 2022, 18, 4342–4351. [Google Scholar] [CrossRef]
  5. Kermani, M.; Adelmanesh, B.; Shirdare, E. Intelligent energy management based on SCADA system in a real microgrid for smart building applications. Renew. Energy 2021, 171, 1115–1127. [Google Scholar] [CrossRef]
  6. Eini, R.; Linkous, L.; Zohrabi, N.; Abdelwahed, S. Smart building management system: Performance specifications and design requirements. J. Build. Eng. 2021, 39, 102222. [Google Scholar] [CrossRef]
  7. Li, C.; Dong, Z.; Li, J.; Li, H. Optimal control strategy of distributed energy storage cluster for prompting renewable energy accomodation in distribution network. Automat. Electron. Power Syst. 2018, 45, 76–83. [Google Scholar]
  8. Hu, W.Q.; Wu, Z.X.; Lv, X.X.; Dinavahi, V. Robust secondary frequency control for virtual synchronous machine-based microgrid cluster using equivalent modeling. IEEE Trans. Smart Grid. 2021, 12, 2879–2889. [Google Scholar] [CrossRef]
  9. Yu, S.; Liu, N.; Zhao, B. Multi-agent Classified Voltage Regulation Method for Photovoltaic User Group. Automat. Electron. Power Syst. 2022, 46, 20–41. [Google Scholar]
  10. Ding, M.; Liu, X.F.; Bi, R.; Hu, D. Method for cluster partition of high-penetration distributed generators based on comprehensive performance index. Automat. Electron. Power Syst. 2018, 42, 47–52. [Google Scholar]
  11. Kyriakou, D.G.; Kanellos, F.D. Optimal Operation of Microgrids Comprising Large Building Prosumers and Plug-in Electric Vehicles Integrated into Active Distribution Networks. Energies 2022, 15, 6182. [Google Scholar] [CrossRef]
  12. Kyriakou, D.G.; Kanellos, F.D. Energy and power management system for microgrids of large-scale building prosumers. IET Energy Syst. Integr. 2023, 5, 228–244. [Google Scholar] [CrossRef]
  13. Farinis, G.Κ.; Kanellos, F.D. Integrated energy management system for microgrids of building prosumers. Electr. Power Syst. Res. 2021, 198, 107357. [Google Scholar] [CrossRef]
  14. Bian, X.Y.; Sun, M.Q.; Dong, L.; Yang, X.W. Distributed source-load coordinated dispatching considering flexible aggregated power. Automat. Electron. Power Syst. 2021, 45, 89–98. [Google Scholar]
  15. Li, Z.H.; Li, T.; Wu, W.C. Minkowski sum based flexibility aggregating method of load dispatching for heat pumps. Automat. Electron. Power Syst. 2019, 43, 14–21. [Google Scholar]
  16. Zhao, W.M.; Huang, H.J.; Zhu, J.Q. Flexible resource cluster response based on inner approximate and constraint space integration. Power Syst. Technol. 2023, 47, 2621–2629. [Google Scholar]
  17. Jian, J.; Li, P.; Ji, H.R.; Bai, L.Q. DLMP-based quantification and analysis method of operational flexibility in flexible distribution networks. IEEE Trans. Sustain. Energy 2022, 13, 2353–2369. [Google Scholar] [CrossRef]
  18. Shao, C.Z.; Ding, Y.; Wang, J.H.; Song, Y.H. Modeling and integration of flexible demand in heat and electricity integrated energy system. IEEE Trans. Sustain. Energy 2018, 19, 361–370. [Google Scholar] [CrossRef]
  19. Majidi, M.; Zare, K. Integration of smart energy hubs in distribution networks under uncertainties and demand response concept. IEEE Trans. Power Syst. 2019, 34, 566–574. [Google Scholar] [CrossRef]
  20. Xue, J.R.; Cao, Y.J.; Shi, X.H.; Zhang, Z. Coordination of multiple flexible resources considering virtual power plants and emergency frequency control. Appl. Sci. 2023, 13, 6390. [Google Scholar] [CrossRef]
  21. Li, W.; Wei, J.; Cao, Y.; Ma, R.; Zhang, H.; Tian, X. A quantitative assessment method for power system flexibility based on probabilistic optimal power flow. Mod. Electr. Power 2023, 40, 1–12. [Google Scholar]
  22. Lin, Z.; Li, H.; Su, Y. Evaluation and expansion planning method of a power system considering flexible carrying capacity. Power Syst. Prot. Control 2021, 49, 46–57. [Google Scholar]
  23. Sun, H.; Fan, X.; Hu, S. Internal and external coordinated bidding strategy of virtual power plants participating in the day ahead power market. Power Syst. Technol. 2022, 46, 1248–1262. [Google Scholar]
  24. Cao, Y.; Zhou, B.; Chung, C.Y.; Shuai, Z.; Hua, Z.; Sun, Y. Dynamic modelling and mutual coordination of electricity and watershed networks for spatio-temporal operational flexibility enhancement under rainy climates. IEEE Trans Smart Grid 2023, 14, 3450–3464. [Google Scholar] [CrossRef]
  25. Zhang, C.; Liu, Q.; Zhou, B.; Chung, C.Y.; Li, J.; Zhu, L.; Shuai, Z. A central limit theorem-based method for DC and AC power flow analysis under interval uncertainty of renewable power generation. IEEE Trans. Sustain. Energy 2023, 14, 563–575. [Google Scholar] [CrossRef]
  26. Wang, X.; Zhang, H.; Zhang, S. Game model of virtual power plant composed of wind power and electric vehicles participating in power market. Autom. Electr. Power Syst. 2019, 43, 155–162. [Google Scholar]
  27. Zhang, W.; Song, J.; Guo, M. Load balancing management strategy of virtual power plant considering electric vehicle charging demand. Autom. Electr. Power Syst. 2022, 46, 118–126. [Google Scholar]
  28. Zhu, Y.W.; Zhang, T.; Ma, Q.S.; Fukuda, H. Thermal performance and optimizing of composite trombe wall with temperature-controlled dc fan in winter. Sustainability 2022, 14, 3080. [Google Scholar] [CrossRef]
  29. Sun, Z.Q.; Si, W.G.; Luo, S.J.; Zhao, J. Real-time demand response strategy of temperature-controlled load for high elastic distribution network. IEEE Access 2021, 9, 69418–69425. [Google Scholar] [CrossRef]
  30. Huang, H.; Nie, S.L.; Jin, L.; Wang, Y.Y.; Dong, J. Optimization of peer-to-peer power trading in a microgrid with distributed pv and battery energy storage systems. Sustainability 2020, 12, 923. [Google Scholar] [CrossRef]
  31. Deng, Y.S.; Jiao, F.S.; Zhang, J.; Li, Z.G. A Short-Term Power Output forecasting model based on correlation analysis and elm-lstm for distributed PV system. J. Electr. Comput. Eng. 2020, 2020, 2051232. [Google Scholar]
  32. Li, J.Y.; Khodayar, M.E.; Wang, J.; Zhou, B. Data-driven distributionally robust co-optimization of P2P energy trading and network operation for interconnected microgrids. IEEE Trans. Smart Grid 2021, 12, 5172–5184. [Google Scholar] [CrossRef]
  33. Dhaifallah, M.; Alaas, Z.; Rezvani, A.; Le, B.N.; Samad, S. Optimal day-ahead economic/emission scheduling of renewable energy resources based microgrid considering demand side management. J. Build. Eng. 2023, 76, 107070. [Google Scholar] [CrossRef]
  34. Pandžić, H.; Marales, J.; Conejo, A.; Kuzle, I. Offering model for a virtual power plant based on stochastic programming. Appl. Energy 2013, 105, 282–292. [Google Scholar] [CrossRef]
  35. Khani, H.; El-Taweel, N.; Farag, H.E.Z. Supervisory scheduling of storage-based hydrogen fueling stations for transportation sector and distributed operating reserve in electricity markets. IEEE Trans. Ind. Inform. 2020, 16, 1529–1538. [Google Scholar] [CrossRef]
Figure 1. Typical modalities of a smart-building system.
Figure 1. Typical modalities of a smart-building system.
Buildings 13 02854 g001
Figure 2. VB model and the grid-connected operation status of BIEV.
Figure 2. VB model and the grid-connected operation status of BIEV.
Buildings 13 02854 g002
Figure 3. VB model principle of TCL in a smart-building system.
Figure 3. VB model principle of TCL in a smart-building system.
Buildings 13 02854 g003
Figure 4. Principle of second-order cone relaxation.
Figure 4. Principle of second-order cone relaxation.
Buildings 13 02854 g004
Figure 5. Flowchart of the proposed method.
Figure 5. Flowchart of the proposed method.
Buildings 13 02854 g005
Figure 6. Modified IEEE 33-bus distribution system with building-integrated flexible DERs.
Figure 6. Modified IEEE 33-bus distribution system with building-integrated flexible DERs.
Buildings 13 02854 g006
Figure 7. Scenario generation and reduction results of PV outputs at node 24: (a) initial scenarios; (b) reduced scenarios.
Figure 7. Scenario generation and reduction results of PV outputs at node 24: (a) initial scenarios; (b) reduced scenarios.
Buildings 13 02854 g007
Figure 8. Building-integrated resource cluster division index function curves.
Figure 8. Building-integrated resource cluster division index function curves.
Buildings 13 02854 g008
Figure 9. Optimal cluster division result for the distribution network with building-integrated flexible DERs.
Figure 9. Optimal cluster division result for the distribution network with building-integrated flexible DERs.
Buildings 13 02854 g009
Figure 10. Flexible regulation capacity of the clusters of the building-integrated flexible DERs.
Figure 10. Flexible regulation capacity of the clusters of the building-integrated flexible DERs.
Buildings 13 02854 g010
Figure 11. Comparative results with Schemes 1–3. (a) The 24 h network loss; (b) The node voltage at noon.
Figure 11. Comparative results with Schemes 1–3. (a) The 24 h network loss; (b) The node voltage at noon.
Buildings 13 02854 g011aBuildings 13 02854 g011b
Table 1. Technical specifications of building-integrated equipment and distribution network.
Table 1. Technical specifications of building-integrated equipment and distribution network.
BIES P c h , i , t , m a x E S   =   300 kW P d i s , i , t , m a x E S   =   300 kW η c h E S   =   0.95
E i , t , m a x E S   =   1200 kWh E i , t , m i m E S   =   200 kWh η d i s E S   =   0.95
BIEV P c h , i , t , m a x E V   =   400 kW P c h , i , t , m a x E V   =   400 kW η c h E V   =   0.95
E i , t 0 E V   =   800 kWh η d i s E V   =   0.95
BITCL A 1   =   8.16   m 2 A 2   =   35   m 2 ρ c V   =   1.25 kJ/°C
η T C L   =   0.80 η W a l l   =   0.85 η C O P   =   0.95
Unit prices λ k P V   =   1.2 USD/kW λ k E S   =   0.5 USD/kW λ k E V   =   0.8 USD/kW
λ k T C L   =   0.8 USD/kW λ k L o s s   =   1.1 USD/kW λ k d e f i   =   0.9 USD/kW
Table 2. Comparative economic performance results.
Table 2. Comparative economic performance results.
SchemeOperation Cost (USD)Network Loss Cost (USD)Flexibility Deficiency Cost (USD)
15697.83089.62608.2
26027.13254.32872.8
36452.63438.13014.5
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, Y.; Hao, S.; Zhang, J.; Yu, L.; Luo, Y.; Zhang, K. Optimal Cluster Scheduling of Active–Reactive Power for Distribution Network Considering Aggregated Flexibility of Heterogeneous Building-Integrated DERs. Buildings 2023, 13, 2854. https://doi.org/10.3390/buildings13112854

AMA Style

Fu Y, Hao S, Zhang J, Yu L, Luo Y, Zhang K. Optimal Cluster Scheduling of Active–Reactive Power for Distribution Network Considering Aggregated Flexibility of Heterogeneous Building-Integrated DERs. Buildings. 2023; 13(11):2854. https://doi.org/10.3390/buildings13112854

Chicago/Turabian Style

Fu, Yu, Shuqing Hao, Junhao Zhang, Liwen Yu, Yuxin Luo, and Kuan Zhang. 2023. "Optimal Cluster Scheduling of Active–Reactive Power for Distribution Network Considering Aggregated Flexibility of Heterogeneous Building-Integrated DERs" Buildings 13, no. 11: 2854. https://doi.org/10.3390/buildings13112854

APA Style

Fu, Y., Hao, S., Zhang, J., Yu, L., Luo, Y., & Zhang, K. (2023). Optimal Cluster Scheduling of Active–Reactive Power for Distribution Network Considering Aggregated Flexibility of Heterogeneous Building-Integrated DERs. Buildings, 13(11), 2854. https://doi.org/10.3390/buildings13112854

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