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 New General-Purpose Algorithm for Mixed-Integer Bilevel Linear Programs

International audienceBilevel optimization problems are very challenging optimization models arising in many important practical contexts, including pricing mechanisms in the energy sector, airline and telecommunication industry, transportation networks, critical infrastructure defense, and machine learning. In this paper, we consider bilevel programs with continuous and discrete variables at both levels, with linear objectives and constraints (continuous upper level variables, if any, must not appear in the lower level problem). We propose a general-purpose branch-and-cut exact solution method based on several new classes of valid inequalities, which also exploits a very effective bilevel-specific preprocessing procedure. An extensive computational study is presented to evaluate the performance of various solution methods on a common testbed of more than 800 instances from the literature and 60 randomly generated instances. Our new algorithm consistently outperforms (often by a large margin) alternative state-of-the-art methods from the literature, including methods exploiting problem-specific information for special instance classes. In particular, it solves to optimality more than 300 previously unsolved instances from the literature. To foster research on this challenging topic, our solver is made publicly available online

A New General-Purpose Algorithm for Mixed-Integer Bilevel Linear Programs

Bilevel optimization problems are very challenging optimization models arising in many important practical contexts, including pricing mechanisms in the energy sector, airline and telecommunication industry, transportation networks, critical infrastructure defense, and machine learning. In this paper, we consider bilevel programs with continuous and discrete variables at both levels, with linear objectives and constraints (continuous upper level variables, if any, must not appear in the lower level problem). We propose a general-purpose branch-and-cut exact solution method based on several new classes of valid inequalities, which also exploits a very effective bilevel-specific preprocessing procedure. An extensive computational study is presented to evaluate the performance of various solution methods on a common testbed of more than 800 instances from the literature and 60 randomly generated instances. Our new algorithm consistently outperforms (often by a large margin) alternative state-of-the-art methods from the literature, including methods exploiting problem-specific information for special instance classes. In particular, it solves to optimality more than 300 previously unsolved instances from the literature. To foster research on this challenging topic, our solver is made publicly available online

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)

Fuzzy Bilevel Optimization

In the dissertation the solution approaches for different fuzzy optimization problems are presented. The single-level optimization problem with fuzzy objective is solved by its reformulation into a biobjective optimization problem. A special attention is given to the computation of the membership function of the fuzzy solution of the fuzzy optimization problem in the linear case. Necessary and sufficient optimality conditions of the the convex nonlinear fuzzy optimization problem are derived in differentiable and nondifferentiable cases. A fuzzy optimization problem with both fuzzy objectives and constraints is also investigated in the thesis in the linear case. These solution approaches are applied to fuzzy bilevel optimization problems. In the case of bilevel optimization problem with fuzzy objective functions, two algorithms are presented and compared using an illustrative example. For the case of fuzzy linear bilevel optimization problem with both fuzzy objectives and constraints k-th best algorithm is adopted.:1 Introduction 1 1.1 Why optimization . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 1 1.2 Fuzziness as a concept . . . . . . . . . . . . . . . . . . . . .. . . . . . . 2 1.3 Bilevel problems . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6 2 Preliminaries 11 2.1 Fuzzy sets and fuzzy numbers . . . . . . . . . . . . . . . . . . . . . 11 2.2 Operations . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15 2.3 Fuzzy order . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16 2.4 Fuzzy functions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .17 3 Optimization problem with fuzzy objective 19 3.1 Formulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19 3.2 Solution method . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20 3.3 Local optimality . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 24 3.4 Existence of an optimal solution . . . . . . . . . . . . . . . . . . . . 25 4 Linear optimization with fuzzy objective 27 4.1 Main approach . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28 4.2 Example . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30 4.3 Optimality conditions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 33 4.4 Membership function value . . . . . . . . . . . . . . . . . . . . . . . . 34 4.4.1 Special case of triangular fuzzy numbers . . . . . . . . . . . . 36 4.4.2 Example . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .39 5 Optimality conditions 47 5.1 Differentiable fuzzy optimization problem . . . . . . . . . . .. . . . 48 5.1.1 Basic notions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 48 5.1.2 Necessary optimality conditions . . . . . . . . . . . . . . . . . . .. 49 5.1.3 Suffcient optimality conditions . . . . . . . . . . . . . . . . . . . . . . 49 5.2 Nondifferentiable fuzzy optimization problem . . . . . . . . . . . . 51 5.2.1 Basic notions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 51 5.2.2 Necessary optimality conditions . . . . . . . . . . . . . . . . . . . 52 5.2.3 Suffcient optimality conditions . . . . . . . . . . . . . . . . . . . . . . 54 5.2.4 Example . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55 6 Fuzzy linear optimization problem over fuzzy polytope 59 6.1 Basic notions . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 60 6.2 The fuzzy polytope . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .63 6.3 Formulation and solution method . . . . . . . . . . . . . . . . . . .. . 65 6.4 Example . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 69 7 Bilevel optimization with fuzzy objectives 73 7.1 General formulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 74 7.2 Solution approaches . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .74 7.3 Yager index approach . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 76 7.4 Algorithm I . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 77 7.5 Membership function approach . . . . . . . . . . . . . . . . . . . . . . .78 7.6 Algorithm II . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . .80 7.7 Example . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 81 8 Linear fuzzy bilevel optimization (with fuzzy objectives and constraints) 87 8.1 Formulation . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87 8.2 Solution approach . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 88 8.3 Algorithm . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90 8.4 Example . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 91 9 Conclusions 95 Bibliography 9

Intersection disjunctions for reverse convex sets

We present a framework to obtain valid inequalities for optimization problems constrained by a reverse convex set, which is defined as the set of points in a polyhedron that lie outside a given open convex set. We are particularly interested in cases where the closure of the convex set is either non-polyhedral, or is defined by too many inequalities to directly apply disjunctive programming. Reverse convex sets arise in many models, including bilevel optimization and polynomial optimization. Intersection cuts are a well-known method for generating valid inequalities for a reverse convex set. Intersection cuts are generated from a basic solution that lies within the convex set. Our contribution is a framework for deriving valid inequalities for the reverse convex set from basic solutions that lie outside the convex set. We begin by proposing an extension to intersection cuts that defines a two-term disjunction for a reverse convex set. Next, we generalize this analysis to a multi-term disjunction by considering the convex set's recession directions. These disjunctions can be used in a cut-generating linear program to obtain disjunctive cuts for the reverse convex set.Comment: 24 page

On SOCP-based disjunctive cuts for solving a class of integer bilevel nonlinear programs

We study a class of integer bilevel programs with second-order cone constraints at the upper-level and a convex-quadratic objective function and linear constraints at the lower-level. We develop disjunctive cuts (DCs) to separate bilevel-infeasible solutions using a second-order-cone-based cut-generating procedure. We propose DC separation strategies and consider several approaches for removing redundant disjunctions and normalization. Using these DCs, we propose a branch-and-cut algorithm for the problem class we study, and a cutting-plane method for the problem variant with only binary variables. We present an extensive computational study on a diverse set of instances, including instances with binary and with integer variables, and instances with a single and with multiple linking constraints. Our computational study demonstrates that the proposed enhancements of our solution approaches are effective for improving the performance. Moreover, both of our approaches outperform a state-of-the-art generic solver for mixed-integer bilevel linear programs that is able to solve a linearized version of our binary instances.Comment: arXiv admin note: substantial text overlap with arXiv:2111.0682

Models and algorithms for decomposition problems

This thesis deals with the decomposition both as a solution method and as a problem itself. A decomposition approach can be very effective for mathematical problems presenting a specific structure in which the associated matrix of coefficients is sparse and it is diagonalizable in blocks. But, this kind of structure may not be evident from the most natural formulation of the problem. Thus, its coefficient matrix may be preprocessed by solving a structure detection problem in order to understand if a decomposition method can successfully be applied. So, this thesis deals with the k-Vertex Cut problem, that is the problem of finding the minimum subset of nodes whose removal disconnects a graph into at least k components, and it models relevant applications in matrix decomposition for solving systems of equations by parallel computing. The capacitated k-Vertex Separator problem, instead, asks to find a subset of vertices of minimum cardinality the deletion of which disconnects a given graph in at most k shores and the size of each shore must not be larger than a given capacity value. Also this problem is of great importance for matrix decomposition algorithms. This thesis also addresses the Chance-Constrained Mathematical Program that represents a significant example in which decomposition techniques can be successfully applied. This is a class of stochastic optimization problems in which the feasible region depends on the realization of a random variable and the solution must optimize a given objective function while belonging to the feasible region with a probability that must be above a given value. In this thesis, a decomposition approach for this problem is introduced. The thesis also addresses the Fractional Knapsack Problem with Penalties, a variant of the knapsack problem in which items can be split at the expense of a penalty depending on the fractional quantity