城市水文学之由流量资料推求设计洪水

,Click to edit Master text styles,Second Level,Third Level,Fourth Level,Fifth Level,Click to edit Master title style,城市水文学-8,由流量资料推求设计洪水,由流量资料推求设计洪水,第一节 概述,一、问题的提出,在河流上兴建水库，目的在于兴利除害。为了兴利（灌溉、发电等），需要设置一定的兴利库容，调节年、月径流，使之符合人们的要求。,但水库单单有兴利库容是否就行了呢？,为水利水库自身安全和下游防护区的安全，还必须设置一定的库容拦蓄洪水。,设计洪水,拦洪库容,设计洪水位；,校核洪水,调洪库容,校核洪水位,;,水库泄洪,泄洪建筑物；,死水位,Z,死,和死库容,V,死,；正常蓄水位,Z,蓄,和兴利库容,V,兴,；,防洪限制水位,Z,限,和结合库容,V,结,；防洪高水位,Z,防,和防洪库容,V,防,；,设计洪水位,Z,设,和拦洪库容,V,拦,；校核洪水位,Z,校,与调洪库容,V,调,；,水库总库容：,V,总,= V,死,+ V,兴,+ V,调,- V,结,三峡工程，正常蓄水位,175m,，防洪限制水位,145m,，枯季消落最低水位,155m,，,100,年一遇洪水位,166.9m,，设计洪水位（,1000,年一遇）,175m,，校核洪水位,180.4m,，坝顶高程,185m,。总库容,393,亿,m,3,（,175m,以下），兴利库容,165m,3,，防洪库容,221.5m,3,，水库库面面积,1084km,2,。,二、设计洪水的涵义和设计标准,水库防洪设计的依据。,两类水库防洪问题：,1.,水库本身安全防洪问题，是确定在某一特大,Q,t,情况下，为了不使洪水漫溢坝顶造成毁坝灾害，所需要的坝顶高程等工程规模数据。,如何设计调洪库容和泄洪建筑物？,水工建筑物的设计洪水,2.,下游地区防洪问题，一般是水库下游河道要求水库下泄流量不超过某一流量值。,如何设计防洪库容？,防护对象的设计洪水。,设计洪水定义,：为解决各类防洪问题，所提供的作为规划设计依据的各种设计标准的洪水。,如何选择水工建筑物的设计洪水，涉及一个标准问题，即,设计标准,。,设计标准定得,过高,，工程投资增大而不经济，但工程比较安全；,设计标准定得,过低,，工程造价降低，但工程遭受破坏的风险增大。,确定设计标准是一个非常复杂的问题。,我国：,SDJ12-78,水利水电枢纽工程等级划分及设,计标准（山区、丘陵区部分）（试行）,GB50201-94,防洪标准,水工建筑物的防洪标准：,正常运用标准,设计洪水：确定水库的设计洪水位、设计泄洪流量等。不超过这种标准的洪水来临时，水库枢纽一切工作维持正常状态。,非常运用标准,校核洪水：确定水库的校核洪水位。这种标准的洪水来临时，水库枢纽的某些正常工作可以暂时破坏，次要建筑物允许损毁，但主要建筑物必须确保安全。,防护对象的防洪标准：,三、设计洪水计算的基本方法和内容,1.,我国推求设计洪水的发展,（,1,）历史最大洪水加成法,以历史上发生过的最大洪水再加上一个安全值作为设计洪水。,缺点：, 对未来洪水超过历史最大洪水的可能性考虑不足，降低了工程的安全程度；, 对大小不同，重要性不同的工程采用同一个标准，显然不合理。,（,2,）频率计算法,以符合某一频率的洪水作为设计洪水，如百年一遇、千年一遇等。,此法将洪水作为随机事件，根据概率理论由已发生的洪水来推估未来可能发生的符合某一频率标准的洪水作为设计洪水。,此法克服了历史加成法存在的缺点，根据工程的重要性和工程规模选择不同的标准，适用面较宽，在我国水利、电力、交通设计中应用广泛。,（,3,）水文气象法,因频率计算缺乏成因概念，如果资料太短，用于推求稀遇洪水根据就很不足。且近年来，我国一再出现超标准的特大洪水，设计标准一再提高。水文气象法从物理成因入手，根据水文气象要素推求一个特定流域在现代气候条件下，可能发生的最大洪水作为设计洪水。,2.,设计洪水的内容,设计洪水三要素：,设计洪峰流量、设计洪量、设计洪水过程线,对于桥梁、涵洞、调节性能小的水库，一般可只推求设计洪峰流量，如葛洲坝电站，其泄洪闸以设计洪峰流量控制（,Qm=110000m,3,/s),。,对于大型水库，调节性能高，可以洪量控制，即库容大小主要由洪水总量决定。如三峡水库，拦洪库容,300.2,亿,m,3,。,一般水库都以峰和量同时控制。,第二节 设计洪峰流量及设计洪量的推求,一、资料审查,“,三性,”,审查： 可靠性、一致性、代表性,1.,资料可靠性的审查与改正,实测洪水资料：,对测验和整编进行检查，重点放在观测与整编质量较差的年份。包括水位观测、流量测验、水位流量关系等。,历史洪水资料：,一是调查计算的洪峰流量可靠性；二是审查洪水发生的年份的准确性。,2.,资料一致性的审查与还原,所谓洪水资料的一致性，就是产生各年洪水的流域产流和汇流条件在调查观测期中应基本相同。,如果发生了较大的变化，需要将变化后的资料还原到原先天然状态的基础上，以保证抽样的随机性（减少人为的干扰），和能与历史资料组成一个具有一致性的系列。,例如上游建了比较大的水库，则应把建库后的资料通过水库调洪计算，修正为未建库条件下的洪水。,3.,资料代表性的审查与插补延长,当洪水资料的频率分布能近似反映洪水的总体分布时，则认为具有代表性；否则，则认为缺乏代表性。实际工作中要求连续实测的洪水年数一般不少于,20,30,年，并有特大洪水加入。,当实测洪水资料缺乏代表性时，应插补延长和补充历史特大洪水，使之满足代表性的要求。插补延长主要是采用相关分析的方法。,二、样本选取,河流上一年内要发生多次洪水，每次洪水具有不同历时的流量变化过程，如何从历年洪水系列资料中选取表征洪水特征值的样本，是洪水频率计算的首要问题。,目前采用年,最大值法选样,：,即从资料中逐年选取一个最大流量和固定时段的最大洪水总量，组成洪峰流量和洪量系列。,固定时段一般采用,1,、,3,、,5,、,7,、,15,、,30,天。大流域，调洪能力大的工程，设计时段可以取得长些；小流域、调洪能力小的工程，可以取得短一些。,T=1,天,T=5,天,T=3,天,t(d),Q(m,3,/s),Q,m,W,1,W,5,W,3,在洪水资料审查中，样本的代表性要求洪水系列长,20,30,年，并有特大洪水加入。,那么下面主要讲什么是特大洪水、为什么要加入特大洪水、加入特大洪水进入后如何进行处理等问题,三、特大洪水的处理,1.,概述,(1),什么是特大洪水,?,特大洪水,是指实测系列和调查到的历史洪水中，比一般洪水大得多的稀遇洪水。,历史上的一般洪水是没有文字记载和留下洪水痕迹,只有特大洪水才有文献记载和洪水痕迹可供查证,所以调查到的历史洪水一般就是特大洪水,.,特大洪水可以发生在实测流量期间之内,也可以发生在实测流量期之外,前者称,资料内特大洪水,后者称,资料外特大洪水,(,历史特大洪水,).,历史调查期,实测期,Q,N,Q,N,资料内特大洪水,资料外特大洪水,(,历史特大洪水,),一般 时，,Q,N,可以考虑作为特大洪水处理。,历史调查期,实测期,（,2,） 特大洪水重现期,重现期,是指某随机变量的取值在长时期内平均多少年出现一次，又称多少年一遇。,要准确地定出特大洪水的重现期是相当困难的，目前，一般是根据历史洪水发生的年代来大致推估。,从发生年代至今为最大,N=,设计年份,-,调查期发生年份,+ 1,从调查考证的最远年份至今为最大,N=,设计年份,-,文献考证期最远年份,+ 1,例,1992,年长江重庆宜昌河段洪水调查,同治九年（,1870,年）川江发生特大洪水，沿江调查到石刻,91,处，推算得宜昌洪峰流量,Q,m,110000m,3,/s,。,N,n,1992,1870,Q,m,110000m,3,/s,如此洪水为,1870,年以来为最大，则,N=1992-1870+1,123,（年）。这么大的洪水平均,130,年就发生一次，可能性不大。,又经调查，在四川忠县长江北岸,2km,处的选溪山洞中调查到宋绍兴,23,年（南宋赵构年号）即,1153,年一次大洪水。,N,n,1992,1870,Q,m,110000m,3,/s,1153,该洪水小于,1870,年洪水，通过调查还可以肯定自,1153,年以来,1870,年洪水为最大，则,1870,年洪水的重现期为,N,1992-1153+1,840,（年）。,这样确定特大洪水的重现期具有相当大的不稳定性，要准确地确定重现期就要追溯到更远的年代，但追溯的年代愈远，河道情况与当前差别越大，记载愈不详尽，计算精度愈差，一般以明、清两代六百年为宜。,（,3,）为什么要考虑特大洪水？,目前我们所掌握的样本系列不长，系列愈短，抽样误差愈大，若用于推求千年一遇、万年一遇的稀遇洪水，根据就很不足。,如果能调查到,N,年（,Nn,）中的特大洪水，就相当于把,n,年资料展延到了,N,年，提高了系列的代表性，使计算结果更合理、准确。,例：河北省滹沱河黄壁庄水库设计洪水计算：,1955,年设计，资料,n=18,年，,Q,0.1%,=12600m,3,/s,；,1956,年发生特大洪水,Q=13100m,3,/s,，直接加入资料系列（,n=19,），未做特大洪水处理，,Q,0.1%,=25900m,3,/s,；,将,1956,年洪水做特大洪水处理，但不加历史特大洪水，,Q,0.1%,=19700m,3,/s,；,再加入历史特大洪水（,1794,、,1853,、,1917,、,1939,），,Q,0.1%,=,22600m,3,/s,；,1963,年又发生了一次特大洪水,Q=12000m,3,/s,，加入并做特大洪水处理，,Q,0.1%,=23300m,3,/s,。,由此可见加入特大洪水有助于提高样本的代表性和设计洪水的可靠性。但应注意的是，年代越久，由于河流演变等原因，推算的洪峰流量可能存在较大误差，必须尽可能的从多方面考察、论证。,2.,考虑特大洪水时经验频率的估算,加入特大洪水后，资料系列的特征：,（,1,）,连序系列和不连序系列：,缺测,所谓“连序”与“不连序”，不是指时间上连续与否，只是说所构成的样本中间有无空位。,所以特大洪水加入系列后，样本成为不连序系列，其经验频率和统计参数的计算与连序系列不同。这样就要研究有特大洪水时的经验频率和统计参数的计算方法，称为,特大洪水处理,。,考虑特大洪水时经验频率的计算基本上是采用将特大洪水的经验频率与一般洪水的经验频率分别计算的方法。,目前国内有两种计算特大洪水与一般洪水经验频率的方法：独立样本法、统一样本法。,设：,N,历史调查期年数：,n,实测系列的年数；,l,n,年中的特大洪水项数；,a,N,年中能够确定排位的特大洪水项数（含资 料内特大洪水,l,项）；,m,实测系列在,n,中由大到小排列的序号，,m,=,l,+1,，,l,+2,，,.,，,n,；,P,m,实测系列第,m,项的经验频率；,P,M,特大洪水第,M,序号的经验频率，,M,=1,2,.,a,N,n,a,项特大洪水,M=1,2,.,a,实测期内特大洪水，,l,项,.,.,T,Q(m,3,/s),.,.,实测一般洪水，,n-,l,项,m=,l,+1,l,+2,.,n,缺测,（,2,）独立样本法,把实测一般洪水系列与特大洪水系列都看作是从总体中独立抽出的两个随机连序样本，各项洪水可分别在各个系列中进行排位，,实测系列,的经验频率仍按连序系列经验频率公式计算：,特大洪水系列,的经验频率计算公式为：,当实测系列中含有特大洪水时，虽然这些特大洪水提到与历史特大洪水一起排序，但这些特大洪水亦应在实测系列中占序号，即实测系列的排序为,m=,l,+1,l,+2,.,n,。,（,2,）统一样本法,将实测系列与特大值系列共同组成一个不连序系列，作为代表总体的一个样本，不连序系列各项可在历史调查期,N,年内统一排位。,特大洪水,的经验频率仍采用下式,(,n-,l,),项实测,一般洪水,的经验频率计算公式为：,a,项特大洪水,M,=1,2,.,a,实测期内特大洪水，,l,项,.,.,P,Q(m,3,/s),.,P,m,P,M,实测一般洪水，,n-,l,项,m=,l,+1,l,+2,.,n,.,P,Ma,1-P,Ma,上述两种方法，我国目前都在使用。一般说，独立样本法把特大洪水与实测一般洪水视为相互独立，这在理论上有些不合理，但比较简单。在特大洪水排位可能有错漏时，因不互相影响，这方面讲则是比较合适的。当特大洪水排位比较准确时，理论上说，用统一样本法更好一些。,例,6-1,某站自,1935,1972,年的,38,年中，有,5,年因战争缺测，故实有洪水资料,33,年。其中,1949,年为最大，并考证应从实测系列中抽出作为特大值处理。另外，查明自,1903,年以来的,70,年间，为首的三次大洪水，其大小排位为,1921,、,1949,、,1903,年，并能判断在这,70,年间不会遗漏掉比,1903,年更大的洪水。同时，还调查到在,1903,年以前，还有三次大于,1921,年的特大洪水，其序位是,1867,、,1852,、,1832,年，但因年代久远，小于,1921,年洪水则无法查清。现按上述两种方法估算各项经验频率。,n=33,N,1,=70,N,2,=141,1921,1949,1903,1949,1921,1867,1852,1832,1972,1903,1832,1935,独立样本法：,1852,年,1832,年,1921,年,1867,年,调查期,N,2,=141,：,统一样本法：,同独立样本法,独立样本法：,1949,年,1903,年,1921,年 已抽到上面排序,调查期,N,1,=70,：,统一样本法：,独立样本法：,1940,年,1968,年,1949,年 已抽到上面排序,实测期,n=33,：,统一样本法：,.,.,.,3.,考虑特大洪水时统计参数的确定,（,1,）初步估计参数,矩法,假设系列中,n-l,年的一般洪水的均值为,x,n-l,、均方差为,n-l,，它们与除去特大洪水后的,N-a,年总的一般洪水系列的均值,x,N-a,、均方差,N-a,相等，即,:,x,j,特大洪水；,x,i,一般洪水,则可导出：,（,2,）三点法初估参数,如已知某曲线的数学方程，用选点法可确定方程中的参数。,例,直线方程,y=a+bx,在直线上任选两点，有：,y,1,=a+bx,1,y,2,=a+bx,2,解之可得,a,，,b,我们知道,P,-,型曲线有三个统计参数，,x,，,Cv,，,Cs,，于是在经验频率曲线上任选三点,:,根据,P,-,型曲线的性质有：,x,P,x,P,1,x,P,2,x,P,3,P,1,P,2,P,3,（,P,1,x,P1,),，,(,P,2,x,P2,),，,(,P,3,x,P3,),联解得：,令：,S,是,Cs,的函数，称偏度系数。计算时，可由计算的,S,值，查,S,Cs,关系表，求,Cs,。,再查,Cs,值表，得,(,P,1,Cs,),，,(,P,2,Cs,),，,(,P,3,Cs,),。,再联解下式,得；,三点的选取一般为：,1,50,99% 3,50,97%,5,50,95% 10,50,90%,例,6-2,已知某水库坝址处共有,1972,1954,年,19,年的实测洪峰流量资料。另外，通过历史洪水调查得知，,1922,年发生过一次大洪水，是,1922,年来最大，,1963,年洪水则为第二大洪水。试根据点绘的经验频率曲线，用三点法初估参数。,P(%),Q,m,(m,3,/s),2080,760,260,5%,50%,95%,选点：,Q,5%,=2080m,3,/s,Q,50%,= 760m,3,/s,Q,95%,= 260m,3,/s,计算,S,值：,由,S,查得：,由,Cs,查附表,1,值表，得：,计算得：,例题：,某水文站实测有,1938,年至,1992,年最大洪峰流量资料，其中最大的五年洪峰流量依次为,28400m,3,/s,，,13200m,3,/s,，,9850 m,3,/s,，,8560 m,3,/s,，,8450 m,3,/s,。另外，调查到,1927,年发生一次洪峰为,32000 m,3,/s,是,1856,年以来最大一次洪水，,1856,年至,1938,年间其余洪水的洪峰流量均在,15000 m,3,/s,以下，试用统一样本法计算上述六项洪峰流量的经验频率。,Cs,值：,对于,Cv,0.5,的地区，,Cs,=(3,4,）,Cv,；,对于0.5,1.0,的地区，,Cs,=(2,3,）,Cv,；,此外，还可以采用,权函数法,来估计,Cs,。,四、频率曲线线型选择,P,型频率曲线,五、频率曲线参数估计,1.,经验适线法（目估适线法）,1),尽量照顾点群的趋势,2,）曲线通过点群中心,3,）侧重考虑中上部的较大洪水点据,4,）对特大洪水作具体分析,2.,优化适线法,六、推求设计洪峰、洪量,（,1,）,在频率格纸上点绘经验点据；,（,2,）,选定频率分布线型（一般选用,P,型）；,（,3,）,参数估计,x,、,Cv,、,Cs,；,（,4,）,根据,x,、,Cv,、,Cs,，,查,附表,1,或,附表,2,，计算,x,p,值，以,x,p,为纵坐标，,p,为横坐标，即可得到频率曲线；,（,5,）,根据频率曲线与经验点据的配合情况，从中选择一条与经验点据配合较好的曲线作为采用曲线；,（,6,）,求指定频率的水文变量设计值。,取,Cv,=0.7,，令,Cs,=3,Cv,。,查附表,2,，得不同频率,P,的,Kp,值。则,Q,p,=,K,p,Q,或查附表,1,，,得不同频率,P,的,值，,Q,p,=,Q(1+Cv ),计算,Qp,。,为避免修改参数的盲目性，要了解参数 对频率曲线形状的影响：,a),值愈大，频率曲线位置愈高；,x,P,P,P,x,P,P,x,P,c),C,S,值愈大，频率曲线上段变陡，下段变缓，中部向左偏。,b),C,V,值愈大，频率曲线愈陡；,求,200,年一遇的洪峰流量？,P=0.5%,，,Cv=0.8,，,Cs=3.5Cv=2.8,查附表,2,，得,Kp=4.87,Q,0.5%,=4.87 587=2859m,3,/s,或,查附表,1,，得,=4.84,，代入,有,Q,0.5%,=587,（,1+0.84.84,）,= 2859m,3,/s,七、设计洪水估计值的抽样误差,频率计算是用有限样本估算总体的参数，必然存在误差。,统计参数的误差与所选的频率曲线线型有关。当总体分布为,P,型，对于,n,年连序序列，用矩法估计参数时，样本的均方误的计算公式为,：,均值的相对误差为：,设计洪水,x,p,的均方误近似公式为：,在实际工作中，对大型或重要的中型工程，应计算校核标准洪水的抽样误差。如果成果偏小，应加安全修正，一般不超过计算值的,20%,。,八、计算成果的合理性检验,（,1,）检查洪峰、各时段洪量的统计参数与历时之间的关系；,历时,增长，,均值,增大，,Cv,、,Cs,一般,减小。,P,Q W,7,d,5,d,3,d,(2),根据上下游、干支流及邻近地区各河流洪水频率分析成果进行比较。,从上游到下游，均值增大、模数减小,P,W,3,d,（下游站、干流站）,3,d,（上游站、支流站,）,（,3,）与暴雨频率分析结果进行比较,一般洪水径流深应小于相应天数的暴雨深，而洪水的,Cv,值大于相应暴雨量的,Cv,值。,第三节 设计洪水过程线的推求,目的：,调洪演算，以确定防洪建筑物的规模尺寸等。,定义：,设计洪水过程线是指具有某一标准的洪水过程线。,方法：,类似推求 设计年径流的年内分配，从实测资料中选取典型洪水过程线，按设计洪峰、洪量放大，即得设计洪水过程线。,同频率放大、同倍比放大,一、典型洪水过程线的选择,影响因素：,1,）洪水（尤其是特大洪水）的形成规律和天气条件；,2,）洪水过程，如大洪水出现的时间、季节、峰型（单峰、双峰或连续峰）、主峰位置、洪水上涨历时、洪量集中程度等。,原则：,1,）峰高量大，接近设计条件下的稀遇洪水；,2,）代表性（季节、地区组成、洪峰次数、峰量关系）；,3,）对过程不利（峰型集中、主峰靠后、与下游洪水遭遇）。,若符合条件的洪水过程可能不止一条，也可以选择几个典型过程，分别推求设计洪水过程线，供调洪计算时分析比较。,例如丹江口水库设计时，就曾选择了,1935,年和,1964,年洪水两个典型，前者是夏季洪水的典型，后者是秋季洪水的典型。,二、放大方法,1.,同频率放大（峰量同频率放大法）,洪峰放大倍比：,最大,1,天洪量放大倍比：,最大,3,天洪量除最大,1,天以外，其余两天的放大倍比：,特点：,放大后的过程线，其洪峰流量和各时段的洪量都符合同一设计频率。,注意：,1.,时段划分，不宜过多，一般以,3,段或,4,段为宜；,2.,对于放大后过程线的不连续现象，可徒手修匀，修匀后仍应保持洪峰和各时段洪量等于设计值。,2.,同倍比放大,1,）按洪峰控制的放大倍比：,2,）按洪量控制的放大倍比：,特点：,用同一放大倍比放大典型过程。,注意：,1.,用峰控制还是用量控制，要看峰、量哪个其主要作用；,2.,设计洪水过程线的峰或量偏离设计值。,“,以峰控制,”,，则洪峰等于设计值，洪量不一定等于设计值；,“,以量控制,”,，则时段洪量等于设计值，而洪峰不一定等于设计值。,第四节 分期设计洪水,目的：,为了满足不同时期的防洪需要；,定义：,一年中某个时段所拟定的设计洪水。,方法：,把整个年内划分为若干个分期，然后在分期的时段内，按年最大值法选样，进行频率计算。,一、洪水季节性变化规律分析和分期划分,依据：,设计流域洪水季节变化规律、工程要求,气候、均值、,Cv,、,Cs,的变化等,原则：,尽可能根据不同成因的洪水，把全年划分为若干分期。 分期不宜过短，一般以不短于,1,个月为宜。,选样：,跨期选样、不跨期选样。,跨期选样的日期不宜超过,5,10,日。,二、分期设计洪水的计算方法,1,）选样,2,）特大洪水处理；,经验频率计算；,参数估计；,配线；,3,）分期设计洪水过程线,施工初期围堰,设计洪峰流量,大坝合龙,设计洪峰和洪量,调洪,设计洪水工程线,4,）与全年最大洪水的峰量频率曲线对比,要点,特大洪水的确定与经验频率选择。,系列特征值的计算（均值、离差系数）,设计洪峰的计算（查表）,设计洪水过程（同频率、同倍比）,综合：了解由流量推求设计洪水的过程。,练习,1,经验频率推求：,某站具有,1952,年到,1983,年实测洪峰资料，其中,1955,年洪水最大，,Qm= 13600m,3,/s,。另调查到期,1878,年洪峰流量,Qm= 18200m,3,/s,，,1935,年的洪峰流量,Qm= 12500m,3,/s,，从,1878,年以来，调查考证可知未遗漏大于,1935,年洪水。,（,1,）、试分析各次特大洪水经验频率？,（,2,）如实测期次大洪水为,1967,年，洪峰,，问该洪峰流量如何排位？其经验频率为多少？,练习,2.,设计洪峰推,已求得某河流断面年最大洪峰流量频率计算结果为,Q=365 m,3,/s,、,Cv=0.72,，,Cs=3Cv,。试推求该河流断面,50,年一遇设计洪峰流量。,（查表：离均系数、模比系数）,