压实黄土在固结和剪切过程中的微结构演化

肖涛; 李萍; 邵生俊

doi:10.11779/CJGE2021S1018

压实黄土在固结和剪切过程中的微结构演化 English Version

肖涛^1,,
李萍^{1, 2, ,},
邵生俊²

1.
西北大学大陆动力学国家重点试验室地质学系,陕西　西安　710069
2.
西安理工大学陕西省黄土力学与工程重点试验室,陕西　西安　710048

基金项目:

黄土力学与工程陕西省重点试验室开放基金项目 LME201803

中国博士后基金项目 2019M653883XB

详细信息

作者简介:
肖涛（1994— ）,男,硕士研究生,从事黄土微结构与力学特性研究。Email：672590412@qq.com

通讯作者:
李萍, E-mail：lp19881028@163.com

中图分类号: TU43
计量
- 文章访问数: 179
- HTML全文浏览量: 28
- PDF下载量: 221
出版历程
- 收稿日期: 2020-12-28
- 网络出版日期: 2022-12-05
- 刊出日期: 2021-06-30

Microstructural evolution of compacted loess during consolidation and shearing

XIAO Tao^1,,
LI Ping^{1, 2, ,},
SHAO Sheng-jun²

1.
State Key Laboratory of Continental Dynamics, Department of Geology, Northwest University, Xi'an 710069, China
2.
Shaanxi Key Laboratory of Loess Mechanics and Engineering, Xi'an 710048, China

摘要

摘要: 为了解压实黄土在固结和剪切过程中微结构的演化过程,运用三轴试验和MIP技术,得出不同应力应变状态的压实黄土试样的孔径分布曲线,通过分析压实土样的孔径分布特征,研究压实黄土在固结和剪切过程中的微结构演化规律。结果显示,该压实黄土在剪应力作用下表现为强烈的剪缩特性,固结和剪切阶段的孔径分布演化规律相似,粒内孔隙没有受到影响,只有粒间孔隙被压缩。另外,基于孔径分布曲线,利用分形理论计算得出的试样孔隙体积分形维数变化不大,表明在固结和剪切过程中,孔隙表面的粗糙程度变化不大,且试样的分形维数和孔隙比、轴向应变之间无明显相关关系。研究结果说明了压实黄土的力学响应是集粒间相互作用的结果,而非颗粒与颗粒之间相互作用的结果。
- 压实黄土 /
- 微结构演化 /
- 孔径分布曲线 /
- 力学响应 /
- 分形维数
Abstract: To understand the microstructural evolution of compacted loess during consolidation and shearing, the triaxial tests and the MIP technology are used to acquire the pore-size distribution curve of compacted loess specimens with different stress and strain states.By analyzing the pore-size distribution characteristics of compacted loess specimens, the microstructural evolution of compacted loess during consolidation and shearing is studied.The results show that the compacted loess tested is strongly contractive during shearing, and the evolution of the pore-size distribution of compacted loess during consolidation is similar to that during shearing.The intra-aggregate pores are not affected, only the inter-aggregate pores are compressed.In addition, based on the pore-size distribution curve, the fractal dimension of each specimen is calculated by the fractal theory, which changes a little from specimen to specimen.It is indicated that the roughness of pore surface changes little during the process of consolidation and shearing, and there is no obvious correlation between fractal dimension, void ratio and axial strain.The research results indicate that the mechanical response of compacted loess is due to the interaction between aggregates, while not between particles.
- compacted loess /
- evolution of microstructure /
- pore-size distribution curve /
- mechanical response /
- fractal dimension

HTML全文

0. 引言

山岭隧道的洞口段作为隧道施工的“咽喉”,常常分布着节理裂隙较发育的严重风化岩体,施工中极易出现围岩失稳坍塌、地表塌陷等工程灾害^[1-2]。鉴于位移和变形是反映工程安全控制最核心的元素,因此,对洞口段地表沉降进行预测研究,对于保障隧道进洞安全及塌方灾害防治具有重要现实意义。

近年来,计算机和人工智能技术飞速发展,智能预测方法在路基沉降^[3-4]和滑坡位移^[5-6]预测中已有较多研究。而对于地表沉降的预测研究相对较少,主要采用的方法有神经网络模型^[7-8]、贝叶斯模型^[9]、粒子群-最大似然模型^[10]等。显然,对地表沉降的预测主要集中在地铁盾构隧道领域,而对山岭隧道洞口段地表沉降的预测缺乏深入研究。虽然上述智能预测方法有了一些应用,但仍存在某些方面的不足：①地表沉降监测值多为非线性和非等距时序数据,由于其复杂性,此时直接基于实测值的预测不易得到良好的结果；②预测模型（如SVM、ELM等）多为“静态”模型,只能学习并使用当前时刻的监测值,并且是一次性学习全部训练数据,以至于不能很好地反映地表沉降的长期变化趋势；③智能优化算法易陷入局部极小值,模型泛化能力变差,且“过拟合”现象较为严重,从而导致预测结果出现较大误差。

综上,本文结合非等距时序分析理论,并引入一种新的位移分解方法——变分模态分解（variational mode decomposition,VMD）^[11],对沉降实测值进行分解,然后采用具有动态预测能力的在线贯序极限学习机（online sequential extreme learning machine,OSELM）^[12]分别预测,针对模型存在的输入权值及隐含层阈值确定较为困难的缺陷,引入一种新型的元启发式智能优化算法——灰狼优化算法（grey wolf optimizer,GWO）^[13]对模型进行优化。将优化后的模型应用于重庆市兴隆隧道洞口段实测地表沉降预测中,并进行对比分析,进一步验证所提模型的合理性和实用性。

1. 地表沉降时序预测模型

1.1 监测数据预处理

根据《公路隧道施工技术规范》^[14],当隧道埋深小于2倍的洞室开挖宽度时,地表沉降的量测是十分重要的,其量测频率见表1。因此,现场实测值往往是非等间距的,不利于时序分析原理的应用,这就需要采用数学手段将其转换为等间距时序数据。3次样条函数插值法作为一种应用十分广泛的插值方法,将其应用于地表沉降监测值的等距处理中,可得到较为连续光滑的曲线,并且能够最大程度的保留原始数据的变化趋势。

表 1 地表沉降量测频率

Table 1. Measurement frequencies of ground settlement

量测项目	d	量测频率
地表沉降	d < 2B	（1~2次）/d
	2b < d < 5B	1次/（2~3 d）
	d > 5B	1次/7 d
注：d为开挖断面距量测断面的距离,B为隧道开挖宽度。

下载: 导出CSV

| 显示表格

现有研究表明,地表沉降是受基础地质条件（如地形地貌和岩性等）和外部因素（如降雨、风荷载、测量条件等）等共同作用形成的^[15]。在前者作用下表现为随时间变化的趋势项位移,反映变形的主要规律,为非平稳序列；在后者作用下则表现为随机变化的随机项位移,为平稳噪声序列,该序列的存在会降低预测真实性及精度。基于此,为了更好的反映各位移分量变化情况,利用时间序列原理,对地表沉降进行分解,

$s (t) = α (t) + β (t)$ ,

(1)

式中, $s (t)$ 为隧道地表沉降总位移, $α (t)$ 为趋势项位移, $β (t)$ 为随机项位移。

目前较为典型的位移序列分解方法主要有EMD^[16]和EEMD^[17]等,他们各有优缺点。本文尝试引入一种新的分解方法——变分模态分解^[11],该方法理论基础完善,能根据位移序列自身特征尺度将其转化为非递归VMD模式,分解所得各位移分量均具有明确物理意义,可有效克服传统EMD和EEMD存在的模态混叠问题。

实现VMD的核心是构造和求解变分问题。

（1）构造变分问题。假设原始位移序列s(t)经分解后得到了K个分量,VMD通过在变分框架内寻找模型的最优解来实现位移的分解,其变分约束表达式为

$min_{{δ_{k}}, {ω_{k}}} {{\sum_{k} ‖ \partial_{t} [(μ (t) + \frac{i}{π n}) * δ_{k} (t)] e^{- i ω_{k} t} ‖}_{2}^{2}}$ ,

(2)

$s . t . \sum_{k}^{K} δ_{k} = s (t)$ ,

(3)

式中,K为分解所得的模态个数（正整数）, ${δ_{k}}$ , ${ω_{k}}$ 分别为分解后的第K个模态分量及对应中心频率, $μ (t)$ 为狄拉克函数, $*$ 为卷积运算符。

（2）求解变分问题。引入Lagrange乘法算子,将上式转变为非约束变分问题,得到增广Lagrange表达式为

$\begin{matrix} L ({δ_{k}}, {ω_{k}}, λ) = σ {\sum_{k} ‖ \partial_{t} [(μ (t) + i / π t) * δ_{k} (t)] e^{- i ω_{k} t} ‖}_{2}^{2} + \\ {‖ s (t) - \sum_{k} δ_{k} (t) ‖}_{2}^{2} + 〈 λ (t), s (t) - \sum_{k} δ_{k} (t) 〉 \end{matrix}$ ,

(4)

式中, $σ$ 为二次惩罚因子。

利用交替方向乘子法将上述问题转化为寻找增广Lagrange方程的“鞍点”问题,也即最优解,限于篇幅,具体过程见文献[18]。本文将其引入地表沉降分解中,得到不同特征尺度的位移分量,更能反映其自身信息。

1.2 地表沉降动态预测的GWO-OSELM模型

（1）在线贯序极限学习机模型

传统神经网络模型（如SVM、ELM等）多为“静态”的,预测时只能学习并使用当前时刻监测值,并且一次性学习全部训练数据,而地表沉降监测值往往是非线性动态变化的,并且彼此之间存在一定联系。鉴于此,将具有在线学习机制的OSELM模型^[12]引入地表沉降预测中,它是对传统ELM模型的改进,能够逐批学习监测值并根据最新监测值不断更新网络,实现动态预测。

OSELM模型主要通过利用递推计算对ELM模型的输出层权值 $\hat{β}$ 进行更新,主要包括两部分。

a）初始化阶段

给定激励函数g(x),隐含层节点数n,以及初始训练子集 $ψ_{0} = {(x_{i}, y_{i}) | i = 1, 2, \dots, N_{0}}$ ,并且 $N_{0} \geq n$ ,则初始隐含层输出权值可由下式求得,

$β_{0} = P_{0} H_{0}^{T} T_{0}$ ,

(5)

式中, $P_{0} = {(H_{0}^{T} H_{0})}^{- 1}$ , $T_{0} = {[\begin{matrix} y_{1} & y_{2} & \cdot \cdot \cdot & y_{N 0} \end{matrix}]}^{T}$ , $H_{0}$ 为隐含层初始输出矩阵。

b）在线学习阶段

令m=0,m为加入训练数据的批次,对m+1批数据,计算其隐层输出矩阵H_m+1,并按下式更新隐层输出权值：

$P_{m + 1} = P_{m} - \frac{P_{m} h_{m + 1}^{T} h_{m + 1} P_{m}}{1 + h_{m + 1} P_{m} h_{m + 1}^{T}}$ ,

(6)

$β_{m + 1} = β_{m} + P_{m + 1} H_{m + 1}^{T} (T_{m + 1}^{T} - H_{m + 1} β_{m})$ ,

(7)

式中, $h_{m} = [\begin{matrix} G (a_{1}, b_{1}, x_{m}) & \cdot \cdot \cdot & G (a_{n}, b_{n}, x_{m}) \end{matrix}]$ , $G (a_{n},$ $b_{n}, x_{m}) =$ $g_{n} (x_{m})$ , $a_{n}$ ,b_n分别是第n个隐层节点的输入权值和阈值。

（2）灰狼优化算法

GWO由澳大利亚学者Mirjalili等^[13]在2014年首次提出,它模拟了自然界中灰狼种群的社会等级制度和群体狩猎行为,灰狼群体按社会等级制度高低依次分为α狼、β狼、γ狼和 $ω$ 狼,其群体狩猎行为如下：

Step 1狼群跟踪并接近猎物。按下式计算距离：

$D = | C \cdot X_{p} (u) - X (u) |$ ,

(8)

$C = 2 \times r_{1}$ ,

(9)

式中,D为灰狼与猎物间的距离, $u$ 为当前迭代次数,X_p( $u$ )为第 $u$ 代猎物位置向量,X( $u$ )为第 $u$ 代灰狼个体位置向量,C为系数向量,r₁为[0,1]间均匀分布的随机向量。

Step 2狼群追捕并包围猎物。按下式更新位置：

$X (u + 1) = X_{p} (u) - η \cdot D$ ,

(10)

$η = 2 b \times r_{2} - b$ ,

(11)

式中, $η$ 为系数向量, $r_{2}$ 为b的随机向量,b的初始值为2,随最大迭代次数增加而线性递减至0。

Step 3狼群攻击并捕杀猎物。依据狼群中 $α$ 狼寻找猎物方向,按下式计算 $ω$ 狼与各等级狼的距离：

$D_{l} = | C_{j} \cdot X_{l} (u) - X (u) |$ ,

(12)

$X_{j} = X_{l} (u) - η_{j} \cdot D_{l}$ ,

(13)

$X (u + 1) = \frac{X_{1} + X_{2} + X_{3}}{3}$ 。

(14)

式中 l=α,β,γ；j=1,2,3。

（3）GWO-OSELM模型动态预测流程

虽然OSELM模型存在诸多优点,但它同样是通过随机产生输入权值及隐层阈值来进行计算,存在诸多缺点,故引入GWO算法对其输入权值及隐层阈值进行优化,以提高模型预测精度,并改善模型泛化能力。如图1所示,基于前述原理建立了该模型,具体步骤如下：

图 1 GWO-OSELM动态预测流程

Figure 1. Flow chart of dynamic prediction of GWO-OSELM

下载: 全尺寸图片幻灯片

a）判断洞口段地表沉降监测值是否等距,若是,则利用VMD将其分解为趋势项及随机项,若否,则先利用3次样条插值将其转化为等距后再进行分解。

b）为便于后续进行最优权值及阈值的搜寻及消除数据量纲的影响,使算法更容易收敛,分别对分解得到的时序数据进行归一化操作,使其位于[-1,1]。

c）初始化算法相关参数,主要包括狼群数量、最大迭代次数、搜索变量最大值及最小值、狼群个体初始位置及初始目标函数值。

d）以模型测试集预测的均方根误差作为个体适应度计算值,将适应度最优的个体依次作为 $α$ 狼、 $β$ 狼和 $γ$ 狼,然后更新狼群位置,并更新 $α$ 狼、 $β$ 狼和 $γ$ 狼。

e）判断是否达到最大迭代次数,若否,则继续执行步骤d）,若是,则输出OSELM最优权值及阈值。

f）将结果用于模型训练和预测中,得位移分量预测结果,由时间序列原理叠加得到地表沉降总预测值。

1.3 模型预测结果评价指标

相关性系数R、均方根误差R_MSE、平均绝对误差M_AE 3项指标常用来衡量位移预测精度^[17]。因此,本文采用它们来评价模型预测效果,

$R = \frac{\sum_{j = 1}^{N} (d_{j} - \bar{d}) ({\hat{d}}_{j} - \bar{\hat{d}})}{\sqrt{\sum_{j = 1}^{N} {(d_{j} - \bar{d})}^{2} {({\hat{d}}_{j} - \bar{\hat{d}})}^{2}}}$ ,

(15)

$R_{MSE} = \sqrt{\frac{1}{N} \sum_{j = 1}^{N} {({\hat{d}}_{j} - d_{j})}^{2}}$ ,

(16)

$M_{AE} = \frac{1}{N} \sum_{j = 1}^{N} | {\hat{d}}_{j} - d_{j} |$ ,

(17)

式中, $d_{j}$ 为沉降监测数据, ${\hat{d}}_{j}$ 为沉降预测数据, $\bar{d}$ 为沉降监测值均值, $\bar{\hat{d}}$ 为沉降预测值均值,3项指标均为相对越小时预测结果越好。

2. 工程实例分析

2.1 工程概况

在建的兴隆隧道位于重庆市渝北区木耳镇,隧道左线全长2553 m,右线全长2539 m。前期地勘结果表明,隧址区域地质较稳定,总体呈单斜构造,主要穿越侏罗系中统上沙溪庙组地层。隧道出洞口段地形较平缓,缓坡主要分布第四系残坡积层粉质黏土,下伏基岩主要为泥岩和砂岩（图2）,岩体较破碎,洞口段山体覆盖层较薄,围岩自稳能力差。因此,施工过程中及时预测地表沉降变化,对于分析隧道进洞稳定性、保障施工安全及塌方灾害防治具有重要意义。

图 2 洞口段地质情况

Figure 2. Geological conditions of portal section of tunnel

下载: 全尺寸图片幻灯片

本文以隧道左线出口段ZK38+020和ZK38+015断面地表沉降实测值为研究对象,隧道采用三台阶法施工,初期支护为钢拱架+系统锚杆+27 cm厚喷射混凝土,监测断面平面位置见图3,沿隧道横断面每隔2～5 m布置监测点位置见图4,选择距隧道拱顶较近的D1-4,D1-5监测点前50 d沉降实测值进行分析。

图 3 监测点位置平面图

Figure 3. Location plan of monitoring arrangement

下载: 全尺寸图片幻灯片

图 4 监测点位置横断面图

Figure 4. Cross section of monitoring points

下载: 全尺寸图片幻灯片

2.2 洞口段监测数据预处理

将监测数据分为2组,第一组为ZK38+020断面D1-4,D1-5监测点,第二组为ZK38+015断面D1-4,D1-5监测点,以进一步验证模型的合理性。实际预测时,以前36 d的监测值为训练集,后14 d监测值为测试集,利用前4 d的监测值依次预测后1 d的监测值。

研究区段地表沉降监测值为非等距时序数据,因此利用前述原理对监测数据进行等距处理,将其转化为间距均为1 d的等距时序数据（图5）。在此基础上,利用前述时间序列分析及VMD分解原理对地表沉降总位移进行分解,首先设定模态数K=2,以位移分解残余项（即趋势项位移）的分解效果为目标,经多次试算后确定二次惩罚因子 $σ$ =2400,上升步长 $τ$ =0.03。各监测点沉降监测值分解结果如图6所示。

图 5 非等距时序数据的等距化

Figure 5. Equidistance of non-equidistant time series data

下载: 全尺寸图片幻灯片

图 6 地表沉降数据时序分解

Figure 6. Time series decomposition of ground settlement data

下载: 全尺寸图片幻灯片

2.3 模型预测结果及分析

GWO-OSELM在传统ELM基础上引入了“时间”概念,通过不断更新网络输出权值及阈值,实现动态预测。本文进行预测时,设定算法参数为：狼群数量10只,最大迭代次数300次,搜索变量上下限分别为1和0；经多次试算后确定隐层神经元个数为10,激励函数为“Sigmoid”,初始训练样本集为15,每次学习3个数据。各位移分量及累计位移预测结果分析如下：

（1）趋势项位移预测及分析

分别采用本文模型及OSELM模型和传统的ELM模型对各监测点的地表沉降趋势项进行预测,预测结果见图7,ZK38+020断面及ZK38+015断面D1-4,D1-5监测点趋势项的预测精度及误差分别见表2,3。

图 7 趋势项位移预测结果

Figure 7. Prediction results of trend terms displacement

下载: 全尺寸图片幻灯片

表 2 ZK38+020趋势项位移预测精度及误差

Table 2. Prediction accuracies and errors of trend term displacement of ZK38+020

模型	R				R_MSE/ mm				M_AE/ mm
	D1-4		D1-5		D1-4		D1-5		D1-4		D1-5
	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集
GWO-OSELM	0.9999	0.9994	0.9999	0.9998	0.0089	0.0106	0.0127	0.0051	0.0077	0.0088	0.0105	0.0044
OSELM	0.9998	0.9996	0.9998	0.9998	0.0484	0.1889	0.0566	0.2349	0.0401	0.1803	0.0506	0.2282
ELM	0.9998	0.9989	0.9996	0.9997	0.0668	0.3581	0.0881	0.4870	0.0553	0.3392	0.0725	0.4724

下载: 导出CSV

| 显示表格

表 3 ZK38+015趋势项位移预测精度及误差

Table 3. Prediction accuracies and errors of trend term displacement of ZK38+015

模型	R				R_MSE/ mm				M_AE/ mm
	D1-4		D1-5		D1-4		D1-5		D1-4		D1-5
	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集
GWO-OSELM	0.9999	0.9991	0.9999	0.9998	0.0123	0.0540	0.0261	0.0166	0.0109	0.0405	0.0217	0.0155
OSELM	0.9994	0.9990	0.9993	0.9998	0.0825	0.2426	0.0711	0.2696	0.0732	0.2363	0.0618	0.2547
ELM	0.9994	0.9992	0.9995	0.9996	0.0945	0.5241	0.0754	0.6082	0.0790	0.5060	0.0475	0.5733

下载: 导出CSV

| 显示表格

由图7中的结果对比可知,在测试集的预测中,具有动态预测能力的GWO-OSELM模型及OSELM模型的预测效果要明显优于传统的“静态”ELM模型,这主要是由于前者具有在线学习能力,能够逐批学习监测数据,并不断的更新网络,而引入GWO后的OSELM模型预测结果也要优于传统OSELM模型,原因主要是灰狼算法具有较强的寻优能力,能够得到更为准确的模型输入权值及阈值,相较于传统模型而言,优势更加明显。

再结合表2,3进一步分析,以ZK38+020断面D1-5监测点测试集预测为例,R达到了0.9998,R_MSE为0.0051 mm,M_AE为0.0044 mm,3项误差指标计算结果均以本文模型更优,这也进一步验证了模型用于趋势项位移预测的可靠性,能得到较为满意的预测结果。

（2）随机项位移预测及分析

与趋势项位移相比,随机项位移表现出了一定程度的波动及随机现象,并且具有一定的平稳性,因此,应用所提模型对随机项位移进行预测,更易体现出模型较强的预测能力。分别利用3种模型对两断面D1-4,D1-5监测点随机项位移进行预测,预测结果见图8,两断面各监测点预测精度和误差分别见表4,5。

图 8 随机项位移预测结果

Figure 8. Predicted results of random term displacement

下载: 全尺寸图片幻灯片

表 4 ZK38+020随机项位移预测精度及误差

Table 4. Prediction accuracies and errors of random term displacement of ZK38+020

模型	R				R_MSE/ mm				M_AE/ mm
	D1-4		D1-5		D1-4		D1-5		D1-4		D1-5
	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集
GWO-OSELM	0.9995	0.9839	0.9994	0.9731	0.0205	0.0162	0.0258	0.0114	0.0168	0.0137	0.0222	0.0093
OSELM	0.9938	0.9887	0.9961	0.9774	0.0748	0.1202	0.0717	0.0868	0.0523	0.1189	0.0570	0.0865
ELM	0.9982	0.9913	0.9989	0.9792	0.0440	0.1317	0.0553	0.1245	0.0369	0.1309	0.0438	0.1243

下载: 导出CSV

| 显示表格

表 5 ZK38+015随机项位移预测精度及误差

Table 5. Prediction accuracies and errors of random term displacement of ZK38+015

模型	R				R_MSE/ mm				M_AE/ mm
	D1-4		D1-5		D1-4		D1-5		D1-4		D1-5
	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集
GWO-OSELM	0.9994	0.9673	0.9995	0.9912	0.0206	0.0124	0.0163	0.0259	0.0118	0.0104	0.0115	0.0241
OSELM	0.9974	0.9721	0.9953	0.9862	0.0449	0.1272	0.0515	0.1444	0.0331	0.1268	0.0328	0.1439
ELM	0.9901	0.9881	0.9991	0.9802	0.0931	0.2190	0.0994	0.1837	0.0666	0.2186	0.0769	0.1822

下载: 导出CSV

| 显示表格

由图8可知,在对各监测点随机项位移预测时,本文模型仍然表现出了较高精度,以ZK38+015断面D1-5监测点表现最为明显,传统ELM预测时,偏差较大,而经GWO优化后的OSELM预测效果明显更优,模型泛化能力也较好,算法整体性能较稳定。

再分析表4,5,同样以ZK38+020断面D1-5监测点测试集预测为例,R达到了0.9731,R_MSE为0.0114 mm,M_AE为0.0093 mm,模型预测精度同样都远优于传统模型,表明所提模型应用于随机项位移预测中的适用性。另外,值得注意的是,在各监测点训练集的拟合中,传统ELM在拟合前期效果较好,但后期性能下降较快,出现了一定程度的“过拟合”现象,使得模型泛化能力变差,出现较大预测误差,这主要是由于传统的ELM为“静态”模型,它一次性对全部训练数据进行学习并拟合,而动态预测模型则是将训练数据分为多个子集分别进行学习并拟合,因此不易出现过拟合现象,这也是本文模型在地表沉降预测中的优势所在。

（3）累积位移预测及分析

利用前述时序分析原理,将趋势项及随机项位移进行叠加,便得到了各监测点累积位移预测结果,如图9所示,各监测点预测精度及误差见表6,7。

图 9 累积位移预测结果

Figure 9. Predicted results of cumulative displacement

下载: 全尺寸图片幻灯片

表 6 ZK38+020累积位移预测精度及误差

Table 6. Prediction accuracies and errors of cumulative displacement of ZK38+020

模型	R				R_MSE/ mm				M_AE/ mm
	D1-4		D1-5		D1-4		D1-5		D1-4		D1-5
	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集
GWO-OSELM	0.9999	0.9994	0.9999	0.9997	0.0230	0.0167	0.0340	0.0093	0.0203	0.0136	0.0281	0.0071
OSELM	0.9989	0.9998	0.9991	0.9997	0.0976	0.3085	0.1125	0.3201	0.0762	0.3000	0.0919	0.3143
ELM	0.9991	0.9986	0.9990	0.9993	0.0976	0.4865	0.1255	0.6100	0.0697	0.4700	0.0962	0.5979

下载: 导出CSV

| 显示表格

表 7 ZK38+015累积位移预测精度及误差

Table 7. Prediction accuracies and errors of cumulative displacement of ZK38+015

模型	R				R_MSE/ mm				M_AE/ mm
	D1-4		D1-5		D1-4		D1-5		D1-4		D1-5
	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集	训练集	测试集
GWO-OSELM	0.9999	0.9975	0.9998	0.9996	0.0274	0.0567	0.0302	0.0246	0.0200	0.0443	0.0241	0.0207
OSELM	0.9989	0.9984	0.9983	0.9986	0.0895	0.3679	0.0989	0.4063	0.0766	0.3643	0.0806	0.3971
ELM	0.9962	0.9951	0.9960	0.9980	0.1736	0.7358	0.1629	0.7763	0.1338	0.7250	0.1125	0.7550

下载: 导出CSV

| 显示表格

由图9对比可得,GWO-OSELM模型的预测结果与地表沉降实测值基本一致,而传统ELM模型预测结果偏差较大,且预测误差表现出了随着监测时间的延长而逐渐增大的特点,这也是传统“静态”模型在沉降长期预测中所表现出的不足。而OSELM模型在两断面D1-5监测点的预测中也出现了较大误差,其原因主要是两断面监测点随机项位移的预测精度均不是特别高,这也说明了随机项位移对累积位移预测精度有一定影响。

再结合表6,7进一步分析,同样以ZK38+020断面D1-5监测点测试集预测为例,R达到了0.9997,R_MSE为0.0093 mm,M_AE为0.0071 mm,模型预测精度也都远优于传统模型,具有较高精度及泛化能力,算法整体稳定性较好。因此,本文所提模型更适合应用于洞口段地表沉降长期预测中,具有较好的推广应用价值。

3. 讨论

3.1 激励函数的选择

由于神经网络的上一层输出需通过激励函数转换之后才能输入到下一层,因此,激励函数的选择十分重要^[19]。为了比较不同激励函数对模型预测性能的影响,选取常用的“径向基函数（RBF）”、“Sigmoid函数”“Sine函数”,并以ZK38+020断面D1-5监测点的实测值为例进行分析。表8为选择不同激励函数时的预测精度和误差,图10为不同激励函数时的预测结果。

表 8 模型预测精度及误差

Table 8. Prediction accuracies and errors of model

类别	R		R_MSE/ mm		M_AE/ mm
类别	训练集	测试集	训练集	测试集	训练集	测试集
RBF	0.9996	0.9813	0.0746	0.5205	0.0625	0.4852
Sigmoid	0.9999	0.9997	0.0340	0.0093	0.0281	0.0071
Sine	0.9999	0.9991	0.0314	0.2221	0.0266	0.2077

下载: 导出CSV

| 显示表格

图 10 不同激励函数时的预测结果

Figure 10. Predicted results of different activation functions

下载: 全尺寸图片幻灯片

由表8可知,对于训练集拟合结果而言,以Sine函数稍优一些,但测试集的预测精度则以Sigmoid函数最优。再结合图10可得,当模型参数设置相同时,采用RBF函数,测试集预测结果偏差最大,而采用Sigmoid函数,预测结果与实测值基本一致,采用Sine函数,预测结果介于两者之间。综上,本文模型预测时采用Sigmoid函数整体要更优,模型泛化能力较好。

值得注意的是,对于新建工程,往往只有前期监测数据,然而对于同一断面而言,其工程地质特征及条件往往是固定的,前期所获监测数据一定程度上也可以反映出其未来变化趋势,因此实际应用时,激励函数可通过前期监测数据的训练和测试进行选择。

3.2 随机项位移优势影响因素分析

事实上,对于随机项位移的预测,若以影响因素作为模型输入变量,则更利于实际应用,但鉴于目前监测手段有限,很多随机性因素无法获得。因此,本文仍以分解所得随机项位移进行模型训练和预测,并根据现有监测手段所获数据,以ZK38+020断面D1-5监测点为例,对其部分影响因素（单日降雨量、双日降雨量、拱顶下沉速率及地表沉降速率,如图11所示）进行分析。

图 11 随机项位移影响因素

Figure 11. Impact factors of random term displacement

下载: 全尺寸图片幻灯片

本文借鉴生物化学领域相关性分析新方法（GRaMM）^[20],尝试将其引入随机项位移优势影响因素分析中,具体流程如图12所示。它是由贾伟教授团队于2019年首次提出的一种相关性分析方法,整合了线性回归、最大信息系数（MIC）等多种方法,能够较好的处理线性和非线性关系,最初应用于代谢组与微生物组的相关性分析中。为更好的验证其合理性,与文献[18]中的MIC分析结果进行对比,见表9,当某一影响因素与随机项位移相关性越强时,数值越接近0,此时可认为该影响因素是相对主要的。

图 12 GRaMM相关性计算流程

Figure 12. Calculation flowchart of correlation with GRaMM

下载: 全尺寸图片幻灯片

表 9 相关性对比结果

Table 9. Comparitive results of correlation

方法	单日降雨量	双日降雨量	拱顶下沉速率	地表沉降速率
GRaMM	0.4968	0.4519	0.7114	0.7664
MIC	0.4716	0.4584	0.8642	0.9600

下载: 导出CSV

| 显示表格

分析表9可知,在现有监测手段所获数据中,地表沉降速率影响程度最大,其次是拱顶下沉速率,而降雨量影响稍小一些,原因是该段隧道施工期间采取了较为得当的地表防排水措施,并且地层沉降较稳定。另外,GRaMM的分析结果与文献[18]中MIC的分析结果基本一致,且前者结果整体偏小一些,主要是由于它可以有效的对数据中存在的异常值进行处理。因此,利用GRaMM分析随机项位移影响因素具备一定的合理性和可行性。随着监测手段日趋完善,将会有更多的影响因素被采集,然后利用该方法分析出优势影响因素,并以此作为模型输入变量,即可更精确地预测随机项位移。

4. 结论

（1）采用三次样条函数插值法对地表沉降实测值进行等距处理,并结合时序分析理论和VMD分解方法对其进行分解,可以有效降低原始数据复杂性,提高预测准确性,并且可以反映出位移分量自身尺度特征。

（2）将GWO-OSELM模型应用于重庆兴隆隧道洞口段ZK38+020和ZK38+015断面D1-4,D1-5监测点实测地表沉降预测中,对比分析结果表明了该模型在趋势项、随机项及累积位移预测中均要优于OSELM模型及传统ELM模型,模型预测精度高且误差小,具有很好的泛化能力,验证了该模型的合理性及工程实用性。

（3）深入探讨了激励函数的选择对模型预测性能的影响,发现Sigmoid激励函数更适合本文模型,此时模型预测性能最佳,且具有很好的泛化能力。

（4）通过引入GRaMM相关性分析新方法,对地表沉降随机项部分影响因素进行定量分析,发现地表沉降速率和拱顶下沉速率对随机项位移影响程度要高于降雨量,可为后续更精确的预测随机项位移提供借鉴。

图 1 原状黄土与压实黄土的SEM图像

Figure 1. SEM photographs of intact loess and compacted loess

下载: 全尺寸图片幻灯片

图 2 具有双层孔隙结构的土的理想孔径分布曲线

Figure 2. Idealized pore-size distribution curve of soils with double-level structures

下载: 全尺寸图片幻灯片

图 3 传统方法测得的孔隙比与MIP测得的孔隙比

Figure 3. Comparison of void ratios measured by traditional method and MIP

下载: 全尺寸图片幻灯片

图 4 不同固结压力下的PSDD

Figure 4. PSDD under different consolidation pressures

下载: 全尺寸图片幻灯片

图 5 三轴试验结果

Figure 5. Results of triaxial tests

下载: 全尺寸图片幻灯片

图 6 不同轴向应变下的PSDD

Figure 6. PSDD under different axial strains

下载: 全尺寸图片幻灯片

图 7 分形维数与孔隙比之间的关系

Figure 7. Relationship between fractal dimension and e_in

下载: 全尺寸图片幻灯片

图 8 分形维数与轴向应变的关系

Figure 8. Relationship between fractal dimension and axial strain

下载: 全尺寸图片幻灯片

图 9 固结和剪切过程中土体微结构变化概念图

Figure 9. Conceptual model for change of soil microstructure during consolidation and shearing

下载: 全尺寸图片幻灯片

表 1 土样基本物理指标

Table 1 Physical properties of soil specimens

含水率/%	天然密度/(g·cm^-3)	孔隙比	土粒相对质量密度	液限/%	塑限/%
15.0~20.0	1.44~1.46	0.94~1.15	2.70	34.2	18.6

下载: 导出CSV

表 2 基于Menger模型方法得出的分形维数表

Table 2 Fractal dimensions based on Menger model method

试样	D₁	D₂	D₃	D₄
Remolded	3.3990	4.9880	2.6382	2.1776
100-0%	2.8176	7.1840	2.7800	2.0711
100-5%	2.4976	7.4127	2.0212	—
100-10%	2.2121	8.4766	2.7554	2.0197
100-15%	2.5802	10.6195	2.8180	1.2609
100-20%	2.5841	7.1734	2.6205	0.7869
200-0%	2.4513	6.8217	2.4318	0.1943
200-5%	2.6717	8.9774	2.7851	1.2439
200-10%	2.6678	7.3717	2.8752	1.9276
200-15%	2.4132	10.7486	2.6901	2.2662
200-20%	2.1522	7.1398	2.8394	2.2585
300-5%	2.5734	9.1076	2.8492	1.9301
300-10%	2.4648	8.9810	2.8211	1.3640
300-15%	2.5172	10.4048	2.8172	1.4712
300-20%	2.2914	6.2945	2.8790	2.4006
400-0%	2.6769	7.2853	2.8671	—
400-5%	2.4841	7.5416	2.8461	2.5499
400-10%	2.3108	6.1133	2.8823	2.7861
400-15%	2.7839	9.1771	2.8287	1.8220
400-20%	2.4452	8.7693	2.8017	2.0617
注：表中的横线表示无数据,即试样只有3个分形区间。

下载: 导出CSV

参考文献(21)

[1]	JUANG C H, DIJKSTRA T, WASOWSKI J, et al. Loess geohazards research in China: Advances and challenges for mega engineering projects[J]. Engineering Geology, 2019, 251: 1-10. doi: 10.1016/j.enggeo.2019.01.019
[2]	LI X, LI L C, SONG Y X, et al. Characterization of the mechanisms underlying loess collapsibility for land-creation project in Shaanxi Province, China—a study from a micro perspective[J]. Engineering Geology, 2019, 249: 77-88. doi: 10.1016/j.enggeo.2018.12.024
[3]	曹杰, 张继文, 郑建国, 等. 黄土地区平山造地岩土工程设计方法浅析[J]. 岩土工程学报, 2019, 41(S1): 109-112. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2019S1029.htm CAO Jie, ZHANG Ji-wen, ZHENG Jian-guo, et al. Design of high fill reclamation projects in loess areas[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(S1): 109-112. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC2019S1029.htm
[4]	谢定义, 齐吉琳. 土结构性及其定量化参数研究的新途径[J]. 岩土工程学报, 1999, 21(6): 651-656. doi: 10.3321/j.issn:1000-4548.1999.06.003 XIE Dingyi, QI Ji-lin. Soil structure characteristics and new approach in research on its quantitative parameter[J]. Chinese Journal of Geotechnical Engineering, 1999, 21(6): 651-656. (in Chinese) doi: 10.3321/j.issn:1000-4548.1999.06.003
[5]	胡再强, 沈珠江, 谢定义. 结构性黄土的变形特性[J]. 岩石力学与工程学报, 2004, 23(24): 4142-4146. doi: 10.3321/j.issn:1000-6915.2004.24.008 HU Zai-qiang, SHEN Zhu-jiang, XIE Ding-yi. Deformation properties of structural loess[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(24): 4142-4146. (in Chinese) doi: 10.3321/j.issn:1000-6915.2004.24.008
[6]	方祥位, 申春妮, 汪龙, 等. Q₂黄土浸水前后微观结构变化研究[J]. 岩土力学, 2013, 34(5): 1319-1324. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201305014.htm FANG Xiang-wei, SHEN Chun-ni, WANG Long, et al. Collapsible deformation properties of Q₂ loess[J]. Rock and Soil Mechanics, 2013, 34(5): 1319-1324. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201305014.htm
[7]	LI P, XIE W L, PAK R Y S, et al. Microstructural evolution of loess soils from the Loess Plateau of China[J]. CATENA, 2019, 173: 276-288. doi: 10.1016/j.catena.2018.10.006
[8]	JIANG M J, SHEN Z J, ADACHI T, et al. Microscopic analysis on artificially-prepared structured collapsible loess[J]. Chinese Journal of Geotechnical Engineering, 1999, 21(4): 486-491. doi: 10.3321/j.issn:1000-4548.1999.04.019
[9]	NG C W W, SADEGHI H, HOSSEN S B, et al. Water retentionandvolumetriccharacteristicsofintactand recompacted loess[J]. Canadian Geotechnical Journal, 2016, 53(8): 1258-1269. doi: 10.1139/cgj-2015-0364
[10]	SHAO X X, ZHANG H Y, TAN Y. Collapse behavior and microstructural alteration of remolded loess under graded wetting tests[J]. Engineering Geology, 2018, 233: 11-22. doi: 10.1016/j.enggeo.2017.11.025
[11]	蒲毅彬, 陈万业, 廖全荣. 陇东黄土湿陷过程的CT结构变化研究[J]. 岩土工程学报, 2000, 22(2): 49-54. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200001008.htm PU Yi-bin, CHEN Wan-ye, LIAO Quan-rong. Research on CT structure changing for damping process of loess in Longdong[J]. Chinese Journal of Geotechnical Engineering, 2000, 22(2): 49-54. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC200001008.htm
[12]	雷胜友, 唐文栋. 黄土在受力和湿陷过程中微结构变化的CT扫描分析[J]. 岩石力学与工程学报, 2004, 23(24): 4166-4169. doi: 10.3321/j.issn:1000-6915.2004.24.013 LEI Sheng-you, TANG Wen-dong. Analysis of microstructure change for loess in the process of loading collapse with CT scanning[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(24): 4166-4169. (in Chinese) doi: 10.3321/j.issn:1000-6915.2004.24.013
[13]	朱元青, 陈正汉. 原状Q₃黄土在加载和湿陷过程中细观结构动态演化的CT-三轴试验研究[J]. 岩土工程学报, 2009, 31(8): 1219-1228. doi: 10.3321/j.issn:1000-4548.2009.08.011 ZHU Yuan-qing, CHEN Zhenghan. Experimental study on dynamic evolution of mesostructure of intact Q₃loess during loading and collapse using CT andtriaxialapparatus[J]. ChineseJournalof Geotechnical Engineering, 2009, 31(8): 1219-1228. (in Chinese) doi: 10.3321/j.issn:1000-4548.2009.08.011
[14]	CASINI F, VAUNAT J, ROMERO E, et al. Consequences on water retention properties of double-porosity features in a compacted silt[J]. Acta Geotechnica, 2012, 7(2): 139-150. doi: 10.1007/s11440-012-0159-6
[15]	JIANG M J, ZHANG F G, HU H J, et al. Structural characterization of natural loess and remolded loess under triaxial tests[J]. Engineering Geology, 2014, 181: 249-260. doi: 10.1016/j.enggeo.2014.07.021
[16]	CHEN Z F, CHEN H E, LI J F, et al. Study on the changing rules of silty clay's pore structure under freeze-thaw cycles[J]. Advances in Civil Engineering, 2019, 2019(1): 1-11.
[17]	DELAGE P, LEFEBVRE G. Study of the structure of a sensitiveChamplainclayandofitsevolutionduring consolidation[J]. Canadian Geotechnical Journal, 1984, 21(1): 21-35. doi: 10.1139/t84-003
[18]	GRIFFITHS F J, JOSHI R C. Change in pore size distribution due to secondary consolidation of clays[J]. Géotechnique, 1989, 39(1): 159-167. doi: 10.1680/geot.1989.39.1.159
[19]	胡冉, 陈益峰, 周创兵. 基于孔隙分布的变形土土水特征曲线模型[J]. 岩土工程学报, 2013, 35(8): 1451-1462. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201308013.htm HU Ran, CHEN Yi-feng, ZHOU Chuang-bing. A water retention curve model for deformable soils based on pore size distribution[J]. Chinese Journal of Geotechnical Engineering, 2013, 35(8): 1451-1462. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201308013.htm
[20]	FRIESEN W I, MIKULA R J. Fractal dimensions of coal particles[J]. Journal of Colloid and Science, 1987, 120(1): 263-271.
[21]	SUN HQ, MAŠÍND , NAJSERJ , etal. Fractal characteristics of pore structure of compacted bentonite studied by ESEM and MIP methods[J]. Acta Geotechnica, 2020, 15(6): 1655-1671. doi: 10.1007/s11440-019-00857-z

施引文献(32)

期刊类型引用(18)

1.	柴永进，渠根启，潘童. 考虑降雨作用下公路高边坡变形预测的ELM方法及工程应用. 科技创新与应用. 2024(01): 98-102 . 百度学术
2.	崔靖奇，吴顺川，程海勇，王涛，姜关照，浦仕江，任子健. 滇中引水软岩隧洞围岩位移时序预测. 清华大学学报(自然科学版). 2024(07): 1215-1225 . 百度学术
3.	马恩临，赖金星，王立新，汪珂，雷升祥，李储军，邱军领. 基于控制区间牵引算法的地下施工变形预测. 岩土力学. 2023(02): 577-594 . 百度学术
4.	哈吉章，杨良义，肖旺槟，李晶生，彭乙芹. 基于组合预测模型的地铁车站地表沉降研究. 广东土木与建筑. 2023(02): 49-53 . 百度学术
5.	王述红，董福瑞. 基于变形预测和参数反演的山岭隧道围岩稳定性分析. 岩土工程学报. 2023(05): 1024-1035 . 本站查看
6.	成睿，李素敏，韩追，毛嘉骐，李彦臣. 时序InSAR与GWO-VMD相结合的地表沉降预测. 贵州大学学报(自然科学版). 2023(03): 78-85 . 百度学术
7.	李永靖，王松，刘维青，文成章，高航. 某隧道穿越富水断层破碎带围岩大变形成因分析及稳定性预测. 矿业研究与开发. 2023(05): 84-90 . 百度学术
8.	侯明华，袁颖，杨丛铭，李云鹏，黄虎城. 基于麻雀搜索算法优化Elman残差自校正地面沉降预测模型. 科学技术与工程. 2023(13): 5470-5480 . 百度学术
9.	付宏渊，刘忠伟，邱祥，罗震宇. 红黏土路堤填挖交界段工后沉降变形研究. 交通科学与工程. 2023(02): 1-7+15 . 百度学术
10.	陈城，史培新，王占生，贾鹏蛟. 基于融合多注意力机制的深度学习的盾构荷载预测方法. 东北大学学报(自然科学版). 2023(11): 1631-1637+1646 . 百度学术
11.	尹宏，王述红，董卓然，侯钦宽. 引入因子分析的结构面粗糙度RBF复合参数模型. 岩土工程学报. 2022(04): 721-730 . 本站查看
12.	林广东，何军，申小军，徐龙飞，裴莉莉，余婷. 基于随机森林的隧道建成初期累计沉降量预测. 计算技术与自动化. 2022(01): 160-163 . 百度学术
13.	宋光浩. 基于数据插值预测路基沉降及误差分析. 地理空间信息. 2022(04): 175-177 . 百度学术
14.	叶勇超，闫超德，罗先学，张瑞峰，袁观杰. 时序InSAR郑州地铁沿线地面沉降分析. 遥感学报. 2022(07): 1342-1353 . 百度学术
15.	胡仕明，杨伟红，李涛，李昕堃. 公路隧道洞口顺层边坡变形规律分析. 岩土工程技术. 2022(06): 477-482 . 百度学术
16.	王述红，董福瑞，朱宝强，刘欢，张泽. 山岭隧道围岩参数智能反演及稳定性分析. 应用基础与工程科学学报. 2021(05): 1171-1185 . 百度学术
17.	朱宝强，王述红，张泽，王鹏宇，董福瑞. 基于时间序列与DEGWO-SVR模型的隧道变形预测方法. 浙江大学学报(工学版). 2021(12): 2275-2285 . 百度学术
18.	柳明. 基于划分沉降区域的基坑沉降预测方法. 测绘标准化. 2021(04): 35-38 . 百度学术