首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
In this study, we attempted to analyse a drawdown pattern around a pumping well in an unconfined sandy gravelly aquifer constructed in a laboratory tank by means of both experimental and numerical modelling of groundwater flow. The physical model consisted of recharge, aquifer and discharge zones. Permeability and specific yield of the aquifer material were determined by Dupuit approximation under steady‐state flow and stepwise gravitational drainage of groundwater, respectively. The drawdown of water table in pumping and neighbouring observation wells was monitored to investigate the effect of no‐flow boundary on the drawdown pattern during pumping for three different boundary conditions: (i) no recharge and no discharge with four no‐flow boundaries (Case 1); (ii) no recharge and reservoir with three no‐flow boundaries (Case 2); (iii) recharge and discharge with two no‐flow boundaries (Case 3). Based on the aquifer parameters, numerical modelling was also performed to compare the simulated drawdown with that observed. Results showed that a large difference existed between the simulated drawdown and that observed in wells for all cases. The reason for the difference could be explained by the formation of a curvilinear type water table between wells rather than a linear one due to a delayed response of water table in the capillary fringe. This phenomenon was also investigated from a mass balance study on the pumping volume. The curvilinear type of water table was further evidenced by measurement of water contents at several positions in the aquifer between wells using time domain reflectometry (TDR). This indicates that the existing groundwater flow model applicable to an unconfined aquifer lacks the capacity to describe a slow response of water table in the aquifer and care should be taken in the interpretation of water table formation in the aquifer during pumping. Copyright © 2001 John Wiley & Sons, Ltd.  相似文献   

2.
This study used a two‐dimensional steady‐state finite‐element groundwater flow model to simulate groundwater flow in two Newfoundland blanket peat complexes and to examine flow system sensitivity to changes in water table recharge and aquifer properties. The modelling results were examined within the context of peat‐forming processes in the two complexes. Modelled flow compared favourably with observed flow. The sensitivity analyses suggested that more highly decomposed bog peat along bog margins probably has/had a positive impact on net peat accumulation within bog interiors. Peat with lower hydraulic conductivity along bog margins effectively impedes lateral drainage, localizes water table drawdown to extreme bog margins, and elevates water tables along bog interiors. Peat formation and elevated water tables in adjacent poor fens/laggs currently rely on placic and ortstein horizons impeding vertical drainage and water flow inputs from adjacent bogs. Modest reductions in atmospheric recharge were found to govern bog‐flow‐system geometries in a way that would adversely affect paludification processes in adjacent fens/laggs. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

3.
This study addresses the influence of landslide dams on surface water drainage and groundwater flow. In the study area of Scanno Lake and Sagittario River (Central Italy), a limestone rockslide‐avalanche formed a lake, which has an outlet that is occasionally active, showing infiltration into the rockslide dam. Several springs are present at the lake's base and are partly fed by seepage through the rockslide debris. Piezometric surveys, discharge measurements, pumping tests and chemical analyses are tools used to build a conceptual model of the groundwater flow and to evaluate the flow through the rockslide debris. Seasonal water isotopic signatures validate the assumed model, showing a mixing of infiltration recharge and groundwater seepage throughout the rockslide debris. Various recharge areas have been found for springs, pointing out those directly fed by the rockslide debris aquifer. Hypotheses about seasonal groundwater mixing between the regional carbonate aquifer and the rockslide debris aquifer are supported by isotope results. Seasonal changes in groundwater table level due to recharge and surface losses from seasonal outlet have been correlated with isotopic groundwater composition from the rockslide debris aquifer and the downstream springs; this relationship highlights the role of the rockslide dam body on the hydrodynamics of the studied area. Relationships between surface waters and groundwater in the area have been completely understood on the basis of water isotopic fingerprinting, finally obtaining a complete evaluation of groundwater renewable resources and its regimen. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

4.
A chalky aquifer located near Liege (Belgium), is recharged by infiltration through the overlying loess. Wells and collecting tunnels produce a daily flow of 60 000 m3 out of this aquifer when hydrogeological balances have shown that an average yield of 100 000 m3/day should be possible.Finite element modelling has been developed to foresee the evolution of the water table to get some additional information especially about the main drainage axis.The transient flow constitutive laws are recalled in confined and unconfined aquifers. A new law is proposed to model the water table surface in transient conditions and with a fixed meshing network. Using the local flow equilibrium and the virtual power principle, the FEM formulation is set up. The time integration and the iteration technique are shortly discussed.The 3D discretization and the modelling of the entire aquifer has been realized. The problem requires about 3600 DOF and 2670 8-nodes isoparametric brick finite elements. The modelling has been quite delicate because of the geometric complexity of the different geological layers. This complexity justifies fully the use of the finite element method; there is indeed a great diversity of geological characteristics and the range of the different permeabilities is very wide. Sensitivity of the model to permeabilities and storage coefficient variations has been studied and various numerical problems have been notified.The calibration procedure is described in its main steps, and the most significant results are presented.  相似文献   

5.
The water budget in clay shale terrain is controlled by a complex interaction between the vertisol soil layer, the underlying fractured rock, land use, topography, and seasonal trends in rainfall and evapotranspiration. Rainfall, runoff, lateral flow, soil moisture, and groundwater levels were monitored over an annual recharge cycle. Four phases of soil–aquifer response were noted over the study period: (1) dry‐season cracking of soils; (2) runoff initiation, lateral flow and aquifer recharge; (3) crack closure and down‐slope movement of subsurface water, with surface seepage; (4) a drying phase. Surface flow predominated within the watershed (25% of rainfall), but lateral flow through the soil zone continued for most of the year and contributed 11% of stream flow through surface seepage. Actual flow through the fractured shale makes up a small fraction of the water budget but does appear to influence surface seepage by its effect on valley‐bottom storage. When the valley soil storage is full, lateral flow exits onto the valley‐bottom surface as seasonal seeps. Well response varied with depth and hillslope position. FLOWTUBE model results and regional recharge estimates are consistent with an aquifer recharge of 1·6% of annual precipitation calculated from well heights and specific yield of the shale aquifer. Copyright © 2005 John Wiley & Sons, Ltd.  相似文献   

6.
One of the most important issues for water resource management is developing strategies for groundwater modelling that are adaptable to data scarcity. These strategies are particularly important in arid and semi‐arid areas where access to data is poor and data collection is difficult, such as the Lake Chad Basin in Africa. In the present study, we establish a numerical groundwater flow model and evaluate the effects of dry and wet periods on groundwater recharge in the Chari–Logone area (96 000 km2) of the Lake Chad Basin. Boundary conditions, flow direction, sources, and sinks for the Chari–Logone local model were obtained by revising and remodelling the Lake Chad Basin regional hydrogeological model (508 400 km2) developed by the BRGM (Bureau de Recherches Géologiques et Minières) in the 1990s. The simulated aquifer water level showed good agreement with observed levels. Aquifer recharge is primarily determined by river–aquifer interactions and mostly occurs in the southern section of the study area. In wet years, groundwater recharge also occurs in the N'Djamena area. The approach we adopted provided relevant results and was useful as an initial step in more detailed modelling of the area. It also proved to be a useful method for groundwater modelling in large semi‐arid and arid regions where available data are scarce. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

7.
A conceptual model of eolian transport is proposed to address the widely distributed, high concentrations of hexavalent chromium (Cr+6) observed in ground water in the Emirate of Abu Dhabi, United Arab Emirates. Concentrations (30 to more than 1000 μg/L Cr+6) extend over thousands of square kilometers of ground water systems. It is hypothesized that the Cr is derived from weathering of chromium-rich pyroxenes and olivines present in ophiolite sequence of the adjacent Oman (Hajar) Mountains. Cr+3 in the minerals is oxidized to Cr+6 by reduction of manganese and is subsequently sorbed on iron and manganese oxide coatings of particles. When the surfaces of these particles are abraded in this arid environment, they release fine, micrometer-sized, coated particles that are easily transported over large distances by wind and subsequently deposited on the surface. During ground water recharge events, the readily soluble Cr+6 is mobilized by rain water and transported by advective flow into the underlying aquifer. Chromium analyses of ground water, rain, dust, and surface (soil) deposits are consistent with this model, as are electron probe analyses of clasts derived from the eroding Oman ophiolite sequence. Ground water recharge flux is proposed to exercise some control over Cr+6 concentration in the aquifer.  相似文献   

8.
ABSTRACT

The estimation of recharge and boundary flux is an important problem in deterministic groundwater modelling since these quantities are often difficult to measure directly in the field. A new method (inferred recharge) is proposed for this estimation, based upon water level and transmissivity observations; it is particularly applicable to those semiarid regions where a steady state formulation can be used and recharge is confined to certain known areas. The inferred recharge method which assumes that recharge occurs in these areas only is compared with an alternative method which does not use this information, and is found to be superior. In the example of the inferred recharge method applied to the Oman coastal plain, random normal errors have been added to both water levels and transmissivities in simulation experiments to assess the effect of errors in the data. The results of the simulations have been used to test the reliability of variance estimates derived from the theory. In addition, statistical tests have been used to examine the inferred recharge results.  相似文献   

9.
Xi Chen  Xunhong Chen   《Journal of Hydrology》2003,280(1-4):246-264
During a flood period, stream-stage increases induce infiltration of stream water into an aquifer; subsequent declines in stream stage cause a reverse motion of the infiltrated water. This paper presents the results of the water exchange rate between a stream and aquifer, the storage volume of the infiltrated stream water in the surrounding aquifer (bank storage), and the storage zone. The storage zone is the part of aquifer where groundwater is replaced by stream water during the flood. MODFLOW was used to simulate stream–aquifer interactions and to quantify rates of stream infiltration and return flow. MODPATH was used to trace the pathlines of the infiltrated stream water and to determine the size of the storage zone. Simulations were focused on the analyses of the effects of the stream-stage fluctuation, aquifer properties, the hydraulic conductivity of streambed sediments, regional hydraulic gradients, and recharge and evapotranspiration (ET) rates on stream–aquifer interactions. Generally, for a given stream–aquifer system, larger flow rates result from larger stream-stage fluctuations; larger storage volumes and storage zones are produced by larger and longer-lasting fluctuations. For a given stream-stage hydrograph, a lower-permeable streambed, an aquitard, or an anisotropic aquifer of low vertical hydraulic conductivity can significantly reduce the rate of infiltration and limit the size of the storage zone. The bank storage solely caused by the stage fluctuation differs slightly between gaining and losing streams. Short-term rainfall recharge and ET loss in the shallow groundwater slightly influence on the flow rate, but their effects on bank storage in a larger area for a longer period can be considerable.  相似文献   

10.
On the basis of one-dimensional theoretical water flow model, we demonstrate that the groundwater level variation follows a pattern similar to recharge fluctuation, with a time delay that depends on the characteristics of aquifer, recharge pattern as well as the distance between the recharge and observation locations. On the basis of a water budget model and the groundwater flow model, we propose an empirical model that links climatic variables to groundwater level. The empirical model is tested using a partial data set from historical records of water levels from more than 80 wells in a monitoring network for the carbonate rock aquifer, southern Manitoba, Canada. The testing results show that the predicted groundwater levels are very close to the observed ones in most cases. The overall average correlation coefficient between the predicted and observed water levels is 0.92. This proposed empirical statistical model could be used to predict variations in groundwater level in response to different climate scenarios in a climate change impact assessment.  相似文献   

11.
Storage of water in aquifers using injection wells is an efficient way for utilizing excess desalinated water in arid regions. In this investigation we estimate the benefits of optimally recharging seasonal surplus desalinated water into a strategic coastal aquifer already benefitting from natural recharge of flash-floods water by a recharge dam. Since, usually the buyers of desalinated water commit to purchase surplus desalinated water under take-or-pay contracts, any attempt in utilizing the paid water is beneficial. Coastal cities are observing an increased urbanization leaving limited space for aquifer recharge infrastructure. In order to determine the optimal location of wells and maximize the use of surplus desalinated water available in winter period, a decision tool combining a numerical groundwater flow simulation model (MODFLOW) with an optimization model is developed. The results of this study show that increasing the number of wells from the existing 45 wells to 173 would allow storing 31.4 million cubic meter per year of excess desalinated water into the aquifer that can be used during later during summer months. The net benefit would reach US$55 million/year while the cost of drilling the new wells is US$5.11 million.  相似文献   

12.
We infer information about the nature of groundwater flow within a karst aquifer from the physical and chemical response of a spring to storm events. The spring discharges from the Maynardville Limestone in Bear Creek Valley, Tennessee. Initially, spring discharge peaks approximately 1–2 h from the midpoint of summer storms. The initial peak is likely due to surface loading, which pressurizes the aquifer and results in water moving out of storage. All of the storms monitored exhibited recessions that follow a master recession curve very closely, indicating that storm response is fairly consistent and repeatable, independent of the time between storms and the configuration of the rain event itself. Electrical conductivity initially increases for 0.5–2.9 days (longer for smaller storms), the result of moving older water out of storage. This is followed by a 2.1–2.5 day decrease in conductivity, resulting from an increasing portion of low conductivity recharge water entering the spring. Stable carbon isotope data and the calcite saturation index of the spring water also support this conceptual model. Spring flow is likely controlled by displaced water from the aquifer rather than by direct recharge through the soil zone.  相似文献   

13.
Regional estimates of aquifer recharge are needed in data-scarce regions such as the Basin of Mexico, where nearly 20 million people are located and where the Basin’s aquifer system represents the main water source. In order to develop the spatio-temporal estimates of aquifer recharge and to analyze to what extent urban growth has affected aquifer recharge, this work presents a daily soil water balance which uses different vegetation and soil types as well as the effect of topography on climatological variables and evapotranspiration. The soil water balance was applied on a daily time step in the Basin of Mexico for the period 1975–1986, obtaining an annually-lumped potential recharge flow of 10.9–23.8 m3/s (35.9–78.1 mm) in the entire Basin, while the monthly values for the year with the largest lumped recharge value (1981 = 78.1 mm) range from 1 m3/s (0.3 mm) in December to 87.9 m3/s (23.7 mm) in June. As aquifer recharge in the Basin mainly occurs by subsurface flow from its enclosing mountains as Mountain Block Recharge, urban growth has had a minimal impact on aquifer recharge, although it has diminished recharge in the alluvial plain.  相似文献   

14.
Intermontane basin aquifers worldwide, particularly in the Himalayan region, are recharged largely by the adjoining mountains. Recharge in these basins can occur either by water infiltrating from streams near mountain fronts (MFs) as mountain front recharge (MFR) or by sub-surface mountain block infiltration as mountain block recharge (MBR). MFR and MBR recharge are challenging to distinguish and are least quantified, considering the lack of extensive understanding of the hydrological processes in the mountains. This study used oxygen and hydrogen isotopes (δ18O and δ2H), electrical conductivity (EC) data, hydraulic head, and groundwater level data to differentiate MFR and MBR. Groundwater level data provide information about the groundwater-surface water interactions and groundwater flow directions, whereas isotopes and EC data are used to distinguish and quantify different recharge sources. The present methodology is tested in an intermontane basin of the Himalayan region. The results suggest that karst springs (KS) and deep groundwater (DGW) recharge are dominated by snowmelt (47% ± 10% and 46% ± 9%) as MBR from adjacent mountains, insignificantly affected by evaporation. The hydraulic head data and isotopes indicate Quaternary shallow groundwater (SGW) aquifer system recharge as MFR of local meteoric water with significant evaporation. The results indicate several flow paths in the aquifer system, a local flow for KS, intermediate flow for SGW, and regional flow for DGW. The findings will significantly impact water resource management in the area and provide vital baseline knowledge for sustainable groundwater management in other Himalayan intermontane basins.  相似文献   

15.
Seawater intrusion into fresh groundwater formations generally results inadvertently from human activities, such as over‐abstraction from coastal aquifers. This article describes the data analysis to quantify drain–aquifer interactions in a low‐lying pump‐drained coastal aquifer, which is subject to saline intrusion due to widespread land drainage, and the resulting development and application of a numerical groundwater model to understand the spatial groundwater system behaviour (including groundwater salinity fluxes). Without measured flow data in this pump‐drained catchment, a novel groundwater head‐dependent approach to hydrograph separation is described. Time‐variant and time‐invariant MODFLOW analyses are utilised to examine the flow processes. A new approach to calculate drain coefficients, which represent the extensive network of drainage ditches in the regional model, using field information, is described; the sum of the drainage coefficients are close to the values independently estimated from the head‐dependent hydrograph separation. Results show that (1) the groundwater flows into the drainage systems are well reproduced using the new drain coefficients, (2) particle tracking of fresh and saline water can explain observed spatial salinity distribution within drainage networks and (3) the modelled flow of seawater across the coast is approximately 25% greater than that discharged by the pumps, demonstrating the need for drainage management to be aware of the slow response of groundwater systems to past drainage system changes. The article demonstrates that numerical groundwater modelling can produce the improved understanding needed to inform management decisions in such complex environments. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

16.
Water table fluctuation (WTF) methods are a primary and well-established way to determine groundwater recharge based on the direct response of the water table to precipitation input. An emerging complexity of recharge is whether it occurs as an episodic and transient process, or a continuous steady-state process, however, most studies have not focused on these short-term vs. long-term timescales, in part because of a lack of data resolution. Here, high-resolution (subhourly) precipitation and water-level data are analyzed for wells in the suburbs of New York City using two contrasting WTF approaches, with a common mathematical basis, that are suited to episodic and continuous processes. The resulting hourly recharge results, like the individual water-level records from comparable wells, are sensitive indicators of subtle differences in aquifer conditions such as thickness of the unsaturated zone, position in the flow system and localized preferential flow. While the episodic, transient approach excludes diffuse recharge by design, the continuous, steady-state approach is influenced by short-term precipitation events, and therefore integrates transient processes to some extent. However, the continuous, steady-state approach is subject to its own limitations relating to position in the aquifer system, and may overestimate recharge if aquifer conditions are not well understood. More widespread use of higher resolution data as well as understanding aquifer conditions and refining aquifer parameters would improve WTF recharge estimation.  相似文献   

17.
To improve the knowledge of the regionally important Continental Terminal 3 (CT3) aquifer in south-western Niger, fifteen magnetic resonance soundings (MRS) were carried out in December 2005 in the vicinity of wells and boreholes. The output MRS geophysical parameters, i.e. water content and decay constants versus depth, were compared to hydrogeological characteristics, i.e. water table depth, total porosity, specific yield and transmissivity estimated from direct measurements, pumping tests and transient groundwater modelling. The MRS-determined parameters were then used to estimate the rates of groundwater recharge.Contained in poorly consolidated Tertiary sandstones, the CT3 aquifer's water table has continuously risen by 4 m in total over the past four decades. Additionally, a significant portion of this increase has occurred in the past decade alone, with an annual rise now ranging between 0.1 and 0.3 m depending on the monitored well. Increase in groundwater recharge due to land clearance and deforestation explains this situation. According to previous estimations, the pre-clearing recharge ranged from 1 to 5 mm per year in 1950–60 s, while more recent recharge rates (1990s–2000s) range from 20 to 50 mm per year. These recharge values are directly affected by estimated aquifer specific yield value, while the spatial variation of rates of water table rise can be attributed to large scale hydrodynamic heterogeneities in the aquifer. However, few field measurements were available to confirm these assumptions.The main results of this study are: (1) The water table depth and aquifer transmissivity are estimated from MRS output parameters with an average accuracy of ± 10% and ± 9% respectively. (2) The MRS-determined water content is linked to both the total porosity and the specific yield of the aquifer, but no quantitative formulation can be proposed as yet. (3) Using the average MRS-determined water content over the investigated area, i.e. 13%, the groundwater recharge rates can be estimated to be ~ 2 mm per year in the 1950–1960s (pre-clearing period), and ~ 23 mm per year for the last decade. (4) The variations in specific yield and transmissivity cannot explain by themselves the spatial variability of the rise of the water table. (5) The ranges in transmissivity and water content obtained from MRS are more realistic than the groundwater modelling outputs. Therefore, MRS could be used to better constrain the aquifer parameters in groundwater modelling with a dense site network.Finally, this work illustrates how MRS can successfully improve characterisation and transient multi-year groundwater balance of commonly found sedimentary aquifers, particularly when integrated with well observations and pumping tests.  相似文献   

18.
Patterns and Age Distribution of Ground-Water Flow to Streams   总被引:2,自引:0,他引:2  
Simulations of ground-water flow in a generic aquifer system were made to characterize the topology of ground-water flow in the stream subsystem and to evaluate its relation to deeper ground-water flow. The flow models are patterned after hydraulic characteristics of aquifers of the Atlantic Coastal Plain and are based on numerical solutions to three-dimensional, steady-state, unconfined flow. The models were used to evaluate the effects of aquifer horizontal-to-vertical hydraulic conductivity ratios, aquifer thickness, and areal recharge rates on flow in the stream subsystem. A particle tracker was used to determine flow paths in a stream subsystem, to establish the relation between ground-water seepage to points along a simulated stream and its source area of flow, and to determine ground-water residence time in stream subsystems. In a geometrically simple aquifer system with accretion, the source area of flow to streams resembles an elongated ellipse that tapers in the down gradient direction. Increased recharge causes an expansion of the stream subsystem. The source area of flow to the stream expands predominantly toward the stream headwaters. Base flow gain is also increased along the reach of the stream. A thin aquifer restricts ground-water flow and causes the source area of flow to expand near stream headwaters and also shifts the start-of-flow to the drainage basin divide. Increased aquifer anisotropy causes a lateral expansion of the source area of flow to streams. Ground-water seepage to the stream channel originates both from near- and far-recharge locations. The range in the lengths of flow paths that terminate at a point on a stream increase in the downstream direction. Consequently, the age distribution of ground water that seeps into the stream is skewed progressively older with distance downstream. Base flow ia an integration of ground water with varying age and potentially different water quality, depending on the source within the drainage basin. The quantitative results presented indicate that this integration can have a wide and complex residence time range and source distribution.  相似文献   

19.
Water budget analyses are important for the evaluation of the water resources in semiarid and arid regions. The lack of observed data is the major obstacle for hydrological modelling in arid regions. The aim of this study is the analysis and calculation of the natural water resources of the Western Dead Sea subsurface catchment, one which is highly sensitive to rainfall resulting in highly variable temporal and spatial groundwater recharge. We focus on the subsurface catchment and subsequently apply the findings to a large‐scale groundwater flow model to estimate the groundwater discharge to the Dead Sea. We apply a semidistributed hydrological model (J2000g), originally developed for the Mediterranean, to the hyperarid region of the Western Dead Sea catchment, where runoff data and meteorological records are sparsely available. The challenge is to simulate the water budget, where the localized nature of extreme rainstorms together with sparse runoff data results in few observed runoff and recharge events. To overcome the scarcity of climate input data, we enhance the database with mean monthly rainfall data. The rainfall data of 2 satellites are shown to be unsuitable to fill the missing rainfall data due to underrepresentation of the steep hydrological gradient and temporal resolution. Hydrological models need to be calibrated against measured values; hence, the absence of adequate data can be problematic. Therefore, our calibration approach is based on a nested strategy of diverse observations. We calculate a direct surface runoff of the Western Dead Sea surface area (1,801 km2) of 3.4 mm/a and an average recharge (36.7 mm/a) for the 3,816 km2 subsurface drainage basin of the Cretaceous aquifer system.  相似文献   

20.
Across equatorial Africa, increasing demand for groundwater has raised concerns about resource sustainability and has highlighted the need for reliable estimates of groundwater recharge. Recharge investigations in this environment are typically inhibited by a shortage of good quality meteorological and hydrogeological records. Moreover, when recharge studies are attempted they tend to rely on a single technique and frequently lack corroborating evidence to substantiate recharge predictions. In recent studies undertaken in the Aroca catchment of the Victoria Nile basin in central Uganda, the timing and magnitude of recharge determined by a soil moisture balance approach are supported by stable isotope data and groundwater flow modelling. The soil moisture balance study reveals that recharge averages in the order of 200 mm year−1 and is more dependent on the number of heavy (more than 10 turn day−1) rainfall events than the total annual volume of rainfall. Stable isotope data suggest independently that recharge occurs during the heaviest rains of the monsoons, and further establish that recharge stems entirely from the direct infiltration of rainfall, an assumption implicit in the soil moisture balance approach. Deforestation over the last 30 years is shown to have more than doubled the recharge estimate. Aquifer flow modelling supports the recharge estimates but demonstrates that the vast majority (over 99%) of recharging waters must be transmitted by the aquifer in the regolith rather the underlying bedrock fractures which have traditionally been developed for rural water supplies.  相似文献   

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

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