CN108897940A - 基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 - Google Patents

基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 Download PDF

Info

Publication number
CN108897940A
CN108897940A CN201810652670.XA CN201810652670A CN108897940A CN 108897940 A CN108897940 A CN 108897940A CN 201810652670 A CN201810652670 A CN 201810652670A CN 108897940 A CN108897940 A CN 108897940A
Authority
CN
China
Prior art keywords
model
dimentional
hydrodynamic model
raster data
scale
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.)
Granted
Application number
CN201810652670.XA
Other languages
English (en)
Other versions
CN108897940B (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.)
Institute of Geographic Sciences and Natural Resources of CAS
Original Assignee
Institute of Geographic Sciences and Natural Resources of CAS
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 Institute of Geographic Sciences and Natural Resources of CAS filed Critical Institute of Geographic Sciences and Natural Resources of CAS
Priority to CN201810652670.XA priority Critical patent/CN108897940B/zh
Publication of CN108897940A publication Critical patent/CN108897940A/zh
Application granted granted Critical
Publication of CN108897940B publication Critical patent/CN108897940B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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；
步骤2，采用插值法将小时尺度栅格数据转化为与二维水动力模型时间步长一致的秒尺度栅格数据
步骤3，降尺度后的栅格数据即可作为驱动条件，输入二维水动力模型中质量守恒方程的源项进行水动力计算，即源项在数值上等于
步骤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相同栅格的数据文件；根据每个栅格的土地利用类型，可查经验糙率系数表得到糙率系数(即曼宁系数)，糙率系数是二维水动力模型的唯一参数。
然后根据所述的基础数据构建全流域水文模型，包括：采用已有的分布式水文模型计算全流域的产流过程即降雨、下渗、蒸散发，得到每个计算栅格中的净雨量。在本步骤中，计算全流域的产流过程即降雨、下渗、蒸散发，得到每个计算栅格中的净雨量，包括：净雨量＝降雨量-下渗量-蒸散发量。本步骤可采用不同的分布式水文模型。虽然不同的分布式水文模型在功能上有所差别，但只要具有计算下渗、蒸散发功能的分布式水文模型即可使用，目的是得到每个计算栅格中的净雨量(小时/天尺度)。
步骤2，数据降尺度。采用插值方法将小时尺度的栅格数据ri,j转化为与二维水动力模型时间步长一致的秒尺度栅格数据
步骤3，输入计算。将降尺度后的栅格数据作为驱动条件，输入二维水动力模型中质量守恒方程的源项进行水动力计算。
步骤4，根据上述方法的具体计算方法如下：
其中，ri,j表示分布式水文模型计算所得的小时尺度净雨量；i，j为栅格序号；n代表当前计算步骤；Δtn为水动力模型当前计算步骤的时间步长；3600代表一小时包含3600秒；表示秒尺度净雨量；为水动力模型质量守恒方程的源项。
附图2本发明方法具体实例中出口断面模拟流量序列和实测流量序列的对比验证的示意图。
与现有技术相比，本发明具有以下优点和有益效果：
(1)基于栅格尺度完成分布式水文模型和二维水动力模型的耦合，保证了栅格尺度上对水循环模拟的连贯性；
(2)结合二维水动力模型的自适应时间步长，采用插值法对净雨量栅格数据进行降尺度计算，确保了二维水动力模型计算的稳定性并降低了误差；
(3)纵向净雨量直接输入二维水动力模型进行计算，保证了系统整体上的质量守恒。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求及其等同限定。

Claims (3)

1.一种基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法，其特征在于，包括如下步骤：
步骤1，获取分布式水文模型计算所得的小时尺度净雨量栅格数据；
步骤2，采用插值法将小时尺度栅格数据转化为与二维水动力模型时间步长一致的秒尺度栅格数据；
步骤3，降尺度后的栅格数据即可作为驱动条件，输入二维水动力模型中质量守恒方程的源项进行水动力计算；
步骤4，根据上述方法的具体计算方法如下：
式(1)为插值法公式(步骤2)，其中，ri,j表示分布式水文模型计算所得的小时尺度净雨量；i，j为栅格序号；n代表当前计算步骤；Δtn为水动力模型当前计算步骤的时间步长；3600代表一小时包含3600秒；表示秒尺度净雨量，式(2)为每个时间步长质量守恒方程源项的计算公式(步骤3)，其中，为水动力模型质量守恒方程的源项。
2.如权利要求1所述的基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法，其特征在于，在所述步骤1中,此处获得的净雨量栅格数据的前提是分布式水文模型与二维水动力模型的构建基于同一个矩形网格，旨在保证两者在空间离散上的一致性。
3.如权利要求1所述的基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法，其特征在于，在所述步骤2中，此处二维水动力模型的时间步长为自适应时间步长，即二维水动力模型完成每一个时间步长后将根据水深和速度的关系重新计算下一个时间步长的大小，因此此处的时间步长为一个变量。
CN201810652670.XA 2018-06-22 2018-06-22 基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法 Active CN108897940B (zh)

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 true CN108897940A (zh) 2018-11-27
CN108897940B 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)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112561212A (zh) * 2021-02-22 2021-03-26 航天宏图信息技术股份有限公司 山洪风险预测方法、装置、设备及计算机存储介质
CN112733344A (zh) * 2020-12-30 2021-04-30 北京大学 一种人工化流域产流量确定方法及系统
CN112904458A (zh) * 2021-01-26 2021-06-04 广东省水文局梅州水文分局 一种超渗-蓄满混合产流模式的水文预报方法及系统
CN116362146A (zh) * 2023-02-01 2023-06-30 重庆交通大学 一种水动力联合计算模型构建方法、系统、设备及介质

Citations (2)

* Cited by examiner, † Cited by third party
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 成都信息工程大学 基于稳定水流场的流域地表径流汇流模拟方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
Title
万民等: "飞来峡流域基于栅格DEM的分布式水文模拟", 《武汉大学学报(工学版)》 *
包红军等: "基于混合产流与二维运动波汇流分布式水文模型", 《水电能源科学》 *
宫兴龙等: "基于格子玻尔兹曼方法的流域汇流数值模型", 《水力发电学报》 *
朱冬冬等: "城市雨洪径流模型研究概述", 《水资源与水工程学报》 *
杨小柳等: "蓄滞洪区洪水资源化的益损定量分析——以大黄堡洼为例", 《水科学进展》 *
黄玉新等: "二、三维耦合水动力模型研究I:模型的建立", 《水道港口》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112733344A (zh) * 2020-12-30 2021-04-30 北京大学 一种人工化流域产流量确定方法及系统
CN112904458A (zh) * 2021-01-26 2021-06-04 广东省水文局梅州水文分局 一种超渗-蓄满混合产流模式的水文预报方法及系统
CN112904458B (zh) * 2021-01-26 2021-11-19 广东省水文局梅州水文分局 一种超渗-蓄满混合产流模式的水文预报方法及系统
CN112561212A (zh) * 2021-02-22 2021-03-26 航天宏图信息技术股份有限公司 山洪风险预测方法、装置、设备及计算机存储介质
CN116362146A (zh) * 2023-02-01 2023-06-30 重庆交通大学 一种水动力联合计算模型构建方法、系统、设备及介质
CN116362146B (zh) * 2023-02-01 2024-04-05 重庆交通大学 一种水动力联合计算模型构建方法、系统、设备及介质

Also Published As

Publication number Publication date
CN108897940B (zh) 2022-07-01

Similar Documents

Publication Publication Date Title
CN108897940A (zh) 基于矩形网格的分布式水文模型和二维水动力模型的单向耦合方法
CN102289570B (zh) 基于降雨-径流-洪水演进计算的洪水预报方法
CN107832931A (zh) 一种平原水网地区内涝风险的模块化分析方法
CN108022047A (zh) 一种海绵城市水文计算方法
CN109101706B (zh) 集总式水文模型和二维水动力模型的耦合方法
CN107704592A (zh) 一种基于WebGIS的洪水预报服务构建方法
CN109145376A (zh) 一种基于cad地形图的水库淹没范围确定方法
CN103886135B (zh) 基于二维非恒定流数值模型的电力工程选址方法
CN109918821A (zh) 一种迎风守恒型河道漫溢出流数值模拟方法
CN108021780A (zh) 一种基于无规则非结构网格模型的山洪动态仿真方法
CN113505546A (zh) 一种洪水风险预测系统
CN112651118B (zh) 一种气候-陆面-水文过程全耦合模拟方法
CN113569438A (zh) 基于多源降雨融合和实时校正的城市洪涝模型构建方法
CN109858666A (zh) 一种湖泊水资源量评估与预测方法
Poretti et al. An approach for flood hazard modelling and mapping in the medium Valtellina
Li et al. Development of 1D and 2D coupled model to simulate urban inundation: an application to Beijing Olympic Village
CN115186040A (zh) 一种高效简便界定山丘区中小河流洪水淹没范围的方法
Adib et al. Flood inundation modeling for Kota Tinggi catchment by combination of 2D hydrodynamic model and flood mapping approach
CN109190213A (zh) 基于土地利用类型的地表径流流域空间离散化方法
CN108920799B (zh) 基于正方形自适应网格的二维水文-水动力单向耦合方法
Malekani et al. Application of GIS in modeling Zilberchai Basin runoff
CN109063306A (zh) 一种网格化河北模型的土壤下渗能力空间离散方法
CN114459571B (zh) 一种基于初态水位曲线簇河道初始水位确定方法
CN112948915B (zh) 一种用于数值模拟试验中竖直涉水建筑物的概化处理方法
CN117852850B (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