CN106026191B - 用于最小化弃光电量光伏电站发电权置换电量评估的方法 - Google Patents
用于最小化弃光电量光伏电站发电权置换电量评估的方法 Download PDFInfo
- Publication number
- CN106026191B CN106026191B CN201610338888.9A CN201610338888A CN106026191B CN 106026191 B CN106026191 B CN 106026191B CN 201610338888 A CN201610338888 A CN 201610338888A CN 106026191 B CN106026191 B CN 106026191B
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- msubsup
- photovoltaic
- power
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000010248 power generation Methods 0.000 title claims abstract description 151
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000005611 electricity Effects 0.000 title claims abstract description 21
- 238000006073 displacement reaction Methods 0.000 title abstract 3
- 230000003287 optical effect Effects 0.000 title abstract 2
- 230000009194 climbing Effects 0.000 claims description 41
- 230000008901 benefit Effects 0.000 claims description 11
- 238000011156 evaluation Methods 0.000 claims description 10
- 238000012795 verification Methods 0.000 claims description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 8
- 239000003245 coal Substances 0.000 claims description 4
- 238000006477 desulfuration reaction Methods 0.000 claims description 4
- 230000023556 desulfurization Effects 0.000 claims description 4
- 238000007599 discharging Methods 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 239000003344 environmental pollutant Substances 0.000 claims description 3
- 239000003550 marker Substances 0.000 claims description 3
- 231100000719 pollutant Toxicity 0.000 claims description 3
- 230000008569 process Effects 0.000 abstract description 7
- 238000005096 rolling process Methods 0.000 abstract description 3
- 238000006467 substitution reaction Methods 0.000 description 6
- 238000011160 research Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000002860 competitive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
- H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P80/00—Climate change mitigation technologies for sector-wide applications
- Y02P80/20—Climate change mitigation technologies for sector-wide applications using renewable energy
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Photovoltaic Devices (AREA)
- Engine Equipment That Uses Special Cycles (AREA)
Abstract
本发明公开了一种用于最小化弃光电量光伏电站发电权置换电量评估的方法，通过对光伏电站和火电厂、水电站日前发电权置换进行评估。从而最大化新能源发电量，最小化弃光。并简化交易电量滚动修正的过程，只进行日前的置换电量评估，简化计算过程。从而达到提高新能源利用率以及提高发电企业利润的目的。
Description
技术领域
本发明涉及大规模新能源发电运行评估技术，具体地，涉及一种用于最小化弃光电量光伏电站发电权置换电量评估的方法。
背景技术
发电权是发电厂商发电运行并将电能出售给电网或通过电网出售给终端用户的权利。在管制电力系统中，它对应于发电厂商在特定机制下获得的发电配额；而在自由竞争的电力市场中，它是指发电厂商在远期电力市场或当前电力市场中通过报价竞标获得的电量。
发电权置换是指上述配额或中标电量的交换。发电权置换本质上是电力系统在调度、发电资源重新配置的过程。积极合理的发电权置换能够使经济高效的发电形式替代高能耗、高污染的发电形式承担负荷，提升电力系统运行的经济性和可靠性、保护环境、促进可持续发展，同时使参与发电权置换的发电集团获益。发电权置换交易完成后，各发电厂的出力计划会相应调整，这也意味着电力系统运行方式的改变。
发电权置换的研究背景是一方面，当前我国部分地区尤其是“三北”风电、光伏等新能源装机增长较快，用电负荷增长较慢，造成大量弃风、弃光；另一方面，由于火电厂消耗大量煤炭等非可再生资源，同时造成环境污染严重，雾霾天气频发。发电权置换的主要对象是常规电源和新能源企业，通过将常规电源的发电量置换成新能源电源的发电量，从而实现降低发电成本、提高企业利润，同时节能减排、保护环境。发电权置换的形式包括中长期发电权置换、短期发电权置换、以及日前发电权置换等置换方式，需要指出的是日前发电权置换需要考虑线路传输容量等运行问题。
目前，发电权置换在发电权重新配置，特别是日前电力系统运行方式调整等方面的效益并未得到充分显示。已开展的相关研究也相对较少。
专利201510375285.1提出了一种“一种促进风电消纳的多时间尺度发电权交易方法”，研究的是不同发电厂商之间的发电权交易，涉及电量交易和资金交易，研究目标是最大化社会效益总和，并通过逐小时修正发电权交易计划最终确定发电权交易电量，因为风电预测的不确定性、负荷预测的不确定性以及机组开停机的变化性等因素，该过程本身就是在一定概率意义下的评估，进行逐小时滚动修正会增加大量的工作量，其效果并不显著。
发明内容
本发明的目的在于，针对上述问题，提出一种用于最小化弃光电量光伏电站发电权置换电量评估的方法，以实现提高新能源利用率以及提高发电企业利润的优点。
为实现上述目的，本发明采用的技术方案是：
一种用于最小化弃光电量光伏电站发电权置换电量评估的方法，
对于同一发电集团旗下的光伏电站与火电厂、水电站进行日前发电权置换，光伏电站与火电厂在t时刻的发电权置换功率为Zw,h(t)，光伏电站与水电站在t时刻的发电权置换功率为Zw,a(t)，t为时间，t∈[0,24)；
通过发电权置换最小化弃光电量，即最大化发电权置换电量的目标函数为：
其中，T为24；
该方法具体包括以下步骤：
步骤1：设置光伏出力下达曲线是否经过校验标志位为False；
步骤2：确定光伏电站出力上下限约束；
其中，Pw(t)为光伏电站在t时刻的原计划上网出力，为第i台光伏逆变器的出力下限，为第i台光伏逆变器的出力上限，Ua，i(t)Uw,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,M，M为光伏电站A的光伏逆变器台数；
步骤3：确定火电厂及水电站出力上下限约束，火电厂出力上下限约束为：
其中，Ph(t)为火电厂在t时刻的原计划上网出力，为第i台火电机组的出力下限，为第i台火电机组的出力上限，Uh,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂B的火电机组台数；
水电站出力上下限约束为：
其中，Pa(t)为水电厂A在t时刻的原计划上网出力，为第i台水电机组的出力下限，为第i台水电机组的出力上限，Uw,i(t)为第i台水电机组在在t时刻的启停状态，i＝1,2,3,...,K，K为水电厂水电机组的台数；
步骤4：确定调度部门给光伏电站预留的送出通道容量约束Pc(t)，根据光伏电站发电规则，有：
Pw(t)+Zw,h(t)+Zw,a(t)≤Pw,c(t)，
其中，Pw(t)为光伏计划出力值，t为时间，t∈[0,24)，Pc(t)为调度部门给光伏电站预留的光伏送出通道容量；
步骤5：确定火电厂B的最小经济技术出力约束t∈[0,24)，及确定水电厂的不弃水最小出力约束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中确定电源-负荷有功功率平衡约束具体为：
根据负荷预测结果，得到置信水平α的预测条带的下限再结合火电厂的最小经济技术出力曲线和水电厂的不弃水最小出力约束，确定给光伏电站预留的出力曲线：
其中，为置信水平α的预测条带的下限，为火电厂最小经济技术出力曲线，为水电厂的不弃水最小出力曲线。
优选的，所述步骤7中确定光伏电站的日前光伏功率预测曲线具体为：
首先利用测光站的实测辐照强度，结合全天空成像仪，修正数值天气预报的预测数据，用修正后的数值天气预报的预测数据结合光伏电站的日前开机计划通过光伏功率预测系统预测明天0-24小时内光伏电站的出力曲线，取置信水平α的预测条带的下限
优选的，所述步骤8生成发电权置换后光伏电站的日前光伏出力下达计划具体为：
优选的，所述步骤9中，光伏组件最小启停时间约束为：
其中，为第i台光伏逆变器的连续运行时间及连续停止时间，为第i台光伏逆变器的最小开停机时间。
优选的，所述步骤11中，火电机组最小启停时间约束为：
其中，为第i台火电机组的连续运行时间及连续停止时间， 为第i台风电机组的最小开停机时间；
所述步骤11中，水电机组最小启停时间约束为：
其中，为第i台水电机组连续运行时间及连续停止时间， 为第i台水电机组的最小开停机时间。
优选的，所述步骤13中，光伏逆变器爬坡速率约束为：
其中，Pw(t)为光伏电站在t时刻的正常上网出力，为第i台光伏逆变器的爬坡率下限，为第i台光伏逆变器的爬坡率上限，Uw,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,M，M为光伏电站的光伏逆变器台数；
所述步骤15中，火电机组爬坡速率约束为：
其中，Ph(t)为火电厂在t时刻的原计划上网出力，为第i台火电机组的上下爬坡率，Uh,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂的火电机组台数；
所述步骤15中，水电机组爬坡速率约束为：
其中，Pa(t)为水电厂t时刻的原计划上网出力，为第i台水电机组的上下爬坡率，Ua,i(t)为第i台水电机组在t时刻的启停状态，i＝1,2,3,...,K，K为水电机组台数。
优选的，所述步骤20计算光伏电站通过日前发电权置换所减少的弃光电量具体为：
其中，Zw,h(t)和Za,h(t)为光伏电站通过发电权置换得到的功率，T为24，表示电量统计的时间为明天一天，为发电权置换后光伏电站的日前光伏出力，Pw(t)为光伏计划出力值。
优选的，所述步骤21计算发电集团增加的利润具体为：
其中，Zw,h(t)和Za,h(t)为t时刻发电权置换的电量，λb(t)为国家每度电的新能源补贴，λz(t)为火电机组的脱硫标杆上网电价，T＝24表示一天的置换电量所增加的利润。
优选的，所述步骤22计算发电权置换社会效益具体为：
其中，Zw,h(t)为t时刻发电权置换的电量，η为火电厂每度电需要消耗的标准煤，γ为火电厂每度电产生的污染物，a、b为加权系数。
本发明的技术方案具有以下有益效果：
本发明的技术方案同一发电集集团旗下不同发电形式的发电厂之间(如水电站、火电厂和光伏电站)的发电权置换，只涉及电量置换，不涉及资金交易，为了最大化新能源发电量，最小化弃光从而简化交易电量滚动修正的过程，只进行日前的置换电量评估，优化简化计算过程，达到提高新能源利用率以及提高发电企业利润的目的。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明实施例所述的用于最小化弃光电量光伏电站发电权置换电量评估的方法的流程图；
图2为典型日负荷曲线计算出由于系统平衡约束给光伏电站预留的出力曲线示意图；
图3为光伏电站的日前光伏功率预测曲线示意图；
图4为经发电权置换后的光伏出力计划下达曲线示意图。
结合附图，本发明实施例中附图标记如下：
1-负荷预测条带下限；2-火电最小经济技术出力；3-给光伏电站的预留出力；4-水电最小经济技术出力。
具体实施方式
以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
如图1所示，一种用于最小化弃光电量光伏电站发电权置换电量评估的方法，
对于同一发电集团旗下的光伏电站与火电厂、水电站进行日前发电权置换，光伏电站与火电厂在t时刻的发电权置换功率为Zw,h(t)，光伏电站与水电站在t时刻的发电权置换功率为Zw,a(t)，t为时间，t∈[0,24)；
通过发电权置换最小化弃光电量，即最大化发电权置换电量的目标函数为：
其中，T为24；
该方法具体包括以下步骤：
步骤1：设置光伏出力下达曲线是否经过校验标志位为False；
步骤2：确定光伏电站出力上下限约束；
其中，Pw(t)为光伏电站在t时刻的原计划上网出力，为第i台光伏逆变器的出力下限，为第i台光伏逆变器的出力上限，Uw,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,M，M为光伏电站A的光伏逆变器台数；
步骤3：确定火电厂及水电站出力上下限约束，火电厂出力上下限约束为：
其中，Ph(t)为火电厂在t时刻的原计划上网出力，为第i台火电机组的出力下限，为第i台火电机组的出力上限，Uh,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂B的火电机组台数；
水电站出力上下限约束为：
其中，Pa(t)为水电厂A在t时刻的原计划上网出力，为第i台水电机组的出力下限，为第i台水电机组的出力上限，Ua，i(t)为第i台水电机组在在t时刻的启停状态，i＝1,2,3,...,K，K为水电厂水电机组的台数；
步骤4：确定调度部门给光伏电站预留的送出通道容量约束Pc(t)，根据光伏电站发电规则，有：
Pw(t)+Zw,t(t)+Zw,a(t)≤Pw,c(t)，
其中，Pw(t)为光伏计划出力值，t为时间，t∈[0,24)，Pc(t)为调度部门给光伏电站预留的光伏送出通道容量；
步骤5：确定火电厂B的最小经济技术出力约束t∈[0,24)，及确定水电厂的不弃水最小出力约束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万千瓦装机的典型火电厂、5万水电站进行发电权置换，光伏电站有40台500kW光伏逆变器，火电厂有1台300MW火电机组，水电站有2台25MW机组，根据上述计算流程。
目标函数：通过发电权置换最小化弃光电量，即最大化发电权置换电量
其中，Zw,h(t)为光伏电站与火电厂发电权置换得到的功率，Zw,a(t)为光伏电站与水电站发电权置换得到的功率，T为24，表示电量统计的时间为明天一天。
步骤1：设置光伏出力下达曲线是否经过校验标志位为False。
步骤2：确定光伏电站出力上下限约束：
其中，Pw(t)为光伏电站在t时刻的原计划上网出力，Zw,h(t)为光伏电站与火电厂发电权置换得到的功率，Zw,a(t)为光伏电站与水电站发电权置换得到的功率，为第i台光伏逆变器的出力下限，为第i台光伏逆变器的出力上限，Uw,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,M，M为光伏电站的光伏逆变器台数。原计划上网出力根据调度部门下发的日前计划结合天前光伏功率预测结果进行确定，启停状态根据光伏电站日前开机计划获得。
根据上述约束，光伏电站正常上网出力+置换出力应大于0，小于所有开机运行的光伏逆变器总容量，根据当前运行情况，明日光伏电站有2台逆变器检修停运，因此光伏电站正常上网出力+置换出力应小于38*500kW＝190MW。
步骤3：确定火电厂及水电站出力上下限约束：
其中，Ph(t)为火电厂在t时刻的正常上网出力，Zw,h(t)为火电厂发电权置换掉的功率，为第i台火电机组的出力下限，为第i台火电机组的出力上限，Uh,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂的火电机组台数。火电厂原计划上网出力由调度部门下发的日前计划确定。
其中，Pa(t)为水电厂在t时刻的原计划上网出力，Zw,a(t)为水电厂置换掉的功率，为第i台水电机组的出力下限，为第i台水电机组的出力上限，Uw,i(t)为第i台水电机组在在t时刻的启停状态，i＝1,2,3,...,K，K为水电厂水电机组的台数。水电厂原计划上网出力由调度部门下发的日前计划确定。
根据上述约束，火电机组正常上网出力减去置换出力应大于0，小于300MW。水电机组正常上网出力减去置换出力应大于0，小于50MW。
步骤4：确定调度部门给光伏电站预留的送出通道容量约束Pc(t)，根据光伏电站发电规则，则有
Pw(t)+Zw,h(t)+Zw,a(t)≤Pw,c(t)
其中，Pw(t)为光伏计划出力值，t为时间，本发明研究日前置换电量评估，则有t∈[0,24)，Zw,h(t)为光伏电站与火电厂发电权置换得到的功率，Zw,a(t)为光伏电站与水电厂发电权置换得到的功率，Pc(t)为调度部门给光伏电站预留的光伏送出通道容量。
调度部门为光伏电站明日预留的送出通道容量为180MW。光伏电站正常上网出力+置换出力应小于180MW。
步骤5：确定火电厂的最小经济技术出力约束t∈[0,24)。
确定水电厂的不弃水最小出力约束t∈[0,24)。
目前，新型火电机组最小技术出力已经可以达到50％额定功率，因此，本算例中为150MW。水电厂不弃水最小出力10MW。
步骤6：确定电源-负荷有功功率平衡约束：
根据负荷预测结果，得到置信水平α(一般可以取α＝90％)的负荷预测条带的下限再结合火电厂的最小经济技术出力曲线和水电厂的不弃水最小出力约束，确定给光伏电站预留的出力曲线
其中，为置信水平α(一般可以取α＝90％)的预测条带的下限，为火电厂最小经济技术出力曲线，为水电厂的不弃水最小出力曲线。
根据典型日负荷曲线计算出由于系统平衡约束给光伏电站预留的出力如图2所示。
步骤7：确定光伏电站的日前光伏功率预测曲线。首先利用测光站的实测辐照强度，结合全天空成像仪，修正数值天气预报(NWP)的预测数据，用修正后的NWP数据结合光伏电站的日前开机计划通过光伏功率预测系统预测明天0-24小时内光伏电站的出力曲线，取置信水平α(一般可以取α＝90％)的预测条带的下限
光伏电站的日前光伏功率预测曲线如图3所示。
步骤8：生成发电权置换后光伏电站的日前光伏出力下达计划
经发电权置换后的光伏出力计划下达曲线如图4所示。
步骤9：是否满足光伏逆变器最小启停时间约束：
其中，为第i台光伏逆变器的连续运行时间及连续停止时间，为第i台光伏逆变器的最小开停机时间。
取光伏逆变器最小启停时间为0。
步骤10：增加光伏逆变器启停限制
将不满足启停约束的机组的启停状态设置为启动并不得停止，本算例经检验不需要进行该设置。
步骤11：是否满足火电机组最小启停时间约束：
其中，为第i台火电机组的连续运行时间及连续停止时间， 为第i台风电机组的最小开停机时间。
其中，为第i台水电机组连续运行时间及连续停止时间， 为第i台水电机组的最小开停机时间。
取火电机组最小启停时间为5小时，水电机组最小开停机时间为0。
如果是，进入步骤13；如果否，进入步骤12；
步骤12：增加火电机组及水电机组启停限制
将不满足启停约束的机组的启停状态设置为启动并不得停止，本算例经检验不需要进行该设置。
步骤13：是否满足光伏逆变器爬坡速率约束：
其中，Pw(t)为光伏电站在t时刻的正常上网出力，Zw,t(t)为光伏电站与火电厂发电权置换在t时刻得到的功率，Zw,h(t)为光伏电站与水电厂发电权置换在t时刻得到的功率，为第i台光伏逆变器的爬坡率下限，为第i台光伏逆变器的爬坡率上限，Uw,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,M，M为光伏电站的光伏逆变器台数。
根据上述约束光伏逆变器正常上网出力+置换出力应变化率应介于最大负爬坡速率与最大正爬坡速率之间。此处，取40％/min，即80MW/min。
如果是，进入步骤15；如果否，进入步骤14；
步骤14：增加光伏逆变器爬坡率限制
从不满足机组爬坡率的时刻开始，逐点减少发电权置换功率，本算例中经校验不需要进行该项设置。
步骤15：是否满足火电机组及水电机组爬坡速率约束：
其中，Ph(t)为火电厂在t时刻的原计划上网出力，为第i台火电机组的上下爬坡率，Uh,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂的火电机组台数。
其中，Pa(t)为水电厂t时刻的原计划上网出力，Za,h(t)为水电厂发电权置换掉的功率，为第i台水电机组的上下爬坡率，Ua,i(t)为第i台水电机组在t时刻的启停状态，i＝1,2,3,...,K，K为水电机组台数。
根据上述约束火电机组正常上网出力减去置换出力应变化率应介于最大负爬坡速率与最大正爬坡速率之间。此处，火电机组取2％/min，即6MW/min，水电机组取100％/min，即50MW/min。
如果是，进入步骤17；如果否，进入步骤16；
步骤16：增加火电机组、水电机组爬坡率限制
从不满足机组爬坡率的时刻开始，逐点减少发电权置换功率，本算例中经校验不需要进行该项设置。
步骤17：光伏出力下达曲线是否经过校验？
如果是，转步骤19。
步骤19：得到经过校验的光伏发电权置换后的日出力计划曲线。
步骤20：计算光伏电站通过日前发电权置换所减少的弃光电量
其中，Zw,h(t)和Za,h(t)为光伏电站通过发电权置换得到的功率，T为24，表示电量统计的时间为明天一天，为发电权置换后光伏电站的日前光伏出力，Pw(t)为光伏计划出力值。
经计算，该光伏电站当日减少弃光电量40MWh。
步骤21：计算发电集团增加的利润
其中，Zw,h(t)和Za,h(t)为t时刻发电权置换的电量，λb(t)为国家每度电的新能源补贴，λz(t)为火电机组的脱硫标杆上网电价，T＝24表示一天的置换电量所增加的利润。
根据本算例，脱硫标杆电价加国家新能源补贴合计为1.0元/kWh，经计算，该发电集团当天由发电权置换带来的利润为40000元。
步骤22：计算发电权置换社会效益
其中，Zw,h(t)为t时刻发电权置换的电量，η为火电厂每度电需要消耗的标准煤，γ为火电厂每度电产生的污染物(CO2、SOX、NOX、烟尘等)。a、b为加权系数。
结合实际情况，设置合理加权系数后，可以算得社会效益的定量结果。
最后应说明的是：以上所述仅为本发明的优选实施例而已，并不用于限制本发明，尽管参照前述实施例对本发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
Claims (10)
1.一种用于最小化弃光电量光伏电站发电权置换电量评估的方法，
对于同一发电集团旗下的光伏电站与火电厂、水电站进行日前发电权置换，光伏电站与火电厂在t时刻的发电权置换功率为Zw,h(t)，光伏电站与水电站在t时刻的发电权置换功率为Zw,a(t)，t为时间，t∈[0,24)；
通过发电权置换最小化弃光电量，即最大化发电权置换电量的目标函数为：
其中，T为24；
其特征在于，该方法具体包括以下步骤：
步骤1：设置光伏出力下达曲线是否经过校验标志位为False；
步骤2：确定光伏电站出力上下限约束；
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>min</mi>
</msubsup>
<msub>
<mi>U</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>max</mi>
</msubsup>
<msub>
<mi>U</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中，Pw(t)为光伏电站在t时刻的原计划上网出力，为第i台光伏逆变器的出力下限，为第i台光伏逆变器的出力上限，Uw,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,M，M为光伏电站A的光伏逆变器台数；
步骤3：确定火电厂及水电站出力上下限约束，火电厂出力上下限约束为：
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>min</mi>
</msubsup>
<msub>
<mi>U</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mi>h</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>max</mi>
</msubsup>
<msub>
<mi>U</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中，Ph(t)为火电厂在t时刻的原计划上网出力，为第i台火电机组的出力下限，为第i台火电机组的出力上限，Uh,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂B的火电机组台数；
水电站出力上下限约束为：
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>min</mi>
</msubsup>
<msub>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<msub>
<mi>P</mi>
<mi>a</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msubsup>
<mi>P</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>max</mi>
</msubsup>
<msub>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中，Pa(t)为水电厂A在t时刻的原计划上网出力，为第i台水电机组的出力下限，为第i台水电机组的出力上限，Ua，i(t)为第i台水电机组在在t时刻的启停状态，i＝1,2,3,...,K，K为水电厂水电机组的台数；
步骤4：确定调度部门给光伏电站预留的送出通道容量约束Pc(t)，根据光伏电站发电规则，有：
Pw(t)+Zw,h(t)+Zw,a(t)≤Pw,c(t)，
其中，Pw(t)为光伏计划出力值，t为时间，t∈[0,24)，Pc(t)为调度部门给光伏电站预留的光伏送出通道容量；
步骤5：确定火电厂B的最小经济技术出力约束t∈[0,24)，及确定水电厂的不弃水最小出力约束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所述的用于最小化弃光电量光伏电站发电权置换电量评估的方法，其特征在于，所述步骤6中确定电源-负荷有功功率平衡约束具体为：
根据负荷预测结果，得到置信水平α的预测条带的下限再结合火电厂的最小经济技术出力曲线和水电厂的不弃水最小出力约束，确定给光伏电站预留的出力曲线：
<mrow>
<msubsup>
<mi>P</mi>
<mi>p</mi>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>d</mi>
<mo>,</mo>
<mi>&alpha;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mi>h</mi>
<mi>min</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msubsup>
<mi>P</mi>
<mi>a</mi>
<mi>min</mi>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
</mrow>
其中，为置信水平α的预测条带的下限，为火电厂最小经济技术出力曲线，为水电厂的不弃水最小出力曲线。
3.根据权利要求2所述的用于最小化弃光电量光伏电站发电权置换电量评估的方法，其特征在于，所述步骤7中确定光伏电站的日前光伏功率预测曲线具体为：
首先利用测光站的实测辐照强度，结合全天空成像仪，修正数值天气预报的预测数据，用修正后的数值天气预报的预测数据结合光伏电站的日前开机计划通过光伏功率预测系统预测明天0-24小时内光伏电站的出力曲线，取置信水平α的预测条带的下限
4.根据权利要求3所述的用于最小化弃光电量光伏电站发电权置换电量评估的方法，其特征在于，所述步骤8生成发电权置换后光伏电站的日前光伏出力下达计划具体为：
<mrow>
<msubsup>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>&alpha;</mi>
</mrow>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>P</mi>
<mi>w</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mo>{</mo>
<msub>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mi>w</mi>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>,</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>&alpha;</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mo>}</mo>
<mo>.</mo>
</mrow>
5.根据权利要求4所述的用于最小化弃光电量光伏电站发电权置换电量评估的方法，其特征在于，所述步骤9中，光伏组件最小启停时间约束为：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>S</mi>
</msubsup>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>o</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>D</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中，为第i台光伏逆变器的连续运行时间及连续停止时间，为第i台光伏逆变器的最小开停机时间。
6.根据权利要求5所述的用于最小化弃光电量光伏电站发电权置换电量评估的方法，其特征在于，所述步骤11中，火电机组最小启停时间约束为：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>T</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>S</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>T</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>o</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>D</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中，为第i台火电机组的连续运行时间及连续停止时间， 为第i台风电机组的最小开停机时间；
所述步骤11中，水电机组最小启停时间约束为：
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>o</mi>
<mi>n</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>S</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msubsup>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>o</mi>
<mi>f</mi>
<mi>f</mi>
</mrow>
</msubsup>
<mo>&GreaterEqual;</mo>
<msubsup>
<mi>T</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mi>D</mi>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
</mrow>
其中，为第i台水电机组连续运行时间及连续停止时间， 为第i台水电机组的最小开停机时间。
7.根据权利要求6所述的用于最小化弃光电量光伏电站发电权置换电量评估的方法，其特征在于，所述步骤13中，光伏逆变器爬坡速率约束为：
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msubsup>
<mi>R</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msubsup>
<msub>
<mi>U</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>&le;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>w</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>w</mi>
</msub>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>a</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
<mo>&le;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>M</mi>
</munderover>
<msubsup>
<mi>R</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msubsup>
<msub>
<mi>U</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>,</mo>
</mrow>
其中，Pw(t)为光伏电站在t时刻的正常上网出力，为第i台光伏逆变器的爬坡率下限，为第i台光伏逆变器的爬坡率上限，Uw,i(t)为第i台光伏逆变器在t时刻的启停状态，i＝1,2,3,...,M，M为光伏电站的光伏逆变器台数；
所述步骤15中，火电机组爬坡速率约束为：
<mrow>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>R</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msubsup>
<msub>
<mi>U</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>&le;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>h</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>h</mi>
</msub>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<msubsup>
<mi>R</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msubsup>
<msub>
<mi>U</mi>
<mrow>
<mi>h</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>,</mo>
</mrow>
其中，Ph(t)为火电厂在t时刻的原计划上网出力，为第i台火电机组的上下爬坡率，Uh,i(t)为第i台火电机组在t时刻的启停状态，i＝1,2,3,...,N，N为火电厂的火电机组台数；
所述步骤15中，水电机组爬坡速率约束为：
<mrow>
<mo>-</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msubsup>
<mi>R</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>d</mi>
<mi>o</mi>
<mi>w</mi>
<mi>n</mi>
</mrow>
</msubsup>
<msub>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>&le;</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>a</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>P</mi>
<mi>a</mi>
</msub>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mi>t</mi>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>t</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>&le;</mo>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>K</mi>
</munderover>
<msubsup>
<mi>R</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msubsup>
<msub>
<mi>U</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>,</mo>
</mrow>
其中，Pa(t)为水电厂t时刻的原计划上网出力，为第i台水电机组的上下爬坡率，Ua,i(t)为第i台水电机组在t时刻的启停状态，i＝1,2,3,...,K，K为水电机组台数。
8.根据权利要求7所述的用于最小化弃光电量光伏电站发电权置换电量评估的方法，其特征在于，所述步骤20计算光伏电站通过日前发电权置换所减少的弃光电量具体为：
<mrow>
<msub>
<mi>Q</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>=</mo>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>T</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>=</mo>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>T</mi>
</munderover>
<mrow>
<mo>(</mo>
<msubsup>
<mi>P</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>&alpha;</mi>
</mrow>
<mrow>
<mi>u</mi>
<mi>p</mi>
</mrow>
</msubsup>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>-</mo>
<msub>
<mi>P</mi>
<mi>w</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>,</mo>
</mrow>
其中，Zw,h(t)和Za,h(t)为光伏电站通过发电权置换得到的功率，T为24，表示电量统计的时间为明天一天，为发电权置换后光伏电站的日前光伏出力，Pw(t)为光伏计划出力值。
9.根据权利要求8所述的用于最小化弃光电量光伏电站发电权置换电量评估的方法，其特征在于，所述步骤21计算发电集团增加的利润具体为：
<mrow>
<msubsup>
<mi>C</mi>
<mi>w</mi>
<mi>z</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>&lambda;</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>b</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&lambda;</mi>
<mi>z</mi>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>T</mi>
</munderover>
<mrow>
<mo>(</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>a</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>,</mo>
</mrow>
其中，Zw,h(t)和Za,h(t)为t时刻发电权置换的电量，λb(t)为国家每度电的新能源补贴，λz(t)为火电机组的脱硫标杆上网电价，T＝24表示一天的置换电量所增加的利润。
10.根据权利要求9所述的用于最小化弃光电量光伏电站发电权置换电量评估的方法，其特征在于，所述步骤22计算发电权置换社会效益具体为：
<mrow>
<msubsup>
<mi>S</mi>
<mi>w</mi>
<mi>z</mi>
</msubsup>
<mrow>
<mo>(</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>a</mi>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>T</mi>
</munderover>
<msub>
<mi>&eta;Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
<mo>+</mo>
<mi>b</mi>
<munderover>
<mo>&Integral;</mo>
<mn>0</mn>
<mi>T</mi>
</munderover>
<msub>
<mi>&gamma;Z</mi>
<mrow>
<mi>w</mi>
<mo>,</mo>
<mi>h</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mi>d</mi>
<mi>t</mi>
</mrow>
其中，Zw,h(t)为t时刻发电权置换的电量，η为火电厂每度电需要消耗的标准煤，γ为火电厂每度电产生的污染物，a、b为加权系数。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610338888.9A CN106026191B (zh) | 2016-05-20 | 2016-05-20 | 用于最小化弃光电量光伏电站发电权置换电量评估的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610338888.9A CN106026191B (zh) | 2016-05-20 | 2016-05-20 | 用于最小化弃光电量光伏电站发电权置换电量评估的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106026191A CN106026191A (zh) | 2016-10-12 |
| CN106026191B true CN106026191B (zh) | 2018-06-01 |
Family
ID=57096161
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610338888.9A Expired - Fee Related CN106026191B (zh) | 2016-05-20 | 2016-05-20 | 用于最小化弃光电量光伏电站发电权置换电量评估的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106026191B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106558884B (zh) * | 2017-01-20 | 2019-03-19 | 许继集团有限公司 | 区域能源互联网控制系统有功功率平衡优化互补调节方法 |
| CN107294088B (zh) * | 2017-06-28 | 2023-08-01 | 中国能源建设集团甘肃省电力设计院有限公司 | 用于电改的增量配电网构建应用方法与多电源直供电系统 |
| CN110378787B (zh) * | 2019-06-27 | 2022-11-11 | 国电南瑞科技股份有限公司 | 一种交易电量的分解方法、装置及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103530696A (zh) * | 2013-09-27 | 2014-01-22 | 国家电网公司 | 计入风电不确定度的火电机组中长期时序递进发电协调算法 |
| CN103824235A (zh) * | 2014-03-18 | 2014-05-28 | 国家电网公司 | 一种基于竞价上网的新能源电站电量消纳方案制定方法 |
| JP2015167461A (ja) * | 2014-03-04 | 2015-09-24 | 日本電信電話株式会社 | 太陽光発電システムの制御方法 |
| CN105470987A (zh) * | 2015-12-23 | 2016-04-06 | 华北电力大学 | 荷-源协调控制的可再生能源日前计划修正方法 |
-
2016
- 2016-05-20 CN CN201610338888.9A patent/CN106026191B/zh not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103530696A (zh) * | 2013-09-27 | 2014-01-22 | 国家电网公司 | 计入风电不确定度的火电机组中长期时序递进发电协调算法 |
| JP2015167461A (ja) * | 2014-03-04 | 2015-09-24 | 日本電信電話株式会社 | 太陽光発電システムの制御方法 |
| CN103824235A (zh) * | 2014-03-18 | 2014-05-28 | 国家电网公司 | 一种基于竞价上网的新能源电站电量消纳方案制定方法 |
| CN105470987A (zh) * | 2015-12-23 | 2016-04-06 | 华北电力大学 | 荷-源协调控制的可再生能源日前计划修正方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106026191A (zh) | 2016-10-12 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104463378B (zh) | 用于水电资源丰富地区省级电网风力发电消纳能力的算法 | |
| Fang et al. | Look-ahead bidding strategy for concentrating solar power plants with wind farms | |
| CN105006844A (zh) | 一种间歇式发电并网条件下的电力系统日前鲁棒调度系统 | |
| CN105932668B (zh) | 用于最小化弃风弃光电量的新能源发电权置换评估的方法 | |
| CN111445107A (zh) | 冷热电联供型微电网多目标优化配置方法 | |
| CN114676991B (zh) | 基于源荷双侧不确定的多能互补系统优化调度方法 | |
| CN110956344A (zh) | 一种计及绿证与碳交易制度的“源-荷”协调优化调度方法 | |
| CN106026191B (zh) | 用于最小化弃光电量光伏电站发电权置换电量评估的方法 | |
| CN115115087B (zh) | 一种考虑氢燃料汽车与氢储能的虚拟电厂协调调度方法 | |
| CN110932261A (zh) | 一种基于全局效益最大化的多能源系统联合装机规划方法 | |
| CN106257792B (zh) | 一种基于两阶段机组组合的新能源优先调度方法 | |
| CN105896535B (zh) | 用于最小化弃风电量的风电场发电权置换电量评估的方法 | |
| CN116599148A (zh) | 一种面向新能源消纳的氢电混合储能两阶段协同规划方法 | |
| CN105939013B (zh) | 最小化弃风电量的风电场发电权置换电量评估方法 | |
| Su et al. | Capacity configuration optimization for green hydrogen generation driven by solar-wind hybrid power based on comprehensive performance criteria | |
| CN114266136A (zh) | 基于时序模拟和风险度量的清洁能源消纳能力评估方法 | |
| CN106712105A (zh) | 一种新能源接纳空间离散概率序列计算方法 | |
| CN116914732A (zh) | 基于深度强化学习的热电联产系统低碳调度方法及系统 | |
| CN106339777A (zh) | 一种发电机组检修计划优化的最小失负荷期望模型计算方法 | |
| Yao et al. | Determination of a dispatch strategy to maximize income for a wind turbine-BESS power station | |
| CN105894134B (zh) | 最小化弃风弃光电量的新能源发电权置换电量评估方法 | |
| CN114764652A (zh) | 一种考虑中长期调度的多周期协调电力平衡系统及方法 | |
| Song et al. | Capacity investment decisions of energy storage power stations supporting wind power projects | |
| CN105976267A (zh) | 最小化弃光电量的光伏电站发电权置换电量评估方法 | |
| CN113052396A (zh) | 一种适用于大规模风电和配置储热的热电联产机组的能量共享方法 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20180601 Termination date: 20200520 |