chu修正的四种径流计算公式的对比分析.docxVIP

chu修正的四种径流计算公式的对比分析 混合产流模型 降水入渗产流是引起土壤侵蚀和生态环境影响的原因之一。许多研究人员对土壤入渗和生产过程进行了研究。1911年Green-Ampt以毛管理论为基础,提出具有相同初始含水量的均质土壤的下渗方程。1973年Mein 和Larson在Green-Ampt公式的基础上确定出了稳定降雨条件下,地表积水时入渗方程(GAML),该方程为: f=dFdtΚ(1+SΜ/F)(1)f=dFdtK(1+SM/F)(1) 式中:f为下渗能力(mm/min);K为土壤饱和导水率(mm/min);SM为有效土壤水吸力(mm);t为历时(min);F为累积下渗量(mm)。 Chu考虑了不稳定降雨条件下的产流特性,利用GAML方程,提出了不稳定降雨条件下的产流计算方法。Chu认为在非稳定降雨情况下,地表时而有积水发生,时而又无积水发生,因此一场降雨有多个积水点。当将降雨过程看作是许多时段降雨的累积,则小时段内可近似认为降雨强度是恒定的。在此假定条件下,可推求出时段净雨。Philip1975年在Rechilds方程的基础上,运用Bolman变换,并结合一定的边界条件和初始条件,得到了方程的级数解,其二项式入渗方程为: fp=12S?t-0.5+A(2)fp=12S?t?0.5+A(2) 式中:fp为下渗能力(mm/min);t为历时(min);S为吸渗率(mm/min-0.5);A为稳渗率(mm/min)。 Horton认为当降雨持续进行时,下渗率逐渐减小。下渗过程是一个消退过程,消退速率与剩余量成正比。即消退速率df/dt与给定时间t内有待消退的量(从该时刻的下渗率f到稳渗率fc的改变量)成正比。根据以上假定得到下渗公式为: fp=fc+(f0-fc)e-kt(3)fp=fc+(f0?fc)e?kt(3) 式中:f0为初始下渗率(mm/min);fc为稳定下渗率(mm/min);k为衰减指数。 美国土壤保持研究局为了计算无降雨过程记录地区的径流量,研制了径流曲线数(Curve Number)法。该计算方法中只有一个反映流域特征的一个综合参数CN,它与流域土壤类型和土地利用现状等有关。该产流计算法具有结构简单、计算方便的特点,同时能考虑土地利用变化等对产流的影响。因此,它在美国及其它一些国家已得到广泛应用。美国许多土壤侵蚀模型如EPIC、CREAMS、AGNPS等模型中都运用该法进行产流量的计算。其计算公式为: Q=(Ρ-0.2S)2Ρ+0.8S(Ρ0.2S)(4a)Q=0(Ρ≤0.2S)(4b)S=25?400/CΝ-254(5)Q=(P?0.2S)2P+0.8S(P0.2S)(4a)Q=0(P≤0.2S)(4b)S=25?400/CN?254(5) 式中:Q为径流量(mm);P为降雨量(mm);S为吸水参数;CN为径流曲线数,与植被覆盖、水文土壤组以及前期影响雨量等因素有关,取值范围为0～100。 在应用较为广泛的Green-Ampt入渗方程、Philip入渗方程和Horton入渗方程的基础上,国内学者针对各地的实际情况,又提出了若干入渗产流模型。如包为民等根据蓄满产流和超渗产流的产流机制,结合流域产流的实际情况,提出了混合产流模型。郭瑛利用Horton入渗方程,结合流域入渗分配曲线,提出了非饱和产流的产流计算方法。郝振纯在饱和入渗理论的基础上,考虑新入渗雨水驱替旧土壤水分的作用,根据力学原理提出了一个简明入渗模型。但是到目前为止,对北京山区的入渗模型研究还比较少,本文用北京山区坡面径流小区的降雨径流资料,对Chu修正的GAML入渗方程,Phillip入渗方程,Horton入渗方程和SCS径流曲线数法进行比较分析,找出适合北京山区的径流计算方法。本研究结果可用于水土保持规划和评价以及水土资源管理。 1 学习方法 1.1 径流、土壤条件 选取了房山蒲洼两个坡耕地坡面径流小区(坡度分别为11°和23°,作物为玉米)1992～1995年30场实测降雨径流资料和密云石匣休闲地坡面径流小区(坡度为14.6°)1993～1998年13场实测降雨径流资料。蒲洼坡面径流小区位于北京市西南石质山区,地理坐标为115°35′～115°46′E, 39°43′～39°49′N,小区于1992年建成。土壤为棕壤,土壤质地粗,土层薄。小区长10 m,宽5 m。密云石匣休闲地径流试验小区位于北京市密云水库东北,地理坐标为117°01′～117°07′E,42°32′～42°38′N,小区于1993年建成,土壤为粗骨褐土,土层薄。小区坡长为10 m,坡宽为5 m。 1.2 illip入渗方程及horton入渗方程 本文选用常用的入渗方程,即:CHU修正的GAML方程,Phillip入渗方程,Horton入渗方程进行计较。同时考虑到无资