CN108133104A - 一种长期跨流域多梯级水电优化运行模拟方法 - Google Patents
一种长期跨流域多梯级水电优化运行模拟方法 Download PDFInfo
- Publication number
- CN108133104A CN108133104A CN201711397656.1A CN201711397656A CN108133104A CN 108133104 A CN108133104 A CN 108133104A CN 201711397656 A CN201711397656 A CN 201711397656A CN 108133104 A CN108133104 A CN 108133104A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- hydropower
- long
- unit
- term
- 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.)
- Granted
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 68
- 230000007774 longterm Effects 0.000 title claims abstract description 59
- 238000005457 optimization Methods 0.000 title claims abstract description 35
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000006243 chemical reaction Methods 0.000 claims abstract description 12
- 238000012545 processing Methods 0.000 claims abstract description 10
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 4
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 93
- 238000003860 storage Methods 0.000 claims description 38
- 239000011159 matrix material Substances 0.000 claims description 26
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000010248 power generation Methods 0.000 claims description 11
- 238000011144 upstream manufacturing Methods 0.000 claims description 10
- 238000012546 transfer Methods 0.000 claims description 8
- 239000008239 natural water Substances 0.000 claims description 6
- 238000001704 evaporation Methods 0.000 claims description 4
- 230000008020 evaporation Effects 0.000 claims description 4
- 238000010521 absorption reaction Methods 0.000 claims description 3
- 230000008901 benefit Effects 0.000 description 15
- 230000005611 electricity Effects 0.000 description 4
- 238000010276 construction Methods 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000012886 linear function Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000002351 wastewater Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A20/00—Water conservation; Efficient water supply; Efficient water use
- Y02A20/152—Water filtration
-
- 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
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E40/00—Technologies for an efficient electrical power generation, transmission or distribution
- Y02E40/70—Smart grids as climate change mitigation technology in the energy generation sector
-
- 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
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Human Resources & Organizations (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Public Health (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种长期跨流域多梯级水电优化运行模拟方法,其步骤:构建由目标函数及约束条件组成的长期水电群运行模拟模型;建立长期水电优化运行模拟模型的决策变量;建立长期水电群运行模拟模型的目标函数;建立长期水电群运行模拟模型的约束条件;对长期水电群运行模拟模型进行简化与求解:对水电出力转换函数进行非线性处理;处理水库与水库之间、水库与机组之间的空间拓扑联系;采用已有混合整数规划算法求解长期水电群运行模拟模型,实现模拟长期跨流域多梯级水电优化运行。本发明可为水电比例较大的电力系统中水电群的优化调度提供技术支撑。
Description
技术领域
本发明涉及一种梯级水电的电力系统优化运行模拟领域,特别是关于一种长期跨流域多梯级水电优化运行模拟方法。
背景技术
我国河川众多,水能资源丰富。20世纪70年代末普查统计,水能资源理论总蕴藏量为5.92万亿千瓦时/年,居世界第一。本世纪以来,随着我国经济的高速发展,能源需求日益增加,水能资源得到了快速的开发建设,截止至2016年,我国水电装机总容量为3.32亿kW,居世界第一位。我国快速的水电建设增长给电网对水电的合理消纳、梯级水电站群的合理调度运行提出了挑战。
在电力系统中,与其他化石能源相比,水电资源不仅是碳排放低、可再生的清洁能源,更由于其调节能力灵活、机组启停快、发电成本低廉等特点,因此是电力系统电源结构中重要的组成部分。在电力系统中,水电系统不仅承担着部分电力负荷的供应,还参与了系统的调频调峰、系统备用等任务。在电力系统生产运行中,如何利用水电的特点,优化水电站群的调度运行方式,尽量减少水库弃水,充分发挥水电的经济效益,进而一方面有效降低全系统的运行成本,另一方面提高电力系统的安全可靠性具有重大意义。
在以全系统发电运行成本最小化的优化运行模拟背景下,水电的经济效益取决于以下几方面因素:防止弃水、调峰效益、水库水头效益、同一水头下水电机组的出力效率效益。在四个因素中防止弃水的效益最大,是优化运行计划中首先考虑的因素,其次是调峰效益,然后水头效益次之,效率效益又次之。单独强调系统中比如火电机组的效益或者单独强调水电机组的效益,在经济上都得不到最优运行状态,尤其当系统电源结构中水电比例相对较大时,该问题更加明显。因此,在考虑大规模水电群的全系统优化运行模拟时,更需要合理协调水电机组与其他类型机组的效益,才能实现整个电力系统最大经济效益。
发明内容
针对上述问题,本发明的目的是提供一种长期跨流域多梯级水电优化运行模拟方法,其能实现减小发电弃水、提高发电平均水头,最终使得水电群的综合效益最大化。
为实现上述目的,本发明采取以下技术方案:一种长期跨流域多梯级水电优化运行模拟方法,其特征在于包括以下步骤:1)构建由目标函数及约束条件组成的长期水电群运行模拟模型;1.1)建立长期水电优化运行模拟模型的决策变量;1.2)建立长期水电群运行模拟模型的目标函数;1.3)建立长期水电群运行模拟模型的约束条件;2)对长期水电群运行模拟模型进行简化与求解:2.1)对水电出力转换函数进行非线性处理;2.2)处理水库与水库之间、水库与机组之间的空间拓扑联系;2.3)采用已有混合整数规划算法求解长期水电群运行模拟模型,实现模拟长期跨流域多梯级水电优化运行。
进一步,所述步骤1.1)中,长期水电群运行模拟模型中的决策变量分为五种:Pi t为水电机组i在t时段的出力,单位MW;为水电机组i在t时段的发电用水量,单位m3;为水库k在t时段的库容,单位m3;为水库k在t时段的用水量,单位m3;为水库k在t时段的弃水量,单位m3;
进一步,所述步骤1.2)中,目标函数为:
其中,E为水电系统在模拟优化运行期内所有水电机组发电量,单位MW;Δt为系统每个计算时段的长度,单位小时;L为目标优化时段数,模型以月为计算时段;N为水电系统的水电机组数。
进一步,所述步骤1.3)中,具体包括:1.3.1)建立水电机组运行约束:1.3.1.1)确定水电出力上下限约束:
Pi,min≤Pi t≤Pi,max
其中,Pi,min为水电机组i最小出力,单位MW;Pi,max为水电机组i最大出力,单位MW;1.3.1.2)确定水电转换函数:
其中,水电机组i处于水库k中;1.3.1.3)确定水电系统的水电消纳能力约束:长期水电群运行模拟模型设定系统消纳水电出力的总量不能超过系统负荷预设比例范围:
其中,γ表示水电出力占系统t计算时段全系统负荷的最大比例、λ表示水电出力占系统t计算时段全系统负荷的最小比例;Dt为t时段全系统社会负荷,单位MWh;1.3.2)建立水库运行约束。
进一步,所述步骤1.3.2)中,具体包括:1.3.2.1)水电机组用水与水库用水平衡约束:
其中,Ωk是水库k包含的所有水电机组的集合;1.3.2.2)水库库容上下限约束:Vk,min≤Vk≤Vk,max,其中,Vk,min为水库k库容下限,单位m3;Vk为水库k的库容,单位m3;Vk,max为水库k库容上限,单位m3;1.3.2.3)水库库容始末约束:
其中,Vk,begin为水库k初始计算时段库容边界,单位m3;Vk,end为水库k终止计算时段库容边界,单位m3;1.3.2.4)水库用水上下限约束:Qk,min≤Qt k≤Qk,max,其中,Qk,min为水库k用水量下限,单位m3;Qk,max为水库k用水量上限,单位m3;Qk,min主要取决于下游航运、正常供水的限制;Qk,max主要取决于水库的泄流能力;1.3.2.5)水库弃水量约束:1.3.2.6)水量平衡约束:
其中,为水库k在t时段自然来水量,单位m3;Πk为水库k包含的所有上游水库的集合;为水库k时段t的蒸发率;λj-k表示水库j用水下泄至水库k的有效比例。
进一步,所述步骤2.1)中,采用离散库容的方法处理水电转换关系:在水库库容上下限范围内,将库容分成若干段,每段库容范围内用该段的平均库容进行计算;然后将机组出力转换为用水量的一元非线性函数,再进一步将出力-用水曲线线性化。
进一步,对库容的离散化处理如下:
表示第k个水库,离散为L段,每段的状态变量有且仅有一个状态变量为1,其他均为0:
水库库容离散为L段,第l段的范围为:[Vk,l-1,Vk,l],Vk,l-1为水库库容离散的第l-1段;Vk,l为水库库容离散的第l段;则有:
Vk,0=Vk,min
Vk,L=Vk,max
当库容处于第l段时,即时需满足如下约束:
得到离散库容约束为:
进一步,处于第k个水库的水电机组i的水电转换函数线性简化处理如下:
其中,λi,l表示第i台机组在第l段离散库容的发电效率;
设表示第i台机组在第l段离散库容的发电出力,则有:
其中M为数值很大的常量,当时,当时,
进一步,所述步骤2.2)中,采用矩阵的方式处理水库与机组之间、水库与水库之间的空间约束关系,设水库-机组关联矩阵AR-U表示水库与机组之间的空间拓扑关系,若机组i,j处于水库r上,则关联矢量为:
则水库-机组关联矩阵为:AR-U=[a1,a2,...,aR]T;
其中R表示水库,U表示机组;设水库-水库关联矩阵AR-R表示水库与水库之间的空间拓扑关系;若水库i是水库r的下游水库,则关联矢量为:
则水库-水库关联矩阵为:AR-R=[a1,a2,...,aR];
设则关联矢量a* R中非零元素所在位置就表示水库r的上游水库编号;则水库-机组用水平衡约束为:
其中,矩阵λR-R为元素组成的方阵;为水库R在t时段的用水量,单位m3;为水库R在t时段的弃水量,单位m3;
水量平衡约束为:
其中,Ik t为水库k在t时段自然来水量,单位m3;矩阵λR-R为元素组成的方阵。
本发明由于采取以上技术方案,其具有以下优点:1、本发明中水电群的长期优化运行模拟模型以月为计算时段,目标是给出具有月以上调节能力水库的优化运行过程,使得水电系统能尽量减小弃水,同时提高平均发电水头,以期使水电系统总发电量最大。2、本发明将水电群的复杂非线性模型与电力系统运行模拟模型相结合,对具有较多水电接入的电力系统的长期运行模拟具有重要价值。3、本发明针对具有调节能力的水库,在满足水库与水电机组的运行约束,保证正常航运、社会用水、防洪等条件下,对水库的来水和用水过程进行调节,实现尽量减小发电弃水、提高发电平均水头,最终使得水电群的综合效益最大化。
附图说明
图1是本发明的水电转换函数关系图;
图2是本发明的分段线性化的水电转换关系图;其中,λ1、λ2、λ3分别表示在第1段、第2段和第3段的离散库容发电效率。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明提供一种长期跨流域多梯级水电优化运行模拟方法,其包括以下步骤:
1)构建由目标函数及约束条件组成的长期水电群运行模拟模型;
长期水电群运行模拟模型的构建过程具体包括以下步骤:
1.1)建立长期水电优化运行模拟模型的决策变量;
该长期水电群运行模拟模型中的决策变量分为五种:Pi t为水电机组i在t时段的出力,单位MW;为水电机组i在t时段的发电用水量,单位m3;为水库k在t时段的库容,单位m3;为水库k在t时段的用水量,单位m3;为水库k在t时段的弃水量,单位m3;
1.2)建立长期水电群运行模拟模型的目标函数;
目标函数为考虑水电系统在模拟优化运行期内所有水电机组发电量最大化:
其中,E为水电系统在模拟优化运行期内所有水电机组发电量,单位MW;Δt为系统每个计算时段的长度,单位小时;L为目标优化时段数,模型以月为计算时段;N为水电系统的水电机组数。
1.3)建立长期水电群运行模拟模型的约束条件,具体包括:
1.3.1)建立水电机组运行约束,具体包括:
1.3.1.1)确定水电出力上下限约束:
Pi,min≤Pi t≤Pi,max (2)
其中,Pi,min为水电机组i最小出力,单位MW;Pi,max为水电机组i最大出力,单位MW。
1.3.1.2)确定水电转换函数:
水电机组的水电转换函数是一个复杂的多维非线性函数,在本实施例中,长期水电群运行模拟模型水电转换函数关系考虑为:水电出力是水电机组用水与水库库容的二元非线性关系:
其中,水电机组i处于水库k中。
1.3.1.3)确定水电系统的水电消纳能力约束:
在全系统水火优化运行模拟的背景下,考虑到火电机组的调峰能力,长期水电群运行模拟模型设定系统消纳水电出力的总量不能超过系统负荷预设比例范围:
其中,γ表示水电出力占系统t计算时段全系统负荷的最大比例、λ表示水电出力占系统t计算时段全系统负荷的最小比例;Dt为t时段全系统社会负荷,单位MWh。
1.3.2)建立水库运行约束,具体包括:
1.3.2.1)水电机组用水与水库用水平衡约束:
水库用水为处于该水库的所有水电机组的用水量和水库弃水量之和:
其中,Ωk是水库k包含的所有水电机组的集合。
1.3.2.2)水库库容上下限约束:
Vk,min≤Vk≤Vk,max (6)
其中,Vk,min为水库k库容下限,单位m3;Vk为水库k的库容,单位m3;Vk,max为水库k库容上限,单位m3。
1.3.2.3)水库库容始末约束:
其中,Vk,begin为水库k初始计算时段库容边界,单位m3;Vk,end为水库k终止计算时段库容边界,单位m3。
1.3.2.4)水库用水上下限约束:
其中,Qk,min为水库k用水量下限,单位m3;Qk,max为水库k用水量上限,单位m3;Qk,min主要取决于下游航运、正常供水的限制;Qk,max主要取决于水库的泄流能力。
1.3.2.5)水库弃水量约束:
1.3.2.6)水量平衡约束:
水库水量平衡为:水库当前时段库容等于上一时段库容与当前时段水库来水量、当前时段上游水库用水量之和,然后减去当前时段水库用水。
其中,为水库k在t时段自然来水量,单位m3;Πk为水库k包含的所有上游水库的集合。由于长期每个计算时段较长,因此需要考虑库容的蒸发量,上一时段库容对当前时段库容的水量贡献需要乘以系数其中为水库k时段t的蒸发率;还需要考虑上游水库用水下泄到下游水库的损失量,水量平衡约束中上游水库用水量对当前时段库容的水量贡献需要乘以系数λj-k,λj-k表示水库j用水下泄至水库k的有效比例;不需要考虑上下游水库之间的流水延迟。
2)对长期水电群运行模拟模型进行简化与求解:
2.1)对水电出力转换函数进行非线性处理;
鉴于精确考虑水电转换函数的复杂性,将长期水电群运行模拟模型水电转换函数关系视为:水电出力是水电机组用水与水库库容的二元非线性关系。
水电转换函数关系通常如图1、图2所示,本发明采用离散库容的方法处理水电转换关系。其主要方法为:在水库库容上下限范围内,将库容分成若干段,每段库容范围内用该段的平均库容进行计算。如此处理之后,将机组出力转换为用水量的一元非线性函数,然后进一步将出力-用水曲线线性化。其中:
①长期水电群运行模拟模型中对库容的离散化处理如下:
表示第k个水库,离散为L段,每段的状态变量显然有且仅有一个状态变量为1,其他均为0:
水库库容离散为L段,第l段的范围为:[Vk,l-1,Vk,l],Vk,l-1为水库库容离散的第l-1段;Vk,l为水库库容离散的第l段;显然有:
Vk,0=Vk,min (13)
Vk,L=Vk,max
当库容处于第l段时,即时需满足如下约束:
可以得到离散库容约束为:
②处于第k个水库的水电机组i的水电转换函数线性简化处理如下:
其中,λi,l表示第i台机组在第l段离散库容的发电效率。
由于涉及变量与的乘积,因此需要使用分离不等式法消除非线性。
设表示第i台机组在第l段离散库容的发电出力,则有:
其中M为数值很大的常量,当时,当时,由于目标函数包含发电量最大部分,因此能够保证约束取到等号。
2.2)处理水库与水库之间、水库与机组之间的空间拓扑联系;
长期水电群运行模拟模型中的水电机组用水与水库用水平衡约束中需要对集合Ωk进行处理。长期水电群运行模拟模型中的水量平衡约束中需要对集合Πk进行处理。
采用矩阵的方式处理水库与机组之间、水库与水库之间的空间约束关系。设水库-机组关联矩阵AR-U表示水库与机组之间的空间拓扑关系,若机组i,j处于水库r上,则关联矢量为:
则水库-机组关联矩阵为,R表示水库,U表示机组:
AR-U=[a1,a2,...,aR]T (19)
设水库-水库关联矩阵AR-R表示水库与水库之间的空间拓扑关系;若水库i是水库r的下游水库,则关联矢量为:
则水库-水库关联矩阵为:
AR-R=[a1,a2,...,aR] (21)
设则关联矢量a* R中非零元素所在位置就表示水库r的上游水库编号。
如上处理后,水库-机组用水平衡约束即可描述为:
其中,矩阵λR-R为元素组成的方阵;为水库R在t时段的用水量,单位m3;为水库R在t时段的弃水量,单位m3;
水量平衡约束即可描述为:
其中,Ik t为水库k在t时段自然来水量,单位m3;矩阵λR-R为元素组成的方阵。
2.3)求解长期水电群运行模拟模型,实现模拟长期跨流域多梯级水电优化运行;
长期水电群运行模拟线性化模型由式(1)、(2)、(4)、(5)、(6)、(7)、(8)、(9)、(12)、(13)、(15)、(17)、(22)和式(23)共同组成,其为大规模混合整数线性规划模型,可以用cplex软件中混合整数规划算法进行求解。
上述各实施例仅用于说明本发明,各个步骤都是可以有所变化的,在本发明技术方案的基础上,凡根据本发明原理对个别步骤进行的改进和等同变换,均不应排除在本发明的保护范围之外。
Claims (9)
1.一种长期跨流域多梯级水电优化运行模拟方法,其特征在于包括以下步骤:
1)构建由目标函数及约束条件组成的长期水电群运行模拟模型;
1.1)建立长期水电优化运行模拟模型的决策变量;
1.2)建立长期水电群运行模拟模型的目标函数;
1.3)建立长期水电群运行模拟模型的约束条件;
2)对长期水电群运行模拟模型进行简化与求解:
2.1)对水电出力转换函数进行非线性处理;
2.2)处理水库与水库之间、水库与机组之间的空间拓扑联系;
2.3)采用已有混合整数规划算法求解长期水电群运行模拟模型,实现模拟长期跨流域多梯级水电优化运行。
2.如权利要求1所述的一种长期跨流域多梯级水电优化运行模拟方法,其特征在于:所述步骤1.1)中,长期水电群运行模拟模型中的决策变量分为五种:Pi t为水电机组i在t时段的出力,单位MW;为水电机组i在t时段的发电用水量,单位m3;为水库k在t时段的库容,单位m3;为水库k在t时段的用水量,单位m3;为水库k在t时段的弃水量,单位m3。
3.如权利要求1所述的一种长期跨流域多梯级水电优化运行模拟方法,其特征在于:所述步骤1.2)中,目标函数为:
其中,E为水电系统在模拟优化运行期内所有水电机组发电量,单位MW;Δt为系统每个计算时段的长度,单位小时;L为目标优化时段数,模型以月为计算时段;N为水电系统的水电机组数。
4.如权利要求1所述的一种长期跨流域多梯级水电优化运行模拟方法,其特征在于:所述步骤1.3)中,具体包括:
1.3.1)建立水电机组运行约束:
1.3.1.1)确定水电出力上下限约束:
Pi,min≤Pi t≤Pi,max
其中,Pi,min为水电机组i最小出力,单位MW;Pi,max为水电机组i最大出力,单位MW;
1.3.1.2)确定水电转换函数:
其中,水电机组i处于水库k中;
1.3.1.3)确定水电系统的水电消纳能力约束:长期水电群运行模拟模型设定系统消纳水电出力的总量不能超过系统负荷预设比例范围:
其中,γ表示水电出力占系统t计算时段全系统负荷的最大比例、λ表示水电出力占系统t计算时段全系统负荷的最小比例;Dt为t时段全系统社会负荷,单位MWh;
1.3.2)建立水库运行约束。
5.如权利要求4所述的一种长期跨流域多梯级水电优化运行模拟方法,其特征在于:所述步骤1.3.2)中,具体包括:
1.3.2.1)水电机组用水与水库用水平衡约束:
其中,Ωk是水库k包含的所有水电机组的集合;
1.3.2.2)水库库容上下限约束:
Vk,min≤Vk≤Vk,max
其中,Vk,min为水库k库容下限,单位m3;Vk为水库k的库容,单位m3;Vk,max为水库k库容上限,单位m3;
1.3.2.3)水库库容始末约束:
其中,Vk,begin为水库k初始计算时段库容边界,单位m3;Vk,end为水库k终止计算时段库容边界,单位m3;
1.3.2.4)水库用水上下限约束:
Qk,min≤Qt k≤Qk,max
其中,Qk,min为水库k用水量下限,单位m3;Qk,max为水库k用水量上限,单位m3;Qk,min主要取决于下游航运、正常供水的限制;Qk,max主要取决于水库的泄流能力;
1.3.2.5)水库弃水量约束:
1.3.2.6)水量平衡约束:
其中,为水库k在t时段自然来水量,单位m3;Πk为水库k包含的所有上游水库的集合;为水库k时段t的蒸发率;λj-k表示水库j用水下泄至水库k的有效比例。
6.如权利要求1所述的一种长期跨流域多梯级水电优化运行模拟方法,其特征在于:所述步骤2.1)中,采用离散库容的方法处理水电转换关系:在水库库容上下限范围内,将库容分成若干段,每段库容范围内用该段的平均库容进行计算;然后将机组出力转换为用水量的一元非线性函数,再进一步将出力-用水曲线线性化。
7.如权利要求6所述的一种长期跨流域多梯级水电优化运行模拟方法,其特征在于,对库容的离散化处理如下:
表示第k个水库,离散为L段,每段的状态变量有且仅有一个状态变量为1,其他均为0:
水库库容离散为L段,第l段的范围为:[Vk,l-1,Vk,l],Vk,l-1为水库库容离散的第l-1段;Vk,l为水库库容离散的第l段;则有:
Vk,0=Vk,min
Vk,L=Vk,max
当库容处于第l段时,即时需满足如下约束:
得到离散库容约束为:
8.如权利要求6或7所述的一种长期跨流域多梯级水电优化运行模拟方法,其特征在于,处于第k个水库的水电机组i的水电转换函数线性简化处理如下:
其中,λi,l表示第i台机组在第l段离散库容的发电效率;
设表示第i台机组在第l段离散库容的发电出力,则有:
其中M为数值很大的常量,当时,当时,
9.如权利要求1所述的一种长期跨流域多梯级水电优化运行模拟方法,其特征在于:所述步骤2.2)中,采用矩阵的方式处理水库与机组之间、水库与水库之间的空间约束关系,设水库-机组关联矩阵AR-U表示水库与机组之间的空间拓扑关系,若机组i,j处于水库r上,则关联矢量为:
则水库-机组关联矩阵为:AR-U=[a1,a2,...,aR]T;
其中R表示水库,U表示机组;设水库-水库关联矩阵AR-R表示水库与水库之间的空间拓扑关系;若水库i是水库r的下游水库,则关联矢量为:
则水库-水库关联矩阵为:AR-R=[a1,a2,...,aR];
设则关联矢量a* R中非零元素所在位置就表示水库r的上游水库编号;则水库-机组用水平衡约束为:
其中,矩阵λR-R为元素组成的方阵;为水库R在t时段的用水量,单位m3;为水库R在t时段的弃水量,单位m3;
水量平衡约束为:
其中,Ik t为水库k在t时段自然来水量,单位m3;矩阵λR-R为元素组成的方阵。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711397656.1A CN108133104B (zh) | 2017-12-21 | 2017-12-21 | 一种长期跨流域多梯级水电优化运行模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711397656.1A CN108133104B (zh) | 2017-12-21 | 2017-12-21 | 一种长期跨流域多梯级水电优化运行模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108133104A true CN108133104A (zh) | 2018-06-08 |
CN108133104B CN108133104B (zh) | 2022-03-08 |
Family
ID=62391308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711397656.1A Active CN108133104B (zh) | 2017-12-21 | 2017-12-21 | 一种长期跨流域多梯级水电优化运行模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108133104B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109117592A (zh) * | 2018-09-14 | 2019-01-01 | 中国长江电力股份有限公司 | 虚拟水库模型及虚拟水库水位计算方法 |
CN109799732A (zh) * | 2019-02-21 | 2019-05-24 | 三峡大学 | 一种数据驱动的梯级水电和新能源联合运行实验系统 |
CN109802445A (zh) * | 2018-12-28 | 2019-05-24 | 国网山西省电力公司电力科学研究院 | 一种考虑输电区段容量限制的多机组协调调度方法 |
CN109936164A (zh) * | 2019-03-31 | 2019-06-25 | 东北电力大学 | 基于电源互补特性分析的多能源电力系统优化运行方法 |
CN111431173A (zh) * | 2020-03-30 | 2020-07-17 | 国网经济技术研究院有限公司 | 考虑负荷互补特性的中压配电网环网结构确定方法及系统 |
CN111950982A (zh) * | 2020-08-10 | 2020-11-17 | 中国电力科学研究院有限公司 | 梯级水电水平衡模型构建方法、系统、设备及存储介质 |
CN112907078A (zh) * | 2021-02-20 | 2021-06-04 | 中国电力科学研究院有限公司 | 水电耦合的市场出清方法、系统、设备及可读存储介质 |
CN113780628A (zh) * | 2021-08-16 | 2021-12-10 | 西安交通大学 | 梯级水电调度模型构建方法、装置、电子设备及存储介质 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101705671A (zh) * | 2009-11-19 | 2010-05-12 | 武汉大学 | 黄河上游梯级水电站运行设计与优化调度方法及设备 |
CN102682409A (zh) * | 2012-05-10 | 2012-09-19 | 中国水利水电科学研究院 | 一种基于gams非线性规划梯级水库群优化调度方法 |
CN102817335A (zh) * | 2012-07-12 | 2012-12-12 | 贵州乌江水电开发有限责任公司 | 一种梯级水库群联合防洪优化调度的方法及系统 |
CN102855591A (zh) * | 2012-08-14 | 2013-01-02 | 贵州乌江水电开发有限责任公司 | 梯级水库群短期联合发电优化调度方法及系统 |
CN102867275A (zh) * | 2012-08-14 | 2013-01-09 | 贵州乌江水电开发有限责任公司 | 梯级水库群中长期联合发电优化调度方法及系统 |
CN104123589A (zh) * | 2014-06-24 | 2014-10-29 | 华中科技大学 | 一种梯级水电站短期优化调度方法 |
CN105869070A (zh) * | 2016-04-06 | 2016-08-17 | 大连理工大学 | 一种跨流域梯级水电站群效益均衡的合作优化调度方法 |
CN106682810A (zh) * | 2016-11-17 | 2017-05-17 | 三峡大学 | 巨型水电站动态投产下跨流域梯级水电站群长期运行方法 |
US9811061B1 (en) * | 2001-05-18 | 2017-11-07 | The Energy Authority, Inc. | Method for management and optimization of hydropower generation and consumption |
-
2017
- 2017-12-21 CN CN201711397656.1A patent/CN108133104B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9811061B1 (en) * | 2001-05-18 | 2017-11-07 | The Energy Authority, Inc. | Method for management and optimization of hydropower generation and consumption |
CN101705671A (zh) * | 2009-11-19 | 2010-05-12 | 武汉大学 | 黄河上游梯级水电站运行设计与优化调度方法及设备 |
CN102682409A (zh) * | 2012-05-10 | 2012-09-19 | 中国水利水电科学研究院 | 一种基于gams非线性规划梯级水库群优化调度方法 |
CN102817335A (zh) * | 2012-07-12 | 2012-12-12 | 贵州乌江水电开发有限责任公司 | 一种梯级水库群联合防洪优化调度的方法及系统 |
CN102855591A (zh) * | 2012-08-14 | 2013-01-02 | 贵州乌江水电开发有限责任公司 | 梯级水库群短期联合发电优化调度方法及系统 |
CN102867275A (zh) * | 2012-08-14 | 2013-01-09 | 贵州乌江水电开发有限责任公司 | 梯级水库群中长期联合发电优化调度方法及系统 |
CN104123589A (zh) * | 2014-06-24 | 2014-10-29 | 华中科技大学 | 一种梯级水电站短期优化调度方法 |
CN105869070A (zh) * | 2016-04-06 | 2016-08-17 | 大连理工大学 | 一种跨流域梯级水电站群效益均衡的合作优化调度方法 |
CN106682810A (zh) * | 2016-11-17 | 2017-05-17 | 三峡大学 | 巨型水电站动态投产下跨流域梯级水电站群长期运行方法 |
Non-Patent Citations (3)
Title |
---|
杜尔顺 等: "适应低碳电源发展的低碳电网规划模型", 《电网技术》 * |
杨柳 等: "南方电网节能与经济运行评价方法", 《电力系统自动化》 * |
贾江涛 等: "考虑水头影响的梯级水电站群短期优化调度", 《电力系统自动化》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109117592A (zh) * | 2018-09-14 | 2019-01-01 | 中国长江电力股份有限公司 | 虚拟水库模型及虚拟水库水位计算方法 |
CN109802445A (zh) * | 2018-12-28 | 2019-05-24 | 国网山西省电力公司电力科学研究院 | 一种考虑输电区段容量限制的多机组协调调度方法 |
CN109799732A (zh) * | 2019-02-21 | 2019-05-24 | 三峡大学 | 一种数据驱动的梯级水电和新能源联合运行实验系统 |
CN109936164A (zh) * | 2019-03-31 | 2019-06-25 | 东北电力大学 | 基于电源互补特性分析的多能源电力系统优化运行方法 |
CN109936164B (zh) * | 2019-03-31 | 2022-02-01 | 东北电力大学 | 基于电源互补特性分析的多能源电力系统优化运行方法 |
CN111431173A (zh) * | 2020-03-30 | 2020-07-17 | 国网经济技术研究院有限公司 | 考虑负荷互补特性的中压配电网环网结构确定方法及系统 |
CN111431173B (zh) * | 2020-03-30 | 2021-08-17 | 国网经济技术研究院有限公司 | 考虑负荷互补特性的中压配电网环网结构确定方法及系统 |
CN111950982A (zh) * | 2020-08-10 | 2020-11-17 | 中国电力科学研究院有限公司 | 梯级水电水平衡模型构建方法、系统、设备及存储介质 |
CN111950982B (zh) * | 2020-08-10 | 2024-04-19 | 中国电力科学研究院有限公司 | 梯级水电水平衡模型构建方法、系统、设备及存储介质 |
CN112907078A (zh) * | 2021-02-20 | 2021-06-04 | 中国电力科学研究院有限公司 | 水电耦合的市场出清方法、系统、设备及可读存储介质 |
CN113780628A (zh) * | 2021-08-16 | 2021-12-10 | 西安交通大学 | 梯级水电调度模型构建方法、装置、电子设备及存储介质 |
CN113780628B (zh) * | 2021-08-16 | 2024-04-12 | 西安交通大学 | 梯级水电调度模型构建方法、装置、电子设备及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN108133104B (zh) | 2022-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108133104B (zh) | 一种长期跨流域多梯级水电优化运行模拟方法 | |
CN107153885B (zh) | 考虑火电机组深度调峰的实时发电计划优化方法 | |
CN109242350B (zh) | 计及可平移负荷的冷热电联供系统容量优化配置方法 | |
CN107276122B (zh) | 适应大规模可再生能源并网的调峰资源调用决策方法 | |
CN108599269B (zh) | 一种考虑风险成本的大电网自适应鲁棒备用优化方法 | |
WO2019006733A1 (zh) | 一种跨省互联水电站群长期联合调峰调度方法 | |
CN102182634B (zh) | 基于改进粒子群的海岛风柴蓄发电功率优化设计方法 | |
CN104795846A (zh) | 一种抽水蓄能电站与风电的联合系统的优化运行方法 | |
CN116667395B (zh) | 基于梯级水电改造的水风光蓄互补泵站容量配置方法 | |
CN110994606B (zh) | 一种基于复杂适应系统理论的多能源电源容量配置方法 | |
CN107104467A (zh) | 一种考虑核电调峰及其安全约束的机组组合优化方法 | |
CN110350512A (zh) | 一种智能园区新能源发电站调度优化方法及系统 | |
CN103795088B (zh) | 一种基于负荷曲线量化的抽水蓄能电站优化调度方法 | |
CN110690729B (zh) | 一种调节电力系统峰谷差的抽蓄优化调度方法 | |
CN104037814B (zh) | 一种分段平滑法与邻域搜索算法相结合的多电网调峰负荷分配方法 | |
CN115081700A (zh) | 基于综合储能技术的数据中心多能协同优化方法及系统 | |
CN114362255A (zh) | 一种源网荷储电力系统多目标日前调度优化方法及系统 | |
CN110571868B (zh) | 微电网的优化配置方法 | |
Li et al. | Research on optimizing operation of the single reservoir of hybrid pumped storage power station | |
CN110649592A (zh) | 含蓄热式电锅炉的电热互联能源系统多场景经济调度方法 | |
Ma et al. | Two-stage optimal dispatching based on wind-photovoltaic-pumped storage-thermal power combined power generation system | |
CN105244879A (zh) | 考虑核电机组检修的核电抽水储能电站联合运行优化方法 | |
Shao et al. | An MILP Model for Annual Power Balance Based on Monthly Multiple Typical Daily Load Scenarios | |
Gumilar et al. | Load Factor Improvement on Daily Load Curve Using Pump Storage to SavingProduction Cost | |
Xiao et al. | Complementary coordinated dispatch of multi-energy systems considering power flexibility margin |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |