水利工程论文-长江中游散乱水位流量关系单值化处理方法研究.doc

水利工程论文-长江中游散乱水位流量关系单值化处理方法研究摘要：本研究在重点分析荆江洞庭湖防洪重要水文站监利、城陵矶(七里山)、螺山、武汉关水位流量长系列实测资料的基础上，初步建立了考虑河段水面比降因子模拟水位流量的单一关系模式：即J=J0+JQ。该式待定系数意义明确，简便易求，并可借此对河道水流特性进行深入分析，利用该方法所求结果与实测值完全吻合。关键词：长江水位单值化1问题的提出长江荆江河段水位流量关系变化既受荆江裁弯河床演变、洞庭湖分流入汇变化、上游葛洲坝兴建，又将受未来三峡水库建库的影响，变得十分复杂。荆江包括洞庭湖区行洪水道，是长江中游防洪最严峻的河段，水位流量的变化，尤其是城陵矶(洞庭湖出湖七里山站)、监利水位流量关系是荆江与洞庭湖防洪研究的重要特征参数。近来研究水位流量关系特性采用的方法主要有：(1)综合落差指数法：对城陵矶七里山站流量进行单值化处理，公式形式为Q=Qm/H=f(H)(1)H=1H1+2H2+3H3(2)式中Qm为实测流量，H为七里山站水位，H1、H2、H3分别为湖区鹿角至七里山、七里山至上游监利、七里山至下游螺山的落差，1、2、3为相应的距离(km)，为待定指数。选取典型年作各大时段单一水位流量关系。(2)幂函数拟合法：依据裁弯前后的一部分水文资料，分别建立在某一顶托流量下七里山站、监利站水位流量关系，采用幂函数拟合，公式形式为Z=Q(3)式中Z为水位，Q研究者认为，上述方法一般相关关系较好，除个别点外，相关系数可达到0.9以上。本文研究对相应河段乃至下游螺山站、武汉关站水位流量关系进行了深入分析，提出关系模式成果，将克服上述两个方法不足。并(1)为仅需一站水位资料。(2)变分级流量关系式为统一的关系式。2分析方法及基本关系式确定2.1水位随流量变化特性点绘监利站1980年至1987年逐日水文关系，如图1所示。由图可知，监利同流量下水位变幅最大可达45m,这种散乱关系，是受下游洞庭湖出流影响的结果，同流量条件下，下游水位高，该站水位就高，反之亦然。既使在关系点相对集中的流量范围，如15000m3/s流量以下，流量由4000m3/s增加到13000m3/s时，水位相应由21.2m抬高至27.8m，变化幅度达6.6m。这种变化的主要因素也是下游水位抬高引起的。下游水位变化小时，流量的变化引起上游水位抬高的变化要小的多，以下游螺山水位2021m为例，流量变化在3960m3/s至12900m3/s，监利水位由23.4m增至25.8m，变化仅2.4m。如将下游螺山站的水位按大小划分若干组，则可知每一组的监利Q、H均存在着图2所示特性。图11980至1987年监利站逐日水位流量关系RelationbetweenwaterlevelanddischargeatJianlistationfrom1980to1987图2监利站水位流量关系(H螺山=2021m)RelationbetweenwaterlevelanddischargeatJianlistation(Hluoshan=2021m)分析认为，监利的水位由以下两部分组成，一是图2中的点延长线至Q=0的交点，记为H0；二是因流量的变化引起的水位变化，记为HQ，即H上=HQ+H0(4)2.2HQ、H0的确定将监利站统计年份逐日流量与监利至螺山水面比降关系点绘在图3中，点群是很散乱的，同样采取下游水位分组方法，如H下=2021m组点，见图4“+”点所示。H下=2728m组的点也一并点绘于图4中，用“*”点表示。图3监利流量与监利至螺山水面比降关系RelationbetweendischargeatJianlistationandwatersurfaceslopeinreachfromJianlitoLuoshan图4不同螺山水位监利流量与水面比降关系RelationbetweendischargeandwatersurfaceslopeatJianlistationunderdifferentwaterlevelsatLuoshanstation3中很散乱的点是由下游水位不同组数据汇合而成的，实际每组的QJ关系是很有规律的。同组条件下，比降随流量的增大而单调增加。进一步分析可知，J的增加与Q可由线性关系表示。如将Q=0交点的比降值用J0表示，流量引起的比降变化用JQ表示，即可得到J=JQ+J0(5)对比(4)式可知，监利水位HQ、H0分别可表示为JQL、J0L+H下(L为两站间距，H下为下站水位)。即有H上=JQL+J0L+H下(6)2.3经验拟合J与QJQ=(H下)Q(7)分析数站实测水文资料，研究后发现(H下)、J表达式分别为(H下)=10aH下b+c(8)J0=d+eH下(9)将(8)、(9)式代入(5)、(6)式中，即为本文建立的水位流量关系基本表达式H=(10aH下b+c)QL+(d+eH下)L+H下(10)或J=(10aH下b+c)Q+(d+eH下)(11)其中a、b、c、d、e五个为待定系数，且a、b、c与d、e数均有明确的数学含义和其简单的计算方法，后文中将详细介绍。其他符号意义同前。3实例计算3.1系数确定水文资料按监利19801987年下游螺山站水位大小分组，计算用两次最小二乘法，分别对每组数据点拟合求J、(H下)，结果点绘于图5及图6中。图5螺山水位与监利J关系RelationbetweenwaterlevelatLuoshanandJ0atJianlistation图6螺山水位与监利(H下)关系RelationbetweenwaterlevelatLuoshanand(H下)atJianlistation由上述图点关系可分别求得a、b、c、d、e(比降J用1/1000Q用103值，距离L用公里数)。系数值列于表1。以同样方法，对七里山站(下站螺山)、螺山站(下站龙口)、武汉关站(下站黄石)1980年1987年逐日水文资料进行分析，确定各站系数并一一列入表1中。表1系数取值(19801987)Valuesofcoefficient(19801987)站名监利(螺山)七里山(螺山)螺山(龙口)武汉关(黄石)a4.223.76344.6267-0.85460a2-5.333-5.2259-6.1998-1.9122a30.0003080.0000433040.00021897-0.000026073a4-0.00061278-0.00071889-0.000621140.00055556a50.0316780.0496500.0286960.00066673.2验算结果将统计年份的点计算对比结果点绘于图7(a、b、c、d)中，可见，计算与实测点完全吻合。综上分析，得出以下认识，(1)无论水文站流量水位点分散或单一，该站水位高低主要受下游站水位变化的影响；(2)利用上述(10)式关系，当已知本站流量、水位、下站水位三个参数中的任意两个，均可计图7实测水位与计算值对比Comparisonofobservedandcomputedstages图8监利站实测流量与计算值对比ComparisonofobservedandcomputeddischargeatJianlistation算出第三个参数值。仍以监利站为例，由监利站水位及下游水位计算得到监利流量,并与实测流量对比，点绘于图8中，点子落在45线上；(3)天然河道均匀流概念很难定义，而由均匀流基本公式J=Q2n2/A2R4/3；仍看不出公式中各参数之间的因果关系。实际上，Q、n、J均是独立量。如果边壁条件不一样，水深或水位变化将引起n值的变化，反映河道沿程阻力特性，可认为流量变化仅引起水深(水位)变化，但对n值不应有所影响。本文推荐的水位流量公式中包括与流量无关的比降项，即J0(=d+eh)。天然河道中均匀流时，Q、n值无法确定，仍有待今后研究，但该项能够具体计算，确实对该问题的解决大大前进一步，这方面研究有待进一步深入。4监利水位年代变化分析建立了上述水位流量关系式，即可很方便地对比分析不同年代各种水文组合条件下的变化特点。4.1水文频率分析受实测资料所限，与80年代相对比时段选为19541959年系列，对监利站水位流量及下游螺山站水位进行频率分析，结果列于表2。由表中数据可看出，50年代与80年代相比，同频率监利流量增大，水位反而有所下降。其原因之一是螺山水位同频率水位也是在下降的。4.2监利流量同频率条件水位变化对于表2中监利各频率流量在不同下游水位条件下用水位流量关系计算出本站水位变化，如图9(a,b),图10(a,b)所示，由图可知：(1)不同年代出现最大流量时，下游同水位条件下80年代水位具有高于50年代水位的趋势，但实际上因下游螺山水位在80年代低于50年代水位，因此监利水位在80年代低于50年代。表2监利流量水位、螺山水位频率分析对比AnalysesoffrequencyofdischargeandwaterlevelatJianliandthewaterlevelatLuoshan站年代频率Max1%5%10%50%min监利19801987Q3520029400232001880078103140H34.733.8932.3131.2625.8722.0819541959Q4570037900292002600090803150H34.5133.3631.9531.2026.0021.16螺山1980198731.3231.2728.8427.9021.2813.991954195930.9829.8128.1027.3121.8915.32图9监利站80年代与50年代水位变化趋势TendencyofstagesatJianlistationin80sand50s图100.10频率流量下监利站水位变化趋势与实测值对比ComparisonofcomputedandobservedstagesatJianlistationataflowfrequencyof0.10(2)其它频率时，监利水位变化参见图9(b)及图10(a，b)。80年代为27.028.0m,相应监利水位33.5m，50年代下游水位31m左右，监利水位达到34.6m,高出1.1m。由图10(a,b)可知，将实测点绘