一种适用于多地貌类型区的降雨输沙模型构建方法与流程

1.本发明涉及工程泥沙技术领域，特别涉及一种适用于多地貌类型区的降雨输沙模型构建方法。


背景技术：

2.流域产沙是降雨和下垫面结合的产物,无论下垫面还是降雨，一旦发生变化，都会产生不同的沙量。在分析流域产沙量变化原因时，需要分析人类活动影响较小的天然下垫面产沙量，在论证新建水库、淤地坝等工程时，需要分析流域现状下垫面的产沙量。在分析不同下垫面产沙量时，一般需要分析建模时段并构建降雨输沙模型。
3.黄土高原地貌复杂多样，包括黄土丘陵沟壑区、河源梁墹区、黄土高塬沟壑区、风沙区、黄土阶地区和黄土丘陵林区等多种地貌类型区，不同地貌类型区的产沙模数差别很大。对于包含多种地貌类型区的某一流域，如果在计算产沙量时不考虑地貌类型的差别，只采用一个降雨输沙模型计算，势必导致沙量计算结果偏大或者偏小。在确定建模时段时，由于水沙关系的复杂性，常用的双累积曲线法分析得出的突变点往往具有不确定性，会出现多个不同的突变点，在确定流域突变点时可能存在较多的人为因素。另外，在建立降雨输沙模型时，如果设定的参数初始值不当，可能导致计算收敛速度慢或者出现指数为负数等不符合物理概念的情况。


技术实现要素：

4.为解决上述技术问题，本发明提供一种适用于多地貌类型区的降雨输沙模型构建方法，用以计算不同下垫面输沙量，包括以下步骤：
5.采集流域雨量站降雨数据、水文站实测输沙量数据、不同地貌类型区面积数据；
6.根据所述流域雨量站降雨数据、水文站实测输沙量数据和不同地貌类型区面积数据分析出流域不同地貌类型区降雨特征指标及分区输沙量；
7.采用非参数统计检验方法分析建模代表时段，获取分析结果；
8.根据所述分析结果与流域不同地貌类型区降雨特征指标及分区输沙量构建不同地貌类型区的单指标降雨输沙模型，分析判别并筛选有效降雨指标；
9.基于筛选的有效降雨指标，构建分区多指标降雨输沙模型，对所述分区多指标降雨输沙模型中参数的取值范围进行量化约束，获得分区最优模型；
10.采用所述分区最优模型计算全流域输沙量，获取计算结果，根据所述计算结果对所述分区最优模型的合理性进行分析。
11.优选的，所述根据所述流域雨量站降雨数据、水文站实测输沙量数据和不同地貌类型区面积数据分析出流域不同地貌类型区降雨特征指标及分区输沙量，包括：
12.根据所述流域雨量站降雨数据，分析不同地貌类型区的降雨特征指标，根据所述水文站实测输沙量数据，分析不同地貌类型区的分区输沙量；
13.其中，所述不同地貌类型区指同一流域包含两种及两种以上地貌类型；
14.所述降雨特征指标包括时段降雨量p、时段平均雨强i、时段降雨笼罩面积f、不同等级降雨量p
n
、不同等级雨强i
n
、不同等级降雨笼罩面积f
n
及最大n日雨量p
maxn
；
15.所述分区输沙量包括不同地貌类型区的年输沙量w
s分区
，根据各地貌分区内水文站实测沙量w
s水文站
，采用面积比法进行计算：
[0016][0017]
其中，f
分区
为不同地貌类型区的面积，f
水文站
为各地貌分区内水文站的集水面积。
[0018]
优选的，所述采用非参数统计检验方法分析建模代表时段，获取分析结果，包括：
[0019]
基于mann
‑
whitney
‑
pettitt法分析不同地貌类型区的产沙系数(产沙系数＝输沙量/(面积*降水量))统计量ut，检测出产沙系数序列的一级突变年份，然后以该突变年份为界将原沙量序列一分为二，继续检测新的突变年份，可检测出二级突变年份，从而获得所述建模代表时段；
[0020]
将所述建模代表时段作为所述分析结果。
[0021]
优选的，所述根据所述分析结果与流域不同地貌类型区降雨特征指标及分区输沙量构建不同地貌类型区的单指标降雨输沙模型，分析判别并筛选有效降雨指标，包括：
[0022]
根据所述分区输沙量与降雨特征指标，采用非线性回归方法，分别建立分区沙量时段降雨量p、时段平均雨强i、时段降雨笼罩面积f、不同等级降雨量p
n
、不同等级雨强i
n
、不同等级降雨笼罩面积f
n
及最大n日雨量p
maxn
等单个降雨指标的关系，从中筛选出与输沙量相关性大于等于预设阈值的多组单个降雨指标，作为有效降雨指标，单指标降雨输沙模型构建方法如下：
[0023]
w
s分区
＝αp
β
[0024]
其中，α为系数，β为指数。
[0025]
优选的，所述基于筛选的有效降雨指标，构建分区多指标降雨输沙模型，对所述分区多指标降雨输沙模型中参数的取值范围进行量化约束，选择分区最优模型，包括：
[0026]
基于筛选的有效降雨指标，进行两两降雨指标组合，组合方式为时段降雨量与平均雨强、时段降雨量与笼罩面积、时段降雨量与最大n日雨量、不同等级降雨量与雨强、不同等级降雨量与笼罩面积进行组合，获取多组组合后的降雨指标；
[0027]
根据每组组合的降雨指标，构建每组组合的降雨指标对应的分区多指标降雨输沙模型：
[0028]
w
s分区
＝αp
β
i
δ
[0029]
其中，α为系数，β、δ为指数；
[0030]
根据黄河流域主要产沙区数十组降雨输沙模型总结分析结果对所述分区多指标降雨输沙模型中的参数的取值范围进行量化约束；
[0031]
在多组多指标降雨输沙模型中选择相关度最高的一组作为采用的分区最优模型。
[0032]
优选的，所述采用所述分区最优模型计算全流域输沙量，获取计算结果，根据所述计算结果对所述分区最优模型的合理性进行分析，包括：
[0033]
利用所述分区最优模型计算各分区输沙量，累加后计算全流域输沙量；
[0034]
将计算的全流域输沙量与流域把口水位站的实测输沙量进行对比分析，若二者的
相似度大于等于第三预设阈值，确认所述分区最优模型的合理性为优，否则重新选降雨输沙模型进行计算。
[0035]
优选的，所述方法还包括：对所述分区最优模型进行稳定性检测，检测步骤包括：
[0036]
设定预设坐标系，在所述流域中固定设置多个不同含沙量的沙坑作为特征点；
[0037]
将各个特征点输入到所述分区最优模型中，获得所述分区最优模型输出的特征数据；
[0038]
利用预设图像采集装置在所述特征数据中采集每个沙坑的图像数据，对每个沙坑的图像数据进行处理获得不同特征点的实时坐标值；
[0039]
根据不同特征点的实时坐标值获得所述流域的第一数字化地貌曲面；
[0040]
对所述流域进行模拟间接性降雨，检测所述模拟间接性降雨过程中的第一降雨序列；
[0041]
获取所述第一降雨序列中的第一分位数，根据所述分位数将所述第一降雨序列分为第一极小值、第一常规值和第一极大值；
[0042]
确定流域内雨水的实测值，根据所述实测值推算出流域实测降雨过程中的第二降雨序列；
[0043]
获取所述第二降雨序列中的第二分位数，根据所述分位数将所述第二降雨序列分为第二极小值、第二常规值和第二极大值；
[0044]
根据所述第二极小值、第二常规值和第二极大值与第一极小值、第一常规值和第一极大值之间的偏差值确定第一降雨序列与第二降雨序列之间的修正因子；
[0045]
基于所述修正因子，对所述第一降雨序列进行修正，获得修正后的第一降雨序列；
[0046]
将所述修正后的第一降雨序列中的多个序列值转化为在所述预设坐标系中的目标坐标值；
[0047]
根据多个目标坐标值获得所述流域的第二数字化地貌曲面；
[0048]
比较所述第一数字化地貌曲面与第二数字化地貌曲面，获得所述流域内不同地貌的形变幅度和速度；
[0049]
获取流域进行模拟间接性降雨的降雨参数，根据所述流域内不同地貌的形变幅度和速度以及所述降雨参数构建流域在不同降雨条件下的稳定性评估模型；
[0050]
获取第一目标数量个预设降雨量，将所述预设降雨量同时输入到所述分区最优模型和稳定性评估模型中，获得第一目标数量个第一输出结果和第二输出结果；
[0051]
确认第一目标数量个第一输出结果与第二输出结果的相似度；
[0052]
统计第一输出结果与第二输出结果的相似度大于等于第四预设阈值的第二目标数量，确认所述第二目标数量是否大于等于预设数量，若是，确认所述分区最优模型的稳定性为优，否则，确认所述分区最优模型的稳定性为差。
[0053]
优选的，所述采用所述分区最优模型计算全流域输沙量，获取计算结果，根据所述计算结果对所述分区最优模型的合理性进行分析，包括：
[0054]
确定所述分区最优模型对应的目标有效降雨指标；
[0055]
确定所述目标有效降雨指标对应的目标降雨强度；
[0056]
计算出目标降雨强度相对于预设基准降雨强度的地面侵蚀强度：
[0057][0058]
其中，k表示为目标降雨强度相对于预设基准降雨强度的地面侵蚀强度，d
i
表示为在目标降雨强度下流域的地面侵蚀量，d
i
表示为在预设基准降雨强度下流域的地面侵蚀量，d1表示为在目标降雨强度与预设基准降雨强度差值的绝对值下流域的地面侵蚀量，δ表示为分区最优模型的精度；
[0059]
根据目标降雨强度相对于预设基准降雨强度的地面侵蚀强度计算出全流域输沙量：
[0060][0061]
其中，q表示为全流域输沙量，n表示为降雨的天数，s
i
表示为第i天的有效降雨量，f
i
表示为第i天的降雨衰减因子，m表示为流域内地貌的数量，g
j
表示为第j个地貌的土地稀松度，x
j
表示为第j个地貌的输沙模数，θ1表示为流域内的地下产沙比重，θ2表示为流域内的地表产沙比重，m表示为流域内的下垫面的总质量；
[0062]
将所述全流域输沙量确认为计算结果；
[0063]
根据所述全流域输沙量与实际检测的输沙量进行比较，若二者的偏差在预设范围内，确认所述分区最优模型的合理性为合理，若二者的偏差在所述预设范围外，确认所述分区最优模型的合理性为不合理。
[0064]
与现有技术相比，本发明具有以下优点和有益效果：
[0065]
本发明为一种适用于多地貌类型区的降雨输沙模型构建方法，充分考虑了不同地貌类型区的产沙特点。在分析建模时段时，采用了非参数统计检验方法，基于原始数据所获得的秩进行统计分析，也不受少数异常值和缺失数据的干扰和影响，分析的建模时段更符合实际。在构建模型时，首先通过构建分区单指标降雨输沙模型，优选分区有效降雨指标，再通过指标合理组合，构建出分区多指标降雨输沙模型，并对模型中的指数取值范围进行量化约束，提高了回归分析的准确性及高效性，选择最优降雨输沙模型，为不同下垫面输沙量计算提供了一种技术解决方案，步骤简明、易于操作、成果可靠，是一种易于为基层科技工作者掌握和使用的简便方法。
[0066]
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
[0067]
附图用来提供对本发明的进一步理解，并且构成说明书的一部分，与本发明的实施例一起用于解释本发明，并不构成对本发明的限制。在附图中：
[0068]
图1为本发明提供的适用于多地貌类型区的降雨输沙模型构建方法的流程图；
[0069]
图2为适应于本发明的典型多地貌类型流域示意图；
[0070]
图3为本发明提供的建模时段划分示例图；
[0071]
图4为本发明提供的某一地貌类型区单指标降雨输沙模型示例图；
[0072]
图5为本发明提供的某一地貌类型区多指标降雨输沙模型示例图；
[0073]
图6为本发明提供的适用于多地貌类型区的降雨输沙模型构建系统的结构示意
图。
具体实施方式
[0074]
以下结合附图对本发明的优选实施例进行说明，应当理解，此处所描述的优选实施例仅用于说明和解释本发明，并不用于限定本发明。
[0075]
另外，在本发明中如涉及“第一”、“第二”等的描述仅用于描述目的，并非特别指称次序或顺位的意思，亦非用以限定本发明，其仅仅是为了区别以相同技术用语描述的组件或操作而已，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外，各个实施例之间的技术方案以及技术特征可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本发明要求的保护范围之内。
[0076]
流域产沙是降雨和下垫面结合的产物,无论下垫面还是降雨，一旦发生变化，都会产生不同的沙量。在分析流域产沙量变化原因时，需要分析人类活动影响较小的天然下垫面产沙量，在论证新建水库、淤地坝等工程时，需要分析流域现状下垫面的产沙量。在分析不同下垫面产沙量时，一般需要分析建模时段并构建降雨输沙模型。
[0077]
黄土高原地貌复杂多样，包括黄土丘陵沟壑区、河源梁墹区、黄土高塬沟壑区、风沙区、黄土阶地区和黄土丘陵林区等多种地貌类型区，不同地貌类型区的产沙模数差别很大。对于包含多种地貌类型区的某一流域，如果在计算产沙量时不考虑地貌类型的差别，只采用一个降雨输沙模型计算，势必导致沙量计算结果偏大或者偏小。在确定建模时段时，由于水沙关系的复杂性，常用的双累积曲线法分析得出的突变点往往具有不确定性，会出现多个不同的突变点，在确定流域突变点时可能存在较多的人为因素。另外，在建立降雨输沙模型时，如果设定的参数初始值不当，可能导致计算收敛速度慢或者出现指数为负数等不符合物理概念的情况为了解决上述问题，本实施例公开了一种适用于多地貌类型区的降雨输沙模型构建方法。
[0078]
一种适用于多地貌类型区的降雨输沙模型构建方法，如图1所示，包括以下步骤：
[0079]
步骤1：采集流域雨量站降雨数据、水文站实测输沙量数据、不同地貌类型区面积数据；
[0080]
步骤2：根据所述流域雨量站降雨数据、水文站实测输沙量数据和不同地貌类型区面积数据分析出流域不同地貌类型区降雨特征指标及分区输沙量；
[0081]
步骤3：采用非参数统计检验方法分析建模代表时段，获取分析结果；
[0082]
步骤4：根据所述分析结果与流域不同地貌类型区降雨特征指标及分区输沙量构建不同地貌类型区的单指标降雨输沙模型，分析判别并筛选有效降雨指标；
[0083]
步骤5：基于筛选的有效降雨指标，构建分区多指标降雨输沙模型，对所述分区多指标降雨输沙模型中参数的取值范围进行量化约束，获得分区最优模型；
[0084]
步骤6：采用所述分区最优模型计算全流域输沙量，获取计算结果，根据所述计算结果对所述分区最优模型的合理性进行分析；
[0085]
在本实施例中，上述步骤1可以从大数据库中调取到上述数据并且进行下载以利用，适用于本方法的典型多地貌类型流域示意图如图2所示。
[0086]
上述步骤2包括：根据所述流域雨量站降雨数据，分析不同地貌类型区的降雨特征指标，根据所述水文站实测输沙量数据，分析不同地貌类型区的分区输沙量；
[0087]
其中，所述不同地貌类型区指同一流域包含两种及两种以上地貌类型；
[0088]
所述降雨特征指标包括时段降雨量p、时段平均雨强i、时段降雨笼罩面积f、不同等级降雨量p
n
、不同等级雨强i
n
、不同等级降雨笼罩面积f
n
及最大n日雨量p
maxn
；
[0089]
所述分区输沙量包括不同地貌类型区的年输沙量w
s分区
，根据各地貌分区内水文站实测沙量w
s水文站
，采用面积比法进行计算：
[0090][0091]
其中，f
分区
为不同地貌类型区的面积，f
水文站
为各地貌分区内水文站的集水面积。
[0092]
上述步骤3包括：基于mann
‑
whitney
‑
pettitt法分析不同地貌类型区的产沙系数(输沙量/面积/降水量)统计量ut，检测出产沙系数序列的一级突变年份，然后以该突变年份为界将原沙量序列一分为二，继续检测新的突变年份，可检测出二级突变年份，从而获得所述建模代表时段；
[0093]
将所述建模代表时段作为所述分析结果；
[0094]
建模代表时段的划分示例图如图3所示。
[0095]
上述步骤4包括：根据所述分区输沙量与降雨特征指标，采用非线性回归方法，分别建立分区沙量时段降雨量p、时段平均雨强i、时段降雨笼罩面积f、不同等级降雨量p
n
、不同等级雨强i
n
、不同等级降雨笼罩面积f
n
及最大n日雨量p
maxn
等单个降雨指标的关系，从中筛选出与输沙量相关性大于等于预设阈值的多组单个降雨指标，作为有效降雨指标，单指标降雨输沙模型构建方法如下：
[0096]
w
s分区
＝αp
β
[0097]
其中，α为系数，β为指数；
[0098]
图4为根据本发明方法构建的某一地貌类型区单指标降雨输沙模型示例图。
[0099]
上述步骤5包括：基于筛选的有效降雨指标，进行两两降雨指标组合，组合方式为时段降雨量与平均雨强(pi)、时段降雨量与笼罩面积(pf)、时段降雨量与最大n日雨量(pp
maxn
)、不同等级降雨量与雨强(p
n
i
n
)、不同等级降雨量与笼罩面积进行组合(p
n
i
n
)，获取多组组合后的降雨指标；
[0100]
根据每组组合的降雨指标，构建每组组合的降雨指标对应的分区多指标降雨输沙模型：
[0101]
w
s分区
＝αp
β
i
δ
[0102]
其中，α为系数，β、δ为指数；
[0103]
根据黄河流域主要产沙区数十组降雨输沙模型总结分析结果对所述分区多指标降雨输沙模型中的参数的取值范围进行量化约束；
[0104]
在多组多指标降雨输沙模型中选择相关度最高的一组作为采用的分区最优模型；
[0105]
图5为根据本发明方法构建的某一地貌类型区多指标降雨输沙模型示例图。
[0106]
上述步骤6包括：利用所述分区最优模型计算各分区输沙量，累加后计算全流域输沙量；
[0107]
将计算的全流域输沙量与流域把口水位站的实测输沙量进行对比分析，若二者的相似度大于等于第三预设阈值，确认所述分区最优模型的合理性为优，否则重新选降雨输沙模型进行计算。
[0108]
上述技术方案的工作原理为：采集流域雨量站降雨数据、水文站实测输沙量数据、不同地貌类型区面积数据，根据所述流域雨量站降雨数据、水文站实测输沙量数据和不同地貌类型区面积数据分析出流域不同地貌类型区降雨特征指标及分区输沙量，采用非参数统计检验方法分析建模代表时段，获取分析结果，根据所述分析结果与流域不同地貌类型区降雨特征指标及分区输沙量构建不同地貌类型区的单指标降雨输沙模型，分析判别并筛选有效降雨指标，基于筛选的有效降雨指标，构建分区多指标降雨输沙模型，对所述分区多指标降雨输沙模型中参数的取值范围进行量化约束，获得分区最优模型，采用所述分区最优模型计算全流域输沙量，获取计算结果，根据所述计算结果对所述分区最优模型的合理性进行分析。
[0109]
上述技术方案的有益效果为：充分考虑了不同地貌类型区的产沙特点。在分析建模时段时，采用了非参数统计检验方法，基于原始数据所获得的秩进行统计分析，也不受少数异常值和缺失数据的干扰和影响，分析的建模时段更符合实际。在构建模型时，首先通过构建分区单指标降雨输沙模型，优选分区有效降雨指标，再通过指标合理组合，构建出分区多指标降雨输沙模型，并对模型中的指数取值范围进行量化约束，提高了回归分析的准确性及高效性，选择最优降雨输沙模型，为不同下垫面输沙量计算提供了一种技术解决方案，步骤简明、易于操作、成果可靠。
[0110]
在一个实施例中，所述方法还包括：对所述分区最优模型进行稳定性检测，检测步骤包括：
[0111]
设定预设坐标系，在所述流域中固定设置多个不同含沙量的沙坑作为特征点；
[0112]
将各个特征点输入到所述分区最优模型中，获得所述分区最优模型输出的特征数据；
[0113]
利用预设图像采集装置在所述特征数据中采集每个沙坑的图像数据，对每个沙坑的图像数据进行处理获得不同特征点的实时坐标值；
[0114]
根据不同特征点的实时坐标值获得所述流域的第一数字化地貌曲面；
[0115]
对所述流域进行模拟间接性降雨，检测所述模拟间接性降雨过程中的第一降雨序列；
[0116]
获取所述第一降雨序列中的第一分位数，根据所述分位数将所述第一降雨序列分为第一极小值、第一常规值和第一极大值；
[0117]
确定流域内雨水的实测值，根据所述实测值推算出流域实测降雨过程中的第二降雨序列；
[0118]
获取所述第二降雨序列中的第二分位数，根据所述分位数将所述第二降雨序列分为第二极小值、第二常规值和第二极大值；
[0119]
根据所述第二极小值、第二常规值和第二极大值与第一极小值、第一常规值和第一极大值之间的偏差值确定第一降雨序列与第二降雨序列之间的修正因子；
[0120]
基于所述修正因子，对所述第一降雨序列进行修正，获得修正后的第一降雨序列；
[0121]
将所述修正后的第一降雨序列中的多个序列值转化为在所述预设坐标系中的目
标坐标值；
[0122]
根据多个目标坐标值获得所述流域的第二数字化地貌曲面；
[0123]
比较所述第一数字化地貌曲面与第二数字化地貌曲面，获得所述流域内不同地貌的形变幅度和速度；
[0124]
获取流域进行模拟间接性降雨的降雨参数，根据所述流域内不同地貌的形变幅度和速度以及所述降雨参数构建流域在不同降雨条件下的稳定性评估模型；
[0125]
获取第一目标数量个预设降雨量，将所述预设降雨量同时输入到所述分区最优模型和稳定性评估模型中，获得第一目标数量个第一输出结果和第二输出结果；
[0126]
确认第一目标数量个第一输出结果与第二输出结果的相似度；
[0127]
统计第一输出结果与第二输出结果的相似度大于等于第四预设阈值的第二目标数量，确认所述第二目标数量是否大于等于预设数量，若是，确认所述分区最优模型的稳定性为优，否则，确认所述分区最优模型的稳定性为差。
[0128]
上述技术方案的有益效果为：通过利用数据坐标检测流域内的第一数字化地貌曲面可实现对流域的相对高频检测，以达到精确监测的目的，提高数据的准确性，进一步地，通过对流域进行模拟降雨获得流域在降雨过程中的第二数字化地貌曲面可以快速地与第一数字化地貌曲面进行比较以确定流域在降雨情况下的地形变化，进一步地，通过利用降雨参数和流域在降雨情况下的地形变化参数构建稳定性评估模型可以保证稳定性评估模型与分区最优模型的数据关联性，进而通过对比二者的相似度即可确定分区最优模型的稳定性，避免了人为主观因素对于分区最优模型的稳定性判断结果的影响，保证最终判定结果的客观性。
[0129]
在一个实施例中，所述采用所述分区最优模型计算全流域输沙量，获取计算结果，根据所述计算结果对所述分区最优模型的合理性进行分析，包括：
[0130]
确定所述分区最优模型对应的目标有效降雨指标；
[0131]
确定所述目标有效降雨指标对应的目标降雨强度；
[0132]
计算出目标降雨强度相对于预设基准降雨强度的地面侵蚀强度：
[0133][0134]
其中，k表示为目标降雨强度相对于预设基准降雨强度的地面侵蚀强度，d
i
表示为在目标降雨强度下流域的地面侵蚀量，d
i
表示为在预设基准降雨强度下流域的地面侵蚀量，d1表示为在目标降雨强度与预设基准降雨强度差值的绝对值下流域的地面侵蚀量，δ表示为分区最优模型的精度；
[0135]
根据目标降雨强度相对于预设基准降雨强度的地面侵蚀强度计算出全流域输沙量：
[0136][0137]
其中，q表示为全流域输沙量，n表示为降雨的天数，s
i
表示为第i天的有效降雨量，f
i
表示为第i天的降雨衰减因子，m表示为流域内地貌的数量，g
j
表示为第j个地貌的土地稀松度，x
j
表示为第j个地貌的输沙模数，θ1表示为流域内的地下产沙比重，θ2表示为流域内的
地表产沙比重，m表示为流域内的下垫面的总质量；
[0138]
将所述全流域输沙量确认为计算结果；
[0139]
根据所述全流域输沙量与实际检测的输沙量进行比较，若二者的偏差在预设范围内，确认所述分区最优模型的合理性为合理，若二者的偏差在所述预设范围外，确认所述分区最优模型的合理性为不合理。
[0140]
上述技术方案的有益效果为：通过计算目标降雨强度相对于预设基准降雨强度的地面侵蚀强度可以预先评估在目标降雨强度下对应地面的侵蚀强度以及产沙的效率，进一步地，通过利用目标降雨强度相对于预设基准降雨强度的地面侵蚀强度计算出全流域输沙量可以根据流域内的不同地貌的的参数以及降雨量来准确地计算出全流域输沙量，保证了数据的精确性。
[0141]
优选的，本实施例还公开了一种适用于多地貌类型区的降雨输沙模型构建方法的系统，如图6所示，该系统包括：
[0142]
采集模块601，用于采集流域雨量站降雨数据、水文站实测输沙量数据、不同地貌类型区面积数据；
[0143]
第一分析模块602，用于根据所述流域雨量站降雨数据、水文站实测输沙量数据和不同地貌类型区面积数据分析出流域不同地貌类型区降雨特征指标及分区输沙量；
[0144]
第二分析模块603，采用非参数统计检验方法分析建模代表时段，获取分析结果；
[0145]
筛选模块604，用于根据所述分析结果与流域不同地貌类型区降雨特征指标及分区输沙量构建不同地貌类型区的单指标降雨输沙模型，分析判别并筛选有效降雨指标；
[0146]
构建模块605，用于基于筛选的有效降雨指标，构建分区多指标降雨输沙模型，对所述分区多指标降雨输沙模型中参数的取值范围进行量化约束，获得分区最优模型；
[0147]
计算模块606，用于采用所述分区最优模型计算全流域输沙量，获取计算结果，根据所述计算结果对所述分区最优模型的合理性进行分析。
[0148]
上述技术方案的工作原理及有益效果在方法权利要求中已经说明，此处不再赘述。