共查询到20条相似文献，搜索用时 922 毫秒
1.
An analytical solution to 1D coupled water infiltration and deformation in layered soils is derived using a Laplace transformation. Coupling between seepage and deformation, and initial conditions defined by arbitrary continuous pore‐water pressure distributions are considered. The analytical solutions describe the transient pore‐water pressure distributions during 1D, vertical infiltration toward the water table through two‐layer unsaturated soils. The nonlinear coupled formulations are first linearized and transformed into a form that is solvable using a Laplace transformation. The solutions provide a reliable means of comparing the accuracy of various numerical methods. Parameters considered in the coupled analysis include the saturated permeability (ks), desaturation coefficient (α), and saturated volumetric water content (θs) of each soil layer, and antecedent and subsequent rainfall infiltration rates. The analytical solution demonstrates that the coupling of seepage and deformation plays an important role in water infiltration in layered unsaturated soils. A smaller value of α or a smaller absolute value of the elastic modulus of the soil with respect to a change in soil suction (H) for layered unsaturated soils means more marked coupling effect. A smaller absolute value of H of the upper layer soil also tends to cause more marked coupling effect. A large difference between the saturated coefficients of permeability for the top and bottom soil layers leads to reduced rainfall infiltration into the deep soil layer. The initial conditions also play a significant role in the pore‐water pressure redistribution and coupling effect. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
2.
Analytical solution to 1D coupled water infiltration and deformation in unsaturated soils 总被引：1，自引：0，他引：1
An analytical solution to 1D coupled water infiltration and deformation is derived using a Fourier integral transform. Exponential functional forms are used to represent the hydraulic conductivity–pore‐water pressure relationship and the soil‐water characteristic curve. Fredlund's incremental‐linear constitutive model for unsaturated soils is adopted. The analytical solution considers arbitrary initial pore‐water pressure distributions and flux and pressure boundary conditions. The corresponding analytical solutions to coupled steady‐state problems are also obtained. The analytical solutions demonstrate that the coupling of seepage and deformation plays an important role in water infiltration in unsaturated soils. In the early stages of infiltration, the difference between uncoupled and coupled conditions becomes marked over time, and in late stages, the difference caused by the coupling effects diminishes toward the steady state. The difference between the uncoupled and coupled conditions increases with decreasing desaturation coefficient (α). Pore‐water pressure or deformation changes caused by the coupling effects are mainly controlled by the degree of soil volume change due to a change in soil suction (H). The smaller the absolute value of H, the greater the effect of coupling on the infiltration and deformation. The ratio of rainfall intensity to saturated permeability (q/ks) also has a strong influence on the coupled seepage and deformation. Copyright © 2008 John Wiley & Sons, Ltd. 相似文献
3.
In this paper, we analyze the impact of physical and chemical heterogeneity on solute travel time to a pumping well. We consider
a solute undergoing reversible linear instantaneous equilibrium sorption. Both the distribution coefficient, K
d
, and the transmissivity field, T, are considered spatially variable, and are modeled as partially correlated spatial random functions. Groundwater flow and
solute transport are then solved within the context of a numerical Monte Carlo framework. The results are analyzed on the
basis of dimensional analysis techniques. Simple and compact expressions characterizing the dependence of the target travel
time moments on relevant dimensionless groups are proposed. The functional form of these expressions is inspired by, and is
consistent with, the previous works of Sanchez-Vila and Rubin (Water Resour. Res. 39(4):1086, 2003) and Riva et al. (J. Contam. Hydrol. 82:23–43, 2006) A key result is that the effects of the chemical and physical heterogeneities on the mean travel time can be decoupled consistently
with existing analytical results. The relative role of physical and geochemical heterogeneities in travel time variance is
more complex, and such a decoupling is not observed. Potential uses of this work include the assessment of aquifer reclamation
time by means of a single pumping well. 相似文献
4.
Technetium (99Tc) spreads predominantly laterally through US Department of Energy Hanford site sediments. Lateral transport implies that at higher tensions, h, in the unsaturated zone, the effective hydraulic conductivity, K, may be strongly anisotropic. A modeling procedure has been developed to predict characteristics of the Tc plumes. The procedure consists of: (1) Adapting existing numerical techniques based on critical path analysis to calculate K(h), (2) Statistically correlating predicted K at various h values with texture, (3) Seeking value of h, for which anisotropy and horizontal K values are both sufficiently large to accommodate multi-kilometer horizontal spreading, (4) Predicting the distribution of K values for vertical flow as a function of system support volume, (5) Comparing the largest likely K value in the vertical direction with the expected K in the horizontal direction, (6) Finding the length scale at which the two K values are roughly equal, and (7) Comparing that length scale with horizontal spreading of the plume. Predictions of the typical value of h at which spreading is occurring compares well with inference. However, the length scale at which Tc transport in the unsaturated zone changes from predominantly horizontal to more nearly vertical appears underestimated. 相似文献
5.
6.
S. A. Mirbagheri S. A. Hashemi Monfared H. R. Kazemi 《Environmental Earth Sciences》2009,59(2):287-296
A one-dimensional mathematical model based on convection–dispersion equation in unsaturated porous media is presented to compute
inorganic total solid concentration in the soil column under the Shiraz landfill. In addition, a dynamic mathematical model
is formulated to simulate concentrations of ions such as Ca2+
, Mg2+, Fe2+, K+, Na+, Cl−, SO4
2− and HCO3
− as well as PH and EC in soil profile under the Shiraz landfill. Leachw model was applied to simulate water flow, water content
and hydraulic conductivity in soil depth. The model was calibrated and verified by using different sets of data collected
from several segments of soil depth in the study area. The numerical solution obtained using finite element method. The simulated
values for the parameters were compared with measured values as well as analytical solution. The simulated results are in
good agreement with measured values. This model could be applied to field scale problems for the landfill management. 相似文献
7.
An approximate analytical solution is presented for the coupled seepage and deformation problem of unsaturated soils. Because of the matric suction dependence of both saturation and permeability coefficient, the coupled governing equations are strongly nonlinear. To obtain an analytical solution, these coupled governing equations are linearized and analytically solved for a specified saturation using the eigenfunction method. Then, the obtained analytical solutions are extended to the entire saturation range. Comparison between the current solution and the previous theoretical solution indicates that the proposed solution yields excellent results. Due to its analytical nature, the proposed procedure can be effectively used to obtain the solution of the coupled seepage and deformation of unsaturated soils. 相似文献
8.
Amir M. Abdolahzadeh Benoit Lacroix Vachon Alexandre R. Cabral 《Geotechnical and Geological Engineering》2011,29(5):783-802
One important step in the design of inclined covers with capillary barrier effect (CCBE) is the determination of the water
diversion length (DL). Numerical simulations can predict the DL more precisely than steady-state analytical solutions. Nevertheless,
as simplified methods have always been part of engineering design, the application of analytical solutions with conservative
boundary conditions, may allow engineers to make reasonable predictions, particularly during the pre-feasibility stage of
a project. In this study, a CCBE was designed, constructed and instrumented at the Saint-Tite-des-Caps landfill, Quebec, Canada.
This CCBE included a seepage control layer superimposing a sand-gravel capillary barrier. The seepage control layer was made
up of deinking by-products (DBP), an industrial by-product that was previously disposed of as waste. The capillary barrier
was designed using an adaptation of the Ross analytical solution and the scenario considered was that of steady-state flow
during constant seepage flow applied uniformly at the top of the sand-gravel capillary barrier. Although these conditions
appear simplistic, they were deemed reasonable because placement of the seepage control layer on the top of the capillary
barrier led to very low suctions at the interface, thereby allowing uniform downward seepage rates, limited by the saturated
hydraulic conductivity of the DBP. In this paper, a discussion about the behaviour of the cover system based on 4 years of
field data from several instruments is presented. The challenge of using DBP, more precisely the settlement of the DBP layer
and its impact on k
sat
, is also assessed. The DL was reassessed considering the new k
sat
. A discussion on the validity of employing analytical solutions to determine DL is also presented. This paper illustrates
how certain variables affect the design of inclined CCBEs that include a highly compressible material as seepage control layer. 相似文献
9.
将数值计算区域用三角形单元进行离散,并为每个单元构建局部坐标系。局部坐标系的X轴为三角形单元某一条边的方向,局部坐标系的原点为该边的其中一个端点。在局部坐标系下,基于“格林公式”及达西定律推导了单元压力梯度及单元流速的解析表达式,给出了流经单元各棱及各节点的流量计算方法。形成了类似固体弹簧系统的渗流管道网络,建立了管道压差与流量的函数关系。将各单元局部坐标系下求得的流速及流量转换至整体坐标系,并在节点上进行凝聚。通过引入流体体积模量实现了节点渗透压力的显式求解,通过引入节点饱和度实现了非饱和问题的求解。基于局部坐标系的方法具有物理意义明确、求解过程简单等特点。通过在局部坐标系下构建管道压差与管道流量的对应关系,将有限元的渗透刚度矩阵简化为两个管道的渗透刚度值,从而节省了内存,提高了计算效率。4个数值算例的计算结果与理论解基本一致,表明了该方法在求解稳态、非稳态、饱和、非饱和渗流问题时的精度。 相似文献
10.
从华南几个具代表性地点的下三叠统灰岩中,获得了有意义的古地磁学初步结果。这些结果支持华南存在分离板块的设想。根据新近获得的地质与古地磁学证据,推断华南存在4个岩石图板块。它们是扬子(D=232.4°,I=-11.8°,K=28.2,α95=5.5°),湘桂(D=198.0°,I=32.2°,α95=12.9°),华夏(D=88.0°,I=13.0°,K=20.1,α95=9.1°)和海南(上二叠统D=338.2°,I=15.9°,K=10.2,α95=14.3°)板块。4个板块的古地磁极位分散,纬度差明显,并存在碰撞后的局部与区域性的旋转作用,尤其是顺时针旋转。华南早中生代的大地构造可理解为4个来自古特提斯和冈瓦纳的岩石圈碎块与欧亚大陆聚合作用的产物。 相似文献
11.
Tian Kanliang Yang Aoqiu Nie Kangyi Zhang Huili Xu Jin Wang Xiaodong 《Acta Geotechnica》2020,15(9):2681-2689
The complexity of unsaturated soil means that its permeability coefficient and seepage law are important topics in soil mechanics and geotechnical engineering. In this study, a new type of equipment for measuring the seepage of unsaturated soil was used to study the steady seepage of unsaturated loess soil collected from northern Shaanxi province in China. The entire process of the steady seepage of unsaturated soil under different seepage hydraulic gradients was obtained, and the permeability coefficient in the steady seepage state was measured. When the soil water content was low, the soil moisture remained only in small pores or at the junction of soil particles, the effective seepage area of the sample was reduced, and the seepage path was extended. Therefore, the lower the water content, the smaller the permeability coefficient of unsaturated soil and the longer the time required for the seepage to reach a steady state. The matric suction of the unsaturated loess soil was measured under different water contents. With decreasing water content, water remained only in the smaller pores, which provided shrink films with a smaller radius of curvature. This allowed the shrink films to withstand a larger pore air pressure and transmit a smaller pore water pressure, so that the matric suction of the loess soil increased with decreasing soil water content. The Gardner model provided an accurate fit of how the matric suction varied with the volumetric water content of the unsaturated soil. The seepage velocity and seepage hydraulic gradient were correlated well with a linear relationship, indicating that the seepage law of unsaturated loess soil accords with Darcy’s law. The experimental results show that the new unsaturated soil permeameter has good performance and controllability, thereby providing support for studying the steady seepage of unsaturated soil.
相似文献12.
在使用有限元方法求解非饱和土渗流问题时,土-水特征曲线和渗透率函数的强烈非线性经常会造成计算中出现迭代不收敛、计算误差大等问题。基于变量变换的思想,结合时间步长自适应技术提出了一种求解非饱和渗流问题的新方法--欠松弛RFT变换方法（ATUR1）。ATUR1方法通过变量变换,大大降低了Richards方程中未知数在空间和时间上的非线性程度,从而改善这种非线性所带来的计算收敛困难和精度差等问题。欠松弛技术的引入减少了迭代过程中的振荡现象,进一步提高了非线性迭代计算的效率。时间步长自适应技术则有效地控制整个计算过程的误差。数值算例结果说明,ATUR1可以有效地提高计算效率和精度,是一种准确有效的计算方法。 相似文献
13.
非饱和-非稳定渗流条件下的边坡临界滑动场 总被引：1，自引：0，他引：1
受季节性降雨或水库运行的影响,岸坡外水位及坡内孔隙水压力场的变化较大,不利于岸坡的稳定性。在水位变化过程中,利用非饱和-非稳定渗流有限元计算得到孔隙水压力场,基于非饱和土的渗流和抗剪强度理论,对水位变化过程中的边坡临界滑动场法进行改进,提出可考虑水位变化与岸坡非饱和-非稳定渗流过程的边坡临界滑动场数值模拟方法。将改进后的水位变化过程中的边坡临界滑动场法分别应用于黏土、粉土岸坡在水位升降过程中的稳定性分析,研究了水位升降速率及基质吸力对岸坡稳定性的影响,并揭示了边坡在水位变化过程中的稳定性变化历程。研究表明,该方法计算结果合理、可靠,更适用于涉水边坡的稳定性计算,且岸坡稳定性变化历程受水位升降速率、基质吸力等多种因素共同影响,只有在考虑非稳定渗流的基础上同时考虑基质吸力的作用才能正确得出水位变化过程中岸坡稳定性变化规律和实质。 相似文献
14.
Existing solutions for analyzing one-dimensional (1-D) consolidation of unsaturated soil are only derived to cater to two extreme drainage conditions (fully drained and undrained). This study presents a new explicit solution for 1-D consolidation of unsaturated soil with semi-permeable drainage boundary. Based on the assumptions of two independent stress variables and the governing equations proposed by Fredlund, the eigenfunction expansion method is adopted to develop an explicit analytical solution to calculate excess pore-water and pore-air pressures in an unsaturated soil when it is subjected to external loads. The developed general solutions are expressed in terms of depth, z, and time, t. For the semi-permeable drainage boundary, eigenvalues and eigenfunctions in the space domain are developed. The technique of Laplace transform is used to solve the coupled ordinary differential equations in the time domain. The newly derived explicit solution is verified with the existing semi-analytical method in the literature, and an excellent agreement is obtained. Compared with the semi-analytical solution, the newly derived analytical solution is more straightforward and explicit so that this solution is relatively easier to be implemented into a computer program to carry out a preliminary assessment of 1-D consolidation of unsaturated soil. 相似文献
15.
A design procedure is proposed to minimize water infiltration into landfills by optimizing the water diversion length of inclined
covers with capillary barrier effect (CCBE). This design procedure is based on a conceptual, mathematical and numerical approach
and aims at selecting materials and optimizing layer thickness. Selection among candidate materials is made based on their
hydraulic conductivity functions and on a threshold infiltration rate imposed on the designer. The capillary break layer (CBL;
bottom layer) is characterized by a weak capillarity, while the moisture retention layer (MRL; upper layer) is characterized
by a compromise between strong capillarity and high hydraulic conductivity. The thickness of the CBL corresponds to the height
where suction reaches its maximum value for a given infiltration rate. This height can be calculated using the Kisch [Géotechnique 9 (1959)] model. The optimal thickness of the MRL is determined by applying an adaptation of the Ross [Water Resources Research 26 (1990)] model. The results obtained using the proposed design procedure were compared to those obtained from numerical simulations
performed using a finite element unsaturated seepage software. The procedure was applied for two cover systems; one where
deinking by-products (DBP) were used as MRL and sand as CBL and another where sand was used as MRL and gravel as CBL. Using
this procedure, it has been shown that an infiltration control system composed of thin layers of sand over gravel is highly
efficient in terms of diversion length and that its efficiency can be enhanced by placing a hydraulic barrier – such as a
layer of DBP – above the MRL. 相似文献
16.
在非饱和多孔介质渗流分析中,近饱和条件下物理模型与数值模型之间的差异会导致数值不稳定问题。为解决这一问题,并保证模拟结果的可靠性,提出了3种方法,并在有限元分析程序U-DYSAC2中分别进行了程序代码的实施。通过数值试验与试验数据的比较,证实了在近饱和条件下土-水特征曲线和水力传导函数的高非线性可引起数值收敛性、稳定性和精度问题,而且在不同条件下含水率和基质吸力的预测结果差异明显。在3种方法中,修正的Van Genuchten模型（MVGM）方法对含水率的预测较为准确,而Line方法对基质吸力的预测较为合理。因此,解决在分析近饱和条件下非饱和多孔介质渗流问题时,为获得接近真实的模拟结果,采用合适的数值方法进行预测是非常关键的。 相似文献
17.
Numerical evaluation of multicomponent cation exchange reactive transport in physically and geochemically heterogeneous porous media 总被引：1，自引：0，他引：1
Experimental evidence and stochastic studies strongly show that the transport of reactive solutes in porous media is significantly
influenced by heterogeneities in hydraulic conductivity, porosity, and sorption parameters. In this paper, we present Monte
Carlo numerical simulations of multicomponent reactive transport involving competitive cation exchange reactions in a two-dimensional
vertical physically and geochemically heterogeneous medium. Log hydraulic conductivity, log K, and log cation exchange capacity (log CEC) are assumed to be random Gaussian functions with spherical semivariograms. Random
realizations of log K and log CEC are used as input data for the numerical simulation of multicomponent reactive transport with CORE2D, a general purpose reactive transport code. Longitudinal features of the fronts of reactive and conservative species are
computed from the temporal and spatial moments of depth-averaged concentrations. Monte Carlo simulations show that: (1) the
displacement of reactive fronts increases with increasing variance of log K, while it decreases with the variance of log CEC; (2) second-order spatial moments increase with increasing variances of
log K and log CEC; (3) uncertainties in the mean arrival time are largest (smallest) for negatively (positively) correlated log
K and Log CEC; (4) cations undergoing competitive cation exchange exhibit different apparent velocities and retardation factors
due to both physical and geochemical heterogeneities; and (5) the correlation between log K and log CEC affects significantly apparent cation retardation factors in heterogeneous aquifers. 相似文献
18.
This study investigates the influence of the water level fluctuation on the stability of soil slopes using coupled seepage and slope stability analysis. A simulation framework was proposed and implemented seamlessly using Python code to seek insights into three factors that have not been thoroughly studied for this issue: soil unit weight variation in the unsaturated zone, unsaturated shear strength models, and velocity of water drawdown. For this purpose, the seepage analysis was carried out by discretizing a numerical seepage analysis model using a finite element analysis platform, FEniCS. The output of the seepage analysis, i.e., pore water pressure distribution, was used as input for the slope stability analysis. Limit equilibrium methods including the Bishop Simplified method and the Ordinary Method of Slices were modified to take into consideration the unsaturated shear strength, unit weight variation in the unsaturated zone, and hydrostatic pressure changes in response to the water level fluctuation of a reservoir. Both seepage and slope analysis modules were validated against commercial programs. Analysis results obtained with the validated framework clearly revealed the distinct influences of the three factors in representative silty and sandy slopes.
相似文献19.
在水电、石油、交通等部门的地下工程中,经常会遇到地下水的渗流问题,如果处理不当,可能诱发应力恶化而影响结构稳定。本文以非饱和介质的渗流模型为基础,实现了地下洞室开挖过程中的非稳定渗流场解析。
按分期施工方案,在保持渗流场变化连续性的前提下,对各施工阶段进行了准稳定渗流分析,并在此基础上建立了地下工程非稳定渗流数值模型。其中,对于洞室系统,基于水量平衡原理,采用了子结构模型;对于防渗帷幕,采用了隔水层模型;对于排水井列,采用了杆单元模型。
根据上述模型,采用有限差分法作为数值模拟方法,编写了地下工程非稳定渗流场分析的有限元程序,实现了复杂渗流区域的渗流场的动态模拟。工程实例计算分析表明,用该模型计算的结果是合理的。 相似文献
按分期施工方案,在保持渗流场变化连续性的前提下,对各施工阶段进行了准稳定渗流分析,并在此基础上建立了地下工程非稳定渗流数值模型。其中,对于洞室系统,基于水量平衡原理,采用了子结构模型;对于防渗帷幕,采用了隔水层模型;对于排水井列,采用了杆单元模型。
根据上述模型,采用有限差分法作为数值模拟方法,编写了地下工程非稳定渗流场分析的有限元程序,实现了复杂渗流区域的渗流场的动态模拟。工程实例计算分析表明,用该模型计算的结果是合理的。 相似文献
20.
This paper proposed a numerical formulation for unsaturated flow problems with nonlinear boundaries of seepage face and soil–atmosphere interface via the concept of parabolic variational inequality (PVI) method. A unified unilateral boundary condition was first proposed to represent the conditions on the seepage face and soil–atmosphere interface boundaries within the partial differential equation (PDE) formulation. A PVI formulation mathematically equivalent to the PDE formulation was then proposed, which automatically transforms the flux part of the unified unilateral boundary condition into the natural boundary condition and eliminates the singularity at seepage points. By discretizing the PVI formulation, a finite element procedure together with an iterative algorithm was suggested. An existing experiment of unsaturated flow in a layered hillside and a laboratory test of unsaturated flow through sand flume performed in this study were used to validate the proposed method, with a good agreement between the measured and computed results and a satisfactory balance of mass being maintained during the simulations. The numerical results also indicated that the problem of mesh dependence associated with unsaturated flow simulations is well addressed with the proposed numerical method. Finally, the process of unsaturated flow in a soil slope with layers of horizontal drains subjected to rainfall/evaporation was further examined. The numerical results reveal that the deployment of drains in a soil slope can significantly lower the pore water pressure around the drains, with the bottom layer drains being most effective in controlling the seepage flow. 相似文献