一种风电场群与火力发电厂联合发电系统的平滑出力方法

摘要

本发明提供一种风电场群与火力发电厂联合发电系统的平滑出力方法，包括：针对风电场群中的每一风电场，获取由该风电场的风力发电功率预测值组成的数据集合，将所有风电场对应的数据集合中的风力发电功率预测值相加，得到由风电场群总发电功率组成的综合数据集合；利用多项式拟合算法对综合数据集合进行拟合，得到风电平滑出力公式，并计算风电平滑出力值；计算风电平滑出力值与火力发电厂基准出力值之和，得到总平滑出力值；根据总平滑出力值、风电场群总发电功率、以及火力发电厂最小出力值、最大出力值、基准出力值，确定风电场群和火力发电厂的实际出力情况。本发明得到的平滑出力曲线不会滞后延时，且确保耗能低、污染少，输出电能平稳。

权利要求

1.一种风电场群与火力发电厂联合发电系统的平滑出力方法，其特征在于，包括：步骤A，针对风电场群中的每一风电场，获取由该风电场的风力发电功率预测值组成的数据集合，将所有风电场对应的数据集合中的风力发电功率预测值相加，得到由风电场群总发电功率组成的综合数据集合；

步骤B，利用多项式拟合算法对所述综合数据集合进行拟合，得到风电平滑出力公式，并根据所述风电平滑出力公式计算风电平滑出力值；

步骤C，计算所述风电平滑出力值与火力发电厂基准出力值之和，得到总平滑出力值；

步骤D，根据所述总平滑出力值、所述风电场群总发电功率、以及火力发电厂最小出力值、最大出力值、基准出力值，确定风电场群和火力发电厂的实际出力情况；

所述步骤A具体为：

步骤A1，针对风电场群中的每一风电场，获取由该风电场的风力发电功率预测值组成的数据集合Pj：Pj＝{(pji,ti)|j＝1,2...,k；i＝1,2...,m}步骤A2，将所有风电场对应的数据集合中的风力发电功率预测值相加，得到由风电场群总发电功率组成的综合数据集合P：P＝{(pi,ti)|i＝1,2...,m}其中，k为风电场群中的风电场总数，j为风电场序列号，m为各风电场对应的数据集合、所述综合数据集合的样本个数，i为样本序列号，Pj为第j个风电场对应的数据集合，pji为第j个风电场的风力发电功率预测值，P为综合数据集合，pi为风电场群总发电功率，ti为pji、pi对应的时间；

所述步骤B具体包括：

步骤B1，根据所述综合数据集合P中风电场群总发电功率pi的波动趋势，确定所述风电平滑出力公式的阶数n，其中n为自然数；

步骤B2，拟合具有所述阶数n的多项式：

antin+an-1tin-1+…+a1ti+a0；

其中，a0～an为多项式系数；

步骤B3，计算所述多项式antin+an-1tin-1+…+a1ti+a0与所述风电场群总发电功率pi的差值平方和Err：步骤B4，利用最小二乘法计算所述差值平方和Err为最小值时，多项式系数a0～an对应的具体值α0～αn；

步骤B5，利用所述具体值α0～αn构建风电平滑出力公式X(t)：X(t)＝αntn+αn-1tn-1+…+α1t+α0；

其中，t为时间；

步骤B6，计算当t＝ti时，所述风电平滑出力公式X(t)的值：X(ti)＝αntin+αn-1tin-1+…+α1ti+α0其中，X(ti)为风电平滑出力值；

所述步骤C具体为：

计算风电平滑出力值X(ti)与火力发电厂基准出力值之和，得到总平滑出力值：Y(ti)＝X(ti)+Pdefault

其中，Y(ti)为总平滑出力值，Pdefault为火力发电厂基准出力值；

所述步骤D具体包括：

步骤D1，计算所述总平滑出力值Y(ti)与风电场群总发电功率pi的差值Δpi：Δpi＝Y(ti)-pi；

步骤D2，令风电场群按照所述风电场群总发电功率pi出力；

步骤D3，若所述差值Δpi小于或等于火力发电厂最小出力值Pmin时，令火力发电厂按照其最小出力值Pmin出力；若所述差值Δpi大于或等于火力发电厂最大出力值Pmax时，令火力发电厂按照其最大出力值Pmax出力；若所述差值Δpi大于火力发电厂最小出力值Pmin且小于火力发电厂最大出力值Pmax时，令火力发电厂按照所述差值Δpi出力。

2.根据权利要求1所述的方法，其特征在于，所述步骤B1具体包括：根据综合数据集合P中风电场群总发电功率pi的波动趋势，确定平滑出力曲线波形；

根据所述平滑出力曲线波形，确定所述风电平滑出力公式的阶数n。

3.根据权利要求2所述的方法，其特征在于，当所述平滑出力曲线波形为直线时，确定所述风电平滑出力公式的阶数n＝1。

4.根据权利要求2所述的方法，其特征在于，当所述平滑出力曲线波形为抛物线时，确定所述风电平滑出力公式的阶数n＝2。

5.根据权利要求1所述的方法，其特征在于，所述步骤B4具体包括：分别对多项式系数a0～an求偏导数，得到如下方程组：求解以上方程组，得到多项式系数a0～an对应的具体值α0～αn。

6.根据权利要求1所述的方法，其特征在于，所述火力发电厂最小出力值Pmin、最大出力值Pmax、基准出力值Pdefault具有如下关系：

说明书

一种风电场群与火力发电厂联合发电系统的平滑出力方法

技术领域

[0001] 本发明涉及电力技术领域，具体地，涉及一种风电场群与火力发电厂联合发电系统的平滑出力方法。

背景技术

[0002] 火力发电是目前最重要、技术最为成熟的发电方式，当前火力发电在我国发电领域中仍占有非常大的比例，但火力发电存在能耗高、污染严重等缺点。

[0003] 近年来，风能等可再生清洁能源已越来越多被应用于发电领域，风电场群与火力发电厂联合发电系统是一种综合了风力发电和火力发电的发电系统，这种发电系统在风能资源充足时增加风力供电的比重，充分发挥风力发电技术能耗低、污染小的优势，在风能资源不够充足时，利用火力发电满足电网运行的需要。图1为风电场群与火力发电厂联合发电系统的结构示意图，多个风电场与一个火力发电厂接入公共电网。

[0004] 风能资源存在着随机性和波动性的问题，导致风力发电的功率具有波动性，提供给电网的电能就不够平稳，影响电网的平稳运行。为了减少这种冲击，需要对风力发电进行功率平抑，以减小功率波动性对电网的影响。

[0005] 东北电力大学宇航、严干贵等人利用一阶低通滤波算法实现风电功率波动平抑的控制策略。该控制策略主要对风电场运行中的高频分量进行滤除，减小风电功率的变化率，为电力系统提供较为稳定的功率输出，而储能系统则是通过其充放电来改变输出功率的幅值，使注入电网的电能更加平稳。但实际应用中发现这种利用一阶低通滤波算法得到的出力平滑曲线存在一定的延时作用，如图2所示，较细的线为风电单独出力曲线，较粗的线为利用这种控制策略得到的风储平滑出力曲线，从图2可明显看出，风储平滑出力曲线滞后于风电单独出力曲线。这是因为这种利用一阶低通滤波算法实现风电功率波动平抑的控制策略是采用本次采样值与上次滤波输出值进行加权得到本次滤波输出值，具体公式为：

[0006] Y(n)＝αX(n)+(1-α)Y(n-1)

[0007] 以上公式中，α为滤波系数；X(n)为本次采样值；Y(n-1)为上次滤波输出值；Y(n)为本次滤波输出值。

[0008] 可见，这种利用一阶低通滤波算法实现风电功率波动平抑的控制策略还存在不足之处。

发明内容

[0009] 本发明实施例的主要目的在于提供一种风电场群与火力发电厂联合发电系统的平滑出力方法，以解决现有技术利用一阶低通滤波算法平抑风电功率波动所得到的出力平滑曲线存在延时现象的问题。

[0010] 为了实现上述目的，本发明实施例提供一种风电场群与火力发电厂联合发电系统的平滑出力方法，包括：

[0011] 步骤A，针对风电场群中的每一风电场，获取由该风电场的风力发电功率预测值组成的数据集合，将所有风电场对应的数据集合中的风力发电功率预测值相加，得到由风电场群总发电功率组成的综合数据集合；

[0012] 步骤B，利用多项式拟合算法对所述综合数据集合进行拟合，得到风电平滑出力公式，并根据所述风电平滑出力公式计算风电平滑出力值；

[0013] 步骤C，计算所述风电平滑出力值与火力发电厂基准出力值之和，得到总平滑出力值；

[0014] 步骤D，根据所述总平滑出力值、所述风电场群总发电功率、以及火力发电厂最小出力值、最大出力值、基准出力值，确定风电场群和火力发电厂的实际出力情况；

[0015] 所述步骤A具体为：

[0016] 步骤A1，针对风电场群中的每一风电场，获取由该风电场的风力发电功率预测值组成的数据集合Pj：

[0017] Pj＝{(pji,ti)|j＝1,2...,k；i＝1,2...,m}

[0018] 步骤A2，将所有风电场对应的数据集合中的风力发电功率预测值相加，得到由风电场群总发电功率组成的综合数据集合P：

[0019] P＝{(pi,ti)|i＝1,2...,m}

[0020]

[0021] 其中，k为风电场群中的风电场总数，j为风电场序列号，m为各风电场对应的数据集合、所述综合数据集合的样本个数，i为样本序列号，Pj为第j个风电场对应的数据集合，pji为第j个风电场的风力发电功率预测值，P为综合数据集合，pi为风电场群总发电功率，ti为pji、pi对应的时间；

[0022] 所述步骤B具体包括：

[0023] 步骤B1，根据所述综合数据集合P中风电场群总发电功率pi的波动趋势，确定所述风电平滑出力公式的阶数n，其中n为自然数；

[0024] 步骤B2，拟合具有所述阶数n的多项式：

[0025]

[0026] 其中，a0～an为多项式系数；

[0027] 步骤B3，计算所述多项式 与所述风电场群总发电功率pi的差值平方和Err：

[0028]

[0029] 步骤B4，利用最小二乘法计算所述差值平方和Err为最小值时，多项式系数a0～an对应的具体值α0～αn；

[0030] 步骤B5，利用所述具体值α0～αn构建风电平滑出力公式X(t)：

[0031] X(t)＝αntn+αn-1tn-1+…+α1t+α0；

[0032] 其中，t为时间；

[0033] 步骤B6，计算当t＝ti时，所述风电平滑出力公式X(t)的值：

[0034]

[0035] 其中，X(ti)为风电平滑出力值；

[0036] 所述步骤C具体为：

[0037] 计算风电平滑出力值X(ti)与火力发电厂基准出力值之和，得到总平滑出力值：

[0038] Y(ti)＝X(ti)+Pdefault

[0039] 其中，Y(ti)为总平滑出力值，Pdefault为火力发电厂基准出力值；

[0040] 所述步骤D具体包括：

[0041] 步骤D1，计算所述总平滑出力值Y(ti)与风电场群总发电功率pi的差值Δpi：

[0042] Δpi＝Y(ti)-pi；

[0043] 步骤D2，令风电场群按照所述风电场群总发电功率pi出力；

[0044] 步骤D3，若所述差值Δpi小于或等于火力发电厂最小出力值Pmin时，令火力发电厂按照其最小出力值Pmin出力；若所述差值Δpi大于或等于火力发电厂最大出力值Pmax时，令火力发电厂按照其最大出力值Pmax出力；若所述差值Δpi大于火力发电厂最小出力值Pmin且小于火力发电厂最大出力值Pmax时，令火力发电厂按照所述差值Δpi出力。

[0045] 借助于上述技术方案，本发明通过对整个计划出力区间的风电场群总发电功率进行多项式拟合，最终得到的平滑出力曲线不会滞后延时，相比于利用一阶低通滤波算法进行风电功率波动平抑的方法，本发明具有更优化的平滑出力效果，同时，本发明将对风电场群进行功率波动平抑的结果与火力发电厂的出力加在一起作为整体来控制联合发电的实际出力情况，以确保整个系统耗能低、污染少，且提供给电网的电能平稳。

附图说明

[0046] 为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

[0047] 图1是本发明背景技术提供的风电场群与火力发电厂联合发电系统的结构示意图；

[0048] 图2是本发明背景技术提供的利用一阶低通滤波算法进行风电波动平抑前后的功率曲线对比示意图；

[0049] 图3是本发明提供的风电场群与火力发电厂联合发电系统的平滑出力方法流程示意图；

[0050] 图4是本发明提供的平滑出力前后功率曲线对比示意图。

具体实施方式

[0051] 下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0052] 本发明提供一种风电场群与火力发电厂联合发电系统的平滑出力方法，如图3所示，该方法包括：

[0053] 步骤S1，针对风电场群中的每一风电场，获取由该风电场的风力发电功率预测值组成的数据集合，将所有风电场对应的数据集合中的风力发电功率预测值相加，得到由风电场群总发电功率组成的综合数据集合。

[0054] 步骤S2，利用多项式拟合算法对综合数据集合进行拟合，得到风电平滑出力公式，并根据风电平滑出力公式计算风电平滑出力值。

[0055] 步骤S3，计算风电平滑出力值与火力发电厂基准出力值之和，得到总平滑出力值。

[0056] 步骤S4，根据总平滑出力值、风电场群总发电功率、以及火力发电厂最小出力值、最大出力值、基准出力值，确定火力发电厂的实际出力情况。

[0057] 下面以公式形式对上述方法的具体过程进行说明：

[0058] 步骤S1具体包括如下两步骤：

[0059] 步骤S11，针对风电场群中的每一风电场，获取由该风电场的风力发电功率预测值组成的数据集合Pj：

[0060] Pj＝{(pji,ti)|j＝1,2...,k；i＝1,2...,m}

[0061] 该步骤可以从每一风电场的功率预测系统SCADA中获取风力发电功率预测值，即功率预测系统SCADA对该风电场在计划出力区间单独出力预测的功率值。

[0062] 步骤S12，将所有风电场对应的数据集合中的风力发电功率预测值相加，得到由风电场群总发电功率组成的综合数据集合P：

[0063] P＝{(pi,ti)|i＝1,2...,m}

[0064]

[0065] 其中，k为风电场群中的风电场总数，j为风电场序列号，m为各风电场对应的数据集合、综合数据集合的样本个数，i为样本序列号，Pj为第j个风电场对应的数据集合，pji为第j个风电场的风力发电功率预测值，P为综合数据集合，pi为风电场群总发电功率，ti为pji、pi对应的时间。

[0066] 其中，风电场群总发电功率pi为ti时刻风电场群中各个风电场的风力发电功率预测值pji之和，综合数据集合P包含了整个计划出力区间的风电场群总发电功率pi。

[0067] 步骤S2具体包括如下各步骤：

[0068] 步骤S21，根据综合数据集合P中风电场群总发电功率pi的波动趋势，确定风电平滑出力公式的阶数n，其中n为自然数。

[0069] 较佳的，该步骤按照如下具体步骤执行：

[0070] 步骤S211，根据综合数据集合P中风电场群总发电功率pi的波动趋势，确定平滑出力曲线波形；

[0071] 步骤S212，根据平滑出力曲线波形，确定风电平滑出力公式的阶数n。

[0072] 例如，当平滑出力曲线波形为直线时，确定风电平滑出力公式的阶数n＝1；当平滑出力曲线波形为抛物线时，确定风电平滑出力公式的阶数n＝2。

[0073] 步骤S22，拟合具有阶数n的多项式：

[0074] antin+an-1tin-1+…+a1ti+a0；

[0075] 其中，a0～an为多项式系数。

[0076] 步骤S23，计算多项式antin+an-1tin-1+…+a1ti+a0与风电场群总发电功率pi的差值平方和Err：

[0077]

[0078] 步骤S24，利用最小二乘法计算差值平方和Err为最小值时，多项式系数a0～an对应的具体值α0～αn。

[0079] 较佳的，该步骤按照采用如下方式计算：

[0080] 分别对多项式系数a0～an求偏导数，得到如下方程组：

[0081]

[0082] 求解以上方程组，得到多项式系数a0～an对应的具体值α0～αn。

[0083] 步骤S25，利用具体值α0～αn构建风电平滑出力公式X(t)：

[0084] X(t)＝αntn+αn-1tn-1+…+α1t+α0；

[0085] 其中，t为时间；

[0086] 步骤S26，计算当t＝ti时，风电平滑出力公式X(t)的值：

[0087]

[0088] 其中，X(ti)为风电平滑出力值。

[0089] 具体的，该步骤利用多项式拟合算法对整个计划出力区间的风电场群总发电功率数据进行拟合，得到的平滑出力值是对风电场群的功率波动进行平抑后的结果，由于不是如现有的利用一阶低通滤波算法那样采用相邻滤波输出值来计算当前的滤波输出值，因此本发明得到的平滑出力曲线(平滑出力公式对应的曲线)不会滞后延时，平滑效果更加优化。

[0090] 图4中虚线所示为某风电场群在整个计划出力区间的风电场群总发电功率pi组成的曲线，利用本发明提供的方法进行平滑出力后得到如图4中实线所示的平滑出力曲线，通过对比可知，平滑出力后得到的平滑出力曲线减小了功率波动，且不存在延时现象。

[0091] 由于本发明研究的是由风电场群与火力发电厂共同组成的联合发电系统，该系统应尽量降低火力发电而使用风力发电，以保证耗能低、污染少，但火力发电厂的最小出力又不可能降为零，这就需要将对风电场群进行功率波动平抑的结果(平滑出力值)与火力发电厂的出力加在一起作为整体来控制联合发电的实际出力情况，以确保整个系统耗能低、污染少，且提供给电网的电能平稳。

[0092] 步骤S3具体为：

[0093] 计算风电平滑出力值X(ti)与火力发电厂基准出力值之和，得到总平滑出力值：

[0094] Y(ti)＝X(ti)+Pdefault

[0095] 其中，Y(ti)为总平滑出力值，Pdefault为火力发电厂基准出力值，是火力发电厂默认情况下的出力值，该基准出力值可根据火力发电厂的具体情况(需要考虑火力发电厂的机组特性)设定。

[0096] 步骤S4具体包括如下步骤：

[0097] 步骤S41，计算总平滑出力值Y(ti)与风电场群总发电功率pi的差值Δpi：

[0098] Δpi＝Y(ti)-pi；

[0099] 步骤S42，令风电场群按照风电场群总发电功率pi出力；

[0100] 步骤S43，若差值Δpi小于或等于火力发电厂最小出力值Pmin时，令火力发电厂按照其最小出力值Pmin出力；

[0101] 若差值Δpi大于或等于火力发电厂最大出力值Pmax时，令火力发电厂按照其最大出力值Pmax出力；

[0102] 若差值Δpi大于火力发电厂最小出力值Pmin且小于火力发电厂最大出力值Pmax时，令火力发电厂按照差值Δpi出力。

[0103] 以上步骤中，火力发电厂最大出力值Pmax和最小出力值Pmin分别是火力发电厂能输出的最大功率和最小功率，可根据火力发电厂的具体情况(需要考虑火力发电厂的机组特性)设定。一般情况下，Pmin＜Pdefault＜Pmax。

[0104] 较佳的，火力发电厂最小出力值Pmin、最大出力值Pmax、基准出力值Pdefault具有如下关系：

[0105]

[0106] 例如，当火力发电厂最小出力值Pmin＝500MW，最大出力值Pmax＝1000MW，基准出力值Pdefault＝750MW。

[0107] 通过步骤4，整个风电场群与火力发电厂联合发电系统将按照如下策略出力：

[0108] 当总平滑出力值Y(ti)与风电场群总发电功率pi的差值Δpi小于或等于火力发电厂最小出力值Pmin时，说明火力发电厂只要按照其最小功率输出电能，再加上与风电场群输出的电能之和，就可以满足总平滑出力值；

[0109] 当总平滑出力值Y(ti)与风电场群总发电功率pi的差值Δpi大于或等于火力发电厂最大出力值Pmax时，说明火力发电厂需要按照其最大功率输出电能，再加上与风电场群输出的电能之和，才能基本满足总平滑出力值；

[0110] 当总平滑出力值Y(ti)与风电场群总发电功率pi的差值Δpi大于火力发电厂最小出力值Pmin且小于火力发电厂最大出力值Pmax时，说明火力发电厂适当出力(差值Δpi)，再加上与风电场群输出的电能之和，即可满足总平滑出力值；

[0111] 以上各种情况下，风电场群都是按照风电场群总发电功率pi出力，即按照其实际能输出的最大功率出力，其主要原因是：风电平滑出力值X(ti)是对风电场群总发电功率pi进行波动平抑后的结果，X(ti)与pi相差并不大，而总平滑出力值Y(ti)＝X(ti)+Pdefault，火力发电厂基准出力值Pdefault＞0，因此，总平滑出力值Y(ti)总会大于风电场群总发电功率pi，这种情况下，单独靠风电场群输出电能(即便按照其实际能输出的最大功率出力)并不能满足总平滑出力值Y(ti)的要求，还需要火力发电厂输出电能来补充，基于耗能低、污染少的原则，尽量降低火力发电的比例，就需要风电场群按照其实际能输出的最大功率出力，即按照风电场群总发电功率pi出力。

[0112] 以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。