CN104142812A - 一种分布式水文模型并行运算方法 - Google Patents
一种分布式水文模型并行运算方法 Download PDFInfo
- Publication number
- CN104142812A CN104142812A CN201410371523.7A CN201410371523A CN104142812A CN 104142812 A CN104142812 A CN 104142812A CN 201410371523 A CN201410371523 A CN 201410371523A CN 104142812 A CN104142812 A CN 104142812A
- Authority
- CN
- China
- Prior art keywords
- simulation
- sub
- streams
- territory
- node
- 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
Links
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明提供了一种分布式水文模型并行运算方法,通过确定循环时段初始状态,主要包括上游依赖子流域数目和各子流域至流域出口线路上所有子流域所需模拟时间的累积和;选择累积运行时间最长且上游依赖河段数等于0的河段,建立节点子进程进行并行模拟;通过节点模拟并行子进程完成整个子流域产汇流运算模拟,结束后将其直接下游子流域的上游依赖子流域数目减1;循环执行动态分配过程和节点模拟过程,直到所有子流域均完成模拟,然后进入下一循环时段进行模拟。本发明提高了分布式水文模型并行效率,能够确保消耗时间最长的子流域优先模拟。
Description
技术领域
本发明属于水文模型领域,尤其涉及一种分布式水文模型并行运算方法。
背景技术
随着分布式水文模型应用尺度的增大,基本计算单元个数越来越多,模型运算所需运行时间也越来越多。多线程并行运算成为提高模型运算速率的重要方法。就模型并行运算方面,许多学者进行了研究,并提出不同的并行运算方法,涉及静态/动态并行任务分配方案、时间/空间并行分配方案等。现有的并行运算方法假设每个子流域具有基本相同的运行时间,对不同子流域具有不同的运行时间情况的优化效率不高,在遇到变时间步长问题时,即不同的子流域之间以及在不同的模拟循环中都具有不同的模拟运行时间,动态分配效率低。
发明内容
本发明的目的在于提供一种分布式水文模型并行运算方法,旨在解决现有的并行运算动态分配技术只能处理同一个模型单个模拟循环中所有子流域以及同一子流域在不同模拟循环中都具有相同的模拟运算时间的情况。
本发明是这样实现的,一种分布式水文模型并行运算方法,包括以下步骤:
S1、确定循环时段初始状态,主要包括上游依赖子流域数目(初始值等于上游子流域数目)和各子流域至流域出口线路上所有子流域所需模拟时间的累积和;
S2、选择累积运行时间最长且上游依赖河段数等于0的河段,建立节点子进程进行并行模拟,即动态分配过程;
S3、通过节点模拟并行子进程完成整个子流域产汇流运算模拟,结束后将其直接下游子流域的上游依赖子流域数目减1,即节点模拟过程;
S4、循环执行动态分配过程和节点模拟过程,直到所有子流域均完成模拟,然后进入下一循环时段进行模拟。
相比于现有技术的缺点和不足,本发明具有以下有益效果:
(1)本发明能够适用于对子流域采用变时间步长模拟的分布式水文模型并行运算;本发明充分考虑了不同模拟循环中各子流域具有不同的模拟时间的情景进行动态分配,提高分布式水文模型并行效率。
(2)本发明根据不同模拟时段各子流域实际所需运行时间进行动态分配,能够确保累积消耗时间最长的子流域优先模拟。
附图说明
图1是本发明分布式水文模型并行运算方法的步骤流程图;
图2是本发明实施例中子流域示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
一种分布式水文模型并行运算方法,如图1所示,包括以下步骤:
S1、确定循环时段初始状态,主要包括上游依赖子流域数目(初始值等于上游子流域数目)和各子流域至流域出口线路上所有子流域所需模拟时间的累积和
在步骤S1中,首先根据子流域内河网拓扑关系,获取子流域上下游拓扑关系信息,如图2和表1所示,其中,图2为子流域示意图,表1为子流域上下游拓扑信息(其中子流域编号0表示没有)。
表1子流域上下游拓扑信息
这部分信息属于基础信息,由外部工具提供,只需计算一次即可。该步骤主要是将上游依赖子流域数目信息进行初始化,初始值等于子流域上下游拓扑关系信息中的上游子流域数目。此外,根据当此模拟循环的具体情况(例如在日尺度模拟情况下,有的子流域按日模拟需要1个单位的基本时间,而有的子流域按小时模拟需要24个单位的基本时间),确定所有子流域所需的模拟时间后,根据子流域上下游拓扑关系计算各子流域所需的累积模拟时间,例如在图2中,1号子流域的累积模拟时间等于1号、5号、6号这三个子流域模拟时间之和。当计算出所有子流域累积模拟时间后,按值从大到小进行排列,即累积模拟时间越大的越靠前。至此,完成当前循环的初始化过程。
S2、选择累积运行时间最长且上游依赖河段数等于0的河段,建立节点子进程进行并行模拟,即动态分配过程
在步骤S2中,动态任务分配主线程采用循环的结构进行检测当前模拟循环过程是否所有子流域均已经完成,如果完成则进入步骤S4;否则检测是否有空闲节点进程,如果没有则继续循环检测;如果有空闲节点进程,则进行动态节点并行任务分配。根据排序后的子流域序列从累积模拟时间最大的子流域开始依次往下,找到第一个上游依赖河段数等于0的子流域节点编号,建立节点并行运算进程,然后主进程继续循环检测。
S3、通过节点模拟并行子进程完成整个子流域产汇流运算模拟,结束后将其直接下游子流域的上游依赖子流域数目减1,即节点模拟过程
在步骤S3中,主要进行子流域产汇流模拟运算工作,相关工作可有具体的水文模型进行工作。节点模拟过程的关键在于,当子流域产汇流模拟结束后,需要将其下游子流域的上游依赖子流域数目减1。例如,在图2中,1号子流域节点进行运行完毕,则需要将5号子流域的上游依赖子流域数目减1(假设2、3号子流域均未模拟,则该步骤就是将5号子流域的上游依赖子流域数由3变成2)。接着2号子流域模拟完成,则5号子流域上游依赖子流域则由2变成1,依次直到等于0,从而可以被主进行选中进程任务分配。
S4、循环执行动态分配过程和节点模拟过程,直到所有子流域均完成模拟,然后进入下一循环时段进行模拟
在步骤S4中,循环执行动态分配过程和节点模拟过程,直到所有子流域均完成模拟。如果所有模拟时段均已完成,则程序介绍,否则进入步骤S1开始下一个模拟循环过程。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (1)
1.一种分布式水文模型并行运算方法,其特征在于包括以下步骤:
S1、确定循环时段初始状态,主要包括上游依赖子流域数目和各子流域至流域出口线路上所有子流域所需模拟时间的累积和;
S2、选择累积运行时间最长且上游依赖河段数等于0的河段,建立节点子进程进行并行模拟,即动态分配过程;
S3、通过节点模拟并行子进程完成整个子流域产汇流运算模拟,结束后将其直接下游子流域的上游依赖子流域数目减1,即节点模拟过程;
S4、循环执行动态分配过程和节点模拟过程,直到所有子流域均完成模拟,然后进入下一循环时段进行模拟。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410371523.7A CN104142812B (zh) | 2014-07-30 | 2014-07-30 | 一种分布式水文模型并行运算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410371523.7A CN104142812B (zh) | 2014-07-30 | 2014-07-30 | 一种分布式水文模型并行运算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104142812A true CN104142812A (zh) | 2014-11-12 |
CN104142812B CN104142812B (zh) | 2017-02-01 |
Family
ID=51851996
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410371523.7A Expired - Fee Related CN104142812B (zh) | 2014-07-30 | 2014-07-30 | 一种分布式水文模型并行运算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104142812B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105160121A (zh) * | 2015-09-17 | 2015-12-16 | 天津市水文水资源勘测管理中心 | 一种有限元控制的分布式水文模型的建模方法 |
CN106897529A (zh) * | 2017-03-02 | 2017-06-27 | 中国水利水电科学研究院 | 基于污水管网拓扑关系的子流域汇流演算顺序计算方法 |
CN109190160A (zh) * | 2018-07-27 | 2019-01-11 | 华中科技大学 | 一种分布式水文模型的矩阵化模拟方法 |
CN109753362A (zh) * | 2019-01-14 | 2019-05-14 | 中国水利水电科学研究院 | 一种分布式水文模型的汇流并行调度方法 |
CN109902366A (zh) * | 2019-02-14 | 2019-06-18 | 中国水利水电科学研究院 | 一种分布式水文模型汇流并行方法 |
CN109918741A (zh) * | 2019-02-13 | 2019-06-21 | 北京科技大学 | 一种适用于大规模水文模拟参数率定方法及率定装置 |
CN114707325A (zh) * | 2022-03-30 | 2022-07-05 | 洛阳师范学院 | 一种基于物理模型的水文过程串并耦合动态模拟方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103092572A (zh) * | 2013-01-11 | 2013-05-08 | 中国科学院地理科学与资源研究所 | 一种集群环境下分布式水文模拟的并行化方法 |
CN103164190A (zh) * | 2013-03-02 | 2013-06-19 | 中国科学院对地观测与数字地球科学中心 | 一种全分布式流域生态水文模型的快速并行化方法 |
-
2014
- 2014-07-30 CN CN201410371523.7A patent/CN104142812B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103092572A (zh) * | 2013-01-11 | 2013-05-08 | 中国科学院地理科学与资源研究所 | 一种集群环境下分布式水文模拟的并行化方法 |
CN103164190A (zh) * | 2013-03-02 | 2013-06-19 | 中国科学院对地观测与数字地球科学中心 | 一种全分布式流域生态水文模型的快速并行化方法 |
Non-Patent Citations (3)
Title |
---|
HAO WANG: ""Maximum speedup ratio curve(MSC) in parallelcomputing of the binary-tree-based drainage network"", 《COMPUTERS & GEOSCIENCES》 * |
TIEJIAN LI: ""Dynamic parallelization of hydrological model simulations"", 《ENVIRONMENTAL MODELLING & SOFTWARE》 * |
刘军志: ""基于栅格分层的逐栅格汇流算法并行化研究"", 《国防科技大学学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105160121A (zh) * | 2015-09-17 | 2015-12-16 | 天津市水文水资源勘测管理中心 | 一种有限元控制的分布式水文模型的建模方法 |
CN105160121B (zh) * | 2015-09-17 | 2018-07-13 | 天津市水文水资源勘测管理中心 | 一种有限元控制的分布式水文模型的建模方法 |
CN106897529A (zh) * | 2017-03-02 | 2017-06-27 | 中国水利水电科学研究院 | 基于污水管网拓扑关系的子流域汇流演算顺序计算方法 |
CN106897529B (zh) * | 2017-03-02 | 2020-03-31 | 中国水利水电科学研究院 | 基于污水管网拓扑关系的子流域汇流演算顺序计算方法 |
CN109190160A (zh) * | 2018-07-27 | 2019-01-11 | 华中科技大学 | 一种分布式水文模型的矩阵化模拟方法 |
CN109753362A (zh) * | 2019-01-14 | 2019-05-14 | 中国水利水电科学研究院 | 一种分布式水文模型的汇流并行调度方法 |
CN109753362B (zh) * | 2019-01-14 | 2020-03-24 | 中国水利水电科学研究院 | 一种分布式水文模型的汇流并行调度方法 |
CN109918741A (zh) * | 2019-02-13 | 2019-06-21 | 北京科技大学 | 一种适用于大规模水文模拟参数率定方法及率定装置 |
CN109918741B (zh) * | 2019-02-13 | 2020-11-27 | 北京科技大学 | 一种适用于大规模水文模拟的参数率定方法 |
CN109902366A (zh) * | 2019-02-14 | 2019-06-18 | 中国水利水电科学研究院 | 一种分布式水文模型汇流并行方法 |
CN109902366B (zh) * | 2019-02-14 | 2020-05-05 | 中国水利水电科学研究院 | 一种分布式水文模型汇流并行方法 |
CN114707325A (zh) * | 2022-03-30 | 2022-07-05 | 洛阳师范学院 | 一种基于物理模型的水文过程串并耦合动态模拟方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104142812B (zh) | 2017-02-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104142812A (zh) | 一种分布式水文模型并行运算方法 | |
CN106201718A (zh) | 一种基于负载预测的云计算资源动态伸缩方法 | |
Bhuiyan et al. | Energy-efficient parallel real-time scheduling on clustered multi-core | |
CN105589736B (zh) | 基于网表分割和多线程并行的硬件描述语言仿真加速方法 | |
GB201013074D0 (en) | Executing a utility in a distributed computing system based on an integrated model | |
Campbell et al. | Water supply network sectorization based on social networks community detection algorithms | |
CN103294599A (zh) | 一种基于云的嵌入式软件交叉测试方法 | |
TWI507989B (zh) | 資源導向之嵌入式系統功率消耗分析方法 | |
GB2516198A (en) | Business process analytics | |
Baki et al. | Tools for energy footprint assessment in urban water systems | |
Charr et al. | Dynamic frequency scaling for energy consumption reduction in synchronous distributed applications | |
Liu et al. | Estimation of theoretical maximum speedup ratio for parallel computing of grid-based distributed hydrological models | |
CN110059094A (zh) | 主动配电网分布式状态估计的分区方法、介质及设备 | |
Hosseinabady et al. | Energy optimization of FPGA-based stream-oriented computing with power gating | |
CN106201655B (zh) | 虚拟机分配方法和虚拟机分配系统 | |
CN103413196B (zh) | 改进的基于关系链的作业车间调度关键外协工序辨识方法 | |
Uddin et al. | Signature-based high-level simulation of microthreaded many-core architectures | |
CN103218249B (zh) | 一种虚拟化集群控制方法及设备、虚拟化集群系统 | |
CN105956363A (zh) | 一种河流生态需水分段分析、补充以及分流方法 | |
CN109902366A (zh) | 一种分布式水文模型汇流并行方法 | |
Jo et al. | Development of the non-revenue water analysis equation through the statistical analysis of main parameter in waterworks system in Incheon City | |
CN102855339A (zh) | 集成电路版图验证并行处理解决方案 | |
Karimzadeh et al. | Developing a Diurnal Pattern of Sewage Flow at Inlet Nodes Using Limited Measured Nodes–A Case Study of Greater Tehran Sewer System | |
Wang et al. | Composite dispatching rule design for photolithography area scheduling in wafer manufacturing system with multiple objectives | |
Shaker Ardakani et al. | Estimating The Energy Demand of Natural Gas, Electricity and Petroleum Products in Manufacturing Industries Based on SURE Method |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170201 Termination date: 20170730 |