Application of comprehensive survey methods in the gypsum mine goaf for Mengxi-Huazhong railway route selection
-
摘要: 以蒙华铁路荆沙段石膏矿采空区勘察为例,介绍了采用遥感解译、调查访问、勘探验证相结合的采空区综合勘察方法。针对铁路穿越采空区选线重点关注的问题,提出了采空区勘察的工作思路,阐述了卫片、航片、多时相分析的采空区遥感解译方法,提出采空区调查访问时应重点关注的问题及如何开展综合分析,并提出应根据区域地质条件、采空区发育特征、勘探方法的适用性、场地条件等,合理选择勘探验证方法。最后,通过实例分析,介绍了金陵、宝安石膏矿采空区的综合勘察,可为类似采空区勘察提供借鉴和参考。Abstract: Taking the investigation of the mined-out area of gypsum mine in Mengxi-Huazhong railway Jingmen-Yueyang section as an example, this paper introduces the comprehensive survey method of the mined-out area by the combination of remote sensing interpretation, investigation and visit and exploration verification. In view of the problems of railway route selection through goaf, the working idea of goaf survey is put forward, and the remote sensing interpretation method of goaf based on satellite photograph, aerial photograph and multi-temporal analysis is expounded. This paper puts forward the problems that should be paid attention to in the investigation and visit of goaf and how to carry out comprehensive analysis. It is also proposed that the exploration verification method should be reasonably selected according to the regional geological conditions, the development characteristics of the goaf, the applicability of exploration methods, site conditions and so on. Finally, the comprehensive survey of goaf in Jinling and Baoan gypsum mines is introduced, which can provide reference for similar goaf investigation.
-
Keywords:
- gypsum mine /
- goaf /
- remote sensing interpretation /
- geophysical exploration /
- comprehensive survey
-
0. 引言
滑坡是威胁山西地区人民生命和财产安全的主要地质灾害之一,广布的黄土层为滑坡灾害的发生提供了物质基础,尤其是汾河流域多为非自重的湿陷性黄土,具有显著的结构疏松与压缩性高的性质[1 − 2],而土质边坡以易发生、危害强、难预测等特点危害性更为突出[3 − 4]。滑坡现象多为内部因素和外部因素共同作用触发,内部因素包括坡体本身的性质弱化以及较差的地质构造环境等,降雨、工程开挖等外部因素是诱发土质边坡失稳的直接原因[5 − 6]。
众多学者围绕土质边坡稳定性问题开展了一系列卓有成效的研究,尤其是20世纪70年代后,有限元法借助计算机开始被广泛应用于边坡稳定分析。然而,在实际应用中,有限元法仍存在一些局限性,如较难确定边坡的初始应力状态以及把握边坡临近破坏时的弹塑性本构关系等。极限平衡法逐渐形成了垂直条分法和滑移线法两个独立的分支[7 − 8]。我国学者在滑动机理及监测方面也取得了一系列研究成果,彭建兵等[9]研究了黄土高原滑坡灾害形成的动力学机制问题,认为区域构造应力、边坡构造应力及易灾特性是黄土滑坡的“三大元凶”;董震等[1]聚焦山西省黄土地质灾害问题,认为黄土地质灾害的发生受其特有的物质成分和工程性质控制,并具有地域性分布的特点;支泽民等[10]研究发现,边坡稳定性影响因子相互作用会促进滑坡发育;蔡欣育等[11]开展了不同降雨类型对边坡稳定性的影响研究,重点分析了不同降雨类型条件下边坡孔隙水压力的分布特征;侯天顺等[12]也进行了人工降雨作用下的滑坡试验,模拟长期连续降雨边坡稳定系数的变化规律;马蓓青等[13]探讨了降雨条件下边坡裂缝发育规律及其对稳定性的影响;苗青等[14]通过LASSO算法建立了能够预测滑坡危险区域的模型。
然而,在工程扰动下土质边坡稳定性愈加复杂,很难用一种模式概括其滑动机理与影响因素。调查表明,多数土质边坡的滑坡与降雨有直接或间接的关系。在近年极端天气频发的背景下,探讨工程场地客观地质条件下极端降雨对土质边坡稳定性的影响是十分必要的。本文在山西太原古交市马兰镇营立村拟建水泥厂地质调研基础上,基于模糊数学理论获取该土质边坡稳定的主要影响因素,确定权值并计算稳定系数,结合FLAC3D软件模拟得到不同降雨强度下边坡稳定系数变化特征,进而得到现有条件下维持该边坡稳定所能承受的极端降雨条件,从而综合判断边坡的稳定状态。
1. 工程地质条件
1.1 地形与岩性特征
研究区位于山西省太原市西部山区的古交市马兰镇营立村水泥厂,属于地形起伏较大的山间河谷地貌单元区,拟建场地地形起伏较大,各钻孔孔口标高介于1182.43~1194.64 m,相对高差为12.21 m。区内地层发育不完全,据勘察及钻孔揭露的地层情况,勘察深度范围内地层主要由全新统粉土、圆砾与中奥陶统石灰岩组成(图1)。粉土(Qh)呈褐黄色,稍密,稍湿,土质较松散;场地局部出露,可见零星煤屑、云母、氧化铁及植物根系,摇振反应中等,无光泽反应,干强度和韧性低。圆砾(
)为杂色,中密,稍湿,母岩成分主要为灰岩和砂岩,呈椭圆状,磨圆度较好,一般粒径为3~10 cm,最大粒径为18 cm,以粉土、中砂充填其中,含大量漂石,中奥陶统马家沟组石灰岩(O2m)为灰白色,强风化-中风化,岩体节理裂隙较发育,岩芯较完整,呈短柱状或柱状,一般在10.0~30.0 cm,敲击声较清脆,不易击碎,为较硬岩。1.2 水文地质与气候
场地内水文地质条件简单,在勘察深度内未发现稳定地下水位。研究区属于典型的季风气候,当地历年降雨资料显示年平均降雨量为426.1 mm,降雨主要集中在7—8月,其降雨量约在60~90 mm,但近年极端天气明显增多[14],多次出现了日降水量超过100 mm的大暴雨,极端天气频现给区域内土质边坡稳定性带来了挑战,也增加了研究的复杂性和难度。
2. 边坡稳定性模糊数学分析
2.1 边坡稳定性评判模型
边坡稳定性受到内因和外因的双重影响。在现场调研基础上,采用层次分析法确定了影响该水泥厂边坡稳定性的因素及其分级标准(表1),确定了地层岩性、地貌特征、风化作用、降雨及人类活动[15]等 5 项影响因素,构建了边坡稳定性的模糊综合评价模型
表 1 稳定性影响因素与量化取值表Table 1. Stability influencing factors and quantified values评价指标 滑坡稳定性级别及其分级标准 因素 因子 稳定(1) 基本稳定(2) 欠稳定(3) 不稳定(4) 降雨U1 降雨量/
(mm·d−1)25~<50 50~<100 100~250 >250 人类活动U2 工程活动 弱 一般 较强 强 地貌特征U3 坡面形态 平直形 凸状 凹状 阶状 坡度(θ)/(°) 0~<10 10~<30 30~<50 50~90 坡高(h)/m 0~<25 25~<35 35~<45 45~55 地层岩性U4 岩土性质 坚石 软石 硬土 普通土 风化作用U5 风化程度 微风化 中风化 强风化 全风化 根据构建的边坡稳定性模糊综合评价模型[16],建立评价对象的关联因素集,表示为:
(1) 根据边坡实际情况,结合半定量取值法,确定不同因素权重集,表示为:
(2) 通过评价对象可能出现的结果建立备择集,表示为:
(3) 根据单个影响因素,联立式(2)与(3),确定各离散值对备择集的隶属程度,得出评判矩阵,表示为:
(4) 根据式(2)与(4),依据权重集与影响因素评判关系矩阵的乘积可以有效地表达所有因素对判定对象的综合影响指数,表示为:
(5) 2.2 边坡稳定影响因素关系矩阵
基于地质资料及调查结果,对各影响因素在不同稳定系数下的影响程度进行评价。基于所构建边坡稳定性评价模型的降雨量、人类活动、地貌特征、地层岩性与风化程度等5项因素建立关系矩阵,矩阵主对角线上数据表示各影响因素在不同稳定系数下的影响程度。矩阵中各行表示该因素对其它因素在不同稳定系数下的影响程度,设该因素为A。矩阵中各列表示其它因素对其在相同稳定系数下的影响程度,设其为B。不同因素对不同稳定系数的影响具有差异性,因而选用编码方式加以区分。按影响程度由弱到强分别赋值为0,1,2,3,4。其中,Aa为关系矩阵中不同稳定系数的稳定性影响程度,该值越大表示影响越大;Bb为关系矩阵中相同稳定系数的稳定性重要程度,该值越大表示越重要。
基于以上原则建立了关系矩阵,计算公式如下:
(6) (7) (8) 式中:Aa——稳定性影响程度;
Bb——稳定性重要程度;
Hi——影响因素的权重(表2)。
表 2 营立村水泥厂切坡稳定性影响因素赋值表Table 2. Impact factors of artificial slope cutting stability in Yingli Village quarryP A A+B A-B Hi/% Aa Bb P1 3 2 1 1 11 23 −1 26 16.26 −0.707 3 P2 2 2 1 10 19 1 22 13.44 0.707 2 1 P3 1 1 7 15 −1 17 10.06 −0.707 2 2 1 P4 3 10 17 3 19 12.02 2.120 1 1 1 1 P5 6 14 −2 16 9.90 −1.410 B 12 9 8 7 8 88 根据地质调查结果及相关地质资料与规范,由式(6)—(8)计算,最终确定各影响因素的权重,降雨为26%、人类活动为22%、地貌特征为17%、地层岩性为19%、风化程度为16%。确定权重集为:
一般土质边坡的稳定系数区间为[1.20, 2.00],区间为0.25均匀分布,建立备择集表示如下:
由此,计算得出综合评价结果为:
最后,对综合评价结果加权平均求得边坡稳定系数,计算结果如下:
式中:Fs——边坡稳定系数。
综上可知,基于模糊数学求得Fs =1.429>1.15,根据《滑坡防治工程勘察规范》(GBT 32864—2016) [17]第13. 3. 4 条滑坡稳定状态划分标准,该土质边坡天然条件下处于稳定状态。
3. 边坡稳定性FLAC3D分析
3.1 模型构建与物理力学参数获取
基于模糊数学分析已知降雨是营立村水泥厂土质边坡稳定性的主要影响因素,为进一步量化分析该土质边坡稳定性,使用FLAC3D 软件模拟不同降雨条件下该边坡的稳定性演变特征[18]。模型主要依据边坡的坡高、坡角等几何尺寸建立,由于边坡开挖为阶梯状,z轴分别为8.5 m、13.5 m、25.3 m,x轴长度为42 m,y轴长度为40 m。模型初始边界条件为:流体边界设置边坡表层为渗透面,坡脚和边坡底部为渗漏面;位移边界固定边坡底面,约束模型两侧的y轴位移。此外,为便于分析坡体在不同降雨条件下各层岩土体弹性模量、容量、剪应力、摩擦角等物理力学参数的变化情况[19],本次模拟在所建模型每级台阶坡脚处在对应台阶内部提取位移监测点与应力监测点(图2)。
岩土体物理力学参数基于地质调查结果及土工试验确定(表3),所建边坡模型采用Mohr-Coulomb准则和强度折减法求取稳定系数[20]。潜在滑动面主要计算原理如下:
表 3 边坡岩土体物理参数指标Table 3. Physical parameter indicators of slope rock and soil类型 容重
/(kg·m−3)泊松比 黏聚力
/MPa内摩擦角
/(°)抗剪强度
/MPa饱和粉土 1960 0.3 0.02 18 0.35 粉土 1820 0.25 0.04 25 0.51 砂岩 2550 0.25 8.2 35 8.91 石灰岩 2800 0.18 11.4 38 12.18 (9) (10) 式中:c、c′——边坡土体初始状态和临界状态下的黏 聚力/kPa;
φ、φ′——边坡土体初始状态和临界状态下的内摩 擦角/(°);
α——最危险滑动面与最小主应力夹角/(°);
δ1、δ3——土体单元的第1和第3主应力/kPa。
3.2 天然状态下边坡稳定模拟分析
基于表3的岩土体力学参数采用FLAC3D数值模拟得到营立村水泥厂土质边坡在天然状态下边坡位移和剪应力空间分布图(图3)。图3(a)表明,天然状态下边坡的位移集中在边坡顶部区域,最大位移为29.3 cm,自坡顶向坡脚位移量具有减小趋势,直至边坡中下部无位移量。图3(b)表明,边坡内部应力变化最为显著,坡体在自重及岩性影响下剪应力集中分布于石灰岩及其上覆第四系沉积层界面附近,剪切力增量向边坡底部集中,最大剪切力累积在坡脚位置,说明在没有断裂的情况下马家沟组灰岩与其上覆第四系沉积层界面存在明显的渗透系数差,是最主要的结构面。模拟分析结果显示,此状态下边坡Fs=1.426,与模糊数学分析法得出的Fs=1.429几乎一致,综合判定该边坡天然状态下处于稳定状态。
3.3 降雨状态下边坡稳定性变化
降雨会引起土体黏聚力和内摩擦角改变,进而影响边坡稳定性[21 − 23]。为了简化研究,将均匀雨型日降雨量设为基础分析条件。FLAC3D 模拟不同降雨强度条件下该边坡的稳定性,结果显示,随着日降水量增加稳定系数小幅升高后下降,整体呈现单峰型曲线形态(图4)。由图4可知,较小的降雨强度对边坡渗流场影响有限,基本无法引起地下水位变化,而雨水渗透引起边坡粉土的黏聚力增大超过了内摩擦角减小的效应,因此边坡稳定性出现了小幅升高。随着日降雨量继续增加,土体黏聚力与内摩擦角均为下降趋势,导致稳定系数快速下降,当降雨量增长至75 mm/d时Fs=1.15,边坡处于稳定状态;当降雨量增长至120 mm/d时Fs=1.05,边坡处于基本稳定状态;当降雨量增长至150 mm/d时,稳定系数为1.004,此时边坡已处于极限平衡的欠稳定状态,存在发生滑坡的潜在风险。
3.4 极端降雨条件下边坡稳定状态分析
150 mm/d的极端降雨引起了潜在滑动带由坡体表面逐渐向坡体深部转移,直至滑动带沟通了岩土界面,计算模型中出现了一条明显的塑性剪切带,见图5(a)。降雨约为0.5 h时,雨水快速入渗引起土体自重增大,变形区域主要集中在表土层,局部具有轻微向外隆起的趋势,Y2、Y3监测点应力随之上升,分别为178.2 kN、132.4 kN;降雨约为3 h时,因入渗土体中水不能及时排出,孔隙水压力抵消一部分土体压力,内部的应力则呈现减小的趋势,Y2、Y3监测点应力分别为123.5 kN、53.5 kN;此后,Y3监测点应力基本不变,Y2监测点经历快速下降(3~9 h)后进入缓慢下降阶段,约至20 h时,Y2监测点应力稳定为59.8 kN。处于坡脚岩土界面附近的Y1监测点应力随降雨历时整体呈上升趋势,变化趋势大致呈三个阶段,6 h内为应力较快上升阶段,6~16 h为应力缓慢上升阶段,16 h后应力增长达到最大限度,最大应力增量约为78 kN,见图5(b)。随降雨历时增加,边坡顶面浅表层首先产生了拉裂缝,裂缝产生的位移量由顶部表层向内部减小,即相同降雨历时条件下,监测点位移量W4>W3>W2>W1。降雨约18 h时,各监测点位移量达到最大值,见图5(c),边坡已基本处于极限平衡状态。
4. 结论
综上可知,降雨对土质边坡稳定性有着显著的影响,雨水入渗导致土体自重增大、强度减小,使边坡内部发生局部剪切作用。随着降雨强度增大,边坡岩土层剪切面贯通形成完整的滑动面,最终导致土质边坡失稳。
(1) 基于地质资料及野外调研基础上,构建了边坡稳定性模糊综合评价模型,结合关系矩阵确定降雨为营立村水泥厂边坡稳定性最主要的影响因素(26%),并求得Fs =1.429>1.15,表明该土质边坡天然条件下处于稳定状态。
(2) 使用FLAC3D 模拟分析结果显示天然状态下边坡Fs=1.426;当降雨量75~120 mm/d时,Fs由1.15降至1.05;降雨量120~150 mm/d时,Fs由1.05降至1.004,此时边坡已处于极限平衡的欠稳定状态,存在发生滑坡的潜在风险。
(3) 模糊数学与FLAC3D 评价营立村水泥厂土质边坡天然状态下稳定性结果几乎一致,联合运用两种方法相互佐证能够提高评价准确性,使结果更符合边坡实际情况。
-
表 1 采空区调查访问内容
Table 1 Investigation and interview contents
访问对象 访问内容 矿上工作人员 开采方法
矿体分布特征,包括矿石类型、矿体埋藏深度、倾向、倾角、是否有断层切割、开采厚度等
开采方向
地下采空分布特征及处理方法
抽排地下水情况附近居民 是否有爆破震动影响
房屋是否有开裂现象
井水是否有减少或干涸现象
地表是否有采空塌陷
耕地灌溉是否有影响 -
[1] 张立琴, 李国和, 苗庆库. 古窑及近代小型采空区勘察方法选择与应用[J]. 铁道勘察,2014,40(4):58 − 60. [ZHANG Liqin, LI Guohe, MIAO Qingku. The surveying methods and application for the ancient and modern small-scale goaf[J]. Railway Investigation and Surveying,2014,40(4):58 − 60. (in Chinese with English abstract) DOI: 10.3969/j.issn.1672-7479.2014.04.017 [2] 张纪星, 师修昌. 浅埋采空区大采高条件下覆岩破坏规律[J]. 中国地质灾害与防治学报,2019,30(5):92 − 97. [ZHANG Jixing, SHI Xiuchang. Failure of overburden rock under large mining height in shallow buried goaf area[J]. The Chinese Journal of Geological Hazard and Control,2019,30(5):92 − 97. (in Chinese with English abstract) [3] 管永伟, 魏义强, 杨振江, 等. 基于数值模拟的房柱式采空区关键层-岩柱系统破坏机理分析[J]. 中国地质灾害与防治学报,2017,28(4):59 − 63. [GUAN Yongwei, WEI Yiqiang, YANG Zhenjiang, et al. Failure mechanism of the key strata-pillar system in room and pillar goaf based on a numerieal simulation[J]. The Chinese Journal of Geological Hazard and Control,2017,28(4):59 − 63. (in Chinese with English abstract) [4] 李国和, 李桂芳. 采空区铁路工程地质选线研究[J]. 铁道工程学报,2012,29(10):15 − 20. [LI Guohe, LI Guifang. Research on railway alignment in goal[J]. Journal of Railway Engineering Society,2012,29(10):15 − 20. (in Chinese with English abstract) DOI: 10.3969/j.issn.1006-2106.2012.10.004 [5] 邓睿, 杜宇本, 蒋良文, 等. 玉磨铁路采空区段地质选线研究[J]. 铁道工程学报,2017,34(2):21 − 25. [DENG Rui, DU Yuben, JIANG Liangwen, et al. Research on the geologic route selection in goafs of Yuxi-mohan railway[J]. Journal of Railway Engineering Society,2017,34(2):21 − 25. (in Chinese with English abstract) DOI: 10.3969/j.issn.1006-2106.2017.02.005 [6] 祝佳. Landsat8卫星遥感数据预处理方法[J]. 国土资源遥感,2016,28(2):21 − 27. [ZHU Jia. Analysis of Landsat8 satellite remote sensing data preprocessing[J]. Remote Sensing for Land & Resources,2016,28(2):21 − 27. (in Chinese with English abstract) [7] 张刚, 李忠权, 李佳奇, 等. 基于OLI青海文巴地区遥感图像增强及构造解译[J]. 矿产勘查,2015,6(5):571 − 576. [ZHANG Gang, LI Zhongquan, LI Jiaqi, et al. Remote sensing image enhance and structural interpretation based on OLI in the Wenba area, Qinghai[J]. Mineral Exploration,2015,6(5):571 − 576. (in Chinese with English abstract) DOI: 10.3969/j.issn.1674-7801.2015.05.011 [8] 强建华, 于浩. 新疆矿山环境遥感监测成果综述[J]. 中国地质调查,2016,3(5):28 − 34. [QIANG Jianhua, YU Hao. Review on remote sensing monitoring results of mine geological environment in Xinjiang[J]. Geological Survey of China,,2016,3(5):28 − 34. (in Chinese with English abstract) [9] 陈文平, 范英霞, 韩小明, 等. 中巴资源卫星影像HR数据在煤矿矿区地面塌陷调查中的应用[J]. 测绘与空间地理信息,2012,35(2):80 − 83. [CHEN Wenping, FAN Yingxia, HAN Xiaoming, et al. Application of CBERS HR image data in investigation of surface collapse in coal mine zone[J]. Geomatics & Spatial Information Technology,2012,35(2):80 − 83. (in Chinese with English abstract) DOI: 10.3969/j.issn.1672-5867.2012.02.023 [10] 国家铁路局. 中华人民共和国铁路运输行业标准: 铁路工程地质遥感技术规程 TB 10041—2018[S]. 北京: 中国铁道出版社, 2018. Railway & Train Standard of the People's Republic of China: Code for technique on remote sensing of railway engineering geology TB 10041—2018[S]. Beijing: China Railway Publishing House, 2018.(in Chinese)
[11] GORELICK N, HANCHER M, DIXON M, et al. Google Earth Engine: Planetary-scale geospatial analysis for everyone[J]. Remote Sensing of Environment,2017,202:18 − 27. DOI: 10.1016/j.rse.2017.06.031
[12] 王宾, 韩晓南, 王康东. 综合物探方法在采空塌陷区地质灾害勘查中的应用[J]. 工程地球物理学报,2013,10(5):725 − 729. [WANG Bin, HAN Xiaonan, WANG Kangdong. The application of integrated geophysical methods to investigating geological hazards in mined-out collapse areas[J]. Chinese Journal of Engineering Geophysics,2013,10(5):725 − 729. (in Chinese with English abstract) DOI: 10.3969/j.issn.1672-7940.2013.05.025 [13] 房浩, 李巧灵, 雷晓东, 等. 平原区深层隐伏岩溶塌陷主控因子地球物理调查方法适用性分析[J]. 水文地质工程地质,2020,47(1):153 − 160. [FANG Hao, LI Qiaoling, LEI Xiaodong, et al. Applicability of geophysical survey methods for the main controlling factors of deep covered karst collapse in plain areas[J]. Hydrogeology and Engineering Geology,2020,47(1):153 − 160. (in Chinese with English abstract) [14] 张淑坤, 张向东, 李永靖. 高速公路下伏老采空区探测及验证技术[J]. 中南大学学报(自然科学版),2015,46(9):3361 − 3367. [ZHANG Shukun, ZHANG Xiangdong, LI Yongjing. Old goaf detection and verification techniques under highway[J]. Journal of Central South University (Science and Technology),2015,46(9):3361 − 3367. (in Chinese with English abstract) DOI: 10.11817/j.issn.1672-7207.2015.09.028 [15] 张扬, 王璐, 刘智慧, 等. 多物探方法在胶东某金矿采空区勘查中的应用研究[J]. 地质与勘探,2019,55(3):809 − 817. [ZHANG Yang, WANG Lu, LIU Zhihui, et al. Application of multi-geophysical methods in exploration of a worked-out area of gold mining in the Jiaodong region[J]. Geology and Exploration,2019,55(3):809 − 817. (in Chinese with English abstract) -
期刊类型引用(0)
其他类型引用(5)