Relaxations of the Steady Optimal Gas Flow Problem for a Non-Ideal Gas

TL;DR

This paper develops tight relaxations for the steady optimal gas flow problem in pipeline networks, enabling fast and near-optimal solutions for non-ideal gas models, significantly improving operational efficiency.

Contribution

It introduces novel polyhedral relaxations for the nonlinear constraints of the OGF problem, scalable to large networks and capable of proving optimality in most cases.

Findings

Relaxations solve within seconds on standard laptops.

Proven optimality in 92% of tested instances.

Effective for large-scale realistic networks.

Abstract

Natural gas ranks second in consumption among primary energy sources in the United States. The majority of production sites are in remote locations, hence natural gas needs to be transported through a pipeline network equipped with a variety of physical components such as compressors, valves, etc. Thus, from the point of view of both economics and reliability, it is desirable to achieve optimal transportation of natural gas using these pipeline networks. The physics that governs the flow of natural gas through various components in a pipeline network is governed by nonlinear and non-convex equality and inequality constraints and the most general steady-flow operations problem takes the form of a Mixed Integer Nonlinear Program. In this paper, we consider one example of steady-flow operations -- the Optimal Gas Flow (OGF) problem for a natural gas pipeline network that minimizes the production cost subject to the physics of steady-flow of natural gas. The ability to quickly determine global optimal solution and a lower bound to the objective value of the OGF for different demand profiles plays a key role in efficient day-to-day operations. One strategy to accomplish this relies on tight relaxations to the nonlinear constraints of the OGF. Currently, many nonlinear constraints that arise due to modeling the non-ideal equation of state either do not have relaxations or have relaxations that scale poorly for realistic network sizes. In this work, we combine recent advancements in the development of polyhedral relaxations for univariate functions to obtain tight relaxations that can be solved within a few seconds on a standard laptop. We demonstrate the quality of these relaxations through extensive numerical experiments on very large scale test networks available in the literature and find that the proposed relaxation is able to prove optimality in 92% of the instances.