基于VUMA-3D的深井矿山热环境模拟研究.pdf

返回 相似 举报
基于VUMA-3D的深井矿山热环境模拟研究.pdf_第1页
第1页 / 共5页
基于VUMA-3D的深井矿山热环境模拟研究.pdf_第2页
第2页 / 共5页
基于VUMA-3D的深井矿山热环境模拟研究.pdf_第3页
第3页 / 共5页
基于VUMA-3D的深井矿山热环境模拟研究.pdf_第4页
第4页 / 共5页
基于VUMA-3D的深井矿山热环境模拟研究.pdf_第5页
第5页 / 共5页
亲,该文档总共5页,全部预览完了,如果喜欢就下载吧!
资源描述:
第 27 卷增刊 2 2018 年 10 月 中 国 矿 业 CHINA MINING MAGAZINE Vol . 27 , Suppl 2 Oct . 2018 收稿日期 2018-08-12 责任编辑 赵奎涛 基金项 目 “ 十 三 五 ” 国 家 重 点 研 发 计 划 项 目 资 助 (编 号 2016YFC0600803) 第一作者简介 朱瑞军(1971 - ) , 男 , 山西阳泉人 , 高级工程师 , 从事 矿山工程咨询设计及科研工作 。 引用格式 朱瑞军 , 李少辉 , 杨卓明 . 基于 VUMA-3D 的深井矿山热 环境模拟研究 [J] . 中国矿业 , 2018 , 27 (S2) 234-237 , 262 . doi 10 . 12075/j . issn . 1004-4051 . 2018 . S2 . 024 矿业纵横 基于 VUMA-3D 的深井矿山热环境模拟研究 朱瑞军 1 , 2 ,李少辉 1 , 2 ,杨卓明 1 , 2 (1 .中国恩菲工程技术有限公司 ,北京 100038 ; 2 .中国矿业信息化协同创新北京市工程研究中心 ,北京 100038) 摘 要 深井高温热害已经成为制约矿山安全高效开采的重大问题之一 , 需要对深井矿山的热环境 进行三维模拟才能有针对性地采取降温措施 。 本文介绍了三维通风模拟软件 VUMA-3D 的功能特点 , 对 深井热源及其影响因素进行了分析 , 针对某深井矿山进行了热环境模拟 , 并对模拟结果进行了对比分析 , 提出了按需降温的理念 。 关键词 VUMA-3D ;深井矿山 ;热环境模拟 ;按需降温 中图分类号 TD727 +.2 文献标识码 A 文章编号 1004-4051(2018)S2-0234-04 Study on thermal environment simulation of deep mine based on VUMA-3D ZHU Ruijun1 , 2,LI Shaohui1 , 2,YANG Zhuoming1 , 2 (1 .China ENFI Engineering Corporation ,Beijing 100038 ,China ; 2 .China Mining Innovation Center ,Beijing 100038 ,China) Abstract Hot disaster of high temperature has become one of the biggest issues with safe and efficient mining .The 3D simulation of the thermal environment for deep mine can be targeted to take cooling measure . This paper introduces the functions and characteristics of the 3D ventilation simulation software VUMA-3D , and analyzes the heat source and its influencing factors .As a case ,simulates the thermal environment simulation of one deep mine , and analyzes the result of the simulation . At last , puts forward the concept of cooling on demand . Keywords VUMA-3D ;deep mine ;thermal environment simulation ;cooling on demand 随着地下矿山开采深度的不断增大 , 井下温度 也随之升高 , 深井高温热害已经成为制约矿山安全 高效开采的重大问题之一 。 美国 、 德国 、 南非等矿业 发达国家很早就开展了矿井热害治理研究工作 , 这 些国家的研究人员在矿井风温预测 、岩石热物性参 数测定 、 地温观测 、系数计算 、井下热源散热量以及 需冷量的计算方面做了大量的研究 , 为矿井热害治 理和井下空调应用奠定了基础 [1-4] 。 当前 , 我国对深 部矿井热环境的研究 , 多集中于理论研究 、数值计 算 、 室内试验等方面 , 还不能对井下热环境的温度场 及其规律进行准确预测 [5-7] 。 对矿井热害控制的前提是对井下热环境进行预 测 。 井下热环境是指由气温 、 周围物体表面温度 、 相 对湿度与气流速度等物理因素组成的作用于人并影 响人的冷热感和健康的环境 。 由于导致金属矿山矿 井热害的热源种类较多 , 不同热源对井下作业环境 的影响程度不同 。 解决金属矿山深井热害问题一般 要查找与分析矿井存在高温的主要原因 , 并利用计 算机软件进行热模拟 , 有针对性地提出解决问题的 技术方案 。 1 VUMA-3D 软件系统简介 VUMA-3D 是国际上比较先进和流行的三维 通风软件 , 能够对地下矿山通风系统 、 热环境 、污染 物环境等进行模拟 , 进行空气动力学 、热力学解算 。 VUMA-3D 软件系统是由专门从事矿井通风降温 技术研究和工程设计施工的南非布卢姆 伯顿工程 万方数据 增刊 2朱瑞军 ,等 基于 VUMA-3D 的深井矿山热环境模拟研究 有限公司(Bluhm Burton Engineering , BBE)基于最 新研究成果和多年的工程实践经验开发 。 在 VUMA-3D 软件系统中 , 通风网络是由多种 分支(元素)连接而成的 , 这些分支包括竖井 、 井底车 场 、 巷道 、 掘进工作面 、 采场 、 风机 、 通风构筑物等 , 每 一个分支都包含特定的空气动力学和热力学特性 。 通过空气动力学及热力学模拟 , 可用于确定风机机 站及构筑物布置方案 、井下制冷方案并进行能耗计 算等 。 VUMA-3D 还可基于在线监控和实时数据 分析技术 , 对生产矿山的通风系统进行优化 , 这对于 能耗控制以及降低运营成本有很大用处 。 VUMA-3D 软件系统的功能特点包括 ① 可进 行 dxf 文件的导入导出 , 方便快速建立通风网络 ; ② 友好易用的用户界面和智能提示功能 , 可进行通 风网络三维设计 , 方便地进行创建 、编辑 、 分类输入 参数等操作 ;③ 基于解算结果自动生成详细的分类 报告 ;④ 为地下采矿不同工作面的散热量计算提供 了完整算法 , 且这些算法都经过了实际工程验证 ; ⑤ 建立了岩石 、 设备和风机的数据库 , 方便快速调用 和计算 。 2 深井矿山热环境影响因素分析 2 . 1 围岩散热 围岩散热量的影响因素有原岩温度 、 岩石热导 率 、 暴露面积 、揭露时间 、表面湿润程度 、支护形 式等 。 岩石的热导率取决于岩石的成分 、 结构 、 形成条 件 、 含水状况 、 温度和压力等 [8] 。 一般情况下 , 岩石 的热导率随压力 、 密度和湿度的加大而增高 , 随温度 的增高而减小 。 岩石热导率一般应采用实测法确 定 , 不具备条件时 , 则参考类似岩石的热导率进行选 取 。 巷道的揭露时间会影响岩壁的散热量 , 揭露时 间短 , 则单位时间内的散热量较大 , 揭露一定时间 后 , 其单位时间内的散热量逐渐趋于稳定 。 岩壁的 湿润程度 、 巷道的支护形式和支护厚度也会影响围 岩的散热量 , 支护厚度大 , 则保温效果较好 , 围岩散 热量也会较小 。 2 . 2 设备散热 由于矿山开采强度的增大 , 井下的机械化程度 不断提升 , 从而使一些机械设备的装机容量及数量 相应增大 , 导致井下热量增多 , 矿井机械设备的散热 已成为深井矿山一个重要的热源 。 设备散热量的影 响因素包括设备类型 、设备数量 、单台功率 、负载情 况 、 设备利用率以及工作环境的潮湿程度等 。 柴油 设备与电动设备的散热量也有较大差别 , 同样功率 及工况条件下 , 柴油设备的散热量一般会大于电动 设备的散热量 。 2 . 3 空气压缩 风流由地表经井筒进入井下的过程中 , 位能转 换为焓 , 温度升高 , 且由于井筒内一般是潮湿的 , 井 巷围岩与风流会发生热湿交换的现象 , 从而使水分 蒸发带走热量 , 风流温度降低 , 导致自压缩对温升的 影响降低 , 但散热量会进一步加大 。 空气压缩热的 影响因素包括空气密度 、 空气流量 、 井筒深度等 , 由 于深井矿山井筒深度大 , 因此其空气压缩热成为影 响井下热环境的重要因素 。 2 . 4 水沟 、 管路散热量计算参数 水沟中的水主要包括地下涌水 , 充填 、 凿岩回水 等 。 水温及流量根据不同位置情况进行设定 , 井下 涌水在地温较高的地段中也是一个重要的热源 。 井 巷中的风水管路也会对热环境产生一定影响 , 其影 响因素包括管路断面尺寸 、 流量或流速 、流体温度 、 管路覆盖情况等 。 2 . 5 爆破作业及矿堆放热 炸药爆炸产生的能量一部分用来破坏矿岩结 构 , 另一部分则以热量的形式向矿内空气释放 , 同时 也使采下的矿石温度升高 。 在原岩温度已知的情况 下 , 岩石爆破放热量跟爆破后矿堆温度有关 , 影响因 素还包括炸药种类 、 炸药量 、 工作面尺寸 、工作面推 进速度 、 涌水情况 、 围岩情况 、 支护情况等 , 而爆破后 矿堆温度跟通风条件 、 矿石在井下停留时间 、 淋水等 因素有关 。 2 . 6 其他 其他热源包含矿石氧化 、木材氧化 、照明 、人体 放热等 , 一般情况下放热量较小 。 3 某深井矿山热环境模拟实例 3 . 1 矿山简介 瑞海金矿地处渤海海湾 , 位于三山岛北部海域 , 矿山生产规模为 12 000 t/d 。 矿床赋存在渤海海面 以下 , 主矿体埋藏深度超过 1 200 m , 设计采用上向 水平分层充填法开采 。 矿井需风量除满足向井下供 给足够新鲜风流 , 有效排出炮烟和粉尘外 , 还要排出 井下各种热源放出的热量 , 从而保证井下作业面的 温度满足安全规程的要求 。 使用 VUMA-3D 软件 对瑞海金矿的井下放热量进行了分析 , 对井下风流 、 空气温度 、 湿度等进行了模拟 。 3 . 2 输入参数 在 VUMA-3D 通风网络系统中 , 分支元素的类 型分为巷道(硐室)、 竖井 、 井底车场 、 采场 、 掘进工作 面 、 风机机站 、通风构筑物等 , 其中巷道(硐室)、竖 井 、 井底车场 、 采场 、 掘进工作面为可能的热源 , 在进 532 万方数据 中 国 矿 业第 27 卷 行热环境模拟时 , 需要对每个分支确定散热量计算 参数 , 进行热力学解算 。 这些参数大部分是可以量 化的 , 而对于那些难以量化的指标 , 按照分级的方法 进行确定 , 这些参数包括岩石表面湿润程度 、 管道绝 缘条件 、 柴油设备工作时间利用率 、热源湿度等级 、 采场涌水量等 。 VUMA-3D 对各个等级进行了描 述性说明 , 方便用户根据实际情况来选取 。 3 . 2 . 1 空气动力学参数 空气动力学模拟为热力学模拟的基础 , 主要巷 道的输入参数见表 1 。 表 1 主要巷道输入数据一览表 巷道名称 摩擦阻力系数/ (Ns2/m4) 直径/ m 断面/ m2 | 周长/ m 主井0 . 0356 . 331. 17219 . 792 副井0 . 0357 . 544. 15623 . 562 进风井0. 026 . 533. 16620 . 41 回风井0. 027 . 544. 15623 . 55 副井石门0 . 01213. 69214. 5 风井石门0 . 012277. 4719 . 678 分段巷道0 . 01213. 69214. 5 有轨运输巷道0 . 0128 . 53411 . 164 胶带斜井0. 0211. 96313 . 141 辅助斜坡道0 . 01213. 69214. 5 表 2 围岩热物理参数表 岩石热导率/ (W/(m K)) 比热/ (J/(kg ℃ )) 密度/ (kg/m3) 上盘围岩7. 9909 2 650 下盘围岩7. 9909 2 650 矿体8. 1909 2 790 表 3 主要巷道热力学参数表 巷道名称揭露时间/a湿润程度支护形式 主井6 Very wet砼/钢筋砼 副井6 Very wet砼/钢筋砼 进风井6 Very wet砼/钢筋砼 进风盲斜井6 Very wet喷砼/不支护 回风井6 Very wet砼/钢筋砼 倒段回风井5 Very wet砼 - 1 060 m 副井石门5 wet/damp喷砼/不支护 进风井石门5 wet/damp喷砼/不支护 回风井石门5 wet/damp喷砼/不支护 分段巷道4 wet/damp喷砼/不支护 有轨运输巷道4 wet/damp喷砼/不支护 胶带斜井4 wet/damp砼/钢筋砼 1 #斜坡道 4 ~ 5 wet/damp喷砼/不支护 2 #斜坡道 4 wet/damp喷砼/不支护 3 . 2 . 2 围岩散热参数 根据详查报告提供的矿区 5 个钻孔井温测量结 果 , 该区地温梯度在 1 . 45 ℃ /100 m ~ 1 . 94 ℃ /100 m 之间 , 1 412 m 深度温度为41 . 7 ℃ 。 在软件模拟中 , 设定地表原岩温度 20 ℃ , 地温梯度 0 . 015 ℃ /m 。 围岩热物理参数见表 2 , 主要巷道热力学参数见 表 3 。 3 . 2 . 3 设备参数 各柴油设备根据其使用地点情况(表 4) , 将其 散热量分配至相应巷道中 。 此外 , 还包括采区变电 所 、 牵引变电所中的变压器等设备(表 5) , 水泵房中 的水泵及排泥泵等设备 。 表 4 柴油设备输入参数表 设备 数量/ 台 单台功率/ kW 负载情况 工作环境 潮湿程度 柴油铲运机4186 Typical Mainly wet/ Intermediate 铲运机(掘进)7102 Typical Mainly wet/ Intermediate 服务车14 96Typical Mainly wet/ Intermediate 坑内卡车5130 Typical Mainly wet/ Intermediate 撬毛台车2120 Typical Mainly wet/ Intermediate 表 5 电动设备输入参数表 设备 数量/ 台 单台功率/ kW 负载情况 工作环境 潮湿程度 凿岩台车14 63Typical Mainly wet/ Intermediate 电动铲运机6132 Typical Mainly wet/ Intermediate 锚杆台车4115 Typical Mainly wet/ Intermediate 电移动空压机6132 Typical Mainly wet/ Intermediate 电机车10 30Typical Mainly wet/ Intermediate 风机42 450ATypical Mainly wet/ Intermediate 无轨维修设备1 套50Typical Mainly wet/ Intermediate 3 . 2 . 4 采场输入参数 VUMA-3D 系统针对不同的采场方法设置了 不同的计算参数 , 支持的采矿方法包括 进路法 、分 层充填法 、 分段空场法 、房柱法 、分段崩落法等 。 所 632 万方数据 增刊 2朱瑞军 ,等 基于 VUMA-3D 的深井矿山热环境模拟研究 需的计算参数包括采场结构参数 、 工作面推进速度 、 设备运行情况 、 涌水情况 、 围岩情况 、 支护情况等 。 本案例中矿山采用上向水平分层充填采矿法 , 以盘区为回采单元组织生产 。 盘区垂直走向布置 , 长 100 m , 宽度为矿体厚度 , 中段高度 60 m 。 每个 中段设 3 个分段 , 分段高度 20 m , 每条分段巷道承 担 6 ~ 7 个分层的回采工作 。 每个盘区内划分矿房 、 矿柱进行两步骤回采 。 矿房 、矿柱均垂直矿体走向 布置 , 宽度均为 10 m , 长度为矿体厚度 。 盘区生产 能力为 1 200 t/d , 10 个盘区同时生产 。 3 . 2 . 5 水沟 、 水管参数 按照井下正常涌水量 6 000 t/d 设置涌水参数 。 坑内凿岩 、 防尘用水点分布于回采和掘进工作面 , 总 用水量为 2 000 m 3 /d 。 主供水管敷设于副井井筒 内 , 由地表到达各中段 , 坑内供水管沿各中段或分段 敷设 , 接至各生产用水点 , 供水支管选用φ159 mm 8 mm无缝钢管 。 3 . 3 模拟结果 在 VUMA-3D 软件系统中建立通风网络三维 模型 , 并输入空气动力学及热力学参数后进行热环 境模拟 , 图 1 所示为井下温度分布情况 。 在进风空气的干球温度为 12 . 5 ℃ , 相对湿度为 80% 的条件下 , 井下热负荷计算结果见表 6 , 从模拟 结果中看出 , 井下放热量主要来自于空气压缩 、 围岩 散热 、 设备散热 、 采场爆破散热及地下水散热 。 局部 通风不畅地段会由于热量集中而导致温度较高 , 如 图 2 所示 。 在模拟结果中 , 可查看每段巷道的风速 、 风量 、 进风温度 、 出风温度 、 空气湿度 、 散热量等 。 图 1 井下温度分布情况图 图 2 温度较高地段示意图 表 6 井下热负荷计算结果 热源散热量/kW 围岩8 960 设备3 424 空气压缩9 711 裂隙水2 199 采场爆破等2 730 生产用水等吸热- 475˙ 其他295 合计26 844 在不同的进风温度情况下 , 井下发热量也不相 同 , 模拟了在大气压强 101 kPa 条件下 , 不同进风温 度和湿度情况下井下总散热量 , 结果见表 7 。 表 7 不同环境温度下的井下总散热量 地表进风参数模拟结果 干球温度/ ℃ 相对湿度/ % 湿球温度/ ℃ 焓值/ (kJ/kg) 含湿量/ (g/kg) 总散热量/ kW 工作面平均 干球温度/℃ 工作面平均 湿球温度/℃ 12 . 58010. 6330E. 877o. 2626 844J28p. 326f. 3 12 . 59011. 5833E. 198o. 1726 330J28p. 726f. 7 15 67 . 411. 6333E. 267o. 1926 331J28p. 726f. 7 15 8012. 9836Z. 78o. 5525 632J29p. 227f. 2 15 9014. 0139E. 459o. 6425 034J29p. 727f. 7 根据模拟结果 , 在进风空气的干球温度为12 .5 ℃ , 相对湿度 80% , 含湿量 7 . 26 g/kg , 焓值 30 . 87 kJ/kg 情况下 , 井下总发热量为 26 844 kW , 此时采掘工作 面的平均干球温度约为 28 ℃ , 达到安全规程所规定 的上限 。 因此 , 当进风空气的焓值大于此情况时 , 则 需要考虑采取制冷措施 。 从模拟结果中看出 , 在不 同的地表气象参数下 , 空气温度 、 湿度不同但焓值相 同时 , 总散热量和工作面温度几乎相同 。 4 结 语 利用 VUMA-3D 系统可有效地对深井矿山热 环境进行模拟 , 可掌握井下整体或局部地区的温度 、 湿度分布情况 , 并进行热源分析 , 方便对模拟矿井采 取针对性的制冷降温措施 。利用VUMA-3D系统 (下转第 262 页) 732 万方数据 中 国 矿 业第 27 卷 区内地下水补给 、 径流及排泄条件 , 严格受地形 地貌及岩性构造因素控制 , 具有典型的山地丘陵及 滨海平原区的特点 。 区内地下水化学动态及特征受 气象 、 水文 、 地形地貌 、 径流 、 排泄条件及人为因素的 控制 , 地下水水化学类型主要取决于其所处地理位 置 、 水文地质单元 、 补 、 径 、 排条件及补给区含水层岩 性 、 污染状况等多种因素 [5] 。 5 矿山地质条件 烟台市主要矿种矿床有金 、 煤 、铁 、铜 、钼钨 、水 泥用灰岩 、 水泥用大理岩 、滑石 、 菱美矿 、 长石 、 饰面 用花岗岩 , 以金矿最富盛名 , 黄金储量占全国的四分 之一 , 位居第一 [6] 。 矿产资源分布图见图 3 。 6 结 语 烟台市矿产资源丰富 , 多种矿产资源储量在全 国占有重要地位 , 其中金矿为优势矿产 , 储量位列全 国第一 。 由于成矿地质条件不同 , 全市矿种分布呈 现出明显的地域差异 。 龙口具有煤炭和油页岩资源 特色 ; 莱州拥有菱镁矿 、高岭土 、长石 、饰面花岗岩 、 天然卤水等资源优势 ; 金矿分布范围广而相对集中 , 查明的储量主要分布于招远 、 莱州 、 牟平三地 ; 有色 金属矿主要分布于栖霞 、福山 ;查明储量的铁矿 85% 分布于莱州 ; 查明储量的滑石 75% 分布于栖 霞 ; 查明的水泥用灰岩与大理岩主要分布于栖霞 、 福 山 、 蓬莱 ; 建筑用大理岩 、 花岗岩多分布于莱州 、 牟平 等地 ; 膨润土 、 沸石 、 珍珠岩分布于莱阳市 。 参考文献 [1 ] 高琴 , 于波 . 数字烟台地理空间基础框架建设初探[J] . 电脑知 识与技术 , 2007 , 3(13) 23-24 . [2 ] 孙媛媛 , 李爽 , 张祖陆 . 烟台市人居环境研究[J] . 云南地理环 境研究 , 2012 , 24(1) 75-80 . [ 3 ] 姜珊 , 于秀芹 , 李来民 , 等 . 1980 ~ 2010 年烟台市降水变化规 律及时空分布特征分析 [J] . 现代农业科技 , 2014 (3) 263- 265 . [ 4 ] 邵金花 , 刘贤赵 , 李德一 . 烟台水资源与社会经济可持续发展 协调度分析[J] . 经济地理 , 2007 , 27(4) 599-602 . [5 ] 王燕 , 王巍萍 , 王慧婷 . 烟台市水资源状况及对策研究[J] . 山 东水利 , 2007(4) 40-41 . [ 6 ] 张云吉 . 烟台市矿产资源的开发利用对策[J] . 资源开发与市 场 , 1998(4) 163-165 . (上接第 237 页) 进行热环境模拟需要掌握较全面矿山热源相关参 数 , 包括地温梯度 、 岩石热导率 、 地下水情况 、 不同巷 道的支护情况 、 围岩揭露时间 、设备运行状况等 , 输 入参数的准确性和全面性直接影响了热环境模拟的 准确性 。 深井矿山通风降温的能耗较高 , 应采取按需降 温的理念 , 对井下温度进行分区控制 , 最大限度降低 能耗 , 例如对于局部通风不畅地段 , 会由于热量集中 而导致温度较高 , 可采取喷雾洒水 、 加强局部通风等 措施 。 参考文献 [1 ] 连海瑛 , 严鹏 , 冯福康 , 等 . 二道沟矿深井开采热环境研究[J] . 黄金 , 2012 , 33(6) 27-30 . [2 ] 赵小稚 , 崔嵛 , 王敬志 . 基于 Ventsim 的深井采矿热环境分析 [J] . 铜业工程 , 2014(3) 52-54 . [ 3 ] 李丽峰 . 矿井热环境预测模型研究及应用[D] . 阜新 辽宁工 程技术大学 , 2013 . [ 4 ] 刘成敏 . 深井铜矿山热环境分析与降温方法研究[D] . 赣州 江西理工大学 , 2015 . [ 5 ] 张习军 , 王长元 , 姬建虎 . 矿井热害治理技术及其发展现状 [J] . 煤矿安全 , 2009(3) 33-37 . [6 ] 姚韦靖 , 庞建勇 . 我国深部矿井热环境研究现状与进展[J] . 矿 业安全与环保 , 2018(1) 113-117 . [ 7 ] 王君恒 , 李淑玲 , 姚长利 . 固体地球物理学概论[M ] . 北京 地 质出版社 , 2014 . [8 ] 马俊生 , 张爱民 , 叶勇 . 深井矿山井下热负荷计算研究[J] . 中 国矿山工程 , 2018(2) 8-12 , 20 . (上接第 257 页) 5 结 论 本文采用 3DMine 矿业工程软件 , 以国内某金 矿地质数据为基础 , 建立了矿体模型以及井巷模型 , 完成了矿体品位估值 , 探讨了软件在采矿设计中的 应用 , 通过以上工作为智能矿山建设提供了数据和 软件基础 。 参考文献 [ 1 ] 吴健飞 , 叶义成 , 王其虎 , 等 . 某多层复杂矿床开采优化的三维 地质建模[J] . 金属矿山 , 2012(9) 124-128 . [2 ] 柳波 , 陈广平 . 基于 3DMine 的贾家堡铁矿储量估算研究[J] . 现代矿业 , 2010 , 26(11) 1-4 . [ 3 ] 朱汉朝 , 张小友 . 赵振江 . 3DMine 软件的矿床模型构建及矿 石质量估算[J] . 科技创新与应用 , 2016(21) 46-47 . 262 万方数据
展开阅读全文

资源标签

最新标签

长按识别或保存二维码,关注学链未来公众号

copyright@ 2019-2020“矿业文库”网

矿业文库合伙人QQ群 30735420