马斯京根法(比赛终).ppt

1第三节第三节 马斯京根法马斯京根法(Muskingum Method)水文学原理水文学原理（Principle of Hydrology）2主要内容主要内容1 基本原理和概念基本原理和概念2 参数和演算时段参数和演算时段3 分段流量演算法分段流量演算法 基本原理基本原理4 非线性流量演算非线性流量演算演算时段的确定演算时段的确定参数的推求参数的推求本节课本节课讲解的内容讲解的内容 马斯京根流量演算马斯京根流量演算3几个问题几个问题1 为什么要学习马斯京根法？为什么要学习马斯京根法？2 什么是什么是“马斯京根法马斯京根法”？3 如何推求马斯京根法的参数？如何推求马斯京根法的参数？4第二章第二章 河段洪水预报河段洪水预报 1 基本概念和原理(concept and theory)u马斯京根法是由马斯京根法是由G.T.麦卡锡于麦卡锡于1938年年提出，因首先应提出，因首先应用于美国马斯京根河而得名。用于美国马斯京根河而得名。u马斯京根法是进行河道流量演算的水文学方法。马斯京根法是进行河道流量演算的水文学方法。u方法虽然简单但有一定的实用价值，所以至今仍被广泛方法虽然简单但有一定的实用价值，所以至今仍被广泛的应用。的应用。5流量演算法是利用水量平衡原理与流量演算法是利用水量平衡原理与河段中的蓄（蓄水量）泄（出流量）河段中的蓄（蓄水量）泄（出流量）关系，将河段内上断面的入流过程关系，将河段内上断面的入流过程演算成下断面的出流过程。演算成下断面的出流过程。水量平衡方程水量平衡方程(water balance equation)槽蓄方程槽蓄方程(channel storage equation)（1）河道洪水演算）河道洪水演算(river flood routing)连续方程连续方程动力方程动力方程Q(t)W Q1Q2t1t2I1I2tQ,I0t图图4-1 河段水量平衡图河段水量平衡图I(t)1 基本概念和原理(concept and theory)6流量演算法本质流量演算法本质流量演算法联解上两式，后者因不同洪水受附加比降的影响各异，相应流量演算法联解上两式，后者因不同洪水受附加比降的影响各异，相应的蓄泄关系也不相同。如果蓄泄关系呈单值线性函数形式，流量演算可的蓄泄关系也不相同。如果蓄泄关系呈单值线性函数形式，流量演算可大为简化。因此，寻求槽蓄曲线呈单值线性函数是河道洪水演算水文学大为简化。因此，寻求槽蓄曲线呈单值线性函数是河道洪水演算水文学方法讨论的主要内容。方法讨论的主要内容。水量平衡方程水量平衡方程槽蓄方程槽蓄方程（1）河道洪水演算）河道洪水演算(river flood routing)1 基本概念和原理(concept and theory)7基本假定：基本假定：(1)河段中的水面是一直线，流量沿程变化是线性的。河段中的水面是一直线，流量沿程变化是线性的。Q与与I，Q线性关系线性关系 W是是Q的线性函数的线性函数I为入流为入流蓄泄关系曲线的坡度蓄泄关系曲线的坡度即传播时间即传播时间Q为出流为出流（1）河道洪水演算）河道洪水演算(river flood routing)QQ/基本假定：基本假定：(2)河段中的某断面流量河段中的某断面流量Q/与河段的蓄水量与河段的蓄水量W之间为单一的线性关系。之间为单一的线性关系。流量比重系数流量比重系数1 基本概念和原理(concept and theory)8（2）流量演算公式）流量演算公式(flow routing)水量平衡方程水量平衡方程槽蓄方程槽蓄方程C0C1C2+=11 基本概念和原理(concept and theory)9（2）流量演算公式）流量演算公式(flow routing)输入输入参数参数K、x计算时段计算时段t入流入流I和和初始出流量初始出流量输出输出1 基本概念和原理(concept and theory)10（3）流量演算实例）流量演算实例已知长江万县已知长江万县宜昌河段的参数宜昌河段的参数x=0.25，K=t=18h，要求，要求将万县流量演算到宜昌流量过程。将万县流量演算到宜昌流量过程。1 基本概念和原理(concept and theory)11已知长江万县已知长江万县宜昌河段的参数宜昌河段的参数x=0.25，K=t=18h，要求，要求将万县流量演算到宜昌流量过程。将万县流量演算到宜昌流量过程。（3）流量演算实例）流量演算实例1 基本概念和原理(concept and theory)12（3）流量演算实例）流量演算实例时间时间 (月，日，时月，日，时)万县实测入万县实测入流量流量I I(m(ms)s)0.260.26I I2 20.480.48I I1 10.260.26Q Q1 1宜昌演算出宜昌演算出流量流量Q Qc c(m(ms)s)（1）（2）（3）（4）（5）（6）7.1.1419900 22800 7.2.824300 7.3.238800 7.3.2050000 7.4.1453800 7.5.850800 7.6.243400 7.6.2035100 7.7.1426900 7.8.822400 7.9.219600 7.9.2017900 631895525928100881300013988132081128491266994582450964654217891166418624240002582424384208321075216848129129408566771291007612497133971275772341110690866001230822006327419387534806451529490654271534948278221 基本概念和原理(concept and theory)13（3）流量演算实例）流量演算实例1 基本概念和原理(concept and theory)142 参数和演算时段(Parameter and Calculus duration)（1）参数的推求）参数的推求已知参数已知参数K、x流量演算流量演算流量演算的流量演算的关键在于推关键在于推求参数求参数K、x 试算法(trial method)分析法(analytical method)最小二乘法(least square method)最常用最常用假设一个假设一个x x值值按方程按方程Q=Qr+x(I-Qr)求出相应的一组求出相应的一组Q值值按方程按方程 求逐时段的槽蓄变量求逐时段的槽蓄变量W累加累加W得各时段末的得各时段末的W(t)做QW(t)的关系曲线的关系曲线WQQWQ/W(t)为线性关系？为线性关系？x为所求：为所求：K=W/QYN推求推求x,K计算流程图计算流程图（1）参数的推求）参数的推求2 参数和演算时段(Parameter and Calculus duration)16时间时间 万县实测万县实测入流量入流量I宜昌实测宜昌实测出流量出流量Q区间径区间径流量流量q修正后的修正后的实测出流实测出流量量Qr=Q-qQ=I-QrWW(ms)(月，日，时月，日，时)(m(ms)s)(m(ms)s)(m(ms)s)(m(ms)s)(m(ms)s)(m(ms)s)18h18h(m(ms)s)18h18hx=0.1=0.1x=0.250.25x=0.150.15（1）（2）（3）（4）（5）（6）（7）（8）（9）（10）（11）7.1.1419900 7.2.824300 237006007.3.238800 2700016007.3.2050000 3780012007.4.1453800 484009007.5.850800 519005007.6.243400 496004007.6.2035100 430004007.7.1426900 356004007.8.822400 293003007.9.219600 242003007.9.2021300200（1）参数的推求）参数的推求Q=Qr+x(I-Qr)2 参数和演算时段(Parameter and Calculus duration)17时间时间 万县实万县实测入流测入流量量I宜昌实宜昌实测出流测出流量量Q区间径区间径流量流量q修正后修正后的实测的实测出流量出流量Qr=Q-qQ=I-QrWW(ms)(月，日，月，日，时时)(m(ms)s)(m(mss)(m(ms)s)(m(ms)s)(m(ms)s)(m(ms)s)18h18h(m(ms)s)18h18hx=0.1=0.1 x=0.250.25 x=0.150.15（1）（2）（3）（4）（5）（6）（7）（8）（9）（10）（11）7.1.14 19900 7.2.8 24300 23700 6007.3.2 38800 27000 16007.3.20 50000 37800 1200（1）参数的推求）参数的推求Q=Qr+x(I-Qr)2310025400366001200 13400 13400 730013400 9850 07300 20700 2322026740379402340028750399502328027410386102 参数和演算时段(Parameter and Calculus duration)18时间时间 万县实测万县实测入流量入流量I宜昌实测宜昌实测出流量出流量Q区间径区间径流量流量q修正后的修正后的实测出流实测出流量量Qr=Q-qQ=I-QrWW(ms)(月，日，时月，日，时)(m(ms)s)(m(ms)s)(m(ms)s)(m(ms)s)(m(ms)s)(m(ms)s)18h18h(m(ms)s)18h18hx=0.1=0.1x=0.250.25x=0.150.15（1）（2）（3）（4）（5）（6）（7）（8）（9）（10）（11）7.1.1419900 7.2.824300 23700600231001200730002322023400232807.3.238800 270001600254001340013400 7300 2674028750274107.3.2050000 37800120036600134009850 20700 3794039950386107.4.1453800 484009004750063002850 30550 4813049075484457.5.850800 5190050051400-600-3200 33400 5134051250513107.6.243400 4960040049200-5800-6650 30200 4862047750483307.6.2035100 4300040042600-7500-7900 23550 4185040725414757.7.1426900 3560040035200-8300-7450 15650 3437033125339557.8.822400 2930030029000-6600-5450 8200 2834027350280107.9.219600 2420030023900-4300-4300 2750 2347022825232557.9.202130020021100（1）参数的推求）参数的推求Q=Qr+x(I-Qr)2 参数和演算时段(Parameter and Calculus duration)19QWK=W/Q1(即即K=18h)2 参数和演算时段(Parameter and Calculus duration)20时间时间 万县实测入流量万县实测入流量I 宜昌实测出流量宜昌实测出流量Q 区间径流量区间径流量q 宜昌演算出流量宜昌演算出流量Qc(月，日，时月，日，时)(m(ms)s)(m(ms)s)(m(ms)s)(m(ms)s)（1）（2）（3）（4）（5）7.1.1419900 217897.2.824300 23700600274197.3.238800 270001600387537.3.2050000 378001200480647.4.1453800 48400900515297.5.850800 51900500490657.6.243400 49600400427157.6.2035100 43000400349487.7.1426900 35600400278227.8.822400 29300300230827.9.219600 242003007.9.2021300200（1）参数的推求）参数的推求2 参数和演算时段(Parameter and Calculus duration)212 参数和演算时段(Parameter and Calculus duration)（1）参数的推求）参数的推求假定本河段的假定本河段的x、K值已选定，即值已选定，即x=0.15，K=18h以及以及t=18h,带入下式：带入下式：2 参数和演算时段(Parameter and Calculus duration)23KK、x x受原始资料干扰受原始资料干扰受原始资料干扰受原始资料干扰河段内水面曲线非线性河段内水面曲线非线性河段内水面曲线非线性河段内水面曲线非线性 问问问问河段内水力特性非线性河段内水力特性非线性河段内水力特性非线性河段内水力特性非线性讨讨讨讨题题题题各次洪水参数不同各次洪水参数不同各次洪水参数不同各次洪水参数不同论论论论参参参参数数数数的的的的物物物物理理理理意意意意义义义义不不不不明明明明确确确确（1）参数的推求）参数的推求p 存在问题：存在问题：QW曲线虽然合拢，但非直线，成曲线虽然合拢，但非直线，成“8”字形字形2 参数和演算时段(Parameter and Calculus duration)24（1）参数的推求）参数的推求(1)用本法确定的参数，有时因所取计算时段不同而有差别，宜作分析比较。(2)应选择区间径流尽可能少的洪水作为分析对象，以减小区间水对参数确定的影响。(3)作为一次完整洪水，应注意始、末流量基本相等。(4)在计算区间洪水总量时，应考虑河段汇流时间，在具体处理时可参照洪峰传播时间或其他方法进行估算。(5)在区间来水分配时，应考虑河段区间面积汇流特性。本例区间来水是以单位线为基础进行分配，是比较合理的。试算法应注意的问题2 参数和演算时段(Parameter and Calculus duration)25中国知网中国知网http:/ http:/ 参考资料参考资料(Refrence)26Exercise根据某河段一次实测洪水资料，用马斯京根法进行河段洪水演算。根据某河段一次实测洪水资料，用马斯京根法进行河段洪水演算。时间时间 上上 下下 q区区 下下-q区区 Q上上+q区区-Q下下 W W Q/(m3/s)（月（月.日日.时）时）（m3/s）(m3/s)(m3/s)(m3/s)(m3/s)(m3/s)12h(m3/s)12h x=0.2 x=0.3 1 2 3 4 5678910.75 75 0.370 80 37.1620 440 73.2210 1680 110.2290 2150 73.1830 2280 37.1220 1680 0.830 1270 .610 880 .480 680 .390 550 .330 450 .300 400 .260 340 .230 290 .200 250 .180 220 .160 200 作作 业业27Thank you！三亚租车三亚租车 http:/