水利工程论文-朗梨站洪峰水位多维混合回归预报模型.doc_第1页
水利工程论文-朗梨站洪峰水位多维混合回归预报模型.doc_第2页
水利工程论文-朗梨站洪峰水位多维混合回归预报模型.doc_第3页
水利工程论文-朗梨站洪峰水位多维混合回归预报模型.doc_第4页
水利工程论文-朗梨站洪峰水位多维混合回归预报模型.doc_第5页
已阅读5页,还剩1页未读 继续免费阅读

下载本文档

版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领

文档简介

水利工程论文-朗梨站洪峰水位多维混合回归预报模型摘要:根据浏渭河流域水文特性对多维混合回归系统模型的结构进行了改进,并在系统输入非线性化处理基础上建立了朗梨水文站洪峰水位多维混合回归系统预报模型,经检验评定,效果良好。关键词:洪峰水位预报多维混合回归预报模型渭河流域0前言流域水文系统因其输入输出和系统输送既受确定性因素的作用又受随机性因素的影响,因此变化十分复杂,其精确描述显然是十分困难的。在当前一种既考虑水文过程形成的物理背景又不苛求预报量与预报因子间确切的运行机制的系统建模思想比较符合人们对客观水文规律认识的现实,因而得到了众多学者和专家的认可15。多维混合回归系统模型是将回归与自回归有机结合起来的混合回归模型,它既考虑了预报量自身的演变规律又充分利用了预报量对预报因子的依赖关系,预报模型具有一定的物理基础又不过分强调每一个参数的物理意义。其建模方式灵活,使用方便,精度较高6。1多维混合回归系统水文模型的结构和参数识别1.1模型结构回归模型是利用因变量与自变量的统计相关关系建立起来的,并可预报因变量的未来值。自回归模型则是利用自变量自身取值的相互依赖的统计关系建立起来的,并可用它预报时间序列的未来值。前者没有利用因变量自身变化的统计规律,而后者没有利用其它自变量对因变量的统计依赖关系,成因概念不十分明确。将回归和自回归结合起来的混合模型,取两家之长,补两家之短,在理论和具体方法上没有实质性的障碍,但在实际应用上却可以拓宽范围,丰富模型种类1。具体到流域水文系统,设其为多输入、单输出线性系统,系统的输出为,输入为Z,y1,y2,ys,则经过流域输送作用后,系统的混合回归描述可表示为式(1)即为多维混合回归系统水文的结构。式中y1,t,y2,t,ys,t分别为系统输入y1,y2,ys在t时刻的数值,Zt为系统输入Zt在t时刻的数值,b(0),b(1),b(2),b(s)为系统响应函数。1.2参数识别式(1)中包含了s+1个变量,是多维的。参数共有(p0+p1+p2+ps+1)个。实际上p0为系统自回归阶数,p1,p2,ps为因子y1,y2,ys的阶数。模型阶数和参数识别算法如下。设已知实测数据Zt,y1,t,y2,t,ys,t,其中t=1,2,m。令p=maxp0,p1,ps,在式(1)中分别取t=p+1,p+2,m,由此可得m-p个等式,其缩写形式为式中根据最小二乘原理,B的最小二乘估计为相应的残差平方和表达式为可按赤池信息准则(AIC准则)确定阶数p0,p1,ps。AIC准则函数为对给定某个备选的最大延迟量时,在所有可能的符合以下不等式的值中找出使式(5)中的AIC最小的p0,p1,ps。0p0,p1,psp(6)如此选取出的p0,p1,ps即为式(1)中各变量的阶数。当s较大时,参数个数很多,识别计算工作量很大。实际上参数估算和阶数确定是相互联系,互为前提的,为减少计算工作量,提高模型的统计稳定性,实际应用时可采用疏系数方法求解6。2朗梨站洪峰水位多维混合回归预报模型2.1流域基本情况朗梨站系湘江一级支流浏渭河的总控制站,控制流域面积3815km2,干流长度110km,平均坡降0.27。该站为平原地形,上下游附近两岸均为堤垸,上游15km浏阳市境内为丘陵地形,河谷仍较开阔,河槽发育。浏阳双江口以上为山间盆地,山势平缓,有二级阶地,在双江口设有双江口水文站,控制流域面积2067km2,区间有江背报汛雨量站。浏渭河在朗梨站下游约25km处汇入湘江。长沙站位于湘江干流浏渭河汇合处上游约2km处。朗梨站洪峰主要由上游流域的暴雨形成,同时也受湘江洪水顶托。但浏渭河洪峰一般先于湘江干流洪峰到达长沙,故朗梨站洪水历时较长,全年最高水位多发生在6-7月。由于洪水波运动的非线性特征加之下游严重的回水顶托影响,使得朗梨站的水位流量关系异常复杂,相似的上游来水和区间降水在朗梨站形成的洪峰水位相距甚远。用传统的上下游水位相关法建立洪峰水位预报方案合格率不高,方案不能正常投入运行,作业预报的难度很大。为此考虑建立该站洪峰水位多维混合回归预报模型。2.2模型结构的改进根据朗梨站洪峰水位预报的实际需要,系统以朗梨站洪峰水位为输出,以双江口站出峰前相关站的水位和区间降水为输入。其混合回归预报模型的结构在式(1)的基础上略加改进,结构形式为式中朗Zt为t时刻朗梨站的水位;为洪峰从双江口至朗梨站的传播时间,也是本模型的预见期;双Zm,t为t时刻双江口站的水位,下标加注m表示其值同时为洪峰水位;长Zt为t时刻长沙站的水位;为(t-1)t时刻双江口站的降水量非线性化因子;为(t-1)t时刻江背站的降水量非线性化因子;si为第i个河段的水力特性因子;bj(k),ai为系统响应。2.3系统输入的非线性化处理由于降水径流和河段水力特性的非线性特征,线性系统的模拟能力有限,对系统的降水和河段水力特性输入因子均作非线性处理。即:根据现行概念性降雨径流模型,降雨径流关系均为非线性关系,且接近指数型曲线,故对降水输入取=p,为简便起见分别取0.1,0.2,0.3,2.0后,按AIC准则确定采用值;对河段水力特性因子主要考虑双江口站出峰时双江口朗梨,朗梨长沙的水面落差,由曼宁公式不难推知河道流量与水面比降的平方根成正比,故模型固定取:其中和为基面换算系数,以保证s1和s2大于零。2.4模型的建立因式(7)较之式(1)有较大改进,在资料处理上已将连续水文系统截断为以单次洪峰为单位的多个序列资料,以之分别代入式(7)采用疏系数方法并结合AIC准则,即可得到最优预报模型。采用朗梨站1965,19731995年共34次洪峰资料,并取近30小时的所有有效信息即给定最大延迟量p=5,得到洪峰水位最优预报模型为:朗Zm,t+=45.76040.6304朗Zt0.6913朗Zt-1+0.2188朗Zt-2+1.8674双Zm,t+0.2970双Zt-1-0.5683双Zt-2+0.0007双Zt-3+0.0409长Zt+0.3240长Zt-1+0.1854长Zt-2+0.0588双+0.0458双+0.0713江-7.8302s1+0.5204s2式中:峰现时间的预报采用同期资料建立双江口站洪峰水位与洪峰传播历时相关图,即(双Zm,t)经验相关图查得。限于篇幅(双Zm,t)关系图从略。2.5模型的检验与评价根据水文情报预报规范7,模型的评定指标主要有合格率和有效性两项。其中:合格率=合格场次/总场次100%模型率定期内34场洪水洪峰水位变幅6.37m,模型拟合均方差为1.61m,评定时允许误差最大值取1.00m,最小值取0.30m,单次洪水洪峰水位允许预报误差据预见期水位变化幅度(i)在0.301.00m间内插确定,允许误差为允许=0.30+(0.70i)/6.37作业预报的评定按每次预报误差的大小分四个等级,即:优级:预报误差/许可误差25%;良级:25%预报误差/许可误差50%;合格:50%预报误差/许可误差合格率70%,0.91效率系数0.71;丙等方案:70%合格率60%,0.71效率系数0.50.对照上表,本方案可综合评定为甲等方案。2.6实际预报作业模型用于1997年至1999年四次洪水预报,预报结果列于表

温馨提示

  • 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

评论

0/150

提交评论