CN109933892B - 一种水温模型电站下泄流量边界条件的修正方法 - Google Patents
一种水温模型电站下泄流量边界条件的修正方法 Download PDFInfo
- Publication number
- CN109933892B CN109933892B CN201910182654.3A CN201910182654A CN109933892B CN 109933892 B CN109933892 B CN 109933892B CN 201910182654 A CN201910182654 A CN 201910182654A CN 109933892 B CN109933892 B CN 109933892B
- Authority
- CN
- China
- Prior art keywords
- sum
- flow
- day
- power station
- calculation
- 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.)
- Active
Links
Images
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
- Other Liquid Machine Or Engine Such As Wave Power Use (AREA)
Abstract
本发明公开了一种水温模型电站下泄流量边界条件的修正方法,包括:步骤1,统计电站可能出现的发电流量并按从小到大排序,依次为Q1,Q2,…;步骤2,定义流量累积变量Qsum并赋予初始值0,赋予儒略日编号n初始值1;步骤3,对于qn作如下计算,并在每一天的计算之后进入步骤4:Qsum=Qsum+qn;步骤4,若Qsum<Q1,则该日流量qn=0,进入步骤6;否则i=2并进入步骤5;步骤5,若Qsum>Qi+0.5×Qi‑1,则该日流量qn=Qi,Qsum=Qsum-Qi,进入步骤6;否则i=i+1,再次执行步骤5;步骤6,若计算至儒略日的最后一天,则计算结束;否则n=n+1,返回步骤3。本发明在不改变电站下泄发电的前提下,最大程度地还原了电站发电可能出现的真实流量工况,显著提高了模型精度。
Description
技术领域
本发明涉及工程领域含有隔水幕墙的宽度平均立面二维水库水温模型建模计算,特别涉及一种水温模型电站下泄流量边界条件的修正方法。
背景技术
高坝大库易在春夏形成水库水温分层的现象,水温较高的水***于上层,温度较低的水***于下层。而电站取水口多高程较低,所取水温为低温水,造成春夏季下泄低温水的问题。电站春夏季下泄低温水,将会影响下游河道的鱼类繁殖和农作物生长,造成严重的生态影响。
隔水幕墙是提升电站下泄水温的工程措施之一,其原理是在取水口上游修建一道位于水下的隔水幕墙,挡住底层低温水,使上层温度较高的水体进入电站进水口,达到提升电站下泄水温的目的。
隔水幕墙在设计时,需建立全库、全年尺度的宽度平均立面二维水库水温模型,模拟计算不同设计参数、不同断面位置的隔水幕墙实施后,水库各月水温结构的变化,以及下泄水温和幕墙的水温改善效果。
通常宽度平均立面二维水库水温模型在输入模型边界条件时,均采用日均值。但是对于含有隔水幕墙的宽度平均立面二维水库水温模型,若发电流量仍然采用日均值,将导致较大的计算误差。其原因在于,水库底部多崎岖不平,而隔水幕墙的修建为水下施工,现有技术水平无法做到底部完全封闭,会留有1m~2m的间隙。库区底层低温水则会由此进入电站进水口,影响下泄水温。
根据坝前局部三维水库水温模型的研究成果,电站发电时,分流比(从隔水幕墙顶部进入电站的高温水流量与从隔水幕墙底部进入电站的低温水流量之比)非固定值,与电站的发电流量相关。发电流量越大,分流比越大,即进入电站的高温水占比越大。
由于电站发电机组在修建后,其发电流量多为固定值。即可能出现的发电流量为少数固定值。电站在实际运行时,多数情况下,一天之内仅有一台机组运行几个小时,其余时候并无下泄流量。在这种情况下,电站的日均发电流量将低于单台机组的发电流量。若将该日均流量作为宽度平均立面二维水库水温模型的发电流量边界条件。那么模型内将会计算为一台发电流量为日均流量的机组运行了24小时。因此,在模型中虽然一日之内总的下泄水量,但机组发电流量低于实际的机组发电流量,机组发电时长长于实际的发电时长。由于不同发电流量下,隔水幕墙的分流比不同,那么在宽度平均立面二维水库水温模型中,将有比实际情况更多的低温水从库区流出。这将导致宽度平均立面二维水库水温模型模拟所得的下泄水温严重偏低,库区水温结构严重失真。
针对上述问题,本发明提出一种针对含有隔水幕墙的宽度平均立面二维水库水温模型计算时,电站下泄水温边界条件的修正方法。
发明内容
本发明的目的在于,针对上述现有技术的不足,提供一种水温模型电站下泄流量边界条件的修正方法,在保证电站年内总发电水量不变的前提下,还原电站出流现状的下泄流量,达到让隔水幕墙分流比接近实际情况的目的。
为解决上述技术问题,本发明所采用的技术方案是:
一种水温模型电站下泄流量边界条件的修正方法,其特点是包括以下步骤:
步骤1,根据电站机组的发电流量,统计电站可能出现的发电流量,并按从小到大排序,依次为Q1,Q2,…;
步骤2,定义流量累积变量Qsum,并赋予初始值:Qsum=0,赋予儒略日编号n初始值n=1;
步骤3,对于qn作如下计算,并在每一天的计算之后进入步骤4:
Qsum=Qsum+qn,
式中,qn为第n天的日均发电流量;
步骤4,若Qsum<Q1,则该日流量qn=0,进入步骤6;否则i=2并进入步骤5;
步骤5,若Qsum>Qi+0.5×Qi-1,则该日流量qn=Qi,Qsum=Qsum-Qi,进入步骤6;否则i=i+1,再次执行步骤5;
步骤6,若计算至儒略日的最后一天,则计算结束;否则n=n+1,返回步骤3。
与现有技术相比,本发明通过算法对水库下泄流量过程进行调整,在不改变电站下泄发电的前提下,最大程度地还原了电站发电可能出现的真实流量工况,极大地降低了数学模型边界条件采用日均值计算时带来的隔水幕墙分流比计算误差,显著提高了模型精度。
附图说明
图1是发明的方法流程示意图。
具体实施方式
本发明的计算流程如图1所示,计算步骤如下:
步骤1,根据电站机组的发电流量,统计电站可能出现的发电流量,并按从小到大排序,依次为Q1,Q2,…;
步骤2,定义流量累积变量Qsum,并赋予初始值:Qsum=0,赋予儒略日编号n初始值n=1;
步骤3,对于qn作如下计算,并在每一天的计算之后进入步骤4:
Qsum=Qsum+qn,
式中,qn为第n天的日均发电流量;
步骤4,若Qsum<Q1,则该日流量qn=0,进入步骤6;否则i=2并进入步骤5;
步骤5,若Qsum>Qi+0.5×Qi-1,则该日流量qn=Qi,Qsum=Qsum-Qi,进入步骤6;否则i=i+1,再次执行步骤5;
步骤6,若计算至儒略日的最后一天,则计算结束;否则n=n+1,返回步骤3。
以三板溪隔水幕墙全库宽度平均立面二维水温模型计算为例,选择连续10个模拟的儒略日,实施步骤如下:
步骤1,三板溪电站共4台机组,可能出现的发电流量为220m3/s、440m3/s、660m3/s、870m3/s;
步骤2,定义流量累积变量Qsum,并赋予初始值:Qsum=0,赋予儒略日编号n初始值n=1;
步骤4~步骤6为逐日的循环计算,计算过程见表1:
表1发电流量修正表单位:m3/s
儒略日 | 日均流量(修正前) | <![CDATA[Q<sub>sum</sub>]]> | 日均流量(修正后) |
1 | 34 | 34 | 0 |
2 | 145 | 179 | 220 |
3 | 280 | 239 | 220 |
4 | 576 | 595 | 660 |
5 | 364 | 299 | 220 |
6 | 683 | 762 | 660 |
7 | 278 | 380 | 440 |
8 | 43 | -17 | 0 |
9 | 124 | 107 | 0 |
10 | 34 | 141 | 220 |
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是局限性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护范围之内。
Claims (1)
1.一种水温模型电站下泄流量边界条件的修正方法,其特征在于,包括以下步骤:
步骤1,根据电站机组的发电流量,统计电站可能出现的发电流量,并按从小到大排序,依次为Q1,Q2,…;
步骤2,定义流量累积变量Qsum,并赋予初始值:Qsum=0,赋予儒略日编号n初始值n=1;
步骤3,对于qn作如下计算,并在每一天的计算之后进入步骤4:
Qsum=Qsum+qn,
式中,qn为第n天的日均发电流量;
步骤4,若Qsum<Q1,则第n天的日均发电流量qn=0,进入步骤6;否则i=2并进入步骤5;
步骤5,若Qsum>Qi+0.5×Qi-1,则第n天的日均发电流量qn=Qi,Qsum=Qsum-Qi,进入步骤6;否则i=i+1,再次执行步骤5;
步骤6,若计算至儒略日的最后一天,则计算结束;否则n=n+1,返回步骤3。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910182654.3A CN109933892B (zh) | 2019-03-12 | 2019-03-12 | 一种水温模型电站下泄流量边界条件的修正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910182654.3A CN109933892B (zh) | 2019-03-12 | 2019-03-12 | 一种水温模型电站下泄流量边界条件的修正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109933892A CN109933892A (zh) | 2019-06-25 |
CN109933892B true CN109933892B (zh) | 2023-04-07 |
Family
ID=66986848
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910182654.3A Active CN109933892B (zh) | 2019-03-12 | 2019-03-12 | 一种水温模型电站下泄流量边界条件的修正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109933892B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112699521B (zh) * | 2020-11-12 | 2023-01-24 | 中国电建集团中南勘测设计研究院有限公司 | 一种隔水幕墙顶部过流高度年逐日变化数据的分析方法 |
CN112580230B (zh) * | 2020-11-12 | 2022-05-17 | 中国电建集团中南勘测设计研究院有限公司 | 一种水库水温结构分析方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102968674A (zh) * | 2012-12-05 | 2013-03-13 | 湖南省电力公司 | 一种基于水位控制的日调节水电站发电优化调度方法 |
CN104636831A (zh) * | 2015-02-12 | 2015-05-20 | 华中科技大学 | 一种面向多电网的水电站短期调峰特征值搜索方法 |
WO2017071230A1 (zh) * | 2015-10-30 | 2017-05-04 | 南京南瑞集团公司 | 一种多Agent的水电站群短期优化调度方法 |
CN106845712A (zh) * | 2017-01-20 | 2017-06-13 | 中国水利水电科学研究院 | 面向大***水库群的聚合分解调度规则的提取方法 |
CN108797533A (zh) * | 2018-06-14 | 2018-11-13 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种泥石流拦砂坝圆形泄流孔的设计方法 |
CN109272229A (zh) * | 2018-09-14 | 2019-01-25 | 中国长江电力股份有限公司 | 基于虚拟水库的电站偏差电量计算方法 |
-
2019
- 2019-03-12 CN CN201910182654.3A patent/CN109933892B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102968674A (zh) * | 2012-12-05 | 2013-03-13 | 湖南省电力公司 | 一种基于水位控制的日调节水电站发电优化调度方法 |
CN104636831A (zh) * | 2015-02-12 | 2015-05-20 | 华中科技大学 | 一种面向多电网的水电站短期调峰特征值搜索方法 |
WO2017071230A1 (zh) * | 2015-10-30 | 2017-05-04 | 南京南瑞集团公司 | 一种多Agent的水电站群短期优化调度方法 |
CN106845712A (zh) * | 2017-01-20 | 2017-06-13 | 中国水利水电科学研究院 | 面向大***水库群的聚合分解调度规则的提取方法 |
CN108797533A (zh) * | 2018-06-14 | 2018-11-13 | 中国科学院、水利部成都山地灾害与环境研究所 | 一种泥石流拦砂坝圆形泄流孔的设计方法 |
CN109272229A (zh) * | 2018-09-14 | 2019-01-25 | 中国长江电力股份有限公司 | 基于虚拟水库的电站偏差电量计算方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109933892A (zh) | 2019-06-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109933892B (zh) | 一种水温模型电站下泄流量边界条件的修正方法 | |
Liu et al. | Deriving optimal refill rules for multi-purpose reservoir operation | |
CN107491635B (zh) | 一种基于嵌套降维算法的梯级水库水沙联合优化调度方法 | |
CN110543731A (zh) | 一种高含沙河流水库水沙分置调度模拟方法及*** | |
CN112257352A (zh) | 一维水动力模型和二维水动力模型的耦合方法及*** | |
CN105574605A (zh) | 一种梯级泵站输水***日优化调控方法 | |
CN104318058B (zh) | 基于雨量监测的泥石流预警方法 | |
CN111080157B (zh) | 一种梯级水电站排磷量的调度方法和*** | |
CN110210109B (zh) | 一种河网中堰闸工程反向水流的数值模拟方法及*** | |
CN109101732A (zh) | 基于地形特征界线的无汊河道二维结构网格剖分方法 | |
CN111340350B (zh) | 长距离输水渠道实现渠道局部水力冲刷的闸门群调控方法 | |
CN116628915A (zh) | 一种库区沿程流量预测方法 | |
CN114819753B (zh) | 一种促进鱼类产卵的梯级水库生态调度方法 | |
CN107862408B (zh) | 一种水电厂弃水最小预警协调滚动优化方法 | |
CN103510487A (zh) | 一种具有季节性河流水质改善功能的景观跌水结构 | |
CN109636015B (zh) | 一种梯级水电虚拟抽蓄电站调度方法 | |
CN115186857A (zh) | 一种基于集成学习的神经网络水库水位预测方法 | |
CN111882116A (zh) | 梯级水库泥沙实时预报方法 | |
CN106651011A (zh) | 一种基于粒子群算法的渠系优化配水方法 | |
Zheng et al. | The effects of intra-annual variability of river discharge on the spatio-temporal dynamics of saltmarsh vegetation at river mouth bar: Insights from an ecogeomorphological model | |
CN113779671A (zh) | 一种基于时空步长自适应技术的明渠调水工程水动力实时计算方法 | |
CN110578317B (zh) | 一种水文模型水库泄水量模拟方法 | |
CN112580230B (zh) | 一种水库水温结构分析方法 | |
CN113361818B (zh) | 一种水电机组中长期发电优化方法及*** | |
CN116415824A (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 |