天津承压含水层粉土渗透特性试验研究

郑 刚，王佳琳，佟婧博，张天奇

天津承压含水层粉土渗透特性试验研究

郑 刚1, 2，王佳琳2，佟婧博1, 2，张天奇1, 2

(1. 天津大学建筑工程学院，天津 300072；2. 滨海土木工程结构与安全教育部重点实验室(天津大学)，天津 300072)

天津市地下工程施工主要集中在粉土层处，这类土层极易随地下水的流失而流失，在该类地区的地下结构在施工及运营期间也易出现渗漏现象．其中土体的工程性质与渗透特性是这种渗流侵蚀灾害的重要影响因素．因此很有必要对天津粉土的工程性质，尤其是渗透特性进行标定．为确定天津地下承压含水层中粉土的渗透特性，首先开展室内试验研究，适当地加工改造《土工试验方法标准》中规定的标准试验仪器，通过比重试验、颗粒分析试验、液塑限试验对天津粉土进行比重、级配、液塑限等基本参数进行标定，详细汇总整理各个试验数据，应用相关公式及数据处理方法对大量试验数据展开分析得出试验结果；然后通过常水头渗透试验与变水头渗透试验相结合的方式对粉土渗透系数的变化特征进行分析．研究结果表明：天津粉土的渗透系数随孔隙率的增大呈现增大的趋势，随时间的增长而减小且最终趋于稳定．从天津粉土的土体特性出发，将既无法传递净水压，也不能产生渗流的由结合水占领的无效孔隙排除，剩下的孔隙比定义为有效孔隙比．利用有效孔隙比的概念对经典的太沙基渗流理论进行修正，提出了适用于天津粉土的渗透系数估算公式．该公式能很好地计算天津粉土的渗透系数，可为天津市承压含水层中地下施工的灾害防治工作提供参考．

天津粉土；颗粒级配；比重；液塑限；渗透试验；太沙基

在以天津为代表的我国华北沿海地区，地下结构工程大都为建设在土质软弱、地下水位高并有粉土分布的基坑工程以及盾构法施工的隧道工程．粉土相对于砂土渗透系数小、颗粒小，易随地下水渗漏而发生流失土质，此类地质情况下，地下结构物极易因施工活动破坏承压含水层原有的应力平衡状态，导致承压含水层易随地下水流失的土体流动，引发严重的地下工程灾害[1]．因此，准确描述粉土的渗透特性具有重要意义．

目前对土的渗透性研究方面，砂土已经相对完备和成熟，前人分别从土的颗粒级配参数、密实度、固结系数、自由水动力黏滞系数等方面研究了其与土的渗透系数之间的关系，并提出了许多关于粗粒土渗透系数的经验公式[2]，如：太沙基渗透系数公式、柯森卡门渗透系数公式、斯托克斯孔隙流渗透系数公式[3]、中国水利水电科学研究院渗透系数公式、用固结度表示渗透系数的公式以及达西渗透系数公式等．这些公式都清楚地反映了渗透系数与各个影响因素之间的关系，使我们可以将试验得到的物理参数直接带入经验公式，从而得到较为准确的砂性土渗透系数．

对于细颗粒的黏土来说，国内外的研究也比较成熟．曾玲玲等[4]通过对压缩过程中重塑黏土渗透系数的试验数据进行多元线性回归分析，明确了压缩过程中重塑黏土渗透系数的主要影响因素是液限和孔隙比，并提出了相应的重塑黏土渗透系数定量表达式．

由于天津粉土的颗粒形状、矿物组成、颗粒级配、孔隙比等特征与上述细粒土有明显差别，以上基于室内试验对黏土的研究与分析，并不能满足粉土渗透特性的计算要求．除此之外，因为土体中存在没有连通或是半连通这两种水无法穿过的和连通但全部是结合水膜的无效孔径，但关于粗粒土的渗透系数经验公式并没有考虑无效孔径的问题．鉴于以上原因，本文基于天津粉土基本物理参数的测定，借助常水头试验和变水头试验结合的方法，系统研究其渗透系数变化规律，揭示地下工程灾害土的渗透特性．结合其工程性质建立一个适用于天津承压含水层粉土的渗透系数计算公式，以期为天津粉土渗透系数的估算搭建更加准确便捷的途径．

1 土样基本参数

天津市地下工程施工主要集中在粉土和粉质黏土层，在这种土层里发生的大多数工程灾害都是由于水土流失导致的．这种直接或间接地水土流失与天津市典型工程地质条件、天津市典型水文地质条件有着直接的关系，具体地质条件如图1、图2所示．取地质属于天津典型工程及水文地质的津南区浯水道与先锋河交口北侧的先锋河调蓄水池及配套管网工程事故场地的土体作为试验土体．在该场地埋深60.00m的范围内，地基土可以按照其成因不同被分为9层，按力学性质可进一步将其划分为15个亚层．工程影响深度内地基土自上而下为：①1杂填土、②2素填土、④1黏土、⑥淤泥质黏土、⑥4粉质黏土、⑦粉质黏土、⑧1粉质黏土、⑧2砂性大粉质黏土、⑨1粉质黏土、⑪1粉质黏土、⑪2粉砂、⑪4粉砂、⑫1粉质黏土、⑬1粉质黏土．根据地基土的岩性分层以及室内渗透试验结果，场地埋深60.00m以上可划分为：潜水含水层、第1微承压含水层和第2微承压含水层．本文研究的是标高-17～-22m地层(粉质黏土、砂性大的粉质黏土，地勘报告显示该土层有黏性)的粉土，即天津地下结构所处的典型土层的土体．

图1 天津市典型工程地质条件

图2 天津市典型水文地质条件

按照《土工试验方法标准》[8]，采用比重瓶法测定目标粉土颗粒比重为2.520．采用密度计法对过0.5mm筛、人工压碎以及机械搅拌3种方式处理过的⑧2粉土进行颗粒级配测定，由计算分析试验结果得到土体的颗分曲线(图3)可知，试验土样级配良好，不均匀系数u为1.28，曲率系数c为1.163．

对粉土进行液、塑限试验，测得其物理力学参数见表1．根据土样各种物性参数并结合塑性图得出结论：本试验采用土的种类为级配良好的低液限粉土．

图3 土样的颗分曲线

表1 试验土样的基本物理力学参数

Tab.1 Basic physical and mechanical parameters of soil

2 渗透试验

2.1 变水头渗透试验

本课题试验对象为粉土，故根据《土工试验方法标准》[8]应采用变水头渗透试验．渗透试验采用南京土壤仪器厂生产的TST-55型变水头渗透仪，见图4.为满足试验环境和操作简便性的需求，在保证试验结果准确性的前提下对传统渗透试验装置做出改造， 得到渗透试验使用的主要仪器设备有：变水头管、渗透容器、供水瓶、接水源管、进水管夹、排气管和出水管.

图4 试验装置

本试验共将土样按照质量由小到大分为7组，整个试验方案共涉及平行试验下7组孔隙率的试样在10个不同下降水头高度时的渗透系数测量与计算，具体试验方案如表2所示．

平行试验分组进行，设置空调温度，保证室温25℃的试验环境．使用上述工程现场取来的地勘报告中编号为⑧2的粉土层，因现场所取扰动试样有结块、含水量不均匀、掺混杂质等问题，故先对所取土样进行烘干处理，之后使用筛分法进行初筛去除较大杂质．再采取人工轻轻碾碎的方法对扰动土样进行处理，使颗粒大小均匀，用电子秤分别称取110g、115g、120g、125g、130g、135g、140g的人工处理土样．按照《土工试验方法标准》[8]进行试验，记录向供水装置中加水直到水头位置分别达到目标高度10cm、20cm、30cm、40cm、50cm、60cm、70cm、80cm和90cm时所需的时间．按照以上流程连续试验、记录4次，每次结束后使水头管的水位回升至初始目标高度，试验后取4次平均值做出记录．对试验数据的处理按照以下公式计算：

表2 变水头试验所对应孔隙率和孔隙比

Tab.2 Porosity n and pore ratio e corresponding to the variable head test

表3 试验时间记录

Tab.3 Test time record s

2.2 常水头渗透试验

本文旨在探究天津粉土在发生地下水渗流时的渗透性变化规律，且常水头试验更符合土体由于处于地下而水头不会发生变化的渗流情况．为更加系统全面地探究其渗透系数变化规律，揭示地下工程灾害土的渗透特性，采用常水头试验和变水头试验[9]配合测量土体渗透系数的方法，以求对天津承压含水层粉土做出全面细致的渗透性质测定．本试验采用与变水头试验相同的渗透仪以及滴定管进行组装后测定，设置空调温度，保证室温25℃的试验环境．使用同样的试样处理方法对土样进行烘干以及一系列处理，然后按照质量由小到大的顺序将试验土样分为7组，分别是125g、130g、135g、140g、145g、150g、160g，作为对照试验．按照《土工试验方法标准》[8]进行试验，记录在0min、10min、20min、30min、40min、50min、60min、70min、80min时经过300s出水管的出水量．对试验数据的处理按照以下公式计算：

表4 常水头试验数据

Tab.4 Data for constant head test

续表4

3 试验结果分析

3.1 孔隙率对渗透系数的影响

将变水头渗透试验中的两组平行试验分别设为A-1、A-2，图5揭示了对于下降水头高度相同的天津粉土，渗透系数随着孔隙率[10]的增大呈现增大的趋势．这一结果符合无黏性土渗透系数的基本理论，在渗透仪容积不变的前提下，粉土质量越大，孔隙率越小，即土体骨架中毛细通道越小，下降相同水头高度所需要的时间越长，计算所得的渗透系数越小．

图5 渗透系数随孔隙率变化曲线

3.2 时间对渗透系数的影响

试样土样进行常水头试验时，相同的孔隙率下渗透系数随时间的变化结果见图6．试验结果表明：在孔隙率固定不变的情况下，在短时间内渗透系数会随着时间的增加有减小的趋势，最终趋于稳定．这是因为无黏性细颗粒土体骨架内形成了多条毛细通道[11]，在渗透水流作用下颗粒迅速重新排布，细颗粒会在水流作用下将大颗粒组成的颗粒骨架中毛细通道填满，随时间的增长土体形成新的骨架并最终稳定[12].从数据上表现为渗透系数随着时间的增长而减小，最终当土体内结构基本稳定时，渗透系数为常数．

除此之外，图6显示时间由10s增长至70s，质量为125g的粉土试样渗透系数从1.7×10-4cm/s减小至1.4×10-4cm/s，减小了17.65%，质量为160g的粉土试样渗透系数从5.83×10-5cm/s减小至3.72×10-5cm/s，减小了36.2%．可知时间的增长对干密度较低的试样影响更为显著[13]．

图6 渗透系数随时间变化曲线

4 天津粉土渗透系数公式研究

4.1 无效孔隙的影响

事实上，相同相对密度的粗粒土和细粒土，细粒土的容重要远远小于粗粒土，也就是粗粒土的孔隙比要远远小于细粒土．这说明相同条件下，粗粒土的孔隙比更小，但渗透系数却更大．所以，除了土粒的大小和级配、土的密实度、水温和封闭气体的含量这些影响渗透系数的主要因素以外，本文认为细粒土的孔隙特征也与之有直接的关系．

土体中的孔径分为有效孔径和无效孔径[14]．有效孔径指的是可以发生渗流的孔径，无效孔径则指的是没有连通，或是半连通这种水无法穿过的和连通但全部是结合水膜的孔径．但关于粗粒土的渗透系数经验公式并没有考虑土体中无效孔径的问题．因此，简单地把粗粒土的渗透系数计算公式应用至天津粉土会造成很大的计算误差．

4.2 经典渗流系数的公式修正

土体中的结合水分为强结合水和弱结合水．强结合水位于最靠近土颗粒表面的地方，水化离子和水分子排列十分紧密；弱结合水不能传递静水压力，位于强结合水外层，与土体结合引力仅次于强结合水，可能会发生从浓度较高处缓慢地向浓度较低处迁移的运动，这种运动与重力无关．结合水占据了粉土孔隙的大部分，却无法传递净水压，也不能产生渗流[15].所以，将这部分由结合水占领的孔隙称为无效孔隙，用f来表示．

将用于计算粗粒土渗透系数的经典渗流公式中孔隙比修正为有效孔隙比，即可找到适用于天津粉土的渗透系数经验公式．其中，计算粗粒土渗透系数的经典公式有柯森-卡门公式、斯托克斯孔隙流渗透系数公式、中国水利水电科学研究院渗透系数公式、用固结度表示渗透系数的公式、太沙基渗透系数公式[16]以及达西渗透系数公式等．出于对已有经验公式计算所需试验参数的明确意义、误差大小以及试验操作难度的考虑，本文选取其中最简便、误差最小的经验公式——太沙基渗透系数公式进行修正．太沙基公式为

式中：为粗粒土的渗透系数，cm/s；10为粒径分布曲线上纵坐标为10%时所对应的土颗粒有效粒径，cm．

由于结合水占据的孔隙部分为无效孔隙，故将传统公式中的无效孔隙去除，得到适用于天津粉土的经验公式

4.3 粉土有效孔隙比的计算方法

根据第2节土样基本参数中的试验研究，天津粉土在含水率低于塑限时，土体内大部分是强结合水和空气；当粉土含水率高于塑限但低于液限时，土体内水大部分是弱结合水和少量的自由水；当含水量高于液限时，土体内出现大量自由水．

于是可以做出3个假定：①土体中结晶水的含量忽略不计；②土颗粒之间的连接只有结合水方式；③土颗粒和水不可压缩．

强结合水含量的最大值或弱结合水含量的最小值等于粉土的塑限，弱结合水含量的最大值比上换算系数等于粉土的液限[3]，进而可以通过已经测得的粉土的临界含水率求得其有效孔隙比．

天津粉土的无效孔隙比所对应的含水率为

式中：fw为无效孔隙比所对应的土颗粒的质量，kg；s为土颗粒的质量，kg；为结合水质量与土体刚好处于液限时含水量的比，也称换算系数；L为粉土液限.

可知粉土中结合水的质量为

粉土的无效孔隙比为

将式(7)、(8)代入太沙基公式(4)中，可得

故修正之后适用于粉土的太沙基公式为

4.4 算例验证

当粉土处于固体和半固体状态时，粉土中的强结合水会在电场的作用下与粉土颗粒紧密结合，形成强结合水膜．之后当粉土处于可塑状态时，含水量增加，强结合水膜外形成弱结合水膜．含水量继续增大，粉土易发生液化，自由水逐渐增多，变成流动状态．可知，强弱结合水的比值以及强弱结合水相对于自由水的多少会直接影响粉土的宏观物理状态.于是可以认为，粉土在可塑和软塑状态的临界点处结合水含量与含水率几乎相等，可以得到换算系数近似为1.4，即为液塑限之比，而换算系数代入修正之后的太沙基渗透系数公式中可得渗透系数值如图7所示．

将换算系数＝1.4代入第2节渗透试验不同孔隙率和时间的试验结果进行误差对比分析，不难发现，修正后的计算结果和室内试验结果相近，均在同一数量级．由图8中修正前后相对误差对比可知，在不考虑无效孔隙影响的情况下，所得天津粉土的渗透系数计算结果与实际渗透系数数量级相差40～80倍，说明没有经过修正的经验公式并不适用于天津粉土．由于太沙基公式应用范围广，所需试验参数简单易得，相对不易产生试验误差，故选择修正后的太沙基公式作为天津粉土的渗透系数经验公式．

图7 粉土的渗透系数

图8 相对误差随孔隙比变化曲线

5 结 论

(1) 采用液塑限联合测定仪对天津粉土进行界限含水率测定，结果显示天津粉土为低液限粉土．低液限粉土较一般土质毛细作用更加剧烈，水稳定性不佳，会对地下工程产生很大的影响．

(2) 通过对TST-55型渗透仪进行改装，完成土体的室内渗透试验，包括常水头试验和变水头试验．其中，由变水头试验得知，对于天津粉土而言，当下降水头高度相同时，随着孔隙率增大，渗透系数呈现增大的趋势．也就是粉土密度越大，孔隙率越小，下降相同水头高度所需要的时间越长，渗透系数也就越小．同一种土，当孔隙率固定不变时，渗透系数随时间的增长总体呈现减小的趋势，最终趋于稳定．由常水头试验可知，同种土在孔隙率相同的情况下，渗透系数会随着时间的增加呈减小的趋势，但最后趋于一个定值．同种土在同一时刻，渗透系数随孔隙率或孔隙比的增大而增大．渗透系数的数量级为 10-5～10-4cm/s．试验数据呈现出明显的规律性，除了部分数据有较小偏差外，其他数据基本符合预期结果.

(3) 根据适用于粗粒土的经典渗流系数公式，对比选择粗粒土渗透系数经验公式，即太沙基渗透系数经典公式．在此基础上对其进行修正．提出无效孔隙比的概念，解释了粉土相对粗粒土孔隙比大而渗透系数小的原因．

(4) 为了使修正后的太沙基公式可以应用到工程分析中，提出换算系数的概念，将换算系数与液塑限联系起来．经过比较分析取＝1.4，得到修正后的太沙基渗透系数公式，确定影响粉土渗透特性的主要因素．

(5) 用天津粉土对修正后的太沙基渗透系数公式进行验证，将试验测定的渗透系数与修正前后的太沙基公式计算结果对比之后发现，修正前的太沙基公式更适用于粗粒土；修正后的太沙基公式更适用于天津粉土，且计算所得渗透系数在同一数量级．总结得出适合于天津粉土的渗透系数估算公式，能够对后续的工程分析提供全面、合理的依据．

［1］ 张天奇，葛隆博，郑 刚. 砂土隧道开挖引起的地表及深层土体变形研究[J]. 天津大学学报(自然科学与工程技术版)，2019，52(增1)：113-119.

Zhang Tianqi，Ge Longbo，Zheng Gang. Deformation of surface and subsurface ground due to tunnel excavation in sand[J]. Journal of Tianjin University (Science and Technology)，2019，52(Suppl1)：113-119(in Chinese).

［2］ 谢康和，庄迎春，李西斌. 萧山饱和软黏土的渗透性试验研究[J]. 岩土工程学报，2005，27(5)：591-594.

Xie Kanghe，Zhuang Yingchun，Li Xibin. Laboratory investigation of permeability characteristics of Xiaoshan clay[J]. Chinese Journal of Geotechnical Engineering，2005，27(5)：591-594(in Chinese).

［3］ 宋林辉，黄 强，闫 迪，等. 水力梯度对黏土渗透性影响的试验研究[J]. 岩土工程学报，2018，40(9)：1635-1641.

Song Linhui，Huang Qiang，Yan Di，et al. Experimental study on effect of hydraulic gradient on permeability of clay[J]. Chinese Journal of Geotechni-cal Engineering，2018，40(9)：1635-1641(in Chinese).

［4］ 曾玲玲，洪振舜，陈福全. 压缩过程中重塑黏土渗透系数的变化规律[J]. 岩土力学，2012，33(5)：1286-1292.

Zeng Lingling，Hong Zhenshun，Chen Fuquan. A law of change in permeability coefficient during compression of remolded clays[J]. Rock and Soil Mechanics，2012，33(5)：1286-1292(in Chinese).

［5］ 刘维正，石名磊，缪林昌. 天然沉积饱和黏土渗透系数试验研究与预测模型[J]. 岩土力学，2013，34(9)：2501-2507，2514.

Liu Weizheng，Shi Minglei，Miu Linchang. Experimental study of permeability coefficient of natural saturated clay and its prediction model[J]. Rock and Soil Mechanics，2013，34(9)：2501-2507，2514(in Chinese).

［6］ 梁健伟，房营光. 极细颗粒黏土渗流特性试验研究[J]. 岩石力学与工程学报，2010，29(6)：1222-1230.

Liang Jianwei，Fang Yingguang. Experimental study of seepage characteristics of tiny-particle clay[J]. Journal of Rock Mechanics and Engineering，2010，29(6)：1222-1230(in Chinese).

［7］ Reddi L N，Thangavadivelu S. Representation of compacted clay minifabric using random networks[J]. Journal of Geotechnical Engineering，1996，122(11)：906-913.

［8］ GB/T50123—1999土工试验方法标准[S]. 北京：中国计划出版社，1999.

GB/T50123—1999 Standard for Soil Test Method [S]. Beijing：China Planning Press，1999(in Chinese).

［9］ 余良贵，周 建，温晓贵，等. 利用HCA研究黏土渗透系数的标准探索[J]. 岩土力学，2019，40(6)：2293-2302.

Yu Lianggui，Zhou Jian，Wen Xiaogui，et al. Standard exploration of permeability coefficient test for clay by HCA[J]. Rock and Soil Mechanics，2019，40(6)：2293-2302(in Chinese).

［10］ Li Gang，Li Xiaosen，Li Chao. Measurement of permeability and verification of Kozeny-Carman equation using statistic method[J]. Energy Procedia，2017，42(10)：4104-4109.

［11］ Xu Peng，Yu Boming. Developing a new form of permeability and Kozeny-Carman constant for homogeneous porous media by means of fractal geometry[J]. Advances in Water Resources，2007，31(1)：1-57.

［12］ 刘一飞，郑东生，杨 兵，等. 粒径及级配特性对土体渗透系数影响的细观模拟[J]. 岩土力学，2019，40(1)：403-412.

Liu Yifei，Zheng Dongsheng，Yang Bing，et al. Study on influence of particle size and gradation on permeability coefficient of soil by meso-scale simulation [J]. Rock and Soil Mechanics，2019，40(1)：403-412(in Chinese).

［13］ 王新志，王 星，胡明鉴，等. 吹填人工岛地基钙质粉土夹层的渗透特性研究[J]. 岩土力学，2017，38(11)：3127-3135.

Wang Xinzhi，Wang Xing，Hu Mingjian，et al. Study of permeability of calcareous silty layer of foundation at an artificial reclamation island[J]. Rock and Soil Mechanics，2017，38(11)：3127-3135(in Chinese).

［14］ 党发宁，刘海伟，王学武，等. 基于有效孔隙比的黏性土渗透系数经验公式研究[J]. 岩石力学与工程学报，2015，34(9)：1909-1917.

Dang Faning，Liu Haiwei，Wang Xuewu，et al. Empirical formulas of permeability of clay based on effective pore ratio[J]. Journal of Rock Mechanics and Engineering，2015，34(9)：1909-1917(in Chinese).

［15］ 高文静. 细粒土渗透特性的影响因素试验研究[D]. 西安：西安理工大学，2018.

Gao Wenjing. Influence of Clay Particles on Soil Permeability[D]. Xi’an：Xi’an University of Technol-ogy，2018(in Chinese).

［16］ 孔令伟，李新明，田湖南. 砂土渗透系数的细粒效应与其状态参数关联性[J]. 岩土力学，2011，32(增2)：21-26，41.

Kong Lingwei，Li Xinming，Tian Hunan. Effect of fines content on permeability coefficient of sand and its correlation with state parameters[J]. Rock and Soil Mechanics，2011，32(Suppl2)：21-26，41(in Chi-nese).

Experimental Study of Silt Permeability Characteristics of Tianjin Compressed Aquifer

Zheng Gang1, 2，Wang Jialin2，Tong Jingbo1, 2，Zhang Tianqi1, 2

(1. School of Civil Engineering，Tianjin University，Tianjin 300072，China；2. Key Laboratory of Coast Civil Structure Safety(Tianjin University)，Ministry of Education，Tianjin 300072，China)

The underground engineering construction in Tianjin is mainly concentrated in the silt soil layer，which is easily lost with the loss of groundwater. In addition，the underground structure in this area is prone to leakage during construction and operation.The engineering properties and seepage characteristics of soil are important factors affecting the seepage erosion hazard. Thus，it is necessary to calibrate the engineering properties of Tianjin silt，especially the permeability. To determine the permeability characteristics of silt in the underground confined aquifer of Tianjin，first，laboratory tests were carried out to properly process and transform the standard test instruments specified in the“Standard for Soil Test Method”. Basic parameters such as specific gravity，grade，and liquid plastic limit of Tianjin silt were calibrated through specific gravity test，particle analysis experiment，and liquid plastic limit test. Each test data was then summarized and sorted in detail. A large number of test data were analyzed using relevant formulas and data processing methods to obtain the test results. The variation characteristics of the silt permeability coefficient were then analyzed by the combination of constant head permeability test and variable head permeability test. Results show that the permeability coefficient of Tianjin silt increases with the increase of porosity，decreases with the increase of time，and finally tends to be stable. Based on the characteristics of Tianjin silt soil，the invalid pores that are occupied by bound water that cannot transfer water purification pressure and cannot produce seepage are excluded，and the remaining pore ratio is defined as the effective pore ratio. This paper uses the concept of effective pore ratio to modify the classical percolation theory of the Terzaghi foundation and proposes an applicable formula for estimating the permeability coefficient of Tianjin silt. The formula more precisely calculates the permeability coefficient of Tianjin silt. Results of this work could provide a reference for disaster prevention and control of underground construction of the confined aquifer in Tianjin.

Tianjin silt；particle gradation；specific gravity；liquid plastic limit；penetration test；Terzaghi

TU441+.33

A

0493-2137(2022)01-0077-08

10.11784/tdxbz202102034

2021-02-25；

2021-04-22.

郑 刚（1967—  ），男，博士，教授，zhenggang1967@163.com．

张天奇，tianqizhang@tju.edu.cn.

国家自然科学基金重点资助项目(41630641)；国家自然科学基金青年基金资助项目(51808387).

Supported by the Key Program of the National Natural Science Foundation of China(No. 41630641)，the Key Program of the National Natural Science Foundation for Young Scientists of China(No. 51808387).

(责任编辑：樊素英)