CN114113457B - 一种基于双源声信号特征的煤温监测实验装置及方法 - Google Patents

一种基于双源声信号特征的煤温监测实验装置及方法 Download PDF

Info

Publication number
CN114113457B
CN114113457B CN202111428024.3A CN202111428024A CN114113457B CN 114113457 B CN114113457 B CN 114113457B CN 202111428024 A CN202111428024 A CN 202111428024A CN 114113457 B CN114113457 B CN 114113457B
Authority
CN
China
Prior art keywords
coal
coal sample
acoustic wave
receiving
sound
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202111428024.3A
Other languages
English (en)
Other versions
CN114113457A (zh
Inventor
郭军
郑学召
金永飞
蔡国斌
王凯旋
刘荫
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shaanxi Jinchuang Ante Technology Co ltd
Xian University of Science and Technology
Original Assignee
Shaanxi Jinchuang Ante Technology Co ltd
Xian University of Science and Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shaanxi Jinchuang Ante Technology Co ltd, Xian University of Science and Technology filed Critical Shaanxi Jinchuang Ante Technology Co ltd
Priority to CN202111428024.3A priority Critical patent/CN114113457B/zh
Publication of CN114113457A publication Critical patent/CN114113457A/zh
Application granted granted Critical
Publication of CN114113457B publication Critical patent/CN114113457B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N31/00—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods
    • G01N31/12—Investigating or analysing non-biological materials by the use of the chemical methods specified in the subgroup; Apparatus specially adapted for such methods using combustion
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
    • G01K11/22—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of acoustic effects
    • G01K11/24—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using measurement of acoustic effects of the velocity of propagation of sound
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/04—Preparation or injection of sample to be analysed
    • G01N30/06—Preparation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N30/00—Investigating or analysing materials by separation into components using adsorption, absorption or similar phenomena or using ion-exchange, e.g. chromatography or field flow fractionation
    • G01N30/02—Column chromatography
    • G01N30/26—Conditioning of the fluid carrier; Flow patterns
    • G01N30/28—Control of physical parameters of the fluid carrier
    • G01N30/30—Control of physical parameters of the fluid carrier of temperature

Abstract

本发明公开了一种基于双源声信号特征的煤温监测实验装置及方法，该实验装置包括立方体型煤样箱、滤网、电加热丝、进气管和出气管；该方法包括步骤：一、煤样填充并检测立方体型煤样箱的气密性；二、给收发式声波换能器编号并对煤样进行区块划分；三、设置煤样的加热温度阈值并开启电加热丝和气相色谱仪；四、初始化多个收发式声波换能器；五、判断煤样的加热温度是否达到加热温度阈值；六、收发式声波换能器在第一工作模式下的煤温监测；七、收发式声波换能器在第二工作模式下的煤温监测。本发明能够完成外加声波与燃烧音的产生、发声、接收、采集全过程，实时监测煤样温度，实现煤样温度三维可视化，能够快速、准确地重建煤自燃过程中温度场。

Description

一种基于双源声信号特征的煤温监测实验装置及方法
技术领域
本发明属于基于双源声信号特征的煤温监测实验技术领域，具体涉及一种基于双源声信号特征的煤温监测实验装置及方法。
背景技术
煤炭作为我国工业发展的重要物质基础，一直以来占据着能源主体地位。煤层自燃火灾不仅是矿井灾害的主要形式之一，也是诱发矿井瓦斯、煤尘爆炸等次生灾害的主要原因。煤自燃精准探测是高效防控的关键，但是限于采空区空间范围大、煤岩体导热性能差、松散煤体自燃高温区域具有隐蔽性、移动性等特点，以及受限空间火源探测技术瓶颈等因素，火源位置精确探测及防控难度极大。现有火源探测方法在很大程度上解决了火区范围的圈定问题，但仍难以实现采空区煤温的快速量化识别。近年来，声学法测温技术在农业、工业以及大气空间领域的精准测温已广泛应用，为采空区隐蔽火源的精准探测提供了新的思路。然而，在实际运用中，环境中存在大量的噪声信号，且煤体自燃过程中会产生燃烧音，井下受限复杂环境中大范围遗煤自燃产生的燃烧音以及环境噪声在很大程度上制约了声波测温技术发展。其根本原因是：松散煤体自燃历程中燃烧音的演变规律、燃烧音与声发射复合声波的传播机制及感温机理等基础科学问题尚未解决。因此，研究松散煤体声发射信号特征识别与提取算法，分析松散煤体自燃进程燃烧音产生机制及演变规律，探究升温松散煤体内“双源”复合声波信号传播衰减特性，构建基于“双源”复合声波的煤温监测实验装置及方法，以期揭示松散煤体的声波温敏特性及感温机制，研究成果对采空区隐蔽火源探测和煤火灾害精准防控具有重要的价值。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足，提供一种基于双源声信号特征的煤温监测实验装置，其设计新颖合理，能够完成外加声波与燃烧音的产生、发声、接收、采集全过程，实时监测煤样温度，实现煤样温度三维可视化，能够快速、准确地重建煤自燃过程中温度场，便于推广使用。
为解决上述技术问题，本发明采用的技术方案是：一种基于双源声信号特征的煤温监测实验装置，其特征在于：包括立方体型煤样箱，立方体型煤样箱下侧设置有滤网，立方体型煤样箱四个侧壁位于滤网以上的部位设置有电加热丝，滤网、立方体型煤样箱顶板和四侧的电加热丝之间围成用于放置煤样的正方体空腔，立方体型煤样箱下部侧壁设置有进气管，进气管通过立方体型煤样箱底板和滤网之间的间隙伸入至立方体型煤样箱，立方体型煤样箱顶部设置有伸入至立方体型煤样箱内的出气管，出气管远离立方体型煤样箱的一端通过软管与气相色谱仪连通，多个收发式声波换能器穿过立方体型煤样箱外侧壁与立方体型的煤样连通，收发式声波换能器的数量为(8+12a)个，(8+12a)个收发式声波换能器的探测端分别分布在立方体型的煤样的八个顶点和十二个棱边上，其中，a为正整数；
立方体型的所述煤样的四个竖向棱边上的收发式声波换能器均通过声波导管与煤样连通，声波导管由外至内向下倾斜，声波导管穿过立方体型煤样箱外侧壁的部位外侧设置有耐火塞；
所述气相色谱仪的输出端与计算机的输入端连接，电加热丝和多个收发式声波换能器均由所述计算机控制，煤样内预埋有多个温度传感器。
上述的一种基于双源声信号特征的煤温监测实验装置，其特征在于：所述煤样箱包括箱体和罐盖，箱体内煤样顶部与罐盖之间设置有绝热盖。
上述的一种基于双源声信号特征的煤温监测实验装置，其特征在于：所述箱体内侧壁和电加热丝之间由外至内依次设置有保温棉和铜板。
同时，本发明还公开了一种方法步骤简单、设计合理的双源声信号特征的煤温监测的方法，其特征在于：该方法包括以下步骤：
步骤一、煤样填充并检测立方体型煤样箱的气密性：将破碎、筛分好的煤样装入立方体型煤样箱中，将出气管通过橡胶软管连接到气相色谱仪的进样管，并检查煤温监测实验装置的气密性，保证密封效果；
步骤二、给多个收发式声波换能器编号并对正方体煤样进行区块划分；
步骤三、设置煤样的加热温度阈值并开启电加热丝和气相色谱仪；
步骤四、初始化多个收发式声波换能器；
步骤五、判断煤样的加热温度是否达到加热温度阈值：利用多个温度传感器的测量平均值作为煤样的加热温度，当煤样的加热温度未达到加热温度阈值时，多个收发式声波换能器启动第一工作模式，执行步骤六；当煤样的加热温度达到加热温度阈值时，多个收发式声波换能器启动第二工作模式，执行步骤七；
其中，多个收发式声波换能器的第一工作模式是指煤样的加热温度未达到加热温度阈值，煤样未燃烧，立方体型煤样箱内没有声源，利用多个收发式声波换能器制造外加声波进行煤温监测；
多个收发式声波换能器的第二工作模式是指煤样的加热温度达到加热温度阈值，煤样开始燃烧，立方体型煤样箱内产生燃烧音的声源，利用多个收发式声波换能器只接收燃烧音进行煤温监测；
步骤六、多个收发式声波换能器在第一工作模式下的煤温监测，过程如下：
步骤601、依次控制(8+12a)个收发式声波换能器单独工作，其中，控制任一收发式声波换能器单独工作的过程均相同；
控制第q个收发式声波换能器单独工作时，控制第q个收发式声波换能器进行声波发射，并控制剩余收发式声波换能器进行声波接收，其中，当第q个收发式声波换能器进行声波发射后，只有与第q个收发式声波换能器不在一个平面上的收发式声波换能器能接收到声波信号，得到不同传播路径上的声波飞渡时间；
其中，q为正整数且q＝1，2，...，(8+12a)；
步骤602、构建第一病态矩阵
Figure BDA0003377993730000041
其中，M为外加声波下声波传播路径的总数，N为正方体煤样划分的区块总数，Am,n为第一病态矩阵中第m个声波传播路径上第n个区块中的算子且/>
Figure BDA0003377993730000042
lm为外加声波下第m个声波传播路径，(xn,yn,zn)为第n个区块的中心坐标，/>
Figure BDA0003377993730000043
为第一径向基函数，m为正整数且m＝1,2，...，M，n为正整数且n＝1,2，...，N；
步骤603、根据公式tA＝AεA，获取第一待定系数矩阵
Figure BDA0003377993730000044
其中，tA为外加声波下M条声波传播路径的声波飞渡时间矩阵，/>
Figure BDA0003377993730000045
为第一待定系数矩阵中第n个元素；
步骤604、根据公式
Figure BDA0003377993730000046
计算外加声波的声波速度分布函数vA(x,y,z)；
步骤605、根据公式
Figure BDA0003377993730000047
计算外加声波下的温度分布函数TA(x,y,z)，其中，γA为外加声波下气相色谱仪测的气体绝热指数，RA为外加声波下气相色谱仪测的普适气体常量，mA为外加声波下气相色谱仪测的气体的摩尔质量；
步骤606、根据外加声波下的温度分布函数TA(x,y,z)重建煤样整个测量区域内的温度场分布；
步骤七、多个收发式声波换能器在第二工作模式下的煤温监测，过程如下：
步骤701、控制(8+12a)个收发式声波换能器全部为接收状态，根据燃烧音的频率、幅度和接收时间确定其产生位置；
步骤702、构建第二病态矩阵
Figure BDA0003377993730000051
其中，W为燃烧音的声波传播路径的总数，Bw,n为第二病态矩阵中第w个声波传播路径上第n个区块中的算子且
Figure BDA0003377993730000052
lw为燃烧音中第w个声波传播路径，φ(xn,yn,zn)为第二径向基函数，w为正整数且w＝1,2，...，W；
步骤703、根据公式tB＝BεB，获取第二待定系数矩阵
Figure BDA0003377993730000053
其中，tB为燃烧音下W条声波传播路径的声波飞渡时间矩阵，/>
Figure BDA0003377993730000054
为第二待定系数矩阵中第n个元素；
步骤704、根据公式
Figure BDA0003377993730000055
计算燃烧音的声波速度分布函数vB(x,y,z)；
步骤705、根据公式
Figure BDA0003377993730000056
计算燃烧音下的温度分布函数TB(x,y,z)，其中，γB为燃烧音下气相色谱仪测的气体绝热指数，RB为燃烧音下气相色谱仪测的普适气体常量，mB为燃烧音下气相色谱仪测的气体的摩尔质量；
步骤706、根据燃烧音下的温度分布函数TB(x,y,z)重建煤样整个测量区域内的温度场分布。
上述的方法，其特征在于：所述煤样破碎、筛分后粒径选取3mm～5mm；所述煤样为褐煤、长焰煤、瘦煤、贫煤、不粘煤、弱粘煤、无烟煤、气煤、肥煤或焦煤；褐煤的加热温度阈值为270℃～310℃，长焰煤的加热温度阈值为275℃～320℃，瘦煤的加热温度阈值为350℃～380℃，贫煤的加热温度阈值为360℃～385℃，不粘煤的加热温度阈值为280℃～305℃，弱粘煤的加热温度阈值为310℃～350℃，无烟煤的加热温度阈值为370℃～420℃，气煤的加热温度阈值为300℃～350℃，肥煤的加热温度阈值为320℃～360℃，焦煤的加热温度阈值为350℃～370℃。
上述的方法，其特征在于：所述第一径向基函数
Figure BDA0003377993730000057
为Multiquadric径向基函数；所述第二径向基函数φ(xn,yn,zn)为sigmoid径向基函数。
本发明与现有技术相比具有以下优点：
1、本发明采用的实验装置，在收发式换能器发射与接收端加装声波导管，有效降低收发式换能器所在位置的辐射温度以及箱体内煤灰等颗粒物的影响，减少声波在进入煤样内之前的扩散效应，保证了声源信号的初始强度，便于推广使用。
2、本发明采用的实验装置，将装煤箱体内部空间设计为正方体，合理布置收发式声波换能器，得到的外加声波以及燃烧音传播路径清晰有效，使用效果好。
3、本发明根据煤自燃前无燃烧音，自燃后产生燃烧音的特征，利用切换收发式换能器收发状态，分别在煤自燃前利用外加声波测温，自燃后利用燃烧音测温，两者可以起到对照作用，增加测量精度。
4、本发明采用的方法，步骤简单，能够完成外加声波与燃烧音的产生、发声、接收、采集全过程，实时监测煤样温度，实现煤样温度三维可视化，能够快速、准确地重建煤自燃过程中温度场，在三维温度场重建技术中，选用基于径向基拟合和奇异值分解的三维温度场重建方式，综合了径向基函数和奇异值分解的优点，实现待测区域温度场重建，便于推广使用。
综上所述，本发明设计新颖合理，确定煤自燃过程中外加声波和燃烧音的路径和规律，重构出其三维温度场，基于此进行煤自燃规律研究；有利于提高松散煤体自燃进程复合声波煤温感知机制研究，可为采空区隐蔽火源的精准监测提供参考，便于推广使用。
下面通过附图和实施例，对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明采用的实验装置的结构连接示意图。
图2为本发明收发式声波换能器与声波导管的连接示意图。
图3为本发明收发式声波换能器在箱体上的布置位置示意图。
图4为本发明方法的流程框图。
附图标记说明:
1—煤样； 2—箱体； 3—罐盖；
4—绝热盖； 5—保温棉； 6—铜板；
7—耐火塞； 8—电加热丝； 9—进气管；
10—滤网； 11—收发式声波换能器； 12—声波导管；
20—出气管。
具体实施方式
如图1至图3所示，本发明所述的一种基于双源声信号特征的煤温监测实验装置，包括立方体型煤样箱，立方体型煤样箱下侧设置有滤网10，立方体型煤样箱四个侧壁位于滤网10以上的部位设置有电加热丝8，滤网10、立方体型煤样箱顶板和四侧的电加热丝8之间围成用于放置煤样1的正方体空腔，立方体型煤样箱下部侧壁设置有进气管9，进气管9通过立方体型煤样箱底板和滤网10之间的间隙伸入至立方体型煤样箱，立方体型煤样箱顶部设置有伸入至立方体型煤样箱内的出气管20，出气管20远离立方体型煤样箱的一端通过软管与气相色谱仪连通，多个收发式声波换能器11穿过立方体型煤样箱外侧壁与立方体型的煤样1连通，收发式声波换能器11的数量为(8+12a)个，(8+12a)个收发式声波换能器11的探测端分别分布在立方体型的煤样1的八个顶点和十二个棱边上，其中，a为正整数。
本实施例中，所述煤样箱包括箱体2和罐盖3，箱体2内煤样1顶部与罐盖3之间设置有绝热盖4。
本实施例中，所述箱体2内侧壁和电加热丝8之间由外至内依次设置有保温棉5和铜板6。
本实施例中，立方体型的所述煤样1的四个竖向棱边上的收发式声波换能器11均通过声波导管12与煤样1连通，声波导管12由外至内向下倾斜，声波导管12穿过立方体型煤样箱外侧壁的部位外侧设置有耐火塞7。
本实施例中，还包括计算机，所述气相色谱仪的输出端与所述计算机的输入端连接，电加热丝8和多个收发式声波换能器11均由所述计算机控制，煤样1内预埋有多个温度传感器。
需要说明的是，在收发式换能器发射与接收端加装声波导管，有效降低收发式换能器所在位置的辐射温度以及箱体内煤灰等颗粒物的影响，减少声波在进入煤样内之前的扩散效应，保证了声源信号的初始强度；将装煤箱体内部空间设计为正方体，合理布置收发式声波换能器，得到的外加声波以及燃烧音传播路径清晰有效；根据煤自燃前无燃烧音，自燃后产生燃烧音的特征，利用切换收发式换能器收发状态，分别在煤自燃前利用外加声波测温，自燃后利用燃烧音测温，两者可以起到对照作用，增加测量精度；能够完成外加声波与燃烧音的产生、发声、接收、采集全过程，实时监测煤样温度，实现煤样温度三维可视化，能够快速、准确地重建煤自燃过程中温度场，在三维温度场重建技术中，选用基于径向基拟合和奇异值分解的三维温度场重建方式，综合了径向基函数和奇异值分解的优点，实现待测区域温度场重建；确定煤自燃过程中外加声波和燃烧音的路径和规律，重构出其三维温度场，基于此进行煤自燃规律研究；有利于提高松散煤体自燃进程复合声波煤温感知机制研究，可为采空区隐蔽火源的精准监测提供参考。
如图4所示的一种基于双源声信号特征的煤温监测实验方法，包括以下步骤：
步骤一、煤样填充并检测立方体型煤样箱的气密性：将破碎、筛分好的煤样1装入立方体型煤样箱中，将出气管20通过橡胶软管连接到气相色谱仪的进样管，并检查煤温监测实验装置的气密性，保证密封效果；
步骤二、给多个收发式声波换能器编号并对正方体煤样进行区块划分；
步骤三、设置煤样的加热温度阈值并开启电加热丝和气相色谱仪；
步骤四、初始化多个收发式声波换能器；
步骤五、判断煤样的加热温度是否达到加热温度阈值：利用多个温度传感器的测量平均值作为煤样的加热温度，当煤样的加热温度未达到加热温度阈值时，多个收发式声波换能器启动第一工作模式，执行步骤六；当煤样的加热温度达到加热温度阈值时，多个收发式声波换能器启动第二工作模式，执行步骤七；
其中，多个收发式声波换能器的第一工作模式是指煤样的加热温度未达到加热温度阈值，煤样未燃烧，立方体型煤样箱内没有声源，利用多个收发式声波换能器制造外加声波进行煤温监测；
多个收发式声波换能器的第二工作模式是指煤样的加热温度达到加热温度阈值，煤样开始燃烧，立方体型煤样箱内产生燃烧音的声源，利用多个收发式声波换能器只接收燃烧音进行煤温监测；
步骤六、多个收发式声波换能器在第一工作模式下的煤温监测，过程如下：
步骤601、依次控制(8+12a)个收发式声波换能器单独工作，其中，控制任一收发式声波换能器单独工作的过程均相同；
控制第q个收发式声波换能器单独工作时，控制第q个收发式声波换能器进行声波发射，并控制剩余收发式声波换能器进行声波接收，其中，当第q个收发式声波换能器进行声波发射后，只有与第q个收发式声波换能器不在一个平面上的收发式声波换能器11能接收到声波信号，得到不同传播路径上的声波飞渡时间；
其中，q为正整数且q＝1，2，...，(8+12a)；
步骤602、构建第一病态矩阵
Figure BDA0003377993730000091
其中，M为外加声波下声波传播路径的总数，N为正方体煤样划分的区块总数，Am,n为第一病态矩阵中第m个声波传播路径上第n个区块中的算子且/>
Figure BDA0003377993730000101
lm为外加声波下第m个声波传播路径，(xn,yn,zn)为第n个区块的中心坐标，/>
Figure BDA0003377993730000102
为第一径向基函数，m为正整数且m＝1,2，...，M，n为正整数且n＝1,2，...，N；
步骤603、根据公式tA＝AεA，获取第一待定系数矩阵
Figure BDA0003377993730000103
其中，tA为外加声波下M条声波传播路径的声波飞渡时间矩阵，/>
Figure BDA0003377993730000104
为第一待定系数矩阵中第n个元素；
步骤604、根据公式
Figure BDA0003377993730000105
计算外加声波的声波速度分布函数vA(x,y,z)；
步骤605、根据公式
Figure BDA0003377993730000106
计算外加声波下的温度分布函数TA(x,y,z)，其中，γA为外加声波下气相色谱仪测的气体绝热指数，RA为外加声波下气相色谱仪测的普适气体常量，mA为外加声波下气相色谱仪测的气体的摩尔质量；
步骤606、根据外加声波下的温度分布函数TA(x,y,z)重建煤样整个测量区域内的温度场分布；
步骤七、多个收发式声波换能器在第二工作模式下的煤温监测，过程如下：
步骤701、控制(8+12a)个收发式声波换能器全部为接收状态，根据燃烧音的频率、幅度和接收时间确定其产生位置；
步骤702、构建第二病态矩阵
Figure BDA0003377993730000107
其中，W为燃烧音的声波传播路径的总数，Bw,n为第二病态矩阵中第w个声波传播路径上第n个区块中的算子且
Figure BDA0003377993730000108
lw为燃烧音中第w个声波传播路径，φ(xn,yn,zn)为第二径向基函数，w为正整数且w＝1,2，...，W；
步骤703、根据公式tB＝BεB，获取第二待定系数矩阵
Figure BDA0003377993730000111
其中，tB为燃烧音下W条声波传播路径的声波飞渡时间矩阵，/>
Figure BDA0003377993730000112
为第二待定系数矩阵中第n个元素；
步骤704、根据公式
Figure BDA0003377993730000113
计算燃烧音的声波速度分布函数vB(x,y,z)；
步骤705、根据公式
Figure BDA0003377993730000114
计算燃烧音下的温度分布函数TB(x,y,z)，其中，γB为燃烧音下气相色谱仪测的气体绝热指数，RB为燃烧音下气相色谱仪测的普适气体常量，mB为燃烧音下气相色谱仪测的气体的摩尔质量；
步骤706、根据燃烧音下的温度分布函数TB(x,y,z)重建煤样整个测量区域内的温度场分布。
本实施例中，所述煤样1破碎、筛分后粒径选取3mm～5mm；所述煤样1为褐煤、长焰煤、瘦煤、贫煤、不粘煤、弱粘煤、无烟煤、气煤、肥煤或焦煤；褐煤的加热温度阈值为270℃～310℃，长焰煤的加热温度阈值为275℃～320℃，瘦煤的加热温度阈值为350℃～380℃，贫煤的加热温度阈值为360℃～385℃，不粘煤的加热温度阈值为280℃～305℃，弱粘煤的加热温度阈值为310℃～350℃，无烟煤的加热温度阈值为370℃～420℃，气煤的加热温度阈值为300℃～350℃，肥煤的加热温度阈值为320℃～360℃，焦煤的加热温度阈值为350℃～370℃。
本实施例中，所述第一径向基函数
Figure BDA0003377993730000115
为Mul t iquadr ic径向基函数；所述第二径向基函数φ(xn,yn,zn)为s igmoid径向基函数。
以上所述，仅是本发明的较佳实施例，并非对本发明作任何限制，凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化，均仍属于本发明技术方案的保护范围内。

Claims (5)

1.一种利用基于双源声信号特征的煤温监测实验装置进行双源声信号特征的煤温监测的方法，其特征在于：所述基于双源声信号特征的煤温监测实验装置包括立方体型煤样箱，立方体型煤样箱下侧设置有滤网(10)，立方体型煤样箱四个侧壁位于滤网(10)以上的部位设置有电加热丝(8)，滤网(10)、立方体型煤样箱顶板和四侧的电加热丝(8)之间围成用于放置煤样(1)的正方体空腔，立方体型煤样箱下部侧壁设置有进气管(9)，进气管(9)通过立方体型煤样箱底板和滤网(10)之间的间隙伸入至立方体型煤样箱，立方体型煤样箱顶部设置有伸入至立方体型煤样箱内的出气管(20)，出气管(20)远离立方体型煤样箱的一端通过软管与气相色谱仪连通，多个收发式声波换能器(11)穿过立方体型煤样箱外侧壁与立方体型的煤样(1)连通，收发式声波换能器(11)的数量为(8+12a)个，(8+12a)个收发式声波换能器(11)的探测端分别分布在立方体型的煤样(1)的八个顶点和十二个棱边上，其中，a为正整数；
立方体型的所述煤样(1)的四个竖向棱边上的收发式声波换能器(11)均通过声波导管(12)与煤样(1)连通，声波导管(12)由外至内向下倾斜，声波导管(12)穿过立方体型煤样箱外侧壁的部位外侧设置有耐火塞(7)；
所述气相色谱仪的输出端与计算机的输入端连接，电加热丝(8)和多个收发式声波换能器(11)均由所述计算机控制，煤样(1)内预埋有多个温度传感器；
该方法包括以下步骤：
步骤一、煤样填充并检测立方体型煤样箱的气密性：将破碎、筛分好的煤样(1)装入立方体型煤样箱中，将出气管(20)通过橡胶软管连接到气相色谱仪的进样管，并检查煤温监测实验装置的气密性，保证密封效果；
步骤二、给多个收发式声波换能器编号并对正方体煤样进行区块划分；
步骤三、设置煤样的加热温度阈值并开启电加热丝和气相色谱仪；
步骤四、初始化多个收发式声波换能器；
步骤五、判断煤样的加热温度是否达到加热温度阈值：利用多个温度传感器的测量平均值作为煤样的加热温度，当煤样的加热温度未达到加热温度阈值时，多个收发式声波换能器启动第一工作模式，执行步骤六；当煤样的加热温度达到加热温度阈值时，多个收发式声波换能器启动第二工作模式，执行步骤七；
其中，多个收发式声波换能器的第一工作模式是指煤样的加热温度未达到加热温度阈值，煤样未燃烧，立方体型煤样箱内没有声源，利用多个收发式声波换能器制造外加声波进行煤温监测；
多个收发式声波换能器的第二工作模式是指煤样的加热温度达到加热温度阈值，煤样开始燃烧，立方体型煤样箱内产生燃烧音的声源，利用多个收发式声波换能器只接收燃烧音进行煤温监测；
步骤六、多个收发式声波换能器在第一工作模式下的煤温监测，过程如下：
步骤601、依次控制(8+12a)个收发式声波换能器单独工作，其中，控制任一收发式声波换能器单独工作的过程均相同；
控制第q个收发式声波换能器单独工作时，控制第q个收发式声波换能器进行声波发射，并控制剩余收发式声波换能器进行声波接收，其中，当第q个收发式声波换能器进行声波发射后，只有与第q个收发式声波换能器不在一个平面上的收发式声波换能器(11)能接收到声波信号，得到不同传播路径上的声波飞渡时间；
其中，q为正整数且q＝1，2，...，(8+12a)；
步骤602、构建第一病态矩阵
Figure FDA0004254671970000021
其中，M为外加声波下声波传播路径的总数，N为正方体煤样划分的区块总数，Am,n为第一病态矩阵中第m个声波传播路径上第n个区块中的算子且/>
Figure FDA0004254671970000031
lm为外加声波下第m个声波传播路径，(xn,yn,zn)为第n个区块的中心坐标，/>
Figure FDA0004254671970000032
为第一径向基函数，m为正整数且m＝1,2，...，M，n为正整数且n＝1,2，...，N；
步骤603、根据公式tA＝AεA，获取第一待定系数矩阵
Figure FDA0004254671970000033
其中，tA为外加声波下M条声波传播路径的声波飞渡时间矩阵，/>
Figure FDA0004254671970000034
为第一待定系数矩阵中第n个元素；
步骤604、根据公式
Figure FDA0004254671970000035
计算外加声波的声波速度分布函数vA(x,y,z)；
步骤605、根据公式
Figure FDA0004254671970000036
计算外加声波下的温度分布函数TA(x,y,z)，其中，γA为外加声波下气相色谱仪测的气体绝热指数，RA为外加声波下气相色谱仪测的普适气体常量，mA为外加声波下气相色谱仪测的气体的摩尔质量；
步骤606、根据外加声波下的温度分布函数TA(x,y,z)重建煤样整个测量区域内的温度场分布；
步骤七、多个收发式声波换能器在第二工作模式下的煤温监测，过程如下：
步骤701、控制(8+12a)个收发式声波换能器全部为接收状态，根据燃烧音的频率、幅度和接收时间确定其产生位置；
步骤702、构建第二病态矩阵
Figure FDA0004254671970000037
其中，W为燃烧音的声波传播路径的总数，Bw,n为第二病态矩阵中第w个声波传播路径上第n个区块中的算子且
Figure FDA0004254671970000038
lw为燃烧音中第w个声波传播路径，φ(xn,yn,zn)为第二径向基函数，w为正整数且w＝1,2，...，W；
步骤703、根据公式tB＝BεB，获取第二待定系数矩阵
Figure FDA0004254671970000041
其中，tB为燃烧音下W条声波传播路径的声波飞渡时间矩阵，/>
Figure FDA0004254671970000042
为第二待定系数矩阵中第n个元素；
步骤704、根据公式
Figure FDA0004254671970000043
计算燃烧音的声波速度分布函数vB(x,y,z)；
步骤705、根据公式
Figure FDA0004254671970000044
计算燃烧音下的温度分布函数TB(x,y,z)，其中，γB为燃烧音下气相色谱仪测的气体绝热指数，RB为燃烧音下气相色谱仪测的普适气体常量，mB为燃烧音下气相色谱仪测的气体的摩尔质量；
步骤706、根据燃烧音下的温度分布函数TB(x,y,z)重建煤样整个测量区域内的温度场分布。
2.按照权利要求1所述的双源声信号特征的煤温监测的方法，其特征在于：所述煤样箱包括箱体(2)和罐盖(3)，箱体(2)内煤样(1)顶部与罐盖(3)之间设置有绝热盖(4)。
3.按照权利要求2所述的双源声信号特征的煤温监测的方法，其特征在于：所述箱体(2)内侧壁和电加热丝(8)之间由外至内依次设置有保温棉(5)和铜板(6)。
4.按照权利要求1所述的双源声信号特征的煤温监测的方法，其特征在于：所述煤样(1)破碎、筛分后粒径选取3mm～5mm；所述煤样(1)为褐煤、长焰煤、瘦煤、贫煤、不粘煤、弱粘煤、无烟煤、气煤、肥煤或焦煤；褐煤的加热温度阈值为270℃～310℃，长焰煤的加热温度阈值为275℃～320℃，瘦煤的加热温度阈值为350℃～380℃，贫煤的加热温度阈值为360℃～385℃，不粘煤的加热温度阈值为280℃～305℃，弱粘煤的加热温度阈值为310℃～350℃，无烟煤的加热温度阈值为370℃～420℃，气煤的加热温度阈值为300℃～350℃，肥煤的加热温度阈值为320℃～360℃，焦煤的加热温度阈值为350℃～370℃。
5.按照权利要求1所述的双源声信号特征的煤温监测的方法，其特征在于：所述第一径向基函数
Figure FDA0004254671970000051
为Multiquadric径向基函数；所述第二径向基函数φ(xn,yn,zn)为sigmoid径向基函数。
CN202111428024.3A 2021-11-26 2021-11-26 一种基于双源声信号特征的煤温监测实验装置及方法 Active CN114113457B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111428024.3A CN114113457B (zh) 2021-11-26 2021-11-26 一种基于双源声信号特征的煤温监测实验装置及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111428024.3A CN114113457B (zh) 2021-11-26 2021-11-26 一种基于双源声信号特征的煤温监测实验装置及方法

Publications (2)

Publication Number Publication Date
CN114113457A CN114113457A (zh) 2022-03-01
CN114113457B true CN114113457B (zh) 2023-07-14

Family

ID=80370686

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111428024.3A Active CN114113457B (zh) 2021-11-26 2021-11-26 一种基于双源声信号特征的煤温监测实验装置及方法

Country Status (1)

Country Link
CN (1) CN114113457B (zh)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109184800A (zh) * 2018-10-15 2019-01-11 西安科技大学 一种煤自燃预警分级方法及预警方法
CN112647923A (zh) * 2020-12-24 2021-04-13 山东科技大学 大尺度型煤煤炭地下气化相似材料模拟试验装置及方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0774348B2 (ja) * 1993-07-20 1995-08-09 バブコック日立株式会社 石炭ガス化炉スラグ流下監視方法及び装置
US6959761B2 (en) * 2000-04-24 2005-11-01 Shell Oil Company In situ thermal processing of a coal formation with a selected ratio of heat sources to production wells
CN111208555B (zh) * 2020-01-14 2023-03-14 山东科技大学 地下煤火危险声波主被动探测及定位方法
CN111141866B (zh) * 2020-01-14 2022-11-11 山东科技大学 模拟煤层自燃过程及高温异常区域探测的试验装置及方法
CN113295769B (zh) * 2021-05-20 2023-05-30 山东科技大学 煤自燃高温异常区域声电信号传输特性测试的模拟装置及方法
CN113419027A (zh) * 2021-07-28 2021-09-21 西安科技大学 基于声学测温的煤自燃高温点运移路径可视化装置及方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109184800A (zh) * 2018-10-15 2019-01-11 西安科技大学 一种煤自燃预警分级方法及预警方法
CN112647923A (zh) * 2020-12-24 2021-04-13 山东科技大学 大尺度型煤煤炭地下气化相似材料模拟试验装置及方法

Also Published As

Publication number Publication date
CN114113457A (zh) 2022-03-01

Similar Documents

Publication Publication Date Title
CN100456035C (zh) 基于爆破源的大坝无损检测系统
CN108802825B (zh) 一种次声波监测煤岩动力灾害定位方法及定位系统
Blair et al. Attenuation of explosion‐generated pulse in rock masses
Lu et al. Case study of blast-induced shock wave propagation in coal and rock
CN103235335A (zh) 一种强声脉冲测井系统
Zhang et al. Experimental study on intensity and energy evolution of deep coal and gas outburst
CN105588783A (zh) 一种超声波监测管道输送膏体充填料浆浓度的方法
CN114113457B (zh) 一种基于双源声信号特征的煤温监测实验装置及方法
CN114412573A (zh) 一种矿井井壁破裂监测预警系统及方法
Chen et al. Correlation mechanism between the law of ultrasonic propagation in coal samples and the migration of water
Ren et al. Sound absorption characteristics of loose bituminous coal porous media with different metamorphic degrees
Vivek et al. Laboratory scale investigation of stress wave propagation and vibrational characteristics in sand when subjected to air-blast loading
Manchao et al. Laboratory study of unloading rate effects on rockburst
Dong et al. Influence of temperature on acoustic emission source location accuracy in underground structure
CN109187761A (zh) 一种气固热耦合的煤岩体滞弹性各向异性探测装置及方法
Matoza et al. One hundred years of advances in volcano seismology and acoustics
CN110018062B (zh) 一种直剪试验中岩石结构面剪切破坏位置定位方法
US20040177692A1 (en) System and method for inspecting an industrial furnace or the like
Zhang et al. Anisotropic characteristics of ultrasonic transmission velocities and stress inversion during uniaxial compression process
CN216386914U (zh) 一种煤自燃全过程声信号与温度映射关系构建实验装置
CN108167090B (zh) 一种固体火箭发动机试验三维燃面退移测试系统与方法
CN102298158A (zh) 断层探查装置以及断层探查方法
CN113295769B (zh) 煤自燃高温异常区域声电信号传输特性测试的模拟装置及方法
Jia et al. The analysis of energy characteristics of infrasonic signal in coal sample with different strength under loading process
CN113419027A (zh) 基于声学测温的煤自燃高温点运移路径可视化装置及方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant