A New Approach to Electricity Market Clearing With Uniform Purchase Price and Curtailable Block Orders

The European market clearing problem is characterized by a set of heterogeneous orders and rules that force the implementation of heuristic and iterative solving methods. In particular, curtailable block orders and the uniform purchase price (UPP) pose serious difficulties. A block is an order that spans over multiple hours, and can be either fully accepted or fully rejected. The UPP prescribes that all consumers pay a common price, i.e., the UPP, in all the zones, while producers receive zonal prices, which can differ from one zone to another. The market clearing problem in the presence of both the UPP and block orders is a major open issue in the European context. The UPP scheme leads to a non-linear optimization problem involving both primal and dual variables, whereas block orders introduce multi-temporal constraints and binary variables into the problem. As a consequence, the market clearing problem in the presence of both blocks and the UPP can be regarded as a non-linear integer programming problem involving both primal and dual variables with complementary and multi-temporal constraints. The aim of this paper is to present a non-iterative and heuristic-free approach for solving the market clearing problem in the presence of both curtailable block orders and the UPP. The solution is exact, with no approximation up to the level of resolution of current market data. By resorting to an equivalent UPP formulation, the proposed approach results in a mixed-integer linear program, which is built starting from a non-linear integer bilevel programming problem. Numerical results using real market data are reported to show the effectiveness of the proposed approach. The model has been implemented in Python, and the code is freely available on a public repository.Comment: 15 pages, 7 figure

A Consensus-ADMM Approach for Strategic Generation Investment in Electricity Markets

This paper addresses a multi-stage generation investment problem for a strategic (price-maker) power producer in electricity markets. This problem is exposed to different sources of uncertainty, including short-term operational (e.g., rivals' offering strategies) and long-term macro (e.g., demand growth) uncertainties. This problem is formulated as a stochastic bilevel optimization problem, which eventually recasts as a large-scale stochastic mixed-integer linear programming (MILP) problem with limited computational tractability. To cope with computational issues, we propose a consensus version of alternating direction method of multipliers (ADMM), which decomposes the original problem by both short- and long-term scenarios. Although the convergence of ADMM to the global solution cannot be generally guaranteed for MILP problems, we introduce two bounds on the optimal solution, allowing for the evaluation of the solution quality over iterations. Our numerical findings show that there is a trade-off between computational time and solution quality

Achieving an optimal trade-off between revenue and energy peak within a smart grid environment

We consider an energy provider whose goal is to simultaneously set revenue-maximizing prices and meet a peak load constraint. In our bilevel setting, the provider acts as a leader (upper level) that takes into account a smart grid (lower level) that minimizes the sum of users' disutilities. The latter bases its decisions on the hourly prices set by the leader, as well as the schedule preferences set by the users for each task. Considering both the monopolistic and competitive situations, we illustrate numerically the validity of the approach, which achieves an 'optimal' trade-off between three objectives: revenue, user cost, and peak demand

Evaluating Resilience of Electricity Distribution Networks via A Modification of Generalized Benders Decomposition Method

This paper presents a computational approach to evaluate the resilience of electricity Distribution Networks (DNs) to cyber-physical failures. In our model, we consider an attacker who targets multiple DN components to maximize the loss of the DN operator. We consider two types of operator response: (i) Coordinated emergency response; (ii) Uncoordinated autonomous disconnects, which may lead to cascading failures. To evaluate resilience under response (i), we solve a Bilevel Mixed-Integer Second-Order Cone Program which is computationally challenging due to mixed-integer variables in the inner problem and non-convex constraints. Our solution approach is based on the Generalized Benders Decomposition method, which achieves a reasonable tradeoff between computational time and solution accuracy. Our approach involves modifying the Benders cut based on structural insights on power flow over radial DNs. We evaluate DN resilience under response (ii) by sequentially computing autonomous component disconnects due to operating bound violations resulting from the initial attack and the potential cascading failures. Our approach helps estimate the gain in resilience under response (i), relative to (ii)