1. 项目概述当距离保护遇上新能源并网距离保护这个在电力系统继电保护领域服役了数十年的“老兵”其核心逻辑清晰而优雅测量保护安装处的电压与电流计算其比值得到视在阻抗这个阻抗值理论上与故障点到保护安装处的距离成正比。只要这个阻抗落在预先设定的“动作区”内继电器就判定为区内故障迅速跳闸隔离故障。这套逻辑在由同步发电机主导的传统电网中经过精心整定和补偿表现一直相当可靠。然而随着以光伏、风电为代表的逆变器并网资源大规模接入电网的“底色”正在发生根本性变化。IBR通过电力电子变流器与电网连接其故障响应特性与旋转的同步发电机截然不同。同步机在故障瞬间能提供巨大的、受其暂态电抗限制的短路电流且相位相对明确。而IBR受其电力电子器件过流能力限制通常在故障后会迅速进入“电流限幅”模式输出一个幅值受限、相位可能快速变化的电流。这个电流的幅值可能远低于同步机的短路电流且其等效内阻抗也不再是简单的固定电抗而是一个受控源模型。这对传统距离保护构成了严峻挑战。保护的动作边界即特性原本是基于确定的系统阻抗和故障电流相位假设来整定的。当IBR提供的故障电流幅值小、相位不确定时计算出的视在阻抗就会严重偏离真实故障距离导致两种危险情况一是保护拒动该跳不跳故障无法快速切除二是保护误动不该跳却跳造成非故障线路停电扩大事故范围。问题的根源在于传统方法的视在阻抗计算严重依赖于故障时刻系统的“运行点”——即全网各处的电压和电流状态而这些状态在含有大量IBR的电网中变得难以预测。正是在此背景下“增量量”技术重新进入了保护工程师的视野。其核心思想堪称“以不变应万变”既然故障后的电气量是故障分量与正常负荷分量的叠加而令人头疼的不确定性主要来自后者那么何不直接将其剔除增量量就是故障后某一时刻的电气量电压或电流减去故障前一个工频周期的电气量。通过这个简单的减法理论上我们就得到了纯由故障引起的“故障分量”。这个分量有一个美妙的性质它主要取决于电网的网络结构线路阻抗、拓扑和故障本身的位置与过渡电阻而与故障前系统带了多大负荷、各电源输出多少功率这个“运行点”几乎无关。这就为构建一个更稳健、更适应新型电力系统的距离保护特性提供了可能。本文将深入拆解一种基于增量量构建距离保护特性的新方法。我们不仅会阐释其背后的数学模型和物理意义更会聚焦于工程实现中的关键细节、参数整定的实际考量以及面对IBR时该方法的优势与仍需注意的局限性。无论你是从事继电保护整定计算的工程师还是研究新型电力系统稳定性的技术人员理解这套方法都将有助于你应对当前电网保护领域最前沿的挑战。2. 增量量距离保护的数学模型与核心思想拆解要理解基于增量量的距离保护必须从其数学模型入手看清它是如何巧妙地“剥离”运行点影响的。整个过程可以看作一个严谨的推导其目标是将继电器测量到的视在阻抗用仅与网络结构和故障参数相关的增量量来表达。2.1 增量量的定义与基本假设首先我们明确定义增量量。设电力系统工频周期为 δ。对于网络中任意母线k其在时间窗口 [t, tδ) 内的三相电压和电流相量记为 v_k(t) 和 i_k(t)。那么在故障发生时刻t我们假设故障发生在t时刻对应的增量电压和增量电流定义为˜v_k v_k(t) - v_k(t-δ)˜i_k i_k(t) - i_k(t-δ)这里我们选取故障前一个周期的电气量作为基准。理论上p可以取任意正整数但p1是最常用且合理的因为它能最大程度保证故障前后网络拓扑未变同时获取最新的“正常状态”信息。整个推导建立在一条关键假设之上所有电源在故障发生前后短暂时间内至少涵盖[t-δ, tδ]这个窗口保持周期性。这意味着对于同步发电机SG其机端电压 v_S(t) v_S(t-δ)。这符合物理实际因为同步机转子惯性大故障后几个周期内其电动势幅值和相位不会突变。对于逆变器并网资源IBR建模为诺顿等效电路一个电流源并联一个导纳则假设其注入电流 i_C(t) i_C(t-δ)。这是该方法面对IBR时最主要的近似也是其局限性的来源。我们后文会详细讨论。这个假设是整个数学魔术的“开关”。它意味着在故障前后所有电源的“源项”SG的电压源电动势、IBR的电流源输出在相量意义上没有变化。那么故障前后电气量的任何变化就只能归因于网络中出现了新的扰动——即故障。2.2 网络建模与增量方程推导考虑一个典型的保护场景一条输电线路其正序串联阻抗为z。保护继电器安装于线路的本侧母线L测量电压v_L和流入线路的电流i_L。对侧为远端母线R。故障发生在线上故障点距本侧的距离占线路全长的比例为 m_T (0 ≤ m_T ≤ 1)故障过渡电阻为 m_F * R_F其中R_F是最大可能过渡电阻m_F ∈ [0, 1]。为了分析远端电流i_R对保护测量的影响必须对全网进行建模。将网络中的所有母线分类C: IBR母线建模为电流源 i_C 并联导纳 Y_CS: 同步发电机母线建模为电压源 v_SJ: 联络节点和负荷母线建模为导纳 Y_J对纯联络节点Y_J0F: 虚拟的故障点母线根据基尔霍夫定律可以列写故障后瞬间(t时刻)全网络的节点电压方程Y(m_T) * [v_F(t); v_J(t); v_C(t); v_S(t)]^T [-Y_F^η(m_F)v_F(t); -Y_J v_J(t); i_C(t)Y_C v_C(t); i_S(t)]其中Y(m_T) 是包含了位于m_T处的虚拟故障节点的系统节点导纳矩阵。方程右侧是注入各节点的电流向量其中故障点注入电流为-Y_F^η v_Fη代表故障类型如A相接地ag、AB相间短路ab等。同样可以列写故障前一个周期(t-δ时刻)的方程。关键的一步来了将故障后的方程减去故障前的方程。由于假设了电源的周期性v_S(t)v_S(t-δ) i_C(t)i_C(t-δ)方程右侧的源项全部被消去同时故障前故障点注入电流为0。相减后我们得到Y(m_T) * [˜v_F; ˜v_J; ˜v_C; 0]^T [-Y_F^η(m_F)v_F(t); -Y_J ˜v_J; Y_C ˜v_C; ˜i_S]经过一系列矩阵变换和推导核心是利用矩阵求逆将故障点电压用继电器测量量表示我们可以得到一个极其重要的关系式——增量远端电流表达式σ_η(m) ˜i_R - ˜v_F / ((1-m_T)z) Ω_η(m) * [v_L(t-δ); i_L(t-δ)]^T这个公式(3)是理解整个方法的关键。它表明故障后线路对侧的增量电流σ_η(m)可以通过一个仅与网络结构、故障类型η、故障位置m_T和过渡电阻m_F相关的矩阵Ω_η(m)乘以故障前一个周期继电器本地测量到的电压和电流来计算。Ω_η(m)矩阵可以离线预先计算好。这意味着在故障发生后继电器无需知道系统实时的运行状态只需要利用故障前一刻这是已知的历史数据的本地信息结合预存的网络参数矩阵就能推算出故障后对侧电流的增量分量。实操心得矩阵Ω_η(m)的离线计算与存储这是工程实现中的核心步骤。Ω_η(m)依赖于网络拓扑和参数Y矩阵以及故障位置m_T。在实际系统中拓扑相对固定但运行方式多变。我们需要针对一条线路可能发生的各种故障类型η以及一系列离散的m_T和m_F取值例如m_T从0到1以0.05为步长m_F取0, 0.5, 1等几个典型值预先计算出对应的Ω_η(m)矩阵并存储在继电保护装置中。这需要保护装置具备一定的存储和快速检索能力。计算Ω_η时需要完整的系统正序、零序参数这依赖于前期的系统参数普查和建模工作。2.3 从增量量到视在阻抗以A相接地故障为例有了增量量的表达我们就可以重新审视距离保护测量的视在阻抗。以最常见的A相接地故障ag为例。传统距离保护中继电器计算的视在阻抗为z_A^ag v_A^ag(t) / i_A^ag(t)其中v_A^ag(t) v_L^a(t)A相电压i_A^ag(t) i_L^a(t) k * i_L^0(t)A相电流加零序补偿电流k为零序补偿系数。根据故障回路分析从继电器安装处到故障点的电压降可以写成v_L^a(t) m_T * z * i_A^ag(t) m_F * r_F * (i_L^a(t) i_R^a(t))其中i_L^a(t) i_R^a(t)就是流过故障过渡电阻的故障电流。将故障电流用增量量表示i_L^a(t) i_R^a(t) (˜i_L^a ˜i_R^a) (i_L^a(t-δ) i_R^a(t-δ))。由于故障前线上该相电流和为零假设是三相对称负荷第二项为0。因此故障电流简化为˜i_L^a ˜i_R^a。而我们已经知道增量远端电流˜i_R^a可以通过σ_ag(m)和故障点增量电压˜v_F^a来表达。经过代入和整理具体推导涉及将˜v_F也用本地测量量表示最终可以得到基于增量量的视在阻抗表达式z_A^ag(m) m_T * z m_F * r_F * (˜i_L^a ψ_ag * σ_ag(m)) / (i_L^a(t) k * i_L^0(t))这个公式(4)看起来比传统公式复杂但其物理意义非常深刻。公式右边第一项m_T * z就是真实的线路阻抗。第二项是过渡电阻引入的附加项。最关键的是这个附加项中的所有变量要么是本地可测的增量量˜i_L^a, ˜i_L^0要么是通过离线矩阵Ω和故障前本地量计算出的σ_ag(m)。整个表达式z_A^ag(m)不再显式依赖于故障时刻全网未知的运行状态如远端电源的电压、其他负荷的大小等而只依赖于故障位置m_T和过渡电阻m_F未知待判定线路参数z,r_F已知本地故障前后电气量v_L(t-δ), i_L(t-δ), v_L(t), i_L(t)已知网络结构矩阵Ω_η(m)离线已知这就实现了设计目标将距离保护的动作特性即视在阻抗z_A随m_T和m_F变化的轨迹从对系统运行点的依赖中解放出来使其主要取决于网络本身和故障参数。对于AB相间短路ab故障可以采用完全类似的思路推导出对应的公式(5)。3. 保护特性的构建从理论公式到实用判据上一节我们得到了视在阻抗z_A^η(m)的表达式它是一个关于故障位置m_T和过渡电阻m_F的复函数。距离保护的动作特性理论上就是当m_T在保护区段内例如0到0.8倍线路全长且m_F在可能范围内0到某个最大值变化时z_A^η(m)在复阻抗平面上扫过的区域。如果实际测量计算出的视在阻抗落在这个区域内就判定为区内故障。3.1 精确特性集的困境与近似必要性定义特性集Z_η { z_A^η(m) | m_T ∈ [0, 1], m_F ∈ [0, 1] }。理想情况下我们希望得到Z_η的精确形状。但观察公式(4)和(5)会发现由于σ_η(m)本身是m_T和m_F的复杂函数通过矩阵Ω_η(m)z_A^η(m)与m的关系是非线性的。这意味着Z_η很可能是一个形状不规则的区域很难用一个简单的几何图形如传统的多边形、圆形或苹果形特性来精确描述。在保护装置中实时计算并判断一个点是否落在这样一个复杂区域内计算量过大难以满足速动性要求。因此论文中提出了两种实用的近似方法旨在用简单的几何形状来覆盖或近似这个复杂的特性集Z_η。3.2 点估计法平行四边形近似第一种方法称为“点估计法”。其思路非常直观既然σ_η(m)随m变化导致非线性我们干脆用一个固定的、代表性的m值来计算它从而简化表达式。我们选择一个标称的故障参数估计值ˆm (ˆm_T, ˆm_F)。一个合理的选择是ˆm_T 0.5故障在线路中点和ˆm_F 1过渡电阻为最大可能值。然后用这个固定的ˆm来计算σ_η(ˆm)代入视在阻抗公式得到简化后的表达式ẑ_A^η(m, ˆm)。以A相接地故障为例ẑ_A^ag(m, ˆm) m_T * z m_F * r_F * (˜i_L^a ψ_ag * σ_ag(ˆm)) / (i_L^a(t) k * i_L^0(t))此时ẑ_A^η(m, ˆm)是m_T和m_F的线性函数。当m_T从0变到1时对应复平面上沿着向量z方向的一条线段。当m_F从0变到1时对应另一条沿着向量(˜i_L^a ψ_ag * σ_ag(ˆm)) / (i_L^a(t) k * i_L^0(t))方向的线段。这两个线段的闵可夫斯基和即一个线段上的每个点与另一个线段上的每个点相加所得到的全部点的集合就构成了一个平行四边形。这个平行四边形就是近似后的保护特性ˆZ_η(ˆm)。它的优点是计算极其简单快速仅需一次σ_η(ˆm)计算和向量加法形状规则易于在继电器中实现判据。只要实际的精确特性集Z_η大部分落在这个平行四边形内且平行四边形本身没有过多地超出线路实际阻抗范围这个近似就是可用的。注意事项点估计法的参数选择与适应性选择ˆm_T0.5和ˆm_F1是一种折中。ˆm_F1考虑了最严重的过渡电阻情况保证了在高阻故障时特性有足够的覆盖范围。ˆm_T0.5则试图用一个“平均”位置来代表全线。然而当实际故障点靠近线路两端或者过渡电阻很小时用中点计算出的σ_η可能与实际情况有偏差导致平行四边形特性与真实特性存在误差。这可能会略微影响保护的灵敏度或选择性。在实际整定时可能需要根据具体线路参数和系统结构通过仿真来验证这个标称点的选择是否合理或者针对不同保护区段如第I段、第II段采用不同的ˆm值。3.3 凸包近似法更精确的拟合第二种方法旨在获得更接近Z_η真实形状的近似。其思路是既然无法解析描述整个区域我们就用足够多的采样点来“描绘”它。我们在(m_T, m_F)的定义域[0,1]×[0,1]内选取一系列离散点ˆm_kk1, 2, ..., n。对于每一个点我们计算其对应的精确视在阻抗z_A^η(ˆm_k)。这样就在复阻抗平面上得到了n个点。然后我们计算这n个点的凸包。凸包是包含这组点的最小凸多边形。如果采样点ˆm_k选得足够多、足够密特别是能覆盖m_T和m_F变化时对z_A^η影响最大的区域那么这个凸包Z_η_ch就能很好地逼近真实的特性集Z_η。论文中提到一个简单的采样方案是取四个角点(0,0), (0,1), (1,0), (1,1)再加上在区域内均匀布点例如一个4×5的网格。凸包近似的精度高于平行四边形近似因为它更直接地基于精确公式采样。计算凸包有成熟的算法如礼品包装算法、Graham扫描法等对于几十个点的集合计算速度在现代保护装置的处理器上是可以接受的论文中提到约3毫秒。实操心得凸包算法的工程实现与优化在嵌入式保护装置中实现凸包算法需要考虑实时性和可靠性。礼品包装算法虽然直观但最坏情况下的时间复杂度为O(nh)其中n是点数h是凸包上的点数。对于点数不多如20-30个的情况完全可行。为了进一步加速可以离线预计算对于固定的网络拓扑和线路参数Ω_η(m)矩阵是固定的。我们可以离线计算出所有采样点ˆm_k对应的z_A^η(ˆm_k)并将这些点的坐标直接存储在保护装置中。在线应用时只需要根据故障类型η调用对应的点集然后在线计算这些点在当前故障量(˜i_L, i_L(t))作用下的实际位置这通常是一个简单的复数乘法或加法再计算这些新点的凸包。这大大减少了在线计算量。采样点优化不必在整个区域均匀采样可以重点采样m_T在保护区段内如0-0.8和m_F在常见高阻范围内的点这样可以用更少的点获得对动作区更好的拟合。判据简化不一定需要实时生成完整的凸包多边形。可以预先存储凸包多边形的边向量方程。在线计算得到测量阻抗z_meas后判断z_meas是否在凸包内可以转化为判断z_meas是否在所有边向量的同一侧例如都在左侧这是一个快速的线性判断。4. 面对逆变器并网资源的优势、挑战与工程考量基于增量量的距离保护方法其最大的吸引力在于应对高比例IBR电网的保护难题。让我们深入分析其优势所在以及在实际工程应用中必须面对的挑战和需要进行的适应性调整。4.1 应对IBR故障特性的核心优势传统距离保护在IBR附近失灵根源在于IBR的故障响应与传统电源不同导致传统的电源模型固定电压源串联阻抗和故障电流相位假设失效。增量量方法从原理上规避了这个问题消除运行点依赖性这是最根本的优势。保护特性不再依赖于故障时刻全网未知的电压电流分布运行点而只依赖于网络结构和故障参数。无论IBR在故障前输出多少功率故障后进入何种控制模式只要其输出电流在故障后短期内满足或近似满足周期性假设该方法都能建立一个相对稳定的阻抗-距离关系。这极大地减少了因系统运行方式变化而需要频繁修改保护定值的需求。减弱对故障电流幅值与相位的敏感度传统距离保护中过渡电阻附加项ΔR R_F * (I_L I_R) / I_L受对侧助增电流I_R影响巨大。IBR提供的助增电流幅值小且相位不确定会严重扭曲这个附加项导致阻抗测量不准。在增量量方法中对应的项变为m_F * r_F * (˜i_L ψη * σ_η(m)) / i_A(t)。其中σ_η(m)是通过网络结构矩阵和故障前量计算得出的它反映的是网络对故障分量的响应而不是对侧电源的绝对出力。这在一定程度上“规范化”了助增电流的影响使其更可预测。对不对称故障的适应性增量量方法分别处理不同类型的短路故障通过不同的η。对于IBR其在不对称故障下的正序、负序、零序输出特性可能更加复杂与传统同步机差异更大。增量量方法通过分别建模可以更精细地处理不同序分量的影响而不是依赖于传统的基于同步机假设的复合序网连接。4.2 关键假设的挑战与应对策略该方法的核心假设——电源周期性——在IBR面前显得尤为脆弱这也是该方法最主要的挑战。幅值突变IBR在检测到故障后会迅速通常在几个毫秒内进入低电压穿越模式或电流限幅模式。其输出电流的幅值会从额定值突然增加到限幅值。这意味着i_C(t)很可能不等于i_C(t-δ)。这破坏了增量量中抵消源项的基础。影响分析这会导致推导中˜i_S项不为零从而使得σ_η(m)的表达式中混入了未抵消的源项不确定性。最终视在阻抗z_A^η(m)将再次受到对侧IBR故障后输出电流变化的影响。应对思路论文中提到一种补救措施是将这个不确定性显式地纳入模型即认为i_L(t)本身也包含一个未知但有界的扰动。这样保护特性Z_η会从一个二维区域随m_T,m_F变化扩展为一个更高维的、更复杂的区域计算和判断将变得困难。更实际的工程思路可能是利用IBR控制响应的时间尺度。IBR的电流控制环响应极快但其电流指令值的突变并非瞬时完成。如果保护的动作速度足够快例如采用基于暂态量的高速保护算法在故障后半个周期内做出判断有可能在IBR电流发生大幅变化之前就完成测量和判断此时周期性假设近似成立。这要求保护装置具有极高的数据采样率和处理速度。相位与频率变化故障可能导致并网点电压相位跳变IBR的锁相环需要跟踪这个跳变期间其输出电流的相位也可能发生快速变化。这同样违反了周期性假设。影响分析相量的计算依赖于对信号频率和相位的准确估计。如果故障前后工频频率或相位基准发生变化那么直接相减得到的“增量量”就不再是纯故障分量还包含了频率/相位变化引入的误差。应对思路需要采用频率自适应的相量提取算法。例如在计算v_k(t)和v_k(t-δ)时不能假设两者频率完全相同。可以采用基于DFT或卡尔曼滤波的动态相量估计算法实时跟踪系统频率并确保在相减时使用的是基于同一频率基准计算的相量。或者考虑使用时域差分等受频率影响更小的算法来近似故障分量。4.3 工程实现中的具体问题与解决方案数据窗与速动性的权衡该方法需要故障前一个周期20ms的数据。这意味着保护至少需要在故障发生20ms后才能获得完整的增量量数据窗。对于要求瞬时动作如距离I段通常要求20-30ms内动作的保护来说这会引入固有延时。一种改进思路是使用更短的故障前数据窗例如半周波并结合预测或插值算法来估计一个完整周期的分量但这会引入误差。对CT/PT测量误差的敏感性增量量是两个大数故障后量减去一个相近的大数故障前量得到一个小数故障分量。这个过程会放大测量误差的相对影响。特别是电流互感器在故障大电流下的饱和可能导致故障前后电流测量均出现畸变其差值增量量的误差可能非常大。解决方案必须采用高精度、抗饱和的测量互感器。同时在算法中加入数据有效性校验和抗干扰模块。例如可以计算故障前后电压电流的突变量启动判据只有突变量超过一定门槛值才启用增量量计算避免在正常轻微波动下进行无意义的相减放大噪声。系统拓扑变化与矩阵Ω的更新该方法依赖于准确的网络导纳矩阵Y来计算Ω_η(m)。如果系统运行方式发生重大变化如重要线路或变压器投退导致网络拓扑改变Ω_η(m)矩阵必须相应更新。解决方案这需要保护装置具备与广域测量系统或调度主站通信的能力接收网络拓扑变化信息并调用预存的不同运行方式下的Ω_η矩阵集。或者研究基于本地测量实时辨识等效系统阻抗的方法动态修正保护特性但这在算法上更具挑战性。高阻故障的灵敏度与所有距离保护一样高阻故障m_F接近1时过渡电阻上的压降占主导测量阻抗z_A的电阻分量会大幅增加可能使阻抗点落在保护特性区域之外导致拒动。增量量方法并未从根本上解决高阻故障检测难题但因其特性更稳定可能为优化高阻故障判据如结合零序功率方向、谐波分量等提供一个更好的基础平台。5. 总结与展望一种有潜力的方向基于增量量的距离保护特性为解决IBR并网带来的保护适应性问题提供了一条富有潜力的新路径。它通过巧妙的数学处理将保护特性的影响因素从难以预测的系统全局运行点转移到了相对固定的网络结构和本地可测或可推算的量上从而增强了保护的鲁棒性。我个人认为这种方法的价值不仅在于其本身更在于它启发了我们思考保护问题的新角度从“绝对量”的保护转向“变化量”的保护。在系统惯量降低、电源特性多元化的未来电网中故障前的稳态运行点可能千变万化但故障瞬间产生的“突变”所遵循的网络物理规律相对而言更为稳定。增量量正是抓住了这个“突变”特征。然而将其投入工程实用仍有很长的路要走。当前的研究和实验多集中在理论推导和离线仿真验证。下一步的关键在于实验室动模与RTDS测试在实时数字仿真器中接入详细的IBR电磁暂态模型全面测试该方法在各种故障类型、不同IBR渗透率、不同控制策略下的性能特别是验证其对于IBR周期性假设破坏的容忍度。与现有保护体系的融合它不太可能完全取代成熟的传统距离保护更可能作为一种附加判据或自适应模块。例如在传统距离保护因IBR影响而动作模糊时启用增量量判据进行辅助决策或者利用增量量特性来动态调整传统保护的动作边界。算法优化与硬件实现研究更快速的凸包生成与判据算法优化Ω_η矩阵的存储与计算开发适应频率变化的增量量计算模块使其能够嵌入到现有保护装置的硬件平台中满足实时性要求。总之基于增量量的距离保护是一种面向未来的有益探索。它提醒我们在新型电力系统背景下继电保护不能止步于沿用旧模型和旧定值必须深入理解新元件的物理本质从基本原理出发创新保护原理与算法。这条路充满挑战但也正是技术发展的魅力所在。