CN112711825B - 一种地表与管网分布式的直接双向耦合方法 - Google Patents

一种地表与管网分布式的直接双向耦合方法 Download PDF

Info

Publication number
CN112711825B
CN112711825B CN202110323521.0A CN202110323521A CN112711825B CN 112711825 B CN112711825 B CN 112711825B CN 202110323521 A CN202110323521 A CN 202110323521A CN 112711825 B CN112711825 B CN 112711825B
Authority
CN
China
Prior art keywords
pipe network
coupling
model
drainage
water
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
Application number
CN202110323521.0A
Other languages
English (en)
Other versions
CN112711825A (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.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
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 China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN202110323521.0A priority Critical patent/CN112711825B/zh
Publication of CN112711825A publication Critical patent/CN112711825A/zh
Application granted granted Critical
Publication of CN112711825B publication Critical patent/CN112711825B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/10—Geometric CAD
    • G06F30/18—Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • E—FIXED CONSTRUCTIONS
    • E03—WATER SUPPLY; SEWERAGE
    • E03F—SEWERS; CESSPOOLS
    • E03F1/00—Methods, systems, or installations for draining-off sewage or storm water
    • E—FIXED CONSTRUCTIONS
    • E03—WATER SUPPLY; SEWERAGE
    • E03F—SEWERS; CESSPOOLS
    • E03F3/00—Sewer pipe-line systems
    • E03F3/02—Arrangement of sewer pipe-lines or pipe-line systems
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Public Health (AREA)
  • Health & Medical Sciences (AREA)
  • Hydrology & Water Resources (AREA)
  • Water Supply & Treatment (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Sewage (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供一种地表与管网分布式的直接双向耦合方法，包括获取地表水文水动力数据，还包括以下步骤：依据所述地表水文水动力数据生成地表模型和管网模型，所述地表模型为地表分布式水文水动力模型；在所述地表模型和所述管网模型的基础上，根据城市排水管网实际组成以及数据覆盖情况，构建地表模型与管网模型之间的耦合模式，所述耦合模式至少包括物理耦合和概念耦合。本发明提出了一种新的分布式的直接双向耦合模式，即降雨（Rainfall）—地表（Overland）—管网（Sewer）的耦合模式，在该模式下，地表采用分布式的水文水动力学模型，降雨直接作用在地表栅格上，降雨产汇流和积水淹没均基于统一的地表栅格进行计算，同时地表模型与管网模型进行双向耦合。

Description

一种地表与管网分布式的直接双向耦合方法
技术领域
本发明涉及市政工程的技术领域，特别是一种地表与管网分布式的直接双向耦合方法。
背景技术
城市暴雨洪涝模拟过程中，地表模型与管网模型的耦合计算一直是关键点和技术难点之一。现有的城市暴雨洪水模拟模型，主要采用降雨（Rainfall）—子流域（Catchment）—管网（Sewer）—地表（Overland）（以下简称RCSO）的模式。按照这种模式，首先须根据排水管网系统中检查井的位置划分子流域，然后将降雨直接作用在子流域上，再通过子流域降雨产汇流模型计算子流域上的产流量，子流域上的产流量直接通过检查井加入管网模型。当管网检查井出现溢流后，驱动地表淹没模型，并将检查井溢流量加入到地表淹没模型中，再通过地表淹没模型计算地表洪水淹没情况。
上述RCSO模式，地表和管网之间不是直接的双向的耦合，这对于评估管网的运行效率是可行且有效的，但对于评估地表洪水淹没却存在明显的缺点。首先，该模式将城市降雨和淹没分割成两个独立的过程：子流域上的降雨产汇流过程与地表洪水淹没过程，这与城市降雨引起局地积水的实际过程是不相符的。实际上，城市降雨产流过程和积水淹没过程，是作用在城市地表上相互影响密不可分的统一的物理过程。其次，在RCSO模式下，只有当管网检查井出现溢流后，才驱动地表淹没模型进行地表洪水淹没计算，因此洪水淹没仅可能发生在检查井周围区域，这与城市暴雨洪水发生的实际过程也是不符的。实际上，城市降雨积水，除了管网检查井溢流造成淹没外，更多的情况是，由于地表地形坡度较大或者管网排水能力不足等，导致雨水快速向低洼地汇集，甚至受管网下游顶托，雨水尚未进入管网前就积滞在低洼区域，从而形成局地洪水淹没。
2017年《水力发电学报》杂志的第36卷第5期公开了曾照洋、、王兆礼、吴旭树、赖成光、陈晓宏的文章《基于SWMM 和LISFLOOD 模型的暴雨内涝模拟研究》。该文章将二维地表模型（LISFLOOD-FP）与一维管网模（SWMM模型）进行耦合，对东莞市典型区域进行暴雨内涝模拟。该方法是典型的RCSO耦合模式，管网模型和地表模型之间是基于检查井的单向耦合，其缺点就是上所述RCSO模式共同的缺点，即将降雨产汇流和地表淹没过程分隔，地表淹没仅发生在管网检查井溢流处，因此模拟的地表淹没过程与实际情况不符。
2020年7月《Chinese Journal of Hydrodynamics》公开了赵仕霖、金生、杨宁的《地表水流与地下管流耦合的数值模拟研究》。该文章提出了二维地表模型与一维管网模型动态耦合的分区匹配算法，并通过试验算了进行了检验。该文章中，地表与排水管网之间的耦合是通过检查井进行流量交换来完成的，并假设检查井处于打开状态。该文章展示的地表与管网的耦合方法，不足之处在于，它没有区分实际排水管网系统中检查井与雨水篦子的不同作用，实际上城市排水系统中的检查井都被盖板覆盖的，并非出于打开状态，地表与排水管网之间的水量交换主要是通过雨水篦子完成的，雨水篦子和检查井须分别处理。另外，该文章所展示的地表与管网的耦合方法不适用于管网资料不足的地区。
发明内容
为了解决上述的技术问题，本发明提出的一种地表与管网分布式的直接双向耦合方法，提出了一种新的直接双向耦合模式，即降雨（Rainfall）—地表（Overland）—管网（Sewer）的耦合模式，在该模式下，地表采用分布式的水文水动力学模型，降雨直接作用在地表栅格上，降雨产汇流和积水淹没均基于统一的地表栅格进行计算，同时地表模型与管网模型进行双向耦合。
本发明的目的是提供一种地表与管网分布式的直接双向耦合方法，包括获取地表水文水动力数据，还包括以下步骤：
步骤1：依据所述地表水文水动力数据生成地表模型和管网模型，所述地表模型为地表分布式水文水动力模型；
步骤2：在所述地表模型和所述管网模型的基础上，根据城市排水管网实际组成以及数据覆盖情况，构建地表模型与管网模型之间的耦合模式，所述耦合模式至少包括物理耦合和概念耦合。
优选的是，所述地表水文水动力数据包括地表水深、地表流速、降雨强度、下渗强度、排水强度、摩阻、管道过流面积、管道内流量、管道内流速、水头和水力半径。
在上述任一方案中优选的是，基于所述地表分布式水文水动力模型，降雨产汇流和地表积水均基于同一地表栅格进行计算，计算公式为：
Figure 113170DEST_PATH_IMAGE001
Figure 81126DEST_PATH_IMAGE002
Figure 990307DEST_PATH_IMAGE003
其中，
Figure 701911DEST_PATH_IMAGE004
是偏微分符号，t是时间，h是水深，u和v是x和y方向的流速，z是水位， r、f和c分别为降雨强度、下渗强度和排水强度，g为重力加速度，n是糙率系数，
Figure 629416DEST_PATH_IMAGE005
为速度的模，且
Figure 768273DEST_PATH_IMAGE006
。
在上述任一方案中优选的是，所述管网模型为一维水动力模型，计算公式为：
Figure 161821DEST_PATH_IMAGE007
Figure 411537DEST_PATH_IMAGE008
其中，A是管道过流面积，Q是流量，U是流速，Z是水头，且
Figure 396810DEST_PATH_IMAGE009
，R是水力半径，H为水深。
在上述任一方案中优选的是，所述物理耦合表示地表和管网之间根据雨水汇集的实际物理过程进行耦合。
在上述任一方案中优选的是，所述物理耦合适用于排水管网资料完备较好的区域。
在上述任一方案中优选的是，所述物理耦合的步骤包括以下子步骤：
步骤201：将代表雨水篦子和检查井的管网节点与它所在的地表网格单元之间进行一对一的空间位置耦合；
步骤202：根据城市雨水排水系统中雨水篦子和检查井所承担的实际功能，对雨水篦子和检查井进行区别计算。
在上述任一方案中优选的是，所述雨水篦子的计算方法包括水流方向和大小采用具有物理意义的堰流公式计算，
Figure 768886DEST_PATH_IMAGE010
Figure 839610DEST_PATH_IMAGE011
其中，Q s 为流量，Z g 为地表水位，Z s 为管网节点水头，m为堰流系数，
Figure 706066DEST_PATH_IMAGE012
为淹没系数，L 1 为节点孔口周长，H 1 和H 2为地表与节点水头差，A 1 为节点孔口面积。
在上述任一方案中优选的是，所述检查井的计算方法包括只有当管网水头足够大顶开盖板，管网水流从检查井溢流进入地表，计算公式为：
Figure 545846DEST_PATH_IMAGE013
。
在上述任一方案中优选的是，所述概念耦合是指地表网格单元与管网节点之间基于某一种经验关系而非物理过程进行耦合。
在上述任一方案中优选的是所述概念耦合适用于排水管网资料缺失的区域。
在上述任一方案中优选的是，所述概念耦合的步骤包括以下子步骤：
步骤211：根据实地调研或者空间就近关系，确定每个区域的排水出口；
步骤212：根据排水出口节点的位置，指定地表网格单元对应的出口节点；
步骤213：根据该区域的设计排水能力，指定地表网格的排水能力C i ，同时考虑地表网格单元的水位和管网节点的水位，计算每个网格的实际排水速度c i 。
在上述任一方案中优选的是，所述每个网格的实际排水速度c i 的计算公式为：
Figure 88823DEST_PATH_IMAGE014
其中，h i 为第i个地表网格单元的水深；C i 为第i个地表网格单元的排水能力，
Figure 381264DEST_PATH_IMAGE015
为时间步长。
在上述任一方案中优选的是，每个概念耦合的节点接收概念排水的速度等于与其概念耦合的所有地表单元的排水速度和Q s ，计算公式为：
Figure 51411DEST_PATH_IMAGE016
其中，A i 为第i个网格面积，c i 为第i个网格的实际排水速度。
本发明提出了一种地表与管网分布式的直接双向耦合方法，在实际城市暴雨洪水模拟过程中，将物理耦合与概念耦合相结合，对于管网资料详实的区域，采用物理耦合，对于管网资料不足的区域，采用概念耦合，可以更好地模拟地表洪水发生情况，评估地表洪水淹没风险。
附图说明
图1为按照本发明的地表与管网分布式的直接双向耦合方法的一优选实施例的流程图。
图2为按照本发明的地表与管网分布式的直接双向耦合方法的地表网格单元与管网节点之间的物理耦合的一实施例的示意图。
图3为按照本发明的地表与管网分布式的直接双向耦合方法的地表网格单元与管网节点之间的概念耦合的一实施例的示意图。
图4为按照本发明的地表与管网分布式的直接双向耦合方法的某城市区域地表模型与管网模型的另一实施例的示意图。
图5为按照本发明的地表与管网分布式的直接双向耦合方法的概念耦合的另一实施例的示意图。
图6为按照本发明的地表与管网分布式的直接双向耦合方法的基于分布式耦合的一实施例的洪水淹没模拟效果示意图。
具体实施方式
下面结合附图和具体的实施例对本发明做进一步的阐述。
实施例一
如图1所示，执行步骤100，获取地表水文水动力数据。所述地表水文水动力数据包括地表水深、地表流速、降雨强度、下渗强度、排水强度、摩阻、管道过流面积、管道内流量、管道内流速、水头和水力半径。
执行步骤110，依据所述地表水文水动力数据生成地表模型和管网模型。所述地表模型为地表分布式水文水动力模型。基于所述地表分布式水文水动力模型，降雨产汇流和地表积水均基于同一地表栅格进行计算，计算公式为：
Figure 745697DEST_PATH_IMAGE017
Figure 459575DEST_PATH_IMAGE018
Figure 317941DEST_PATH_IMAGE019
其中，
Figure 713150DEST_PATH_IMAGE004
是偏微分符号，t是时间，h是水深，u和v是x和y方向的流速，z是水位， r、f和c分别为降雨强度、下渗强度和排水强度，g为重力加速度，n是糙率系数，
Figure 261943DEST_PATH_IMAGE005
为速度的模，且
Figure 412302DEST_PATH_IMAGE006
。
所述管网模型为一维水动力模型，计算公式为：
Figure 679335DEST_PATH_IMAGE007
Figure 711792DEST_PATH_IMAGE008
其中，A是管道过流面积，Q是流量，U是流速，Z是水头，且
Figure 115092DEST_PATH_IMAGE009
，R是水力半径，H为水深，n为糙率。
执行步骤120，在所述地表模型和所述管网模型的基础上，根据城市排水管网实际组成以及数据覆盖情况，构建地表模型与管网模型之间的耦合模式，所述耦合模式至少包括物理耦合和概念耦合。
执行步骤130，判断管网数据是否完整。如果管网数据完整，则执行步骤140，构建地表模型与管网模型之间的物理耦合。所述物理耦合表示地表和管网之间根据雨水汇集的实际物理过程进行耦合。所述物理耦合适用于排水管网资料完备的区域。在步骤140中，执行步骤141，将代表雨水篦子和检查井的管网节点与它所在的地表网格单元之间进行一对一的空间位置耦合；执行步骤142，根据城市雨水排水系统中雨水篦子和检查井所承担的实际功能，对雨水篦子和检查井进行区别计算。所述雨水篦子的计算方法包括水流方向和大小采用具有物理意义的堰流公式计算，
Figure 249401DEST_PATH_IMAGE010
Figure 3731DEST_PATH_IMAGE011
其中，Q s 为流量，Z g 为地表水位，Z s 为管网节点水头，m为堰流系数，
Figure 803059DEST_PATH_IMAGE012
为淹没系数，L 1 为节点孔口周长，H 1 和H 2为地表与节点水头差，A 1 为节点孔口面积。
所述检查井的计算方法包括只有当管网水头足够大顶开盖板，管网水流可以从检查井溢流进入地表，计算公式为：
Figure 326445DEST_PATH_IMAGE013
。
如果管网数据不完整，则执行步骤150，构建地表模型与管网模型之间的概念耦合。所述概念耦合是指地表网格单元与管网节点之间基于某一种经验关系而非物理过程进行耦合。所述概念耦合适用于排水管网资料缺失的区域。在步骤150中，执行步骤151，根据实地调研或者空间就近关系，确定每个区域的排水出口；执行步骤152，根据排水出口节点的位置，指定地表网格单元对应的出口节点；执行步骤153，根据该区域的设计排水能力，指定地表网格的排水能力C i ，同时考虑地表网格单元的水位和管网节点的水位，计算每个网格的实际排水速度c i 。所述每个网格的实际排水速度c i 的计算公式为：
Figure 569338DEST_PATH_IMAGE014
其中，h i 为第i个地表网格单元的水深；C i 为第i个地表网格单元的排水能力，
Figure 545384DEST_PATH_IMAGE015
为时间步长。
每个概念耦合的节点接收概念排水的速度等于与其概念耦合的所有地表单元的排水速度和Q s ，计算公式为：
Figure 148404DEST_PATH_IMAGE016
其中，A i 为第i个网格面积，c i 为第i个网格的实际排水速度。
实施例二
本发明提出了一种新的直接双向耦合模式，即降雨（Rainfall）—地表（Overland）—管网（Sewer）的耦合模式，以下简称ROS模式。该模式下，地表采用分布式的水文水动力学模型，降雨直接作用在地表栅格上，降雨产汇流和积水淹没均基于统一的地表栅格进行计算，同时地表模型与管网模型进行双向耦合。该模式符合城市降雨造成积水的实际物理过程，从而可以更合理地模拟城市暴雨引起的地表洪水淹没过程。
在ROS模式中，首先地表模型采用的分布式的水文水动力模型，如公式（1）、（2）和（3）所示。这是地表模型与管网模型进行直接双向耦合的基础。基于地表分布式水文水动力模型，降雨产汇流和地表积水均基于同一地表栅格进行计算，从而避免了RCSO中降雨产汇流和地表洪水淹没人为分割的弊端。
Figure 791875DEST_PATH_IMAGE020
（1）
Figure 127041DEST_PATH_IMAGE021
（2）
Figure 669012DEST_PATH_IMAGE022
（3）
其中，
Figure 747827DEST_PATH_IMAGE004
是偏微分符号，t是时间，h是水深，u和v是x和y方向的流速，z是水位， r、f和c分别为降雨强度、下渗强度和排水强度，g为重力加速度，n是糙率系数，
Figure 308121DEST_PATH_IMAGE005
为速度的模，且
Figure 814189DEST_PATH_IMAGE006
。
管网模型采用一维水动力模型，如公式（4）、（5）所示。
Figure 30407DEST_PATH_IMAGE007
（4）
Figure 723032DEST_PATH_IMAGE008
（5）
其中，A是管道过流面积，Q是流量，U是流速，Z是水头，且
Figure 75515DEST_PATH_IMAGE009
，
Figure 80381DEST_PATH_IMAGE023
是水力半径。
本发明在上述地表分布式水文水动力模型和管网一维水动力模型的基础上，根据城市排水管网实际组成以及数据覆盖情况，构建地表模型与管网模型之间两种耦合模式，即物理耦合和概念耦合。
（1）物理耦合，即地表和管网之间根据雨水汇集的实际物理过程进行耦合。该耦合方式适用于管网数据质量较好的区域。一般来说，城市公共区域，如主干道路、广场等区域的地下排水管网，归市政管理部门管理，数值质量较好，且易进行现场勘测，因此，能够收集到雨水篦子和检查井等详细的管网数据，适用于构建物理耦合方式。在这种物理耦合模式下，首先将代表雨水篦子和检查井的管网节点与它所在的地表网格单元之间进行一对一的空间位置耦合，如图2所示。在完成空间位置耦合的基础上，对雨水篦子和检查井进行区别计算，这是因为在实际城市雨水排水系统中，雨水篦子和检查井承担着完全不同的作用。雨水篦子是雨水收集的终端，承担着直接收集地表雨水的功能，而检查井是不同管段的连接点且一般都有井盖覆盖，一般不承担直接收集地表雨水的功能。因此，在物理耦合模式下，对于雨水篦子，地表与管网之间可以双向交换水流，其水流方向和大小采用具有物理意义的堰流公式计算，如公式（6）和（7）所示。而对于检查井，由于其不直接收集地表的雨水，只有当管网水头足够大顶开盖板，管网水流可以从检查井溢流进入地表，如（7）所示。
Figure 518315DEST_PATH_IMAGE010
（6）
Figure 938932DEST_PATH_IMAGE011
（7）
（2）概念耦合，是指地表网格单元与管网节点之间基于某一种经验关系而非物理过程进行耦合。该耦合方式，适用于排水管网资料不足的区域，例如城市社区或者私密区域。这些区域一般较难获得内部详细的管网数据。首先根据实地调研或者空间就近关系，确定每个区域的排水出口，然后根据排水出口节点的位置，指定地表网格单元对应的出口节点，如图3所示，从而建立地表与管网之间概念上的耦合关系。根据该区域的设计排水能力，指定地表网格的排水能力
Figure 224551DEST_PATH_IMAGE024
，同时考虑地表网格单元的水位和管网节点的水位，计算每个网格的实际排水速度
Figure 72421DEST_PATH_IMAGE025
，如下公式（8）所示。
Figure 59969DEST_PATH_IMAGE014
（8）
每个概念耦合的节点，接收概念排水的速度，等于与其概念耦合的所有地表单元的排水速度和，如下公式（9）所示，
Figure 284277DEST_PATH_IMAGE016
（9）
通过上述描述可知，基于地表分布式的水文水动力模型与管网一维水动力模型，构建了地表模型与管网模型之间直接双向耦合模式。同时考虑到管网数据覆盖情况，提出了物理耦合和概念耦合两者耦合模式。其中物理耦合模式，考虑了城市排水管网中雨水篦子和检查井的不同作用，对雨水篦子和检查井采用不同的水流计算，可以更加合理地模拟地表与管网之间的水流交换，从而合理地模拟地表洪水淹没过程。物理耦合模式，符合城市降雨产汇流的实际物理过程，从而可以合理地模拟管网资料较好的区域的洪水淹没情况。而概念耦合方式，虽然基于一定的概念关系，但依然能够较好地评估管网不足地区地表发生积水的风险，是物理耦合模式的一种有效的补充。在实际城市暴雨洪水模拟过程中，将物理耦合与概念耦合相结合，对于管网资料详实的区域，采用物理耦合，对于管网资料不足的区域，采用概念耦合，可以更好地模拟地表洪水发生情况，评估地表洪水淹没风险。
实施例三
某一城市区域的地表模型和管网模型分布示意如图4所示，主干道路上的管网数据详实，雨水篦子和检查井分类且连接关系清晰可靠。而主干道路以外的社区内部管网数据缺失。针对现有的管网数据情况，将主干道路上的地表网格单元，与排水管网的雨水篦子和检查井，根据位置关系，构建一一对应的物理耦合关系。而对于居民区和建筑物的地表网格单元，根据空间就近关系，分配网格单元的排水出口，如图5所示，每个地表网格单元上所标识的数字，代表该单元网格的排水出口即检查井的编号。通过这种概念上的排水关系，居民区所在的区域，虽然缺乏实际的管网资料，但依然可以按照设计排水能力将该区域的降雨径流排泄到附近的管网里。通过上述物理耦合和概念耦合两者耦合方式相互补充，从而可以更加合理地表达出地表淹没的情况。如图6所示，模拟的淹没点S1和S3，位于管网资料详实的主干道路上，淹没点S1位于管网资料缺失的居民区。由此可见，本发明提出的地表和管网之间的直接双向耦合，不仅可以准确模拟管网资料覆盖区域的积水淹没，还可以合理评估管网资料不足地区的地表淹没风险，在城市暴雨洪水模拟应用中，具有较大的实用价值。
为了更好地理解本发明，以上结合本发明的具体实施例做了详细描述，但并非是对本发明的限制。凡是依据本发明的技术实质对以上实施例所做的任何简单修改，均仍属于本发明技术方案的范围。本说明书中每个实施例重点说明的都是与其它实施例的不同之处，各个实施例之间相同或相似的部分相互参见即可。对于系统实施例而言，由于其与方法实施例基本对应，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。

Claims (4)

1.一种地表与管网分布式的直接双向耦合方法，包括获取地表水文水动力数据，其特征在于，还包括以下步骤：
步骤1：依据所述地表水文水动力数据生成地表模型和管网模型，所述地表模型为地表分布式水文水动力模型；基于所述地表分布式水文水动力模型，降雨产汇流和地表积水均基于同一地表网格单元进行计算，计算公式为：
Figure FDA0003077371090000011
Figure FDA0003077371090000012
Figure FDA0003077371090000013
其中，
Figure FDA0003077371090000014
是偏微分符号，t是时间，h是地表水深，u和v是地表x和y方向的流速，z是地表水位，r、f和c分别为降雨强度、下渗强度和排水强度，g为重力加速度，n是糙率系数，|V|为速度的模，且
Figure FDA0003077371090000015
所述管网模型为一维水动力模型，计算公式为：
Figure FDA0003077371090000016
Figure FDA0003077371090000017
其中，A是管道过流面积，Q是管道流量，U是管道流速，Z是管道水头，且Z＝Z0+H，R是管道水力半径，H为管道水深；
步骤2：在所述地表模型和所述管网模型的基础上，根据城市排水管网实际组成以及数据覆盖情况，构建地表模型与管网模型之间的耦合模式，所述耦合模式至少包括物理耦合和概念耦合；
所述物理耦合表示地表和管网之间根据雨水汇集的实际物理过程进行耦合，所述物理耦合的步骤包括以下子步骤：
步骤201：将代表雨水篦子和检查井的管网节点与它所在的地表网格单元之间进行一对一的空间位置耦合；
步骤202：根据城市雨水排水系统中雨水篦子和检查井所承担的实际功能，对雨水篦子和检查井进行区别计算；
所述概念耦合是指地表网格单元与管网节点之间基于某一种经验关系而非物理过程进行耦合，所述概念耦合的步骤包括以下子步骤：
步骤211：根据实地调研或者空间就近关系，确定每个区域的排水出口；
步骤212：根据排水出口节点的位置，指定地表网格单元对应的出口节点；
步骤213：根据该区域的设计排水能力，指定地表网格单元的排水能力Ci，同时考虑地表网格单元的水位和管网节点的水位，计算每个地表网格单元的实际排水速度ci。
2.如权利要求1所述的地表与管网分布式的直接双向耦合方法，其特征在于，所述雨水篦子的计算方法包括水流方向和大小采用具有物理意义的堰流公式计算，
Figure FDA0003077371090000021
for zg＞zs
Figure FDA0003077371090000022
for zg≤zs
其中，Qs为流量，Zg为管网节点所在的地表网格单元的水位，Zs为管网节点的水头，m为堰流系数，σ为淹没系数，L1为管网节点的孔口周长，H1为管网节点所在的地表网格单元的水深，H2为管网节点的水头与其所在的地表网格单元的水位之间的差，即H2＝Zs-Zg，A1为管网节点的孔口面积。
3.如权利要求2所述的地表与管网分布式的直接双向耦合方法，其特征在于，所述检查井的计算方法包括只有当管网水头足够大顶开盖板，管网水流从检查井溢流进入地表，计算公式为：
Figure FDA0003077371090000031
4.如权利要求3所述的地表与管网分布式的直接双向耦合方法，其特征在于，所述每个地表网格单元的实际排水速度ci的计算公式为：
Figure FDA0003077371090000032
其中，hi为第i个地表网格单元的水深；Ci为第i个地表网格单元的排水能力，Δt为时间步长。
CN202110323521.0A 2021-03-26 2021-03-26 一种地表与管网分布式的直接双向耦合方法 Active CN112711825B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110323521.0A CN112711825B (zh) 2021-03-26 2021-03-26 一种地表与管网分布式的直接双向耦合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110323521.0A CN112711825B (zh) 2021-03-26 2021-03-26 一种地表与管网分布式的直接双向耦合方法

Publications (2)

Publication Number Publication Date
CN112711825A CN112711825A (zh) 2021-04-27
CN112711825B true CN112711825B (zh) 2021-07-30

Family

ID=75550278

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110323521.0A Active CN112711825B (zh) 2021-03-26 2021-03-26 一种地表与管网分布式的直接双向耦合方法

Country Status (1)

Country Link
CN (1) CN112711825B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113204927B (zh) * 2021-06-08 2021-11-19 中国水利水电科学研究院 一种针对城市降雨洪涝过程的不同尺度分级嵌套模拟方法
CN113344305B (zh) * 2021-08-05 2021-10-29 中国水利水电科学研究院 一种暴雨内涝事件的快速预测方法
CN113627817B (zh) * 2021-08-19 2024-05-28 东莞理工学院 一种基于城市内涝条件下道路连通性分析方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107832931A (zh) * 2017-10-31 2018-03-23 上海市政工程设计研究总院（集团）有限公司 一种平原水网地区内涝风险的模块化分析方法
CN108388714A (zh) * 2018-02-07 2018-08-10 杭州师范大学 流域水系和城市管网耦合的平原河网城市洪水模拟方法
CN109255185A (zh) * 2018-09-11 2019-01-22 中国水利水电科学研究院 一种基于城市地表地下管网的一二维水动力学耦合分析方法
CN110298076A (zh) * 2019-05-27 2019-10-01 广州奥格智能科技有限公司 一种基于gis和swmm的城市内涝智能建模及分析方法
CN111046551A (zh) * 2019-12-10 2020-04-21 河海大学 一种城市群排水过程模拟方法
CN111369059A (zh) * 2020-03-05 2020-07-03 山东大学 基于内涝快速模拟耦合模型的城市内涝预测方法及系统
CN112052545A (zh) * 2020-08-25 2020-12-08 水利部交通运输部国家能源局南京水利科学研究院 一种基于元胞自动机的城市地表径流与管网汇流耦合方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107832931A (zh) * 2017-10-31 2018-03-23 上海市政工程设计研究总院（集团）有限公司 一种平原水网地区内涝风险的模块化分析方法
CN108388714A (zh) * 2018-02-07 2018-08-10 杭州师范大学 流域水系和城市管网耦合的平原河网城市洪水模拟方法
CN109255185A (zh) * 2018-09-11 2019-01-22 中国水利水电科学研究院 一种基于城市地表地下管网的一二维水动力学耦合分析方法
CN110298076A (zh) * 2019-05-27 2019-10-01 广州奥格智能科技有限公司 一种基于gis和swmm的城市内涝智能建模及分析方法
CN111046551A (zh) * 2019-12-10 2020-04-21 河海大学 一种城市群排水过程模拟方法
CN111369059A (zh) * 2020-03-05 2020-07-03 山东大学 基于内涝快速模拟耦合模型的城市内涝预测方法及系统
CN112052545A (zh) * 2020-08-25 2020-12-08 水利部交通运输部国家能源局南京水利科学研究院 一种基于元胞自动机的城市地表径流与管网汇流耦合方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"Research and development of flood control";Chai Fuxin 等;《EPiC Series in Engineering》;20181231;第3卷;第731-738页 *

Also Published As

Publication number Publication date
CN112711825A (zh) 2021-04-27

Similar Documents

Publication Publication Date Title
CN112711825B (zh) 一种地表与管网分布式的直接双向耦合方法
CN109492299B (zh) 基于swmm与modflow耦合的水资源模拟方法
CN107832931B (zh) 一种平原水网地区内涝风险的模块化分析方法
CN109657841B (zh) 一种城市暴雨内涝积水深度提取方法
Aronica et al. Drainage efficiency in urban areas: a case study
CN114240119A (zh) 基于数字孪生的国土空间全域全要素防洪防涝系统及预警方法
CN115391712A (zh) 一种城市洪涝风险预测方法
CN106372328A (zh) 一种基于排水能力评估的排水管道修复方法及系统
CN114647881B (zh) 顾及建筑物微观水文过程的城市内涝建模方法
CN111581767B (zh) 一种管网-河流耦合模型校核特征参数率定方法
CN111475950B (zh) 一种下凹式立交桥下雨洪的模拟方法
CN106779232A (zh) 一种城市内涝的建模预测方法
CN110222427A (zh) 一种基于数学模型城市内涝的分析方法
CN110532641B (zh) 一种地表网格分层建模方法及系统
CN106529739A (zh) 基于网格化的洪水预测方法
Shahapure et al. Coastal urban flood simulation using FEM, GIS and remote sensing
CN112084671A (zh) 城市时变增益降雨-径流过程模拟计算方法
CN110187413A (zh) 一种城市内涝预报方法、电子设备及存储介质
Filipova Urban flooding in Gothenburg-A MIKE 21 study
CN110287595B (zh) 一种城市不同下垫面减灾效应分析方法
CN115510771A (zh) 一种基于水文水动力耦合模型的社区雨涝数值模拟方法
Zeleňáková et al. Impact of land use changes on surface runoff in urban areas-Case study of Myslavsky Creek Basin in Slovakia.
Boonya-aroonnet et al. Modeling of urban flooding in Bangkok
Méndez Hydraulic analysis of urban drainage systems with conventional solutions and sustainable technologies: Case study in Quito, Ecuador
CN113887053A (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