【正文】
附 錄附錄一:外文原文Sensitivity Analysis of the CCHE1D Channel Network Model Weiming Wu (1), Dalmo A. Vieira (2), Abdul Khan (3) and Sam S. Y. Wang (4) (1), (2), (3) and (4), National Center for Computational Hydroscience and Engineering, School of Engineering, The University of Mississippi, MS 38677。 PH (662) 9155673 / (662) 9157788。 FAX (662) 9157796。 Email: wuwm@ AbstractThe CCHE1D model was designed to simulate longterm flow and sediment transport in channel networks to support the DEC project. It uses either the dynamic wave or the diffusive wave model to pute unsteady flows in channel networks with pound cross sections, taking into account the effects of instream hydraulic structures, such as culverts, weirs, drop structures, and bridge crossings. It simulates nonuniform sediment transport using a nonequilibrium approach, and calculates bank toe erosion and mass failure due to channel incision. The CCHE1D model decouples the flow and sediment transport calculations but couples the calculations of nonuniform sediment transport, bed changes and bed material sorting in order to enhance the numerical stability of the model. In this paper, the sensitivity of CCHE1D to parameters such as the nonequilibrium adaptation length of sediment transport and the mixing layer thickness is evaluated in cases of channel aggradation and degradation in laboratory flumes as well as in a natural channel network. In the case of channel degradation, the simulated scour process is not sensitive to variation in values of the nonequilibrium adaptation length, but the determination of the mixing layer thickness is important to the putations of the equilibrium scour depth and of the bedmaterial size distribution at the armoring layer. The simulated bed profiles in the case of channel aggradation and the calculated sediment yield in the case of natural channel network are insensitive to the prescription of both the nonequilibrium adaptation length and the mixing layer thickness. The CCHE1D model can provide reliable results even when these two parameters are given a wide range of values. Introduction The CCHE1D model was designed to simulate longterm flow and sediment transport in channel networks to support the Demonstration Erosion Control (DEC) project, which is an interagency cooperative effort among the US Army Corps of Engineers (COE), the Natural Resources Conservation Service (NRCS) and the Agricultural Research Service (ARS) of the US Department of Agriculture. The CCHE1D version was based on the unsteady flow model DWAVNET (Diffusion WAVe model for channel NETworks, Langendeon, 1996) and the sediment transport model BEAMS (Bed and Bank Erosion Analysis Model for Streams, Li et al., 1996). It was significantly improved by implementing the dynamic wave model and the nonequilibrium sediment transport model (Wu, Vieira and Wang 2000). The CCHE1D was integrated with the landscape analysis tool TOPAZ (Garbrecht and Martz, 1995) and with the watershed models AGNPS (Bosch et al., 1998) and SWAT (Arnold et al., 1993), through an ArcView GISbased graphical user interface (Vieira and Wu, 2000). The CCHE1D has been successfully tested in various experimental and field cases. Because several parameters in CCHE1D must be prescribed empirically, it is very important to know the response of the model to the uncertainty of these parameters. In this study, the sensitivity of CCHE1D to model parameters such as the nonequilibrium adaptation length of sediment transport and the mixing layer thickness is analyzed in cases of channel aggradation and degradation in laboratory flumes as well as in a natural channel network. Description of the CCHE1D Channel Network Model Hydrodynamic Model. The CCHE1D flow model simulates unsteady flow in channel networks with pound crosssections using either the diffusive wave model or the dynamic wave model. The dynamic wave model solves the full St. Venant equations. The Preissmann implicit,fourpoint, finite difference scheme is used to discretize the governing equations. Linearized iteration schemes for the discretized governing equations are established and solved using a double sweep algorithm. The influence of hydraulic structures such as culverts, measuring flumes, bridge crossings and drop structures has been considered in the CCHE1D model. Stagedischarge relations for hydraulic structures are derived so that the hydraulic structures bee an intrinsic part of the numerical algorithm. Sediment Transport Model. The CCHE1D model calculates nonuniform sediment transport in rivers using a nonequilibrium approach. The governing equation for the nonequilibrium transport of nonuniform total load is (1)where A is the flow area。 Ctk is the depthaveraged totalload concentration of size class k。 Qtk is the actual totalload transport rate。 Qt*k is the totalload transport capacity。 Ls is the adaptation length of nonequilibrium sediment transport。 and qlk is the side sediment discharge from banks or tributaries per unit channel length, with the contribution from banks being simulated by CCHE1D bank erosion and bank failure module, and the contribution from upland erosion being simulated by SWAT or AGNPS. The sediment transport capacity can be written as a general form (2)where pbk is the bed material gradation。 Q*tk is the potential sediment transport rate, which is determined with SEDTRA module (Garbrecht et al., 1995), Wu, Wang and Jia’s formula (2000), the modified Ackers and White’s 1973 formula (Proffitt and Sutherland, 1983), or the modified Engelund and Hansen’s 1967 formula (with Wu, Wang and Jia’s correction factor, 2000). The bed deformation due to size class k is determined with (3)where p′ is the bed material porosity, which is calculated with t