分布式水文模型的子流域-基本单元双层并行计算方法.pdfVIP

致谢 博士四年，最要感谢的是我的导师朱阿兴研究员！朱老师学识渊博，思维活 跃，极具学者风度，在学术和生活上都是学生的榜样。朱老师对我们既严格要求， 又给予充分自由的发挥空间；传授给我们更多的是方法，而不仅仅是知识。跟朱 老师学习这四年，我受益匪浅。 感谢加拿大圭尔夫大学的刘永波研究员。我从博士阶段才开始学习水文模型， 刘老师是我这方面的启蒙老师。每次我向刘老师请教问题，总能得到刘老师的耐 心讲解，我非常感激。 感谢中国科学院大学的赵英时教授和宋现锋教授。两位老师从我硕士期间就 一直给予我关心和帮助，我一定继续努力，不辜负你们的教诲。 感谢精细地理过程模拟研究组的李宝林、裴韬、秦承志和杨琳老师四年来对 我的关心和帮助。 感谢国防科技大学的景宁、陈荦、吴秋云、刘露老师和杨岸然博士在 863 项目集成时给予的帮助。 感谢马志鹏、姚一鸣、王雪梅和王淑强老师为我们提供了良好的科研环境。 感谢2415 办公室的兄弟姐妹们。你们让简单重复的科研生活变得不那么枯 燥。 感谢二十多年来一直为我默默付出的家人，你们的亲情和关心是我不断前进 的动力。 最后，感谢百忙之中抽出宝贵时间审阅本文的各位评委老师！ 摘 要 大流域、高分辨率、多过程耦合的分布式水文模拟计算量巨大，传统串行计算技术不能 满足其对计算能力的需求，迫切需要并行计算的支持。当前分布式水文模型的并行计算研究 大多以子流域作为调度单元。基于子流域进行并行计算优点是：子流域之间只存在松散耦合 关系，通信开销小。但是，很多分布式水文模型的基本模拟单元往往采用比子流域更精细的 单元（如栅格），这些模型除了在子流域层次具有可并行性，在基本单元层次也往往具有很 好的可并行性。基于子流域的并行计算方法没有利用这些较细层次上的可并行性，因而不能 充分发挥这些模型的并行潜力。 针对上述问题，本文提出一种能同时利用子流域层次和更细的基本单元层次上可并行性 的分布式水文模型并行计算方法——子流域-基本单元双层并行计算方法，并以作者开发的 能体现水流空间运动的分布式水文模型（Spatially Explicit Integrated Modeling, SEIM）为例， 对该并行方法进行了研究。 子流域-基本单元双层并行计算方法由两种相互独立且互补的并行计算方法组成，分别 是： (1) 子流域层次消息传递方式的并行计算方法 在SEIM 模型中，子流域只在出口与相邻子流域存在数据依赖关系，子流域层次的并行 计算任务间通信量小，因此采用消息传递的方式在计算机集群平台上进行并行计算。对于该 层次的并行计算，各个子流域的坡面过程计算不存在相互依赖关系，具有很高的可并行性； 而河道汇流过程的计算必须按照从上游到下游的顺序依次进行，不存在上下游关系的子流域 可以并行执行，存在上下游关系的子流域只能串行进行，可并行性较差。 针对这一特点，本文在借鉴前人研究成果的基础上，提出了一种坡面和河道过程分离的 主从式并行计算方法，该方法中包括一个管理进程、多个计算进程和一个水文数据库，任务 分配采用METIS 图分解算法。 (2) 基本单元层次共享内存方式的并行计算方法 在基本单元层次，SEIM 模型中的坡面汇流等水文子过程的并行计算需要较大的通信量， 因此采用通信代价小的共享内存并行模型。根据计算特点，可以将基本单元层次的子过程模 拟方法分为计算独立型和顺序依赖型两类。对于计算独立型方法，采用简单分块的方式进行 并行计算；对于顺序依赖型方法，本文提出一种基于模拟单元分层的并行计算方法。 在上述两种方法中，子流域层次的并行计算可运行于计算机集群环境，在计算架构上具 有很好的可伸缩性，但是没有利用子流域内部更细粒度的可并行性，不能充分发挥分布式水 文模型的并行潜力。而基本单元层次的并行计算充分利用了基本单元层次的可并行性，但是 依赖于共享内存硬件平台，可伸缩性较差。针对两种方法的优缺点，本文将上述两种方法进 行组合，得到子流域-基本单元双层并行计算方法，以充分利用两种方法的优点提高并行计 算性能。 i 对两种方法进行组合的思路为：以多核集群为运行平台，首先将子流域的计算任务分配 到集群的多个CPU （Central Processing Unit ）进行消息传递