自燃钼镍矿矿渣内部温度分布特征及拟合研究.pdf

Series No.566August 2023 金 属 矿 山METAL MINE 总 第566 期2023 年第 8 期收稿日期 2023-04-24基金项目 贵州省地矿局地质科研项目(编号:黔地矿科201904 号,202220 号);遵义市地质环境治理恢复工程项目(编号:GZBS-2022-03)。作者简介 江攀和(1994),男,工程师,注册岩土工程师。通信作者 梁劲松(1982),男,高级工程师。自燃钼镍矿矿渣内部温度分布特征及拟合研究江攀和 陈 磊 张明思 唐 广 魏泽权 梁劲松(贵州省地矿局一 0 二地质大队,贵州 遵义 563000)摘 要 利用 K 型高温热电偶测试了自燃钼镍矿矿渣不同深度的温度,基于反距离权重法(IDW)和统计学方法分别绘制了不同深度温度平面分布图及竖向温度曲线,分析了水平及竖向温度分布特征,查明了高温区分布及自燃趋势性最大的区域,初步确定了火源中心;建立了竖向温度一元三次函数拟合模型,根据相关系数、方差分析及测量误差理论,确定了研究区竖向温度最优拟合模型,采用 ZK5 实测温度值进行了验证。以上研究结果为研究区自燃钼镍矿矿渣的修复治理提供了技术依据。关键词 自燃 钼镍矿矿渣 温度分布 拟合 中图分类号TD75+.1,X753 文献标志码A 文章编号1001-1250(2023)-08-291-06DOI 10.19614/ki.jsks.202308038Study on Temperature Distribution Characteristics and Fitting of Molybdenum-nickel Slag in Spontaneous CombustionJIANG Panhe CHEN Lei ZHANG Mingsi TANG Guang WEI Zequan LIANG Jinsong(102 Geological Brigade,Guizhou Bureau of Geology and Mineral Resources,Zunyi 563000,China)Abstract The temperature of spontaneous combustion molybdenum-nickel slag at different depths was measured by K-type high temperature thermocouple.Based on IDW and statistical methods,the plane and vertical temperature curves of tem-perature at different depths were drawn respectively.The horizontal and vertical temperature distribution characteristics were analyzed,and the high temperature distribution and the largest spontaneous combustion trend area were identified.Based on the correlation coefficient,analysis of variance and measurement error theory,the optimal fitting model of vertical temperature in the study area was established,and verified by the measured temperature of ZK5.The above results provide a technical basis for the remediation of moly nickel slag in the study area.Keywords spontaneous combustion,molybdenum-nickel slag,temperature distribution,fit 镍钼矿属黑色岩系多金属共生矿,成分较复杂,含有多种金属及非金属元素1-3,主要分布在贵州遵义、湘西北、江西都昌、云南和浙江富阳等地。钼镍矿矿层主要产于寒武系下统牛蹄塘组底部黑色炭质泥页岩中,具有明显的硫化金属矿特征,矿体及附近常发育伴生黄铁矿,这种矿藏由于有机碳含量高而被称为“石煤”4-10。钼镍矿开采方式主要为井下开采,由于矿层薄,采掘比大,导致开采过程中钼镍矿矿渣较多,矿渣的发热量有 4 606.58 kJ,是标煤发热量的 15.73%,热值较高;以往不规范的开采,导致大量的钼镍矿矿渣露天堆积于斜坡或冲沟中,受到风吹雨淋,含有黄铁矿的钼镍矿矿渣极易发生氧化放热反应,使矿渣内部温度升高最终导致矿渣自燃,自燃的钼镍矿矿渣将会释放大量 H2S、SO2、CO 等有毒气体,威胁矿渣堆周边民众的生命健康安全11-12。杜永吉等13以 280 煤矸石山为例,研究得出自燃煤矸石山温度水平分布无规律性,垂直方向由表及里先升高后降低;陈兵等14基于一元二次函数和指数函数的曲线拟合,建立了温度拟合模型,揭示出自燃煤矸石山内部温度的变化趋势;杜玉玺等15基于一元三次函数曲线拟合,建立了煤矸石山深部温度拟合最优模型,实现了自燃煤矸石山的火源定位及温度分布可视化;杨娜等16实测了煤矸石山不同深度的温度,建立了浅深层温度间的一元二次拟合模型,发现高温区主要分布在靠近边坡的一侧,表面 6.0 m 处192自燃倾向性最高;杨娜17对荫营煤矿自燃煤矸石成分进行了测试分析,煤矸石含有 Si、Al、Fe、Ca、Mg、S等氧化物和 C,主要成分是氧化钙、二氧化硅以及氧化铁。由于我国目前对钼镍矿矿渣内部自燃温度场的相关研究较少,经采集研究区不同自燃状态的矿渣样化验分析,其矿物成分及元素含量与煤矸石主成分相近,同时含有提供自燃初期热量的黄铁矿。本文将借鉴煤矸石温度场研究相关内容,通过实测钼镍矿矿渣地表及内部温度,分析内部温度分布特征,基于最小二乘法建立竖向温度拟合模型,划定自燃趋势性较大的区域及初步确定火源位置,为后续钼镍矿矿渣自燃的防治提供技术依据。1 研究区背景研究区位于遵义市汇川区某钼镍金属矿矿渣堆北侧矿渣区域,该矿渣堆 2007 年开始堆填,2011 年停止堆填,纵向长约 150 m,横向宽约 60 m,总面积为7 420 m2,其中平台面积为 3 597 m2,斜坡面积 3 823 m2,斜坡平均坡度约 26。由于地形条件限制,钼镍矿矿渣自然堆放于冲沟中,平台区轻度压实,斜坡区未经压实,这种堆积方式使得斜坡区矿渣堆受氧面积增大,容易发生自燃。2010 年在研究区发现钼镍矿矿渣存在自燃现象,矿渣堆表面出现臭鸡蛋气味的黑烟,2019 年底对该矿渣堆采取了覆盖法及注浆法处理,自燃基本得到控制,目前斜坡区及靠近斜坡的平台区域矿渣已完成自燃,但平台区内部矿渣温度仍在升高,依然存在自燃的可能性。2 现场温度数据采集方法及结果为真实准确地测量和预测钼镍矿矿渣内部温度,在矿渣平台区、斜坡区及其交界处共布设了 10 个钻孔,钻孔位置分布见图 1。温度测量仪器主要使用 K型高温热电偶,由铠装线连接温度传感器及就地温度显示仪,测量温度范围为-101 000。各孔深度有所不同,温度测量间隔均为 1 m,本次实验有 ZK1 至ZK10 共 10 个测温点,首先测量每个测温点表层温度(0.1 m,矿渣表层平均覆土厚度约 0.5 m),为保证孔内温度数据测量的准确性,采用干钻法机械成孔,每钻进一米,间隔 30 min(降低干钻法钻头摩擦升温的影响),然后将温度传感器置入孔内测温点,停留约30 s 待温度数值稳定后读数,温度测量结果见表 1。已完全自燃矿渣岩芯呈暗红色或紫红色,自燃初期矿渣岩芯呈黑色。图 1 钻孔位置分布Fig.1 Distribution of boreholes表 1 钼镍矿矿渣温度测量数据Table 1 Temperature measurement data of molybdenum-nickel slag钻孔编号不同钻孔测温深度下测点温度/0.1 m1 m2 m3 m4 m5 m6 m7 m8 m9 m10 m11 m12 m13 m14 mZK1133841475457565762636876ZK21959110141168186194220208ZK3203349646584ZK420435570687383120110ZK5194351556878ZK61829488287ZK71719495361758687ZK8162136464853556164697175797370ZK9184258596159605859ZK102225272838404044485370725351 注:表示该孔矿渣处于自燃初期;表示该孔矿渣已完全自燃。3 温度分布特征3.1 水平方向温度分布特征温度是衡量钼镍矿矿渣自燃的重要指标,为了研究温度水平分布规律及划定自燃趋势性较大的区域,选择了测温点较多的 0.18.0 m 深度范围内的温度数据,经统计分析,0.18.0 m 温度数据变异系数分292总第 566 期 金 属 矿 山 2023 年第 8 期别为 13.66%、34.84%、42.41%、47.43%、50.47%、54.54%、63.40%、66.75%、66.22%,一般变异系数小于 15%称为弱变异,15%36%为中等变异,大于36%为强变异18,反距离权重法(IDW)对较强变异性的数据插值效果最佳19。本文对钼镍矿矿渣每一层深度的散点温度利用 ArcGis10.2 空间分析中的反距离权重法(IDW)进行插值分析,得到了每一层钼镍矿矿渣的温度平面分布图(图 2),该图可以直观地显示出水平层的温度分布特征及温度变化趋势。由图 2 可知,整体而言,温度变化趋势趋于一致,图 2 自燃钼镍矿矿渣各深度温度平面分布Fig.2 Temperature distribution at different depths of molybdenum nickel slag in spontaneous combustion均以 ZK2 为中心向周边发散。钼镍矿矿渣各层温度高温区域主要分布在研究区中部平台区 ZK2 及附近,研究区北侧斜坡区温度明显较低。经现场调查,研究区北侧为斜坡区,亦为迎风面,与平台区相比,有两个面(斜坡面及坡顶平面)的氧气供给,矿渣中的黄铁矿可以充分氧化放热,目前斜坡区及坡顶交界处矿渣已发生自燃,现处于自燃温度衰减末期,故温度较低;研究区中部为平台区,钼镍矿矿渣上覆黏土平均厚度约 0.5 m,具备一定的聚热条件,矿渣中黄铁矿不断氧化放热,由于氧气供给不充分(氧气供给充足约是供给不充足时释放热量的 3.5 倍),致使平台区中部温度升温较缓,尚未达到碳类物质的燃点 280,目前所测最高温度为 ZK2 孔中 7.0 m 深度处的220。现目前以 ZK2 为中心,四周外扩约 10 m 范392 江攀和等:自燃钼镍矿矿渣内部温度分布特征及拟合研究 2023 年第 8 期围内钼镍矿矿渣的自燃趋势性最大。3.2 竖直方向温度分布特征为研究钼镍矿矿渣内部温度竖向分布特征,选取测量深度较大的 6 个钻孔绘制温度曲线(图 3)。由图 3 可知,在 0.16.0 m 深度范围内,大部分测点温度随着深度变化趋势大体相同,竖向温度变化非线性,随着深度增加先快速升高后缓慢增高或趋于平缓。在 0.12.0 m 深度范围内温度快速升高,温度高的监测点与温度低的监测点相比,升温变化率较大,2.0 m 以后温度升温率逐渐降低。在接近矿渣与下伏基岩界面之上约 1.0 m 位置(所有钻孔均钻进至矿渣下伏炭质泥岩地层 0.51.5 m),出现了 2 种情况:一是在 2.0 m 以后,矿渣与下伏基岩界面之上约 1.0 m 位置,温度达到最高,随着深度继续增加,温度开始逐渐降低,如 ZK2、ZK4、ZK8、ZK10;二是在2.0 m 以后,随着深度增加,温度升温速率逐渐降低后温度趋于缓慢增加,如 ZK1、ZK9。图 3 自燃钼镍矿矿渣各深度温度曲线Fig.3 Temperature curves at different depths of moly-bdenum nickel slag in spontaneous combustion经现场调查及勘查,ZK2、ZK4、ZK8、ZK10出现温度先增加后降低的现象,说明钼镍矿矿渣在竖直方向上具有火源中心的特征,火源中心主要集中在矿渣与下伏基岩界面之上约 1.0 m 位置,表层不利于聚热,且越接近地表受到大气温度的影响越大,深层又缺乏氧气,不利于钼镍矿矿渣中黄铁矿的氧化放热或自燃,虽 ZK8、ZK10 比 ZK2、ZK4 的火源中心位置深度大 4.05.0 m,主要原因是 ZK8、ZK10 分别为平台区边界孔与边坡孔,虽然矿渣厚度较大,但是深部氧气供给与平台区内部相比较充足。ZK1、ZK9 出现温度缓慢增加的现象,可能是 ZK1、ZK9 位于平台区内部,聚热条件较好,但是供氧条件不足,导致钼镍矿矿渣中的黄铁矿一直处于初期缓慢氧化放热阶段,无明显的高温区或燃点,故整体温度差异不明显,预计温度随着深度的增加受到氧气浓度的限制而降低。研究区温度最高点位于 ZK2 地下 7.0 m 深度处,为 220;由于斜坡区及斜坡区与平台区交界处属于迎风面,氧气供给相对充足,目前该区域内钼镍矿矿渣已自燃。4 竖向温度曲线拟合4.1 拟合模型的建立为了直观地了解竖向温度的变化趋势,同时对深部温度进行预测,本文在最小二乘法20的基础上,运用 EXCEL 对研究区内 9 个钻孔温度数据基于指数模型和多项式模型进行拟合(表 2),保留 ZK5 作为最优拟合模型的验证孔。模型拟合优度用相关系数 R2进行判断,R2越接近 1,说明各个测温点越靠近拟合曲线,残差越小,拟合效果越好。表 2 结果表明,一元三次函数拟合模型要明显优于指数函数拟合模型。表 2 竖向温度拟合曲线结果Table 2 Fitting curve results of vertical temperature钻孔编号一元三次函数拟合拟合模型R2指数函数拟合拟合模型R2ZK1y=0.142 7x3-2.698 9x2+17.9x+16.2160.976 8y=27.914e0.102 1x0.634 1ZK2y=0.018 5x3-3.514 4x2+51.208x+17.1620.993 8y=44.014e0.248 9x0.723 0ZK3y=0.333 3x3-3.035 7x2+19.56x+18.9290.978 8y=24.075e0.270 7x0.920 2ZK4y=0.277 8x3-3.456 7x2+21.876x+22.1820.927 0y=31.228e0.179 4x0.832 5ZK6y=-3.083 3x3+18.714x2-8.273 8x+18.7290.990 3y=19.340e0.419 0 x0.958 9ZK7y=-0.128 8x3+0.703 5x2+11.845x+14.7270.966 4y=20.568e0.242 8x0.854 8ZK8y=-0.000 6x3-0.360 7x2+9.191 1x+16.2850.978 1y=27.096e0.091 2x0.745 3ZK9y=0.367x3-5.826 1x2+28.262x+18.8280.985 4y=33.919e0.096 4x0.438 1ZK10y=-0.080 4x3+1.397 9x2-2.385 7x+24.9010.871 4y=24.284e0.080 8x0.837 44.2 模型显著性分析模型显著性分析可以检验自变量与因变量之间的函数关系在总体上是否显著成立,常用的统计分析方法有 t 检验、卡方检验和方差分析等,本研究的模型显著性分析采用方差分析法,利用 EXCEL 表中的数据分析工具,主要通过各个数据之间显示的偏差与F 分布函数统计数据中认为是属于误差范围内的偏差进行比较,来测验各组数据之间是否存在显著性差492总第 566 期 金 属 矿 山 2023 年第 8 期异14,若统计量 FF 临界值,则模型回归效果显著,反之,不显著,统计分析结果见表 3。表 3 方差分析法模型显著性检验结果Table 3 Significance test results of variance analysis model钻孔编号自由度FF 临界值ZK12394.011 27.945 4ZK21736.913 58.530 9ZK31127.302 510.044 3ZK41741.957 88.530 9ZK6913.466 111.258 6ZK71529.333 58.861 6ZK82990.329 97.635 6ZK917101.747 08.530 9ZK102774.490 47.721 2 由表 3 可知,在选取显著性水平 a=0.01 的情况下,采用方差分析法进行显著性检验,结果显示 F 值远大于 1,说明组间差异具备统计学意义,同时 F 值均不同程度大于 F 临界值,因此拟合模型结果是可信的。4.3 最优拟合模型的选择与验证从竖向温度曲线拟合模型得知,所有点的最优拟合模型均是基于一元三次函数,根据测量误差理论,将每一个钻孔的实测温度值分别与 9 个钻孔的一元三次函数拟合模型预测值相比取差值,计算其偏差值的平方和,取平方和最小的拟合模型为研究区最优拟合模型,计算结果表明,基于一元三次函数的最优拟合模型为 y=0.333 3x3-3.035 7x2+19.56x+18.929,R2=0.978 8。通过对 ZK5 实测温度值与最优拟合模型对 ZK5的预测温度值进行统计分析,结果见表 4。结果表明,预测和实测温度值偏差较小,最优拟合模型可信度较高。表 4 最优拟合模型验证数据统计Table 4 Statistics of optimal fitting model validation data钻孔深度/mZK5 预测值/ZK5 实测值/0.118.9319135.7943248.5751359.2955469.9368582.50785 结 论(1)基于 IDW,对钼镍矿矿渣不同深度的温度数据进行插值处理,得到了不同深度温度平面分布图,整体而言,温度变化趋于一致,高温区主要分布在研究区中部平台区 ZK2 及附近,研究区北侧斜坡温度明显较低。目前以 ZK2 为中心向四周外扩约 10 m范围内的钼镍矿矿渣自燃趋势性最大。(2)在竖直方向上,钼镍矿矿渣受到聚热及供氧条件的影响,在 0.12.0 m 深度范围内温度升温率较高,2.0 m 以后温度升温率逐渐降低,一是 2.0 m以后,在矿渣与下伏基岩界面之上约 1.0 m 位置,温度达到最高,随着深度继续增加,温度开始逐渐降低,初步确定其为火源中心位置,二是在 2.0 m 以后,随着深度增加,温度升温速率逐渐降低后温度趋于缓慢增加。(3)一元三次函数模型比指数函数模型的拟合效果更佳,回归效果显著,可以作为研究区内竖向温度预测的经验回归公式。根据相关系数及测量误差理论,确定研究区竖向温度最优模型为 y=0.333 3x3-3.035 7x2+19.56x+18.929,R2=0.978 8。ZK5 实测温度值与预测温度值偏差较小,拟合模型可信度较高。参 考 文 献1 Beate ORBERGER,Jan PASAVA,Jean Paul GALLIEN,et al.Se,As,Mo,Ag,Cd,In,Sb,Pt,Au,Tl,Re traces in biogenic and abio-genic sulfides from Black Shales(Selwyn Basin,Yukon territories,Canada):A nuclear microprobe studyJ.Nuclear Inst and Methods in Physics Research,B,2003,210(1):441-448.2 ORBERGER B,VYMAZALOVA A,WAGNER C,et al.Biogenic ori-gin of intergrown Mo-sulphide and carbonaceous matter in Lower Cambrian black shales Zunji Formation,S ChinaJ.Chemical Ge-ology,2007,238(3-4):213-231.3 刘俊场,杨大锦,李怀仁,等.钼镍矿综合利用工艺现状及趋势J.云南冶金,2013,42(5):31-36.LIU Junchang,YANG Dajin,LI Huairen,et al.Comperehensive utili-zation process present situation and trend of Ni and Mo depositJ.Yunnan Metallurgy,2013,42(5):31-36.4 李坤,刘波,卢贸达.贵州纳雍地区牛蹄塘组钼镍多金属矿床地质特征及找矿前景J.中国钼业,2018,42(4):24-30.LI Kun,LIU Bo,LU Maoda.Geological characteristic and prospec-ting perspective of Mo-Ni polymetallic deposit host in niutitang for-mation,in Nayong of Guizhou provinceJ.China Molybdenum In-dustry,2018,42(4):24-30.5 杜小全,赵远由,张志强.贵州遵义地区钼镍矿规模远景及开发中存在的问题J.矿产勘查,2011,2(4):376-382.DU Xiaoquan,ZHAO Yuanyou,ZHANG Zhiqiang.Scale prospect and the problems existing in the development of Ni-Mo deposit in Zhunyi area,GuizhouJ.Mineral Exploration,2011,2(4):376-382.6 杨光龙,崔登伟.黔北牛蹄塘组黑色岩系中的矿产特征及找矿远景分析J.贵州地质,2010,27(3):208-212.YANG Guanglong,CUI Dengwei.Mineral character and prospecting potential in the black rock measures of Niutitang formation in north GuizhouJ.Guizhou Geology,2010,27(3):208-212.592 江攀和等:自燃钼镍矿矿渣内部温度分布特征及拟合研究 2023 年第 8 期7 杨剑.黔北地区下寒武统黑色岩系形成环境与地球化学研究D.西安:长安大学,2009.YANG Jian.Study on the Formation Environment and Geochemistry of Lower Cambrian Black Shale Series,Northern Guizhou Province,ChinaD.Xian:Changan University,2009.8 杨剑,易发成.黔北下寒武统黑色岩系元素赋存状态及富集模式J.矿物学报,2012,32(2):281-287.YANG Jian,YI Facheng.A study on occurrence modes and enrich-ment patterns of Lower Cambrian black shale series in Northern Guizhou Province,China J.Acta Mineralogica Sinica,2012,32(2):281-287.9 魏怀瑞,杨瑞东.贵州天柱早寒武世黑色岩系自燃现象的发现J.地质论评,2009,55(2):189-190.WEI Huairui,YANG Ruidong.Discovery of spontaneous combustion of Early Cambrian black rock series in Tianzhu,GuizhouJ.Geo-logical Review,2009,55(2):189-190.10 董允杰,缪加坦.我国钼镍矿及生产现状J.中国钼业,2008(2):60.DONG Yunjie,MIAO Jiatan.Molybdenum nickel ore and produc-tion status in ChinaJ.China Molybdenum Industry,2008(2):60.11 孙亚莉,王砚耕.贵州钼镍钒矿资源开发及其环境问题J.矿床地质,2012(S1):951-952.SUN Yali,WANG Yangeng.Molybdenum nickel vanadium resource development and environmental problems in GuizhouJ.Mineral Deposits,2012(S1):951-952.12 梁成,卢益平,程胜高.遵义松林镇冉村沟钼镍矿矿渣化学成分分析J.环境污染与防治,2015,37(8):124-124.LIANG Cheng,LU Yiping,CHENG Shenggao.Chemical composi-tion analysis of slag from Rancun Gou Molybdenum Nickel Mine in Songlin Town,ZunyiJ.Environmental Pollution&Control,2015,37(8):124-124.13 杜永吉,张成梁,李温雯.自燃煤矸石山温度特征研究J.矿业研究与开发,2010(1):92-95.DU Yongji,ZHANG Chengliang,LI Wenwen.Study on temperature characteristics of spontaneous combustion of gangue pilesJ.Min-ing Research and Development,2010(1):92-95.14 陈兵,胡振琪,赵艳玲,等.自燃煤矸石山内部温度拟合研究J.能源环境保护,2011,25(2):17-20,24.CHEN Bing,HU Zhenqi,ZHAO Yanling,et al.Temperature fitting study of combustion spontaneous coal waste pilesJ.Energy Envi-ronmental Protection,2011,25(2):17-20,24.15 杜玉玺,苏未曰,吴媛婧,等.自燃煤矸石山内部温度拟合与可视化研究J.矿业安全与环保,2018,45(5):32-36.DU Yuxi,SU Weiyue,WU Yuanjing,et al.Study on fitting and vi-sualization of internal temperature of spontaneous combustion gangue hillJ.Mining Safety&Environmental Protection,2018,45(5):32-36.16 杨娜,张永波,牛金荣.荫营煤矿自燃矸石山温度场分布及深部温度拟合J.矿业安全与环保,2021,48(4):23-27,32.YANG Na,ZHANG Yongbo,NIU Jinrong.Temperature field distri-bution and deep temperature fitting of spontaneous combustion gangue dump in Yinying Coal MineJ.Mining Safety&Environ-mental Protection,2021,48(4):23-27,32.17 杨娜.荫营煤矿自燃矸石山温度分布及动态演变规律研究D.太原:太原理工大学,2021.YANG Na.Study on Temperature Distribution and Dynamic Evolu-tion Law of Spontaneous Combustion Gangue Hill in Yinying Coal MineD.Taiyuan:Taiyuan University of Technology,2021.18 范俊楠,贺小敏,陆泗进,等.湖北省重点行业企业周边土壤重金属污染现状及潜在生态危害评价J.华中农业大学学报,2018,37(5):74-80.FAN Junnan,HE Xiaomin,LU Sijin,et al.Status of heavy metal pollution and potential ecological hazard assessment in surrounding soil of enterprises in Hubei Province,ChinaJ.Journal of Hua-zhong Agricultural University,2018,37(5):74-80.19 陈思萱,邹滨,汤景文.空间插值方法对土壤重金属污染格局识别的影响J.测绘科学,2015,40(1):63-67.CHEN Sixuan,ZOU Bin,TANG Jingwen.Impact of spatial interpo-lation methods on identifying structure of heavy metal polluted soilJ.Science of Surveying and Mapping,2015,40(1):63-67.20 高秋英,王丽丽,王荣忠.最小二乘法曲线拟合及优化算法研究J.工业控制计算机,2021,34(11):100-101.GAO Qiuying,WANG Lili,WANG Rongzhong.Research on least square curve fitting and optimization algorithmJ.Industrial Con-trol Computer,2021,34(11):100-101.692总第 566 期 金 属 矿 山 2023 年第 8 期