瞬时单位线法求雨水管网系统入流流量过程线的数值计算方法_百度文库.DOCVIP

第30卷　第5期哈　尔　滨　建　筑　大　学　学　报Vol.30No.5　　　　　瞬时单位线法求雨水管网系统入流 流量过程线的数值计算方法 周玉文　赵洪宾　　李玉华 (给水排水系统研究室)　(给水排水教研室) 　　摘　要　用瞬时单位线法求S曲线,的计算方法,并给出了计算实例。 关键词　入流流量过程线;雨水管网系统;;分类号　TU992 线其中,、数学推导严密、参数可通过统计计算求得便于。本文旨在介绍用瞬时单位线法求城市雨水管网。 1　基本情况 1.1　瞬时单位线法的基本概念 在特定流域上,单位时段内均匀分布的单位净雨所形成的流域出口地表径流过程线为单位线。当净雨历时趋向于无限小时所求得的单位线叫瞬时单位线,通常用u(0,t)或u(t)表示。 u(0,t)=()N-1e- K#(N)Kk(1) 式中:#(N)为以N为参数的伽马函数;N为统计参数,相当于线性水库数或调节次数;K为统计参数,相当于流域汇流时间的参数,(min)。 112　基本假设 采用单位线法或瞬时单位线法计算流量过程线均假设流域为线性系统,符合倍比迭加原理。即 (1)单位时段内净雨深不同,它所形成的地表径流过程的总历时相同; (2)如果单位时段内净雨深不是一个单位,而是n个单位,那么它形成的地表径流过程线的流量值为单位线的n倍; (3)如果净雨历时不是一个时段,而是m个时段,那么各时段净雨量所形成的流量过程线之间互不干扰,出口断面的流量为各时段净雨量所形成的流量之和。 收稿日期:1996-11-25 周玉文　男　副教授 哈尔滨建筑大学市政环境工程学院(150008) 国家自然科学基金项目 42哈　尔　滨　建　筑　大　学　学　报第30卷113　城市雨水管网系统地表径流的瞬时单位线参数 根据北京市百万庄小区雨水管网系统降雨径流观测资料求得瞬时单位线统计参数为: 非铺砌地表　　非铺砌地表　　　　　　　　　　　　　　　(4) 　K=14115F-0102i0101(5) 式中:F为汇水面积,104m2;i为净雨强度,mm min; 　　以上公式适用于F≤115×103m2,　　综合地表NK(min)　N=2113F0138i-　　　　　　　　　　　　　　　(2)　K=10177F-0102i-0106　　　　　　　　　　　　　　(3)　N=1177F0138i6) 2,0=010015～0105不透水地表占50%左右的城市雨水2　数值计算方法 根据瞬时单位线理论,在确定了流域瞬时单位线参数后,即可以用数值方法计算S曲线、时段单位线和雨水口流域流量过程线。进行数值计算时首先应根据资料的精度条件和对计算结果的精度要求确定计算时间步长?t,根据城市雨水管网系统的特点,?t可以取1,2,5,10min。 211　S曲线 S(t)曲线为瞬时单位线的积分 S(t)=u(01t)dt∫0t =0K()N-1e-#NKt Kd()K =#(N)K0(N-1-)eKt Kd()K(7) 　　#(N)的值可以查表或用计算机求近似值,常用的方法有斯特令(Stirling)渐近展开公式或台劳展开式。Stirling公式: 1 2#(N)=NN-1e-N(2ΠN){1+++…}(8)2-3-12N288N51840N2488320N4 　　台劳展开式: 第5期周玉文等:瞬时单位线法求雨水管网系统入流流量过程线的数值计算方法43 (9)#(N)={A0+A1N+A2N2+…} [(N+1)N] 　　其中N≤1,Ai(i=0,1,2…)为系数,详见计算程序,当取前11项时,可以保证前6位有效数字准确。 积分式∫(K)0t KN-1-t Ked()为N的不完全伽马函数,可以采用34节点高斯积分K 公式求值。对于每一个ti=i?t,有: S(ti)=1-2#(N)∑N-1-wj[+]e22j=134[(40-t K)U40+t K+]22(10) 式中:Wj,Uj分别为高斯积分系数,详见计算程序。 用以上方法计算的S曲线其精度均高于现有S(t),,计算速度较快。 212　时段单位线 可以采用S,u(?t,t?t[(t)-S(t-?t)](11) 则 u(?t,ti)=?t[S(ti)-S(ti-1)](12) 其中:S(t0)=0 213　雨水管网系统入流流量过程线 雨水管网系统入流流量过程线就是各个雨水口流域的出流过程线。根据单位线的三条假设,可以用下式计算雨水口流域流量过程线。 Q(t)=∫u(?t,Σ)idΣ =∫[S(Σ)-S(Σ-0t-ΣΣmΣm 0?Σ)]?tdΣ(13) 则 m QK,j=∑H i=1j-i+1u(?t,ti)FK 60?t u(?t,ti)=m 60?t∑H i=1j-i+1(14) 式中:Qk,j为第K个雨水口第j个时段的流量值,L s; Hj-i+1