CN108897940B - 基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 - Google Patents
基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 Download PDFInfo
- Publication number
- CN108897940B CN108897940B CN201810652670.XA CN201810652670A CN108897940B CN 108897940 B CN108897940 B CN 108897940B CN 201810652670 A CN201810652670 A CN 201810652670A CN 108897940 B CN108897940 B CN 108897940B
- Authority
- CN
- China
- Prior art keywords
- model
- grid
- dimensional hydrodynamic
- hydrodynamic model
- dimensional
- 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
Images
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
-
- 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
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提出了一种基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法,包括以下步骤:步骤1,获取分布式水文模型计算所得的小时尺度净雨量栅格数据;步骤2,采用插值法将小时尺度栅格数据转化为与二维水动力模型时间步长一致的秒尺度栅格数据;步骤3,降尺度后的栅格数据即可作为驱动条件,输入二维水动力模型中质量守恒方程的源项进行水动力计算。本发明方法中基于栅格尺度建立从分布式水文模型到二维水动力模型的单向耦合,保证了水循环过程的连贯性,解决了传统方法中由空间尺度不统一带来的质量不守恒问题。本发明方法还具有计算过程简单、好操作、易推广等优势。
Description
技术领域
本发明涉及水文模型和水动力模型耦合计算技术领域,特别涉及一种基于矩形网格的分布式水文模型和二维水动力模型的耦合方法。
背景技术
分布式水文模型是目前研究流域洪水预报的主要工具。在模型应用中,一般利用DEM生成数字流域,在每个小的子流域(或DEM网格)上应用水文模型模拟产流过程来推求径流值,再进行汇流演算,最后求得每个子流域(或网格)出口断面处的流量过程、峰值流量及洪水到达时间等洪水预报数据。
基于二维浅水动力学方程的水动力模型,常用于洪水过程模拟,可模拟复杂流态的动力过程,更加适用于具有突发性、水量集中、流速大等特点的暴雨洪水模拟,可反映出暴雨洪水特征值在空间和时间上的分布和动态变化。
为综合分布式水文模型和二维水动力模型的优点,采用水文模型和水动力模型耦合的方法进行流域洪水预报。常用的耦合方法可分为单向耦合法和双向耦合法两类。其中,单向耦合法是利用水文模型模拟流域的降雨产流、区域汇流和出流过程,提供区域出流过程作为水动力学模型的输入条件,采用水动力学模型模拟河道或重点区域的洪水过程,即两者之间通过边界条件的衔接,实现水文模型与水动力学模型的耦合。
实际应用和深入研究表明,通过边界条件连接存在一些缺点:(1)两者在空间尺度上存在巨大差别,前者(水文模型)的计算区域为整个流域,后者(水动力模型)的计算区域为河道或流域中的某个区域;(2)两者的空间差异导致纵向和横向汇入后者的净雨量无法完全通过边界条件体现,导致整体质量不守恒;(3)两者机制上的差异导致流域水循环过程模拟上的不连贯性。
发明内容
针对现有方法的不足,本发明针对水文-水动力单向耦合方法在空间尺度、质量守恒和连贯性这三个难点问题,提出一种基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法。
为解决上述技术问题,本发明采用如下技术方案:
一种基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法,包括以下步骤:
步骤1,获取分布式水文模型计算所得的小时尺度净雨量栅格数据ri,j;
步骤4,根据上述方法的具体计算方法如下:
其中,ri,j表示分布式水文模型计算所得的小时尺度净雨量;i,j为栅格序号;n代表当前计算步骤;Δtn为水动力模型当前计算步骤的时间步长;3600代表一小时包含3600秒;表示秒尺度净雨量;为水动力模型质量守恒方程的源项。
本发明的特征在于:分布式水文模型和二维水动力模型基于同一个矩形网格;采用插值法对栅格数据进行降尺度处理;降尺度处理后的净雨量栅格数据可直接输入到二维水动力模型质量守恒方程的源项。
与现有技术相比,本发明具有以下优点和有益效果:
(1)基于栅格尺度完成分布式水文模型和二维水动力模型的耦合,保证了栅格尺度上对水循环模拟的连贯性;
(2)结合二维水动力模型的自适应时间步长,采用插值法对净雨量栅格数据进行降尺度计算,确保了二维水动力模型计算的稳定性并降低了误差;
(3)纵向净雨量直接输入二维水动力模型进行计算,保证了系统整体上的质量守恒。
本发明方法中基于栅格尺度建立从分布式水文模型到二维水动力模型的单向耦合,保证了水循环过程的连贯性,解决了传统方法中由空间尺度不统一带来的质量不守恒问题。本发明方法还具有计算过程简单、好操作、易推广等优势。
附图说明
附图1为本发明基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法的流程图;
附图2为本发明方法具体实例中出口断面模拟流量序列和实测流量序列的对比验证示意图。
具体实施方式
以下通过实例,并结合附图对本发明的技术方案做进一步具体说明。
如附图1所示,基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法,包括如下步骤:
步骤1,数据准备。确定研究区域、所用模型和参数取值后,运行分布式水文模型得到小时尺度净雨量栅格数据ri,j。
具体的,下面对获得小时尺度净雨量栅格数据ri,j的过程进行说明:
首先,获取分布式水文模型和二维水动力模型所需基础数据,基于与地面高程数据DEM栅格数据文件相一致的规则笛卡尔计算网格构建覆盖全流域的水文模型和水动力模型,基于此计算网格提取土地利用类型数据;并根据每个计算栅格的土地利用类型,查表得到各栅格的糙率系数。其中,分布式水文模型所需基础数据,包括:地面高程数据(DEM)、土地利用数据、土壤类型数据、气象数据,其中,气象数据为驱动条件,其它数据为初始条件。
二维水动力模型所需基础数据,包括:地面高程数据、土地利用数据、净雨量,其中净雨量为驱动条件,其它数据为初始条件,且净雨量通过分布式水文模型计算得到的,它是两个模型耦合连接的纽带,也是水文模型和水动力模型耦合的关键条件。
然后,设置水文模型和水动力模型采用与地面高程数据DEM栅格数据文件相一致的规则笛卡尔计算网格,利用ArcGIS将土地利用数据处理为与地面高程数据DEM相同栅格的数据文件;根据每个栅格的土地利用类型,可查经验糙率系数表得到糙率系数(即曼宁系数),糙率系数是二维水动力模型的唯一参数。
然后根据所述的基础数据构建全流域水文模型,包括:采用已有的分布式水文模型计算全流域的产流过程即降雨、下渗、蒸散发,得到每个计算栅格中的净雨量。在本步骤中,计算全流域的产流过程即降雨、下渗、蒸散发,得到每个计算栅格中的净雨量,包括:净雨量=降雨量-下渗量-蒸散发量。本步骤可采用不同的分布式水文模型。虽然不同的分布式水文模型在功能上有所差别,但只要具有计算下渗、蒸散发功能的分布式水文模型即可使用,目的是得到每个计算栅格中的净雨量(小时/天尺度)。
步骤4,根据上述方法的具体计算方法如下:
其中,ri,j表示分布式水文模型计算所得的小时尺度净雨量;i,j为栅格序号;n代表当前计算步骤;Δtn为水动力模型当前计算步骤的时间步长;3600代表一小时包含3600秒;表示秒尺度净雨量;为水动力模型质量守恒方程的源项。
附图2本发明方法具体实例中出口断面模拟流量序列和实测流量序列的对比验证的示意图。
与现有技术相比,本发明具有以下优点和有益效果:
(1)基于栅格尺度完成分布式水文模型和二维水动力模型的耦合,保证了栅格尺度上对水循环模拟的连贯性;
(2)结合二维水动力模型的自适应时间步长,采用插值法对净雨量栅格数据进行降尺度计算,确保了二维水动力模型计算的稳定性并降低了误差;
(3)纵向净雨量直接输入二维水动力模型进行计算,保证了系统整体上的质量守恒。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。
Claims (3)
1.一种基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法,其特征在于,包括如下步骤:
步骤1,获取分布式水文模型计算所得的小时尺度净雨量栅格数据;
步骤2,采用插值法将小时尺度栅格数据转化为与二维水动力模型时间步长一致的秒尺度栅格数据;
步骤3,降尺度后的栅格数据即可作为驱动条件,输入二维水动力模型中质量守恒方程的源项进行水动力计算;
步骤4,根据上述方法的具体计算方法如下:
2.如权利要求1所述的基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法,其特征在于,在所述步骤1中,此处获得的净雨量栅格数据的前提是分布式水文模型与二维水动力模型的构建基于同一个矩形网格,旨在保证两者在空间离散上的一致性。
3.如权利要求1所述的基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法,其特征在于,在所述步骤2中,此处二维水动力模型的时间步长为自适应时间步长,即二维水动力模型完成每一个时间步长后将根据水深和速度的关系重新计算下一个时间步长的大小,因此此处的时间步长为一个变量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810652670.XA CN108897940B (zh) | 2018-06-22 | 2018-06-22 | 基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810652670.XA CN108897940B (zh) | 2018-06-22 | 2018-06-22 | 基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108897940A CN108897940A (zh) | 2018-11-27 |
CN108897940B true CN108897940B (zh) | 2022-07-01 |
Family
ID=64345813
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810652670.XA Active CN108897940B (zh) | 2018-06-22 | 2018-06-22 | 基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108897940B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112733344B (zh) * | 2020-12-30 | 2021-11-02 | 北京大学 | 一种人工化流域产流量确定方法及系统 |
CN112904458B (zh) * | 2021-01-26 | 2021-11-19 | 广东省水文局梅州水文分局 | 一种超渗-蓄满混合产流模式的水文预报方法及系统 |
CN112561212A (zh) * | 2021-02-22 | 2021-03-26 | 航天宏图信息技术股份有限公司 | 山洪风险预测方法、装置、设备及计算机存储介质 |
CN116362146B (zh) * | 2023-02-01 | 2024-04-05 | 重庆交通大学 | 一种水动力联合计算模型构建方法、系统、设备及介质 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7136756B1 (en) * | 2004-11-02 | 2006-11-14 | Vieux And Associates, Inc. | Method for determining runoff |
CN107590354A (zh) * | 2017-10-26 | 2018-01-16 | 成都信息工程大学 | 基于稳定水流场的流域地表径流汇流模拟方法 |
-
2018
- 2018-06-22 CN CN201810652670.XA patent/CN108897940B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7136756B1 (en) * | 2004-11-02 | 2006-11-14 | Vieux And Associates, Inc. | Method for determining runoff |
CN107590354A (zh) * | 2017-10-26 | 2018-01-16 | 成都信息工程大学 | 基于稳定水流场的流域地表径流汇流模拟方法 |
Non-Patent Citations (6)
Title |
---|
二、三维耦合水动力模型研究I:模型的建立;黄玉新等;《水道港口》;20130828(第04期);全文 * |
城市雨洪径流模型研究概述;朱冬冬等;《水资源与水工程学报》;20110615(第03期);全文 * |
基于格子玻尔兹曼方法的流域汇流数值模型;宫兴龙等;《水力发电学报》;20140430;第33卷(第2期);全文 * |
基于混合产流与二维运动波汇流分布式水文模型;包红军等;《水电能源科学》;20161125(第11期);全文 * |
蓄滞洪区洪水资源化的益损定量分析——以大黄堡洼为例;杨小柳等;《水科学进展》;20140930;第25卷(第9期);全文 * |
飞来峡流域基于栅格DEM的分布式水文模拟;万民等;《武汉大学学报(工学版)》;20101028(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN108897940A (zh) | 2018-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108897940B (zh) | 基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 | |
CN102289570B (zh) | 基于降雨-径流-洪水演进计算的洪水预报方法 | |
CN109101706B (zh) | 集总式水文模型和二维水动力模型的耦合方法 | |
Xia et al. | Development of distributed time-variant gain model for nonlinear hydrological systems | |
Xiong et al. | A two-parameter monthly water balance model and its application | |
CN106202790B (zh) | 一种分布式河北模型构建方法及其应用 | |
WO2022032872A1 (zh) | 一种基于大数据的水文预报方法 | |
CN107239657A (zh) | 一种面向对象的水动力学建模要素管理方法 | |
CN108898250B (zh) | 一种基于D藤copula函数的月径流模拟方法 | |
CN113569438B (zh) | 基于多源降雨融合和实时校正的城市洪涝模型构建方法 | |
CN112765912A (zh) | 基于气候模式集合的洪涝灾害社会经济暴露度的评估方法 | |
CN105930635A (zh) | 一种内陆河干旱绿洲区干旱评价系统 | |
Chang et al. | Simulated wave-driven ANN model for typhoon waves | |
CN110765694B (zh) | 一种基于简化型浅水方程组的城市地表水流数值模拟方法 | |
CN112651118A (zh) | 一种气候-陆面-水文过程全耦合模拟方法 | |
CN115983158B (zh) | 一种地下水模型与二维水动力模型松散耦合的方法 | |
CN108920799B (zh) | 基于正方形自适应网格的二维水文-水动力单向耦合方法 | |
CN103870699A (zh) | 基于双层异步迭代策略的水动力学洪水演进模拟方法 | |
Adib et al. | Evaluation of synthetic outlet runoff assessment models | |
Wolfs et al. | Flood probability mapping by means of conceptual modeling | |
Yao et al. | Evaluation of the effects of underlying surface change on catchment hydrological response using the HEC-HMS model | |
Shekar | Rainfall-Runoff Modelling of a River Basin Using HEC HMS: A Review Study | |
Malekani et al. | Application of GIS in modeling Zilberchai Basin runoff | |
CN108345738A (zh) | 一种中小流域暴雨洪水汇流产流模型参数自我率定方法 | |
CN115796599B (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 |