首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Hysteresis is a common feature exhibited in hydraulic properties of an unsaturated soil. The movement of wetting front and the hysteresis effect are important factors which impact the shear strength of the unsaturated soil and the mechanics of shallow landslides. These failures are mainly triggered by the deepening of the wetting front accompanied by a decrease in matric suction induced by infiltration. This research establishes a method for determining a stability analysis of unsaturated infinite soil slopes, integrating the influence of infiltration and the water retention curve hysteresis. Furthermore, the present stability analysis method including the infiltration model and the advanced Mohr–Coulomb failure criterion calculates the variations of the safety factor (FS) in accordance with different slope angle, depth and hydrological processes. The experimentally measured data on the effect of hysteresis are also carried out for comparison. Numerical analyses, employing both wetting and drying hydraulic behaviour of unsaturated soil, are performed to study the difference in soil‐water content as observed in the experiments. The simulating approximations also fully responded to the experimental data of sand box. The results suggest that the hysteresis behaviour affect the distribution of soil‐water content within the slope indeed. The hysteresis made the FS values a remarkable recovery during the period of non‐rainfall in a rainfall event. The appropriate hydraulic properties of soil (i.e. wetting or drying) should be used in accordance with the processes that unsaturated soil actually experience. This method will enable us to acquire more accurate matric suction head and the unsaturated soil‐shear strength as it changes with the hysteretic flow, in order to calculate into the stability analysis of shallow landslides. An advanced understanding of the process mechanism afforded by this method is critical to realizing a reliable and appropriate design for slope stabilization. It also offers some immediate reference information to the disaster reduction department of the government. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

2.
This study examined the influence of tidally‐induced oscillations of the beach water table in regulating beach surface moisture dynamics. A series of laboratory experiments were conducted to assess the influence of hysteresis and transient flow effects on surface moisture variability. The experimental apparatus utilized a column of well‐sorted fine sand partially immersed in a reservoir of water. The water level in the reservoir was raised and lowered via a diaphragm‐metering pump to simulate tidally induced fluctuations of the water table, and the moisture content profile within the column was monitored using an array of Delta‐T probes. Moisture contents at specific elevations within the column were utilized as proxies to represent various ‘surface’ elevations (relative to the high water table). Results indicate that surface moisture content behaves in a distinctly hysteretic manner. Examination of water flow scanning curves illustrated that for all surface elevations considered, higher moisture contents for a given pressure head occurred during the drying cycle than during the wetting cycle. This observation is particularly evident with shallow surface elevations (i.e. water table close to the surface) where the Haines Jump phenomenon was found to have a significant influence on moisture content dynamics. Additionally, an assessment of the accuracy of hysteretic and non‐hysteretic models to predict the measured moisture contents demonstrated that hysteretic simulations consistently provide a better representation of the observed moisture contents than non‐hysteretic simulations. A time lag was found between the respective maxima and minima in water table elevation surface moisture content. At the near surface water table positions the time lag ranged between 30 and 100 minutes, and it increased to 240 minutes (four hours) with the high water table at 60 cm below the surface. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

3.
Hu Liu  Wenzhi Zhao  Zhibin He  Jintao Liu 《水文研究》2015,29(15):3328-3341
A combination of field measurements, continuous monitoring and numerical modelling was used to evaluate soil moisture regimes at four sites across a landscape gradient of the Heihe River Basin. Recorded data of precipitation, irrigation and floods were used to build the model, and an optimization technique was employed to calibrate the parameters. Based on the optimized parameters and estimates of future scenarios, the modelling structure was employed to predict the changes in the growing season soil moisture regimes due to climate change and intensive management. The results suggest that the upper‐reach Yeniugou and Xishui sites will become wetter because of alterations in the precipitation regime, and this trend could be strengthened by the expected amplified interannual variability. Precipitation features at middle‐reach Linze and lower‐reach Ejina, however, are not expected to change in the future. We assumed that a more water‐saving irrigation system will replace the current traditional one, and hence, the soil moisture probability density function at the Linze site would tend to be narrowed to ranges around either the wilting point or the point of incipient water stress, depending on how the intervention point and target level are settled. Ejina is expected to experience the most extreme ecological conversion effects in the future, and soil moisture would be more frequently recharged by water delivery. However, the soil moisture regime would not change much because of the poor water‐holding capacity and intensive evaporation. The revealed patterns and predicted shifts in soil moisture dynamics could provide a useful reference for identifying robust long‐term water resource management strategies for the Heihe River Basin. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

4.
Infiltration into frozen soil is a key hydrological process in cold regions. Although the mechanisms behind point‐scale infiltration into frozen soil are relatively well understood, questions remain about upscaling point‐scale results to estimate hillslope‐scale run‐off generation. Here, we tackle this question by combining laboratory, field, and modelling experiments. Six large (0.30‐m diameter by 0.35‐m deep) soil cores were extracted from an experimental hillslope on the Canadian Prairies. In the laboratory, we measured run‐off and infiltration rates of the cores for two antecedent moisture conditions under snowmelt rates and diurnal freeze–thaw conditions observed on the same hillslope. We combined the infiltration data with spatially variable data from the hillslope, to parameterise a surface run‐off redistribution model. We used the model to determine how spatial patterns of soil water content, snowpack water equivalent (SWE), and snowmelt rates affect the spatial variability of infiltration and hydrological connectivity over frozen soil. Our experiments showed that antecedent moisture conditions of the frozen soil affected infiltration rates by limiting the initial soil storage capacity and infiltration front penetration depth. However, shallow depths of infiltration and refreezing created saturated conditions at the surface for dry and wet antecedent conditions, resulting in similar final infiltration rates (0.3 mm hr?1). On the hillslope‐scale, the spatial variability of snowmelt rates controlled the development of hydrological connectivity during the 2014 spring melt, whereas SWE and antecedent soil moisture were unimportant. Geostatistical analysis showed that this was because SWE variability and antecedent moisture variability occurred at distances shorter than that of topographic variability, whereas melt variability occurred at distances longer than that of topographic variability. The importance of spatial controls will shift for differing locations and winter conditions. Overall, our results suggest that run‐off connectivity is determined by (a) a pre‐fill phase, during which a thin surface soil layer wets up, refreezes, and saturates, before infiltration excess run‐off is generated and (b) a subsequent fill‐and‐spill phase on the surface that drives hillslope‐scale run‐off.  相似文献   

5.
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.  相似文献   

6.
A periodic vertical movement of the groundwater table results in a subsequent cyclic response of the water content and pressure profiles in the vadose zone. The sequence of periodic wetting and drying processes can be affected by hysteresis effects in this zone. A one-dimensional saturated/unsaturated flow model based on Richards’ equation and the Mualem (Soil Sci. 137 (1984) 283) hysteresis model is formulated which can take into account multi-cycle hysteresis effects in the relation between capillary pressure and water content. The numerical integration of the unsaturated flow equation is based on a Galerkin-type finite element method. The flow domain is discretised by finite elements with linear shape functions. Simulations start with static water content and pressure profiles, which correspond to either a boundary drying or wetting retention curve. To facilitate the numerical solution of the hysteretic case an implicit non-iterative procedure was chosen for the solution of the nonlinear differential equation. Laboratory experiments were performed with a vertical sand column by imposing a high frequency periodic pressure head at the lower end of the column. The total water volume in the column, and the periodic water content profile averaged over time were measured. The boundary drying and wetting curves of the relation between water content and capillary pressure were determined by independent experiments. The simulations of the experimental conditions show a clear effect of the hysteresis phenomenon on the water content profile. The simulations with hysteresis agree well with the measurements. Computed dimensionless water content profiles are presented for different oscillation frequencies with and without consideration of hysteresis.  相似文献   

7.
Often the soil hydraulic parameters are obtained by the inversion of measured data (e.g. soil moisture, pressure head, and cumulative infiltration, etc.). However, the inverse problem in unsaturated zone is ill‐posed due to various reasons, and hence the parameters become non‐unique. The presence of multiple soil layers brings the additional complexities in the inverse modelling. The generalized likelihood uncertainty estimate (GLUE) is a useful approach to estimate the parameters and their uncertainty when dealing with soil moisture dynamics which is a highly non‐linear problem. Because the estimated parameters depend on the modelling scale, inverse modelling carried out on laboratory data and field data may provide independent estimates. The objective of this paper is to compare the parameters and their uncertainty estimated through experiments in the laboratory and in the field and to assess which of the soil hydraulic parameters are independent of the experiment. The first two layers in the field site are characterized by Loamy sand and Loamy. The mean soil moisture and pressure head at three depths are measured with an interval of half hour for a period of 1 week using the evaporation method for the laboratory experiment, whereas soil moisture at three different depths (60, 110, and 200 cm) is measured with an interval of 1 h for 2 years for the field experiment. A one‐dimensional soil moisture model on the basis of the finite difference method was used. The calibration and validation are approximately for 1 year each. The model performance was found to be good with root mean square error (RMSE) varying from 2 to 4 cm3 cm?3. It is found from the two experiments that mean and uncertainty in the saturated soil moisture (θs) and shape parameter (n) of van Genuchten equations are similar for both the soil types. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

8.
Accurate estimation of the soil water balance (SWB) is important for a number of applications (e.g. environmental, meteorological, agronomical and hydrological). The objective of this study was to develop and test techniques for the estimation of soil water fluxes and SWB components (particularly infiltration, evaporation and drainage below the root zone) from soil water records. The work presented here is based on profile soil moisture data measured using dielectric methods, at 30‐min resolution, at an experimental site with different vegetation covers (barley, sunflower and bare soil). Estimates of infiltration were derived by assuming that observed gains in the soil profile water content during rainfall were due to infiltration. Inaccuracies related to diurnal fluctuations present in the dielectric‐based soil water records are resolved by filtering the data with adequate threshold values. Inconsistencies caused by the redistribution of water after rain events were corrected by allowing for a redistribution period before computing water gains. Estimates of evaporation and drainage were derived from water losses above and below the deepest zero flux plane (ZFP), respectively. The evaporation estimates for the sunflower field were compared to evaporation data obtained with an eddy covariance (EC) system located elsewhere in the field. The EC estimate of total evaporation for the growing season was about 25% larger than that derived from the soil water records. This was consistent with differences in crop growth (based on direct measurements of biomass, and field mapping of vegetation using laser altimetry) between the EC footprint and the area of the field used for soil moisture monitoring. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

9.
Infiltration systems are widely used as an effective urban stormwater control measure. Most design methods and models roughly approximate the complex physical flow processes in these systems using empirical equations and fixed infiltration rates to calculate emptying times from full. Sophisticated variably saturated flow models are available, but rarely applied owing to their complexity. This paper describes the development and testing of an integrated one‐dimensional model of flow through the porous storage of a typical infiltration system and surrounding soils. The model accounts for the depth in the storage, surrounding soil moisture conditions and the interaction between the storage and surrounding soil. It is a front‐tracking model that innovatively combines a soil‐moisture‐based solution of Richard's equation for unsaturated flow with piston flow through a saturated zone as well as a reservoir equation for flow through a porous storage. This allows the use of a simple non‐iterative numerical solution that can handle ponded infiltration into dry soils. The model is more rigorous than approximate stormwater infiltration system models and could therefore be valuable in everyday practice. A range of test cases commonly used to test soil water flow models for infiltration in unsaturated conditions, drainage from saturation and infiltration under ponded conditions were used to test the model along with an experiment with variable depth in a porous storage over saturated conditions. Results show that the model produces a good fit to the observed data, analytical solutions and Hydrus. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

10.
11.
To study the effects of biological soil crusts (BSCs) on hydrological processes and their implications for disturbance in the Mu Us Sandland, the water infiltration, evaporation and soil moisture of high coverage (100% BSCs), middle coverage (40% BSCs) and low coverage (0% BSCs, bare sand) of moss‐dominated crusts were conducted in this study, respectively. The conclusions are as follows: (1) the main effects of moss‐dominated crusts in the Mu Us Sandland on the infiltration of rainwater were to reduce the infiltration depths and to retain the limited rainwater in shallow soil; (2) moss‐dominated crusts have no significant effects on daily evaporation when the volumetric water content at 4 cm depth in 100% BSCs (VWC4) was over 24.7%, on enhanced daily evaporation when the VWC4 ranged from 6.5% to 24.7% and on reduced daily evaporation when the VWC4 was less than 6.5%; and (3) decreasing the coverage of moss‐dominated crusts (from 100% to 40%) did not significantly change its effects on infiltration, evaporation and soil moisture. Our results demonstrated that for the growth and regeneration of shrubs, which were dominated by Artemisia ordosica in the Mu Us Sandland, high coverage of moss‐dominated crusts has negative effects on hydrological processes, and these negative effects could not be significantly reduced by decreasing the coverage of moss‐dominated crusts from 100% to 40%. Therefore, for the sustained and healthy development of shrub communities in the Mu Us Sandland, it is necessary to take appropriate measures for the well‐developed BSCs in the sites with high vegetation coverage in the rainy season. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

12.
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.  相似文献   

13.
This study presents a versatile index for the quantification of hysteretic loops between hydrological variables at the runoff event timescale. The conceptual development of the index is based on a normalization of the input data and the computation of definite integrals at fixed intervals of the independent variable. The sum, the minimum and the maximum of the differences between integrals computed for the rising and the falling curves provide information on the direction, the shape and the extent of the loop. The index was tested with synthetic data and field data from experimental catchments in Northern Italy. Hysteretic relations between streamflow (the independent variable) and soil moisture, depth to water table, isotopic composition and electrical conductivity of stream water (dependent variables) were correctly identified and quantified by the index. The objective quantification of hysteresis by the index allows for the automatic classification of hysteretic loops and thus the determination of differences in hydrological responses during different events. The index was also used to examine the seasonal dynamics in the relation between streamflow and soil moisture and captured the switch in the direction of the loop with changes in event size and antecedent wetness conditions. The sensitivity of the index to the temporal resolution of the measurements and measurement errors was also tested. The index can successfully quantify hysteresis, except for very noisy data or when the temporal resolution of the measurements is not well suited to study hysteresis between the variables. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

14.
One‐dimensional flow simulations were conducted at four locations of the shallow alluvial aquifer of the upper Rhine River (at the Erstein polder) to quantify the time‐dependent moisture distribution, the water flux and the water volume infiltrated in the unsaturated zone as a function of soil heterogeneities during a five‐day‐long flooding event. Three methods of estimating the hydraulic parameters of soil in the vadose zone were tested. They are based on the following: (1) experimental data, (2) soil particle‐size distribution and (3) pedology information on soils. Water fluxes calculated from modelling approaches 2 and 3 were compared with those of the experiment‐based values and the effect of these differences on the arrival time and velocity of water at the water table were analysed. Major differences in water fluxes were found among the methods of estimating the hydrodynamic parameters. At the Terrace location, the groundwater recharge predicted using soil data from methods 1 and 2 are approximately 4500 and 2400 mm, respectively. Flow simulations using soil data and the experiment‐based method show the highest velocities of infiltrating water at the soil surface and largest volume of groundwater infiltration but result in the lowest centres of the moisture content mass. The results obtained using soil data based on the pedological method are similar to those calculated using soil parameters based on the particle‐size distribution of extracted soil samples. Water pressure profiles calculated on Terrace and Channel location, 3 and 7 days after the inundation event agreed reasonably well with those observed when using hydrodynamic parameters from the experiment‐based method. However, the flow model using the pedology‐based parameters largely underestimates the time needed to achieve hydrostatic conditions of the soil water profile once water flooding at the soil surface stops. This can be mainly attributed to the low values of estimated van Genuchten parameter α. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

15.
16.
Han Yu  Craig C. Douglas 《水文研究》2015,29(6):1225-1237
On the basis of unsaturated Darcy's law, the Talbot–Ogden method provides a fast unconditional mass conservative algorithm to simulate groundwater infiltration in various unsaturated soil textures. Unlike advanced reservoir modelling methods that compute unsaturated flow in space, it only discretizes the moisture content domain into a suitable number of bins so that the vertical water movement is estimated piecewise in each bin. The dimensionality of the moisture content domain is extended from one dimensional to two dimensional in this study, which allows us to distinguish pore shapes within the same moisture content range. The vertical movement of water in the extended model imitates the infiltration phase in the Talbot–Ogden method. However, the difference in this extension is the directional redistribution, which represents the horizontal inter‐bin flow and causes the water content distribution to have an effect on infiltration. Using this extension, we mathematically analyse the general relationship between infiltration and the moisture content distribution associated with wetting front depths in different bins. We show that a more negatively skewed moisture content distribution can produce a longer ponding time, whereas a higher overall flux cannot be guaranteed in this situation. It is proven on the basis of the water content probability distribution independent of soil textures. To illustrate this analysis, we also present numerical examples for both fine and coarse soil textures. Copyright © 2014 John Wiley & Sons, Ltd.  相似文献   

17.
Evaporation of soil moisture is one of the most important processes affecting water availability in semiarid ecosystems. Biological soil crusts, which are widely distributed ground cover in these ecosystems, play a recognized role on water processes. Where they roughen surfaces, water residence time and thus infiltration can be greatly enhanced, whereas their ability to clog soil pores or cap the soil surface when wetted can greatly decrease infiltration rate, thus affecting evaporative losses. In this work, we compared evaporation in soils covered by physical crusts, biological crusts in different developmental stages and in the soils underlying the different biological crust types. Our results show that during the time of the highest evaporation (Day 1), there was no difference among any of the crust types or the soils underlying them. On Day 2, when soil moisture was moderately low (11%), evaporation was slightly higher in well‐developed biological soil crusts than in physical or poorly developed biological soil crusts. However, crust removal did not cause significant changes in evaporation compared with the respective soil crust type. These results suggest that the small differences we observed in evaporation among crust types could be caused by differences in the properties of the soil underneath the biological crusts. At low soil moisture (<6%), there was no difference in evaporation among crust types or the underlying soils. Water loss for the complete evaporative cycle (from saturation to dry soil) was similar in both crusted and scraped soils. Therefore, we conclude that for the specific crust and soil types tested, the presence or the type of biological soil crust did not greatly modify evaporation with respect to physical crusts or scraped soils. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

18.
Large-scale simulation of the soil-derived dust emission in semi-arid regions needs to account for the influence of the soil moisture on the wind erosion threshold. Soil water retention consists of molecular adsorption on the soil grain surface and capillary forces between the grain. Interparticle capillary forces (characterized by the moisture tension) are the main factor responsible for the increase of the wind erosion threshold observed when the soil moisture increases. When the soil moisture content is close to but smaller than the maximum amount of adsorbed water, w′ (depending on the soil texture), these capillary forces are considered as not strong enough to significantly increase the erosion threshold. An expression of the moisture tension as a function of soil moisture and w′ is derived from retention curves. From this expression, a parametrization of the ratio of the wet to dry erosion thresholds has been developed as a function of soil moisture and soil texture. The coefficients of this parametrization have been determined by using experimental data from the literature. An empirical relationship between w′ and soil clay content has been established. The erosion threshold ratios simulated for different soil textures were found to be in good agreement with the experimental data.  相似文献   

19.
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.  相似文献   

20.
Soil moisture is widely recognized as a fundamental variable governing the mass and energy fluxes between the land surface and the atmosphere. In this study, the soil moisture modelling at sub‐daily timescale is addressed by using an accurate representation of the infiltration component. For that, the semi‐analytical infiltration model proposed by Corradini et al. (1997) has been incorporated into a soil water balance model to simulate the evolution in time of surface and profile soil moisture. The performances of this new soil moisture model [soil water balance module‐semi‐analytical (SWBM‐SA)] are compared with those of a precedent version [SWBM‐Green–Ampt (GA)] where the GA approach was employed. Their capability to reproduce in situ soil moisture observations at three sites in Italy, Spain and France is analysed. Hourly observations of quality‐checked rainfall, temperature and soil moisture data for a 2‐year period are used for testing the modelling approaches. Specifically, different configurations for the calibration and validation of the models are adopted by varying a single parameter, that is, the saturated hydraulic conductivity. Results indicate that both SWBMs are able to reproduce satisfactorily the hourly soil moisture temporal pattern for the three sites with root mean square errors lower than 0.024 m3/m3 both in the calibration and validation periods. For all sites, the SWBM‐SA model outperforms the SWBM‐GA with an average reduction of the root mean square error of ~20%. Specifically, the higher improvement is observed for the French site for which in situ observations are measured at 30 cm depth, and this is attributed to the capability of the SA infiltration model to simulate the time evolution of the whole soil moisture profile. The reasonable models performance coupled with the need to calibrate only a single parameter makes them useful tools for soil moisture simulation in different regions worldwide, also in scarcely gauged areas. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

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

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