《山地小流域降雨及输沙过程的计算模型及其数值模拟》由会员分享，可在线阅读，更多相关《山地小流域降雨及输沙过程的计算模型及其数值模拟（9页珍藏版）》请在金锄头文库上搜索。
1、1山地小流域降雨及输沙过程的计算模型 及其数值模拟摘要：本文依据山地小流域产流及汇流特征，采用数理分析方法，探讨了其降雨强度、出口流量及输沙率随时间变化过程的数学表达关系，建立了从降雨至径流输沙的系统模型，经编程计算，模拟值与某小流域降雨和汇流资料吻合较好。 关键词：小流域 降雨过程 输沙过程 非线性关系 开展小流域的治理及其环境生态恢复，是日前西部大开发，水土保持及山川秀美工程的主要任务。本文以山地小流域为研究对象，利用数理分析的方法，探讨降雨、汇流、输沙过程的数值计算关系，拟建立相应的系统计算模型，为进一步研究江河泥沙灾害的预测及其防治，打下初步基础。1 小流域单元划分及流量过程的计算图
2、1 小流域示意图Sketch of Small Watershed图 2 小流域简化图Predigesting Sketch of Small Watershed2本文以流域内部的脊线及分水岭作为分界边线，将其划分成 n 个小单元，其面积从流域出口处开始往上游编号，它们依次为 F1、F2、Fn，如图 1 所示。为研究的方便，图 1 可进一步简化为图 2 的形式。设两相邻单元出口位置(含小流域出口断面)的主沟长度为 Lij(i=0,2.n-1;j=1,2n)，次降雨过程是由多峰叠加组成的，其强度 I 服从半周期正弦函数分布，其与时间 t 的关系可表示为(1)式中 Imj 为降雨强度中第 j 个峰
3、值；t0j 为第 i 个峰值起始时间；T0j 为第 j 个雨峰历时；l 为峰值个数土壤入渗强度 f 可分段表达为1f=I(t) ttp(2.1)f=mt-n+fc tpttc(2.2)f=fc ttc(2.3)式中 tp 为初始产流历时；tc 为达到稳渗率历时；fc 为稳渗率；m、n 分别为系数和3指数。山地小流域产流模式多为超渗产流，由于坡陡流急，可忽略降雨过程中的蒸发及小单元的填洼水量，则 t 时刻 i 单元的径流量 Qi(t）为当 ttp,f=I(t)， Qi(t)=0(3)当 t(tp,tc,f=mt-n+fc(4)当 ttc,f=fc(5)令小流域出口断面为 0 号断面。假设流域内每
4、一小单元的水量在向 0 号断面汇集时，其相互影响忽略不计，则 t 时刻每一小单元水量汇聚到 0 断面的时间和数值分别是：当 ttp 时4Q0(t)0.0(6)当 ttp 时，出口流量有可能不为 0，对于 t 时刻的降雨，流至出口断面的流量分别是Q01(t+t01)=Q1(t);t01=L01A01Q1(t)(7)式中 L01 为 0 断面至 1 号小单元出口的距离；A01 为 Q01 对应的主沟道过水面积。同理Q02(t+t02)=Q1(t+t02-t01)Q2(t);t02L02A02/Q2(t)(8)Q0j(t+t0j)=Q1(t+t0j-t01)Q2(t+t0j-t02)Qj-1(t+t
5、0j-t0j-1)Qj(t); t0jL0jA0jQj(t)(9)Q0n(t+t0n)=Q1(t+t0n-t01)Q2(t+t0n-t02)Qn-1(t+t0n-t0n-1)Qn(t); t0n=L0nA0nQn(t)(10)式中 L0j 为 0 断面至 j(j=2,3n)号小单元出口的距离；A0j 为 Q0j 对应的主5沟道过水面积。出口断面的总流量应为同一时刻每一小单元单独产生的流量过程对应的出口流量之迭加Q0(T)=Q01(T)+Q02(T)+Q0n(T)(11)(11)式即为小流域出口断面流量变化过程关系，由此可模拟出其径流过程。 2 小流域输沙过程的计算坡面上的雨滴溅蚀及漫流侵蚀的泥
6、沙，均向沟底汇集，流向小单元的出口，当此处的水流挟沙能力小于来沙量时，将有部分泥沙沉积；反之，水流将从沟底摄取部分泥沙，输向下游，造成沟道冲刷。因此，可以认为，小单元出口处的水流含沙量接近于水流的挟沙能力 S，其值可由下式计算2S*4h0.5J/n2w(12)式中 h 为小单元出口位置处的水深：，J 为坡度；n 为糙率；w 为泥沙的沉速，利用上述流域主沟内水量迭加原理，可得小流域出口断面输沙量的迭代关系：当ttp 时，各小单元的沙量输至 0 断面的时间及数值分别是Qs01(t+t01)=Q1(t)S*1(t) t01=L01A01/Q1(t)6(13)Qs02(t+t02)=Q1(t+t02-
7、t01)S*1(t+t02-t01)Q2(t)S*2t）t02=L02A02/Q2t）(14)Qs0j(t+t0j)=Q1(t+t0j-t01)S*1(t+t0j-t01)+Q2(t+t0j-t02)S*2(t+t0j-t02) +Qj-1(t+t0j-t0j-1)S*j-1(t+t0j-t0j-1)+Qj(t)S*j(t);t0j=L0jA0j/Qj(t)(15)Qs0n(t+t0n)Q1(t+t0n-t01)S*1(t+t0n-t01)+Q2(t+t0n-t02)S*2(t+t0n- t02）+Qn-1(t+t0n-t0n-1)S*n-1(t+t0n-t0n-1)+Qn(t)S*n(t);
8、 t0n=L0nA0n/Qn(t)(16)式中 L0j 为 0 断面至 j(j=2,3n)号小单元出口的距离；A0j 为 Qoj 对应的主沟道过水面积。小流域出口断面的输沙过程，应为同一时刻每一小单元产生的输沙过程对应的输沙量之迭加Qs0(T)=Qso1(T)Qs02(T)+Qson(T) (17)小流域出口断面的总输沙量为(18) 3 验证为了检验上述计算关系是否合理，笔者利用湖北省某径流试验场的观测资料，进行了验证。该径流场集水面积 31540m2，长度 432m，平均宽度 73m，平均坡度7607；坡向朝西北，坡面上淤积物及基岩占 65.7%和 27.1%，其余为坡残积物。流域内土质为中
9、砂土、重砾质轻壤土和重砾质砂壤土，其比例分别为 29.7%、34.1%和 9.1%。依据地形条件，将其划分成 7 个小单元；又由 1966 年 9 月降雨观测资料，将整个降雨强度过程描述成 4 个单峰小过程的迭加。具体计算数据如下表所示2表 1 小流域汇流产沙系统模型计算的有关数据表The data related to the calculation of system model in a small watershed流域小单元面积(m2)降雨峰值(mm/min)渗透参数F1F2F3F48F5F6F7Im1Im2Im3Im4mnfc3046321840435721730739454259
10、1.5391.821.00.43.950.50.4利用上述推导的基本关系及流域计算数据，笔者进行了编程计算。降雨过程如图 3 所示，从图中可以看出，模拟曲线与实际降雨曲线基本吻合。在此基础上，计算了小流域出口断面的输沙率过程，其值绘于图 4 中。图 3 降雨强度过程曲线Intensity process of a rainfall 图 4 小流域出口断面输沙率过程Sediment load process of outlet section in a watershed由于野外观测数据无法绘成连续曲线，故仅给出了数据点；从图 4 中可以看出，模拟曲线与观测点据吻合较好。对比降雨及输沙过程，可明
11、显看出，小流域出口断10面输沙率过程，迟后于降雨过程；降雨出现两次较大的峰值，则输沙过程也出现两次峰值；当降雨强度小于土壤入渗强度，地表就不产生径流，但下渗的水体，会在重力作用下，向低处流动，导致下游土层含水量增高，甚至超出入渗强度，冒出地表，此现象常造成短历时降雨后，流域出口仍有小量的水体流出，使汇流过程拉长。4 结语本文依据山地小流域降雨、产流、汇流、输沙的基本规律，建立了降雨强度过程模拟关系(1)；小流域出口断面的流量过程的计算关系式(7)(11)；和小流域出口断面的输沙率过程的计算关系式(13)(18)。由于上述公式中参数较多，计算时尚需根据当地的具体资料率定相关系数。参 考 文 献1 蔡强国，王贵平，陈永宗。黄土高原小流域侵蚀产沙过程与模拟。科学出版社，1998.2 钱宁，万兆惠。泥沙运动力学。科学出版社，北京，1986.3 朱蕾。潮湿山区小流域产流规律的分析与模拟。四川联合大学硕士学位论文，1998.4.
