工程科学学报，第38卷，第10期：1350-1358,2016年10月 Chinese Journal of Engineering,Vol.38,No.10:1350-1358,October 2016 D0l:10.13374/j.issn2095-9389.2016.10.002:http://journals..ustb.edu.cn 非均质孔隙率采空区氧化升温规律四维动态模拟 周佩玲”，张英华”，黄志安区，袁飞2》，高玉坤)，王辉”，孙倩” 1)北京科技大学金属矿山高效开采与安全教有部重点实验室，北京1000832)北京科技大学治金与生态工程学院，北京100083 3)北京科技大学钢铁治金新技术国家重点实验室，北京100083 ☒通信作者，E-mail:huang-za@qg.com 摘要运用Fuet动网格模型实现采空区的四维动态变化，并用用户自定义函数将煤低温氧化动力学机理及非均质孔隙 率函数编入Fut中，结合时间和空间，对U+L型通风系统采空区升温规律进行四维动态模拟研究.研究表明：非均质孔隙 率四维动态模型能更真实地反应孔隙率的空间与时间变化，空间某一位置的孔隙率随时间呈负指数递减：工作面推进速度越 大，采空区升温速率越小，推进速度为3.6m"d1时平均升温速率仅为推进速度为1.2m·d时的1/5：然而，推进速度越大，高 温点的深度越大，不利于自燃的预防：尾巷的存在使得温度场范围扩大，温度升高，C0主要从尾巷流出，尾巷释放的C0量是 回风巷C0释放量的10倍.最后利用现场实测的数据对结果进行验证，表明模拟结果是正确可信的. 关键词采空区：非均质孔隙率：煤燃烧：氧化升温：四维：动态模拟 分类号TD75·2 4D dynamic simulation of coal oxidation heating law in gobs with heterogeneous porosity ZHOU Pei--Hing》,ZHANG Ying-hua》,HUANG Zhi--ane,YUAN Fei2,3》,GA0Yu-hun',WANG Hui,SUN Qian” 1)Key Laboratory of High-Efficient Mining and Safety of Metal Mines (the Ministry of Education),University of Science and Technology Beijing, Beijing 100083,China 2)School of Metallurgical and Ecological Engineering,University of Science and Technology Beijing,Beijing 100083,China 3)State Key Laboratory of Advanced Metallurgy,University of Science and Technology Beijing,Beijing 100083,China Corresponding author,E-mail:huang_za@qq.com ABSTRACT Using the Fluent dynamic mesh model to realize the 4D dynamic movement of a gob and inputting the kinetic mecha- nism of coal low-temperature oxidation and the dynamic change of heterogeneous porosity into Fluent through the user defined function, 4D dynamic simulation is performed on the spontaneous heating law in a gob for U+L ventilation of a certain mine.The results show that unsteady heterogeneous porous media can be more realistic in response to the change of porosity,which exponentially decreases with time.The greater the advancing speed of the working face,the smaller the heating rate,and the average heating rate at an advan- cing speed of 3.6md is 1/5 of that at an advancing speed of 1.2m'd.However,the higher the advancing speed,the deeper the depth of the high temperature region,which is unfavorable for the prevention from spontaneous combustion of coal.Because of the existence of a tail roadway in the gob,the temperature field expands,the temperature rises,the main way to release CO is the tail roadway,and the amount of CO from the tail roadway is 10 times that from the outlet.Finally,the results are verified by the field test data,indicating that the simulation results are correct and reliable. KEY WORDS gobs;heterogeneous porosity:coal combustion:oxidation heating:four-dimensional:dynamic simulation 收稿日期：201603-29 基金项目：国家自然科学基金资助项目(51474017)：新疆维吾尔族自治区基金资助项目(2014211013)工程科学学报，第 38 卷，第 10 期: 1350--1358，2016 年 10 月 Chinese Journal of Engineering，Vol． 38，No． 10: 1350--1358，October 2016 DOI: 10． 13374 /j． issn2095--9389． 2016． 10． 002; http: / /journals． ustb． edu． cn 非均质孔隙率采空区氧化升温规律四维动态模拟 周佩玲1) ，张英华1) ，黄志安1) ，袁 飞2，3) ，高玉坤1) ，王 辉1) ，孙 倩1) 1) 北京科技大学金属矿山高效开采与安全教育部重点实验室，北京 100083 2) 北京科技大学冶金与生态工程学院，北京 100083 3) 北京科技大学钢铁冶金新技术国家重点实验室，北京 100083  通信作者，E-mail: huang_za@ qq． com 摘 要 运用 Fluent 动网格模型实现采空区的四维动态变化，并用用户自定义函数将煤低温氧化动力学机理及非均质孔隙 率函数编入 Fluent 中，结合时间和空间，对 U + L 型通风系统采空区升温规律进行四维动态模拟研究． 研究表明: 非均质孔隙 率四维动态模型能更真实地反应孔隙率的空间与时间变化，空间某一位置的孔隙率随时间呈负指数递减; 工作面推进速度越 大，采空区升温速率越小，推进速度为 3. 6 m·d － 1时平均升温速率仅为推进速度为 1. 2 m·d － 1时的 1 /5; 然而，推进速度越大，高 温点的深度越大，不利于自燃的预防; 尾巷的存在使得温度场范围扩大，温度升高，CO 主要从尾巷流出，尾巷释放的 CO 量是 回风巷 CO 释放量的 10 倍． 最后利用现场实测的数据对结果进行验证，表明模拟结果是正确可信的． 关键词 采空区; 非均质孔隙率; 煤燃烧; 氧化升温; 四维; 动态模拟 分类号 TD75 + 2 收稿日期: 2016--03--29 基金项目: 国家自然科学基金资助项目( 51474017) ; 新疆维吾尔族自治区基金资助项目( 2014211B013) 4D dynamic simulation of coal oxidation heating law in gobs with heterogeneous porosity ZHOU Pei-ling1) ，ZHANG Ying-hua1) ，HUANG Zhi-an1)  ，YUAN Fei2，3) ，GAO Yu-kun1) ，WANG Hui1) ，SUN Qian1) 1) Key Laboratory of High-Efficient Mining and Safety of Metal Mines ( the Ministry of Education) ，University of Science and Technology Beijing， Beijing 100083，China 2) School of Metallurgical and Ecological Engineering，University of Science and Technology Beijing，Beijing 100083，China 3) State Key Laboratory of Advanced Metallurgy，University of Science and Technology Beijing，Beijing 100083，China  Corresponding author，E-mail: huang_za@ qq． com ABSTＲACT Using the Fluent dynamic mesh model to realize the 4D dynamic movement of a gob and inputting the kinetic mecha￾nism of coal low-temperature oxidation and the dynamic change of heterogeneous porosity into Fluent through the user defined function， 4D dynamic simulation is performed on the spontaneous heating law in a gob for U + L ventilation of a certain mine． The results show that unsteady heterogeneous porous media can be more realistic in response to the change of porosity，which exponentially decreases with time． The greater the advancing speed of the working face，the smaller the heating rate，and the average heating rate at an advan￾cing speed of 3. 6 m·d － 1 is 1 /5 of that at an advancing speed of 1. 2 m·d － 1 ． However，the higher the advancing speed，the deeper the depth of the high temperature region，which is unfavorable for the prevention from spontaneous combustion of coal． Because of the existence of a tail roadway in the gob，the temperature field expands，the temperature rises，the main way to release CO is the tail roadway，and the amount of CO from the tail roadway is 10 times that from the outlet． Finally，the results are verified by the field test data，indicating that the simulation results are correct and reliable． KEY WOＲDS gobs; heterogeneous porosity; coal combustion; oxidation heating; four-dimensional; dynamic simulation