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

Preprint

Article

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

Altmetrics

Downloads

Views

Comments

A peer-reviewed article of this preprint also exists.

Yu Fu,Shuqing Hao,Junhao Zhang,Liwen Yu,

Yuxin Luo,

Kuan Zhang^*

Yu Fu,Shuqing Hao,Junhao Zhang,Liwen Yu,

Yuxin Luo,

Kuan Zhang^*

This version is not peer-reviewed

Submitted:

17 September 2023

Posted:

18 September 2023

You are already at the latest version

Alerts

Abstract

This paper proposes an active-reactive power collaborative scheduling model with cluster division for the flexible distributed energy resources (DERs) of smart building system to resolve the high complexity of 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 penalty cost for the flexibility deficiency, and the second-order cone-based branch flow method is utilized to convert the power flow equations into the linearized cone constraints, reducing the nonlinearity and heavy computation burden of the scheduling model. Customized virtual battery models for the 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 flexible balance contribution index is formulated for cluster division to achieve the information exchange and energy interaction in the distribution network with high proportion of building-integrated flexible DERs. Comparative studies have demonstrated the superior performance of the proposed methodology on the economic merits and voltage regulation.

Keywords:

Subject: Engineering - Electrical and Electronic Engineering

1. Introduction

Construction sector is the third largest energy consumer, accounting for 33% of the global energy demand and 16% of global CO₂ 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 will be increased by 40% in 2030 and even 50% in 2050 [2]. The aggregation of flexible DERs in the buildings is recognized as a promising alternative method to reduce CO₂ emissions and energy consumption, promoting the profound integration of distributed energy resources (DERs) and building energy supply [3]. The smart building system equipped with power access points will achieve the flexible resources aggregation and exploit the regulation capacity of these resources sufficiently, intensifying the interdependency between the buildings and electricity network [3,4]. Nevertheless, the location of building-integrated photovoltaic (BIPV), building-integrated electric vehicle (BIEV), building-integrated energy storage (BIES), building-integrated temperature-controlled load (BITC) and other comprehensive energy production and storage equipment in the smart building system is 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 distribution network with cluster division for flexible resources integrated into the smart building system will be ideal for the high-penetrated deployment of DERs to meet the carbon-neutral targets.

One typical method studied in [7] for the 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. A chromosome encoding method considering modular index based on the electrical designed to satisfy the demand of distribution network programmer [8]. A K-means clustering algorithm based on the variation sensitivity of voltage-to-power between nodes to solve the voltage overlimit and voltage fluctuations [9]. In addition, in order to tackle the difficulty in the uncertainty and centralized control, a multi-time rolling dynamic cluster division method with the equilibrium degree of reactive power for active distribution network [10]. However, the cluster division in these studies only involve a single index, without considering the comprehensive indexes including the flexibility regulation capacity and electrical connection.

Integration and quantification for the regulation capacity of flexible resources is regarded as an effective means to take advantage of the flexibility of the smart building system [4,11,12,13,14,15]. Lots of inspiring works have been reported on the virtual battery model [11], virtual synchronous machine model [12], flexible supply and demand balance model [13], and node power model [14] for the aggregation of flexible resources. A generic and scalable approach considering zonotopic sets for flexible energy systems is studied in [15] to quantitatively describe their flexibility. A constraint space superposition method based on Minkowski Sum is designed to seek the maximal inner cube constraints for the virtual synchronous machine of the heat pump cluster [16]. Moreover, a day-ahead optimal scheduling model based on the generalized virtual battery model is developed to determine the energy schedule and reserve capacity of electric vehicle and air conditioning [17]. After the flexible resources aggregation, the coordination and optimization of 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 the building-integrated flexible DERs for distribution networks are not involved.

In this paper, the customized virtual battery models are proposed for building-integrated flexible DERs to aggregate the power characteristics of flexible resources while quantifying the regulation capacity of flexible resources with time-shifting power and energy boundaries. Moreover, a cluster division algorithm considering the structural and functional cluster division indexes for building-integrated flexible DERs is developed to automatically divide the distribution network with high proportion renewables and flexible DERs into the optimal clusters. Furthermore, an optimal active-reactive power collaborative scheduling model for the clusters of the building-integrated flexible DERs is formulated to minimize the operational cost, the power loss cost, and penalty cost for flexibility deficiency while handling the uncertainties from the renewable energy in the distribution network.

2. Flexible regulation capacity of flexible DERs in smart building system

2.1. Smart building system

Figure 1 depicts the typical modality of smart building system among the high renewable penetrated distribution network. The smart building system is generally equipped with cogeneration systems, distributed photovoltaics, electric vehicle, heat pumps, energy storage to satisfy the demand of electricity, cold and heat load in the building [5]. In order to fully utilize the space and time flexibility of the building-integrated flexible DERs, the establishment of 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 the 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 the smart building system (1), consisting of charge power constraints, discharge power constraints, and energy state constraints is developed to mathematically describe the operating characteristic,

P_{ch, i, t, \min}^{ES} \leq P_{ch, i, t}^{ES} \leq P_{ch, i, t, \max}^{ES}

(1a)

P_{dis, i, t, \min}^{ES} \leq P_{dis, i, t}^{ES} \leq P_{dis, i, t, \max}^{ES}

(1b)

E_{i, t, \min}^{ES} \leq E_{i, t}^{ES} \leq E_{i, t, \max}^{ES}

(1c)

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} - \frac{P_{dis, i, t}^{ES}}{η_{dis}^{ES}}) Δ t

(1d)

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 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 coefficient of BIES, and

∆ t

is the time interval.

The flexible regulation capacity of energy storage in the smart building system with energy time-shift capacity and power regulation capacity can be expressed as,

F_{i, t}^{ES, s, up} = P_{i, t, \max}^{ES} - P_{i, t}^{ES}

(2a)

F_{i, t}^{ES, s, dn} = P_{i, t}^{ES} - P_{i, t, \min}^{ES}

(2b)

where

F_{i, t}^{E S, s, u p}

and

F_{i, t}^{E S, s, d n}

denote the upgraded flexible energy supply and 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.2. Virtual battery (VB) model in the smart building system

The VB model can accurately describe the operation and regulation characteristics of flexible DERs by constraining the energy state, energy boundary and power boundary, which can be expressed as,

P_{i, t, \min}^{G} \leq P_{i, t}^{G} \leq P_{i, t, \max}^{G}

(3a)

E_{i, t, \min}^{G} \leq E_{i, t}^{G} \leq E_{i, t, \max}^{G}

(3b)

E_{i, t + 1}^{G} = E_{i, t}^{G} + P_{i, t}^{G} Δ t + Δ E_{i, t}^{G}

(3c)

where

P_{i, t}^{G}

and

E_{i, t}^{G}

are the regulation power and flexibility reserve energy of the flexible DERs in the 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}

E_{i, t, m a x}^{G}

denote the lower and upper bounds of the regulation power and flexibility reserve energy; (3c) represents the state of the flexibility reserve energy of the flexible DERs and

∆ E_{i, t}^{G}

is the effect of other factors on the electric energy of the VB model.

2.3. VB model of EV in the 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} \leq P_{ch, i, t}^{EV} \leq P_{ch, i, t, \max}^{EV}

(4a)

P_{di s, i, t, \min}^{EV} \leq P_{di s, i, t}^{EV} \leq P_{di s, i, t, \max}^{EV}

(4b)

E_{i, t_{0}}^{EV} - P_{dis, i, t, \max}^{EV} (t - t_{0}) \leq E_{i, t}^{EV} \leq E_{i, t_{0}}^{EV} + P_{ch, i, t, \max}^{EV} (t - t_{0})

(4c)

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} - \frac{P_{dis, i, t}^{EV}}{η_{dis}^{EV}}) Δ t

(4d)

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 threshold 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}

E_{i, t}^{E V}

denote the 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 coefficient of BIEV, (4c) and (4d) represent the boundary of operation capacity and operation state.

Then, considering the power boundary

P_{i, t, m a x}^{E V}

and

P_{i, t, m i n}^{E V}

of BIEV is affected by the operation capacity, the upgraded flexible energy supply

F_{i, t}^{E V, s, u p}

and downgraded flexible energy supply

F_{i, t}^{E V, s, d n}

of EV in the smart building system can be expressed as,

F_{i, t}^{EV, s, up} = \min {P_{i, t, \max}^{EV} - P_{i, t}^{EV}, \frac{E_{i, t, \max}^{EV} - E_{i, t}^{EV}}{Δ t}}

(5a)

F_{i, t}^{EV, s, dn} = \min {P_{i, t}^{EV} - P_{i, t, \min}^{EV} P, \frac{E_{i, t}^{EV} - E_{i, t, \min}^{EV}}{Δ t}}

(5b)

2.4. VB model of TCL in the smart building system

The equivalent thermal parameter model is adopted to describe the thermal dynamics characteristics of TCL, which equivalents internal environment, external environment, and heating (cooling) capacity to circuit devices to analyze the dynamic change relationship of indoor temperature and power [18,19]. 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}

. When the ambient temperature

T_{o u t, i, t}^{T C L}

is higher than the indoor temperature

T_{i n, i, t}^{T C L}

, the temperature of the heating equipment

T_{h e a t, i, t}^{T C L}

and the ambient temperature combine to raise the indoor temperature, increasing the energy of the VB model, as follows,

{\begin{cases} 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} = \frac{P_{i, t}^{TCL} Δ t}{Q_{i, t}^{1}} \end{cases}

(6a)

P_{i, t, \min}^{TCL} \leq P_{i, t}^{TCL} \leq P_{i, t, \max}^{TCL}

(6b)

T^{base} - σ \leq T_{in, i, t}^{TCL} \leq T^{base} + σ

(6c)

T_{in, i, t + 1}^{TCL} = (1 - \frac{A_{2} η^{Wall}}{ρ c V}) T_{in, i, t}^{TCL} + \frac{A_{2} η^{Wall}}{ρ c V} T_{out, i, t}^{TCL} + \frac{P_{i, t}^{TCL}}{ρ c V η^{COP}} Δ t

(6d)

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

V

are the indoor air density, air heat capacity, and room volume. (6b) and (6c) represent the boundary of actual power and the indoor temperature, where

T^{b a s e}

and

σ

are the indoor temperature baseline and dead zone value.

The upgraded flexible energy supply

F_{i, t}^{T C L, s, u p}

and downgraded flexible energy supply

F_{i, t}^{T C L, s, d n}

of TCL in the smart building system can be expressed as,

F_{i, t}^{TCL, s, up} = P_{i, t, \max}^{TCL} - P_{i, t}^{TCL}

(7a)

F_{i, t}^{TCL, s, dn} = P_{i, t}^{TCL} - P_{i, t, \min}^{TCL}

(7b)

3. Active-reactive power collaborative optimization scheduling of distribution network cluster with the aggregation of flexible DERs in smart building system

3.1. Cluster division index of building-integrated flexible DERs

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

ξ = \frac{1}{2 m} \sum_{i} \sum_{j} (μ_{i j} - \frac{r_{i} r_{j}}{2 m}) δ (i, j)

(8a)

{\begin{cases} μ_{i j} = 1 - \frac{L_{i j}}{L_{i j, \max}} \\ m = \sum_{i} \sum_{j} μ_{i j} \\ r_{i} = \sum_{i} μ_{i j} \end{cases}

(8b)

L_{i j} = \sqrt{{(d_{i 1} - d_{j 1})}^{2} + {(d_{i 2} - d_{j 2})}^{2} + ... + {(d_{i n} - d_{j n})}^{2}}

(8c)

d_{i j} = (S_{i i}^{P} - S_{i j}^{P}) + (S_{i i}^{Q} - S_{i j}^{Q})

(8d)

Δ U = S^{P} Δ P + S^{Q} Δ Q

(8e)

where

μ_{i j}

m

, and

r_{i}

denote the edge weight connecting node

i

and

j

, sum of the edge weight in the network, edge weight connecting node

i

;

δ (i, j)

, a 0-1 variable, equals to 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 of node

i

and

j

as shown in (8c) and (8d);

L_{i j, m a x}

is the maximum threshold of the electrical distance; (8e) indicates the sensitivity relationship between node voltage and node injected power;

∆ U

∆ P

, and

∆ Q

are the variation of the 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 the building-integrated flexible DERs refers to its flexible energy supply capacity to satisfy the flexible demands of net load in the smart building system, which can be divided into upgraded flexible energy supply capacity and downgraded flexible energy supply capacity. Flexible balance contribution index

φ

can measure the matching degree between distributed power consumption demands and the supply of the building-integrated flexible DERs through calculating the equilibrium proportion of the flexible DERs bearing the node flexibility demands at a certain moment,

φ = \frac{1}{n} (\sum_{i = 1}^{n} α_{i, t}^{G, up} + \sum_{i = 1}^{n} α_{i, t}^{G, dn})

(9a)

{\begin{cases} α_{i, t}^{G, up} = \frac{F_{i, t}^{G, s, up}}{F_{i, t}^{dem}} \\ α_{i, t}^{G, dn} = \frac{F_{i, t}^{G, s, dn}}{F_{i, t}^{dem}} \end{cases} G \in {ES, EV, TCL}

(9b)

F_{i, t}^{dem} = P_{i, t}^{Load} - P_{i, t}^{PV}

(9c)

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 the smart building system, deriving from the difference between the load demands

P_{i, t}^{L o a d}

and BIPV generation.

Then, the cluster division index

γ

of the 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 formula (8a) and formula (9a),

γ = a_{1} ξ + a_{2} φ

(10)

The critical procedures of the cluster division process of the 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 the building-integrated flexible DERs based on formula (10);
Node $j$ is randomly selected from the remaining other 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 repeating the second procedure. The division procedures will be stopped until the nodes in the network cannot merge and the cluster division index $γ$ of the 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 the building-integrated flexible DERs to handle the uncertainties from the renewable energy. Specifically, each cluster

k

aims to minimize the operational cost

f_{k}^{1}

, power loss cost

f_{k}^{2}

, and penalty cost for flexibility deficiency

f_{k}^{3}

, as shown in (11a), (12a), and (13a).

f_{k}^{1}

is composed of the maintenance cost

C_{k, i, t}^{P V}

of BIPVs due to the renewable generation curtailment, lifetime degradation cost

C_{k, i, t}^{E S}

of BIESs, compensation cost

C_{k, i, t}^{C L}

of BIEVs and BITCs stemmed from involving power regulation, and the operation cost

C_{k, i, t}^{S V C}

of static var compensator (SVC) from the uncertainties considering the renewable generation.

f_{k}^{2}

can be derived from the power flow calculation during the scheduling circulation, where

U_{k, i, t}

G_{k, i, t}

B_{k, i, t}

δ_{k, i j}

represent the nodal voltage magnitudes, conductance, susceptance, and phase angle of line

i j

at time

t

within cluster

k

;

λ_{k}^{L o s s}

is the compensation price of network loss. The penalty cost

f_{k}^{3}

for flexibility deficiency is introduced when the flexibility supply cannot satisfy the load demands within the clusters, including the upgraded flexibility deficiency

F_{k, t}^{d e f i, u p}

and downgraded flexibility deficiency

F_{k, t}^{d e f i, d n}

, as follows,

f_{k}^{1} = \sum_{t \in Γ} (\sum_{i \in Ω^{PV}}^{N_{k}} C_{k, i, t}^{PV} + \sum_{i \in Ω^{ES}}^{N_{k}} C_{k, i, t}^{ES} + \sum_{i \in Ω^{CL}}^{N_{k}} C_{k, i, t}^{CL} + \sum_{i \in Ω^{SVC}}^{N_{k}} C_{k, i, t}^{SVC}) k = 1, 2, 3..., K

(11a)

C_{k, i, t}^{PV} = λ_{k}^{PV} P_{k, i, t}^{PV} = λ_{k}^{PV} (P_{k, i, t}^{PV 0} - Δ P_{k, i, t}^{PV})

(11b)

C_{k, i, t}^{ES} = λ_{k}^{ES} | P_{k, i, t}^{ES} |

(11c)

C_{k, i, t}^{CL} = λ_{k}^{EV} | P_{k, i, t}^{EV} | + λ_{k}^{TCL} | P_{k, i, t}^{TCL} |

(11d)

C_{k, i, t}^{SVC} = λ_{k}^{SVC} | Q_{k, i, t}^{SVC} |

(11e)

f_{k}^{2} = \sum_{t \in Γ} (λ_{k}^{Loss} \sum_{(i, j) \in N_{k}} G_{k, i j} ({(U_{k, i, t})}^{2} + {(U_{k, j, t})}^{2} - 2 U_{k, i, t} U_{k, j, t} \cos δ_{k, i j})) k = 1, 2, 3..., K

(12a)

f_{k}^{3} = \sum_{t \in Γ} (λ_{k}^{defi} (F_{k, t}^{defi, up} + F_{k, t}^{defi, dn})) k = 1, 2, 3..., K

(13a)

F_{k, t}^{defi, up} = \sum_{i \in N_{k}} Δ F_{k, i, t}^{up} = {\begin{cases} \sum_{i \in N_{k}} | F_{k, i, t}^{G, s, up} - F_{k, i, t}^{dem} | Δ F_{k, i, t}^{up} \leq 0 \\ 0 Δ F_{k, i, t}^{up} \geq 0 \end{cases}

(13b)

F_{k, t}^{defi, dn} = \sum_{i \in N_{k}} Δ F_{k, i, t}^{dn} = {\begin{cases} \sum_{i \in N_{k}} | F_{k, i, t}^{G, s, dn} - F_{k, i, t}^{dem} | Δ F_{k, i, t}^{dn} \leq 0 \\ 0 Δ F_{k, i, t}^{dn} \geq 0 \end{cases}

(13c)

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 total number of divided clusters;

λ_{k}^{P V}

λ_{k}^{E S}

λ_{k}^{E V}

λ_{k}^{T C L}

, and

λ_{k}^{S V C}

represent unit operation prices of BIPV, BIES, BIEV, BITC, and SVC connected at node

i

and time

t

within cluster

k

;

P_{k, i, t}^{P V}

P_{k, i, t}^{P V 0}

, and

Δ P_{k, i, t}^{P V}

denote the actual generation, predicted generation, and generation curtailment of BIPV. Due to the frequent charging/discharging behaviors of the building-integrated flexible DERs, the aging of these resources is inevitable, where

P_{k, i, t}^{E S}

P_{k, i, t}^{E V}

P_{k, i, t}^{T C L}

represent the actual power of BIES, BIEV, BITC connected at node

i

and time

t

within cluster

k

;

Q_{k, i, t}^{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}^{u p}

and

Δ F_{k, i, t}^{d n}

of the negative value, and

λ_{k}^{d e f i}

is the penalty factor for the flexibility deficiency.

Considering the 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}

(14)

3.2.2. Operation constraints

Newton-Laphson method is adopted to describe the 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}

{\begin{cases} Δ P_{k, i, t} - U_{k, i, t} \sum_{(i, j) \in N_{k}} U_{k, j, t} (G_{k, i j} \cos δ_{k, i j} + B_{k, i j} \sin δ_{k, i j}) = 0 \\ Δ Q_{k, i, t} - U_{k, i, t} \sum_{(i, j) \in N_{k}} U_{k, j, t} (G_{k, i j} \sin δ_{k, i j} - B_{k, i j} \cos δ_{k, i j}) = 0 \end{cases}

(15a)

{\begin{cases} Δ P_{k, i, t} = P_{k,, i, t}^{PV} + P_{k, i, t}^{ES} - P_{k, i, t}^{EV} - P_{k, i, t}^{TCL} - P_{k, i, t}^{Load} \\ Δ Q_{k, i, t} = Q_{k, i, t}^{PV} + Q_{k, i, t}^{ES} + Q_{k, i, t}^{SVC} \end{cases}

(15b)

U_{k, i, t, \min} \leq U_{k, i, t} \leq U_{k, i, t, \max}

(15c)

The power output curtailment of BIPV is supposed to be lower than the specified maximum thresholds of the active and reactive power,

Δ P_{k, i, t, m a x}^{P V}

and

Δ Q_{k, i, t, m a x}^{P V}

. Moreover, sufficient power factor

c o s θ

for the operation of photovoltaic inverter is needed to limit the reactive power circulation, where the amount of reactive power injected or absorbed is constrained by the bound (16a) and (16b),

（ P_{k, i, t}^{PV} ）^{2} + （ Q_{k, i, t}^{PV} ）^{2} \leq （ S_{k, i, t}^{PV} ）^{2}

(16a)

\cos θ \sqrt{（ P_{k, i, t}^{PV} ）^{2} + （ Q_{k, i, t}^{PV} ）^{2}} \leq P_{k, i, t}^{PV}

(16b)

0 \leq Δ P_{k, i, t}^{PV} \leq Δ P_{k, i, t, \max}^{PV}

(16c)

0 \leq Δ Q_{k, i, t}^{PV} \leq Δ Q_{k, i, t, \max}^{PV}

(16d)

Also, the state 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} \leq Q_{k, i, t}^{SVC} \leq Q_{k, i, t, \max}^{SVC}

(17)

In addition, the sum of the regulation capacity of the 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}^{F l e x}

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}^{Flex} = P_{k, t}^{ES} \oplus P_{k, t}^{EV} \oplus P_{k, t}^{TCL} = \sum_{i \in N_{k}} P_{k, i, t}^{ES} + \sum_{i \in N_{k}} P_{k, i, t}^{EV} + \sum_{i \in N_{k}} P_{k, i, t}^{TCL}

(18a)

P_{k, t, \min}^{Flex} \leq P_{k, t}^{Flex} \leq P_{k, t, \max}^{Flex}

(18b)

{\begin{cases} P_{k, t, \max}^{Flex} = \sum_{i \in N_{k}} P_{i, t, \max}^{ES} + \sum_{i \in N_{k}} P_{i, t, \max}^{EV} + \sum_{i \in N_{k}} P_{i, t, \max}^{TCL} \\ P_{k, t, \min}^{Flex} = \sum_{i \in N_{k}} P_{i, t, \min}^{ES} + \sum_{i \in N_{k}} P_{i, t, \min}^{EV} + \sum_{i \in N_{k}} P_{i, t, \min}^{TCL} \end{cases}

(18c)

3.2.3. Solution methodology

Due to the nonlinearity and nonconvexity of the power flow equations increasing the heavy computation burden of 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}

Y_{k, i, t}

Z_{k, i, t}

and cone constraints, as shown in (19),

X_{k, i, t} = {(U_{k, i, t})}^{2}

(19a)

Y_{k, i, t} = U_{k, i, t} U_{k, j, t} \cos δ_{k, i j}

(19b)

Z_{k, i, t} = U_{k, i, t} U_{k, j, t} \sin δ_{k, i j}

(19c)

\sqrt{{(Y_{k, i, t})}^{2} + {(Z_{k, i, t})}^{2} + {(\frac{X_{k, i, t} - X_{k, j, t}}{2})}^{2}} \leq \frac{X_{k, i, t} + X_{k, j, t}}{4}

(19d)

Then the power loss cost function (12a), the power flow equations (15a), and the nodal voltage magnitudes constraint (15c) can be reformulated as follows,

f_{k}^{2} = \sum_{t \in Γ} (ρ_{k}^{Loss} \sum_{(i, j) \in N_{k}} G_{k, i j} (X_{k, i, t} + X_{k, j, t} - 2 Y_{k, i, t})) k = 1, 2, 3..., K

(20)

{\begin{cases} Δ P_{k, i, t} - (G_{k . i i} X_{k, i, t} + \sum_{\begin{array}{l} (i, j) \in N_{k} \\ j \neq 1 \end{array}}^{n} (G_{k, i j} Y_{k, i, t} + B_{k, i j} Z_{k, i, t})) = 0 \\ Δ Q_{k, i, t} - (- B_{k, i i} X_{k, i, t} - \sum_{\begin{array}{l} (i, j) \in N_{k} \\ j \neq 1 \end{array}}^{n} (B_{k, i j} Y_{k, i, t} - G_{k, i j} Z_{k, i, t})) = 0 \end{cases}

(21)

{(U_{k, i, t, \min})}^{2} \leq X_{k, i, t} \leq {(U_{k, i, t, \max})}^{2}

(22)

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 4. The installed capacities of BIPVs in nodes 1, 6, 13, 17, 18, 24, 29 are set as 0.5MW, 1.0MW, 0.2MW, 1.5MW, 0.8MW, 1.2MW, 0.4MW to simulate the diverse PV endowments. Figure 5 illustrates the predicted energy outputs of BIPVs derived from the actual values in [6,20,21] and the estimated 24-h loads with a system standard load of 3.175+j2.3MVA [22,23]. The prediction errors of BIPVs generation and node loads follow Gaussian distributions with zero means and 10% standard deviations of the prediction values [21,22,24]. The detailed technical specifications of building-integrated equipment and distribution network are summarized in Table 1 [3,4,24].

The appropriate weights for the cluster division index of the building-integrated DERs and the multi-object function are set as 0.41, 0.29; 0.63, 0.48, and 0.25 after a large number of simulation tests, respectively. The optimization scheduling of distribution network with cluster division for flexible DERs integrated into the smart building system is implemented on an hourly basis over a 24-h horizon, and solved by the CPLEX Optimizer on MATLAB platform.

4.2. Comparative Results and Analysis

Three comparative schemes are performed for in depth analysis on the effectiveness of the proposed methodology.

Scheme 1 performs the proposed optimization scheduling of distribution network with cluster division for flexible DERs integrated into the smart building system in Section 3.
Scheme 2 adopts the centralized optimization scheduling of distribution network without considering the cluster division.
Scheme 3 is the initial distribution network before the optimization scheduling.

4.2.1. Cluster division of the building-integrated flexible DERs

Consider the differences in node voltage amplitude, power injection, and load demands at different times, the detailed data of each node at 12:00 noon is exploited for cluster division. Figure 6 depicts the function curves for the cluster division index of the 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 7,

4.2.2. Flexible regulation capacity results of the building-integrated flexible DERs

Figure 8 illustrates the flexible regulation capacity of the building-integrated flexible DERs based on the flexibility constraint space superimposed by Minkowski sum. It can be observed from Figure 8 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 in cluster 1. In addition, the sum of the total amount of BIPV generation and the regulation capacity of the building-integrated flexible DERs can satisfy the total load demands within the cluster 2, indicating that cluster 2 achieves the complete self-consistency of the active power. On the contrary, the total load demands within the cluster 1 is greater than the sum of the total amount of BIPV generation and the regulation capacity of the building-integrated flexible DERs at time

t = 4

and

t = 20

, thus give rise to the flexibility power deficiency.

4.2.3. Optimization scheduling results with scheme 1-3

The power loss and the node voltage amplitude with schemes 1-3 are illustrated in Figure 9, where the node voltage at 12:00 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 BIPVs 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, remaining the node voltage overlimit. However, since the 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 schemes1 and 2, including the total cost, operational cost, power loss cost, and flexibility deficiency cost. Compared to scheme 2, scheme 1 considering 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 schemes 3 on the economic merits and voltage regulation voltage regulation.

5. Conclusions

In this paper, an optimal active-reactive power collaborative scheduling model for the building flexible DER clusters is proposed to resolve the uncertainties from the renewable energy in 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 voltage regulation than the centralized optimization scheduling of 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 in demand side, further on research will focus on the collaborate flexibility regulation capacity of building-integrated flexible DERs in the supply side and demand side.

Author Contributions

Conceptualization, Y. Fu and K. Zhang; methodology, S. Hao; software, J. Zhang; validation, L. Yu, Y. Fu and K. Zhang; formal analysis, S. Hao; investigation, Y. Luo; resources, J. Zhang; data curation, L. Yu; writing—original draft preparation, Y. Luo; writing—review and editing, Y. Luo and K. Zhang; visualization, Y. Luo; supervision, K. Zhang; project administration, Y. Fu; funding acquisition, Y. Fu, S. Hao, J. Zhang and L. Yu. All authors have read and agreed to the published version of the manuscript.

Funding

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

Conflicts of Interest

The authors declare no conflict of interest.

References

IEA, World Energy Outlook 2022: An updated roadmap to Net Zero Emissions by 2050. Available online: https://iea.blob.core.windows.net/assets/75cd37b8-e50a-4680-bfd7-0424e 04a1968/WorldEnergyOutlook2022.pdf.
Xue, Q.W.; Wang, Z.J.; Chen, Q.Y. Multi-objective optimization of building design for life cycle cost and CO₂ emissions: A case study of a low-energy residential building in a severe cold climate. Build. Simul. 2022, 15, 83–98. [Google Scholar] [CrossRef]
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]
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]
Kermani, M.; Adelmanesh, B.; Shirdare, E. Intelligent energy management based on SCADA system in a real microgrid for smart building applications. Renew. Energ. 2021, 171, 1115–1127. [Google Scholar] [CrossRef]
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]
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]
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]
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]
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]
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]
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]
Zhao, W.M.; Huang, H.J.; Zhu, J.Q. Flexible resource cluster response based on inner approximate and constraint space integration. Power Syst. Tenchno. 2023, 47, 2621–2629. [Google Scholar]
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]
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]
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]
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]
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]
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]
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]
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. [Google Scholar]
Li, J.Y.; Khodayar, M.E.; Wang, J.; Zhou, B. Data-driven distributionally robust co-optimization of P₂P energy trading and network operation for interconnected microgrids. IEEE Trans. Smart Grid 2021, 12, 5172–5184. [Google Scholar] [CrossRef]
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]
Khani, H.; El-Taweel, N.; Farag, H.E.Z. Supervisory scheduling of storagebased hydrogen fueling stations for transportation sector and distributed operating reserve in electricity markets. IEEE Trans. Ind. Informat. 2020, 16, 1529–1538. [Google Scholar] [CrossRef]

Figure 1. Typical modalities of smart building system.

Figure 2. VB model and the grid-connected operation status of BIEV.

Figure 3. VB model principle of TCL in the smart building system.

Figure 4. Modified IEEE 33-bus distribution system with the building-integrated flexible DERs.

Figure 5. Prediction curve of BIPVs generation and node loads. (a) The predicted energy outputs of BIPVs generation; (b) The estimated 24-h loads of the distribution network.

Figure 6. Building-integrated resource cluster division index function curves.

Figure 7. Optimal cluster division result for the distribution network with building-integrated flexible DERs.

Figure 8. Flexible regulation capacity of the clusters of the building-integrated flexible DERs.

Figure 9. Comparative results with scheme 1-3. (a) The 24-h network loss; (b) The node voltage at 12:00 noon.

Table 1. Technical specifications of building-integrated equipment and distribution network.

BIES	$P_{c h, i, t, m a x}^{E S}$ =300kW	$P_{d i s, i, t, m a x}^{E S}$ =300kW	$η_{c h}^{E S}$ =0.95
BIES	$E_{i, t, m a x}^{E S}$ =1200kWh	$E_{i, t, m i m}^{E S}$ =200kWh	$η_{d i s}^{E S}$ =0.95
BIEV	$P_{c h, i, t, m a x}^{E V}$ =400kW	$P_{c h, i, t, m a x}^{E V}$ =400kW	$η_{c h}^{E V}$ =0.95
	$E_{i, t_{0}}^{E V}$ =800kWh		$η_{d i s}^{E V}$ =0.95
BITCL	$A_{1}$ =8.16 $m^{2}$	$A_{2}$ =35 $m^{2}$	$ρ c V$ =1.25kJ/℃
	$η^{T C L}$ =0.80	$η^{W a l l}$ =0.85	$η^{C O P}$ =0.95
Unit prices	$λ_{k}^{P V}$ =1.2 $/kW	$λ_{k}^{E S}$ =0.5 $/kW	$λ_{k}^{E V}$ = 0.8 $/kW
	$λ_{k}^{T C L}$ =0.8 $/kW	$λ_{k}^{L o s s}$ =1.1 $/kW	$λ_{k}^{d e f i}$ =0.9 $/kW

Table 2. Comparative economic performance results.

Scheme	Operation cost ($)	Network loss cost ($)	Flexibility deficiency cost ($)
1	5697.8	3089.6	2608.2
2	6027.1	3254.3	2872.8
3	6452.6	3438.1	3014.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.

© 2023 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.

MDPI Initiatives

Important Links

Choose an area of interest and we will send you notifications of new preprints at your preferred frequency.

Disclaimer

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

Abstract

1. Introduction

2. Flexible regulation capacity of flexible DERs in smart building system

2.1. Smart building system

2.2. Battery model of ES in the smart building system

2.2. Virtual battery (VB) model in the smart building system

2.3. VB model of EV in the smart building system

2.4. VB model of TCL in the smart building system

3. Active-reactive power collaborative optimization scheduling of distribution network cluster with the aggregation of flexible DERs in smart building system

3.1. Cluster division index of building-integrated flexible DERs

3.1.1. Module degree index

3.1.2. Flexible balance contribution index

3.2. Active-reactive power collaborative optimization scheduling within clusters

3.2.1. Optimization objective

3.2.2. Operation constraints

3.2.3. Solution methodology

4. Discussion

4.1. System data

4.2. Comparative Results and Analysis

4.2.1. Cluster division of the building-integrated flexible DERs

4.2.2. Flexible regulation capacity results of the building-integrated flexible DERs

4.2.3. Optimization scheduling results with scheme 1-3

5. Conclusions

Author Contributions

Funding

Conflicts of Interest

References

MDPI Initiatives

Important Links

Subscribe