徐 磊, 林茂林, 李 聪, 张亚倩
(1.无锡市水利设计研究院有限公司, 江苏 无锡 214000; 2.无锡太湖学院商学院, 江苏 无锡 214100)
苏南农村经济发达,居民富足,地少人多,自然村落密集,千亩以上圩区均有居民居住区和工业企业坐落,无概念上的纯农业圩区。江阴市共有大小圩区96个,圩区面积1.8万hm2,共建有大小圩区96个,其中千亩以上圩区有53个,万亩以上圩区有8个。圩内河网纵横,水面率超过8%,圩区内排涝泵站规模不大,单泵流量很少超过4 m3/s,且具有以下社会经济特点:
(1)工业企业多,且有一定的规模。苏南是乡镇企业发展地,改革开放以来,各镇、村均兴办了大量的工业企业。经过40多年的发展,催生了一大批大型企业乃至龙头企业。江阴市月城镇马甲圩为江阴市的一座万亩圩区,8.19 km2的圩区内有30多家企业,年工业生产总值超过20亿元,一旦受淹,经济损失巨大。
(2)经济作物多,传统农业种植少。随着农业产业结构的调整及土地流转,农业规模经营发展加快,一大批企业从工业转型为农业,高效农业、设施农业发展迅速,高经济附加值的水果及有机蔬菜的种植面积已超过了水稻等传统农业。
(3)居民生活条件好,但住房抗灾能力不高。2020年江阴市农民人均可支配收入38 225元,农民收入的连年提高,必然使包括住房在内的固定资产增加。但由于农村向城镇集中的政策,从20世纪90年代开始,苏南一些地区已控制农村翻建房屋,农民居住房屋多为20世纪80年代建造,房屋抗灾能力普遍较低。
(4)各级政府对防洪排涝等抗灾工程设施建设十分重视。由于经济发展,政府财政收入增加,对水利等基础设施投入也逐年增加。农民由于受淹后损失较大,质疑政府不作为,极易产生社会不稳定因素。不受涝、抗大灾的意识在各级政府中形成共识。苏南农村混合型圩区的排涝能力应达到短时受涝或不受涝的要求,这也是各级政府及群众的意愿所在。
《江阴市水系规划》(2011—2030年)提出,农业圩区执行20年一遇24 h暴雨雨后1日排除标准。圩区排涝流量的计算一般采用平均排除法。计算式如下。
(1)
其中:R水田=P-h田蓄-E
R旱=αP
R水=P-h水蓄-E
式中:Q为设计排涝流量;F水田为圩区水田面积;t为规定的排涝时间;P为设计暴雨量;h田蓄、h水蓄为水田和水面的滞蓄水深;E为历时为t的水田田间腾发量;α为径流系数
此方法考虑了各种种植结构及水塘调蓄等综合情况,目前普遍采用此法计算排涝流量。但对于调蓄能力较差、面积较大的混合型圩区,按此法计算的流量偏小。在规划圩区排涝流量时,往往根据硬化区域大小,对按平均排除法计算出的排涝流量进行适当增加,随意性较大,缺乏科学合理的依据。
苏南混合型圩区经济总量大,受淹后损失大,规划排涝流量应以满足一定频率暴雨情况不受淹为原则。考虑以河道允许最高水位作为控制条件,计算设计排涝流量。
河道控制水位涉及调蓄能力的计算,控制水位越低,排涝流量越高,应选择科学合理、确保区域不受淹的河道水位作为控制水位。一般以地面最远或最低处作为圩区的典型点,控制河道水位典型点高程以下0.2 m。
江阴市圩区执行的排涝标准为20年一遇24 h暴雨1日排出,因此,设计暴雨历时取24 h,以1 h、6 h为控制时段进行雨量分配。设计暴雨可采用实测暴雨资料。若无实测暴雨资料,可采用《江苏省洪水暴雨图集》。
圩区硬化区域与经济作物区由于不能受淹,在计算设计净雨过程中,采用后损1 mm/h计算。对于传统农业区域可受淹1 d以上区域,采用次降雨径流相关法计算。
采用瞬时单位线法进行汇流计算。瞬时单位线法是依据水量平衡方程和线性蓄槽方程,经过连续演算和数学处理导出的数学方程式:
(2)
式中:Γ为伽马函数;n为水库个数或调节次数的参数;K为水库滞时,相当于流域汇流的时间参数。
参数m1与m2可以通过数学方法进行优选。考虑苏南圩区缺乏比降资料,采用m1=2.25F0.38(F为圩区的面积)。参数m2比较稳定,平原区取0.5。
由m1及m2可以得到区域2 h单位线,乘以时段总径流量,并错开相加,可以得到设计洪水过程线。
假定河道的初始水位一般为预降后的水位,根据水量平衡计算出抽排后水位计算式。
(3)
式中:V1、V2分别为计算时段始末的河道蓄水量;Q1、Q2分别为计算时段始末的入流流量;q为假设的圩区排涝模数。
控制各时段的最高水位不超过控制水位,可计算出设计排涝流量。
本文选取江阴市月城镇马甲圩作为计算实例。马甲圩面积8.19 km2,圩区内乡镇企业多,农民大量种植经济作物,为典型的混合型圩区。
设计暴雨采用图集法,由《江苏省暴雨洪水图集》查得月城片区暴雨统计资料,时段分为1 h、6 h及24 h,结果如表1。
表1 24 h暴雨统计参数
表2 设计暴雨成果表
以1 h和6 h为控制时段,计算出10年一遇及20年一遇1 d暴雨时程分配。
马甲圩区水果种植面积约333.33 hm2,由于不能受淹,视作建成区,则生态区面积约为32%。生态区净雨采用降雨径流相关图法,可采用下式计算:
Pα=αImax=0.65×95=61.8 mm
(4)
式中:R为净雨深,P为设计雨量,Cp、Ci为产流参数。
建成区净雨按每小时扣1 mm来计算,计算出设计净雨深。
由瞬时单位线参数m1=2.25×F0.38=5.003,得到该区2 h单位线,将时段净雨利用公式
(5)
得出时段总径流量,再乘以时段单位线可得到20年一遇洪水过程线。
根据圩区的平均地面高程确定河道调蓄水深为0.5m,将产汇流作为河道的输入流量,假定排涝模数进行推演计算,得到河道的蓄水量以及相应的水位,将计算得到的河道最高水位与控制水位进行比较,河道最高水位与控制水位相等时的排涝模数就是所求的排涝模数。用式(6)进行调节计算。
经计算,马甲圩圩区需要的排涝流量为15.6 m3/s。本计算是针对整个圩区均匀排涝的计算,可作为排涝规划的参考。具体排涝设计时,由于分区排涝不同,排涝河道标准不一,应分别对分区排涝进行计算。
以达到河道控制水位为标准而计算出的排涝流量能确保一定频率标准下暴雨该地区不受涝,但该方法未考虑排水河道中水位差因素,因此较适合排水河道不长的中小圩区,苏南混合型圩区均为这类圩区。对于有多座排涝泵站的圩区,可以按每座排涝站的集水面积单独进行计算,确保各泵站设计流量的科学合理。