分布式水文模型并行计算的研究进展.pdf

第32 卷 第4 期 地 理 科 学 进 展 Vol.32, No.4 2013 年04 月 PROGRESSINGEOGRAPHY Apr., 2013 分布式水文模型的并行计算研究进展 1,2 1,3 1 4 1,2 1,2 刘军志 ，朱阿兴 ，秦承志 ，陈腊娇 ，吴 辉 ，江净超 (1. 中国科学院地理科学与资源研究所资源与环境信息系统国家重点实验室，北京 100101；2. 中国科学院大学，北京 100049；3. 威斯康星大学麦迪逊分校地理系，麦迪逊53706 ，美国；4. 中国科学院遥感与数字地球研究所，北京 100094) 摘 要：大流域、高分辨率、多过程耦合的分布式水文模拟计算量巨大，传统串行计算技术不能满足其对计算能力 的需求，因此需要借助于并行计算的支持。本文首先从空间、时间和子过程三个角度对分布式水文模型的可并行 性进行了分析，指出空间分解的方式是分布式水文模型并行计算的首选方式，并从空间分解的角度对水文子过程 计算方法和分布式水文模型进行了分类。然后对分布式水文模型的并行计算研究现状进行了总结。其中，在空 间分解方式的并行计算方面，现有研究大多以子流域作为并行计算的基本调度单元；在时间角度的并行计算方 面，有学者对时空域双重离散的并行计算方法进行了初步研究。最后，从并行算法设计、流域系统综合模拟的并 行计算框架和支持并行计算的高性能数据读写方法3 个方面讨论了当前存在的关键问题和未来的发展方向。 关 键 词：分布式水文模型；并行计算；研究进展；空间分解 doi: 10.11820/dlkxjz.2013.04.006 1 引言 随着并行计算软硬件技术的发展，并行计算的 门槛不断降低。在硬件方面，多核处理器、图形处 分布式水文模型考虑了气候和下垫面因子的 理器(Graphic Processing Unit, GPU)和计算机集群 空间异质性，与集总式模型相比能更准确地模拟流 等并行计算设备已成为普通用户容易获取的计算 域水文过程，从而能更好地预测自然和人为因素对 工具(陈国良等, 2009; 卢风顺等, 2011) 。在软件方 流域水资源和水环境的影响，为流域科学管理提供 面，MPI(Message Passing Interface) 、OpenMP(Open 决策依据(Borah et al, 2004; 王中根等, 2002) 。随着 Multi-Processing) 、CUDA(Compute Unified Device 地理信息系统和遥感技术的发展，分布式水文模型 Architecture) 、OpenCL(Open Computing Language) 所需的地形、植被、土壤等空间数据的获取日益方 等并行编程标准和编程库降低了并行编程的难度 便，并且分辨率越来越高，促使分布式水文模拟的 (Rauber et al, 2010; 贾海鹏等, 2012) ，提高了并行程 空间范围越来越大，所采用的时空分辨率越来越精 序的可移植性。在信息基础设施的利用方面，网格 细(徐宗学等, 2010) 。与此同时，模拟涉及的流域过 计算技术在科学研究中得到越来越广泛的应用 程也越来越多，从传统的单一水文过程模拟逐渐扩 (Wang et al, 2009; Zhao et al, 2013) 。这些并行计算 展到生态、土壤侵蚀和非点源污染等多过程耦合的 技术的发展为开发并行化的分布式水文模型，解决 流域系统综合模拟( 陈腊娇等, 2011; 袁再健等, 大范围、高分辨率分布式水文模拟的性能瓶颈问题 2006; 朱阿兴等, 2012) ，其计算过程更趋复杂。分 提供了良好的条件。 布式水文模拟的这些新的