一种基于横向或纵向扫描线处理的填洼算法

针对传统PlanchouDarbous算法采用递归执行方式导致处理效率低下等问题,该文提出了一种适合基于CUDA的GPU并行执行的DEM洼地填充算法。此算法采用淹没-排水思路,多次推进式扫描完成排水;在串行执行和并行执行两种方式下,算法处理有所差异,其中串行算法代码更短;串行算法的时间复杂度为O(n),测试结果表明串行执行时快于PlanchonDarboux算法,且在CUDA环境下的并行执行用时不到串行执行用时的一半。算法思路简单,易于编程实现。     

海河流域近60a降水极值的频率分析及时空分布特征

以海河流域为研究区域,利用MK和F检验对站点的年最大日降雨序列进行趋势、突变和跳跃分析。基于水文气象分区线性矩法进行一致区的划分、最优分布的选择和降水极值的频率估计值计算,并分析其空间分布特征。结果表明:MK趋势和突变检验显示只有8个站点(22.9%)呈现显著下降趋势,通过了信度为0.1的显著性水平检验,特别是京津唐区域具有显著的下降趋势,突变时段主要发生在1980-1990年。均值和方差的跳跃性显示大部分站点都呈现出显著的向下跳跃,主要分布在流域的滦河子流域和北三河水系;趋势和跳跃的综合分析能够对降水极值有可能引起的旱涝灾害进行更全面合理的判断。不同重现期下的频率估计值的空间分布总体趋势一致,从东南到西北、从沿海到内陆逐渐减少,并与地形表现出很好的一致性;降水极值空间分布的中心主要集中在滦河子流域的遵化和青龙附近;大部分站点50 a一遇的估计值和AMP序列的最大观测值能够保持较好的一致性,间接反映了估计值的合理准确性。     

利用区域气候模式RIEMS产品分析日蒸散量及其影响

利用区域气候模式RIEMS输出的各种气象参数,采用了BEF等4种不同方法计算了沂沭河上游流域的潜在蒸散量,并与该流域6个气象站实测蒸发数据计算的陆面潜在蒸散量进行了比较。结果表明,根据平均偏差、平均绝对偏差、均方根差和相关系数指标的综合判断,该4种方法的估测精度从高到低依次为双线性曲面回归经验函数法(BEF)、Hargreaves-Samani(Harg)法、Pristley-Tayler(P-T)法和Penman-Monteith(P-M)法。在时间序列上,4种方法计算的逐日蒸散量与观测值呈相同的变化趋势,但计算值在蒸散发最强、最弱和降水最多、气温最高的7-9月有较大差异。BEF法估测的精度最高,与观测值最接近,Harg法、P-M法和P-T法都有明显的偏高现象。BEF法只需要较少的参数就能得到较高的估测精度,因此可作为利用区域气候模式RIEMS产品计算沂沭河流域蒸散量的首选方法,进而为RIEMS模式中耦合的陆面水文过程模型TOPX提供满足精度要求的日蒸散量驱动参数。     

地区线性矩法估算的暴雨频率设计值空间连续性问题探讨

地区线性矩法(Regional L-moments Analysis,RLMA)是目前最先进的频率估算方法之一,自20世纪90年代起被广泛应用于美国及其他一些国家的防洪设计标准的分析计算中;我国近年也在太湖流域暴雨频率分析中成功地采用推广了这一水文气象途径的频率计算方法,积累了实际经验。但是,由于地面雨量站点有限、资料长度有限以及站点空间分布的不均匀等问题,研究区中相邻水文气象一致区之间暴雨频率估算值会出现空间不连续性的现象,这个空间不一致的问题尚没有得到合理的解决。在简要介绍地区线性矩法推求频率估计值的基础上,着重讨论并介绍一种概念清楚、简单易操作的基于空间往返二次内插的校正方法,来解决频率估算值空间不连续的问题。这种空间内插平差主要通过构造一个与研究区实测站点控制面积相应分辨率的虚拟网格站网来实现。研究结果表明,经过平差后校正的研究区内各站点的频率估计值的经验频率与理论概率更加接近,频率估计值的空间分布也更加合理。     

多普勒雷达降水反演及其与改进的区域气候模式RIEMS的同化研究

在前人构建的陆面水文过程模型TOPX和区域气候模式RIEMS的耦合模式的框架基础上,针对耦合模式中RIEMS对降水和蒸散发的模拟精度较低,RIEMS和TOPX模式之间尺度不匹配等进行研究,结果如下:(1)多普勒雷达资料的定量估测降水方面。首先,获取研究区最佳的雷达反射率-雨强(Z-I)关系。选择改进     

BP神经网络在多普勒雷达降水量的估测中的应用

利用2005年4次降雨过程的多普勒雷达体扫的回波强度资料及相应的雨量计观测资料, 通过BP神经网络方法来估测临沂地区的降雨量, 同时以改进的最佳窗概率配对法建立的Z-R关系估测的降雨量为对照, 进一步验证BP神经网络方法的优越性。根据各个站点的平均相对误差、 均方根差、 相关系数和相关曲线斜率4个指标的比较, 小时雨量和累计降雨量估测结果表明： BP神经网络估测精度要明显优于Z-R关系式, 训练样本的精度高于检验样本的精度, BP神经网络估测的降雨量与站点实测雨量吻合性较好, 能够较真实地反映地面降雨情况; Z-R关系式估测的降雨量随着雨强的不同表现为不同程度的低估现象。     

淮河流域极值降雨频率估计时空分布一致性检验和纠正

由于地面雨量站点空间不均及站点资料长度有限、相邻时段和相邻一致区的最优分布函数选择不同,造成降雨频率估计值出现交叉和空间不连续的现象,时空不一致的问题尚没有得到合理的解决。为了获取更加完善合理、时空一致的频率估计值,利用水文气象分区线性矩法计算出的频率估计值,通过误差分摊和二次空间插值及平滑技术对时空分布一致性进行检验和纠正,消除分布曲线的交叉和边界梯度的问题。研究结果表明:通过误差分摊纠正后,消除了交叉异常曲线,各个时段的频率估计值相互制约更加合理;通过空间二次插值纠正后,边界梯度和误差得到明显改善,各站点的频率估计值的经验频率与理论概率更接近,空间分布也更合理。多时段不同重现期下的频率估计值空间分布总体趋势一致,最大值主要分布在淮河流域的东北沂蒙山区,最小值分布在流域的北部平原地区,与站点的实际观测相一致,可为工程设计暴雨和洪水提供基础资料和重要水文依据。     

基于贝叶斯三角帽法的多源降水数据融合分析及应用

目前的降水产品依然存在较大的不确定性,采用多源降水数据融合可以更准确地估计降水量和空间分布情况。为实现无资料地区的数据融合,本文在不使用任何先验信息的前提下,通过整合站点插值、卫星遥感和再分析的降水产品,基于贝叶斯三角帽(Bayesian-Three Cornered Hat, BTCH)法,融合多源降水数据,探究不同输入数量的降水产品对于融合数据精度的影响以及每个降水产品对于融合数据精度的贡献率,并在黄河源区进行应用。结果表明：在月尺度上,融合数据性能优于原始降水产品;在日尺度上,融合数据性能明显高于卫星遥感和再分析降水产品,但低于基于站点的降水产品CHM＿PRE;2套基于站点的降水产品CN05.1和CHM＿PRE对于融合数据有最大的贡献率。在黄河源区的应用表明,该数据融合方法确实能够更准确地估计降水量,可应用于无实测降水资料地区,为数据融合分析及应用提供参考。     

统计权重矩阵法在雷达估测降水集成中的应用

Z～R关系法和6种雷达雨量计联合法反演的区域降水量与雨量计观测得到的降水场存在较大的误差,将这7种降水估测结果作为信息源,采用统计权重矩阵法对上述7种反演结果进行集成分析,提出了一种改进雷达估测降水的方法。结果表明:在被集成资料中,Z～R关系法估测的降雨场具有明显的偏低现象,精度最差,6种雷达雨量计联合法的估测精度明显优于Z～R关系法。通过统计权重矩阵集成后,精度比集成前所有方法均有明显提高,尤其是降水场分布形势和降水中心的强度都与雨量计场非常吻合。集成得到的降水空间分布场能够较真实地反映地面的降水情况,可以在估测区域降水量中进行业务试用。     

地区线性矩法在太湖流域暴雨频率分析中的应用

介绍了基于水文气象途径的地区线性矩法的概念,通过基于次序统计量的线性矩进行参数估计与基于水文气象一致区的地区分析法相结合,以太湖流域1d时段的年极值降雨资料为例,进行暴雨频率分析。应用水文气象一致区的判别准则,将太湖流域划分为8个水文气象一致区;综合考虑三种拟合优度检测方法,选择1~8区的最优分布线型分别为:GEV、GLO、GEV、GEV、GNO、GNO、GEV、GNO;根据地区分析法原理,估算各雨量站的暴雨频率设计值。分析表明:太湖流域各重现期下的年极值降雨空间分布形态基本一致,西南山区是太湖流域的暴雨高风险区,应该在地区防洪规划中引起重视。结果表明:地区线性矩法具有很高的学术和实用价值,建议在全国范围内推广,作为防洪规划的顶层设计和基础工作,以满足工程防洪设计、地区防洪规划、山洪预警和城市防涝防洪规划等方面的需求。