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

三星堆祭祀区地表干缩开裂病害发育后期发展演化特征

赵凡, 姚雪, 胡芮, 李思凡, 黎俊豪

赵凡, 姚雪, 胡芮, 李思凡, 黎俊豪. 三星堆祭祀区地表干缩开裂病害发育后期发展演化特征[J]. 岩土工程学报, 2025, 47(5): 1025-1035. DOI: 10.11779/CJGE20231288
引用本文: 赵凡, 姚雪, 胡芮, 李思凡, 黎俊豪. 三星堆祭祀区地表干缩开裂病害发育后期发展演化特征[J]. 岩土工程学报, 2025, 47(5): 1025-1035. DOI: 10.11779/CJGE20231288
ZHAO Fan, YAO Xue, HU Rui, LI Sifan, LI Junhao. Late-stage development and evolution characteristics of surface desiccation cracking in sacrifice area of Sanxingdui site[J]. Chinese Journal of Geotechnical Engineering, 2025, 47(5): 1025-1035. DOI: 10.11779/CJGE20231288
Citation: ZHAO Fan, YAO Xue, HU Rui, LI Sifan, LI Junhao. Late-stage development and evolution characteristics of surface desiccation cracking in sacrifice area of Sanxingdui site[J]. Chinese Journal of Geotechnical Engineering, 2025, 47(5): 1025-1035. DOI: 10.11779/CJGE20231288

三星堆祭祀区地表干缩开裂病害发育后期发展演化特征  English Version

基金项目: 

国家重点研发计划项目 2020YFC1522900

四川省科技计划项目 2023YFS0463

教育部人文社会科学研究项目 22YJCZH219

西南民族大学中央高校基本科研业务费专项资金项目 2022SCXTD01

详细信息
    作者简介:

    赵凡(1987—),男,博士研究生,副研究员,主要从事岩土文物劣化机理与保护技术方面的研究工作。E-mail: zhaofan_163@163.com

    通讯作者:

    姚雪, E-mail: yaozhibushixue0726@126.com

  • 中图分类号: TU449

Late-stage development and evolution characteristics of surface desiccation cracking in sacrifice area of Sanxingdui site

  • 摘要: 针对黏性土干缩开裂发育进程问题,以三星堆祭祀区考古遗址地表干缩开裂病害为对象,通过病害发育及赋存环境现场原位长时序监测,分析病害发育后期干缩裂隙发展演化特征,探讨病害发育与赋存环境的关联性、裂隙尺寸变化的作用方式和裂隙灌浆方法治理病害的可行性。结果表明,裂隙形态和尺寸变化主要表现为部分主干裂隙不断扩展及少量次级裂隙间歇性延伸或萌生,主干裂隙宽度和深度总体不断增大并趋于相对稳定,其尺寸越大,增大越明显。在主控环境因素地下土壤含水率时空变化驱动下,不同深度土体经历不同收缩阶段,加之裂块自身应力调整,抑制尺寸较小裂隙发育,导致裂隙尺寸变化呈现出时空分异性。采用裂隙灌浆治理病害从方法层面基本可行,关键在于选择裂隙发育相对稳定时间。
    Abstract: To explore the development process of desiccation cracking in clay soil, the desiccation cracking deterioration on the surface of the archaeological site in sacrifice area of Sanxingdui Ruins site is taken as the research object. The long-term and in-situ monitoring of desiccation cracking development is conducted along the concurrent environmental data collection. The data analysis focuses on studying the evolution characteristics of desiccation cracking at the late stage of deterioration development, exploring the correlation between desiccation cracking development and the environment, and investigating the feasibility of the grouting method based on the temporal characteristics of desiccation cracking. The results indicate that the development of desiccation cracking includes both morphological evolution and dimensional changes. Specifically, some primary main fractures demonstrate continuous expansion slowly, while a minority of the secondary fractures intermittently extend or emerge. Overall, the width and depth of the primary main fractures gradually increase and become stable over time. Larger fractures exhibit more obvious size increments during development. The moisture content of the soil is the predominant environmental factor influencing the desiccation cracking. Change in the controlling factor results in different development stages of fractures at varying depths within the soil, leading to intersecting fractures that form fracture blocks, and under the stress adjustment, smaller fractures are suppressed from developing, resulting in spatial and temporal heterogeneity in change of fractures sizes. Practical experience has demonstrated the feasibility of relieving desiccation cracking deterioration through the grouting method at stable periods of desiccation cracking development.
  • 地震作用引起的饱和砂土液化是导致地基失效等地震灾害的主要原因,因此预测土壤在地震作用下是否发生液化是液化灾害防治的首要任务。液化判别可以评估土体性能,对场地是否液化进行预测,以便及时应对可能发生的液化危害。而建立合理有效的砂土地震液化判别模型,是解决这一问题的重要手段[1]

    液化判别主要分为室内试验和现场测试。标准贯入试验是目前应用最广泛的现场测试方法,有大量现场试验记录的数据资料,其测试结果也比较可靠。近年来,以标准贯入试验为依据的地震液化判别方法的研究不断有成果出现,但是大部分的研究成果都是以Seed等[2]提出的“简化方法”为基础,通过采用更新的现场试验液化案例数据资料来对“简化方法”进行修正,再根据某一土层深度处地震引起的循环剪应力比与土层的循环抗力比的大小的对比结果来判断地震液化是否发生。Youd等[3]对Seed等提出的“简化方法”进行了更新和方法总结,该研究也成为后来学者们对地震液化判别方法进行改进的重要参考。例如:Idriss等[4]重新整理了重大地震液化场地的标准贯入试验资料,并基于更新数据资料对标贯试验的N值进行了校正;陈国兴等[5]合并了谢君斐和Idris等的数据资料,并基于该扩展数据集对Idriss等改进的“简化方法”进行了修正。

    与“简化方法”具有同样代表性并在国内广泛应用的液化判别方法是中国《建筑抗震设计规范》中的液化判别方法[6]。国内液化判别法自1970年试行版规范的提出到现行的2010版规范被广泛使用以来,进行了几次修改。由于原始资料的局限性,1989年版规范中的线性判别公式只适用于土层深度在10 m以内的情况;随着数据资料的丰富,2001年版规范判别式扩展到适用于土层深度在20 m内的分段线性公式,而现行2016年版《建筑抗震设计规范》采用的是2010年版规范中的对数型连续公式[7]

    由于采用中国建规方法在地震烈度为Ⅷ度和Ⅸ时对10~20 m深度处砂土的液化判别结果过于保守;而采用Seed等的“简化方法”判别液化时会出现标准贯入击数临界值随土层深度增加先增加后减小的不合理现象。因此,孙锐等[8]提出了双曲线型的液化判别模型,弥补了中国建规法和“简化方法”在砂土地震液化判别中的缺陷,但这个原始双曲线模型仅仅适用于不考虑细粒含量的影响的洁净砂地震液化判别。

    根据地震液化的震害调查结果发现,地震液化的砂土场地中通常都含有一定量的细粒土,为此研究者们通过室内试验来研究分析细粒含量对砂土液化的影响。马利超[9]通过动三轴试验和对砂土中细粒分布状态的微观观察发现随着细粒含量的增加,含细粒砂土中的细颗粒会逐渐分布在砂粒间的接触部位,增加了土体的抗液化强度。陈宇龙等[10]通过对含非塑性细粒的饱和砂土进行单调加载和循环扭剪试验发现随着细粒含量的增大,饱和砂土达到液化时的循环次数在减少。由于砂土中细粒的存在会影响土颗粒的分布状态和颗粒间的接触,从而影响土体的抗液化能力,因此对含细粒砂土的地震液化评估进行研究具有重要实际意义和工程价值。

    自1976年唐山地震中出现了含细粒砂土发生液化的案例以来,国内外学者展开了一系列针对含细粒砂土液化判别的研究。含细粒砂土中的“细粒”通常用黏粒含量(粒径小于0.005 mm)或细粒含量(粒径小于0.075 mm)来表示,研究中一般根据不同的土体类别选用不同的细粒指标[11]。关于这方面的研究,钟龙辉[12]基于唐山地震实际液化资料提出了轻亚黏土的液化判定方法,按黏粒含量的大小修正实际标贯击数,最终通过修正的标贯击数和实测标贯击数的对比分析来判定土壤是否液化。Seed等[13]在其1971年对“简化方法”的研究的基础上根据实测资料建立了细粒含量分别为15%和35%的土层的循环抗力比CRR与修正标贯击数(N1)60的液化临界曲线。刘颖[14]认为标准贯入击数与黏粒含量百分数的倒数(Cv=1/ρc)之间存在一定的函数关系,由此建议了建规法中考虑黏粒含量的修正项。董林等[11]通过分析集集地震数据资料,对“简化方法”和建规法中的含细粒影响项进行了改进,指出在对含细粒砂性土进行液化判别时应综合考虑黏粒含量和细粒含量的影响。这些学者提出的针对含细粒砂土的液化判别方法都是以建规方法和Seed等提出的“简化方法”中的洁净砂模型为基础展开研究,在考虑细粒含量影响的砂土液化判别模型中也存在前述相同的问题。

    本文基于原始的双曲线模型,根据大陆以往地震液化资料和台湾集集地震的标准贯入试验数据资料,建立了针对含细粒砂性土液化判别的计算公式,并利用国内外大地震液化数据集对提出的计算公式进行检验来验证所建判别模型的正确性,为含细粒砂土的地震液化判别研究的发展提供参考。

    双曲线模型的建立是基于大陆建规法形成和检验所采用的国内数据资料,来源于20世纪六七十年代大陆几次大地震的现场记录,主要包括唐山地震、通海地震和海城地震的液化资料[15]。筛除了实测标准贯入击数大于50的数据资料后,共有现场数据150例,其中液化样本95例,非液化样本55例,且该组成数据在各地震烈度区的分布较为均衡(Ⅶ度、Ⅷ度和Ⅸ度的样本数分别为50例、47例和53例),但其中90%的数据资料的土层深度在10 m以内,浅层土的数据资料占多数。

    1999年台湾集集地震数据是近年来记录的发生地震规模较大,液化资料记录较为详尽的数据资料,包含有许多含细粒砂土在不同深度处液化的记录情况[16]。该地震资料中的水平峰值加速度是由地震台记录的,较为准确且不需要估计;该数据资料包含了许多土层深度大于10 m的土壤液化资料,应用该数据资料可以有效地检验双曲线模型对不同烈度区中土层深度在20 m以内的液化判别情况的适用性。另外,台湾集集地震液化资料是为数不多的同时记录了场地土的黏粒含量ρc和细粒含量FC的数据资料,对含细粒砂性土的研究提供了宝贵的资料。

    Cetin等[17]将世界范围内大地震液化资料总结形成了Cetin数据集,该数据资料收集了各地区的原始记录情况,对各地区的原始资料进行了详细的评定,并修正了不同地区现场测试的标准贯入锤击数。该地震液化数据集包括113个液化场地和95个非液化场地数据以及2个边缘场地案例。

    本文分别采用这3种地震液化资料作为样本数据集对双曲线模型进行改进。首先采用原始数据资料作为基础数据集修正了原始双曲线模型,并对原始双曲线模型和修正双曲线模型的判别情况进行了对比分析;然后采用数据资料可靠且同时记录了场地的细粒含量和黏粒含量的集集地震数据集来确定砂土含细粒时液化判别的修正系数,建立了含细粒砂土的液化判别双曲线模型;最后采用丰富性和多样性较强的,并且对不同地区的液化资料进行了统一修正的Cetin数据集作为检验样本来检验双曲线型的液化判别模型在不同土层条件或不同地震烈度区的适应性。

    在孙锐等[8]一文中,利用大陆地震液化数据来验证所建立的双曲线模型的回判成功率,从其回判结果上看,用于砂土液化判别的双曲线模型相对于2010版建规判别法在液化和非液化的判别成功率上均有提高。另外,文献中收集了日本阪神地震和台湾集集地震数据资料,并对该模型的判别结果进行验证,结果显示所建立的双曲线模型保持了建规法对液化样本判别成功率较高的优点,还大大提高了非液化样本的判别成功率。图1显示了孙锐双曲线模型对各数据集的判别成功率。由于世界各国所使用的标准贯入试验标准不同,所测得的标贯击数N值需进行的修正也有所不同。中国地震原始资料中并未记录标准贯入试验的修正参数,故采用标贯试验的实测值进行分析;而台湾集集地震数据资料和Cetin数据集采用的是修正后的标准贯入击数值(N1)60作为标贯试验的实测值来进行判别对比。

    图  1  原始双曲线模型对各数据集的判别结果
    Figure  1.  Evaluated results of original hyperbolic model for data sets

    图1中原始液化判别双曲线模型对3种数据集的检验结果可以看出,该模型在评估大陆数据和集集地震液化案例时的效果较好,且其液化样本的判别成功率和非液化样本的判别成功率大致相等,均在80%以上;而在对Cetin综合数据集的判别分析中,非液化案例的判别结果良好,液化案例的判别成功率相对较低。分析原因有以下两点,一是综合数据集中有许多早期的液化数据,不同地区的标准贯入试验结果存在较大差异,N值修正后这种差异仍然存在;二是原始液化判别的双曲线模型选取的临界标准贯入击数初值偏小,也就是各烈度区的临界标贯击数基准值偏小,这样会导致原始双曲线模型对地表附近砂土的判别结果偏危险。本文针对这个问题通过对不同烈度区的实测标贯击数散点图的统计分析重新确定了临界标贯击数基准值的大小,并对液化判别双曲线模型中的土层深度和地下水位深度这两个液化判别影响项采用了先独立再组合的分解分析法,在此基础上提出了针对洁净砂地震液化判别的修正双曲线模型。

    本文采用与原始双曲线模型相同的大陆地震液化资料作为回归总样本,并基于原始双曲线模型[8]提出了修正的液化判别双曲线模型为

    Ncr=N"(1k3dw)(k2+dsds+k1)γ (1)

    式中 Ncr为标贯击数临界值;N'为临界标准贯入击数基准值;ds为标准贯入点的土层深度;dw为地下水位深度;k1,k2,k3为回归系数;γ为组合系数。

    修正的液化判别双曲线模型的确定主要分为三步:第一步是确定临界标准贯入基准值N',基于原始数据共150个样本点的地震烈度和实测标准贯入击数的统计关系,直观地确定液化和非液化之间的临界标贯击数值基准值N' [15];第二步是确定dwds对标贯击数临界值Ncr的影响,通过对原始样本点的回归计算,分别考虑dwds这两个因素对液化判别的影响,确定各影响项的回归系数;第三步是将这两个影响项组合,通过对整体样本的回归计算,确定了组合系数γ的值为1。从而得到了修正的砂土液化判别双曲线模型如下式所示:

    Ncr=N"(10.035dw)(0.7+dsds+6) (2)

    式中,确定临界标准贯入基准值N'时,参考建规法中确定临界标准贯入击数的方法[18],基于地表附近(ds=3 m,dw=2 m)的砂土发生液化时的临界标准贯入击数,并在采用的原始数据中增加了海城地震和唐山地震液化的资料数据,通过各样本点地震烈度和实测标准贯入击数的统计关系,最终确定临界标准贯入击数基准值N'在地震烈度为Ⅶ度、Ⅷ度、Ⅸ度时分别为9,12,17。

    为了验证修正双曲线模型的判别效果,将国内地震液化原始数据代入修正的双曲线模型来计算其回判成功率,并与2010版建规判别法和原始双曲线模型的回判成功率进行比较,得到了不同判别方法的回判结果如表1所示。通过分析这3种液化判别方法的回判结果可以发现,原始双曲线模型和修正双曲线模型对原始样本的回判结果都优于2010版建规判别法,并且从全部样本的回判结果上看,修正双曲线模型在保持非液化样本判别成功率的基础上,提高了对液化样本的判别成功率。

    表  1  液化回判成功率
    Table  1.  Success rates of liquefaction estimating  (%)
    烈度场地情况2010建规法原始双曲线模型修正双曲线模型
    Ⅶ度样本液化5974100
    非液化918778
    整体748090
    Ⅷ度样本液化938996
    非液化808590
    整体888893
    Ⅸ度样本液化859090
    非液化727278
    整体818486
    全部样本液化808595
    非液化828282
    整体818490
    下载: 导出CSV 
    | 显示表格

    一般情况下,含细粒砂是指由无黏性砂粒、粉粒和有黏性的黏粒组成的砂土,其中粉粒和砂粒的物理性质比较接近,而且含细粒砂土的液化性能主要受黏粒含量的影响[19],因此本文首先是将黏粒含量作为含细粒指标来判别液化。20世纪后全球发生的多起强震中都发现了含有非塑性细粒的深层土液化的案例。为了扩大该修正双曲线模型的应用范围,使其可适用于非塑性细粒土的液化判别,本文引入含细粒影响范围更广的细粒含量作为细粒指标,并将修正双曲线模型的含细粒影响项从按黏粒含量进行修正扩展到根据细粒含量进行修正。以此来判别含细粒砂土的地震液化情况。

    本文首先考虑黏粒含量对含细粒砂土液化判别的影响,并以黏粒含量的百分比值作为分类标准。参考中国2010版《建筑抗震设计规范》中判别饱和砂土液化的地震烈度要求,在地震烈度为Ⅵ度时砂土液化现象出现的可能性较低,危害性较小,故不考虑地震烈度为Ⅵ度时的砂土液化样本。表2给出了集集地震资料中272例土体样本按黏粒含量进行分类的样本统计结果。

    表  2  集集地震数据中按黏粒含量分类统计的样本统计结果
    Table  2.  Statistical results of Chi-Chi earthquake sample data by clay content
    黏粒含量ρc/%液化数非液化数合计
    ≤35557112
    4172037
    5231235
    614519
    712214
    85611
    9628
    10617
    11303
    12437
    >1213619
    总数158114272
    下载: 导出CSV 
    | 显示表格

    表2中,当砂土中的黏粒含量ρc ≤ 3%时,通常把该砂土作为“洁净砂”考虑,不对砂土的临界标准贯入击数值进行修正。而当土体中的黏粒含量ρc>12%时,黏粒的增加对标贯击数N值无明显影响,所以仍按ρc=12%考虑。由于黏粒含量ρc>6%的各分类组下的数据样本较少且不均衡,为了保证结果的稳定性,针对非液化样本中实测标贯击数N值普遍偏大的现象,选择筛去ρc=8%和ρc=12%的分类下偏离临界状态曲线较远的点。利用该统计分类方法可以分别计算各黏粒含量分类组所对应的临界标准贯入击数修正系数。

    关于临界标准贯入击数修正系数的确定[12],本文依据建规判别法的建立过程,通过对标贯击数临界值进行含细粒修正来确定临界标贯击数等效值,首先基于国内样本数据得到考虑黏粒含量的临界标准贯入击数修正系数,然后根据黏粒含量修正系数确定相应的细粒含量修正系数,最终形成考虑细粒含量的砂土液化判别的修正双曲线模型,具体过程如下:

    (1)首先,用α表示修正双曲线模型中计算的临界标准贯入击数的修正系数,以此来考虑黏粒含量大小对砂土液化的影响。Ncr等效表示按黏粒含量修正后的标准贯入击数临界值的等效值,Ncr表示由修正的双曲线模型计算得到的临界标贯击数,两者的关系如下式所示:

    Ncr=αNcr (3)

    具体地,利用MATLAB计算程序进行统计回归,以实测的标贯击数N值作为Ncr等效的目标值,由此得到各黏粒含量下的α值。

    (2)分析得出不同黏粒含量下的α值整体呈下降趋势。将不同黏粒含量下的修正系数α与黏粒含量的变化关系用图2表示。从图2可以看出,临界标准贯入击数的修正系数α随着黏粒含量ρc的增加而逐渐减小。假设αρc之间存在负相关的线性关系:

    α=t×ρc+h, (4)
    图  2  黏粒含量和修正系数线性关系图
    Figure  2.  Linear relation between clay content and correction coefficient

    式中,t,h为线性函数系数。由线性回归可得t=0.0555,h=1.2537。根据该线性关系即可得到各黏粒含量ρc下的修正系数α的线性结果值。

    (3)分别针对砂土中不同黏粒含量分类下对应的细粒含量值进行分析,得到它们各自对应的细粒含量的平均数,中位数,众数,以此来确定各黏粒含量对应的细粒含量范围,由于ρc值较大的部分数据样本较少,故不分析其众数。图3给出了砂土中不同黏粒含量下对应的细粒含量平均值和中位数的变化趋势,根据细粒含量统计参数随黏粒含量变化的变化趋势,直观地确定各黏粒含量对应的细粒含量值,即图中星号点所表示的对应值。

    图  3  细粒含量变化趋势图
    Figure  3.  Trends of fine grain content

    (4)根据图3砂土的细粒含量随黏粒含量的变化趋势图,确定黏粒含量与细粒含量的对应关系,结果如表3所示。当砂土中细粒含量的实际记录值非整数时,用插值法得到修正系数α[12]

    表  3  标贯击数的修正系数确定值
    Table  3.  Deterministic coefficients of SPT
    细粒含量FC/%≤1530354045
    黏粒含量ρc/%≤367912
    修正系数α1.090.920.870.750.59
    下载: 导出CSV 
    | 显示表格

    综上,将洁净砂液化判别的双曲线模型和含细粒砂土中细粒影响项以及其他因素影响项进行组合,得到了含细粒砂性土液化判别双曲线模型表示如下:

    Ncr=αN"(10.035dw)(0.7+dsds+6), (5)

    式中,Ncr为考虑细粒修正后的临界标贯击数等效值,α为临界标贯击数修正系数,按表3取得。当砂层实测的标准贯入值小于临界标贯击数等效值时,即N < Ncr时,判定为砂土发生液化,反之则判别为砂土不液化。

    为了检验修正的双曲线模型对洁净砂地震液化判别结果,将原始双曲线模型和修正双曲线模型对集集地震液化样本的判别情况进行了对比,对比结果如图4所示。图4把双曲线模型中地下水位dw=2 m时的结果作为典型的代表曲线,对原始双曲线模型和修正双曲线模型的评估结果进行了对比分析,并将样本点按烈度区(Ⅶ度、Ⅷ度、Ⅸ度)分别绘制在图4(a),(b),(c)3个坐标图中,图中实心点表示液化点,空心点表示非液化点。从临界标准贯入击数与土层深度的关系图中可以明显地看出,在地震烈度为Ⅶ度时,修正双曲线模型对液化样本的判别结果要优于原始双曲线模型;而在地震烈度为Ⅷ度时,修正双曲线模型从一定程度上改善了原始双曲线模型对非液化样本的判别结果;地震烈度为Ⅸ度时由于高烈度区样本案例不足,两种双曲线模型的对比结果不是很明显,但从总体上看修正双曲线模型不仅保持了原始双曲线模型判别结果良好的优势,也改善了原始双曲线模型在浅层处的液化临界线的偏于危险的情况,提高了液化判别结果的安全程度。

    图  4  临界标准贯入击数与土层深度的关系
    Figure  4.  Relationship between critical SPT number and soil depth

    为了分析修正的双曲线模型中增加含细粒修正项之后的对砂土地基的液化评估效果,将修正的含细粒砂土液化判别双曲线模型分别代入台湾集集地震液化样本和Cetin地震液化数据集样本进行判别分析,并将结果分别绘制在图5(a)5(b)中。图中以液化场地样本点的临界标贯击数Ncr为样本点的纵坐标,以实测标贯击数N值为样本点的横坐标来确定图中的基础样本点;以标准贯入试验的判断标准为依据,并确定N实测=α·Ncr的关系线为临界线,根据临界标贯击数修正系数α的范围确定临界线的位置区间来反映含细粒修正后的判别结果。从图5(a)5(b)的结果可以发现,在临界线区间范围内有许多液化点和非液化点,临界线区间范围左侧主要是液化点,右侧主要是非液化点。

    图  5  含细粒修正双曲线模型判别结果的检验
    Figure  5.  Results of evaluation for Cetin’s data set by corrected hyperbolic model

    图5的结果也反映了本文提出的考虑含细粒影响的修正双曲线模型可以有效地应用于含细粒砂土的地震液化判别。随着细粒含量的增加,临界线逐渐左移使得非液化区范围增大,非液化样本判别成功的数量就会增多,这主要是因为砂土中细粒的存在会贡献于土的骨架,使得土体的抗液化能力增强,利用现场试验测得的标贯击数就会增大,而采用修正双曲线模型得到的等效临界标准贯入击数值就会相对减小,这样有利于非液化样本的判别。

    同时,将该Cetin数据集样本代入修正的砂土液化判别双曲线模型中得到的判别结果与实际液化情况进行对比,由于数据集中只有细粒含量FC的记录资料,故只计算考虑和不考虑细粒含量FC影响的修正双曲线模型这两种情况的液化判别成功率。通过计算得到,不考虑细粒含量FC影响的修正的砂土液化判别双曲线模型对液化和非液化样本的判别成功率分别为71%(80/113)和85%(81/95),而图1中原始双曲线模型对液化和非液化样本的判别成功率分别为71%和82%,从对比结果可以看出,修正的含细粒砂土液化判别双曲线模型在保证液化样本的判别成功率的基础上,通过采用对临界标贯击数进行含细粒修正的方法,提高了Cetin综合数据资料中非液化样本的判别成功率;当对修正双曲线模型计算得到的标贯击数临界值采用细粒含量修正系数修正后,修正的含细粒砂土液化判别双曲线模型对液化和非液化样本的判别成功率分别为67%(76/113)和88%(84/95)。与图5(a)5(b)中反映的规律相同,修正双曲线模型对液化样本的判别成功率稍有降低,而其对于非液化样本的判别成功率结果有一定提高,但从总体样本的判别成功率上看,考虑含细粒修正后的双曲线模型对场地整体的液化或非液化判别的总成功率是保持稳定的。

    由于本文计算中所采用的砂土细粒含量和黏粒含量的对应关系是通过统计方法得到的,且该计算模型假设对于黏粒含量ρc≤3%的砂土,规范视其为纯净砂。而“Seed简化法”的改进方法“NCEER 推荐法”则规定细粒含量FC≤5%的砂土为纯净砂。而在ρc≤ 3%,FC > 5%的区域,台湾集集地震数据资料存在大量的样本点共110个,占总样本数的40%,其中液化案例53个,非液化案例57个。如果采用建规法和“简化方法”同时对该区域样本点进行判别,会出现是否需要对该区域样本点进行“纯净砂简化”的争议问题,影响所采用的液化判别方法对该区域样本点判别结果的可信度。而本文所建立的修正双曲线模型是基于统计分析的方法得到的,对于黏粒含量ρc≤ 3%和细粒含量FC≤5%的界定并非“洁净砂”的定义,只作为修正双曲线模型中考虑黏粒含量或细粒含量影响的统计边界。本文的修正双曲线模型如果按砂土黏粒含量修正系数考虑该区域的样本点时,修正系数α=1.09,此时修正双曲线模型的液化判别成功率为91%(48/53),而非液化判别成功率为91%(52/57);如果按砂土细粒含量修正系数考虑该区域的样本点时,修正双曲线模型的液化判别成功率为85%(45/53),非液化判别成功率为95%(54/57)。两种情况的判别结果基本一致。由此表明了本文提出的修正双曲线模型不受两种细粒指标关于“洁净砂”定义的限制,并不会因为选用指标不同得出较大差异的判别结果。

    本文基于砂土液化判别的原始双曲线模型,根据大陆地震液化资料和台湾集集地震的标准贯入试验数据资料,建立了适用于含细粒砂性土的液化判别双曲线模型,提出了关于细粒含量影响项的修正系数,并对修正双曲线模型的有效性和适用性进行了检验,主要结论如下:

    (1)从修正的砂土液化判别双曲线模型对于中国大陆地震液化资料的回判结果表明,考虑细粒含量影响的修正双曲线模型能有效提高液化的判别成功率,且能适应不同地震烈度区土层深度在20 m内的判别情况。

    (2)将液化判别双曲线模型进行含细粒影响项的修正弥补了原始双曲线模型只适用于洁净砂的局限性,对于含细粒的砂性土的液化判别结果表现较为稳定。

    (3)本文提出了利用地震液化砂土中细粒含量和黏粒含量的对应关系建立修正系数的方法,该方法可以有效地消除砂土中细粒含量和黏粒含量应用过程中的转换误差。

  • 图  1   三星堆祭祀区考古遗址平面图

    Figure  1.   Archeological excavation site in sacrifice pits area Ruins site of Sanxingdui

    图  2   裂隙形态

    Figure  2.   Morphologies of fractures

    图  3   病害监测区和测点布置

    Figure  3.   Monitoring area and points for desiccation cracking

    图  4   裂隙图像处理流程

    Figure  4.   Process of fracture image processing

    图  5   裂隙宽度和深度测量方法示意图

    Figure  5.   Measuring method for width and depth of fracture

    图  6   测区3裂隙网络正射影像变化对比

    Figure  6.   Orthophoto comparison of fracture network change in measuring area 3

    图  7   测区3裂隙网络数字图像随时间变化

    Figure  7.   Digital image of development progress of fracture network in measuring area 3

    图  8   裂隙网络特征指标及其增加率随时间变化

    Figure  8.   Change of characteristic indexes of fracture network and their increase rates with times

    图  9   裂隙宽度和深度增加量随时间变化

    Figure  9.   Characteristics of increment of width and depth of fractures with time

    图  10   裂隙宽度和深度空间分布特点

    Figure  10.   Spatial distribution characteristics of width and depth of fractures

    图  11   裂隙宽度和深度之间关系

    Figure  11.   Correlation between width and depth of fractures

    图  12   赋存环境相关指标随时间变化

    Figure  12.   Characteristics of related environment factors with time

    图  13   主干裂隙尺寸与赋存环境Pearson相关系数

    Figure  13.   Pearson correlation coefficients of main fractures and different environment factors

    图  14   干缩裂隙尺寸变化的作用方式示意图

    Figure  14.   Variation in size of desiccation cracking at various depths

    图  15   第一次裂隙灌浆试验

    Figure  15.   First grouting experiment

    图  16   第二次裂隙灌浆试验

    Figure  16.   Second grouting experiment

  • [1] 王旭东. 潮湿环境土遗址保护理念探索与保护技术展望[J]. 敦煌研究, 2013(1): 1-6, 125. doi: 10.3969/j.issn.1000-4106.2013.01.001

    WANG Xudong. Exploration of conservation philosophy for earthen sites in humid environments and an outlook on future conservation technology[J]. Dunhuang Research, 2013(1): 1-6, 125. (in Chinese) doi: 10.3969/j.issn.1000-4106.2013.01.001

    [2]

    YAO X, ZHAO F. A comprehensive evaluation of the development degree and internal impact factors of desiccation cracking in the Sanxingdui archaeological site[J]. Heritage Science, 2023, 11(1): 93. doi: 10.1186/s40494-023-00934-3

    [3] 孙满利, 陈彦榕, 沈云霞. 土遗址病害研究新进展与展望[J]. 敦煌研究, 2022(2): 136-148.

    SUN Manli, CHEN Yanrong, SHEN Yunxia. New progress and prospects in research on earthen site deterioration[J]. Dunhuang Research, 2022(2): 136-148. (in Chinese)

    [4] 曲瑾, 马建林, 杨柏. 三星堆城墙干缩裂缝开裂与扩展机理[J]. 工程地质学报, 2020, 28(3): 610-618.

    QU Jin, MA Jianlin, YANG Bai. Crack initiation and propagation mechanism of earth wall at sanxingdui city[J]. Journal of Engineering Geology, 2020, 28(3): 610-618. (in Chinese)

    [5] 孙满利, 张景科. 文物保护学的理论探讨[J]. 西北大学学报(自然科学版), 2022, 52(2): 192-198.

    SUN Manli, ZHANG Jingke. Theoretical discussion on conservation of cultural heritage[J]. Journal of Northwest University (Natural Sciences), 2022, 52(2): 192-198. (in Chinese)

    [6] 唐朝生, 施斌, 崔玉军. 土体干缩裂隙的形成发育过程及机理[J]. 岩土工程学报, 2018, 40(8): 1415-1423. doi: 10.11779/CJGE201808006

    TANG Chaosheng, SHI Bin, CUI Yujun. Behaviors and mechanisms of desiccation cracking of soils[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(8): 1415-1423. (in Chinese) doi: 10.11779/CJGE201808006

    [7] 赵贵刚, 黄英, 张浚枫, 等. 干湿循环作用下云南红土裂缝发展研究[J]. 水土保持学报, 2017, 31(2): 157-165.

    ZHAO Guigang, HUANG Ying, ZHANG Junfeng, et al. Investigation on the development of cracks of laterite of Yunnan under wetting-drying cycles[J]. Journal of Soil and Water Conservation, 2017, 31(2): 157-165. (in Chinese)

    [8] 王峰, 原俊红, 吴图那胜. 玄武岩纤维对黏土干缩开裂特征的影响[J]. 岩土工程学报, 2023, 45(增刊 1): 128-131. doi: 10.11779/CJGE2023S10003

    WANG Feng, YUAN Junhong, WU Tunasheng. Influences of basalt fibers on characteristics of shrinkage cracking of clay[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(S1): 128-131. (in Chinese) doi: 10.11779/CJGE2023S10003

    [9]

    LI J H, ZHANG L M. Study of desiccation crack initiation and development at ground surface[J]. Engineering Geology, 2011, 123(4): 347-358. doi: 10.1016/j.enggeo.2011.09.015

    [10] 赵凡, 姚雪, 胡芮, 等. 三星堆祭祀区地表干缩开裂病害程度与土的性质的关联性[J/OL]. 土木与环境工程学报(中英文): 1-12. DOI: 10.11835/j.issn.2096-6717.2023.094.

    ZHAO Fan, YAO Xue, HU Rui, et al. Correlation between desiccation cracking and soil properties in Sanxingdui sacrifice archeology site[J/OL]. Journal of Civil and Environmental Engineering: 1-12. DOI: 10.11835/j.issn.2096-6717.2023.094. (in Chinese)

    [11] 唐朝生, 王德银, 施斌, 等. 土体干缩裂隙网络定量分析[J]. 岩土工程学报, 2013, 35(12): 2298-2305. http://cge.nhri.cn/article/id/15610

    TANG Chaosheng, WANG Deyin, SHI Bin, et al. Quantitative analysis of soil desiccation crack network[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(12): 2298-2305. (in Chinese) http://cge.nhri.cn/article/id/15610

    [12]

    TANG C S, SHI B, LIU C, et al. Experimental characterization of shrinkage and desiccation cracking in thin clay layer[J]. Applied Clay Science, 2011, 52(1/2): 69-77.

    [13] 林銮, 唐朝生, 程青, 等. 基于数字图像相关技术的土体干缩开裂过程研究[J]. 岩土工程学报, 2019, 41(7): 1311-1318. doi: 10.11779/CJGE201907016

    LIN Luan, TANG Chaosheng, CHENG Qing, et al. Desiccation cracking bebavior of soils based on digital image correlation technique[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(7): 1311-1318. (in Chinese) doi: 10.11779/CJGE201907016

    [14]

    DASOG G S, SHASHIDHARA G B. Dimension and volume of cracks in a vertisol under different crop covers[J]. Soil Science, 1993, 156(6): 424-428. doi: 10.1097/00010694-199312000-00007

    [15]

    XU J, DING T Y, WU D, et al. Extracting the 3D spatial structure of soil cracks in the North China Plain from paraffin casting and exploring the developmental patterns of vertical cross-sectional morphology of the cracks[J]. Geoderma, 2023, 436: 116554. doi: 10.1016/j.geoderma.2023.116554

    [16]

    SANCHEZ M, ATIQUE A, KIM S, et al. Exploring desiccation cracks in soils using a 2D profile laser device[J]. Acta Geotechnica, 2013, 8(6): 583-596. doi: 10.1007/s11440-013-0272-1

    [17]

    CHENG Q, TANG C S, ZHU C, et al. Drying-induced soil shrinkage and desiccation cracking monitoring with distributed optical fiber sensing technique[J]. Bulletin of Engineering Geology and the Environment, 2020, 79(8): 3959-3970. doi: 10.1007/s10064-020-01809-8

    [18]

    TANG C S, WANG D Y, ZHU C, et al. Characterizing drying-induced clayey soil desiccation cracking process using electrical resistivity method[J]. Applied Clay Science, 2018, 152: 101-112.

    [19] 潘斌, 曾召田, 莫红艳, 等. 胀缩性土收缩特性的温度效应试验[J]. 岩土工程学报, 2022, 44(增刊 1): 115-120. doi: 10.11779/CJGE2022S1021

    PAN Bin, ZENG Zhaotian, MO Hongyan, et al. Temperature effects on shrinkage properties of swell-shrink soils[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(S1): 115-120. (in Chinese) doi: 10.11779/CJGE2022S1021

    [20] 唐朝生, 崔玉军, TANG A M, 等. 土体干燥过程中的体积收缩变形特征[J]. 岩土工程学报, 2011, 33(8): 1271-1279. http://cge.nhri.cn/article/id/14163

    TANG Chaosheng, CUI Yujun, TANG A M, et al. Volumetric shrinkage characteristics of soil during drying[J]. Chinese Journal of Geotechnical Engineering, 2011, 33(8): 1271-1279. (in Chinese) http://cge.nhri.cn/article/id/14163

    [21]

    BRAUDEAU E, COSTANTINI J M, BELLIER G, et al. New device and method for soil shrinkage curve measurement and characterization[J]. Soil Science Society of America Journal, 1999, 63(3): 525-535. doi: 10.2136/sssaj1999.03615995006300030015x

    [22]

    BOIVIN P, GARNIER P, TESSIER D. Relationship between clay content, clay type, and shrinkage properties of soil samples[J]. Soil Science Society of America Journal, 2004, 68(4): 1145-1153.

    [23] 杜长城, 祝艳波, 苗帅升, 等. 三趾马红土失水收缩裂缝演化规律研究[J]. 岩土力学, 2019, 40(8): 3019-3027, 3036.

    DU Changcheng, ZHU Yanbo, MIAO Shuaisheng, et al. The evolution of cracks in the dewatering shrinkage process of hipparion red soil[J]. Rock and Soil Mechanics, 2019, 40(8): 3019-3027, 3036. (in Chinese)

  • 期刊类型引用(15)

    1. 李强强,韩赟,史贞,郭永强. 二维周期排桩隔振技术应用于兰州中川国际机场的初探研究. 噪声与振动控制. 2024(02): 249-255 . 百度学术
    2. 余虔,曹正龙,寇璟媛,高妍婕,马新岩. 不同区域粉土累积塑性应变预测模型. 公路交通科技. 2024(07): 56-64 . 百度学术
    3. 童俊,刘勤峰,谭聆言,马胜浩,王旭宾,钱建固. 飞机与地铁动载下道基动力响应现场监测及分析. 结构工程师. 2024(05): 122-127 . 百度学术
    4. 李飞龙,姜昌山,蔡国庆,余虔,韩进宝,张合青. 飞机滑行荷载对水泥混凝土道面及下穿通道的动力响应影响. 土木工程学报. 2024(S2): 80-87 . 百度学术
    5. 董倩,张献民,包伊婷,程少锋,张宇辉. 飞机-道面耦合作用下飞机随机动荷载研究. 振动与冲击. 2023(05): 287-294+312 . 百度学术
    6. 雷宇,刘希重,宣明敏,余虔,叶新宇,张升. 基于服役需求的机场粉土道基临界动应力研究. 铁道科学与工程学报. 2023(03): 950-960 . 百度学术
    7. 罗其奇,张升,李强,叶新宇,张兴胜,余虔. 飞机荷载作用下湿化对粉土道基应力响应的影响. 工程科学与技术. 2023(03): 87-99 . 百度学术
    8. 易文妮,刘津丞,余虔,宣明敏,刘希重,叶新宇,张升. 循环荷载作用下非饱和盐化粉土动力特性. 哈尔滨工业大学学报. 2023(06): 125-133 . 百度学术
    9. 黄博,王宇,盛文军,何淳健,凌道盛. 飞机制动滑行作用下跑道的动力响应. 中南大学学报(自然科学版). 2022(08): 3052-3061 . 百度学术
    10. 冯多,徐林荣,蔡雨,苏娜. 移动荷载作用下横观各向同性层状地基–薄板结构动力响应半解析研究. 岩土工程学报. 2021(05): 858-866 . 本站查看
    11. 罗其奇,张升,叶新宇,李强,马新岩,张兴中. 飞机荷载作用下湿化粉土道基变形特性研究. 中南大学学报(自然科学版). 2021(07): 2188-2199 . 百度学术
    12. 魏川尧,黄忠凯,陈颂,杜一鸣. 下穿通道建设对滑行道变形影响的数值模拟. 佳木斯大学学报(自然科学版). 2021(06): 1-5+14 . 百度学术
    13. 田小芳. 地下结构下穿机场设计中飞机荷载计算研究. 现代隧道技术. 2019(S2): 532-537 . 百度学术
    14. 蔡学彬. 机场复合道面关键力学响应研究. 建筑技术开发. 2018(06): 1-4 . 百度学术
    15. 范鹏贤,王贾博,王德荣. 人工填筑岛礁机场的中长期沉降问题. 防护工程. 2018(04): 70-78 . 百度学术

    其他类型引用(10)

图(16)
计量
  • 文章访问数:  194
  • HTML全文浏览量:  18
  • PDF下载量:  44
  • 被引次数: 25
出版历程
  • 收稿日期:  2023-12-27
  • 网络出版日期:  2024-08-20
  • 刊出日期:  2025-04-30

目录

/

返回文章
返回