CN114614495A - 一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法 - Google Patents

一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法 Download PDF

Info

Publication number
CN114614495A
CN114614495A CN202210297367.9A CN202210297367A CN114614495A CN 114614495 A CN114614495 A CN 114614495A CN 202210297367 A CN202210297367 A CN 202210297367A CN 114614495 A CN114614495 A CN 114614495A
Authority
CN
China
Prior art keywords
power
irrigation system
pumping
storage
power generation
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.)
Pending
Application number
CN202210297367.9A
Other languages
English (en)
Inventor
陈冲
喻飞
张鹏
刘凯
李成子
张旭
陈昕
袁博
张发印
徐峰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Changjiang Institute of Survey Planning Design and Research Co Ltd
Original Assignee
Changjiang Institute of Survey Planning Design and Research Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Changjiang Institute of Survey Planning Design and Research Co Ltd filed Critical Changjiang Institute of Survey Planning Design and Research Co Ltd
Priority to CN202210297367.9A priority Critical patent/CN114614495A/zh
Publication of CN114614495A publication Critical patent/CN114614495A/zh
Pending legal-status Critical Current

Links

Images

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/381—Dispersed generators
    • 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/003—Load forecast, e.g. methods or systems for forecasting future load demand
    • 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/004—Generation forecast, e.g. methods or systems for forecasting future energy generation
    • 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/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
    • H02J3/0075—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources for providing alternative feeding paths between load and source according to economic or energy efficiency considerations, e.g. economic dispatch
    • 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]
    • 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
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/22—The renewable source being solar energy
    • H02J2300/24—The renewable source being solar energy of photovoltaic origin
    • 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
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy
    • 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
    • H02J2310/00—The network for supplying or distributing electric power characterised by its spatial reach or by the load
    • H02J2310/50—The network for supplying or distributing electric power characterised by its spatial reach or by the load for selectively controlling the operation of the loads
    • H02J2310/56—The network for supplying or distributing electric power characterised by its spatial reach or by the load for selectively controlling the operation of the loads characterised by the condition upon which the selective controlling is based
    • H02J2310/62—The condition being non-electrical, e.g. temperature
    • H02J2310/64—The condition being economic, e.g. tariff based load management

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明提出了一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法，首先调研获取电网第二天的风力、日照强度和常规负荷用电功率，建立环境因子多维度影响约束，计算得到风电发电功率和光电发电功率，预测计算第二天的灌溉用水量；继而，分析抽蓄灌溉系统的运行特性，建立抽蓄灌溉系统的运行模型，构建常规发电机组运行模型和电网功率平衡约束；在上述基础上，以系统日运行成本最小为目标，对抽蓄灌溉系统进行运行控制。本发明的运行控制方法充分考虑了风力和日照强度气象因子对光伏发电、风电、灌溉用水量以及常规发电机组的多维度影响，显著提高了系统运行经济性，在工程应用方面也具有丰富的参考价值。

Description

一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法
技术领域
本发明属于电力与水利工程技术领域，涉及一种抽蓄灌溉系统的运行控制方法。
背景技术
在能源危机和环保需求的激励下，高渗透率的风电和光伏发电并入电网。考虑到可再生能源发电的随机性、波动性和反调峰特性，部分地区已经建成抽水蓄能电站参与电网削峰填谷。可见，如何优化抽水蓄能电站工作，提高结合风电和光伏发电的抽蓄灌溉系统的运行经济性，具有重要的现实意义。
目前，对结合风电和光伏发电的抽蓄灌溉系统的运行控制中，鲜有考虑当地居民从水库汲水的灌溉需求，居民从水库汲水将导致水库的水量降低，影响抽水蓄能电站发电。并且，现有的运行控制方法中普遍忽视了日照强度、风力等气象因子的多维度影响，日照强度会影响光伏发电量，风力会影响风电发电量，进而还会影响灌溉用水量，进而影响抽水蓄能电站工作。
发明内容
为解决背景技术中所述的问题，本发明提出了一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法。
本发明的技术方案包括以下步骤：
获取电网第二天的风力、日照强度和常规负荷用电功率；
依据风力和日照强度计算得到风电发电功率和光电发电功率；
通过过往同一时节的灌溉用水量与日照强度的调研数据，预测计算第二天的灌溉用水量；
依据抽蓄灌溉系统的运行特性，建立抽蓄灌溉系统的工作模型、常规发电机组运行模型以及电网发电功率和用电功率的平衡约束，建立以系统日运行成本最小为目标函数的抽蓄灌溉系统运行控制模型；
对抽蓄灌溉系统运行控制模型进行求解，得到常规发电机组日前运行功率值、抽水蓄能电站日前运行功率值以及抽水负荷日前运行功率值；
依据得到的常规发电机组日前运行功率值、抽水蓄能电站日前运行功率值以及抽水负荷日前运行功率值，对抽蓄灌溉系统进行运行控制。
进一步地，所述的电网第二天的风力、日照强度和常规负荷用电功率表示为集合A：
Figure BDA0003562094670000021
其中，v1,v2,...,vt表示第二天各个时刻的风速，G1,G2,...,Gt表示第二天各个时刻的日照强度，P1 L,P2 L,...,Pt L表示第二天各个时刻的常规负荷用电功率。
进一步地，所述的风电发电功率和光电发电功率为：
Figure BDA0003562094670000031
Pt PV＝PSTCGt[1+ε(Tc,t-Tr)]/GSTC，
其中，Pt WT是风电t时刻发电功率；υt是t时刻风速；
Figure BDA0003562094670000032
υr和
Figure BDA0003562094670000033
分别代表切入风速、额定风速和切出风速；Pt PV是光伏发电t时刻发电功率；PSTC为标准测试条件下单位光伏装机容量下输出的最大测试功率；GSTC为标准测试条件下的日照强度；ε为功率温度系数；Tc,t为t时刻电池板的工作温度；Tr为参考温度。
进一步地，所述的第二天的灌溉用水量的预测计算方法为：
①调取过往同一时节的灌溉用水量和日照强度的数据，在二维坐标系上绘制相关曲线；
②基于MATLAB平台利用多项式拟合日照强度-灌溉用水量曲线，确定系数ai(i＝1,2,...,5)；
③依据日照强度-灌溉用水量的拟合曲线和获取的电网第二天的日照强度，得到第二天的灌溉用水量为：
Qt＝a1Gt+a2(Gt)2+a3(Gt)3+a4(Gt)4+a5(Gt)5，
其中，Qt为第二天t时刻的灌溉用水量，Gt为第二天t时刻的日照强度。
进一步地，所述的抽蓄灌溉系统的工作模型中，抽蓄灌溉系统的抽水蓄能运行模式对应的功率约束为：
Figure BDA0003562094670000034
抽蓄灌溉系统的发电运行模式对应的功率约束为：
Figure BDA0003562094670000041
抽蓄灌溉系统在任意时刻只能工作在抽水蓄能运行模式或发电运行模式的约束为：
Figure BDA0003562094670000042
抽蓄灌溉系统的爬坡速率约束为：
-ΔPWP≤Pt WP,c-Pt WP,d≤ΔPWP；
抽蓄灌溉系统的水库储水值与抽蓄电站工作功率和灌溉用水量的关系为：
Figure BDA0003562094670000043
抽蓄灌溉系统的水库容量限制约束为：
Figure BDA0003562094670000044
其中，Pt WP,c和Pt WP,d分别为抽水蓄能电站t时刻的抽水蓄能功率和发电功率，
Figure BDA0003562094670000045
表示抽蓄电站额定工作功率，布尔型变量
Figure BDA0003562094670000046
和
Figure BDA0003562094670000047
分别表示抽水蓄能电站t时刻的工作模式，ΔPWP表示抽水蓄能电站爬坡速率，ηc和ηd分别是抽蓄电站的能水转化和水能转换效率，Wt Lake、
Figure BDA0003562094670000048
和
Figure BDA0003562094670000049
分别表示水库t时刻储水、最小储水和最大储水值。
进一步地，所述的常规发电机组运行模型，其中常规发电机组的工作范围为：
Figure BDA00035620946700000410
常规发电机组的爬坡速率约束为：
Figure BDA00035620946700000411
其中，Pt MT表示常规发电机组t时刻运行功率，
Figure BDA0003562094670000051
表示常规发电机组最大允许运行功率，ΔPMT表示常规发电机组额定爬坡速率；
所述的电网发电功率和用电功率的平衡约束为：
Pt WT+Pt PV+Pt MT+Pt WP,d-Pt WP,c＝Pt L+Pt SL，
其中，Pt SL表示t时刻灌溉负荷，
Figure BDA0003562094670000052
Figure BDA0003562094670000053
表示灌溉负荷额定运行功率，
Figure BDA0003562094670000054
ESL表示灌溉负荷日总用电量。
进一步地，所述的抽蓄灌溉系统运行控制模型以系统日运行成本最小为优化目标，目标函数可表示为：
Figure BDA0003562094670000055
Figure BDA0003562094670000056
其中，
Figure BDA0003562094670000057
表示t时刻常规发电机组运行成本，b3、b2、b1和b0为常规发电机发电特性参数。
进一步地，所述的抽蓄灌溉系统运行控制模型的求解方法为：
①定义函数f(Pt MT)满足下式：
f(Pt MT)＝b3(Pt MT)3+b2(Pt MT)2+b1(Pt MT)+b0，
其中，Pt MT表示常规发电机组t时刻运行功率；
②根据常规发电机组运行功率范围
Figure BDA0003562094670000058
等间距设置K+1个断点功率，其满足如下关系：
Figure BDA0003562094670000059
③对于任意功率Pt MT，可以用不同权重下的断点功率叠加组合表示：
Figure BDA0003562094670000061
其中，
Figure BDA0003562094670000062
表示断点k对应的功率权重；
④常规发电机组的运行成本可以表示为：
Figure BDA0003562094670000063
然后基于MATALB平台编码刻画优化问题，采用yalmip软件调用商业求解器Gurobi求解得到发电机组运行成本最小时的运行功率值。
进一步地，所述的断点k对应的功率权重
Figure BDA0003562094670000064
满足以下约束：
①常规发电机组实际运行功率只能落在单一功率段之内：
Figure BDA0003562094670000065
②断点功率点权重之和等于1：
Figure BDA0003562094670000066
③线性组合时只能允许相邻固定功率点参与：
Figure BDA0003562094670000067
Figure BDA0003562094670000068
Figure BDA0003562094670000069
Figure BDA0003562094670000072
其中，布尔型变量
Figure BDA0003562094670000071
用于标识常规发电机组具体的运行功率区间。
本发明与现有技术相比，所提出的抽蓄灌溉系统的运行控制方法中，首先调研获取电网第二天的风力、日照强度和常规负荷用电功率，建立环境因子多维度影响约束，计算得到风电发电功率和光电发电功率，预测计算第二天的灌溉用水量；继而，量化分析抽蓄灌溉系统的运行特性，建立抽蓄灌溉系统的运行模型；此外，构建发电机组运行模型；在上述基础上，以系统日运行成本最小为目标，对抽蓄灌溉系统进行运行控制。
本发明的运行控制方法充分考虑了风力和日照强度气象因子对光伏发电、风电、灌溉用水量以及常规发电机组的多维度影响，以系统日运行成本最小为目标，对抽蓄灌溉系统进行运行控制，显著提高了抽蓄灌溉系统运行经济性，在工程应用方面也具有丰富的参考价值。
附图说明
图1为本发明实施例中获取的电网第二天的风力-时间曲线图、光照强度-时间曲线图；
图2为本发明实施例中灌溉负荷功率-时间曲线图。
图3为本发明抽蓄灌溉系统的各运行功率-时间变化图。
具体实施方式
下面结合附图详细说明本发明的实施情况，但它们并不构成对本发明的限定，仅做举例而已，同时通过说明，将更加清楚地理解本发明的优点。本领域的普通的技术人员能从本发明公开的内容直接导出或联想到的所有变形，均应认为是本发明的保护范围。实施例中所述的位置关系均与附图所示一致，实施例中其他未详细说明的部分均为现有技术。
1、调研获取电网第二天的风力、日照强度和常规负荷用电功率，并用集合A表示：
Figure BDA0003562094670000081
其中，v1,v2,...,vt表示第二天各个时刻的风速，G1,G2,...,Gt表示第二天各个时刻的日照强度，P1 L,P2 L,...,Pt L表示第二天各个时刻的常规负荷用电功率。
在本实施例中，获取的电网第二天的风力-时间曲线图、光照强度-时间曲线图如图1所示。
2、量化分析气象因子作用，建立环境因子多维度影响约束。
首先，根据风力发电和光伏发电公式，计算得到日前可再生能源发电功率为：
Figure BDA0003562094670000082
Pt PV＝PSTCGt[1+ε(Tc,t-Tr)]/GSTC ③，
其中，Pt WT是风电t时刻发电功率；υt是t时刻风速；
Figure BDA0003562094670000083
υr和
Figure BDA0003562094670000084
代表切入风速、额定风速和切出风速；Pt PV是光伏t时刻发电功率；PSTC为标准测试条件下单位光伏装机容量下输出的最大测试功率；GSTC为标准测试条件下的日照强度；ε为功率温度系数；Tc,t为t时刻电池板的工作温度；Tr为参考温度。
3、预测计算第二天的灌溉用水量。
首先调研过往同一时节的当地居民灌溉用水量和日照强度的数据，在二维坐标系上绘制相关曲线。
然后基于MATLAB平台利用多项式拟合日照强度-灌溉用水量曲线，确定系数ai(i＝1,2,...,5)。
最后依据日照强度-灌溉用水量的拟合曲线和获取的电网第二天的日照强度，得到第二天的灌溉用水量为：
Qt＝a1Gt+a2(Gt)2+a3(Gt)3+a4(Gt)4+a5(Gt)5 ④，
其中，Qt为第二天t时刻的灌溉用水量，Gt为第二天t时刻的日照强度。
4、依据抽蓄灌溉系统的运行特性，建立抽蓄灌溉系统的工作模型。
考虑到目前抽蓄灌溉系统的运行模型较为复杂，给求解带来了巨大的麻烦，为此本发明提出一个简化的抽蓄灌溉系统工作模型，如式(5)～(10)所示。其中，式(5)和(6)表示抽蓄蓄能电站的抽水蓄能和发电两种运行模式对应的功率约束，式(7)限制了抽水蓄能电站任意时刻只能工作在抽蓄状态或发电状态，式(8)表示抽水蓄能电站爬坡速率约束，式(9)表示水库储水值与抽蓄电站工作功率和灌溉用水量的关系，式(10)表示水库容量限制。
Figure BDA0003562094670000091
Figure BDA0003562094670000092
Figure BDA00035620946700001013
-ΔPWP≤Pt WP,c-Pt WP,d≤ΔPWP ⑧，
Figure BDA00035620946700001014
Figure BDA0003562094670000101
其中，Pt WP,c和Pt WP,d分别为抽水蓄能电站t时刻的抽水蓄能功率和发电功率，
Figure BDA0003562094670000102
表示抽蓄电站额定工作功率，布尔型变量It WP,c和It WP,d分别表示抽水蓄能电站t时刻的工作模式，ΔPWP表示抽水蓄能电站爬坡速率，ηc和ηd分别是抽蓄电站的能水转化和水能转换效率，Wt Lake、
Figure BDA0003562094670000103
和
Figure BDA0003562094670000104
分别表示水库t时刻储水、最小储水和最大储水值。
5、构建常规发电机组运行模型，以及电网发电功率和用电功率的平衡约束。
首先，构建常规发电机组运行模型，如式
Figure BDA0003562094670000108
和
Figure BDA0003562094670000109
所示，式
Figure BDA00035620946700001010
表示常规发电机组要在允许范围内工作，式
Figure BDA00035620946700001011
表示常规发电机组的爬坡速率约束：
Figure BDA0003562094670000105
Figure BDA0003562094670000106
式中，Pt MT表示常规发电机组t时刻运行功率，
Figure BDA0003562094670000107
表示常规发电机组最大允许运行功率，ΔPMT表示常规发电机组额定爬坡速率。
继而，建立电网发电功率和用电功率平衡约束：
Pt WT+Pt PV+Pt MT+Pt WP,d-Pt WP,c＝Pt L+Pt SL
Figure BDA00035620946700001012
其中，Pt SL表示t时刻灌溉负荷，其可以进行需求侧响应参与电网削峰填谷，其相应的运行约束可以表示为：
Figure BDA0003562094670000111
Figure BDA0003562094670000112
其中，Pt SL表示t时刻灌溉负荷，ESL表示灌溉负荷日总用电量。
本实施例中灌溉负荷功率-时间曲线图如图2所示。
6、建立以系统日运行成本最小为目标函数的抽蓄灌溉系统运行控制模型。
系统日运行成本只考虑常规发电机组的燃料成本，系统日运行成本最小即为常规发电机组运行成本最小。目前往往使用二次函数拟合常规发电机组的燃料消耗，本发明采用三次多项式，拟合常规发电机组的运行成本特性，得到抽蓄灌溉系统运行控制模型的目标函数为：
Figure BDA0003562094670000113
Figure BDA0003562094670000114
其中，
Figure BDA0003562094670000115
表示t时刻常规发电机组运行成本，Pt MT表示常规发电机组t时刻运行功率，b3、b2、b1和b0为常规发电机发电特性参数。
7、对目标函数
Figure BDA0003562094670000118
和
Figure BDA0003562094670000119
进行求解，得到抽蓄灌溉系统日运行成本最小时的常规发电机组日前运行功率值、抽水蓄能电站日前运行功率值以及抽水负荷日前运行功率值。
首先，定义函数f(Pt MT)满足下式：
f(Pt MT)＝b3(Pt MT)3+b2(Pt MT)2+b1(Pt MT)+b0
Figure BDA00035620946700001110
在此基础上，根据常规发电机组运行功率范围
Figure BDA0003562094670000116
等间距设置K+1个断点功率，其满足如下关系：
Figure BDA0003562094670000117
对于任意功率Pt MT，可以用不同权重下的断点功率叠加组合表示：
Figure BDA0003562094670000121
其中，
Figure BDA0003562094670000122
表示断点k对应的功率权重，其满足约束
Figure BDA0003562094670000123
式
Figure BDA00035620946700001211
表示常规发电机组实际运行功率只能落在单一功率段之内，式
Figure BDA00035620946700001212
表示断点功率点权重之和等于1，式
Figure BDA00035620946700001213
限制了线性组合时只能允许相邻固定功率点参与。
Figure BDA0003562094670000124
Figure BDA0003562094670000125
Figure BDA0003562094670000126
Figure BDA0003562094670000127
Figure BDA0003562094670000128
Figure BDA00035620946700001214
其中，布尔型变量
Figure BDA0003562094670000129
用于标识常规发电机组具体的运行功率区间。
最后，常规发电机组的运行成本可以表示为：
Figure BDA00035620946700001210
然后基于MATALB平台编码刻画优化问题，采用yalmip软件调用商业求解器Gurobi求解得到抽蓄灌溉系统日运行成本最小时的常规发电机组日前运行功率值、抽水蓄能电站日前运行功率值以及抽水负荷日前运行功率值。
在本实施例中，通过上述公式
Figure BDA0003562094670000131
可求解得到系统日运行成本最小值为82176万元。
8、抽蓄灌溉系统日运行成本最小时的常规发电机组日前运行功率值、抽水蓄能电站日前运行功率值以及抽水负荷日前运行功率值，对抽蓄灌溉系统进行运行控制。
采用上述所述的运行控制方法后，该抽蓄灌溉系统的各运行功率-时间变化图如图3所示，常规负荷功率与灌溉负荷功率之和即为抽蓄灌溉系统的总输出功率，抽蓄灌溉系统的总输入功率为风电和光电得到的可再生能源功率与常规发电机组功率之和，通过图3可得知，常规发电机组在一天中的运行时间极少，该抽蓄灌溉系统的运行功率来源主要依靠于运行成本极低的风电和光电得到的可再生能源功率。该运行控制方法充分考虑了风力和日照强度气象因子对光伏发电、风电、灌溉用水量以及常规发电机组的多维度影响，在系统日运行成本最小时的前提下实现了抽蓄灌溉系统的正常运行，显著提高了抽蓄灌溉系统运行经济性。
以上结合附图及具体实施例详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种简单变型，这些简单变型均属于本发明的保护范围。

Claims (9)

1.一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法，其特征在于，包括以下步骤：
获取电网第二天的风力、日照强度和常规负荷用电功率；
依据风力和日照强度计算得到风电发电功率和光电发电功率；
通过过往同一时节的灌溉用水量与日照强度的调研数据，预测计算第二天的灌溉用水量；
依据抽蓄灌溉系统的运行特性，建立抽蓄灌溉系统的工作模型、常规发电机组运行模型以及电网发电功率和用电功率的平衡约束，建立以系统日运行成本最小为目标函数的抽蓄灌溉系统运行控制模型；
对抽蓄灌溉系统运行控制模型进行求解，得到常规发电机组日前运行功率值、抽水蓄能电站日前运行功率值以及抽水负荷日前运行功率值；
依据得到的常规发电机组日前运行功率值、抽水蓄能电站日前运行功率值以及抽水负荷日前运行功率值，对抽蓄灌溉系统进行运行控制。
2.根据权利要求1所述的一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法，其特征在于：所述的电网第二天的风力、日照强度和常规负荷用电功率表示为集合A：
Figure FDA0003562094660000011
其中，v1,v2,...,vt表示第二天各个时刻的风速，G1,G2,...,Gt表示第二天各个时刻的日照强度，
Figure FDA0003562094660000021
表示第二天各个时刻的常规负荷用电功率。
3.根据权利要求2所述的一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法，其特征在于：所述的风电发电功率和光电发电功率为：
Figure FDA0003562094660000022
Pt PV＝PSTCGt[1+ε(Tc,t-Tr)]/GSTC，
其中，Pt WT是风电t时刻发电功率；υt是t时刻风速；
Figure FDA0003562094660000023
υr和
Figure FDA0003562094660000024
分别代表切入风速、额定风速和切出风速；Pt PV是光伏发电t时刻发电功率；PSTC为标准测试条件下单位光伏装机容量下输出的最大测试功率；GSTC为标准测试条件下的日照强度；ε为功率温度系数；Tc,t为t时刻电池板的工作温度；Tr为参考温度。
4.根据权利要求3所述的一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法，其特征在于：所述的第二天的灌溉用水量的预测计算方法为：
①调取过往同一时节的灌溉用水量和日照强度的数据，在二维坐标系上绘制相关曲线；
②基于MATLAB平台利用多项式拟合日照强度-灌溉用水量曲线，确定系数ai(i＝1,2,...,5)；
③依据日照强度-灌溉用水量的拟合曲线和获取的电网第二天的日照强度，得到第二天的灌溉用水量为：
Qt＝a1Gt+a2(Gt)2+a3(Gt)3+a4(Gt)4+a5(Gt)5，
其中，Qt为第二天t时刻的灌溉用水量，Gt为第二天t时刻的日照强度。
5.根据权利要求4所述的一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法，其特征在于：所述的抽蓄灌溉系统的工作模型中，抽蓄灌溉系统的抽水蓄能运行模式对应的功率约束为：
Figure FDA0003562094660000031
抽蓄灌溉系统的发电运行模式对应的功率约束为：
Figure FDA0003562094660000032
抽蓄灌溉系统在任意时刻只能工作在抽水蓄能运行模式或发电运行模式的约束为：
Figure FDA0003562094660000033
抽蓄灌溉系统的爬坡速率约束为：
-ΔPWP≤Pt WP,c-Pt WP,d≤ΔPWP；
抽蓄灌溉系统的水库储水值与抽蓄电站工作功率和灌溉用水量的关系为：
Figure FDA0003562094660000034
抽蓄灌溉系统的水库容量限制约束为：
Figure FDA0003562094660000035
其中，Pt WP,c和Pt WP,d分别为抽水蓄能电站t时刻的抽水蓄能功率和发电功率，
Figure FDA0003562094660000036
表示抽蓄电站额定工作功率，布尔型变量
Figure FDA0003562094660000037
和
Figure FDA0003562094660000038
分别表示抽水蓄能电站t时刻的工作模式，ΔPWP表示抽水蓄能电站爬坡速率，ηc和ηd分别是抽蓄电站的能水转化和水能转换效率，Wt Lake、
Figure FDA0003562094660000039
和
Figure FDA00035620946600000310
分别表示水库t时刻储水、最小储水和最大储水值。
6.根据权利要求5所述的一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法，其特征在于：所述的常规发电机组运行模型中，其工作范围为：
Figure FDA0003562094660000041
常规发电机组的爬坡速率约束为：
Figure FDA0003562094660000042
其中，Pt MT表示常规发电机组t时刻运行功率，
Figure FDA0003562094660000043
表示常规发电机组最大允许运行功率，ΔPMT表示常规发电机组额定爬坡速率；
电网发电功率和用电功率的平衡约束为：
Pt WT+Pt PV+Pt MT+Pt WP,d-Pt WP,c＝Pt L+Pt SL，
其中，Pt SL表示t时刻灌溉负荷，
Figure FDA0003562094660000044
Figure FDA0003562094660000045
表示灌溉负荷额定运行功率，
Figure FDA0003562094660000046
ESL表示灌溉负荷日总用电量。
7.根据权利要求6所述的一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法，其特征在于：所述的抽蓄灌溉系统运行控制模型的目标函数为：
Figure FDA0003562094660000047
Figure FDA0003562094660000048
其中，
Figure FDA0003562094660000049
表示t时刻常规发电机组运行成本，b3、b2、b1和b0为常规发电机发电特性参数。
8.根据权利要求7所述的一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法，其特征在于：所述的抽蓄灌溉系统运行控制模型的求解方法为：
①定义函数f(Pt MT)满足下式：
f(Pt MT)＝b3(Pt MT)3+b2(Pt MT)2+b1(Pt MT)+b0，
其中，Pt MT表示常规发电机组t时刻运行功率；
②根据常规发电机组运行功率范围
Figure FDA0003562094660000051
等间距设置K+1个断点功率，其满足如下关系：
Figure FDA0003562094660000052
③对于任意功率Pt MT，可以用不同权重下的断点功率叠加组合表示：
Figure FDA0003562094660000053
其中，
Figure FDA0003562094660000054
表示断点k对应的功率权重；
④常规发电机组的运行成本可以表示为：
Figure FDA0003562094660000055
然后基于MATALB平台编码刻画优化问题，采用yalmip软件调用商业求解器Gurobi求解得到系统日运行成本最小时的常规发电机组、抽水蓄能电站以及抽水负荷日前运行功率值。
9.根据权利要求8所述的一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法，其特征在于：所述的断点k对应的功率权重
Figure FDA0003562094660000061
满足以下约束：
①常规发电机组实际运行功率只能落在单一功率段之内：
Figure FDA0003562094660000062
②断点功率点权重之和等于1：
Figure FDA0003562094660000063
③线性组合时只能允许相邻固定功率点参与：
Figure FDA0003562094660000064
Figure FDA0003562094660000065
Figure FDA0003562094660000066
Figure FDA0003562094660000067
其中，布尔型变量
Figure FDA0003562094660000068
用于标识常规发电机组具体的运行功率区间。
CN202210297367.9A 2022-03-24 2022-03-24 一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法 Pending CN114614495A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210297367.9A CN114614495A (zh) 2022-03-24 2022-03-24 一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210297367.9A CN114614495A (zh) 2022-03-24 2022-03-24 一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法

Publications (1)

Publication Number Publication Date
CN114614495A true CN114614495A (zh) 2022-06-10

Family

ID=81863960

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210297367.9A Pending CN114614495A (zh) 2022-03-24 2022-03-24 一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法

Country Status (1)

Country Link
CN (1) CN114614495A (zh)

Similar Documents

Publication Publication Date Title
CN109508499B (zh) 多时段多场景分布式电源最优接入位置与容量研究方法
CN103023035B (zh) 一种多能互补发电系统的优化配置方法
CN113343478B (zh) 考虑不确定性和需求响应的独立微网容量优化配置方法
CN102419394B (zh) 一种预测分辨率可变的风光功率预测方法
CN101515722A (zh) 用于提高功率受限的风力发电厂的性能的系统和方法
CN108879793A (zh) 一种风光蓄能水电站互补的离网式混合能源系统优化方法
CN110601260B (zh) 一种限定联络线上功率波动的光-蓄系统容量优化方法
CN112600209A (zh) 一种含潮流能的海岛独立微电网容量多目标优化配置方法
Tian et al. A coordinated multi-time scale robust scheduling framework for isolated power system with ESU under high RES penetration
CN106600022B (zh) 基于多目标优化的风-光-燃气-海水抽水蓄能孤立电力系统容量优化配置方法
CN112769156B (zh) 一种计及大规模海上风电并网的源网荷储协调运行方法
CN115313508A (zh) 一种微电网储能优化配置方法、装置及存储介质
CN115333161A (zh) 一种绿色水厂供电系统的容量优化配置方法
CN113536581B (zh) 一种计及运行策略的储能系统多状态可靠性建模方法
CN110768303A (zh) 一种海岛型能源系统设备容量的优化配置方法
CN116914860B (zh) 多能互补发电系统电源容量配置优化方法及系统
CN110571868B (zh) 微电网的优化配置方法
CN117592690A (zh) 一种水电站群迎峰度夏顶峰能力分析方法
CN117410958A (zh) 一种基于urfc的微电网容量配置优化方法
Hammad et al. Numerical Analysis Of PV Water Pumping System In Different Geographical Locations In Egypt.
CN108062606A (zh) 一种基于黎曼积分的虚拟电厂调度优化方法
CN116562554A (zh) 一种考虑灵活性供需平衡的长短期储能联合规划方法
CN114614495A (zh) 一种配置风电和光伏发电的抽蓄灌溉系统的运行控制方法
Han et al. Analysis of economic operation model for virtual power plants considering the uncertainties of renewable energy power generation
CN110544958A (zh) 一种电力系统消纳随机出力电源能力的确定方法及装置

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination