资源描述:
收稿日期2020-02-24 作者简介王明 (1994) , 男, 硕士研究生。 某矿山硫化矿堆自燃过程温度场时效分析 王明张建华叶永喜黄刚 1 (武汉理工大学资源与环境学院, 湖北 武汉 430070) 摘要针对某铜矿山硫化矿石的堆放, 解决其硫化矿堆自燃的安全问题。考虑堆积体内部硫化矿石自身的氧 化生成热、 绝热温升2种升温模式, 利用ANSYS建立矿山常用的靠侧壁三角形和平坦地基上梯形2种堆积体的温度场 模型, 研究2种堆积方式同体积下硫化矿石堆内部温度场及自燃区域随时间的变化规律。通过计算模拟, 获得堆积体 内部不同时刻温度场分布图、 自燃区域分布图。研究结果表明采取防治措施应在日升温为6 ℃以前的阶段; 相同体积 下表面积更大的矿堆形状, 内部升温速度更慢、 自燃时间更短。根据研究结果推荐采用梯形堆放方式, 以及相关防治 措施。 关键词硫化矿堆自燃温度场升温曲线自燃区域 中图分类号TD714文献标志码A文章编号1001-1250 (2020) -06-210-06 DOI10.19614/ki.jsks.202006033 Time-effect Analysis of Temperature Field during Spontaneous Combustion of Sulfide Ore Reactor in a Mine Wang MingZhang JianhuaYe YongxiHuang Gang2 (School of Resources and Environmental Engineering, Wuhan University of Technology, Wuhan 430070, China) AbstractIn view of sulfide ore stacking in a copper mine, the safety problem on spontaneous combustion of sulfide ore stack should be solved. Considering the two heating modes of oxidation heat and adiabatic temperature rise of sulfide ore in the deposits,and with use of ANSYS,the temperature field models for two kinds of deposits are established,that is triangle near the side wall and trapezoid on flat foundation,which are commonly used in mines,to study the variation of temperature field and spontaneous combustion zone with time in the same volume of sulfide ore stack under two kinds of deposit modes. Through the calculation and simulation, the temperature field and the spontaneous combustion area distribution map at differ⁃ ent time inside the deposit are obtained. The results show that the prevention and control measures should be taken before the daily temperature rise of 6 ℃. The shape of ore heap with larger surface area under the same volume has slower internal tem⁃ perature rise and shorter spontaneous combustion time. According to the research results,the trapezoidal stacking mode and its related prevention measures are recommended. KeywordsSpontaneous combustion, Temperature field, Heating curve, Spontaneous combustion area 硫化矿石的自燃发火问题是开采含硫矿床的矿 山企业生产中需要预防和可能遭受的灾害问题。通 过相关数据统计分析发现, 国内具有含硫矿石自燃 安全隐患的金属矿山在所有矿山企业中占比高达 20, 其中大部分矿山都出现过不同程度的含硫矿石 自燃火灾事故 [1-3]。硫化矿石发生自燃火灾会导致矿 山企业的正常生产运作系统产生中断, 造成自燃区 域内大量矿石的损失, 并在氧化自燃时生成大量具 有腐蚀性的气体。因此, 自燃火灾给矿山企业造成 的巨额经济损失的同时产生众多不可忽视的安全隐 患及环境问题 [4-8]。 关于硫化矿石的自燃问题 [7-10], 国内外学者从多 种学科研究角度做了细致的探究。多种研究方法都 被研究人员用于硫化矿石自燃问题研究。例如统计 学在硫化矿石自燃危险性评价 [11-12]中被广泛应用, 饶 运章等 [13]建立硫化矿石堆氧化自热温度的GRNN神 经网络模型, 以含硫量、 矿石块度、 升温梯度作为试 验影响因素来研究硫化矿石堆氧化自热温升速率, 总第 528 期 2020 年第 6 期 金属矿山 METAL MINE Series No. 528 June 2020 210 ChaoXing 最终预测误差为3.51; 文虎等[14]将复杂的传热模式 利用有效导热系数法简化为矿石间的导热问题, 利 用COMSOL Multiphysics数值软件对硫化矿石堆孔隙 间具有不同流速的空气对硫化矿石堆温度场所产生 的影响规律进行数值仿真; 潘伟等 [15]运用 ANSYS、 MATLAB数值软件分别对硫化矿石堆自热过程、 采场 硫化矿石堆风流场分布情况及SO2、 O2浓度随时间的 变化趋势进行了模拟研究并建立了硫化矿石堆的二 维多孔介质模型。 上述研究对影响硫化矿石自燃的各种因素已经 做出较多论述, 但实际矿山生产过程中硫化矿石的 自燃问题最有效的解决方法, 是在矿石到达自燃点 之前升温较小阶段采取防治措施。故针对矿山实际 状况, 研究硫化矿堆温度场随时间的变化, 对于相似 矿山的生产安全工作具有指导意义。本研究根据矿 山实际生产状况建立2种最为常见的硫化矿堆模型, 利用ANSYS对其内部温度场的状况以及自燃区域进 行模拟, 研究矿堆升温曲线以及自燃区域动态变化, 探明堆积过程变化规律。 1工程概况 缅甸莱比塘铜矿是一座位于缅甸联邦西北部实 皆省南部的大型露天矿, 矿区面积约32.73 km2。矿 山年计划采剥量7 000万t, 台阶高度为15 m, 坡面角 70, 钻孔直径为250 mm。在L45-2088炮区施工时, 炮孔附近矿堆陆续出现冒烟现象, 导致炸药自燃, 炮 区装填重铵油炸药的炮孔共计41个, 发生炸药自燃 的共计20个。发生炸药自燃的区域属于含黄铁矿区 域; 发生炸药自燃炮孔与其他含黄铁矿区域炮孔对 比发现该区域硫化矿石Fe2、 S离子含量相差不大, 但 其pH值相对更低, 最小值为2.89。初步判定炸药自 燃原因为该爆破区域内矿堆的黄铁矿矿石接触空 气、 水发生预氧化或氧化反应释放出一定的热量, 导 致附近装药炮孔内的温度上升。当装入重铵油炸药 时, 热量积累, 导致炸药自燃、 炮孔冒烟。由此可见 该矿山硫化矿石堆其自燃倾向性很大, 模拟研究其 矿堆堆积时温度场变化规律对矿山生产有重要的意 义。 2岩石热传导理论 研究主体为多孔散体介质的硫化矿石堆, 且矿 石堆自身含有内热源。矿石堆内部热源为硫化矿石 自身发生氧化时所生成的热量。在硫化矿石堆放的 环境下, 其内部与周围堆积岩体及环境之间矿石材 料氧化热的等效热传导方程满足 [16] ∂T ∂t α ∂2T ∂x2 ∂2T ∂y2 ∂2T ∂z2 Q ˉ ρc ,(1) 式中,ρ为矿石材料密度;c为矿石材料比热容;T为外 界环境温度;t为时间; α λ ρc 为有效导热系数;λ为导 热系数;Q ˉ为单位时间单位体积矿石释放出的热量。 在式 (1) 中Q ˉ只考虑硫化矿石自身氧化时所释放 的热量, 即 Q ˉ ρc ∂θ t ∂t ρcθ0 ∂f t ∂t ,(2) 式中,θ t是硫化矿石堆积体的中心区域绝热温升负 指数函数方程;θ0为硫化矿石最终绝热升温时温度。 将式 (2) 代入式 (1) 中, 可得在仅考虑矿石材料氧化 热的等效热传导方程, 即 ∂T ∂t α ∂2T ∂x2 ∂2T ∂y2 ∂2T ∂z2 θ0 ∂f t ∂t .(3) 在此基础上则将硫化矿石堆积体温度场计算分 析问题转变为在设定的初始条件及边界条件下对函 数 (3) 求值的硫化矿石堆积体温度场计算分析问题 。 由设定的边界条件, 利用有限元数值分析法将式 (3) 所在求解区域内进行离散化处理, 最终得出该热 传导问题的有限元求解方程。 []K{ }T []C { } ∂T ∂τ { }F ,(4) []C ∫ Ω ρc[]N T[ ]N dΩ,(5) []K ∫ Ω[ ]B T[ ]D[]B dΩ ∫ Γ β[]N T[ ]N dΓ, (6) { }F ∫ Ω Q ˉ[ ]N TdΩ - ∫Γqq[ ]N TdΓ q∫ΓβTα[ ]N TdΓ, (7) []B ■ ■ ■ ■ ■ ■ ■ ■■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■ ■■ ■ ■ ■ ■ ■ ∂ ∂x ∂ ∂y ∂ ∂z []N ,(8) []D λI,(9) 式中,[]N []Ni,Nj,,Nr为单元型函数。 式 (4) 中对任意时间τ均成立,设 ΔTn Tn 1- Tn Δτn ■ ■ ■ ■ ■ ■ ■ ■ 1 - s ∂T ∂τ n s ∂T ∂τ n 1 . (10) 王明等 某矿山硫化矿堆自燃过程温度场时效分析2020年第6期 211 ChaoXing 可得 []K 1 sΔτn []C{}Tn 1 1 - s s []K - 1 sΔτn []C{ }Tn 1 - s s { }Fn{}Fn 1, (11) 式中,{ }Tn、{ }Fn、{}Fn 1均为已知量,{}Tn 1为未知 量, 所以, 式(11)为关于{}Tn 1的线性方程组, 求解 可得出各节点在ττn1时的温度{}Tn 1。 3模型建立 3. 1建立矿堆模型 3. 1. 1 硫化松散矿堆模型理论化假设 对研究矿山的硫化矿堆, 进行现场调研收集相 关资料。分析相关数据, 为了模拟的可实现性, 在 现场实际数据的基础上, 对建立的模型做如下假 设 ①所建立的硫化矿石堆模型的含水率为零; ②硫 化矿石的密度及比热容为某一固定常量; ③硫化矿 石材料均匀连续; ④在对模型进行温度场分析处理 时, 考虑堆积体内部硫化矿石自身的氧化生成热、 绝 热温升2种升温模式。 3. 1. 2建立模型 硫化矿石在采场中的堆积型式较为简单, 一般 状况下采用靠侧壁三角形和平坦地基上梯形2种堆 积型式, 在保证2种堆积形式其体积相同的前提下进 行比较分析。2种堆积形式几何模型如图1所示。 在进行数值仿真分析时将基于2种硫化矿石堆 积体的横截面建立二维模型对硫化矿石堆进行瞬态 热分析。其中, 2种模型中与基岩接触边界均设置为 温度类型, 设初始地基及硫化矿石堆温度设定为 27 ℃, 矿石自燃点为442 ℃; 对2种模型进行网格划 分时, 基岩及硫化矿石梯形堆积体均采用四边形网 格进行划分, 三角形堆积体则采用三角形网格划分, 2种模型网格划分如图2所示。利用有限元法对模型 进行计算求解, 进行2次计算。将第1次计算求解时 间步长设置为 86 400 s。第 2次计算调整 ANSYS计 算分析时间步长为7 200 s; 仅对不同堆积形式数值 模型中的硫化矿石堆部分的温度场重新进行数值计 算。利用第1次数值计算结果, 将三角形硫化矿石堆 积体从矿石堆积300 h时开始计算至矿石堆积480 h; 梯形硫化矿石堆积体从矿石堆积350 h时开始计算 至500 h。 3. 2定义材料参数 研究成果表明硫化矿石的氧化生热过程受到多 种因素的影响, 包括矿石中铁离子含量、 温度、 含水 量以及空气中氧含量等的影响, 复杂的变量使得学 者无法根据现有成果对某一类矿石给出合适的温升 曲线。故需对现场矿石进行勘察, 将所有的样品参 数进行处理。在进行模型的建立时, 忽略堆积形状 带来的一些模型参数的细微差别, 使梯形矿堆和三 角形矿堆模型材料参数保持一致, 构建该矿山常用 形状的硫化矿石堆数值模型时, 选用最常见的黄铁 矿石矿堆的堆积参数。模型的含水率为零, 比热容 及密度等参数为固定常数。 最终取矿石密度ρ0为2 700 kg/m3; 基岩密度ρ1为 2 700 kg/m3; 矿石堆密度ρK为2 310 kg/m3; 矿石比热容 金属矿山2020年第6期总第528期 212 ChaoXing c0为2.222 J/ (kg K);基岩比热容c1为0.7 J/ (kg K) ; 基 岩有效导热系数λh为4.0 W/ (m K) ; 基岩表面传热系 数β1为5.85 W/ (m2 K) ; 矿石堆孔隙率n为32; 有效 导热系数λh0.001 3 T0.037 3。内源热强度 qv ■ ■ ■ ■ ■■ ■ ■ ■■ ■ 838.96 - 807.95 1 T*97.26 8.03, T * ≤ 125 -2.93 T*- 384.84, T* 125 ; (12) 其中,T *的单位为℃。 其中, 矿石堆内热源强度是一个关于温度的函 数方程表达式; 将该函数用LM算法拟合成一个关于 时间的函数 [17], 硫化矿石内热源强度是一个关于时 间的分段函数, 分为自热低温氧化、 快速氧化、 持续 散热3个阶段。其函数关系如图3所示。 在硫化矿石氧化自热低温氧化阶段, 矿石氧化 热生成速率低, 此时内热源强度处于一个相对平稳 阶段; 在经过低温氧化阶段后矿石堆积聚了一定的 热量, 硫化矿石堆开始进入快速氧化阶段, 该阶段硫 化矿石氧化速率迅速提升, 其自热幅度大大增加, 当 达到矿石自燃点时, 硫化矿石便会发生自燃, 其内热 源强度显著升高; 在硫化矿石经过持续数天或更长 时间的燃烧后, 硫化矿石堆基本全部被氧化, 此时矿 堆放热速率迅速下降, 其内热源强度同样迅速降低。 4模拟结果分析 4. 1堆积体内部温度场仿真模拟过程分析 硫化矿石堆的自燃点, 取决于硫化矿堆内温度 最高点的大小, 故研究矿石堆的温度场时, 对矿石堆 内温度各个时刻最高点着重进行分析。利用ANSYS 后处理器提取出2种堆积形式的硫化矿石堆数值模 型各个时刻最高温度节点, 所得2种硫化矿石堆积体 内部最高温度随时间变化曲线如图4所示。 图4为60 d内不同堆积形式硫化矿石堆内部最 高温度随时间变化曲线, 2种堆积体的温度变化曲线 趋势大致相同, 变化规律如下 (1) 矿石堆积的前8 d内, 硫化矿石氧化速率较 小, 生成热较少, 温度变化平稳而缓慢, 堆积体内部 温度场分布无显著变化。2种堆积体的变化过程分 别为 ①三角形堆积体堆积前8 d内, 日均升温量为 6 ℃左右, 到第8 d后温度变化速率开始增加, 此时温 度场分布如图5 (a) ; ②梯形堆积体堆积前10 d内, 日 均升温量为2.5 ℃左右, 到第10 d后温度变化速率开 始增加, 此时温度场分布如图5 (b) 。 (2) 堆积第8 d至第20 d, 矿石氧化热量在内部不 断聚集提升温度, 氧化速率迅速提升, 放热速率不断 增大, 最终到达矿石自燃点, 发生自燃, 随后温度降 至自燃点下方。2种堆积体的变化过程分别为 ①三 角形堆积体从第8 d开始, 内部升温幅度变大, 第13 d 时达到自燃点, 到第16 d时内部温度开始下降, 当矿 石堆积体到第18 d其内部最高温度开始低于硫化矿 石自燃点, 期间硫化矿堆内部最高温度为578 ℃, 此 时堆积体内部温度场分布如图5 (c) ; ②梯形堆积体 从第10 d开始, 内部升温幅度变大, 第15 d时达到硫 化矿石的自燃点, 到第18 d时内部温度开始下降, 当 矿石堆积体到第19 d其内部最高温度开始低于硫化 矿石自燃点, 硫矿堆积体内部最高温度达到563℃, 此时堆积体内部温度场分布如图5 (d) 。 (3) 在此之后, 硫化矿石与氧气的化学反应进入 末期, 其放热速率逐渐降低, 使得硫化矿石堆内部温 度呈现下降趋势如图5 (e) 、 图5 (f) , 最终硫化矿堆放 热反应基本完成, 硫矿堆内部温度场不断与周围环 境发生热交换, 直至温度降至与周围环境温度大致 相同。 4. 2矿堆模型自燃区域变化规律。 由温度场数值模拟结果确定硫化矿堆的自燃过 程所经历的时间段, 对该时间段内堆积模型自燃区 域动态变化进行数值模拟, 利用积分公式求出该时 段内自燃区域的叠加面积变化曲线如图6所示。 由图 6 可知, 三角形硫矿堆自燃开始时刻为第 300 h, 其自燃面积占模型面积比例为 0.02; 在 456 h 王明等 某矿山硫化矿堆自燃过程温度场时效分析2020年第6期 213 ChaoXing [1] [2] [3] 时, 硫化矿石自燃面积所占模型总面积比例达到最 大值0.907, 如图7 (a) 。随后自燃区域叠加面积总占 比稳定在该数值, 硫化矿堆自燃持续时间约为156 h。 梯形硫矿堆在348 h时发生自燃, 其自燃面积占模型 面积比例为0.055; 在448 h时, 硫化矿石自燃面积所 占模型总面积比例达到最大值0.812, 如图7 (b) 。随 后自燃区域叠加面积总占比稳定在该数值; 硫化矿 堆自燃持续时间约为100 h。2种堆积体的自燃开始 区域均在矿堆中央。 5结论 利用 ANSYS 对某铜矿硫化矿石堆进行数值模 拟, 根据该矿山的矿石种类以及其他相关参数建立 模型计算后获得了该铜矿山硫化矿石堆积体在堆积 时间内的温度分布云图, 自燃区域面积大小及位置。 数值模拟分析结果表明, 矿山常用的2种堆积方式的 温度场变化规律大致相同, 采取自燃防治措施的最 佳时间在堆积前 8 d 之内; 自燃发生区域在矿堆中 心, 因而对长期堆放的矿石需进行中心部位温度定 期测定; 使用梯形堆积这种表面积大的堆积方式有 利于延缓硫化矿堆到达自燃点的时间, 降低自燃最 高温度, 减少自燃时间和自燃矿石损失量。 参考文献 郭虎. 硫化矿石自燃影响因素分析与实验研究 [D] . 重庆重 庆大学, 2016. Guo Hu. Analysis and Experimental Study on Factors Affecting Spontaneous Combustion of Sulfide Ore [D] . ChongqingChongqing University, 2016. 赵国彦, 古德生, 吴超. 硫化矿床内因火灾综合防治措施研究 [J] .矿业研究与开发, 2001 (1) 17-19. Zhao Guoyan,Gu Desheng,Wu Chao. Research on integrated fire prevention and control measures in sulfide deposits [J] . Mining Re⁃ search and Development, 2001 (1) 17-19. 胡汉华. 铜山铜矿采场防灭火试验研究 [J] .金属矿山, 2001 (5) 48-51. Hu Hanhua. Experimental study on stope fire prevention and extin⁃ 金属矿山2020年第6期总第528期 214 ChaoXing [4] [5] [6] [7] [8] [9] [10] [11] [12] [13] [14] [15] [16] [17] guishing in Tongshan Copper Mine [J] . Metal Mine,2001(5) 48- 51. 石振东. 硫对设备高温腐蚀产物自燃性及自燃预防措施的研究 [D] .沈阳 东北大学, 2009. Shi Zhendong. Study on the Spontaneous Combustion of High Tem⁃ perature Corrosion Products of Equipment by Sulfur and the Preven⁃ tion Measures of Spontaneous Combustion [D] .ShenyangNortheast⁃ ern University, 2009. 阳富强. 金属矿山硫化矿自然发火机理及其预测预报技术研究 [D] .长沙中南大学, 2011. Yang Fuqiang. Study on Spontaneous Combustion Mechanism and Prediction Technology of Sulfide Ore in Metal Mines[D] . Chang⁃ sha Central South University, 2011. 饶运章, 徐水太, 张见.含硫矿山酸性废水及其金属离子污染 防治方法的探讨 [J] .中国钨业, 2004, 19 (4) 33-35. Rao Yunzhang,Xu Shuitai,Zhang Jian. Discussion on the preven⁃ tion and control s of acid wastewater and metal ion pollution in sulfur-containing mines [J] . China Tungsten Industry,2004, 19 (4) 33-35. 王张辉. 含硫矿石自燃过程及机理的热重实验研究 [D] .西安 西安科技大学, 2010. Wang Zhanghui. Thermogravimetric Experimental Study on Sponta⁃ neous Combustion Process and Mechanism of Sulfur-Bearing Ores [D] . Xi′an Xi′an University of Science and Technology, 2010. 刘辉, 吴超, 阳富强, 等.红外热像技术探测硫化矿石自燃 火源的影响因素及其解决方法 [J] .科技导报, 2010, 28 (2) 91- 95. Liu Hui,Wu Chao,Yang Fuqiang,et al. Influence factors and so⁃ lutions of infrared thermography in detecting spontaneous combus⁃ tion of sulfide ore[J] . Science and Technology Herald,2010, 28 (2) 91-95. 王国利.句容硫铁矿矿石自燃倾向和炸药自爆危险性评价 [J] . 有色金属, 1999 (1) 8-10. Wang Guoli. Spontaneous combustion tendency of Jurong pyrite ore and risk assessment of explosive self explosion[J] . Nonferrous Metals, 1999(1) 8-10. 李孜军, 陈天丰, 李明, 等.化学泡沫阻化剂对硫化矿石自燃 阻化效果研究 [J] .中国安全科学报, 2017, 27 (10) 44-49. Li Zijun,Chen Tianfeng, Li Ming, et al. Study on sulfide ore spon⁃ taneous combustion inhibition effect of chemical foam inhibitor [J] . China Safety Science Journal, 2017, 27(10) 44-49. 蔡逸伦, 阳富强, 刘晓霞.硫化矿石自燃灾害预警的RBF神经网 络模型及应用 [J] .有色金属工程, 2019, 9 (7) 72-78. Cai Yilun,Yang Fuqiang,Liu Xiaoxia. RBF neural network model and application for early warning of sulfide ore spontaneous com⁃ bustion disaster[J] . Nonferrous Metals Engineering,2019, 9(7) 72-78. 阳富强, 吴超, 贾永军, 等. GM (1, 1) 模型在含硫矿堆温度场 预测中的应用 [J] .矿业快报, 2007 (6) 13-15. Yang Fuqiang,Wu Chao,Jia Yongjun,et al. Application of GM ( 1, 1 )Model in predicting temperature of sulphide ore piles [J] . Express Ination of Mining Industry, 2007 (6) 13-15. 饶运章, 袁博云, 吴卫强, 等. 基于GRNN模型的硫化矿石堆氧 化自热温度预测 [J] . 金属矿山, 2016 (6) 149-152. Rao Yunzhang, Yuan Boyun, Wu Weiqiang, et al. Prediction of self heating temperature of sulfide ore heap oxidation based on GRNN model[J] . Metal Mine, 2016 (6) 149-152. 文虎. 煤自燃过程的实验及数值模拟研究 [D] .西安西安科 技大学, 2003. Wen Hu. Experimental and Numerical Simulation of Coal Spontane⁃ ous Combustion Process[D] . Xi′an Xi′an University of Science and Technology, 2003. 潘伟. 硫化矿石堆自热过程的非线性及数值仿真研究 [D] .长 沙中南大学, 2011. Pan Wei. Nonlinear and Numerical Simulation of Self-Heating Pro⁃ cess of Sulfide Ore Heap[D] . Changsha Central South University, 2011. 赵军.不同堆积型式硫化矿石堆积体温度场数值模拟 [J] .中 国安全生产科学技术, 2012, 8 (5) 29-34. Zhao Jun. Numerical simulation of temperature field of sulfide ore deposits of different stacking types[J] . China Safety Production Science and Technology, 2012, 8(5) 29-34. 朱建芳, 郭文杰, 段嘉敏, 等. 基于活化能的煤自燃倾向性鉴定 方法探讨 [J] .华北科技学院学报, 2017, 14 (1) 6-10. Zhu Jianfang,Guo Wenjie,Duan Jiamin,et al. Discussion on the identification of coal spontaneous combustion tendency based on activation energy[J] . Journal of North China Institute of Science and Technology, 2017, 14(1) 6-10. 王明等 某矿山硫化矿堆自燃过程温度场时效分析2020年第6期 215 ChaoXing
展开阅读全文