海河流域下垫面变化对蓄滞洪区最高滞洪水位影响分析

2015-12-15 01:11谭桂秋
海河水利 2015年1期
关键词:文安滞洪区下垫面

齐 静,谭桂秋

(海河水利委员会科技咨询中心,天津300170)

海河流域下垫面变化对蓄滞洪区最高滞洪水位影响分析

齐 静,谭桂秋

(海河水利委员会科技咨询中心,天津300170)

近年来,海河流域下垫面情况发生了剧烈变化,加上气候变化的影响,流域总体洪涝水明显减少,局部城市沥涝水增加。这些变化对海河流域内蓄滞洪区的滞洪水位将会产生一定影响。选取东淀、文安洼、贾口洼3个蓄滞洪区进行数值模拟分析,比较不同频率情况下垫面变化对蓄滞洪区最高滞洪水位的影响。

蓄滞洪区;水力学模型;滞洪水位;流域

1 问题提出

海河流域曾是洪涝灾害比较严重的流域,但从1963年以来流域内一直未发生流域性的大洪水,水资源尤其是地下水资源的过度开发、土地利用方式的改变、城市化率的增加、工程建设等因素导致流域下垫面发生剧烈的变化,加之气候变化的影响,流域总体洪涝水明显减少,局部城市沥涝水增加。由此引起的洪水情势变化,会对流域内蓄滞洪区的相应蓄滞洪特征指标产生一定的影响。笔者以海河流域大清河中下游的东淀、文安洼、贾口洼3个蓄滞洪区为例,通过二维水动力学数值模拟,分析下垫面变化对以上3个蓄滞洪区最高滞洪水位的影响。

2 三大蓄滞洪区基本情况

2.1 东淀

东淀位于大清河下游河北省霸州市、文安县和天津市静海县、西青区境内,为大清河南北支洪水和清南、清北沥水的总归宿地,洪水由独流减河和海河干流入海。东淀属于纯行洪区,设计超过3年一遇洪水就要漫滩行洪。东淀遇50年一遇洪水滞洪水位(第六堡)6.44 m时,淹没村庄77个,涉及24 105户9.4万人;淹没土地376 km2,耕地1.58万hm2,房屋258万m2,GDP16亿元,固定资产15.4亿元;淹没区人口达蓄滞洪区总人口的98%。

2.2 文安洼

文安洼地处大清河下游,位于河北省文安县、大城县、任丘市及天津市静海县境内,承纳白洋淀及东淀分洪洪水和清南沥水,西、北靠千里堤、隔淀堤,东依子牙河左堤,南至津保公路南线。文安洼有3处进洪口门,分别是王村分洪闸、滩里分洪口门及小关分洪口门。退水主要靠滩里口门,经东淀由独流减河和海河干流下泄。文安洼规划分两区滞洪,Ⅰ区位于小白河以北和滩里新渠、赵磨支渠、任河大干渠以东,其余为Ⅱ区。当白洋淀水位超过9.0 m、出流大于2 700 m3/s时,超量洪水由赵王新渠上的王村闸向文安洼分洪;当东淀第六堡水位超过6.44 m时,向文安洼或贾口洼分洪。文安洼滞洪控制水位为大赵站5.94 m,当Ⅰ区、贾口洼、团泊洼已经运用,东淀第六堡水位超过6.44 m且继续上涨时,分别在马武营附近和文安县城防洪围堤以西扒开口门分洪入Ⅱ区。文安洼遇50年一遇洪水滞洪水位5.94 m时,淹没村庄529个,涉及153 671户58万人;淹没土地1 310 km2,耕地9.33万hm2,房屋1 489万m2,GDP101亿元,固定资产95.8亿元;淹没土地面积为本蓄滞洪区总土地面积的84%,95%的耕地被淹没,淹没区人数达蓄滞洪区总人口数的92%。

2.3 贾口洼

贾口洼位于东淀南侧,处于黑龙港河下游、子牙河与南运河之间,地跨天津市的静海县、河北省的大城县和青县,西北依子牙河右堤,东靠南运河左岸,南邻子牙新河左堤北侧自然高地。贾口洼兼具行洪与滞洪的功能,当遇到20年一遇洪水时启用。贾口洼遇50年一遇洪水滞洪水位5.94 m时,淹没村庄292个,涉及72 639户27.3万人;淹没土地815 km2,耕地7.32万hm2,房屋689万m2,GDP33.1亿元,固定资产53.0亿元;淹没土地面积为本蓄滞洪区总土地面积的89%,96%的耕地被淹没,淹没区人数达蓄滞洪区总人口数的93%。

三蓄滞洪区位置,如图1所示。

图1 三蓄滞洪区位置

3 设计洪水成果

3.1 现行设计洪水成果

现行的洪水设计成果采用河北省水利水电勘测设计院1990年12月完成的《大清河流域设计洪水分析报告》成果,该成果经海委审查后颁布使用至今。新镇站为大清河洪水控制断面。大清河流域控制站设计洪水成果,见表1。

表1 大清河流域控制站设计洪水成果

文安洼蓄滞洪区还涉及到平原沥水的蓄滞问题。大清河清南及清北平原的设计沥水,采用水利部天津水利水电勘测设计研究院1997年完成的《大清河永定河下游地区超标准洪水安排规划》中的分析成果。大清河流域洪沥总成果以第六堡为大清河流域总控制断面,30 d洪量系列计算采用山区洪水控制,平原采取相应的取样方法,将大清河流域30 d洪量与平原相应涝水相加,组成大清河全流域30 d洪量系列进行频率分析,得出大清河流域30 d洪量地区组成成果,见表2。

表2 大清河流域30 d洪量地区组成成果亿m3

3.2 设计洪水下垫面一致性修订成果

在水利部公益性行业科研项目“下垫面变化条件下设计洪水修订技术研究”中对下垫面变化对大清河主要控制站设计洪水影响进行了修订,提出下垫面影响一致性修订设计洪水成果。新镇及第六堡控制断面设计洪水修订成果与现有成果比较情况如下:50年一遇洪水标准时30 d洪量,新镇控制断面比现行成果下降14.2%,第六堡控制断面比现行成果下降25.1%;100年一遇洪水标准30 d洪量,新镇控制断面比现行成果下降11.7%,第六堡控制断面比现行成果下降20.8%。

4 蓄滞洪区洪水模拟模型

本次采用的数学模型为平面二维非恒定流水力学模型,已在南水北调中线设计及多项铁路、公路工程的防洪评价中采用。

描述二维非恒定流运动规律的是圣维南方程组,圣维南方程组为:

式中:h为柱面平均水深(m);H为水面高程(m);U、V分别为X、Y方向的平均流速(m/s);τxb、τyb分别为X、Y方向的床面剪切应力(N/m2);M、N分别为X、Y方向的单宽流量(m3/s)。

由式(1)—(3)推导,得出以下微分方程:

将上述微分方程省略对流项后,得到如下求解形式的差分方程组:

式中水深与流速(单宽流量)在空间和时间上呈交错分布,如图2所示。

具体计算步骤为:空间上,在(i+1/2,j+1/2)位置上计算h,在(i+1,j+1/2)、(i,j+1/2)位置上计算X方向上的M,U,在(i+1/2,j+1)、(i+1/2,j)处计算Y方向上的N,V。时间上,在n+1、n+3的时间层上计算h,在n、n+2时间层上计算M、N、U、V。确定初始及边界条件后,可由式(7)—(9)逐时段模拟洪水过程。

图2 差分格式示意

5 模拟范围及洪水数值模拟分析

由于大清河中下游的东淀、文安洼、贾口洼为联合调度运用洼地,文安洼和贾口洼的分洪运用对东淀流态、流势均有一定影响,所以需将上述三洼均纳入数值模拟范围。具体范围为:西部边界,清南地区为白洋淀千里堤,东淀区域为刘家铺村;南部边界为津保公路;北部边界为东淀北大堤;东部边界为南运河左堤、东淀内西河右堤。模拟区总面积3 228 km2,覆盖了任大公路以北的所有清南地区以及整个东淀、贾口洼蓄滞洪区。大清河三洼洪水数值模拟区地形剖分立体图,如图3所示。

图3 大清河三洼洪水数值模拟区地形剖分立体图

按大清河洪水调度方案,当东淀第六堡水位达到6.44 m且继续上涨威胁天津市区安全时,在保证河道充分泄洪的情况下,视洪水情况向文安洼或贾口洼分洪。当发生北支为主的洪水时,若文安洼尚未分洪运用,则向贾口洼分洪;当贾口洼八堡站水位达到5.94 m且继续上涨时,视洪水情况向文安洼分洪。

针对大清河流域主要控制站修订前、后的设计洪水成果,按照现状工程条件和现行调度运用方案,采用二维非恒定流模型对东淀、文安洼、贾口洼3个蓄滞洪区的蓄滞洪过程进行模拟分析。三洼蓄滞洪区典型位置最高洪水位成果,见表3。

6 结论

由于流域下垫面变化、气候变化等的影响,流域总体洪涝水明显减少。在相同的运用方式和条件下,经过模拟计算,50年一遇洪水时,东淀第六堡最高滞洪水位6.44 m、与现行设计洪水成果一致,文安洼大赵最高滞洪水位降低0.63 m,贾口洼锅底闸最高滞洪水位降低0.45 m;100年一遇洪水时,东淀第六堡最高滞洪水位降低0.10 m,文安洼大赵最高滞洪水位降低0.29 m,贾口洼锅底闸最高滞洪水位降低0.55 m。

TV122+.3;TV873

A

1004-7328(2015)01-0041-04

10.3969/j.issn.1004-7328.2015.01.015

2014-10-11

齐静(1979-),女,硕士,工程师,主要从事水利规划、水资源方面研究工作。

猜你喜欢
文安滞洪区下垫面
城市下垫面渗蓄性能量化模拟试验研究
大陆泽蓄滞洪区滏阳河老油召桥至老北定桥段工程地质问题概述
胖头泡蓄滞洪区的工程管理制度与职责探析
北京与成都城市下垫面闪电时空分布特征对比研究
雨姑娘快下凡
2018高文安设计新航向AIDG聚创国际新启航
将不动
香港室内设计之父高文安设计作品探析
潖江蓄滞洪区建设与管理初探
下垫面变化条件下的官厅山峡段流域洪水管理