关于水库洪水预报方法的探究.pdfVIP

关于水库洪水预报方法的探究

摘要：莱西市境内有大型水库一座，中型水库两座，汛期能否进行正确的洪

水调度，关系到我市的工农业生产和人民群众的生命财产安全，及时准确的洪水

预报，是进行洪水调度的前提和基础。分析降雨径流，提出水库洪水预报，要探

讨两大问题，一是产流，二是汇流。产流是解决由毛雨过程求净雨过程的问题，

其中确定前期影响雨量Pa为主要前提；汇流是解决由净雨过程求入库流量过程

的问题，精度较高的计算方法是经验单位线法。本文结合我市洪水预报的实际情

况，对前期影响雨量Pa、净雨R、经验单位线法进行分析探究。

1前期影响雨量Pa

前期影响雨量Pa是反映雨前土壤蓄水量变化的参数，即雨前土壤蓄水量。

其主要影响因素是表层土壤的组成成分和厚度、降雨量的大小、蒸发量的多少、

土壤入渗量的快慢等，较适用的计算方法大致分两种：

1.1扣损法

无雨日：Pa(t+1)=Pa(t)-E；有雨不产流日：Pa(t+1)=Pa(t)+P-E；

有雨产流日：Pa(t+1)=Pa(t)+P-R-E≤Im；Pa(t)，Pa(t+1)—第t天，第(t+1)天的

Pa值。

P—流域平均降雨；R—净雨量；E—第t日的陆面蒸发量，Im—土壤最大持

水量。

我市大型水库水文站有E601蒸发皿观测资料，就是陆面蒸发量。在汛期，

蒸发量变化范围1—7毫米，高温天气5—7毫米，晴天3—5毫米，阴天2—3

毫米，雨天1—2毫米，据统计，6月份平均日蒸发量为4.7毫米，7月份为4.5

毫米,8月份为3.9毫米，9月份为3.6毫米，逐月呈递减趋势。

有雨产流日公式中P—R—E指土壤入渗的部分降雨量。严格的讲P—R—E

的部分降雨量不能全部入渗土壤，总有一部分降雨量被坑、洼、库、塘拦蓄，但

计入拦蓄部分，使Pa的计算更接近实际，当土壤达到饱和，即土壤含水量达到

Im，可以认为后期降雨量不再入渗土壤，全为产流净雨，这时的Pa就达到最大

值Im。

扣损法看起来较麻烦，实际运用还是可行的，关键在于确定土壤的最大持水

量Im，Im的计算可选用流域久旱之后发生的一场大暴雨，用公式Im=P-R-E计

算，久旱指Pa=0，大暴雨是指土壤含水量能够达到饱和，R取实际的径流深，

选择实际发生的若干次暴雨资料，即可确定Im的值。Im与流域的土壤状况和径

流坡度密切相关。经分析，我市产芝水库土壤植被最好，高格庄水库次之，北墅

水库最差，径流坡度分别为0.0014；0.0027；0.0032。经测算，产值水库的Im为

100—105毫米，高格庄水库Im为93—98毫米，北墅水库Im为85—90毫米。

Im值在实际取用上不能采用同一固定值，这是因为土壤干燥时，坑、洼、库、

塘要截留一部分雨量，降雨较缓时土壤有充分时间吸收更多的雨量使之深层也达

到饱和，故针对一次降雨过程，当Pa值较小，降雨量缓慢Im取值相对稍大，

反之较小，这样可以提高预报精度。

实际上，流域的植被情况（主要是农作物、采木树林的影响）和流域的蓄水

情况（坑、洼、库、塘的影响）都会给洪水预报带来偏差，这些情况有没有固定

的公式可以计算，而扣损法计算过程中Im列入了这些因素的影响，使理论与实

际更为相符。

1.2系数法

Pa(t+1)=K【Pa(t)+p】但须控制Pa(t+1)≤Im

式中，K—流域土壤蓄水的日消退系数。我市三座水库预报均采用此法，精

确与否取决于选用的系数K，K值取用范围大致在0.85—1.0，当降雨小于5毫

米时K=0.9，当降雨在5—15毫米时K=0.95，当降雨大于15毫米时K=1.。

此法是近三十几年来采用较多的计算方法，简单适用性好。系数K=1-E/Im，

按每日的蒸发量逐日推算，将提高预报精度。多数水库没有蒸发观测，可近似的

按每月平均蒸发能力E推算，系数K取值可以凭借一定的工作经验，相对灵活

取用。

2净雨量R

净雨量是指毛降雨扣除损失后的产流雨量，损失量包括蒸发、植被截留、填

洼、土壤入渗等，其中土壤入渗对净雨量R影响较大，R的计算方法大致分为三

种：

2.1径流系数法：R=ap,a—径流系数法，选择水库流域若干次降雨径流资料