基于多参数耦合的煤堆压实防治自燃不适风速的数值模拟与实验研究

基于多参数耦合的煤堆压实防治自燃不适风速的数值模拟与实验研究一、引言1.1研究背景与意义煤炭作为一种重要的基础能源，在全球能源结构中占据着举足轻重的地位。我国是煤炭生产和消费大国，煤炭在一次能源生产和消费结构中的占比长期保持在较高水平。然而，煤炭在储存和堆放过程中，煤堆自燃问题却一直是困扰煤炭行业的一个难题。煤堆自燃不仅会造成煤炭资源的大量浪费，还会引发一系列严重的后果。从经济角度来看，煤堆自燃导致煤炭发热量降低，品质下降，直接影响煤炭的销售价格和使用价值，给企业带来巨大的经济损失。据相关统计，每年因煤堆自燃造成的煤炭损耗在数百万吨以上，经济损失高达数十亿元。从安全角度而言，煤堆自燃容易引发火灾，火势一旦失控，可能会对周边的建筑物、设备以及人员生命安全构成严重威胁。例如，2018年某大型煤场发生煤堆自燃引发的火灾事故，造成了周边多座仓库被烧毁，直接经济损失超过千万元，同时还导致了数名消防人员受伤。此外，煤堆自燃过程中会释放出大量的有害气体，如二氧化硫（SO_2）、一氧化碳（CO）、硫化氢（H_2S）等，这些气体不仅会污染空气，危害周边居民的身体健康，还会对大气环境造成长期的破坏，加剧酸雨、雾霾等环境问题。为了有效防治煤堆自燃，目前常用的方法有多种。周期性压实通过减小煤堆孔隙率，降低氧气进入量，从而抑制煤堆氧化；表面覆盖利用覆盖材料隔绝氧气与煤体接触；细煤覆盖在煤堆表面形成一层保护屏障；惰化向煤堆中注入惰性气体，降低氧气浓度；风障改变风流方向，减少煤堆表面风速，降低氧气供应；喷水增加煤堆湿度，带走部分热量；灌浆填充煤堆孔隙，隔绝氧气。然而，这些方法在实际应用中都存在一定的局限性。例如，使用阻化剂虽然能抑制煤炭氧化，但可能会降低煤质，且成本较高；表面喷水受环境温湿度及风速影响较大，有时反而会加速煤堆自燃。在众多防治方法中，煤堆压实技术因其操作相对简单、成本较低，在工程实践中得到了广泛应用。然而，大量的实践经验和研究表明，煤堆压实技术并非在所有情况下都能有效防治自燃，在高风速条件下，压实后的煤堆有时反而更容易发生自燃。这是因为风速的变化会显著影响煤堆内部的漏风速度和进气量。当风速较大时，压实后的煤堆孔隙率降低，透气性变差，使得煤堆内部产生的热量难以散发出去，同时氧气又能通过漏风不断进入煤堆，从而导致煤堆内部温度不断升高，加速煤的氧化和自燃进程。因此，深入研究煤堆压实防治自燃与不适风速之间的关系，对于提高煤堆自燃防治的有效性和科学性具有至关重要的意义。一方面，通过对不同堆放参数（如煤堆形状、高度、顶底部宽度、角度、孔隙率等）下煤堆压实防治自燃的不适风速进行模拟研究，可以准确掌握煤堆在不同条件下的自燃特性，为煤堆的合理堆放和压实操作提供科学依据。例如，明确在何种风速范围内压实煤堆是安全有效的，以及在不同风速条件下应如何调整煤堆的堆放参数，以避免因压实不当而引发自燃。另一方面，这一研究成果还可以为煤场的规划设计和日常管理提供重要参考，有助于制定更加完善的煤堆自燃防治措施和应急预案，从而有效降低煤堆自燃的风险，保障煤炭资源的安全储存和利用，减少经济损失和环境污染，促进煤炭行业的可持续发展。1.2国内外研究现状煤堆自燃是一个复杂的物理化学过程，受到多种因素的综合影响。国内外学者针对煤堆自燃的机理、影响因素以及防治措施等方面开展了大量的研究工作，取得了一系列有价值的成果。在煤堆自燃机理研究方面，国外学者较早开展了相关工作。例如，Alazmi和Hooman研究发现风流对煤堆的影响存在从自然对流到混合对流再到强制对流的转变过程，并且在这个过程中，风速的变化对煤堆自燃有着显著影响。F.Akguna等建立了煤堆自热自燃的二维非稳态计算模型，通过该模型深入研究了煤堆高度、堆放时间、角度、粒径孔隙率等因素对煤堆自热自燃的影响，指出孔隙率在0.15-0.2时，由于煤堆透气性降低，自热升温较慢。Ejlali在多风速条件下对煤堆自燃进行研究，得出进风量是影响煤堆最高温度变化较为显著的因素。A.Ejlali等进一步强调了孔隙率是影响煤堆内部自热-自燃和最高温度变化的关键参数，并提出煤堆压实并非在所有情况下都能有效保护煤堆，有时可能会起到相反的作用。国内学者在煤堆自燃机理研究方面也做出了重要贡献。徐精彩等认为煤堆压实是防治煤堆自燃的首选方法，通过压实可以减小煤堆孔隙率，降低氧气进入量，从而抑制煤堆氧化。董希琳等则指出，任何可以影响煤堆自热-自燃的参数发生微小变化，都可能导致煤堆内部系统发生较大的非连续性变化，进而引起煤堆的不安全。张瑞新、谢和平、谢之康等通过试验证明煤堆自燃受环境影响较大，应优先采取预防措施，如合理选择煤堆的堆放位置、控制煤堆的高度和形状等。在煤堆堆放参数对自燃影响的研究中，国内外学者主要关注煤堆的形状、高度、顶底部宽度、角度以及孔隙率等参数。国外研究中，部分学者通过实验和数值模拟，分析了不同堆放参数下煤堆内部的气流分布和温度场变化，发现煤堆的高度和坡度会影响煤堆表面的风速和气流方向，进而影响煤堆的散热和氧气供应。例如，研究表明较高的煤堆在相同风速下，内部热量更难散发，更容易达到自燃温度；而较大的坡度会使煤堆表面风速增大，增加氧气进入量，也可能促进煤堆自燃。国内学者则更侧重于研究这些参数之间的相互关系以及它们对煤堆自燃的综合影响。有研究通过构建多参数耦合的煤堆自燃模型，模拟不同堆放参数组合下煤堆的自燃过程，发现孔隙率与煤堆高度、形状等参数相互作用，共同影响煤堆的自燃特性。例如，在相同孔隙率下，不同形状的煤堆其内部的气流通道和散热路径不同，导致自燃倾向性也有所差异。关于不适风速与煤堆压实防治自燃关系的研究，目前相关报道相对较少。国外一些研究尝试通过实验测量不同风速下压实煤堆的温度变化和气体成分，初步探索不适风速的范围，但尚未形成系统的理论和方法。国内学者董子文等基于煤堆自燃理论，使用COMSOLMultiphysics数值仿真软件建立煤堆自燃模型，对不同堆放规格尺寸煤堆进行压实后温度变化和自然发火期变化进行了分析，从温度变化和自然发火期两个角度研究了煤堆最低不适用压实风速，建立了与煤堆孔隙率相关的最低不适用压实的风速方程。然而，该研究仅考虑了有限的几个堆放参数，对于煤堆形状、角度等参数对不适风速的影响研究不够深入，且在实际应用中，该模型的准确性和普适性还有待进一步验证。尽管国内外在煤堆自燃研究方面取得了一定进展，但仍存在一些不足之处。现有研究对于煤堆自燃机理的认识还不够全面和深入，特别是在多因素耦合作用下煤堆自燃的微观过程和宏观表现之间的联系，尚未完全明确。在煤堆堆放参数对自燃影响的研究中，大多数研究仅考虑单一或少数几个参数的变化，缺乏对多个参数同时变化时相互作用机制的深入研究。对于不适风速与煤堆压实防治自燃关系的研究还处于起步阶段，缺乏系统的理论和实验研究，无法准确地确定不同堆放参数下煤堆压实防治自燃的不适风速范围，难以满足实际工程的需求。1.3研究内容与方法1.3.1研究内容本研究旨在深入探究堆放参数对煤堆压实防治自燃的不适风速影响，具体研究内容包括以下几个方面：不同堆放参数下煤堆压实防治自燃的不适风速研究：系统分析煤堆形状（如圆锥体、长方体、棱台体等）、高度（设置不同高度梯度，如5m、10m、15m等）、顶底部宽度（改变顶底部宽度比例，如1:2、1:3、1:4等）、角度（调整坡面角度，如30°、45°、60°等）、孔隙率（通过压实程度控制，设定孔隙率范围为0.1-0.4）等堆放参数，对煤堆压实防治自燃的不适风速进行全面研究。利用数值模拟软件，建立不同堆放参数组合下的煤堆模型，模拟不同风速条件（设定风速范围为0-20m/s，以1m/s为步长）下煤堆内部的气流分布、温度场变化以及氧气浓度分布等情况，确定不同堆放参数下煤堆压实防治自燃的不适风速范围。煤堆堆放参数与不适风速关系模型的构建：基于模拟结果，运用数据拟合和回归分析等方法，构建煤堆堆放参数与不适风速之间的定量关系模型。通过对模型的验证和优化，提高模型的准确性和可靠性，使其能够准确预测不同堆放参数下煤堆压实防治自燃的不适风速。在构建模型过程中，考虑各堆放参数之间的相互作用和耦合关系，例如煤堆高度和孔隙率对不适风速的协同影响，以及煤堆形状和角度对气流分布的综合作用等。模型验证与分析：采用实验数据对构建的关系模型进行验证，对比模拟结果与实验数据，评估模型的准确性和有效性。对模型进行灵敏度分析，确定各堆放参数对不适风速的影响程度和敏感性，找出影响不适风速的关键参数。例如，通过改变某一参数的值，观察不适风速的变化情况，从而确定该参数对不适风速的影响程度。根据模型验证和分析结果，提出针对性的煤堆自燃防治建议，为实际工程应用提供科学依据。例如，对于敏感性较高的参数，在煤堆堆放和压实过程中应严格控制其取值范围，以降低煤堆自燃的风险。1.3.2研究方法本研究综合采用数值模拟和实验研究相结合的方法，确保研究结果的准确性和可靠性。具体研究方法如下：数值模拟方法：选用COMSOLMultiphysics、ANSYSFluent等专业数值模拟软件，建立煤堆自燃的三维物理模型。在模型中，考虑煤堆的热传导、对流换热、氧气扩散、煤氧化学反应等物理过程，通过设置合理的边界条件和初始条件，模拟不同堆放参数和风速条件下煤堆内部的物理场变化。利用软件的后处理功能，获取煤堆内部的温度分布、气流速度分布、氧气浓度分布等数据，为分析煤堆压实防治自燃的不适风速提供依据。例如，通过绘制不同时刻煤堆内部的温度云图，直观地展示煤堆温度随时间和空间的变化情况。实验研究方法：设计并搭建煤堆自燃实验平台，该平台应具备模拟不同风速、控制煤堆堆放参数、监测煤堆内部温度和气体成分等功能。采用实际煤样，按照不同的堆放参数（如不同的高度、形状、孔隙率等）进行煤堆堆放实验。在实验过程中，通过改变风速，利用热电偶、气体传感器等设备实时监测煤堆内部不同位置的温度、氧气浓度、一氧化碳浓度等参数的变化。对比不同风速下煤堆的温度变化和气体成分变化，确定煤堆压实防治自燃的不适风速。例如，当煤堆温度在某一风速下迅速升高，且一氧化碳浓度明显增加时，可判断该风速可能为不适风速。将实验结果与数值模拟结果进行对比分析，验证数值模拟模型的准确性，为进一步优化模型提供实验依据。二、煤堆自燃及压实防治原理2.1煤堆自燃机理煤堆自燃是一个复杂的物理化学过程，其机理至今尚未完全明确。目前，煤氧复合作用学说被广泛接受，该学说认为煤堆自燃是煤与氧气相互作用的结果。在常温下，煤中的碳（C）、氢（H）等元素会与空气中的氧气发生反应，生成一氧化碳（CO）、甲烷（CH_4）及其他烷烃等可燃物。这些氧化反应是放热过程，如果产生的热量不能及时散发出去，就会在煤堆内部积聚，导致煤堆内部温度逐渐升高。而温度的升高又会进一步加速煤的氧化速度，使其放出更多的可燃物质和热量。当热量不断聚集，温度上升到一定值时，就会引发可燃物质的燃烧，从而导致煤堆自燃。煤堆自燃过程通常可分为三个阶段，各阶段具有不同的特征和所需条件。潜伏期：在这一阶段，煤的低温氧化反应处于缓慢状态。煤表面与空气中的氧气相互作用，主要以物理吸附的方式在煤的表面形成过氧络合物。此过程中生成的热量较少，煤的温度变化也较小，几乎没有明显的宏观现象。潜伏期的长短主要取决于煤的性质，如煤化程度、硫铁矿含量等。一般来说，煤化程度越低、硫铁矿含量越高，潜伏期越短。此外，环境温度和湿度也会对潜伏期产生影响，较高的环境温度和湿度会缩短潜伏期。自热期：经过潜伏期后，被活化的煤炭能够更快地吸收氧气，煤氧化的速度加快。氧化产生的热量大幅增加，如果热量不能及时散发出去，煤堆温度就会逐渐上升。当达到煤的自热临界温度（一般为60-80℃）时，煤的吸氧能力将快速加大，致使煤的氧化速度急剧上升，煤温同步快速上升。同时，会生成一氧化碳（CO）、二氧化碳（CO_2）、甲烷（CH_4）及其他烷烃等可燃物质。自热期的发展受到多种因素的影响，除了煤的性质外，煤堆的透气性和散热条件起着关键作用。如果煤堆孔隙率较大，透气性好，氧气容易进入煤堆内部，同时热量也容易散发出去，自热期就会相对较长；反之，如果煤堆孔隙率小，透气性差，热量积聚快，自热期就会缩短。燃烧期：当煤温达到着火温度（一般为210-350℃，不同煤种着火温度有所差异）时，就会引发煤炭自燃。此时，有煤烟和明火冒出，煤堆进入燃烧期。在燃烧期，煤的燃烧反应剧烈，会释放出大量的热量和有害气体。影响燃烧期的主要因素是氧气的供应和可燃物质的浓度。充足的氧气供应和较高的可燃物质浓度会使燃烧更加剧烈。此外，煤堆的几何形状和尺寸也会对燃烧期产生影响，较大的煤堆在燃烧时可能会形成更复杂的燃烧区域和火势发展情况。2.2煤堆压实防治自燃原理煤堆压实是防治煤堆自燃的一种重要方法，其原理主要基于对煤堆内部孔隙结构和氧气进入量的控制。煤堆在自然堆放状态下，内部存在大量的孔隙，这些孔隙相互连通，形成了气体流通的通道。空气中的氧气可以通过这些孔隙自由进入煤堆内部，与煤发生氧化反应。同时，煤堆内部由于氧化反应产生的热量也会通过孔隙与外界进行热交换。当对煤堆进行压实操作时，煤颗粒之间的排列更加紧密，煤堆的孔隙率显著降低。孔隙率的减小使得煤堆内部的气体流通通道变得狭窄和曲折，从而增加了氧气进入煤堆的阻力。根据气体扩散定律，气体的扩散速率与扩散系数成正比，与扩散距离成反比。在压实后的煤堆中，氧气需要经过更长的扩散路径才能到达煤堆内部的反应区域，这大大减缓了氧气的扩散速度，降低了煤与氧气的接触机会，从而抑制了煤的氧化反应速率。例如，当煤堆孔隙率从0.3降低到0.2时，氧气的扩散系数可能会降低30%-50%，使得煤的氧化速度明显下降。此外，压实还能减少煤堆内部的漏风现象。漏风是指空气在煤堆内部的不合理流动，它会加速氧气的进入和热量的传递。在未压实的煤堆中，由于孔隙分布不均匀，容易形成局部的漏风通道，使得氧气能够快速进入煤堆内部，同时也会带走部分热量，不利于热量的积聚和自燃的发生。而压实后的煤堆，孔隙结构更加均匀，漏风现象得到有效抑制，煤堆内部的氧化反应更加均匀和缓慢，减少了局部热点的形成，降低了煤堆自燃的风险。从热量积聚的角度来看，压实后的煤堆内部孔隙率减小，空气的流动性变差，这使得煤堆内部因氧化反应产生的热量难以通过空气对流的方式快速散发出去。在一定程度上，这有助于保持煤堆内部的低温状态，因为缓慢的热量积聚不易使煤堆温度升高到自燃的临界温度。然而，当风速较大时，这种热量积聚的特性可能会带来负面影响。由于压实后煤堆透气性差，外界大风带来的热量和氧气难以有效排出，导致煤堆内部温度迅速升高，反而增加了自燃的可能性。这也正是研究煤堆压实防治自燃与不适风速关系的重要原因之一。2.3风速对煤堆自燃的影响机制风速作为影响煤堆自燃的重要环境因素之一，其对煤堆自燃的影响机制较为复杂，主要通过改变煤堆的进气量、氧气量、热散失和积聚等方面来影响煤堆的自燃进程。从进气量角度来看，风速的变化会直接影响煤堆表面的风压，进而改变煤堆内部的漏风速度和进气量。当风速较低时，煤堆表面的风压较小，煤堆内部的漏风速度较慢，进气量相对较少。在这种情况下，煤与氧气的接触机会有限，煤的氧化反应速率较慢，产生的热量也相对较少。同时，由于进气量少，煤堆内部因氧化反应产生的热量难以被快速带走，热量逐渐积聚，可能会使煤堆内部温度升高。然而，如果风速过低，煤堆内部的热量积聚速度可能会超过热量散失速度，导致煤堆温度持续上升，当达到一定温度时，就会引发煤堆自燃。随着风速的增加，煤堆表面的风压增大，煤堆内部的漏风速度加快，进气量显著增加。充足的氧气供应使得煤与氧气的接触更加充分，煤的氧化反应速率加快，产生的热量也随之增多。在一定范围内，增加的风速可以及时带走煤堆内部产生的热量，使得煤堆内部温度保持在较低水平，抑制煤堆自燃的发生。例如，当风速在2-5m/s时，煤堆内部的热量能够得到有效散发，煤堆温度相对稳定，自燃风险较低。但当风速过高时，情况则会发生变化。过高的风速虽然能带来更多的氧气，但也会使煤堆内部的热量难以积聚。对于压实后的煤堆，由于其孔隙率降低，透气性变差，过高的风速使得外界的热量和氧气快速进入煤堆，而煤堆内部产生的热量却无法及时排出。这就导致煤堆内部温度迅速升高，加速煤的氧化和自燃进程。研究表明，当风速超过10m/s时，对于一些孔隙率较低的压实煤堆，自燃的风险会明显增加。风速对煤堆内部的氧气分布也有重要影响。在低风速条件下，氧气在煤堆内部的扩散主要依靠分子扩散，扩散速度较慢。这使得煤堆内部不同位置的氧气浓度分布相对不均匀，靠近煤堆表面的区域氧气浓度较高，而煤堆内部深处的氧气浓度较低。这种氧气浓度的差异会导致煤的氧化反应主要集中在煤堆表面附近，煤堆内部的氧化反应相对较弱。随着风速的增大，氧气在煤堆内部的扩散方式逐渐从分子扩散转变为对流扩散，扩散速度大大加快。这使得煤堆内部的氧气浓度分布更加均匀，煤与氧气能够在更广泛的范围内接触，从而加速煤的氧化反应。然而，当风速过大时，虽然氧气供应充足，但由于热量散失过快，煤堆内部的温度难以升高到自燃的临界温度，在一定程度上会抑制煤堆自燃。但对于压实煤堆，过高风速带来的大量氧气会加剧煤堆内部的氧化反应，即便热量有部分散失，仍可能因氧化反应过于剧烈而导致温度升高引发自燃。在热散失和积聚方面，风速起着关键的调节作用。当风速较小时，煤堆内部的热量主要通过热传导和少量的自然对流方式与外界进行热交换。由于自然对流较弱，热量散失速度较慢，容易在煤堆内部积聚。例如，在无风或微风环境下，煤堆内部的热量积聚速度可能会超过热量散失速度，导致煤堆温度逐渐升高。当风速适中时，强制对流成为煤堆内部热量传递的主要方式。此时，风速能够有效地将煤堆内部因氧化反应产生的热量带出煤堆，使得煤堆内部的温度保持在相对稳定的状态。例如，在风速为3-6m/s的情况下，煤堆内部的热量能够及时散发，煤堆温度不会出现明显的上升，从而降低了煤堆自燃的风险。而当风速过高时，虽然热量散失速度加快，但同时也会带来大量的外界热量。对于压实煤堆，由于其散热能力有限，过多的外界热量进入煤堆后无法及时排出，会导致煤堆内部温度升高。此外，过高的风速还可能破坏煤堆表面的覆盖层或保护层（如果有），进一步加剧煤堆内部的氧化和热量积聚，增加煤堆自燃的可能性。三、研究方法与模型建立3.1数值模拟软件介绍本研究选用COMSOLMultiphysics软件进行数值模拟分析。COMSOLMultiphysics是一款功能强大的多物理场仿真软件，在科研和工程领域有着广泛的应用。它基于有限元方法，能够精确地求解各种复杂的物理问题。COMSOLMultiphysics的核心优势在于其卓越的多物理场耦合模拟能力。在煤堆自燃的研究中，涉及到多种物理过程的相互作用，如热传导、对流换热、氧气扩散以及煤氧化学反应等。该软件能够将这些不同的物理场进行直接耦合，全面、准确地模拟煤堆内部的真实物理现象。例如，在模拟煤堆内部的温度分布时，它可以同时考虑热传导和对流换热对热量传递的影响，以及氧气扩散和煤氧化学反应对热量产生的作用。通过这种多物理场耦合的模拟方式，能够得到更符合实际情况的结果，为研究提供更可靠的依据。软件提供了直观易用的图形用户界面（GUI），即使对于没有深厚编程基础的研究人员，也能轻松上手。在建模过程中，用户可以通过GUI方便地定义煤堆的几何形状、设置材料属性、添加物理场方程以及设定边界条件和初始条件等。例如，在定义煤堆的几何形状时，用户只需在GUI中选择相应的几何模型（如圆锥体、长方体、棱台体等），并输入相关的尺寸参数，即可快速创建煤堆的几何模型。在设置材料属性时，用户可以直接在GUI中输入煤的密度、比热容、导热系数等参数，软件会自动将这些参数应用到模型中。COMSOLMultiphysics拥有丰富的物理模型库，涵盖了电磁学、结构力学、声学、流体流动、传热以及化学反应等多个领域。在煤堆自燃模拟中，可以直接调用其中的传热、流体流动和化学反应等物理模型，大大节省了建模时间和工作量。以传热模型为例，软件提供了多种传热方式的模拟选项，包括热传导、对流换热和辐射换热等，用户可以根据实际情况选择合适的传热模型进行模拟。在化学反应模型方面，软件支持多种反应动力学模型，能够准确地描述煤氧化学反应的过程。软件还具备强大的后处理功能，能够对模拟结果进行全面、深入的分析和可视化展示。用户可以通过后处理模块生成各种数据图表（如温度随时间变化曲线、氧气浓度分布曲线等）、云图（如煤堆内部温度云图、氧气浓度云图等）和动画（如煤堆温度随时间变化的动画），直观地展示煤堆内部物理量的分布和变化情况。这些可视化结果有助于研究人员更清晰地理解模拟结果，发现其中的规律和趋势。例如，通过观察煤堆内部温度云图，研究人员可以直观地看到煤堆内部高温区域的分布位置和范围，以及温度随时间的变化情况，从而分析煤堆自燃的发展趋势。在计算效率方面，COMSOLMultiphysics采用了先进的数值算法和高效的求解器，能够快速、准确地求解大规模的复杂模型。对于煤堆自燃这种涉及多物理场耦合和复杂几何形状的模拟问题，该软件能够在较短的时间内得到精确的结果。同时，软件还支持并行计算，能够充分利用多核处理器的计算资源，进一步提高计算效率。例如，在模拟大型煤堆的自燃过程时，通过并行计算可以将计算时间缩短数倍，大大提高了研究效率。3.2煤堆模型的建立与参数设定为了准确模拟不同堆放参数下煤堆压实防治自燃的不适风速，本研究建立了三维煤堆模型。在实际的煤场中，煤堆的形状多种多样，常见的有圆锥体、长方体、棱台体等。考虑到实际情况以及模型的代表性，本研究选择了圆锥体和棱台体作为煤堆的几何形状进行建模。对于圆锥体煤堆模型，其底面半径r设置为5m，高度h分别设置为5m、10m、15m，以研究不同高度对煤堆压实防治自燃的不适风速的影响。圆锥体的母线与底面的夹角\theta通过公式\theta=\arctan(\frac{h}{r})计算得出，在不同高度设置下，夹角也相应变化。例如，当高度为5m时，夹角\theta=\arctan(1)\approx45^{\circ}；当高度为10m时，夹角\theta=\arctan(2)\approx63.4^{\circ}；当高度为15m时，夹角\theta=\arctan(3)\approx71.6^{\circ}。这种角度的变化会影响煤堆表面的风速分布和气流方向，进而影响煤堆的散热和氧气供应。棱台体煤堆模型的上底面长a_1和宽b_1分别设置为5m和3m，下底面长a_2和宽b_2分别设置为10m和6m，高度h同样设置为5m、10m、15m。棱台体的侧面与底面的夹角\alpha和\beta可通过三角函数关系计算得出。例如，对于侧面与长方向底面的夹角\alpha，可根据公式\alpha=\arctan(\frac{h}{\frac{a_2-a_1}{2}})计算；对于侧面与宽方向底面的夹角\beta，可根据公式\beta=\arctan(\frac{h}{\frac{b_2-b_1}{2}})计算。在不同高度设置下，这些夹角也会发生变化，从而对煤堆内部的气流和温度分布产生影响。煤堆的孔隙率\varepsilon是一个重要的参数，它直接影响煤堆内部的气体流通和热量传递。本研究通过控制压实程度来设定不同的孔隙率，孔隙率范围设定为0.1-0.4，以0.05为步长进行取值。在实际操作中，孔隙率的控制可以通过不同的压实设备和压实工艺来实现。例如，使用重型压路机对煤堆进行多次碾压，可以有效降低煤堆的孔隙率。在模拟中，孔隙率的取值依据相关的实验数据和工程经验。研究表明，当孔隙率较低时，煤堆内部的氧气扩散速度较慢，有利于抑制煤的氧化反应；但当孔隙率过低且风速较大时，煤堆内部的热量难以散发，容易导致自燃。煤的物理化学参数对于模拟结果的准确性至关重要。煤的密度\rho取值为1300kg/m^3，这是根据常见煤种的密度范围以及相关的实验测量数据确定的。不同煤种的密度会有所差异，但在本研究中，选取的这个典型值能够代表大多数实际情况。比热容c_p取值为1260J/(kg\cdotK)，该值反映了煤在吸收或释放热量时温度变化的特性，也是基于大量的实验研究和相关文献资料确定的。导热系数\lambda取值为0.26W/(m\cdotK)，它决定了煤堆内部热量传导的速度，同样是根据实际测量和相关研究成果来确定的。煤的氧化反应动力学参数也是模拟中的关键参数。反应活化能E取值为80000J/mol，指的是煤与氧气发生氧化反应时，分子从常态转变为容易发生化学反应的活跃状态所需要的能量。这个数值是通过对不同煤种的氧化反应进行实验研究，并结合相关的化学反应动力学理论得出的。频率因子A取值为1\times10^7s^{-1}，表示单位时间内反应物分子的有效碰撞次数，它反映了反应进行的难易程度。该值同样是基于实验数据和理论分析确定的。这些参数的取值在一定程度上会受到煤种、煤的变质程度等因素的影响，但在本研究中，选取的这些典型值能够满足对煤堆自燃过程的模拟需求。3.3边界条件与初始条件设置在数值模拟中，边界条件和初始条件的合理设置对于准确模拟煤堆自燃过程至关重要。本研究根据实际情况和相关理论，对边界条件和初始条件进行了如下设置。对于风速边界条件，考虑到实际煤场中可能出现的不同风速情况，设定入口风速v的范围为0-20m/s，以1m/s为步长进行模拟。在模型的入口边界，采用速度入口条件，即给定风速的大小和方向。例如，当模拟风速为5m/s时，在入口边界设置风速大小为5m/s，方向垂直于煤堆的迎风面。在出口边界，采用压力出口条件，设定出口压力为标准大气压，即101325Pa。这是因为在实际情况中，煤堆周围的大气压力接近标准大气压，采用压力出口条件可以更真实地模拟煤堆与外界大气的相互作用。温度边界条件方面，煤堆表面与周围环境进行热交换。根据实际的环境温度，设定煤堆表面的环境温度T_{env}为25℃。这是一个常见的环境温度值，在大多数情况下能够代表煤堆所处的实际环境温度。在模拟过程中，考虑煤堆与环境之间的对流换热和辐射换热。对流换热系数h根据经验公式计算，取值为10W/(m^2\cdotK)。该值反映了煤堆表面与空气之间的对流换热能力，是基于大量的实验研究和工程经验确定的。辐射换热采用斯蒂芬-玻尔兹曼定律进行计算，煤堆表面的发射率\varepsilon_{rad}取值为0.9，该值表示煤堆表面发射和吸收辐射能的能力，也是根据相关的实验数据和理论分析确定的。氧气浓度边界条件上，在煤堆表面，氧气浓度C_{O_2}设定为空气中的氧气含量，即0.21mol/m^3。这是因为煤堆表面直接与空气接触，氧气浓度与空气中的氧气含量相同。随着向煤堆内部深入，氧气浓度逐渐降低，其降低的规律遵循扩散定律。在模拟中，通过设置氧气的扩散系数来描述氧气在煤堆内部的扩散过程，氧气的扩散系数D_{O_2}取值为2\times10^{-5}m^2/s，该值反映了氧气在煤堆内部的扩散能力，是根据实验测量和相关研究确定的。在初始条件设置方面，煤堆的初始温度T_0设置为25℃，这与环境温度相同，符合煤堆在初始堆放时的实际情况。煤堆内部的初始氧气含量与表面相同，均为0.21mol/m^3。这是因为在初始状态下，煤堆内部的氧气尚未与煤发生明显的反应，氧气浓度与表面一致。同时，假设煤堆内部的初始水分含量均匀分布，水分含量w_0取值为5%，该值是根据常见煤种的水分含量范围以及实际煤场中的煤堆水分情况确定的。水分在煤堆自燃过程中会影响煤的氧化反应和热量传递，因此合理设置初始水分含量对于准确模拟煤堆自燃过程非常重要。3.4模型的验证与可靠性分析为了验证所建立的煤堆模型的准确性和可靠性，本研究采用了与实验数据对比的方法。在实验方面，设计并搭建了专门的煤堆自燃实验平台。实验平台主要由实验煤堆装置、风速调节系统、温度监测系统和气体成分监测系统等部分组成。实验煤堆装置采用有机玻璃材质制作，以便于观察煤堆内部的情况。煤堆的尺寸和形状尽量与数值模拟中的模型保持一致，以确保对比的有效性。例如，制作了底面半径为5m的圆锥体煤堆和上底面长5m、宽3m，下底面长10m、宽6m的棱台体煤堆，高度均设置为5m。通过不同的压实工艺来控制煤堆的孔隙率，使其分别达到0.1、0.2、0.3和0.4。风速调节系统采用风机和整流装置，能够精确地调节风速，使其在0-20m/s的范围内稳定运行。通过风速传感器实时监测风速，确保实验过程中风速的准确性。在实验中，分别设置风速为3m/s、6m/s、9m/s、12m/s、15m/s和18m/s，模拟不同风速条件下煤堆的状态。温度监测系统由多个热电偶组成，均匀分布在煤堆内部不同位置，能够实时监测煤堆内部的温度变化。热电偶的测量精度为±0.1℃，能够准确地捕捉煤堆温度的细微变化。在实验过程中，每隔1小时记录一次煤堆内部各测点的温度。气体成分监测系统采用气相色谱仪，用于检测煤堆内部的氧气浓度、一氧化碳浓度等气体成分的变化。气相色谱仪能够精确地分析气体成分，其检测精度可以达到ppm级别。通过对气体成分的监测，可以了解煤堆内部氧化反应的进行程度。在数值模拟方面，使用COMSOLMultiphysics软件对相同条件下的煤堆进行模拟。将实验中设置的煤堆尺寸、孔隙率、风速、温度和气体成分等参数输入到模拟模型中，确保模拟条件与实验条件一致。运行模拟程序，得到煤堆内部的温度分布、气流速度分布和氧气浓度分布等模拟结果。对比实验数据与模拟结果，以验证模型的准确性。在温度对比方面，选取煤堆内部几个关键位置的温度进行对比。例如，对于圆锥体煤堆，选取底面圆心正上方1m、2m和3m处的温度；对于棱台体煤堆，选取上底面中心正下方1m、下底面中心正上方1m以及侧面中点处的温度。以风速为6m/s，孔隙率为0.2的圆锥体煤堆为例，实验测得底面圆心正上方2m处的温度在24小时后为32℃，而模拟结果为31.5℃，两者相对误差为1.56%。通过对多个位置和不同工况下的温度对比，发现模拟结果与实验数据的相对误差大部分在5%以内，说明模型在温度模拟方面具有较高的准确性。在氧气浓度对比方面，同样选取煤堆内部不同位置进行对比。例如，在棱台体煤堆中，选取上底面四个角点下方0.5m处和下底面四个角点上方0.5m处的氧气浓度进行对比。当风速为9m/s，孔隙率为0.3时，实验测得下底面一个角点上方0.5m处的氧气浓度为0.18mol/m³，模拟结果为0.175mol/m³，相对误差为2.78%。通过对多个位置和不同工况下的氧气浓度对比，发现模拟结果与实验数据的相对误差大部分在6%以内，表明模型在氧气浓度模拟方面也具有较好的准确性。此外，还将模拟得到的煤堆内部气流速度分布与实验中通过粒子图像测速（PIV）技术测量得到的气流速度分布进行对比。在不同风速和孔隙率条件下，模拟结果与PIV测量结果在气流速度大小和方向上都具有较好的一致性。例如，在风速为12m/s，孔隙率为0.25的情况下，模拟得到的煤堆迎风面中部的气流速度为8.5m/s，PIV测量结果为8.2m/s，相对误差为3.66%。通过与实验数据的全面对比，验证了所建立的煤堆模型在模拟煤堆内部温度分布、氧气浓度分布和气流速度分布等方面具有较高的准确性和可靠性。这为后续深入研究堆放参数对煤堆压实防治自燃的不适风速影响提供了坚实的基础。同时，本研究还将模拟结果与相关文献中的数据进行了对比分析。在查阅了多篇关于煤堆自燃数值模拟和实验研究的文献后，选取了与本研究工况相近的文献数据进行对比。例如，文献[具体文献]中研究了某特定尺寸和孔隙率的煤堆在不同风速下的温度变化情况，将其与本研究中相同参数设置下的模拟结果进行对比，发现温度变化趋势和数值大小都具有较好的一致性。这进一步证明了本研究模型的可靠性和有效性，能够准确地反映煤堆在不同堆放参数和风速条件下的物理现象。四、堆放参数对不适风速影响的模拟结果与分析4.1孔隙率对不适风速的影响4.1.1不同孔隙率下煤堆温度场和氧气浓度场分布利用数值模拟软件COMSOLMultiphysics，对不同孔隙率下煤堆的温度场和氧气浓度场进行了模拟分析。以圆锥体煤堆模型为例，设定底面半径为5m，高度为10m，风速为8m/s，分别模拟了孔隙率为0.1、0.2、0.3和0.4时煤堆内部的物理场分布情况。当孔隙率为0.1时，从煤堆温度场云图（图1）可以看出，煤堆内部温度分布相对较为均匀，整体温度较低，最高温度出现在煤堆中心偏下位置，约为32℃。这是因为孔隙率较低，煤堆内部的气体流通通道狭窄，氧气进入量少，煤与氧气的接触机会有限，氧化反应速率较慢，产生的热量也较少。同时，由于气体流通不畅，热量难以在煤堆内部积聚，大部分热量通过煤堆表面与外界进行热交换而散失。在氧气浓度场云图（图2）中，煤堆表面的氧气浓度接近空气中的氧气含量，为0.21mol/m^3，随着向煤堆内部深入，氧气浓度迅速降低。在距离煤堆表面约1m处，氧气浓度已降至0.1mol/m^3以下，在煤堆中心位置，氧气浓度接近0。这表明在低孔隙率条件下，氧气在煤堆内部的扩散受到极大限制，难以到达煤堆深部，从而抑制了煤的氧化反应。当孔隙率增加到0.2时，温度场分布发生了明显变化（图3）。煤堆内部温度有所升高，最高温度达到38℃，且高温区域范围有所扩大，不仅在煤堆中心偏下位置温度较高，煤堆中上部的部分区域温度也明显上升。这是因为孔隙率的增大使得煤堆内部的透气性有所改善，氧气进入量增加，煤与氧气的接触更加充分，氧化反应速率加快，产生的热量增多。同时，气体流通能力的增强也使得热量在煤堆内部的传递更加容易，导致高温区域的范围扩大。氧气浓度场方面（图4），煤堆内部的氧气浓度分布相对0.1孔隙率时有所改善，氧气能够向煤堆内部更深处扩散。在距离煤堆表面2m处，氧气浓度仍能保持在0.15mol/m^3左右，在煤堆中心位置，氧气浓度约为0.05mol/m^3。这说明随着孔隙率的增加，氧气在煤堆内部的扩散阻力减小，能够为煤的氧化反应提供更多的氧气，促进氧化反应的进行。孔隙率为0.3时，煤堆温度场（图5）显示最高温度进一步升高至45℃，高温区域几乎覆盖了煤堆的大部分区域，煤堆表面温度也有明显上升。此时，煤堆内部的氧化反应更加剧烈，产生的大量热量使得煤堆整体温度大幅升高。由于孔隙率较大，气体流通顺畅，热量在煤堆内部迅速传递，难以在局部积聚，导致高温区域分布较为广泛。氧气浓度场（图6）表明，煤堆内部的氧气浓度分布较为均匀，在距离煤堆表面3m处，氧气浓度仍能达到0.18mol/m^3，煤堆中心位置氧气浓度为0.1mol/m^3左右。这充分说明在高孔隙率条件下，氧气能够在煤堆内部自由扩散，与煤充分接触，极大地加速了煤的氧化反应。当孔隙率达到0.4时，煤堆温度场（图7）显示最高温度已超过50℃，整个煤堆几乎都处于较高温度状态，表明煤堆内部的氧化反应极为剧烈，已经接近自燃的临界状态。此时，大量的氧气进入煤堆，与煤充分反应，产生的热量远远超过了煤堆的散热能力，导致煤堆温度急剧上升。氧气浓度场（图8）显示煤堆内部氧气浓度在大部分区域都能保持在0.2mol/m^3左右，与煤堆表面的氧气浓度相差无几。这说明在极高孔隙率下，煤堆内部的气体流通几乎不受阻碍，氧气能够迅速扩散到煤堆的各个部位，使得煤堆内部的氧化反应在高氧气浓度条件下快速进行。通过对不同孔隙率下煤堆温度场和氧气浓度场分布的分析可以看出，孔隙率对煤堆内部的氧化反应和热量传递有着显著的影响。随着孔隙率的增加，煤堆内部的氧气供应逐渐充足，氧化反应速率加快，产生的热量增多，同时热量传递也更加容易，导致煤堆温度升高，高温区域范围扩大。当孔隙率达到一定程度时，煤堆内部的氧化反应会变得极为剧烈，容易引发自燃。4.1.2孔隙率与不适风速的定量关系为了深入探究孔隙率与不适风速之间的定量关系，对不同孔隙率和风速条件下的模拟数据进行了详细分析。以棱台体煤堆模型为例，上底面长5m、宽3m，下底面长10m、宽6m，高度为15m，孔隙率分别取0.1、0.15、0.2、0.25、0.3、0.35和0.4，风速在0-20m/s范围内以1m/s为步长进行模拟。通过模拟得到不同孔隙率和风速组合下煤堆内部的最高温度数据，将这些数据进行整理和分析。以孔隙率为横坐标，不适风速（定义为煤堆内部最高温度超过40℃时对应的风速）为纵坐标，绘制散点图（图9）。从散点图中可以看出，随着孔隙率的增加，不适风速呈现出先减小后增大的趋势。为了确定孔隙率与不适风速之间的具体函数关系，采用最小二乘法对散点数据进行拟合。经过多次尝试和分析，发现二次函数能够较好地拟合两者之间的关系，拟合函数为：v=-25.36\varepsilon^2+15.78\varepsilon+6.85，其中v为不适风速（m/s），\varepsilon为孔隙率。该拟合函数的决定系数R^2=0.93，表明拟合效果较好，能够较好地反映孔隙率与不适风速之间的定量关系。当孔隙率较低时，如\varepsilon=0.1，代入拟合函数可得不适风速v=-25.36\times0.1^2+15.78\times0.1+6.85=8.17m/s。这意味着在孔隙率为0.1时，当风速超过8.17m/s，煤堆内部最高温度可能会超过40℃，存在自燃风险。随着孔隙率增加到\varepsilon=0.2，不适风速v=-25.36\times0.2^2+15.78\times0.2+6.85=9.02m/s。此时，虽然孔隙率增大使得煤堆透气性改善，但由于风速的影响，热量散失和积聚的平衡发生变化，不适风速有所增加。当孔隙率继续增大到\varepsilon=0.3，不适风速v=-25.36\times0.3^2+15.78\times0.3+6.85=9.34m/s。然而，当孔隙率进一步增大到\varepsilon=0.4时，不适风速v=-25.36\times0.4^2+15.78\times0.4+6.85=8.83m/s，出现了下降趋势。这是因为孔隙率过大时，煤堆内部气体流通过于顺畅，热量容易被带走，但同时也会带来更多的氧气，加速氧化反应，当风速达到一定程度时，氧化反应产生的热量超过了热量散失的速度，导致煤堆内部温度升高，不适风速降低。通过对孔隙率与不适风速定量关系的分析可知，在煤堆压实防治自燃过程中，孔隙率是一个关键参数。合理控制孔隙率可以有效改变不适风速的范围，从而降低煤堆自燃的风险。例如，在实际工程中，如果预计风速较高，可以适当降低煤堆孔隙率，以提高不适风速的阈值；反之，如果风速较低，可以适当增加孔隙率，以保证煤堆内部的通风散热，防止热量积聚。4.2堆放高度对不适风速的影响4.2.1不同高度煤堆的温度变化和自燃进程为了深入探究堆放高度对煤堆温度变化和自燃进程的影响，以圆锥体煤堆模型为例，设定底面半径为5m，孔隙率为0.25，风速为10m/s，分别对高度为5m、10m、15m的煤堆进行了数值模拟。通过模拟得到了不同高度煤堆在不同时刻的温度场分布情况。对于高度为5m的煤堆，在初始阶段，煤堆整体温度与环境温度相同，为25℃。随着时间的推移，煤堆表面开始发生氧化反应，温度逐渐升高。在第3天，煤堆表面温度最高达到30℃，但煤堆内部大部分区域温度仍在28℃以下，温度升高较为缓慢。这是因为煤堆高度较低，热量更容易通过表面散失到周围环境中，抑制了煤堆内部温度的上升。同时，较低的煤堆使得氧气在内部的扩散距离较短，煤与氧气的反应相对较为均匀，不易形成局部高温区域。当煤堆高度增加到10m时，温度变化趋势明显不同。在第3天，煤堆表面温度最高达到35℃，而煤堆内部中心位置温度也上升到32℃。此时，煤堆内部温度升高速度加快，高温区域范围逐渐扩大。这是由于煤堆高度增加，热量散失相对困难，内部积聚的热量增多，加速了煤的氧化反应。同时，氧气在较高的煤堆内部扩散距离变长，导致煤堆内部不同位置的氧气浓度分布差异增大，使得氧化反应在不同区域的进行程度不同，容易形成局部高温区域。对于高度为15m的煤堆，温度变化更为显著。在第3天，煤堆表面温度最高达到40℃，煤堆内部中心位置温度已经超过35℃，高温区域几乎覆盖了煤堆的大部分区域。随着时间继续推移，煤堆温度持续上升，在第5天，煤堆内部最高温度已经接近50℃，表明煤堆已经处于自燃的边缘。这是因为煤堆高度进一步增加，热量积聚效应更加明显，煤堆内部的氧化反应在高温和充足氧气的条件下迅速进行，产生大量热量，使得煤堆温度急剧上升，加速了自燃进程。通过对不同高度煤堆温度随时间变化曲线（图10）的分析，可以更直观地看出高度对煤堆自燃进程的影响。随着煤堆高度的增加，煤堆内部温度上升速度加快，达到自燃临界温度的时间缩短。例如，高度为5m的煤堆在模拟时间内温度上升较为平缓，未达到自燃临界温度；高度为10m的煤堆温度上升速度适中，在模拟后期接近自燃临界温度；而高度为15m的煤堆温度上升迅速，在较短时间内就达到了接近自燃的温度。这表明煤堆高度是影响煤堆自燃进程的重要因素之一，较高的煤堆更容易发生自燃。4.2.2高度与不适风速的关联分析为了研究堆放高度与不适风速之间的关联关系，对不同高度和风速条件下的模拟数据进行了详细分析。以棱台体煤堆模型为例，上底面长5m、宽3m，下底面长10m、宽6m，孔隙率为0.3，高度分别设置为5m、10m、15m，风速在0-20m/s范围内以1m/s为步长进行模拟。通过模拟得到不同高度和风速组合下煤堆内部的最高温度数据，将这些数据进行整理和分析。以煤堆高度为横坐标，不适风速（定义为煤堆内部最高温度超过45℃时对应的风速）为纵坐标，绘制散点图（图11）。从散点图中可以看出，随着煤堆高度的增加，不适风速呈现出逐渐减小的趋势。采用线性回归分析方法对散点数据进行拟合，得到高度与不适风速之间的线性关系方程为：v=-1.2h+12.5，其中v为不适风速（m/s），h为煤堆高度（m）。该线性回归方程的决定系数R^2=0.88，表明拟合效果较好，能够较好地反映高度与不适风速之间的线性关系。当煤堆高度为5m时，代入线性关系方程可得不适风速v=-1.2×5+12.5=6.5m/s。这意味着在高度为5m时，当风速超过6.5m/s，煤堆内部最高温度可能会超过45℃，存在自燃风险。随着煤堆高度增加到10m，不适风速v=-1.2×10+12.5=0.5m/s。此时，煤堆高度的增加使得煤堆内部热量积聚更容易，即使风速较低，也可能导致煤堆内部温度升高，超过自燃临界温度。当煤堆高度进一步增加到15m时，不适风速v=-1.2×15+12.5=-5.5m/s（实际情况中不适风速不可能为负，此处仅为方程计算结果，说明在该高度下，较低的风速就可能引发煤堆自燃）。这充分说明煤堆高度对不适风速有着显著的影响，高度越高，煤堆对风速的耐受性越低，即使是较低的风速也可能导致煤堆内部温度升高，引发自燃。通过对高度与不适风速关联关系的分析可知，在煤堆堆放和管理过程中，合理控制煤堆高度是降低煤堆自燃风险的重要措施之一。根据不同的风速条件，选择合适的煤堆高度，可以有效提高煤堆的安全性。例如，在风速较大的地区，应适当降低煤堆高度，以提高不适风速的阈值，避免煤堆在大风条件下自燃；反之，在风速较小的地区，可以适当增加煤堆高度，但也需要密切关注煤堆内部温度变化，防止因热量积聚而引发自燃。4.3堆放角度对不适风速的影响4.3.1不同堆放角度煤堆的气流分布和热量传递利用数值模拟软件，对不同堆放角度的煤堆进行模拟，以探究其气流分布和热量传递特性。以棱台体煤堆模型为例，上底面长5m、宽3m，下底面长10m、宽6m，高度为10m，孔隙率为0.3，分别模拟了堆放角度为30°、45°、60°时煤堆在风速为8m/s条件下的情况。当堆放角度为30°时，从煤堆气流速度矢量图（图12）可以看出，气流在煤堆迎风面受到阻挡后，部分气流沿着煤堆表面向上流动，在煤堆顶部形成一个气流汇聚区域，此处气流速度较大。而在煤堆背风面，由于气流的绕流作用，形成了一个气流漩涡区域，气流速度相对较小且方向不稳定。在煤堆内部，气流主要沿着孔隙通道流动，但由于孔隙结构较为曲折，气流速度在不同位置存在较大差异。在靠近煤堆表面的区域，气流速度相对较大，随着向煤堆内部深入，气流速度逐渐减小。在热量传递方面，从煤堆温度场云图（图13）可以看出，煤堆表面温度较高，这是因为表面与外界空气直接接触，受到气流的加热作用明显。热量主要通过传导和对流的方式向煤堆内部传递，但由于气流在煤堆内部的分布不均匀，导致热量传递也不均匀。在气流速度较大的区域，热量能够较快地被带走，温度相对较低；而在气流速度较小的区域，热量积聚，温度较高。在煤堆中心位置，由于气流难以到达，热量传递缓慢，温度相对较高。当堆放角度增加到45°时，气流分布发生了明显变化（图14）。在煤堆迎风面，气流沿着表面向上流动的趋势更加明显，在煤堆顶部形成的气流汇聚区域范围扩大，气流速度进一步增大。在煤堆背风面，气流漩涡区域的范围和强度也有所增加。在煤堆内部，孔隙结构相对30°时更加规则，气流在孔隙内的流动更加顺畅，气流速度分布相对更加均匀。在热量传递方面（图15），煤堆表面温度仍然较高，但由于气流分布的变化，热量传递更加均匀。在煤堆内部，由于气流速度分布的均匀性提高，热量能够更有效地在煤堆内部传递，使得煤堆内部温度分布相对更加均匀，高温区域和低温区域的温差减小。当堆放角度达到60°时，气流分布又呈现出不同的特点（图16）。在煤堆迎风面，气流几乎垂直向上流动，在煤堆顶部形成强烈的气流汇聚区域，气流速度达到最大值。在煤堆背风面，气流漩涡区域的范围进一步扩大，强度增强。在煤堆内部，孔隙结构更加有利于气流的流动，气流速度在整个煤堆内部都相对较大且分布较为均匀。在热量传递方面（图17），由于气流速度大且分布均匀，热量能够迅速地在煤堆内部传递，煤堆内部温度分布非常均匀，几乎不存在明显的高温区域和低温区域差异。然而，由于气流带来的热量过多，且煤堆内部的散热能力有限，导致煤堆整体温度升高。通过对不同堆放角度煤堆的气流分布和热量传递分析可知，堆放角度对煤堆内部的气流和热量传递有着显著影响。随着堆放角度的增加，煤堆表面的气流速度增大，气流在煤堆内部的分布更加均匀，但同时也会导致煤堆整体温度升高，热量积聚更加明显，从而增加了煤堆自燃的风险。4.3.2角度与不适风速的内在联系通过理论分析和模拟结果进一步探讨堆放角度与不适风速之间的内在联系。从理论上讲，堆放角度的变化会改变煤堆表面的气流形态和速度分布。当堆放角度较小时，煤堆表面相对平缓，气流在表面的流动受到的阻力较小，气流速度相对较小，能够带走的热量也相对较少。随着堆放角度的增大，煤堆表面变得更加陡峭，气流在表面流动时受到的阻力增大，气流速度增大，能够带走更多的热量。然而，当气流速度过大时，会带来过多的氧气，加速煤的氧化反应，同时也会使煤堆内部的热量积聚难以散发，导致煤堆温度升高，增加自燃风险。以圆锥体煤堆模型为例，底面半径为5m，高度为15m，孔隙率为0.2，对不同堆放角度（30°、45°、60°）和风速（0-20m/s）条件下煤堆内部的最高温度进行模拟分析。模拟结果显示，当堆放角度为30°时，在较低风速（0-5m/s）范围内，煤堆内部最高温度随着风速的增加而缓慢升高；当风速超过5m/s后，煤堆内部最高温度升高速度加快，在风速为10m/s时，煤堆内部最高温度达到40℃，此时可认为10m/s为该堆放角度下的不适风速。当堆放角度增加到45°时，在风速为0-3m/s范围内，煤堆内部最高温度变化较为缓慢；当风速超过3m/s后，煤堆内部最高温度迅速升高，在风速为8m/s时，煤堆内部最高温度达到40℃，不适风速降低到8m/s。这是因为45°的堆放角度使得煤堆表面气流速度增大，在较低风速下就能为煤堆内部带来较多的氧气，加速氧化反应，同时热量积聚也更快，导致不适风速降低。当堆放角度达到60°时，在风速为0-2m/s范围内，煤堆内部最高温度略有上升；当风速超过2m/s后，煤堆内部最高温度急剧升高，在风速为6m/s时，煤堆内部最高温度就达到了40℃，不适风速进一步降低到6m/s。60°的堆放角度使得煤堆表面气流速度更大，在较小的风速下就能引发煤堆内部的快速氧化和热量积聚，从而降低了不适风速。通过对不同堆放角度下不适风速的模拟分析可知，随着堆放角度的增大，煤堆压实防治自燃的不适风速呈现逐渐降低的趋势。这是因为堆放角度的增大改变了煤堆表面的气流特性，使得煤堆在较低风速下就容易出现氧气供应过多、热量积聚难以散发的情况，从而增加了煤堆自燃的风险。在实际煤堆堆放过程中，应根据预计的风速条件，合理选择堆放角度，以降低煤堆自燃的可能性。五、实验研究与验证5.1实验方案设计5.1.1实验煤样选取为确保实验结果的可靠性和代表性，实验煤样选取自某大型煤矿的实际开采煤。该煤矿的煤种为烟煤，具有典型的煤堆自燃特性。在采样过程中，严格遵循相关的采样标准和规范，采用多点采样法，在煤矿的不同开采区域选取了多个采样点，以保证采集到的煤样能够全面反映该煤矿煤的特性。将采集到的煤样进行混合、破碎和筛分处理，使其粒度分布符合实验要求。经过筛分，得到粒度主要分布在0-5mm、5-10mm、10-20mm三个区间的煤样，各区间煤样的质量占比分别为30%、40%和30%。这样的粒度分布能够较好地模拟实际煤堆中煤的粒度情况，因为在实际煤堆中，煤的粒度往往是不均匀的，存在一定的粒度分布范围。5.1.2实验装置搭建实验装置主要包括实验煤堆系统、风速调节系统、温度监测系统和气体成分监测系统。实验煤堆系统：实验煤堆采用有机玻璃材质制作的模型，以便于观察煤堆内部的情况。模型设计为圆锥体和棱台体两种形状，以模拟实际煤堆的常见形状。圆锥体煤堆模型底面半径设置为0.5m，高度分别设置为0.5m、1m、1.5m，通过调整高度来研究不同高度对煤堆压实防治自燃的不适风速的影响。棱台体煤堆模型上底面长0.5m、宽0.3m，下底面长1m、宽0.6m，高度同样设置为0.5m、1m、1.5m。在搭建煤堆模型时，通过不同的压实工艺来控制煤堆的孔隙率。例如，使用小型压实设备对煤堆进行分层压实，每层压实后测量孔隙率，直至达到预定的孔隙率值。孔隙率分别设置为0.1、0.2、0.3和0.4，以研究孔隙率对不适风速的影响。风速调节系统：风速调节系统采用小型风机和整流装置，能够精确地调节风速。风机的风速调节范围为0-20m/s，可以满足实验对不同风速条件的需求。整流装置用于使风机吹出的气流更加均匀稳定，以确保实验过程中风速的准确性和稳定性。在实验过程中，通过风速传感器实时监测风速，并根据需要调整风机的转速，使风速保持在设定值。风速传感器的测量精度为±0.1m/s，能够准确地测量风速的变化。温度监测系统：温度监测系统由多个热电偶组成，均匀分布在煤堆内部不同位置。热电偶的测量精度为±0.1℃，能够准确地捕捉煤堆温度的细微变化。在圆锥体煤堆模型中，热电偶布置在底面圆心正上方0.1m、0.3m、0.5m（对应高度为0.5m的煤堆），以及高度方向上每隔0.2m的位置（对应高度为1m和1.5m的煤堆）。在棱台体煤堆模型中，热电偶布置在上底面中心正下方0.1m、下底面中心正上方0.1m以及侧面中点处等关键位置。通过数据采集系统实时采集热电偶测量的温度数据，并将数据传输到计算机进行分析处理。气体成分监测系统：气体成分监测系统采用便携式气相色谱仪，用于检测煤堆内部的氧气浓度、一氧化碳浓度等气体成分的变化。气相色谱仪能够精确地分析气体成分，其检测精度可以达到ppm级别。在实验过程中，定期从煤堆内部不同位置采集气体样本，通过气相色谱仪进行分析，以了解煤堆内部氧化反应的进行程度。例如，在煤堆不同高度和径向位置设置气体采样点，每隔一定时间采集一次气体样本，分析其中氧气和一氧化碳的浓度变化。5.1.3变量控制与测量方法在实验过程中，严格控制各个变量，以确保实验结果的准确性和可靠性。实验的自变量包括煤堆形状（圆锥体、棱台体）、高度（0.5m、1m、1.5m）、孔隙率（0.1、0.2、0.3、0.4）和风速（0-20m/s）。通过精确设置实验装置的参数来控制这些自变量。例如，在搭建煤堆模型时，使用精确的测量工具确保煤堆的形状和尺寸符合预定要求；通过不同的压实工艺和测量手段控制孔隙率；利用风速调节系统和风速传感器准确控制和测量风速。因变量为煤堆内部的温度和气体成分（氧气浓度、一氧化碳浓度）。采用上述的温度监测系统和气体成分监测系统对因变量进行测量。在测量温度时，每隔1小时记录一次各热电偶测量的温度数据。对于气体成分的测量，每隔3小时采集一次气体样本进行分析。在实验过程中，保持环境温度和湿度相对稳定，环境温度控制在25±2℃，湿度控制在50±5%。通过空调和加湿器等设备来调节环境温度和湿度，以减少环境因素对实验结果的影响。同时，对实验过程中的各种数据进行详细记录，包括实验时间、实验条件（煤堆形状、高度、孔隙率、风速等）、测量数据（温度、气体成分等）。在每次实验结束后，对实验数据进行整理和初步分析，及时发现数据中的异常情况，并采取相应的措施进行处理。例如，如果发现某个热电偶测量的温度数据异常偏高或偏低，检查热电偶是否损坏或安装位置是否正确，必要时重新进行测量。5.2实验过程与数据采集在完成实验方案设计与装置搭建后，严格按照实验流程开展实验，以获取准确且具有代表性的数据。实验开始前，再次检查实验装置的各个部分，确保风速调节系统、温度监测系统和气体成分监测系统均正常运行，各传感器的校准准确无误。同时，对实验环境进行全面检查，确保环境温度和湿度稳定在设定范围内，避免环境因素对实验结果产生干扰。将选取好的煤样按照预定的孔隙率和形状要求进行堆放。对于圆锥体煤堆，使用模具辅助成型，确保底面半径和高度符合设计尺寸。在堆放过程中，采用分层压实的方法，每层铺设一定厚度的煤样后，使用小型压实设备进行压实，通过测量压实后煤样的体积和质量，计算并调整孔隙率，直至达到预定的孔隙率值。例如，在构建孔隙率为0.2的圆锥体煤堆时，先将煤样分层铺设，每层厚度控制在10-15cm，每次压实后使用孔隙率测量仪进行测量，根据测量结果调整压实力度和铺设厚度，最终达到目标孔隙率。棱台体煤堆的搭建同样遵循类似的方法，确保上底面和下底面的尺寸以及高度准确，通过精确控制压实工艺实现不同孔隙率的设置。完成煤堆搭建后，启动风速调节系统，将风速设置为初始值0m/s，保持1小时，以便煤堆内部的温度和气体成分达到稳定状态。随后，以1m/s为步长逐渐增加风速，每调整一次风速，稳定运行2小时，以确保煤堆内部的物理过程充分发展并达到相对稳定状态。在每个风速条件下稳定运行期间，按照预定的数据采集频率进行数据采集。温度数据采集方面，温度监测系统中的热电偶将实时测量的温度数据传输至数据采集系统。数据采集系统每隔10分钟记录一次各热电偶测量的温度值。例如，在圆锥体煤堆中，位于底面圆心正上方0.1m处的热电偶每10分钟将测量的温度数据传输并记录，同时记录该时刻的实验时间、风速等相关信息。对于棱台体煤堆，布置在上底面中心正下方0.1m、下底面中心正上方0.1m以及侧面中点处的热电偶同样按照此频率进行数据采集。气体成分数据采集时，每隔30分钟从煤堆内部不同位置采集气体样本。使用气体采样器在预定的采样点抽取气体，将采集到的气体样本迅速注入便携式气相色谱仪进行分析。例如，在煤堆高度方向上每隔0.3m设置一个气体采样点，在径向方向上每隔0.2m设置一个采样点，按照预定时间间隔进行采样分析。气相色谱仪能够精确分析气体样本中的氧气浓度、一氧化碳浓度等成分，并将分析结果传输至计算机进行记录和存储。在整个实验过程中，安排专人负责记录实验数据和观察实验现象。记录的数据包括实验时间、风速、各热电偶测量的温度、气相色谱仪分析的气体成分以及实验过程中出现的任何异常情况。例如，若发现某个热电偶的测量数据出现异常波动，立即检查热电偶的连接是否松动或损坏，同时对该位置附近的煤堆状态进行观察，记录是否存在明显的气体泄漏或其他异常现象。实验人员还需密切关注煤堆表面是否有冒烟、温度异常升高等可能预示自燃的现象，一旦发现异常，及时采取相应措施，并详细记录异常发生的时间、风速条件以及煤堆的具体状态。当风速达到20m/s并完成数据采集后，逐渐降低风速至0m/s，关闭风速调节系统。继续监测煤堆内部的温度和气体成分变化1小时，以观察煤堆在停止通风后的恢复情况。最后，对实验装置进行清理和维护，为下一次实验做好准备。通过严格按照上述实验过程进行操作和数据采集，确保了实验数据的准确性和可靠性，为后续的实验结果分析和与模拟结果的对比提供了坚实的基础。5.3实验结果与模拟结果对比分析将实验测得的不同堆放参数下煤堆压实防治自燃的不适风速数据与数值模拟结果进行对比，结果如表1所示。以圆锥体煤堆为例，当孔隙率为0.2，高度为10m时，实验测得的不适风速为9.5m/s，而模拟结果为9.2m/s，相对误差为3.16%。在棱台体煤堆中，当孔隙率为0.3，高度为15m时，实验得到的不适风速为7.8m/s，模拟结果为7.5m/s，相对误差为3.85%。煤堆形状孔隙率高度(m)实验不适风速(m/s)模拟不适风速(m/s)相对误差(%)圆锥体0.2109.59.23.16圆锥体0.3158.27.93.66棱台体0.3157.87.53.85棱台体0.4106.56.24.62从整体对比情况来看，实验结果与模拟结果具有较好的一致性，相对误差大部分在5%以内。这表明数值模拟模型能够较为准确地预测不同堆放参数下煤堆压实防治自燃的不适风速，验证了模拟模型的可靠性和有效性。然而，在部分数据对比中仍存在一定的差异。例如，在圆锥体煤堆孔隙率为0.3，高度为15m的情况下，虽然相对误差在可接受范围内，但模拟结果略低于实验结果。分析其原因，一方面可能是实验过程中存在一些难以完全控制的因素，如煤样的不均匀性、实验装置的微小误差以及环境因素的细微变化等。尽管在实验过程中采取了一系列措施来控制这些因素，但仍无法完全消除其影响。另一方面，数值模拟模型虽然考虑了多种物理过程和参数，但在实际情况中，煤堆内部的物理化学过程可能更加复杂，模型可能无法完全准确地描述所有细节。例如，煤堆内部的水分迁移、化学反应的复杂性以及煤颗粒之间的相互作用等，都可能对煤堆的自燃过程产生影响，而模拟模型在这些方面的描述可能存在一定的简化。通过对实验结果与模拟结果的对比分析，进一步验证了数值模拟在研究煤堆压实防治自燃不适风速方面的可行性和准确性。同时，也明确了实验和模拟过程中存在的不足之处，为后续进一步优化模拟模型和改进实验方法提供了方向。在未来的研究中，可以进一步提高实验的精度和控制水平，减少实验误差。例如，采用更加先进的实验设备和技术，对煤样进行更精细的处理和分析，以提高煤样的均匀性和实验数据的准确性。在模拟模型方面，可以进一步完善模型的物理过程描述，考虑更多的影响因素，如煤堆内部的水分迁移、化学反应动力学的复杂性等，以提高模型的准确性和可靠性。六、基于模拟与实验结果的煤堆自燃防治策略6.1根据堆放参数优化煤堆压实方案基于上述模拟与实验结果，针对不同的堆放参数，可制定如下优化的煤堆压实方案。孔隙率方面：根据模拟得到的孔隙率与不适风速的定量关系，在实际煤堆堆放中，应根据预计的风速情况合理控制孔隙率。当预计风速较低时，可以适当增加煤堆的孔隙率，以保证煤堆内部有良好的通风散热条件。例如，在风速常年低于5m/s的地区，可将煤堆孔隙率控制在0.3-0.4之间。这样的孔隙率能够使煤堆内部的氧气供应相对充足，同时又能保证热量及时散发，避免因热量积聚导致自燃。而当预计风速较高时，应降低煤堆孔隙率，以提高不适风速的阈值。比如在风速经常超过10m/s的区域，将煤堆孔隙率控制在0.1-0.2之间。较低的孔隙率可以减少氧气进入量，降低煤的氧化反应速率，即使在大风条件下，也能有效抑制煤堆自燃。在实际操作中，可以通过不同的压实设备和压实工艺来控制孔隙率。使用重型压路机进行多次碾压，可使煤堆孔隙率降低；而采用较为松散的堆放方式，则可适当提高孔隙率。堆放高度方面：由模拟结果可知，煤堆高度与不适风速呈负相关关系，即煤堆高度越高，不适风速越低，煤堆越容易自燃。因此，在煤堆堆放过程中，应严格控制煤堆高度。对于风速较大的地区，应显著降低煤堆高度。在风速超过15m/s的情况下，煤堆高度不宜超过5m。较低的煤堆高度可以使热量更容易散发到周围环境中，减少热量积聚，降低自燃风险。而在风速较小的区域，煤堆高度可以适当增加，但也不宜过高。比如在风速低于5m/s的地方，煤堆高度可控制在10m左右。同时，在煤堆高度较高时，可以采取分层压实的方法，在每层煤之间铺设一定厚度的惰性材料（如沙子、粉煤灰等），以进一步抑制氧气的进入和热量的传递，降低自燃可能性。堆放角度方面：模拟结果表明，随着堆放角度的增大，不适风速逐渐降低，煤堆自燃风险增加。在实际煤堆堆放时，应根据风速条件合理选择堆放角度。当风速较大时，应选择较小的堆放角度。在风速超过10m/s的情况下，堆放角度不宜超过30°。较小的堆放角度可以使煤堆表面的气流速度相对较小，减少氧气进入量，同时也有利于热量的散发。而当风速较小时，可以适当增大堆放角度，但一般也不宜超过45°。比如在风速低于5m/s的情况下，堆放角度可选择40°左右。这样既能保证煤堆有一定的稳定性，又能在一定程度上利用气流带走热量，降低自燃风险。此外，在煤堆表面可以铺设一些防风材料（如防风网、土工布等），以进一步减小风速对煤堆的影响，特别是对于堆放角度较大的煤堆，防风材料的作用更为明显。6.2考虑风速因素的煤堆储存管理措施根据模拟与实验结果，在不同风速条件下，应采取不同的煤堆储存和管理措施。低风速条件下（风速小于5m/s）：在低风速环境中，煤堆内部的热量散失相对较慢，因此重点应放在加强通风散热和合理控制氧气供应上。可以适当增大煤堆的孔隙率，以促进空气流通，提高散热效率。在堆放煤堆时，可以采用较为松散的堆放方式，使煤堆孔隙率保持在0.3-0.4之间。这样既能保证有足够的氧气供应，维持煤堆内部一定的氧化反应速率，又能通过空气的流通及时带走氧化反应产生的热量，避免热量积聚导致自燃。同时，定期对煤堆进行翻动，使煤堆内部的热量能够均匀散发。例如，每隔一周对煤堆进行一次翻动，将煤堆表层的煤翻到内部，内部的煤翻到表层，促进热量交换。在煤堆表面设置一些通风管道或通风孔，进一步增强通风效果。通风管道可以采用直径为10-15cm的PVC管，按照一定的间距（如1-2m）均匀分布在煤堆表面，深入煤堆内部一定深度（如0.5-1m），使空气能够直接进入煤堆内部，加速热量的散失。中