马斯京根洪水演算方程参数确定的最小二乘法.pdfVIP

维普资讯 水文计算 · 马斯京根供水演算方程参数确定的最小二乘法 马 斯 京 根 洪 水 演 算 方 程 参 数 确 定 的 最 小 二 乘 法 黄 才 安，- (扬州大学水利学院 摘要 本文根据线性 回归的最小二乘法原理 ．提 出了一什由安测湛水 资料 直接确定 马斯京丰量 潜水漓算 方程参数 c C 、C 的方法 畦方法与 目前广泛苯用 的试错法杷比 ，具有计革的参数精度 较 高、计算过程可 以程序化等优点 。 关键词 辨水演算 马斯京根法 参敷计算 最小二柬法 马斯京根 (Muskingum)洪水演算法 由于具有原理简单，计算简便 ．而且计算结果的精度可 以满足实际要求等优点 因此 ，它是 目前应用较为普遍的洪水演算方法之一。马靳京根洪水演 算方程为 “： 0。+1一CoL+1+C1Ii C：0． (1) 式 中：of、Q+为时段始 、末的河段下游断面出流量 ，I．．L 为时段始 、末的河段上辨断面人流量 ， C。、C。、C：为洪水演算系数，可 由下列各式确定： 0．5△ t—kx … ， 、 L 一 —k-- kx+—O SA t ’ ’ ‘ ． r1 一 kx+。_垒 r¨ 一 k--kx+0．5△ t … k--kx--O．5△ t … L 一 —k -- kx+—0 SA t ． 式中，△t为时段长，k为槽蓄系数 ，x为流量 比重因子。 为 了利用式(2)～ (4)求得洪水演算系数 c。、c 、c：，必然首先确定其中的△t、k、x为此 ，除 演算时段长△t可根据问题的要求确定外 ，槽蓄系数 k和流量 比重因子 x则需要根据实测洪 水资料进行估算。目前 ．试错法。是估算k和X值的主要方法。在应用试错法求参数的计算过 程 中，需要数次假定x值，计算相应的示储流量Q =xI+ (1一x)O，然后点绘槽蓄量 s= (I— O)与Q 的关系曲线 ．由目估法判断该 曲线是否为单值 曲线 ，以确定假定之x值是否正确 。若假 定之x正确，由该线的坡度定出k值 。再 由所确定的k、x值，利用式 (2)～(4)求洪水演算系数 C…C 、C：，从而确定了马斯京根洪水演算方程 。显然 ，这样的参数求解过程存在着：①计算工作 量较大，而且不容易程序化，难 以由计算机完成计算工作；②参数计算的精度较差 ，特别是存在 人为因素对参数精度韵影响。 有鉴于此，本文利用最小二乘法原理，直接 由实测洪水资料确定方程 (1)中的参数 c。、c 、 C ， 克服试错法 中的不足 。 ·30 · 维普资讯 第i3卷 1906 《第 2甥 2