瓦斯抽放条件下采空区流场数值模拟

煤矿瓦斯抽放掘进法的数值模拟研究煤矿瓦斯爆炸事故一直是我国煤矿安全的“绊脚石”，为了避免这类事故的发生，煤矿瓦斯抽放掘进法成为了煤矿安全的重要手段之一。

该方法通过采取抽放设备，将瓦斯通道内的瓦斯抽出进行处理，起到了防止煤矿瓦斯爆炸的作用。

近年来，数值模拟技术在煤矿瓦斯抽放掘进法中的应用越来越广泛，本文将探讨煤矿瓦斯抽放掘进法的数值模拟研究。

一、数值模拟的基础理论数值模拟是指将实际问题通过数学建模、计算机模拟的方法，得到与实际问题尽可能相近的数值解。

其中，数值模拟的基础理论包括：计算流体力学（CFD）、数值计算方法等。

CFD方法是指通过数值方法求解流体运动的物理方程，其基本原理是将流体分割为无数个小单元，在每个小单元内求解流体状态。

而数值计算方法包括基于有限元、有限差或者其他数值算法的方法，它们的基本原理是将物理问题离散化，通过数字计算得到问题的数值解。

这两个基础理论在数值模拟中是重要的支撑，也是煤矿瓦斯抽放掘进法的数值模拟研究的基础。

二、模拟方法综述煤矿瓦斯抽放掘进法的数值模拟方法主要包括多相流模拟、离散相模拟和耦合模拟等。

多相流模拟是指将流体和固体或多个流体进行数值模拟的方法，通过将瓦斯、颗粒和煤尘等不同相进行数值计算来分析瓦斯抽放的可行性。

而离散相模拟则是通过计算颗粒分布、运动轨迹等参数，来分析其对瓦斯抽放的影响，同时也可以描绘出煤尘分布的情况。

耦合模拟则是将多相流和离散相模拟相结合，综合考虑不同相对瓦斯抽放的影响，以达到更加精确的模拟效果。

三、模拟结果分析煤矿瓦斯抽放掘进法的数值模拟研究在其内部排风系统的设计、瓦斯抽放量、抽放速度等方面起到了重要的作用。

通过数值模拟，可以得到模拟结果，并分析出其对实际工程中的瓦斯抽放的影响。

同时，模拟结果还可以用于模型优化、方案比较、风险分析等方面。

瓦斯抽放掘进法的数值模拟研究还需要结合工程实际，对实际工程环境进行准确的模拟。

目前的数值模拟技术还存在一些问题，如计算精度不高、模型参数固定等方面的问题，需要进一步的改进和优化。

Y型通风采空区瓦斯流场数值模拟研究的开题报告
一、研究背景及意义
煤矿采矿过程中，由于地质构造、煤体构造、采矿工艺等原因，难免会产生瓦斯，如果不能及时有效地排放，就会产生矿井瓦斯事故。

因此，瓦斯治理对于安全生产至关重要。

而煤矿瓦斯治理的核心是通风系统，通风系统的设计和优化是保障矿井安全的必要条件。

Y型通风采空区是煤矿通风系统中常见的一种形式，研究其瓦斯流场规律，对于提高煤矿通风系统的运行效率和安全性具有重要意义。

二、研究目的
本文旨在通过数值模拟方法研究Y型通风采空区的瓦斯流场规律，探究其影响因素和优化策略，为煤矿通风系统的运行和瓦斯治理提供理论依据和技术支持。

三、研究内容
1.建立Y型通风采空区瓦斯流场数值模型；
2.探究Y型通风采空区瓦斯流场规律，分析其影响因素；
3.优化Y型通风采空区通风系统，提高通风效率和安全性；
4.分析Y型通风采空区瓦斯爆炸危险性，提出瓦斯治理建议。

四、研究方法
本文采用计算流体力学（CFD）数值模拟方法，利用Fluent等软件建立Y型通风采空区瓦斯流场模型，并进行数值计算分析。

根据数值模拟结果，结合现场实验数据和文献资料，分析影响因素和优化策略，提出瓦斯治理建议。

五、预期结果
通过数值模拟研究，预计可以得出以下结论：
1.Y型通风采空区的内部瓦斯流场规律及其影响因素；
2.Y型通风采空区通风系统的优化方案，提高通风效率和安全性；
3.Y型通风采空区中瓦斯爆炸的危险性分析及瓦斯治理建议。

六、研究意义
本研究可为煤矿通风系统的建设和瓦斯治理提供理论基础和实践支持，有助于提高煤矿采空区的安全性和运行效率，对于促进我国煤炭行业的健康可持续发展具有重要意义。

基于COMSOL的采空区瓦斯抽采数值模拟研究胡延伟;孙路路;江城浩;黄腾瑶;陈连军【摘要】为了研究采空区内瓦斯达到稳定后的分布规律,从而确定瓦斯抽采巷道的位置,结合孔庄煤矿7433工作面实例,基于“O型圈”理论,采用分块赋值孔隙率的方法,通过COMSOL有限元分析软件对采空区瓦斯分布规律进行了数值模拟.模拟结果表明:工作面漏风不断流入采空区与瓦斯持续解吸涌出形成了1个动态平衡结果;7433工作面回采至180 m处时瓦斯富集,可以确定瓦斯抽采巷处于裂隙带上,瓦斯抽采巷道内错距离在10～30 m范围时抽采效果达到最优.【期刊名称】《煤矿安全》【年(卷),期】2018(049)005【总页数】5页(P167-170,174)【关键词】采空区;瓦斯分布;瓦斯抽采;瓦斯涌出;数值模拟【作者】胡延伟;孙路路;江城浩;黄腾瑶;陈连军【作者单位】山东科技大学矿业与安全工程学院,山东青岛266590;山东科技大学矿业与安全工程学院,山东青岛266590;山东科技大学矿山灾害预防控制-省部共建国家重点实验室培育基地,山东青岛266590;山东科技大学矿业与安全工程学院,山东青岛266590;山东科技大学矿业与安全工程学院,山东青岛266590;山东科技大学矿业与安全工程学院,山东青岛266590;山东科技大学矿山灾害预防控制-省部共建国家重点实验室培育基地,山东青岛266590【正文语种】中文【中图分类】TD712随着我国煤矿开采逐渐向深部发展，地质环境受多种因素影响，瓦斯已经成为威胁煤矿安全生产及工人生命安全的最重要因素[1-3]。

瓦斯抽采是防治瓦斯灾害的主要技术措施之一，但由于相应的瓦斯抽采理论的匮乏以及研究手段的制约，导致采空区瓦斯抽采这一手段的利用相对滞后，尤其是对位于工作面后方的采空区，其具有人员不能进入且很难被常规方法处理的特殊性及困难性，应用模拟仿真的方法来解决这一问题已经成为当前研究领域一种主流研究方法[4-6]，除此以外采空区内的瓦斯涌出运移过程极其其复杂，往往涉及到高阶非线性的偏微分方程，针对以上所提主要问题，采用多物理场仿真模拟软件COMSOL Multiphysics来模拟[7]采空区内瓦斯动态平衡后的分布规律，为确定采空区瓦斯治理提供依据使采空区瓦斯分布规律特征的理论讨论更趋向于完善。

煤矿采空区瓦斯渗流规律及其数值模拟研究煤矿采空区上覆岩层结构和移动规律分析综放工艺在开采高含量瓦斯厚煤层的推广应用中之所以遇到困难，往往是由于综放面上隅角瓦斯易超限，从而被迫断电撤人、中断生产所导致的。

上隅角瓦斯的主要来源一是工作面煤壁释放出的瓦斯，二是采煤工作面新采落下来的煤炭中散发出来的瓦斯，三是从采空区涌出的瓦斯，其中采空区涌出瓦斯是主要的来源。

由于采动影响在采动断裂带形成的破断裂隙和离层裂隙，采动裂隙网络与采空区相连通形成采动断裂带，由于瓦斯的升浮、扩散和渗透作用，在采动断裂带形成瓦斯富集区，这是瓦斯抽采的重点区域。

因此，要研究采空区内瓦斯的渗流规律，有必要先研究采空区岩体的垮落特征，按照采场覆岩横向采动特征，将采空区按照自然堆积区、载荷影响区和压实稳定区在横向进行划分，弄清各区碎胀系数、空隙率的分布特点；研究采空区上覆岩层采动断裂带的高度、碎胀系数及空隙率等特征，以便较全面地分析和研究采空区内空气—瓦斯混合气体在冒落带和采动断裂带内的渗流规律。

采空区瓦斯流场数学模型研究煤矿采空区内的瓦斯流动情况，建立起瓦斯流场的数学模型，对于认识采空区内瓦斯的真实流动状况以及对于进行数值模拟都有重要的基础意义。

垮落带之上的采动断裂带，在存在破断裂隙和离层裂隙相互贯通的同时，煤岩体内的裂隙还会与综放采场和采空区连通。

研究瓦斯在采动断裂带内的渗流、升浮和扩散原理，可以为解释采动断裂带是瓦斯聚集带，为其内布置钻孔抽采、巷道排放等瓦斯治理技术提供科学依据。

求解方法的选择FLUENT提供三种不同的解格式：分离解；隐式耦合解；显式耦合解。

三种解法都可以在很大流动范围内提供准确的结果，但是它们也各有优缺点。

分离解和耦合解方法的区别在于，连续性方程、动量方程、能量方程以及组分方程的解的步骤不同，分离解是按顺序解，耦合解是同时解。

两种解法都是最后解附加的标量方程（比如：湍流或辐射）。

隐式解法和显式解法的区别在于线化耦合方程的方式不同。

论文题目：采空区瓦斯运移规律数值模拟专业：应用数学硕士生：李书兵（签名）_____________ 指导教师：曹根牛（签名）_____________摘要长期的生产实践证明，矿井瓦斯运移规律预测是瓦斯防治不可缺少的重要技术环节。

矿井瓦斯是煤矿生产过程中主要的不安全因素，瓦斯灾害是威胁我国煤矿安全生产的最主要灾害，所以加强对瓦斯灾害的治理是保障矿井安全、高效生产的必要前提，瓦斯的运移、分布规律对矿井的设计、建设和开采都有重要影响。

随着开采深度和产量的增加，瓦斯潜在的影响更加显著，其模拟结果的正确与否，将直接影响矿井的技术经济指标。

为此，提出了研究矿井“采空区瓦斯运移规律的数值模拟”研究课题。

本文通过对陈家山矿416工作面采空区的长期实验观测，得到了采空区瓦斯分布的数据，首先分析了该采空区瓦斯主要来源于邻近层涌入的瓦斯和本煤层开采涌出的瓦斯，并对影响采空区瓦斯运移的因素进行层次分析法分析，得出漏入采空区的氧气是影响采空区中瓦斯运移最重要的因素，其次运用数学中的迭代方法对采空区中的非线性渗流方程进行改进，划分出该采空区中不同的流态区域，最后在总结目前采空区瓦斯运移特点、运移规律预测及数值模拟的研究的基础上，建立采空区与巷道风流流动场方程的数学模型，运用FLUENT软件进行数值模拟，以图形的方式展示了采空区瓦斯的分布规律，为采空区瓦斯治理及工作面瓦斯治理提供了可靠的理论依据。

关键词：采空区；瓦斯；层次分析法；迭代；数值模拟研究类型：应用研究Subject ：Numerical Simulation of Gas Migration Rule in Goaf Specialty ：Applied MathematicsName ：Li Shubing (Signature) ______________ Instructor ：Cao Genniu (Signature) ______________ABSTRACTThe long-term production proved that the prediction of migration rule of gas in mine prediction is an important and indispensable technology for mine's prevention and control link. mine gas is the main factors of insecurity in the process of coal production, Gas disaster is the major disaster of threat to China's mine safety production.so it is a necessary prerequisite to strengthen the governance of gas disaster for protecting the mine safety and efficient production and it is a significant impact to gas transport and distribution for mine design, construction and mining. With the increase in mining depth and production,gas potential impact will have more significant and the simulation results of the right or wrong, it will directly affect the technical and economic indicators of mine.Toward this end,we made a study of “Numerical Simulation of Gas Migration Rule in Goaf”.This article have goaf gas distribution data through the long-term experimental observations in the Chenjiashan Mine Goaf 416 face.First of all we analysis the gob gas mainly from the adjacent layer of the influx of gas and the coal seam gas emission and use AHP to analysis the factors affecting goaf gas migration，and found oxygen gas in goaf is the most important factor to affect the gas transport Second, we use the iterative method in mathematics to improve the nonlinear flow equations in goaf and carved out the different flow pattern of regional.Finally,we establish Merry mathematical model of flow field equations in the goaf and the roadway based on the study of summarizing the current characteristics of goaf gas migration, migration laws of prediction and numerical simulation. This article use FLUENT software for numerical simulation in order to display goaf gas flow distribution in graphical. Provides a reliable theoretical basis for gas governance and management.Keywords:Goaf Gas Analytic Hierarchy Process Iteration Numerical Simulation Thesis : Application Research1 绪论1.1 选题背景煤炭是我国一次能源的主体，煤炭工业承载着经济发展、社会进步和民族振兴的历史重任，是现代工业的血液，同时煤矿安全工作是全国工业安全工作的重中之重。

抽采条件下综放采空区漏风流场数值模拟
苏来旺
【期刊名称】《中国科技信息》
【年(卷),期】2014(000)009
【摘要】通过现场实测抽采条件下综放采空区瓦斯浓度、氧气浓度、工作面压力、采空区遗煤耗氧速度、采空区遗煤瓦斯释放速度、瓦斯抽采量.假定综放采空区为
各项同性的多孔介质,建立综放采空区漏风流场数学模型.本文使用了ANSYS软件
和FLUENT软件,其中ANSYES软件是先用DesignModeler进行三维建模(针对采空区),之后对综放采空区有关参数进行说明,这里运用的是UDF用户自定义函数.采用FLUENT软件(ANSYS环境下)对抽采条件下综放采空区漏风流场进行数值模拟,得出了抽采条件下综放采空区漏风流场规律.
【总页数】4页(P105-108)
【作者】苏来旺
【作者单位】安徽理工大学能源与安全学院,安徽淮南232001;山西汾西矿业集团
有限责任公司,山西介休032000
【正文语种】中文
【相关文献】
1.黄岩汇矿综放面采空区三维抽采条件下氧浓度数值模拟研究 [J], 朱兴攀;文虎;王凯;乔懿麟
2.三源两汇综放面采空区漏风与“三带”分布模拟研究 [J], 王康明;周生国;杨丁丁
3.综放采空区漏风流场数值模拟研究 [J], 刘伟华
4.新峪矿6114综放面回采期间采空区漏风测定 [J], 李斌;张景钢;刘海忠;康怀宇
5.大倾角坚硬顶板综放面采空区漏风数值模拟 [J], 谢振华
因版权原因，仅展示原文概要，查看原文内容请购买。

高瓦斯孤岛工作面采空区流场分布规律数值模拟及应用谷守生【摘要】为了研究高瓦斯煤层孤岛工作面采空区的流场分布规律，以天池矿15102孤岛工作面为对象，采用FLUENT数值模拟软件，对15102工作面采空区的风流场、氧浓度场、自燃三带以及瓦斯浓度场分布进行了模拟，并根据模拟结果进行了现场应用。

结果表明：与孤岛工作面进风侧相邻的采空区局部存在自燃危险性，通过采用FR-1阻化泡沫灌注采空区，可有效预防15102孤岛工作面采空区浮煤自然发火，保证工作面的安全回采。

%15102 island coal seam is taken as the object in order to study flow field distribution law of face goaf at high gas island. By use of FLUENT numerical simulation software, wind flow field, oxygen concentration field, spontaneous combustion three zones and gas concentra-tion distribution of the goaf of 15102 face are simulated, and the simulation results are put into applications. The results show that the goaf ad-jacent to the air inlet side is of spontaneous combustion risks. Use of FR-1 resistance foam to fill the goaf can effectively prevent 15102 island face from spontaneous combustion, so as to ensure the safe recovery at the mining face.【期刊名称】《煤矿现代化》【年(卷),期】2014(000)005【总页数】3页(P88-90)【关键词】孤岛工作面;自燃;采空区漏风;氧气浓度;阻化泡沫【作者】谷守生【作者单位】山西和顺天池能源有限责任公司，山西晋中 032700【正文语种】中文【中图分类】TD75+2.2通过大量研究表明:自然发火预测预报和防治(如合理配风量、注氮参数的确定)的基础工作之一就是划分采空区自燃"三带"[1-2].目前,国内外针对单一煤层采空区自燃"三带"划分已经较为完善,但对于高瓦斯煤层开采后形成的孤岛工作面采空区的流场分布以及自燃"三带"的合理划分还需要进一步的研究,这对于保证孤岛工作面安全、高效回采具有重要意义,同时对于探索高瓦斯矿井综放面采空区自燃"三带"分布特征也有重要的理论价值[3-6].采煤工作面采空区的自燃"三带"划分目前通常可以使用理论计算(数值模拟)以及现场仪器监测等方法来确定,现场仪器监测方法通常情况下周期较长,实施工艺也较为复杂,采空区冒落情况也会影响测试结果;然而数值模拟方法计算方便,输入参数后,便可模拟计算不同情况下的结果,能够预先或及时指导工作面的实际开采情况[7-8].本文通过数值模拟方法,以天池煤矿15102工作面为模型,来研究孤岛工作面采空区的流场分布规律以及自燃"三带"分布特征,并根据模拟进行现场应用.研究表明:煤层开采后,采空区内的岩体和松散的煤体均可被视为多孔均匀介质,因此,在整个采空区的计算域内,可以假设流体的密度不变,采空区中的气体流动符合达西定律.采空区中的氧气与浮煤发生复杂的化学反应产生CO2和CO等气体,同时消耗一部分浮煤.另外由于产生与消耗的气体一样,空气质量没有变化,但由于空气中存在浓度差,从而使不同组分的气体按照Fick定律扩散.煤自燃氧化过程通常来说比较缓慢,于是在不考虑能量的条件下,可以假设采空区的气体渗流、扩散与化学反应是一个稳态的过程,符合如下公式[9-13]:式中:x、y、z分别是不同方向的轴坐标,Q是工作面的漏风强度(即单位面积煤岩通过的漏风量),Qx、Qy、Qz分别是漏风强度在x、y、z三个方向上的分量,H是气体的压力(Pa),是绝对渗透率,μ是空气粘性系数,常温下取μ=1.7894X10-5kg/(m.s),是煤的耗氧速度.2.115102 工作面概况天池煤矿15102综放工作面位于15101和15103工作面中间,15101及15103综放工作面先于该工作面开采,15102综放工作面开采时就成为了孤岛工作面.该工作面面长160m,推进方向900m,煤厚4.6~6.2m,平均5.3m.采用综合机械化放顶煤开采技术,采高2.6m,采放比为1:1.04.瓦斯相对涌出量是20.95m3/t,瓦斯绝对涌出量是66.6m3/min,煤层为自燃煤层.2.2 数值模型的建立正确的数值模拟结果需要与现场相同条件的参数设置,因此,通过15102工作面相关参数的测试(结果见表1),结合15102工作面配风量、大气参数、瓦斯涌出情况以及顶板冒落情况等条件,构建采空区的数学模型,并对模型各个参数进行赋值,从而模拟15101、15102和15103工作面的采空区流场分布特征.根据天池煤矿15102采煤工作面的实际情况,同时为了便于模型的建立和合理参数的选取,在尽量与15102工作面实际情况相符合的原则下,数值模型利用FLUENT 里的GAMBIT建模软件来建立.模型分为采空区多孔介质空间渗流部分和工作面自由空间紊流部分.依据15102工作面的实际供风量和模型解算需要,速度入口边界是进风巷,风量2800m3/min,速度3.46m/s;自由出口边界是回风巷及尾巷,其出口风量分别是1500m3/min和1150m3/min,速度分别是2.65m/s和2.03m/s.通过数值模拟可以得到如图1~4所示的采空区气体风速矢量及等值线、氧浓度等值线、自燃三带和瓦斯浓度等值线等分布.由图1可知,由于采动影响,15101采空区、15103采空区与15102工作面之间的煤柱发生破坏存在漏风通道,导致15103和15101采空区有少量来自于15102工作面的漏风流,漏风流速度最大0.00003m/s.从图2可知,由于15102进风巷风压较大,加之采动导致煤柱之间的漏风通道,15102工作面进风巷附近的15103工作面采空区的局部地段氧气浓度超过8%,此处浮煤自燃机率较大;15102工作面回风巷附近的15101工作面采空区氧气浓度较低,最大为5%,此处浮煤不太可能自燃.靠近15102采煤工作面处在瓦斯尾巷和回风巷之间有两道联络巷已打开,上隅角的大部分瓦斯通过这两道联络巷经瓦斯尾巷排出,上隅角瓦斯浓度低于1%.研究表明:一般可以根据氧气浓度等值线(8% ~18%为自燃带),同时结合风速等值线图(0.004~ 0.001%为自燃带)来确定采空区自燃"三带".则由图3可知,回风巷一侧采空区的自燃带与工作面间距是32m~49m,进风巷一侧是27m~49m.15102采煤工作面是一个孤岛面,因此其自燃三带有其较为特殊的方面,15103采空区临近进风巷4m~5m区域的氧气浓度均高于8%,说明该处采空区遗煤有可能发生自燃,而15101采空区内的氧气浓度最大仅为5%,说明该采空区不存在自燃的可能.从4图可知,瓦斯浓度在进风巷和回风巷较近处变化幅度较大,而中间一般较小.瓦斯浓度在5~15%范围的位置,在进风巷一侧与工作面距离为15m~30m回风巷一侧与工作面间距在30m之后.综合图1~4可得到15102工作面采空区遗煤自燃与瓦斯爆炸界限范围如表2所示.针对15102工作面的现场条件,结合数值模拟结果,选用FR-1阻化泡沫来进行堵漏风,同时作为阻化剂防止采空区的遗煤发生自燃,如图5所示.分别从15102采煤工作面两侧巷道,每隔30m向邻近的15103与15101采空区施工φ42mm的钻孔,同时在工作面向采空区施工与进风巷的间距是40m和80m的φ42mm的两组钻孔.所有钻孔均要进入采空区,每个钻孔日注泡沫不小于2h.15102工作面在2009年12月底到2010年12月共推进900m,平均日推进度2.5m,月推进度75m.在该工作面开采期间,通过应用采空区气体监测、控制风量、采空区两顺槽砌筑隔离档墙、注液体阻化泡沫等综合预防措施,15102回风巷的瓦斯浓度在0.11%~0.8%之间,平均为0.48%,回风巷、尾巷的CO浓度在0~0.0014%,说明采用FR-1阻化泡沫灌注采空区的防灭火技术能有效预防15102孤岛工作面采空区浮煤自燃.由于现场测定工艺复杂、周期较长,通过数值模拟可以较为准确的模拟孤岛工作面采空区的流场分布规律,为工作面防灭火工作和瓦斯治理提供指导.(1)15102孤岛工作面的相邻两工作面都会出现漏风现象,但漏风量不大.15103采空区内靠近15102进风顺槽的部分氧气浓度在8%以上,这部分采空区浮煤有出现自燃的可能,15101采空区遗煤没有自燃危险.(2)工作面进风巷一侧采空区的自燃带与工作面间距是27m~49m,回风巷一侧是32m~49m.瓦斯浓度在5%~15%范围的位置,在进风巷一侧与工作面距离为15m~30m回风巷一侧与工作面间距在30m之后.(3)采用FR-1阻化泡沫灌注采空区的防灭火技术能有效预防15102孤岛工作面采空区浮煤自然发火和可能引发的瓦斯爆炸,从而保证了15102采煤工作面的安全开采.【相关文献】[1]杜志刚.大规模采空区煤炭自燃与防治关键技术研究[J].中国安全生产科学技术,2012,8(增):15-18.[2]李宗翔,吴强,肖亚宁.采空区瓦斯涌出与自燃耦合关系基础研究[J].中国矿业大学学报,2008,37(1):38-42.[3]李化敏,王文,熊祖强.采动围岩活动与工作面瓦斯涌出关系[J].采矿与安全工程学报,2008,25(1):11-16.[4]曹代勇,樊新杰,吴查查,等.内蒙古乌达煤田火区相关裂隙研究[J].煤炭学报,2009,34(8):1009-1014.[5]赵聪,陈长华.基于模糊渗流理论的采场自然发火[J].辽宁技术工程大学学报(自然科学版),2009,28(S):31-33.[6]褚廷湘,余明高,杨胜强,等.瓦斯抽采对U+Ⅱ型近距离煤层自燃的耦合关系[J].煤炭学报,2010,35(12):2082-2087.[7]徐精彩.煤自燃危险区域判定理论[M],北京:煤炭工业出版社,2001.[8]崔凯,张东海,杨胜强.采空区遗煤自燃带确定及风流场数值模拟[J].山东科技大学学报(自然科学版)2002,21(4): 88-92.[9]文虎.综放工作面采空区煤自燃过程的动态数值模拟[J].煤炭学报,2002,27(1):54-58.[10]张辛亥,刘灿,周金生,等.综放面采空区流场模拟及自燃危险区域划分[J].西安科技大学,2006(26)1:6-9.[11]吕文陵,杨胜强,徐全,等.高瓦斯矿井孤岛综放采空区遗煤自燃综合防治技术[J].中国安全生产科学技术,2010,6 (5):60-64.[12]谢军,薛生.综放采空区空间自燃三带划分指标及方法研究[J].煤炭科学技术,2011,39(1):65-68.。