首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 31 毫秒
1.
Deep basin aquifers are increasingly used in water‐stressed areas, though their potential for sustainable development is inhibited by overlying aquitards and limited recharge rates. Long open interval wells (LOIWs)—wells uncased through multiple hydrostratigraphic units—are present in many confined aquifer systems and can be an important mechanism for deep basin aquifers to receive flow across aquitards. LOIWs are a major control on flow in the deep Cambrian–Ordovician sandstone aquifers of the upper Midwest, USA, providing a source of artificial leakage from shallow bedrock aquifers and equilibrating head within the sandstone aquifers despite differential pumpage. Conceptualizing and quantifying this anthropogenic flow has long been a challenge for groundwater flow modellers, particularly on a regional scale. Synoptic measurements of active production wells and well completion data for northeast Illinois form the basis for a transient, head‐specified MODFLOW model that determines mass balance contributions to the region and estimates LOIW leakage to the aquifers. Using this insight, transient LOIW leakage was simulated using transiently changing KV zones in a traditional, Q‐specified MODFLOW‐USG model, a novel approach that allows the KV in a cell containing a LOIW to change transiently by use of the time‐variant materials (TVM) package. With this modification, we achieved a consistent calibration through time, averaging 19.9 m root mean squared error. This model indicates that artificial leakage via LOIWs contributed a minimum of 10–13% of total flow to the sandstone aquifers through the entire history of pumping, up to 50% of flow around 1930. Removal from storage exceeds 40% of flow during peak withdrawals, much of this flow sourced from units other than the primary sandstone aquifers via LOIWs. As such, understanding the timing and magnitude of LOIW leakage is essential for predicting future water availability in deep basin aquifers.  相似文献   

2.
This work evaluated the spatial variability and distribution of heterogeneous hydraulic conductivity (K) in the Choushui River alluvial fan in Taiwan, using ordinary kriging (OK) and mean and individual sequential Gaussian simulations (SGS). A baseline flow model constructed by upscaling parameters was inversely calibrated to determine the pumping and recharge rates. Simulated heads using different K realizations were then compared with historically measured heads. A global/local simulated error between simulated and measured heads was analysed to assess the different spatial variabilities of various estimated K distributions. The results of a MODFLOW simulation indicate that the OK realization had the smallest sum of absolute mean simulation errors (SAMSE) and the SGS realizations preserved the spatial variability of the measured K fields. Moreover, the SAMSE increases as the spatial variability of the K field increases. The OK realization yields small local simulation errors in the measured K field of moderate magnitude, whereas the SGS realizations have small local simulation errors in the measured K fields, with high and low values. The OK realization of K can be applied to perform a deterministic inverse calibration. The mean SGS method is suggested for constructing a K field when the application focuses on extreme values of estimated parameters and small calibration errors, such as in a simulation of contaminant transport in heterogeneous aquifers. The individual SGS realization is useful in stochastically assessing the spatial uncertainty of highly heterogeneous aquifers. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

3.
The hydraulic properties of aquitards are not easily obtained because monitoring wells are usually installed in aquifers for groundwater resources management. Earthquake‐induced crust stress (strain) triggers groundwater level variations over a short period of time in a large area. These groundwater anomalies can be used to investigate aquifer systems. This study uses a poroelastic model to fit the postseismic variations of groundwater level triggered by the Chi‐Chi earthquake to evaluate the hydraulic properties of aquitards in the Jhoushuei River alluvial fan (JRAF), Taiwan. Six of the adopted eight wells with depths of 70 to 130 m showed good agreement with the recovery theory. The mean hydraulic conductivities (K) of the aquifers for the eight wells are 1.62 × 10?4 to 9.06 × 10?4 m/s, and the thicknesses are 18.8 to 46.1 m. The thicknesses of the aquitards are 11.3 to 42.0 m. Under the isotropic assumption for K, the estimated values of K for the aquitards are 3.0 × 10?8 to 2.1 × 10?6 m/s, corresponding to a silty medium. The results match the values obtained for the geological material of the drilling core and those reported in previous studies. The estimated values were combined with those given in previous studies to determine the distribution of K in the first two aquitards in the JRAF. The distribution patterns of the aquitards reflect the sedimentary environments and fit the geological material. The proposed technique can be used to evaluate the K value of aquitards using inverse methods. The inversion results can be used in hydrogeological analyses, contaminant modeling, and subsidence evaluation.  相似文献   

4.
A process-based methodology was used to compare the vulnerability of public supply wells tapping seven study areas in four hydrologically distinct regional aquifers to volatile organic compound (VOC) contamination. This method considers (1) contributing areas and travel times of groundwater flowpaths converging at individual supply wells, (2) the oxic and/or anoxic conditions encountered along each flowpath, and (3) the combined effects of hydrodynamic dispersion and contaminant- and oxic/anoxic-specific biodegradation. Contributing areas and travel times were assessed using particle tracks generated from calibrated regional groundwater flow models. These results were then used to estimate VOC concentrations relative to an unspecified initial concentration (C/C0) at individual public supply wells. The results show that the vulnerability of public supply wells to VOC contamination varies widely between different regional aquifers. Low-recharge rates, long travel times, and the predominantly oxic conditions characteristic of Basin and Range aquifers in the western United States leads to lower vulnerability to VOCs, particularly to petroleum hydrocarbons such as benzene and toluene. On the other hand, high recharge rates and short residence times characteristic of the glacial aquifers of the eastern United States leads to greater vulnerability to VOCs. These differences lead to distinct patterns of C/C0 values estimated for public supply wells characteristic of each aquifer, information that can be used by resource managers to develop monitoring plans based on relative vulnerability, to locate new public supply wells, or to make land-use management decisions.  相似文献   

5.
为了研究抚仙湖紫外辐射(UVR)和光合有效辐射(PAR)衰减的时空特征及其与有色可溶性有机物(CDOM)、悬浮物(SS)、浮游植物(叶绿素a表征)等因子的关系,于2014年10月(秋季)、2015年1月(冬季)开展现场调查,结果显示:秋季不同波长(段)的漫射衰减系数Kd(305)、Kd(340)和Kd(PAR)分别为1.27±0.12、0.68±0.11和0.32±0.13 m-1,冬季分别为1.13±0.10、0.63±0.07和0.36±0.07 m-1;秋季CDOM的不同波长吸收系数ag(254)、ag(305)和ag(340)分别为4.09±0.26、1.18±0.09和0.57±0.05 m-1,冬季分别为2.95±0.24、0.61±0.11和0.11±0.07 m-1,秋季ag(254)、ag(305)和ag(340)显著高于冬季;秋季Kd(305)显著大于冬季,这与秋季(雨季)较高的CDOM丰度、浮游植物生物量(及SS浓度)有关.秋季ag(305)/Kd(305)、ag(340)/Kd(340)均显著高于冬季;秋季及秋冬季整体而言,ag(254)与Kd(305)、Kd(340)呈显著正相关,各多元逐步回归方程中均包含ag(254),说明CDOM吸收对UVR的衰减有重要贡献.空间差异方面,秋季北部的ag(254)、Kd(305)和Kd(340)显著高于南部,冬季南北部无明显差异,或与雨旱季北岸河流输入的CDOM和SS的情况有关.此外,浮游植物对UV-B衰减的影响和SS(与CDOM的交互作用)对UV-A衰减的影响更在于季节变化方面,而影响UVR、PAR衰减的各因子的相对贡献有待进一步量化.  相似文献   

6.
Aquifer hydraulic parameters are commonly inferred from constant-rate pumping tests, while variable pumping rates are frequently encountered in actual field conditions. In this study, we propose a generally applicable dimensionless form of the analytical solution for variable-rate pumping tests in confined aquifers. In particular, we adopt a piecewise-linear fitting of variable pumping rates and propose a new type-curve method for estimating the hydraulic conductivity (K ) and specific storage (Ss ) of the investigated confined aquifer. For each test, a series of type curves, which depend on the variable pumping rates, the location of observation wells and the introduced first dimensionless inflection time, need to be provided for matching the observed drawdown data on a log-log graph. We first demonstrate the applicability and robustness of this method through a synthetic pumping test. Subsequently, we apply this method to analyze drawdown data from four pumping tests conducted within a multilayered aquifer/aquitard system in Wuxi city, Jiangsu Province, China. The parameter estimates are then compared with those reported by PEST. The K and Ss values estimated by the new type-curve method are found to be quite close to PEST-based estimates. Parameter estimation results demonstrate the difference in K and Ss values between observation wells. The difference could be attributed to the spatial heterogeneity in K and Ss . A future research topic may focus on the characterization of K and Ss heterogeneity with the currently available drawdown data from variable-rate pumping tests.  相似文献   

7.
Mathematical relationships have been developed for reaeration rate coefficient (Ka) by various researchers. These relationships have a number of variables such as depth, velocity, width, slope, Froud number, molecular diffusion coefficient, kinematic viscosity and the gas‐transfer Reynolds number. From these variables, 29 relations have been developed and divided into four groups. To evaluate their predictive capability for highly variable flow rivers receiving high pollution loads form large cities, these relationships have been used to model dissolved oxygen (DO) in the River Ravi. Such rivers are either saturated with DO during high flows or anaerobic during critical low‐flow conditions. The evaluation is based on the agreement between model DO values calculated using Ka obtained from the available equations and the measured DO concentrations in the river samples in terms of sum of square of residuals (SSR) and coefficient of determination (R2). It has been found that in general, the group of equations containing depth and velocity as the only two variables affecting Ka performed better than the equations in other groups as reflected by lower SSR and higher R2 values. The study results also reveal that the turbulence‐based reaeration rate coefficient equation containing additional variables also resulted in close agreement between DO model results and the measured values. The study results identify the most important parameters affecting the reaeration rate coefficient and the suitability of various Ka relationships as well for rivers with highly variable flows. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

8.
Flow‐through column experiments were carried out to investigate the influence of pH on the sorption of three phenols (2‐methyl‐4, 6‐dinitrophenol, 2, 4, 6‐trichlorophenol, pentachlorophenol) onto a natural sandy aquifer material collected from a bank filtration site of River Elbe, Germany. For the phenols investigated, an increase in sorption (retardation) with decreasing pH is observed indicating a stronger sorption of the neutral species in comparison to that of the anions formed by dissociation. The anions of 2‐methyl‐4, 6‐dinitrophenol and 2, 4, 6‐trichlorophenol do not show significant sorption. On the contrary, pentachlorophenol showed sorption not only in neutral form but also in ionic form significantly which should be taken into account while assessing the fate and transport of such compound. A linear model based on the degree of protonation (calculated from pH and pKa) can be used to resolve the apparent (observed) sorption coefficient (Kd, app) into its neutral (Kd, n) and ionised (Kd, i) components. Knowing pKa, Kd, n, and Kd, i the apparent sorption coefficient for pH values other than experimentally investigated can be predicted.  相似文献   

9.
In granite aquifers, fractures can provide both storage volume and conduits for groundwater. Characterization of fracture hydraulic conductivity (K) in such aquifers is important for predicting flow rate and calibrating models. Nuclear magnetic resonance (NMR) well logging is a method to quickly obtain near-borehole hydraulic conductivity (i.e., KNMR) at high-vertical resolution. On the other hand, FLUTe flexible liner technology can produce a K profile at comparable resolution but requires a fluid driving force between borehole and formation. For three boreholes completed in a fractured granite, we jointly interpreted logging NMR data and FLUTe K estimates to calibrate an empirical equation for translating borehole NMR data to K estimates. For over 90% of the depth intervals investigated from these boreholes, the estimated KNMR are within one order of magnitude of KFLUTe. The empirical parameters obtained from calibrating the NMR data suggest that “intermediate diffusion” and/or “slow diffusion” during the NMR relaxation time may occur in the flowing fractures when hydraulic aperture are sufficiently large. For each borehole, “intermediate diffusion” dominates the relaxation time, therefore assuming “fast diffusion” in the interpretation of NMR data from fractured rock may lead to inaccurate KNMR estimates. We also compare calibrations using inexpensive slug tests that suggest reliable KNMR estimates for fractured rock may be achieved using limited calibration against borehole hydraulic measurements.  相似文献   

10.
Characterization of hydraulic conductivity (K) in aquifers is critical for evaluation, management, and remediation of groundwater resources. While estimates of K have been traditionally obtained using hydraulic tests over discrete intervals in wells, geophysical measurements are emerging as an alternative way to estimate this parameter. Nuclear magnetic resonance (NMR) logging, a technology once largely applied to characterization of deep consolidated rock petroleum reservoirs, is beginning to see use in near‐surface unconsolidated aquifers. Using a well‐known rock physics relationship—the Schlumberger Doll Research (SDR) equation—K and porosity can be estimated from NMR water content and relaxation time. Calibration of SDR parameters is necessary for this transformation because NMR relaxation properties are, in part, a function of magnetic mineralization and pore space geometry, which are locally variable quantities. Here, we present a statistically based method for calibrating SDR parameters that establishes a range for the estimated parameters and simultaneously estimates the uncertainty of the resulting K values. We used co‐located logging NMR and direct K measurements in an unconsolidated fluvial aquifer in Lawrence, Kansas, USA to demonstrate that K can be estimated using logging NMR to a similar level of uncertainty as with traditional direct hydraulic measurements in unconsolidated sediments under field conditions. Results of this study provide a benchmark for future calibrations of NMR to obtain K in unconsolidated sediments and suggest a method for evaluating uncertainty in both K and SDR parameter values.  相似文献   

11.
A Python Script to Compute Isochrones for MODFLOW   总被引：1，自引：0，他引：1       下载免费PDF全文
  相似文献   

12.
In this methods note, we present a simple analytical formula to quantify the steady‐state leakage flux (Q) over a perforated aquitard. The flux depends on the aquitard thickness (D), the radius of the perforation (R), the hydraulic conductivity of the material inside the perforation (kfill), the conductivities of the overlying and underlying aquifers (k1 and k2, respectively), and the head difference between the two aquifers (ΔH): This equation assumes an aquitard separating two homogeneous and infinite aquifers (R ? aquifer thickness) in which radial flow to and from the perforation occurs, with no other recharge or discharge boundaries near the perforation. The flux through a perforation in a hypothetical case study with D = 10 m, k1 = 10 m/d, k2 = 20 m/d, R = 0.072 m, and ΔH = 1 m ranges between less than 1 mL/d if the hole is backfilled with bentonite (kfill = 10?4 m/d), to several liters per day if the perforation is backfilled with sand from the overlying aquifer (kfill = 10 m/d), to several m3/d if the perforation forms an open conduit (kfill = 105 m/d). The leakage fluxes calculated with this model agree well with those calculated using a numerical model (MODFLOW).  相似文献   

13.
The vertical hydraulic conductivity (Kv) of a stream or lake sediment is often determined in the field using standpipe tests. Calculation of Kv is based on the assumption that the hydraulic head in the pipe is equal to that of the stream or lake stage. In this work, a modified equation for Kv is developed for the standpipe test which is applicable when this assumption is not valid. The equation involves not only the hydraulic head at different times but also the difference in the hydraulic head (a) between the groundwater level and river stage. The effects of certain factors on Kv, such as the ratio of the hydraulic head at different times (h1/h2), the difference a, and the initial water table height (h0), are also discussed. The results show that when h1/h2 is constant, the relative error (Er) in Kv increases with the ratio a/h2. Furthermore, if a/h2 < 0.05, then for any value of h1/h2, Er is less than 5% using the modified equation. Also, if a/h2 is large, hydraulic head readings with larger h1/h2 ratios must be used to avoid large Er values. The results of a field test also indicate that the error in Kv decreases as the value of h0 increases. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

14.
The time required at a field site to obtain a few measurements of saturated hydraulic conductivity (Ks) will allow for many measurements of soil air permeability (ka). This study investigates if ka measured in situ (ka, in situ) can be a substitute for measurement of Ks in relation to infiltration and surface runoff modelling. Measurements of ka, in situ were carried out in two small agricultural catchments. A spatial correlation of the log‐transformed values existed having a range of approximately 100 m. A predictive relationship between Ks and ka measured on 100‐cm3 soil samples in the laboratory was derived for one of the field slopes and showed good agreement with an earlier suggested predictive Ks–ka relationship. In situ measurements of Ks and ka suggested that the predictive relationships also could be used at larger scale. The Ks–ka relationships together with the ka, in situ data were applied in a distributed surface runoff (DSR) model, simulating a high‐intensity rainfall event. The DSR simulation results were highly dependent on whether the geometric average of ka, in situ or kriged values of ka, in situ was used as model input. When increasing the resolution of Ks in the DSR model, a limit of 30–40 m was found for both field slopes. Below this limit, the simulated runoff and hydrograph peaks were independent of resolution scale. If only a few randomly chosen values of Ks were used to represent the spatial variation within the field slope, very large deviations in repeated DSR simulation results were obtained, both with respect to peak height and hydrograph shape. In contrast, when using many predicted Ks values based on a Ks–ka relationship and measured ka, in situ data, the DSR model generally captured the correct hydrograph shape although simulations were sensitive to the chosen Ks–ka relationship. As massive measurement efforts normally will be required to obtain a satisfactory representation of the spatial variability in Ks, the use of ka, in situ to assess spatial variability in Ks appears a promising alternative. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

15.
Abstract

Time-domain reflectometry (TDR) is an electromagnetic technique for measurements of water and solute transport in soils. The relationship between the TDR-measured dielectric constant (Ka ) and bulk soil electrical conductivity ([sgrave]a) to water content (θW) and solute concentration is difficult to describe physically due to the complex dielectric response of wet soil. This has led to the development of mostly empirical calibration models. In the present study, artificial neural networks (ANNs) are utilized for calculations of θw and soil solution electrical conductivity ([sgrave]w) from TDR-measured Ka and [sgrave]a in sand. The ANN model performance is compared to other existing models. The results show that the ANN performs consistently better than all other models, suggesting the suitability of ANNs for accurate TDR calibrations.  相似文献   

16.
Abstract

Abstract Recently, substantial progress has been made in detection and observation of non-aqueous phase liquids (NAPLs) in the subsurface using different experimental techniques. However, there is still a lack of appropriate direct methods to measure the saturation of NAPL (θNAPL). This paper provides a guide for estimating θNAPL and water content (θ w ) in unsaturated and saturated sand based on direct measurements of soil dielectric constant (Ka ) and electrical conductivity (σ a ) using time domain reflectometry (TDR). The results show that the previously used dielectric mixing model fails to predict θNAPL in the case of a four-phase system. A new methodology is suggested and exemplified by showing that the measured Ka gives accurate estimation of θNAPL for a three-phase system while in a four-phase system, both θ a and Ka need to be measured. The results show that using the suggested methodology, accurate predictions of θ w (R 2 = 0.9998) and θNAPL lower than 0.20 m3 m-3 (average R 2 = 0.9756) are possible.  相似文献   

17.
Pressure–volume–temperature relations have been measured to 32 GPa and 2073 K for natural magnesite (Mg0.975Fe0.015Mn0.006Ca0.004CO3) using synchrotron X-ray diffraction with a multianvil apparatus at the SPring-8 facility. A least-squares fit of the room-temperature compression data to a third-order Birch–Murnaghan equation of state (EOS) yielded K0 = 97.1 ± 0.5 GPa and K′ = 5.44 ± 0.07, with fixed V0 = 279.55 ± 0.02 Å3. Further analysis of the high-temperature compression data yielded the temperature derivative of the bulk modulus (∂KT/∂T)P = −0.013 ± 0.001 GPa/K and zero-pressure thermal expansion α = a0 + a1T with a0 = 4.03 (7) × 10−5 K−1 and a1 = 0.49 (10) × 10−8 K−2. The Anderson–Grüneisen parameter is estimated to be δT = 3.3. The analysis of axial compressibility and thermal expansivity indicates that the c-axis is over three times more compressible (KTc = 47 ± 1 GPa) than the a-axis (KTc = 157 ± 1 GPa), whereas the thermal expansion of the c-axis (a0 = 6.8 (2) × 10−5 K−1 and a1 = 2.2 (4) × 10−8 K−2) is greater than that of the a-axis (a0 = 2.7 (4) × 10−5 K−1 and a1 = −0.2 (2) × 10−8 K−2). The present thermal EOS enables us to accurately calculate the density of magnesite to the deep mantle conditions. Decarbonation of a subducting oceanic crust containing 2 wt.% magnesite would result in a 0.6% density reduction at 30 GPa and 1273 K. Using the new EOS parameters we performed thermodynamic calculations for magnesite decarbonation reactions at pressures to 20 GPa. We also estimated stability of magnesite-bearing assemblages in the lower mantle.  相似文献   

18.
A new MODFLOW package (Nonlinear Flow Process; NLFP) simulating nonlinear flow following the Forchheimer equation was developed and implemented in MODLFOW‐2005. The method is based on an iterative modification of the conductance calculated and used by MODFLOW to obtain an effective Forchheimer conductance. The package is compatible with the different layer types, boundary conditions, and solvers as well as the wetting capability of MODFLOW. The correct implementation is demonstrated using four different benchmark scenarios for which analytical solutions are available. A scenario considering transient flow in a more realistic setting and a larger model domain with a higher number of cells demonstrates that NLFP performs well under more complex conditions, although it converges moderately slower than the standard MODFLOW depending on the nonlinearity of flow. Thus, this new tool opens a field of opportunities to groundwater flow simulation with MODFLOW, especially for core sample simulation or vuggy karstified aquifers as well as for nonlinear flow in the vicinity of pumping wells.  相似文献   

19.
  相似文献   

20.
Streambed horizontal hydraulic conductivity (Kh) has a substantial role in controlling exchanges between stream water and groundwater. We propose a new approach for determining Kh of the shallow streambed sediments. Undisturbed sediment samples were collected using tubes that were horizontally driven into streambeds. The sediment columns were analysed using a permeameter test (PT) on site. This new test approach minimizes uncertainties due to vertical flow in the vicinity of test tube and stream stage fluctuations in the computation of the Kh values. Ninety‐eight PTs using the new approach were conducted at eight sites in four tributaries of the Platte River, east‐central Nebraska, USA. The Kh values were compared with the nondirectional hydraulic conductivity values (Kg) determined from 12 empirical grain‐size analysis methods. The grain‐size analysis methods used the same sediment samples as Kh tests. Only two methods, the Terzaghi and Shepherd methods, yielded Kg values close to the Kh values. Although the Sauerbrei method produced a value relatively closer to Kh than other nine grain‐size analysis methods, the values from this method were not as reliable as the methods of Terzaghi and Shepherd due to the inconsistent fluctuation of the average estimates at each of the test sites. The Zunker, Zamarin, Hazen, Beyer, and Kozeny methods overestimated Kh, while the Slichter, US Bureau of Reclamation (USBR), Harleman, and Alyamani and Sen methods underestimated Kh. Any of these specific grain‐size methods might yield good estimates of streambed Kh at some sites, but give poor estimates at other sites, indicating that the relationship between Kg and Kh is significantly site dependent in our study. Copyright © 2011 John Wiley & Sons, Ltd.  相似文献   

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

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