水-岩及风化作用下层状炭质板岩巴西劈裂试验研究

李二强; 冯吉利; 张龙飞; 张洪昌; 朱天宇

doi:10.11779/CJGE202102013

水-岩及风化作用下层状炭质板岩巴西劈裂试验研究 English Version

李二强^{1, 2,},
冯吉利^{1, 2, ,},
张龙飞^{1, 2},
张洪昌^{1, 2},
朱天宇^{1, 2}

1.
深部岩土力学与地下工程国家重点实验室,北京　100083
2.
中国矿业大学(北京）力学与建筑工程学院,北京　100083

基金项目:

国家重点研发计划项目 2016YFC0600901

国家自然科学基金项目 41172116

国家自然科学基金项目 U1261212

详细信息

作者简介:
李二强（1987— ）,男,博士研究生,主要从事岩石力学与地下工程的研究工作。E-mail：lerqiang@163.com

通讯作者:
冯吉利, E-mail: fjl@cumtb.edu.cn

中图分类号: TU452
计量
- 文章访问数: 259
- HTML全文浏览量: 20
- PDF下载量: 139
出版历程
- 收稿日期: 2020-04-20
- 网络出版日期: 2022-12-04
- 刊出日期: 2021-01-31

Brazilian tests on layered carbonaceous slate under water-rock interaction and weathering

LI Er-qiang^{1, 2,},
FENG Ji-li^{1, 2, ,},
ZHANG Long-fei^{1, 2},
ZHANG Hong-chang^{1, 2},
ZHU Tian-yu^{1, 2}

1.
State Key Laboratory for Geomechanics and Deep Underground Engineering, Beijing 100083, China
2.
School of Mechanics and Civil Engineering, China University of Mining and Technology Beijing, Beijing 100083, China

摘要

摘要: 横观各向同性层状板岩在地质环境及施工裸露后形成的水-岩和自然风化作用下,其物理力学性质将发生退化,由此常引起岩土隧道、边坡等工程产生抗拉失稳破坏。以甘肃岷县木寨岭隧道炭质板岩为研究对象,对其开展即时烘干试样（即试样加工后烘干处理）和静置风化试样（加工后室内常温通风放置60 d）巴西劈裂试验,分析了两种条件下炭质板岩的抗拉力学性能。结果表明：两种炭质板岩巴西劈裂试验均表现为脆性破坏,但力学响应及峰值荷载差异极大；炭质板岩抗拉强度均受层理影响显著,但即时烘干试样和静置风化试样的劈裂破坏形态受层理影响程度不同；以即时烘干试样为基础,静置风化试样软化系数为0.11～0.13,体现了水岩作用及其静置风化对岩样具有显著劣化效应,但对层理形成的横观各向同性抗拉性能影响不大。此外,借助SEM对水岩静置风化机理进行分析,具体为水侵入后形成物理化学反应改变矿物微观成分、结构进而结合风化作用表现为矿物骨架承载的累计劣化损伤。
- 炭质板岩 /
- 层理 /
- 水-岩作用 /
- 风化 /
- 抗拉强度 /
- 破坏模式 /
- 劣化机理
Abstract: The physical-mechanical properties of transversely isotropic layered slates may be degraded due to water-rock interaction (WRI) and natural weathering, which potentially leads to the instability or collapse of tunnels, slopes and mining. In this study, based on the specimens sampled from the carbonaceous slates of Muzhailing Tunnel, the corresponding drying samples (drying in oven after fabrication) and weathering samples (for 60 days’ natural weathering after fabrication) are prepared, respectively. The Brazilian tensile tests are conducted. The test results show that the two kinds of carbonaceous slates are both characterized by brittle failure, but the mechanical response covering failure displacement and the peak load are rather different. Further, the tensile strength are significantly affected by bedding while the cleavage failure patterns are also dominated by bedding. The softening coefficient of the static weathering samples is 0.11~0.13, which implies that the WRI and natural weathering play a vital role in the course of rock failure, but they have few influences on the transverse isotropic tensile properties of bedding. Moreover, the mechanisms of specimen failure are explained by the SEM technique to analyze the micro-components and observe the process of specimen deterioration due to physico-chemical reaction.
- carbonaceous slate /
- bedding /
- water-rock interaction /
- weathering /
- tensile strength /
- failure pattern /
- degradation mechanism

HTML全文

0. 引言

持水特征曲线指土体持水状态（可用体积含水率、质量含水率或饱和度表达）与土中水的能量状态（可用吸力、势能或相对湿度表达）之间的关系^[1],描述了土中水（液相）与土（固相）间的能量平衡状态^[2]。

土体持水的物理化学机制有毛细、吸附与渗透。毛细机制与弯液面上的表面张力有关,受Young- Laplace方程控制（土体孔径越小其毛细吸力越大）,当土中水能量状态低至其空化点以下,毛细水汽化,相应毛细机制消失。吸附机制指（黏）土颗粒表面由于van der Waals力、双电层力、阳离子水化等造成的短程吸附作用导致的土中水能量降低,可用Kelvin方程描述。通常将毛细与吸附机制导致的土中水相对自由水状态而言势能的降低称为基质势（用负值表达）；将对应土体固相的持水/吸水能力（用正值表达）称为基质吸力；二者数值相等,符号相反。此外,渗透机制指溶解盐会降低土中水的能量状态,其能量降低程度与溶质类型和浓度相关。

SWRC表征土体持水能力,原本用于农业领域,随后的研究表明非饱和状态下的渗流、强度与变形等土力学行为^[1]乃至土的成分与结构等土质学特征^[2]均与SWRC密切相关,因而成为描述非饱和状态下土体性状的关键指标与工具,在非饱和土力学研究与工程应用中具有基础的、特殊的、关键的作用与位置,其研究成果亦非常丰富^[1-3],但大部分SWRC^[1-3]由轴平移方法获得,受陶土板进气值限制,其基质吸力目前仅能控制在1500 kPa以内,根据Kelvin方程（式（1））可知其对应的相对湿度RH为98.9%以上（环境温度25℃时）。然而,自然界中RH变化范围远比轴平移方法所能控制的RH范围大；例如中国主要城市各月RH统计^[4]显示,最低14%（拉萨,2月）,最高90%（长沙,10月与11月；海口,1月）；因而高吸力（RH低于95%或吸力高于7.1 MPa）下的土体持水特征研究具有特别的现实意义。

高吸力下土体持水特征测试常用方法主要有：滤纸法^[5-8]（适用于10～100 MPa）、冷镜湿度计法^[9]（适用于0～300 MPa）、饱和盐溶液等压湿度控制法^[10-12]（适用于3～400 MPa）、干湿空气混合湿度控制法^[13]（适用于7～600 MPa）等。前两种为吸力量测方法,难以控制吸力；后两种为吸力控制方法,适于测定完整的SWRC,但蒸汽平衡历时漫长,如文献[14]采用饱和盐溶液等压湿度控制法测定3.3～359.1 MPa吸力范围的膨胀土SWRC耗时319 d,阻碍了高吸力下土体持水特征的研究与应用。

近年来出现的动态露点等温线方法^[15-16]使得上述情况有所改观,该方法能够在较短时间（几天内）获得相对湿度在3%～95%范围的SWRC。为了给高吸力下相关土力学工程实践提供依据和参考,本文采用该方法：①对从黏土质砂到极高塑性黏土的9种土样开展高吸力下吸脱湿过程持水特征测试；②采用RH = 95%下的吸湿质量含水率w₉₅量化高吸力下持水能力,采用吸脱湿SWRC间滞回区域面积HHA量化高吸力下SWRC滞回效应强弱,采用同一吸力下吸脱湿SWRC含水率之差最大值Δw_max量化水力滞回导致的含水率偏差范围；③探讨w₉₅、HHA、Δw_max与液限w_L、塑性指数I_p、阳离子交换量CEC、比表面积SSA等土性参数间量化关系；④构建修正的Fredlund-Xing模型描述高吸力下吸脱湿过程SWRC。

1. 研究方案

1.1 试验土样

试验采用Wyoming膨润土、宁明膨胀土、荆门黄褐色膨胀土、Denver黏土岩、荆门棕褐色膨胀土、武汉黏土、三门峡粉质黏土、郑州粉土和开封黄河砂共9种土样,其物理性质指标测试结果见表1。

表 1 试验用土的物理性质指标

Table 1. Physical property indexes of soils

土样名称	颗粒相对质量密度^[17]	颗粒所占比例^[17]/%				液限^[17]/%	塑限^[17]/%	塑性指数	塑性^[18]	USCS 定名^[19]	比表面积EGME^[20]/(m²·g^-1)	CEC(NH₄⁺)^[17]/( mmol·kg^-1)	自由膨胀率^[21]/%	膨胀潜势^[21]
土样名称	颗粒相对质量密度^[17]	>0.075 mm	0.005~0.075 mm	0.002~0.005 mm	<0.002 mm	液限^[17]/%	塑限^[17]/%	塑性指数	塑性^[18]	USCS 定名^[19]	比表面积EGME^[20]/(m²·g^-1)	CEC(NH₄⁺)^[17]/( mmol·kg^-1)	自由膨胀率^[21]/%	膨胀潜势^[21]
Wyoming膨润土	2.70	5.1	14.2	3.8	77.0	217.7	33.1	184.5	极高	CH	500.1	621	775	强
宁明膨胀土	2.73	5.8	2.9	0.3	91.1	164.6	37.6	127.0	极高	CH	393.0	508	655	强
荆门黄褐色膨胀土	2.75	2.1	47.4	21.5	29.0	62.9	25.5	37.4	高	CH	236.5	309	75	中
Denver黏土岩	2.72	8.5	44.2	20.6	26.7	46.0	23.2	22.8	中等	CL	88.7	160	38	—
荆门棕褐色膨胀土	2.72	11.6	42.4	26.7	19.3	41.8	20.7	21.1	中等	CL	110.4	254	40	弱
武汉黏土	2.73	5.5	63.9	9.6	21.0	40.3	18.2	22.1	中等	CL	112.5	154	39	—
三门峡粉质黏土	2.72	0.5	82.8	5.7	11.1	35.2	19.1	16.1	中等	CL	108.5	149	30	—
郑州粉土	2.70	8.7	77.4	2.5	11.4	26.1	14.8	11.3	低	CL	27.7	50	9	—
开封黄河砂	2.68	62.6	36.0	0.2	1.3	—	—	—	—	SC	7.5	16	—	—

下载: 导出CSV

| 显示表格

1.2 蒸汽吸附分析仪工作原理

美国Meter公司生产的蒸汽吸附分析仪（Vapor Sorption Analyzer, VSA）见图1,用于测试动态的或平衡的水分吸附等温线（moisture sorption isotherms）,即给定温度下的SWRC。

图 1 蒸汽吸附分析仪

Figure 1. Vapor sorption analyzer

下载: 全尺寸图片幻灯片

VSA采用动态露点等温线（dynamic dew point isotherm, DDI）方法获得动态水分吸附等温线；采用控制湿度下测试样质量的方法（dynamic vapor sorption, DVS）获得静态（平衡）水分吸附等温线。其工作原理是：DDI方法既不控制含水率亦不控制RH,而是在用湿气流湿化/干气流脱湿过程中通过高精度天平测试样质量变化、通过冷镜露点法测RH得到动态水分吸附等温线。DVS方法是在控制RH条件下跟踪试样质量变化以获得静态（平衡）水分吸附等温线,每级湿度下的稳定（平衡）标准可用预设时间间隔或预设试样质量变化幅度。Likos等^[15]与Arthur等^[22]试验证实了DDI方法的精度与可靠性,系统试验研究表明DDI方法与DVS方法获得的SWRC无明显偏差。

VSA相对湿度控制范围是3%～95%（25℃下对应吸力范围为7.1～482.9 MPa）,试样质量建议范围为500～5000 mg,试样温度可控范围为15℃～60 ℃。

1.3 高吸力下持水特征试验方案

采用VSA测定9种土样高吸力下SWRC。将土样风干碾散过0.5 mm筛,置于烘箱中烘至衡量（105 ℃下烘干48 h）,取出置于干燥器内冷却至室温,每次取约1000 mg试验。试验方案与步骤如下：

（1）仪器准备

将VSA在工作台上调平,安装干燥管,水箱中注入去离子水,检查天平,标定相对湿度。

（2）设定SWRC量测步骤

鉴于采用DDI方法与DVS方法测土样SWRC具有很好的一致性^[22],且前者测试历时约为后者的1/25^[15],因此采用DDI方法。由于试样是完全烘干后的,质量含水率w为0,第1步设置初始RH为3%,最终RH为95%,用DDI方法采集吸湿过程w数据RH,即土样SWRC吸湿边界线；第2步设置初始RH为95%、最终RH为3%,目的是获得土样SWRC的脱湿曲线；由于脱湿曲线起始点为RH=95%对应点,所以并非脱湿边界线。此外,设定样品室温度为25℃。

（3）测定SWRC

将试样置入样品杯,放入试样室开始测试,测试由VSA自动执行。试验结束后,取出试样杯称试样质量,烘干后复测干土质量,并在计算机终端下载SWRC试验数据。

2. 试验结果分析

由于Kelvin方程（式(1)）确定了吸力ψ与RH间关系,VSA试验结果既可用w-RH关系亦可用w-ψ关系表达；鉴于非饱和土力学中常用w-ψ关系,因此将9种土样w-ψ关系表达的SWRC绘制在图2中,将每种土样试验历时,试验数据点数,吸湿起始点、吸湿结束点（脱湿起始点）、脱湿结束点对应的RH,ψ,w值列在表2中。

图 2 高吸力下的SWRC（RH = 3%~95%）

Figure 2. Soil-water retention curves under high suctions

下载: 全尺寸图片幻灯片

表 2 高吸力下SWRC的吸湿/脱湿起始点与结束点

Table 2. Start point/end point of adsorption/desorption of SWRCs under high suctions

土样名称	试验历时/min	试验数据点数	吸湿起始点			吸湿结束点（脱湿起始点）			脱湿结束点
土样名称	试验历时/min	试验数据点数	RH	ψ/kPa	w/%	RH	ψ/kPa	w₉₅/%	RH	ψ/kPa	w/%
Wyoming膨润土	3135	203	0.0375	451647	0.67	0.9457	7686	21.46	0.0350	461106	0.87
宁明膨胀土	3777	205	0.0324	471755	0.76	0.9472	7468	19.45	0.0343	463909	1.02
荆门黄褐色膨胀土	2982	204	0.0297	484119	1.28	0.9456	7700	13.40	0.0331	468823	1.96
Denver黏土岩	2113	201	0.0313	476506	0.79	0.9476	7410	6.79	0.0338	465943	1.14
荆门棕褐色膨胀土	2523	212	0.0299	483240	0.89	0.9474	7439	9.63	0.0245	510367	1.28
武汉黏土	2343	202	0.0294	485114	0.44	0.9455	7715	7.37	0.0326	470929	0.91
三门峡粉质黏土	2186	203	0.0244	510892	0.43	0.9478	7380	6.19	0.0327	470498	0.97
郑州粉土	1787	197	0.0246	509791	0.42	0.9500	7061	2.34	0.0245	510347	0.52
开封黄河砂	1295	148	0.0245	510383	0.65	0.9519	6786	1.28	0.0247	509237	0.67

下载: 导出CSV

| 显示表格

$ψ = - \frac{R T}{ν_{w}} \ln RH$ ,

(1)

式中 ψ为吸力（kPa）；R = 8.314 J/(K·mol),为通用气体常数；T为热力学温度,样品室温度25℃对应T = 298.15 K； $ν_{w}$ = 1.8×10^-5 m³/mol,为水的摩尔体积；RH为相对湿度,量纲为1。

由表2可见,9种土样中宁明膨胀土测试历时最长,为3777 min,即2.6 d,且有205个数据点。一方面说明试验历时为文献[14]中饱和盐溶液等压湿度控制法测试历时（319 d）的1/122,极大缩短了测试时间；另一方面说明数据点丰富,能够更好描述土样的持水特征。

从以下3方面探讨高吸力下的SWRC：①反映土样亲水性强弱的持水能力特征点w₉₅；②描述水力滞回效应强弱的SWRC滞回区域面积HHA；③反映水力滞回导致的含水率偏差范围的SWRC滞回效应最大值Δw_max。

2.1 高吸力下持水能力特征点w₉₅

由图2与表2可见：9种土样从完全干燥状态吸湿结束时（控制RH为95%）对应的质量含水率w₉₅差别相当大,高至Wyoming膨润土的21.46%,宁明膨胀土的19.45%,表明持水能力/亲水性非常强,低至开封黄河砂的1.28%,表明其持水能力/亲水性非常弱。由此可见,w₉₅作为RH = 95%下的吸湿质量含水率,是表征土样持水（吸水）能力的特征点,其数值反映了土样亲水性强弱。

将9种土样的w₉₅与表1给出的CEC对应绘制在图3中,可见w₉₅与CEC具有非常高的线性相关性,其线性拟合结果见式（2）,决定系数R²高达0.99。

$w_{95} = 0.035 CEC + 1.3$ 。

(2)

图 3 9种土的w₉₅-CEC关系

Figure 3. Relationship between w₉₅ and CEC of 9 kinds of soils

下载: 全尺寸图片幻灯片

鉴于《膨胀土地区建筑技术规范》^[21]采用CEC分别为340,260,170 mmol/kg作为判断土样膨胀潜势强（CEC > 340 mmol/kg）、中（CEC =260～340 mmol/kg）、弱（CEC = 170～260 mmol/kg）的界限,将其代入式（2）得到采用w₉₅判别土样膨胀潜势的相应界限分别为13.0%,10.2%,7.1%。即：w₉₅ > 13.0%,判别为强膨胀土；w₉₅=10.2%～13.0%,判别为中膨胀土；w₉₅=7.1%～10.2%,判别为弱膨胀土；w₉₅<7.1%,判别为非膨胀土。

下面对上述方法的可行性进行验证：表3为根据《膨胀土地区建筑技术规范》^[21]测定的Denver膨润土、信阳黏土、商用高岭土的CEC、自由膨胀率及确定的膨胀潜势。采用VSA获得的w₉₅亦列在表3中,Denver膨润土w₉₅为22.47%,大于13.0%,根据w₉₅判别为强膨胀土；信阳黏土w₉₅为8.77%,在7.1%～10.2%,根据w₉₅判别为弱膨胀土；商用高岭土w₉₅为1.03%,小于7.1 %,根据w₉₅判别为非膨胀土。表明采用w₉₅作为膨胀潜势判别标准与采用CEC的判别结果具有很好的一致性,至于如何提高其判别精度,则应收集更多的试验结果,确定更为可靠的强、中、弱膨胀潜势的w₉₅界限值。

表 3 3种土样膨胀潜势的判别指标

Table 3. Eestimation indexes of swelling potential of 3 kinds of soils

土样名称	CEC(NH₄⁺)^[17]/(mmol·kg^-1)	自由膨胀率^[21]/%	膨胀潜势^[21]	w₉₅/%
Denver膨润土	976	111	强	22.47
信阳黏土	206	40	弱	8.77
商用高岭土	4	—	—	1.03

下载: 导出CSV

| 显示表格

2.2 高吸力下SWRC滞回区域面积HHA

由图2可见,高吸力段(RH = 3%～95%),9种土样SWRC均呈现或强或弱的水力滞回行为,即相同吸力下脱湿路径上的w高于吸湿路径上的w。

SWRC水力滞回效应强弱可用脱湿与吸湿曲线间滞回区域面积HHA (hydraulic hysteresis area)描述。根据图2中的脱湿与吸湿曲线数据在Origin软件中数值积分获得9种土样滞回区域面积见表4。

表 4 高吸力下SWRC滞回效应特征参数（R_H = 3%~95%）

Table 4. Hysteresis properties of SWRCs under high suctions

土样名称	滞回区域面积HHA/kPa	HHA 相对值	Δw_max/%	ψ(Δw_max)/kPa	RH(Δw_max)	w_d(Δw_max)/%	w_w(Δw_max)/%	w相对误差/%
土样名称	滞回区域面积HHA/kPa	HHA 相对值	Δw_max/%	ψ(Δw_max)/kPa	RH(Δw_max)	w_d(Δw_max)/%	w_w(Δw_max)/%	$\frac{Δ w_{\max}}{w_{d}}$	$\frac{Δ w_{\max}}{w_{w}}$
Wyoming膨润土	13641	1.00	6.27	68737	0.6069	13.08	6.81	47.9	92.1
宁明膨胀土	10717	0.79	4.18	72744	0.5895	11.00	6.81	38.0	61.4
荆门黄褐色膨胀土	4512	0.33	1.77	37559	0.7612	9.87	8.10	17.9	21.9
Denver黏土岩	2367	0.17	0.69	37184	0.7633	4.80	4.11	14.4	16.8
荆门棕褐色膨胀土	3520	0.26	1.38	30993	0.7984	7.41	6.03	18.6	22.9
武汉黏土	3412	0.25	1.15	31077	0.7979	5.13	3.98	22.4	28.9
三门峡粉质黏土	3439	0.25	1.04	89267	0.5228	3.55	2.51	29.3	41.4
郑州粉土	1002	0.07	0.37	155335	0.3235	1.22	0.85	30.3	43.5
开封黄河砂	196	0.01	0.16	221700	0.1999	0.88	0.72	18.2	22.2

下载: 导出CSV

| 显示表格

由图2与表4可见：9种土样HHA差别显著,对于强膨胀土,Wyoming膨润土高达13641 kPa,宁明膨胀土亦高达10717 kPa；对于中膨胀土,荆门黄褐色膨胀土为4512 kPa；对于弱膨胀土,荆门棕褐色膨胀土为3520 kPa；武汉黏土、三门峡粉质黏土亦达到3412,3439 kPa；而郑州粉土低至1002 kPa；开封黄河砂仅有196 kPa,是Wyoming膨润土的1/70,滞回效应可忽略不计。

以上数据表明：HHA与土样膨胀潜势正相关,将9种土样的HHA与表1给出的CEC、SSA对应绘制在双y轴坐标系中,见图4,可见HHA与CEC、SSA均具有非常高的线性相关性,其过坐标原点的线性拟合结果为

图 4 9种土的CEC-SSA-HHA关系

Figure 4. CEC-SSA-HHA relationship of 9 kinds of soils

下载: 全尺寸图片幻灯片

$CEC = 0 .049 HHA$ ,

(3)

$SSA = 0 .037 HHA$ ,

(4)

式中,HHA与CEC线性相关的R²达到0.94,HHA与SSA线性相关的R²高达0.97。

CEC与SSA是强烈相关但有明显区别的土性指标。CEC指单位质量土样吸附的可交换阳离子总量,反映由于颗粒带负电荷导致的土样吸附阳离子的能力,更大的CEC意味着更强的膨胀潜势。SSA指单位质量土样具有的表面积,反映土样的颗粒级配与矿物成分,更大的SSA意味着更大的塑性,更强的膨胀潜势,更强的亲水性。HHA与CEC、SSA均正相关表明,滞回效应大小既受土样吸附阳离子能力的强烈影响,亦受土样比表面积控制,且从R²上看,其受土样比表面积影响更大。

将8种细粒土HHA与其液限w_L、塑性指数I_p绘制在双y轴坐标系中,见图5。

图 5 8种细粒土的w_L-I_p-HHA关系

Figure 5. w_L-I_p-HHA relationship of 8 kinds of fine-grained soils

下载: 全尺寸图片幻灯片

可见HHA与w_L、I_p均正相关,过坐标原点的线性拟合结果分别为

$w_{L} = 0 .015 HHA$ ,

(5)

$I_{p} = 0.012 HHA$ 。

(6)

式中,HHA与w_L线性相关的R²高达0.97,HHA与I_p线性相关的R²亦达到0.92,表明细粒土HHA与w_L具有更高的线性相关性。由于w_L、I_p代表土样塑性大小,表明HHA受土样塑性大小控制。

2.3 高吸力下SWRC滞回效应最大值Δw_max

SWRC滞回效应指水力滞回（脱/吸湿路径）对SWRC（此处是w-ψ关系）的影响,此处SWRC滞回效应最大值定义为相同ψ下,脱湿曲线对应含水率w_d与吸湿曲线对应含水率w_w之差的最大值Δw_max。根据该定义,在图2中搜索到每种土的Δw_max,及其对应的吸力ψ、RH、脱湿曲线含水率w_d、吸湿曲线含水率w_w,列入表4中。

由表4与图2可见,9种土样Δw_max存在较大差异：强膨胀土Δw_max相当大,如Wyoming膨润土为6.27%,宁明膨胀土为4.18%；中膨胀土,荆门黄褐色膨胀土为1.77 %；弱膨胀土,荆门棕褐色膨胀土为1.38%；中等塑性黏土,武汉黏土为1.15%,三门峡粉质黏土为1.04%,Denver黏土岩为0.69%；低塑性黏土,郑州粉土为0.37%；开封黄河砂仅为0.16%。对比表1可见,土样Δw_max值与w_L,I_p,CEC,SSA均正相关,反映了水力滞回导致的含水率偏差范围。取Δw_max与相应w_d、w_w的比值作为SWRC的含水率相对误差,计算9种土样Δw_max/w_d、Δw_max/w_w结果列于表4与图6。

图 6 由于水力滞回导致的含水率相对误差

Figure 6. Relative errors of water content due to hydraulic hysteresis

下载: 全尺寸图片幻灯片

表4与图6中,同一土样w_w小于w_d,因此采用吸湿曲线含水率w_w作为分母的相对误差比采用脱湿曲线含水率w_d作为分母的相对误差要大。由相对误差计算结果可见：在高吸力范围,不考虑水力滞回会导致持水特征曲线SWRC上质量含水率w最高92.1 %（Wyoming膨润土）的相对误差。

另一方面,如果工程中需要用当前含水率估算吸力,则同一w下,脱湿路径上吸力与吸湿路径上吸力的比值最高可达2.71（图2中Wyoming膨润土,w为6%,脱湿路径上吸力205959 kPa,吸湿路径上吸力76119 kPa,比值为2.71）,印证了高吸力下强膨胀土SWRC滞回效应之强烈。

2.4 高吸力下SWRC滞回行为机制探讨

由图2及2.2节、2.3节对SWRC滞回区域面积HHA与滞回效应最大值Δw_max的分析,可知高吸力下,9种土样中,高塑性（液限高、塑性指数大）、阳离子交换量大、比表面积大的土样滞回效应大,在此对高吸力下SWRC滞回行为机制展开探讨。

水力滞回行为机制解释主要有墨水瓶效应^[1]与接触角效应^[13],前者由于土中孔隙尺寸与形状（可用孔径r概化）的不均匀性所致,后者由于吸湿接触角 $α_{w}$ 通常大于脱湿接触角 $α_{d}$ 所致^[23]。墨水瓶效应与接触角效应均属毛细范畴,受Young-Laplace方程控制：

$ψ_{c} = u_{a} - u_{w} = \frac{2 T_{s} \cos α}{r}$ 。

(7)

式中 ψ_c为毛细吸力,可用孔隙气压力u_a与孔隙水压力u_w之差表达；T_s为水的表面张力； $α$ 为接触角；r为毛细半径。当吸力增大、土中水能量状态低到其空化点以下,毛细水汽化,（土中毛细水比自由水空化点更低,但其对应的ψ_c或RH尚无定论,一般而言RH低于80%,毛细水难以存在^[2]）,墨水瓶效应与接触角效应难以解释高吸力下的滞回行为。

Lu等^[24]给出了高吸力下SWRC滞回机制的一种解释：非膨胀黏土,脱湿与吸湿主要发生在颗粒表面,因而滞回效应较小。膨胀性土,脱湿与吸湿亦可发生于存在可交换阳离子的晶层之间,吸湿时阳离子水化吸附极性水分子,水分子需克服较大的,van der Waals力方能进入层间；脱湿时由于层间距离较大van der Waals力较弱,因此相同ψ下,吸湿路径对应较小的w,脱湿路径对应较大的w。在此结合文中数据探讨该解释的合理性：

高吸力下Wyoming膨润土HHA=13641 kPa,Δw_max=6.27 %,滞回效应最大；自由膨胀率为775 %,膨胀性最大；SSA高达500.1 m²/g,表明其包含了较大的层间比表面积^[25],CEC高达621 mmol/kg,表明其可交换阳离子总量大；印证了吸湿状况下膨胀土层间阳离子水化受van der Waals力阻碍致其含水率较低的机制解释。

相对而言,郑州粉土,SSA为27.7 m²/g,CEC为50 mmol/kg,比表面积小,可交换阳离子总量小。一方面其自由膨胀率（9%）为Wyoming膨润土的1.2 %,表明其为非膨胀土；另一方面,其HHA（1002 kPa）为Wyoming膨润土的7%,Δw_max（0.37%）为Wyoming膨润土的6%,滞回效应很小。印证了非膨胀土脱湿吸湿主要发生在颗粒表面（而非层间）致其滞回效应较小的机制解释。

以上讨论明晰了2.2节中的论述：高吸力下滞回效应既受CEC的强烈影响,亦受SSA控制,且从R²上看,其受SSA影响更大。至于是否考虑滞回效应,应视土样与工况等具体情况而定,鉴于Wyoming膨润土具有最大的HHA,可采用其作为标准,其它土样的HHA相对值亦列入表4中,可见从开封黄河砂的0.01到宁明膨胀土的0.79,HHA相对值可作为实际工况中是否考虑滞回效应的依据。

3. 高吸力下持水特征的本构描述

基于对以上9种土样高吸力下SWRC规律性的认知,经尝试,建议了Fredlund-Xing（1994）模型^[1]的修正形式描述高吸力下SWRC：

。

$\begin{matrix} w (ψ) = C (ψ) \frac{w_{u}}{{\ln [\exp (1) + {(\frac{(ψ - ψ_{b})}{a})}^{n}]}^{m}} \\ = [1 - \frac{\ln (1 + ψ / ψ_{r})}{\ln (1 + 10^{6} / ψ_{r})}] \frac{w_{u}}{{\ln [\exp (1) + {(\frac{(ψ - ψ_{b})}{a})}^{n}]}^{m}} 。 \end{matrix}$

(8)

式中 w为质量含水率（%）；ψ为吸力（kPa）；w(ψ)为质量含水率表达的SWRC；C(ψ)为修正系数函数,能够迫使SWRC通过零含水率对应吸力10⁶ kPa的点。4个参数a,n,m,ψ_r为Fredlund-Xing(1994)模型已有参数,a为与拐点位置相关的土性参数（kPa）；n为与拐点处斜率相关的土性参数,量纲为1,m为与曲率相关的土性参数,量纲为1,ψ_r为残余吸力（kPa）。

为描述高吸力下SWRC,将Fredlund-Xing模型中饱和质量含水率w₀用w_u替代,w_u为高吸力下SWRC质量含水率上边界（最大值）,相应吸力为ψ_b。增加土性参数ψ_b是式（8）与Fredlund-Xing模型形式上的主要区别,ψ_b单位为kPa,物理意义是高吸力下SWRC吸力下边界。由此可见,修正的Fredlund-Xing模型共有6个参数：w_u,ψ_b,ψ_r,a,n,m。为区分吸湿与脱湿过程,相应参数增加上标w（wetting,表示吸湿）或d（drying,表示脱湿）。

采用式（8）分别对9种土样高吸力下SWRC进行拟合,流程如下：

（1）采用表2给出的吸湿结束点（脱湿起始点）含水率w₉₅值作为w_u值,对应吸力值作为ψ_b值；w_u与ψ_b为吸脱湿过程共用参数值。

（2）对于吸湿过程,采用表2给出的吸湿起始点吸力值作为ψ_r值；至此,吸湿过程所需确定的参数仅余a,n,m。依据式（8）在Excel中通过规划求解功能对吸湿过程SWRC数据进行非线性拟合获得a,n,m最优解以及拟合结果的决定系数R²。

（3）对于脱湿过程,采用表2给出的脱湿结束点吸力值作为ψ_r值；依据式（8）在Excel中通过规划求解功能对脱湿过程SWRC数据进行非线性拟合获得a,n,m最优解以及R²。

前8种土样采用上述流程获得的模型参数与R²见表5,可见8种土样,无论脱湿吸湿,修正的Fredlund-Xing模型再现结果的R²均在0.99以上。为节约篇幅,其修正Fredlund-Xing模型数值再现结果亦放在图2中,可见式（8）相当好地再现了高吸力下8种土样脱湿与吸湿过程SWRC。

表 5 高吸力下修正Fredlund-Xing模型参数

Table 5. Modified Fredlund-Xing model parameters for 9 kinds of soils under high suctions

土样名称	共用参数		吸湿过程					脱湿过程
土样名称	w_u/%	ψ_b/kPa	ψ_r^w/kPa	a^w/kPa	n^w	m^w	R²	ψ_r^d/kPa	a^w/kPa	n^d	m^d	R²
Wyoming膨润土	21.46	7686	451647	248967	0.94	10.54	1.00	461106	90593	1.06	2.09	0.99
宁明膨胀土	19.45	7468	471755	245967	0.82	8.02	1.00	463909	90230	0.94	2.13	1.00
荆门黄褐色膨胀土	13.40	7700	484119	237692	0.61	4.34	1.00	468823	91452	0.82	1.90	1.00
Denver黏土岩	6.79	7410	476506	267853	0.57	4.14	1.00	465943	91075	0.66	1.79	1.00
荆门棕褐色膨胀土	9.63	7439	483240	87263	0.64	2.81	1.00	510367	86362	0.80	1.96	0.99
武汉黏土	7.37	7715	485114	81572	0.67	3.77	1.00	470929	80590	0.71	2.39	0.99
三门峡粉质黏土	6.19	7380	510892	213682	0.67	4.50	1.00	470498	91368	0.82	1.56	1.00
郑州粉土	2.34	7061	509791	30152	0.78	1.17	0.99	510347	31576	0.58	0.86	1.00
开封黄河砂	1.28	6786	—	76542	0.66	1.04	0.96	—	76542	0.66	1.04	0.96

下载: 导出CSV

| 显示表格

对第9种土样开封黄河砂,由图2可见其SWRC滞回行为非常小,因此采用同一组模型参数描述其SWRC,此外不考虑式（8）中修正系数函数C(ψ),即C(ψ) = 1,对吸脱湿SWRC数据进行非线性拟合获得a,n,m最优解见表5,同时描述其吸脱湿SWRC的R²达到0.96,由图2的再现结果可见式（8）非常好地对其SWRC进行了本构描述。

此外,对比表5中吸脱湿过程a,n,m值可见,对同种土样,相对脱湿过程,吸湿过程a值较大说明其拐点位置靠后,n值较低说明其拐点处斜率相对较小,m值较大说明其曲率较大。

由此可见,文中建议的Fredlund-Xing模型修正形式有6个参数：w_u,ψ_b,ψ_r值直接从实测SWRC数据确定,a,n,m值通过非线性拟合获得。Fredlund- Xing(1994)模型中4个参数（a,n,m,ψ_r）需通过非线性拟合获得,还有一个参数饱和质量含水率w₀需试验测定（从实测SWRC数据确定）,实际上有5个参数。文中建议的修正形式仅增加了一个参数ψ_b,且ψ_b与w_u、ψ_r直接从实测SWRC数据确定,只需非线性拟合3个参数（a,n,m）,且该修正形式能够很好地分别描述吸湿与脱湿过程中SWRC强烈的非线性行为。

4. 结论

（1）9种土样高吸力下吸湿与脱湿过程持水特征测试最长历时为2.6 d,有205个数据点,说明动态露点法极大缩短了高吸力下SWRC测试时间,且数据点丰富,能够更好描述土样持水特征。

（2）w₉₅作为RH = 95%下的吸湿质量含水率,是表征土样持水能力的特征点,反映土样亲水性强弱,w₉₅与CEC具有非常高的线性相关性,w₉₅作为膨胀潜势的判别标准是可行的。

（3）SWRC水力滞回效应强弱可用吸脱湿曲线间滞回区域面积HHA描述,HHA与CEC、SSA均有非常高的线性相关性,表明滞回效应强弱既受土样吸附阳离子能力的强烈影响,亦受土样SSA控制,且从R²看,其受SSA影响更大。HHA与w_L,I_p均正相关,表明HHA受土样塑性大小控制。

（4）SWRC滞回效应最大值可用同一ψ下脱湿与吸湿曲线对应w之差最大值Δw_max量化,土样Δw_max值与w_L,I_p,CEC,SSA均正相关,反映了水力滞回导致的w偏差范围。高吸力范围,不考虑水力滞回会导致SWRC上w最高92.1 %的相对误差。如工程中需要用当前w估算ψ,则同一w下,脱湿路径上ψ与吸湿路径上ψ的比值最高可达2.71。

（5）建议采用Fredlund-Xing模型的修正形式描述高吸力下SWRC,方程有6个参数,w_u,ψ_b,ψ_r值直接从实测SWRC数据确定,a,n,m值通过非线性拟合获得。对9种土样SWRC的修正Fredlund -Xing模型数值再现结果表明该方程可分别很好地描述吸湿与脱湿过程中SWRC强烈的非线性行为。

图 1 巴西圆盘试样示意图

Figure 1. Illustration for carbonaceous slate samples with bedding angles of 0°, 30°, 45°, 60° and 90°

下载: 全尺寸图片幻灯片

图 2 加载仪器及试件放置

Figure 2. Experimental setup and placement of samples

下载: 全尺寸图片幻灯片

图 3 炭质板岩不同层理倾角下巴西劈裂典型位移-荷载曲线

Figure 3. Load-displacement curves of Brazilian carbonaceous slate tests under different bedding angles

下载: 全尺寸图片幻灯片

图 4 炭质板岩即时烘干试样巴西劈裂破坏形态图

Figure 4. Brazilian splitting failure modes of drying carbonaceous slate samples

下载: 全尺寸图片幻灯片

图 5 炭质板岩静置风化试样巴西劈裂典型破坏模式

Figure 5. Typical Brazilian failure modes of weathering carbonaceous slate samples

下载: 全尺寸图片幻灯片

图 6 炭质板岩抗拉强度对比曲线图

Figure 6. Comparative curves of tensile strength of two kinds of carbonaceous slates

下载: 全尺寸图片幻灯片

图 7 即时烘干炭质板岩SEM照片

Figure 7. SEMs of drying carbonaceous slate samples

下载: 全尺寸图片幻灯片

图 8 静置风化炭质板岩SEM照片

Figure 8. SEMs of weathering carbonaceous slate samples

下载: 全尺寸图片幻灯片

图 9 炭质板岩巴西劈裂受力及劣化示意图

Figure 9. Illustration for degradation of a Brazilian carbonaceous slate sample

下载: 全尺寸图片幻灯片

表 1 炭质板岩抗拉强度统计表

Table 1 Tensile strengths of two kinds of carbonaceous slate samples

层理倾角/(°)	即时					静置
层理倾角/(°)	试样抗拉强度			均值	标准差	试样抗拉强度					均值	标准差
0	12.19	10.81	11.60	11.53	0.21	1.33	1.72	1.43	1.58	1.19	1.45	0.21
30	9.70	10.80	8.43	9.64	1.19	1.12	0.73	1.35	1.16	—	1.09	0.26
45	7.45	7.56	8.17	7.73	0.47	0.97	0.73	0.56	1.30	—	0.89	0.32
60	5.13	7.30	7.18	6.54	1.22	0.76	0.99	0.97	0.41	—	0.78	0.27
90	4.88	6.57	7.61	6.35	1.38	0.59	0.71	0.68	0.78	—	0.69	0.08
注：①抗拉强度单位均为MPa；②考虑静置试样易受人为影响产生数据离散,故总计完成21例试样。

下载: 导出CSV

参考文献(21)

[1]	CHO J W, KIM H, JEON S, et al. Deformation and strength anisotropy of Asan gneiss, Boryeong shale, and Yeoncheon schist[J]. International Journal of Rock Mechanics and Mining Sciences, 2012, 50(2): 158-169.
[2]	DAN D Q, KONIETZKY H, HERBST M. Brazilian tensile strength tests on some anisotropic rocks[J]. International Journal of Rock Mechanics and Mining Sciences, 2013, 58(2): 1-7.
[3]	吕坤, 高旭, 吴峥, 等. 弱黏结复合层状顶板冒顶隐患分级方法与应用[J]. 矿业科学学报, 2018, 3(3): 253-259. https://www.cnki.com.cn/Article/CJFDTOTAL-KYKX201803006.htm LÜ Kun, GAO Xu, WU Zheng, et al. Weakly bonding composite layered roof risk classification method and its application[J]. Journal of Mining Science and Technology, 2018, 3(3): 253-259. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KYKX201803006.htm
[4]	黄琪嵩, 程久龙. 层状底板采动应力场的解析计算模型研究[J]. 矿业科学学报, 2017, 2(6): 559-565. https://www.cnki.com.cn/Article/CJFDTOTAL-KYKX201706006.htm HUANG Qi-song, CHENG Jiu-long.Research on analytic calculation model for mining induced stress in multi-layered floor rock[J]. Journal of Mining Science and Technology, 2017, 2(6): 559-565. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KYKX201706006.htm
[5]	LI Z G, XU G L, HUANG P, et al. Experimental study on anisotropic properties of Silurian silty slates[J]. Geotechnical and Geological Engineering, 2017, 35(2): 1-12.
[6]	朱思尘, 李江腾. 干燥和饱水状态下含层理构造板岩巴西劈裂实验能量研究[J]. 中南大学学报(自然科学版), 2018, 49(8): 2024-2030. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201808024.htm ZHU Si-chen, LI Jiang-teng. Experimental energy study of Brazil test of bedded structural slate under dry and saturated conditions[J]. Journal of Central South University (Science and Technology), 2018, 49(8): 2024-2030. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD201808024.htm
[7]	张娜, 赵方方, 张毫毫, 等. 岩石气态水吸附特性及其影响因素实验研究[J]. 矿业科学学报, 2017, 2(4): 34-45. https://www.cnki.com.cn/Article/CJFDTOTAL-KYKX201704004.htm ZHANG Na, ZHAO Fang-fang, ZHANG Hao-hao, et al. Experimental study on water vapor adsorption of rock and its influencing factors[J]. Journal of Mining Science and Technology, 2017, 2(4): 336-347. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KYKX201704004.htm
[8]	张娜, 王水兵, 何枭, 等. 深部煤系页岩吸水及软化效应微观机理研究[J]. 矿业科学学报, 2019, 2(4): 308-317. https://www.cnki.com.cn/Article/CJFDTOTAL-KYKX201904004.htm ZHANG Na, WANG Shui-bing, HE Xiao, et al. Study on micromechanism of water absorption and its softening effect on shale rock in deep coal measures[J]. Journal of Mining Science and Technology, 2019, 2(4): 308-317. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KYKX201904004.htm
[9]	滕腾, 杜玉冰, 陈朋飞, 等. 砂岩变形率与水理效应的力学特性研究[J]. 矿业科学学报, 2020, 5(3): 342-352. https://www.cnki.com.cn/Article/CJFDTOTAL-KYKX202003013.htm TENG Teng, DU Yu-bing, CHEN Peng-fei, et al. Effects of deformation rate and hydrated condition on the mechanical property of sandstone[J]. Journal of Mining Science and Technology, 2020, 5(3): 342-352. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KYKX202003013.htm
[10]	郑明雄, 李保珠, 申岳龙. 贵州泥堡金矿凝灰岩的风化机理[J]. 昆明理工大学学报(自然科学版), 2019, 44(6): 31-38. ZHENG Ming-xiong, LI Bao-zhu, SHEN Yue-long. Weathering mechanism of tuff in Guizhou Nibao gold mine[J]. Journal of Kunming University of Science and Technology (Natural Science), 2019, 44(6): 31-38. (in Chinese)
[11]	WONG L N Y, MARUVANCHERY V, LIU G. Water effects on rock strength and stiffness degradation[J]. Acta Geotechnica, 2016, 11(4): 713-737.
[12]	TAVALLALI A, VERVOORT A. Effect of layer orientation on the failure of layered sandstone under Brazilian test conditions[J]. International Journal of Rock Mechanics and Mining Sciences, 2010, 47(2): 313-322.
[13]	NEZHAD M M, FISHER Q J, GIRONACCI E, et al. Experimental study and numerical modeling of fracture propagation in shale rocks during brazilian disk test[J]. Rock Mechanics and Rock Engineering, 2018, 51(6): 1755-1775.
[14]	KHANLARI G, RAFIEI B, ABDILOR Y. An experimental investigation of the Brazilian tensile strength and failure patterns of laminated sandstones[J]. Rock Mechanics and Rock Engineering, 2015, 48(2): 843-852.
[15]	李德建, 祁浩, 李春晓, 等. 含层理面煤试样的巴西圆盘劈裂实验及数值模拟研究[J]. 矿业科学学报, 2020, 5(2): 150-159. https://www.cnki.com.cn/Article/CJFDTOTAL-KYKX202002003.htm LI De-jian, QI Hao, LI Chun-xiao, et al. Brazilian disc splitting tests and numerical simulations on coal samples containing bedding planes[J]. Journal of Mining Science and Technology, 2020, 5(2): 150-159. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KYKX202002003.htm
[16]	班宇鑫, 傅翔, 谢强, 等. 页岩巴西劈裂裂缝形态评价及功率谱特征分析[J]. 岩土工程学报, 2019, 41(12): 2307-2315. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201912023.htm BAN Yu-xin, FU Xiang, XIE Qiang, et al. Evaluation of fracture morphology of shale in Brazilian tests andanalysis of power spectral characteristics[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(12): 2307-2315. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201912023.htm
[17]	刘运思, 傅鹤林, 伍毅敏, 等. 基于单弱面理论对板岩巴西劈裂试验研究[J]. 煤炭学报, 2013, 38(10): 1775-1780. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201310012.htm LIU Yun-si, FU He-lin, WU Yi-min, et al. Study on Brazilian splitting test for slate based on single weak plane theory[J]. Journal of China Coal Society, 2013, 38(10): 1775-1780. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201310012.htm
[18]	DEBECKER B, VERVOORT A. Experimental observation of fracture patterns in layered slate[J]. International Journal of Fracture, 2009, 159(1): 51-62.
[19]	丁长栋, 胡大伟, 周辉, 等. 考虑三维片理的板岩巴西劈裂试验研究[J]. 岩石力学与工程学报, 2019, 38(2): 90-101. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201902008.htm DING Chang-dong, HU Da-wei, ZHOU Hui, et al. Brazil test of slate considering three-dimensional schistosity effect[J]. Chinese Journal of Rock Mechanics and Engineering, 2019, 38(2): 90-101. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201902008.htm
[20]	GHOLAMI R, RASOULI V. Mechanical and elastic properties of transversely isotropic slate[J]. Rock Mechanics and Rock Engineering, 2014, 47(5): 1763-1773.
[21]	ISRM. Suggested methods for determining tensile strength of rock materials[J]. International Journal of Rock Mechanics and Mining Sciences and Geomechanics, 1978, 15(3): 99-103.

施引文献(23)

期刊类型引用(7)

1.	王东，王智佼，陈伟祥，陈党辉. 不同锚孔出水量下预应力树脂锚索锚固剂组合试验研究. 四川建筑. 2024(01): 180-182 . 百度学术
2.	彭亚雄，周子霈，姚颖康，刘运思，左清军. 爆破作用下层状围岩隧道突变失稳判据研究. 中国安全科学学报. 2024(01): 171-178 . 百度学术
3.	刘晓颖，郭青林，张博，尚东娟，余静，裴强强. 降雪、温差作用下大尺寸砂岩风化模拟试验研究. 岩石力学与工程学报. 2024(06): 1519-1534 . 百度学术
4.	刘金芳，屈玉龙. 风化作用对侵蚀构造地貌碳质砂岩高边坡的影响. 科技与创新. 2023(05): 119-121+124 . 百度学术
5.	尹大伟，丁屹松，汪锋，江宁，谭毅，李宗蓄. 压力水浸泡下煤岩抗拉特性试验研究. 岩石力学与工程学报. 2023(S1): 3178-3191 . 百度学术
6.	王世鸣，熊咸瑞，王嘉琪，颜世军，吴秋红，翁磊. TSL材料防护作用下砂岩动态抗拉性能的试验研究. 振动与冲击. 2023(21): 192-199+218 . 百度学术
7.	陈伟，万文，赵延林，王卫军，吴秋红，吴小凡，谢森林. 高湿环境平行双裂隙砂岩单轴压缩破坏及裂纹扩展特性. 岩土工程学报. 2021(11): 2094-2104 . 本站查看