Damage evolution mechanism of coal rock under long-term soaking
-
摘要: 富水状态下煤岩力学性能的持续劣化是导致煤矿水害的重要原因,为探究长期浸润作用下煤岩的渐进损伤特征,以神木地区弱黏煤为研究对象,首先以无压力水浸泡的方式模拟煤岩的长期浸润(15,30,60,120,240 d);然后对干燥及浸润试样进行渐进破坏试验,辅以声发射(AE)监测,定量分析煤岩的软化特性;最后利用扫描电镜(SEM)和核磁共振(NMR)研究不同浸润时长下煤岩孔隙结构的劣化规律。研究表明:①煤岩饱和度随浸润时长的增加而增加,增加速度先快后慢,最后达到稳定。②随饱和度的增加,煤岩力学性能劣化;饱和后,随浸润时长的持续增加,煤岩力学性能仍持续劣化,表现为抗压强度、弹性模量的减小和累计振铃数的增加;③干燥处理后,原状煤岩表面干燥收缩,浸水后吸水膨胀,随着浸润时间的增加,小孔隙逐渐增多,孔隙结构趋于复杂。研究结果可为浸水煤岩失稳破坏问题提供有益参考。Abstract: The deterioration of the mechanical properties of water-rich coal rock is a predisposition to many mine geologic hazards. To explore the progressive damage characteristics of coal rock under long-term soaking, the weakly caking coal of the Shenmu region is selected as the research object. Firstly, the long-term soaking tests on coal rock are carried out without pressures for 15, 30, 60, 120 and 240 days. Then, the uniaxial compression tests with acoustic emission (AE) monitoring are carried out on the dry and then soaked coal rock to analyze its water-softening characteristics. Finally, the evolution of pore structure of coal rock under different soaking durations is studied using the scanning electron microscope (SEM) and nuclear magnetic resonance (NMR). The results show that: ①The water saturation of coal rock increases with the duration of soaking, and the increase rate continuously decreases, and tends to a constant. ②In the process of saturation, the mechanical properties of coal rock deteriorate, and those of saturated coal rock still deteriorate with the increase of the soaking duration, characterized by the decrease of the compressive strength and elastic modulus, and an increase in the cumulative AE ringing numbers. ③A surface shrink may occur after the drying process, and during the soaking tests, the small pores of samples gradually increase, and the pore structure becomes more complex with the increase of the soaking duration. The research results can provide useful prevention for the instability and failure of coal rock under water-rich conditions.
-
Keywords:
- coal rock /
- soaking duration /
- mechanical property /
- multi-scale /
- damage mechanism
-
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为隧道开挖宽度。现有研究表明,地表沉降是受基础地质条件(如地形地貌和岩性等)和外部因素(如降雨、风荷载、测量条件等)等共同作用形成的[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}{∑k‖∂t[(μ(t)+iπn)∗δk(t)]e−iωkt‖22}, (2) s.t. K∑kδk=s(t), (3) 式中,K为分解所得的模态个数(正整数),
{δk} ,{ωk} 分别为分解后的第K个模态分量及对应中心频率,μ(t) 为狄拉克函数,∗ 为卷积运算符。(2)求解变分问题。引入Lagrange乘法算子,将上式转变为非约束变分问题,得到增广Lagrange表达式为
L({δk},{ωk},λ)=σ∑k‖∂t[(μ(t)+i/πt)∗δk(t)]e−iωkt‖22+‖s(t)−∑kδk(t)‖22+〈λ(t),s(t)−∑kδk(t)〉, (4) 式中,
σ 为二次惩罚因子。利用交替方向乘子法将上述问题转化为寻找增广Lagrange方程的“鞍点”问题,也即最优解,限于篇幅,具体过程见文献[18]。本文将其引入地表沉降分解中,得到不同特征尺度的位移分量,更能反映其自身信息。
1.2 地表沉降动态预测的GWO-OSELM模型
(1)在线贯序极限学习机模型
传统神经网络模型(如SVM、ELM等)多为“静态”的,预测时只能学习并使用当前时刻监测值,并且一次性学习全部训练数据,而地表沉降监测值往往是非线性动态变化的,并且彼此之间存在一定联系。鉴于此,将具有在线学习机制的OSELM模型[12]引入地表沉降预测中,它是对传统ELM模型的改进,能够逐批学习监测值并根据最新监测值不断更新网络,实现动态预测。
OSELM模型主要通过利用递推计算对ELM模型的输出层权值
ˆβ 进行更新,主要包括两部分。a)初始化阶段
给定激励函数g(x),隐含层节点数n,以及初始训练子集
ψ0={(xi,yi)|i=1, 2, ⋯, N0} ,并且N0≥n ,则初始隐含层输出权值可由下式求得,β0=P0HT0T0, (5) 式中,
P0=(HT0H0)−1 ,T0=[y1y2⋅⋅⋅yN0]T ,H0 为隐含层初始输出矩阵。b)在线学习阶段
令m=0,m为加入训练数据的批次,对m+1批数据,计算其隐层输出矩阵Hm+1,并按下式更新隐层输出权值:
Pm+1=Pm−PmhTm+1hm+1Pm1+hm+1PmhTm+1, (6) βm+1=βm+Pm+1HTm+1(TTm+1−Hm+1βm), (7) 式中,
hm=[G(a1,b1,xm)⋅⋅⋅G(an,bn,xm)] ,G(an, bn,xm)= gn(xm) ,an ,bn分别是第n个隐层节点的输入权值和阈值。(2)灰狼优化算法
GWO由澳大利亚学者Mirjalili等[13]在2014年首次提出,它模拟了自然界中灰狼种群的社会等级制度和群体狩猎行为,灰狼群体按社会等级制度高低依次分为α狼、β狼、γ狼和
ω 狼,其群体狩猎行为如下:Step 1狼群跟踪并接近猎物。按下式计算距离:
D=|C⋅Xp(u)−X(u)|, (8) C=2×r1, (9) 式中,D为灰狼与猎物间的距离,
u 为当前迭代次数,Xp(u )为第u 代猎物位置向量,X(u )为第u 代灰狼个体位置向量,C为系数向量,r1为[0,1]间均匀分布的随机向量。Step 2狼群追捕并包围猎物。按下式更新位置:
X(u+1)=Xp(u)−η⋅D, (10) η=2b×r2−b, (11) 式中,
η 为系数向量,r2 为b的随机向量,b的初始值为2,随最大迭代次数增加而线性递减至0。Step 3狼群攻击并捕杀猎物。依据狼群中
α 狼寻找猎物方向,按下式计算ω 狼与各等级狼的距离:Dl=|Cj⋅Xl(u)−X(u)|, (12) Xj=Xl(u)−ηj⋅Dl, (13) X(u+1)=X1+X2+X33。 (14) 式中 l=α,β,γ;j=1,2,3。
(3)GWO-OSELM模型动态预测流程
虽然OSELM模型存在诸多优点,但它同样是通过随机产生输入权值及隐层阈值来进行计算,存在诸多缺点,故引入GWO算法对其输入权值及隐层阈值进行优化,以提高模型预测精度,并改善模型泛化能力。如图1所示,基于前述原理建立了该模型,具体步骤如下:
a)判断洞口段地表沉降监测值是否等距,若是,则利用VMD将其分解为趋势项及随机项,若否,则先利用3次样条插值将其转化为等距后再进行分解。
b)为便于后续进行最优权值及阈值的搜寻及消除数据量纲的影响,使算法更容易收敛,分别对分解得到的时序数据进行归一化操作,使其位于[-1,1]。
c)初始化算法相关参数,主要包括狼群数量、最大迭代次数、搜索变量最大值及最小值、狼群个体初始位置及初始目标函数值。
d)以模型测试集预测的均方根误差作为个体适应度计算值,将适应度最优的个体依次作为
α 狼、β 狼和γ 狼,然后更新狼群位置,并更新α 狼、β 狼和γ 狼。e)判断是否达到最大迭代次数,若否,则继续执行步骤d),若是,则输出OSELM最优权值及阈值。
f)将结果用于模型训练和预测中,得位移分量预测结果,由时间序列原理叠加得到地表沉降总预测值。
1.3 模型预测结果评价指标
相关性系数R、均方根误差RMSE、平均绝对误差MAE 3项指标常用来衡量位移预测精度[17]。因此,本文采用它们来评价模型预测效果,
R=N∑j=1(dj−ˉd)(ˆdj−ˉˆd)√N∑j=1(dj−ˉd)2(ˆdj−ˉˆd)2, (15) RMSE=√1NN∑j=1(ˆdj−dj)2, (16) MAE=1NN∑j=1|ˆdj−dj|, (17) 式中,
dj 为沉降监测数据,ˆdj 为沉降预测数据,ˉd 为沉降监测值均值,ˉˆd 为沉降预测值均值,3项指标均为相对越小时预测结果越好。2. 工程实例分析
2.1 工程概况
在建的兴隆隧道位于重庆市渝北区木耳镇,隧道左线全长2553 m,右线全长2539 m。前期地勘结果表明,隧址区域地质较稳定,总体呈单斜构造,主要穿越侏罗系中统上沙溪庙组地层。隧道出洞口段地形较平缓,缓坡主要分布第四系残坡积层粉质黏土,下伏基岩主要为泥岩和砂岩(图2),岩体较破碎,洞口段山体覆盖层较薄,围岩自稳能力差。因此,施工过程中及时预测地表沉降变化,对于分析隧道进洞稳定性、保障施工安全及塌方灾害防治具有重要意义。
本文以隧道左线出口段ZK38+020和ZK38+015断面地表沉降实测值为研究对象,隧道采用三台阶法施工,初期支护为钢拱架+系统锚杆+27 cm厚喷射混凝土,监测断面平面位置见图3,沿隧道横断面每隔2~5 m布置监测点位置见图4,选择距隧道拱顶较近的D1-4,D1-5监测点前50 d沉降实测值进行分析。
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所示。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。
表 2 ZK38+020趋势项位移预测精度及误差Table 2. Prediction accuracies and errors of trend term displacement of ZK38+020模型 R RMSE/ mm MAE/ 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 表 3 ZK38+015趋势项位移预测精度及误差Table 3. Prediction accuracies and errors of trend term displacement of ZK38+015模型 R RMSE/ mm MAE/ 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 由图7中的结果对比可知,在测试集的预测中,具有动态预测能力的GWO-OSELM模型及OSELM模型的预测效果要明显优于传统的“静态”ELM模型,这主要是由于前者具有在线学习能力,能够逐批学习监测数据,并不断的更新网络,而引入GWO后的OSELM模型预测结果也要优于传统OSELM模型,原因主要是灰狼算法具有较强的寻优能力,能够得到更为准确的模型输入权值及阈值,相较于传统模型而言,优势更加明显。
再结合表2,3进一步分析,以ZK38+020断面D1-5监测点测试集预测为例,R达到了0.9998,RMSE为0.0051 mm,MAE为0.0044 mm,3项误差指标计算结果均以本文模型更优,这也进一步验证了模型用于趋势项位移预测的可靠性,能得到较为满意的预测结果。
(2)随机项位移预测及分析
与趋势项位移相比,随机项位移表现出了一定程度的波动及随机现象,并且具有一定的平稳性,因此,应用所提模型对随机项位移进行预测,更易体现出模型较强的预测能力。分别利用3种模型对两断面D1-4,D1-5监测点随机项位移进行预测,预测结果见图8,两断面各监测点预测精度和误差分别见表4,5。
表 4 ZK38+020随机项位移预测精度及误差Table 4. Prediction accuracies and errors of random term displacement of ZK38+020模型 R RMSE/ mm MAE/ 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 表 5 ZK38+015随机项位移预测精度及误差Table 5. Prediction accuracies and errors of random term displacement of ZK38+015模型 R RMSE/ mm MAE/ 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 由图8可知,在对各监测点随机项位移预测时,本文模型仍然表现出了较高精度,以ZK38+015断面D1-5监测点表现最为明显,传统ELM预测时,偏差较大,而经GWO优化后的OSELM预测效果明显更优,模型泛化能力也较好,算法整体性能较稳定。
再分析表4,5,同样以ZK38+020断面D1-5监测点测试集预测为例,R达到了0.9731,RMSE为0.0114 mm,MAE为0.0093 mm,模型预测精度同样都远优于传统模型,表明所提模型应用于随机项位移预测中的适用性。另外,值得注意的是,在各监测点训练集的拟合中,传统ELM在拟合前期效果较好,但后期性能下降较快,出现了一定程度的“过拟合”现象,使得模型泛化能力变差,出现较大预测误差,这主要是由于传统的ELM为“静态”模型,它一次性对全部训练数据进行学习并拟合,而动态预测模型则是将训练数据分为多个子集分别进行学习并拟合,因此不易出现过拟合现象,这也是本文模型在地表沉降预测中的优势所在。
(3)累积位移预测及分析
利用前述时序分析原理,将趋势项及随机项位移进行叠加,便得到了各监测点累积位移预测结果,如图9所示,各监测点预测精度及误差见表6,7。
表 6 ZK38+020累积位移预测精度及误差Table 6. Prediction accuracies and errors of cumulative displacement of ZK38+020模型 R RMSE/ mm MAE/ 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 表 7 ZK38+015累积位移预测精度及误差Table 7. Prediction accuracies and errors of cumulative displacement of ZK38+015模型 R RMSE/ mm MAE/ 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 由图9对比可得,GWO-OSELM模型的预测结果与地表沉降实测值基本一致,而传统ELM模型预测结果偏差较大,且预测误差表现出了随着监测时间的延长而逐渐增大的特点,这也是传统“静态”模型在沉降长期预测中所表现出的不足。而OSELM模型在两断面D1-5监测点的预测中也出现了较大误差,其原因主要是两断面监测点随机项位移的预测精度均不是特别高,这也说明了随机项位移对累积位移预测精度有一定影响。
再结合表6,7进一步分析,同样以ZK38+020断面D1-5监测点测试集预测为例,R达到了0.9997,RMSE为0.0093 mm,MAE为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 RMSE/ mm MAE/ 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 由表8可知,对于训练集拟合结果而言,以Sine函数稍优一些,但测试集的预测精度则以Sigmoid函数最优。再结合图10可得,当模型参数设置相同时,采用RBF函数,测试集预测结果偏差最大,而采用Sigmoid函数,预测结果与实测值基本一致,采用Sine函数,预测结果介于两者之间。综上,本文模型预测时采用Sigmoid函数整体要更优,模型泛化能力较好。
值得注意的是,对于新建工程,往往只有前期监测数据,然而对于同一断面而言,其工程地质特征及条件往往是固定的,前期所获监测数据一定程度上也可以反映出其未来变化趋势,因此实际应用时,激励函数可通过前期监测数据的训练和测试进行选择。
3.2 随机项位移优势影响因素分析
事实上,对于随机项位移的预测,若以影响因素作为模型输入变量,则更利于实际应用,但鉴于目前监测手段有限,很多随机性因素无法获得。因此,本文仍以分解所得随机项位移进行模型训练和预测,并根据现有监测手段所获数据,以ZK38+020断面D1-5监测点为例,对其部分影响因素(单日降雨量、双日降雨量、拱顶下沉速率及地表沉降速率,如图11所示)进行分析。
本文借鉴生物化学领域相关性分析新方法(GRaMM)[20],尝试将其引入随机项位移优势影响因素分析中,具体流程如图12所示。它是由贾伟教授团队于2019年首次提出的一种相关性分析方法,整合了线性回归、最大信息系数(MIC)等多种方法,能够较好的处理线性和非线性关系,最初应用于代谢组与微生物组的相关性分析中。为更好的验证其合理性,与文献[18]中的MIC分析结果进行对比,见表9,当某一影响因素与随机项位移相关性越强时,数值越接近0,此时可认为该影响因素是相对主要的。
表 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 分析表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 煤岩的工业分析和元素分析
Table 1 Proximate analysis and ultimate analysis of coal rock
样品 工业分析/% 元素分析/(%daf) Mad Aad Vdaf C H O N S C/H H/O 煤岩 5.65 5.23 32.65 78.21 4.77 15.32 1.65 0.43 1.37 0.226 -
[1] 中华人民共和国国家统计局. 中国统计年鉴[M]. 北京: 中国统计出版社, 2022. National Bureau of Statistics of the People's Republic of China. China Statistical Yearbook[M]. Beijing: China Statistics Press, 2022. (in Chinese)
[2] ZHANG Anlin, XIE Heping, ZHANG Ru, et al. Mechanical properties and energy characteristics of coal at different depths under cyclic triaxial loading and unloading[J]. International Journal of Rock Mechanics and Mining Sciences, 2023, 161(1): 105271.
[3] 张玉军. 控水采煤技术原理、关键技术及在砂岩含水层下综放开采实践[J]. 煤炭学报, 2020, 45(10): 3380-3388. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB202010002.htm ZHANG Yujun. Principle and key technologies of controlled water mining and practice of fully-mechanized mining under soft sandstone aquifer[J]. Journal of China Coal Society, 2020, 45(10): 3380-3388. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB202010002.htm
[4] 武强, 申建军, 王洋. "煤-水"双资源型矿井开采技术方法与工程应用[J]. 煤炭学报, 2017, 42(1): 8-16. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201701002.htm WU Qiang, SHEN Jianjun, WANG Yang. Mining techniques and engineering application for "Coal-Water" dual-resources mine[J]. Journal of China Coal Society, 2017, 42(1): 8-16. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201701002.htm
[5] CHEN Tian, YAO Qiangling, WEI Fei. Effects of water intrusion and loading rate on mechanical properties of and crack propagation in coal–rock combinations [J]. Journal of Central South University, 2017, 24(2): 183-191.
[6] YAO Qiangling, CHEN Tian, XIA Ze, et al. Mechanisms of failure in coal samples from underground water reservoir[J]. Engineering Geology, 2020, 267: 105494. doi: 10.1016/j.enggeo.2020.105494
[7] ZHONG Chunlin, ZHANG Zhenyu, RANJITH P G, et al. The role of pore water plays in coal under uniaxial cyclic loading[J]. Engineering Geology, 2019, 257: 105125. doi: 10.1016/j.enggeo.2019.05.002
[8] WEN Guangcai, YANG Shuo, LIU Yanbao, et al. Influence of water soaking on swelling and microcharacteristics of coal[J]. Energy Science and Engineering, 2020, 8, 50-60. doi: 10.1002/ese3.508
[9] YANG Shuo, WEN Guangcai, YAN Fazhi, et al. Swelling characteristics and permeability evolution of anthracite coal containing expansive clay under water-saturated conditions [J]. Fuel, 2020, 279: 118501. doi: 10.1016/j.fuel.2020.118501
[10] 刘驰, 刘晓丽, 张东, 等. 软岩软化的水岩界面动力学模型及其演化规律[J]. 岩土工程学报, 2022, 44(12): 2280-2289. doi: 10.11779/CJGE202212015 LIU Chi, LIU Xiaoli, ZHANG Dong, et al. Dynamic model for water-rock interface of softening of soft rock and its evolution law[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(12): 2280-2289. (in Chinese) doi: 10.11779/CJGE202212015
[11] MENG Fanzhen, ZHOU Hui, ZHANG Chuanqing, et al. Evaluation methodology of brittleness of rock based on post-peak stress-strain curves[J]. Rock Mechanics and Rock Engineering, 2015, 48(5): 1787-1805. doi: 10.1007/s00603-014-0694-6
[12] VISHAL V, RANJITH P G, SINGH T N. An experimental investigation on behaviour of coal under fluid saturation, using acoustic emission[J]. Journal of Natural Gas Science and Engineering, 2015, 22: 428-436. doi: 10.1016/j.jngse.2014.12.020
[13] 夏冬, 杨天鸿, 王培涛, 等. 干燥及饱和岩石循环加卸载过程中声发射特征试验研究[J]. 煤炭学报, 2014, 39(7): 1243-1247. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201407008.htm XIA Dong, YANG Tianhong, WANG Peitao, et al. Experimental study of acoustic emission characteristics of dry and saturated rocks during cyclic loading and unloading process[J]. Journal of China Coal Society, 2014, 39(7): 1243-1247. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB201407008.htm
[14] 夏冬, 杨天鸿, 徐涛, 等. 浸水时间对饱水岩石损伤破坏过程中声发射特征影响的试验[J]. 煤炭学报, 2015, 40(增刊2): 337-345. https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB2015S2005.htm XIA Dong, YANG Tianhong, XU Tao, et al. Experimental study on AE properties during the damage process of water saturated rock specimens based on time effect[J]. Journal of China Coal Society, 2015, 40(S2): 337-345. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTXB2015S2005.htm
[15] ZHENG Chao, MA Dongmin, CHEN Yue, et al. Pore structure of different macroscopically distinguished components within low-rank coals and its methane desorption characteristics[J]. Fuel, 2021, 293: 120465. doi: 10.1016/j.fuel.2021.120465
[16] 张慧, 李小彦. 扫描电子显微镜在煤岩学上的应用[J]. 电子显微学报, 2004, 23(4): 467. doi: 10.3969/j.issn.1000-6281.2004.04.169 ZHANG Hui, LI Xiaoyan. Application of scanning electron microscope in coal petrology[J]. Journal of Chinese Electron Microscopy Society, 2004, 23(4): 467. (in Chinese) doi: 10.3969/j.issn.1000-6281.2004.04.169
[17] 张慧, 吴静, 袁立颖, 等. 煤中气孔的发育特征与影响因素浅析[J]. 煤田地质与勘探, 2019, 47(1): 78-85, 91. doi: 10.3969/j.issn.1001-1986.2019.01.011 ZHANG Hui, WU Jing, YUAN Liying, et al. Analysis on the development characteristics and influencing factors of gas pores in coal[J]. Coal Geology & Exploration, 2019, 47(1): 78-85, 91. (in Chinese) doi: 10.3969/j.issn.1001-1986.2019.01.011
[18] XU Yong, CHEN Xuexi, ZHAO Wei, et al. Effect of water intrusion on the characteristics of surface morphology and pore fracture spaces in argillaceous meagre coal[J]. Journal of Natural Gas Science and Engineering, 2020, 81: 103404. doi: 10.1016/j.jngse.2020.103404
[19] 张村, 贾胜, 王方田, 等. 水岩作用下煤体孔裂隙结构演化特征及驱动机制实验研究[J]. 应用基础与工程科学学报, 2023, 31(1): 185-196. https://www.cnki.com.cn/Article/CJFDTOTAL-YJGX202301014.htm ZHANG Cun, JIA Sheng, WANG Fangtian, et al. An experimental research on the pore and fracture evolution characteristics and its driving mechanism for coal samples under water-rock interaction[J]. Journal of Basic Science and Engineering, 2023, 31(1): 185-196. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YJGX202301014.htm
[20] 常会珍, 李贵山, 张蒙, 等. 松软低渗区不同煤体结构煤的孔裂隙发育特征[J]. 煤炭技术, 2023, 42(1): 133-137. https://www.cnki.com.cn/Article/CJFDTOTAL-MTJS202301027.htm CHANG Huizhen, LI Guishan, ZHANG Meng, et al. Pore and fissure development characteristics of coal with different coal structure in soft and low permeability area[J]. Coal Technology, 2023, 42(1): 133-137. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-MTJS202301027.htm
[21] XU Hao, TANG Dazhen, CHEN Yanpeng, et al. Effective porosity in lignite using kerosene with low-field nuclear magnetic resonance[J]. Fuel, 2018, 213: 158-163. doi: 10.1016/j.fuel.2017.10.112
[22] LIU Zhen, YANG He, WANG Wenyu, et al. Experimental Study on the pore structure fractals and seepage characteristics of a coal sample around a borehole in coal seam water infusion[J]. Transport in Porous Media, 2018, 125(2): 1-21.
[23] JIA Hailiang, DING Shun, WANG Yan, et al. An NMR-based investigation of pore water freezing process in sandstone [J]. Cold Regions Science and Technology, 2019, 168, 102893. doi: 10.1016/j.coldregions.2019.102893
[24] XIE Jingna, NI Guanhua, HAO Jian, et al. The pore characteristics of high temperature coal mine and the impact of invasive water on its development[J] Arabian Journal of Geosciences, 2019, 12(21): 657. doi: 10.1007/s12517-019-4856-4
[25] 王龙飞, 蒋仲安, 陈举师, 等. 低压注水对煤体孔隙特征及渗透率的影响[J]. 中国安全生产科学技术, 2018, 14(6): 108-113. https://www.cnki.com.cn/Article/CJFDTOTAL-LDBK201806017.htm WANG Longfei, JIANG Zhongan, CHEN Jushi, et al. Influence of low pressure water injection on pore characteristics and permeability of coal. Journal of Safety Science and Technology, 2018, 14(6): 108-113. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-LDBK201806017.htm
[26] 荣腾龙, 刘克柳, 周宏伟, 等. 采动应力下深部煤体渗透率演化规律研究[J]. 岩土工程学报, 2022, 44(6): 1106-1114. doi: 10.11779/CJGE202206015 RONG Tenglong, LIU Keliu, ZHOU Hongwei, et al. Permeability evolution of deep coal under mining stress[J]. Chinese Journal of Geotechnical Engineering, 2022, 44(6): 1106-1114. (in Chinese) doi: 10.11779/CJGE202206015
[27] 冯国瑞, 文晓泽, 郭军, 等. 含水率对煤样声发射特征和碎块分布特征影响的试验研究[J]. 中南大学学报(自然科学版), 2021, 52(8): 2910-2918. https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD202108035.htm FENG Gourui, WEN Xiaoze, GUO Jun, et al. Study on influence of moisture content on coal sample AE properties and fragment distribution characteristics[J]. Journal of Central South University (Science and Technology), 2021, 52(8): 2910-2918. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-ZNGD202108035.htm
[28] 张慧梅, 王云飞. 冻融红砂岩损伤演化多尺度分析[J]. 岩土力学, 2022, 43(8): 2103-2114. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202208006.htm ZHANG Huimei, WANG Yunfei. Multi-scale analysis of damage evolution of freezing-thawing red sandstones[J]. Rock and Soil Mechanics, 2022, 43(8): 2103-2114. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202208006.htm
[29] 楚亚培, 张东明, 王满, 等. 基于核磁共振技术和压汞法的液氮冻融煤体孔隙结构损伤演化规律试验研究[J]. 岩石力学与工程学报, 2022, 41(9): 1820-1831. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202209005.htm CHU Yapei, ZHANG Dongming, WANG Man, et al. Experiment study on influence of liquid nitrogen freeze-thaw on pore structure of coal based on nuclear magnetic resonance technology and mercury intrusion methods[J]. Chinese Journal of Rock Mechanics and Engineering, 2022, 41(9): 1820-1831. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX202209005.htm
-
期刊类型引用(0)
其他类型引用(5)
-
其他相关附件
-
本文英文pdf
20240609-G23-0229英文发布稿 点击下载(4305KB)
-