[发明专利]一种基于db6小波的风电汇聚效应时频分析方法

摘要

本发明涉及一种基于db6离散小波变换的风电集群功率汇聚效应的时频分析方法，针对当前缺乏从时频域角度分析风电集群功率汇聚效应的问题，本发明以db6小波的特性及其离散变换原理为基础，综合考虑系统的典型调峰、调频时间尺度，确立了面向电网运行的风电频段划分，建立时频分析指标，以此在各个频段内对风电功率的汇聚效应进行时频特性分析。本发明方法通过分析风电汇聚效应时频特性，对于大规模风电集群并网的规划、安全运行及电网调度具有重要的指导意义。

主权项

1.一种基于db6小波的风电汇聚效应时频分析方法，其特征是：分析db6小波的特性及其离散变换原理，分析结合系统的典型调峰、调频时间尺度，建立面向电网运行的风电频段划分，在此基础上建立时频特性分析指标，并在各个频段内对风电汇聚效应的时频特性进行分析，具体包括以下步骤：1)db6小波及其离散变换原理分析风电功率是典型的不规则信号，为了把握各频段内风电功率的时频特性，利用dbN小波对风电功率进行多层分解，在不同的频率范围内对风电功率进行时频分析；db6小波为紧支撑正交实小波，正则性好、消失矩大，能够用于风电场出力的时频域特性分析；将db6小波作为基本小波，通过改变伸缩因子、平移因子的大小，使得风电功率信号通过不同的带通滤波器，从而将不同频段的风电功率信号进行分离，进而在各个频段内对风电汇聚效应进行时频分析；2)面向电网运行的风电频段划分对风电功率依靠db6小波进行二进制离散变换，随着小波分解层数的增多，各层细节信号的频率以1/2的幂指数递减，结合系统典型的调峰、调频尺度，对时间尺度为1s的信号进行db6小波11层分解，将分解得到的前5层细节信号合并作为风电的高频信号，影响系统的一次调频，第6～9层细节信号合并作为次高频信号，影响系统的二次调频，第10～11层细节信号合并作为中频信号，影响系统的短时间尺度调峰，第11层逼近信号作为低频信号，影响系统的中长期调峰；3)汇聚效应时频特性分析指标体系随着集群汇聚规模的增大，不同位置风电机组的出力波动存在平抑效果，集群风电的输出功率会逐渐趋于平缓，称为“汇聚效应”，为了分析风电汇聚效应的时频特性，需要在不同频段内对不同集群规模的风电功率进行分析，因此建立如下时频指标进行分析；①各频段功率能量占比E_i^*小波变换能够将信号无泄露、不重叠的正交分解到各个频段内，并且满足能量守恒定律，根据Parsecal能量积分公式能够得到，信号f(t)在时域内的能量E表示为：其中：R为信号f(t)积分区间，t为时间序列；离散信号在小波域用小波分解得到的各个频段的信号序列，某频段序列内各个节点系数的平方和等于对应频段在时间域内的能量：(2)式中E_i表示频段i在时间域内的能量，N表示总节点数，C_i,j表示频段i内第j个节点的系数；为了方便将各个频段内的数据进行对比，对各个频段内的能量进行归一化处理：能量占比越大，此频段所占的能量比重越高；②各频段功率极值差ΔP_i,max经过离散小波变换后的各频段功率存在波动，将各频段内功率的最大值与最小值之差定义为该频段功率的极值差ΔP_max；极值差越大，说明在此频段的风电功率波动范围越广，波动越明显；四个频段的功率序列{P_i}能够进行如下表示：{P_i}＝{p_i1,p_i2,...,p_in}(i＝1,2,3,4)             (4)其中：p_in为频段i中第n个数据点的功率值，n为功率信号序列的数据个数；第i个频段功率的极值差ΔP_i,max表示为：ΔP_i,max＝p_i,max‑p_i,min(i＝1,2,3,4)           (5)其中：p_i,max为第i个频段功率的最大值，p_i,min为第i个频段功率的最小值；③各频段功率波动率V_i为了反映各个频段功率的波动情况，定义各频段功率波动率指标V_i，步骤为：将风电频段i内第j+1个数据点与第j个数据点之间的功率变化标幺值记为X_ij^*：X_ij^*＝(P_i,j+1‑P_i,j)/P_N(j＝1,2,···n‑1)               (6)其中：P_i,j+1表示第i个频段第j+1个数据点的功率，P_i,j表示第i个频段第j个数据点的功率，n为数据点个数，P_N为装机容量；功率变化标幺值序列的均值记为功率变化序列的标准差记为σ_i：其中：X_in‑1^*为风电频段i内第n个数据点与第n‑1个数据点之间的功率变化标幺值；风电在第i个频段的波动率记为V_i：④各频段功率平滑系数S_i为了衡量各频段风电出力的平滑效果，定义频段i的平滑系数指标S_i：其中S_i(i＝1,2,3,4)分别表示低、中、次高、高频段相应的平滑系数，平滑系数越大，说明集群风电相对于单机波动的改善效果越明显；R_i,单机为单机在频段i内的风电波动置信区间，R_i,机群为机群在频段i内的风电波动置信区间，波动置信区间R的计算如下：P(|X|≤R)＝p                    (11)X＝(P(t+1)‑P(t))/P_N                  (12)其中：P表示概率，X表示风电差额功率标幺值，P(|X|≤R)表示差额功率标幺值的绝对值不超过R的概率，p为给定概率值，P(t+1)表示第t+1个数据点的功率，P(t)表示第t个数据点的功率，P_N为额定装机容量；⑤汇聚效应时频特征因子F极值差、波动率、平滑系数同为描述波动性的指标，将三者进行加权得到综合波动指标q_i：q_i＝λ_iΔP_i,max+β_iV_i‑γ_iS_i               (13)其中：λ_i,β_i,γ_i分别为极值差ΔP_i,max、波动率V_i、平滑系数S_i对应的权重系数，采用熵权法进行客观定权，q_i越大，表示频段i的波动性对系统造成的影响越大；数值越小，表明对系统影响越小；理论上该指标允许出现负值，即在加权后的极值差与波动率之和小于加权后的平滑系数时，综合波动指标q_i为负；在相同综合波动指标下，某频段能量占比E_i^*越大，则对系统的影响也越大，因此将能量占比与综合波动指标的乘积，作为特征因子f_i以评判频段i风电对系统的影响程度，f_i定义为：各个频段特征因子之和定义为汇聚效应的时频特征因子F：F值越大，该风电集群功率各频段总体波动特性对系统影响越大。