首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 109 毫秒
1.
The Soil Conservation Service curve number (CN) method commonly uses three discrete levels of soil antecedent moisture condition (AMC), defined by the 5‐day antecedent rainfall depth, to describe soil moisture prior to a runoff event. However, this way may not adequately represent soil water conditions of fields and watersheds in the Loess Plateau of China. The objectives of this study were: (1) to determine the effective soil moisture depth to which the CN is most related; (2) to evaluate a discrete and a linear relationship between AMC and soil moisture; and (3) to develop an equation between CN and soil moisture to predict runoff better for the climatic and soil conditions of the Loess Plateau of China. The dataset consisted of 10 years of rainfall, runoff and soil moisture measurements from four experimental plots cropped with millet, pasture and potatoes. Results indicate that the standard CN method underestimated runoff depths for 85 of the 98 observed plot‐runoff events, with a model efficiency E of only 0·243. For our experimental conditions, the discrete and linear approaches improved runoff estimation, but still underestimated most runoff events, with E values of 0·428 and 0·445 respectively. Based on the measured CN values and soil moisture values in the top 15 cm of the soil, a non‐linear equation was developed that predicted runoff better with an E value of 0·779. This modified CN equation was the most appropriate for runoff prediction in the study area, but may need adjustments for local conditions in the Loess Plateau of China. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

2.
Runoff estimations based on the standard USDA–NRCS curve number (CN) table without calibration have a tendency to give inaccurate results when the CN values are applied in South Korea which has many high slope watersheds and that has a continental monsoon climate. Particularly for the design flood estimation, accurately calibrated CN values are required because the estimated peak flow is very sensitive to the selection of CN. However, the lack of flood data makes it difficult to calibrate and assign runoff CNs to Korean watersheds. Even if sufficient data are available to estimate CN values, it is also difficult to obtain the direct flows by separating base flows from total runoff hydrographs due to the temporal irregularity of rainfall events and the resulting complex pattern of runoff. Therefore, an alternative method for estimating CNs needs to be developed to overcome these issues. The purpose of this study is to present a method for estimating runoff CNs using the soil and water assessment tool (SWAT) model which can take into account watershed heterogeneities such as climate conditions, land use and soil types. The proposed CN estimation method uses the simulated flow data by SWAT instead of using measured flow data. This method has advantages in estimating CN values spatially for each subbasin division considering watershed characteristics. The use of daily data can reduce the sensitivity to the abnormality that is commonly involved in flow data with a small time scale. The SWAT‐based CN estimation method, combined with the asymptotic CN method, was applied to the Chungju dam watershed in South Korea. A regression equation was then developed from this approach, which was used to estimate CN values that decrease exponentially as rainfall amounts increase and that converge to 60·6 and 79·4 without and with considering subsurface lateral flow, respectively. Furthermore, the CN values for the antecedent moisture conditions were determined using the probabilistic approach. The CN associated with the 50% probability for the Chungju dam watershed is 87·8 which can be taken to be representative of antecedent moisture condition (AMC) II. The CNI and CNIII associated with 90% and 10% probabilities are 78·9 and 94·1, respectively. The estimated CNII = 87·8 differs markedly from the geographic information system (GIS)‐based CN 65·0, which implies that the standard USDA–NRCS CN method should be calibrated to the studied area of interest. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

3.
This paper describes the results of benchmark testing of land use change impact on direct runoff using Soil Conservation Service-Curve Number (SCS-CN) model in two ungauged neighbouring urban watersheds (Çınar and Kadıyakuplu) in Istanbul, Turkey. To examine this impact, the model was applied to daily rainfall data using three different dated (1982, 1996 and 2012) hydrological soil groups and land use of the two ungauged urban watersheds. Finally, the impact of land use change and model performance were evaluated with the rainfall-runoff regression, the coefficient of determination and the NSE test using benchmark runoff data based on 1982 land use conditions. The results of the analysis indicate that the changing of land use types from natural surfaces to impervious surfaces has a significant impact on surface runoff. Additionally, remarkable spatial variations of the land use changes and their impact on the runoff in 1996 and 2012 were more detected in the Çınar watershed compared with the Kadıyakuplu watershed. The planning decision on land use of the watersheds, has vital role in these differences. The results of this research also reveal that change to intensive land use in urban watersheds has a significantly larger impact on runoff generation than those rainfall.  相似文献   

4.
Spatial and seasonal variations of curve number (CN) and initial abstraction ratio (λ) in a watershed can result in inaccurate runoff volume estimations when using the US Natural Resources Conservation Service (SCS-CN) method with constant values for these parameters. In this paper, parameters of CN and λ are considered as calibration parameters and the sensitivity of estimated runoff to these parameters using the SCS-CN method is scrutinized. To incorporate the uncertainty associated with CN and λ, fuzzy linear regression (FLR) is applied to derive the relationships of CN and λ with rainfall depth (P) by employing a large dataset of storm events from four watersheds in Iran. Results indicate that the proposed approach provides more accuracy in estimation of runoff volume compared to the SCS method with constant values of CN and λ, and gives a straightforward technique for evaluating the hydrological effects of CN, λ, and P on runoff volume.  相似文献   

5.
This paper investigates the variation of the popular curve number (CN) values given in the National Engineering Hand Book–Section 4 (NEH‐4) of the Soil Conservation Service (SCS) with antecedent moisture condition (AMC) and soil type. Using the volumetric concept, involving soil, water, and air, a significant condensation of the NEH‐4 tables is achieved. This leads to a procedure for determination of CN for gauged as well as ungauged watersheds. The rainfall‐runoff events derived from daily data of four Indian watersheds exhibited a power relation between the potential maximum retention or CN and the 5‐day antecedent rainfall amount. Including this power relation, the SCS‐CN method was modified. This modification also eliminates the problem of sudden jumps from one AMC level to the other. The runoff values predicted using the modified method and the existing method utilizing the NEH‐4 AMC criteria yielded similar results. Copyright © 2006 John Wiley & Sons, Ltd.  相似文献   

6.
Exploring the chemical characterization of dissolved organic matter (DOM) is important for understanding the fate of laterally transported organic matter in watersheds. We hypothesized that differences in water-extractable organic matter (WEOM) in soils of varying land uses and rainfall events may significantly affect the quality and the quantity of stream DOM. To test our hypotheses, characteristics of rainfall-runoff DOM and WEOM of source materials (topsoil from different land uses and gullies, as well as typical vegetation) were investigated at two adjacent catchments in the Loess Plateau of China, using ultraviolet–visible absorbance and excitation emission matrix fluorescence with parallel factor analysis (PARAFAC). Results indicated that land-use types may significantly affect the chemical composition of soil WEOM, including its aromaticity, molecular weight, and degree of humification. The PARAFAC analysis demonstrated that the soils and stream water were dominated by terrestrial/allochthonous humic-like substances and microbial transformable humic-like fluorophores. Shifts in the fluorescence properties of stream DOM suggested a pronounced change in the relative proportion of allochthonous versus autochthonous material under different rainfall patterns and land uses. For example, high proportions of forestland could provide more allochthonous DOM input. This study highlights the relevance of soils and hydrological dynamics on the composition and fluxes of DOM issuing from watersheds. The composition of DOM in soils was influenced by land-use type. Precipitation patterns influenced the proportion of terrestrial versus microbial origins of DOM in surface runoff. Contributions of allochthonous, terrestrially derived DOM inputs were highest from forested landscapes.  相似文献   

7.
Knowledge of the effective impervious area (EIA) or the degree to which impervious surfaces are hydraulically connected to the drainage system is useful for improving hydrological and environmental models and assessing the effectiveness of green stormwater infrastructure in urban watersheds. The goal of this research is to develop a method to estimate EIA fraction in urban watersheds using readily available data. Since EIA is dependent on rainfall–runoff response and cannot be solely determined based on the physical characteristics of a watershed, the EIA is linked with the asymptotic curve number (CN), a watershed index that represents runoff characteristics. In order for the method to be applicable to ungauged watersheds, the asymptotic CN is predicted using land cover and soil data from 35 urban catchments in Minnesota and Texas, USA. Similar data from 11 other urban catchments in Wisconsin and Texas, USA, are used to validate the results. A set of runoff depth versus EIA fraction curves is also developed to assess the impact of EIA reduction on discharge from an urban watershed in land-use planning studies.  相似文献   

8.
Linking quickflow response to subsurface state can improve our understanding of runoff processes that drive emergent catchment behaviour. We investigated the formation of non-linear quickflows in three forested headwater catchments and also explored unsaturated and saturated storage dynamics, and likely runoff generation mechanisms that contributed to threshold formation. Our analyses focused on two reference watersheds at the Coweeta Hydrologic Laboratory (CHL) in western North Carolina, USA, and one reference watershed at the Susquehanna Shale Hills Critical Zone Observatory (SHW) in Central Pennsylvania, USA, with available hourly soil moisture, groundwater, streamflow, and precipitation time series over several years. Our study objectives were to characterise (a) non-linear runoff response as a function of storm characteristics and antecedent conditions, (b) the critical levels of shallow unsaturated and saturated storage that lead to hourly flow response, and (c) runoff mechanisms contributing to rapidly increasing quickflow using measurements of soil moisture and groundwater. We found that maximum hourly rainfall did not significantly contribute to quickflow production in our sites, in contrast to prior studies, due to highly conductive forest soils. Soil moisture and groundwater dynamics measured in hydrologically representative areas of the hillslope showed that variable subsurface states could contribute to non-linear runoff behaviour. Quickflow generation in watersheds at CHL were dominated by both saturated and unsaturated pathways, but the relative contributions of each pathway varied between catchments. In contrast, quickflow was almost entirely related to groundwater fluctuations at SHW. We showed that co-located measurements of soil moisture and groundwater supplement threshold analyses providing stronger prediction and understanding of quickflow generation and indicate dominant runoff processes.  相似文献   

9.
Over the last thirty years, French hydrologists have observed flow characteristics at the outlets of small tropical watersheds, characteristics which vary between arid zones, humid zones and equatorial rainforest zones. The main physical phenomena related to the various types of runoff generation have been observed and analysed. The generation of runoff and subsurface runoff following tropical rainstorms depends on soil surface formations, texture and structure as well as vegetation cover, geological substrata, weathered material and climate.

Analysis is also made of experimental rainfall-runoff relationships, the shape of flood hydrographs and the conditions necessary for the occurrence of steady base flow. Scale effects between experimental plot and watershed have been studied to identify any effect this may have upon runoff generation and flow. Tables detail watershed characteristics and relationships between these and rainfall and runoff.  相似文献   


10.
Temporally weighted average curve number method for daily runoff simulation   总被引:1,自引:0,他引:1  
Nam Won Kim  Jeongwoo Lee 《水文研究》2008,22(25):4936-4948
The modified Soil Conservation Service curve number (CN) method is widely used in long‐term continuous models to predict daily surface runoff. However, it has been shown that this method gives poor results in reproducing peak flows in high rainfall periods. This is because there is an inaccuracy stemming from the model algorithm as it adjusts the daily runoff curve number as a function of soil moisture content at the end of the previous day. This paper proposes an alternative daily based curve number technique that can provide better prediction of daily runoff during the high flow season. The proposed method uses the temporally weighted average curve number (TWA‐CN) to estimate daily surface runoff, while considering the effect of rainfall during a given day as well as the antecedent soil moisture condition. To test the applicability of the TWA‐CN method, it was incorporated with the long‐term, continuous simulation watershed models SWAT and SWAT‐G. Simulations were conducted for the Miho River watershed located in the middle of South Korea. The graphical displays and statistics of the determination coefficient (R2) and the Nash–Sutcliffe model efficiency (NSE) of the observed and simulated daily runoff indicated that the modified SWAT with the TWA‐CN method may provide better runoff prediction (R2 = 0·837, NSE = 0·833) than the original SWAT (R2 = 0·815, NSE = 0·824). Likewise, the determination coefficient (R2 = 0·816) and the Nash–Sutcliffe efficiency (NSE = 0·834) for the modified SWAT‐G are also higher than the original version (R2 = 0·782, NSE = 0·825). It is expected that the improved capability in predicting surface runoff using the suggested CN estimate method will provide a sound contribution to the accurate simulations of water yield. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

11.
Because the traditional Soil Conservation Service curve‐number (SCS‐CN) approach continues to be used ubiquitously in water quality models, new application methods are needed that are consistent with variable source area (VSA) hydrological processes in the landscape. We developed and tested a distributed approach for applying the traditional SCS‐CN equation to watersheds where VSA hydrology is a dominant process. Predicting the location of source areas is important for watershed planning because restricting potentially polluting activities from runoff source areas is fundamental to controlling non‐point‐source pollution. The method presented here used the traditional SCS‐CN approach to predict runoff volume and spatial extent of saturated areas and a topographic index, like that used in TOPMODEL, to distribute runoff source areas through watersheds. The resulting distributed CN–VSA method was applied to two subwatersheds of the Delaware basin in the Catskill Mountains region of New York State and one watershed in south‐eastern Australia to produce runoff‐probability maps. Observed saturated area locations in the watersheds agreed with the distributed CN–VSA method. Results showed good agreement with those obtained from the previously validated soil moisture routing (SMR) model. When compared with the traditional SCS‐CN method, the distributed CN–VSA method predicted a similar total volume of runoff, but vastly different locations of runoff generation. Thus, the distributed CN–VSA approach provides a physically based method that is simple enough to be incorporated into water quality models, and other tools that currently use the traditional SCS–CN method, while still adhering to the principles of VSA hydrology. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

12.
Abstract

The normalized antecedent precipitation index (NAPI) model by Heggen for the prediction of runoff yield is analytically derived from the water balance equation. Heggen's model has been simplified further to a rational form and its performance verified with the Soil Conservation Service Curve Number (SCS-CN) model. The simplified model has three coefficients specific to a watershed, and requires two inputs: rainfall and the derived parameter, NAPI. The characteristic behaviour of the NAPI has resonance with the curve number (CN) of the SCS model. The proposed NAPI model was applied to three watersheds in the semi-arid region of India to simulate runoff yield. The model showed improved correlation between the observed and predicted runoff data compared to the SCS-CN model. The F test and paired t test also confirmed the reliability of the model with significance levels of 0.01 and 0.001%, respectively. The proposed model could be used successfully for rainfall–runoff modelling in a watershed.

Citation Ali, S., Ghosh, N. C. & Singh, R. (2010) Rainfall–runoff simulation using a normalized antecedent precipitation index. Hydrol. Sci. J. 55(2), 266–274.  相似文献   

13.
Since 1999, large-scale ecosystem restoration has been implemented in the Loess Plateau, effectively increasing regional vegetation coverage. Vegetation restoration has significantly elevated the saturated hydraulic conductivity (Ks) of the near-surface soil layers and increased the vertical heterogeneity of the Ks profile. Many studies have examined the change of runoff due to revegetation, yet the impacts of Ks profile on the soil moisture distribution and runoff generation processes were less explored. In this study, numerical simulations were conducted to investigate how changes in the Ks profile caused by vegetation restoration influenced the hydrological responses at event scale. The numerical simulation results show that the increase of surface Ks caused by vegetation restoration can effectively reduce runoff at event scale. Moreover, the enhancement of vertical heterogeneity of Ks profiles can significantly change the vertical profile of soil water content, prompting more water to percolate into the deep soil layer. When rainfall exceeds a threshold, the accumulation of soil water above the relatively less permeable layer can cause short-term saturation in shallow soil layers, resulting in a transient perched water table. As a result, after the vegetation restoration in the Loess Plateau, though Horton overland flow is still the main runoff generation mechanism, there is a possibility of the emergence of Dunne overland flow under the high vegetation coverage (e.g., NDVI larger than 0.5). This emergence of new runoff generation mechanism, saturation excess runoff, in the Loess Plateau due to the vegetation restoration could provide scientific guidance for water and sediment movement, soil and water conservation practices, and desertification control in the Loess Plateau.  相似文献   

14.
ABSTRACT

When discharge measurements are not available, design of water structures relies on using frequency analysis of rainfall data and applying a rainfall–runoff model to estimate a hydrograph. The Soil Conservation Service (SCS) method estimates the design hydrograph first through a rainfall–runoff transformation and next by propagating runoff to the basin outlet via the SCS unit hydrograph (UH) method. The method uses two parameters, the Curve Number (CN) and the time of concentration (Tc). However, in data-scarce areas, the calibration of CN and Tc from nearby gauged watersheds is limited and subject to high uncertainties. Therefore, the inherent uncertainty/variability of the SCS parameters may have considerable ramifications on the safety of design. In this research, a reliability approach is used to evaluate the impact of incorporating the uncertainty of CN and Tc in flood design. The sensitivity of the probabilistic outcome against the uncertainty of input parameters is calculated using the First Order Reliability Method (FORM). The results of FORM are compared with the conventional SCS results, taking solely the uncertainty of the rainfall event. The relative importance of the uncertainty of the SCS parameters is also estimated. It is found that the conventional approach, used by many practitioners, might grossly underestimate the risk of failure of water structures, due to neglecting the probabilistic nature of the SCS parameters and especially the Curve Number. The most predominant factors against which the SCS-CN method is highly uncertain are when the average rainfall value is low (less than 20 mm) or its coefficient of variation is not significant (less than 0.5), i.e. when the resulting rainfall at the design return period is low. A case study is presented for Egypt using rainfall data and CN values driven from satellite information, to determine the regions of acceptance of the SCS-CN method.
EDITOR D. Koutsoyiannis; ASSOCIATE EDITOR A. Efstratiadis  相似文献   

15.
Evaluating the benefits of sediment and runoff reduction in different vegetation types is essential for studying the mechanisms of soil and water conservation on the Loess Plateau.The experiment was conducted in shrub-grass plots with nine levels of mixed vegetation coverage from 0%to 70%,three slopes(10,15,and 20)and two rainfall intensities(1.0 and 2.5 mm/min).The results showed that the vegetation coverage and slope gradient significantly affect runoff and sediment yield.Shrub-grass vegetation coverage had a significant effect on the runoff start-time,runoff flow velocity,runoff rate,and soil erosion rate on hillslopes.Mixed vegetation coverage could effectively delay the runoff starttime and decrease the runoff flow velocity.However,the effects of the slope gradient on runoff and sediment yield are opposite to those of vegetation coverage.Shrub-grass vegetation coverage could effectively increase runoff and sediment yield reduction benefits,while their benefits were affected by the rainfall intensity.At the 1.0 mm/min rainfall intensity,the reduction in the sediment production rate was greater than that under the 2.5 mm/min intensity.However,when the shrub-grass vegetation coverage exceeded 42%,the runoff reduction benefit was more obvious at higher rainfall intensities.The cumulative sediment yield increased with increasing cumulative runoff,and the rate of increase in the cumulative runoff was greater than that of the cumulative sediment yield with increasing of shrub-grass vegetation coverage.Moreover,there was a power function relationship between cumulative sediment yield and cumulative runoff yield(P<0.05).Our paper is expected to provide a good reference on the ecological environment and vegetation construction on the Loess Plateau.  相似文献   

16.
The Soil Conservation Service Curve Number (CN) method is routinely used to estimate the effects of forest fires on hydrological response. However, despite recent efforts, CN values are still not well known for burned conditions. A major forest fire in Attica, Greece, which affected the Lykorrema stream experimental watershed, provided an opportunity for the estimation of post-fire CN variation using detailed pre-fire and post-fire rainfall–runoff datasets. The CN values for both periods were estimated and compared using a wide range of available methods. Methods considering the spatial variability of soil-cover complexes were also used to investigate the effect of spatial heterogeneity. The post-fire watershed response changed from complacent to standard. Direct runoff depths and peak flows increased by a factor of more than 7.7 and 11.8, respectively. On average, the estimated post-fire CN values for the studied soil-cover complexes increased by about 25 units. This study may assist the improvement of existing post-fire hydrological assessment tools.  相似文献   

17.
Abstract

The objective of this study is to develop a Modified Rational Equation (MoRE) that combines the advantages of the Rational Equation (e.g. simplicity and global acceptance) and those of the standard US Department of Agriculture (USDA) Soil Conservation Service (SCS) curve number (CN) method (e.g. easy parameterization and extensive verification across the world). Herein, the hypothesis is that the MoRE is more accurate, consistent and robust than the SCS-CN method and its improved versions in predicting runoff in watersheds with limited data. The MoRE was designed to have a simple structure that is described by four intrinsic parameters: CN, permanent wilting point, field capacity and saturation soil moisture, and does not include initial abstraction as a variable. An evaluation of 77 USDA small agricultural watersheds indicated that CN of the MoRE has different physical meanings from CN of the SCS-CN method. The MoRE (mean Nash-Sutcliffe coefficient, E > 0.73) performed better than the SCS-CN (mean E < 0.32) and the four improved models (mean E < 0.56) in reproducing the runoff of the study watersheds. Performance of all six models varied greatly between watersheds, as well as between events, but was independent of watershed drainage area. However, the model performances tend to be better for watersheds and/or events with a runoff-to-rainfall ratio of between 0.1 and 0.3 than for those with a ratio outside this range. The MoRE has the most consistent and robust performance.

Editor D. Koutsoyiannis; Associate editor I. Nalbantis

Citation Wang, X., Liu, T., and Yang, W., 2012. Development of a robust runoff-prediction model by fusing the rational equation and a modified SCS-CN method. Hydrological Sciences Journal, 57 (6), 1118–1140.  相似文献   

18.
The impacts of climate change on storm runoff and erosion in Mediterranean watersheds are difficult to assess due to the expected increase in storm frequency coupled with a decrease in total rainfall and soil moisture, added to positive or negative changes to different types of vegetation cover. This report, the second part of a two‐part article, addresses this issue by analysing the sensitivity of runoff and erosion to incremental degrees of change (from ? 20 to + 20%) to storm rainfall, pre‐storm soil moisture, and vegetation cover, in two Mediterranean watersheds, using the MEFIDIS model. The main results point to the high sensitivity of storm runoff and peak runoff rates to changes in storm rainfall (2·2% per 1% change) and, to a lesser degree, to soil water content (?1·2% per 1% change). Catchment sediment yield shows a greater sensitivity than within‐watershed erosion rates to both parameters: 7·8 versus 4·0% per 1% change for storm rainfall, and ? 4·9 versus ? 2·3% per 1% change for soil water content, indicating an increase in sensitivity with spatial scale due to changes to sediment connectivity within the catchment. Runoff and erosion showed a relatively low sensitivity to changes in vegetation cover. Finally, the shallow soils in one of the catchments led to a greater sensitivity to changes in storm rainfall and soil moisture. Overall, the results indicate that decreasing soil moisture levels caused by climate change could be sufficient to offset the impact of greater storm intensity in Mediterranean watersheds. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

19.
In the Soil Conservation Service Curve Number (SCS‐CN) method, the three levels of antecedent moisture condition (AMC) permit unreasonable sudden jumps in curve numbers, which result into corresponding jumps in the estimated runoff. A few recently developed SCS‐CN‐based models obviate this problem, yet they have several limitations. In this study, such a model incorporating a continuous function for antecedent moisture has been presented. It has several advantages over the other existing SCS‐CN‐based models. Its application to a large dataset from US watersheds showed to perform better than the existing SCS‐CN method and the others based on it. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

20.
The analysis of the physical processes involved in a conceptual model of soil water content balance is addressed with the objective of its application as a component of rainfall–runoff modelling. The model uses routinely measured meteorological variables (rainfall and air temperature) and incorporates a limited number of significant parameters. Its performance in estimating the soil moisture temporal pattern was tested through local measurements of volumetric water content carried out continuously on an experimental plot located in central Italy. The analysis was carried out for different periods in order to test both the representation of infiltration at the short time‐scale and drainage and evapotranspiration processes at the long time‐scale. A robust conceptual model was identified that incorporated the Green–Ampt approach for infiltration and a gravity‐driven approximation for drainage. A sensitivity analysis was performed for the selected model to assess the model robustness and to identify the more significant parameters involved in the principal processes that control the soil moisture temporal pattern. The usefulness of the selected model was tested for the estimation of the initial wetness conditions for rainfall–runoff modelling at the catchment scale. Specifically, the runoff characteristics (runoff depth and peak discharge) were found to be dependent on the pre‐event surface soil moisture. Both observed values and those estimated by the model gave good results. On the contrary, with the antecedent wetness conditions furnished by two versions of the antecedent precipitation index (API), large errors were obtained. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

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

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