·31·
水 土 工 程中国农业文摘·农业工程 2021年第1期水源涵养是区域生态系统服务功能的重要构成，是水土保持工程效益的重要组成，也是区域生态环境直接影响人类生产生活的重要生态功能[1]。尤其在重要水源区和干旱缺水地区，水源涵养功能是保证区域植物落存活并充分发挥效益的基础条件，也是改善区域小气候、修复区域生态系统、保障区域生态安全、提升区域环境质量的重要指标。甘肃省是全国水土流失最严重的省区之一，水土流失面积大、范围广、类型多、强度高、危害重，严重制约着区域经济、社会发展。尤其是陇中区域干旱缺水，对于提升区域水源涵养的要求更为迫切。水源涵养功能是区域生态系统重要服务功能之一，是区域生态系统状况的重要指示器[2]。影响区域水源涵养功能发挥的因素很多，主要包括区域干旱程度、植被拦蓄能力、土壤入渗程度、降雨强度等，以及水分在各因素之间运移、转化的自然过程；水源涵养功能的变化也能反作用于区域气候、水文、土壤、植被等自然因素；直接影响区域生态服务功能的发挥。对于水源涵养的研究成果很多，各类评估方法也较为成熟，但对小流域尺度利用遥感数据结合实地测量，展开小流域尺度的水源
涵养动态评估的研究尚不多见。甘肃省中部的黄土丘陵沟壑区属于典型的半干旱区，区域内降水时空分布不均，水资源贫乏，缺水成为区域生态环境改善、社会及经济发展的主要制约因素，在小流域尺度上定量的评估甘肃省黄土丘陵沟壑区水土保持工程水源涵养功能的发挥，分析其影响因素及空间特征，对科学认识黄土丘陵沟壑区水土保持工程的水源涵养功能，合理保护和利用区域水土资源，展开黄河上中
游区域大保护，稳定和增强甘肃省生态屏障作用，制定区域生态环境保护规划、设计、政策等，具有十分重要的意义。
1 研究区概况1.1 区域概况
本文研究区位于国家水土保持重点治理区的甘肃省通渭县西北部，黄河二级支流渭河一级支流牛谷河流域上中游，属黄土丘陵沟壑区第三副区。牛谷河上游流域共包含有马营、水岔、李家大河、朱家营滩、段家峡、万家岔、蒋家川7条主要支流，研究区域内高差为747m，干沟长37.5km，干沟比降2.6%，主支沟长64.6km，主支沟平均比降3.7%，其它支毛沟长491.91km，沟壑密度1.43km/km 2。研究区地理位置介于东经104°54′41″-105°17′20″，
北纬35°10′08″-35°23′34″之间，海拔1 774- 2 521m，地貌为黄土梁峁沟壑地貌，沟壑占比达9.0%，小支毛沟多为“V”型沟，主支沟多为“U”型沟。研究区属中温带季风大陆性气候，多年平均
基于GIS的水土保持工程水源涵养评估
——以通渭县牛谷河项目区为例
徐剑春，王 博，罗进选
（甘肃省水土保持科学研究所，甘肃兰州 730020）
摘要：国家生态文明战略和黄河大保护战略全面实施，甘肃省水土保持工程建设发展迅速，在工程建设和水土保持监测高标准的要求下，研究基于GIS的监测评价方法在小流域尺度的应用具有很强的现实指导意义。现场调查结合遥感影像及成果，以水量平衡方程为理论指导，选择降雨、径流、蒸散发数据构建水源涵养价值评估体系，以文献分析法获取经验径流系数，利用ArcGIS计算不同因子水源涵养指数及其面积，分析水源涵养量。以通渭县牛谷河研究区为例，根据项目实施前后的对比，探索小尺度水土保持工程GIS监测评价的方法，评估区域水源涵养价值及其动态变化；分析结果表明牛谷河项目区水源涵养量从6 982.55万m 3增加至8 131.58万m 3，增加了16.46%。
关键词：水土保持工程；水源涵养；GIS；评估
课题项目：重大生态工程水土保持监测技术研究，甘水
科外【2017】76号，甘水水保发【2016】204号。
作者简介：徐剑春，硕士，高级工程师，主要从事黄土丘
陵沟壑区水土流失综合治理及水土保持监测技术研究。
·32·水 土 工 程中国农业文摘·农业工程 2021年第1期
≥10℃的活动积温1 226℃，年均气温6.6℃，年平均太阳辐射总量为129kj/（cm2.a），多年平均日照时数2 239h，多年平均无霜期131d，多年平均大风日数为32天。研究区内包含42个行政村，人口密度166人/ km2，人均耕地0.37hm2，涉及通渭县的华岭镇、马营镇、北城铺镇、三铺乡、平襄镇、陇阳镇共5镇1乡。
1.2 建设概况
土加于近年来研究区内新增建设各类水土保持措施150km2，其中建设梯田4 007hm2，水土保持林2 085.4hm2，水土保持种草1 664hm2，封禁治理 7 017.11hm2，新建土谷坊210座，水窖96眼，沟头防护4.24km。
在地形坡度小于15°距村庄较近的坡耕地中建设水平梯田；坡度15°-25°距村庄较远立地条件和水分条件较好的坡耕地中退耕种植紫花苜蓿；坡度＞25°的陡坡耕地及15°-25°植被稀疏的荒山荒沟，水平阶整地后种植侧柏及沙棘；陡坡耕地及荒坡，鱼鳞坑整地后种植山杏、柠条、沙棘；对于坡度＞35°且水土流失严重的荒坡荒沟以及不宜布设其它防治措施的荒地采用封禁措施；根据各流域具体地形条件和沟道特征，选择适宜的地点全面建设谷坊及沟头防护措施。
2 评价方法
2.1 评价指标及模型
水源涵养是陆地生态系统对降雨进行分配的一个过程，通过土壤的入渗、微地形的截留、植物拦蓄与蒸腾、地表蒸发等，再利用区域生态系统特定结构，实现降雨、水流、水循环的调控来影响水源涵养量。水源涵养的生态功能主要体现在消减降雨侵蚀力，减轻土壤溅蚀，缓和径流形成，拉长行洪时间，降低洪水危害，增加土壤含水量，补充地下水，滞洪补枯、降低径流波动，保证水质等方面。陆地生态系统中影响水源涵养的因素很多，主要有地表覆盖因子、降雨强度因子、地形地貌因子、农业耕作方式、土壤植被类型等。在本文中以水源涵养量作为水土保持工程水源涵养功能的评估指标，采用水量平衡方程的基本原理来计算区域水土保持工程的水源涵养量。评估模型为：
（1）
式中：TQ为水源涵养量（m3），P i为第i区域年均降雨量（mm），R i为第i区域年均地表径流量（mm），ET i为第i区域年均蒸散发量（mm），A i为第i区域面积（km2），i为研究区不同影响因子综合分析形成数据集中的第i个单元，j为研究区影响因子数据集的单元数。
2.2 数据来源
根据上述水量平衡方程以及各类影响因子的分析，进行水土保持工程水源涵养评估所需数据主要有：气象数据集（年均降雨量、雨强、气象站分布资料等）、蒸散发数据集（研究区年均蒸散发量分布栅格图）、土壤数据集（土壤类型及分布、主要土壤的理化性质、实测土壤绘制的栅格图、现场调查的
土地利用类型分布）、高程数据集（DEM数据及治理完成区域高程测量资料）、遥感影像等，具体计算中在收集现有遥感资料的基础上，对遥感数据尺度较大的采用现场实测数据加密后进行插值计算。
（1）降雨量因子
研究区附近有2处气象观测站，华家岭气象站纬度35.38°经度104.83°海拔2 450.6m，通渭气象站纬度35.22°经度105.23°海拔1 768.2m，系列资料均能够满足评价所需。降雨量因子，利用2处气象站观测的年均数据，采用ArcGIS软件，插值计算绘制降雨量因子（P值），栅格分布图（见图1）。
图1 降雨量因子（P）插值分布图
（2）地表径流因子
地表径流因子利用地表径流量来进行分析，在本文中采用经验值计算径流量，即由降雨量乘以地表径流系数获得，计算公式如下：
R=P×α （2）
式中：R为地表径流量（mm），P
为多年平均降
·33·
水 土 工 程中国农业文摘·农业工程 2021年第1期
雨量（mm），α为该区域平均地表径流系数。
平均地表径流系数主要依据现场调查的治理前后植被类型区、土地利用类型，在ArcGIS中绘制形成要素图，结合研究区NDVI数据，筛选通过文献分析法获取的各类地表径流经验值（见表1），将经验值赋值在绘制的要素图中，转制形成地表径流栅格图。
表1 各类型生态系统地表径流系数均值表
生态系统类型1生态系统类型2平均地表径流系数（%）
森林
常绿阔叶林 2.67常绿针叶林
3.02针阔混交林 2.29落叶阔叶林 1.33落叶针叶林0.88稀疏林19.20灌丛
常绿阔叶灌丛
4.26落叶阔叶灌丛 4.17针叶灌丛 4.17稀疏灌丛19.20草地
草甸
8.20草原 4.78草丛9.37稀疏草地
18.27
在ArcGIS中利用栅格计算器，根据公式（2）及降雨量栅格图，绘制治理前后的地表径流因子（R 值）栅格分布图。
图2 治理前径流系数（α）分布图
图3 治理后径流系数（α）分布图
（3）蒸散发因子
蒸散发因子根据国家生态系统观测研究网络科技资源服务系统网站提供的产品数据。原始数据空间分辨率为1km，通过ArcGIS 软件对栅格图像重采样后得到小流域尺度的蒸散发因子（ET）栅格分布图（见图4）。
图4 项目区蒸散发量（ET）分布图
（4）面积因子
面积因子利用治理前后的水源涵养指数（P-R-ET）栅格分布图，在ArcGIS中计算获取不同指数单元的面积。
2.3 水源涵养能力计算
将降雨量因子（P ）、地表径流因子（R ）、蒸散发因子（ET ）3类评价因子栅格图，统一形成10m 分辨率的栅格数据，在ArcGIS中利用栅格计算器，根据评价模型P-R-ET计算得到水源涵养指数（P-R-ET）栅格分布图（见图5、图6）。在计算得出的栅格图中提取不同模块的面积数据，根据评估模型加权计算治理前后的水源涵养量。3 评价结果
根据以上得出的治理前后水源涵养量（TQ）分
布图（见图7、图8），分析牛谷河项目区水土保持工程实施前后水源涵养能力的变化（见表2），造成项目区水源涵养量变化的主要原因为通过各类水土保持措施的实施，改变了地表起伏度（坡度、坡长），增加了项目区植被覆盖系数，从而降低项目区径流系数，提高水源涵养力。
从表2得出牛谷河区域水源涵养量从6 982.55万m 3增加至8 131.58万m 3，增加了1 149.03万m 3，增加了16.46%。
4 结语
基于遥感及GIS技术、物联网监测技术进行生态
·34·水 土 工 程中国农业文摘·农业工程 2021年第1期
Ⲵ䶒〟ᮠᦞˈṩᦞ䇴ՠ⁑ර࣐ᵳ䇑㇇⋫⨶ࡽਾⲴ≤Ⓚ⏥ޫ䟿DŽ
图5 治理前水源涵养指数（P-R-ET）分布图 图6 治理后水源涵养指数（P-R-ET）分布图
图7 治理前水源涵养量（TQ）分布图 图8 治理后水源涵养量（TQ）分布图
表2 治理前后项目区水源涵养变化分析表
分析因子
治理前治理后
栅格数（个）面积（km2）百分比（%）栅格数（个）面积（km2）百分比（%）
地形起伏度
≤0.2585 70058.577.53 1 479 700147.9719.04 0.2-0.4 2 637 600263.7633.93 2 028 200202.8226.09 0.4-0.5 2 087 100208.7126.85 1 806 100180.6123.23 0.5-0.6 1 366 800136.6817.58 1 364 400136.4417.55≥0.6 1 096 200109.6214.10 1 095 100109.5114.09
植被覆盖度
≥0.6 3 616 200361.6246.52 3 646 200364.6246.91 0.4-0.5 2 281 700228.1729.35 2 854 300285.4336.72 0.3-0.4 1 596 200159.6220.53 1 163 500116.3514.97 0.2-0.368 200 6.820.8831 200 3.120.40≤0.2211 10021.11 2.7278 3007.83 1.01
水源涵养指数
≤50 1 440 500144.0518.53904 40090.4411.63 50-100 3 729 600372.9647.98 2 988 200298.8238.44 100-150 1 978 500197.8525.45 2 932 900293.2937.73 150-200617 70061.777.95880 50088.0511.33≥20072000.720.0945 200 4.520.58
水源涵养量（万m3） 6 982.558 131.58
·35·
水 土 工 程中国农业文摘·农业工程 2021年第1期
项目监测，是水土保持监测行业未来发展的方向，探索小流域尺度的GIS监测评价体系及方法，研究现场实测结合遥感分析的数据获取及处理技术，具有很强的现实指导意义。本文尝试利用GIS技术针对具体治理项目展开监测，但还存在遥感数据精度不高，现场气象监测等数据来源较少，通用计算模型与区域实际不完全匹配等一些问题，还需后续进一步展开相关研究，逐步完善甘肃省水土保持生态建设项目的监测体系和方法。参考文献
[1] 郭灵辉，高江波，王天平，等.太行山典型森林水源
涵养分层特征与差异研究[J].人民黄河，2018，40（8）：91-95.
[2] 龚诗涵，肖洋，郑华，等.中国生态系统水源涵养空
间特征及其影响因素[J].生态学报，2017，37（7）：2455-2462.
[3] 贾立，胡光成，郑超磊，等.中国-东盟1km分辨率
地表蒸散发数据集（2013）[J].全球变化数据学
报，2017，1（3）：282-289. DOI：10.3974/geodp. 2017.03.05.
[4] 郑超磊，贾立，胡光成.全球陆表实际蒸散发数据
集（2013-2014）.国家青藏高原科学数据中心，2019. doi：10.11888/Hydro.tpdc.270298.[5] 莫小娥，张维江，梁如心.海子流域土壤物理特性及
其水源涵养能力分析[J].人民黄河，2018，40（1）：71-75.
[6] 郑飞龙，谭军利，王西娜，等.黄土高原水源涵
养林截蓄降雨研究进展[EB/OL].人民黄河，2019-08-19. http：//knski/kcms/detail/41.1128.TV.20190819.1117.006.html.[7] 包玉斌，李婷，柳辉，等.基于InVEST模型的陕北
黄土高原水源涵养功能时空变化[J].地理研究，2016，35（4）：664-676.
[8] 罗娅，杨胜天，刘晓燕，等.黄土高原次降雨径流产
沙变化及其与雨强和植被盖度变化的相关性[J].干旱区研究，2015，32（4）：698-709.
氟环唑2 000倍液喷雾。以上药物都可以有效的这三种病症。当然经济林木幼年树的虫害问题也是较为严重的，像红蜘蛛、锈壁虱、潜叶蛾等。这些害虫会严重危害经济林木嫁接的成功率，所以需要重点防治。种植户在以虫治虫，以菌防虫，防治结合的基础上，红蜘蛛可以使用24%螺螨酯4 000倍液喷雾防治；而潜叶蛾则是使用1.8%阿维菌素1 000倍液喷雾防治[6]。5 总结
总的来说，我国的经济林木嫁接以及培育种植技术十分成熟，只要技术管理到位，经济林木的嫁接成活率以及幼年树的生长发育都会较为喜人，当然这一过程也是需要做好嫁接处理以及嫁接后的防虫防病，而且在相关的培植过程中，还可以借助多种栽培手段保障经济林木幼年树的成长和品质，进而促进我
国经济林木产业的健康发展和进步。参考文献
[1] 陈文汉.我国经济林木种植存在的问题及对策[J].
江西农业，2018（2）：95.
[2] 栾艳波.经济林木嫁接技术及管理[J].现代农业，
2016（8）：75-76.
[3] 周波.经济林木嫁接技术及管理[J].科技展望，
2016，26（3）：88.
[4] 闫淑英，王瑗，迟卫京.收益现值法在经济林木价值
损失鉴定中的应用探讨[J].内蒙古林业调查设计，2019，42（4）：86-87.
[5] 姚永泉，刘永辉.老桑树截干嫁接复壮及快速成林
管理技术[J].蚕桑茶叶通讯，2015（4）：19-20.[6] 刘润元.甘肃庆阳市干果经济林优良品种引种研究
报告[J].中国园艺文摘，2018，34（1）：211-212.
（上接第16页）
版权声明：本站内容均来自互联网，仅供演示用，请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系QQ:729038198，我们将在24小时内删除。
发表评论