基于栅格分层的逐栅格汇流算法并行化研究.pdfVIP

第３５卷 第１期 国　防　科　技　大　学　学　报 Ｖｏｌ．３５Ｎｏ．１ ２０１３年２月 ＪＯＵＲＮＡＬＯＦＮＡＴＩＯＮＡＬＵＮＩＶＥＲＳＩＴＹＯＦＤＥＦＥＮＳＥＴＥＣＨＮＯＬＯＧＹ Ｆｅｂ．２０１３ 基于栅格分层的逐栅格汇流算法并行化研究 １，２ １，３ ４ １ １，２ １，２ １ 刘军志 ，朱阿兴 ，刘永波 ，秦承志 ，陈腊娇 ，吴　辉 ，杨　琳 （１．中国科学院 资源与环境信息系统国家重点实验室，北京　１００１０１； ２．中国科学院大学，北京　１０００４９； ３．ＤｅｐａｒｔｍｅｎｔｏｆＧｅｏｇｒａｐｈｙ，ＵｎｉｖｅｒｓｉｔｙｏｆＷｉｓｃｏｎｓｉｎＭａｄｉｓｏｎ，Ｍａｄｉｓｏｎ，ＷＩ５３７０６，ＵＳＡ； ４．ＤｅｐａｒｔｍｅｎｔｏｆＧｅｏｇｒａｐｈｙ，ＵｎｉｖｅｒｓｉｔｙｏｆＧｕｅｌｐｈ，５０ＳｔｏｎｅＲｏａｄＥａｓｔ，Ｇｕｅｌｐｈ，ＯｎｔａｒｉｏＮ１Ｇ２Ｗ１，Ｃａｎａｄａ） 摘　要：分布式水文模型中的逐栅格汇流算法计算量大，需要借助并行计算以满足大流域长历时模拟的 要求。针对目前鲜有对基于隐式有限差分的逐栅格汇流算法进行并行计算研究的情况，基于栅格分层的思 想提出一种适用于共享内存并行计算环境的逐栅格汇流并行算法。该算法首先根据流向进行栅格分层，使 同一层中栅格的计算相互独立，然后将同一层中栅格的计算任务分配到多个计算单元并行计算。采用 Ｃ＋＋ 编程语言与ＯｐｅｎＭＰ并行编程库实现了该算法，并选择河北省清水河流域为实验区，在不同数据规模（３０ｍ、 ９０ｍ、２７０ｍ分辨率）、不同核数（２～２０个）以及不同栅格分层方法的情况下对算法性能进行了测试。实验结 果表明本文提出的并行算法具有较好的加速比和并行效率，且并行效率随数据规模的增大而增大。栅格分 层方法对算法并行性能有明显影响，从上到下的分层方法比从下到上的方法具有更高的并行效率。 关键词：逐栅格汇流；隐式有限差分；栅格分层；并行计算；ＯｐｅｎＭＰ 中图分类号：Ｐ３３３．９　　文献标志码：Ａ　　 文章编号：１００１－２４８６（２０１３）０１－０１２３－０７ Ｐａｒａｌｌｅｌｉｚａｔｉｏｎｏｆａｇｒｉｄｔｏｇｒｉｄｒｏｕｔｉｎｇａｌｇｏｒｉｔｈｍｂａｓｅｄ ｏｎｇｒｉｄｓｌａｙｅｒｉｎｇ １，２ １，３ ４ １ １，２ １，２ １ ＬＩＵＪｕｎｚｈｉ ，ＺＨＵＡＸｉｎｇ ，ＬＩＵＹｏｎｇｂｏ，ＱＩＮＣｈｅｎｇｚｈｉ，ＣＨＥＮＬａｊｉａｏ ，ＷＵＨｕｉ ，ＹＡＮＧＬｉｎ （１．ＳｔａｔｅＫｅｙＬａｂｏｆＲｅｓｏｕｒｃｅｓａｎｄＥｎｖｉｒｏｎｍｅｎｔａｌＩｎｆｏｒｍａｔｉｏｎＳｙｓｔｅｍ， ＩｎｓｔｉｔｕｔｅｏｆＧｅｏｇｒａｐｈｉｃＳｃｉｅｎｃｅｓａｎｄＮａｔｕｒａｌＲｅｓｏｕｒｃｅｓＲｅｓｅａｒｃｈ，ＣＡＳ，Ｂｅｉｊｉｎｇ１００１０１，Ｃｈｉｎａ； ２．ＵｎｉｖｅｒｓｉｔｙｏｆＣｈｉｎｅｓｅＡｃａｄｅｍｙｏｆＳｃｉｅｎｃｅｓ，ＣＡＳ，Ｂｅｉｊｉｎｇ１０００４９，Ｃｈｉｎａ； ３．ＤｅｐａｒｔｍｅｎｔｏｆＧｅｏｇｒａｐｈｙ，ＵｎｉｖｅｒｓｉｔｙｏｆＷｉｓｃｏｎｓｉｎＭａｄｉｓｏｎ，ＭａｄｉｓｏｎＷＩ５３７０６，ＵＳＡ； ４．ＤｅｐａｒｔｍｅｎｔｏｆＧｅｏｇｒａｐｈｙ，ＵｎｉｖｅｒｓｉｔｙｏｆＧｕｅｌｐｈ，Ｇｕｅｌｐｈ，ＯｎｔａｒｉｏＮ１Ｇ２Ｗ１，Ｃａｎａｄａ）