朗梨站洪峰水位多维混合回归预报模型Word格式文档下载.docx
《朗梨站洪峰水位多维混合回归预报模型Word格式文档下载.docx》由会员分享,可在线阅读,更多相关《朗梨站洪峰水位多维混合回归预报模型Word格式文档下载.docx(3页珍藏版)》请在冰豆网上搜索。
前者没有利用因变量自身变化的统计规律,而后者没有利用其它自变量对因变量的统计依赖关系,成因概念不十分明确。
将回归和自回归结合起来的混合模型,取两家之长,补两家之短,在理论和具体方法上没有实质性的障碍,但在实际应用上却可以拓宽范围,丰富模型种类[1]。
具体到流域水文系统,设其为多输入、单输出线性系统,系统的输出为,输入为{Z},{y1},{y2},…,{ys},则经过流域输送作用后,系统的混合回归描述可表示为式即为多维混合回归系统水文的结构。
式中y1,t,y2,t,…,ys,t分别为系统输入{y1},{y2},…,{ys}在t时刻的数值,Zt为系统输入{Zt}在t时刻的数值,{b(0)},{b
(1)},{b
(2)},…,{b(s)}为系统响应函数。
参数识别式中包含了s+1个变量,是多维的。
参数共有个。
实际上p0为系统自回归阶数,p1,p2,…,ps为因子{y1},{y2},…,{ys}的阶数。
模型阶数和参数识别算法如下。
设已知实测数据{Zt},{y1,t},{y2,t},…,{ys,t},其中t=1,2,…,m。
令p=max{p0,p1,…,ps},在式中分别取t=p+1,p+2,…,m,由此可得m-p个等式,其缩写形式为式中根据最小二乘原理,B的最小二乘估计为相应的残差平方和表达式为可按赤池信息准则确定阶数p0,p1,…,ps。
AIC准则函数为对给定某个备选的最大延迟量时,在所有可能的符合以下不等式的值中找出使式中的AIC最小的p0,p1,…,ps。
0≤p0,p1,…,ps≤p如此选取出的p0,p1,…,ps即为式中各变量的阶数。
当s较大时,参数个数很多,识别计算工作量很大。
实际上参数估算和阶数确定是相互联系,互为前提的,为减少计算工作量,提高模型的统计稳定性,实际应用时可采用疏系数方法求解。
2朗梨站洪峰水位多维混合回归预报模型流域基本情况朗梨站系湘江一级支流浏渭河的总控制站,控制流域面积3815km2,干流长度110km,平均坡降‰。
该站为平原地形,上下游附近两岸均为堤垸,上游15km浏阳市境内为丘陵地形,河谷仍较开阔,河槽发育。
浏阳双江口以上为山间盆地,山势平缓,有二级阶地,在双江口设有双江口水文站,控制流域面积2067km2,区间有江背报汛雨量站。
浏渭河在朗梨站下游约25km处汇入湘江。
长沙站位于湘江干流浏渭河汇合处上游约2km处。
朗梨站洪峰主要由上游流域的暴雨形成,同时也受湘江洪水顶托。
但浏渭河洪峰一般先于湘江干流洪峰到达长沙,故朗梨站洪水历时较长,全年最高水位多发生在6-7月。
由于洪水波运动的非线性特征加之下游严重的回水顶托影响,使得朗梨站的水位流量关系异常复杂,相似的上游来水和区间降水在朗梨站形成的洪峰水位相距甚远。
用传统的上下游水位相关法建立洪峰水位预报方案合格率不高,方案不能正常投入运行,作业预报的难度很大。
为此考虑建立该站洪峰水位多维混合回归预报模型。
模型结构的改进根据朗梨站洪峰水位预报的实际需要,系统以朗梨站洪峰水位为输出,以双江口站出峰前相关站的水位和区间降水为输入。
其混合回归预报模型的结构在式的基础上略加改进,结构形式为式中朗Zt为t时刻朗梨站的水位;
τ为洪峰从双江口至朗梨站的传播时间,也是本模型的预见期;
双Zm,t为t时刻双江口站的水位,下标加注m表示其值同时为洪峰水位;
长Zt为t时刻长沙站的水位;
为(t-1)~t时刻双江口站的降水量非线性化因子;
为(t-1)~t时刻江背站的降水量非线性化因子;
si为第i个河段的水力特性因子;
bj(k),ai为系统响应。
系统输入的非线性化处理由于降水径流和河段水力特性的非线性特征,线性系统的模拟能力有限,对系统的降水和河段水力特性输入因子均作非线性处理。
即根据现行概念性降雨径流模型,降雨径流关系均为非线性关系,且接近指数型曲线,故对降水输入取=pα,为简便起见α分别取,,,…,后,按AIC准则确定采用值;
对河段水力特性因子主要考虑双江口站出峰时双江口~朗梨,朗梨~长沙的水面落差,由曼宁公式不难推知河道流量与水面比降的平方根成正比,故模型固定取其中α和β为基面换算系数,以保证s1和s2大于零。
模型的建立因式较之式有较大改进,在资料处理上已将连续水文系统截断为以单次洪峰为单位的多个序列资料,以之分别代入式采用疏系数方法并结合AIC准则,即可得到最优预报模型。
采用朗梨站1965,1973~1995年共34次洪峰资料,并取近30小时的所有有效信息即给定最大延迟量p=5,得到洪峰水位最优预报模型为朗Zm,t+τ= –朗Zt– 朗Zt-1+朗Zt-2+双Zm,t+双Zt-1-双Zt-2+双Zt-3+长Zt+长Zt-1+长Zt-2+双+双+江- s1+s2式中:
峰现时间的预报采用同期资料建立双江口站洪峰水位与洪峰传播历时相关图,即经验相关图查得。
限于篇幅关系图从略。
模型的检验与评价根据《水文情报预报规范》,模型的评定指标主要有合格率和有效性两项。
其中合格率=合格场次/总场次×
100%模型率定期内34场洪水洪峰水位变幅,模型拟合均方差σ△为,评定时允许误差最大值取,最小值取,单次洪水洪峰水位允许预报误差据预见期水位变化幅度在~间内插确定,允许误差为△允许= +(△i)/作业预报的评定按每次预报误差的大小分四个等级,即优级:
预报误差/许可误差25%;
良级:
25%≤预报误差/许可误差50%;
合格:
50%≤预报误差/许可误差100%;
不合格:
预报误差/许可误差≥100%。
模型的有效性评定采用效率系数衡量,效率系数公式为其中;
为预报误差;
△i为预见期内水位变幅;
为△i的均值;
n为点据数。
以1966~1971年共10场洪水对模型作验证预报。
率定期34场洪水回顾预报和验证期10场洪水的检验预报精度的综合评价详见表1。
表1模型预报精度的综合评价单次预报误差幅度确定性系数单次预报优等率单次预报良好率单次预报合格率资料分段|△|min~|△|maxd△(%)(%)(%)率定期0~ 验证期~ 按《水文情报预报规范》预报方案的等级划分标准和我国现行通用的水文预报精度评定方法:
甲等方案:
合格率≥85%,效率系数≥;
乙等方案:
85%合格率≥70%,效率系数≥;
丙等方案:
70%合格率≥60%,效率系数≥对照上表,本方案可综合评定为甲等方案。
实际预报作业模型用于1997年至1999年四次洪水预报,预报结果列于表2。
表2模型作业预报精度评定洪号实测洪峰水位预报洪峰水位(m)预报误差(m)允许误差(m)评定970608 - 优秀980614 良好980627 优秀990718 - 良好四次预报二优二良,可见预报效果是令人满意的。
3存在问题和解决办法流域区间有官庄、友谊、乌川、东庄等水利工程,由于报汛站网布局不尽合理,难于收集水库调蓄和开闸泄洪的信息,对部份洪峰水位的预报精度造成一定影响,有待报汛站网的补充和调整。
在双江口出现洪峰后,如流域内特别是流域上游继续降大到暴雨,本模型对此类洪水洪峰水位的预报精度不高。
对这类洪水应考虑过程预报方案。
鉴于本模型使用的水文资料序列不长,很难完全反映浏渭河流域的水文系统特征。
用信息论的观点看,属于信息不完全系统。
在正式投产预报作业时,应及时补充新的水文资料,重新识别模型参数,以求尽可能地接受信息,以包容各种水情状况,增强模拟和实际预报效果。
参考文献[1].李纪人,刘德平.水文时间序列模型及预报方法[M].南京:
河海大学出版社,~226.李正最,吴雅琴.混合门限自回归模型在长期水文预报中的应用[J].江西水利科技,1993,19:
309~314.李正最,吴雅琴.非方程灰色系统在长期水文预报中的应用初探[J].四川水力发电,1994,13
(1):
16~22.徐向阳.具有物理基础的太湖水位多元回归模型[J].河海大学学报,1990,18:
48~54.葛守西.现代洪水预报技术[M].北京:
中国水利水电出版社,1999:
30~129.吴雅琴.多维混合回归系统水文预报模型及其应用[A].迈向新世纪[C]。
长沙:
湖南科学技术出版社,1999:
251~255.中华人民共和国行业标准.水文情报预报规范[S].北京:
中国水利水电出版社,~14.长江水利委员会.水文预报方法[M].北京:
水利电力出版社,~362AMultidimensionalMixedRegressionModeltoForecastPeakStageatLangliStationAbstract:
Themultidimensionalmixedregressionmodelisacomprehensivemodelcombinedwithregressionandautoregression.AccordingtothehydrologiccharacteristicsofLiuweiheRiverBasin,thispaperestablishesmultidimensionalmixedregressionmodeltoforecastpeakstageatLanglihydrometricstation.Thetestforecastresultsaresatisfactory.Keywords:
forecastingofpeakstage;
multidimensionalmixedregressionforecastmodel;
LiuweiheRiverBasin