首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 453 毫秒
1.
A method is presented to evaluate the storm runoff contributions from different land‐use class areas within a river basin using the geographical information system‐based hydrological model WetSpa. The modelling is based on division of the catchment into a grid mesh. Each cell has a unique response function independent of the functioning of other cells. Summation of the flow responses from the cells with the same land‐use type results in the storm runoff contribution from these areas. The model was applied on the Steinsel catchment in the Alzette river basin, Grand Duchy of Luxembourg, with 52 months of meteo‐hydrological measurements. The simulation results show that the direct runoff from urban areas is dominant for a flood event compared with runoff from other land‐use areas in this catchment, and this tends to increase for small floods and for the dry‐season floods, whereas the interflow from forested, pasture and agricultural field areas contributes to recession flow. It is demonstrated that the relative contribution from urban areas decreases with flow coefficient, that cropland relative contribution is nearly constant, and that the relative contribution from grassland and woodland increases with flow coefficient with regard to their percentage of land‐use class areas within the study catchment. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

2.
We began a study, in 1996, to compare ground water quality under irrigated and nonirrigated agriculture, sewered and nonsewered residential developments, industrial, and nondeveloped land uses. Twenty-three monitoring wells were completed in the upper meter of an unconfined sand aquifer. Between 1997 and 2000, sampling occurred quarterly for major ions, trace inorganic chemicals, volatile organic compounds (VOCs), herbicides, and herbicide degradates. On single occasions, we collected samples for polynuclear aromatic hydrocarbons (PAHs), perchlorate, and coliform bacteria. We observed significant differences in water chemistry beneath different land uses. Concentrations of several trace inorganic chemicals were greatest under sewered urban areas. VOC detection frequencies were 100% in commercial areas, 52% in sewered residential areas, and <10% for other land uses. Median nitrate concentrations were greatest under irrigated agriculture (15,350 microg/L) and nonsewered residential areas (6080 microg/L). Herbicides and degradates of acetanilide and triazine herbicides were detected in 86% of samples from irrigated agricultural areas, 68% of samples from nonirrigated areas, and <10% of samples from other land uses. Degradates accounted for 96% of the reported herbicide mass. We did not observe seasonal differences in water chemistry, but observed trends in water chemistry when land use changes occurred. Our results show land use is the dominant factor affecting shallow ground water quality. Trend monitoring programs should focus on areas where land use is changing, while resource managers and planners must consider potential impacts of land use changes on ground water quality.  相似文献   

3.
The impact of Hurricane Ivan on water quality in Pensacola Bay was investigated by MODIS 250 m remote sensing of chlorophyll-a concentrations at different time slots before and after the hurricane event. Before the hurricane, the mean chlorophyll-a in the Bay was 5.3 μg/L. Heavy rainfall occurred during the hurricane landfall. The 48 h rainfall reached 40 cm and the peak storm surge reached 3 m on 9/16. After the rainstorm and during the storm surge on 9/17/2004, the mean chlorophyll-a concentration substantially increased to 14.7 μg/L. 26.3% water area was in the poor-water-quality condition (chl-a > 20 μg/L). This indicates that heavy nutrient loads from urban stormwater runoff and storm-surge inundation simulated chlorophyll bloom. After the end of the storm surge on 9/18/2004, the mean chlorophyll dropped to 2.0 μg/L, suggesting the effective flushing of polluted water from the bay to the Gulf of Mexico by the storm-surge. The good water quality condition lasted for almost several weeks after the storm surge. The peak river flow, arriving on the 4th day after the peak storm surge, did not alter the good water quality situation in the bay. This indicates that urban stormwater runoff rather than the river inflow is the major pollutant source for water quality in Pensacola Bay during the hurricane.  相似文献   

4.
The water quality of urban drainage ditches in lowlands in the Rhine‐Meuse delta was analysed with principal component analysis (PCA) during a dry period and a rain storm, and related to the seepage of polluted river water and effective impervious area (EIA). This was done in order to test the hypothesis that seepage of river water and storm water runoff from impervious areas strongly determine the water quality of urban drainage systems along large lowland rivers. Our analysis revealed that upward seepage of groundwater originating from rivers Rhine and Meuse was positively correlated with nitrate, potassium, sodium and chloride and negatively correlated with alkalinity, calcium, magnesium and iron. EIA was correlated with very few environmental variables (i.e. phosphate, pH and iron in the dry period and iron during the rain storm). Nickel and zinc concentrations generally exceeded the maximum allowable concentrations (MAC), while lead and phosphorus concentrations were just above the nutrient standards and MAC in a few locations during the rain storm. To optimize water quality in urban water systems, attention should be paid to all sources of pollution and not only to EIA. The impact of local groundwater seepage originating from large rivers in lowlands on the chemistry of urban water systems is often underestimated and should be taken into account when assessing water quality and improving water quality status. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

5.
Effective impervious area for runoff in urban watersheds   总被引：2，自引：0，他引：2       下载免费PDF全文
Effective impervious area (EIA), or the portion of total impervious area (TIA) that is hydraulically connected to the storm sewer system, is an important parameter in determining actual urban runoff. EIA has implications in watershed hydrology, water quality, environment, and ecosystem services. The overall goal of this study is to evaluate the application of successive weighted least square (WLS) method to urban catchments with different sizes and various hydrologic conditions to determine EIA fraction. Other objectives are to develop insights on the data selection issues, EIA fraction, EIA/TIA ratio, and runoff source area patterns in urban catchments. The successive WLS method is applied to 50 urban catchments with different sizes from less than 1 ha to more than 2000 ha in Minnesota, Wisconsin, Texas, USA as well as Europe, Canada, and Australia. The average, median, and standard deviation of EIA fractions for the 42 catchments with residential land uses are found to be 0.222, 0.200, and 0.113, respectively. These values for the EIA/TIA ratio in the same 42 catchments are 0.50, 0.48, and 0.21, respectively. While the EIA/TIA results indicate the importance of EIA, 95% prediction interval of the mean EIA/TIA is found to be 0.07 to 0.93, which shows that using an average value for this ratio in each land use to determine EIA from TIA in ungauged urban watersheds can be misleading. The successive WLS method was robust and is recommended for determining EIA in gauged urban catchments. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

6.
Research on the water quality of urban runoff has so far focussed on the post‐development phase of urban catchments, whereas water quality in developing areas under construction has remained less understood. The construction phase, however, may constitute a considerable source of diffuse pollution in urban areas. This study investigated the mechanisms affecting water quality in residential areas during the construction and post‐construction phases. Water quality (suspended solids, total phosphorus, total nitrogen and chemical oxygen demand) and runoff were monitored over a period of 5 years in three catchments located in the city of Espoo in Southern Finland. The catchments included an urbanizing area under construction and low‐density and high‐density residential areas. The water quality was quantified in terms of event mean concentrations and loads. The key influential variables explaining the water quality in a multiple linear regression analysis included hydrological variables (event volume and intensity), antecedent conditions and a variable describing the ongoing construction projects. Construction activities in the developing catchment had a profound impact on water quality. Inclusion of the variables describing activities, such as earth moving works, paving, house construction and temporary wastewater discharges, was necessary to explain water quality variations in the developing catchment. The importance of antecedent conditions as an explanatory variable depended on the site and the pollutant in question. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

7.
Hydrological modelling of mesoscale catchments is often adversely affected by a lack of adequate information about specific site conditions. In particular, digital land cover data are available from data sets which were acquired on a European or a national scale. These data sets do not only exhibit a restricted spatial resolution but also a differentiation of crops and impervious areas which is not appropriate to the needs of mesoscale hydrological models. In this paper, the impact of remote sensing data on the reliability of a water balance model is investigated and compared to model results determined on the basis of CORINE (Coordination of Information on the Environment) Land Cover as a reference. The aim is to quantify the improved model performance achieved by an enhanced land cover representation and corresponding model modifications. Making use of medium resolution satellite imagery from SPOT, LANDSAT ETM+ and ASTER, detailed information on land cover, especially agricultural crops and impervious surfaces, was extracted over a 5-year period (2000–2004). Crop-specific evapotranspiration coefficients were derived by using remote sensing data to replace grass reference evapotranspiration necessitated by the use of CORINE land cover for rural areas. For regions classified as settlement or industrial areas, degrees of imperviousness were derived. The data were incorporated into the hydrological model GROWA (large-scale water balance model), which uses an empirical approach combining distributed meteorological data with distributed site parameters to calculate the annual runoff components. Using satellite imagery in combination with runoff data from gauging stations for the years 2000–2004, the actual evapotranspiration calculation in GROWA was methodologically extended by including empirical crop coefficients for actual evapotranspiration calculations. While GROWA originally treated agricultural areas as homogeneous, now a consideration and differentiation of the main crops is possible. The accuracy was determined by runoff measurements from gauging stations. Differences in water balances resulting from the use of remote sensing data as opposed to CORINE were analysed in this study using a representative subcatchment. Resulting Nash–Sutcliff model efficiencies improved from 0.372 to 0.775 and indicate that the enhanced model can produce thematically more accurate and spatially more detailed local water balances. However, the proposed model enhancements by satellite imagery have not exhausted the full potential of water balance modelling, for which a higher temporal resolution is required.  相似文献   

8.
For the appropriate management of water resources in a watershed, it is essential to calculate the time distribution of runoff for the given rainfall event. In this paper, a kinematic‐wave‐based distributed watershed model using finite element method (FEM), geographical information systems (GIS) and remote‐sensing‐based approach is presented for the runoff simulation of small watersheds. The kinematic wave equations are solved using FEM for overland and channel flow to generate runoff at the outlet of the watershed concerned. The interception loss is calculated by an empirical model based on leaf area index (LAI). The Green‐Ampt Mein Larson (GAML) model is used for the estimation of infiltration. Remotely sensed data has been used to extract land use (LU)/land cover (LC). GIS have been used to prepare finite element grid and input files such as Manning's roughness and slope. The developed overland flow model has been checked with an analytical solution for a hypothetical watershed. The model has been applied to a gauged watershed and an ungauged watershed. From the results, it is seen that the model is able to simulate the hydrographs reasonably well. A sensitivity analysis of the model is carried out with the calibrated infiltration parameters, overland flow Manning's roughness, channel flow Manning's roughness, time step and grid size. The present model is useful in predicting the hydrograph in small, ungauged watersheds. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

9.
A three‐dimensional numerical modelling system is developed to study transformation processes of water resources in alluvial fan and river basin along the middle reaches of the Heihe River Basin, Northwest China, an arid and semi‐arid region. Integrating land utilization, remote sensing and geographic information systems, we have developed a numerical modelling system that can be used to quantify the effects of land use and anthropogenic activities on the groundwater system as well as to investigate the interaction between surface water and groundwater. Various hydraulic measurements are used to identify and calibrate the hydraulic boundary conditions and spatial distributions of hydraulic parameters. In the modelling study, various water exchanges and human effects on the watershed system are considered. These include water exchange between surface water and groundwater, groundwater pumping, lateral water recharges from mountain areas, land utilization, and infiltration and evaporation in the irrigation and non‐irrigation areas. The modelling system provides a quantitative method to describe spatial and temporal distributions and transformations between various water resources, and it has application to other watersheds in arid and semi‐arid areas. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

10.
Large urban areas are typically characterized by a mosaic of different land uses, with contrasting mixes of impermeable and permeable surfaces that alter “green” and “blue” water flux partitioning. Understanding water partitioning in such heterogeneous environments is challenging but crucial for maintaining a sustainable water management during future challenges of increasing urbanization and climate warming. Stable isotopes in water have outstanding potential to trace the partitioning of rainfall along different flow paths and identify surface water sources. While isotope studies are an established method in many experimental catchments, surprisingly few studies have been conducted in urban environments. Here, we performed synoptic sampling of isotopes in precipitation, surface water and groundwater across the complex city landscape of Berlin, Germany, for a large -scale overview of the spatio-temporal dynamics of urban water cycling. By integrating stable isotopes of water with other hydrogeochemical tracers we were able to identify contributions of groundwater, surface runoff during storm events and effluent discharge on streams with variable degrees of urbanization. We could also assess the influence of summer evaporation on the larger Spree and Havel rivers and local wetlands during the exceptionally warm and dry summers of 2018 and 2019. Our results demonstrate that using stable isotopes and hydrogeochemical data in urban areas has great potential to improve our understanding of water partitioning in complex, anthropogenically-affected landscapes. This can help to address research priorities needed to tackle future challenges in cities, including the deterioration of water quality and increasing water scarcity driven by climate warming, by improving the understanding of time-variant rainfall-runoff behaviour of urban streams, incorporating field data into ecohydrological models, and better quantifying urban evapotranspiration and groundwater recharge.  相似文献   

11.
The urban environment modifies the hydrologic cycle resulting in increased runoff rates, volumes, and peak flows. Green infrastructure, which uses best management practices (BMPs), is a natural system approach used to mitigate the impacts of urbanization onto stormwater runoff. Patterns of stormwater runoff from urban environments are complex, and it is unclear how efficiently green infrastructure will improve the urban water cycle. These challenges arise from issues of scale, the merits of BMPs depend on changes to small‐scale hydrologic processes aggregated up from the neighborhood to the urban watershed. Here, we use a hyper‐resolution (1 m), physically based hydrologic model of the urban hydrologic cycle with explicit inclusion of the built environment. This model represents the changes to hydrology at the BMP scale (~1 m) and represents each individual BMP explicitly to represent response over the urban watershed. Our study varies both the percentage of BMP emplacement and their spatial location for storm events of increasing intensity in an urban watershed. We develop a metric of effectiveness that indicates a nonlinear relationship that is seen between percent BMP emplacement and storm intensity. Results indicate that BMP effectiveness varies with spatial location and that type and emplacement within the urban watershed may be more important than overall percent.  相似文献   

12.
In the face of increasing urbanization, there is growing interest in application of microscale hydrologic solutions to minimize storm runoff and conserve water at the source. In this study, a physically based numerical model was developed to understand hydrologic processes better at the urban residential scale and the interaction of these processes among different best management practices (BMPs). This model simulates hydrologic processes using an hourly interval for over a full year or for specific storm events. The model was applied to treatment and control single‐family residential parcels in Los Angeles, California. Data collected from the control and treatment sites over 2 years were used to calibrate and validate the model. Annual storm runoff to the street was eliminated by 97% with installation of rain gutters, a driveway interceptor, and lawn retention basin. Evaluated individually, the driveway interceptor was the most effective BMP for storm runoff reduction (65%), followed by the rain gutter installation (28%), and lawn converted to retention basin (12%). An 11 m3 cistern did not substantially reduce runoff, but provided 9% of annual landscape irrigation demand. Simulated landscape irrigation water use was reduced 53% by increasing irrigation system efficiency, and adjusting application rates monthly based on plant water demand. The model showed that infiltration and surface runoff processes were particularly sensitive to the soil's physical properties and its effective depth. Replacing the existing loam soil with clay soil increased annual runoff discharge to the street by 63% when climate and landscape features remained unchanged. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

13.
  相似文献   

14.
In the modeling of urban storm water runoff temperatures, the contribution of rooftops to the heating of rainfall runoff is usually neglected or not mentioned in the literature. In this paper we examine the accuracy of this assumption (a) by analyzing temperature data that we recorded on a residential rooftop, a commercial rooftop, and a concrete driveway, and (b) by simulating temperature profiles within rooftops and pavements, and estimating heat transfer amounts from these surfaces to rainfall runoff (‘heat export’). Analysis of both wet‐ and dry‐weather temperature data recorded in the Minneapolis/St. Paul metropolitan area (north‐central USA) over periods of several months leads to the conclusion that (a) a concrete driveway has a far greater capacity for heat storage and release than a shingled residential rooftop, and (b) an insulated commercial rooftop is able to store and release more heat than the residential rooftop. Unexpectedly, the rainfall events with the highest dew point (rainfall) and surface temperatures often occurred during late night or early morning hours, and not during daylight hours. The analysis of three rainfall events showed that the heat export from the commercial rooftop was roughly three times that of the residential rooftop, but only 30%–90% of the heat export from a concrete driveway. Maximum (potential) heat export was significantly higher for the driveway than for either rooftop. In conclusion, the results of the data analysis and heat export simulations support the assumption that residential rooftops contribute very little heating to runoff from rainfall, while commercial rooftops may have a thermal impact on rainfall runoff because of their greater thermal storage capacity. Commercial rooftops, in addition to asphalt and concrete pavements, should be considered when estimating water temperature of rainfall runoff from urbanized areas and the associated impact on the thermal regime of streams and fish habitat. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

15.
The quality of stormwater runoff from seaports can be an important source of pollution to the marine environment. Currently, little knowledge exists with regards to the pollutant generation capacity specific to seaports as they do not necessarily compare well with conventional urban land use. The research project focussed on the assessment of pollutant build-up and wash-off. The study was undertaken using rainfall simulation and small impervious plots for different port land uses with the results obtained compared to typical urban land uses.The study outcomes confirmed that the Port land uses exhibit comparatively lower pollutant concentrations. However, the pollutant characteristics varied across different land uses. Hence, the provision of stereotypical water quality improvement measures could be of limited value. Particle size <150 μm was predominant in suspended solids. Therefore, if suspended solids are targeted as the surrogate parameter for water quality improvement, this particle size range needs to be removed.  相似文献   

16.
Experimental research in the Ethiopian highlands found that saturation excess induced runoff and erosion are common in the sub‐humid conditions. Because most erosion simulation models applied in the highlands are based on infiltration excess, we, as an alternative, developed the Parameter Efficient Distributed (PED) model, which can simulate water and sediment fluxes in landscapes with saturation excess runoff. The PED model has previously only been tested at the outlet of a watershed and not for distributed runoff and sediment concentration within the watershed. In this study, we compare the distributed storm runoff and sediment concentration of the PED model against collected data in the 95‐ha Debre Mawi watershed and three of its nested sub‐watersheds for the 2010 and 2011 rainy seasons. In the PED model framework, the hydrology of the watershed is divided between infiltrating and runoff zones, with erosion only taking place from two surface runoff zones. Daily storm runoff and sediment concentration values, ranging from 0.5 to over 30 mm and from 0.1 to 35 g l?1, respectively, were well simulated. The Nash Sutcliffe efficiency values for the daily storm runoff for outlet and sub‐watersheds ranged from 0.66 to 0.82, and the Nash–Sutcliffe efficiency for daily sediment concentrations were greater than 0.78. Furthermore, the model uses realistic fractional areas for surface and subsurface flow contributions, for example between saturated areas (15%), degraded areas (30%) and permeable areas (55%) at the main outlet, while close similarity was found for the remaining hydrology and erosion parameter values. One exception occurred for the distinctly greater transport limited parameter at the actively gullying lower part of the watershed. The results suggest that the model based on saturation excess provides a good representation of the observed spatially distributed runoff and sediment concentrations within a watershed by modelling the bottom lands (as opposed to the uplands) as the dominant contributor of the runoff and sediment load. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

17.
Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.  相似文献   

18.
During the last two decades, remote sensing data have led to tremendous progress in advancing flood inundation modelling. In particular, low‐cost space‐borne data can be invaluable for large‐scale flood studies in data‐scarce areas. Various satellite products yield valuable information such as land surface elevation, flood extent and water level, which could potentially contribute to various flood studies. An increasing number of research studies have been dedicated to exploring those low‐cost data towards building, calibration and evaluation, and remote‐sensed information assimilation into hydraulic models. This paper aims at reviewing these recent scientific efforts on the integration of low‐cost space‐borne remote sensing data with flood modelling. Potentials and limitations of those data in flood modelling are discussed. This paper also introduces the future satellite missions and anticipates their likely impacts in flood modelling. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

19.
Heejun Chang 《水文研究》2007,21(2):211-222
This study investigates changes in streamflow characteristics for urbanizing watersheds in the Portland Metropolitan Area of Oregon for the period from 1951 to 2000. The objective of this study was to assess how mean annual runoff ratio, mean seasonal runoff ratio, annual peak runoff ratio, changes in streamflow in response to storm amount, the fraction of time that the daily mean flow exceeds the annual mean flow, 3‐day recession constants, and dry/wet flow ratio vary among watersheds with different degrees of urban development. There were no statistically significant changes in annual runoff ratio and annual peak runoff ratio for the mixed land‐use watershed (Tualatin River watershed) and the urban watershed (Johnson Creek watershed) during the entire study period. The Tualatin River watershed, where most of the urban development occurred in a lower part of the watershed, showed a statistically significant increase in annual peak runoff ratio during the 1976 and 2000 period. The Upper Tualatin River watershed illustrated a significant decrease in annual peak runoff ratio for the entire study period. With significant differences in seasonal runoff ratio, only Johnson Creek exhibited a significant increase in both wet and dry season runoff ratios. Streamflow during storm events declined rapidly in the urban watershed, with a high 3‐day recession constant. At an event storm scale, streamflow in Fanno Creek, which is the most urbanized watershed, responded quickly to precipitation input. The fraction of time that the daily mean flow exceeded the annual mean flow and dry/wet flow ratio are all lower in Johnson Creek. This suggests a shorter duration of storm runoff and lower baseflow in the urbanized watershed when compared to the mixed land use watershed. The findings of this study demonstrate the importance of spatial and temporal scale, climate variability, and basin physiographic characteristics in detecting the hydrologic effects of urbanization in the Pacific Northwest of the USA. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

20.
Urbanization strongly changes natural catchment by increasing impervious coverage and by creating a need for efficient drainage systems. Such land cover changes lead to more rapid hydrological response to storms and change distribution of peak and low flows. This study aims to explore and assess how gradual hydrological changes occur during urban development from rural area to a medium‐density residential catchment. The Stormwater Management Model (SWMM) is utilized to simulate a series of scenarios in a same developing urban catchment. Sub‐hourly hydro‐meteorological data in warm season is used to calibrate and validate the model in the fully developed catchment in 2006. The validated model is then applied to other cases in development stage and runoff management scenarios. Based on the simulations and observations, three key problems are solved: (1) how catchment hydrology changes with land cover change, (2) how urban development changes pre‐development flows, and (3) how stormwater management techniques affect catchment hydrology. The results show that the low‐frequency flow rates had remarkably increased from 2004 to 2006 along with the increase of impervious areas. Urbanization in the residential catchment expands the runoff contributing area, accelerates hydrological response, raises peak flows in an order of magnitude of over 10, and more than doubles the total runoff volume. The effects of several LID controls on runoff hydrograph were simulated, and the techniques were able to reduce flows towards the pre‐development levels. However, the partly restored flow regime was still clearly changed in comparison to the pre‐development flow conditions. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

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

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