发明名称 一种风电场机群划分方法
摘要 本发明属于电力系统仿真领域,尤其涉及一种风电场机群划分方法,该方法基于风电场实测运行数据对场内机组进行机群划分。在数据的采集过程中,由于测量系统的缺陷或执行错误等因素,实测数据可能会含有噪声数据,为减小噪声数据的干扰,本发明首先基于样本点处的势值对风电场实测运行数据进行孤立点数据处理。在机群划分时,当两个机群中心初始位置较近时,其包含的冗余信息较多,划分结果容易陷入局部最优,本发明针对此将欧氏距离最小的样本组向均值点移动,并用移动后的样本组的均值替换原来的样本组,通过该方法获得具有多样性的初始机群中心的位置,提高了全局搜索能力。采用本发明提出的风电场机群划分方法能够将具有相近运行点的风电机组划分为同一机群,优化了风电场等值建模方法。
申请公布号 CN103955521B 申请公布日期 2017.04.12
申请号 CN201410193280.2 申请日期 2014.05.08
申请人 华北电力大学 发明人 林俐;潘险险;赵双;张凌云;李丹;李亮玉;李凯;吴聪聪;邹兰青;周鹏;李诗童
分类号 G06F17/30(2006.01)I;G06F17/50(2006.01)I 主分类号 G06F17/30(2006.01)I
代理机构 北京众合诚成知识产权代理有限公司 11246 代理人 朱琨
主权项 一种风电场机群划分方法,其特征在于,所述方法包括以下步骤:步骤1、风电场实测数据孤立点数据处理;步骤2、初始机群中心优化过程;步骤3、机群划分,判断是否满足终止条件;如果是,则机群划分结束;否则执行步骤4;步骤4、更新机群中心;所述步骤1中,风电场实测数据孤立点数据处理过程为:步骤101、假设某时段内风电场中有n台风电机组并网运行,选取该时段内并网运行风电机组的实测运行数据,包括风速、有功功率或无功功率作为样本,建立实测运行数据样本矩阵:<img file="FDA0001188202590000011.GIF" wi="1582" he="263" />x<sub>ij</sub>表示第i台风电机组在第j个时刻的实测运行数据,n为场内风电机组台数,t为该时段内实测运行数据的采集时间点的个数;步骤102、对于某一时刻的所有机组实测数据X<sub>j</sub>而言,可定义样本点x<sub>ij</sub>(i=1,2,…,n)处的势函数为<maths num="0001"><math><![CDATA[<mrow><mi>&Phi;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>a</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mi>e</mi><mrow><mo>-</mo><mfrac><mn>4</mn><msup><mi>R</mi><mn>2</mn></msup></mfrac><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>a</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001188202590000012.GIF" wi="1203" he="159" /></maths>其中<maths num="0002"><math><![CDATA[<mrow><mi>R</mi><mo>=</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msqrt><mrow><mfrac><mn>1</mn><mrow><mi>n</mi><mrow><mo>(</mo><mi>n</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><munderover><mo>&Sigma;</mo><mrow><mi>a</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>-</mo><msub><mi>x</mi><mrow><mi>a</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow></msqrt><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001188202590000013.GIF" wi="1289" he="150" /></maths>Φ(x<sub>ij</sub>)为样本点x<sub>ij</sub>处的势函数,R表示邻域半径,x<sub>ij</sub>表示第i台风电机组在第j个时刻的实测运行数据,x<sub>aj</sub>表示第a台风电机组在第j个时刻的实测运行数据,当样本点x<sub>ij</sub>处的势越大,表明x<sub>ij</sub>附近聚集的样本点越多;步骤103、计算x<sub>ij</sub>处的势占该时刻所有数据样本点处势总和的比值,定义为势比,<maths num="0003"><math><![CDATA[<mrow><mi>&eta;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>&Phi;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><munderover><mo>&Sigma;</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>&Phi;</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0001188202590000021.GIF" wi="1086" he="191" /></maths>根据此方法计算每个数据样本点的势比,并将其降序排列η<sub>1</sub>≥η<sub>2</sub>≥…≥η<sub>n×t</sub>≥0,计算该排列中相邻势比的差值,定义为势比差Δη<sub>g</sub>,Δη<sub>g</sub>=η<sub>g</sub>‑η<sub>g+1</sub>(g=1,2,…,n×t‑1)   (5)若Δη<sub>b</sub>(b∈g)最大,则可确定势比差阈值为η<sub>b</sub>,即当η(x<sub>ij</sub>)&lt;η<sub>b</sub>(i=1,2,…,n;j=1,2,…,t)时,可确定x<sub>ij</sub>为孤立点数据,为不影响实测数据在时间上的连续性,对于孤立点数据,设其为某台机组在某一时刻的实测数据,采用该台机组在其前后时刻的实测数据平均值替换该孤立点数据,若存在某几个连续时刻测得的数据均为孤立点数据,则选择其最邻近的前后时刻的正常数据进行等分取值;为表达方便,经过孤立点数据处理后的风电场实测运行数据仍用式(1)表示,样本组换成机组样本组<img file="FDA0001188202590000022.GIF" wi="1557" he="262" />x<sub>ij</sub>表示经过孤立点数据处理后第i台风电机组在第j个时刻的实测运行数据,n为场内风电机组台数,t为该时段内实测运行数据的采集时间点的个数,X<sub>i</sub>为一台风电机组所有时刻的实测数据,即一个机组样本组。
地址 102206 北京市昌平区朱辛庄北农路2号