本发明涉及土壤侵蚀与水土保持、农田面源污染控制领域，具体地说是涉及一种浅沟侵蚀输沙能力确定的系统和方法。
背景技术：
浅沟是坡地侵蚀产沙的主要来源和输沙通道，发育有浅沟的坡地在降雨过程中迅速汇集径流，形成冲刷力很强的股流，侵蚀动能剧增，成为切沟侵蚀和沟头前进的动力源泉，使土地资源遭到极大的破坏。由于浅沟侵蚀对坡面总侵蚀输沙量的显著贡献，近十几年来引起了国内外对浅沟侵蚀模拟和预报研究的极大兴趣。目前可以模拟浅沟侵蚀过程的模型有wepp(watererosionpredictionproject)，但由于试验手段和计算方法的局限，产沙方程中输沙能力等尚不能经试验直接测定与计算，限制了模型在浅沟侵蚀预报中的应用。已开发的浅沟侵蚀模型egem(ephemeralgullyerosionmodel)可用于预报浅沟年均侵蚀量，但不能反应浅沟侵蚀产沙的时空变化过程，只能预报有限空间和时间点的土壤侵蚀信息。因缺乏合理的参数计算方法，当模拟值与实测值出现差异时，很难辨别误差来自模型本身还是参数选择有误。为了实现土壤侵蚀预报完全建立在物理过程基础上，模型参数应具有明确的物理意义，并能用力学表达式描述或由试验直接测量。
输沙能力(sedimenttransportcapacity)的概念是在物理过程分析的基础上得来的。在浅沟产生以后水流的侵蚀、搬运和沉积作用是一个不断转化的过程，它被定义为水流可携带的最大含沙量。如果实际水流含沙量小于这个值，那么水流继续侵蚀搬运土壤颗粒；如果含沙量大于输沙能力，那么多余的泥沙就会沉积。从能量的角度而言，浅沟中的水流具有剥离土壤和搬运泥沙的能力。在浅沟产生初期径流中的泥沙含量较少，水流的剥离能力占有主导地位，随着水流不断的剥蚀土壤泥沙量逐渐增多，用于剥离的能力也逐渐转换为输送泥沙的能量，直到水流的全部能量都用于输送泥沙，即不再剥离土壤。与此同时，水流中的泥沙含量处于饱和状态并达到稳定值，出现暂时的平衡状态。平衡状态下的稳定值对应于水流的输沙能力，表征由剥离作用进入水流中的泥沙量与水流中泥沙的沉积量相对平衡、净剥离量为零的一种状态。
技术实现要素：
本发明所要解决的技术问题是针对现有技术的不足，提供一种浅沟侵蚀输沙能力确定的定量表达方法。
输沙能力(sedimenttransportcapacity)的概念是在物理过程分析的基础上得来的。在浅沟产生以后水流的侵蚀、搬运和沉积作用是一个不断转化的过程，它被定义为水流可携带的最大含沙量。如果实际水流含沙量小于这个值，那么水流继续侵蚀搬运土壤颗粒；如果含沙量大于输沙能力，那么多余的泥沙就会沉积。从能量的角度而言，浅沟中的水流具有剥离土壤和搬运泥沙的能力。在浅沟产生初期径流中的泥沙含量较少，水流的剥离能力占有主导地位，随着水流不断的剥蚀土壤泥沙量逐渐增多，用于剥离的能力也逐渐转换为输送泥沙的能量，直到水流的全部能量都用于输送泥沙，即不再剥离土壤。与此同时，水流中的泥沙含量处于饱和状态并达到稳定值，出现暂时的平衡状态。平衡状态下的稳定值对应于水流的输沙能力，表征由剥离作用进入水流中的泥沙量与水流中泥沙的沉积量相对平衡、净剥离量为零的一种状态。在以上理论分析的基础上，我们提出了一种浅沟侵蚀输沙能力确定的系统和方法。
从能量耗散率基本理论出发，在给定的坡度和流量下，浅沟侵蚀的含沙量随径流对泥沙的剥蚀速率增加而增加，直到径流所能携带的输沙量达到最大值。应用初始条件，整合关于浅沟侵蚀下坡距离的输沙量方程，我们得出了一种解决方案，以计算在浅沟侵蚀中集中水流的输沙能力。在给定稳定水流的条件下，测定多个坡段不同浅沟长度水流的泥沙含量，借助数学方法得到泥沙含量与浅沟长度的对应关系。当浅沟长度达到临界值时，浅沟水流中的泥沙含量达到稳定或不再显着增加时，理论上水流中的泥沙含量为对应于浅沟水流输沙能力的泥沙含量。在坡度较缓坡长较短的情况下，即使浅沟长度达到试验水槽的长度时，水流中的泥沙含量仍未达到稳定，可以借助数学表达式，由试验测得的泥沙含量与浅沟长度的函数关系进行拟合，用拟合函数的极限值来计算浅沟侵蚀输沙能力。该测定系统和计算方法的关键是设计侵蚀水槽并测定浅沟侵蚀不同沟长对应的泥沙含量，得到浅沟侵蚀含沙量与沟长的对应关系。
本发明解决上述技术问题的技术方案如下：
本发明提供一种浅沟侵蚀输沙能力的定量表达方法，包括如下步骤：
s1，测定不同水动力条件下浅沟随不同沟长对应的泥沙含量，获取浅沟侵蚀含沙量与沟长的对应关系；
s2，根据s1中得到的浅沟侵蚀含沙量与沟长的对应关系，计算不同水动力条件下的浅沟侵蚀的单位时间水流的输沙量；
s3，根据s2中的数据得到不同水动力条件下单位时间水流的输沙量随沟长的变化复合非线性数学关系；
s4，根据s3中的输沙量随沟长的变化复合非线性数学关系，计算输沙能力。
进一步，在步骤s1之前，制备可精确控制水动力条件的侵蚀水槽，并使用所述侵蚀水槽进行步骤s1中的测定。
进一步，所述步骤s1中浅沟侵蚀含沙量与沟长的对应关系为：在预设的水动力条件条件下，浅沟沟长距离径流入水口xm-1的含沙量为cm-1，则在浅沟沟长距离径流入水口xm的含沙量为cm，表达式为
xm＝xm-1+δx(m∈[1，12])(1)
其中cm为浅沟从入水口至沟长xm处的含沙量，kg/m3；m是段数，因次；xm是浅沟从入水口至点m的沟长，m。
进一步，所述步骤s2中单位时间水流的输沙量g的计算公式为：
其中，其中g为水流中输沙量，kg/m·s；g＝qc，q是单位宽度径流量，m2/s。
进一步，所述步骤s3中不同水动力条件下单位时间水流的输沙量随沟长的变化复合非线性数学关系为：
g(x)＝a(1-e-βx)(5)
其中g为输沙量(kg/m·s)，β为衰减系数，x是沟长(m)，a是回归系数。
进一步，所述步骤s4中输沙能力的计算公式为：
α＝tcβ(6.1)
tc＝a(6.2)
其中，tc输沙能力(kgm-1s-1)，α为系数。
进一步，所述水动力条件包括径流量和坡度。
本发明的有益效果是：本发明的浅沟侵蚀输沙能力的定量表达方法，只需要建立浅沟输沙量和坡长的响应关系，就可以计算得到任意水动力条件下的浅沟侵蚀输沙能力；并且本发明的方法不完全基于简单的数学推导，数学式能明确表达泥沙剥蚀输移的反馈关系，为坡地剥蚀输沙的动态过程及作用机理研究提供了新手段。
附图说明
图1为本发明的浅沟侵蚀输沙能力的定量表达方法的流程示意图。
具体实施方式
以下结合附图对本发明的原理和特征进行描述，所举实例只用于解释本发明，并非用于限定本发明的范围。
一种浅沟侵蚀输沙能力的定量表达方法，该方法包括以下步骤：
浅沟侵蚀含沙量测定
从能量耗散率基本理论出发，在给定的坡度和径流量下，浅沟侵蚀的含沙量随径流对泥沙的剥蚀速率增加而增加，直到径流所能携带的输沙量达到最大值。应用初始条件，整合关于浅沟侵蚀下坡距离的输沙量方程，我们得出了一种解决方案，以计算在浅沟侵蚀中集中水流的输沙能力。在给定稳定水流的条件下，测定多个坡段不同浅沟长度水流的泥沙含量，借助数学方法得到泥沙含量与浅沟长度的对应关系。当浅沟长度达到临界值时，浅沟水流中的泥沙含量达到稳定或不再显着增加时，理论上水流中的泥沙含量为对应于浅沟水流输沙能力的泥沙含量。在坡度较缓坡长较短的情况下，即使浅沟长度达到试验水槽的长度时，水流中的泥沙含量仍未达到稳定，可以借助数学表达式，由试验测得的泥沙含量与浅沟长度的函数关系进行拟合，用拟合函数的极限值来计算浅沟侵蚀输沙能力。该测定系统和计算方法的关键是设计侵蚀水槽并测定浅沟侵蚀不同沟长对应的泥沙含量，得到浅沟侵蚀含沙量与沟长的对应关系。
为了精确地描述浅沟侵蚀产流产沙过程及因此造成的坡面形态变化过程，必须对研究对象抽象或简化，将影响浅沟侵蚀和发育的许多复杂因素尽可能地分离出来，在可控的条件下进行专门和系统地观测。从动力学的角度影响浅沟侵蚀产沙量的动力因子主要是径流强度、坡面径流运动速度、水力比降和沟长。由于坡面径流运动速度，尤其是点速度尚难以观测，因而很难建立起径流侵蚀产沙与流速的关系，所以在设计时，主要考虑了径流量和坡度对侵蚀产沙量的影响。因此，设计了室内浅沟侵蚀模拟试验，选取了径流量和坡度为影响浅沟侵蚀发生发展的主要水动力影响因子，测定不同水动力条件下浅沟随不同沟长对应的泥沙含量，获取浅沟侵蚀含沙量与沟长的对应关系。设在某一特定坡度和径流量条件下，浅沟沟长距离径流入水口xm-1的含沙量为cm-1，则在浅沟沟长距离径流入水口xm的含沙量为cm，表达式为
xm＝xm-1+δx(m∈[1，12])(1)
其中cm为浅沟从入水口至沟长xm处的含沙量，kg/m3；m是段数，因次；xm是浅沟从入水口至点m的沟长，m。
浅沟侵蚀输沙量计算
根据不同沟长所测得的含沙量结果，计算不同坡度和径流量下的浅沟侵蚀的输沙量。单位时间水流的输沙量g的公式为
其中g为水流中输沙量，kg/m·s；g＝qc，q是单位宽度流量，m2/s。
用公式(4)计算在不同径流量下在不同浅沟沟长的输沙量。在浅沟沟道径流引入口处，径流中的所有能量都可用于剥蚀土壤颗粒，因此，含沙量从入水口处随着沿浅沟沟长的增加而增加，随着含沙量的增加需要分配更多的能量用于输送剥蚀土壤颗粒，从而可用于剥蚀土壤颗粒的能量相应减少，含沙量沿沟长的增加速率趋于减少。
浅沟侵蚀输沙能力计算
不同坡度及流量下输沙量随沟长的变化复合非线性数学关系：
g(x)＝a(1-e-βx)(5)
其中,g为输沙量(kg/m·s)，β为衰减系数，x是沟长(m)，a是回归系数。
输沙能力计算
输沙量与坡长的函数所表达的物理意义为输沙量随着沟长的增加逐渐增加，但增加的幅度越来越小，逐渐趋近极限值a。β是输沙量随距离增加其增长幅度衰减的指数，a值实际上代表的是水流达到饱和时的输沙量，也就是一定坡度和径流量下水流的输沙能力：
α＝tcβ(6.1)
tc＝a(6.2)
其中，tc输沙能力(kgm-1s-1)，α为系数。
实施例1
在径流量为：32lmin-1，坡度为：5°的条件下，测量浅沟随不同沟长对应的泥沙含量，具体数值如下表所示：
表1浅沟随不同沟长对应的泥沙含量
由上表可推算：在坡度为5°和径流量为32lmin-1的条件下，浅沟沟长距离径流入水口xm-1的含沙量为cm-1，则在浅沟沟长距离径流入水口xm的含沙量为cm，表达式为
xm＝xm-1+δx(m∈[1，12])(1)
其中，cm为浅沟从入水口至沟长xm处的含沙量，kg/m3；m是段数，因次；xm是浅沟从入水口至点m的坡长，m。
根据上述的浅沟侵蚀含沙量与沟长的对应关系，浅沟侵蚀的单位时间水流的输沙量，计算公式如下：
其中，其中g为单位时间水流的输沙量，kg/m·s；g＝qc，q是单位宽度流量，m2/s。得到在坡度为5°和径流量为32lmin-1的条件下，不同坡段单位时间水流的输沙量见表1。
并可以推算出，单位时间水流的输沙量随沟长的变化复合非线性数学关系为：
g(x)＝a(1-e-βx)(5)
其中g为输沙量(kg/m·s)，β为衰减系数，x是沟长距离(m)，a是回归系数，表征一定坡度和流量下水流所能携带的最大输沙量。得到在坡度为5°和径流量为32lmin-1的条件下，含沙量最大值为440.12kg/m3,单位时间水流的最大输沙量为4.69kg/m·s。
输沙量与坡长的函数所表达的物理意义为输沙量随着沟长的增加逐渐增加，但增加的幅度越来越小，逐渐趋近极限值a。β是输沙量随距离增加其增长幅度衰减的指数，a值实际上代表的是水流达到饱和时的输沙量，也就是一定坡度和流量下水流的输沙能力，即：
α＝tcβ(6.1)
tc＝a(6.2)
其中，tc输沙能力(kgm-1s-1)，α为系数。
由此可以得出，在坡度为**和径流量为**的条件下浅沟的输沙能力tc为：4.69(kgm-1s-1)。
以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。