基于数值模拟的综放开采采空区液氮降温防灭火效能剖析

基于数值模拟的综放开采采空区液氮降温防灭火效能剖析一、引言1.1研究背景煤炭作为我国重要的基础能源，在国家能源结构中占据着举足轻重的地位。随着煤炭开采技术的不断进步，综放开采技术因其高产高效、巷道掘进率低、对地质条件适应性强等显著优势，在我国厚煤层开采中得到了广泛应用。据统计，我国厚煤层储量丰富，约占全部煤层储量的44%，而综放开采技术的应用使得厚煤层开采的效率和安全性得到了极大提升。然而，综放开采过程中也面临着诸多挑战，其中采空区火灾隐患是制约其安全生产的关键问题之一。由于综放开采时采空区丢煤较多，这些浮煤呈破碎状态，与空气接触面积大，且采空区内存在漏风现象，为煤炭自燃提供了充足的氧气和热量积聚条件。相关研究表明，全国重点煤矿中，有煤层自然发火倾向的矿井占比高达54.19%，而综放工作面采空区的煤炭自燃问题尤为突出。一旦发生火灾，不仅会造成煤炭资源的大量浪费，还会产生大量有毒有害气体，如一氧化碳、二氧化硫等，严重威胁井下作业人员的生命安全，导致人员中毒伤亡；同时，火灾还可能引发瓦斯爆炸、煤尘爆炸等次生灾害，对矿井设施造成严重破坏，使矿山生产中断，给企业带来巨大的经济损失。目前，针对采空区火灾防治，虽然已经有多种技术手段，如灌浆防灭火、均压防灭火、阻化剂防灭火等，但这些传统技术在实际应用中存在一定的局限性。例如，灌浆防灭火技术存在浆液流失、固结效果差等问题；均压防灭火技术对通风系统要求较高，实施难度较大；阻化剂防灭火技术则可能对环境造成一定的污染。液氮作为一种新型的防灭火介质，近年来在煤矿防灭火领域逐渐得到应用。液氮具有温度极低（沸点为-196℃）、性质稳定、不燃不助燃等特点，在防灭火过程中具有冷却降温、隔绝氧气、抑制爆炸等多重作用。当液氮注入采空区后，迅速汽化吸收大量热量，可使高温煤体温度急剧下降，降低煤炭自燃的可能性；同时，汽化后的氮气能够排挤采空区内的空气，降低氧气浓度，使煤炭氧化反应无法持续进行；此外，氮气还能稀释可燃气体浓度，有效抑制爆炸的发生。数值模拟作为一种先进的研究手段，能够对综放开采采空区液氮降温防灭火过程进行可视化模拟分析，深入研究其作用机理和影响因素，为实际工程应用提供科学依据和技术支持。通过建立合理的数学模型和物理模型，利用专业的数值模拟软件，可以准确地模拟液氮在采空区内的流动、扩散、汽化以及与煤体、空气之间的热质交换过程，预测采空区温度场、氧气浓度场、氮气浓度场等参数的变化规律，从而优化液氮防灭火方案，提高防灭火效果。综上所述，开展综放开采采空区液氮降温防灭火数值模拟研究具有重要的现实意义和工程应用价值，不仅有助于解决综放开采过程中的火灾隐患问题，保障煤矿安全生产，还能为液氮防灭火技术的进一步推广应用提供理论支持和技术指导。1.2研究目的与意义本研究旨在通过数值模拟的方法，深入探究综放开采采空区液氮降温防灭火的作用机理和影响因素，为实际工程应用提供科学依据和技术支持。具体而言，主要有以下几个方面的目的和意义：保障煤矿安全生产：综放开采采空区火灾是威胁煤矿安全生产的重大隐患，通过对液氮降温防灭火过程的数值模拟研究，可以准确掌握液氮在采空区内的流动、扩散规律以及对煤体温度、氧气浓度等参数的影响，从而优化液氮防灭火方案，有效预防和控制采空区火灾的发生，保障井下作业人员的生命安全，减少火灾事故对矿井设施的破坏，避免因火灾引发的瓦斯爆炸、煤尘爆炸等次生灾害，确保煤矿生产的安全稳定进行。提高煤炭资源利用率：采空区煤炭自燃不仅会造成煤炭资源的浪费，还会影响矿井的正常生产秩序。采用液氮降温防灭火技术，能够抑制煤炭自燃，减少煤炭损失，提高煤炭资源的回收率，使有限的煤炭资源得到更充分的利用，对于缓解我国能源紧张局面，保障国家能源安全具有重要意义。丰富和完善煤矿防灭火理论：目前，虽然液氮防灭火技术在煤矿领域有了一定的应用，但相关的理论研究还不够深入和系统。本研究通过数值模拟手段，对液氮降温防灭火过程中的热质交换、流动扩散等复杂物理现象进行定量分析，有助于揭示其内在作用机理，丰富和完善煤矿防灭火理论体系，为今后该技术的进一步发展和应用提供坚实的理论基础。同时，研究成果也可为其他类似矿山火灾防治提供参考和借鉴，推动整个矿山安全领域的技术进步。促进环保与可持续发展：传统的防灭火技术如灌浆、阻化剂等可能会对环境造成一定的污染，而液氮作为一种清洁、无污染的防灭火介质，在使用过程中不会产生有害物质，对环境友好。开展液氮降温防灭火数值模拟研究，推广液氮防灭火技术的应用，有利于减少煤矿开采对环境的负面影响，实现煤炭行业的绿色、可持续发展。1.3国内外研究现状1.3.1综放开采采空区火灾防治研究现状煤炭自燃是一个复杂的物理化学过程，国内外众多学者围绕其展开了深入研究。在自燃机理方面，国外早在20世纪初就开始了相关探索，提出了诸如活性基团学说、煤氧复合学说等理论。煤氧复合学说认为，煤炭自燃是煤与氧气发生物理吸附、化学吸附以及化学反应的过程，在这个过程中会产生热量，当热量积聚到一定程度且散热条件不佳时，就会引发煤炭自燃。随着研究的不断深入，国内学者在这些理论的基础上，结合我国煤矿的实际情况，进一步完善和发展了煤炭自燃理论。例如，通过热分析技术、量子化学计算等手段，对煤自燃过程中的微观结构变化、化学反应动力学等进行了研究，揭示了不同煤种自燃倾向性差异的本质原因。在综放开采采空区火灾防治技术方面，国内外已取得了丰硕的成果。传统的防治技术主要包括灌浆防灭火、均压防灭火、阻化剂防灭火等。灌浆防灭火技术是将泥浆、粉煤灰等灌浆材料与水混合后，通过管道输送到采空区，以包裹浮煤、隔绝氧气，达到预防和灭火的目的。均压防灭火技术则是通过调整通风系统，使采空区内外压力平衡，减少漏风，从而抑制煤炭自燃。阻化剂防灭火技术是向煤体喷洒或注入阻化剂，降低煤的氧化活性，延缓煤炭自燃进程。近年来，一些新型的防灭火技术也不断涌现，如凝胶防灭火、三相泡沫防灭火、氮气防灭火等。凝胶防灭火技术是利用凝胶材料在采空区形成具有一定强度和稳定性的凝胶体，包裹浮煤、封堵漏风通道，同时凝胶中的水分蒸发还能吸收热量，起到降温灭火的作用。三相泡沫防灭火技术是将气、液、固三相物质混合形成三相泡沫，注入采空区后，泡沫能在浮煤表面形成一层保护膜，隔绝氧气，且泡沫中的水分和固体颗粒能起到降温、阻化的作用。氮气防灭火技术是将氮气注入采空区，降低氧气浓度，抑制煤炭氧化反应。其中，液氮作为氮气的一种特殊形式，因其具有低温、快速降温等特性，在煤矿防灭火领域逐渐受到关注。1.3.2液氮防灭火数值模拟研究现状数值模拟技术在煤矿防灭火领域的应用越来越广泛，为研究液氮防灭火过程提供了有力的工具。国外学者率先开展了液氮防灭火数值模拟研究，利用CFD（计算流体动力学）软件对液氮在采空区的流动、扩散以及热质交换过程进行了模拟分析。通过建立物理模型和数学模型，考虑了液氮的汽化潜热、采空区的多孔介质特性、漏风等因素，模拟结果能够较为准确地反映液氮在采空区的作用效果。国内学者在液氮防灭火数值模拟方面也取得了显著进展。他们在借鉴国外研究成果的基础上，结合我国煤矿的实际地质条件和开采工艺，对数值模拟模型进行了优化和改进。例如，考虑了采空区遗煤分布的不均匀性、煤体的热物理性质随温度变化等因素，使模拟结果更加符合实际情况。同时，利用数值模拟研究了液氮注入量、注入位置、注入方式等因素对防灭火效果的影响，为液氮防灭火方案的优化提供了理论依据。然而，目前液氮防灭火数值模拟研究仍存在一些不足之处。一方面，模型中对一些复杂物理现象的描述还不够准确，如液氮与煤体之间的复杂化学反应、采空区漏风的动态变化等，这些因素可能会影响模拟结果的准确性和可靠性；另一方面，数值模拟结果与实际工程应用之间还存在一定的差距，需要进一步加强现场实测和验证，以提高数值模拟的实用性和指导性。1.4研究内容与方法1.4.1研究内容本研究围绕综放开采采空区液氮降温防灭火展开，主要内容包括：液氮降温防灭火机理研究：深入剖析液氮在采空区的物理特性变化，如汽化过程中的吸热原理，以及其与煤体、空气之间复杂的热质交换机制，明确其在降低煤体温度、隔绝氧气、抑制煤炭氧化等方面的作用路径和关键影响因素，从理论层面揭示液氮降温防灭火的内在本质。数值模型建立与验证：结合综放开采采空区的实际地质条件，包括煤层厚度、倾角、顶板岩性等，以及开采工艺参数，如采煤速度、放煤方式等，构建准确反映采空区液氮降温防灭火过程的数学模型和物理模型。运用专业的数值模拟软件进行模拟计算，并通过与现场实测数据或实验室实验数据对比，验证模型的准确性和可靠性。影响因素分析：系统研究液氮注入量、注入位置、注入时间以及采空区漏风强度、遗煤分布等因素对液氮降温防灭火效果的影响。通过改变模型中的相关参数，模拟不同工况下采空区温度场、氧气浓度场、氮气浓度场等参数的变化，分析各因素与防灭火效果之间的定量关系。方案优化：基于上述研究成果，提出针对综放开采采空区液氮降温防灭火的优化方案。确定最佳的液氮注入参数，如注入量、注入速度、注入频率等，以及合理的采空区通风控制策略，实现对采空区火灾隐患的有效防控，为实际工程应用提供科学、可行的技术指导。1.4.2研究方法为实现上述研究内容，本研究将综合运用以下方法：数值模拟方法：选用专业的CFD数值模拟软件Fluent进行研究。该软件具有强大的计算能力和丰富的物理模型库，能够准确模拟流体流动、传热传质等复杂物理过程。在Fluent软件中，通过设置合适的边界条件，如采空区进风口的风速、温度、氧气浓度，出风口的压力等，以及定义相关物理参数，如煤体的热导率、比热容，液氮的密度、汽化潜热等，对综放开采采空区液氮降温防灭火过程进行数值模拟。模拟过程中，采用有限体积法对控制方程进行离散求解，得到采空区各物理量的分布和变化情况。实验研究方法：开展实验室实验，搭建小型的综放开采采空区物理模型，模拟实际开采过程中的采空区环境。在模型中设置不同的液氮注入条件，测量煤体温度、氧气浓度、氮气浓度等参数的变化，获取实验数据。实验数据一方面用于验证数值模拟结果的准确性，另一方面为深入理解液氮降温防灭火机理提供依据。理论分析方法：运用传热学、流体力学、燃烧学等相关理论知识，对液氮降温防灭火过程中的热质交换、流动扩散等现象进行理论分析。建立数学模型，推导相关计算公式，从理论上分析各因素对防灭火效果的影响，为数值模拟和实验研究提供理论支持。二、综放开采采空区特性及火灾成因分析2.1综放开采工艺特点综放开采即综合机械化放顶煤开采，是厚煤层开采的一种高效采煤方法。其具体流程为：在沿煤层（或分段）底部布置的综采工作面中，首先采煤机进行割煤作业，将煤层底部的煤采出。割煤后，液压支架及时跟进并支护顶板，防止顶板垮落，随后推移工作面前部输送机至煤帮，完成煤炭的初步运输。接着，操作后部输送机专用千斤顶，将后部输送机相应前移。在采过1-3刀后，按照规定的放煤工艺要求，打开液压支架后部的放煤窗口，使上方已被矿山压力或辅以松动爆破等方法破碎成散体的顶煤，依靠自重经放煤窗口落下，并通过后部输送机运出。当放出煤炭中的矸石含量超过一定限度后，及时关闭放煤口，至此完成一个采煤工艺循环。综放开采具有诸多显著特点。从开采强度来看，由于将割煤和放煤工序在空间和时间上进行了有效分离，使得采煤效率大幅提高，能够实现高产高效开采。与普通采煤工艺相比，综放开采多出了放煤这一工序，增加了出煤量。例如，在一些赋存条件良好的厚煤层矿井中，综放工作面的日产煤炭量可达上万吨，年产可达数百万吨甚至更高。然而，高产带来的通风需求也更为突出，随着开采强度的增大，工作面所需通风量显著增加，以满足人员呼吸、稀释瓦斯等需求，这对矿井通风系统的能力和稳定性提出了更高的要求。遗煤量方面，虽然综放开采提高了煤炭采出率，但相较于其他开采工艺，采空区内仍会遗留较多的浮煤。在放煤过程中，由于顶煤的破碎程度、放煤工艺以及放煤口的控制等因素，难以将顶煤全部放出，总会有部分煤炭残留。这些遗煤呈破碎状态，与空气接触面积大，为煤炭自燃提供了物质基础。同时，综放开采时，采空区的冒落空间较高。随着采煤工作面的推进，顶板冒落形成较大的空间，采空区内的气体流动和热量传递较为复杂，不利于热量的散失，为煤炭自燃创造了蓄热条件。而且，在综放工作面的两道（进风顺槽和回风顺槽），由于开采过程中的顶煤垮落、煤壁片帮等原因，也会堆积较多的松散煤体，增加了煤炭自燃的风险。2.2采空区特征分析采空区作为综放开采过程中形成的特定空间，其内部结构复杂，具有一系列独特的物理特征，这些特征对煤炭自燃以及防灭火措施的实施有着至关重要的影响。从结构上看，采空区由垮落带、裂隙带和弯曲下沉带构成。垮落带位于采空区最下部，是随着采煤工作面推进，直接顶岩层在矿山压力作用下垮落形成的。垮落的岩块呈不规则堆积，碎胀系数较大，一般在1.3-1.5左右，随着时间推移和上覆岩层压力的作用，部分压实后碎胀系数可降至1.03左右。裂隙带处于垮落带上方，由于“砌体梁”结构的存在，岩体并未完全垮落，主体以岩块形式存在，在横向上连续分布，具有离层裂隙；在竖向上，由于岩体垮落下沉的不一致性，产生竖向的破断裂隙，这些裂隙为气体流动和热量传递提供了通道。弯曲下沉带在裂隙带上方，该区域岩层受开采影响较小，主要在重力作用下发生弯曲下沉，竖向破断裂隙发育较少，以横向离层裂隙为主。采空区的孔隙率分布呈现出明显的不均匀性。在自然堆积区，由于垮落岩石压实程度小，呈自然堆积状态，孔隙率最大，一般可达30%-40%，漏风风流的通过能力最强，风速较大，气体流动状态多为紊流。随着向载荷影响区过渡，采动岩层垮落，冒矸及遗煤受压，堆积体孔隙率逐渐减小，一般在15%-30%之间，漏风量也相应减小。到了压实稳定区，采空区上方岩层应力逐渐趋于稳定，煤岩堆积体的孔隙率最小，通常小于15%，且变化缓慢，不易漏风。漏风情况是采空区的一个关键特征。由于采煤工作面两端的进风口和回风口存在压力差，致使工作面后方采空区产生漏风。漏风强度与压力差、采空区孔隙率、漏风通道的长度和粗糙度等因素密切相关。在靠近工作面的区域，漏风较为严重，随着远离工作面，漏风强度逐渐减弱。据研究，在一些综采工作面，靠近工作面入口处和移架时支架交错地点的漏风量可占总漏风量的30%-50%。漏风不仅为煤炭自燃提供了充足的氧气，还影响着采空区内的热量传递和温度分布。当漏风量较小时，煤炭氧化反应速度较慢，但产生的热量不易被带走，容易积聚引发自燃；当漏风量大时，虽然能带走部分热量，但也会加速煤炭的氧化进程，增加自燃风险。气体流动在采空区内十分复杂。除了漏风引起的气体流动外，煤炭氧化、瓦斯涌出等也会导致气体成分和流动状态的变化。在自然堆积区，由于孔隙率大、风速高，气体混合较为充分；而在压实稳定区，气体流动缓慢，容易形成气体分层现象。例如，在一些采空区中，靠近顶板处可能会积聚高浓度的瓦斯，而靠近底板处则氧气浓度相对较高。温度分布方面，采空区的温度受到多种因素的综合影响。在散热带，靠近采煤作业面，漏风较大，氧气充足，煤氧反应产生的热量能及时被风流带走，温度较低，一般接近进风温度。随着向自燃带深入，漏风速度和体积逐渐减小，遗煤氧化产生的热量难以散失，温度逐渐升高，当达到煤炭自燃点时，就可能引发自燃。在窒息带，由于氧气浓度极低，煤炭氧化反应微弱，温度相对较低且变化不大。2.3火灾形成原因及危害采空区火灾的形成是一个复杂的过程，涉及多种因素的相互作用，其中遗煤氧化和漏风是两个最为关键的因素。煤炭自燃是一个复杂的物理化学过程，煤与空气中的氧气发生氧化反应是其起始点。在这个过程中，煤中的有机成分与氧气结合，产生一系列化学反应，如煤中的碳与氧气反应生成一氧化碳和二氧化碳，同时释放出热量。不同煤种的自燃倾向性存在差异，这主要与煤的变质程度、水分、灰分、挥发分以及硫分等因素密切相关。一般来说，变质程度较低的煤，其挥发分含量较高，化学活性较强，更容易发生自燃。例如，褐煤的自燃倾向性通常比无烟煤强，因为褐煤的挥发分含量可高达40%-60%，而无烟煤的挥发分含量一般在10%以下。水分在煤炭自燃过程中也起着重要作用，适量的水分可以促进煤的氧化反应，因为水可以作为催化剂，加速煤与氧气的化学反应；但过多的水分则会起到抑制作用，因为水分蒸发会吸收热量，降低煤体温度。在综放开采过程中，采空区内会遗留大量的浮煤。这些遗煤呈破碎状态，与空气接触面积大幅增加，为煤炭氧化提供了充足的物质基础。随着时间的推移，遗煤不断与氧气发生氧化反应，产生的热量逐渐积聚。当热量积聚到一定程度，且散热条件不佳时，煤体温度就会不断升高，从而引发煤炭自燃。漏风是采空区火灾形成的另一个重要因素。由于采煤工作面向前推进时，采空区顶板冒落不严实，会形成漏风通道，同时通风系统不完善、通风设施设置不合理或管理不善，也会造成风流短路或漏风。漏风为遗煤自燃提供了源源不断的氧气，加速了煤炭的氧化进程。当漏风量较小时，煤炭氧化反应速度相对较慢，但产生的热量不易被带走，容易积聚引发自燃；当漏风量大时，虽然能带走部分热量，但也会使遗煤氧化速度加快，增加自燃风险。例如，在一些采空区中，靠近工作面入口处和移架时支架交错地点的漏风较为严重，这些区域的煤炭自燃风险也相对较高。采空区火灾一旦发生，会带来诸多严重危害，对人员、资源和环境都造成巨大的威胁。对人员安全而言，火灾发生时，会消耗大量氧气，导致采空区内氧气浓度急剧下降，人员呼吸不到足够的空气，容易造成窒息。同时，燃烧会产生高温和大量有毒、有害气体，如一氧化碳、二氧化硫、氮氧化物等。一氧化碳是一种无色无味的剧毒气体，它与人体血红蛋白的结合能力比氧气强得多，一旦人体吸入一氧化碳，会迅速与血红蛋白结合，形成碳氧血红蛋白，阻碍氧气的运输和利用，导致人体缺氧中毒，严重时可危及生命。据统计，在煤矿火灾事故中，因一氧化碳中毒伤亡的人数占总伤亡人数的比例高达70%-80%。此外，火灾还可能引发瓦斯爆炸、煤尘爆炸等次生灾害，爆炸产生的强大冲击波和高温火焰，会对井下作业人员的生命安全造成直接威胁，导致人员伤亡惨重。从资源角度来看，采空区火灾会烧毁大量的煤炭资源，造成资源的巨大浪费。煤炭是不可再生资源，其储量有限，火灾导致的煤炭损失，不仅影响了煤矿企业的经济效益，也对国家的能源供应产生不利影响。同时，为了扑灭火灾，需要投入大量的人力、物力和财力，增加了企业的生产成本。例如，一些大型煤矿采空区火灾的扑救工作，可能需要持续数月甚至数年，耗费大量的资金用于灭火材料、设备和人员投入。在环境方面，采空区火灾会对周边环境造成严重污染。火灾产生的浓烟和粉尘会排放到大气中，其中含有大量的有害物质，如颗粒物、二氧化硫、氮氧化物等，这些物质会对空气质量造成严重影响，危害人体健康，引发呼吸道疾病、心血管疾病等。此外，火灾还可能导致地下水污染，燃烧产生的有害物质渗入地下，污染地下水，影响周边居民的生活用水和农业灌溉用水。而且，火灾对生态环境也会造成破坏，烧毁植被，破坏生态平衡，影响生物多样性。三、液氮降温防灭火原理及优势3.1液氮的物理性质液氮，即液态的氮气，英文名为Liquidnitrogen，常缩写为LN2。在常温常压下，它呈现为无色透明、无味无臭的液体状态，且具有良好的流动性。液氮无毒且无刺激性，化学性质十分稳定，只有在高温和放电等特殊条件下，才能够与金属或非金属发生化合反应，生成氮化物。液氮的温度极低，在常压下，其沸点为-196℃，熔点更是低至-210℃。这使得液氮成为一种极为高效的冷源，在许多领域都有着广泛的应用，如医疗领域中用于保存活体组织、外科手术中的冷冻治疗，以及食品行业的冷冻保鲜等。其密度为0.808g/cm³（-196℃），相比水的密度（1g/cm³）略小。当液氮处于-196℃时，其质量体积比为1.24m³/kg。液氮的汽化潜热为5.56kJ/mol，这意味着在液氮汽化过程中，每摩尔液氮会吸收5.56千焦的热量。这种特性使得液氮在汽化时能够吸收大量的热量，从而起到显著的冷却降温作用。例如，在火灾现场，液氮的汽化可以迅速降低周围环境的温度，有效抑制火势的蔓延。当液氮从液态转变为气态时，其体积会发生急剧膨胀。在标准状况下（0℃，101.325kPa），1m³的液氮汽化后可形成约643m³的气氮；若温度升高至25℃，1m³液氮则可形成约700m³的气氮。这种体积的大幅膨胀，使得液氮在防灭火过程中能够快速排挤周围的空气，降低氧气浓度，从而达到窒息灭火的效果。3.2降温防灭火作用机制3.2.1汽化吸热降温液氮在常压下的沸点为-196℃，当液氮注入到温度相对较高的采空区后，会迅速吸收周围环境的热量而发生汽化。这一过程遵循热力学中的汽化潜热原理，即物质从液态转变为气态时需要吸收大量的热量。每千克液氮汽化时大约需要吸收199.2kJ的热量，这些热量的来源主要是高温煤体以及周围的空气。从微观角度来看，液氮分子在获得足够的能量后，其热运动加剧，分子间的距离增大，从而实现从液态到气态的转变。在这个过程中，液氮分子与周围的煤体分子和空气分子频繁碰撞，将煤体和空气中的热量传递给自身，使得煤体和空气的温度降低。这种降温作用可以有效地降低煤炭自燃的可能性，因为煤炭自燃需要一定的温度条件，当煤体温度降低到其自燃点以下时，煤炭氧化反应的速率会显著减缓，甚至停止。例如，在某煤矿采空区火灾防治中，通过向采空区注入液氮，使得高温煤体的温度在短时间内从100℃降低到50℃以下，有效地抑制了煤炭的自燃。而且，由于液氮的汽化速度快，能够在短时间内吸收大量热量，形成局部的低温环境，对火灾的蔓延起到了很好的遏制作用。3.2.2稀释氧气浓度随着液氮的汽化，大量的氮气迅速充斥在采空区内。氮气是一种惰性气体，化学性质稳定，不支持燃烧。在采空区中，氮气的增加会稀释氧气的浓度，从而破坏煤炭自燃所需要的氧气条件。根据理想气体状态方程PV=nRT（其中P为压强，V为体积，n为物质的量，R为普适气体常量，T为温度），在采空区体积和温度相对稳定的情况下，氮气的增加会导致氧气的物质的量分数降低。当氧气浓度降低到一定程度时，煤炭氧化反应将无法持续进行。研究表明，当氧气浓度低于12%时，煤炭的氧化反应会受到明显抑制；当氧气浓度低于8%时，煤炭基本不会发生自燃。在实际应用中，通过合理控制液氮的注入量和注入位置，可以有效地降低采空区内氧气的浓度。例如，在一些综放工作面采空区，通过在采空区氧化带附近注入液氮，使得该区域的氧气浓度从原来的20%左右降低到10%以下，成功地预防了煤炭自燃。而且，由于氮气的密度与空气相近，在采空区内能够较为均匀地分布，从而实现对氧气的有效稀释。3.2.3隔绝氧气液氮汽化后形成的氮气不仅能够稀释氧气浓度，还能在煤体表面形成一层气体保护膜，起到隔绝氧气的作用。这层保护膜可以阻止氧气与煤体的直接接触，从而抑制煤炭的氧化反应。从分子层面来看，氮气分子会吸附在煤体表面，占据煤体表面的活性位点，使得氧气分子难以与煤体发生化学反应。而且，由于氮气的存在，即使有少量氧气分子能够扩散到煤体表面，也会因为氮气分子的阻挡而难以与煤体进行充分的接触和反应。这种隔绝氧气的作用类似于在煤体表面涂覆了一层防护涂层，有效地延缓了煤炭的氧化进程。例如，在实验室模拟实验中，将煤样暴露在空气中，其氧化速率较快，而在煤样表面通入氮气后，煤样的氧化速率明显降低。在实际的煤矿采空区中，通过注入液氮，在遗煤表面形成的氮气保护膜能够长期有效地隔绝氧气，减少煤炭自燃的风险。而且，随着时间的推移，即使部分氮气扩散流失，由于液氮的持续汽化补充，仍然能够维持一定的氮气浓度，保证隔绝氧气的效果。3.3与传统防灭火方法对比在煤矿采空区防灭火领域，传统的方法如灌浆防灭火、注惰气（主要为氮气）防灭火等应用历史较为悠久，各有其特点和适用场景，与液氮降温防灭火技术相比，存在诸多差异。灌浆防灭火是较为常见的传统防灭火方法之一，它是将泥浆、粉煤灰等固体材料与水混合形成的浆液，通过管道输送到采空区，以包裹浮煤、隔绝氧气，从而达到预防和灭火的目的。在实际应用中，浆液的输送需要较大的管径和较高的压力，以确保其能够到达采空区的各个部位。但由于浆液的流动性较差，在输送过程中容易出现堵管现象，影响防灭火效果。而且，灌浆后浆液的固结时间较长，一般需要数天甚至数周才能达到较好的固结效果，在这段时间内，浆液可能会流失，导致防灭火效果不佳。从降温速度来看，灌浆的降温作用主要依靠水的蒸发吸热，降温速度相对较慢。以某煤矿为例，采用灌浆防灭火时，高温区域的温度在一周内才从100℃降低到80℃左右。相比之下，液氮注入采空区后，能在短时间内迅速汽化，大量吸收热量，使高温煤体温度急剧下降。在相同的高温区域，注入液氮后，煤体温度在数小时内就可从100℃降低到50℃以下，降温速度明显更快。注惰气防灭火技术主要是利用氮气等惰性气体，通过向采空区注入惰气，降低氧气浓度，抑制煤炭氧化反应。传统的注惰气方式一般是通过制氮设备产生氮气后，经管道输送至采空区。这种方式虽然能在一定程度上降低氧气浓度，但由于氮气在常温下为气态，其在采空区的扩散速度相对较慢，难以在短时间内形成有效的惰化区域。而液氮在注入采空区后，不仅能迅速汽化产生大量氮气，稀释氧气浓度，而且由于其汽化时体积急剧膨胀，能够快速在采空区内扩散，更有效地排挤空气，降低氧气含量。从灭火效率方面来看，传统注惰气防灭火在面对已经发生的火灾时，由于其作用速度相对较慢，往往难以迅速控制火势。而液氮具有快速降温、稀释氧气和隔绝氧气的多重作用，能够更快速地扑灭火灾，提高灭火效率。例如，在某煤矿火灾事故中，采用传统注惰气防灭火方法，经过数小时才将火势初步控制，而采用液氮灭火，在短时间内就使火势得到了有效遏制。在环保方面，灌浆防灭火中使用的泥浆等材料可能会对环境造成一定的污染，如浆液中的重金属等物质可能会渗入地下，污染土壤和地下水。而液氮作为一种清洁、无污染的防灭火介质，在使用过程中不会产生有害物质，对环境友好。在操作便利性上，灌浆防灭火需要准备大量的固体材料和水，并且要对浆液进行搅拌、输送等一系列操作，过程较为繁琐。注惰气防灭火则需要配备专门的制氮设备，设备投资较大，且设备的维护和运行成本也较高。相比之下，液氮可以通过槽车运输到现场，直接注入采空区，操作相对简单便捷。综上所述，液氮降温防灭火技术在降温速度、灭火效率、环保性和操作便利性等方面相较于传统的灌浆防灭火和注惰气防灭火方法具有明显的优势，为综放开采采空区的防灭火工作提供了一种更为高效、可靠的选择。四、数值模拟理论与模型建立4.1数值模拟基本理论数值模拟作为一种强大的研究工具，在综放开采采空区液氮降温防灭火研究中发挥着关键作用，其理论基础涵盖计算流体力学、传热学等多个领域，这些理论相互交织，共同为准确模拟复杂的物理过程提供了支撑。计算流体力学（CFD）是一门通过数值计算方法求解流体流动控制方程，以研究流体流动现象的学科。在综放开采采空区液氮降温防灭火的数值模拟中，涉及到采空区内的气体流动，包括漏风、氮气的扩散以及空气与氮气的混合等过程。这些过程需要依据CFD理论进行精确描述。流体流动需遵循质量守恒方程、动量守恒方程和能量守恒方程。质量守恒方程，也被称作连续性方程，其数学表达式为\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{v})=0，其中\rho代表流体密度，t为时间，\vec{v}是速度矢量。该方程表明在单位时间内，微元体中流体质量的增加等于同一时间间隔内流入该微元体的净质量，确保了流体在流动过程中质量的守恒。动量守恒方程，即纳维-斯托克斯方程（NS方程），在笛卡尔坐标系下的表达式为：\rho(\frac{\partialu}{\partialt}+u\frac{\partialu}{\partialx}+v\frac{\partialu}{\partialy}+w\frac{\partialu}{\partialz})=-\frac{\partialp}{\partialx}+\mu(\frac{\partial^2u}{\partialx^2}+\frac{\partial^2u}{\partialy^2}+\frac{\partial^2u}{\partialz^2})+\rhof_x\rho(\frac{\partialv}{\partialt}+u\frac{\partialv}{\partialx}+v\frac{\partialv}{\partialy}+w\frac{\partialv}{\partialz})=-\frac{\partialp}{\partialy}+\mu(\frac{\partial^2v}{\partialx^2}+\frac{\partial^2v}{\partialy^2}+\frac{\partial^2v}{\partialz^2})+\rhof_y\rho(\frac{\partialw}{\partialt}+u\frac{\partialw}{\partialx}+v\frac{\partialw}{\partialy}+w\frac{\partialw}{\partialz})=-\frac{\partialp}{\partialz}+\mu(\frac{\partial^2w}{\partialx^2}+\frac{\partial^2w}{\partialy^2}+\frac{\partial^2w}{\partialz^2})+\rhof_z式中，u、v、w分别为速度在x、y、z方向上的分量，p是压力，\mu为动力粘度，f_x、f_y、f_z是质量力在x、y、z方向上的分量。此方程体现了微元体中流体动量的增加率等于作用在微元体上各种力之和，涵盖了压力、粘性力和质量力等对流体动量变化的影响，全面描述了流体动量的传输和转化过程。能量守恒方程描述了流体微团内能变化率与流入微团的净热流量以及体积力和表面力对流体微团做功功率之间的关系。其表达式为\rhoc_p(\frac{\partialT}{\partialt}+\vec{v}\cdot\nablaT)=\nabla\cdot(k\nablaT)+S，其中c_p是定压比热容，T为温度，k是热导率，S为源项，包含了化学反应热、辐射热等各种热源。该方程反映了在采空区中，热量通过传导、对流等方式进行传递，以及各种热源对温度分布的影响，对于研究液氮汽化吸热导致的采空区温度变化至关重要。在实际的采空区环境中，流体流动通常处于湍流状态，这使得流动更加复杂。为了准确模拟湍流流动，常用的湍流模型有标准k-\epsilon模型、RNGk-\epsilon模型和Realizablek-\epsilon模型等。标准k-\epsilon模型通过求解湍动能k和湍动耗散率\epsilon的输运方程来封闭NS方程，能够较好地模拟一般的湍流流动。其湍动能k方程为\frac{\partial(\rhok)}{\partialt}+\frac{\partial(\rhoku_i)}{\partialx_i}=\frac{\partial}{\partialx_j}[(\mu+\frac{\mu_t}{\sigma_k})\frac{\partialk}{\partialx_j}]+G_k-\rho\epsilon，湍动耗散率\epsilon方程为\frac{\partial(\rho\epsilon)}{\partialt}+\frac{\partial(\rho\epsilonu_i)}{\partialx_i}=\frac{\partial}{\partialx_j}[(\mu+\frac{\mu_t}{\sigma_{\epsilon}})\frac{\partial\epsilon}{\partialx_j}]+\frac{C_{1\epsilon}\epsilon}{k}G_k-C_{2\epsilon}\rho\frac{\epsilon^2}{k}，其中G_k是由平均速度梯度引起的湍动能产生项，\mu_t是湍动粘度，C_{1\epsilon}、C_{2\epsilon}、\sigma_k、\sigma_{\epsilon}是模型常数。RNGk-\epsilon模型在标准k-\epsilon模型的基础上，考虑了湍流的旋转和曲率效应，对复杂流动的模拟精度更高。Realizablek-\epsilon模型则在湍动粘度的计算和湍动耗散率方程中进行了改进，使其能够更好地模拟具有强旋流和分离流动的情况。传热学是研究热量传递规律的学科，在综放开采采空区液氮降温防灭火过程中，涉及到多种传热方式。热传导是指由于物体内部分子、原子和电子等微观粒子的热运动而产生的热量传递现象。在采空区的煤体和岩石中，热传导起着重要作用。其遵循傅里叶定律，表达式为q=-k\nablaT，其中q是热流密度，负号表示热量传递方向与温度梯度方向相反。热对流是指流体各部分之间发生相对位移而引起的热量传递过程，在采空区内，由于漏风、氮气的流动等，热对流是热量传递的主要方式之一。热对流可分为自然对流和强制对流，自然对流是由于流体内部温度不均匀导致密度差异而引起的对流，而强制对流则是在外部动力（如风机等）作用下产生的对流。热辐射是物体通过电磁波传递能量的过程，在高温环境下，热辐射的作用不可忽视。在采空区火灾发生时，高温煤体和火焰会向外辐射热量，影响周围环境的温度分布。斯蒂芬-玻尔兹曼定律描述了黑体的热辐射能力，表达式为E_b=\sigmaT^4，其中E_b是黑体的辐射力，\sigma是斯蒂芬-玻尔兹曼常数，T是黑体的绝对温度。对于实际物体，其辐射力还需要考虑发射率\epsilon的影响，即E=\epsilon\sigmaT^4。在综放开采采空区液氮降温防灭火的数值模拟中，这些计算流体力学和传热学的理论相互关联、相互影响。例如，液氮在采空区的流动会带动热量的传递，而热量传递又会影响流体的密度和粘度，进而改变流体的流动状态。因此，在建立数值模型时，需要综合考虑这些理论，准确描述各种物理现象，以获得可靠的模拟结果。4.2采空区物理模型构建为了准确模拟综放开采采空区液氮降温防灭火过程，需依据某煤矿综放工作面的实际开采条件，构建三维几何模型。该工作面的开采煤层厚度为6-8m，平均厚度7m，采用综放开采工艺，采煤机割煤高度为3m，放煤高度为4m。工作面长度为200m，推进长度为1500m。采空区顶板为砂岩，厚度约10m，底板为泥岩，厚度约5m。在构建模型时，进行了适当的简化处理。忽略采空区内一些微小的地质构造和局部不规则性，将采空区视为一个规则的长方体空间。同时，将采煤工作面简化为一个平面，不考虑采煤机、液压支架等设备的具体结构，仅考虑其对风流和气体分布的宏观影响。这样的简化处理既能突出主要物理过程，又能减少计算量，提高模拟效率。模型的边界条件设置对于模拟结果的准确性至关重要。在进风口边界，设定为速度入口边界条件，根据实际通风情况，将风速设置为4m/s，温度为25℃，氧气浓度为21%。在出风口边界，采用压力出口边界条件，出口压力设定为101325Pa。采空区的四周和顶部设置为壁面边界条件，壁面无滑移，即速度为零。对于液氮注入口，设置为质量流量入口边界条件，根据不同的模拟工况，调整液氮的注入质量流量。例如，在研究液氮注入量对防灭火效果的影响时，分别设置液氮注入质量流量为50kg/s、100kg/s、150kg/s等不同数值。在采空区与煤层的交界面，考虑煤体的透气性和瓦斯涌出等因素，设置相应的源项和边界条件。例如，根据煤体的瓦斯含量和涌出规律，设置瓦斯涌出源项，以准确模拟采空区内瓦斯浓度的分布和变化。通过合理设置这些边界条件，能够更真实地反映采空区的实际物理过程，为后续的数值模拟分析提供可靠的基础。4.3数学模型建立在综放开采采空区液氮降温防灭火的数值模拟中，建立准确的数学模型是核心任务之一，这涉及到多个物理过程的数学描述，包括质量守恒、动量守恒、能量守恒以及组分输运等，同时还需考虑煤氧反应动力学模型，以全面反映采空区内复杂的物理化学现象。质量守恒方程，即连续性方程，是描述流体流动过程中质量守恒的基本方程。在采空区的复杂环境中，流体包括空气、氮气以及其他可能存在的气体成分。对于可压缩流体，其质量守恒方程的一般形式为\frac{\partial\rho}{\partialt}+\nabla\cdot(\rho\vec{v})=0。其中，\rho表示流体的密度，它会随着温度、压力以及气体成分的变化而改变。在采空区中，由于液氮的注入和煤氧反应等过程，气体的温度和成分会发生动态变化，从而导致密度的改变。例如，液氮汽化后，氮气的密度与空气不同，会影响整个采空区内的密度分布。t代表时间，反映了流体流动过程的动态特性。随着时间的推移，采空区内的气体流动状态、温度分布以及组分浓度都会发生变化。\vec{v}是速度矢量，描述了流体在空间中的运动方向和速度大小。在采空区，由于漏风、液氮注入等因素，气体的速度分布十分复杂，在不同区域和方向上都存在差异。动量守恒方程，即纳维-斯托克斯方程（NS方程），是描述流体动量变化的关键方程。在笛卡尔坐标系下，其表达式为：\rho(\frac{\partialu}{\partialt}+u\frac{\partialu}{\partialx}+v\frac{\partialu}{\partialy}+w\frac{\partialu}{\partialz})=-\frac{\partialp}{\partialx}+\mu(\frac{\partial^2u}{\partialx^2}+\frac{\partial^2u}{\partialy^2}+\frac{\partial^2u}{\partialz^2})+\rhof_x\rho(\frac{\partialv}{\partialt}+u\frac{\partialv}{\partialx}+v\frac{\partialv}{\partialy}+w\frac{\partialv}{\partialz})=-\frac{\partialp}{\partialy}+\mu(\frac{\partial^2v}{\partialx^2}+\frac{\partial^2v}{\partialy^2}+\frac{\partial^2v}{\partialz^2})+\rhof_y\rho(\frac{\partialw}{\partialt}+u\frac{\partialw}{\partialx}+v\frac{\partialw}{\partialy}+w\frac{\partialw}{\partialz})=-\frac{\partialp}{\partialz}+\mu(\frac{\partial^2w}{\partialx^2}+\frac{\partial^2w}{\partialy^2}+\frac{\partial^2w}{\partialz^2})+\rhof_z其中，u、v、w分别为速度在x、y、z方向上的分量，它们决定了流体在三维空间中的运动轨迹和速度分布。在采空区，由于顶板垮落形成的不规则空间、漏风通道的复杂性以及液氮注入的扰动，不同方向上的速度分量会相互影响，使得流体运动呈现出复杂的形态。p是压力，它在采空区内的分布与流体的流动状态密切相关。压力差会驱动流体流动，例如，采煤工作面进风口和出风口之间的压力差会导致采空区内的漏风现象。同时，压力的变化也会影响气体的密度和温度。\mu为动力粘度，它反映了流体内部的粘性力，影响着流体的流动阻力。在采空区，气体的动力粘度会受到温度和成分的影响，例如，温度升高时，气体分子的热运动加剧，动力粘度会发生变化。f_x、f_y、f_z是质量力在x、y、z方向上的分量，常见的质量力如重力。在采空区，重力会对气体的流动产生一定的影响，特别是在一些较大的采空区空间中，气体在重力作用下会发生自然对流。能量守恒方程用于描述采空区内能量的传递和转化。其表达式为\rhoc_p(\frac{\partialT}{\partialt}+\vec{v}\cdot\nablaT)=\nabla\cdot(k\nablaT)+S。其中，c_p是定压比热容，它表示单位质量的流体在定压条件下温度升高1℃所吸收的热量。不同气体的定压比热容不同，在采空区中，随着液氮的注入和煤氧反应的进行，气体成分发生变化，定压比热容也会相应改变。T为温度，它是采空区状态的重要参数，温度的分布和变化直接反映了能量的传递和积聚情况。例如，液氮汽化时吸收大量热量，会导致周围区域温度降低，而煤氧反应产生的热量则会使煤体和周围气体温度升高。k是热导率，它决定了热量在介质中的传导能力。在采空区，煤体和岩石的热导率与气体不同，这会影响热量在不同介质之间的传递速度。S为源项，包含了化学反应热、辐射热等各种热源。煤氧反应会释放大量的化学反应热，这是采空区温度升高的重要原因之一；而在高温环境下，热辐射也会成为能量传递的重要方式，例如，高温煤体和火焰会向外辐射热量，影响周围环境的温度分布。组分输运方程用于描述采空区内各种气体组分浓度的变化。对于第i种组分，其组分输运方程为\frac{\partial(\rhoY_i)}{\partialt}+\nabla\cdot(\rho\vec{v}Y_i)=\nabla\cdot(\rhoD_{i,m}\nablaY_i)+S_i。其中，Y_i是第i种组分的质量分数，它反映了该组分在混合气体中的相对含量。在采空区，随着液氮的注入和煤氧反应的进行，氧气、氮气、一氧化碳等气体组分的质量分数会发生变化。D_{i,m}是第i种组分在混合气体中的扩散系数，它决定了该组分在混合气体中的扩散速度。不同气体组分的扩散系数不同，这会影响它们在采空区内的分布和混合情况。S_i是第i种组分的源项，它包括化学反应产生或消耗该组分的速率。例如，在煤氧反应中，氧气被消耗，一氧化碳等产物生成，这些过程都会通过源项S_i反映在组分输运方程中。煤氧反应动力学模型是描述煤炭与氧气发生化学反应的模型，它对于理解煤炭自燃过程至关重要。常用的煤氧反应动力学模型基于Arrhenius定律，该定律表明化学反应速率与温度之间存在指数关系。煤氧反应的速率方程可以表示为R=Ae^{-\frac{E}{RT}}C_O^n。其中，R是煤氧反应速率，它决定了煤炭氧化的快慢程度。反应速率受到多种因素的影响，如温度、氧气浓度和煤的性质等。A是指前因子，它与煤的化学结构和反应活性有关，不同煤种的指前因子不同，反映了它们氧化反应活性的差异。E是活化能，它表示化学反应发生所需要克服的能量障碍。活化能越高，反应越难以进行，煤的自燃倾向性就越低。R是气体常数，它是一个固定值。T是绝对温度，温度对煤氧反应速率的影响非常显著，随着温度升高，分子热运动加剧，反应速率会迅速增加。C_O是氧气浓度，氧气是煤氧反应的反应物，其浓度的高低直接影响反应速率。n是反应级数，它反映了反应速率对氧气浓度的敏感程度。一般来说，煤氧反应的反应级数在0.5-2之间，不同煤种和反应条件下反应级数会有所不同。这些控制方程和煤氧反应动力学模型相互关联、相互影响，共同构成了综放开采采空区液氮降温防灭火数值模拟的数学模型。在实际模拟过程中，需要根据采空区的具体条件和物理过程，对这些方程进行合理的简化和求解，以获得准确的模拟结果。4.4模型求解与验证在完成采空区物理模型和数学模型的构建后，选用Fluent软件进行数值求解。Fluent是一款功能强大的CFD软件，广泛应用于各种复杂流动和传热问题的模拟。它采用有限体积法对控制方程进行离散求解，将计算区域划分为一系列控制体积，使每个控制体积都包围一个节点，通过对控制方程在每个控制体积上进行积分，得到离散方程，进而求解得到各物理量在节点上的值。在Fluent软件中，首先对构建好的采空区三维几何模型进行网格划分。网格划分的质量直接影响计算结果的准确性和计算效率，因此采用非结构化网格对模型进行划分，以更好地适应采空区复杂的几何形状。在网格划分过程中，对液氮注入口、采空区氧化带等关键区域进行局部加密，确保这些区域的计算精度。例如，在液氮注入口附近，将网格尺寸设置为0.1m，而在远离注入口的区域，网格尺寸可适当增大至0.5m。通过这种局部加密的方式，既能保证关键区域的计算精度，又能有效控制计算量。对于湍流模型，经过对比分析，选用Realizablek-\epsilon模型。该模型在湍动粘度的计算和湍动耗散率方程中进行了改进，能够更好地模拟采空区复杂的湍流流动。在采空区中，由于漏风、液氮注入等因素，流体流动处于湍流状态，Realizablek-\epsilon模型能够准确地描述湍流的特性，如湍动能的产生和耗散、湍流的扩散等。例如，在模拟采空区漏风时，该模型能够准确地预测漏风的速度分布和湍流强度，为研究漏风对煤炭自燃和液氮防灭火效果的影响提供了可靠的依据。在求解过程中，采用SIMPLE算法进行压力-速度耦合求解。SIMPLE算法是一种常用的求解不可压缩流体流动问题的算法，它通过迭代求解压力修正方程和动量方程，逐步逼近真实的压力和速度分布。在采空区液氮降温防灭火的数值模拟中，压力和速度的准确求解对于描述气体的流动和传热过程至关重要。例如，在模拟液氮注入采空区的过程中，SIMPLE算法能够准确地计算液氮的喷射速度和压力分布，以及液氮与周围气体的相互作用。同时，为了提高计算精度和收敛速度，对动量方程、能量方程和组分输运方程等采用二阶迎风差分格式进行离散。二阶迎风差分格式能够更准确地描述物理量在空间上的变化，减少数值误差，提高计算结果的准确性。例如，在计算采空区内的温度分布时，采用二阶迎风差分格式能够更精确地捕捉温度梯度的变化，使模拟结果更接近实际情况。为了验证所建立模型的准确性，将数值模拟结果与某煤矿综放工作面的现场监测数据进行对比。该煤矿在采空区布置了多个温度传感器和气体浓度传感器，实时监测采空区内的温度和氧气浓度变化。选取监测数据中液氮注入前后的时间段，将模拟得到的温度场和氧气浓度场与实际监测数据进行对比分析。在温度对比方面，以距离工作面不同位置处的温度为例。在距离工作面50m处，监测数据显示液氮注入前采空区温度为35℃，注入液氮一段时间后，温度降至25℃。模拟结果显示，液氮注入前该位置温度为34.5℃，注入后降至24.8℃。模拟值与监测值的相对误差在合理范围内，温度变化趋势也与实际情况相符。从温度分布云图来看，模拟结果能够准确地反映出液氮注入后采空区内温度的降低情况，以及低温区域的扩展范围。在氧气浓度对比方面，同样选取距离工作面不同位置进行分析。在距离工作面80m处，监测数据表明液氮注入前氧气浓度为18%，注入后降至10%。模拟结果显示，注入前氧气浓度为17.8%，注入后降至9.6%。模拟值与监测值的相对误差较小，氧气浓度的变化趋势也与实际监测结果一致。通过对不同位置处温度和氧气浓度的对比分析，验证了所建立的数值模型能够较为准确地模拟综放开采采空区液氮降温防灭火过程，为后续的研究和工程应用提供了可靠的基础。五、数值模拟结果与分析5.1采空区流场分布特征利用Fluent软件对综放开采采空区液氮降温防灭火过程进行数值模拟，得到了采空区流场的分布情况，包括风流速度和压力分布。图1展示了采空区风流速度分布云图，从图中可以清晰地看出，在进风口附近，风流速度较大，达到了设定的入口风速4m/s。这是因为新鲜风流从进风口进入采空区，尚未受到太多阻碍，能够保持较高的速度。随着风流向采空区内部流动，由于受到垮落岩石和遗煤的阻挡，以及与周围介质的摩擦作用，风流速度逐渐减小。在靠近工作面的区域，风流速度相对较大，这是因为该区域通风条件较好，漏风现象较为明显。而在采空区深部，风流速度明显降低，部分区域的风速甚至接近于零。这是由于深部区域距离进风口较远，风流在流动过程中能量逐渐消耗，且受到采空区复杂地形的影响，通风阻力增大，导致风速减小。在采空区的边缘部分，由于与周围岩体的摩擦和边界效应，风流速度也相对较低。为更直观地展示风流速度在采空区内的变化规律，选取了一条沿采空区长度方向的中心线，绘制了该中心线上的风流速度变化曲线，如图2所示。从曲线中可以看出，在进风口处，风流速度迅速下降，这是由于风流进入采空区后，突然从相对开阔的巷道进入到充满垮落岩石和遗煤的采空区空间，通风断面减小，风流受到较大的阻力，导致速度急剧降低。随着向采空区内部深入，风流速度逐渐趋于稳定，但仍然呈现出缓慢下降的趋势。在距离进风口约50m处，风流速度下降到约1m/s。在采空区深部，风流速度下降更为明显，在距离进风口150m处，风流速度已降至0.2m/s左右。这表明采空区内部的通风状况较差，漏风主要集中在靠近工作面和进风口的区域。图3为采空区压力分布云图，从图中可以看出，进风口处压力最高，设定为相对压力0Pa。这是因为进风口是新鲜风流的入口，风机提供的能量使得风流具有较高的压力。随着风流在采空区内流动，由于克服通风阻力做功，压力逐渐降低。在出风口处，压力最低，为负压状态。这是由于出风口与大气相通，风流流出采空区时，压力逐渐降低至大气压力以下。在采空区内部，压力分布呈现出不均匀的特点。在靠近进风口和出风口的区域，压力变化较为明显，而在采空区中部，压力变化相对较小。这是因为中部区域受到进风口和出风口的影响相对较小，通风阻力分布较为均匀。同样选取沿采空区长度方向的中心线，绘制该中心线上的压力变化曲线，如图4所示。从曲线中可以看出，进风口处压力迅速下降，这是由于风流进入采空区后，通风阻力迅速增大，导致压力快速降低。在距离进风口约30m处，压力下降了约100Pa。随着向采空区内部深入，压力下降速度逐渐减缓，在采空区中部，压力变化较为平缓。在距离进风口100-120m区域，压力基本保持稳定，变化幅度在10Pa以内。在靠近出风口处，压力再次迅速下降，这是因为风流在流出采空区时，需要克服较大的阻力，导致压力急剧降低。在出风口处，压力降至约-150Pa。采空区的漏风情况与风流速度和压力分布密切相关。由于进风口和出风口之间存在压力差，使得采空区内形成了漏风通道。从风流速度分布云图和压力分布云图可以看出，漏风主要集中在靠近工作面和进风口的区域。这些区域风流速度较大，压力变化明显，说明漏风较为严重。而在采空区深部和边缘部分，漏风相对较小。漏风不仅为煤炭自燃提供了充足的氧气，还会影响采空区内的温度分布和氮气的扩散。因此，准确掌握采空区流场分布特征，对于分析煤炭自燃风险和优化液氮防灭火方案具有重要意义。5.2温度场变化规律在对采空区流场分布特征进行分析后，进一步研究采空区温度场的变化规律。图5展示了注液氮前采空区的温度分布云图，从图中可以看出，在靠近工作面的区域，由于煤炭开采活动以及风流的影响，温度相对较低，约为25℃，这是因为新鲜风流不断带入低温空气，带走了部分热量。随着向采空区内部深入，温度逐渐升高，在采空区深部，温度最高可达40℃左右。这是由于采空区内遗煤与氧气发生氧化反应，产生热量，且深部区域通风条件较差，热量难以散失，导致温度积聚升高。在采空区的边缘部分，温度也相对较高，这是因为边缘部分与周围岩体的热交换相对较弱，热量容易积累。图6为注液氮1h后采空区的温度分布云图。可以明显观察到，在液氮注入口附近，温度急剧下降，形成了一个低温区域，最低温度降至-100℃以下。这是由于液氮注入后迅速汽化，吸收大量热量，使得周围区域温度大幅降低。随着距离注入口的增加，温度逐渐升高，但与注氮前相比，整体温度仍有显著下降。在距离注入口50m范围内，温度基本都降至20℃以下。在采空区的其他区域，由于氮气的扩散和热传递作用，温度也有所降低，不过降低幅度相对较小。为了更直观地了解温度随时间和空间的变化趋势，选取了距离工作面不同位置的三个监测点，分别为A点（距离工作面30m）、B点（距离工作面80m）和C点（距离工作面130m），绘制了这三个点的温度随时间变化曲线，如图7所示。从曲线中可以看出，在注液氮前，三个监测点的温度基本保持稳定，A点温度约为28℃，B点温度约为32℃，C点温度约为38℃。当开始注液氮后，A点温度迅速下降，在1h内就降至10℃以下。这是因为A点距离注入口较近，受到液氮的冷却作用明显。B点温度在注氮后也逐渐下降，但下降速度相对较慢，在2h后降至20℃左右。C点温度下降更为缓慢，在3h后才降至30℃左右。这表明随着距离注入口的增加，液氮的冷却效果逐渐减弱。从空间分布来看，沿着采空区长度方向，温度呈现出逐渐降低的趋势。在注氮初期，低温区域主要集中在注入口附近，随着时间的推移，低温区域逐渐向采空区深部扩展。在注氮3h后，距离注入口100m范围内的温度都得到了有效控制，基本降至25℃以下。这说明液氮的冷却作用能够在一定范围内有效地降低采空区温度，抑制煤炭自燃。从温度分布的不均匀性来看，在注氮过程中，采空区内温度分布存在明显的梯度变化。在低温区域与高温区域的交界处，温度梯度较大，这是由于液氮的冷却作用与煤炭氧化放热之间的相互作用导致的。随着注氮时间的延长，温度梯度逐渐减小，说明采空区内的温度分布逐渐趋于均匀。采空区温度场的变化规律与流场分布密切相关。风流的流动会影响热量的传递和扩散，而液氮的注入也会改变流场的分布，进而影响温度场。例如，在风流速度较大的区域，热量能够更快地被带走，温度相对较低；而在风流速度较小的区域，热量容易积聚，温度较高。液氮注入后，由于其汽化产生的气流作用，会改变采空区内的风流方向和速度，从而影响温度的分布和变化。5.3氧气浓度分布图8为采空区氧气浓度分布云图，从中可以清晰地看出采空区氧气浓度的分布特征。在进风口区域，氧气浓度保持在初始设定的21%，这是因为新鲜风流不断从进风口涌入，带来充足的氧气。随着风流在采空区内流动，氧气浓度逐渐发生变化。在靠近工作面的区域，由于通风条件相对较好，风流速度较大，氧气能够及时补充，氧气浓度仍维持在较高水平，约为18%-20%。然而，在采空区深部，由于风流速度减小，氧气的补充受到限制，同时遗煤与氧气发生氧化反应，消耗了大量氧气，使得氧气浓度显著降低。在距离工作面100m以外的深部区域，氧气浓度降至10%以下，部分区域甚至低于5%。为进一步分析氧气浓度的变化规律，选取了沿采空区长度方向的中心线，绘制该中心线上的氧气浓度变化曲线，如图9所示。从曲线中可以看出，在进风口处，氧气浓度为21%，随着向采空区内部推进，氧气浓度逐渐下降。在距离进风口50m处，氧气浓度降至15%左右。在距离进风口100m处，氧气浓度已降至8%左右。这表明随着距离进风口距离的增加，氧气的供给逐渐减少，而煤炭氧化对氧气的消耗相对稳定，导致氧气浓度不断降低。氧气浓度对煤炭氧化和自燃有着至关重要的影响。煤炭氧化是一个复杂的化学反应过程，其反应速率与氧气浓度密切相关。当氧气浓度较高时，煤炭与氧气的接触机会增多，氧化反应速率加快，产生的热量也相应增加。研究表明，在氧气浓度为21%的条件下，煤炭的氧化反应速率相对较快，煤体温度上升明显。随着氧气浓度的降低，煤炭氧化反应速率逐渐减缓。当氧气浓度降至12%以下时，煤炭氧化反应受到明显抑制。这是因为氧气浓度降低，煤炭与氧气的接触概率减小，反应活性降低，从而使得氧化反应难以持续进行。当氧气浓度低于8%时，煤炭基本不会发生自燃。在采空区深部，由于氧气浓度较低，煤炭氧化反应微弱，煤体温度相对稳定，不易发生自燃。而在靠近工作面和进风口的区域，氧气浓度相对较高，煤炭氧化反应较为活跃，若热量不能及时散发，就容易引发煤炭自燃。因此，通过控制采空区氧气浓度，可以有效预防煤炭自燃的发生。5.4液氮运移特性通过数值模拟，深入研究了液氮在采空区的运移特性，包括其扩散、气化和分布规律，这对于理解液氮降温防灭火的作用机制具有重要意义。从扩散规律来看，液氮注入采空区后，在初始阶段，由于注入压力的作用，液氮以较高的速度向周围扩散。随着扩散距离的增加，注入压力逐渐减弱，液氮的扩散速度也逐渐减小。同时，采空区内的风流对液氮的扩散也有显著影响。在风流速度较大的区域，液氮会被风流携带，向风流方向扩散。例如，在进风口附近，风流速度较大，液氮会随着风流迅速向采空区内部扩散，使得该区域的氮气浓度迅速增加。而在风流速度较小的区域，液氮的扩散主要依靠分子扩散作用，扩散速度相对较慢。在采空区深部，风流速度较小，液氮的扩散范围相对较窄。液氮的气化过程是其发挥降温防灭火作用的关键环节。在注入采空区后，液氮迅速吸收周围环境的热量而发生气化。气化过程主要集中在液氮注入口附近，随着距离注入口的增加，气化速率逐渐降低。这是因为注入口附近的温度相对较高，液氮与周围介质的温差较大，热量传递速度快，有利于液氮的气化。在注入口周围10m范围内，液氮的气化速率较快，大部分液氮在这个区域内完成气化。随着气化的进行，液氮吸收大量热量，使得周围区域温度急剧下降，形成一个低温区域。在低温区域内，液氮的气化速率会受到一定影响，因为温度降低后，热量传递速度减缓，气化所需的热量供应相对不足。在分布规律方面，液氮注入采空区后，氮气在采空区内呈现出不均匀的分布状态。在注入口附近，氮气浓度最高，随着距离注入口的增加，氮气浓度逐渐降低。在距离注入口30m处，氮气浓度可达到50%以上；而在距离注入口100m处，氮气浓度降至20%以下。在采空区的不同高度上，氮气浓度也存在差异。由于氮气的密度略小于空气，在重力作用下，氮气有向上运动的趋势，因此在采空区的上部，氮气浓度相对较高。在距离顶板5m范围内，氮气浓度比距离底板5m范围内的氮气浓度高10%-20%。影响液氮运移的因素众多。注入压力是一个重要因素，注入压力越大，液氮的初始扩散速度越快，能够更快地到达采空区的各个部位。当注入压力从0.5MPa增加到1MPa时，液氮在1h内的扩散距离从30m增加到50m。采空区的孔隙率也对液氮运移有显著影响。孔隙率越大，采空区内的气体流动阻力越小，液氮的扩散和气化速度越快。在孔隙率为30%的区域，液氮的扩散速度比孔隙率为15%的区域快50%左右。漏风强度同样不容忽视，漏风会改变采空区内的风流场，从而影响液氮的运移。当漏风强度增大时，风流对液氮的携带作用增强，液氮会随着风流更快地扩散，但同时也可能导致液氮在某些区域的浓度分布不均匀。在漏风强度为0.5m/s的情况下，液氮在风流方向上的扩散距离比漏风强度为0.2m/s时增加了20m。综上所述，液氮在采空区的运移特性受到多种因素的综合影响，深入研究这些特性和影响因素，对于优化液氮防灭火方案，提高防灭火效果具有重要的指导意义。六、影响液氮降温防灭火效果的因素分析6.1液氮注入参数的影响液氮注入参数包括注入量、注入速度和注入位置，这些参数对综放开采采空区液氮降温防灭火效果有着显著的影响，通过数值模拟对这些因素进行深入分析，有助于确定最优的液氮注入方案，提高防灭火效果。6.1.1注入量对防灭火效果的影响为研究注入量的影响，在其他条件保持不变的情况下，设置不同的液氮注入量进行数值模拟。分别设定注入量为500kg/h、1000kg/h和1500kg/h。模拟结果表明，注入量对采空区温度和氧气浓度有着明显的影响。当注入量为500kg/h时，采空区内的温度下降速度相对较慢。在注氮1h后，距离注入口30m范围内的温度可降至30℃左右，而在距离注入口较远的区域，温度下降幅度较小。这是因为注入量较少，液氮汽化吸收的热量有限，无法对整个采空区起到有效的降温作用。同时，氧气浓度的降低也较为缓慢，在注氮2h后，采空区内大部分区域的氧气浓度仍在15%以上。当注入量增加到1000kg/h时，采空区的降温效果明显增强。注氮1h后，距离注入口50m范围内的温度可降至20℃以下，且随着时间的推移，低温区域逐渐扩大。在注氮3h后，距离注入口80m范围内的温度都得到了有效控制。氧气浓度也下降得更快，在注氮2h后，采空区内大部分区域的氧气浓度降至12%以下，有效地抑制了煤炭氧化反应。当注入量达到1500kg/h时，采空区的温度迅速下降。注氮1h后，距离注入口70m范围内的温度可降至10℃以下，整个采空区的温度都得到了显著降低。氧气浓度在短时间内大幅下降，注氮2h后，采空区内大部分区域的氧气浓度降至8%以下，煤炭基本不会发生自燃。通过对比不同注入量下的模拟结果，可以看出注入量越大，采空区的降温速度越快，氧气浓度降低越明显，防灭火效果越好。然而，注入量并非越大越好，过大的注入量不仅会增加成本，还可能对采空区的其他参数产生不利影响。在实际应用中，需要根据采空区的具体情况，如采空区体积、遗煤量、漏风情况等，综合考虑确定合适的液氮注入量。例如，对于体积较大、遗煤较多且漏风严重的采空区，需要适当增加注入量，以确保防灭火效果；而对于体积较小、遗煤较少的采空区，过大的注入量可能会造成资源浪费。6.1.2注入速度对防灭火效果的影响注入速度也是影响液氮降温防灭火效果的重要因素之一。在模拟过程中，固定液氮注入量为1000kg/h，分别设置注入速度为0.5m/s、1m/s和1.5m/s进行模拟。当注入速度为0.5m/s时，液氮在采空区内的扩散速度较慢。由于注入速度低，液氮在注入口附近积聚，难以快速扩散到采空区的各个部位。在注氮1h后，低温区域主要集中在注入口周围20m范围内，对采空区其他区域的降温效果不明显。同时，氧气浓度的降低也局限在注入口附近，采空区大部分区域的氧气浓度仍然较高。当注入速度提高到1m/s时，液氮的扩散速度加快，能够在较短时间内扩散到更广泛的区域。注氮1h后，低温区域扩大到注入口周围40m范围内，采空区内更多区域的温度得到了有效降低。氧气浓度的降低范围也相应扩大，采空区内大部分区域的氧气浓度有所下降，煤炭氧化反应受到一定程度的抑制。当注入速度达到1.5m/s时，液氮迅速扩散到采空区的各个角落。注氮1h后，低温区域覆盖了注入口周围60m的范围，整个采空区的温度都有明显下降。氧气浓度在采空区内均匀降低，注氮2h后，大部分区域的氧气浓度降至10%以下，有效地遏制了煤炭自燃。由此可见，注入速度越快，液氮在采空区内的扩散速度越快，能够更快地覆盖采空区，使采空区的温度和氧气浓度得到更均匀的降低，防灭火效果更好。但注入速度过大也可能导致液氮在采空区内的停留时间过短，无法充分发挥其降温防灭火作用。在实际应用中，需要根据采空区的尺寸、形状以及通风条件等因素，合理选择注入速度。例如，对于通风条件较好、采空区空间较大的情况，可以适当提高注入速度，以加快液氮的扩散；而对于通风条件较差、采空区较为狭窄的区域，过高的注入速度可能会使液氮迅速流出采空区，此时应适当降低注入速度。6.1.3注入位置对防灭火效果的影响注入位置的选择对液氮降温防灭火效果同样至关重要。在模拟中，设置了三种不同的注入位置：靠近工作面的进风口、采空区中部和采空区深部。当注入位置靠近工作面的进风口时，由于进风口处风流