废弃采空区地基稳定性极简评价指标体系构建-矿业114网 
首页 >> 文献频道 >> 矿业论文 >> 正文
废弃采空区地基稳定性极简评价指标体系构建
2019-09-29
废弃采空区不良地基残余变形威胁其上方拟建工程的运营安全,对其开发利用前需进行地基稳定性 评价,而当前地基稳定性评价指标体系中因素纷繁复杂且数量庞大,极易导致评价结果失去科学性和合理性。以 河南古汉山矿和吴村矿废弃采空区上方的武云高速建设为例,基于研究区域的实际地质采矿环境,首先遵从唯一 性、目的性、可行性和可观测性四大原则对初选指标进行定性筛选;其次,基于改进灰色关联度算法对定性筛选出 的结果再次进行定量筛选,确定了煤层埋深、建(构)筑物荷载扰动深度与导水裂隙带垂向间距以及停采时间作为评 价指标,进而构建了适用于武云高速建设地基稳定性评判的极简评价指标体系,通过引入效度系数和相关性系数 分析了...
Series No. 519 金 属 矿 山 总第 519 期 September 2019 METAL MINE 2019 年第 9 期 废弃采空区地基稳定性极简评价指标体系构建 1 ,2 3 1 11 鑫 郭庆彪 李英明 王 亮 吕 ( 1. 安徽理工大学测绘学院,安徽 淮南 232001;2. 煤炭资源与安全开采国家重点实验室,江苏 徐州 221116; 3 . 安徽理工大学能源与安全学院,安徽 淮南 232001) 摘 要 废弃采空区不良地基残余变形威胁其上方拟建工程的运营安全,对其开发利用前需进行地基稳定性 评价,而当前地基稳定性评价指标体系中因素纷繁复杂且数量庞大,极易导致评价结果失去科学性和合理性。以 河南古汉山矿和吴村矿废弃采空区上方的武云高速建设为例,基于研究区域的实际地质采矿环境,首先遵从唯一 性、目的性、可行性和可观测性四大原则对初选指标进行定性筛选;其次,基于改进灰色关联度算法对定性筛选出 的结果再次进行定量筛选,确定了煤层埋深、建(构)筑物荷载扰动深度与导水裂隙带垂向间距以及停采时间作为评 价指标,进而构建了适用于武云高速建设地基稳定性评判的极简评价指标体系,通过引入效度系数和相关性系数 分析了该指标体系的有效性和可靠性。研究表明:运用所构建的极简指标体系评价武云高速建设地基处于基本稳 定状态,评价结果与其他指标体系的评价结果较一致,且由于所构建的指标体系仅含有 3 个主要指标,可极大减轻 系统的计算负担,简化评估过程。 关键词 废弃采空区 地基稳定性 指标筛选 改进灰联度算法 中图分类号 TD353 文献标志码 A 文章编号 1001-1250(2019)-09-179-06 DOI 10.19614/j.cnki.jsks.201909029 Construction of Minimalist Foundation Stability Evaluation Index System for Abandoned Goaf 1,2 3 1 12 Guo Qingbiao Li Yingming Wang Liang Lu Xin (1. Anhui University of Science and Technology,Huainan 232001,China;2. State Key Laboratory of Coal Resources and Safe Mining,Xuzhou 221116,China;3. School of Energy and Safety,Anhui University of Science and Technology, Huainan 232001,China) Abstract The residual deformation of foundation in an abandoned goaf threatens the safety operation of the proposed project,and the foundation stability evaluation is required before its development and utilization. However,the current foun⁃ dation stability evaluation index system is complicated,which can easily lead to unscientific and unreasonable evaluation re⁃ sults. In view of this,this paper takes the Wuyun expressway construction as an example,which is located above the aban⁃ doned goaf of Guhanshan Mine and Wucun Mine in Henan Province.First of all,according to the four principles of uniqueness, purpose,feasibility and observability,qualitative screening is conducted.Secondly,screened quantitatively is conducted based on the theory of improved grey relational degree,and then a minimalist stability evaluation system suitable for Wuyun expressway construction foundation is established.The effectiveness and reliability of established minimalist stability evalua⁃ tion system is analyzed by using the validity coefficient and correlation coefficient.Wuyun expressway construction foundation is convicted in a basically stable state by using the proposed minimalist stability evaluation system,and the evaluation result is consistent with those of other evaluation index systems,confirming its accuracy. Because the minimalist stability evaluation system only contains three main indicators,it can greatly reduce the calculation burden of the system and simplify the evalua⁃ tion process. Keywords Abandoned goaf,Foundation stability,Indexes screening,Improved grey relational degree algorithm 收稿日期 2019-07-24 基金项目 国家自然科学基金项目(编号:51574006),煤炭资源与安全开采国家重点实验室开放研究基金项目(编号:SKLCRSM19KF005),安徽 理工大学国家自然科学基金预研项目(编号:YY201802),安徽理工大学博士科研启动基金项目(编号:11953)。 作者简介 郭庆彪(1990—),男,讲师,博士。通讯作者 李英明(1975—),男,教授,博士,博士研究生导师。 · 179 · 总第519期 金 属 矿 山 2019年第9期 废弃采空区内赋存大量空洞、空隙,在内外荷载 耦合扰动下,其内部空洞、空隙被压密,由此释放出 移动变形并向上逐层传递,直至地表,引起地表(建 筑地基)产生残余变形,威胁其上方拟建工程安全运 营。为此,在进行废弃采空区土地资源开发利用前, 有必要进行地基稳定性评价。废弃采空区地基稳定 性与煤层埋深、采厚、地质构造、覆岩岩性、开采方法 等多因素密切相关,使得地基稳定性评价是一项极 的实际地质采矿环境,开展研究区域地基稳定性评 价主要指标的筛选研究,进而构建适用于武云高速 建设地基稳定性评价的极简评价体系。 1 工程概况 武云高速位于河南省焦作市修武县、武陟县境 内,是河南省高速公路“十一五”规划的重点工程。 由现场调研资料可知,武云高速 K29+335~K33+000 路 段 穿 越 古 汉 山 煤 矿 和 吴 村 煤 矿 ,其 中 K30 + 800~K32+900路段位于老采空区上方,二者之间平面 相对位置关系如图1所示。 [1-2] 其复杂的系统工程,其稳定状态难以精确评定 [3] 。 吴启红等 运用多级模糊评判理论,选取 15 个评判 因子,建立了矿山采空区稳定性的II级模糊综合评价 模型,并将该模型应用于具体工程,其评价结果与实 研究区地表海拔 90~115 m,平均坡度 0.542 1%。 下伏采空区均为全部垮落法管理顶板。浅部地层主 要为第四系地层,深部地层主要为二叠系地层,地层 平均倾角为15°,区内基本构造轮廓为向南东倾斜的 单 斜 构 造 ,仅 在 K30 + 520~K29 + 650 处 和 K29 + 520~K29+650段下方赋存团相断层和油坊蒋村断层, 断层数量少。 [ 4] 际情况相吻合;高峰等 根据采空区稳定性评价中评 价指标间交互作用机制,提出了基于 RES-云模型的 采空区稳定性二维评价方法,并将其应用于云南某 [4] [5] 锡矿,验证了该模型的可行性 ;邓红卫等 建立了 基于不确定型层次分析法(UAHP)的采空区稳定性模 糊评价方法,有效解决了评价指标权重确定难题;董 2 武云高速建设地基稳定性评价指标筛选 2. 1 定性筛选 [ 6] 宪久 选取了工程地质环境、岩体质量、采空区几何参 数和围岩应力等因素构建了采空区稳定性评价指标 体系,基于功效系数法原理,采用三标度改进层次分 析法确定了各评价指标的权重系数,建立了采空区稳 定性评价分级新模型,其适用性和可靠性得到了工程 结合废弃采空区地基稳定性评价的研究成 15-16] ,梳理并归纳了与武云高速建设地基稳定性 [ 果 有关的评价指标,如表1所示。 由表1可知:与武云高速建设地基稳定性有关的 评价指标众多,部分指标之间具有一定的相关性,表 达信息冗余,如此庞大的指标集不仅大幅增加计算工 作量,还会淡化主要评价指标的重要程度。为此,本 研究遵循唯一性、目的性、可行性和可观测性四大原 则,首先进行定性筛选,以保证定性筛选后的指标可 进一步进行量化分析。其中唯一性、目的性、可行性 和可观测性原则可表述为:①唯一性,删除信息重叠、 冗余和无用指标 ②目的性,删除无法切实反应废弃 采空区地基稳定性特征的指标;③可行性,筛选出的 指标数据能够保质保量获取;④可观测性,无论是直 接还是间接观测,筛选出的指标必须含义清晰、明确。 基于上述原则,对表 1 中的指标进行了定性筛 选,筛选后的结果如表2所示。 [ 7] 实际的有效验证;康钦容等 采用综合指数方法对大 面积石膏矿采空区失稳的危险性进行了综合评价。 上述研究表明:将非线性科学理论应用于评价采空区 [ 8-11] 稳定性近年来得到了广大学者们的深入研究 。 反观上述非线性评价模型,学者们为全面反映 废弃采空区地基稳定性的真实状态,从水文地质、采 空区自身及外荷载扰动等方面选取较多的评价指 标,构建多维量化评价指标体系。然而,在实际应用 中,评价指标并非越多越好,庞大的指标体系极易导 [ 12-14] 致评价结果失去科学性和合理性 ,存在的不足主 要有:①庞大的指标体系中不可避免地存在着相关 性指标,表达信息冗余,增加计算成本;②过多的评 价指标容易引起评价者判断上的错觉和混乱,导致 关键指标权重减小,造成评价结果失真;③由于评价 者对于定性指标的理解各不相同,在采用同一指标 体系评价相同目标时,会产生相差较大的评价结果, 即指标体系的稳定性和可靠性较差。 2. 2 定量筛选 本研究采用改进的灰色关联度算法对上述定性 筛选后的指标进行二次筛选。 2. 2. 1 改进的灰色关联度模型 为避免上述问题,有必要对稳定性评价指标进 行有效筛选,在庞大的指标体系中筛选出少量具有 重要参考价值的指标,进而构建废弃采空区地基稳 定性评价的极简化指标体系。本研究以河南武云高 速穿越废弃采空区的具体工程为例,根据研究区域 灰色关联是将样本数据分为参考数据和评价数 据,分别计算各评价数据与参考数据间的关联度,关 [ 17-18] 联度越大,则该评价数据越重要 。传统灰色关联 度计算易受分辨系数取值的影响,计算结果具有主 观不确定性,影响最终决策。本研究采用改进的灰 · 180 · 郭庆彪等:废弃采空区地基稳定性极简评价指标体系构建 2019年第9期 色关联度对指标权重进行求解,具体过程如下。 1)邀请m位行业内专家同时对n个指标进行权 重赋值,得到表征指标权重的矩阵 D ( éd11 d d21 d22 ... d2m 12 ... d ù 1m ú ê D = , (1) ê ê ú ú . .. ... ... ... dn1 dn2 ... dnm ë û 式中,dnm 为第 m 个专家对第 n 个待筛选指标权重的 评定值。 (2)从矩阵 D 中每列挑选出最大值作为参考值, 组成参考向量 D0 · 181 · 总第519期 金 属 矿 山 2019年第9期 D0 = { d01,d02,d03,⋯,d0m} , (2) 式,d0m 为第m位专家确定的参考权重。 ( 3)求取各个指标与参考指标数据之间的距离, m 2 D0i = (d0k - dik) , (3) ∑ k = 1 式中,D0i 为第 i 个指标与参考指标之间的距离;dik 为第k位专家对第i个指标赋予的权值。 (4)计算各指标的权重。公式为 wi = 1/(1 + D0i) , (4) n - w = wi / wi , ∑ i = 1 (5) i - 系进行评价时,对问题的认识是否趋向一致,越趋向 一致,则指标体系的有效性越高,也就越能反应实际 情况。本研究采用效度系数来判定筛选后指标的有 式中,w 为第i个指标的权重。 为第 i 个指标与参考指标数据的关联度;wi i 2 . 2. 2 武云高速建设地基稳定性评价指标定量筛选 邀请业内 5 名专家结合武云高速建设地基的实 [19-20] 效性 。 设评价体系 F ={ f1,f2,⋯,fi,⋯fn} ,参加评价的专 家 人 数 为 S,专 家 j 对 评 价 目 标 的 评 分 集 为 际地质采矿环境,对表 2 中的评价指标进行权重赋 值,形成了指标权重矩阵D é0.25 0.18 0.20 0.20 0.24ù X = { x 1j,x2j,⋯,xij,⋯,xnj } S ,定义指标 f 的效度系数为 i ê0.10 0.12 0.18 0.15 0.13ú ê ú ú ú . — ê0.06 0.10 0.08 0.07 0.08ú ê βi = |xi - xij| SM , (8) ∑ ê0.10 0.13 0.09 0.10 0.10ú ê j = 1 D = (6) ê0.05 0.08 0.05 0.04 0.05ú ê ú ú — ê0.04 0.04 0.05 0.04 0.05ú ê 式中,M为指标 fi 的评分集中评分最优值;xi 为评价 ê0.18 0.15 0.15 0.18 0.15ú ê ú ë0.22 0.20 0.20 0.22 0.20û 由此组成参考向量 D0 D0 = 0.25,0.20,0.20,0.22,0.24 指标 fi 得分平均值,可由下式计算出 S — xi = xij S . (9) ∑ { } . (7) j = 1 则由筛选后指标体系F的效度系数为 联合式(3)~式(5),计算的定性筛选后的指标权 重如表3所示。 n β = βi n . (10) ∑ i = 1 效度系数越小,表明指标体系的有效性越高。 一般情况下,当指标体系效度小数小于 0.1 时,则该 指标体系的有效性较高。结合式(6)与式(8)~(10), 计算的图 2 指标评价体系的效度系数为 0.076,即本 研究建立的武云高速建设地基稳定性评价指标体系 有效性较高。 3 . 2 可靠性分析 若多个专家采用相同的评价指标体系进行稳定 性评价时,评价结果差异较小,则认为该指标体系可 靠性强。本研究采用相关系数分析武云高速建设地 [ 19-20] 分析表3可知:与武云高速建设地基稳定性关联 度最高的前 3 个指标分别是指标 1、指标 8 和指标 7, 即煤层埋深、采空区停采时间和建(构)筑物荷载扰 动深度与导水裂隙带垂向间距这 3 个指标最能代表 原始的评价指标体系,如图2所示。 基稳定性评价指标体系的可靠性 。设专家组评 分的平均数据组 Y ={y1,y2,⋯,yi,⋯,yn} ,其中 S yi = xij S . (11) (12) ∑ j = 1 则指标体系可靠性系数为 S 3 筛选指标的有效性和可靠性 . 1 有效性分析 ρ = ρj S , ∑ 3 j = 1 筛选指标的有效性是指不同专家采用该指标体 式中,ρj 可由下式计算出 · 182 · 郭庆彪等:废弃采空区地基稳定性极简评价指标体系构建 2019年第9期 n — xij - x ÷çyi - y j — ì 废弃采空区地基稳定性评价指标的筛选流程,包括 指标初选、定性筛选和定量筛选。在对初选指标定 性筛选后,采用改进灰色关联度算法进一步进行定 量筛选,最终确定煤层埋深、采空区停采时间和建 (构)筑物荷载扰动深度与导水裂隙带垂向间距为主 要指标,来评价武云高速建设地基的稳定性。 æ öæ ö ÷ ç è ï ï ∑ øè ø i = 1 ïρj = ï n 2 n 2 ï ï — — æ ö æ ç yi - y ∑ è i = 1 ö ç xij - x ÷ ÷ ï ï ∑ j è ø ø i = 1 . (13) í n — ï ï xj = xij n ï ï ∑ i = 1 n ï ï — ï y = yi n ï ∑ (2)利用筛选出的3个指标建立了适用于武云高 î i = 1 速建设地基稳定性评价的极简指标体系,引入效度 系数和相关性系数,分析得到该评价指标体系具有 较高的有效性和可靠性。 一般情况下,当 ρ ∈ 价体系可靠性较高;当 ρ ∈ 评价体系可靠性一般;当 ρ ∈ ( 0.90,0.95 ) 时,可认为该评 ( 0.80,0.90 ) 时,可认为该 ( 0,0.80 时,可认为该 ) (3)运用本研究构建的极简评价指标体系判定 评价体系可靠性较差。结合上式计算的可靠性系数 为 0.906,即本研究构建的武云高速建设地基稳定性 评价指标体系可靠性较高。 武云高速建设地基处于基本稳定状态,通过在路基 中加铺土工格栅的方式即可保证武云高速公路运营 安全。该评价结果与其他评价指标体系的评价结果 一致,表明该评价指标体系能够有效评估武云高速 建设地基的稳定状态,且由于该指标体系仅含有3个 指标,有助于大幅减轻系统的计算负担,简化评价过 程。 4 武云高速建设地基稳定状态评价 以模糊综合评判理论为基础,采用本研究构建 的极简评价指标体系对武云高速建设地基的稳定性 进行了评价。对表3数据进行了标准化处理,得到评 价体系中的指标权重,如表4所示。 参 考 文 献 [1] 郭庆彪.煤矿老采空区上方高速公路建设安全性评价及其关键 技术研究[D].徐州:中国矿业大学,2017. Guo Qingbiao.Research on the Safety Evaluation and Key Technolo⁃ gies for the Expressway Construction on Old Goaf of Coal Mine[D]. Xuzhou:China University of Mining and Technology,2017. 参考文献[20]中上述3个指标的隶属度取值(式 [2] 杨 锋,郭广礼,万战胜,等.煤矿老采空区上方高速公路建设场 ( 14)),选取 M(⊕,·) 算子进行模糊合成运算,最终判 地稳定性评价[J].金属矿山,2017(11):137-144. Yang Feng,Guo Guangli,Wan Zhansheng,et al.Stability evalua⁃ tion of the expressway construction site on the old goaf of coal mine [J].Metal Mine,2017(11):137-144. 定武云高速建设地基处于基本稳定状态,通过采取 简易的抗变形设计方案即可保证武云高速公路运营 安全。 [ 3] 吴启红,彭振斌,陈科平,等.矿山采空区稳定性二级模糊综合评 é0.080 0.429 0.383 0.108ù ë0.100 0.700 0.200 0.000úû R = êê0.300 0.500 0.150 0.050ú . (14) 判[J].中南大学学报:自然科学版,2010,41(2):661-667. Wu Qihong,Peng Zhenbin,Chen Keping,et al.Synthetic judgment on two-stage fuzzy of stability of mine gob area[J].Journal of Cen⁃ tral South University:Science and Technology,2010,41(2):661- 根据评价结果,提出采用在路基中加铺土工格 栅的方式以提高武云高速的抗变形能力,规避潜在 威胁,如图3所示。 6 67. [4] 高 峰,高宇旭,周科平.基于RES-云模型的采空区稳定性评价 研究[J].灾害学,2019,34(1):17-26. Gao Feng,Gao Yuxu,Zhou Keping.Economic loss assessment of ty⁃ phoon based on BP neural network and VSM[J].Journal of Catastro⁃ phology,2019,34(1):17-26. [5] 邓红卫,贾 明.基于UAHP的采空区稳定性模糊综合评判[J]. 中国安全科学学报,2012,22(3):24-29. Deng Hongwei,Jia Ming.Fuzzy comprehensive evaluation of gob ar⁃ ea stability based on uncertain analytic hierarchy process[J]. Chi⁃ na Safety Science Journal,2012,22(3):24-29. 本研究评价结果与文献[2]的评价结果一致,证 实了所构建的极简评价体系的准确性,并反映出该 评价体系的优越性。 [ 6] 董宪久.基于功效系数与改进AHP法的地下矿山采空区稳定性 5 结 论 评判[J].有色金属:矿山部分,2018,70(4):24-29. (1)结合武云高速建设的具体工程实例,分析了 Dong Xianjiu.Stability evaluation study of underground mine goaf · 183 · 总第519期 金 属 矿 山 2019年第9期 based on efficacy coefficient method and improved AHP[J]. Non⁃ ferrous Metals:Mine Section,2018,70(4):24-29. Command Control,2018,43(2):37-42. [14] 李随成,陈敬东,赵海刚.定性决策指标体系评价研究[J].系统 工程理论与实践,2001(9):22-27. [ 7] 康钦容,张卫中,倪小山,等.采空区稳定性评价的综合指数法 [ J].金属矿山,2017(3):162-166. Li Suicheng,Chen Jingdong,Zhao Haigang.Studying on the method of appraising qualitative decision indication system[J].System En⁃ gineering Theory and Practice,2001(9):22-27. Kang Qinrong,Zhang Weizhong,Ni Xiaoshan,et al. Comprehen⁃ sive Index Method in Goaf Stability Evaluation[J]. Metal Mine, 2 017(3): 162-166. [15] 王正帅,刘冰晶,邓喀中.老采空区稳定性的模糊可拓评价模型 [J].地下空间与工程学报,2016,12(3):553-559. Wang Zhengshuai,Liu Bingjing,Deng Kazhong.Fuzzy extension as⁃ sessment model of old goaf stability[J].Chinese Journal of Under⁃ ground Space and Engineering,2016,12(3):553-559. [16] 王正帅.老采空区残余沉降非线性预测理论及应用研究[D].徐 州:中国矿业大学,2011. [ 8] 车安刚,郭庆彪,郭广礼,等.基于随机介质理论的武云高速老空 区段建设场地稳定性评价[J].公路,2016(12):53-57. Che Angang,Guo Qingbiao,Guo Guangli,et al.Stability evaluation of construction site for goaf section of Wuzhi- Yuntaishan express⁃ way based on stochastic medium theory[J].Highway,2016(12):53- 5 7. [ 9] 罗东海.模糊评判法在煤矿采空区稳定性评价中的应用[J].水 Wang Zhengshuai.Study on the Non- linear Prediction Theory of Old Goaf Residual Subsidence and Its Application[D].Xuzhou:Chi⁃ na University of Mining and Technology,2011. 土保持通报,2015,35(1):176-186. Luo Donghai.Application of Fuzzy evaluation method in stability evaluation of mined-out area[J].Bulletin of Soil and Water Conser⁃ vation,2015,35(1):176-186. [17] 张 磐,汪宇倩,姜 宁,等.基于改进灰色关联度的配电网调控 水平评价指标筛选方法[J].电力系统及其自动化学报,2018 (11):64-69. [10] 杨继华,路新景,闫长斌,等.基于模糊综合评判的护盾式 TBM 施工围岩稳定性分类[J].吉林大学学报:地球科学版,2018,48 Zhang Pan,Wang Yuqian,Jiang Ning,et al.Screening method for evaluation indicators of distribution network control level based on improved grey correlation degree[J].Proceedings of the CSU- EP⁃ SA,2018(11):64-69. (5):1596-1602. Yang Jihua,Lu Xinjing,Yan Changbin,et al.Stability classification of tunnel surrounding rock during shield TBM construction based on fuzzy comprehensive evaluation[J].Journal of Jilin University: Earth Science Edition,2018,48(5):1596-1602. [18] 方 齐.膏体带状充填开采复合支撑体稳定性模拟研究[D].徐 州:中国矿业大学,2016. [ 11] 葛鹏飞,郭广礼,郭庆彪,等.基于层次分析法的采空区稳定性评 Fang Qi.Simulation Study on Stability of Composite Support Pillar for Strip Mining with Paste Backfilling[D].Xuzhou:China Universi⁃ ty of Mining and Technology,2016. 价[J].煤矿安全,2017,48(7):228-231. Ge Pengfei,Guo Guangli,Guo Qingbiao,et al.Stability evaluation of goaf based on analytic hierarchy process[J].Safety in Coal Mines,2017,48(7):228-231. [19] 娄云永.电力通信网可靠性评价指标体系的研究[D].保定:华北 电力大学,2009. [ 12] 刘雨萌,李战江,尹 伟,等.非参数下贝叶斯判别与聚类分析的 Lou Yunyong.The Reliability Evaluation Index System Research of Eletric Power Communication Network[D].Baoding:North China Electric Power University,2009. 信用指标筛选模型[J].统计与决策,2018(22):5-10. Liu Yumeng,Li Zhanjiang,Yin Wei,et al.Credit index screening model of nonparametric Bayesian discriminant and cluster analysis [20] 郭广礼.高速公路下伏采空区沉陷机理及治理关键技术研究 [R].徐州:中国矿业大学,2017. [J].Statistics & Decision,2018(22):5-10. [ 13] 路云飞,李琳琳,张 壮.基于灰色粗糙集的指标筛选方法[J]. Guo Guangli.Study on subsidence mechanism and key technology of goaf under expressway[R].Xuzhou:China University of Mining and Technoloy,2017. 火力与指挥控制,2018,43(2):37-42. Lu Yunfei,Li Linlin,Zhang Zhuang.Research on approach of eval⁃ uation index screening based on grey rough set[J].Fire Control & (责任编辑 王小兵) · 184 ·
  • 中矿传媒与您共建矿业文档分享平台下载改文章所需积分:  5
  • 现在注册会员立即赠送 10 积分


皖公网安备 34050402000107号