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

频发微小地震下顺层岩质边坡累积损伤及稳定性分析

刘新荣, 许彬, 刘永权, 王继文, 蔺广义

刘新荣, 许彬, 刘永权, 王继文, 蔺广义. 频发微小地震下顺层岩质边坡累积损伤及稳定性分析[J]. 岩土工程学报, 2020, 42(4): 632-641. DOI: 10.11779/CJGE202004005
引用本文: 刘新荣, 许彬, 刘永权, 王继文, 蔺广义. 频发微小地震下顺层岩质边坡累积损伤及稳定性分析[J]. 岩土工程学报, 2020, 42(4): 632-641. DOI: 10.11779/CJGE202004005
LIU Xin-rong, XU Bin, LIU Yong-quan, WANG Ji-wen, LIN Guang-yi. Cumulative damage and stability analysis of bedding rock slope under frequent microseisms[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(4): 632-641. DOI: 10.11779/CJGE202004005
Citation: LIU Xin-rong, XU Bin, LIU Yong-quan, WANG Ji-wen, LIN Guang-yi. Cumulative damage and stability analysis of bedding rock slope under frequent microseisms[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(4): 632-641. DOI: 10.11779/CJGE202004005

频发微小地震下顺层岩质边坡累积损伤及稳定性分析  English Version

基金项目: 

国家自然科学基金项目 41972266

国家自然科学基金项目 41772319

重庆市研究生科研创新项目 CYS18026

详细信息
    作者简介:

    刘新荣(1969—),男,教授,博士生导师,主要从事岩土工程的教学与科研工作。E-mail: liuxrong@126.com

    通讯作者:

    许彬, E-mail: geotechnicale2016@163.com

  • 中图分类号: TU432

Cumulative damage and stability analysis of bedding rock slope under frequent microseisms

  • 摘要: 三峡库区蓄水后诱发的高频度微小地震对边坡稳定性产生了重大影响。采用振动台模型试验和UDEC离散元数值计算方法,深入地探究了库区典型顺倾层状岩质边坡在高频次微小地震下的累积损伤和稳定性。研究表明:①地震持续作用下,边坡自振频率、阻尼比、损伤度和损伤速率依次降低、增大、累积变大和不断提升,各测点PGA响应表现为“高程效应”和“趋表效应”,且PGA放大系数均呈降低趋势;②高、低动荷载振幅阶段边坡岩体非线性累积损伤模型可分别用指数函数和三次函数描述,其演化曲线分别呈初期轻微降低、中期线性递增、后期平缓微增的“S”型特征和急速增长的“陡升”型特征;③边坡累积损伤–失稳破坏呈现为起伏体爬坡–啃断–磨平、次级节理(层面)起裂–扩展–贯通、坡体沿复合破坏面发生整体滑移、失稳破坏后岩体以破碎–大型–巨型块状堆积于坡脚,且含起伏体边坡整体稳定性更优;④动荷载振幅、动荷载频率、坡高、坡角增大而层面厚度减小时,边坡临界失稳微震作用次数减小、累积永久位移增大、稳定性系数减小,且层面出露边坡更易发生失稳破坏。
    Abstract: Micro-seism frequently occurs in the Three Gorges Reservoir after impoundment and has a significant impact on slope stability. To investigate the cumulative damage and stability of typical bedding rock slope under frequent microseisms, the shaking table model tests and UDEC discrete element numerical analyses are conducted. The results indicate that: (1) The damping ratio, damage degree and damage rate of slope increase with the increase of seismic duration, while the natural frequency and PGA amplification coefficient decrease. The PGA response of slope shows "elevation effect" and "surface effect". (2) The non-linear cumulative damage model under high and low seismic amplitudes can be represented by exponential and cubic functions respectively, and the evolution curve shows "S" characteristics of slight decline in initial stage, linear increase in intermediate stage and gentle increase in later stage, and "steep" characteristics of rapid rise. (3) The gradual process of cumulative damage and instability is as follows: undulant body climbs, shear and grinds, cracks initiate, propagate and run through, the slope slips wholly along the composite sliding surface, and rock mass accumulates at the slope foot in the form of broken, large and giant blocks. The stability of slope with undulant body is better. (4) The critical times of microseism decrease and the cumulative permanent displacement increases with the increase of seismic amplitude, frequency, slope height and slope angle and the decrease of bedding thickness, and accordingly the stability factor decreases. The slope with outcropped bedding is more prone to fail.
  • 在中国建设地铁的过程中,人工冻结法是针对沿海地区地层进行加固的一种有效工法。具有冻融敏感特性的淤泥质黏土作为沿海地区冻结施工土层,其本身地下水的渗流速度很小,多数情况下无需考虑渗流因素。然而随着沿海城市及滨海地区的地下空间不断向更深处进一步开发,当工程为越江(跨海)隧道[1]、工程开挖周围存在溶洞、工程建设在地下水排泄通道附近[2],或工程开挖含水层的隔水板时,遭遇的渗流就会产生无法及时顺利交圈和不均匀沉降等不利影响[3]

    近几年,在一些特殊工程实践中发现渗流对冻结法施工存在着影响,但研究集中在冻结效果的数值模拟或室内样品试验及施工现场应对措施等方面。赖远明等[4]以温度场和渗流的有限元计算公式、透水介质的渗流体力计算方法、冻胀荷载的估算方法以及耦合方程的求解方法为基础,提出了寒区岩石隧道温度场,渗流和应力场非线性耦合理论。程桦等[5]基于Harlan水热耦合模型,使用数值模拟的方法对饱和砂土在竖井冻结时冻结壁的扩展过程进行了模拟;周晓敏等[6]使用温度,渗流,应力耦合理论,数值模拟了在低流速水平渗流的作用下,立井冻结时的模型。Hu等[7]使用原位冻胀模型数值模拟了砂性土在水平渗流的作用下二维的冻结效果。Vitel等[8]建立了一个完全热力学一致的耦合热-水力数值模型,并模拟了渗流条件下饱和不可变形多孔介质的地层冻结效果。Ahmed等[9]建立了一种水热耦合有限元模型,用以在渗流情况下通过寻找最佳位置来优化隧道砂层冻结的冻结管排布。Pimentel等[10]设计了针对测试高速渗流作用下砂性土冻结法试验装置并进行了冻结效果的模拟试验;崔灏[11]、陈湘生等[12]、于琳琳等[13]也做过类似的模拟试验研究。可以看出针对渗流作用下冻结法的研究大多是针对自身渗透性高的砂性土层。软黏土冻结时受分凝作用影响明显,从而体现出与砂性土完全不同的冻融特性[14-15]。同时渗流问题的存在条件也区别很大。软土地区冻结法面临的渗流问题,往往是冻结软黏土时存在较大渗流边界,或者渗流排泄区。此时渗流不仅对冻结效果产生影响,比如温度场、冻结时间、冻结速率等;还会增大对周围软土环境的不确定性影响,不均匀沉降问题、冻胀力分布问题会更加显著。目前国内外针对软黏土下伏存在较大渗流砂层,渗流对上覆冻融敏感特性软黏土冻结施工及周围软土环境带来影响的可参考性研究成果相对较少。

    本文通过沿海地区典型的上黏下砂的组合地层设计,分析了渗流对软黏土冻结效果及周围环境效应的影响,研究结果可以帮助工程优化冻结方案,为控制不均匀沉降、防止冻胀力过大提供指导;从而提高沿海复杂地质环境下冻结工程的安全性。

    随着沿海城市地下空间进一步扩大化建设,越来越多的大型工程交叉影响,使得人工地层冻结法面临的水文地质条件更具复杂化。图1是当隧道联络通道周围存在大型基坑降水引发的强渗流环境的工况(典型案例为上海复兴东路越江通道冻结工程)。本模型试验针对沿海软土地区典型的地层单元,设计了软黏土下伏边界砂层不同渗流速度下,人工地层冻结法的模型工况,组合地层渗流冻结法模型设计如图1所示,上海市《旁通道冻结法技术规程(DG/TJ 08— 902—2016)》要求渗流速度超过2 m·d-1的工况需要进行特殊处理,故设计冻结法可能遇到的地下水渗流速度设计了3,1.2,0.5 m·d-1 3种渗流速度的模型试验,渗流循环水体在砂土中水平流过。其中3 m·d-1的渗流会致使冻结帷幕无法扩展到设计厚度(DG/TJ 08—902—2016)。根据第1批模型试验结果,在能达到冻结帷幕设计要求的渗流区间(≤1.2 m·d-1)内补充了0.75,0.3,0 m·d-1的3种渗流速度下的试验。

    图  1  工况设计
    Figure  1.  Design of working conditions

    本试验主要研究土体宏观的应力应变特性,在宏观角度针对土体应力应变特性进行相似比设计,而土体材料微观颗粒特性在试验中作为次要问题尽可能弱化,因此模型土体没有根据相似材料的缩尺设计,而是采用相同的土体材料进行重塑土制样。重塑过程严格遵守土工试验方法标准(GB/T 50123—2019),软黏土采用泥浆固结法,而砂土采用分层击实法。

    土是有记忆性的,且软黏土具有触变性的特点,本试验通过排水固结及静置的方式,使重塑模型土排水固结至设定含水率并逐渐恢复强度,尽可能接近原状土(见表1);砂性土使用分层击实法,通过控制含水率和孔隙比进行重塑,强度亦尽量达到未扰动前。这也是目前模型试验都采用重塑土样的基本前提:在尽可能(或部分)恢复土样基础物理指标及强度的同时,又能保证模型土体的均匀性,尽量减少交叉影响,从而更好了揭示研究因素的规律特征。

    表  1  模型试验土体特性表
    Table  1.  Soil properties of model soils
    类型内摩擦角/(°)黏聚力/kPa重度/(kN·m-3)含水率/%孔隙比
    原型(软黏土)11.514.016.9049.861.38
    模型(软黏土)11.012.716.5550.001.38
    原型(砂土)34.1018.9026.900.97
    模型(砂土)27.000.97
    注:ac,φ数据来自正常固结的快剪试验,模型(砂土)通过控制孔隙比进行重塑。
    下载: 导出CSV 
    | 显示表格

    试验首先使用分层击实法对砂土层进行装填。然后装填黏土层。当软黏土层装填至规定高度后,使用千斤顶对整个组合地层进行固结加载。固结应力与现场冻结区域的竖向应力230 kPa相匹配,尽量使得模拟冻结区域的软黏土能够具有原状软黏土围压下的强度。固结预压21 d[16]且每天沉降量不超过0.01 mm时,取下千斤顶,土体回弹达到7 d且每天回弹量不超过0.01 mm时,预压过程结束。固结预压完成后,取模型箱角落处的黏土样进行正常固结的快剪试验,结果显示重塑软黏土基本恢复冻结区软黏土的物理指标与强度。

    模型试验基本物理量量纲及常规相似常数见表2

    表  2  模型试验相似特性表
    Table  2.  Similar characteristics of model tests
    物理量量纲相似常数
    长度lL1/N
    应力σML-1T-21/N
    位移lL1/N
    时间tT1/N
    应变e 1
    下载: 导出CSV 
    | 显示表格

    在冻结法施工的工况下,地下水流速较小,符合达西定律:

    v=Ki, (1)

    式中,v为渗流速度(m·d-1),K为渗透系数(m·d-1),i为水力梯度。当原型和模型使用同一种土体以及流体时:

    KmKP=11 (2)

    渗流速度v、渗流时间t的相似比为

    vmvP=KmimKPiP=11, (3)
    tmtP=lm/vmlP/vP=1N, (4)

    式中,角标为m的参量代表模型参量,角标为p的参量代表现场原型参量。

    研究表明[12],对直径不超过4 mm的多孔介质,在GrPr<1000时(Gr为格拉晓夫数,Pr普朗特数,这些是表征对流换热及辐射换热在传热过程中的相对重要性的无量纲参数),对流换热和辐射贡献微小,故本试验设计仅考虑热传导和相变传热。

    q=λdTdl, (5)

    式中,q为热流密度(W·m-2);λ为热传导介质的导热系数(W·m-1·K-1)。

    当原型和模型使用同一种土体以及流体时,材料导热系数相同。为保证试验温度比尺为1∶1,模型热流密度的相似比为

    qmqP=λmdTm/dlmλpdTp/dlp=dlpdlm=N1 (6)

    固体换热与相变对流换热、相变物质相变过程中所产生的潜热为

    ΔQ=Cphase1mphase1+Cphase2mphase2, (7)

    式中,Cphase1,mphase1分别为相变前物质的比热容和质量,Cphase2,mphase2为相变后物质的比热容和质量。当原型和模型使用同一种土体时,

    ΔQmΔQP=1N3 (8)

    模型试验装置一般是通过以几何相似常数为基础进行协调统一,但当子系统因为相似统一存在实现困难时,以主要控制变量为主,次要控制变量为辅进行协调设计。根据实际大小和模型箱大小几何相似常数1∶30为宜,但工程实践冻结管尺寸一般为89~105 mm,而模型试验在工业上能够制作模拟的冻结管直径最小为8 mm,无法满足几何相似常数的1∶30。但是像以往部分学者[17-18]采用的将冻结管作为冷源处理,只要在相似设计里满足傅里叶准则和科索维奇准则计算出相应的盐水温度即认为满足冻结管周围土体内部温度场的相似,存在着不合理性[19-20]。因为冻结管除了循环流体的温度,当流量不一样时,其吸热能力存在区别,也就是冻土的降温速度就不一样[21-22]。因此本文综合考虑冷源热对流相似和冻结管周围土体热传导相似,通过循环流体温度、流量及冻结管管径尺寸综合决定模型与原型的相似性,从而保证温度场变化规律的一致性,忽略冻结管系统本身的绝对几何相似。过程如下:

    根据旁通道冻结法技术规程(DG/TJ 08—902— 2016),冻结管吸热能力QT按下式计算:

    QT=Lπdq, (9)

    式中,L为冻结管长度(m);d为冻结管的直径(m);

    盐水泵流量W应按下式确定

    W=QΔtγc (10)

    式中 Q为1个冷冻站所需的冷量(kg/h);γ为盐水密度(kg/m3);c为盐水比热(kJ·kg-1·K);Δt为去回路盐水温差(K)。

    冷冻站制冷能力Q由下式确定:

    Q=mQT, (11)

    式中,m为冷量损失系数,取1.1~1.2。

    在不考虑热量损失的情况下,从能量守恒的角度来看,周围土体降温热传导放出热量的速率与冻结管中吸热能力相当,才能保证整个循环系统能量平衡,那么原型及模型的能量比值应该相等:

    QmQP=qmAmqPAm=qmπdmLmqpπdPLP=WmΔtmγmcmWpΔtpγpcp (12)

    冻结管的直径d,冻结管的根数n,盐水泵流量W,应满足:

    WmWp=QmQP/ΔtmγmcmΔtpγpcp=qmdmLmΔtpγpcpqpdPLPΔtmγmcm=N1dmLmΔtpγpcpdPLPΔtmγmcm (13)

    模型与原型的相似比采用1∶30。

    本模型试验原型为上海地区最为常见的埋置深度为12 m、冻结宽度为6 m的地铁隧道联络通道。经几何相似计算得冻结管埋深40 cm,冻结系统设计冻结管直径8 mm,数量为7根。冻结盐水泵流量经热流密度相似计算为30 L·min-1,渗流速度与原型相同。模型箱体两端使用带孔的有机玻璃隔板布设了进出水槽,进水水槽外墙顶部与出水槽底部分别布置3排不同高度的可拆卸进水阀门及进出水导管。箱体四周铺满保温隔热材料(见图2),以防止散热。现场试验与模型试验冻结管的流量设计都是依据冻结总管或工作泵。模型试验使用的冻结总管直径为40 mm,远大于冻结管直径,可以尽量保证冻结管处于满流状态。箱体两端有机玻璃隔板设置密集过水圆孔使水槽与试验土体进行水力联系。恒温水箱通过可拆卸控水阀门及进出水导管与进出水水槽相连,流量控制式清水泵置于恒温水箱内部,形成一定温度一定流量的渗流循环水体;一端出水导管连接有集水容量杯,从而精确复核试验土体的渗流循环流量。试验装置俯视图、传感器布置图如图2,3所示:温度传感器以相同深度处沿水平渗流方向以4或者5个串联布置一排测线,冻结管、冻结帷幕上下边缘及砂层处一共布设4排(见图3):土压力传感器在模型箱中心处沿深度在冻结管、冻结帷幕上缘及砂层处各布置1个,冻结帷幕下沿处由于渗流影响较大横向布置3个;位移传感器沿渗流方向在模型箱中线位置处分别布置5个,温度传感器采用热电偶串,尺寸为0.2 mm2精度0.01℃;土压传感器采用的应变式压力计,大小为0.6 cm2,精度0.2~0.4 kPa;位移传感器采用的差动式LVDT,探头仅0.1 cm2,精度0.01 mm。相比整个模型土体尺寸,监测探头尺寸效应不明显。试验土体初始温度为21℃(现场隧道附近监测温度),试验渗流水进入模型装置时的温度恒定为18℃(地下土体温度),试验环境使用空调设定为20℃(隧道通风温度)。

    图  2  试验装置俯视图
    Figure  2.  Top view of test devices
    图  3  传感器布置图
    Figure  3.  Layout of sensors

    试验的填土、传感器等布设完成后,对土体进行预固结。试验阶段给模型土体制冷后,实时监测土体内部各传感器的监测指标变化规律,当冻结帷幕扩展至规定厚度后土体监测温度全部降低至规定平均温度(-10℃),达到冻结设计要求;或当渗流速度特别大时,即使冻结足够长时间,仍然达不到-10℃,则当温度稳定后(1 h内温度降低不超过0.1℃),关闭制冷系统,开始融化。当温度全部回升到17℃左右且地表沉降不再增加时,试验结束。

    通过4排温度测线的布设,可以对试验土体重要温度变化区温度场进行严格的实时监测。渗流对冻结效果存在临界点:当渗流小于1.2 m·d-1时,冻结帷幕能够完成。渗流速度越大,冻结时间越长;当渗流大于1.2 m·d-1时,冻结将无法达到要求。由于0.5 m·d-1的流速适中,且冻结效果较有代表性,以0.5 m·d-1的工况为例,对温度分布特性进行分析;图4表明在没有渗流时模型土体的等温线相对水平,冻结帷幕形状近似为矩形;而图5显示,有渗流存在时,等温线呈一定程度的不均匀化偏移,上游等温线密集,下游等温线稀疏;且冻结帷幕的形状发生明显变化,冻结帷幕形状偏椭圆,呈入水端尖,出水端扁。渗流循环过程中,循环水体会吸收冻结土体的热量,对冻结帷幕的形成起到削弱作用。从监测的温度分布来看,上游的温度普遍高于下游,渗流对冻结壁整体都有削弱,但对上游削弱更明显。然而等温线并没有与渗流方向大角度相交,故渗流引入的渗流水体与周围软黏土的热对流作用仅仅是增加了其对土体冻结的影响程度,土体传热依然以热传导为主。

    图  4  冻结稳定时等温线分布图(0 m·d-1
    Figure  4.  Isotherm plot at freezing stability (0 m·d-1)
    图  5  冻结稳定时等温线分布图(0.5 m·d-1
    Figure  5.  Isotherm plot at freezing stability (0.5 m·d-1)

    冻结帷幕下缘是软黏土层和砂层的边界处,8号传感器位于中心位置受边界效应较小,因此以8号温度传感器为例,图6为不同渗流速度下8号传感器温度监测结果。土体在冻结阶段的温度变化曲线基本符合指数型特征;融沉阶段,除去未能降低至零下的大渗流速度的工况,温度阶梯型回升特性较明显。

    图  6  8号传感器温度监测图
    Figure  6.  Temperature monitoring of sensor No. 8

    除去渗流致使冻结帷幕无法冻结的工况,冻结帷幕下沿处8号传感器渗流速率与土体相变持续时间所对应关系如图7所示,可知渗流致使冻结土融化速率加快,这是因为区域内的渗流水体和周围土体存在热对流效应,且渗流速度越大,潜热释放的热量通过热对流被带走的越多,致使相变平衡时间随渗流速度的增大而线性减小。土体冻结降温阶段,相变平衡时间较短,主要原因是土体降温受冻结管影响较大,渗流对冻结所产生的影响不容易被传感器捕捉到。

    图  7  渗流速度与相变时间关系图
    Figure  7.  Relationship between seepage velocity and phase transition time

    图8统计了各渗流速度冻结帷幕上下缘冻结温度及冻结所需时间结果。渗流直接影响区域处的冻结帷幕下缘,当组合地层渗流速度大于1.2 m·d-1以上时,无法完成冻结设计需求,而在组合地层渗流速度1.2 m·d-1以内的,冻结帷幕可扩展到设计厚度,但渗流均延长了其达到冻结设计温度所需的冻结时间;而渗流间接影响区域处的冻结帷幕上缘,渗流的影响效果并不完全一致,尤其是最终稳定温度。渗流在延长冻结时间的同时,反而没有增强该区域的降温效果。

    图  8  上冻结帷幕下缘稳定温度及所需的冻结时间
    Figure  8.  Stable temperature and freezing time required at different positions

    综合以上分析,对于工程实际而言,在冻结帷幕能够达到设计的冻结效果时,并不是越大的渗流速度就越要加大注浆治理的投入。当组合地层渗流速度小于1.2 m·d-1时,若想缩短冻结时间,在冻结帷幕冻结管的内部调整优化即可。适当调整冻结管的排布方案,加密靠近渗流边界区域的冻结管,尽量使冻结液流动循环方向与渗流流动方向相同。当组合地层渗流大于1.2 m·d-1时,若想达到设计要求,需增加下缘冻结管数量,或在冻结帷幕下缘人为地降低地下水流速。上缘仍可在允许范围内增大间距。

    选取0.5 m·d-1各深度土压力监测结果(图9)进行对比分析。总的来说冻结管附近冻胀力集中效应明显,冻胀力下游明显大于上游,冻结帷幕上缘大于冻结帷幕下缘。在靠近渗流边界的区域,由于渗流循环水体和周围冻结土体的热交换,直接影响了冻结效果。渗流对冻胀力存在普遍削弱作用。

    图  9  0.5 m·d-1最大土压力(冻胀力)分布图
    Figure  9.  Distribution of maximum earth pressure (frost-heave force)

    除了冻胀力数值的区别,冻结帷幕上下缘冻胀力发展模式还存在较大不同,如图10所示。对于远离渗流边界的区域:冻结帷幕上缘以及冻结管附近区域(见土压力传感器2-4,2-8),冻结阶段具有一个突变增长的过程。冻结开始后,短时间内冻胀力猛增,增长至最大值后趋近于平稳。当进入融沉阶段时,土压力极速降低。对于靠近渗流边界的区域(见土压力传感器2-3,2-6),冻结和融解阶段土压力呈渐变型缓慢增加或减小。原因还是在靠近渗流边界的区域,渗流水体与冻结土体时刻处于直接接触状态,哪怕是降温到零度,该部分的土体也并不是处于完全冻结状态,而是依然存在热量交换的动态平衡,所以冻结帷幕下缘的冻胀力始终小于冻结帷幕上缘,且增长缓慢。

    图  10  0.5 m·d-1上冻结帷幕下缘土压力监测曲线
    Figure  10.  Monitoring curves of earth pressure (0.5 m·d-1)

    在实际工程施工中,冻胀力是监测预警的关键指标。对于冻结帷幕上缘部分,应该重点关注冻结前冻胀力的突变增长,应加强此阶段的预警及措施。而冻结帷幕下沿因为前期冻胀力增长较缓慢且数值较小,冻结后期应该是关注的重点。在融沉初始阶段,冻结帷幕上缘的快速消散应是工程注意的重点。

    对于冻融对地表位移的影响,以地表中心位置传感器1-6位移监测曲线为例,如图11。对于同一位置的地表变形来说,组合地层的渗流速度越大,冻胀和融沉的速率越小;地表最终冻结隆起量和融化沉降量也越小。对于融化后地表沉降平面形状,如图12,可以看出沿着渗流方向,地表沉降越来越小。且渗流速度越大,这种不均匀沉降差异也越明显。渗流速度越大,最终沉降量反而有所降低。实际工程中,对于地表环境而言,不均匀沉降也是需要关注的灾害。

    图  11  不同渗流速度下位移传感器1-6位移监测图
    Figure  11.  Monitoring of displacement sensors 1-6 under different seepage velocities
    图  12  地表最终变形曲线
    Figure  12.  Curves of surface deformation

    渗流与融沉速率的影响关系如图13可知,在小渗流的情况下,渗流速度增大将降低融沉速率,也就大大地增长了融沉处理的时间,这对实际工程也是不利的。产生这种现象的原因是渗流速度越大的土体冻结后渗透系数越小,这就在一定程度上增长了排水固结的时间。

    图  13  渗流与融沉速率关系图
    Figure  13.  Relationship between seepage and thawing rate

    渗流可以对组合地层冻结时间产生重大影响。为降低模型试验不可避免的偶然误差与系统误差,本文将可以达到设计要求的渗流工况下的模型试验温度监测曲线进行拟合与处理[23],模型采用热传导的一维解模型(拟合结果及方差见图14表3,表明拟合效果较好)。冻结帷幕下缘设计冻结厚度处的土体,有渗流与无渗流工况下温度降低到同一温度时,温度降低速率的比值关系为

    dTvdt/dT0dt=1/1.25v, (14)
    图  14  8号温度传感器温度时程曲线拟合结果
    Figure  14.  Fitting results of time-history curve of temperature sensor No. 8
    表  3  温度拟合方差表
    Table  3.  Variances of temperature fitting
    渗流速度/(m·d-1)方差
    00.907
    0.300.935
    0.500.878
    0.750.901
    1.200.954
    3.000.789
    下载: 导出CSV 
    | 显示表格

    式中 dTv/dt为含有渗流工况冻结帷幕下缘降低到某一温度时的温度降低速率(K·h-1);dT0/dt为无渗流工况冻结帷幕下缘降低到某一温度时的温度降低速率(K·h-1)。

    对温度监测曲线进行拟合,得有、无渗流的工况下,完成冻结施工所需时间的关系为

    tv=t0e0.803v, (15)

    式中tv为含有渗流工况冻结帷幕下缘扩展到设计厚度所需的时间(h);t0为无渗流工况冻结帷幕下缘扩展到设计厚度所需的时间(h),为22.6 h。

    将模型试验数据根据相似比进行放大以反映实际工程,可对组合地层中具有渗流工况下现场冻结施工所需的时间进行预测,如表4

    表  4  渗流工况下现场冻结施工所需的时间预测表
    Table  4.  Time required for on-site freezing construction
    渗流速度/(m·d-1)00.10.20.30.40.50.60.70.80.81.01.1
    时间/d22.624.526.528.731.233.636.639.743.046.750.654.8
    下载: 导出CSV 
    | 显示表格

    针对不同渗流的现场冻结施工工况,本文给出建议列于表5

    表  5  现场冻结施工建议措施表
    Table  5.  Suggestions for freezing construction
    渗流速度/(m·d-1)建议
    0~0.5正常施工并做好延长工期的计划与安排
    0.5~1.2可打防渗帷幕或适当采用抽水的方式降低流速
    1.2~须采取特殊措施降低冻结区域地下水流速,如打竖排的冻结管,制作防渗帷幕,人工降水等措施
    下载: 导出CSV 
    | 显示表格

    在靠近渗流边界的冻结帷幕下沿以及砂层区域,渗流致使冻结阶段所产生的冻胀力随渗流速度的增大而减小,故随着渗流速度的增大隧道所需要抵抗冻胀所造成的弯矩将变小,隧道变形与整体位移将会减小。渗流对隧道变形控制来说是一个有利的因素[24-25]。隧道变形示意图,如图15所示。

    图  15  隧道变形示意图
    Figure  15.  Tunnel deformation

    对于冻结区域隧道每延米受到的冻胀力合力:

    E=topbottomPK0dl (16)

    式中 E为冻结区域每延米隧道预计会受到的冻胀力的合力(MPa﹒m);P为竖直向的冻胀力,为竖直向土压力计监测到的土压力值减初始土压力监测值(Pa);l为隧道上某点相对于隧道底部的垂直高度(m);K0为冻结土体主动土压力系数。

    根据模型试验数据与相似比可以对现场工况进行预测。将数据拟合并计算整理得渗流工况下冻结区域每延米隧道所受到的冻胀合力预测值如表6所示。

    表  6  渗流速度对隧道所受冻胀合力影响预测
    Table  6.  Effects of seepage velocity on frost-heave force on tunnel
    渗流速度/(m·d-1)00.250.500.751.0
    冻胀合力/(MPa·m)24.722.220.218.817.0
    下载: 导出CSV 
    | 显示表格

    根据模型试验数据与相似比可对现场地表沉降量进行预测[26-27]。针对工程实际的地面沉降预测与治理建议见表7,8

    表  7  渗流速度与地表变形关系表
    Table  7.  Relationship between seepage velocity and surface deformation
    渗流速度/(m·d-1)最大隆起/mm最大沉降/mm平均沉降梯度(0.001)
    06.368.71.53
    0.255.453.21.72
    0.504.547.61.93
    0.753.843.72.03
    1.03.541.52.10
    下载: 导出CSV 
    | 显示表格
    表  8  渗流工况下地表变形防治注意事项表
    Table  8.  Precautions for prevention of surface deformation under seepage conditions
    渗流速度/(m·d-1)建议
    0应注意防治较大的地表沉降量
    其他注意防治不均匀沉降,以及长期沉降
    下载: 导出CSV 
    | 显示表格

    (1)渗流速度在3 m·d-1以内时,砂层渗流边界增加了渗流循环水体对周围软黏土的热对流作用,改变了等温线分布形状,但软黏土层的热传导依然为土体主要传热方式。

    (2)渗流会使组合地层冻结时间增长,会减小渗流边界处的冻胀力,会减小地表因冻融所产生的隆起与沉降。

    (3)组合地层渗流对软黏土冻结帷幕上下缘的冻结效果差别较大,尤其存在渗流速度临界点。组合地层渗流处于临界速度以内时,若想缩短冻结时间,在冻结帷幕冻结管的内部调整优化即可,无需盲目增大冻结管投入。而超过临界速度时,需增加下缘冻结管数量,或在冻结帷幕下缘人为地降低地下水流速。上缘仍可在允许范围内增大间距。

    (4)渗流削弱了直接影响区域冻结软黏土的潜热影响效果,即缩短了相变平衡时间。相变平衡时间随渗流速度增大呈线性减小。

    (5)组合渗流对冻结帷幕上下缘内部的冻胀力分布及其发展模式的影响存在差异。对于冻结帷幕上缘,冻结前期就需要重点关注,及时采取预警及措施;而冻结帷幕下沿,更应关注冻结后期的冻胀力主要增长期。

    (6)渗流速度越大,沿渗流方向地表的不均匀沉降差异也越明显。对于实际工程的地表环境而言,不均匀沉降也是需要关注的灾害。

    (7)渗流会致使隧道所需承受的冻胀合力减小,这有利于冻结工程隧道的变形控制。

  • 图  1   试验设备

    Figure  1.   Test apparatus

    图  2   边坡模型概化设计示意图

    Figure  2.   Generalized design diagrams of slope models

    图  3   试验边坡模型

    Figure  3.   Slope models for tests

    图  4   试验地震波加速度时程

    Figure  4.   Time histories of acceleration of seismic wave for tests

    图  5   加速度绝对传递函数虚部

    Figure  5.   Imaginary part of absolute transfer function of acceleration

    图  6   自振频率和阻尼比随白噪声序列号的变化

    Figure  6.   Variation of natural frequency and damping ratio with white noise sequence numbers

    图  7   各阶段测点PGA放大系数随h/Hb/B的变化

    Figure  7.   Variation of PGA amplification factor at monitoring points with h/H and b/B at different stages

    图  8   边坡岩体非线性累积损伤演化曲线

    Figure  8.   Evolution curves of nonlinear cumulative damage of slope rock mass

    图  9   含起伏体边坡累积损伤-失稳破坏演化过程

    Figure  9.   Evolution of cumulative damage and instability failure of slope with undulant body

    图  10   无起伏体边坡累积损伤-失稳破坏演化过程

    Figure  10.   Evolution of cumulative damage and instability failure of slope without undulant-body

    图  11   边坡失稳破坏滑面形态

    Figure  11.   Sliding surface shape of slope after instability failure

    图  12   数值计算模型

    Figure  12.   Numerical model

    图  13   边坡累积永久位移和稳定性系数的变化

    Figure  13.   Variation of cumulative permanent displacement and stability factor of slope

    图  14   边坡累积损伤-失稳破坏演化过程

    Figure  14.   Evolution of cumulative damage and instability failure of slope

    图  15   多影响因素下边坡累积永久位移和稳定性系数的变化

    Figure  15.   Variation of cumulative permanent displacement and stability factor under different influencing factors

    表  1   关键物理量及其相似系数

    Table  1   Key physical quantities and their similitude coefficients

    关键物理量量纲相似关系相似系数
    密度ML3 Cρ 1
    弹性模量ML1T2 CE 32.6
    泊松比无量纲Cμ 1
    黏聚力ML1T2 Cc=CECε 22.8
    内摩擦角无量纲Cφ 1
    应力ML1T2 Cσ=CECε 22.8
    应变无量纲Cε=CρCgClC1E 0.7
    几何尺寸LCl=C0.5EC0.5ρCt 22.8
    位移LCu=ClCε 16
    时间TCt 4
    频率T1 Cf=C1t 0.25
    速度LT1 Cv=CuC1t 4
    运动加速度LT2 Ca=CuC2t 1
    重力加速度LT2 Cg=CuC2t 1
    阻尼比无量纲Cζ 1
    下载: 导出CSV

    表  2   相似材料力学特性参数

    Table  2   Mechanical property parameters of equivalent materials

    项目密度/(g·cm-3)抗压强度/MPa弹性模量/MPa泊松比黏聚力/kPa内摩擦角/(°)
    三峡库区范围2.0~2.530~700.15~4.00.18~0.32~1030~45
    典型原岩2.4934.142000.3200041.8
    相似系数122.832.6132.61
    目标值2.491.50128.80.361.341.8
    实际值2.501.401800.2512042
    下载: 导出CSV

    表  3   试验加载工况

    Table  3   Loading conditions for tests

    加载工况动荷载
    加载类型振幅/g频率/Hz持时/s加载次数
    1白噪声10.08301
    2正弦波0.043071
    3天然波0.043271
    4正弦波0.04571
    5正弦波0.041071
    6正弦波0.041571
    7正弦波0.043071
    8天然波0.043271
    9天然波0.063271
    10天然波0.093271
    11天然波0.153271
    12白噪声20.08301
    13~62(微震)天然波0.0432750
    63白噪声30.08301
    重复加载工况13~63,直至微震作用500次
    523(震中高烈度)天然波0.35326.51
    524(小震)白噪声130.08301
    525~574正弦波0.1530750
    重复加载工况524~574,直至小震作用500次
    1035白噪声230.08301
    持续施加振幅为0.2g~0.4g、频率为15 Hz和持时为7 s的正弦波直至边坡失稳破坏
    下载: 导出CSV

    表  4   数值计算工况

    Table  4   Numerical conditions

    工况层面倾角β/(°)层面间距Ld/m坡角α/(°)坡高H/m动荷载振幅/g动荷载频率/Hz
    #102.36022.90.0430
         0.045
    #2302.36022.90.1520
         0.230
    #3602.36022.90.0430
    #4301.56022.90.0430
    #5300.56022.90.0430
    #6302.33022.90.0430
    #7302.34522.90.0430
    #8302.36050.00.0430
    #9302.360150.00.0430
    下载: 导出CSV

    表  5   岩体与边界参数

    Table  5   Parameters of rock mass and boundary

    项目密度/(g·cm-3)弹性模量/MPa泊松比内摩擦角/(°)黏聚力/kPa抗拉强度/MPa
    岩体2.4942000.341.840002
    边界2.5042000.3
    下载: 导出CSV

    表  6   层面与节理参数

    Table  6   Parameters of bedding surface and joint

    项目内摩擦角/(°)黏聚力/kPa抗拉强度/MPa法向刚度/(GPa·m-1)切向刚度/ (GPa·m-1)
    层面32300.08010.26.4
    节理2060.0038.25.1
    下载: 导出CSV
  • [1] 黄润秋. 岩石高边坡稳定性工程地质分析[M]. 北京: 科学出版社, 2012.

    HUANG Run-qiu. Engineering Geology for High Rock Slopes[M]. Beijing: Science Press, 2012. (in Chinese)

    [2] 王儒述. 三峡水库诱发地震的监测与预报[J]. 三峡大学学报(自然科学版), 2017, 39(1): 30-37. https://www.cnki.com.cn/Article/CJFDTOTAL-WHYC201701025.htm

    WANG Ru-shu. Induced earthquake monitoring and forecast of Three Gorges reservoir[J]. Journal of China Three Gorges University (Natural Sciences), 2017, 39(1): 30-37. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-WHYC201701025.htm

    [3] 李兴旺. 频发微震作用下顺层岩质边坡的稳定性研究[D]. 重庆: 重庆大学, 2015.

    LI Xing-wang. The Analysis of the Stability of Bedding Slope Under the Effect of Frequent Seismic[D]. Chongqing: Chongqing University, 2015. (in Chinese)

    [4] 王思敬. 岩石边坡动态稳定性的初步探讨[J]. 地质科学, 1977(4): 372-376. https://www.cnki.com.cn/Article/CJFDTOTAL-DZKX197704005.htm

    WANG Si-jing. Preliminary notes on the dynamic stability of rock slopes[J]. Chinese Journal of Geology (Scientia Geologica Sinica), 1977(4): 372-376. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DZKX197704005.htm

    [5]

    GRIFFITHS D W, BOLLINGER G A. The effect of appalachian mountain topography on seismic waves[J]. Bulletin of the Seismological Society of America, 1979, 69(4): 1081-1105. doi: 10.1785/BSSA0690041081

    [6]

    HONG Y S, CHEN R H, WU C S, et al. Shaking table tests and stability analysis of steep nailed slopes[J]. Canadian Geotechnical Journal, 2005, 42(5): 1264-1279. doi: 10.1139/t05-055

    [7] 黄润秋, 李果, 巨能攀. 层状岩体斜坡强震动力响应的振动台试验[J]. 岩石力学与工程学报, 2013, 32(5): 865-875. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201305003.htm

    HUANG Run-qiu, LI Guo, JU Neng-pan. Shaking table test on strong earthquake response of stratified rock slopes[J]. Chinese Journal of Rock Mechanics and Engineering, 2013, 32(5): 865-875. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201305003.htm

    [8] 董金玉, 杨国香, 伍法权, 等. 地震作用下顺层岩质边坡动力响应和破坏模式大型振动台试验研究[J]. 岩土力学, 2011, 32(10): 2977-2982, 2988. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201110015.htm

    DONG Jin-yu, YANG Guo-xiang, WU Fa-quan, et al. The large-scale shaking table test study of dynamic response and failure mode of bedding rock slope under earthquake[J]. Rock and Soil Mechanics, 2011, 32(10): 2977-2982, 2988. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201110015.htm

    [9] 罗璟, 裴向军, 黄润秋, 等. 强震作用下滑坡岩体震裂损伤程度影响因素研究[J]. 岩土工程学报, 2015, 37(6): 1105-1114. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201506020.htm

    LUO Jing, PEI Xiang-jun, HUANG Run-qiu, et al. Influencing factors for damage degree of shattered landslide rock mass under high seismic action[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(6): 1105-1114. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201506020.htm

    [10] 许强, 刘汉香, 邹威, 等. 斜坡加速度动力响应特性的大型振动台试验研究[J]. 岩石力学与工程学报, 2010, 29(12): 2420-2428. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201012009.htm

    XU Qiang, LIU Han-xiang, ZOU Wei, et al. Large-scale shaking table test study of acceleration dynamic responses characteristics of slopes[J]. Chinese Journal of Rock Mechanics and Engineering, 2010, 29(12): 2420-2428. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201012009.htm

    [11]

    LIU X R, HE C M, LIU S L, et al. Dynamic response and failure mode of slopes with horizontal soft and hard interbeddings under frequent microseisms[J]. Arabian Journal for Science and Engineering, 2018, 43(10): 5397-5411.

    [12] 杨忠平, 刘树林, 刘永权, 等. 反复微震作用下顺层及反倾岩质边坡的动力稳定性分析[J]. 岩土工程学报, 2018, 40(7): 1277-1286. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201807018.htm

    YANG Zhong-ping, LIU Shu-lin, LIU Yong-quan, et al. Dynamic stability analysis of bedding and toppling rock slopes under repeated micro-seismic action[J]. Chinese Journal of Geotechnical Engineering, 2018, 40(7): 1277-1286. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201807018.htm

    [13] 刘永权. 频发微震下库区顺层岩质边坡累积损伤演化机理及稳定性研究[D]. 重庆: 重庆大学, 2017.

    LIU Yong-quan. Study on Cumulative Damage Evolution Mechanism and Stability of Bedding Rock Slope in Reservoir Area Under Frequent Microseismic[D]. Chongqing: Chongqing University, 2017. (in Chinese)

    [14]

    LIU X R, LIU Y Q, HE C M, et al. Dynamic stability analysis of the bedding rock slope considering the vibration deterioration effect of the structural plane[J]. Bulletin of Engineering Geology and the Environment, 2018, 77(1): 87-103.

    [15] 江洎洧. 频发微震作用下三峡库首段典型滑坡变形机制及动力响应研究[D]. 武汉: 中国地质大学, 2012.

    JIANG Ji-wei. Research on the Deformation Mechanism and Dynamic Response of Typical Landslides in Three Gorges Reservoir in Case of Frequent Microseisms[D]. Wuhan: China University of Geosciences, 2012. (in Chinese)

    [16] 顾大钊. 相似材料和相似模型[M]. 徐州: 中国矿业大学出版社, 1995.

    GU Da-zhao. Equivalent Materials and Similitude Models[M]. Xuzhou: China University of Mining & Technology Press, 1995. (in Chinese)

    [17] 杨静黎. 三峡库区滑坡岩土力学参数分区研究[D]. 重庆: 重庆交通大学, 2014.

    YANG Jing-li. Study on Rock and Soil Mechanics Parameters Partition of Landslides in the Three Gorges Reservoir[D]. Chongqing: Chongqing Jiaotong University, 2014. (in Chinese)

    [18] 陈剑, 李晓, 杨志法. 三峡库区滑坡的时空分布特征与成因探讨[J]. 工程地质学报, 2005, 13(3): 305-309. https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ200503003.htm

    CHEN Jian, LI Xiao, YANG Zhi-fa. On the distribution and mechanism of landslides in the three gorges reservoir area[J]. Journal of Engineering Geology, 2005, 13(3): 305-309. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCDZ200503003.htm

    [19] 曹树谦, 张文德, 萧龙翔. 振动结构模态分析——理论、实验与应用[M]. 2版.天津: 天津大学出版社, 2014.

    CAO Shu-qian, ZHANG Wen-de, XIAO Long-xiang. Modal Analysis of Vibration Structures-Theory, Experiments and Applications[M]. 2nd ed. Tianjin: Tianjin University Press, 2014. (in Chinese)

    [20] 朱传云, 喻胜春. 爆破引起岩体损伤的判别方法研究[J]. 工程爆破, 2001, 7(1): 12-16. https://www.cnki.com.cn/Article/CJFDTOTAL-GCBP200101002.htm

    ZHU Chuan-yun, YU Sheng-chun. Study on the criterion of rock mass damage caused by blasting[J]. Engineering Blasting, 2001, 7(1): 12-16. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GCBP200101002.htm

    [21]

    NAKATA N, SNIEDER R, KURODA S, et al. Monitoring a building using deconvolution interferometry: I Earthquake- data analysis[J]. Bulletin of the Seismological Society of America, 2013, 103(3): 1662-1678.

  • 期刊类型引用(20)

    1. 杨哲,蔡海兵,王彬,李孟凯,庞昌强. 高流速富水砂层地铁隧道液氮冻结温度场演化规律研究. 铁道科学与工程学报. 2025(01): 307-319 . 百度学术
    2. 谭智勇,王超林,龙安发. 外部水源作用下岩石液氮冻结试验研究. 岩土工程学报. 2024(02): 415-425 . 本站查看
    3. 刘爽,李晓康,李旭,聂雯,林旸,盛志刚. 渗流作用下粗粒土冻结壁交圈规律及预测模型探索. 冰川冻土. 2024(01): 247-259 . 百度学术
    4. 刘欣,冯利华,任全军,沈宇鹏,韩昀希,刘越,韩风雷. 渗流地层地铁联络通道冻结壁形成过程及其影响因素分析. 中南大学学报(自然科学版). 2024(04): 1463-1476 . 百度学术
    5. 孙立强,时鹏,郎瑞卿,商安策. 渗流作用下人工冻结特性室内模型试验研究. 岩石力学与工程学报. 2024(S1): 3530-3542 . 百度学术
    6. 孙立强,商安策,郎瑞卿,苗雨. 渗流地层人工冻结壁交圈时间计算方法. 岩石力学与工程学报. 2023(S1): 3663-3673 . 百度学术
    7. 樊文虎,杨平,王升福. 粉质黏土冻融前后细观结构试验研究. 森林工程. 2023(03): 182-190 . 百度学术
    8. 宋帅,姜艳,尹文纲. 冻结法施工中渗流场对温度场的影响分析. 徐州工程学院学报(自然科学版). 2023(04): 68-76 . 百度学术
    9. 荣传新,王彬,程桦,董艳宾,杨凡. 大流速渗透地层人工冻结壁形成机制室内模型试验研究. 岩石力学与工程学报. 2022(03): 596-613 . 百度学术
    10. 陈盈盈,周桂云,余长青. 渗流对人工冻结温度场影响的研究. 四川水泥. 2022(03): 61-63 . 百度学术
    11. 刘利. 冻结法施工在某隧道扩挖中的应用. 市政技术. 2022(04): 144-148+192 . 百度学术
    12. 李忠超,白天麒,梁荣柱,肖铭钊,蔡兵华,叶超,吴文兵. 富水粉细砂层水平冻结效果试验及数值模拟. 工业建筑. 2022(03): 1-9 . 百度学术
    13. 黄建华,严耿明,覃少杰. 液氮冻结加固冻结管内换热机制及对流换热系数研究. 岩土力学. 2022(09): 2624-2633 . 百度学术
    14. 徐陈明,胡俊,熊辉,林小淇,周禹暄,冯继超,王志鑫. 一种基坑坑底加固结构温度场变化规律研究. 海南大学学报(自然科学版). 2022(04): 434-440 . 百度学术
    15. 丁飞,张鲁鲁,曹新刚,李远荣,汪亦显. 微冻结法盾构隧道衬砌管片冻结性能参数试验研究. 工程与建设. 2022(05): 1379-1381+1416 . 百度学术
    16. 赵象卓,王春林,王宇,李彦民,姬强,闫永乐. 富水岩层井筒超长冻结自然解冻趋势与注浆技术. 煤炭工程. 2021(08): 43-48 . 百度学术
    17. 刘娟. 江底强渗透地层冻结法可行性分析. 江西建材. 2021(08): 193-195 . 百度学术
    18. 周洁,李泽垚,田万君. 渗流作用对人工冻结土体特性的影响. 铁道工程学报. 2021(07): 12-17 . 百度学术
    19. 汪恩良,任志凤,韩红卫,田雨,胡胜博,刘兴超. 超低温冻结黏土单轴抗压力学性质试验研究. 岩土工程学报. 2021(10): 1851-1860 . 本站查看
    20. 周扬,武子寒,许程,卢萌盟,周国庆. 高温下饱和冻土一维融化热固结模型及解答. 岩土工程学报. 2021(12): 2190-2199 . 本站查看

    其他类型引用(11)

图(15)  /  表(6)
计量
  • 文章访问数:  409
  • HTML全文浏览量:  26
  • PDF下载量:  223
  • 被引次数: 31
出版历程
  • 收稿日期:  2019-05-19
  • 网络出版日期:  2022-12-07
  • 刊出日期:  2020-03-31

目录

/

返回文章
返回