Hierarchical Scheduling Scheme for AC / DC Hybrid Active Distribution Network Based on Multi-Stakeholders

This paper presents a hierarchical multi-stage scheduling scheme for the AC/DC hybrid active distribution network (ADN). The load regulation center (LRC) is considered in the developed scheduling strategy, as well as the AC and DC sub-network operators. They are taken to be different stakeholders. To coordinate the interests of all stakeholders, a two-level optimization model is established. The flexible loads are dispatched by LRC in the upper-level optimization model, the objective of which is minimizing the loss of the entire distribution network. The lower-level optimization is divided into two sub-optimal models, and they are carried out to minimize the operating costs of the AC/DC sub-network operators respectively. This two-level model avoids the difficulty of solving multi-objective optimization and can clarify the role of various stakeholders in the system scheduling. To solve the model effectively, a discrete wind-driven optimization (DWDO) algorithm is proposed. Then, considering the combination of the proposed DWDO algorithm and the YALMIP toolbox, a hierarchical optimization algorithm (HOA) is developed. The HOA can obtain the overall optimization result of the system through the iterative optimization of the upper and lower levels. Finally, the simulation results verify the effectiveness of the proposed scheduling scheme.


Introduction
The active distribution network (ADN) is considered an effective way to address the issues caused by the large-scale integration of distributed energy resources (DERs).For this reason, it has aroused great attention in both research and industry fields, and a great number of studies have been done [1][2][3][4].With the large-scale integration of distributed generation (DG), flexible loads, energy storage and new types of power electronic control devices, a large amount of DC equipment has been installed in AC distribution networks.Due to the large number of electrical energy conversion links, using AC distribution networks for power supply only will influence both the cost and the efficiency [5].The best topology that can accommodate both AC and DC technologies with less need for such conversion is a hybrid one [6].The idea is to merge the AC and DC distribution networks through a bidirectional converter and establish an AC/DC hybrid ADN in which AC-or DC-type energy sources and loads can be flexibly integrated into the distribution networks such that power can smoothly flow between them [7].Furthermore, flexible DC technology is applied to construct the AC/DC hybrid ADN, which can achieve energy scheduling over a wide area.The AC/DC hybrid ADN is controllable and allows DC loads and various DERs to plug-and-play easily; it will be an important implementation form of the future distribution network [8].The characteristics of the AC/DC hybrid ADN provide an efficient means for the integration of DERs with minimum modifications of the current distribution grid [9].At present, most research on AC/DC hybrid ADNs is focused on the optimal power flow (OPF) problem and coordinated control methods.An efficient OPF algorithm for hybrid AC/DC grids with discrete control devices is presented in [10].On this basis, the corrective security constraint is taken into account in OPF for a meshed AC/DC power transmission network [11].Meanwhile, second-order cone programming (SOCP) is applied, and [12] presents a SOCP formulation of the OPF problem for AC/DC systems with voltage source converter (VSC) technology.In terms of coordinated control methods, Ref. [13] presents an overview of various control schemes used for voltage and frequency regulation in standalone and transition mode operation of the hybrid micro-grid.Ref. [14] proposes an instantaneous power-based current control scheme for reactive power compensation in hybrid AC/DC networks, while Ref. [15] presents a novel three-phase reactive power and voltage distributed control method for AC/DC hybrid ADNs based on model predictive control.Furthermore, in order to ensure the safety of the converters and the grid facility, a flexible control strategy is proposed for the AC/DC hybrid grid in [16].These studies provide strong support for the optimal scheduling and operation of the AC/DC hybrid ADN.
In general, the core technologies for operation of AC/DC hybrid ADN contain flexible DC technology and scheduling technology [17].Flexible DC technology can provide a better interface for DGs, energy storage, electric vehicles, and other devices [18][19][20].Meanwhile, the ADN scheduling technology can flexibly dispatch various types of flexible resources in the distribution network, and ultimately realize the optimal operation of the entire system.Using a hybrid-type distribution system can eliminate several conversion stages and thus improve efficiency and reduce investment costs [21].
Recently, several studies have involved the optimization of the scheduling strategies of the AC/DC hybrid ADN.Ref. [22] builds a system structure for an AC/DC hybrid ADN and then establishes a multi-time scale optimal dispatch model based on a multi-agent system (MAS).The model is designed to reduce the operation cost and energy loss on the premise of absorbing the most DERs.As the AC and DC sub-networks perform power interaction through an interlink convertor (ILC), a hierarchical planning method is introduced into the scheduling process.Ref. [23] focuses on power management in the AC/DC microgrid, and the combination of PV, fuel cells, wind, and battery storage with adjustable parameters is analyzed.An optimization model is investigated to reduce the implementation costs by using a multi-objective particle swarm optimization (MOPSO) algorithm.For an ADN with a single DC section and multiple AC sections, [24] proposes a hierarchical and distributed coordinated multi-source optimal scheduling strategy which divides the dispatch area into a local level and an area level.In [25], a two-stage stochastic centralized dispatch scheme is presented, where the first stage produces day-ahead dispatch decisions for the dispatchable DG units, while the second stage determines appropriate corrective decisions for a set of possible scenarios.On this basis, a robust optimal method is introduced.Reference [26] proposes a bi-level two-stage robust optimal scheduling model for AC/DC hybrid multi-microgrids.In this model, the system is divided into the utility level and the supply level, and two-stage robust optimization is carried out.A multi-interval-uncertainty constrained robust dispatch model is proposed in [27], in order to deal with the uncertainties of renewable energy generation and load power in an AC/DC hybrid microgrid.
The above studies made great efforts towards developing optimal operation models for AC/DC hybrid ADN; however, there are still some persisting problems that need to be addressed.Firstly, as different stakeholders, the AC and DC subgrids should have their own scheduling objectives.Meanwhile, the distribution network operator (DNO) needs to regulate the power grid globally.Few studies consider their scheduling strategy based on multi-stakeholders of an AC/DC hybrid ADN.Secondly, most of the present studies pay little attention to the coordination of various schedulable resources in an AC/DC hybrid ADN.In addition, the interaction between the AC and DC sides still deserves further study.Thus, this paper constructs a coordinated optimization framework for AC/DC hybrid ADN based on multi-stakeholders, and a hierarchical multi-stage scheduling model is established.The upper level is optimized by the load regulation center (LRC), taking the minimum loss of the entire distribution network as the scheduling objective.Meanwhile, the lower level is optimized by the AC and DC sub-networks.Therefore, lower-level optimization is divided into two suboptimal problems, which have objective functions related to their own stakeholders.The advantages of using the above hierarchical optimization scheduling model include: (a) avoiding the difficulty of solving multi-objective optimization; and (b) ensuring that all stakeholders pay attention to their own benefits, and obtain the optimal operating state of the system through mutual cooperation and coordination.Meanwhile, considering the combination of the proposed DWDO algorithm and the YALMIP toolbox, the model is solved by a hierarchical optimization algorithm.The simulation results show that the proposed strategy can minimize the operating costs of AC/DC sub-networks while reducing system losses.
In summary, the contributions of this paper are as follows: (a) a hierarchical multi-stage scheduling model for AC/DC hybrid and is proposed, to avoid the difficulty of solving multi-objective optimization; (b) AC and DC sub-networks are regarded as different stakeholders so that their interaction characteristics are analyzed in a scheduling problem; and (c) a hierarchical optimization algorithm comprising the DWDO algorithm and YALMIP toolbox is developed to solve the hierarchical optimization model.
The rest of this paper is organized into the following sections.Section 2 describes the network structure of the AC/DC hybrid ADN and proposes a hierarchical multi-stage scheduling framework based on multi-stakeholders.Section 3 details the mathematical formulation of the hierarchical optimization scheduling model.In Section 4, a hierarchical optimization algorithm is proposed and researched.In Section 5, case studies and discussions are carried out.Finally, Section 6 put forwards the conclusions that can be drawn.

System Framework
As mentioned above, various types of distributed scheduling resources are integrated into the AC/DC hybrid ADN.When the distribution network is in operation, the schedulable resources will interact with each other.Scheduling only a single or partial resource makes it difficult to achieve overall system optimization, and easily leads to redundancy and waste of resources.Thus, global optimization of AC/DC hybrid ADN is necessary.On the other hand, each of the AC and DC areas constitutes a sub-network, which can be scheduled as an independent stakeholder.In view of this, the network structure is shown in Figure 1.The arrows in the figure show the positive power direction.Meanwhile, the lower level is optimized by the AC and DC sub-networks.Therefore, lower-level optimization is divided into two suboptimal problems, which have objective functions related to their own stakeholders.The advantages of using the above hierarchical optimization scheduling model include: (a) avoiding the difficulty of solving multi-objective optimization; and (b) ensuring that all stakeholders pay attention to their own benefits, and obtain the optimal operating state of the system through mutual cooperation and coordination.Meanwhile, considering the combination of the proposed DWDO algorithm and the YALMIP toolbox, the model is solved by a hierarchical optimization algorithm.The simulation results show that the proposed strategy can minimize the operating costs of AC/DC sub-networks while reducing system losses.In summary, the contributions of this paper are as follows: (a) a hierarchical multi-stage scheduling model for AC/DC hybrid and is proposed, to avoid the difficulty of solving multi-objective optimization; (b) AC and DC sub-networks are regarded as different stakeholders so that their interaction characteristics are analyzed in a scheduling problem; and (c) a hierarchical optimization algorithm comprising the DWDO algorithm and YALMIP toolbox is developed to solve the hierarchical optimization model.
The rest of this paper is organized into the following sections.Section 2 describes the network structure of the AC/DC hybrid ADN and proposes a hierarchical multi-stage scheduling framework based on multi-stakeholders.Section 3 details the mathematical formulation of the hierarchical optimization scheduling model.In Section 4, a hierarchical optimization algorithm is proposed and researched.In Section 5, case studies and discussions are carried out.Finally, Section 6 put forwards the conclusions that can be drawn.

System Framework
As mentioned above, various types of distributed scheduling resources are integrated into the AC/DC hybrid ADN.When the distribution network is in operation, the schedulable resources will interact with each other.Scheduling only a single or partial resource makes it difficult to achieve overall system optimization, and easily leads to redundancy and waste of resources.Thus, global optimization of AC/DC hybrid ADN is necessary.On the other hand, each of the AC and DC areas constitutes a sub-network, which can be scheduled as an independent stakeholder.In view of this, the network structure is shown in Figure 1.The arrows in the figure show the positive power direction.In Figure 1, the AC/DC hybrid ADN is divided into an AC section and a DC section.The main components of the DC section include photovoltaics (PVs), wind turbines (WT), fuel cells (FC), compressed air energy storage (CAES), and DC loads, which are connected to the DC bus.The AC section mainly constitutes a micro-gas turbine (MT) and AC loads, with the upper-level power grid In Figure 1, the AC/DC hybrid ADN is divided into an AC section and a DC section.The main components of the DC section include photovoltaics (PVs), wind turbines (WT), fuel cells (FC), compressed air energy storage (CAES), and DC loads, which are connected to the DC bus.The AC section mainly constitutes a micro-gas turbine (MT) and AC loads, with the upper-level power grid (UPG) connected to it.The upper-level power grid supplies electricity based on time-of-use price.Meanwhile, the DC section interacts with the AC section through ILC, which can transmit power in both directions.Additionally, in order to reflect the operating characteristics of the ADN, flexible loads are considered to participate in the scheduling.
In general, the power on both the AC and DC sides influences each other, and it is difficult to directly obtain the optimal operating status of the entire system.Thus, a hierarchical multi-stage scheduling framework based on multi-stakeholders is proposed, as shown in Figure 2. It is necessary to consider both active and reactive power in the optimal power flow of the distribution network [28].However, in this paper, the main purpose of the lower-level optimization is to determine the active economic dispatch of the power system.Meanwhile, considering the complexity of this model, only active power is analyzed in both the AC and DC sections.
Energies 2018, 11, 2830 4 of 17 (UPG) connected to it.The upper-level power grid supplies electricity based on time-of-use price.
Meanwhile, the DC section interacts with the AC section through ILC, which can transmit power in both directions.Additionally, in order to reflect the operating characteristics of the ADN, flexible loads are considered to participate in the scheduling.
In general, the power on both the AC and DC sides influences each other, and it is difficult to directly obtain the optimal operating status of the entire system.Thus, a hierarchical multi-stage scheduling framework based on multi-stakeholders is proposed, as shown in Figure 2. It is necessary to consider both active and reactive power in the optimal power flow of the distribution network [28].However, in this paper, the main purpose of the lower-level optimization is to determine the active economic dispatch of the power system.Meanwhile, considering the complexity of this model, only active power is analyzed in both the AC and DC sections.In Figure 2, the stakeholder of the upper-level optimization is the load regulation center (LRC).All flexible loads of DC and AC sections are managed by the LRC, in order to reduce the losses of the distribution network.Then, the scheduling results of the flexible loads are transferred to the lower-level optimization, which is comprised of the internal optimization of both the DC network and the AC network.The DGs and CAES are scheduled by the lower-level optimization, as well as the transmission power of the UPG.The results are also transferred to the upper level, and the system performs repeated iterative optimization until it reaches a dynamic balance.Through the hierarchical optimization method, optimal scheduling results for the entire AC/DC hybrid ADN are expected to finally be obtained.The results not only reflect the global optimization of the power grid, but also reflect the local optimization of the AC and DC power grids.
In the next section, we will establish the hierarchical optimization model, and give a detailed solution process for it.

Optimization Model of the Upper Level
In the upper-level optimization, flexible loads of the entire network are scheduled.The scheduling objective is to reduce system losses, which mainly include line loss and loss of ILC.Therefore, the objective function is expressed in Equation (1).
where T is the scheduling period.
The detailed expression of the above power loss can be expressed as Equations ( 2) and (3).In Figure 2, the stakeholder of the upper-level optimization is the load regulation center (LRC).All flexible loads of DC and AC sections are managed by the LRC, in order to reduce the losses of the distribution network.Then, the scheduling results of the flexible loads are transferred to the lower-level optimization, which is comprised of the internal optimization of both the DC network and the AC network.The DGs and CAES are scheduled by the lower-level optimization, as well as the transmission power of the UPG.The results are also transferred to the upper level, and the system performs repeated iterative optimization until it reaches a dynamic balance.Through the hierarchical optimization method, optimal scheduling results for the entire AC/DC hybrid ADN are expected to finally be obtained.The results not only reflect the global optimization of the power grid, but also reflect the local optimization of the AC and DC power grids.
In the next section, we will establish the hierarchical optimization model, and give a detailed solution process for it.

Optimization Model of the Upper Level
In the upper-level optimization, flexible loads of the entire network are scheduled.The scheduling objective is to reduce system losses, which mainly include line loss and loss of ILC.Therefore, the objective function is expressed in Equation (1).
where T is the scheduling period.The detailed expression of the above power loss can be expressed as Equations ( 2) and (3). (2) where l 1 and l 2 are lines of the AC and DC networks, respectively, and N AC line and N DC line are the total numbers of lines in AC and DC networks, respectively.
The active and reactive power of the network lines are obtained using the existing power flow calculation method.In the optimization process, the node voltage amplitude constraint of the network is considered to be as in Equation (4).
where N node is total network node of the system.ε 1,s,t and ε 2,s,t are relaxation variables.In practical calculation, voltage constraints are added as penalty terms to the objective function to avoid the occurrence of voltage over-limit conditions.In this model, the flexible loads are mainly transferable loads (TLs), which can be expressed as Equation (5).
where α j,t ∈ {0, 1} are the actual operating state of transferable load j at time t, 1 is operation and 0 is outage.N TL is the total number of transferable loads in the system.Equation ( 5) is the constraint condition for the flexible load scheduling.Through the upper-level optimization, the flexible load scheduling result is obtained.Then it is passed to the lower level.Thus, the actual loads of DC and AC sections can be expressed as Equation (6).

Optimization Model of the Lower Level
The lower-level optimization is divided into two parts: the DC section and the AC section, which have their own optimization objectives.As different stakeholders, they are committed to minimizing operating costs in their respective section.

Optimization Model of the DC Section
For the DC section, the scheduling costs mainly include the maintenance cost of generators and energy storage, the fuel cost of FC and CAES, and the electricity transaction cost.The objective function is shown in Equation (7).
The detailed expressions of the costs are expressed in Equation (8).
where N PV , N WT , N CAES , N FC are the total number of PV, WT, CAES and FC in the system.The second items in the right of C PV,t and C WT,t are the penalty cost of PV and wind power curtailment.These can ensure that the power of PV and WT is curtailed as little as possible.
Since CAES that has been applied on a large scale at this stage is mostly supplementarily fired, it needs to consume natural gas to supplement combustion during discharge stage.This cost has been considered in Equation (8).The simplified power model of CAES can be described as Equation ( 9). .
where k CAES_C and k CAES_G are the operation coefficients for compression and generation conditions, respectively, and their detailed expressions are presented in literature [29].
As CAES usually uses the working pressure of the gas storage chamber to indicate its charge and discharge capacity, the equivalent state of charge (SOC) can be described as Equation (10).
u CAES_C,t + u CAES_G,t ≤ 1 (13) Equation ( 11) is the power constraint of CAES for the charge and discharge stages.Equation ( 12) is the equivalent SOC constraint of CAES.Equation ( 13) is the charge and discharge status constraint of CAES; this constraint ensures that CAES can only be charging, discharging or in a stop state at any given time.Equation ( 14) is the equivalent SOC balance constraint of CAES; this constraint ensures that CAES maintains the pressure of the gas storage chamber before and after a scheduling cycle.
where Equation ( 15) is the power balance constraint of the DC section.Equation ( 16) is the output power constraint of the FC.Equation ( 17) is the output power constraint of the ILC.

Optimization Model of the AC Section
Similarly, the objective function of the AC section is shown in Equation (18). min The constraints are shown in Equations ( 19)- (22).
Equation ( 19) is the power balance constraint of the AC section, Equations ( 20) and ( 21) are the power constraints of UPG and MT respectively.According to the positive power direction shown in Figure 1, Equation ( 22) ensures the power balance of both the DC and AC sections.
Thus, the hierarchical scheduling model is established.It is summarized in Equation (23).

Model Solution
The optimization model of the AC/DC hybrid ADN is a bi-level optimization problem.The'network loss information in the upper-level optimization needs to be obtained according to the lower-level optimization result.Conversely, the lower-level optimization requires the output information of flexible loads determined in the upper-level optimization.Based on existing research, this chapter will present a DWDO algorithm, and combine it with the YALMIP toolkit to solve the bi-level optimization problem.

Discrete Wind-Driven Optimization Algorithm
Through the treatment of flexible loads in this paper, the upper-level optimization is a nonlinear integer optimization problem.At present, most nonlinear optimization problems are solved by intelligent algorithms, such as particle swarm optimization (PSO), genetic algorithm (GA), and so on.However, most of them have the problems that they find it easy to fall into local optima and are computationally inefficient.Thus, a wind-driven optimization (WDO) technique was proposed by [30], which is also an iterative heuristic global optimization algorithm.WDO is used to describe the motion of air parcels within the earth's atmosphere, and Newton's second law of motion is used to describe the N-dimensional search space.As a conclusion, [31] points out that WDO is well-suited for problems with both discrete and continuous-valued parameters.The algorithm process of WDO is as follows.
(1) Initialize the population size, set the maximum number of iterations, related parameters, and search boundaries, and define pressure functions.(2) Initialize the air mass points and assign the starting speeds and locations randomly.(3) Calculate the pressure value of the air particles in the current iteration and rearrange the population according to the pressure values.(4) Update the speeds of the air particles using Equation (24).
where u new is the updated speed vector of air mass points, u cur is current speed vector, x cur is current location vector, x opt is the global optimum location so far, and u otherdim cur is the speed vector of any other dimension.i is the arrangement number of the air mass points in step 3. α, g, RT and c are parameters which will be changed according to the optimization problems.
(5) Update the positions of the air particles using Equation (25).
where x new is the updated location vector of the air mass points.
(6) If the termination condition is not reached, go to step 3.
The purpose of the upper-level optimization is to determine the output state of the flexible loads at each time.This is a discrete nonlinear optimization problem, which cannot be solved directly using the WDO algorithm.Therefore, the original WDO algorithm needs to be discretized, which is called DWDO in this paper.In DWDO, the starting position parameters should be randomly initialized as 0-1 integer variables.The speed update formula is the same as Equation ( 24), but the positions need to be updated based on the probability of taking 1.Now, introduce the sigmoid function, which is shown in Equation (26).
In this place, s(u new ) indicates the probability of taking 1 for x new .So the locations of the air particles can be updated by Equation (27).
The remaining optimization steps of DWDO are similar to the original WDO algorithm.Using the DWDO algorithm, the upper-level optimization can obtain the optimized output for flexible loads.

Solution Process of the Hierarchical Scheduling Model
The lower-level optimization is a mixed-integer linear programming (MILP) problem.As a mature optimization toolkit, YALMIP has great advantages in solving MILP problems, and has been widely used [32][33][34].Furthermore, it can easily interact with MATLAB.Therefore, YALMIP is used to solve lower-level optimization and iterates with the DWDO algorithm in the upper-level optimization.The solution process is shown in Figure 3.
As this paper focuses on active economic dispatch of ADN, only active power is allowed to be exchanged between the ADN and the UPG.The DC and AC sections exchange active power through ILC.In addition to this, the AC section also exchanges active power with the UPG through the boundary line.Reactive power is generated locally by suitable means to satisfy the reactive demand and losses for both DC and AC sections.Thus, the solution process of the hierarchical scheduling model is as follows.
(1) Initialize the values of P CAES , P FC , P MT , P ILC using operation data of another similar day, and preset the planned output of the flexible loads.(2) Solve the upper-level optimization using DWDO, in order to obtain the optimized output of the flexible loads (P * TL_DC ,P * TL_AC ).In this process, we expect that PV and wind power can be absorbed as much as possible.In addition, as the AC side is connected to the upper-level power grid, its internal power adjustment ability is strong.Therefore, the optimization starts from the DC side in this step.(3) Solve the optimization of the DC and AC sections using YALMIP, in order to obtain P * CAES , P * FC , P * ILC , P * MT .(4) Determine whether the algorithm reaches a defined number of iterations or convergence accuracy.
If the change of the system cost is lower than a certain constant between two iterations, the algorithm converges and then the result is presented.If not, go to step 2.
Energies 2018, 11, 2830 9 of 17 The remaining optimization steps of DWDO are similar to the original WDO algorithm.Using the DWDO algorithm, the upper-level optimization can obtain the optimized output for flexible loads.

Solution Process of the Hierarchical Scheduling Model
The lower-level optimization is a mixed-integer linear programming (MILP) problem.As a mature optimization toolkit, YALMIP has great advantages in solving MILP problems, and has been widely used [32][33][34].Furthermore, it can easily interact with MATLAB.Therefore, YALMIP is used to solve lower-level optimization and iterates with the DWDO algorithm in the upper-level optimization.The solution process is shown in Figure 3.As this paper focuses on active economic dispatch of ADN, only active power is allowed to be exchanged between the ADN and the UPG.The DC and AC sections exchange active power through ILC.In addition to this, the AC section also exchanges active power with the UPG through the boundary line.Reactive power is generated locally by suitable means to satisfy the reactive demand and losses for both DC and AC sections.Thus, the solution process of the hierarchical scheduling model is as follows.).In this process, we expect that PV and wind power can be absorbed as much as possible.In addition, as the AC side is connected to the upper-level power grid, its internal power adjustment ability is strong.Therefore, the optimization starts from the DC side in this step.(4) Determine whether the algorithm reaches a defined number of iterations or convergence accuracy.If the change of the system cost is lower than a certain constant between two iterations, the algorithm converges and then the result is presented.If not, go to step 2.

Case Study and Discussion
The hierarchical scheduling model is tested in a modified 38-node test system, which is shown in Figure 4.The modified test system is divided into a DC section and an AC section.The DC section includes two PVs, one WT, one FC and one CAES, while the AC section includes one MT

Case Study and Discussion
The hierarchical scheduling model is tested in a modified 38-node test system, which is shown in Figure 4.The modified test system is divided into a DC section and an AC section.The DC section includes two PVs, one WT, one FC and one CAES, while the AC section includes one MT and is connected to the upper-level power grid.The configuration of loads is also shown in Figure 4.The line parameters are given in literature [35].In this system, the forecast output of the PVs and WT is shown in Figure 5, while the total rigid load prediction curve of the whole system is shown in Figure 6.The data for transferable loads is shown in Table 1.The data for CAES is shown in Table 2.The time-of-use price (TOU) of UPG is shown in Table 3.The maintenance cost of the PVs, WT, FC and MT are 0.0096 CNY/kWh, 0.0296 In this system, the forecast output of the PVs and WT is shown in Figure 5, while the total rigid load prediction curve of the whole system is shown in Figure 6.The data for transferable loads is shown in Table 1.The data for CAES is shown in Table 2.The time-of-use price (TOU) of UPG is shown in Table 3.The maintenance cost of the PVs, WT, FC and MT are 0.0096 CNY/kWh, 0.0296 CNY/kWh, 0.088 CNY/kWh and 0.088 CNY/kWh, respectively.The fuel cost of the FC and MT are 0.525 CNY/kWh and 0.641 CNY/kWh, respectively.The penalty cost coefficients of PV and wind power curtailment is 10 CNY/kWh.The power limits of the ILC and the UPG are ±1000 kW and ±1500 kW, respectively.
The parameters of DWDO are α = 0.4, g = 0.2, RT = 3, c = 0.4.The number of air mass points is 30, the number of iterations is 500.All programs are coded and tested in MATLAB, and the lower-level optimization is solved by YALMIP.The programs run on an Intel(R) Core(TM) i7-4710MQ @ 2.50 GHz, RAM 8 GB system.In this system, the forecast output of the PVs and WT is shown in Figure 5, while the total rigid load prediction curve of the whole system is shown in Figure 6.The data for transferable loads is shown in Table 1.The data for CAES is shown in Table 2.The time-of-use price (TOU) of UPG is shown in Table 3.The maintenance cost of the PVs, WT, FC and MT are 0.0096 CNY/kWh, 0.0296 CNY/kWh, 0.088 CNY/kWh and 0.088 CNY/kWh, respectively.The fuel cost of the FC and MT are 0.525 CNY/kWh and 0.641 CNY/kWh, respectively.The penalty cost coefficients of PV and wind power curtailment is 10 CNY/kWh.The power limits of the ILC and the UPG are 1000kW  and 1500  kW, respectively.The parameters of DWDO are =0.4 . The number of air mass points is 30, the number of iterations is 500.All programs are coded and tested in MATLAB, and the lower-level optimization is solved by YALMIP.The programs run on an Intel(R) Core(TM) i7-4710MQ @ 2.50 GHz, RAM 8 GB system.To analyze the effectiveness of the proposed hierarchical scheduling strategy, we consider the following two case studies: CASE 1: Hierarchical optimization of two levels, but the lower-level is optimized as a whole.CASE 2: Hierarchical optimization of two levels, and the lower-level is divided into two stakeholders.

Economic Comparison
CASE 1: This case is the conventional scheduling strategy that already exists.In this case, the AC/DC hybrid ADN is regarded as one stakeholder, and the ILC power limitation is not considered.The scheduling result is shown in Figure 7.In the depth of the night, the CAES charges, as both the system load and the electricity price of the UPG are low.At noon, from 12:00 to 14:00, the outputs of the PVs are high, and the power supply for the ADN is sufficient.The electricity price of the UPG reaches its maximum, which is higher than the MT and FC.Therefore, all DGs run in order to sell electricity to the UPG in order to gain substantial profits.At 15:00-17:00, although the rigid load of the system reaches peak value, the power of the PVs and WT is able to meet most load requirements.At 18:00-21:00, the electricity price of the UPG reaches maximum again.Now the outputs of the PVs are almost zero, and the FC and MT are running to sell electricity to the UPG.However, at this time, the sales power is less than that at noon.As a result, the operating costs of the AC and DC sections are 16,764 CNY and −2967 CNY respectively, and the cost of system loss is 4651 CNY.

Economic Comparison
CASE 1: This case is the conventional scheduling strategy that already exists.In this case, the AC/DC hybrid ADN is regarded as one stakeholder, and the ILC power limitation is not considered.The scheduling result is shown in Figure 7.In the depth of the night, the CAES charges, as both the system load and the electricity price of the UPG are low.At noon, from 12:00 to 14:00, the outputs of the PVs are high, and the power supply for the ADN is sufficient.The electricity price of the UPG reaches its maximum, which is higher than the MT and FC.Therefore, all DGs run in order to sell electricity to the UPG in order to gain substantial profits.At 15:00-17:00, although the rigid load of the system reaches peak value, the power of the PVs and WT is able to meet most load requirements.At 18:00-21:00, the electricity price of the UPG reaches maximum again.Now the outputs of the PVs are almost zero, and the FC and MT are running to sell electricity to the UPG.However, at this time, the sales power is less than that at noon.As a result, the operating costs of the AC and DC sections are 16,764 CNY and −2967 CNY respectively, and the cost of system loss is 4651 CNY.
CASE 2: This case is the proposed hierarchical scheduling strategy in this paper.The lower level is divided into two stakeholders, and the transmission power limit between the two stakeholders is taken into account.The scheduling result is shown in Figure 8.Compared to CASE 1, the ADN will also sell electricity to the UPG at 12:00-14:00 and 18:00-21:00 in CASE 2, but the sales power is less than CASE 1.As a result, the operating costs of the AC and DC sections are 16,673 CNY and 2580 CNY, respectively, and the cost of the system loss is 4643 CNY.Obviously, the total cost in CASE 2 is more than that in CASE 1, due to the transmission power limitation of ILC.From Figure 7, we know that the maximum transmission power of ILC is about 1800 kW, which is far more than the 1000 kW in CASE 2. This means more investment is needed in the construction of the distribution network.Thus, it is inappropriate to regard the AC/DC hybrid ADN as one stakeholder without considering the constraints between the AC and DC sub-networks.CASE 2: This case is the proposed hierarchical scheduling strategy in this paper.The lower level is divided into two stakeholders, and the transmission power limit between the two stakeholders is taken into account.The scheduling result is shown in Figure 8.Compared to CASE 1, the ADN will also sell electricity to the UPG at 12:00-14:00 and 18:00-21:00 in CASE 2, but the sales power is less than CASE 1.As a result, the operating costs of the AC and DC sections are 16,673 CNY and 2580 CNY, respectively, and the cost of the system loss is 4643 CNY.Obviously, the total cost in CASE 2 is more than that in CASE 1, due to the transmission power limitation of ILC.From Figure 7, we know that the maximum transmission power of ILC is about 1800 kW, which is far more than the 1000 kW in CASE 2. This means more investment is needed in the construction of the distribution network.Thus, it is inappropriate to regard the AC/DC hybrid ADN as one stakeholder without considering the constraints between the AC and DC sub-networks.
Figure 7, we know that the maximum transmission power of ILC is about 1800 kW, which is far more than the 1000 kW in CASE 2. This means more investment is needed in the construction of the distribution network.Thus, it is inappropriate to regard the AC/DC hybrid ADN as one stakeholder without considering the constraints between the AC and DC sub-networks.

Scheduling Result of Loads
The scheduling result of the loads in CASE 2 is shown in Figure 9.The planned loads are low at 0:00-6:00 and high at 15:00-17:00.This may lead to excess or insufficient power supply for the ADN.By means of the optimization scheduling, flexible loads are transferred to achieve peak load shifting.The scheduled load curve is more stable, and the load peak-to-valley difference is smaller.

Scheduling Result of Loads
The scheduling result of the loads in CASE 2 is shown in Figure 9.The planned loads are low at 0:00-6:00 and high at 15:00-17:00.This may lead to excess or insufficient power supply for the ADN.By means of the optimization scheduling, flexible loads are transferred to achieve peak load shifting.The scheduled load curve is more stable, and the load peak-to-valley difference is smaller.

ILC Power Analysis
The power balance of ILC is important in the optimization process.Usually, in optimization scheduling, decision-makers mainly focus on the active power.The power profiles of ILC on the DC side and the AC side are shown in Figure 10.According to the previous analysis, we know that a large amount of load was transferred to 0:00-2:00.Now, the electricity price of UPG is low and the AC sub-network buys electricity from the UPG.Part of the purchased electricity is transmitted to the DC sub-network through the ILC.Between around 10:00 and 12:00, the outputs of WT and PVs are high, and ILC delivers power from the DC sub-network to the AC sub-network at maximum power.At 18:00-21:00, the electricity price of UPG is high, and the output of WT is high.

ILC Power Analysis
The power balance of ILC is important in the optimization process.Usually, in optimization scheduling, decision-makers mainly focus on the active power.The power profiles of ILC on the DC side and the AC side are shown in Figure 10.According to the previous analysis, we know that a large amount of load was transferred to 0:00-2:00.Now, the electricity price of UPG is low and the AC sub-network buys electricity from the UPG.Part of the purchased electricity is transmitted to the DC sub-network through the ILC.Between around 10:00 and 12:00, the outputs of WT and PVs are high, and ILC delivers power from the DC sub-network to the AC sub-network at maximum power.At 18:00-21:00, the electricity price of UPG is high, and the output of WT is high.Thus, power is also delivered from the DC sub-network to the AC sub-network during this time.It can be seen from Figure 10 that the power of the ILC on the AC and DC side is balanced in a scheduling period.

ILC Power Analysis
The power balance of ILC is important in the optimization process.Usually, in optimization scheduling, decision-makers mainly focus on the active power.The power profiles of ILC on the DC side and the AC side are shown in Figure 10.According to the previous analysis, we know that a large amount of load was transferred to 0:00-2:00.Now, the electricity price of UPG is low and the AC sub-network buys electricity from the UPG.Part of the purchased electricity is transmitted to the DC sub-network through the ILC.Between around 10:00 and 12:00, the outputs of WT and PVs are high, and ILC delivers power from the DC sub-network to the AC sub-network at maximum power.At 18:00-21:00, the electricity price of UPG is high, and the output of WT is high.Thus, power is also delivered from the DC sub-network to the AC sub-network during this time.It can be seen from Figure 10 that the power of the ILC on the AC and DC side is balanced in a scheduling period.

Rejected Power Analysis
Due to the existence of system constraints, the active power may be rejected by other stakeholders.The actual output powers of WT and PVs are shown in Figure 11.The actual outputs of PVs are the same as the prediction curve, but there is a reduction in the power of WT.In a scheduling period, the maximum total output of wind power is 19,811 kWh, but the total output is actually 19,501 kWh.The transmission power constraints of ILC and UPG will affect it, as well as the capacity of CAES.

Rejected Power Analysis
Due to the existence of system constraints, the active power may be rejected by other stakeholders.The actual output powers of WT and PVs are shown in Figure 11.The actual outputs of PVs are the same as the prediction curve, but there is a reduction in the power of WT.In a scheduling period, the maximum total output of wind power is 19,811 kWh, but the total output is actually 19,501 kWh.The transmission power constraints of ILC and UPG will affect it, as well as the capacity of CAES.For the DC section, the rejected power of WT and PVs is mainly related to the ILC power limitation and the capacity of the energy storage system.Meanwhile, for the whole AC/DC hybrid ADN, the rejected power is mainly affected by the power limitation of the boundary line with UPG.For instance, when the transmission power limit of the UPG is increased to 3000 kW, no rejected power will appear.Detailed quantitative analysis will be carried out in future research.

Conclusions and Prospects
As an AC/DC hybrid, ADN is becoming an important implementation forms of the future distribution network, its global optimization scheduling is of great significance.For this issue, a hierarchical multi-stage scheduling model is established in this paper.The upper level optimizes the total system losses by scheduling the flexible loads, while the lower level minimizes the scheduling costs of the AC and DC sections by scheduling the controllable DGs and CASE.The For the DC section, the rejected power of WT and PVs is mainly related to the ILC power limitation and the capacity of the energy storage system.Meanwhile, for the whole AC/DC hybrid ADN, the rejected power is mainly affected by the power limitation of the boundary line with UPG.For instance, when the transmission power limit of the UPG is increased to 3000 kW, no rejected power will appear.Detailed quantitative analysis will be carried out in future research.

Conclusions and Prospects
As an AC/DC hybrid, ADN is becoming an important implementation forms of the future distribution network, its global optimization scheduling is of great significance.For this issue, a hierarchical multi-stage scheduling model is established in this paper.The upper level optimizes the total system losses by scheduling the flexible loads, while the lower level minimizes the scheduling costs of the AC and DC sections by scheduling the controllable DGs and CASE.The DWDO algorithm and the YALMIP toolbox are combined to solve the scheduling model.The case study shows that the proposed scheduling scheme can achieve lower system loss with low operating costs of both AC and DC sections.Thus, not only is the overall system optimized, but the profits of multi-stakeholders are also taken into account.The research in this paper can provide several theoretical references for the scheduling of the AC/DC hybrid distribution network.

Figure 1 .
Figure 1.Network structure of the AC/DC hybrid ADN.

Figure 1 .
Figure 1.Network structure of the AC/DC hybrid ADN.

Figure 2 .
Figure 2. Hierarchical multi-stage scheduling framework based on multi-stakeholders.

Figure 2 .
Figure 2. Hierarchical multi-stage scheduling framework based on multi-stakeholders.

Figure 3 .
Figure 3. Solution process of the hierarchical scheduling model.

( 1 )
Initialize the values of CAES P , FC P , MT P , ILC P using operation data of another similar day, and preset the planned output of the flexible loads.(2) Solve the upper-level optimization using DWDO, in order to obtain the optimized output of the flexible loads (

( 3 )
Solve the optimization of the DC and AC sections using YALMIP, in order to obtain * CAES

Figure 3 .
Figure 3. Solution process of the hierarchical scheduling model.

Figure 5 .
Figure 5. Prediction curve of PVs and WT.Figure 5. Prediction curve of PVs and WT.

Figure 10 .
Figure 10.Power of ILC on both the DC and the AC side.Figure 10.Power of ILC on both the DC and the AC side.

Figure 10 .
Figure 10.Power of ILC on both the DC and the AC side.Figure 10.Power of ILC on both the DC and the AC side.

Figure 11 .
Figure 11.Actual output power and MPPT power of WT, PV1 and PV2.

Figure 11 .
Figure 11.Actual output power and MPPT power of WT, PV1 and PV2.
,PV P PV,i PV ,t + c aban,PV (P PV,i PV ,t − P PV,i PV ,t )]∆t C WT,t = ∑ i WT ∈N WT [c o,WT P WT,i WT ,t + c aban,WT (P WT,i WT ,t − P WT,i WT ,t )]∆t C CAES,t = ∑ i CAES ∈N CAES c CAES P CAES,i CAES,t u CAES_G,t + c o,CAES P CAES,i CAES,t ∆t C FC,t = ∑ i FC ∈N FC c FC P FC,i FC,t + c o,FC P FC,i FC,t ∆t C ILC_DC,t = c g,t P ILC_DC,t ∆t

Table 1 .
Data of transferable loads.

Table 1 .
Data of transferable loads.

Table 3 .
Time-of-use price of UPG.