初始主应力方向角对饱和珊瑚砂液化特性影响的试验

马维嘉; 陈国兴; 秦悠; 吴琪

doi:10.11779/CJGE202003022

初始主应力方向角对饱和珊瑚砂液化特性影响的试验 English Version

马维嘉^1,,
陈国兴^{1, 2, ,},
秦悠¹,
吴琪¹

1.
南京工业大学岩土工程研究所,江苏　南京　210009
2.
江苏省土木工程防震技术研究中心,江苏　南京　210009

基金项目:

国家自然科学基金项目 51678299

详细信息

作者简介:
马维嘉（1991— ）,男,博士研究生,主要从事土动力学的研究工作。E-mail：njtechmwj@163.com

通讯作者:
陈国兴, E-mail：gxc6307@163.com

中图分类号: TU443
计量
- 文章访问数: 322
- HTML全文浏览量: 33
- PDF下载量: 151
出版历程
- 收稿日期: 2019-03-24
- 网络出版日期: 2022-12-07
- 刊出日期: 2020-02-29

Experimental studies on effects of initial major stress direction angles on liquefaction characteristics of saturated coral sand

MA Wei-jia^1,,
CHEN Guo-xing^{1, 2, ,},
QIN You¹,
WU Qi¹

1.
Institute of Geotechnical Engineering, Nanjing Tech University, Nanjing 210009, China
2.
Civil Engineering and Earthquake Disaster Prevention Center of Jiangsu Province, Nanjing 210009, China

摘要

摘要: 为探究不同的初始主应力方向角α₀条件下南沙岛礁饱和珊瑚砂的液化特性,对饱和珊瑚砂进行了一系列不排水循环扭剪试验。试验结果表明：①α₀对饱和珊瑚砂液化特性的影响显著；α₀≠45°时的超静孔压u_e不会到达初始有效平均主应力 ${p^{'}}_{0}$ ；但α₀ = 45°时u_e会达到 ${p^{'}}_{0}$ 。②不同α₀条件下试样各应变分量发展形式的差异显著,但均会产生较大的广义剪应变幅值 $γ_{ga}$ ,且超静孔压比u_e/ ${p^{'}}_{0}$ 与 $γ_{ga}$ 之间的相关性可表示为正切函数。③以惯用的循环应力比CSR作为应力水平指标时,对于给定的CSR,到达液化标准所需的循环次数N_L随着α₀的增大而减小；采用本文定义的修正循环应力比CSR_M作为应力水平指标时,不同α₀条件下CSR_M与N_L的关系曲线规准化为一条曲线,且文献中3种饱和砂土的原始试验数据也验证了CSR_M的适用性。
- 饱和珊瑚砂 /
- 初始主应力方向角 /
- 超静孔隙水压力 /
- 变形特性 /
- 循环强度 /
- 修正循环应力比
Abstract: In order to investigate the liquefaction characteristics of saturated coral sand from Nansha Islands under different initial major stress direction angles (α₀), a series of undrained cyclic torsional shear tests for saturated coral sand are carried out by using the GDS hollow cylinder torsional apparatus. A remarkable finding of this study is that, α₀ has significant influence on the liquefaction characteristics of coral sand, the excess pore water pressure (u_e) for case α₀ ≠ 45° cannot reach the initial mean principal stress ${p^{'}}_{0}$ , but the u_e for case α₀ = 45° can eventually reach do. The developments of strain components for various α₀ cases are different, but the specimens for all α₀ cases will produce larger general shear strain amplitude (γ_ga), and the correlation between the u_e to ${p^{'}}_{0}$ and γ_ga can be expressed approximately as tangent function. Based on the experimental data, the number of cycles (N_L) required to induce initial liquefaction decreases with the increasing α₀ for a given cyclic stress ratio (CSR). However, by defining the modified cyclic stress ratio (CSR_M) as a new index for cyclic stress level, a virtually unique relationship between CSR_M and N_L can be established for all α₀ cases. The applicability of CSR_M is verified by the original test data of saturated sands from the literatures.
- saturated coral sand /
- initial major stress direction angle /
- excess pore water pressure /
- deformation characteristic /
- cyclic resistance /
- modified cyclic stress ratio

HTML全文

0. 引言

非饱和渗透系数是非饱和土水力特性研究的重要参数之一^[1]，它是基质吸力（或体积含水量）的函数，变化范围可达数个数量级。目前，室内试验测定非饱和渗透系数，耗时费力，尤其是在高吸力阶段，测量周期较长。因此，学界建立了不同的模型来间接预测非饱和土的渗透系数。这些模型大致可分为经验模型^[2]、宏观模型^[3]和统计模型^[4-5]等。经验模型和宏观模型一般采用相对简单的数学公式对渗透系数进行预测，需要较多试验数据以获得最佳拟合参数，不便于复杂条件下的非饱和土水力特性理论分析。统计模型往往利用土-水特征曲线（SWCC）确定非饱和土渗透系数，是目前比较流行的方法，常见的基本统计模型有Childs Collis-George（CCG）模型^[4]，Mualem^[5]，Burdine模型^[6]等，许多学者对模型进行了修正，以便更准确地预测非饱和渗透系数。Van Genuchten（VG）^[7]建立了土-水特征曲线模型，并结合Burdine和Mualem提出的传导率模型建立了相对渗透系数方程。Fredlund等^[8]提出了一种新的SWCC模型，并结合CCG模型推导了整个吸力范围内的相对渗透系数函数方程。陶高梁等^[9-10]将SWCC视为反映土体不同尺寸孔隙渗透通道的间接指标，建立了考虑微观渗透通道起始水力梯度的饱和黏性土非线性渗透模型；基于渗透通道分形特性，建立了不同水力梯度下的饱和黏性土非线性渗透分形新模型。

除了基于SWCC模型的方法外，陈正汉等^[11]、姚志华等^[12]结合自主设计的原状黄土取样设备对原状和重塑黄土的渗水特性进行了一维水平土柱入渗试验，形成了成套的直接取样测量方法，该方法过程简明，原理清晰，所得结果较为可靠。

值得一提的是，直接测量方法一般试样尺寸较大，试验过程耗时相对较长，工作量较大；对于基于SWCC的间接预测方法，由于需要获得SWCC，同样需要较长的测量时间与较大的工作量。因此，提出一种可快速准确预测非饱和渗透系数的方法具有重要的实际意义和工程应用价值。

核磁共振（NMR）技术是一种灵敏度高、测试速度快、无破坏性的先进水分探测技术，核磁共振技术可应用于岩土工程领域中土料水分和孔隙变化的准确测量，具有测量时间快的优点，需要指出的是，NMR方法设备昂贵，相关成本较高。

目前利用NMR技术探究土壤中孔隙水分布和迁移的微观机理的研究已有一定基础^[13]。在多孔介质渗透特性预测方面，Seevers^[14]推导了基于弛豫时间原理的砂岩饱和渗透率方程。在此基础上，Kleinberg等^[15]提出了Schlumberger-Doll Research（SDR）方程，并采用平均弛豫时间作为预测岩土饱和渗透系数的参数。然而，现有的NMR试验经典模型主要通过连通孔隙的宏观总体积来描述岩石的渗透率，并未充分考虑核磁共振中横向弛豫时间分布几何形态等信息。综上所述，上述模型大多依赖于经验模型来推导渗透率函数，适用范围较小，并且预测对象都是饱和土，对于非饱和土渗透系数的预测较少，而本文的研究重点是通过NMR技术提出一种能够快速预测非饱和渗透系数的新方法。

结合NMR理论和渗流理论，本文提出了一种基于NMR曲线快速预测饱和/非饱和渗透系数的方法。为验证该方法的准确性，以湖南黏土为例，制备不同初始孔隙比土样，采用瞬时剖面法得到了非饱和渗透系数与体积含水量的关系；利用核磁共振技术测量了不同初始孔隙比试样脱湿及加湿过程的NMR曲线，分析了基于试样不同状态下NMR曲线的非饱和渗透系数预测结果，并与已有基于SWCC的预测模型预测结果进行了比较。

1. 基于NMR曲线的饱和/非饱和渗透系数新模型

假设土样由不同孔径的连通孔隙通道组成，土料中水的渗流运动皆发生于这些连通的孔隙通道中。饱和土的渗透系数即为海量连通孔隙通道饱和渗透系数的叠加值，可表示为^[16]

${k}_{\text{s}}=\frac{Q}{AJ}=\frac{{\displaystyle \sum _{i=1}^{n}\frac{\gamma {d}_{i}{}^{2}}{32\mu }\times {A}_{i}J}}{AJ}=\frac{{\displaystyle \sum _{i=1}^{{}_{n}}\frac{\gamma {d}_{i}{}^{2}}{32\mu }\times {A}_{i}}}{A}\text{ }。$

(1)

式中：Q为单位时间流量（cm³/s）；J为水力梯度（cm）；A为分析土样的横截面总面积（cm²）；μ为黏度（Pa·s）； $\gamma$ 为水的重度（N/m³）；d_i（cm）为第i级孔隙通道的直径；A_i（cm²）为第i级孔隙通道的横截面积。

假设第i级孔隙通道的实际长度与土样长l（cm）之比为p_i，则孔隙通道实际长度可记为p_il。若第i级孔隙通道的所含水分对应的体积含水量为θ_i，则相应的水分总体积为θ_iV_T，其中V_T（m³）为土样的总体积。则第i级孔隙通道截面积为

${A_i} = \frac{{{\theta _i}{V_{\text{T}}}}}{{{p_i}l}} 。$

(2)

将式（2）代入式（1）中可得

${k_{\text{s}}} = \sum\limits_{i = 1}^{_n} {\frac{{\gamma {d_i}^2}}{{32\mu }} \times \frac{{{\theta _i}{V_{\text{T}}}}}{{{p_i}lA}}} = \sum\limits_{i = 1}^{_n} {\frac{{\gamma {d_i}^2}}{{32\mu }} \times \frac{{{\theta _i}}}{{{p_i}}}} 。$

(3)

通过核磁共振（NMR）技术可获得土样的弛豫时间（T₂）分布曲线，该曲线反映了信号幅度和弛豫时间的对应关系。已有研究^[17]表明，弛豫时间T₂与土料孔隙特性直接相关：

$\frac{1}{{{T_2}}} \approx {\rho _2}\frac{S}{V} = \frac{{\lambda {\rho _2}}}{d} 。$

(4)

式中：S，V分别为水分所处孔隙的表面积（cm²）与体积（cm³）； ${\rho _2}$ 为横向弛豫率（Hz），其与土颗粒表面的物理化学性质有关；d为孔隙直径（cm），与式（3）中d_i都表示孔径， $\lambda$ 为与孔隙形状相关的因子，柱状孔隙为4，球形孔隙为6。

将式（4）简化后可得

$d = \lambda {\rho _2}{T_2} 。$

(5)

由上式可知，横向弛豫时间T₂与孔隙直径d成正比，因此弛豫时间T₂分布曲线实际上反应了孔隙分布特性。将T₂分布曲线的横坐标从右至左依次划分为n，n-1，n-2， $\cdots$ ，2级，1级，假设大小相近的第i级孔隙组成连续孔隙通道，其等效孔径为d_i（n级孔隙最大、1级孔隙最小），通过标定T₂分布曲线总信号幅度为土样总体积含水量，可得到该级通道对应的体积含水量θ_i，如图 1所示，图中编号表示了孔隙通道的分级方式。

图 1 孔隙通道分级示意图

Figure 1. Schematic diagram of classification of pore channels

下载: 全尺寸图片幻灯片

将式（5）代入式（3），得到

${k_{\text{s}}} = \sum\limits_{i = 1}^n {\frac{{\gamma {\lambda ^2}\rho _2^2{\theta _i}{T_{2i}}^2}}{{32\mu {p_i}}}} 。$

(6)

陶高梁等^[16]认为对于同一种土而言，p_i可假设为常数，于是上式可简化为

${k_{\text{s}}} = \sum\limits_{i = 1}^{_n} {{k_{\text{c}}}{\theta _i}T_{2i}^2} 。$

(7)

式中： ${k_{\text{c}}}$ = $\gamma {\lambda ^2}\rho _2^2$ /(32μ ${p_i}$ )，对于同一种土为常数。

假设总孔隙通道分为n级，现只有1~m级通道充满水（m＜n），此时试样总体积含水量为θ（θ＜θ_s，θ_s为试样饱和体积含水量），根据式（7）可得到此时非饱和渗透系数k(θ)：

$k(\theta ) = \sum\limits_{i = 1}^m {{k_{\text{c}}}{\theta _i}{T_{2i}}^2} 。$

(8)

结合式（7），（8），非饱和相对渗透系数k_r可以表示为

${k_{\text{r}}} = {{\sum\limits_{i = 1}^m {{\theta _i}{T_{2i}}^2} } \mathord{\left/ {\vphantom {{\sum\limits_{i = 1}^m {{\theta _i}{T_{2i}}^2} } {\sum\limits_{i = 1}^n {{\theta _i}} }}} \right. } {\sum\limits_{i = 1}^n {{\theta _i}} }}{T_{2i}}^2 。$

(9)

综上所述，基于NMR曲线，本文建立的饱和土渗透系数预测模型如（7）式所示，建立的非饱和土绝对渗透系数和相对渗透系数模型分别如式（8），（9）所示。

2. 湖南黏土饱和/非饱和渗透试验与NMR试验

2.1 饱和土渗透试验

本试验研究对象为湖南邵阳某地非饱和黏性土，属于原生黏土经过再次搬运堆积形成的次生黏土，土粒相对质量密度为2.76，液限w_L为46.3%，塑限w_p为27.8%。

将土料风干后碾碎，过2 mm圆孔筛，配置含水率一定的土料，选用直径为45 mm，高度为2 mm的环刀制样，制备初始孔隙比为1.12，1.04，0.97，0.90，0.84的试样，抽真空饱和，采用变水头法测量其饱和渗透系数。为提高试验精度，经多次重复试验取得平均值，进行温度修正后，得到20℃条件下不同初始孔隙比土样的饱和渗透系数^[18-19]，如表 1所示。

表 1 不同初始孔隙比试样的饱和渗透系数（20℃）

Table 1. Saturated permeability coefficients of soil samples with different initial void ratios (20℃)

初始孔隙比	1.12	1.04	0.97	0.90	0.84
渗透系数k_s/(cm·s^-1)	7.72×10⁻⁴	4.15×10⁻⁴	2.49×10⁻⁴	1.73×10⁻⁴	7.63×10⁻⁵

下载: 导出CSV

| 显示表格

2.2 非饱和土渗透试验

采用瞬时剖面法测量非饱和渗透系数。试验装置为课题组自主设计的有机玻璃桶，直径23 cm，高1 m，沿桶周均匀钻打五列竖向排列的圆形孔，孔与孔之间的间距为5 cm，直径为1 cm，并用A、B、C、D、E进行标记。试验开始前将这些小孔封闭以防漏水，本试验采用土仍为2.1节所述湖南黏土。在该试验的准备阶段，配置含水率为19%的土样待用，待水分迁移均匀后复测含水率，采用特制击实器分层击实，并在土柱顶部铺上8 cm厚的细砂层，随后从土柱顶部持续加水20 min，共计1500 mL，使水分在整个界面下渗均匀，观察湿润锋动态变化特点；分不同时间间隔于打好的孔中取土烘干以实测A、B、C、D、E列不同深度土样的含水率。结合所测试样的土-水特征曲线^[18]，计算得到非饱和渗透系数和体积含水量的关系^[18-19]。

将图 2中非饱和渗透系数实测值除以表 1中饱和渗透系数实测值，得到非饱和相对渗透系数实测值并用于后续验证。值得指出，瞬态剖面法试验中，上部土料以脱湿为主，下部土料以吸湿为主，研究表明土料脱湿与吸湿过程非饱和渗透系数存在一定差异，但相对于含水率的影响，该因素可以忽略不计^[8]。故一般情况下本文不考虑非饱和土渗透系数的滞回效应。

图 2 湖南红黏土非饱和渗透系数实测数据

Figure 2. Measured data of unsaturated permeability coefficient of Hunan clay

下载: 全尺寸图片幻灯片

2.3 核磁共振试验

（1）试验设备及原理

本文核磁共振试验于中科院武汉岩土所进行，所用仪器为苏州纽迈公司研发的PQ-001低磁场核磁共振仪。该试验过程中，试样置于试管中，其有效测试范围是60 mm×Φ60 mm。

核磁共振（NMR）的信号来源主要为氢质子，氢质子越多，表明含水率越大，反之则越低。因此通过信号量定标的方法，该技术可以用来测量物质中水的质量。NMR通过测定土样弛豫时间与信号强度的关系，计算峰面积，进而获得土样的水分尺度分布及孔隙分布特性。为了消除磁体磁场的非均匀性对T₂测量的影响，运用NMR测量分析软件及CPMG序列采集样品T₂衰减曲线，运用反演软件得到NMR试验数据。

（2）试验方案

本试验同样采用2.1节所述湖南黏土，制作两组平行试样，每组试样分别包含初始孔隙比为1.12，1.04，0.97，0.90，0.84的5个试样，用于脱湿与吸湿试验，分别模拟上层土料脱湿与下层土料吸湿过程。脱湿过程中，对饱和试样进行称重并记录，将不同初始孔隙比试样放在小铝盒中，自然风干。开始时，需每小时一次对试样进行称重，避免结果差异过大。土样脱水速度逐渐变慢后，可逐渐增大测量间隔。吸湿过程中，将初始与饱和含水率的间隔分为十级，每级之间的差相同以确定每级加湿目标含水率。由初始含水率开始，用滴水球分级向试样加水，密封72 h使水分均匀分布，最终使每个试样达到该初始含水率下的饱和含水率。利用核磁共振测试系统分别测定不同初始孔隙比土样在吸湿和脱湿过程中不同含水率试验的NMR曲线，NMR试验对应试样的含水率如表 2所示，NMR试验共计95次。

表 2 脱湿及吸湿过程NMR试验方案及其对应的试样含水率

Table 2. Moisture contents of samples corresponding to NMR tests of desorption and absorption process

孔隙比	脱湿过程					吸湿过程
孔隙比	1.12	1.04	0.97	0.90	0.84	1.12	1.04	0.97	0.90	0.84
体积含水量/%	53.16	52.35	50.39	49.18	48.36	19.85	20.61	21.38	22.14	22.91
	46.49	45.62	43.69	42.62	41.91	22.87	23.45	24.36	24.48	25.62
	40.89	40.16	38.44	36.74	36.53	29.17	26.43	27.15	27.49	28.07
	36.49	35.59	34.38	32.58	32.33	29.77	29.85	30.30	30.01	30.84
	36.15	35.21	34.08	32.23	31.98	32.75	32.85	33.28	32.70	33.65
	30.67	29.09	28.10	26.74	26.75	36.50	36.38	36.34	35.44	36.36
	29.69	28.05	27.13	25.69	25.62	40.57	39.62	39.40	38.21	39.23
	25.62	23.92	23.25	22.07	21.90	43.33	42.63	42.30	40.86	41.78
	14.35	13.64	14.36	13.17	14.13	46.72	45.56	45.26	43.50	44.52
	11.91	11.64	12.22	11.49	12.00

下载: 导出CSV

| 显示表格

（3）试验结果

a）不同初始孔隙比试样脱湿过程NMR曲线

按上节试验方案进行脱湿过程的核磁共振试验，获得不同初始孔隙比核磁共振试验反演后获得的核磁共振曲线，由于篇幅限制，本节选用e=1.12与e=0.84两种初始孔隙比结果作为代表，如图 3所示。

图 3 脱湿过程中核磁共振分布曲线试验数据

Figure 3. Experimental data of NMR distribution curves during desorption process

下载: 全尺寸图片幻灯片

由图 3可知，从饱和含水率到残余含水率过程中，曲线所围成的峰面积逐渐减小，表明含水率逐渐减少。不同初始孔隙比试样的NMR曲线变化规律存在相似之处，即均存在两个峰值，第一个峰值在0.5~0.8 ms处，第二个峰值在30~80 ms处。在脱湿过程中，随着含水率的降低，第二峰逐渐消失，说明在脱湿过程前期，大孔隙优先排水；而后第一峰值逐渐降低，表明在脱湿过程后期以小孔隙排水为主。

b）不同初始孔隙比试样吸湿过程NMR曲线

按上节试验方案进行吸湿过程的核磁共振试验，同样以e=1.12与e=0.84两种初始孔隙比反演后获得的核磁共振曲线为代表，如图 4所示。

图 4 吸湿过程中核磁共振分布曲线试验数据

Figure 4. Experimental data of NMR distribution curves in absorption process

下载: 全尺寸图片幻灯片

可以看出，与脱湿过程类似，不同初始孔隙比的NMR曲线所包围的峰面积随着含水率的增加而增加，不同初始孔隙比试样的NMR曲线同样有两个峰值，且第一个峰值的位置与脱湿过程大致相同，第二峰值点在20~30 ms。随着含水率的增加，第一峰值逐渐升高；然而不同初始孔隙比试样的NMR曲线在后五级吸湿过程中，其曲线峰值几乎不变。第二峰在前期吸湿过程中并未出现，直到含水率增加至较大含水率时，峰值略有起势，之后随含水率上升逐渐升高。这一试验结果表明，在吸湿过程中，水分优先进入小孔隙，由于小孔隙数量远多于大孔隙，因此小孔隙含水率总和大于大孔隙。NMR曲线表现为：第一峰面积明显大于第二峰面积，随着初始孔隙比的减小，大孔隙体积压缩量高于小孔隙，使得二者间差值增大。

c）不同初始孔隙比试样饱和NMR曲线

图 5为不同初始孔隙比下饱和试样的核磁共振反演曲线。

图 5 饱和状态试样NMR曲线

Figure 5. NMR curves under saturated state

下载: 全尺寸图片幻灯片

随着初始孔隙比的减小，NMR曲线所围成的总面积逐渐减小，即饱和土样的含水率减少。饱和试样的NMR曲线同样有两个峰值，大约分别位于0.8 ms和30~50 ms。随着初始孔隙比的减小，第二峰峰面积明显减小，而第一峰峰面积略有增加。上述现象表明，随着土样初始孔隙比的降低，大孔隙体积急剧减小，导致孔隙总体积减小，孔隙水质量随之减少。

3. NMR非饱和渗透系数模型验证

利用NMR试验预测渗透系数，关键是通过弛豫时间T₂分布曲线获得各孔隙通道对应的含水率。将不同初始孔隙比试样的总含水率除以各自的总信号幅度便可确定不同初始孔隙比试样中单位信号幅度对应的含水率, 随后便可将各弛豫时间 $T_{2i}^{}$ 对应的信号幅度转变为相应的体积含水量 ${\theta _i}$ 。

将中不同初始孔隙比试样NMR曲线的纵坐标分别求和，得到总信号幅度，随后将不同初始孔隙比试样的饱和含水率分别除以相对应的总信号幅度，得到单位信号幅度对应的含水率，信号幅度与单位信号幅度对应的含水率的乘积即为每个纵坐标对应的含水率，如此便将每级横向弛豫时间T_2i对应的的信号幅度转换为对应 ${\theta _i}$ 。基于式（7），将每级横向弛豫时间 $T_{2i}^{}$ 平方后乘以对应 ${\theta _i}$ 求和，便得到 $\sum\limits_{i = 1}^n {{\theta _i}} T_{2i}^2$ ，将其作为横坐标，表 1的饱和渗透系数实测值作为纵坐标进行线性拟合，设置截距为0，结果如图 6所示。

图 6 不同初始孔隙比湖南黏性土k_c线性拟合

Figure 6. Linear fitting of k_c of Hunan clay with different initial void ratios

下载: 全尺寸图片幻灯片

对于湖南黏土，其计算结果k_c=0.1046 cm/s³，均方根误差为7.151×10^-5，相关系数为0.9311，说明拟合结果良好，k_c对于同种类土样来说可视为定值。图 6证明了式（7）的合理性，式（9）是基于式（7）推导而来，因而证明了式（9）的理论基础是科学的。

利用脱湿、吸湿过程及饱和试样的NMR曲线，基于式（9）分别预测非饱和相对渗透系数。值得说明的是，本文为了准确刻画土料脱、吸湿路径不同含水率状态水分分布特性，共进行了95次核磁共振试验，在利用式（9）预测非饱和相对渗透系数时，计算方法有3种，即：吸湿NMR曲线预测方法、脱湿NMR曲线预测方法和饱和NMR曲线预测方法。3种方法的具体预测计算过程分别如下：

（a）对于脱湿、吸湿核磁共振方法而言，分母部分 $\sum\limits_{i = 1}^n {{\theta _i}} T_{2i}^2$ 按照土料饱和状态NMR曲线（）进行计算，采取的方法与中横坐标的获取方式一致，分子部分 $\sum\limits_{i = 1}^m {{\theta _i}} T_{2i}^2$ 的计算方法与分母相同，只不过需要分别选取非饱和相对渗透系数对应含水率的NMR曲线（图 3与）计算得到的 ${\theta _i}$ 和 $T_{2i}^{}$ 。

（b）对于饱和NMR曲线预测方法，同样是利用式（9）进行计算，分母 $\sum\limits_{i = 1}^n {{\theta _i}} T_{2i}^2$ 仍然按照饱和状态的NMR曲线进行同样的计算，不同的是，在计算分子 $\sum\limits_{i = 1}^m {{\theta _i}} T_{2i}^2$ 时，假设土料排水为大孔先排出，再排小孔的水（图 3脱湿NMR曲线也证明了这一现象），因此在计算某一特征含水率θ对应的非饱和相对渗透系数时，以饱和状态NMR曲线为基础，如图 1所示，根据NMR数据所得的T₂弛豫时间分布将孔隙通道分为1~n级，若θ值与第1~i级孔隙的累计体积含水量相等（或处于第1~（i-1）级孔隙的累计体积含水量与第1~i级孔隙的累计体积含水量之间），则分子 $\sum\limits_{i = 1}^m {{\theta _i}} T_{2i}^2$ 计算时只计算饱和NMR曲线1~i级孔隙的对应数值，如此便得到利用饱和NMR曲线的不同特征含水率下的非饱和相对渗透系数预测值。

将脱湿、吸湿及饱和状态NMR曲线的非饱和相对渗透系数式（9）预测值和实测值进行比较，见图 7。

图 7 式（9）非饱和渗透系数预测值与实测值对比

Figure 7. Comparison between predicted values from three kinds of NMR curves and experimental results

下载: 全尺寸图片幻灯片

由图 7可知，基于脱湿、吸湿过程及饱和状态NMR预测值整体上与实测值基本吻合，吸湿和脱湿过程NMR曲线预测值整体低于实测值。其中，当含水率大于35%时，吸湿过程NMR曲线预测误差最大，其预测值低于实测值、脱湿过程NMR曲线预测值和饱和NMR曲线预测值；当含水率达到35%左右时，脱湿过程NMR曲线预测值要低于实测值和饱和状态NMR曲线预测值，并与吸湿过程NMR曲线预测值相近，与实测值相差较远，此时饱和状态NMR曲线预测值更接近实测值。整体上看，饱和状态NMR曲线预测值误差最小，其预测值与实测值总体较为吻合，尤其在初始孔隙比为1.04时基本重合。

在含水率大于35%时，瞬态剖面法中土料一般处于上部土体脱湿过程中，因此脱湿过程NMR预测方法所得误差较小。在吸湿过程中，滴水球滴入试样中的水分不能完全被试样吸收，部分水分会在空气中挥发；而在密封迁移水分的过程中，部分水分会残留这密封膜上，导致试样的实际含水率偏小。在脱湿过程中，随着含水率的降低，试样慢慢收缩，逐渐与环刀脱离，试验过程中需要对试样进行频繁称重，导致部分土颗粒脱落，形成质量误差。对于饱和状态NMR曲线预测方法而言，该方法并没有上述试验操作过程，因此避免了试验操作带来的各项误差，对于其深层次的机理原因，有待进一步深入研究与分析。

4. 分析与讨论

4.1 核磁共振曲线缺级的影响

由2.3节（3）中的T₂分布曲线可知，较为密实的饱和试样及低含水率的非饱和试样容易出现孔隙水含量缺级，如图 8所示，连续NMR曲线出现中断的现象被称为缺级。

图 8 核磁共振曲线缺级示意图

Figure 8. Schematic diagram of missing grade of NMR curves

下载: 全尺寸图片幻灯片

此时，大孔隙水分总含量较少，大孔隙间很难形成直接的联通孔隙，需借道小孔隙形成联通孔隙通道，因此连通孔隙通道的等效孔径由小孔隙孔径控制。本文将这部分水分所对应的孔隙通道称为“非连续孔隙通道”，其对土样的渗透系数不起控制作用。若按照式（7）~（9）计算，会产生较大误差，应舍弃这些数据。图 9分别分析了在上述e=1.12试样中，脱湿NMR曲线和吸湿NMR曲线预测非饱和相对渗透系数时，NMR曲线缺级的影响。

图 9 NMR曲线缺级对脱湿和吸湿过程非饱和渗透系数预测的影响

Figure 9. Effects of lack of level in NMR curves on predicted results during desorption and absorption processes

下载: 全尺寸图片幻灯片

可以明显看出，若不舍弃非连续孔隙通道（NMR曲线中缺级部分）的数据，预测值将非常不稳定，并且从10^-5至10²数量级之间变化，总体上远远偏离实测值；当舍弃非连续孔隙通道数据时，预测结果均处于10^-4~1数量级区间，与实测值吻合较好。因此在模型验证部分，非连续孔隙通道数据，本文均予以舍弃。

4.2 脱湿、吸湿及饱和状态NMR曲线预测方法的比较

非饱和渗透系数可通过脱湿、吸湿和饱和试样的NMR曲线来预测。虽然这3种情况下的预测值都与试验值吻合较好，但整体上来看基于饱和状态下的NMR曲线总体预测效果更好。另一方面，采用脱湿或吸湿过程的NMR曲线对某一含水率条件下的非饱和渗透系数进行预测时，需要先对相应含水率试样进行NMR试验测量相应NMR曲线，预测不同含水率下的非饱和渗透系数则需要进行多次试验，试验周期长，工作量大，成本高；而饱和NMR曲线预测方法仅需要一次核磁共振试验便可预测全含水率区间的所有非饱和相对渗透系数，其工作量和成本远低于其他两种方法。因此，本文建议采用饱和NMR曲线直接预测饱和/非饱和渗透系数。

4.3 NMR曲线预测方法与已有方法对比

目前为止，基于土-水特征曲线（SWCC）预测非饱和渗透系数是最常用的方法，然而，SWCC本身测量周期较长，获取全吸力范围的试验数据较为困难，不利于实际工程应用。值得一提的是，核磁共振（NMR）方法具有快速方便的优点，值得进一步研究。前文已对NMR方法的有效性进行了论证，其中饱和NMR曲线预测方法是最值得推广的方法。为进一步比较SWCC预测方法与NMR预测方法，本文选取了两种SWCC预测方法：CCG（Childs Collis-George）模型和文献[16]模型，对湖南黏土非饱和相对渗透系数进行预测。图 10给出了饱和NMR预测法、CCG模型预测法和文献[16]模型预测法结果与非饱和相对渗透系数实测值。

图 10 NMR预测方法与已有SWCC预测方法对比

Figure 10. Comparison of prediction methods based on NMR curves and SWCCs

下载: 全尺寸图片幻灯片

由图 10可知，CCG模型预测值均高于实测值，且预测误差较大，土样含水率越小，预测误差越大，当含水率接近30%时，预测值偏离较大，可见CCG模型对湖南黏土适用性不佳；在含水率高于35%时，文献[16]与NMR预测方法预测结果接近，且与实测值吻合较好，而在体积含水量低于35%时，文献[16]模型预测值远低于实测值，而NMR预测方法仍具有较好预测效果。

本文提出的饱和NMR曲线预测方法能够快速与较为准确地预测湖南黏土的非饱和相对渗透系数，但对于其他土类的适用性仍有待验证。同时，对存在缺级现象的土料，直接应用时误差可能较大，应舍弃缺级部分。值得一提的是，在应用式（9）预测非饱和相对渗透系数时不需要获得k_c，若需要应用式（7）或式（8）获得非饱和/饱和（绝对）渗透系数，需要对不同初始孔隙比的试样进行多次饱和渗透试验与核磁共振试验以确定k_c，增加了试验的工作量。此外，由图 10可以看出，SWCC预测方法（CCG、文献[16]模型）无法给出低含水率区间的非饱和相对渗透系数预测结果，其原因在于，SWCC通常由压力板仪测得，一般无法给出高基质吸力（即低含水率）区间的实测结果，而饱和NMR预测方法可以同时给出全含水率范围内的渗透系数预测结果。尽管如此，饱和NMR预测方

法在低含水率区间的预测效果仍需进一步验证。

5. 结论

结合核磁共振（NMR）理论和渗流理论建立了基于核磁共振曲线的饱和/非饱和渗透系数预测模型。为验证该方法的合理性，以湖南黏土为研究对象，进行了渗透试验及核磁共振试验，利用试验数据验证了该模型的合理性，得到以下4点结论。

（1）首次提出了通过核磁共振曲线直接预测饱和/非饱和渗透系数的模型（式（7）~（9）），并结合湖南黏土的饱和/非饱和渗透系数实测值进行了验证，证明该模型可靠有效。

（2）脱湿、吸湿及饱和NMR曲线均能较好预测湖南黏性土非饱和相对渗透系数，其中饱和NMR曲线预测效果最佳。此外，饱和NMR曲线预测方法具有成本低、工作量小、速度更快等优点，本文建议采用饱和NMR曲线预测方法预测非饱和渗透系数。

（3）对于缺级的NMR曲线，非连续孔隙通道数据对预测结果具有较大影响，计算时应舍弃NMR曲线缺级之后的大孔隙试验数据，即“非连续孔隙通道”数据。

（4）采用CCG模型、文献[16]模型两种基于土-水特征曲线的预测方法对湖南黏性土非饱和渗透系数进行了预测，并与本文提出的饱和NMR曲线预测方法的预测结果、湖南黏土非饱和渗透系数实测结果进行对比，发现本文提出的饱和NMR法的预测结果与实测值吻合更好。

图 1 空心圆柱扭剪仪（HCA）

Figure 1. Hollow cylinder torsional apparatus

下载: 全尺寸图片幻灯片

图 2 珊瑚砂颗粒电镜扫描图像

Figure 2. Electronic microscope scanning graphs of coral sand

下载: 全尺寸图片幻灯片

图 3 珊瑚砂级配曲线

Figure 3. Grain-size distribution curve of coral sand

下载: 全尺寸图片幻灯片

图 4 空心圆柱试样受力状态

Figure 4. Hollow cylinder sample under loads

下载: 全尺寸图片幻灯片

图 5 试验中应力路径示意图

Figure 5. Diagram of stress paths in tests

下载: 全尺寸图片幻灯片

图 6 不同 $α_{0}$ 时饱和珊瑚砂超静孔压的发展

Figure 6. Development tendency of excess pore water pressure at various $α_{0}$ for saturated coral sand

下载: 全尺寸图片幻灯片

图 7 不同 $α_{0}$ 时饱和珊瑚砂应变分量的发展

Figure 7. Development tendency of strain components at various $α_{0}$ for saturated coral sand

下载: 全尺寸图片幻灯片

图 8 不同 $α_{0}$ 时饱和珊瑚砂广义剪应变幅值的发展

Figure 8. Development tendency of generalized shear strain amplitude γ_ga at various $α_{0}$ for saturated coral sand

下载: 全尺寸图片幻灯片

图 9 不同 $α_{0}$ 时饱和珊瑚砂广义剪应变与超静孔压比的关系

Figure 9. Relationships between amplitudes of generalized shear strains and excess pore water pressure ratios at various $α_{0}$ for saturated coral sand

下载: 全尺寸图片幻灯片

图 10 参数a,b与 $α_{0}$ 之间的关系

Figure 10. Relationships among parameters a, b and $α_{0}$

下载: 全尺寸图片幻灯片

图 11 不同 $α_{0}$ 时饱和珊瑚砂的动强度曲线

Figure 11. Curves of cyclic resistance at various $α_{0}$ for saturated coral sand

下载: 全尺寸图片幻灯片

图 12 不同 $α_{0}$ 时N_L与CSR_M之间的关系

Figure 12. Relationships between CSR_M and N_L at various $α_{0}$ for saturated coral sand

下载: 全尺寸图片幻灯片

图 13 3种砂土的CSR、CSR_M与N_L之间的关系

Figure 13. CSR and CSR_M of three sands vs. N_L

下载: 全尺寸图片幻灯片

表 1 珊瑚砂的基本物理性质指标

Table 1 Index properties of coral sand

平均粒径d₅₀	不均匀系数C_u	曲率系数C_c	相对密度G_s	最大孔隙比e_max	最小孔隙比e_min
0.31	4.67	0.86	2.80	1.72	0.99

下载: 导出CSV

表 2 试验方案

Table 2 Summary of test schemes

试验编号	初始主应力方向角α₀	内压 ${p^{'}}_{i}$ /kPa	外压 ${p^{'}}_{o}$ /kPa	轴向应力 ${σ^{'}}_{z}$ /kPa	扭矩M/(N·m)	CSR
1	0	106.67	96.00	120.00	0	0.25
2						0.30
3						0.35
4	22.5	104.71	97.17	114.14	2.9027	0.25
5						0.30
6						0.35
7	45	100.00	100.00	100	4.1050	0.25
8						0.30
9						0.35
10	67.5	95.29	102.83	85.86	2.9027	0.20
11						0.25
12						0.30
13	90	93.33	104.00	80.00	0	0.20
14						0.25
15						0.30

下载: 导出CSV

表 3 4种砂土在不同固结状态下 ${σ^{'}}_{1 (0)}$ , ${σ^{'}}_{3 (0)}$ ,β₁,β₂的取值

Table 3 Values of ${σ^{'}}_{1 (0)}$ , ${σ^{'}}_{3 (0)}$ , β₁ and β₂ under different consolidation conditions for four sands

土类（平均大主应力）	m	$α_{0}$ /(°)	R	b	${σ^{'}}_{1 (0)}$ /kPa	${σ^{'}}_{3 (0)}$ /kPa	β₁	β₂
南沙岛礁珊瑚砂（本文）（ ${p^{'}}_{0}$ = 100 kPa）	0.77	0	1.5	0.5	130.00	100.00	1.11	0.93
		22.5			128.28	92.93	1.08	0.92
		45			120.00	80.00	1.00	0.89
		67.5			107.07	71.72	0.92	0.90
		90			100.00	70.00	0.88	0.91
南京细砂^[10]（ ${p^{'}}_{0}$ = 100 kPa）	1.42	—	1.0	—	100.00	100.00	1.00	1.00
		0	2.0	0.5	150.00	100.00	1.37	0.89
		0	1.5	0.5	130.00	100.00	1.22	0.93
		30	1.5	0.5	126.37	88.63	1.11	0.91
		45	1.5	0.5	120.00	80.00	1.00	0.89
		60	1.5	0.5	111.37	73.63	0.90	0.89
		90	1.5	0.5	100.00	70.00	0.79	0.91
南沙岛礁珊瑚砂^[20]（ ${p^{'}}_{0}$ = 100 kPa）	0.98	0	1.66	0.5	137.22	100.00	1.18	0.92
		30	1.66	0.5	132.72	85.89	1.09	0.89
		45	1.66	0.5	124.81	75.19	1.00	0.87
		60	1.66	0.5	114.11	67.28	0.91	0.86
		90	1.66	0.5	100.00	62.78	0.82	0.88
		0	1.66	0	154.10	100.00	1.26	0.89
		0	1.66	1.0	122.92	100.00	1.11	0.95
福建砂^[10]（ ${p^{'}}_{0}$ = 100 kPa）	2.49	0	1.66	0.5	137.22	100.00	1.53	0.92
		45	1.66	0.5	124.81	75.19	1.00	0.87
		90	1.66	0.5	100.00	62.78	0.60	0.88

下载: 导出CSV

参考文献(21)

[1]	朱长歧, 陈海洋, 孟庆山, 等. 钙质砂颗粒内孔隙的结构特征分析[J]. 岩土力学, 2014, 35(7): 1831-1836. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201407003.htm ZHU Chang-qi, CHEN Hai-yang, MENG Qing-shan, et al. Microscopic characterization of intra-pore structures of calcareous sands[J]. Rock and Soil Mechanics, 2014, 35(7): 1831-1836. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201407003.htm
[2]	陈海洋, 汪稔, 李建国, 等. 钙质砂颗粒的形状分析[J]. 岩土力学, 2005, 26(9): 1389-1392. doi: 10.3969/j.issn.1000-7598.2005.09.008 CHEN Hai-yang, WANG Ren, LI Jian-guo, et al. Grain shape analysis of calcareous soil[J]. Rock and Soil Mechanics, 2005, 26(9): 1389-1392. (in Chinese) doi: 10.3969/j.issn.1000-7598.2005.09.008
[3]	HYODO M, HYDE A F L, ARAMAKI N. Liquefaction of crushable soils[J]. Géotechnique, 1998, 48(4): 527-543. doi: 10.1680/geot.1998.48.4.527
[4]	MAO X, FAHEY M. Behaviour of calcareous soils in undrained cyclic simple shear[J]. Géotechnique, 2003, 53(8): 715-727. doi: 10.1680/geot.2003.53.8.715
[5]	SHARMA S S, ISMAIL M A. Monotonic and cyclic behavior of two calcareous soils of different origins[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2006, 132(12): 1581-1591. doi: 10.1061/(ASCE)1090-0241(2006)132:12(1581)
[6]	李建国. 波浪荷载作用下饱和钙质砂动力特性的试验研究[D]. 武汉: 中国科学院武汉岩土力学研究所, 2005. LI Jian-guo. Experimental Research on Dynamic Behavior of Saturated Calcareous Sand Under Wave Loading[D]. Wuhan: Institute of Rock & Soil Mechanics China Academy of Sciences, 2005. (in Chinese)
[7]	虞海珍. 复杂应力条件下饱和钙质砂动力特性的试验研究[D]. 武汉: 华中科技大学, 2006. YU Hai-zhen. Experimental Research on Dynamic Behavior of Saturated Calcareous Sand Under Complex Stress Conditions[D]. Wuhan: Huazhong University of Science and Technology, 2006. (in Chinese)
[8]	SATO K, YASUHARA K, YOSHIDA N. Effect of pre-shearing with drainage on undrained cyclic shear behaviour of dense sand[C]//Proeeedings of the Ninth International Offshore and Polar Engineering Conferenee. France, 1999: 542-547.
[9]	付磊, 王洪瑾, 周景星, 等. 主应力偏转角对砂砾料动力特性影响的试验研究[J]. 岩土工程学报, 2000, 22(4): 435-440. doi: 10.3321/j.issn:1000-4548.2000.04.010 FU Lei, WANG Hong-jin, ZHOU Jing-xing, et al. Effect of initial rotation angle of principal stress on the dynamic properties of soil[J]. Chinese Journal of Geotechnical Engineering, 2000, 22(4): 435-440. (in Chinese) doi: 10.3321/j.issn:1000-4548.2000.04.010
[10]	郭莹, 栾茂田, 许成顺, 等. 主应力方向变化对松砂不排水动强度特性的影响[J]. 岩土工程学报, 2003, 25(6): 666-670. doi: 10.3321/j.issn:1000-4548.2003.06.004 GUO Ying, LUAN Mao-tian, XU Chen-shun. et al. Effect of variation of principal stress orientation on undrained dynamic strength behavior of loose sand[J]. Chinese Journal of Geotechnical Engineering, 2003, 25(6): 666-670. (in Chinese) doi: 10.3321/j.issn:1000-4548.2003.06.004
[11]	潘华. 复杂应力条件下饱和南京细砂动力特性试验研究[D]. 南京: 南京工业大学, 2011. PAN Hua. Experimental Research on Dynamic Behavior of Saturated Nanjing Fine Sand under Complex Stress Conditions[D]. Nanjing: Nanjing University of Technology, 2011. (in Chinese)
[12]	HIGHT D W, GENS A, SYMES M J. The development of a new hollow cylinder apparatus for investigating the effects of principal stress rotation in soils[J]. Géotechnique, 1983, 33(4): 355-383. doi: 10.1680/geot.1983.33.4.355
[13]	YANG Z X, LI X S, YANG J. Undrained anisotropy and rotational shear in granular soil[J]. Géotechnique, 2007, 57(4): 371-384. doi: 10.1680/geot.2007.57.4.371
[14]	CHEN G X, WU Q, ZHOU Z L, et al. Undrained anisotropy and cyclic resistance of saturated silt subjected to various patterns of principal stress rotation[J]. Géotechnique, in press.
[15]	ISHIHARA K, YAMAZAKI A. Analysis wave-induced liquefaction in seabed deposits of sand[J]. Soils and Foundations, 1984, 24(3): 85-100. doi: 10.3208/sandf1972.24.3_85
[16]	HUANG B, CHEN X, ZHAO Y. A new index for evaluating liquefaction resistance of soil under combined cyclic shear stresses[J]. Engineering Geology, 2015, 199: 125-139. doi: 10.1016/j.enggeo.2015.10.012
[17]	SEED H B, IDRISS I M, ARANGO I. Evaluation of liquefaction potential using field performance data[J]. Journal of Geotechnical Engineering, 1983, 109(3): 458-482. doi: 10.1061/(ASCE)0733-9410(1983)109:3(458)
[18]	CASAGRANDE A. Charateristics of cohensionless soils affecting the stability of slopes and earth fill[J]. Journal of Boston Society of Civil Engineering, 1936: 257-276.
[19]	俞茂宏. 岩土类材料的统一强度理论及其应用[J]. 岩土工程学报, 1994, 16(2): 1-10. doi: 10.3321/j.issn:1000-4548.1994.02.001 YU Mao-hong. Unified strength theory for geomaterials and its applications[J]. Chinese Journal of Geotechnical Engineering, 1994, 16(2): 1-10. (in Chinese) doi: 10.3321/j.issn:1000-4548.1994.02.001
[20]	李建国, 汪稔, 虞海珍, 等. 初始主应力方向对钙质砂动力特性影响的试验研究[J]. 岩土力学, 2005, 26(5): 723-727. doi: 10.3969/j.issn.1000-7598.2005.05.009 LI Jian-guo, WANG Ren, YU Hai-zhen, et al. Experimental research on effect of initial principal stress orientation on dynamic properties of calcareous sand[J]. Rock and Soil Mechanics, 2005, 26(5): 723-727. (in Chinese) doi: 10.3969/j.issn.1000-7598.2005.05.009
[21]	CHEN G X, WU Q, SUN T, et al. Cyclic behaviors of saturated sand-gravel mixtures under undrained cyclic triaxial loading[J]. Journal of Earthquake Engineering, in press.