首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 62 毫秒
1.
I. MUZIK 《水文研究》1996,10(10):1401-1409
The concept of a spatially distributed unit hydrograph is based on the fact that the unit hydrograph can be derived from the time–area curve of a watershed by the S-curve method. The time–area diagram is a graph of cumulative drainage area contributing to discharge at the watershed outlet within a specified time of travel. Accurate determination of the time–area diagram is made possible by using a GIS. The GIS is used to describe the connectivity of the links in the watershed flow network and to calculate distances and travel times to the watershed outlet for various points within the watershed. Overland flow travel times are calculated by the kinematic wave equation for time to equilibrium; channel flow times are based on the Manning and continuity equations. To account for channel storage, travel times for channel reaches are increased by a percentage depending on the channel reach length and geometry. With GIS capability for rainfall mapping, the assumption of a uniform spatial rainfall distribution is no longer necessary; hence the term, spatially distributed unit hydrograph. An example of the application for the Waiparous Creek in the Alberta Foothills is given. IDRISI is used to develop a simple digital elevation model of the 229 km2 watershed, using 1 km × 1 km grid cells. A grid of flow directions is developed and used to create an equivalent channel network. Excess rainfall for each 1 km × 1 km cell is individually computed by the Soil Conservation Service (SCS) runoff curve method and routed through the equivalent channel network to obtain the time–area curve. The derived unit hydrograph gave excellent results in simulating an observed flood hydrograph. The distributed unit hydrograph is no longer a lumped model, since it accounts for internal distribution of rainfall and runoff. It is derived for a watershed without the need for observed rainfall and discharge data, because it is essentially a geomorphoclimatic approach. As such, it allows the derivation of watershed responses (hydrographs) to inputs of various magnitudes, thus eliminating the assumption of proportionality of input and output if needed. The superposition of outputs is retained in simulating flood hydrographs by convolution, since it has been shown that some non-linear systems satisfy the principle of superposition. The distributed unit hydrograph appears to be a very promising rainfall runoff model based on GIS technology.  相似文献   

2.
ROGER MOUSSA 《水文研究》1997,11(5):429-449
Recently, several attempts have been made to relate the hydrological response of a catchment to its morphological and topographical features using different hypotheses to model the effect of the drainage network. Several transfer functions were developed and some of these are based on the theory of a linear model, the geomorphological unit hydrograph. The aim of this paper is to present a methodology to automatically identify the transfer function, using digital elevation models for applications in distributed hydrological modelling. The transfer function proposed herein is based on the Hayami approximation solution of the diffusive wave equation especially adapted for the routing hydrograph through a channel network. The Gardon d’Anduze basin, southern France, was retained for applications. Digital elevation models were used to extract the channel network and divide the basin into subcatchments. Each subcatchment produces, at its own outlet, an impulse response which is routed to the outlet of the whole catchment using the diffusive wave model described by two parameters: celerity and diffusivity functions of geometrical characteristics of the channel network. Firstly, a geomorphological unit hydrograph obtained by routing a homogeneous effective rainfall was compared with the unit hydrograph identified by a lumped model scheme, then the distributed model was applied to take into account the spatial variability of effective rainfall in the catchment. Results show that this new method seems to be adapted for distributed hydrological modelling; it enables identification of a transfer function response for each hydrological unit, here subcatchments, and then simulation of the contribution of each unit to the hydrograph at the outlet. © 1997 by John Wiley & Sons, Ltd.  相似文献   

3.
Distributed, continuous hydrologic models promote better understanding of hydrology and enable integrated hydrologic analyses by providing a more detailed picture of water transport processes across the varying landscape. However, such models are not widely used in routine modelling practices, due in part to the extensive data input requirements, computational demands, and complexity of routing algorithms. We developed a two‐dimensional continuous hydrologic model, HYSTAR, using a time‐area method within a grid‐based spatial data model with the goal of providing an alternative way to simulate spatiotemporally varied watershed‐scale hydrologic processes. The model calculates the direct runoff hydrograph by coupling a time‐area routing scheme with a dynamic rainfall excess sub‐model implemented here using a modified curve number method with an hourly time step, explicitly considering downstream ‘reinfiltration’ of routed surface runoff. Soil moisture content is determined at each time interval based on a water balance equation, and overland and channel runoff is routed on time‐area maps, representing spatial variation in hydraulic characteristics for each time interval in a storm event. Simulating runoff hydrographs does not depend on unit hydrograph theory or on solution of the Saint Venant equation, yet retains the simplicity of a unit hydrograph approach and the capability of explicitly simulating two‐dimensional flow routing. The model provided acceptable performance in predicting daily and monthly runoff for a 6‐year period for a watershed in Virginia (USA) using readily available geographic information about the watershed landscape. Spatial and temporal variability in simulated effective runoff depth and time area maps dynamically show the areas of the watershed contributing to the direct runoff hydrograph at the outlet over time, consistent with the variable source area overland flow generation mechanism. The model offers a way to simulate watershed processes and runoff hydrographs using the time‐area method, providing a simple, efficient, and sound framework that explicitly represents mechanisms of spatially and temporally varied hydrologic processes. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

4.
In this paper a very general rainfall-runoff model structure (described below) is shown to reduce to a unit hydrograph model structure. For the general model, a multi-linear unit hydrograph approach is used to develop subarea runoff, and is coupled to a multi-linear channel flow routing method to develop a link-node rainfall-runoff model network. The spatial and temporal rainfall distribution over the catchment is probabilistically related to a known rainfall data source located in the catchment in order to account for the stochastic nature of rainfall with respect to the rain gauge measured data. The resulting link node model structure is a series of stochastic integral equations, one equation for each subarea. A cumulative stochastic integral equation is developed as a sum of the above series, and includes the complete spatial and temporal variabilities of the rainfall over the catchment. The resulting stochastic integral equation is seen to be an extension of the well-known single area unit hydrograph method, except that the model output of a runoff hydrograph is a distribution of outcomes (or realizations) when applied to problems involving prediction of storm runoff; that is, the model output is a set of probable runoff hydrographs, each outcome being the results of calibration to a known storm event.  相似文献   

5.
Attempts to reduce the number of parameters in distributed rainfall–runoff models have not yet resulted in a model that is accurate for both natural and anthropogenic hillslopes. We take on the challenge by proposing a distributed model for overland flow and channel flow based on a combination of a linear response time distribution and the hillslope geomorphologic instantaneous unit hydrograph (GIUH), which can be calculated with only a digital elevation model and a map with field boundaries and channel network as input. The spatial domain is subdivided into representative elementary hillslopes (REHs) for each of which we define geometric and flow velocity parameters and compute the GIUH. The catchment GIUH is given by the sum of all REH responses. While most distributed models only perform well on natural hillslopes, the advantage of our approach is that it can also be applied to modified hillslopes with for example a rectangular drainage network and terrace cultivation. Tests show that the REH‐GIUH approach performs better than classical routing functions (exponential and gamma). Simulations of four virtual hillslopes suggest that peak flow at the catchment outlet is directly related to drainage density. By combining the distributed flow routing model with a lumped‐parameter infiltration model, we were also able to demonstrate that terrace cultivation delays the response time and reduces peak flow in comparison to the same hillslope, but with a natural stream network. The REH‐GIUH approach is a first step in the process of coupling distributed hydrological models to erosion and water quality models at the REH (associated with agricultural management) and at the catchment scale (associated with the evaluation of the environmental impact of human activities). It furthermore provides a basis for the development of models for large catchments and urban or peri‐urban catchments. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

6.
In this paper a very general rainfall-runoff model structure (described below) is shown to reduce to a unit hydrograph model structure. For the general model, a multi-linear unit hydrograph approach is used to develop subarea runoff, and is coupled to a multi-linear channel flow routing method to develop a link-node rainfall-runoff model network. The spatial and temporal rainfall distribution over the catchment is probabilistically related to a known rainfall data source located in the catchment in order to account for the stochastic nature of rainfall with respect to the rain gauge measured data. The resulting link node model structure is a series of stochastic integral equations, one equation for each subarea. A cumulative stochastic integral equation is developed as a sum of the above series, and includes the complete spatial and temporal variabilities of the rainfall over the catchment. The resulting stochastic integral equation is seen to be an extension of the well-known single area unit hydrograph method, except that the model output of a runoff hydrograph is a distribution of outcomes (or realizations) when applied to problems involving prediction of storm runoff; that is, the model output is a set of probable runoff hydrographs, each outcome being the results of calibration to a known storm event.  相似文献   

7.
The physical basis of the linkage between magnitude and timing of channel flow hydrographs and drainage network morphometry is reviewed. Small Hortonian and structurally Hortonian networks are analysed using numerical runoff simulation. For Hortonian networks the variability of the geometry of individual channels and subcatchments within each Strahler order has generally little effect upon the overall character of the hydrograph in channels of higher order. If the network is also structurally Hortonian, the analysis of the simultaneous formation, travel, and concentration of the hydrographs in all channels of the network can be simplified to a sequence of one representative hydrograph per channel order. This approach is used in this study. Three major runoff processes control the flow hydrograph characteristics: the overland flow process which determines the water supply to the drainage network; the channel flow process which translates the hydrograph in space and time; and the drainage network process which concentrates and magnifies the flow at the junctions of the drainage network. Functional relations for the hydrograph peak, timing, and flow velocity are presented. For a given uniform rainfall and infiltration rate, the peak of the channel flow hydrograph is shown to increase geometrically with channel order, and its magnitude is directly related to the bifurcation ratio. The travel time of the peak also increases geometrically with channel order, and it is directly related to the channel length ratio over velocity ratio. The flow velocity of the peak changes in a downstream direction as a function of the bifurcation and slope ratio. It was also found that for negligible channel storage the channel flow and drainage network processes do not contribute significantly to the observed nonlinear response of a watershed to precipitation.  相似文献   

8.
针对我国在西南构造活跃区修建的水库蓄水与地震活动的关系,本文对蓄水已长达7年并在高水位运维多个周期的云南小湾水库,采用结合波形互相关技术的双差地震定位法对水库库区及周边地区2005年7月至2014年12月发生的M≥1.0级地震进行了精定位处理,结果显示出明显的地震成丛活动特征,库区内外的地震震源深度差别较大.对地震震源深度、地震活动与水库蓄水水位及b值分析结果表明:小湾水库蓄水后地震活动明显增多,有水库触发地震发生,触发地震主要分布在沿黑惠江(A)和澜沧江流域(B、C)的3组地震丛中,且3个区的触发地震类型均为快速响应型;在水库蓄水响应活动最明显的地震丛集区A,展现出明显的随水库蓄水水体渗透发生地震"迁移"活动的现象;但库区内也存在着与蓄水关系不大的可能属于正常构造地震的活动,而库区外的地震活动与水库蓄水没有什么相关性,很可能是属于正常的构造地震.综合断层展布、岩性分布及震源深度分析,认为水库蓄水引起的溶岩作用和渗透作用及断层活动可能是小湾水库触发地震的主控因素.  相似文献   

9.
A geomorphological instantaneous unit hydrograph (GIUH) rainfall‐runoff model was applied in a 31 km2 montane catchment in Scotland. Modelling was based on flow path length distributions derived from a digital terrain model (DTM). The model was applied in two ways; a single landscape unit response based on the DTM alone, and a two‐landscape unit response, which incorporated the distribution of saturated areas derived from field‐validated geographic information system (GIS) analysis based on a DTM and soil maps. This was to test the hypothesis that incorporation of process‐information would enhance the model performance. The model was applied with limited multiple event calibration to produce parameter sets which could be applied to a spectrum of events with contrasting characteristics and antecedent conditions. Gran alkalinity was used as a tracer to provide an additional objective measure for assessing model performance. The models captured the hydrological response dynamics of the catchment reasonably well. In general, the single landscape unit approach produced the best individual model performance statistics, though the two‐landscape unit approach provided a range of models, which bracketed the storm hydrograph response more realistically. There was a tendency to over‐predict the rising limb of the hydrograph, underestimate large storm event peaks and anticipate the hydrograph recession too rapidly. Most of these limitations could be explained by the simplistic assumptions embedded within the GIUH approach. The modelling also gave feasible predictions of stream water chemistry, though these could not be used as a basis for model rejection. Nevertheless, the study suggested that the approach has potential for prediction of hydrological response in ungauged montane headwater basins. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

10.
The separation of baseflow is an important issue in hydrology. The objective of this paper is to develop a new baseflow separation method based on the Horton infiltration capacity curve. For this purpose, the mathematical expressions of three parameters for the Horton infiltration capacity curve were derived in terms of rainfall and runoff data, and the lag time of the effective rainfall in the unsaturated zone and the groundwater flow routing equation are also presented. With these equations, the baseflow hydrographs at the outlet of the basin can be separated. The flow chart of the proposed method for baseflow hydrograph separation is given. Three recent baseflow separation methods, i.e. digital filter, non‐linear reservoir and the Boussinesq equation, were chosen as parallel schemes to compare with the proposed method. Rainfall‐runoff data from four watersheds located in different climatic regions in China were selected and used as case studies. Test and application results indicate that the proposed baseflow hydrograph separation method is in accordance with the hydrological physical process. The proposed method is comparable with current available methods and reduces some of the subjective aspects for the rising limb of the baseflow hydrograph, and it is useful for unit hydrograph analysis and for the study of the rainfall‐runoff relationship. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

11.
Based on the measuring data and Digital Elevation Data (DEM) in a typical watershed--Hemingguan Watershed, Nanbu County, Sichuan Province of China, a GIS-based distributed soil erosion model was developed particularly for the purple soil type. It takes 20 m × 20 m grid as calculating unit and operates at 10-minute time interval. The required input data to the model include DEM, soil, land use, and time-series of precipitation and evaporation loss. The model enables one to estimate runoff, erosion and sediment yield for each grid cell and route the flow along its flow path to the watershed outlet. Furthermore, the model is capable of calculating the total runoff; erosion and sediment yield for the entire watershed by recursion algorithm. The validation of the model demonstrated that it could quantitatively simulate the spatial distribution of hydrological variables in a watershed, such as runoff, vegetation entrapment, soil erosion, the degree of soil and water loss. Moreover, it can evaluate the effect of land use change on the runoff generation and soil erosion with an accuracy of 80% and 75% respectively. The application of this model to a neighboring watershed with similar conditions indicates that this distributed model could be extended to other similar regions in China.  相似文献   

12.
In watershed modelling, the traditional practice of arbitrarily filling topographic depressions in digital elevation models has raised concerns. Advanced high‐resolution remote sensing techniques, including airborne scanning laser altimetry, can identify naturally occurring depressions that impact overland flow. In this study, we used an ensemble physical and statistical modelling approach, including a 2D hydraulic model and two‐point connectivity statistics, to quantify the effects of depressions on high‐resolution overland flow patterns across spatial scales and their temporal variations in single storm events. Computations for both models were implemented using graphic processing unit‐accelerated computing. The changes in connectivity statistics for overland flow patterns between airborne scanning laser altimetry‐derived digital elevation models with (original) and without (filled) depressions were used to represent the shifts of overland flow response to depressions. The results show that depressions can either decrease or increase (to a lesser degree and shorter duration) the probability that any two points (grid locations) are hydraulically connected by overland flow pathways. We used macro‐connectivity states (Φ) as a watershed‐specific indicator to describe the spatiotemporal thresholds of connectivity variability caused by depressions. Four states of Φ are identified in a studied watershed, and each state represents different magnitudes of connectivity and connectivity changes (caused by depressions). The magnitude of connectivity variability corresponds to the states of Φ, which depend on the topological relationship between depressions, the rising/recession limb, and the total rainfall amount in a storm event. In addition, spatial distributions of connectivity variability correlate with the density of depression locations and their physical structures, which cause changes in streamflow discharge magnitude. Therefore, this study suggests that depressions are “nontrivial” in watershed modelling, and their impacts on overland flow should not be neglected. Connectivity statistics at different spatial scales and time points within a watershed provide new insights for characterizing the distributed and accumulated effects of depressions on overland flow.  相似文献   

13.
Shin‐Jen Cheng 《水文研究》2010,24(20):2851-2870
This study explores the hydrograph characteristics of quick and slow runoffs in watershed outlet hydrographs. The quick and slow runoffs were modelled using a conceptual model of three linear cascade reservoirs that have exponential decay response expressions. Mean rainfall for model inputs was estimated using the block Kriging method. The 107 storms during the 1966–2008 events were classified as large, medium and small events according to the observed streamflow. The optimal hydrograph parameters for 61 rainfall‐runoff events were calibrated using the shuffled complex evolution optimal algorithm. The efficacy of the model was verified using the seven averaged parameters of three types of 46 events and was compared with three evaluation criteria resulting from the Nash model. The 61 calibrations were used to analyse and compare the characteristics of quick and slow flows in early and later periods (1966–2002 and 2003–2008). Finally, the following five conclusions were obtained: (1) The base time of a slow runoff hydrograph is the same as that of a total runoff hydrograph. (2) A quick runoff with a long period occurs when soil antecedent moisture is low and with a short period under a high value. (3) The time to peak of hydrograph components is directly proportional to peak time of a hyetograph; the time to peak of quick and slow flows is about 0·97 and 1·12 times the peak time of a hyetograph, respectively. (4) The peak of hydrograph components is relative to a total runoff hydrograph; the percentages for quick runoff are approximately 71% and 13% for slow flow. Finally, (5) the volume of a quick runoff component is 49% of a total runoff volume and 37% for a slow runoff volume. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

14.
For snowmelt-driven flood studies, snow water equivalent (SWE) is frequently estimated using snow depth data. Accurate measurements of snow depth are important in providing data for continuous hydrologic simulations of such watersheds. A new hydrologic fidelity metric is proposed in this study to evaluate the potential contribution of particular snow depth datasets to flow characteristics using observed data and hydrologic modeling using the Variable Infiltration Capacity (VIC) model. Data-based hydrologic fidelity of snow depth measurements is defined as a categorical skill score between the snow depth in the watershed and the hydrograph peak or volume at the watershed outlet. Similarly, model-based hydrologic fidelity is defined as a categorical skill score between the model-simulated snow depth and the model-simulated hydrograph peak or volume. The proposed framework is illustrated using the Pecatonica River watershed in the USA, indicating which sites have a higher hydrologic fidelity, which is preferred in hydrologic studies.  相似文献   

15.
This study was designed to develop a physically based hydrological model to describe the hydrological processes within forested mountainous river basins. The model describes the relationships between hydrological fluxes and catchment characteristics that are influenced by topography and land cover. Hydrological processes representative of temperate basins in steep terrain that are incorporated in the model include intercepted rainfall, evaporation, transpiration, infiltration into macropores, partitioning between preferential flow and soil matrix flow, percolation, capillary rise, surface flow (saturation‐excess and return flow), subsurface flow (preferential subsurface flow and baseflow) and spatial water‐table dynamics. The soil–vegetation–atmosphere transfer scheme used was the single‐layer Penman–Monteith model, although a two‐layer model was also provided. The catchment characteristics include topography (elevation, topographic indices), slope and contributing area, where a digital elevation model provided flow direction on the steepest gradient flow path. The hydrological fluxes and catchment characteristics are modelled based on the variable source‐area concept, which defines the dynamics of the watershed response. Flow generated on land for each sub‐basin is routed to the river channel by a kinematic wave model. In the river channel, the combined flows from sub‐basins are routed by the Muskingum–Cunge model to the river outlet; these comprise inputs to the river downstream. The model was applied to the Hikimi river basin in Japan. Spatial decadal values of the normalized difference vegetation index and leaf area index were used for the yearly simulations. Results were satisfactory, as indicated by model efficiency criteria, and analysis showed that the rainfall input is not representative of the orographic lifting induced rainfall in the mountainous Hikimi river basin. Also, a simple representation of the effects of preferential flow within the soil matrix flow has a slight significance for soil moisture status, but is insignificant for river flow estimations. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

16.
C. Fleurant  B. Kartiwa  B. Roland 《水文研究》2006,20(18):3879-3895
The rainfall‐runoff modelling of a river basin can be divided into two processes: the production function and the transfer function. The production function determines the proportion of gross rainfall actually involved in the runoff. The transfer function spreads the net rainfall over time and space in the river basin. Such a transfer function can be modelled using the approach of the geomorphological instantaneous unit hydrograph (GIUH). The effectiveness of geomorphological models is actually revealed in rainfall‐runoff modelling, where hydrologic data are desperately lacking, just as in ungauged basins. These models make it possible to forecast the hydrograph shape and runoff variation versus time at the basin outlet. This article is an introduction to a new GIUH model that proves to be simple and analytical. Its geomorphological parameters are easily available on a map or from a digital elevation model. This model is based on general hypotheses on symmetry that provide it with multiscale versatile characteristics. After having validated the model in river basins of very different nature and size, we present an application of this model for rainfall‐runoff modelling. Since parameters are determined relying on real geomorphological data, no calibration is necessary, and it is then possible to carry out rainfall‐runoff simulations in ungauged river basins. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

17.
The influence of the method of identification of the drainage network on its geomorphological characteristics and on its hydrological response is analysed. Blue lines, photo-interpreted networks and networks generated from digital elevation models (DEMs) by an automatic algorithm are compared with field observations for two small alpine catchments. The comparisons are carried out in quantitative terms by using several geomorphological indices and functions and by calculating the hydrological response of the networks as represented by their geomorphologic instantaneous unit hydrograph (GIUH). The results show that the effect of the identification method on the geomorphological indices and on the hydrological response is significant, and that the threshold area for channel initiation is not constant. Moreover, the available data show a poor correlation between local slope and threshold area. Finally, the influence of the threshold area on the shape of the GIUH is larger when the residence time on the hillslopes is of the same order as the residence time in the network. In the opposite case, the variability of the flow velocity along the network seems to play an important role. © 1997 John Wiley & Sons, Ltd.  相似文献   

18.
Abstract

The snow and rain in the Himalayas are the main sources of supply for the rivers in the Indo-Gangetic plains. Irrigation, hydropower generation, and water supply are very much dependent on the availability of water in the Himalaya rivers. Mathematical models serve as important aids for the estimation of water availability in rivers. In the present study the SLURP watershed model is applied to a rainfed area of the Satluj catchment located in the western Himalayas, India. The SLURP model developed at NHRI, Canada, is a distributed conceptual model which simulates the behaviour of a watershed by carrying out vertical water balances for each element of a matrix of landcovers and subareas of a watershed and then routing the resulting runoff between subareas. The ILWIS geographic information system was used to prepare the input data required for SLURP and land use data were obtained from the IRS satellite LISS II visible and near infrared sensors. The simulated flows at the Bhakhra Dam outlet of the Satluj catchment were computed and found to compare well with the observed flows.  相似文献   

19.
The design storm approach, where the subject criterion variable is evaluated by using a synthetic storm pattern composed of identical return frequencies of storm pattern input, is shown to be an effective approximation to a considerably more complex probabilistic model. The single area unit hydrograph technique is shown to be an accurate mathematical model of a highly discretized catchment with linear routing for channel flow approximation, and effective rainfalls in subareas which are linear with respect to effective rainfall output for a selected “loss” function. The use of a simple “loss” function which directly equates to the distribution of rainfall depth-duration statistics (such as a constant fraction of rainfall, or a ?-index model) is shown to allow the pooling of data and thereby provide a higher level of statistical significance (in estimating T-year outputs for a hydrologic criterion variable) than use of an arbitrary “loss” function. The above design storm unit hydrograph approach is shown to provide the T-year estimate of a criterion variable when using rainfall data to estimate runoff.  相似文献   

20.
Book Review     
Abstract

The instantaneous unit hydrograph (IUH) of a watershed is the result of one instantaneous unit of rainfall excess distributed uniformly over the watershed. Although the geomorphological characteristics of the basin remain relatively constant, the variable characteristics of storms cause variations in the shape of the resulting hydrographs. It is, therefore, inadequate to use one typical IUH to represent the hydrological response generated from any specific storm. In this study, a variable IUH was derived that directly reflects the time-varying rainfall intensity during storms. The rainfall intensity used to generate the variable IUH at time t is the mean rainfall intensity occurring from the time t—T c to t in which T c is the watershed time of concentration. Hydrological records from three watersheds in Taiwan were used to demonstrate the applicability of the proposed model. The results show that better simulations can be obtained by using the proposed model than by using the conventional unit hydrograph method, especially for concentrated rainstorm cases.  相似文献   

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

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