Next Article in Journal
Round-Trip Wireless Charging Infrastructure for Heterogeneous Electric Vehicles on Highways: Modelling and Optimization
Previous Article in Journal
Personalized Collision Avoidance Control for Intelligent Vehicles Based on Driving Characteristics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Mixed Linear Model of a Safety Dispatch Model in an Active Distribution Network for Source–Grid–Load Interactions

Department of Economic Management, North China Electric Power University, Beijing 102206, China
*
Author to whom correspondence should be addressed.
World Electr. Veh. J. 2023, 14(6), 159; https://doi.org/10.3390/wevj14060159
Submission received: 29 April 2023 / Revised: 23 May 2023 / Accepted: 31 May 2023 / Published: 14 June 2023

Abstract

:
There are a large number of plug-and-play loads in an active distribution network, such as EVs(Electric Vehicles), energy storage, solar power, etc. Due to the lack of security control methods for each terminal node and the lack of distributed power voltage control methods, the large number of loads brings a huge challenge to the security of the distribution network. At present, some regional distribution networks dominated by new energy in China have long-standing problems, such as high voltage impact and frequency flickering, which are extremely harmful to electric equipment, and the resulting load-side accidents have brought huge economic losses. Therefore, research on an optimization model of the source–grid–load interaction in the active distribution network considering the safety characteristics, especially the voltage of the system, will help to improve the quality of the grid dispatch. In this paper, the safety limits of the independent operation of a source network loaded on three sides are used as the operating constraints of the system, and the social welfare of the interaction is maximized as the goal. A joint optimization modeling after the independent solution of the three sides is used as the core means, a heuristic algorithm is used to solve the overall optimization of the whole system, a scheduling optimization model that meets the system security goals is constructed, and this model is used to guide the operation strategy of each node in the system. The Lagrangian relaxation factor is introduced for structural optimization, and finally, the simplified 36-node model of the actual power grid is used for verification. The results show that under the goal of ensuring the economy of the system, the system voltage is controlled within the specified range of the safe operation of the system, which can meet the safety needs of the interaction.

1. Introduction

Power distribution systems around the world are experiencing a large-scale deployment of distributed energy (DER), such as renewable and nonrenewable distributed power generation (DG), distributed energy storage systems (DESS), plug-in electric vehicles (PEV), and micro-combined heat and power (CHP) plants [1]. In recent years, due to the encouragement and stimulus of governments around the world, the penetration rate of renewable energy has continued to rise. This rapid onset of DERs has brought about a paradigm shift in the way that conventionally passive distribution networks have been planned and operated [2], leading to the advent of active distribution networks (ADNs). ADNs are essentially driven by advanced information and communication technologies (ICT) and active network management (ANM), which facilitate the DG control and optimization [3], the distributed storage utilization [4], multi-energy integration [5,6], the coordinated control of various system elements, and the demand response programs [7,8].
The conventional DNP is incapable of providing an adequate planning solution for active distribution networks. This is because the DG interconnections transform the distribution network operation in terms of bidirectional power flows, voltage rises and fluctuations, aggravated fault levels, lower power losses, and reliability and stability problems [9]. Moreover, the intermittent generation of the nondispatchable renewable DGs imposes operational uncertainties on the power availability. Numerous other factors including load variability, demand growth, and electricity market prices introduce further operational uncertainties to distribution system planning. Note that operational uncertainty is the level of accuracy with which an input parameter of the planning problem can be forecasted [10]. The compounded effect of these uncertainties can bring about several operational and control challenges: for instance, protection degradation and stability issues [11]. The typical thorough treatment of uncertain input parameters of the planning problem can either lead to computational inefficiency or intractability; hence, this problem demands the application of accurate and efficient uncertainty modelling techniques [12]. Security deals with network responses to all kinds of disturbances, and so is concerned with network dynamics and transient operating conditions, such as voltage instability in the event of major generation and transmission outages [13]. The persistent expansion of smart grids makes the reliability assessment more complex and challenging. Increasing distributed energy resources improves the power supply adequacy due to their inexhaustible nature, but it can reduce the security as a result of power injection and the power flow increase in the lines in the distribution network [14]. To achieve acceptable and optimal reliability, the old-fashioned ‘fit and forget’ strategy can no longer be applied in modern smart networks. The DERs’ output powers must be frequently changed in different time horizons and operating conditions, where the network must be continuously reconfigured, i.e., network reconfiguration, to address such changes when needed [15]. This poses a huge challenge to distribution network security [16].
Currently, source–grid–load interaction has become a better method for active distribution network dispatch [17], but due to the large difference between the power dispatching method and the main grid, the unit generally suffers from a lack of adjustability. Under normal circumstances, there will be many safety issues, such as machine terminal voltage rise, system frequency fluctuations, voltage flickering, and other related power quality issues [18]. Among them, the problem of voltage lifting is particularly serious. In order to maximum the output of the generator, the unit will raise the voltage, leading to the high voltage of the whole grid, and at the same time raising the safety risk [19,20]. At present, some areas in China are mostly new energy sources, such as some small hydropower-rich areas in Sichuan, where this problem is particularly prominent. The problem of burning user-side related equipment due to voltage rise has repeatedly appeared. Therefore, in an active distribution network, the priority goal of source–grid–load interaction should be the safety of the entire grid, ensuring power quality and reducing overall failure rates [21].
This paper is based on the basic calculation and analysis tools of active distribution networks; it analyzes and calculates through the three-sided independent safety operation standard of the source network load, and then uses a heuristic algorithm for joint analysis and calculation to construct an interactive strategy that considers the security characteristics of each node in the entire network. This analysis method takes the voltage stability and operating limits of the whole network as the optimal solution target for the joint calculation of the whole network. According to the node voltage control target issued by the power grid, each node performs independent analysis and calculation to ensure that the voltage fluctuation is fixed within a certain range. The method mainly depends on the internal combination and external coordination of the generator set, and the demand side capacity response coordination of the load node. Finally, a joint operation optimization is carried out through a heuristic algorithm to realize the optimal safety control by the matching of supply and demand.
Taking into account the complexity of the three-sided joint calculation, the Lagrangian relaxation method is used to further optimize the solution process to ensure the feasibility of the calculation. Finally, the model is verified through the simulation calculation of the calculation example.
This algorithm is suitable for active distribution networks with a large number of independent power stations and independent micro-grids, especially for the main body of those with power generation characteristics such as small hydropower, centralized wind farms, and photovoltaic power stations in river basins. The overall applicability is strong, and it can realize the self-safety restraint of the micro-grid in each grid-connected node. The independent optimization model used in this chapter after analysis of node characteristics has achieved good results in the security scheduling strategy of a certain distribution network in Sichuan.

2. Independent Grid Nodes Safety Control Model

2.1. System Safety Dispatch Control Target Model

In order to ensure the overall safety of the system, under the premise that the grid-connected node voltage is the basic constraint, the entire network control target is optimized, and the node voltage target control strategy is issued in consideration of the grid supply and demand forecast. Based on the node compactness matrix, the goal of adjusting the power flow distribution of the active distribution network is achieved by optimizing distributed power output, adjustable unit gears, reactive power compensation devices, and network topology reconstruction, thereby achieving the lowest overall network loss and the most secure grid operation dispatching method [22].
The objective function optimization objective is to maximize the system safety index, and the node voltage deviation degree model is used as the final function, which needs to meet the minimum deviation degree constraint:
max t = 1 T [ i = 1 N b r n = 1 N n π ( n ) ( V i n V r , i ) 2 / V r , i 2 ]
where the node voltage has the minimum deviation from the square difference, T is a daily voltage measurement period, Nbr is the number of nodes, π ( n ) is the node running state quantity, V i n is the actual value of the node voltage, and V r , i is the standard voltage value of the node.
In addition to the system power balance, system backup, node power limit, and voltage limit constraints in the unit calculation formula, the constraint function also requires frequency constraints, load adjustment value constraints, energy storage system constraints, energy saving constraints, and gateway area control error (ACE) constraints. Inequality constraints are shown below:
{ F min f ( t ) F max P min , C l P C l ( t ) P max , C l lim A C E ( t ) = 0 t T 0 P E X ( t ) P max , E X 0 Q E X ( t ) Q max , E X
In these equations, f is the frequency value at time t, P C l ( t ) is the power control adjustment value, ACE is the power exchange deviation of the contact line, and PEX/ Q E X are active/reactive variables for contact line exchange.
The energy conservation constraint on the energy storage nodes is as follows [23]:
E s ( t ) = E s ( t 1 ) P s ( t 1 ) Δ t
In the objective function, the loss cost of the node’s topological reconstruction can be expressed by the node’s active power consumption cost and the reactive support cost.
{ C n c = i = 1 N N c α i P N C i α i = C N C j Δ t P N C i = ( P D N i + 1 P R , N C i + 1 ) + ( P R , N C i P D N i ) C Q = j = 1 N Q β j Q j , β j = C Q j Δ t E r l = i = 1 N b r δ P r l i , δ = C r l Δ t P r l i = i = 1 N b r R i ( m , n ) G i ( m , n ) [ V m 2 + V n 2 2 V m V n cos ( θ m θ n ) ]
where C n c is the active loss of the reconstructed node, α i is the number of nodes, C N C j is the adjustment cost of node power, and P N C i , P R , N C i , and P D N i are the value of node power scheduling, the current value of node I, and the control value of node i, respectively. P r l i is the active power loss of node i, Nbr is the total number of branch nodes of the distribution network in the system, and Ri(m,n) and Gi(m,n) are the tightness value and conductance of the m and n branch of node i, respectively.
Finally, the power grid current constraint is as follows:
f l max i = 1 N P i , t G l , i j = 1 N D j , t G l , j f l max
where G is the transfer distribution factor of line l to the power generation nodes, and f is the current line value.

2.2. Combination Optimization Model of Safety Control of Power Generation Nodes

The security control of power generation nodes is based on the node association matrix. In accordance with the optimal goal of ensuring the safety of the system, the model must control the operating limits of the voltage and frequency of the parallel gateway. Independent consideration is given to the unit combination in the microgrid node to deal with system voltage rise and other issues. In the new energy-related node output combination mathematical model for active distribution networks, various constraints need to be considered, and the objective function is mainly formulated based on the system operation limit [24]. Current main functions include voltage limit, frequency limit, gate grid-connected power limit, standby reliability, non-new energy unit pollution emission indicators, etc. In addition, the voltage value is used as the optimal objective function of the distribution network unit combination. Due to the existence of randomly fluctuating power sources in active distribution networks, such as solar sources and wind turbines [25,26], the uncertainty discrete model is established for the load and power sources, and the normal distributed method is adopted when the predicted value is known [27].
{ D ^ t s = D t s ( 1 + e ^ s ) G ^ t s = G ^ t s ( 1 + e ^ s )
Voltage control is directly related to reactive power, so the lowest reactive power loss is used as the solution goal to ensure the constraints of voltage and the power factor. The voltage limit and power factor are used as constraints, and the overall network loss is calculated by dividing the load among 96 periods of a day, which is formulated as follows:
P l o s s , t = min t = 1 96 i 1 n V j , t G i j , t cos θ i j , t C s i
where t is the time period, P l o s s , t is the reactive power loss at moment t, and C s i is the exponential function of the shutdown function T i , t o f f of the system unit at the cut-off time t, which can be represented as:
C s i = U i , t ( 1 U i , t 1 ) [ α i + β i ( 1 exp T i , t o f f τ i ) ]
where U i , t is the operation state of the system at time t.
System constraints are as follows:
System power balance constraints:
i = 1 N P i t + j = 1 N w P w j t + k = 1 N s P s k t = P l t
where P i t is the power of conventional adjustable unit i in period t; P w j t is the power of j, a new energy random output unit, in period t; P s k t is the power during period t of the kth energy storage related facilities; P l t is the sum of system load power; and N, Nw, and Ns are the number of various facilities.
Node power balance constraint:
{ Δ P i , t = P G i , t P D i , t V i , t j = 1 n V j , t ( G i j , t cos θ i j , t + B i j , t sin θ i j , t ) = 0 Δ Q i , t = Q G i , t + Q C i , t Q D i , t V i , t j = 1 n V j , t ( G i j , t sin θ i j , t B i j , t cos θ i j , t ) = 0
where P G i , t and Q G i , t are the active and reactive inputs of unit i at time t, P D i , t and Q D i , t are the active and reactive outputs of unit i at time t, Q C i , t is the reactive input of the reactive compensation device, V is the voltage amplitude, θ is the phase angle, and G and B represent nodal conductance and susceptance.
System back-up capacity constraints:
i N a U i , t P l max P t l + P s max + R t P t w
where U i , t is the operating state of the system at time t, P l m a x is the maximum load power of the system, P t l is the load power of the system at time t, P s m a x is the maximum operating power of energy storage, and R t is the reserve capacity at time t. Consider that the system is an active distribution system, which is a small power system, and this work takes 5% of the maximum load as a constant.
Minimum open/down time constraint for the unit:
{ ( U i , t 1 U i , t ) ( S i , t o n O i o n ) 0 ( U i , t U i , t 1 ) ( S i , t o f f O i o f f ) 0
where S i , t o n and S i , t o f f represent the minimum startup and shutdown time, and O i o n and O i o f f represent the unit startup and shutdown time.
Variable inequality constraint:
{ P i t max { P i min , P i t 1 r i d × Δ t } P i t min { P i max , P i t 1 + r i u × Δ t } P G i , min P G i , t P G i , max Q G i , min Q G i , t Q G i , max 0 Q Ci , t Q Ci , max V i , min V i , t V i , max
{ T l , min T l T l , max t = 1 24 | T l , t T l , t 1 | T max t = 1 24 | C i , t C i , t 1 | C max
In these inequalities, P i max and P i min are the upper and lower output limits of unit i, r i u and r i d are the upper and lower climbing powers of unit i, and P i t and P i t 1 are the load power of unit i at time t and time t − 1, respectively. P G i , t / Q G i , t is the compensation value of the active/reactive power capacitor of unit i at time t, P G i , min and P G i , max     are the maximum and minimum values of the active capacitance compensation value of unit i, Q G i , min and Q G i , max     are the maximum and minimum values of the reactive capacitance compensation value of unit i, Q c i , t is the reactive power input of the reactive power compensation device, and V i , t is the voltage amplitude of unit i at time t. T l , t is the unit gear at time t, T l , max and T l , min are the maximum and minimum values of the unit gear, and   C i , t is the number of capacitor switches of unit i at time t.
Consider the irreducible constraint in the inequality T and C, using the system deformation equation below:
t = 1 n ( x t x t 1 ) F ( x t x t 1 ) X max F ( x t x t 1 ) = { 1   x t > x t 1 1 x t < x t 1
Continuous accessibility functions can be obtained by:
F ( x t x t 1 ) = exp ( ( x t x t 1 ) / T ) exp ( ( x t x t 1 ) / T ) exp ( ( x t x t 1 ) / T ) + exp ( ( x t x t 1 ) / T )
Then, we can find the derivative function of the two functions T and C.
Based on the unit safety constraints for various new energy sources connected to the grid, the discrete hierarchical calculation method is adopted to decompose the overall solution into the node optimal control equation without considering the safety constraints, the node limit analysis of the unit start, or the node control equation for introducing the penalty function.
The mathematical model for the node limit analysis is:
min a A s P s t T [ i N g C i P i , t U i , t + ε l L s R ( χ ) ]
R ( χ ) is the quadratic penalty function, which is given by:
R ( χ ) = ( x a | x a 1 1 2 F ( x a ) x a | x a 1 + 1 2 F ( x a ) )
F ( x a ) is a function that represents the relationship between transformer gear step size and capacitor switching times. It is a non-continuously varying discrete quantity. In the process of optimisation, the adjustable parameters move closer to the centre of the definition interval in their definition domain. A quadratic penalty function is added as follows.
ϕ ( x a ) = 1 2 μ a ( x a x a 1 ) 2
where μ a is the punishment factor.
The penalty variable reflects the restriction of the current calculation on the unit output; when the penalty variable is nonzero, the safety constraint will not be satisfied. The calculation needs to be revised under the nonsafety constraint in the first step, and cyclic calculations need to be performed until for lim R ( χ ) = 0 x N , N is the number of iterative constraints; if the constraint number is exceeded or cannot converge, the approximation value is adopted as the optimal value.
The network constraint is the grid operation constraint. As one of the input data of the source network charge interaction, the three-sided joint calculation is completed.

2.3. Virtual Power Generation Safety Dispatching Model of Load Nodes

At present, the loads participate in the interaction primarily by way of virtual power generation or demand response. Therefore, the highest active node power response mechanism or the lowest cost is the direct goal of flexible control, power quality, such as voltage and harmonics, can be used as a constraint, and system reliability indicators need to be considered [28].
Special carriers, such as energy storage, are defined as flexible reverse loads whose power fluctuates between the maximum capacity and the optimal SOC state to participate in system optimization.
First, we defined the system node response value formula as follows:
V ( t ) = t = 1 T i = 1 N r c C v i , t P r i ( t ) Δ t
where C v i , t is the node response price of node i at moment t and P r i is the adjustable power for the response price of the node at moment t, which is calculated by:
P r i ( t ) = P i , t P i , t 1
The value of a node can be expressed as:
V ( i , t ) = V ( 0 , t ) + C 0 , t [ V ( i , t ) V ( 0 , t ) ] × { 1 + V ( i , t ) V ( 0 , t ) 2 γ ( t ) V ( 0 , t ) }
where C 0 , t is the electricity price before the node period t and γ ( t ) is the node elasticity coefficient at time t, which can be calculated by the node power limit. This is an open-down parabolic function, and all node accumulation requires segmented optimization.
The active reactive and voltage relationships of the nodes are also considered as follows:
V ( i , t ) = V ( 0 , t ) { 1 + i = 1 , i t T γ ( t , i ) V ( i , t ) V ( i , 0 ) + C i , t D ( i , t ) max { D ( i , t ) } V ( i , 0 ) }
At the same time, consider the active and reactive power and voltage relationships of the node as follows:
{ P i = ( α i V i 2 Z N i + β i V i I N i + δ i S N i ) cos ( φ i ) Q i = ( α i V i 2 Z N i + β i V i I N i + δ i S N i ) sin ( φ i )
where α , β , δ are the load-type proportional coefficients, which are calculated from the node. V, Z, I, and S are node voltage, impedance, power, and power, respectively. Through calculation, it can be found that the voltage fluctuation of the node has a proportional effect on the power, that is, the relationship of Y = k X + b . Under the requirements of the system node voltage limit, voltage fluctuations need to be controlled. According to the PQ adjustment value given by the load multi-agent, the instantaneous voltage fluctuation of the adjustment is calculated in the reverse direction, and it needs to be controlled within the voltage limit as one of the main constraint conditions [29]. The probability of load loss is mainly caused by grid faults and the insufficient capacity of the unit. Loss of load will cause the load to be completely unable to participate in the response. The constructed load loss probability formula is.
{ Pr ( λ t ) = { i = 1 N g ( 1 ε i , t ( 1 K i ) ) | ( 1 γ i , t M i ) & ( R i , t i = 1 N l P i , t + s = 1 N s P s , t ( 1 ϕ s o c ) ) } ( 1 P t × L s ( t ) ) P t = i = 1 N G α i , t δ i , t 1 P i , t 1 + i = 1 N G j > i N G α i , t δ i j , t 2 P i j , t 2 δ i j , t 2 = ( P i , t + R i , t ) p i ( P i , t + R i , t ) p i + ( P j , t + R j , t ) p j L s ( t ) = i = 1 N G α i , t δ i , t 1 P i , t 1 ( P i , t + R i , t S R i , t ) + i = 1 N G j > i N G α i , t δ i j , t 2 P i j , t 2 ( P i , t + R i , t + P j , t + R j , t S R t )
where ε i , t and γ i , t are the node fault state of the power grid and the unit; K i is the reliability index value of the power grid, which is generally calculated by the local power grid company according to the annual fault time; and M i is the unplanned shutdown rate of the unit which is an annual statistical value. Additionally, R t is the system capacity, P i , t is the reduced power value for the faulty unit, P s , t is the energy storage power, ϕ s o c is the optimal SOC threshold for the energy storage battery, δ i , t 1 is the contribution factor to unit load loss, and R is the rotating standby capacity.
In this paper, in nonspecial cases, the independent node value is calculated by the reference node electricity price method, and independent node pricing is realized based on the principle of power economic dispatching [30].
The objective function of the system is the maximum load response value, i.e.:
max ( V ( t ) t = 1 T P l l , t C l l Δ t L s ( t ) )
Reliability constraints:
| 1 max { ( P s e t P i ) k p + V i , ( Q s e t Q i ) k q + V i } V s t d | α
where α is the normal voltage operation limit, generally 5–10%; P s e t is the set and adjusts the active power value, and k p is the slope of the active power and voltage (the reactive power is the same).

3. Optimization Model of Mixed Linear Multi-Node Combination

A joint optimization model based on a heuristic algorithm.
The optimal problem is typically a nonlinear programming problem, which contains multiple input variables and solving constraints, with a general mathematical model described as follows:
{ max ( min ) f ( x ) s . t h ( x ) = 0 g min g ( x ) g max
It can also be expressed as a form of P = ( S , H , F ) , where f ( x ) R 1 is a constraint equation, h ( x ) R i is a combination of i equality constraints, and g ( x ) R m is a combination of inequality constraints with m containing upper and lower limits.
Supposing there is a set of contact variables Xi, i =1,2…,n, the problem of combinatorial optimization is the global optimal solution that satisfies the H constraints. For any of s* ∈ S, if f(s*) ≤ f(s), then s* is a value more approximating the global optimal solution and provides at least one global optimal solution for F.
Heuristic algorithms are intuitively or empirically constructed algorithms that give a feasible solution for each instance of the combinatorial optimisation problem to be solved at an acceptable cost (in terms of computational time and space), which deviates from the optimal solution to an extent that cannot generally be anticipated. Therefore, based on the existing “source–grid–load” tripartite independent optimization model, this paper selects a heuristic algorithm as the basic method for joint optimization. Furthermore, by using a self-control model for each of the nodes of the grid, we have already changed the grid optimization problem from a nonlinear problem to a mixed linear problem.
Currently, the relatively common heuristic algorithms are the simulated annealing algorithm (SA), the genetic algorithm (GA), ant colony optimization (ACO), artificial neural networks (ANN), etc.
The ant colony system (ACO) was first proposed by Italian scholars Dorigo and Maniezzo, etc. It was first proposed in the 1990s, and it is an approximate optimization model based on the study of real ant colony behavior. The basic idea of the application of the ant colony algorithm in optimization problems is as follows: the process of the problem to be optimized is represented by the walking path of the ants, and all the path components of the ant colony are regarded as the feasible solutions and space of the problem to be optimized. The ant with the shortest path contains more information. According to the passage of time, the information concentration accumulated by the shorter path gradually increases, and the number of ants selected by the path gradually increases. In the end, the ant colony will concentrate on this optimal path under the effect of positive feedback, and the corresponding solution at this time is the optimal solution of the problem to be optimized. The basic calculation process of ACO is shown in Figure 1.
The basic calculation process of the ACO is:
(1)
Set up multiple ants according to the specific optimization goals, and at the same time make three initial groups for the three sides of the source network load, and search separately;
(2)
Initialize an equal amount of pheromones on each path, as:
τ i j ( 0 ) = C Δ τ i j ( 0 ) = 0
(3)
Take an ant, calculate the transition probability, select the next optimized calculation node according to the roulette method, and update the taboo table [16]. After each ant completes an optimized output, it releases pheromone on the path of the optimized combination. The amount of pheromone is proportional to the quality of the solution. Considering the correlation between nodes, a random local search strategy is adopted. It can optimize the operation of two adjacent nodes to the best one, and the amount of pheromone on the better node is also the largest. Later, the probability of ant selection increases;
(4)
Each ant takes the legal system node optimization path and retains the pheromone amount Δ τ i j k that the ant did not release between the nodes, and the ant dies;
(5)
Repeat steps 3–4 until all ants complete the optimization process;
(6)
Compute pheromone increments Δ τ i j k and pheromones for all optimization modes τ i j ( t + n ) ;
(7)
Record this iteration path and update the current optimum combination to empty the taboo table;
(8)
Reach a predetermined number of iteration steps or stagnation (all ants choose the same path and the solution no longer changes) [17]; the algorithm ends with the current optimal solution as the optimal solution of the problem, otherwise the iteration continues.
The pheromones and transfer probabilities are modeled as follows:
(1)
Transfer probabilistic calculation formula:
p i j k ( t ) = { [ [ τ i j ( t ) ] ] [ η i j ( t ) ] β s J k ( i ) [ τ i s ( t ) ] α [ η i s ( t ) ] β , i f j J k ( i ) 0 , e l e s
Inspired factor:
η i j = 1 d i j
where α and β are two constants, representing the weighted values of pheromones and visibility, respectively.
(2)
Pheromone calculation formula:
τ i j ( t ) = ( 1 ρ ) τ i j + k = 1 m Δ τ i j k Δ τ i j k = { Q L k , i f   a n t   k   g o e s   t h r o u g h   t h e   o p t i m i z a t i o n   n o d e ( i , j )   i n   t h i s   t r a v e r s a l 0 , e l s e
where Q is normal and Lk is all of the ants in this traversal.
(3)
Pheromone update method
τ 0 = ( n L n n ) 1 τ i ( t + 1 ) = ( 1 ρ ) τ i ( t ) + ρ , Δ τ i ( t + 1 ) + ε ( t + 1 ) , ( i = 1 , 2 , 3 , n )
where ρ is the evaporation rate, which is an additional function: ε = e Δ τ i j b , s .
After updating the rules, the relationship can be obtained as follows:
Δ τ i j = s S m [ p ( S m ) s S m | s < i , j > s F ( s ) ]
The combined problem p, which handles the transfer probability in traditional ACO, increases the difficulty of solving and uses the continuous probability distribution function instead as a solution, and the three-dimensional independent optimization result adopts the high-dimensional normal distribution, namely the multi-dimensional Gaussian function.
The multidimensional Gaussian function formula is:
Ν ( x | μ , Σ ) = 1 ( 2 π ) D / 2 1 | Σ | 1 / 2 exp { 1 2 ( x μ ) T Σ 1 ( x μ ) }
For the assumed dataset X = (x1,…,xn)T, each sample is obtained independently from the Gaussian distribution using a maximum likelihood estimation method with its likelihood function of:
ln p ( X | μ , Σ ) = N D 2 ln ( 2 π ) N 2 ln | Σ | 1 2 ( x n μ ) T Σ 1 ( x n μ )
Using matrix derivative properties, we can obtain the Gaussian probability function to obtain the mean μ derivative form:
u ln p ( X | μ , Σ ) = n = 1 N Σ 1 ( x n μ )
With derivative Σ , the following can be obtained:
u ln p | Σ | = ( Σ 1 ) T = Σ 1 Σ i j { ( x n μ ) T Σ 1 ( x n μ ) } = ( x n μ ) T Σ 1 Σ i j ( x n μ )
Set up:
B = Σ 1 ( x n μ ) B T = { Σ 1 ( x n μ ) } T = ( x n μ ) T Σ 1
Finally, the calculation is made by matrix multiplication and the following is obtained:
μ ln p ( X | μ , Σ ) = N 2 Σ 1 + 1 2 Σ 1 { n = 1 N ( x n μ ) ( x n μ ) T } Σ 1
Taking the multi-dimensional Gaussian ACO function, the probability density function is calculated for each step of the ant solution, and finally the complex weight combination value of multiple Gaussian functions is obtained, i.e.:
N i ( x | μ , Σ ) = k = 1 k ω k N ( x | μ , Σ ) = k = 1 k 1 δ k 2 π e ( k 1 ) 2 2 δ 2 k 2 N ( x | μ , Σ )
where δ is the global optimal solution computational equilibrium parameter, the smaller it is the higher the computational complexity, but the higher the probability of selecting the optimal solution. In each iteration, each ant selects a solution based on the weight in the probability density function, and Ni samples the solution through the Gaussian function, according to its covariance matrix.
For the joint scheduling optimization model in this model, there are three possible system operation states, SG(i,n), SN(i,n) and SL(i,n), representing the operating amount of the i node of the n topological island in the system. The information mainly includes the node power value, voltage value, current, and switching signal volume, and a 3D Gaussian function is used as the density function. The solution of the model is the security scheduling optimization result with the lowest joint operation cost of the ant on the topological connection path with the equilibrium node as the initial search node in any topological island.

4. Solving the Calculation Example

4.1. Basic Scenario

We use a simplified 36-node system structure as the main case. Due to the confidentiality needs of the power grid, the simplified node parameters are not provided temporarily [18], and the simplified topology graph is shown in Figure 2.
The 36 nodes simplified by this model are used as the basic model of distribution network calculation, including 11 independent power nodes with a total installed capacity of 270.34 MW; 16 independent load nodes with a 197.7 MW capacity for final assembly; and 9 micro-network nodes with an elastic capacity of 91.14 MW and a loading capacity of 84.63 MW. It can be seen that this distribution network is in the basic supply-and-demand balance and full delivery state.
According to this method of dividing the model, the node types within the system are shown in Table 1.
The formed division of distribution network family groups are shown in Table 2, and the calculated correlation tightness matrix is too large in dimension, and thus will not be shown here.
The network contains three voltage levels, 110 kV, 35 kV, and 10 kV, Each user takes different grid voltages according to different characteristics; furthermore, each user interacts with other distribution network supply areas in the regional power grid through 110 kV. The overall supply and demand situation in the system has certain fluctuations along with the season. The supply–demand ratio is calculated by the annual monthly electricity average index. The analysis contains four dimensions: grid power, network load power, online power, and offline power, as shown in Figure 3. The relationship between the four key quantities is:
P g e n + P o f f = P l d + P o n
where P g e n is the power generation amount, P l d is the load power amount, P o f f is the off-grid power amount, and P o n is the on-grid power amount; this formula is also established for the whole grid and the power exchange of the contact line at any time.
In a historical operation situation, an auto part manufacturer burned the relevant production equipment due to high voltage in the network. The voltage of the production line data monitoring points was analyzed in Figure 4, where the average maximum daily voltage fluctuation reached 4.52%, and the maximum voltage fluctuation was 12.68%, seriously exceeding the limit. Therefore, it can be seen that the overall distribution network has had high voltage for a long time, and the terminal voltage of some generating nodes has significantly deviated from the grid-connection safety license.

4.2. Distribution Network Operation Scenario

Based on the 36-node distribution network model, the resilient node devices within it are decomposed in this paper in order to analyse the impact of the combination of security policies within each generation node. The key elastic node parameters are shown in Table 3, with the power generation ranking shown first, and the load ranking second (where 2 × 3.5 means two 3.5 MW installed units).
The overall operation curve of the whole network’s typical day is shown in Figure 5, where the relationship between the four key quantities is:
P l d + P o u t = P r e + P t r
where P o u t is the external power, P l d is the load power, P r e is the new energy power generation (including scenery power), and P t r is the traditional power generation (water and fire).
Initial generation input power can be seen in Table 4.
In the optimization analysis, each node prediction is used to superposition load predictions and power generation predictions on the full network. According to the outsourcing plan of the power grid (the 96 points of the optimization data, as shown in Figure 5), as the basic setting of the distribution network solution, new energy is designed according to the principle of priority power generation, and gas and biomass are designed as power regulation power supplies.
The first three pictures in Figure 6 show the network summary data load, power generation, and outgoing data; the last three pictures show the power generation and a single day’s predictions for the typical load of elastic users, including wind power, using a super short-term 5 min forecast for a cycle, runoff small hydropower, adopting a super short-term 15 min for a cycle, and a user, adopting a super short-term 15 min forecast for a cycle. Other users of the whole network are not detailed here.

4.3. Model Calculation Results

4.3.1. Grid Voltage Control Target Value

The 96-point voltage control target of the 36-node network is shown in Table 5 and Figure 7, showing the voltage and phase angles of the grid at 0:00. All are displayed in unitary units.
It can be seen from the voltage control target curve that, due to the uncontrollable power supply, to ensure the system supply, the highest voltage point has been relaxed to nearly 1.07, while the minimum voltage of the whole network is about 0.95, and the overall day solution target is 0.063752189, which fully meets the voltage control target of the power grid.

4.3.2. Optimization Result of Safety Control of Power Generation Nodes

All family combinations of units in the network are shown in Figure 7. The shutdown status identification is −1 and the boot identification is 1, to be distinguished by the line of x = 0.
In Figure 8, the first figure shows the output curve at the same time to show the cascade scheduling combination and the voltage fluctuation diagram under the combination of each family unit.
The analysis of each relevant unit combination is as follows:
The role of energy storage in the power generation unit can be intuitively displayed from the diagram, especially with wind power and photovoltaic outputs, as a power regulation device with a high enabling frequency, which can effectively make up for the characteristics of scenery fluctuations;
Nodes 4, 16, 18, and 29 are connected with unadjustable hydropower and photovoltaic power; of these, photovoltaic power has reactive compensation, so all start-up strategies are adopted. Because the families in which the small hydropower is located all have other adjustable power sources, the results for the small hydropower are all units fully open, the same as for node 4;
Nodes 19 and 26 are high-controllable generator sets, and the basic unit combination strategies are the same. It can be found from the unit combinations that the combination of power outputs is not the minimum open number combination completed by the algebraic superposition mode;
The optimization combination strategy of a storage capacity adjustable hydropower unit is relatively complex. Since most of them are cascading hydropower combination models with storage capacity plus runoff, the active and reactive situation of their overall outputs and the utilization of water conservancy resources should be considered. Therefore, adopt all power-up strategies for cascading hydropower to limit the power generation and drainage to achieve the voltage control target;
Under the reactive control of the combined voltage of independent units, the time limit of current trend analysis will occur even with the target constraints. In particular, the randomness of the hydropower + wind power combination is always high, resulting from the reactive output of two types of units. Therefore, it is seen that the reactive control of an ADN cannot completely rely on the power generation side, and requires the coordination of load-side and power grid reactive devices.

4.3.3. Optimization Result of Load Node Safety Control

The virtual power generation combination of whole network load flexible resources in the ideal state is shown in Figure 9; this is the output result under the unified adjustment mode of the whole network load, which is not consistent with the actual output according to family status, as shown in Figure 10. It is obvious that the number of load nodes involved is more than in the ideal state.

4.3.4. Combined Optimization

According to the combination of nodes, the initial point section, and the N-1 safety analysis, we obtained the results shown in Figure 11; N-1 shows the N-1 scene number and the vertical coordinates and the gray scale indicates the load rate of each branch in different N-1 scenarios. The overall visible voltage control condition is good between 0.98 and 1.07, but there is a serious branch limit. The most serious branch, 13, reaches 1.349, and cannot meet the system safety operation requirements.
The initial point current flow and N-1 safety analysis after the ant colony optimization whole network solution is shown in Figure 12. When the voltage is in a stable state, the overload problem is completely eliminated and meets the operational requirements of the power grid to optimize the load; the maximum branch load is 0.9487. There is still some room for optimization. The overall joint optimization results solved by the heuristic algorithm are shown in Figure 13.
After the optimization and solution of the whole network, the new energy combination of the power generation nodes has not changed. Node 26 is readjusted, and the user-side energy storage devices, heat load, and charging piles are adjusted greatly. Combined with the maximum load response elastic space, the power exchange is adjusted to meet the balance of supply and demand in the network and ensure the maximum demand response value.
The above calculation example’s analysis is as follows:
With the maximum load of 96 points, N-1 is 0.9621. After multiple optimizations and iterative solutions, it is found that the calculation has reached the current optimal value, and the equipment transformation should be conducted for this branch;
The whole network voltage is maintained between 0.98 and 1.04, with good results. The number of changes in the optimization of the power generation combinations is small, mainly because it is the power output end, and the appropriate increase in the reactive compensation devices on the power grid side and the user side is helpful to further stabilize the whole network voltage;
The installed capacity and the inverter power determine the actual start-stop state, which may be quite different from the theoretical calculation;
The independent power supply family in an ADN is about equal to 0 with current technology, especially the distribution grids that rely mainly upon new energy, which needs to be adjusted by conventional power supplies, energy storage, contact lines, etc.

5. Conclusions

This paper mainly analyzes an interactive model of source network load scheduling optimization with operational safety as the main goal. In the safety optimization analysis, considering the different operating characteristics of the distribution network and the transmission network, the method of joint optimization after independent analyses and calculations of the source network and the loads on the three sides is adopted.
(1)
For power generation nodes, the mode of source-source interaction is adopted, the unit combination mode of the output node is considered, the safety control model of the power generation node is constructed, and the operating boundary conditions of the generator set are used as constraints to obtain the operating safety combination of the units in the network.
(2)
For the load node, this paper uses a flexible control model to calculate and analyze the safety scheduling strategy of the load; the model of node voltage deviation in one day is used as the ultimate objective function to solve the safety scheduling of the distribution network. After calculating the independent safety scheduling model, a joint optimization model is established according to a heuristic algorithm model to analyze the overall network security characteristics of the active distribution network, and a Lagrangian relaxation calculation is introduced to ensure that the optimization target is obtained in the active distribution network.
(3)
The model availability and accuracy of this section are verified through the actual grid calculation examples, and it is found via the actual comparison that the power quality in the network can be effectively optimized and adapted to the optimization analysis of the active distribution grid of various new energy scenarios.
(4)
The overall safety requirements of the source–network–load interaction, especially voltage stability, are the basis for ensuring that the source–network–load interaction can be implemented. Good distribution network planning and internal family division of the distribution network can maximize the internal reliability indicators of the distribution network. Therefore, in the design and planning of the new ADNs, a good family analysis is required. Some current grid plans can already meet this demand to a great extent.
However, this article does not consider distribution network scheduling optimization strategies in the context of AC–DC hybrid distribution networks, electric vehicles, or large-scale energy storage applications. From the overall optimization principle, the scheduling unit of the distribution network family can be further refined. It has some similarities with current load agents, microgrids, etc., but its operating mechanism and scheduling mechanism are better. Taking the family as the basic scheduling unit is more conducive to optimization analysis and calculation of the plug-and-play distribution network. There is no need to consider the characteristics of the access node too much, as it will increase the complexity of the model. Therefore, the model of the internal family division of distribution networks should be an important research direction for the optimization and analysis of the scheduling of active distribution networks.

Author Contributions

Methodology, J.D.; Software, Y.Z.; Resources, Y.Z.; Data curation, P.J.; Writing—original draft, P.J.; Funding acquisition, J.D. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the State Grid Economic and Technological Research Institute (PN:SGTYHT/21-JS-223).

Data Availability Statement

Not applicable.

Acknowledgments

The completion of this paper has been helped by many teachers and classmates. We would like to express our gratitude to them for their help and guidance.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jiayi, H.; Chuanwen, J.; Rong, X. A review on distributed energy resources and MicroGrid. Renew. Sustain. Energy Rev. 2008, 12, 2472–2483. [Google Scholar] [CrossRef]
  2. Ehsan, A.; Yang, Q. State-of-the-art techniques for modelling of uncertainties in active distribution network planning: A review. Appl. Energy 2019, 239, 1509–1523. [Google Scholar] [CrossRef]
  3. ALKaabi, S.S.; Zeineldin, H.H.; Member, S.; Khadkikar, V. Planning active distributionnetworks considering multi-DG configurations. IEEE Trans. Power Syst. 2014, 29, 785–793. [Google Scholar] [CrossRef]
  4. Xing, H.; Cheng, H.; Zhang, Y.; Zeng, P. Active distribution network expansion planningintegrating dispersed energy storage systems. IET Gener. Transm Distrib. 2016, 10, 638–644. [Google Scholar] [CrossRef]
  5. Mashayekh, S.; Stadler, M.; Cardoso, G.; Heleno, M.; Madathil, S.C.; Nagarajan, H.; Bent, R.; Mueller-Stoffels, M.; Lu, X.; Wang, J. Security-constrained design of isolated multi-energy microgrids. IEEE Trans. Power Syst. 2018, 33, 2452–2462. [Google Scholar] [CrossRef]
  6. Koltsaklis, N.E.; Kopanos, G.M.; Georgiadis, M.C. Design and operational planning of energy networks based on combined heat and power units. Ind. Eng. Chem. Res. 2014, 53, 16905–16923. [Google Scholar] [CrossRef]
  7. Perković, L.; Mikulčić, H.; Pavlinek, L.; Wang, X.; Vujanović, M.; Tan, H.; Baleta, J.; Duić, N. Coupling of cleaner production with a day-ahead electricity market: A hypothetical case study. J. Clean. Prod. 2017, 143, 1011–1020. [Google Scholar] [CrossRef]
  8. Kumar, S.; Sushama, M. Sushama, Strategic demand response framework for energy management in distribution system based on network loss sensitivity. Energy Environ. 2020, 31, 1385–1402. [Google Scholar] [CrossRef]
  9. Tan, W.-S.; Hassan, M.Y.; Majid, M.S.; Abdul Rahman, H. Optimal distributed renewable generation planning: A review of diferent approaches. Renew. Sustain. Energy Rev. 2013, 18, 626–645. [Google Scholar] [CrossRef]
  10. Santos, S.F.; Fitiwi, D.Z.; Bizuayehu, A.W.; Shafie-khah, M.; Asensio, M.; Contreras, J.; Cabrita, C.M.; Catalão, J.P. Novel multi-stage stochastic DG investment planning with recourse. IEEE Trans. Sustain. Energy 2016, 8, 164–178. [Google Scholar] [CrossRef]
  11. Yang, Q.; Barria, J.A.; Green, T.C. Communication infrastructures for distributed control of power distribution networks. IEEE Trans. Ind. Inform. 2011, 7, 316–327. [Google Scholar] [CrossRef]
  12. Ehsan, A.; Yang, Q.; Cheng, M. A scenario-based robust investment planning model for multi-type distributed generation under uncertainties. IET Gener. Transm. Distrib. 2018, 12, 4426–4434. [Google Scholar] [CrossRef]
  13. Billinton, R.; Allan, R.N. Reliability Evaluation of Power Systems; Springer Science & Business Media: Boston, MA, USA, 2013. [Google Scholar]
  14. Heidari, A.; Agelidis, V.G.; Kia, M.; Pou, J.; Aghaei, J.; Shafie-Khah, M.; Catalão, J.P. Reliability optimization of automated distribution networks with probability customer interruption cost model in the presence of dg units’. IEEE Trans. Smart Grid. 2017, 8, 305–315. [Google Scholar]
  15. Niknam, T.; Kavousifard, A.; Aghaei, J. Scenario-based multiobjective distribution feeder reconfiguration considering wind power using adaptive modified particle swarm optimization. IET Renew. Power Gener. 2012, 6, 236–247. [Google Scholar] [CrossRef]
  16. Ghadi, M.J.; Ghavidel, S.; Rajabi, A.; Azizivahed, A.; Li, L.; Zhang, J. A review on economic and technical operation of active distribution systems. Renew. Sustain. Energy Rev. 2019, 104, 38–53. [Google Scholar] [CrossRef]
  17. Li, Y.; Yang, Z.; Li, G.; Mu, Y.; Zhao, D.; Chen, C.; Shen, B. Optimal scheduling of isolated microgrid with an electric vehicle battery swapping station in multi-stakeholder scenarios: A bi-level programming approach via real-time pricing. Appl. Energy 2018, 232, 54–68. [Google Scholar] [CrossRef]
  18. Mahmood, S.A.; Zedan, M.J.M. Network Load Balancing in Teleconferencing Systems. In Proceedings of the 2022 8th International Engineering Conference on Sustainable Technology and Development (IEC), Erbil, Iraq, 23–24 February 2023; pp. 12–16. [Google Scholar]
  19. 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]
  20. Weishan, D.; Yongjuan, L. Research on grid-connected distributed power grids based on solar power generation systems. Power Technol. 2017, 1052–1054. [Google Scholar]
  21. Li, Y.; Wei, X.; Li, Y.; Dong, Z.; Shahidehpour, M. Detection of false data injection attacks in smart grid: A secure federated deep learning approach. IEEE Trans. Smart Grid 2022, 13, 4862–4872. [Google Scholar] [CrossRef]
  22. Siwu, L.; Sichang, X.; Fangmei, B.; Yulin, Y.; Xiaodi, M.; Chan, P. The analysis of business scenarios and implementation path of “5G+ Source-network-load-storage” multi-station integration. E3S Web Conf. 2021, 248, 02031. [Google Scholar]
  23. Li, Y.; Feng, B.; Wang, B.; Sun, S. Joint planning of distributed generations and energy storage in active distribution networks: A Bi-Level programming approach. Energy 2022, 245, 123226. [Google Scholar] [CrossRef]
  24. Fan, H.; Yu, Z.; Xia, S.; Li, X. Review on coordinated planning of source-network-load-storage for integrated energy systems. Front. Energy Res. 2021, 9, 641158. [Google Scholar] [CrossRef]
  25. Chen, X.U.E.; Jing, R.E.N.; Zhang, X.; Wang, P.; Zhou, X.; Liu, Y.; Kuang, H. An optimal dispatch method for high proportion new energy power grid based on source-network-load-storage interaction. In Proceedings of the 2021 4th International Conference on Electron Device and Mechanical Engineering (ICEDME), Guangzhou, China, 19–21 March 2021; pp. 119–122. [Google Scholar]
  26. Li, Y.; Zhang, M.; Chen, C. A deep-learning intelligent system incorporating data augmentation for short-term voltage stability assessment of power systems. Appl. Energy 2022, 308, 118347. [Google Scholar] [CrossRef]
  27. Xu, R.; Chen, C.; Chen, D.; Zhou, G.; Jin, Y.; Wu, X.; Lin, Z. Maximum openable capacity optimization method of active distribution network considering multiple users access. Energy Rep. 2022, 8, 43–50. [Google Scholar] [CrossRef]
  28. Piao, Z.; Peng, T.; Lu, S.; Zhang, Y.; Cheng, D.; Wang, J. Efficient collaborative utilization mode of “Source-network-load-storage” to improve energy-self-balance ability of microgrid. In Proceedings of the 2022 4th International Conference on Power and Energy Technology (ICPET), Beijing, China, 28–31 July 2022; pp. 498–503. [Google Scholar]
  29. Li, Y.; Yang, Z.; Li, G.; Zhao, D.; Tian, W. Optimal scheduling of an isolated microgrid with battery storage considering load and renewable generation uncertainties. IEEE Trans. Ind. Electron. 2019, 66, 1565–1575. [Google Scholar] [CrossRef]
  30. Wang, J.; Huo, S.; Yan, R.; Cui, Z. Leveraging heat accumulation of district heating network to improve performances of integrated energy system under source-load uncertainties. Energy 2022, 252, 124002. [Google Scholar] [CrossRef]
Figure 1. Basic solution process of ACO.
Figure 1. Basic solution process of ACO.
Wevj 14 00159 g001
Figure 2. 36-node active distribution network structure diagram.
Figure 2. 36-node active distribution network structure diagram.
Wevj 14 00159 g002
Figure 3. The 36-node grid characteristic analyses.
Figure 3. The 36-node grid characteristic analyses.
Wevj 14 00159 g003aWevj 14 00159 g003b
Figure 4. Voltage fluctuation situations of a special user.
Figure 4. Voltage fluctuation situations of a special user.
Wevj 14 00159 g004
Figure 5. The whole net power curve of a typical day.
Figure 5. The whole net power curve of a typical day.
Wevj 14 00159 g005
Figure 6. The 96 data points for a day’s optimization. Each color in the figure represents the predicted terminal load of different branches.
Figure 6. The 96 data points for a day’s optimization. Each color in the figure represents the predicted terminal load of different branches.
Wevj 14 00159 g006
Figure 7. Targets of voltage control for the grid. Each color in the figure represents the voltage of a different node.
Figure 7. Targets of voltage control for the grid. Each color in the figure represents the voltage of a different node.
Wevj 14 00159 g007
Figure 8. Generator combinations of the network. Each color in the figure represents the output power of different units.
Figure 8. Generator combinations of the network. Each color in the figure represents the output power of different units.
Wevj 14 00159 g008aWevj 14 00159 g008bWevj 14 00159 g008c
Figure 9. Load response combinations of the ideal state.
Figure 9. Load response combinations of the ideal state.
Wevj 14 00159 g009
Figure 10. Load response combinations of safe strategies.
Figure 10. Load response combinations of safe strategies.
Wevj 14 00159 g010aWevj 14 00159 g010b
Figure 11. Original combination safety check results. Each color in the figure represents different branches.
Figure 11. Original combination safety check results. Each color in the figure represents different branches.
Wevj 14 00159 g011
Figure 12. Optimization of union calculation results.
Figure 12. Optimization of union calculation results.
Wevj 14 00159 g012
Figure 13. Comparison of the starting state of each node before and after optimization.
Figure 13. Comparison of the starting state of each node before and after optimization.
Wevj 14 00159 g013aWevj 14 00159 g013bWevj 14 00159 g013cWevj 14 00159 g013d
Table 1. Classification of nodes in this case.
Table 1. Classification of nodes in this case.
Node NumberNode TypeResponse CharacteristicsMain PropertiesElastic Index
1LoadRigidCommercial electricity consumption15.21%
2LoadRigidResidential areas4.36%
3Power generationElasticityCascading hydropower40.93%
4Power generationRigidBasin runoff is a small hydropower group12.67%
5LoadRigidChemical plant20.05%
6LoadRigidWater works8.32%
7LoadElasticitySchool, energy storage36.69%
8Power generationElasticityWind power, energy storage28.96%
9LoadRigidMechanical plant9.13%
10LoadRigidFood processing plant10.14%
11LoadElasticityCommercial and self-built photovoltaic energy storage46.68%
12LoadRigidGovernment office31.24%
13Power generationRigidWind force16.21%
14LoadElasticityTextile factory, belt energy storage16.46%
15LoadRigidPaper mill8.83%
16Power generationRigidPhotovoltaic power5.49%
17LoadElasticityCotton textile mill30.13%
18Power generationRigidBasin runoff is a small hydropower group10.25%
19Power generationElasticityBiomass power generation60.14%
20LoadRigidSmall textile mills14.32%
21LoadRigidAutomobile manufacturers10.93%
22LoadRigidEnergy storage battery manufacturer4.16%
23LoadRigidSteel processing plant2.86%
24LoadRigidCommercial office building20.57%
25LoadElasticityAuto parts manufacturer, PV29.12%
26Power generationElasticityGas-fired power generation63.17%
27Power generationElasticityComprehensive energy demonstration park79.23%
Table 2. Classification of families in this case.
Table 2. Classification of families in this case.
Family MembersList of Load NodesList of Power Generation Nodes
11, 2, 7, 9, 11, 174, 8
25, 1519
36, 10, 1213, 32
414, 21, 23, 33, 3536
520, 22, 2518, 26, 29
624, 30, 313
727, 28, 3416
Table 3. Flexible node parameters.
Table 3. Flexible node parameters.
Node NumberNode TypeVoltage LevelMain Parameters (Installed Capacity, in MW)
3Cascading hydropower110 kV2 × 3.5 + 3 × 1.25 + 3 × 0.63 + 2 × 0.32 + 4 × 0.5
8Wind farm35 kV16 × 2 + 1 × 1.5 + 5, Storage: 16
19Biomass power generation35 kV8 × 1.5 + 4 × 3 + 3 × 5
26Gas-fired power generation10 kV4 × 4.5 + 6 × 1.5
32Photovoltaic power station10 kV1.2 × 5, Storage: 3.6
36Storage capacity and power station110 kV6 × 4.6 + 4 × 2.5 + 8 × 1.25
7Class II load10 kVStorage: 1.5. Installation capacity: 3
11Class II load35 kVLight: 1.2 × 3, Storage: 2, Installed capacity: 6.3
14Class I load110 kVStorage: 3.5. Installation capacity: 4.6
25Class I load110 kVLight: 1 × 4, installed capacity: 7.83
27Class II load35 kVLight: 1 × 14, Storage: 12, Gas: 8 × 1.5, Heat: 0.4 × 6
Installed capacity: 34
28Class II load35 kVLight 1.2 × 6, Storage: 10. Pile: 10 × 0.18 + 12 × 0.035
Installed capacity: 16
30Three types of load10 kVStorage: 2. Installation capacity: 2.4
34Class II load35 kVStorage: 6. Installation capacity: 6.3
35Class II load10 kVStorage: 3.2, Pile: 6 × 0.36 + 14 × 0.24, installed capacity: 4.2
Table 4. Original values of generators.
Table 4. Original values of generators.
Generation NodeActive Power (MW)Capacity (MW)
313.201915.28
416.674518.48
817.902538.5
136.887630
16−0.06620
1810.520612.48
1919.96839
269.93627
29−0.147816
320.04436
3640.317247.6
Table 5. Targets of node voltage control.
Table 5. Targets of node voltage control.
Node NameVoltagePhase AngleNode NameVoltagePhase Angle
BUS110BUS191.011764051−31.1464
BUS100.989535501−19.7508BUS201.011727813−31.1951
BUS111.019591909−23.7891BUS211.001133166−38.2493
BUS21.011749759−31.1677BUS260.991429899−41.5284
BUS221.001028601−38.277BUS271−27.8
BUS230.988518931−30.1358BUS281.019514092−16.4535
BUS240.988722183−30.0071BUS291.01950393−16.4518
BUS250.991470488−41.5135BUS341.019493038−16.4671
BUS30.967180274−41.2169BUS40.991399466−41.513
BUS510.994779656−41.5325BUS51.009353618−22.8737
BUS90.994544278−28.8754BUS501.033333767−31.8241
BUS120.990260647−11.9921BUS520.994779656−41.5325
BUS131.002626092−37.9659BUS60.989535501−19.7508
BUS140.998759999−31.5197BUS300.988722183−30.0071
BUS151.039171437−30.079BUS311−41.2368
BUS160.977696715−20.1593BUS331−10.7673
BUS171.046380497−5.50664BUS71−10.6603
BUS181.019633416−23.8103BUS80.989868176−20.2014
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

Jiang, P.; Dong, J.; Zhu, Y. Mixed Linear Model of a Safety Dispatch Model in an Active Distribution Network for Source–Grid–Load Interactions. World Electr. Veh. J. 2023, 14, 159. https://doi.org/10.3390/wevj14060159

AMA Style

Jiang P, Dong J, Zhu Y. Mixed Linear Model of a Safety Dispatch Model in an Active Distribution Network for Source–Grid–Load Interactions. World Electric Vehicle Journal. 2023; 14(6):159. https://doi.org/10.3390/wevj14060159

Chicago/Turabian Style

Jiang, Peng, Jun Dong, and Yuan Zhu. 2023. "Mixed Linear Model of a Safety Dispatch Model in an Active Distribution Network for Source–Grid–Load Interactions" World Electric Vehicle Journal 14, no. 6: 159. https://doi.org/10.3390/wevj14060159

APA Style

Jiang, P., Dong, J., & Zhu, Y. (2023). Mixed Linear Model of a Safety Dispatch Model in an Active Distribution Network for Source–Grid–Load Interactions. World Electric Vehicle Journal, 14(6), 159. https://doi.org/10.3390/wevj14060159

Article Metrics

Back to TopTop