• 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

温度梯度作用下非饱和土中铅离子迁移机理研究

周凤玺, 罗玚, 马强

周凤玺, 罗玚, 马强. 温度梯度作用下非饱和土中铅离子迁移机理研究[J]. 岩土工程学报, 2023, 45(9): 1822-1830. DOI: 10.11779/CJGE20220741
引用本文: 周凤玺, 罗玚, 马强. 温度梯度作用下非饱和土中铅离子迁移机理研究[J]. 岩土工程学报, 2023, 45(9): 1822-1830. DOI: 10.11779/CJGE20220741
ZHOU Fengxi, LUO Yang, MA Qiang. Migration mechanism of lead ions in unsaturated loess under temperature gradient[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(9): 1822-1830. DOI: 10.11779/CJGE20220741
Citation: ZHOU Fengxi, LUO Yang, MA Qiang. Migration mechanism of lead ions in unsaturated loess under temperature gradient[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(9): 1822-1830. DOI: 10.11779/CJGE20220741

温度梯度作用下非饱和土中铅离子迁移机理研究  English Version

基金项目: 

国家自然科学基金项目 11962016

国家自然科学基金项目 51978320

甘肃省基础研究创新群体项目 20JR5RA478

甘肃省教育厅高校教师创新基金项目 2023A-164

详细信息
    作者简介:

    周凤玺(1979—),男,博士,教授,主要从事岩土工程方面的教学和科研工作。E-mail: geolut@163.com

  • 中图分类号: TU431

Migration mechanism of lead ions in unsaturated loess under temperature gradient

  • 摘要: 基于多孔介质理论及连续介质热弹性力学,研究了温度梯度作用下重金属离子在非饱和黄土中的迁移机理。首先,通过非等温吸附试验确定了黄土颗粒对铅离子的吸附模型;其次,以固、液、气三相系统中各组分质量、能量以及动量守恒方程为基础,建立了表征温度梯度下非饱和土中污染物迁移的多场耦合数学模型,模型考虑了重金属离子的吸附、孔隙流体的迁移及相变、温度梯度下土体的变形等因素,以温度、孔隙率、孔隙水压、孔隙气压、重金属离子浓度、吸附浓度、土体骨架位移为基本未知量,较全面反映了非饱和土体中热质迁移规律及变形特点;最后,开展非饱和黄土中铅离子迁移土柱试验,验证理论模型的可靠性,并分析了各因素对铅离子迁移规律的影响。结果表明:Langmuir非线性吸附模型比较符合黄土颗粒对铅离子的吸附特征,最大吸附量随温度升高而降低,温度梯度对非饱和黄土中污染物迁移有较明显促进作用;不计土体湿陷效应时,温度梯度及吸湿膨胀使土柱产生拉伸变形,重力、孔隙水压及孔隙气压共同作用使土柱产生压缩变形,后者影响大于前者,故土柱整体为压缩变形;污染物迁移随时间增长而减弱,吸附作用是长期过程,时间增长对其影响不显著。
    Abstract: Based on the theory of porous media and the continuum thermoelasticity, the migration mechanism of heavy metal ions in unsaturated soils under temperature gradient is studied. Firstly, the adsorption model for the lead ion by soil particles is determined through the non-isothermal adsorption tests. Secondly, based on the mass conservation, energy conservation and momentum conservation equations for various components in the three-phase system of solid, liquid and gas, a multi-physical field coupling mathematical model characterizing the migration of pollutants in unsaturated soils under temperature gradient is established. The model takes into account the adsorption effects of heavy metal ions, the migration of pore fluid, the compressibility of phase change soil skeleton and pore fluid phase and the deformation of soils regarding the temperature, porosity, pore water pressure, pore gas pressure, heavy metal ion concentration, adsorption concentration and soil skeleton displacement as the basic unknowns, which comprehensively reflect the characteristics of heat and mass migration and deformation in unsaturated soils. Finally, a soil column test on the migration of lead ions in unsaturated loess is carried out to verify the effectiveness of the theoretical model and the numerical results, and the influences of various factors on the migration laws of lead ions are analyzed. The results show that the Langmuir nonlinear adsorption model is more consistent with the adsorption characteristics of unsaturated loess for lead ions, and the maximum adsorption capacity decreases with the increase of temperature. The temperature gradient can obviously promote the migration of pollutants in unsaturated loess. While the pollutants migrate in unsaturated loess soil column, the whole soil column is in tensile deformation state. The migration of pollutants decreases with time. The adsorption is a long-term process, and the effects of time on it are not significant.
  • 土壤液化是国际上具有普遍意义的典型震害现象,是土木工程防灾减灾中的复杂难题之一。真实场地的液化响应规律及特征,仍是当今液化研究中没有解决好的突出问题之一,尚需深入研究[1]

    试验是液化研究中不可或缺的手段。然而,迄今尚没有一种试验,能够完全复现真实场地的液化响应过程及行为特征。室内液化试验,具有大量且突出的优点。可是,无论动三轴等土单元试验,还是振动台或动力离心机等模型试验,目前都不能同时、全面地对可液化土所处的现场条件和影响因素进行模拟,与真实液化场地存在差距。作为室内试验的一种补充,直接开展现场液化试验或现场液化监测,逐渐成为当今土壤液化研究的必需手段[2]。采用室内和现场液化试验两种不同的研究手段,通过不断地“相互对比、验证、补充”,可以进一步完善、修正当前的液化分析方法。

    相比室内液化试验方面的丰富成果,人工激振下现场液化试验的研究成果尚少,目前国际上总体处于探索阶段[2]。中国虽有现场液化试验方面的研究成果并取得一些开创性的进展[3-7],但是总体仍处于起步阶段,亟待进一步发展。

    值得注意的是,目前国内外已有的人工激振下现场液化试验,无论是基于原状土还是重塑土,可液化土的上覆压力基本不超过20 kPa,绝大多数在10 kPa以内,大致对应1~2 m的土层埋深。根据对1962—2010年中国大陆地区地震液化实例的统计,液化土层的埋深集中在0~6 m,约占总数的80%,其中2~4 m最为集中。而其地下水位埋深集中在0~3 m,约占总数的90%,其中1~2 m最为集中[8]。以此粗略推算,真实地震中液化土层所对应的上覆压力范围应集中在0~120 kPa,最为集中的均值约50 kPa左右。众所周知,上覆压力(或埋深)是控制土壤液化响应的主要因素之一。显然,当今开展的现场液化试验中可液化土的上覆压力较低,与真实液化场地存在较大差距。而且,由于试验操作的扰动,有效围压50 kPa以下的室内土单元试验成果相对很少。这也导致了已有现场液化试验成果难以和室内液化试验,尤其是动三轴等土单元试验,在相近的动、静应力水平下进行对比分析。可以说,过低的上覆压力,使得目前通过现场液化试验研究土壤液化问题的价值大打折扣。

    现有的振动台或动力离心机液化试验中,模型土的上覆压力和真实液化场地也会存在差距,往往通过堆载重物的方式来改善[9-10]。动力离心机试验,当应力水平无法满足时,在ng重力环境下,采用比土的密度大几倍的材料(如钢材)作为堆载重物,通过较小的堆载高度就会达到理想效果。然而,在振动台或现场试验等1g重力环境下,当所需上覆压力较大时,通过堆载重物的方式并不理想,甚至不可行。毕竟,试验中堆载重物的高度不可能无限制增加,而且实际加压效果、可操作性等都会遇到问题。若采用传统的“千斤顶+反力架”的方式,似乎非常容易实现对模型土施加上覆压力,且能完全满足液化试验需求的压力上限值。然而,对于液化试验而言,这种传统的加压方式并不可行。因为液化发生前,土体一定会率先出现沉降变形,微小的沉降都会导致通过千斤顶施加的荷载无法传递,造成上覆压力“瞬间消失”。因此,需专门发展试验技术,才能开展考虑上覆压力影响的现场液化试验。

    鉴于此,笔者在已有现场液化试验研究基础上[6-7],围绕能够开展“多工况、可重复、易操作、低成本”试验的基本设想,专门研制施加上覆压力的设备,发展考虑上覆压力影响的现场液化试验技术,设计完成2次不同上覆压力(25,50 kPa)的饱和砂、干砂对比试验,分析其动力响应规律及特征,并与已有认识对比,验证试验技术的可靠性。本文研究成果,可为当今人工激振下现场液化试验的发展提供技术借鉴,也为在1g重力环境下,以更接近真实场地的模拟条件开展土壤液化问题研究,提供一种简单实用、可操作性强的试验新方法。

    为评价考虑上覆压力影响的现场液化试验是否成功实施,本文首先提出试验孔压比的概念,并阐述其必要性和合理依据。

    由于现场液化试验是在1g重力环境下进行,上覆压力必须通过特定技术手段施加给试验土体。在动荷载开始前,预先施加给土体的上覆压力能够较容易地保持平稳和恒定状态。但是,在动荷载作用下,人为施加的上覆压力会发生很大变化,尤其是采用机械设备以“接触力传递”方式施加的上覆压力。理论上,液化的土体会软化和变形,通过机械设备一直维持初始施加的上覆压力的平稳,几乎不可能,因为液化的土体无法提供良好的反力支撑。但是,在液化发展阶段,尤其是超静孔压急剧上升时,如果预先施加的上覆压力下降(损失)过快、过大,则必然会对现场液化试验结果产生较大影响,甚至失去考虑和控制上覆压力的意义。因此,在现场液化试验中,施加给土体的上覆压力应保持大致的稳定。只有预先施加的上覆压力损失有限,以此得到的现场液化试验结果才有价值。而且,需要对现场液化试验中上覆压力的损失程度,进行定量评价和判断。

    一方面,根据已有试验结果和液化认识,当孔压比达到0.4~0.6左右,是液化变形及对地震动产生影响等动力响应特征的阈值[6, 11-12]。另一方面,从液化机理上,当土壤液化的条件均能充分满足时,超静残余孔压的发展是不可能大于其上覆压力的。也就是说,如果设备预先施加的上覆压力下降(损失)过快、过大,超静孔压一定不会发展过高。

    因此,在土壤液化的条件均能满足的情况下,只要现场液化试验中的超静孔压与预先施加的初始上覆压力的最大比值能够达到0.4~0.6及以上,就可获取该上覆压力作用下土壤液化的主要响应特征。本文将该最大比值,定义为试验孔压比,以区别通常描述液化特征的孔压比。因此,试验孔压比,可初步判断预先施加的上覆压力是否大致稳定,从而评价可考虑上覆压力影响的现场液化试验是否成功实施(有价值)。

    当然,当土体抗液化强度较高或动荷载较小时,液化发展程度本身有限,就不会发生较大的土体软化或变形,预先施加的上覆压力自然也会保持较平稳状态。此时,则无需采用试验孔压比进行评价。

    本文现场液化试验,所用的砂土全部采用福建标准砂。通过室内土工试验实测,其物理指标如表 1所示,其颗粒级配曲线如图 1所示。

    表  1  试验砂物理指标
    Table  1.  Physical indexes of sands used in tests
    最大干密度ρdmax/(g·cm-3) 最小干密度ρdmin/(g·cm-3) 有效粒径d10/mm 特征粒径d30/mm 控制粒径d60/mm 不均匀系数Cu
    1.953 1.538 0.309 0.475 0.889 2.877
    下载: 导出CSV 
    | 显示表格
    图  1  试验用砂颗粒级配曲线
    Figure  1.  Grain-size distribution curves of sands used in the tests

    本文现场液化试验所用的动力加载设备,与笔者以往开展的现场液化试验相同[6-7]。该设备是自主研发,主要有由激振装置、刚体基础、变频控制器、电动机等组成,能在试验场地中激励产生水平向的等幅正弦波地震动。加速度传感器型号为LC0113M,量程±2g,灵敏度2500 mv/g左右。孔压传感器型号为HC-25,量程0~50 kPa,灵敏度0.2%FS;土压力传感器型号HC-16,量程0~150 kPa,灵敏度0.2%FS。位移传感器型号为HCW-20,量程0~15 cm,灵敏度0.1%FS。动态数据采集设备,采用日本东京测器研究所生产的TMR多功能记录仪,最高数据采样率可达100 kHz。部分试验设备及传感器如图 2所示。由图可知,所用的加速度、孔压、土压力传感器体积小,能够满足土体“点测量”的基本要求。

    图  2  现场试验用的部分试验设备及传感器
    Figure  2.  Photos of experimental apparatus and sensors used in in-situ tests

    施加上覆压力的设备,是本文现场液化试验能否成功实施的关键硬件。笔者自主研发了适用于现场液化试验的施加上覆压力的设备。该设备主要由钢筋混凝土反力架、液压千斤顶、加载钢板及持压装置组成,其示意图及现场实景如图 3所示。该设备,是在传统的“反力架+千斤顶”模式的基础上,增加了自主研发的持压装置,解决了因土体沉降导致上覆压力“瞬间消失”的问题。持压装置的工作原理,是通过合理刚度的弹簧的回弹来补偿沉降,从而实现“伺服”控制。通过该设备,现场液化试验所需的上覆压力,不但可以被定量施加,而且能够在液化发展阶段保持大致稳定的状态。该设备是无源工作,操作简单,尤其适合野外试验环境。

    图  3  施加上覆压力设备的示意图及现场实景
    Figure  3.  Sketch and photo of apparatus for applying overburden pressure in in-situ liquefaction tests

    本文试验的主要目的,是在已有重塑土现场液化试验基础上,围绕“多工况、可重复、易操作、低成本”的基本设想,发展考虑上覆压力影响的现场液化试验技术并验证其可靠性。考虑到本文试验的探索性,并为后续针对性地开展具体液化问题的研究奠定认识基础,本次试验方案设计遵循以下基本原则:①设置同工况的干砂试验,作为对比参照;②仅把上覆压力作为试验变量,其他影响液化因素暂不考虑。如前文所述,考虑当前上覆压力基本现状,优先设计上覆压力分别为25,50 kPa的两次工况的现场液化试验。

    基于上述设想,在试验场地直接开挖2个相同尺寸的试坑,然后在坑内分别制备干砂、饱和砂。人工激励产生的水平向震动,通过刚体基础周围的土层传递,以“地震动”的形式直接对试坑内的砂土施加动荷载,因此不存在振动台试验的动力边界效应问题。由于试坑尺寸的大小不会对试验技术可行性产生实质影响,经综合考虑,试坑的尺寸确定为40 cm×40 cm×50 cm。试坑和刚体基础的空间相对位置关系如图 4所示。

    图  4  现场液化试验的试坑布置图
    Figure  4.  Sketch of test pits used in in-situ liquefaction tests

    本文设计试坑内砂土的初始相对密实度Dr为40%。砂土的制备技术及操作方法,与笔者之前进行的现场液化试验相同[6-7]。饱和砂,仍采用水沉法制备,其内部的传感器布置方案如图 5所示。干砂内部,除不放置孔压传感器外,其他布置均与饱和砂相同。考虑到试验的探索性,砂土内部的传感器均采用“备份”原则布置。即在砂土中心和左边的不同深度测点,逐一对应,布置相同类型的传感器。在干砂试坑一侧的地表,与砂土表面中心对应(均距离基础边缘70 cm),单独埋设1个加速度传感器,其实测加速度作为试验的输入地震动。位移传感器固定于反力架上端,通过加载钢板来测量砂土的地表沉降。

    图  5  饱和砂土的传感器布置图
    Figure  5.  Sketch of sensor distribution for saturated sands

    为表述方便,本文上覆压力为25 kPa工况的试验命名为Test-L试验,上覆压力为50 kPa工况的试验命名为Test-H试验。实际中,施加至设计上覆压力后,砂土均有一定程度的沉降。根据沉降值估算,在Test-L试验中,干砂的密实度变为49%,饱和砂密实度变为46%。在Test-H试验中,干砂的密实度为54%,饱和砂变为50%。砂土密实度的变化程度,均在可接受的范围。动荷载频率,与笔者之前的现场液化试验相同,仍为15 Hz。输入地震动的加速度幅值为0.11g左右,数据采样频率为100 Hz。两个工况现场试验分别加载。

    根据笔者观察,加载期间,无论Test-L试验,还是Test-H试验,饱和砂土均出现明显的喷砂冒水现象。图 6为Test-H试验后的饱和砂土表面实景,可以看出大量的砂堆积在加载钢板的一角,是加载期间喷出。相比笔者以前的现场液化试验,本文试验中饱和砂土的喷水高度、喷砂量均有大幅增加。这应该是,较大的上覆压力,导致砂土液化时超静孔压明显升高的结果。试验宏观现象初步表明,两次现场液化试验均取得了成功。加载期间,除去Test-L试验干砂中的土压力时程数据因电源线意外故障导致采集失败外,其他测点均成功采集到了数据。

    图  6  现场液化试验后的饱和砂实景(Test-H工况)
    Figure  6.  Photo of saturated sands in in-situ liquefaction tests after dynamic loadings (case of Test-H)

    图 7分别给出了两个工况现场试验中试坑外侧的地表加速度时程。经对比,试坑内砂土底部的加速度与其一致,限于篇幅不再给出,并将该数据视为试验的地震动输入加速度。可以看出,两次试验的输入加速度的幅值几乎一致,均在0.11g~0.12g左右。

    图  7  试验的输入加速度时程
    Figure  7.  Time histories of input acceleration in tests

    超静孔压是液化发展程度最直接的反映和标志。由于试坑内饱和砂土中的孔压传感器距离真实地表仅15 cm,如前文所述,倘若没有额外施加的上覆压力,超静的残余孔压是不可能超过1.5 kPa的(即大致15 cm高度饱和砂的有效自重)。图 8给出了两次现场液化试验饱和砂中心测点的超静孔压时程。结合图 7的地表加速度时程可知,在地震动开始的最初3 s内(约45个荷载循环),超静孔压迅速发展并达到最大值,之后孔压开始逐渐消散,这与以往的现场液化试验孔压响应规律和认识一致[6-7]。所不同的是,超静孔压到达最大值后,便表现为消散趋势且速率较快,这应该与实际排水路径过短(此处仅15 cm)有关。从图 8中可以看出,Test-L试验中超静孔压的最大值约为17 kPa,与其初始上覆压力23 kPa的比值,即试验孔压比,为0.75。Test-H试验中超静孔压的最大值约为25 kPa,与初始上覆压力50 kPa的比值为0.5。超静孔压的响应表明,现场液化试验中,设备施加给饱和砂的上覆压力能够保持大致稳定。这也说明两次现场液化试验均是成功的,可以结合土压力、地表沉降等响应做进一步分析。

    图  8  现场液化试验的超静孔压时程(中心测点)
    Figure  8.  Time histories of excess pore pressure in in-situ liquefaction tests(at central measuring point)

    另外,从理论角度,上覆压力是影响液化响应的重要因素。相比Test-L试验,在相同的加速度幅值(0.11g左右)和相对密度(50%左右)下,Test-H试验中的试验孔压比相对较小,应该还与饱和砂本身的上覆压力过大(50 kPa)有较大关系。后文,结合上覆压力的响应,这点也被证实。为后续分析需要,图 9分别给出两次现场液化试验中饱和砂中心和左边测点的超静孔压时程。为图示清晰,仅保留残余孔压。由图 9可知,同一工况试验中,两测点的孔压时程几乎完全一致。

    图  9  现场液化试验饱和砂中心和左边测点超静孔压时程
    Figure  9.  Time histories of excess pore pressure at central and left measuring points in the in-situ liquefaction tests

    砂土中的上覆压力通过土压传感器实测得到。上覆压力,是本文试验非常关心的问题,直接关系试验技术的可行性,进行重点分析。

    图 10给出了两次现场试验的砂土中心和左边测点的土压力时程,也即上覆压力时程。低上覆压力工况Test-L试验中,因未能成功采集到加载期间的干砂的土压力时程数据,故无法给出。根据试验记录,加载前,Test-L试验干砂中心和左边测点实测上覆压力分别为23,20 kPa。从图 10可以看出,饱和砂中心和左边测点的上覆压力为25,21 kPa,两者基本一致。在高上覆压力工况Test-H试验中,干砂中心和左边测点的上覆压力分别为60,40 kPa,对应的饱和砂上覆压力约为61,42 kPa,中心和左边测点的上覆压力值差别较大。

    图  10  现场液化试验的上覆压力时程
    Figure  10.  Time histories of overburden pressure in in-situ liquefaction tests

    砂土的上覆压力,是由液压千斤顶通过加载钢板中心位置施加给砂土。笔者分析认为,中心和边侧压力值差别较大,主要原因应是加载板的刚度不够大导致了应力集中。只是荷载较小的情况下,加载板的刚度影响尚不明显。从理论角度,施加的上覆压力越大时,加载板的受力就越接近于点荷载模式。根据弹性力学的Boussinesq解答,砂土中应力分布的规律就应该是中间大,两侧小。定性上,实测结果与理论解答相符。由于土压力、孔压传感器体积很小,能够实现“点”测量,上覆压力和孔压均可实测得到,因此考虑上覆压力对液化响应的影响时,可以做到精确分析。基于简化考虑,本文施加上覆压力时,将中心和左边测点的上覆压力平均值视为最终的上覆压力。实际中,Test-L试验干砂、饱和砂的最终上覆压力分别为22,23 kPa,Test-H试验干砂、饱和砂的最终上覆压力分别为50,51 kPa,达到试验设计预期。

    为更好地分析上覆压力响应,略去实测数值的动态波动,并同时给出孔压比时程,如图 11所示,并将孔压急剧上升阶段前后的5~6 s时程详细给出。图 11给出的孔压比是每一时刻的残余孔压与对应的平均上覆压力的比值,这与通常液化分析中的孔压比的概念完全一致。

    图  11  现场液化试验的上覆压力与孔压比时程
    Figure  11.  Time histories of overburden pressure and excess pore pressure ratio in in-situ liquefaction tests

    图 11可以看出,地震动开始时(两次试验的时程上分别为20.5,13.0 s处),上覆压力均有短暂的小幅度增加,推测仍与加载钢板刚度不够大引起的弹性回弹有关。在动荷载作用下,干砂只会发生沉降,不存在土体软化问题。从其上覆压力响应可以看出,在整个加载过程中,Test-H试验中干砂的平均上覆压力缓慢下降,最终从50 kPa下降至40 kPa左右,而且中心和左边测点处的上覆压力下降趋势和程度基本一致。根据试验记录,Test-L试验中干砂的平均上覆压力,最终从22 kPa下降至14 kPa左右。干砂的上覆压力响应表明,施加上覆压力设备的工作性能是较为理想的。相比之下,饱和砂的上覆压力,随着超静孔压的急剧上升和土体软化,在短时间内出现了较为明显的下降,尤其是饱和砂的中心测点处。对于液化试验,重点关注孔压急剧上升阶段。当超静孔压发展到最大值时,Test-L试验中的平均上覆压力从23 kPa跌至18 kPa左右,孔压比最大值约为0.85;Test-H试验中的平均上覆压力从51 kPa跌至35 kPa左右,孔压比最大值约为0.68。由于Test-L和Test-H工况试验中,仅存在上覆压力上的差别,即使考虑加载期间上覆压力的下降,从图 11中可知,当孔压比发展到最大时,Test-H试验中的上覆压力最低值也大于Test-L试验中的初始上覆压力。因此,这充分说明,Test-H试验中孔压发展程度较低(即孔压比相对较小),是上覆压力较高导致。也就是说,Test-H试验中的试验孔压比虽然为0.5,若在更大幅值的加速度作用下,还会增大的。这也表明,其他条件相同时,上覆压力越大,饱和砂的抗液化强度越高,这与目前对液化的理论认识和实际震害调查结果一致。

    图 11还可以看出,同一工况试验中,饱和砂中心和左边测点的上覆压力响应不同,在初始压力值、下降程度两方面均有差异。尤其是,在孔压上升阶段,中心测点处上覆压力的下降程度明显大于左边测点,这应该与边侧的饱和砂受到非液化土的约束有关。为进一步分析上覆压力对液化响应的影响,将两测点处的孔压比分别计算,即为该测点处每一时刻的残余孔压与对应的上覆压力的比值,如图 12所示。为了图示清晰,仅给出孔压急剧上升阶段前后的5 s时程。由图 12可知,Test-L试验饱和砂中心和左边测点的最大孔压比分别为1和0.7,Test-H试验饱和砂中心和左边测点的最大孔压比分别为0.75和0.6。Test-L试验中,砂土中心测点处的孔压比为1,这说明此处的砂土已经完全液化。

    图  12  现场液化试验饱和砂中心和左边测点的孔压比时程
    Figure  12.  Time histories of excess pore pressure ratio at central and left measuring points in in-situ liquefaction tests

    值得注意的是,Test-L和Test-H工况试验中,两者上覆压力的差别本身并不大,仅有25 kPa,但是在地震动加速度幅值0.11g下,饱和砂的液化程度(孔压比)却有较大差异。而且,对于同一工况试验而言,饱和砂中心和两侧的液化发展程度也不一样。从图 12中也可看到,Test-L和Test-H试验中,中心和左边测点的孔压比时程曲线均存在交叉现象。由此可知,液化发展初期(边侧上覆压力较小),边侧的孔压比是大于中心处的,随着发展(边侧上覆压力变大),边侧的孔压比是小于中心处的。这说明,在此工况条件下,上覆压力对液化的影响较为明显。尤其值得注意的是,Test-L工况试验中,中心和左侧处的上覆压力差值仅有4 kPa(相当于40 cm厚的饱和土层),但是时程曲线上仍表现出交叉现象,定性上的响应规律与Test-H工况是一致的。这更表明,在上覆压力水平较低(至少50 kPa以下)时,上覆压力这一因素对液化响应的影响是敏感的。

    上覆压力50 kPa,大致对应真实场地中埋深2.5~5 m土层的埋藏条件。可以说,真实地震中,普遍且大量存在着上覆压力为50 kPa及以下的液化场地。但是,目前液化研究中,上覆压力(有效围压)50 kPa以下的液化试验成果相对很少,很多问题尚缺乏清晰认识。而且,当前一些实用的液化分析方法,考虑埋深影响的修正,也多是以上覆压力(有效围压)100 kPa作为基准参照。因此,低上覆压力水平(0~50 kPa)下的液化响应问题,应需重视并进行深入探究。

    本文现场液化试验的反力架是稳定可靠的固定端,因此能够方便、准确地获取整个加载期间的砂土地表沉降时程。

    图 13给出了两次现场试验的砂土地表沉降位移和孔压比时程,并将孔压急剧上升阶段前后的5~6 s时程详细给出。从图 13中可以看出,无论Test-L试验还是Test-H试验中,孔压迅速发展阶段,饱和砂发生的地表沉降了约为3.5~4.5 mm,将该沉降值换算成体应变约为1%(砂土高度300 mm)。孔压刚开始发展(触发)时的沉降太小,实测的位移难以准确甄别,但是可以确定当地表沉降在0.3~0.5 mm时,即体应变约为0.1%时,孔压比发展到约0.1左右。由此可知,当体应变在0.1%以下时,孔压就已经开始发展。当体应变在0.1%~1%时,孔压快速地大幅度发展。本文得到的孔压与应变关系的定量结果,与笔者之前的现场液化结果基本一致[6],也与目前国际上普遍认同的Dobry等人有关应变与孔压关系的定量认识基本一致[13]。需说明的是,Dobry等试验结果的剪应变γ,是根据土样轴向应变ε,通过γ=(1+ν)ε的理论关系换算得到,其泊松比ν取为0.5。由于前面的系数仅为1.5,因此剪应变和侧限条件下体应变的量级不会发生变化。因此,定量认识上两者是一致的。这表明,本文试验中获取的地表沉降响应是可靠的,能够反映液化的主要变形特征,并可做进一步分析。

    图 13中还可以发现,地震动开始阶段,同等工况下的干砂、饱和砂的沉降量和沉降速率几乎完全一致。但是,当饱和砂的孔压比发展到0.6左右时,两者的沉降响应开始出现明显差异。此后一段时间,饱和砂的沉降量和沉降速率要明显小于同等工况下的干砂,这可能与饱和砂的超静孔隙水无法及时排出有较大关系,而干砂内的气体却能随时排出,从而表现出这种沉降差异现象。随着孔压的逐渐消散,饱和砂的沉降量和沉降速率明显增大,并最终大于干砂。由此可知,饱和砂液化后的沉降主要发生在孔压消散阶段。这一沉降规律,与他人的动力离心机液化试验一致,也与真实地震现场观察到的现象一致[14-15]。这也说明,本文现场液化试验,能够较理想地再现真实液化场地的地表沉降响应特征。

    图  13  现场液化试验的地表沉降与孔压比时程
    Figure  13.  Time histories of soil surface settlement and excess pore pressure ratio in in-situ liquefaction tests

    通过对宏观试验现象以及超静孔压、土压力(上覆压力)、地表沉降等响应的分析,并与已有的液化认识对比,本文设计完成的两个不同上覆压力工况的现场液化试验,均可得到液化的主要响应特征和规律,并能够以此获取新认识。这说明,本文提出的考虑上覆压力影响的现场液化试验技术,总体上是可行和可靠的,能够为在1g重力环境下,以更接近真实场地的模拟条件开展液化相关问题研究,提供一种可能的新手段。

    本文在已有现场液化试验工作基础上,围绕“多工况、可重复、易操作、低成本”的基本设想,提出了一套考虑上覆压力影响的现场液化试验技术,设计完成了2次不同上覆压力(25,50 kPa)的饱和砂、干砂对比现场试验,分析了其动力响应规律及特征,并与已有认识对比,验证了试验技术的可靠性。主要得到以下3点结论。

    (1)新提出的试验孔压比这一概念,可作为评价考虑上覆压力影响的现场液化试验是否成功的基本指标。两个工况的现场试验表明,新提出的考虑上覆压力影响的现场液化试验技术可行,能够得到土壤液化的主要响应特征。这为1g重力环境下,以更接近真实场地的模拟条件开展液化相关问题研究,提供了一种可能的新手段。

    (2)在低上覆压力水平(0~50 kPa)下,大致对应真实场地0~5 m的埋深,上覆压力是影响土壤液化响应的敏感因素。实际地震中,普遍且大量存在上覆压力50 kPa及以下条件的液化场地,而目前与该静应力水平一致的液化研究成果并不丰富,相关认识也不清晰。低上覆压力水平(0~50 kPa)下的土壤液化响应问题,应需重视并进行深入探究。

    (3)当饱和砂的体应变在0.1%以下时,孔压已经开始发展;当饱和砂体应变在0.1%~1%时,孔压迅速且大幅度发展。当孔压比发展到0.6左右时,同等工况下的干砂、饱和砂的地表沉降响应出现差异。并且,液化砂土的沉降主要发生在孔压消散阶段。

  • 图  1   Langmuir吸附模型拟合曲线

    Figure  1.   Fitting curves of Langmuir adsorption model

    图  2   土水特征曲线及试验装置

    Figure  2.   Soil-water characteristic curves and test device

    图  3   浓度、温度对比图

    Figure  3.   Comparison of concentration and temperature

    图  4   浓度对比图

    Figure  4.   Comparison of concentration

    图  5   孔隙水压、基质吸力、水蒸气压强及水蒸气密度变化曲线

    Figure  5.   Variation curves of pore water pressure, matrix suction, water vapor pressure and water vapor density

    图  6   各点位移、有效应力变化曲线

    Figure  6.   Variation curves of displacement and effective stress at each point

    图  7   高度方向浓度、吸附浓度变化曲线

    Figure  7.   Variation curves of concentration and adsorption concentration in height direction

    表  1   铅离子吸附试验数值

    Table  1   Test values of Pb ion adsorption

    温度/
    K
    初始浓度C0/(mg·L-1) 液相浓度C/(mg·L-1) 吸附Cs/(mg·kg-1)
    278 100 0.36 1244.51
    120 6.48 1619.98
    150 13.93 1741.25
    180 22.98 1915.12
    200 32.34 2021.25
    283 100 0.42 1216.86
    120 7.28 1548.64
    150 14.86 1646.25
    180 24.64 1865.46
    200 33.48 1986.46
    293 100 0.56 1186.88
    120 8.64 1508.24
    150 15.58 1604.25
    180 26.12 1768.56
    200 34.64 1868.36
    313 100 0.78 1108.46
    120 10.16 1428.62
    150 16.88 1516.26
    180 28.12 1648.36
    200 36.46 1726.54
    下载: 导出CSV

    表  2   吸附系数

    Table  2   Adsorption coefficients

    温度
    T/K
    非线性吸附常数
    α/(mg·L-1)
    最大吸附量β/
    (mg·g-1)
    278 0.55 2.32
    283 0.50 2.01
    293 0.56 1.82
    313 0.37 1.79
    下载: 导出CSV

    表  3   计算参数表

    Table  3   Parameters for calculation

    参数名称 取值 参数名称 取值
    土体试样高度H/m 0.3 土体密度ρs/(kg·m-3) 1560
    初始时刻污染物
    浓度C0/(kg·m-3)
    0.5
    污染物扩散
    系数D0/(m2·s-1)
    5×10-9
    土体试样半径r/m 0.1 土体纵向
    弥散度αL/m
    0.35
    n0 0.423 干燥气体摩尔
    质量Ma(kg·mol-1)
    2.88×10-2
    气体摩尔常数R/(J·mol·K-1) 8.214 孔隙度n0
    渗透率K0/(1/m2)
    1×10-16
    最大饱和度Ss 0.9 残余饱和度Sr 0.1
    吸湿膨胀系数βs 1.5×10-4 弹性模量Es/MPa 25
    泊松比μ 0.3 非线性吸附常数α/ (mg·L-1) 0.55
    最大吸附量β/(mg·g-1) 2.32 温度T0/K 278.15
    正反应速率
    常数λ1/(1·s-1)
    2.5×10-6 P0/Pa 0.9×105
    逆反应速率
    常数λ2/(1·s-1)
    1.5×10-6 βm/(1·K-1) -2.1×10-3
    土体弯曲因子τ 0.7 液态水质量热容Cw/
    (J·kg·K-1)
    4180
    空气质量热容Ca/
    (J·kg·K)
    1000 土体骨架质量热容Cp/ (J·kg·K-1) 1000
    水蒸气质量热容Cg/(J·kg·K) 1900 混合气相导热系数λg/ (W·m·K-1) 0.024
    液态水导热系数λl/(W·m·K-1) 0.56 土体骨架导热系数λs/ (W·m·K-1) 2.93
    液相水相变潜热Lwg/(J·kg-1) 2.4×106 水蒸汽摩尔质量Mw/
    (kg·mol-1)
    0.018
    热膨胀系数
    βsT/(1·K-1)
    7.8×10-6 液态水密度ρw/
    (kg·m-3)
    1000
    下载: 导出CSV

    表  4   边界条件

    Table  4   Boundary conditions

    初始条件 上边界条件 下边界条件
    孔隙率n 0.432 孔隙率n 孔隙率n
    孔隙水压Pl/kPa -190 孔隙水压Pl/kPa 孔隙水压
    Pl
    0
    孔隙气压Pg/kPa 102.1 孔隙气压Pg/kPa 102.1 孔隙气压
    Pg
    0
    位移u/
    m
    0 位移u/
    m
    位移u 0
    浓度C/
    (kg·m-3)
    0 浓度C/
    (kg·m-3)
    浓度C/
    (kg·m-3)
    0.5
    吸附浓度
    Cs/(kg·kg-1)
    0 吸附浓度
    Cs/(kg·kg-1)
    吸附浓度
    Cs/(kg·kg-1)
    温度
    T/K
    278.15 温度
    T/K
    278.15 温度
    T/K
    303.15
    下载: 导出CSV
  • [1] 罗春泳. 黏土的环境土工特性及垃圾填埋场衬垫性状研究[D]. 杭州: 浙江大学, 2004.

    LUO Chunyong. Study on Geoenvironmental Characteristics of Clay and the Behavior of Clay Liner in MSW Landfills[D]. Hangzhou: Zhejiang University, 2004. (in Chinese)

    [2] 蔡国庆, 赵成刚, 田辉. 高放废物地质处置库中非饱和缓冲层的热-水-力耦合数值模拟[J]. 岩土工程学报, 2013, 35(增刊1): 1-8. http://www.cgejournal.com/cn/article/id/15204

    CAI Guoqing, ZHAO Chenggang, TIAN Hui. Numerical simulation of coupled thermo-hydro- mechanical behavior for engineered barriers in high-level radioactive waste disposal[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(S1): 1-8. (in Chinese) http://www.cgejournal.com/cn/article/id/15204

    [3] 滕继东, 贺佐跃, 张升, 等. 非饱和土水气迁移与相变: 两类"锅盖效应"的发生机理及数值再现[J]. 岩土工程学报, 2016, 38(10): 1813-1821. doi: 10.11779/CJGE201610010

    TENG Jidong, HE Zuoyue, ZHANG Sheng, et al. Moisture transfer and phase change in unsaturated soils: physical mechanism and numerical model for two types of "canopy effect"[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(10): 1813-1821. (in Chinese) doi: 10.11779/CJGE201610010

    [4] 秦冰, 陈正汉, 方振东, 等. 基于混合物理论的非饱和土的热-水-力耦合分析模型Ⅰ[J]. 应用数学和力学, 2010, 31(12): 1476-1488. https://www.cnki.com.cn/Article/CJFDTOTAL-YYSX201012010.htm

    QIN Bing, CHEN Zhenghan, FANG Zhendong, et al. Analysis of coupled thermo-hydro- mechanical behavior of unsaturated soils based on theory of mixtures Ⅰ[J]. Applied Mathematics and Mechanics, 2010, 31(12): 1476-1488. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YYSX201012010.htm

    [5] 陈正汉, 郭楠. 非饱和土与特殊土力学及工程应用研究的新进展[J]. 岩土力学, 2019, 40(1): 1-54. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201901002.htm

    CHEN Zhenghan, GUO Nan. New developments of mechanics and application for unsaturated soils and special soils[J]. Rock and Soil Mechanics, 2019, 40(1): 1-54. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201901002.htm

    [6] 郭志光. 多孔介质中的水-热作用效应及污染物迁移特征研究[D]. 北京: 北京交通大学, 2018.

    GUO Zhiguang. The Hydro-thermal Effect and Transport Characteristics of Contaminants in Porous Media[D]. Beijing: Beijing Jiaotong University, 2018. (in Chinese)

    [7]

    ADALA M, BENNACER R, SAMMOUDA H, et al. Heat and mass transport in the unsaturated porous media: application to the soil dry drainage method[J]. Defect and Diffusion Forum, 2009, 283: 589-593.

    [8] 张嘉熙. 重金属污染物与悬浮颗粒耦合迁移过程及迁移距离研究[D]. 北京: 北京交通大学, 2019.

    ZHANG Jiaxi. Experimental Study on the Transport Process and Distance of Coupled Suspended Particles and Heavy Metal Ion[D]. Beijing: Beijing Jiaotong University, 2019. (in Chinese)

    [9] 王洪涛. 多孔介质污染物迁移动力学[M]. 北京: 高等教育出版社, 2008.

    WANG Hongtao. Dynamics of Fluid Flow and Contaminant Transport in Porous Media[M]. Beijing: Higher Education Press, 2008. (in Chinese)

    [10]

    HEITMAN J L, HORTON R, REN T, et al. A test of coupled soil heat and water transfer prediction under transient boundary temperatures[J]. Soil Science Society of America Journal, 2008, 72(5): 1197-1207. doi: 10.2136/sssaj2007.0234

    [11] 陈正汉. 重塑非饱和黄土的变形、强度、屈服和水量变化特性[J]. 岩土工程学报, 1999, 21(1): 82-90. http://www.cgejournal.com/cn/article/id/10260

    CHEN Zhenghan. Deformation, strength, yield and moisture change of a remolded unsaturated loess[J]. Chinese Journal of Geotechnical Engineering, 1999, 21(1): 82-90. (in Chinese) http://www.cgejournal.com/cn/article/id/10260

    [12] 秦冰, 陈正汉, 孙发鑫, 等. 高吸力下持水曲线的温度效应及其吸附热力学模型[J]. 岩土工程学报, 2012, 34(10): 1877-1886. http://www.cgejournal.com/cn/article/id/14878

    QIN Bing, CHEN Zhenghan, SUN Faxin, et al. Temperature effect on water retention curve under high suction and its modeling based on thermodynamics of sorption[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(10): 1877-1886. (in Chinese) http://www.cgejournal.com/cn/article/id/14878

    [13] 周凤玺, 周立增, 王立业, 等. 温度梯度作用下非饱和盐渍土水盐迁移及变形特性研究[J]. 岩石力学与工程学报, 2020, 39(10): 2115-2130. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202010015.htm

    ZHOU Fengxi, ZHOU Lizeng, WANG Liye, et al. Study on water and salt migration and deformation properties of unsaturated saline soil under temperature gradient[J]. Chinese Journal of Rock Mechanics and Engineering, 2020, 39(10): 2115-2130. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202010015.htm

    [14]

    RUTQVIST J, BÖRGESSON L, CHIJIMATSU M, et al. Thermohydromechanics of partially saturated geological media: governing equations and formulation of four finite element models[J]. International Journal of Rock Mechanics and Mining Sciences, 2001, 38(1): 105-127.

    [15] 陈宗芳. 吉林西部盐渍土冻结温度预测及水盐迁移规律研究[D]. 长春: 吉林大学, 2020.

    CHEN Zongfang. Study on Freezing Temperature Prediction and Water-Salt Migration of Saline Soil in Western Jilin[D]. Changchun: Jilin University, 2020. (in Chinese)

    [16] 陈星欣, 白冰, 蔡奇鹏. 考虑吸附-解吸效应的多孔介质中颗粒加速迁移问题解析解[J]. 岩土力学, 2015, 36(6): 1698-1706. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201506023.htm

    CHEN Xingxin, BAI Bing, CAI Qipeng. Analytical solutions for accelerated transport of particles in porous media with considering adsorption-desorption effect[J]. Rock and Soil Mechanics, 2015, 36(6): 1698-1706. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201506023.htm

    [17]

    TREMOSA J, GAILHANOU H, CHIABERGE C, et al. Effects of smectite dehydration and illitisation on overpressures in sedimentary basins: a coupled chemical and thermo-hydro-mechanical modelling approach[J]. Marine and Petroleum Geology, 2020, 111: 166-178.

    [18] 周凤玺, 高国耀. 非饱和土中热湿盐多场耦合过程分析[J]. 岩土工程学报, 2019, 41(5): 813-820. doi: 10.11779/CJGE2022S1010

    ZHOU Fengxi, GAO Guoyao. Multi-field coupling process of heat-moisture-salt in unsaturated soil[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(5): 813-820. (in Chinese) doi: 10.11779/CJGE2022S1010

  • 期刊类型引用(1)

    1. 梁越,李越,许彬,杨牛虎,张斌,潘剑,尹朝恒. 考虑围压影响的间断级配散粒土潜蚀规律试验研究. 河海大学学报(自然科学版). 2025(01): 80-86 . 百度学术

    其他类型引用(1)

图(7)  /  表(4)
计量
  • 文章访问数:  336
  • HTML全文浏览量:  49
  • PDF下载量:  156
  • 被引次数: 2
出版历程
  • 收稿日期:  2022-06-12
  • 网络出版日期:  2023-03-04
  • 刊出日期:  2023-08-31

目录

/

返回文章
返回