Fissure evolution and variation of pyrolysis kinetics parameters of tar-rich coal during heat treatment under different atmosphere
-
摘要:
富油煤的清洁高效开发有助于保障我国的能源安全。选择柠条塔矿富油煤,开展了N2和欠氧气氛条件下的加热试验,利用体视显微镜获取不同温度处理后富油煤的表观形貌,在氩气气氛中,对富油煤在0~600 ℃下热重(TG)、微商热重(DTG)、差示扫描量热(DSC)进行测试,基于激光闪射法监测富油煤热导率变化,并分析了不同升温速率下活化能、频率因子等热解动力学参数的变化规律,探讨富油煤热处理后裂隙演化和热解动力学参数变化。研究结果表明:随着加热温度的升高,富油煤表面裂隙数目明显增多,逐渐贯通,且镜煤裂隙发育程度更高。相比N2气氛,同温度下欠氧气氛中煤样裂隙率更高。质量损失率变化趋势与裂隙率相似,由于欠氧环境中煤样发生氧化反应,导致其增速更快。N2气氛中富油煤质量损失率与裂隙率迅速提升的温度点明显滞后于欠氧气氛。室温至350 ℃范围内,由于热膨胀影响,富油煤热导率随着温度升高有所增加。根据热重分析,富油煤随温度变化可分为室温至300 ℃、300~500 ℃和500~600 ℃ 3个阶段,在阶段Ⅰ和Ⅱ分别由于吸附水析出和有机质的分解吸热,频率因子和活化能逐渐增大;在阶段Ⅲ由于无机矿物的分解放热,频率因子和活化能降低。研究结果有助于理解富油煤高温开发过程中的结构演化机制和传热机理。
Abstract:The clean and efficient development of tar-rich coal is helpful in ensuring China's energy security. The tar-rich coal from Ningtiaota Mine was selected and heated under N2 and hypoxia atmosphere. Using a stereomicroscope to obtain the apparent morphology of tar-rich coal treated at different temperatures, the thermal gravimetric (TG), differential scanning calorimetry (DSC), and differential scanning calorimetry (DTG) of tar-rich coal were tested at 0-600 ℃ in an argon atmosphere. Based on the laser flash method, the thermal conductivity changes of tar-rich coal were monitored, and the changes in pyrolysis kinetics parameters such as activation energy and frequency factor under different heating rates were analyzed. The crack evolution and pyrolysis kinetics parameter variation of tar-rich coal after heat treatment were explored. The results show that with the increase of heating temperature, the number of cracks on the surface of the tar-rich coal sample increases significantly and gradually. Additionally, the degree of fracture development in vitrain is higher. Compared with the N2 atmosphere, the crack rate of coal samples in the hypoxia atmosphere at the same temperature is higher. The changing trend of mass loss ratio is similar to the trend of fracture rate, and the growth rate is faster due to the oxidation reaction of coal samples in the hypoxia environment. The temperature at which the mass loss rate and fracture rate of tar-rich coal rapidly increase in the N2 atmosphere is clearly lagging behind that in the hypoxia atmosphere. In the range of room temperature to 350 ℃, due to the effect of thermal expansion, the thermal conductivity of tar-rich coal increases with the increase of temperature. According to TG analysis, the tar-rich coal can be divided into three stages with temperature changes: 25-300 ℃, 300-500 ℃ and 500-600 ℃. In stages I and II, due to the precipitation of adsorbed water and the heat absorption of organic matter decomposition, the frequency factor and activation energy gradually increase; In stage III, due to the exothermic decomposition of inorganic minerals, the frequency factor and activation energy decrease. The results are help understand the structural evolution mechanism and heat transfer mechanism of tar-rich coal in the process of high temperature development.
-
Keywords:
- tar-rich coal /
- high temperature /
- cracks /
- thermal analysis /
- pyrolysis kinetics
-
0. 引言
煤矿智能化开采的3要素包括智能感知、智能决策和智能控制[1]。其中,智能感知是实现智能化开采的第1步,也是智能决策和智能控制的基础和重要依据。带式输送机作为煤矿生产过程中的重要运输设备[2],其智能化监测是实现煤矿智能化开采的必要环节[3]。目前,常用的带式输送机故障监测装置为电子类开关量传感器[4-5],缺乏对故障演化过程的实时定量监控,不利于数据的深度融合与分析,难以及时指导参数调整和停机维修[6]。同时,传统的电子类传感器在复杂煤矿环境中的抗干扰能力和远距离传输能力不足,难以满足实时定量监测对信号传输的高精度要求[7]。因此,寻求一种安全可靠、可实时定量测量的带式输送机故障监测方式非常必要。
方崇全[8]提出利用巡检机器人收集带式输送机运行状态音、视频信息并分析健康状态,降低了人工巡检成本和工作量,但单点巡检的方式在故障实时全面监测方面仍然存在一定的局限性。汪楚森[9]与韩涛等[10]利用机器视觉技术对带式输送机进行故障监测,具有成本小、非接触式测量等优势,但视觉监测的精度易受光线影响,且煤尘污染镜头可能会直接导致监测功能失效。马宏伟等[11]采用红外热成像技术监测带式输送机的关键部件温度,克服了传统点式温度测量法的监测区域小、布线复杂等难题,但该技术亦存在抗干扰能力弱、分辨率低等缺点。而光纤型传感器可实现对应变、温度等物理量的精确测量[12],相较于传统传感器具有:体积小、抗电磁干扰、灵敏度高、耐腐蚀、可定量测量、本质安全等优点[13],在航空、医疗、机械生产等领域已有广泛应用。近年来,光纤传感器在煤矿领域的应用也逐渐得到了发展。NARUSE等[14]构建了一个基于布里渊光时域反射原理的分布式光纤地下矿井应变传感系统,现场实验证明该系统能够有效监测煤矿大规模开采区域周围应力分布不均而导致的回风巷道变形。柴敬等[15]研发的钻孔植入式光纤光栅传感系统,为煤岩变形的监测提供了一种有效方法。梁敏富等[16]提出了一种感知锚杆杆体轴力的光纤光栅传感器,并比较测试了影响应变传递效率的因素,发现黏结剂种类对应变传递效率影响最为明显。李丽君等[17]提出了一种光纤光栅矿压传感器,可以埋入煤、岩层用于压力测量。张登攀等[18]针对煤层瓦斯压力不能多点同步测量的问题,提出一种基于弹性膜片的光纤光栅瓦斯压力传感器,能够有效避免安装误差和迟滞带来的影响,同时便于组链开展多点分布式测量。梁敏富等[19]基于光纤光栅传感原理及弹性膜片结构,设计了一种新型的压力传感器,实现了温度补偿的压力测量并成功用于锚杆支护质量监测。孙丽等[20]以连通器为基本原理,设计了一种基于等强度悬臂梁的新型双光纤光栅静力水准仪,用于解决建筑结构不均匀沉降监测问题。于涛等[21]针对煤矿巷道顶板离层监测设备可靠性低、准确度差及无法实现实时监测等问题,设计了一种基于等强度悬臂梁的光纤光栅顶板离层仪,具有灵敏度高、量程大、安装方便等优点。此外,光纤传感技术在液压支架阻力监测[22]、刮板机直线度感知[23]、设备结构件变形测量[24]等方面也有广泛应用。上述研究利用光纤传感器的应变敏感特性,实现了对煤矿环境参数和设备状态的监测,为煤矿的安全高效生产提供了技术支持。然而,带式输送机作为煤矿生产的主要运输设备,利用光纤传感技术对其进行故障监测的研究却鲜有报道。
为了寻求安全可靠、可定量测量的带式输送机故障监测与保护方案,笔者提出了一种基于FBG的带式输送机综合保护系统,明确了光纤型传感器阵列的设计是系统搭建的重要前提。在分析输送机常见故障特征的基础上,设计了多种光纤型传感器并组成了带式输送机综合保护系统的传感器阵列。作为光纤型传感器阵列的核心敏感元件,笔者对FBG–等强度悬臂梁的结构尺寸设计和材料选择进行了深入研究,借助理论分析、数值模拟及实验室试验等手段对其结构有效性、灵敏度和稳定性进行了验证。所设计的光纤型传感器阵列解决了现有带式输送机运行故障无法实时定量感知的问题,为实现设备智能化感知与控制提供了坚实的数据基础。
1. 带式输送机结构及故障分析
1.1 带式输送机结构分析
在众多带式输送机类型中,托辊式带式运输机凭借其良好的稳定性、适应性和承载能力成为煤矿运输领域应用最广的运输设备[25],故本文以托辊式带式输送机为故障监测对象(下文简称带式输送机)。
如图1所示,带式输送机主要由5部分构成:输送带、驱动装置、滚筒、托辊及机身机架。其中,输送带承载并运输物料;驱动装置(电机与减速器)提供动力;滚筒分为主动滚筒和改向滚筒,前者连接驱动装置驱动输送带,后者改变输送带方向;托辊数量众多,支撑并导向输送带以保证其稳定运行;机架为输送机各部分提供支撑。此外,还包括头尾轮、导向轮、张紧装置和联结装置等辅助部件,共同构成完整的带式输送机系统。
1.2 常见故障分析
带式输送机在运行过程中,常会出现诸如跑偏、打滑、堆料、过热和撕带等故障。
1)跑偏故障。跑偏故障是指输送带偏离期望轨道的现象,其主要原因是张紧装置异常导致的输送带两侧牵引力不均,从而发生跑偏;其次,输送带的老化、损坏、连接处松动,滚筒和托辊磨损严重以及机架变形引起的倾斜也可能引发跑偏故障。
2)打滑故障。打滑故障是指输送带与驱动滚筒之间因摩擦力不足发生打滑而无法正常传递动力的现象。张紧装置调节不当或故障,会导致输送带张力不足,无法与主动滚筒产生足够摩擦力;此外,输送带卡阻,主动滚筒变形和包胶材料磨损老化导致的摩擦因数减小也会诱发打滑故障。
3)堆料故障。堆料故障是指物料在输送带或输送路径上堆积的现象。常见原因包括:驱动装置异常引起的输送速度与上料速度不匹配;张紧装置异常会导致输送带松弛,物料在输送过程中易滑落堆积;此外,滚筒和托辊的磨损、变形、卡滞等问题,也会影响输送机正常运行,造成带上物料分布不均,增加堆料风险。
4)过热故障。过热故障是指带式输送机组成部件温度异常升高的现象。输送带磨损老化、托辊卡滞和滚筒打滑等可能导致输送带摩擦过热严重,进而引发输送带或物料起火;前述跑偏和堆料故障的未及时纠正会加剧摩擦升温;此外,驱动装置故障或过载,也是温度异常的重要原因。
5)撕带故障。撕带故障是指输送带出现的撕裂或破损现象。输送带的磨损与老化会直接导致其发生撕裂或断裂;滚筒、托辊表面磨损严重,大块物料落入或卡入会施加异常摩擦力或冲击力于输送带上,从而引发撕带故障。此外张紧装置异常,导致输送带张力偏离正常范围,也会增加撕带风险。
带式输送机运行故障与部件健康状态的关系,如图2所示。这些故障不仅能反映组成部件的损坏与性能劣化,更是衡量整个输送机健康状态的重要指标[26]。通过对这些故障的实时定量监测与综合分析,能够帮助工作人员及时确定故障部件,把握输送机的整体健康状态,提高设备运维效率,确保高效生产和设备稳定运行。
2. 带式输送机综保系统传感器阵列
2.1 FBG带式输送机综合保护系统
传统的带式输送机综合保护系统的监测形式主要表现为开关量信号的采集,这种简单的监测形式缺乏对故障演化过程的实时定量监控,在数据融合与分析方面存在明显局限性,难以支撑现代工业对智能化、精细化保护系统的迫切需求。为了推动带式输送机向更智能化的方向迈进,亟需引入实时、定量的故障监测方式,为设备的潜在故障提前预警、故障根源精准定位和自调节提供数据基础。
针对上述问题,笔者提出了一种基于FBG的带式输送机综合保护系统,如图3所示。该系统主要由光纤型传感器阵列、光纤数据传输系统、数据处理与分析单元、报警与预警系统、人机交互界面和故障应急控制系统六大功能部分组成。
光纤型传感器阵列实时监测带式输送机的运行状态,运行故障会触发相应光纤型传感器动作并产生承载故障信息的光信号;借助光纤数据传输系统,故障信息得以在复杂的煤矿环境中,以传输安全、抗干扰能力强的光信号形式传输至主机室;经光信号解调主机解调得到的故障信息,一方面被传输到地面控制室存储显示,方便工作人员随时了解输送机健康状态;另一方面被传输至数据处理与分析单元,经报警(预警)逻辑分析后,生成控制指令并传递给报警与预警系统和故障应急控制系统;报警与预警系统接受控制信号后,立即发出声光报警提醒工作人员故障发生;故障应急控制系统接受控制信号后,自动执行诸如调整设备运行参数或停机等操作[27],最小化故障对设备和生产的影响。
作为系统运行、决策、动作的依据,故障监测传感器是带式输送机综合保护系统的核心单元之一。相较传统监测系统,所提出的基于FBG的带式输送机综合保护系统的优势均得益于其采用了光纤传感技术对带式输送机进行故障监测。然而,尽管光纤传感技术在工业监测领域已经展现出了巨大潜力,但目前其在带式输送机故障监测领域的研究仍处于初级阶段。因此,光纤型传感器阵列的设计是构建光纤型带式输送机综合保护系统的重要前提。
2.2 光纤型传感器阵列
光纤型传感器凭借其实时定量测量的特性,可为搭建自动化、智能化的带式输送机综合保护系统提供数据基础,从而实现报警条件自定义、故障提前预警和健康状态综合分析等功能;此外,在煤矿井下这类复杂环境中,光纤型传感器具有相较传统电子类传感器更强的抗电磁干扰和远距离传输能力,更适用于有实时定量要求的精准测量。
在分析带式输送机故障原理的基础上,基于光纤光栅的应变敏感特性设计了以FBG–等强度悬臂梁为核心敏感元件的光纤型跑偏传感器、光纤型堆料传感器、光纤型撕带传感器和光纤型拉绳开关;基于光纤光栅的温度敏感特性设计了光纤型温度传感器。这些传感器协同工作,组成了带式输送机综合保护系统的光纤型传感器阵列。
2.2.1 应变敏感类故障传感器
应变敏感类故障传感器以FBG–等强度悬臂梁为核心敏感元件,通过相应的传感结构将故障监测量转换为等强度悬臂梁的应变;粘贴在等强度悬臂梁表面的光纤光栅由于应变敏感特性,内部反射光中心波长发生偏移;经解调仪解调和温度补偿后得到由单一应变引起的中心波长偏移量,最后根据相应的结构动作关系建立的传感模型即可解算出故障监测量。该类传感器的温度补偿措施通常为内置一个温补光纤光栅,粘贴在无应变的敏感基材上(文中为尼龙6)或搭配柔性导热材料使用,在避免应变干扰的同时,又具有良好的导热性,以达到理想的温度补偿效果。
1)跑偏传感器。跑偏传感器被成对安装在输送带张力变化较大区域的输送带两侧,尤其是输送机的机头和机尾处,能够实时监测到输送带的跑偏情况,便于系统及时采取应对措施。如图4所示,跑偏传感器主要由跑偏传动机构、凸轮、压力传动机构、等强度悬臂梁、光纤光栅及外壳等组成。其中,跑偏传动机构由拨杆、连接柄、扭力弹簧、转轴组成;压力传动机构由滑轮、滑轮轴、弹簧板1、压力弹簧、限位筒、弹簧板2组成。
带式输送机发生跑偏时,输送带推动拨杆绕转轴摆动,进一步通过跑偏传动机构驱使凸轮旋转;压力传动机构将凸轮旋转角转换为施加在等强度悬臂梁上的作用力,最后使粘贴在等强度悬臂梁表面的光纤光栅产生应变。利用光纤光栅的应变敏感特性和传感模型分析,即可得到输送带的跑偏量。
2)堆料传感器。堆料传感器被安装在输送机上料点、下料点和两条输送机重叠点处的输送带中线上方,安装高度据平均物料流量和堆料裕量设定。如图5a所示,堆料传感器主要由拨杆、球阀、压力传动机构、等强度悬臂梁、光纤光栅及外壳等组成。其中压力传动机构由上顶板、下顶板、压力弹簧组成。测量过程如图5b和图5c所示。
当带式输送机发生堆料故障时,堆积的物料推动拨杆,进一步驱使球阀转动;压力传动机构将球阀转动角转换为施加在等强度悬臂梁上的作用力,最后使粘贴在等强度悬臂梁表面的光纤光栅产生应变。利用光纤光栅的应变敏感特性和传感模型分析,即可得到堆料高度。
3)撕带传感器。撕带传感器被安装在输送带上料点向物流方向上的10~15 m处,且介于上下输送带之间。通常撕带故障一旦出现即需立刻停机,而利用光纤光栅实现定量测量的意义在于可根据现场情况自定义报警阈值,有效减少因物料飞溅误触发而对正常生产造成的妨碍。此外,撕带传感器还能够实现对物料飞溅异常增多的感知,为故障的提前预警分析提供数据基础。如图6所示,撕带传感器主要由载料板、侧挡板、下底壳、限位柱、承重簧、压力传动机构、等强度悬臂梁和光纤光栅等组成。
撕带故障的测量过程包括3个阶段:预加载与稳定支撑阶段,承重簧首先执行第1压缩行程,用于支撑传感器自身质量和正常洒落的物料。撕裂测量阶段,当输送带发生纵向撕裂故障或物料飞溅异常增多时,超过一定质量的物料落在载料板上,承重簧进入第2压缩行程;同时,压力传动机构也开始工作,产生压力作用在等强度悬臂梁的自由端;最后使粘贴在等强度悬臂梁表面的光纤光栅产生应变;利用光纤光栅的应变敏感特性和传感模型分析,即可得到物料飞溅情况和撕带情况。过载保护阶段,载料板上方物料超出承重上限,承重簧和压力传动结构停止动作,限位柱接管承重任务避免传感器过载损坏。
4)拉绳开关。工作人员发现异常情况,可通过拉绳开关手动控制停机。为确保在输送机沿线均能及时控制停机,拉绳开关被间隔安装在输送机的两侧机架上和机头、机尾处,两两之间用钢绳连接,钢绳中间用托环支撑。拉绳开关的定量化增加了触发条件的可自定义功能,降低了对连接钢绳的安装要求;同时,通过合理设置触发动作范围和触发时间阈值,可减少意外接触引起的误触发,避免妨碍正常生产。如图7所示,堆料传感器与跑偏传动机构相似,主要由拉力传动机构、凸轮、压力传动机构、等强度悬臂梁、光纤光栅及外壳等组成。其中,拉力传动机构由拉杆、转轴、扭力弹簧组成;压力传动机构由滑轮、滑轮轴、弹簧板1、压力弹簧、限位筒、弹簧板2组成。
工作人员发现故障需要紧急停机带式输送机时,可拉动沿线安装的拉绳开关;拉力传动机构驱动凸轮旋转,凸轮旋转角经压力传动机构转换为施加在等强度悬臂梁上的作用力,最后使粘贴在等强度悬臂梁表面的光纤光栅产生应变。利用光纤光栅的应变敏感特性和传感模型分析,即可感知到拉绳动作。
2.2.2 温度传感器
温度传感器被铺设在温度异常高风险的区域内,用于监测过热故障,如输送机机头和机尾的易高温区域、沿线机架上、滚筒和托辊附近等。FBG温度传感器可以实现温度的准分布式测量,同时,根据所在通道数和原始中心波长可以实现温度异常点的定位功能。如图8所示,温度传感器主要由光缆、光纤光栅、金属导热管、光缆接头等组成,金属导热管套设在光纤光栅栅区外,方便导热以及防止光纤光栅受应变影响。
FBG温度传感器被按需铺设在输送机上,测温点间隔可以通过调整光纤光栅栅区间隔按需定制;同一光纤通道,支持不同中心波长的多个光纤光栅串联同时测量;光纤光栅的反射光中心波长受温度影响产生偏移;经解调仪解调后,可以得到由温度变化产生的中心波长偏移量。根据光纤光栅的温度敏感特性,可进一步解算出温度情况。
所设计的光纤型传感器阵列考虑了如下措施提高传感精度和灵敏度:① 光纤光栅具有一定的栅区长度,采用等强度悬臂梁而非普通悬臂梁,使测量区的应变接近相同值;② 考虑黏胶长度与应变传递效率的关系[28-29],选用光纤光栅的栅区长度为1 cm,对应黏胶长度2 cm;③ 考虑温度补偿,通过设置温补光纤,消除温度影响。
2.3 磁敏式打滑传感器
目前,光纤型敏感元件在转速测量方面的实现难度颇高,相较之下,磁敏式转速传感器凭借其灵敏度高、抗干扰能力强、热稳定性好以及非接触式测量等优点,在低转速测量领域得到了广泛应用[30]。如图9所示,将磁敏式转速传感器作为打滑故障监测的解决方案。此类传感器以霍尔元件或磁敏电阻为敏感元件,当被测物体旋转带动磁场变化时,磁性敏感元件内部的电子受洛伦兹力作用,会在导体中形成电动势或电流。电信号和方向与磁场的强度和变化率有关,通过对该电信号的测量即可实现对物体转速的感知,进而为打滑故障的监测提供可靠依据。
打滑故障程度可以通过比较目标驱动速度和实际运行速度衡量。目标驱动速度可以通过测量驱动滚筒的转速,结合其半径计算得到;而由于张紧滚筒、托辊的转动完全由输送带摩擦驱动,故输送带的实际运行速度可以通过测量张紧滚筒、托辊的转速,并结合其半径计算得到。
3. FBG–等强度悬臂梁
在所设计的光纤型传感器阵列中,包括凸轮结构、压力弹簧和扭力弹簧等经典机械元件已经能够实现精确设计与加工,且动作传递关系明确。而等强度悬臂梁虽在传感器领域已经有应用,但其在光纤传感领域的适用性以及如何实现预期传感效果的研究尚不成熟,主要体现在结构尺寸设计和材料选择两个方面。下文以跑偏传感器的测量要求和传感模型为例,对FBG–等强度悬臂梁的结构设计和材料选择进行深入研究,并分析传感精度和灵敏度的影响因素。研究成果同样适用于光纤型传感器阵列中的其他传感器。
3.1 FBG感知及补偿原理
光纤光栅常用于测量温度和应变2个基本物理量,其基本感知理论[16]如下:
$$ \frac{{\Delta \lambda }}{\lambda } = (1 - {P_{\mathrm{e}}})\varepsilon + {K_{\mathrm{T}}}\Delta T $$ (1) 式中:$\lambda $,$\Delta \lambda $,${P_{\mathrm{e}}}$,$\varepsilon $,${K_{\mathrm{T}}}$,$\Delta T$分别为光纤光栅的初始中心波长,波长偏移量,弹光系数(${P_{\mathrm{e}}} = 0.22$),应变,温度灵敏度和温度变化量。
当光纤光栅工作在无应变条件下时($\varepsilon = 0$),反射光中心波长与温度的变化关系如式(2)所示;用光纤光栅测量应变时,只需在相同测量环境中再设置一温补光纤光栅,即可实现温度补偿,光纤光栅中心波长与应变关系如式(3)所示:
$$ \frac{{\Delta {\lambda _{\rm{T}}}}}{{{\lambda _{\rm{T}}}}} = {K_{\rm{T}}}\Delta T $$ (2) $$ \frac{{\Delta {\lambda _{\rm{p}}}}}{{{\lambda _{\rm{p}}}}} - \frac{{\Delta {\lambda _{\rm{T}}}}}{{{\lambda _{\rm{T}}}}} = (1 - {P_{\rm{e}}})\varepsilon $$ (3) 式中:${\lambda _{\rm{T}}}$,$\Delta {\lambda _{\rm{T}}}$,${\lambda _{\text{p}}}$,$\Delta {\lambda _{\text{p}}}$,$\Delta {T_{\rm{T}}}$分别为温度敏感光纤光栅的初始中心波长和波长偏移量,应变敏感光纤光栅的初始中心波长和中心波长偏移量,温度变化量。
根据式(2)和式(3),可分别实现对温度和应变的测量。
3.2 尺寸和结构设计
FBG–等强度悬臂梁设计的关键在于:根据传感器的一般工况,控制其受力和应变范围,在不超过光纤光栅应变测量范围前提下,尽可能控制应变范围更大,以弥补传感器加工精度的不足,提高传感精度。跑偏传感器的测量要求为跑偏角变化范围$({0^ \circ },{60^ \circ })$,对应作用力变化范围$(0,100\;{\mathrm{N}})$,在此前提下,对FBG–等强度悬臂梁进行尺寸和结构设计。
根据式(3),初始中心波长${\lambda _{\rm{p}}}$的确定由最大中心波长偏移量$\Delta {\lambda _{{\rm{p}}\max }}$和最大应变${\varepsilon _{\max }}$决定。$\Delta {\lambda _{{\rm{p}}\max }}$和${\varepsilon _{\max }}$的增大,能提高FBG–等强度悬臂梁及传感器的精度和灵敏度。当应变$\varepsilon \leqslant $
2000 ×10−6,波长变化量$\Delta {\lambda _{\rm{p}}} \leqslant 2\;{\text{nm}}$时,光纤光栅的测量稳定性最佳[31-32],为了达到最高精度和灵敏度,理想状态下的参数选择:${\varepsilon _{\max }} = $2000 ×10−6,$\Delta {\lambda _{{\rm{p}}\max }} = 2\;{\text{nm}}$,${\lambda _{\rm{p}}} = 1\;282\;{\text{nm}}$。根据光纤光栅解调仪常见解调范围,调整参数为${\lambda _{\rm{p}}} = 1\;530\;{\text{nm}}$,$\Delta {\lambda _{{\rm{p}}\max }} = 2\;{\text{nm}}$,${\varepsilon _{\max }} = $1675 ×10−6。再结合跑偏传感器的一般工况,确定FBG–等强度悬臂梁的设计目标为:作用力F的变化范围(0,100 N),对应应变$\varepsilon $的变化范围(0,1675 ×10−6)。所设计的等强度悬臂梁外形如图10a所示,其中,h为截面高度;截面宽度$b = b(x),x \in (0,L)$,x为截面到固定端的距离;L为悬臂梁的最大长度。
应力$ \sigma $与弯矩$M(x)$的关系为
$$ \sigma = \frac{{M(x)}}{{W(x)}} $$ (4) 等强度悬臂梁允许的最大应力${\sigma _{\max }}$等于所用材料的许用应力$\left[ \sigma \right]$,即:
$$ \frac{{M(x)}}{{W(x)}} = \left[ \sigma \right] $$ (5) 弯矩$M(x)$为
$$ M(x) = F(L - x) $$ (6) 抗弯截面系数$W(x)$为
$$ W(x) = \frac{{b(x){h^2}}}{6} $$ (7) 联立式(5)、式(6)和式(7),可得截面宽度函数$b(x)$为
$$ b(x) = \frac{{6{F_0}(L - x)}}{{{h^2}\left[ \sigma \right]}} $$ (8) 式中:${F_0}$为预设力,即设计允许的自由端的最大受力。
分析式(8)可知,等强度悬臂梁的截面宽度$b(x)$是关于x且斜率为负的一次函数。特别地,当$x = L$时,$b(x) = 0$,显然不能满足剪切强度要求,因此需要根据剪切强度要求确定自由端附近的最小截面宽度${b_{\min }}$。
悬臂梁的切应力强度条件为
$$ {\tau _{\max }} = \frac{{3{F_{\rm{s}}}}}{{2bh}} \leqslant \left[ \tau \right] $$ (9) 其中,$\left[ \tau \right]$为许用切应力,${F_{\rm{s}}}$为剪力。悬臂梁的剪力为
$$ {F_{\rm{s}}}(x) = F $$ (10) 联立式(9)和式(10)可得:
$$ {b_{\min }} = \frac{{3{F_0}}}{{2\left[ \tau \right]h}} $$ (11) 且当$x = L - \dfrac{{h\left[ \sigma \right]}}{{4\left[ \tau \right]}}$时,$b(x) = {b_{\min }}$。
综上,FBG–等强度悬臂梁的尺寸设计依据如下:
$$ b(x) = \frac{{6{F_0}(L - x)}}{{{h^2}\left[ \sigma \right]}} $$ (12) $$ {b_{\min }} = \frac{{3{F_0}}}{{2\left[ \tau \right]h}}\text{,}x = {L_0} = L - \frac{{h\left[ \sigma \right]}}{{4\left[ \tau \right]}} $$ (13) $$ {b_{\max }} = \frac{{6{F_0}L}}{{{h^2}\left[ \sigma \right]}}\text{,}x = 0 $$ (14) 分析式(12)、式(13)和式(14)可知,等强度悬臂梁的尺寸设计影响因素包括:① 材料特性。许用应力$\left[ \sigma \right]$和许用切应力$\left[ \tau \right]$;② 预设力${F_0}$。其中,材料特性难以控制,应将材料选择作为设计考虑的主要因素;预设力${F_0}$的设定相对灵活,将其作为设计考虑的辅助因素,补偿材料特性引起的尺寸和应变范围的不合理。
传感器中的悬臂梁常选用弹性材料或弹塑性材料,对于弹塑性材料,以屈服极限${\sigma _{\rm{s}}}$为基 准,许用应力$\left[ \sigma \right]$和许用切应力$\left[ \tau \right]$分别为
$$ \left[\sigma \right]=\frac{{\sigma }_{{\mathrm{s}}}}{{n}_{1}}\text{,}{n}_{1}\in (1.5,2.5) $$ (15) $$ \left[ \tau \right] = \left[ \sigma \right]{n_2},{n_2} \in (0.6,1.0) $$ (16) 经后续的理论分析和有限元仿真发现,当${n_1} = 2.5$,${n_2} = 0.6$,选用材料尼龙6,并将预设力调整到${F_0} = 1\;000\;{\text{N}}$时,FBG–等强度悬臂梁的尺寸和应变表现最好。尼龙6的屈服极限${\sigma _{\rm{s}}} = 1.036 \times {10^{ - 8}}\;{\text{N/}}{{\text{m}}^{\text{2}}}$,弹性模量$E = 2.62 \times {10^9}\;{\text{N/}}{{\text{m}}^{\text{2}}}$。代入计算,得到等强度悬臂梁的最合理尺寸:截面高度$h = 12.06\;{\text{mm}}$,总长$L = 30\;{\text{mm}}$,最大截面宽度${b_{\max }} = 24.86\;{\text{mm}}$,当${L_0} = 24.975\;{\text{mm}}$时取到最小截面宽度$ {b_{\min }} = 5\;{\text{mm}} $。
为了方便受力、固定和粘贴光纤光栅,本文优化了等强度悬臂梁结构,分别设置了受力区域、固定区域和光纤光栅尾纤槽,如图10b所示。
3.3 应变传递模型
FBG–等强度悬臂梁在外力作用下能够产生均匀应变并进行测量,这是光纤型传感器阵列工作的重要一环,其转换关系值得深入研究。
应变$ \varepsilon $和悬臂梁应力$\sigma $的计算公式分别为
$$ \varepsilon = \frac{\sigma }{E} $$ (17) $$ \sigma = \frac{{Mh}}{{2{I_{\rm{z}}}}} $$ (18) ${I_{\rm{z}}}$为悬臂梁的惯性矩:
$$ {I_{\rm{z}}} = \frac{{{h^3}}}{{12}}b(x) $$ (19) 联立式(17)、式(18)和式(19)可得:
$$ \varepsilon = \frac{{6M(x)}}{{E{h^2}b(x)}} $$ (20) 联立式(8)、式(20)可得应变$\varepsilon $与力$F$的传递关系:
$$ \varepsilon = \frac{{F\left[ \sigma \right]}}{{{F_0}E}} $$ (21) 式(21)中,$\left[ \sigma \right]$、${F_0}$、$E$均为常量,因此,所设计的等强度悬臂梁的应变$\varepsilon $随力$F$变化而均匀变化。
同时,式(21)也提供了目标工况下等强度悬臂梁的尺寸设计思路:可以通过材料的$\dfrac{{\left[ \sigma \right]}}{E}$来初步调节应变$\varepsilon $与作用力$F$的变化关系;在确定材料后,可以通过调整预设力${F_0}$来进一步调节应变$\varepsilon $达到合理范围。
3.4 传感模型分析
为了探究FBG–等强度悬臂梁的相关设计因素对传感器精度和灵敏度的影响,根据设计的跑偏传感器建立了跑偏量传感模型,并基于此模型进行分析。跑偏传感器的监测量与实际测得的中心波长偏移量的映射关系可概括为:跑偏量−角度−力−应变−中心波长偏移量。带式输送机和传感器安装的相对位置决定跑偏量–角度的转换关系;力与应变的转换关系如式(21)所示;应变–波长偏移量的转换关系如式(3)所示;角度–力的转换关系如下。
如图11所示,跑偏传感器的压力传动机构内部主要为压力弹簧,当凸轮转动带动滑轮杆上下移动时,弹簧压缩并产生弹力作用在等强度悬臂梁的自由端。
通过合适的凸轮外形设计和弹簧选择,可以将转动角度$\theta $与压力$F$的对应关系控制为
$$ F = K\theta $$ (22) 式中:K为传动系数,根据跑偏传感器的一般工况:$ \theta \in ({0\text{°} },{60\text{°} }) \to F \in (0,100\;{\text{N}}) $,得到$K = 1.67\;{\text{N}}{{\text{/}}({\text{°} })}$。
综合式(3)、式(12)、式(21)和式(22),得到跑偏传感模型为
$$ F = \frac{{{F_0}E}}{{(1 - {P_{\rm{e}}})\left[ \sigma \right]}}\left(\frac{{\Delta {\lambda _{\rm{p}}}}}{{{\lambda _{\rm{p}}}}} - \frac{{\Delta {\lambda _{\rm{T}}}}}{{{\lambda _{\rm{T}}}}}\right) $$ (23) $$ \theta = \frac{{{F_0}E}}{{(1 - {P_{\rm{e}}})\left[ \sigma \right]K}}\left(\frac{{\Delta {\lambda _{\rm{p}}}}}{{{\lambda _{\rm{p}}}}} - \frac{{\Delta {\lambda _{\rm{T}}}}}{{{\lambda _{\rm{T}}}}}\right) $$ (24) 跑偏传感器的传感灵敏度S1和FBG–等强度悬臂梁的灵敏度S2分别为
$$ {S _1} = \frac{{{F_0}E}}{{(1 - {P_{\rm{e}}})\left[ \sigma \right]{\lambda _{\rm{p}}}K}} $$ (25) $$ {S _2} = \frac{{{F_0}E}}{{(1 - {P_{\rm{e}}})\left[ \sigma \right]{\lambda _{\rm{p}}}}} $$ (26) $$ {S _2} = K{S _1} $$ (27) 由式(25)可知,跑偏传感器的灵敏度与材料特性$\dfrac{E}{{\left[ \sigma \right]}}$、预设力${F_0}$正相关,与初始中心波长${\lambda _{\rm{p}}}$、传动系数K反相关;其中$\dfrac{E}{{\left[ \sigma \right]}}$、${F_0}$、${\lambda _{\rm{p}}}$是FBG–等强度悬臂梁的主要设计因素。在设计传感器时,应该协调各参数关系,满足尺寸合理性,并达到目标应变范围。尼龙6的$E = 2.62 \times {10^9}\;{\text{N/}}{{\text{m}}^{\text{2}}}$,$\left[ \sigma \right] = 4.144 \times {10^7}\;{\text{N/}}{{\text{m}}^{\text{2}}}$,代入计算得到跑偏传感器的理论灵敏度${S _1} = {31.723\;3{(\text{°} )}}{\text{/nm}}$,FBG–等强度悬臂梁的理论灵敏度${S _2} = 52.978\;0\;{\text{N/nm}}$。在跑偏传感器最大量程$\theta = {60{\text{°}}}$时,对应波长变化量$\Delta {\lambda _{\rm{p}}} \approx 1.89\;{\text{nm}}$(目标$\Delta {\lambda _{\rm{p}}} = 2\;{\text{nm}}$),符合设计要求。
4. 数值分析及有限元仿真
4.1 材料分析
测试分析了多种材料设计的等强度悬臂梁,其中金属材料包括:304不锈钢、易切削黄铜、3003–铝合金、1345铝合金等;塑料材料包括:尼龙101、尼龙6/10、尼龙6等,通过计算和ANSYS有限元仿真得到各自最合理尺寸及最大量程应变,结果见表1。
表 1 材料选择与尺寸设计Table 1. Table of material selection and dimension design材料 ${F_0}$/N ${F_{\max }}$/N 最优尺寸 最大量程应变 评价 bmax/mm bmin/mm h/mm ${L_0}$/mm 304不锈钢 100 100 93 2 1.5 29.4 $4.5{ \times }{10^{ - 4}}$ 尺寸不合理 易切削黄铜 100 100 55.14 2 2.52 28.9 $4.0{ \times }{10^{ - 4}}$ 尺寸不合理 3003–铝合金 100 100 39.5 3 5.1 27.9 $2.55{ \times }{10^{ - 4}}$ 应变过小 1345铝合金 100 100 25.6 3 7.5 26.9 $1.8{ \times }{10^{ - 4}}$ 应变过小 尼龙101 100 100 25.6 2 5.2 27.8 $2.3{ \times }{10^{ - 2}}$ 应变过大 尼龙6/10 100 100 62.05 2 2.25 29.1 $6.3{ \times }{10^{ - 3}}$ 尺寸不合理 尼龙6 100 100 45.73 2 3 28.7 $1.5{ \times }{10^{ - 2}}$ 应变较大 由表1可知,金属材料的屈服极限和弹性模量较大,制作的等强度悬臂梁在尺寸合理性和传感精度方面均表现不佳。以304不锈钢制作的等强度悬臂梁为例,最大截面宽度为93 mm,需要在29.4 mm的距离内过渡到2 mm的最小截面宽度,尺寸表现显然不合理;同时,其最大量程应变均值仅有450×10−6(目标
1675 ×10−6),传感精度也表现不佳。使用铝合金1345虽然能够获得合理的设计尺寸,但同样其最大量程应变仅有180×10−6,难以满足传感精度要求。多数尼龙类材料虽然能够获得合理的设计尺寸,但其弹性模量较小,最大量程应变过大,超出了光纤光栅的应变测量范围,因此可以根据式(21)调整预设力${F_0}$,从而获得更合理的应变范围和外形参数。材料尼龙6具有强度高、耐磨和耐腐蚀等优点,常被用于制造各种机械零件。综合考虑尺寸和应变范围的合理性以及材料的易获得性等因素,确定尼龙6作为等强度悬臂梁的材料。最终设计参数见表2。
表 2 最终材料选择与尺寸设计Table 2. Table of final material selection and dimension design材料 ${F_0}$/N ${F_{\max }}$/N 最优尺寸 最大量程应变 评价 bmax/mm bmin/mm h/mm ${L_0}$/mm 尼龙6 1 000 100 24.8 5 12.1 24.9 ${\text{2}}{{.2\times }}{10^{ - 3}}$ 尺寸、应变均合理 4.2 有限元仿真与分析
在确定等强度悬臂梁的尺寸设计和材料选择后,利用有限元分析软件ANSYS对等强度悬臂梁作进一步的力学分析。网格单元尺寸设置为:0.2 mm。划分网格后,生成单元513 164个,节点301 221个。对等强度悬臂梁的固定区域施加固定载荷,模拟现实中对悬臂梁的固定约束。设置作用力变化范围$ (0,100\;{\text{N}}) $。从0开始,按照10 N的梯度加载到100 N,并记录应变表现。加载区域为受力端圆环区域,加载方向垂直于圆环区域表面。仿真分析得到等强度悬臂梁的应变情况如图12、图13所示,测量区的应变均值随压力变化情况如图14所示。
由图12可知,所设计的等强度悬臂梁在上表面中轴线附近处具有良好的等值应变表现,对于提升光纤光栅的栅区测量精度意义重大。由图13进一步得知,在确定的压力作用下,等强度悬臂梁在上表面中线区域的应变波动较小,即在误差可接受的情况下,所设计的等强度悬臂梁应变与距离无关,满足光纤光栅的栅区应变近似相等的设计目标。由图14可知,所设计的等强度悬臂梁应变均值与自由端受力F的大小成线性正相关关系,与理论推导式(21)相符。
5. 性能测试与结果分析
5.1 传感器测试系统及装置
通过实验室试验进一步对FBG–等强度悬臂梁的结构有效性、灵敏度和稳定性进行了验证。测试系统由图15中的光纤光栅、等强度悬臂梁、夹具、光纤光栅解调仪、客户端、应变片、应变仪和图16中的万能试验机等组成。光纤光栅粘贴在等强度悬臂梁上表面中轴线位置,连接光纤光栅解调仪,测量应变引起的反射光中心波长偏移量;客户端连接光纤光栅解调仪,显示和记录波长变化情况;应变片粘贴在等强度悬臂梁下表面,连接应变仪,用于测量实际应变并与波长变化量进行实时对应;等强度悬臂梁固定在夹具上,方便万能试验机加压,实验室试验如图16所示。
5.2 试验及结果分析
试验采用逐步加载、分级记录的试验形式。加载增速为$1\;{\text{N/s}}$,加载终点设置为105 N。考虑万能试验机的加载特点,以$F = 2.4\;{\text{N}}$作为初始点,此时中心波长作为初始波长,应变仪调零。记录作用力在0 ~100 N内的波长和应变的变化情况,记录间隔5 N。
5.2.1 灵敏度及重复性试验
在室温10~12 ℃的试验环境下,对试制的3个FBG–等强度悬臂梁试样(试样1,2,3)分别进行压力应变试验,得到各试样在不同压力作用下,等强度悬臂梁受力与应变变化曲线,如图17所示;等强度悬臂梁受力与光纤光栅中心波长偏移量变化曲线,如图18所示。选取灵敏度表现较好的试样1,再进行重复性加载试验,得到等强度悬臂梁受力与应变重复性试验变化曲线,如图19所示。
据图17可知,所设计的等强度悬臂梁应变随压力增大而线性增大。3个试样均能够有效地将压力转换为应变,且转换关系表现出良好的线性关系。受到最大作用力100 N时,3个试样的实际应变均在
2000 ×10−6附近,符合设计目标(目标2000 ×10−6)。由图18可知,光纤光栅波长偏移量随着压力增大而线性增大,利用光纤光栅这点特性可以测量等强度悬臂梁的应变,即所设计的FBG–等强度悬臂梁结构能够实现力–应变–光纤光栅波长偏移量的转换,且拟合结果(表3)表现出良好的线性关系。试验所得3个试样的平均灵敏度为$S _2^* = 38.115\;7\;{\text{N/nm}}$,代入跑偏传感模型得到实际的平均灵敏度为$S _1^* = {22.823\;8({\text{°}} )}{\text{/nm}}$。此外,受最大作用力100 N时,3个试样产生的光纤光栅波长偏移量均在2.0~2.5 nm,符合设计要求(目标2 nm)。据图19可知,所设计的FBG–等强度悬臂梁试样在重复加载的情况下,3次试验中力–应变表现基本重合,平均重复性误差仅为1.002%,表现出良好的重复性和稳定性。表 3 等强度悬臂梁受力与光纤光栅中心波长偏移量数据拟合Table 3. Table of data fitting between the applied force on an equal-strength cantilever beam and the central wavelength shift of fiber bragg gratings试样 拟合函数 拟合系数 传感灵敏度/(N·nm−1) 平均灵敏度$S_2^*$/(N·nm−1) 试样1 y=0.023 3x−0.060 7 0.998 2 42.918 5 38.115 7 试样2 y=0.028 0x−0.103 1 0.998 4 35.714 3 试样3 y=0.028 0x−0.041 9 0.998 1 35.714 3 试验测得的跑偏传感器的灵敏度$S _1^* = {22.823\;8{(\text{°} )}}{\text{/nm}}$,FBG–等强度悬臂梁的灵敏度$S _2^* = 38.115\;7\;{\text{N/nm}}$,相比理论灵敏度平均灵敏度$S _1^{} = {31.723\;3{(\text{°} )}}{\text{/nm}}$,$S _2^{} = 52.978\;0\;{\text{N/nm}}$偏小,分析误差原因:① 胶层应变传递效率和光纤光栅应变测量效率难以达到100%,这是误差形成的主要原因;② 光纤光栅的手工粘贴位置和胶层情况难以控制;③ 等强度悬臂梁尺寸计算设计及加工有误差。但在加工传感器时,需要根据实际表现调整标定关系,可以接受小范围误差的存在。
5.2.2 温度敏感性及温度补偿试验
1)温度敏感性试验
在不加载的条件下,对试样1和试样2的光纤光栅进行温度补偿敏感性试验,测试试样在不同温度下的中心波长变化情况,温度测试间隔5 ℃,测试范围10~50 ℃。试验结果如图20所示。
对实验结果进行线性拟合分析,试样1的线性拟合系数为0.997 8,试样2的线性拟合系数为0.998 3,可见光纤光栅对温度变化敏感,且中心波长随着温度变化而线性变化。这表明:用光纤光栅测量应变时,需要进行温度补偿以消除温度影响,提高传感精度;此外,还可利用该特性设计温度传感器,实时测量带式输送机相关部件的温度。
2)温度补偿试验
在3个试样中选择压力应变试验中灵敏度表现最佳的试样1,将其波长变化量代入式(3)得到应变测量计算值,并参考应变仪读出的实际应变,对比其在温度补偿前后的表现,试验结果如图21所示。
据图21可知,在温度补偿前后,试样1通过光纤光栅测量得到的应变都随压力变化而线性变化,与应变片测量得到的实际应变趋势表现基本一致。而实际应变与温度补偿前后通过光纤光栅测量得到的应变存在一定差异,分析其主要原因为胶层应变传递效率和光纤光栅应变测量效率难以达到100%,在加工传感器时,尤其要考虑该误差因素,并根据实际情况进行标定。
温度补偿试验数据拟合情况见表4,相比温度补偿前FBG–等强度悬臂梁的应变测量灵敏度
19.5212 ×10−6 N−1,温度补偿后的应变测量灵敏度为19.6413 ×10−6 N−1,更接近实际应变拟合灵敏度21.5961 ×10−6 N−1。可见,即便在微小温度变化环境下(10~12 ℃),补偿机制仍能有效提高灵敏度0.6%。因此,在实际应用过程中,尤其是在温度变化复杂的环境下,需要考虑温度补偿。表 4 温度补偿试验数据拟合Table 4. Table of data fitting for temperature compensation experiments试验内容 拟合函数 拟合系数 灵敏度/(10−6 N−1) 实际应变 y=21.596 1x−52.406 9 0.998 6 21.596 1 温补前测量应变 y=19.521 2x−50.864 5 0.998 3 19.521 2 温补后测量应变 y=19.641 3x−49.575 7 0.998 4 19.641 3 6. 结 论
1)为了寻求更加安全有效的带式输送机监测与保护方案,本文分析了带式输送机常见故障原理与表现形式,基于光纤传感技术,提出了一种基于FBG的带式输送机综合保护系统,明确了光纤型传感器阵列的设计是系统搭建的重要前提。以FBG–等强度悬臂梁为核心传感结构设计了多种传感器,组成了带式输送机综合保护系统的传感器阵列。
2)对FBG–等强度悬臂梁的结构设计和材料选择进行了深入研究,优选了尼龙6作为等强度悬臂梁的制作材料。通过理论分析、数值模拟及实验室试验等手段对FBG–等强度悬臂梁的结构有效性、灵敏度和稳定性进行了验证。理论分析和试验表明:FBG–等强度悬臂梁的理论灵敏度为$S = 52.978\;0\;{\text{N/nm}}$,实际灵敏度为$S_{}^* = 38.115\;7\;{\text{N/nm}}$;光纤型传感器对温度变化敏感,通过温度补偿能够提高传感精度。所设计的光纤型传感器阵列解决了现有带式输送机运行故障无法实时定量感知的问题,为实现设备智能化感知与控制提供了坚实的数据基础。
3)以视觉、激光、红外热成像等新兴技术为代表的非接触式测量具有适用性强、测量范围广等优点,但易受环境波动的影响、测量精度有限。包括所设计的光纤型传感器阵列在内的接触式测量,虽然测量精度高,但在高强度恶劣工作环境中难免出现精度衰减或损坏的情况。鉴于此,未来应该着眼于多种传感器的综合应用。具体而言,对于监测范围广、数量多且精度要求不高的部件采用非接触式测量。而对于精度要求高,且测量点位规律或集中的部件,采用接触式测量。从而,有效提高故障监测效率,降低设备故障率,推动煤矿行业向更高水平的自动化、智能化方向迈进。
-
表 1 富油煤分析测试结果
Table 1 Characterization of tar-rich coal sample
工业分析/% 元素分析/% 发热量Qgr,ad/(MJ·kg−1) 焦渣特征 Mad Ad Vdaf FCdaf Cdaf Hdaf Odaf Ndaf St,d 3.70 6.15 40.71 59.29 81.35 5.00 12.14 1.00 0.49 29.80 2 表 2 欠氧气氛热处理后煤样表面裂隙发育
Table 2 Development of surface cracks in coal samples after heat treatment under air atmosphere
温度 100 ℃ 150 ℃ 200 ℃ 225 ℃ 显微图 二值化 温度 250 ℃ 275 ℃ 300 ℃ 325 ℃ 显微图 二值化 表 3 N2气氛热处理后煤样表面裂隙发育
Table 3 Development of surface cracks in coal samples after heat treatment under N2 atmosphere
温度 200 ℃ 250 ℃ 300 ℃ 350 ℃ 显微图 二值化 温度 400 ℃ 500 ℃ 550 ℃ 600 ℃ 显微图 二值化 -
[1] 王双明,孙 强,乔军伟,等. 论煤炭绿色开采的地质保障[J]. 煤炭学报,2020,45(1):8−15. WANG Shuangming,SUN Qiang,QIAO Junwei,et al. Geological guarantee of coal green mining[J]. Journal of China Coal Society,2020,45(1):8−15.
[2] 王双明,师庆民,王生全,等. 富油煤的油气资源属性与绿色低碳开发[J]. 煤炭学报,2021,46(5):1365−1377. WANG Shuangming,SHI Qingmin,WANG Shengquan,et al. Resource property and exploitation concepts with green and low-carbon of tar-rich coal as coal-based oil and gas[J]. Journal of China Coal Society,2021,46(5):1365−1377.
[3] 薛香玉,王长安,邓 磊,等. 基于全生命周期的富油煤原位热解碳排放[J]. 煤炭学报,2023,48(4):1773−1781. XUE Xiangyu,WANG Changan,DENG Lei,et al. Carbon emissions from in-situ pyrolysis of tar-rich coal based on full life cycle analysis method[J]. Journal of China Coal Society,2023,48(4):1773−1781.
[4] 于海洋,许永彬,陈智明,等. 双碳目标下煤炭深部流态化开采及前景[J]. 洁净煤技术,2023,29(1):15−32. YU Haiyang,XU Yongbin,CHEN Zhiming,et al. Deep fluidized coal mining and its prospect under the target of carbon peak and carbon neutralization[J]. Clean Coal Technology,2023,29(1):15−32.
[5] 王双明,王 虹,任世华,等. 西部地区富油煤开发利用潜力分析和技术体系构想[J]. 中国工程科学,2022,24(3):49−57. doi: 10.15302/J-SSCAE-2022.03.006 WANG Shuangming,WANG Hong,REN Shihua,et al. Potential analysis and technical conception of exploitation and utilization of tar-rich coal in western China[J]. Strategic Study of CAE,2022,24(3):49−57. doi: 10.15302/J-SSCAE-2022.03.006
[6] 申艳军,王 旭,师庆民,等. 榆神府矿区富油煤煤相及孔隙结构特征试验研究[J]. 煤矿安全,2021,52(10):30−37,44. SHEN Yanjun,WANG Xu,SHI Qingmin,et al. Coal facies characteristics and pore structure response of oil-rich coal in Yushenfu Mining Area[J]. Safety in Coal Mines,2021,52(10):30−37,44.
[7] 张 旭,王利斌,裴贤丰,等. 煤热解提高焦油产率及品质关键技术与研究进展[J]. 煤炭科学技术,2019,47(3):227−233. ZHANG Xu,WANG Libin,PEI Xianfeng,et al. Research progress and key technology of improving coal tar yield and quality by coal pyrolysis[J]. Coal Science and Technology,2019,47(3):227−233.
[8] HAN Lei,SHEN Jian,WANG Jinyue,et al. Characteristics of pore evolution and its maceral contributions in the huolinhe lignite during coal pyrolysis[J]. Natural Resources Research,2021,30:2195−2210. doi: 10.1007/s11053-020-09792-3
[9] MENG Dexi,YUE Chengyan,WANG Tian,et al. Evolution of carbon structure and functional group during Shenmu lump coal pyrolysis[J]. Fuel,2021,287:119538. doi: 10.1016/j.fuel.2020.119538
[10] 胡 鑫,王超勇,孙 强,等. 核磁共振技术对热解煤孔隙结构分形特征研究[J]. 矿业研究与开发,2021,41(11):67−75. HU Xin,WANG Chaoyong,SUN Qiang,et al. Study on fractal characteristics of pore structure of pyrolysis coal by nuclear magnetic resonance technology[J]. Mining Research and Development,2021,41(11):67−75.
[11] 张 蕾,韩智坤,舒 浩,等. 陕北富油煤低温热解提油基础特性[J]. 煤炭工程,2022,54(9):124−128. ZHANG Lei,HAN Zhikun,SHU Hao,et al. Basic characteristics of tar extraction in low temperature pyrolysis of tar-rich coal from Northen Shaanxi[J]. Coal Engineering,2022,54(9):124−128.
[12] JU Yang,ZHU Yan,ZHOU Hongwei,et al. Microwave pyrolysis and its applications to the in situ recovery and conversion of oil from tar-rich coal:an overview on fundamentals,methods,and challenges[J]. Energy Reports,2021,7:523−536. doi: 10.1016/j.egyr.2021.01.021
[13] 单鹏飞,杨 攀,来兴平,等. 微波–水交互作用下富油煤岩渐进性破坏规律试[P]. 中国: CN202211277067.0,2023-03-28. [14] JU Yang,ZHU Yan,ZHANG Yuwei,et al. Effects of high-power microwave irradiation on tar-rich coal for realising in situ pyrolysis,fragmentation,and low-carbon utilisation of tar-rich coal[J]. International Journal of Rock Mechanics and Mining Sciences, 2022, 157: 105165.
[15] 王双明,石智军,王生全,等. 一种过热水蒸汽原位热解富油煤高效提油方法[P]. 中国:CN110541695A,2019−12−06. [16] 李利峰. 褐煤热损伤过程中的孔裂隙细观特性研究[J]. 煤矿安全,2020,51(3):12−15. LI Lifeng. Study on mesoscopic characteristics of pores and cracks during thermal damage of lignite[J]. Safety in Coal Mines,2020,51(3):12−15.
[17] 李海琪. 过热水蒸汽热解长焰煤的孔裂隙结构演化规律研究[D]. 太原:太原理工大学,2021. LI Haiqi. Study on the evolution law of pore and fissure structure of long flame coal after superheated steam pyrolysis[D]. Taiyuan:Taiyuan University of Technology,2021.
[18] 毛崎森,王长安,侯育杰,等. 富油煤原位热解对流加热过程传热规律数值模拟[J]. 洁净煤技术,2023,29(8):48−58. MAO Qisen,WANG Changan,HOU Yujie,et al. Numerical simulation of waste heat utilization process after in-situ pyrolysis of tar-rich coal[J]. Clean Coal Technology,2023,29(8):48−58.
[19] 董光顺,朱超凡,厉家宗,等. 黄陵矿区富油煤对流加热原位转化开发效果数值模拟[J]. 煤田地质与勘探,2023,51(4):57−67. DONG Guangshun,ZHU Chaofan,LI Jiazong,et al. Numerical simulation on development effect of tar-rich coal through in-situ conversion by convective heating in Huangling Mining Area[J]. Coal Geology & Exploration,2023,51(4):57−67.
[20] 陈美静,漆博文,王长安,等. 富油煤地下原位热解余热利用过程数值模拟[J]. 洁净煤技术,2023,29(1):48−58. CHEN Meijing,QI Bowen,WANG Changan,et al. Numerical simulation of waste heat utilization process after in-situ pyrolysis of tar-rich coal[J]. Clean Coal Technology,2023,29(1):48−58.
[21] JIU Bo,HUANG Wenhui,HE Mingqian,et al. Quantitative analysis of micron-scale and nano-scale pore throat characteristics of tight sandstone using matlab[J]. Applied Sciences,2018,8(8):1272 doi: 10.3390/app8081272
[22] YANG Xiaohan,REN Ting,TAN Lihai. Size distribution measurement of coal fragments using digital imaging processing[J]. Measurement,2020,160:107867. doi: 10.1016/j.measurement.2020.107867
[23] 王登科,魏 强,魏建平,等. 煤的裂隙结构分形特征与分形渗流模型研究[J]. 中国矿业大学学报,2020,49(1):103−−109, 122. WANG Dengke,WEI Qiang,WEI Jianping,et al. Fractal characteristics of fracture structure and fractal seepage model of coal[J]. Journal of China University of Mining & Technology,2020,49(1):103−109, 122.
[24] HU Xin,SUN Qiang,SHI Qingmin,et al. Radon exhalation characteristics after pyrolysis of long flame coal[J]. Science of the Total Environment,2023,904:167228. doi: 10.1016/j.scitotenv.2023.167228
[25] 刘 芸,张 源,郭娅男,等. 不同煤阶煤热解气化特性及机理研究[J]. 科学技术与工程,2021,21(18):7672−7677. LIU Yun,ZHANG Yuan,GUO Yanan,et al. Pyrolysis and gasification of coal with different ranks:characteristics and mechanism[J]. Science Technology and Engineering,2021,21 (18): 7672−7677.
[26] 尹振勇,许 浩,汤达祯,等. 不同煤阶煤热解过程中孔隙结构变化规律研究[J]. 煤炭科学技术,2019,47(9):74−79. YIN Zhenyong,XU Hao,TANG Dazhen,et al. Study on pore structure change during different coal grade pyrolysis[J]. Coal Science and Technology,2019,47(9):74−79.
[27] WEN Hu,LU Junhui,XIAO Yang,et al. Temperature dependence of thermal conductivity,diffusion and specific heat capacity for coal and rocks from coalfield[J]. Thermochimica Acta,2015,619:41−47. doi: 10.1016/j.tca.2015.09.018
[28] DENG Jun,LI Qingwei,XIAO Yang,et al. Predictive models for thermal diffusivity and specific heat capacity of coals in huainan mining area,China[J]. Thermochimica Acta,2017,656:101−111. doi: 10.1016/j.tca.2017.09.005
[29] 邓 军,任帅京,肖 旸,等. 煤低温氧化与热解过程的传热特性对比研究[J]. 煤炭学报,2019,44(S1):171−177. DENG Jun,REN Shuaijing,XIAO Yang,et al. Comparative study on heat transfer characteristics of coal during low temperature oxidation and pyrolysis[J]. Journal of China Coal Society,2019,44(S1):171−177.
[30] ZHUMAGULOV M G. Experimental study of thermophysical properties of Shubarkol coal[J]. Chemistry and Technology Fuels Oils,2013,49:100−108. doi: 10.1007/s10553-013-0419-7
[31] 肖 旸,尹 岚,马 砺,等. 不同预氧化温度下煤样热物性参数的实验研究[J]. 西安科技大学报,2018,38(3):383−388. XIAO Yang,YIN Lan,MA Li,et al. Experimental study on coal thermo-physical parameters under the different peroxidation temperature[J]. Journal of Xi’an University of Science and Technology,2018,38(3):383−388.
[32] BERRYMAN J G. Thermal conductivity of porous media[J]. Applied Physics Letters,2005,86(3):143.
[33] YIN Lan,XIAO Yang,LI Qingwei W,et al. Temperature effect on thermal conductivity of oxidised coal associated with its predictive model during coal pyrolysis[J]. Energy,2023,274:127417. doi: 10.1016/j.energy.2023.127417
[34] SOLOMON P R,FLETCHER T H,PUGMIRE R J. Progress in coal pyrolysis[J]. Fuel,1993,72(5):587−597. doi: 10.1016/0016-2361(93)90570-R
[35] LIU Quanrun,HU Haoquan,ZHOU Qiang,et al. Effect of inorganic matter on reactivity and kinetics of coal pyrolysis[J]. Fuel,2004,83(6):713−718. doi: 10.1016/j.fuel.2003.08.017
[36] SERIO M A,HAMBLEN D G,MARKHAM J R,et al. Kinetics of volatile product evolution in coal pyrolysis:experiment and theory[J]. Energy & Fuels,1987,1(2):138−152.
[37] 温 昶,徐明厚,何兰兰,等. 褐煤热解与燃烧时矿物转化和细灰形成的CCSEM研究[J]. 工程热物理学报,2012,33(4):703−706. WEN Chang,XU Minghou,HE Lanlan,et al. CCSEM Investigation on the mineral transformation and fine ash formation during lignite pyrolysis and combustion[J]. Journal of Engineering Thermophysics,2012,33(4):703−706.
[38] 赵世永,刘美美,徐 婕,等. 神府煤热解前后矿物质变化的研究[J]. 煤炭技术,2016,35(3):278−280. ZHAO Shiyong,LIU Meimei,XU Jie,et al. Study on changes of minerals before and after pyrolysis of Shenfu coal. Coal Technology,2016,35(3):278−280.
[39] COATS A W,REDFERN J. Kinetic parameters from thermogravimetric data[J]. Nature,1964,201:68−69. doi: 10.1038/201068a0
[40] 楚希杰,赵丽红,李 文,等. 神华煤及其直接液化残渣热解动力学试验研究[J]. 煤炭科学技术,2010,38(5):121−124. CHU Xijie,ZHAO Lihong,LI Wen,et al. Research on thermolysis kinetics test of Shenhua coal and residue from direct liquefaction[J]. Coal Science and Technology,2010,38(5):121−124.
[41] 董志浩. 煤炭地下气化覆岩高温损伤与评估研究[D]. 徐州: 中国矿业大学,2020. DONNG Zhihao. Study on thermal damage and evaluation of overburden rock in underground coal gasification[D]. Xuzhou: China University of Mining and Technology,2020.
[42] 畅志兵,初 茉,孙任晖,等. 煤直接液化残渣与褐煤共热解动力学研究[J]. 煤炭科学技术,2015,43(3):138−141,39. CHANG Zhibing,CHU Mo,SUN Renhui,et al. Study on co-pyrolysis kinetics of coal direct iquefaction residue and lignite[J]. Coal Science and Technology,2015,43(3):138−141,39.