CN101013809A - 高速采样的间断角闭锁差动保护数值方法 - Google Patents

Abstract

高速采样的间断角闭锁差动保护数值方法，采样速率大于或等于2400Hz，步骤如下：交流采样、含角差转换的差流计算、计算差流微分、对差流微分值进行滤波并取绝对值、浮动门槛计算、波宽和间断角计算、小扰动保持措施、补偿角差值计算、差动判断。本发明从理论上给出对应数值算法的动作门槛，利用涌流的物理特征，直接、快速、正确地区分涌流与故障电流，达到在变压器空投、内部故障、外部故障切除，或涌流情况下能快速动作，提高了变压器保护在选择性、可靠性和灵敏度等方面的性能，同时增加了小扰动保持措施，使本方法可以在小电流信号情况下有效防止采样噪声和数值算法舍入误差影响，进一步保证了差动保护可靠动作。

Description

高速采样的间断角闭锁差动保护数值方法

技术领域

本发明涉及一种变压器继电保护方法，特别是一种高速采样的间断角闭锁差动保护数值方法。

背景技术

变压器继电保护装置是变压器和电网安全运行的重要组成部分，使用谐波闭锁或波形鉴别闭锁励磁涌流的方法是变压器数字保护常用的方法。

目前上述方法存在的问题有：

1)由于现代变压器工作磁通密度的提高，励磁涌流中二次谐波分量减少，使得二次谐波闭锁的灵敏性、选择性和可靠性降低。

2)在含有串补电容或电力电缆的高压系统中的变压器如果发生内部故障，谐波含量会超过15％，引起差动拒动。

3)采样速率较低的间断角数值方法不能满足间断角原理上的精度要求。

4)提高采样速率来判别间断角的方法，将使计算机负荷增加量巨大，而目前高压变压器保护一般采用主后备一体化，因此很难实时完成高速间断角计算。

5)变压器轻微故障时，由于故障信号幅值较低，由一般的间断角数值算法难免受到采样噪声和计算过程中舍入误差的影响，会造成保护拒动。

6)没有精确的理论推导来得到在不同间断角数值算法情况下的动作门槛。

间断角原理的变压器差动保护率先由我国于60年代提出，间断角是励磁涌流固有的物理特征，相对于抽象的谐波制动原理，具有物理意义明确的特点，据此原理实现的模拟式保护装置至今已得到广泛应用，参见王祖光.间断角原理的变压器差动保护.电力系统自动化，1979，3(1)。随着微机保护技术的发展也有一些尝试将间断角闭锁原理应用于微机变压器保护，参见朱亚明，等.间断角原理的变压器差动保护的性能特点及微机实现.电力系统自动化，1996，20(11)，并提出一种采样频率较低的间断角闭锁数值方法，其计算过程是将采样值进行插值后，与浮动门槛比较并计数，累加得到间断角和波宽，在采样频率为1800Hz时，其固有的误差为10°，因此其判别效果具有天生的局限性。为了使间断角数值算法实用化，一个基本方法是提高采样频率，但采样频率的提高，首先增加了运算量，另外最重要的是采样频率提高的同时增加了高频采样噪声的负面影响，因此，并非如其他学者预测的那样，简单的提高采样频率就能够提高间断角判别的准确性。因此对于间断角的应用仍然停留在理论分析的层次，很多研究成果都对间断角原理数值实现的复杂性做出了正确的判断。参见赵永彬，等.多原理的数字式变压器保护的应用.电力系统自动化学报，2004，16(5)；唐跃中，等.几种变压器励磁涌流判别方法的特点及其内在联系的分析.电力系统自动化，1995，19(9)。

发明内容

本发明所要解决的技术问题是，克服现有技术的缺点，提供一种高速采样的间断角闭锁差动保护数值方法。

本发明解决其技术问题所采用的技术方案如下：一种高速采样的间断角闭锁差动保护数值方法，采样速率大于或等于2400Hz，对应采样频率每周波至少48点，其包括如下步骤：

(1)、交流采样：获得变压器各电压等级侧的电流采样；

(2)、含角差转换的差流计算：n采样点的差电流 $i_n=\underset{j=1}{\overset{N}{\mathrm{\Σ}}}i(j,n),$ 其中，N为变压器所有分支的个数，i(j，n)表示n采样点经过角差转换的第j个变压器分支的采样电流；

(3)、计算差流微分：n采样点的差电流微分值id_n＝i_n-i_n-1；

(4)、对差流微分值进行滤波并取绝对值： ${{\mathrm{id}}_n}^{\′}=\left|\underset{i=0}{\overset{3}{\mathrm{\Σ}}}{\mathrm{id}}_{n-1}\right|;$

(5)、浮动门槛计算：浮动门槛i_TH(n)＝max(λ·max(id_m′)，i_THF)，m∈[1，n]，且m为整数，其中λ：固定系数，取值范围为[0.1，0.4]，i_THF：固定门槛的标么值，其取值范围为[0.1，0.3]，每个id_n′对应一浮动门槛i_TH(n)；

(6)、波宽和间断角计算：一个采样波内，设每周波的采样点数为N_sam，连续的id_n′值大于对应浮动门槛i_TH(n)的点数分别计为a₁，a₂，...，a_i，与其对应的各个波宽分别为

连续的id_n′值小于对应浮动门槛i_TH(n)的点数分别计为b₁，b₂，...，b_j与其对应的各个间断角分别为

上述式子中i、j均为正整数；

(7)小扰动保持措施：设置一抗扰动角度值θ_min，θ_min的取值范围为[5°，15°]，设连续的id_n′值小于对应浮动门槛i_TH(n)的点的个数为c，若c满足不等式(I)：

则将这c个点视为小扰动，并将这c个点以及这c个点的后一波宽的数值计入这c个点的前一波宽的数值中去；设连续的id_n′值大于对应浮动门槛i_TH(n)的点的个数为d，若d满足不等式(II)：则这d个点视为小扰动，并将这d个点以及d个点的后一间断角的数值计入前一间断角的数值中去，上述式子中，N_sam为每周波采样点数；从波形的起始位置开始，在满足不等式(I)、(II)的采样点依次交替出现的情况下，在上述的采样点范围内，如果不等式(I)先成立，则将上述满足不等式(I)、(II)的采样点全部计入前一波宽的数值中去，否则将这些采样点全部计入前一间断角的数值中去；

(8)、补偿角差值计算：相邻的两个采样点id_n-1′、id_n′，分别满足不等式id_n-1′＞i_TH(n-1)和id_n′＜i_TH(n)，并且它们分别被计入波宽和间断角的数值中，则通过插值计算得当前波宽补偿角 ${\mathrm{\θ}}_{\mathrm{comB}}={\mathrm{\θ}}_{\mathrm{sam}}\·\frac{|{{\mathrm{id}}_{n-1}}^{\′}-i_{\mathrm{TH}(n-1)}|}{{{|\mathrm{id}}_{n-1}}^{\′}-{{\mathrm{id}}_n}^{\′}|},$ 将θ_comB加入当前波宽的数值中去；如果id_n-1′、id_n′分别满足不等式id_n-1′＜i_TH(n-1)和id_n′＞i_TH(n)，并且它们分别被计入间断角和波宽的数值中，通过插值计算得当前间断角补偿角 ${\mathrm{\θ}}_{\mathrm{comJ}}={\mathrm{\θ}}_{\mathrm{sam}}\·\frac{{{|\mathrm{id}}_{n-1}}^{\′}-i_{\mathrm{TH}(n-1)}|}{|{{\mathrm{id}}_{n-1}}^{\′}-{{\mathrm{id}}_n}^{\′}|},$ 将θ_comJ加入当前间断角的数值中去，上述等式中，

N_sam为每周波采样点数；

(9)、差动判断：判断一周波内的最大波宽值是否大于115°，如果不大于115°，则闭锁差动，如果大于115°，则判断该周波内的最大间断角是否小于70°，如果该间断角不小于70°，则闭锁差动，如果小于70°，则开放差动。

本发明的有益效果如下：本发明提供了在变压器差动保护领域内使用高速采样间断角闭锁原理的数值方法，其对差流进行微分计算，滤除了直流分量对差电流的影响，保证了间断角计算的正确性，对差流微分值进行4点的平滑滤波，去除高频信号的干扰，从而可以精确计算波宽和间断角，并从理论上给出对应数值算法的动作门槛，以确定开放或闭锁差动保护；该方法利用涌流的物理特征，直接、快速、正确地区分涌流与故障电流，达到在变压器空投、内部故障、外部故障切除，或涌流情况下能快速动作，可靠闭锁差动保护，提高了变压器保护在选择性、可靠性和灵敏度等方面的性能，同时增加了小扰动保持措施，使本方法可以在小电流信号的情况下有效防止采样噪声和数值算法舍入误差的影响，进一步保证了差动保护可靠动作。

附图说明

图1为本发明高速采样的间断角闭锁差动保护数值方法的流程框图。

图2为利用计算机仿真软件按照采样频率为每周波48点采样直接生成差电流微分(滤波后的)波形图。

图3为对图2中的滤波后的差流微分值(纵坐标值)取绝对值后形成的波形图。

具体实施方式

下面参照附图并结合实施例对本发明作进一步详细描述。但是本发明不限于所给出的例子。

本发明高速采样的间断角闭锁差动保护数值方法，采样速率大于或等于2400Hz，对应采样频率每周波至少48点，其包括如下步骤：

(1)、交流采样：获得变压器各电压等级侧的电流采样；

(2)、含角差转换的差流计算：n采样点的差电流 $i_n=\underset{j=1}{\overset{N}{\mathrm{\Σ}}}i(j,n),$ 其中，N为变压器所有分支的个数，i(j，n)表示n采样点经过角差转换的第j个变压器分支的采样电流；

(3)、计算差流微分：n采样点的差电流微分值id_n＝i_n-i_n-1；

(4)、对差流微分值进行滤波并取绝对值： ${{\mathrm{id}}_n}^{\′}=\left|\underset{i=0}{\overset{3}{\mathrm{\Σ}}}{\mathrm{id}}_{n-i}\right|;$

其中，步骤(1)至步骤(4)为目前通常采用的技术手段，无需进行详细描述，具体计算公式也已经给出。

发明人利用计算机仿真软件按照采样频率为每周波48点采样直接生成差电流微分(滤波后的)波形，并选取该波形中

个波形(60个连续的采样点)，前 个波形(前12个采样点)为正常稳态状态下的波形，具有一定的抖动，但是在允许的范围之内，后1个波形(后48个采样点)为故障波形，是本是实施例的实验对象，前

个波形引入的主要目的在于直观的比较稳态波形与故障波形的差别，如图2所示为

个波形的差电流微分(滤波后的)波形图，图3为对图2中滤波后的差流微分值取绝对值后得到的波形图，图3中，曲线f为浮动门槛曲线，对应的数据参见表1。

步骤(5)、浮动门槛计算，i_TH(n)＝max(λ·max(id_m′)，i_THF)，m∈[1，n]，且m为整数，其中λ：固定系数，取值0.26，i_THF：固定门槛的标么值，其取值为0.2，每个id_n′对应一浮动门槛i_TH(n)，见表1中，以第9个采样点为例，m∈[1，9]，且m为整数，max(id_m′)＝0.1571，i_TH(9)＝max(0.26·0.1571，0.2)＝max(0.0408，0.2)＝0.2；在第19个采样点为例，m∈[1，19]，且m为整数，max(id_m′)＝0.8388，i_TH(19)＝max(0.26·0.8388，0.2)＝max(0.2181，0.2)＝0.2181；再在以表1中第29个采样点为例，m∈[1，29]，且m为整数，max(id_m′)＝1.1180，i_TH(19)＝max(0.26·1.1180，0.2)＝max(0.2907，0.2)＝0.2907。

步骤(6)、波宽和间断角计算：见表1，一个采样波内，每周波的采样点数为48，连续的id_n′值大于对应浮动门槛i_TH(n)的点数分别计为20、1、1、19，与其对应的各个波宽分别为：

即分别为：150°，7.5°，7.5°，142.5°，连续的id_n′值小于对应浮动门槛i_TH(n)的点数分别计为1、1、3、1、1与其对应的各个间断角为

即分别为：7.5°、7.5°、22.5°、7.5°、7.5°。

步骤(7)、小扰动保持措施，设置一抗扰动角度值θ_min，其取值范围应为[5°，15°]，本实施例中θ_min的取值为10°，设连续的id_n′值小于对应浮动门槛i_TH(n)的点的个数为c，若c满足不等式(I)：

即c＜1.67，则这c个点视为小扰动，并将这c个点以及这c个点的后一波宽的数值计入这c个点的前一波宽的数值中去，如表1所示的实施例中，第34个采样点满足不等式(I)，将其视为小扰动，并将第34个点以及该采样点后一波宽的数值(7.5°)计入第34个采样点的前一波宽的数值(150°)中去，这样第一个波宽的值为

设连续的id_n′值大于对应浮动门槛i_TH(n)的点的个数为d，若d满足不等式(II)： 即d＜1.67，则这d个点视为小扰动，并将这d个点以及d个点的后一间断角的数值计入前一间断角的数值中去，如表1所示的实施例中，第39个采样点满足不等式(II)，将其视为小扰动，并将第39个点以及该采样点后一间断角的数值(7.5°)计入第39个采样点的前一间断角的数值(22.5°)中去，这样第三个间断角的值为

从波形的起始位置开始，在满足不等式(I)、(II)的采样点依次交替出现的情况下，在上述的采样点范围内，如果不等式(I)先成立，则将上述满足不等式(I)、(II)的采样点全部计入前一波宽的数值中去，否则将这些采样点全部计入前一间断角的数值中去；本实施例中，如果将θ_min取值为20。(实际操作中θ_min的取值范围应为[5°，15°]，这里将θ_min的取值为20°仅仅用以说明本发明中小扰动保持的方法)，则采样点34，36、37、38，40满足不等式(I)，采样点35，39满足不等式(II)，满足不等式(I)、(II)的采样点依次交替出现，即连续出现小抖动，此时不等式(I)首先被满足，因此这些连续抖动的满足不等式(I)、(II)的采样点全部计入前一波宽(第1个波宽)的数值中去，计算得

同样的如果首先被满足的不是不等式(I)，那么这些点就应计入间断角的数值中去。采用小扰动措施是考虑到故障电流和差电流不能突变，在小电流信号下由于判别门槛相应地降低，使得采样噪声的影响变得更加显著，但差电流的基本变化趋势不变，在噪声的影响下计算的波宽和间断角如果以小角度交替出现，那么将较小角度的波宽或间断保持计入上一段计算的间断或波宽。

步骤(8)、补偿角差值计算：相邻的两个采样点id_n-1′、id_n′，分别满足不等式id_n-1′＞i_TH(n-1)和id_n′＜i_TH(n)，即它们位于各自对应浮动门槛的上、下，并且它们分别被计入波宽和间断角的数值中，也就是说在这两个点不是小抖动，那么通过插值计算得波宽补偿角 ${\mathrm{\θ}}_{\mathrm{comB}}={\mathrm{\θ}}_{\mathrm{sam}}\·\frac{{|{\mathrm{id}}_{n-1}}^{\′}-i_{\mathrm{TH}(n-1)}|}{|{{\mathrm{id}}_{n-1}}^{\′}-{{\mathrm{id}}_n}^{\′}|}$ (上述等式中，

N_sam为每周波采样点数)，并将θ_comB加入当前波宽的数值中去，从表1可知，采样点35、36满足上述的条件，采样点59、60也满足上述条件，以采样点35、36为例说明，  id₃₅′＝0.3637，i_TH(35)＝0.2907，id₃₆′＝0.0722，i_TH(36)＝0.2907，因此id₃₅′＞i_TH(35)，id₃₆′＜i_TH(36)，并且采样35被计入波宽，采样点36被计入间断角，因此通过插值计算得波宽补偿角

将其加入当前波宽的数值中去，即得当前的波宽(第一个波宽)＝165°+1.878°＝166.878°；同理，根据采样点59、60算得的波宽补偿角算 ${\mathrm{\θ}}_{\mathrm{comB}}={\mathrm{\θ}}_{\mathrm{sam}}\·\frac{|{{\mathrm{id}}_{59}}^{\′}-i_{\mathrm{TH}\left(59\right)}|}{|{{\mathrm{id}}_{59}}^{\′}-{{\mathrm{id}}_{60}}^{\′}|}$

得第4个波宽值＝142.5°+2.742°＝145.242°；如果id_n-1′、id_n′分别满足不等式id_n-1′＜i_TH(n-1)和id_n′＞i_TH(n)，并且它们分别被计入间断角和波宽的数值中，通过插值计算得间断角补偿角 ${\mathrm{\θ}}_{\mathrm{comJ}}={\mathrm{\θ}}_{\mathrm{sam}}\·\frac{|{{\mathrm{id}}_{n-1}}^{\′}-i_{\mathrm{TH}(n-1)}|}{|{{\mathrm{id}}_{n-1}}^{\′}-{{\mathrm{id}}_n}^{\′}|}$ (上述等式中，

N_sam为每周波采样点数)，将θ_comJ加入间断角的数值中去；从表1可知，采样点13、14满足上述的条件，采样点40、41也满足上述条件，以采样点40、41为例说明，id₄₀′＝0.1867，i_TH(40)＝0.2907，id₄₁′＝0.3548，i_TH(41)＝0.2907，因此id₄₀′＜i_TH(40)，id₄₁′＞i_TH(41)，并且采样40被计入间断角，采样点41被计入波宽，因此通过插值计算得波宽补偿角

将其加入当前间断角的数值中去，即得当前的间断角(第三个间断角)＝37.5°+4.640°＝42.140°；同理，根据采样点13、14算得的间断角补偿角算

得第1个间断角＝7.5°+1.253°＝8.753°；

(9)、差动判断：判断一周波内的最大波宽值是否大于115°，如果不大于115°，则闭锁差动，如果大于115°，则判断该周波内的最大间断角是否小于70°，如果该间断角不小于70°，则闭锁差动，如果小于70°，则开放差动，本实施例中的最大波宽为166.878°，最大间断角为42.140°，最大波宽大于115°，则比较最大间断角是否小于70°，最大间断角42.140°＜70°，则开放差动。

表1

采样点	差流微分值	取绝对值后的差流微分值	浮动门槛值	波宽	间断	补偿角	小扰动	备注
									1	0.1516	0.1516	0.2000	0	1	0
2	0.0197	0.0197	0.2000	0	1	0
									3	-0.0864	0.0864	0.2000	0	1	0
4	-0.1571	0.1571	0.2000	0	1	0
									5	0.0908	0.0908	0.2000	0	1	0
6	0.0981	0.0981	0.2000	0	1	0
									7	-0.1462	0.1462	0.2000	0	1	0
8	0.1071	0.1071	0.2000	0	1	0
									9	-0.1219	0.1219	0.2000	0	1	0
10	-0.1634	0.1634	0.2000	0	1	0
									11	0.1072	0.1072	0.2000	0	1	0
12	0.1043	0.1043	0.2000	0	1	0
									13	0.1801	0.1801	0.2000	0	1	1.253		故障开始时刻
14	0.2992	0.2992	0.2000	1	0	0
									15	0.2025	0.2025	0.2000	1	0	0
16	0.5063	0.5063	0.2000	1	0	0
									17	0.3625	0.3625	0.2000	1	0	0
18	0.7392	0.7392	0.2000	1	0	0
									19	0.8388	0.8388	0.2181	1	0	0
20	0.5956	0.5956	0.2181	1	0	0
									21	0.8901	0.8901	0.2314	1	0	0
22	0.9879	0.9879	0.2568	1	0	0
									23	1.0402	1.0402	0.2705	1	0	0
24	1.0454	1.0454	0.2718	1	0	0
									25	1.0427	1.0427	0.2718	1	0	0
26	0.9548	0.9548	0.2718	1	0	0
									27	1.1180	1.1180	0.2907	1	0	0
28	1.1086	1.1086	0.2907	1	0	0
									29	0.8826	0.8826	0.2907	1	0	0
30	0.6454	0.6454	0.2907	1	0	0
									31	0.5334	0.5334	0.2907	1	0	0
32	0.7030	0.7030	0.2907	1	0	0

33	0.6359	0.6359	0.2907	1	0	0
									34	0.2335	0.2335	0.2907	0	1	0	√
35	0.3637	0.3637	0.2907	1	0	1.878
									36	0.0722	0.0722	0.2907	0	1	0
37	-0.1868	0.1868	0.2907	0	1	0
									38	-0.1965	0.1965	0.2907	0	1
39	-0.4070	0.4070	0.2907	1	0		√
									40	-0.1867	0.1867	0.2907	0	1	4.640
41	-0.3548	0.3548	0.2907	1	0	0
									42	-0.7363	0.7363	0.2907	1	0	0
43	-0.6792	0.6792	0.2907	1	0	0
									44	-0.9495	0.9495	0.2907	1	0	0
45	-0.9600	0.9600	0.2907	1	0	0
									46	-0.9988	0.9988	0.2907	1	0	0
47	-0.9044	0.9044	0.2907	1	0	0
									48	-0.8579	0.8579	0.2907	1	0	0
49	-0.8622	0.8622	0.2907	1	0	0
									50	-0.8679	0.8679	0.2907	1	0	0
51	-0.9349	0.9349	0.2907	1	0	0
									52	-1.0663	1.0663	0.2907	1	0	0
53	-0.8621	0.8621	0.2907	1	0	0
									54	-0.9197	0.9197	0.2907	1	0	0
55	-0.6914	0.6914	0.2907	1	0	0
									56	-0.5917	0.5917	0.2907	1	0	0
57	-0.4803	0.4803	0.2907	1	0	0
									58	-0.3476	0.3476	0.2907	1	0	0
59	-0.4195	0.4195	0.2907	1	0	2.742
									60	0.0672	0.0672	0.2907	0	1	0

注：表1中的对于一个采样点，其“波宽”标记为“0”，“间断角”标记为“1”，则表示该采样点的“取绝对值后的差流微分值”小于对应的“浮动门槛值”；对于一个采样点，其“波宽”标记为“1”，“间断角”标记为“0”，则表示该采样点的“取绝对值后的差流微分值”大于对应的“浮动门槛值”；“小扰动”一栏中的“√”表示该点被判为小扰动。

本实施例仅以每周波48点的采样频率、对应采样速率为2400Hz为例说明本发明的方法，实际操作中，采样速率大于或等于2400Hz，采样频率每周波至少48点，即可满足，采样频率可以是每周波60点、对应采样速率为3000Hz或者100点、对应采样速率为5000Hz，其原理、方法与48点采样、对应采样速率为2400Hz的原理、方法一致。

Claims (3)

1、高速采样的间断角闭锁差动保护数值方法，采样速率大于或等于2400Hz，对应采样频率每周波至少48点，其包括如下步骤：

(1)、交流采样：获得变压器各电压等级侧的电流采样；

(2)、含角差转换的差流计算：n采样点的差电流 $i_n=\underset{j=1}{\overset{N}{\mathrm{\Σ}}}i(j,n),$ 其中，N为变压器所有分支的个数，i(j，n)表示n采样点经过角差转换的第j个变压器分支的采样电流；

(3)、计算差流微分：n采样点的差电流微分值id_n＝i_n-i_n-1；

(4)、对差流微分值进行滤波并取绝对值： ${\mathrm{id}}_n^{\′}=\left|\underset{i=0}{\overset{3}{\mathrm{\Σ}}}{\mathrm{id}}_{n-1}\right|;$

(5)、浮动门槛计算：浮动门槛i_TH(n)＝max(λ·max(id_m′)，i_THF)，m∈[1，n]，且m为整数，其中λ：固定系数，取值范围为[0.1，0.4]，i_THF：固定门槛的标么值，其取值范围为[0.1，0.3]，每个id_n′对应一浮动门槛i_TH(n)；

(6)、波宽和间断角计算：一个采样波内，设每周波的采样点数为N_sam，连续的id_n′值大于对应浮动门槛i_TH(n)的点数分别计为a₁，a₂，...，a_i，与其对应的各个波宽分别为 连续的id_n′值小于对应浮动门槛i_TH(n)的点数分别计为b₁，b₂，...，b_j与其对应的各个间断角分别为上述式子中i、j均为正整数；

(7)小扰动保持措施：设置一抗扰动角度值θ_min，θ_min的取值范围为[5°，15°]，设连续的id_n′值小于对应浮动门槛i_TH(n)的点的个数为c，若c满足不等式(I)：

则将这c个点视为小扰动，并将这c个点以及这c个点的后一波宽的数值计入这c个点的前一波宽的数值中去；设连续的id_n′值大于对应浮动门槛i_TH(n)的点的个数为d，若d满足不等式(II)：

则这d个点视为小扰动，并将这d个点以及d个点的后一间断角的数值计入前一间断角的数值中去，上述式子中，N_sam为每周波采样点数；从波形的起始位置开始，在满足不等式(I)、(II)的采样点依次交替出现的情况下，在上述的采样点范围内，如果不等式(I)先成立，则将上述满足不等式(I)、(II)的采样点全部计入前一波宽的数值中去，否则将这些采样点全部计入前一间断角的数值中去；

(8)、补偿角差值计算：相邻的两个采样点id_n-1′、id_n′，分别满足不等式id_n-1′＞i_TH(n-1)和id_n′＜i_TH(n)，并且它们分别被计入波宽和间断角的数值中，则通过插值计算得当前波宽补偿角 ${\mathrm{\θ}}_{\mathrm{comB}}={\mathrm{\θ}}_{\mathrm{sam}}\·\frac{{|{\mathrm{id}}_{n-1}}^{\′}-i_{\mathrm{TH}(n-1)}|}{|{{\mathrm{id}}_{n-1}}^{\′}-{{\mathrm{id}}_n}^{\′}|},$ 将θ_comB加入当前波宽的数值中去；如果i_dn-1′、id_n′分别满足不等式id_n-1′＜i_TH(n-1)和id_n′＞i_TH(n)，并且它们分别被计入间断角和波宽的数值中，通过插值计算得当前间断角补偿角 ${\mathrm{\θ}}_{\mathrm{comJ}}={\mathrm{\θ}}_{\mathrm{sam}}\·\frac{{|{\mathrm{id}}_{n-1}}^{\′}-i_{\mathrm{TH}(n-1)}|}{|{{\mathrm{id}}_{n-1}}^{\′}-{{\mathrm{id}}_n}^{\′}|},$ 将θ_comJ加入当前间断角的数值中去，上述等式中，

N_sam为每周波采样点数；

(9)、差动判断：判断一周波内的最大波宽值是否大于115°，如果不大于115°，则闭锁差动，如果大于115 °，则判断该周波内的最大间断角是否小于70°，如果该间断角不小于70°，则闭锁差动，如果小于70°，则开放差动。

2、根据权利要求1所述的高速采样的间断角闭锁差动保护数值方法，其特征是：其采样速率为3000Hz，对应采样频率为每周波60点。

3、根据权利要求1所述的高速采样的间断角闭锁差动保护数值方法，其特征是：其采样速率为5000Hz，对应采样频率为每周波100点。

Images