高 宁 朱宝训 李建刚
(1.河南城建学院测绘工程学院,河南 平顶山467036;2.矿山空间信息技术国家测绘地理信息局重点实验室,河南 焦作454003;3.河北联合大学迁安学院矿业与建工系,河北 唐山064400)
近年来,随着我国煤炭开采深度的增加,矿井软岩巷道的支护、维护等问题日益成为制约煤炭资源安全开采和向纵深发展的关键因素[1-5]。软岩[1-5,7-8](又称松软岩层)是指在工程外力作用下,能产生显著变形的工程岩体,此类岩体通常是不连续面、非均质的一类岩体,其力学变形机理复杂,受到开采方式、爆破、地质环境等众多因素的影响,具有大变形、大地压、难支护等特点。目前,基于现有的岩体力学理论来解决软岩体未来力学行为问题困难较大。为此,从灰色系统理论出发,避开软岩体力学的诸多复杂因素,将软岩巷道围岩变形时序看作是在一定范围内的灰色量[1-3,6],通过对灰色过程的数据开发,建立描述其动态变化规律的数学模型,并据此对软岩巷道围岩的变化进行定量预测。
(1)巷道围岩变形过程是一灰色过程。白色系统(信息完全)与灰色系统(信息不完全的系统)的区别在于,系统内各因素间是否具有确定性的函数关系。随着煤炭等矿产资源的开采,围岩应力随之发生改变,当其累积到一定量时,便出现岩体的塌落与破坏,从而影响矿山的安全生产。巷道的安全性、稳定性受到多种地质因素(煤柱宽度、地应力、围岩强度、支护形式等)影响,且各种地质因素均对岩体变形产生影响,影响关系复杂,难以建立明确的函数关系。
(2)巷道围岩变形数据符合灰度数据特征。采用各类监测设备对巷道围岩进行变形监测,所得观测数据序列实际上是巷道围岩变形系统多种影响因素总效应的综合反映,其中有众多难以描述的随机因素影响,此外由于观测过程受到各类误差的影响,使得观测数据具有某种程度上的不确定性特征。
式中,a 为灰发展系数(反映围岩变形原始观测数据序列和AGO 生成序列的发展态势);u 为灰色作用量(反映数据的内在变化)。
设待求参数序列为[a,u]T,由最小二乘法可得
设巷道围岩变形原始观测数据序列为x(0)(i)(i= 1,2,…,n),在建模过程中,可根据矿山开采进度等因素对围岩变形的影响,从x(0)(i)中选取不同的样本数据进行建模,通常不同的样本数据所构建的模型的预测结果差异较大。因此,可根据围岩变形发展态势,在预测建模中灵活选取建模样本,建立不同的GM(l,1)模型,在建模过程中所建立的不同GM(l,1)模型,称之为灰色模型群[1-3,6,9-11]。
设建模序列为x(0)= {x(0)(1),x(0)(2),…,x(0)(n)},将x(0)(n)定义为时间坐标轴的原点,则称t <n 为过去,t = n 为现在,t >n 时为未来。当t≤n 时,称为模型模拟值;t >n 时,称为模型的预测值。利用x(0)(i)预测时,为了提高预测结果的精度,应确保有足够高的模拟精度,特别是当t = n 时的模拟精度,根据软岩巷道围岩变形的实际情况,提出了基于3 类数据的灰色模型群:
(1)全数据预测模型。设观测数据序列为x(0)= {x(0)(1),x(0)(2),…,x(0)(n)}时,利用全部观测数据建立的GM(l,1)模型。
(2)新信息预测模型。设在建模过程中,获取最新观测数据x(0)(n +1),将x(0)(n +1)置入x(0),则利用包含最新数据建立的模型。
(3)新陈代谢预测模型。在建模过程中,置入最新观测数据x(0)(n + 1),同时删除最陈旧的历史数据x(0)(1),而建立x(0)= {x(0)(1),x(0)(2),…,x(0)(n),x(0)(n +1)}的模型。
某矿巷道位于松软破碎的变质岩中,对该巷道某变形严重段的顶底板1#点,两帮2#点分别进行围岩变形预测。建立了基于3 类数据的灰色模型群:①模型1(全数据模型),假定在巷道围岩监测过程中,仅获取1 ~7 期数据,建立GM(l,1)模型;②模型2(新信息模型),在建模过程中,又获取最新围岩变形监测量,即第8 期数据,在模型1 的建模数据样本中加入最新信息x(0)(8),即利用1 ~8 期观测数据建模;③模型3(新陈代谢模型),在模型1 的原始序列中删除最陈旧的信息x(0)(1),置入最新信息x(0)(8),即利用2 ~8 期观测数据建模。以上3 个模型的建模值、预测值与原始监测值对比结果分别见图1、图2。
由图1、图2 可知,利用包含不同信息的样本数据进行建模,其预测精度存有较大差异,新陈代谢GM(l,1)模型预测精度最高,新信息GM(l,1)模型次之,全数据GM(l,1)模型预测精度最差。说明在利用GM(l,1)建模预测时,应着重考虑对新信息的利用。从图1、图2 还可以看出,随着GM(l,1)模型预测步长的增大,误差也逐渐变大,因此,该灰色模型组不适合对巷道围岩变形进行长期预测。
图1 1#点变形监测值与灰色模型群的建模值、预测值比较Fig.1 Comparison of the monitored values with the modeling value and predicted value of the grey model groups at 1# point
图2 2#点变形监测值与灰色模型群的建模值、预测值比较Fig.2 Comparison of the monitored values and the modeling value and predicted value of grey model groups at 2# point
顾及了新信息的GM(l,1)预测模型之所以能够取得较好的预测精度,是因为:
(1)建模过程顾及了变形系统的各扰动因素对预测精度的影响。巷道围岩变形系统可看作1 个发展变化的灰动力学系统,随着煤炭资源的开采以及外部各类地质环境的不断改变,使得岩体变形时刻处于动态变化过程中。为此,在预测建模过程中,及时将最新扰动因素融入到建模数据序列中,能较好地反映变形系统的发展态势。
(2)建模样本顾及了信息量对预测精度的影响。新陈代谢GM(l,1)模型预测精度最高,对于巷道围岩变形系统而言,随着矿产资源的开采及外部环境的改变,历史数据的信息意义逐渐降低,新信息的作用明显增强。在建模过程中,需要不断地融入最新观测数据,同时删除陈旧信息,使得建模序列更好地反映出变形系统当前的变化特征。
(1)软岩巷道围岩变形是一复杂的岩体力学问题,将其看作灰色系统,利用相应模型进行预测,能避开岩体力学的复杂条件,仅利用巷道围岩变形的原始观测数据进行变形预测,建模过程简单方便。
(2)利用GM(l,1)建模,新信息的利用程度对预测精度有较大影响,新陈代谢GM(l,1)模型能充分利用新信息,及时摒弃历史冗余数据,建模过程符合巷道围岩变形的规律,具有较好的预测效果。此外,利用GM(l,1)模型群对软岩巷道围岩变形预测时,短期预测精度较高,但不适合进行中长期预测。
[1] 王永岩.软岩巷道变形与压力分析控制及预测[D]. 阜新:辽宁工程技术大学,2001.
Wang Yongyan. Analysis,Control and Forecast of Deformation and Pressure in Soft Rock Tunnel[D]. Fuxin:Liaoning Technical Universtiy,2001.
[2] 王永岩,马士进,高 菲.软岩巷道围岩变形时序预测方法的研究[J].辽宁工程技术大学学报:自然科学版,2001,20(4):505-506.
Wang Yongyan,Ma Shijin,Gao Fei.Study on time sequence prediction method of the deformation of surrounding rock in soft rock roadway[J]. Journal of Liaoning Technical University:Natural Science Edition,2001,20(4):505-506.
[3] 沈建波,屈建国,詹召伟.基于灰色理论的软岩巷道围岩变形预测研究[J].煤矿开采,2011,16(1):22-23.
Shen Jianbo,Qu Jianguo,Zhan Zhaowei. Prediction of surrounding rock deformation of roadway with soft rock based on grey theory[J].Coal Mining Technology,2011,16(1):22-23.
[4] 吴和平,陈建宏,张 涛,等. 高应力软岩巷道变形破坏机理与控制对策研究[J].金属矿山,2007(9):50-54.
Wu Heping,Chen Jianhong,Zhang Tao,et al. Study on deformation destruction mechanism of high stress weak-rock tunnel and control measures[J].Metal Mine,2007(9):50-54.
[5] 杨利平,张凤杰. 深部软岩巷道支护围岩结构的变化规律研究[J].煤炭技术,2014,33(6):19-21.
Yang Liping,Zhang Fengjie. Deformation regular study on bearing structure of surrounding rock in deep soft rock roadway[J]. Coal Technology,2014,33(6):19-21.
[6] 刘祖强,张正禄,邹启新,等. 工程变形监测分析预报的理论与实践[M].北京:中国水利水电出版社,2008.
Liu Zuqiang,Zhang Zhenglu,Zou Qixin,et al.Theory and Practice of Deformation Monitoring Analysis and Prediction for Engineering Constructions[M].Beijing:China Water & Power Press,2008.
[7] 张向东,吕金伟.特厚煤层大断面巷道支护设计及分析[J]. 金属矿山,2014(8):29-32.
Zhang Xiangdong,Lu Jinwei. Support design and analysis of large section roadway in extremely thick coal seam[J].Metal Mine,2014(8):29-32.
[8] 李男男,顾士坦,于凤海,等. 复杂条件下深井软岩开拓巷道变形破坏分析[J].金属矿山,2012(7):45-49.
Li Nannan,Gu Shitan,Yu Fenghai,et al.Analysis on the deformation of the surrounding rock of development tunnel under the complex conditions on deep tunnels[J].Metal Mine,2012(7):45-49.
[9] 刘思峰,党耀国.灰色系统理论及其应用.[M].5 版.北京:科学出版社,2010.
Liu Sifeng,Dang Yaoguo. Grey System Theory and Its Application[M].The Fifth Edition.Beijing:Science Press,2010.
[10] 党耀国,刘思峰,王正新,等.灰色预测与决策模型研究.[M].北京:科学出版社,2009.
Dang Yaoguo,Liu Sifeng,Wang Zhengxin,et al. Grey Prediction and Decision[M].Beijing:Science Press,2009.
[11] 高 宁. 贫信息小样本条件下时空动态变形预报方法研究[D].北京:中国矿业大学(北京),2013.
Gao Ning.Research on Forecasting Method of Spatial Temporal Dynamic Deformation Under the Condition of Poor Information and Small Sample Data Set[D]. Beijing:China University of Mining and Technology(Beijing),2013.
[12] 霍玉华.隧道围岩变形量预测的灰色模型应用比较研究[J].北京交通大学学报,2006,30(4):42-45.
HuoYuhua.Study on the grey model and its application in forecasting deformation of tunnel surrounding rock[J]. Journal of Beijing Jiaotong University,2006,30(4):42-45.
[13] 吴益平,李亚伟.灰色-进化神经网络模型在深埋隧道围岩变形预测中的应用[J].岩土力学,2008(S):263-266.
Wu Yiping,Li Yawei.Application of Grey-ENN model to prediction of wall-rock deformation in deep buried tunnels[J].Rock and Soil Mechanics,2008(S):263-266.
[14] 高 宁,潘传姣,李建刚.深埋隧道围岩变形预测的非线性组合模型[J].金属矿山,2014(8):45-48.
Gao Ning,Pan Chuanjiao,Li Jiangang.Non-linear combination forecast model to prediction surrounding rock deformation in deep buried tunnels[J].Metal Mine,2014(8):45-48.
[15] 高文华,朱建群,黄自永,等.隧道围岩变形动态预测的灰色自适应模型及其参数智能辨识[J].公路交通科技,2012,29(1):114-120.
Gao Wenhua,Zhu Jianqun,Huang Ziyong,et al.Grey self-adaptive model of dynamic prediction of surrounding rock deformation of tunnel and intelligent identification of parameters[J].Journal of Highway and Transportation Research and Development,2012,29(1):114-120.