Natural Gas Maximal Load Delivery for Multi-contingency Analysis

TL;DR

This paper develops a convex relaxation approach to determine maximal load delivery in damaged natural gas pipeline networks, aiding risk assessment and operational decision-making during multi-contingency disruptions.

Contribution

It introduces a mixed-integer convex relaxation for the complex maximal load delivery problem, enabling efficient capacity evaluation of large-scale damaged gas networks.

Findings

Relaxation provides tight bounds on transport capacity.

Method converges quickly for networks up to 4197 junctions.

Effective in assessing impacts of multi-contingency disruptions.

Abstract

As the use of renewable generation has increased, electric power systems have become increasingly reliant on natural gas-fired power plants as fast ramping sources for meeting fluctuating bulk power demands. This dependence has introduced new vulnerabilities to the power grid, including disruptions to gas transmission networks from natural and man-made disasters. To address the operational challenges arising from these disruptions, we consider the task of determining a feasible steady-state operating point for a damaged gas pipeline network while ensuring the maximal delivery of load. We formulate the mixed-integer nonconvex maximal load delivery (MLD) problem, which proves difficult to solve on large-scale networks. To address this challenge, we present a mixed-integer convex relaxation of the MLD problem and use it to determine bounds on the transport capacity of a gas pipeline system. To demonstrate the effectiveness of the relaxation, the exact and relaxed formulations are compared across a large number of randomized damage scenarios on nine natural gas pipeline network models ranging in size from 11 to 4197 junctions. A proof of concept application, which assumes network damage from a set of synthetically generated earthquakes, is also presented to demonstrate the utility of the proposed optimization-based capacity evaluation in the context of risk assessment for natural disasters. For all but the largest network, the relaxation-based method is found to be suitable for use in evaluating the impacts of multi-contingency network disruptions, often converging to the optimal solution of the relaxed formulation in less than ten seconds.