共查询到20条相似文献，搜索用时 20 毫秒
1.
Melina G. Sidiropoulou Konstadinos N. Moutsopoulos Vassilios A. Tsihrintzis 《水文研究》2007,21(4):534-554
This study focuses on the determination of the Forchheimer equation coefficients a and b for non‐Darcian flow in porous media. Original theoretical equations are evaluated and empirical relations are proposed based on an investigation of available data in the literature. The validity of these equations is checked using existing experimental data, and their accuracy versus existing approaches is studied. On the basis of this analysis, some insight into the physical background of the phenomenon is also provided. The dependence of the coefficients a and b on the Reynolds number is also detected, and potential future research areas, e.g. investigation of inertial effects for consolidated porous media, are pointed out. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
2.
Surface‐wetting properties are an important cause of changing the groundwater and two‐phase fluid flows. Various factors affecting the surface wettability were investigated in a parallel‐walled glass fracture with non‐aqueous phase liquid (NAPL) (gasoline, diesel, trichloroethylene, and creosote) wetted surfaces. First, the effect of the duration of NAPL exposure on wettability change was considered at pre‐wet fracture surfaces using the various NAPL species, and the result showed that the surface became hydrophobic after the exposure time of NAPL exceeded 2000 min. Second, the initial wetting state of the surface affected the timing when the wettability change begins as well as the extent of the wettability change in an NAPL‐wetted rock fractures. Under the dry condition, the wettability change was completed within a very short time of exposure to NAPL (~5 min), and then it finally reached the intermediate and weakly NAPL wetting (contact angle of 118°). Under the pre‐wet condition, a relatively long time of exposure (~5000 min) was needed to observe the obvious change of the surface wettability, which was changed up to strongly NAPL wetting (contact angle of 142°). Third, the wettability changed by NAPL exposure was stable and maintained for a long time, regardless of water flushing rate and temperature. Finally, the wettability change by the exposure of NAPL on parallel fracture surfaces was evaluated at various groundwater flow velocities. Result showed that groundwater flow velocity has an important impact upon measured contact angle. Although fracture surfaces were exposed to NAPL at the low groundwater flow velocity, the wettability was not changed from hydrophilic to hydrophobic when the contact time between NAPL and mineral surfaces was not sufficient owing to the pulse‐type movement of NAPL. This implies that the variation of exposure pattern due to groundwater flow on the wettability change can be an important factor affecting the wettability change of fracture surface and migration behaviour at natural fractured rock aquifers in case of NAPL spill. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
3.
This paper is focused on the lag time of flow confluence. Lag is described empirically by the precipitation factor in hydrological modelling, but the traditional way to establish the relationship between lag and intensity is not very satisfactory in arid and semi‐arid regions. A total 215 rainfall–runoff experiments were conducted to reveal the effects of net rain duration on lag and intensity. The results show that a correlation between lag and rain intensity exists under certain conditions. There is a critical value in net rain intensity (0·8–1 mm min?1) in the correlation curve of lag and net rain intensity in the given experimental watershed in the laboratory. Similarly, a critical value of net rain duration also exists. This value is the total confluence time. The features of lag time change significantly if intensity or duration is less or greater than the critical value. The study also explores the joint effects of net rain intensity and duration on lag. The formula established among lag, intensity and duration resulted in a better fit. Therefore, the two‐parameter (intensity and duration) empirical formula for lag is better than the traditional single‐parameter (intensity) method. This two‐parameter correlation formula can also be applied to a temporally and spatially uneven runoff processes. A typical field watershed is selected to test the results of the experiments. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
4.
Peter F. Germann 《水文研究》2018,32(9):1166-1172
Preferential flow, a term that includes macropore flow, non‐equilibrium flow, and finger flow, stands in well known conflict with Richards' ( 1931 ) capillary flow. Acoustic velocity experiments demonstrate that preferential flow moves independently from, faster than, and before capillary flow during gravity‐driven infiltration. Viscous flow in permeable media is briefly introduced to the point where Richards' ( 1931 ) particular treatment of viscosity turns out as the hydro‐mechanical bifurcation from general laminar flow. Preferential flow is expected during significant infiltration, however, spatio‐temporarily limited according to the viscous‐flow regime. Two ways of delineating capillary flow from viscous flow reveal minimum path widths of preferential flow in the range of about 20 μm. 相似文献
5.
Barret L. Kurylyk Dylan J. Irvine Sean K. Carey Martin A. Briggs Dale D. Werkema Mariah Bonham 《水文研究》2017,31(14):2648-2661
Groundwater flow advects heat, and thus, the deviation of subsurface temperatures from an expected conduction‐dominated regime can be analysed to estimate vertical water fluxes. A number of analytical approaches have been proposed for using heat as a groundwater tracer, and these have typically assumed a homogeneous medium. However, heterogeneous thermal properties are ubiquitous in subsurface environments, both at the scale of geologic strata and at finer scales in streambeds. Herein, we apply the analytical solution of Shan and Bodvarsson ( 2004 ), developed for estimating vertical water fluxes in layered systems, in 2 new environments distinct from previous vadose zone applications. The utility of the solution for studying groundwater‐surface water exchange is demonstrated using temperature data collected from an upwelling streambed with sediment layers, and a simple sensitivity analysis using these data indicates the solution is relatively robust. Also, a deeper temperature profile recorded in a borehole in South Australia is analysed to estimate deeper water fluxes. The analytical solution is able to match observed thermal gradients, including the change in slope at sediment interfaces. Results indicate that not accounting for layering can yield errors in the magnitude and even direction of the inferred Darcy fluxes. A simple automated spreadsheet tool (Flux‐LM) is presented to allow users to input temperature and layer data and solve the inverse problem to estimate groundwater flux rates from shallow (e.g., <1 m) or deep (e.g., up to 100 m) profiles. The solution is not transient, and thus, it should be cautiously applied where diel signals propagate or in deeper zones where multi‐decadal surface signals have disturbed subsurface thermal regimes. 相似文献
6.
New methods for obtaining and quantifying spatially distributed subsurface moisture are a high research priority in process hydrology. We use simple linear regression analyses to compare terrain electrical conductivity measurements (EC) derived from multiple electromagnetic induction (EMI) frequencies to a distributed grid of water‐table depth and soil‐moisture measurements in a highly instrumented 50 by 50 m hillslope in Putnam County, New York. Two null hypotheses were tested: H0(1), there is no relationship between water table depth and EC; H0(2), there is no relationship between soil moisture levels and EC. We reject both these hypotheses. Regression analysis indicates that EC measurements from the low frequency EM31 meter with a vertical dipole orientation could explain over 80% of the variation in water‐table depth across the test hillslope. Despite zeroing and sensitivity problems encountered with the high frequency EM38, EC measurements could explain over 70% of the gravimetrically determined soil‐moisture variance. The use of simple moisture retrieval algorithms, which combined EC measurements from the EM31 and EM38 meters in both their vertical and horizontal orientations, helped increase the r2 coefficients slightly. This first hillslope hydrological analysis of EMI technology in this way suggests that it may be a promising method for the collection of a large number of distributed soilwater and groundwater depth measurements with a reasonable degree of accuracy. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
7.
A tuned mass damper (TMD) system consists of an added mass with properly functioning spring and damping elements for providing frequency‐dependent damping in a primary structure. The advantage of a friction‐type TMD, that is, a nonlinear TMD, is its energy dissipation via a friction mechanism. In contrast, the disadvantages of a passive friction TMD (PF‐TMD) are its fixed and predetermined slip load and loss of tuning and energy dissipation capabilities when it is in a stick state. A semi‐active friction TMD (SAF‐TMD) is used to overcome these disadvantages. The SAF‐TMD can adjust its slip force in response to structure motion. To verify its feasibility, a prototype SAF‐TMD was fabricated and tested dynamically using a shaking table test. A nonsticking friction control law was used to keep the SAF‐TMD activated and in a slip state in earthquakes at varying intensities. The shaking table test results demonstrated that: (i) the experimental results are consistent with the theoretical results; (ii) the SAF‐TMD is more effective than the PF‐TMD given a similar peak TMD stroke; and (iii) the SAF‐TMD can also prevent a residual TMD stroke in a PF‐TMD system. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
8.
The relationship between the longitudinal dispersion (DL) and Peclet number (Pe) is crucial for predicting and simulating tracer through the variable‐aperture fracture. In this study, the roughness of the self‐affine fracture wall was decomposed into primary roughness (relatively large‐scale waviness) and secondary roughness (relatively small‐scale waviness) by a multiscaled wavelet analysis technique. Based on the complete dispersion mechanism (diffusion, macrodispersion, and Taylor dispersion) in the variable‐aperture fracture, three relationships (second‐order, power‐law, and linear relationships) between the DL and Pe were investigated at large and small scales, respectively. Our results showed that the primary roughness mostly controlled the Taylor dispersion mechanism, whereas the secondary roughness was a dominant factor for the macrodispersion mechanism. Increasing the Hurst exponent and removing the secondary roughness led to the decreasing range of Pe where macrodispersion mechanism dominated the solute transport. It was found that estimating the DL from the power‐law relationship based on Taylor dispersion theory resulted in considerable errors, even in the range of Pe where the Taylor dispersion mechanism dominated. The exponent of the power‐law relationship increased as the secondary roughness was removed. Analysing the linear relationship between the DL and Pe revealed that the longitudinal dispersivity αL increased linearly. However, this linear increase became weak as the Taylor dispersion mechanism dominated. In the range of Pe where the macrodispersion mechanism dominated, increasing the Hurst exponent caused the increase of αL and the secondary roughness played a significant role in enhancing the αL. As the Taylor dispersion mechanism dominated, the αL was insensitive to the influence of multiscale roughness in variable‐aperture fractures. 相似文献
9.
This work attempted to locate clean and safe groundwater for irrigation use in the Choushui River alluvial fan. Multiple‐variable indicator kriging (MVIK) was adopted to evaluate numerous hydrochemical parameters for a standard of water quality for irrigation in Taiwan. Many hydrochemical parameters in groundwater were distinguished into three main categories—salinity/sodium hazard, nitrogen hazard and heavy metal hazard. Safe and potential hazardous regions of groundwater for irrigation were delineated according to different probabilities estimated by MVIK. The probabilistic results of the classifications gave an opportunity to explore the spatial uncertainty of the hazards and helped government administrators establish a sound policy associated with the development and management of groundwater resources. Analysis of the results indicate that the central distal‐fan and mid‐fan aquifers are the best places to extract clean and safe groundwater for irrigation, and the deep aquifer (exceeding 200 m depth) has wider regions with clean and safe groundwater for irrigation than shallow aquifers. The northern and southern aquifers, with multiple hazards, limit groundwater use for irrigation. Although the proximal‐fan aquifer is a zone of groundwater recharge, the high nitrogen content seriously affects the environment and is not suitable for irrigation use. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
10.
Based on the theory of gravity‐driven groundwater flow systems, we have developed a complex Flow System Sand‐Box Model (FSM). It enables the visual observations of the development and characteristics and temporal evolution of complex Tóthian flow systems in the laboratory. The configuration of the regional, intermediate and local flow systems can be controlled and observed; hydraulic head, flow direction and travel time can be measured; and the scale and shape of the sub‐flow systems as well as the path lines and flow lines can be observed directly. The experiments demonstrate the Tóthian flow systems in a small basin with multiple sources and sinks. Greater local topographic (water table) undulation will lead to larger local flow systems. Greater regional and less local topographic undulation will enhance the development of intermediate and regional flow systems. In homogeneous media, increasing fluid‐potential differences between source and sink increase the spatial scale of the generated flow systems. The FSM is a useful teaching aid and experimental device to study and develop an intuitive insight into gravity‐driven groundwater flow systems. It helps to visualize and understand the hydraulic properties and controlling factors of Tóthian flow systems and may be used to study problems related to the chemical and temperature characteristics of the flow systems as well. Copyright © 2010 John Wiley & Sons, Ltd. 相似文献
11.
Florian Kobierska Tobias Jonas Nena Griessinger Christian Hauck Stephan Huxol Stefano M. Bernasconi 《水文研究》2015,29(6):817-827
We implemented multiple independent field techniques to determine the direction and velocity of groundwater flow at a specific stream reach in a glacier forefield. Time‐lapse experiments were conducted using two electrical resistivity tomography (ERT) lines installed in a cross pattern. A circular array of groundwater tubes was also installed to monitor groundwater flow via discrete salt injections. Both inter‐borehole and ERT results confirmed this stream section as a losing reach and enabled quantification of the flow direction. Both techniques yielded advection velocities varying between 5.7 and 21.8 m/day. Estimates of groundwater flow direction and velocity indicated that groundwater infiltrates from the stream nearby and not from the adjacent lateral moraine. Groundwater age estimated from radon concentration measurements supported this hypothesis. Despite uncertainties inherent to each of the methods deployed, the combination of multiple field techniques allowed drawing consistent conclusions about local groundwater flow. We thus regard our multi‐method approach as a reliable way to characterize the two‐dimensional groundwater flow at sites where more invasive groundwater investigation techniques are difficult to carry out and local heterogeneities can make single measurements unreliable. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
12.
Surface water and groundwater in the Heihe river basin of China are interconnected and the pattern of water resources exploitation has a direct effect on the interaction of groundwater and surface water, especially on a downstream oasis. A three‐dimensional groundwater flow simulation model with eight model layers was established to simulate the regional groundwater flow in the multilayered aquifer system and the interaction among the rivers, springs, and groundwater. The model was calibrated not only with historical water levels but also with the investigated baseflow and spring flux. The simulation results of the numerical model match reasonably well with the observed groundwater levels, baseflow to rivers, and spring flux. The numerical simulation also demonstrates that the hydraulic connection between the river and the aquifers has transferred from the coupling to decoupling at some reaches. It is suggested that there is a vital need to reduce groundwater withdrawal and to rationalize the use of both groundwater and surface water in order to maintain sustainable development in the study area. Copyright © 2007 John Wiley & Sons, Ltd. 相似文献
13.
A two‐dimensional variable‐density groundwater flow and transport model was developed to provide a conceptual understanding of past and future conditions of nitrate (NO3) transport and estimate groundwater nitrate flux to the Gulf of Mexico. Simulation results show that contaminant discharge to the coast decreases as the extent of saltwater intrusion increases. Other natural and/or artificial surface waters such as navigation channels may serve as major sinks for contaminant loading and act to alter expected transport pathways discharging contaminants to other areas. Concentrations of NO3 in the saturated zone were estimated to range between 30 and 160 mg?L?1 as NO3. Relatively high hydraulic vertical gradients and mixing likely play a significant role in the transport processes, enhancing dilution and contaminant migration to depth. Residence times of NO3 in the deeper aquifers vary from 100 (locally) to about 300 years through the investigated aquifer system. NO3 mass fluxes from the shallow aquifers (0 to 5.7 × 104 mg?m?2?day?1) were primarily directed towards the navigation channel, which intersects and captures a portion of the shallow groundwater flow/discharge. Direct NO3 discharge to the sea (i.e. Gulf of Mexico) from the shallow aquifer was very low (0 to 9.0 × 101 mg · m?2?day?1) compared with discharge from the deeper aquifer system (0 to 8.2 × 103 mg?m?2?day?1). Both model‐calibrated and radiocarbon tracer‐determined contaminant flux estimates reveal similar discharge trends, validating the use of the model for density‐dependent flow conditions. The modelling approach shows promise to evaluate contaminant and nutrient loading for similar coastal regions worldwide. Copyright © 2015 John Wiley & Sons, Ltd. 相似文献
14.
David L. Dunkerley 《地球表面变化过程与地形》2003,28(5):475-491
Friction coefficients in overland flows are customarily estimated from mean flow properties (depth, velocity, slope) that subsume spatial variations in flow arising from two major causes: microtopography and obstacles. This paper uses laboratory experiments in shallow flumes to examine the extent of non‐uniformity in flow conditions associated with each cause. Randomly placed emergent obstacles in a flume with a shallow axial channel generally yielded higher hydraulic roughness than the same pattern of obstacles on a planar flume, as well as greater variation in roughness as the obstacle locations were altered. In both flumes, hydraulic roughness fell with increasing Reynolds number for 10% obstacle cover, showed a flattening trend at 20% cover, and exhibited a convex‐downward trend at 30% obstacle cover. These results indicate the progressive onset of flow controls at narrow gaps in the obstacle field. In such flows, the use of mean flow properties conceals the existence of two main subdivisions of flow: flow filaments and backwater flows. In the experiments, flow filaments involved velocities more than twice the overall mean, whereas backwater flows were much slower than the mean. The existence of fast‐moving flow filaments may be significant in understanding soil transport in surface runoff, and backwater depths may modify splash detachment. Similarly, friction coefficients that fail to reflect these important non‐uniform flow components may not be optimal for hydraulic calculations or in erosion models. It is concluded that new approaches to observing and processing flow data may be required, in order to avoid the loss of important flow detail that is entailed in assuming uniform flow conditions. Copyright © 2003 John Wiley & Sons, Ltd. 相似文献
15.
This study uses a numerical model to investigate the groundwater flow and salt transport mechanisms below islands in the Okavango Delta. Continuous evapotranspiration on the islands results in accumulation of solutes and the formation of a saline boundary layer, which may eventually become unstable. A novel Lagrangian method is employed in this study and compared to other numerical methods. The numerical results support the geophysical observations of density fingering on Thata Island. However, the process is slow and it takes some hundreds of years until density fingering is triggered. The results are sensitive to changes of the hydraulic gradient and the evapotranspiration rate. Small changes may lead to different plume developments. Results further demonstrate that density effects may be entirely overridden by lateral flow on islands embedded in a sufficiently high regional hydraulic gradient. 相似文献
16.
Stable isotope variations are extremely useful for flow partitioning within the hydrologic cycle but remain poorly understood throughout the tropics, particularly in watersheds with rapidly infiltrating soils, such as Andisols in Central America. This study examines the fluctuations of stable isotope ratios (δ18O and δ2H) in the hydrologic components of a tropical coffee agroforestry watershed (~1 km2) with Andisol soils in Costa Rica. Samples were collected in precipitation, groundwater, springs, and stream water over 2 years. The local meteoric water line for the study site was δ2H = 8.5 δ18O + 18.02 (r2 = 0.97, n = 198). The isotope ratios in precipitation exhibited an enriched trend during the dry season and a notable depletion at the beginning of the wet season. The δ18O compositions in groundwater (average = ?6.4‰, σ = 0.7) and stream water (average = ?6.7‰, σ = 0.6) were relatively stable over time, and both components exhibited more enriched values in 2013, which was the drier year. No strong correlation was observed between the isotope ratios and the precipitation amount at the event or daily time‐step, but a correlation was observed on a monthly scale. Stream water and base flow hydrograph separations based on isotope end‐member estimations showed that pre‐event water originating from base flow was prevalent. However, isotope data indicate that event water originating from springs appears to have been the primary driver of initial rises in stream flow and peak flows. These results indicate that isotope sampling improves the understanding of water balance components, even in a tropical humid location, where significant variations in rainfall challenge current modelling efforts. Further research using fine‐scale hydrometric and isotopic data would enhance understanding the processes driving spring flow generation in watersheds. 相似文献
17.
18.
Anisotropy and heterogeneity of hydraulic conductivity (K) are seldom considered in models of mire hydrology. We investigated the effect of anisotropy and heterogeneity on groundwater flow in bog peat using a steady‐state groundwater model. In five model simulations, four sets of K data were used. The first set comprised measured K values from an anisotropic and heterogeneous bog peat. These data were aggregated to produce the following simplified data sets: an isotropic and heterogeneous distribution of K; an isotropic and homogeneous distribution; and an anisotropic and homogeneous distribution. We demonstrate that, where anisotropy and heterogeneity exist, groundwater flow in bog peat is complex. Fine‐scale variations in K have the potential to influence patterns and rates of groundwater flow. However, for our data at least, it is heterogeneity and not anisotropy that has the greater influence on producing complex patterns of groundwater flow. We also demonstrate that patterns and rates of groundwater flow are simplified and reduced when measured K values are aggregated to create a more uniform distribution of K. For example, when measured K values are aggregated to produce isotropy and homogeneity, the rate of modelled seepage is reduced by 28%. We also show that when measured K values are used, the presence of a drainage ditch can increase seepage through a modelled cross‐section. Our work has implications for the accurate interpretation of hydraulic head data obtained from peat soils, and also the understanding of the effect of drainage ditches on patterns and rates of groundwater flow. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
19.
Hossein Afzalimehr 《水文研究》2010,24(3):331-341
Non‐uniform flows encompassing both accelerating and decelerating flows over a cobble‐bed flume have been experimentally investigated in a flume at a scale of intermediate relative submergence. Measurements of mean longitudinal flow velocity u, and determinations of turbulence intensities u′, v′, w′, and Reynolds shear stress ?ufwf have been made. The longitudinal velocity distribution was divided into the inner zone close to the bed and the outer zone far from the bed. In the inner zone of the boundary layer (near the bed) the velocity profile closely followed the ‘Log Law’; however, in the outer zone the velocity distribution deviated from the Log Law consistently for both accelerating and decelerating flows and the changes in bed slopes ranging from ?2% to + 2% had no considerable effect on the outer zone. For a constant bed slope (S = ±0·015), the larger the flow rate, the smaller the turbulence intensities. However, no detectable pattern has been observed for u′, v′ and w′ distributions near the bed. Likewise, for a constant flow rate (Q = 0·040 m3/s), with variation in bed slope the longitudinal turbulent intensity profile in the longitudinal direction remained concave for both accelerating and decelerating flows; whereas vertical turbulent intensity (w′) profile presented no specific form. The results reveal that the positions of maximum values of turbulence intensities and the Reynolds shear stress depend not only on the flow structure (accelerating or decelerating) but also on the intermediate relative submergence scale. Copyright © 2009 John Wiley & Sons, Ltd. 相似文献
20.
Compared to hydrograph recession analysis, which is widely applied in engineering hydrology, the quantitative assessment of stream salinity with time (i.e. the salinograph) has received significantly less attention. In particular, while in many previous hydrological studies an inverse relationship between hydrograph and salinograph responses is apparent, the concept of salinity accession (the inversely related salinity counterpart to hydrograph recession) has not been introduced nor quantitatively evaluated in previous literature. In this study, we conduct a mathematical analysis of salinograph accession, and determine new quantitative relationships between salinity accession and hydrograph recession parameters. An equation is formulated that reproduces the general trend in salinity accession. A salinity accession parameter kc is then introduced and is shown to be the ratio of direct runoff to total stream flow recession parameters: kr/k. The groundwater recession parameter kg was estimated using a simple and rapid method that uses both salinograph and hydrograph data. Salinity accession type‐curves illustrate that under certain conditions, the relative steepness of individual salinographs is dependent upon the ratio of groundwater salinity to direct runoff salinity: Cg/Cr. The salinity accession algorithms are applied to two contrasting field settings: Scott Creek, South Australia and Sandy Creek, northern Queensland, Australia. It was found that kg > k during periods of obvious stream flow recession, for the events analysed. Salinograph accession behaviour was fairly similar for both sites, despite contrasting environments. Using assumed end‐member salinities for groundwater and direct runoff based upon field observations, the behaviour of kc from the Scott Creek site was approximately reproduced by varying the initial groundwater to runoff flow ratio: Qg0/Qr0, within reasonable parameter ranges. The use of salinograph information when used in addition to standard hydrograph analyses provided useful information on recession characteristics of stream components. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献