基于数值模拟探究电厂烟气注入采空区对遗煤自燃的影响

基于数值模拟探究电厂烟气注入采空区对遗煤自燃的影响一、绪论1.1研究背景与意义在全球能源结构中，煤炭作为重要的一次能源，在电力生产、工业制造等领域占据着举足轻重的地位。我国是煤炭生产和消费大国，煤炭在我国能源结构中一直占据主导地位。尽管近年来随着能源结构的调整，煤炭占比有所下降，但在未来较长一段时间内，煤炭仍将是我国的主要能源来源之一。然而，在煤炭开采和利用过程中，采空区遗煤自燃以及电厂碳排放等问题，给能源利用、安全生产和环境保护带来了严峻挑战。采空区遗煤自燃是矿井安全生产的重大隐患，严重威胁着工作人员的生命安全。据相关统计，我国每年因采空区遗煤自燃导致的煤炭资源损失高达数千万吨，同时还会诱发瓦斯、煤尘燃烧或爆炸等连锁事故，造成巨大的经济损失和环境污染。煤炭自燃的发生需要同时具备煤具有自燃倾向性、有连续不断的供氧条件以及有热量易于积聚的客观环境这三个条件。其过程是一个极其复杂的物理化学变化过程，除了取决于煤炭本身内在的物理、化学、力学性质外，还与地质条件、开拓条件、通风条件等因素密切相关。目前，普遍采用以注氮为主的综合防灭火措施来预防采空区煤炭自燃。然而，制氮过程能耗高、成本大，造成了能源的浪费。近年来，随着我国煤电一体化建设的推进，大量坑口电厂应运而生。电厂烟道气体作为燃烧产物，主要成分是CO₂等惰性气体，含量达80%以上，氧含量在10%以下。若将电厂烟气通过管路直接注入井下采空区，不仅能有效预防和控制煤炭自燃火灾，还可避免制氮造成的能源浪费，为企业节省大量资金。与此同时，全球气候变化问题日益严峻，减少温室气体排放已成为国际社会的共识。CO₂作为主要的温室气体，其减排和封存备受关注。将电厂烟气注入井下采空区，能够实现CO₂等气体的封存，有助于缓解温室效应，对环境保护具有重要意义。我国高度重视节能减排工作，国务院印发的《节能减排综合性工作方案》明确提出，在重点行业推广重大节能减排技术，并出台了一系列有利于节能减排的财税政策。利用电厂烟气预防遗煤自燃和实现碳封存，符合国家节能减排的战略要求，具有广阔的应用前景和重大的现实意义。通过深入研究电厂烟气预防遗煤自燃的机理及其碳封存效果，不仅可以为矿井防灭火工作提供新的技术手段和理论依据，有效降低采空区遗煤自燃的风险，保障矿井安全生产，还能为电厂烟气的资源化利用和CO₂减排提供创新思路，实现煤炭资源的清洁高效利用，促进能源与环境的协调发展，具有显著的经济效益、社会效益和环境效益。1.2国内外研究现状在自然发火防治方面，国内外学者进行了大量研究。国外对于煤炭自燃的研究起步较早，在自燃机理、预测方法以及防治技术等方面积累了丰富的经验。例如，美国、澳大利亚等煤炭资源丰富的国家，通过对不同煤种的自燃特性进行深入研究，建立了较为完善的煤炭自燃预测模型。在防治技术上，注氮防灭火技术在国外得到了广泛应用，其技术成熟度较高，设备先进，能够实现对采空区氧气浓度的有效控制。国内在自然发火防治领域也取得了显著进展。研究人员通过实验研究和数值模拟等手段，对煤炭自燃的机理和防治技术进行了深入研究。有学者利用自主研发的实验装置，模拟采空区遗煤自燃过程，研究发现煤的自燃倾向性与煤的变质程度、煤岩成分、水分含量等因素密切相关。还有学者通过数值模拟，分析了采空区漏风规律、氧气分布和温度场变化，为自然发火防治提供了理论依据。在惰气防灭火方面，注氮防灭火技术是目前应用最广泛的惰气防灭火技术之一。国内外学者对注氮防灭火的机理、工艺和设备进行了大量研究。研究表明，注氮能够降低采空区内氧气浓度，抑制煤炭氧化反应，从而达到预防自燃的目的。此外，还有学者研究了其他惰气如二氧化碳、氮气与二氧化碳混合气体等在防灭火中的应用。在采空区遗煤自燃数值模拟方面，随着计算机技术的发展，数值模拟方法在采空区遗煤自燃研究中得到了广泛应用。国内外学者建立了多种数学模型来模拟采空区遗煤自燃过程，包括传热传质模型、化学反应动力学模型、多场耦合模型等。通过数值模拟，可以深入研究采空区遗煤自燃的影响因素，如漏风风速、氧气浓度、煤体特性等，为防灭火措施的制定提供科学依据。然而，当前研究仍存在一些不足。在电厂烟气防治遗煤自燃方面，对于不同煤质、不同地质条件下电厂烟气的最佳注入参数和注入方式的研究还不够深入，缺乏系统的理论和技术指导。在碳封存方面，虽然开展了多个示范项目，但碳封存的成本较高，技术的稳定性和安全性仍有待进一步提高，且对于碳封存后长期环境影响的研究相对较少。此外，将电厂烟气防治遗煤自燃与碳封存相结合的研究还处于起步阶段，两者协同作用的机理和效果还需要进一步深入探究。1.3研究内容与方法本文主要研究内容如下：电厂烟气防治遗煤自燃的理论研究：分析电厂烟气的成分和特性，研究电厂烟气中各成分对煤炭氧化自燃的影响机制，从化学反应动力学、传热传质等角度探讨电厂烟气防治遗煤自燃的原理。采空区遗煤自燃的数学模型构建：综合考虑采空区漏风、氧气扩散、煤炭氧化放热、热量传递等因素，建立采空区遗煤自燃的多场耦合数学模型，并对模型进行求解和验证。电厂烟气注入采空区的数值模拟研究：利用建立的数学模型，对电厂烟气注入采空区的过程进行数值模拟，分析不同注入参数（如注入量、注入位置、注入时间等）对采空区内气体浓度分布、温度场变化和煤炭自燃进程的影响，优化电厂烟气注入方案。现场实测与验证：选择合适的煤矿现场，进行电厂烟气注入采空区的现场试验，监测采空区内气体成分、温度等参数的变化，验证数值模拟结果的准确性，为实际工程应用提供数据支持。本文采用的研究方法主要包括：文献研究法：查阅国内外相关文献，了解煤炭自燃防治、惰气防灭火以及碳封存等领域的研究现状和发展趋势，为本文的研究提供理论基础和参考依据。数值模拟法：运用CFD（计算流体力学）软件和自编程序，对采空区遗煤自燃过程和电厂烟气注入采空区的效果进行数值模拟，通过模拟结果分析影响因素和作用规律，为方案优化提供依据。实验研究法：设计并开展实验室实验，模拟采空区遗煤自燃和电厂烟气注入过程，研究电厂烟气对煤炭自燃的抑制作用，获取相关实验数据，验证数值模拟结果的可靠性。现场实测法：在煤矿现场进行电厂烟气注入采空区的试验，安装监测设备，实时监测采空区内气体成分、温度等参数的变化，对数值模拟和实验研究结果进行实际验证，为工程应用提供实践经验。二、电厂烟气防灭火理论及技术2.1电厂烟气来源及成分分析在我国，随着煤电一体化建设的大力推进，坑口电厂成为了电厂烟气的主要来源。坑口电厂依托煤矿资源，将煤炭在锅炉内进行燃烧，产生高温高压的蒸汽，驱动汽轮机发电。在这一过程中，煤炭中的可燃成分与空气中的氧气发生剧烈的氧化反应，生成了大量的烟气。这些烟气从锅炉尾部排出，经过一系列的处理后，一部分被排放到大气中，而另一部分则可作为潜在的防灭火资源，用于井下采空区的遗煤自燃防治。电厂烟气是一种成分复杂的混合气体，其主要成分包括CO₂、N₂、O₂、SO₂、NOx以及少量的CO、H₂O等。其中，CO₂等惰性气体含量达80%以上，这是由于煤炭燃烧过程中，碳元素与氧气充分反应生成CO₂，而N₂则主要来自于空气中未参与反应的部分。以某典型坑口电厂为例，其烟气中CO₂含量约为16%，N₂含量约为75%，两者之和超过了80%。而O₂含量在10%以下，这是因为在燃烧过程中，大部分氧气被消耗，剩余的氧气主要是由于燃烧不完全或空气过量引入的。SO₂和NOx是电厂烟气中的主要污染物，它们的产生与煤炭中的硫、氮元素含量以及燃烧条件密切相关。当煤炭中的硫元素在高温下与氧气反应时，会生成SO₂；而燃料型NOx则是由煤中氮元素在燃烧过程中热分解并进一步氧化产生，热力型NOx则是在高温下空气中氮气与氧气反应生成。这些污染物不仅对环境造成危害，如导致酸雨、光化学烟雾等问题，还可能对井下设备和人员健康产生影响。因此，在将电厂烟气注入采空区之前，需要对其进行严格的净化处理，以降低污染物的含量。此外，电厂烟气中还含有少量的CO和H₂O。CO是由于煤炭不完全燃烧产生的，其含量与燃烧的充分程度有关；H₂O则主要来自于煤炭中的水分蒸发以及燃烧过程中氢气与氧气反应生成。虽然这些成分的含量相对较少，但在某些情况下，它们也可能对采空区的环境和遗煤自燃产生一定的影响。2.2电厂烟气防灭火原理电厂烟气之所以能够起到防灭火的作用，主要是基于以下几个方面的原理：降低氧浓度：氧气是煤炭氧化自燃的必要条件之一。电厂烟气中氧含量在10%以下，远低于空气中21%的氧浓度。当将电厂烟气注入采空区后，烟气会与采空区内的空气混合，从而稀释了采空区内的氧气浓度。随着氧气浓度的降低，煤炭氧化反应的速率会显著减缓。因为氧气作为氧化剂，其浓度的降低使得煤炭与氧气接触并发生反应的机会减少，从而抑制了煤炭的氧化过程，降低了煤炭自燃的可能性。抑制氧化反应：电厂烟气中的CO₂、N₂等惰性气体，化学性质稳定，不易与其他物质发生化学反应。它们在采空区内形成一种惰性环境，包裹在煤炭表面，阻碍了氧气与煤炭的接触。这就如同在煤炭与氧气之间建立了一道屏障，使得煤炭无法获得足够的氧气来维持氧化反应的进行。例如，CO₂气体的密度比空气大，会在采空区内逐渐下沉，聚集在煤炭周围，有效地隔离了氧气，从而抑制了煤炭的氧化反应。吸附作用：煤是一种具有复杂表面结构的非晶体物质，对CO₂气体具有较强的吸附能力。研究表明，煤对CO₂气体的吸附量比煤对O₂及N₂的吸附量大2-3倍左右。当电厂烟气注入采空区后，煤会优先吸附CO₂。这种吸附作用不仅减少了煤炭表面可与氧气接触的活性位点，而且吸附在煤炭表面的CO₂分子形成了一层保护膜，进一步阻止了氧气与煤炭的接触，从而有效地阻止了煤的氧化自燃。降温作用：电厂烟气在注入采空区的过程中，由于其温度相对较低（一般低于采空区内的温度），会与采空区内的高温空气和煤炭进行热量交换，带走一部分热量。这有助于降低采空区内的温度，使煤炭的温度难以升高到自燃的临界温度，从而抑制了煤炭自燃的发生。2.3电厂烟气注入工艺系统及技术指标注入工艺系统：电厂烟气注入采空区的工艺系统主要包括烟气收集、净化处理、输送和注入四个部分。在烟气收集环节，通过烟道将电厂锅炉排出的烟气收集起来。由于电厂烟气中含有SO₂、NOx、粉尘等污染物，若直接注入采空区，会对井下环境和设备造成损害，因此需要进行净化处理。常用的净化处理方法包括脱硫、脱硝和除尘。脱硫工艺如石灰石-石膏法，通过向烟气中加入石灰石浆液，使SO₂与石灰石反应生成亚硫酸钙，再氧化为石膏，从而去除烟气中的SO₂；脱硝工艺如选择性催化还原法（SCR），在催化剂的作用下，向烟气中喷入氨气，使NOx与氨气反应生成氮气和水，实现脱硝；除尘则可采用电除尘器、袋式除尘器或电袋复合除尘器等设备，去除烟气中的粉尘。净化后的烟气通过专门的输送管路输送到井下采空区。输送管路一般采用耐腐蚀、耐压的管道，以确保烟气能够安全、稳定地输送。在注入环节，根据采空区的具体情况，选择合适的注入方式，将烟气注入到采空区内。2.注入方式：电厂烟气注入采空区的注入方式主要有埋管注气和钻孔注气两种。埋管注气是在采煤工作面推进过程中，沿采空区回风侧预埋注气管路。随着工作面的推进，采空区逐渐形成，烟气通过预埋管路注入到采空区内。这种方式的优点是注气均匀，能够覆盖较大的采空区范围；缺点是管路容易受到采空区顶板垮落等因素的影响，导致管路损坏，影响注气效果。钻孔注气则是在采空区周围的巷道或地面上，向采空区内打钻孔，通过钻孔将烟气注入采空区。这种方式的灵活性较高，可根据采空区的具体情况选择钻孔位置和深度；但注气范围相对较小，需要合理布置钻孔才能达到较好的注气效果。3.技术指标：电厂烟气防灭火惰化指标是衡量电厂烟气注入效果的重要依据。一般来说，采空区内空气中氧浓度值应＜7%，当氧浓度在这个范围内时，能够有效预防和控制煤炭自燃。此外，还需要考虑烟气的注入量、注入压力等指标。烟气注入量应根据采空区的大小、漏风情况以及煤炭自燃的危险程度等因素进行确定。例如，对于一个较大的采空区且漏风严重的情况，需要较大的烟气注入量才能保证采空区内的氧气被充分稀释。注入压力则应保证烟气能够顺利地输送到采空区内，并在采空区内均匀分布，但压力也不能过高，以免对采空区的顶板和巷道造成破坏。2.4电厂烟气防灭火技术评价成本优势：与传统的制氮防灭火技术相比，电厂烟气防灭火技术具有显著的成本优势。制氮过程需要专门的制氮设备，如空气分离设备等，这些设备投资巨大，运行成本高，包括设备的能耗、维护费用等。而电厂烟气是发电过程中的副产品，将其用于防灭火，无需额外的气体制备成本，只需对烟气进行简单的净化处理和输送，大大降低了防灭火的成本。效果显著：从防灭火效果来看，电厂烟气能够有效地降低采空区内的氧气浓度，抑制煤炭的氧化自燃。通过合理的注入工艺和参数控制，可以使采空区内的氧气浓度保持在较低水平，从而达到预防煤炭自燃的目的。同时，烟气中的惰性气体和吸附作用也进一步增强了防灭火的效果。环保效益：电厂烟气中含有CO₂等温室气体，将其注入采空区实现了CO₂的封存，减少了温室气体的排放，有助于缓解全球气候变化。此外，通过对电厂烟气的净化处理，降低了SO₂、NOx等污染物的排放，对改善大气环境质量具有积极意义。存在的不足：然而，电厂烟气防灭火技术也存在一些不足之处。首先，电厂烟气的成分和流量会受到电厂发电负荷、煤炭质量等因素的影响，导致烟气的稳定性较差。这就需要建立完善的监测和调控系统，实时监测烟气的成分和流量，并根据实际情况进行调整。其次，在烟气输送过程中，可能会出现管路腐蚀、堵塞等问题，需要定期对管路进行维护和清理。此外，对于一些特殊地质条件的采空区，如采空区顶板破碎、漏风严重等，电厂烟气的注入效果可能会受到一定的影响，需要进一步研究和优化注入工艺。三、工作面自燃规律的现场测定3.1工作面概况及自然发火特点本次研究选取了某煤矿的1101综采工作面作为研究对象。该工作面位于井田一采区东翼，所采煤层为3号煤层，煤层平均厚度为4.5m，倾角在8°-12°之间，属于近水平煤层。煤层赋存稳定，结构简单，无夹矸。顶底板岩性主要为泥岩和砂岩，顶板具有一定的稳定性，但在采动影响下，容易出现垮落现象。该工作面采用综采放顶煤采煤工艺，全部垮落法管理顶板。工作面长度为200m，推进长度为1500m。采煤机割煤速度为4-6m/min，日推进度为8-10m。在采煤过程中，会产生大量的浮煤，这些浮煤堆积在采空区内，为煤炭自燃提供了物质基础。从自然发火的潜在风险来看，该工作面具有以下特点：一是煤层具有自燃倾向性，根据煤样的自燃倾向性鉴定结果，3号煤层为Ⅱ类自燃煤层，自燃倾向性较强；二是采空区遗煤量大，综采放顶煤工艺使得采空区内遗留大量煤炭，这些煤炭在适宜的条件下容易发生自燃；三是采空区漏风严重，由于工作面采用全部垮落法管理顶板，采空区与工作面之间存在漏风通道，风流中的氧气能够进入采空区，为煤炭氧化提供了供氧条件；四是工作面推进速度相对较慢，导致采空区内遗煤氧化时间较长，热量容易积聚，增加了煤炭自燃的风险。3.2工作面漏风量测定采空区漏风原理及危害：采空区漏风是指在通风压力的作用下，风流从通风系统中通过采空区的裂隙、孔洞等通道，直接渗透到回风道或地面的现象。其形成的根本原因是采空区两端存在压差，以及采空区内部存在漏风通道。在该工作面，采空区与工作面之间的漏风通道主要包括采空区顶板垮落形成的裂隙、支架间的缝隙以及工作面上下隅角等。采空区漏风会带来诸多危害。风流中携带的氧气能够促进遗煤的氧化自燃，引发火灾。当漏风强度较小时，煤炭氧化产生的热量能够积蓄，且有连续的氧气供给，采空区遗煤容易自燃；而漏风强度较大时，虽然能带走部分热量，但也会增加煤炭与氧气的接触机会，在一定程度上促进煤炭氧化。漏风还可能导致瓦斯积聚，引发瓦斯灾害。因为漏风会改变采空区内的气体流动状态，使瓦斯难以排出，从而在局部区域积聚，当瓦斯浓度达到爆炸界限时，遇到火源就会发生爆炸。此外，采空区漏风会使有效风流减少，影响工作面的正常通风，威胁井下的安全生产及工人的生命安全。2.测定方法：为了准确测定工作面的漏风量，采用了示踪气体法。该方法的原理是在进风流中释放一定量的示踪气体（如六氟化硫SF₆），示踪气体随风流进入采空区，由于采空区漏风，示踪气体在回风流中的浓度会发生变化，通过测定进、回风流中示踪气体的浓度以及进风量，就可以计算出采空区的漏风量。具体操作步骤如下：首先，在工作面进风巷距离工作面50m处设置示踪气体释放点，利用气体释放装置以恒定的速度释放SF₆气体，释放时间为30min，释放量为1000mL。然后，在工作面回风巷距离工作面50m处以及采空区回风隅角设置气体采样点，使用气相色谱仪每隔5min采集一次气样，分析气样中SF₆气体的浓度。同时，使用风表在进风巷和回风巷的测风站测量进风量和回风量，测量时间与气样采集时间同步，每次测量持续10min，取平均值作为测量结果。3.漏风数据及分析：经过连续3天的测量，得到的漏风数据如表1所示：测量日期进风量(m³/min)回风巷SF₆浓度(ppm)采空区回风隅角SF₆浓度(ppm)计算漏风量(m³/min)第一天25000.51.2300第二天24500.451.1280第三天25500.551.3320从表中数据可以看出，该工作面的漏风量较大，平均漏风量达到300m³/min左右。这主要是由于采空区顶板垮落形成的裂隙较多，以及工作面上下隅角封堵不严，导致漏风通道较多。漏风量的大小与进风量、采空区的通风阻力以及漏风通道的分布等因素密切相关。进风量越大，在相同的通风阻力下，漏风量也会相应增大。而采空区通风阻力的大小则取决于顶板垮落情况、支架布置以及漏风通道的复杂程度等。3.3工作面自燃“三带”测定束管监控系统布置：为了准确划分工作面自燃“三带”，在该工作面布置了束管监控系统。束管监控系统主要由束管、采样泵、气体分析仪、数据采集与处理系统等组成。束管选用16芯聚乙烯管，具有良好的柔韧性和耐腐蚀性。在工作面进风巷、回风巷以及采空区内共布置了10个测点，其中进风巷2个测点，分别位于距离工作面50m和100m处；回风巷3个测点，分别位于距离工作面50m、100m和150m处；采空区内5个测点，沿着采空区走向方向，每隔50m布置一个测点，距离工作面最近的测点距离为50m。每个测点都安装了粉尘过滤器和滤水器，以防止粉尘和水分进入束管，影响测量结果。束管从测点引出后，沿着巷道壁敷设，通过分路箱连接到采样泵。采样泵将采空区内的气体抽到地面的气体分析仪中进行分析。气体分析仪采用气相色谱仪，能够准确测量气体中的氧气、一氧化碳、二氧化碳、甲烷等成分的浓度。数据采集与处理系统实时采集气体分析仪的测量数据，并进行存储、分析和显示。2.自燃“三带”划分及分析：根据实测的氧浓度数据，采用氧气浓度划分法对采空区自燃“三带”进行划分。目前，国内外广泛认可的划分标准是将氧气浓度的8%和18%作为采空区氧浓度的极限值，当氧气浓度超过18%时，为散热带；8%-18%为氧化自燃带；低于8%则为窒息带。通过对各测点氧浓度数据的分析，得到该工作面采空区自燃“三带”的分布情况如图1所示：[此处插入采空区自燃“三带”分布示意图]从图中可以看出，散热带主要分布在距离工作面较近的区域，宽度约为50-80m。在这个区域内，由于漏风较大，风流能够及时带走煤炭氧化产生的热量，使得煤炭温度难以升高，氧化反应速度较慢。氧化自燃带位于散热带的后方，宽度约为80-150m。在这个区域内，氧气浓度适中，漏风强度相对较小，煤炭氧化产生的热量能够积聚，具备煤炭自燃的条件，是煤炭自燃的危险区域。窒息带位于氧化自燃带的后方，氧气浓度低于8%，由于缺乏氧气，煤炭氧化反应受到抑制，一般不会发生自燃。该工作面自燃“三带”的分布受到多种因素的影响，如采空区漏风、顶板垮落情况、煤炭自燃倾向性等。采空区漏风是影响自燃“三带”分布的关键因素，漏风强度和方向决定了氧气的供给和热量的传递，从而影响煤炭氧化自燃的进程。顶板垮落情况会影响采空区的通风阻力和漏风通道的分布，进而影响自燃“三带”的范围和位置。而煤炭自燃倾向性则决定了煤炭氧化自燃的难易程度，自燃倾向性越强的煤炭，在相同的条件下越容易发生自燃。四、采空区遗煤自燃三维多场耦合模型构建4.1三维多场耦合解算模型流场控制方程：采空区内的风流运动满足质量守恒定律和动量守恒定律，其连续性方程和动量方程是描述流场的基础。连续性方程表示在单位时间内，流入和流出控制体的流体质量之差等于控制体内流体质量的变化率，对于不可压缩流体，其数学表达式为\nabla\cdot\vec{u}=0，其中\vec{u}为流体速度矢量。动量方程则描述了流体动量的变化率等于作用在流体上的外力之和，包括压力梯度力、粘性力和重力等，其表达式为\rho(\frac{\partial\vec{u}}{\partialt}+\vec{u}\cdot\nabla\vec{u})=-\nablap+\nabla\cdot(\mu\nabla\vec{u})+\rho\vec{g}，其中\rho为流体密度，p为压力，\mu为动力粘度，\vec{g}为重力加速度。多组分扩散方程：氧气、一氧化碳等气体在采空区内的扩散过程对遗煤自燃起着关键作用。根据费克定律，多组分扩散方程可表示为\frac{\partial(\rhoY_i)}{\partialt}+\nabla\cdot(\rho\vec{u}Y_i)=\nabla\cdot(\rhoD_{i,m}\nablaY_i)+S_i，其中Y_i为第i种组分的质量分数，D_{i,m}为第i种组分的扩散系数，S_i为第i种组分的源项，主要包括煤炭氧化反应产生或消耗的气体量。温度场控制方程：温度场的变化是遗煤自燃过程的重要特征，其控制方程基于能量守恒定律。在采空区内，能量的传递包括热传导、对流和煤炭氧化放热等过程，温度场控制方程可表示为\rhoc_p(\frac{\partialT}{\partialt}+\vec{u}\cdot\nablaT)=\nabla\cdot(k\nablaT)+Q_{coal}，其中c_p为比热容，T为温度，k为热导率，Q_{coal}为煤炭氧化反应产生的热量。煤炭氧化反应是一个复杂的化学反应过程，会释放出大量的热量，其产热速率与煤的氧化反应速率、反应热等因素有关。湍流模型：采空区内的风流通常处于湍流状态，为了准确描述湍流对各物理场的影响，选择合适的湍流模型至关重要。标准k-\epsilon湍流模型是工程中常用的湍流模型之一，它通过求解湍动能k和湍流耗散率\epsilon的输运方程来封闭雷诺应力项。湍动能k的方程为\frac{\partial(\rhok)}{\partialt}+\nabla\cdot(\rho\vec{u}k)=\nabla\cdot[(\mu+\frac{\mu_t}{\sigma_k})\nablak]+G_k-\rho\epsilon，其中G_k为湍动能的生成项，主要由平均速度梯度产生；湍流耗散率\epsilon的方程为\frac{\partial(\rho\epsilon)}{\partialt}+\nabla\cdot(\rho\vec{u}\epsilon)=\nabla\cdot[(\mu+\frac{\mu_t}{\sigma_{\epsilon}})\nabla\epsilon]+\frac{\epsilon}{k}(C_{1\epsilon}G_k-C_{2\epsilon}\rho\epsilon)，其中\mu_t为湍流粘度，\sigma_k、\sigma_{\epsilon}为湍流普朗特数，C_{1\epsilon}、C_{2\epsilon}为经验常数。该模型在处理高雷诺数湍流流动时具有较好的精度和稳定性，能够较好地模拟采空区内的湍流特性。4.2物理模型及参数确定物理模型构建与网格划分：为了准确模拟采空区遗煤自燃过程，根据某煤矿1101综采工作面的实际尺寸和地质条件，构建了三维物理模型。模型包括采煤工作面、进风巷、回风巷和采空区四个部分，其中采煤工作面长度为200m，宽度为8m，高度为4.5m；进风巷和回风巷的长度均为200m，宽度为5m，高度为3m；采空区沿工作面推进方向长度为300m，宽度为200m，高度为4.5m。在构建物理模型后，对其进行网格划分。采用结构化网格划分方法，在采煤工作面、进风巷和回风巷等区域采用较密的网格，以提高计算精度；在采空区内部采用相对较疏的网格，以减少计算量。经过多次调试和优化，最终确定网格数量为500万个，网格质量满足计算要求。通过对网格无关性的验证，确保不同网格数量下的计算结果差异在合理范围内，从而保证网格划分的合理性和计算结果的准确性。2.遗煤自燃影响参数确定：煤的物性参数：煤的物性参数是影响遗煤自燃的重要因素，通过对该煤矿3号煤层的煤样进行实验测定，得到了煤的密度为1350kg/m³，比热容为1.26kJ/(kg・K)，热导率为0.25W/(m・K)。这些参数对于准确模拟煤的氧化放热、热量传递等过程至关重要，直接影响到温度场和气体浓度场的计算结果。边界条件：边界条件的设置直接影响到模拟结果的准确性。在进风巷入口处，设置为速度入口边界条件，根据实际测量的进风量，确定入口风速为4m/s，同时给定入口处的氧气浓度为0.21，温度为293K；在回风巷出口处，设置为压力出口边界条件，出口压力为101325Pa；在采煤工作面和采空区的壁面处，设置为无滑移壁面边界条件，即速度为0，同时考虑壁面与周围介质的热交换。这些边界条件的设置基于实际的通风系统和地质条件，能够真实地反映采空区内的物理过程。初始条件：初始条件是模拟计算的起点，对于采空区遗煤自燃模拟，初始时采空区内的气体浓度和温度均匀分布，氧气浓度为0.21，一氧化碳浓度为0，温度为293K。这些初始条件的设定基于采空区开采初期的实际情况，为后续的模拟计算提供了合理的基础。4.3自燃危险区域模拟分析漏风流场模拟结果分析：通过对采空区漏风流场的模拟，得到了采空区内风流速度的分布情况。从模拟结果可以看出，在进风巷和回风巷附近，风流速度较大，这是由于风流在巷道内流动时受到通风压力的作用，流速较快。而在采空区内部，风流速度逐渐减小，且分布不均匀。在采空区靠近工作面的区域，由于顶板垮落形成的裂隙较多，漏风通道较为复杂，风流速度相对较大；而在采空区深部，漏风通道相对较少，风流速度较小。漏风流场的分布对氧气的输送和热量的传递有着重要影响。在风流速度较大的区域，氧气能够更快地输送到遗煤表面，促进煤炭的氧化反应；同时，风流也能够带走部分热量，抑制煤炭温度的升高。而在风流速度较小的区域，氧气供应相对不足，煤炭氧化反应速率减缓，但热量容易积聚，增加了煤炭自燃的风险。2.氧浓度场模拟结果分析：氧浓度是煤炭自燃的关键因素之一，通过模拟得到了采空区内氧浓度的分布情况。在进风巷入口处，氧浓度为0.21，随着风流进入采空区，氧浓度逐渐降低。在采空区靠近工作面的区域，由于漏风较大，氧浓度相对较高，一般在0.15-0.20之间；而在采空区深部，氧浓度迅速降低，在距离工作面150m以外的区域，氧浓度低于0.08，进入窒息带。氧浓度的分布直接影响着煤炭氧化反应的速率和范围。在氧浓度较高的区域，煤炭氧化反应较为剧烈，容易产生热量积聚，从而引发自燃；而在氧浓度较低的区域，煤炭氧化反应受到抑制，一般不会发生自燃。因此，通过对氧浓度场的模拟分析，可以确定采空区内煤炭自燃的危险区域，为防灭火措施的制定提供依据。3.温度场模拟结果分析：温度场的变化是煤炭自燃的重要标志，模拟结果显示，在采空区初始阶段，温度分布较为均匀，与进风温度相近。随着时间的推移，在采空区靠近工作面的区域，由于煤炭氧化放热，温度逐渐升高。在距离工作面50-100m的范围内，温度升高较为明显，最高温度可达323K左右。而在采空区深部，由于氧气浓度较低，煤炭氧化反应较弱，温度升高不明显。温度场的分布与漏风流场和氧浓度场密切相关。在漏风较大、氧浓度较高的区域，煤炭氧化反应剧烈，产生的热量较多，温度升高较快；而在漏风较小、氧浓度较低的区域，煤炭氧化反应受到抑制，温度升高较慢。通过对温度场的模拟分析，可以直观地了解采空区内煤炭自燃的发展趋势，及时发现潜在的自燃危险区域。4.自燃危险区域确定：综合漏风流场、氧浓度场和温度场的模拟结果，采用氧气浓度和温度作为指标，确定采空区自燃危险区域。当氧浓度在0.08-0.18之间，且温度高于313K时，定义为自燃危险区域。根据这一标准，确定该采空区的自燃危险区域主要分布在距离工作面50-150m的范围内，呈带状分布。自燃危险区域的确定对于制定合理的防灭火措施具有重要意义。针对自燃危险区域，可以采取加强通风、封堵漏风通道、注入惰性气体等措施，降低氧浓度，带走热量，从而有效预防煤炭自燃的发生。同时，通过对自燃危险区域的监测和分析，可以及时调整防灭火措施，确保采空区的安全。五、电厂烟气注入对遗煤自燃影响的模拟分析5.1对流场及扩散影响区域分析通过数值模拟，详细研究了电厂烟气注入后采空区内流场的变化情况。在注入初期，电厂烟气从注入口迅速喷出，形成高速射流，对周围的空气产生强烈的扰动。随着烟气的扩散，采空区内的风流结构发生显著改变，原本较为规则的风流流线变得紊乱，这是由于烟气的注入打破了采空区内原有的风流平衡。在采空区的不同位置，流场的变化存在明显差异。靠近注入口的区域，烟气的流速较高，形成一个强扰动区域，该区域内的气体混合较为剧烈。随着距离注入口距离的增加，烟气的流速逐渐降低，在采空区的深部，流速趋于稳定，但仍低于未注入烟气时的风速。这是因为烟气在扩散过程中，与周围空气不断进行动量交换，导致其流速逐渐衰减。通过追踪烟气中特定示踪粒子的运动轨迹，清晰地揭示了烟气在采空区内的运移路径。烟气从注入口出发，首先沿着阻力较小的通道向前运移，遇到采空区顶板垮落形成的裂隙或空洞时，会发生分流和绕流现象。部分烟气会沿着顶板裂隙向上扩散，而另一部分则会沿着采空区底部向深部运移。在运移过程中，烟气与周围空气不断混合，其浓度逐渐降低。随着注入时间的增加，烟气的扩散范围逐渐扩大。在注入后的第1天，烟气主要集中在注入口附近半径约50m的范围内；到第3天，扩散范围扩大到半径100m左右；第7天时，烟气已基本覆盖了距离注入口200m范围内的采空区区域。这表明随着时间的推移，烟气能够逐渐扩散到采空区的更远处，对更大范围的采空区进行惰化。为了更直观地展示烟气的惰化区域，绘制了不同时刻采空区内氧气浓度的分布云图。在未注入烟气时，采空区内氧气浓度较为均匀，约为0.21。注入烟气后，在烟气扩散到的区域，氧气浓度明显降低。以注入后第7天为例，在距离注入口150m范围内，氧气浓度已降至0.1以下，形成了一个明显的惰化区域。这说明电厂烟气的注入能够有效地降低采空区内氧气浓度，抑制煤炭氧化自燃。5.2对遗煤自燃“三带”的影响为了深入探究电厂烟气含氧量以及注入位置对自燃“三带”范围和分布的影响，进行了多组模拟实验。首先，固定注入位置，改变电厂烟气的含氧量，分别设置含氧量为5%、7%、9%进行模拟。结果表明，随着电厂烟气含氧量的降低，自燃“三带”中的氧化自燃带范围明显减小。当含氧量为5%时，氧化自燃带的宽度约为50m；含氧量提高到7%时，氧化自燃带宽度增加到70m；含氧量为9%时，氧化自燃带宽度进一步扩大到90m。这是因为较低的含氧量能够更有效地抑制煤炭的氧化反应，减少氧化生热，从而缩小了氧化自燃带的范围。在含氧量较低的情况下，煤炭与氧气的接触机会减少，氧化反应速率减缓，热量产生量降低，使得氧化自燃带难以向外扩展。而含氧量较高时，煤炭氧化反应相对剧烈，产生的热量较多，有利于氧化自燃带的发展。接着，固定电厂烟气含氧量为7%，改变注入位置进行模拟。当注入位置靠近工作面时，氧化自燃带在靠近工作面一侧的范围明显减小，但在采空区深部的范围略有增加。这是因为靠近工作面注入时，烟气能够迅速稀释工作面附近的氧气浓度，抑制该区域煤炭的氧化自燃。然而，由于烟气在向采空区深部扩散过程中，受到漏风等因素的影响，使得深部区域的氧气浓度降低相对较慢，导致氧化自燃带在深部的范围有所增加。当注入位置远离工作面时，氧化自燃带在靠近工作面一侧的范围变化不大，但在采空区深部的范围明显减小。这是因为远离工作面注入时，烟气对工作面附近的氧气浓度影响较小，但在深部区域，烟气能够更有效地降低氧气浓度，抑制煤炭氧化，从而减小了氧化自燃带在深部的范围。为了量化分析电厂烟气注入对自燃“三带”的影响，计算了不同工况下氧化自燃带的面积占比。结果显示，当电厂烟气含氧量为5%，注入位置靠近工作面时，氧化自燃带面积占采空区总面积的比例为15%；含氧量提高到9%，注入位置不变时，该比例增加到25%。当含氧量固定为7%，注入位置远离工作面时，氧化自燃带面积占比为20%，而注入位置靠近工作面时，该比例为18%。这些数据进一步表明，电厂烟气含氧量和注入位置对自燃“三带”的范围和分布具有显著影响，合理调整这两个参数能够有效控制氧化自燃带的范围，降低煤炭自燃的风险。5.3对温度场影响模拟结果显示，在注入电厂烟气后，采空区内的温度场发生了明显变化。在注入初期，由于烟气的温度相对较低，注入点附近区域的温度迅速下降。随着时间的推移，这种降温效应逐渐向周围扩散。在注入后的第1天，注入点周围半径30m范围内的温度下降了约5K；到第3天，降温范围扩大到半径50m，温度下降了约8K；第7天时，半径80m范围内的温度下降了约10K。这表明电厂烟气的注入能够有效地降低采空区内的温度，抑制煤炭氧化自燃的发展。通过分析不同时刻采空区内的温度分布云图，可以清晰地看到温度场的变化趋势。在未注入烟气时，采空区内温度分布相对均匀，随着煤炭氧化反应的进行，在靠近工作面的氧化自燃带内，温度逐渐升高。注入烟气后，在烟气扩散到的区域，温度明显降低，形成一个低温区域。在低温区域内，煤炭氧化反应受到抑制，因为温度是影响煤炭氧化反应速率的重要因素之一，降低温度能够减缓氧化反应速率，减少热量产生。为了研究电厂烟气含氧量对温度的影响，设置了不同含氧量的模拟工况。当含氧量为5%时，在注入后的第7天，采空区内最高温度为305K；含氧量提高到9%时，最高温度上升到315K。这说明含氧量越高，煤炭氧化反应越剧烈，产生的热量越多，导致采空区内温度升高。因为氧气是煤炭氧化反应的氧化剂，含氧量的增加会提供更多的氧化反应动力，使反应更加剧烈，从而释放更多的热量。电厂烟气注入位置对温度也有显著影响。当注入位置靠近工作面时，工作面附近区域的温度下降更为明显。在注入后的第7天，靠近工作面注入时，工作面附近区域的平均温度为300K；而注入位置远离工作面时，该区域平均温度为305K。这是因为靠近工作面注入时，烟气能够直接作用于煤炭氧化反应较为剧烈的区域，迅速带走热量，降低温度。而远离工作面注入时，烟气在扩散过程中，热量传递效率相对较低，对工作面附近区域的降温效果不如靠近工作面注入时明显。5.4注入参数确定综合考虑防灭火效果、成本以及对采空区环境的影响等因素，确定电厂烟气注入的最佳参数。在防灭火效果方面，通过对不同注入参数下采空区内氧气浓度、温度以及自燃“三带”范围的模拟分析，发现当电厂烟气含氧量为6%-7%，注入位置距离工作面50-80m时，能够有效地降低采空区内氧气浓度和温度，缩小氧化自燃带范围，防灭火效果最佳。在该参数下，采空区内氧气浓度可降至0.08以下，温度可降低10-15K，氧化自燃带范围缩小30%-40%。从成本角度考虑，电厂烟气的净化处理成本和输送成本是主要的成本构成。降低烟气含氧量需要更复杂的净化工艺，会增加成本。因此，在保证防灭火效果的前提下，选择相对容易实现的含氧量范围。同时，注入位置的选择也会影响输送成本，距离过远会增加管路铺设成本和输送能耗。综合考虑，确定距离工作面50-80m的注入位置较为经济合理。对采空区环境的影响也是确定注入参数的重要因素。过高的注入压力可能会导致采空区顶板垮落，破坏采空区的稳定性；而注入量过大可能会影响采空区内的气体平衡，导致有害气体积聚。因此，在确定注入参数时，需要结合采空区的地质条件和实际情况，合理控制注入压力和注入量。通过模拟分析，确定注入压力为0.2-0.3MPa，注入量为1000-1500m³/h时，既能保证烟气能够顺利注入采空区，又不会对采空区环境造成不利影响。通过多因素综合分析，确定电厂烟气注入的最佳参数为：含氧量6%-7%，注入位置距离工作面50-80m，注入压力0.2-0.3MPa，注入量1000-1500m³/h。在实际应用中，可根据具体的矿井条件和需求，对这些参数进行适当调整，以达到最佳的防灭火效果和经济效益。六、结论与展望6.1研究结论通过对电厂烟气注入采空区对遗煤自燃影响的数值模拟研究，得出以下主要结论：流场及扩散影响：电厂烟气注入后显著改变了采空区内的流场结构。在注入初期，烟气形成高速射流，强烈扰动周围空气，使风流流线变得紊乱。随着烟气扩散，采空区内不同位置的流速发生明显变化，靠近注入口区域流速高、扰动强，气体混合剧烈；远离注入口处流速逐渐降低，在采空区深部趋于稳定，但低于未注入烟气时的风速。烟气从注入口出发，沿阻力小的通道运移，遇顶板裂隙或空洞会分流绕流，部分向上扩散，部分向深部运移。注入时间增加，烟气扩散范围扩大，注入第1天主要集中在注入口附近半径约50m范围，第7天基本覆盖距离注入口200m范围内区域。在烟气扩散到的区域，氧气浓度明显降低，注入后第7天，距离注入口150m范围内氧气浓度降至0.1以下，形成有效惰化区域，抑制煤炭氧化自燃。对遗煤自燃“三带”的影响：电厂烟气含氧量和注入位置对自燃“三带”范围和分布影响显著。降低烟气含氧量，氧化自燃带范围明显减小，含氧量为5%时氧化自燃带宽度约50m，含氧量提高到9%时宽度扩大到90m。靠