1、北京市城市雨水排除规划设计标准研究北京市降雨变化趋势研究,北京市气象局气候中心 北京市城市规划设计研究院二零零九年十二月,一、概述二、气候变化事实分析三、北京市历史降雨变化趋势四、气候变化对北京降雨的影响五、北京地区未来降雨变化趋势预估六、北京市未来不同历时最大降水量重现期分析七、结论与展望,主 要 内 容,1.1 研究背景,城市雨水排除系统是现代城市的重要基础设施，对城市的经济发展、城市品味、人们的生产和日常生活等有着极其重要的影响。 北京市城市雨水排除系统规划设计标准研究。其中，北京市气象局气候中心主要承担本项目的子课题2北京市降雨变化趋势研究。,1.2 研究目标与内容,按照北京市防洪排水
2、安全的要求，研究适合北京市实际情况的雨水排除系统规划设计标准。本项目的总目标：提高城市雨水基础设施的运行能力，减少洪涝灾害。 北京市降雨变化趋势的研究目标是，基于北京现有最新的雨量观测记录数据，研究北京地区降雨时空变化规律，预估未来降雨趋势及不同历时重现期的降雨量，为北京城市雨水排除规划设计标准的修订提供参考。,1.2 研究目标与内容,按照任务书的要求，主要完成两部分的工作：1、资料的整理和提供根据北京市暴雨分区研究结论，整理和提供北京地区代表气象站逐分钟精细化降雨记录资料。每站每年选取1820场最大降雨样本。2、未来降水趋势预估 2.1气候变化的趋势：分析北京地区近五十年来降水的空间和时间变
3、化趋势。时间变化趋势分析的变量包括：年总降水量、短历时降雨强度。暴雨强度的空间分布。极端降雨事件的变化趋势。分析北京地区代表站（观象台）近百年降水量变化趋势。 2.2未来降水趋势预估：分析未来（2020、2050年）北京地区的降雨趋势（年总降水量、最大日降雨量、最大过程降雨总量）。分析不同历时（5、10、15、20、30、45、60、90、120、150、180、240、360、720、1440分钟）重现期（1、2、3、5、10、20、50、100年）降雨量。,1.4 研究方案,（1）北京地区20个气象站资料 20个气象站：海淀、朝阳、观象台、丰台、顺义、汤河口、密云、怀柔、上甸子、平谷、通州
4、、大兴、斋堂、门头沟、房山、霞云岭、延庆、昌平、佛爷顶、石景山。其中，平原地区6站包括：观象台、丰台、顺义、密云、房山、昌平。（2）北京地区121个雨量站资料 北京市潮白、蓟运、北运、永定、大清河五个水系共121个雨量站1956-2008年降雨资料，包括年、汛期、10分钟、30分钟、60分钟、1440分钟等六个不同时间尺度。,1.资料的收集整理,1.4 研究方案,（3）短历时代表气象站逐分钟雨量资料短历时（24小时）资料取代表站：观象台、密云、海淀，自建站起2008年12月31日的自记降水纸记录资料。应用国家气象局预测减灾司提供的“降水自记纸数字化处理软件”（2004）对降水记录纸进行扫描、扫
5、描检查、降水曲线提取、降水强度数据转换和质量检查，提取逐分钟雨量资料。除了降水自记纸预处理系人工作业外，降水自记纸数字化处理软件由四部分组成：扫描、扫描检查、降水曲线提取、降水强度数据转换和质量检查。,1.资料的收集整理,1.4 研究方案,1.资料的收集整理,（3）短历时代表气象站逐分钟雨量资料,降水自记纸数字化处理业务流程如下：1.降水自记纸预处理：对降水自记纸进行扫描输入和数字化处理之前，首先要对降水自记纸进行一系列技术处理，为降水自记纸图像扫描和图形数字化的顺利和准确进行做好准备。这种技术处理，称为降水自记纸预处理。2.降水自记纸图形扫描输入，形成扫描图像文件。3.扫描图形检查（扫描质量
6、、相关信息和补遗）。4.降水曲线自动提取，人机对话认可或修正，形成降水曲线文件。5.计算降水分钟强度，按规定格式存储。按站形成降水自记数据（降水分钟强度）文件和标准（降水分钟强度、降水小时强度）文件。,1.4 研究方案,1.资料的收集整理,（3）短历时代表气象站逐分钟雨量资料,1.4 研究方案,分析北京市近百年和近五十年来的降雨变化趋势，包括多个历时（年、1440分钟、720分钟、360分钟、120分钟、60分钟和5分钟）的时间变化趋势，并应用arcgis绘制空间分布图，探讨空间变化规律。重点分析1981年至今的降雨趋势，北京市极端降雨事件频率和强度的变化趋势。,2.北京市降雨变化趋势研究,1
7、.4 研究方案,应用国家气象局发布的中国地区气候变化预估数据集和天气发生器方法预估北京地区未来（2020、2050年）的降雨趋势（日降雨量极值、次降雨量极值）。,4.北京市未来降雨趋势分析,应用极大似然法和概率加权法（PWM），对中国东部地区日降水广义极值分布（GEV）和广义帕雷托分布（GPD）模拟的情况进行比较。利用概率散点图和百分位散点图对北京短历时降水的情况进行拟合效果检验，最终GPD对降水的拟合效果最好，并利用GPD方法预估了北京城区代表站不同历时重现期降雨量。,5.北京市未来不同历时最大降雨量重现期分析,1.5 主要研究成果,收集整理并提供了北京地区代表气象站逐分钟降雨信息化资料。观
8、象台站为19412008年，海淀为19752008年，密云为19612008年。分析了北京市短历时降雨变化的趋势，包括时间变化趋势和空间分布规律。分析了极端降雨的时间变化趋势，并探讨了影响北京地区强降雨的因子。分析了气候变化对北京地区降雨的影响，包括水汽输送总量和水汽输送路径的变化规律。预估了北京地区未来（2020、2050年）降雨趋势，并预测了各个历时（5、10、15、20、30、45、60、90、120、150、180、240、360、720、1440分钟）的重现期（1、2、3、5、10、20、50、100年）降雨量。,二、气候变化事实分析,2.1 中国及华北地区气候变化事实2.2 北京市
9、气候变化事实,2.1 中国及华北地区气候变化事实,降水,华北区域暴雨频数与暴雨雨量均表现出显著下降趋势。1956-2007 年间，华北区域降雨频数下降是导致降雨量显著下降的首要原因，而其最主要原因是暴雨频数下降导致暴雨雨量下降。,气候变化国家评估报告华北地区年降水量极值的面积覆盖率表现出明显的下降趋势，与此同时，年降水日数极端偏多的范围表现出更为明显的下降趋势，这就造成了平均降水强度极端偏强范围的明显增加。,气候变化对北京市的影响体现在诸多方面，诸如农业、水资源、森林和其它自然生态系统、城市运行保障、人体健康等。就对水资源的影响而言，气候变暖加剧了北京市水资源的供需矛盾。随着气温的升高，降雨量
10、也呈现出逐年降低的趋势，对北京市水资源供需平衡造成了不利影响。而随着社会经济活动和居民生活水平的提高，加剧了水资源的供需矛盾。另外，极端气候事件频发，局地强降雨频现。,2.2 北京市气候变化事实,三、北京市历史降雨变化趋势分析,3.1 一百多年来北京降雨量的趋势3.2 极端降雨事件的趋势分析3.3 短历时强降雨的演变趋势3.4 影响北京地区强降雨的因子分析,3.1 一百多年来北京降雨量的趋势,（1）一百多年来，北京地区有两个多雨时段，一段是1880年代中期至1890年代中期，另外一段是1950年代；三个少雨期是1850年代中期至1860年代中期、1940年代、1970年代中期至今。（2）195
11、0年代以来降雨量呈减少趋势，1960年以后年降雨距平百分率基本为负值，即降雨量少于常年平均值，平均每10年降雨量减少16毫米。（3）19992007年连续9年降雨偏少，造成水资源总量、地表水和地下水资源量持续偏少。北京每年的降雨很容易产生重大的偏离，以致出现旱涝无常的状况。据统计，北京地区旱多于涝，但大涝远多于大旱。,时间变化趋势,近五十年北京平原地区年降水量呈减少趋势，减少速率为31毫米/10年。尤其是近9年（19992007年）连续为偏少年，年平均降水量为462.1毫米，比近30年平均（556.0毫米）偏少17%。目前北京处于20世纪50年代后的相对少雨时段，干旱现象十分严重。近五十年降水
12、最少的三个年份分别为1999年、1965年、2006年。降水日数也呈减少趋势。,降水日数,降水距平百分率,3.1 一百多年来北京降雨量的趋势,时间变化趋势,3.1 一百多年来北京降雨量的趋势,位于怀柔、密云附近的大于600mm降雨量等值线的范围随年代发生明显改变。1980年代降雨中心主要位于怀柔、密云一带。1990年代降雨中心逐渐向城区移动，城区北部的降雨量逐渐增加。另外，位于房山的霞云岭附近出现一个降雨高值区，北京的主要降雨落区呈东北西南向分布。2000年之后降雨中心量值减少非常明显，中心位置回撤到密云附近。,空间变化趋势,1980年代,1990年代,2000年代,3.1 一百多年来北京降雨
13、量的趋势,北京地区20个气象观测站近三十年降雨量变化趋势： 全市各站降雨均呈减少趋势，东部地区减少最为明显，每十年减少4050毫米。,空间变化趋势,3.2 极端降雨时间的趋势分析,19712007年，北京年极端降雨发生的频次在减小，1970年代极端降水事件的频次为6.4次/年，到近十年（19992007年）仅为2.6次/年。但日降雨强度有增大的趋势。这表明，北京的暴雨集中度在增加，极易出现强降雨。,极端降水事件近年我国各地频繁出现多个破历史记录的极端气候事件。,3.3 短历时强降雨的演变趋势,5分钟,北京地区代表站（观象台、海淀和密云）近50年短历时5分钟年最大雨强，19412008年5分钟雨
14、强年际间差异较大，最大值主要集中在1950年代和1960年代。其中，1959年、1963年、1976和1951年最高，5分钟雨量在19.221.5毫米之间。从1980年代中期1990年代末，5分钟雨强呈减小的趋势。但从2001年至今，北京地区5分钟雨强呈增加的趋势。,3.3 短历时强降雨的演变趋势,60分钟,19412008年北京地区60分钟雨强的年际差异明显，从18.2129.7毫米/60分钟。1950年代和1960年代60分钟雨强最大，从1985年至2002年，60分钟雨强呈减小的趋势。2003年至今，海淀和密云的60分钟雨强呈增加的趋势，观象台变化趋势不大,3.3 短历时强降雨的演变趋势
15、,60分钟,利用北京地区121个水文站60分钟最大降雨量资料，从1980年代中期到2001年，北京地区60分钟时间间隔内发生的强降雨极值呈明显减小的趋势，而从20022008年60分钟雨强极值呈增加的趋势，尤其是近3年来，60分钟强降雨为近五十年来最大值。各年代的60分钟雨强极值从大到小依次为：1980年代1970年代2000年代1990年代。1980年代最大，1990年代1小时雨强最小，20012008年雨强比1990年代增大。,3.3 短历时强降雨的演变趋势,120分钟,360分钟,3.3 短历时强降雨的演变趋势,3.3 短历时强降雨的演变趋势,1980年代后，一方面北京地区的降雨明显偏少
16、，另一方面，发生区域性大雨的频次明显减小，近年来北京暴雨多呈现出局地性的强对流天气，它是由尺度很小的积雨云单体或单个雨团所引起的，这种天气的产生具有一定随机性，并且生命史极短，所以生消发展都很迅速，这种特点大大增加了预报难度，一旦出现，极易造成局地积涝，使得交通瘫痪、市民出行受到影响。另外，由于现代城市的发展，不透水路面增多，地表容易形成径流，即使遇到一般暴雨也会出现局部地区积水，影响城市交通。,1960和1970年代发生区域性大雨的频次较高，1970年代最为明显，一次大雨过程可以导致区域内大范围台站24小时之内同时出现大于25毫米的降雨。,2.广义帕雷托分布: GPD的最大优点在于，它直接由
17、原始资料数据以给定门限值为标准来抽取超过该门限值的极大(或极小)值，其所需资料年限大大节省，从而增加了极值的样本量，该分布更加符合实际。广义Pareto分布函数为：分布密度函数(PDF)为：GPD的重现期值公式为：,6.1 方法介绍,2.广义帕雷托分布:广义帕雷托GPD分布参数估计方法，采用概率加权矩法(PWM)，不但有较高的精度可与公认的“极大似然法”相媲美，而且其计算简便易行，不需要利用繁琐的迭代计算。PMW定义为其通式可写为根据前三个参数无偏估计式，GPD分布参数的估计通式为,6.1 方法介绍,2.广义帕雷托分布:（4）广义帕雷托分布拟合及其效果检验选取我国不同区域的几个代表站，用par
18、eto 函数拟合1050 年不同长度样本资料的极端降雨。该方法能达到较高的拟合度，且基本上资料年份越长，效果越好。这也说明极端降雨采用pareto分布拟合是可行的。,6.1 方法介绍,检验方法：概率散点图（P-P图）和百分位散点图（Q-Q图） 概率散点图： 概率散点图描述了样本的次序统计量与经验分布函数的拟合程度。对于极值模型最关心的就是数据取大值时模型是否合适，即关心较低次序统计量 (相对于较小的k)的拟合度。 百分位散点图：极值理论认为，当X的分布函数为 F(X) 时，P-P图和Q-Q图应近似为直线。这是检验模型合理性的主要根据。,6.2 广义极值分布和广义帕雷托分布对北京降雨模拟,5分钟
19、,6.2 广义极值分布和广义帕雷托分布对北京降雨模拟,（2）广义帕雷托分布对北京观象台短历时降雨拟合检验,120分钟,6.2 广义极值分布和广义帕雷托分布对北京降雨模拟,（2）广义帕雷托分布对北京观象台短历时降雨拟合检验,K-S检验是一种拟合优度检验，研究样本观察值的分布和设定的理论分布间是否吻合，通过对两个分布差异的分析，确定是否有理由认为样本的观察结果来自所设定的理论分布总体。观象台不同历时极端降雨分布拟合效果检验,6.2 广义极值分布和广义帕雷托分布对北京降雨模拟,（4）科斯（K-S）检验方法对GEV和GPD模拟效果的检验,6.3 广义帕雷托分布的重现期结果,计算步骤：（1）门限值的选定
20、利用年平均交叉率接近于1的方法计算出不同历时的门限值。,6.4 北京市各历时降雨量的重现期结果,七、结论与展望,北京市历史降雨变化趋势研究（1）北京地区18412008年年降雨量变化曲线显示：北京有两个多雨时段和三个少雨时段。近三十年，全市各站降雨均呈下降趋势，东部地区下降最为明显。同时，降雨日数也呈减少趋势。（2）19712008年，北京地区年极端降雨发生的频次在减小，但日降雨强度有增大的趋势，北京的暴雨集中度在增加。（3）北京地区近五十年短历时强降雨呈波动变化的规律，最大值主要集中在1970年代之前，以1950和1960年代强降雨量最高，19412008年5分钟雨强变化趋势不大，而60分钟
21、、120分钟和360分钟年最大雨强从1980年代后期至1990年代末呈明显减小的趋势。2001年至今，短历时降雨雨强极值呈增加的趋势，近十年短历时雨强极值处于由偏小转向偏大的过渡时期。（4）19701980年代，北京地区强降雨主要为全市区域性强降雨，1990年代之后北京的短历时强降雨呈现出局地性的特征，降雨分布不均，强降雨中心大致成东北西南向带状分布。近三十年，全市短历时雨强1980年代最高，1990年代最小，20012008年雨强比1990年代增大。尤其是近3年来，1小时强降雨为近五十年来最高值。（5）1960和1970年代北京市发生区域性大雨的频次较高，1970年代最为明显，1980年代之
22、后发生区域性大雨的频次明显减小。,七、结论与展望,北京市未来不同历时最大降雨量重现期分析（1）选取中国东部地区78个站50年(1951-2000)逐日降雨资料作为研究对象，用广义帕雷托分布函数拟合1050 年不同长度样本资料的极端降雨，经检验，该方法能达到较高的拟合度，且基本上资料年份越长，效果越好。说明极端降雨采用广义帕雷托分布拟合是可行的。（2） 对北京城区代表站观象台5分钟和120分钟降雨极值分布进行GEV和GPD拟合，检验结果表明两种分布函数均能拟合出观象台降雨变化。（3）应用广义极值分布和广义帕雷托分布方法对极端降雨的拟合都非常好，两种方法对不同重现期下不同历时降雨量的估算结果可信。通过K-S检验，广义帕雷托分布GPD对极值降雨的模拟效果要好于广义极值GEV分布，因此采用GPD方法预估不同历时的重现期变化。（4）通过对实测资料的对比拟合，重现期模拟结果可靠。,谢谢！,
