多点入流汇流计算法在北京城市洪水计算中的应用与研究.docx

多点入流汇流计算法在北京城市洪水计算中的应用与研究文立道范本贤徐晓峰韦明杰(北京城市规划设计研究院)摘要根据北京城市化水平(下垫面情况及排水条件)、城市汇流特点(雨水管道及明渠排水),通过用多点入流汇流计算法模拟实测洪水过程线,分析了北京城市河道的汇流特点与规律,初步提出了一种适合于北京城市洪水计算的方法,为城市防洪排水工程规划、设计提供了依据。关键词汇流计算产流分析规划流量城市排水系统北京市护暗沟和南、北扩城河等。通惠河乐家花园站以上水系及站网分布详见图1。1前言随着北京市经济建设的发展,人口不断向城市地区集中,城市规模在不断扩大。大规模地建造房屋、铺砌道路、衬砌河道、绿化环境等,迅速地改变着原有下垫面条件;城市排水河道的大规模整治与改建、雨水管道的系统化,改变了原有降雨径流的形成规律。由于城市暴雨洪水的产汇流条件极其复杂,北京市以往采用的等流时线法、推理公式法、排涝模数法等,干支流或上下游断面之间洪水不便组合,使得计算成果都不很理想。为确保城市防洪安全,并为城市防洪排水工程规划提供依据,需要进行城市洪水计算方法的研究。本文选择了水文资料较多、建设区面积较大的市区主要排水河道通惠河乐家花园站以上流域为代表。从分析实测雨洪资料入手,研究城市化对洪水特征值的影响;用多点入流汇流计算法模拟实测洪水过程线,研究该方法的适用性,初步分析出了不同城市化水平的汇流规律,并修正了设计暴雨;最后计算出了通惠河乐家花园站以上河道的规划流量。通惠河自西向东汇入北运河,乐家花园站以上流域面积约190km2,上游支流有南旱河、永引渠、京引渠、长河、前三门暗沟、东护暗沟、西图1通惠河乐家花园站以上水系及站网分布示意图2多点入流汇流计算法2.1多维线性定常汇流系统多点入流汇流计算法为多维线性定常汇流系统,它考虑了净雨空间分布的不均匀性,运用河槽汇流曲线进行以等流时块为单元区的流域汇流计算或以支流流域为单元区的河网汇流演算,以及考虑区间支流的河道洪水演算等问题。虽然情况各异,但具有共同的特点:一是入22流非单一;二是出流由各入流通过各自相应流程的河槽汇流曲线,转换为单元出流过程,并叠加而得到出口断面的出流过程。为便于研究,将这种汇流问题归纳为多点入流汇流问题,其系统框图如图2所示。当输入(净雨或入流流量时,其表达式为:CknF·Rki-j+1∑∑Pk,j·Qi=(3)△tk=1j=1(1≤i-j+1≤r)式中:Ck为第k等流时面积的汇流曲线时段数;F与R为等流时面积及其时段净雨深;n为等流时块数;P为长办汇流曲线系数。2.2长办汇流曲线长办汇流曲线是于1962年,在г.п.加里宁河槽汇流曲线概念基础上,研究提出的一种入流为矩形的河槽汇流曲线,其含义是:一历时为tk,流量为1m3/s的矩形入流,流经任意n个河段后产生的出流过程,其计算公式如下:n-1im1-e-M·∑i!i(0≤m≤mk)图2多点入流系统框图(4)=0Pn=n-1n-j-i过程)为阶梯型时,则有:rjMjimk∑j!·e-M-mk·∑i!)(1-e(mk≤m)j=0i=0mk;m=t/τ;mk=tk/τ;n与τ为Qj(t)=∑qji·p(t-(1)式中:M=m-(i-1)△t)njj=1流经河段数与汇流时间;tk为入流历时;Pn系∑Qj(t),故其多维线系统总输出为Qj(t)=河段数为n的长办汇流曲线系数。j=1性定常汇流系统方程为:3产流分析rjnQj(t)=∑∑qji·pj(t-(i-1)△t)(2)3.1下垫面的分析采用航空遥感照片进行下垫面分析,将通惠河乐家花园站以上流域分为楼房区、老平房区、新平房区、学校和体育场馆区、工厂区、农田绿地区、河湖区等7类地区。对各类地区的不透水面积,采用典型块法分析计算,在比例尺为1:2000及1:6000的典型块航空遥感照片上分析和辅以现场调查,进而计算出1959年、1983年两个不同时期乐家花园站以上各地区的下垫面组成情况。分析结果是:1959年全流域不透水面积占33%,其中城区站(指玉渊潭出口站及松林闸站至乐家花园的区间流域)不透水面占61%;1983年全流域不透水面积占47%,其中城区站占77%。1983年各区间流域的不透水面积百分比均比1959年有所增加,反映了城市建设的发展对下垫面的影响。23j=1i=1(t-(i-1)△t)≥0式中:t与△t分别为时程与阶梯型输入的单位时段、汇流曲线的入流历时或计算时段;n为子系统数(等流时块数);q为阶梯形输入(等流时块上的净雨);p为子系统的传递函数(汇流曲线);r为阶梯形输入情况的入流时段数;Qj与Q为第j子系统与系统的输出。从上述系统方程可以看,每个子系统的汇流曲线不同,这就反映了此方法用于河道汇流计算时,考虑了干流上游断面和区间各支流入流点至下游断面的距离不等应受的不同调节作用。用于流域或河网汇