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

第三章 单位线分析计算 本章着重介绍如何由净雨量过程预报流域出口的流量过程。净雨量经过流域汇流形成出口的流量过程线，流域汇流历时是降雨径流预报预见期的来源，流域汇流物理过程是编制预报方案的理论依据。 3.1 舍尔曼时段单位线 3.1.1 基本原理 舍尔曼(L. K. Sherman)于1932年提出了单位线的概念。其定义为：流域上分布均匀的1个单位净雨直接径流产流量，所形成的直接径流过程线，即为单位线，记为。 1个单位净雨是指单位时段内单位净雨深。单位时段长可以任取，例如2h、3h、6h，等。而单位净雨深通常取为10mm。 而实际发生的净雨，常常不是1个时段，也不是1个单位，应用于分析单位线时，有一些假定。这些假定可归纳为以下两点： 如果单位时段内净雨深不是一个单位，而是个，它所形成的出流过程线，总历时与UH相同，流量则是UH的倍。 如果净雨历时不是一个时段，而是个，则各个时段净雨所形成的出流过程之间互不干扰，出流过程的流量过程等于个流量过程之和。 由以上假定，净雨、出流与的纵坐标之间的关系如下： (2-1) 式中：，为净雨时段数。和的单位为，则用单位净雨深的倍来表示。 如果已知，根据上式，可由净雨转换为出流，计算十分简便。关键是如何求得。可以根据流域的实测水文资料，分析出净雨及直接径流过程后，依据上式推求出，为式(2-1)的逆过程。 3.1.2 单位线的推求 推求是使用次洪时段净雨深及相隔为计算时段长的直接径流时序过程。前者由次洪降雨量经过扣损后得到，后者由径流过程线分割地下水后得到。这里需要补充说明一下具体问题： 由扣损方案求得的次洪净雨深，常不等于过程线分割得到的实测值，为了不把扣损的误差带入汇流计算，需要将计算值改正，或谓平差。在平差前应分析误差的来源，作出较为合理的修正，以不改变原来的雨型为原则。 (2) 计算时段长的确定，主要考虑峰形的控制，一般取涨洪历时的1/3～1/4，常与流域大小成正比 对于多时段净雨而言，传统的单位线分析方法有以下一些： (1) 分析法 由 ，得到 因此，上式组成的方程组为一个多元线性代数方程组，求解该代数方程组可以得到、的数值。最简单的解法是逐一消去法。由上面的方程组可以看出，由式(3-2)，已知和，解出： (3-5) 已知，将其带入式(3-3)，可得到 (3-6) 如此递推下去，得到 (3-7) 计算结果见表3-1和图3-1所示。 直接代数法对于降雨与径流实测资料没有误差，流域汇流符合线性时不变系统时，能得出正确的唯一解。但实际情况是：提供分析的实测资料存在着观测与分析误差，流域汇流为非线性系统，因此直接分析解法不但不能得出唯一解，由于误差累积，其解常很不合理。如分析得到的呈锯齿形，如表3-1和图3-1的退水段，有的甚至出现负值，无法继续计算下去。所以很少有人使用。 图3-1 分析法得到的单位线 表3-1 单位线分析法结果 日 时 时序 计算值 修正值 计算值 7 6 0 0 0 0 0 12 1 186 24.5 76 76 186 18 2 667 20.3 209 210 669 8 0 3 1935 616 617 1938 6 4 2450 489 490 2453 12 5 1900 370 355 1864 18 6 1280 216 240 1309 9 0 7 850 168 155 867 6 8 560 89 105 572 12 9 400 89 73 392 18 10 277 39 52 276 10 0 11 202 50 38 199 6 12 142 16 22 131 12 13 80 19 12 74 18 14 40 1 0 24 11 0 15 0 0 0 0 (2) 试错法 这个方法是假设单位线，目估对比推流的与实测的。当两者最接近时，所假设的即为所求。 初始的，可用其它洪水已分析得来的成果，或用斜线分割法的结果，也可任意假定。 试错法应用比较广泛，有的用科林法试错，有的凭经验试错。但主要确定是单位线过程的初始假定比较困难，试错过程有时也会出现不合理现象，对其过程进行修正也不易做得完好。 (3) 系统识别方法 既然是线性系统的单位响应，进一步可以应用线性系统鉴别的方法推求的最优解，如有约束的最小二乘法和回归法等。 (4) 各种方法的特点 如上所述，推求单位线的方法主要有分析法、试错法和系统识别法。试错法又可分为目估试错法和科伦试错法，