一种基于负阻尼贡献的风电场次同步振荡切除方法技术

一种基于负阻尼贡献的风电场次同步振荡切除方法，该方法通过采集风电场各条35kV集电线以及升压站出线三相相电压和三相相电流获取时域序列，然后对时域序列通过选择性滤波，滤除工频分量、高频分量，获取其中次同步/超同步振荡分量序列，对次同步/超同步振荡分量序列进行特征信息计算，得到次同步/超同步等效阻抗，接着对每一条集电线的阻抗按照电阻特性进行排序，按照排序切除负阻尼贡献最大的集电线，直至各条集电线都提供正电阻为止。该方法能够根据各条集电线对次同步/超同步振荡贡献的大小，精细化有步骤地切除产生振荡的风电机组，适合对风电机组大规模并网产生次同步/超同步振荡时进行有效防控。

【技术实现步骤摘要】
一种基于负阻尼贡献的风电场次同步振荡切除方法
本专利技术属于电力系统监测与控制
，具体涉及一种基于负阻尼贡献的风电场次同步振荡切除方法。
技术介绍
与传统发电机不同，风、光等新能源发电均通过电力电子装置并入电网，其多时间尺度的控制特性与电网自身特征相互作用，将可能引发次同步到谐波频段内，从而带来传统电力系统中没有的控制不稳定和振荡问题。近年来，这类问题已经在电网中逐步凸显，且波及范围越来越广，造成的后果也越来越严重。2014年德国北部海上风电场经直流输电送出系统发生200多Hz谐波振荡，谐波电流达到基波的40％以上，导致高压直流整流器的滤波电容爆炸，造成整个风场关断10个月之久，震动了整个风电和高压直流产业界。2011年以来，我国河北沽原地区风电场发生了上百次由风电机群与串补电网相互作用而引发的次同步谐振，其频率在3～10Hz内变化，曾造成变压器异常振动和大量风机脱网。2015年我国新疆哈密风电汇集地区多次发生电网次/超同步振荡现象，甚至引发距离200km以外的多台高压直流配套火电机组扭振保护动作跳机。事后事故分析表明，造成此次次同步振荡的原因是由于大量新能源汇集在电网中引入了大量的次同步谐波，该谐波频率与发电机轴系的自然振荡频率互补时，导致电气-机械扭振互作用现象发生，即次同步振荡。更重要的是，该次同步振荡的形成机理与传统电力系统不一致，是由于大量新能源汇集引入的次同步谐波造成的，已有的监测、保护控制方法已经不再适用。本专利技术针对大规模风电场汇集引入次同步谐波问题，通过对风电场升压站35kV集电线的负电阻特性提取，精准、有序地切除产生次同步振荡的风电机组，有效隔离激励源，保障系统安全稳定运行。
技术实现思路
本专利技术的目的是提供一种基于负阻尼贡献的风电场次同步振荡切除方法，能够准确有序地切除产生次同步振荡的风电机组，快速平息风电场的次同步振荡。一种基于负阻尼贡献的风电场次同步振荡切除方法，该方法主要应用于风电场升压站，通过对集电线的负电阻进行特征排序，来选择要切除的集电线路。首先，采集各集电线以及升压站出线的三相电压和电流，计算升压站出线瞬时功率并以此判断是否发生振荡，若振荡发生，根据振荡频率滤波得到集电线的谐振电压、谐振电流、谐振阻抗；其次，判断谐振阻抗是否呈现负电阻特性，若存在负电阻则选取负电阻最大的集电线路切除，根据切除前后出线瞬时功率的衰减系数变化率，决定后续需要切除的集电线路；再次，若剩余集电线路均不具有负电阻特性而次同步振荡仍持续，则选择次同步振荡电流幅值最大的集电线路进行切除；最后，切除后延时稳定一段时间，重复以上步骤继续选择切除线路，直至次同步振荡消失或线路全部切除。具体方案包括以下步骤：(1)以采样频率fs采集第i条集电线(升压站进线)三相相电压uai、ubi、uci，三相相电流iai、ibi、ici，获取三相相电压uai、ubi、uci，三相相电流iai、ibi、ici的时间序列uai(k)、ubi(k)、uci(k)和iai(k)、ibi(k)、ici(k)，同时以采样频率fs采集风场升压站出线的三相相电压uaoj、uboj、ucoj，三相相电流iaoj、iboj、icoj，获取对应时间序列uaoj(k)、uboj(k)、ucoj(k)和iaoj(k)、iboj(k)、icoj(k)，计算风场升压站出线的瞬时功率pj(k)；(2)对pj(k)序列截取数据长度为N的有限长度数据窗，对数据窗内数据进行FFT计算得到频谱数据，寻找pj(k)频谱的极值点Pj(k1)并且Pj(k1)≥TH，即可认为存在振荡频率为的次同步/超同步振荡，其中，TH为阈值，同时记录瞬时功率数据pj(k)，寻找振荡起始阶段瞬时功率的振荡包络线，并求取衰减系数dj，其中j表示次同步振荡切除线路的轮次，j＝0，1，…，n；(3)对满足条件的频点频率信号进行滤波提取，得到该频率下的实时序列和(4)在时间窗T内通过比较实时序列，分别获取序列的峰值以及同相电压电流之间的相角得到三相相阻抗的计算方式；以及(5)分别计算的实部，得到从集电线看到的风电机组外部电阻特性，(6)在风电场升压站计算集电线次同步频率下的相阻抗，实部分别为：对于各条集电线的按照从小到大进行排序；(7)如果存在实部为负，则第一轮切除操作先切除各条集电线的中电阻最小的一条集电线，最小值为Rjmin，同时记录风场出线的瞬时功率为pj+1(k)，计算衰减系数dj+1；(8)计算Δdj＝dj+1-dj，获取关联系数cj＝Δdj/Rjmin；(9)去除切除掉的线路后，延长固定时间T1，重复步骤(2)～(6)，重新生成按照从小到大排列的各条线路阻抗实部集合，如果集合中存在负电阻，则从集合中从小到大选取负电组的线路，使得选择的线路负电阻之和R(j+1)min满足dj+1≤cj|R(j+1)min|，切除选择的负电阻特性线路，令j＝j+1，记录风场出线的瞬时功率为pj+1(k)，计算衰减系数dj+1；(10)重复步骤(8)～(9)，直至不存在负电阻特性线路为止；(11)判断计算出来的瞬时功率pj+1(k)，如果持续振荡，切除集电线中振荡电流幅值最大的线路；(12)去除切除掉的线路后，延长固定时间T1，重复步骤(11)，直至振荡消除。其中，在上述步骤(2)中：寻找pj(k)频谱的极值点PPj(k1)并且PPj(k1)≥TH，其中TH以集电线额定功率和的0.5％来选取，对于瞬时功率的包络线的选择采用滑动数据窗内pj(k)最大值点来选择，滑动数据窗的时间宽度为100ms，以500ms内最大值点来进行衰减系数的计算，衰减系数的计算方式如下：(1)对于500ms内pj(k)的最大值序列为pjmax(n)，分别求取对数，即lnpjmax(n)；(2)对于对数序列数组[lnpjmax(n)，tn]，采用最小二乘法拟合，求取线性表达式y＝ax+b中的a，即衰减系数，其中，y对应的数据序列为lnpjmax(n)，x对应的数据序列为tn。其中，在上述步骤(4)中：在时间窗T内通过比较实时序列，分别获取序列的峰值以及同相电压电流之间的相角时间窗T通常取1s，电压和电流的峰值通过判断极值点的方法获得，同相电压电流之间的相角计算方法如下：(1)以A相为例进行说明，B、C相采用同样的计算过程，相角的计算基于DFT算法；(2)对1s数据窗内的数据序列进行离散傅里叶变换，得到电流序列在频率的实部和虚部具体计算公式如下：其中，N为数据窗长度；(3)根据实部和虚部，得到相角把相角进行360度取余，使得相角在0～360°范围之间；(4)对电压序列重复(2)～(3)，得到电压序列的相角(5)A相电压电流之间的相角单位为度。本专利技术具有以下有益效果：本专利技术提出的一种基于负阻尼贡献的风电场次同步振荡切除方法，相对于现阶段大规模风电场并网引发次同步振荡常规依据振荡幅值进行线路切除来说，基于负阻尼贡献的切除方法能够精确有序的切除振荡线路，快速隔离激励源，为风电场并网引发次同步振荡的保护与控制方案提供了一种解决方式。附图说明图1为本专利技术方法的实现流程图；图2为瞬时功率衰减系数拟合计算示意图；图3为电压电流相角示意图；图4为从原始采样数据中分离出来的电压电流次同步(超同步)振荡分量以及相角及电阻图。具体实施方式下面本文档来自技高网...

【技术保护点】
1.一种基于负阻尼贡献的风电场次同步振荡切除方法，其特征在于，该方法包括以下步骤：首先，采集各集电线以及升压站出线的三相电压和电流，计算升压站出线瞬时功率并以此判断是否发生振荡，若振荡发生，根据振荡频率滤波得到集电线的谐振电压、谐振电流、谐振阻抗；其次，判断谐振阻抗是否呈现负电阻特性，若存在负电阻则选取负电阻最大的集电线路切除，根据切除前后出线瞬时功率的衰减系数变化率，决定后续需要切除的集电线路；再次，若剩余集电线路均不具有负电阻特性而次同步振荡仍持续，则选择次同步振荡电流幅值最大的集电线路进行切除；最后，切除后延时稳定一段时间，重复以上步骤继续选择切除线路，直至次同步振荡消失或线路全部切除。

【技术特征摘要】
1.一种基于负阻尼贡献的风电场次同步振荡切除方法，其特征在于，该方法包括以下步骤：首先，采集各集电线以及升压站出线的三相电压和电流，计算升压站出线瞬时功率并以此判断是否发生振荡，若振荡发生，根据振荡频率滤波得到集电线的谐振电压、谐振电流、谐振阻抗；其次，判断谐振阻抗是否呈现负电阻特性，若存在负电阻则选取负电阻最大的集电线路切除，根据切除前后出线瞬时功率的衰减系数变化率，决定后续需要切除的集电线路；再次，若剩余集电线路均不具有负电阻特性而次同步振荡仍持续，则选择次同步振荡电流幅值最大的集电线路进行切除；最后，切除后延时稳定一段时间，重复以上步骤继续选择切除线路，直至次同步振荡消失或线路全部切除。2.一种基于负阻尼贡献的风电场次同步振荡切除方法，其特征在于，所述方法具体步骤如下：(1)以采样频率fs采集第i条集电线(升压站进线)三相相电压uai、ubi、uci，三相相电流iai、ibi、ici，获取三相相电压uai、ubi、uci，三相相电流iai、ibi、ici的时间序列uai(k)、ubi(k)、uci(k)和iai(k)、ibi(k)、ici(k)，同时以采样频率fs采集风场升压站出线的三相相电压uaoj、uboj、ucoj，三相相电流iaoj、iboj、icoj，获取对应时间序列uaoj(k)、uboj(k)、ucoj(k)和iaoj(k)、iboj(k)、icoj(k)，计算风场升压站出线的瞬时功率pj(k)；(2)对pj(k)序列截取数据长度为N的有限长度数据窗，对数据窗内数据进行FFT计算得到频谱数据，寻找pj(k)频谱的极值点Pj(k1)并且Pj(k1)≥TH，即可认为存在振荡频率为fk1＝k1fs/N的次同步/超同步振荡，其中，TH为阈值，同时记录瞬时功率数据pj(k)，寻找振荡起始阶段瞬时功率的振荡包络线，并求取衰减系数dj，其中j表示次同步振荡切除线路的轮次，j＝0,1，…,n；(3)对满足条件的频点频率信号进行滤波提取，得到该频率下的实时序列和(4)在时间窗T内通过比较实时序列，分别获取序列的峰值以及同相电压电流之间的相角得到三相相阻抗的计算方式；以及(5)分别计算的实部，得到从集电线看到的风电机组外部电阻特性，(6)在风电场升压站计算集电线次同步频率下的相阻抗，实部分别为：对于各条集电线的按照从小到大进行排序；(7)如果存在实部为负，则第一轮切除操作先切除各条集电线的或者或者中电阻最小的一条集电线，最小值为Rjmin，同时记录风场出线的瞬时功率为pj+1(k)，计算衰减系数dj+1；(8)计算△dj＝dj+1-dj，获取关联系数cj＝△dj/Rjmin；(9)去除切除掉的线路后，延长固定时间T1，重复步骤(2)～(6)，重新生成按照从小到大排列的...

【专利技术属性】
技术研发人员：王春华，宋晓喆，刘志君，郭雷，赵伟华，高培生，孙福寿，常富杰，孙铭徽，曲振军，姜旭，刘家岩，丁浩，
申请(专利权)人：国网吉林省电力有限公司，北京四方继保自动化股份有限公司，
类型：发明
国别省市：吉林,22