山地小流域降雨及输沙过程的计算模型及其数值模拟.docVIP

山地小流域降雨及输沙过程的计算模型及其数值模拟 目录 TOC \o 1-9 \h \z \u 目录 1 正文 1 文1：山地小流域降雨及输沙过程的计算模型及其数值模拟 1 文2：山地小流域降雨及输沙过程的计算模型及其数值模拟 9 1 小流域单元划分及流量过程的计算 9 2 小流域输沙过程的计算 11 3 验证 12 4 结语 14 参考文摘引言： 14 原创性声明（模板） 14 文章致谢（模板） 15 正文 山地小流域降雨及输沙过程的计算模型及其数值模拟 文1：山地小流域降雨及输沙过程的计算模型及其数值模拟 开展小流域的治理及其环境生态恢复，是日前西部大开发， 水土保持 及山川秀美工程的主要任务。本文以山地小流域为研究对象，利用数理分析的方法，探讨降雨、汇流、输沙过程的数值计算关系，拟建立相应的系统计算模型，为进一步研究江河泥沙灾害的预测及其防治，打下初步基础。 1 小流域单元划分及流量过程的计算 图 1 小流域示意图 Sketch of Small Watehed 图 2 小流域简化图 Predigesting Sketch of Small Watehed 本文以流域内部的脊线及分水岭作为分界边线，将其划分成 n 个小单元，其面积从流域出口处开始往上游编号，它们依次为 F 1 、 F 2 、… F n ，如图 1 所示。为研究的方便，图 1 可进一步简化为图 2 的形式。 设两相邻单元出口位置 ( 含小流域出口断面 ) 的主沟长度为 L ij (i=0,2...n-1;j=1,2 …… n) ，次降雨过程是由多峰叠加组成的，其强度 I 服从半周期正弦函数分布，其与时间 t 的关系可表示为 (1) 式中 I mj 为降雨强度中第 j 个峰值； t 0j 为第 i 个峰值起始时间； T 0j 为第 j 个雨峰历时； l 为峰值个数土壤入渗强度 f 可分段表达为 [1] f=I(t) tt p f=mt -n +f c tp≤t≤t c f=f c tt c 式中 t p 为初始产流历时； t c 为达到稳渗率历时； f c 为稳渗率； m 、 n 分别为系数和指数。 山地小流域产流模式多为超渗产流，由于坡陡流急，可忽略降雨过程中的蒸发及小单元的填洼水量，则 t 时刻 i 单元的径流量 Q i (t ）为 当 t≤t p ,f=I(t) ， Q i (t)=0 (3) 当 t∈(t p ,t c ］ ,f=mt -n +f c (4) 当 tt c ,f=f c (5) 令小流域出口断面为 0 号断面。假设流域内每一小单元的水量在向 0 号断面汇集时，其相互影响忽略不计，则 t 时刻每一小单元水量汇聚到 0 断面的时间和数值分别是：当 t ≤ t p 时 Q 0 (t) ＝ (6) 当 tt p 时，出口流量有可能不为 0 ，对于 t 时刻的降雨，流至出口断面的流量分别是 Q 01 (t+Δt 01 )=Q 1 (t);Δt 01 =L 01 A 01 ／ Q 1 (t) (7) 式中 L 01 为 0 断面至 1 号小单元出口的距离； A 01 为 Q 01 对应的主沟道过水面积。同理 Q 02 (t+Δt 02 )=Q 1 (t+Δt 02 -Δt 01 ) ＋ Q 2 (t);Δt 02 ＝ L 02 A 02 /Q 2 (t) (8) Q 0j (t+Δt 0j )=Q 1 (t+Δt 0j -Δt 01 ) ＋ Q 2 (t+Δt 0j -Δt 02 ) ＋ ……Q j-1 (t+Δt 0j -Δt 0j-1 ) ＋ Q j (t);Δt 0j ＝ L 0j A 0j ／ Q j (t) (9) Q 0n (t+Δt 0n )=Q 1 (t+Δt 0n -Δt 01 ) ＋ Q 2 (t+Δt 0n -Δt 02 ) ＋ ……Q n-1 (t+Δt 0n -Δt 0n-1 ) ＋ Q n (t);Δt 0n =L 0n A 0n ／ Q n (t) (10) 式中 L 0j 为 0 断面至 j(j=2,3 …… n) 号小单元出口的距离； A 0j 为 Q 0j 对应的主沟道过水面积。出口断面的总流量应为同一时刻每一小单元单独产生的流量过程对应的出口流量之迭加 Q 0 (T)=Q 01 (T)+Q 02 (T)+……Q 0n (T) (11) (11) 式即为小流域出口断面流量变化过程关系，由此可模拟出其径流过程。 2 小流域输沙过程的 计算 坡面上的雨滴溅蚀及漫流侵蚀的泥沙，均向沟底汇集，流向小单元的出口，当此处的水流挟沙能力小于来沙量时，将有部分泥沙沉积；反之，水流将从沟底摄