CN105894134A - 最小化弃风弃光电量的新能源发电权置换电量评估方法 - Google Patents

Abstract

本发明公开了一种最小化弃风弃光电量的新能源发电权置换电量评估方法，通过对风电场、光伏电站和火电厂日前发电权置换进行评估。从而最大化新能源发电量，最小化弃风及弃光。并简化交易电量滚动修正的过程，只进行日前的置换电量评估，简化计算过程。从而达到提高新能源利用率以及提高发电企业利润的目的。

Description

最小化弃风弃光电量的新能源发电权置换电量评估方法

技术领域

本发明涉及新能源发电量运行评估技术，具体地，涉及一种最小化弃风弃光电量的新能源发电权置换电量评估方法。

背景技术

发电权是发电厂商发电运行并将电能出售给电网或通过电网出售给终端用户的权利。在管制电力系统中，它对应于发电厂商在特定机制下获得的发电配额；而在自由竞争的电力市场中，它是指发电厂商在远期电力市场或当前电力市场中通过报价竞标获得的电量。

发电权置换是指上述配额或中标电量的交换。发电权置换本质上是电力系统在调度、发电资源重新配置的过程。积极合理的发电权置换能够使经济高效的发电形式替代高能耗、高污染的发电形式承担负荷，提升电力系统运行的经济性和可靠性、保护环境、促进可持续发展，同时使参与发电权置换的发电集团获益。发电权置换交易完成后，各发电厂的出力计划会相应调整，这也意味着电力系统运行方式的改变。

发电权置换的研究背景是一方面，当前我国部分地区尤其是“三北”风电、光伏等新能源装机增长较快，用电负荷增长较慢，造成大量弃风、弃光；另一方面，由于火电厂消耗大量煤炭等非可再生资源，同时造成环境污染严重，雾霾天气频发。发电权置换的主要对象是常规电源和新能源企业，通过将常规电源的发电量置换成新能源电源的发电量，从而实现降低发电成本、提高企业利润，同时节能减排、保护环境。发电权置换的形式包括中长期发电权置换、短期发电权置换、以及日前发电权置换等置换方式，需要指出的是日前发电权置换需要考虑线路传输容量等运行问题。

目前，发电权置换在发电权重新配置，特别是日前电力系统运行方式调整等方面的效益并未得到充分显示。已开展的相关研究也相对较少。

专利201510375285.1提出了一种“一种促进风电消纳的多时间尺度发电权交易方法”，研究的是不同发电厂商之间的发电权交易，涉及电量交易和资金交易，研究目标是最大化社会效益总和，并通过逐小时修正发电权交易计划最终确定发电权交易电量，因为风电预测的不确定性、负荷预测的不确定性以及机组开停机的变化性等因素，该过程本身就是在一定概率意义下的评估，进行逐小时滚动修正会增加大量的工作量，其效果并不显著。

发明内容

本发明的目的在于，针对上述问题，提出一种最小化弃风弃光电量的新能源发电权置换电量评估方法，以实现提高新能源利用率以及提高发电企业利润的优点。

为实现上述目的，本发明采用的技术方案是：

一种最小化弃风弃光电量的新能源发电权置换电量评估方法，

对于同一发电集团旗下的风电场、光伏电站与火电厂进行日前发电权置换，风电场与火电厂在t时刻的发电权置换功率为Z_w,h(t)，光伏电站与火电厂在t时刻的发电权置换功率为Z_p,h(t)，t为时间，t∈[0,24)；

通过发电权置换最小化弃风弃光电量，即最大化发电权置换电量的目标函数为：

$\begin{array}{cc}max& Q_{w,h}=\underset{0}{\overset{T}{\∫}}\left({\mathrm{\αZ}}_{w,h}\right(t)+{\mathrm{\βZ}}_{p,h}(t\left)\right)dt\end{array}$

其中，T为24；

该方法具体包括以下步骤：

步骤1：设置新能源出力下达曲线是否经过校验标志位为False，所述新能源包括风电场和光伏电站；

步骤2：确定风电场出力上下限约束以及光伏电站的出力上下限约束：

所述风电场出力上下限约束为：

$\underset{i=1}{\overset{M}{\Σ}}{P}_{w,i}^{\min }{U}_{w,i}\left(t\right)\≤P_w\left(t\right)+Z_{w,h}\left(t\right)\≤\underset{i=1}{\overset{M}{\Σ}}{P}_{w,i}^{\max }{U}_{w,i}\left(t\right)$

其中，P_w(t)为风电场在t时刻的原计划上网出力，为第i台风电机组的出力下限，为第i台风电机组的出力上限，U_w,i(t)为第i台风电机组在t时刻的启停状态，i＝1,2,3,...,M，M为风电场的风电机组台数；

所述光伏电站的出力上下限约束为：

$\underset{i=1}{\overset{K}{\Σ}}{P}_{p,i}^{\min }{U}_{p,i}\left(t\right)\≤P_p\left(t\right)+Z_{p,h}\left(t\right)\≤\underset{i=1}{\overset{K}{\Σ}}{P}_{p,i}^{\max }{U}_{p,i}\left(t\right)$

其中，P_p(t)为光伏电站在t时刻的原计划上网出力，为第i台光伏逆变器的出力下限，为第i台光伏逆变器的出力上限，U_p,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,K，K为光伏电站的光伏逆变器台数；

步骤3：确定火电厂出力上下限约束：

$\underset{i=1}{\overset{N}{\Σ}}{P}_{h,i}^{\min }{U}_{h,i}\left(t\right)\≤P_h\left(t\right)-Z_{w,h}\left(t\right)\≤\underset{i=1}{\overset{N}{\Σ}}{P}_{h,i}^{\max }{U}_{h,i}\left(t\right),$

其中，P_h(t)为火电厂在t时刻的原计划上网出力，为第i台火电机组的出力下限，为第i台火电机组的出力上限，U_h,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂的火电机组台数；

步骤4：确定调度部门给风电场预留的送出通道容量约束P_w,c(t)，以及确定调度部门给光伏电站预留的送出通道容量约束P_p,c(t)；

步骤5：确定火电厂的最小经济技术出力约束t∈[0,24)；

步骤6：确定风电场的日前风电功率预测曲线；

步骤7：确定电源-负荷有功功率平衡约束；

步骤8：生成发电权置换后风电场和光伏电站的日前出力下达计划；

步骤9：判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足风电机组及光伏逆变器最小启停时间约束，如果是，进入步骤11，如果否，进入步骤10；

步骤10：增加风电机组及光伏逆变器启停限制，将不满足启停约束的机组的启停状态设置为启动并不得停止，转步骤9；

步骤11：判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足火电机组最小启停时间约束，如果是，进入步骤13，如果否，进入步骤12；

步骤12：增加火电机组启停限制，将不满足启停约束的机组的启停状态设置为启动并不得停止，转步骤11；

步骤13：判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足风电机组及光伏逆变器爬坡速率约束，如果是，进入步骤15，如果否，进入步骤14；

步骤14：增加风电机组、光伏电站爬坡率限制，从不满足机组爬坡率的时刻开始，逐点减少发电权置换功率，并转步骤13；

步骤15：判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足火电机组爬坡速率约束，如果是，进入步骤17，如果否，进入步骤16；

步骤16：增加火电机组爬坡率限制，从不满足机组爬坡率的时刻开始，逐点减少发电权置换功率，并转步骤15；

步骤17：风电出力下达曲线是否经过校验，如果是，转步骤19；如果否，转步骤18；

步骤18：设置风电出力下达曲线是否经过校验标志位为True，并转步骤2；

步骤19：得到经过校验的新能源发电权置换后的日出力计划曲线；

步骤20：计算风电场及光伏电站通过日前发电权置换所减少的弃风弃光电量；

步骤21：计算发电集团增加的利润；

步骤22：计算发电权置换社会效益；

步骤23：完成最小化弃风弃光电量的新能源发电权置换电量评估。

优选的，所述步骤4中确定调度部门给风电场预留的送出通道容量约束P_w,c(t)具体为：

根据风电场发电规则，有

P_w(t)+Z_w,h(t)≤P_w,c(t)，

其中，P_w(t)为风电计划出力值，t为时间，t∈[0,24)，Z_w,t(t)为风电场通过发电权置换得到的功率，P_w,c(t)为调度部门给风电场预留的风电送出通道容量；

所述步骤4中确定调度部门给光伏电站预留的送出通道容量约束P_p,c(t)具体为：

根据光伏电站发电规则，则有：

P_p(t)+Z_p,h(t)≤P_p,c(t)，

其中，P_p(t)为光伏电站计划出力值，t为时间，t∈[0,24)，P_p,c(t)为调度部门给光伏电站预留的风电送出通道容量。

优选的，所述步骤6中确定风电场的日前风电功率预测曲线具体为：

首先利用风资源监测网络的测风塔的实测风速，修正数值天气预报的预测数据，用修正后的数值天气预报的预测数据结合风电场的日前开机计划通过风电功率预测系统预测明天0-24小时内风电场的出力曲线，取置信水平α的风电预测条带的下限

确定光伏电站的日前光伏功率预测曲线，首先利用光伏资源监测网络的测光站的实测辐照度，修正数值天气预报的预测数据，用修正后的数值天气预报的预测数据结合光伏电站的日前开机计划通过光伏功率预测系统预测明天0-24小时内光伏电站的出力曲线，取置信水平α的光伏预测条带的下限

优选的，所述步骤7中确定电源-负荷有功功率平衡约束具体为：

根据负荷预测结果，得到置信水平α的负荷预测条带的下限再结合火电厂的最小经济技术出力曲线，确定给风电场和光伏电站预留的出力曲线

$P_w^{up}\left(t\right)+P_p^{up}\left(t\right)=P_{d,\mathrm{\α}}^{down}\left(t\right)-P_h^{\min }\left(t\right)$

为功率平衡约束下风电的出力上限，为功率平衡约束下光伏的出力上限；

风电场和光伏电站的出力分配遵循调度优先级排序原则，即根据预先设定的优先级，选择优先调度风电还是优先调度光伏；

即根据电源-负荷平衡约束压制风电，让光伏电站尽量满发，

据此，则有：

$P_p^{up}\left(t\right)=P_{d,\mathrm{\α}}^{down}\left(t\right)-P_h^{\min }\left(t\right)$

$P_w^{up}\left(t\right)=P_{d,\mathrm{\α}}^{down}\left(t\right)-P_h^{min}\left(t\right)-P_{p,\mathrm{\α}}^{up}\left(t\right)$

其中，为功率平衡约束下光伏的出力上限，为置信水平α的负荷预测条带的下限，为火电厂的最小经济技术出力。为发电权置换后的光伏出力。

优选的，所述步骤8：生成发电权置换后风电场和光伏电站的日前出力下达计划中，光伏电站应满足约束为：

$P_{p,\mathrm{\α}}^{up}\left(t\right)=P_p\left(t\right)+Z_{p,h}\left(t\right)=min\{P_{p,c}\left(t\right),P_p^{up}\left(t\right),P_{p,\mathrm{\α}}^{down}\left(t\right)\},$

其中，P_p,c(t)为调度部门给光伏电站预留的送出通道容量约束，为风电预测条带的下限，为对光伏电站的电源-负荷功率平衡约束；

风电场应满足约束为：

$P_{w,\mathrm{\α}}^{up}\left(t\right)=P_w\left(t\right)+Z_{w,h}\left(t\right)=min\{P_{w,c}\left(t\right),P_w^{up}\left(t\right),P_{w,\mathrm{\α}}^{down}\left(t\right)\},$

其中，P_w,c(t)为调度部门给风电场预留的送出通道容量约束，为风电预测条带的下限，为对风电场的电源-负荷功率平衡约束。

优选的，所述步骤9判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足风电机组及光伏逆变器最小启停时间约束中，风电机组最小启停时间约束为：

$\left\{\begin{array}{c}{T}_{w,i}^{on}\≥T_{w,i}^S\\ T_{w,i}^{off}\≥T_{w,i}^D\end{array}\right.,$

其中，为第i台风电机组的连续运行时间及连续停止时间，为第i台风电机组的最小开停机时间；

光伏逆变器最小启停时间约束为：

$\left\{\begin{array}{c}{T}_{p,i}^{on}\≥T_{p,i}^S\\ T_{p,i}^{off}\≥T_{p,i}^D\end{array}\right.$

其中，为第i台光伏逆变器的连续运行时间及连续停止时间，为第i台光伏逆变器的最小开停机时间。

优选的，所述步骤11判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足火电机组最小启停时间中，火电机组最小启停时间约束具体为：

$\left\{\begin{array}{c}{T}_{h,i}^{on}\≥T_{h,i}^S\\ T_{h,i}^{off}\≥T_{h,i}^D\end{array}\right.,$

其中，为第i台火电机组的连续运行时间及连续停止时间， 为第i台风电机组的最小开停机时间。

优选的，所述步骤13中判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足风电机组及光伏逆变器爬坡速率约束，风电机组爬坡速率约束为：

$\begin{array}{c}-\underset{i=1}{\overset{M}{\Σ}}{R}_{w,i}^{down}{U}_{w,i}\left(t\right)\mathrm{\Δ}t\≤\left(P_w\right(t)+Z_{w,h}(t\left)\right)-\\ \left(P_w\right(t-\mathrm{\Δ}t)+Z_{w,h}(t-\mathrm{\Δ}t\left)\right)\≤\underset{i=1}{\overset{M}{\Σ}}{R}_{w,i}^{up}{U}_{w,i}\left(t\right)\mathrm{\Δ}t\end{array},$

其中，P_w(t)为风电场在t时刻的原计划上网出力，Z_w,h(t)为风电场通过发电权置换在t时刻得到的功率，为第i台风电机组的爬坡率下限，为第i台风电机组的爬坡率上限，U_w,i(t)为第i台风电机组在t时刻的启停状态，i＝1,2,3,...,M，M为风电场的风电机组台数；

光伏逆变器爬坡速率约束为：

$\begin{array}{c}-\underset{i=1}{\overset{K}{\Σ}}{R}_{p,i}^{down}{U}_{p,i}\left(t\right)\mathrm{\Δ}t\≤\left(P_p\right(t)+Z_{p,h}(t\left)\right)-\\ \left(P_p\right(t-\mathrm{\Δ}t)+Z_{p,h}(t-\mathrm{\Δ}t\left)\right)\≤\underset{i=1}{\overset{K}{\Σ}}{R}_{p,i}^{up}{U}_{p,i}\left(t\right)\mathrm{\Δ}t\end{array},$

P_p(t)为光伏电站在t时刻的原计划上网出力，Z_p,h(t)为光伏电站通过发电权置换在t时刻得到的功率，为第i台光伏逆变器的爬坡速率下限，为第i台光伏逆变器的爬坡上限，U_p,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,K，K为光伏电站的光伏逆变器台数；

所述步骤15中，火电机组爬坡速率约束为：

$\begin{array}{c}-\underset{i=1}{\overset{N}{\Σ}}{R}_{h,i}^{down}{U}_{h,i}\left(t\right)\mathrm{\Δ}t\≤\left(P_h\right(t)-Z_{w,h}(t\left)\right)-\\ \left(P_h\right(t-\mathrm{\Δ}t)-Z_{w,h}(t-\mathrm{\Δ}t\left)\right)\≤\underset{i=1}{\overset{N}{\Σ}}{R}_{h,i}^{up}{U}_{h,i}\left(t\right)\mathrm{\Δ}t\end{array}$

其中，P_h(t)为火电厂在t时刻的正常上网出力，Z_w,h(t)为火电厂发电权置换掉的功率，为第i台火电机组的上下爬坡率，U_h,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂B的火电机组台数。

优选的，所述步骤20计算风电场及光伏电站通过日前发电权置换所减少的弃风弃光电量中：

减少的弃风电量为：

$Q_{w,h}=\underset{0}{\overset{T}{\∫}}{Z}_{w,h}\left(t\right)dt=\underset{0}{\overset{T}{\∫}}\left(P_{w,\mathrm{\α}}^{up}\right(t)-P_w(t\left)\right)dt,$

减少的弃光电量为：

$Q_{p,h}=\underset{0}{\overset{T}{\∫}}{Z}_{p,h}\left(t\right)dt=\underset{0}{\overset{T}{\∫}}\left(P_{p,\mathrm{\α}}^{up}\right(t)-P_p(t\left)\right)dt,$

其中，Z_w,h(t)为风电场通过发电权置换得到的功率，T为24，表示电量统计的时间为明天一天，为发电权置换后风电场的日前风电出力，P_w(t)为风电原计划出力值，

Z_p,h(t)为光伏电站通过发电权置换得到的功率，T为24，表示电量统计的时间为明天一天，为发电权置换后光伏电站的日前出力，P_p(t)为光伏原计划出力值。

优选的，所述步骤21中计算发电集团增加的利润具体为：

$C_w^z\left(Z_{w,h}\right)=\left({\mathrm{\λ}}_{w,b}\right(t)+{\mathrm{\λ}}_z(t\left)\right)\underset{0}{\overset{T}{\∫}}{Z}_{w,h}\left(t\right)dt+\left({\mathrm{\λ}}_{p,b}\right(t)+{\mathrm{\λ}}_z(t\left)\right)\underset{0}{\overset{T}{\∫}}{Z}_{p,h}\left(t\right)dt$

其中，Z_w,h(t)为t时刻风电发电权置换的电量，Z_p,h(t)为t时刻光伏发电权置换的电量，λ_w,b(t)为国家每度风电的新能源补贴，λ_z(t)为火电机组的脱硫标杆上网电价，λ_p,b(t)为国家每度光伏发电的新能源补贴，T＝24表示一天的置换电量所增加的利润；

所述步骤22中计算发电权置换社会效益具体为：

$S_w^z\left(Z_{w,h}\right)=a\mathrm{\η}\underset{0}{\overset{T}{\∫}}\left(Z_{w,h}\right(t)+Z_{p,h}(t\left)\right)dt+b\mathrm{\γ}\underset{0}{\overset{T}{\∫}}\left(Z_{w,h}\right(t)+Z_{p,h}(t\left)\right)dt$

其中，Z_w,h(t)为t时刻风电发电权置换的电量，Z_p,h(t)为t时刻光伏发电权置换的电量，η为火电厂每度电需要消耗的标准煤，γ为火电厂每度电产生的污染物，a、b为加权系数。

本发明的技术方案具有以下有益效果：

本发明的技术方案同一发电集集团旗下不同发电形式的发电厂之间(如火电厂和风电场、火电厂和光伏电站)的发电权置换，只涉及电量置换，不涉及资金交易，为了最大化新能源发电量，最小化弃风、弃光从而简化交易电量滚动修正的过程，只进行日前的置换电量评估，优化简化计算过程，达到提高新能源利用率以及提高发电企业利润的目的。

下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。

附图说明

图1为本发明实施例所述的最小化弃风弃光电量的新能源发电权置换电量评估方法的流程图；

图2为风电场的日前风电功率预测曲线和光伏电站的日前光伏功率预测曲线示意图；

图3为根据典型日负荷曲线计算出由于系统平衡约束给风电场及光伏电站预留的出力曲线示意图；

图4为经发电权置换后的风电出力计划下达曲线示意图。

结合附图，本发明实施例中附图标记如下：

1-负荷预测条带下限；2-给风电场的预留出力；3-风电送出容量约束；4-火电最小经济技术出力；5-光伏送出容量约束；6-带约束的光伏预测出力；7-光伏出力预测。

具体实施方式

以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。

如图1所示，一种最小化弃风弃光电量的新能源发电权置换电量评估方法，

对于同一发电集团旗下的风电场、光伏电站与火电厂进行日前发电权置换，风电场与火电厂在t时刻的发电权置换功率为Z_w,h(t)，光伏电站与火电厂在t时刻的发电权置换功率为Z_p,h(t)，t为时间，t∈[0,24)；

通过发电权置换最小化弃风弃光电量，即最大化发电权置换电量的目标函数为：

$\begin{array}{cc}max& Q_{w,h}=\underset{0}{\overset{T}{\∫}}\left({\mathrm{\αZ}}_{w,h}\right(t)+{\mathrm{\βZ}}_{p,h}(t\left)\right)dt\end{array}$

其中，T为24；

该方法具体包括以下步骤：

步骤1：设置新能源出力下达曲线是否经过校验标志位为False，所述新能源包括风电场和光伏电站；

步骤2：确定风电场出力上下限约束以及光伏电站的出力上下限约束：

所述风电场出力上下限约束为：

$\underset{i=1}{\overset{M}{\Σ}}{P}_{w,i}^{\min }{U}_{w,i}\left(t\right)\≤P_w\left(t\right)+Z_{w,h}\left(t\right)\≤\underset{i=1}{\overset{M}{\Σ}}{P}_{w,i}^{\max }{U}_{w,i}\left(t\right)$

其中，P_w(t)为风电场在t时刻的原计划上网出力，为第i台风电机组的出力下限，为第i台风电机组的出力上限，U_w,i(t)为第i台风电机组在t时刻的启停状态，i＝1,2,3,...,M，M为风电场的风电机组台数；

所述光伏电站的出力上下限约束为：

$\underset{i=1}{\overset{K}{\Σ}}{P}_{p,i}^{\min }{U}_{p,i}\left(t\right)\≤P_p\left(t\right)+Z_{p,h}\left(t\right)\≤\underset{i=1}{\overset{K}{\Σ}}{P}_{p,i}^{\max }{U}_{p,i}\left(t\right)$

其中，P_p(t)为光伏电站在t时刻的原计划上网出力，为第i台光伏逆变器的出力下限，为第i台光伏逆变器的出力上限，U_p,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,K，K为光伏电站的光伏逆变器台数；

步骤3：确定火电厂出力上下限约束：

$\underset{i=1}{\overset{N}{\Σ}}{P}_{h,i}^{\min }{U}_{h,i}\left(t\right)\≤P_h\left(t\right)-Z_{w,h}\left(t\right)\≤\underset{i=1}{\overset{N}{\Σ}}{P}_{h,i}^{\max }{U}_{h,i}\left(t\right),$

其中，P_h(t)为火电厂在t时刻的原计划上网出力，为第i台火电机组的出力下限，为第i台火电机组的出力上限，U_h,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂的火电机组台数；

步骤4：确定调度部门给风电场预留的送出通道容量约束P_w,c(t)，以及确定调度部门给光伏电站预留的送出通道容量约束P_p,c(t)；

步骤5：确定火电厂的最小经济技术出力约束t∈[0,24)；

步骤6：确定风电场的日前风电功率预测曲线；

步骤7：确定电源-负荷有功功率平衡约束；

步骤8：生成发电权置换后风电场和光伏电站的日前出力下达计划；

步骤9：判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足风电机组及光伏逆变器最小启停时间约束，如果是，进入步骤11，如果否，进入步骤10；

步骤10：增加风电机组及光伏逆变器启停限制，将不满足启停约束的机组的启停状态设置为启动并不得停止，转步骤9；

步骤11：判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足火电机组最小启停时间约束，如果是，进入步骤13，如果否，进入步骤12；

步骤12：增加火电机组启停限制，将不满足启停约束的机组的启停状态设置为启动并不得停止，转步骤11；

步骤13：判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足风电机组及光伏逆变器爬坡速率约束，如果是，进入步骤15，如果否，进入步骤14；

步骤14：增加风电机组、光伏电站爬坡率限制，从不满足机组爬坡率的时刻开始，逐点减少发电权置换功率，并转步骤13；

步骤15：判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足火电机组爬坡速率约束，如果是，进入步骤17，如果否，进入步骤16；

步骤16：增加火电机组爬坡率限制，从不满足机组爬坡率的时刻开始，逐点减少发电权置换功率，并转步骤15；

步骤17：风电出力下达曲线是否经过校验，如果是，转步骤19；如果否，转步骤18；

步骤18：设置风电出力下达曲线是否经过校验标志位为True，并转步骤2；

步骤19：得到经过校验的新能源发电权置换后的日出力计划曲线；

步骤20：计算风电场及光伏电站通过日前发电权置换所减少的弃风弃光电量；

步骤21：计算发电集团增加的利润；

步骤22：计算发电权置换社会效益；

步骤23：完成最小化弃风弃光电量的新能源发电权置换电量评估。

下文结合具体应用对技术方案进行具体说明：

考虑一个20.1万千瓦装机的典型风电场、一个5万装机的光伏电站与一个30万千瓦装机的典型火电厂进行发电权置换，风电场有134台1.5MW风电机组，光伏电站有100台500kW的逆变器，火电厂有1台300MW火电机组，根据上述计算流程。

目标函数：通过发电权置换最小化弃风弃光电量，即最大化发电权置换电量

$\begin{array}{cc}max& Q_{w,h}=\underset{0}{\overset{T}{\∫}}\left({\mathrm{\αZ}}_{w,h}\right(t)+{\mathrm{\βZ}}_{p,h}(t\left)\right)dt\end{array}$

其中，Z_w,h(t)为风电场通过发电权置换得到的功率，Z_p,h(t)为光伏电站通过发电权置换得到的功率，T为24，表示电量统计的时间为明天一天。

步骤1：设置新能源出力下达曲线是否经过校验标志位为False。

步骤2：确定风电场A出力上下限约束：

$\underset{i=1}{\overset{M}{\Σ}}{P}_{w,i}^{\min }{U}_{w,i}\left(t\right)\≤P_w\left(t\right)+Z_{w,h}\left(t\right)\≤\underset{i=1}{\overset{M}{\Σ}}{P}_{w,i}^{\max }{U}_{w,i}\left(t\right)$

其中，P_w(t)为风电场A在t时刻的原计划上网出力，为第i台风电机组的出力下限，为第i台风电机组的出力上限，U_w,i(t)为第i台风电机组在t时刻的启停状态，i＝1,2,3,...,M，M为风电场A的风电机组台数。原计划上网出力根据调度部门下发的日前计划结合天前风电功率预测结果进行确定，启停状态根据风电场日前开机计划获得。

确定光伏电站B的出力上下限约束：

$\underset{i=1}{\overset{K}{\Σ}}{P}_{p,i}^{\min }{U}_{p,i}\left(t\right)\≤P_p\left(t\right)+Z_{p,h}\left(t\right)\≤\underset{i=1}{\overset{K}{\Σ}}{P}_{p,i}^{\max }{U}_{p,i}\left(t\right)$

其中，P_p(t)为光伏电站B在t时刻的原计划上网出力，Z_p,h(t)为光伏电站通过发电权置换得到的功率，为第i台光伏逆变器的出力下限，为第i台光伏逆变器的出力上限，U_p,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,K，K为光伏电站B的光伏逆变器台数。原计划上网出力根据调度部门下发的日前计划结合天前光伏发电功率预测结果进行确定，启停状态根据光伏电站日前开机计划获得。

根据上述约束，风电正常上网出力+置换出力应大于0，小于所有开机运行的风电机组总容量，根据当前运行情况，明日风电场有4台机组检修停运，因此风电场正常上网出力+置换出力应小于130*1.5MW＝195MW。

根据上述约束，光伏正常上网出力+置换出力应大于0，小于所有运行光伏逆变器的总容量，根据当前运行情况，明日光伏电站有10台光伏逆变器检修停运，因此光伏电站正常上网出力+置换出力应小于40MW。

步骤3：确定火电厂B出力上下限约束：

$\underset{i=1}{\overset{N}{\Σ}}{P}_{h,i}^{\min }{U}_{h,i}\left(t\right)\≤P_h\left(t\right)-Z_{w,h}\left(t\right)\≤\underset{i=1}{\overset{N}{\Σ}}{P}_{h,i}^{\max }{U}_{h,i}\left(t\right)$

其中，P_h(t)为火电厂B在t时刻的正常上网出力，Z_w,h(t)为火电厂发电权置换掉的功率，为第i台火电机组的出力下限，为第i台火电机组的出力上限，U_h,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂B的火电机组台数。火电厂原计划上网出力由调度部门下发的日前计划确定。

根据上述约束，火电机组正常上网出力-置换出力应大于0，小于300MW。

步骤4：确定调度部门给风电场A预留的送出通道容量约束P_w,c(t)，根据风电场发电规则，则有

P_w(t)+Z_w,h(t)≤P_w,c(t)

其中，P_w(t)为风电计划出力值，t为时间，本发明研究日前置换电量评估，则有t∈[0,24)，P_w,c(t)为调度部门给风电场A预留的风电送出通道容量。

确定调度部门给光伏电站B预留的送出通道容量约束P_p,c(t)，根据光伏电站发电规则，则有

P_p(t)+Z_p,h(t)≤P_p,c(t)

其中，P_p(t)为光伏电站计划出力值，t为时间，本发明研究日前置换电量评估，则有t∈[0,24)，P_p,c(t)为调度部门给光伏电站B预留的风电送出通道容量。

调度部门为风电场明日预留的送出通道容量为180MW。风电场正常上网出力+置换出力应小于180MW。

调度部门为光伏电站明日预留的送出通道容量为30MW。光伏电站正常上网出力+置换出力应小于30MW。

步骤5：确定火电厂B的最小经济技术出力约束t∈[0,24)。

目前，新型火电机组最小技术出力已经可以达到50％额定功率，因此，本算例中为150MW。

步骤6：确定风电场A的日前风电功率预测曲线。首先利用风资源监测网络的测风塔的实测风速，修正数值天气预报(NWP)的预测数据，用修正后的NWP数据结合风电场的日前开机计划通过风电功率预测系统预测明天0-24小时内风电场A的出力曲线，取置信水平α(一般可以取α＝90％)的风电预测条带的下限

确定光伏电站B的日前光伏功率预测曲线。首先利用光伏资源监测网络的测光站的实测辐照度，修正数值天气预报(NWP)的预测数据，用修正后的NWP数据结合光伏电站的日前开机计划通过光伏功率预测系统预测明天0-24小时内光伏电站B的出力曲线，取置信水平α(一般可以取α＝90％)的光伏预测条带的下限

风电场A的日前风电功率预测曲线和光伏电站B的日前光伏功率预测曲线如图2所示，

步骤7：确定电源-负荷有功功率平衡约束：

根据负荷预测结果，得到置信水平α(一般可以取α＝90％)的负荷预测条带的下限再结合火电厂C的最小经济技术出力曲线，确定给风电场A和光伏电站B预留的出力曲线

$P_w^{up}\left(t\right)+P_p^{up}\left(t\right)=P_{d,\mathrm{\α}}^{down}\left(t\right)-P_h^{\min }\left(t\right)$

为功率平衡约束下风电的出力上限，为功率平衡约束下光伏的出力上限。

风电场A和光伏电站B的出力分配遵循调度优先级排序原则，即根据预先设定的优先级，选择优先调度风电还是优先调度光伏。例如，考虑发电集团利润时，可以选择优先调度光伏，即根据电源-负荷平衡约束压制风电，让光伏电站尽量满发。

据此，则有

$P_p^{up}\left(t\right)=P_{d,\mathrm{\α}}^{down}\left(t\right)-P_h^{\min }\left(t\right)$

$P_w^{up}\left(t\right)=P_{d,\mathrm{\α}}^{down}\left(t\right)-P_h^{min}\left(t\right)-P_{p,\mathrm{\α}}^{up}\left(t\right)$

其中，为功率平衡约束下光伏的出力上限，为置信水平α(一般可以取α＝90％)的负荷预测条带的下限，为火电厂C的最小经济技术出力。为发电权置换后的光伏出力。

根据典型日负荷曲线计算出由于系统平衡约束给风电场、光伏电站预留的出力如图3所示。

步骤8：生成发电权置换后风电场A和光伏电站B的日前出力下达计划

对于光伏电站B应满足约束

$P_{p,\mathrm{\α}}^{up}\left(t\right)=P_p\left(t\right)+Z_{p,h}\left(t\right)=min\{P_{p,c}\left(t\right),P_p^{up}\left(t\right),P_{p,\mathrm{\α}}^{down}\left(t\right)\}$

其中，P_p,c(t)为调度部门给光伏电站B预留的送出通道容量约束，为风电预测条带的下限，为对光伏电站的电源-负荷功率平衡约束。

对于风电场A应满足约束：

$P_{w,\mathrm{\α}}^{up}\left(t\right)=P_w\left(t\right)+Z_{w,h}\left(t\right)=min\{P_{w,c}\left(t\right),P_w^{up}\left(t\right),P_{w,\mathrm{\α}}^{down}\left(t\right)\}$

其中，P_w,c(t)为调度部门给风电场A预留的送出通道容量约束，为风电预测条带的下限，为对风电场的电源-负荷功率平衡约束。

经发电权置换后的风电出力计划下达曲线如图4所示。

步骤9：是否满足风电机组最小启停时间约束：

$\left\{\begin{array}{c}{T}_{w,i}^{on}\≥T_{w,i}^S\\ T_{w,i}^{off}\≥T_{w,i}^D\end{array}\right.$

其中，为第i台风电机组的连续运行时间及连续停止时间， 为第i台风电机组的最小开停机时间。

$\left\{\begin{array}{c}{T}_{p,i}^{on}\≥T_{p,i}^S\\ T_{p,i}^{off}\≥T_{p,i}^D\end{array}\right.$

其中，为第i台光伏逆变器的连续运行时间及连续停止时间，为第i台光伏逆变器的最小开停机时间。

取风电机组最小启停时间为0，光伏逆变器最小启停时间为0。

步骤10：增加风电机组、光伏逆变器启停限制

将不满足启停约束的机组的启停状态设置为启动并不得停止，转步骤9；

本例中均满足该约束。

步骤11：是否满足火电机组最小启停时间约束：

$\left\{\begin{array}{c}{T}_{h,i}^{on}\≥T_{h,i}^S\\ T_{h,i}^{off}\≥T_{h,i}^D\end{array}\right.,$

其中，为第i台火电机组的连续运行时间及连续停止时间， 为第i台风电机组的最小开停机时间。

取火电机组最小启停时间为5小时。

如果是，进入步骤13；如果否，进入步骤12；

步骤12：增加火电机组启停限制

将不满足启停约束的机组的启停状态设置为启动并不得停止，本算例经检验不需要进行该设置。

步骤13：是否满足风电机组、光伏逆变器爬坡速率约束：

$\begin{array}{c}-\underset{i=1}{\overset{M}{\Σ}}{R}_{w,i}^{down}{U}_{w,i}\left(t\right)\mathrm{\Δ}t\≤\left(P_w\right(t)+Z_{w,h}(t\left)\right)-\\ \left(P_w\right(t-\mathrm{\Δ}t)+Z_{w,h}(t-\mathrm{\Δ}t\left)\right)\≤\underset{i=1}{\overset{M}{\Σ}}{R}_{w,i}^{up}{U}_{w,i}\left(t\right)\mathrm{\Δ}t\end{array},$

$\begin{array}{c}-\underset{i=1}{\overset{K}{\Σ}}{R}_{p,i}^{down}{U}_{p,i}\left(t\right)\mathrm{\Δ}t\≤\left(P_p\right(t)+Z_{p,h}(t\left)\right)-\\ \left(P_p\right(t-\mathrm{\Δ}t)+Z_{p,h}(t-\mathrm{\Δ}t\left)\right)\≤\underset{i=1}{\overset{K}{\Σ}}{R}_{p,i}^{up}{U}_{p,i}\left(t\right)\mathrm{\Δ}t\end{array},$

其中，P_w(t)为风电场在t时刻的原计划上网出力，Z_w,h(t)为风电场通过发电权置换在t时刻得到的功率，为第i台风电机组的爬坡率下限，为第i台风电机组的爬坡率上限，U_w,i(t)为第i台风电机组在t时刻的启停状态，i＝1,2,3,...,M，M为风电场A的风电机组台数。

P_p(t)为光伏电站在t时刻的原计划上网出力，Z_p,h(t)为光伏电站通过发电权置换在t时刻得到的功率，为第i台光伏逆变器的爬坡速率下限，为第i台光伏逆变器的爬坡上限，U_p,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,K，K为光伏电站的光伏逆变器台数。

根据上述约束新能源正常上网出力+置换出力应变化率应介于最大负爬坡速率与最大正爬坡速率之间。此处，风电取20％/min，即40MW/min，光伏50％/min，即25MW/min。

如果是，进入步骤15；如果否，进入步骤14；

步骤14：增加风电机组、光伏电站爬坡率限制

从不满足机组爬坡率的时刻开始，逐点减少发电权置换功率，本算例中经校验不需要进行该项设置。

步骤15：是否满足火电机组爬坡速率约束：

$\begin{array}{c}-\underset{i=1}{\overset{N}{\Σ}}{R}_{h,i}^{down}{U}_{h,i}\left(t\right)\mathrm{\Δ}t\≤\left(P_h\right(t)-Z_{w,h}(t\left)\right)-\\ \left(P_h\right(t-\mathrm{\Δ}t)-Z_{w,h}(t-\mathrm{\Δ}t\left)\right)\≤\underset{i=1}{\overset{N}{\Σ}}{R}_{h,i}^{up}{U}_{h,i}\left(t\right)\mathrm{\Δ}t\end{array}$

其中，P_h(t)为火电厂在t时刻的正常上网出力，Z_w,h(t)为火电厂发电权置换掉的功率，为第i台火电机组的上下爬坡率，U_h,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂B的火电机组台数。

根据上述约束火电机组正常上网出力-置换出力应变化率应介于最大负爬坡速率与最大正爬坡速率之间。此处，取2％/min，即6MW/min。

如果是，进入步骤17；如果否，进入步骤16；

步骤16：增加火电机组爬坡率限制

从不满足机组爬坡率的时刻开始，逐点减少发电权置换功率，并转步骤15；

步骤17：风电出力下达曲线是否经过校验？

如果是，转步骤19。

步骤19：得到经过校验的风电发电权置换后的日出力计划曲线。

步骤20：计算风电场、光伏电站通过日前发电权置换所减少的弃风弃光电量

$Q_{w,h}=\underset{0}{\overset{T}{\∫}}{Z}_{w,h}\left(t\right)dt=\underset{0}{\overset{T}{\∫}}\left(P_{w,\mathrm{\α}}^{up}\right(t)-P_w(t\left)\right)dt$

$Q_{p,h}=\underset{0}{\overset{T}{\∫}}{Z}_{p,h}\left(t\right)dt=\underset{0}{\overset{T}{\∫}}\left(P_{p,\mathrm{\α}}^{up}\right(t)-P_p(t\left)\right)dt$

其中，Z_w,h(t)为风电场通过发电权置换得到的功率，T为24，表示电量统计的时间为明天一天，为发电权置换后风电场A的日前风电出力，P_w(t)为风电原计划出力值。

Z_p,h(t)为光伏电站通过发电权置换得到的功率，T为24，表示电量统计的时间为明天一天，为发电权置换后光伏电站B的日前出力，P_p(t)为光伏原计划出力值。

经计算，该风电场当日减少弃风电量96MWh，光伏电站减少弃光5MWh。

步骤21：计算发电集团增加的利润

$C_w^z\left(Z_{w,h}\right)=\left({\mathrm{\λ}}_{w,b}\right(t)+{\mathrm{\λ}}_z(t\left)\right)\underset{0}{\overset{T}{\∫}}{Z}_{w,h}\left(t\right)dt+\left({\mathrm{\λ}}_{p,b}\right(t)+{\mathrm{\λ}}_z(t\left)\right)\underset{0}{\overset{T}{\∫}}{Z}_{p,h}\left(t\right)dt$

其中，Z_w,h(t)为t时刻风电发电权置换的电量，Z_p,h(t)为t时刻光伏发电权置换的电量，λ_w,b(t)为国家每度风电的新能源补贴，λ_z(t)为火电机组的脱硫标杆上网电价，λ_p,b(t)为国家每度光伏发电的新能源补贴，T＝24表示一天的置换电量所增加的利润。

根据本算例，采用脱硫标杆电价加国家新能源补贴合计，即风电上网电价0.54元/kWh，光伏上网电价0.9元/kWh，经计算，该风电场当天由发电权置换带来的利润为51840元，光伏电站当天由发电权置换带来的利润为4500元，。

步骤22：计算发电权置换社会效益

$S_w^z\left(Z_{w,h}\right)=a\mathrm{\η}\underset{0}{\overset{T}{\∫}}\left(Z_{w,h}\right(t)+Z_{p,h}(t\left)\right)dt+b\mathrm{\γ}\underset{0}{\overset{T}{\∫}}\left(Z_{w,h}\right(t)+Z_{p,h}(t\left)\right)dt$

其中，Z_w,h(t)为t时刻风电发电权置换的电量，Z_p,h(t)为t时刻光伏发电权置换的电量，η为火电厂每度电需要消耗的标准煤，γ为火电厂每度电产生的污染物(CO₂、SO_X、NO_X、烟尘等)。a、b为加权系数。

结合实际情况，设置合理加权系数后，可以算得社会效益的定量结果。

最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。

Claims (10)

1.一种最小化弃风弃光电量的新能源发电权置换电量评估方法，

对于同一发电集团旗下的风电场、光伏电站与火电厂进行日前发电权置换，风电场与火电厂在t时刻的发电权置换功率为Z_w,h(t)，光伏电站与火电厂在t时刻的发电权置换功率为Z_p,h(t)，t为时间，t∈[0,24)；

通过发电权置换最小化弃风弃光电量，即最大化发电权置换电量的目标函数为：

$\begin{array}{cc}\max & Q_{w,h}=\underset{0}{\overset{T}{\∫}}\left({\mathrm{\αZ}}_{w,h}\right(t)+{\mathrm{\βZ}}_{p,h}(t\left)\right)dt\end{array}$

其中，T为24；

其特征在于，该方法具体包括以下步骤：

步骤1：设置新能源出力下达曲线是否经过校验标志位为False，所述新能源包括风电场和光伏电站；

步骤2：确定风电场出力上下限约束以及光伏电站的出力上下限约束：

所述风电场出力上下限约束为：

$\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{P}_{w,i}^{\min }{U}_{w,i}\left(t\right)\≤P_w\left(t\right)+Z_{w,h}\left(t\right)\≤\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{P}_{w,i}^{\max }{U}_{w,i}\left(t\right)$

其中，P_w(t)为风电场在t时刻的原计划上网出力，为第i台风电机组的出力下限，为第i台风电机组的出力上限，U_w,i(t)为第i台风电机组在t时刻的启停状态，i＝1,2,3,...,M，M为风电场的风电机组台数；

所述光伏电站的出力上下限约束为：

$\underset{i=1}{\overset{K}{\Σ}}{P}_{p,i}^{\min }{U}_{p,i}\left(t\right)\≤P_p\left(t\right)+Z_{p,h}\left(t\right)\≤\underset{i=1}{\overset{K}{\Σ}}{P}_{p,i}^{\max }{U}_{p,i}\left(t\right)$

其中，P_p(t)为光伏电站在t时刻的原计划上网出力，为第i台光伏逆变器的出力下限，为第i台光伏逆变器的出力上限，U_p,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,K，K为光伏电站的光伏逆变器台数；

步骤3：确定火电厂出力上下限约束：

$\underset{i=1}{\overset{N}{\Σ}}{P}_{h,i}^{\min }{U}_{h,i}\left(t\right)\≤P_h\left(t\right)-Z_{w,h}\left(t\right)\≤\underset{i=1}{\overset{N}{\Σ}}{P}_{h,i}^{\max }{U}_{h,i}\left(t\right),$

其中，P_h(t)为火电厂在t时刻的原计划上网出力，为第i台火电机组的出力下限，为第i台火电机组的出力上限，U_h,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂的火电机组台数；

步骤4：确定调度部门给风电场预留的送出通道容量约束P_w,c(t)，以及确定调度部门给光伏电站预留的送出通道容量约束P_p,c(t)；

步骤5：确定火电厂的最小经济技术出力约束t∈[0,24)；

步骤6：确定风电场的日前风电功率预测曲线；

步骤7：确定电源-负荷有功功率平衡约束；

步骤8：生成发电权置换后风电场和光伏电站的日前出力下达计划；

步骤9：判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足风电机组及光伏逆变器最小启停时间约束，如果是，进入步骤11，如果否，进入步骤10；

步骤10：增加风电机组及光伏逆变器启停限制，将不满足启停约束的机组的启停状态设置为启动并不得停止，转步骤9；

步骤11：判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足火电机组最小启停时间约束，如果是，进入步骤13，如果否，进入步骤12；

步骤12：增加火电机组启停限制，将不满足启停约束的机组的启停状态设置为启动并不得停止，转步骤11；

步骤13：判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足风电机组及光伏逆变器爬坡速率约束，如果是，进入步骤15，如果否，进入步骤14；

步骤14：增加风电机组、光伏电站爬坡率限制，从不满足机组爬坡率的时刻开始，逐点减少发电权置换功率，并转步骤13；

步骤15：判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足火电机组爬坡速率约束，如果是，进入步骤17，如果否，进入步骤16；

步骤16：增加火电机组爬坡率限制，从不满足机组爬坡率的时刻开始，逐点减少发电权置换功率，并转步骤15；

步骤17：风电出力下达曲线是否经过校验，如果是，转步骤19；如果否，转步骤18；

步骤18：设置风电出力下达曲线是否经过校验标志位为True，并转步骤2；

步骤19：得到经过校验的新能源发电权置换后的日出力计划曲线；

步骤20：计算风电场及光伏电站通过日前发电权置换所减少的弃风弃光电量；

步骤21：计算发电集团增加的利润；

步骤22：计算发电权置换社会效益；

步骤23：完成最小化弃风弃光电量的新能源发电权置换电量评估。

2.根据权利要求1所述的最小化弃风弃光电量的新能源发电权置换电量评估方法，其特征在于，所述步骤4中确定调度部门给风电场预留的送出通道容量约束P_w,c(t)具体为：

根据风电场发电规则，有

P_w(t)+Z_w,h(t)≤P_w,c(t)，

其中，P_w(t)为风电计划出力值，t为时间，t∈[0,24)，Z_w,t(t)为风电场通过发电权置换得到的功率，P_w,c(t)为调度部门给风电场预留的风电送出通道容量；

所述步骤4中确定调度部门给光伏电站预留的送出通道容量约束P_p,c(t)具体为：

根据光伏电站发电规则，则有：

P_p(t)+Z_p,h(t)≤P_p,c(t)，

其中，P_p(t)为光伏电站计划出力值，t为时间，t∈[0,24)，P_p,c(t)为调度部门给光伏电站预留的风电送出通道容量。

3.根据权利要求2所述的最小化弃风弃光电量的新能源发电权置换电量评估方法，其特征在于，所述步骤6中确定风电场的日前风电功率预测曲线具体为：

首先利用风资源监测网络的测风塔的实测风速，修正数值天气预报的预测数据，用修正后的数值天气预报的预测数据结合风电场的日前开机计划通过风电功率预测系统预测明天0-24小时内风电场的出力曲线，取置信水平α的风电预测条带的下限

确定光伏电站的日前光伏功率预测曲线，首先利用光伏资源监测网络的测光站的实测辐照度，修正数值天气预报的预测数据，用修正后的数值天气预报的预测数据结合光伏电站的日前开机计划通过光伏功率预测系统预测明天0-24小时内光伏电站的出力曲线，取置信水平α的光伏预测条带的下限

4.根据权利要求3所述的最小化弃风弃光电量的新能源发电权置换电量评估方法，其特征在于，所述步骤7中确定电源-负荷有功功率平衡约束具体为：

根据负荷预测结果，得到置信水平α的负荷预测条带的下限再结合火电厂的最小经济技术出力曲线，确定给风电场和光伏电站预留的出力曲线

$P_w^{up}\left(t\right)+P_p^{up}\left(t\right)=P_{d,\mathrm{\α}}^{down}\left(t\right)-P_h^{\min }\left(t\right)$

为功率平衡约束下风电的出力上限，为功率平衡约束下光伏的出力上限；

风电场和光伏电站的出力分配遵循调度优先级排序原则，即根据预先设定的优先级，选择优先调度风电还是优先调度光伏；

即根据电源-负荷平衡约束压制风电，让光伏电站尽量满发，

据此，则有：

$P_p^{up}\left(t\right)=P_{d,\mathrm{\α}}^{down}\left(t\right)-P_h^{\min }\left(t\right)$

$P_w^{up}\left(t\right)=P_{d,\mathrm{\α}}^{down}\left(t\right)-P_h^{min}\left(t\right)-P_{p,\mathrm{\α}}^{up}\left(t\right)$

其中，为功率平衡约束下光伏的出力上限，为置信水平α的负荷预测条带的下限，为火电厂的最小经济技术出力。为发电权置换后的光伏出力。

5.根据权利要求4所述的最小化弃风弃光电量的新能源发电权置换电量评估方法，其特征在于，所述步骤8：生成发电权置换后风电场和光伏电站的日前出力下达计划中，光伏电站应满足约束为：

$P_{p,\mathrm{\α}}^{up}\left(t\right)=P_p\left(t\right)+Z_{p,h}\left(t\right)=min\{P_{p,c}\left(t\right),P_p^{up}\left(t\right),P_{p,\mathrm{\α}}^{down}\left(t\right)\},$

其中，P_p,c(t)为调度部门给光伏电站预留的送出通道容量约束，为风电预测条带的下限，为对光伏电站的电源-负荷功率平衡约束；

风电场应满足约束为：

$P_{w,\mathrm{\α}}^{up}\left(t\right)=P_w\left(t\right)+Z_{w,h}\left(t\right)=min\{P_{w,c}\left(t\right),P_w^{up}\left(t\right),P_{w,\mathrm{\α}}^{down}\left(t\right)\},$

其中，P_w,c(t)为调度部门给风电场预留的送出通道容量约束，为风电预测条带的下限，为对风电场的电源-负荷功率平衡约束。

6.根据权利要求5所述的最小化弃风弃光电量的新能源发电权置换电量评估方法，其特征在于，所述步骤9判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足风电机组及光伏逆变器最小启停时间约束中，风电机组最小启停时间约束为：

$\left\{\begin{array}{c}{T}_{w,i}^{on}\≥T_{w,i}^S\\ T_{w,i}^{off}\≥T_{w,i}^D\end{array}\right.,$

其中，为第i台风电机组的连续运行时间及连续停止时间，为第i台风电机组的最小开停机时间；

光伏逆变器最小启停时间约束为：

$\left\{\begin{array}{c}{T}_{p,i}^{on}\≥T_{p,i}^S\\ T_{p,i}^{off}\≥T_{p,i}^D\end{array}\right.$

其中，为第i台光伏逆变器的连续运行时间及连续停止时间，为第i台光伏逆变器的最小开停机时间。

7.根据权利要求6所述的最小化弃风弃光电量的新能源发电权置换电量评估方法，其特征在于，所述步骤11判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足火电机组最小启停时间中，火电机组最小启停时间约束具体为：

$\left\{\begin{array}{c}{T}_{h,i}^{on}\≥T_{h,i}^S\\ T_{h,i}^{off}\≥T_{h,i}^D\end{array}\right.,$

其中，为第i台火电机组的连续运行时间及连续停止时间， 为第i台风电机组的最小开停机时间。

8.根据权利要求7所述的最小化弃风弃光电量的新能源发电权置换电量评估方法，其特征在于，所述步骤13中判断发电权置换后风电场和光伏电站的日前出力下达计划是否满足风电机组及光伏逆变器爬坡速率约束，风电机组爬坡速率约束为：

$\begin{array}{c}-\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{R}_{w,i}^{down}{U}_{w,i}\left(t\right)\mathrm{\Δ}t\≤\left(P_w\left(t\right)+Z_{w,h}\left(t\right)\right)-\\ \left(P_w\left(t-\mathrm{\Δ}t\right)+Z_{w,h}\left(t-\mathrm{\Δ}t\right)\right)\≤\underset{i=1}{\overset{M}{\mathrm{\Σ}}}{R}_{w,i}^{up}{U}_{w,i}\left(t\right)\mathrm{\Δ}t\end{array},$

其中，P_w(t)为风电场在t时刻的原计划上网出力，Z_w,h(t)为风电场通过发电权置换在t时刻得到的功率，为第i台风电机组的爬坡率下限，为第i台风电机组的爬坡率上限，U_w,i(t)为第i台风电机组在t时刻的启停状态，i＝1,2,3,...,M，M为风电场的风电机组台数；

光伏逆变器爬坡速率约束为：

$\begin{array}{c}-\underset{i=1}{\overset{K}{\mathrm{\Σ}}}{R}_{p,i}^{down}{U}_{p,i}\left(t\right)\mathrm{\Δ}t\≤\left(P_p\left(t\right)+Z_{p,h}\left(t\right)\right)-\\ \left(P_p\left(t-\mathrm{\Δ}t\right)+Z_{p,h}\left(t-\mathrm{\Δ}t\right)\right)\≤\underset{i=1}{\overset{K}{\mathrm{\Σ}}}{R}_{p,i}^{up}{U}_{p,i}\left(t\right)\mathrm{\Δ}t\end{array},$

P_p(t)为光伏电站在t时刻的原计划上网出力，Z_p,h(t)为光伏电站通过发电权置换在t时刻得到的功率，为第i台光伏逆变器的爬坡速率下限，为第i台光伏逆变器的爬坡上限，U_p,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,K，K为光伏电站的光伏逆变器台数；

所述步骤15中，火电机组爬坡速率约束为：

$\begin{array}{c}-\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{R}_{h,i}^{down}{U}_{h,i}\left(t\right)\mathrm{\Δ}t\≤\left(P_h\left(t\right)-Z_{w,h}\left(t\right)\right)-\\ \left(P_h\left(t-\mathrm{\Δ}t\right)-Z_{w,h}\left(t-\mathrm{\Δ}t\right)\right)\≤\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{R}_{h,i}^{up}{U}_{h,i}\left(t\right)\mathrm{\Δ}t\end{array}$

其中，P_h(t)为火电厂在t时刻的正常上网出力，Z_w，h(t)为火电厂发电权置换掉的功率，为第i台火电机组的上下爬坡率，U_h，i(t)为第i台火电机组在t时刻的启停状态，i＝1，2，3，...，N，N为火电厂B的火电机组台数。

9.根据权利要求8所述的最小化弃风弃光电量的新能源发电权置换电量评估方法，其特征在于，所述步骤20计算风电场及光伏电站通过日前发电权置换所减少的弃风弃光电量中：

减少的弃风电量为：

$Q_{w,h}=\underset{0}{\overset{T}{\∫}}{Z}_{w,h}\left(t\right)dt=\underset{0}{\overset{T}{\∫}}\left(P_{w,\mathrm{\α}}^{up}\right(t)-P_w(t\left)\right)dt,$

减少的弃光电量为：

$Q_{p,h}=\underset{0}{\overset{T}{\∫}}{Z}_{p,h}\left(t\right)dt=\underset{0}{\overset{T}{\∫}}\left(P_{p,\mathrm{\α}}^{up}\right(t)-P_p(t\left)\right)dt,$

其中，Z_w，h(t)为风电场通过发电权置换得到的功率，T为24，表示电量统计的时间为明天一天，为发电权置换后风电场的日前风电出力，P_w(t)为风电原计划出力值，

Z_p，h(t)为光伏电站通过发电权置换得到的功率，T为24，表示电量统计的时间为明天一天，为发电权置换后光伏电站的日前出力，P_p(t)为光伏原计划出力值。

10.根据权利要求9所述的最小化弃风弃光电量的新能源发电权置换电量评估方法，其特征在于，所述步骤21中计算发电集团增加的利润具体为：

$C_w^z\left(Z_{w,h}\right)=\left({\mathrm{\λ}}_{w,b}\right(t)+{\mathrm{\λ}}_z(t\left)\right)\underset{0}{\overset{T}{\∫}}{Z}_{w,h}\left(t\right)dt+\left({\mathrm{\λ}}_{p,b}\right(t)+{\mathrm{\λ}}_z(t\left)\right)\underset{0}{\overset{T}{\∫}}{Z}_{p,h}\left(t\right)dt$

其中，Z_w，h(t)为t时刻风电发电权置换的电量，Z_p，h(t)为t时刻光伏发电权置换的电量，λ_w，b(t)为国家每度风电的新能源补贴，λ_z(t)为火电机组的脱硫标杆上网电价，λ_p，b(t)为国家每度光伏发电的新能源补贴，T＝24表示一天的置换电量所增加的利润；

所述步骤22中计算发电权置换社会效益具体为：

$S_w^z\left(Z_{w,h}\right)=a\mathrm{\η}\underset{0}{\overset{T}{\∫}}\left(Z_{w,h}\right(t)+Z_{p,h}(t\left)\right)dt+b\mathrm{\γ}\underset{0}{\overset{T}{\∫}}\left(Z_{w,h}\right(t)+Z_{p,h}(t\left)\right)dt$

其中，Z_w，h(t)为t时刻风电发电权置换的电量，Z_p，h(t)为t时刻光伏发电权置换的电量，η为火电厂每度电需要消耗的标准煤，γ为火电厂每度电产生的污染物，a、b为加权系数。