CN114580910A - 一种考虑道路影响的流域泥沙连通性评价方法 - Google Patents
一种考虑道路影响的流域泥沙连通性评价方法 Download PDFInfo
- Publication number
- CN114580910A CN114580910A CN202210211774.3A CN202210211774A CN114580910A CN 114580910 A CN114580910 A CN 114580910A CN 202210211774 A CN202210211774 A CN 202210211774A CN 114580910 A CN114580910 A CN 114580910A
- Authority
- CN
- China
- Prior art keywords
- grid
- tool
- road
- connectivity
- calculating
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 239000013049 sediment Substances 0.000 title claims abstract description 41
- 238000011156 evaluation Methods 0.000 title claims abstract description 8
- 238000000034 method Methods 0.000 claims abstract description 63
- 238000004364 calculation method Methods 0.000 claims abstract description 32
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 51
- 230000008569 process Effects 0.000 claims description 33
- 238000009825 accumulation Methods 0.000 claims description 18
- 230000005540 biological transmission Effects 0.000 claims description 15
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000011160 research Methods 0.000 claims description 11
- 230000001174 ascending effect Effects 0.000 claims description 9
- 238000012732 spatial analysis Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 7
- 230000000007 visual effect Effects 0.000 claims description 4
- 230000003746 surface roughness Effects 0.000 claims description 3
- 239000002689 soil Substances 0.000 abstract description 9
- 230000007547 defect Effects 0.000 abstract description 5
- 230000008859 change Effects 0.000 description 5
- 238000004162 soil erosion Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 239000004576 sand Substances 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 230000008021 deposition Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 238000012854 evaluation process Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 239000003309 forest litter Substances 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000001556 precipitation Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
- 238000009418 renovation Methods 0.000 description 1
- 238000005211 surface analysis Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Educational Administration (AREA)
- Economics (AREA)
- Development Economics (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Primary Health Care (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明属于生态水文技术领域，公开一种考虑道路影响的流域泥沙连通性评价方法，该方法通过改进连通性指数IC的计算，即在计算过程中以能更好地反映山区森林区域的泥沙连通性的相对平滑度指数RS作为权重因子，同时，将道路作为掩膜图层参与IC的计算过程，从而考虑路网影响下河网掩膜图层的生成，进而反映道路网络对地表径流方向的影响、对汇流路径长度的改变等，以弥补道路作为排水通道对泥沙连通性影响的不足，使其适用于道路影响的山区森林区域。本发明可以弥补道路作为排水通道对泥沙连通性影响的不足，使得计算结果既能反映山区森林流域的土壤覆盖与地形特征，又能充分反映道路网络对流域泥沙连通性的影响。
Description
技术领域
本发明属于生态水文技术领域，尤其涉及一种考虑道路影响的流域泥沙连通性评价方法。
背景技术
泥沙连通性是指流域不同地貌或景观单元间的泥沙级联关系，用于表征流域内泥沙克服区域阻力通过分离和输移从源到汇区(从坡面到沟道)的运输程度，从而衡量泥沙在地貌单元运动的潜能。这既取决于地貌或景观单元的空间格局，也受控于气候、地貌、生态、水文、土壤侵蚀等地表过程及其时空变异。
泥沙连通性是近十年来国际上的研究热点，与地形地貌、生态水文、土壤侵蚀及水质保护等多个研究领域密切相关，它强调不同地貌或景观单元间的泥沙输移、分离和沉积等过程，是控制流域景观演变的最重要属性之一，可用于理解土壤侵蚀与泥沙输移过程及其时空分布。泥沙连通性重点关注流域内不同地貌或景观单元间的泥沙级联关系，是决定流域泥沙输移方式、效率和规模的重要因素，是探究泥沙来源、识别水土保持重点区域的基础，在径流产生、泥沙运移相关的水文过程中起着重要作用。因此，开展流域泥沙连通性研究将有助于揭示流域水沙运移过程变化的影响机制，进而为流域生态修复提供科学依据。
近年来，由于人类活动强度加剧，原生地带性稳定群落退化严重、森林资源锐减、水土流失问题日趋严峻。特别是，受森林采伐、旅游开发、矿产开采等因素的影响，道路网络持续扩张，如对原有林区采伐废旧道路的恢复、现有道路的翻修和新规划道路的建设等，进一步加剧了水土流失、改变了流域原有的水沙输移过程。泥沙连通性作为识别与量化泥沙输移路径及其时空变化的基础，在探究泥沙来源热点、识别水土流失治理重点区域等方面发挥着重要作用，且与地貌演变、水文过程、景观生态、土壤侵蚀、水质演变及水资源管理等密切相关。因此，研究泥沙连通性具有重要的理论意义与应用价值。
野外调查法作为定性的研究方法，需要在发生明显的次降水或侵蚀事件后在流域内调查土壤侵蚀及泥沙沉积等情况，且需要详细的遥感影像、航片和土地利用等资料，但定性的判断只能用来评价或解释定量化指标的计算结果。
近年来，随着科技手段的不断进步，能表征连通性的数据的分辨率得到大大的提升。基于高精度数据，量化泥沙连通性的方法亦逐渐被提出，大致分为图论法、指标法和模型模拟法。如Heckmann等人应用图论法分析了阿尔卑斯山脉中部流域泥沙级联及其路径的空间分布。高海东等人基于图论的原理研究了淤地坝如何对径流造成影响。这种方法给泥沙通道的网络结构以及上下坡泥沙连通性的分析提供了极大的帮助。然而，它只能适用于泥沙连通性的空间格局分析，对输沙能力的分析较为缺乏。指标法中Broselli等人开发的连通性指数(IC)应用最为广泛，如Foerster等将其用于评估西班牙土地利用变化对连通性的影响。模型模拟法中，Lesschen等人应用LAPSUS模型发现植被斑块和农业梯田的空间分布是影响流域泥沙连通性的主要影响因素。但是这些方法在针对道路影响下泥沙连通性的空间分布及其变化的影响尚不明确。尤其是，道路网络在改变泥沙传输方向的同时捕获部分泥沙和径流，使其沿道路延伸方向向下游传输，并在道路与河流交汇点向流域出水口传输，进而改变原有的泥沙传输路径。
综上，图论法针对不同的研究对象需要选择不同的定量指标进行表达，虽能直观展示不同景观单元间泥沙级联的空间分布与动态变化，但适用的空间尺度有限，在较大尺度的流域上对小区、坡面、切沟和细沟等较小景观单元的连通性表达能力有限。指标法和模型模拟法都可用于模拟流域侵蚀泥沙源和汇，定量化表征泥沙连通性，尤其是基于连通性指数的指标法是目前应用最为广泛的研究方法。其中，连通性指数能表征流域泥沙源-汇间的潜在联系。但上述技术在评估流域泥沙连通性时均较少考虑道路作为传输通道对泥沙传输方向、汇流路径长度、有效汇流面积等的影响。
发明内容
本发明针对泥沙连通性评价过程中，现有方法未考虑适应山区森林区域的土地覆盖和地表特征，以及道路作为传输通道对流域流向及径流与泥沙传输路径的影响，从而导致泥沙连通性评价过程中道路影响的准确性差的问题，提出一种考虑道路影响的流域泥沙连通性评价方法，基于指标法中的连通性指数(IC)，在计算过程中考虑适用于山区森林流域的相对平滑度指数(RS)作为权重因子；考虑道路作为传输路径对流向的影响，并在影响流域河网的基础上影响汇流路径长度、上坡有效集水区面积等，使得计算结果既能反映山区森林流域的土壤覆盖与地形特征，又能充分反映道路网络对流域泥沙连通性的影响。
为了实现上述目的，本发明采用以下技术方案：
一种考虑道路影响的流域泥沙连通性评价方法，包括：
步骤1、基于ArcGIS平台，进行数据的获取和处理；所述数据包括数字高程模型DEM数据、谷歌卫星影像数据、根据谷歌卫星影像数据目视解译得到的流域内道路网络与河流网络；
步骤2、利用ArcGIS平台，基于流域内DEM数据，运用空间分析中的坡度工具计算得到研究区坡度因子；
步骤3、计算相对平滑度指数作为权重因子，以反映山地森林地区土地覆盖和地表特征对泥沙连通性的影响；
步骤4、以谷歌卫星影像数据矢量化得到的流域内道路网络为基础，基于ArcGIS平台，运用要素转栅格工具和重分类工具区分流域内的道路和其它土地利用类型，并将道路栅格赋值为0，其它土地利用类型赋值为1，进而生成道路掩膜图层；
步骤5、以流域内DEM数据为基础，基于ArcGIS平台，获取水流流向与汇流累积量图层；
步骤6、通过空间分析中的水流长度工具计算沿流路径从每个像元到栅格边上的汇点或出水口的下坡距离，并以流域泥沙传输的最终流向作为水流流向的栅格数据、
作为权重栅格数据，从而得到下坡组分；
步骤7、借助汇流累积量工具分别以相对平滑度指数和坡度为权重栅格数据，计算不同权重栅格数据下的汇流累积量，根据不同权重栅格数据下的汇流累积量，通过栅格计算器工具计算上坡平均权重因子和上坡平均坡度，从而得到上坡组分；
步骤8、通过栅格计算器工具中的Log10函数，计算得到泥沙连通性指数。
进一步地，所述步骤1中按照如下方式得到流域内道路网络与河流网络：
利用ArcGIS平台，基于2m分辨率的谷歌卫星影像数据，在Catalog下新建要素类型为折线的矢量图层，调用编辑器在编辑命令下描绘谷歌卫星影像数据中呈现的所有道路，描绘完成后点击编辑器下的“保存编辑内容”与“停止编辑”，即生成道路网络；重复上述步骤，生成河流网络。
进一步地，所述步骤2还包括：
通过栅格计算器工具为坡度＝0的栅格字段赋值为0.005。
进一步地，所述步骤3包括：
以流域内DEM为基础，通过Spatial Analyst Tools中的表面分析-坡度，计算流域内的坡度图层；随后，借助栅格计算器，通过地面粗糙度的计算公式，即
得到粗糙度图层；然后，根据粗糙度指数与曼宁系数的关系式，计算得到曼宁系数：
Y＝2.4234×e0.3005ln(n)
式中：Y为地表粗糙度，n为曼宁系数；
随后，按照下式得到相对平滑度指数：
式中：nmin为曼宁系数的最小值；RS为相对平滑度指数，是一个无量纲数字，从0.0～1.0不等。
进一步地，所述步骤5包括：
步骤5.1、以流域DEM数据为基础，利用ArcGIS平台的填洼工具先对DEM进行填洼处理，随后运用空间分析工具中的水流流向工具提取DEM填洼后的水流方向，即FlowDir1；
步骤5.2、通过栅格计算器工具，将道路掩膜图层叠加到起始FlowDir1上，即道路掩膜*FlowDir1，从而得到路网影响下的水流流向，即FlowDir2；
步骤5.3、基于路网影响下的FlowDir2，运用空间分析中的汇流累积量工具计算栅格单元所流经的水量大小，即汇流累积量FlowACC1；
步骤5.4、根据FlowACC1，通过栅格计算器工具，结合矢量化河网确定阈值，从而提取流域河网掩膜图层；再次通过栅格计算器，将河网掩膜图层叠加在FlowDir2上，得到流域泥沙传输的最终流向，即FlowDir3；
步骤5.5、利用汇流累积量工具得到水流流向FowDir3的汇流累积量，即FlowACC2；利用栅格计算器，使FlowACC2每个栅格的汇流累积量值+1，得到最终的汇流累积量图层，即FlowACC3。
进一步地，所述步骤6包括：
通过空间分析中的水流长度工具计算沿流路径从每个像元到栅格边上的汇点或出水口的下坡距离，并以FlowDir3作为水流流向的栅格数据、
作为权重栅格数据，得到下坡组分，具体公式如下：
式中：Ddn表示径流和泥沙沿路线到达指定流域出水口或沟道的可能性；di为研究单元到泥沙汇的距离；Wi是第i个单元的权重因子，此处用相对平滑度指数代替；Si是第i个单元的坡度因子；
借助栅格计算器工具为上述计算结果为0的栅格单元赋值。
进一步地，所述步骤7包括：
借助汇流累积量工具分别以相对平滑度指数和坡度为权重栅格数据，计算不同权重栅格数据下的汇流累积量，即FlowRS和FlowS；随后根据FlowRS和FlowS，通过栅格计算器工具计算上坡平均权重因子和上坡平均坡度，并根据下式、借助栅格计算器，计算得到上坡组分：
式中：Dup表示上坡的连通分量，决定了上坡产生的径流和泥沙向下游流动的潜力；
是上坡贡献面积的平均权重因子，代表山区森林流域径流和泥沙在输移过程中可能遇到的阻抗；
为上坡贡献区域面积的平均坡度；
为上坡的贡献面积。
与现有技术相比，本发明具有的有益效果：
本发明通过改进连通性指数(IC)的计算，即在计算过程中以能更好地反映山区森林区域的泥沙连通性的相对平滑度指数(RS)作为权重因子，同时，将道路作为“掩膜”图层参与IC的计算过程，从而考虑路网影响下“河网掩膜”图层的生成，进而反映道路网络对地表径流方向的影响、对汇流路径长度的改变等，以弥补道路作为排水通道对泥沙连通性影响的不足，使其适用于道路影响的山区森林区域。本发明可以弥补道路作为排水通道对泥沙连通性影响的不足，使得计算结果既能反映山区森林流域的土壤覆盖与地形特征，又能充分反映道路网络对流域泥沙连通性的影响。
附图说明
图1为本发明实施例一种考虑道路影响的流域泥沙连通性评价方法的整体流程图；
图2为本发明实施例计算生成的坡度因子图；
图3为本发明实施例计算生成的相对平滑度指数图；
图4为本发明实施例计算生成的“道路掩膜”图；
图5为本发明实施例计算生成的道路影响下的“水流流向”图；
图6为本发明实施例计算生成的道路影响下的“汇流累积量”图；
图7为本发明实施例计算生成的下坡组分图；
图8为本发明实施例计算生成的平均权重因子图；
图9为本发明实施例计算生成的平均坡度因子图；
图10为本发明实施例计算生成的上坡组分图；
图11为本发明实施例计算生成的泥沙连通性指数图。
具体实施方式
下面结合附图和具体的实施例对本发明做进一步的解释说明：
如图1所示，一种考虑道路影响的流域泥沙连通性评价方法，包括：
步骤1、基于ArcGIS平台，进行数据的获取和处理；所述数据包括数字高程模型DEM数据、谷歌卫星影像数据、根据谷歌卫星影像数据目视解译得到的流域内道路网络与河流网络；
步骤2、利用ArcGIS平台，基于流域内DEM数据，运用空间分析中的坡度工具计算得到研究区坡度因子；
步骤3、计算相对平滑度指数作为权重因子，以反映山地森林地区土地覆盖和地表特征对泥沙连通性的影响；
步骤4、以谷歌卫星影像数据矢量化得到的流域内道路网络为基础，基于ArcGIS平台，运用要素转栅格工具和重分类工具区分流域内的道路和其它土地利用类型，并将道路栅格赋值为0，其它土地利用类型赋值为1，进而生成道路掩膜图层；
步骤5、以流域内DEM数据为基础，基于ArcGIS平台，获取水流流向与汇流累积量图层；
步骤6、通过空间分析中的水流长度工具计算沿流路径从每个像元到栅格边上的汇点或出水口的下坡距离，并以流域泥沙传输的最终流向作为水流流向的栅格数据、
作为权重栅格数据，从而得到下坡组分；
步骤7、借助汇流累积量工具分别以相对平滑度指数和坡度为权重栅格数据，计算不同权重栅格数据下的汇流累积量，根据不同权重栅格数据下的汇流累积量，通过栅格计算器工具计算上坡平均权重因子和上坡平均坡度，从而得到上坡组分；
步骤8、通过栅格计算器工具中的Log10函数，计算得到泥沙连通性指数。
进一步地，步骤1中数据的获取与处理包括：
本发明所使用的数据包括：15m分辨率的数字高程模型(DEM)数据、2m谷歌卫星影像数据、根据谷歌卫星影像数据目视解译得到的流域内道路网络与河流网络。
流域内道路网络与河流网络的提取步骤包括：利用ArcGIS平台，基于2m分辨率的谷歌卫星影像数据，在Catalog下新建要素类型为折线的矢量图层，调用“编辑器”在编辑命令下手工描绘谷歌卫星影像数据中呈现的所有道路，描绘完成后点击“编辑器”下的“保存编辑内容”与“停止编辑”，即生成“道路”网络。重复上述步骤，生成“河流”网络。
进一步地，步骤2中坡度因子(S)的计算过程包括：
以流域内高程数据DEM为基础，基于ArcGIS平台，运用空间分析中的坡度工具计算得到研究区“坡度”因子。为避免最小值为0时引起的连通性指数在计算过程中出现不收敛的情况，通过栅格计算器工具为“坡度”＝“0”的栅格字段赋值为0.005，即[(“坡度”＝＝0)*0.05+“坡度”](如图2)。
进一步地，步骤3中权重因子(W)的计算过程包括：
权重因子是反映流域地表状况和土地利用的加权系数，直接决定泥沙输送的过程和效率。本发明选择相对平滑度指数(RS)作为权重因子，以反映山地森林地区土地覆盖和地表特征对泥沙连通性的影响。该指数由曼宁系数推导出来的，同时保留泥沙连通性指数(IC)中权重因子(W)在计算过程的无量纲特性。
计算过程以流域DEM为基础，通过“Spatial Analyst Tools”中的“表面分析”-“坡度”，计算流域内的“坡度图层”。随后，借助“栅格计算器”，通过地面粗糙度的计算公式，即
得到“粗糙度图层”。然后，根据粗糙度指数与曼宁系数的关系式，计算得到曼宁系数。具体公式如下：
Y＝2.4234×e0.3005ln(n)
式中：Y为地表粗糙度，n为曼宁系数。
随后，借助相对平滑度指数与曼宁系数之间的关系式，计算得到相对平滑度指数RS(如图3)。具体公式如下：
式中：n为曼宁系数；nmin为曼宁系数的最小值；RS是一个无量纲数字，从0.0～1.0不等，保留了IC计算的无量纲性。这里使用RS作为权重因子，能够反映森林土壤表面积累的森林凋落物对径流与泥沙传输的影响。
综上，本发明对权重因子进行了改进：原有泥沙连通性指数在计算过程中侧重于考虑植被覆盖管理因子(C因子)对农业和森林环境中径流和泥沙通量的量化，本发明基于对山区森林环境的考虑，选择相对平滑度指数(RS)对权重因子加以改进，改进后的权重因子能更好地反映土地覆盖和地表特征对山区森林流域泥沙连通性的影响。
进一步地，步骤4中“道路掩膜”图层的获取过程包括：
以谷歌卫星影像数据矢量化得到的道路网络为基础，基于ArcGIS平台，运用要素转栅格工具和重分类工具区分流域内的道路和其它土地利用类型，并给道路栅格赋值为0，其它土地利用类型赋值为1，进而生成“道路掩膜”图层(图4)。
进一步地，步骤5中“水流流向”与“汇流累积量”图层的获取过程包括：
步骤5.1、以流域数字高程DEM为基础，为避免计算过程中伪地形对径流路径造成的影响，在数据处理时，利用ArcGIS平台的填洼工具先对DEM进行填洼处理。随后运用空间分析工具中的水流流向工具提取DEM填洼后的水流方向，即“FlowDir1”。
步骤5.2、通过栅格计算器工具，将“道路掩膜”图层叠加到起始“FlowDir1”上，即“道路掩膜*FlowDir1”，从而得到路网影响下的水流流向，即“FlowDir2”。
步骤5.3、基于路网影响下的“FlowDir2”，运用空间分析中的汇流累积量工具计算栅格单元所流经的水量大小，即汇流累积量“FlowACC1”。
步骤5.4、根据“FlowACC1”，通过栅格计算器工具，结合矢量化河网确定阈值，从而提取流域“河网掩膜”图层。再次通过栅格计算器，将“河网掩膜”图层叠加在“FlowDir2”上，得到流域泥沙传输的最终流向，即“FlowDir3”(图5)。
步骤5.5、利用汇流累积量工具得到水流流向“FowDir3”的“汇流累积量”，即“FlowACC2”。为避免在后续计算过程中出现零和无穷大的情况，利用栅格计算器，使“FlowACC2”每个栅格的“汇流累积量”值做“+1”处理，得到最终的“汇流累积量”图层，即“FlowACC3”(图6)。
综上，本发明对径流与泥沙传输通道进行了改进：通过“道路掩膜”图层的叠加，改变径流与泥沙的传输方向，在此基础之上生成“河网掩膜”图层，影响径流与泥沙从源到汇的传输路径长度，使得到达各个汇区的汇流累积量有所不同。在此过程中，道路的添加如何改变流向、汇流路径长度、汇流累积量，进而影响河网是泥沙连通性计算过程的关键点。这主要包括道路网络的栅格化，并在计算过程中将其考虑到道路影响下径流和泥沙沿着道路方向传输，以及在道路-河流交汇处向流域出水口的输移。
进一步地，步骤6中下坡组分(Ddn)的计算过程包括：
通过空间分析中的水流长度工具计算沿流路径从每个像元到栅格边上的汇点或出水口的下坡距离，并以“FlowDir3”作为水流流向的栅格数据、
作为栅格数据的权重栅格数据，得到“Flowlength”，此时即得到下坡组分，具体公式如下：
式中：Ddn表示径流和泥沙沿路线到达指定流域出水口或沟道的可能性；di为研究单元到泥沙汇的距离；Wi是第i个单元的权重因子，此处用相对平滑度指数(RS因子)代替；Si是第i个单元的坡度因子。
为避免后续计算过程中出现无穷大等不收敛的情况以及更准确的表达道路所经处的实际情况，借助栅格计算器工具为上述计算结果为0的栅格单元赋值，最终计算得到下坡组分(Ddn)(图7)。
进一步地，步骤7中上坡组分(Dup)的计算过程包括：
为计算上坡平均权重因子
和平均坡度
借助汇流累积量工具分别以“相对平滑度指数”和“坡度”为权重栅格数据，计算不同权重栅格数据下的“汇流累积量”，即FlowRS和FlowS，此时的“汇流累积量”分别受“相对平滑度指数”以及“坡度”的影响。随后根据FlowRS和FlowS，通过栅格计算器工具计算
(图8)和
(图9)，并根据上坡组分的计算公式(具体公式如下)借助“栅格计算器”，计算得到Dup(图10)。
式中：Dup表示上坡的连通分量，决定了上坡产生的径流和泥沙向下游流动的潜力；
是上坡贡献面积的平均权重因子，代表山区森林流域径流和泥沙在输移过程中可能遇到的阻抗；
为上坡贡献区域面积的平均坡度；
为上坡的贡献面积。
进一步地，步骤8中泥沙连通性指数(IC)的计算过程包括：
通过栅格计算器工具中的“Log10”函数，计算得到泥沙连通性指数(图11)，具体公式如下：
综上，本发明通过改进连通性指数(IC)的计算，即在计算过程中以能更好地反映山区森林区域的泥沙连通性的相对平滑度指数(RS)作为权重因子，同时，将道路作为“掩膜”图层参与IC的计算过程，从而考虑路网影响下“河网掩膜”图层的生成，进而反映道路网络对地表径流方向的影响、对汇流路径长度的改变等，以弥补道路作为排水通道对泥沙连通性影响的不足，使其适用于道路影响的山区森林区域。本发明可以弥补道路作为排水通道对泥沙连通性影响的不足，使得计算结果既能反映山区森林流域的土壤覆盖与地形特征，又能充分反映道路网络对流域泥沙连通性的影响。
以上所示仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (7)
1.一种考虑道路影响的流域泥沙连通性评价方法，其特征在于，包括：
步骤1、基于ArcGIS平台，进行数据的获取和处理；所述数据包括数字高程模型DEM数据、谷歌卫星影像数据、根据谷歌卫星影像数据目视解译得到的流域内道路网络与河流网络；
步骤2、利用ArcGIS平台，基于流域内DEM数据，运用空间分析中的坡度工具计算得到研究区坡度因子；
步骤3、计算相对平滑度指数作为权重因子，以反映山地森林地区土地覆盖和地表特征对泥沙连通性的影响；
步骤4、以谷歌卫星影像数据矢量化得到的流域内道路网络为基础，基于ArcGIS平台，运用要素转栅格工具和重分类工具区分流域内的道路和其它土地利用类型，并将道路栅格赋值为0，其它土地利用类型赋值为1，进而生成道路掩膜图层；
步骤5、以流域内DEM数据为基础，基于ArcGIS平台，获取水流流向与汇流累积量图层；
步骤6、通过空间分析中的水流长度工具计算沿流路径从每个像元到栅格边上的汇点或出水口的下坡距离，并以流域泥沙传输的最终流向作为水流流向的栅格数据、
作为权重栅格数据，从而得到下坡组分；
步骤7、借助汇流累积量工具分别以相对平滑度指数和坡度为权重栅格数据，计算不同权重栅格数据下的汇流累积量，根据不同权重栅格数据下的汇流累积量，通过栅格计算器工具计算上坡平均权重因子和上坡平均坡度，从而得到上坡组分；
步骤8、通过栅格计算器工具中的Log10函数，计算得到泥沙连通性指数。
2.根据权利要求1所述的一种考虑道路影响的流域泥沙连通性评价方法，其特征在于，所述步骤1中按照如下方式得到流域内道路网络与河流网络：
利用ArcGIS平台，基于2m分辨率的谷歌卫星影像数据，在Catalog下新建要素类型为折线的矢量图层，调用编辑器在编辑命令下描绘谷歌卫星影像数据中呈现的所有道路，描绘完成后点击编辑器下的“保存编辑内容”与“停止编辑”，即生成道路网络；重复上述步骤，生成河流网络。
3.根据权利要求1所述的一种考虑道路影响的流域泥沙连通性评价方法，其特征在于，所述步骤2还包括：
通过栅格计算器工具为坡度＝0的栅格字段赋值为0.005。
4.根据权利要求1所述的一种考虑道路影响的流域泥沙连通性评价方法，其特征在于，所述步骤3包括：
以流域内DEM为基础，通过Spatial Analyst Tools中的表面分析-坡度，计算流域内的坡度图层；随后，借助栅格计算器，通过地面粗糙度的计算公式，即
得到粗糙度图层；然后，根据粗糙度指数与曼宁系数的关系式，计算得到曼宁系数：
Y＝2.4234×e0.3005ln(n)
式中：Y为地表粗糙度，n为曼宁系数；
随后，按照下式得到相对平滑度指数：
式中：nmin为曼宁系数的最小值；RS为相对平滑度指数，是一个无量纲数字，从0.0～1.0不等。
5.根据权利要求1所述的一种考虑道路影响的流域泥沙连通性评价方法，其特征在于，所述步骤5包括：
步骤5.1、以流域DEM数据为基础，利用ArcGIS平台的填洼工具先对DEM进行填洼处理，随后运用空间分析工具中的水流流向工具提取DEM填洼后的水流方向，即FlowDir1；
步骤5.2、通过栅格计算器工具，将道路掩膜图层叠加到起始FlowDir1上，即道路掩膜*FlowDir1，从而得到路网影响下的水流流向，即FlowDir2；
步骤5.3、基于路网影响下的FlowDir2，运用空间分析中的汇流累积量工具计算栅格单元所流经的水量大小，即汇流累积量FlowACC1；
步骤5.4、根据FlowACC1，通过栅格计算器工具，结合矢量化河网确定阈值，从而提取流域河网掩膜图层；再次通过栅格计算器，将河网掩膜图层叠加在FlowDir2上，得到流域泥沙传输的最终流向，即FlowDir3；
步骤5.5、利用汇流累积量工具得到水流流向FowDir3的汇流累积量，即FlowACC2；利用栅格计算器，使FlowACC2每个栅格的汇流累积量值+1，得到最终的汇流累积量图层，即FlowACC3。
6.根据权利要求5所述的一种考虑道路影响的流域泥沙连通性评价方法，其特征在于，所述步骤6包括：
通过空间分析中的水流长度工具计算沿流路径从每个像元到栅格边上的汇点或出水口的下坡距离，并以FlowDir3作为水流流向的栅格数据、
作为权重栅格数据，得到下坡组分，具体公式如下：
式中：Ddn表示径流和泥沙沿路线到达指定流域出水口或沟道的可能性；di为研究单元到泥沙汇的距离；Wi是第i个单元的权重因子，此处用相对平滑度指数代替；Si是第i个单元的坡度因子；
借助栅格计算器工具为上述计算结果为0的栅格单元赋值。
7.根据权利要求1所述的一种考虑道路影响的流域泥沙连通性评价方法，其特征在于，所述步骤7包括：
借助汇流累积量工具分别以相对平滑度指数和坡度为权重栅格数据，计算不同权重栅格数据下的汇流累积量，即FlowRS和FlowS；随后根据FlowRS和FlowS，通过栅格计算器工具计算上坡平均权重因子和上坡平均坡度，并根据下式、借助栅格计算器，计算得到上坡组分：
式中：Dup表示上坡的连通分量，决定了上坡产生的径流和泥沙向下游流动的潜力；
是上坡贡献面积的平均权重因子，代表山区森林流域径流和泥沙在输移过程中可能遇到的阻抗；
为上坡贡献区域面积的平均坡度；
为上坡的贡献面积。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210211774.3A CN114580910B (zh) | 2022-03-04 | 2022-03-04 | 一种考虑道路影响的流域泥沙连通性评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202210211774.3A CN114580910B (zh) | 2022-03-04 | 2022-03-04 | 一种考虑道路影响的流域泥沙连通性评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN114580910A true CN114580910A (zh) | 2022-06-03 |
| CN114580910B CN114580910B (zh) | 2022-09-20 |
Family
ID=81778354
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202210211774.3A Active CN114580910B (zh) | 2022-03-04 | 2022-03-04 | 一种考虑道路影响的流域泥沙连通性评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114580910B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115081069A (zh) * | 2022-06-21 | 2022-09-20 | 河南大学 | 一种道路沉沙池规划方法及系统 |
| CN115860458A (zh) * | 2022-11-08 | 2023-03-28 | 广东省科学院广州地理研究所 | 滑坡横向泥沙连通性的评估方法、系统、设备和介质 |
| CN117172142A (zh) * | 2023-06-27 | 2023-12-05 | 长江水利委员会水文局 | 考虑地形影响的水文模型水流沿程再分配方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190354873A1 (en) * | 2018-02-16 | 2019-11-21 | Lucas Pescarmona | Analysis system and hydrology management for basin rivers |
| CN113627817A (zh) * | 2021-08-19 | 2021-11-09 | 东莞理工学院 | 一种基于城市内涝条件下道路连通性分析方法 |
-
2022
- 2022-03-04 CN CN202210211774.3A patent/CN114580910B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190354873A1 (en) * | 2018-02-16 | 2019-11-21 | Lucas Pescarmona | Analysis system and hydrology management for basin rivers |
| CN113627817A (zh) * | 2021-08-19 | 2021-11-09 | 东莞理工学院 | 一种基于城市内涝条件下道路连通性分析方法 |
Non-Patent Citations (2)
| Title |
|---|
| SHILIANG LIU ET AL.: "The Road Impacts on the River Landscape in a Typical Valley in Yunnnan Province", 《IEEE XPLORE》 * |
| 张润莲等: "规则网格数字高程模型中基于距离与坡度的路径规划算法", 《计算机应用》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115081069A (zh) * | 2022-06-21 | 2022-09-20 | 河南大学 | 一种道路沉沙池规划方法及系统 |
| CN115860458A (zh) * | 2022-11-08 | 2023-03-28 | 广东省科学院广州地理研究所 | 滑坡横向泥沙连通性的评估方法、系统、设备和介质 |
| CN117172142A (zh) * | 2023-06-27 | 2023-12-05 | 长江水利委员会水文局 | 考虑地形影响的水文模型水流沿程再分配方法 |
| CN117172142B (zh) * | 2023-06-27 | 2024-05-17 | 长江水利委员会水文局 | 考虑地形影响的水文模型水流沿程再分配方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN114580910B (zh) | 2022-09-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN114580910B (zh) | 一种考虑道路影响的流域泥沙连通性评价方法 | |
| Maksimović et al. | Overland flow and pathway analysis for modelling of urban pluvial flooding | |
| Kumar et al. | Surface runoff estimation of Sind river basin using integrated SCS-CN and GIS techniques | |
| Ahmad et al. | Estimation of a unique pair of Nash model parameters: an optimization approach | |
| CN109492259A (zh) | 一种城市水文模拟系统 | |
| Omran et al. | New ArcGIS tools developed for stream network extraction and basin delineations using Python and java script | |
| CN115391712A (zh) | 一种城市洪涝风险预测方法 | |
| Arthur et al. | Potential for small hydropower development in the Lower Pra River Basin, Ghana | |
| Kayembe et al. | Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape | |
| Amatya et al. | Application of LiDAR data for hydrologic assessments of low-gradient coastal watershed drainage characteristics | |
| Muluneh et al. | Morphometric analysis of Didessa river catchment in blue nile basin, western Ethiopia | |
| CN114528672A (zh) | 基于3s技术的城市水文站网布设方法及系统 | |
| Minh et al. | Application of GIS technology to establish a drainage density hierarchical map for flood hazard zoning in Lam river basin | |
| Athira et al. | Flood modelling and inundation mapping of Meenachil River using HEC-RAS and HEC-HMS software | |
| Jin et al. | Corrected GCM data through CMFD data to analysis future runoff changes in the source region of the Yangtze River, China | |
| CN114385959A (zh) | 一种近坝区子流域单元划分方法、装置及存储介质 | |
| Hamad | Surface runoff estimation of Wadi Ba Al-Arid watershed NE Libya using SCS-CN, GIS and RS data | |
| Šamanović et al. | Influence of pit removal algorithms on surface runoff simulation | |
| Mitterer | TATOO–Python Topographic Analysis Tool Library for semi-automated setup of high-resolution integrated hydrologic models | |
| Khemiri et al. | DEM-Based GIS Algorithms and 3D Geospatial Mapping for Creation of Hydrogeological Models Data in Foussana Basin (Central Tunisia) | |
| Shamsi | GIS and water resources modeling: state-of-the-art | |
| Upadhyay et al. | Delineation of flood risk zones in Jodhpur District, Rajasthan, India using remote sensing, GIS and multi-criteria decision making techniques | |
| Shahzad | Assessing, Mapping and Optimizing the Potential Dam Site Selection Using ArcGIS | |
| Muhaisen et al. | Prediction of Surface Runoff for Mosul Dam Reservoir from Different Regional Catchment Areas Using Arc SWAT Model | |
| Kumar et al. | Study of Rising Surface Water Levels on Land Submergence and Groundwater |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |