水文学原理课件(第十章).ppt

'第十章洪水演算(FloodRouting) 本章主要内容引言(Introduction)1线性扩散波演算法(Lineardiffusionwaveroutingmethod)2线性特征河长演算法(Linearcharacteristiclengthroutingmethod)3线性运动波演算法(Linearkinematicwaveroutingmethod)4 Part1Introduction1、定义(Definition)洪水演算：根据河段上断面出现的洪水过程来推算其下断面洪水过程的方法称洪水演算。2、分类(Classification)具有物理基础的洪水演算法—主要用于河道汇流基于“黑箱子”分析的洪水演算法—主要用于流域汇流 Part1Introduction3、具有物理基础的洪水演算法概述(Physically-basedfloodroutingmethod)水力学水文学水力学方法的关键是如何求解方程式。水文学方法的关键是如何确定槽蓄方程。 Part2Lineardiffusionwaveroutingmethod定解问题的构成(Compositionofsolutionproblems)1入流过程的处理(Processingofinflowhydrograph)2基本解(Basicsolution)3出流过程的计算(Derivationofoutflowhydrograph)4参数的确定方法(Determinationofparameters)5 Part2Lineardiffusionwaveroutingmethod2-1定解问题的构成(Compositionofsolutionproblems)洪水波—扩散波（Diffusionwave）初始条件—稳定流（Steadyflow）上边界过程—上断面洪水过程（inflowhydrographatuppersection)下边界条件—半无限长，自由下边界（semi—infinitelong,freelowerboundary） Part2Lineardiffusionwaveroutingmethod2-2Processingofinflowhydrograph(1)简单入流过程(Simpleinflowhydrograph)a)单位入流过程(UnitInflowhydrograph)：一直维持1个单位流量的入流过程(例如1m3/s)延迟单位入流 Part2Lineardiffusionwaveroutingmethod2-2Processingofinflowhydrographb)单位矩形入流过程(UnitRectangularPulseInput)在单位时间(△t)内维持1个单位流量的入流称之.延迟单位矩形入流(1)简单入流过程(Simpleinflowhydrograph) Part2Lineardiffusionwaveroutingmethod2-2Processingofinflowhydrographc)单位瞬时脉冲入流(UnitinstantaneousPulseInput)时段很短（），强度很大（），但入流总量为有限数（等于1个单位）的入流称之.(1)简单入流过程(Simpleinflowhydrograph) Part2Lineardiffusionwaveroutingmethod2-2Processingofinflowhydrograph(2)入流过程的数学处理(ProcessingofInflowHydrograph)用柱状图代替光滑的I(t)曲线.用减小△t的方法提高近似精度.第一个矩形：第二个矩形：以此类推。 Part2Lineardiffusionwaveroutingmethod2-3Basicsolution简单入流的解称为基本解，又称为(河段)汇流曲线(flowconcentrationcurve)。如果河段入流为单位入流函数H(t)，则下断面出流为S曲线，又称为单位过程线(UnitHydrograph)如果某河段上断面入流为单位瞬时脉冲入流δ(t)，则河段下断面出流称为瞬时单位线u(0,t)，简写为u(t)如果河段上断面入流为单位矩形入流，则下断面出流称为时段单位线u(△t,t) Part2Lineardiffusionwaveroutingmethod2-3Basicsolution一般情况下,先求u(0,t)，再求S(t)，最后求u(△t,t)u(0,t)S(t)u(△t,t) Part2Lineardiffusionwaveroutingmethod2-3Basicsolution应用Laplace变换，容易求得瞬时单位线的表达式为： Part2Lineardiffusionwaveroutingmethod2-3Basicsolution若：应用Laplace变换，求得S(t)曲线的表达式：式中—余误差函数(complementaryerrorfunction) Part2Lineardiffusionwaveroutingmethod2-4Derivationofoutflowhydrograph首先把入流过程离散化(Discretizinginflowhydrograph)用单位矩形入流表示 Part2Lineardiffusionwaveroutingmethod2-4Derivationofoutflowhydrograph计算出流量过程的主要步骤：1)确定计算时段△t。2)按在上断面入流过程线上从起涨时刻开始摘取时刻的入流量3)计算时段平均入流量或时段平均入流量的增量。4)按求定的河段的时段单位线或S—曲线，计算出流过程。 Part2Lineardiffusionwaveroutingmethod2-4Derivationofoutflowhydrograph现以按u(△t,t)推流为例，说明计算Q(t)-t过程。 Part2Lineardiffusionwaveroutingmethod2-4Derivationofoutflowhydrographt-τ为汇流时间离散化(Discretization)卷积公式(Convolutionformula) Part2Lineardiffusionwaveroutingmethod2-5Determinationofparameters(1)根据参数的物理公式计算 Part2Lineardiffusionwaveroutingmethod2-5Determinationofparameters(2)根据实测流量资料确定（试算法）直到假定的使得推算的与实测足够接近。只要满足方程组：但由于的表达式很复杂，直接解出较困难。最优化方法： Part3Linearcharacteristiclengthroutingmethod描述洪水波运动的基本微分方程式(Basicdifferentialequationsoffloodwavemovement)1汇流曲线的确定(Determinationofflowconcentrationcurve)2参数确定方法(Determinationofparameters)3 Part3Linearcharacteristiclengthroutingmethod3-1Basicdifferentialequationsoffloodwavemovement在每个子河段内，满足:每个子河段： Part3Linearcharacteristiclengthroutingmethod3-1Basicdifferentialequationsoffloodwavemovement先不对每个河段写方程，而对整个河段写方程，则由上式可得：消去所有中间变量得到: Part3Linearcharacteristiclengthroutingmethod3-1Basicdifferentialequationsoffloodwavemovement从而 Part3Linearcharacteristiclengthroutingmethod3-2Determinationofflowconcentrationcurve假设每一特征河长段传播时间相等，则则方程变为：设初始条件为零，求入流为时的出流,此情况如下：用Laplace变换，得： Part3Linearcharacteristiclengthroutingmethod3-3Determinationofparameters Part4Linearkinematicwaveroutingmethod运动波差分方程的建立(Differenceequationofkinematicwave)1数值扩散的概念(Numericaldiffusion)2汇流系数的计算(Calculationofflowconcentrationcoefficient)3连续演算问题(successiverouting)4讨论(Discussions)5 Part4Linearkinematicwaveroutingmethod4-1Differenceequationofkinematicwave当运动波为线性运动波（即波速C恒定时），运动波既不变形，也不坦化，像“刚体运动”一般。即运动波在实际河流中很少发生，其解析解意义不大，但微分改成差分后带入了误差项，正好可模拟扩散项，所以其差分解法更有实际意义。,运动波运动特征：运动波 Part4Linearkinematicwaveroutingmethod4-1Differenceequationofkinematicwave对j断面：对j-1断面： Part4Linearkinematicwaveroutingmethod4-1Differenceequationofkinematicwave采取折中办法加权平均表示和 Part4Linearkinematicwaveroutingmethod4-2Numericaldiffusion以差商代替微商带来了截断误差，若能将误差考虑进来，则近似程度可得到提高。截断误差可由泰勒公式(Taylorformula)分析得到。则有： Part4Linearkinematicwaveroutingmethod4-2Numericaldiffusion若取二阶精度的截断误差，则有：虽然运动波方程本身不能描述具有衰减的洪水波，但具有二阶精度的运动波差分方程在一定条件下可达到这个目的，条件为： Part4Linearkinematicwaveroutingmethod4-2Numericaldiffusion扩散波方程： Part4Linearkinematicwaveroutingmethod4-2Numericaldiffusion若已知上断面出流，求下断面出流：(Cungemethod) Part4Linearkinematicwaveroutingmethod4-3Calculationofflowconcentrationcoefficient1)入流过程离散化(seeppt)(连续Q~t曲线称为过程，离散化的Q~t称为时间序列)延迟单位序列：单位序列所形成的下断面的出流称为汇流系数(flowconcentrationcoefficient)单位序列：只在某时刻有一个单位强度(1m3/s)的入流称单位序列，用δ(n)表示。 Part4Linearkinematicwaveroutingmethod4-3Calculationofflowconcentrationcoefficient首先确定然后求当n=0时，=1……P:汇流系数 Part4Linearkinematicwaveroutingmethod4-3CalculationofflowconcentrationcoefficientTimeI(n)P1(n)P.IQ(n)I(1)I(2)I(3)183I0P1(0)I0P1(0)I0P1(0)6I1P1(1)I0P1(1)I1P1(0)I0P1(1)+I1P1(0)9I2P1(2)I0P1(2)I1P1(1)I2P1(0)I0P1(2)+I1P1(1)+I2P1(0)12P1(3)I0P1(3)I1P1(2)I2P1(1)I0P1(3)+I1P1(2)+I2P1(1)……………… Part4Linearkinematicwaveroutingmethod4-4Successiverouting隐式差分—无条件稳定显式差分—有条件稳定—满足柯朗条件，即河段长正好为因此，为了保证差分格式的稳定性,对于河长大于Δχ的河段,需要按Δχ长度将河段划分为若干子河段.差分格式的稳定性问题： Part4Linearkinematicwaveroutingmethod4-5Discussions假定槽蓄方程的形式为： Part4Linearkinematicwaveroutingmethod4-5Discussions则上式可改变为：上式与运动波的二阶精度的差分方程等价.以上就是1932年提出的Muskingum流量算法。美国陆军工程兵团(USArmyCorpsofEngineers)首次在美国的马斯京根河上应用.1969年Frenchengineer(法国工程师)Cunge(康吉)给出了Muskingum的理论基础，并给出了x的物理意义，即： 水文学原理PrincipleofHydrology本章结束，谢谢！'