CN109446599A - 应用于大区域水流的计算方法和计算系统 - Google Patents
应用于大区域水流的计算方法和计算系统 Download PDFInfo
- Publication number
- CN109446599A CN109446599A CN201811188575.5A CN201811188575A CN109446599A CN 109446599 A CN109446599 A CN 109446599A CN 201811188575 A CN201811188575 A CN 201811188575A CN 109446599 A CN109446599 A CN 109446599A
- Authority
- CN
- China
- Prior art keywords
- parameter
- flow
- equation
- water information
- grid
- 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.)
- Pending
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 30
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 151
- 239000002689 soil Substances 0.000 claims abstract description 51
- 239000002352 surface water Substances 0.000 claims abstract description 39
- 238000004590 computer program Methods 0.000 claims abstract description 10
- 239000003673 groundwater Substances 0.000 claims description 29
- 238000012545 processing Methods 0.000 claims description 16
- 230000008595 infiltration Effects 0.000 claims description 12
- 238000001764 infiltration Methods 0.000 claims description 12
- 239000010410 layer Substances 0.000 claims description 9
- 238000005259 measurement Methods 0.000 claims description 8
- 229920006395 saturated elastomer Polymers 0.000 claims description 7
- 239000011229 interlayer Substances 0.000 claims description 6
- 230000000694 effects Effects 0.000 abstract description 6
- 238000004088 simulation Methods 0.000 abstract description 6
- 238000000034 method Methods 0.000 description 14
- 238000010586 diagram Methods 0.000 description 8
- 241000208340 Araliaceae Species 0.000 description 5
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 5
- 235000003140 Panax quinquefolius Nutrition 0.000 description 5
- 235000008434 ginseng Nutrition 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 239000013535 sea water Substances 0.000 description 4
- 230000004888 barrier function Effects 0.000 description 3
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical compound O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000006870 function Effects 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 238000006424 Flood reaction Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000005206 flow analysis Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
Classifications
-
- 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
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/14—Pipes
-
- 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)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Networks & Wireless Communication (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种应用于大区域水流的计算方法和计算系统,基于数据高程模型通过计算机程序以计算大区域水流。计算方法包括:由数据高程模型对大区域水流取得多个栅格,其中,各栅格相连成大区域水流图像;划分各栅格为多个子栅格,通过栅格计算方程式来计算各子栅格的地表水信息、土壤水信息或地下水信息;根据多个子栅格中的地表水信息、土壤水信息或地下水信息通过一维流动方程式及/或二维流动方程式来计算多个参数以模拟大区域水流。本发明提供的计算方法,能够提高模拟大区域水流的准确性,进而提高洪涝模拟的计算效果。
Description
技术领域
本发明涉及水流计算技术领域,具体而言,涉及一种应用于大区域水流的计算方法和计算系统。
背景技术
水流类型的划分是水流学研究的基础工作。数据高程模型(DEM,digitalelevation module)为水流类型划分提供了大量的基础数据,数字水流分析和水流类型的自动划分已成为研究热点。现有技术上,DEM设计实现了针对水流类型的提取方法和流程。
一般来说,大面积范围的洪涝风险分析或洪涝灾情预报是必要的,洪涝灾情预报可以给民众预防灾害。原因是不少极端大量降雨事件影响的地域范围很大,例如,台风或暴雨带来的强降雨过程。由于流域是个有机整体,需将整个流域作为研究对象,才利于集成的计算与分析。例如,流域的上游地区与下游地区是密切相关的,上游地区的洪水会直接影响下游地区的汛情。下游地区的河流水位高程或出海口海平面高程,会通过回水效应影响上游地区的河流流动。
发明内容
鉴于上述问题,本发明提供了一种应用于大区域水流的计算方法和计算系统,采用高分辨率的栅格计算单元对大区域水流进行预测的研究,能够提高水流预测的准确性。
为了实现上述目的,本发明采用如下的技术方案:
第一方面,本发明提供了一种应用于大区域水流的计算方法,基于数据高程模型通过计算机程序以计算大区域水流,所述计算方法包括:
由所述数据高程模型对大区域水流取得多个栅格,其中各栅格相连成大区域水流图像;
划分各栅格为多个子栅格,通过栅格计算方程式来计算各子栅格的地表水信息、土壤水信息或地下水信息;
根据所述多个子栅格中的所述地表水信息、所述土壤水信息或所述地下水信息通过一维流动方程式及/或二维流动方程式来计算多个预设参数以模拟大区域的水流。
作为一种可选的实施方式,在根据所述多个子栅格中的所述地表水信息、所述土壤水信息或所述地下水信息通过一维流动方程式及/或二维流动方程式来计算多个预设参数以模拟大区域水流中,所述计算方法包括:
计算地表产流参数,根据所述子栅格中的所述地表水信息和所述土壤水信息来计算所述地表产流参数,所述地表产流参数包括蓄满产流因子和超渗产流因子;或
计算坡面汇流参数,通过所述一维流动方程式比较下游水位因子及本地水位因子来计算所述坡面汇流参数。
作为一种可选的实施方式,在根据所述多个子栅格中的所述地表水信息、所述土壤水信息或所述地下水信息通过一维流动方程式及/或二维流动方程式来计算多个预设参数以模拟大区域水流中,所述计算方法包括:
计算地表积水流动参数,通过至少一个所述一维流动方程式来计算所述地表积水流动参数;或
计算排水管网汇流参数,通过所述一维流动方程式或所述二维流动方程式来计算所述排水管网汇流参数。
作为一种可选的实施方式,在根据所述多个子栅格中的所述地表水信息、所述土壤水信息或所述地下水信息通过一维流动方程式及/或二维流动方程式来计算多个预设参数以模拟大区域水流中,所述计算方法包括:
计算地下水流动参数,根据所述子栅格中的土壤水信息及所述地下水信息来计算所述地下水流动参数,其中,所述地下水流动参数为计算地表下的上部非饱和层及下部饱和层间的变化因子;或
计算河道汇流参数,通过所述二维流动方程式来计算所述河道汇流参数,其中,所述河道汇流参数为计算河道的地表水及地下水的交换因子。
作为一种可选的实施方式,在根据所述多个子栅格中的所述地表水信息、所述土壤水信息或所述地下水信息通过一维流动方程式及/或二维流动方程式来计算多个预设参数以模拟大区域水流中,所述计算方法包括:
计算海水倒灌参数,通过所述二维流动方程式来比较海平面高程因子及沿海平面高程因子以计算所述海水倒灌参数。
第二方面,本发明提供了一种应用于大区域水流的计算系统,基于数据高程模型通过计算机程序以计算大区域水流,所述计算系统包括:
测量模块,用于通过所述数据高程模型对大区域水流取得多个栅格,其中各栅格相连成大区域水流图像;
划分模块,用于划分各栅格为多个子栅格,通过栅格计算方程式来计算各子栅格的地表水信息、土壤水信息或地下水信息;
处理模块,根据所述多个子栅格中的所述地表水信息、所述土壤水信息或所述地下水信息通过一维流动方程式及/或二维流动方程式来计算多个预设参数以模拟大区域水流,其中,所述划分模块连接所述测量模块及所述处理模块。
作为一种可选的实施方式,所述处理模块包括:
地表入渗产流单元,用于计算地表产流参数,所述地表入渗产流单元根据所述子栅格中的所述地表水信息和所述土壤水信息来计算所述地表产流参数,所述地表产流参数包括蓄满产流因子和超渗产流因子;或
坡面汇流单元,用于计算坡面汇流参数,所述坡面汇流单元通过所述一维流动方程式比较下游水位因子及本地水位因子来计算所述坡面汇流参数。
作为一种可选的实施方式,所述处理模块包括:
地表积水流动单元,用于计算地表积水流动参数,所述地表积水流动单元通过至少一个所述一维流动方程式来计算所述地表积水流动参数;或
排水管网汇流单元,用于计算排水管网汇流参数,排水管网汇流单元通过所述一维流动方程式或所述二维流动方程式来计算所述排水管网汇流参数。
作为一种可选的实施方式,所述处理模块包括:
地下水流动单元,用于计算地下水流动参数,所述地下水流动单元根据所述子栅格中的所述土壤水信息及所述地下水信息来计算所述地下水流动参数,其中,所述地下水流动参数为计算地表下的上部非饱和层及下部饱和层间的变化因子;或
河道汇流单元,用于计算河道汇流参数,所述河道汇流单元通过所述二维流动方程式来计算所述河道汇流参数,其中,所述河道汇流参数为计算河道的地表水及地下水的交换因子。
作为一种可选的实施方式,所述处理模块包括:
海水倒灌单元,用于计算海水倒灌参数,所述海水倒灌单元通过所述二维流动方程式来比较海平面高程因子及沿海平面高程因子以计算所述海水倒灌参数。
根据本发明提供的应用于大区域水流的计算方法和计算系统;由所述数据高程模型对大区域水流取得多个栅格;划分各栅格为多个子栅格;根据所述多个子栅格中的地表水信息、土壤水信息或地下水信息来模拟大区域水流的趋势信息。可见,实施本发明的技术方案能通过同时考虑地表入渗产流参数、坡面汇流参数、地表积水流动参数、排水管网汇流参数、地下水流动参数、河道汇流参数及高程变异系数海水倒灌参数等等的预设参数组合来提高预测水流趋势的准确性,从而提高水流环境预测的实时性和准确性。为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,应当理解,以下附图仅示出了本发明的某些实施例,因此不应被看作是对本发明范围的限定。
图1a是本发明实施例1提供应用于大区域水流的计算方法的流程示意图。
图1b是本发明实施例1提供大区域水流图像示意图。
图1c是本发明实施例1提供地下水深度的栅格示意图。
图2是本发明实施例2提供应用于大区域水流的计算系统的方块示意图。
主要元件符号说明:
200-应用于大区域水流的计算系统;210-测量模块;220-划分模块;230-处理模块;231-地表入渗产流单元;232-坡面汇流单元;233-地表积水流动单元;234-排水管网汇流单元;235-地下水流动单元;236-河道汇流单元;237-海水倒灌单元。
具体实施方式
下面将结合本发明实施例中附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此,以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围,而是仅仅表示本发明的选定实施例。基于本发明的实施例,本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
针对现有技术中的问题,本发明提供了一种应用于大区域水流的计算方法和计算系统;由数据高程模型对大区域水流取得多个栅格,其中各栅格相连成大区域的水流图像,通过栅格计算方程式来计算各栅格的地表水信息、土壤水信息或地下水信息;划分各栅格为多个子栅格;根据多个子栅格中的地表水信息、土壤水信息或地下水信息通过一维流动方程式及/或二维流动方程式来计算多个参数以模拟大区域水流。可见,实施本发明的技术方案能够提高预测水流趋势的准确性,同时考虑多个关于水流参数的组合,从而提高水流环境预测的实时性和准确性。并且,该技术可以采用相关的软件或硬件实现,下面通过实施例进行描述。
实施例1
请同时参阅图1a、图1b及图1c,图1a是本发明实施例1提供应用于大区域水流的计算方法(以下简称“计算方法”)的流程示意图。图1b是本发明实施例1提供大区域水流图像示意图。图1c是本发明实施例1提供地下水深度的栅格示意图。计算方法基于数据高程模型(DEM,digital elevation module)通过计算机程序以计算大区域水流。如图1所示,该计算方法包括以下步骤:
S101、由数据高程模型对大区域水流取得多个栅格,其中各栅格相连成大区域水流图像。
本一实施例中,并且本发明实施例中对栅格的获取可以从网络中进行获取,或是从近地卫星上高分辨率的测量系统内置的计算机程序中获取。进一步来说,数据高程模型通过栅格计算方程式来计算各栅格的地表水信息、土壤水信息或地下水信息,通过近地卫星上的测量系统来测量每个栅格具有多少的地表水、多少的土壤水及多少的地下水。在图1b中,由多个栅格相连成大区域水流图像,例如,图中灰阶的深浅可以表示不同的河道高程,深灰阶表示较高的河道高程,浅灰阶表示较低的河道高程。
S103、划分各栅格为多个子栅格,通过栅格计算方程式来计算各子栅格的地表水信息、土壤水信息或地下水信息。
在一实施例中,栅格为大栅格,子栅格为小栅格,由数据高程模型通过近地卫星上高分辨率的测量系统内置的计算机程序来将各大栅格划分为小栅格。例如,图1c代表栅格的地下水深度,负数越小代表地下水深度越深,负数越靠近零代表地下水深度越浅,将每栅格代表30米的量尺划分为每子栅格代表10米的量尺,每子栅格代表10米的量尺相对于每栅格代表30米的量尺可以显示出更精密的下水信息。在图1c中,整体栅格的平均地下水深度值为-5,中间子栅格的地下水深度值相对于邻近子栅格的地下水深度值,可以影响子栅格间的水流流向。
S105、根据多个子栅格中的地表水信息、土壤水信息或地下水信息通过模拟一维流动方程式及/或二维流动方程式来计算多个预设参数以大区域水流。
在一实施例中,通过近地卫星上高分辨率的测量系统内置的计算机程序来模拟大区域水流可以考虑降雨、地表入渗、地表产流、坡面汇流、河道汇流、地下水、海水潮汐、台风等等影响因素。例如,一维流动方程式可以为圣维南方程式。一维流动方程式可以包括一个水量平衡方程式和一个动量方程式。二维流动方程式可以包括一个水量平衡方程式和两个动量方程式。动量方程式可以包括运动波方法、扩散波方法、惯性方法或动力波方法。简单来说,水量平衡方程式需考虑流域中输入的水及输出的水,水量平衡方程式可针对为大气领域、水流领域、土壤领域、地下水领域来进行计算。其中,数据高程模型根据多个子栅格中的地表水信息、土壤水信息或地下水信息,接着通过一维流动方程式及/或二维流动方程式计算多个参数以模拟大区域水流,计算多个预设参数可以包括计算地表产流参数、坡面汇流参数、地表积水流动参数、排水管网汇流参数、地下水流动参数、河道汇流参数或海水倒灌参数。以下为计算各预设参数的说明。
在另一实施例中,计算地表产流参数,数据高程模型根据子栅格中的地表水信息和土壤水信息来计算地表产流参数,地表产流参数包括蓄满产流因子和超渗产流因子。蓄满产流因子代表地表土壤的饱和蓄水度,超渗产流因子代表地表土壤的蓄水度(已超过地表土壤的饱和蓄水度)。例如,当超渗产流因子超过一阀值时,代表降水强度超过地表土壤的饱和蓄水度。当降水强度较大而发生洪涝灾害时,容易发生超渗产流。例如,同时参照蓄满产流因子和超渗产流因子入渗到土壤里的水量,能够更符合实际模拟的蓄满产流和超渗产流。另外,当地表积水时,积水可能通过地表入渗进入到土壤,因此,洪涝过程与地表水信息和土壤水信息是高度正相关的。
在又一实施例中,计算坡面汇流参数,通过一维流动方程式比较下游水位因子及本地水位因子来计算坡面汇流参数。
例如,在计算坡面汇流参数时,坡面汇流参数相关于水位因子,当下游水位因子高于或等于本地水位因子时,回水效应明显;当下游水位因子低于本地水位因子时,回水效应不明显。例如,坡面汇流参数可以为坡面汇流时间,当坡面水流流程为数十公尺时,坡面汇流时间为数分钟。
在又一实施例中,计算地表积水流动参数,通过至少一个一维流动方程式来计算地表积水流动参数。
进一步来说,当地表排水能力不佳且发生淹水情况时,地势高区域与地势低区域之间就会发生地表积水流动。一般来说,计算地表积水流动参数是通过两维流动方程式,然而,两维流动方程式的计算负荷量较高。在本实施例中,数据高程模型通过两个一维流动方程式来计算地表积水流动参数可以有效降低计算量,提高计算效率,增加计算稳定性。
在又一实施例中,通过一维流动方程式或二维流动方程式来计算排水管网汇流参数。
例如,排水管网汇流参数可以通过一维流动方程式或二维流动方程式来计算,数据高程模型通过运动波方式的一维流动方程式计算排水管网里的水流。排水管网汇流参数可以同时参考计算坡面汇流参数的回水效应,提高水流模拟的准确性。例如,当雨水降落地面后,需要通过排水管网来进行排水,排水管网的密集程度相关于排水管网汇流参数。因此,在城市区域建设中,排水管网的设计需要受到高度重视,以防止大量淹水。
在又一实施例中,计算地下水流动参数,数据高程模型根据子栅格中的土壤水信息及地下水信息来计算地下水流动参数,其中,地下水流动参数为计算地表下的上部非饱和层及下部饱和层间的变化因子。
进一步来说,根据子栅格中的土壤水信息及地下水信息来计算地表下的上部非饱和层及下部饱和层间的变化因子。例如,地表与岩石之间的土壤可以分为上部非饱和层及下部饱和层,上部非饱和层及下部饱和层的界面即是地下潜水水位,此地下潜水水位能够动态变化而产生变化因子。地下潜水可以在相邻栅格或子栅格间进行流动,流动方向取决于地下潜水水位的高低关系,地下潜水的流动为由高潜水水位流向低潜水水位。
在又一实施例中,计算河道汇流参数,数据高程模型通过二维流动方程式来计算河道汇流参数,其中,河道汇流参数为计算河道的地表水及地下水的交换因子。
河道汇流参数是通过二维流动方程式来计算。相较于一维流动方程式,二维流动方程式具有更高的精密度,更准确计算河道单元的地表水及河道两岸的地表水及地下水的交换因子。例如,河道水位高程超过堤防高程时,河道两岸发生的淹没动态过程即需要通过二维流动方程式来模拟计算。例如,河道汇流参数的影响要素可以为降雨、地形高程差、流道障碍物。降雨越大,河道汇流越急;降雨越小,河道汇流越缓。地形高程差越大,河道汇流越急;地形高程差越小,河道汇流越缓。流道障碍物越少,河道汇流越急;流道障碍物越多,河道汇流越缓。
在又一实施例中,计算海水倒灌参数,数据高程模型通过二维流动方程式来比较海平面高程因子及沿海平面高程因子以计算海水倒灌参数。
进一步来说,数据高程模型可以通过二维流动方程式根据海平面高程的变化来实现对海水倒灌参数的模拟。如果海平面高程增加,海平面高程因子高于沿海平面高程因子,海水即倒灌进入沿海或内河。例如,当台风引起海平面高程升高,海水即流入沿海低地,可以通过数据高程模型对海水倒灌参数的模拟。
实施例2
请参阅图2,图2是本发明实施例2提供应用于大区域水流的计算系统(以下简称“计算系统”)的方块示意图。基于数据高程模型通过计算机程序以计算大区域水流。计算系统200包括:
测量模块210,用于通过数据高程模型对大区域水流取得多个栅格,其中各栅格相连成大区域水流图像。
划分模块220,用于划分各栅格为多个子栅格,通过栅格计算方程式来计算各子栅格的地表水信息、土壤水信息或地下水信息。
处理模块230,用于根据多个子栅格中的地表水信息、土壤水信息或地下水信息通过一维流动方程式及/或二维流动方程式来计算多个参数以模拟大区域水流,其中,划分模块220连接测量模块210及处理模块230。
处理模块230进一步包括地表入渗产流单元231、坡面汇流单元232、地表积水流动单元233、排水管网汇流单元234、地下水流动单元235、河道汇流单元236或海水倒灌单元237。以下为分别说明:
地表入渗产流单元231,用于计算地表产流参数,地表入渗产流单元231根据子栅格中的地表水信息和土壤水信息来计算地表产流参数,地表产流参数包括蓄满产流因子和超渗产流因子。
坡面汇流单元232,用于计算坡面汇流参数,坡面汇流单元232通过一维流动方程式比较下游水位因子及本地水位因子来计算坡面汇流参数。
地表积水流动单元233,用于计算地表积水流动参数,地表积水流动单元233通过至少一个一维流动方程式来计算地表积水流动参数。
排水管网汇流单元234,用于计算排水管网汇流参数,排水管网汇流单元234通过一维流动方程式或二维流动方程式来计算排水管网汇流参数。
地下水流动单元235,用于计算地下水流动参数,地下水流动单元235根据子栅格中的土壤水信息及地下水信息来计算地下水流动参数,其中,地下水流动参数为计算地表下的上部非饱和层及下部饱和层间的变化因子。
河道汇流单元236,用于计算河道汇流参数,河道汇流单元236通过二维流动方程式来计算河道汇流参数,其中,河道汇流参数为计算河道的地表水及地下水的交换因子。
海水倒灌单元237,用于计算海水倒灌参数,海水倒灌单元237通过二维流动方程式来比较海平面高程因子及沿海平面高程因子以计算海水倒灌参数。
可见,图2所描述的计算系统,能够精确地计算得到洪涝模拟,从而提高了计算效率及精密度。另一方面,由测量模块测量并取得多个栅格。由划分模块划分各栅格为多个子栅格。由处理模块根据多个子栅格中的地表水信息、土壤水信息或地下水信息通过一维流动方程式及/或二维流动方程式来模拟大区域水流的趋势信息。可见,实施本发明的技术方案能够提高预测洪涝及大区域水流的准确性,同时考虑地表入渗产流参数、坡面汇流参数、地表积水流动参数、排水管网汇流参数、地下水流动参数、河道汇流参数及高程变异系数海水倒灌参数等等的预设参数,并通过一维流动方程式及/或二维流动方程式来模拟大区域水流的趋势信息,从而提高预测洪涝及大区域水流的准确性。
另外,在本发明各个实施例中的各功能模块或单元可以集成在一起形成一个独立的部分,也可以是各个模块单独存在,也可以两个或更多个模块集成形成一个独立的部分。
所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时,可以存储在一个计算机可读取存储介质中。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来,该计算机软件产品存储在一个存储介质中,包括若干指令用以使得一台计算机设备(可以是智能手机、个人计算机、服务器、或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括:U盘、移动硬盘、只读存储器(ROM,Read-Only Memory)、随机存取存储器(RAM,Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (10)
1.一种应用于大区域水流的计算方法,基于数据高程模型通过计算机程序以计算大区域水流,其特征在于,所述计算方法包括:
由所述数据高程模型对大区域水流取得多个栅格,其中各栅格相连成大区域水流图像;
划分各栅格为多个子栅格,通过栅格计算方程式来计算各子栅格的地表水信息、土壤水信息或地下水信息;
根据所述多个子栅格中的所述地表水信息、所述土壤水信息或所述地下水信息通过一维流动方程式及/或二维流动方程式来计算预设参数以模拟大区域水流。
2.根据权利要求1所述的计算方法,其特征在于,在根据所述多个子栅格中的所述地表水信息、所述土壤水信息或所述地下水信息通过一维流动方程式及/或二维流动方程式来计算多个预设参数以模拟大区域水流中,所述计算方法包括:
计算地表产流参数,根据所述子栅格中的所述地表水信息和所述土壤水信息来计算所述地表产流参数,所述地表产流参数包括蓄满产流因子和超渗产流因子;或
计算坡面汇流参数,通过所述一维流动方程式比较下游水位因子及本地水位因子来计算所述坡面汇流参数。
3.根据权利要求1所述的计算方法,其特征在于,在根据所述多个子栅格中的所述地表水信息、所述土壤水信息或所述地下水信息通过一维流动方程式及/或二维流动方程式来计算多个预设参数以模拟大区域水流中,所述计算方法包括:
计算地表积水流动参数,通过至少一个所述一维流动方程式来计算所述地表积水流动参数;或
计算排水管网汇流参数,通过所述一维流动方程式或所述二维流动方程式来计算所述排水管网汇流参数。
4.根据权利要求1所述的计算方法,其特征在于,在根据所述多个子栅格中的所述地表水信息、所述土壤水信息或所述地下水信息通过一维流动方程式及/或二维流动方程式来计算多个预设参数以模拟大区域水流中,所述计算方法包括:
计算地下水流动参数,根据所述子栅格中的土壤水信息及所述地下水信息来计算所述地下水流动参数,其中,所述地下水流动参数为计算地表下的上部非饱和层及下部饱和层间的变化因子;或
计算河道汇流参数,通过所述二维流动方程式来计算所述河道汇流参数,其中,所述河道汇流参数为计算河道的地表水及地下水的交换因子。
5.根据权利要求1所述的计算方法,其特征在于,在根据所述多个子栅格中的所述地表水信息、所述土壤水信息或所述地下水信息通过一维流动方程式及/或二维流动方程式来计算多个预设参数以模拟大区域水流中,所述计算方法包括:
计算海水倒灌参数,通过所述二维流动方程式来比较海平面高程因子及沿海平面高程因子以计算所述海水倒灌参数。
6.一种应用于大区域水流的计算系统,基于数据高程模型通过计算机程序以计算大区域水流,其特征在于,所述计算系统包括:
测量模块,用于通过所述数据高程模型对大区域水流取得多个栅格,其中各栅格相连成大区域水流图像,通过栅格计算方程式来计算各栅格的地表水信息、土壤水信息或地下水信息;
划分模块,用于划分各栅格为多个子栅格;
处理模块,用于根据所述多个子栅格中的所述地表水信息、所述土壤水信息或所述地下水信息通过一维流动方程式及/或二维流动方程式来计算多个预设参数以模拟大区域水流,其中,所述划分模块连接所述测量模块及所述处理模块。
7.根据权利要求6所述的计算系统,其特征在于,所述处理模块包括:
地表入渗产流单元,用于计算地表产流参数,所述地表入渗产流单元根据所述子栅格中的所述地表水信息和所述土壤水信息来计算所述地表产流参数,所述地表产流参数包括蓄满产流因子和超渗产流因子;或
坡面汇流单元,用于计算坡面汇流参数,所述坡面汇流单元通过所述一维流动方程式比较下游水位因子及本地水位因子来计算所述坡面汇流参数。
8.根据权利要求6所述的计算系统,其特征在于,所述处理模块包括:
地表积水流动单元,用于计算地表积水流动参数,所述地表积水流动单元通过至少一个所述一维流动方程式来计算所述地表积水流动参数;或
排水管网汇流单元,用于计算排水管网汇流参数,所述排水管网汇流单元通过所述一维流动方程式或所述二维流动方程式来计算所述排水管网汇流参数。
9.根据权利要求6所述的计算系统,其特征在于,所述处理模块包括:
地下水流动单元,用于计算地下水流动参数,所述地下水流动单元根据所述子栅格中的土壤水信息及所述地下水信息来计算所述地下水流动参数,其中,所述地下水流动参数为计算地表下的上部非饱和层及下部饱和层间的变化因子;或
河道汇流单元,用于计算河道汇流参数,所述河道汇流单元通过所述二维流动方程式来计算所述河道汇流参数,其中,所述河道汇流参数为计算河道的地表水及地下水的交换因子。
10.根据权利要求6所述的计算系统,其特征在于,所述处理模块包括:
海水倒灌单元,用于计算海水倒灌参数,所述海水倒灌单元通过所述二维流动方程式来比较海平面高程因子及沿海平面高程因子以计算所述海水倒灌参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811188575.5A CN109446599A (zh) | 2018-10-12 | 2018-10-12 | 应用于大区域水流的计算方法和计算系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811188575.5A CN109446599A (zh) | 2018-10-12 | 2018-10-12 | 应用于大区域水流的计算方法和计算系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109446599A true CN109446599A (zh) | 2019-03-08 |
Family
ID=65544789
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811188575.5A Pending CN109446599A (zh) | 2018-10-12 | 2018-10-12 | 应用于大区域水流的计算方法和计算系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109446599A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111914396A (zh) * | 2020-06-30 | 2020-11-10 | 南京大学 | 基于高分辨率dem数据的次网格地形三维地表太阳辐射强迫效应快速参数化方法 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102034001A (zh) * | 2010-12-16 | 2011-04-27 | 南京大学 | 一种以栅格为模拟单元的分布式水文模型设计方法 |
-
2018
- 2018-10-12 CN CN201811188575.5A patent/CN109446599A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102034001A (zh) * | 2010-12-16 | 2011-04-27 | 南京大学 | 一种以栅格为模拟单元的分布式水文模型设计方法 |
Non-Patent Citations (2)
Title |
---|
万民等: "飞来峡流域基于栅格DEM的分布式水文模拟", 《武汉大学学报(工学版)》 * |
包红军等: "基于混合产流与二维运动波汇流分布式水文模型", 《水电能源科学》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111914396A (zh) * | 2020-06-30 | 2020-11-10 | 南京大学 | 基于高分辨率dem数据的次网格地形三维地表太阳辐射强迫效应快速参数化方法 |
CN111914396B (zh) * | 2020-06-30 | 2023-09-08 | 南京大学 | 基于高分辨率dem数据的次网格地形三维地表太阳辐射强迫效应快速参数化方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Fujii et al. | Hydrological roles of the Cambodian floodplain of the Mekong River | |
CN113723024B (zh) | 一种适用于滨海地区的“溪流”-“河道”-“河口”分布式洪水过程模拟方法 | |
Okedu et al. | Harnessing the potential of small hydro power in Cross River state of Southern Nigeria | |
Duong et al. | Assessing climate change impacts on the stability of small tidal inlets: Part 1-Data poor environments | |
Luo et al. | Numerical modelling of hydrodynamics and sand transport in the tide-dominated coastal-to-estuarine region | |
Wu et al. | Hydrological response to climate change and human activities: A case study of Taihu Basin, China | |
CN115935732A (zh) | 一种复杂感潮河网的多源复合灾害淹没风险的模拟方法 | |
Baldoni et al. | A comprehensive study of the river plume in a microtidal setting | |
Ton et al. | Longshore sediment transport by large-scale lake circulations at low-energy, non-tidal beaches: A field and model study | |
Sundborg et al. | Erosion and sedimentation by water: problems and prospects | |
CN109299428A (zh) | 应用运动波的水流计算方法及系统 | |
CN109446596A (zh) | 应用网格自动机的水流计算方法及系统 | |
Thomas et al. | Improving coastal flooding predictions by switching meshes during a simulation | |
Hung et al. | Coastline and river mouth evolution in the central part of the Red River Delta | |
Duong | Climate change induced coastline change adjacent to small tidal inlets | |
CN109446599A (zh) | 应用于大区域水流的计算方法和计算系统 | |
Wang et al. | Impacts of different dynamic factors on the saltwater intrusion in the northern branch of the yangtze estuary | |
Elmoustafa et al. | Defining the degree of flood hazard using a hydrodynamic approach, a case study: Wind turbines field at west of Suez Gulf | |
Schallenberg et al. | A temperate, tidal lake‐wetland complex 1. Water balance and ecological implications | |
Wilkens | Medium scale morphodynamics of the Central Dithmarschen Bight | |
Goodwin et al. | Hydrodynamic simulation of small-scale tidal wetlands | |
CN108984999A (zh) | 应用概念模型的水流计算方法及系统 | |
Hubbert et al. | Modelling storm surges and coastal ocean flooding | |
Huang et al. | Morphodynamic evolution of the Xiaoqing River mouth: a Huanghe River-derived mixed energy estuary | |
Asplund et al. | Coastal erosion in the region of Thu Bon River mouth, Vietnam |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190308 |