第三章单位线分析计算.doc

第三章单位线分析计算本章着重介绍如何由净雨量过程预报流域出口的流量过程净雨量经过流域汇流形成出口的流量过程线，流域汇流历时是降雨径流预报预见期的来源，流域汇流物理过程是编制预报方案的理论依据舍尔曼时段单位线基本原理舍尔曼(L.K.Sherman)于1932年提出了单位线的概念其定义为：流域上分布均匀的1个单位净雨直接径流产流量，所形成的直接径流过程线，即为单位线，记为UH1个单位净雨是指单位时段内单位净雨深单位时段长可以任取，例如2h、3h、6h,等而单位净雨深通常取为10mm而实际发生的净雨，常常不是1个时段，也不是1个单位，应用于分析单位线时,有一些假定这些假定可归纳为以下两点：(1) 如果单位时段内净雨深不是一个单位，而是n个，它所形成的出流过程线，总历时与UH相同，流量则是UH的n倍2) 如果净雨历时不是一个时段，而是m个，则各个时段净雨所形成的出流过程之间互不干扰，出流过程的流量过程等于m个流量过程之和由以上假定，净雨rd、出流Qd与UH的纵坐标q之间的关系如下：mQd,trd,iqti1(2-1)i1式中：i1,2,3,m，为净雨时段数Qd和q的单位为m3/s，rd则用单位净雨深的n倍来表示。

如果UH已知，根据上式，可由净雨转换为出流，计算十分简便关键是如何求得UH可以根据流域的实测水文资料，分析出净雨及直接径流过程后，依据上式推求出UH，为式(2-1)的逆过程单位线的推求推求UH是使用次洪时段净雨深及相隔为计算时段长的直接径流时序过程前者由次洪降雨量经过扣损后得到，后者由径流过程线分割地下水后得到这里需要补充说明一下具体问题：(1) 由扣损方案求得的次洪净雨深，常不等于过程线分割得到的实测值，为了不把扣损的误差带入汇流计算，需要将计算值改正，或谓平差在平差前应分析误差的来源，作出较为合理的修正，以不改变原来的雨型为原则2) 计算时段长t的确定，主要考虑峰形的控制，一般取涨洪历时的1/3〜1/4，常与流域大小成正比对于多时段净雨而言，传统的单位线分析万法有以下一些：(1)分析法m由Qd,trd,iqti1i1，得到Qd,1刼1(3-2)Qd,2“2"1(3-3)Qd,3“3rd,2q2rd,3q1(3-4)因此，上式组成的方程组为一个多元线性代数方程组，求解该代数方程组可以得到q1、q2的数值最简单的解法是逐一消去法由上面的方程组可以看出，由式(3-2)，已知Qd,1和S，解出q1:rd,1(3-5)q2rd,1q1已知，将其带入式(3-3)，可得到(3-6)如此递推下去，得到,t,trd,iqti1(3-7)rd,1计算结果见表3-1和图3-1所示。

直接代数法对于降雨与径流实测资料没有误差，流域汇流符合线性时不变43线位单43线位单修正前单位线过程修正后单位线过程系统时，能得出正确的唯一解但实际情况是：提供分析UH的实测资料存在着观测与分析误差，流域汇流为非线性系统，因此直接分析解法不但不能得出唯一解，由于误差累积，其解常很不合理如分析得到的q~t呈锯齿形，如表3-1和图3-1的退水段，有的甚至出现负值，无法继续计算下去所以很少有人使用图3-1分析法得到的单位线表3-1单位线分析法结果日时q时序Qdrdq计算值q修正值Qd计算值760000012118676761861826672092106698031935616617193864245048949024531251900370355186418612802162401309907850168155867685608910557212940089733921810277395227610011202503819961214216221311213801912741814401024110150000(2) 试错法这个方法是假设单位线，目估对比推流的与实测的Qd~to当两者最接近时，所假设的UH即为所求。

初始的UH，可用其它洪水已分析得来的成果，或用斜线分割法的结果，也可任意假定试错法应用比较广泛，有的用科林法试错，有的凭经验试错但主要确定是单位线过程的初始假定比较困难，试错过程有时也会出现不合理现象，对其过程进行修正也不易做得完好3) 系统识别方法UH既然是线性系统的单位响应，进一步可以应用线性系统鉴别的方法推求UH的最优解，如有约束的最小二乘法和回归法等4) 各种方法的特点如上所述，推求单位线的方法主要有分析法、试错法和系统识别法试错法又可分为目估试错法和科伦试错法，后者是一种迭代性质的试错法分析法当降雨为1个时段时能取得很好的效果，但当降雨时段数大于1个时段，分析法推求的结果可能出现锯齿型，甚至出现负值，主要原因是误差的累积目估试错法是通过人工目估，逐步试错，使计算的流量过程与实测的流量过程吻合，避免了误差的累积，但任意性大，吻合程度缺乏客观标准，并且要求调试者具有相当的经验，费时费力科伦试错法是首先拟定初始单位线，然后按初始单位线计算出除最大降雨时段外的其它时段降雨产生的流量过程，并从实测流量过程中将其减去，差值即当成最大降雨时段产生的流量过程，据此推求新的单位线，将新单位线作为初始单位线，重复上述步骤，直至新旧单位线之间的差值达到给定的误差要求。

优点是当初始单位线拟定后，不必目估试错，能自动逐步修正，适合计算机编程缺点是计算结果及迭代的收敛性与初始单位线的拟定有关，推求的结果不一定是最佳结果另外当降雨时段较多时，迭代难以继续下去因此，科伦试错法仅当降雨较集中在最大时段或降雨时段数较小时，效果才较好单位线的综合为求得流域的单位线预报方案，需分析多次洪水多数流域，各次洪水所得到的UH会有差别这时，首先要检查原始资料的观测误差，及根据原始资料计算用于分析UH的rd~t和Qd~t过程中的操作误差，作出可能的改正，然后对UH的变化根据UH原理上存在的问题进行分析及处理如果各次洪水的UH变化不大，可求出平均UH平均UH的绘制须注意，平均纵标值不宜取同时流量的平均值应先根据各个UH，计算出平均峰值及峰现时间，确定峰点位置然后根据各个UH的外形，初绘光滑的平均UH，最后修正并校核总量为单位雨深单位线时段转换UH是有一定长度的，同一次洪水，如单位时段长度不同，实际雨强不等，UH不相同如原来UH时段长为T,现推求2T的UH,根据线性假定，可将UH滞后Th,与原UH相加，将该过程线的纵坐标除以2即得2T的UH该UH的峰值应等于上述两个UH的交点，比原来UH峰值低并滞后，这是由于历时增长，雨强降低所致。

因此，对一个给定的流域综合UH方案，或将不同地区的UH方案综合，需将不同时段的UH换算为同一时段UH是作为一个经验方法提出来的，但方法简易，使用精度不错，生产上直沿用至今系统概念引入后，建立了UH的理论，明确了方法的物理实质，这对于推求UH及正确的应用都有益UH是个黑箱模型，并和流域上实际的水力状态没有关系，使用时必须有实测水文资料，这是它的应用受到限制的主要原因纳须单位线(Nash)纳须模型是现行概念性流域汇流模型中得到比较广泛应用的一种，该模型结构简单，包含的参数较少，且具有较好的适应性本文介绍纳须模型的结构及其参数率定方法纳须瞬时单位线1945年克拉克(C.O.Clark)首先提出瞬时单位线的概念所谓瞬时单位线，就是流域上均匀分布的、历时趋于无限小、强度趋于无穷大、但净雨总量为1个单位的净雨所形成的流域出口断面过程线通常用u(t)表示1957年，纳须(Nash)把流域看作是一连串的n个相同的“线性水库”，如图3-2所示可以推导出一(3-8)u(t)个单位的瞬时入流进入水库系统后，其对应的出流即瞬时单位线的数学方程式为式中：为伽马函数，当n为整数时，其值为(n1)!；n为相当于水库个数；kk(n)k为一个线性水库的蓄泄系数。

图3-2纳须模型示意图纳须时段单位线在实际应用中，需要将瞬时单位线转换成时段单位线，一般用S(t)曲线按照S(t)曲线的疋义，S(t)等于瞬时单位线的积分，即S(t)u(t)dt，对该式进0行积分，可以很方便地导出S(t)曲线，再利用线性叠加原理，u(t,t)S(t)S(tt)，从而得到时段单位线即为2!n1tt.nu(tit}ek[ekiiRhT)nini](3-9)式(3-9)为当n为自然数时，计算步长为t时的纳须时段单位线的计算公式纳须单位线参数n、k的确定如何确定参数n和k是纳须模型应用中的重要问题之一纵观目前所有的确定该模型的方法，按依据的资料分，则可归纳为3类：一是依据降雨和径流一一对应资料的方法，如矩法、累积量法、最优化方法等；二是依据流域地形和地貌资料的方法，如经验公式法等；三是依据流域出口断面流量资料的方法本文着重介绍应用较为广泛的矩法根据矩法，可以求出参数:[MQ1)MP]2nQ2)n(2)(3-10)nQ2)n(2)mQ1〉m『(3-11)式中：Mi⑴、MQ1)分别为入流量(净雨量)及出流量的一阶原点矩；N⑵、NQ2)分别为入流量和出流量的二阶中心矩入流量和出流量的一阶原点矩的计算公式分别为(1)I0I(t)tdt0I(t)dtn_「i(t)tii1n_「i(t)i1li(t)miii1Ii(t)i1(3-12)(1)Q0Q(t)tdtQi(t)tii1Qi(t)miii1(3-13)0Q(t)dtQi(t)Qi(t)式中：l(t)、l[(t)分别为t时刻的净雨量和净雨量的时段平均值；Q(t)、Qi(t)分别为t时刻的出流量和出流量的时段平均值。

mi1,3,5,7,,2n3,2n1同理，可以计算出净雨量和出流量的二阶原点矩为Mi⑵n「i(t)m2i1n「i(t)i1(3-14)nQ(t)m2T2mQ2)-i1(3-15)nQi(t)i1由于入流量和出流量的二阶中心矩可用原点矩来表示，即n(2)M|⑵[M(1).2I](3-16)nQ2)MQ2)[M(1)]2QJ(3-17)在推求某一流域纳须时段单位线时，首先选择几场有代表性的洪水资料，根据每场次洪水的相应入流和出流过程，分别求出它们的单位线参数n和k，然后把n和k平均(如果相差较大，可以分成几组)概化出流域的单位线参数，便可用纳须时段单位线计算公式计算出流域时段单位线应用实例某站的一次降雨径流过程下表，依据资料利用矩法求纳须时段单位线时间净雨(mm)实测出流(m3/s)基流(m3/s)地面径流m3/s7,9,210810808516109407141000。