流域水文模型的研究大约始于上世纪50年代,70年代至80年代中期是其蓬勃发展时期,一些比较著名的模型,如美国的斯坦福
1.1预报效果不理想
现有流域水文模型都是概念性模型,水文现象十分复杂,使得人们至今还不能用数学物理方程严格地描述其中每一个子过程。因此,现有流域水文模型在许多结构环节上仍主要借助于概念性元素模拟或经验函数关系描述。这样的模拟往往只涉及现象的表面而不涉及现象的本质或物理机制,因此,现有流域水文模型包含的许多参数都缺乏明确的物理意义,只能主要依据实测降雨和径流资料来反求,而这样求得的模型参数必然带有经验统计性,只能反映有关影响因素对流域径流形成过程的平均作用。这是现有流域水文模型拟合一组资料中的大部分虽可达到令人满意的程度,但对该组资料中的个别特殊情况,或者该组资料以外的另一些资料却不一定能获得令人满意的拟合结果的症结所在。
1.2用最优化方法确定模型参数对实测降雨径流资料的依赖性很大
现有流域水文模型一般都包含有两个以上的参数需要由实测降雨和径流资料来反求。在数学上处理这类
流域水文模型的输入是流域上各点的降雨过程,输出是流域出口断面的流量或水位过程,因此,它是一种输入具有分散性和输出具有集中性的模型。现有流域水文模型在结构上一般与此并不匹配,在实际应用中考虑这一问题时,几乎无一例外地采用将全流域按雨量站划分成若干个单元面积的方法,认为当单元面积的尺度小到一定程度时,即可作集中输入和集中输出的流域水文模型来模拟该单元面积的径流形成,最后将各单元面积对全流域出口断面输出的贡献迭加起来作为其出口断面的输出。显然,这种实用的处理方法仍是不完善的。
1.4有效的计算方法问题
与计算机技术发展不可分割的流域水文模型,总是以离散方式来模拟流域径流形成过程,因此采用有效的计算方法是十分重要的。流域产、汇流过程是一个连续过程,当用数学模型模拟时,必然因离散带来离散误差。计算机的数字有效位数总量有限,这又必然带来误差。如果离散化方法不适当,计算方法不够好,则非物理的数值效应可能会导致背离物理图景的计算结果。这就表明,流域水文模型是否成功,除了合理的结构和参数外,还有赖于采用合适的离散化格式和有效的计算方法,换言之,流域水文模拟存在数值计算是否满足重系数:
Wj=〔WW∧
通过加权后可得到总输出Q
Qj=WX+§=WX+
2.2黑箱模型预报原理
在实际应用中,影响一场洪水大小的因素,可归纳为降雨量、当前流域断面出口流量、蒸散发能力三个方面。在黑箱模型中,可将三个影响因素作为模型的输入,计算流域出口的产流量Q。
现假设有某流域,需要预报该流域t时刻后v个时段(一般采用一小时作为一个时段)的出流过程,且流域面用泰森多边形划分为f个单元面积,同时具有t时刻的降雨量m、蒸发量e
E=
X=
假设权重系数:
其中θ为计算阀值
其中每一行系数,表示一个预报时段的权重系数。
则流域出口断面t时刻出流过程为:Qv==WvX
2.2 权重系数的求解方法
权重系数的求解过程,实际上是记录洪水涨落规律的过程,这些规律的表现形式就是权重系数。对于一个流域的权重系数求解,必须使用该流域的多次洪水进行求解,一般预报时段数就是需要使用的洪水场数。
为清晰分析权重系数的求解过程,现假设流域有n场洪水资料,其中包括降雨M、蒸发E、流域断面总出流量Q。且已知输入向量
X=MER=mmKmee∧er
Qv=
Wv=其中θ为计算阀值
根据黑箱数学模型
Qv=X Wv 可得到j个线性方程组
其中
j=1,2,3, ┄
k=经验常数(一般k=3);
f=流域单元数;
v=预报时段数;
n=洪水资料场次数
4 结语
洪水预报的黑箱模型还不成熟,但对大量无实测资料的中小水库进行洪水预报有作较大的意义,可以一劳永逸的解决洪水预报模型问题,为标准化的编制预报软件打下理论基础,是很有价值的一种思路。
参考文献
赵人俊,流域水文模拟,水利电力出版社,1984
苑希民 李鸿雁 刘树坤 崔广涛,神经网络和遗传算法在水科学领域的应用,2002
朱华,水情自动测报系统,1992
