0.   引言

石灰岩是地壳中分布最广的岩体之一，广泛存在于中国水利和岩土工程中。水库蓄、放水、地下工程开挖等过程产生的加载、卸荷效应对裂隙岩体渗流的影响不容忽视^[1-2]。因裂隙岩体渗流问题发生大坝失事案例不少^[3]，如法国Malpasset双曲拱坝、Vajoint拱坝及中国梅山连拱坝等^[3]。研究加卸载作用下石灰岩粗糙单裂隙渗流规律具有重要理论意义。

Snow等^[4]通过光滑平行板裂隙渗流试验，提出了单裂隙渗流立方定律。据此，国内外学者^[5-8]针对粗糙裂隙，提出了许多立方定律修正公式，但其大多没有考虑裂隙外部荷载变化的影响，描述加卸载作用下粗糙单裂隙渗流特性时存在局限。石灰岩裂隙岩体受地震、构造作用影响，易形成加卸载应力环境，加载、卸荷过程对岩体裂隙渗流影响有所不同。梁宁慧等^[9]指出卸荷过程岩体裂隙渗透系数变化存在滞后性，验证了卸荷、加载过程裂隙渗透系数变化的差异性，提出了卸荷量与渗透系数的关系式。孙可明等^[10]开展了不同粗糙裂隙岩体加卸载渗流试验，分析试验前后裂隙面形貌特征，揭示了粗糙裂隙渗流加载、卸荷作用下渗透率、平均宽度等参数变化规律。杨秀荣等^[11]研究了渗流-应力耦合作用下石灰岩渗透破坏模式，揭示了破裂岩样渗透率随净围压、渗压加卸载变化规律。Tian等^[12]等以细砂岩、砂岩和复合岩体为试验岩样，研究了裂隙水压力和围压对裂隙渗透系数的影响，表明裂隙渗透系数随裂隙水压、净围压增加，分别呈线性增加和指数下降的变化趋势。Zhao等^[13]研究了3种不同裂隙石灰岩岩样在应力加载、卸荷环境下的裂隙渗流特性和形貌特征变化规律，提出了一种考虑裂隙面形貌特征和弹塑性变形影响的渗透率分析模型。Fu等^[14]研究了围压作用下砂质泥岩粗糙裂隙渗透性变化规律，指出裂隙渗透系数受裂隙面粗糙度和围压共同影响，其中围压加载的影响最显著，二者呈指数函数关系。Zhang等^[15]发现围压首次加卸载过程对裂隙岩体渗透率影响最显著，围压卸荷过程中，渗透率恢复存在滞后性。Chao等^[16]基于围压循环加卸载条件下渗透试验成果，提出了裂隙渗透张量的数学计算模型。

综上研究表明，考虑裂隙面粗糙度、加卸载作用对岩体裂隙渗流特性的研究成果已有不少，但同时针对不同节理粗糙度裂隙面、研究围压和渗压加卸载作用下石灰岩单粗糙裂隙渗流规律的成果较少，多因素影响下裂隙渗流特性演化机理尚不清楚。鉴于此，选定5组不同JRC石灰岩裂隙试样，开展围压、渗压加卸载条件下石灰岩单裂隙试样渗流试验，揭示围压、渗压加卸载过程中粗糙单裂隙渗流特性演化规律。

1.   JRC和渗透率计算

裂隙面粗糙程度对裂隙岩体渗流特性影响显著，研究粗糙裂隙渗流特性时，需先明确裂隙面形貌特征表征形式^[17-19]。挪威学者Barton^[5]于1973年提出了节理粗糙度系数JRC的概念。Barton等^[20]通过136个岩样剪切试验，绘制了10条标准轮廓曲线，明确了JRC取值范围为0~20，提出了基于10条标准轮廓曲线的JRC计算方法，受到国际岩石力学学会的认可。诸多学者^[21-24]结合先进测量技术和方法，从维度、表征参数和计算公式等方面进行裂隙粗糙程度表征，其中Tse^[21]提出的JRC计算公式应用广泛。Blelem等^[22]提出平均梯度模的Z_2s概念。王珂等^[25]将Tse公式的表征维度提升至三维，本文采用该方法^[25]进行JRC计算，具体公式如下：

式中：${Z_{{\text{2S}}}}$为平均梯度模；N_x，N_y分别为X轴，Y轴方向上采样点的数量；∆x，∆y分别为X轴，Y轴方向上采样点的间隔；Z_i，Z_i+1分别为第i和第i+1裂隙离散点的Z轴方向坐标，i，j分别为沿X轴和Y轴方向上细观平面的序号。

由于岩体基质渗透性远小于裂隙，裂隙岩体渗透性主要由裂隙渗透率决定。本文采取立方定律反算裂隙渗透率^[26]，计算公式具体如下：

式中：Q为通过裂隙面的渗流量（m³/s）；$A$为裂隙面的横断面面积（m²）；$\upsilon$为通过裂隙面流体的流速（m/s）；$K$为裂隙面的渗透系数（m/s）；$D$为裂隙试样的直径（m）；${b_{\text{h}}}$为裂隙的等效水力隙宽（m）；$\Delta P$为裂隙试样的两端压差（Pa）；$\rho$为水流的密度（kg/m³）；$L$为裂隙试样的长度（m）；$\mu$为水流的动力黏滞系数（Pa·s）；$k$为石灰岩裂隙面的渗透率（m²）；$\gamma$为水流的重度（N/m³）。

2.   试验概况

2.1   试样制备

石灰岩试样取自济宁嘉祥纸坊镇矿山，质地较密，其基本物理性质如表 1所示。采取水钻法在完整岩块上钻取直径50 mm，高100 mm的圆柱体试样。采用微机控制电液伺服万能试验机（WAW-1000）劈裂圆柱体试样^[27]，设定不同的加载速率获取20组不同JRC劈裂试样，通过南京Wiiboox公司提供的Reeyee 3M三维光学扫描仪采集裂隙表面几何特征，将扫描数据导入自主编写的JRC计算程序中计算，选取5个典型JRC裂隙面试样进行渗流试验，典型JRC裂隙面单侧壁面数字化图像如图 1所示。

表  1  石灰岩物理性质

Table  1.  Physical properties of limestone

物理量	计量单位	平均值
密度	g/cm3	2.70
吸水率	%	0.21
抗压强度	MPa	87.50
抗拉强度	MPa	8.20
抗弯强度	MPa	9.11

下载: 导出CSV 

| 显示表格

图  1  不同JRC的裂隙面数字化图像

Figure  1.  Digital images of fracture surfaces with different JRCs

下载: 全尺寸图片 幻灯片

为保证围压和渗压加卸载试验中试样裂隙面JRC统一，以劈裂岩样裂隙面为模板，采用激光雕刻技术复刻形貌特征与之相同的裂隙岩样。试样雕刻步骤如下：①对5个不同JRC裂隙岩样劈裂面进行扫描，获得坐标点云数据；②将数据导入岩石雕刻设备所依托的维宏软件重构劈裂面模型；③依据重构模型，对石灰岩进行雕刻；④将雕刻好的两半裂隙面对齐拼接，钻芯取样。将雕刻所得两半岩样贴紧，涂抹704硅胶密封侧壁，以避免试验过程中发生侧壁渗漏，图 2为经过劈裂、雕刻、胶装处理后的试样形态图。

图  2  制备阶段的试样形态图

Figure  2.  Shapes of samples during preparation process

下载: 全尺寸图片 幻灯片

2.2   试验装置

渗流试验采用自主研制的岩体裂隙渗流试验装置，包括轴向、围压和渗压加卸载装置，流量、轴压、围压和渗压数据采集系统等，试验系统主要装置见图 3。

图  3  试验系统主要装置图

Figure  3.  Diagram of main device of test system

下载: 全尺寸图片 幻灯片

2.3   试验过程和方案

选取不同JRC（2.64，7.78，11.40，15.47，19.51）石灰石单裂隙试样进行渗流试验，试验步骤如下：①将试样放置于压力室上、下压头中间，施加轴压固定；②烘紧热缩膜密封试样表面及上、下压头接触部分；③升高保温仓密封压力室，打开排气阀，待仓内满油后关闭；④打开温控系统，维持舱内温度恒定；⑤预加围压2.0 MPa后，施加渗压，转应力控制模式维持稳定，完成试验前准备工作。操作过程中施加初始轴压5.0 MPa固定试样，保持稳定，开展各工况渗透试验时，轴向采用位移伺服控制模式，轴压随围压与渗压协调变化，保持轴压等于净围压。

围压和渗压加卸载过程分为10级，以围压加卸载试验为例，试验开始前设置围压为2 MPa，启动试验，保持2 MPa围压2 h后以0.005 MPa/s速率加压至下一级，继续保持2 h，依次类推，直至增加至10 MPa后以相同速率卸压，降至2 MPa后保持2 h，记录各级压力下的数据后结束试验。渗压加卸载过程与围压相同，起始压力为0.2 MPa，加载速率为0.001 MPa，最大渗透压设置为1.5 MPa。分级加卸载试验方案见表 2。

表  2  不同加卸载条件下渗流试验方案

Table  2.  Schemes of seepage test under different loading and unloading conditions

试样编号	JRC	围压/MPa	渗压/MPa	加载速率/(MPa·s-1)
M-1	2.64	2，3，4，6，8，10，8，6，4，3，2	0.5	0.005
M-2	7.78
M-3	11.40
M-4	15.47
M-5	19.51
H-1	2.64	3	0.2，0.5，0.8，1.0，1.2，1.5，1.2，1.0，0.8，0.5，0.2	0.001
H-2	7.78
H-3	11.40
H-4	15.47
H-5	19.51

下载: 导出CSV 

| 显示表格

3.   结果分析

3.1   围压加卸载作用的影响

图 4为石灰岩裂隙渗流量随围压加卸载变化曲线。

图  4  渗流量随围压变化曲线

Figure  4.  Curves of seepage flow with confining pressure

下载: 全尺寸图片 幻灯片

由图 4可知，不同JRC裂隙渗流量随围压增大逐渐下降，下降率逐渐减小并趋于稳定，JRC越小，初始渗流量越大，下降量越大。低围压条件下，裂隙平均开度较大，裂隙面接触面积较小，渗流量大；围压加载至4.0 MPa时，岩样M-1~M-5渗流量显著下降，下降量占总下降量比值相近，最小为70.01%。这是由于侧向围压作用下，试样裂隙被压密，隙宽减小，裂隙左、右粗糙面接触点局部应力增大，超过岩体强度而发生变形，转为面接触或相互嵌入，接触面积增大，对流面积减小，渗流量下降；随着围压加大，渗流量持续下降，但下降量明显减少，这是因为裂隙面变化已趋于稳定，裂隙接触面应力不均匀性减小，隙宽和对流面积已达到下降阈值。

围压卸载过程中，不同JRC裂隙渗流量变化呈现先减小后增大的趋势，渗流量回升存在滞后性^[9]；围压降至2.0 MPa时，渗流量和试验初始值存在较大差距，说明在围压作用下，岩样裂隙变化规律符合弹塑性变形理论^[13]。卸荷过程中，裂隙隙宽在水压力作用下有所增大，接触面变形恢复，但存在不可逆转的塑性损伤，导致相同围压下裂隙渗透性较之加载时显著下降。围压加卸载过程，裂隙面JRC越大，相同围压下渗流量越小，这是由于JRC越大，裂隙面凹凸起伏度越大，凹凸体数量偏多，渗液在裂隙内流动路径延长，导致渗流量减小^[27]。

从裂隙渗流演化机理上分析，在发展期，压力作用下裂隙面凸起会嵌入另一面，接触面增大，隙宽减少，过流通道减小；此外，裂隙内破碎颗粒增多，阻碍了水流在裂隙内流动。因此，初期渗流量下降最为迅速。围压徐变作用下，裂隙隙宽持续减小，左、右裂隙面贴合更紧密。进入过渡期，裂隙左、右裂隙面由点接触变为面接触，凹凸体破损和变形过程减缓，移动颗粒因间隙变小被逐渐固定，渗流量下降趋势减缓。当裂隙面接触面积达到一定程度，裂隙面之间的反作用力将平衡围压施加的作用力，此时裂隙面凹凸体变形和隙宽稳定，渗流量不再下降，裂隙渗流进入稳定期。在整个裂隙渗流过程中，同时期裂隙面JRC越小，渗流量越大，这是因为过流面积和路径在裂隙渗流中起主要作用。JRC越小，裂隙面凹凸体数量越少，凹凸体接触面积越小，渗液在裂隙面过流路径越短，过流面积较大。因此，围压作用下，JRC越小的裂隙，其渗流量随时间不断减小，但其值仍高于同时期JRC大的裂隙。

依据式（3）~（5）计算裂隙渗透率，绘制裂隙渗透率随围压加卸载的变化曲线，如图 5所示。结合图 4，5可知，裂隙渗透率变化规律与渗流量一致。相同试样及试验条件下，单位时间内裂隙渗流量与k^3/2成正比，故在后续试验结果分析中，选取渗透率k进行详细陈述。

图  5  渗透率与围压变化曲线

Figure  5.  Curves of permeability and confining pressure

下载: 全尺寸图片 幻灯片

为定量分析围压加载过程裂隙渗透率的变化规律，对岩样M-1~M-5各级围压渗透率进行拟合，拟合关系曲线如图 6所示。由图 6可知，围压加载过程中，裂隙渗透率和围压呈良好的幂函数关系，不同JRC岩样裂隙渗透率拟合公式的R²最小为0.9583，拟合公式如下：

图  6  渗透率与围压P_a的拟合关系曲线

Figure  6.  Fitting curves between permeability and confining pressure

下载: 全尺寸图片 幻灯片

式中：P_a为围压（MPa）；a，b为与JRC相关的试验参数。

选取不同JRC岩样对a，b参数值进行回归分析，拟合曲线如图 7所示，拟合公式如下：

图  7  试验参数a，b与JRC的拟合关系曲线

Figure  7.  Fitting curves of test parameters a, b and JRC

下载: 全尺寸图片 幻灯片

围压卸荷过程裂隙渗透率与围压不存在明显相关性。卸荷过程裂隙渗透率变化存在明显的滞回效应，围压减小初期，裂隙渗透率维持稳定或轻微下降，当围压卸载到一定阈值后，渗透率开始回升。不同JRC裂隙渗透率卸荷过程的滞后程度不同，JRC越极端，滞后程度越小，围压阈值越大。如岩样M-1渗透率回升阈值为8.0 MPa，岩样M-3围压阈值仅为3.0 MPa，为直观分析JRC对卸荷过程裂隙渗透率滞回效应的影响，提取滞回阶段裂隙渗透率下降量k_l和回升量k_f，k_l和k_f随JRC的拟合关系曲线如图 8所示。由图 8可知，滞回阶段裂隙渗透率下降量和回升量随JRC增大而减小，在2.0~19.0的JRC范围内，k_l和k_f与JRC的拟合公式R²分别为0.8633和0.9805，公式如下：

图  8  k_l和k_f随JRC变化的关系曲线

Figure  8.  Curves of k_l and k_f with JRC

下载: 全尺寸图片 幻灯片

综上所述，JRC对卸荷过程裂隙渗透率的影响主要体现在滞回阶段渗透率下降回升程度及回升的围压阈值；JRC越大，渗透率下降和回升程度越明显。

3.2   渗压加卸载作用的影响

图 9为不同JRC裂隙渗流量随渗透压加卸载的变化曲线。由图 9可知，不同JRC岩样裂隙的初始渗流量接近，最大差值为0.98×10^-8 m³/s，裂隙面JRC与渗流量之间不存在明显对应关系；渗压加卸载与围压加卸载渗透试验初期渗流量与JRC之间的相关性受裂隙面净围压影响较大。由表 2可知，渗压加卸载条件下的渗流试验初始净围压值为2.8 MPa，远大于围压加卸载条件下的1.5 MPa；净围压增大导致裂隙隙宽减小，裂隙壁面贴合更紧密，渗流量下降。此时不同JRC裂隙面渗流量均处于较低水平，JRC对渗流量影响不明显。渗流量随渗压增大而上升，在0.2 MPa至0.8 MPa加载阶段，各岩样裂隙渗流量上升较为显著，增长量占整体增量的69.65%~78.99%，此外，相同渗压下裂隙渗流量随JRC增大而减小，其规律可与围压加卸载试验结果相互佐证，说明在0.5 MPa至1.5 MPa渗压和2.0 MPa至10.0 MPa围压变化范围内，裂隙渗流量与JRC负相关。

图  9  裂隙渗流量随渗压加卸载变化曲线

Figure  9.  Curves of seepage flow of fracture with loading and unloading of seepage pressure

下载: 全尺寸图片 幻灯片

渗压卸荷过程，渗流量逐渐减小，待降至初始值0.2 MPa时，岩样H-1~H-5裂隙渗流量相近，最大差值为0.46×10^-8 m³/s，差值较小；与初始值相比，各岩样裂隙渗流量下降了2.60×10^-8~3.10×10^-8 m³/s，降幅较大。为研究初次渗压加卸载后裂隙渗流量与JRC的变化关系，绘制渗流量下降量$\Delta Q$随JRC变化关系曲线如图 10所示。由图 10可知，$\Delta Q$随JRC增大呈现先下降后增大的变化趋势，建立二次函数拟合关系，其R²为0.9116，具体公式如下：

图  10  JRC与渗流量下降量拟合关系

Figure  10.  Fitting relationship between JRC and decrease of seepage discharge

下载: 全尺寸图片 幻灯片

初始渗流量${Q_0}$和裂隙面JRC也呈二次函数关系，拟合结果见图 11所示，其R²为0.813，具体公式如下：

图  11  JRC与初始渗流量的拟合关系

Figure  11.  Fitting relationship between JRC and initial seepage rate

下载: 全尺寸图片 幻灯片

将监测所得初始渗流量代入式（12）反算裂隙JRC，再将JRC代入式（11），即可得到0.2~1.6 MPa范围内渗压加卸载后裂隙$\Delta Q$。

为直观分析渗压加卸载过程裂隙岩体渗透性变化规律，绘制裂隙渗透率随渗压加卸载的变化曲线图 12。由图 12可知，渗压加载过程，不同JRC裂隙渗透率逐渐降低，与图 9呈现的渗流量变化趋势相反，在3.0 MPa围压和0.2~1.6 MPa渗压范围内，随着渗压加载，裂隙内水压力增大，渗透水流速增大，冲蚀作用增强，裂隙面岩体强度下降，在围压挤压和渗透水冲蚀作用下，裂隙内壁薄弱处会发生二次破坏 ^[27]；此外，碳酸钙等矿物质易溶解于渗液而析出，渗流速度增大会促进析出过程，凹凸体数量及起伏度增加，渗流路径延长、左右壁面接触面积增大；因此在应力、冲蚀和溶解等作用影响下，渗流量随渗压增加而增大的程度被削弱，裂隙等效水力隙宽和渗透率呈下降趋势。

图  12  裂隙渗透率随渗压加卸载变化曲线

Figure  12.  Curves of fracture permeability with seepage pressure under loading and unloading conditions

下载: 全尺寸图片 幻灯片

结合式（3）~（5）可进一步印证渗透率变化规律，渗压加载过程中，渗流量Q增加（图 9），但其增幅小于渗透压力差$\Delta P$的增幅，故Q/$\Delta P$比值呈下降趋势；裂隙渗透率与Q/$\Delta P$正相关，故渗压加载过程中，裂隙渗透率亦呈下降趋势。在0.2~0.5 MPa的渗压加载初期，JRC越大，渗透率下降越显著，渗透率大小与JRC相关性不显著；在0.5~1.5 MPa的渗压加载中后期，不同JRC裂隙渗透率下降率波动较大，相同渗压下，裂隙渗透率随JRC增大而减小。渗压卸荷过程，不同JRC裂隙渗透率变化趋势与渗流量相似，呈下降趋势；说明裂隙内水压力对水力隙宽减小存在延缓作用，故随着渗压减小，裂隙内水压力下降，隙宽减小，裂隙渗透率逐级下降。综上可知，在3.0 MPa围压和0.2~1.6 MPa渗压加卸载过程中，渗压变化对裂隙面渗透率主要起削减作用。

3.3   形貌特征分析

围压和渗压加卸载过程中，作用于裂隙面上的应力作用和冲蚀作用均会影响裂隙面形貌特征。图 13是JRC值为11.40的裂隙面渗透试验前后形貌特征对比图，其中左图为裂隙面实物图，右图是其凸起高程分布图（以裂隙面最低点为基准参考面）。

图  13  试验前后试样裂隙（JRC = 11.40）形貌特征图

Figure  13.  Characteristics of morphologies of fracture surface of samples with JRC of 11.40 before and after tests

下载: 全尺寸图片 幻灯片

由图 13可知，裂隙表面存在因凹凸体嵌入产生的深色点状印记，对比试验前后裂隙面形貌特征可知，渗流试验后点状印记数量显著增多，相较于试验前，分布更加均匀，说明围压作用下，壁面贴合更加紧密，壁面接触面积增大，相互挤压的凹凸体数量增多、嵌入程度增加；历经围压加卸载或渗压加卸载的渗透试验后，裂隙表面凸起高度均减小，前者降幅较大，说明作用于裂隙壁面的法向应力和渗流作用对裂隙面JRC起削弱作用，但随着试验前裂隙面JRC减小，削弱效果逐渐下降，当裂隙面JRC为2.64时，出现增强效果，其原因为裂隙中流动的偏中性水溶液对裂隙表面矿物质溶解作用强度较低。JRC较大时，裂隙表面凹凸起伏程度大，凹凸体数量偏多，此时应力和渗流冲蚀产生的削弱作用占据主导；而JRC较小时，裂隙表面整体趋于平滑，此时溶解作用对裂隙面凹凸体数量及起伏程度的增强作用占据主导，试验后JRC有所增加。

4.   讨论

围压加载过程，不同JRC裂隙渗流量和渗透率均与围压呈负相关，卸荷过程存在滞回阶段，其程度与JRC相关，相同围压下，JRC越大，渗流量和渗透率越大；裂隙渗流量和渗透率随围压变化趋势与孙可明等^[10]研究成果基本一致，但由于试验岩样不同，孙可明等^[10]试验得出的渗透率和流量在法向荷载卸荷过程中未明显下降。渗压加卸载过程中，不同JRC裂隙渗流量与渗压正相关，渗透率则呈持续降低趋势。在3.0 MPa围压和0.2~1.6 MPa渗压作用下，渗压增大对石灰岩裂隙渗透率的影响以削弱为主。低渗压（0.2~0.5 MPa）下裂隙渗流量、渗透率与JRC相关性不强。本文渗压加卸载试验成果仅适用于单一贯通裂隙或断层的渗流分析，当石灰岩破裂，裂隙发展成网状分布时，这一变化规律适用性下降^[11]。围压或渗压初次加卸载完成后，裂隙渗流量和渗透率均有所降低，说明围压和渗压加卸载过程会导致裂隙粗糙面形貌特征发生不可扭转的塑性变形。

为避免岩体内部裂隙网络，尤其是贯穿裂隙或断层造成基岩发生大面积渗漏，在大坝选址前的地质勘探过程中，宜对易发生大面积渗漏的裂隙或断层进行渗流量监测和应力分析，获得其初始围压与渗流量。通过式（3）~（8）反算裂隙初始JRC，结合该区域裂隙岩体渗透率设计标准，可估算围压增量；结合式（9），（10）可估算50 m水头下2.0~10.0 MPa范围内围压首次加卸载作用对岩体裂隙渗透率的削减程度。联立式（11），（12），基于当前渗流量监测值，估算对应位置裂隙面JRC，并推算历经单次0.2~1.6 MPa幅度渗压加卸载作用后的裂隙渗流量，在监测点故障时可作为补充数据；此外，通过对比分析与实际监测数据差异，可提前预知相关裂隙所处渗压环境变化。

由于实际工程中，岩体裂隙的应力环境极为复杂，围压、渗压等变化幅度及频率随机性大，天然裂隙或断层的倾角和尺寸与实验室样品存在较大差异，将本文提出的经验公式用于首次加卸载后渗流参数预测计算时，应结合现场压水试验成果对经验公式进行修正。且本文试验围压比渗压至少大1.5 MPa，试验成果未能体现渗压加载接近围压或超过围压后的粗糙裂隙渗流特性，也未考虑轴压、岩性、裂隙开度、充填物、温度和渗透介质酸碱性等诸多因素对裂隙渗流规律的影响。故本文提出的裂隙渗流量及渗透率计算公式在多因素交叉影响下的适用性有待进一步论证。

5.   结论

（1）石灰岩裂隙渗流量和渗透率随围压增大而增大；围压加载时不同JRC裂隙渗透率均与围压呈幂函数关系，R²最小值为0.9583；围压卸荷阶段裂隙渗透率变化呈现先下降后回升的滞回效应；不同JRC裂隙渗透率与围压无明确的定量关系。

（2）石灰岩裂隙渗流量和渗透率随渗压加载、卸荷的变化趋势存在差异，裂隙渗流量变化与渗压正相关，渗透率呈下降趋势；初始渗流量及渗压加载、卸荷后渗流量下降量与裂隙JRC均呈二次函数关系。

（3）2.0~10 MPa围压加卸载范围内，JRC越小，渗流量和渗透率越大；围压卸荷时渗透率下降量和回升量均与JRC相关，JRC越大，卸荷时渗透率滞回效应越明显。在0.5~1.5 MPa范围内的渗压加卸载作用下，JRC越大，渗流量和渗透率越大；渗压降至0.5 MPa以下，渗流量和渗透率与JRC无明显的相关性。

（4）本文提出了一种围压加卸载作用下粗糙裂隙渗透率估算模型及渗压加卸载作用下渗流量预测模型，建立了裂隙渗透率与JRC、围压和渗压等参数之间的定量关系，研究成果可为加卸载环境下岩体裂隙渗透参数和渗漏量估算提供理论依据。