Anisotropic frost heaving coefficient of saturated permafrost considering moisture migration process based on standards
-
摘要:
正冻土冻胀是寒区工程产生冻害的关键因素,其冻胀过程是水热力相互耦合的动态作用结果,在开放系统中,温度、温度梯度、含水率、水分补给强度等都是影响正冻土冻胀变形的重要因素。冻土冻胀是水分迁移产生的竖直方向分凝冻胀和原位冻胀的共同作用,其冻胀力学特性属于各向异性。本研究参考规范内土体含冰量随冻结的变化过程,考虑泊松比、地下水位深度、降温速率等因素,得到正冻土的在冻结过程中水平与竖直方向的冻胀系数的计算方法,通过对比粉土和粉质黏土的冻胀系数,与试验结果吻合较好。案例中粉土在−0.2~−3 °C、0.2~1 m范围内竖向冻胀系数为−1.37×10−3~−7.67×10−3,水平向冻胀系数为−0.81×10−3~−4.85×10−3,差值百分比为10.4%~77.7%,说明考虑分凝冻胀产生的各向异性是必要的。本研究提出的水平与竖直方向的冻胀系数计算方法,可以为科研和设计工作提供参考依据。
Abstract:Frost heave in positive permafrost is a critical factor leading to frost damage of engineering in the cold region. The frost heave process results from the dynamic interaction of hydro-thermal forces. In an open system, factors such as temperature, temperature gradient, water content, and water supply intensity significantly influence the deformation of positive permafrost frost heave. Frost heave in permafrost is the combined effect of vertical segregated frost heave due to moisture migration and in-situ frost heave, and its mechanical characteristics of frost heave are anisotropic. This study, referencing the standard process of ice content variation in soils during freezing, takes into account factors such as Poisson's ratio, groundwater table depth, and cooling rate. the calculation methods for the frost heave coefficients in the horizontal and vertical directions during the freezing process of positive permafrost were derived. By comparing the frost heave coefficients of silt and silty clay, the results show good agreement with experimental data. In the case study, the vertical frost heave coefficient of silt ranges from −1.37×10−3 to −7.67×10−3, and the horizontal frost heave coefficient ranges from −0.81×10−3 to −4.85×10−3 within the temperature range of −0.2 °C to −3 °C and a depth range of 0.2 m to 1 m. The percentage difference ranges from 10.4% to 77.7%, indicating the necessity of considering the anisotropy resulting from segregated frost heave. The calculation methods for horizontal and vertical frost heave coefficients in this study can provide valuable information for frost damage of engineering in the cold region.
-
Keywords:
- permafrost /
- moisture migration /
- temperature /
- depth of water table /
- frost heaving coefficient
-
冻土是一种温度在0 °C以下,且土颗粒被冰所胶结的特殊岩土[1],是由气体、水、冰、土颗粒组成的各向异性的非均质四相复合体。在开放系统中土的冻胀是土体中的液态水、气态水以及迁移水分相变成冰体积膨胀的过程,土颗粒、冰以及未冻水在水、热、力等各项因素的作用下相互影响[2 − 3],在这个过程中水分迁移是影响土体冻胀量的重要因素[4 − 5]。冻胀系数是描述土体在冻结过程中膨胀程度的参数,是材料的基本属性。
冻土按照温度、含冰状况和水热输运特征,可分为已冻土、正冻土、融土等三个互相联系的子系统[6]。正冻土是指土体温度处于0 °C以下,且在冷却过程中有冰晶形成和生长(冻融界面移动)的冻土[7]。
20世纪30年代,Taber等[8]用苯代替水在试验中证明了水分迁移是土体冻胀的主要因素,否定了早期对冻胀机理的理解,开辟了冻胀研究的新方向。在此基础上,国内外学者对土体冻胀过程中的水分迁移量以及影响水分迁移的因素进行了研究,例如毛细冻胀理论[9]将冻胀过程中毛细吸力作为水分迁移的主要驱动力,对水分迁移引起的冻胀量作出了定量的解释;分凝冻胀理论[10]认为温度梯度是未冻水向已冻区迁移的主要原因;分凝势理论[11]将热状况达到稳定的正冻土外界水分迁移入流通量与温度梯度的比值定义为分凝势,认为水分迁移速率与温度梯度成正比。
基于不同的冻胀理论模型,考虑不同因素对土体冻胀的影响,郭富强等[12]通过分析不同地下水位冻胀观测数据提出地下水位降低0.5~1.0 m,冻胀率减小71%~83.8%;杨自友等[13]考虑冻土冻融循环及不均匀冻胀特性,研究了非静水压力作用下不均匀冻胀系数、体积冻胀率对应力场的影响。由于分凝冻胀中的水分迁移方向以及土体围压作用,土体冻胀变形为各向异性,蔡海兵等[14]基于冻土正交各向异性变形特征,得到了考虑各向异性的热应变分量,建立了地层冻胀的弹塑性热力耦合数学模型,并验证了考虑冻胀各向异性的必要性;王贺等[15]仅考虑原位冻胀产生的冻胀变形,通过弹性力学方法引入土的泊松比,得到了正冻土的冻胀系数与冻胀率之间的关系。
开放系统中正冻土冻结过程由于地下水位深度、水分迁移强度、冻胀的各向异性等影响因素,目前仍无法对冻胀量进行准确的定量表示,由于土体参数的复杂性,这些理论中涉及的土体参数较难确定,且计算过程复杂,在工程实际中应用较为局限。大部分理论模型在计算冻胀量时仅考虑分凝冻胀,影响了计算结果的准确性。为解决以上问题,本文基于《冻土地区建筑地基基础设计规范》(JGJ 118—2011)[16]引入地下水位深度、降温速率等影响因素,通过对开放系统中饱和正冻土的原位冻胀量和分凝冻胀量累加得到较为准确的冻胀量计算公式,从而推导出开放系统中饱和正冻土冻胀系数与冻胀率、地下水位深度的关系,便于工程应用。
1. 开放系统中饱和正冻土冻胀计算方法
在开放系统中,土的冻胀根据水分来源可以分为原位冻胀和分凝冻胀,未冻水在温度梯度作用下会产生可观的水分迁移,从而形成分凝冻胀[17 − 19],水分迁移过程中水由高势能处向低势能处转移,水势差越大,水分迁移的速率越大[20 − 21]。天然条件下由于外界水分补给及水分迁移的作用,分凝冻胀是构成土体冻胀的主要分量,很多学者在研究冻胀时通常会忽略原位冻胀引起的体积变化。研究表明在单向冻结模式下,有外界水源补给的冻胀率是无水源补给的3~15倍,原位冻胀量占总冻胀量的6.25%~25%[22],是一个区间,因此在某些情况下原位冻胀并不能被忽略,为了计算更准确的冻胀量,应同时考虑原位水分冻结及迁移水分冻结引起的体积变化。假设土体为各向同性的线弹性体,图1为土体发生原位冻胀和分凝冻胀的微元体示意图。
土体在开放系统中z方向的变形量
,其中 为原位冻胀引起的变形量, 为外界水分补给后相变引起的变形量。即土体的冻胀率等于土体的原位冻胀率和分凝冻胀率的总和:(1) 式中:
——土体的总冻胀率/%; ——土体原2位冻胀产生的冻胀率/%; ——分凝冻胀产生的冻胀率/%。根据冻胀系数的定义
, 为温度变化量,引入泊松比,得到考虑围压的原位冻胀率与冻胀系数的关系。分凝冻胀量的计算是建立在外界水源流入后全部相变为冰的假定上,考虑外界水分渗透强度和地下水深度得到单位时间的外界水源的入流通量,对入流通量进行累积得到分凝冻胀引起的土体冻胀变化量。将原位冻胀量和分凝冻胀量累加得到土体总冻胀量,根据定义得到冻胀系数与冻胀率的关系式。选取文献中的试验进行同参数计算,将计算值与试验值进行对比,验证关系式的可靠性。2. 开放系统中饱和正冻土冻胀系数与冻胀率间的关系
冻土冻胀系数
,量纲为 ,是衡量冻土状态及性质的主要参数之一。冻胀系数是指土体在含水率不变的情况下,冻结时体积膨胀的比例,是土体体积增大率与温度下降率之间的比值,与温度变化率呈正比。冻胀率为土体冻胀总量与冻结深度(不包括冻胀量)的比值[23],是土体冻结前后体积之差与冻结前土体体积之比,以百分数表示。一般以土体试样冻结前后的高度差(∆z)与冻结前试样高度(h)之比表示,按照式(2)计算[16]:
(2) 式中:
——土的平均冻胀率/%; ——地表冻胀量/mm; ——冻层厚度/mm。在封闭系统中,即没有地下水补给的条件下,
与含水率( )间的关系可以按下式计算[16]:(3) 式中:
——土的干密度,取1.5 g/cm3; ——水的密度,取1.0 g/cm3; ——塑限含水率/%。按照实际情况,土的微元体在xz、yz平面受到围压约束,只在z方向发生变形,根据弹性力学中按位移求解温度应力的平面应变问题,推导出
与 之间的关系[15]:(4) 式中:
——初始冻结温度/°C; ——完全冻结温度/°C;α——冻胀系数;
——泊松比。相对含冰率指冻土中冰与全部水的质量比,对不同温度(
)下的相对含冰率进行线性插值得到不同温度下原位冻胀系数与冻胀率的关系:(5) 式中:
——温度Ti下的原位冻胀系数; —— 下的相对含冰率。表1为不同类型正冻土在不同温度下的相对含冰量和温度修正系数值[16],在土体冻结过程中可以通过相对含冰率判断正冻土的冻结状态,根据此表中的温度修正系数可以计算土体冻结过程中不同温度下的相对含冰率:
表 1 不同土质不同温度下的温度修正系数Table 1. The temperature correction coefficient under different soil types and temperatures土质 塑性指标 温度修订系数 −0.2 °C −0.5 °C −1.0 °C −2.0 °C −3.0 °C −5.0 °C −10 °C 砂土 — 0.35 0.22 0.15 0.08 0.07 0.05 0.02 粉土 Ip≤10 0.70 0.50 0.30 0.20 0.15 0.15 0.10 粉质黏土 10<Ip≤13 0.90 0.65 0.50 0.40 0.35 0.30 0.25 13<Ip≤17 1.00 0.80 0.70 0.60 0.50 0.45 0.40 黏土 17<Ip 1.10 0.90 0.80 0.70 0.60 0.55 0.50 注:Ip为塑性指数;—表示不取值。 (6) 式中:
——不同温度对应的温度修正系数,取值见表1。表2为《冻土地区建筑地基基础设计规范》(JGJ 118—2011)[16]中的季节冻土冻胀性分类表,根据此表可以计算不同含水率对应的冻胀率,例如可以采用线性插值法计算得到:
表 2 季节冻土冻胀性分类Table 2. Frost heaving classification of seasonal frozen soil土质 hw/m 冻胀等级 冻胀类别 黏性土 w0≤wp+2 >2.0 ≤1.0Ⅰ 不冻胀 ≤2.0 1.0< ≤3.5Ⅱ 弱冻胀 wp+2<w0≤wp+5 >2.0 ≤2.0 3.5< ≤6.0Ⅲ 冻胀 wp+5<w0≤wp+9 >2.0 ≤2.0 6.0< ≤12.0Ⅳ 强冻胀 wp+9< w0≤wp+15 >2.0 ≤2.0 >12.0Ⅴ 特强冻胀 粉土 w0≤19 >1.5 ≤1.0Ⅰ 不冻胀 ≤1.5 1.0< ≤3.5Ⅱ 弱冻胀 19< w0≤22 >1.5 1.0< ≤3.5Ⅱ 弱冻胀 ≤1.5 3.5< ≤6.0Ⅲ 冻胀 22< w0≤26 >1.5 ≤1.5 6.0< ≤12.0Ⅳ 强冻胀 26< w0≤30 >1.5 ≤1.5 >12.0Ⅴ 特强冻胀 注:w0为冻前天然含水率。 (7) 式中:
——含水率上下限; ——冻胀率上下限。分凝冻胀量为水分迁移产生的聚冰作用引起的体积增大:
(8) (9) 式中:
——分凝冻胀系数; ——弹性微元体在温度变化为 时,z方向产生的形变分量; ——入流通量/(m3∙s−1); ——水分迁移时间/s。(10) 式中:
——孔隙率。徐学祖等[6]通过统计分析方法得到入流通量计算公式:
(11) 式中:
——冻前地下水位距设计冻深的最小距离/m; ——设计冻深/m; ——冻结温度Tu~Tf的历时/s; ——渗透强度,定义为渗透强度与压差的 乘积,可以取最大值 ,P为外荷载/KPa。qt可以不考虑外荷载作用、时间因素,hf(t)取设计冻深。
入流通量:
(12) 入流量:
(13) 联立式(7)(8)(11)(12)得:
(14) 根据冻胀系数的定义:
(15) 将式(14)代入式(13)得:
(16) 式中:
, 为降温速度/(°C·h−1)。联立式(3)—(5)、(14)—(15),得到不同温度下原位冻胀及水分迁移作用下x、y、z三个方向冻胀系数与冻胀量之间的关系:
(17) 3. 算例分析
根据文献[24]中开放系统中庄河灌区土样单向冻结试验1,冻结在低温冻融循环机中进行,土样中大部分为细粒土,属于低液限粉土,采取底部补水模拟地下水补给,土样下部铺设透水石。试验2[12]为室外试验,土样为粉质黏土,采用不同高度的台座模拟不同地下水位补给。采用试验中的土样参数及边界条件计算理论冻胀系数,与试验值进行比较验证公式的可靠性,土体参数及边界条件见下表3。
表 3 土体参数及边界条件Table 3. Soil parameters and boundary conditions参数 试验1 试验2 土质 粉土 粉质黏土 降温速率/(°C∙h−1) −2 −0.07 初始含水率/% 22 23.45 泊松比 0.25 0.33 孔隙率 0.3 0.2 塑限含水率/% 21.7 18.5 孔隙率
,由式(10)可以得到饱和土在完全冻结下的原位冻胀量 。根据表1可知粉土在−10~0 °C发生冻胀,将温度修正系数代入式(6)得到完全冻结时的相对含冰率 。采用试验1的参数进行试算:当地下水位深度0.5 m,温度
=−0.2 °C时,相对含冰率 =1−0.7=0.3,在温度 =0 °C时,土体处于临界状态,相对含冰率 =0,冻胀系数:(18) 温度
=−0.5 °C时,相对含冰率 =1−0.5=0.5,冻胀系数:(19) 同理,可以通过试验1参数得到不同地下水位深度、不同温度下的z方向的冻胀系数(αz)结果,计算结果见表4。
表 4 不同地下水位深度不同温度下的Table 4. αz at different water table depths and temperatures温度
/°Chw=0.2 m hw=0.4 m hw=0.5 m hw=0.6 m hw=0.8 m hw=1 m 0 0 0 0 0 0 0 −0.2 −7.67 −6.26 −5.98 −5.79 −5.55 −5.41 −0.5 −6.05 −4.64 −4.36 −4.17 −3.93 −3.79 −1 −6.05 −4.64 −4.36 −4.17 −3.93 −3.79 −2 −4.44 −3.03 −2.75 −2.56 −2.32 −2.18 −3 −3.63 −2.22 −1.94 −1.75 −1.51 −1.37 −5 — — — — — — −10 −3.63 −2.22 −1.94 −1.75 −1.51 −1.37 根据试验1,地下水位深度hw=0 m,初始含水率
粉土在开放系统单向冻结试验中不同温度下的冻胀率 及冻胀系数 见表5。表 5 不同温度下z方向冻胀系数试验值Table 5. Experimental values of frost heaving coefficient in z direction at different temperatures温度/°C 冻胀率/% 冻胀系数/10−3 0 0 0 −0.5 0.2 −4 −1.8 1.6 −8.89 −3.5 2.2 −6.29 −5 2.9 −5.80 −7.5 3.4 −4.53 −9.4 3.6 −3.83 −11.5 3.8 −3.30 3.1 冻胀系数计算值与试验值对比
冻结试验1在低温冻融循环机中进行,降温速率为2 °C/h,初始含水率22%,采取底部补水模拟地下水补给,采用与表3中相同的土体试验参数计算得到理论冻胀系数值,对hw=0.2 m冻胀系数计算值与hw=0 m试验值对比,验证公式的有效性。对比值见图2,计算值与试验值随温度的变化趋势一致,地下水位hw=0.2 m计算值小于hw=0 m试验值,随着地下水位的降低,冻胀系数减小,说明地下水位距设计冻深的距离hw与冻胀系数呈反比,与实际相符。
冻胀系数计算值在−0.2~0 °C变化率较大,呈上升趋势,温度T=0.2 °C达到峰值,−3~−0.2 °C呈下降趋势,−10~−3 °C随着温度的下降冻胀系数值趋于平缓。在开放系统单向冻结试验中,水源距试样距离hw=0 m,冻胀系数试验值大于hw=0.2 m计算值,地下水位降低,渗透强度随之减小,冻胀系数减小。冻胀系数试验值在−1.8~0 °C变化率较大,呈上升趋势,温度T=−1.8 °C达到峰值,−3.5~−1.8 °C呈下降趋势, −11.5~−3.5 °C随着温度的下降冻胀系数值趋于平缓,试验值与计算值总体趋势相同。冻胀系数试验峰值滞后于计算值,这是因为试验采用土样大部分为细粒土,存在部分粗粒土,粗粒土冻胀速率相对较慢,并且在试验初期制冷机启动时间较短导致土体冻结速率较慢。
由于未冻水含量与温度变化有着密切的联系,正冻土在冻结过程中,不同温度下冻结强度不同,崔托维奇[25]将正冻土划分为剧烈相变区(下降1 °C含冰量变化≥1%)、过渡相变区(下降1 °C含冰量变化0.1%~1%)、已冻结区(下降1 °C含冰量变化≤1%)。
冻胀系数随温度变化的三个阶段符合崔托维奇按照温度划分的剧烈相变区、过渡相变区、已冻结区。
试验2在室外进行,采用0.5,1.0 m台座模拟0.5,1.0 m地下水位高度,土样为粉质黏土,降温速率通过日平均气温计算得到,为0.07 °C/d,具体计算参数见表3中试验2参数。通过公式(17)计算得到hw=0.5,1.0 m冻胀系数值与试验值对比见图3。
从图中可以看到在相同地下水位条件下,计算值与试验值变化趋势一致,hw=0.5 m冻胀系数值在−0.2 °C达到峰值
,试验值在−0.35 °C达到峰值 ,峰值相差 。hw=1.0 m冻胀系数值在−0.2 °C达到峰值 ,试验值在−0.35 °C达到峰值 ,峰值相差 。冻胀系数试验值整体低于计算值,这是因为计算值采用日平均气温,而室外日气温处于昼夜波动状态,冻胀相对较弱。冻胀系数变化趋势与试验1相同,符合剧烈相变区、过渡相变区、已冻结区的划分。随着地下水位高度的降低,冻胀系数值减小,hw=0.5 m与hw=1.0 m冻胀系数计算值差值达到72.7%,试验值差值为75.9%,进一步验证了不同地下水位条件公式的有效性。3.2 冻胀系数各向异性必要性
对比不同温度不同地下水位深度各方向通过试验1参数计算得到的冻胀系数值,定义冻胀系数差值百分比为
,来描述不同方向冻胀系数值的差异情况。不同温度不同地下水位深度各向异性冻胀系数计算对比值见表6。表 6 不同温度不同地下水位深度各向异性冻胀系数比较Table 6. Comparison of anisotropic frost heave coefficients at different temperatures and different groundwater depthshw/m 冻胀系数 温度/°C −0.2 −0.5 −2 −3 0.2 αz/10−3 −7.67 −6.05 −4.44 −3.63 αx,y/10−3 −4.85 −3.23 −1.62 −0.81 差值百分比 36.8% 46.6% 63.5% 77.7% 0.4 αz/10−3 −6.26 −4.64 −3.03 −2.22 αx,y/10−3 −4.85 −3.23 −1.62 −0.81 差值百分比 22.5% 30.4% 46.5% 63.5% 0.5 αz/10−3 −5.98 −4.36 −2.75 −1.94 αx,y/10−3 −4.85 −3.23 −1.62 −0.81 差值百分比 18.9% 25.9% 41.1% 58.2% 0.6 αz/10−3 −5.79 −4.17 −2.56 −1.75 αx,y/10−3 −4.85 −3.23 −1.62 −0.81 差值百分比 16.2% 22.5% 36.7% 53.7% 0.8 αz/10−3 −5.55 −3.93 −2.32 −1.51 αx,y/10−3 −4.85 −3.23 −1.62 −0.81 差值百分比 12.6% 17.8% 30.2% 46.4% 1.0 αz/10−3 −5.41 −3.79 −2.18 −1.37 αx,y/10−3 −4.85 −3.23 −1.62 −0.81 差值百分比 10.4% 14.8% 25.7% 40.9% 分析不同温度不同地下水位深度的冻胀系数计算值,得到冻胀系数各向异性差值百分比,给定温度条件,随地下水位距设计冻深距离hw的增大,冻胀系数差值百分比减小,在−3~−0.2 °C,0.2~1 m范围内冻胀系数差值百分比为10.4%~77.7%,这说明了冻胀系数各向异性的必要性,不同方向的差值并不能被忽略;温度为−0.2 °C时,冻胀系数差值百分比为10.4%,这说明在剧烈相变区,地下水位较深,水分补给较弱时,αx,y不能被忽略。
图4为不同温度下冻胀系数差值百分比随地下水位深度变化图,给定地下水位深度,随温度的降低,冻胀系数差值百分比增大,这是因为在剧烈相变区孔隙内原有水分相变成冰,原位冻胀作用明显,冻胀系数差值百分比较小,αx,y被忽略后对结果准确性影响较大;随着温度的降低,原位冻胀作用减弱,外界水分依然在温度梯度的作用下向上迁移相变成冰不同方向冻胀系数差值百分比增大,这与Konrad[11]提出的当正冻土的热状况达到稳定时,冻结缘附近水分迁移入流通量与温度梯度呈一个比值相符合。
3.3 冻胀系数与地下水位深度的关系
水源补给是冻结时水分迁移及由此所产生冻胀的持续性与强弱的直接原因,而地下水位深度是水源补给强度的重要因素。大量试验数据和文献表明[25 − 28],在相同土体参数、边界条件下,地下水位深度与冻胀率呈负指数关系:
(20) 式中:a、b——与当地气候、土质等因素有关的经验系数。
根据本文推导的冻胀系数计算公式采用试验1参数分别计算不同温度下地下水位深度hw=0.2~1.0 m、z方向冻胀系数
,hw步距为0.2,不同地下水位冻胀系数对比见图5。随着地下水位深度的增加,冻胀系数值减小,由于室内冻胀试验水位深度为0,试验值大于计算值。这是由于当地下水埋深较浅时,土体在冻胀变形过程中水分补给量较大,整体冻胀较为剧烈。
对比冻胀趋于稳定时相邻地下水位深度的冻胀系数值,0.2~0.4 m冻胀系数值减小38.8%;0.4~0.6 m冻胀系数值减小21.2%;0.6~0.8 m冻胀系数值减小13.7%;0.8~1.0 m冻胀系数值减小9.3%。
冻胀系数计算值
随地下水位深度变化见图6,随着地下水位的降低,冻胀系数 的变化并不是线性的,随着地下水位的降低,冻胀系数 减小且变化率呈降低趋势,这表明地下水位降低到一定高度后,水分补给变化减小,趋于稳定,地下水位深度对冻胀系数的影响削弱,冻胀系数 与地下水位呈负指数关系,这与式(20)得到的结论相同,进一步印证了本文提出公式的准确性。通过线性插值得到0,−1,−2,−3,−5,−10 °C冻胀系数试验值,与hw=1 m的冻胀系数计算值进行差值比较。土的冻胀稳定后,地下水位高度降低1 m冻胀系数减小62.8%~80.6%,不同温度下冻胀系数减小百分率见表7,与郭富强等[8]通过试验得到当地下水位高度降低0.5~1.0 m后,冻胀率将减小71%~83.8%的结论基本吻合。
表 7 不同温度下地下水位深度对冻胀系数影响Table 7. Influence of groundwater depth on frost heave coefficient at different temperatures with hw = 1 m温度/°C 冻胀系数减小百分率/% 0 0 0 0 −1 −5.88 −3.79 35.5 −2 −8.58 −2.18 74.6 −3 −7.05 −1.37 80.6 −5 −5.80 — — −10 −3.68 −1.37 62.8 4. 结论
(1)通过试验对照,本文冻胀系数理论计算结果与试验值在不同温度下变化趋势相同,在计算值中也可以体现土体冻胀剧烈相变区、过渡相变区、已冻结区,不同地下水位深度对冻胀系数的影响与试验值基本一致,验证了关系式的准确性。
(2)通过不同方向冻胀系数差值百分比对比,不同方向冻胀系数的差值不能被忽略,冻胀系数各向异性具有必要性。不同方向冻胀系数的差值百分比与地下水位距设计冻深的最小距离hw呈反比,越靠近剧烈相变区、地下水位距离设计冻深越远,不同方向冻胀系数差值越小,原位冻胀在冻胀过程中的占比就越大,产生的冻胀作用明显,
不能被忽略。(3)z方向冻胀系数
与地下水位距设计冻深的最小距离hw呈负指数关系,土体冻胀稳定后,地下水位降低1m,z方向冻胀系数 减小62.8%~80.6%。在开放系统中地下水深度是影响土体冻胀率、冻胀系数的重要因素。由于土体冻胀是多物理场耦合的结果,且土体冻胀影响因素具有复杂性,进一步的研究还需要增加补水条件中的外部压力条件和外界降水条件,考虑在水分迁移过程中水热力相互耦合的影响。
-
表 1 不同土质不同温度下的温度修正系数
Table 1 The temperature correction coefficient under different soil types and temperatures
土质 塑性指标 温度修订系数 −0.2 °C −0.5 °C −1.0 °C −2.0 °C −3.0 °C −5.0 °C −10 °C 砂土 — 0.35 0.22 0.15 0.08 0.07 0.05 0.02 粉土 Ip≤10 0.70 0.50 0.30 0.20 0.15 0.15 0.10 粉质黏土 10<Ip≤13 0.90 0.65 0.50 0.40 0.35 0.30 0.25 13<Ip≤17 1.00 0.80 0.70 0.60 0.50 0.45 0.40 黏土 17<Ip 1.10 0.90 0.80 0.70 0.60 0.55 0.50 注:Ip为塑性指数;—表示不取值。 表 2 季节冻土冻胀性分类
Table 2 Frost heaving classification of seasonal frozen soil
土质 hw/m 冻胀等级 冻胀类别 黏性土 w0≤wp+2 >2.0 ≤1.0Ⅰ 不冻胀 ≤2.0 1.0< ≤3.5Ⅱ 弱冻胀 wp+2<w0≤wp+5 >2.0 ≤2.0 3.5< ≤6.0Ⅲ 冻胀 wp+5<w0≤wp+9 >2.0 ≤2.0 6.0< ≤12.0Ⅳ 强冻胀 wp+9< w0≤wp+15 >2.0 ≤2.0 >12.0Ⅴ 特强冻胀 粉土 w0≤19 >1.5 ≤1.0Ⅰ 不冻胀 ≤1.5 1.0< ≤3.5Ⅱ 弱冻胀 19< w0≤22 >1.5 1.0< ≤3.5Ⅱ 弱冻胀 ≤1.5 3.5< ≤6.0Ⅲ 冻胀 22< w0≤26 >1.5 ≤1.5 6.0< ≤12.0Ⅳ 强冻胀 26< w0≤30 >1.5 ≤1.5 >12.0Ⅴ 特强冻胀 注:w0为冻前天然含水率。 表 3 土体参数及边界条件
Table 3 Soil parameters and boundary conditions
参数 试验1 试验2 土质 粉土 粉质黏土 降温速率/(°C∙h−1) −2 −0.07 初始含水率/% 22 23.45 泊松比 0.25 0.33 孔隙率 0.3 0.2 塑限含水率/% 21.7 18.5 表 4 不同地下水位深度不同温度下的
Table 4 αz at different water table depths and temperatures
温度
/°Chw=0.2 m hw=0.4 m hw=0.5 m hw=0.6 m hw=0.8 m hw=1 m 0 0 0 0 0 0 0 −0.2 −7.67 −6.26 −5.98 −5.79 −5.55 −5.41 −0.5 −6.05 −4.64 −4.36 −4.17 −3.93 −3.79 −1 −6.05 −4.64 −4.36 −4.17 −3.93 −3.79 −2 −4.44 −3.03 −2.75 −2.56 −2.32 −2.18 −3 −3.63 −2.22 −1.94 −1.75 −1.51 −1.37 −5 — — — — — — −10 −3.63 −2.22 −1.94 −1.75 −1.51 −1.37 表 5 不同温度下z方向冻胀系数试验值
Table 5 Experimental values of frost heaving coefficient in z direction at different temperatures
温度/°C 冻胀率/% 冻胀系数/10−3 0 0 0 −0.5 0.2 −4 −1.8 1.6 −8.89 −3.5 2.2 −6.29 −5 2.9 −5.80 −7.5 3.4 −4.53 −9.4 3.6 −3.83 −11.5 3.8 −3.30 表 6 不同温度不同地下水位深度各向异性冻胀系数比较
Table 6 Comparison of anisotropic frost heave coefficients at different temperatures and different groundwater depths
hw/m 冻胀系数 温度/°C −0.2 −0.5 −2 −3 0.2 αz/10−3 −7.67 −6.05 −4.44 −3.63 αx,y/10−3 −4.85 −3.23 −1.62 −0.81 差值百分比 36.8% 46.6% 63.5% 77.7% 0.4 αz/10−3 −6.26 −4.64 −3.03 −2.22 αx,y/10−3 −4.85 −3.23 −1.62 −0.81 差值百分比 22.5% 30.4% 46.5% 63.5% 0.5 αz/10−3 −5.98 −4.36 −2.75 −1.94 αx,y/10−3 −4.85 −3.23 −1.62 −0.81 差值百分比 18.9% 25.9% 41.1% 58.2% 0.6 αz/10−3 −5.79 −4.17 −2.56 −1.75 αx,y/10−3 −4.85 −3.23 −1.62 −0.81 差值百分比 16.2% 22.5% 36.7% 53.7% 0.8 αz/10−3 −5.55 −3.93 −2.32 −1.51 αx,y/10−3 −4.85 −3.23 −1.62 −0.81 差值百分比 12.6% 17.8% 30.2% 46.4% 1.0 αz/10−3 −5.41 −3.79 −2.18 −1.37 αx,y/10−3 −4.85 −3.23 −1.62 −0.81 差值百分比 10.4% 14.8% 25.7% 40.9% 表 7 不同温度下地下水位深度对冻胀系数影响
Table 7 Influence of groundwater depth on frost heave coefficient at different temperatures with hw = 1 m
温度/°C 冻胀系数减小百分率/% 0 0 0 0 −1 −5.88 −3.79 35.5 −2 −8.58 −2.18 74.6 −3 −7.05 −1.37 80.6 −5 −5.80 — — −10 −3.68 −1.37 62.8 -
[1] 马巍,王大雁. 冻土力学[M]. 北京:科学出版社,2014. [MA Wei,WANG Dayan. Mechanics of frozen soil[M]. Beijing:Science Press,2014. (in Chinese)] MA Wei, WANG Dayan. Mechanics of frozen soil[M]. Beijing: Science Press, 2014. (in Chinese)
[2] 陈新瑞,宋玲,孙雯,等. “开敞系统”下单向冻融试验装置的研制与应用[J]. 水文地质工程地质,2020,47(3):69 − 78. [Chen Xinrui,Song Ling,Sun Wen,et al. Development and application of unidirectional freeze-thawing test device under the “open-system conditions”[J]. Hydrogeology & Engineering Geology,2020,47(3):69 − 78. (in Chinese with English abstract)] Chen Xinrui, Song Ling, Sun Wen, et al. Development and application of unidirectional freeze-thawing test device under the “open-system conditions”[J]. Hydrogeology & Engineering Geology, 2020, 47(3): 69 − 78. (in Chinese with English abstract)
[3] 温智,邓友生,冯文杰,等. 冻土水分迁移机理研究:评述与展望[J]. 冰川冻土,2023,45(2):588 − 598. [WEN Zhi,DENG Yousheng,FENG Wenjie,et al. Study on the mechanism of moisture migration in freezing soils:Review and prospect[J]. Journal of Glaciology and Geocryology,2023,45(2):588 − 598. (in Chinese with English abstract)] WEN Zhi, DENG Yousheng, FENG Wenjie, et al. Study on the mechanism of moisture migration in freezing soils: Review and prospect[J]. Journal of Glaciology and Geocryology, 2023, 45(2): 588 − 598. (in Chinese with English abstract)
[4] 胡建华,汪稔,胡明鉴,等. 多年冻土路堤水平排水板的作用机理与试验研究[J]. 中南大学学报(自然科学版),2009,40(5):1451 − 1456. [HU Jianhua,WANG Ren,HU Mingjian,et al. Test and mechanism of horizontal plastic drain in permafrost regions embankment[J]. Journal of Central South University (Science and Technology),2009,40(5):1451 − 1456. (in Chinese with English abstract)] HU Jianhua, WANG Ren, HU Mingjian, et al. Test and mechanism of horizontal plastic drain in permafrost regions embankment[J]. Journal of Central South University (Science and Technology), 2009, 40(5): 1451 − 1456. (in Chinese with English abstract)
[5] 程国栋,金会军. 青藏高原多年冻土区地下水及其变化[J]. 水文地质工程地质,2013,40(1):1 − 11. [CHENG Guodong,JIN Huijun. Groundwater in the permafrost regions on the Qinghai-Tibet Plateau and it changes[J]. Hydrogeology & Engineering Geology,2013,40(1):1 − 11. (in Chinese with English abstract)] CHENG Guodong, JIN Huijun. Groundwater in the permafrost regions on the Qinghai-Tibet Plateau and it changes[J]. Hydrogeology & Engineering Geology, 2013, 40(1): 1 − 11. (in Chinese with English abstract)
[6] 徐学祖,王家澄,张立新. 冻土物理学[M]. 北京:科学出版社,2001. [XU Xuezu,WANG Jiacheng,ZHANG Lixin. Permafrost physics [M]. Beijing:Science Press,2001. (in Chinese)] XU Xuezu, WANG Jiacheng, ZHANG Lixin. Permafrost physics [M]. Beijing: Science Press, 2001. (in Chinese)
[7] 肖东辉,马巍,张泽,等. 冻融与静荷载双重作用下土体中水分场变化规律模型试验研究[J]. 水文地质工程地质,2016,43(5):155 − 160. [XIAO Donghui,MA Wei,ZHANG Ze,et al. Model test on moisture’s change law in the soil subjected to freezing-thawing and static load actions[J]. Hydrogeology & Engineering Geology,2016,43(5):155 − 160. (in Chinese with English abstract)] XIAO Donghui, MA Wei, ZHANG Ze, et al. Model test on moisture’s change law in the soil subjected to freezing-thawing and static load actions[J]. Hydrogeology & Engineering Geology, 2016, 43(5): 155 − 160. (in Chinese with English abstract)
[8] TABER S. The mechanics of frost heaving[J]. the Journal of Geology,1930,38(4):303 − 317. DOI: 10.1086/623720
[9] EVERETT D H. The thermodynamics of frost damage to porous solids[J]. Transactions of the Faraday Society,1961,57:1541 − 1551. DOI: 10.1039/tf9615701541
[10] TAKAGI S. Segregation freezing as the cause of suction force for ice lens formation[J]. Engineering Geology,1979,13(1/2/3/4):93 − 100.
[11] KONRAD J M,MORGENSTERN N R. The segregation potential of a freezing soil[J]. Canadian Geotechnical Journal,1981,18(4):482 − 491. DOI: 10.1139/t81-059
[12] 郭富强,史海滨,程满金,等. 不同地下水位下渠基冻胀规律与保温板适宜厚度确定[J]. 农业工程学报,2018,34(19):95 − 103. [GUO Fuqiang,SHI Haibin,CHENG Manjin,et al. Law of frost heave of canal foundation and appropriate thickness of insulation board under different groundwater levels[J]. Transactions of the Chinese Society of Agricultural Engineering,2018,34(19):95 − 103. (in Chinese with English abstract)] DOI: 10.11975/j.issn.1002-6819.2018.19.013 GUO Fuqiang, SHI Haibin, CHENG Manjin, et al. Law of frost heave of canal foundation and appropriate thickness of insulation board under different groundwater levels[J]. Transactions of the Chinese Society of Agricultural Engineering, 2018, 34(19): 95 − 103. (in Chinese with English abstract) DOI: 10.11975/j.issn.1002-6819.2018.19.013
[13] 杨自友,高鹏,尚阳光,等. 考虑冻融循环效应的非静水压寒区隧道应力场弹性分析[J]. 隧道建设(中英文),2023,43(增刊2):48 − 58. [YANG Ziyou,GAO Peng,SHANG Yangguang,et al. Elastic analysis of tunnel stress field in cold zone under non−hydrostatic pressure considering freeze−thaw cycle effect[J]. Tunnel Construction,2023,43(Sup2):48 − 58. (in Chinese with English abstract)] YANG Ziyou, GAO Peng, SHANG Yangguang, et al. Elastic analysis of tunnel stress field in cold zone under non−hydrostatic pressure considering freeze−thaw cycle effect[J]. Tunnel Construction, 2023, 43(Sup2): 48 − 58. (in Chinese with English abstract)
[14] 蔡海兵,程桦,姚直书,等. 基于冻土正交各向异性冻胀变形的隧道冻结期地层位移数值分析[J]. 岩石力学与工程学报,2015,34(8):1667 − 1676. [CAI Haibing,CHENG Hua,YAO Zhishu,et al. Numerical analysis of ground displacement due to orthotropic frost heave of frozen soil in freezing period of tunnel[J]. Chinese Journal of Rock Mechanics and Engineering,2015,34(8):1667 − 1676. (in Chinese with English abstract)] CAI Haibing, CHENG Hua, YAO Zhishu, et al. Numerical analysis of ground displacement due to orthotropic frost heave of frozen soil in freezing period of tunnel[J]. Chinese Journal of Rock Mechanics and Engineering, 2015, 34(8): 1667 − 1676. (in Chinese with English abstract)
[15] 王贺,郭春香,吴亚平,等. 基于弹性力学考虑冰水相变过程下多年冻土冻胀系数与冻胀率之间的关系[J]. 岩石力学与工程学报,2018,37(12):2839 − 2845. [WANG He,GUO Chunxiang,WU Yaping,et al. Relationship between the frost−heaving coefficient and the frost−heaving rate of permafrost soils considering the ice water phase transformation based on elastic mechanics[J]. Chinese Journal of Rock Mechanics and Engineering,2018,37(12):2839 − 2845. (in Chinese with English abstract)] WANG He, GUO Chunxiang, WU Yaping, et al. Relationship between the frost−heaving coefficient and the frost−heaving rate of permafrost soils considering the ice water phase transformation based on elastic mechanics[J]. Chinese Journal of Rock Mechanics and Engineering, 2018, 37(12): 2839 − 2845. (in Chinese with English abstract)
[16] 中华人民共和国住房和城乡建设部. 冻土地区建筑地基基础设计规范:JGJ 118—2011[S]. 北京:中国建筑工业出版社,2012. [Ministry of Housing and Urban−Rural Development of the People's Republic of China. Code for design of soil foundation of building in frozen soil region:JGJ 118—2011[S]. Beijing:China Architecture & Building Press,2012. (in Chinese)] Ministry of Housing and Urban−Rural Development of the People's Republic of China. Code for design of soil foundation of building in frozen soil region: JGJ 118—2011[S]. Beijing: China Architecture & Building Press, 2012. (in Chinese)
[17] 周家作,韦昌富,李东庆,等. 饱和粉土冻胀过程试验研究及数值模拟[J]. 岩石力学与工程学,2017,36(2):485 − 495. [ZHOU Jiazuo,WEI Changfu,LI Dongqing,et al. Experimental study and numerical simulation to the process of frost heave in saturated silt[J]. Chinese Journal of Rock Mechanics and Engineering,2017,36(2):485 − 495. (in Chinese with English abstract)] ZHOU Jiazuo, WEI Changfu, LI Dongqing, et al. Experimental study and numerical simulation to the process of frost heave in saturated silt[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(2): 485 − 495. (in Chinese with English abstract)
[18] 刘泉声,黄诗冰,康永水,等. 低温饱和岩石未冻水含量与冻胀变形模型研究[J]. 岩石力学与工程学报,2016,35(10):2000 − 2012. [LIU Quansheng,HUANG Shibing,KANG Yongshui,et al. Study of unfrozen water content and frost heave model for saturated rock under low temperature[J]. Chinese Journal of Rock Mechanics and Engineering,2016,35(10):2000 − 2012. (in Chinese with English abstract)] LIU Quansheng, HUANG Shibing, KANG Yongshui, et al. Study of unfrozen water content and frost heave model for saturated rock under low temperature[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(10): 2000 − 2012. (in Chinese with English abstract)
[19] 周洁,李泽垚,唐益群,等. 软黏土三维水−热−力耦合冻融模型及应用[J]. 哈尔滨工业大学学报,2022,54(2):117 − 125. [ZHOU Jie,LI Zeyao,TANG Yiqun,et al. Three−dimensional thermal−hydraulic−mechanical coupling frost thaw model of soft clay and its application[J]. Journal of Harbin Institute of Technology,2022,54(2):117 − 125. (in Chinese with English abstract)] DOI: 10.11918/202011033 ZHOU Jie, LI Zeyao, TANG Yiqun, et al. Three−dimensional thermal−hydraulic−mechanical coupling frost thaw model of soft clay and its application[J]. Journal of Harbin Institute of Technology, 2022, 54(2): 117 − 125. (in Chinese with English abstract) DOI: 10.11918/202011033
[20] 于钱米,邰博文,牛吉强,等. 细粒土空间不均匀分布对水分迁移的影响[J]. 中南大学学报(自然科学版),2020,51(12):3503 − 3514. [YU Qianmi,TAI Bowen,NIU Jiqiang,et al. The effect of unevenly distributed fine−grained soil in space on moisture migration[J]. Journal of Central South University(Science and Technology),2020,51(12):3503 − 3514. (in Chinese with English abstract)] DOI: 10.11817/j.issn.1672-7207.2020.12.024 YU Qianmi, TAI Bowen, NIU Jiqiang, et al. The effect of unevenly distributed fine−grained soil in space on moisture migration[J]. Journal of Central South University(Science and Technology), 2020, 51(12): 3503 − 3514. (in Chinese with English abstract) DOI: 10.11817/j.issn.1672-7207.2020.12.024
[21] 曹文炳,万力,陈康,等. 季节冻土对包气带水分迁移的影响[J]. 水文地质工程地质,2014,41(6):1 − 5. [CAO Wenbing,WAN Li,CHEN Kang,et al. Influence of seasonal frozen soil on moisture migration in the aerated zone[J]. Hydrogeology & Engineering Geology,2014,41(6):1 − 5. (in Chinese with English abstract)] CAO Wenbing, WAN Li, CHEN Kang, et al. Influence of seasonal frozen soil on moisture migration in the aerated zone[J]. Hydrogeology & Engineering Geology, 2014, 41(6): 1 − 5. (in Chinese with English abstract)
[22] 王天亮,王海航,宋宏芳,等. 人工冻结粉质黏土力学性能演化规律研究[J]. 中国铁道科学,2019,40(1):1 − 7. [WANG Tianliang,WANG Haihang,SONG Hongfang,et al. Evolution laws of mechanical properties of artificially frozen silty clay[J]. China Railway Science,2019,40(1):1 − 7. (in Chinese with English abstract)] DOI: 10.3969/j.issn.1001-4632.2019.01.01 WANG Tianliang, WANG Haihang, SONG Hongfang, et al. Evolution laws of mechanical properties of artificially frozen silty clay[J]. China Railway Science, 2019, 40(1): 1 − 7. (in Chinese with English abstract) DOI: 10.3969/j.issn.1001-4632.2019.01.01
[23] 唐益群,洪军,杨坪,等. 人工冻结作用下淤泥质黏土冻胀特性试验研究[J]. 岩土工程学报,2009,31(5):772 − 776. [TANG Yiqun,HONG Jun,YANG Ping,et al. Frost−heaving behaviors of mucky clay by artificial horizontal freezing method[J]. Chinese Journal of Geotechnical Engineering,2009,31(5):772 − 776. (in Chinese with English abstract)] DOI: 10.3321/j.issn:1000-4548.2009.05.021 TANG Yiqun, HONG Jun, YANG Ping, et al. Frost−heaving behaviors of mucky clay by artificial horizontal freezing method[J]. Chinese Journal of Geotechnical Engineering, 2009, 31(5): 772 − 776. (in Chinese with English abstract) DOI: 10.3321/j.issn:1000-4548.2009.05.021
[24] 孙忠臣. 庄河灌区渠基土冻胀试验研究及保温板防冻厚度数值模拟分析[D]. 沈阳:沈阳农业大学,2020. [SUN Zhongchen. Experimental study on frost heave of canal subsoil in Zhuanghe Irrigation Area and numerical simulation analysis of anti−frost thickness of insulation board[D]. Shenyang:Shenyang Agricultural University,2020. (in Chinese with English abstract)] SUN Zhongchen. Experimental study on frost heave of canal subsoil in Zhuanghe Irrigation Area and numerical simulation analysis of anti−frost thickness of insulation board[D]. Shenyang: Shenyang Agricultural University, 2020. (in Chinese with English abstract)
[25] 崔托维奇. 冻土力学[M]. 张长庆,朱元林,译. 北京:科学出版社,1985. [TSYTOVICH H A. The mechanics of frozen ground[M]. ZHANG Changqing,ZHU Yuanlin,trans. Beijing:Science Press,1985. (in Chinese)] TSYTOVICH H A. The mechanics of frozen ground[M]. ZHANG Changqing, ZHU Yuanlin, trans. Beijing: Science Press, 1985. (in Chinese)
[26] 陈肖柏,刘建坤,刘鸿绪,等. 土的冻结作用与地基[M]. 北京:科学出版社,2006. [CHEN Xiaobai,LIU Jiankun,LIU Hongxu,et al. Frost action of soil and foundation engineering[M]. Beijing:Science Press,2006. (in Chinese)] CHEN Xiaobai, LIU Jiankun, LIU Hongxu, et al. Frost action of soil and foundation engineering[M]. Beijing: Science Press, 2006. (in Chinese)
[27] 肖旻,王正中,刘铨鸿,等. 考虑冻土与结构相互作用的梯形渠道冻胀破坏弹性地基梁模型[J]. 水利学报,2017,48(10):1229 − 1239. [XIAO Min,WANG Zhengzhong,LIU Quanhong,et al. Elastic foundation beam model of frost heave damage of trapezoidal canal considering interaction between frozen soil and lining stucture[J]. Journal of Hydraulic Engineering,2017,48(10):1229 − 1239. (in Chinese with English abstract)] XIAO Min, WANG Zhengzhong, LIU Quanhong, et al. Elastic foundation beam model of frost heave damage of trapezoidal canal considering interaction between frozen soil and lining stucture[J]. Journal of Hydraulic Engineering, 2017, 48(10): 1229 − 1239. (in Chinese with English abstract)
[28] 何鹏飞,候光亮,董建华,等. 梯形渠道衬砌冻胀破坏弹性地基板模型[J]. 农业工程学报,2022,38(23):91 − 100. [HE Pengfei,HOU Guangliang,DONG Jianhua,et al. Elastic foundation plate model for the frost heave damage of trapezoidal canal lining[J]. Transactions of the Chinese Society of Agricultural Engineering,2022,38(23):91 − 100. (in Chinese with English abstract)] DOI: 10.11975/j.issn.1002-6819.2022.23.010 HE Pengfei, HOU Guangliang, DONG Jianhua, et al. Elastic foundation plate model for the frost heave damage of trapezoidal canal lining[J]. Transactions of the Chinese Society of Agricultural Engineering, 2022, 38(23): 91 − 100. (in Chinese with English abstract) DOI: 10.11975/j.issn.1002-6819.2022.23.010
[29] 王正中,孙涛,郑艾磊,等. 水−热−力耦合作用下寒区弧底梯形渠道结构优化设计[J]. 中南大学学报(自然科学版),2023,54(5):1916 − 1929. [WANG Zhengzhong,SUN Tao,ZHENG Ailei,et al. Structural optimization design for trapezoidal canal with arc−bottom during thermal−moisture−mechanical coupling in cold regions[J]. Journal of Central South University (Science and Technology),2023,54(5):1916 − 1929. (in Chinese with English abstract)] DOI: 10.11817/j.issn.1672-7207.2023.05.025 WANG Zhengzhong, SUN Tao, ZHENG Ailei, et al. Structural optimization design for trapezoidal canal with arc−bottom during thermal−moisture−mechanical coupling in cold regions[J]. Journal of Central South University (Science and Technology), 2023, 54(5): 1916 − 1929. (in Chinese with English abstract) DOI: 10.11817/j.issn.1672-7207.2023.05.025