首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 78 毫秒
1.
Quantifying groundwater flow at seepage faces is crucial because seepage faces influence the hydroecology and water budgets of watersheds, lakes, rivers and oceans, and because measuring groundwater fluxes directly in aquifers is extremely difficult. Seepage faces provide a direct and measurable groundwater flux but there is no existing method to quantitatively image groundwater processes at this boundary. Our objective is to determine the possibilities and limitations of thermal imagery in quantifying groundwater discharge from discrete seeps. We developed a conceptual model of temperature below discrete seeps, observed 20 seeps spectacularly exposed in three dimensions at an unused limestone quarry and conducted field experiments to examine the role of diurnal changes and rock face heterogeneity on thermal imagery. The conceptual model suggests that convective air‐water heat exchange driven by temperature differences is the dominant heat transfer mechanism. Thermal imagery is effective at locating and characterizing the flux of groundwater seeps. Areas of active groundwater flow and ice growth can be identified from thermal images in the winter, and seepage rates can be differentiated in the summer. However, the application of thermal imagery is limited by diverse factors including technical issues of image acquisition, diurnal changes in radiation and temperature, and rock face heterogeneity. Groundwater discharge rates could not be directly quantified from thermal imagery using our observations but our conceptual model and experiments suggest that thermal imagery could quantify groundwater discharge when there are large temperature differences, simple cliff faces, non‐freezing conditions, and no solar radiation.  相似文献   

2.
In cases when an equivalent porous medium assumption is suitable for simulating groundwater flow in bedrock aquifers, estimation of seepage into underground mine workings (UMWs) can be achieved by specifying MODFLOW drain nodes at the contact between water bearing rock and dewatered mine openings. However, this approach results in significant numerical problems when applied to simulate seepage into an extensive network of UMWs, which often exist at the mine sites. Numerical simulations conducted for individual UMWs, such as a vertical shaft or a horizontal drift, showed that accurate prediction of seepage rates can be achieved by either applying grid spacing that is much finer than the diameter/width of the simulated openings (explicit modeling) or using coarser grid with cell sizes exceeding the characteristic width of shafts or drifts by a factor of 3. Theoretical insight into this phenomenon is presented, based on the so-called well-index theory. It is demonstrated that applying this theory allows to minimize numerical errors associated with MODFLOW simulation of seepage into UMWs on a relatively coarse Cartesian grid. Presented examples include simulated steady-state groundwater flow from homogeneous, heterogeneous, and/or anisotropic rock into a vertical shaft, a horizontal drift/cross-cut, a ramp, two parallel drifts, and a combined system of a vertical shaft connected to a horizontal drift.  相似文献   

3.
The quantification of groundwater recharge is an important but challenging task in groundwater flow modeling because recharge varies spatially and temporally. The goal of this study is to present an innovative methodology to estimate groundwater recharge rates and zone structures for regional groundwater flow models. Here, the unknown recharge field is partitioned into a number of zones using Voronoi Tessellation (VT). The identified zone structure with the recharge rates is associated through a simulation‐optimization model that couples MODFLOW‐2000 and the hybrid PSOLVER optimization algorithm. Applicability of this procedure is tested on a previously developed groundwater flow model of the Tahtal? Watershed. Successive zone structure solutions are obtained in an additive manner and penalty functions are used in the procedure to obtain realistic and plausible solutions. One of these functions constrains the optimization by forcing the sum of recharge rates for the grid cells that coincide with the Tahtal? Watershed area to be equal to the areal recharge rate determined in the previous modeling by a separate precipitation‐runoff model. As a result, a six‐zone structure is selected as the best zone structure that represents the areal recharge distribution. Comparison to results of a previous model for the same study area reveals that the proposed procedure significantly improves model performance with respect to calibration statistics. The proposed identification procedure can be thought of as an effective way to determine the recharge zone structure for groundwater flow models, in particular for situations where tangible information about groundwater recharge distribution does not exist.  相似文献   

4.
Variations in lake seepage were studied along a 130 m shoreline of Mirror Lake NH. Seepage was downward from the lake to groundwater; rates measured from 28 seepage meters varied from 0 to ?282 cm/d. Causes of this variation were investigated using electrical resistivity surveys and lakebed sediment characterization. Two‐dimensional (2D) resistivity surveys showed a transition in lakebed sediments from outwash to till that correlated with high‐ and low‐seepage zones, respectively. However, the 2D survey was not able to predict smaller scale variations within these facies. In the outwash, fast seepage was associated with permeability variations in a thin (2 cm) layer of sediments at the top of the lakebed. In the till, where seepage was slower than that in the outwash, a three‐dimensional resistivity survey mapped a point of high seepage associated with heterogeneity (lower resistivity and likely higher permeability). Points of focused flow across the sediment–water interface are difficult to detect and can transmit a large percentage of total exchange. Using a series of electrical resistivity geophysical methods in combination with hydrologic data to locate heterogeneities that affect seepage rates can help guide seepage meter placement. Improving our understanding of the causes and types of heterogeneity in lake seepage will provide better data for lake budgets and prediction of mass transfer of solutes or contaminants between lakes and groundwater.  相似文献   

5.
The PULSE analytical model, which calculates daily groundwater discharge on the basis of user‐specified recharge, was originally developed for calibration using streamflow data. This article describes a model application in which groundwater level data constitute the primary control on model input. As a test case, data were analyzed from a small basin in central Pennsylvania in which extensive groundwater level data are available. The timing and intensity of daily water‐level rises are used to ascertain temporal distribution of recharge, and the simulated groundwater discharge hydrograph has shape features that are similar to the streamflow hydrograph. This article does not include details about calibration, but some steps are illustrated and general procedures are described for calibration in specific hydrologic studies. The PULSE model can be used to assess results of fully automated base flow methods and can be used to define groundwater recharge and discharge at a relatively small time scale.  相似文献   

6.
Regional groundwater flow in high mountainous terrain is governed by a multitude of factors such as geology, topography, recharge conditions, structural elements such as fracturation and regional fault zones as well as man‐made underground structures. By means of a numerical groundwater flow model, we consider the impact of deep underground tunnels and of an idealized major fault zone on the groundwater flow systems within the fractured Rotondo granite. The position of the free groundwater table as response to the above subsurface structures and, in particular, with regard to the influence of spatial distributed groundwater recharge rates is addressed. The model results show significant unsaturated zones below the mountain ridges in the study area with a thickness of up to several hundred metres. The subsurface galleries are shown to have a strong effect on the head distribution in the model domain, causing locally a reversal of natural head gradients. With respect to the position of the catchment areas to the tunnel and the corresponding type of recharge source for the tunnel inflows (i.e. glaciers or recent precipitation), as well as water table elevation, the influence of spatial distributed recharge rates is compared to uniform recharge rates. Water table elevations below the well exposed high‐relief mountain ridges are observed to be more sensitive to changes in groundwater recharge rates and permeability than below ridges with less topographic relief. In the conceptual framework of the numerical simulations, the model fault zone has less influence on the groundwater table position, but more importantly acts as fast flow path for recharge from glaciated areas towards the subsurface galleries. This is in agreement with a previous study, where the imprint of glacial recharge was observed in the environmental isotope composition of groundwater sampled in the subsurface galleries. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

7.
A hydrogeophysical survey is performed at small earthen dam that overlies a confined aquifer. The structure of the dam has not shown evidence of anomalous seepage internally or through the foundation prior to the survey. However, the surface topography is mounded in a localized zone 150 m downstream, and groundwater discharges from this zone periodically when the reservoir storage is maximum. We use self‐potential and electrical resistivity tomography surveys with seismic refraction tomography to (1) determine what underlying hydrogeologic factors, if any, have contributed to the successful long‐term operation of the dam without apparent indicators of anomalous seepage through its core and foundation; and (2) investigate the hydraulic connection between the reservoir and the seepage zone to determine whether there exists a potential for this success to be undermined. Geophysical data are informed by hydraulic and geotechnical borehole data. Seismic refraction tomography is performed to determine the geometry of the phreatic surface. The hydro‐stratigraphy is mapped with the resistivity data and groundwater flow patterns are determined with self‐potential data. A self‐potential model is constructed to represent a perpendicular profile extending out from the maximum cross‐section of the dam, and self‐potential data are inverted to recover the groundwater velocity field. The groundwater flow pattern through the aquifer is controlled by the bedrock topography and a preferential flow pathway exists beneath the dam. It corresponds to a sandy‐gravel layer connecting the reservoir to the downstream seepage zone.  相似文献   

8.
Pore water radon (222Rn) distributions from Indian River Lagoon, Florida, are characterized by three zones: a lower zone where pore water 222Rn and sediment-bound radium (226Ra) are in equilibrium and concentration gradients are vertical; a middle zone where 222Rn is in excess of sediment-bound 226Ra and concentration gradients are concave-downward; and an upper zone where 222Rn concentration gradients are nearly vertical. These 222Rn data are simulated in a one-dimensional numerical model including advection, diffusion, and non-local exchange to estimate magnitudes of submarine groundwater discharge components (fresh or marine). The numerical model estimates three parameters, fresh groundwater seepage velocity, irrigation intensity, and irrigation attenuation, using two Monte Carlo (MC) simulations that (1) ensure the minimization algorithm converges on a global minimum of the merit function and the parameter estimates are consistent within this global minimum, and (2) provide 90% confidence intervals on the parameter estimates using the measured 222Rn activity variance. Model estimates of seepage velocities and discharge agree with previous estimates obtained from numerical groundwater flow models and seepage meter measurements and show the fresh water component decreases offshore and varies seasonally by a factor of nine or less. Comparison between the discharge estimates and precipitation patterns suggests a mean residence time in unsaturated and saturated zones on the order of 5 to 7 months. Irrigation rates generally decrease offshore for all sampling periods. The mean irrigation rate is approximately three times greater than the mean seepage velocity although the ranges of irrigation rates and seepage velocities are the same. Possible mechanisms for irrigation include density-driven convection, wave pumping, and bio-irrigation. Simulation of both advection and irrigation allows the separation of submarine groundwater discharge into fresh groundwater and (re)circulated lagoon water.  相似文献   

9.
地下水的空间分布及运移规律是建立地下水流系统和分析渗流场规律的重要因素.具有直接找水特性的地面核磁共振方法(SNMR)可以用于确定地下水的空间分布,而自然电场法(SP)是一种可以探测地下水流动方向的物探方法,这两种方法的联合使用为地下水研究提供了一种较好的工作模式.本文以滑动带水体探测为例,阐述该工作模式的应用效果.首...  相似文献   

10.
Knowledge of groundwater seepage to and from lakes can be an important component of scientific investigations involving water and geochemical budgets. Measurements of groundwater seepage at Lake Anna, a man-made lake in central Virginia, show that inflow to the lake occurs even under dry summer conditions. Seepage rates were found to be correlated with the elevation of the near-shore water table, which responded rapidly to rainfall events in the fractured rock terrane in which Lake Anna is located. Seepage rates did not decline uniformly with distance offshore. This result contrasts with those for lakes underlain by relatively homogeneous porous media where measurements generally confirm the prediction that seepage rates drop off exponentially with distance from shore. The along shore variability of seepage rates in Lake Anna was related to a topographic index that is used to describe drainage from hillslopes. This suggests that seepage in impoundments such as Lake Anna may be strongly controlled by drainage pathways that pre-date the lake.  相似文献   

11.
This paper describes the formulation and application of a coupled unsaturated/saturated model framework developed to investigate the impact of mining on catchment water yield and groundwater dynamics. The model conceptualization was implemented in both a finite‐element (SUSCAT) and finite‐difference (WEC‐C) solution scheme and found to give similar results. The model framework simulates a coupled surface‐water and groundwater system in which a physically based solution scheme was used to simulate one‐dimensional movement through the unsaturated zone, and a distributed model was used to simulate two‐dimensional saturated groundwater flow. Each soil column comprises a series of layers, each layer being connected to adjacent cells. Subsurface lateral flow is considered when any cell within a layer develops a saturated thickness. Simulation results presented are based on a catchment in the Darling Range, Western Australia that was progressively mined and subsequently rehabilitated. The results predicted the groundwater system beneath the mine areas to have a peak rise owing to mining of between 2 and 4 m. Six years after mining, and following vegetation rehabilitation, the groundwater rise had reduced to 1 m above simulated unmined levels. The corresponding streamflow increase as a result of mining was estimated to peak at 21 mm/year and declined to 7·4 mm/year eight years after revegetation of the mined areas. The simulated groundwater response and streamflow results derived from both models were found to be consistent with observed data. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

12.
Anisotropy and heterogeneity of hydraulic conductivity (K) are seldom considered in models of mire hydrology. We investigated the effect of anisotropy and heterogeneity on groundwater flow in bog peat using a steady‐state groundwater model. In five model simulations, four sets of K data were used. The first set comprised measured K values from an anisotropic and heterogeneous bog peat. These data were aggregated to produce the following simplified data sets: an isotropic and heterogeneous distribution of K; an isotropic and homogeneous distribution; and an anisotropic and homogeneous distribution. We demonstrate that, where anisotropy and heterogeneity exist, groundwater flow in bog peat is complex. Fine‐scale variations in K have the potential to influence patterns and rates of groundwater flow. However, for our data at least, it is heterogeneity and not anisotropy that has the greater influence on producing complex patterns of groundwater flow. We also demonstrate that patterns and rates of groundwater flow are simplified and reduced when measured K values are aggregated to create a more uniform distribution of K. For example, when measured K values are aggregated to produce isotropy and homogeneity, the rate of modelled seepage is reduced by 28%. We also show that when measured K values are used, the presence of a drainage ditch can increase seepage through a modelled cross‐section. Our work has implications for the accurate interpretation of hydraulic head data obtained from peat soils, and also the understanding of the effect of drainage ditches on patterns and rates of groundwater flow. Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

13.
地下水动态受水文因素影响较大,对地震和构造活动具有较灵敏的响应。判别并排除各种水文干扰,确认地下水在地震前的异常变化,对提高地震分析预报能力,具有重要作用。马家沟矿井水位动态观测层与地下水开采层为同一含水层,井水位于2010年出现破年变异常,加速持续上升,截至2015年,最大上升幅度约30 m。依据该井水文地质环境特征,根据唐山市区2001-2015年地下水位、降雨量、地下水开采量实测资料,建立合理多元回归模型和三维地下水流动模型,发现地下水开采量减少应为影响马家沟矿井水位动态的可能因素。文中采取的异常识别与分析方法,可为其他类似井孔的地下水动态异常识别及判定提供一定借鉴。  相似文献   

14.
Coastal areas are usually the preferred place of habitation for human beings. Anthropogenic activities such as the construction of high‐rise buildings and underground transport systems usually require extensive deep foundations and ground engineering works, which may unintentionally modify the coastal groundwater system because the construction materials of foundations are usually of low hydraulic conductivity. In this paper, the impact of these building foundations on the groundwater regime is studied using hypothetical flow and transport models. Various possible realizations of foundation distributions are generated using stochastic parameters derived from a topographical map of an actual coastal area in Hong Kong. The effective hydraulic conductivity is first calculated for different realizations and the results show that the effective hydraulic conductivity can be reduced significantly. Then a hypothetical numerical model based on FEFLOW is set up to study the change of hydraulic head, groundwater discharge, and saltwater‐fresh water interface. The groundwater level and flow are modified to various degrees, depending on the foundations percentage and the distribution pattern of the buildings. When the foundations percentage is high and the building foundations are aggregated, the hydraulic head is raised significantly and the originally one‐dimensional groundwater flow field becomes complicated. Seaward groundwater discharge will be reduced and some groundwater may become seepage through the ground surface. The transport model shows that, after foundations are added, overall the seawater and fresh groundwater interface moves landward, so extensive foundations may induce seawater intrusion. It is believed that the modification of the coastal groundwater system by building foundations may have engineering and environmental implications, such as submarine groundwater discharge, foundation corrosion, and slope stability. Copyright © 2007 John Wiley & Sons, Ltd.  相似文献   

15.
Protection of fens–wetlands dependent on groundwater discharge–requires characterization of groundwater sources and stresses. Because instrumentation and numerical modeling of fens is labor intensive, easy-to-apply methods that model fen distribution and their vulnerability to development are desirable. Here we demonstrate that fen areas can be simulated using existing steady-state MODFLOW models when the unsaturated zone flow (UZF) package is included. In cells where the water table is near land surface, the UZF package calculates a head difference and scaled conductance at these “seepage drain” cells to generate average rates of vertical seepage to the land. This formulation, which represents an alternative to blanketing the MODFLOW domain with drains, requires very little input from the user because unsaturated flow-routing is inactive and results are primarily driven by easily obtained topographic information. Like the drain approach, it has the advantage that the distribution of seepage areas is not predetermined by the modeler, but rather emerges from simulated heads. Beyond the drain approach, it takes account of intracell land surface variation to explicitly quantify multiple surficial flows corresponding to infiltration, rejected recharge, recharge and land-surface seepage. Application of the method to a basin in southeastern Wisconsin demonstrates how it can be used as a decision-support tool to first, reproduce fen distribution and, second, forecast drawdown and reduced seepage at fens in response to shallow pumping.  相似文献   

16.
Groundwater-flow models depend on hydraulic head and flux observations for evaluation and calibration. A different type of observation—change in storage measured using repeat microgravity—can also be used for parameter estimation by simulating the expected change in gravity from a groundwater model and including the observation misfit in the objective function. The method is demonstrated using new software linked to MODFLOW input and output files and field data from the vicinity of the All American Canal in southeast California, USA. Over a 10-year period following lining of the previously highly permeable canal with concrete, gravity decreased by over 100 μGal (equivalent to about 2.5 m of free-standing water) at some locations as seepage decreased and the remnant groundwater mound dissipated into the aquifer or was removed by groundwater pumping. Simulated gravity from a MODFLOW model closely matched observations, and repeat microgravity data proved useful for constraining both hydraulic conductivity and specific yield estimates. Specific yield estimated using the infinite-horizontal slab approximation agreed well with model-derived values, and the departure from the linear, flat-water-table approximation was small, less than 2%, despite relatively large and dynamic water-table slope. First-order second-moment parameter uncertainty analysis shows reduction in uncertainty for all hydraulic conductivity and specific yield parameter estimates with the addition of repeat microgravity data, as compared to drawdown data alone.  相似文献   

17.
Self-potential (SP) and electrical resistivity measurements are used to investigate seepage at a remote moraine dam in the Sierra Nevada of California. The site is a small terminal moraine impounding roughly 300,000 m3 of water at ~ 3400 m a.s.l. Suspicious fine sediment in a small lake at the dam's downstream toe prompted initial concerns that anomalous seepage may be eroding matrix material from the moraine. 235 individual SP measurements covering the surface of the dam were collected in order to investigate electrokinetic current sources resulting from seepage, while resistivity soundings probed moraine stratigraphy and suggest that the till contains interstitial ice. Contoured SP data reveal a non-uniform voltage distribution over the moraine dam and two distinct negative SP anomalies. The first, located in the central area of the moraine, shows a broad negative SP zone around the crest and increasingly positive SP moving downhill towards both the upstream and downstream toes. This anomaly can be explained by shallow gravitational groundwater flow in the near subsurface combined with upward groundwater flux through evapotranspiration; numerical simulation of the combined effect matches field data well. The second SP anomaly has a tightly localized distribution and can be explained by vertically descending flow into a bedrock fault conduit. Our conceptual seepage model suggests that flow travels from Dana Lake first at the boundary of ice-filled moraine and bedrock before converging on a concentrated channel in the subvertical fault zone. Positive SP near the dam abutments results from groundwater inflow from adjacent hillslopes. Combined analyses suggest that seepage erosion is not currently affecting the moraine dam, and that the sediment observed on the bed of the downstream toe lake is likely a remnant of past outflow events.  相似文献   

18.
The link between groundwater and surface hydrology in a small headwater drainage basin in the zone of glacial deposition of southern Ontario south of the Precambrian Shield was examined for two years. The basin is situated in a discharge zone of a regional aquifer and contains a small treed spring-fed swamp. The swamp exists because of the groundwater and has little effect on the maintenance of streamflow. Groundwater input to the swamp is an order of magnitude larger than precipitation. Groundwater of local and regional origin passes through the swamp by two routes: surface streamlets, where groundwater that emerges at specific seepage points in the swamp is conveyed over the ground surface with little interaction with the swamp itself, and by diffuse seepage in the swamp and through the bed of the stream. While the diffuse seepage input is the smaller component of groundwater it maintains the swamp's saturation. Groundwater input to the swamp from the specific seepage points and diffuse flow varies little over a year; therefore the saturation of the swamp and baseflow from the basin display little seasonal variation compared to other wetland types. The existence of the valley bottom in the headwater basin alters the seasonal and storm hydrology and is important to biogeochemical transformation of emerging groundwater.  相似文献   

19.
This paper, based on a real world case study (Limmat aquifer, Switzerland), compares inverse groundwater flow models calibrated with specified numbers of monitoring head locations. These models are updated in real time with the ensemble Kalman filter (EnKF) and the prediction improvement is assessed in relation to the amount of monitoring locations used for calibration and updating. The prediction errors of the models calibrated in transient state are smaller if the amount of monitoring locations used for the calibration is larger. For highly dynamic groundwater flow systems a transient calibration is recommended as a model calibrated in steady state can lead to worse results than a noncalibrated model with a well-chosen uniform conductivity. The model predictions can be improved further with the assimilation of new measurement data from on-line sensors with the EnKF. Within all the studied models the reduction of 1-day hydraulic head prediction error (in terms of mean absolute error [MAE]) with EnKF lies between 31% (assimilation of head data from 5 locations) and 72% (assimilation of head data from 85 locations). The largest prediction improvements are expected for models that were calibrated with only a limited amount of historical information. It is worthwhile to update the model even with few monitoring locations as it seems that the error reduction with EnKF decreases exponentially with the amount of monitoring locations used. These results prove the feasibility of data assimilation with EnKF also for a real world case and show that improved predictions of groundwater levels can be obtained.  相似文献   

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

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

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