计及山地风资源特性的风电场无功控制方法
【专利摘要】本发明提供一种计及山地风资源特性的风电场无功补偿的控制方法，包括以下步骤：对目标山地风电场各风机的实测历史风速序列进行特征分解，获得风速序列的空间特征向量的各分量及各分量对应的方差贡献率；按照各空间特征向量方差贡献率由大到小排列，选出第一个方差贡献率最高的空间特征向量；根据所属方差贡献率最高的空间特征向量的各分量相对大小安排风机的出力，并对各个风电场出力水平进行离线优化计算，得到各风机的无功出力值；根据各风机的无功出力值对风机及无功补偿设备进行控制。
【专利说明】计及山地风资源特性的风电场无功控制方法
【技术领域】
[0001]本发明涉及一种风电场无功的控制方法，尤其涉及一种计及山地风资源特性的风电场无功的控制方法，属于电力系统新能源发电领域。
【背景技术】
[0002]风电是一种间歇式能源，出力具有随机性和波动性的特征，在传统的普通风电场中由于集电线路很短，这种波动性所带来的风电场内的电压问题并不明显。
[0003]然而随着山地风电场的发展，山地风电场具有集电线路长、风资源分散性强、变化快的特点，这导致了山地风电场内有较明显的无功电压问题，已经发生过多起由于电压不合格导致的部分风机脱网事故，因此充分考虑山地风电场内地形和风资源的差异，建立合适的无功电压控制方法和策略，对山地风电的开发、指导实际山地风电场运行具有重大意义。
【发明内容】
[0004]综上所述，确有必要提供一种能够充分考虑山地风电场内地形和风资源的差异的风电场无功控制方法。
[0005]一种计及山地风资源特性的风电场无功补偿的控制方法，包括以下步骤:对目标山地风电场各风机的实测历史风速序列进行特征分解，获得风速序列的空间特征向量的各分量及各分量对应的方差贡献率；按照各空间特征向量方差贡献率由大到小排列，选出第一个方差贡献率最高的空间特征向量；根据所属方差贡献率最高的空间特征向量的各分量相对大小安排风机的出力，并对各个风电场出力水平进行离线优化计算，得到各风机的无功出力值；根据各风机的无功出力值对风机及无功补偿设备进行控制。
[0006]相对于现有技术，本发明提供的风电场无功控制方法，以山地风资源特性为基础，利用特征分解在挖掘资源共性与特性方面的优势，进行山地风资源特性提取，适应存在明显地形起伏导致风资源分散性强、风速变化快的山地风电场，改善了山地风电场内的电压分布，提高了优化控制的速度。
【专利附图】
【附图说明】
[0007]图1为本发明提供的计及山地风资源特性的风电场无功控制方法的流程图。
【具体实施方式】
[0008]下面根据说明书附图并结合具体实施例对本发明的技术方案进一步详细表述。
[0009]请参阅图1，本发明提供的计及山地风资源特性的风电场无功控制方法，包括如下步骤:
[0010]步骤S10，对目标山地风电场各风机的实测历史风速序列进行特征分解，获得风速序列的空间特征向量及与各空间特征向量对应的方差贡献率；
[0011]步骤S20，按照各空间特征向量方差贡献率由大到小排列，选出第一个方差贡献率最闻的空间特征向量；
[0012]步骤S30，按照空间特征向量的相对大小安排风机的出力，并对各个风电场出力水平进行离线优化计算，得到各风机的无功出力值；
[0013]步骤S40，根据各风机的无功出力值对风机及无功补偿设备进行控制。
[0014]在步骤SlO中，所述空间特征向量及其对应的方差贡献率的获取具体包括如下步骤:
[0015]步骤S11，对风电场的风机进行编号，并得到风机的实测风速序列。
[0016]对风电场的所有风机进行编号，假设风电场有m台风机，则某时刻t所述风电场内所有风机的风速可以用一个m维向量Xt来表示，设已知有η个时刻的实测风速序列:
[0017](Xi t，Xgt，...，Χ?」.，...，Xmt)，t_ I j 2 j...j π (I)
[0018]式中，Xij为i号风机在时刻j时的风速。
[0019]步骤S12，根据风机的实测风速序列，构建实测风速序列的协方差矩阵Σ。
[0020]所述协方差矩阵Σ为一个m行m列的矩阵，它的第i行第j列元素具体表示如下:
[0021 ]Xij = XilXj^Xi2Xj2+...+XinXjn (2)
[0022]步骤S13，计算协方差矩阵Σ的特征值λ。
[0023]所述特征值λ可通过解以下方程式得到:
[0024]Σ-λ I = O (3)
[0025]式中，I为单位矩阵。设解式⑶可得到k个解，记为λ i，i = 1，2，…k。
[0026]步骤S14，针对每个特征值λ i，求特征值λ j对应的单位特征向量Vp
[0027]所述单位特征向量Vi可通过以下公式获取:
[0028](Σ-λ I)V = O (4)
[0029]k个特征值λ可以得到k个单位特征向量，记为Vi, i = 1，2，…k，每个单位特征向量Vi又称为描述空间特征的空间特征向量。
[0030]步骤S15，计算每个单位特征向量Vi的方差贡献率Qi。
[0031]所述单位特征向量Vi的方差贡献率Qi可通过以下公式计算:
【权利要求】
1.一种计及山地风资源特性的风电场无功控制方法，包括以下步骤: 对目标山地风电场各风机的实测历史风速序列进行特征分解，获得风速序列的空间特征向量及与个空间特征向量对应的方差贡献率； 按照各空间特征向量的方差贡献率由大到小排列，选出第一个方差贡献率最高的空间特征向量； 按照空间特征向量的相对大小安排风机的出力，并对各个风电场出力水平进行离线优化计算，得到各风机的无功出力值； 根据各风机的无功出力值对风机及无功补偿设备进行控制。
2.如权利要求1所述的计及山地风资源特性的风电场无功控制方法，其特征在于，风速序列的空间特征向量及对应的方差贡献率包括如下步骤: 对风电场的风机进行编号，并得到风机的实测风速序列； 根据风机的实测风速序列，构建实测风速序列的协方差矩阵Σ ； 计算协方差矩阵Σ的特征值Ai; 针对每个特征值Xi，求特征值对应的多个单位特征向量Vi ;以及 计算每个单位特征向量Vi的方差贡献率Q”
3.如权利要求2所述的计及山地风资源特性的风电场无功控制方法，其特征在于，设风电场有m台风机，以m维向量Xt来表示某时刻t所述风电场内所有风机的风速，则η个时刻的实测风速序列:
Xt (Xit，^2t?...，Xij，...，Xmt)，t_I，2，...，η ; 式中，Xu为i号风机在时刻j时的风速。
4.如权利要求2所述的计及山地风资源特性的风电场无功控制方法，其特征在于，所述协方差矩阵Σ为一个m行m列的矩阵，它的第i行第j列元素具体表示如下:
Σ Jj 一 XiiXji+Xi2Xj2+...+XinXjn。
5.如权利要求4所述的计及山地风资源特性的风电场无功控制方法，其特征在于，所述特征值λ通过解以下方程式得到:
Σ-λ 11 = O, 式中，I为单位矩阵，解方程式得到k个解，记为Xi，i = l，2r-k。
6.如权利要求5所述的计及山地风资源特性的风电场无功控制方法，其特征在于，所述单位特征向量V可通过以下公式获取:
(Σ-λ I) V = O, 其中，k个特征值λ可以得到k个单位特征向量,记为Vi, i = 1，2，…k,每个单位特征向量Vi即为描述空间特征的空间特征向量。
7.如权利要求6所述的计及山地风资源特性的风电场无功控制方法，其特征在于，所述单位特征向量Vi的方差贡献率Qi通过以下公式计算:
8.如权利要求7所述的计及山地风资源特性的风电场无功控制方法，其特征在于，按照方差贡献率Qi的大小将单位特征向量Vi由大到小进行排列，并选出方差贡献率Qi最高的单位特征向量Vmax，将其余单位特征向量舍去不考虑。
9.如权利要求1所述的计及山地风资源特性的风电场无功控制方法，其特征在于，各个风电场出力水平进行离线优化计算包括如下步骤: 将风电场可能的出力水平范围(O~Pmax)进行离散化处理: ? = 0，1，2，"％?_(单位:丽)， 其中P为风电场的整体出力水平，Pmax为风电场的最大出力； 将所选出的空间特征向量表示为:
Vnlax = (Vl，V2.…，Vm)T， 根据每个风电场的整体出力水平P计算对应的各风机的假想出力Pi为:
Pi = P*Vi/ Ol+V2 +…+Vm)， 式中Pi为第i台风机的假想出力； 在各风机假想有功出力情况下，分别对风机的无功出力进行离线优化，优化模型为:
式中g为风机的台数，Ui为节点i电压实测值，Uref为机端电压参考值，Pei和Qei分别是节点i风机发出的有功功率和无功功率，Pu和Qu分别是是节点i的有功负荷和无功负荷，Qci为节点i的无功补偿装置发出的无功补偿值，Gij和Bu分别为节点i和节点j间的节点电导和节点电纳，S u是节点i和节点j之间电压的相角差，Qcd min和Qcd max分别为无功补偿装置补偿量的上下限，Qei min和Qei _分别为双馈发电机无功出力的上下限。
10.如权利要求9所述的计及山地风资源特性的风电场无功控制方法，其特征在于，风机及无功补偿设备进行控制包括如下步骤: 测得风电场实时有功出力，找到与实时有功出力对应的离散出力水平P ; 调取离散出力水平P下的离散优化结果，以该离散优化结果作为初值，按风机的无功出力的离线优化模型进行计算，得到风机及无功补偿设备的无功出力值； 根据风机及无功补偿设备的无功出力值对风机及无功补偿设备进行控制。
【文档编号】H02J3/16GK104201686SQ201410356118
【公开日】2014年12月10日 申请日期:2014年7月24日 优先权日:2014年7月24日
【发明者】乔颖, 鲁宗相, 吴林林 申请人:清华大学, 国家电网公司, 华北电力科学研究院有限责任公司