共查询到20条相似文献，搜索用时 796 毫秒
1.
The vertical hydraulic conductivity (Kv), elastic (Sske), and inelastic (Sskv) skeletal specific storage of aquitards are three of the most critical parameters in land subsidence investigations. Two new analytic methods are proposed to estimate the three parameters. The first analytic method is based on a new concept of delay time ratio for estimating Kv and Sske of an aquitard subject to long‐term stable, cyclic hydraulic head changes at boundaries. The second analytic method estimates the Sskv of the aquitard subject to linearly declining hydraulic heads at boundaries. Both methods are based on analytical solutions for flow within the aquitard, and they are jointly employed to obtain the three parameter estimates. This joint analytic method is applied to estimate the Kv, Sske, and Sskv of a 34.54‐m thick aquitard for which the deformation progress has been recorded by an extensometer located in Shanghai, China. The estimated results are then calibrated by PEST (Doherty 2005), a parameter estimation code coupled with a one‐dimensional aquitard‐drainage model. The Kv and Sske estimated by the joint analytic method are quite close to those estimated via inverse modeling and performed much better in simulating elastic deformation than the estimates obtained from the stress‐strain diagram method of Ye and Xue (2005). The newly proposed joint analytic method is an effective tool that provides reasonable initial values for calibrating land subsidence models. 相似文献
2.
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. 相似文献
3.
Inverse modelling using PS‐InSAR data for improved land subsidence simulation in Las Vegas Valley,Nevada 
下载免费PDF全文
下载免费PDF全文 Las Vegas Valley has had a long history of groundwater development and subsequent surface deformation. InSAR interferograms have revealed detailed and complex spatial patterns of subsidence in the Las Vegas Valley area that do not coincide with major pumping regions. This research represents the first effort to use high spatial and temporal resolution subsidence observations from InSAR and hydraulic head data to inversely calibrate transmissivities (T), elastic and inelastic skeletal storage coefficients (Ske and Skv) of the developed‐zone aquifer and conductance (CR) of the basin‐fill faults for the entire Las Vegas basin. The results indicate that the subsidence observations from InSAR are extremely beneficial for accurately quantifying hydraulic parameters, and the model calibration results are far more accurate than when using only groundwater levels as observations, and just a limited number of subsidence observations. The discrepancy between distributions of pumping and greatest levels of subsidence is found to be attributed to spatial variations in clay thickness. The Eglington fault separates thicker interbeds to the northwest from thinner interbeds to the southeast and the fault may act as a groundwater‐flow barrier and/or subsidence boundary, although the influence of the groundwater barrier to this area is found to be insignificant. Copyright © 2016 John Wiley & Sons, Ltd. 相似文献
4.
Assessing radial transmissivity variation in heterogeneous aquifers using analytical techniques 下载免费PDF全文
下载免费PDF全文 Pumping tests are one of the most commonly used in situ testing techniques for assessing aquifer hydraulic properties. Numerous researches have been conducted to predict the effects of aquifer heterogeneity on the groundwater levels during pumping tests. The objectives of the present work were as follows: (1) to predict drawdown conditions and to estimate aquifer properties during pumping tests undertaken in radially symmetric heterogeneous aquifers, and (2) to identify a method for assessing the transmissivity field along the radial coordinate in radially symmetric and fully heterogeneous transmissivity fields. The first objective was achieved by expanding an existing analytical drawdown formulation that was valid for a radially symmetric confined aquifer with two concentric zones around the pumping well to an N concentric zone confined aquifer having a constant transmissivity value within each zone. The formulation was evaluated for aquifers with three and four concentric zones to assess the effects of the transmissivity field on the drawdown conditions. The specific conditions under which aquifer properties could be identified using traditional methods of analysis were also evaluated. The second objective was achieved by implementing the inverse solution algorithm (ISA), which was developed for petroleum reservoirs to groundwater aquifer settings. The results showed that the drawdown values are influenced by a volumetric integral of a weighting function and the transmissivity field within the cone of depression. The weighting function migrates in tandem with the expanding cone of depression. The ability of the ISA to predict radially symmetric and log‐normally distributed transmissivity fields was assessed against analytical and numerical benchmarks. The results of this investigation indicated that the ISA method is a viable technique for evaluating the radial transmissivity variations of heterogeneous aquifer settings. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
5.
Ground subsidence is a critical problem in populated areas, which requires attention and routine monitoring. Specifically, subsidence related to the exploitation of aquifer systems is of particular interest in terms of risk for the possible damage to buildings and availability of water resources. We propose to realise a survey of the aquifer systems near the city of Lorca, southern Spain, using Synthetic Aperture Radar Interferometry (InSAR) with a light method necessitating a small number of interferograms rather than analysing hundreds of radar images. We computed nine interferograms with good coherence over the years 2004 and 2005. We analysed them by a two-step filtering process in order to characterise the atmospheric contribution and the ground deformation. We are then able to identify two deforming areas and quantify subsidence rates at 0.20 ± 0.02 mm/day south of Lorca and southeast to the city of Alhama de Murcia with the same efficiency as the studies using sophisticated InSAR methods. Using continuous records of water levels in boreholes, we deduced a mean storage coefficient (Sk) at 4.1 × 10− 3 which is very close to the previous published values, and coincides in this case to an inelastic deformation mode of all the system. With this value, we determine the groundwater level variations from the measurements of the ground subsidence. 相似文献
6.
Javad Ashjari 《Ground water》2013,51(1):122-127
An aquifer test is used mostly to determine the storage coefficient and transmissivity. Although residual drawdown data are widely used in estimating the transmissivity of aquifers, the estimation of storage coefficients with recovery data is controversial. Some researchers have proposed methods to estimate storage coefficients with recovery data by assuming equality of storage coefficients for the recovery and pumping periods (S = S′). The aim of this study is to determine storage coefficients without such an assumption, that is, S≠S′. The method is a modified version of Banton‐Bangoy's method without considering drawdown data due to pumping. Drawdown is plotted vs. the logarithmic ratio (t′/t) or time since pumping stopped to the duration of pumping and the ratio of storage coefficient during recovery to the storage coefficient from the pumping period (S′/S). The method is verified with one case study and two synthetic examples. Thus, it is possible to determine storage coefficient of pumping period accurately without any data from pumping period by recovery data. 相似文献
7.
Carl Keller 《Ground water》2017,55(2):244-254
This study describes a new technique for measuring the head profile in a geologic formation. The technique provides rapid, low cost information on the depth of water‐producing zones and aquitards in heterogeneous aquifers, yielding estimates of hydraulic heads in each zone while identifying any potential for cross contamination between zones. The measurements can be performed in a typical borehole in just a few hours. The procedure uses both the continuous transmissivity profile obtained by the installation (eversion) of a flexible borehole liner into an open borehole and the subsequent removal (inversion) of the same liner from the borehole. The method is possible because of the continuous transmissivity profile (T profile described by Keller et al. 2014) obtained by measuring the rate of liner eversion under a constant driving head. The hydraulic heads of producing zones are measured using the reverse head profile (RHP) method (patent no. 9,008,971) based on a stepwise inversion of the borehole liner. As each interval of the borehole is uncovered by inversion of the liner, the head beneath the liner is allowed to equilibrate to a steady‐state value. The individual hydraulic heads contributing to each measurement are calculated using the measured transmissivity for each zone. Application of the RHP method to a sedimentary bedrock borehole in New Jersey verified that it reproduced the head distribution obtained the same day in the same borehole instrumented with a multilevel sampling system. 相似文献
8.
In this paper, we perform an inverse method to simultaneously estimate aquifer parameters, initial condition, and boundary conditions in groundwater modelling. The parameter estimation is extended to a complete inverse problem that makes the calibrated groundwater flow model more realistic. The adjoint state method, the gradient search method, and the least square error algorithm are combined to build the optimization procedure. Horizontal two‐dimensional groundwater flow in a confined aquifer is exemplified to demonstrate the correlation between unknowns, the contribution of observation, as well as the suitability of applying the inverse method. The correlation analysis shows the connection between storage coefficient and initial condition. Besides, transmissivity and boundary conditions are also highly correlated. More observations at different location and time are necessary to provide sufficient information. A time series of unsteady head is requested for estimation of storage coefficient and initial condition. Observation near boundary is very effective for boundary condition estimation. The observation at pumping well mostly contributes to the estimation of transmissivity. According to all observations, it is possible to identify parameters, initial condition, and boundary condition simultaneously. Furthermore, the results not only illustrate the traditional assumption of known boundary condition but also initial condition, which may cause an incorrect estimation. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
9.
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. 相似文献
10.
A conceptual model of anisotropic and dynamic permeability is developed from hydrogeologic and hydromechanical characterization of a foliated, complexly fractured, crystalline rock aquifer at Gates Pond, Berlin, Massachusetts. Methods of investigation include aquifer‐pumping tests, long‐term hydrologic monitoring, fracture characterization, downhole heat‐pulse flow meter measurements, in situ extensometer testing, and earth tide analysis. A static conceptual model is developed from observations of depth‐dependent and anisotropic permeability that effectively compartmentalizes the aquifer as a function of foliation intensity. Superimposed on the static model is dynamic permeability as a function of hydraulic head in which transient bulk aquifer transmissivity is proportional to changes in hydraulic head due to hydromechanical coupling. The dynamic permeability concept is built on observations that fracture aperture changes as a function of hydraulic head, as measured during in situ extensometer testing of individual fractures, and observed changes in bulk aquifer transmissivity as determined from earth tides during seasonal changes in hydraulic head, with higher transmissivity during periods of high hydraulic head, and lower transmissivity during periods of relatively lower hydraulic head. A final conceptual model is presented that captures both the static and dynamic properties of the aquifer. The workflow presented here demonstrates development of a conceptual framework for building numerical models of complexly fractured, foliated, crystalline rock aquifers that includes both a static model to describe the spatial distribution of permeability as a function of fracture type and foliation intensity and a dynamic model that describes how hydromechanical coupling impacts permeability magnitude as a function of hydraulic head fluctuation. This model captures important geologic controls on permeability magnitude, anisotropy, and transience and therefor offers potentially more reliable history matching and forecasts of different water management strategies, such as resource evaluation, well placement, permeability prediction, and evaluating remediation strategies. 相似文献
11.
12.
We present a novel pilot-point-based hydraulic tomography (HT) inversion procedure to delineate preferential flow paths and estimate hydraulic properties in a fractured aquifer. Our procedure considers a binary prior model developed using a randomized algorithm. The randomized algorithm involves discretizing the domain into grid cells, assigning a binary label to each cell, traversing the grid randomly, and choosing the optimal grid configuration cell-by-cell. This binary prior model is used to guide the placement of pilot points and to constrain aquifer parameters during pilot-point-based HT inversion. A two-dimensional fractured granite rock block was considered to test our methodology under controlled laboratory conditions. Multiple pumping tests were conducted at selected ports and the pressure responses were monitored. The pumping datasets thus obtained were preprocessed using median filters to remove random noise, and then analyzed using the proposed procedure. The proposed binary prior algorithm was implemented in C++ by supplying the forward groundwater model, HydroGeoSphere (HGS). Pilot-point-assisted HT inversion was performed using the parameter-estimation tool, coupled to HGS. The resulting parameter distributions were assessed by: (1) a visual comparison of the K- and Ss-tomograms with the known topology of the fractures and (2) comparing model predictions with measurements made at two validation ports that were not used in calibration. The performance assessment revealed that HT with the proposed randomized binary prior could be used to recover fracture-connectivity and to predict drawdowns in fractured aquifers with reasonable accuracy, when compared to a conventional pilot-point inversion scheme. 相似文献
13.
A new method is presented to search for hydraulically transmissive features in open boreholes in bedrock. A flexible borehole liner made of a watertight, nylon fabric is filled with water to create a constant driving head to evert (reverse of invert) the liner down the hole so that the liner pushes the borehole water out into transmissive fractures or other permeable features. The descent rate is governed by the bulk transmissivity of the remaining permeable features below the liner. Initially, the liner descent rate or velocity is a measure of transmissivity (T) of the entire hole. As the everting liner passes and seals each permeable feature, changes in the liner velocity indicate the position of each feature and an estimate of T using the Thiem equation for steady radial flow. This method has been performed in boreholes with diameters ranging from 96 to 330 mm. Profiling commonly takes a few hours in holes 200‐ to 300‐m long. After arrival of the liner at the bottom of the hole, the liner acts as a seal preventing borehole cross connection between transmissive features at different depths. Liner removal allows the hole to be used for other purposes. The T values determined using this method in a dolostone aquifer were found to be similar to the values from injection tests using conventional straddle packers. This method is not a replacement for straddle‐packer hydraulic testing of specific zones where greater accuracy is desired; however, it is effective and efficient for scanning entire holes for transmissive features. 相似文献
14.
The HySuf‐FEM code (Hydrodynamic of Subsurface Flow by Finite Element Method) is proposed in this article in order to estimate the spatial variability of the transmissivity values of the Berrechid aquifer (Morocco). The calibration of the model is based on the hydraulic head, hydraulic conductivity and recharge. Three numerical tests are used to validate the model and verify its convergence. The first test case consists in using the steady analytical solution of the Poisson equation. In the second, the model has been compared with the hydrogeological system which is characterized by an unconfined monolayer (isotropic layer) and computed by using PMWIN‐MODFLOW software. The third test case is based on the comparison between the results of HySuf‐FEM and the multiple cell balance method in the aquifer system with natural boundaries case. Good agreement between the Hydrodynamic of Subsurface Flow, the numerical tests and the spatial distribution of the thickening of the hydrogeological system is deduced from the analysis and the interpretations of hydrogeological wells. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
15.
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. 相似文献
16.
Actual pumping tests may involve continuously decreasing rates over a certain period of time, and the hydraulic conductivity (K) and specific storage (Ss) of the tested confined aquifer cannot be interpreted from the classical constant‐rate test model. In this study, we revisit the aquifer drawdown characteristics of a pumping test with an exponentially decreasing rate using the dimensionless analytical solution for such a variable‐rate model. The drawdown may decrease with time for a short period of time at intermediate pumping times for such pumping tests. A larger ratio of initial to final pumping rate and a smaller radial distance of the observation well will enhance the decreasing feature. A larger decay constant results in an earlier decrease, but it weakens the extent of such a decrease. Based on the proposed dimensionless transformation, we have proposed two graphical methods for estimating K and Ss of the tested aquifer. The first is a new type curve method that does not employ the well function as commonly done in standard type curve analysis. Another is a new analytic method that takes advantage of the decreasing features of aquifer drawdown during the intermediate pumping stage. We have demonstrated the applicability and robustness of the two new graphical methods for aquifer characterization through a synthetic pumping test. 相似文献
17.
This paper presents a geostatistical approach to multi-directional aquifer stimulation in order to better identify the transmissivity
field. Hydraulic head measurements, taken at a few locations but under a number of different steady-state flow conditions,
are used to estimate the transmissivity. Well installation is generally the most costly aspect of obtaining hydraulic head
measurements. Therefore, it is advantageous to obtain as many informative measurements from each sampling location as possible.
This can be achieved by hydraulically stimulating the aquifer through pumping, in order to set-up a variety of flow conditions.
We illustrate the method by applying it to a synthetic aquifer. The simulations provide evidence that a few sampling locations
may provide enough information to estimate the transmissivity field. Furthermore, the innovation of, or new information provided
by, each measurement can be examined by looking at the corresponding spline and sensitivity matrix. Estimates from multi-directional
stimulation are found to be clearly superior to estimates using data taken under one flow condition. We describe the geostatistical
methodology for using data from multi-directional simulations and address computational issues. 相似文献
18.
Determination of hydraulic head, H, as a function of spatial coordinates and time, in ground water flow is the basis for aquifer management and for prediction of contaminant transport. Several computer codes are available for this purpose. Spatial distribution of the transmissivity, T(x,y), is a required input to these codes. In most aquifers, T varies in an erratic manner, and it can be characterized statistically in terms of a few moments: the expected value, the variance, and the variogram. Knowledge of these moments, combined with a few measurements, permits one to estimate T at any point using geostatistical methods. In a review of transmissivity data from 19 unconsolidated aquifers, Hoeksema and Kitanidis (1985) identified two types of the logtransmissivity Y= ln(T) variations: correlated variations with variance sigma2Yc and correlation scale, I(Y), on the order of kilometers, and uncorrelated variations with variance sigma2Yn. Direct identification of the logtransmissivity variogram, Gamma(Y), from measurements is difficult because T data are generally scarce. However, many head measurements are commonly available. The aim of the paper is to introduce a methodology to identify the transmissivity variogram parameters (sigma2Yc, I(Y), and sigma2Yn) using head data in formations characterized by large logtransmissivity variance. The identification methodology uses a combination of precise numerical simulations (carried out using analytic element method) and a theoretical model. The main objective is to demonstrate the application of the methodology to a regional ground water flow in Eagle Valley basin in west-central Nevada for which abundant transmissivity and head measurements are available. 相似文献
19.
Forced and free oscillations of water level were recorded in the YuZ-5 well, Kamchatka due to the passage of seismic waves from the Sumatra-Andaman earthquake of December 26, 2004, M w = 9.3, hypocentral distance 8250 km. The greatest amplitude of water level oscillations, at least 5 cm, was observed during the onset of seismic surface waves with a typical period of 20–50 s. The total duration of the forced and free water level oscillations was about ten hours. The available theoretical models that describe oscillations of water level in a well due to seismic waves and rapid injection of water were used to estimate the transmissivity of the aquifer. The values obtained exceed by at least two orders of magnitude the transmissivity derived from pumping test measurements. A hypothesis was proposed to explain the temporary increase in aquifer transmissivity during the passage of seismic waves by invoking disturbances in the structure of the crack-pore space and a sharp increase in aquifer rock permeability. 相似文献
20.
We introduce a simple correction to coastal heads for constant‐density groundwater flow models that contain a coastal boundary, based on previous analytical solutions for interface flow. The results demonstrate that accurate discharge to the sea in confined aquifers can be obtained by direct application of Darcy's law (for constant‐density flow) if the coastal heads are corrected to ((α + 1)/α)hs ? B/2α, in which hs is the mean sea level above the aquifer base, B is the aquifer thickness, and α is the density factor. For unconfined aquifers, the coastal head should be assigned the value . The accuracy of using these corrections is demonstrated by consistency between constant‐density Darcy's solution and variable‐density flow numerical simulations. The errors introduced by adopting two previous approaches (i.e., no correction and using the equivalent fresh water head at the middle position of the aquifer to represent the hydraulic head at the coastal boundary) are evaluated. Sensitivity analysis shows that errors in discharge to the sea could be larger than 100% for typical coastal aquifer parameter ranges. The location of observation wells relative to the toe is a key factor controlling the estimation error, as it determines the relative aquifer length of constant‐density flow relative to variable‐density flow. The coastal head correction method introduced in this study facilitates the rapid and accurate estimation of the fresh water flux from a given hydraulic head measurement and allows for an improved representation of the coastal boundary condition in regional constant‐density groundwater flow models. 相似文献