基于回归分析的概率积分法参数获取

  • 格式:pdf
  • 大小:491.26 KB
  • 文档页数:3

第29卷第5期 2015年9月 山东理工大学学报(自然科学版) 
Journal of Shandong University of Technology(Natural Science Edition) 
V01.29 No.5 

Sep.2015 

文章编号:1672—6197(2015)05—0061—03 
基于回归分析的概率积分法参数获取 

李 强 
(山东科技大学测绘科学与工程学院,山东青岛266510) 

摘 要:利用回归分析的方法对某煤矿7269开采工作面进行预计参数的确定,并根据实际开采情 
况进行参数修正,以较为快速准确地得到开采面的预计参数,并将预计结果与实际观测数据进行比 
较分析,验证了该方法的预测精度. 
关键词:概率积分法;回归分析;参数计算;误差分析 
中图分类号:P258 文献标志码:A 

The parameters of probability integral method 
acquisition based on regression analysis 

LI Qiang 
(College of Geomatics,Shandong University of Science and Technology,Qingdao 266510,China) 

Abstract:Regression analysis method was performed in a coal mine mining for 7269 to determine 
the parameters,and according to the actual mining parameters were modified to get the expected 
parameters of surface mining more quickly and accurately,and the prediction accuracy of the 
method was verified by comparing the expected results and the actual observation data. 
Key words:probability integral method;regression analysis;parameter calculation;error analysis 

地下矿物被采出后,会打破岩体力学原有的平 
衡状态,应力会重新分布以达到新的平衡状态,当开 
采面积达到一定程度后,这种变形会扩展到地表,使 
地表发生形变,就造成了开采沉陷问题.沉陷现象不 
仅会对矿区的正常开采造成危害,而且会威胁周围 
村庄建筑物的安全,所以对开采可能引起的沉陷作 
出预计是十分重要的.目前,在矿山开采预计的各种 
方法中,概率积分法是较为常用的一种方法,但是如 
果只是根据经验值对预计参数进行选取会存在较大 
的误差,所以应选用较合适的模型进行处理已得到 
更好的参数值[1].本文采用回归分析的方法,获得概 
率积分法的参数值,并与曲线拟合得到的参数值进 
行比较. 

收稿日期:20l4—1l—O8 
作者简介:李强,男,liqiangxin@163.com 

1概率积分法 
概率积分法是由于其所用的变形预计公式里含 
有概率积分而得名的,它主要是基于随机介质理论 
而建立的.其主要思想是:将岩体看作是随机介质的 
颗粒,当矿物被采出时,其上方颗粒介质的运动用颗 
粒的随机运动来表示,大量颗粒介质的移动是符合 
随机过程的.因此可以将整个采区的开采划分为无 
数个“单元开采”,整个采区的沉陷可以看作是无数 
个“小单元”的塌陷所造成的.从统计学的角度看,各 
单元开采影响之和等于整个开采对地表的影响,这 
个叠加和计算过程可以用概率分布密度曲线的积分 
来完成.使用概率积分法对矿山开采沉陷进行预计主 
要需完成以下参数的求取:下沉影响系数g、水平移动 
系数6、开采影响传播系数k和主要影响角 引. 
62 山东理工大学学报(自然科学版) 
2 回归分析模型的选取 
在矿山开采过程中,地表的形变程度与很多因 
素有关,如煤层厚度 、煤层倾角a、采深H、开采方 
法、地表的覆岩情况等.利用文献[3]中对观测数据 
相对较为完整的100个典型工作面作为数值分析的 
样本进行回归分析的结果,可以得到各个参数与地 
质采矿因素之间的函数关系. 
2.1最大下沉值W。 
最大下沉值主要与上覆岩层的性质和深厚比有 
关,其关系式为 
w。 55+ +0.04( sa 。 

式中:H为采深; 为煤层采厚;a表示煤层倾角; E 大小为3 600MPa,E表示该采区的变形模量. 2.2最大水平移动值U。 最大水平位移值主要受采厚m和煤层倾角a 的影响,通过拟合分析得到最大水平位移的计算公 式为 U 一0.06m 。+0.36e一。。 +0.003 式中,v为煤层倾角的弧度值. 2.3 影响半径r 影响半径可以反映出开采的影响范围,它主要 与开采深度H有关,得到关系式为 r一1.95H。・ 2.4开采影响传播角 开采影响传播角0主要取决于覆岩的岩性和煤 层的倾角,在该开采区选取的是中硬(覆岩强度介于 30MPa与60MPa之间)覆岩,计算公式为 0—221.8a-0- 一10】.7 3 矿区概况 某煤矿井田地貌属黄淮冲积平原,为第四系地 层覆盖地区,地势较平坦,陆地地面高程33.54~ 37.47m,东部湖区湖底高程为3O.O0~33.00m,该 采区的地质状况见表1.该矿区7269工作面自2010 年10月25日开始回采,2011年9月13日回采结 束.工作面走向长度852 ̄876m,平均864m,倾向长 度170m,煤层倾角7。~12。,平均9。,煤层采厚 5.0m,平均采深780m.根据该开采区的地质状况报 告,可知预计时应选用的变形模量为3 900MPa. 表1采区7269工作面地质情况 3.1采区参数计算 根据上述采区已有的资料和回归分析模型中各 参数的计算公式,可计算得到该采区充分采动的预 计系数,结果见表2、表3. 
表2基于回归分析的计算值 

表3 采区7269工作面的预计参数的计算结果 
根据表3可以得到7269工作面充分开采条件 
下的各个预计参数值,但是根据实际的开采状况可 
知,该采区为不充分采动,所以要对得到的预计值进 
行一定的修正【4].下沉系数修正采用的公式为 
g1一q(0.97n 一0.07n+0.39) 
水平移动系数修正采用的公式为 
b1一b(0.9125n 一0.5375n+0.836) 
在修正式中q,b为计算出的充分采动下的下沉系 

数, = ̄/n1 2, 1===D1/(2r), z=D2/(2r), 
D ,D 为实际开采的走向和倾向长度,r为主要影 
响半径[5].经过计算得到修正后的下沉系数q=== 
0.30,水平移动系数6—0.42,这样就可以得到完整 
的一套预计参数. 
3.2预计成果分析 
根据该采区布设的实测观测线资料,采用曲线 
拟合的方法对7269开采工作面进行参数确定,得到 
的预计参数与回归分析结果见表4.