CN105976267A - 最小化弃光电量的光伏电站发电权置换电量评估方法 - Google Patents

Abstract

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

Description

最小化弃光电量的光伏电站发电权置换电量评估方法

技术领域

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

背景技术

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

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

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

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

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

发明内容

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

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

一种最小化弃光电量的光伏电站发电权置换电量评估方法，

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

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

$\begin{array}{cc}\max & Q_{w,h}=\underset{0}{\overset{T}{\∫}}{Z}_{w,h}\left(t\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为光伏电站的光伏逆变器台数；

步骤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时刻的原计划上网出力，Z_w,h(t)为火电厂发电权置换掉的功率，为第i台火电机组的出力下限，为第i台火电机组的出力上限，U_h,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂的火电机组台数；

步骤4：确定调度部门给光伏电站预留的送出通道容量约束P_c(t)，根据光伏电站发电规则，则有：

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

其中，P_w(t)为光伏计划出力值，t为时间，t∈[0,24)，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：完成最小化弃光电量的光伏电站发电权置换电量评估，结束。

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

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

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

优选的，所述步骤7确定光伏电站的日前光伏功率预测曲线具体为：

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

优选的，所述步骤8生成发电权置换后光伏电站的日前光伏出力下达计划具体为：

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

优选的，所述步骤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台光伏逆变器的最小开停机时间。

优选的，所述步骤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时刻的正常上网出力，为第i台光伏逆变器的爬坡率下限，为第i台光伏逆变器的爬坡率上限，U_w,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,M，M为光伏电站的光伏逆变器台数；

所述步骤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时刻的正常上网出力，为第i台火电机组的上下爬坡率，U_h,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂的火电机组台数。

优选的，所述步骤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,$

其中，T为24，表示电量统计的时间为明天一天，为发电权置换后光伏电站的日前光伏出力，P_w(t)为光伏计划出力值。

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

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

其中，λ_b(t)为国家每度电的新能源补贴，λ_z(t)为火电机组的脱硫标杆上网电价，T＝24表示一天的置换电量所增加的利润。

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

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

其中，η为火电厂每度电需要消耗的标准煤，γ为火电厂每度电产生的污染物，a、b为加权系数。

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

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

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

附图说明

图1为本发明实施例所述的最小化弃光电量的光伏电站发电权置换电量评估方法示意图；

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

图3为光伏电站A的日前光伏功率预测曲线示意图；

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

具体实施方式

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

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

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

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

$\begin{array}{cc}\max & Q_{w,h}=\underset{0}{\overset{T}{\∫}}{Z}_{w,h}\left(t\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为光伏电站的光伏逆变器台数；

步骤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时刻的原计划上网出力，Z_w,h(t)为火电厂发电权置换掉的功率，为第i台火电机组的出力下限，为第i台火电机组的出力上限，U_h,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂的火电机组台数；

步骤4：确定调度部门给光伏电站预留的送出通道容量约束P_c(t)，根据光伏电站发电规则，则有：

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

其中，P_w(t)为光伏计划出力值，t为时间，t∈[0,24)，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万千瓦装机的典型光伏电站与一个30万千瓦装机的典型火电厂进行发电权置换，光伏电站有40台500kW光伏逆变器，火电厂有1台300MW火电机组，根据上述计算流程。

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

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

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

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

步骤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_c(t)，根据光伏电站发电规则，则有

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

其中，P_w(t)为光伏计划出力值，t为时间，本发明研究日前置换电量评估，则有t∈[0,24)，Z_w,h(t)为光伏电站通过发电权置换得到的功率，P_c(t)为调度部门给光伏电站A预留的光伏送出通道容量。

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

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

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

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

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

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

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

步骤7：确定光伏电站A的日前光伏功率预测曲线。首先利用测光站的实测辐照强度，结合全天空成像仪，修正数值天气预报(NWP)的预测数据，用修正后的NWP数据结合光伏电站的日前开机计划通过光伏功率预测系统预测明天0-24小时内光伏电站A的出力曲线，取置信水平α(一般可以取α＝90％)的预测条带的下限光伏电站A的日前光伏功率预测曲线如图3所示。

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

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

经发电权置换后的光伏出力计划下达曲线如图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台光伏逆变器的最小开停机时间。

取光伏逆变器最小启停时间为0。

步骤10：增加光伏逆变器启停限制

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

步骤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}$

其中，P_w(t)为光伏电站在t时刻的正常上网出力，Z_w,h(t)为光伏电站通过发电权置换在t时刻得到的功率，为第i台光伏逆变器的爬坡率下限，为第i台光伏逆变器的爬坡率上限，U_w,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,M，M为光伏电站A的光伏逆变器台数。

根据上述约束光伏逆变器正常上网出力+置换出力应变化率应介于最大负爬坡速率与最大正爬坡速率之间。此处，取40％/min，即80MW/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：增加火电机组爬坡率限制

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

步骤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$

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

经计算，该光伏电站当日减少弃光电量40MWh。

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

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

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

根据本算例，脱硫标杆电价加国家新能源补贴合计为1.0元/kWh，经计算，该发电集团当天由发电权置换带来的利润为40000元。

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

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

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

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

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

Claims (10)

1.一种最小化弃光电量的光伏电站发电权置换电量评估方法，

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

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

$\begin{array}{cc}\max & Q_{w,h}=\underset{0}{\overset{T}{\∫}}{Z}_{w,h}\left(t\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为光伏电站的光伏逆变器台数；

步骤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时刻的原计划上网出力，Z_w,h(t)为火电厂发电权置换掉的功率，为第i台火电机组的出力下限，为第i台火电机组的出力上限，U_h,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂的火电机组台数；

步骤4：确定调度部门给光伏电站预留的送出通道容量约束P_c(t)，根据光伏电站发电规则，则有：

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

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

步骤5：确定火电厂的最小经济技术出力约束

步骤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所述的最小化弃光电量的光伏电站发电权置换电量评估方法，其特征在于，所述步骤6中确定电源-负荷有功功率平衡约束具体为：

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

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

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

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

4.根据权利要求3所述的最小化弃光电量的光伏电站发电权置换电量评估方法，其特征在于，所述步骤8生成发电权置换后光伏电站的日前光伏出力下达计划具体为：

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

5.根据权利要求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台光伏逆变器的最小开停机时间。

6.根据权利要求5所述的最小化弃光电量的光伏电站发电权置换电量评估方法，其特征在于，所述步骤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台光伏逆变器的最小开停机时间。

7.根据权利要求6所述的最小化弃光电量的光伏电站发电权置换电量评估方法，其特征在于，所述步骤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时刻的正常上网出力，为第i台光伏逆变器的爬坡率下限，为第i台光伏逆变器的爬坡率上限，U_w,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,M，M为光伏电站的光伏逆变器台数；

所述步骤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时刻的正常上网出力，为第i台火电机组的上下爬坡率，U_h,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂的火电机组台数。

8.根据权利要求7所述的最小化弃光电量的光伏电站发电权置换电量评估方法，其特征在于，所述步骤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,$

其中，T为24，表示电量统计的时间为明天一天，为发电权置换后光伏电站的日前光伏出力，P_w(t)为光伏计划出力值。

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

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

其中，λ_b(t)为国家每度电的新能源补贴，λ_z(t)为火电机组的脱硫标杆上网电价，T＝24表示一天的置换电量所增加的利润。

10.根据权利要求9所述的最小化弃光电量的光伏电站发电权置换电量评估方法，其特征在于，所述步骤22计算发电权置换社会效益具体为：

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

其中，η为火电厂每度电需要消耗的标准煤，γ为火电厂每度电产生的污染物，a、b为加权系数。