首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
All river engineering schemes require flood discharge estimates as part of the design and appraisal process. Unfortunately, continuous measurement of flood discharges is limited to those river sites with instrumented gauging stations, which constitute only a small proportion of channel reaches where information is required. Therefore, considerable research effort has been devoted to the development of reliable indirect techniques of flood discharge estimation. Research on the interrelationship of stream channel geometry and river discharge has provided the basis for an indirect method of flood estimation – the channel-geometry method – which employs river channel dimensions alone to estimate discharge characteristics at ungauged river sites. Channel-geometry equations are developed empirically by relating streamflow data from gauging stations and channel dimensions measured from natural river reaches in the vicinity of the gauge, and take the form of power function relations. Once regional channel-geometry equations have been defined, a channel width or channel capacity measurement is the only variable needed to estimate the flood flow characteristics at a specified river site. The method is useful as an alternative to traditional catchment-based approaches or as a rapid reconnaissance technique. In addition to the application for flood discharge prediction, channel-geometry equations could prove helpful in the management of river channels, first, by providing a basis for assessing local deviations in the channel form–discharge relation, deviations which could be employed as indicators of the sensitivity of particular stretches of river channel to change, and secondly, in the computation of natural channel dimensions for use in river channel design and river restoration.  相似文献   

2.
The braiding intensity and dynamics in large braiding rivers are well known to depend on peak discharges, but the response in braiding and channel–floodplain transformations to changes in discharge regime are poorly known. This modelling study addresses the morphodynamic effects of increasing annual peak discharges in braiding rivers. The study site is a braiding reach of the Upper Yellow River. We estimated the effects on the larger‐scale channel pattern, and on the smaller‐scale bars, channel branches and floodplains. Furthermore, we determined the sensitivity of the channel pattern to model input parameters. The results showed that the dominant effect of a higher peak discharge is the development of chute channels on the floodplains, formed by connecting head‐cut channels and avulsive channels. Widening of the main channel by bank erosion was found to be less dominant. In addition, sedimentation on the bars and floodplains increased with increasing peak discharge. The model results also showed that the modelled channel pattern is especially sensitive to parametrization of the bed slope effect, whereas the effect of median grain size was found to be relatively small. Copyright © 2018 John Wiley & Sons, Ltd.  相似文献   

3.
Nicholas Pinter 《水文研究》2010,24(8):1088-1093
This study tests the hypothesis that historical float‐based discharge measurements on the Mississippi River systematically over‐stated actual flood flows by 10% to > 30% relative to measurements using current meters. This assertion has been repeated over the past 25 years and recently has been used to adjust historical discharges used for flood‐frequency analysis. This study tests the hypothesis above using 2150 historical discharge measurements digitized from the three principal gauging stations on the Middle Mississippi River (MMR): data that include 626 float‐based discharges and 1516 meter‐based discharges, including 122 paired measurements. Multiple comparative tests show that the hypothesis above cannot be supported; if anything, the float‐based measurements slightly underestimate flows (not over‐estimate) over a broad range of discharges up to large floods. In response to the purported data bias above (‘changing history’; Dieckmann RJ, Dyhouse GR. 1998. Changing history at St. Louis—adjusting historic flows for frequency analysis. First Federal Inter‐Agency Hydrologic Modeling Conference, April 20–22, 1998. Las Vegas, NV; 4·31–4·36), historical flood discharges on the MMR have been modified, most by 10–20% and several by > 30%. These altered discharges are now being promulgated, in particular, through the Upper Mississippi River System Flow Frequency Study (UMRSFFS). New flow frequencies, flood profiles, and new flood maps from the UMRSFFS may significantly underestimate the actual flood hazard on the MMR if the original hydrologic data have been erroneously altered on the basis of an assumption of data bias. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

4.
ABSTRACT

An accurate comprehension of celerity (flood wave speed) dynamics is a key step for understanding flood wave propagation in rivers. We present the results of empirically estimated celerity values in 12 Brazilian rivers, and analyse the behaviour of celerity–discharge relationships (CxQ). Celerity was estimated with a reach-scale (RS) method, based on the peak travel time between stations; and with a local-scale (LS) method, based on the derivative of discharge–cross-section area relationships surveyed at gauging stations. The results indicate that the magnitudes of celerity values obtained by the methods are reasonably comparable, and can rarely be considered constant, varying with river discharge. Three reaches presented differing CxQ relationships at local and reach scales, which suggests that in situ cross-sections at gauging stations should not be extrapolated as representative of the whole reach for flood routing studies, and that CxQ relationship assessments might provide relevant insights for hydrological modelling.  相似文献   

5.
A well-developed subglacial drainage system consisting of large cavities developed in the lee of bedrock steps connected together by a network of Nye channels is exposed on an area of recently deglaciated limestone bedrock in front of Glacier de Tsanfleuron, Switzerland. This system covers some 51 per cent of the bedrock surface area, and is believed to have transported the bulk of supraglacially-derived meltwaters through the glacier. Using the cavity hydraulics model of Kamb (1987), it is shown that the geometry of the system rendered it stable against collapse by meltback of channel roofs into a tunnel-dominated system. For likely combinations of glacier geometry and meltwater discharge, the steady state water pressure in this system would have been only a small fraction of that required for flotation, and for discharges of less than about 0·5–5 m3 s?1 water would have flowed at atmospheric pressure. The system appears to have adjusted to varying discharges by a combination of varying water pressure and changing the total cross-sectional area of flow by altering the number of active channels connecting cavities. Glacier sliding velocity would have been independent of meltwater discharge for discharges at which water flowed at atmospheric pressure, but would have risen with increasing discharge for higher flows. Velocities on the order of 0·1 m d?1 are predicted for a realistic range of discharges and effective pressures, and these are believed to be plausible. Episodes of enhanced sliding in glaciers with similar drainage systems could be triggered by a rise in meltwater discharge across the threshold between flows at atmospheric pressure and flow under pressure from the glacier.  相似文献   

6.
Bankfull discharge was identified in some 30 gravel-bed rivers representing in total c. 40 gauging stations. The catchment sizes cary from 4km2 to nearly 2700km2. Bankfull discharge value increases with basin size. In the case of gravel-bed rivers developed on an impermeable substratum, the following equation emerges: Qb=0·087 A1·044. Bankfull discharge recurrence interval was determined by fitting maximum annual floods (Ta) into Gumbel's distribution and then using the partial duration series (Tp) in this same distribution. Recurrence interval is below 0·7 years (Tp) for small pebble-bed rivers developed on an impermeable substratum; it reaches 1·1 to 1·5 years when the catchment size of these rivers exceeds 250km2. Rivers incised in the soft schists of the Famenne show larger channel capacity at bankfull stage, a small width/depth ratio and thus higher recurrence intervals (1·4–5·3 years with Ta and 1–4·4 years with Tp). Baseflow-dominated gravel-bed streams and sandy or silty rivers experience less frequent bankfull discharges, with a recurrence interval higher than 2 or even 3 years (Tp). © 1997 John Wiley & Sons, Ltd.  相似文献   

7.
This paper reports on an evaluation of the use of artificial neural network (ANN) models to forecast daily flows at multiple gauging stations in Eucha Watershed, an agricultural watershed located in north‐west Arkansas and north‐east Oklahoma. Two different neural network models, the multilayer perceptron (MLP) and the radial basis neural network (RBFNN), were developed and their abilities to predict stream flow at four gauging stations were compared. Different scenarios using various combinations of data sets such as rainfall and stream flow at various lags were developed and compared for their ability to make flow predictions at four gauging stations. The input vector selection for both models involved quantification of the statistical properties such as cross‐, auto‐ and partial autocorrelation of the data series that best represented the hydrologic response of the watershed. Measured data with 739 patterns of input–output vector were divided into two sets: 492 patterns for training, and the remaining 247 patterns for testing. The best performance based on the RMSE, R2 and CE was achieved by the MLP model with current and antecedent precipitation and antecedent flow as model inputs. The MLP model testing resulted in R2 values of 0·86, 0·86, 0·81, and 0·79 at the four gauging stations. Similarly, the testing R2 values for the RBFNN model were 0·60, 0·57, 0·58, and 0·56 for the four gauging stations. Both models performed satisfactorily for flow predictions at multiple gauging stations, however, the MLP model outperformed the RBFNN model. The training time was in the range 1–2 min for MLP, and 5–10 s for RBFNN on a Pentium IV processor running at 2·8 GHz with 1 MB of RAM. The difference in model training time occurred because of the clustering methods used in the RBFNN model. The RBFNN uses a fuzzy min‐max network to perform the clustering to construct the neural network which takes considerably less time than the MLP model. Results show that ANN models are useful tools for forecasting the hydrologic response at multiple points of interest in agricultural watersheds. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

8.
《Journal of Hydrology》2003,270(1-2):105-134
Freshwater availability has been recognized as a global issue, and its consistent quantification not only in individual river basins but also at the global scale is required to support the sustainable use of water. The WaterGAP Global Hydrology Model WGHM, which is a submodel of the global water use and availability model WaterGAP 2, computes surface runoff, groundwater recharge and river discharge at a spatial resolution of 0.5°. WGHM is based on the best global data sets currently available, and simulates the reduction of river discharge by human water consumption. In order to obtain a reliable estimate of water availability, it is tuned against observed discharge at 724 gauging stations, which represent 50% of the global land area and 70% of the actively discharging area. For 50% of these stations, the tuning of one model parameter was sufficient to achieve that simulated and observed long-term average discharges agree within 1%. For the rest, however, additional corrections had to be applied to the simulated runoff and discharge values. WGHM not only computes the long-term average water resources of a country or a drainage basin but also water availability indicators that take into account the interannual and seasonal variability of runoff and discharge. The reliability of the modeling results is assessed by comparing observed and simulated discharges at the tuning stations and at selected other stations. The comparison shows that WGHM is able to calculate reliable and meaningful indicators of water availability at a high spatial resolution. In particular, the 90% reliable monthly discharge is simulated well. Therefore, WGHM is suited for application in global assessments related to water security, food security and freshwater ecosystems.  相似文献   

9.
Eroding channels can usually be characterized by a power relationship between channel width (W) and channel discharge (Q). This paper examines the WQ relation using a recently developed channel junction approach to extend the validity of the WQ relation and to develop a procedure for estimating the WQ exponent and proportionality coefficient. Rill and gully channel data from the literature, and new data collected in different badland areas and in a few forest mountain streams, are analysed. Analysis shows that the WQ relation for channel width collected in badlands and forests agrees with trends observed for cropland. The exponent increases with increasing channel width in a continuous fashion rather than in a step‐like way and tends to a maximum whose value ranges between 0·5 and 0·6. The proportionality coefficient can be split into two terms, one expressing the case in which an eroding channel can broaden, the other reflecting the difficulties in removing the less erodible clods or rock fragments from the channel bed. Its splitting allows the development of a more correct form of the WQ relation in agreement with modern approaches of channel geometry: one part has the dimension of a discharge and makes the power base dimensionless, while the other brings the dimension of a length, needed for the channel width, into the WQ relation. The interpretation of the two constants is supported by data collected in rainfall‐runoff simulation experiments conducted in the field. Values characterizing the two constants in some environments are also given. Nevertheless the approach is not sufficiently parameterized yet to be of practical use (e.g. in models or for estimating peak discharge in areas where rill channels have formed). Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

10.
Alpine glacial basins are a significant source and storage area for sediment exposed by glacial retreat. Recent research has indicated that short‐term storage and release of sediment in proglacial channels may control the pattern of suspended sediment transfer from these basins. Custom‐built continuously recording turbidimeters installed on a network of nine gauging sites were used to characterize spatial and temporal variability in suspended sediment transfer patterns for the entire proglacial area at Small River Glacier, British Columbia, Canada. Discharge and suspended sediment concentration were measured at 5 min intervals over the ablation season of 2000. Differences in suspended sediment transfer patterns were then extracted using multivariate statistics (principal component and cluster analysis). Results showed that each gauging station was dominated c. 80% of days by diurnal sediment transfer patterns and ‘low’ suspended sediment concentrations. ‘Irregular’ transfer patterns were generally associated with ‘high’ sediment concentrations during snowmelt and rainfall events, resulting in the transfer of up to 70% of the total seasonal suspended sediment load at some gauging stations. Suspended sediment enrichment of up to 600% from channel storage release and extrachannel inputs occurred between the glacial front and distal proglacial boundary. However, these patterns differed significantly between gauging stations as determined by the location of the gauging station within the catchment and meteorological conditions. Overall, the proglacial area was the source for up to 80% of the total suspended sediment yield transferred from the Small River Glacier basin. These results confirmed that sediment stored and released in the proglacial area, in particular from proglacial channels, was controlling suspended sediment transfer patterns. To characterize this control accurately requires multiple gauging stations with high frequency monitoring of suspended sediment concentration. Accurate characterization of this proglacial control on suspended sediment transfer may therefore aid interpretation of suspended sediment yield patterns from glacierized basins. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

11.
Stream response to surface coal mining and reclamation was studied in 29 small (0·13 to 5·72 km2) watersheds located in the bituminous coal fields of Central Pennsylvania. These basins, up to 82 per cent mined, were selected from 176 first-order tributaries of Beech Creek with similar vegetation, soil, lithology, and basin characteristics. Measurements were made at 262 cross-sections (an average of nine cross-sections per stream) of channel cross-section area, bankfull width, mean bankfull depth, dimensions of the largest moving blocks, stream slope, valley-side slope, basin area, and mined area. Observed differences in channel morphology were related to differences in extent of mining by means of scatter plots, correlation, cluster analysis, and bivariate regression. Stream response to increased peak discharge and channel shear stress produced by increased surface runoff from regraded mine spoil takes the form of enlarged channels and increases in the size of moving blocks. Large basin areas appear to dampen the effect of mining, resulting in limited channel enlargement with greater extent of mining. In contrast, where peak discharges and associated shear stresses exceed the combined erosional resistance of floodplain vegetation, colluvial blocks, and channel banks, streams adjust extensively to higher levels of mining, causing an abrupt increase in the size of transported blocks and eroded channels. In the first-order basins studied, this stepped response occurs at approximately 0·45 km2 mined area and 50 per cent of the total basin area mined. For streams that have exceeded both threshold levels, disequilibrium is demonstrated by a strong, positive correlation between local stream slope and basin area. Where both threshold levels of mining are exceeded, steep channel slopes reinforce the tendency of stream cross-sections to increase with greater disturbance by mining, necessitating that these streams rapidly adjust their morphology in order to attain a new equilibrium which is compatible with the conditions imposed by mining and reclamation.  相似文献   

12.
A system has been installed to automatically monitor rainfall, streamflow, bedload discharge and suspended sediment concentration in the arid to hyper‐arid setting of Nahal Rahaf, Southern Judean Desert in Israel. The Rahaf gauging station is located in a relatively steep, wide channel with an unsteady bed driven by flash floods. It is an attempt to deploy modern automatic equipment for continuous sediment transport monitoring in harsh, arid fluvial environments. Unit bedload discharges are the highest recorded hitherto, suggesting they may represent an upper end member in the worldwide climate–bedload discharge relationship. Suspended sediment concentration is much higher than is typical of perennial fluvial humid environments. There is high correlation between suspended sediment concentration and water discharge on an event scale, with diverse intra‐event relations. The sediment yield of individual events is large, but the small number of floods limits the mean annual sediment yield to low values in this arid environment. This also has environmental implications, as large‐scale quarrying requires a long period of self‐restoration in such an arid fluvial setting. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

13.
Adjustment of stream-channel shape to hydrologic regime   总被引：1，自引：0，他引：1  
Bed-load channels tend to adjust their cross-sections so that given slope, roughness and sediment load, channel shape approaches the optimum for bed-load transport. The extent to which any one shape is the optimum varies with discharge, so four Cumberland Basin stream channels have been investigated to determine the discharges at which their present cross-sections represent the optimum for bed-load transport. These discharges have return periods ranging from 1.1 to 1.5 years on the annual series. The return periods closely correspond with return periods for the discharge at which, over a period of time, the most bed-load is transported. These return periods vary from 1.15 to 1.45 years when the same bed-load equation is used. The close correspondence between sets of return periods suggests that bed-load channels tend to adjust their cross-sections to become the optimum shape for bed-load transport at or close to the discharge at which the most bed-load transport is accomplished.  相似文献   

14.
Springs along the south rim of the Grand Canyon, Arizona, are important ecological and cultural resources in Grand Canyon National Park and are discharge points for regional and local aquifers of the Coconino Plateau. This study evaluated the applicability of electrical resistance (ER) sensors for measuring diffuse, low-stage (<1.0 cm) intermittent and ephemeral flow in the steep, rocky spring-fed tributaries of the south rim. ER sensors were used to conduct a baseline survey of spring flow timing at eight sites in three spring-fed tributaries in Grand Canyon. Sensors were attached to a nearly vertical rock wall at a spring outlet and were installed in alluvial and bedrock channels. Spring flow timing data inferred by the ER sensors were consistent with observations during site visits, with flow events recorded with collocated streamflow gauging stations and with local precipitation gauges. ER sensors were able to distinguish the presence of flow along nearly vertical rock surfaces with flow depths between 0.3 and 1.0 cm. Laboratory experiments confirmed the ability of the sensors to monitor the timing of diffuse flow on impervious surfaces. A comparison of flow patterns along the stream reaches and at springs identified the timing and location of perennial and intermittent flow, and periods of increased evapotranspiration.  相似文献   

15.
The seventy-kilometre-long Herbert Gorge of northeastern Australia preserves a record of past floods in slackwater deposits and palaeostage indicators. Step-backwater modelling of water-surface profiles indicates that discharges ranging from 11000 to 17000 m3s?1 have occurred six times in the gorge during the last 900 years. These flood reconstructions provide insight into the role of extreme flows in shaping bedrock channel morphology. In particular, the hydraulics of extreme flows can be related to boulder transport, and to the location of large boulder bars. Large boulder bars occur throughout the Herbert Gorge, being best developed at loci of stream power minima along the inside of bends, at tributary junctions, and at obstructions in the channel caused by bedrock highs. Only the flows exceeding approximately 8000 m3 s?1 are competent to transport the boulders which constitute the bars. In the straight channel reaches, the boulder accumulations and bedrock highs have a fairly regular spacing which appears to be independent of lithologic or structural controls. The bars provide an efficient means of energy dissipation, and they are interpreted as a result of the inherent high turbulence of flow in a steep channel. The regular spacing of the bars, and their correspondence with the hydraulics of large flows, suggest that the bars and associated bedrock highs may represent a self-regulating mechanism akin to the pool-riffle sequence of alluvial channels. It may therefore be appropriate to view bedrock channels as deformable on the timescale of extreme discharges.  相似文献   

16.
The potential impacts of climate change can alter the risk to critical infrastructure resulting from changes to the frequency and magnitude of extreme events. As well, the natural environment is affected by the hydrologic regime, and changes in high flows or low flows can have negative impacts on ecosystems. This article examines the detection of trends in extreme hydrological events, both high and low flow events, for streamflow gauging stations in Canada. The trend analysis involves the application of the Mann–Kendall non‐parametric test. A bootstrap resampling process has been used to determine the field significance of the trend results. A total of 68 gauging stations having a nominal record length of at least 50 years are analysed for two analysis periods of 50 and 40 years. The database of Canadian rivers investigated represents a diversity of hydrological conditions encompassing different extreme flow generating processes and reflects a national scale analysis of trends. The results reveal more trends than would be expected to occur by chance for most of the measures of extreme flow characteristics. Annual and spring maximum flows show decreasing trends in flow magnitude and decreasing trends in event timing (earlier events). Low flow magnitudes exhibit both decreasing and increasing trends. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

17.
Parsimonious stage–fall–discharge rating curve models for gauging stations subject to backwater complications are developed from simple hydraulic theory. The rating curve models are compounded in order to allow for possible shifts in the hydraulics when variable backwater becomes effective. The models provide a prior scientific understanding through the relationship between the rating curve parameters and the hydraulic properties of the channel section under study. This characteristic enables prior distributions for the rating curve parameters to be easily elicited according to site‐specific information and the magnitude of well‐known hydraulic quantities. Posterior results from three Norwegian and one American twin‐gauge stations affected by variable backwater are obtained using Markov chain Monte Carlo simulation techniques. The case studies demonstrate that the proposed Bayesian rating curve assessment is appropriate for developing rating procedures for gauging stations that are subject to variable backwater. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

18.
This study investigated the effect of urbanization on runoff from the On-Cheon Stream watershed in Pusan, Korea. This watershed has been experiencing considerable urbanization since the 1960s. There are two gauging stations in the watershed. For one of the stations there are recent flow data and for the other flow data were observed in the past. A linear reservoir model was chosen and runoff was analysed for several flood events. The linear reservoir model has been found to generate flood hydrographs accurately for both gauging stations, and its applicability to the study area has also been established. Using two methods of computing effective rainfall or rainfall excess (ϕ-index and constant percentage method), the results of runoff analyses were investigated. The ϕ-index method yielded better results than the constant percentage method. A comparison of hydrographs observed in the past with the simulation results at the Ie-Sup bridge site revealed that the peak discharge increased and the mean lag time of the study area decreased owing to urbanization over the past two decades. It is also possible to evaluate the effect of urbanization quantitatively. © 1998 John Wiley & Sons, Ltd.  相似文献   

19.
This paper compares hydrologic records and geomorphic effects of several historic floods in the central Appalachian region of the eastern United States. The most recent of these, occurring in November 1985, was the largest ever recorded in West Virginia, with peak discharges exceeding the estimated 500-year discharge at eight of eleven stations in the South Branch Potomac River and Cheat River basins. Geomorphic effects on valley floors included some of the most severe and widespread floodplain erosion ever documented and exceeded anything seen in previous floods, even though comparable or greater rainfall and unit discharge have been observed several times in the region over the past 50 years. Comparison of discharge-drainage area plots suggests that the intensity and spatial scale of the November 1985 flood were optimal for erosion of valley floors along the three forks of the South Branch Potomac River. However, when a larger geographic area is considered, rainfall totals and discharge-drainage area relationships are insufficient predictors of geomorphic effectiveness for valley floors at drainage areas of 250 to 2500 km2. Unit stream power was calculated for the largest recorded flood discharge at 46 stations in the central Appalachians. Maximum values of unit stream power are developed in bedrock canyons, where the boundaries are resistant to erosion and the flow cross-section cannot adjust its width to accommodate extreme discharges. The largest value was 2570 W m?2; record discharge at most stations was associated with unit stream power values less than 300 W m?2, but more stations exceeded this value in the November 1985 flood than in the other floods that were analysed. Unit stream power at indirect discharge measurement sites near areas experiencing severe erosion in this and other central Appalachian floods generally exceeded 300 W m?2; reach-average values of 200-500 W m?2 were calculated for valleys where erosion damage was most widespread. Despite these general trends, unit stream power is not a reliable predictor of geomorphic change for individual sites. Improved understanding of flood impacts will require more detailed investigation of interactions between local site characteristics and patterns of flood flow over the valley floor.  相似文献   

20.
After its formation, a rill may remain in the field for months, often receiving lower flow rates than the formative discharge. The objective of this work was to evaluate the rill flow transport capacity of soil aggregates at discharges unable to erode the rill, and to analyse the influence of the rill macro‐roughness on this transport process. A non‐erodible rill was built in which roughness was reproduced in detail. In order to assess only the rill macro‐roughness, a flat channel with a similar micro‐roughness to that in the rill replica was built. Rill and channel experiments were carried out at a slope of 8 and at six discharges (8·3 × 10?5 to 5·2 × 10?4 m3 s?1) in the rill, and eight discharges (1·6 × 10?5 to 5·2 × 10?4 m3 s?1) in the channel. Non‐erodible aggregates of three sizes (1–2, 3–5 and 5–10 mm) were released at the inlet of the rill/channel. The number of aggregates received at the outlet was registered. The number and position of the remaining aggregates along the rill/channel were also determined. The rill flow was a major sediment transport mechanism only during the formation of the rill, as during that period the power of the flow was great enough to overcome the influence of the macro‐roughness of the rill bed. At lower discharges the transport capacity in the previously formed rill was significantly less than that in the flat channel under similar slope and discharge. This was determined to be due to local slowing of flow velocities at the exit of rill pools. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号