Geometry of Distance Protection

基本信息

• 论文ID: 2510.04379
• 标题: Geometry of Distance Protection
• 作者: Joshua A. Taylor (New Jersey Institute of Technology), Alejandro D. Domínguez-García (University of Illinois at Urbana-Champaign)
• 分类: math.OC cs.IT cs.SY eess.SY math.IT
• 发表时间: 2025年10月（arXiv预印本）
• 论文链接: https://arxiv.org/abs/2510.04379

摘要

距离继电器用于检测输电线路故障，但面临故障位置、故障电阻以及远端电流等不确定性。本文使用Minkowski和来聚合这些不确定性，从而能够显式建模继电器线路周围的电网，并适应同步机和逆变器资源的任意组合。为了简化继电器的任务，逆变器可以注入扰动或辅助信号（如负序电流）。作者使用Farkas引理构建了设计逆变器辅助信号的优化问题。

研究背景与动机

问题描述

1. 传统距离保护面临的挑战：距离继电器通过测量本地电压和电流计算到故障点的阻抗来检测故障，但面临故障位置、故障电阻和远端电流的不确定性。
2. 逆变器资源带来的新问题：随着逆变器基础资源(IBR)在电网中占比增加，传统基于同步发电机假设的距离保护特性可能出现误动作。
3. 现有方法的局限性：
   • 传统四边形特性基于物理启发的临时调整来处理远端电流
   • 这些调整假设电源为同步发电机，在IBR主导的电网中可能失效
   • 缺乏系统性的方法来设计IBR辅助信号

研究动机

作者旨在通过几何方法（Minkowski和）形式化距离保护中的不确定性聚合，并开发优化方法设计IBR辅助信号，以提高距离保护在现代电网中的可靠性。

核心贡献

1. 理论贡献：首次使用Minkowski和形式化距离继电器特性中的不确定性聚合，为传统四边形特性提供了数学基础。
2. 建模创新：提出了能够显式建模电网和电源的框架，可适应同步机和IBR的任意组合。
3. 优化方法：基于Farkas引理构建了设计IBR辅助信号的优化问题，确保继电器能够明确区分不同故障场景。
4. 计算方法：开发了高效的计算方法，包括zonotope投影算法和zonogon松弛，使得新特性能够在线计算。

方法详解

任务定义

输入：电网拓扑、电源不确定性、故障参数范围 输出：继电器特性区域、最优辅助信号 目标：最小化辅助信号幅度，同时确保所有故障场景可分离

物理模型

网络建模

• 将同步机建模为电压源：$v_S = v_S^{\circ}$
• 将IBR建模为电流源：$i_C = i_C^{\circ} + \Delta$（包含辅助信号）
• 负载采用阻抗模型或RLC等效电路

故障环路分析

对于不同故障类型，建立基尔霍夫电压定律(KVL)：

接地故障(ag)： $v_L^a = m_z z(i_L^a + ki_L^0) + m_r r_F(i_L^a + i_R^a)$

相间故障(ab)： $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)$

其中$m_z \in [0,1]$为标准化故障距离，$m_r \in [0,1]$为标准化故障电阻。

不确定性建模

电源不确定性

定义不确定性集合： $U = \left\{\begin{bmatrix} i_C^{\circ} \\ v_S^{\circ} \end{bmatrix} + \Sigma\lambda \mid \lambda \in \Lambda\right\}$

其中$\Lambda$可以是：

• 范数有界：$\Lambda = \{\lambda \mid \|\lambda\|_2 \leq 1\}$
• 多面体(zonotope)：$\Lambda = \{\lambda \mid P\lambda \leq 1\}$

Minkowski和聚合

对于故障场景$\eta$，表观阻抗的不确定性集合为： $\tilde{Z}_A^{\eta} = Z_T \oplus \tilde{Z}_F^{\eta}$

其中$Z_T = \{m_z z \mid m_z \in [\underline{m_z}, 1]\}$是线路阻抗集合，$\tilde{Z}_F^{\eta}$是故障阻抗集合。

测试前后不确定性

测试后不确定性（继电器特性）

测试后已知本地电流$i_L$，不确定性集合变为凸集，可高效计算投影到复平面。

测试前不确定性（辅助信号设计）

测试前$i_L$未知，包含双线性项$m_z\lambda$和$m_r\lambda$，需要凸松弛：

松弛方法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\}$

松弛方法2：分离变量$\lambda_{m_z}$和$\lambda_{m_r}$，去除耦合约束。

辅助信号优化

分离条件

定义模型对$(η_1, η_2)$被辅助信号$\Delta$分离，当且仅当： $\breve{P}^{η_1,η_2}(\Delta) = \emptyset$

即两个测试前不确定性集合的交集为空。

优化问题

\begin{align} \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}

Farkas引理应用

使用Farkas引理将空集约束转换为对偶系统的可行性： $\breve{P}^{η_1,η_2}(\Delta) = \emptyset \Leftrightarrow \mathcal{S}^{η_1,η_2}(\Delta) \neq \emptyset$

最终优化问题变为双凸规划，使用ADMM算法求解。

实验设置

测试系统

• 基础网络：IEEE 14节点系统
• 电源配置：
  • 同步机节点：$\mathcal{S} = \{1,3,5,14\}$
  • IBR节点：$\mathcal{C} = \{2,4,7,8,12\}$
  • 负载节点：$\{6,10,11,13\}$

参数设置

• 同步机：单位幅值电压源，相角均匀分布在$[0, 2π/10]$
• IBR：单位幅值电流源，相角均匀分布在$[0, 2π/3]$
• 噪声：每个源添加$\lambda_k[1, α, α^2]^T$，$λ_k \in \Lambda_k$（正多边形）
• 故障参数：$k=2$（零序补偿系数），$r_F=1$，$\underline{m_z}=0.15$

计算实现

• 语言：Python
• 可视化：Matplotlib
• 优化器：CVXPy + Clarabel求解器
• 计算时间：
  • 所有11种故障的$\tilde{Z}_A^{\eta}$：约0.9秒
  • Zonogon松弛$\hat{Z}_A^{\eta}$：约1毫秒

实验结果

测试后不确定性集合

图4展示了ag和ab故障的精确不确定性集合$\tilde{Z}_A^{\eta}$及其zonogon松弛$\hat{Z}_A^{\eta}$。结果显示：

• 松弛能够很好地逼近精确集合
• 添加切割约束$\text{Im}[z]\text{Re}[z_A^{\eta}] \geq \text{Re}[z]\text{Im}[z_A^{\eta}]$可以收紧松弛

辅助信号设计结果

可视化分析

图5展示了分离不同故障对的辅助信号集合：

• 分离(N,ag)、(N,ab)和(ag,ab)的最小辅助信号：$δ = j0.2$，$\|\Delta\|_2 ≈ 0.53$
• 分离所有LG和LL故障的最小辅助信号：$δ = 1.6+j0.6$，$\|\Delta\|_2 ≈ 4.52$

优化结果

使用ADMM算法优化得到： $\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$

最终幅值$\|\Delta\|_2 ≈ 4.35$，相比网格搜索结果提升4%。

与传统方法对比

• 四边形特性：基于临时调整，假设同步机主导
• 本文方法：精确的Minkowski和，适应任意电源组合
• 计算效率：zonogon松弛可在1毫秒内完成，支持在线计算

相关工作

距离保护基础

• 传统距离保护基于阻抗测量和四边形特性
• 在IBR主导电网中面临误动作问题

IBR对保护的影响

• IBR的控制特性与同步机不同，导致故障电流特性变化
• 现有研究主要关注控制策略调整

辅助信号方法

• 谐波注入、负序电流注入等方法
• IEEE Standard 2800规定IBR应注入90°超前的负序电流

几何方法在保护中的应用

• 本文首次将Minkowski和应用于距离保护
• 提供了系统性的不确定性聚合框架

结论与讨论

主要结论

1. 理论贡献：Minkowski和为距离保护不确定性聚合提供了严格的数学框架
2. 实用价值：新方法能适应现代电网的复杂电源组合
3. 计算可行：zonogon松弛使得在线计算成为可能
4. 优化有效：基于Farkas引理的辅助信号设计方法能够确保故障可分离性

局限性

1. 模型简化：
   • 电源建模相对粗糙（电压源/电流源）
   • 使用标称故障参数计算网络矩阵
   • 忽略了动态过程
2. 计算复杂性：
   • 测试前不确定性集合的双线性特性需要松弛
   • ADMM算法的收敛性依赖于初始化
3. 实用性考虑：
   • 需要准确的网络模型和电源参数
   • 辅助信号的实际实现可能面临约束

未来方向

1. 电磁暂态仿真验证：在详细的电力系统仿真中验证方法有效性
2. 网络矩阵改进：考虑多个标称故障参数或自适应更新
3. 多继电器协调：扩展到多条线路的协调保护
4. 其他保护方案：将核心思想应用于过流保护、时域保护等

深度评价

优点

1. 理论创新性强：
   • 首次将Minkowski和引入距离保护领域
   • 提供了严格的数学框架替代临时调整方法
   • Farkas引理的应用巧妙且有效
2. 实用价值高：
   • 解决了IBR主导电网中的实际问题
   • 计算方法高效，支持在线应用
   • 辅助信号设计系统性强
3. 方法完整性：
   • 涵盖了11种故障类型的完整分析
   • 提供了多种凸松弛方法
   • 算法实现详细且开源

不足

1. 模型假设限制：
   • 电源模型过于简化，可能影响实际应用
   • 忽略了电网动态特性
   • 对中间馈线的处理不够充分
2. 验证不足：
   • 缺乏实际电网数据验证
   • 未与商业继电器进行对比
   • 电磁暂态仿真验证有待完成
3. 计算复杂性：
   • 双线性优化问题的全局最优性难以保证
   • 大规模网络的计算可扩展性未充分验证

影响力评估

1. 学术影响：
   • 为电力系统保护领域引入了新的数学工具
   • 可能启发更多几何方法在保护中的应用
   • 对IBR保护研究具有重要参考价值
2. 工程价值：
   • 为现代电网保护设备升级提供了理论基础
   • 辅助信号设计方法具有实际应用潜力
   • 有助于制定相关标准和规范

适用场景

1. IBR渗透率高的电网：特别适用于可再生能源比例较高的系统
2. 新建输电线路：可在设计阶段考虑新的保护策略
3. 保护设备升级：为现有继电器的智能化升级提供方案
4. 电网规划：在电网规划中考虑保护性能优化

参考文献

论文引用了26篇重要文献，主要包括：

• 距离保护经典教材：Ziegler (2011), Horowitz等 (2022)
• IBR对保护影响的最新研究：Kasztenny (2022), Baeckeland等 (2022)
• 辅助信号相关工作：Saleh等 (2020), Banaiemoqadam等 (2019)
• 数学工具相关：Minkowski和理论、Farkas引理、zonotope理论

总体评价：这是一篇高质量的学术论文，在理论创新和实用价值方面都有重要贡献。作者巧妙地将几何方法引入传统的电力系统保护领域，为解决现代电网面临的保护挑战提供了新的思路。尽管在模型假设和验证方面还有改进空间，但其核心思想和方法框架具有重要的学术价值和工程应用前景。