CN114547531B - 一种城市不透水面有效性量化方法 - Google Patents
一种城市不透水面有效性量化方法 Download PDFInfo
- Publication number
- CN114547531B CN114547531B CN202210182297.2A CN202210182297A CN114547531B CN 114547531 B CN114547531 B CN 114547531B CN 202210182297 A CN202210182297 A CN 202210182297A CN 114547531 B CN114547531 B CN 114547531B
- Authority
- CN
- China
- Prior art keywords
- effectiveness
- urban
- grids
- impervious surface
- runoff
- 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.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- Y—GENERAL 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
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A30/00—Adapting or protecting infrastructure or their operation
- Y02A30/60—Planning or developing urban green infrastructure
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Processing Of Solid Wastes (AREA)
Abstract
本发明提供了一种城市不透水面有效性量化方法,属于城市不透水面量化技术领域。该方法包括:基于城市不透水面的空间分布,分别获取高程数据和土地利用数据;根据所述高程数据和土地利用数据,计算得到水流方向;根据所述水流方向和流径下垫面类型,计算得到研究区域出口最终径流量;根据所述研究区域出口最终径流量,计算得到城市不透水面有效性。利用该方法进行城市下垫面空间特征对地表产汇流的影响研究,能更进一步分析影响程度,推动相关研究的进展,实现城市不透水面有效性的高效且精准的计算,将有助于提高常规城市不透水面的计算的精度。
Description
技术领域
本发明属于城市不透水面量化技术领域,尤其涉及一种城市不透水面有效性量化方法。
背景技术
作为城市下垫面的最主要组成部分,城市不透水面是引发城市热岛效应的主要因素,城市不透水面通常定义为阻止地表水入渗裸土的人造地表,主要由建设用地等城市用地类型构成,由于城市建设用地包括房屋、道路和广场等多种形式,且随着城市的发展不断扩张,城市不透水面逐步演绎成多种空间分布形式。研究不透水面空间分布对产汇流影响机理的重要基础是对其进行量化表征,所以一种新型的城市不透水面有效性量化方法的提出十分有必要。在过去城市径流量计算,通常是采用降雨量与径流系数的乘积形式。在城市中,不透水面是径流系数的主要影响因素,所以研究不透水面的有效性至关重要。由于过去多使用不透水面的比例作为不透水率,但是不透水面在空间上的分布各不相同,可能会直接导致不透水面的有效性有所不同。因此,不透水面可分为“形式不透水面”和“有效不透水面”。通过本发明可以解决如何量化不透水面的有效性问题,对以后解决相关研究问题和运用到模型计算中都具有重要意义。
发明内容
针对现有技术中的上述不足,本发明提供的一种城市不透水面有效性量化方法,能够实现城市不透水面有效性的高效且精准的计算,将有助于提高常规城市不透水面的计算的精度。
为了达到以上目的,本发明采用的技术方案为:
本方案提供一种城市不透水面有效性量化方法,包括以下步骤:
S1、基于城市不透水面的空间分布,分别获取高程数据和土地利用数据;
S2、根据所述高程数据和土地利用数据,计算得到水流方向;
S3、根据所述水流方向和流径下垫面类型,计算得到研究区域出口最终径流量;
S4、根据所述研究区域出口最终径流量,计算得到城市不透水面有效性。
本发明的有益效果是:本发明首先基于城市不透水面的空间分布,获取高程数据和土地利用数据,并了解每个部分的降雨量,计算水流方向,依据最陡坡度确定唯一流向,计算得到不同部分的水流流向,再根据流经下垫面类型确定研究区域出口最终径流量,最后计算不透水面有效性。利用该方法进行城市下垫面空间特征对地表产汇流的影响研究,能更进一步分析影响程度,推动相关研究的进展,实现城市不透水面有效性的高效且精准的计算,将有助于提高常规城市不透水面的计算的精度。
再进一步地,所述步骤S2包括以下步骤:
S201、根据所述高程数据和土地利用数据,针对每个栅格对周围八个方向均有高程差,且水流方向最陡时,计算中心栅格与相邻八个方向的栅格间的坡度;
S202、根据所述中心栅格与相邻八个方向的栅格间的坡度,确定每个栅格的最陡下降方向;
S203、判断所述每个栅格的最陡下降方向是否均相同,若是,则进入步骤S204,否则,进入步骤S205;
S204、扩大相邻栅格范围,并返回步骤S201;
S205、得到最大下降方向,并根据所述最大下降方向得到水流方向。
上述进一步方案的有益效果是:本发明的创新在于根据地理信息技术提出了一种新的计算方法,可以更加精确的分析城市不透水面上的产流过程,通过此算法可以准确的计算出水流过程,以便于简化出口最终径流量的计算过程。
再进一步地,所述步骤S201中栅格间的坡度的表达式如下:
其中,J表示栅格间的坡度,Δh表示栅格间的水平距离差,Δs表示栅格间的高程差。
上述进一步方案的有益效果是:本发明通过计算栅格间的坡度,为后续判断最大下降方向提供了前提条件。
再进一步地,所述步骤S3中总的径流量的表达式如下:
其中,F表示研究区域出口最终径流量,Pi表示研究区域降雨量,a表示透水面径流系数,m表示水流经过透水面的个数,i表示研究区域的编号,i=1,2,3,...,n,n表示研究区域的总编号。
上述进一步方案的有益效果是:本发明通过上述方法可以区别原有的径流量计算,并且准确的计算出区域出口最终径流量,以便于进行城市不透水面有效性的分析。
再进一步地,所述步骤S4中城市不透水面有效性的表达式如下:
其中,S表示城市不透水面有效性,F表示研究区域出口最终径流量,F0表示基于不透水率得到的出流量。
上述进一步方案的有益效果是:本发明提出了一种原理明确、操作简便一种城市不透水面有效性量化计算方法,利用该方法进行城市下垫面空间特征对地表产汇流的影响研究,能更进一步分析影响程度,推动相关研究的进展。
附图说明
图1为本发明的方法流程图。
图2为本实施例中城市下垫面复杂空间特征对地表产汇流过程示意图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
实施例
如图1所示,本发明提供了一种城市不透水面有效性量化方法,其实现方法如下:
S1、基于城市不透水面的空间分布,分别获取高程数据和土地利用数据;
本实施例中,高程数据:获取研究区DEM数据,数字高程模型(Digital ElevationModel),简称DEM,是通过有限的地形高程数据实现对地面地形的数字化模拟(即地形表面形态的数字化表达),它是用一组有序数值阵列形式表示地面高程的一种实体地面模型,是数字地形模型的一个分支,其它各种地形特征值均可由此派生。土地利用数据:根据土地利用数据获取不同高程网格的土地利用类型,根据土地利用类型获取相应的下渗系数。
S2、根据所述高程数据和土地利用数据,计算得到水流方向,其实现方法如下:
S201、根据所述高程数据和土地利用数据,针对每个栅格对周围八个方向均有高程差,且水流方向最陡时,计算中心栅格与相邻八个方向的栅格间的坡度;
S202、根据所述中心栅格与相邻八个方向的栅格间的坡度,确定每个栅格的最陡下降方向;
S203、判断所述每个栅格的最陡下降方向是否均相同,若是,则进入步骤S204,否则,进入步骤S205;
S204、扩大相邻栅格范围,并返回步骤S201;
S205、得到最大下降方向,并根据所述最大下降方向得到水流方向。
本实施例中,本发明提供了一种对于下垫面水流流向的新型算法,根据以获取的DEM数据,分别对每个栅格计算流向,假设每个栅格对周围8个方向均有高程差,且假设水流向最陡,计算中心栅格与相邻8个方向的栅格间的坡度,水流向为来自每个区域栅格的最陡下降方向,最陡下降方向计算如下:
首先,用下式计算栅格之间的坡度:
其中,J表示栅格间的坡度,△h表示栅格间的水平距离差,Δs表示栅格间的高程差。
S3、根据所述水流方向和流径下垫面类型,计算得到研究区域出口最终径流量;
本实施例中,采用径流系数法,计算研究区域出口最终径流量:
其中,F表示研究区域出口最终径流量,Pi表示研究区域降雨量,a表示透水面径流系数,m表示水流经过透水面的个数,i表示研究区域的编号,i=1,2,3,...,n,n表示研究区域的总编号。
S4、根据所述研究区域出口最终径流量,计算得到城市不透水面有效性:
其中,S表示城市不透水面有效性,F表示研究区域出口最终径流量,F0表示基于不透水率得到的出流量。
本实施例中,在过去计算城市径流量,通常采用降雨量乘以径流系数,也就是说仅仅根据土地利用数据得到的不透水率来计算径流,而不考虑空间分布的影响,因此,本发明通过上述方法计算得到的径流量,在上述公式中定义为F0和F。
本实施例中,在之前的研究中,一部分不透水面不直接与排水***相连,其产流在汇流到排水口之前,会流经其他透水面,称为非有效透水面(在本发明中为形式不透水面),而有效不透水面直接与排水***相连,为了精确的分析城市地表产汇流本发明提出了一种新的计算不透水面有效性的方法,由于两种径流量的计算考虑因素不同,其得到的值可能会大于1,也充分说明了研究的重要性。
本发明首先基于城市不透水面的空间分布,对研究区域按高程、下垫面类型划分为多个部分,并了解每个部分的降雨量,计算水流方向,依据最陡坡度确定唯一流向,计算得到不同部分的水流流向,再根据流经下垫面类型确定总的径流量,最后计算不透水面有效性,利用该方法进行城市下垫面空间特征对地表产汇流的影响研究,能更进一步分析影响程度,推动相关研究的进展,实现城市不透水面有效性的高效且精准的计算,将有助于提高常规城市不透水面的计算的精度。
下面以一个简明算例说明城市下垫面复杂空间特征对地表产汇流过程的影响。
如图2所示,图2(a)给出了一个九宫格的高程图,假设其为某城市区域的下垫面高程,计算水流方向,依据最陡坡度确定唯一流向,计算得到不同高程格点的水流流向,从图2(a)中可以看出,区域最终出水口为高程为240的格点,最后进行出口径流量和透水面有效性计算。图2(b)~(f)共给出了5种下垫面情景,图2(b)为100%透水面情景,图2(c)100%不透水面情景,图2(d)、(e)、(f)分别给出了2/3不透水面(不透水率66.67%)条件下的三种不同空间分布情景。
在仅考虑地表径流且忽略时间因素的情况下,假设各网格点均在同一时间获得1个单位降雨量,透水面径流系数为0.5,不透水面径流系数为1.0,采用径流系数法计算出水口最终径流量。计算结果为:
首先计算水流过程,根据图2(a)可以看出,最终的出口为高程240的格点,根据栅格间的坡度的表达式求得图2(a)中的9个格点与周围小于自身高程的格点之间的坡度。
例如:高程为324与303的格点之间的坡度为21:
高程为324与297的格点之间的坡度为19.09:
同理高程为324与312的格点之间坡度为12,所以最大下降方向为324→303,其他的格点之间同理依次求出最大下降方向,得出如图2(a)的水流过程。
本实施例中,100%透水面情景下,9个格点均为透水面,且每个格点的降水量相等,水流方向如图2(b)可知高程为324、350的格点经过4次透水面,高程为312、300、303的格点经过3次透水面,高程为297、260、322的格点经过2次透水面,高程为240的格点经过1次不透水面。根据径流系数法公式,可得:
F=2×α4+3×α3+3×α2+1×α
其中,α为0.5,最终求得径流量为1.75个单位。
同理可得100%不透水面情景下,径流量为9个单位;而不透水率为66.67%的三种不同空间分布情景下,径流量分别为:情景1(图2(d))为7.5个单位,情景2(图2(e))为3.375个单位,情景3(图2(f))为6.5个单位。
则三个情景的不透水面有效性计算分别为:
S1=(F1)/(F0)*100%=1
S2=(F2)/(F0)*100%=0.45
S3=(F3)/(F0)*100%=0.87
上从计算结果可知,地表径流量随着不透水率的增加而显著增加。同时,在同一不透水率下,下垫面空间分布特征不同,也会导致径流量差异显著。本算例中,在不透水率同为66.67%的三种不透水面空间分布情景下,最大径流量差异高达55%,对比100%不透水率和66.67%不透水率下的情景1,二者的径流量差异仅为16.67%,一定程度上表明,下垫面空间特征分布不同引起的径流量差异,甚至有可能比不透水率自身变化导致的差异更加剧烈。
Claims (3)
1.一种城市不透水面有效性量化方法,其特征在于,包括以下步骤:
S1、基于城市不透水面的空间分布,分别获取高程数据和土地利用数据;
S2、根据所述高程数据和土地利用数据,计算得到水流方向;
S3、根据所述水流方向和流径下垫面类型,计算得到研究区域出口最终径流量;
S4、根据所述研究区域出口最终径流量,计算得到城市不透水面有效性;
所述步骤S3中研究区域出口最终径流量的表达式如下:
其中,F表示研究区域出口最终径流量,Pi表示研究区域降雨量,a表示透水面径流系数,m表示水流经过透水面的个数,i表示研究区域的编号,i=1,2,3,...,n,n表示研究区域的总编号;
所述步骤S4中城市不透水面有效性的表达式如下:
其中,S表示城市不透水面有效性,F表示研究区域出口最终径流量,F0表示基于不透水率得到的出流量。
2.根据权利要求1所述的城市不透水面有效性量化方法,其特征在于,所述步骤S2包括以下步骤:
S201、根据所述高程数据和土地利用数据,针对每个栅格对周围八个方向均有高程差,且水流方向最陡时,计算中心栅格与相邻八个方向的栅格间的坡度;
S202、根据所述中心栅格与相邻八个方向的栅格间的坡度,确定每个栅格的最陡下降方向;
S203、判断所述每个栅格的最陡下降方向是否均相同,若是,则进入步骤S204,否则,进入步骤S205;
S204、扩大相邻栅格范围,并返回步骤S201;
S205、得到最大下降方向,并根据所述最大下降方向得到水流方向。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210182297.2A CN114547531B (zh) | 2022-02-25 | 2022-02-25 | 一种城市不透水面有效性量化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210182297.2A CN114547531B (zh) | 2022-02-25 | 2022-02-25 | 一种城市不透水面有效性量化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114547531A CN114547531A (zh) | 2022-05-27 |
CN114547531B true CN114547531B (zh) | 2022-09-06 |
Family
ID=81680350
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210182297.2A Active CN114547531B (zh) | 2022-02-25 | 2022-02-25 | 一种城市不透水面有效性量化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114547531B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116611249B (zh) * | 2023-05-23 | 2023-11-21 | 中国水利水电科学研究院 | 一种城市区域雨水入渗量的计算方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109086500A (zh) * | 2018-07-19 | 2018-12-25 | 武汉大学 | 基于空间分布式径流系数的未控区径流计算方法 |
CN109492259A (zh) * | 2018-10-15 | 2019-03-19 | 华北水利水电大学 | 一种城市水文模拟*** |
CN109490209A (zh) * | 2018-11-15 | 2019-03-19 | 中国科学院遥感与数字地球研究所 | 一种基于城市不透水面和dem的城市地表径流检测方法 |
CN109657841A (zh) * | 2018-11-27 | 2019-04-19 | 杭州师范大学 | 一种城市暴雨内涝积水深度提取方法 |
CN110232737A (zh) * | 2019-05-13 | 2019-09-13 | 杭州师范大学 | 一种城市汇水区划分方法 |
CN110457771A (zh) * | 2019-07-19 | 2019-11-15 | 河海大学 | 一种基于高程偏差传递的dem水流方向计算方法 |
CN112163346A (zh) * | 2020-10-21 | 2021-01-01 | 中国科学院生态环境研究中心 | 一种基于透水-不透水格局的面源污染迁移能力估算方法 |
-
2022
- 2022-02-25 CN CN202210182297.2A patent/CN114547531B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109086500A (zh) * | 2018-07-19 | 2018-12-25 | 武汉大学 | 基于空间分布式径流系数的未控区径流计算方法 |
CN109492259A (zh) * | 2018-10-15 | 2019-03-19 | 华北水利水电大学 | 一种城市水文模拟*** |
CN109490209A (zh) * | 2018-11-15 | 2019-03-19 | 中国科学院遥感与数字地球研究所 | 一种基于城市不透水面和dem的城市地表径流检测方法 |
CN109657841A (zh) * | 2018-11-27 | 2019-04-19 | 杭州师范大学 | 一种城市暴雨内涝积水深度提取方法 |
CN110232737A (zh) * | 2019-05-13 | 2019-09-13 | 杭州师范大学 | 一种城市汇水区划分方法 |
CN110457771A (zh) * | 2019-07-19 | 2019-11-15 | 河海大学 | 一种基于高程偏差传递的dem水流方向计算方法 |
CN112163346A (zh) * | 2020-10-21 | 2021-01-01 | 中国科学院生态环境研究中心 | 一种基于透水-不透水格局的面源污染迁移能力估算方法 |
Non-Patent Citations (4)
Title |
---|
不同空间布局对城市降雨径流的影响;金牧青 等;《水土保持通报》;20180430;第38卷(第02期);第33-39页 * |
城市典型下垫面产流过程模拟实验;刘慧娟 等;《资源科学》;20151130;第37卷(第11期);第2219-2227页 * |
基于DEM的城市降雨入渗条件调控模型的设计;李鹏 等;《滁州学院学报》;20091031;第11卷(第05期);第112-113页 * |
基于SWMM模型的不透水率与产汇流关系研究;李阳 等;《水电能源科学》;20170228;第35卷(第02期);第34-37页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114547531A (zh) | 2022-05-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106884405B (zh) | 一种无资料地区溃堤型山洪灾害分析评价方法 | |
CN111369059B (zh) | 基于雨洪模拟耦合模型的城市内涝快速预测方法及*** | |
Aher et al. | Quantification of morphometric characterization and prioritization for management planning in semi-arid tropics of India: a remote sensing and GIS approach | |
CN110232737B (zh) | 一种城市汇水区划分方法 | |
CN109871621A (zh) | 城市暴雨内涝汇水区分析方法 | |
CN110598290A (zh) | 考虑气候变化的流域未来水电发电能力预测方法和*** | |
CN111062125B (zh) | 海绵型综合管廊水文效应评估方法 | |
CN113011685A (zh) | 一种无径流资料地区内陆湖泊水位变化模拟预测方法 | |
Rao et al. | A distributed model for real-time flood forecasting in the Godavari Basin using space inputs | |
CN103093114A (zh) | 一种基于地形和土壤特性的分布式流域缺水量测算方法 | |
CN111538798B (zh) | 一种顾及dsm和dlg的城市汇水区精细化提取方法 | |
Wu et al. | Hydrological response to climate change and human activities: A case study of Taihu Basin, China | |
CN108009753A (zh) | 城市内涝预测方法、装置、终端及计算机可读存储介质 | |
CN114547531B (zh) | 一种城市不透水面有效性量化方法 | |
CN116305933B (zh) | 一种基于dem数据的坡面产汇流简易计算方法及装置 | |
CN106156949A (zh) | 一种基于内涝风险评估的雨水低影响开发规划方法 | |
CN114139258A (zh) | 城市内涝建模评估方法、***及计算机可读存储介质 | |
CN110377868B (zh) | 一种基于实时雨情的动态水系提取方法 | |
CN106096129B (zh) | 一种基于山地汇水计算的山脚水面规模分析方法 | |
CN104750936A (zh) | 一种基于gis的城市水文径流计算方法及*** | |
CN117195603B (zh) | 基于高分辨率遥感要素的洪涝灾害推演方法、设备及介质 | |
CN108920571B (zh) | 基于分布式水文模型的室内水文实验模型比例尺确定方法 | |
CN110109194A (zh) | 基于栅格陆面模型与向量化流径的河道栅格入流计算方法 | |
CN111178372B (zh) | 基于遥感影像和地形数据的大区域尺度黄土塬面提取方法 | |
Salem et al. | Quantification and evaluation of water erosion: Application of the Model SDR–InVEST in the Ziz Basin in South-East Morocco |
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 |