CN103942736B - 一种风电场多机等值建模方法 - Google Patents
一种风电场多机等值建模方法 Download PDFInfo
- Publication number
- CN103942736B CN103942736B CN201410194177.XA CN201410194177A CN103942736B CN 103942736 B CN103942736 B CN 103942736B CN 201410194177 A CN201410194177 A CN 201410194177A CN 103942736 B CN103942736 B CN 103942736B
- Authority
- CN
- China
- Prior art keywords
- group
- mrow
- planes
- msub
- cluster centre
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Landscapes
- Wind Motors (AREA)
Abstract
本发明属于电力系统仿真领域,尤其涉及一种风电场多机等值建模方法。该方法在机群划分过程中,从优化机群中心初始位置出发,综合考虑邻域半径和机群密集度因素进行初始机群中心的优化选择过程,并在后续的机群聚类中心更新过程中,基于目标函数在当前解处的伪梯度信息,对机群聚类中心进行优化搜索,提高了机群划分的精确性,优化了风电场多机等值建模方法,降低机群划分结果对初始机群聚类中心位置的敏感性,建立的风电场多机等值模型可以较完整地反映风电场的真实情况,在实测数据准确可信的前提下,建立的风电场多机等值模型可以更加准确地反映风电场的实际响应特性,模型的精确性较高。
Description
技术领域
本发明属于电力系统仿真领域,尤其涉及一种风电场多机等值建模方法。
背景技术
随着能源危机和环境污染的日益严重,风力发电受到世界各国的高度重视,随着风力发电和并网技术的发展,国内外学者对风电的研究也越来越关注。准确合理的风电机组模型是风电在电力系统中仿真研究的关键,也是其他深入研究的基础。
在风电场建模研究领域中,部分学者提出多机表征模型,其主要思想是以风电机组具有相近运行点为机群划分原则,采用某些聚类算法进行机群划分,将同群内的机组合并为一台等值机。目前,对机群划分最常用的方法是k-means聚类算法,该方法对处理大数据集具有快速、高效以及可伸缩性好等优点。但是k-means聚类算法对初始聚类中心非常敏感,如果机群聚类中心的初始位置选择不当,划分结果容易陷入局部。另外,采用传统的k-means聚类算法进行机群划分,后续的机群聚类中心更新方法一般是将同群内的机组样本组求均值,当同群内机组样本组分布不均匀时,将会导致机群聚类中心远离数据密集区,在很大程度上降低聚类质量。
发明内容
本发明的目的在于,针对目前风电场多机等值建模质量不高问题,提出了一种风电场多机等值建模方法,主要包括以下步骤:
步骤1、建立风电场多机等值模型的数学模型;
步骤2、初始机群聚类中心的优化选择过程;
步骤3、机群划分过程;
步骤4、更新机群聚类中心;
步骤5、算法迭代过程。
所述步骤1中,建立风电场多机等值模型的数学模型过程包括:
步骤101、假设风电场中风电机组共有n+m台,若某个时段内有n台风电机组并网运行,m台风电机组由于某些原因而与电网脱离,则将m台离网的风电机组数据剔除,选取风电场在该时段内并网运行的风电机组的实测风速数据作为样本,在该时段内实测风速数据的采样点数为t个,将n台风电机组的实测风速数据建立样本矩阵V:
其中vij表示第i台风电机组在第j个时刻的实测风速数据,把向量vi看作空间V的一个样本组,则V中样本组数量为n,维数为t;
步骤102、建立风电场多机等值模型的数学模型:
先确定各机群聚类中心,其当前的位置表示为一个解向量X,
X=(x1,x2,…,xa,…,xk) (2)
k(k<n)为机群划分的个数,xa为第a(a=1,2,…,k)个机群的聚类中心。
风电场机群划分目标函数定义为:
其中Ga代表第a(a=1,2,…,k)个机群,xa为第a个机群的聚类中心,vi为空间V中的第i个样本组,vi∈Ga表示vi是属于第a个机群内的样本组。
所述步骤2中,初始机群聚类中心的优化选择过程包括:
步骤201、定义两个样本组之间的欧氏距离db,c,
xbj表示第b个机群的聚类中心的第j个分量,xcj表示第c个机群的聚类中心的第j个分量,对于第i台风电机组样本组vi来说,其邻域半径为
其中di,max=max(di,e|e=1,2,…,n),表示V中与第i台风电机组样本组vi欧氏距离最远的样本组与vi的欧氏距离,μ为调节系数,取值范围在[0,1],初始值设为1,第i台风电机组样本组vi的邻域为以vi为中心,半径为ri的区域,则该机组的机群密集度为
Ni为在第i台风电机组样本组vi的邻域范围内的机组台数,ri为vi的邻域半径,选取ρi最大的样本组作为第一个初始机群中心;
步骤202、将第一个初始机群中心及其邻域内的所有样本组除去,并降低调节系数
重复执行步骤2中的以上步骤,找出第二个初始机群中心,依此类推,直到找到k个初始机群中心,
设根据以上步骤找到的k个机群的初始位置向量空间为
其中k为机群划分的个数,t为实测数据的采样点数,xa,g表示第a个机群聚类中心在第g时刻对应的分量。
所述步骤3中,机群划分过程包括:
根据欧氏距离进行机群划分,对于V中的任意一个样本组vi(i=1,2,…,n)而言,总是存在一个机群聚类中心xa(a=1,2,…,k)使d(vi,xa)达到最小,d(vi,xa)表示vi与xa之间的欧氏距离,则将vi划分到机群Ga中,Ga表示第a个机群,遍历完所有样本组,完成一次机群划分。
所述步骤4中,更新机群聚类中心包括:
更新机群的聚类中心,为不失一般性,以第a(a=1,2,…,k)个机群中心xa为例说明机群聚类中心优化更新过程,随机生成一个向量
Δxa=(Δxa,1,Δxa,2,…,Δxa,g,…,Δxa,t) (9)
其中
g(g=1,2,…,t),参数β(0<β<1)为机群聚类中心更新过程的步长,取值范围在[0,1],计算
Δxa,g为式(10)中的量,向量f′a,g(xa)=(f′a,1(xa,1),f′a,2(xa,2),…,f′a,g(xa,g),…,f′a,t(xa,t))为目标函数在当前解xa处的伪梯度,f′a,g(xa,g)为目标函数在当前解xa处的伪梯度的第g(g=1,2,…,t)个分量;令
x′a=(x′a,1,x′a,2,…x′a,g,…x′a,t) (12)
x′a,g=xa,g+β×sign(f′a,g(xa,g)) (13)
xa,g第a个机群聚类中心的g(g=1,2,…,t)个分量,x′a,g为x′a的第g(g=1,2,…,t)个分量,sign(y)为符号函数,其功能是取某个数的符号,
将xa和x′a分别代入式(3),计算目标函数值,当f(x′a)≥f(xa)时,重复执行该步骤,直到出现x′a使f(x′a)<f(xa),将xa更新为x′a,完成第a个机群聚类中心的更新;分别对k个机群聚类中心进行上述的更新过程,更新完之后转至步骤3进行机群划分。
所述步骤5中,算法迭代过程包括:
重复执行步骤4,直到相邻两次迭代的目标函数变化小于等于0.02或达到指定的循环次数。
本发明提出的风电场多机等值建模方法,其特点和效果是,本发明建立的风电场多机等值模型是基于风电机组实测运行数据的,能够反映风电机组的实际运行状态;本发明所采用的聚类方法可以降低机群划分结果对初始机群聚类中心位置的敏感性,并在后续的机群聚类中心更新过程中基于目标函数在当前解处的伪梯度信息,对机群聚类中心进行优化搜索,提高了机群划分的精确性,建立的风电场多机等值模型可以较完整地反映风电场的真实情况,在实测数据准确可信的前提下,建立的风电场多机等值模型可以更加准确地反映风电场的实际响应特性,模型的精确性较高。
附图说明
图1为本发明中风电场多机等值建模方法的基本步骤框图。
图2为使用本发明的实际风电场机组位置平面图。
图3为使用本发明的IEEE9风电场并网仿真系统。
图4为系统侧短路故障时,风电场采用传统的单机等值模型(single)、基于无监督谱聚类算法建立的多机表征模型(un-supervised)、本文基于半监督谱聚类算法建立的多机表征模型(semi-supervised)以及详细模型(detail)时有功出力的对比图。
具体实施方式
下面结合附图,具体说明本发明的实施方式。
本发明针对目前风电场建模聚类算法质量低的问题,提出了一种风电场多机等值建模方法。在机群划分过程中,该发明从优化机群中心初始位置出发,综合考虑邻域半径和机群密集度因素进行初始机群中心的优化选择过程。并在后续的机群聚类中心更新过程中,基于目标函数在当前解处的伪梯度信息,对机群聚类中心进行优化搜索,提高了机群划分的精确性,优化了风电场多机等值建模方法。
本发明提出了一种风电场多机等值建模方法的,主要包括以下步骤:
步骤1、建立风电场多机等值模型的数学模型;
步骤2、初始机群聚类中心的优化选择过程;
步骤3、机群划分过程;
步骤4、更新机群聚类中心;
步骤5、算法迭代过程。
步骤1中,建立风电场多机等值模型的数学模型过程为:
步骤101、假设风电场中风电机组共有n+m台,若某个时段内有n台风电机组并网运行,m台风电机组由于某些原因而与电网脱离,则将m台离网的风电机组数据剔除,选取风电场在该时段内并网运行的风电机组的实测风速数据作为样本,在该时段内实测风速数据的采样点数为t个,将n台风电机组的实测风速数据建立样本矩阵V:
其中vij表示第i台风电机组在第j个时刻的实测风速数据,把向量vi看作空间V的一个样本组,则V中样本组数量为n,维数为t。
步骤102、建立风电场多机等值模型的数学模型:
1、解的表述
在风电场机群划分过程中,需要先确定各机群聚类中心,其当前的位置表示为一个解向量X,
X=(x1,x2,…,xa,…,xk) (2)
k(k<n)为机群划分的个数,xa为第a(a=1,2,…,k)个机群的聚类中心。
2、目标函数
风电场机群划分的目的是将具有相近运行点的风电机组划分为同一机群,即同群内的机组与该机群的聚类中心的平方误差之和最小,故目标函数可定义为:
其中Ga代表第a(a=1,2,…,k)个机群,xa为第a个机群的聚类中心,vi为空间V中的第i个样本组,vi∈Ga表示vi是属于第a个机群内的样本组。
步骤2中,初始机群聚类中心的优化选择过程:
定义两个样本组之间的欧氏距离db,c,
xbj表示第b个机群的聚类中心的第j个分量,xcj表示第c个机群的聚类中心的第j个分量,对于第i台风电机组样本组vi来说,其邻域半径为
其中di,max=max(di,e|e=1,2,…,n),表示V中与第i台风电机组样本组vi欧氏距离最远的样本组与vi的欧氏距离,μ为调节系数,取值范围在[0,1],初始值设为1。第i台风电机组样本组vi的邻域为以vi为中心,半径为ri的区域,则该机组的机群密集度为
Ni为在第i台风电机组样本组vi的邻域范围内的机组台数,ri为vi的邻域半径。ρi越大,表明vi在邻域半径越小的情况下,邻域范围内聚集的机组样本组越多,即该样本组位于高密度区域,因此可选取ρi最大的样本组作为第一个初始机群中心。
将第一个初始机群中心及其邻域内的所有样本组除去,并降低调节系数
重复执行步骤2中的以上步骤,找出第二个初始机群中心,依此类推,直到找到k个初始机群中心,
设根据以上步骤找到的k个机群的初始位置向量空间为
其中k为机群划分的个数,t为实测数据的采样点数,xa,g表示第a个机群聚类中心在第g时刻对应的分量。
步骤3中,机群划分过程:
根据欧氏距离进行机群划分,对于V中的任意一个样本组vi(i=1,2,…,n)而言,总是存在一个机群聚类中心xa(a=1,2,…,k)使d(vi,xa)达到最小,d(vi,xa)表示vi与xa之间的欧氏距离,则将vi划分到机群Ga中,Ga表示第a个机群。遍历完所有样本组,完成一次机群划分。
步骤4中,更新机群聚类中心:
更新机群的聚类中心,为不失一般性,以第a(a=1,2,…,k)个机群中心xa为例说明机群聚类中心优化更新过程,随机生成一个向量
Δxa=(Δxa,1,Δxa,2,…,Δxa,g,…,Δxa,t) (9)
其中
g(g=1,2,…,t),参数β(0<β<1)为机群聚类中心更新过程的步长,取值范围在[0,1]。计算
Δxa,g为式(10)中的量,向量f′a,g(xa)=(f′a,1(xa,1),f′a,2(xa,2),…,f′a,g(xa,g),…,f′a,t(xa,t))为目标函数在当前解xa处的伪梯度,f′a,g(xa,g)为目标函数在当前解xa处的伪梯度的第g(g=1,2,…,t)个分量。令
x′a=(x′a,1,x′a,2,…x′a,g,…x′a,t) (12)
x′a,g=xa,g+β×sign(f′a,g(xa,g)) (13)
xa,g第a个机群聚类中心的g(g=1,2,…,t)个分量,x′a,g为x′a的第g(g=1,2,…,t)个分量,sign(y)为符号函数,其功能是取某个数的符号,
将xa和x′a分别代入式(3),计算目标函数值,当f(x′a)≥f(xa)时,重复执行该步骤的以上过程,直到出现x′a使f(x′a)<f(xa),将xa更新为x′a,完成第a个机群聚类中心的更新。分别对k个机群聚类中心进行上述的更新过程,更新完之后转至步骤3进行机群划分。
步骤5中,算法迭代过程:
重复执行步骤4,直到相邻两次迭代的目标函数变化不明显(前一次与后一次的目标函数差值为某一很小的正实数,一般取0.02)或达到指定的循环次数(人为预先指定该循环次数)。
下面通过一个实际风电场来说明本发明提出的风电场动态等值建模方法。
以某地区实际风电场为例进行分析,该风电场内共有33台风电机组,其中机组类型为GE1.5MW,风电场总装机容量为49.5MW,机组位置平面图如图2所示。
选取2011年9月1日至2011年9月30日的实测风速数据进行分析,对场内的风电机组采用本发明提出的风电场动态等值建模方法进行机群划分,根据机群划分结果,可将风电场等值为三台风电机组,采用容量加权法计算等值风电机组的参数。将上述等值的3台风电机组并联接在图3所示的IEEE9风电场并网仿真系统的母线3上。
图4为系统侧母线9发生三相接地短路故障下,风电场采用传统的k-means聚类算法建立的多机等值模型、本发明提出的方法建立的多机等值模型以及详细模型时,风电场有功功率动态响应曲线。
从图4可以看出,在系统侧发生三相接地短路故障下,与传统的k-means聚类算法建立的多机等值模型相比,风电场采用本发明提出的方法建立的多机等值模型时,其动态响应特性与详细模型更接近,模型的精确性有了一定的提高。
以上所述仅为本发明的较佳实施例而已,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内所作的任何修改、等同替换、改进等,均包含在本发明的保护范围内。
Claims (5)
1.一种风电场多机等值建模方法,其特征在于,所述方法主要包括:
步骤1、建立风电场多机等值模型的数学模型;
步骤2、初始机群聚类中心的优化选择过程;
步骤3、机群划分过程;
步骤4、更新机群聚类中心;
步骤5、算法迭代过程;
所述步骤2中,初始机群聚类中心的优化选择过程包括:
步骤201、定义两个样本组之间的欧氏距离db,c,
<mrow>
<msub>
<mi>d</mi>
<mrow>
<mi>b</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>d</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>b</mi>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>c</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msqrt>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>j</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>t</mi>
</munderover>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>b</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mrow>
<mi>c</mi>
<mi>j</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>,</mo>
<mrow>
<mo>(</mo>
<mi>b</mi>
<mo>,</mo>
<mi>c</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>n</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
n为风电机组台数,V为将n台风电机组的实测风速数据建立样本矩阵,xbj表示第b个机群的聚类中心的第j个分量,xcj表示第c个机群的聚类中心的第j个分量,对于第i台风电机组样本组vi来说,其邻域半径为
<mrow>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>d</mi>
<mrow>
<mi>i</mi>
<mo>,</mo>
<mi>m</mi>
<mi>a</mi>
<mi>x</mi>
</mrow>
</msub>
<mrow>
<mn>2</mn>
<mo>&times;</mo>
<mi>k</mi>
<mo>&times;</mo>
<mi>&mu;</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
其中di,max=max(di,e|e=1,2,…,n),表示V中与第i台风电机组样本组vi欧氏距离最远的样本组与vi的欧氏距离,μ为调节系数,取值范围在[0,1],初始值设为1,第i台风电机组样本组vi的邻域为以vi为中心,半径为ri的区域,则该机组的机群密集度为
<mrow>
<msub>
<mi>&rho;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>N</mi>
<mi>i</mi>
</msub>
<msub>
<mi>r</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
Ni为在第i台风电机组样本组vi的邻域范围内的机组台数,ri为vi的邻域半径,选取ρi最大的样本组作为第一个初始机群中心;
步骤202、将第一个初始机群中心及其邻域内的所有样本组除去,并降低调节系数
<mrow>
<msup>
<mi>&mu;</mi>
<mo>&prime;</mo>
</msup>
<mo>=</mo>
<mi>&mu;</mi>
<mo>-</mo>
<mfrac>
<mn>1</mn>
<mi>k</mi>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
重复执行步骤2中的以上步骤,找出第二个初始机群中心,依此类推,直到找到k个初始机群中心,
设根据以上步骤找到的k个机群的初始位置向量空间X为
其中k为机群划分的个数,t为实测数据的采样点数,xa,g表示第a个机群聚类中心在第g时刻对应的分量。
2.根据权利要求1所述方法,其特征在于,所述步骤1中,建立风电场多机等值模型的数学模型过程包括:
步骤101、假设风电场中风电机组共有n+m台,若某个时段内有n台风电机组并网运行,m台风电机组由于某些原因而与电网脱离,则将m台离网的风电机组数据剔除,选取风电场在该时段内并网运行的风电机组的实测风速数据作为样本,在该时段内实测风速数据的采样点数为t个,将n台风电机组的实测风速数据建立样本矩阵V:
其中vij表示第i台风电机组在第j个时刻的实测风速数据,把向量vi看作空间V的一个样本组,则V中样本组数量为n,维数为t;
步骤102、建立风电场多机等值模型的数学模型:
先确定各机群聚类中心,其当前的位置表示为一个解向量X,
X=(x1,x2,…,xa,…,xk) (2)
k(k<n)为机群划分的个数,xa为第a(a=1,2,…,k)个机群的聚类中心;
风电场机群划分目标函数定义为:
<mrow>
<mi>min</mi>
<mi> </mi>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mn>1</mn>
</msub>
<mo>,</mo>
<msub>
<mi>x</mi>
<mn>2</mn>
</msub>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>a</mi>
</msub>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<msub>
<mi>x</mi>
<mi>k</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>a</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>k</mi>
</munderover>
<munder>
<mo>&Sigma;</mo>
<mrow>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mo>&Element;</mo>
<msub>
<mi>G</mi>
<mi>a</mi>
</msub>
</mrow>
</munder>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>v</mi>
<mi>i</mi>
</msub>
<mo>-</mo>
<msub>
<mi>x</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
其中Ga代表第a(a=1,2,…,k)个机群,xa为第a个机群的聚类中心,vi为空间V中的第i个样本组,vi∈Ga表示vi是属于第a个机群内的样本组。
3.根据权利要求1所述方法,其特征在于,所述步骤3中,机群划分过程包括:
根据欧氏距离进行机群划分,对于V中的任意一个样本组vi(i=1,2,…,n)而言,总是存在一个机群聚类中心xa(a=1,2,…,k)使d(vi,xa)达到最小,d(vi,xa)表示vi与xa之间的欧氏距离,则将vi划分到机群Ga中,Ga表示第a个机群,遍历完所有样本组,完成一次机群划分,n为风电机组台数,V为将n台风电机组的实测风速数据建立样本矩阵,k为机群划分的个数。
4.根据权利要求1所述方法,其特征在于,所述步骤4中,更新机群聚类中心包括:
更新机群的聚类中心,为不失一般性,以第a(a=1,2,…,k)个机群中心xa为例说明机群聚类中心优化更新过程,随机生成一个向量
Δxa=(Δxa,1,Δxa,2,…,Δxa,g,…,Δxa,t) (9)
其中
g(g=1,2,…,t),参数β(0<β<1)为机群聚类中心更新过程的步长,取值范围在[0,1],计算
<mrow>
<msubsup>
<mi>f</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>g</mi>
</mrow>
<mo>&prime;</mo>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfrac>
<mrow>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>a</mi>
</msub>
<mo>+</mo>
<msub>
<mi>&Delta;x</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mi>f</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>x</mi>
<mi>a</mi>
</msub>
<mo>-</mo>
<msub>
<mi>&Delta;x</mi>
<mi>a</mi>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>2</mn>
<msub>
<mi>&Delta;x</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>g</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
Δxa,g为式(10)中的量,向量f′a,g(xa)=(f′a,1(xa,1),f′a,2(xa,2),…,f′a,g(xa,g),…,f′a,t(xa,t))为目标函数在当前解xa处的伪梯度,f′a,g(xa,g)为目标函数在当前解xa处的伪梯度的第g(g=1,2,…,t)个分量;令
xa′=(x′a,1,x′a,2,…x′a,g,…x′a,t) (12)
x′a,g=xa,g+β×sign(f′a,g(xa,g)) (13)
xa,g第a个机群聚类中心的g(g=1,2,…,t)个分量,x′a,g为xa′的第g(g=1,2,…,t)个分量,sign(y)为符号函数,其功能是取某个数的符号,
<mrow>
<mi>s</mi>
<mi>i</mi>
<mi>g</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
<mtd>
<mrow>
<mi>y</mi>
<mo>></mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>y</mi>
<mo><</mo>
<mn>0</mn>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
将xa和xa′分别代入式(3),计算目标函数值,当f(xa′)≥f(xa)时,重复执行该步骤,直到出现xa′使f(xa′)<f(xa),将xa更新为xa′,完成第a个机群聚类中心的更新;
分别对k个机群聚类中心进行上述的更新过程,更新完之后转至步骤3进行机群划分,k为机群划分的个数。
5.根据权利要求1所述方法,其特征在于,所述步骤5中,算法迭代过程为:重复执行步骤4,直到相邻两次迭代的目标函数变化小于等于0.02或达到指定的循环次数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410194177.XA CN103942736B (zh) | 2014-05-08 | 2014-05-08 | 一种风电场多机等值建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410194177.XA CN103942736B (zh) | 2014-05-08 | 2014-05-08 | 一种风电场多机等值建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103942736A CN103942736A (zh) | 2014-07-23 |
CN103942736B true CN103942736B (zh) | 2017-10-10 |
Family
ID=51190388
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410194177.XA Expired - Fee Related CN103942736B (zh) | 2014-05-08 | 2014-05-08 | 一种风电场多机等值建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103942736B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105184683A (zh) * | 2015-10-10 | 2015-12-23 | 华北电力科学研究院有限责任公司 | 一种基于风电场运行数据的概率聚类分群方法 |
CN105468867A (zh) * | 2015-12-21 | 2016-04-06 | 华北电力大学(北京) | 一种基于clarans聚类的风电场机组聚合方法 |
CN105703396B (zh) * | 2016-03-21 | 2018-01-16 | 广西电网有限责任公司电力科学研究院 | 一种双馈风力发电场多机表征等值建模方法 |
CN108664671B (zh) * | 2017-03-28 | 2021-01-12 | 中国电力科学研究院有限公司 | 一种风电场多机聚合模型参数辨识方法及装置 |
CN107272643B (zh) * | 2017-06-19 | 2019-03-29 | 哈尔滨工业大学 | 新能源场站的单机等值方法 |
CN108460228B (zh) * | 2018-03-21 | 2021-06-04 | 电子科技大学 | 一种基于多目标优化算法进行风电场等值的方法 |
CN109409575A (zh) * | 2018-09-27 | 2019-03-01 | 贵州电网有限责任公司 | 基于Gap Statistic的风电场机群划分方法 |
CN111950131B (zh) * | 2020-07-20 | 2022-06-10 | 东南大学 | 一种考虑限电因素的风电场出力等值聚合模型构建方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103400009A (zh) * | 2013-08-07 | 2013-11-20 | 华北电力大学 | 基于分裂层次半监督谱聚类算法的风电场动态等值方法 |
-
2014
- 2014-05-08 CN CN201410194177.XA patent/CN103942736B/zh not_active Expired - Fee Related
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103400009A (zh) * | 2013-08-07 | 2013-11-20 | 华北电力大学 | 基于分裂层次半监督谱聚类算法的风电场动态等值方法 |
Non-Patent Citations (1)
Title |
---|
基于运行数据的风电场等效建模研究;陈迎;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20130315(第3期);8,13-16 * |
Also Published As
Publication number | Publication date |
---|---|
CN103942736A (zh) | 2014-07-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103942736B (zh) | 一种风电场多机等值建模方法 | |
CN104978608B (zh) | 一种风电功率预测装置及预测方法 | |
Shi et al. | Hybrid forecasting model for very-short term wind power forecasting based on grey relational analysis and wind speed distribution features | |
CN103268366A (zh) | 一种适用于分散式风电场的组合风电功率预测方法 | |
CN103400009B (zh) | 基于分裂层次半监督谱聚类算法的风电场动态等值方法 | |
CN106228278A (zh) | 光伏功率预测系统 | |
WO2024051524A1 (zh) | 一种水风光发电功率联合预测方法、装置 | |
CN110264002B (zh) | 基于聚类分析的风电场微观选址方案评价方法 | |
CN102012956A (zh) | 一种考虑风电场输入风速和风向随机波动的风电场等值方法 | |
CN105069521A (zh) | 一种基于加权fcm聚类算法的光伏电站输出功率预测方法 | |
CN103489046A (zh) | 风电场短期功率预测方法 | |
CN109272139A (zh) | 一种基于非线性集成深度学习的短期风速预测方法 | |
CN115425680B (zh) | 一种多能源联合发电系统的功率预测模型构建及预测方法 | |
US20240054267A1 (en) | Method for planning a layout of a renewable energy site | |
Hu et al. | Short-term photovoltaic power prediction based on similar days and improved SOA-DBN model | |
CN106849066A (zh) | 一种区域风电功率预测方法 | |
CN103440405A (zh) | 基于实测数据的风电场稳态建模方法及系统 | |
CN103955521A (zh) | 一种风电场机群划分方法 | |
Chen et al. | Research on wind power prediction method based on convolutional neural network and genetic algorithm | |
CN110909310A (zh) | 一种基于模型参数优化的光伏短期发电量预测方法及系统 | |
Sarkar et al. | Proficiency assessment of adaptive neuro-fuzzy inference system to predict wind power: A case study of Malaysia | |
Liu et al. | A weight-varying ensemble method for short-term forecasting PV power output | |
Das et al. | Optimized support vector regression-based model for solar power generation forecasting on the basis of online weather reports | |
CN103530508B (zh) | 一种建立风速功率转换概率模型的方法 | |
Sun et al. | Short-term photovoltaic power prediction modeling based on AdaBoost algorithm and Elman |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171010 Termination date: 20180508 |