CN111046574B - 一种平原湖区湖河闸泵***防洪排涝演算方法 - Google Patents
一种平原湖区湖河闸泵***防洪排涝演算方法 Download PDFInfo
- Publication number
- CN111046574B CN111046574B CN201911323888.1A CN201911323888A CN111046574B CN 111046574 B CN111046574 B CN 111046574B CN 201911323888 A CN201911323888 A CN 201911323888A CN 111046574 B CN111046574 B CN 111046574B
- Authority
- CN
- China
- Prior art keywords
- lake
- flood
- gate
- water level
- river
- 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
- 238000004364 calculation method Methods 0.000 title claims abstract description 37
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 148
- 238000000034 method Methods 0.000 claims abstract description 43
- 230000008569 process Effects 0.000 claims abstract description 33
- 238000005086 pumping Methods 0.000 claims abstract description 18
- 238000001363 water suppression through gradient tailored excitation Methods 0.000 claims abstract description 3
- 238000003860 storage Methods 0.000 claims description 31
- 238000007599 discharging Methods 0.000 claims description 6
- 238000011144 upstream manufacturing Methods 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- 230000001133 acceleration Effects 0.000 claims 1
- 230000005484 gravity Effects 0.000 claims 1
- 238000009826 distribution Methods 0.000 abstract description 9
- 238000005094 computer simulation Methods 0.000 abstract description 6
- 238000011160 research Methods 0.000 description 6
- 125000000217 alkyl group Chemical group 0.000 description 5
- 238000010276 construction Methods 0.000 description 5
- 238000013461 design Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 241000532412 Vitex Species 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 235000009347 chasteberry Nutrition 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000002265 prevention Effects 0.000 description 2
- 238000009991 scouring Methods 0.000 description 2
- 241000317410 Arisaema dracontium Species 0.000 description 1
- 241000195493 Cryptophyta Species 0.000 description 1
- 101100136092 Drosophila melanogaster peng gene Proteins 0.000 description 1
- 238000006424 Flood reaction Methods 0.000 description 1
- 244000248021 Vitex negundo Species 0.000 description 1
- 235000010363 Vitex negundo Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- HVCNNTAUBZIYCG-UHFFFAOYSA-N ethyl 2-[4-[(6-chloro-1,3-benzothiazol-2-yl)oxy]phenoxy]propanoate Chemical compound C1=CC(OC(C)C(=O)OCC)=CC=C1OC1=NC2=CC=C(Cl)C=C2S1 HVCNNTAUBZIYCG-UHFFFAOYSA-N 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000002262 irrigation Effects 0.000 description 1
- 238000003973 irrigation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000013486 operation strategy Methods 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
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
- 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
-
- 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
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Barrages (AREA)
Abstract
本发明公开了一种平原湖区湖河闸泵***防洪排涝演算方法,包括下述步骤:根据平原湖区水文测站的逐日降雨、逐日水位资料,采用降雨径流模型设计入湖洪水过程和入河道洪水过程;利用有侧向入流的明渠非恒定流基本方程计算两级串联水闸群首级闸闸下水位、末级闸闸上水位以及泵站站前水位、河道测点水位;建立水闸不同流态下的闸门出流模型计算水闸泄流能力;建立分洪区流量计算模型计算分洪流量;建立泵站排水模型计算泵站的排水能力。本发明能有效地由已知条件并指定调度策略,求解防洪排涝***状态结果的模型;并根据防洪排涝***状态计算机模拟结果提出大型排涝泵站调度和主要湖泊控制运用意见。
Description
技术领域
本发明涉及水资源与环境技术领域,尤其涉及一种平原湖区湖河闸泵***防洪排涝演算方法。
背景技术
我国南方湖区系江湖冲积平原,地形平坦,土地肥沃,雨量充沛,河网密布、湖泊众多,是重要的农副业生产基地。但由于地面高程大多低于外江洪水位,防洪安全依赖堤防保护。汛期降雨量大而集中,湖区洪水受外江持续高水位顶托,自流排水困难,易形成溃涝灾害。经过治理,大部分湖区已形成有调蓄区(含调蓄湖泊和备蓄区等)、河网、排水闸和泵站(包括外排泵站和内排泵站)等工程构成的湖河闸泵防洪排涝***。开展平原湖区湖河闸泵***防洪排涝演算模型研究对大型排涝泵站调度和主要湖泊控制运用意见编制、降低泵站群运行成本、提升湖区防洪排涝管理水平具有重要的理论意义和实用价值。
本发明将提出一种平原湖区湖河闸泵***防洪排涝演算模型。图1是较为常见的带分蓄洪区一湖两串联水闸一泵站型防洪排涝***概化图。***运行时各外排设施(外排闸、外排泵站)优先排出其区域内的田间涝水(称为排田);没有外排设施区域的排田流量进入湖泊调蓄;外排闸(末级水闸)关闭期间,当外排泵站站前水位高于泵站起排水位时外排泵站开机排水(称为排湖),低于停机水位时泵站停机;若排水河道监测点水位达到河道分洪水位ZfR且预报有雨,则执行入河道洪水分洪调蓄方案;若外排泵站全部开机,湖泊水位仍不能稳定在湖泊分洪水位ZfL且预报有雨,为保证湖堤及下游人民生命财产安全,执行湖泊汛期控制运用分洪调蓄方案。
目前常见的排水排涝***有:平原圩区除涝***、城市港渠闸站排涝***、灌区内部排水除涝***,对于这些类似的水资源***调度的研究多采用***分析的方法。***分析方法是一种用于研究复杂***的组合构成及运行规律的概念性方法,通常来讲它包括了确定***的构成要素、组织架构、对象关系、功能效用等基本要素,以及运行过程和联系等动态要素。***分析方法已经被广泛的用于平原湖区的除涝排水优化运行和城市河湖闸泵群排涝***优化调度问题研究,并且采用了多种计算方法,比如线性规划、非线性规划、动态规划、大***分解协调、多目标规划等等,建立了较为科学准确的数学模型或运算程序,提高了生产效益同时,为未来的规划建设,运行策略也给予了科学合理的指导。
研究水资源***的规划建设、运营方式、调度方案等优化问题的目的是为了更好的指导其未来和现在的运行,使***的“硬件设施”能够充分发挥其功能效用。国内外的经验也表明,一个出色的水资源***,通常是科学合理的规划建设和优化的调度运行策略相辅相成的。到目前为止,国内外无论是对平原湖区除涝***还是城市的河湖闸泵群排涝***优化调度的研究大多侧重于建设规划或项目论证,对于平原湖区防洪排涝设施组成的大排水***的优化排涝调度运行方面研究较少。
发明内容
本发明要解决的技术问题在于针对现有技术中的缺陷,提供一种平原湖区湖河闸泵***防洪排涝演算方法,能有效地由已知条件并指定调度操作,求解排水***状态结果的模型,并根据排水***状态计算机模拟结果提出大型排涝泵站调度和主要湖泊控制运用意见修编建议。
本发明解决其技术问题所采用的技术方案是:
本发明提供一种平原湖区湖河闸泵***防洪排涝演算方法,该方法包括以下步骤:
步骤1:收集平原湖区水文测站的逐日降雨、逐日水位资料,采用降雨径流模型设计入湖洪水过程Q(t)和入河道洪水过程Q2(t);
步骤2:输入湖泊汛限水位Z0;泵站起排水位Zyp;湖泊分洪水位ZfL;河道分洪水位ZfR;湖泊分洪量W1max;河道分洪量W2max;t=0天;
步骤3:令t=t+1开始迭代计算;
步骤4:输入入湖洪水过程Q(t)和入河道洪水过程Q2(t)和末级闸闸下外江水位Zd (2)(t);
步骤5:假设t时段末湖泊水位ZL(t),由时段初的湖泊水位开始试算,时段初的湖泊水位取湖泊汛限水位Z0;
步骤6:假设t时段末首级闸排湖流量Q1 ′(t);
步骤7:由防洪排涝***概化模型确定首级闸排湖流量Q1 ′(t)、河道侧向入流q=Q2(t)-q(t);再根据有侧向入流的明渠非恒定流基本方程计算首级闸闸下河道水位Zd (1)(t)、末级闸闸上水位Zu (1)(t)、泵站站前水位Zp(t)、河道测点水位ZR(t);
步骤8:根据水闸不同流态下的闸门出流模型计算首级闸泄流能力Q1(t);
步骤9:若|Q1(t)-Q′1(t)|≤ε,ε为给定的收敛精度,则迭代收敛,首级闸泄流能力Q1(t)即为所求;否则,重新设Q′1(t)=Q1(t)+Δ1,返回步骤6,进行下次迭代,Δ1为给定增量;若Q′1(t)>Q1(t),则Δ1>0;否则Δ1<0;
步骤10:由湖泊水量平衡模型计算时段末湖泊容积变化:ΔV’=V′(t)-V′(t-1);
步骤11:由时段初的湖泊水位,然后根据湖泊水位-容积关系曲线求得时段初湖泊容积V′(t-1),进而得到时段末湖泊容积V′(t);再据湖泊水位-容积关系曲线模型求得时段末湖泊水位ZL′(t);
步骤12:若ZR(t)≥ZfR且W2(t)<W2max,ZfR为分洪调蓄方案确定的河道水位测点分洪水位,W2max为分洪调蓄方案确定的河道侧向入流分洪量,则执行分洪调蓄方案,由分洪区流量计算模型确定河道侧向入流分洪流量q(t+1),并更新计算W2(t),重新计算河道侧向入流q=Q2(t)-q(t),返回步骤7;否则往下继续执行;
步骤13:若ZL′(t)≥ZfL且W1(t)<W1max,ZfL为分洪调蓄方案确定的湖泊分洪水位,W1max为湖泊分洪量,则执行分洪调蓄方案,由分洪区流量计算模型确定湖泊分洪流量q’(t),并更新计算W1(t),返回步骤10;否则往下继续执行;
步骤14:若|ZL(t)-Z'L(t)|≤ε,ε为给定的收敛精度,则迭代收敛,时段末湖泊水位ZL′(t)即为所求;否则,重新设ZL(t)=ZL′(t)+Δ2,返回步骤4,进行下次迭代,Δ2为给定增量;若ZL(t)>ZL′(t),则Δ2>0;否则Δ2<0;
步骤15:计算末级闸泄流过程q1(t)、泵站排水能力Dp(t);
步骤16:若t≥tmax,则时段末湖泊水位ZL′(t)、首级闸排湖流量Q1(t)即为所求;否则,返回步骤3;
步骤17:输出湖泊水位过程ZL′(t)、外排闸末级泄流过程q1(t)、外排泵站排水过程Dp(t)。
进一步地,本发明的步骤7中,有侧向入流的明渠非恒定流基本方程为:
其中,Q1′(t)表示t时段末首级闸排湖流量,m3/s;q表示河道侧向入流,q=Q2(t)-q(t),m3/s;Q2(t)表示时段末入河洪水过程,m3/s;q(t)表示时段末河道侧向入流分洪流量,m3/s;A表示河道过水断面面积,m2;Z表示水位,m;R表示水力半径,m;C表示谢才系数;α表示动量修正系数。
进一步地,本发明的步骤8中,判断水闸不同流态的流程为:
平板闸门出流流态的判别按时段内闸上水头H、计入行近流速水头的闸上水头H0、闸下水深hs和闸门开闸高度e的相互关系来界定:
1)e/H≤0.65且hs/e≤1,为闸孔自由出流;
2)e/H≤0.65且1<hs/e≤5.1,为闸孔不完全淹没出流;
3)e/H≤0.65且hs/e>5.1,为闸孔完全淹没出流;
4)e/H>0.65且hs/H0≤0.80,为宽顶堰自由出流;
5)e/H>0.65且0.80<hs/H0<0.90,为宽顶堰低淹没出流;
6)e/H>0.65且0.90≤hs/H0,为宽顶堰高淹没出流。
进一步地,本发明的步骤8中,不同流态下的闸门出流模型具体为:
其中,H(1)表示首级闸上水头,m;行近流速小于一定阈值时,略去行近流速水头,有:e(1)、B(1)分别为首级闸开闸高度、闸门宽度,m;表示t时段末首级闸闸下水深,m;表示t时段末首级闸上游水位,m;表示首级闸底板高程,为已知条件,m;表示t时段末首级闸下游水位,m;ZL(t)表示t时段末湖泊水位,m。
进一步地,本发明的步骤10中,湖泊水量平衡模型的方程为:
其中,Q(t-1)、Q(t)分别为t时段初、末入湖洪水过程,m3/s;Q1(t-1)、Q1(t)分别为t时段初、末首级闸泄流量,m3/s;q′(t-1)、q′(t)为t时段初、末湖泊分洪流量,m3/s;V'(t-1)、V'(t)为t时段初、末湖泊容积,m3;t表示时间,以日计,Δt=86400s。
进一步地,本发明的步骤11中,湖泊水位-容积关系曲线为:
Z'L(t)=F(V'(t)) (5)
其中,Z'L(t)为t时段末湖泊的水位,m;V'(t)为t时段末湖泊容积,m3;F(V'(t))为湖泊的水位容积关系曲线。
进一步地,本发明的步骤12中更新计算W2(t)的方法为:W2(t)=W2(t)+q(t)*24*0.36,单位万m3;步骤13中更新计算W1(t)的方法为:W1(t)=W1(t)+q’(t)*24*0.36,单位万m3。
进一步地,本发明的步骤12中,分洪区流量计算模型为:
其中,B1为分蓄洪区分洪口门宽度,m;H1为分蓄洪区分洪口门上游水头,m;行近流速小于一定阈值时,略去行近流速水头,有:H10≈H1;Z1f分蓄洪区分洪口门底板高程,m。
进一步地,本发明的步骤15中,通过泵站排水模型计算泵站排水能力,泵站排水模型为:
其中,Dp(t)、Qp(t)分别为泵站在t时段末的总抽排流量和来流量,m3/s;Dpb为泵站的抽排能力,m3/s;Zp(t-1)为t时段初的泵站站前水位,m;Zyp为泵站的起排水位,m。
本发明产生的有益效果是:本发明的平原湖区湖河闸泵***防洪排涝演算方法,相对于现有平原湖区湖河闸泵***防洪排涝演算模型,能有效地由已知条件并指定调度操作,求解排水***状态结果的模型,并根据排水***状态计算机模拟结果提出大型排涝泵站调度和主要湖泊控制运用意见修编建议。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1为本发明实施例的流程图;
图2为本发明实施例中典型平原湖区湖河闸泵防洪排涝***示意图;
图3为本发明实施例中闸门出流流态判别流程图;
图4为本发明实施例中四湖流域在湖北省的位置示意图;
图5为本发明实施例中长湖流域防洪排涝***现状示意图;
图6为本发明实施例中长湖流域湖河闸泵防洪排涝***概化图;
图7为本发明实施例中长湖流域30日最大暴雨频率曲线图;
图8为本发明实施例中长湖水位-容积关系曲线;
图9为本发明实施例中计算机模拟2016年6月23日-7月26日长湖水位与实测水位的比较图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅用以解释本发明,并不用于限定本发明。
本实施例以湖北长湖流域为例对平原湖区湖河闸泵防洪排涝***联合调度计算机模拟模型进行阐述。
四湖流域在湖北省的位置见图4。地处长江中游荆江北岸、汉江及其东荆河以南,为江汉平原心腹地带。因境内原有四个大型湖泊(长湖、三湖、白露湖、洪湖)而得名。目前仅保留长湖、洪湖两个湖泊。在行政区划上包括荆州市的荆州区、沙市区、江陵县、监利县、洪湖市和石首市的一部分及荆门市的沙洋县、潜江市一部分,流域总面积11547.5平方公里,内垸面积10375平方公里,洲滩民垸面积1172.5平方公里。
根据水系地理特征及等高截流、分散治理的原则,四湖流域分为上、中、下三区。三个区既相互独立、又相互联系,统筹兼顾,分区负担调蓄及防洪排涝任务。长湖、田关河以上的地区,称为上区,又称为长湖流域,流域面积3240km2,其中丘陵山区面积2360km2,占上区的72.8%。长湖、田关河以下,洪湖、下新河、排涝河以上的地区为中区,汇流面积5980km2,均为平原。洪湖、下新河、排涝河以下的地区为下区,汇流面积1155km2。
四湖流域常见的自然灾害有洪、涝、旱、渍等灾害,但最严重的以洪涝最为频繁。2016年7月23日长湖水位达到33.45m,超过历史最高水位0.15米,接近长湖50年一遇洪水位(33.50m)。长湖内垸9.7km2内垸全部分洪,对长湖沿岸的城镇和农业水产业造成了极大损失,也增加了长湖围堤的防汛劳务用工和费用。亟待开展对长湖防洪排涝调度决策支持***研究。
长湖流域防洪排涝***现状示意图见图5。长湖流域又分为长湖及田北两片,长湖片面积为2265.5km2,田北片面积为974.5km2。其中沙洋境内西荆河和大港河来水可经双店渠汇入长湖,亦可直接汇入田关河;潜江境内来水则直接汇入田关河。经过多年的规划和建设,长湖流域目前已建成内蓄与外排相结合、自排与提排相结合、内垸分洪与外部分洪相结合的防御工程体系。
长湖调蓄工程:(1)长湖相当于一座大型平原水库,自然汇流面积为2265.5km2,人工汇流面积为974.5km2。根据《湖北省长湖湖泊保护规划》得到的最新湖泊详查资料,长湖水位33.00m,对应水面面积为130km2,相应湖容为5.64亿m3。(2)长湖库堤西起沙市沙桥门东至荆门市蝴蝶咀,堤防总长49.391km;滨湖内垸防洪堤,包括太湖港南堤,太湖港北堤和纪南防洪堤,桥河堤、后港围堤、蛟尾围堤和藻湖围堤等,堤防总长约62.27km。
长湖外排工程:长湖常规外排工程主要为田关河排水***工程,包括长湖主要出湖闸、排水干渠和外排闸站。(1)长湖主要出湖闸为刘岭闸,位于田关河渠首,建于1965年,设计流量为229m3/s。(2)田关河是长湖的主要排水干渠,过流能力达250m3/s,全长约30.46km,目前渠道底宽120m,堤顶高程达34.50m。(3)田关闸和田关泵站位于田关河尾端,是长湖流域的主要外排出口。田关闸位于东荆河堤右岸,于1965年重建,设计排水流量250m3/s;田关泵站位于田关河尾端田关闸右侧,装机6×2800kW,设计排水流量220m3/s,于1989年建成。
内分洪工程:(1)长湖内垸包括乔子湖外垸、马子湖、胜利垸、外六台、幸福垸等,分洪面积9.7平方公里,分洪量0.19亿立方米。(2)彭塚湖分蓄洪区:位于荆门市境,包括北湖(彭湖)、南湖(潘家湖)和宋湖等,计划蓄洪面积7.6km2,分蓄洪2280万m3,目前采取扒口分洪方式。(3)借粮湖分蓄洪区:蓄洪水位33.0m,蓄洪面积53.0km2,可分蓄洪1~1.2亿m3。目前采取窑场扒口分洪方式。
外分洪工程:在非常规情况下,长湖的洪水可通过习家口闸和高场南闸、张义咀闸、中沙河闸等向四湖流域中下区分洪。分洪的途径主要有以下两条:(1)总干分洪工程:首端建筑物为习家口闸,位于总干渠渠首,设计流量107m3/s。(2)东干分洪工程:首端建筑物为高场南闸、张义咀闸、中沙河闸等,主要向中下区的东干渠分洪,分洪总流量可达120m3/s。
以下通过2016年防洪调度验证本发明的正确性。
采用平原湖区湖河闸泵防洪排涝***联合调度计算机模拟模型,再根据2016年实际防洪调度情况进行长湖流域的湖河闸泵防洪排涝***联合调度计算机模拟,计算2016年计算期的长湖水位;通过与计算期长湖实测水位的比较,验证上述湖河闸泵防洪排涝***联合调度计算机模拟的可靠性。
长湖流域30日最大暴雨频率曲线见图7。2016年6月23日~7月22日长湖流域30d累计面雨量499.7mm。
表1长湖流域2016年6月23日~7月22日30d降雨
日期(月/日) | 6/23 | 6/24 | 6/25 | 6/26 | 6/27 | 6/28 | 6/29 | 6/30 |
面雨量(mm) | 6.0 | 29.2 | 12.8 | 16.7 | 0.8 | 0.0 | 0.0 | 93.3 |
日期(月/日) | 7/1 | 7/2 | 7/3 | 7/4 | 7/5 | 7/6 | 7/7 | 7/8 |
面雨量(mm) | 35.5 | 7.5 | 42.7 | 1.3 | 9.5 | 7.5 | 0.0 | 0.0 |
日期(月/日) | 7/9 | 7/10 | 7/11 | 7/12 | 7/13 | 7/14 | 7/15 | 7/16 |
面雨量(mm) | 0.0 | 2.7 | 0.0 | 0.0 | 14.2 | 43.3 | 7.2 | 0.0 |
日期(月/日) | 7/17 | 7/18 | 7/19 | 7/20 | 7/21 | 7/22 | ||
面雨量(mm) | 0.8 | 43.7 | 115.2 | 9.8 | 0.0 | 0.0 |
以日为时段查暴雨径流相关曲线得出相应净雨见表。
表2长湖流域2016年6月23日~7月22日30d净雨
日期(月/日) | 6/23 | 6/24 | 6/25 | 6/26 | 6/27 | 6/28 | 6/29 | 6/30 |
面雨量(mm) | 0.0 | 12.2 | 6.9 | 9.3 | 0.5 | 0.0 | 0.0 | 56.1 |
日期(月/日) | 7/1 | 7/2 | 7/3 | 7/4 | 7/5 | 7/6 | 7/7 | 7/8 |
面雨量(mm) | 22.5 | 4.8 | 27.7 | 0.9 | 6.2 | 4.9 | 0.0 | 0.0 |
日期(月/日) | 7/9 | 7/10 | 7/11 | 7/12 | 7/13 | 7/14 | 7/15 | 7/16 |
面雨量(mm) | 0.0 | 1.8 | 0.0 | 0.0 | 9.4 | 29.0 | 4.8 | 0.0 |
日期(月/日) | 7/17 | 7/18 | 7/19 | 7/20 | 7/21 | 7/22 | ||
面雨量(mm) | 0.6 | 29.7 | 80.0 | 6.9 | 0.0 | 0.0 |
再根据长湖片和田北片经验单位线分别进行长湖片和田北片汇流计算,分别得到长湖片和田北片2016年洪水过程,见表3。
表3长湖片和田北片2016年设计洪水过程,单位:m3/s
序号 | 长湖片 | 田北片 |
0 | 0 | 0 |
1 | 50 | 11 |
2 | 161 | 41 |
3 | 191 | 82 |
4 | 177 | 117 |
5 | 96 | 121 |
6 | 43 | 87 |
7 | 252 | 90 |
8 | 716 | 157 |
9 | 626 | 233 |
10 | 446 | 262 |
11 | 463 | 240 |
12 | 291 | 175 |
13 | 208 | 149 |
14 | 144 | 142 |
15 | 70 | 109 |
16 | 33 | 72 |
17 | 18 | 30 |
18 | 26 | 5 |
19 | 13 | 0 |
20 | 43 | 22 |
21 | 223 | 80 |
22 | 398 | 140 |
23 | 261 | 162 |
24 | 116 | 136 |
25 | 183 | 80 |
26 | 686 | 80 |
27 | 1109 | 150 |
28 | 659 | 212 |
29 | 278 | 207 |
30 | 140 | 141 |
31 | 79 | 54 |
32 | 39 | 8 |
33 | 3 | 0 |
请见图1,本发明提供的一种平原湖区湖河闸泵***防洪排涝演算方法,包括以下步骤:
步骤1:收集四湖流域上区四个水文测站的逐日降雨资料。采用降雨径流模型设计长湖片洪水过程Q(t)和田北片洪水过程Q2(t)。
步骤2:
1)输入长湖汛限水位Z0=30.50m;
2)田关泵站起排水位Zyp=31.00m;
3)长湖分洪水位ZfL(根据《湖北省大型排涝泵站调度及主要湖泊控制运用意见(鄂政发〔2011〕74号)》调度规则:33.00m,但2016年7月22日8时长湖水位33.35m,超过历史最高水位0.05m才分洪。);
4)西荆河向彭冢湖分洪水位ZfR=33.00m(根据《湖北省大型排涝泵站调度及主要湖泊控制运用意见(鄂政发〔2011〕74号)》调度规则:当西荆河流量超过250m3/s,高场水位接近33.00m,预报后续洪水流量更大,高场水位超过33.00m,则***西荆河右岸(青龙闸附近)堤防,向彭冢湖分洪。);
5)西荆河向彭冢湖分洪量W2max(彭冢湖分蓄洪区分洪面积7.6km2,蓄洪量0.228亿m3。);
W2max=0.228亿m3 (8)
6)湖泊分洪量W1max(长湖内垸包括乔子湖外垸、马子湖、胜利垸(含外六台)、幸福垸等,分洪面积9.7km2,分洪量0.19亿m3;彭冢湖分蓄洪区分洪面积7.6km2,蓄洪量0.228亿m3;借粮湖分蓄洪区分洪面积53km2,分洪量1至1.2亿m3。)
t=0天。
步骤3:令t=t+1开始迭代计算。
步骤4:输入长湖片洪水过程Q(t)和田北片洪水过程Q2(t)和田关闸闸下东荆河水位Zd (2)(t)。
步骤5:假设t时段末长湖水位ZL(t)(2016年6月22日长湖水位30.28m)。
步骤6:假设t时段末刘岭闸排湖流量Q1′(t)。
步骤7:由防洪排涝***概化图6确定刘岭闸排湖流量Q1′(t)、田关河侧向入流q=Q2(t)-q(t)。再根据有侧向入流的明渠非恒定流基本方程计算刘岭闸闸下水位Zd (1)(t)、田关闸闸上水位Zu (1)(t)、田关泵站站前水位Zp(t)、高场水位ZR(t)。
步骤8:根据水闸不同流态下的闸门出流模型计算刘岭闸泄流能力Q1(t)。刘岭闸出流流态判断流程图见图2。
式中:
e(1)、B(1)——分别为刘岭闸开闸高度、闸门宽度,m。
ZL(t)——假定t时段末长湖水位,m。
步骤9:若|Q1(t)-Q′1(t)|≤ε(ε为给定的收敛精度),则迭代收敛,刘岭闸泄流能力Q1(t)即为所求;否则,重新设Q1′(t)=Q1(t)+Δ1,返回步骤6,进行下次迭代,Δ1为给定增量。若Q1′(t)>Q1(t),则Δ1>0;否则Δ1<0。
步骤10:由长湖水量平衡模型计算时段末长湖容积变化ΔV’=V′(t)-V′(t-1)。
式中:
Q(t-1)、Q(t)——分别为t时段初、末长湖片洪水过程,m3/s。
Q1(t-1)、Q1(t)——分别为t时段初、末刘岭闸泄流量,m3/s。
q′(t-1)、q′(t)——为t时段初、末长湖分洪流量,m3/s。
V'(t-1)、V'(t)——为t时段初、末长湖容积,m3。
V'(t-1)、V'(t)——为t时段初、末长湖容积,m3。
步骤11:由2016年6月22日长湖水位30.28m,根据长湖水位-容积关系曲线求得时段初长湖容积V′(t-1),进而得到时段末长湖容积V′(t)。再据长湖水位-容积关系曲线模型求得时段末湖泊水位ZL′(t)。
长湖水位-容积关系曲线见图7。
步骤12:若ZR(t)≥ZfR且W2(t)<W2max,则执行分洪调蓄方案,由分洪区流量计算模型确定西荆河向彭冢湖分洪流量q(t+1),计算W2(t)=W2(t)+q(t)*24*0.36(万m3),重新计算田关河侧向入流q=Q2(t)-q(t),返回步骤7;否则往下继续执行(ZfR为高场分洪水位33.00m,W2max为分洪调蓄方案确定的西荆河向彭冢湖分洪量0.228亿m3)。
步骤13:若ZL′(t)≥ZfL且W1(t)<W1max,则执行分洪调蓄方案,由分洪区流量计算模型确定湖泊分洪流量q’(t),计算W1(t)=W1(t)+q’(t)*24*0.36(万m3),返回步骤10;否则往下继续执行(根据《湖北省大型排涝泵站调度及主要湖泊控制运用意见(鄂政发〔2011〕74号)》调度规则:33.00m,但2016年7月22日8时长湖水位33.35m,超过历史最高水位0.05m才分洪。W1max为长湖分洪量)。
步骤14:若|ZL(t)-Z'L(t)|≤ε(ε为给定的收敛精度),则迭代收敛,时段末长湖水位ZL′(t)即为所求;否则,重新设ZL(t)=ZL′(t)+Δ2,返回步骤5,进行下次迭代,Δ2为给定增量。若ZL(t)>ZL′(t),则Δ2>0;否则Δ2<0。
步骤15:计算田关闸泄流过程q1(t)、田关泵站排水能力Dp(t)。
田关闸出流流态判断流程图见图2,根据水闸不同流态下的闸门出流模型计算田关闸泄流能力q1(t)。
式中:
e(2)、B(2)——分别为田关闸开闸高度、闸门宽度,m。
ZL(t)——假定t时段末长湖水位,m。
步骤17:输出长湖水位过程ZL′(t)。
根据《湖北省大型排涝泵站调度及主要湖泊控制运用意见(鄂政发〔2011〕74号)》调度规则,只运用长湖内垸、彭冢湖分洪,不运用借粮湖分洪,计算2016年计算期的长湖水位与长湖实测水位的比较见图6。
1)2016年洪水,没有运用借粮湖分洪,长湖水位排涝演算结果与实测水位结果比较见图8。由图8知:对2016年洪水,通过运用长湖内垸分洪、彭冢湖分洪等措施,不能使长湖水位控制在33.00m(保证水位),最高可达33.46m,超保证水位0.46m。
2)上区2016年6月23日~7月22日防洪调度模拟结果见表4。
表4上区2016年6月23日~7月22日防洪调度模拟结果
3)本发明的平原湖区湖河闸泵***防洪排涝计算机模拟模型可以用于长湖流域调度运用方案的论证。
应当理解的是,对本领域普通技术人员来说,可以根据上述说明加以改进或变换,而所有这些改进和变换都应属于本发明所附权利要求的保护范围。
Claims (8)
1.一种平原湖区湖河闸泵***防洪排涝演算方法,其特征在于,该方法包括以下步骤:
步骤1:收集平原湖区水文测站的逐日降雨、逐日水位资料,采用降雨径流模型设计入湖洪水过程Q(t)和入河道洪水过程Q2(t);
步骤2:输入湖泊汛限水位Z0;泵站起排水位Zyp;湖泊分洪水位ZfL;河道分洪水位ZfR;湖泊分洪量W1max;河道分洪量W2max;t=0天;
步骤3:令t=t+1开始迭代计算;
步骤4:输入入湖洪水过程Q(t)和入河道洪水过程Q2(t)和末级闸闸下外江水位Zd (2)(t);
步骤5:假设t时段末湖泊水位ZL(t),由时段初的湖泊水位开始试算,时段初的湖泊水位取湖泊汛限水位Z0;
步骤6:假设t时段末首级闸排湖流量Q1′(t);
步骤7:由防洪排涝***概化模型确定首级闸排湖流量Q1′(t)、河道侧向入流q=Q2(t)-q(t);再根据有侧向入流的明渠非恒定流基本方程计算首级闸闸下河道水位Zd (1)(t)、末级闸闸上水位Zu (1)(t)、泵站站前水位Zp(t)、河道测点水位ZR(t);
步骤8:根据水闸不同流态下的闸门出流模型计算首级闸泄流能力Q1(t);
步骤9:若|Q1(t)-Q′1(t)|≤ε,ε为给定的收敛精度,则迭代收敛,首级闸泄流能力Q1(t)即为所求;否则,重新设Q1′(t)=Q1(t)+△1,返回步骤6,进行下次迭代,△1为给定增量;若Q1′(t)>Q1(t),则△1>0;否则△1<0;
步骤10:由湖泊水量平衡模型计算时段末湖泊容积变化:ΔV’=V′(t)-V′(t-1);
步骤11:由时段初的湖泊水位,然后根据湖泊水位-容积关系曲线求得时段初湖泊容积V′(t-1),进而得到时段末湖泊容积V′(t);再据湖泊水位-容积关系曲线模型求得时段末湖泊水位ZL′(t);
步骤12:若ZR(t)≥ZfR且W2(t)<W2max,ZfR为分洪调蓄方案确定的河道水位测点分洪水位,W2max为分洪调蓄方案确定的河道侧向入流分洪量,则执行分洪调蓄方案,由分洪区流量计算模型确定河道侧向入流分洪流量q(t+1),并更新计算W2(t),重新计算河道侧向入流q=Q2(t)-q(t),返回步骤7;否则往下继续执行;
步骤13:若ZL′(t)≥ZfL且W1(t)<W1max,ZfL为分洪调蓄方案确定的湖泊分洪水位,W1max为湖泊分洪量,则执行分洪调蓄方案,由分洪区流量计算模型确定湖泊分洪流量q’(t),并更新计算W1(t),返回步骤10;否则往下继续执行;
步骤14:若|ZL(t)-Z'L(t)|≤ε,ε为给定的收敛精度,则迭代收敛,时段末湖泊水位ZL′(t)即为所求;否则,重新设ZL(t)=ZL′(t)+△2,返回步骤4,进行下次迭代,△2为给定增量;若ZL(t)>ZL′(t),则△2>0;否则△2<0;
步骤15:计算末级闸泄流过程q1(t)、泵站排水能力Dp(t);
步骤16:若t≥tmax,则时段末湖泊水位ZL′(t)、首级闸排湖流量Q1(t)即为所求;否则,返回步骤3;
步骤17:输出湖泊水位过程ZL′(t)、外排闸末级泄流过程q1(t)、泵站排水能力Dp(t)。
3.根据权利要求1所述的平原湖区湖河闸泵***防洪排涝演算方法,其特征在于,步骤8中,判断水闸不同流态的流程为:
平板闸门出流流态的判别按时段内闸上水头H、计入行近流速水头的闸上水头H0、闸下水深hs和闸门开闸高度e的相互关系来界定:
1)e/H≤0.65且hs/e≤1,为闸孔自由出流;
2)e/H≤0.65且1<hs/e≤5.1,为闸孔不完全淹没出流;
3)e/H≤0.65且hs/e>5.1,为闸孔完全淹没出流;
4)e/H>0.65且hs/H0≤0.80,为宽顶堰自由出流;
5)e/H>0.65且0.80<hs/H0<0.90,为宽顶堰低淹没出流;
6)e/H>0.65且0.90≤hs/H0,为宽顶堰高淹没出流。
6.根据权利要求1所述的平原湖区湖河闸泵***防洪排涝演算方法,其特征在于,步骤11中,湖泊水位-容积关系曲线为:
Z'L(t)=F(V'(t))
其中,Z'L(t)为t时段末湖泊的水位,m;V'(t)为t时段末湖泊容积,m3;F(V'(t))为湖泊的水位容积关系曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911323888.1A CN111046574B (zh) | 2019-12-20 | 2019-12-20 | 一种平原湖区湖河闸泵***防洪排涝演算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201911323888.1A CN111046574B (zh) | 2019-12-20 | 2019-12-20 | 一种平原湖区湖河闸泵***防洪排涝演算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111046574A CN111046574A (zh) | 2020-04-21 |
CN111046574B true CN111046574B (zh) | 2022-03-11 |
Family
ID=70237339
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201911323888.1A Active CN111046574B (zh) | 2019-12-20 | 2019-12-20 | 一种平原湖区湖河闸泵***防洪排涝演算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111046574B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111709134B (zh) * | 2020-06-12 | 2021-05-04 | 河海大学 | 一种基于闸泵站群控制的多水源跨流域调度方法及装置 |
CN112101818B (zh) * | 2020-10-13 | 2024-04-23 | 南昌工程学院 | 一种适用于复杂水力联系的海绵城市洪涝优化调度方法 |
CN112431772A (zh) * | 2020-11-23 | 2021-03-02 | 王程 | 圩垸区防汛排涝泵站群优化调度运行方法 |
CN113204819B (zh) * | 2021-04-28 | 2024-02-09 | 浙江省水利河口研究院(浙江省海洋规划设计研究院) | 感潮水闸闸前橄榄型调蓄湖的布置方法 |
CN114819322B (zh) * | 2022-04-18 | 2024-05-14 | 武汉大学 | 湖泊入湖流量的预报方法 |
CN115423346B (zh) * | 2022-09-20 | 2023-04-07 | 中国水利水电科学研究院 | 一种明渠调水工程串联闸群汛期可调节能力实时测算方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103455725A (zh) * | 2013-09-06 | 2013-12-18 | 武汉大学 | 管网***非恒定流模拟方法 |
CN103543639A (zh) * | 2013-10-30 | 2014-01-29 | 武汉大学 | 河湖水体量质耦合自优化模拟调控方法 |
CN104978451A (zh) * | 2015-08-12 | 2015-10-14 | 河海大学 | 一种基于ArcGIS的复杂平原河网地区水文模型处理方法 |
US10094091B1 (en) * | 2015-09-02 | 2018-10-09 | John A. Tesvich | Sediment suction sink and method for sediment control in rivers, streams, and channels |
CN110276145A (zh) * | 2019-06-26 | 2019-09-24 | 天津神州海创科技有限公司 | 城市排水***仿真建模与调度方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AR109623A1 (es) * | 2018-02-16 | 2019-01-09 | Pescarmona Enrique Menotti | Proceso y sistema de análisis y gestión hidrológica para cuencas |
-
2019
- 2019-12-20 CN CN201911323888.1A patent/CN111046574B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103455725A (zh) * | 2013-09-06 | 2013-12-18 | 武汉大学 | 管网***非恒定流模拟方法 |
CN103543639A (zh) * | 2013-10-30 | 2014-01-29 | 武汉大学 | 河湖水体量质耦合自优化模拟调控方法 |
CN104978451A (zh) * | 2015-08-12 | 2015-10-14 | 河海大学 | 一种基于ArcGIS的复杂平原河网地区水文模型处理方法 |
US10094091B1 (en) * | 2015-09-02 | 2018-10-09 | John A. Tesvich | Sediment suction sink and method for sediment control in rivers, streams, and channels |
CN110276145A (zh) * | 2019-06-26 | 2019-09-24 | 天津神州海创科技有限公司 | 城市排水***仿真建模与调度方法 |
Non-Patent Citations (3)
Title |
---|
Physical Fractional Steps for Numerical Solution of Unsteady 2D Open Channel Flows;Sun Guangcai 等;《2010 International Conference on Computational and Information Sciences》;20110204;第854-856页 * |
小浪底水库射流冲吸式清淤设备研究;陆宏圻 等;《人民黄河》;20110430;第33卷(第4期);第15-18页 * |
近20年青海湖水量变化遥感分析;张洪源 等;《地理科学进展》;20180630;第823-832页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111046574A (zh) | 2020-04-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111046574B (zh) | 一种平原湖区湖河闸泵***防洪排涝演算方法 | |
CN110046469B (zh) | 多约束条件下水电站坝前河床冲淤变形的计算方法 | |
Feldman | HEC models for water resources system simulation: theory and experience | |
CN110298076A (zh) | 一种基于gis和swmm的城市内涝智能建模及分析方法 | |
CN113723024B (zh) | 一种适用于滨海地区的“溪流”-“河道”-“河口”分布式洪水过程模拟方法 | |
CN106777959A (zh) | 人工干扰无水文资料地区河流环境流量分区界定计算方法 | |
Nagy et al. | Hydrological dimensioning and operation of reservoirs: Practical design concepts and principles | |
CN107885958A (zh) | 一种平原感潮河网区纳污能力计算方法 | |
CN112101818B (zh) | 一种适用于复杂水力联系的海绵城市洪涝优化调度方法 | |
Wu et al. | Hydrological response to climate change and human activities: A case study of Taihu Basin, China | |
CN112052545B (zh) | 一种基于元胞自动机的城市地表径流与管网汇流耦合方法 | |
CN106284239B (zh) | 一种大型坝式水电站初步选址方法 | |
CN114997591A (zh) | 一种基于水环境数学模型的入河污染物削减及河道断面水质预测方法及装置 | |
Li et al. | Development of 1D and 2D coupled model to simulate urban inundation: an application to Beijing Olympic Village | |
CN110096670B (zh) | 调蓄池削减峰值流量的一种动态计算方法 | |
CN112926786B (zh) | 一种基于关联规则模型和数值模拟的浅水湖泊目标水位逆向预测方法及*** | |
CN109537520A (zh) | 一种建设项目防排洪设小水库留空库容削峰调洪***及方法 | |
CN111090711B (zh) | 河流水文地貌变化关键因素贡献度分析方法及*** | |
CN116305452A (zh) | 一种流域尺度淤地坝拦沙量的动态计算方法 | |
CN112818438B (zh) | 一种基于swmm的河道模型及其调度模拟概化方法 | |
CN110889185A (zh) | 一种小流域的洪峰流量分析方法及其应用 | |
CN111340649B (zh) | 一种水系结构连通性的量测方法 | |
Lu et al. | Modeling Monthly Fluctuations in Submersion Area of a Dammed River Reservoir: A Case Study 1 | |
CN113047212A (zh) | 一种城市山洪防治与内河生态补水结合的治水方法 | |
Nagy et al. | Planning for dams and reservoirs: Hydrologic design elements and operational characteristics of storage reservoirs |
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 |