首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 437 毫秒
1.
Abstract

Kanchanapally watershed covering an area of about 11 km2 in Nalgonda district, Andhra Pradesh, India is located in granitic terrain. Groundwater recharge has been estimated from a water balance model using hydrometeorological data from 1978–1994. The monthly recharge estimates obtained from the water balance model formed input for the groundwater flow model during transient model testing. The groundwater flow model has been prepared to simulate steady state groundwater conditions of 1977 using the nested squares finite difference method. The transient groundwater flow model has been tested during 1977–1994 using the estimated recharge values. The present study helped verify the usefulness of monthly recharge estimates for accounting dynamic variations in recharge as reflected in water level fluctuations in hydrographs.  相似文献   

2.
Nitrate transport in the unsaturated zone of a riverbank filtration (RBF) system in Karany, Czech Republic, was studied. Previous study of the system estimated RBF recharge as 60% riverbank filtrate and 40% local groundwater contaminated by nitrates. Nitrate concentrations observed in RBF recently cannot be explained by simple groundwater contamination and a new conception of groundwater recharge is suggested. A two‐component model based on water 18O data modelled recharge of local groundwater. One component of groundwater recharge is rainfall and irrigation water moving through the unsaturated zone of the Quaternary sediments in piston flow. The second component is groundwater from the Cretaceous deposits with a free water table. Both the components of groundwater recharge have different nitrate concentrations, and resulting contamination of groundwater depends on the participation of water from Quaternary and Cretaceous deposits. Nitrates' origins and their mixing in the subsurface were traced by 15N data. Nitrate transport from the unsaturated zone is important and time variable source of groundwater contamination. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

3.
Chahardouly basin is located in the western part of Iran and is characterized by semi‐arid climatic conditions and scarcity in water resources. The main aquifer systems are developed within alluvial deposits. The availability of groundwater is rather erratic owing to the occurrence of hard rock formation and a saline zone in some parts of the area. The aquifer systems of the area show signs of depletion, which have taken place in recent years due to a decline in water levels. Groundwater samples collected from shallow and deep wells were analysed to examine the quality characteristics of groundwater. The major ion chemistry of groundwater is dominated by Ca2+ and HCO3?, while higher values of total dissolved solids (TDS) in groundwater are associated with high concentrations of all major ions. An increase in salinity is recorded in the down‐gradient part of the basin. The occurrence of saline groundwater, as witnessed by the high electrical conductivity (EC), may be attributed to the long residence time of water and the dissolution of minerals, as well as evaporation of rainfall and irrigation return flow. Based on SAR values and sodium content (%Na), salinity appears to be responsible for the poor groundwater quality, rendering most of the samples not suitable for irrigation use. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

4.
Abstract

Reports on the occurrence of fluoride in natural water resources and the associated health hazards due to human consumption have been made from many parts of India during the last decade. With the objective of organizing a systematic scientific programme to understand the behaviour of fluoride in natural water resources in relation to the local hydrogeological and climatic conditions and agricultural use, a typical area constituting the lower Vamsadhara River basin was chosen for a detailed study. High fluoride concentrations in the groundwater reaching a maximum of 3.4 mg 1?1 were observed to be associated with weathered formations of pyroxene amphibolites and pegmatites. The groundwater in the clayey soils contained much less fluoride as compared to the sandy soils. The complex depositional pattern of these sandy and clayey soils plays an important role in the uneven spatial distribution of fluoride in the groundwater. The contribution of fluoride from geological formations is far greater than that from agriculture: the maximum yield of fluoride by superphosphate fertilizer to irrigation water is observed to be 0.34 mg 1?1. The fluoride concentration is expected to be increased in the future as the groundwater is subsaturated with respect to fluorite. An inverse relationship between F and Ca and positive relationships of F with Na, HCO3, PO4 and electrical conductivity were observed. Best relationships were obtained in the fluoride range of 1.0–3.4 mg 1?1.  相似文献   

5.
The effects of surface water flow system changes caused by constructing water‐conservation areas and canals in southeast Florida on groundwater quality under the Atlantic Coastal Ridge was investigated with numerical modeling. Water quality data were used to delineate a zone of groundwater with low total dissolved solids (TDS) within the Biscayne aquifer under the ridge. The delineated zone has the following characteristics. Its location generally coincides with an area where the Biscayne aquifer has high transmissivities, corresponds to a high recharge area of the ridge, and underlies a part of the groundwater mound formed under the ridge prior to completion of the canals. This low TDS groundwater appears to be the result of pre‐development conditions rather than seepage from the canals constructed after the 1950s. Numerical simulation results indicate that the time for low TDS groundwater under the ridge to reach equilibrium with high TDS surface water in the water‐conservation areas and Everglades National Park are approximately 70 and 60 years, respectively. The high TDS groundwater would be restricted to the water‐conservation areas and the park due to its slow eastward movement caused by small hydraulic gradients in Rocky Glades and its mixing with the low TDS groundwater under the high‐recharge area of the ridge. The flow or physical boundary conditions such as high recharge rates or low hydraulic conductivity layers may affect how the spatial distribution of groundwater quality in an aquifer will change when a groundwater flow system reaches equilibrium with an associated surface water flow system.  相似文献   

6.
The deeply buried river‐connected Xishan karst aquifer (XKA) in western Beijing, China, has been suffering from diminishing recharge for several decades, which in turn leads to the disappearing of spring water outflows and continuously lowering of groundwater level in the area. Thus, it is important to correctly recognize the groundwater recharge and flow paths for the sustainable development of the XKA. To investigate these issues, the hydrochemical and isotopic compositions are analysed for both surface water and groundwater samples collected over an area of about 280 km2. Results show that (a) the river water is characterized by high Na contents; (b) the δ2H and δ18O values in the river water are distinctively higher than those of groundwater samples, after experiencing the long‐time evaporative enrichment in the upstream reservoir; (c) the Sr concentrations and 87Sr/86Sr ratios of groundwater clearly indicated the interaction between water and carbonate minerals but excluded the water–silicate interaction; and (d) the groundwater samples in the direct recharge area of the XKA have the lowest Na concentrations and the δ2H and δ18O values. Based on the large differences in the Na contents and 18O values of groundwater and surface water, a simple two‐component mixing model is developed for the study area and the fractions of the river water are estimated for groundwater samples. We find that the distribution pattern of the river water fractions in the XKA clearly shows a change of directions in the preferential flow path of the groundwater from its source zone to the discharge area. Overall, our results suggest that the recharged surface water can be a useful evidence for delineating the groundwater flow path in river‐connected karst aquifer. This study improves our understanding of the heterogeneity in karst groundwater systems.  相似文献   

7.
We evaluated sources and pathways of groundwater recharge for a heterogeneous alluvial aquifer beneath an agricultural field, based on multi‐level monitoring of hydrochemistry and environmental isotopes of a riverside groundwater system at Buyeo, Korea. Two distinct groundwater zones were identified with depth: (1) a shallow oxic groundwater zone, characterized by elevated concentrations of NO3? and (2) a deeper (>10–14 m from the ground surface) sub‐oxic groundwater zone with high concentrations of dissolved Fe, silica, and HCO3?, but little nitrate. The change of redox zones occurred at a depth where the aquifer sediments change from an upper sandy stratum to a silty stratum with mud caps. The δ18O and δ2H values of groundwater were also different between the two zones. Hydrochemical and δ18O? δ2H data of oxic groundwater are similar to those of soil water. This illustrates that recharge of oxic groundwater mainly occurs through direct infiltration of rain and irrigation water in the sandy soil area where vegetable cropping with abundant fertilizer use is predominant. Oxic groundwater is therefore severely contaminated by agrochemical pollutants such as nitrate. In contrast, deeper sub‐oxic groundwater contains only small amounts of dissolved oxygen (DO) and NO3?. The 3H contents and elevated silica concentrations in sub‐oxic groundwater indicate a somewhat longer mean residence time of groundwater within this part of the aquifer. Sub‐oxic groundwater was also characterized by higher δ18O and δ2H values and lower d‐excess values, indicating significant evaporation during recharge. We suggest that recharge of sub‐oxic groundwater occurs in the areas of paddy rice fields where standing irrigation and rain water are affected by strong evaporation, and that reducing conditions develop during subsequent sub‐surface infiltration. This study illustrates the existence of two groundwater bodies with different recharge processes within an alluvial aquifer. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

8.
This paper characterizes a seasonally inundated Danish floodplain wetland in a state close to naturalness and includes an analysis of the major controls on the wetland water and nitrogen balances. The main inputs of water are precipitation and percolation during ponding and unsaturated conditions. Lateral saturated subsurface flow is low. The studied floodplain owes its wetland status to the hydraulic properties of its sediments: the low hydraulic conductivity of a silt–clay deposit on top of the floodplain maintains ponded water during winter, and parts of autumn and spring. A capillary fringe extends to the soil surface, and capillary rise from groundwater during summer maintains near‐saturated conditions in the root zone, and allows a permanently very high evapotranspiration rate. The average for the growing season of 1999 is 3·6 mm day?1 and peak rate is 5·6 mm day?1. In summer, the evapotranspiration is to a large degree supplied by subsurface storage in a confined peat layer underlying the silt–clay. The floodplain sediments are in a very reduced state as indicated by low sulphate concentrations. All nitrate transported into the wetland is thus denitrified. However, owing to modest water exchange with surrounding groundwater and surface water, denitrification is low; 71 kg NO3–N ha?1 during the study period of 1999. Reduction of nitrate diffusing into the sediments during water ponding accounts for 75% of nitrate removal. Biomass production and nitrogen uptake in above‐ground vegetation is high—8·56 t dry matter ha?1 year?1 and 103 kg N ha?1 year?1. Subsurface ammonium concentrations are high, and convective upward transport into the root zone driven by evapotranspiration amounted to 12·8 kg N ha?1year?1. The floodplain wetland sediments have a high nitrogen content, and conditions are very favourable for mineralization. Mineralization thus constitutes 72% of above‐ground plant uptake. The study demonstrates the necessity of identifying controlling factors, and to combine surface flow with vadose and groundwater flow processes in order to fully comprehend the flow and nitrogen dynamics of this type of wetland. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

9.
Abstract

Groundwater is an important water resource and its management is vital for integrated water resources development in semiarid catchments. The River Shiyang catchment in the semiarid area of northwestern China was studied to determine a sustainable multi-objective management plan of water resources. A multi-objective optimization model was developed which incorporated water supplies, groundwater quality, ecology, environment and economics on spatial and temporal scales under various detailed constraints. A calibrated groundwater flow model was supplemented by grey simulation of groundwater quality, thus providing two lines of evidence to use in the multi-objective water management. The response matrix method was used to link the groundwater simulation models and the optimization model. Multi-phase linear programming was used to minimize and compromise the objectives for the multi-period, conjunctive water use optimization model. Based on current water demands, this water use optimization management plan was able to meet ecological, environmental and economic objectives, but did not find a final solution to reduce the overall water deficit within the catchment.  相似文献   

10.
Abstract

Abstract Water balances for a re-vegetated xerophyte shrub (Caragana korshinskii) area were compared to that of a bare surface area by using auto-weighing type lysimeters during the 1990–1995 growing seasons at the southeast Tengger Desert, Shapotou, China. The six-year experiment displayed how major daily water balance components might vary for a bare and a re-vegetated sand dune area. Evapotranspiration from the C. korshinskii lysimeter represented a major part of the water balance. The average annual ET/P ratios varied between 69 and 142%. No seepage was observed for the vegetated lysimeter. For the bare lysimeter, on the other hand, 48 mm or 27% of observed rainfall per year occurred as seepage. These results suggest that re-vegetating large sandy areas with xerophytic shrubs could reduce soil water storage by transpiration. Also, the experimental results indicate that re-vegetating large sandy areas could significantly change groundwater recharge conditions. However, from a viewpoint of desert ecosystem reconstruction, it appears that natural rainfall can sustain xerophytic shrubs such as C. korshinskii which would reduce erosion loss of sand. However, re-vegetation has to be balanced with recharge/groundwater needs of local populations.  相似文献   

11.
ABSTRACT

The temporal dynamics of groundwater–surface water interaction under the impacts of various water abstraction scenarios are presented for hydraulic fracturing in a shale gas and oil play area (23 984.9 km2), Alberta, Canada, using the MIKE-SHE and MIKE-11 models. Water-use data for hydraulic fracturing were obtained for 433 wells drilled in the study area in 2013 and 2014. Modelling results indicate that water abstraction for hydraulic fracturing has very small (<0.35%) negative impacts on mean monthly and annual river and groundwater levels and stream and groundwater flows in the study area, and small (1–4.17%) negative impacts on environmental flows near the water abstraction location during low-flow periods. The impacts on environmental flow depend on the amount of water abstraction and the daily flow over time at a specific river cross-section. The results also indicate a very small (<0.35%) positive impact on mean monthly and annual groundwater contributions to streamflow because of the large study area. The results provide useful information for planning long-term seasonal and annual water abstractions from the river and groundwater for hydraulic fracturing in a large study area.  相似文献   

12.
Stable isotopes, 2Hwater, 18Owater as well as 18Osulphate and 34Ssulphate, were used to study the flow system of shallow groundwater and soil water at the base area of a former leaching heap at the uranium mining area of Ronneburg, Germany. The flow paths and water‐retention times were estimated by comparison of δ2H and δ18O values in groundwater and soil water to the δ2H and δ18O signature of precipitation, giving distinctive inputs of summer or winter precipitation. The points of measuring the groundwater were divided into three categories with different flow conditions: rapid flow, stagnant conditions and a transition zone by hierarchical cluster analysis of δ2H and δ18O values of groundwater. The transit time of groundwater in the rapid flow area is less than 6 months, whereas water in the stagnant zone is stored for at least 1 year. In soil water, a clear response to different input signals is detectable only in the 30‐cm horizon (retention time is about 6 months), whereas at deeper levels a mixing with older water is taking place. The isotopic composition of the dissolved sulphate was used to identify oxidation of sulphides as the source of sulphate. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

13.
Small‐scale heterogeneities and large changes in hydraulic gradient over short distances can create preferential groundwater flow paths that discharge to lakes. A 170 m2 grid within an area of springs and seeps along the shore of Shingobee Lake, Minnesota, was intensively instrumented to characterize groundwater‐lake interaction within underlying organic‐rich soil and sandy glacial sediments. Seepage meters in the lake and piezometer nests, installed at depths of 0·5 and 1·0 m below the ground surface and lakebed, were used to estimate groundwater flow. Statistical analysis of hydraulic conductivity estimated from slug tests indicated a range from 21 to 4·8 × 10?3 m day?1 and small spatial correlation. Although hydraulic gradients are overall upward and toward the lake, surface water that flows onto an area about 2 m onshore results in downward flow and localized recharge. Most flow occurred within 3 m of the shore through more permeable pathways. Seepage meter and Darcy law estimates of groundwater discharge agreed well within error limits. In the small area examined, discharge decreases irregularly with distance into the lake, indicating that sediment heterogeneity plays an important role in the distribution of groundwater discharge. Temperature gradients showed some relationship to discharge, but neither temperature profiles nor specific electrical conductance could provide a more convenient method to map groundwater–lake interaction. These results suggest that site‐specific data may be needed to evaluate local water budget and to protect the water quality and quantity of discharge‐dominated lakes. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

14.
There has been a great deal of research interest regarding changes in flow path/runoff source with increases in catchment area. However, there have been very few quantitative studies taking subscale variability and convergence of flow path/runoff source into account, especially in relation to headwater catchments. This study was performed to elucidate how the contributions and discharge rates of subsurface water (water in the soil layer) and groundwater (water in fractured bedrock) aggregate and change with catchment area increase, and to elucidate whether the spatial variability of the discharge rate of groundwater determines the spatial variability of stream discharge or groundwater contribution. The study area was a 5‐km2 forested headwater catchment in Japan. We measured stream discharge at 113 points and water chemistry at 159 points under base flow conditions. End‐member mixing analysis was used to separate stream water into subsurface water and groundwater. The contributions of both subsurface water and groundwater had large variability below 1 km2. The contribution of subsurface water decreased markedly, while that of groundwater increased markedly, with increases in catchment area. The specific discharge of subsurface water showed a large degree of variability and decreased with catchment area below 0.1 km2, becoming almost constant above 0.1 km2. The specific discharge of groundwater showed large variability below 1 km2 and increased with catchment area. These results indicated that the variabilities of stream discharge and groundwater contribution corresponded well with the variability of the discharge rate of groundwater. However, below 0.1 km2, it was necessary to consider variations in the discharge rates of both subsurface water and groundwater. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

15.
The present study was carried out as part of a complex survey of urban groundwater quality and quantity in Szeged, southeast Hungary. The concentrations of 12 inorganic contaminants in 28 shallow groundwater monitoring wells were determined over a 2-year period (2010–2012). The evaluation of concentrations indicates remarkable contamination all over the city. Discriminant analysis (DA) was used to evaluate the spatial changes of groundwater quality. The groundwater levels were measured over a 14-year period (2000–2013). The fractal properties of water level fluctuations time series characterizing the groundwater system and Tisza River were investigated using rescaled range (R/S) analysis. The resulting Hurst exponents clearly showed the persistency and thus long memory effects of both the groundwater and the river flow. Comparison of the results of DA with the results of R/S analysis thus implies that the geological conditions and the changing groundwater quantities are not related to groundwater quality.
EDITOR A. Castellarin

ASSOCIATE EDITOR A. Fiori  相似文献   

16.
《水文科学杂志》2013,58(6):1230-1240
Abstract

This study aims to identify the hydrogeochemical processes influencing the high fluoride concentrations in groundwater of the Nairobi area, Kenya. For this purpose 16 groundwater samples were collected and analysed. Fluoride concentrations above the WHO standard are found in the downstream areas. The high F? concentrations are correlated with high sodium and pH and low Ca2+ concentrations. Weathering of sodium-rich alkaline igneous rocks causes a pH increase resulting in an increase in HCO3 ? and CO3 2- by dissolution of CO2. Groundwater becomes oversaturated compared to calcite and calcite precipitation occurs, leading to a decrease in Ca2+. This causes a sub-saturation with respect to fluorite and dissolution of fluorite increases the F? concentration. These reactions were modelled using the PHREEQC model and the results showed a good agreement with the measured groundwater quality, indicating that the proposed reactions are plausible for explaining the observed concentrations in groundwater.  相似文献   

17.
Groundwater-surface water interactions (GSI) connect rivers and streams with riparian areas and the adjacent aquifer. Although these interactions exert a substantial control of quantity and quality of both groundwater and surface water, knowledge on GSI along rivers at the regional scale, particularly for inland waterways, is still limited. We investigated GSI along the river Moselle, an important federal inland waterway in Germany, by using radon and tritium to identify gaining (water flux from the aquifer to the surface water) and losing (water flux from the surface water to the aquifer) stream conditions, respectively. Gaining stream conditions were identified by continuously measuring radon along the river during boat surveys with a high spatial resolution (every 2 km) during intermediate (October 2020) and near low flow conditions (August/September 2021). The tritium concentrations in surface water and groundwater and the resulting tritium inventories were used to characterize losing stream conditions Monthly tritium inventories from 2017 to 2022 revealed a mean loss for the whole period of 20.3 % and a mean gain of 21.8%. Both were probably triggered by a combination of losing stream conditions and flood-induced mass transfer of water from the aquifer back into the river as well as discharge fluctuations. At the investigated site Lehmen there were direct indications of an influence of surface water due to elevated tritium concentrations in the groundwater (up to 13.3 Bq L−1). Using radon mass balance modelling, good agreements of simulated versus measured radon data with respect to two groundwater end-member scenarios were obtained during intermediate flow (Spearman's ρ: 0.97 and 0.99; MAE: 10.1 and 3.4 Bq L−1) and near low flow (Spearman's ρ: 0.97 and 0.99; MAE: 11 and 6.5 Bq L−1). Considerable groundwater inflow was limited to the meander of Detzem, where cumulated groundwater inflow of about 19 m3 s−1 (9.5% of total discharge) and 4.2 m3 s−1 (3.8% of total discharge) was simulated during intermediate and near low flow, respectively. However, the groundwater inflow was relatively low compared to alpine streams, for example. The study will help to better identify and quantify GSI at the regional scale and provide methodological guidance for future studies focusing on inland waterways.  相似文献   

18.
Understanding groundwater–surface water (GW–SW) interactions is vital for water management in karstic catchments due to its impact on water quality. The objective of this study was to evaluate and compare the applicability of seven environmental tracers to quantify and localize groundwater exfiltration into a small, human-impacted karstic river system. Tracers were selected based on their emission source to the surface water either as (a) dissolved, predominantly geogenic compounds (radon-222, sulphate and electrical conductivity) or (b) anthropogenic compounds (predominantly) originating from wastewater treatment plant (WWTP) effluents (carbamazepine, tramadol, sodium, chloride). Two contrasting sampling approaches were compared (a) assuming steady-state flow conditions and (b) considering the travel time of the water parcels (Lagrangian sampling) through the catchment to account for diurnal changes in inflow from the WWTP. Spatial variability of the concentrations of all tracers indicated sections of preferential groundwater inflow. Lagrangian sampling techniques seem highly relevant for capturing dynamic concentration patterns of WWTP-derived compounds. Quantification of GW inflow with the finite element model FINIFLUX, based on observed in-stream Rn activities led to plausible fluxes along the investigated river reaches (0.265 m3 s−1), while observations of other natural or anthropogenic environmental tracers produced less plausible water fluxes. Important point sources of groundwater exfiltration can be ascribed to locations where the river crosses geological fault lines. This indicates that commonly applied concepts describing groundwater–surface water interactions assuming diffuse flow in porous media are difficult to transfer to karstic river systems whereas concepts from fractured aquifers may be more applicable. In general, this study helps selecting the best suited hydrological tracer for GW exfiltration and leads to a better understanding of processes controlling groundwater inflow into karstic river systems.  相似文献   

19.
Abstract

Abstract Land development often results in adverse environmental impact for surface and subsurface water systems. For areas close to the coast, land changes may also result in seawater intrusion into coastal aquifers. Due to this, it is important to evaluate potential adverse effects in advance of any land development. For evaluation purposes a combined groundwater recharge model is proposed with a quasi three-dimensional unconfined groundwater flow equation. The catchment water balance for a planned new campus area of Kyushu University in southern Japan, was selected as a case study to test the model approach. Since most of the study area is covered with forest, the proposed groundwater recharge model considers rainfall interception by forest canopy. The results show that simulated groundwater and surface runoff agree well with observations. It is also shown that actual evapotranspiration, including rainfall interception by forest canopy, is well represented in the proposed simulation model. Several hydrological components such as direct surface runoff rate, groundwater spring flow rate to a ground depression, trans-basin groundwater flow etc., were also investigated.  相似文献   

20.
Subsurface flow and heat transport near Freienbrink, NE Germany, was simulated in order to study groundwater–surface water exchange between a floodplains aquifer and a section of the lowland River Spree and an adjacent oxbow. Groundwater exfiltration was the dominant process, and only fast surface water level rises resulted in temporary infiltration into the aquifer. The main groundwater flow paths are identified based on a 3D groundwater flow model. To estimate mass fluxes across the aquifer–surface water interfaces, a 2D flow and heat transport modelling approach along a transect of 12 piezometers was performed. Results of steady‐state and transient water level simulations show an overall high accuracy with a Spearman coefficient ρ = 0.9996 and root mean square error (RMSE) = 0.008 m. Based on small groundwater flow velocities of about 10?7 to 10?6 ms?1, mean groundwater exfiltration rates of 233 l m?2 d?1 are calculated. Short periods of surface water infiltration into the aquifer do not exceed 10 days, and the infiltration rates are in the same range. The heat transport was modelled with slightly less accuracy (ρ = 0.8359 and RMSE = 0.34 °C). In contrast to the predominant groundwater exfiltration, surface water temperatures determine the calculated temperatures in the upper aquifer below both surface water bodies down to 10 m during the whole simulation period. These findings emphasize prevailing of heat conduction over advection in the upper aquifer zones, which seems to be typical for lowland streams with sandy aquifer materials and low hydraulic gradients. Moreover, this study shows the potential of coupled numerical flow and heat transport modelling to understand groundwater–surface water exchange processes in detail. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

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

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