[论文解读] Natural Gas Flow Solutions with Guarantees: A Monotone Operator Theory Approach

Krishnamurthy Dvijotham, Marc Vuffray+2 more|arXiv (Cornell University)|Jun 19, 2015

Gas Dynamics and Kinetic Theory被引用 13

一句话总结

本文提出了一种单调算子理论方法，用于求解天然气流动方程，并在解的存在性、唯一性和计算方面提供可证明的保证。通过使用线性矩阵不等式（LMIs）识别单调性区域，该方法确保每个区域内至多存在一个解，并实现多项式时间内的计算，同时通过流量比的简化边界实现高效的区域选择。

ABSTRACT

We consider balanced flows in a natural gas transmission network and discuss computationally hard problems such as establishing if solution of the underlying nonlinear gas flow equations exists, if it is unique, and finding the solution. Particular topologies, e.g. trees, are known to be easy to solve based on a variational description of the gas flow equations, but these approaches do not generalize. In this paper, we show that the gas flow problem can be solved efficiently using the tools of monotone operator theory, provided that we look for solution within certain monotonicity domains. We characterize a family of monotonicity domains, described in terms of Linear Matrix Inequalities (LMI) in the state variables, each containing at most one solution. We also develop an efficient algorithm to choose a particular monotonicity domain, for which the LMI based condition simplifies to a bound on the flows. Performance of the technique is illustrated on exemplary gas networks.

研究动机与目标

• 解决大规模输气管网中非线性天然气流动方程的计算挑战。
• 建立天然气流动方程解的存在性与唯一性的可证明条件。
• 开发一种高效算法，用于在保证的单调性区域内计算解。
• 将现有可解情形（如树状网络、无环网络中的损失最小化）推广至任意网络拓扑。
• 提供一种框架，利用算子理论确保解的唯一性与多项式时间内的可计算性。

提出的方法

• 将天然气流动方程表述为关于流量和势变量的多元二次方程组。
• 通过单调算子理论将问题重新表述为变分不等式，其中解对应于单调算子的零点。
• 利用状态变量（压力和流量）的线性矩阵不等式（LMIs）定义单调性区域，确保每个区域内至多存在一个解。
• 通过矩阵值参数引入一族参数化的单调性区域，实现系统的区域选择。
• 对于特定参数选择，将LMI条件简化为网络中最大与最小流量之比的边界。
• 使用拟凸优化高效计算满足简化流量比条件的单调性区域参数。

实验结果

研究问题

• RQ1 天然气流动方程能否在保证存在性与唯一性的前提下高效求解？
• RQ2 何种条件可确保天然气流动算子为单调，从而支持多项式时间求解方法？
• RQ3 如何表征与计算单调性区域，以保证每个区域内至多一个解？
• RQ4 基于LMIs的单调性条件能否简化为适用于实际网络的可计算边界？
• RQ5 识别适合的单调性区域以实现高效解计算的计算复杂度是多少？

主要发现

• 在特定LMIs定义的区域内，天然气流动算子是单调的，确保每个区域内至多存在一个解。
• 在每个单调性区域内，天然气流动方程的解可实现多项式时间计算。
• 对于特定参数选择，LMI条件可简化为网络中最大与最小流量之比的边界。
• 通过求解拟凸优化问题，可高效计算单调性区域的参数。
• 该方法将已知的可解情形（如树状网络、无环网络中的损失最小化）推广至任意拓扑结构。
• 在测试网络上的数值评估验证了该方法的高效性与实用性。