2025-11-25T00:07:17.537013

Geometry of Distance Protection

Taylor, Domínguez-García
Distance relays detect faults on transmission lines. They face uncertainty from the fault's location and resistance, as well as the current from the line's remote terminal. In this paper, we aggregate this uncertainty with the Minkowski sum. This allows us to explicitly model the power grid surrounding the relay's line, and in turn accommodate any mix of synchronous machines and inverter-based resources. To make the relay's task easier, inverters can inject perturbations, or auxiliary signals, such as negative-sequence current. We use Farkas' lemma to construct an optimization for designing inverter auxiliary signals.
academic

Geometria della Protezione di Distanza

Informazioni Fondamentali

  • ID Articolo: 2510.04379
  • Titolo: Geometry of Distance Protection
  • Autori: Joshua A. Taylor (New Jersey Institute of Technology), Alejandro D. Domínguez-García (University of Illinois at Urbana-Champaign)
  • Classificazione: math.OC cs.IT cs.SY eess.SY math.IT
  • Data di Pubblicazione: Ottobre 2025 (preprint arXiv)
  • Link Articolo: https://arxiv.org/abs/2510.04379

Riassunto

I relè di distanza vengono utilizzati per rilevare i guasti nelle linee di trasmissione, ma affrontano incertezze relative alla posizione del guasto, alla resistenza di guasto e alla corrente del lato remoto. Questo articolo utilizza la somma di Minkowski per aggregare queste incertezze, consentendo di modellare esplicitamente la rete elettrica intorno al relè e adattarsi a combinazioni arbitrarie di risorse sincrone e invertitori. Per semplificare il compito del relè, gli invertitori possono iniettare segnali di disturbo o ausiliari (come corrente di sequenza negativa). Gli autori utilizzano il lemma di Farkas per costruire problemi di ottimizzazione per la progettazione di segnali ausiliari degli invertitori.

Contesto di Ricerca e Motivazione

Descrizione del Problema

  1. Sfide della protezione di distanza tradizionale: I relè di distanza rilevano i guasti misurando la tensione e la corrente locali per calcolare l'impedenza fino al punto di guasto, ma affrontano incertezze relative alla posizione del guasto, alla resistenza di guasto e alla corrente del lato remoto.
  2. Nuovi problemi introdotti dalle risorse invertitori: Con l'aumento della percentuale di risorse basate su invertitori (IBR) nella rete, le caratteristiche tradizionali di protezione di distanza basate su ipotesi di generatori sincroni potrebbero presentare malfunzionamenti.
  3. Limitazioni dei metodi esistenti:
    • Le caratteristiche quadrilatere tradizionali si basano su aggiustamenti temporanei ispirati dalla fisica per gestire la corrente del lato remoto
    • Questi aggiustamenti presuppongono fonti sincrone, che potrebbero fallire in reti dominate da IBR
    • Manca un approccio sistematico per progettare segnali ausiliari IBR

Motivazione della Ricerca

Gli autori mirano a formalizzare l'aggregazione delle incertezze nella protezione di distanza attraverso metodi geometrici (somma di Minkowski) e sviluppare metodi di ottimizzazione per progettare segnali ausiliari IBR, al fine di migliorare l'affidabilità della protezione di distanza nelle reti moderne.

Contributi Principali

  1. Contributi Teorici: Primo utilizzo della somma di Minkowski per formalizzare l'aggregazione delle incertezze nelle caratteristiche dei relè di distanza, fornendo una base matematica per le caratteristiche quadrilatere tradizionali.
  2. Innovazione nella Modellazione: Propone un framework in grado di modellare esplicitamente la rete e le fonti di energia, adattandosi a combinazioni arbitrarie di macchine sincrone e IBR.
  3. Metodo di Ottimizzazione: Costruisce problemi di ottimizzazione per la progettazione di segnali ausiliari IBR basati sul lemma di Farkas, garantendo che il relè possa distinguere chiaramente diversi scenari di guasto.
  4. Metodi Computazionali: Sviluppa metodi computazionali efficienti, inclusi algoritmi di proiezione zonotope e rilassamenti zonogon, consentendo il calcolo online delle nuove caratteristiche.

Dettagli del Metodo

Definizione del Compito

Input: Topologia della rete, incertezze delle fonti, intervalli di parametri di guasto Output: Regione caratteristica del relè, segnale ausiliario ottimale Obiettivo: Minimizzare l'ampiezza del segnale ausiliario garantendo la separabilità di tutti gli scenari di guasto

Modello Fisico

Modellazione della Rete

  • Modellazione della macchina sincrona come fonte di tensione: vS=vS∘v_S = v_S^{\circ}
  • Modellazione dell'IBR come fonte di corrente: iC=iC∘+Δi_C = i_C^{\circ} + \Delta (incluso il segnale ausiliario)
  • Carico modellato con impedenza o circuiti equivalenti RLC

Analisi del Circuito di Guasto

Per diversi tipi di guasto, stabilire la legge di Kirchhoff sulla tensione (KVL):

Guasto a terra (ag): vLa=mzz(iLa+kiL0)+mrrF(iLa+iRa)v_L^a = m_z z(i_L^a + ki_L^0) + m_r r_F(i_L^a + i_R^a)

Guasto tra fasi (ab): vLa−vLb=mzz(iLa−iLb)+mrrF2(iLa−iLb+iRa−iRb)v_L^a - v_L^b = m_z z(i_L^a - i_L^b) + \frac{m_r r_F}{2}(i_L^a - i_L^b + i_R^a - i_R^b)

dove mz∈[0,1]m_z \in [0,1] è la distanza di guasto normalizzata e mr∈[0,1]m_r \in [0,1] è la resistenza di guasto normalizzata.

Modellazione dell'Incertezza

Incertezza della Fonte

Definire l'insieme di incertezza: U={[iC∘vS∘]+Σλ∣λ∈Λ}U = \left\{\begin{bmatrix} i_C^{\circ} \\ v_S^{\circ} \end{bmatrix} + \Sigma\lambda \mid \lambda \in \Lambda\right\}

dove Λ\Lambda può essere:

  • Limitato in norma: Λ={λ∣∥λ∥2≤1}\Lambda = \{\lambda \mid \|\lambda\|_2 \leq 1\}
  • Poliedrico (zonotope): Λ={λ∣Pλ≤1}\Lambda = \{\lambda \mid P\lambda \leq 1\}

Aggregazione mediante Somma di Minkowski

Per lo scenario di guasto η\eta, l'insieme di incertezza dell'impedenza apparente è: Z~Aη=ZT⊕Z~Fη\tilde{Z}_A^{\eta} = Z_T \oplus \tilde{Z}_F^{\eta}

dove ZT={mzz∣mz∈[mz‾,1]}Z_T = \{m_z z \mid m_z \in [\underline{m_z}, 1]\} è l'insieme di impedenza della linea e Z~Fη\tilde{Z}_F^{\eta} è l'insieme di impedenza di guasto.

Incertezza Prima e Dopo il Test

Incertezza Dopo il Test (Caratteristica del Relè)

Dopo il test, la corrente locale iLi_L è nota, l'insieme di incertezza diventa un insieme convesso, calcolabile efficientemente mediante proiezione nel piano complesso.

Incertezza Prima del Test (Progettazione del Segnale Ausiliario)

Prima del test, iLi_L è sconosciuta, contenendo termini bilineari mzλm_z\lambda e mrλm_r\lambda, richiedendo rilassamento convesso:

Metodo di Rilassamento 1: UmRel={m[iC∘vS∘]+Σλ∣λ∈Λ,m∈[m‾,1]}U_m^{\text{Rel}} = \left\{m\begin{bmatrix} i_C^{\circ} \\ v_S^{\circ} \end{bmatrix} + \Sigma\lambda \mid \lambda \in \Lambda, m \in [\underline{m}, 1]\right\}

Metodo di Rilassamento 2: Separare le variabili λmz\lambda_{m_z} e λmr\lambda_{m_r}, eliminando i vincoli di accoppiamento.

Ottimizzazione del Segnale Ausiliario

Condizione di Separazione

Definire una coppia di modelli (η1,η2)(η_1, η_2) separata dal segnale ausiliario Δ\Delta se e solo se: P˘η1,η2(Δ)=∅\breve{P}^{η_1,η_2}(\Delta) = \emptyset

ovvero l'intersezione dei due insiemi di incertezza prima del test è vuota.

Problema di Ottimizzazione

\min_{\Delta} &\quad \Delta^{\dagger}Q\Delta \\ \text{s.t.} &\quad \breve{P}^{η_1,η_2}(\Delta) = \emptyset, \quad \forall (η_1,η_2) \in \mathcal{F}^2 \end{align}$$ #### Applicazione del Lemma di Farkas Utilizzare il lemma di Farkas per convertire il vincolo di insieme vuoto in fattibilità del sistema duale: $$\breve{P}^{η_1,η_2}(\Delta) = \emptyset \Leftrightarrow \mathcal{S}^{η_1,η_2}(\Delta) \neq \emptyset$$ Il problema di ottimizzazione finale diventa una programmazione biconvessa, risolvibile con l'algoritmo ADMM. ## Configurazione Sperimentale ### Sistema di Test - **Rete di Base**: Sistema IEEE a 14 nodi - **Configurazione delle Fonti**: - Nodi con macchine sincrone: $\mathcal{S} = \{1,3,5,14\}$ - Nodi IBR: $\mathcal{C} = \{2,4,7,8,12\}$ - Nodi di carico: $\{6,10,11,13\}$ ### Impostazione dei Parametri - Macchine sincrone: Fonti di tensione con ampiezza unitaria, angoli uniformemente distribuiti in $[0, 2π/10]$ - IBR: Fonti di corrente con ampiezza unitaria, angoli uniformemente distribuiti in $[0, 2π/3]$ - Rumore: Aggiungere $\lambda_k[1, α, α^2]^T$ a ogni fonte, $λ_k \in \Lambda_k$ (poligono regolare) - Parametri di guasto: $k=2$ (coefficiente di compensazione di sequenza zero), $r_F=1$, $\underline{m_z}=0.15$ ### Implementazione Computazionale - **Linguaggio**: Python - **Visualizzazione**: Matplotlib - **Ottimizzatore**: CVXPy + risolutore Clarabel - **Tempi di Calcolo**: - $\tilde{Z}_A^{\eta}$ per tutti gli 11 tipi di guasto: circa 0,9 secondi - Rilassamento Zonogon $\hat{Z}_A^{\eta}$: circa 1 millisecondo ## Risultati Sperimentali ### Insiemi di Incertezza Dopo il Test La Figura 4 mostra gli insiemi di incertezza precisi $\tilde{Z}_A^{\eta}$ e i loro rilassamenti zonogon $\hat{Z}_A^{\eta}$ per guasti ag e ab. I risultati mostrano: - Il rilassamento approssima bene l'insieme preciso - L'aggiunta di vincoli di taglio $\text{Im}[z]\text{Re}[z_A^{\eta}] \geq \text{Re}[z]\text{Im}[z_A^{\eta}]$ può stringere il rilassamento ### Risultati della Progettazione del Segnale Ausiliario #### Analisi Visiva La Figura 5 mostra gli insiemi di segnali ausiliari che separano diverse coppie di guasti: - Segnale ausiliario minimo per separare (N,ag), (N,ab) e (ag,ab): $δ = j0.2$, $\|\Delta\|_2 ≈ 0.53$ - Segnale ausiliario minimo per separare tutti i guasti LG e LL: $δ = 1.6+j0.6$, $\|\Delta\|_2 ≈ 4.52$ #### Risultati dell'Ottimizzazione Utilizzando l'algoritmo ADMM si ottiene: $$\Delta ≈ [1.55+j0.57, 1.56+j0.56, 1.59+j0.50, 1.50+j0.64, 1.50+j0.63, 1.52+j0.60, 1.48+j0.66]^T$$ L'ampiezza finale $\|\Delta\|_2 ≈ 4.35$, con un miglioramento del 4% rispetto ai risultati della ricerca su griglia. ### Confronto con Metodi Tradizionali - **Caratteristiche Quadrilatere**: Basate su aggiustamenti temporanei, presuppongono dominanza sincrona - **Metodo Proposto**: Somma di Minkowski precisa, adatta a combinazioni arbitrarie di fonti - **Efficienza Computazionale**: Il rilassamento zonogon può essere completato in 1 millisecondo, supportando il calcolo online ## Lavori Correlati ### Fondamenti della Protezione di Distanza - La protezione di distanza tradizionale si basa sulla misurazione dell'impedenza e sulle caratteristiche quadrilatere - Affronta problemi di malfunzionamento in reti dominate da IBR ### Impatto degli IBR sulla Protezione - Le caratteristiche di controllo degli IBR differiscono dalle macchine sincrone, causando cambiamenti nelle caratteristiche della corrente di guasto - La ricerca esistente si concentra principalmente su aggiustamenti della strategia di controllo ### Metodi di Segnali Ausiliari - Metodi di iniezione armonica, iniezione di corrente di sequenza negativa, ecc. - Lo Standard IEEE 2800 specifica che gli IBR devono iniettare corrente di sequenza negativa in anticipo di 90° ### Applicazione di Metodi Geometrici nella Protezione - Primo utilizzo della somma di Minkowski nella protezione di distanza - Fornisce un framework sistematico per l'aggregazione delle incertezze ## Conclusioni e Discussione ### Conclusioni Principali 1. **Contributi Teorici**: La somma di Minkowski fornisce un framework matematico rigoroso per l'aggregazione delle incertezze nella protezione di distanza 2. **Valore Pratico**: Il nuovo metodo si adatta a combinazioni complesse di fonti nelle reti moderne 3. **Fattibilità Computazionale**: Il rilassamento zonogon rende il calcolo online possibile 4. **Efficacia dell'Ottimizzazione**: Il metodo di progettazione del segnale ausiliario basato sul lemma di Farkas garantisce la separabilità dei guasti ### Limitazioni 1. **Semplificazioni del Modello**: - La modellazione delle fonti è relativamente grossolana (fonte di tensione/corrente) - La matrice di rete viene calcolata utilizzando parametri di guasto nominali - Vengono ignorati i processi dinamici 2. **Complessità Computazionale**: - La natura bilineare dell'insieme di incertezza prima del test richiede rilassamento - La convergenza dell'algoritmo ADMM dipende dall'inizializzazione 3. **Considerazioni Pratiche**: - Richiede modelli di rete e parametri di fonte accurati - L'implementazione pratica dei segnali ausiliari potrebbe affrontare vincoli ### Direzioni Future 1. **Verifica mediante Simulazione Elettromagnetica Transitoria**: Verificare l'efficacia del metodo in simulazioni dettagliate di sistemi di potenza 2. **Miglioramento della Matrice di Rete**: Considerare parametri di guasto nominali multipli o aggiornamento adattivo 3. **Coordinamento Multi-Relè**: Estendere a protezione coordinata di più linee 4. **Altri Schemi di Protezione**: Applicare le idee fondamentali a protezione di sovracorrente, protezione nel dominio del tempo, ecc. ## Valutazione Approfondita ### Punti di Forza 1. **Forte Innovazione Teorica**: - Primo utilizzo della somma di Minkowski nel campo della protezione di distanza - Fornisce un framework matematico rigoroso in sostituzione dei metodi temporanei - L'applicazione del lemma di Farkas è ingegnosa ed efficace 2. **Alto Valore Pratico**: - Risolve problemi reali nella protezione di reti dominate da IBR - I metodi computazionali sono efficienti e supportano applicazioni online - La progettazione del segnale ausiliario è sistematica e robusta 3. **Completezza del Metodo**: - Analisi completa che copre 11 tipi di guasto - Fornisce molteplici metodi di rilassamento convesso - Implementazione algoritmica dettagliata e open source ### Carenze 1. **Limitazioni delle Ipotesi del Modello**: - Il modello delle fonti è eccessivamente semplificato, potrebbe influenzare l'applicazione pratica - Ignora le caratteristiche dinamiche della rete - Il trattamento dei circuiti intermedi non è sufficientemente approfondito 2. **Verifica Insufficiente**: - Mancanza di verifica con dati di reti reali - Nessun confronto con relè commerciali - La verifica mediante simulazione elettromagnetica transitoria è ancora in sospeso 3. **Complessità Computazionale**: - Difficile garantire l'ottimalità globale per problemi di ottimizzazione bilineare - La scalabilità computazionale per reti di grandi dimensioni non è stata sufficientemente verificata ### Valutazione dell'Impatto 1. **Impatto Accademico**: - Introduce nuovi strumenti matematici nel campo della protezione dei sistemi di potenza - Potrebbe ispirare più applicazioni di metodi geometrici nella protezione - Ha importante valore di riferimento per la ricerca sulla protezione IBR 2. **Valore Ingegneristico**: - Fornisce una base teorica per l'aggiornamento dei dispositivi di protezione nelle reti moderne - Il metodo di progettazione del segnale ausiliario ha potenziale di applicazione pratica - Aiuta a formulare standard e normative correlate ### Scenari Applicabili 1. **Reti con Alta Penetrazione di IBR**: Particolarmente adatto a sistemi con elevata percentuale di energie rinnovabili 2. **Nuove Linee di Trasmissione**: Può considerare nuove strategie di protezione nella fase di progettazione 3. **Aggiornamento dei Dispositivi di Protezione**: Fornisce soluzioni per l'intelligentizzazione dei relè esistenti 4. **Pianificazione della Rete**: Considera l'ottimizzazione delle prestazioni di protezione nella pianificazione della rete ## Bibliografia L'articolo cita 26 importanti riferimenti, principalmente includenti: - Testi classici sulla protezione di distanza: Ziegler (2011), Horowitz et al. (2022) - Ricerche recenti sull'impatto degli IBR sulla protezione: Kasztenny (2022), Baeckeland et al. (2022) - Lavori correlati ai segnali ausiliari: Saleh et al. (2020), Banaiemoqadam et al. (2019) - Strumenti matematici correlati: Teoria della somma di Minkowski, lemma di Farkas, teoria dei zonotopi --- **Valutazione Complessiva**: Questo è un articolo accademico di alta qualità con importanti contributi sia in termini di innovazione teorica che di valore pratico. Gli autori hanno abilmente introdotto metodi geometrici nel tradizionale campo della protezione dei sistemi di potenza, fornendo nuove prospettive per affrontare le sfide di protezione affrontate dalle reti moderne. Sebbene vi sia ancora spazio per miglioramenti nelle ipotesi del modello e nella verifica, il concetto fondamentale e il framework metodologico possiedono importante valore accademico e prospettive di applicazione ingegneristica.