第三章 单位线分析计算
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
第三章 单位线分析计算
本章着重介绍如何由净雨量过程预报流域出口的流量过程。净雨量经过流域汇流形成出口的流量过程线,流域汇流历时是降雨径流预报预见期的来源,流域汇流物理过程是编制预报方案的理论依据。 舍尔曼时段单位线 3.1.1 基本原理
舍尔曼(L. K. Sherman)于1932年提出了单位线的概念。其定义为:流域上分布均匀的1个单位净雨直接径流产流量,所形成的直接径流过程线,即为单位线,记为UH 。
1个单位净雨是指单位时段内单位净雨深。单位时段长可以任取,例如2h 、3h 、6h ,等。而单位净雨深通常取为10mm 。
而实际发生的净雨,常常不是1个时段,也不是1个单位,应用于分析单位线时,有一些假定。这些假定可归纳为以下两点:
(1) 如果单位时段内净雨深不是一个单位,而是n 个,它所形成的出流过程线,总历时与UH 相同,流量则是UH 的n 倍。
(2) 如果净雨历时不是一个时段,而是m 个,则各个时段净雨所形成的出流过程之间互不干扰,出流过程的流量过程等于m 个流量过程之和。
由以上假定,净雨d r 、出流d Q 与UH 的纵坐标q 之间的关系如下:
(
∑=+-=m
i i t i d t d q r Q 1
1,, (2-1)
式中:m i ,3,2,1 =,为净雨时段数。d Q 和q 的单位为s /m 3,d r 则用单位净雨深的n 倍来表示。
如果UH 已知,根据上式,可由净雨转换为出流,计算十分简便。关键是如何求得UH 。可以根据流域的实测水文资料,分析出净雨及直接径流过程后,依据上式推求出UH ,为式(2-1)的逆过程。 3.1.2 单位线的推求
推求UH 是使用次洪时段净雨深及相隔为计算时段长的直接径流时序过程。
前者由次洪降雨量经过扣损后得到,后者由径流过程线分割地下水后得到。这里需要补充说明一下具体问题:
(1) 由扣损方案求得的次洪净雨深,常不等于过程线分割得到的实测值,为了不把扣损的误差带入汇流计算,需要将计算值改正,或谓平差。在平差前应分析误差的来源,作出较为合理的修正,以不改变原来的雨型为原则。
(2) 计算时段长t ∆的确定,主要考虑峰形的控制,一般取涨洪历时的1/3~1/4,常与流域大小成正比
对于多时段净雨而言,传统的单位线分析方法有以下一些: (1) 分析法
由 ∑=+-=m
i i t i d t d q r Q 1
1,,,得到
《
4)-(3 3)-(3 2)-(3 13,22,31,3,12,21,2,11,1,q r q r q r Q q r q r Q q r Q d d d d d d d d d ++=+==
因此,上式组成的方程组为一个多元线性代数方程组,求解该代数方程组可以得到1q 、2q 的数值。最简单的解法是逐一消去法。由上面的方程组可以看出,由式(3-2),已知1,d Q 和1,d r ,解出1q :
1
,1,1d d r Q q =
(3-5)
1q 已知,将其带入式(3-3),可得到
1
,1
2,2,2d d d r q r Q q -=
(3-6)
如此递推下去,得到
1
,2
1
,,d m
i i t i d t d t r q r Q q ∑=+--=
(3-7)
计算结果见表3-1和图3-1所示。
直接代数法对于降雨与径流实测资料没有误差,流域汇流符合线性时不变系统时,能得出正确的唯一解。但实际情况是:提供分析UH 的实测资料存在着
观测与分析误差,流域汇流为非线性系统,因此直接分析解法不但不能得出唯一
解,由于误差累积,其解常很不合理。如分析得到的t
q~呈锯齿形,如表3-1和图3-1的退水段,有的甚至出现负值,无法继续计算下去。所以很少有人使用。
图3-1 分析法得到的单位线
/
表3-1 单位线分析法结果
(2) 试错法
这个方法是假设单位线,目估对比推流的与实测的t
Q
~。当两者最接近时,
d
所假设的UH即为所求。
初始的UH,可用其它洪水已分析得来的成果,或用斜线分割法的结果,也可任意假定。
试错法应用比较广泛,有的用科林法试错,有的凭经验试错。但主要确定是单位线过程的初始假定比较困难,试错过程有时也会出现不合理现象,对其过程进行修正也不易做得完好。
(3) 系统识别方法
—
UH既然是线性系统的单位响应,进一步可以应用线性系统鉴别的方法推求UH的最优解,如有约束的最小二乘法和回归法等。
(4) 各种方法的特点
如上所述,推求单位线的方法主要有分析法、试错法和系统识别法。试错法又可分为目估试错法和科伦试错法,后者是一种迭代性质的试错法。
分析法当降雨为1个时段时能取得很好的效果,但当降雨时段数大于1个时段,分析法推求的结果可能出现锯齿型,甚至出现负值,主要原因是误差的累积。目估试错法是通过人工目估,逐步试错,使计算的流量过程与实测的流量过程吻合,避免了误差的累积,但任意性大,吻合程度缺乏客观标准,并且要求调试者具有相当的经验,费时费力。科伦试错法是首先拟定初始单位线,然后按初始单位线计算出除最大降雨时段外的其它时段降雨产生的流量过程,并从实测流量过程中将其减去,差值即当成最大降雨时段产生的流量过程,据此推求新的单位线,将新单位线作为初始单位线,重复上述步骤,直至新旧单位线之间的差值达到给定的误差要求。优点是当初始单位线拟定后,不必目估试错,能自动逐步修正,适合计算机编程。缺点是计算结果及迭代的收敛性与初始单位线的拟定有关,推求的结果不一定是最佳结果。另外当降雨时段较多时,迭代难以继续下去。