Journal of Modern Power Systems and Clean Energy

ISSN 2196-5625 CN 32-1884/TK

网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

A Projection and Decomposition Approach for Multi-agent Coordinated Scheduling in Power Systems  PDF

  • Haifeng Qiu 1 (Senior Member, IEEE)
  • Zhigang Li 2 (Senior Member, IEEE)
  • Hongjun Gao 3 (Senior Member, IEEE)
  • Hung Dinh Nguyen 1 (Member, IEEE)
  • Veerapandiyan Veerasamy 1 (Member, IEEE)
  • Hoay Beng Gooi 1 (Fellow, IEEE)
1. School of Electrical and Electronic Engineering, Nanyang Technological University, 639798 Singapore; 2. School of Electric Power Engineering, South China University of Technology, Guangzhou 510641, China; 3. College of Electrical Engineering, Sichuan University, Chengdu 610065, China

Updated:2024-05-20

DOI:10.35833/MPCE.2023.000422

  • Full Text
  • Figs & Tabs
  • References
  • Authors
  • About
CITE
OUTLINE

Abstract

Aiming at multi-agent coordinated scheduling problems in power systems under uncertainty, a generic projection and decomposition (P&D) approach is proposed in this letter. The canonical min-max-min two-stage robust optimization (TSRO) model with coupling constraints is equivalent to a concise robust optimization (RO) model in the version of mixed-integer linear programming (MILP) via feasible region projection. The decentralized decoupling of the non-convex MILP problem is realized through a dual decomposition algorithm, which ensures the fast convergence to a high-quality solution in the distributed optimization. Numerical tests verify the superior performance of the proposed P&D approach over the existing distributed TSRO method.

I. Introduction

WITH the growing penetration of renewables, its strong uncertainty brings considerable challenges to power system scheduling. In addition, due to the coexistence of multiple agents in power systems, the efficient coordination of independent stakeholders considering privacy-preserving has become the focus of the current research [

1]. Distributed two-stage robust optimization (TSRO) models have been investigated for multi-agent coordinated scheduling under uncertainty [2].

At present, the traditional distributed TSRO method deals with uncertainty by constructing a centralized min-max-min RO model. The augmented Lagrangian method, e.g., the alternating direction method of multipliers (ADMM) [

3] or the analysis target cascading (ATC) method [4], further decouples the centralized TSRO model in a distributed manner. The above method has been applied to a few coordinated scheduling problems [2]-[5]. Nonetheless, some limitations exist in the application of the traditional distributed TSRO method.

1) TSRO is generally solved by an iterative algorithm, such as the column and constraint generation (C&CG) algorithm or the Benders decomposition. The iterative calculations could hardly satisfy the computational requirements in practical cases due to their large calculation expense. Such a burden grows particularly heavy if the conventional iterative solution algorithm is embedded in distributed decomposition [

6]. The TSRO and the iterative solution algorithm would be therefore inappropriate to execute online distributed calculation for scheduling.

2) The representative augmented Lagrangian-based distributed methods can equivalently decouple convex centralized models, and such methods could derive the identical results as the centralized optimization [

7]. However, the TSRO model features the programming in a min-max-min formulation, and the computational optimality and convergence cannot be warranted for such complex models using conventional distributed methods [8]. Although some techniques have been devised to address the convergence and optimality issues [9], they are not universally efficient, especially when coupling constraints exist in both stages of a TSRO model. As the computational performance of the traditional distributed TSRO model could hardly fulfill the requirements of practical applications, there arises a need to bring a concise yet strong approach for distributed decomposition of TSRO scheduling problems in power systems.

In view of the above research gaps, this letter proposes a novel projection and decomposition (P&D) approach for multi-agent coordinated scheduling problems under uncertainty. Compared to the traditional distributed TSRO method, the proposed P&D approach identifies a high-quality solution with low computational cost and strong scalability.

II. General Model of Centralized TSRO

The centralized TSRO model for multi-agent scheduling with uncertainty in power systems is written in matrices as (1). The compact model M1 minimizes the total operation objective for system scheduling in the predicted nominal scenario, and enforces the power balance under uncertainty. The feasibility of these basic scheduling plans in the nominal scenario is ensured for all uncertainty scenarios, thus keeping the safe operation of the system [

10].

M1:minx1,x2,,xmΦ(x1,x2,,xm)i=1mciTxis.t.  Bixiei    i       i=1mCixif       Φ(x1,x2,,xm)=maxu={u1,u2,,um}uiΩi,iminy={y1,y2,,ym}ξ={ξ1,ξ2,,ξm}i=1m1Tξi       Giyi+ξihi-Liui-Mixi    i       i=1mKiyig (1)

where xi denotes the first-stage variables of agent i before uncertainty, that is, the basic scheduling plans, including the binary commitments and the continuous unit outputs; ciTxi denotes the local operation objective of agent i in the nominal scenario, and the total operation objective of all agents i=1,2,,m is minimized in the first-stage optimization. The second line of M1 represents the local constraints of agent i in the nominal scenario, mainly containing the unit operation and network security within its jurisdiction, and the third line denotes the coupling constraints among agents, such as the consistency for shared line power. In the fourth line,Φ(x1,x2,,xm) denotes the second-stage optimization model; and yi denotes the continuous recourse variables of agent i after uncertainty. The second-stage recourse plans modify xi to hedge the uncertainty. ξi denotes the non-negative slack variables introduced to characterize power imbalance with two parts (ξi+,ξi-). Note that when there is an inflow power imbalance, ξi->0 and ξi+=0; otherwise, when an outflow power imbalance emerges, ξi+>0 and ξi-=0. ui denotes the uncertainty decision variables of agent i, e.g., the uncertainty in source-load power; and Ωi denotes the uncertainty set, such as the boxed or polyhedral set. The fifth and sixth lines represents the local constraints of agent i and the coupling constraints between multiple agents under uncertainty, respectively. Bi, ei, Ci, f, Gi, hi, Li, Mi, Ki, and g are the constant matrices for the first- and second-stage constraints.

III. Proposed P&D Approach

Considering that a centralized TSRO model with coupling constraints in both stages cannot be decoupled with guaranteed optimality in distributed optimization, this section proposes a P&D approach to tackle this problem.

A. Feasible Region Projection of TSRO Model

For M1, the decision-making of the second stage intrinsically checks the scheduling feasibility under uncertainty, indicating that with the first-stage solution xi, there exists yi for any uncertainty scenario that ensures ξi=0. Hence, M1 is recast as:

M2:minx1,x2,,xmΦ(x1,x2,,xm)i=1mciTxis.t.  Bixiei    i       i=1mCixif       Φ(x1,x2,,xm)=uiΩi,i,(y1,y2,,ym):             Mixi+Giyi+Liuihi,Kiyig (2)

Compared to M1, the slack variables ξi are omitted in M2. The max-min bi-level model in the second stage is described as a multi-dimensional feasible region of the recourses, which is defined by the inner-layer linear programming constraints with respect to uncertainty. According to the Fourier-Motzkin elimination [

11], such a region in the max-min optimization for each agent represented by the constraints in the fifth line of M2 is projected to the domain formed by xi in the first stage and ui, so as to eliminate the second-stage variables yi [12]. After eliminating yi, the linear constraints can be generated as:

AiXxi+AiUuiki (3)

where AiX, AiU, and ki are the parameter matrices produced in the elimination; and i identifies the constraint set generated by the corresponding agent. The general procedures are given in [

13] to clarify how to derive such kinds of constraints in (3) by the Fourier-Motzkin elimination. Some accelerating policies such as redundant constraint filtering could be employed to enhance the computational performance of the Fourier-Motzkin elimination [14]. Furthermore, considering that the recourses and second-stage constraints are built by series, the Fourier-Motzkin elimination can uniformly handle each constraint series, thereby further improving its performance.

It is worth mentioning that the last term in (2) essentially represents the coupling constraints among agents only on their coupling variables in yi under uncertainty, and it does not involve any internal variables of other agents. Such types of constraints are shared among the neighboring agents, e.g., the consistency constraints on tie-line variables between different areas in power systems. Therefore, as agent i knows these coupling variables and constraints with its neighbors, the projection for the second-stage optimization could be conducted independently by agent i without the privacy of its neighboring agents, and each agent can obtain the linear constraints in (3). A simple multi-agent optimization model with coupling constraints is provided in [

13] to clarify the elimination procedures.

According to the property of the Fourier-Motzkin elimination, we can get all feasible solutions of a linear programming model from the solution results fitting these linear constraints generated by the Fourier-Motzkin elimination with two presuppositions [

15]: ① the optimization of the original mathematical programming has a nonempty set; and ② the constraint coefficients of the original programming are not all zero. The specific meaning of the feasible region projection for this model is that for all agents, if (xi, ui) denotes a solution set filling (3), yi certainly exists so that (xi, yi, ui) is a set of solutions falling into the feasible region characterized in the third and fourth lines in M2. After the projection, M2 is equivalently expressed as:

M3:minx1,x2,,xmΦ(x1,x2,,xm)i=1mciTxis.t.  Bixiei    i       i=1mCixif       Φ(x1,x2,,xm)={uiΩi,i,AiXxi+AiUuiki} (4)

M3 is still intractable because Φ(x1,x2,,xm) contains infinite uncertainty scenarios. As ui is exogenous, to ensure (3) valid for uiΩi, the only need is to meet:

AiXxi+maxuiΩiAiUuiki (5)

In the above formulation, the rth dimensional maximum optimization for ui is derived as:

maxui[AiU]rui s.t.  Ωi={ui|CiUuidiU:[δi]r0} (6)

where CiU and diU denote the coefficient matrices related to Ωi of agent i; and [δi]r denotes dual variables for the rth dimensional constraint. In practice, the max objectives AiUui can be directly determined by setting ui as the boundaries of the uncertainty interval. For example, if AiU>0 and ui takes the upper boundaries, (5) can be satisfied, then (3) must hold for uiΩi. Without loss of generality, the duality expression for each dimension in (5) identifies the same max values.

After substituting the dual formulation of the max optimization for each dimension in (6), M3 is transformed as:

M4:minx1,x2,,xm,δ1,δ2,,δmi=1mciTxis.t.  Bixiei    i       i=1mCixif       AiXxi+δiTdiUki    i       CiUδiT=AiU    i (7)

where δi collects [δi]r by columns for all dimensions. As the Fourier-Motzkin elimination is an exact equivalent procedure, the original min-max-min TSRO model, i.e., M1, is equivalent to a mixed-integer linear programming (MILP) problem M4 after feasible region projection. The solution of M4 does not need any iteration compared with M1, which is convenient for the engineering implementation.

B. Dual Decomposition of Equivalent MILP Problem

The MILP problem M4 after projection is refined into a compact formulation P1.

P1:minx1,x2,,xmi=1mciTxis.t.  i=1mAixib        xiXi ,Xi={xi: Dixidi}    i (8)

where δi in M4 is integrated as the decision variables of agent i to form the new variable vector xi for each agent in P1; Ai, b, Di, and di are the constant matrices or vectors for the coupling and local constraints in the refined formulation; Xi is the compact set of all local constraints for agent i in M4, i.e., the first, third, and fourth constraints in (7). The second line of (8) represents the corresponding coupling constraints, that is, the second constraint in (7), where the dimension of b is q.

The traditional ADMM and ATC method are less effective for an MILP because such a non-convex problem no longer meets the optimality and convergence conditions of the traditional decomposition approaches [

8]. Hence, a dual decomposition algorithm is deployed for the decoupling of the projected MILP problem to guarantee the solution convergence and quality [16]. The main steps include dualizing the coupling constraints using the Lagrange multipliers and handling the Lagrangian dual programming to obtain the dual variables. Then, the primal solution is recovered by solving the local MILP of each agent given the dual variables. This guarantees the recovered primal solution satisfies the local constraints but does not ensure that the obtained primal solution satisfies the coupling constraints because of the non-convexity of MILP. Hence, the dual subgradient algorithm iteratively generates the tentative solution and updates the dual variables by each agent in a decentralized manner. By suitably averaging these tentative solutions during iterations, we could identify one solution that fits both the coupling and local constraints via the iteration between the primal and dual problems. Besides, to enforce the feasibility of the obtained dual solution, the constraint tightening strategy is adopted to restore the feasible primary solution, which shrinks the duality gap and guarantees the solution quality. A tightened version for P1 is defined as P2 after introducing the non-negative tightening coefficients ρ.

P2:minx1,x2,,xmi=1mciTxis.t.  i=1mAixib-ρ    ρ0,xiXi,Xi={xi: Dixidi},i (9)

Its Lagrangian dual expression is formulated as:

maxλ0(-λT(b-ρ))+i=1mminxiXi(ciT+λTAi)xi (10)

where λ denotes the q-dimensional dual variables; and the definition of ρ is as follows: finding the maximal difference in the contributions of feasible solution to the coupling constraints.

[ρ]j=qmaxi{1,2,,m}maxxiXi[Ai]jxi-minxiXi[Ai]jxi (11)

where j is the index of the matrix line in the coupling constraints.

For traditional distributed scheduling, a coordination center is set up to collect the optimal results of all shared variables among agents. Their contributions to the coupling constraints help update ρ and λ, which are broadcast to all agents for the next iteration. In fact, it is a burden to set a coordination center for the whole network in the real-world cases. A max consensus algorithm is further devised to assign λ and ρ to avert a coordination center. The detailed steps of the decentralized algorithm are listed as below.

Algorithm: decentralized dual decomposition algorithm for P1
Step 1: Initialize k=0, set λi(0)=0, ρi(0)=0, si+(0)=-, si-(0)=+, for all i=1,2,,m
Step 2: Repeat
Step 3: For i=1,2,,m do
Step 4: ψi(k)=jNiαij(k)λj(k)
Step 5: xi(k+1)=arg minxiXi(ciT+ψi(k)TAi)xi
Step 6: ϕi(k)=maxjNi{ρj(k)}
Step 7: si+(k+1)=max{si+(k),Aixi(k+1)}
Step 8: si-(k+1)=min{si-(k),Aixi(k+1)}
Step 9: ρi(k+1)=max{ϕi(k),q(si+(k+1)-si-(k+1))}
Step 10: λi(k+1)=[ψi(k)+β(k)(Aixi(k+1)-(b-ρi(k+1))/m]+
Step 11: k=k+1
Step 12: Until x1(k),x2(k),,xm(k) satisfy the coupling constraints

In Step 4, agent i gathers λj(k) of its connected agent j, and the average ψi(k) for λj(k) is constructed by using the weighted factor αij(k); and Ni denotes the group of the neighboring agents for agent i. The dual variables of agent i are fixed as ψi(k), and the min optimization is performed in Step 5 to reap a tentative primary solution xi(k+1). In Steps 6-9, each agent refines the coefficient ρi(k+1) via the max and min optimizations si+(k+1) and si-(k+1) of components ϕi(k), respectively. Specifically, Steps 7 and 8 obtain the worst contributions of Aixi(k+1) to the coupling constraints by the tentative solution xi(k+1), and a tightening coefficient ρi(k+1) is derived as the maximum one of ρi(k) and ρj(k) in Steps 6 and 9. In Step 10, the dual variables are updated along the gradients via xi(k+1); β(k) denotes the step size whose selection rules can be found in [

17] to ensure the convergence property of this algorithm; and [·]+ denotes a projection to the non-negative quadrant with q-dimensions. Note that Steps 4 and 6 realize the decentralization of the dual decomposition, and each agent does not need to disclose its privacy information, but only the tentative λi(k) and ρi(k).

Note that the decentralized solution converges to the dual optimal solution in finite iterations with the dual subgradient algorithm. Going into detail, the dual variable sequence λi(k) created by Step 10 guarantees the convergence to the dual optimal solution. Such a dual finite-time convergence proposition for MILP models is presented in [

17]. Since a duality gap for the solutions of the primal and dual programming will exist due to non-convexity, the coupling constraint tightening in P2 recovers a feasible primal solution using the dual solution of the tightened primal problem [16]. The primal finite-time feasibility property has been proven (see Theorem 1 in [18]). Meanwhile, Theorem 2 in [18] further indicates the optimal performance property. Therefore, the dual decomposition algorithm could derive a decentralized feasible solution for the primal programming that fills tight performance bounds within finite iterations. Interested researchers can refer to [18] for elaborated mathematical proofs on the algorithmic optimal and convergent properties of MILP problems.

C. Remarks on Proposed P&D Approach

In this study, we consider the source-load power uncertainty shown in (1) for multi-agent coordinated scheduling of power systems. Note that there would be uncertainty in the availability of power generation or consumption resources in practice, while the proposed P&D approach is not limited considering certain kinds of uncertainty. In other words, the scheduling model in (1) is general for multiple uncertainty. For instance, when we further consider the N-k output uncertainty of power units [

19], which is a typical uncertainty in the availability of generation, we only need to incorporate the uncertainty set of unit outage into ui and Ωi, and this does not affect the essential formulation of the mathematical model in (1). Hence, the proposed P&D approach is still applicable under this condition.

It is noteworthy that this paper focuses on how to achieve an effective decentralized solution to a centralized RO problem and ensure that the derived decentralized results are as close as possible to the optimal results of the centralized method. From the perspective of game theory, the relationship of the multi-agent objectives in this paper belongs to the cooperative game. There also exists the condition of multiple agents with conflicting objectives in scheduling. Currently, some effective models or methods have been investigated to address this problem, such as the Stackelberg game model [

20] and the Nash equilibrium model, but this topic goes beyond the scope of this work.

IV. Numerical Tests

Three cases of multi-area scheduling are evaluated to verify the superiority and scalability of the proposed P&D approach. Case 1 is a small-scale 2-area 12-bus system with one wind farm integrated into each area and an interconnected tie-line. Case 2 and Case 3 both adopt a 1441-bus system with practical size, which are divided into 4 areas and 8 areas, respectively. The real 1441-bus system contains 130 wind farms, 320 conventional generators, 2057 internal lines, and 23 tie-lines connecting multiple areas with a voltage level of 500 kV. Detailed data of the two test systems are available in [

9] and [2], respectively. For each case, the proposed P&D approach (A0), the centralized TSRO method (A1) [10], and the distributed TSRO method based on ADMM (A2) [4] are used for comparison. The detailed centralized TSRO model before distributed decomposition is formulated in [13]. All TSRO models in compared methods are tackled by the C&CG algorithm with strong duality. Besides, 1000 stochastic scenarios are generated for Case 1 and 20000 for Cases 2 and 3, to check the feasibility of their basic plans, and those scenarios with power balance are recorded. The optimization horizon and the period interval are 24 hours and 1 hour, respectively. Numerical tests are implemented on a laptop under 2.30 GHz i5-CPU and 8 GB RAM, and CPLEX 12.8 solves all MILPs. The optimal results of different cases using various robust scheduling methods are listed in Table I.

TABLE I  Optimal Results of Different Cases
CaseMethodCost ($)IterationTime (s)Number of scenariosProblem scale
Number of constraintsNumber of variables

Case 1

(2-area 12-bus system)

A1 150243 N/A 4.9 1000 22646 10638
A0 150280 23 37.6 1000 28126 7352†
A2 152326 78 1146.5 1000 22694 10710‡

Case 2

(4-area 1441-bus system)

A1 29376412 N/A 100.4 20000 1248106 465184
A0 29382876 27 95.2 20000 1512704 303374†
A2 29768991 132 3880.6 20000 1248298 465568‡

Case 3

(8-area 1441-bus system)

A1 29376412 N/A 100.4 20000 1248106 465184
A0 29381995 30 63.3 20000 1512822 303402†
A2 29617300 146 3433.7 20000 1248490 465952‡

Note:   † represents the total constraints and variables of the models after elimination for all agents; and ‡ represents the total constraints and variables of the decoupled models after ADMM for all agents.

In Case 1, it is noteworthy that the total operation cost of A0 is approximately equal to the referring centralized solution in A1, and the relative optimality gap is 0.025%, whereas the gap for A2 is 1.39%. The dual decomposition algorithm realizes the effective decoupling of the non-convex MILP problem, and the constraint tightening contracts the duality gap that ensures the solution quality. The number of distributed iterations and the total solution time of A2 are 3.4 and 30.5 times of those in A0, respectively. The reason is that the dual decomposition algorithm updates the dual variables in the direction of gradients, which accelerates the convergence of the problem solution. The update of the multipliers and the penalty factors in ADMM depends on the empirical values, resulting in a low convergence speed. In addition, each agent directly solves the MILP-type robust scheduling model after projection, which greatly reduces the computational time for each distributed iteration.

For larger-scale testing systems in Case 2 and Case 3, the optimality gap between A0 and A1 is further reduced. In particular, the optimality gaps for Case 2 and Case 3 are 0.022% and 0.019%, respectively, indicating that the optimal cost is closer to the centralized objective with more agents. Besides, the iteration of A0 increases slowly, hence its computational performance is relatively stable. The solution within minutes in A0 fully meets the scheduling requirements, while the iteration number of A2 is 4.8 times of that by A0, and the total run time even reaches 40-54 times. It can also be observed that the increase of the modeling scale leads to a significant rise in the iterations and the run time of A2, and the hour-level computational time hinders its practical application. Although the total number of constraints after elimination in A0 increases compared with A1 and A2, the dual decomposition further divides the entire problem into smaller local scheduling problems for each agent, hence it will not further pose huge computational challenges to the decentralized implementation by agents. Therefore, A0 enjoys stable computational performance and high scalability in terms of agent number and system size, and these strengths are more prominent compared with A2. Besides, for all cases, the three methods (A0, A1, and A2) derive the optimal power plans that guarantee 100% feasibility to all stochastic scenarios, indicating the robustness of their solution results.

V. Conclusion

This letter proposes a novel P&D approach for the robust coordinated scheduling of multiple agents under uncertainty. Numerical tests denote that the proposed P&D approach overcomes the defects of low computational efficiency, slow convergence, and suboptimality in the traditional distributed TSRO method. The merits of the proposed P&D approach in optimality and computational performance are more notable with the enlargement of the system scale and agent quantity. Note that the proposed P&D approach is a universal distributed TSRO method that can be easily extended to numerous coordinated scheduling problems such as multi-area unit commitment, transmission-distribution coordination, energy management of networked microgrids, peer-to-peer trading in prosumers, and integrated energy system dispatch. Future research will consider extending the P&D approach for distributed TSRO methods with binary recourses and endogenous uncertainty. Besides, further research is needed on distributed scheduling of multiple agents with conflicting objectives under uncertainty.

References

1

J. Wei, Y. Zhang, J. Wang et al., “Decentralized demand management based on alternating direction method of multipliers algorithm for industrial park with CHP units and thermal storage,” Journal of Modern Power Systems and Clean Energy, vol. 10, no. 1, pp. 120-130, Jan. 2022. [Baidu Scholar] 

2

Z. Li, W. Wu, M. Shahidehpour et al., “Adaptive robust tie-line scheduling considering wind power uncertainty for interconnected power systems,” IEEE Transactions on Power Systems, vol. 31, no. 4, pp. 2701-2713, Jul. 2016. [Baidu Scholar] 

3

H. Gao, J. Liu, L. Wang et al., “Decentralized energy management for networked microgrids in future distribution systems,” IEEE Transactions on Power Systems, vol. 33, no. 4, pp. 3599-3610, Jul. 2018. [Baidu Scholar] 

4

H. Qiu and F. You, “Decentralized-distributed robust electric power scheduling for multi-microgrid systems,” Applied Energy, vol. 269, p. 115146, Jul. 2020. [Baidu Scholar] 

5

Z. Chen, Z. Li, C. Guo et al., “Fully distributed robust reserve scheduling for coupled transmission and distribution systems,” IEEE Transactions on Power Systems, vol. 36, no. 1, pp. 169-182, Jan. 2021. [Baidu Scholar] 

6

Y. Ji, Q. Xu, and Y. Xia, “Distributed robust energy and reserve dispatch for coordinated transmission and active distribution systems,” Journal of Modern Power Systems and Clean Energy, vol. 11, no. 5, pp. 1494-1506, Sept. 2023. [Baidu Scholar] 

7

X. Liang, Z. Li, W. Huang et al., “Relaxed alternating direction method of multipliers for hedging communication packet loss in integrated electrical and heating system,” Journal of Modern Power Systems and Clean Energy, vol. 8, no. 5, pp. 874-883, Sept. 2020. [Baidu Scholar] 

8

N. Michelena, H. Park, and P. Papalambros, “Convergence properties of analytical target cascading,” AIAA Journal, vol. 41, no. 5, pp. 897-905, May 2003. [Baidu Scholar] 

9

Z. Li, M. Shahidehpour, W. Wu et al., “Decentralized multiarea robust generation unit and tie-line scheduling under wind power uncertainty,” IEEE Transactions on Sustainable Energy, vol. 6, no. 4, pp. 1377-1388, Oct. 2015. [Baidu Scholar] 

10

N. G. Cobos, J. M. Arroyo, N. Alguacil et al., “Robust energy and reserve scheduling considering bulk energy storage units and wind uncertainty,” IEEE Transactions on Power Systems, vol. 33, no. 5, pp. 5206-5216, Sept. 2018. [Baidu Scholar] 

11

J. Imbert, “Fourier’s elimination: which to choose?” in Proceedings of Principles and Practice of Constraint Programming, Newport, USA, Jun. 1993, pp. 117-129. [Baidu Scholar] 

12

J. Zhen, D. Hertog, and M. Sim, “Adjustable robust optimization via Fourier-Motzkin elimination,” Operations Research, vol. 66, no. 4, pp. 1086-1100, Jul. 2018. [Baidu Scholar] 

13

H. Qiu. (2023, Oct.). Supplementary. [Online]. Available: https://entuedu-my.sharepoint.com/:b:/g/personal/haifeng_qiu_staff_main_ntu_edu_sg/ETahe6g9sy9Eoy005oaY7CEBZyF4B65aDjrfRYD_jGIv4w?e=tdkKI8 [Baidu Scholar] 

14

Y. Liu, L. Wu, and J. Li, “A fast LP-based approach for robust dynamic economic dispatch problem: a feasible region projection method,” IEEE Transactions on Power Systems, vol. 35, no. 5, pp. 4116-4119, Sept. 2020. [Baidu Scholar] 

15

G. B. Dantzig and B. C. Eaves, “Fourier-Motzkin elimination and its dual,” Journal of Combinatorial Theory (A), vol. 14, no. 3, pp. 288-297, May 1973. [Baidu Scholar] 

16

H. Qiu, A. Vinod, S. Lu et al., “Decentralized mixed-integer optimization for robust integrated electricity and heat scheduling,” Applied Energy, vol. 350, p. 121693, Nov. 2023. [Baidu Scholar] 

17

A. Falsone and M. Prandini, “A distributed dual proximal minimization algorithm for constraint-coupled optimization problems,” IEEE Control Systems Letters, vol. 5, no. 1, pp. 259-264, Jan. 2021. [Baidu Scholar] 

18

A. Falsone, K. Margellos, and M. Prandini, “A distributed iterative algorithm for multi-agent MILPs: finite-time feasibility and performance characterization,” IEEE Control Systems Letters, vol. 2, no. 4, pp. 563-568, Oct. 2018. [Baidu Scholar] 

19

H. Qiu, W. Gu, W. Sheng et al., “Resilience-oriented multistage scheduling for power grids considering nonanticipativity under tropical cyclones,” IEEE Transactions on Power Systems, vol. 38, no. 4, pp. 3254-3267, Jul. 2023. [Baidu Scholar] 

20

H. Gao, H. Pan, R. An et al., “Bi-level multi-leader multi-follower stackelberg game model for multi-energy retail package optimization,” Journal of Modern Power Systems and Clean Energy, vol. 12, no. 1, pp. 225-237, Jan. 2024. [Baidu Scholar]