CN103530462A - 面向山洪演进数值模拟的计算网格流出率的修正方法 - Google Patents

面向山洪演进数值模拟的计算网格流出率的修正方法 Download PDF

Info

Publication number
CN103530462A
CN103530462A CN201310478469.1A CN201310478469A CN103530462A CN 103530462 A CN103530462 A CN 103530462A CN 201310478469 A CN201310478469 A CN 201310478469A CN 103530462 A CN103530462 A CN 103530462A
Authority
CN
China
Prior art keywords
grid
per unit
unit width
water
discharge per
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.)
Withdrawn
Application number
CN201310478469.1A
Other languages
English (en)
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.)
Nanjing Xiaozhuang University
Original Assignee
Nanjing Xiaozhuang 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 Nanjing Xiaozhuang University filed Critical Nanjing Xiaozhuang University
Priority to CN201310478469.1A priority Critical patent/CN103530462A/zh
Publication of CN103530462A publication Critical patent/CN103530462A/zh
Withdrawn legal-status Critical Current

Links

Images

Classifications

    • 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

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明涉及一种面向山洪演进数值模拟的计算网格流出率的修正方法，包括：通过修正系数对负水深网格进行修正，获得该网格的实际水深，克服了在以往的移动边界计算中，经常发生水面标高低于地面标高的情况，也就是计算后的水深变为负值，即出现负水深网格，从而导致计算域内的质量不能保证守恒，稳定性变差，甚至计算发散而得不到结果，直接影响到山洪演进数值模拟计算无法顺利进行的技术问题，通过引入所述修正系数λ，使当出现水面标高低于地面标高时，山洪演进数值模拟计算得已顺利进行。

Description

面向山洪演进数值模拟的计算网格流出率的修正方法
技术领域
本发明涉及一种面向山洪演进数值模拟的计算网格流出率的修正方法。
背景技术
山洪灾害是指山丘地区在强降雨影响下，短时间内形成具有较大洪峰流量的洪水。我国地处东亚季风区，山区和丘陵地区占国土地面积的三分之二。其中，山洪灾害的优先预防面积达97万 ，影响人口1.3亿。近年山洪灾害造成的死亡人数占全国洪涝灾害死亡人数的比例超过70%，成为造成人员伤亡的主要灾种。随着社会经济的发展，山洪灾害的防治工作越来越被重视。
早先，国际通用的山洪灾害的预测技术是对沟道、沟口进行实地采样，根据有可能的灾害种类和等级确定危险指数。最具代表性的是Aulitzky提出的荒溪分类及危险区制图指数法，通过收集9种指标51个具体因子划分出不同等级的危险区。随着地理信息系统、数字高程模型、遥感和卫星遥测等现代科学技术高速发展，基于平面浅水波方程的模拟方法被广泛应用于对流域内洪水、山洪及泥石流等灾害现象的预测和定量分析中。该方法不受模型实验相似性理论的限制，可快速、精确的揭示灾害发生的原因及过程，从而大大提高山洪等突发洪水的预见期。最常用的平面二维模拟数值方法包括有限元法、有限体积法和有限差分法。针对洪水模拟，Jin等提出了使用二阶迎风格式离散动量方程的非线性对流项和二阶leap-frog格式离散线性对流项模拟Nakdong 流域的泛滥过程。Roger 等提出一种MacCormack+TVD格式的边界拟合数值模型，通过检测每个时间步长下的水深是否到达干枯临界值判定动边界范围。丹麦DHI水资源与环境研究院采用隐式交替方向算法开发了水动力学模拟软件Mike21。
在面向山洪泥石流灾害二维数值模拟中，由于流域地形陡峭，水流量急速变化，在较大计算时间步长下的动边界处理过程中,对流出量简单的归零处理往往导致负水深，造成模拟过程中的质量与动量不守恒，最终导致计算数值不稳定甚至计算发散而得不到结果。近年来，Satofuka和Mizuyama（2005）采用一维数值模型计算了明渠山洪对山区河床、水坝造成的影响。而Nakatani等（2008）采用二维数值模拟模型，基于山洪淹没深度和沉降的变化对其形成的冲积扇进行了模拟计算。
由于山洪发生地的地形具有陡峭且凹凸不平的特点，因此洪水的边界范围、流速和水深都急剧变化。在以往的移动边界计算中，经常发生水面标高低于地面标高的情况，也就是计算后的水深变为负值。导致计算域内的质量不能保证守恒，稳定性变差，甚至计算发散而得不到结果，从而直接影响到山洪演进数值模拟计算无法顺利进行。
发明内容
本发明所要解决的技术问题是提供一种面向山洪演进数值模拟的计算网格流出率的修正方法，该修正方法解决了当出现水面标高低于地面标高的情况，山洪演进数值模拟计算无法顺利进行的技术问题。
 为了解决上述问题，本发明提供了一种面向山洪演进数值模拟的计算网格流出率的修正方法，包括：
①建立计算域内地形的数字高程模型，即，所述计算域按照一定空间步长（                                                
Figure 21577DEST_PATH_IMAGE001
、
Figure 81413DEST_PATH_IMAGE002
）进行分割后得到二维网格模型，定义（i，j）为该二维网格模型中一网格，其中，设定该二维网格模型中任一网格的水深为H，且与该网格按X轴同方向的单宽流量M，与Y轴同方向的单宽流量N；同时设定所述山洪演进数值模拟计算的时间步n，以及该时间步n对应的时间步长
Figure 66686DEST_PATH_IMAGE003
。
②通过单宽流量公式获得所述二维网格模型中各网格在n+1时间步的单宽流量M和单宽流量N，其方法，包括：
预设初始条件，即，在n时间步，且设定位于所述计算域的边界上的流量输入起始界的各边界网格的初值参数，该初值参数包括：初始水深
Figure 438762DEST_PATH_IMAGE004
，所述单宽流量M的对应流速矢量
Figure 322535DEST_PATH_IMAGE005
，以及所述单宽流量N的对应流速矢量
Figure 375942DEST_PATH_IMAGE006
；在所述单宽流量公式中位于所述计算域外的网格的初值参数与该网格相邻的一边界网格的初值参数相同，且位于所述计算域内的网格的初值参数的相应取值为0。
以及，所述单宽流量M和单宽流量N的初始值的计算公式：
Figure 278039DEST_PATH_IMAGE007
     （1）
Figure 758699DEST_PATH_IMAGE008
     （2）
所述单宽流量公式：
Figure 2013104784691100002DEST_PATH_IMAGE009
（3） 
Figure 926506DEST_PATH_IMAGE010
  （4） 
③根据所述各网格在n+1时间步的单宽流量M和单宽流量N，建立所述各网格的水深计算式，即
Figure 2013104784691100002DEST_PATH_IMAGE011
 （5）
其中， n+2时间步作为水深计算的起始时间步，H表示所述网格水深，表示在所述起始时间步时网格（i，j）的区域内的水深；
Figure 2013104784691100002DEST_PATH_IMAGE013
表示在n时间步的所述网格的区域内的水深，且该时间步的水深为预设值；
Figure 353256DEST_PATH_IMAGE014
表示在n+1时间步的所述网格在X轴方向上流体的流入单宽流量；表示相邻网格在X轴方向上的流入单宽流量，即，所述网格在X轴方向上的流出单宽流量； 
Figure 922571DEST_PATH_IMAGE016
表示在n+1时间步的所述网格在Y轴方向上流体的流入单宽流量；
Figure 317780DEST_PATH_IMAGE017
表示相邻网格在Y轴方向上的流入单宽流量，即，所述网格在Y轴方向上的流出单宽流量。
④若步骤③计算出一网格的水深H小于0，即该网格为负水深网格，则引入一修正系数
Figure 679622DEST_PATH_IMAGE018
重新对该负水深网格的水量进行修正计算，直到该负水深网格修正为实际水深后，再转入步骤②通过所述单宽流量公式获得所述二维网格模型中各网格在下一时间步的单宽流量M和单宽流量N。
其中，当一网格出现负水深时，根据负水深网格在当前时间步n的单宽流量M和单宽流量N计算出该负水深网格的流入总量和流出总量。
所述流入总量公式：     （6）
所述流出总量公式：
Figure 512263DEST_PATH_IMAGE022
   （7）
所述计算修正系数λ的公式：
        （8） 
并根据公式（1）、（2）、（3）、（4）、（6）、（7）修正当前负水深网格在n+1时间步的流入量
Figure 400902DEST_PATH_IMAGE024
和流出量
Figure 137913DEST_PATH_IMAGE025
，得出所述负水深网格在n+2时间步的实际水深：
Figure 723616DEST_PATH_IMAGE027
      （9）。
本发明相对于现有技术具有积极的效果：本发明克服了在以往的移动边界计算中，经常发生水面标高低于地面标高的情况，也就是计算后的水深变为负值，即出现负水深网格，从而导致计算域内的质量不能保证守恒，稳定性变差，甚至计算发散而得不到结果，直接影响到山洪演进数值模拟计算无法顺利进行的技术问题，通过引入所述修正系数λ，使当出现水面标高低于地面标高（负水深网格）时，山洪演进数值模拟计算得已顺利进行。
附图说明
为了清楚说明本发明的创新原理及其相比于现有产品的技术优势，下面借助于附图通过应用所述原理的非限制性实例说明一个可能的实施例。在图中：
图1为本发明的一种面向山洪演进数值模拟的计算网格流出率的修正方法的流程图；
图2为本发明中公式（6）和（7）的水流方向示意图。 
具体实施方式
参照附图，对本发明的实施方式和实施例进行详细说明。
图1示出本发明的面向山洪演进数值模拟的计算网格流出率的修正方法的流程图，现结合图1所述流程图对该修正方法的各步骤作具体说明。
在步骤S001中，建立计算域内地形的数字高程模型，即，所述计算域按照一定空间步长（
Figure 153460DEST_PATH_IMAGE001
、
Figure 939626DEST_PATH_IMAGE002
）进行分割后得到二维网格模型，定义（i，j）为该二维网格模型中一网格，其中，设定该二维网格模型中任一网格的水深为H，且与该网格按X轴同方向的单宽流量M，与Y轴同方向的单宽流量N；同时设定所述山洪演进数值模拟计算的时间步n，以及该时间步n对应的时间步长
Figure 542645DEST_PATH_IMAGE003
；其中，单宽流量M和单宽流量 N均为矢量，也可以用于表示相应流量的方向。
其中，所述数字高程模型，即 DEM，用一组有序数值阵列形式表示地面高程的一种实体地面模型，定义（i，j）为该二维网格模型中一网格，所述i和j的取值i=1,2,3,……，j=1,2,3,……，利用i和j来限定网格在所述二维网格模型中的具体位置，即网格（1,1），网格（1,2）类似的表示形式，所述网格为矩形网格，即可以是长方形，也可以是正方形，正方形网格也可以称为栅格DEM。建立计算域内地形的数字高程模型的方法，以及从数字高程模型提取出网格数据的方法在现有技术的文献中均有上述记载，这里不再重复。
所述单宽流量：单位宽度上河流或输水管的输水流量，这里的单位宽度即网格。
所述空间步长也可以简称步长，针对DEM空间的分辨力，也就是网格的精度，一般为30M、90M两种，就是用30*30或90*90的DEM网格来表示地形，当然也可以根据计算需要，另外设置相应步长，或者使用长方形网格（矩形网格）,其中
Figure 186116DEST_PATH_IMAGE001
、分别表示一个网格的长、宽距离。
所述时间步n为整个计算过程中的时刻间隔，即计算步，所述时间步n的取值n=1,2,3,……；所述时间步长：相邻两计算步的时间间隔，对应的时间步长，一般可以去0.1s或0.01s，也可以根据计算设置任意时间。
在步骤S002中，通过单宽流量公式获得所述二维网格模型中各网格在n+1时间步的M和N，其方法，包括：
预设初始条件，即，在n时间步，且设定位于所述计算域的边界上的流量输入起始界的各边界网格的初值参数，该初值参数包括：初始水深
Figure 17434DEST_PATH_IMAGE004
，所述单宽流量M的对应流速矢量
Figure 515412DEST_PATH_IMAGE005
，以及所述单宽流量N的对应流速矢量
Figure 83796DEST_PATH_IMAGE006
；在所述单宽流量公式中位于所述计算域外的网格的初值参数与该网格相邻的一边界网格的初值参数相同，且位于所述计算域内的网格的初值参数的相应取值为0。
以及，所述单宽流量M和单宽流量N的初始值的计算公式：
Figure 113063DEST_PATH_IMAGE028
     （1）
     （2）
所述单宽流量公式：
（3）
Figure 163037DEST_PATH_IMAGE031
 （4）
其中，所述初值参数可以通过多参数水文监测仪来获得，获得方式属于现有技术，这里不再详细叙述。
在步骤S003中，根据所述各网格在n+1时间步的单宽流量M和单宽流量N，建立所述各网格的水深计算式，即
Figure 600972DEST_PATH_IMAGE032
 （5）
其中， n+2时间步作为水深计算的起始时间步，H表示所述网格水深，表示在所述起始时间步时网格（i，j）的区域内的水深；
Figure 431842DEST_PATH_IMAGE013
表示在n时间步的所述网格的区域内的水深，且该时间步的水深为预设值；
Figure 279712DEST_PATH_IMAGE014
表示在n+1时间步的所述网格在X轴方向上流体的流入单宽流量；
Figure 17992DEST_PATH_IMAGE015
表示相邻网格在X轴方向上的流入单宽流量，即，所述网格在X轴方向上的流出单宽流量； 
Figure 304617DEST_PATH_IMAGE016
表示在n+1时间步的所述网格在Y轴方向上流体的流入单宽流量；
Figure 444742DEST_PATH_IMAGE017
表示相邻网格在Y轴方向上的流入单宽流量，即，所述网格在Y轴方向上的流出单宽流量。
在步骤S004中，判断n+2时间步是否有负水深网格存在。
所述网格流出率修正具体步骤包括：若步骤S003中计算出一网格的水深H小于0，即该网格为负水深网格，则引入一修正系数
Figure 729093DEST_PATH_IMAGE018
重新对该负水深网格的水量进行修正计算，直到该负水深网格修正为实际水深后，再转入步骤S002通过所述单宽流量公式获得所述二维网格模型中各网格在下一时间步的单宽流量M和单宽流量N，该下一时间步表示一种循环关系，即开始下一个计算周期。
在步骤S005中，结合图2，其中，当一网格出现负水深时，根据负水深网格在当前时间步n的单宽流量M和单宽流量N计算出该负水深网格的流入总量
Figure 203937DEST_PATH_IMAGE033
和流出总量
Figure 776476DEST_PATH_IMAGE034
；
所述流入总量公式：
Figure 20376DEST_PATH_IMAGE035
     （6）
所述流出总量公式：
Figure 475628DEST_PATH_IMAGE036
   （7）。
在步骤S006中，用于所述计算修正系数
Figure 188500DEST_PATH_IMAGE018
，该计算修正系数λ的公式：        （8）。
在步骤S007中，根据公式（1）、（2）、（3）、（4）、（6）、（7）修正当前负水深网格在n+1时间步的流入量
Figure 853017DEST_PATH_IMAGE038
和流出量
Figure 292219DEST_PATH_IMAGE039
，即，按照公式（1）、（2）、（3）、（4）获得的在n+1时间步的单宽流量M和单宽流量N分别代入公式（6）、（7）。
在步骤S008中，根据上述步骤S007和S006中的相应计算，得出所述负水深网格在n+2时间步的实际水深：
Figure 741655DEST_PATH_IMAGE040
      （9）
在步骤S009中，输出网格的实际水深。
在步骤S010中，转入步骤S002通过所述单宽流量公式获得所述二维网格模型中各网格在下一时间步的单宽流量M和单宽流量N。
显然，上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而这些属于本发明的精神所引伸出的显而易见的变化或变动仍处于本发明的保护范围之中。

Claims (1)

1.一种面向山洪演进数值模拟的计算网格流出率的修正方法，包括：
①建立计算域内地形的数字高程模型，即，所述计算域按照一定空间步长（                                                
Figure 288387DEST_PATH_IMAGE001
、
Figure 205528DEST_PATH_IMAGE002
）进行分割后得到二维网格模型，定义（i，j）为该二维网格模型中一网格，其中，设定该二维网格模型中任一网格的水深为H，且与该网格按X轴同方向的单宽流量M，与Y轴同方向的单宽流量N；同时设定所述山洪演进数值模拟计算的时间步n，以及该时间步n对应的时间步长
Figure 47582DEST_PATH_IMAGE003
；
②通过单宽流量公式获得所述二维网格模型中各网格在n+1时间步的单宽流量M和单宽流量N，其方法，包括：
预设初始条件，即，在n时间步，且设定位于所述计算域的边界上的流量输入起始界的各边界网格的初值参数，该初值参数包括：初始水深
Figure 442791DEST_PATH_IMAGE004
，所述单宽流量M的对应流速矢量
Figure 70213DEST_PATH_IMAGE005
，以及所述单宽流量N的对应流速矢量
Figure 892675DEST_PATH_IMAGE006
；在所述单宽流量公式中位于所述计算域外的网格的初值参数与该网格相邻的一边界网格的初值参数相同，且位于所述计算域内的网格的初值参数的相应取值为0；
以及，所述单宽流量M和单宽流量N的初始值的计算公式：
Figure 222025DEST_PATH_IMAGE007
     （1）
Figure 420926DEST_PATH_IMAGE008
     （2）
所述单宽流量公式：
Figure 902854DEST_PATH_IMAGE010
（3）；
Figure 958534DEST_PATH_IMAGE012
  （4）；
③根据所述各网格在n+1时间步的单宽流量M和单宽流量N，建立所述各网格的水深计算式，即
Figure 525913DEST_PATH_IMAGE014
 （5）
其中， n+2时间步作为水深计算的起始时间步，H表示所述网格水深，
Figure 325242DEST_PATH_IMAGE015
表示在所述起始时间步时网格（i，j）的区域内的水深；表示在n时间步的所述网格的区域内的水深，且该时间步的水深为预设值；表示在n+1时间步的所述网格在X轴方向上流体的流入单宽流量；
Figure 330216DEST_PATH_IMAGE018
表示相邻网格在X轴方向上的流入单宽流量，即，所述网格在X轴方向上的流出单宽流量； 
Figure 667657DEST_PATH_IMAGE019
表示在n+1时间步的所述网格在Y轴方向上流体的流入单宽流量；
Figure 124177DEST_PATH_IMAGE020
表示相邻网格在Y轴方向上的流入单宽流量，即，所述网格在Y轴方向上的流出单宽流量； 
④若步骤③计算出一网格的水深H小于0，即该网格为负水深网格，则引入一修正系数重新对该负水深网格的水量进行修正计算，直到该负水深网格修正为实际水深后，再转入步骤②通过所述单宽流量公式获得所述二维网格模型中各网格在下一时间步的单宽流量M和单宽流量N；
其中，当一网格出现负水深时，根据负水深网格在当前时间步n的单宽流量M和单宽流量N计算出该负水深网格的流入总量
Figure 250582DEST_PATH_IMAGE022
和流出总量
Figure 142446DEST_PATH_IMAGE023
；
所述流入总量公式：
Figure 702740DEST_PATH_IMAGE024
     （6）；
所述流出总量公式：
Figure DEST_PATH_IMAGE025
     （7）；
所述计算修正系数λ的公式：
        （8）；
并根据公式（1）、（2）、（3）、（4）、（6）、（7）修正当前负水深网格在n+1时间步的流入量和流出量
Figure 993017DEST_PATH_IMAGE028
，得出所述负水深网格在n+2时间步的实际水深：
Figure DEST_PATH_IMAGE029
      （9）。
CN201310478469.1A 2013-10-14 2013-10-14 面向山洪演进数值模拟的计算网格流出率的修正方法 Withdrawn CN103530462A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310478469.1A CN103530462A (zh) 2013-10-14 2013-10-14 面向山洪演进数值模拟的计算网格流出率的修正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310478469.1A CN103530462A (zh) 2013-10-14 2013-10-14 面向山洪演进数值模拟的计算网格流出率的修正方法

Publications (1)

Publication Number Publication Date
CN103530462A true CN103530462A (zh) 2014-01-22

Family

ID=49932469

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310478469.1A Withdrawn CN103530462A (zh) 2013-10-14 2013-10-14 面向山洪演进数值模拟的计算网格流出率的修正方法

Country Status (1)

Country Link
CN (1) CN103530462A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107451372A (zh) * 2017-08-09 2017-12-08 中国水利水电科学研究院 一种运动波与动力波相结合的山区洪水过程数值模拟方法
CN114580283A (zh) * 2022-03-07 2022-06-03 水利部交通运输部国家能源局南京水利科学研究院 山溪性强潮河口分汊段中长期动力地貌演变数值模拟方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107451372A (zh) * 2017-08-09 2017-12-08 中国水利水电科学研究院 一种运动波与动力波相结合的山区洪水过程数值模拟方法
CN114580283A (zh) * 2022-03-07 2022-06-03 水利部交通运输部国家能源局南京水利科学研究院 山溪性强潮河口分汊段中长期动力地貌演变数值模拟方法
CN114580283B (zh) * 2022-03-07 2023-08-01 水利部交通运输部国家能源局南京水利科学研究院 山溪性强潮河口分汊段中长期动力地貌演变数值模拟方法

Similar Documents

Publication Publication Date Title
Ramly et al. Application of HEC-GeoHMS and HEC-HMS as rainfall–runoff model for flood simulation
Yao et al. Bank erosion and accretion along the Ningxia–Inner Mongolia reaches of the Yellow River from 1958 to 2008
Timbadiya et al. A 1D–2D coupled hydrodynamic model for river flood prediction in a coastal urban floodplain
CN109543275B (zh) 一种城区地表径流二维数值模拟方法
CN103559411A (zh) 一种基于数字高程模型的洪水演进数值模拟计算修正方法
JP2008050903A (ja) 洪水予測方法および洪水予測システム
CN103530461A (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
Formánek et al. Two-dimensional model of Ciliwung river flood in DKI Jakarta for development of the regional flood index map
Efthimiou The role of sediment rating curve development methodology on river load modeling
Danielson Utilizing a high resolution Digital Elevation Model (DEM) to develop a Stream Power Index (SPI) for the Gilmore creek watershed in Winona County, Minnesota
Khalil et al. Floodplain mapping for Indus river: Chashma–Taunsa reach
Huang An effective alternative for predicting coastal floodplain inundation by considering rainfall, storm surge, and downstream topographic characteristics
Wahyudi et al. Simulating on water storage and pump capacity of “Kencing” river polder system in Kudus regency, Central Java, Indonesia
CN103530462A (zh) 面向山洪演进数值模拟的计算网格流出率的修正方法
Aniskin et al. Potential impact of land-use changes on river basin hydraulic parameters subjected to rapid urbanization
CN103559403A (zh) 洪水演进数值模拟计算中网格流出率的修正方法
Dahm et al. Next generation flood modelling using 3Di: A case study in Taiwan
Krishna et al. Flood frequency analysis of Prakasam barrage reservoir Krishna district, Andhra Pradesh using Weibull, Gringorten and L-moments formula
Malekani et al. Application of GIS in modeling Zilberchai Basin runoff
Ackerman et al. Hydrologic and hydraulic models for performing dam break studies
Schoenbaechler et al. TxBLEND model calibration and validation for the Nueces estuary
Namadi et al. The numerical investigation of the effect of subsequent check dams on flood peaks and the time of concentration using the MIKE 11 modeling system (Case study: golabdareh catchment, Iran)
Andrássy et al. Numerical modeling of groundwater flow close to drinking water resources during flood events
CN117852850B (zh) 基于虚拟河道替代闸孔泄流计算的蓄滞洪区洪水调度方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C04 Withdrawal of patent application after publication (patent law 2001)
WW01 Invention patent application withdrawn after publication

Application publication date: 20140122