主权项 |
一种风电场机群划分方法,其特征在于,所述方法包括以下步骤:步骤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>Φ</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow><mo>=</mo><munderover><mo>Σ</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>Σ</mo><mrow><mi>a</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><munderover><mo>Σ</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>η</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>Φ</mi><mrow><mo>(</mo><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>)</mo></mrow></mrow><mrow><munderover><mo>Σ</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><mi>Φ</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>)<η<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>为一台风电机组所有时刻的实测数据,即一个机组样本组。 |