CN104182634A - 一种梯级水电站联合运行水位控制断面优选方法 - Google Patents
一种梯级水电站联合运行水位控制断面优选方法 Download PDFInfo
- Publication number
- CN104182634A CN104182634A CN201410416309.9A CN201410416309A CN104182634A CN 104182634 A CN104182634 A CN 104182634A CN 201410416309 A CN201410416309 A CN 201410416309A CN 104182634 A CN104182634 A CN 104182634A
- Authority
- CN
- China
- Prior art keywords
- section
- water level
- upstream
- water
- overbar
- 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
Abstract
本发明公开了一种梯级水电站联合运行水位控制断面优选方法。本发明利用梯级水电站两坝间河道的水下地形及河道断面数据,建立一维水动力学模型,并进行上游水电站不同下泄流量和下游水电站不同坝前水位组合工况下沿程水位过程线计算,生成两坝间河道各待选断面的水位集合,对各待选断面的水位集合,计算其断面水位稳定性、断面水位对于上下游水位的敏感性、断面水位与上下游水位的相关性等特征值,并结合梯级水电站两坝间河道特点,优选一个合适的控制断面。本发明能够量化两坝间河道不同断面水位在不同来水频率下的变化特征,优选出的控制断面的水位作为上下梯级水电站联合发电优化调度快捷建模的决策变量。
Description
技术领域
本发明属于水力学与统计学的交叉技术领域,更具体地,涉及一种梯级水电站联合发电优化运行水位控制断面优选方法。该方法能够为流域管理机构进行梯级联合发电调度快捷建模提供决策变量。
背景技术
梯级水电站间存在水力和电力联系,不仅存在水流迟滞,某些梯级水电站间还存在水头衔接关系,因此,梯级间水位控制对梯级水电站联合运行、充分利用河流水位落差十分重要。通常,需分别采用上下游水电站坝前水位作为各水电站的决策变量。但由于梯级间洪水传播存在延时,还需进行两坝间河道洪水演算或通过流量平移进行简化处理,计算量十分庞大。因此,两坝间某断面水位若能正确反映不同来水情势下两水电站坝前水位的关系,将会使梯级联合发电优化调度计算快捷和高效。如何合理的选择梯级间河道某一断面水位来代表两坝前水位,直接关系到梯级上下游水电站发电水头和出力的计算精度、以及水电站实际生产中工况的快速调节性能,对实现梯级联合优化调度、提高梯级总出力具有重要意义。
发明内容
针对现有技术的以上缺陷或改进需求,本发明提供一种梯级水电站联合运行水位控制断面优选方法,通过梯级水电站两坝间沿程水面线水力学模拟生成了待选断面水位集合,建立了断面水位特征值计算公式,采用梯级水电站两坝间一维水动力学模拟,生成两坝间待选断面在不同流量方案下的水位集合,据此计算两坝间河道不同断面水位的稳定性、敏感性、相关性等三种特征值,并结合梯级水电站两坝间河道特点,优选一个合适的控制断面。该断面水位作为梯级水电站联合发电优化调度快捷建模的决策变量,能够实现梯级电站联合发电优化调度的实时响应和控制。
本发明提供一种梯级水电站联合运行水位控制断面优选方法,包括以下步骤:
步骤1在梯级水电站间河道上选择多个待选断面a1~aN,N表示所述待选断面个数;
步骤2根据所述待选断面a1~aN以及上游入口断面水位和下游出口断面水位,选择m个不同的组合工况,其中每个组合工况包括一个上游水电站下泄流量和一个下游水电站坝前水位,针对每个组合工况进行沿程水位过程线计算,得到所述待选断面a1~aN的水位集合h1~hN、所述上游入口断面的水位集合hup和所述下游出口断面的水位集合hdown:
hi=[hi,1,hi,2,…,hi,m],i=1,…,N
hup=[hup_1,hup_2,…,hup_m]
hdown=[hdown_1,hdown_2,…,hdown_m]
其中,hi,j、hup_j和hdown_j分别表示待选断面ai、所述上游入口断面和所述下游出口断面在第j种组合工况下的水位,j=1,2,…,m;m表示所述组合工况的数量;
步骤3在获得所述水位集合h1~hN后,进行所述待选断面a1~aN水位的特征值计算,其中,所述特征值包括以下三种:
(1)断面水位稳定性:
其中,Εi表示所述待选断面ai水位的方差,指示所述待选断面ai水位的稳定性;表示所述待选断面ai根据所述步骤2得到的水位集合hi中各水位的平均值;
(2)断面水位对于上下游水位的敏感性:
其中,ηi为所述待选断面ai水位的敏感性;
(3)断面水位与上下游水位的相关性:
ri=rup_i+rdown_i,(i=1,…,N)
其中:
其中,ri表示所述待选断面ai与所述上下游水位的相关性;rup_i表示所述待选断面ai与所述上游入口断面水位的相关系数;rdown_i表示所述待选断面ai与所述下游出口断面水位的相关系数;和分别表示根据所述步骤2计算得到的所述上游入口断面水位集合hup和所述下游出口断面水位集合hdown中各水位的平均值;
步骤4:根据各个待选断面的特征值(Ei,ηi,ri),选出合适的控制断面,其中,Ei表示所述待选断面ai的稳定性;ηi表示所述待选断面ai对于所述上下游水位的敏感性;ri表示所述待选断面ai与所述上下游水位的相关性。
总体而言,通过本发明所构思的以上技术方案与现有技术相比,具有以下有益效果:
本发明运用水力学方法建立梯级电站两坝间水力联系,生成梯级水电站两坝间河道待选断面的水位集合,提出评价待选断面水位集合的三种特征值,并根据两坝间河道特点优选出一个合适的控制断面。本发明提出的断面水位的特征值及其评价方法考虑了断面水位与上下游电站运行水位间的联系,优选出断面的水位可作为梯级水电站联合发电优化调度的水位控制对象,具有实用性。
附图说明
图1为本发明梯级水电站联合运行水位控制断面优选实施流程图;
图2为本发明待选断面划分示意图;
图3为本发明水力学计算点布置图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
图1所示为本发明梯级水电站运行水位控制断面位置优选实施流程图,具体包括以下步骤:
步骤1:根据河道地形数据和流域长系列水文资料,在梯级水电站间河道上选择多个待选断面a1~aN,N表示待选断面个数。图2所示为本发明待选断面划分示意图,在上游水电站和下游水电站间,选择多个待选断面a1~aN,其中待选断面的数量N无具体限制。待选断面应是有实测数据的河道断面和汊点(干支流交汇点)断面,与现有水文站点和拟建水文站点所在的断面尽量接近,以保证待选断面具有较高的工程应用意义。
步骤2:根据待选断面a1~aN以及上游入口断面和下游出口断面水位,选择m个不同的组合工况(每个组合工况包括一个上游水电站下泄流量和一个下游水电站坝前水位,应包含水库所有常用工况),针对这m个组合工况进行沿程水位过程线计算,得到待选断面a1~aN的水位集合h1~hN、上游入口断面的水位集合hup和下游出口断面的水位集合hdown:
hi=[hi,1,hi,2,…,hi,m],i=1,…,N
hup=[hup_1,hup_2,…,hup_m]
hdown=[hdown_1,hdown_2,…,hdown_m]
其中,hi,j、hup_j和hdown_j分别表示待选断面ai、上游入口断面和下游出口断面在第j种组合工况下的水位,j=1,2,…,m;m表示组合工况的数量。
在本发明实施例中,每个组合工况下的梯级间河道待选断面a1~aN的水位以及上游入口断面和下游出口断面的水位,都通过圣维南方程组描述,并对其离散后进行求解:
其中,Q表示断面流量;A表示断面过水面积;x表示距水道参考断面沿流程的距离;t表示计算时段序号;q表示单位长度河流的侧向入流或出流流量;α表示动量分布系数,默认取1.0;g表示重力加速度;C表示谢才系数;R表示断面的阻力半径;h表示断面水位。
为求解上述圣维南方程组,本发明采用图3所示的水力学计算点布置方式。在本发明实施例中,依次布置h点(水位节点)与Q点(流量节点)。具体布置方法为,在河道的入口断面、出口断面、有实测数据的河道断面以及汊点断面各设置一个h点,即h点完全包含待选断面a1~aN以及上游入口断面和下游出口断面;在两个h点的中点处布置Q点,如此得到n个h点和n-1个Q点。接下来,应用上述圣维南方程组对这2n-1个点采用水动力计算常用的Abbott格式进行离散,得到由2n-1个方程构成的线性方程组,对于每个节点的线性方程如下所示:
其中,k表示节点序号,k=1,2,…,2n-1;t表示计算时段序号;对于h点,和分别表示第k个h点前后2个Q点的流量值,表示第k个h点的水位值;对于Q点,和分别表示第k个Q点前后2个h点的水位值,表示第k个Q点的流量值;αj、βj、γj、δj分别表示圣维南方程组采用Abbott格式离散所求得的线性方程系数。
本发明采用追赶法对上述线性方程组进行求解,最终得到各h点处的水位值和各Q点处的流量值,由于h点包含了待选断面a1~aN以及上游入口断面和下游出口断面,那么在此组合工况下的断面a1~aN以及上游入口断面和下游出口断面的水位即已求出。接下来将每种组合工况都进行上述计算,即可求得待选断面a1~aN的水位集合h1~hN,以及上游入口断面的水位集合hup和下游出口断面的水位集合hdown。
步骤3:对各个待选断面a1~aN的水位集合进行特征值计算。在本发明实施例中采用了三种特征值,如下所示:
(1)断面水位稳定性,值越小表示其稳定性越好:
其中,Εi表示待选断面ai水位的方差,指示断面水位的稳定性;hi,j表示待选断面ai在第j种组合工况下的水位;表示待选断面ai根据上述方法得到的水位集合hi中各水位的平均值。
(2)断面水位对于上下游水位的敏感性,值越大表示其敏感性越好:
其中,ηi为待选断面ai水位的敏感性。
(3)断面水位与上下游水位的相关性,其值越接近2表示其相关性越好:
ri=rup_i+rdown_i,(i=1,…,N)
其中:
其中,ri表示待选断面ai与上下游水位的相关性;rup_i表示待选断面ai与上游入口断面水位的相关系数;rdown_i表示待选断面ai与下游出口断面水位的相关系数;hup_j和hdown_j分别表示第j种组合工况下的上游入口断面水位和下游出口断面水位;hi,j表示待选断面ai在第j种组合工况下的水位;和分别表示根据上述水动力学方法计算得到的上游入口断面的水位集合hup和下游出口断面的水位集合hdown中各水位的平均值。
步骤4:根据各个待选断面的特征值,选出合适的控制断面。由以上步骤得到的待选断面ai的特征值形式为:
(Ei,ηi,ri)
其中,Ei表示待选断面ai的稳定性;ηi表示待选断面ai的敏感性;ri表示待选断面ai的相关性。
在本发明实施例中,为了兼顾上述3种特征,并保证一定的选择空间,当梯级水电站两坝间河道地形复杂、流量较大时,选择稳定性Ei数值由低到高排序前20%的断面作为备选断面集;当梯级水电站两坝间河道较长、支流较多时,选择敏感性ηi数值由高到低排序排名前20%的断面作为备选断面集。然后,在备选断面集中,选择具有相关性ri数值最高的断面作为运行水位控制断面。
本发明在梯级水电站间河道上选择多个待选断面,利用梯级水电站间河道的水下地形及河道断面数据,建立了两坝间河道一维水动力学模型,并运用该模型生成各待选断面在不同来水频率下的水位集合。在此基础上,计算各待选断面水位的稳定性、敏感性和相关性等三种特征值。然后,结合两坝间河道特点,依据稳定性和敏感性特征值生成备选断面集,再依据相关性特征值选出了一个合适的控制断面。该断面水位能够合理、充分地反映梯级水电站上游水电站尾水位和下游水电站坝前水位之间的对应关系,为流域梯级联合发电调度生产运行提供一个实时便捷的控制对象。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种梯级水电站联合运行水位控制断面优选方法,其特征在于,包括:
步骤1在梯级水电站间河道上选择多个待选断面a1~aN,N表示所述待选断面个数;
步骤2根据所述待选断面a1~aN以及上游入口断面水位和下游出口断面水位,选择m个不同的组合工况,其中每个组合工况包括一个上游水电站下泄流量和一个下游水电站坝前水位,针对每个组合工况进行沿程水位过程线计算,得到所述待选断面a1~aN的水位集合h1~hN、所述上游入口断面的水位集合hup和所述下游出口断面的水位集合hdown:
hi=[hi,1,hi,2,…,hi,m],i=1,…,N
hup=[hup_1,hup_2,…,hup_m]
hdown=[hdown_1,hdown_2,…,hdown_m]
其中,hi,j、hup_j和hdown_j分别表示待选断面ai、所述上游入口断面和所述下游出口断面在第j种组合工况下的水位,j=1,2,…,m;m表示所述组合工况的数量;
步骤3在获得所述水位集合h1~hN后,进行所述待选断面a1~aN水位的特征值计算,其中,所述特征值包括以下三种:
(1)断面水位稳定性:
其中,Εi表示所述待选断面ai水位的方差,指示所述待选断面ai水位的稳定性;表示所述待选断面ai根据所述步骤2得到的水位集合hi中各水位的平均值;
(2)断面水位对于上下游水位的敏感性:
其中,ηi为所述待选断面ai水位的敏感性;
(3)断面水位与上下游水位的相关性:
ri=rup_i+rdown_i,(i=1,…,N)
其中:
其中,ri表示所述待选断面ai与所述上下游水位的相关性;rup_i表示所述待选断面ai与所述上游入口断面水位的相关系数;rdown_i表示所述待选断面ai与所述下游出口断面水位的相关系数;和分别表示根据所述步骤2计算得到的所述上游入口断面水位集合hup和所述下游出口断面水位集合hdown中各水位的平均值;
步骤4:根据各个待选断面的特征值(Ei,ηi,ri),选出合适的控制断面,其中,Ei表示所述待选断面ai的稳定性;ηi表示所述待选断面ai对于所述上下游水位的敏感性;ri表示所述待选断面ai与所述上下游水位的相关性。
2.如权利要求1所述的方法,其特征在于,所述步骤2中计算不同组合工况下的梯级间河道沿程水位过程线,采用如下圣维南方程组描述两坝间的水力学过程:
其中,Q表示断面流量;A表示断面过水面积;x表示距水道参考断面沿流程的距离;t表示计算时段序号;q表示单位长度河流的侧向入流或出流;α表示动量分布系数;g表示重力加速度;C表示谢才系数;R表示断面的阻力半径;h表示断面水位。
3.如权利要求2所述的方法,其特征在于,在河道的入口断面、出口断面、有实测数据的河道断面以及汊点断面各设置一个水位节点,在两个所述水位节点的中点处布置流量节点。
4.如权利要求3所述的方法,其特征在于,应用所述圣维南方程组对所述水位节点和所述流量节点进行离散化,得到线性方程组,对于每个节点的线性方程如下所示:
其中,k表示节点序号,k=1,2,…,2n-1;t表示计算时段序号;对于所述水位节点,和分别表示第k个水位节点前后2个流量节点的流量值,表示第k个水位节点的水位值;对于所述流量节点,和分别表示第k个流量节点前后2个水位节点的水位值,表示第k个流量节点的流量值;αj、βj、γj、δj分别表示所述圣维南方程组离散化所求得的线性方程系数。
5.如权利要求1-4中任一项所述的方法,其特征在于,当梯级水电站两坝间河道地形复杂、流量较大时,根据所述稳定性Ei数值在所述待选断面a1~aN中选择多个断面作为备选断面集;当梯级水电站两坝间河道较长、支流较多时,根据所述敏感性ηi数值在所述待选断面a1~aN中选择多个断面作为所述备选断面集,然后,在所述备选断面集中,根据所述相关性ri数值选择一个断面作为所述控制断面。
6.如权利要求1-4中任一项所述的方法,其特征在于,当梯级水电站两坝间河道地形复杂、流量较大时,选择所述稳定性Ei数值由低到高排序前20%的断面作为备选断面集;当梯级水电站两坝间河道较长、支流较多时,选择所述敏感性ηi数值由高到低排序排名前20%的断面作为所述备选断面集,然后,在所述备选断面集中,选择所述相关性ri数值最高的断面作为所述控制断面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410416309.9A CN104182634B (zh) | 2014-08-21 | 2014-08-21 | 一种梯级水电站联合运行水位控制断面优选方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410416309.9A CN104182634B (zh) | 2014-08-21 | 2014-08-21 | 一种梯级水电站联合运行水位控制断面优选方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104182634A true CN104182634A (zh) | 2014-12-03 |
CN104182634B CN104182634B (zh) | 2015-07-22 |
Family
ID=51963669
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410416309.9A Active CN104182634B (zh) | 2014-08-21 | 2014-08-21 | 一种梯级水电站联合运行水位控制断面优选方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104182634B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107288092A (zh) * | 2016-09-22 | 2017-10-24 | 武汉大学 | 冲积通航河道枯水河槽形态调整预测模型的构建方法 |
CN109837862A (zh) * | 2019-02-27 | 2019-06-04 | 长江水利委员会水文局 | 基于生态水位衔接的河流梯级开发设计和调度运用方法 |
CN110929929A (zh) * | 2019-11-19 | 2020-03-27 | 天津大学 | 一种考虑下游河道水力安全的水电站泄洪优化调控方法 |
CN111126847A (zh) * | 2019-12-24 | 2020-05-08 | 华中科技大学 | 耦合河道水动力过程的梯级水库短期优化调度方法和系统 |
WO2021129039A1 (zh) * | 2019-12-25 | 2021-07-01 | 华中科技大学 | 一种跨水库流域河道径流演算方法和系统 |
CN113887154A (zh) * | 2021-10-21 | 2022-01-04 | 中国电建集团成都勘测设计研究院有限公司 | 一种梯级水电站群中长期发电能力评估方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009139622A1 (en) * | 2008-05-13 | 2009-11-19 | Van Den Noort Innovations Bv | Self closing flood barrier and method for protecting a hinterland using same |
CN102817335A (zh) * | 2012-07-12 | 2012-12-12 | 贵州乌江水电开发有限责任公司 | 一种梯级水库群联合防洪优化调度的方法及系统 |
CN102982221A (zh) * | 2012-12-26 | 2013-03-20 | 北京奥特美克科技股份有限公司 | 一种河流断面动态模拟方法及系统 |
-
2014
- 2014-08-21 CN CN201410416309.9A patent/CN104182634B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009139622A1 (en) * | 2008-05-13 | 2009-11-19 | Van Den Noort Innovations Bv | Self closing flood barrier and method for protecting a hinterland using same |
CN102817335A (zh) * | 2012-07-12 | 2012-12-12 | 贵州乌江水电开发有限责任公司 | 一种梯级水库群联合防洪优化调度的方法及系统 |
CN102982221A (zh) * | 2012-12-26 | 2013-03-20 | 北京奥特美克科技股份有限公司 | 一种河流断面动态模拟方法及系统 |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107288092A (zh) * | 2016-09-22 | 2017-10-24 | 武汉大学 | 冲积通航河道枯水河槽形态调整预测模型的构建方法 |
CN109837862A (zh) * | 2019-02-27 | 2019-06-04 | 长江水利委员会水文局 | 基于生态水位衔接的河流梯级开发设计和调度运用方法 |
CN110929929A (zh) * | 2019-11-19 | 2020-03-27 | 天津大学 | 一种考虑下游河道水力安全的水电站泄洪优化调控方法 |
CN111126847A (zh) * | 2019-12-24 | 2020-05-08 | 华中科技大学 | 耦合河道水动力过程的梯级水库短期优化调度方法和系统 |
CN111126847B (zh) * | 2019-12-24 | 2022-08-02 | 华中科技大学 | 耦合河道水动力过程的梯级水库短期优化调度方法和系统 |
WO2021129039A1 (zh) * | 2019-12-25 | 2021-07-01 | 华中科技大学 | 一种跨水库流域河道径流演算方法和系统 |
CN113887154A (zh) * | 2021-10-21 | 2022-01-04 | 中国电建集团成都勘测设计研究院有限公司 | 一种梯级水电站群中长期发电能力评估方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104182634B (zh) | 2015-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104182634B (zh) | 一种梯级水电站联合运行水位控制断面优选方法 | |
Wu et al. | Projected climate change impacts on the hydrology and temperature of Pacific Northwest rivers | |
Chang et al. | Efficiency evaluation of hydropower station operation: A case study of Longyangxia station in the Yellow River, China | |
CN105893763B (zh) | 南方感潮河网区水流水质耦合模拟系统 | |
CN109992909B (zh) | 树状河网梯级水库水动力水质泥沙耦合模拟方法及系统 | |
CN105243502B (zh) | 一种基于径流区间预测的水电站调度风险评估方法及系统 | |
Quaranta et al. | CFD simulations to optimize the blade design of water wheels | |
Li et al. | Identifying explicit formulation of operating rules for multi-reservoir systems using genetic programming | |
CN103729550A (zh) | 基于传播时间聚类分析的多模型集成洪水预报方法 | |
CN104047258B (zh) | 一种面向环境保护的联合调度水库群的生态库容确定方法 | |
CN108614915B (zh) | 基于情景驱动的水文模型自由组建策略方法 | |
CN103455725A (zh) | 管网系统非恒定流模拟方法 | |
Chen et al. | Understanding flow dynamics and density currents in a river-reservoir system under upstream reservoir releases | |
CN108647778A (zh) | 一种城市雨水系统排水口排水流量的动态预测方法 | |
Dorchies et al. | Climate change impacts on multi-objective reservoir management: case study on the Seine River basin, France | |
Baláž et al. | On the use of the Muskingum method for the simulation of flood wave movements | |
Shouyu et al. | Fuzzy Optimization Neural Network Approach for Ice Forecast in the Inner Mongolia Reach of the Yellow River/Approche d'Optimisation Floue de Réseau de Neurones pour la Prévision de la Glace Dans le Tronçon de Mongolie Intérieure du Fleuve Jaune | |
Müller et al. | Flow field in a reservoir subject to pumped-storage operation–in situ measurement and numerical modeling | |
Bao et al. | Dynamic correction of roughness in the hydrodynamic model | |
CN106284239B (zh) | 一种大型坝式水电站初步选址方法 | |
CN103088783B (zh) | 一种面向生态的水库调度方案生成方法 | |
Davidsen et al. | The cost of ending groundwater overdraft on the North China Plain | |
Attab | Estimation of Manning’s roughness coefficient for Tigris River by using HEC-RAS model | |
Szydłowski et al. | Hydropower potential of the lower Vistula | |
CN106228277A (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 | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |