首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
基于FEFLOW的海水入侵数值模拟   总被引:3,自引:1,他引:2  
卢薇  朱照宇  刘卫平 《地下水》2010,32(3):19-21,129
在系统分析珠江口东岸地区地质及水文地质条件的基础上,建立了研究区海水入侵三维溶质模型,利用基于有限元原理的FEFLOW软件对建立的模型进行求解,用地下水位和浓度的动态观测资料对模型进行了识别和校验。运用识别后的模型预测了在不同开采条件下,研究区海水入侵的趋势。研究结果表明减少地下水开采量,海水入侵面积可逐步减少,入侵可得到有效控制;在减少开采量的情况下,相较于基岩裂隙含水层,第四系含水层海水入侵面积收缩较慢。  相似文献   

2.
3.
渗透系数空间变异性研究   总被引:20,自引:6,他引:20       下载免费PDF全文
水文地质参数的空间变异性是随机理论研究的基础,而渗透系数是最为重要的水文地质参数。国外有关渗透参数空间变异性的研究工作已开展很多,但渗透系数究竟服从什么分布目前尚无确切答案。利用Borden含水层试验数据,对渗透系数的空间变异性进行探讨,结果表明若处理方法得当,渗透系数应服从对数正态分布。同时,还对今后野外开展含水层渗透系数空间变异性试验研究提出几点建议。  相似文献   

4.
Sedimentological processes often result in complex three-dimensional subsurface heterogeneity of hydrogeological parameter values. Variogram-based stochastic approaches are often not able to describe heterogeneity in such complex geological environments. This work shows how multiple-point geostatistics can be applied in a realistic hydrogeological application to determine the impact of complex geological heterogeneity on groundwater flow and transport. The approach is applied to a real aquifer in Belgium that exhibits a complex sedimentary heterogeneity and anisotropy. A training image is constructed based on geological and hydrogeological field data. Multiple-point statistics are borrowed from this training image to simulate hydrofacies occurrence, while intrafacies permeability variability is simulated using conventional variogram-based geostatistical methods. The simulated hydraulic conductivity realizations are used as input to a groundwater flow and transport model to investigate the effect of small-scale sedimentary heterogeneity on contaminant plume migration. Results show that small-scale sedimentary heterogeneity has a significant effect on contaminant transport in the studied aquifer. The uncertainty on the spatial facies distribution and intrafacies hydraulic conductivity distribution results in a significant uncertainty on the calculated concentration distribution. Comparison with standard variogram-based techniques shows that multiple-point geostatistics allow better reproduction of irregularly shaped low-permeability clay drapes that influence solute transport.  相似文献   

5.
Numerical models encompassing source zones and receptors, based on representative conceptual models and accounting for aquifer heterogeneity, are needed to understand contaminant migration and fate; however, aquifer characterization seldom provides the necessary data. This study aimed to develop a workflow for field characterization and data integration, which could: (1) be adapted to the definition of subwatershed-scale aquifer heterogeneity (over 10 km2) and (2) adequately support mass transport model development. The study involved the field investigation of a shallow granular aquifer in a 12-km2 subwatershed in Saint-Lambert-de-Lauzon, Canada, in which a decommissioned landfill is emitting a leachate plume managed by natural attenuation. Using proven field methods, the characterization sequence was designed to optimize each method in terms of location, scale of acquisition, density and quality. The emphasis was on the acquisition of detailed indirect geophysical data that were integrated with direct hydraulic and geochemical data. This report focuses on the first qualitative and geostatistical data integration steps of the workflow leading to the development of a hydrogeological conceptual model. This is a prerequisite for further integration steps: prediction of hydrofacies and hydraulic conductivity (K), geostatistical simulations of K, studies of geochemical processes and numerical modeling.  相似文献   

6.
Numerical simulation of sea water intrusion near Beihai, China   总被引:6,自引:0,他引:6  
 A leaky aquifer system occurs in the coastal plain near Beihai, China. Seawater intrusion into the confined aquifer took place along the northern coast. Chloride concentrations at some observation wells increased steadily from 1988 and were at their peak in 1993. A quasi-three-dimensional element model has been developed to simulate the spatial and temporal evolution of hydraulic heads and chloride concentrations of the groundwater near the northern coast. The simulation model was based on the transition zone approach, which requires simultaneous solution of the governing water flow and solute transport equations. An irregular grid of a quadrangle was used to discretize the flow domain. Various aquifer parameters were verified with the numerical model in order to obtain satisfactory matches between computed values and observed data from an investigation. Three pumpage schemes were designed to use the calibrated model for prediction of future changes in water levels and chloride concentrations in groundwater in the study area. Results show that seawater intrusion would worsen in the confined aquifer if the current rates of groundwater pumpage continue. The alternative, to eliminate pumpage in the intruded area and to moderate pumpage rates from water supply wells far from the seashore, may limit seawater intrusion significantly and is considered attractive in the area. Received: 27 September 1999 · Accepted: 27 December 1999  相似文献   

7.
In the framework of various studies to characterize the aquifer at the groundwater experimental field near Montalto Uffugo, Italy, the present work estimates the spatial distribution of hydraulic conductivity of the unconsolidated deposits that underlie the area, by applying the geostatistical technique of kriging with external drift to electrical-resistivity and hydraulic-conductivity data. The reliability of the estimation method was tested by implementing a model, based on the method of cells, that simulates groundwater flow, with the estimated values of hydraulic conductivity. The results obtained indicate that the estimation method used has an acceptable degree of reliability.  相似文献   

8.
Quantitative evaluation of management strategies for long-term supply of safe groundwater for drinking from the Bengal Basin aquifer (India and Bangladesh) requires estimation of the large-scale hydrogeologic properties that control flow. The Basin consists of a stratified, heterogeneous sequence of sediments with aquitards that may separate aquifers locally, but evidence does not support existence of regional confining units. Considered at a large scale, the Basin may be aptly described as a single aquifer with higher horizontal than vertical hydraulic conductivity. Though data are sparse, estimation of regional-scale aquifer properties is possible from three existing data types: hydraulic heads, 14C concentrations, and driller logs. Estimation is carried out with inverse groundwater modeling using measured heads, by model calibration using estimated water ages based on 14C, and by statistical analysis of driller logs. Similar estimates of hydraulic conductivities result from all three data types; a resulting typical value of vertical anisotropy (ratio of horizontal to vertical conductivity) is 104. The vertical anisotropy estimate is supported by simulation of flow through geostatistical fields consistent with driller log data. The high estimated value of vertical anisotropy in hydraulic conductivity indicates that even disconnected aquitards, if numerous, can strongly control the equivalent hydraulic parameters of an aquifer system.  相似文献   

9.
Ongoing hydrogeological research aims to develop a correct management model for the Plio-Pleistocene multi-aquifer system of the Albegna River coastal plain (southern Tuscany, Italy); overexploitation of this aquifer for irrigation and tourism has caused seawater intrusion. The conceptual model is based on field and laboratory data collected during the 1995–2003 period. Meteoric infiltration and flows from the adjoining carbonate aquifer recharge the aquifer. Natural outflow occurs through a diffuse flow into the sea and river; artificial outflow occurs through intensive extraction of groundwater from wells. Water exchanges in the aquifer occur naturally (leakage, closing of aquitard) and artificially (multiscreened wells). The aquifer was represented by a three-dimensional finite element model using the FEFLOW numerical code. The model was calibrated for steady-state and transient conditions by matching computed and measured piezometric levels (February 1995–February 1996). The model helped establish that seawater intrusion is essentially due to withdrawals near the coast during the irrigation season and that it occurs above all in the Osa-Albegna sector, as well as along the river that at times feeds the aquifer. The effects of hypothetical aquifer exploitation were assessed in terms of water budget and hydraulic head evolution.  相似文献   

10.
Korba aquifer is one of the most typical examples of overexploited coastal aquifer in the Mediterranean countries. In fact, from 1985, a considerable piezometric level drop, water salinization, and seawater intrusion were registered in the aquifer. In December 2008, Tunisian authorities initiated a general plan to groundwater management in order to augment groundwater resources, restore the piezometric levels, and improve water quality. The plan consists of artificial recharge of groundwater used treated wastewater through three infiltration basins. During the first 4 years (from December 2008 to December 2012), 1.41 Mm3 of treated wastewater was injected to the Korba aquifer. This study presents a hydrogeological assessment of groundwater evolution during the recharge processes. In this study, 32 piezometric and chemical surveys of 70 piezometers and observed wells are used to present hydrogeological investigation and water quality evolution of wastewater reuse through artificial recharge in Korba coastal aquifer. The piezometric evolution maps are used to specify the positive effect in groundwater level that exceeding 1.5 m in some regions. The interpretation of salinity evolution maps are used to indicate the improving of groundwater quality.  相似文献   

11.
A numerical assessment of seawater intrusion in Gaza, Palestine, has been achieved applying a 3-D variable density groundwater flow model. A two-stage finite difference simulation algorithm was used in steady state and transient models. SEAWAT computer code was used for simulating the spatial and temporal evolution of hydraulic heads and solute concentrations of groundwater. A regular finite difference grid with a 400 m2 cell in the horizontal plane, in addition to a 12-layer model were chosen. The model has been calibrated under steady state and transient conditions. Simulation results indicate that the proposed schemes successfully simulate the intrusion mechanism. Two pumpage schemes were designed to use the calibrated model for prediction of future changes in water levels and solute concentrations in the groundwater for a planning period of 17 years. The results show that seawater intrusion would worsen in the aquifer if the current rates of groundwater pumpage continue. The alternative, to eliminate pumpage in the intruded area, to moderate pumpage rates from water supply wells far from the seashore and to increase the aquifer replenishment by encouraging the implementation of suitable solutions like artificial recharge, may limit significantly seawater intrusion and reduce the current rate of decline of the water levels.  相似文献   

12.
A transient finite difference groundwater flow model has been calibrated for the Nasia sub-catchment of the White Volta Basin. This model has been validated through a stochastic parameter randomization process and used to evaluate the impacts of groundwater abstraction scenarios on resource sustainability in the basin. A total of 1500 equally likely model realizations of the same terrain based on 1500 equally likely combinations of the data of the key aquifer input parameters were calibrated and used for the scenario analysis. This was done to evaluate model non-uniqueness arising from uncertainties in the key aquifer parameters especially hydraulic conductivity and recharge by comparing the realizations and statistically determining the degree to which they differ from each other. Parameter standard deviations, computed from the calibrated data of the key parameters of hydraulic conductivity and recharge, were used as a yardstick for evaluating model non-uniqueness. All model realizations suggest horizontal hydraulic conductivity estimates in the range of 0.03–78.4 m/day, although over 70 % of the area has values in the range of 0.03–14 m/day. Low standard deviations of the horizontal hydraulic conductivity estimates from the 1500 solutions suggest that this range adequately reflects the properties of the material in the terrain. Lateral groundwater inflows and outflows appear to constitute significant components of the groundwater budgets in the terrain, although estimated direct vertical recharge from precipitation amounts to about 7 % of annual precipitation. High potential for groundwater development has been suggested in the simulations, corroborating earlier estimates of groundwater recharge. Simulation of groundwater abstraction scenarios suggests that the domain can sustain abstraction rates of up to 200 % of the current estimated abstraction rates of 12,960 m3/day under the current recharge rates. Decreasing groundwater recharge by 10 % over a 20-year period will not significantly alter the results of this abstraction scenario. However, increasing abstraction rates by 300 % over the period with decreasing recharge by 10 % will lead to drastic drawdowns in the hydraulic head over the entire terrain by up to 6 m and could cause reversals of flow in most parts of the terrain.  相似文献   

13.
Seawater intrusion is a major problem to freshwater resources especially in coastal areas where fresh groundwater is surrounded and could be easily influenced by seawater. This study presents the development of a conceptual and numerical model for the coastal aquifer of Karareis region (Karaburun Peninsula) in the western part of Turkey. The study also presents the interpretation and the analysis of the time series data of groundwater levels recorded by data loggers. The SEAWAT model is used in this study to solve the density-dependent flow field and seawater intrusion in the coastal aquifer that is under excessive pumping particularly during summer months. The model was calibrated using the average values of a 1-year dataset and further verified by the average values of another year. Five potential scenarios were analyzed to understand the effects of pumping and climate change on groundwater levels and the extent of seawater intrusion in the next 10 years. The result of the analysis demonstrated high levels of electrical conductivity and chloride along the coastal part of the study area. As a result of the numerical model, seawater intrusion is simulated to move about 420 m toward the land in the next 10 years under “increased pumping” scenario, while a slight change in water level and TDS concentrations was observed in “climate change” scenario. Results also revealed that a reduction in the pumping rate from Karareis wells will be necessary to protect fresh groundwater from contamination by seawater.  相似文献   

14.
The overexploitation of groundwater in coastal aquifers is often accompanied by seawater intrusion, intensified by climate change and sea level rise. Heading long-term water quality safety and thus the determination of vulnerable zones to seawater intrusion becomes a significant hydrogeological task for many coastal areas. Due to this background, the present study focussed the established methodology of the GIS-based GALDIT model to assess the aquifer vulnerability to seawater intrusion for the Algerian example of the Quaternary coastal Collo aquifer. According to the result analysis overall, more than half of the total surface of the northern study area can be classified as highly vulnerable. Besides the coastline, the areas nearby the local wadis of Guebli and Cherka occur to be the most vulnerable in the region. In view of further map removal performance as well as single-parameter sensitivity analyses from a coupled perspective respectively the GALDIT parameters, distance from the shore (D) and aquifer hydraulic conductivity (A) have been found to be of key significance regarding the model results (mean effective weightings ~?18–19%). Overall, the study results provide a good approximation basis for future management decisions of the Collo aquifer region, including various perspectives such as identification of suitable settings for prospective groundwater pumping wells.  相似文献   

15.
Overexploitation of shallow aquifers on the Syrian coast, north of Latakia (Damsarkho) for irrigation and tourism has caused an intrusion of seawater. The seawater intrusion into this aquifer has been presented by a three-dimensional finite element model using the FEFLOW numerical code. This conceptual model is based on field and laboratory data collected during the period 1966–2003. Meteoric infiltration and flows from the adjoining carbonate aquifer recharge the aquifer; natural outflow occurs through a diffuse flow into the sea; and artificial outflow occurs through intensive extraction of groundwater from wells. Water exchanges in the aquifer occur naturally (leakage) and artificially (multi-screened wells). The model was calibrated for transient conditions. The model helped establish that seawater intrusion is essentially due to withdrawals near the coast during the irrigation season and that it mainly occurs in the Damsarkho plain. The effects of hypothetical aquifer exploitation were assessed in terms of salt budget.  相似文献   

16.
为评价滨海地区海(咸)水入侵区不同类型地层渗透性,提出振荡试验适宜模型选择方法。利用Geoprobe系统对不同平面位置地层实施了气压式振荡试验,部分点位进行了分层测试;利用配线法计算地层渗透系数。结果表明:当试验目标层为承压含水层时,利用Butler模型分析欠阻尼振荡试验,KGS模型分析过阻尼振荡试验,均可获得良好的计算结果;当试验目标层为弱透水层时,KGS模型可用于过阻尼振荡试验的数据分析,但计算结果会受到相邻高渗透性含水层以及钻孔编录不准确的影响。地层渗透性等级未知条件下,需根据地层类型、介质岩性、测试井结构和水位恢复过程曲线特征选择适宜的振荡试验分析模型,以准确获取地层渗透系数值。  相似文献   

17.
基于不同地质统计方法的渗透系数场对污染物运移的影响   总被引:1,自引:0,他引:1  
渗透系数场的空间变异性是影响污染物运移结果的决定因素,而地质统计方法是解决渗透系数空间变异性的主要技术手段。本文利用野外场地实测数据,采用普通克里格法和指示克里格法、顺序高斯模拟法和顺序指示模拟法四种地质统计方法,插值估测和模拟再现随机渗透系数场,进而对比研究四种渗透系数场对大尺度污染物运移的影响。研究结果表明,污染羽的质心位置(一阶矩)主要由渗透系数的平均值来决定;污染羽在空间上的展布范围(二阶矩)主要受渗透系数空间变异方差的影响;条件模拟克服了估计法的平滑效果,较好地再现真实曲线的波动性,渗透系数( lnK)估计方差与污染羽空间二阶矩随着条件模拟次数的增加而减小,并且顺序指示模拟程度更加明显。  相似文献   

18.
Geostatistical estimations of the hydraulic conductivity field (K) in the Carrizo aquifer, Texas, are performed over three regional domains of increasing extent: 1) the domain corresponding to a three-dimensional groundwater flow model previously built (model domain); 2) the area corresponding to the 10 counties encompassing the model domain (County domain), and; 3) the full extension of the Carrizo aquifer within Texas (Texas domain). Two different approaches are used: 1) an indirect approach where transmissivity (T) is estimated first and K is retrieved through division of the T estimate by the screen length of the wells, and; 2) a direct approach where K data are kriged directly. Due to preferential well screen emplacement, and scarcity of sampling in the deeper portions of the formation (> 1 km), the available data set is biased toward high values of hydraulic conductivities. Kriging combined with linear regression, simple kriging with varying local means, kriging with an external drift, and cokriging allow the incorporation of specific capacity as secondary information. Prediction performances (assessed through cross-validation) differ according to the chosen approach, the considered variable (log-transformed or back-transformed), and the scale of interest. For the indirect approach, kriging of log T with varying local means yields the best estimates for both log-transformed and back-transformed variables in the model domain. For larger regional scales (County and Texas domains), cokriging performs generally better than other kriging procedures when estimating both (log T) and T. Among procedures using the direct approach, the best prediction performances are obtained using kriging of log K with an external drift. Overall, geostatistical estimation of the hydraulic conductivity field at regional scales is rendered difficult by both preferential well location and preferential emplacement of well screens in the most productive portions of the aquifer. Such bias creates unrealistic hydraulic conductivity values, in particular, in sparsely sampled areas.  相似文献   

19.
Management of groundwater resources can be improved by using groundwater models to perform risk analyses and to improve development strategies, but a lack of extensive basic data often limits the implementation of sophisticated models. Dar es Salaam in Tanzania is an example of a city where increasing groundwater use in a Pleistocene aquifer is causing groundwater-related problems such as saline intrusion along the coastline, lowering of water-table levels, and contamination of pumping wells. The lack of a water-level monitoring network introduces a problem for basic data collection and model calibration and validation. As a replacement, local water-supply wells were used for measuring groundwater depth, and well-top heights were estimated from a regional digital elevation model to recalculate water depths to hydraulic heads. These were used to draw a regional piezometric map. Hydraulic parameters were estimated from short-time pumping tests in the local wells, but variation in hydraulic conductivity was attributed to uncertainty in well characteristics (information often unavailable) and not to aquifer heterogeneity. A MODFLOW model was calibrated with a homogeneous hydraulic conductivity field and a sensitivity analysis between the conductivity and aquifer recharge showed that average annual recharge will likely be in the range 80–100 mm/year.  相似文献   

20.
多尺度非均质多孔介质中溶质运移的蒙特卡罗模拟   总被引:4,自引:0,他引:4       下载免费PDF全文
探讨了将蒙特卡罗(Monte Carlo)方法应用于多尺度非均质含水层中溶质运移模拟的方法。所研究的含水层由两种具有不同渗透系数统计特征的多孔介质所组成,每一种多孔介质是非均质的,且其渗透系数场符合平稳假设,而整个模拟区的渗透系数是非平稳的。Monte Carlo方法要求参数是平稳的,因此,分别对两种多孔介质产生若干随机渗透系数场后,用两种方法进行组合,并进行溶质运移的模拟计算。通过对计算结果的分析,综合考虑计算精度、计算时间等因素,得出了处理多尺度非均质多孔介质中溶质运移问题的较好方法。  相似文献   

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

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