首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 31 毫秒
1.
Cross-borehole flowmeter tests have been proposed as an efficient method to investigate preferential flowpaths in heterogeneous aquifers, which is a major task in the characterization of fractured aquifers. Cross-borehole flowmeter tests are based on the idea that changing the pumping conditions in a given aquifer will modify the hydraulic head distribution in large-scale flowpaths, producing measurable changes in the vertical flow profiles in observation boreholes. However, inversion of flow measurements to derive flowpath geometry and connectivity and to characterize their hydraulic properties is still a subject of research. In this study, we propose a framework for cross-borehole flowmeter test interpretation that is based on a two-scale conceptual model: discrete fractures at the borehole scale and zones of interconnected fractures at the aquifer scale. We propose that the two problems may be solved independently. The first inverse problem consists of estimating the hydraulic head variations that drive the transient borehole flow observed in the cross-borehole flowmeter experiments. The second inverse problem is related to estimating the geometry and hydraulic properties of large-scale flowpaths in the region between pumping and observation wells that are compatible with the head variations deduced from the first problem. To solve the borehole-scale problem, we treat the transient flow data as a series of quasi-steady flow conditions and solve for the hydraulic head changes in individual fractures required to produce these data. The consistency of the method is verified using field experiments performed in a fractured-rock aquifer.  相似文献   

2.
Zheng C  Gorelick SM 《Ground water》2003,41(2):142-155
Several recent studies at the Macrodispersion Experiment (MADE) site in Columbus, Mississippi, have indicated that the relative preferential flowpaths and flow barriers resulting from decimeter-scale aquifer heterogeneities appear to have a dominant effect on plume-scale solute transport. Numerical experiments are thus conducted in this study to explore the key characteristics of solute transport in two-dimensional flow fields influenced by decimeter-scale preferential flowpaths. A hypothetical but geologically plausible network of 10 cm wide channels of high hydraulic conductivity is used to represent the relative preferential flowpaths embedded in an otherwise homogeneous aquifer. When the hydraulic conductivity in the channels is 100 times greater than that in the remaining portion of the aquifer, the calculated concentration distributions under three source configurations all exhibit highly asymmetrical, non-Gaussian patterns. These patterns, with peak concentrations close to the source and extensive spreading downgradient, resemble that observed at the MADE site tracer tests. When the contrast between the channel and nonchannel hydraulic conductivities is reduced to 30:1 from 100:1, the calculated mass distribution curve starts to approach a Gaussian one with the peak concentration near the central portion of the plume. Additional analysis based on a field-scale model demonstrates that the existence of decimeter-scale preferential flowpaths can have potentially far-reaching implications for ground water remediation. Failure to account for them in numerical simulation could lead to overestimation of the effectiveness of the remedial measure under consideration.  相似文献   

3.
《水文科学杂志》2013,58(6):1106-1124
Abstract

Natural soils very often contain micro- and macropores, having different hydraulic properties. At the macroscopic scale, the unsaturated flow in such soils can be described with various models, depending on the hydraulic diffusivity ratio of the components and the connectivity of the most conductive component. Three macroscopic models recently derived by the homogenization method are discussed. The limit passages between the models are studied. A unified model suitable for the entire range of the hydraulic diffusivity ratio is proposed. A numerical example shows the application of the model to macroscopically one-dimensional infiltration in a porous medium containing inclusions. A parametric study for varying conductivity (diffusivity) ratio is performed.  相似文献   

4.
Litang Hu 《水文科学杂志》2013,58(15):2694-2703
ABSTRACT

Estimation of hydraulic properties in the field is usually small-scale and not cost-effective. This paper proposes an innovative method for estimating hydraulic diffusivity at regional scale. Monthly groundwater storage change over the period from 2003 to 2013 is first estimated from GRACE-derived terrestrial water storage (TWS). Assuming that the aquifer system is unconfined and the hydraulic properties are uniform in a geographical cell, the water balance principle and Darcy’s law are used to establish a relation between groundwater storage and hydraulic diffusivity. The value of hydraulic diffusivity is then adjusted using the generalized least squares and linear correlation method. This GRACE-derived hydraulic diffusivity estimation method, or GHDE method for short, is first verified with a hypothetical case and then applied in the Beishan area with available field-measured hydraulic conductivity data. The hypothetical case study demonstrates that the method works perfectly if the TWS data are error free. The Beishan case study illustrates that the estimated hydraulic diffusivities using the GHDE method correlate reasonably well with field test results, suggesting that this method is applicable. The accuracy of this method is constrained by the resolution of the GRACE-derived TWS data and is most suitable for very large scale groundwater problems due to the current accuracy of the GRACE data.
EDITOR A. Castellarin ASSOCIATE EDITOR N. Verhoest  相似文献   

5.
Most published solutions for aquifer responses to ocean tides focus on the one-sided attenuation of the signal as it propagates inland. However, island aquifers experience periodic forcing from the entire coast, which can lead to integrated effects of different tidal signals, especially on narrow high-permeability islands. In general, studies disregard a potential time lag as the tidal wave sweeps around the island. We present a one-dimensional analytical solution to the ground water flow equation subject to asynchronous and asymmetric oscillating head conditions on opposite boundaries and test it on data from an unconfined volcanic aquifer in Maui. The solution considers sediment-damping effects at the coastline. The response of Maui Aquifers indicate that water table elevations near the center of the aquifer are influenced by a combination of tides from opposite coasts. A better match between the observed ground water head and the theoretical response can be obtained with the proposed dual-tide solution than with single-sided solutions. Hydraulic diffusivity was estimated to be 2.3 × 107 m2/d. This translates into a hydraulic conductivity of 500 m/d, assuming a specific yield of 0.04 and an aquifer thickness of 1.8 km. A numerical experiment confirmed the hydraulic diffusivity value and showed that the y -intercepts of the modal attenuation and phase differences estimated by regression can approximate damping factors caused by low-permeability units at the boundary.  相似文献   

6.
The study demonstrates spectral relationships in the time–frequency domain for one‐dimensional groundwater flow in aquifers bounded by fluctuating boundaries. By nature, the solutions of spectral equations are non‐linear complex functions. To determine hydraulic diffusivity in the governing equations, it is required that the data are collected from the spectra of water levels at the fluctuating boundaries and observation wells. Hydraulic diffusivity thus can be obtained by an iterative inverse approach. This paper presents an application in Pingtung County of Taiwan to determine the hydraulic diffusivity of a sandy aquifer under confined conditions. Spectral density function of water level obtained from tidal boundaries and observation wells are used to approximate hydraulic diffusivity, which yields an averaged value of 1·26 × 106 m2/h. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

7.
Fractured rock aquifers cover much of Earth's surface and are important mountain sites for groundwater recharge but are poorly understood. To investigate groundwater systematics of a fractured-dominated aquifer in Baja California Sur, Mexico, we examined the spatial patterns of aquifer recharge and connectivity using the geochemistry of springs. We evaluate a range of geochemical data within the context of two endmember hypotheses describing spatial recharge patterns and fracture connectivity. Hypothesis 1 is that the aquifer system is segmented, and springs are fed by local recharge. Hypothesis 2 is that the aquifer system is well connected, with dominant recharge occurring in the higher elevations. The study site is a small <15 km2 catchment. Thirty-four distinct springs and two wells were identified in the study area, and 24 of these sites were sampled for geochemical analyses along an elevation gradient and canyon transect. These analyses included major ion composition, trace element and strontium isotopes, δ18O and δ2H isotopes, radiocarbon, and tritium. δ18O and δ2H isotopes suggest that the precipitation feeding the groundwater system has at least two distinct sources. Carbon isotopes showed a change along the canyon transect, suggesting that shorter flowpaths feed springs in the top of the transect, and longer flowpaths discharge near the bottom. Geochemical interpretations support a combination of the two proposed hypotheses. Understanding of the connectivity and provenance of these springs is significant as they are the primary source of water for the communities that inhabit this region and may be impacted by changes in recharge and use.  相似文献   

8.
运用加卸载响应比方法,计算和分析了山西省2008年以来数字化地下水位固体潮观测数据,并选取大于2倍加卸载响应比值作为异常判定指标,将山西及周边地区发生的ML≥4.0地震的参数与井水位固体潮响应比异常变化对比分析。结果表明,山西地区9口数字水位观测井的地下水位固体潮加卸载响应比异常变化与附近的地震活动存在不同程度的对应关系。从而说明,利用井水位固体潮观测的加卸载响应比是一种能有效识别和提取地震前兆的方法。  相似文献   

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

10.
Analyses of independent laboratory- and field-scale measurements from two sites on Sapelo Island, Georgia reveal heterogeneity in hydraulic parameters across the upland–estuary interface. Regardless of the method used (short-duration pumping tests, amplitude attenuation of tidal pumping data, sediment grain size distributions, and falling head permeameter tests), we obtain hydraulic conductivity of 10−4 m s−1 for the fine-grained, well-sorted, clean sands that make up the upland areas. Proximal to the upland–estuary boundary, the tidal pumping analyses and permeameter tests suggest that hydraulic conductivities decrease by more than two orders of magnitude, a result consistent with the presence of a clogging layer. Such a clogging layer may arise due to a variety of physical, chemical, or biological processes. The extent and orientation of the layers of reduced hydraulic conductivity near the upland–estuary boundary influence the nature of the aquifer's response to tidal forcing. Where the lower conductivity layer forms a relatively flat creek bank, tidal pumping produces a primarily mechanical response in the adjacent aquifer. Where the creek bank is nearly vertical, there is a more direct hydraulic connection between the tidal creek and the adjacent aquifer. The clogging layer likely contributes to the development of complicated flow pathways across the upland–estuary boundary. Effective flow paths calculated from tidal pumping data terminate within the marsh, beyond the boundary of the upland aquifer, suggesting a diffuse regime of groundwater discharge in the marsh. We postulate that, in many settings, submarsh flow may be as important as seepage faces for groundwater discharge into the marsh–estuary complex.  相似文献   

11.
According to the fact that the Xinfengjiang reservoir earthquakes are caused mainly by water seepage, this paper using the data ofM s?2. 0 earthquakes, studies the hydraulic diffusivity of the mainshock zone by the expansion of the distribution area of epicenters. It is indicated thatin-situ hydraulic diffusivity during the preshock activity of the Xinfengjiang reservoir region was about 6. 2 m2/s. However, after the main shock, thein-situ hydraulic diffusivity in the main shock region increased by fifty percent, that is to say, to 9. 7 m2/s. During the long period after main shock occurrence thein-situ hydraulic diffusivity was affected by significant anisotropy of the medium and fluctuation of water level. No regularity can be found. In addition, we compare the diffusivity found by experiment with rock samples collected with thein-situ hydraulic diffusivity estimated. It is shown that the diffusivity of rock samples with fractures is about the same as the diffusivity estimatedin-situ. However, the diffusivity of whole rock samples is 3 orders of magnitudes smaller than that estimatedin-situ. Finally, we discuss the limits on the method by the expansion of distribution area of epicenters in the study of reservoir induced seismicity.  相似文献   

12.
Tidal and wind-driven surface currents in the German Bight between shallow mudflats of the North Frisian islands and the island of Helgoland are studied using coastal high-frequency radar (HFR) observations and hindcasts from a primitive equation numerical model. The setup of the observational system is described, and estimates of expected measurement errors are given. A quantitative comparison of numerical model results and observations is performed. The dominant tidal components are extracted from the two data sources using tidal harmonic analysis and the corresponding tidal ellipses are defined. Results show that the spatial patterns of different tidal ellipse parameters are consistent in the two data sets. Model sensitivity studies with constant and variable salinity and temperature distributions are used to study density-related mechanisms of circulation. Furthermore, the role of the surface wind field in driving the German Bight circulation is investigated using the complex correlation between wind and surface current vectors. The observed change of the respective correlation patterns from the coastal to open ocean is shown to be due to a combination of density effects, the coastline and topography. The overall conclusion is that HFR observations resolve the small-scale and rapidly evolving characteristics of coastal currents well in the studied area and could present an important component for regional operational oceanography when combined with numerical modelling. Some unresolved issues associated with the complex circulation and large instability of circulation in front of the Elbe River Estuary justify further considerations of this area using dedicated surveys and modelling efforts.  相似文献   

13.
According to the fact that the Xinfengjiang reservoir earthquakes are caused mainly by water seepage, this paper using the data ofM s⩾2. 0 earthquakes, studies the hydraulic diffusivity of the mainshock zone by the expansion of the distribution area of epicenters. It is indicated thatin-situ hydraulic diffusivity during the preshock activity of the Xinfengjiang reservoir region was about 6. 2 m2/s. However, after the main shock, thein-situ hydraulic diffusivity in the main shock region increased by fifty percent, that is to say, to 9. 7 m2/s. During the long period after main shock occurrence thein-situ hydraulic diffusivity was affected by significant anisotropy of the medium and fluctuation of water level. No regularity can be found. In addition, we compare the diffusivity found by experiment with rock samples collected with thein-situ hydraulic diffusivity estimated. It is shown that the diffusivity of rock samples with fractures is about the same as the diffusivity estimatedin-situ. However, the diffusivity of whole rock samples is 3 orders of magnitudes smaller than that estimatedin-situ. Finally, we discuss the limits on the method by the expansion of distribution area of epicenters in the study of reservoir induced seismicity. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,13, 364–371, 1991. This subject is sponsored by the Chinese Joint Seismological Science Foundation.  相似文献   

14.
Unstructured mesh models can resolve the model domain with a variable and very fine mesh resolution. Nevertheless, tuning the model setup is still required (for example because of parametrized sub-grid processes). Adjoint models are commonly used to calculate sensitivities of ocean models and optimize their parameters so that better agreement is achieved between model simulations and observations. One major obstacle in developing an adjoint model is the need to update the reverse code after each modification of the forward code, which is not always straightforward. Automatic differentiation is a tool to generate the adjoint model code without user input. So far this method has mainly been used for structured mesh ocean models. We present here an unstructured mesh, adjoint, tidal model using this technique, and discuss the sensitivities of the misfit between simulated and observed elevations with respect to open boundary values, the bottom friction coefficient and the bottom topography. The forward model simulates tides on the European Continental Shelf and we show that the tidal model dynamics in the adjoint simulations can be used to define regions, where parameters or mesh has to be optimized. We analyze the dependence of the sensitivities on the wave type and mesh resolution to specify whether the model misfit originates from physical or numerical model deficiencies. In the sensitivity patterns, it is possible to identify islands not resolved in the mesh. We propose to refine the mesh prior to the parameter optimization.  相似文献   

15.
Monitoring runoff generation processes in the field is a prerequisite for developing conceptual hydrological models and theories. At the same time, our perception of hydrological processes strongly depends on the spatial and temporal scale of observation. Therefore, the aim of this study is to investigate interactions between runoff generation processes of different spatial scales (plot scale, hillslope scale, and headwater scale). Different runoff generation processes of three hillslopes with similar topography, geology and soil properties, but differences in vegetation cover (grassland, coniferous forest, and mixed forest) within a small v‐shaped headwater were measured: water table dynamics in wells with high spatial and temporal resolution, subsurface flow (SSF) of three 10 m wide trenches at the bottom of the hillslopes subdivided into two trench sections each, overland flow at the plot scale, and catchment runoff. Bachmair et al. ( 2012 ) found a high spatial variability of water table dynamics at the plot scale. In this study, we investigate the representativity of SSF observations at the plot scale versus the hillslope scale and vice versa, and the linkage between hillslope dynamics (SSF and overland flow) and streamflow. Distinct differences in total SSF within each 10 m wide trench confirm the high spatial variability of the water table dynamics. The representativity of plot scale observations for hillslope scale SSF strongly depends on whether or not wells capture spatially variable flowpaths. At the grassland hillslope, subsurface flowpaths are not captured by our relatively densely spaced wells (3 m), despite a similar trench flow response to the coniferous forest hillslope. Regarding the linkage between hillslope dynamics and catchment runoff, we found an intermediate to high correlation between streamflow and hillslope hydrological dynamics (trench flow and overland flow), which highlights the importance of hillslope processes in this small watershed. Although the total contribution of SSF to total event catchment runoff is rather small, the contribution during peak flow is moderate to substantial. Additionally, there is process synchronicity between spatially discontiguous measurement points across scales, potentially indicating subsurface flowpath connectivity. Our findings stress the need for (i) a combination of observations at different spatial scales, and (ii) a consideration of the high spatial variability of SSF at the plot and hillslope scale when designing monitoring networks and assessing hydrological connectivity. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

16.
In tidal environments, the response of suspended sediment concentration (SSC) to the current velocity is not instantaneous, the SSC lagging behind the velocity (phase lag), and the amplitude of SSC variation decreasing with height above the bed (amplitude attenuation). In order to quantitatively describe this phenomenon, a one-dimensional vertical advection–diffusion equation of SSC is derived analytically for uniform unsteady tidal flow by defining a concentration boundary condition using a constant vertical eddy diffusivity and sediment settling velocity. The solution, in simple and straightforward terms, shows that the vertical phase lag increases linearly with the height above the bed, while the amplitude of the SSC variation decreases exponentially with the height. The relationship between the SSC and the normalized current velocity can be represented by an ellipse or a line, depending on the phase lag. The lag of sediment movement or “diffusion/settling lag” is the mechanism generating the phase lag effect. Field observations used for validation show that the theoretically predicted and the observed curves of the vertical SSC phase lag and amplitude attenuation show reasonable agreement. The procedure proposed in this paper substantially simplifies the modeling of suspended matter transport in tidal flows.  相似文献   

17.
In watershed modelling, the traditional practice of arbitrarily filling topographic depressions in digital elevation models has raised concerns. Advanced high‐resolution remote sensing techniques, including airborne scanning laser altimetry, can identify naturally occurring depressions that impact overland flow. In this study, we used an ensemble physical and statistical modelling approach, including a 2D hydraulic model and two‐point connectivity statistics, to quantify the effects of depressions on high‐resolution overland flow patterns across spatial scales and their temporal variations in single storm events. Computations for both models were implemented using graphic processing unit‐accelerated computing. The changes in connectivity statistics for overland flow patterns between airborne scanning laser altimetry‐derived digital elevation models with (original) and without (filled) depressions were used to represent the shifts of overland flow response to depressions. The results show that depressions can either decrease or increase (to a lesser degree and shorter duration) the probability that any two points (grid locations) are hydraulically connected by overland flow pathways. We used macro‐connectivity states (Φ) as a watershed‐specific indicator to describe the spatiotemporal thresholds of connectivity variability caused by depressions. Four states of Φ are identified in a studied watershed, and each state represents different magnitudes of connectivity and connectivity changes (caused by depressions). The magnitude of connectivity variability corresponds to the states of Φ, which depend on the topological relationship between depressions, the rising/recession limb, and the total rainfall amount in a storm event. In addition, spatial distributions of connectivity variability correlate with the density of depression locations and their physical structures, which cause changes in streamflow discharge magnitude. Therefore, this study suggests that depressions are “nontrivial” in watershed modelling, and their impacts on overland flow should not be neglected. Connectivity statistics at different spatial scales and time points within a watershed provide new insights for characterizing the distributed and accumulated effects of depressions on overland flow.  相似文献   

18.
Non-perennial streams comprise over half of the global stream network and impact downstream water quality. Although aridity is a primary driver of stream drying globally, surface flow permanence varies spatially and temporally within many headwater streams, suggesting that these complex drying patterns may be driven by topographic and subsurface factors. Indeed, these factors affect shallow groundwater flows in perennial systems, but there has been only limited characterisation of shallow groundwater residence times and groundwater contributions to intermittent streams. Here, we asked how groundwater residence times, shallow groundwater contributions to streamflow, and topography interact to control stream drying in headwater streams. We evaluated this overarching question in eight semi-arid headwater catchments based on surface flow observations during the low-flow period, coupled with tracer-based groundwater residence times. For one headwater catchment, we analysed stream drying during the seasonal flow recession and rewetting period using a sensor network that was interspersed between groundwater monitoring locations, and linked drying patterns to groundwater inputs and topography. We found a poor relationship between groundwater residence times and flowing network extent (R2 < 0.24). Although groundwater residence times indicated that old groundwater was present in all headwater streams, surface drying also occurred in each of them, suggesting old, deep flowpaths are insufficient to sustain surface flows. Indeed, the timing of stream drying at any given point typically coincided with a decrease in the contribution from near-surface sources and an increased relative contribution of groundwater to streamflow at that location, whereas the spatial pattern of drying within the stream network typically correlated with locations where groundwater inputs were most seasonally variable. Topographic metrics only explained ~30% of the variability in seasonal flow permanence, and surprisingly, we found no correlation with seasonal drying and down-valley subsurface storage area. Because we found complex spatial patterns, future studies should pair dense spatial observations of subsurface properties, such as hydraulic conductivity and transmissivity, to observations of seasonal flow permanence.  相似文献   

19.
Connectivity describes the efficiency of material transfer between geomorphic system components such as hillslopes and rivers or longitudinal segments within a river network. Representations of geomorphic systems as networks should recognize that the compartments, links, and nodes exhibit connectivity at differing scales. The historical underpinnings of connectivity in geomorphology involve management of geomorphic systems and observations linking surface processes to landform dynamics. Current work in geomorphic connectivity emphasizes hydrological, sediment, or landscape connectivity. Signatures of connectivity can be detected using diverse indicators that vary from contemporary processes to stratigraphic records or a spatial metric such as sediment yield that encompasses geomorphic processes operating over diverse time and space scales. One approach to measuring connectivity is to determine the fundamental temporal and spatial scales for the phenomenon of interest and to make measurements at a sufficiently large multiple of the fundamental scales to capture reliably a representative sample. Another approach seeks to characterize how connectivity varies with scale, by applying the same metric over a wide range of scales or using statistical measures that characterize the frequency distributions of connectivity across scales. Identifying and measuring connectivity is useful in basic and applied geomorphic research and we explore the implications of connectivity for river management. Common themes and ideas that merit further research include; increased understanding of the importance of capturing landscape heterogeneity and connectivity patterns; the potential to use graph and network theory metrics in analyzing connectivity; the need to understand which metrics best represent the physical system and its connectivity pathways, and to apply these metrics to the validation of numerical models; and the need to recognize the importance of low levels of connectivity in some situations. We emphasize the value in evaluating boundaries between components of geomorphic systems as transition zones and examining the fluxes across them to understand landscape functioning. © 2018 John Wiley & Sons, Ltd.  相似文献   

20.
GPS Height and Gravity Variations Due to Ocean Tidal Loading Around Taiwan   总被引:2,自引:1,他引:1  
This study presents predicts ocean tidal loading (OTL) effects using a Green’s function approach and validates a novel tidal model for Taiwan. Numerical integration of OTL is performed using the Gauss quadrature method and a local tidal model for the inner zone and a global model for the outer zone. Observed time series of GPS-derived vertical displacements and gravity variations (3–7 days) at five co-located GPS-gravimeter stations along the South East China and Taiwan coasts were utilized to assess the accuracy of the proposed models and two other models. The OTL-induced gravity variations are 3–16 μgal and the vertical site displacements are 13–27 mm. Generally, an OTL model using a mixed global and local tidal model generates better agreement with the observations than an OTL model using a global tidal model only. However, containing a local model inside a global model does not always produce a good agreement with the observations. The relatively large discrepancies between modeled and observed OTL values at some stations indicate that there is a need for an improved local tidal model in the study area. An erratum to this article can be found at  相似文献   

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

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