基于神经网络的水沙运动预报模型与回归模型比较及应用(李义天,李荣,黄伟)

热度136票 浏览172次 【共0条评论】【我要评论 时间:2001年1月01日 19:14

基于神经网络的水沙运动预报模型与回归模型比较及应用

|.\"G~0z6i6W0

李义天1,李 荣2, 黄 伟2
u6zV4Wl0(1.武汉大学 水沙科学教育部重点实验室,2.浦东新区水利

:]E [n1Q2|3UO0

摘要:河道水沙运动过程是一个复杂的非线性动力学过程,鉴于神经网络具有很强的处理大规模复杂非线性动力学系统的能力,本文将神经网络理论应用于河道水沙运动规律的模拟与预报,以期识别水流泥沙运动变化过程与其影响因子之间的复杂非线性关系,为河道水沙模拟和预报提供一条新的方法和途径。在此基础上将神经网络模型与传统回归模型进行了对比分析并用于螺山站流量大于50000m3/s的水沙规律模拟,两模型模拟结果与实际水沙规律一致。神经网络模型在水沙因素选择方面较回归模型简单,有较成熟的数学理论基础,而回归模型数学表达直观,易于理解。水利论文4{E.j b A:\iP2tb:x s'P

关键词:神经网络模型; 泥沙灾害; 回归模型水利论文.cO@E$T

基金项目国家自然科学基金委和水利部联合资助重大项目(59890200)

$g5P:w0W!\&]8\ pt0

作者简介:李义天(1957-),男,武汉大学教授。

8OxoF s?6B6rd0

1 引言水利论文hRi}.^0E;j}V+Q&t#CS

  90年代以来,我国长江、珠江、松花江、嫩江、辽河及淮河等洪水灾害此起彼伏,损失严重并成逐步加重趋势。特别1998年长江流域和松花江流域发生的特大洪水引起严重的洪水灾害,直接经济损失超过2000亿元。近年来的水灾加重与河道泥沙淤积加剧、蓄滞洪区围垸造地等有关。准确地揭示河道及湖区水流泥沙运动机理、进行水沙过程的高精度预报是提高防洪能力与减轻洪水灾害的有效途径。

!Z_ U w iu:cC0

  目前对河道水沙变化过程的研究方法主要有水动力学方法和水文学方法。水动力学方法从完整的圣维南方程组和泥沙运动基本方程出发,详细考虑水流泥沙沿河道传输中的各水力要素,计算精度高。但由于江湖关系非常复杂,缺乏完整的地形资料、水沙资料和区间产流资料,因此水动力学方法的应用还存在许多需要解决的问题。水文学方法是一种简化的计算方法,如相应水位(流量)法,流量(水位)时段涨差法,出流与槽蓄关系法等。这类方法经验半经验性较强,当河道地形发生改变和外界条件发生突然改变时,用以往资料建立的河段槽蓄量值函数关系不再适用,该方法的应用受到一定的限制[1]。河道水沙灾害形成机制非常复杂,受造成水沙变化过程的动力因子(如流量、上游来沙、河道地形和下游控制基准面等)的影响,这些因子的共同作用使得水沙演变过程成为一个复杂的非线性动力学过程。因此,寻求一种对河道地形资料要求不高又能准确地反映水沙运动及河床演变,特别是洪水演进变化规律的研究方法具有重要的理论意义和实用价值。

2{y5eL0_8v*q6p0

  人工神经网络是借鉴大脑和神经系统存储和处理信息的某些特征抽象出来的一种数学模型。目前应用最多,研究比较成熟的多层前馈网络误差反传算法模型,即BP模型是一种较特殊的非线性映射方法,它是通过一元函数的多次复合来逼近多元函数的映射方法。鉴于神经网络具有很强的处理大规模复杂非线性动力学系统的能力,神经网络理论在河道水流预报中得到初步应用[2,3],由于模型中忽略了河道冲淤变化对水流演进的影响,这类模型对河床冲淤变化不大的河道来说具有较高的预报精度;但对于河床冲淤变形比较大的河道,现有模型很难反映实际水流泥沙之间的相互作用情况,其长期预报精度普遍较低。为此,本文利用神经网络理论,考虑河床变化对河道水流运动的影响,在此基础上建立了河道水沙非线性动力学模型,以期识别河道水情变化过程与其影响因子之间的复杂非线性关系,为河道水沙运动变化和洪水预报提供一条新的途径。水利论文&CbLmMA0q/L"l

  多元回归拟合是传统的资料处理方法,通过给定拟合模型方程式的结构,利用资料率定公式系数,用于自然规律的预测预报。多元回归模型在水文、泥沙等方面已得到广泛的应用,具有成熟的理论和应用基础。

/V!O zd9r`)v?0

  人工神经网络的非线性映射方法与常用回归拟合方法在许多方面有着相似之处,但也有差别。基于神经网络理论的水沙预报模型(简称RFS-ANN模型)和回归拟合预报模型是河道水沙预报的两种不同的研究手段和方法。两种模型在模型原理预报精度等方面进行充分比较是今后模型选择的基本依据也是模型进一步研究发展的前提。

Q*} B-R*rRi0

2 模型建立与比较水利论文T1@8R:q8Y^

2.1 模型建立水利论文.@2C1iIz(u z

  对于复杂河道某一断面,其水沙规律一般受上游来水来沙条件,下游断面水沙控制条件、河道地形、分流分沙比等复杂条件的影响。写成表达式为水利论文.R fO5]I5r

水利论文IYaFB!x3W

(河道地)………

(1)

M&n \zg^?s0

式中φ代表所求断面的水沙变量;Qup表示上游来流条件;SupZdo为下游水位控制条件;λQ、λS分别表示河道的分流特性和分沙特性量;χ表示河道地形特性;ψ为非线性函数。水利论文 T,AW4P m ?

  RFS-ANN模型:一般来讲,ψ为复杂非线性函数,其表达式的确定存在相当大的困难或根本无法表达。而BP神经网络可以通过简单线性函数的复合来达到处理复杂非线性函数的目的。将影响水沙变量φ的影响要素Qup、Sup、Zdo、λQ、λS、χ以及其它条件作为网络输入节点变量,水沙变量φ作为网络输出节点变量,以三层网络结构为例建立水沙变量φ与影响要素之间的非线性映射神经网络模

水利论文+?dv(w7y`'G6I
图1 网络结构示意图
xl%[*`MnM0K,p1C1ql0Fig.1 Draft of net structure水利论文+S2|1p+FBN0[

型(如图1)。通过对历史资料的学习调整网络权重达到输入输出量之间的匹配,网络调整后的权重矩阵也就是非线性函数ψ的隐含形式。用调整后的网络就可以对未来水沙变量φ进行预测和分析。水利论文.d2D:Q6^ATa&p

  回归拟合模型:在详细分析因变量和各个自变量之间定性关系的前提下,确定非线性函数ψ的公式结构,通过历史资料确定函数中各项的系数,达到对历史水沙运动的正确模拟。水沙变量φ和影响要素Qup、Sup、Zdo、λQ、λS、χ等之间的回归模型可以写为

5o7W9S,ek0d.NI3X:O*X0

水利论文|f trNWQ

(2)

O6F S(C1A@F0

ni,j为任意的实数。上式表示水沙变量φ可以表达成影响要素Qup、Sup、Zdo、λQ、λS、χ等之间不同组合形式的表达。系数率定后的公式反映了水沙运动的基本规律,可以用于水沙运动的预测与分析。

B#DsGj:N0`0

2.2 模型原理比较

8O"fSO ]m(Nx0

  RFS-ANN模型和回归拟合模型都是基于对水沙规律的充分认识基础上,提取影响水沙运动的主要动力因子,建立影响因子与水沙因变量之间的关系,通过历史水沙资料来率定模型的方法。因此,影响因子选择的正确与否关系到两模型对实际水沙现象模拟的精度,历史资料的选择和资料的完整性也是决定两模型能否充分反映水沙运动规律、正确进行水沙变化预报的关键。尽管两种方法有很多的相似处,但RFS-ANN模型和回归拟合模型在建立模型方法和“学习”算法等方面存在着许多不同,主要在于:水利论文!VI'Y EM{#s

  (1) 模型形式不同:RFS-ANN模型是通过简单函数的复合来建立描述水沙运动的非线性映射函数,模型没有明显的数学表达式,其因变量和自变量之间的函数关系由神经网络的连接权重来隐含确定;而回归拟合模型有直观的函数关系式,因变量和自变量之间的函数关系是通过函数式显式表达出来的。水利论文mUF9ft+\Nu@3K

  (2) 学习方式差异:两种模型都存在历史资料的学习。RFS-ANN模型通过调整网络连接权重来满足计算输出量与实际要求输出量之间误差的要求,利用误差反传算法(即BP算法)来实现网络权重的自动调节;而回归拟合模型是利用最小二乘法来率定公式系数,最小二乘法与梯度下降算法比较相对简单些。

.k gNZa)};|*H2w4[0

  (3) 模型建立难易程度不同:当影响水沙运动的动力因子选择以后,影响因子就是RFS-ANN模型的输入节点变量,可以通过连接权重直接建立因变量和自变量之间的非线性关系,具有成熟的数学理论基础;而回归拟合不同,它需要考虑数学表达式中各项因子之间的组合和因子的指数,当选择不当,会造成较大的计算误差,回归拟合模型的主要困难是数学表达形式的很难选择,缺乏理论指导,经验和半经验性较强。水利论文J^ U2}%H

  (4) 模型容错性不同:容错性是神经网络的一个基本特征,网络某一连接权的误差对整体非线性映射的影响不大,而回归模型的某一项系数出现误差,对模拟精度将会造成较大误差。 (5) 计算速度上的差异:由于RFS-ANN模型利用的是混合算法(最速下降和共轭梯度结合算法),其模型率定速度与最小二乘法相比明显慢,但模型预测速度差别不很明显。

6d b _ny5NE}c0

3 模型应用比较

9k9V6R2L|XXz0

  本文以长江中游及洞庭湖为研究区域,将本文所研究的RFS-ANN模型和回归拟合模型应用于城陵矶站水位的模拟,对水位与各个影响因子之间的关系进行详细分析并与传统回归模型进行了比较。

\.s*v}ZA:[0

3.1 基本概况

.Sd%uJJg{S0

  洞庭湖位于长江中游荆江河段南岸,湖南省北部。北面有松滋、太平、藕池和调弦(1958年堵口)四口分泄长江来水,西南有湘、资、沅、澧四水入汇。洪水经洞庭湖调蓄后在城陵矶注入长江,与下荆江来水来沙汇合注入螺山汉口河段,构成了非常复杂的江湖关系。

1w`8O s9R#TA#J0

  近年来,长江中游特别是荆江和洞庭湖区水位普遍上涨,这种流量不大而水位抬高的根本原因是河道和湖区的不断淤积造成江湖关系的恶化。城陵矶作为洞庭湖的出口水文站近年的水位抬高尤其明显,城陵矶水位变化作为洞庭湖区水位变化的重要特征量,对其进行深入研究具有重要的意义。城陵矶水情变化受湖区来水、长江干流的顶托等因素的影响,直接研究城陵矶站水情变化存在一定的困难。已有资料表明,城陵矶站和螺山站的水位相关性比较好[5],而影响螺山站水情变化的因子相对比较简单,故我们可以用对螺山站水位的研究来代替对城陵矶站水位变化的分析。水利论文 \eW({7`$I5Jj3H

  螺山断面水位变化受螺山进口流量,下游汉口水位和螺山汉口河段地形变化的影响。在汉口断面水情稳定的条件下[4,5],同流量下螺山水位的抬高与汉口螺山段河床的多年累计淤积必然存在着一定的联系。尽管螺山至汉口段累积淤积是螺山、城陵矶、洞庭湖区和荆江河段水位抬高的根本原因,这一点已经基本形成共识[6-11],但定量研究螺山至汉口河段淤积与湖区和长江干流水位抬高之间的关系很少。到目前为止,研究螺山水位变化的手段主要有水动力学方法、水文学方法、河床演变分析等,尽管研究手段有差异,其目的都是通过对历史水流泥沙运动资料的分析,预测未来水流泥沙运动的变化。因此,提出符合水沙演变规律的非线性模型、研究河道淤积对水位抬高的影响对推动整个江湖关系的研究、减少自然灾害损失具有一定的理论意义和实用价值。螺山站水位~流量关系出现同流量下水位抬高的现象,一般水文学模型和不考虑河道淤积影响的模型很难正确反映这一规律。考虑到水力学模型对地形资料的严格要求和水文学模型对泥沙淤积变化影响的反应不足,本文试图通过多元回归方法和神经网络思想建立水沙预报模型对螺山至汉口段水流泥沙进行分析并对两模型进行对比。水利论文xcnob%F6^!u[

3.2 模型率定水利论文g.fsd5WK.u J+?wb%Z

  螺山水位RFS-ANN模型:影响螺山水位变化的主要因素有汉口水文站的水位、螺山至汉口河段的地形变化、螺山上游的来流情况等。由于每年都对河段进行地形测量或及时获得详细地形资料存在相当大的困难。在缺乏精确地形资料的条件下,本文用河段累积淤积量近似代替地形资料描述该河段地形的变化。水利论文tb FG!?}E,[

  建立螺山断面水位Z与该断面来流量Q、下游汉口站水位Z汉以及河段累积淤积量V等因素之间的关系,即

9C$H/V{W.bK @0

Z=f(Q,Z,V)

/_NKI&}.K0

(3)

%cr$[l.KZwd(T/ec w0

与一般神经网络模型相比,该模型,综合考虑各个影响因子对水位的影响,特别是模型考虑了河道变形V对上游水位的影响。

Iumd f1pq0

  BP网络的输入、输出层的结点数是根据实际问题的需要来决定的。输入层的节点数相应为非线性映射函数的自变量数。对螺山汉口河段来讲,影响螺山水位变化的主要动力因子有螺山站的来流量、汉口水位、河道地形等,因此模型中输入结点数取3(分别对应螺山站来流量Q、汉口站水位Z以及螺山汉口河段累积淤积量V),输出结点取1(对应螺山站水位Z)。中间层结点数的选取具有一定的任意性,目前在理论上还没有定论。研究表明网络的输入结点数与隐层结点数相等时,网络性能比较稳定[12,13],故本文模型网络结构为3-3-1型。

.Q-z ^9A.\ Q0NC p0
  建立人工神经网络预报模型,实质上是要获得反映训练样本数据中输入与输出映射关系的权值。如何获得预报权重是建立预报模型的关键所在。在选定网络结构(输入结点数、隐层数、输出结点数等)以及初始权重后,用大量样本训练神经网络,训练达到精度要求后,得到网络权重,以此网络来对未来进行预测。


Y5wQ?TCoS0bF0图2 网络结构示意图
2`-l)fJ%L&B0Fig.2 Draft of net structure
水利论文]lW ~E N \aS

  选取螺山水文站1955~1992年38年汛期50000ms/s以上流量系列训练网络。为了消除各因子由于量纲和单位不同的影响,对于输入样本进行归一化处理,使比较大的输入仍落在神经元转换函数梯度最大的区域。归一化公式如下

7r(w9j$K6wG.@Qq0

%qOj9M"Z0

(4)水利论文nD}.g _ @jF

式中 Ximax、Ximin分别表示第i个影响因子在学习样本资料中的最大、最小值,Xi表示归一化后的影响因子值。经过上述变换后,影响因子都转换到[0,1]范围内。

!boqI$s?S{ t0

  为了防止神经元达到过饱和,将输出状态作如下变换水利论文8eoY}gaA

{oS%H%{'e0`KE0

(5)水利论文V$l8Uf"n B\{2]

式中 yi,yimax,yimin分别表示因变量向量第i个分量及其最大、最小值;d1,d2分别取0.998和0.001。

7d8FyW~7v*x}8G/V&l0

  螺山水位变化回归拟合公式:文献[14]建立了螺山水位(Zl)与汉口水位(Zh)、螺山流量(Q)、螺山汉口河段54年以后的累积淤积量(V)之间的经验关系。即水利论文(r/Oh|8y!Z5Mn

Zl=f(Zh,Q,V)

7B"i wA8T otQ$l0

(6)

u3{1Y {)]{wF0

选取1955~1992年38年汛期50000ms/s以上流量系列率定公式结构,由水流运动方程可知,上游水位与下游控制水位之间是线性关系,与流量Q和淤积量V之间是复杂的非线性关系。将公式(6)Tailor展开,可写成Zh、QiVj(i=0,1,2…,j=0,1,2…)的线性组合形式,由于Qi(i≥3)或Vj(j≥2)的各项其偏相关系数普遍小于0.5,因此,在率定公式时取如下公式结构

]mP?q7o{-V0

vlow%[/G)p0

(7)水利论文)e9x a[(QKU \m ?F

式中 ai(i=1,2…6)为常数,分别为a1=5.96165,a2=0.70677,a3=2.278×10-4,a4=-0.147,a5=-1.898×10-9和a6=3.059×10-6,公式中各项偏相关系数见表1。

3n rSnp;h0

表1 式(7)中各项的偏相关数水利论文*`q&Ubz,A a/~6x
Table 1 Correlative value of different parts in the equation (7)

y9d;E"Sd0

项目名称Zh QV Q2 QV

偏相关系数0.99930.99860.9963

0.994水利论文%z-x }e:{2f4| |g

0.9973

  分别用RFS-ANN模型和公式(7)计算螺山1955~1992各年流量大于50000m3/s的系列,计算值与实测值比较见图3。两模型平均误差都控制为0.02m,RFS-ANN模型螺山水位计算值与实际值之间的最大误差为0.15m;式(7)计算结果与实测值最大误差为0.20m。在同样控制误差的条件下,回归拟合模型由于经验公式各个项的选择经验性半经验性较强,其选择具有一定的任意性,模型拟合精度比神经网络模型稍差。本文两模型在拟合精度上尽管存在一定的差异,但两模型对历史资料的模拟精度满足计算的要求,模型能够较好地反映螺山水位变化的规律,可以用于水位变化规律的分析和预测。

&t}Aaq*YUd0

水利论文1o.IyC/j0j+z z0jFz,I
图3 计算洪水位与实测值比较水利论文b#v P!nu[J)H {S
Fig.3 Comparison of the computed
水利论文 ]-PZ!Q grt+q.e4NS
and measured flood level

a;x*MVq0c'a:j(X0
       

3.3 资料完整条件下的模型预测比较 水利论文c-T7W0v?HN~R&l_

  图4为RFS-ANN模型和回归模型对93年流量在50000m3/s以上的水位进行预测的结果。由图可知RFS-ANN模型预测结果与实测值之间最大误差小于25cm,平均误差为3cm。回归模型预测结果与实测值之间最大误差小于28cm, 平均误差为3.2 cm。回归拟合模型对93年水位预测精度比神经网络模型预测精度稍差。水利论文)Z.b*Bw3x!Mj

水利论文Bl Z$h3aU&gQ5UV(v
图4 93年螺山水位预测值与实测值比较(考虑淤积)水利论文u1|wM7[m l&d]jj-n
Fig.4 Comparison of computed and measured flood level at Luoshan in flood period 1993

^8}UV&MI1sz0
     


1ZR'{3S4x|5[0图5 93年螺山水位预测值与实测值比较(不考虑淤积)水利论文u6P6phMy TUh~
Fig.5 Comparison of computed and measured flood level at Luoshan in flood period 1993
水利论文F1^P!CU$d?}

  图5为RFS-ANN模型和回归拟合模型在不考虑河床变形条件下对93年水位的预测结果与实际水位比较。由图可知:不考虑河床累积淤积情况下两模型预测结果均比实测值小,分析出现这种情况的根源是如果不考虑河床变化对水位的影响,其实质是在汉口水位流量关系保持不变时将螺山水位表达成流量的单一函数,在下游河床抬高后,同一流量下的水位仍为河床变化之前的水位,忽略了螺山汉口河段累积淤积对水位壅高的影响,因此会出现计算结果明显比实测值小的情况。随着淤积量的增加,水位流量关系不是单一曲线,水位随着淤积量的增加而抬升。因此,尽管两模型计算原理不相同,其对水沙运动规律的反映基本一致,反映了泥沙淤积对水位影响的客观自然规律。

5U-W)y3jdJ]0

  近年来汉口水位比较稳定,螺山站出现小流量、高水位的主要原因在于螺山汉口段的不断累积淤积所造成的河床形态的改变,在建立螺山水位变化预测模型时应加以考虑。

0Y~[p{ {%`f7p0

图6、图7分别为RFS-ANN模型对不同流量级下(分别为60000m3/s、65000m3/s、70000m3/s)螺山水位随河段累积淤积的变化模拟结果。水利论文?8sH?4P"{7WRC

水利论文.V:r/^XQ;weA
图6 河床淤积与水位的关系(RFS-ANN模型)
~y~ k0ko)Xc:}(E"h0Fig.6 Relationship between deposition and flood level
水利论文-C3W7GkD/~z7a

水利论文4} ma1t$qm#xw
图7 水位随累积淤积量变化(回归模型)
G$U.vI.|5uc8G0Fig.7 Variation of water level with the
水利论文|l)bp1j'DI6A
accumulation of deposition

(u3~(L[#ZB7a I;p"L0

  由计算结果可知,流量在60000m3/s和65000m3/s时,河道每淤积100×10t泥沙,RFS-ANN模型模拟的螺山水位分别增加6.2cm和6.3cm;传统回归模型水位抬高分别为3.63cm和5.2cm。回归模型模拟结果比神经网络模拟结果略小。流量在70000m3/s时,水位受河床淤积影响明显,河道每淤积100×10t泥沙,水位抬高分别约7cm(RFS-ANN模型)和6.7cm(回归拟合模型),两模型在水位抬高数量上基本相近,即螺山汉口段每淤积10亿t泥沙将引起螺山断面水位抬升近70cm。正是这种大量泥沙的不断淤积导致了近年来水位上升、水情的不断恶化。如果不对螺山汉口段的泥沙采取适当的措施,这样的小水大灾会进一步恶化,甚至威胁武汉的安全。

-U/R{ZeNe0

  螺山在60000~70000m3/s之间出现流量增加水位随淤积量变化增加的原因在于泥沙不同淤积部位(分河槽淤积和滩地淤积)对不同流量下的水位影响存在一定的差异。根据资料统计[8],1954~1993年间螺山至宜昌河段以及洞庭湖区年平均泥沙淤积180×103t。水利论文 ?%gI {5D

根据长江科学院的研究成果[10],用地形法计算1959~1993年得到的50000m3/s流量(相当于平滩水位下的流量)下莲花塘至汉口淤积量为0.624×10m3,泥沙干密度以1.4t/m3计算,折合为0.87×10t;按输沙率法计算,螺山至汉口段在1959~1993年间共淤积量为1.96×10t。可见,螺山至汉口段淤积滩地淤积与主槽淤积量基本相同。因此,汛期出现随着流量的增加,水位的抬高值增加是主槽淤积和滩地淤积共同作用的结果。水利论文;`5p:B\kDJ

3.4 资料不完整条件下的模型预测比较水利论文IQVqU,l

  1998年长江中游发生了大洪水,经济损失高达1000多亿元。为了减轻洪水灾害,许多水利学者对1998洪水成因提出了见解。与1954年长江特大洪水相比,1998洪水主,水位很高等特点,螺山站的水位抬高了近1.8m。研究和探讨螺山水位变化与螺山至汉口河段累积淤积量之间的关系对荆江及洞庭湖的治理、减轻长江中游洪水灾害具有重要的研究和指导意义。水利论文X5b;M"C`\%]+C Q

  由于资料原因,1994年以后的螺山至汉口段泥沙淤积量只有通过以前该河段以及长江干流泥沙运动变化特征进行估算。在长江干流宜昌及四水来水来沙变化不大、下游汉口站输出水沙量和水位基本保持不变的条件下,汉口以上长江干流及洞庭湖区的总淤积量也不应该有较大的变化。在洞庭湖区淤积量逐年减少的同时,必然导致螺山至汉口河段淤积量的相应增加,即泥沙淤积发生“搬家”现象。

Ke,YV;qf4t0

  通过整理长江中游和洞庭湖水沙资料可得螺山至汉口河段和洞庭湖区泥沙淤积比λV与三口分沙比λS之间的定量关系[8]

Q w[4O_N~/V0

λV=(1/2.14λS+0.123)-1

*Yh`Yn1{E BM5a0

(8)

J%j,m7|,n/cj+J0

  由式(8)计算结果可知各个时期湖区和螺山至汉口河段多年平均泥沙淤积总量基本保持在180×106t。分沙比的减少只是调整了螺山至汉口河段与湖区之间的淤积关系,使得本应在湖区淤积的泥沙转移到螺山至汉口河段,增加了螺山至汉口河段的淤积量。螺山至汉口河段淤积量的增加,抬高了洞庭湖区和荆江河段下游的侵蚀基准面,导致该地区的防洪形势恶化,这也是近年来频繁出现小流量高水位的直接原因。

Tu4aWY _3gfn0

自然和人为双重因素的作用,三口分流分沙比在1954~1982年以前呈明显的减少趋势,1982年以后仍有减少趋势但不明显。为了研究问题的方便,本文将1988年以后的分沙比取为1983~1988年间的平均分沙比,即近似地假定1983年以后的分沙比基本保持在0.197,以此计算1983~1993年间螺山至汉口河段年平均淤积量为80.67×106t,与实测值80.45×106t相近。如果1993年以后的宜昌和四水来沙条件不变,洞庭湖区淤积率保持在74%,则螺山至汉口河段1993年后平均年淤积80.67×106t,1994~1998年共淤积403×106t。

tZ.A^"l$QJ0u0

水利论文+U!Ra-a*a3Sk}Y.T~)Y
图8 1998螺山洪水位预测(RFS-ANN模型)水利论文Q&TF&h~)O`
Fig.8 Observed versus predicted water level at Luoshanby Model 2 for
flood period 1998'水利论文0f8E%N dA

水利论文z}OTry
图9 泥沙淤积对1998洪水位的影响水利论文k(`1Pk dAT
Fig.9 Influence of sediment deposition in 1998' flood level

8\2IUzENF3b0

  为了更明显地反映河段淤积对螺山水位的影响,本文用RFS-ANN模型在假定1954年地形和考虑1954年后河道累积淤积两种情况下对1998年螺山水位进行模拟。计算结果与1998年实测值比较如图8所示,由图可知:考虑1954年后的淤积模型计算结果与实测值之间平均误差为0.008m,说明模型在对水位进行预测时必须充分考虑泥沙淤积对水位的影响。同时计算结果还表明,如果螺山至汉口河段1954年以后不发生累积淤积,1998年洪水水位最大减低1.8 m,防洪压力会大大减轻,受灾范围也相应会得到控制。

!\b)uovx#O0

  图9为回归拟合模型在考虑1954年以后河道淤积和不考虑1954年后河道淤积的两种情况下对1998年汛期水情模拟的结果与实测值比较。考虑河道淤积计算结果基本与实测值吻合,一方面说明螺山至汉口河段在近5年淤积量在403×106t左右,另一方面模型对洪水预测的实用性和准确性。如果不考虑1954年后的地形变化,即以1954年河段地形,1998年来水条件计算1998年洪水过程,水位明显下降,汛期平均降低0.92m,最大降低近1.4m。

$]6wua2o7A,D0

  两模型都表明,河床的不断累积淤积抬高了螺山站的洪水位,加重了洞庭湖区和长江中游的洪水灾害,增加了长江中游湖北、湖南两省的防洪困难。因此,目前从防洪利益出发,对螺山汉口河段进行疏浚和遏止该河段的累积淤积将有利于该段河道行洪,有利于减少和减轻长江中游的洪水灾害。

^!J],L,j'BBH5kY0

  从两模型对1993年和1998年螺山水位模拟预报情况分析,神经网络模型与回归拟合模型对1993年的预测差别不大,但1998年预测差别比较明显;另外,两模型在河道淤积对水位抬高影响程度的体现上也存在一定的差别。神经网络模型只要选择好影响因子,选择好网络结构及误差控制,其模拟精度就由网络自身确定,而经验回归方法在恰当选择好影响因子后还必须选择好公式的结构,如果公式结构选择不当,将会在很大程度上影响预测精度。因此,如果对水沙规律了解比较全面,对水沙相互之间的作用了解得比较清楚,建立直观、明了的回归拟合模型不会影响模型的预报效果,反之如果对河段水沙之间的关系了解不是非常清楚,利用神经网络的非线性映射功能可以较好地达到预测的目的。另外尽管两模型在泥沙淤积与水位抬高之间的相互影响数量上有一点差别,但两模型模拟和预报结果都反映在建立模型时河道的淤积对水位的抬高必须考虑,否则模型精度和预报性能将得不到保证。水利论文|+CVr(vi:BgJo s

  本文通过模型对螺山水位的模拟,还得到这样的一个认识:长江中游和洞庭湖区洪水灾害的核心是泥沙淤积和泥沙淤积部位的问题,解决长江中游和洞庭湖水患灾害的关键必须处理好泥沙来量和泥沙淤积分布。从治沙必须控制泥沙来源的角度,上游通过水土保持、建拦沙坝等减少宜昌上游泥沙来量应该成为我国的一个长期策略,是解决长江中游泥沙问题的根本办法。但水土保持等措施需要一段较长的时间,目前可以通过在长江中游低洼区域沉积一部分泥沙,尽量减少螺山至汉口段的淤积来有效降低洞庭湖区和荆江河段的洪水位。

-L&N&o4Nq1s"G'_,K^0

4 结论

;^+WZ qjhc*L0

  RFS-ANN模型与传统回归拟合模型相比,RFS-ANN模型在模型结构确定、模型容错性能等方面优于回归模型,在模型学习速度、算法复杂程度、实用性和直观性等方面逊于回归模型。当两模型在充分考虑影响水沙变化因子的前提下,都具有较高的模拟和预测精度。本文将两模型用于螺山站水位的模拟和预测,结果与实测值吻合理想。同时,两模型对螺山至汉口河段淤积与螺山水位抬高关系进行了计算和分析,结果表明,下游淤积是螺山和湖区水位抬高的根本原因,是近年来长江中游小水大灾的主要影响因素,在建立该河段水沙模型时必须考虑河道泥沙淤积对水沙运动的影响。水利论文1y:z:z yDkw

参考文献

-b%? fr'S0

[1]吴道喜。复杂河网洪水调度数学模型研究。武汉水利电力大学硕士论文,1997.水利论文Ta6`g3ts

[2]N. Karunanithi, et al.(1994). Neural Network for River Flow Prediction. J. of Computer Civil Engneering. 8(2): 201~220.

(dL8hn DKH:qz0

[3]Knanina, E. D. (1990). A Simple Procedure for Pruning Back-Propagation Trained Neural Networks. IEEE Trans. On Neural Networks.1(2):239~242.

7Y]4HJa6dN#D%M0t0

[4]李义天,倪晋仁。泥沙输移对长江中游水位抬升的影响。应用基础与工程科学学报。Vol.6, 1998.

LP{v+y8M| F7L{J0

[5]李义天,李荣,邓金运。长江中游泥沙输移规律及对防洪影响研究。泥沙研究,2000,(3)。

F _+RY0l/r[.p0

[6]蒋学东。洞庭湖口及长江干流城螺河段水情变化主要研究成果综述。水利规划研究,1994,(1).

s9j2soR3S4g0G0

[7]段文忠。下荆江裁弯与城陵矶水位抬高的关系。泥沙研究,1993,(1).

*R)y$xV}:O:a,E.P0

[8]施修端等。城陵矶水位流量关系及其水沙变化分析。长江水利委员会水文局科研报告,1997年12月。

y!|1` t4T ~SF@(u0

[9]Li Rong and Li Yitian. Applied Researches of Neural Network in Flood Prediction”。International. Journal. of Sediment Research. No.2,1999.

D$h7K@)[0

[10]Li Yitian, Gan Minghui and Li Rong.“Sediment deposition and flood disaster in Dongting Lake” Proceedings of the Seventh International Symposium on River Sedimentation.1999.水利论文D}6x a s(^9}

[11]欧阳履泰。初论荆江和洞庭湖及其治理。湖南水利,1981,(3).水利论文cP)uDS)I+yQP"{!~ p

[12]焦李成。神经网络系统理论,西安电子科技大学出版社,1990.

P.SR!U1H(Ejqsd2l;e0

[13]胡铁松。神经网络预测与优化。大连海事大学出版社,1997.

!T{ h?p!t!E5F/Y0

[14]李义天,李荣,邓金运,孙昭华。泥沙输移变化对荆江及洞庭湖区洪水位抬高的影响。武汉水利电力大学科研报告。1999.

gh [m0pI2Mr0
TAG: 黄伟 李荣 李义天 模型 神经网络
顶:21 踩:17
【已经有99人表态】
16票
极差
9票
很差
7票
较差
12票
稍差
11票
稍好
15票
较好
11票
很好
18票
极好
下一篇:中国河口最大浑浊带刍议(沈焕庭,贺松林,茅志昌,李九发)
上一篇:1977~1996年黄河下游水文断面反映的河床演变(梁志勇)
查看全部回复【已有0位网友发表了看法】

广告投放

广告投放