标题 | 流域下垫面变化影响下入库设计洪水的修订 |
范文 | 丛娜 冯平 摘要:受基础设施建设、工农业生产、生态环境建设等人类活动影响,海河流域下垫面发生了一些变化,导致设计洪水也随之发生相应的变化。现根据王快水库实测降雨径流资料,对降雨和洪峰进行趋势与变异分析,识别出入库洪水序列的变异点,然后采用降雨径流相关法、峰量相关法对洪量和洪峰序列进行了修订,并对修订前后的洪水序列进行了频率分析。结果表明,下垫面变化的影响会使入库设计洪水减小,相应水库校核防洪标准10 000年一遇的洪峰流量减少了2.38%。探讨下垫面变化后入库设计洪水的修订问题,对于调整流域防洪对策和防洪工程的布局具有重要的现实意义。 关键词:下垫面变化;设计洪水;洪水修订;王快水库 中图分类号:TV124文献标识码:A文章编号: 16721683(2014)02000605 Revision on Reservoir Inflow Design Flood under the Variation Impact of Underlying Surface CONG Na,FENG Ping (School of Civil Engineering,Tianjin University,Tianjin 300072,China) Abstract:Under the influence of infrastructure construction,industrial and agricultural production,ecological environment construction,and other human activities,the underlying surface of the Haihe River Basin has been changed,which leads to corresponding variation in the design flood.In this paper,based on the data of measured rainfall and runoff in the Wangkuai Reservoir,the trend and aberrance point of rainfall and flood peak were analyzed and used to identify the aberrance point in the series of reservoir inflow flood.The rainfallrunoff relationship method and the peakvolume correlation method were used to modify the series of flood volume and flood peak in the Wangkuai Reservoir,and flood frequency analysis was performed on the flood series before and after modification.The results showed that the variation of underlying surface can cause the reduction in the reservoir inflow design flood,and the corresponding flood control standard of the worst flood once in ten thousand years has decreased by 2.38%.The study of the revision on reservoir inflow design flood under then variation of underlying surface is of important significance for the adjustment of flood control measures and layout of flood control engineering. Key words:variation of underlying surface;design flood;flood revision;Wangkuai Reservoir 由于人类活动对下垫面的扰动,流域水文序列的一致性发生了一定的变化,这将影响到防洪工程布局与防洪对策的调整。因此,估算下垫面变化条件下的设计洪水(即对原洪峰序列和洪量序列进行修订,具有重要的理论和现实意义)。国内外许多学者针对下垫面变化对洪水、径流(包括设计洪水)的影响问题进行了大量研究[14]。 当水文站以上流域的人类活动频繁且有水量耗用时,必须对观测值进行修订,这一过程即一致性修订[5]。目前国内外水文序列修订方法主要有逐项调查法、蒸发插值法、模型法等[6],但这些方法往往存在计算精度不高、还原项目难以考虑完全、天然径流还原基础较难一致等问题。 根据流域内的降雨、径流资料,可以建立降雨―径流相关关系[7],对流域径流序列进行修订。测站以上的现状洪水径流量与天然洪水径流量的差值,就是流域洪水径流量在下垫面变化后引起的减少量,利用这一减少量便可以对变化前洪量序列进行修订。因此。本文采用降雨径流相关法,将洪水径流量“还现”到现状下垫面条件下的径流量,即“向后还原”,来实现对洪量序列的修订;然后再采用洪水峰量相关法,进行洪峰序列的修订,从而达到对入库设计洪水的修订目的,为完善水库的防洪调度与管理提供支持。 1研究方法 1.1水文序列的趋势和变异分析 对水文序列进行一致性修订前,首先需要对水文序列的变化趋势进行分析,即趋势性和变异性,并通过一定的方法检验。检验趋势性与变异性的方法众多,其中非参数的方法运用较多。熊立华等[8]用Spearman 和MannKendlla秩次相关检验法检验了长江宜昌站年径流趋势,王荣恩等[9]用方差分析法分析了黑龙江省主要江河水文序列,张爱静[10]用双累计曲线法分析了东北地区流域的径流序列,谢平等[11]采用成因分析法与统计分析法分别对水文序列确定性成分和随机性成分进行识别与检验。本文将采用Spearman[12]、MannKendall秩次相关检验法[13]、滑动平均法[14]、滑动F检验法[15],双累积曲线法[16],对海河流域大清河水系主要水库站暴雨量和洪峰流量序列进行趋势性和变异性分析。 1.1.1趋势检验方法 (1)Spearman秩次相关检验法。Spearman秩次相关检验法是通过分析时间序列与其时序的相关性来检验趋势性的,它是采用t检验法检验序列秩次相关系数是否异于零,来判断序列的趋势变化。其假定统计量T服从自由度为(n2)的t分布,并用统计量T作为时间序列趋势性大小的衡量标度。T越大,则序列的趋势性变化越显著。 T=r(n-4 1-r2)1/2,r=1-6∑n i=1(Ri-i)2 n3-n(1) 式中:r为秩次相关系数;n为序列长度;i为序列时序;Ri为秩次。 (2)MannKendall秩次相关检验法(MK方法)。MK方法趋势检验是在序列平稳的条件下,首先假定时间序列无趋势,计算其统计量S和MK检验统计量Z。当Z为正值(负值)时,表明时间序列呈上升(下降)趋势。 S=∑n-1 i=1∑n j=i+1sgn(Xj-Xi)(2) Z=(S-1)/var(S)S>0 0S=0 (S+1)/var(S)S<0(3) 当n≥10时,统计量S近似服从正态分布,其均值E(S)和方差var(S)分别为: E(S)=0,var(S)=n(n-1)(2n+5)/18(4) 式中:n为序列长度。 1.1.2变异检验方法 (1)双累积曲线法。双累积曲线法是在直角坐系中绘制同时期内一个参考变量的连续累积值与另一个被检验变量的连续累积值的关系线,是检验变量变化的常用方法。采用降雨―径流双累积曲线法,能便捷直观地确定径流序列的突变时间。 (2)滑动F检验法。假设变异点m前后,两序列总体的分布函数分别为F1(x)和F2(x),从总体F1(x)和F2(x)中分别抽取容量为n1和n2的两个样本,要求检验原假设:F1(x)=F2(x)。 设x1,x2,…,xn1与y1,y2,…,yn2分别代表变异点m前后两个样本系列,则样本均值和方差分别为: X=1 n∑n1 i=1Xi,y=1 n∑n1 i=1yi(5) S21=1 n1∑n1 i=1(Xi-X)2,S22=1 n2∑n2 i=1(yi-y)2(6) 令M=n1S12/(n1-1),N=n2S22/(n2-1),若M>N则F=M/N,自由度为V1=n1-1,V2=n2-1;若M 本文采用降雨―径流相关法“还现”径流量。首先根据实测的场次降雨量和洪水流量资料,分别计算与降雨、径流相关的水文要素,包括前期影响雨量Pa、平均降雨量和地面径流深;然后以洪水序列突变点为分界,分别点绘两个序列的P+Pa~R相关图,据此便可分析流域下垫面变化对次洪量的影响。即根据下垫面变化前某次洪水所对应的降雨量,从相关图中分别查出下垫面变化前后的径流深,由其变化值就是下垫面一致性修订值,以此进行逐次洪水修订,即可实现洪量系列的修订。 通过上述降雨―径流相关法还可以对入库洪量序列进行修订。因此只要再建立不同时段洪量与洪峰的相关图,就可以根据不同时段洪量的变化幅度,并估算确定洪峰的变化量,从而得到修订后的洪峰流量序列,达到对入库洪峰流量进行“还现”的目的。 2实例分析 2.1流域基本情况 大清河水系位于海河流域中部,西起太行山,东临渤海湾,北邻永定河,南界子牙河,东西长约275 km,南北宽约200 km,流域面积34 060 km2。流域范围内降雨量春冬季较少,夏季较多,主要集中在7月、8月。流域内主要有王快、西大洋、龙门、横山岭等水库。本文以王快水库为例,来探讨下垫面变化影响下设计洪水修订问题。 王快水库位于河北省曲阳县郑家庄西、大清河水系沙河上游,控制流域面积为3 770 km2,总库容13.89亿m3,是一座以防洪为主,结合灌溉、发电等综合利用的大(Ⅰ)型水利枢纽工程。水库始建于1958年6月,2002年开始进行以扩建溢洪道和大坝加固为主的水库除险加固工程,工程完成后,水库校核防洪标准达10 000年一遇。20世纪80年代以来,水保工程、退耕还林等活动使其流域下垫面发生了明显变化。 2.2降雨和流量序列的趋势与变异分析 从1956年―2005年王快水库实测降雨径流资料中,选取年最大洪峰流量和相应的次降雨量,分别采用MK检验法、滑动平均法和Spearman法对该降雨和洪峰流量序列进行趋势分析,结果见表1。 从表1可以看出,王快水库的次降雨量和洪峰流量均有下降趋势,其中洪峰序列的检验结果为下降显著,而对降雨序列的MK和Spearman检验结果分别为显著和不显著。由图1可以看出,位于序列首端的1963年出现了特大暴雨,这可能导致了降雨量序列的趋势统计结果为显著或不显著的下降趋势。 图1和图2为降雨和洪峰流量的5年滑动平均结果,从中可看出,均呈下降趋势。由图3的双累积曲线可知,洪峰流量序列变异点为1979年。并且滑动F检验法计算得到的洪峰流量序列变异点也为1979年。综上,王快水库次降雨量、入库洪峰均呈现下降趋势,变异点可确定为1979年。 图1降雨量序列滑动平均图 Fig.15year running mean of the precipitation series 图2洪峰流量序列滑动平均图 Fig.25year running mean of the flood peak series 图3降雨径流双累积曲线 Fig.3Double mass curve of the rainfall and runoff 2.3洪量的修订 选取了王快水库以上流域1956年-2004年实测降雨径流资料中42场暴雨洪水过程,分别计算每场洪水相应的降 图4降雨-径流相关图 Fig.4Relationship between (P+Pa) and R 雨总量P、洪水总量R和每场洪水的前期影响雨量Pa,建立1980年前后的(P+Pa)~R相关图(图4)。结果表明,在王快水库以上流域,相同降雨量情况下,1980年后的次洪流量与1980年前相比有所减小,说明在降雨量相同的情况下,流域下垫面条件的变化使产流量变小。由其变化值可得到降雨量P的相应次洪量的下垫面一致性修订值,因此可实现洪量系列的修订。 2.4洪峰的修订 利用上述42场暴雨洪水资料,分别建立王快水库不同时段入库洪量与洪峰的相关图 (图5),即可根据不同时段洪量的修订值,来计算确定洪峰的变化量,从而实现对洪峰序列的修订。 图5峰量相关图 Fig.5Correlation diagram of flood peak and flood volume 2.5设计洪水计算 通过上述分析,可以得到考虑下垫面变化影响后“还现”的王快水库入库洪峰和不同时段的洪量序列,这样便可在考虑历史特大洪水情况下,对修订前和修订后洪水序列进行频率分析计算,结果见表2。图6和图7分别是是修订前后洪峰序列的频率曲线。 图6王快水库洪峰序列频率曲线 Fig.6Frequency curve of the original flood peak series of the Wangkuai Reservoir 图7修正后的王快水库洪峰序列频率曲线 Fig.7Frequency curve of the revised flood peak series of the Wangkuai Reservoir 通过对修订后与修订前的洪水序列频率分析结果比较(表2),发现其均值和不同重现期的设计值均有减小,其中洪峰流量均值减小了10.82%,1日、3日和6日洪量均值分别减小了12.50%、12.73%和11.84%。相应水库校核防洪标准10 000年一遇的洪峰流量减少了238%。说明由于下垫面变化的影响,特别是水保工程、退耕还林等形式的下垫面的变化,王快水库的洪峰流量和洪量均有有减小的趋势,并且小洪水的洪峰和洪量比对大洪水减少的更多。这种结果符合流域对洪水的调蓄作用,下垫面的变化对小洪水比对大洪水的调蓄作用更加明显。 3结语 本文以大清河流域王快水库为研究对象,对1956年―2005年的天然年径流系列进行一致性分析,确定了水文序列的变异点为1979年。然后用降雨径流相关法,修订了洪量序列,根据峰量相关进一步修订洪峰值,实现了下垫面影响较大水文序列的一致性修订。通过对王快修订前后水库入库洪水序列的频率分析计算,水保工程、退耕还林等形式的下垫面的变化导致流域的产流量减少,设计洪水均呈变小的趋势。相应水库校核防洪标准10 000年一遇的洪峰流量减少了2.38%。 基于传统的相关分析法进行下垫面变化影响下设计洪水的修订,其本质是对除降水以外影响径流的各种因素进行的统一修订,即把除降水以外的其他因素都归于下垫面变化,其计算精度还需进一步讨论。 参考文献(References): [1]陈民,谢悦波,冯宇鹏.人类活动对海河流域径流系列一致性影响的分析[J].水文,2007,(3):5759.(CHEN Min,XIE Yuebo,FENG Yupeng.Influence of Human Activity on Consistency of Runoff Series in Haihe River Basin[J].Journal of China Hydrology,2007,27(3):5759.(in Chinese)) [2]姚治君,管彦平,高迎春.潮白河径流分布规律及人类活动对径流的影响分析[J].地理科学进展,2003,22(6):599606.(YAO Zhijun,GUAN Yanping,GAO Yingchun.Analysis of Distribution Regulation of Annual Runoff and Affection to Annual Runoff by Human Activity in the Chaobaihe River[J].Progress in Geography,2003,22(6):599606.(in Chinese)) [3]THERESIA P,BRUNO M.Trends in Flood Magnitude,Frequency and Seasonality in Germany in the Period 19512002[J].Journal of Hydrology,2009,371:129141. [4]张建中,毛慧慧.下垫面变化条件下设计洪水的修订方法研究[J].水利水电技术,2012,43(12):14,8.(ZHANG Jianzhong,MAO Huihui.Study on Method for Modification of Design Flood under Condition of Change of Underlying Surface[J].Water Resources and Hydropower Engineering,2012,43(12):14,8.(in Chinese)) [5]周鹏,张中华.天然径流一致性修正方法应用及问题探讨[J].东北水利水电,2008,26(12):1112,71.(ZHOU Peng,ZHANG Zhonghua.Discussion and Application of Natural Runoff Consistency Correction Method[J].Water Resources & Hydropower of Northeast China,2008,26(12):1112,71.(in Chinese)) [6]沈宏.天然径流还原计算方法初步探讨[J].水利规划与设计,2003,(3):1518,47.(SHEN Hong.Probing on the Reducing Calculation Method of Natural Runoff[J].Water Resources Planning and Design,2003,(3):1518,47.(in Chinese)) [7]张强,宋智明.于桥水库降雨径流相关图预报方案修订分析[J].海河水利,2005,(5):5051.(ZHANG Qiang,SONG Zhiming.Analysis of Yuqiao Reservoir RainfallRunoff Correction Graph Forecast Program Revision[J].Haihe Water Resources,2005,(5):5051.(in Chinese)) [8]XIONG Lihua,GUO Shenglian.Trend Test and ChangePoint Detection for the Annual Discharge Series of the Yangtze River at the Yichang Hydrological Station[J].Hydrological Sciences,2004,49(1):99112. [9]王恩荣,耿鸿江.黑龙江省主要江河水文要素的周期分析[J].水文,1995,(1):4247,53.(WANG Rongen,GENG Hongjiang.The Main Elements of the Heilongjiang River Hydrological Cycle Analysis[J].Journal of China Hydrology,1995,(1):4247,53.(in Chinese)) [10]张爱静.东北地区流域径流对气候变化与人类活动的响应特征研究[D].大连:大连理工大学,2013:3345.(ZHANG Aijing.Research on the Streamflow Responses to Climate Variability and Human Activities in Northest China[D].Dalian:Dalian University of Technology,2013:3345.(in Chinese)) [11]谢平,陈广才,夏军.变化环境下非一致性年径流序列的水文频率计算原理[J].武汉大学学报(工学版),2005,38(6):69,15.(XIE Ping,CHENG Guangcai,XIA Jun.Hydrological Frequency Calculation Principle of Inconsistent Annual Runoff Series under Changing Environments [J].Engineering Journal of Wuhan University,2005,38(6):69,15.(in Chinese)) [12]许全喜,石国钰,陈泽方.长江上游近期水沙变化特点及其趋势分析[J].水科学进展,2004,(4):420426.(XU Quanxi,SHI Guoyu,CHEN Zefang.Analysis of Recent Changing Characteristics and Tendency Runoff and Sediment Transport in the Upper Reach of Yangtze River[J].Advances in Water Science,2004,(4):420426.(in Chinese)) [13]刘叶玲,翟晓丽,郑爱勤.关中盆地降水量变化趋势的MannKendall分析[J].人民黄河,2012,34(2):2830,33.(LIU Yeling,ZHAI Xiaoli,ZHANG Aiqin.Analysis of Precipitation Trend in the Guanzhong Basin Based on the Mann Kendall Method[J].Yellow River,2012,34(2):2830,33.(in Chinese)) [14]郭鹏,陈晓玲,刘影.鄱阳湖湖口、外洲、梅港三站水沙变化及趋势分析(19552001年)[J].湖泊科学,2006,18(5):458463.(GUO Peng,CHEN Xiaoling,LIU Ying.Analysis on the Runoff and Sediment Transportation in the Houkou,Waizhou and Meigang Stations of Lake Poyang during 19552001[J].Journal of Lake Sciences,2006,18(5):458463.(in Chinese)) [15]陈广才,谢平.水文变异的滑动F识别与检验方法[J].水文,2006,26(2):5760.(CEN Guangcai,XIE Ping.Slide F test of ChangePoint Analysis[J].Journal of China Hydrology,2006,26(2):5760.(in Chinese)) [16]穆兴民,张秀勤,高鹏,王飞.双累积曲线方法理论及在水文气象领域应用中应注意的问题[J].水文,2010,30(4):4751.(MU Xingmin,ZHANG Xiuqin,GAO Peng,WANG Fei.Theory of Double Mass Curves and Its Applications in Hydrology and Meteorology[J].Journal of China Hydrology,2010,30(4):4751.(in Chinese)) |
随便看 |
|
科学优质学术资源、百科知识分享平台,免费提供知识科普、生活经验分享、中外学术论文、各类范文、学术文献、教学资料、学术期刊、会议、报纸、杂志、工具书等各类资源检索、在线阅读和软件app下载服务。