云南山地中小河流众多,但许多流域内水文站较少,且很少进行泥沙观测,缺乏充足的泥沙资料,给水利工程设计和建设带来一定的困难。目前中小河流水库设计工作中泥沙大多以土壤侵蚀量近似代表河道输沙量进行泥沙计算,此方法明显较实际输沙量偏大不少,增加了无谓的工程成本,因此急需探讨一种精度较高的无资料地区河流输沙量估算方法。本研究基于流域的产沙机制,考察云南20个水文站的多年实测悬移质泥沙数据和流域内的土地利用、土壤分布等流域环境因素的关系,建立了计算河流输沙量的多元回归方程,为无水文站监测的小流域地区水库设计提供可信度较高的泥沙数据,亦为类似区域的输沙量估算提供参考。 There are many medium and small rivers in Yunnan mountainous areas. However, there are few hydrological stations in many basins, sediment observation is rarely carried out in these hydro-logical stations, which leads to the lack of sufficient sediment data, and brings some difficulty to the design and construction of water conservancy projects. At present, most of sediment data are approximately represented by soil erosion when design medium and small reservoir, which are larger than actual sediment discharge and increase the unnecessary engineering cost. Therefore, it’s urgent to find a method with high accuracy to estimate river sediment discharge in medium and small rivers. This study is based on sediment yield mechanism of basin, which analyzed the relationship between suspended sediment data from 20 hydrological stations in Yunnan and land-use, soil distribution and some other environmental factors, and established suspended sediment calculation model. This study provides reliable sediment data for reservoir design in small basin without hydrological station monitoring, and also provides reference for sediment discharge estimation in similar areas.
云南山地中小河流众多,但许多流域内水文站较少,且很少进行泥沙观测,缺乏充足的泥沙资料,给水利工程设计和建设带来一定的困难。目前中小河流水库设计工作中泥沙大多以土壤侵蚀量近似代表河道输沙量进行泥沙计算,此方法明显较实际输沙量偏大不少,增加了无谓的工程成本,因此急需探讨一种精度较高的无资料地区河流输沙量估算方法。本研究基于流域的产沙机制,考察云南20个水文站的多年实测悬移质泥沙数据和流域内的土地利用、土壤分布等流域环境因素的关系,建立了计算河流输沙量的多元回归方程,为无水文站监测的小流域地区水库设计提供可信度较高的泥沙数据,亦为类似区域的输沙量估算提供参考。
河流输沙量,无泥沙资料地区,中小河流,云南
Kaidao Fu1, Mei Ning1, Yating Dai1, Zhengtao Shi2, Tianli Zhang3, Shixiang Gu3*
1Institute of International Rivers and Eco-Security, Yunnan University, Kunming Yunnan
2College of Tourism and Geography Science, Yunnan Normal University, Kunming Yunnan
3Yunnan Institute of Water & Hydropower Engineering Investigation, Design and Research, Kunming Yunnan
Received: Jan. 26th, 2021; accepted: Mar. 10th, 2021; published: Mar. 18th, 2021
There are many medium and small rivers in Yunnan mountainous areas. However, there are few hydrological stations in many basins, sediment observation is rarely carried out in these hydrological stations, which leads to the lack of sufficient sediment data, and brings some difficulty to the design and construction of water conservancy projects. At present, most of sediment data are approximately represented by soil erosion when design medium and small reservoir, which are larger than actual sediment discharge and increase the unnecessary engineering cost. Therefore, it’s urgent to find a method with high accuracy to estimate river sediment discharge in medium and small rivers. This study is based on sediment yield mechanism of basin, which analyzed the relationship between suspended sediment data from 20 hydrological stations in Yunnan and land-use, soil distribution and some other environmental factors, and established suspended sediment calculation model. This study provides reliable sediment data for reservoir design in small basin without hydrological station monitoring, and also provides reference for sediment discharge estimation in similar areas.
Keywords:River sediment Discharge, Rivers without Sediment Data, Medium and Small Rivers, Yunnan
Copyright © 2021 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
伴随经济迅速发展,资源需求逐年增加。在化石能源日渐枯竭的今天,能源的转换和利用已成为许多国内科学家的研究课题,其中水利发电作为能源转换的重要组成,也广受关注 [
云南省有怒江、澜沧江等6大水系,水力资源丰富,水利水电工程建设较多。在干流或较大的支流上规划建设大型水电工程时,会在坝前或坝后建立水文站,进行泥沙观测。但在中小河流流域内建成的水文站较少,导致工程建设时缺乏充足的泥沙观测资料进行参考,给工程建设带来一定的阻碍 [
国内外关于流域输沙量计算的研究大多是计算流域土壤侵蚀量,再通过泥沙输移比 [
针对现存无资料河流来沙量估算方法的误差和不足,本文从流域产沙机制出发,分析流域内坡耕地面积、土壤性质、植被面积等因素对流域内悬移质输沙量的影响,构建计算悬移质输沙量的模型,用以推求无资料中小河流的输沙量 [
由于建造水利水电工对河流的径流量及下游输沙量产生一定影响 [
图1. 云南省水文站分布图
80年代土地利用数据缺乏,选择年份最接近水文数据的2000年土地利用类型图进行耕地和植被面积提取。而流域坡度、土壤类型在短期内变化不大,因此选择精度和分辨率较高的2009年DEM数据与1995年土壤类型数据(表1)。
数据名称 | 日期(年) | 来源 | 格式 | 分辨率 | 用途 |
---|---|---|---|---|---|
数字高程图 | 2009 | 地理空间数据云 | TIFF | 30 m | 提取集水面积 |
土地利用类型图 | 2000 | 中国西部环境与生态科学数据中心 | TIFF | 1:10万 | 提取流域坡耕地面积、植被面积 |
土壤类型图 | 1995 | 国家地球系统科学数据共享服务平台 | shp | 500 m | 计算土壤系数 |
泥沙资料 | 1981~1985,1987 | 藏南滇西河流水文资料 | Excel表格 | — | 建立回归模型 |
表1. 研究区基础数据
选取拟建的永平岔河水库对模型进行验证,水库建坝断面有2014~2016年实测输沙量数据。同时,采用《云南省土壤侵蚀遥感调查报告》(云南省水利水电厅2004年8月新编)的土壤侵蚀图表,结合卫星图片(2013年4月10日现场)显示的地表植被状况进行综合估算,得到岔河水库多年平均总侵蚀量为18.39万吨(表2)。
数据来源 | 输沙量(万吨) |
---|---|
水文站测量 | 5.97 |
侵蚀模数估算 | 18.39 |
表2. 岔河流域输沙量数据
单位径流深输沙模数是指流域内单位径流深在每单位面积上对土壤产生侵蚀并输移至河流中的年含沙量 [
m s = M s / y 。
式中:Ms——流域多年平均输沙模数(t/km2)。
Y——流域多年平均年径流深(mm)。
ms——单位径流深输沙模数(t/(km2∙mm))。
单位径流深输沙模数反映的是流域内地表的可侵蚀程度,它的大小是由流域内可侵蚀松散土的提供量以及侵蚀表土在流域内的输运条件决定的。当前许多研究采用流域内植被覆盖率来衡量流域的侵蚀程度。植被能有效减少土壤侵蚀,降低流域产沙 [
云南是高原山地广布的省份,气候条件、地质、地貌、土壤和植被类型复杂多变。山地河流地区虽构造活跃,地层破碎,但多为基岩,土壤层较薄,难以提供河流悬移质泥沙物质。而高原面上人口聚居,流域面上耕地面积大 [
根据分析,单位径流深输沙模数的函数表达为如下函数关系式 [
m s = f ( F , P , P K , K , V ) 。
其中:ms——单位径流深输沙模数(t/(km2∙mm))。
F——流域面积(km2)。
P——流域内坡耕地面积百分数(%)。
PK——坡度大于25˚耕地占坡耕地面积百分数(%)。
K——流域内土壤可蚀性(t∙hm2∙h/(MJ∙mm∙hm2))。
V——流域内植被面积百分数(%)。
在ArcGIS10.2环境下,利用云南省数字高程图、土壤类型图、土地利用类型图等数据,提取各个流域的影响因子。利用SPSS软件对20个站的各影响因子数据进行回归分析,建立单位径流深输沙模数方程。
选取拟建的永平岔河水库对构建的模型精度验证。通过ArcGIS提取永平流域的各影响因子数据并代入单位径流深输沙模数方程,计算得出永平流域悬移质输沙量。将计算结果与水文站实测输沙量、侵蚀模数估算结果进行对比,验证模型精度。
由于影响单位径流深输沙模数的影响因子包括集水面积、流域内坡耕地面积百分比、流域内大于25˚坡耕地面积百分比、土壤可蚀性和流域内植被面积百分比。根据单位径流深输沙模数的表达式,在ArcGIS10.2环境下,利用云南省数字高程图、土壤类型图、土地利用类型图等,提取各个流域的影响因子,计算出集水面积、坡耕地面积、大于25˚耕地面积和土壤系数。
将提取的集水面积与水文站给出的数据进行对比,发现各流域的F因子与实际集水面积误差仅为1%左右,分析可能是由于投影和DEM数据源选择不同造成的。将流域集水面积图与中国西部环境与生态科学数据中心获得的流域土地利用类型图、国家地球系统科学数据共享服务平台下载的土壤类型图进行叠加分析,得到其他影响因子数据,计算结果较为精确。
利用SPSS软件将云南境内20个水文站控制流域80年代多年实测河流泥沙资料及相应流域内土地利用资料,将5个影响因子数据与单位径流深输沙模数进行多元回归,得出如下经验关系:
m = K P 0.435 P 0.225 F − 0.092 V − 0.05 。
式中:ms——单位径流深输沙模数(t/(km2∙mm))。
K——流域内土壤可蚀性(t∙hm2∙h/(MJ∙mm∙hm2))。
P——流域内坡耕地面积百分数(%)。
PK——坡度大于25˚耕地占坡耕地面积百分数(%)。
F——流域面积(km2)。
V——流域内植被面积百分数(%)
单位径流深输沙模数方程计算得出的ms与实测值进行比较,发现:
1) 计算值与实测值相对误差小于30%的流域占总数的60%。
2) 当实测值小于2时,计算值与实测值相对误差小于30%的点占据86%,即实测单位径流深输沙模数小于2 t/(km2∙mm)时,模型计算结果更接近实测值,能更好拟合实测输沙量。
3) 当实测单位径流深输沙模数大于2 t/(km2∙mm)时,方程的计算结果误差较大,且实测单位径流深输沙模数越大,误差越大。
4) 滇西地区的龙川江、怒江、澜沧江3个流域中的12个水文站点,计算值与实测值相对误差小于30%的点占据75%。
站点 | 计算值 | 实测值 | 误差率 |
---|---|---|---|
景谷 | 1.1311 | 1.6127 | −29.87% |
勐海 | 0.4387 | 0.4665 | −5.96% |
勐省 | 0.7644 | 0.8619 | −11.32% |
孟连 | 0.9149 | 1.2058 | −24.13% |
羊庄坪 | 0.9021 | 0.8363 | 7.87% |
塘上 | 0.8273 | 0.9876 | −16.23% |
拉贺练 | 0.5365 | 0.4749 | 12.96% |
麻栗坝 | 0.6104 | 1.1979 | −49.04% |
等戛 | 0.7530 | 0.6830 | 10.25% |
大文 | 1.1626 | 1.6422 | −29.21% |
龙潭寨 | 1.5520 | 1.6445 | −5.62% |
南溪街 | 1.0219 | 1.0679 | −4.31% |
牛街 | 1.0936 | 0.7553 | 44.78% |
西山 | 0.5638 | 0.4511 | 24.98% |
表3. 实测值小于2 t/(km2∙mm)水文站点
由表3可以看出,除了麻栗坝和牛街两个水文站控制的流域,其余12个流域的计算值和实测值相对误差均在30%以内。分析模型各个因子数据,发现麻栗坝站的流域内坡耕地面积百分数仅为9.98,大于25˚耕地占坡耕地面积百分数为9.10,均在10%以下,与其他水文站的P、Pk值相比较小。而单位径流深输沙模数方程的构建是根据流域产沙机制建立,坡耕地面积百分数是泥沙来源的重要参数,麻栗坝站的P值、Pk值较低也导致计算结果比实测值更低。
牛街站位于云南省东北边界,集水流域超出了云南境内,而研究所用的数据只包括云南境内,没有境外数据,这是导致牛街站计算值误差较大的主要原因。此外,牛街站仅有84年的水文数据,而流域的输沙量、径流量均有年间差异,无法保证84年的水文数据是表征其多年平均水平,这也是导致牛街站的计算结果与实测单位径流深输沙模数差异较大的原因。
由于输沙模数 [
M S = R S ⋅ F = y S [
由此可得 S = M s / y
又由于 m s = M s / y ,故有
S = m s 。
其中,Ms——流域多年平均输沙模数(t/km2)。
R——多年平均径流量(m3/s)。
S——河流多年平均含沙量(kg/m3)。
F——流域面积(km2)。
y——流域多年平均径流深(mm)。
ms——单位径流深输沙模数(t/(km2∙mm)或kg/m3)。
根据输沙模数和单位径流深输沙模数的定义,发现河流多年平均含沙量与单位径流深输沙模数在数值上是相等的,因此ms的值为河流多年平均含沙量,其单位为kg/m3。再结合径流量数据,就可求得河流多年平均悬移质输沙量。
岔河水库位于大理市永平县杉阳镇境内倒流河上,岔河水库规模为中型水库,主要功能为农村生活和灌溉供水。岔河水库干流方案坝址位于倒流河干流上,距下游永和村约2 km,地理坐标东经99˚27′41.6″,北纬25˚15′20″ (1954年北京坐标系),坝址处高程1524 m,坝址以上控制径流面积79.5 km2,多年平均径流量3285万m3,P = 75%径流量2430万m3。
在ArcGIS环境下,以云南省DEM数据为基础,使用水文分析工具提取流域集水面积图,可得出岔河水库的集水面积。将集水面积图与DEM数据进行叠加裁剪,得到流域高程图(图2)。
图2. 岔河流域水系高程图
将流域高程图与云南省土地利用图进行叠加分析,可提取出流域内坡耕地和流域内大于25˚坡耕地分布图(图3),根据计算得出流域内坡耕地面积百分比及大于25˚耕地占坡耕地面积百分比。
图3. 岔河流域坡耕地分布图
将土壤分布图与流域集水面积图叠加分析,得到流域内不同类型土壤的分布面积,参考RUSLE模型中K因子计算方法,得到各因子的数据如表4。
影响因子 | F | P | PK | K | V |
---|---|---|---|---|---|
计算值 | 78.2 | 28.88 | 24.25 | 0.352 | 53.93 |
表4. 岔河水库泥沙淤积影响因子
将算得的F,P,PK,K,V因子的值带入多元方程 m s = K P 0.435 P K 0.025 F − 0.092 V − 0.05 中,求得ms (单位径流深输沙模数)等于1.71 t/(km2∙mm),根据前文分析,河流多年平均含沙量在数值上等于单位径流深输沙模数,因此求得干流方案中控制流域多年平均含沙量为2.04 kg/m3。根据实测资料,干流方案坝址多年平均来水量3285万m3,求得该坝址处多年平均悬移质输沙量为6.7万吨。泥沙推悬比参考多沙河流经验取值为30%,则该坝址处多年平均推移质输沙量2.01万吨,总输沙量为8.71万吨。
岔河站有2014~2016年实测泥沙资料,3年平均悬移质数沙量为4.59万吨,泥沙推悬比取经验值30%,总输沙量为5.97万吨,考虑区域土壤侵蚀研究期间(2014~2016)云南省处于旱季,与3年径流量对比,监测结果较多年平均值较小。
根据单位径流深输沙模数计算得出干流坝址多年平均总输沙量为8.93万吨,比这三年平均值大,符合实际情况。同时,通过对比周边水文站以及根据岔河水库断面的径流输沙能力进行判断,岔河水库多年平均总输沙量8.93万吨这个计算结果可信度较高。
建立多元回归方程所使用的P、Pk这两个因子数据是通过云南省DEM数据和云南省土地利用类型图叠加分析提取得到。由于DEM数据分辨率为30m,而实际生活中一块完整耕地面积很少有900 m2,因此提取出的数据与实际相比可能存在一定程度的误差。由于能够获取的数据有限,只能使用20个流域的数据建立多元回归模型,可以添加一定数量的水文站,使模型能够适用于整个云南省境内的河流输沙量计算。
将计算值与实测值进行对比,当实测单位径流深输沙模数小于2时,单位径流深输沙模数方程计算值与实测值误差小于30%的点占86%。而四川地区研究构建的模型 [
金沙江流域内的牛街站位于云南省东北,其集水区域超出了云南省境内。而土地利用、土壤类型等数据只有云南省境内的,因此只计算了这个流域在云南省境内的数据,导致计算结果误差较大。
1) 经过研究,得出计算云南省中小河流单位径流深输沙模数方程为 m s = K P 0.435 P K 0.025 F − 0.092 V − 0.05 。本方程由云南省境内资料建立,原则上只适用于云南地区。
2) 将模型计算结果与实测泥沙资料相比较,发现该方程在实测ms小于2的水文站的控制流域计算结果更精确,且误差小于30%的水文站占据86%。当单位径流深输沙模数大于2时,ms的实测值越高,计算结果误差越大。
3) 使用单位径流深输沙模数方程对永平流域进行输沙量计算,得到该流域输沙量为8.71万吨,与永平流域多年平均输沙量相近。
国家自然科学基金项目(编号:42061006,41561144012)。
傅开道,宁 梅,戴雅婷,史正涛,张天力,顾世祥. 云南省中小河流悬移质泥沙计算模型研究Study on Suspended Sediment Calculation Model of Medium and Small Rivers in Yunnan[J]. 地球科学前沿, 2021, 11(03): 247-256. https://doi.org/10.12677/AG.2021.113021