首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
ABSTRACT: Precipitation and streamflow data from three nested subwatersheds within the Little Washita River Experimental Watershed (LWREW) in southwestern Oklahoma were used to evaluate the capabilities of the Soil and Water Assessment Tool (SWAT) to predict streamflow under varying climatic conditions. Eight years of precipitation and streamflow data were used to calibrate parameters in the model, and 15 years of data were used for model validation. SWAT was calibrated on the smallest and largest sub‐watersheds for a wetter than average period of record. The model was then validated on a third subwatershed for a range in climatic conditions that included dry, average, and wet periods. Calibration of the model involved a multistep approach. A preliminary calibration was conducted to estimate model parameters so that measured versus simulated yearly and monthly runoff were in agreement for the respective calibration periods. Model parameters were then fine tuned based on a visual inspection of daily hydrographs and flow frequency curves. Calibration on a daily basis resulted in higher baseflows and lower peak runoff rates than were obtained in the preliminary calibration. Test results show that once the model was calibrated for wet climatic conditions, it did a good job in predicting streamflow responses over wet, average, and dry climatic conditions selected for model validation. Monthly coefficients of efficiencies were 0.65, 0.86, and 0.45 for the dry, average, and wet validation periods, respectively. Results of this investigation indicate that once calibrated, SWAT is capable of providing adequate simulations for hydrologic investigations related to the impact of climate variations on water resources of the LWREW.  相似文献   

2.
Forty‐five flood control reservoirs, authorized in the Watershed Protection and Flood Prevention Act 1954, were installed by United States Department of Agriculture (USDA) between 1969 and 1982 in the Little Washita River Experimental Watershed (LWREW), located in central Oklahoma. Over time, these reservoirs have lost sediment and flood storage capacity due to sedimentation, with rates dependent on upstream land use and climate variability. In this study, sedimentation rates for 12 reservoirs representing three major land use categories within LWREW were measured based on bathymetric surveys that used acoustic profiling system. Physiographic and climate attributes of drainage area of surveyed reservoirs were extracted from publicly available data sources including topographic maps, digital elevation models, USDA Natural Resource Conservation Service soils, and weather station databases. Correlation, principal component analysis, and stepwise regression were utilized to analyze the relationship between normalized reservoir sedimentation rates (ReSRa) and the drainage area characteristics to determine the major variables controlling sedimentation within the LWREW. Percent of drainage area with extreme slopes, saturated hydraulic conductivity, and maximum daily rainfall event recorded in spring explained most of the variability in ReSRa. It was also found that percent reduction in reservoir surface area can be used as a surrogate for estimating ReSRa. The implications of the results are discussed.  相似文献   

3.
We have enhanced the ability of a widely used watershed model, Hydrologic Simulation Program — FORTRAN (HSPF), to predict low flows by reconfiguring the algorithm that simulates groundwater discharge. During dry weather periods, flow in most streams consists primarily of base flow, that is, groundwater discharged from underlying aquifers. In this study, HSPF's groundwater storage‐discharge relationship is changed from a linear to a more general nonlinear relationship which takes the form of a power law. The nonlinear algorithm is capable of simulating streamflow recession curves that have been found in some studies to better match observed dry weather hydrographs. The altered version of HSPF is implemented in the Chesapeake Bay Program's Phase 5 Model, an HSPF‐based model that simulates nutrient and sediment loads to the Chesapeake Bay, and is tested in the upper Potomac River basin, a 29,950 km2 drainage area that is part of the Bay watershed. The nonlinear relationship improved median Nash‐Sutcliffe efficiencies for log daily flows at the model's 45 calibration points. Mean absolute percent error on low‐flow days dropped in five major Potomac River tributaries by up to 12 percentage points, and in the Potomac River itself by four percentage points, where low‐flow days were defined as days when observed flows were in the lowest 5th percentile range. Percent bias on low‐flow days improved by eight percentage points in the Potomac River, from ?11 to ?3%.  相似文献   

4.
Abstract: This article describes the development of a calibrated hydrologic model for the Blue River watershed (867 km2) in Summit County, Colorado. This watershed provides drinking water to over a third of Colorado’s population. However, more research on model calibration and development for small mountain watersheds is needed. This work required integration of subsurface and surface hydrology using GIS data, and included aspects unique to mountain watersheds such as snow hydrology, high ground‐water gradients, and large differences in climate between the headwaters and outlet. Given the importance of this particular watershed as a major urban drinking‐water source, the rapid development occurring in small mountain watersheds, and the importance of Rocky Mountain water in the arid and semiarid West, it is useful to describe calibrated watershed modeling efforts in this watershed. The model used was Soil and Water Assessment Tool (SWAT). An accurate model of the hydrologic cycle required incorporation of mountain hydrology‐specific processes. Snowmelt and snow formation parameters, as well as several ground‐water parameters, were the most important calibration factors. Comparison of simulated and observed streamflow hydrographs at two U.S. Geological Survey gaging stations resulted in good fits to average monthly values (0.71 Nash‐Sutcliffe coefficient). With this capability, future assessments of point‐source and nonpoint‐source pollutant transport are possible.  相似文献   

5.
Abstract: Storm‐flow transients (i.e., hydrograph rise and fall dynamics) may represent an important aspect of understanding streamflow dynamics. However, little is known about how temporal resolution of transient data and climate variability may color these potential indicators of hydrologic pattern or condition. Warm‐season stream stage and rainfall were monitored continuously (5 min) during the 2002 water year in eight tributaries of the Little Miami River (Ohio), which drain 17‐58 km2 catchments. Rise rates generated using 5‐min data were different than those generated with mean daily data [calculated with the Indicators of Hydrologic Alteration (IHA) software], though fall rates were similar for fine and coarse temporal data. This result suggests that data with low temporal resolution may not be adequate to fully represent the dynamics of storm rise rates. Conversely, fall rates based on daily stage data (via IHA) were similar to those based on the 5‐min data, and so daily mean data may be appropriate for characterizing fall rates. We next analyzed the possible correlations between rainfall variability and storm‐flow stage dynamics. We derived rise and recession rates from storm stage hydrographs by assuming exponential rise and decay of a runoff peak. We found that raw rise rates (Rraw) were correlated with both the maximum rainfall rate and the time to the centroid of a rain event. We subsequently removed the trend based on these rainfall characteristics, which yielded new representations of rise rates abbreviated as Rrate and Rtcent, respectively, and that had lower variability than the uncorrected (raw) data. Fall rates were found to be independent of rainfall characteristics. Due to the predominant influence of stream hydrology upon aquatic biota and nutrient fluxes, our work suggests that these stage data analysis protocols can refine or otherwise reduce variability in these indices by accounting for relevant factors such as rainfall forcing. These protocols for derivation of transient indices should be tested for their potential to improve correlations between stream hydrology and temporally aligned biotic data and dissolved nutrient fluxes in streams.  相似文献   

6.
One of the major factors contributing to surface water contamination in agricultural areas is the use of pesticides. The Soil and Water Assessment Tool (SWAT) is a hydrologic model capable of simulating the fate and transport of pesticides in an agricultural watershed. The SWAT model was used in this study to estimate stream flow and atrazine (2-chloro-4-(ethylamino)-6-(isopropylamino)-s-triazine) losses to surface water in the Cedar Creek Watershed (CCW) within the St. Joseph River Basin in northeastern Indiana. Model calibration and validation periods consisted of five and two year periods, respectively. The National Agricultural Statistics Survey (NASS) 2001 land cover classification and the Soil Survey Geographic (SSURGO) database were used as model input data layers. Data from the St. Joseph River Watershed Initiative and the Soil and Water Conservation Districts of Allen, Dekalb, and Noble counties were used to represent agricultural practices in the watershed which included the type of crops grown, tillage practices, fertilizer, and pesticide application rates. Model results were evaluated based on efficiency coefficient values, standard statistical measures, and visual inspection of the measured and simulated hydrographs. The Nash and Sutcliffe model efficiency coefficients (E(NS)) for monthly and daily stream flow calibration and validation ranged from 0.51 to 0.66. The E(NS) values for atrazine calibration and validation ranged from 0.43 to 0.59. All E(NS) values were within the range of acceptable model performance standards. The results of this study indicate that the model is an effective tool in capturing the dynamics of stream flow and atrazine concentrations on a large-scale agricultural watershed in the midwestern USA.  相似文献   

7.
Abstract: Both ground rain gauge and remotely sensed precipitation (Next Generation Weather Radar – NEXRAD Stage III) data have been used to support spatially distributed hydrological modeling. This study is unique in that it utilizes and compares the performance of National Weather Service (NWS) rain gauge, NEXRAD Stage III, and Tropical Rainfall Measurement Mission (TRMM) 3B42 (Version 6) data for the hydrological modeling of the Middle Nueces River Watershed in South Texas and Middle Rio Grande Watershed in South Texas and northern Mexico. The hydrologic model chosen for this study is the Soil and Water Assessment Tool (SWAT), which is a comprehensive, physical‐based tool that models watershed hydrology and water quality within stream reaches. Minor adjustments to selected model parameters were applied to make parameter values more realistic based on results from previous studies. In both watersheds, NEXRAD Stage III data yields results with low mass balance error between simulated and actual streamflow (±13%) and high monthly Nash‐Sutcliffe efficiency coefficients (NS > 0.60) for both calibration (July 1, 2003 to December 31, 2006) and validation (2007) periods. In the Middle Rio Grande Watershed NEXRAD Stage III data also yield robust daily results (time averaged over a three‐day period) with NS values of (0.60‐0.88). TRMM 3B42 data generate simulations for the Middle Rio Grande Watershed of variable qualtiy (MBE = +13 to ?16%; NS = 0.38‐0.94; RMSE = 0.07‐0.65), but greatly overestimates streamflow during the calibration period in the Middle Nueces Watershed. During the calibration period use of NWS rain gauge data does not generate acceptable simulations in both watersheds. Significantly, our study is the first to successfully demonstrate the utility of satellite‐estimated precipitation (TRMM 3B42) in supporting hydrologic modeling with SWAT; thereby, potentially extending the realm (between 50°N and 50°S) where remotely sensed precipitation data can support hydrologic modeling outside of regions that have modern, ground‐based radar networks (i.e., much of the third world).  相似文献   

8.
ABSTRACT: The performance of two popular watershed scale simulation models — HSPF and SWAT — were evaluated for simulating the hydrology of the 5,568 km2 Iroquois River watershed in Illinois and Indiana. This large, tile drained agricultural watershed provides distinctly different conditions for model comparison in contrast to previous studies. Both models were calibrated for a nine‐year period (1987 through 1995) and verified using an independent 15‐year period (1972 through 1986) by comparing simulated and observed daily, monthly, and annual streamflow. The characteristics of simulated flows from both models are mostly similar to each other and to observed flows, particularly for the calibration results. SWAT predicts flows slightly better than HSPF for the verification period, with the primary advantage being better simulation of low flows. A noticeable difference in the models' hydrologic simulation relates to the estimation of potential evapotranspiration (PET). Comparatively low PET values provided as input to HSPF from the BASINS 3.0 database may be a factor in HSPF's overestimation of low flows. Another factor affecting baseflow simulation is the presence of tile drains in the watershed. HSPF parameters can be adjusted to indirectly account for the faster subsurface flow associated with tile drains, but there is no specific tile drainage component in HSPF as there is in SWAT. Continued comparative studies such as this, under a variety of hydrologic conditions and watershed scales, provide needed guidance to potential users in model selection and application.  相似文献   

9.
ABSTRACT: A grid based daily hydrologic model for a watershed with paddy fields was developed to predict the stream discharge. ASCII formatted elevation, soil, and land use data supported by the GRASS Geographic Information System are used to generate distributed results such as surface runoff and subsurface flow, soil water content, and evapotranspiration. The model uses a single flow path algorithm and simulates a water balance at each grid element. A linear reservoir assumption was used to predict subsurface runoff components. The model was applied to a 75.6 km2 watershed located in the middle of South Korea, and observed stream flow hydrographs from 1995 and 1996 were compared to model predictions. The stream flow predictions of 1995 and 1996 generally agreed with the observed flow, resulting in a Nash‐Sutcliffe efficiency R2 of 0.60 and 0.62, respectively. The hydraulic conductivity for percolating water through the saturated layer affected baseflow generation. The levee height of the paddy influenced the time and magnitude of the surface runoff, depending on irrigation management. The model will be used for making low flow management decisions by evaluating the role of each land use to stream flow, especially in case of paddy decrease by gradual urbanization of a watershed.  相似文献   

10.
This article presents SWATMOD‐Prep, a graphical user interface that couples a SWAT watershed model with a MODFLOW groundwater flow model. The interface is based on a recently published SWAT‐MODFLOW code that couples the models via mapping schemes. The spatial layout of SWATMOD‐Prep guides the user through the process of importing shape files (sub‐basins, hydrologic response units [HRUs], river network) from an existing SWAT model, creating a grid, performing necessary geo‐processing operations to link the models, writing out SWAT‐MODFLOW files, and running the simulation. The option of creating a new single‐layer MODFLOW model for near‐surface alluvial aquifers is available, with the user prompted to provide groundwater surface elevation (through a digital elevation model), aquifer thickness, and necessary aquifer parameter values. The option of simulating nitrate transport in the aquifer also is available, using the reactive transport model RT3D. The interface is in the public domain. It is programmed in Python, with various software packages used for geo‐processing operations (e.g., selection, intersection of rasters) and inputting/outputting data, and is written for Windows. The use of SWATMOD‐Prep is demonstrated for the Little River Experimental Watershed, Georgia. SWATMOD‐Prep and SWAT‐MODFLOW executables are available with an accompanying user's manual at: http://swat.tamu.edu/software/swat-modflow/ . The user's manual also accompanies this article as Supporting Information.  相似文献   

11.
Abstract: A stochastic, spatially explicit method for assessing the impact of land cover classification error on distributed hydrologic modeling is presented. One‐hundred land cover realizations were created by systematically altering the North American Landscape Characterization land cover data according to the dataset’s misclassification matrix. The matrix indicates the probability of errors of omission in land cover classes and is used to assess the uncertainty in hydrologic runoff simulation resulting from parameter estimation based on land cover. These land cover realizations were used in the GIS‐based Automated Geospatial Watershed Assessment tool in conjunction with topography and soils data to generate input to the physically‐based Kinematic Runoff and Erosion model. Uncertainties in modeled runoff volumes resulting from these land cover realizations were evaluated in the Upper San Pedro River basin for 40 watersheds ranging in size from 10 to 100 km2 under two rainfall events of differing magnitudes and intensities. Simulation results show that model sensitivity to classification error varies directly with respect to watershed scale, inversely to rainfall magnitude and are mitigated or magnified by landscape variability depending on landscape composition.  相似文献   

12.
Historically, many watershed studies have been based on using the streamflow flux, typically from a single gauge at the basin's outlet, to support calibration. In this setting, there is great potential for equifinality of parameters during the optimization process, especially for parameters that are not directly related to streamflow. Therefore, some of the optimal parameter values achieved during the autocalibration process may be physically unrealistic. In recent decades a vast array of data from land surface models and remote sensing platforms can help to constrain hydrologic fluxes such as evapotranspiration (ET). While the spatial resolution of these ancillary datasets varies, the continuous spatial coverage of these gridded datasets provides flux measurements across the entire basin, in stark contrast to point‐based streamflow data. This study uses Global Land Evaporation: the Amsterdam Model data to constrain Soil and Water Assessment Tool parameter values associated with ET to a more physically realistic range. The study area is the Little Washita River Experimental Watershed, in southern Oklahoma. Traditional objective metrics such as the Nash‐Sutcliffe coefficients record no performance improvement after application of this method. However, there is a dramatic increase in the number of days with receding flow where simulations match observed streamflow.  相似文献   

13.
Abstract: The main objective of the study is to examine the accuracy of and differences among simulated streamflows driven by rainfall estimates from a network of 22 rain gauges spread over a 2,170 km2 watershed, NEXRAD Stage III radar data, and Tropical Rainfall Measuring Mission (TRMM) 3B42 satellite data. The Gridded Surface Subsurface Hydrologic Analysis (GSSHA), a physically based, distributed parameter, grid‐structured, hydrologic model, was used to simulate the June‐2002 flooding event in the Upper Guadalupe River watershed in south central Texas. There were significant differences between the rainfall fields estimated by the three types of measurement technologies. These differences resulted in even larger differences in the simulated hydrologic response of the watershed. In general, simulations driven by radar rainfall yielded better results than those driven by satellite or rain‐gauge estimates. This study also presents an overview of effects of land cover changes on runoff and stream discharge. The results demonstrate that, for major rainfall events similar to the 2002 event, the effect of urbanization on the watershed in the past two decades would not have made any significant effect on the hydrologic response. The effect of urbanization on the hydrologic response increases as the size of the rainfall event decreases.  相似文献   

14.
Abstract: The hydrologic performance of DRAINMOD 5.1 was assessed for the southern Quebec region considering freezing/thawing conditions. A tile drained agricultural field in the Pike River watershed was instrumented to measure tile drainage volumes. The model was calibrated using water table depth and subsurface flow data over a two‐year period, while another two‐year dataset served to validate the model. DRAINMOD 5.1 accurately simulated the timing and magnitude of subsurface drainage events. The model also simulated the pattern of water table fluctuations with a good degree of accuracy. The R2 between the observed and simulated daily WTD for calibration was >0.78, and that for validation was 0.93. The corresponding coefficients of efficiency (E) were >0.74 and 0.31. The R2 and E values for calibration/validation of subsurface flow were 0.73/0.48 and 0.72/0.40, respectively. DRAINMOD simulated monthly subsurface flow quite accurately (E > 0.82 and R2 > 0.84). The model precisely simulated daily/monthly drain flow over the entire year, including the winter months. Thus DRAINMOD 5.1 performed well in simulating the hydrology of a cold region.  相似文献   

15.
ABSTRACT: With the increasing availability of digital and remotely sensed data such as land use, soil texture, and digital elevation models (DEMs), geographic information systems (GIS) have become an indispensable tool in preprocessing data sets for watershed hydrologic modeling and post processing simulation results. However, model inputs and outputs must be transferred between the model and the GIS. These transfers can be greatly simplified by incorporating the model itself into the GIS environment. To this end, a simple hydrologic model, which incorporates the curve number method of rainfall‐runoff partitioning, the ground‐water base‐flow routine, and the Muskingum flow routing procedure, was implemented on the GIS. The model interfaces directly with stream network, flow direction, and watershed boundary data generated using standard GIS terrain analysis tools; and while the model is running, various data layers may be viewed at each time step using the full display capabilities. The terrain analysis tools were first used to delineate the drainage basins and stream networks for the Susquehanna River. Then the model was used to simulate the hydrologic response of the Upper West Branch of the Susquehanna to two different storms. The simulated streamflow hydrographs compare well with the observed hydrographs at the basin outlet.  相似文献   

16.
ABSTRACT: The routing of flood waves through the Central Basin of the Passaic River in New Jersey is complex because of flat gradients and flow reversals. The one-dimensional unsteady flow program DWOPER, developed by the National Weather Service, was used to simulate flood wave movement through the Basin. A historical event was used for calibration and two synthetic events were simulated. Boundary conditions consisted of discharge hydrographs at inflow points to the study area, local flow hydrographs at interior points, and a stage discharge relation for flow over the crest of a diversion dam at the basin outlet. Manning's n values were adjusted based on stage and discharge data for the historical event; however, verification data were not available for events comparable in magnitude to the synthetic events. Aspects of the investigation reported include techniques for characterizing the flow system, model calibration, techniques for representing a tunnel diversion, and simulation results.  相似文献   

17.
Abstract: Studies to regionalize conceptual hydrologic models generally require rainfall and river flow data from multiple watersheds. Besides the considerable time (cost) to assemble and process rainfall data for many watersheds, investigators often need to choose from a number of candidate gauges, subjectively weighing the relative importance of proximity and elevation to select a representative rainfall dataset. The Unified Raingauge Dataset (URD) is a gridded daily rainfall dataset that covers the conterminous United States at 0.25 × 0.25 degrees spatial resolution and is available from 1948 to present. The objective of this study was to determine whether uncertainty in daily river flow predictions using the conceptual hydrologic model IHACRES in small to moderate size watersheds (50‐400 km2) in southern California would increase if URD gridded rainfall data were used in place of single rain gauge data to calibrate the model. Rain gauge data were obtained from the gauge nearest the watershed centroid and the gauge closest in elevation to the watershed mean elevation. Results from 20 randomly selected watersheds indicated that IHACRES calibration performance was similar using rainfall data from the URD grids and rain gauge data. There was some evidence of greater uncertainties associated with the URD calibrations in areas where topography may affect rainfall amounts. In contrast to the URD data, monthly gridded data produced by the Parameter‐Elevation Regressions on Independent Slopes Model (PRISM) includes adjustments for elevation and produces gridded values at a finer spatial resolution (4 km2). A limited test on two watersheds demonstrated that scaling the URD daily rainfall estimates to match the PRISM monthly values may improve rainfall estimates and model simulation performance.  相似文献   

18.
The Watershed Flow and Allocation model (WaterFALL®) provides segment‐specific, daily streamflow at both gaged and ungaged locations to generate the hydrologic foundation for a variety of water resources management applications. The model is designed to apply across the spatially explicit and enhanced National Hydrography Dataset (NHDPlus) stream and catchment network. To facilitate modeling at the NHDPlus catchment scale, we use an intermediate‐level rainfall‐runoff model rather than a complex process‐based model. The hydrologic model within WaterFALL simulates rainfall‐runoff processes for each catchment within a watershed and routes streamflow between catchments, while accounting for withdrawals, discharges, and onstream reservoirs within the network. The model is therefore distributed among each NHDPlus catchment within the larger selected watershed. Input parameters including climate, land use, soils, and water withdrawals and discharges are georeferenced to each catchment. The WaterFALL system includes a centralized database and server‐based environment for storing all model code, input parameters, and results in a single instance for all simulations allowing for rapid comparison between multiple scenarios. We demonstrate and validate WaterFALL within North Carolina at a variety of scales using observed streamflows to inform quantitative and qualitative measures, including hydrologic flow metrics relevant to the study of ecological flow management decisions.  相似文献   

19.
Accurate records of high‐resolution rainfall fields are essential in urban hydrology, and are lacking in many areas. We develop a high‐resolution (15 min, 1 km2) radar rainfall data set for Charlotte, North Carolina during the 2001‐2010 period using the Hydro‐NEXRAD system with radar reflectivity from the National Weather Service Weather Surveillance Radar 1988 Doppler weather radar located in Greer, South Carolina. A dense network of 71 rain gages is used for estimating and correcting radar rainfall biases. Radar rainfall estimates with daily mean field bias (MFB) correction accurately capture the spatial and temporal structure of extreme rainfall, but bias correction at finer timescales can improve cold‐season and tropical cyclone rainfall estimates. Approximately 25 rain gages are sufficient to estimate daily MFB over an area of at least 2,500 km2, suggesting that robust bias correction is feasible in many urban areas. Conditional (rain‐rate dependent) bias can be removed, but at the expense of other performance criteria such as mean square error. Hydro‐NEXRAD radar rainfall estimates are also compared with the coarser resolution (hourly, 16 km2) Stage IV operational rainfall product. Stage IV is adequate for flood water balance studies but is insufficient for applications such as urban flood modeling, in which the temporal and spatial scales of relevant hydrologic processes are short. We recommend the increased use of high‐resolution radar rainfall fields in urban hydrology.  相似文献   

20.
Accurate discharge simulation is one of the most common objectives of hydrological modeling studies. However, a good simulation of discharge is not necessarily the result of a realistic simulation of hydrological processes within the catchment. We propose an evaluation framework that considers both discharge and water balance components as evaluation criteria for calibration of the Soil and Water Assessment Tool (SWAT). In this study, we integrated average annual values of surface runoff, groundwater flow, and evapotranspiration in the model evaluation procedure to constrain the selection of good model runs for the Little River Experimental Watershed in Georgia, United States. For evaluating water balance and discharge dynamics, the Nash‐Sutcliffe efficiency (NSE) and percent bias (PBIAS) were used. In addition, the ratio of root mean square error and standard deviation of measured data (RSR) was calculated for individual segments of the flow duration curve to identify the best model runs in terms of discharge magnitude. Our results indicate that good statistics for discharge do not guarantee realistic simulations of individual water balance components. Therefore, we recommend constraining the ranges of water balance components to achieve a more realistic simulation of the entire hydrological system, even if tradeoffs between good statistics for discharge simulations and reasonable amounts of the water balance components are unavoidable. Editor's note : This paper is part of the featured series on SWAT Applications for Emerging Hydrologic and Water Quality Challenges. See the February 2017 issue for the introduction and background to the series.  相似文献   

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

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