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

考虑蠕变特性的高应力软岩隧道变形预测方法与实践

郭新新, 汪波, 王振宇, 于家武

郭新新, 汪波, 王振宇, 于家武. 考虑蠕变特性的高应力软岩隧道变形预测方法与实践[J]. 岩土工程学报, 2023, 45(3): 652-660. DOI: 10.11779/CJGE20220058
引用本文: 郭新新, 汪波, 王振宇, 于家武. 考虑蠕变特性的高应力软岩隧道变形预测方法与实践[J]. 岩土工程学报, 2023, 45(3): 652-660. DOI: 10.11779/CJGE20220058
GUO Xinxin, WANG Bo, WANG Zhenyu, YU Jiawu. Methods and practices for deformation prediction in high-stress soft rock tunnels considering creep characteristics[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(3): 652-660. DOI: 10.11779/CJGE20220058
Citation: GUO Xinxin, WANG Bo, WANG Zhenyu, YU Jiawu. Methods and practices for deformation prediction in high-stress soft rock tunnels considering creep characteristics[J]. Chinese Journal of Geotechnical Engineering, 2023, 45(3): 652-660. DOI: 10.11779/CJGE20220058

考虑蠕变特性的高应力软岩隧道变形预测方法与实践  English Version

基金项目: 

国家自然科学基金项目 U2034205

甘肃省科技计划项目 19ZD2GA005

中铁隧道集团二处有限公司科技研发计划项目 2020-05

详细信息
    作者简介:

    郭新新(1990—),男,博士,副研究员,主要从事隧道工程方面的教学和科研工作。E-mail: zj_gxinxin@163.com

    通讯作者:

    汪波, E-mail: ahbowang@163.com

  • 中图分类号: TU431;U45

Methods and practices for deformation prediction in high-stress soft rock tunnels considering creep characteristics

  • 摘要: 为研究高应力软岩蠕变特性对隧道围岩变形预测的影响,以木寨岭公路隧道为依托,首先采用三维计算模型与多元线性回归相结合的方法分析初始地应力场,并结合围岩段落划分,选择典型计算断面;其次,提出基于[BQ]值的围岩参数取值方法,确定典型计算断面的围岩参数;而后,开展基于M-C模型和Cvsic模型的断面变形计算,剖析岩体蠕变特性对围岩变形的影响;最终,对比了预测结果与实际围岩变形。结果表明:①岩体蠕变特性对围岩变形具有明显增大效应,围岩位移增长量与横断面平均主应力呈正相关;②围岩条件越差,蠕变增大效应越显著;横断面平均主应力越大,蠕变增大效应中位移增长量越大,而位移增长率变化不明显;③蠕变特性对围岩变形等级预测有明显影响,M-C模型预测结果弱于Cvisc模型,与实际围岩变形情况存在较大差异。研究结果为在高应力软岩隧道变形预测中引入岩体蠕变效应奠定了实践基础。
    Abstract: To study the influences of high-stress soft rock creep characteristics on deformation prediction of surrounding rock of tunnels, the Muzhailing highway tunnel is taken as the support project. Firstly, the initial ground stress field is analyzed using a three-dimensional simulation model combined with multiple linear regression. The typical calculation sections are selected in conjunction with the section division of the surrounding rock. Secondly, the relevant method based on the [BQ] value is proposed to determine the parameters of the surrounding rock in the typical calculation section. Then, the section deformation is calculated based on the M-C model and the Cvsic model to analyze the influences of the creep characteristics of rock on the deformation of the surrounding rock. Finally, the predicted results are compared with the actual ones. The results show that: (1) The creep characteristics of rock mass have significant increase effects on the deformation, and the displacement growth rate is positively correlated with the average principal stress of cross-section. (2) The worse the surrounding rock condition, the more pronounced the creeping increase effects. The greater the average principal stress of the cross-section, the larger the displacement growth rate in the creep increase effects, while the change in the displacement growth rate is less noticeable. (3) The creep characteristics have significant effects on the prediction of the deformation level of the surrounding rock. The predicted results by the M-C model are weaker than those by the Cvisc model, which are pretty different from the actual deformations of the surrounding rock. The results lay a practical foundation for introducing the creep effects into the deformation prediction of high-stress soft rock tunnels.
  • 非饱和土在世界各地的工程实践中均有涉及,土体的抗剪强度随含水率或吸力的变化对岩土结构的稳定性至关重要[1-3]。基质吸力与含水率或者饱和度之间的关系被称之为土水特征曲线(soil-water characteristic curve, SWCC)或者土体持水曲线(soil-water retention curve, SWRC)。SWRC是体现非饱和土力学本构关系特征的重要响应曲线,其对土体的非饱和抗剪强度有显著影响[4-7]

    不同非饱和土的微细观测试分析表明[8],SWRC受到孔径分布特征的影响,大致可区分为单峰SWRC和双峰SWRC。Cai等[9]针对桂林红黏土联合轴平移技术、滤纸法和盐溶液蒸汽平衡法进行了宽广吸力范围内的持水特性试验,结果表明压实样呈现双孔隙结构,并具有双峰SWRC,而预固结试样呈现单一孔隙和单峰SWRC。Satyanaga等[10]针对砂土和高岭土的混合物进行了一系列水力特性的测定试验,指出具有双孔隙结构的土体不一定呈现双峰SWRC,但双峰SWRC出现的必要条件是具备双孔隙结构。已有数据表明间断级配土体、粗砂-砾石混合物以及压实黏土等易呈现双峰SWRC。

    过往的大量试验研究了抗剪强度随吸力的演化规律,其中土样普遍为具有单一孔径结构或单峰SWRC的土体。如张俊然等[11]利用轴平移技术和盐溶液蒸汽平衡法对吸力范围在38 MPa内的南阳膨胀土进行了一系列常吸力三轴剪切试验,分析了峰值抗剪强度的演化规律并提出了预测模型;张季如等[12]对不同应力路径下钙质砂抗剪强度进行了试验研究与机理分析;徐筱等[13]对高吸力下黏性土的抗剪强度进行了三轴试验研究;Gao等[14-15]研究了不同制样方法和不同初始干密度对非饱和土应力应变关系的影响。试验结果和理论分析表明吸力和强度之间往往呈现非线性关系[4, 16-17]。对于粗粒土,利用SWRC将吸力乘以饱和度,孔隙应力扩为宏观应力时,强度包络线呈现直线特征,这是因为此时的宏观应力是有效应力的一部分。对于细粒土,毛细部分对剪切强度起主导作用[6-7, 16, 18-19],故而吸力乘以毛细饱和度才是骨架应力,或者有效应力的一部分。已有大量研究表明,结合上述有效应力,极限剪应力的临界状态线可以由唯一的强度包络线表征,即存在一组适用于饱和与非饱和状态的相同强度参数[5-7, 20]。以上结论仅在具有单峰SWRC的土体中得到过验证。

    近期对于具有双孔隙结构,特别是具有双峰SWRC土体力学特性的相关试验在陆续开展。Zhao等[21]对间断级配粗粒土进行了低吸力段的三轴试验,并将其应力应变行为分为4个阶段;Ng等[22]对具有双峰SWRC的压实黄土试样进行了一系列控制吸力的直剪试验。但上述土体的吸力-强度特性和具有单峰SWRC的土体能否统一还未有相关研究。

    本文研究对象为同时具有双孔隙结构及双峰SWRC的土体(不仅限于压实膨胀土)。联合轴平移技术和盐溶液蒸汽平衡法在宽广吸力范围内对内乡膨胀土压实样进行了一系列的非饱和三轴剪切试验,探讨不同吸力对土体应力-应变行为及抗剪强度的影响特性。采用双峰SWRC模型区别毛细和吸附的影响效应,进而在毛细吸应力空间(毛细饱和度乘以吸力)分析土体的非饱和峰值抗剪强度,并利用双孔隙结构土体水力特性试验成果对峰值强度分析方法作评价与验证。

    本试验选用膨胀土取自河南内乡县张岗村以西200 m处,是内乡220 KV变电站拟选站址,取土深度为地下5 m范围,天然干密度为1.69 g/cm3,天然含水率为22.74%。经过X射线衍射分析,内乡膨胀土的主要矿物成分为石英、黏土矿物,其次为斜长石和微斜长石。砂粒的含量为9.97%,粉粒含量为56.77%,黏粒含量为33.26%。粒径大于0.075 mm的颗粒质量为总质量的9.97%,小于50%,按照土的工程分类标准可判断该土样为细粒土,其它物理特性指标可见表 1。内乡膨胀土颗粒分布曲线如图 1所示,根据颗粒分布曲线及塑性指数,可知其为低液限黏土(CL)。

    表  1  内乡膨胀土基本物理指标
    Table  1.  Index properties of Neixiang expansive soil
    颗粒相对质量密度 液限/
    %
    塑限/
    %
    最大干密度/
    (g·cm-3)
    最优含水率/% 自由膨胀率/%
    2.74 45.1 25.3 1.70 21.2 59
    下载: 导出CSV 
    | 显示表格
    图  1  内乡膨胀土的颗粒分布曲线
    Figure  1.  Grain-size distribution curve of Neixiang expansive soil

    较早开展的压汞试验表明[23],内乡膨胀土的重塑预固结试样具有单峰孔径分布;与之不同的是,重塑压实样的微观结构为双孔隙结构,由集聚体间的大孔隙和集聚体内的小孔隙组成(如图 2(a))。联合轴平移技术(0~200 kPa)、滤纸法(200~8 MPa)和饱和盐溶液蒸汽平衡法(3.29~198.14 MPa)进行了宽广吸力范围的持水特性试验[23]。结果显示,预固结样的持水曲线呈现单峰形态。压实样具有双峰持水曲线,即由两段斜率陡峭的曲线组成(双“S”分布),呈对数坐标下的两段半正态分布形式,中间有明显过渡区或者平台段。同时曲线上包括两个进气值AEV1和AEV2(分别对应空气进入到集聚体间和集聚体内孔隙时的孔隙气压力)以及两个残余饱和度R1R2(如图 2(b))。本文进行三轴试验的试样为内乡膨胀土压实样,具有双孔隙结构和双峰持水曲线。

    图  2  重塑压实样和重塑预固结样的孔径分布密度和持水曲线
    Figure  2.  PSDs and SWRCs of compacted and reconstituted specimens

    将自然风干土通过碎土机打碎,过2 mm筛,配制目标含水率为24%的湿土,放入密闭容器中静置两周待水分均匀后取用。试验采用击实法在三开模中制样,分为6层击实以保证均匀性,每一层通过175 g的落锤从300 mm的高度落下。试样直径38 mm,高76 mm,并控制其初始干密度为1.40 g/cm3

    三轴试验采用HKUST双压力室的GDS UNSAT非饱和土三轴仪进行(如图 3),其压力室的底座嵌着进气值为500 kPa的陶土板。差压传感器(DPT)和HKUST双压力室(double-cell)系统量测非饱和土试样的体积变化。试验前进行标定,以消除非饱和土三轴试验系统中体变量测的系统误差。

    图  3  GDS双压力室非饱和土三轴试验仪
    Figure  3.  GDS double-cell triaxial apparatus for testing unsaturated soils with HKUST

    为了实现全吸力范围抗剪强度的测定,联合采用轴平移技术法(0~500 kPa)和饱和盐溶液蒸汽平衡法(3.29~198.14 MPa)进行宽广范围吸力的控制。低吸力试验中,利用饱和陶土板的过水不过气的功能控制低吸力,同时轴平移技术避免孔隙水气化。应力路径如图 4所示,从初始状态开始,首先施加σ3n=20 kPa,s=100,200 kPa进行吸力平衡(ABi),吸力平衡之后施加σ3n=200 kPa进行固结平衡(BiCi),变形稳定之后在恒定围压和吸力值下进行控制应变速率的三轴排水剪切(CiDi)。

    图  4  宽广吸力范围内三轴试验应力路径示意图
    Figure  4.  Stress paths for triaxial tests over a wide suction range

    高吸力试验中,利用饱和盐溶液来控制密封干燥皿内空气的相对湿度[13, 24-26],从而控制土体的总吸力,不同的饱和盐溶液对应的吸力值见表 2[9, 25-26],在不考虑渗透吸力的情况下,总吸力近似等于基质吸力。依据每小时体积变化量小于0.01 cm3的平衡标准[11, 13-14],该过程耗时3~4个月。吸力平衡后把试样装在非饱和土三轴试验仪上,进行常含水率下的固结平衡和三轴剪切,此时可认为是等吸力试验[11, 13, 26]。试样顶部通过孔隙气压控制管路和大气连通,试验中保证孔隙气压为0。采用不透水塑料薄膜将陶土板和试样隔开,保证含水率不变[11, 13]。应力路径如图 4所示(仅展示控制吸力为3.29,21.82 MPa的路径),从初始状态开始,首先利用饱和盐溶液进行高吸力平衡(AXi),吸力平衡之后施加围压σ3n=200 kPa进行固结平衡(XiYi),变形稳定之后在恒定围压和常含水率下进行控制应变速率的不排水(排气)剪切(Yi Zi)。

    表  2  饱和盐溶液及对应的吸力值
    Table  2.  Saturated salt solutions and corresponding suctions
    饱和盐溶液 RH/% 总吸力/MPa
    K2SO4 97.6 3.29
    Na2SO3·7H2O 90.8 13.10
    KCl 85.1 21.82
    NaCl 75.5 38.00
    NaBr 59.1 71.20
    K2CO3 43.2 113.50
    CH3COOK 23.1 198.14
    下载: 导出CSV 
    | 显示表格

    吸力平衡和固结平衡的标准根据进每2 h出水量小于0.01 cm3和每2 h体积变化量小于0.02 cm3确定。采取轴向应变控制三轴剪切加载模式。相关学者对宽广吸力范围内不同三轴剪切速率对峰值强度和体变的影响做了相关研究[11, 27-28],本试验选取轴向剪切速度约为0.00192 mm/min,保证剪切中孔隙水压完全消散。

    本文对重塑压实膨胀土样进行了不同吸力水平下的9组试验,图 5显示不同常吸力下,三轴剪切过程中的应力应变关系和体变结果(压缩为正,剪胀为负)。结果表明,低吸力条件下(s=100 MPa,s=200 kPa)以及中吸力条件下(s=3.29 MPa,s=13.1 MPa),应力应变曲线呈现硬化现象,同时有明显的剪缩行为。随着吸力增加,剪缩体变逐渐减小。

    图  5  三轴剪切过程中应力应变关系
    Figure  5.  Relationship between deviator stress, volumetric strain and axial strain from triaxial shear tests

    高吸力条件下(s为21.82,38,71.2,113.5,198.14 MPa),应力应变曲线呈现峰值后软化现象,试样开始在经历了1%~3%的体积收缩变形之后开始出现明显剪胀。在高吸力状态下(s=113.5 MPa,s=198.14 MPa),试样在较小应变的情况下很快达到峰值强度,随后软化至大应变时偏应力急剧降低。

    随着吸力的增加,剪胀性逐渐增加(或剪缩性减弱)。应力应变曲线的峰后软化现象及体积剪胀现象可以归因于试样失水后孔隙比减小明显(特别是大孔隙收缩[27])或过度固结的应力历史[11]

    图 5中应力应变关系存在峰值,则峰值抗剪强度取偏应力最大值;若应力应变关系不存在峰值,则峰值抗剪强度取轴向应变为15%处的偏应力作为峰值抗剪强度(如图 6)。在吸力s小于21.82 MPa时,峰值和残余偏应力均增加,特别是在吸力s介于200 kPa至21.82 MPa之间增幅较大。在吸力s大于21.82 MPa时,峰值偏应力继续增加,残余偏应力开始下降。这是由于在恒定围压条件下,超过峰值剪应力的连续剪切形成的剪切带严重破坏了土体颗粒接触周围的空气-水弯液面系统,因此吸力对峰值强度的贡献往往大于残余状态。

    图  6  不同吸力下内乡膨胀土峰值强度和残余强度
    Figure  6.  Peak and residual deviator stresses at different suctions of Neixiang expansive soil

    除了外部载荷外,吸力的增加往往会在土体颗粒的接触点上施加表面张力以及吸附作用,从而在土体颗粒之间产生额外的黏结,使土壤结构变得更加坚硬,在剪切过程中表现出脆性。以脆性指数评估土体脆性大小的方法如下[29-30]

    IB=qpeakqrqr
    (1)

    式中:qpeak为峰值抗剪强度;qr为残余抗剪强度。

    图 67可以清楚地看到试件的脆性随吸力的增加而增长,表明剪切破坏后破坏区域逐渐过渡,即在中低吸力下呈桶形或中心鼓形的延性破坏模式,在高吸力下发生应变局部化现象,并产生明显的剪切破坏面。此外,脆性增加了峰后软化的幅度,表现为图 6中峰值偏应力与残余破坏状态之间的差值增大。

    图  7  不同吸力下试样剪切后状态
    Figure  7.  Typical failures of specimens at different suction states

    Bishop[31]提出了经典的非饱和土有效应力公式:

    σij=σijuaδij+χsδij
    (2)

    式中:σij为总应力张量;ua为孔隙气压力;χ为有效应力参数,s为基质吸力;δij为Kronecker符号。

    式(2)的有效应力参数χ最初采用饱和度Sr加以量度,这种做法由于高估了吸力对骨架应力的作用,忽视了吸附吸力的存在[5-7, 16, 18-20],导致预测的峰值抗剪强度往往大于实际强度,特别是对于细粒土以及低饱和度的情况。为此,近来一些研究发现,可将总饱和度Sr的量度替换成毛细饱和度Scapr[5-7, 20, 32-33]以区分毛细和吸附不同作用的影响:

    σij=σijuaδij+Scaprsδij
    (3)

    研究表明,结合上述有效应力,极限剪应力的临界峰值强度状态线可以由唯一的强度包络线表征[5-7, 20](如图 8(a),即存在一组适用于饱和与非饱和状态的相同强度参数。唯一的强度包络线不仅在低吸力适用,也在高吸力范围内得到了验证,同时对于粗粒土和细粒土显示具有线性特征[34]

    图  8  Madrid黏性砂土和内乡膨胀土峰值剪切强度特征
    Figure  8.  Peak shear strength behaviors for Madrid clayey sand and Neixiang expansive soil

    需要注意的是,以上研究仅针对大量具有单峰持水曲线的土体,而没有涉及具有双峰持水曲线的土体。这里利用笔者[5, 20, 23, 32]提出的可以区分毛细和吸附饱和度的双峰持水曲线模型来计算毛细饱和度:

    Scapr=(1α)A+(ααC(1B)1C(B1))B
    (4)
    A=12erfc(ln(s/sm1)2ζ1)
    (5)
    B=12erfc(ln(s/sm2)2ζ2)
    (6)
    C=β(1lnslnsd)
    (7)

    式中:sm1为集聚体间中值孔径对应的基质吸力(kPa),sm2为集聚体内中值孔径对应的基质吸力(kPa);ζ1为表征集聚体间孔径分布宽度的无量纲参数,ζ2为表征集聚体内孔径分布宽度的无量纲参数,erf( )为误差函数;α=θs2/θs1θs1为第一进气值对应的含水率,即土材料完全饱和时的体积含水率,θs2为第二进气值对应的含水率,即集聚体内的饱和含水率;β为表征吸附作用大小的材料参数;sd为干燥情况下的最大吸力,取sd=106(kPa)。

    参数α代表两种大小孔隙结构相对大小,技术上可以根据压汞试验进行标定。参数β控制了吸附水在最大吸力处的变化速率,可以通过检测吸附水含量的变化确定,技术上可以通过核磁共振技术实现[35]。参数smsm1sm2)代表中等孔径对应的吸力值,ζζ1ζ2)是对数坐标下孔径的标准差,这些参数可以通过微观试验,例如压汞试验进行测定[36]。以上参数的标定方法都是基于相应的微观试验或者这些参数直接控制的测试结果来确定。然而,在未进行微观试验的状况下,模型参数可以通过对持水数据的拟合获得,即基于非线性拟合(详细的讨论见文献[37])的方法进行,相同的参数标定流程见笔者对双峰持水曲线的建模过程[23]

    图 8(b)表明,利用相同的表征方法利用毛细饱和度将吸力从孔隙应力扩为宏观应力时(参数见表 3,其中αβζ1ζ2为无量纲参数),强度包络线由两条直线段组成,这一强度特征的微细观机理可解释为:由于具有双峰持水曲线的土体常存在集聚体结构,集聚体间大孔隙以毛细水占主导,集聚体内小孔隙赋存吸附水。对于毛细水,来源于空气-水弯液面提供的表面张力,能够在骨架间进行应力传递,完全贡献了第一段直线的强度;吸附作用主要源于3种不同的物理化学机制[38-40],即离子水化、内表面水化及颗粒表面水化,离子水化来源于原子尺度的电子力作用,内表面水化和颗粒表面水化受到范德华力影响。此种物理化学作用虽不能形成骨架应力,往往会产生额外的胶结作用,从而增加土体的第二段抗剪强度[27]

    表  3  双峰持水曲线模型参数
    Table  3.  Parameters of bimodal water-retention curve
    土体 sm1/
    kPa
    sm2/
    kPa
    ζ1 ζ2 β α
    膨胀土 41.21 14500 1.05 0.74 0.95 0.70
    黄土[17, 41] 9.69 2097 1.00 1.24 0.20 0.59
    M3[10] 9.35 121.9 0.24 0.86 0.05 0.85
    M5[10] 14.22 135 0.19 0.77 0.05 0.87
    M6[10] 11.74 131.7 0.10 0.91 0.30 0.85
    M8[10] 0.53 28.12 0.23 1.34 0.38 0.84
    红黏土[42-43] 68.74 9430 1.22 0.74 0.48 0.79
    下载: 导出CSV 
    | 显示表格

    Ng等[22, 27]对具有双峰SWRC的压实黄土试验进行了宽广范围的直剪试验;Satyanaga等[10]对砂土和高岭土不同配比的混合物M3(70%砂土、30%高岭土)、M5(50%砂土、50%高岭土)、M6(40%砂土、60%高岭土)以及M8(20%砂土、80%高岭土)利用轴平移技术进行了三轴剪切试验;孙德安等[42-43]对桂林红黏土进行了一系列土水特性试验和控制吸力的直剪试验,测试土体的基本物理指标见表 4,双峰SWRC的模型参数见表 3。基于土体持水特性和剪切试验的试验结果,在毛细吸应力空间重新整理峰值抗剪强度值(如图 9),结果显示强度随毛细吸应力的变化均由双段直线组成,第二段直线位于第一段直线的上方。这种现象的产生机理可同样归因于:对于具有双孔隙结构的土体,其持水曲线呈现双峰形态,高吸力下的峰值抗剪强度不仅由毛细作用形成的空气-水弯液面提供,还受到吸附作用引起的胶结作用增强。

    表  4  具有双峰持水曲线土体的基本物理指标
    Table  4.  Physical properties of soils with bimodal SWRC
    土体 液限/
    %
    塑限/
    %
    最大干密
    度/(g·cm-3)
    c/kPa φ/(°)
    黄土[17, 41] 36.0 19.0 1.30 0 27.7
    M3[10] 36.8 21.6 1.85 5.00 34.0
    M5[10] 46.7 27.4 1.75 9.00 26.0
    M6[10] 48.9 28.2 1.59 12.00 26.0
    M8[10] 57.8 33.6 1.45 16.00 25.0
    红黏土[42, 43] 77.8 42.1 1.52 1.57 32.1
    下载: 导出CSV 
    | 显示表格
    图  9  具有双孔隙结构和双峰SWRC土体的非饱和峰值强度特性(毛细吸应力-剪应力/偏应力空间)
    Figure  9.  Unsaturated shear strengths of soils with dual-porosity structure and bimodal SWRC

    试验数据和理论研究表明,对于具有双孔隙结构的土体,在集聚体间(低吸力范围),用毛细饱和度乘以吸力将吸力从孔隙应力扩为宏观应力时,能更为合理地量度吸力对有效应力的贡献;在集聚体内(高吸力范围),此时除毛细吸应力贡献峰值抗剪强度之后,还受到结构重组和胶结作用的影响(铁和锰氧化物的胶结[27])。后续研究需要考虑双孔隙结构以及胶结作用对抗剪强度的影响特性。

    本文针对具有双孔隙结构及双峰SWRC的非饱和重塑膨胀土的抗剪强度开展了试验及理论研究。联合轴平移技术和盐溶液蒸汽平衡法对宽广吸力范围内膨胀土压实样进行了一系列的非饱和三轴剪切试验。同时对峰值抗剪强度在毛细吸应力空间(毛细饱和度乘以吸力与偏应力)进行了分析讨论。得到以下4点结论。

    (1)中低吸力下的应力应变关系为应变硬化,伴随着剪缩行为,试样呈现中心鼓形的延性破坏模式;高吸力下为应变软化,伴随着剪胀现象,试样有明显的破坏剪切带。

    (2)脆性随吸力的增加而增长,表现为峰后软化的幅度增大,即峰值抗剪强度与残余抗剪强度之间的差值增大。

    (3)毛细吸应力能更为合理地揭示吸力对剪切强度的贡献,在毛细吸应力-强度空间,双孔隙结构非饱和土(不仅限于压实膨胀土)均呈现双段线性的强度特征。

    (4)低吸力范围内,双孔隙结构非饱和土的抗剪强度由饱和强度与毛细吸应力贡献;高吸力范围内,抗剪强度应由饱和强度、毛细吸应力与胶结作用提供。

  • 图  1   实测应力分布

    Figure  1.   Distribution of measured stress

    图  2   隧址区三维地质模型

    Figure  2.   Three-dimensional geological model

    图  3   模型加载示意图

    Figure  3.   Model loading

    图  4   初始应力场分布规律

    Figure  4.   Distribution laws of initial stress field

    图  5   隧道计算模型图

    Figure  5.   Model for tunnel calculation

    图  6   最大位移

    Figure  6.   Maximum displacement

    图  7   Cvisc模型最大位移增长量、增长率

    Figure  7.   Maximum displacement increments and growth rates by Cvisc model

    图  8   蠕变增大效应与应力场关系

    Figure  8.   Relationship between creep increase effect and stress field

    图  9   不同岩性时,蠕变增大效应与XZ向平均应力关系曲线

    Figure  9.   Relation curve between creep increase effect and XZ average stress under different lithologies

    图  10   不同围岩本构模型变形等级预测结果

    Figure  10.   Predicted results of deformation grade of different constitutive models for surrounding rock

    图  11   典型断面位移

    Figure  11.   Displacements of typical section

    表  1   各岩层的单轴饱和抗压强度

    Table  1   Uniaxial saturated compressive strengths of rock strata

    统计项目 单轴饱和极限抗压强度
    统计指标岩土名称 统计个数n 范围值 算术平均值fm 标准差σf 变异系数δ 修正系数 标准值
    炭质板岩
    (P1
    30 11.23~45.66 28.56 7.99 0.28 0.91 26.0
    炭质板岩
    (C1
    6 16.37~41.24 25.48 8.53 0.33 0.75 19.2
    断层压碎岩 8 12.12~15.94 14.00 1.40 0.10 0.93 13.00
    注:已剔除异常值。
    下载: 导出CSV

    表  2   主要岩层物理力学参数

    Table  2   Main physical and mechanical parameters of rock strata

    岩体类型 变形模量E/GPa 泊松比 重度/(kN·m-3)
    中风化炭质板岩(P1 2.0 0.35 27.0
    断层压碎岩 1.5 0.40 27.0
    中风化灰岩(C1 6.0 0.30 27.0
    中风化炭质板岩(C1 3.0 0.35 27.0
    下载: 导出CSV

    表  3   N2钻孔主应力反演与实测对比

    Table  3   Comparison between inversion and measurement of N2 principal stress

    钻孔编号 测点 SH/MPa Sh/MPa Sv/MPa
    实测值 回归值 绝对误差 相对误差 实测值 回归值 绝对误差 相对误差 实测值 回归值 绝对误差 相对误差
    N2 1 25.7 25.1 0.6 2.3% 16.8 20.5 3.7 22.2% 10.7 12.8 2.1 19.2%
    2 24.1 25.4 1.3 5.7% 17.0 20.9 3.9 23.0% 11.1 13.3 2.2 19.1%
    下载: 导出CSV

    表  4   围岩弹塑性参数

    Table  4   Elastic-plastic parameters of surrounding rock

    岩体类型 界限 变形模量E/GPa 泊松比ν 黏聚力c/MPa 内摩擦角φ/(°) 重度γ/(kN·m-3) 备注
    炭质板岩 上限 2.00 0.35 0.80 28 27.0 [BQ]=214
    下限 1.20 0.38 0.50 25 27.0 [BQ]=54
    断层压碎岩 1.20 0.39 0.45 24 27.0
    下载: 导出CSV

    表  5   木寨岭隧道围岩蠕变参数

    Table  5   Creep parameters of surrounding rock in Muzhailing tunnel

    围岩岩性 Maxwell剪切模量GM/GPa Maxwell黏度ηM/(GPa·h-1) Kelvin剪切模量GK/GPa Kelvin黏度ηK/(GPa·h-1)
    炭质板岩 1.00 8300 2.00 138
    断层压碎岩 0.75 6225 1.50 103.50
    注:灰岩蠕变参数按炭质板岩参数选取。
    下载: 导出CSV

    表  6   木寨岭公路隧道变形预测相关计算资料(左线)

    Table  6   Data of deformation prediction of Muzhailing highway tunnel (left line)

    段落起讫里程
    (ZK)
    主要岩性 [BQ]均值 计算断面1 计算断面2 段落起讫里程
    (ZK)
    主要岩性 [BQ]均值 计算断面1 计算断面2
    212+000—+185 190.5 211+900 212+100 216+005—217+100 171 216+100 217+100
    212+185—+385 54 212+200 212+300 217+100—+300 54.3 217+100 217+300
    212+385—+585 212+400 212+500 217+300—+520 217+300 217+500
    212+585—+785 54 212+600 212+700 217+520—+720 54.3 217+600 217+700
    212+785—213+355 189 212+800 213+200 217+720—219+030 214 217+700 219+000
    213+355—+555 92.8 213+400 213+500 219+030—+230 84.3 219+100 219+200
    213+555—+715 213+600 213+700 219+230—219+610 219+300 219+600
    213+715—+915 92.8 213+800 213+900 219+610—810 84.3 219+700 219+800
    213+915—214+285 190 214+000 214+200 219+810—221+340 203 219+900 221+100
    214+285—+385 114 214+300 221+340—+540 92.5 221+400 221+500
    214+385—+505 214+400 214+500 221+540—221+690 221+600
    214+505—+705 93 214+600 214+700 221+690—+890 92.5 221+700 221+800
    214+705—+830 133 214+800 221+890—222+810 196.8 221+900 222+800
    214+830—215+030 54.3 214+900 215+000 222+810—223+010 92.5 222+900 223+000
    215+030—+150 215+100 223+010—+105 223+100
    215+150—+635 54.3 215+400 215+600 223+105—+305 92.5 223+200 223+300
    215+635—+805 215+700 215+800 223+305—224+000 190.6 223+400 224+000
    215+805—216+005 54.3 215+900 216+000
    注:岩性①为炭质板岩,岩性②为断层压碎岩;里程K212前、及K224后,因垂直应力小,地应力场模拟出现明显偏差,予以剔除。
    下载: 导出CSV
  • [1]

    TIAN Siming, WANG Wei, GONG Jiangfeng. Development and prospect of railway tunnels in China(including statistics of railway tunnels in China by the end of 2020)[J]. Tunnel Construction, 2021, 41(2): 308-325.

    [2]

    LI Zhijun, GUO Xinxin, MA Zhenwang, et al. Research status and high-strength pre-stressed primary (type) support system for tunnels with large deformation under squeezing conditions[J]. Tunnel Construction, 2020, 40(6): 755-782.

    [3]

    SINGH B, JETHWA J L, DUBE A K, et al. Correlation between observed support pressure and rock mass quality[J]. Tunnelling and Underground Space Technology, 1992, 7(1): 59-74. doi: 10.1016/0886-7798(92)90114-W

    [4]

    GOEL R K, JETHWA J L, PAITHANKAR A G. Tunnelling through the young Himalayas—a case history of the Maneri-Uttarkashi power tunnel[J]. Engineering Geology, 1995, 39(1/2): 31-44.

    [5]

    HOEK E, MARINOS P. Predicting squeeze (Part 2) [J]. Tunnels and Tunnelling International, 2000, 32(12): 33-36.

    [6] 陈卫忠, 田云, 王学海, 等. 基于修正[BQ]值的软岩隧道挤压变形预测[J]. 岩土力学, 2019, 40(8): 3125-3134.

    CHEN Wei-zhong, TIAN Yun, WANG Xue-hai, et al. Squeezing prediction of tunnel in soft rocks based on modified [BQ] [J]. Rock and Soil Mechanics, 2019, 40(8): 3125-3134. (in Chinese)

    [7] 骆顺天, 杨凡杰, 周辉, 等. 基于统计分析的地下厂房边墙最大收敛变形多指标预测方法[J]. 岩土力学, 2020, 41(10): 3415-3424. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202010028.htm

    LUO Shuntian, YANG Fanjie, ZHOU Hui, et al. Multi-index prediction method for maximum convergence deformation of underground powerhouse side wall based on statistical analysis[J]. Rock and Soil Mechanics, 2020, 41(10): 3415-3424. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202010028.htm

    [8] 赵阳升. 岩体力学发展的一些回顾与若干未解之百年问题[J]. 岩石力学与工程学报, 2021, 40(7): 1297-1336. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202107001.htm

    ZHAO Yangsheng. Retrospection on the development of rock mass mechanics and the summary of some unsolved centennial problems[J]. Chinese Journal of Rock Mechanics and Engineering, 2021, 40(7): 1297-1336. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202107001.htm

    [9] 王道远, 刘佳, 张逴, 等. 高地应力深埋隧道断裂破碎带大变形控制方法现场试验研究[J]. 岩土工程学报, 2020, 42(4): 658-666. doi: 10.11779/CJGE202004008

    WANG Daoyuan, LIU Jia, ZHANG Chuo, et al. Field tests on large deformation control method for surrounding rock of deep tunnel in fault zone with high geostress[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(4): 658-666. (in Chinese) doi: 10.11779/CJGE202004008

    [10] 徐国文, 何川, 汪耀, 等. 流变荷载作用下隧道裂损二次衬砌结构安全性能研究[J]. 土木工程学报, 2016, 49(12): 114-123. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201612013.htm

    XU Guowen, HE Chuan, WANG Yao, et al. Study on the safety performance of cracked secondary lining under action of rheological load[J]. China Civil Engineering Journal, 2016, 49(12): 114-123. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC201612013.htm

    [11] 白国权. 隧道二衬严重开裂段病害综合治理研究[J]. 铁道标准设计, 2014, 58(8): 133-138. https://www.cnki.com.cn/Article/CJFDTOTAL-TDBS201408035.htm

    BAI Guoquan. Research on comprehensive treatment for severe cracking of tunnel secondary lining[J]. Railway Standard Design, 2014, 58(8): 133-138. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TDBS201408035.htm

    [12] 孟陆波, 潘皇宋, 李天斌, 等. 鹧鸪山隧道二次衬砌开裂机理及支护时机探讨[J]. 现代隧道技术, 2017, 54(2): 129-136. https://www.cnki.com.cn/Article/CJFDTOTAL-XDSD201702020.htm

    MENG Lubo, PAN Huangsong, LI Tianbin, et al. Secondary lining cracking mechanism and the best time for supporting the zhegushan tunnel[J]. Modern Tunnelling Technology, 2017, 54(2): 129-136. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-XDSD201702020.htm

    [13] 汪波, 李天斌, 何川, 等. 强震区软岩隧道大变形破坏特征及其成因机制分析[J]. 岩石力学与工程学报, 2012, 31(5): 928-936. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201205010.htm

    WANG Bo, LI Tianbin, HE Chuan, et al. Analysis of failure properties and formatting mechanism of soft rock tunnel in meizoseismal areas[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(5): 928-936. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201205010.htm

    [14] 汪波, 何川, 吴德兴, 等. 苍岭特长公路隧道地应力场反演分析[J]. 岩土力学, 2012, 33(2): 628-634. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201202051.htm

    WANG Bo, HE Chuan, WU Dexing, et al. Inverse analysis of in situ stress field of Cangling super-long highway tunnel[J]. Rock and Soil Mechanics, 2012, 33(2): 628-634. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201202051.htm

    [15] 代聪, 何川, 陈子全, 等. 超大埋深特长公路隧道初始地应力场反演分析[J]. 中国公路学报, 2017, 30(10): 100-108. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL201710013.htm

    DAI Cong, HE Chuan, CHEN Ziquan, et al. Inverse analysis of initial ground stress field of deep embedded and extra long highway tunnel[J]. China Journal of Highway and Transport, 2017, 30(10): 100-108. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGGL201710013.htm

    [16] 李建伟, 雷胜友, 李振, 等. 木寨岭隧道炭质板岩流变力学特性研究[J]. 隧道建设, 2012, 32(1): 36-40. https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201201011.htm

    LI Jianwei, LEI Shengyou, LI Zhen, et al. Investigation on rheologic properties of carbonaceous slate in muzhailing tunnel[J]. Tunnel Construction, 2012, 32(1): 36-40. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-JSSD201201011.htm

    [17] 宋勇军, 雷胜友, 邹翀, 等. 干燥与饱水状态下炭质板岩蠕变特性研究[J]. 地下空间与工程学报, 2015, 11(3): 619-625, 664. https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201503015.htm

    SONG Yongjun, LEI Shengyou, ZOU Chong, et al. Study on creep characteristics of carbonaceous slates under dry and saturated states[J]. Chinese Journal of Underground Space and Engineering, 2015, 11(3): 619-625, 664. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201503015.htm

    [18] 陶志刚, 罗森林, 康宏伟, 等. 公路隧道炭质板岩变形规律及蠕变特性研究[J]. 中国矿业大学学报, 2020, 49(5): 898-906. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD202005010.htm

    TAO Zhigang, LUO Senlin, KANG Hongwei, et al. Analysis of deformation law and creep characteristics of carbonaceous slate in highway tunnel[J]. Journal of China University of Mining & Technology, 2020, 49(5): 898-906. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD202005010.htm

    [19] 徐慧臣. 木寨岭深埋隧道板岩吸水强度软化结构效应实验研究[D]. 北京: 中国矿业大学(北京), 2019.

    XU Huichen. Experimental Study on Strength Softening of Slate after Water Absorption Due to Structure Effects for Muzhailing Deep-Buried Tunnel[D]. Beijing: China University of Mining & Technology, Beijing, 2019. (in Chinese)

    [20]

    WANG M N, WANG Z L, TONG J J, et al. Support pressure assessment for deep buried railway tunnels using BQ-index[J]. Journal of Central South University, 2021, 28(1): 247-263.

  • 期刊类型引用(2)

    1. 张村,方尚鑫,贾胜,王永乐,王方田,白庆升. 基于CT扫描的三维重构煤体加载损伤演化特征及尺寸效应. 矿业科学学报. 2024(03): 413-425 . 百度学术
    2. 李文新. 基于非饱和渗流特征的土体抗剪强度衰减预测研究. 矿产勘查. 2024(10): 1914-1920 . 百度学术

    其他类型引用(6)

  • 其他相关附件

图(11)  /  表(6)
计量
  • 文章访问数:  370
  • HTML全文浏览量:  55
  • PDF下载量:  138
  • 被引次数: 8
出版历程
  • 收稿日期:  2022-01-12
  • 网络出版日期:  2023-03-15
  • 刊出日期:  2023-02-28

目录

/

返回文章
返回