基于数值模拟的综掘工作面巷道通风风流与粉尘分布规律深度剖析

基于数值模拟的综掘工作面巷道通风风流与粉尘分布规律深度剖析一、引言1.1研究背景与意义煤炭作为我国重要的能源资源，在能源结构中占据着关键地位。随着采煤机械化水平的持续提升，综掘工作面在煤炭开采过程中扮演着愈发重要的角色。然而，综掘工作面在掘进作业时会产生大量粉尘，这些粉尘不仅对煤矿安全生产构成严重威胁，还会对作业人员的身体健康造成极大危害。在煤矿生产中，粉尘是引发尘肺病的主要原因之一。尘肺病是一种不可逆的肺部疾病，会导致肺部组织纤维化，严重影响呼吸功能，降低作业人员的生活质量，甚至危及生命。据相关统计数据显示，尘肺病在我国职业病中所占比例高达80%以上，而煤矿行业又是尘肺病的高发行业。除了对人体健康的危害，粉尘还可能引发煤尘爆炸事故，给煤矿企业带来巨大的经济损失和人员伤亡。当空气中的煤尘浓度达到一定范围，遇到火源时，就可能发生爆炸，爆炸产生的高温、高压和冲击波会摧毁巷道、设备，造成严重的破坏。此外，高浓度的粉尘还会降低巷道内的能见度，影响作业人员的视线，增加误操作的风险，进而影响生产效率和安全性。在粉尘浓度较高的环境中，作业人员难以清晰地观察设备运行状态和周围环境，容易导致操作失误，引发设备故障或安全事故。为了有效降低综掘工作面粉尘浓度，保障煤矿安全生产和作业人员的身体健康，深入研究巷道内通风风流与粉尘分布规律至关重要。通风作为降低粉尘浓度的重要手段之一，其风流的运动状态直接影响着粉尘的扩散和分布。通过合理的通风设计，可以有效地稀释和排出粉尘，降低工作面粉尘浓度。因此，研究通风风流与粉尘分布规律，对于优化通风系统、提高通风效率、降低粉尘浓度具有重要的指导意义。传统的研究方法主要依赖于现场实测和物理模型试验。现场实测虽然能够获取真实的工况数据，但存在测量难度大、成本高、周期长等问题，且在某些复杂工况下难以进行全面测量。物理模型试验则受到模型制作、相似准则等因素的限制，难以准确模拟实际情况，且试验过程较为繁琐，成本也较高。随着计算机技术和计算流体力学（CFD）的飞速发展，数值模拟方法在煤矿通风与粉尘防治领域得到了广泛应用。数值模拟方法能够通过建立数学模型，对综掘工作面巷道内的通风风流和粉尘分布进行精确模拟，克服了传统研究方法的局限性。通过数值模拟，可以在计算机上快速、准确地分析不同通风条件、掘进工艺等因素对风流和粉尘分布的影响，为通风系统的优化设计和粉尘防治措施的制定提供科学依据。同时，数值模拟还可以减少现场试验的次数，降低研究成本，提高研究效率。因此，采用数值模拟方法研究综掘工作面巷道内通风风流与粉尘分布规律具有重要的现实意义和应用价值。1.2国内外研究现状在煤矿开采领域，综掘工作面巷道通风风流与粉尘分布规律一直是研究的重点和热点。国内外众多学者和研究机构运用多种方法，从不同角度对这一课题展开了深入研究，取得了一系列有价值的成果。在国外，早期的研究主要侧重于通风系统的设计与优化，以满足井下作业的基本通风需求。随着对粉尘危害认识的加深，研究逐渐转向通风与粉尘控制的协同作用。美国学者[具体姓名1]通过现场实测，分析了不同通风方式下巷道内风流速度和粉尘浓度的变化情况，发现压入式通风在稀释粉尘方面具有一定优势，但也容易导致粉尘在巷道内的扩散范围扩大。德国的研究团队[具体团队名称1]利用物理模型试验，研究了通风参数对粉尘分布的影响，提出了通过优化通风参数来降低工作面粉尘浓度的方法。他们发现，合理调整风筒出口位置和风速，可以有效改善风流场，减少粉尘的积聚。近年来，随着数值模拟技术的飞速发展，国外在该领域的研究取得了新的突破。英国学者[具体姓名2]运用CFD软件对综掘工作面巷道内的通风风流和粉尘运动进行了数值模拟，详细分析了不同工况下粉尘的扩散轨迹和浓度分布，为通风系统的优化设计提供了重要依据。他们的研究表明，数值模拟能够准确预测粉尘的分布情况，为实际生产中的粉尘防治提供了科学指导。澳大利亚的研究人员[具体姓名3]则将人工智能技术引入到通风与粉尘控制的研究中，通过建立智能模型，实现了对通风系统的自动调控，提高了粉尘防治的效果。他们利用机器学习算法，对大量的通风和粉尘数据进行分析，建立了通风参数与粉尘浓度之间的关系模型，实现了根据粉尘浓度实时调整通风参数的功能。在国内，相关研究起步相对较晚，但发展迅速。早期，国内主要借鉴国外的研究成果和经验，开展现场实测和物理模型试验研究。中国矿业大学的[具体姓名4]通过对多个煤矿综掘工作面的现场实测，总结了通风风流与粉尘分布的基本规律，提出了一些针对性的防尘措施，如增加通风量、合理布置通风设施等。这些措施在一定程度上降低了工作面粉尘浓度，但也存在一些局限性，如通风能耗增加、防尘效果不稳定等。随着计算机技术的普及，数值模拟方法在国内得到了广泛应用。山东科技大学的[具体姓名5]采用计算流体力学的离散相模型，对长压短抽式通风掘进巷道中的粉尘运动规律进行了数值模拟，得出了工作面附近区域粉尘质量浓度较高，然后沿程逐渐降低；距离压入式风筒另一侧煤壁越近，回风侧空间内粉尘质量浓度越高的分布规律。这一研究成果为优化通风系统和制定粉尘防治措施提供了重要参考。西安科技大学的[具体姓名6]利用流场模拟试验方法，分析获取出风口风流参数变化下粉尘浓度的回归数据，以掘进机司机位置处呼尘及回风侧全尘质量浓度最低为优化目标，构建风流调控下的粉尘浓度双目标优化模型，基于粒子群算法求解该模型，并选取理想解作为最优风流调控方案。通过实际应用，验证了该方案的有效性，为实现综掘面安全高效通风降尘提供了依据。尽管国内外在综掘工作面巷道通风风流与粉尘分布规律的研究方面取得了显著进展，但仍存在一些不足之处。部分研究在建立模型时，对实际工况的考虑不够全面，如未充分考虑巷道的复杂形状、设备的干扰以及粉尘的二次飞扬等因素，导致模拟结果与实际情况存在一定偏差。在多因素耦合作用下，通风风流与粉尘分布规律的研究还不够深入，例如通风方式、风速、粉尘粒径、湿度等因素之间的相互影响机制尚未完全明确。此外，目前的研究主要集中在常规条件下的通风与粉尘控制，对于特殊工况，如高瓦斯矿井、深部开采等情况下的通风风流与粉尘分布规律研究较少，难以满足复杂地质条件下煤矿安全生产的需求。综上所述，进一步深入研究综掘工作面巷道内通风风流与粉尘分布规律，完善数值模拟模型，加强多因素耦合作用和特殊工况下的研究，对于提高煤矿通风效率、降低粉尘浓度、保障安全生产具有重要意义。1.3研究内容与方法1.3.1研究内容本文旨在通过数值模拟的方法，深入研究综掘工作面巷道内通风风流与粉尘分布规律，具体研究内容如下：建立数值模型：依据某煤矿综掘工作面的实际地质条件、巷道尺寸、通风系统以及设备布置等参数，运用专业的建模软件，构建精准的三维巷道物理模型。充分考虑巷道的形状、坡度、转弯等因素，以及掘进机、刮板输送机、风筒等设备对风流和粉尘运动的影响，确保模型能够真实反映实际工况。在此基础上，根据计算流体力学和多相流理论，确定合适的数学模型和控制方程，如采用标准k-ε湍流模型来描述风流的湍流特性，利用离散相模型（DPM）来追踪粉尘颗粒的运动轨迹，同时考虑粉尘颗粒与风流之间的相互作用，如曳力、重力、浮力等，为后续的数值模拟提供可靠的理论基础。模拟通风风流分布规律：在建立的数值模型基础上，设定不同的通风工况，包括不同的通风方式（如压入式通风、抽出式通风、长压短抽式通风等）、通风量以及风筒出口位置和角度等参数。运用数值模拟软件对不同工况下巷道内的通风风流场进行模拟计算，得到风流速度、压力、流线等参数的分布情况。通过对模拟结果的分析，揭示通风风流在巷道内的流动特性和分布规律，研究不同通风参数对风流场的影响机制，如通风方式如何影响风流的有效射程和覆盖范围，通风量的变化如何改变风流速度和压力分布，以及风筒出口位置和角度对风流扩散和涡流形成的影响等，为优化通风系统设计提供理论依据。模拟粉尘分布规律：在模拟通风风流场的基础上，考虑粉尘的产生源、粒径分布、初始浓度等因素，对不同工况下巷道内的粉尘分布进行数值模拟。通过离散相模型追踪粉尘颗粒在风流中的运动轨迹，计算不同时刻、不同位置处的粉尘浓度，得到粉尘浓度场的分布情况。分析粉尘在巷道内的扩散、沉降、积聚等规律，研究通风参数、粉尘特性以及巷道结构等因素对粉尘分布的影响，如通风方式和通风量如何影响粉尘的扩散范围和浓度降低速度，粉尘粒径大小对其沉降速度和扩散距离的影响，以及巷道的弯道、交叉口等结构对粉尘积聚的影响等，为制定有效的粉尘防治措施提供科学依据。分析影响因素：综合考虑通风方式、通风量、粉尘粒径、湿度、巷道壁面粗糙度等多种因素，通过数值模拟和对比分析，深入研究各因素对通风风流与粉尘分布规律的单独影响和耦合作用。建立各因素与通风风流和粉尘分布参数之间的定量关系，如通风量与风流速度、粉尘浓度之间的函数关系，粉尘粒径与沉降速度、扩散距离之间的数学模型等。通过敏感性分析，确定影响通风风流与粉尘分布的关键因素，为针对性地优化通风系统和粉尘防治措施提供理论支持，明确在实际生产中应重点控制哪些因素，以达到最佳的通风降尘效果。优化通风与降尘方案：根据模拟结果和影响因素分析，以降低巷道内粉尘浓度、提高通风效率为目标，提出优化的通风与降尘方案。对通风系统的参数进行优化调整，如确定最佳的通风方式、通风量、风筒出口位置和角度等，以改善风流场分布，增强对粉尘的稀释和排出能力。结合其他降尘措施，如喷雾降尘、除尘器除尘等，研究不同降尘措施与通风系统的协同作用，提出综合的降尘方案。通过数值模拟对优化后的方案进行验证和评估，对比优化前后巷道内通风风流和粉尘分布的变化情况，分析优化方案的降尘效果和经济效益，确保优化方案的可行性和有效性，为煤矿综掘工作面的实际生产提供技术指导。1.3.2研究方法数值模拟软件：选用计算流体力学（CFD）软件FLUENT作为主要的数值模拟工具。FLUENT具有强大的物理模型库和求解器，能够精确模拟各种复杂的流体流动和传热传质现象，在航空航天、汽车工程、能源等领域得到了广泛应用。在煤矿通风与粉尘防治领域，FLUENT也被众多学者用于研究通风风流和粉尘的运动规律。它可以方便地处理各种复杂的几何模型，提供多种湍流模型、多相流模型和离散相模型，能够满足本文对综掘工作面巷道内通风风流与粉尘分布规律的研究需求。通过FLUENT软件，能够对建立的巷道物理模型进行网格划分、边界条件设定、求解控制和结果后处理等操作，得到准确的模拟结果。理论基础：以计算流体力学（CFD）理论为核心，结合多相流理论和粉尘运动学理论，对综掘工作面巷道内的通风风流和粉尘分布进行数值模拟。计算流体力学是通过计算机数值计算和图像显示，对包含有流体流动和热传导等相关物理现象的系统所做的分析，它基于质量守恒方程、动量守恒方程和能量守恒方程，通过数值方法求解这些方程，得到流场的各种参数分布。在本文中，运用CFD理论来描述巷道内风流的运动特性，求解风流的速度、压力、温度等参数。多相流理论则用于处理风流与粉尘颗粒的相互作用，考虑了颗粒相在连续相（风流）中的运动、扩散和相互作用力。粉尘运动学理论主要研究粉尘颗粒在风流中的受力情况，如曳力、重力、浮力、Saffman升力等，通过建立颗粒的运动方程，追踪粉尘颗粒的运动轨迹，从而得到粉尘的分布规律。这些理论的综合应用，为数值模拟提供了坚实的理论基础，确保模拟结果的准确性和可靠性。二、综掘工作面巷道通风与粉尘的相关理论2.1综掘工作面巷道通风理论基础2.1.1通风的目的与任务综掘工作面巷道通风是保障煤矿安全生产和作业人员身体健康的关键环节，其目的与任务具有重要意义。通风的首要任务是为掘进工作面的工作人员提供足够的新鲜空气。在井下狭窄且封闭的空间中，人员的呼吸会不断消耗氧气并产生二氧化碳。正常情况下，人体每分钟需要吸入约0.5-1升的氧气，以维持正常的生理代谢。若通风不畅，氧气含量会迅速下降，二氧化碳浓度则会急剧上升。当氧气浓度低于18%时，人体会出现头晕、乏力、呼吸急促等症状；当氧气浓度低于12%时，可能会导致昏迷甚至死亡。因此，通过有效的通风系统，持续向工作面输送新鲜空气，确保氧气浓度保持在20%左右，是保障人员生命安全的基本要求。冲淡、排除掘进工作面的有害气体也是通风的重要任务之一。在综掘过程中，会产生多种有害气体，如一氧化碳（CO）、硫化氢（H₂S）、二氧化氮（NO₂）等。一氧化碳是一种无色、无味、无臭的剧毒气体，它与人体血红蛋白的亲和力比氧气大200-300倍，一旦吸入，会迅速与血红蛋白结合，形成碳氧血红蛋白，使血红蛋白失去携氧能力，导致人体组织缺氧。当空气中一氧化碳浓度达到0.016%时，数小时后会出现轻微头痛；当浓度达到0.128%时，1-2小时内会使人感到剧烈头痛、眩晕、呕吐，甚至昏迷；当浓度达到1.28%时，短时间内即可致人死亡。硫化氢具有臭鸡蛋气味，是一种剧毒气体，它对人体的中枢神经系统、呼吸系统和眼部等都有强烈的刺激和损害作用。二氧化氮是一种棕红色、有刺激性气味的气体，它会对呼吸道黏膜产生强烈的刺激，引发咳嗽、气喘、肺水肿等症状。这些有害气体若不能及时排出，会在巷道内积聚，对作业人员的生命安全构成严重威胁。通风系统通过引入新鲜空气，将有害气体稀释并排出到地面，使巷道内有害气体浓度低于安全标准，保障作业环境的安全。此外，通风还承担着创造良好工作环境的重要职责。在掘进过程中，机械设备的运转会产生大量的热量，同时，煤岩体的氧化、人员的活动等也会释放热量。若热量不能及时散发，会导致巷道内温度升高。当温度超过30℃时，人体会感到不适，工作效率明显下降；当温度超过35℃时，可能会引发中暑等疾病，严重影响作业人员的身体健康。通风系统通过空气的流动，将热量带走，使巷道内温度保持在适宜的范围内，一般要求井下作业地点的空气温度不得超过26℃。同时，通风还可以调节巷道内的湿度，避免因湿度过高或过低对设备和人员造成不利影响。适宜的湿度范围一般在40%-65%之间，湿度过高会使设备容易生锈，影响其使用寿命，还会使人感到闷热不适；湿度过低则会导致呼吸道黏膜干燥，引发咳嗽、咽干等症状。通过通风系统的调节，为作业人员创造一个温度适宜、湿度适中的工作环境，有助于提高工作效率和保障人员的身体健康。2.1.2常见通风方式及特点在综掘工作面巷道通风中，常见的通风方式主要有压入式通风、抽出式通风和压抽混合式通风，它们各自具有独特的工作原理、适用场景及优缺点。压入式通风：压入式通风是利用通风机产生的压力，将新鲜空气经风筒压入掘进工作面，使工作面的污浊空气沿巷道排出。其工作原理是通风机将外界的新鲜空气吸入，通过刚性或柔性风筒，将空气以一定的速度和压力输送到掘进工作面的迎头。新鲜空气在工作面与污浊空气混合后，沿着巷道向回风方向流动，最终排出巷道。这种通风方式在高瓦斯浓度的煤巷道和半煤岩巷中应用较为广泛。其优点显著，安置便捷，通风设备相对简单，成本低廉，只需将通风机和风筒布置在合适位置即可投入使用。而且有效射程远，风流从风筒末端射向工作面，有效射程一般可达7-8m，能够快速排出工作面的乏风和粉尘，通风效果较好。同时，局部通风机安设在新鲜风流中，安全性能较高，降低了因通风设备引发瓦斯爆炸等事故的风险。然而，压入式通风也存在一些缺点。在掘进过程中，综掘面需风量大，而压入式通风容易使瓦斯等有害气体在巷道内的扩散距离大、范围广。当巷道内采用压入式通风方式时，在送风筒射流有效射程之外的空间会出现循环涡流区，这会导致部分污浊空气在该区域积聚，难以排出。为了尽可能多排出巷道内污风和有害气体，送风筒出口与迎头面距离应在有效射程之内，否则会影响通风效果。抽出式通风：抽出式通风是利用抽风机产生的负压，将掘进工作面的含尘风流经风筒抽出，使新鲜空气从巷道进入工作面。抽风机工作时，在风筒内形成负压，将工作面的含尘风流吸入风筒，然后排出到巷道外。这种通风方式通常适用于低瓦斯矿井或对空气质量要求较高的掘进工作面。其优点在于能保持巷道内空气新鲜，因为含尘风流是通过风筒抽出，减少了对巷道内空气的污染，为人员创造了较好的作业环境。然而，抽出式通风也存在一些局限性。抽风筒一般为刚性材质，如铁皮等，成本较高，且安装和维护相对困难。此外，倘若吸走的污风中含有瓦斯和粉尘，当其通过发生故障的通风机时，容易发生爆炸事故，存在较大的安全隐患，故在高瓦斯井巷中不宜使用抽出式通风。压抽混合式通风：压抽混合式通风综合了压入式通风和抽出式通风的特点，在巷道内同时布置压风筒和抽风筒。压风筒为综掘巷道供风，将新鲜空气压入工作面；抽风筒在回风侧吸走含尘风流，将污浊空气排出。混合式通风可分为长抽短压、长压短抽和长抽长压三种，其中长压短抽通风方式使用最为广泛。长压短抽通风方式中，通常以压入式通风为主、抽出式通风为辅，且压入风量大于抽吸风量。这种通风方式的优点明显，相较于单独使用压入式或抽出式通风，它在保证综掘面供风的同时，还能缩短排出含尘风流的时间，缩小了粉尘在巷道内的扩散范围。在实际应用中，长压短抽通风方式中的压风筒采用柔性材质制成，成本较低，安装方便；抽出式一般采用铁皮等硬质材料制成，以保证其强度和密封性。然而，该种通风方式管理较为困难，需要合理调整压入风量和抽出风量的比例，以及风筒的布置位置和长度，确保通风系统的稳定运行。2.2综掘工作面巷道粉尘产生及危害2.2.1粉尘产生的主要工序与来源在综掘工作面巷道的掘进过程中，多种工序会产生大量粉尘，这些粉尘的来源广泛且复杂。掘进机切割煤岩是最主要的产尘工序。掘进机截割头在高速运转过程中，截割齿与煤岩壁面相接触，对煤岩进行强力破碎。在这个过程中，煤岩受到截齿的挤压、冲击和摩擦作用，其内部结构被破坏，原本紧密结合的颗粒被分离，从而产生大量粉尘。由于切割过程的高速性和冲击性，这些粉尘会以较高的速度向外飞溅，在掘进机附近形成高浓度的粉尘区域。据相关研究和实际测量，在未采取任何除尘措施时，掘进机切割区域附近的粉尘浓度可高达4800mg/m³以上，该工序产生的粉尘占巷道内粉尘总量的比例可达80%以上。煤岩塌落也是重要的产尘来源之一。在掘进作业时，部分煤岩块由于受到掘进机的震动、自身重力以及地质条件变化等因素的影响，会从巷道顶部或侧壁塌落至巷道底板。煤岩块在塌落过程中与掘进机外壳、巷道底板发生碰撞，以及铲板、转载机对煤块的进一步摩擦破碎作用，都会使煤岩块表面的细小颗粒脱落，从而产生粉尘。此外，在煤岩装载过程中，铲斗与物料的接触、矿车与运输设备的碰撞等，也会因摩擦和振动产生一定量的粉尘。地质构造的复杂性对粉尘产生也有显著影响。在煤层形成过程中，由于地质作用，煤岩层会发生断裂和错位现象，在这些错位煤层的空隙中会留存大量粉尘。当掘进机进行作业时，原本封闭的煤层被打开，这些空隙中的粉尘暴露于外界，在掘进机的扰动和风流的作用下，迅速扩散到巷道空间中，成为巷道粉尘的重要组成部分。在地质构造复杂的区域，如断层、褶皱等附近，掘进开采时的产尘量通常会明显增大。二次扬尘也是不可忽视的粉尘来源。在巷道掘进过程中，已经沉降在巷道底板和掘进机外壳上的粉尘，会受到风流的作用而再次扬起，重新进入空气中，形成二次扬尘。当通风风速过大时，风流对沉降粉尘的扰动作用增强，更容易引发二次扬尘，从而进一步提高巷道内的粉尘浓度。巷道内的其他作业活动，如人员走动、设备移动等，也可能导致沉降粉尘的再次飞扬。2.2.2粉尘对人体健康和安全生产的危害粉尘对人体健康和煤矿安全生产都有着严重的危害，其影响不容忽视。在人体健康方面，长期吸入粉尘会对呼吸系统造成极大的损害，尘肺病是最为典型的职业病。尘肺病是由于在职业活动中长期吸入生产性粉尘，并在肺内潴留而引起的以肺组织弥漫性纤维化为主的全身性疾病。当人体吸入的粉尘进入肺部后，巨噬细胞会试图吞噬这些粉尘颗粒，但一些难以被消化的粉尘，如含有游离二氧化硅的粉尘，会导致巨噬细胞死亡。随着时间的推移，大量的粉尘在肺部积聚，引发肺部组织的炎症反应，促使成纤维细胞增生，逐渐形成胶原纤维结节，导致肺组织纤维化。尘肺病一旦发生，是不可逆的，患者会逐渐出现咳嗽、咳痰、胸痛、呼吸困难等症状，肺功能不断下降，严重影响生活质量，甚至危及生命。据统计，尘肺病在我国职业病中所占比例高达80%以上，煤矿行业是尘肺病的高发行业，许多煤矿工人因长期暴露在高浓度粉尘环境中而患病。粉尘对安全生产也构成严重威胁，煤尘爆炸是最为严重的安全事故之一。煤尘具有可燃性，当空气中的煤尘浓度达到一定范围，一般为45-2000g/m³，遇到火源时，就可能发生爆炸。煤尘爆炸产生的高温可达2000℃以上，高压可达7-8个大气压，爆炸产生的冲击波会以极高的速度传播，摧毁巷道、设备，造成严重的人员伤亡和财产损失。煤尘爆炸还可能引发连锁反应，导致瓦斯爆炸等更严重的事故，进一步扩大灾害范围。在煤矿生产历史上，发生过多起因煤尘爆炸导致的重大事故，给煤矿企业和社会带来了巨大的伤痛和损失。此外，高浓度的粉尘还会降低巷道内的能见度，严重影响作业人员的视线。在粉尘浓度较高的环境中，光线被大量散射和吸收，作业人员难以清晰地观察设备运行状态和周围环境，增加了误操作的风险。这不仅容易导致设备故障，影响生产效率，还可能引发其他安全事故，如人员碰撞、滑倒等。粉尘还会对设备造成磨损和腐蚀，缩短设备的使用寿命，增加设备维护成本，进一步影响煤矿生产的正常进行。2.3数值模拟相关理论与方法2.3.1气固两相流理论气固两相流理论是研究气体与固体颗粒混合流动现象的重要理论，在描述综掘工作面巷道内通风风流与粉尘运动中发挥着关键作用。在综掘工作面巷道中，通风风流可视为连续相，而粉尘则是离散分布于风流中的分散相，二者共同构成了气固两相流体系。当风流在巷道中流动时，会与粉尘颗粒发生相互作用。一方面，风流为粉尘颗粒提供了运动的动力，通过曳力的作用，使粉尘颗粒跟随风流一起运动。曳力是气固两相流中最为重要的作用力之一，其大小与风流速度、粉尘颗粒的形状和尺寸、以及流体的粘性等因素密切相关。根据斯托克斯定律，对于球形颗粒在层流状态下，曳力与颗粒直径的平方、风流速度以及流体粘性成正比。另一方面，粉尘颗粒也会对风流的流动特性产生影响。大量的粉尘颗粒会增加流体的有效密度和粘度，改变风流的速度分布和压力场。当粉尘浓度较高时，颗粒之间的相互碰撞和摩擦会消耗风流的能量，导致风流速度降低，压力损失增大。此外，在气固两相流中，还需要考虑其他多种作用力对粉尘颗粒运动的影响。重力是不可忽视的作用力，它使粉尘颗粒有向下沉降的趋势。对于粒径较大、密度较高的粉尘颗粒，重力的作用更为明显。在巷道中，若通风风速较小，不足以克服重力的影响，粉尘颗粒就会逐渐沉降到巷道底部。浮力则是由于颗粒与周围流体密度不同而产生的作用力，当粉尘颗粒的密度小于风流的密度时，浮力会使颗粒向上运动。Saffman升力是在非均匀流场中，由于速度梯度的存在而产生的对颗粒的作用力，它会使颗粒在垂直于速度梯度的方向上发生偏移。此外，当颗粒带有电荷时，还会受到静电力的作用，静电力会影响颗粒之间的相互作用以及颗粒与巷道壁面的粘附和分离。在实际的综掘工作面巷道中，气固两相流的运动非常复杂，受到多种因素的综合影响。巷道的几何形状、通风方式、通风量、粉尘的产生源和粒径分布等因素都会对气固两相流的特性产生显著影响。在弯曲的巷道中，风流会产生二次流，从而改变粉尘颗粒的运动轨迹；不同的通风方式会导致风流的速度场和压力场分布不同，进而影响粉尘的扩散和沉降。因此，在研究综掘工作面巷道内通风风流与粉尘分布规律时，需要综合考虑这些因素，运用气固两相流理论建立准确的数学模型，以深入理解和预测其运动行为。2.3.2计算流体力学（CFD）方法计算流体力学（CFD）方法是一种基于计算机数值计算和图像显示，对包含流体流动和热传导等相关物理现象的系统进行分析的数值模拟技术，在模拟巷道内复杂流场和粉尘扩散方面具有独特的优势。CFD方法的核心是基于质量守恒方程、动量守恒方程和能量守恒方程这三大基本守恒定律。质量守恒方程，也称为连续性方程，它表明在一个封闭的控制体内，流体的质量不会凭空产生或消失，流入控制体的质量流量等于流出控制体的质量流量与控制体内质量变化率之和。对于不可压缩流体，其密度为常数，连续性方程可简化为速度散度为零。在综掘工作面巷道通风中，通过连续性方程可以确保风流在巷道内的流动满足质量守恒，准确计算风流在不同位置的流量和速度分布。动量守恒方程，即纳维-斯托克斯方程（N-S方程），它描述了流体动量的变化与作用在流体上的力之间的关系。这些力包括压力梯度力、粘性力、重力等。在巷道内，风流受到通风机提供的压力作用，克服巷道壁面的摩擦阻力以及重力等其他作用力，从而实现流动。通过求解N-S方程，可以得到风流的速度、压力等参数的分布情况，进而分析风流在巷道内的流动特性，如是否存在涡流、流速的大小和方向变化等。能量守恒方程则考虑了流体的内能、动能和势能的变化，以及热量的传递和做功等因素。在巷道通风中，虽然热量传递和做功的影响相对较小，但在一些特殊情况下，如高温矿井或通风系统中存在热交换设备时，能量守恒方程对于准确描述风流的状态变化至关重要。在实际应用CFD方法模拟巷道内通风风流和粉尘扩散时，首先需要对巷道进行几何建模，将实际的巷道形状和尺寸转化为计算机可处理的几何模型。然后，对几何模型进行网格划分，将巷道空间离散为众多的小单元，这些小单元被称为网格。网格的质量和数量对模拟结果的准确性和计算效率有很大影响。高质量的网格应具有合理的形状和尺寸分布，能够准确地描述巷道的几何特征和流场的变化。一般来说，在流场变化剧烈的区域，如通风机出口、巷道弯道和交叉口等位置，需要加密网格，以提高计算精度；而在流场变化相对平缓的区域，可以适当放宽网格尺寸，以减少计算量。接着，根据实际情况设置边界条件，包括入口边界条件、出口边界条件和壁面边界条件等。入口边界条件通常指定风流的速度、温度、压力等参数，以及粉尘的浓度和粒径分布等信息；出口边界条件则根据巷道的实际情况，选择合适的压力或速度条件；壁面边界条件主要考虑巷道壁面对风流的粘性作用，一般采用无滑移边界条件，即认为壁面处的流体速度为零。在完成几何建模、网格划分和边界条件设置后，选择合适的数值求解方法对控制方程进行求解。常见的数值求解方法有有限差分法、有限体积法和有限元法等。有限体积法由于其在守恒性和计算效率方面的优势，在CFD模拟中得到了广泛应用。通过数值求解，可以得到巷道内不同位置处风流的速度、压力、温度等参数的分布，以及粉尘颗粒的运动轨迹和浓度分布。最后，对模拟结果进行后处理和分析。利用CFD软件提供的后处理工具，可以将模拟结果以图形化的方式展示出来，如速度矢量图、压力云图、流线图和粉尘浓度分布云图等。这些直观的图形有助于研究人员深入理解巷道内通风风流和粉尘的分布规律，分析不同因素对其的影响，从而为通风系统的优化设计和粉尘防治措施的制定提供科学依据。2.3.3常用数值模拟软件介绍（如FLUENT）FLUENT是一款功能强大的计算流体力学（CFD）软件，在处理通风与粉尘问题时展现出诸多显著的功能和优势，成为众多科研人员和工程师在该领域进行数值模拟研究的首选工具。FLUENT具备丰富的物理模型库，涵盖了多种适用于通风与粉尘模拟的模型。在处理通风风流问题时，它提供了多种湍流模型，如标准k-ε模型、RNGk-ε模型、Realizablek-ε模型以及k-ω系列模型等。标准k-ε模型是应用最为广泛的湍流模型之一，它基于湍动能k和湍流耗散率ε这两个变量来描述湍流特性，能够较好地模拟一般工程中的湍流流动，在巷道通风风流模拟中，对于大部分工况都能给出较为准确的结果。RNGk-ε模型则在标准k-ε模型的基础上，考虑了湍流的旋流效应和流线弯曲等因素，对于一些复杂的流场，如存在涡流和弯曲巷道的情况，具有更好的模拟精度。Realizablek-ε模型引入了新的湍流粘性公式和耗散率方程，能够更准确地预测一些具有强旋流和分离流动的流场，在处理通风风流中的复杂流动现象时具有独特的优势。k-ω系列模型则在近壁区域的模拟中表现出色，对于研究巷道壁面附近的风流特性具有重要意义。在处理粉尘问题方面，FLUENT提供了离散相模型（DPM）来追踪粉尘颗粒的运动轨迹。DPM模型基于拉格朗日坐标系，通过求解单个颗粒的运动方程，考虑颗粒受到的各种作用力，如曳力、重力、浮力、Saffman升力等，来计算颗粒在流场中的运动轨迹和速度变化。同时，该模型还可以考虑颗粒与风流之间的相互作用，以及颗粒之间的碰撞和聚并等现象。在模拟粉尘扩散时，DPM模型能够准确地预测粉尘在巷道内的分布情况，为研究粉尘的传播规律和制定有效的降尘措施提供了有力的支持。此外，FLUENT还提供了多相流模型，如欧拉-拉格朗日模型和欧拉-欧拉模型等，这些模型可以更全面地考虑风流与粉尘颗粒之间的相互作用，以及颗粒相内部的相互作用，适用于处理一些更为复杂的气固两相流问题。FLUENT的求解器具有高效性和稳定性，能够快速准确地求解复杂的流体力学方程。它采用了先进的数值算法，如有限体积法、SIMPLE算法及其改进算法等，能够有效地处理各种复杂的边界条件和非线性问题。在模拟综掘工作面巷道内通风风流和粉尘分布时，即使面对复杂的巷道几何形状和多物理场耦合问题，FLUENT的求解器也能够稳定收敛，得到准确的计算结果。同时，FLUENT还支持并行计算，通过利用多核心处理器或集群计算资源，可以显著提高计算效率，缩短模拟时间，使得大规模的数值模拟研究成为可能。在模型建立和网格划分方面，FLUENT拥有强大的前处理功能。它可以直接导入多种常见的三维建模软件生成的几何模型，如ANSYSDesignModeler、SolidWorks、Pro/E等，方便用户根据实际巷道情况进行精确建模。同时，FLUENT自身也提供了丰富的几何建模工具，能够满足一些简单模型的创建需求。在网格划分方面，FLUENT支持结构化网格、非结构化网格以及混合网格等多种网格类型。对于复杂的巷道模型，非结构化网格能够更好地适应几何形状的变化，提高网格质量和计算精度；而结构化网格则在一些规则形状的区域具有计算效率高的优势，混合网格则结合了两者的优点。FLUENT还提供了多种网格生成方法和网格加密技术，用户可以根据流场的特点和计算精度的要求，灵活地进行网格划分和优化，以获得最佳的模拟效果。此外，FLUENT还具备强大的后处理功能，能够以多种方式展示模拟结果。它可以生成各种直观的图形，如速度矢量图、压力云图、流线图、粉尘浓度分布云图等，帮助用户清晰地了解通风风流和粉尘在巷道内的分布情况。同时，FLUENT还支持数据的输出和处理，用户可以将模拟结果导出为各种格式的数据文件，以便进行进一步的数据分析和处理。例如，通过对不同工况下的模拟结果进行对比分析，可以研究通风参数、巷道结构等因素对通风风流和粉尘分布的影响规律，为通风系统的优化设计和粉尘防治措施的制定提供科学依据。三、数值模拟模型的建立3.1模型的假设与简化为了在合理范围内简化数值模拟过程，同时确保模拟结果能有效反映综掘工作面巷道内通风风流与粉尘分布的关键特征，对巷道、设备、风流和粉尘等进行了以下合理的假设与简化处理。在巷道模型方面，实际的综掘工作面巷道往往存在一定的起伏、不规则形状以及复杂的地质构造，但为了便于建模和计算，将巷道简化为规则的几何形状。例如，通常将巷道视为直墙半圆拱形或矩形，忽略巷道壁面的微小凹凸不平和局部地质变化对风流和粉尘运动的细微影响。尽管实际巷道可能存在一些拐弯和分支，但在本次模拟中，先重点研究直巷道的情况，将巷道视为单一的、连续的通道，不考虑巷道分支和复杂连接结构对风流和粉尘分布的影响，以突出主要的通风和粉尘扩散规律。同时，假设巷道壁面为刚性、光滑表面，不考虑壁面的粗糙度和孔隙度对风流的摩擦阻力和吸附作用，虽然实际壁面存在一定粗糙度会增加风流阻力并影响粉尘的附着和再扬起，但在初步模拟中简化这一因素，可使问题更容易求解，后续再通过添加修正系数等方式考虑壁面粗糙度的影响。对于设备模型，掘进机、刮板输送机、风筒等设备在实际工作中形状和结构都较为复杂。为简化模型，将掘进机简化为规则的长方体，忽略其截割头、悬臂、履带等复杂部件的细节，仅保留其主要的几何尺寸和占据空间的位置信息。刮板输送机则简化为扁平的长方体，不考虑刮板的具体形状和运动对物料和风流的影响。风筒假设为理想的圆形管道，忽略风筒的接头、褶皱以及可能存在的漏风情况，将风筒壁视为光滑的刚性壁面，这样可以方便地确定风筒内风流的边界条件和流动特性。虽然这些设备在实际工作中会对风流和粉尘产生复杂的扰动和阻挡作用，但通过简化模型可以先获得基本的通风和粉尘分布规律，后续再逐步细化模型以考虑设备的详细结构和工作状态对风流和粉尘的影响。在风流模型方面，假设风流为不可压缩的牛顿流体，即风流的密度和粘性系数在模拟过程中保持不变。尽管实际风流在井下复杂环境中可能会受到温度、压力变化以及气体成分变化的影响，导致其密度和粘性发生改变，但在通常的通风条件下，风流的压缩性和非牛顿特性对通风风流与粉尘分布规律的影响相对较小，通过这一假设可以大大简化计算过程。同时，忽略风流中的热传递和质量扩散现象，不考虑风流与巷道壁面、设备以及粉尘之间的热量交换和物质交换，仅关注风流的流动特性和对粉尘的携带作用，这样可以将问题聚焦于通风风流与粉尘的气固两相流运动，简化模型的复杂性。在粉尘模型方面，考虑到粉尘颗粒的形状和性质差异较大，为简化计算，将粉尘颗粒视为球形，且忽略粉尘颗粒之间的相互碰撞和聚并作用。尽管实际粉尘颗粒形状不规则，且在高浓度情况下颗粒之间的相互作用较为明显，但在初步模拟中，将粉尘颗粒理想化处理，可以先获得粉尘在风流中的基本运动规律和分布特征。同时，假设粉尘颗粒的物理性质，如密度、粒径等，在整个模拟区域内均匀分布，不考虑粉尘颗粒在产生源处的初始速度和方向差异，以及由于不同工序产生的粉尘特性差异，这样可以减少模型的参数数量，降低计算难度，后续可通过进一步细化模型来考虑这些因素对粉尘分布的影响。3.2几何模型的构建以某煤矿综掘工作面巷道为研究对象，该巷道位于煤层中，为满足矿井生产和运输需求而设计。巷道设计长度为1000m，在本次模拟中，选取其中具有代表性的200m一段进行建模分析。巷道采用直墙半圆拱形断面，这种断面形状在煤矿巷道中应用广泛，具有较好的稳定性和承载能力。其净宽为4.5m，净高为3.8m，其中直墙部分高度为1.8m，半圆拱半径为1.35m。运用专业的三维建模软件ANSYSDesignModeler进行模型构建。首先，创建一个长度为200m的长方体，其宽度和高度分别对应巷道的净宽4.5m和净高3.8m，以此作为巷道的基本轮廓。然后，利用软件的布尔运算功能，对长方体的顶部进行处理，通过拉伸和切割操作，构建出半径为1.35m的半圆拱形，使其与直墙部分组合形成完整的直墙半圆拱形巷道断面。在巷道模型中，掘进机是重要的设备之一。该掘进机型号为EBZ200，其外形近似长方体，长度为12m，宽度为3.2m，高度为2.5m。在建模软件中，通过创建相应尺寸的长方体来模拟掘进机，并将其放置在巷道内距离迎头20m的位置，该位置是根据实际掘进作业中的常见工况确定的，能够较好地反映掘进机工作时对通风风流和粉尘分布的影响。刮板输送机用于运输掘进过程中产生的煤岩，其沿着巷道中轴线布置，贯穿整个模拟巷道。刮板输送机可简化为一个长度为200m，宽度为0.8m，高度为0.3m的长方体。在建模时，将其放置在巷道底部中央位置，确保其与巷道和其他设备的相对位置关系准确，以模拟其在实际工作中的运行状态和对风流、粉尘的影响。风筒是通风系统的关键部件，负责将新鲜空气输送到掘进工作面。该风筒采用直径为0.8m的柔性风筒，悬挂在巷道顶板一侧，距离巷道左帮0.5m，距离顶板0.3m。在建模软件中，通过创建圆柱体来模拟风筒，并根据实际悬挂位置进行定位。风筒入口位于巷道开口处，出口则靠近掘进工作面迎头，距离迎头5m，这样的布置方式符合该煤矿综掘工作面的实际通风设计，能够准确模拟新鲜空气的输送路径和风流在巷道内的分布情况。在构建几何模型的过程中，严格按照实际尺寸和相对位置进行建模，确保模型的准确性和真实性。同时，对模型进行了仔细的检查和修正，避免出现几何错误和不合理的结构，为后续的网格划分和数值模拟提供可靠的基础。通过精确的几何模型构建，能够更准确地模拟综掘工作面巷道内通风风流与粉尘的分布规律，为研究和解决实际工程问题提供有力的支持。3.3网格划分将在ANSYSDesignModeler中构建好的综掘工作面巷道几何模型导入到专业的网格划分软件ICEMCFD中，进行网格划分操作。在网格划分过程中，选择合适的网格类型和划分方法是确保计算精度和效率的关键。考虑到巷道模型的复杂几何形状，尤其是直墙半圆拱形的巷道断面以及其中布置的掘进机、刮板输送机、风筒等设备，采用非结构化四面体网格进行划分。非结构化四面体网格能够更好地适应复杂的几何边界，对模型的细节部分进行准确的网格离散，从而提高网格质量和计算精度。在确定网格尺寸时，需要综合考虑计算精度和计算资源的平衡。对于巷道整体区域，初始设置全局网格尺寸为0.3m。这样的尺寸既能保证对巷道整体结构和主要风流流动特征的准确模拟，又不会导致计算量过大。然而，在一些对计算精度要求较高的关键区域，如掘进机周围、风筒出口附近以及巷道壁面附近等，需要对网格进行局部加密处理。掘进机在工作过程中会对风流和粉尘的运动产生强烈的扰动，其周围的流场变化复杂，因此将掘进机周围区域的网格尺寸加密至0.1m，以更精确地捕捉风流在该区域的速度变化、涡流形成以及粉尘颗粒与掘进机表面的相互作用。风筒出口是风流的主要入口，风流在此处的速度和方向变化较大，对整个巷道内的风流场分布起着关键作用，所以将风筒出口附近区域的网格尺寸也加密至0.1m，以准确模拟风流从风筒射出后的扩散和衰减过程。巷道壁面附近存在边界层效应，风流速度在壁面处会急剧降低，为了准确描述这一现象，对巷道壁面附近一定厚度范围内的网格进行加密，网格尺寸同样设置为0.1m。在网格划分完成后，对生成的网格质量进行严格检查。通过ICEMCFD软件提供的网格质量检查工具，查看网格的各项质量指标，如最小内角、纵横比、雅克比行列式等。确保网格的最小内角大于15°，以避免出现过于尖锐的网格单元，防止在计算过程中产生数值不稳定的情况。纵横比尽量控制在10以下，保证网格单元的形状较为规则，提高计算精度。雅克比行列式的值在合理范围内，一般要求大于0.2，以确保网格的变形不会过大，保证计算结果的可靠性。经过检查和调整，最终生成的网格质量良好，能够满足数值模拟的要求。整个模型共划分了约[X]万个网格单元，这些网格单元能够准确地离散巷道模型的几何形状和物理场，为后续在FLUENT软件中的数值计算提供了坚实的基础，确保能够准确地模拟综掘工作面巷道内通风风流与粉尘的分布规律。3.4边界条件的设定在FLUENT软件中对模型进行边界条件的设定，以准确模拟综掘工作面巷道内通风风流与粉尘的实际运动情况。对于进风口，即风筒入口，采用速度入口边界条件。根据该煤矿综掘工作面的实际通风设计，风筒出口风速设定为15m/s。这一风速值是基于矿井通风系统的风量需求以及风筒的规格和阻力计算得出的，能够保证向掘进工作面提供足够的新鲜空气，满足作业人员的呼吸需求并有效稀释有害气体和粉尘。同时，设置湍流强度为5%，水力直径为0.8m，以准确描述风流进入风筒时的湍流特性。湍流强度表示风流中脉动速度的大小，5%的湍流强度符合实际通风系统中风流的湍流情况，能够合理地模拟风流在风筒内的流动和进入巷道后的扩散过程。水力直径则是用于表征风筒内流体流动通道的特征尺寸，对于圆形风筒，其水力直径等于风筒的直径，即0.8m。出风口采用压力出口边界条件，表压力设置为0Pa。这是因为出风口与大气相通，其压力与大气压力相等，在相对压力的计算中，表压力为0。同时，设置回流湍流强度为5%，回流湍流粘度比为10。回流湍流强度和回流湍流粘度比用于描述出风口处回流现象的湍流特性，5%的回流湍流强度和10的回流湍流粘度比能够较好地反映实际情况下出风口处风流的回流和混合情况，确保模拟结果的准确性。壁面边界条件方面，巷道壁面、掘进机表面、刮板输送机表面和风筒壁面均采用无滑移壁面边界条件。这意味着在这些壁面处，风流的速度为零，符合实际情况中流体与固体壁面之间的粘附作用。同时，考虑到壁面对风流的摩擦阻力，设置壁面粗糙度高度为0.001m，粗糙度常数为0.5。壁面粗糙度高度和粗糙度常数用于描述壁面的粗糙程度对风流的影响，0.001m的粗糙度高度和0.5的粗糙度常数是根据实际巷道壁面和设备表面的粗糙情况确定的，能够合理地模拟壁面对风流的摩擦阻力，使模拟结果更接近实际情况。在粉尘源的设置上，考虑到掘进机截割头在工作过程中是主要的粉尘产生源，在截割头表面设置面源。根据实际测量和相关研究，确定粉尘的初始质量流量为0.1kg/s，这一数值是基于掘进机的工作参数、煤岩的性质以及现场的粉尘浓度测量数据综合确定的，能够准确反映截割头在单位时间内产生的粉尘量。同时，设置粉尘颗粒的直径范围为1-100μm，平均粒径为10μm。粉尘颗粒的粒径分布对其在风流中的运动和扩散规律有重要影响，1-100μm的粒径范围涵盖了实际生产中常见的粉尘粒径大小，平均粒径为10μm则是为了在模拟中能够代表大部分粉尘颗粒的粒径特征，以便更准确地模拟粉尘的运动轨迹和浓度分布。此外，粉尘颗粒的密度设置为1350kg/m³，这是根据所开采煤岩的实际密度确定的，能够合理地反映粉尘颗粒在重力作用下的沉降特性。通过合理设定以上边界条件，能够使数值模拟模型更加真实地反映综掘工作面巷道内通风风流与粉尘的实际运动和分布情况，为后续的模拟分析提供可靠的基础。3.5求解器与计算参数的选择在FLUENT软件中，选择分离式求解器进行数值计算。分离式求解器适用于求解低速、不可压缩或弱可压缩流动问题，通过将压力和速度的耦合求解过程分离，采用迭代的方式逐步逼近真实解。在每一次迭代中，先求解动量方程得到速度场，然后根据速度场求解压力修正方程，以修正压力和速度，使得连续性方程得到满足。这种求解方式在计算资源消耗和计算精度之间取得了较好的平衡，对于本文研究的综掘工作面巷道内通风风流与粉尘分布问题，能够稳定且高效地求解。在压力-速度耦合算法方面，选用SIMPLEC算法。SIMPLEC算法是SIMPLE（Semi-ImplicitMethodforPressure-LinkedEquations）算法的改进版本，它通过对压力修正方程的推导和处理，减少了迭代过程中的压力速度耦合误差，提高了收敛速度。在SIMPLEC算法中，对压力修正方程中的系数进行了简化和近似处理，使得计算过程更加简便，同时保持了较好的收敛特性。与其他压力-速度耦合算法相比，SIMPLEC算法在处理复杂流场和多相流问题时，能够更快地收敛到稳定解，尤其适用于本文中通风风流与粉尘气固两相流的模拟计算。对于离散格式的选择，动量方程采用二阶迎风离散格式。二阶迎风离散格式在处理对流项时，考虑了上游节点的信息，能够更准确地捕捉流体的流动特性。与一阶迎风离散格式相比，二阶迎风离散格式具有更高的精度，能够减少数值扩散误差，更真实地反映风流速度在巷道内的分布情况。在处理复杂的巷道几何形状和风流的复杂流动，如涡流、回流等现象时，二阶迎风离散格式能够提供更准确的模拟结果。湍动能k和湍流耗散率ε方程同样采用二阶迎风离散格式，以保证对湍流特性的准确模拟。二阶迎风离散格式在处理湍动能和湍流耗散率的输运过程中，能够更好地捕捉湍流的产生、发展和耗散机制，从而准确地描述风流的湍流特性，为研究通风风流与粉尘的相互作用提供可靠的基础。在迭代计算过程中，设置收敛残差为1×10⁻⁶。收敛残差是判断计算是否收敛的重要指标，它表示在迭代过程中，当前迭代步的计算结果与上一步计算结果之间的差异程度。当收敛残差小于设定值时，认为计算已经收敛，得到的结果是稳定可靠的。1×10⁻⁶的收敛残差设置能够保证计算结果具有较高的精度，确保模拟结果能够准确反映综掘工作面巷道内通风风流与粉尘的分布规律。同时，为了加快计算收敛速度，采用多重网格技术。多重网格技术通过在不同分辨率的网格上进行迭代计算，利用粗网格上的计算结果为细网格的迭代提供初始值，从而加速收敛过程。在粗网格上，计算量较小，能够快速得到一个大致的解，然后将这个解插值到细网格上，作为细网格迭代的初始值，这样可以减少细网格上的迭代次数，提高计算效率，缩短计算时间。四、通风风流分布规律的模拟结果与分析4.1不同通风方式下风流分布模拟结果在对综掘工作面巷道内通风风流分布规律的研究中，采用数值模拟方法对压入式通风、抽出式通风和压抽混合式通风这三种常见通风方式进行模拟，以深入探究不同通风方式下巷道内风流的分布特性。图4-1展示了压入式通风方式下巷道内风流速度矢量图和流线图。从速度矢量图中可以清晰地看到，新鲜空气由风筒高速射出，在巷道内形成一股强劲的射流。风筒出口附近的风流速度较大，可达15m/s，这是由于通风机提供的动力使空气具有较高的动能。随着风流在巷道内的传播，射流逐渐扩散，速度逐渐降低。在射流的作用下，巷道内形成了明显的回流区，回流区的存在使得部分污浊空气在巷道内循环，不利于空气的排出。在巷道壁面附近，由于壁面的摩擦阻力，风流速度明显减小，形成了边界层。从流线图中可以看出，风流从风筒出口射出后，沿着巷道壁面流动，然后在巷道中部形成回流，最终排出巷道。这种风流分布方式使得工作面附近的空气能够得到一定程度的更新，但也容易导致粉尘和有害气体在巷道内的扩散。【此处插入压入式通风风流速度矢量图和流线图】图4-2为抽出式通风方式下的模拟结果。在抽出式通风中，抽风机在风筒内形成负压，将工作面的含尘风流吸入风筒。从速度矢量图中可以观察到，风筒吸入口附近的风流速度较大，形成了较强的吸力。然而，由于抽风筒的有效吸程有限，在远离吸入口的区域，风流速度较小，导致部分区域的空气无法及时被抽出，容易造成粉尘和有害气体的积聚。在巷道的入口处，新鲜空气的进入受到一定的阻碍，风流速度相对较低。流线图显示，新鲜空气从巷道入口缓慢流入，经过工作面后被抽风机吸入风筒排出，整个巷道内的风流速度分布不均匀，通风效果受到一定影响。【此处插入抽出式通风风流速度矢量图和流线图】图4-3呈现了压抽混合式通风方式下的风流分布情况。在这种通风方式下，压风筒和抽风筒共同作用，使得巷道内的风流分布更加复杂。压风筒将新鲜空气压入工作面，抽风筒则在回风侧吸走含尘风流。从速度矢量图中可以看到，压风筒出口的风流速度较大，能够有效地将新鲜空气输送到工作面，同时，抽风筒吸入口附近的风流速度也较大，能够快速地将含尘风流排出。在压风筒和抽风筒之间的区域，风流速度相对较小，但整体上风流分布较为均匀。流线图表明，风流在压风筒和抽风筒的作用下，形成了较为规则的流动路径，减少了粉尘和有害气体在巷道内的积聚，通风效果相对较好。【此处插入压抽混合式通风风流速度矢量图和流线图】通过对三种通风方式下风流分布模拟结果的对比分析，可以看出不同通风方式对巷道内风流分布产生了显著影响。压入式通风有效射程远，能够快速排出工作面的乏风和粉尘，但容易使粉尘和有害气体在巷道内扩散；抽出式通风能保持巷道内空气新鲜，但有效吸程小，通风效果有限；压抽混合式通风综合了前两种通风方式的优点，在保证供风的同时，能缩短排出含尘风流的时间，缩小粉尘在巷道内的扩散范围，通风效果相对最佳。这些模拟结果为综掘工作面巷道通风系统的优化设计提供了重要依据，在实际应用中，应根据具体的地质条件、巷道布置和生产需求，合理选择通风方式，以提高通风效率，保障安全生产。4.2风流速度、压力等参数的分布特征在压入式通风方式下，巷道内风流速度呈现出明显的不均匀分布特征。图4-4展示了沿巷道中轴线方向的风流速度变化曲线。在风筒出口附近，风流速度急剧下降。风筒出口处的初始风速为15m/s，在距离风筒出口5m的位置，风流速度迅速降至8m/s左右，这是由于风流从风筒射出后，与周围空气发生强烈的掺混和动量交换，能量逐渐耗散，导致速度降低。随着距离的进一步增加，风流速度下降趋势逐渐变缓，在距离风筒出口20m处，风流速度降至4m/s左右，之后速度变化趋于平稳，在巷道末端风流速度稳定在1m/s左右。【此处插入压入式通风沿巷道中轴线风流速度变化曲线】在巷道横断面上，风流速度分布也存在显著差异。图4-5为距离风筒出口10m处的巷道横断面风流速度矢量图。可以看出，巷道中心区域的风流速度较大，而靠近巷道壁面处的风流速度较小。在巷道中心位置，风流速度可达5m/s，而在距离巷道壁面0.5m的范围内，风流速度降至1m/s以下。这是因为巷道壁面的摩擦阻力对风流产生了阻碍作用，使得靠近壁面的风流速度减小，形成了边界层。边界层的厚度随着风流在巷道内的流动逐渐增加，在巷道入口处边界层较薄，随着风流的发展，边界层逐渐增厚，对风流速度分布的影响也越来越大。【此处插入压入式通风距离风筒出口10m处巷道横断面风流速度矢量图】巷道内的风流压力分布同样呈现出一定的规律。在压入式通风中，风筒出口处的风流压力较高，随着风流在巷道内的传播，压力逐渐降低。图4-6为沿巷道中轴线方向的风流压力变化曲线。风筒出口处的静压为1050Pa，随着风流向巷道深处流动，静压逐渐减小，在距离风筒出口20m处，静压降至980Pa左右，之后静压继续缓慢下降，在巷道末端静压降至950Pa左右。风流压力的降低主要是由于风流在流动过程中克服巷道壁面的摩擦阻力以及与周围空气的掺混等因素，导致能量损失，压力降低。【此处插入压入式通风沿巷道中轴线风流压力变化曲线】在抽出式通风方式下，巷道内风流速度和压力分布与压入式通风有明显不同。图4-7展示了抽出式通风沿巷道中轴线方向的风流速度变化曲线。在抽风筒吸入口附近，风流速度迅速增大，在吸入口处风流速度可达12m/s，这是由于抽风机产生的负压作用，使得风流被快速吸入风筒。随着距离吸入口的增加，风流速度逐渐减小，在距离吸入口20m处，风流速度降至3m/s左右，之后速度继续缓慢下降，在巷道入口处风流速度降至0.5m/s左右。【此处插入抽出式通风沿巷道中轴线风流速度变化曲线】图4-8为距离抽风筒吸入口10m处的巷道横断面风流速度矢量图。可以看出，抽风筒吸入口附近的风流速度较大，形成了明显的抽吸区域。在吸入口周围一定范围内，风流速度方向指向吸入口，而在远离吸入口的区域，风流速度较小且分布较为均匀。在距离吸入口0.5m的范围内，风流速度可达8m/s，而在巷道其他区域，风流速度大多在1m/s以下。【此处插入抽出式通风距离抽风筒吸入口10m处巷道横断面风流速度矢量图】抽出式通风时，巷道内的风流压力分布也有其特点。图4-9为沿巷道中轴线方向的风流压力变化曲线。在抽风筒吸入口处，风流压力最低，为900Pa，这是由于抽风机的抽吸作用导致此处形成负压。随着距离吸入口的增加，风流压力逐渐升高，在距离吸入口20m处，压力升高至960Pa左右，之后压力继续缓慢升高，在巷道入口处压力升高至980Pa左右。风流压力的升高是因为风流从低压区域流向高压区域，在流动过程中，压力逐渐恢复。【此处插入抽出式通风沿巷道中轴线风流压力变化曲线】在压抽混合式通风方式下，巷道内风流速度和压力分布更为复杂。图4-10展示了压抽混合式通风沿巷道中轴线方向的风流速度变化曲线。在压风筒出口附近，风流速度较大，随着风流的传播，速度逐渐降低。在抽风筒吸入口附近，风流速度又迅速增大。压风筒出口处的初始风速为15m/s，在距离压风筒出口5m的位置，风流速度降至8m/s左右，之后速度缓慢下降。在距离抽风筒吸入口5m处，风流速度迅速增大至10m/s左右，之后随着距离吸入口的增加，速度逐渐减小，在巷道末端风流速度稳定在1m/s左右。【此处插入压抽混合式通风沿巷道中轴线风流速度变化曲线】图4-11为距离压风筒出口10m、抽风筒吸入口5m处的巷道横断面风流速度矢量图。可以看出，在压风筒出口附近，风流速度较大且方向指向巷道深处；在抽风筒吸入口附近，风流速度较大且方向指向吸入口；在两者之间的区域，风流速度相对较小且分布较为均匀。在压风筒出口中心位置，风流速度可达6m/s，在抽风筒吸入口周围0.5m范围内，风流速度可达8m/s，而在两者之间的区域，风流速度大多在2m/s左右。【此处插入压抽混合式通风距离压风筒出口10m、抽风筒吸入口5m处巷道横断面风流速度矢量图】图4-12为压抽混合式通风沿巷道中轴线方向的风流压力变化曲线。在压风筒出口处，风流压力较高，随着风流的传播，压力逐渐降低。在抽风筒吸入口处，风流压力最低，形成负压区域。压风筒出口处的静压为1050Pa，在距离压风筒出口20m处，静压降至980Pa左右。在抽风筒吸入口处，静压降至920Pa左右，之后随着距离吸入口的增加，压力逐渐升高，在巷道末端静压升高至950Pa左右。风流压力的这种变化是由于压风筒和抽风筒的共同作用，使得风流在巷道内的能量分布发生改变。【此处插入压抽混合式通风沿巷道中轴线风流压力变化曲线】综上所述，不同通风方式下巷道内风流速度、压力等参数的分布存在显著差异。通风方式对风流参数分布具有重要影响，在实际应用中，应根据具体情况合理选择通风方式，并优化通风参数，以获得良好的通风效果，降低粉尘浓度，保障矿井安全生产。4.3通风参数对风流分布的影响在综掘工作面巷道通风中，通风参数的变化对风流分布有着显著影响，深入研究这些影响对于优化通风系统设计、提高通风效果具有重要意义。风量是影响风流分布的关键参数之一。通过数值模拟，分别设置不同的风量工况，当风量从30m³/s增加到60m³/s时，巷道内的风流速度明显增大。在风筒出口附近，风流速度随着风量的增加而线性增加，从10m/s提升至20m/s左右。这是因为风量的增大意味着单位时间内进入巷道的空气量增多，在风筒出口处形成的射流能量更强，从而使风流速度增大。随着风流在巷道内的传播，风量的增加使得风流的有效射程也显著增加，风流能够更深入地到达巷道远端，对巷道内空气的更新和稀释作用增强。同时，较大的风量有助于减少巷道内的涡流和回流区域，使风流分布更加均匀。在风量为30m³/s时，巷道内存在明显的涡流区，涡流区面积占巷道总面积的20%左右；而当风量增加到60m³/s时，涡流区面积减小至10%左右，风流的稳定性得到提高，更有利于将粉尘和有害气体排出巷道。风筒直径的改变也会对风流分布产生重要影响。模拟中采用了直径为0.6m、0.8m和1.0m的风筒。当风筒直径从0.6m增大到0.8m时，风筒出口处的风流速度略有降低，从15m/s降至13m/s左右。这是因为在风量不变的情况下，风筒直径增大，风流的过流面积增加，根据流量公式Q=vA（Q为流量，v为流速，A为过流面积），流速会相应减小。然而，随着风筒直径的增大，风流的能量损失减小，风筒内的压力损失降低，使得风流在巷道内的有效射程增加。在风筒直径为0.6m时，风流的有效射程为15m左右；当风筒直径增大到0.8m时，有效射程增加到20m左右。此外，较大直径的风筒还能使风流在巷道内的分布更加均匀，减少了风流的偏析现象。当风筒直径为1.0m时，巷道横断面上的风流速度分布更加均匀，速度差值减小，有利于提高通风效果。风筒位置的变化同样会对风流分布产生显著影响。分别模拟了风筒距离巷道顶板0.3m、0.5m和0.7m三种位置情况。当风筒距离顶板0.3m时，风筒出口射流靠近顶板，在顶板附近形成高速风流区域，而巷道底部的风流速度相对较小。随着风筒距离顶板距离的增加，如增大到0.5m时，风流在巷道横断面上的分布更加均匀，巷道顶部和底部的风流速度差值减小。当风筒距离顶板0.7m时，风流在巷道底部的速度有所增加，但由于风筒位置较高，风流在到达巷道底部时能量有所损失，导致巷道底部的风流速度仍低于巷道中部和顶部。此外，风筒位置的变化还会影响巷道内涡流的形成和分布。当风筒距离顶板较近时，在风筒出口下方容易形成涡流；而当风筒位置适当调整后，涡流的强度和范围会减小，有利于改善通风效果。综上所述，风量、风筒直径和风筒位置等通风参数对综掘工作面巷道内的风流分布有着重要影响。在实际的通风系统设计和运行中，应根据巷道的具体情况和通风要求，合理调整通风参数，以获得最佳的风流分布，提高通风效率，降低粉尘浓度，保障矿井安全生产。五、粉尘分布规律的模拟结果与分析5.1粉尘浓度分布模拟结果通过数值模拟，得到了不同通风方式和通风参数下综掘工作面巷道内粉尘质量浓度云图和等值线图，这些图直观地展示了粉尘在巷道内的分布情况。图5-1为压入式通风方式下，距离掘进机截割头10s时的粉尘质量浓度云图和等值线图。从云图中可以明显看出，掘进机截割头附近区域的粉尘浓度极高，呈现出明显的高浓度核心区域。这是因为截割头在切割煤岩时，会产生大量的粉尘，这些粉尘在初始阶段还未充分扩散，导致浓度积聚。在距离截割头较近的区域，粉尘质量浓度可达3000mg/m³以上。随着距离的增加，粉尘浓度逐渐降低，但在巷道的大部分区域，粉尘浓度仍然较高，尤其是在风筒射流的影响范围内，由于风流的携带作用，粉尘被输送到较远的位置。从等值线图中可以更清晰地看到粉尘浓度的梯度变化，等值线在截割头附近较为密集，表明此处浓度变化剧烈，而在远离截割头的区域，等值线逐渐稀疏，浓度变化相对平缓。在距离截割头20m处，粉尘质量浓度降至1000mg/m³左右，在巷道末端，粉尘质量浓度仍维持在200mg/m³左右。【此处插入压入式通风10s时粉尘质量浓度云图和等值线图】图5-2展示了抽出式通风方式下，相同时间点的粉尘分布情况。与压入式通风不同，抽出式通风时，抽风筒吸入口附近的粉尘浓度较高，形成了一个高浓度区域。这是因为抽风机的抽吸作用，使得含尘风流向吸入口汇聚。在吸入口周围一定范围内，粉尘质量浓度可达2500mg/m³以上。然而，在远离吸入口的区域，由于抽风筒的有效吸程有限，粉尘浓度下降较为缓慢，导致部分区域的粉尘浓度仍然较高。在巷道入口处，由于新鲜空气的进入，粉尘浓度相对较低，但仍有一定量的粉尘存在，此处粉尘质量浓度约为500mg/m³。从等值线图可以看出，粉尘浓度在吸入口附近呈现出向中心汇聚的趋势，等值线较为密集，而在其他区域，等值线分布相对较为均匀，但浓度值仍然较高，说明抽出式通风在控制粉尘扩散方面存在一定的局限性。【此处插入抽出式通风10s时粉尘质量浓度云图和等值线图】图5-3为压抽混合式通风方式下的粉尘浓度分布。在这种通风方式下，由于压风筒和抽风筒的共同作用，巷道内的粉尘分布得到了明显改善。压风筒将新鲜空气压入工作面，有效地稀释了截割头附近的粉尘浓度，使得高浓度区域的范围明显减小。抽风筒则及时将含尘风流抽出，进一步降低了巷道内的粉尘浓度。从云图中可以看到，截割头附近的粉尘质量浓度虽然仍然较高，但相比压入式通风和抽出式通风，有了显著的降低，最高浓度约为1500mg/m³。在压风筒和抽风筒之间的区域，粉尘浓度相对较低，大部分区域的粉尘质量浓度在500mg/m³以下。从等值线图可以清晰地看到，粉尘浓度在巷道内呈现出较为均匀的分布，等值线的分布也相对较为稀疏，说明压抽混合式通风能够有效地控制粉尘的扩散，降低巷道内的粉尘浓度。【此处插入压抽混合式通风10s时粉尘质量浓度云图和等值线图】随着时间的推移，粉尘在巷道内的分布也会发生变化。图5-4展示了压入式通风方式下，距离掘进机截割头30s时的粉尘质量浓度云图和等值线图。与10s时相比，粉尘浓度的高值区域进一步扩大，这是因为随着时间的增加，粉尘有更多的时间扩散。在截割头附近，粉尘质量浓度仍然维持在较高水平，约为2500mg/m³，但在巷道的更远处，粉尘浓度也有所升高。在距离截割头30m处，粉尘质量浓度从10s时的800mg/m³左右升高到1200mg/m³左右，在巷道末端，粉尘质量浓度也从200mg/m³左右升高到300mg/m³左右。这表明在压入式通风方式下，随着时间的延长，粉尘的扩散范围逐渐增大，浓度分布也更加均匀，但整体浓度水平有所上升。【此处插入压入式通风30s时粉尘质量浓度云图和等值线图】图5-5为抽出式通风方式下，30s时的粉尘分布情况。此时，抽风筒吸入口附近的高浓度区域范围有所扩大，粉尘质量浓度进一步升高，最高可达3000mg/m³以上。这是因为随着时间的增加，更多的粉尘被抽吸到吸入口附近。在巷道的其他区域，粉尘浓度也有所上升，尤其是在远离吸入口的区域，由于抽风效果的减弱，粉尘浓度升高较为明显。在巷道入口处，粉尘质量浓度从10s时的500mg/m³左右升高到700mg/m³左右，说明抽出式通风在长时间运行时，对粉尘的控制效果逐渐变差，粉尘浓度会逐渐升高。【此处插入抽出式通风30s时粉尘质量浓度云图和等值线图】图5-6展示了压抽混合式通风方式下，30s时的粉尘浓度分布。与10s时相比，虽然截割头附近的粉尘浓度仍然较高，但整体浓度水平有所下降，最高浓度约为1200mg/m³。在压风筒和抽风筒之间的区域，粉尘浓度进一步降低，大部分区域的粉尘质量浓度在300mg/m³以下。这表明压抽混合式通风在长时间运行时，仍然能够有效地控制粉尘的扩散，降低巷道内的粉尘浓度，保持较好的通风降尘效果。【此处插入压抽混合式通风30s时粉尘质量浓度云图和等值线图】通过对不同通风方式和不同时间下粉尘浓度分布模拟结果的分析，可以看出通风方式对粉尘分布有着显著的影响。压入式通风容易使粉尘在巷道内扩散，随着时间的增加，粉尘浓度会逐渐升高；抽出式通风在控制粉尘扩散方面存在一定局限性，长时间运行时，粉尘浓度也会升高；而压抽混合式通风能够有效地控制粉尘的扩散，降低巷道内的粉尘浓度，在不同时间下都能保持较好的通风降尘效果。这些模拟结果为优化通风系统和制定有效的粉尘防治措施提供了重要依据。5.2粉尘运动轨迹分析利用离散相模型（DPM）对粉尘颗粒在风流作用下的运动轨迹进行追踪，能够深入了解粉尘在巷道内的扩散和沉降规律，这对于制定有效的粉尘防治措施具有关键意义。在压入式通风方式下，从粉尘颗粒的运动轨迹图（图5-7）可以清晰地看到，在掘进机截割头产生粉尘的瞬间，粉尘颗粒在截割头的高速旋转和风流的共同作用下，以较高的初速度向四周飞溅。由于风筒出口射流的影响，大部分粉尘颗粒被风流携带，沿着射流方向向巷道深处运动。在靠近风筒出口的区域，粉尘颗粒的运动速度较快，轨迹较为集中，这是因为该区域风流速度较大，对粉尘颗粒的携带能力较强。随着粉尘颗粒向巷道远端移动，风流速度逐渐降低，粉尘颗粒之间的相互碰撞和摩擦加剧，导致部分粉尘颗粒的运动方向发生改变，轨迹变得较为分散。在巷道壁面附近，由于壁面的摩擦阻力和边界层效应，粉尘颗粒的运动速度明显减小，部分粉尘颗粒会与壁面发生碰撞，然后反弹或沉降到壁面上。此外，在巷道的顶部和底部，由于风流速度相对较低，粉尘颗粒的沉降现象较为明显，尤其是粒径较大的粉尘颗粒，更容易在重力作用下迅速沉降到巷道底部。【此处插入压入式通风粉尘颗粒运动轨迹图】在抽出式通风方式下，粉尘颗粒的运动轨迹呈现出不同的特征（图5-8）。在截割头附近产生的粉尘，受到抽风机负压的影响，迅速向抽风筒吸入口方向运动。在吸入口周围一定范围内，粉尘颗粒的运动轨迹较为集中，速度较大，这是因为抽风机的抽吸作用使得该区域形成了较强的吸力。然而，在远离吸入口的区域，由于抽风筒的有效吸程有限，风流速度较小，粉尘颗粒的运动速度也随之降低，轨迹变得较为分散。部分粉尘颗粒在风流的作用下，会在巷道内形成环流，难以被及时抽出，导致粉尘在巷道内积聚。与压入式通风类似，在巷道壁面附近，粉尘颗粒也会与壁面发生碰撞和沉降，但由于抽出式通风时巷道内的风流速度分布不均匀，壁面附近的粉尘沉降情况相对较为复杂，不同位置的沉降程度存在差异。【此处插入抽出式通风粉尘颗粒运动轨迹图】对于压抽混合式通风方式，粉尘颗粒的运动轨迹更为复杂（图5-9）。在压风筒出口附近，粉尘颗粒在新鲜风流的推动下，向工作面方向运动，与抽出式通风时向吸入口运动的方向相反。在工作面区域，压入的新鲜风流与抽出的含尘风流相互作用，使得粉尘颗粒的运动轨迹发生弯曲和交织。部分粉尘颗粒被压入的新鲜风流稀释