[发明专利]用于机组阻尼在线分析的实时数据筛选方法

权利要求书

1.用于机组阻尼在线分析的实时数据筛选方法，其特征在于：包括以下步骤：

S01）、采样模块对机组励磁电压U_f、有功功率P、汽机主控指令U_k进行实时采样，将采样数据送入第一判断模块；

S02）、第一判断模块根据采样的励磁电压U_f和有功功率P实时监测有功功率是否出现扰动，当出现扰动时进入第二判断模块；

S03）、第二判断模块根据励磁电压U_f和有功功率P判断有功功率扰动后出现的振荡是否为衰减振荡，若为衰减振荡则将数据送入第三判断模块；

S04）、第三判断模块根据有功功率P和汽机主控指令U_k的变化方向判断该有功功率扰动是否为电气扰动，若为电气扰动，则认为此次扰动引起的功率振荡是有效的，可以用于在线阻尼分析。

2.根据权利要求1所述的用于机组阻尼在线分析的实时数据筛选方法，其特征在于：步骤S04具体为：

S41）、计算有功功率开始变化的时刻，t1时刻的有功功率值记作P(t1)，从t1时刻往前取10个数，包含t1时刻往后取Δt1*f个数,取的数记为p(1)、p(2)…p(Δt1*f+10)，其中f为采样频率，从t1开始的连续时间段记为Δt1；计算相邻数据的差值pd(n)，pd(n)=p(n)-p(n-1)，n=[2，Δt1*f+10]，对pd(11)到pd(Δt1*f+10)之间的数据pd(n)依次进行以下判断：M=max(pd(n-9),pd(n-8)…pd(n-1))，若满足[pd(n)1.5M pd(n+1)1.5M]，则将n的值计入数组a,读取a(1)的值，记为N，N即为功率开始变化的点，对应的时刻t_N= t1+（N-11）*1/f；

S42）、计算功率振荡的周期，从t_N时刻后选用2048个数据，记作p(m),m=1,2…2048，对P(m)数组进行傅里叶变换，识别幅值最大的非零频分量频率值，将此频率记作f_s，对应的周期为T_s=1/f_s；

S43）、寻找有功功率极值点，从t_N时刻判断P(x)是否大于（x-T_s/4，x+T_s/4）时间段内所有的有功功率值，若满足则判断P(x)为极大值点，记录时刻为t_maxi，功率为p_maxi；判断P(x)是否小于（x-T_s/4，x+T_s/4）时间段内所有的有功功率值，若满足则判断P(x)为极小值点,记录时刻为t_mini,功率为p_mini；x等于t_N、t_N+1/f、t_N+2/f……,当找到极值点后x的取值从t_maxi或t_mini直接往后跳过T_s/4继续寻找极值；

S44）、判断有功功率的首次变化方向，依据步骤S43寻找到的有功功率极值点，计算相邻有功功率极值点之间的变化量d1、d2…,dL，确定d1、d2…,dL中的最大值d_max，并将0.2d_max与d_max之前的极值点之间的变化量作比较，如果0.2d_max小于d_max之前的极值点之间的变化量，取d_max更前一级的极值点之间的变化量与0.2d_max比较，直至0.2d_max不小于与其比较的变化量，判定最后一次大于0.2d_max的变化量为初定的首次功率变化方向；将最后一次大于0.2d_max的变化量与第一次不大于0.2d_max的变化量比较，如果最后一次大于0.2d_max的变化量大于第一次不大于0.2d_max的变化量的两倍，则判定最后一次大于0.2d_max的变化量为首次功率变化方向，变化量向上为正，向下为负；

S45）、判断汽机主控指令在有功功率极值点附近的变化趋势，找到有功功率首次变化的时刻极值点对应的主控指令数据值U_K（n），判断[U_K（n-2）, U_K（n+2）]之间的数据是否存在大于或小于所有[U_K（n-3）, U_K（n+3）]之间其它数据的值，如果有最大值，汽机主控变化为正，有最小值，汽机主控变化为负；

S46）、判断有功功率变化方向和汽机主控变化方向是否相反，若方向相反，认为该有功功率振荡是电气扰动引起的，可以用于机组阻尼的在线分析，若方向相同，判断为无效的阻尼分析数据。