自适应时间步长法在土体冻结水热耦合模型中的应用

由于相变的存在,土体冻结过程中的温度传导与水分迁移是一个复杂的物理过程。为了更好地描述冻结过程中水分与温度的变化规律,通过对不饱和土体水分传导方程的研究,考虑冻结过程中的相变,建立了一维冻土水热耦合模型。给出了相应的差分与有限元程序,并对室内冻结实验进行了模拟。提出误差因子的概念,通过对程序计算中时间步长与计算用时、误差关系的分析,论证了进行时间步长优化的必要性。在两种不同数值方法的对比中,体现了有限元计算的稳定性。提出了调整后的自适应时间步长计算方法。计算结果表明,优化时间步的自适应步长法,在不影响模型计算准确度的前提下,可以大幅减少计算用时,提高计算效率。     

基于双曲线模型的冻土中锚杆荷载传递特性研究

为研究不同冻结温度及含水率对冻结粉土中锚杆抗拔性能的影响,根据冻土与混凝土接触面的力学特性和变形规律,采用双曲线模型描述锚杆-冻土接触面的剪切特性。基于荷载传递法,建立考虑温度和含水率影响的锚杆荷载传递方程,采用有限差分法进行求解,得到锚固段的剪应力、轴力及承载力计算式;通过ABAQUS数值模拟验证荷载传递方程的合理性,结合算例分析表明：冻土中锚杆抗拔承载力大于常温土,且温度越低承载力越高;在相同荷载作用下,冻土温度越低,轴力沿深度衰减越快,剪应力分布越不均匀;相同冻结温度下,锚杆承载力随含水率的增大呈先增大后减小的趋势;增加锚固体直径能有效地提高锚杆的极限承载力。     

The effect of seasonally frozen soil on stochastic response of elevated water tank under random excitation

Significant seismic events have occurred around the world during winter months in regions where cold temperatures cause ground freezing. Current seismic design practice does not address the effects of cold temperatures in the seasonally frozen areas. Since many elevated water tank structures in cold regions are located in seismic active zones, determining the effect of seasonally frozen soil on the stochastic response of elevated water tank structures subjected to random seismic excitation is an important structural consideration. A three dimensional finite element model, which considers viscous boundaries, was built up to obtain the stochastic seismic behavior of an elevated water tank–fluid–soil interaction system for frozen soil condition. For this model, the power spectral density function represents random ground motion applied to each support point of the three dimensional finite element model of the elevated water tank–fluid–soil interaction system. Numerical results show that the soil temperature affects the seismic response of the elevated water tank; whereas the variation in the thickness of the frozen soil causes insignificant changes on the response. In addition, the effect of the variation in water tank’s fullness on the stochastic response of the coupled system is investigated in the study. As a result, the seasonal frost changes the foundation soil stiffness and may impact seismic behavior of the water tank.     

离石加筋黄土动弹性模量变化规律研究

针对目前加筋土动弹性模量演化机理不明确、测试过程复杂等问题,结合GDS动三轴,开展不同因素(围压、含水率、加筋材料、加筋层数)影响下加筋土动弹性模量测试试验。测试结果表明:(1)加筋土的动弹性模量受含水率、围压、加筋材料和加筋层数等多种因素影响,其影响的敏感性等级为含水率围压加筋层数加筋材料;(2)加筋土的动弹性模量主要受筋材与土颗粒的联结作用和排列方式影响;(3)提出不同因素影响下动弹性模量衰减的预测模型,其表达式符合双曲线规律。研究结论对路基设计、处治具有一定的指导意义。     

Dynamic shear modulus and damping ratio of frozen compacted sand subjected to freeze–thaw cycle under multi-stage cyclic loading

Frozen soil plays an important role on the stability of railway and highway subgrade in cold regions. However, the dynamic properties of frozen soil subjected to the freeze–thaw cycles have rarely been investigated. In this study, cryogenic cyclic triaxial tests were conducted on frozen compacted sand from Nehe, Heilongjiang Province in China which was subjected to the closed-system freeze–thaw cycles. A modified Hardin hyperbolic model was suggested to describe the backbone curves. Then, dynamic shear modulus and damping ratio versus cyclic shear strain were analyzed under the different freeze–thaw cycles, temperatures, initial water contents, loading frequencies and confining pressures. The results indicate that the freeze–thaw process plays a significant effect on the dynamic shear modulus and damping ratio, which slightly change after one freeze–thaw cycle. Dynamic shear modulus increases with increasing initial water content, temperature, loading frequency and confining pressure. Damping ratio increases with increasing initial water content, while decreases with increasing temperature and loading frequency. The effect of confining pressure on the damping ratio was found not significant. Furthermore, the empirical expressions were formulated to estimate dynamic shear modulus and damping ratio of the frozen compacted sand. The results provide guidelines for evaluating the infrastructures in cold regions.     

冻土蠕变过程中结构的CT分析

ＣＴ为非破坏的持续检测和定量描述不同实验条件下试样的内部结构提供了可能，本文利用这一手段观测分析了冻土蠕变过程结构的变化情况，认为冻土蠕变过程中所进行的发育微裂隙、颗粒集合体的破坏以及其它结构缺陷的增生与扩展，制约着土结构的强化与弱化作用，控制着蠕变过程的形态特征，同时还给出了ＣＴ值与物密度和含水量之间的定量关系。     

Laboratory study of infiltration into two frozen engineered (sandy) soils recommended for bioretention

Infiltration of water into two frozen engineered soils of different gradation was studied in laboratory soil columns 1.2 m long and 0.1 m in diameter. Prior to testing, the soil moisture was adjusted to two levels, described by the gravimetric water content of 5% or 10%, and soils were compacted to about 80–90% of the maximum dry density and refrigerated to temperatures ranging from ?8 to ?2 °C. Water with temperatures 8–9 °C was thereafter fed on the top of columns at a constant head, and the times of water breakthrough in the column and reaching a steady percolation rate, as well as the percolation rate, were recorded. The soil water content was a critical factor affecting the thawing process; during freezing, soil moisture was converted into ice, which blocked pores, and its melting required high amounts of energy supplied by infiltrating water. Hence, the thawing of soils with higher initial water content was much slower than in lower moisture soils, and water breakthrough and the attainment of steady percolation required much longer times in higher moisture soils. Heat transfer between infiltrating water, soil ice, and frozen soil particles was well described by the energy budget equations, which constitute a parsimonious model of the observed processes. The finer grained soil and more compacted soil columns exhibited reduced porosity and required longer times for soil thawing. Practical implications of study results for design of bioretention facilities (BFs) in cold climate include the use of coarse engineered soils and fitting bioretention facilities with a drain facilitating soil drainage before the onset of freezing weather. Copyright © 2015 John Wiley & Sons, Ltd.     

Effects of antecedent moisture and macroporosity on infiltration and water flow in frozen soil

Infiltration into frozen soil plays an important role in soil freeze–thaw and snowmelt-driven hydrological processes. To better understand the complex thermal energy and water transport mechanisms involved, the influence of antecedent moisture content and macroporosity on infiltration into frozen soil was investigated. Ponded infiltration experiments on frozen macroporous and non-macroporous soil columns revealed that dry macroporous soil produced infiltration rates reaching 10³ to 10⁴ mm day⁻¹, two to three orders of magnitude larger than dry non-macroporous soil. Results suggest that rapid infiltration and drainage were a result of preferential flow through initially air-filled macropores. Using recorded flow rates and measured macropore characteristics, calculations indicated that a combination of both saturated flow and unsaturated film flow likely occurred within macropores. Under wet conditions, regardless of the presence of macropores, infiltration was restricted by the slow thawing rate of pore ice, producing infiltration rates of 2.8 to 5.0 mm day⁻¹. Reduced preferential flow under wet conditions was attributed to a combination of soil swelling, due to smectite-rich clay (that reduced macropore volume), and pore ice blockage within macropores. In comparison, dry soil column experiments demonstrated that macropores provided conduits for water and thermal energy to bypass the frozen matrix during infiltration, reducing thaw rates compared with non-macroporous soils. Overall, results showed the dominant control of antecedent moisture content on the initiation, timing, and magnitude of infiltration and flow in frozen macroporous soils, as well as the important role of macropore connectivity. The study provides an important data set that can aid the development of hydrological models that consider the interacting effects of soil freeze–thaw and preferential flow on snowmelt partitioning in cold regions.     

Evaluating model of frozen soil environment change under engineering actions

The change of frozen soil environment is evaluated by permafrost thermal stability, thermal thaw sensibility and surface landscape stability and the quantitatively evaluating model of frozen soil environment is proposed in this paper. The evaluating model of frozen soil environment is calculated by 28 ground temperature measurements along Qinghai-Xizang Highway. The relationships of thermal thaw sensibility and freezing and thawing processes and seasonally thawing depth, thermal stability and permafrost table temperature, mean annual ground temperature and seasonally thawing depth, and surface landscape stability and freezing and thawing hazards and their forming possibility are analyzed. The results show that thermal stability, thermal thaw sensibility and surface landscape stability can be used to evaluate and predict the change of frozen soil environment under human engineering action.     

Using a H∞ filter assimilation procedure to estimate root zone soil water content

Root zone soil water content impacts plant water availability, land energy and water balances. Because of unknown hydrological model error, observation errors and the statistical characteristics of the errors, the widely used Kalman filter (KF) and its extensions are challenged to retrieve the root zone soil water content using the surface soil water content. If the soil hydraulic parameters are poorly estimated, the KF and its extensions fail to accurately estimate the root zone soil water. The H‐infinity filter (HF) represents a robust version of the KF. The HF is widely used in data assimilation and is superior to the KF, especially when the performance of the model is not well understood. The objective of this study is to study the impact of uncertain soil hydraulic parameters, initial soil moisture content and observation period on the ability of HF assimilation to predict in situ soil water content. In this article, we study seven cases. The results show that the soil hydraulic parameters hold a critical role in the course of assimilation. When the soil hydraulic parameters are poorly estimated, an accurate estimation of root soil water content cannot be retrieved by the HF assimilation approach. When the estimated soil hydraulic parameters are similar to actual values, the soil water content at various depths can be accurately retrieved by the HF assimilation. The HF assimilation is not very sensitive to the initial soil water content, and the impact of the initial soil water content on the assimilation scheme can be eliminated after about 5–7 days. The observation interval is important for soil water profile distribution retrieval with the HF, and the shorter the observation interval, the shorter the time required to achieve actual soil water content. However, the retrieval results are not very accurate at a depth of 100 cm. Also it is complex to determine the weighting coefficient and the error attenuation parameter in the HF assimilation. In this article, the trial‐and‐error method was used to determine the weighting coefficient and the error attenuation parameter. After the first establishment of limited range of the parameters, ‘the best parameter set’ was selected from the range of values. For the soil conditions investigated, the HF assimilation results are better than the open‐loop results. Copyright © 2010 John Wiley & Sons, Ltd.     

黄土高原大理河流域水沙耦合模型应用研究

水沙模型是定量描述水沙关系及水沙规律的重要工具,现阶段国内外对于水沙模型的研究大都为基于某个典型流域的经验统计模型或基于流域大量基础资料的物理模型,极大限制了其使用范围及模拟精度.本文建立了结构与参数均具有物理意义的流域水沙耦合物理概念模型,其优点是物理概念清楚,模拟精度高,实用性强,易于深入研究泥沙基本规律.该模型将概念性水文模型和泥沙模型耦合,提出水流挟沙能力和土壤抗侵蚀能力概念,用对数曲线近似描述流域土壤抗侵蚀能力的空间变异性,在拜格诺河道水流悬移质泥沙公式基础上建立概念性沟蚀产沙公式,参照水流汇集相似性建立泥沙汇集演算公式.选取黄河中游大理河流域4个不同流域尺度的实际流域对模型进行应用检验,模拟结果表明,该模型的水流泥沙两部分均有很高的模拟精度,可以很好地模拟黄河中游地区不同流域尺度和年际尺度上的洪水过程和泥沙产生及输移过程,充分证明了该模型结构、参数和计算方法的合理性,可用于定量分析流域内各项水保措施的减水减沙效益及水沙关系变化趋势,对该模型的推广应用可做进一步分析研究.     

寒区黏土与结构接触面冻结强度特性试验研究

寒区土与结构接触面冻结强度可以视为是与土性、温度、含水量、界面粗糙度、法向压力等诸多因素直接相关的复杂函数形式,并直接影响到上部结构的承载能力及稳定性。以青藏高原黏土与不同粗糙度的钢板结构接触面为研究对象,通过不同含水率、不同温度及不同法向压力下冻结黏土与钢板结构接触面的直剪正交试验,研究土与结构接触面冻结强度的影响因素及影响程度大小,并对冻结强度的变化规律进行初步分析和探讨。研究结果表明含水量对抗剪强度的影响最大,在不同影响因素共同作用下,界面抗剪强度最小值为0.13 MPa,最大值为0.45 MPa。界面抗剪强度随含水量的增加、温度的升高而明显降低,随界面法向压力和粗糙度的增大呈明显增大趋势。界面强度基本可通过摩尔-库仑准则,利用界面的黏聚力和内摩擦角进行表示,并在文中给出了界面强度的参考值。     

Study of water tension differences in heterogeneous sandy soils using surface ERT

Herbaceous vegetation in the Sahel grows almost exclusively on sandy soils which preferentially retain water through infiltration and storage. The hydrological functioning of these sandy soils during rain cycles is unknown. One way to tackle this issue is to spatialize variations in water content but these are difficult to measure in the vadose zone. We investigated the use of Electrical Resistivity Tomography (ERT) as a technique for spatializing resistivity in a non-destructive manner in order to improve our knowledge of relevant hydrological processes. To achieve this, two approaches were examined. First, we focused on a possible link between water tension (which is much easier to measure in the field by point measurements than water content), and resistivity (spatialized with ERT). Second, because ERT is affected by solution non-uniqueness and reconstruction smoothing, we improved the accuracy of ERT inversion by comparing calculated solutions with in-situ resistivity measurements. We studied a natural microdune during a controlled field experiment with artificial sprinkling which reproduced typical rainfall cycles. We recorded temperature, water tension and resistivity within the microdune and applied surface ERT before and after the 3 rainfall cycles. Soil samples were collected after the experiment to determine soil physical characteristics. An experimental relationship between water tension and water content was also investigated. Our results showed that the raw relationship between calculated ERT resistivity and water tension measurements in sand is highly scattered because of significant spatial variations in porosity. An improved correlation was achieved by using resistivity ratio and water tension differences. The slope of the relationship depends on the soil solution conductivity, as predicted by Archie's law when salted water was used for the rain simulation. We found that determining the variations in electrical resistivity is a sensitive method for spatializing the differences in water tension which are directly linked with infiltration and evaporation/drainage processes in the vadose zone. However, three factors complicate the use of this approach. Firstly, the relation between water tension and water content is generally non-linear and dependent on the water content range. This could limit the use of our site-specific relations for spatializing water content with ERT through tension. Secondly, to achieve the necessary optimization of ERT inversion, we used destructive resistivity measurements in the soil, which renders ERT less attractive. Thirdly, we found that the calculated resistivity is not always accurate because of the smoothing involved in surface ERT data inversion. We conclude that further developments are needed into ERT image reconstruction before water tension (and water content) can be spatialized in heterogeneous sandy soils with the accuracy needed to routinely study their hydrological functioning.     

黏性土中结合水测定及影响因素分析

土中结合水是形成黏性土的物理、化学、力学等性质的重要因素。通过对一系列人工混合黏性土和2种天然红黏土进行比表面积测试和等温吸附试验,确定各类黏性土的比表面积和土中不同结合水类型及含量,在此基础上探讨黏性土中结合水与矿物成分、比表面积等影响因素之间的定量关系。研究结果表明：(1)黏性土的比表面积与黏土矿物含量线性正相关;(2)等温吸附法中,相对湿度RH=0.90和0.98分别为黏性土中强结合水、弱结合水、自由水的界限划分点;(3)相同矿物成分的混合黏性土吸附结合水含量随蒙脱土含量、比表面积的增加而递增,大致呈线性关系;(4)天然红黏土的吸附结合水含量与其矿物成分总量、比表面积具有正相关关系,黏土矿物含量是最根本的影响因素。     

Coupling of a simultaneous heat and water model with a distributed hydrological model and evaluation of the combined model in a cold region watershed

Snow and frozen soil prevail in cold regions worldwide, and the integration of these processes is crucial in hydrological models. In this study, a combined model was developed by fully coupling a simultaneous heat and water model with a geomorphologically based distributed hydrological model. The combined model simulates vertical and lateral water transfer as well as vertical heat fluxes and is capable of representing the effects of frozen soil and snowmelt on hydrological processes in cold regions. This model was evaluated by using in situ observations in the Binggou watershed, an experimental watershed for cold region hydrology of the Watershed Allied Telemetry Experimental Research Project. Results showed that the model was able to predict soil freezing and thawing, unfrozen soil water content, and snow depth reasonably well. The simulated hydrograph was in good agreement with the in situ observation. The Nash–Sutcliffe coefficient of daily discharge was 0.744 for the entire simulation period, 0.472 from April to June, and 0.711 from June to November. This model can improve our understanding of hydrological processes in cold regions and assess the impacts of global warming on hydrological cycles and water resources. Copyright © 2012 John Wiley & Sons, Ltd.     

Investigating plant root effects on soil electrical conductivity: An integrated field monitoring and statistical modelling approach

Plants have been shown to affect soil water content and temperature. Previous studies were conducted mainly in forestry and agricultural soils, where conditions of soil and vegetation are different from those in an urban landscape. In an urban landscape, the influence of plant roots on electrical conductivity, soil water content and temperature is still not clear. This study aims to investigate the effects of soil water content and temperature on electrical conductivity in vegetated soils through an integrated field monitoring and computational modelling approach. A new relationship between soil electrical conductivity and water content as well as temperature is proposed. Field monitoring was conducted in both vegetated (tree species) and bare soils. The monitoring included measurements of soil water content, soil temperature and soil electrical conductivity. This was followed by response surface regression modelling. Measured results show that soil temperature at shallow depths was lower in vegetated soil than that in the bare soil. This observation was also consistent with the higher soil water content and hence, higher electrical conductivity under tree canopy. The model developed could predict nonlinear relationships between electrical conductivity and soil temperature and water content. Uncertainty analysis indicated normal distribution for electrical conductivity under variation of soil temperature and water content. © 2018 John Wiley & Sons, Ltd.     

大气作用下膨胀土地基的水分迁移与胀缩变形分析

运用土体渗流和蒸发理论,建立了大气-非饱和土相互作用模型;以现场观测的气象数据作为边界条件,进行了地基土中水分迁移的数值模拟,得到了大气作用下地基土体含水量的动态分布规律。计算结果表明,地基土中含水量变化幅度随深度增加而递减,3.5 m深度以下土体的体积含水量基本不变,从而确定了南宁地区膨胀土地基的大气影响层深度为3.5 m。在此基础上,结合已有膨胀土胀缩性指标的干湿循环效应研究成果,提出了一种同时考虑干湿循环效应和1.0 m深处含水量变化的膨胀土地基胀缩变形计算方法,通过算例将该法与传统方法进行比较,结果显示该法更加符合工程实际。     

Infiltration into frozen soil: From core‐scale dynamics to hillslope‐scale connectivity

Infiltration into frozen soil is a key hydrological process in cold regions. Although the mechanisms behind point‐scale infiltration into frozen soil are relatively well understood, questions remain about upscaling point‐scale results to estimate hillslope‐scale run‐off generation. Here, we tackle this question by combining laboratory, field, and modelling experiments. Six large (0.30‐m diameter by 0.35‐m deep) soil cores were extracted from an experimental hillslope on the Canadian Prairies. In the laboratory, we measured run‐off and infiltration rates of the cores for two antecedent moisture conditions under snowmelt rates and diurnal freeze–thaw conditions observed on the same hillslope. We combined the infiltration data with spatially variable data from the hillslope, to parameterise a surface run‐off redistribution model. We used the model to determine how spatial patterns of soil water content, snowpack water equivalent (SWE), and snowmelt rates affect the spatial variability of infiltration and hydrological connectivity over frozen soil. Our experiments showed that antecedent moisture conditions of the frozen soil affected infiltration rates by limiting the initial soil storage capacity and infiltration front penetration depth. However, shallow depths of infiltration and refreezing created saturated conditions at the surface for dry and wet antecedent conditions, resulting in similar final infiltration rates (0.3 mm hr^?1). On the hillslope‐scale, the spatial variability of snowmelt rates controlled the development of hydrological connectivity during the 2014 spring melt, whereas SWE and antecedent soil moisture were unimportant. Geostatistical analysis showed that this was because SWE variability and antecedent moisture variability occurred at distances shorter than that of topographic variability, whereas melt variability occurred at distances longer than that of topographic variability. The importance of spatial controls will shift for differing locations and winter conditions. Overall, our results suggest that run‐off connectivity is determined by (a) a pre‐fill phase, during which a thin surface soil layer wets up, refreezes, and saturates, before infiltration excess run‐off is generated and (b) a subsequent fill‐and‐spill phase on the surface that drives hillslope‐scale run‐off.     

Effect of macropores on soil freezing and thawing with infiltration

An understanding of heat transport and water flow in unsaturated soils experiencing freezing and thawing is important when considering hydrological and thermal processes in cold regions. Macropores, such as cracks, roots, and animal holes, provide efficient conduits for enhanced infiltration, resulting in a unique distribution of water content. However, the effects of macropores on soil freezing and thawing with infiltration have not been well studied. A one‐directional soil‐column freezing and thawing experiment was conducted using unsaturated sandy and silt loams with different sizes and numbers of macropores. During freezing, macropores were found to retard the formation of the frozen layer, depending on their size and number. During thawing, water flowed through macropores in the frozen layer and reached the underlying unfrozen soil. However, infiltrated water sometimes refroze in a macropore. The ice started to form at near inner wall of the macropore, grew to the centre, and blocked flow through the macropore. The blockage ice in the macropore could not melt until the frozen layer disappeared. Improving a soil freezing model to consider these macropore effects is required. Copyright © 2016 John Wiley & Sons, Ltd.