Impermeability characteristics of junctional zone between compacted broadly graded clayey soil and hard surface
-
摘要: 研发了土料与刚性面结合带在相互正交、独立可控的压应力、剪切变形和水力比降作用下抗渗性能试验的大型和中型设备,最大剪切变形量值大于1 m。对压实天然宽级配黏土开展了结合带抗渗特性试验研究,发现了试样渗透系数随剪切变形突然增大后逐步减小并稳定的变化规律,剪切变形启停历史会弱化结合带渗透系数随剪切变形启动而突然增大的现象,初始密度高、应力高、试样尺寸大的土样渗透系数更低。考虑结合带导水系数与初始值的比值和试样渗透系数与初始值的比值之间的关系,分析渗透系数剪切后达峰值时结合带厚度在0.075~1.0 mm之间,结合带厚度随正应力减小、随试样尺寸增大而增厚。提出改进的黏土颗粒不规则形状集簇模型,通过组构变化较好解释了剪切过程中渗透性持续演化的机理,也支撑了结合带厚度的分析结论。Abstract: Large-scale and medium-scale test equipments for impermeability characteristics of junctional zone between soil and hard surface are developed. The mutually orthogonal compressive stress, shear deformation and hydraulic gradient of the junctional zone between the soil and the rigid surface are controllable independently. The maximum shear deformation is greater than 1 m. The experimental studies on a compacted natural broadly graded clayey soil are carried out. The variation law of the permeability coefficient of the samples is discovered, that the permeability coefficient decreases gradually and then becomes stable after the sudden increase along with the start of the shear deformation. The history of shear deformation may weaken the phenomenon of the sudden increase of the permeability coefficient after the start of shear deformation. The permeability coefficient of soil samples with higher initial density, higher stress and larger external size is lower. Based on the corresponding relationship between the ratio of water conductivity coefficient of the junctional zone to its initial value and the ratio of permeability coefficient of the sample to its initial value, the possible thickness of the junctional zone considered to be between 0.075 mm and 1.0 mm is analyzed when the peak value of permeability coefficient is achieved. A revised model for clay particle clusters with irregular shape is proposed. The continuous change mechanism of permeability during shear deformation can be explained by fabric adjustment. The conjecture of the thickness of the junctional zone is also supported.
-
0. 引言
岩体工程多赋存于水环境中,并常处于饱和状态。环境中的水进入岩体孔隙后易发生水岩相互作用,将改变岩石的矿物组成和微细观结构,引起强度、变形等物理力学性能的劣化。为保证水环境中岩体工程的稳定性,研究饱和岩石的水软化特性和机制十分必要。
目前,许多学者对多种饱和岩石的水软化特性进行了研究。普遍认为,岩石饱和后强度降低,弹性模量减小,泊松比增加,摩擦系数减小。如周翠英等[1]发现泥岩饱和后,其强度下降了约90%。Mann等[2]对3种饱和砂岩进行压缩试验,发现弹性模量减小了6%~19%。熊德国等[3]得到砂岩、砂质泥岩和泥岩的抗压强度水软化系数(饱和岩石抗压强度与干燥岩石抗压强度之比)范围为0.58~0.94。Zhu等[4]得到大理岩饱和后,单轴抗压强度减小了29.18%,泊松比增加13.90%。黄弈茗[5]通过剪切试验得到饱和灰岩的摩擦系数降低了25%。Chen等[6]、Erguler等[7]和Vasarhelyi[8]发现凝灰岩饱和后强度降低了2%~88%。工程上常将单轴抗压强度小于20 MPa的岩石称为软岩,20~40 MPa的称为中硬岩,大于40 MPa的称为硬岩。总结来看,目前软岩的水软化特性研究成果丰富,硬岩的水软化特性研究存在不足。
为合理解释饱和岩石力学性能的劣化现象,一些学者提出了诸如黏土矿物软化、矿物颗粒摩擦力下降、孔隙水压力等水软化机制[9-13]。然而,受限于试验手段和分析方法,饱和岩石力学性能劣化的物理力学机制尚未达成共识,试验成果多,理论分析成果少。Zhu等[4]综合已有试验成果发现,硬岩的水软化机制可能包括:机制一,黏土矿物软化;机制二,孔隙水压力作用;机制三,摩擦弱化效应;机制四,应力腐蚀。其中,机制一、四分别存在于含黏土矿物和石英矿物的岩石中,且机制一已被公认。机制二、三(统称水岩作用机制,下同)普遍存在于各类饱和岩石中,受到了广泛关注。值得提出的是,上述水软化机制常常同时存在,并相互影响。因此,排除其他针对某一个或几个特定水软化机制,定性和定量描述其存在条件和作用程度成为亟待解决的问题。
声发射(acoustic emission,AE)是材料加载过程中局部应变能释放所产生的瞬态弹性波,因其可以反映材料内部裂隙发育和破坏特性而被广泛应用于岩石力学中[14-15]。声发射信号的分析方法主要包括参数分析和波形分析。目前一些学者利用声发射特征参数,如声发射事件率、声发射振铃计数、能量等,研究含水岩石的破坏特性[16-18]。虽然声发射参数分析能反映饱和岩石破坏的一些释能特征,但将复杂难懂的波形信号抽象成一系列的参数来加以描述,难免会遗失波形所携带的一些重要的信息。近年来,一些学者开始研究声发射波形信号的频谱特征[14, 19-20]。关于声发射频谱在含水岩石破坏特性研究中的运用,少量学者收获了十分有限的成果。例如,Chugh等[21]、Read等[22]分析了干燥和含水岩石的声发射波形信号的主频差异。张艳博等[23]发现饱和砂岩的声发射信号主频主要分布在35~60 kHz,且主频范围比干燥砂岩集中。然而,基于声发射波形,特别是波形信号的主频统计特征,分析饱和岩石的微观破坏特性进而研究水岩作用机制的成果鲜有报道。
为分析饱和岩石水岩作用机制的存在条件和作用程度,本文选取不含黏土和硅质矿物的典型硬岩,以排除机制一、四的影响,开展硬岩在干燥、饱和状态下的单轴压缩声发射试验,分析压缩强度、弹性模量、破裂形态等力学软化特性,探讨饱和硬岩破坏的声发射波形信号主频统计特征。结合力学软化特性和声发射主频统计特征,研究硬岩水岩作用机制的存在条件和作用程度。研究成果可丰富饱和硬岩力学软化特性的研究内容,并为量化描述水岩作用机制的研究提供新的视角和途径。
1. 试验方法与过程
1.1 试样制备
试验用岩石包括3种:①宝兴大理岩,取样自四川省雅安市宝兴县,②锦屏大理岩,取样自四川省凉山彝族自治州盐源县锦屏一级水电站库岸边坡,③灰岩,取样自四川省凉山彝族自治州雷波县马湖边坡。依据国际岩石力学学会(ISRM)建议方法将岩芯加工成50 mm×H100 mm的标准圆柱体试样。如图 1所示,采用X射线单晶衍射仪测定3种岩石的矿物成分和各矿物占比,发现宝兴大理岩和灰岩均由单一矿物方解石组成,锦屏大理岩由方解石(12%)、白云石(87%)、氯化钾(0.5%)、氧化铁(0.5%)组成。说明3种岩石均不含黏土矿物和石英矿物,黏土矿物软化和应力腐蚀作用可以忽略。电镜扫描(SEM)图像显示3种岩石的矿物胶结形态、孔隙数量存在较大差异。利用核磁共振岩芯分析系统(NMR)测定试样的孔隙度。试样的物理特性参数如表 1所示。
表 1 试样基本物理性质参数Table 1. Basic physical parameters of samples岩石类型 天然密度/(kg·m-3) 孔隙度/% 纵波波速/(m·s-1) 矿物成分 宝兴大理岩 2.69 (±0.04) 6.16 (±0.25) 4262 (±30.73) 方解石 锦屏大理岩 2.98 (±0.17) 2.68 (±0.13) 4689 (±29.32) 方解石、白云石、氯化钾、氧化铁 灰岩 2.63 (±0.11) 4.80 (±0.09) 3826 (±18.67) 方解石 注:括号内数值为试样样本参数的标准差。 为获得干燥、饱和试样,将加工好的试样放入105℃的恒温烘箱内烘48 h,使其达到完全干燥状态。冷却至室温后取一半试样置于干燥箱内保存。采用强制饱和法,将剩余试样置于100 kPa的负压真空抽气机中强制饱和48 h,并静置4 h。随机选取3种岩石的干燥、饱和试样各2个,采用干馏法测定水饱和度,发现其值分别接近0和1,认为获得的干燥、饱和试样满足试验要求。
1.2 试验设备与方案
试验采用MTS815 Flex Test GT岩石力学试验系统和PCI-2声发射采集系统进行。每个试样布置8个Micro30高精度传感器,可以有效采集试验过程的声发射信号[4, 13]。声发射前置放大器门槛值(增益)为40 dB,采样门槛值为30 dB。试样与传感器间涂抹耦合剂,以保证信号数据的有效采集。试验设备与传感器布置如图 2所示。试验采用LVDT进行轴向位移控制,加载速率为0.05 mm/min。每种岩石包括干燥、饱和试验组各1组,每组4个试样,共计24个。试验过程中,保证力学加载系统与声发射采集系统的时间同步。
2. 试验结果
2.1 力学软化特性
3种岩石典型试样的应力应变曲线如图 3所示。可以看出,试样在单轴压缩荷载作用下的应力应变曲线具有明显的阶段性,即加载初期的非线性段(压密段)、线弹性段、峰值应力前的塑性变形段和峰后段。饱和试样的峰值应变(峰值应力对应的应变)比干燥试样的大,表现出应变软化特性。相比较而言,锦屏大理岩力学性能最强,灰岩次之,宝兴大理岩最弱。
干燥、饱和状态下3种岩石的单轴抗压强度、弹性模量E及软化系数f如图 4,5所示。其中,软化系数f为饱和试样组的力学参数均值与干燥试样组的力学参数均值的比值,和分别表示单轴抗压强度软化系数和弹性模量软化系数。
可以看出,3种岩石饱和后,力学性能均出现明显劣化,表现为单轴抗压强度和弹性模量减小。宝兴大理岩、锦屏大理岩、灰岩的单轴抗压强度软化系数分别为0.74,0.94,0.76,弹性模量软化系数分别为0.72,0.92,0.79。相比较而言,宝兴大理岩水软化性能最强,灰岩次之,锦屏大理岩最弱。
干燥、饱和状态下典型试样的宏观破坏形态与声发射事件定位见图 6,图中声发射事件的颜色对应不同的加载时间。可以看出,声发射事件多发生于试样的宏观破坏面或裂纹处,宏观破裂形态与声发射事件定位能较好吻合。相比较而言,锦屏大理岩破坏时释放的声发射信号比宝兴大理岩和灰岩的多。在压缩荷载作用下,干燥岩石多发生剪切破坏,试样可见明显的单斜面或共轭斜面剪切破坏面。岩石饱和后,平行于试样轴向的张拉裂纹和张拉破裂面增多。
2.2 声发射波形主频特性
(1)数据处理
试验采集了加载全过程岩石释放的声发射波形信号。本文采用声发射信号频谱统计分析方法,即将所有通道采集的海量信号数据进行统计分析。首先,通过MATLAB编程,采用快速傅立叶变换(FFT)获得声发射波形信号的主频参数(图 7)。其次,对获得的波形信号主频进行区段划分。室内岩石的声发射频率位于10~500 kHz范围内[24],且一般低于450 kHz。为便于统计分析,参照李林芮[25]、邓建辉等[26]的研究,采用的主频区段划分方法为:主频0~450 kHz,按每10 kHz划分一个主频段,共45个主频段;主频451~500 kHz,划分为一个主频段,合计46个主频段。
(2)声发射主频特征
干燥、饱和状态下3种岩石典型试样的主频占比统计见图 8。可以发现,在压缩荷载作用下,干燥、饱和岩石的声发射波形主频均集中在2个主频带范围内,称之为高主频带和低主频带。作为邓建辉等[26]关于岩石声发射双主频特性的后续研究,说明双主频现象是岩石破坏的一种固有特性,与岩石类型和含水状态无关。岩石饱和后,高、低主频带范围发生变化。以宝兴大理岩为例,干燥宝兴大理岩的高、低主频带分别集中在230~320 kHz(第24~32主频段)和11~80 kHz(第2~8主频段),饱和宝兴大理岩高、低主频带分别集中在230~270 kHz(第24~27主频段)和0~50 kHz(第1~5主频段)。
将占比超过1%的主频段计入样本范围,统计得到干燥、饱和岩石的声发射高、低主频数量占比和能量占比,如图 9所示。干燥、饱和试验组试样主频数量占比均值和能量占比均值如表 2所示。其中能量以声发射信号幅值的平方表征。干燥岩石破坏的声发射以携带高能量的高主频带信号为主,最小高主频信号数量占比均值为69.51%,最小高主频信号能量占比均值为70.70%。3种硬岩饱和后,高主频带波形信号明显减少,低主频带波形信号显著增多。低主频带信号占据主导地位,且携带能量更多。最小低主频信号数量占比均值为58.77%,最小低主频信号能量占比均值为65.01%。
表 2 干燥、饱和岩石的声发射主频参数均值统计Table 2. Average values of dominant frequency parameters of different rocks under dry and saturated conditions岩石类型 状态 N/104 低频带 高频带 /% /% /% /% 宝兴大理岩 干燥 6.89 25.55 29.14 70.29 70.70 饱和 5.56 73.13 73.16 23.59 26.84 锦屏大理岩 干燥 2.42 28.63 25.15 70.84 74.35 饱和 1.31 58.77 65.01 35.18 29.11 灰岩 干燥 5.79 27.93 21.39 69.51 78.24 饱和 3.48 68.65 78.52 29.73 20.60 注:表中N,,分别表示波形数、不同主频带信号数量占比、不同主频带信号能量占比。 3. 水软化机制分析
硬岩的矿物成分决定水软化机制的存在条件。宝兴大理岩、锦屏大理岩和灰岩不含黏土矿物和硅质矿物,可能存在的水软化机制包括孔隙水压力(机制二)和摩擦弱化效应(机制三)。如表 2和图 9所示,上述3类硬岩饱和后,低主频信号数量及携带能量显著增加,并占绝对优势。李林芮[25]、邓建辉等[26]提出并验证,岩石的高主频波形信号对应微观剪切破坏,低主频波形信号对应微观拉伸破坏。换言之,本文研究的硬岩饱和后,发生的微观拉伸破坏明显增多,并占据主导地位。而3类岩石只可能存在孔隙水压力和摩擦弱化效应水软化机制,一种比较直观的推论是两种水软化机制中的一种或全部导致微观拉伸破坏数量的增加。
为论证上述推论,采用ANSYS数值软件对有(无)孔隙压力的单一孔隙周边的应力分布进行模拟分析。将试样孔隙简化为半径为1 mm的圆形孔洞,轴向压力为20 MPa。通过迭代计算获得孔洞周边的应力分布,如图 10所示。对比分析可以得出,有孔隙压力时,孔隙周边的拉应力普遍更大,最大拉应力由19 MPa(无孔隙压力时)增加至21.7 MPa,即孔隙压力导致拉应力增加。因此可以推断,饱和岩石的孔隙水压力作用机制将导致孔隙周边拉应力增加,产生更多微观拉伸破坏。这一推断也从图 6中饱和硬岩的破坏形态得到证实,即饱和硬岩出现更多张拉裂纹和张拉破裂面。而微观拉伸破坏对应声发射低主频波形信号,所以声发射低主频信号数量增加。总结来说,饱和岩石孔隙水压力作用的强弱与声发射低主频信号的多少存在对应关系。
3类硬岩饱和后低主频信号占比的增量与单轴抗压强度软化系数的关系如图 11所示。可以得到,低主频信号占比增量与强度软化系数间存在负相关关系,即低主频信号占比增量越大,强度降低越明显。这佐证了孔隙水压力是上述3种硬岩强度降低的主导因素。
相比较而言,宝兴大理岩和灰岩的低主频信号增量相差不多,强度软化系数也接近,而锦屏大理岩则表现出巨大差异。笔者认为这与岩石的矿物成分和孔隙结构密切相关。宝兴大理岩和灰岩的矿物成分均为方解石,孔隙度分别为6.16%和4.80%,两者矿物成分和孔隙结构相似,因而软化系数接近,分别为0.74和0.76。锦屏大理岩由多种矿物组成,孔隙度为2.68%。由图 1的SEM图像可以看出,矿物颗粒胶结紧密,未见明显孔隙。可以证实,锦屏大理岩结构致密,未提供孔隙水压力作用的良好条件,发生的微观拉伸破坏比宝兴大理岩和灰岩少,因而声发射低主频信号占比增量更小。综上分析,不含黏土、硅质矿物岩石的矿物成分和孔隙结构决定孔隙水压力作用的强弱。
4. 工程应用探讨
大多数岩体工程,如水电站库岸边坡、部分地下工程等,常处于水环境中,易引起岩体物理力学性能劣化,威胁工程的长期稳定和安全。在岩体工程各类分级评估中,包括工程岩体分级、岩体参数评估等,饱和单轴抗压强度作为强制性指标被广泛应用[27]。由于岩石饱和后,强度降低,采用饱和单轴抗压强度进行评估提高了工程的安全系数。但“全盘”要求使用饱和单轴抗压强度使工程稳定性评估有时过于保守,不利于建设的高效和经济。
不同类型饱和岩石的矿物成分、微结构、孔隙度和连通性等差异巨大,这决定了各种水软化机制的存在条件和作用程度。确定岩石矿物成分是水软化机制分析的首要工作。如前言所述,不含黏土矿物和硅质矿物的岩石可能存在孔隙水压力作用和摩擦弱化效应。本文分析得出,饱和碳酸钙类岩石的水软化机制以孔隙水压力为主导,产生更多的微观拉伸破坏,且孔隙水压力作用强弱受到孔隙度的影响。由于岩石抗压剪不抗拉伸[28],岩体工程分级评估时,应以提高安全系数为原则,即由大孔隙度、碳酸钙类矿物组成的岩体工程处于水环境中时,应严格按饱和单轴抗压强度进行设计评估,以保证工程安全。
声发射低主频信号增量与强度软化系数的负相关关系可应用于水环境中岩体工程的设计评估。在压缩荷载作用下,饱和岩石的声发射低主频信号增量与孔隙水压力作用强弱相对应,进而可预测发生微观拉伸破坏的多少。普遍认为,微观拉伸破坏与岩石强度间存在相关性。如图 12所示,笔者进行不同加载速率下宝兴大理岩干燥试样的单轴压缩试验,获得低主频信号(拉伸破坏)占比与单轴抗压强度的拟合曲线。发现两者间存在负相关性,即低主频信号越多,岩石强度越低。这与岩石抗压剪不抗拉伸的事实相符,由于岩石内部发生的微观拉伸破坏(低主频信号)多,岩石强度低。说明对应拉伸破坏的低主频信号可用于表征岩石强度变化。本文试验用3种岩石的强度软化系数随声发射低主频信号增量的增加而减小,说明声发射低主频信号可以表征碳酸钙类饱和岩石强度的劣化程度。声发射低主频信号可作为一种新的特征参数应用于含水岩体工程的设计和安全评估。
5. 结论
对3类硬岩(宝兴大理岩、锦屏大理岩、灰岩)饱和后的力学软化特性和声发射主频统计特征进行了分析,探究了水岩作用机制之孔隙水压力的存在条件和作用程度,取得以下3点结论。
(1)硬岩饱和后力学软化特性明显,单轴抗压强度和弹性模量减小,出现更多张拉裂纹和张拉破坏面。
(2)硬岩破坏释放的声发射信号存在双主频现象,且与岩石类型和含水状态无关。硬岩饱和后,声发射高主频带信号明显减少,低主频带信号显著增加。
(3)水岩作用机制之孔隙水压力是不含黏土、硅质矿物硬岩性能劣化的主导因素,其作用强弱与声发射低主频信号多少存在对应关系,且取决于矿物成分和孔隙结构。
-
表 1 试验土料基本物理性质
Table 1 Physical properties of test soil
液限 /% 塑限 /% 塑性指数 最优含水率 /% 最大干密度 /(g·cm-3) 2.70 34.0 21.0 13.0 13.7 1.88 表 2 试验方案
Table 2 Test schemes
方案编号 土样外径/mm 干密度/(g·cm-3) 结合面正应力/kPa 水力比降 剪切速率/(mm·d-1) 剪切变形及停顿点/mm T01 100 1.88 50 30 200 0-500-1000 T02 100 1.88 350 200 200 0-500-1000 T03 100 1.88 700 500 200 0-500-1000 T04 100 1.88 700 500 400 0-500-1000 T05 100 1.88 700 500 200 0-200-500-800-1000 T06 100 1.88 700 500 200 0-1000 T07 100 1.69 50 30 200 0-500-1000 T08 100 1.69 350 200 200 0-500-1000 T09 100 1.69 700 500 200 0-500-1000 T10 700 1.88 350 66 200 0-500-1000 T11 700 1.88 700 166 200 0-500-1000 -
[1] 邓刚, 丁勇, 张延亿, 等. 土质心墙土石坝沿革及体型和材料发展历程的回顾[J]. 中国水利水电科学研究院学报, 2021, 19(2): 1-13. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGSX202104006.htm DENG Gang, DING Yong, ZHANG Yan-yi, et al. Evolution of earth core embankment dams along with the development of configuration and material[J]. Journal of China Institute of Water Resources and Hydropower Research, 2021, 19(2): 1-13. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGSX202104006.htm
[2] KJAERNSLI B, TORBLAA I. Compaction of moraine in three feet layers[C]//Proceeding of the 7th International Congress on Large Dams, 1961, Rome.
[3] ASAO I. The Miboro Dam[C]//Proceeding of the 8th International Congress on Large Dams, 1964, Edinbergh.
[4] COOKE J B. Design Methods of Construction and Performance of High Rockfill Dams (above or about 80 m)[C]//Proceeding of the 8th International Congress on Large Dams, 1964, Edinburgh.
[5] ALBERRO J, MORENO E. Interaction phenomena in the Chicoasén dam: Construction and first filling[C]//Proceeding of the 14th International Congress on Large Dams, 1982, Rio De Janeiro.
[6] 邓刚, 韩巍巍, 温彦锋, 等. Hyttejuvet坝突然渗漏事故的回顾和心墙堆石坝水力劈裂的反思[C]//土石坝技术2012年论文集, 2012, 北京. DENG Gang, HAN Wei-wei, WEN Yan-feng, et al. A review of the sudden leakage of hyttejuvet dam and a reflection on the hydraulic fracturing of the earth core rockfill dams[C]//Proceeding of the Symposium on Technology of Earth-Rockfill Dam Rockfill Dams, 2012, Beijing. (in Chinese)
[7] 邓刚, 皇甫泽华, 武颖利, 等. 土质心墙土石坝变形协调控制发展与展望[J]. 水力发电学报, 2020, 39(5): 1-16. DENG Gang, HUANGFU Ze-hua, WU Yin-li, et al. Development and prospect of deformation compatibility control of earth core embankment dams[J]. Journal of Hydroelectric Engineering, 2020, 39(5): 1-16. (in Chinese)
[8] 邓刚, 陈辉, 张茵琪, 等. 基于坝壳湿化过程数值模拟的心墙坝初蓄水力劈裂机理研究[J]. 中国水利水电科学研究院学报, 2021, 19(1): 90-98. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGSX202101012.htm DENG Gang, CHEN Hui, ZHANG Yin-qi, et al. Study on mechanism of hydraulic fracture of earth core embankment dam during first filling based on the numerical simulation of wetting deformation development with time of coarse aggregate[J]. Journal of China Institute of Water Resources and Hydropower Research, 2021, 19(1): 90-98. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZGSX202101012.htm
[9] DEWHURST D N, CLENNELL M B, BROWN K M, et al. Fabric and hydraulic conductivity of sheared clays[J]. Géotechnique, 1996, 46(4): 761-768. doi: 10.1680/geot.1996.46.4.761
[10] DEWHURST D N, BROWN K M, CLENNELL M B, et al. A comparison of the fabric and permeability anistrophy of consolidated and sheared silty clay[J]. Engineering Geology, 1996, 42(4): 253-267. doi: 10.1016/0013-7952(95)00089-5
[11] ZHANG S Q, TERRY E T. The effect of fault slip on permeability and permeability anisotropy in quartz gouge[J]. Tectonophysics, 1998, 295: 41-52. doi: 10.1016/S0040-1951(98)00114-0
[12] KIMURA S, KANEKO H, ITO T, et al. Investigation of fault permeability in sands with different mineral compositions (evaluation of gas hydrate reservoir)[J]. Energies, 2015, 8: 7202-7223. doi: 10.3390/en8077202
[13] IMURA S, KANEKO H, NODA S, et al. Shear-induced permeability reduction and shear-zone development of sand under high vertical stress[J]. Engineering Geology, 2018, 238: 86-98. doi: 10.1016/j.enggeo.2018.02.018
[14] WANG G, WEI X, ZOU T. A hollow cylinder radial-seepage apparatus for evaluating permeability of sheared compacted clay[J]. Geotechnical Testing Journal, 2019, 42(5): 1133-1149.
[15] 魏星, 邹婷, 王刚. 压-剪耦合条件下黏土渗透特性的试验研究.[J]. 岩石力学与工程学报, 2017, 36(增刊1): 3561-3568. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2017S1050.htm WEI Xing, ZOU Ting, WANG Gang. Experimental study on permeability of clay during coupled compression and shear[J]. Chinese Journal of Rock Mechanics and Geotechnical Engineering, 2017, 36(S1): 3561-3568. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2017S1050.htm
[16] 王刚, 游克勤, 魏星, 等. 压实黏土剪切带渗透特性试验研究[J]. 岩土工程学报, 2019, 41(8): 1530-1537. WANG Gang, YOU Ke-qin, WEI Xing, et al. Experimental study on permeability of shear bands in compacted clay[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(8): 1530-1537. (in Chinese)
[17] 雷红军, 卞锋, 于玉贞, 等. 黏土大剪切变形中的渗透特性试验研究[J]. 岩土力学, 2010, 31(4): 1130-1133. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201004025.htm LEI Hong-jun, BIAN Feng, YU Yu-zhen, et al. Experimental study of permeability of clayey soil during process of large shear deformation[J]. Rock and Soil Mechanics, 2010, 31(4): 1130-1133. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201004025.htm
[18] LEI H J, WU Y K, YU Y Z, et al. Influence of shear on permeability of clayey soil[J]. International Journal of Geomechanics, 2016, 16(5): 04016010
[19] LIU Q H, WU Y K, LI Q M, et al. Modified model for hydraulic conductivity of clayey soil under shear[J]. International Journal of Geomechanics, 2019, 19(11): 06019015.
[20] 王刚, 韦林邑, 魏星, 等. 压实黏土三轴压缩变形过程中的渗透性变化规律[J]. 岩土力学, 2020, 41(1): 32-38. WANG Gang, WEI Lin-yi, WEI Xing, et al. Permeability evolution of compacted clay during triaxial compression[J]. Rock and Soil Mechanics, 2020, 41(1): 32-38. (in Chinese)
[21] 雷红军, 刘中阁, 于玉贞, 等. 黏土-结构接触面大剪切变形后渗流特性试验研究[J]. 岩土力学, 2011, 32(4): 1040-1044. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201104013.htm LEI Hong-jun, LIU Zhong-ge, YU Yu-zhen, et al. Experimental study of seepage characteristics of clayey soil-structure interface under large shear deformation[J]. Rock and Soil Mechanics, 2011, 32(4): 1040-1044. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201104013.htm
[22] LUO Y L, JIN X, LI X, et al. A new apparatus for evaluation of contact erosion at the soil-structure interface[J]. Geotechnical Testing Journal, 2013, 36(2): 256-263.
[23] 邓刚, 湛正刚, 张幸幸, 等. 一种旋转剪切式接触面抗渗特性试验装置:中国,201610305992.8[P]. 2019-01-25. DENG Gang, ZHAN Zheng-gang, ZHANG Xing-xing, et al. The invention discloses a rotary shear type test device for impermeability of contact surface: China, 201610305992.8[P]. 2019-01-25. (in Chinese)
[24] 邓刚, 湛正刚, 张茵琪, 等. 一种模拟接触面渗流耦合特性的试验装置及试验方法: CN111912760A[P]. 2020-11-10. DENG Gang, ZHAN Zheng-gang, ZHANG Yin-qi, et al. The invention relates to a test device and a test method for simulating the seepage coupling characteristics of contact surfaces: CN111912760A[P]. 2020-11-10. (in Chinese)
[25] COLLINS K, MCGOWN A. The form and function of microfabric features in a variety of natural soils[J]. Géotechnique, 1974, 24(2): 223-254.
[26] MITCHELL J K, SOGA K. Fundamental of soil behavior[M]. 3rd ed. New Jersey: John Wiley and Sons, 2005.
[27] OLSEN H W. Hydraulic flow through saturated clay[C]//Proceedings of the Ninth National Conference on Clays and Clay Minerals, 1962, West Lafayette.
[28] DELAGE P, LEFEBVRE G. Study of the structure of a sensitive Champlain clay and of its evolution during consolidation[J]. Canadian Geotechnical Journal, 1984, 21(1): 21-35.
[29] MORGENSTERN N R, TCHALENKO J S. Microscopic structures in Kaolin subjected to direct shear[J]. Géotechnique, 1967, 17: 309-328.