[論文レビュー] Numerical Polynomial Homotopy Continuation Method to Locate All The Power Flow Solutions

Dhagash Mehta, Hung T. Nguyen+1 more|arXiv (Cornell University)|Aug 12, 2014

Power System Optimization and Stability被引用数 5

ひとこと要約

本稿では、複素ホモトピー続行法と$γ$-トリックを用いて分岐の問題を回避し、電力系統の潮流方程式のすべての解を体系的かつ網羅的に特定するための数値多項式ホモトピー続行法（NPHC）を提案する。この手法は数学的にすべての解を発見することを保証しており、14バスまでの系統で実証され、過渡安定性および電圧安定性解析にスケーラブルかつ並列処理可能なフレームワークを提供する。

ABSTRACT

The manuscript addresses the problem of finding all solutions of power flow equations or other similar nonlinear system of algebraic equations. This problem arises naturally in a number of power systems contexts, most importantly in the context of direct methods for transient stability analysis and voltage stability assessment. We introduce a novel form of homotopy continuation method called the numerical polynomial homotopy continuation (NPHC) method that is mathematically guaranteed to find all the solutions without ever encountering a bifurcation. The method is based on embedding the real form of power flow equation in complex space, and tracking the generally unphysical solutions with complex values of real and imaginary parts of the voltage. The solutions converge to physical real form in the end of the homotopy. The so-called $\\gamma$-trick mathematically rigorously ensures that all the paths are well-behaved along the paths, so unlike other continuation approaches, no special handling of bifurcations is necessary. The method is \ extit{embarrassingly parallelizable} and can be applied to reasonably large sized systems. We demonstrate the technique by analysis of several standard test cases up to the 14-bus system size. Finally, we discuss possible strategies for scaling the method to large size systems, and propose several applications for transient stability analysis and voltage stability assessment.

研究の動機と目的

• 非線形潮流方程式のすべての解を特定する課題に取り組み、過渡安定性および電圧安定性評価に不可欠である。
• 従来の手法がすべての解を特定できない、または分岐に敏感であるといった限界を克服するため、数学的に厳密なアプローチを導入する。
• 14バスまでの実際の電力系統サイズを扱えるスケーラブルかつ並列処理可能な手法を開発する。

提案手法

• NPHC法は実潮流方程式を複素空間に埋め込み、解の軌道を既知の解を持つ出発系から目的の系へ追跡する。
• 全次数ホモトピーと多面体ホモトピーを用いてホモトピー経路を構築し、後者は計算コストを顕著に低減する。
• $γ$-トリックにより、すべての解軌道が滑らかで特異でないことが保証され、分岐の特別な取り扱いが不要になる。
• この手法は、混合体積によるニュートン多面体の体積に基づく Bernstein-Khovanskii-Kushnirenko (BKK) 定理に依拠しており、孤立解の数の上限が保証される。
• MixedVolアルゴリズムを用いてBKK上限を効率的に計算するため、PHCpack や HOM4PS2 といった高度な数値ソルバーを活用する。
• この手法は本質的に並列処理可能であり、複数のプロセッサにわたり経路追跡計算を効率的に分散できる。

実験結果

リサーチクエスチョン

• RQ1 分岐や収束の問題によって解を逃すことがないよう、すべての潮流方程式の解を保証的に特定できるホモトピー続行法を設計できるか？
• RQ2 NPHC法は、従来の続行法やニュートン・ラプソン法と比較して、複数の潮流解を特定する際の効率性および正確性においてどのように差を示すか？
• RQ3 NPHC法はどの程度大きな電力系統にスケーリング可能であり、14バスを超える系統に拡張するにはどのような計算戦略が必要か？
• RQ4 NPHC法を用いて、直接的過渡安定性解析に必要な「最も近い不安定平衡点（CUEPs）」を体系的に計算できるか？
• RQ5 混合体積とニュートン多面体は、潮流方程式における孤立解の正確な数を推定・計算するために果たす役割は何か？

主な発見

• 標準的なテストケースにおいて、14バス系までにすべての解が正しく特定された。これは、現実の電力系統サイズで有効であることを示している。
• 多面体ホモトピーは全次数ホモトピーと比較して計算コストを顕著に低減し、手法のスケーラビリティを向上させた。
• $γ$-トリックにより、解軌道が滑らかで特異でないことが保証され、分岐に起因する失敗が回避されるため、数学的にすべての解が発見可能である。
• 混合体積を用いて計算されたBKK上限は、孤立解の数の正確な上限を提供しており、数値実験で検証された。
• この手法は本質的に並列処理可能であり、複数のプロセッサにわたり経路追跡計算を効率的に分散できる。
• 解の一般形はシステムパラメータにのみ依存するため、将来的にはパラメトリック安定性解析やロバストネス評価への応用が可能になる。