CN109543275B - 一种城区地表径流二维数值模拟方法 - Google Patents
一种城区地表径流二维数值模拟方法 Download PDFInfo
- Publication number
- CN109543275B CN109543275B CN201811359262.1A CN201811359262A CN109543275B CN 109543275 B CN109543275 B CN 109543275B CN 201811359262 A CN201811359262 A CN 201811359262A CN 109543275 B CN109543275 B CN 109543275B
- Authority
- CN
- China
- Prior art keywords
- city
- building
- rainwash
- confluence
- rain pipe
- 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
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- 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
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
- G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
-
- 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
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Civil Engineering (AREA)
- Architecture (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- Computer Graphics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了城区地表径流二维数值模拟方法。获取城区地形、地貌数据后,采用三角形非结构网格离散城区计算区域,建筑物所在区域不参与剖分,其轮廓线作为网格剖分的控制线,每根建筑物雨水管通过汇流连接边与地面相连。先计算每根雨水管的汇流面积,再计算每根雨水管的出流过程。二维网格单元设置初始计算条件后启动二维Godunov模型进行模拟计算,雨水管出流过程作为与连接边对应单元的流量边界条件,地面网格单元内的降雨在模型中当作源项处理。本发明充分考虑城区内建筑物汇流的特点,既能保证城区洪涝计算时整体的质量守恒特性,同时又能考虑雨水管出口水流动量对地面径流的影响,克服现有城市洪涝模型无法考虑建筑物汇流过程的缺点,提高计算精度。
Description
技术领域
本发明属于水利工程技术领域,特别涉及一种考虑建筑物汇流特点的城区地表径流二维数值模拟方法。
背景技术
近几年,城市暴雨洪涝灾害呈常态化现象,成为干扰城市正常运行、威胁城市公共安全的重大问题。在城市洪涝问题的研究工作中,城市洪涝模型成为解决城市规划、城市洪涝预警等领域的重要科学工具。
目前,城市洪涝过程的模拟大致来讲主要有三种模式:水文模型、简化模型和水动力模型。相比于水文模型和简化模型,基于物理机理的水动力模型在模拟具有复杂城市下垫面条件的水流运动时具有更好的准确性,能够提供更加翔实的水力要素信息。这类模型通常把计算区域离散成三角形(四边形)非结构网格或正方形网格,通过浅水方程或简化的动力波、扩散波模型来模拟地表的降雨径流过程。现有的城市洪涝水动力模型方面的分析,较多的是关注一个较大区域的城市洪水淹没趋势,因此在网格离散方面选取了相对较大的尺度,对城市地貌进行了较大的概化处理。然而,由于城市复杂下垫面条件下的水流运动极易受到微地形的影响,因此,精细化是城市洪涝模型发展的必然趋势。
城市地表分布有密集的建筑物,这点与普通的自然流域有显著不同,建筑物的汇流过程与普通城市地面汇流过程有明显不同,建筑物一般是通过雨水管将房顶的雨水收集后集中排到地面,然后再参与到城市地表的汇流运动中来。当前,尚无城市洪涝模型考虑城区建筑物与普通城市地面汇流过程的差异,势必对城市洪涝过程的模拟精度造成影响。
发明内容
本发明的目的是提供一种考虑建筑物汇流特点的城区地表径流二维数值模拟方法,该方法区别于传统的城市洪涝模型方法,详细考虑了建筑物汇流过程的特点,首先定义建筑物与普通地面的单元连接边,再计算每根雨水管对应的集水面积,然后计算每根雨水管出口的流量过程,该过程作为连接边对应三角形单元的入流过程进行处理。该方法在保证整体水量守恒特性的同时,还能够考虑雨水管出口水流动量的影响,从而达到提高城市洪涝模型计算精度的目的。
本发明的目的是通过以下方案实现的:
本发明为一种考虑建筑物汇流特点的城区地表径流二维数值模拟方法,在城市精细化地形、地貌数据的基础上,定义城市建筑物与普通地面的连接边,确定每根雨水管对应的建筑物顶部的集水面积,计算每根雨水管出口处的流量过程,该过程作为单元连接边的流量边界条件加入到计算中去。通过该方法能够很精细的将建筑物的汇流过程和普通城市地面的汇流过程区分和联系起来,提高城市洪涝模型的计算精度。该方法的具体步骤如下:
(1)获取城区详细的地形数据和楼房(含雨水管数据)、街道等地貌几何特征数据。
(2)计算区域离散:采用三角形非结构网格离散城区计算区域,建筑物所在区域不参与剖分,其轮廓线作为网格剖分的控制线。剖分完成后,建筑物轮廓线被离散作为多个新生成的三角形单元的单元边,建筑物雨水管在水平面的投影点落在哪条单元边上,哪条边则被定义为连接建筑物和普通城市地面的汇流连接边。
(3)计算每根雨水管的楼顶集水面积。
(4)计算每条汇流连接边的流量过程:建筑物顶部的降雨径流通过雨水管对应的汇流连接边进入网格单元,每根雨水管末端的出流流量过程如式(1)所示:
Qk(t)=I(t-τ)*Ak (1)
其中,Qk(t)为第k根雨水管出口流量过程;t为时间变量;τ为建筑物顶部到地面的平均汇流时间,表述为τ=τ1+τ2,τ1为雨水管对应集水区的平面汇流时间,τ2为雨水管垂向的汇流时间;I为降雨强度;Ak为第k根雨水管的集水面积。
(5)设置计算初始条件:采用地形数据对三角形网格节点和网格型心进行高程插值,获取网格高程数据;给每个网格设置糙率和初始条件,初始流速和水深均设置为0。
(6)城区二维地表径流计算:采用基于完整二维浅水方程组的Godunov模型模拟城区二维地表径流过程,各条汇流连接边处的Qk(t)作为各对应单元的流量边界条件处理,这样既能保证建筑物所在区域的径流质量守恒特性,又能考虑雨水管末端水流自身的动量影响。普通地面网格单元内的降雨作为源项处理,地表汇流过程直接采用二维Godunov模型进行模拟。
进一步优选的方案为:步骤(1)中获取城区地形数据和楼房、街道地貌几何特征数据,比例尺为1:2000或1:2000以上。
进一步的,步骤(2)计算区域离散中,为保证数值计算精度,计算区域建筑物轮廓线的剖分尺度不宜超过5m。
进一步的,步骤(3)中,如果有详细的建筑物顶部高程数据,则可以根据D8算法计算每根雨水管的集水面积;如果无楼顶高程数据,可将建筑物顶部的集水面积平均分配到每根雨水管。
进一步的,步骤(6)城市二维地表径流计算中,所采用的模型为基于完整二维浅水方程组的Godunov模型,该模型所采用的控制方程的守恒形式如式(2)所示:
其中:
h为水深,u,v分别为x,y方向的流速,t为时间,分别为x,y方向的单位向量;
分别为x,y方向的坡度,Zb为地面高程,g为重力加速度;
分别为x,y方向的摩阻项,其中n为Manning糙率系数,r为净雨源项,由降雨数据经过产流计算得到;
在模型中,采用Roe格式计算界面通量,底坡项直接积分处理,摩阻项采用半隐式的处理方法,该模型能够捕捉城市地表径流运动的细节,如十字路口水流交互发生的流态变化以及坡度较大道路处发生的滚波等。
进一步的,步骤(6)城市二维地表径流计算中,在处理连接边处的流量过程Qk(t)时,采用Roe格式的近似Riemann解直接求解通过汇流连接边处的数值通量,公式如式(3)所示:
其中,F*·n为通过汇流连接边外法线方向的数值通量;h*,u*,v*分别为该连接边中点处的水深、x方向流速和y方向流速;nx,ny分别为汇流连接边外法线单位向量在x方向和y方向的分量;g为重力加速度。
本发明的有益效果:
在本发明中,充分考虑城区建筑物汇流过程与城市普通地面汇流过程的不同特点,在网格剖分时即将建筑物与普通地面区分开来,建筑物所在区域不参与网格剖分,将雨水管所对应的网格单元边定义为连接边(连接建筑物与普通地面),雨水管出口的流量过程作为连接边对应单元的流量边界处理,普通网格内的降雨作为源项处理,网格内的水流运动采用二维Godunov模型进行模拟,这种做法可以充分保证城区洪涝过程计算时整体的水量质量守恒特性,同时又能很好的考虑雨水管出口水流动量对地面径流的影响。在本发明中,有区别的考虑建筑物汇流和城市普通地面汇流过程的不同,可以更好的模拟城区精细化的地表汇流过程,提高现有城市洪涝过程模拟的精度。
附图说明
图1为本发明的城区地表径流二维数值模拟方法流程图;
图2为本发明的城区地表径流数二维值模拟方法示意图。
具体实施方式
下面结合附图1和附图2对本发明作进一步详细说明。
本发明为一种考虑建筑物汇流特点的城区地表径流二维数值模拟方法,充分考虑城区内建筑物汇流与普通地面汇流的区别与联系(即建筑物顶部的降雨径流过程一般是通过雨水管道流到建筑物周围的地面上后,再参与到城区地面汇流的过程中来)来进行更精确的城区地表径流数值过程模拟。该方法主要包括以下步骤:
(1)获取城区详细的地形数据和楼房(含雨水管数据)、街道等地貌几何特征数据,比例尺1:2000以上(含1:2000)。
(2)计算区域离散:采用三角形非结构网格离散城区计算区域,建筑物所在区域不参与剖分,其轮廓线作为网格剖分的控制边界线。剖分完成后,建筑物轮廓线被离散作为多个新生成的三角形单元的单元边,建筑物雨水管在水平面的投影点落在哪条单元边上,哪条边则被定义为连接建筑物和普通城市地面的汇流连接边。为保证数值计算精度,建筑物计算区域轮廓线的剖分尺度不宜超过5m。
(3)计算每根雨水管的楼顶集水面积:如果有详细的建筑物顶部高程数据,则可以根据D8算法计算每根雨水管的集水面积;如果无楼顶高程数据,可将建筑物顶部的集水面积平均分配到每根雨水管。
(4)计算每条汇流连接边的流量过程:建筑物顶部的降雨径流通过雨水管对应的汇流连接边进入网格单元,每根雨水管末端的出流流量过程如式(1)所示:
Qk(t)=I(t-τ)*Ak (1)
其中,Qk(t)为第k根雨水管出口流量过程;t为时间变量;τ为建筑物顶部到地面的平均汇流时间,表述为τ=τ1+τ2,τ1为雨水管对应集水区的平面汇流时间,τ2为雨水管垂向的汇流时间;I为降雨强度;Ak为第k根雨水管的集水面积。
(5)设置计算初始条件:采用地形数据对三角形网格节点和网格型心进行高程插值,获取网格高程数据;给每个网格设置糙率和初始条件,初始流速和水深均设置为0。
(6)城区二维地表径流计算:采用基于完整二维浅水方程组的Godunov模型模拟城区二维地表径流过程,各条汇流连接边处的Qk(t)作为各对应单元的流量边界条件处理,这样既能保证建筑物所在区域的径流质量守恒特性,又能考虑雨水管末端水流自身的动量影响。普通地面网格单元内的降雨作为源项处理,地表汇流过程直接采用二维Godunov模型进行模拟。
基于完整二维浅水方程组的Godunov模型,该模型所采用的控制方程的守恒形式如式(2)所示:
其中:
h为水深,u,v分别为x,y方向的流速,t为时间,分别为x,y方向的单位向量;
分别为x,y方向的坡度,Zb为地面高程,g为重力加速度;
分别为x,y方向的摩阻项,其中n为Manning糙率系数,r为净雨源项,由降雨数据经过产流计算得到;
在模型中,采用Roe格式计算界面通量,底坡项直接积分处理,摩阻项采用半隐式的处理方法,该模型能够捕捉城市地表径流运动的细节,如十字路口水流交互发生的流态变化以及坡度较大道路处发生的滚波等。具体的模型构建过程可以参见如下文献(张大伟,权锦,马建明,向立云.基于Godunov格式的流域地表径流二维数值模拟.水利学报,2018,49(7):787-794.)
在处理连接边处的流量过程Qk(t)时,采用Roe格式的近似Riemann解直接求解通过汇流连接边处的数值通量,公式如式(3)所示:
其中,F*·n为通过汇流连接边外法线方向的数值通量;h*,u*,v*分别为该连接边中点处的水深、x方向流速和y方向流速;nx,ny分别为汇流连接边外法线单位向量在x方向和y方向的分量;g为重力加速度。
上述的实施例仅是本发明的部分体现,并不能涵盖本发明的全部,在上述实施例以及附图的基础上,本领域技术人员在不付出创造性劳动的前提下可获得更多的实施方式,因此这些不付出创造性劳动的前提下获得的实施方式均应包含在本发明的保护范围内。
Claims (6)
1.一种城区地表径流二维数值模拟方法,其特征在于:
包括以下步骤:
(1)获取城区地形数据和楼房、街道地貌几何特征数据,楼房地貌几何特征数据包含雨水管数据;
(2)计算区域离散:采用三角形非结构网格离散城区计算区域,建筑物所在区域不参与剖分,建筑物轮廓线作为网格剖分的控制边界线;剖分完成后,建筑物轮廓线被离散作为多个新生成的三角形单元的单元边,建筑物雨水管在水平面的投影点落在的单元边,被定义为连接建筑物和普通城市地面的汇流连接边;
(3)计算每根雨水管的楼顶集水面积;
(4)计算每条汇流连接边的流量过流:建筑物顶部的降雨径流通过雨水管对应的汇流连接边进入网格单元,每根雨水管末端的出流流量过程如式(1)所示:
Qk(t)=I(t-τ)*Ak (1)
其中,Qk(t)为第k根雨水管出口流量过程;t为时间变量;τ为建筑物顶部到地面的平均汇流时间,表述为τ=τ1+τ2,τ1为雨水管对应集水区的平面汇流时间,τ2为雨水管垂向的汇流时间;I为降雨强度;Ak为第k根雨水管对应的集水面积;
(5)设置计算初始条件:采用地形数据对三角形网格节点和网格型心进行高程插值,获取网格高程数据;给每个网格设置糙率和初始条件,初始流速和水深均设置为0;
(6)城区二维地表径流计算:采用基于完整二维浅水方程组的Godunov模型模拟城区二维地表径流过程,各条汇流连接边处的Qk(t)作为各对应单元流量边界条件处理;普通地面网格单元内的降雨作为源项处理,地表汇流过程直接采用二维Godunov模型进行模拟。
2.根据权利要求1所述的城区地表径流二维数值模拟方法,其特征在于:
步骤(1)中获取城区地形数据和楼房、街道地貌几何特征数据,比例尺为1:2000或1:2000以上。
3.根据权利要求1所述的城区地表径流二维数值模拟方法,其特征在于:步骤(2)计算区域离散中,计算区域建筑物轮廓线的剖分尺度不超过5m。
4.根据权利要求1所述的城区地表径流二维数值模拟方法,其特征在于:如果有建筑物顶部高程数据,则根据D8算法计算每根雨水管的集水面积;如果无楼顶高程数据,则将建筑物顶部的集水面积平均分配到每根雨水管。
5.根据权利要求1所述的城区地表径流二维数值模拟方法,其特征在于:步骤(6)城市二维地表径流计算中,所采用的模型为基于完整二维浅水方程组的Godunov模型,该模型所采用的控制方程组的守恒形式如式(2)所示:
其中:
h为水深,u,v分别为x,y方向的流速,t为时间,分别为x,y方向的单位向量;
分别为x,y方向的坡度,Zb为地面高程,g为重力加速度;
分别为x,y方向的摩阻项,其中n为Manning糙率系数,r为净雨源项,由降雨数据经过产流计算得到。
6.根据权利要求1所述的城区地表径流二维数值模拟方法,其特征在于:步骤(6)城市二维地表径流计算中,在处理连接边处的流量过程Qk(t)时,直接采用Roe格式的近似Riemann解求解通过汇流连接边处的数值通量,公式如式(3)所示:
其中,F*·n为通过汇流连接边法线方向的数值通量;h*,u*,v*分别为该连接边中点处的水深、x方向流速和y方向流速;nx,ny分别为汇流连接边外法线单位向量在x方向和y方向的分量;g为重力加速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811359262.1A CN109543275B (zh) | 2018-11-15 | 2018-11-15 | 一种城区地表径流二维数值模拟方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811359262.1A CN109543275B (zh) | 2018-11-15 | 2018-11-15 | 一种城区地表径流二维数值模拟方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109543275A CN109543275A (zh) | 2019-03-29 |
CN109543275B true CN109543275B (zh) | 2019-07-09 |
Family
ID=65847418
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811359262.1A Active CN109543275B (zh) | 2018-11-15 | 2018-11-15 | 一种城区地表径流二维数值模拟方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109543275B (zh) |
Families Citing this family (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110188483B (zh) * | 2019-06-03 | 2020-10-09 | 中国水利水电科学研究院 | 一种二维水动力水质模型构建方法 |
CN110362925B (zh) * | 2019-07-16 | 2020-05-19 | 中国水利水电科学研究院 | 一种包含库区的土石坝漫顶溃决洪水数值模拟方法 |
CN110378023B (zh) * | 2019-07-19 | 2023-07-11 | 同济大学建筑设计研究院(集团)有限公司 | 屋面降雨径流模拟方法、装置、计算机设备和存储介质 |
CN110765694B (zh) * | 2019-11-21 | 2021-11-05 | 华南理工大学 | 一种基于简化型浅水方程组的城市地表水流数值模拟方法 |
CN110929443B (zh) * | 2019-12-02 | 2021-11-30 | 中国水利水电科学研究院 | 一种基于高精度地形概化的二维洪水模拟方法 |
CN110992653A (zh) * | 2019-12-13 | 2020-04-10 | 软通动力信息技术有限公司 | 城市内涝预警***及方法 |
CN111046567B (zh) * | 2019-12-18 | 2020-07-31 | 中国水利水电科学研究院 | 一种基于Godunov格式的城市排水管网水流数值模拟方法 |
CN111931263B (zh) * | 2020-05-19 | 2021-09-14 | 长江水利委员会水文局 | 一种基于优化求解平面二维浅水方程的恒定流模拟方法 |
CN111931330B (zh) * | 2020-06-19 | 2024-01-05 | 西安理工大学 | 一种无管网资料城区的管网排水过程概化计算方法 |
CN112906316B (zh) * | 2021-02-05 | 2022-02-01 | 浙江省水利河口研究院(浙江省海洋规划设计研究院) | 一种提升平原防洪排涝能力的方法 |
CN113204927B (zh) * | 2021-06-08 | 2021-11-19 | 中国水利水电科学研究院 | 一种针对城市降雨洪涝过程的不同尺度分级嵌套模拟方法 |
CN117150823B (zh) * | 2023-10-31 | 2024-02-02 | 中国海洋大学 | 一种滨海城市洪水灾害的高分辨率快速模拟方法 |
CN118052059B (zh) * | 2024-02-27 | 2024-07-12 | 中国水利水电科学研究院 | 城市降雨径流计算中不透水构筑物雨水收集分配方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106599457A (zh) * | 2016-12-13 | 2017-04-26 | 中国水利水电科学研究院 | 一种基于Godunov格式一、二维耦合技术的山洪数值模拟方法 |
CN107239657A (zh) * | 2017-05-31 | 2017-10-10 | 中国水利水电科学研究院 | 一种面向对象的水动力学建模要素管理方法 |
CN107451372A (zh) * | 2017-08-09 | 2017-12-08 | 中国水利水电科学研究院 | 一种运动波与动力波相结合的山区洪水过程数值模拟方法 |
CN108021780A (zh) * | 2018-01-24 | 2018-05-11 | 国电南瑞科技股份有限公司 | 一种基于无规则非结构网格模型的山洪动态仿真方法 |
CN108182543A (zh) * | 2018-01-17 | 2018-06-19 | 福建四创软件有限公司 | 一种精细化网格内涝水淹预报方法 |
CN108446502A (zh) * | 2018-03-22 | 2018-08-24 | 中国水利水电科学研究院 | 一种利用完整二维浅水方程组获得流域单位线的方法 |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9726782B2 (en) * | 2012-11-30 | 2017-08-08 | International Business Machines Corporation | Methods, systems and computer program storage devices for generating a response to flooding |
CN103886151B (zh) * | 2014-03-27 | 2017-02-08 | 北京工业大学 | 不同暴雨情景下城市下凹式立交桥区简化水力模型构建方法 |
CN107832931B (zh) * | 2017-10-31 | 2022-03-29 | 上海市政工程设计研究总院(集团)有限公司 | 一种平原水网地区内涝风险的模块化分析方法 |
CN108491604A (zh) * | 2018-03-13 | 2018-09-04 | 广州地理研究所 | 一种亚热带地区水土流失耦合模型构建方法 |
CN108793609A (zh) * | 2018-07-03 | 2018-11-13 | 四川理工学院 | 城市雨水管理与处理的方法 |
-
2018
- 2018-11-15 CN CN201811359262.1A patent/CN109543275B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106599457A (zh) * | 2016-12-13 | 2017-04-26 | 中国水利水电科学研究院 | 一种基于Godunov格式一、二维耦合技术的山洪数值模拟方法 |
CN107239657A (zh) * | 2017-05-31 | 2017-10-10 | 中国水利水电科学研究院 | 一种面向对象的水动力学建模要素管理方法 |
CN107451372A (zh) * | 2017-08-09 | 2017-12-08 | 中国水利水电科学研究院 | 一种运动波与动力波相结合的山区洪水过程数值模拟方法 |
CN108182543A (zh) * | 2018-01-17 | 2018-06-19 | 福建四创软件有限公司 | 一种精细化网格内涝水淹预报方法 |
CN108021780A (zh) * | 2018-01-24 | 2018-05-11 | 国电南瑞科技股份有限公司 | 一种基于无规则非结构网格模型的山洪动态仿真方法 |
CN108446502A (zh) * | 2018-03-22 | 2018-08-24 | 中国水利水电科学研究院 | 一种利用完整二维浅水方程组获得流域单位线的方法 |
Non-Patent Citations (1)
Title |
---|
基于 Godunov格式的流域地表径流二维数值模拟;张大伟 等;《水利学报》;20180731;第49卷(第7期);第787-802页 |
Also Published As
Publication number | Publication date |
---|---|
CN109543275A (zh) | 2019-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109543275B (zh) | 一种城区地表径流二维数值模拟方法 | |
CN106884405B (zh) | 一种无资料地区溃堤型山洪灾害分析评价方法 | |
CN108022047A (zh) | 一种海绵城市水文计算方法 | |
CN110298076A (zh) | 一种基于gis和swmm的城市内涝智能建模及分析方法 | |
CN107563139A (zh) | 一种点源突发性水污染源事故溯源贡献程度计算方法 | |
CN109948866A (zh) | 一种基于swwm模拟的北方滨海地区海绵城市规划方法 | |
CN107832931A (zh) | 一种平原水网地区内涝风险的模块化分析方法 | |
CN109815305A (zh) | 一种无资料地区场次洪水径流过程反演的方法 | |
CN105894116A (zh) | 一种流域梯级水库与蓄滞洪区联合调度方法 | |
CN104933268B (zh) | 一种基于一维非恒定流数值模型的洪水分析方法 | |
CN113409550B (zh) | 基于径流汇流模拟的泥石流灾害预警方法及*** | |
Joseph et al. | Assessment of environmental flow requirements using a coupled surface water-groundwater model and a flow health tool: A case study of Son river in the Ganga basin | |
Ibrakhimov | Spatial and temporal dynamics of groundwater table and salinity in Khorezm (Aral Sea Basin), Uzbekistan | |
CN110232737A (zh) | 一种城市汇水区划分方法 | |
Formánek et al. | Two-dimensional model of Ciliwung river flood in DKI Jakarta for development of the regional flood index map | |
CN114942596B (zh) | 一种城市防洪排水智能控制*** | |
CN110187413A (zh) | 一种城市内涝预报方法、电子设备及存储介质 | |
CN112052545A (zh) | 一种基于元胞自动机的城市地表径流与管网汇流耦合方法 | |
CN106096129B (zh) | 一种基于山地汇水计算的山脚水面规模分析方法 | |
CN104750936A (zh) | 一种基于gis的城市水文径流计算方法及*** | |
CN110377868A (zh) | 一种基于实时雨情的动态水系提取方法 | |
CN113469456A (zh) | 一种河道水量预测方法 | |
Liu et al. | Urban flood modelling in Qiqihar city based on MIKE flood | |
CN105740812B (zh) | 一种基于数字表面模型的城市汇水区提取方法 | |
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 |