TY - JOUR A1 - Ambrosius, Mirjam A1 - Grimm, Veronika A1 - Kleinert, Thomas A1 - Liers, Frauke A1 - Schmidt, Martin A1 - Zöttl, Gregor T1 - Endogenous Price Zones and Investment Incentives in Electricity Markets: An Application of Multilevel Optimization with Graph Partitioning JF - Energy Economics N2 - In the course of the energy transition, load and supply centers are growing apart in electricity markets worldwide, rendering regional price signals even more important to provide adequate locational investment incentives. This paper focuses on electricity markets that operate under a zonal pricing market design. For a fixed number of zones, we endogenously derive the optimal configuration of price zones and available transfer capacities on a network in order to optimally govern investment and production decisions in the long run. In a multilevel mixed-integer nonlinear model that contains a graph partitioning problem on the first level, we determine welfare-maximizing price zones and available transfer capacities for a given electricity market and analyze their impact on market outcomes. Using a generalized Benders decomposition approach developed in Grimm et al. (2019) and a problem-tailored scenario clustering for reducing the input data size, we are able to solve the model to global optimality even for large instances. We apply the approach to the German electricity market as an example to examine the impact of optimal zoning on key performance indicators such as welfare, generation mix and locations, or electricity prices. It turns out that even for a small number of price zones, an optimal configuration of zones induces a welfare level that almost approaches the first best. KW - Electricity Markets KW - Price Zones KW - Investment Incentives KW - Multilevel Optimization KW - Graph Partitioning Y1 - 2018 IS - 92 ER - TY - INPR A1 - Liers, Frauke A1 - Martin, Alexander A1 - Merkert, Maximilian A1 - Mertens, Nick A1 - Michaels, Dennis T1 - Solving Mixed-Integer Nonlinear Optimization Problems using Simultaneous Convexification - a Case Study for Gas Networks N2 - Solving mixed-integer nonlinear optimization problems (MINLPs) to global optimality is extremely challenging. An important step for enabling their solution consists in the design of convex relaxations of the feasible set. Known solution approaches based on spatial branch-and-bound become more effective the tighter the used relaxations are. Relaxations are commonly established by convex underestimators, where each constraint function is considered separately. Instead, a considerably tighter relaxation can be found via so-called simultaneous convexification, where convex underestimators are derived for more than one constraint function at a time. In this work, we present a global solution approach for solving mixed-integer nonlinear problems that uses simultaneous convexification. We introduce a separation method that relies on determining the convex envelope of linear combinations of the constraint functions and on solving a nonsmooth convex problem. In particular, we apply the method to quadratic absolute value functions and derive their convex envelopes. The practicality of the proposed solution approach is demonstrated on several test instances from gas network optimization, where the method outperforms standard approaches that use separate convex relaxations. KW - Mixed-Integer Nonlinear Programming KW - Simultaneous Convexification KW - Convex Envelope KW - Gas Network Optimization Y1 - 2020 ER - TY - INPR A1 - Kuchlbauer, Martina A1 - Liers, Frauke A1 - Stingl, Michael T1 - Adaptive bundle methods for nonlinear robust optimization N2 - Currently, there are few theoretical or practical approaches available for general nonlinear robust optimization. Moreover, the approaches that do exist impose restrictive assumptions on the problem structure. We present an adaptive bundle method for nonlinear and non-convex robust optimization problems with a suitable notion of inexactness in function values and subgradients. As the worst case evaluation requires a global solution to the adversarial problem, it is a main challenge in a general non-convex nonlinear setting. Moreover, computing elements of an epsilon-perturbation of the Clarke subdifferential in the l2-norm sense is in general prohibitive for this class of problems. In this article, instead of developing an entirely new bundle concept, we demonstrate how existing approaches, such as Noll's bundle method for non-convex minimization with inexact information (Computational and analytical mathematics 50: 555-592, 2013) can be modified to be able to cope with this situation. Extending the non-convex bundle concept to the case of robust optimization in this way, we prove convergence under two assumptions: Firstly, that the objective function is lower C1 and secondly, that approximately optimal solutions to the adversarial maximization problem are available. The proposed method is hence applicable to a rather general setting of nonlinear robust optimization problems. In particular, we do not rely on a specific structure of the adversary's constraints. The considered class of robust optimization problems covers the case that the worst-case adversary only needs to be evaluated up to a certain precision. One possibility to evaluate the worst case with the desired degree of precision is the use of techniques from mixed-integer linear programming (MIP). We investigate the procedure on some analytic examples. As applications, we study the gas transport problem under uncertainties in demand and in physical parameters that affect pressure losses in the pipes. Computational results for examples in large realistic gas network instances demonstrate the applicability as well as the efficiency of the method. Y1 - 2020 ER - TY - INPR A1 - Aigner, Kevin-Martin A1 - Clarner, Jan-Patrick A1 - Liers, Frauke A1 - Martin, Alexander T1 - Robust Approximation of Chance Constrained DC Optimal Power Flow under Decision-Dependent Uncertainty N2 - We propose a mathematical optimization model and its solution for joint chance constrained DC Optimal Power Flow. In this application, it is particularly important that there is a high probability of transmission limits being satisfied, even in the case of uncertain or fluctuating feed-in from renewable energy sources. In critical network situations where the network risks overload, renewable energy feed-in has to be curtailed by the transmission system operator (TSO). The TSO can reduce the feed-in in discrete steps at each network node. The proposed optimization model minimizes curtailment while ensuring that there is a high probability of transmission limits being maintained. The latter is modeled via (joint) chance constraints that are computationally challenging. Thus, we propose a solution approach based on the robust safe approximation of these constraints. Hereby, probabilistic constraints are replaced by robust constraints with suitably defined uncertainty sets constructed from historical data. The uncertainty sets are calculated by encompassing randomly drawn scenarios using the scenario approach proposed by Margellos et al. (IEEE Transactions on Automatic Control, 59 (2014)). The ability to discretely control the power feed-in then leads to a robust optimization problem with decision-dependent uncertainties, i.e. the uncertainty sets depend on decision variables. We propose an equivalent mixed-integer linear reformulation for box uncertainties with the exact linearization of bilinear terms. Finally, we present numerical results for different test cases from the Nesta archive, as well as for a real network. We consider the discrete curtailment of solar feed-in, for which we use real-world weather and network data. The experimental tests demonstrate the effectiveness of this method and run times are very fast. Moreover, on average the calculated robust solutions lead only to a small increase in curtailment, when compared to nominal solutions. KW - OR in energy KW - optimal power flow KW - chance constrained programming KW - robust optimization KW - decision-dependent uncertainty Y1 - 2020 ER - TY - JOUR A1 - Adelhütte, Dennis A1 - Aßmann, Denis A1 - Gonzàlez Grandòn, Tatiana A1 - Gugat, Martin A1 - Heitsch, Holger A1 - Liers, Frauke A1 - Henrion, René A1 - Nitsche, Sabrina A1 - Schultz, Rüdiger A1 - Stingl, Michael A1 - Wintergerst, David T1 - Joint model of probabilistic/robust (probust) constraints applied to gas network optimization N2 - Optimization tasks under uncertain conditions abound in many real-life applications. Whereas solution approaches for probabilistic constraints are often developed in case the uncertainties can be assumed to follow a certain probability distribution, robust approaches are usually used in case solutions are sought that are feasible for all realizations of uncertainties within some pre-defined uncertainty set. As many applications contain different types of uncertainties that require robust as well as probabilistic treatments, we deal with a class of joint probabilistic/robust constraints as its appears in optimization problems under uncertainty. Focusing on complex uncertain gas network optimization problems, we show the relevance of this class of problems for the task of maximizing free booked capacities in an algebraic model for a stationary gas network. We furthermore present approaches for their solution. Finally, we study the problem of controlling a transient system that is governed by the wave equation. The task consists in determining controls such that a certain robustness measure remains below some given upper bound, with high probability. KW - robust optimization KW - chance constraints KW - optimal control KW - spheric-radial decomposition Y1 - 2017 U6 - https://doi.org/10.1007/s10013-020-00434-y ER - TY - INPR A1 - Liers, Frauke A1 - Schewe, Lars A1 - Thürauf, Johannes T1 - Radius of Robust Feasibility for Mixed-Integer Problems N2 - For a mixed-integer linear problem (MIP) with uncertain constraints, the radius of robust feasibility (RRF) determines a value for the maximal “size” of the uncertainty set such that robust feasibility of the MIP can be guaranteed. The approaches for the RRF in the literature are restricted to continuous optimization problems. We first analyze relations between the RRF of a MIP and its continuous linear (LP) relaxation. In particular, we derive conditions under which a MIP and its LP relaxation have the same RRF. Afterward, we extend the notion of the RRF such that it can be applied to a large variety of optimization problems and uncertainty sets. In contrast to the setting commonly used in the literature, we consider for every constraint a potentially different uncertainty set that is not necessarily full-dimensional. Thus, we generalize the RRF to MIPs as well as to include “safe” variables and constraints, i.e., where uncertainties do not affect certain variables or constraints. In the extended setting, we again analyze relations between the RRF for a MIP and its LP relaxation. Afterward, we present methods for computing the RRF of LPs as well as of MIPs with safe variables and constraints. Finally, we show that the new methodologies can be successfully applied to the instances in the MIPLIB 2017 for computing the RRF. KW - Robust Optimization, Mixed-integer programming, Uncertainty sets, Robust feasibility Y1 - 2019 U6 - https://doi.org/10.1287/ijoc.2020.1030 PB - Informs Journal on Computing ER - TY - INPR A1 - Aßmann, Denis A1 - Liers, Frauke A1 - Stingl, Michael T1 - Decomposable Robust Two-Stage Optimization: An Application to Gas Network Operations Under Uncertainty N2 - We study gas network problems with compressors and control valves under uncertainty that can be formulated as two-stage robust optimization problems. Uncertain data are present in the physical parameters of the pipes as well as in the overall demand. We show how to exploit the special decomposable structure of the problem in order to reformulate the two-stage robust problem as a standard single-stage optimization problem. Since this structure is present in similar problems on e.g., water or direct current electricity networks, we investigate the consequences of the decomposable structure in an abstract setting: The right-hand side of the single-stage problem can be precomputed by solving a series of optimization problems and multiple elements of the right-hand side can be combined into one optimization task. In order to apply our results to gas network problems, we extend piecewise relaxations and preprocessing techniques to incorporate uncertain input data. The practical feasibility and effectiveness of our approach is demonstrated with benchmarks on realistic gas network instances. We observe large speedups due to the described aggregation method together with the developed preprocessing strategies. Furthermore, we are able to solve even comparably large gas network instances quickly for the price of slightly more conservative solutions. KW - robust optimization KW - gas networks KW - relaxations Y1 - 2017 ER - TY - INPR A1 - Aigner, Kevin-Martin A1 - Bärmann, Andreas A1 - Braun, Kristin A1 - Liers, Frauke A1 - Pokutta, Sebastian A1 - Schneider, Oskar A1 - Sharma, Kartikey A1 - Tschuppik, Sebastian T1 - Data-driven Distributionally Robust Optimization over Time N2 - Stochastic Optimization (SO) is a classical approach for optimization under uncertainty that typically requires knowledge about the probability distribution of uncertain parameters. As the latter is often unknown, Distributionally Robust Optimization (DRO) provides a strong alternative that determines the best guaranteed solution over a set of distributions (ambiguity set). In this work, we present an approach for DRO over time that uses online learning and scenario observations arriving as a data stream to learn more about the uncertainty. Our robust solutions adapt over time and reduce the cost of protection with shrinking ambiguity. For various kinds of ambiguity sets, the robust solutions converge to the SO solution. Our algorithm achieves the optimization and learning goals without solving the DRO problem exactly at any step. We also provide a regret bound for the quality of the online strategy which converges at a rate of $ O(\log T / \sqrt{T})$, where $T$ is the number of iterations. Furthermore, we illustrate the effectiveness of our procedure by numerical experiments on mixed-integer optimization instances from popular benchmark libraries and give practical examples stemming from telecommunications and routing. Our algorithm is able to solve the DRO over time problem significantly faster than standard reformulations. KW - distributionally robust optimization KW - learning over time KW - gradient descent Y1 - 2023 ER - TY - INPR A1 - Bernhard, Daniela A1 - Liers, Frauke A1 - Stingl, Michael T1 - A Gradient-Based Method for Joint Chance-Constrained Optimization with Continuous Distributions N2 - The input parameters of an optimization problem are often affected by uncertainties. Chance constraints are a common way to model stochastic uncertainties in the constraints. Typically, algorithms for solving chance-constrained problems require convex functions or discrete probability distributions. In this work, we go one step further and allow non-convexities as well as continuous distributions. We propose a gradient-based approach to approximately solve joint chance-constrained models. We approximate the original problem by smoothing indicator functions. Then, the smoothed chance constraints are relaxed by penalizing their violation in the objective function. The approximation problem is solved with the Continuous Stochastic Gradient method that is an enhanced version of the stochastic gradient descent and has recently been introduced in the literature. We present a convergence theory for the smoothing and penalty approximations. Under very mild assumptions, our approach is applicable to a wide range of chance-constrained optimization problems. As an example, we illustrate its computational efficiency on difficult practical problems arising in the operation of gas networks. The numerical experiments demonstrate that the approach quickly finds nearly feasible solutions for joint chance-constrained problems with non-convex constraint functions and continuous distributions, even for realistically-sized instances. Y1 - ER - TY - INPR A1 - Kuchlbauer, Martina T1 - Outer approximation for generalized convex mixed-integer nonlinear robust optimization problems N2 - We consider nonlinear robust optimization problems with mixed-integer decisions as well as nonconvexities. In detail, we consider cases where objective and constraint functions can be nonsmooth and generalized convex, i.e., f°-quasiconvex or f°-pseudoconvex. We propose an algorithm for such robust optimization problems that does not require a certain structure of the adversarial problem but only requires that approximate worst cases are available. As a result, our algorithm finds a robust optimal solution up to a tolerance. Our method integrates a bundle method into an outer approximation approach where the bundle method is used for the arising continuous subproblems. We rely on methods from the literature, namely a bundle method for nonlinear and nonconvex robust optimization problems and outer approximation approaches for quasiconvex settings. Our contribution is to combine them to one convergent robust optimization method that can cope with inexactness of worst-case evaluations. Further, we propose the gas transport under uncertainties as a relevant application and demonstrate that generalized convexity is fulfilled for a type of a network structure. Y1 - ER -