1 / 19

Hydrologic River Flood Routing

Hydrologic River Flood Routing. Basic Equation. Muskingum Method. For a non-uniform flow , S t  f(O t ) Because I t =f(y 1 ) and O t = f(y 2 ) . Therefore, S t =f(I t ,O t ) . Prism storage = f(O t ); Wedge storage = f(I t , O t ).

sylvia
Download Presentation

Hydrologic River Flood Routing

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. Hydrologic River Flood Routing Basic Equation Routing

  2. Muskingum Method • For a non-uniform flow, Stf(Ot) • Because It=f(y1)and Ot= f(y2). Therefore, St=f(It,Ot). • Prism storage = f(Ot); Wedge storage = f(It, Ot). • Assume: Prism storage = KOt; Wedge storage = XK (It – Ot) • Total storage volume at any time instant, t, is St = KOt + XK (It - Ot) = K [X It + (1 - X)Ot] where K = travel time of flood wave in a reach, [hrs]; X = weighting factor, 0 ~ 0.5 Routing

  3. Continuity Equation in Difference Form • Referring to figure, the continuity equation in difference form can be expressed as Routing

  4. Derivation of Muskingum Routing Equation • By Muskingum Model, at t = t2, S2 = K [X I2 + (1 - X)O2] at t = t1, S1 = K [X I1 + (1 - X)O1] • Substituting S1, S2 into thecontinuity equation and after some algebraic manipulations, one has O2 = Co I2 + C1 I1 + C2 O1 • Replacing subscript 2 by t +1 and 1 by t, the Muskingum routing equation is Ot+1 = Co It+1 + C1 It + C2 Ot, for t = 1, 2, … where ; ; C2 = 1 – Co – C1 Note: K and t must have the same unit. Routing

  5. Muskingum Routing Procedure • Given (knowns): O1; I1, I2, …; t; K; X • Find (unknowns): O2, O3, O4, … • Procedure: (a) Calculate Co, C1, and C2 (b) Apply Ot+1 = Co It+1 + C1 It + C2 Otstarting from t=1, 2, … recursively. Routing

  6. Muskingum Routing - Example Routing

  7. Estimating Muskingum Parameters (1) Graphical Method: • Referring to the Muskingum Model, find X such that the plot of XIt+ (1-X)Ot vs St behaves almost nearly as a single value curve. The corresponding slope is K. Routing

  8. Example (Graphical Method) Routing

  9. Estimating Muskingum Parameters (2) • Gill’s Least Square Method (Gill, 1978, “Flood routing by Muskingum method”, J. of Hydrology, V.36) St = K[XIt+(1-X)Ot] = KX It+K(1-X)Ot = AIt+BOt where A = KX, B = K(1-X) • Model: St’ = So + A It + B Ot, t = 1, 2, …, n where St’ = cumulative relative storage at time t; So = initial storage of channel before flood • Use least squares method to solve for A and B and So. Then, K and X can be obtained as K = A + B, X = A / (A + B) Routing

  10. Example (Gill’s Method) Routing

  11. Estimating Muskingum Parameters (3) • Stephensen’s Least Squares Method (Stephensen, 1979, “Direct optimization of Muskingum routing coefficient, “ J. of Hydrology, V.41) • Ot+1 = Co It+1 + C1 It + C2 Ot, t = 1, 2, …, n-1 • Use least squares method to estimate the values of Co, C1 and C2 directly. (Does not guarantee satisfying Co + C1 + C2 = 1) • Alternatively,  Co + C1 + C2 =1, C2 = 1 – (Co + C1) (Ot+1 - Ot) = Co (It+1 - Ot) + C1 (It - Ot), t = 1, 2, …, n  ; Routing

  12. Example (Stephensen Method) Routing

  13. Disadvantages of Hydrologic Routing (1) Ignore the dynamic effect of flow; (2) Assume stage and storage is a single–valued function of discharge - implying flow is changing slowly with time. Routing

  14. Nonlinear Muskingum Model • In some stream reaches, St and XIt + (1-X)Ot reveals pronounced nonlinear relation. • To accommodate the non-linearity, the Muskingum model can be modified as St = K[XIt + (1-X)Ot]m or St = K[XItp + (1-X)Otq] • Using the nonlinear relationship, the routing becomes more difficult. A procedure using state-variable technique was developed by Tung (1985) to perform channel routing using nonlinear Muskingum model. Routing

  15. Nonlinear St vs XIt + (1-X)Ot Relation Routing

  16. Muskingum Crest Segment Routing • Purposes: • solve for single outflow rate or • route only a segment of inflow hydrograph • Basis: for t =1, 2, … Since By substituting Ot into the routing equation for Ot+1 and, after some algebraic manipulations, one has From recursive substitution of outflow at one time in terms of that at the previous time, one could have where K1=C0, K2 = C0 C2 + C1; Ki = Ki-1 C2, i≥3 Routing

  17. Distributed Flow Routing • Estimate flow rate or water surface elevation at different locations and time simultaneously, rather than separately at different locations, so that the unsteady, non-uniform nature of actual flow phenomena are more closely modeled. • Modeling of flow movement can be in 1-D, 2-D, or 3-D in space and in time, depending on the dominant flow velocity field to be modeled. Routing

  18. Governing Eqs. For Flow Routing(Saint-Venant Equations) Routing

  19. Steady vs Unsteady / Uniform vs Non-uniform Routing

More Related