[发明专利]一种确定流域尺度次降雨泥沙来源的方法

说明书

技术领域

本发明涉及流域泥沙来源判断技术领域，尤其涉及一种确定流域尺度次降雨泥沙来源的方法。

背景技术

水土流失已成为我国严重的环境问题之一。水土流失不仅会降低土地生产力，泥沙进入河道也会破坏河流生态环境和堵塞河道；同时，污染物通常附着在泥沙表面，随泥沙进入河流从而造成地表水的污染。准确的判断泥沙来源及其贡献可以为做好水土保持和水资源管理提供重要的依据。

常用的指纹识别法可以弥补目前土壤侵蚀模型和土壤侵蚀测量方法在大流域尺度预测和切沟侵蚀预测中的不足。二十世纪七十年代，指纹识别法开始被应用到确定流域泥沙来源及其空间分布的研究。目前，此方法已被广泛应用在流域尺度的泥沙来源研究中。

使用指纹识别法，首先需要划分泥沙来源，其次要确定示踪因子，再次使用模型计算泥沙来源相对贡献率。通常指纹识别法采用混合模型计算泥沙来源的相对贡献率。但混合模型存在的问题是当泥沙来源数增加时，其预测结果准确率明显下降。Lees等(1997)研究了混合模型的预测结果随泥沙来源数的变化，他们发现当泥沙来源数从2个增加到6个时，模型预测成功率从75％下降到0。Small等(2002)提出，当泥沙来源数增加时，混合模型会预测某些泥沙来源无贡献或仅给出虚假的数值解。

发明内容

本发明的目的是提供一种确定流域尺度次降雨泥沙来源的方法，该方法可靠性较高，受泥沙来源个数的影响更小，可以研究次降雨空间产沙规律，计算不同降雨强度泥沙来源贡献率，对于水土保持时间和空间的规划具有指导意义。

本发明的目的是通过以下技术方案实现的：

一种确定流域尺度次降雨泥沙来源的方法，该方法包括：

将研究区划分为M个泥沙来源区，并采集每一泥沙来源区的表层土壤样品；

在次降雨径流过程中采集流域出口泥沙样品；

将所述每一泥沙来源区的土壤样品与泥沙样品分别进行预处理后，进行全元素测量；

基于判别分析法，从所有土壤样品的测量结果中筛选出能够区分泥沙来源区的化学元素，并建立对应的判别方程，再计算判别方程分类泥沙样品来源的正确比例；

根据判别方程中所包含的化学元素种类，将每一泥沙来源区的土壤样品以及泥沙样品对应化学元素及相应的浓度分别代入判别方程中，从而计算出每一泥沙来源区的中心位置，以及采集泥沙样品的位置；

根据每一泥沙来源区的中心位置、采集泥沙样品的位置以及判别方程分类泥沙样品来源的正确比例，计算每一泥沙来源区的中心位置与采集泥沙样品的位置之间的距离，从而确定每一泥沙来源区对泥沙样品的相对贡献。

进一步的，所述将所述每一泥沙来源区的土壤样品与泥沙样品分别进行预处理后，进行全元素测量包括：

使用去离子水分别将土壤样品与泥沙样品过湿筛，保留小于预定尺寸的泥沙颗粒；

将过湿筛的土壤样品与泥沙样品在预定温度下烘干后，使用电感耦合等离子体质谱仪进行全元素测量，获得每一土壤样品与泥沙样品所包含的化学元素种类及对应的浓度。

进一步的，所述基于判别分析法，从所有土壤样品的测量结果中筛选出能够区分泥沙来源区的化学元素，并建立对应的判别方程包括：

对比所有土壤样品的测量结果，筛选浓度含量不同的化学元素；

将浓度含量不同的化学元素在不同泥沙来源区的差别进行量化，筛选出量化结果大于预定值的化学元素，并建立对应的判别方程。

进一步的，所述根据判别方程中所包含的化学元素种类，将每一泥沙来源区的土壤样品以及泥沙样品对应化学元素及相应的浓度分别代入判别方程中，从而计算出每一泥沙来源区的中心位置，以及采集泥沙样品的位置包括：

当仅有一个判别方程时，将则将每一土壤样品以及泥沙样品分别代入判别方程中，获得的数值则为对应泥沙来源区的中心位置，以及采集泥沙样品的位置；

当有两个判别方程时，则将每一土壤样品以及泥沙样品分别代入这两个判别方程中，获得两个数值；将这两个数值分别作为二维坐标系中的x值和y值，从而计算出每一泥沙来源区的中心位置，以及采集泥沙样品的位置。

进一步的，所述根据每一泥沙来源区的中心位置、采集泥沙样品的位置以及判别方程分类泥沙样品来源的正确比例，计算每一泥沙来源区的中心位置与采集泥沙样品的位置之间的距离包括：

若判别方程的数量为n个，则第m个泥沙来源区的中心位置与采集泥沙样品的位置之间的距离D_m计算公式为：