共查询到20条相似文献，搜索用时 31 毫秒
1.
The value of subsidence data in ground water model calibration 总被引：2，自引：0，他引：2
The accurate estimation of aquifer parameters such as transmissivity and specific storage is often an important objective during a ground water modeling investigation or aquifer resource evaluation. Parameter estimation is often accomplished with changes in hydraulic head data as the key and most abundant type of observation. The availability and accessibility of global positioning system and interferometric synthetic aperture radar data in heavily pumped alluvial basins can provide important subsidence observations that can greatly aid parameter estimation. The aim of this investigation is to evaluate the value of spatial and temporal subsidence data for automatically estimating parameters with and without observation error using UCODE-2005 and MODFLOW-2000. A synthetic conceptual model (24 separate cases) containing seven transmissivity zones and three zones each for elastic and inelastic skeletal specific storage was used to simulate subsidence and drawdown in an aquifer with variably thick interbeds with delayed drainage. Five pumping wells of variable rates were used to stress the system for up to 15 years. Calibration results indicate that (1) the inverse of the square of the observation values is a reasonable way to weight the observations, (2) spatially abundant subsidence data typically produce superior parameter estimates under constant pumping even with observation error, (3) only a small number of subsidence observations are required to achieve accurate parameter estimates, and (4) for seasonal pumping, accurate parameter estimates for elastic skeletal specific storage values are largely dependent on the quantity of temporal observational data and less on the quantity of available spatial data. 相似文献
2.
气压作用下, 井-含水层系统中地下水流是一类流体力学问题。 本文应用井壁水流通量边界条件和气压作用下井壁内外水(孔)压平衡条件, 提出了一个井水位随气压变化的解析公式。 解析式表明, 气压系数随时段长度增加而增大, 并趋于气压常数; 气压系数随时段长度的变化只依赖于导水系数与井半径平方的比值(T/r2w), 而与气压变化过程无关; 气压常数只与含水层的一维荷载效率(B)有关, 而与导水系数和井半径无关。 解析解所反映的气压系数与时段长度的关系, 与南溪井实测序列数据分析结果具有很好的一致性。 根据气压系数随时段长度变化过程, 提出了一个参数估计方法, 应用于估计南溪井含水层气压常数和导水系数, 并对本文提出的参数估计方法进行了讨论。 相似文献
3.
Measurements of the tide and groundwater levels in coastal zones are of importance in determining the properties of coastal aquifers. The solution to a one‐dimensional unsteady groundwater flow model in a coastal confined aquifer with sinusoidal fluctuation of the tide shows that the tidal efficiency decreases exponentially with distance and the time lag increases linearly with distance from the coast. The aquifer property described by the ratio of storage coefficient to transmissivity is determined if the damping constant of the tidal efficiency or the slope of the time lag with distance are obtained on the basis of tidal measurements. Hourly observations of the tide and groundwater levels at 10 wells on the northern coast near Beihai, China show that with distance from the coast, tidal efficiency decreases roughly exponentially and the time lag increases roughly linearly. The estimated ratio of storage coefficient to transmissivity of the confined aquifer ranges from 1·169 × 10?6 d m?2 to 1·83 × 10?7 d m?2. For a given transmissivity of 750 m2 d?1, the storage coefficient of the aquifer is 8·7675 × 10?4 with the tidal efficiency method and 1·3725 × 10?4 with the time lag method. The damping constant of the tidal efficiency with distance can be defined as the tidal propagation coefficient. The value of the confined aquifer is determined as 0·0018892 m?1. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
4.
Xiuyu Liang You-Kuan Zhang Keith Schilling 《Stochastic Environmental Research and Risk Assessment (SERRA)》2016,30(1):1-8
Spatiotemporal variations of groundwater level due to a white noise recharge time series and a random transmissivity field in a bounded unconfined aquifer was studied. The analytical solutions for the variance and covariance of groundwater level were derived with non-stationary spectral analyses and superposition principle. It was found that the fluctuations of groundwater level are spatially non-stationary due to a fixed head boundary condition and temporal non-stationary at early time but gradually became stationary as time progresses due to effect of the initial condition. The variation in groundwater level is mainly caused by the random source/sink in the case of temporally random recharge and spatially random transmissivity. The effect of heterogeneity is to increase the variation of groundwater level and the maximum effect occurs close to the constant head boundary because of the linear mean hydraulic gradient. The heterogeneity also enhances the correlation of groundwater level, especially at large time intervals and small spatial distances. 相似文献
5.
Groundwater management decisions are often founded upon estimates of aquifer hydraulic properties, recharge and the rate of groundwater usage. Too often hydraulic properties are unavailable, recharge estimates are very uncertain, and usage is unmetered or infrequently metered over only recent years or estimated using numerical groundwater models decoupled from the drivers of drawdown. This paper extends the HydroSight groundwater time-series package ( http://peterson-tim-j.github.io/HydroSight/ ) to allow the joint estimation of gross recharge, transmissivity, storativity, and daily usage at multiple production bores. A genetic evolutionary scheme was extended from estimating time-series model parameters to also estimating time series of usage that honor metered volumes at each production bore and produces (1) the best fit with the observed hydrograph and (2) plausible estimates of actual evapotranspiration and hence recharge. The reliability of the approach was rigorously tested. Repeated calibration of models for four bores produced estimates of transmissivity, storativity, and mean recharge that varied by a factor of 0.22-0.32, 0.13-0.2, and 0.03-0.48, respectively, when recharge boundary effects were low and the error in monthly, quarterly, and biannual metered usage was generally <10%. Application to the 30 observation bores within the Warrion groundwater management area (Australia), produced a coefficient of efficiency of ≥0.80 at 22 bores and ≥0.90 at 12 bores. The aquifer transmissivity and storativity were reasonably estimated, and were consistent with independent estimates, while mean gross recharge may be slightly overestimated. Overall, the approach allows greater insights from the available data and provides opportunity for the exploration of usage and climatic scenarios. 相似文献
6.
Abstract. A method to calculate aquifer transmissivity, storage coefficient, and the leakage coefficient from pumping test data for a leaky aquifer is presented. The method is carried out by a computer program and is based on a minimization of the sum of squares of differences between drawdown in the observation well and the theoretical values from the Hantush and Jacob formula. No user defined starting points are necessary. Random error estimates for the parameters are given. Applications of the method are illustrated using data from pumping tests performed in leaky aquifers at the Cauca River Valley, Colombia. 相似文献
7.
The hydraulic diffusivity gives a measure of diffusion speed of pressure disturbances in groundwater system; large values of hydraulic diffusivity lead to fast propagation of signals in aquifer. This research provides a novel design and derives spectral representation to determine hydraulic diffusivity using spectral analysis of groundwater levels coupled with time-dependent boundary adjacent to marine system and no flow boundary in aquifer system. To validate the proposed method, water levels of fluctuated boundary and groundwater well in a sandy confined aquifer were collected. The hydraulic diffusivity is then obtained by an inverse process in the non-linear complex form of spectral relationship. The method essentially is constructed on the conceptual design of natural forcing transmitted in large aquifer. It is unlike the conventional field pumping test which is only used to determine hydraulic properties of groundwater in small range around the well. Hydraulic diffusivity of the confined aquifer is determined using real observation and then checked by comparing to the published range. It suggests that without local aquifer test to estimate hydraulic diffusivity in a coastal aquifer using spectral representation with its relevant flow system and boundary has become feasible. 相似文献
8.
Experimental hydraulic heads and electrical (self-potential) signals associated with a pumping test were used in an inverse model to estimate the transmissivity distribution of a real aquifer. Several works reported in the literature show that there is a relatively good linear relationship between the hydraulic heads in the aquifer and electrical signals measured at the ground surface. In this experimental test field, first, the current coupling coefficient was determined by the best fit between experimental and modeled self-potential signals at the end of the pumping phase. Soon afterward, with the hydraulic heads obtained from the self-potential signals, the transmissivity distribution of the aquifer was conditioned by means an inverse model based on the successive linear estimator (SLE). To further substantiate the estimated T field from the SLE analysis, we analyzed the drawdown rate, the derivative of the drawdown with respect to the ln(t), because the drawdown rate is highly sensitive to the variability in the transmissivity field. In our opinion, these results show that self-potential signals allow the monitoring of subsurface flow in the course of pumping experiments, and that electrical potentials serve as a good complement to piezometric observations to condition and characterize the transmissivity distribution of an aquifer. 相似文献
9.
Priyanka Bangalore Nagaraj Mohan Kumar Mandalagiri Subbarayappa Vouillamoz Jean-Michel Johan Hoareau 《水文研究》2021,35(10):e14395
Estimation of hydraulic parameters is essential to understand the interaction between groundwater flow and seawater intrusion. Though several studies have addressed hydraulic parameter estimation, based on pumping tests as well as geophysical methods, not many studies have addressed the problem with clayey formations being present. In this study, a methodology is proposed to estimate anisotropic hydraulic conductivity and porosity values for the coastal aquifer with unconsolidated formations. For this purpose, the one-dimensional resistivity of the aquifer and the groundwater conductivity data are used to estimate porosity at discrete points. The hydraulic conductivity values are estimated by its mutual dependence with porosity and petrophysical parameters. From these estimated values, the bilinear relationship between hydraulic conductivity and aquifer resistivity is established based on the clay content of the sampled formation. The methodology is applied on a coastal aquifer along with the coastal Karnataka, India, which has significant clayey formations embedded in unconsolidated rock. The estimation of hydraulic conductivity values from the established correlations has a correlation coefficient of 0.83 with pumping test data, indicating good reliability of the methodology. The established correlations also enable the estimation of horizontal hydraulic conductivity on two-dimensional resistivity sections, which was not addressed by earlier studies. The inventive approach of using the established bilinear correlations at one-dimensional to two-dimensional resistivity sections is verified by the comparison method. The horizontal hydraulic conductivity agrees with previous findings from inverse modelling. Additionally, this study provides critical insights into the estimation of vertical hydraulic conductivity and an equation is formulated which relates vertical hydraulic conductivity with horizontal. Based on the approach presented, the anisotropic hydraulic conductivity of any type aquifer with embedded clayey formations can be estimated. The anisotropic hydraulic conductivity has the potential to be used as an important input to the groundwater models. 相似文献
10.
Marathwada Agricultural University, Pharbani, has developed about 560 hectares of Wagarwadi watershed in Pharbani district since 1987. Groundwater monitoring on 16 observations wells at weekly intervals commenced in January 1992, and rainfall and pan evaporation has been measured daily at a hydrometeorological station situated in the nearby university campus. Aquifer parameters, namely, transmissivity and specific yield, have been estimated by carrying out a pumping test on a large diameter well. Groundwater recharge resulting from rainfall has been estimated using a water balance model of the soil moisture zone considering soil zone thickness and crops grown. The SCS (Soil Conservation Service) curve number method was used for surface runoff estimation. The groundwater flow model has been constructed using the nested squares, finite difference method. Nested square meshes of sizes 160 m×160 m and 80 m×80 m have been used and the steady-state condition of aquifer system was simulated in the model assuming the June 1992 water level configuration under equilibrium conditions. The model has been calibrated for transient conditions incorporating additional seepage from the water harvesting structures and their contribution to the groundwater regime has been assessed. © 1998 John Wiley & Sons, Ltd. 相似文献
11.
A single recovery type curve from Theis'' exact solution 总被引：2，自引：0，他引：2
The Theis type curve matching method and the Cooper-Jacob semilog method are commonly used for estimation of transmissivity and storage coefficient of infinite, homogeneous, isotropic, confined aquifers from drawdown data of a constant rate pumping test. Although these methods are based on drawdown data, they are often applied indiscriminately to analyze both drawdown and recovery data. Moreover, the limitations of drawdown type curve to analyze recovery data collected after short pumping times are not well understood by the practicing engineers. This often may result in an erroneous interpretation of such recovery data. In this paper, a novel but simple method is proposed to determine the storage coefficient as well as transmissivity from recovery data measured after the pumping period of an aquifer test. The method eliminates the dependence on pumping time effects and has the advantage of employing only one single recovery type curve. The method based on the conversion of residual drawdown to recovered drawdown (buildup) data plotted versus a new equivalent time (delta(t) x t(p)/t(p) + delta(t)). The method uses the recovery data in one observation point only, and does not need the initial water level h0, which may be unknown. The accuracy of the method is checked with three sets of field data. This method appears to be complementary to the Cooper-Jacob and Theis methods, as it provides values of both storage coefficient and transmissivity from recovery data, regardless of pumping duration. 相似文献
12.
Abstract In order to calculate the transmissivity from the inverse problem corresponding to the groundwater flow in an isotropic horizontal aquifer, a numerical conservative approach is tested. The method deals with triangulation of the domain and applies the conservation of mass to elements of the mesh using the harmonic mean for internodal transmissivities. An optimal sweeping algorithm is used to evaluate nodal transmissivities from one element to another with a minimal relative error accumulation. The practical importance of the method is demonstrated through two synthetic examples representing those experienced in the field, then through application to a Moroccan aquifer. The computed hydraulic head is well fitted to the reference one, which confirms the validity of the identified transmissivity model. 相似文献
13.
Tensor Hydraulic Conductivity Estimation for Heterogeneous Aquifers under Unknown Boundary Conditions 
下载免费PDF全文
下载免费PDF全文 A physically based inverse method is developed using hybrid formulation and coordinate transform to simultaneously estimate hydraulic conductivity tensors, steady‐state flow field, and boundary conditions for a confined aquifer under ambient flow or pumping condition. Unlike existing indirect inversion techniques, the physically based method does not require forward simulations to assess model‐data misfits. It imposes continuity of hydraulic head and Darcy fluxes in the model domain while incorporating observations (hydraulic heads, Darcy fluxes, or well rates) at measurement locations. Given sufficient measurements, it yields a well‐posed inverse system of equations that can be solved efficiently with coarse grids and nonlinear optimization. When pumping and injection are active, well rates are used as measurements and flux sampling is not needed. The method is successfully tested on synthetic aquifer problems with regular and irregular geometries, different hydrofacies and flow patterns, and increasing conductivity anisotropy ratios. All problems yield stable inverse solutions under increasing head measurement errors. For a given set of observations, inversion accuracy is strongly affected by the conductivity anisotropy ratio. Conductivity estimation is also affected by flow pattern: within a hydrofacies, when Darcy flux component is very small, the corresponding directional conductivity perpendicular to streamlines becomes less identifiable. Finally, inversion is successful even if the location of aquifer boundaries is unknown. In this case, the inversion domain is defined by the location of the measurements. 相似文献
14.
While tomographic inversion has been successfully applied to laboratory- and field-scale tests, here we address the new issue of scale that arises when extending the method to a basin. Specifically, we apply the hydraulic tomography (HT) concept to jointly interpret four multiwell aquifer tests in a synthetic basin to illustrate the superiority of this approach to a more traditional Theis analysis of the same tests. Transmissivity and storativity are estimated for each element of a regional numerical model using the geostatistically based sequential successive linear estimator (SSLE) inverse solution method. We find that HT inversion is an effective strategy for incorporating data from potentially disparate aquifer tests into a basin-wide aquifer property estimate. The robustness of the SSLE algorithm is investigated by considering the effects of noisy observations, changing the variance of the true aquifer parameters, and supplying incorrect initial and boundary conditions to the inverse model. Ground water flow velocities and total confined storage are used as metrics to compare true and estimated parameter fields; they quantify the effectiveness of HT and SSLE compared to a Theis solution methodology. We discuss alternative software that can be used for implementing tomography inversion. 相似文献
15.
The buried-valley aquifers that are common in the glacial deposits of the northern hemisphere are a typical case of the strip aquifers that occur in many parts of the world. Pumping from a narrow strip aquifer leads to much greater drawdown and much more distant drawdown effects then would occur in a sheet aquifer with a similar transmissivity and storage coefficient. Widely used theories for radial flow to wells, such as the Theis equation, are not appropriate for narrow strip aquifers. Previously published theory for flow to wells in semiconfined strip aquifers is reviewed and a practical format of the type curves for pumping-test analysis is described. The drawdown response of strip aquifers to pumping tests is distinctive, especially for observation wells near the pumped well. A case study is presented, based on extensive pumping test experience for the Estevan Valley Aquifer in southern Saskatchewan, Canada. Evaluation of groundwater resources in such buried-valley aquifers needs to take into account the unusually large drawdowns in response to pumping. 相似文献
16.
Accurate estimation of aquifer parameters, especially from crystalline hard rock area, assumes a special significance for management of groundwater resources. The aquifer parameters are usually estimated through pumping tests carried out on water wells. While it may be costly and time consuming for carrying out pumping tests at a number of sites, the application of geophysical methods in combination with hydro-geochemical information proves to be potential and cost effective to estimate aquifer parameters. Here a method to estimate aquifer parameters such as hydraulic conductivity, formation factor, porosity and transmissivity is presented by utilizing electrical conductivity values analysed via hydro-geochemical analysis of existing wells and the respective vertical electrical sounding (VES) points of Sindhudurg district, western Maharashtra, India. Further, prior to interpolating the distribution of aquifer parameters of the study area, variogram modelling was carried out using data driven techniques of kriging, automatic relevance determination based Bayesian neural networks (ARD-BNN) and adaptive neuro-fuzzy neural networks (ANFIS). In total, four variogram model fitting techniques such as spherical, exponential, ARD-BNN and ANFIS were compared. According to the obtained results, the spherical variogram model in interpolating transmissivity, ARD-BNN variogram model in interpolating porosity, exponential variogram model in interpolating aquifer thickness and ANFIS variogram model in interpolating hydraulic conductivity outperformed rest of the variogram models. Accordingly, the accurate aquifer parameters maps of the study area were produced by using the best variogram model. The present results suggest that there are relatively high value of hydraulic conductivity, porosity and transmissivity at Parule, Mogarne, Kudal, and Zarap, which would be useful to characterize the aquifer system over western Maharashtra. 相似文献
17.
Geophysical Assessment of Groundwater Potential: A Case Study from Mian Channu Area,Pakistan 下载免费PDF全文
下载免费PDF全文 An integrated study using geophysical method in combination with pumping tests and geochemical method was carried out to delineate groundwater potential zones in Mian Channu area of Pakistan. Vertical electrical soundings (VES) using Schlumberger configuration with maximum current electrode spacing (AB/2 = 200 m) were conducted at 50 stations and 10 pumping tests at borehole sites were performed in close proximity to 10 of the VES stations. The aim of this study is to establish a correlation between the hydraulic parameters obtained from geophysical method and pumping tests so that the aquifer potential can be estimated from the geoelectrical surface measurements where no pumping tests exist. The aquifer parameters, namely, transmissivity and hydraulic conductivity were estimated from Dar Zarrouyk parameters by interpreting the layer parameters such as true resistivities and thicknesses. Geoelectrical succession of five‐layer strata (i.e., topsoil, clay, clay sand, sand, and sand gravel) with sand as a dominant lithology was found in the study area. Physicochemical parameters interpreted by World Health Organization and Food and Agriculture Organization were well correlated with the aquifer parameters obtained by geoelectrical method and pumping tests. The aquifer potential zones identified by modeled resistivity, Dar Zarrouk parameters, pumped aquifer parameters, and physicochemical parameters reveal that sand and gravel sand with high values of transmissivity and hydraulic conductivity are highly promising water bearing layers in northwest of the study area. Strong correlation between estimated and pumped aquifer parameters suggest that, in case of sparse well data, geophysical technique is useful to estimate the hydraulic potential of the aquifer with varying lithology. 相似文献
18.
R. T. Miller 《Ground water》1984,22(5):532-537
The U.S. Geological Survey is studying the potential for storage of heated water in a sandstone aquifer in St. Paul, Minnesota. The efficiency of the aquifer to store thermal energy is related, in part, to the hydrogeologic characteristics of the aquifer. The movement of heat away from the injection well is directly related to the anisotropy. Aquifer tests were conducted to determine the anisotropy near the heated-water injection well. The maximum and minimum values of transmissivity along the principal directions of the hydraulic conductivity tensors of the Ironton and Galesville Sandstones in St. Paul, Minnesota are approximately 1,090 and 480 feet squared per day. The storage coefficient is 4.5 × 10−5. These values represent the average of four determinations of nonsteady flow to a well in an idealized infinite anisotropic aquifer. Analysis of the values of transmissivity and storage coefficient for hypothetical changes in location of two of the monitoring wells where depth-deviation surveys were not available indicates that computed transmissivities vary less than 5 percent and storage coefficients vary less than ±6 percent. 相似文献
19.
Estimating River Conductance from Prior Information to Improve Surface‐Subsurface Model Calibration 下载免费PDF全文
下载免费PDF全文 Yohann Cousquer Alexandre Pryet Nicolas Flipo Célestine Delbart Alain Dupuy 《Ground water》2017,55(3):408-418
Most groundwater models simulate stream‐aquifer interactions with a head‐dependent flux boundary condition based on a river conductance (CRIV). CRIV is usually calibrated with other parameters by history matching. However, the inverse problem of groundwater models is often ill‐posed and individual model parameters are likely to be poorly constrained. Ill‐posedness can be addressed by Tikhonov regularization with prior knowledge on parameter values. The difficulty with a lumped parameter like CRIV, which cannot be measured in the field, is to find suitable initial and regularization values. Several formulations have been proposed for the estimation of CRIV from physical parameters. However, these methods are either too simple to provide a reliable estimate of CRIV, or too complex to be easily implemented by groundwater modelers. This paper addresses the issue with a flexible and operational tool based on a 2D numerical model in a local vertical cross section, where the river conductance is computed from selected geometric and hydrodynamic parameters. Contrary to other approaches, the grid size of the regional model and the anisotropy of the aquifer hydraulic conductivity are also taken into account. A global sensitivity analysis indicates the strong sensitivity of CRIV to these parameters. This enhancement for the prior estimation of CRIV is a step forward for the calibration and uncertainty analysis of surface‐subsurface models. It is especially useful for modeling objectives that require CRIV to be well known such as conjunctive surface water‐groundwater use. 相似文献
20.
Effectively managing groundwater relies heavily on estimating the amount of precipitation that may infiltrate the subsurface and supply groundwater. In this study, we present a novel estimation method based on a stochastic approach to evaluate the quantity of precipitation that may recharge groundwater. The precipitation recharge coefficient is also investigated based on an unconfined aquifer with an unbound, infinitely extended boundary condition. Moreover, a spectrum's relationship to the precipitation and groundwater level variation is also derived. The precipitation recharge coefficient can be obtained from the solution of the spectrum equation. Furthermore, sensitivity analysis is performed in order to determine the key variable on the precipitation recharge coefficient. Analysis results indicate that the location of an observation well affects the estimated precipitation recharge coefficient. If the precipitation recharge area is large enough, the precipitation recharge coefficient becomes insensitive to the location of the observation well. The spectrum's relationship between the precipitation recharge and groundwater level variation is also applied when estimating the precipitation recharge coefficient upstream of the Cho‐Shui River alluvial fan. According to those results, the precipitation recharge coefficient is 0·03 and the amount of groundwater recharge from precipitation is 35 million tons of water annually upstream of the Cho‐Shui River alluvial fan. Copyright © 2000 John Wiley & Sons, Ltd. 相似文献