CN106682810A - 巨型水电站动态投产下跨流域梯级水电站群长期运行方法 - Google Patents
巨型水电站动态投产下跨流域梯级水电站群长期运行方法 Download PDFInfo
- Publication number
- CN106682810A CN106682810A CN201611022560.2A CN201611022560A CN106682810A CN 106682810 A CN106682810 A CN 106682810A CN 201611022560 A CN201611022560 A CN 201611022560A CN 106682810 A CN106682810 A CN 106682810A
- Authority
- CN
- China
- Prior art keywords
- power station
- basin
- water
- year
- load
- 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
- 238000000034 method Methods 0.000 title claims abstract description 78
- 230000007774 longterm Effects 0.000 title claims abstract description 9
- 238000009434 installation Methods 0.000 claims abstract description 18
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 99
- 230000005611 electricity Effects 0.000 claims description 38
- 230000008569 process Effects 0.000 claims description 31
- 239000008239 natural water Substances 0.000 claims description 27
- 230000005619 thermoelectricity Effects 0.000 claims description 12
- 238000011161 development Methods 0.000 claims description 9
- 230000037361 pathway Effects 0.000 claims description 9
- 238000009826 distribution Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 claims description 5
- 230000002411 adverse Effects 0.000 claims description 4
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000001447 compensatory effect Effects 0.000 claims description 4
- PEDCQBHIVMGVHV-UHFFFAOYSA-N Glycerine Chemical compound OCC(O)CO PEDCQBHIVMGVHV-UHFFFAOYSA-N 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000000205 computational method Methods 0.000 claims description 3
- 238000007796 conventional method Methods 0.000 claims description 3
- 238000005265 energy consumption Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 229910052729 chemical element Inorganic materials 0.000 claims 1
- 230000008901 benefit Effects 0.000 abstract description 4
- 230000008859 change Effects 0.000 abstract description 3
- 230000003442 weekly effect Effects 0.000 abstract 1
- 241000196324 Embryophyta Species 0.000 description 5
- 238000005457 optimization Methods 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000011144 upstream manufacturing Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 240000002853 Nelumbo nucifera Species 0.000 description 1
- 235000006508 Nelumbo nucifera Nutrition 0.000 description 1
- 235000006510 Nelumbo pentapetala Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- 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
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Quality & Reliability (AREA)
- Educational Administration (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Supply And Distribution Of Alternating Current (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及巨型水电站动态投产下跨流域梯级水电站群长期运行方法,它包括以下步骤:步骤1,预测负荷。步骤2,从月平均负荷中扣除不参与计算的电站出力过程。步骤3,计算全网流域不同频率天然来水过程。步骤4,计算各水电站出力过程;步骤5,计算火电站出力过程。该方法结合云南电网实际运行特点,将水电站分为两类:参与计算和不参与计算,参与计算的电站主要包括调节性能较好、装机规模很大的水电站和火电站,不参与计算的电站主要包括风电、光伏、周调节性能以下水电和部分火电,这部分电站长期运行过程基本变化不大,可按照电网的实际运行情况给定其出力过程。
Description
技术领域
本发明涉及水电站工程技术领域,尤其是涉及一种巨型水电站动态投产下跨流域梯级水电站群长期运行方法。
背景技术
在今后20年间,随着金沙江、雅砻江、大渡河、澜沧江、怒江、黄河上游干流梯级上大容量、巨型机组的水电站群大规模动态投产,水电集中化程度将更高。当水电比例增长到70%以上后,极端气候下的天然来水对电网、水电系统的安全运行影响越来越大,包括大规模水电输送、汛期水电大规模消纳、枯期跨区跨省的电网调峰等复杂问题,最近几年云南连续干旱对电网停电、调度计划安排造成重大影响就是一个突出例子。因此,在巨型水电站群动态投产情况下,需要开展不同频率天然来水对电网及水电调度方式影响的研究,提出好的调度预案及对策,防范可能发生的风险。
巨型水电站群动态投产面临的问题属于大规模、不确定性的水电系统优化问题,涉及不确定来水、负荷预测条件下,水电系统优化(包括梯级电站群调度、跨流域水电站群调度等)、多电源协调(特别是水火电联合优化)等关键问题,目前国内外学者已经针对上述部分问题开展研究,刘双全等针对径流不确定性梯级水电站群优化调度问题,以澜沧江梯级水电站群为研究背景,采用考虑径流随机特性的随机优化调度策略来指导澜沧江梯级水电站群中长期运行,但是,该技术方案仍然存在不足,只能适应梯级水电站个数是固定情况,不适应动态投产的跨流域梯级水电站群。原因之一是水电站动态投产对上下游流量演进影响甚大,特别是在两座已建成电站之间重新投产一座或多座调节能力较强水电站,完全打破了原有的上下游流量传播规律,利用历史径流模拟出来的随机径流不能反映水电站动态投产之后梯级上下游水力联系;原因之二是单一梯级水电站群在地理位置比较接近,其水文条件(如降雨、径流等)也相似,而跨流域梯级水电站由于地理位置相距较远,导致它们之间的水文特性差异较大,相同年份的天然来水频率却不相同,难以获得接近历史规律的径流过程。曾筠等针对南方电网跨流域水电系统,提出跨流域水电站群长期优化调度的可变策略搜索求解算法,该方法可以快速求解跨流域水电系统长期优化调度问题,实现跨流域水电站群补偿调度,利于更好地配置水电系统资源,发挥跨流域水电站群的综合效益。但是,该技术方案仍然存在不足,只能适应固定个数的跨流域梯级水电站群,不适应动态投产的梯级水电站群,并且该技术方案仅涉及水电,没有涉及火电等其它能源,不能分析整个电力系统电力电量平衡情况,所以,当大规模跨流域梯级水电站群处于动态投产情况下时,该技术方案不能很好地配置水电系统资源,发挥跨流域水电站群的综合效益。
随着我国西南地区水电规模不断扩大,水电在电网中的装机容量和发电量达到绝对支配地位,研究一种巨型水电站动态投产下跨流域梯级水电站群长期运行方法,对保障水电占绝对支配地位的电力系统安全、稳定、经济运行至关重要。
发明内容
本发明所要解决的技术问题是提供一种巨型水电站动态投产下跨流域梯级水电站群长期运行方法。该方法结合云南电网实际运行特点,将水电站分为两类:参与计算和不参与计算,参与计算的电站主要包括调节性能较好、装机规模很大的水电站和火电站,不参与计算的电站主要包括风电、光伏、周调节性能以下水电和部分火电,这部分电站长期运行过程基本变化不大,可按照电网的实际运行情况给定其出力过程。
为了解决上述技术问题,本发明提出以下技术方案:巨型水电站动态投产下跨流域梯级水电站群长期运行方法,它包括以下步骤:
步骤1,预测负荷,根据云南省规划及经济发展情况,结合历年月平均负荷预测出2012年~2015年的月平均负荷;
步骤2,从月平均负荷中扣除不参与计算的电站出力过程,这部分电站包括风电、光伏、周调节性能以下水电站和部分火电站,从而得到参与计算的水电站和火电站共同承担负荷Ct,t为月时段;
步骤3,计算全网流域不同频率天然来水过程;
步骤4,计算各水电站出力过程。采用最大剩余负荷最小模型对参与计算的水电站进行优化调度,目标函数为:
式中:目标函数F表示电网负荷经水电调节后的剩余负荷最大值最小;Hm,t表示水电站m在时段t时的平均出力;M表示参与计算的水电站总数,T为计算周期。约束条件主要包括:水量平衡方程、发电流量上限、出力流量上限、水库水位上下限、出力上下限、水库始末水位限制等;本发明中提到的水电站m和水库m均指同一对象。
步骤5,计算火电站出力过程。各火电站考虑出力上下限后,按其装机容量比平衡剩余负荷。
所述步骤1中预测负荷方法为:
根据云南省规划及经济发展情况,结合2007-2011年GDP、电力、电量增长率,预测2012年~2015年月平均负荷。根据“十二五”规划,云南省经济发展高方案为:2010~2020年GDP年均增长9.3%;中方案2010~2020年GDP年均增长8.4%;低方案2010~2020年GDP年均增长7.5%;实际2007-2011年电量平均增长率13.4%,最大负荷增长率13.2%,同期GDP年均增长约12%。因此,按单位GDP能耗不变,预测2012年~2015年高方案负荷和电量平均增长率10.5%和10.7%;中方案负荷和电量平均增长率9.6%和9.8%;低方案负荷和电量平均增长率8.7%和8.8%;另一方面,为体现“十二五”期间国家节能降耗水平提高的要求,采用逐年近似等值降低增长率方案来预测2012年~2015年负荷。
所述步骤3中天然来水丰枯评估方法为:
单一省级电站包含的流域往往多达几十条,这些流域受不同的地理气候及大气环流影响,其水文特性并不同步,很难从全网的角度确定各条流域相同频率的天然来水过程。为解决此问题,本发明以年为统计时段,计算各流域不同年份来水量与多年平均来水量的差距——距平值,然后以各大流域的总装机容量为权重,结合各流域历年天然入库水量的距平值,求出全网全流域在不同装机规模下的距平值,最后对距平值进行频率排序,通过线性插值取25%、50%和75%对应的年份分别表示全网全流域的丰水年、平水年和枯水年;
单一梯级水库群天然来水的丰枯年计算方法是将梯级所有水库纳入基本计算单元,统计年内所有水库的天然入库流量与天然来水多年平均值,通过下面公式进行计算:
式中:Dj,k表示第j条流域在第k年天然来水距平值,%;Qm,k和分别表示水库m在k年的平均入库流量和多年平均入库流量,m3/s;Mj表示第j条流域的水库数量;
以上为同一流域的天然来水丰枯评估方法,以此为基础分析全网各大流域水库综合来水丰枯评估方法如下式:
式中:Ak表示全网第k年所有流域综合来水距平值,%;Cj和Gj分别表示第j条流域来水距平权重系数和总装机容量,MW;J表示全网流域总数。
所述步骤4中最大剩余负荷最小模型求解方法为:
最大剩余负荷最小模型是给定调度期内入库流量过程和水库始末水位,在考虑各种约束条件下,确定各个具有长期调节能力的水电站水库的调度过程,使得系统最大剩余平均负荷最小。该模型主要应用于水火电补偿调度的情况下,目的在于加大水电在枯期发电量,增大水电枯期调峰能力,从而减小水电站水库天然径流时空分配不均的不利影响,使留给火电和其他电源的剩余负荷在调度期内尽量均匀,使系统对火电装机需求尽量小。
所述步骤5中计算火电站出力过程为:
本发明利用POA(Progress Optimality Algorithem)算法对最大剩余负荷最小模型进行求解,其中水电站出力过程求解步骤如下:
步骤1:按常规方法得到初始解。
步骤2:将所有参与计算的水电站按流域个数划分J组,即全网流域总数为J,每条流域包含的水电站个数为Mj,j表示流域编号,j≤J。
步骤3:设起始计算时段t=1,计算周期为T,初始搜索步长ε/2。
步骤4:设j=1;
步骤5:将第j组水电站m在时段t末水位和t+2时段末水位固定,用以下步骤求解水电站m在时段t的出库流量
a)若t<T-1,判断前一轮t+1时段各电站末水位是否有变化,同时若t>0,判断本轮次t-1时段末水位是否有变化,若都无变化则t=t+1返回步骤3。
b)每组内按照设定初始步长ε,在水电站m初始出库流量上下得到初始廊道,由sm代表廊道中的三个水位值,即这样,本组内所有的组合在这一时段具有个状态。采用动态规划的方法求解这一个两时段的问题,求得最优的状态使得目标函数最大。
步骤6:设j=j+1,重复步5直到j>J,转入步骤7。
步骤7:如果采用单步收敛模式,但单步没有收敛,则转入步骤4,否则转入步骤8。
步骤8:如果采用单步收敛模式,设ε=ε/2,若未达精度要求,返回步骤5,否则转入步骤9。
步骤9:设t=t+1,若t<T-1,则转入步骤4,否则转入步骤10。
步骤10:若未收敛,转入步骤3;否则若采用单步收敛模式,设ε=ε/2,若未达到精度要求,转入步骤3,否则转入步骤11。
步骤11:寻优结束,输出结果。
本发明有如下有益效果:
(1)本发明通过全流域距平值评估出水文特性不同步的全流域相同频率的天然来水过程,以各流域总装机容量为权重考虑了水电站动态投产情况,进而解决了巨型水电站动态投产下跨流域梯级水电站群长期运行时来水频率不同步的问题,提出的全网全流域丰枯来水评价方法在一定程度上可以反映出不同频率的天然来水对云南水电运行方式的影响。
(2)本发明提出的最大剩余负荷最小模型可以加大水电在枯期发电量,增大水电枯期调峰能力,从而减小水电站水库天然径流时空分配不均的不利影响,使留给火电和其他电源的剩余负荷在调度期内尽量均匀,使系统对火电装机需求尽量小,从而达到节能减排的作用。
附图说明
下面结合附图和实施例对本发明作进一步说明。
图1是云南省大中型水电站分布图。
图2是云南电网装机容量分布图。
图3是云南电网全部流域2012年~2015年丰枯评估结果。
图4-1是2012-2015年云南电网在丰水年月平均负荷平衡图。
图4-2是2012-2015年云南电网在枯水年月平均负荷平衡图。
图4-3是2012-2015年云南电网在平水年月平均负荷平衡图。
具体实施方式
下面结合附图对本发明的实施方式做进一步的说明。
实施例1:
本发明结合云南电网实际运行特点,将云南电网统调电站分为两类:参与计算和不参与计算,参与计算的电站主要包括调节性能较好、装机规模很大的水电站和火电站,不参与计算的电站主要包括风电、光伏、周调节性能以下水电和部分火电,这部分电站其长期运行过程基本变化不大,可按照电网的实际运行情况给定其出力过程。求解巨型水电站动态投产下跨流域梯级水电站群长期运行方法问题应按照下述步骤(1)-(5)予以实现:
(1)预测负荷。根据云南省规划及经济发展情况,结合历年月平均负荷预测出2012年~2015年的月平均负荷。
根据云南省规划及经济发展情况,结合2007-2011年GDP、电力、电量增长率,预测2012年~2015年月平均负荷。根据“十二五”规划,云南省经济发展高方案为:2010~2020年GDP年均增长9.3%;中方案2010~2020年GDP年均增长8.4%;低方案2010~2020年GDP年均增长7.5%;实际2007-2011年电量平均增长率13.4%,最大负荷增长率13.2%,同期GDP年均增长约12%。因此,按单位GDP能耗不变,预测2012年~2015年高方案负荷和电量平均增长率10.5%和10.7%;中方案负荷和电量平均增长率9.6%和9.8%;低方案负荷和电量平均增长率8.7%和8.8%;另一方面,为体现“十二五”期间国家节能降耗水平提高的要求,采用逐年近似等值降低增长率方案来预测2012年~2015年负荷。因此,预测得出的高、中、低负荷方案如表1:
表1:2012-2015年各方案电力电量预测增长率
(2)从月平均负荷中扣除不参与计算的电站出力过程,这部分电站包括风电、光伏、周调节性能以下水电站和部分火电站,从而得到参与计算的水电站和火电站共同承担负荷Ct,t为月时段。
(3)计算全网流域不同频率天然来水过程。
单一省级电站包含的流域往往多达几十条,这些流域受不同的地理气候及大气环流影响,其水文特性并不同步,很难从全网的角度确定各条流域相同频率的天然来水过程。为解决此问题,本发明以年为统计时段,计算各流域不同年份来水量与多年平均来水量的差距——距平值,然后以各大流域的总装机容量为权重,结合各流域历年天然入库水量的距平值,求出全网全流域在不同装机规模下的距平值,最后对距平值进行频率排序,通过线性插值取25%、50%和75%对应的年份分别表示全网全流域的丰水年、平水年和枯水年。
单一梯级水库群天然来水的丰枯年计算方法是将梯级所有水库纳入基本计算单元,统计年内所有水库的天然入库流量与天然来水多年平均值,通过下面公式进行计算:
式中:Dj,k表示第j条流域在第k年天然来水距平值,%;Qm,k和分别表示水库m在k年的平均入库流量和多年平均入库流量,m3/s;Mj表示第j条流域的水库数量;
以上为同一流域的天然来水丰枯评估方法,以此为基础分析全网各大流域水库综合来水丰枯评估方法如下式:
式中:Ak表示全网第k年所有流域综合来水距平值,%;Cj和Gj分别表示第j条流域来水距平权重系数和总装机容量,MW;J表示全网流域总数。
(4)计算各水电站出力过程。采用最大剩余负荷最小模型对参与计算的水电站进行优化调度,目标函数为:
式中:目标函数F表示电网负荷经水电调节后的剩余负荷最大值最小;Hm,t表示水电站m在时段t时的平均出力;M表示参与计算的水电站总数,T为计算周期。约束条件主要包括:水量平衡方程、发电流量上限、出力流量上限、水库水位上下限、出力上下限、水库始末水位限制等;本发明中提到的水电站m和水库m均指同一对象。
最大剩余负荷最小模型是给定调度期内入库流量过程和水库始末水位,在考虑各种约束条件下,确定各个具有长期调节能力的水电站水库的调度过程,使得系统最大剩余平均负荷最小。该模型主要应用于水火电补偿调度的情况下,目的在于加大水电在枯期发电量,增大水电枯期调峰能力,从而减小水电站水库天然径流时空分配不均的不利影响,使留给火电和其他电源的剩余负荷在调度期内尽量均匀,使系统对火电装机需求尽量小。
(5)计算火电站出力过程。各火电站考虑出力上下限后,按其装机容量比平衡剩余负荷。
本发明利用POA(Progress Optimality Algorithem)算法对最大剩余负荷最小模型进行求解,其中水电站出力过程求解步骤如下:
步骤1:按常规方法得到初始解。
步骤2:将所有参与计算的水电站按流域个数划分J组,即全网流域总数为J,每条流域包含的水电站个数为Mj,j表示流域编号,j≤J。
步骤3:设起始计算时段t=1,计算周期为T,初始搜索步长ε/2。
步骤4:设j=1;
步骤5:将第j组水电站m在时段t末水位和t+2时段末水位固定,用以下步骤求解水电站m在时段t的出库流量
a)若t<T-1,判断前一轮t+1时段各电站末水位是否有变化,同时若t>0,判断本轮次t-1时段末水位是否有变化,若都无变化则t=t+1返回步骤3。
b)每组内按照设定初始步长ε,在水电站m初始出库流量上下得到初始廊道,由sm代表廊道中的三个水位值,即这样,本组内所有的组合在这一时段具有个状态。采用动态规划的方法求解这一个两时段的问题,求得最优的状态使得目标函数最大。
步骤6:设j=j+1,重复步5直到j>J,转入步骤7。
步骤7:如果采用单步收敛模式,但单步没有收敛,则转入步骤4,否则转入步骤8。
步骤8:如果采用单步收敛模式,设ε=ε/2,若未达精度要求,返回步骤5,否则转入步骤9。
步骤9:设t=t+1,若t<T-1,则转入步骤4,否则转入步骤10。
步骤10:若未收敛,转入步骤3;否则若采用单步收敛模式,设ε=ε/2,若未达到精度要求,转入步骤3,否则转入步骤11。
步骤11:寻优结束,输出结果。
以云南电网为研究背景,利用预测得到2012-2015年月平均负荷,模拟巨型水电站动态投产下跨流域梯级水电站群长期运行,并以2011年为基准年对模拟结果进行对比分析,以验证本发明提出的求解方法。截至2011年底,云南省总装机容量为40817MW,中调及以上调度装机容量32723MW,其中统调火电站为13座,装机容量为10600MW,占总装机的32.39%;统调水电站为107座,装机容量为21313MW,占总装机容量65.14%;风电789MW,占总装机的2.41%;光伏20MW,占总装机的0.06%,可以看出,云南省调平衡水电站数量非常庞大,但网内具有季调节及以上水库调节能力的水电站只有13座,仅占12%,其余88%水电站为径流式调节水库或调节性能差的日、周调节水库,发电能力受径流影响较大,丰水期无法充分利用水资源,汛期弃水调峰现象时有发生。根据“十二五”规划,国家将大力开发金沙江和澜沧江流域,如图1所示,届时云南水电2012-2015年主要集中分布在金沙江和澜沧江两大干流梯级水电站群,图2表示金沙江和澜沧江干流梯级水电总装机容量比重图;表2表示金沙江和澜沧江流域梯级电站2012-2015年动态投产情况:
表2:金沙江和澜沧江流域梯级动态投产情况:
采用步骤3的详细步骤对云南电网全部流域2012-2015年丰枯评估结果如图3,可以看出,由于不同年份各流域投产装机容量不同,计算出来的权重及全网天然来水距平值出不同,将各年份的距平值从大到小排列,分别取频率为25%、50%和75%对应的年份为丰水年、平水年和枯水年,因此,得到丰平枯水年份也不同。
采用整个求解思路对云南电网2012-2015年水火电进行电力电量平衡计算,其水电计算结果如表3和图4,其中列(3)、(6)和(12)表示该年年发电量;(4)和(7)表示该年枯期(1-5,11,12月)发电量;(5)、(8)和(12)表示枯期电量/汛期电量;(11)表示超出负荷的弃电量。
表3:澜沧江、金沙江及全网水电模拟调度结果 单位:10亿kWh
注:列(3)、(6)和(12)表示该年总发电量;(4)和(7)表示该年枯期(1-5,11,12月)发电量;(5)、(8)和(12)表示枯期电量/汛期电量;(11)表示超出负荷的弃电量;
从表3可以看出:
(1)从第(5)、(8)和(12)列数据可以看出,2012~2015年与2011年相比,澜沧江、金沙江以及全网水电的枯丰电量比均有所提高,如全网水电2012-2014年丰平枯水年平均枯丰电量比分别为0.91、0.85和0.82,均比2011年的0.8有所提高,说明随着巨型梯级电站的动态投产,极大地改善了云南水电系统的发电特性,缓解了丰枯发电不平衡特性,这种改变在丰水年更加突出;但全网枯丰电网量从2012-2015年呈下降趋势,主要原因在于后期投产的水电站调节性能较差。
(2)从第(11)列数据可以看出,说明随着巨型梯级水电站的投产,基本消除了枯期缺电现象,但由于火电有最小出力限制,全网水电出现了电量过剩的现象,因此可以适当加大外送消纳这部分电量;
(3)从第(3)、(6)和(9)列数据可以看出,丰水年、平水年和枯水年发电量呈依次递减趋势,但2013年全网水电优化调度结果中,枯水年的发电量却比丰水年、平水年的大,原因在于金沙江流域来水非常大,而且大多数年份的天然来水集中在6~10月(约占全年的75%以上),但金沙江流域的梯级水电站均为季调节电站,调节能力有限,丰水年的丰水期弃水量较大,造成丰水年的发电量比平水年和枯水年少。
随着中国西南地区水电在电网中的比重迅速提升,巨型水电站动态投产对电网水电系统的运行方式将产生重大影响,包括极端干旱气候对电网安全运行方式的影响,以及水电外送电规模和其它电源规模等系统问题。本发明以云南电网为背景,采取一系列简化问题的策略,重点分析金沙江和澜沧江流域2012年~2015年巨型梯级水电站动态投产下的云南水电系统模拟调度,评估出了不同频率天然来水下云南水电的发电能力,以便调度部门做出科学决策,防范可能发生的弃电与缺电等风险。
通过上述的说明内容,本领域技术人员完全可以在不偏离本项发明技术思想的范围内,进行多样的变更以及修改都在本发明的保护范围之内。本发明的未尽事宜,属于本领域技术人员的公知常识。
Claims (5)
1.巨型水电站动态投产下跨流域梯级水电站群长期运行方法,其特征在于,它包括以下步骤:
步骤1,预测负荷,根据云南省规划及经济发展情况,结合历年月平均负荷预测出2012年~2015年的月平均负荷;
步骤2,从月平均负荷中扣除不参与计算的电站出力过程,这部分电站包括风电、光伏、周调节性能以下水电站和部分火电站,从而得到参与计算的水电站和火电站共同承担负荷Ct,t为月时段;
步骤3,计算全网流域不同频率天然来水过程;
步骤4,计算各水电站出力过程。采用最大剩余负荷最小模型对参与计算的水电站进行优化调度,目标函数为:
式中:目标函数F表示电网负荷经水电调节后的剩余负荷最大值最小;Hm,t表示水电站m在时段t时的平均出力;M表示参与计算的水电站总数,T为计算周期。约束条件主要包括:水量平衡方程、发电流量上限、出力流量上限、水库水位上下限、出力上下限、水库始末水位限制等;本发明中提到的水电站m和水库m均指同一对象。
步骤5,计算火电站出力过程。各火电站考虑出力上下限后,按其装机容量比平衡剩余负荷。
2.根据权利要求1所述的巨型水电站动态投产下跨流域梯级水电站群长期运行方法,其特征在于,所述步骤1中预测负荷方法为:
根据云南省规划及经济发展情况,结合2007-2011年GDP、电力、电量增长率,预测2012年~2015年月平均负荷。根据“十二五”规划,云南省经济发展高方案为:2010~2020年GDP年均增长9.3%;中方案2010~2020年GDP年均增长8.4%;低方案2010~2020年GDP年均增长7.5%;实际2007-2011年电量平均增长率13.4%,最大负荷增长率13.2%,同期GDP年均增长约12%。因此,按单位GDP能耗不变,预测2012年~2015年高方案负荷和电量平均增长率10.5%和10.7%;中方案负荷和电量平均增长率9.6%和9.8%;低方案负荷和电量平均增长率8.7%和8.8%;另一方面,为体现“十二五”期间国家节能降耗水平提高的要求,采用逐年近似等值降低增长率方案来预测2012年~2015年负荷。
3.根据权利要求1所述的巨型水电站动态投产下跨流域梯级水电站群长期运行方法,其特征在于,所述步骤3中天然来水丰枯评估方法为:
单一省级电站包含的流域往往多达几十条,这些流域受不同的地理气候及大气环流影响,其水文特性并不同步,很难从全网的角度确定各条流域相同频率的天然来水过程。为解决此问题,本发明以年为统计时段,计算各流域不同年份来水量与多年平均来水量的差距——距平值,然后以各大流域的总装机容量为权重,结合各流域历年天然入库水量的距平值,求出全网全流域在不同装机规模下的距平值,最后对距平值进行频率排序,通过线性插值取25%、50%和75%对应的年份分别表示全网全流域的丰水年、平水年和枯水年;
单一梯级水库群天然来水的丰枯年计算方法是将梯级所有水库纳入基本计算单元,统计年内所有水库的天然入库流量与天然来水多年平均值,通过下面公式进行计算:
式中:Dj,k表示第j条流域在第k年天然来水距平值,%;Qm,k和分别表示水库m在k年的平均入库流量和多年平均入库流量,m3/s;Mj表示第j条流域的水库数量;
以上为同一流域的天然来水丰枯评估方法,以此为基础分析全网各大流域水库综合来水丰枯评估方法如下式:
式中:Ak表示全网第k年所有流域综合来水距平值,%;Cj和Gj分别表示第j条流域来水距平权重系数和总装机容量,MW;J表示全网流域总数。
4.根据权利要求1所述的巨型水电站动态投产下跨流域梯级水电站群长期运行方法,其特征在于,所述步骤4中最大剩余负荷最小模型求解方法为:
最大剩余负荷最小模型是给定调度期内入库流量过程和水库始末水位,在考虑各种约束条件下,确定各个具有长期调节能力的水电站水库的调度过程,使得系统最大剩余平均负荷最小。该模型主要应用于水火电补偿调度的情况下,目的在于加大水电在枯期发电量,增大水电枯期调峰能力,从而减小水电站水库天然径流时空分配不均的不利影响,使留给火电和其他电源的剩余负荷在调度期内尽量均匀,使系统对火电装机需求尽量小。
5.根据权利要求1所述的巨型水电站动态投产下跨流域梯级水电站群长期运行方法,其特征在于,所述步骤5中计算火电站出力过程为:
本发明利用POA(Progress Optimality Algorithem)算法对最大剩余负荷最小模型进行求解,其中水电站出力过程求解步骤如下:
步骤1:按常规方法得到初始解。
步骤2:将所有参与计算的水电站按流域个数划分J组,即全网流域总数为J,每条流域包含的水电站个数为Mj,j表示流域编号,j≤J。
步骤3:设起始计算时段t=1,计算周期为T,初始搜索步长ε/2。
步骤4:设j=1;
步骤5:将第j组水电站m在时段t末水位和t+2时段末水位固定,用以下步骤求解水电站m在时段t的出库流量
a)若t<T-1,判断前一轮t+1时段各电站末水位是否有变化,同时若t>0,判断本轮次t-1时段末水位是否有变化,若都无变化则t=t+1返回步骤3。
b)每组内按照设定初始步长ε,在水电站m初始出库流量上下得到初始廊道,由sm代表廊道中的三个水位值,即这样,本组内所有的组合在这一时段具有个状态。采用动态规划的方法求解这一个两时段的问题,求得最优的状态使得目标函数最大。
步骤6:设j=j+1,重复步5直到j>J,转入步骤7。
步骤7:如果采用单步收敛模式,但单步没有收敛,则转入步骤4,否则转入步骤8。
步骤8:如果采用单步收敛模式,设ε=ε/2,若未达精度要求,返回步骤5,否则转入步骤9。
步骤9:设t=t+1,若t<T-1,则转入步骤4,否则转入步骤10。
步骤10:若未收敛,转入步骤3;否则若采用单步收敛模式,设ε=ε/2,若未达到精度要求,转入步骤3,否则转入步骤11。
步骤11:寻优结束,输出结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611022560.2A CN106682810B (zh) | 2016-11-17 | 2016-11-17 | 巨型水电站动态投产下跨流域梯级水电站群长期运行方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611022560.2A CN106682810B (zh) | 2016-11-17 | 2016-11-17 | 巨型水电站动态投产下跨流域梯级水电站群长期运行方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106682810A true CN106682810A (zh) | 2017-05-17 |
CN106682810B CN106682810B (zh) | 2021-06-18 |
Family
ID=58867609
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611022560.2A Active CN106682810B (zh) | 2016-11-17 | 2016-11-17 | 巨型水电站动态投产下跨流域梯级水电站群长期运行方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106682810B (zh) |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107563642A (zh) * | 2017-09-01 | 2018-01-09 | 河海大学 | 一种基于投影寻踪的水电站河流水文情势评价方法 |
CN107844871A (zh) * | 2017-12-16 | 2018-03-27 | 博维恩冷冻科技(苏州)有限公司 | 一种电站发电效益优化方法 |
CN108108838A (zh) * | 2017-12-18 | 2018-06-01 | 华电福新能源股份有限公司福建分公司 | 一种高水量利用率的季调节水库优化调度方法 |
CN108133104A (zh) * | 2017-12-21 | 2018-06-08 | 清华大学 | 一种长期跨流域多梯级水电优化运行模拟方法 |
CN109636140A (zh) * | 2018-11-27 | 2019-04-16 | 广东电网有限责任公司韶关供电局 | 考虑径流量的微电网中小水电站弃水电量计算方法 |
CN109687506A (zh) * | 2018-11-27 | 2019-04-26 | 广东电网有限责任公司韶关供电局 | 微电网中小水电站发电量预测方法 |
CN110932337A (zh) * | 2019-12-03 | 2020-03-27 | 昆明电力交易中心有限责任公司 | 一种基于弃水比例均衡的电量调度方法、装置及设备 |
CN111428970A (zh) * | 2020-03-05 | 2020-07-17 | 三峡大学 | 一种大规模水电站群跨省外送能力分析模型及求解方法 |
CN111612268A (zh) * | 2020-05-28 | 2020-09-01 | 国家电网公司西南分部 | 考虑市场交易的龙头水库水电集群运行优化方法 |
CN111612282A (zh) * | 2020-06-30 | 2020-09-01 | 广东电网有限责任公司梅州供电局 | 一种地区性水电站的日负荷预测方法、系统及介质 |
CN111724253A (zh) * | 2020-05-21 | 2020-09-29 | 中国南方电网有限责任公司 | 一种水电交易执行偏差调度方法、系统、装置及存储介质 |
CN114781764A (zh) * | 2022-06-21 | 2022-07-22 | 国能大渡河流域水电开发有限公司 | 一种梯级水电站枯水期调度多模式自适应匹配方法及系统 |
CN117639111A (zh) * | 2024-01-25 | 2024-03-01 | 南京南瑞水利水电科技有限公司 | 一种基于梯级径流式水电的光伏波动平滑控制方法与系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102867275A (zh) * | 2012-08-14 | 2013-01-09 | 贵州乌江水电开发有限责任公司 | 梯级水库群中长期联合发电优化调度方法及系统 |
CN104967126A (zh) * | 2015-07-14 | 2015-10-07 | 华中科技大学 | 一种面向区域电网的跨流域水电站群多电网联合调峰方法 |
CN105869070A (zh) * | 2016-04-06 | 2016-08-17 | 大连理工大学 | 一种跨流域梯级水电站群效益均衡的合作优化调度方法 |
-
2016
- 2016-11-17 CN CN201611022560.2A patent/CN106682810B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102867275A (zh) * | 2012-08-14 | 2013-01-09 | 贵州乌江水电开发有限责任公司 | 梯级水库群中长期联合发电优化调度方法及系统 |
CN104967126A (zh) * | 2015-07-14 | 2015-10-07 | 华中科技大学 | 一种面向区域电网的跨流域水电站群多电网联合调峰方法 |
CN105869070A (zh) * | 2016-04-06 | 2016-08-17 | 大连理工大学 | 一种跨流域梯级水电站群效益均衡的合作优化调度方法 |
Non-Patent Citations (2)
Title |
---|
吴家乐: "流域梯级电站集中控制和联合优化调度浅析", 《水电厂自动化》 * |
申建建 等: "区域电网水火电系统短期协调控制方法及应用实践", 《电力自动化设备》 * |
Cited By (19)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107563642A (zh) * | 2017-09-01 | 2018-01-09 | 河海大学 | 一种基于投影寻踪的水电站河流水文情势评价方法 |
CN107844871A (zh) * | 2017-12-16 | 2018-03-27 | 博维恩冷冻科技(苏州)有限公司 | 一种电站发电效益优化方法 |
CN108108838B (zh) * | 2017-12-18 | 2021-08-27 | 华电福新能源股份有限公司福建分公司 | 一种高水量利用率的季调节水库优化调度方法 |
CN108108838A (zh) * | 2017-12-18 | 2018-06-01 | 华电福新能源股份有限公司福建分公司 | 一种高水量利用率的季调节水库优化调度方法 |
CN108133104A (zh) * | 2017-12-21 | 2018-06-08 | 清华大学 | 一种长期跨流域多梯级水电优化运行模拟方法 |
CN109636140A (zh) * | 2018-11-27 | 2019-04-16 | 广东电网有限责任公司韶关供电局 | 考虑径流量的微电网中小水电站弃水电量计算方法 |
CN109687506A (zh) * | 2018-11-27 | 2019-04-26 | 广东电网有限责任公司韶关供电局 | 微电网中小水电站发电量预测方法 |
CN110932337A (zh) * | 2019-12-03 | 2020-03-27 | 昆明电力交易中心有限责任公司 | 一种基于弃水比例均衡的电量调度方法、装置及设备 |
CN110932337B (zh) * | 2019-12-03 | 2021-05-18 | 昆明电力交易中心有限责任公司 | 一种基于弃水比例均衡的电量调度方法、装置及设备 |
CN111428970A (zh) * | 2020-03-05 | 2020-07-17 | 三峡大学 | 一种大规模水电站群跨省外送能力分析模型及求解方法 |
CN111724253A (zh) * | 2020-05-21 | 2020-09-29 | 中国南方电网有限责任公司 | 一种水电交易执行偏差调度方法、系统、装置及存储介质 |
CN111612268B (zh) * | 2020-05-28 | 2021-11-30 | 国家电网公司西南分部 | 考虑市场交易的龙头水库水电集群运行优化方法 |
CN111612268A (zh) * | 2020-05-28 | 2020-09-01 | 国家电网公司西南分部 | 考虑市场交易的龙头水库水电集群运行优化方法 |
CN111612282A (zh) * | 2020-06-30 | 2020-09-01 | 广东电网有限责任公司梅州供电局 | 一种地区性水电站的日负荷预测方法、系统及介质 |
CN111612282B (zh) * | 2020-06-30 | 2022-12-20 | 广东电网有限责任公司梅州供电局 | 一种地区性水电站的日负荷预测方法、系统及介质 |
CN114781764A (zh) * | 2022-06-21 | 2022-07-22 | 国能大渡河流域水电开发有限公司 | 一种梯级水电站枯水期调度多模式自适应匹配方法及系统 |
CN114781764B (zh) * | 2022-06-21 | 2022-10-18 | 国能大渡河流域水电开发有限公司 | 一种梯级水电站枯水期调度多模式自适应匹配方法及系统 |
CN117639111A (zh) * | 2024-01-25 | 2024-03-01 | 南京南瑞水利水电科技有限公司 | 一种基于梯级径流式水电的光伏波动平滑控制方法与系统 |
CN117639111B (zh) * | 2024-01-25 | 2024-04-09 | 南京南瑞水利水电科技有限公司 | 一种基于梯级径流式水电的光伏波动平滑控制方法与系统 |
Also Published As
Publication number | Publication date |
---|---|
CN106682810B (zh) | 2021-06-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106682810A (zh) | 巨型水电站动态投产下跨流域梯级水电站群长期运行方法 | |
Ming et al. | Optimizing utility-scale photovoltaic power generation for integration into a hydropower reservoir by incorporating long-and short-term operational decisions | |
Syahputra et al. | Renewable energy systems based on micro-hydro and solar photovoltaic for rural areas: A case study in Yogyakarta, Indonesia | |
Yuan et al. | Stochastic optimization model for the short-term joint operation of photovoltaic power and hydropower plants based on chance-constrained programming | |
Yang et al. | Economic-emission dispatch problem in integrated electricity and heat system considering multi-energy demand response and carbon capture Technologies | |
Wang et al. | Coordinated operation of conventional hydropower plants as hybrid pumped storage hydropower with wind and photovoltaic plants | |
CN104268653B (zh) | 基于集束径流预报的梯级水库优化调度方法 | |
CN103683337B (zh) | 一种互联电网cps指令动态分配优化方法 | |
Wen et al. | Optimizing the sizes of wind and photovoltaic plants complementarily operating with cascade hydropower stations: Balancing risk and benefit | |
CN109492861A (zh) | 一种梯级水电站群中期电量交易计划分解方法 | |
CN106485346A (zh) | 一种混联水库提前蓄水多目标优化调度方法 | |
CN105335561A (zh) | 一种基于指标排序的梯级水电站群超短期调度方法 | |
CN105427017A (zh) | 一种水电富集电网特大规模电站群短期计划编制方法 | |
CN104537428A (zh) | 一种计及风电接入不确定性的经济运行评估方法 | |
Helseth et al. | Long-term hydro-thermal scheduling including network constraints | |
CN105260801B (zh) | 一种水电富集电网大规模电站群长期电力电量平衡分析方法 | |
CN115271304A (zh) | 基于气候变化的跨流域调水工程优化调度方法 | |
Zhang et al. | Assessing the integration potential of new energy in river basin clean energy corridors considering energy-power coupled complementary operation modes | |
Wang et al. | Hybrid time-scale optimal scheduling considering multi-energy complementary characteristic | |
Jia et al. | Optimal operation of cascade hydro-wind-photovoltaic complementary generation system with vibration avoidance strategy | |
CN104636831A (zh) | 一种面向多电网的水电站短期调峰特征值搜索方法 | |
Ding et al. | Long-term operation rules of a hydro–wind–photovoltaic hybrid system considering forecast information | |
CN114123313A (zh) | 一种时序生产模拟新能源电力系统消纳方法 | |
Zhang et al. | Self-optimization simulation model of short-term cascaded hydroelectric system dispatching based on the daily load curve | |
CN117013535A (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |