CN110543984B - 一种基于水安全的城市水面率规划方法 - Google Patents
一种基于水安全的城市水面率规划方法 Download PDFInfo
- Publication number
- CN110543984B CN110543984B CN201910793729.1A CN201910793729A CN110543984B CN 110543984 B CN110543984 B CN 110543984B CN 201910793729 A CN201910793729 A CN 201910793729A CN 110543984 B CN110543984 B CN 110543984B
- Authority
- CN
- China
- Prior art keywords
- water
- area
- rainwater
- maximum
- rainfall
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 180
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000004364 calculation method Methods 0.000 claims abstract description 34
- 238000013461 design Methods 0.000 claims abstract description 8
- 238000005192 partition Methods 0.000 claims description 33
- 238000009825 accumulation Methods 0.000 claims description 28
- 238000004458 analytical method Methods 0.000 claims description 18
- 238000007726 management method Methods 0.000 claims description 11
- 244000025254 Cannabis sativa Species 0.000 claims description 6
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 claims description 6
- 230000001174 ascending effect Effects 0.000 claims description 6
- 230000001788 irregular Effects 0.000 claims description 6
- 230000035699 permeability Effects 0.000 claims description 6
- 238000004088 simulation Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 4
- 230000002265 prevention Effects 0.000 claims description 4
- 208000006011 Stroke Diseases 0.000 claims description 3
- 125000004432 carbon atom Chemical group C* 0.000 claims description 3
- 238000000605 extraction Methods 0.000 claims description 3
- 235000021384 green leafy vegetables Nutrition 0.000 claims description 3
- 238000013507 mapping Methods 0.000 claims description 3
- 238000000691 measurement method Methods 0.000 claims description 3
- 238000005457 optimization Methods 0.000 claims description 3
- 238000011160 research Methods 0.000 claims description 3
- 239000011850 water-based material Substances 0.000 claims description 3
- 238000012821 model calculation Methods 0.000 abstract description 2
- 238000010276 construction Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000012502 risk assessment Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 235000006506 Brasenia schreberi Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000000556 factor analysis Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000012876 topography Methods 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/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
-
- 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
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Tourism & Hospitality (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Theoretical Computer Science (AREA)
- Development Economics (AREA)
- General Physics & Mathematics (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Sewage (AREA)
Abstract
本发明公开了一种基于水安全的城市水面率规划方法，包括：测量DEM地形高程数据；划定雨水管控分区；建立水力模型，识别积水区域；通过积水面积比例，定义每个雨水管控分区内涝风险等级；通过地区暴雨参数图集查算暴雨径流量；计算蓄水容积与水面率。本发明充分考虑各雨水管控分区积水区域面积、深度和内涝风险，根据城市暴雨数据和城市蓄水量，计算得出基于水安全的最大水域面积及水面率，可操作性强，结果合理，极具参考性。以模型计算、雨洪计算、水位安全计算为主要手段，作为城市水安全布局规划、规模设计的重要支撑，提高了城市水安全规划的科学性。
Description
技术领域
本发明涉及城市规划的技术领域，尤其是指一种基于水安全的城市水面率规划方法。
背景技术
城市水面率是通常指城市一定区域范围内河流(江、河、渠等)、湖泊(天然或人工湖泊)、湿地(天然或人工湿地)等承载水域功能的区域面积占区域总面积的比率。城市水面具有多种服务功能和经济价值(何俊仕,吴迪,魏国.城市适宜水面率及其影响因素分析[J].干旱区资源与环境)，城市水面率是影响城市生态系统、防洪排涝系统的关键指标，对城市防灾能力、环境品质、经济发展等各方面具有重要影响。目前已有水面率判定方法主要包括：
(1)试做方案法：通过设计多个土地利用方案，在方案对比过程中兼顾、平衡土地开发与城市水域功能之间的空间关系，确定方案后量算城市水面率数值。该方法在城市水面率判定的实践中较多采用，其判定结果存在较大的主观性，受方案对比过程专业融合程度及试做方案者的经验和水平。
(2)经验判断法。主要由相关专家和专业技术人员通过借鉴相关城市经验，考虑城市开发需求、城市景观建设需求、现状水面率等情况直接估算。该方法在城市水面率应用较少且估算的主观性较强，估算的准确性受专业技术人员配置、水平及案例选择影响较大。
(3)规范标准。《城市水系规划规范》(GB 50513-2009)(2016年版))以省/自治区为单位，将全国划为三个分区适宜水面率方向性意见。该判定方法存在问题包括：1)针对性不足，各地区下垫面差异较大，部分地区难以满足规范要求；2)难以确定水面率统计的单元边界，实践中受边界范围线影响较大；3)科学性不足。该方法在实践中应用较少，指导性不强。
此外，赵璧奎,黄本胜,邱静,等(海绵城市建设中区域适宜水面率研究及应用[J].广东水利水电)以海绵城市径流控制的指标计算城市适宜水面率[3]，通过确定河湖水系水位涨幅目标从而计算各个雨水管控分区的适宜水面率。该方法仅从城市防洪排涝及水安全保障的角度明确了雨水管控分区的适宜水面率域值，而对于城市水安全考虑不周。
由上可见，现有技术中城市水面率的计算方法均存在以下问题：
1、水面率受主观和经验的影响很大；
2、现有技术中没有有效手段以城市水安全为目的对城市水面率进行规划。
发明内容
本发明的目的在于克服现有技术不足，提出了一种基于水安全的城市水面率规划方法，通过各雨水管控分区的水域面积的计算为城市水面规划提出解决方案。
为实现上述目的，本发明所提供的技术方案为：一种基于水安全的城市水面率规划方法，包括以下步骤：
1)测量DEM地形高程数据；
2)建立水力模型，识别积水区域；
3)通过积水面积比例，定义每个雨水管控分区内涝风险等级；
4)通过地区暴雨参数图集查算暴雨径流量；
5)计算蓄水容积与水面率。
在步骤1)中，采用大地测量方法测量城市规划区域地形高程数据，测图精度不低于1：200，成图精度不低于1：500；测图高程精度不低于水准V。
在步骤2)中，建立水力模型，识别积水区域，具体过程如下：
以DEM地形高程数据为基础，通过GIS和计算水力学软件对数据进行转换，并建立地形的闭边界，使模型降雨落入边界范围内，避免模型发散；
设计降雨：采用2h降雨或24h降雨；2h降雨主要针对小排水系统规划，用芝加哥法推求；24h降雨主要针对大排水系统规划，用同频率分析方法推求；
根据地区标准设计雨型及降雨数据制作时间序列文件，模型模拟不受雨峰系数的影响，在推求设计降雨时能够做成前峰型降雨，也能够做成后峰型降雨；
设计一种综合糙率的估算方法，使得在计算中快速逼近真实糙率；采用的估算方法包括以下步骤：
2.1)划定雨水管控分区：
2.1.1)利用步骤1)测量得到的DEM地形高程数据，基于GIS的3D分析模块进行不规则三角网(TIN)的建立，并转化为网格型连续高程表面，得到所需的DEM数据；
2.1.2)以DEM数据为基础，采用GIS的水文分析模块，通过洼地填充——水流方向计算——流量计算——汇流提取——汇点识别——流域分区，形成初步的模拟汇流路径与流域分区边界；
2.1.3)输入现状河流水系及已确定的规划河流水系数据，重新优化水文分析模型，重新生成汇流路径，重新生成流域分区边界；
2.1.4)输入不可穿越边界，包括：堤围、铁路、高速公路，定义汇流路径不得穿越上述设施，进一步优化水文分析模型，调整汇流路径与流域分区的边界；其中，所述的进一步优化水文分析模型是指：设置铁路地基的透水率和区域汇流时间及设置高速公路地基的透水率和区域汇流时间；
2.1.5)从行政管理角度，结合道路路网、行政边界、规划管理单元边界调整分区边界，形成最终的雨水管控分区，单个分区面积为1-5平方公里；
2.2)在每个雨水管控分区内获取所有雨水排水口的最大汇水距离，具体是以雨水排水口为圆心及以DEM网格最小尺度为半径，向外侧做圆，直到某个圆与地形高程的交点从上升趋势转为下降趋势位置，选择上升与下降之间的拐点作为该雨水排水口的最大汇水距离；而后测量拐点到雨水排水口的距离和高差计算其坡度；
2.3)计算所有雨水排水口的坡度，取其平均值作为糙率计算坡度Ja；
2.4)概化水力半径Ra的取值，分为两类：绿地和硬化地面；
对于绿地取为植草高度的0.5-1.5倍，对于硬化地面取值为公路路沿高度的一半；
其中居住用地、商业用地、工业用地、公用设施用地均采用其绿化率与硬化地面的加权值，对于道路广场用地采用硬化地面取值，即为公路路沿高度的一半；对于公园绿地取值为公园绿地植草高度的0.5-1.5倍；将计算区域的所有面积加权平均得到概化水力半径；
综合糙率n的计算公式为：
确定综合径流系数：研究区域的流域范围内综合径流系数按照各下垫面种类径流系数加权平均计算得出；
上式中,
为综合径流系数，
为不同下垫面对应径流系数，Fi为不同下垫面面积，F为流域面积；
管道排水能力概化值的计算方法：雨水管道排除能力值通过概化计算方式获取：
上式中，i为雨水管道排除能力，mm/min；s为汇水面积，hm2；Q为排水管渠的设计流量，m3/s；v为水流断面的平均流速，m/s；n为粗糙系数；A为管渠断面面积，m2；R为管渠按满流设计的水力半径，m；I为出水口管道坡度；
依据上述数据，可利用计算水力学软件建立水力学概化模型，并模拟最大积水深度图，提取积水深度大于0.80m的积水面积。
在步骤3)中，积水深度大于0.80m的积水点为严重积水；计算出每个雨水管控分区严重积水区域的面积占比；
上式中,Pi为不同雨水管控分区严重积水区域面积占比，Ji为不同雨水管控分区严重积水区域面积，Si为对应雨水管控分区面积；
定义内涝风险区：
高风险区域：雨水管控分区内严重积水区域占比≥30％；
中风险区域：雨水管控分区内严重积水区域占比为10％-30％；
低风险区域：雨水管控分区内严重积水区域占比＜10％。
在步骤4)中，根据内涝防治标准与地区暴雨统计参数图集，通过地区最大1h、2h、6h点雨量均值等值线图及地区最大1h、2h、6h点雨量变差系数Cv等值线图，获取地区最大1h、2h、6h点雨量均值及对应变差系数Cv，利用皮尔逊Ⅲ型频率曲线获取模比系数Kp值，计算地区最大1h、2h、6h的暴雨量，并根据综合径流系数，计算每个雨水管控分区最大1h、2h、6h的暴雨径流量分别如式(6)、式(7)、式(8)所示；
上式中,L1、L2、L6为最大1h、2h、6h的暴雨径流量，Hp1、Hp2、Hp6为最大1h、2h、6h的暴雨点流量，Kp1、Kp2、Kp6为最大1h、2h、6h的对应模比系数，
为综合径流系数，F为流域面积；
在步骤5)中，蓄水容积的计算方法为：
高风险区：蓄存容积为最大6小时暴雨径流量，即为式(8)计算结果；
中风险区：蓄存容积为最大2小时暴雨径流量，即为式(7)计算结果；
低风险区：蓄存容积为最大1小时暴雨径流量，即为式(6)计算结果；
确定河湖水系水位安全变幅值ΔH：河湖水系水位安全变幅值ΔH受设计雨量、综合径流系数、地区河湖水系常水位以上的水位安全调蓄变幅的相互影响，取值为0.5-0.8m；
根据蓄水总容积V及水位安全变幅值ΔH，计算各管控分区水域面积：
上式中,Ai水为各雨水管控分区水域面积也就是规划水面面积，其分布也就是城市规划水面，Vi为各雨水管控分区积水容积，ΔH为河湖水系水位安全变幅值；
水域面积与雨水管控分区总面积的比值，即为水面率。
本发明与现有技术相比，具有如下优点与有益效果：
1、本发明充分考虑各雨水管控分区积水区域面积、深度和内涝风险，根据城市暴雨数据和城市蓄水量，计算得出基于水安全的最大水域面积及水面率，可操作性强，结果合理，极具参考性。
2、以模型计算、雨洪计算、水位安全计算为主要手段，作为城市水安全布局规划、规模设计的重要支撑，提高了城市水安全规划的科学性。
附图说明
图1为本发明实施例模拟区域地形图。
图2为本发明实施例雨水管控分区图。
图3为本发明实施例5年一遇2h历时降雨时间序列图。
图4为本发明实施例模拟最大积水深度图。
图5为本发明方法的流程示意图。
具体实施方式
下面以一个具体实施案例对本发明的技术方案进行详细说明。
如附图5所示，某地区基于水安全的城市水面率规划包括以下步骤：
步骤1、获取该地区DEM高程数据
采用大地测量方法测量城市规划区域地形高程数据，测图精度不低于1：200，成图精度不低于1：500；测图高程精度不低于水准V。
步骤2、建立水力模型，识别积水区域；
以DEM地形高程数据为基础，通过GIS和计算水力学软件对数据进行转换，并建立地形的闭边界；
设计降雨采用5年一遇2h历时降雨，根据地区标准设计雨型及降雨数据，用芝加哥法推求制作时间序列文件(dfs0)。模型模拟不受雨峰系数的影响，在推求设计降雨时可以做成前峰型降雨，也可以做成后峰型降雨，如附图3。
设计一种综合糙率的估算方法，使得在计算中快速逼近真实糙率；采用的估算方法包括以下步骤：
2.1)划定该地区的雨水管控分区
2.1.1)利用上述测量过程得到的DEM地形高程数据，基于GIS的3D分析模块进行不规则三角网(TIN)的建立，并转化为网格型连续高程表面(DEM)，如附图1；
2.1.2)以DEM数据为基础，采用GIS的水文分析模块，通过洼地填充——水流方向计算——流量计算——汇流提取——汇点识别——流域分区，形成初步的模拟汇流路径与流域分区边界；
2.1.3)地形数据生成的汇流与实际河流有一定的差别，城市规划中河网水系也会有一定的调整。输入现状河流水系或已确定的规划河流水系等重要排水通道数据，在同一流域输入数据替代初步生成的模拟汇流，优化水文分析模型，通过新生成的汇流路径重新生成流域分区；
2.1.4)城市中人工环境的干预对场地雨水汇流有直接的影响，输入堤围、铁路、高速公路(非架空)等数据，定义汇流路径不得穿越上述设施，结合实际情况人工进一步优化水文分析模型，调整汇流路径与流域分区的边界；
输入不可穿越边界，包括：堤围、铁路、高速公路，定义汇流路径不得穿越上述设施，进一步优化水文分析模型，调整汇流路径与流域分区的边界；上述的进一步优化水文分析模型是指：a.设置铁路地基的透水率和区域汇流时间；b.设置高速公路地基的透水率和区域汇流时间；
2.1.5)从行政管理角度，结合道路路网、行政边界、规划管理单元边界适当调整分区边界，形成最终的雨水管控分区，如附图2，单个分区面积为1-5平方公里。
雨水管控分区信息表
| 管控分区序号 | 面积(km<sup>2</sup>) | 综合径流系数 | 管控分区序号 | 面积(km<sup>2</sup>) | 综合径流系数 |
| 1 | 4.26 | 0.67 | 4 | 2.23 | 0.61 |
| 2 | 3.52 | 0.23 | 5 | 1.67 | 0.38 |
| 3 | 1.06 | 0.34 | 6 | 3.26 | 0.41 |
2.2)在每个雨水管控分区内获取所有雨水排水口的最大汇水距离，其具体方法为以雨水排水口为圆心以DEM网格最小尺度为半径，向外侧做圆，直到某个圆与地形高程的交点从上升趋势转为下降趋势位置，选择上升与下降之间的拐点作为该雨水排水口的最大汇水距离；测量拐点到雨水排水口的距离和高差计算其坡度；
2.3)计算所有雨水排水口的坡度，取其平均值作为糙率计算坡度Ja＝0.005；
2.4)概化水力半径Ra的取值分为两类：绿地和硬化地面；
对于绿地取为植草高度的1倍为0.07m，对于硬化地面取值为公路路沿高度的一半为0.125m；
其中居住用地、商业用地、工业用地、公用设施用地均采用其绿化率与硬化地面的加权值，对于道路广场用地采用硬化地面取值，即为公路路沿高度的一半；对于公园绿地取值为公园绿地植草高度的0.5-1.5倍；将计算区域的所有面积加权平均得到概化水力半径Ra＝0.276；
综合糙率n的计算公式为：
确定综合径流系数：
广州市某地区的流域范围内综合径流系数按照各下垫面种类径流系数加权平均计算得出。
上式中,
为综合径流系数，
为不同下垫面对应径流系数，Fi为不同下垫面面积，F为流域面积。建设用地径流系数：中高密度居住用地：0.55；商业用地：0.75；工业用地、公用设施用地：0.65；道路广场用地：0.8；公园绿地0.05。
综合径流系数计算表
管道排水能力概化值的计算方法：
雨水管道排除能力值通过概化计算方式获取：
Q＝vA (3)
上式中，i为雨水管道排除能力，mm/min；s为汇水面积，hm2；Q为排水管渠的设计流量，m3/s；v为水流断面的平均流速，m/s；n为粗糙系数；A为管渠断面面积，m2；R为管渠按满流设计的水力半径，m；I为出水口管道坡度。
采用历史淹水位置处管网参数，v取1m/s,计算雨水排除能力值，i＝0.3mm/min。
依据上述数据，可利用计算水力学软件建立水力学概化模型，并模拟最大积水深度图，提取积水深度大于0.80m的积水面积，见附图4。
步骤3、计算每个雨水管控分区严重积水区域面积占比
积水深度大于0.80m的积水点为严重积水。可计算出每个雨水管控分区严重积水区域的面积占比。
上式中,Pi为不同雨水管控分区严重积水区域面积占比，Ji为不同雨水管控分区严重积水区域面积，Si为对应雨水管控分区面积。
以管控分区1为例：
各雨水管控分区严重积水区域面积占比
步骤4：定义内涝风险区
高风险区域：雨水管控分区内严重积水区域占比≥20％；
中风险区域：雨水管控分区内严重积水区域占比为10％-20％；
低风险区域：雨水管控分区内严重积水区域占比＜10％。
各雨水管控分区内涝风险评估
| 管控分区序号 | 风险评估 | 管控分区序号 | 风险评估 |
| 1 | 低风险区域 | 4 | 低风险区域 |
| 2 | 高风险区域 | 5 | 中风险区域 |
| 3 | 中风险区域 | 6 | 中风险区域 |
步骤5：通过地区暴雨统计参数图集查算暴雨径流量
根据内涝防治标准与地区暴雨统计参数图集，通过地区最大1h、2h、6h点雨量均值等值线图及地区最大1h、2h、6h点雨量变差系数Cv等值线图，获取地区最大1h、2h、6h点雨量均值及对应变差系数Cv，利用皮尔逊Ⅲ型频率曲线获取模比系数Kp值，计算地区最大1h、2h、6h的暴雨量，并根据步骤1计算的径流系数，计算每个雨水管控分区最大1h、2h、6h的暴雨径流量分别如式(6)、式(7)、式(8)所示。
上式中,L1、L2、L6为最大1h、2h、6h的暴雨径流量，Hp1、Hp2、Hp6为最大1h、2h、6h的暴雨点流量，Kp1、Kp2、Kp6为最大1h、2h、6h的对应模比系数，
为综合径流系数，Yi为各雨水管控分区面积。
步骤6：计算蓄水总容积目标
高风险区：蓄存最大6小时暴雨径流量，即为式8计算结果。
中风险区：蓄存最大2小时暴雨径流量，即为式7计算结果。
低风险区：蓄存最大1小时暴雨径流量，即为式6计算结果。
以管控分区1为例：
各雨水管控分区蓄水总容积
步骤7、计算蓄水容积与水面率
河湖水系水位安全变幅值ΔH受设计雨量、综合径流系数、地区河湖水系常水位以上的水位安全调蓄变幅的相互影响，取值为0.5-0.8m。
步骤8、计算每个管控分区最大水域面积
根据蓄水总容积V及水位安全变幅值ΔH，计算各管控分区的最大水域面积。
上式中,Ai水为各雨水管控分区水域面积也就是规划水面面积，Vi为各雨水管控分区积水容积，ΔH为河湖水系水位安全变幅值取0.7。
以管控分区1为例：
各雨水管控分区水面率
由于各管控分区的空间位置是已知的，因此获得城市中不同空间位置的水面分布情况，也就是获得了基于水安全的城市水面规划方案。
以上所述实施例只为本发明之较佳实施例，并非以此限制本发明的实施范围，故凡依本发明之形状、原理所作的变化，均应涵盖在本发明的保护范围内。
Claims (2)
1.一种基于水安全的城市水面率规划方法，其特征在于，包括以下步骤：
1)测量DEM地形高程数据；
2)建立水力学概化模型，识别积水区域，具体过程如下：
以DEM地形高程数据为基础，通过GIS和计算水力学软件对数据进行转换，并建立地形的闭边界，使水力学概化模型降雨落入边界范围内；
设计降雨：采用2h降雨或24h降雨；2h降雨针对小排水系统规划，用芝加哥法推求；24h降雨针对大排水系统规划，用同频率分析方法推求；
根据地区标准设计雨型及降雨数据制作时间序列文件，水力学概化模型模拟不受雨峰系数的影响，在推求设计降雨时，设计为前峰型降雨或后峰型降雨；
设计一种综合糙率的估算方法，采用的估算方法包括以下步骤：
2.1)划定雨水管控分区
2.1.1)利用步骤1)测量得到的DEM地形高程数据，基于GIS的3D分析模块进行不规则三角网TIN的建立，并转化为网格型连续高程表面，得到所需的DEM地形高程数据；
2.1.2)以DEM地形高程数据为基础，采用GIS的水文分析模块，依次通过洼地填充、水流方向计算、流量计算、汇流提取、汇点识别和流域分区，形成初步的模拟汇流路径与流域分区边界；
2.1.3)输入现状河流水系及已确定的规划河流水系数据，重新优化水文分析模型，重新生成汇流路径，重新生成流域分区边界；
2.1.4)输入不可穿越边界，包括堤围、铁路、高速公路设施，定义汇流路径不得穿越上述设施，优化水文分析模型，调整汇流路径与流域分区的边界；其中，优化水文分析模型是指设置铁路地基的透水率和区域汇流时间及设置高速公路地基的透水率和区域汇流时间；
2.1.5)结合道路路网、行政边界、规划管理单元边界调整分区边界，形成最终的雨水管控分区，单个分区面积为1-5平方公里；
2.2)在每个雨水管控分区内获取所有雨水排水口的最大汇水距离，具体是以雨水排水口为圆心及以DEM网格最小尺度为半径，向外侧做圆，直到某个圆与地形高程的交点从上升趋势转为下降趋势位置，选择上升与下降之间的拐点作为该雨水排水口的最大汇水距离；而后测量拐点到雨水排水口的距离和高差计算其坡度；
2.3)计算所有雨水排水口的坡度，取其平均值作为糙率计算坡度Ja；
2.4)概化水力半径Ra的取值，分为两类：绿地和硬化地面；
对于绿地取为植草高度的0.5-1.5倍，对于硬化地面取值为公路路沿高度的一半；
其中，居住用地、商业用地、工业用地、公用设施用地均采用绿化率与硬化地面的加权值，对于道路广场用地采用硬化地面取值，即为公路路沿高度的一半；对于公园绿地取值为公园绿地植草高度的0.5-1.5倍；将计算区域的所有面积加权平均得到概化水力半径；
综合糙率n的计算公式为：
确定综合径流系数：研究区域的流域范围内综合径流系数按照各下垫面种类径流系数加权平均计算得出；
上式中,
为综合径流系数，
为不同下垫面对应径流系数，Fi为不同下垫面面积，F为流域面积；
管道排水能力概化值的计算方法：雨水管道排除能力值通过概化计算方式获取：
Q＝vA (3)
上式中，i为雨水管道排除能力，单位为mm/min；s为汇水面积，单位为hm2；Q为排水管渠的设计流量，单位为m3/s；v为水流断面的平均流速，单位为m/s；A为管渠断面面积，单位为m2；R为管渠按满流设计的水力半径，单位为m；I为出水口管道坡度；
依据上述数据，利用计算水力学软件建立水力学概化模型，并模拟最大积水深度图，提取积水深度大于0.80m的积水面积；
3)通过积水面积比例，定义每个雨水管控分区内涝风险等级，具体如下：
积水深度大于0.80m的积水面积为严重积水，计算出每个雨水管控分区严重积水区域的面积占比；
上式中,Pi为不同雨水管控分区严重积水区域面积占比，Ji为不同雨水管控分区严重积水区域面积，Si为对应雨水管控分区面积；
定义内涝风险区：
高风险区域：雨水管控分区内严重积水区域占比≥30％；
中风险区域：雨水管控分区内严重积水区域占比为10％-30％；
低风险区域：雨水管控分区内严重积水区域占比＜10％；
4)通过地区暴雨参数图集查算暴雨径流量，具体如下：
根据内涝防治标准与地区暴雨统计参数图集，通过地区最大1h、2h、6h点雨量均值等值线图及地区最大1h、2h、6h点雨量变差系数Cv等值线图，获取地区最大1h、2h、6h点雨量均值及对应变差系数Cv，利用皮尔逊Ⅲ型频率曲线获取模比系数Kp值，计算地区最大1h、2h、6h的暴雨量，并根据综合径流系数，计算每个雨水管控分区最大1h、2h、6h的暴雨径流量分别如式(6)、式(7)、式(8)所示；
上式中,L1为最大1h的暴雨径流量，L2为最大2h的暴雨径流量，L6为最大6h的暴雨径流量，Hp1为最大1h的暴雨点流量，Hp2为最大2h的暴雨点流量，Hp6为最大6h的暴雨点流量，Kp1为最大1h的对应模比系数，Kp2为最大2h的对应模比系数，Kp6为最大6h的对应模比系数，
为综合径流系数，F为流域面积；
5)计算蓄水容积与水面率，具体如下：
蓄水容积的计算方法为：
高风险区：蓄存容积为最大6小时暴雨径流量，即为式(8)计算结果；
中风险区：蓄存容积为最大2小时暴雨径流量，即为式(7)计算结果；
低风险区：蓄存容积为最大1小时暴雨径流量，即为式(6)计算结果；
确定河湖水系水位安全变幅值ΔH：河湖水系水位安全变幅值ΔH受设计雨量、综合径流系数、地区河湖水系常水位以上的水位安全调蓄变幅的相互影响，取值为0.5-0.8m；
根据蓄水总容积V及水位安全变幅值ΔH，计算各管控分区水域面积：
上式中,Ai水为各雨水管控分区水域面积也就是规划水面面积，其分布也就是城市规划水面，Vi为各雨水管控分区积水容积，ΔH为河湖水系水位安全变幅值；
水域面积与雨水管控分区总面积的比值，即为水面率。
2.根据权利要求1所述的一种基于水安全的城市水面率规划方法，其特征在于：在步骤1)中，采用大地测量方法测量城市规划区域地形高程数据，测图精度不低于1：200，成图精度不低于1：500；测图高程精度不低于水准V。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910793729.1A CN110543984B (zh) | 2019-08-27 | 2019-08-27 | 一种基于水安全的城市水面率规划方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910793729.1A CN110543984B (zh) | 2019-08-27 | 2019-08-27 | 一种基于水安全的城市水面率规划方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110543984A CN110543984A (zh) | 2019-12-06 |
| CN110543984B true CN110543984B (zh) | 2022-06-14 |
Family
ID=68712132
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910793729.1A Active CN110543984B (zh) | 2019-08-27 | 2019-08-27 | 一种基于水安全的城市水面率规划方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110543984B (zh) |
Families Citing this family (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111046551B (zh) * | 2019-12-10 | 2020-11-20 | 河海大学 | 一种城市群排水过程模拟方法 |
| CN111881537B (zh) * | 2020-08-10 | 2024-02-20 | 哈尔滨工业大学 | 一种基于灰绿融合的海绵建设效果评价方法 |
| CN112287285B (zh) * | 2020-10-29 | 2023-04-04 | 南昌工程学院 | 基于城市内涝基础信息的蓄涝水面率分区处理方法 |
| CN112942230B (zh) * | 2021-01-21 | 2022-12-27 | 中城(广州)城乡规划设计有限公司 | 面向水污染控制的城市人工湖水域面积值域计算方法 |
| CN113023958B (zh) * | 2021-04-09 | 2022-07-26 | 山东速力动力科技有限公司 | 一种建筑工地雨水净化循环处理系统 |
| CN113449983B (zh) * | 2021-06-25 | 2024-01-05 | 中国水利水电科学研究院 | 一种新型城市湿地规划设计方法 |
| CN113592330B (zh) * | 2021-08-06 | 2023-09-12 | 中国民航大学 | 一种基于飞机起降安全的机场跑道评价方法 |
| CN114969659A (zh) * | 2022-05-13 | 2022-08-30 | 中国电建集团贵阳勘测设计研究院有限公司 | 一种基于水文拐点的水华防治生态调控流量计算方法 |
| CN115374638B (zh) * | 2022-08-24 | 2023-05-30 | 中国水利水电科学研究院 | 一种基于径流模拟的城市不透水面空间管控方法 |
| CN116468269B (zh) * | 2023-04-10 | 2023-10-20 | 广州市城市规划勘测设计研究院 | 一种洪涝高风险区识别方法、装置、设备及存储介质 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109242355A (zh) * | 2018-10-24 | 2019-01-18 | 上海市政工程设计研究总院（集团）有限公司 | 控制性详细规划阶段海绵城市低影响开发指标分解方法 |
| CN109409694A (zh) * | 2018-09-30 | 2019-03-01 | 昆山市建设工程质量检测中心 | 基于实测值的地块海绵城市绩效考核指标计算方法 |
| CN109948866A (zh) * | 2019-04-11 | 2019-06-28 | 沈阳建筑大学 | 一种基于swwm模拟的北方滨海地区海绵城市规划方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20180017710A1 (en) * | 2016-07-18 | 2018-01-18 | 2NDNATURE Software Inc. | Systems and Methods for Event-based Modeling of Runoff and Pollutant Benefits of Sustainable Stormwater Management |
-
2019
- 2019-08-27 CN CN201910793729.1A patent/CN110543984B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109409694A (zh) * | 2018-09-30 | 2019-03-01 | 昆山市建设工程质量检测中心 | 基于实测值的地块海绵城市绩效考核指标计算方法 |
| CN109242355A (zh) * | 2018-10-24 | 2019-01-18 | 上海市政工程设计研究总院（集团）有限公司 | 控制性详细规划阶段海绵城市低影响开发指标分解方法 |
| CN109948866A (zh) * | 2019-04-11 | 2019-06-28 | 沈阳建筑大学 | 一种基于swwm模拟的北方滨海地区海绵城市规划方法 |
Non-Patent Citations (3)
| Title |
|---|
| The Response of DEM Reconstruction Algorithm on Simulating Rainfall-runoff Model;Xiaoli Li et al.;《2018 5th IEEE International Conference on Cloud Computing and Intelligence》;20181123;全文 * |
| 典型下垫面径流系数预测的神经网络方法研究;武晟 等;《环境科学与技术》;20070531;第30卷(第5期);全文 * |
| 海绵城市建设中区域适宜水面率研究及应用;赵璧奎 等;《广东水利水电》;20170531(第5期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110543984A (zh) | 2019-12-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110543984B (zh) | 一种基于水安全的城市水面率规划方法 | |
| Hassan et al. | Climate change effect on storm drainage networks by storm water management model | |
| CN107832931A (zh) | 一种平原水网地区内涝风险的模块化分析方法 | |
| CN108984823B (zh) | 一种合流制溢流调蓄池规模的确定方法 | |
| CN110188323B (zh) | 一种适宜于城市水生态保护的水面率值域测算方法 | |
| CN106779232B (zh) | 一种城市内涝的建模预测方法 | |
| CN110232737B (zh) | 一种城市汇水区划分方法 | |
| CN111581767B (zh) | 一种管网-河流耦合模型校核特征参数率定方法 | |
| CN115391712A (zh) | 一种城市洪涝风险预测方法 | |
| CN111062125B (zh) | 海绵型综合管廊水文效应评估方法 | |
| CN112052561A (zh) | 一种排水系统防涝应急预案制定方法 | |
| CN115859676B (zh) | 一种顾及气候要素的多层次城市内涝耦合模拟方法 | |
| CN112242003B (zh) | 一种顾及地类与流向的城市子汇水区划分方法 | |
| CN110287595B (zh) | 一种城市不同下垫面减灾效应分析方法 | |
| CN106156949A (zh) | 一种基于内涝风险评估的雨水低影响开发规划方法 | |
| CN116561957B (zh) | 一种基于数值模拟的地下空间内涝评估方法 | |
| Ben-Daoud et al. | The use of low impact development technologies in the attenuation of flood flows in an urban area: Settat city (Morocco) as a case | |
| CN115238444A (zh) | 耦合swmm与gast模型的城市雨洪模拟方法 | |
| CN104750936A (zh) | 一种基于gis的城市水文径流计算方法及系统 | |
| CN114676473A (zh) | 一种基于人工智能算法的绿色基础设施空间布局优化方法 | |
| CN113742896A (zh) | 一种城市尺度的浅表流雨水排水系统的构建方法 | |
| CN112663760A (zh) | 基于lid措施的城市雨水径流调控方法 | |
| CN109033589A (zh) | 一种基于swmm模型考虑lid净化作用的污染物去除率计算方法 | |
| CN111985129A (zh) | 一种城市暴雨内涝精细化模拟方法 | |
| CN115563740B (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 |