CN110442988A - 一种基于元胞自动机的城市地表径流流向计算方法及装置 - Google Patents

一种基于元胞自动机的城市地表径流流向计算方法及装置 Download PDF

Info

Publication number
CN110442988A
CN110442988A CN201910737873.3A CN201910737873A CN110442988A CN 110442988 A CN110442988 A CN 110442988A CN 201910737873 A CN201910737873 A CN 201910737873A CN 110442988 A CN110442988 A CN 110442988A
Authority
CN
China
Prior art keywords
cellular
water
downstream
neighborhood
level value
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
CN201910737873.3A
Other languages
English (en)
Other versions
CN110442988B (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.)
Nanjing Institute Of Water Conservancy Sciences State Energy Bureau Ministry Of Transportation Ministry Of Water Conservancy
Nanjing Hydraulic Research Institute of National Energy Administration Ministry of Transport Ministry of Water Resources
Original Assignee
Nanjing Institute Of Water Conservancy Sciences State Energy Bureau Ministry Of Transportation Ministry Of Water Conservancy
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 Institute Of Water Conservancy Sciences State Energy Bureau Ministry Of Transportation Ministry Of Water Conservancy filed Critical Nanjing Institute Of Water Conservancy Sciences State Energy Bureau Ministry Of Transportation Ministry Of Water Conservancy
Priority to CN201910737873.3A priority Critical patent/CN110442988B/zh
Publication of CN110442988A publication Critical patent/CN110442988A/zh
Application granted granted Critical
Publication of CN110442988B publication Critical patent/CN110442988B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

本发明公开了一种基于元胞自动机的城市地表径流流向计算方法及装置，该方法包括：建立研究区的元胞基础地理场景；根据水位值大小确定中间元胞的下游元胞，计算中间元胞水量转移到各个下游元胞所需时间；比较中间元胞与其邻域元胞之间水位值标准差和标准差阈值相对大小，选择水量分配权重值，确定各个下游元胞中的水量分配比例；比较模型时间步长与水量实际转移时间的大小，确定一个时间步长内的水量分配比例；根据中间元胞水量与水量分配比例计算中间元胞向各个下游元胞转移的水量，并设置各个下游方向转移上限，修正转移至各个下游元胞的水量。本发明方法更加贴合实际的流向运动过程，能够弥补其他多流向算法的不足，提高城市汇流过程模拟精度。

Description

一种基于元胞自动机的城市地表径流流向计算方法及装置
技术领域
本发明涉及城市水文和环境科学领域，尤其涉及一种基于元胞自动机的城市地表径流流向计算方法及装置。
背景技术
城市水文模型可以实现城市降雨径流过程的精确模拟，为城市内涝预警、雨洪应急疏散方案制定以及城市面源污染防控和治理等提供理论依据。城市地表径流流向能够描述水、沉积物等地表物质在地貌单元之间从高到低的传输与流动路线，决定着汇水面积、地形指数等重要地形和水文参数的计算，在城市水文模型研究中起着重要作用。
流向算法分为单流向算法和多流向算法两大类。单流向算法原理是水流将水力梯度最大的相邻单元作为下游单元，对DEM高程数据比较敏感，不适合模拟缓坡坡面漫散流动。多流向算法则是按照水量分配权重p将水量分配至不同的下游单元中，一般基于两种假设，其一是“越陡下游方向获取的水量越多”，其二则基于“动力系统倾向于往平衡状态发展”的假设，但目前多流向算法中p值为固定值，忽略了地形特征对流场划分的影响，在坡度较陡情况下模拟结果较差。因此，本发明提出了一种基于元胞自动机的城市地表径流流向计算方法，结合多流向算法的两种基本假设，计算3×3单元水位值标准差作为直接反映局部地形变化程度的指标，并将此作为水量分配权重p的取值依据，以弥补其他多流向算法的不足，提高城市汇流过程模拟精度，这也本发明的主要研究课题。
发明内容
发明目的：本发明目的是提供一种基于元胞自动机的城市地表径流流向计算方法及装置，更加贴合实际的流向运动过程，为城市地表径流流向计算提供新方法。
技术方案：为实现上述发明目的，本发明采用如下技术方案：
一种基于元胞自动机的城市地表径流流向计算方法，包括如下步骤：
(1)建立研究区的基于元胞自动机模式的基础地理场景，其中元胞的邻域关系为Moore型八邻域；
(2)比较中间元胞与邻域元胞的水位值大小，排除水位值高于中间元胞的邻域元胞，将剩下的邻域元胞作为水量转移的下游元胞，计算中间元胞水量转移到各个下游元胞所需时间；
(3)计算中间元胞与其八个邻域元胞之间水位值的标准差及标准差阈值，对比计算标准差和标准差阈值相对大小，选择水量分配权重值，确定各个下游元胞中的水量分配比例；
(4)比较模型时间步长与水量实际转移时间的大小，确定一个模型计算时间步长内的水量分配比例；
(5)根据中间元胞水量与水量分配比例计算中间元胞向各个下游元胞转移的水量，并设置各个下游方向转移上限，修正转移至各个下游元胞的水量。
进一步地，步骤(2)中，根据如下公式计算中间元胞水量转移到各个下游元胞所需时间：
其中，Ti为中间元胞水量转移到第i个下游元胞所需时间，depth0为中间元胞水深，si为中间元胞与第i个下游元胞之间的水面坡度，ni为曼宁糙率系数，Di为中间元胞与第i个下游元胞之间的距离。
进一步地，步骤(3)中，根据如下公式确定各个下游元胞中的水量分配比例：
其中，di为中间元胞对第i个下游邻域元胞的水量分配比例，tanβi为中间元胞到第i个下游邻域元胞的坡度比降，Li为对第i个下游邻域元胞的等高线长度加权因子，M为下游邻域元胞的数量，p为水量分配权重，SD和SD0为中间元胞与其八个邻域元胞之间水位值的标准差及标准差阈值。
进一步地，中间元胞与其八个邻域元胞之间水位值的标准差SD及标准差阈值SD0为：
其中，hk为第k个元胞的水位值，为中间元胞与其八个邻域元胞的水位值均值，N＝9，表示八个邻域和一个中间元胞，dCA为元胞边长。
进一步地，步骤(4)中，对比模型时间步长t与步骤(3)中得到的水量转移至各个下游元胞中的时间Ti的大小，若Ti大于时间步长t，只有部分水量可以被转移，此时向下游转移水量比例为如果Ti小于时间步长t，认为水量可以全部转移，此时向下游元胞转移的水量比例为di。
进一步地，步骤(5)中，根据如下公式对中间元胞向各个下游元胞转移的水量进行修正：
其中，Pi为中间元胞向第i个下游元胞转移的水量，hi为第i个下游元胞水位值，为中间元胞与第i个下游元胞的平均水位值，A为元胞面积。
基于相同的发明构思，本发明提供的一种基于元胞自动机的城市地表径流流向计算装置，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述计算机程序被加载至处理器时实现所述的基于元胞自动机的城市地表径流流向计算方法。
有益效果：本发明提供的一种基于元胞自动机的城市地表径流流向计算方法，基于“越陡下游方向获取的水量越多”假设，通过计算中间元胞与其邻域元胞之间水位值的标准差及标准差阈值，对比计算标准差和标准差阈值相对大小，并以此作为选择水量分配权重值的取值依据，确定各个下游元胞中的水量分配比例，再结合模型时间步长与水量实际转移时间的大小，确定一个模型计算时间步长内的最终水量分配比例；并且基于“动力系统倾向于往平衡状态发展”假设，设置各个下游方向转移上限，修正转移至各个下游元胞的水量。与现有技术相比，本发明的计算方法更加贴合实际的流向动态过程，能够弥补其他多流向算法的不足，提高城市汇流过程模拟精度。
附图说明
图1是本发明的一个实施例的流程示意图；
图2是本发明的元胞空间领域关系示意图；
图3是典型城市地表径流流向判断算法流程示意图；
图4是流向计算算法流程示意图。
具体实施方式
如图1所示，本发明实施例公开的一种基于元胞自动机的城市地表径流流向计算方法，包括以下步骤：
(1)获取研究区的下垫面类型、数字高程数据和水深，设置元胞大小、时间步长，建立基于元胞自动机模式的基础地理场景。
该步骤具体包括：
获取研究区的下垫面类型和数字高程数据(DEM)构建元胞基础地理场景，获得每个元胞的下垫面类型(例如：道路、屋面和草地等)和高程，定义每个元胞的大小、曼宁糙率系数、坐标、邻域关系，模型的时间步长和模拟时长，相关参数如表1所示。如图2所示，本发明中元胞的邻域关系为Moore型八邻域。
表1城市元胞自动机模型参数示例表
参数 意义 说明
d<sub>CA</sub> 元胞边长 模型设置
DEM 元胞高程 资料获取
depth 元胞水深/水量 资料获取
n 曼宁糙率 下垫面材质平均糙率
t 时间步长 设置模拟时间的步长
(2)比较中间元胞与邻域元胞的水位值大小，排除水位值高于中间元胞的邻域元胞，将剩下的邻域元胞作为水量转移的下游元胞；利用曼宁公式计算中间元胞水量全部转移到各个下游元胞所需时间。
该步骤具体包括：
(2-1)水量转移下游元胞判断：计算每个元胞对应的水位值h：
h＝DEM+depth
其中，DEM为元胞的高程，m；depth为元胞水深，m；h为每个元胞的水位值，m。比较中间元胞与邻域元胞的水位大小，排除水位值高于中间元胞的邻域元胞，将剩下的邻域元胞作为水量转移的下游元胞。
(2-2)水量转移时间计算：
其中，Ti为中间元胞水量转移到第i个下游元胞所需时间，s；Vi为流速，m/s；depth0为中间元胞水深，m；si为中间元胞与第i个下游元胞之间的水面坡度；ni为曼宁糙率系数，s/m1/3；Di为中间元胞与第i个下游元胞之间的距离，对于垂直邻域元胞而言，D为元胞边长(D＝dCA)，对于对角邻域元胞，D为倍的元胞边长m。
(3)计算中间元胞与其八个邻域元胞之间水位值的标准差，计算标准差阈值，对比计算标准差和标准差阈值相对大小，选择水量分配权重p值，确定各个下游元胞中的水量分配比例。
该步骤具体包括：
(3-1)标准差和标准差阈值计算：根据中间元胞与其八个邻域元胞之间水位值，计算标准差SD和标准差阈值SD0。
其中，SD为中间元胞和邻域八个元胞的标准差计算结果；hk为第k个元胞的水位值，h为中间元胞与其八个邻域元胞的水位值均值，N＝9，表示八个邻域和一个中间元胞；dCA为元胞尺寸，m。
(3-2)中间元胞水量分配比例计算：计算中间元胞对下游元胞的水量分配比例di。
其中，di为中间元胞对第i个下游邻域元胞的水量分配比例；tanβi为中间元胞到第i个下游邻域元胞的坡度比降；Li为对第i个下游邻域元胞的等高线长度加权因子，本实施例中，垂直方向元胞的等高线长度加权权重为0.500，对角方向元胞的加权权重为0.354(即为的近似值)；M为下游邻域元胞的数量，p为水量分配权重。
(4)比较模型时间步长与水量实际转移时间的大小，确定一个模型计算时间步长内的最终水量分配比例。
该步骤具体包括：
对比模型时间步长t与(3)中得到的水量转移至各个下游元胞中的时间Ti的大小，若Ti大于时间步长t，只有部分水量可以被转移，此时向下游转移水量比例为如果Ti小于时间步长t，认为水量可以全部转移，此时向下游元胞转移的水量比例为di。
(5)计算中间元胞的转移水量，并基于“动力系统倾向于往平衡状态发展”假设，设置各个下游方向转移上限，修正转移至各个下游元胞的水量。
该步骤具体包括：
(5-1)转移水量计算：根据中间元胞水量(面积×水深)与中间元胞水量转移比例，乘积即可获得中间元胞向各个下游元胞转移的水量。
Pi＝depth0×A×di
其中，depth0为中间元胞水深，m3；A表示元胞面积，也即dCA×dCA，m2；Pi为中间元胞向各个下游元胞转移的水量。
(5-2)转移水量修正：结合“动力系统倾向于往平衡状态发展”假设，对从中间元胞转移水量至第i个下游元胞设置上限，进行转移水量的修正。
其中，hi表示第i个下游元胞水位值，m；表示中间元胞与第i个下游元胞的平均水位值，m；A表示元胞面积，也即dCA×dCA，m2。
基于相同的发明构思，本发明实施例提供的一种基于元胞自动机的城市地表径流流向计算装置，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述计算机程序被加载至处理器时实现所述的基于元胞自动机的城市地表径流流向计算方法。
以上所揭露的仅为本发明一种较佳实施例而已，不能以此来限定本发明之权利范围，因此依本发明权利要求所作的等同变化，仍属本发明所涵盖的范围。

Claims (7)

1.一种基于元胞自动机的城市地表径流流向计算方法，其特征在于，包括如下步骤：
(1)建立研究区的基于元胞自动机模式的基础地理场景，其中元胞的邻域关系为Moore型八邻域；
(2)比较中间元胞与邻域元胞的水位值大小，排除水位值高于中间元胞的邻域元胞，将剩下的邻域元胞作为水量转移的下游元胞，计算中间元胞水量转移到各个下游元胞所需时间；
(3)计算中间元胞与其八个邻域元胞之间水位值的标准差及标准差阈值，对比计算标准差和标准差阈值相对大小，选择水量分配权重值，确定各个下游元胞中的水量分配比例；
(4)比较模型时间步长与水量实际转移时间的大小，确定一个模型计算时间步长内的水量分配比例；
(5)根据中间元胞水量与水量分配比例计算中间元胞向各个下游元胞转移的水量，并设置各个下游方向转移上限，修正转移至各个下游元胞的水量。
2.根据权利要求1所述的基于元胞自动机的城市地表径流流向计算方法，其特征在于，步骤(2)中，根据如下公式计算中间元胞水量转移到各个下游元胞所需时间：
其中，Ti为中间元胞水量转移到第i个下游元胞所需时间，depth0为中间元胞水深，si为中间元胞与第i个下游元胞之间的水面坡度，ni为曼宁糙率系数，Di为中间元胞与第i个下游元胞之间的距离。
3.根据权利要求1所述的基于元胞自动机的城市地表径流流向计算方法，其特征在于，步骤(3)中，根据如下公式确定各个下游元胞中的水量分配比例：
其中，di为中间元胞对第i个下游邻域元胞的水量分配比例，tanβi为中间元胞到第i个下游邻域元胞的坡度比降，Li为对第i个下游邻域元胞的等高线长度加权因子，M为下游邻域元胞的数量，p为水量分配权重，SD和SD0为中间元胞与其八个邻域元胞之间水位值的标准差及标准差阈值。
4.根据权利要求3所述的基于元胞自动机的城市地表径流流向计算方法，其特征在于，中间元胞与其八个邻域元胞之间水位值的标准差SD及标准差阈值SD0为：
其中，hk为第k个元胞的水位值，为中间元胞与其八个邻域元胞的水位值均值，N＝9，表示八个邻域和一个中间元胞，dCA为元胞边长。
5.根据权利要求1所述的基于元胞自动机的城市地表径流流向计算方法，其特征在于，步骤(4)中，对比模型时间步长t与步骤(3)中得到的水量转移至各个下游元胞中的时间Ti的大小，若Ti大于时间步长t，只有部分水量可以被转移，此时向下游转移水量比例为如果Ti小于时间步长t，认为水量可以全部转移，此时向下游元胞转移的水量比例为di。
6.根据权利要求1所述的基于元胞自动机的城市地表径流流向计算方法，其特征在于，步骤(5)中，根据如下公式对中间元胞向各个下游元胞转移的水量进行修正：
其中，Pi为中间元胞向第i个下游元胞转移的水量，hi为第i个下游元胞水位值，为中间元胞与第i个下游元胞的平均水位值，A为元胞面积。
7.一种基于元胞自动机的城市地表径流流向计算装置，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述计算机程序被加载至处理器时实现根据权利要求1-6任一项所述的基于元胞自动机的城市地表径流流向计算方法。
CN201910737873.3A 2019-08-12 2019-08-12 一种基于元胞自动机的城市地表径流流向计算方法及装置 Active CN110442988B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910737873.3A CN110442988B (zh) 2019-08-12 2019-08-12 一种基于元胞自动机的城市地表径流流向计算方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910737873.3A CN110442988B (zh) 2019-08-12 2019-08-12 一种基于元胞自动机的城市地表径流流向计算方法及装置

Publications (2)

Publication Number Publication Date
CN110442988A true CN110442988A (zh) 2019-11-12
CN110442988B CN110442988B (zh) 2020-07-28

Family

ID=68434436

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910737873.3A Active CN110442988B (zh) 2019-08-12 2019-08-12 一种基于元胞自动机的城市地表径流流向计算方法及装置

Country Status (1)

Country Link
CN (1) CN110442988B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111931440A (zh) * 2020-08-28 2020-11-13 中国水利水电科学研究院 一种考虑地表积水深变化的坡面多流向漫流分析方法
CN112052545A (zh) * 2020-08-25 2020-12-08 水利部交通运输部国家能源局南京水利科学研究院 一种基于元胞自动机的城市地表径流与管网汇流耦合方法
CN112270115A (zh) * 2020-11-25 2021-01-26 同济大学 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102008063453A1 (de) * 2008-12-17 2010-07-15 Siemens Aktiengesellschaft Verfahren zur Verbesserung der Simulation von Objektströmen mittels Flutungsalgorithmus
CN102708227A (zh) * 2012-04-18 2012-10-03 国家测绘局卫星测绘应用中心 基于sph算法的洪水溃坝过程仿真方法与仿真系统
CN103544390A (zh) * 2013-10-18 2014-01-29 西南交通大学 基于元胞自动机的溃决洪水快速演进模拟分析方法
US20170103172A1 (en) * 2015-10-07 2017-04-13 The Arizona Board Of Regents On Behalf Of The University Of Arizona System And Method To Geospatially And Temporally Predict A Propagation Event
CN109190263A (zh) * 2018-09-10 2019-01-11 柳创新 基于全流域降雨径流及水动力模型预测降水流量的方法
CN109190168A (zh) * 2018-07-31 2019-01-11 中国科学院遥感与数字地球研究所 一种实时测量数据支持下的洪水演进模拟动态修正方法
CN109583035A (zh) * 2018-11-05 2019-04-05 水利部交通运输部国家能源局南京水利科学研究院 基于元胞自动机的城市面源污染物累积冲刷过程计算方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102008063453A1 (de) * 2008-12-17 2010-07-15 Siemens Aktiengesellschaft Verfahren zur Verbesserung der Simulation von Objektströmen mittels Flutungsalgorithmus
CN102708227A (zh) * 2012-04-18 2012-10-03 国家测绘局卫星测绘应用中心 基于sph算法的洪水溃坝过程仿真方法与仿真系统
CN103544390A (zh) * 2013-10-18 2014-01-29 西南交通大学 基于元胞自动机的溃决洪水快速演进模拟分析方法
US20170103172A1 (en) * 2015-10-07 2017-04-13 The Arizona Board Of Regents On Behalf Of The University Of Arizona System And Method To Geospatially And Temporally Predict A Propagation Event
CN109190168A (zh) * 2018-07-31 2019-01-11 中国科学院遥感与数字地球研究所 一种实时测量数据支持下的洪水演进模拟动态修正方法
CN109190263A (zh) * 2018-09-10 2019-01-11 柳创新 基于全流域降雨径流及水动力模型预测降水流量的方法
CN109583035A (zh) * 2018-11-05 2019-04-05 水利部交通运输部国家能源局南京水利科学研究院 基于元胞自动机的城市面源污染物累积冲刷过程计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
QI SHAO: "RunCA: A cellular automata model for simulating surface runoff at different scales", 《JOURNAL OF HYDROLOGY》 *
何梦男: "基于元胞自动机的应急疏散最短路径优化算法", 《中国安全科学学报》 *
赖泽辉: "基于元胞自动机的城市地表积水模拟研究", 《水土保持通报》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112052545A (zh) * 2020-08-25 2020-12-08 水利部交通运输部国家能源局南京水利科学研究院 一种基于元胞自动机的城市地表径流与管网汇流耦合方法
CN112052545B (zh) * 2020-08-25 2021-10-08 水利部交通运输部国家能源局南京水利科学研究院 一种基于元胞自动机的城市地表径流与管网汇流耦合方法
CN111931440A (zh) * 2020-08-28 2020-11-13 中国水利水电科学研究院 一种考虑地表积水深变化的坡面多流向漫流分析方法
CN111931440B (zh) * 2020-08-28 2021-03-05 中国水利水电科学研究院 一种考虑地表积水深变化的坡面多流向漫流分析方法
CN112270115A (zh) * 2020-11-25 2021-01-26 同济大学 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法
CN112270115B (zh) * 2020-11-25 2022-06-14 同济大学 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法

Also Published As

Publication number Publication date
CN110442988B (zh) 2020-07-28

Similar Documents

Publication Publication Date Title
CN108108911B (zh) 基于水生态系统健康的大型湖泊分区水质目标制定方法
CN110442988A (zh) 一种基于元胞自动机的城市地表径流流向计算方法及装置
Fang et al. Modelling blowing snow redistribution to prairie wetlands
Chen et al. Load estimation and source apportionment of nonpoint source nitrogen and phosphorus based on integrated application of SLURP model, ECM, and RUSLE: a case study in the Jinjiang River, China
CN107729695A (zh) 一种用于小流域次洪模拟的水文模型率定方法
CN108388714A (zh) 流域水系和城市管网耦合的平原河网城市洪水模拟方法
Elias et al. Hydrodynamic validation of Delft3D with field measurements at Egmond
CN108256177A (zh) 一种河流水沙模型的参数优化方法和系统
CN112270115B (zh) 一种基于元胞自动机的复杂地形洪水淹没进程模拟方法
CN112784502B (zh) 一种水文水力学动态双向耦合的洪水预测方法
Hu A review of the models for Lake Taihu and their application in lake environmental management
CN106875060A (zh) 一种基于全局自动优化算法的洪水实时校正方法
CN115375198B (zh) 区域河湖水系连通联合调度和水质安全保障方法及系统
CN110187413B (zh) 一种城市内涝预报方法、电子设备及存储介质
CN113011685A (zh) 一种无径流资料地区内陆湖泊水位变化模拟预测方法
CN107609713A (zh) 一种通过降雨和径流双要素校正的陆气耦合实时预报方法
Grabs et al. Discharge observation networks in Arctic regions: Computation of the river runoff into the Arctic Ocean, its seasonality and variability
CN108399577A (zh) 一种基于蒸散发的林地植被生态需水的定量测算方法
Chea et al. Flow simulation in an ungauged catchment of Tonle Sap Lake Basin in Cambodia: Application of the HEC-HMS model
CN109583035A (zh) 基于元胞自动机的城市面源污染物累积冲刷过程计算方法
CN103577895B (zh) 一种资料短缺情形下的二次耦合月径流预报方法
CN106320255B (zh) 一种漫滩水流滩槽过流量的计算方法
CN111105061A (zh) 一种河流水质的预测方法、预测装置及终端
CN110287504B (zh) 一种模拟近岸海水交换规律的耦合模型的建立方法
CN106096129B (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