水文分析之水库库容计算

作为成熟的GIS软件,ArcGIS软件在地表水分析方面的研究成果还是相当不错的,目前已经非常成熟;本章节主要介绍ArcGIS水库库容计算,力求让大家了解其运作原理。了解原理之后才能制定其分析流程。


水库库容计算的知识准备


在开始前需要有一定的ArcGIS的水文分析基础,如果不了解原理,换了一个数据,换了一个地点,可能就无法计算出结果。你需要掌握:
1、必须了解相关的水文专业名词,如河流、出水口、流域、分水岭、河流链接等。水文专业人员基本上没问题,但如果其他非水利从业者,估计需要花点时间进行词汇的对接。
2、必须了解一些特定的GIS名词。如DEM数据、栅格数据、矢量数据。GIS常用的数据格式。
3、了解ArcGIS 水文分析工具的原理、流程及每个工具的细节参数对水文分析的影响。
4、ArcGIS 桌面软件是一个非常庞大的体系,除了掌握水文分析的几个工具之外,必须配合其他工具一起使用,例如栅格转面、重分类、条件函数、要素折点转点等一些周边的工具。
5、最好能够掌握modelbuilder建模方法,对于复杂的流程可以制作分析模型来操作。


水库库容计算的基本流程


在水利行业中,计算库容自然有自己的一套公式和方法。也不是说这种方法不准确和不科学,主要是能够有比较直观的、快捷的方式,而且是可调节参数的方式,不妨可以通过新的方法去替换掉原有的方法。
如果基于传统的计算方式,计算水库的库容只是得出一串数字,那在水利的综合应用上没什么大的意义,充其量也就是水利工程上得知具体的数值而后确定其方案的合理性。但跟GIS结合在一起就完全不一样了。
GIS具有非常高的运算能力,有非常良好的地图可视化功能,可以通过调节不同的参数计算出不同水位淹没的区域,从而计算出库容。 淹没的空间范围可以看作水库移民的依据,与土地利用图层叠加分析,还可以计算出哪些地方被淹没后需要作出多少的补偿。如果库容不足,是否需要开挖,挖方量是多少,结合工程施用需要多少泥头车进行拉运泥土。移民、征地、青苗补偿、开挖、运泥等等一系列的问题都跟钱(财政)有着密切的联系。所以,如果将水库库容计算只看作一次性的数值计算,那就相当浪费了。


水库库容的计算方法


通过水库的坝口(也就是说要筑坝的地方,可以用点来表示)作为出水口,计算出该区域的集水区; 利用集水区裁剪出DEM数据; 设定该区域一个最高的蓄水高程(指的是水库蓄水最高时水平面的海拔高度)与集水区的DEM数据进行裁剪。
裁剪后的结果可以利用ArcGIS的表面体积工具计算出其实际库容。


分析流程


第一步:地形预处理,生成无凹陷点的DEM数据。


第二步:利用坝口作为出水口,捕捉到最近的倾泻点,并计算出该倾泻点的流域。其实就是一个流域提取的过程。


将一些必要的信息加入到地图中,例如坝口位置,图中黄色十字位置。坝口位置一般为选定的水坝位置。如下图

水文分析之水库库容计算的图1

通过无凹陷点的DEM数据,进行流向和流量分析。然后利用捕捉倾泻点工具捕捉到坝口的倾泻点栅格。如下图所示

水文分析之水库库容计算的图2

虽然坝口位置即便可以看成是出水口,但也不能直接使用。因为它与原DEM不是同一份数据,位置上跟DEM对应的出水口还是有差别的。所以,必须要进行捕捉,以确保捕捉到的倾泻点是该点流量最大的栅格,否则就会出现那种很小的面。
利用【分水岭】工具计算出坝口的集水区。 这里要求输入的倾斜点就是捕捉到的倾斜点。如下图所示

水文分析之水库库容计算的图3

其结果是一个栅格面数据,如下图

水文分析之水库库容计算的图4

为了方便后续的分析和查看,可以利用【栅格转面】工具将这个栅格的集水区转换为矢量面数据。通过将矢量面设置只显示边界,去除填充部分,就可以非常清晰的看到这个集水区的边界了。下图是该集水区的边界范围

水文分析之水库库容计算的图5

第三步:利用集水区边界裁剪原始的DEM 。要注意的是,这里的DEM是原始的DEM,而不是填洼过的无凹陷点DEM。切记,因为后续要基于这个准确的DEM去计算体积。
利用【数据管理】-【栅格】-【栅格处理】-【裁剪】工具,对该原始DEM进行栅格裁剪。 裁剪的时候一定要注意,必须要勾选上【使用输入要素裁剪几何】这个参数,不然裁剪出来的是集水面的矩形区域。

水文分析之水库库容计算的图6

其实关于水文分析部分到这一步就已经完全完成了。接下来就是利用GIS工具计算一定该区域“蓄水”的计算了,这就是传统GIS部分的计算。
第四步:根据水位的海拔高度计算出水库蓄水时淹没的区域。
这是一个关键的步骤,因为蓄水高度一旦定下来,那水面与下垫面之间的空间就是整个水库的库容了。怎么计算蓄水的边界是个问题。
ArcGIS提供了栅格计算器工具,可以使用setnull工具来计算;或者直接调用【spatial analyst】–【条件分析】–【设为空函数】工具。
笔者强烈推荐使用第二种方法,因为在modelbuilder里面这个工具比较直观。
这个工具的作用是将识别到的像元值设置为NoData. 假如将集水区栅格通过该工具查找象元值大于2000米,将会将大于2000米的所有像元都设置为没有数据,这样的出来的就是一个只有最高高程为2000米的区域。这就好比水面像一把刀子一样将高于2000米的区域削去。
其设置如下图

水文分析之水库库容计算的图7


最终可以得到一个最大水位时的淹没区,该淹没区也是个栅格数据,如下图所示


水文分析之水库库容计算的图8

第五步:计算最终的库容。 使用【3D 分析】–【功能性表面】–【表面体积】工具,计算这个淹没区下垫面的体积。
该工具如下图所示,在【参考平面】参数选择Below ,这表示计算下垫面的体积。

水文分析之水库库容计算的图9

最终计算出来的是一个文本文件,里面包含了一些结果信息。如下图

水文分析之水库库容计算的图10

可以看到该区域的蓄水量就是16370702100 立方米。
假如这个数值还没达到预期的库容量,可以通过调整蓄水面的高度进行计算。当然,如果高度已经是工程评估过的最高水位,仍达不到理想的库容,要么就是放弃该区域,要么就是通过人工挖方来扩容。将理想的库容量减去计算出来的库容量,就是需要挖方的实际数值。
为了方便,笔者特意利用modelbuilder写了一个处理流程,这样比较方便进行数据上的调整。

水文分析之水库库容计算的图11

水文分析之水库库容计算的图12

在前后两处蓄水高度的调整可以直接输入数值,然后对输出结果进行比较。如下图

水文分析之水库库容计算的图13


这样就方便多了。


属于GIS的外延部分


水库库容的计算是给水利工程使用的。 对于这个工程对周边的影响,其计算就是GIS最拿手的地方。
例如通过计算出蓄矢量面,与居民点进行叠加,可以找出需要移民的居民地。如下图高亮选中的居民点属于必须移民的居民点。

水文分析之水库库容计算的图14

而落在黄色区域外红色区域内的居民点属于存在淹没风险的区域。而这两个区域划定之后,基本上可以确定哪些农田是被淹没的,哪些农田是没有淹没,还可以正常耕作,这样就为移民补偿作精准的依据。


总结


利用GIS做水库库容计算,既需要充分利用水文分析的一些知识,也需要数量掌握GIS有效的工具,这样才能做到科学和准确。 同时GIS作为一种可视化的手段,比实际数学公式计算出来的数值结果更加直观,更具说服力,能够在其他环节中为水库的选择作出科学的断定。

文章来源:水利水电资料库

水文分析水库

水文分析之水库库容计算的评论0条

    暂无评论

    水文分析之水库库容计算的相关案例教程

    影响地质灾害形成的自然因素众多,历史地震地质灾害发生的数量、分布范围、活动规模都直接反映了地层岩性、地形地貌、现存新老滑坡以及有关地震动力环境对地震诱发地质灾害的控制作用;此外,土地利用、地下水、植物条件等因素也对震后地质灾害形成起到一定程度的影响。本文采用统计学方法,对研究区山地灾害点与各因子的每个属性进行相对频率组合的定量计算方法,综合天水市秦州区震后地质灾害发育情况,本次危险性区划分析中选用
    影响地质灾害形成的自然因素众多,历史地震地质灾害发生的数量、分布范围、活动规模都直接反映了地层岩性、地形地貌、现存新老滑坡以及有关地震动力环境对地震诱发地质灾害的控制作用;此外,土地利用、地下水、植物条件等因素也对震后地质灾害形成起到一定程度的影响。本文采用统计学方法,对研究区山地灾害点与各因子的每个属性进行相对频率组合的定量计算方法,综合天水市秦州区震后地质灾害发育情况,本次危险性区划分析中选用
    去年,千寻位置发布了入门级航测新势力“千巡翼X1无人机”,身材小巧的它重量仅1237g,轴距397mm,续航时间可达36分钟,标配1英寸CMOS 2000万像素相机,支持实时厘米级定位,六向720°全向避障,让免像控高精度测绘作业得心应手。 自发布以来,业内已使用千巡翼X1进行了多种类型的航测作业。今天让我们一起看看,在古建筑实景三维建模、正射影像、高精度带状正射影
    1、影像格式的转换 例如把jpg格式转换为tiff格式。可以在arctoolbox中的conversiontools-->to Raster-->Raster to Other Format multiple)。 2、矢量化准备 矢量化前建立一个.mxd工程文件,该文件为索引文件,记录的是所有图层的相对路径和图层渲染标注等信息。然后通过arcaltalog建立点线面图层,添加到工程中即可进行矢量化
    还是让人头疼的投影坐标转换。 推荐下面这篇,该文章分了4个点讲解,希望能给你解惑。 本文要点: 1、地理坐标系是什么?常见的地理坐标系有哪些? 2、投影坐标系是什么?常见的投影坐标系有哪些? 3、关于坐标转换。 4、如何在ArcGis里玩转坐标系?中国有哪些另类的坐标系? 在进行空间数据相关的分析时,坐标系是一个绕不过的坎。 建立工程时需要选择坐标系,成果数据有时还需要进行坐标转换、Cad数据导入
    影响力
    粉丝
    内容
    获赞
    收藏
      项目客服
      培训客服
      0 0