采煤塌陷区土地破坏面积和破坏程度的预测.doc

第期矿山测量№ 3 .3年月 20019 MINE SURVEYING Sept.2001 21 采煤塌陷区土地破坏面积和破坏程度一般可用地表移动变形预计法和地表移动范围角量参数法进行预测。 应用地表移动变形预计法预测地表塌陷面积及破1 坏程度 采用山区地表移动预计软件进行预计,该软件对于平地即为平常的概率积分法,对于丘陵和山区则在概率积分法的基础上附加地形引起的采动滑移影响。在地形、地质和开采条件已知的情况下,可以预计开采影响范围内平地或山区地表任意点任意方向的下沉和水平移动及变形值。预计公式如下 1 2 上列各式中的A 、P 、t 为滑移影响函数参数,可根 据矿区地表移动观测资料按多元函数回归方法求取。山西有关矿区求得的A 2π,P , 2t π。上式中的W m 为该地质采矿条件下的最大下沉值 W max 以计。mm W x 、W y 和W x, y 分别为主剖面、X Y 上 x,y 点和平面上 x,y 点因开采引起的地表的下 沉值,可分别按概率积分法公式计算。 地表任意点的下沉和水平移动 W x , y 和 式中W x , y 和U x , y φ 分别为地表任意点x , y 的下沉值和任意点、任意方向φ的水平移动值。α为 x , y 点的地形趋势面倾角,可按修匀后 的地形剖面取值。 φ为 x , y 点的地表倾斜方向 角;φ为计算方向角 ; φ、 φ均由x轴正向逆时针方向计算。式中的D 为 x , y 点的地表特性系 数,可按表取值。式中的1P x 、P y 和 P x,y 分 别为 X 、Y 主剖面上 x 、y 点和平面上 x , y 点的滑移影响函数,可按公式计算3-5 U x , y φ求出后其任意点垂直变形包括倾斜 , i x ,y φ、曲率K x , y φ和水平变形ε x ,y φ可 按给定的点间水平距离d ,按数值导数规则计算点 间的平均垂直变形与水平变形。设相邻点编号为、、i j k,则有 倾斜 mm/m 6 W x , y W x ,y D x , y P x , y W x , y tg 2α′x , y U x , y φ U x , y φD x , y [P x cos φx ,y cos φ P y sin φx , y sin φ]W x , y tg αx , y ′ x , y x , y x , y x , y P x , y P x cos 2 φx , y P y sin 2 φx , y P x P y cos 2 φx , y sin 2 φx , y tg 2 α′ x , y 5 采煤塌陷区土地破坏面积和破坏程度的预测 胡海峰何万龙康建荣 太原理工大学 摘要 文中在已有的山区地表移动规律研究成果和采煤塌陷区土地破坏面积调研的基础上,对采煤塌陷区土地破坏面积和 破坏程度的预测方法作了研究,提出了一套塌陷区土地破坏面积和破坏程度的预测及分类方法。 关健词 采煤塌陷土地破坏 变形值移动角 P x 1A exp[-12 x r P 2 ]W m exp[-t x r P 2 ] 3 P y 1A exp[-12 y r P 2 ]W m exp[-t y r P 2] 4 地表类型表土层与地面植被特征 特 性 系 数 D 凸 形 地 貌凹 形 地 貌 Ⅰ沙质粘土层,厚度02 m , 地面有密集的灌木丛或树林 0.30.5-0.2-0.3 Ⅱ沙质粘土层,厚度25 m , 地面有密集的灌木丛和疏林的荒坡 0.51.0-0.3-0.5 Ⅲ亚粘土或坡积物,厚度大于5 m ,地面为果园或耕地 1.01.5-0.5-0.8 Ⅳ 具有垂直节理的湿陷性黄土或坡积物,厚度大于5 m ,地面为耕地 1.01.5 -0.5-0.8 表山区地表特性系数1 D i x , y φij W x , y j -W x , y i d ij 开采沉陷 开采沉陷 第期矿山测量年月 3 2001922 曲率 10-3 /m m 7 水平变形 mm/m 8 表矿区地表移动预计参数3 由开采沉陷学可知,对于平地一般单一开采 的缓倾斜矩形工作面,产生一个略偏于下山方向的槽形盆地,其地表移动边界和下沉等值线近似于长椭圆形。当多个工作面连续开采时,地表将呈现连续的槽形塌陷盆地和连续的长椭圆形移动边界。在槽形盆地接壤的边界部位,将产生正曲率k 和拉 伸变形ε叠加区,叠加后的正曲率和水平拉伸变形可能大于半无限开采相应的最大变形值但最多不会超过一倍。当单一工作面倾向开采不充分,其槽形盆地中心倾向主断面上的负曲率-k 和压缩变形-ε绝对值可能大于半无限开采相应的变形最大值,但最多也不会超过最大值的一倍。上述半无限开采的最大变形值可按下列公式计算 最大倾斜在处 ; x 0最大曲率在 ; ≈处x 0.4r 最大水平移动在处 ;x 0最大水平变形 ; 在≈处 x 0.4r 式中为最大下沉值为主要影响半径, ; r 见式为水平移动系数。取中硬覆岩条件10; b 下长壁大冒顶开采的地表移动参数下沉系数 q 依此类推,即可算出所有点和点间的垂直与水平变形。按点间水平距离d 计算变形值时,应由表按2采深H 选取d 值。 预计参数应以工作面为单元输入,按矿区部分观测资料和类比法求得的预计参数如表。 3,水平移动系数0.850.90b ,主要影响角正 0.33切β,拐点偏距可算出不tg 2.2 2.4 , 同采深、采厚的水平煤层开采平地的最大移动变形值如表。 4以上计算都是按平地的预计结果。山区由于受地形和采动滑移附加移动变形的影响,其移动变形值较平地可能有所不同,概括起来可能有以下特点⑴山区采动滑移方向总是指向坡体的下坡方 向,因而下坡方向的水平移动量将增加,而上坡方向的水平移动量将减小。 ⑵在凸形地貌和凸形变坡部位,地表因受附 加拉应力影响而产生附加下沉和水平拉伸,使地表的下沉和水平拉伸变形值增大,同时还可能产生较大的正、负曲率变形,从而使土地的破坏程度有所加重。 ⑶在凹形地貌和凹形变坡部位,地表因受地 形附加压应力影响而产生附加上升和水平压缩变形,使地表的下沉值略有减小裂缝宽度和密度也;略有减小,因而土地的破坏程度略有减轻。 预计地区的地表移动与变形值算出后,即可参照表的移动变形指标和地形开采条件参数确定地5表塌陷的面积和破坏等级。表的指标是根据山西5省塌陷区土地破坏等级的有关裂缝指数,以及矿区地表变形与土地裂缝的相关规律确定的。 应用矿区地表移动角量参数预测地表塌陷面积与2 破坏程度 2.1 地表塌陷范围和塌陷面积的确定 矿区地表塌陷范围可按矿区地表移动角量参数中的边界角、移动角或裂缝角圈定。 移动边界角是在主断面上按地表移动边界下沉为的点至开采边界连线在煤柱一侧与水平线 10mm 所夹角度,煤层走向以δ 0表示,下山和上山分别 K x , y φj i x , y φjk -i x , y φij 0.5d ij d jk ε x , y φij U x , y φj -U x , y φi d ij i max W max r k max 1.52 W max r 2U max bW max 开采深度H ,m 500点间水平距离d ,m 5 7.5 10 15 20 25 30 表按数值求导规则计算地表变形值的点间平距2 d 参 数 名 称 及 代 号参数取值表达式下沉系数q q0.80.1水平移动系数b b0.330.08主要影响角正切tg βtg β2.00.4拐点偏移距S i S i 0.10.05H i 开采影响传播角θ θ90-α α≤3 θ90-0.6 α3 εmax 1.52b W max r W max S i 0.1H i 第期胡海峰等采煤塌陷区土地破坏面积和破坏程度的预测年月 3 2001923 以β0和γ0表示。地表综合移动角不分基岩与表土层是在主断面上按临界变形值i 0、 3 mm/mK 00.210-3、ε/m0确定的点至开采边界2mm/m的连线在煤柱一侧与水平线所夹角度,走向为δ,下山和上山分别为β和γ。地表裂缝角是在主断面上地表塌陷最外边缘的地表裂缝位置与开采边界连线在煤柱一侧与水平线所夹角度,走向为δ“,下山和上山分别为β“和γ“。山西除大同以外的石炭二迭纪煤田中硬覆岩矿区采厚的上述角量参2 6 m 数大致如表。 6目前国内划定地表塌陷范围有的按边界角,有 的按移动角,对于低潜水位地区以移动角或裂缝角划定地表塌陷范围较为合适。因为这些矿区已利用土地的表层大多为黄土所覆盖。现场观测表明,黄土在水平拉伸变形达到左右即可出现微小的2mm/m 裂缝,故一般黄土覆盖的平坦地区应按移动角划定 表山西石炭二迭纪煤层开采地表移动角量参数参考值 6 地表塌陷范围。对于其它丘陵地区,一般可按裂缝角划定地表塌陷范围。 丘陵山区因为高差变化较大,故在井上下对照图上按移动角或裂缝角划定地表塌陷范围时,应用趋近法进行,如图所示。先按开采边界在井上下 1 表不同采深采厚半无限开采地表最大移动变形值 4 序号采厚 M m采深H m 最大下沉W max mm最大平移U max mm最大倾斜 i max 10-3 最大曲率k max 10-3/m最大水平变形 εmax 10-3 土地破坏程度1 2.3100195564543.0 1.4421.6中重2 2.6200221072924.30.4112.2轻中3 3.5100297598265.4 2.1932.8中重4 3.72503145103828.90.4014.5轻中5 4.61503910129060.0 1.4030.1中重6 4.62503910129035.90.5018.0中7 5.52004590163356.9 1.0028.5中8 6.510055251823121.5 4.0660.9重9 6.81506120202097.9 2.3849.1重10 6.52005525182360.8 1.0230.5中重11 6.53005525182340.50.4520.3中127.83007020231756.20.6828.2中重13 6.54005525182330.40.2515.2轻中14 6.55005525182324.30.1612.2轻15 8.5 500 7650 2524 33.7 0.22 16.9 轻中 表矿区地表移动变形值与土地破坏等级对照表5 破坏等级地表变形值 深厚比 H/M 宽深比 L/H 微地貌 水平拉伸变形 εmm/m 倾 斜 i mm/m曲 率 K 10-3/mⅠ级轻度 210 320 0.20.5 45 20 40 1.0 1.0 平地、坡地和凸形地貌 开采深度Hm边界角δo 移动角δ 裂缝角δ“ 松散层移动角厚度 m 干燥含水1006065706060 55 300400707575787880400500757678808085 备 注 1.α50时,βδ-0.6α,γδ。 第期矿山测量年月 3 2001924 图开采边界至塌陷边界水平距离计算示意 1 D 对照图上划出主剖面线X 和Y ,然后按各主剖面边 界的开采深度H 0和移动角分别按下列趋近公式确定开采边界至塌陷边界的水平距离D 下山 ;上山 ;走向 ; 式中 求得D 之后,由各主断面的开采边界向外量取各D 值,即可得各主断面的塌陷边界点位置,然后将各边界点连结起来即得到地表塌陷边界。 矿区地表移动观测资料分析和现场调查表明,在一定开采条件下,地表永久性裂缝充分发育区域一般是从拐点至开采边界的中点起,向煤柱一侧O 扩展2.5b 的范围,其水平宽度 d 大致可按下列公式计算虽然总的永久性拉伸或裂缝范围应从拐点位置起算,至按裂缝角确定的裂缝边界点为止,但 0 两拐点之间的部分则应视为开采引起的动态变形拉伸区。如果是平地,动态拉伸区的裂缝可能在开采沉陷稳定后逐渐自行闭合,但在丘陵和山区,由于附加水平拉伸的影响,位于凸形地貌和凸形变坡部位山顶和梁峁的边缘的动态裂缝不但不能闭合,相反其宽度、密度和落差可能会变得更大,而凹形地貌部位沟谷内由于附加压缩变形的影响,其地表塌陷范围划定后,即可用求积仪计算塌 陷范围的面积。每一图斑的面积应量取二次,二次面积互差在允许范围内时,可取平均值作为最后的结果。 2.2 地表塌陷土地破坏程度的确定 山区地表塌陷破坏主要取决于地表裂缝的宽 度、落差和密度状况,而地表裂缝状况则与地质采矿条件及微地貌形态有关。当一个矿区的覆岩 和地表岩土性质基本相同且不存在特殊断裂构造的情况下,水平缓倾斜煤层单一矩形工作面开采引起的地表永久性静态裂缝主要分布在地表移动水平拉伸变形区,如图所示。如不考虑地形和2地貌影响则在充分或接近充分开采时地表最大水,,平拉伸变形距拐点的水平距离为a 在一般情况下,中硬覆岩矿区的主要影响角正切β;拐点偏距,因而最tg 2.20.2 大水平拉伸变形距开采边界外侧的水平距离应 为动态裂缝将不能得到充分发育,其宽度、密度和落差将有所减小。 综上所述,按地质开采条件和微地貌形态确 定塌陷区土地破坏等级时,应遵循下列原则 初次采动永久性拉伸范围内一般平坦地区 1 和凸形地貌耕地的裂缝发育区破坏等级可按地质采矿条件划分,其中深厚比和水平拉伸变形为主要指标,其它为参考指标。 裂缝发育区以外的永久性裂缝区和动态裂 2 缝区的破坏等级应比裂缝发育区降低一级。 同一区域内凸形地貌变坡部位的破坏等级 3 应比该区域的破坏等级提高一级;凹形地貌部位的破坏等级则应降低一级。下转第页 37D β H 0h 1h 2h n βctg β D γ H 0h 1h 2h n γctg γ D δ H 0h 1h 2h n δctg δ h 1 H 1-H 0 ; h 2 H 2-H 0 ;h n H n -H 0 。 a i 0.4r i 。 s i 0.1H i b i b i a i S i 0.4H i tg β -0.1H i 0.0820.2 H i d i 2.5b i 0.20.02H i S 0 b S 0 a d b a d 1 2 1S 0 b S 0 a d b a d 1 2 11永久性裂缝 d永久性裂缝发育区 2动态裂缝区 图2 塌陷区地表裂缝区划分示意 3 20019第期矿山测量年月 2937