共查询到20条相似文献，搜索用时 84 毫秒
1.
The structure of macroporous or aggregated soils and fractured rocks is generally so complex that it is impractical to measure the geometry at the microscale (i.e., the size and the shape of soil aggregates or rock matrix blocks, and the myriad of fissures or fractures), and use such data in geometry-dependent macroscale flow and transport models. This paper analyzes a first-order type dual-porosity model which contains a geometry-dependent coefficient, β, in the mass transfer term to macroscopically represent the size and shape of soil or rock matrix blocks. As a reference, one- and two-dimensional geometry-based diffusion models were used to simulate mass transport into and out of porous blocks of defined shapes. Estimates for β were obtained analytically for four different matrix block geometries. Values for β were also calculated by directly matching analytical solutions of the diffusion models for a number of selected matrix block geometries to results obtained with the first-order model assuming standard boundary conditions. Direct matching improved previous results for cylindrical macropore geometries, especially when relatively small ratios between the outer soil mantle and the radius of the inner cylinder were used. Results of our analysis show that β is closely related to the ratio of the effective surface area available for mass transfer, and the soil matrix volume normalized by the effective characteristic length of the matrix system. Using values of β obtained by direct matching, an empirical function is derived to estimate macroscopic geometry coefficients from medium properties which in principle are measurable. The method permits independent estimates of β, thus allowing the dual-porosity approach eventually to be applied to media with complex and mixed types of structural geometry. 相似文献
2.
Michael Stoelzle Markus Weiler Kerstin Stahl Andreas Morhard Tobias Schuetz 《水文研究》2015,29(6):1301-1313
Previous work has shown that streamflow response during baseflow conditions is a function of storage, but also that this functional relationship varies among seasons and catchments. Traditionally, hydrological models incorporate conceptual groundwater models consisting of linear or non‐linear storage–outflow functions. Identification of the right model structure and model parameterization however is challenging. The aim of this paper is to systematically test different model structures in a set of catchments where different aquifer types govern baseflow generation processes. Nine different two‐parameter conceptual groundwater models are applied with multi‐objective calibration to transform two different groundwater recharge series derived from a soil‐atmosphere‐vegetation transfer model into baseflow separated from streamflow data. The relative performance differences of the model structures allow to systematically improve the understanding of baseflow generation processes and to identify most appropriate model structures for different aquifer types. We found more versatile and more aquifer‐specific optimal model structures and elucidate the role of interflow, flow paths, recharge regimes and partially contributing storages. Aquifer‐specific recommendations of storage models were found for fractured and karstic aquifers, whereas large storage capacities blur the identification of superior model structures for complex and porous aquifers. A model performance matrix is presented, which highlights the joint effects of different recharge inputs, calibration criteria, model structures and aquifer types. The matrix is a guidance to improve groundwater model structures towards their representation of the dominant baseflow generation processes of specific aquifer types. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
3.
This paper deals with the derivation of the hydrological response of a hillslope on the assumption of quick runoff by surface runoff generation. By using the simple non‐linear storage based model, first proposed by Horton, an analytical solution of the overland flow equations over a plane hillslope was derived. This solution establishes a generalization for different flow regimes of Horton's original solution, which is valid for the transitional flow regime only. The solution proposed was compared successfully with that of Horton and, for the turbulent flow regime, to the one derived from kinematic wave theory. This solution can be applied easily to both stationary and non‐stationary rainfall excess events. An analytical solution for the instantaneous response function (IRF) was also derived. Finally, simple expressions to compute peak and time to peak of IRF are proposed. Copyright © 2001 John Wiley & Sons, Ltd. 相似文献
4.
Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation 
下载免费PDF全文
下载免费PDF全文 In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field. 相似文献
5.
《Advances in water resources》1996,19(5):317-327
A new lumped-parameter approach to simulating unsaturated flow processes in dual-porosity media such as fractured rocks or aggregated soils is presented. Fluid flow between the fracture network and the matrix blocks is described by a non-linear equation that relates the imbibition rate of the local difference in liquid-phase pressure between the fractures and the matrix blocks. Unlike a Warren-Root-type equation, this equation is accurate in both the early and late time regimes. The fracture/matrix interflow equation has been incorporated into an existing unsaturated flow simulator, to serve as a source/sink term for fracture gridblocks. Flow processes are then simulated using only fracture gridblocks in the computational grid. This new lumped-parameter approach has been tested on two problems involving transient flow in fractured/porous media, and compared with simulations performed using explicit discretisation of the matrix blocks. The new procedure seems to accurately simulate flow processes in unsaturated fractured rocks, and typically requires an order of magnitude less computational time than do simulations using fully-discretised matrix blocks. 相似文献
6.
Hyperconcentrated flows on a forested alluvial fan of eastern Canada: geomorphic characteristics,return period,and triggering scenarios 下载免费PDF全文
下载免费PDF全文 The assessment of the dominant flow type on alluvial fans usually refers to two categories: debris‐flow fans (i.e. sediment gravity flows) and fluvial fans (i.e. fluid gravity flows). Here we report the results of combined morphometric, stratigraphic and sedimentological approaches which suggest that hyperconcentrated flows, a transitional process rheologically distinct from debris flows and floods and sometimes referred to as debris floods, mud floods, or transitional debris flows, are the dominant fan building process in eastern Canada. These flows produce transitional facies between those of debris flows which consist of a cohesive matrix‐supported diamicton, and those of river flows which display more distinct stratification. The size of the blocks in the channels and the abrasion scars at the base of several trees attest to the high transport capacity of these flows. The fan channels are routed according to various obstacles comprised primarily of woody debris that impede sediment transit. However, these conditions of sediment storage are combined with readily available sediment due to the friable nature of the local lithology. Tree‐ring analysis allowed the reconstruction of eight hydrogeomorphic events which are characterized by a return period of 9.25 years for the period 1934–2008, although most of the analyzed events occurred after 1970. Historical weather data analysis indicates that they were related to rare hydrometeorological events at regional and local scales. This evidence led to the elaboration of weather scenarios likely responsible for triggering flows on the fan. According to these scenarios, two distinct hydrologic regimes emerge: the torrential rainfall regime and the nival regime related to snowmelt processes. Hydrogeomorphic processes occurring in a cold‐temperate climate, and particularly on small forested alluvial fans of north‐eastern North America, should receive more attention from land managers given the hazard they represent, as well as because of their sensitivity to various meteorological parameters. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
7.
Hydrological responses vary spatially and temporally according to watershed characteristics. In this study, the hydrological models that we developed earlier for the Little Miami River (LMR) and Las Vegas Wash (LVW) watersheds in the USA were used for detailed sensitivity analyses. To compare the relative sensitivities of the hydrological parameters of these two models, we used normalized root mean square error (NRMSE). By combining the NRMSE index with the flow duration curve analysis, we derived an approach to measure parameter sensitivities under different flow regimes. Results show that the parameters related to groundwater are highly sensitive in the LMR watershed, whereas the LVW watershed is primarily sensitive to near-surface and impervious parameters. The high and medium flows are more impacted by most of the parameters. The low flow regime was highly sensitive to groundwater-related parameters. Moreover, our approach is found to be useful in facilitating model development and calibration.
EDITOR D. KoutsoyiannisASSOCIATE EDITOR S. Huang 相似文献
8.
Peter A. Waldner Martin Schneebeli Ute Schultze‐Zimmermann Hannes Flühler 《水文研究》2004,18(7):1271-1290
Water flow through a melting snow pack modifies its structure and stability and affects the release of water and nutrients into soils and surface waters. Field and laboratory observations indicate a large spatial variability on various scales of the liquid water content and flow, a dominant system feature currently not included in numerical models. We investigated experimentally water and dye tracer movement through microstructurally different snow pack horizons and the persistence of preferential flow paths. Naturally rounded snow of varying grain size was artificially packed to obtain well known conditions by sieving it into rectangular bins. Surface melt was induced with infrared lamps. The flow paths were visualized with tracers and liquid water content was monitored with time domain reflectometry probes. Vertical cuts through the snow pack were imaged. The dye tracer patterns allowed the two flow regimes ‘matrix flow’ and ‘preferential flow’ to be distinguished. Matrix flow is apparently dominated by film and capillary flow in the unsaturated snow matrix. The capillary barrier effect at a boundary between a fine over a coarse textured layer on matrix flow in snow was confirmed. In contrast, preferential flow appears as well‐defined flow fingers that advance from 0·1 to 1 cm s?1. During a melt phase, the advancing flow fingers enlarge and are only partially time invariant. It remains to be shown whether the continuum concept, including the Darcy–Buckingham law is apt to describe the extremely non‐linear nature of water flow and the travel time of solutes in snow under conditions of melt water percolation. Probably, snow packs that include faceted crystals and large variations in bulk density, feature more pronounced capillary barriers and preferential flow triggering, but also stronger impeding of fingers by lateral dispersion. Further, triggering and persistence of preferential flow is complicated by the usually transient infiltration rate. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
9.
How nature uses overland flow regime to maximize equilibrium detention storage and flood attenuation
Tommy S. W. Wong 《水文研究》2008,22(26):5004-5012
Equilibrium detention storage is an important parameter as it has a proportional effect on flood attenuation. In this paper, based on the kinematic wave theory, a working formula for the equilibrium detention storage of an overland plane with upstream inflow has been derived. Since the flow regime over a concrete plane can vary throughout the entire range laminar to turbulent, this case has been selected to examine the effect of flow regime on the equilibrium detention storage. In the examination, the derived formula has been applied to four flow regimes: (a) laminar, (b) transitional, (c) near turbulent, and (d) turbulent. The examination shows that for planes with a small discharge, laminar flow gives the maximum detention storage. For planes with a medium discharge, transitional flow gives the maximum detention storage, and for planes with a large discharge, near turbulent flow gives the maximum detention storage. The flow regime can cause more than two‐fold increase in detention storage. All these results can be attributed to the respective flow resistance, and have been endorsed with analyses of the water surface profile and the rising limb of the hydrograph. Finally, relating the results to real‐life situations, it shows that the flow regime that gives the maximum detention storage is also the dominating flow regime in nature. Hence, extraordinarily, the flow regimes that exist in nature in fact provide maximum flood attenuation. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
10.
Management of open-channel flow systems requires accurate models of flow transfer. This article presents a simple nonlinear model representative of the flow transfer in a river reach. The model is obtained through linearization of a physical model, simplification using the cumulant matching method and analytic identification of a nonlinear model coinciding with the linear model around equilibrium points, corresponding to the hydraulic permanent regimes. The methodology is illustrated on the diffusive wave equation and the Saint-Venant equations. The obtained nonlinear models are compared in simulation to the initial models. The nonlinear model is shown to ensure mass conservation, despite the variable delay element of the model. The proposed model can reproduce the nonlinear behavior of the time-delay with discharge variations. It is well-suited for fast simulations, flow forecasting, and for controller design. 相似文献
11.
《水文科学杂志》2013,58(6):1106-1124
AbstractNatural 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. 相似文献
12.
Carbonate reservoir is one of the important reservoirs in the world. Because of the characteristics of carbonate reservoir, horizontal well has become a key technology for efficiently developing carbonate reservoir. Establishing corresponding mathematical models and analyzing transient pressure behaviors of this type of well-reservoir configuration can provide a better understanding of fluid flow patterns in formation as well as estimations of important parameters. A mathematical model for a oil–water two-phase flow horizontal well in triple media carbonate reservoir by conceptualizing vugs as spherical shapes are presented in this article. A semi-analytical solution is obtained in the Laplace domain using source function theory, Laplace transformation, and superposition principle. Analysis of transient pressure responses indicates that seven characteristic flow periods of horizontal well in triple media carbonate reservoir can be identified. Parametric analysis shows that water saturation of matrix, vug and fracture system, horizontal section length, and horizontal well position can significantly influence the transient pressure responses of horizontal well in triple media carbonate reservoir. The model presented in this article can be applied to obtain important parameters pertinent to reservoir by type curve matching. 相似文献
13.
The fundamental solutions of axisymmetric elastodynamic problem for the multilayered half-space due to an impulsive ring source acting within a layered elastic media are derived in time domain with the aid of Laplace–Hankel mixed transform and transfer matrix techniques. In addition, an effective numerical procedure, which utilizes the fast Hankel transform algorithm, is also proposed to calculate these solutions. Illustrative examples have been given to demonstrate that the fundamental solutions can be readily evaluated and the numerical results are of high accuracy. The present solutions can be directly applied to determine the transient wave fields caused by a seismic source and show the potential application to the elastodynamic problems solved by the boundary element method. 相似文献
14.
AbstractSteep mountainous areas account for 70% of all river catchments in Japan. To predict river discharge for the mountainous catchments, many studies have applied distributed hydrological models based on a kinematic wave approximation with surface and subsurface flow components (DHM-KWSS). These models reproduce observed river discharge of catchments in Japan well; however, the applicability of a DHM-KWSS to catchments with different geographical and climatic conditions has not been sufficiently examined. This research applied a DHM-KWSS to two river basins that have different climatic conditions from basins in Japan to examine the transferability of the DHM-KWSS model structure. Our results show that the DHM-KWSS model structure explained flow regimes for a wet river basin as well as a large flood event in an arid basin; however, it was unable to explain long-term flow regimes for the arid basin case study. 相似文献
15.
An infiltration model based on flow variability in macropores: development, sensitivity analysis and applications 总被引：1，自引：0，他引：1
Simulating infiltration in soils containing macropores still provides unsatisfactory results, as existing models seem not to capture all relevant processes. Recent studies of macropore flow initiation in natural soils containing earthworm channels revealed a distinct flow rate variability in the macropores depending on the initiation process. When macropore flow was initiated at the soil surface, most of the macropores received very little water while a few macropores received a large proportion of the total inflow. In contrast, when macropore flow was initiated from a saturated or nearly saturated soil layer, macropore flow rate variation was much lower. The objective of this study was to develop, evaluate, and test a model, which combines macropore flow variability with several established approaches to model dual permeability soils. We then evaluate the INfiltration–INitiation–INteraction Model (IN3M) to explore the influence of macropore flow variability on infiltration behavior by performing a sensitivity analysis and applying IN3M to sprinkling and dye tracer experiments at three field sites with different macropore and soil matrix properties. The sensitivity analysis showed that the flow variability in macropores reduces interaction between the macropores and the surrounding soil matrix and thus increases bypass flow, especially for surface initiation of macropore flow and at higher rainfall intensities. The model application shows reasonable agreement between IN3M simulations and field data in terms of water balance, water content change, and dye patterns. The influence of macropore flow variability on the hydrological response of the soil was considerable and especially pronounced for soils where initiation occurs at the soil surface. In future, the model could be applied to explore other types of preferential flow and hence to get a generally better understanding of macropore flow. 相似文献
16.
A model describing the three‐dimensional matrix flow along a slope with rock fragments or impermeable blocks was developed. The model was combined with modified Picard's iteration to ensure mass conservation in the unsaturated flow. We found that rock fragments obstruct water flow along the slope. The groundwater table must be raised to provide a sufficient pore water pressure gradient to facilitate water flow, but higher pore water pressure may induce slope failure. We also conducted a bench‐scale laboratory flume experiment to examine the effects of impermeable blocks on downstream seepage flow. In addition, a numerical experiment was conducted to examine how different arrangements of impermeable blocks affect downstream seepage flow and pore water pressure. This research demonstrated that the hydraulic phenomena were affected when impermeable blocks were present, and pore water pressure increased as the position of impermeable blocks was lowered. Copyright © 2006 John Wiley & Sons, Ltd. 相似文献
17.
Hillslopes have complex three‐dimensional shapes that are characterized by their plan shape, profile curvature of surface and bedrock, and soil depth. To investigate the stability of complex hillslopes (with different slope curvatures and plan shapes), we combine the hillslope‐storage Boussinesq (HSB) model with the infinite slope stability method. The HSB model is based on the continuity and Darcy equations expressed in terms of storage along the hillslope. Solutions of the HSB equation account explicitly for plan shape by introducing the hillslope width function and for profile curvature through the bedrock slope angle and the hillslope soil depth function. The presented model is composed of three parts: a topography model conceptualizing three‐dimensional soil mantled landscapes, a dynamic hydrology model for shallow subsurface flow and water table depth (HSB model) and an infinite slope stability method based on the Mohr–Coulomb failure law. The resulting hillslope‐storage Boussinesq stability model (HSB‐SM) is able to simulate rain‐induced shallow landsliding on hillslopes with non‐constant bedrock slope and non‐parallel plan shape. We apply the model to nine characteristic hillslope types with three different profile curvatures (concave, straight, convex) and three different plan shapes (convergent, parallel, divergent). In the presented model, the unsaturated storage has been calculated based on the unit head gradient assumption. To relax this assumption and to investigate the effect of neglecting the variations of unsaturated storage on the assessment of slope stability in the transient case, we also combine a coupled model of saturated and unsaturated storage and the infinite slope stability method. The results show that the variations of the unsaturated zone storage do not play a critical role in hillslope stability. Therefore, it can be concluded that the presented dynamic slope stability model (HSB‐SM) can be used safely for slope stability analysis on complex hillslopes. Our results show that after a certain period of rainfall the convergent hillslopes with concave and straight profiles become unstable more quickly than others, whilst divergent convex hillslopes remain stable (even after intense rainfall). In addition, the relation between subsurface flow and hillslope stability has been investigated. Our analyses show that the minimum safety factor (FS) occurs when the rate of subsurface flow is a maximum. In fact, by increasing the subsurface flow, stability decreases for all hillslope shapes. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
18.
Ren-Shi Nie Ying-Feng Meng Yong-Lu Jia Jian-Lin Shang Yong Wang Jian-Guo Li 《Acta Geophysica》2012,60(1):232-259
The paper deals with unsteady inter-porosity flow modeling of underground fluid in a multiple media reservoir. Assuming spherical
vugs, symmetrically distributed pressure, negligible inter-porosity flow between matrix and vug systems and centrifugal flow
of the fluid from matrix blocks or vugs to fractures, and treating media directly connected with wellbore as the fracture
system, we establish and solve a model of unsteady inter-porosity flow for dual and triple porosity media reservoirs. We provide
simulated graphs of pressure and pressure derivative log-log type curves, and analyze the transient flow process and characteristics
of type curves affected by different parameters. The new type curves of unsteady inter-porosity flow modeling are evidently
different in shape and characteristics from those of pseudo-steady inter-porosity flow modeling. The location of dimensionless
pressure of unsteady inter-porosity is lower than that of pseudo-steady inter-porosity, which indicates that unsteady inter-porosity
flow accelerates an energy supplement during production. Qualitatively, the unsteady inter-porosity flow modeling reduces
the classical V-shaped response. We also estimated parameters from well test data in real applications using this model. 相似文献
19.
20.
基于孔隙介质的Biot理论,首先利用Laplace变换,给出圆柱坐标系下横观各向同性饱和弹性多孔介质在变换域上的波动方程;将波动方程解耦后,根据方位角的Fourier展开和径向Hankel变换,求解了Biot波动方程,得到以土骨架位移、孔隙水压力和土介质总应力分量的积分形式的一般解;借助一般解,建立了有限厚度饱和土层和饱和半空间的精确动力刚度矩阵,并由土层的层间界面连续条件建立三维非轴对称层状饱和地基的总刚度方程;在此基础上,系统研究了横观各向同性饱和半空间体在内部集中荷载激励下的动力响应,并给出了问题的瞬态解答.该研究为运用边界元法求解饱和地基动力响应奠定了理论基础. 相似文献