[发明专利]基于插值DFT的同步相量数据次同步振荡参数辨识法

权利要求书

1.基于插值DFT的同步相量数据次同步振荡参数辨识法，其特征在于，对相量测量单元获取的相量数据进行频谱分析，从而获得次同步振荡的参数，包括以下步骤：

步骤一、采用广域测量系统内的相量测量单元获取实时高分辨率时域波形数据；

步骤二、在相量测量单元中对时域波形数据进行傅里叶变换得到同步相量数据，采用上报频率f_r，即以间距f_pr＝f_p/f_r对同步相量数据进行重新采样，在主站获得上报的同步相量X_r(m)，f_p为相量测量单元对时域波形数据的采样频率；

步骤三、在定窗宽为Nr的矩形窗内对同步相量X_r(m)进行离散傅里叶变换操作，得到离散傅里叶变换谱F_r(k)；

步骤四、用次同步振荡相量的最大幅值谱确定次同步振荡幅度最大的频谱线k_m；

步骤五、构造指标相量[k1，k2，k3，k4]，引入二阶差分光谱的频谱F_rs(k_i)，获得二阶差分的比值R，利用泰勒级数展开的一阶项对F_rs(k_i)进行重写，利用且k_m≈L_rs，获得简化后的二阶差分比值R；

步骤六、获得次同步振荡分量的频率f_s和次同步振荡分量的阻尼系数α_s的计算公式，计算次同步振荡分量的频率f_s和次同步振荡分量的阻尼系数α_s的值；

步骤七、对同步相量X_r(m)进行Hann窗口截取并执行离散傅里叶变换操作，用次同步振荡相量的最大幅值谱定位第k_m谱线，得到并利用步骤六得到的f_s和α_s计算W_s(k_m)，反推得到次同步振荡分量幅值A_s的计算公式，计算次同步振荡分量幅值A_s的值。

2.如权利要求1所述的基于插值DFT的同步相量数据次同步振荡参数辨识法，其特征在于，所述步骤二获得的同步相量X_r(m)为：

X_r(m)＝[X_p(0)，…，X_p(mf_pr-1)，…，X_p((N_r-1)f_pr-1)]，m＝0，1，…，N_r-1。

3.如权利要求1所述的基于插值DFT的同步相量数据次同步振荡参数辨识法，其特征在于，所述步骤三得到离散傅里叶变换谱F_r(k)的公式为：

其中k为谱数，

其中F_r1是对应于基频分量的频谱，F_rs是对应次同步振荡分量的频谱。

4.如权利要求1所述的基于插值DFT的同步相量数据次同步振荡参数辨识法，其特征在于，所述步骤五获得的二阶差分比值R公式为：

利用泰勒级数展开的一阶项对F_rs进行重写，利用且k_m≈L_rs，

获得简化后的二阶差分比值R公式为：

其中

其中

5.如权利要求1所述的基于插值DFT的同步相量数据次同步振荡参数辨识法，其特征在于，所述步骤六中获得的次同步振荡分量的频率f_s和次同步振荡分量的阻尼系数α_s的计算公式分别为：

其中

Im(·)和Re(·)分别表示复变量的虚部和实部。