2025-11-18T17:07:13.972479

An Augmented Lagrangian Method on GPU for Security-Constrained AC Optimal Power Flow

Pacaud, NurkanoviÄ, Pozharskiy et al.

We present a new algorithm for solving large-scale security-constrained optimal power flow in polar form (AC-SCOPF). The method builds on Nonlinearly Constrained augmented Lagrangian (NCL), an augmented Lagrangian method in which the subproblems are solved using an interior-point method. NCL has two key advantages for large-scale SC-OPF. First, NCL handles difficult problems such as infeasible ones or models with complementarity constraints. Second, the augmented Lagrangian term naturally regularizes the Newton linear systems within the interior-point method, enabling to solve the Newton systems with a pivoting-free factorization that can be efficiently parallelized on GPUs. We assess the performance of our implementation, called MadNCL, on large-scale corrective AC-SCOPFs, with complementarity constraints modeling the corrective actions. Numerical results show that MadNCL can solve AC-SCOPF with 500 buses and 256 contingencies fully on the GPU in less than 3 minutes, whereas Knitro takes more than 3 hours to find an equivalent solution.

academic

GPU पर सुरक्षा-विवश AC इष्टतम शक्ति प्रवाह के लिए एक संवर्धित लैग्रेंजियन विधि

मूल जानकारी

पेपर ID: 2510.13333
शीर्षक: An Augmented Lagrangian Method on GPU for Security-Constrained AC Optimal Power Flow
लेखक: François Pacaud, Armin Nurkanović, Anton Pozharskiy, Alexis Montoison, Sungho Shin
वर्गीकरण: math.OC (अनुकूलन और नियंत्रण)
प्रकाशन समय: 15 अक्टूबर 2025
पेपर लिंक: https://arxiv.org/abs/2510.13333

सारांश

यह पेपर बड़े पैमाने पर सुरक्षा-विवश AC इष्टतम शक्ति प्रवाह (AC-SCOPF) को हल करने के लिए एक नई एल्गोरिदम प्रस्तावित करता है। यह विधि अरैखिक विवश संवर्धित लैग्रेंजियन (NCL) विधि पर आधारित है, जो उप-समस्याओं को हल करने के लिए आंतरिक बिंदु विधि का उपयोग करती है। NCL के बड़े पैमाने पर SC-OPF के लिए दो मुख्य लाभ हैं: पहला, NCL कठिन समस्याओं जैसे अक्षम समस्याओं या पूरक विवशताओं वाले मॉडल को संभाल सकता है; दूसरा, संवर्धित लैग्रेंजियन पद आंतरिक बिंदु विधि में न्यूटन रैखिक प्रणाली को स्वाभाविक रूप से नियमित करता है, जिससे पिवट-मुक्त अपघटन का उपयोग करके न्यूटन प्रणाली को हल किया जा सकता है, जिसे GPU पर कुशलतापूर्वक समानांतर किया जा सकता है। संख्यात्मक परिणाम दिखाते हैं कि MadNCL 500 नोड्स और 256 विफलताओं के साथ AC-SCOPF को 3 मिनट से कम समय में GPU पर पूरी तरह से हल कर सकता है, जबकि Knitro को समकक्ष समाधान खोजने में 3 घंटे से अधिक समय लगता है।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या विवरण

विद्युत संचरण नेटवर्क में, इष्टतम प्रेषण आमतौर पर सुरक्षा-विवश इष्टतम शक्ति प्रवाह (SCOPF) को हल करके गणना की जाती है। यह प्रेषण भौतिक विवशताओं (शक्ति प्रवाह, लाइन प्रवाह सीमाएं) और जनरेटर क्षमता पर विचार करते हुए दिए गए मानदंड (लागत या नेटवर्क हानि) को कम करता है। इसके अलावा, प्रेषण लाइन या जनरेटर विफलताओं के अनुरूप आपातकालीन स्थितियों की एक श्रृंखला में व्यवहार्य रहना चाहिए (N-1 सुरक्षा मानदंड)।

समस्या का महत्व

SCOPF आमतौर पर DC-SCOPF नामक बड़े पैमाने की रैखिक प्रोग्रामिंग समस्या के रूप में तैयार किया जाता है, जिसका आकार विफलताओं की संख्या के साथ रैखिक रूप से बढ़ता है। लेकिन यह अरैखिक भौतिक विवशताओं के रैखिकीकरण की कीमत पर आता है, जो समाधान की सटीकता को प्रभावित करता है। हालांकि, मूल अरैखिक सूत्रीकरण के साथ AC-SCOPF को हल करना अभी भी एक खुली चुनौती है।

मौजूदा विधियों की सीमाएं

अरैखिक सूत्रीकरण में दो समस्याएं हैं:

AC-SCOPF एक विशाल पैमाने की अरैखिक प्रोग्रामिंग समस्या है, जिसका आकार Ipopt या Knitro जैसे अत्याधुनिक अरैखिक अनुकूलन सॉल्वर की क्षमता से अधिक है
AC-SCOPF मॉडल में पूरक विवशताएं संख्यात्मक रूप से संभालने में कठिन हैं, विशेष एल्गोरिदम को अपनाया जाना चाहिए

अनुसंधान प्रेरणा

बड़े पैमाने पर AC-SCOPF की विशेषताएं एल्गोरिदम को सीमा तक धकेलती हैं, क्योंकि पूरक विवशताओं की संख्या विफलताओं की संख्या के साथ रैखिक रूप से बढ़ती है। इस चुनौती का सामना करने के लिए, लेखकों ने AC-SCOPF को हल करने के लिए अरैखिक विवश संवर्धित लैग्रेंजियन (NCL) पर आधारित संवर्धित लैग्रेंजियन विधि का उपयोग करने का प्रस्ताव दिया है।

मुख्य योगदान

संवर्धित लैग्रेंजियन विधि का पहला अनुप्रयोग: पहली बार पूरक विवशताओं के साथ सुधारात्मक SCOPF को हल करने के लिए संवर्धित लैग्रेंजियन-आधारित एल्गोरिदम का उपयोग किया गया है
GPU त्वरण कार्यान्वयन: MadNCL विकसित किया गया, जो GPU त्वरण का समर्थन करने वाला NCL कार्यान्वयन है, जो रैखिक समाधान के लिए NVIDIA cuDSS का उपयोग करता है
पूरक विवशताओं को संभालना: MadNCL AC-SCOPF में पूरक विवशताओं को अच्छी तरह से संभाल सकता है और अक्षम समस्याओं का प्रभावी ढंग से पता लगा सकता है
महत्वपूर्ण प्रदर्शन सुधार: बड़े पैमाने के उदाहरणों पर पारंपरिक विधियों की तुलना में महत्वपूर्ण त्वरण प्राप्त किया गया है, GPU संस्करण CPU संस्करण से 6 गुना से अधिक तेज है

विधि विवरण

कार्य परिभाषा

सुरक्षा-विवश AC इष्टतम शक्ति प्रवाह (AC-SCOPF) समस्या को इस प्रकार परिभाषित किया गया है:

$\min_{x,u} f(x^0, u^0)$

विवशताओं के अधीन:

आधार स्थिति: $g^0(x^0, u^0) = 0$ , $h^0(x^0, u^0) \leq 0$
प्रत्येक विफलता $k \in \{1, \cdots, K\}$ $k \in {1, \dots, K}$ के लिए:
- $g^k(u^0, x^k, u^k) = 0$
- $h^k(x^k, u^k) \leq 0$
- $0 \leq G(x^k, u^k) \perp H(x^k, u^k) \geq 0$

जहां पूरक विवशताएं निम्नलिखित से आती हैं:

स्वचालित पीढ़ी नियंत्रण (AGC): सक्रिय शक्ति समायोजन आवृत्ति को नियंत्रित करता है
PV/PQ स्विचिंग: वोल्टेज नियंत्रण और प्रतिक्रियाशील शक्ति सीमाएं

मॉडल आर्किटेक्चर

MPCC पुनर्निर्माण

AC-SCOPF को गणितीय प्रोग्रामिंग और पूरक विवशताओं (MPCC) में पुनर्निर्मित किया गया है:

$\min_{w \in \mathbb{R}^n} \phi(w) \quad \text{s.t.} \quad \begin{cases} c(w) = 0, \quad w_0 \geq 0 \\ 0 \leq w_1 \perp w_2 \geq 0 \end{cases}$

NCL एल्गोरिदम

NCL दो स्तरों पर काम करता है:

बाहरी पुनरावृत्ति: दंड पैरामीटर $\rho^{(n)}$ और गुणक अनुमान $(λ^{(n)}, ν_0^{(n)})$ को अपडेट करता है
आंतरिक पुनरावृत्ति: विवश अरैखिक उप-समस्या को हल करता है:

$\min_{w,r,t} L_ρ(w, r, t, λ^{(n)}, ν_0^{(n)})$

विवशताओं के अधीन: $c(w) - r = 0, \quad W_1W_2e \leq t, \quad (w_0, w_1, w_2) \geq 0$

न्यूटन प्रणाली संरचना

उप-समस्या की न्यूटन प्रणाली निम्नलिखित संरचना रखती है:

$\begin{bmatrix} A & B^⊤ \\ B & -C \end{bmatrix} \begin{bmatrix} Δw \\ Δy \end{bmatrix} = \begin{bmatrix} r_1 \\ r_2 \end{bmatrix}$

जहां संवर्धित लैग्रेंजियन पद द्वारा प्रदान किया गया नियमितकरण पिवट-मुक्त अपघटन का उपयोग करने की अनुमति देता है।

तकनीकी नवाचार

प्राकृतिक नियमितकरण: संवर्धित लैग्रेंजियन पद न्यूटन रैखिक प्रणाली को स्वाभाविक रूप से नियमित करता है, यहां तक कि जब सख्त पूरकता विफल हो तब भी प्रणाली की गैर-विलक्षणता को बनाए रखता है
पिवट-मुक्त अपघटन: नियमितकरण प्रतीकात्मक Cholesky अपघटन जैसी पिवट-मुक्त विधियों का उपयोग करने की अनुमति देता है, जिन्हें GPU पर कुशलतापूर्वक समानांतर किया जा सकता है
अक्षमता का पता लगाना: जब समस्या अक्षम हो, तो NCL स्वचालित रूप से व्यवहार्यता समस्या में वापस आता है, दंड पैरामीटर $ρ^{(n)}$ को अनंत तक बढ़ाकर

प्रायोगिक सेटअप

डेटासेट

MATPOWER लाइब्रेरी से उदाहरण का उपयोग किया गया:

118ieee, ACTIVSg200, 300ieee, ACTIVSg500
1354pegase, ACTIVSg2000, 2869pegase
विफलताओं की संख्या 2 से 256 तक भिन्न है

मूल्यांकन मेट्रिक्स

समाधान समय: कुल समाधान समय और प्रति पुनरावृत्ति समय
पुनरावृत्ति संख्या: आंतरिक बिंदु विधि पुनरावृत्ति संख्या
उद्देश्य मान: इष्टतम समाधान का उद्देश्य फ़ंक्शन मान
व्यवहार्यता: अक्षम विफलताओं का पता लगाने की क्षमता

तुलना विधियां

Knitro: MPCC का समर्थन करने वाला अत्याधुनिक अनुकूलन सॉल्वर, $\ell_1$ सटीक दंड विधि का उपयोग करता है
MadNCL-CPU: HSL MA57 का उपयोग करने वाला CPU संस्करण
MadNCL-GPU: NVIDIA cuDSS का उपयोग करने वाला GPU संस्करण

कार्यान्वयन विवरण

प्रोग्रामिंग भाषा: Julia 1.11
अभिसरण सहनशीलता: 1e-6
न्यूनतम बाधा पैरामीटर: $μ_{min} = 10^{-7}$
हार्डवेयर: AMD EPYC 7430 प्रोसेसर, NVIDIA A30 GPU (24GB मेमोरी)

प्रायोगिक परिणाम

मुख्य परिणाम

विफलता स्क्रीनिंग प्रदर्शन

विफलता स्क्रीनिंग कार्य में, MadNCL Knitro से काफी बेहतर है:

उदाहरण	Knitro (s)	MadNCL-CPU (s)
118ieee	0.5	0.01
ACTIVSg500	5.4	0.3
2869pegase	238.4	14.1

MadNCL 300 नोड्स से अधिक के उदाहरणों पर कम से कम 10 गुना तेज है।

बड़े पैमाने पर AC-SCOPF समाधान

ACTIVSg500 उदाहरण के लिए, विफलताओं की संख्या बढ़ने के साथ:

K	चर संख्या	Knitro समय(s)	MadNCL-GPU समय(s)	त्वरण अनुपात
64	241,900	2159.59	27.96	77.2×
128	480,300	4852.33	46.40	104.6×
256	957,100	11136.08	170.75	65.2×

विलोपन प्रयोग

GPU बनाम CPU प्रदर्शन

MadNCL-GPU की तुलना में MadNCL-CPU का प्रदर्शन सुधार:

K≥64 के लिए, GPU संस्करण CPU संस्करण से लगभग 6 गुना तेज है
K≥64 के लिए, GPU संस्करण Knitro से लगभग 20 गुना से अधिक तेज है

प्रति पुनरावृत्ति समय विश्लेषण

विफलताओं की संख्या बढ़ने के साथ, MadNCL-GPU के प्रति पुनरावृत्ति समय में सबसे कम वृद्धि होती है, जो अच्छी स्केलेबिलिटी प्रदर्शित करता है।

प्रायोगिक निष्कर्ष

स्केलेबिलिटी: MadNCL उत्कृष्ट स्केलेबिलिटी प्रदर्शित करता है, लगभग दस लाख चर की समस्याओं को संभाल सकता है
मजबूती: NCL स्वचालित रूप से अक्षम समस्याओं का पता लगा सकता है और उन्हें संभाल सकता है
समानांतर दक्षता: GPU कार्यान्वयन समानांतर कंप्यूटिंग के लाभों का पूरी तरह से उपयोग करता है
संख्यात्मक स्थिरता: संवर्धित लैग्रेंजियन का नियमितकरण संख्यात्मक स्थिरता में सुधार करता है

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

MadNCL बड़े पैमाने पर AC-SCOPF समस्याओं को प्रभावी ढंग से हल कर सकता है, लगभग दस लाख चर को संभाल सकता है
GPU त्वरण संस्करण पारंपरिक CPU सॉल्वर की तुलना में दसियों गुना प्रदर्शन सुधार प्राप्त करता है
संवर्धित लैग्रेंजियन विधि पूरक विवशताओं को संभालने के लिए एक मजबूत समाधान प्रदान करती है

सीमाएं

शर्त संख्या समस्या: समस्या के आकार के बढ़ने के साथ, रैखिक प्रणाली की शर्त संख्या बिगड़ जाती है
अभिसरण: कुछ बड़े पैमाने के उदाहरणों पर अभिसरण पर्याप्त रूप से स्थिर नहीं है
मेमोरी सीमाएं: GPU मेमोरी सीमाएं संभाली जा सकने वाली अधिकतम समस्या के आकार को सीमित करती हैं

भविष्य की दिशाएं

आंतरिक बिंदु विधि न्यूटन प्रणाली में खराब शर्त संख्या समस्या को हल करना
बड़े पैमाने के उदाहरणों तक विस्तार (हजारों नोड्स, सैकड़ों विफलताएं)
संख्यात्मक स्थिरता में सुधार के लिए पूर्व-शर्त तकनीकों में सुधार

गहन मूल्यांकन

लाभ

विधि नवाचार: पहली बार NCL को AC-SCOPF पर लागू किया गया, तकनीकी मार्ग नवीन है
कार्यान्वयन गुणवत्ता: उच्च गुणवत्ता का GPU कार्यान्वयन, समानांतर कंप्यूटिंग के लाभों का पूरी तरह से उपयोग करता है
व्यापक प्रयोग: व्यापक प्रायोगिक मूल्यांकन, स्केलेबिलिटी और मजबूती परीक्षण सहित
व्यावहारिक मूल्य: महत्वपूर्ण प्रदर्शन सुधार बड़े पैमाने पर वास्तविक समय अनुप्रयोगों को संभव बनाता है

कमियां

सैद्धांतिक विश्लेषण: SCOPF समस्याओं पर NCL के अभिसरण सिद्धांत विश्लेषण का अभाव
संख्यात्मक स्थिरता: अधिकतम पैमाने के उदाहरणों पर अभी भी संख्यात्मक स्थिरता समस्याएं मौजूद हैं
सामान्यता: विधि की प्रयोज्यता मुख्य रूप से विद्युत प्रणाली अनुकूलन क्षेत्र तक सीमित है

प्रभाव

शैक्षणिक योगदान: बड़े पैमाने पर गैर-उत्तल अनुकूलन के लिए नई समाधान सोच प्रदान करता है
व्यावहारिक मूल्य: विद्युत प्रणाली संचालन और योजना के लिए महत्वपूर्ण व्यावहारिक मूल्य है
तकनीकी प्रचार: GPU त्वरण अनुकूलन एल्गोरिदम का सफल केस स्टडी

लागू परिस्थितियां

विद्युत प्रणाली प्रेषण: वास्तविक समय और दिन-पहले बाजार की सुरक्षा-विवश अनुकूलन
बड़े पैमाने पर गैर-उत्तल अनुकूलन: पूरक विवशताओं वाली अन्य इंजीनियरिंग अनुकूलन समस्याएं
GPU उच्च-प्रदर्शन कंप्यूटिंग: तेजी से समाधान की आवश्यकता वाले अनुकूलन अनुप्रयोग

संदर्भ

पेपर 31 संबंधित संदर्भों का हवाला देता है, जो SCOPF मॉडलिंग, MPCC समाधान विधियों, संवर्धित लैग्रेंजियन सिद्धांत और GPU अनुकूलन सहित कई पहलुओं के महत्वपूर्ण कार्यों को शामिल करता है, जो अनुसंधान के लिए एक ठोस सैद्धांतिक आधार प्रदान करता है।