第二章 降雨径流相关预报

第二章降雨径流相关预报在现代水文预报中，虽然大量使用流域水文模型，例如新安江模型、萨克门托模型、水箱模型和陕北模型等进行流域降雨径流预报。但是，不少生产单位，尤其是一些大型水库的管理单位，他们在长期的工作实践中已建立了一套适合于当地实际情况的经验性降雨径流预报方案。2.1降雨径流相关图的形式降雨径流经验关系曲线有各种形式，一般有产流量R=f（次雨量P，前期影响雨量p,季节，温度）、R=f（前期影响雨量。，洪水起涨流量。°）和考虑雨强的超渗式关系曲线形式。这里介绍国内普遍使用的产流量与降雨量和前期影响雨量三者的关系，即P〜P〜R相关图。使用P使用P〜P〜R洪水初期的P值，把时段雨量序列变成累积雨量序列，用累积雨量查出累积净a雨，由累积净雨再转化成时段净雨量序列；另一种方法是根据时段降雨序列资料直接推求时段净雨序列。第一种方法的缺点是在整个洪水过程中，使用一条P〜R曲线，没有考虑洪水期中P的变化。而后者的不足是，当时段取的过小时，a一般时段雨量不大，推求净雨时的查线计算易集中在曲线的下段。两种方法的结果存在差别，至于何者更接近实际也很难断言。2.2前期影响雨量P的计算aPA由前期雨量计算，也称前期影响雨量，是反映土壤湿度的参数。其计算公式为若前一个时段有降雨量，即P1>0时，则Pa广K(PaI+P/ 3若前一个时段无降雨时，即P1=0，则P广牛" (2-2)式中：K为土壤含水量衰减系数，对于日模型而言，一般地取K^0.85；P和a,t—1Pat分别为前一个时段和本时段的前期影响雨量；P—1为前一个时段降雨量。式(2-1)和(2-2)为连续计算式。由于|Pa,t=K(Pa,,—i+P—i)Pa,t—i=K(Pa,t—2+Pt—2)P,t—2"(Pa,t—3+P—3) (2-3) lPa,t—n+1=K(Pa,t—n+P—D将式(2-3)各行逐一代入得到P=KP+K2P+K3P+…+Kn(P +P) (2-4)a,t t—1 t—2 t—3 a,t—n t—n式(2-4)为向前倒数n天的一次计算式。一般取15天既可满足计算要求。用I表示流域最大损失量，在数值上等于流域蓄水容量。以mm表示，通m常Iw60-100mm。当计算的P>I时，则以I作P值计算，即认为，此后的降雨量P不再补充初损量，全部形成径流R。当计算时段长Atu24h时，土壤含水量衰减系数K应该用下式换算K=KD1/n (2-5)式中：N=24/At，KD为土壤含水量日衰减系数，K为计算时段是At小时的土壤含水量衰减系数。

根据计算出的流域平均降雨量P和P所产生的径流量R，以及相应的前期影响雨量P，便可建立降雨径流相关图。显然，P是影响降雨径流关系最主要的因素，因为流域的产流决定于非饱和带的物理特性，而前期影响雨量的物理含意是土壤含水量，它反映了非饱和带土壤的物理性质，但它不是唯一的因素，在有些情况下，其它因素不可忽略。除了前期影响雨量以外，季节、降雨历时、流域平均雨强等也不同程度地影响着降雨径流关系。由式R=f(P,P)建立起来的三变数降雨径流相关图，如图2-2所示。由于结构简单，使用方便，且能满足精度上的要求，所以被广泛地应用于雨洪径流预报。图2-2 诅河猴子岩流域降雨径流相关图由式R=f(P,P.)建立起来的三变数降雨径流相关图具有下列共同特点：当P一定时，P越大，R也就越大，所以P等值线呈左小右大。P=0线的延线交P的截距为I,P丰0线的延线交P的截距为D(流域土壤缺水量).在P和R取同一比例时，P=七线与横坐标的夹角略大于45度线。由于超渗产流和局部蓄满产流，也就是说，末满足流域平均土壤缺水量就产流，因此曲线下端曲率较大，上端由于土壤渐趋饱和而逐渐趋于直线且与P=I平行。在同一流域平均径流深R下，P越小产流面积就越小，所需的雨量就越a大，因此曲线下端的曲率随着P的减小而增大。a在同一P情况下，P越大，径流系数a越大。a在点绘三变数降雨径流相关图时，应考虑上述特点来定线。2.4相关图推流计算用P〜P〜R相关图作流域降雨径流计算的步骤是：首先摘录每条P〜R曲线的各点坐标，把P〜P〜R曲线坐标和土壤含水量衰减系数K、土壤最大损失量七、计算时段长&以及计算开始时的前期影响雨量P，作为模型参数输给机瞬程序中的相应变量，然后根据时段降雨序列计算出每个时段的净雨量。2.5计算实例已知黄河流域伊河上陆浑全龙门镇区间流域面积为1826km2，现得到其P〜P〜R曲线的坐标点。已知降雨过程为12.327.544.740.034.013.29.46