首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献，搜索用时 453 毫秒
1.
A three‐dimensional, numerical free‐surface‐flow simulation model is developed to investigate hydrodynamics of a lake and the predictive capabilities of the numerical model are validated by comparing them with field results obtained from Lake Naini, Uttar Pradesh, using environmental isotope (δ18O, δD and 3H) techniques. This has been carried out in order to understand the hydrodynamics of a lake (i.e. circulations, mixing and transport of water within the lake). Copyright © 2002 John Wiley & Sons, Ltd.  相似文献   

2.
Stream flow predictions in ungauged basins are one of the most challenging tasks in surface water hydrology because of nonavailability of data and system heterogeneity. This study proposes a method to quantify stream flow predictive uncertainty of distributed hydrologic models for ungauged basins. The method is based on the concepts of deriving probability distribution of model's sensitive parameters by using measured data from a gauged basin and transferring the distribution to hydrologically similar ungauged basins for stream flow predictions. A Monte Carlo simulation of the hydrologic model using sampled parameter sets with assumed probability distribution is conducted. The posterior probability distributions of the sensitive parameters are then computed using a Bayesian approach. In addition, preselected threshold values of likelihood measure of simulations are employed for sizing the parameter range, which helps reduce the predictive uncertainty. The proposed method is illustrated through two case studies using two hydrologically independent sub‐basins in the Cedar Creek watershed located in Texas, USA, using the Soil and Water Assessment Tool (SWAT) model. The probability distribution of the SWAT parameters is derived from the data from one of the sub‐basins and is applied for simulation in the other sub‐basin considered as pseudo‐ungauged. In order to assess the robustness of the method, the numerical exercise is repeated by reversing the gauged and pseudo‐ungauged basins. The results are subsequently compared with the measured stream flow from the sub‐basins. It is observed that the measured stream flow in the pseudo‐ungauged basin lies well within the estimated confidence band of predicted stream flow. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

3.
Although it has been increasingly acknowledged that groundwater flow pattern is complicated in the three‐dimensional (3‐D) domain, two‐dimensional (2‐D) water table‐induced flow models are still widely used to delineate basin‐scale groundwater circulation. However, the validity of 2‐D cross‐sectional flow field induced by water table has been seldom examined. Here, we derive the analytical solution of 3‐D water table‐induced hydraulic head in a Tóthian basin and then examine the validity of 2‐D cross‐sectional models by comparing the flow fields of selected cross sections calculated by the 2‐D cross‐sectional model with those by the 3‐D model, which represents the “true” cases. For cross sections in the recharge or discharge area of the 3‐D basin, even if head difference is not significant, the 2‐D cross‐sectional models result in flow patterns absolutely different from the true ones. For the cross section following the principal direction of groundwater flow, although 2‐D cross‐sectional models would overestimate the penetrating depth of local flow systems and underestimate the recharge/discharge flux, the flow pattern from the cross‐sectional model is similar to the true one and could be close enough to the true one by adjusting the decay exponent and anisotropy ratio of permeability. Consequently, to determine whether a 2‐D cross‐sectional model is applicable, a comparison of hydraulic head difference between 2‐D and 3‐D solutions is not enough. Instead, the similarity of flow pattern should be considered to determine whether a cross‐sectional model is applicable. This study improves understanding of groundwater flow induced by more natural water table undulations in the 3‐D domain and the limitations of 2‐D models accounting for cross‐sectional water table undulation only.  相似文献   

4.
Vic Kelson 《Ground water》2012,50(6):918-926
Groundwater flow models are commonly used to design new wells and wellfields. As the spatial scale of the problem is large and much local‐scale detail is not needed, modelers often utilize two‐dimensional (2D) or quasi three‐dimensional models based on the Dupuit‐Forchheimer assumption. Dupuit models offer a robust set of tools for simulating regional groundwater flow including interactions with surface waters, the potential for well interference, and varying aquifer properties and recharge rates. However, given an assumed operating water level or drawdown at a well screen, Dupuit models systematically overpredict well yields. For design purposes, this discrepancy is unacceptable, and a method for predicting accurate well yields is needed. While published methods exist for vertical wells, little guidance is available for predicting yields in horizontal screens or collector wells. In plan view, a horizontal screen has a linear geometry, and will likely extend over several neighboring cells that may not align with rows or columns in a numerical model. Furthermore, the model must account for the effects of converging three‐dimensional (3D) flow to the well screens and hydraulic interference among the well screens; these all depend on the design of a specific well. This paper presents a new method for simulating the yield of angled or horizontal well screens in numerical groundwater flow models, specifically using the USGS code MODFLOW. The new method is compared to a detailed, 3D analytic element model of a collector well in a field of uniform flow.  相似文献   

5.
The lake levels in Lake Michigan‐Huron have recently fallen to near historical lows, as has the elevation difference between Lake Michigan‐Huron compared to Lake Erie. This decline in lake levels has the potential to cause detrimental impacts on the lake ecosystems, together with social and economic impacts on communities in the entire Great Lakes region. Results from past work suggest that morphological changes in the St Clair River, which is the only natural outlet for Lake Michigan‐Huron, could be an appreciable factor in the recent trends of lake level decline. A key research question is whether bed erosion within the river has caused an increase in water conveyance, therefore, contributed to the falling lake level. In this paper, a numerical modeling approach with field data is used to investigate the possibility of sediment movement in the St Clair River and assess the likelihood of morphological change under the current flow regime. A two‐dimensional numerical model was used to study flow structure, bed shear stress, and sediment mobility/armoring over a range of flow discharges. Boundary conditions for the numerical model were provided by detailed field measurements that included high‐resolution bathymetry and three‐dimensional flow velocities. The results indicate that, without considering other effects, under the current range of flow conditions, the shear stresses produced by the river flow are too low to transport most of the coarse bed sediment within the reach and are too low to cause substantial bed erosion or bed scour. However, the detailed maps of the bed show mobile bedforms in the upper St Clair River that are indicative of sediment transport. Relatively high shear stresses near a constriction at the upstream end of the river and at channel bends could cause local scour and deposition. Ship‐induced propeller wake erosion also is a likely cause of sediment movement in the entire reach. Other factors that may promote sediment movement, such as ice cover and dredging in the lower river, require further investigation. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

6.
Complex flow processes at river bifurcations and the influence of the layout of a bifurcation make it difficult to predict sediment distribution over the downstream branches in case bedload transport dominates. In one‐dimensional models we need a nodal point relationship that prescribes the distribution of sediment over the downstream branches. We have identified which factors need to be included in such a relationship for the division of bedload transport at bifurcations. Next, irrotational flow theory for idealized geometries has been used to derive a simple physics‐based nodal point relationship that accounts for the effects of helical flow in the situation that a channel takes off under an angle from a straight main channel. This first step towards a complete nodal point relationship is applicable to bedload transport situations if the flow is clearly curved and if there is no pronounced bed topography. The relationship has been tested against data from a unique set of laboratory measurements, numerical data and data from a scale model of the Rhine bifurcation at Pannerden in the Netherlands. We find that the derived model yields a reasonable prediction of the sediment division over the downstream branches, and yields better predictions than the Wang et al. model for the situation considered. Considering the relative complexity and limited accuracy of the nodal point relationship for the effect of helical flow alone, however, we conclude thatderiving a practical physics‐based 1‐D relationship including all relevant processes is not feasible. We therefore recommend 2‐D or 3‐D modelling for all cases in general where morphological evolution depends on the division of bedload transport at bifurcations. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

7.
A riparian ecosystem downstream of a small dam in central Texas was instrumented for sap flow, soil moisture content, and stream level from 2001. Stable isotopes in water (D and 18O) were analysed from rainfall, stream, lake, and cored sapwood cellulose from cedar elm (Ulmus crassifolia). The isotope signature of water source to cedar elm was identified by back calculation starting with the water isotopes in cellulose, and accounting for leaf‐water evaporation and biological fractionation during cellulose synthesis. The estimated mean isotope of the source water to cedar elm was enriched above rainfall in similarity to stream water during 2002. Flow paths that may have contributed to estimated variability from regional base flow and recharge water were identified using the variably saturated HYDRUS‐2D model. Copyright © 2004 John Wiley & Sons, Ltd.  相似文献   

8.
The salinity of Lake Kinneret, Israel, is significantly higher than the salinity of the water from surface streams that flow to the lake. The relatively high salinity is a result of the activity of saline springs located at the bottom of the lake.The purpose of this work is to establish a general model for the salinization mechanism of Lake Kinneret. The model is based on the main components of the annual water and solute balance. Changes in time of the solute mass of the lake were described as a differential equation of a linear reservoir on an annual time scale. The model assumes that under any long-term operation policy of the lake, the components of the annual solute and water balance stay nearly constant in time.The model was tested for both steady-state conditions, and during changes in time, against measured lake salinity over the years 1968-2000. It was found that the major changes of lake salinity throughout the years were described well, despite the variety of rainfall amounts. Predictions of the expected lake salinity changes were proposed for the cases of controlled increase or decrease of saline springs discharge to the lake; for the changes of water quantity allowed to flow into or pumped out of the lake; and for various initial salinities. Predictions agree well with previous predictions made by statistical models.  相似文献   

9.
Two‐dimensional (2D) hydrodynamic models have been increasingly used to quantify aquatic habitat and stream processes, such as sediment transport, streambed morphological evolution, and inundation extents. Because river topography has a strong influence on predicted hydraulic conditions, 2D models require accurate and detailed bathymetric data of the stream channel and surrounding floodplains. Besides collection of mass points to construct high‐resolution three‐dimensional surfaces, bathymetries may be interpolated from cross‐sections. However, limited information is available on the effects of cross‐section spacing and the derived interpolated bathymetry on 2D model results in large river systems. Here, we investigated the effects of cross‐section spacing on flow properties simulated with 2D modeling at low, medium and high discharges in two morphologically different reaches, a simple (almost featureless with low sinuosity) and a complex (presenting pools, riffles, runs, contractions and expansions) reach of the Snake River (Idaho, USA), the tenth largest river in the United States in terms of drainage area. We compared the results from 2D models developed with complete channel bathymetry acquired with multibeam sonar data and photogrammetry, with 2D model results that were developed using interpolated topography from uniformly distributed transects. Results indicate that cross‐sections spaced equal to or greater than 2 times the average channel width (W*) smooths the bathymetry and suppresses flow structures. Conversely, models generated with cross‐sections spaced at 0.5 and 1 W* have stream flow properties, sediment mobility and spatial habitat distribution similar to those of the complete bathymetry. Furthermore, differences in flow properties between interpolated and complete topography models generally increase with discharge and with channel complexity. Copyright © 2013 John Wiley & Sons, Ltd.  相似文献   

10.
The air flow and heat transfer in gravel embankment in permafrost areas   总被引：3，自引：0，他引：3  
A comparative numerical investigation of transient temperature profile and pore-air velocities in horizontal rock block embankments are conducted using the "gravels model", in which the embankment is composed of stones and air, and the "porous media model" respectively. As the velocities from the "gravels model" directly reflect the true flow of air and winter-time convection, in this paper it can be concluded that computational results from the "gravels model"are superior to the "porous media model". In addition, the "gravels model" has the advantages of reflecting the effect of the dimensions and collocation of gravels upon the temperature fields.Therefore, the computation of the gravels embankment is mainly based on the gravels model.Simulation results show that in summer, a clockwise circulation of the pore-air extends throughout most of the embankment. However its motion is very weak that results in relatively straight horizontal isotherm lines. And heat transfer is mainly maintained through conduction. But in winter, the pore-air velocities are higher and multiple vortexes are formed in the embankment.Natural convection then becomes the dominant influence on the isotherm shapes within the embankment. The isotherms are complex and alternative upward and downward flowing plumes exist. The winter-time convection can further reduce the temperature of the foundation soil beneath the gravel embankment. In addition, the effects of the gravel dimensions within the embankment have been analyzed and compared in the gravels model. It shows that in winter, large stones, e.g. 200 mm, lead to stronger vortexes than those of small stones, say 60 mm. Consequently, the zone of low-temperature beneath the large-stone embankment extends deeper into the ground.  相似文献   

11.
Growing human pressure and potential change in precipitation pattern induced by climate change require a more efficient and sustainable use of water resources. Hydrological models can provide a fundamental contribution to this purpose, especially as increasing availability of meteorological data and forecast allows for more accurate runoff predictions. In this article, two models are presented for describing the flow formation process in a sub‐alpine catchment: a distributed parameter, physically based model, and a lumped parameter, empirical model. The scope is to compare the two modelling approaches and to assess the impact of hydrometeorological information, either observations or forecast, on water resources management. This is carried out by simulating the real‐time management of the regulated lake that drains the catchment, using the inflow predictions provided by the two models. Copyright © 2010 John Wiley & Sons, Ltd.  相似文献   

12.
This is the first of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory (also known as Dupuit or Boussinesq theory). Here, we examine the effect of lateral flow on the downward fluxes of water and solutes through perched groundwater at steady state. We derive analytical expressions describing the decline in the downward flux rate with depth. Using these, we obtain analytical expressions for water age in a number of cases. The results show that when the permeability field is homogeneous, the spatial structure of water age depends qualitatively on a single dimensionless number, Hi. This number captures the relative contributions to the lateral hydraulic potential gradient of the relief of the lower‐most impermeable boundary (which may be below the weathering front within permeable or incipiently weathered bedrock) and the water table. A “scaled lateral symmetry” exists when Hi is low: age varies primarily in the vertical dimension, and variations in the horizontal dimension x almost disappear when the vertical dimension z is expressed as a fraction z/H(x) of the laterally flowing system thickness H(x). Taking advantage of this symmetry, we show how the lateral dimension of the advection–diffusion‐reaction equation can be collapsed, yielding a 1‐D vertical equation in which the advective flux downward declines with depth. The equation holds even when the permeability field is not homogeneous, as long as the variations in permeability have the same scaled lateral symmetry structure. This new 1‐D approximation is used in the accompanying paper to extend chemical weathering models derived for 1‐D columns to hillslope domains.  相似文献   

13.
Soil erosion by water is the root cause of ecological degradation in the Shiwalik foothills of Northern India. Simulation of runoff and its component processes is a pre‐requisite to develop the management strategies to tackle the problem, successfully. A two‐dimensional physically based distributed numerical model, ROMO2D has been developed to simulate runoff from small agricultural watersheds on an event basis. The model employs the 2‐D Richards equation with sink term to simulate infiltration and soil moisture dynamics in the vadoze zone under variable rainfall conditions, and 2‐D Saint‐Venant equations under the kinematic wave approximation along with Manning's equation as the stage‐discharge equation for runoff routing. The various flow‐governing equations have been solved numerically by employing a Galerkin finite element method for spatial discretization using quadrilateral elements and finite difference techniques for temporal solutions. The ROMO2D computer program has been developed as a class‐based program, coded in C + + in such a way that with minor modifications, the model can be used to simulate runoff on a continuous basis. The model writes output for a runoff hydrograph of each storm. Model development is described in this paper and the results of model testing and field application are to be presented in a subsequent paper. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

14.
This study aims to realistically simulate the seismic responses of typical highway bridges in California with considerations of soil–structure interaction effects. The p‐y modeling approaches are developed and validated for embankments and pile foundations of bridges. The p‐y approach models the lateral and vertical foundation flexibility with distributed p‐y springs and associated t‐z and q‐z springs. Building upon the existing p‐y models for pile foundations, the study develops the nonlinear p‐y springs for embankments based on nonlinear 2D and 3D continuum finite element analysis under passive loading condition along both longitudinal and transverse directions. Closed‐form expressions are developed for two key parameters, the ultimate resistant force pult and the displacement y50, where 0.5pult is reached, of embankment p‐y models as functions of abutment geometry (wall width and height, embankment fill height, etc.) and soil material properties (wall‐soil friction angle, soil friction angle, and cohesion). In order to account for the kinematic and site responses, depth‐varying ground motions are derived and applied at the free‐end of p‐y springs, which reflects the amplified embankment crest motion. The modeling approach is applied to simulate the seismic responses of the Painter Street Bridge and validated through comparisons with the recorded responses during the 1992 Petrolia earthquake. It is demonstrated that the flexibility and motion amplification at end abutments are the most crucial modeling aspects. The developed p‐y models and the modeling approach can effectively predict the seismic responses of highway bridges. Copyright © 2016 John Wiley & Sons, Ltd.  相似文献   

15.
The knowledge of the velocity model in microseismic jobs is critical to achieving statistically reliable microseismic event locations. The design of microseismic networks and the limited sources for calibration do not allow for a full tomographic inversion. We propose optimizing a priori velocity models using a few active shots and a non‐linear inversion, suitable to poorly constrained systems. The considered models can be described by several layers with different P‐ and S‐wave velocities. The velocities may be constant or have 3D gradients; the layer interfaces may be simple dipping planes or more complex 3D surfaces. In this process the P‐ and S‐ wave arrival times and polarizations measured on the seismograms constitute the observed data set. They are used to estimate two misfit functions: i) one based on the measurement residuals and ii) one based on the inaccuracy of the source relocation. These two functions are minimized thanks to a simulated annealing scheme, which decreases the risk of converging to a local solution within the velocity model. The case study used to illustrate this methodology highlights the ability of this technique to constrain a velocity model with dipping layers. This was performed by jointly using sixteen perforation shots recorded during a multi‐stage fracturing operation from a single string of 3C‐receivers. This decreased the location inaccuracies and the residuals by a factor of six. In addition, the retrieved layer dip was consistent with the pseudo‐horizontal trajectories of the wells and the background information provided by the customer. Finally, the theoretical position of each calibration shot was contained in the uncertainty domain of the relocation of each shot. In contrast, single‐stage inversions provided different velocity models that were neither consistent between each other nor with the well trajectories. This example showed that it is essential to perform a multi‐stage inversion to derive a better updated velocity model.  相似文献   

16.
Five tracer experiments have been performed in a coarse‐textured soil near the new main airport at Gardermoen, Norway. In two lysimeter walls, 30 and 40 measuring points form the basis for spatial moment calculations. Although experiments were performed under different meteorological conditions (autumn and snowmelt) and at two different sites, the ratios of centres of vertical mass over cumulative infiltration were of the same order of magnitude, indicating a gravity‐dominated flow. Two‐dimensional transport simulations with SUTRA (Voss, 1984), with a priori estimated input parameters and random fields of soil hydraulic properties revealed a relatively good agreement with the experimental results. Three possible sources of heterogeneity affecting the vertical displacement of solute during snowmelt were identified: variability of soil physical properties, soil surface elevations and variability of ground frost during the melting period. To obtain accurate predictions, soil heterogeneity was the most important factor to characterize for the coarse‐textured soil under consideration. Copyright © 1999 John Wiley & Sons, Ltd.  相似文献   

17.
Long‐term monitoring of changes in dissolved oxygen (DO) and pH is of great importance to quantifying aquatic ecosystem metabolism, particularly for lakes under the changing global environment. During 173 days, diel DO cycles were measured in situ along with the main driving variables of pH, wind speed (WS), and net solar radiation (Rn) in a temperate shallow lake. Best‐fit multiple non‐linear regression (MNLR) models of diel DO time series were built and validated on a monthly basis, with R2 values ranging from 42.4% in September to 95.4% in November for validation. The strong relationship between diel DO and pH (r = 0.6) appeared to be related to the patterns of ecosystem productivity and respiration, and sensitivity of decomposing bacteria to changes in pH. pH‐driven lake metabolism appears to have significant implications for diel and seasonal lake metabolism in a changing global environment.  相似文献   

18.
A cross‐sectional model, based on the two dimensional groundwater flow equation of Edelman, was applied at seven transects distributed over four geological cross sections to estimate groundwater heads and recharge from/or groundwater discharge to Lake Nasser. The lake with a length of 500 km and an average width of 12 km was created over the period 1964–1970, the time for constructing the Aswan High Dam (AHD). The model, constrained by regional‐scale groundwater flow and groundwater head data in the vicinity of the lake, was successfully calibrated to timeseries of piezometeric heads collected at the cross sections in the period 1965–2004. Inverse modeling yielded high values for the horizontal hydraulic conductivity in the range of 6.0 to 31.1 m day?1 and storage coefficient between 0.01 and 0.40. The results showed the existence of a strong vertical anisotropy of the aquifer. The calibrated horizontal permeability is systematically higher than the vertical permeability (≈1000:1). The calibrated model was used to explore the recharge from/or groundwater discharge to Lake Nasser at the seven transects for a 40‐year period, i.e. from 1965 to 2004. The analysis for the last 20‐year period, 1985–2004, revealed that recharge from Lake Nasser reduced by 37% compared to the estimates for the first 20‐year period, 1965–1984. In the period 1965–2004, seepage of Lake Nasser to the surrounding was estimated at 1.15 × 109 m3 year?1. This led to a significant rise of the groundwater table. Variance‐based sensitivity and uncertainty analysis on the Edelman results were conducted applying quasi‐Monte Carlo sequences (Latin Hypercube sampling). The maximum standard deviation of the total uncertainty on the groundwater table was 0.88 m at Toshka (west of the lake). The distance from the lake, followed by the storage coefficient and hydraulic conductivity, were identified as the most sensitive parameters. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

19.
Tension infiltrometers allow water to infiltrate into the soils at various specified pressure heads. The resulting infiltration rates can then be analysed for soil hydraulic properties by either analytical or inverse numerical methods. Tension infiltrometers however are primarily designed to be deployed on horizontal land surfaces, and their applications have been studied widely using an inverse numerical tool HYDRUS‐2D. However, natural landscapes are often nonhorizontal, and infiltration through tension infiltrometers on sloped surfaces is no longer an axisymmetrical two‐dimensional (2D) process but a fully three‐dimensional (3D) one. In addition, minimal research has examined the effect of simplifying the 3D problem to a 2D one on the hydraulic conductivity estimated using tension infiltrometer data from different land slopes of various soil types. Therefore, in this study, tension infiltrometer data on different slopes have been obtained from a catchment located at National University of Singapore. In addition, tension infiltrometer data of six soil types on different slopes and with different initial water content were simulated using HYDRUS‐3D. Combining field measurements, forward and inverse modelling, the influence of applying a 2D approximation on hydraulic property estimations using tension infiltrometer data was examined. The results show that the estimations for soils with high infiltration rates are more sensitive to application of the 2D approximation. The maximum allowable slopes for employing 2D approximation on clayey and sandy soils are 25° and less than 3°, respectively. Furthermore, the maximum allowable slope decreases with decreasing initial water content. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

20.
Short‐circuiting flow, commonly experienced in many constructed wetlands, reduces hydraulic retention times in unit wetland cells and decreases the treatment efficiency. A two‐dimensional (2‐D), physically based, distributed modelling approach was used to systematically address the effects of bathymetry and vegetation on short‐circuiting flow, which previously have been neglected or lumped in one‐dimensional wetland flow models. In this study, a 2‐D transient hydrodynamics with advection‐dispersion model was developed using MIKE 21 and calibrated with bromide tracer data collected at the Orlando Easterly Wetland Cell 7. The estimated topographic difference between short‐circuiting flow zone and adjacent area ranged from 0·3 to 0·8 m. A range of the Manning roughness coefficient at the short‐circuiting flow zone was estimated (0·022–0·045 s m?1/3). Sensitivity analysis of topographical and vegetative heterogeneity deduced during model calibration shows that relic ditches or other ditch‐shaped landforms and the associated sparse vegetation along the main flow direction intensify the short‐circuiting pattern, considerably affecting 2‐D solute transport simulation. In terms of hydraulic efficiency, this study indicates that the bathymetry effect on short‐circuiting flow is more important than the vegetation effect. Copyright © 2009 John Wiley & Sons, Ltd.  相似文献   

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

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