CN110232737A - 一种城市汇水区划分方法 - Google Patents

一种城市汇水区划分方法 Download PDF

Info

Publication number
CN110232737A
CN110232737A CN201910394827.8A CN201910394827A CN110232737A CN 110232737 A CN110232737 A CN 110232737A CN 201910394827 A CN201910394827 A CN 201910394827A CN 110232737 A CN110232737 A CN 110232737A
Authority
CN
China
Prior art keywords
remittance
charge
water
digital elevation
grid
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
Application number
CN201910394827.8A
Other languages
English (en)
Other versions
CN110232737B (zh
Inventor
胡潭高
沈黎达
李瑶
张路
李东臣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hangzhou Normal University
Original Assignee
Hangzhou Normal University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Hangzhou Normal University filed Critical Hangzhou Normal University
Priority to CN201910394827.8A priority Critical patent/CN110232737B/zh
Publication of CN110232737A publication Critical patent/CN110232737A/zh
Application granted granted Critical
Publication of CN110232737B publication Critical patent/CN110232737B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/05Geographic models
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/40Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Tourism & Hospitality (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Educational Administration (AREA)
  • Human Resources & Organizations (AREA)
  • Computer Graphics (AREA)
  • Remote Sensing (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Development Economics (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • General Business, Economics & Management (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Processing Or Creating Images (AREA)

Abstract

本发明为一种城市汇水区划分方法,包括获取精细化数字高程模型,对数字高程模型进行修正,对不同分辨率的数字高程模型进行汇水区划分和确定最佳分辨率。本发明的优点是:对数字高程模型进行精细化处理,充分考虑所有城市地形对汇水区划分的影响;基于城市排水点对数字高程模型进行修正,将排水点对于汇水区划分的影响转换为高程值,提高汇水区划分的准确性;通过确定最佳分辨率,充分考虑不同分辨率的DEM数据对于汇水区划分的影响,提高汇水区划分的合理性。

Description

一种城市汇水区划分方法
技术领域
本发明涉及城市水文领域,尤其涉及一种城市汇水区划分方法。
背景技术
汇水区划分(catchment delineation)是分布式水文模型中汇水区数据输入的基础和前提。分布式水文模型中的每个子汇水区是独立的水力学单元,在这些单元中,地形和排水系统因子使得地表径流直接汇入到一个排出点(O’Callaghan等,1984)。汇水区划分的尺度、准确性及流向真实性等会对模拟精度和效率产生重要影响(周思斯等,2014)。在研究城市暴雨积涝模型时,由于复杂的模型结构,涉及到更多的模型参数以及更大的数据量,汇水区划分的准确性会对模型结果产生更大的影响。同时,由于城市下垫面具有高度复杂、破碎、起伏微小却突兀等特征,诸如道路隔断、河道堤坝、地下空间、建筑结构等城市微地形会显著改变水流方向,对地表径流以及积水扩散过程产生明显影响,进而增加了汇水区划分的难度。因此,在建模前如何科学合理的划分子汇水区成为研究中一项十分重要工作。
目前,子汇水区划分方法大致归纳为以下类型:(1)人工勾画法。传统的汇水区划分一般以遥感影像图为背景,通过人工勾绘得到,但是对于大面积、多管线点的区域,这将是一项非常耗时且繁杂的工作,而且人工勾绘的随机性较大,地形因素考虑不周到,所以会对模型结果造成较大的影响。(2)排水管网为主的汇水区划分方法。该方法主要根据雨水井的空间位置采用泰森多边形的方法进行汇水区划分,使得每个雨水井在理论上处于汇水区的中心。该方法虽然考虑了排水管网的影响,但是未充分考虑地形对径流的影响,难以准确模拟城市地表径流的真实情况。
发明内容
本发明主要解决了现有的汇水区划分方法忽略了城市管网、道路、建筑物等复杂下垫面对于径流方向的影响,缺乏将地形和城市排水管网结合考虑,提供了一种综合考虑城市复杂下垫面及地形和排水网关的城市汇水区划分方法。
本发明解决其技术问题所采用的技术方案是,一种城市汇水区划分方法,包括以下步骤:
S1:获取精细化数字高程模型;
S2:对数字高程模型进行修正;
S3:对不同分辨率的数字高程模型进行汇水区划分;
S4:确定最佳分辨率。
对城市不同的地形进行补偿后获取精细化数字高程模型,提高水流方向判断的准确性,结合城市地下网管数据对精细化数字高程模型进行进一步修正,提高汇水区划分的准确性。
作为上述方案的一种优选方案,所述步骤S1,包括以下步骤:
S11:获取城市的数字高程模型;
S12:利用获取的城市土地利用类型数据对研究区数字高程模型中对应的栅格单元进行补偿。城市土地利用类型包括道路、建筑、绿地、水渠等,根据不同的城市土地利用类型对相应的栅格单元进行高程值增加或减少。
作为上述方案的一种优选方案,所述步骤S2,包括以下步骤:
S21:确定研究区的排水点;
S22:以每个排水点为中心建立缓冲区;
S23:对高程减少计算公式对缓冲区的高程作降低处理,高程减少计算公式如下,
式中h为减少的高程值,d为排水点雨水井的深度,单位为m;S截面为排水点的横截面积,单位为m2;n为缓冲区分级数量;S范围为排水点影响范围的面积,单位为m2。将地下管网对于汇水区划分的影响转化为单元栅格的高程值,提高最终汇水区划分的准确度。
作为上述方案的一种优选方案,所述步骤S3,包括以下步骤:
S31:对精细化数字高程模型进行重采样,得到不同分辨率的数字高程模型数据;
S32:进行汇水区划分。
作为上述方案的一种优选方案,所述步骤S32中汇水区划分,包括以下步骤:
S321:对数字高程模型进行填洼处理,获得无洼地数字高程模型;
S322:提取每个栅格单元的水流方向,获取区域地形的水流方向数据;
S323:根据区域地形的水流方向数据计算每一栅格流过的水量数值,得到该区域的汇流累积量;
S324:确定划分阈值,以汇流累积量大于划分阈值为依据的提取出河网水系;
S325:确定区域出水口,根据水流方向找到该出水口上游的栅格,确定汇水区。
作为上述方案的一种优选方案,步骤S322中水流方向通过以下步骤获取:
S01:计算一个栅格单元与该栅格单元相邻栅格单元的距离权落差,距离权落差为两个栅格单元中心点的高程差与距离的比值;
S02:筛选出距离权落差最大的栅格为该栅格流出栅格;
S03:确定水流方向为流出栅格对应该栅格的方向。
作为上述方案的一种优选方案,所述划分阈值利用对应地区的现有地形图确定。
本发明的优点是:对数字高程模型进行精细化处理,充分考虑所有城市地形对汇水区划分的影响;基于城市排水点对数字高程模型进行修正,将排水点对于汇水区划分的影响转换为高程值,提高汇水区划分的准确性;通过确定最佳分辨率,充分考虑不同分辨率的DEM数据对于汇水区划分的影响,提高汇水区划分的合理性。
附图说明
图1为本发明的一种流程示意图。
图2为本发明中获取精细化数字高程模型的一种流程示意图。
图3为本发明中对数字高程模型进行修正的一种流程示意图。
图4为本发明中对不同分辨率的数字高程模型进行汇水区划分的一种流程示意图。
图5为本发明中汇水区划分的一种流程示意图。
图6为本发明中水流方向获取的一种流程示意图。
具体实施方式
下面通过实施例,并结合附图,对本发明的技术方案作进一步的说明。
实施例1:
本实施例一种城市汇水区划分方法,如图1所示,包括以下步骤:
S1:获取精细化数字高程模型,如图2所示,包括以下步骤:
S11:利用地形图、高分辨率遥感影像、城市行政边界数据、城市土地利用数据等获取城市的数字高程模型;
S12:利用获取的城市土地利用类型数据对研究区数字高程模型中对应的栅格单元进行补偿,得到精细化数字高程模型,城市土地利用类型包括道路、建筑物、绿地和水渠,道路的凸起容易阻止水系的生成,两侧低洼地势会将降雨引入周边河流或者排水口,因此,对道路图层所占栅格处的高程值增加相应的值;密集分布的建筑物常常阻隔自然径流,自然径流常常沿建筑物边界汇入河流。因此,对建筑物所占栅格高程值增加相应的值;绿地作使得降雨到达地表后部分水量耗损于植被截留,并且粗糙系数较大,对地表径流有较大的阻碍作用,会影响地表径流的方向。对绿地所占栅格高程值减少相应的值;水渠的地势较低,对于地表径流起到分流作用。因此,对沟渠图层的高程值减去一个值;不同地形增加或减少的值具体由现场调查获取。
S2:对数字高程模型进行修正,如图3所示,包括以下步骤:
S21:以研究区的地下管网数据为基础,根据主干道路交叉点和建筑物的分布确定研究区的排水点;
S22:以每一个排水点为中心建立三级缓冲区;
S23:对高程减少计算公式对缓冲区的高程作降低处理,高程减少计算公式如下,
式中h为减少的高程值,d为排水点雨水井的深度,单位为m;S截面为排水点的横截面积,单位为m2;n为缓冲区分级数量;S范围为排水点影响范围的面积,单位为m2
S3:对不同分辨率的数字高程模型进行汇水区划分,如图4所示,包括以下步骤:
S31:对精细化数字高程模型进行重采样,得到分辨率为2m、4m、6m、8m的数字高程模型数据;
S32:进行汇水区划分,如图5所示,包括以下步骤:
S321:结合实际地形对数字高程模型在采集时因误差产生的洼地进行填洼处理,获得无洼地数字高程模型,;
S322:提取每个栅格单元的水流方向,获取区域地形的水流方向数据,如图6所示,包括以下步骤:
S01:在3×3的栅格区域上,计算中心栅格与各相邻栅格间的距离权落差,距离权落差为两个栅格单元中心点的高程差与距离的比值;
S02:筛选出距离权落差最大的栅格为该栅格流出栅格;
S03:在3×3的栅格区域上,中心栅格有8个有效输出方向,根据不同方向对8个相连栅格进行编码,从中心栅格正右方的栅格开始,以顺时针方向依次编码为1、2、4、8、16、32、64、128,分别代表中心栅格的东、东南、南、西南、西、西北、北、东北8个方向,每一个中心栅格的水流都向都由这8个值中的某一个值来确定,流出栅格对应的编码即为该中心栅格的水流方向;
S323:数字高程模型每点处均有一个单位的水量,按照水流从高处往低处流的自然规律,根据区域地形的水流方向数据计算每一栅格处所流过的水量数值,便得到了该区域的汇流累积量;
S324:确定划分阈值,以汇流累积量大于划分阈值为依据的提取出河网水系,不同区域对应的阈值不同,所以应根据研究区实际的地物状况,利用现有地形图等其他资料辅助检验的方法来确定合适的划分阈值,从而提取出符合实际的河网水系,以海宁市为例,最佳划分阈值为450;
S325:先确定出水口,即该集水区的最低点,然后结合水流方向,分析搜索该出水口上游所有流过该出水口的栅格,搜索至流域的边界为止。
S4:确定最佳分辨率,经过上述步骤后会得到4种分辨率的汇水区划分结果,对不同分辨率下的海宁市汇水区划分结果进行局部放大进行比较,2m和4m分辨率的DEM数据所划分的汇水区较为细碎,遥感图像上的建筑物区块被分为独立的子汇水区数量过多,不符合实际汇流情况,且过高分辨率DEM划分产生了一些错误的平形状汇水区。8m分辨率的DEM数据所划分的汇水区数量过少,6m分辨率DEM得到的结果沿主干道将建筑物区块划分为四个主要的子汇水区,且将四个相对独立的社区区块划分开来,汇水区边界也较为贴合建筑物和道路,符合实际的汇流情况。因此确定最佳分辨率为6m。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。

Claims (7)

1.一种城市汇水区划分方法,其特征是:包括以下步骤:
S1:获取精细化数字高程模型;
S2:对数字高程模型进行修正;
S3:对不同分辨率的数字高程模型进行汇水区划分;
S4:确定最佳分辨率。
2.根据权利要求1所述的一种城市汇水区划分方法,其特征是:所述步骤S1,包括以下步骤:
S11:获取城市的数字高程模型;
S12:利用获取的城市土地利用类型数据对研究区数字高程模型中对应的栅格单元进行补偿。
3.根据权利要求1所述的一种城市汇水区划分方法,其特征是:所述步骤S2,包括以下步骤:
S21:确定研究区的排水点;
S22:以每个排水点为中心建立缓冲区;
S23:对高程减少计算公式对缓冲区的高程作降低处理,高程减少计算公式如下,
式中h为减少的高程值,d为排水点雨水井的深度,单位为m;S截面为排水点的横截面积,单位为m2;n为缓冲区分级数量;S范围为排水点影响范围的面积,单位为m2
4.根据权利要求1所述的一种城市汇水区划分方法,其特征是:所述步骤S3,包括以下步骤:
S31:对精细化数字高程模型进行重采样,得到不同分辨率的数字高程模型数据;
S32:进行汇水区划分。
5.根据权利要求4所述的一种城市汇水区划分方法,其特征是:所述步骤S32中汇水区划分,包括以下步骤:
S321:对数字高程模型进行填洼处理,获得无洼地数字高程模型;
S322:提取每个栅格单元的水流方向,获取区域地形的水流方向数据;
S323:根据区域地形的水流方向数据计算每一栅格流过的水量数值,得到该区域的汇流累积量;
S324:确定划分阈值,以汇流累积量大于划分阈值为依据的提取出河网水系;
S325:确定区域出水口,根据水流方向找到该出水口上游的栅格,确定汇水区。
6.根据权利要求5所述的一种城市汇水区划分方法,其特征是:所述步骤S322中水流方向通过以下步骤获取:
S01:计算一个栅格单元与该栅格单元相邻栅格单元的距离权落差,距离权落差为两个栅格单元中心点的高程差与距离的比值;
S02:筛选出距离权落差最大的栅格为该栅格流出栅格;
S03:确定水流方向为流出栅格对应该栅格的方向。
7.根据权利要求5所述的一种城市汇水区划分方法,其特征是:所述划分阈值利用对应地区的现有地形图确定。
CN201910394827.8A 2019-05-13 2019-05-13 一种城市汇水区划分方法 Active CN110232737B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910394827.8A CN110232737B (zh) 2019-05-13 2019-05-13 一种城市汇水区划分方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910394827.8A CN110232737B (zh) 2019-05-13 2019-05-13 一种城市汇水区划分方法

Publications (2)

Publication Number Publication Date
CN110232737A true CN110232737A (zh) 2019-09-13
CN110232737B CN110232737B (zh) 2023-04-18

Family

ID=67860544

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910394827.8A Active CN110232737B (zh) 2019-05-13 2019-05-13 一种城市汇水区划分方法

Country Status (1)

Country Link
CN (1) CN110232737B (zh)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111538798A (zh) * 2020-04-09 2020-08-14 武汉大学 一种顾及dsm和dlg的城市汇水区精细化提取方法
CN111680886A (zh) * 2020-05-18 2020-09-18 华南师范大学 一种内涝风险预测方法及其系统
CN112242003A (zh) * 2020-10-19 2021-01-19 中国测绘科学研究院 一种顾及地类与流向的城市子汇水区划分方法
CN114547531A (zh) * 2022-02-25 2022-05-27 中国水利水电科学研究院 一种城市不透水面有效性量化方法
CN116011154A (zh) * 2022-11-30 2023-04-25 天津大学 一种基于道路高程的城市汇水区划分方法
CN116109143A (zh) * 2023-04-07 2023-05-12 南京信息工程大学 面向暴雨洪涝灾害风险分析的地形起伏影响指数确定方法、系统、装置及存储介质

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110219009A1 (en) * 2010-03-03 2011-09-08 Cellguide Ltd. Method of generating an elevation database
WO2013121340A1 (en) * 2012-02-13 2013-08-22 Stellenbosch University Digital elevation model
CN107563019A (zh) * 2017-08-17 2018-01-09 河北工程大学 一种针对复杂下垫面城区分布式水文模型数字化方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20110219009A1 (en) * 2010-03-03 2011-09-08 Cellguide Ltd. Method of generating an elevation database
WO2013121340A1 (en) * 2012-02-13 2013-08-22 Stellenbosch University Digital elevation model
CN107563019A (zh) * 2017-08-17 2018-01-09 河北工程大学 一种针对复杂下垫面城区分布式水文模型数字化方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
左俊杰等: "平原河网地区汇水区的划分方法――以上海市为例", 《水科学进展》 *
左文君等: "基于DEM的城市数字河网提取", 《水电能源科学》 *
薛丰昌等: "面向城市平原地区暴雨积涝汇水区分级划分的方法研究", 《地球信息科学学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111538798A (zh) * 2020-04-09 2020-08-14 武汉大学 一种顾及dsm和dlg的城市汇水区精细化提取方法
CN111538798B (zh) * 2020-04-09 2023-09-19 武汉大学 一种顾及dsm和dlg的城市汇水区精细化提取方法
CN111680886A (zh) * 2020-05-18 2020-09-18 华南师范大学 一种内涝风险预测方法及其系统
CN112242003A (zh) * 2020-10-19 2021-01-19 中国测绘科学研究院 一种顾及地类与流向的城市子汇水区划分方法
CN114547531A (zh) * 2022-02-25 2022-05-27 中国水利水电科学研究院 一种城市不透水面有效性量化方法
CN114547531B (zh) * 2022-02-25 2022-09-06 中国水利水电科学研究院 一种城市不透水面有效性量化方法
CN116011154A (zh) * 2022-11-30 2023-04-25 天津大学 一种基于道路高程的城市汇水区划分方法
CN116109143A (zh) * 2023-04-07 2023-05-12 南京信息工程大学 面向暴雨洪涝灾害风险分析的地形起伏影响指数确定方法、系统、装置及存储介质

Also Published As

Publication number Publication date
CN110232737B (zh) 2023-04-18

Similar Documents

Publication Publication Date Title
CN110232737A (zh) 一种城市汇水区划分方法
CN109543275B (zh) 一种城区地表径流二维数值模拟方法
CN108399309B (zh) 一种大尺度复杂地形区分布式水文模型的子流域划分方法
CN104978451B (zh) 一种基于ArcGIS的复杂平原河网地区水文模型处理方法
CN108388714A (zh) 流域水系和城市管网耦合的平原河网城市洪水模拟方法
CN108643116B (zh) 一种山区性水库库区河道水面宽度的估算方法
CN107239657B (zh) 一种面向对象的水动力学建模要素管理方法
CN111339711B (zh) 一种小流域设计洪水推求方法
CN110298076A (zh) 一种基于gis和swmm的城市内涝智能建模及分析方法
CN110543984B (zh) 一种基于水安全的城市水面率规划方法
CN106920202A (zh) 一种平原城市河网畅流活水方法
CN107133427A (zh) 一种基于2dgis平台的洪水分析模型的构建方法
CN109657841A (zh) 一种城市暴雨内涝积水深度提取方法
CN111475950B (zh) 一种下凹式立交桥下雨洪的模拟方法
CN108763615A (zh) 基于管网和道路双系统swmm对城市涝积水深模拟方法
CN114647881B (zh) 顾及建筑物微观水文过程的城市内涝建模方法
CN111062128B (zh) 一种流域汇流模拟估测方法
CN110532641B (zh) 一种地表网格分层建模方法及系统
CN110188323B (zh) 一种适宜于城市水生态保护的水面率值域测算方法
CN115358650B (zh) 一种洪涝灾害应急避险与物资实时调配方法
CN112711825B (zh) 一种地表与管网分布式的直接双向耦合方法
Formánek et al. Two-dimensional model of Ciliwung river flood in DKI Jakarta for development of the regional flood index map
CN111046551A (zh) 一种城市群排水过程模拟方法
CN108319758A (zh) 一种基于水力模型的城市下穿地道排水系统优化设计方法
CN112242003A (zh) 一种顾及地类与流向的城市子汇水区划分方法

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