技术领域
[0001] 本
发明属于电
力系统领域，特别涉及一种频率变化情况下序分量计算的方法，并涉及应用该方法的继电保护装置和安全自动装置。
背景技术
[0002] 电力系统是基于额定频率（50Hz或者60Hz）运行的，《GB/T15945-2008
电能质量电力系统
频率偏差》第3.1条规定：“电力系统正常运行条件下频率偏差限值为±0.2Hz。当系统容量较小时，偏差限值可以放宽到±0.5Hz。”。序分量计算是电力系统应用领域的基本
算法之一，由于电力系统频率变化很小，因此以往的序分量计算都是基于系统额定频率进行的，在频率偏离额定值较多或者变化时序分量计算将出现错误结果，文献《计算机继电保护原理与技术》（陈德树主编，
水利电力出版社1992年出版）和《
变压器及中低压网络数字式保护》（许正亚编著，中国水利
水电出版社2004年出版）等给出的序分量计算方法的隐含条件均是电力系统在额定频率下运行。但在一些新的应用场合，如抽水蓄能机组和
燃气轮机组的变频启动过程，或系统解列事故导致发
电机组输出功率突降过程，
发电机组仍处于三相对称运行，但其频率将快速变化或在低频下运行，以往基于额定频率的序分量计算不再适用。
[0003] 基于以上分析，为了克服现有序分量计算基于固定频率的问题，需要一种适用于电力系统工作频率变化情况下的序分量计算方法，本案由此产生。
发明内容
[0004] 本发明的目的，在于提供一种频率变化情况下序分量的计算方法，其可提高频率变化时序分量计算的准确度。
[0005] 为了达成上述目的，本发明的解决方案是：
[0006] 一种频率变化情况下序分量的计算方法，包括如下步骤：
[0007] （1）继电保护装置或安全自动装置对被保护对象的
电压或
电流互感器二次侧电气量进行
采样，采用相间短窗测频算法求取电气量当前频率；
[0008] （2）根据步骤（1）求取的电气量当前频率，求取电气量的正、负、零序分量。
[0009] 上述步骤（1）的具体内容是：
[0010] （1.1）继电保护装置或安全自动装置对被保护对象的电压或电流互感器二次侧电气量进行采样，得到电气量二次值数据序列ea(n)、eb(n)、ec(n)，其中，符号n表示序列号，下标a、b、c分别表示电气量A相、B相、C相；由前述三个相电气量数据序列得到三个线电气量数据序列：
[0011]
[0012] 其中，下标ab、bc、ca分别表示电气量AB相、BC相、CA相；
[0013] （1.2）设当前
采样频率为fS，则采样间隔为
选定ex相和ey相进行频率计算，ex和ey从ea(n)、eb(n)、ec(n)、eab(n)、ebc(n)、eca(n)六个电气量中选择，ex领先于ey的
相位差为Δθ，Δθ<2π；设ex相从负到正过零点的前一个采样数据为ex(p-1)，后一个采样数据为ex(p)；设ey相从负到正过零点的前一个采样数据为ey(m-1)，后一个采样数据为ey(m)，p和m分别表示序列号，则电气量当前周期Te和频率fe为
[0014]
[0015] 其中，K表示ex相与ey相过零点之间的采样数。
[0016] 上述步骤（2）的具体内容是：
[0017] （2a）根据当前采样频率fS和电气量当前频率fe，计算当前每周波的采样点数N，然后采用变数据窗傅氏算法求取三相电气量的
实部和
虚部；
[0018] 设X相的实部eXr和虚部eXi，X相表示A相、B相和C相中的任意一相，计算公式为：
[0019]
[0020] 其中，eX(n)表示X相电气量数据序列，[]表示四舍五入取整；
[0021] 因此得到电气量相量 如下：
[0022]
[0023] （2b）根据序分量计算公式，计算电气量零序分量 正序分量 和负序分量 ：
[0024]
[0025] 其中，eY表示Y相的电气量，Y相表示X相的正相序排列的下一相；eZ表示Z相的j120°电气量，Z相表示Y相的正相序排列的下一相；α表示矢量运算算子e 。
[0026] 上述步骤（2）的具体内容是：
[0028]
[0029] 其中，kmin和kmax分别为最小、最大调节速率系数，ΔfS为当前采样频率fS的参考量与电气量当前频率fe偏差，ΔfSmax为最大频率偏差，N为每周波的采样点数，Δf'S为采样频率调节量；
[0030] （2B）调整采样频率，使当前采样频率的参考量与电气量当前频率偏差ΔfS的绝对值不大于固定
门槛：
[0031] |ΔfS|≤ΔfS·set
[0032] 其中，ΔfS·set表示频率调整偏差固定门槛；
[0033] （2C）采样频率调整完成后，根据A、B、C三相采样数据计算各序分量。
[0034] 上述步骤（2C）的具体内容是：采用固定数据窗傅氏算法求取三相电气量的实部和虚部，进而采用序分量计算公式得到各序分量；
[0035] 首先采用固定数据窗傅氏算法求取三相电气量的实部和虚部,设X相的实部eXr和虚部eXi，X相表示A相、B相和C相中的任意一相，计算公式为：
[0036]
[0037] 其中，eX(n)表示X相电气量数据序列；N为常数，表示每周波采样数；
[0038] 从而得到电气量相量 如下：
[0039]
[0040] 最终采用序分量计算公式得到零序分量 正序分量 和负序分量 ：
[0041]
[0042] 其中，eY表示Y相的电气量，Y相表示X相的正相序排列的下一相；eZ表示Z相的j120°电气量，Z相表示Y相的正相序排列的下一相；α表示矢量运算算子e 。
[0043] 上述步骤（2C）的具体内容是：先采用采样值滤序方法滤出电气量的正、负、零序分量采样序列，再采用固定数据窗傅氏算法求取各序分量；
[0044] 首先采用采样值滤序方法滤出电气量的正、负、零序分量采样序列，设X相的正、负、零序分量序列分别为eX1(n)、eX2(n)、eX0(n)，X相表示A相、B相和C相中的任意一相：
[0045]
[0046] 其中，eX(n)、eY(n)、eZ(n)分别表示X相、Y相、Z相电气量数据序列，Y相表示X相的正相序排列的下一相；Z相表示Y相的正相序排列的下一相；N为常数，表示每周波采样数。
[0047] 采用上述方案后，本发明快速准确地计算电气量当前频率，进而采用变数据窗序分量算法或非线性加速频率
跟踪算法消除了序分量计算受频率变化的影响，解决了以往序分量计算因受频率变化影响而计算不准确的问题，提高了序分量计算在频率变化时的准确度，且方法易于实现。
附图说明
[0049] 图2是电力系统正常运行时各相相位关系图；
[0050] 图3是线性拟合过零点算法示意图；
[0051] 图4是相间短窗测频算法的示意图（采用A、B相短窗）；
[0052] 图5是相间短窗测频算法的示意图（采用A、BC相短窗）；
[0053] 图6是相间短窗测频算法的示意图（采用AB、A相短窗）。
具体实施方式
[0054] 以下将结合附图，对本发明的技术方案进行详细说明。
[0055] 本发明提供一种频率变化情况下序分量的计算方法，结合图1所示，包括如下步骤：
[0056] （1）继电保护装置或安全自动装置测量被保护对象的电压或电流互感器的二次侧电气量，并求取电气量当前频率。
[0057] 频率测量目前比较常用的一种算法是线性拟合过零点算法，文献《变压器及中低压网络数字式保护》等对此有详细阐述，可参考图3所示。线性拟合过零点算法需要一个完整电气量周期Te才能计算出正确的测频结果，所需时间较长。由于电力系统为三相交流系统，正常情况下处于三相对称运行，在非
短路故障情况下各相电气量具有固定对称的相位关系（结合图2所示），基于此本案提出一种改进的加速方法，称之为相间短窗测频方法，可综合利用各相电气量
相位差来提高频率计算速度。
[0058] 继电保护装置或安全自动装置对被保护对象的电压或电流互感器二次侧电气量进行采样，得到电气量二次值数据序列ea(n)、eb(n)、ec(n)，其中，符号n表示序列号，下标a、b、c分别表示电气量A相、B相、C相。由三个相电气量数据序列得到三个线电气量数据序列。
[0059]
[0060] 其中，下标ab、bc、ca分别表示电气量AB相、BC相、CA相。
[0061] 根据电力系统基本知识，我们知道被测电气量ea、eb、ec、eab、ebc、eca各相频率相等，且有固定的相位关系，因此我们可以利用不同相别之间的较短相位差来提高频率计算速度。
[0062] 假定当前采样频率为fS，则采样间隔为 根据需要我们选定ex相和ey相进行频率计算，ex和ey从上述六个电气量中选择，ex领先于ey的相位差为Δθ（Δθ<2π）。设ex相从负到正过零点的前一个采样数据为ex(p-1)，后一个采样数据为ex(p)；设ey相从负到正过零点的前一个采样数据为ey(m-1)，后一个采样数据为ey(m)，p和m分别表示序列号，ex相与ey相过零点之间采样数为K。则电气量当前周期Te和频率fe为
[0063]
[0064] 这样频率计算时间缩短为
[0065] 假定我们选择ea和eb相进行计算，如图4所示，其相位差为 则频率计算时间缩短为 选择ea和ebc相进行计算，配合图5所示，其相位差为 则频率计算时间缩短为选择eab和ea相进行计算，配合图6所示，其相位差为 则频率计算时间缩短为[0066] 相间短窗测频算法能够明显缩短频率计算时间，加快计算速度，但数据窗越短频率测量误差相对越大，可根据应用场合要求选择合适数据窗。另外选择较高数据采样频率时，对提高频率测量
精度有较好效果。
[0067] （2）采用以下两种方法中的任一种求取电气量的正、负、零序分量。
[0068] 第一种方法：采用变数据窗傅氏算法求取三相电气量的实部和虚部，进而利用序分量公式计算出电气量的各序分量，具体步骤是：
[0069] （2a）根据当前采样频率fS和电气量当前频率fe，计算当前每周波的采样点数N，也即变数据窗长度，然后采用傅氏算法求取三相电气量的实部和虚部；
[0070] 设X相的实部eXr和虚部eXi，X相可任意选择A相、B相或C相，计算公式为：
[0071]
[0072] 其中，eX(n)表示X相电气量数据序列，[]表示四舍五入取整；
[0073] 因此得到电气量相量 如下：
[0074]
[0075] （2b）根据序分量计算公式，计算电气量零序分量 正序分量 和负序分量 ，X相表示A相、B相和C相中的任意一相，Y相表示X相正相序排列的下一相，Z相表示Y相正相序排列的下一相：
[0076]
[0077] 其中，α表示矢量运算算子ej120°。
[0078] 第二种方法：采用非线性加速频率跟踪法，使当前采样频率fS快速跟踪电气量当前频率fe，保证二者之比为固定值，即每周波的采样点数N为固定值（N为常数），然后利用傅氏算法和序分量算法求取电气量的各序分量。
[0079] （2A）为加快频率跟踪速度，采用以下非线性加速算法：
[0080]
[0081] 其中，kmin和kmax分别为最小、最大调节速率系数，ΔfS为当前采样频率fS的参考量与电气量当前频率fe偏差，ΔfSmax为最大频率偏差，Δf'S为采样频率调节量。
[0082] （2B）调整采样频率，使当前采样频率的参考量与电气量当前频率偏差ΔfS的绝对值不大于固定门槛：
[0083] |ΔfS|≤ΔfS·set
[0084] 其中，ΔfS·set表示频率调整偏差固定门槛，该定值根据序分量计算精度要求确定。
[0085] （2C）采样频率调整完成后，根据A、B、C三相采样数据计算各序分量，通常可采用以下两种方法：一种是采用固定数据窗傅氏算法求取三相电气量的实部和虚部，进而采用序分量计算公式得到各序分量；
[0086] 首先采用固定数据窗傅氏算法求取三相电气量的实部和虚部,设X相的实部eXr和虚部eXi，X相表示A相、B相和C相中的任意一相，计算公式为：
[0087]
[0088] 其中，eX(n)表示X相电气量数据序列；N为常数，表示每周波采样数。
[0089] 因此得到电气量相量 如下：
[0090]
[0091] 最终采用序分量计算公式得到零序分量 正序分量 和负序分量 X相表示A相、B相和C相中的任意一相，Y相表示X相的正相序排列的下一相，Z相表示Y相的正相序排列的下一相：
[0092]
[0093] 其中，α表示矢量运算算子ej120°。
[0094] 另一种是先采用采样值滤序方法滤出电气量的正、负、零序分量采样序列，再采用固定数据窗傅氏算法求取各序分量。
[0095] 首先采用采样值滤序方法滤出电气量的正、负、零序分量采样序列。设X相的正、负、零序分量序列分别为eX1(n)、eX2(n)、eX0(n)，X相表示A相、B相和C相中的任意一相，Y相表示X相的正相序排列的下一相，Z相表示Y相的正相序排列的下一相。
[0096]
[0097] 其中，eX(n)、eY(n)、eZ(n)分别表示X相、Y相、Z相电气量数据序列；N为常数，表示每周波采样数。
[0098] 需要说明的是，这里给出的采样值滤序方法是一种基本方法，尚有多种改进或
变形方法可参见《变压器及中低压网络数字式保护》（许正亚编著，中国水利水电出版社2004年出版）等参考资料。
[0099] 然后采用固定数据窗傅氏算法求取各序分量的实部和虚部。设X相序分量的实部为eXmr，虚部为eXmi，X相表示A相、B相和C相中的任意一相，m表示某序分量（0为零序，1为正序，2为负序）。
[0100]
[0101] 因此得到电气量的各序分量如下：
[0102]
[0103] 为清楚地说明本发明的方法，本
实施例以发电机组功率突降切机装置应用为例，阐明本发明的具体实施方式。
[0104] 火力
发电厂输电通道发生突然解列事故时，汽轮发电机输出功率将突然降低，转速迅速上升，最高可能超过3300转/分钟，对应频率可能超过55Hz，发电机组功率突降切机装置检测事故特征，快速动作于停机，可靠保障发电机组的运行安全。在这种应用环境下，可以按本发明的方法准确计算发电机电气量的各序分量值，以发电机机端电压序分量计算为例说明如下，类似的其他电气量序分量计算均可应用本方法实现。
[0105] 本发明的具体实施步骤是：
[0106] （1）电压互感器二次侧采样
[0107] 装置对发电机机端电压互感器二次侧电压进行采样，得到发电机机端
三相电压（PT电压二次值）数据序列ua(n)、ub(n)、uc(n)，其中，符号n表示序列号，下标a、b、c分别表示发电机机端电压A相、B相、C相。由三相电压数据序列得到三个线电压数据序列。
[0108]
[0109] 其中，下标ab、bc、ca表示发电机机端线电压AB相、BC相、CA相。
[0110] 注：前文采用符号e表示电压或者电流电气量，此处特指电压量，根据电力专业符号使用习惯采用符号u表示。
[0111] （2）电压频率计算
[0112] 采用相间短窗测频算法计算发电机机端电压频率：
[0113]
[0114] 根据需要采用A相电压和B相电压过零点计算发电机机端电压频率；也可以利用A相电压和BC线电压过零点计算发电机机端电压频率；或者利用AB线电压和A相电压过零点计算发电机机端电压频率。
[0115] （3）序分量计算
[0116] 得到当前发电机机端电压频率fe后，采用变数据窗傅氏算法，求取A、B、C三相电压的实部和虚部：
[0117]
[0118] 则得到发电机机端三相电压 为：
[0119]
[0120] 采用序分量计算公式求取发电机机端电压零序分量 正序分量 和负序分量 ，X相表示A相、B相和C相中的任意一相，Y相表示X相的正相序排列的下一相，Z相表示Y相的正相序排列的下一相：：
[0121]
[0122] 其中，α表示矢量运算算子ej120°。也可以采用非线性加速频率跟踪方法，使装置采样频率和发电机机端电压频率比值固定为常数N，然后采用傅氏算法和序分量算法求取发电机机端电压各序分量。
[0123] 装置采样频率调节量按如下公式确定。
[0124]
[0125] 装置调整采样频率，使当前采样参考频率与电气量当前频率偏差的绝对值不大于固定门槛：
[0126] |ΔfS|≤ΔfS·set
[0127] 装置采样频率调节后，采用傅氏算法求取A、B、C三相电压的实部和虚部。
[0128]
[0129] 然后采用前述序分量计算公式求取发电机机端电压正、负、零序分量。
[0130] 以上实施例仅为说明本发明的技术思想，不能以此限定本发明的保护范围，凡是按照本发明提出的技术思想，在技术方案
基础上所做的任何改动，均落入本发明保护范围之内。
