资源描述:
第 卷第 期煤 炭 科 学 技 术 年 月 移动扫码阅读 王 亮,廖晓雪,褚 鹏,等瓦斯抽采穿层钻孔钻扩造穴卸压增透机理研究煤炭科学技术,,() , , , , , ( ) 瓦斯抽采穿层钻孔钻扩造穴卸压增透机理研究 王 亮,廖晓雪,褚 鹏,张晓磊,刘清泉 (中国矿业大学 安全工程学院,江苏 徐州 ;常州大学 环境与安全工程学院,江苏 常州 ) 摘 要随着我国煤矿开采逐渐向深部开拓,煤层瓦斯压力增大、含量增加,煤层渗透率普遍较低,增 加了瓦斯抽采的难度。 钻孔钻扩造穴卸压增透技术能提高煤层渗透率,是增加瓦斯抽采效率的有效 手段。 为解决寺家庄煤矿 号煤层渗透率低,瓦斯难抽采的问题, 以寺家庄煤矿北翼辅助运输巷 号煤层抽采钻孔造穴为工程背景,以弹性力学应变软化模型和扩散渗流煤层瓦斯流动理论模型为 基础,建立了穿层钻孔钻扩造穴后煤层渗透率演化方程和穿层钻孔造穴煤层瓦斯流动方程。 通过 多物理场数值模拟软件对方程进行解算,分析了钻孔造穴的增透机理和渗透率 分布规律,得到了瓦斯抽采量、瓦斯压力分布以及渗透率等关键参数;结合模拟结论与现场条件,确定 了钻孔钻扩造穴的最优造穴半径为 ,最佳布孔间距为 ,为现场施工造穴半径和钻孔 间距的确定提供了指导。 最后,在寺家庄煤矿 号煤层实施普钻钻孔与造穴钻孔,进行了瓦斯抽采 效果对比。 试验结果表明实施水力钻扩造穴技术后,瓦斯抽采率、抽采浓度和抽采纯量相对普钻钻 孔分别提高了约 、、 倍;瓦斯抽采周期、钻孔施工工程量降低了约 倍,提高了瓦斯抽采技术 经济效益,应用结果验证模拟结果可靠,抽采设计可行,可以指导现场施工。 关键词穿层钻孔;钻扩造穴增透;流固耦合; 瓦斯抽采 中图分类号 文献标志码 文章编号() , , , , ( , , ,; , , ,) 收稿日期;责任编辑曾康生 基金项目国家自然科学基金面上资助项目();中央高校基本科研业务资助项目();江苏省自然科学基金青年基金资助 项目();中国博士后科学基金批面上资助项目() 作者简介王 亮(),男,江苏东海人,教授,博士生导师,博士。 通讯作者张晓磊(),男,安徽淮北人,讲师,博士。 , , , , , , , , , , , , 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 年第 期煤 炭 科 学 技 术第 卷 , , , , , ; ; ; 引 言 煤炭是我国的主导能源, 年在我国能源的 消费比例约占 。 煤层瓦斯作为煤的伴生产 物,是储量丰富的优质清洁能源,同时也是威胁煤矿 安全生产的主要灾害之一。 煤与瓦斯突出和瓦斯爆 炸事故常造成重大的人员伤亡和经济损失,而瓦斯 抽采是消除瓦斯事故、保障安全开采的主要方法,同 时能够使瓦斯得到充分利用。 煤层渗透率是决 定瓦斯抽采质量和产量的关键因素,然而,由于 我国煤层地质条件复杂,煤体渗透率普遍偏低,通常 在 以下,导致瓦斯抽采难度大,抽采效 率低。 为提高煤层渗透率,促进瓦斯高效抽采,有学者 提出水力压裂、水力割缝等技术措施,通过水的高压 作用在周围煤体产生人工裂隙,促使煤体内部瓦斯 渗流,达到增透的目的。 然而,水力割缝对软煤 的持久性差,在地应力大的情况下裂缝容易闭合;水 力压裂易导致煤体局部应力集中,影响增透效果,两 者对松软构造煤的瓦斯抽采效果较差。 钻孔钻扩造 穴扩大钻孔是近年来广泛应用的增透技术,钻扩造 穴是在钻刀切割煤体的过程中,同时打开高压水泵, 利用高压水射流动力软化、破碎钻孔周围的煤体,再 将破碎的煤体冲出钻孔,构建大尺寸孔洞,使周围煤 体卸压。 造穴过程中,钻孔周围煤体应力重新分布, 当煤体强度不足以承受重新分布的应力时,煤体就 会发生破坏。 研究表明,煤体破坏后,渗透率能 增加几十至几千倍。 钻孔钻扩造穴技术集钻 孔冲孔为一体,克服了水力压裂、水力割缝和传统 水力冲孔方法的不足,出煤量大,增透有效。 根据瓦斯抽采钻孔形式,钻孔造穴有穿层钻孔 和顺层钻孔 种形式,穿层钻孔是从岩石巷道向煤 层打钻孔穿透煤层。 因穿层钻孔在岩石中比顺层钻 孔更加稳固,抽采服务时间更久。 近年来,水力钻扩 造穴技术,特别是实施穿层钻孔抽采瓦斯时的钻扩 造穴,在国内煤矿得到了广泛应用。 孙四清等对 松软突出煤层进行穿层钻孔造穴瓦斯抽采后,瓦斯 抽采量增加了 倍。 杜昌华等对大倾角松软 厚煤层实施穿层水力扩孔技术后,造穴钻孔瓦斯抽 采量和单孔瓦斯体积分数分别增加了 和 倍, 并且扩孔半径越大,煤中残余瓦斯含量越低。 牟全 斌等对芦岭煤矿 工作面实施了机械造穴 穿层钻孔后,单孔瓦斯抽采浓度和瓦斯抽采纯量相 应提高了 倍、 倍。 以上研究表 明,实施穿层钻孔造穴卸压增透后,煤层瓦斯抽采效 率显著提高。 为了获得最有效的增透效果,降低施 工成本,郝丛猛等根据不同煤层的特点,采用数 值模拟软件对造穴半径、布孔间距等参数进行了合 理设计,为现场工程提供了指导。 于宝种建立了 冲孔造穴过程中煤体的渗透率演化模型,模拟了不 同造穴半径对煤层渗透率的影响,并在新景煤矿进 行了工程验证,试验结果表明冲孔造穴能够减少钻 孔工作量,提高瓦斯抽采效率。 根据现场突出危险性鉴定,寺家庄煤矿 号煤 层为突出煤层。 受强构造应力影响,该煤层渗透率 低,煤质极松软。 为了防止北翼辅助运输大巷掘进 过程发生煤与瓦斯突出,必须预先抽采瓦斯,再掘进 施工。 由于普通钻孔瓦斯抽采应力扰动范围小,抽 采效率低,抽采时间长,因此采用水力钻扩造穴技术 提高煤层的瓦斯抽采效率。 作者以弹性力学软化理 论和双孔介质瓦斯流动理论为基础建立了造穴钻孔 周围煤体渗透率演化模型和瓦斯运移模型,并利用 多物理场数值模拟软件对建 立的模型进行解算,分析了造穴钻孔周围煤体的渗 透率和应力分布规律,对比了不同钻孔半径和钻孔 间距下的煤层瓦斯压力、瓦斯抽采量,并在寺家庄煤 矿 号煤层对模拟结果进行了现场验证。 钻扩造穴卸压增透机理及瓦斯流动模型 钻扩造穴卸压增透机理 钻孔周围煤体应力分布 煤体是一种弹塑性材料,当所受应力超过峰值 应力后会表现出应变软化特性,采取水力造穴措施 可使钻孔周围煤体应力重新分布。 根据煤岩体力学 特性应变软化模型,煤体受力破坏的应力变化过程 可分为 个阶段(图 )弹性阶段、软化阶段和残余 阶段。 钻扩造穴施工后,在钻孔周围的煤体相应形 成弹性区、塑性区和破碎区。 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 王 亮等瓦斯抽采穿层钻孔钻扩造穴卸压增透机理研究 年第 期 图 煤岩体应变软化模型及钻孔周围煤体应力分布特征 等效塑性剪切应变可以作为描述应变软化过程 的软化参数为 ( , , ,) () 式中 为等效塑性剪切应变; ,、 ,、 ,为沿 个主应变方向的塑性主应变。 应变软化过程是在内摩擦角不变的情况下失去 黏聚力的过程,黏聚力随着等效塑性剪切应变的增 加而呈线性减小,在塑性条件下,全应力应变曲线上 的黏聚力可以表示的等效塑性剪切应变分段线性函 数为 ( ) ( ) ( ) () 式中 为黏聚力;为初始黏聚力;为残余黏聚 力; 为临界塑性剪切应变。 数值模拟中,通常使用 ()准 则的六边形与 ()失稳准则外接 圆进行匹配作为岩石的破坏准则。 准则考 虑了静水压对岩石破坏的影响,认为材料的破坏由 偏应力第一不变量和第二不变量共同决定,其表达 式为 () 式中、为偏应力第一不变量、第二不变量, ; 和 为材料常数; 为煤的内摩擦角,()。 渗透率模型 煤层渗透率是控制煤层瓦斯流动的关键参数, 主要受地应力改变引起的煤体裂隙变化影响。 水力 冲孔使煤层应力重新分布,煤体发生变形和破坏,裂 隙宽度增大,数量增多,进而煤层渗透率增加。 煤体 所处的应力状态不同,发生的变形不同,因此渗透率 随地应力状态的变化遵循不同的函数关系。 根据前 人研究, 在弹性阶段 ( )、 塑性软化阶段 ( )、残余阶段( ),水力冲孔钻孔 周围煤体渗透率与体积应力变化量的关系满足下 式,即 () ( ) ( ) () ( ) ( ) () ( ) () 式中 为渗透率,;为初始渗透率,;为裂隙 压缩因子, ;为体积应力,; 为渗透率 跃变系数。 瓦斯流动模型 基质瓦斯扩散 煤基质的瓦斯流动遵循质量守恒方程为 () ( ) () 式中为单位体积煤基质中的瓦斯质量, ; 为单位体积煤基质同裂隙系统的质量交换率, ();为朗格缪尔体积, ;为孔隙 瓦斯压力,;为朗格缪尔压力,;为甲烷 的摩尔质量, ; 为煤体视密度, ;为 气体摩尔体积, ; 为煤基质孔隙 率,; 为理想气体常数, (); 为煤层温 度,; 为吸附时间,等于煤体中 的瓦斯解吸 出来的时间,;为裂隙瓦斯压力,。 将式()代入质量守恒式(),得基质瓦斯流动 方程,即 ( ) ( ) ( ) () 裂隙瓦斯渗流 裂隙中瓦斯流动满足质量守恒方程为 () () 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 年第 期煤 炭 科 学 技 术第 卷 式中为单位体积煤体裂隙瓦斯质量,; 为煤 体裂隙率,; 为裂隙瓦斯密度, ;为煤体 裂隙瓦斯流动速度, ; 为渗透率,; 为甲烷 动力黏度,。 将方程()代入质量守恒式(),得裂隙瓦斯流 动的方程为 ( ) ( ) () 模拟与结果 以上建立了可用于数值模拟的钻扩造穴后煤层 渗透率演化和瓦斯流动方程。 笔者采用 数值模拟软件对以上方程进行了计算 求解,数值模拟中选用了 软件中的 个模 块钻孔开挖选择固体力学模块,基质瓦斯扩散选择 模块,裂隙瓦斯渗流选择达西定律模块。 模型描述和输入参数 根据山西阳煤寺家庄煤矿北翼辅助运输大巷底 板预抽巷道现场穿层水力冲孔情况,本次模拟的几 何模型及边界条件如图 所示。 为消除边界的影 响,建立了 的矩形区域,远大于钻孔的 增透区。 二维模型是真实三维煤层的简化, 模型 不会影响计算的精度,并且利于显示应力和瓦斯压 力分布。 对于固体变形模型,模型的左边和上 边为恒定应力条件,代表煤层的水平地应力,右边和 下边为滚轴边界。 对于瓦斯抽采模型,钻孔内为恒 定压力边界,四周为无流动边界条件。 图 几何模型和边界条件 所有参数均通过实验室试验获得或现场测量得 到。 模拟中所用参数如下 煤层泊松比 煤层弹性模量 初始水平应力 原始瓦斯压力 煤层初始黏聚力 煤层残余黏聚力 煤体视密度 ( ) 煤层内摩擦角 () 煤层软化参数的临界值 原始渗透率 渗透率突变系数 裂隙压缩因子 煤体裂隙率 基质孔隙率 朗格缪尔体积 ( ) 朗格缪尔压力 摩尔质量 ( ) 气体摩尔体积 ( ) 动力黏度 () 气体常数 ( ) 吸附时间 煤层温度 造穴半径对增透效果的影响 以造穴半径 为例,分析钻孔周围煤体的 应力、等效塑性应变分布。 钻孔施工后,钻孔周围煤 的径向应力、切向应力、体积应力和等效塑性剪切应 变 如图 所示。 根据平面应变假设,为常数, 体积应力为 和 之和。 由图 可知,钻孔周围的 ,体积应力很小。 随着距钻孔的距离增加,等 效塑性应变急剧下降,体积应力增加。 时,煤 体处于应变软化阶段; 时,煤体处于残余 阶段,发生塑性破坏。 因此,钻孔周围塑性破坏区为 ,软化区的范围为 ,增透区的厚度为 ,即等效增透半径为 。 当距钻孔中心 图 造穴钻孔周围应力和等效塑性应变分布 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 王 亮等瓦斯抽采穿层钻孔钻扩造穴卸压增透机理研究 年第 期 距离超过 时, ,体积应力为原始值,煤 体处于弹性阶段。 为研究造穴半径对煤体渗透率分布的影响,确 定合理的造穴半径,分别对造穴半径为 、、 、 和 情况下钻孔周围煤体的渗透率演 化特征和瓦斯压力进行了对比分析。 如图 所示, 实施水力造穴后,钻孔周围煤体渗透率得到提高,且 距钻孔越近,煤层渗透率越大,不同造穴半径下煤层 渗透率的最大值均为原始渗透率的 倍。 随着 距离远离钻孔,煤体渗透率与原始煤层渗透率比 ( )逐渐减小,最终趋近于 ,即煤体渗透率等于 原始煤层渗透率。 随着造穴半径的增大,煤层增透 范围增大,表明造穴半径越大,钻孔增透半径越大。 如图 所示,钻孔增透半径随造穴孔半径呈线性增 加,增透半径约为造穴半径的 倍,说明瓦斯抽 采影响半径随着造穴半径的增加而增大。 同时,图 显示了不同造穴半径抽采 后的煤层瓦斯压 力,随着造穴半径增加,煤体瓦斯压力减小,且随着 距离钻孔中心距离的增加瓦斯压力逐渐增大,最终 趋近于恒定值(煤层原始瓦斯压力 ),表明通 过水力钻扩造穴技术增加钻孔的半径可以对钻孔周 围较大范围内的煤体进行卸压增透。 图 不同造穴半径钻孔周围渗透率、增透半径及瓦斯压力分布 , 由图 可知,造穴半径越大,增透范围越大,越 利于瓦斯抽采,但这并不表明现场工程实施时造穴 半径越大越好,一方面造穴半径过大会造成单个钻 孔的成本增加,另一方面造穴半径受技术和钻机的 条件限制。 由于现场施工过程中,机械钻杆刀臂的 最大半径为 ,因此,结合现场条件,施工造穴半 径采用 。 钻孔间距对瓦斯抽采的影响 井下煤层瓦斯抽采极其复杂, 个区域内有大 量的抽采钻孔,因此钻孔间距是影响瓦斯抽采的 个重要因素。 对于同 个抽采区域,钻孔间距越小, 钻孔数量越多,同时抽采效率越高。 但是钻孔数量 的增加意味着成本的增加,因此,确定合理的钻孔间 距是非常重要的。 当钻孔半径为 时,单个钻 孔的增透半径为 。 为了使煤层的增透效果最 优化,钻孔间距应为 。 然而,由于多孔抽采孔 间互扰的影响,钻孔间距过小会降低每个钻孔的抽 采效率。 因此,为确定合理的钻孔间距,模拟了钻孔 间距为 共 种工况下的瓦斯抽采情况。 不同钻孔间距下抽采 后煤层瓦斯压力云 图如图 所示,煤层中间监测线瓦斯压力分布如图 所示。 根据图 可以看出,钻孔间距为 时, 个 钻孔间煤层瓦斯压力最低,随着钻孔间距的增大,瓦斯 图 不同钻孔间距的煤层瓦斯压力分布 压力逐渐增加。 以消除煤层突出危险性的瓦斯压力界 限 为界,钻孔间距小于 均符合要求。 不同钻孔间距下 内的抽采瓦斯累积量(图 )。 相同时间内,瓦斯抽采总量随着钻孔间距的增 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 年第 期煤 炭 科 学 技 术第 卷 加而增加,但增加量逐渐减小。 钻孔间距为 、、 和 时, 内瓦斯抽采总量分别为 、 、 和 ,相比于前一种工 况分别增加 、和 。 钻孔间距 为 、、 和 时, 内瓦斯抽采总量分别为 、 、 和 ,相比于前 一种工况分别增加 、、和 。 因此,为减少钻孔之间互相干扰的影响,若以增长率 为标准,合理的钻孔间距为 ;若以增长率 为标准,合理的钻孔间距为 。 图 不同钻孔间距的瓦斯抽采量 现场应用试验 考察地点位于阳泉矿区寺家庄煤矿 号煤层 北翼辅助运输大巷。 在巷道的掘进过程中,由于煤层 的起伏,需要进行石门揭煤。 揭煤区域全长 , 平均煤厚 ,煤层原始瓦斯含量为 ,煤层具有突出危险性。 在大巷掘进前, 沿着巷道方向分别布置普钻区和造穴区,以穿层钻 孔的方式进行瓦斯抽采。 普通钻孔区域长度为 ,钻孔间排距为 ,布置 共 排瓦斯抽采 钻孔,每排含 个瓦斯抽采钻孔,巷道两侧钻孔布置 范围分别为 和 。 在普钻区之后,布置造穴 钻孔区,钻孔区域长度为 ,巷道两侧钻孔布置 范围不变。 造穴区布置 共 排钻孔,钻孔 排间距为 。 在奇数排(、、、、、 排) 施工 、、、、 号钻孔并进行水力钻扩造穴,在偶 数排(、、、、 排)施工 、、、 号钻孔并 进行水力钻扩造穴。 图 为普钻钻孔与造穴钻孔布 置平面图和剖面图。 造穴钻孔施工过程中,钻杆通过穿层钻孔的方 式深入到巷道周围煤体内,同时打开高压水射流泵, 高压水射流从钻冲两用钻头上的造穴喷嘴喷出,形 成造穴孔洞,并在钻穴孔段中往复钻进回退钻杆, 通过水流的作用将破碎煤冲出孔洞。 当孔洞出煤停 止,冲出的水流变清时,停止造穴,从而实现煤体的 高效卸压增透。 巷道造穴区施工过程中造穴水压为 ,平均水压约 ;单穴造穴时间和 出煤量分别为 、 ;根据出煤量、煤 层厚度和煤的密度计算出实际单穴造穴半径平均值 为 。 然而,第 节确定的最优单穴造穴半 径为 ,这是由于现场施工过程中采用机械和 水力联合造穴,机械钻杆刀臂的最大半径为 , 而高压水射流增大了造穴半径,因此实际造穴半径 均大于设计的 。 图 普钻钻孔与造穴钻孔布置平剖面示意 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 王 亮等瓦斯抽采穿层钻孔钻扩造穴卸压增透机理研究 年第 期 在寺家庄矿 号煤层不同区域的瓦斯抽采过程中, 对普钻区和造穴区的瓦斯抽采数据进行了计量,普 钻区瓦斯抽采时间约 个月,造穴区瓦斯抽采时间 约 个月。 图 记录了瓦斯抽采措施实施 以 来,普钻区和造穴区的瓦斯抽采纯量和抽采浓度的 情况。 图 表明,抽采过程中,普钻区瓦斯抽采纯 量为 ,而造穴区瓦斯抽采纯量为 ,明显高于普钻区,且修正后造穴 区的平均瓦斯抽采纯量约为普钻区的 倍。 图 表明,普钻区瓦斯抽采体积分数为 ,平均 瓦斯抽采体积分数约 ;而造穴区瓦斯抽采体积 分数为 ,平均瓦斯抽采体积分数可达 ,相对普钻区提高了 倍左右。 图 普钻区和造穴区瓦斯抽采纯量和瓦斯抽采浓度对比 北翼辅运底抽巷普通穿层钻孔和水力钻扩造穴 钻孔瓦斯抽采技术的经济效益情况如图 所示。 对 普钻区进行 个月的瓦斯抽采后,共抽采瓦斯 ;而对造穴区进行 个月的瓦斯抽采后, 共抽采瓦斯 ,瓦斯抽采总量约提高了 倍;瓦斯抽采率由原来的 提高到目前的 ,提高了约 倍;瓦斯抽采纯量由原来的 提高到目前的 ,提高了约 倍。 采用新技术之后,瓦斯抽采周期由原来的 降低为目前的 (含钻孔施工过程中的瓦 斯抽采);钻孔施工工程量由原来的 (煤中进 尺 )降低为目前的约 (煤层中进尺 ); 长度区域的钻孔数量由原来的 个 减少为目前的 个。 以上结果均表明,实施水力钻 扩造穴技术后,瓦斯抽采经济效益显著提高。 图 瓦斯抽采技术经济效益对比分析 结 论 )以弹性力学应变软化模型为基础建立了钻孔钻 扩造穴后煤层的渗透率演化方程,分析了钻孔造穴的 增透机理,并以扩散渗流煤层瓦斯流动理论模型为基 础建立了穿层钻孔造穴煤层瓦斯流动方程。 )利用 多物理场数值模 拟软件求解了寺家庄煤矿北翼辅助运输大巷底板预 抽巷道对 号煤层进行钻扩造穴情况,分析了造穴 钻孔周围煤层渗透率分布和瓦斯抽采情况,并得出 了合理的造穴半径 ,钻孔间距 。 )经现场工程试验 号煤层的瓦斯抽采过 程,实施水力钻扩造穴技术后,瓦斯抽采率、抽采浓 度和抽采纯量分别提高了约 、、 倍;普钻区 经 个月的瓦斯抽采后,共抽采瓦斯 ;造 穴区经 个月的瓦斯抽采后,共抽采瓦斯 ,瓦斯抽采总量约提高了 倍,施工工程量和瓦 斯抽采周期降低了约 倍。 参考文献() , ,, ,,() , ,, , , , ,,() , , , 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et 年第 期煤 炭 科 学 技 术第 卷 ,,() 王永革,令狐建设 地面井水力压裂增透技术在立井快速揭煤 中的应用 矿业安全与环保,,() , ,,() 曹建军 超高压水力割缝卸压抽采区域防突技术应用研究 煤炭科学技术,,() ,,() 邓 强 水力压裂消突技术在低透气性煤层瓦斯治理的应用 煤矿安全,,() ,,() 刘明举,崔 凯,刘彦伟,等 深部低透气性煤层水力冲孔措施 防突机理分析 煤炭科学技术,,() , , , ,,() 王峰,陶云奇,刘东 水力冲孔卸压范围及瓦斯抽采规律 研究 煤炭科学技术,,() , , ,,() , , , , ,,() , , ,, 刘东,刘文 水力冲孔压裂卸压增透抽采瓦斯技术研究 煤炭科学技术,,() , ,, () 孙四清,张 俭,安鸿涛 松软突出煤层穿层洞穴完井钻孔瓦 斯抽采实践 煤炭科学技术,,(), , , ,,(), 杜昌华,冯仁俊 低透松软破碎厚煤层水力扩孔增透技术研 究 煤炭科学技术,,() , ,,() 牟全斌,赵继展 基于机械造穴的钻孔瓦斯强化抽采技术研 究 煤炭科学技术,,(), , ,,(), 郝从猛,刘洪永,程远平 穿层水力造穴钻孔瓦斯抽采效果数 值模拟研究 煤矿安全,,() , , ,,() 赵继展 井下瓦斯抽采钻孔机械造穴技术研究 矿业安全 与环保,,(), ,,(), 于宝种 松软低透煤层高压射流造穴强化抽采技术研究 矿业安全与环保,,(), ,,(), 张 浩构造煤层掘进工作面区域性顺层水力造穴强化瓦斯 抽采机理与工程应用徐州中国矿业大学, , , , ,,() , , ,() , ,,() , ,,() , , ,, , , ,() , , , ,,() , , , ,, () , , ,,() 中国煤炭行业知识服务平台 w w w . c h in a c a j . n et
展开阅读全文