柳江盆地作为河北省重要的煤炭资源供应基地,同时作为华北板块及河北省内板块构造的缩影,其构造演化、岩浆侵入过程对华北克拉通破坏的研究具有典型性和重要意义[1-2]。盆地赋存石炭二叠纪和侏罗纪2个煤田,属于典型的双煤系煤田,也代表了华北煤田的2种典型类型煤田。其中的石炭二叠纪煤田是盆地内的主力煤田,其内赋存的优质无烟煤是周边市县的主要煤炭供应基地。侏罗纪以来盆地经历了多个期次的岩浆活动,在盆地中部形成了巨厚层、大面积分布的中性安山岩岩体,同时在岩浆的喷发过程中岩浆熔岩沿着煤层内、煤层其顶底板围岩软弱界面侵入,因此,在煤层及其围岩内形成了一系列的岩浆岩侵入体。经过多年的野外地质调查、地质填图、遥感地质解译、矿井生产地质调查、煤田地质勘查、钻孔岩芯精细描述以及煤层顶底板岩石力学测试等发现:相同岩性的煤层顶板岩石力学性质存在巨大的差异,煤层顶板的岩性组合类型在平面上分布具规律性,且盆地内火成岩的分布尤其是古火山口位置的分布与具有耦合性。
国内外对柳江盆地的研究最早可追溯到20世纪四五十年代,国内学者刘锡文、刘发等先后在盆地内开展了地质构造[3]、地层古生物[4-5]等方面的研究;近年来郝素琴、陈丹玲、李越等[6-8]对盆地内的火成岩、盆地周边火成岩岩性组合、盆地基底构造及其演化过程等进行了研究;马良、刘善德等[9-10]盆地内煤炭煤岩煤质受火成岩的影响以及盆地内的煤与天然焦的划分等进行了研究。此外,国内的学者[11-19]对我国不同煤田侵入煤层对煤岩煤质、孔隙结构、煤层瓦斯吸附性以及坚硬顶板致灾机理和防治措施等[20-21]进行了研究,且取得了理想的效果。但是以往研究多针对某一个矿井,而对盆地尺度的研究尤其是以柳江盆地具为代表且有典型性的盆地的研究较少,尤其缺少对盆地内岩浆的喷发、侵入对盆地内资源赋存及其开采技术条件的系统研究。
岩石力学特征是进行煤层顶、底板稳定性评价,矿井开拓设计和井巷围岩支护以及矿井采掘支护设备选型等的重要依据。研究表明,岩石力学特征决定煤层顶板的稳定性,而岩石力学特征受岩石岩性、成分、胶结物、填充物等多种因素的影响[22-23],后期的构造运动起着强化早期特征的作用[24-27]。而后期的构造运动尤其是煤系地层内岩浆侵入活动究竟对原生顶板力学性质有如何的影响呢?基于上述问题开展本次研究,以柳江盆地石炭二叠纪煤田主采煤层3煤层及其顶板为研究对象,从岩浆侵入煤层及其顶板的模式入手研究其对岩石力学特征的影响,通过建立盆地岩浆侵入煤层及其顶板的模式,对不同岩浆侵入模式下煤层顶板不同岩性、不同位置,与岩浆岩不同距离、厚度下的煤层及其顶板岩石进行系统的样品采集和室内测试分析,以期为柳江盆地内煤层顶板稳定性评价、矿井设计、井巷围岩采掘支护以及我国相同地质条件的煤田等提供参考。
柳江盆地位于河北省秦皇岛市海港区石门寨镇,大地构造单元属于塔里木中朝地台(Ⅰ)-河北板内块体(Ⅱ2)-燕山断褶带燕山南坡断隆区(图1),为燕山期构造运动形成的一个火山-构造盆地[3,9]。盆地总体平面形态为近椭圆形,南北长约17 km,东西宽约13 km,面积约221 km2,亦有学者认为盆地为走向南北向西开口的断陷盆地[10]。
图1 区域构造及其位置
Fig.1 Regional structure and its location
柳江盆地构造演化大致经历了基底形成、稳定盖层沉积和强烈活动3个阶段。盖层沉积阶段沉积了本区含煤地层太原组、山西组、下花园组;强烈活动阶段的中生代由于强烈的构造运动和岩浆活动,使得区内构造复杂化,形成不同类型的褶皱断裂构造,最终形成了盆地中西部近南北向东缓西陡的不对称向斜构造格局[8]。
盆地内的主力煤田为石炭二叠纪煤田,煤系太原组、山西组是主要层系。山西组的3煤层是主力煤层,太原组的4,5煤层为辅煤层。煤系厚度87~200 m,平均170 m;地层倾角5°~15°,平均10°;含煤层数3~5层,其中可采煤层3层,主力煤层3煤层厚度0~11.86 m,平均4.25m。煤层结构受侏罗纪岩浆侵入活动的影响导致煤层结构复杂、厚度变化大;煤类以无烟煤为主;3煤层直接顶板以闪长玢岩为主,局部发育粉砂岩、细粒砂岩、泥岩直接顶板(图2)。
图2 柳江盆地煤系地层综合柱状
Fig.2 Comprehensive histogram of coal rocks measures strata in Liujiang Basin
柳江盆地火成岩包括盆地周缘前寒武纪太古界花岗岩、燕山期花岗岩和盆地内中部大积分布的巨厚层状的中性安山岩体(图2)。其中对煤系影响大的侵入事件发生于侏罗世,即燕山运动第1期安山岩喷出及闪长岩侵入旋回,相当于燕山运动Ⅱ幕未至Ⅲ幕初的火山事件,地质年龄约165±10 Ma;岩浆侵入深度为500~1 200 m,其形成环境介于喷发与侵入岩之间,属次火山岩类[19]。
本次地质调查、钻孔揭露及以往研究[10]均表明侏罗纪髫髻山期火成岩的喷发通道为盆地中部的断裂构造,喷发方式为裂隙式间歇式喷发。钻孔及地面地质剖面测量发现火山的喷发序列大概经历了4~5个期次的喷发,在盆地内形成了由爆发相-喷出-溢流相的火山岩序列。在岩浆喷发过程中岩浆熔岩在煤系地层内顺煤层、构造、裂隙以及软弱面侵入、吞噬煤系从而形成了目前的煤层及其顶底板结构。煤系内侵入岩以浅灰白色、浅灰黑色安山玢岩、玄武安山玢岩和次闪长玢岩为主,显微镜下具有斑状结构,斑晶以斜长石为主含少量的角闪石和黑云母,基质与斑晶相似,细小板状斜长石散乱分布,其间充填暗色矿物和玻璃质;斑晶和基质均发生了强烈的碳酸盐化,呈变余碎屑假象(图3(f))。
盆地内岩浆岩在地表的盆地中部大面积分布在平面上岩性以安山岩、安山玄武岩(图3(e)~(f))、辉石安山岩、角闪安山岩(图3(d))、火山角砾岩、火山熔岩为主,局部可见火山弹(图3(h))、火山熔岩、枕状熔岩(图3(h))等;在纵向上发育火山集块岩(图3(i))、火山角砾岩(图3(a))、粗面安山岩(图3(c))、粗面岩、安山岩等;火成岩的厚度在平面上中部厚向四周变薄,地形高点和地层厚度极值点分布在老君顶、305高地、大石洼、265高地等地。经统计侵入岩的厚度总体上中部较厚、向四周变薄;各地层内侵入岩的厚度与地表喷出岩体的厚度呈正相关性。
按照文献[25-29]中提出的煤层直接顶板划分方法(表1),结合煤层直接顶板各类岩石力学测试结果,采用《煤矿区水文地质、工程地质及环境地质勘查评价标准》(MT/T 1091—2008)中煤层开采后冒落带的高度计算公式,同时参考盆地内现有国有大型煤矿采用的煤层开采工艺、工作面的长度、宽度等布置情况以及生产矿井揭露取4倍的煤层厚度作为垮落带高度即煤直接顶板厚度计算的参数。
图3 岩浆岩特征
Fig.3 Characteristics of magmatic
表1 煤层顶板类型划分方案[28-29]
Table 1 Coal seam roof type division scheme[28-29]
顶板类型相对位置岩性厚度阶段地勘阶段矿井阶段伪顶直接位于煤层之上泥岩、页岩、炭质泥岩≤0.50 m直接顶板直接位于伪顶或直接顶板1层或几层不同岩性的岩层等于垮落带的高度基本顶板直接位于直接顶板之上1层或几层不同岩性的岩层等于导水裂隙带的高度规范或经验给出的公式计算生产矿井或邻近生产实际测量值
按照上述方法结合本次勘查工程部署的62个钻孔,共对盆地内162个煤田地质勘查钻孔的煤层顶底板进行逐孔直接顶板岩性组合、厚度、类型划分和归纳,将盆地内主力煤层3煤层直接顶板类型归纳划分为5种类型(图4):
(1)吞噬型顶板。岩浆岩将煤层顶板炭质泥岩、泥岩或煤层全部吞噬煤层、岩浆沿着断裂通道向上喷出过程中的岩溶通道区域。该区域在盆地内主要为2种形式:① 岩浆岩吞噬煤层区域的无煤区(图4中的2205孔区域);② 3煤与3上煤层之间由于岩浆的顺煤层侵入形成的岩浆岩顶板及煤层分叉,岩浆的吞噬不仅完全改变了煤体及其特性,而且改造煤层宏观煤岩结构、显微煤岩组分(图5,图6(e)~(k))和煤质特征以及煤炭资源的加工利用[10-11]。
(2)侵入型顶板。岩浆岩沿着3煤顶部的泥岩、煤层侵入,但是3煤仍然具有一定的厚度,形成了煤层顶板岩性全部或大部分为侵入的闪长玢岩侵入岩体,厚度为煤层厚度的3~4倍的区域(图4中的3410孔区域)。从图5可以看出,岩浆侵入煤层后煤层的变质改造程度随着距离的增加而减弱。
(3)复合型顶板。在直接顶板范围内存在岩浆侵入体,形成了岩浆侵入体与其它岩性沉积岩的顶板岩性组合如图4中的3405孔区域(图6(a),(c))。
(4)强化改造型顶板。煤层直接顶板范围无岩浆侵入体,但是在其近距离范围内发育侵入岩,同时岩浆对顶板岩性及其力学性质发生了改造作用,强化并改造了顶板的力学性质,导致顶板岩石力学性质远优于原岩的岩石力学性质如图4的3402和3408孔区域(图6(d),(e))。
图4 柳江盆地石炭二叠纪煤田岩浆侵入及煤层直接顶板类型组合模式
Fig.4 Magma intrusion and coal seam direct roof type combination model of carboniferous-Permian coalfield in Liujiang Basin
图5 火成岩对煤层煤体结构的破坏显微煤岩特征(×500倍油侵镜下)
Fig.5 Maceral characteristic diagram of igneous rock intrusion coal seam coal body structure
(5)原生型顶板。岩浆的活动对直接顶板岩石力学性质基本无影响,因此,煤层顶板基本上保持了原始的岩石力学性质(图4的3411孔区域)。
在地面地质填图和矿井地质调查的基础上,根据古火山口的大概位置结合盆地构造复杂程度和煤层稳定程度,部署煤田地质勘查工程,详细开展地质岩芯钻探、岩芯精细描述、工程地质编录,并布设4条工程地质剖面,对剖面上的11个钻孔煤层顶板垮落带范围内的岩芯进行岩石力学样采集和剖面上的16个钻孔进行工程地质编录及工程地质测井(图1,2)。岩石力学样品根据岩浆侵入煤层顶板的模式分层系统采集,所有采样钻孔均采用绳索钻具、φ95 mm的金刚石钻头钻进,共采集175组、1 595块岩石力学样品,所有样品长度均大于120 mm,岩芯直径大于60 mm,在岩芯管提出井口后第一时间内进行采集,并及时密封防止岩芯风化、失水;根据试验结果分析不同岩性的岩石力学参数特征与侵入岩的厚度、与侵入岩的距离的关系。
2.3.1 煤层及其顶板总体特征
根据测试结果重点分析其中的149组煤层顶板岩石样品,对抗压强度、抗剪切强度、抗拉强度、岩石密度、孔隙度、含水率、杨氏模量、泊松比、RQD等关键力学参数进行测试(表2),结果发现:
(1)盆地内石炭二叠纪煤田煤层及其顶底板围岩岩石属于较软弱~坚硬岩类,以较坚硬岩类为主;其中泥岩、粉砂质泥岩、煤层属于较软弱岩类,粉砂岩属于较坚硬岩类,细粒砂岩、中粗粒砂岩属于坚硬岩类;侵入岩闪长玢岩、安山玢岩属于坚硬岩类,岩体完整,岩体质量好;各类顶板岩石完整性劣~完整,以较完整为主。
(2)盆地内煤层受到岩浆侵入活动的影响,导致煤层的岩石力学特征发生了改变,表现为抗压强度的增强,密度变大,抗剪强度增大等。
(3)各类岩石的抗压强度、抗剪切强度、RQD等关键岩石力学参数具有煤层<泥岩类<砂岩类<火成岩类的特征,其中砂岩类岩石具有随着砂岩粒度增加力学性质增强的特征;而岩石的含水率则相反,泊松比、孔隙度相差不大。
图6 火成岩侵入对煤系岩石的影响
Fig.6 Effect of igneous rock intrusion on coal measure rocks
表2 煤层顶板岩石力学总体特征
Table 2 Overall rock mechanics characteristics of coal seam roof
岩层(组数)自然抗压强度Rn/MPa饱和抗压强度Rc/MPa内摩擦角φ/(°)抗剪强度T/MPa弹性模量E/GPaRQD/%泊松比μ含水率/%岩石密度ρ/(g·cm-3)岩石孔隙度φe/%煤(5)31.35~68.5050.3420.74~34.3225.9735.81~36.9736.394.17~5.224.86—18~9746——1.75~1.881.81—泥岩(35)27.59~126.4867.0014.30~80.0937.5634.99~40.1137.724.02~13.507.351.44~5.813.1337~90660.20~0.240.220.15~1.443.132.58~2.912.741.03~8.223.99粉砂质泥岩(26)35.51~68.2250.8813.01~47.5127.9135.81~38.0636.821.87~3.862.791.73~3.252.4342~97780.20~0.260.230.40~1.030.652.58~2.782.700.72~4.552.23粉砂岩(29)28.39~121.5074.2717.02~82.3744.5635.10~39.7337.413.93~13.307.591.43~4.973.2373~96870.17~0.260.220.18~1.030.512.68~2.792.760.72~3.972.22细粒砂岩(30)52.07~113.7088.5130.97~94.0758.2836.48~40.2638.596.44~12.108.763.21~5.404.2873~99920.17~0.260.210.12~1.430.482.67~2.762.720.72~5.342.13中、粗粒砂岩(7)42.46~147.5891.4822.67~106.0360.9136.05~39.3237.965.86~12.778.861.84~6.353.9885~100940.19~0.260.220.12~1.020.382.55~2.742.651.09~7.473.57闪长玢岩(17)79.06~163.23105.5456.17~117.2074.4538.60~39.6738.906.89~13.779.783.73~6.784.9387~100960.16~0.200.180.04~1.070.552.69~2.752.722.52~4.583.38汇总(149)27.59~163.2379.0013.01~117.2050.1134.99~40.2637.861.87~13.778.041.43~6.783.7218~100800.16~0.260.210.04~1.440.521.75~2.912.590.72~8.222.65
注:数据格式下同。
2.3.2 不同侵入模式下岩石力学特征
在煤层顶板岩浆岩侵入煤层顶板类型划分和岩石力学特征总体分类统计、分析的基础上,按照岩性、距离侵入岩的距离、侵入岩的厚度等分类统计,分析岩石力学各参数与上述各因素之间的相关关系(表3,4)。
表3 岩浆对煤层及其顶板不同岩性岩石力学参数的影响
Table 3 Influence of magma on rock mechanical parameters of different lithology in coal seam and its roof
岩层(组数)侵入岩体距离/m侵入岩厚度/m自然抗压强度Rn/MPa饱和抗压强度Rc/MPa内摩擦角φ/(°)抗剪强度T/MPa弹性模量E/GPa泊松比μ含水率w/%岩石密度ρ/(g·cm-3)岩石孔隙度φe/%煤(5)0.051.0068.5034.3237.255.18———1.78~1.901.85—0.104.5051.1722.8436.975.22—————0.661.2231.3520.7435.814.17—————3.250.8728.7215.2334.353.98—————原煤(1)28.7215.2334.353.981.75~1.881.81侵入(3)31.35~68.5050.3420.74~34.3225.9734.35~37.2536.094.17~5.224.86吞噬(1)79.0666.2439.576.873.730.160.042.763.17泥岩(35)04.3069.4535.23—9.34—————0.691.2865.0631.1338.326.893.100.200.462.835.081.323.13126.4880.0940.1110.075.810.220.272.761.083.2011.4927.5914.3034.994.021.440.240.942.718.225.151.5552.6937.6136.306.402.860.201.442.673.99原岩27.5914.3034.994.021.440.240.942.718.22岩浆侵入后52.69~126.4876.6331.13~80.0849.6136.30~40.1138.245.18~13.508.412.86~5.813.920.20~0.240.220.27~1.440.722.67~2.832.751.08~5.083.38粉砂质泥岩(26)0.106.4576.3439.0338.567.543.250.240.152.911.030.265.5068.2247.5135.819.743.250.200.402.750.726.415.5048.9123.2038.065.182.320.221.032.584.5557.961.5535.5113.0136.584.801.730.260.532.781.43原岩35.5113.0136.584.801.730.260.532.781.43岩浆侵入后48.91~76.3464.4923.20~47.5136.5835.81~38.5637.485.18~9.747.492.32~3.252.940.20~0.240.220.15~1.030.532.58~2.912.750.75~4.552.10粉砂岩(29)012.10~12.3812.2498.80~121.50110.15——10.40~13.3011.85—————0.47~3.511.850.70~11.493.8828.39~101.7563.8717.02~72.7539.4135.10~39.7337.363.95~10.686.471.43~4.833.040.17~0.260.230.18~1.030.532.68~2.952.750.72~3.972.646.51~9.608.160.45~10.404.6451.12~118.8182.3625.78~101.7563.8735.32~39.6037.456.55~9.628.272.90~5.914.400.18~0.240.210.26~0.640.422.74~2.782.761.08~2.171.7119.4510.4089.4565.3035.804.461.630.230.832.780.7251.041.5538.1621.9639.168.494.280.190.412.792.47原岩28.0916.9535.803.851.400.251.032.652.93岩浆侵入后28.39~121.5074.2717.02~82.3744.5635.10~39.7337.413.95~13.307.591.43~4.973.230.17~0.260.230.18~1.030.512.68~2.792.760.72~3.972.22细粒砂岩(30)0(2)0.45~12.386.42105.10~113.70109.4094.0739.808.85~10.509.685.280.170.332.752.160.42~2.541.151.40~10.404.5752.07~108.5280.8330.97~81.2452.8536.48~39.1538.196.44~12.108.763.21~5.073.940.17~0.260.210.12~1.430.622.67~2.762.700.72~5.342.235.37~6.485.811.94~10.405.1669.83~108.0993.5640.42~75.9963.1038.48~38.7238.576.46~10.158.723.23~5.404.580.18~0.240.200.17~0.330.242.73~2.812.760.73~2.141.2014.0010.4092.9846.0240.267.094.680.230.332.744.23原岩52.0730.9736.486.443.210.241.432.671.85岩浆侵入后52.07~113.7088.5130.97~94.0759.9036.48~40.2638.596.44~12.108.763.21~5.404.280.17~0.260.210.12~1.430.482.67~2.762.720.72~5.342.13中、粗粒砂岩(7)0.140.4571.1549.5238.507.473.170.200.122.627.472.288.2094.4550.0038.508.453.950.250.182.631.877.750.7042.4622.6736.055.861.840.261.022.692.908.284.30101.7676.3537.449.754.590.200.372.554.5348.9824.10147.58106.0339.3212.776.350.190.212.741.09原岩42.4622.6736.055.861.840.261.022.692.90岩浆侵入后42.46~147.5891.4822.67~106.0360.9136.05~38.5037.965.86~12.778.861.84~6.353.980.19~0.260.220.12~1.020.382.65~2.812.731.09~7.473.57
表4 岩浆活动对岩石力学参数的影响总体特征
Table 4 Statistical table of overall characteristics of influence of magmatic activities on rock mechanical parameters
侵入岩体距离/m侵入岩厚度/m自然抗压强度Rn/MPa饱和抗压强度Rc/MPa内摩擦角φ/(°)抗剪强度T/MPa弹性模量E/GPa泊松比μ含水率/%岩石密度ρ/(g·cm-3)岩石孔隙度φe/%00.45~12.387.6568.50~121.5096.3335.23~94.0764.6539.805.280.170.332.752.160.10~3.511.460.45~11.494.1727.59~126.4871.4914.30~81.2444.4934.99~40.1137.673.95~12.107.381.43~5.813.510.17~0.260.220.12~1.430.522.62~2.832.730.72~8.222.885.15~9.607.060.45~10.404.1942.46~118.8177.8122.67~82.3750.4735.32~39.6037.585.18~10.157.861.84~5.403.670.18~0.260.210.17~1.440.572.55~2.782.700.73~4.552.4014.00~19.4516.850.45~10.407.0876.34~92.9886.2539.03~65.3050.1238.56~40.2639.337.09~8.497.713.25~4.684.070.19~0.240.220.15~0.410.302.74~2.912.811.03~4.232.5824.55~48.9536.762.20~24.1013.15107.41~147.58127.5073.86~106.0389.9539.32~39.5439.438.04~12.7710.414.66~6.355.510.19~0.210.200.21~0.250.232.68~2.742.711.09~1.111.1051.04~57.9654.501.5535.51~38.1636.8413.01~21.9619.4935.80~36.5836.194.46~4.804.631.63~1.731.680.23~0.260.250.53~0.830.682.780.72~1.431.08原岩27.59~52.0735.7413.01~30.9718.8635.35~36.5835.713.85~6.444.831.40~3.211.920.24~0.260.250.53~1.430.991.75-2.782.551.43-8.223.47侵入27.59~147.5876.7513.01~106.0347.7434.99~40.2637.763.95~13.507.901.43~6.353.600.17~0.260.220.12~1.440.512.25~2.992.730.72~8.222.58吞噬79.0666.2439.576.873.730.160.042.763.17
图7 岩浆侵入厚度、侵入距离与岩石力学参数关系
Fig.7 Relationship diagram of magma intrusion thickness,distance and rock mechanical parameters
(1)盆地内煤层及其直接顶板岩石的抗压强度、抗剪强度、弹性模量与岩浆侵入岩层的距离呈负相关性,且为指数关系(图7(e));而与岩石的密度、孔隙度、含水率以及泊松比等呈正相关性,且为对数关系(图7(f));但是各类岩石力学参数与岩层侵入岩层的厚度相关性不明显,可能是多期次、长时间的岩浆活动作用相叠加而掩盖了二者之间的关系。
(2)通过统计盆地内地表喷出中性安山岩体的厚度与煤系地层内侵入岩闪长玢岩、安山玢岩等岩体的在同一个钻孔上的厚度,发现二者具有正相关性(图7(a)),侧面印证了侵入岩与喷出岩是同期的产物、盆地中部的岩浆岩为原地喷发的观点[6]。
(3)在与煤层顶、底板岩浆侵入体距离相同的情况下,各类岩石的抗压强度、抗剪切强度、杨氏模量力学参数与岩浆侵入岩体的厚度呈正相关性,且为线性关系(图7(c),8);而与内摩擦角、泊松比、岩石的密度等相关不明显(图7(b),8);而与岩石的含水率、孔隙度呈负相关性,其中的含水率表现为幂函数相关关系、孔隙度表现为线性相关关系(图7(d),8)。
(4)岩浆活动对煤层及其顶板侵入程度不同对岩石力学参数的影响表现为:抗压强度、内摩擦角、岩石的密度、弹性模量吞噬区>侵入区>原岩区;含水率、泊松比、孔隙度等参数均受岩浆侵入影响均减少,表现为孔隙度、泊松比原岩区>侵入区>吞噬区(图8(g)~(i),8(j)~(l),9)。
根据上述统计分析所得的结论,对盆地内主力煤层3煤进行顶板稳定性评价,根据传统的煤层直接顶板划分方法和前述划分结果将顶板按照吞噬型、侵入型、复合型、改造型、原生型5种岩性组合进行分类统计,最后按照两种划分方法分别进行煤层顶底板稳定性评价,并对比评价结果、分析方法的先进性。
评价方法统一采用RQD,Z,M,BQ值评价指标法结合煤层顶板岩性组合、顶板面积占比、顶板厚度等进行综合评价,评价结果见表5。
采用《煤矿区水文地质、工程地质及环境地质勘查评价标准》(MT/T1091—2008)、《矿区水文地质工程地质勘探规范》(GB12719—91)和《岩土工程勘察规范》(GB50021—2001)中提供的RQD值法(表6)、岩体质量系数法(表7、式(1))和岩体质量指标法(表8、式(2))以及岩体基本质量评价法(式(3))对比评价个岩层岩体质量及其完整性,评价方法见表7~9和式(1)~(3)。
RQD值评价法:岩石质量等级和岩体完整性评价。
岩体质量系数法Z值为
Z=IfS
(1)
岩体质量指标法M值为
(2)
式中,Z为岩体质量系数;I为岩体完整系数(以RQD值代替);f为结构面摩擦因数;S为岩块坚硬系数(S=Rc/10);M为岩体质量指标;Rc为岩石饱和抗压强度,MPa。
BQ值评价法:
BQ=100+3Rc+250Kv
(3)
式中,Kv为岩体完整性系数(以RQD值代替)。
按照上述评价方法,采用揭露点面积占比、平面上的分布的连续性、Z值、M值、BQ值等指标,根据3煤层直接顶板的类型分区分类进行顶板稳定性评价,评价结果表明(表3、图4):
(1)煤层顶板的岩石力学性质受岩浆侵入的影响和控制,形成了顶板类型平面上环条带状分布的特征,其中以侵入型为主,复合型、强化型、原生型次之(表5);相应的顶板稳定性评价的结果亦具有这一特征(图10(a),(b)),2种方法评价的结果均显示了这一规律,但是新的评价方法表现的这种规律性更强、且更为全面,同时也符合地质规律。
图8 岩浆侵入距离、厚度与岩石力学特征关系
Fig.8 Relationship between magma intrusion distance,hickness and rock mechanics characteristics
图9 煤层顶板岩石力学特征总体特征
Fig.9 General characteristics of rock mechanics characteristics of coal seam roof
表5 煤层顶板稳定性评价成果
Table 5 Evaluation results of coal seam roof stability
顶板类型揭露点数面积占比/%Rc/MPa岩体完整性Z值M值BQ值基本质量级别顶板类型存在的工程地质问题吞噬型92105.54完整7.823.34562Ⅱ难冒落(Ⅲ类)坚硬顶板长期悬顶,侵入型523676.75完整5.152.35520Ⅱ难冒落(Ⅲ类)顶板坚硬长期悬顶,需要人工放顶复合型241860.62较完整3.941.78516Ⅲ中等冒落(Ⅱ类)坚硬、软弱岩层相间发育,顶板冒落事故频发地段强化改造型322044.56较完整2.831.29447Ⅲ中等—易于冒落(Ⅱ类)介于复合顶板与原岩型之间原岩型452435.74较完整2.181.01404Ⅲ易冒落(Ⅰ类)顶板较坚硬,易于冒落
表6 岩石坚硬程度分级
Table 6 Rock hardness grading table
岩石坚硬程度极软岩软岩较软岩较坚硬岩坚硬岩单轴饱和抗压强度/MPaRc≤515≥Rc>530≥Rc>1560≥Rc>30Rc>60
表7 岩体Z值范围及其优劣分级
Table 7 Rock mass Z value range and its advantages and disadvantages classification table
岩体质量系数Z<0.10.1~0.30.3~2.52.5~4.5>4.5岩体质量等级极坏坏一般好特好
表8 岩体质量分级
Table 8 Rock mass quality classification table
岩体分类ⅤⅣⅢⅡⅠ岩体质量指标M<0.010.01~0.120.12~1.01.0~3.0>3岩体质量坏差中等良优
表9 岩石质量等级
Table 9 Rock quality grade table
岩石质量等级ⅤⅣⅢⅡⅠRQD值/%<2525~5050~7575~9090~100岩石质量极劣的劣的中等的好的极好的岩体完整性评价破碎差中等较完整完整
图10 煤层顶板稳定性分区评价及其对比
Fig.10 Zoning evaluation chart of roof stability of coal seam
(2)按照传统的直接顶板岩性进行评价,研究区内顶板以侵入岩为主,掩盖了顶板岩性组合类型、岩浆侵入对顶板岩石力学特征的影响且忽略了开采因素、开采引发的顶板工程地质问题(图10(a));而按照新的顶板类型划分方法与岩浆侵入岩层对岩石力学特征的影响程度相结合,进行岩浆侵入区煤层顶板类型划分和稳定性评价,进而揭露潜在的工程地质问题(图10(b))。
(3)盆地内主力煤层顶板大面积受岩浆活动的影响,岩浆对煤层顶板岩石力学性质发生了不同程度的改造作用,总体上强化了岩石力学参数和岩石的完整性;导致盆地内受岩浆活动影响范围内煤层顶板的稳定性评价为中等~难冒落顶板,而原生型顶板区域为易于冒落顶板。
(4)盆地内的生产矿井对于吞噬区的探查、需要采取有效的措施,该区是未来矿井生产能否正常生产的关键,同时可以利用岩浆吞噬区作为矿井设计关键硐室、井筒落地的层位选择区域;侵入型顶板区域顶板坚硬,存在长期悬顶,需要采取人工强制放顶而避免一次性大面积跨落造成矿震动和应力集中形成的矿井冲击地压显现;而对于复合型顶板区域,薄层的坚硬的侵入岩浆岩与软弱的砂、泥岩互层层面是软弱面,易于形成顶板事故,是顶板管理的重点区域;强化区域和原岩区域顶板易冒落。
(5)盆地内的煤层及其直接顶板岩石力学性质受岩浆侵入活动的影响,尤其是对煤层,其本质是煤体结构的改造,导致煤层变为较坚硬-坚硬岩体,因此,在巷道掘进、工作面回采及其顶板支护设备的参数选择中需要重点考虑;同时,吞噬和煤层内的岩浆侵入体将严重影响采掘的进度,首先需要采取有效的技术措施查明吞噬区分布位置、火成岩夹矸在工作面的空间展布特征,然后采区有效的措施进行超前设计和弱化处理。
(1)建立了柳江盆地内侏罗纪岩浆活动对盆地内煤层及其顶板岩石力学特征的侵入模式,盆地侏罗纪经历了多期次喷发和侵入活动,对煤层顶板发生了改造作用,并将煤层顶板划分为5种类型,即顶板泥岩、煤层的吞噬型、侵入型、复合型和强化改造型以及原生型。
(2)岩浆的侵入、高温烘烤、吞噬及侵入和喷发过程中的挤压作用等综合效应强化了煤层及其顶底板岩石的力学和工程地质参数特征,表现为岩石密度增加,孔隙度、含水率等减少;抗压强度、抗剪强度、岩体的完整性等关键力学参数增大,且与侵入岩体的厚度、距离等具有一定的相关性。
(3)在盆地内圈定了5个岩浆吞噬区,吞噬区位于古火山口或岩浆喷发通道相上;以古火山口或火山通道相形成的吞噬区为中心煤层顶板类型向外依次渐变为侵入型、复合型、强化改造型和原生型。吞噬型顶板区可以作为矿井关键硐室、井筒落底的理想选择位置;复合型顶板区域是顶板事故的多发区,是顶板管理的重点区域;侵入型顶板区域在工作面回采过程中存在坚硬顶板长期悬的安全开采地质问题,需要采用人工强制放顶措施,该区也是矿井冲击地压管理的重点区域;强化改造型、原生型顶板区域属于中等-易于冒落顶板。
(4)柳江盆地作为华北板块的典型和缩影,上述评价方法及其岩浆活动对煤层顶板岩石力学特征的影响研究结果不仅可以指导盆地内矿井设计、顶板稳定性评价、顶板支护管理以及开采设备的选型;而且可以为华北煤田具有同或相似地质条件的矿井提供参考和研究思路。
(5)揭示了柳江盆地内侏罗纪岩浆侵入活动对煤层及其顶板岩石力学各参数的影响程度,岩石力学参数与岩浆岩侵入层厚度、距离之间的关系;揭露了盆地内主力煤层顶板类型的环带状分布的特征,同时顶板稳定性评价结果与其具有一致性。
(6)通过对比分析文中采用的顶板类型划分方法与传统的顶板类型划分、评价的结果,说明文中采用的方法可以揭示矿井开采过程中顶板潜在的工程地质问题,适合于岩浆岩侵入煤层类的顶板类型划分和顶板稳定性评价。
(7)盆地内的矿井设备选型中一定要考虑岩浆侵入对煤及其顶板岩石力学特征的强化作用,为了保证矿井正常开采,在工作面回采前需要对工作面内火成岩夹矸、吞噬体及其改造后的坚硬煤体首先需要进行有效的探查,查明吞噬体的规模、火成岩夹矸的空间展布,然后采区有效的技术措施如吞噬区的预先采掘,厚层岩浆岩夹矸的水力致裂、爆破等超前弱化预处理,以保证矿井的安全、正常生产。
[1] 朱日祥,徐义刚.西太平洋板块俯冲与华北克拉通破坏[J].中国科学:地球科学,2019,49(5):836-847.
ZHU Rixiang,XU Yigang.The subduction of the west Pacific plate and the destruction of the North China Craton[J].Science China Earth Sciences,2019,49,49(5):836-847.
[2] 钟福平,钟建华,艾合买提江·阿不都热合曼,等.华北克拉通破坏时间与破坏范围分布特征——来自银根—额济纳旗盆地苏红图坳陷早白垩世火山岩的启示[J].中国地质,2015,42(2):435-456.
ZHONG Fuping,ZHONG Jianhua,AHMATJAN Abdurahman,et al.Timing and scale of the destruction of the north china craton:Revelation from the early cretaceous volcanic rocks in Suhongtu depression of Inggen-Ejin Banner Basin[J].Geology in China,2015,42(2):435-456.
[3] 刘锡文.从秦皇岛柳江盆地岩浆活动论其大地构造性质及成矿特征[J].唐山工程技术学院学报,1995(1):76-79.
LIU Xiwen.On the geotectonic properties and metalorganic characteristics of Liujiang Basin in Qinhuangdao[J].Journal of Tangshan Institute of Engineering and Technology,1995(1):76-79.
[4] 刘发.河北秦皇岛石门寨地区晚石炭世本溪组下部的腕足动物及其古环境[J].古生物学报,1992(2):249-255,283-284.
LIU Fa.Brachiopods and paleo-environment of the lower part of the Late Carboniferous Benxi Formation in Shimenzhai area,Qinhuangdao,Hebei[J].Journal of Paleontology,1992(2):249-255,283-284.
[5] 冯诗海,李红,蒋佳俊,等.华北秦皇岛地区柳江盆地马家沟组碳酸盐岩多期白云化作用[J].沉积学报,2017,35(4):664-680.
FENG Shihai,LI Hong,JIANG Jiajun,et al.Multiphasic dolomitization of carbonate rocks in Majiagou Formation,Liujiang Basin,Qinhuangdao Area,North China[J].Journal of Sedimentation,2017,35(4):664-680.
[6] 郝素琴,曲以秀,杨洪英.柳江火山—构造盆地的特征及演化[J].沈阳黄金学院学报,1997(1):16-22.
HAO Suqin,QU Yixiu,YANG Hongying,Liujiang volcanic-tectonic basin characteristics and evolution[J].Journal of Shenyang Gold Institute,1997(1):16-22.
[7] 陈丹玲,赖绍聪,刘养杰.秦皇岛柳江盆地混合花岗岩的锆石U-Pb定年[J].西北大学学报(自然科学版),2007(2):277-281.
CHEN Danling,LAI Shaocong,LIU Yangjie.Zircon U-Pb dating of mixed granites in Liujiang Basin,Qinhuangdao[J].Journal of Northwest University(Natural Science Edition),2007(2):277-281.
[8] 马良,舒建生,王海军,等.抚宁秦能一矿整合区煤炭资源储量核实报告(2017)[R].西安:中煤科工集团西安研究院有限公司,2018:76-93.
[9] 李越,季建清,涂继耀,等.燕山东部柳江地区构造属性新解与郯庐断裂系活动[J].岩石学报,2009,25(3):675-681.
LI Yue,JI Jianqing,TU Jiyao,et al.New solution of structural attributes and Tanlu fault system activity in Liujiang area,Eastern Yanshan[J].Journal of Petrology,2009,25(3):675-681.
[10] 马良.柳江盆地内岩浆侵入活动对煤层煤质的影响[J].煤炭科学技术,2019,47(8):226-234.
MA Liang.Effect of magmatic intrusion activity on coal bed and quality in Liujiang Basin[J].Coal Science and Technology,2019,47(8):226-234.
[11] 刘善德,等.柳地勘时期无烟煤与天然焦的煤岩学识别[J].中国煤炭地质,2019,31(8):1-7.
LIU Shande.Coal petrology recognition of anthracite and natural coke during geological survey[J].Coal Geology of China,2019,31(8):1-7.
[12] 闫浩,张吉雄,张强,等.巨厚火成岩下采动覆岩应力场-裂隙场耦合演化机制[J].煤炭学报,2016,41(9):2173-2179.
YAN Hao,ZHANG Jixiong,ZHANG Qiang,et al.Coupled evolution mechanism of stress field and fissure field in overlying strata under huge igneous rocks[J].Journal of China Coal Society,2016,41(9):2173-2179.
[13] 赵志根,唐修义,黄文辉,等.祁东煤矿岩浆侵入特征的研究[J].中国煤炭地质,2009,21(1):27-29.
ZHAO Zhigen,TANG Xiuyi,HUANG Wenhui,et al.Research on characters of magmatic intrusion in Qidong Coal mine[J].Coal Geology of China,2009,21(1):27-29.
[14] 蒋静宇,程远平.淮北矿区岩浆岩侵入对煤储层微孔隙特征的影响[J].煤炭学报,2012,37(4):634-640.
JIANG Jingyu,CHENG Yuanping.Influence of igneous intrusion on micropore characteristics of coal reserves in Huaibei coal field[J].Journal of China Coal Society,2012,37(4):634-640.
[15] 蒋静宇,程远平,王海锋,等.岩浆侵入对煤吸附瓦斯特性的影响分析[J].采矿与安全工程学报,2012,29(1):118-123.
JIANG Jingyu,CHENG Yuanping,WANG Haifeng,et al.Effect of igneous intrusion on adsorption characteristics of coal to gas[J].Journal of Mining & Safety Engineering,2012,29(1):118-123.
[16] 卢平,鲍杰,沈兆武.岩浆侵蚀区煤层孔隙结构特征及其对瓦斯赋存之影响分析[J].中国安全科学学报,2001,11(6):41-43.
LU Ping,BAO Jie,SHEN Zhaowu.Pore structure characteristics of coal seam in magmatic rocke rosive area and its effect on gas storage[J].China Safety Science Journal,2001,11(6):41-43.
[17] 李恒乐,张玉贵,秦勇,等.淮北矿区祁东井田中生代岩浆侵入作用对煤层瓦斯赋存的影响[J].煤炭学报,2013,38(11):1982-1987.
LI Hengle,ZHANG Yugui,QIN Yong,et al.Influence of Mesozoic magmatic intrusion on coalbed gas occurrence in Qidong well field,Huaibei mining area e[J].Journal of China Coal Society,2013,38(11):1982-1987.
[18] 李胜利,郝俊杰,徐颖,等.冀南地区燕山期岩浆活动与煤的变质作用研究[J].中国煤炭地质,2017,29(4):1-7.
LI Shengli,HAO Junjie,XU Ying,et al.Yanshanian magmatism and coal metamorphism in southern Hebei[J].Coal geology of China,2017,29(4):1-7.
[19] 杨培举,何烨,郭卫彬.采场上覆巨厚坚硬岩浆岩致灾机理与防控措施[J].煤炭学报,2013,38(12):2106-2112.
YANG Peiju,HE Ye,GUO Weibin.Disaster-causing mechanism and prevention and control measures of super-thick hard magmatic rocks above the stope[J].Journal of China Coal Society,2013,38(12):2106-2112.
[20] 王海军.盘关向斜火烧铺井田瓦斯赋存地质特征研究[J].煤炭科学技术,2019,47(12):90-101.
WANG Haijun.Research on geological characteristics of gas occurrence in Huoshaopu minefield of Panguan syncline[J].Coal Science and Technology,2019,47(12):90-101.
[21] SMITH Jack A.Analysis of the influence of groundwater and the stress regime on bolt behaviour in underground coal mines[J].International Journal of Coal Science & Technology,2019,6(2):286-300.
[22] 孟召平,彭苏萍.煤系泥岩组分特征及其对岩石力学性质的影响[J].煤田地质与勘探,2004,32(2):14-16.
MENG Zhaoping,PENG Suping.Mudstone composition of coal measures and its influence on the mechanical properties[J].Coal Geology & Exploration,2004,32(2):14-16.
[23] 孟召平,彭苏萍,傅继彤.含煤岩系岩石力学性质控制因素探讨[J].岩石力学与工程学报,2002,21(1):102-106.
MENG Zhaoping,PENG Suping,FU Jitong.Study on control factors of rock mechanics properties of coal bearing formation[J].Chinese Journal of Rock Mechanics and Engineering,2002,21(1):102-106.
[24] 蒋景东,陈生水,徐婕,等.不同含水状态下泥岩的力学性质及能量特征[J].煤炭学报,2018,43(8):2217-2224.
JIANG Jingdong,CHEN Shengshui,XU Jie,et al.Mechanical properties and energy characteristics of mudstone under different containing moisture states[J].Journal of China Coal Society,2018,43(8):2217-2224.
[25] 王海军.永陇矿区LYL井田2号煤层顶板沉积环境及其稳定性评价[J].煤田地质与勘探,2016,44(6):38-44.
WANG Haijun.Sedimentary environment and stability evaluation of seam roof of No.2 seam in LYL coal mine of Yonglong coal field[J].Coal Geology & Exploration,2016,44(6):38-44.
[26] 王海军.煤层顶板沉积环境对其稳定性影响研究[J].煤炭科学技术,2017,45(2):178-184.
WANG Haijun.Study on impact of sedimentary environment to stability of seam roof[J].Coal Science and Technology,2017,45(2):178-184.
[27] 王海军,马良.陕北侏罗纪煤田三角洲平原沉积环境及其岩石力学特征[J].煤田地质与勘探,2019,47(3):61-69.
WANG Haijun,MA Liang.Study on sediment environment and rock mechanics characteristics of the delta plain of Jurassic coalfield in northern Shaanxi[J].Coal Geology & Exploration,2019,47(3):61-69.
[28] 王海军.关于煤层顶板划分中岩性及厚度取值的探讨[J].能源与环保,2017,39(1):86-89.
WANG Haijun.Discuss on division of lithology and thickness values of coal roof[J].China Energy and Environmental Protection,2017,39(1):86-89.
[29] 王海军,马良,朱玉英.煤层顶板类型划分新方法及其稳定性评价[J].煤矿安全,2020,50(4):207-212.
WANG Haijun,MA Liang,ZHU Yuying.A new method of coal seam roof classification and its stability evaluation[J].Safety in Coal Mines,2020,50(4):207-212.