藉由Philip入渗公式推求有效降雨进行运动波模式特性法模拟迳流.pdfVIP

海峡雨岸山地炎害舆瑕境保育研究．第三卷 2002，昆明 藉由Philip入渗公式推求有效降雨逛行迤勤波 模式特性法模凝蓬流 隙榭群’、巫仲明2、黄美君3、何智武4 【1中典大掣水土保持早系教授萎中；2中典大挈水土保持早系博士生塞中 3桃圜褓政府壤保局技士桃圆；4中莘防挺单宙祭警理事畏矗中) 摘要本研究乃於均匀降雨之保件下，藉由Philip入渗公式推求有效降雨，韭以特性法演算 逼勤波模式，推蹲出特性栋上之水深及流量方程式，以此模凝漫地流流况，分析水剖面曲觫及流 量塍腺之解析解，韭且不以梭雅的有限差分形式，而是利用简罩的数值方法推尊由遗界憬件横成 之流量膣糠，针算屡於牛解析解的部份。由於本研究中引入陪梯函数表示违檩方程式之侧流瑁， 因此降雨停止後之遇水段不须再假彀初始水位重新演算。 明键嗣逼勤波理徐漫地流特性法Philip入渗模式 l前雷 最常被臆用於求解逼勤波模式的方法有二，其一是有限差分法，男一是特性法。有限 差分法有敷值援散(numeficaldiffusion)的现象，故可能燕法碓寅榆定出最佳之漫地流糙 度、漫地流畏度等遥些集水匿参数，模式桧定踌所割定之辋格，也因模撮之降雨强度不同 而须有所差具，而且凿正常水深之福禄敷燮小峙，属了逵到鞍佳的模凝秸果，初始辋格畴 and 同之筒距须明颢的燮／j、(Woolhiser Liggett，1967)。若所需的值非落在辋格黠上，则必须 再利用其他方法逛行内差种算。因此，搁格的大小常常影警模式的辈磋性及方便性。有箍 於此，本研究掇探用特性法逛行蓬勤波模式的解析计算，特性法除了没有数值方法中敷值 接散的缺黠，其精磋性及方便性亦均较敷值方法龠佳。由於特性法之敷挈推尊技巧比较繁 愎，馑在侧流项有较多限制，但也因此特性法之演算遇程有明碓的物理意羲。特性法臆用 於求解逼勤波理输之解析解，主要在於筒罩畿何地形上规则的降雨和入渗率的案例分析， Henderson and and and CundyTento(1985)以及Luce andWu 稼之颞式方程式，游出了秸合Philip入渗方程式的逗勤波理输之解析解法。Chen (1997)更将逼勤波模式秸合Phillips斜坡入渗公式湮用特性法之补算方式，臆用於地下水次 遥流之演算。 本文主要是虑用Philip入渗公式模凝地表入渗，於有限降雨延睛之均匀降雨绦件下计 算有效降雨，以此属违槽方程式中之侧流项(1ateral inflow)，利用莲勤波特性法分析罩一坡 and 面上漫地流。本研究所探用之方法是CundyTento(1985)的座檩褥换方式：fⅡ用特性波 之遗界傈件舆初始傈件爱敷r及逼勤波波傅畸同S惰成之S—r座檩平面，推尊特性禄上 9 之水深方程h(s，r)式及流量方程式q(s，r)。再以特性波波傅距雕建立S—r平面舆X--t平 面同之座檩搏换方法t臆用此座檩转换阴保将S—r平面上之水深及流量映射至x一，平面 andTento 上，分别求得漫地流水剖面曲腺舆流量塍橡。本研究之退水段演算有别於Cundy 之方法t富降雨停止之後，不再假彀一水面曲腺喾成退水段之初始保件，而是引入陪梯函 敷(unit step 法。 2逼勤波理萧 本研究中之逼勤波模式是屑於一稚漫地流模式，其控制方程式鹞： 皂+丝：吼d扩x t“ r11 s。=s，