黄土坡面发育的细沟水动力学特征的研究docWord格式文档下载.docx

1、2 细沟水力要素的观测计算细沟侵蚀是径流汇集后发生的侵蚀过程，其水流运动方式已经有别于坡面漫流而类似于河流，只是坡面细沟流的沟床比降较大，径流深较小，及其形态演化较为迅速而已。因此，在目前还没有成熟的坡面流理论之时，可以借鉴河流动力学的原理和方法，用其相应的有关公式来计算研究细沟的水动力学特征。径流深、过水断面宽度、径流平均流速、雷诺数、弗氏数以及阻力系数等水力要素是反映水流动力学特征的主要指标，可以用相应的水力学公式对细沟流的上述水力要素进行计算或测定。Darcy weisbach阻力系数由公式计算f=8gRJ/U2（1）式中 f是Darcy weisbach糙度系数；g是重力加速度；R为水

2、力半径，R=A/P,A为过水断面面积，P为湿周。J为水面能波，U为水流平均流速。无量纲参数雷诺数及弗氏数也可用来表征细沟流紊动性程度和流态Re=Uh/v（2）Fr=U（3）Re为雷诺数，Fr为水流弗罗德数，U为细沟断面平均流速，h为断面平均径流深，v是水动粘性系数。g是重力加速度。曼宁糙率系数用曼宁公式计算U=（1/n）R2/3J1/2（4）n=（R2/3J1/2）/U（5）n曼宁糙率系数。R为水力半径，在宽深比较大时可用径流深代替。J为水流能坡，在细沟底坡较大时可近似的用地表坡度代替。U为径流的平均速度。过水断面宽度是根据量测的细沟横断面形态数据而建立的细沟断面上不同深度与相应深度处的断面宽

3、度的关系，经验关系式（6）计算过水断面宽度B=6.148h0.431（6）式中 B为过水断面宽度，cm，h为径流深，cm。径流平均流速不仅是主要的水力学参数，而且也是计算其他许多水力要素所必须知道的参数。试验中用染色法测定了径流平均流速。其方法是多次测定径流流经固定细沟区间的时间，用式U=1/mL/ti计算平均径流流速。m为测验定次数，L为测验定段的细沟长短度，ti为径流流经测定段所需时间。水深是反映水力特征的重要因子。试验中用淤积地形测定仪测定计算了细沟中平均径流深。其方法是多次测定不同细沟断面处的径流深，用式h=1/mhi计算平均径流深。 m为测定次数，hi为每次测定的径流深。3 试验程序

4、及方法据前人的研究3，坡面细沟侵蚀主要受制于细沟中水流的特征及土壤性质，受雨滴打击的影响很小。一场暴雨过程中，自然坡面上发育的细沟，其深宽因分布的部位、土质条件以及降雨强弱的不同，其规模在尺度上变化很大，但其宽度一般为几厘米到几十厘米不等，深度大多也不超过二十厘米（黄土坡面上梨底层的深度）。因此，细沟侵蚀过程可以通过水流冲刷槽试验来模拟，细沟流的水力学特征及其之间的相互关系也就可以通过水流冲刷试验来研究。本试验在宽为50cm，长500cm，深70cm的可调坡钢制冲刷槽内进行。在装填实验土之前，先在实验槽中铺填30cm厚的天然沙，以保持实验土的透水状况接近天然坡面。然后，现场采集的土样经过筛选后

5、，把通过2cm孔径的土样填入试验槽，填土过程中边填边用力压实，然后再将表层20cm左右翻松以模拟耕作层，填土完成后试验槽土壤表层的干容重为1.17g/cm3。供水设备采用定水头控制流量，从试验槽上端按设计要求加注不同的流量。为了保证每次试验的初始条件一致，试验开始前先均匀地在试验土表面撒水，撒水量控制在能使土壤表面充分饱和，但又未发生产流的程度。试验开始后，仔细观察径流冲刷过程，把握细沟形成时刻，一但细沟形成，使用染色法测定细沟中径流的平均流速，用清华大学研制的淤积地形测定仪测量水深的变化。试验过程分为三阶段，每一阶段10分钟，每个阶段之后，在试验槽的不同部位用测针随机地测定多处细沟断面的形态

6、，整个实验持续30分钟。试验流量变化于40500ml/s之间， 每个流量值连续重复试验二次。坡度采用6、10、12和15共四级变化。试验用土表1 试验土壤的粒径组成Particle size distribution and median diameter of tested soil粒径（mm）0.250.050.050.010.010.0050.0050.0010.001百分比（%）3.363.214.312.86.4中数粒径0.0185mm为采自北京官厅水库的黄土母质，土样中数粒径为0.0185mm,其粒径组成如表1所示。4 结果分析4.1 坡面股流的流态分析坡面漫流汇集成股流以后，由于

7、流速的增大，径流侵蚀力大增，其结果在坡面形成细沟。同时，径流的流态发生了很大的变化。表2是根据试验中实测的细沟流的水力要素值计算的水流弗罗德数，表明在试验的坡度及流量范围内，细沟流的弗罗德数都大于1，说明坡面径流汇集成股流后，其流态已属于急流。而且有坡度和流量越大，弗罗德数也越大；在流量相同时，坡度越大弗罗德数也越大的变化趋势。弗罗德数随流量的变化较为复杂，在坡度为6的坡面上，弗罗德数随流量的变化基本上保持不变。而当坡度较大时，弗罗德数随流量的增加有变大的趋势，同时坡度越陡，这种增大的趋势越明显。径流弗罗德数的这种变化是由于流量及坡度的大小对水深和流速影响作用的具体体现。径流流速随流量和坡度的

8、增加都会不断增大，而水深随流量的增加虽然也会变大，但随坡度的增加又会减小。而且，坡度越大，流速增大程度越大，径流深减小变化也越大。因此，径流弗罗德数就有随流量和坡度的增大而变化，且坡度越陡，这种增大越明显的趋势。坡面流汇集成股流后，径流雷诺数随流量及坡度的不同，其大小变化较大。表3是根据试验中实测的坡面细沟水力要素值计算的径流雷诺数大小，表明在试验的坡度及流量范围内，坡面股流的雷诺数变化于12007700之间。而且，雷诺数的大小也是随流量和坡度的不同有所变化；在流量不变坡度增大时，由于径流含沙量变化的复杂性，雷诺数变化的规律性不强，且变化也不显著。而在坡度相同流量增大时，雷诺数则有迅速增大的趋

9、势。说明坡面股流的雷诺数变化，受径流量的影响较大。若按照明渠水流的划分标准4，上述雷诺数的变化范围说明坡面股流基本上属于层流和紊流的过渡流范围，在流量较大时，其流态也进入紊流状态。但考虑到坡面流受边壁扰动大等实际特点，结合其弗罗德数的大小，可以认为坡面股流已属于急紊流范围。表2 坡面股流的弗罗德数Froude number of flow in rills on the slope流量（m/s）坡度（）61012150.610-41.0341.4471.4021.50011.0651.0851.1431.43421.0801.5831.4711.65631.0601.2371.4621.702

10、5-1.4031.1782.326表3 坡面股流的雷诺数Reynolds number of flow in rills on the slope172012401435185827982385246225153628343836603517502344305.4843665276762971004.2 细沟水力要素间的相互关系分析径流在顺坡流下的过程中，不断沿程侵蚀搬运土壤，在坡面上形成一条条细沟。在细沟的形成发展过程中，径流会通过不断地自我调整，力图使其能量消耗达到最小。径流的自我调整，使其各种水力要素此消彼涨，但其间存在着相互联系、相互制约的必然关系。深入认识细沟流各个水力要素间的关系及

11、其表现特征，对于揭示坡面水流的侵蚀动力学及各种水力要素的估算都十分重要。Lepold和Maddock曾指出5，河流断面的水力几何形态可以用建立于平均过水断面宽度B，平均径流深h，以及平均流速与流量之间的指数关系来描述B=aQb（7）h=cQd（8）U=eQm（9）其中 a，b，c，d，e，m都是经验系数，由质量连续定律可知，上述经验系数之间有关系：ace=1,b+d+m=1。也可以用同样的方法来研究细沟的水力学特征。为此，Lane和Foster6，Gilley7以及Govers等8研究了形成于缓坡条件下的细沟的水力学特征，结果表明上述关系式中的经验系数b变化于0.140.48之间，平均为0.3

12、；经验系数m变化于0.230.39，平均值为0.298，也近似地等于0.3，这样也就可以推算出d等于0.4。Govers的结果更表明m不受坡面坡度大小的影响，稳定在0.3左右。最近，Abrahams9研究了形成于草地上的细沟的水力特征后，发现b,d,m分别为0.33，0.34和0.33。我们的试验结果表明（图1，2），在试验的坡度和流量范围内，形成于黄土坡面上的细沟中的径流的水力要素与流量关系也满足指数方程，其中的指数b,d,m分别为0.26，0.48和0.26，即有关系式B=0.39Q0.26J-0.026（10）h=0.68Q0.48J-0.17（11）U=2.12Q0.26J0.25（12）图1 径流深随流量及坡度的变化图2 径流宽度随流量及坡度的变化Depth of flow versus discharge and slopeWidth of flow versus discharge and slope其