CN112816653A - 预混可燃气体爆炸火焰传播速度研究装置和方法 - Google Patents

预混可燃气体爆炸火焰传播速度研究装置和方法 Download PDF

Info

Publication number
CN112816653A
CN112816653A CN202110226048.4A CN202110226048A CN112816653A CN 112816653 A CN112816653 A CN 112816653A CN 202110226048 A CN202110226048 A CN 202110226048A CN 112816653 A CN112816653 A CN 112816653A
Authority
CN
China
Prior art keywords
explosion
combustible gas
pipeline
propagation speed
premixed
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.)
Pending
Application number
CN202110226048.4A
Other languages
English (en)
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.)
Beijing Institute of Petrochemical Technology
Original Assignee
Beijing Institute of Petrochemical 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 Beijing Institute of Petrochemical Technology filed Critical Beijing Institute of Petrochemical Technology
Priority to CN202110226048.4A priority Critical patent/CN112816653A/zh
Publication of CN112816653A publication Critical patent/CN112816653A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/22—Fuels; Explosives
    • G01N33/225—Gaseous fuels, e.g. natural gas

Landscapes

  • Chemical & Material Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Medicinal Chemistry (AREA)
  • Oil, Petroleum & Natural Gas (AREA)
  • General Chemical & Material Sciences (AREA)
  • Food Science & Technology (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Physics & Mathematics (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明公开了一种预混可燃气体爆炸火焰传播速度研究装置和方法，包括预混可燃气体爆炸管道、爆速仪、可燃气体钢瓶和点火系统；预混可燃气体爆炸管道为密闭不透光管道，其一端装有压力表；点火系统包括点火电极，点火电极装于预混可燃气体爆炸管道的一端；爆速仪包括多个光纤传感器，每个光纤传感器分别通过固定调节器装于预混可燃气体爆炸管道的一侧；预混可燃气体爆炸管道的侧壁连接有出气阀门、进气阀门、循环泵和可燃气体钢瓶。可以对预混可燃气体爆炸开展火焰传播速度测量研究，拟合出预混可燃气体爆炸火焰传播速度与管道长度的关系，为解释预混可燃气体爆炸特性、寻求抑制气体爆炸新技术有重要意义。

Description

预混可燃气体爆炸火焰传播速度研究装置和方法
技术领域
本发明专利涉及一种气体爆炸特性研究技术和装置，尤其涉及一种预混可燃气体爆炸火焰传播速度研究装置和方法，可以在难以实际测定情况下，对预混可燃气体爆炸火焰传播速度进行数值模拟计算。
背景技术
随着现代工业的发展，生产工艺越来越复杂，生产集约化水平也越来越高，工业爆炸的危险也随之增高，气体爆炸发生的频率很高，造成的损失也最大。可燃气体爆炸速度、压力和温度是可燃气体爆炸的三个重要特征，加强可燃气体爆炸基础数据的研究，对研究可燃气体爆炸反应和防治技术有重要意义。
可燃气体爆炸燃烧分为爆炸和爆轰阶段，火焰传播速度从数米/每秒到数千米/每秒不等。火焰传播速度取决于火焰阵面前方扰动气体的运动速度，火焰前方扰动气流的速度超过或接近音速时，会产生巨大的冲击波，可燃气体爆炸燃烧从爆炸转向爆轰阶段，爆炸到爆轰是一个量变到质变的过程。可燃气体爆炸火焰传播速度是火焰相对于前方已扰动气体的运动速度，它与反应物质有关，是反应物质的特征量，不同反应物质爆炸燃烧速度也不同。
国内外众多学者用不同仪器测量不同物质的爆炸燃烧火焰速度，如爆速仪、高速摄像仪等等。
爆速仪多用测量炸药爆轰的速度，杨学山等用一种爆速仪测得雷管的爆轰速度为3000m/s左右。郭亚飞等人用高精度多段智能爆速仪测量了雷管爆轰的速度为7000m/s左右；
学者多用高速摄像仪测量预混可燃气体爆炸火焰速度，高建村等人在1m的圆形管道用高速摄像得出了苯的火焰传播特性，苯蒸气的浓度在3.2％时，火焰的平均速度为4.63m/s。杜增晖等用高速摄像测量了260mm*260mm*315mm的方形管道的甲烷-空气燃爆速度为1.103m/s。魏旭星等用高速纹影摄像法测量了直径130mm、长310mm的圆形管道甲烷的爆炸火焰速度为1.162m/s。卢矍然等用高速纹影摄像法测量了直径140mm长度180mm的圆形管道甲烷的平均速度为1.25m/s。高全研究甲烷-空气预混气体前期最大自加速为1.36m/s。实验研究大多是定容爆炸，影响可燃气体爆炸燃烧速度有诸多因素，如实验管道的长度、宽度、直径比和管道内的构造等等，预混可燃气体爆炸燃烧实验研究达到燃烧的爆炸阶段，火焰传播速度远远达不到爆轰的速度。
火焰传播速度是可燃气体爆炸特征的重要参数，火焰传播速度快慢直接反应了燃爆反应的特征。目前还没有一种预混可燃气体爆炸传播速度的测量装置和方法，在难以实际测定情况下，对不同长度管道内预混可燃气体爆炸火焰传播速度进行数值模拟计算和预测。
发明内容
本发明的目的是提供一种预混可燃气体爆炸火焰传播速度研究装置和方法。
本发明的目的是通过以下技术方案实现的：
本发明的预混可燃气体爆炸火焰传播速度研究装置，包括预混可燃气体爆炸管道、爆速仪、可燃气体钢瓶和点火系统；
所述预混可燃气体爆炸管道为密闭不透光管道，其一端装有压力表；
所述点火系统包括点火电极，所述点火电极装于所述预混可燃气体爆炸管道的一端；
所述爆速仪包括多个光纤传感器，每个光纤传感器分别通过固定调节器装于所述预混可燃气体爆炸管道的一侧；
所述预混可燃气体爆炸管道的侧壁的一端连接有出气阀门，另一端连接有进气阀门，出气阀门与进气阀门之间通过管道连接，并在连接管道上设有循环泵，所述可燃气体钢瓶与所述循环泵与进气阀门之间的管道连接。
上述的预混可燃气体爆炸火焰传播速度研究装置进行预混可燃气体爆炸火焰传播速度研究的方法，包括步骤：
A、安装光纤传感器，将光纤传感器安装到固定调节器上，使其成为一体，在实验管道上打孔，将光纤传感器和固定调节器安装到实验管道上，光纤传感器与实验管道齐平或超出实验管道少许部分，伸进实验管道的部分光纤传感器做耐高温处理，不局限于材料和方法，固定调节器起到固定光纤传感器、调节光纤传感器高度和保证实验管道的气密性良好的作用，安装三个光纤传感器和固定调节器到实验管道上，每个光纤传感器的距离是300mm；
B、关闭进气阀，打开出气阀，开启循环泵，抽出管道内空气，用压力表对实验装置进行测试，确保管道负压下具有良好的气密性；
C、确保实验装置气密性良好，用循环泵使管道内形成负压，关闭出气阀，打开进气阀，向管道内准确注入计算量体积的可燃气体，并向管道内补入空气使之形成常压；
D、打开出气阀，进气阀和出气阀均连接循环泵，循环管道内气体，依次关闭出气阀、进气阀和循环泵，静置；
E、将光纤传感器连接好爆速仪，设置爆速仪测量参数，将爆速仪设置成测量状态；
F、接好点火电极，点火系统设置点火能，待静置完成，点火；
G、待爆炸结束后，打开出气阀，尾气进行无害化处理，收集爆炸后可燃气体爆炸燃烧火焰传播速度数据；
H、把数据导入数据处理软件不限于Matlab，拟合爆炸火焰传播速度与管道长度的关系式。
由上述本发明提供的技术方案可以看出，本发明提供的预混可燃气体爆炸火焰传播速度研究装置和方法，可以对预混可燃气体爆炸开展火焰传播速度测量研究，拟合出(不局限于方法)预混可燃气体爆炸火焰传播速度与管道长度的关系，为解释预混可燃气体爆炸特性、寻求抑制气体爆炸新技术有重要意义。
附图说明
图1为本发明实施例提供的预混可燃气体爆炸火焰传播速度研究装置的原理示意图；
图2为本发明实施例Matlab Curve Fitting Tool操作界面；
图3为本发明实施例9.5％甲烷-空气预混可燃气体爆炸火焰传播速度与管道长度的关系曲线；
图4为本发明实施例5％丙烷-空气预混可燃气体爆炸火焰传播速度与管道长度的关系曲线。
图中：
1-实验管道(预混可燃气体爆炸管道)、2-爆速仪、3-光纤传感器、4-固定调节器、5-压力表、6-点火电极、7-点火系统、8-出气阀门、9-循环泵、10-进气阀门、11-可燃气体钢瓶。
具体实施方式
下面将对本发明实施作进一步详细描述。本发明实施例中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
本发明的预混可燃气体爆炸火焰传播速度研究装置，其较佳的具体实施方式是：
包括预混可燃气体爆炸管道、爆速仪、可燃气体钢瓶和点火系统；
所述预混可燃气体爆炸管道为密闭不透光管道，其一端装有压力表；
所述点火系统包括点火电极，所述点火电极装于所述预混可燃气体爆炸管道的一端；
所述爆速仪包括多个光纤传感器，每个光纤传感器分别通过固定调节器装于所述预混可燃气体爆炸管道的一侧；
所述预混可燃气体爆炸管道的侧壁的一端连接有出气阀门，另一端连接有进气阀门，出气阀门与进气阀门之间通过管道连接，并在连接管道上设有循环泵，所述可燃气体钢瓶与所述循环泵与进气阀门之间的管道连接。
所述点火系统为高能点火系统。
所述爆速仪为多段式高精度爆速仪，优选包括8个光纤传感器，8个光纤传感器沿所述预混可燃气体爆炸管道的轴向均布。
所述光纤传感器为全反射感光元件，伸进所述预混可燃气体爆炸管道部分的光纤传感器为耐高温部件。
所述固定调节器为固定和调节所述光纤传感器的高度的部件，其与所述预混可燃气体爆炸管道之间设有密封元件。
上述的预混可燃气体爆炸火焰传播速度研究装置进行预混可燃气体爆炸火焰传播速度研究的方法，包括步骤：
A、安装光纤传感器，将光纤传感器安装到固定调节器上，使其成为一体，在实验管道上打孔，将光纤传感器和固定调节器安装到实验管道上，光纤传感器与实验管道齐平或超出实验管道少许部分，伸进实验管道的部分光纤传感器做耐高温处理，不局限于材料和方法，固定调节器起到固定光纤传感器、调节光纤传感器高度和保证实验管道的气密性良好的作用，安装三个光纤传感器和固定调节器到实验管道上，每个光纤传感器的距离是300mm；
B、关闭进气阀，打开出气阀，开启循环泵，抽出管道内空气，用压力表对实验装置进行测试，确保管道负压下具有良好的气密性；
C、确保实验装置气密性良好，用循环泵使管道内形成负压，关闭出气阀，打开进气阀，向管道内准确注入计算量体积的可燃气体，并向管道内补入空气使之形成常压；
D、打开出气阀，进气阀和出气阀均连接循环泵，循环管道内气体，依次关闭出气阀、进气阀和循环泵，静置；
E、将光纤传感器连接好爆速仪，设置爆速仪测量参数，将爆速仪设置成测量状态；
F、接好点火电极，点火系统设置点火能，待静置完成，点火；
G、待爆炸结束后，打开出气阀，尾气进行无害化处理，收集爆炸后可燃气体爆炸燃烧火焰传播速度数据；
H、把数据导入数据处理软件不限于Matlab，拟合爆炸火焰传播速度与管道长度的关系式。
进行以下任一种实验，及其数据分析：
5-15％V/V甲烷-空气预混气体爆炸火焰传播速度研究；
2.1-9.5％V/V丙烷-空气预混气体爆炸火焰传播速度研究；
任一种在爆炸极限范围内的可燃气体爆炸火焰传播速度研究；
不同强度磁场下在爆炸极限范围内的可燃气体爆炸火焰传播速度研究；
填充不同阻隔防爆材料的可燃气体爆炸火焰传播速度研究；
不同强度磁场下填充不同阻隔防爆材料的可燃气体爆炸火焰传播速度研究。
利用实验加数值模拟，拟合不同条件不同浓度不同种类预混可燃气体爆炸火焰传播速度与管道长度的关系式，根据关系式预测在难以实际测定情况下，对预混可燃气体爆炸火焰传播速度进行数值模拟计算。
预混可燃气体爆炸火焰传播速度与管道长度的关系式如下：
v＝a*lb+c；
其中v代表爆炸火焰传播速度，单位是m/s，I代表管道长度，单位是m，a、b、c是系数；
系数a、b、c通过下式进行所述关系式的拟合：
Figure BDA0002956230450000051
R2代表拟合优度，在[0，1]之间，值越大代表关系式拟合的越好，
Figure BDA0002956230450000052
是回归方程的平均值，
Figure BDA0002956230450000053
代表回归方程的总方差，
Figure BDA0002956230450000054
是所有y值的平均值，
Figure BDA0002956230450000055
是y的总方差，用1减去y对回归方程的方差(未解释离差)与y的总方差的比值，y减去
Figure BDA0002956230450000056
也就是残差，是拟合方程中不能解释的部分，用1减去不能解释的部分，那么剩下的就是解释的部分，也就是说自变量解释了因变量变动的百分比的多少，那么R2的值肯定是越大越好，意味着该模型把y的变动解释得好；
不同浓度不同气体拟合出来的爆炸火焰传播速度和管道长度的关系式的系数a、b、c是不同的，取决于可燃气体所含的能量大小，能量大的可燃气体拟合出来的关系式的系数a就越大，爆炸火焰传播速度随着管道长度越长，产生的击波更大，爆炸火焰传播速度增加越快，爆炸火焰传播速度v与管道长度l是幂函数方程关系，管道长度l是自变量，爆炸火焰传播速度v是幂也就是因变量，b是自变量管道长度l的指数，c是修正方程后为了减小误差的补充数，让调整R2更大，更接近于1，拟合的方程越贴近实验数值。
对于9.5％甲烷-空气预混可燃气体在圆形管道内爆炸火焰速度传播的研究，v＝8.106*l4.528+2.235；
对于5％丙烷-空气预混可燃气体在圆形管道内爆炸火焰速度传播的研究，v＝17.45*l4.211+1.747。
本发明为预混可燃气体爆炸速度研究装置和方法，可以在可燃气体爆炸实验管道中对不同浓度不同预混可燃气体和其他不同条件下，测量不同浓度不同预混可燃气体在不同长度的管道爆炸火焰传播速度，利用数据处理软件不限于Mathlab工具软件拟合爆炸火焰传播速度和管道长度关系式，根据这个关系式可以得到不同管道长度对应的预混可燃气体爆炸火焰传播速度，可以对火焰特征进行更好的研究。
具体实施例：
一、实验准备与步骤：
实验系统主要有预混可燃气体爆炸实验管道、爆速仪、固定调节器和点火系统组成，本实施例以下面为例，但不仅限于此装置和设备。
1、预混可燃气体爆炸实验管道为不限于自行设计加工的密闭不透光圆形管道。
2、爆速仪为高精度多段式爆速仪，可以连接8个光纤传感器但不仅限于此数量，测量多段爆炸火焰燃烧的速度。测量爆炸火焰速度的方法也不仅限于爆速仪。
3、固定调节器为固定光纤传感器和调节光纤传感器高度、保证试验管道气密性良好的元件，固定调节器做耐高温处理，不局限于材料和方法。
4、点火系统由高能点火系统和点火电极组成，点火能量1-1000mJ之间，点火能可控。
实验系统如图1所示。
1.2实验仪器和材料
根据化学反应当量摩尔比计算结果，配制甲烷-空气预混气体浓度为9.5％(V/V)、丙烷-空气预混气体浓度为5％(V/V)，乙烯-空气预混气体浓度为。
CH4+2O2＝2H2O+CO2△H＝-413kJ/mol
C3H8+5O2＝4H2O+3CO2△H＝-2217.8kJ/mol
二、利用Mathlab进行数值模拟
实验数据处理用到Matlab中Curve Fitting Tool，界面如图2，
在命令行窗口输入cftool调出工具界面，选择数据源，管道长度作为横轴变量，爆炸火焰速度作为竖轴变量，在命令窗口输入实验测得的数据，在列表中选择拟合方式，经过对比本实验选择Curve Fitting Tool中的Power函数类型最为接近，运行的结果如图3和图4，爆速仪是根据测量两个光纤传感器距离的时间差得出管道距离的平均爆炸火焰传播速度，采用10组数据的平均值，每组数据为5段距离对应的爆炸火焰传播速度的平均值，拟合出爆炸火焰传播速度与管道长度的关系式，用第6段数据做对比，每次实验数据不一样，拟合得出的关系式用的代码也不一样，会在具体实施例中表示。
三、实施例(以以下实验为例，不限于以下实验)
3.1 9.5％甲烷-空气爆炸火焰传播速度研究实施例
3.1.1 1m管道9.5％甲烷-空气爆炸火焰传播速度
在1m的管道内测量9.5％甲烷-空气爆炸火焰传播速度，测得3组数据，如表1所示：
表11m管道9.5％甲烷-空气爆炸火焰传播速度
Figure BDA0002956230450000071
由上表可知，在1m管道内9.5％甲烷-空气爆炸火焰传播第一段的平均速度为1.12075m/s，在0.75m处处爆炸火焰传播速度为5.23619m/s。
3.1.2 2m管道9.5％甲烷-空气爆炸火焰传播速度
在2m的管道内测量9.5％甲烷-空气爆炸火焰传播速度，测得10组数据做平均值，如表2所示：
表22m管道9.5％甲烷-空气爆炸火焰传播速度
Figure BDA0002956230450000072
由上表可知，在2m管道内9.5％甲烷-空气第一段距离的爆炸火焰传播速度为1.08576m/s，第2段距离爆炸火焰传播速度为5.32914m/s，第3段距离爆炸火焰传播速度为15.94572m/s，第4段距离爆炸火焰传播速度为44.81628m/s，第5段距离爆炸火焰传播速度为104.680985m/s。
采用5组数据运用Matlab拟合出了9.5％甲烷-空气预混可燃气体爆炸火焰传播速度与管道长度的关系式，代码如下：
function[fitresult,gof]＝createFit(PipeLength,FlamePropagatingVelocity)
％％Fit:'untitled fit 1'.
[xData,yData]＝prepareCurveData(PipeLength,FlamePropagatingVelocity)；
％Set up fittype and options.
ft＝fittype('power2')；
opts＝fitoptions('Method','NonlinearLeastSquares')；
opts.Display＝'Off'；
opts.StartPoint＝[14.0533494806809 3.165992480937473.57257035472464]；
％Fit model to data.
[fitresult,gof]＝fit(xData,yData,ft,opts)；
％Plot fit with data.
figure('Name','untitled fit 1')；
h＝plot(fitresult,xData,yData)；
legend(h,'FlamePropagatingVelocity vs.PipeLength','untitled fit 1','Location','NorthEast')；
％Label axes
xlabelPipeLength
ylabelFlamePropagatingVelocity
grid off
由Matlab可得，关系式如下：
v＝8.106*l4.528+2.235
其中v代表爆炸火焰传播速度，单位是m/s，l代表管道长度，单位是m。
Figure BDA0002956230450000081
由Matlab可得，R-square＝0.9983
Adjust R-square＝0.9986
拟合优度为0.9983，接近于1，证明拟合度很好，调整R2为0.9986，让关系式拟合度进一步上升。
3.1.3 3m管道9.5％甲烷-空气爆炸火焰传播速度
在3m的管道内测量9.5％甲烷-空气爆炸火焰传播速度，用实际测得最后一段距离的爆炸火焰传播速度数据与拟合的火焰传播速度做对比，来验证拟合的关系式是否正确，如表3：
表3实际测量的爆炸火焰传播速度与拟合的爆炸火焰传播速度对比
Figure BDA0002956230450000091
用第6段距离的爆炸火焰传播速度验证关系式的拟合度，拟合的爆炸火焰传播速度为189.24741m/s，实际测到的爆炸火焰传播速度为187.57924m/s，由于甲烷-空气预混可燃气体爆炸每一次的火焰传播速度不可能完全相同，根据公式得出的爆炸火焰传播速度与爆速仪实际测出的爆炸火焰传播速度误差很小，误差完全在可接受的范围内，证明拟合的预混可燃气体爆炸火焰传播速度与管道长度的关系式是正确的。
3.25％丙烷-空气爆炸火焰传播速度研究实施例
3.2.1 1m管道5％丙烷-空气爆炸火焰传播速度
测量1m管道5％丙烷-空气爆炸火焰传播速度，测得3组数据，如表4所示：
表4 1m管道5％丙烷-空气爆炸火焰传播速度
Figure BDA0002956230450000092
由上表所知，1m管道5％丙烷-空气预混可燃气体在0.45m处平均爆炸火焰传播速度为1.75797m/s，在0.75m处爆炸火焰传播速度为7.95462m/s。
3.2.2 2m管道5％丙烷-空气爆炸火焰传播速度
测量2m管道5％丙烷-空气爆炸火焰传播速度，测得10组数据做平均值，如表5所示：
表5 2m管道5％丙烷-空气爆炸火焰传播速度
Figure BDA0002956230450000101
由上表可知，在2m管道内5％丙烷-空气第一段爆炸火焰传播速度为1.76819m/s，第2段爆炸火焰传播速度为7.83572m/s，第3段爆炸火焰传播速度为27.43914m/s，第4段爆炸火焰传播速度为85.25016m/s，第5段平均爆炸火焰传播速度为185.92309m/s。
采用5段距离对应的爆炸火焰传播速度数据运用Matlab拟合5％丙烷-空气预混可燃气体爆炸火焰传播速度与管道长度的关系式，代码如下：
function[fitresult,gof]＝createFit(PipeLength,FlamePropagatingVelocity)
％％Fit:'untitled fit 1'.
[xData,yData]＝prepareCurveData(PipeLength,FlamePropagatingVelocity)；
％Set up fittype and options.
ft＝fittype('power2')；
opts＝fitoptions('Method','NonlinearLeastSquares')；
opts.Display＝'Off'；
opts.StartPoint＝[23.9494344566351 3.180559460912748.84502066822748]；
％Fit model to data.
[fitresult,gof]＝fit(xData,yData,ft,opts)；
％Plot fit with data.
figure('Name','untitled fit 1')；
h＝plot(fitresult,xData,yData)；
legend(h,'FlamePropagatingVelocity vs.PipeLength','untitled fit 1','Location','NorthEast')；
％Label axes
xlabelPipeLength
ylabelFlamePropagatingVelocity
grid off
由Matlab可得，关系式如下：
v＝17.45*l4.211+1.747
其中v代表爆炸火焰传播速度，单位是m/s，l代表管道长度，单位是m。
Figure BDA0002956230450000102
由Matlab可得，R-square＝0.9999
Adjust R-square＝0.9999
拟合优度为0.9999，接近于1，证明拟合度很好，调整R2为0.9999。
3.2.33m管道5％丙烷-空气爆炸火焰传播速度
在3m的管道内测量5％丙烷-空气爆炸火焰传播速度，用实际测得最后一段距离的爆炸火焰传播速度数据与拟合的火焰传播速度做对比，来验证拟合的关系式是否正确，如表6：
表65％丙烷-空气实际测量的爆炸火焰传播速度与拟合的爆炸火焰传播速度对比
Figure BDA0002956230450000111
用第6段距离的爆炸火焰传播速度验证关系式的拟合度，拟合的爆炸火焰传播速度为324.91827m/s，实际测到的爆炸火焰传播速度为329.32738m/s，由于丙烷-空气预混可燃气体爆炸每一次的火焰传播速度不可能完全相同，根据公式得出的爆炸火焰传播速度与爆速仪实际测出的爆炸火焰传播速度误差很小，误差完全在可接受的范围内，证明拟合的5％丙烷-空气预混可燃气体爆炸火焰传播速度与管道长度的关系式同样是正确的，也证明此方法同样适用不同浓度的丙烷-空气预混可燃气体爆炸火焰速度研究。
由于实验条件的限制，实验管道不可能无限长，本发明提出的预混可燃气体爆炸火焰传播速度研究技术利用实验加数值模拟可以得到不同浓度不同预混可燃气体的爆炸火焰传播速度与管道长度的关系，也可以预测不同浓度不同预混可燃气体在更长管道的爆炸火焰速度。
以上所述，仅为本发明较佳的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明披露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求书的保护范围为准。

Claims (10)

1.一种预混可燃气体爆炸火焰传播速度研究装置，其特征在于，包括预混可燃气体爆炸管道、爆速仪、可燃气体钢瓶和点火系统；
所述预混可燃气体爆炸管道为密闭不透光管道，其一端装有压力表；
所述点火系统包括点火电极，所述点火电极装于所述预混可燃气体爆炸管道的一端；
所述爆速仪包括多个光纤传感器，每个光纤传感器分别通过固定调节器装于所述预混可燃气体爆炸管道的一侧；
所述预混可燃气体爆炸管道的侧壁的一端连接有出气阀门，另一端连接有进气阀门，出气阀门与进气阀门之间通过管道连接，并在连接管道上设有循环泵，所述可燃气体钢瓶与所述循环泵与进气阀门之间的管道连接。
2.根据权利要求1所述的预混可燃气体爆炸火焰传播速度研究装置，其特征在于，所述点火系统为高能点火系统。
3.根据权利要求1所述的预混可燃气体爆炸火焰传播速度研究装置，其特征在于，所述爆速仪为多段式高精度爆速仪，包括8个光纤传感器，8个光纤传感器沿所述预混可燃气体爆炸管道的轴向均布。
4.根据权利要求1所述的预混可燃气体爆炸火焰传播速度研究装置，其特征在于，所述光纤传感器为全反射感光元件，伸进所述预混可燃气体爆炸管道部分的光纤传感器为耐高温部件。
5.根据权利要求1所述的预混可燃气体爆炸火焰传播速度研究装置，其特征在于，所述固定调节器为固定和调节所述光纤传感器的高度的部件，其与所述预混可燃气体爆炸管道之间设有密封元件。
6.一种权利要去1至5任一项所述的预混可燃气体爆炸火焰传播速度研究装置进行预混可燃气体爆炸火焰传播速度研究的方法，其特征在于，包括步骤：
A、安装光纤传感器，将光纤传感器安装到固定调节器上，使其成为一体，在实验管道上打孔，将光纤传感器和固定调节器安装到实验管道上，光纤传感器与实验管道齐平或超出实验管道少许部分，伸进实验管道的部分光纤传感器做耐高温处理，不局限于材料和方法，固定调节器起到固定光纤传感器、调节光纤传感器高度和保证实验管道的气密性良好的作用，安装三个光纤传感器和固定调节器到实验管道上，每个光纤传感器的距离是300mm；
B、关闭进气阀，打开出气阀，开启循环泵，抽出管道内空气，用压力表对实验装置进行测试，确保管道负压下具有良好的气密性；
C、确保实验装置气密性良好，用循环泵使管道内形成负压，关闭出气阀，打开进气阀，向管道内准确注入计算量体积的可燃气体，并向管道内补入空气使之形成常压；
D、打开出气阀，进气阀和出气阀均连接循环泵，循环管道内气体，依次关闭出气阀、进气阀和循环泵，静置；
E、将光纤传感器连接好爆速仪，设置爆速仪测量参数，将爆速仪设置成测量状态；
F、接好点火电极，点火系统设置点火能，待静置完成，点火；
G、待爆炸结束后，打开出气阀，尾气进行无害化处理，收集爆炸后可燃气体爆炸燃烧火焰传播速度数据；
H、把数据导入数据处理软件不限于Matlab，拟合爆炸火焰传播速度与管道长度的关系式。
7.根据权利要求6所述的预混可燃气体爆炸火焰传播速度研究的方法，其特征在于，进行以下任一种实验，及其数据分析：
5-15％V/V甲烷-空气预混气体爆炸火焰传播速度研究；
2.1-9.5％V/V丙烷-空气预混气体爆炸火焰传播速度研究；
任一种在爆炸极限范围内的可燃气体爆炸火焰传播速度研究；
不同强度磁场下在爆炸极限范围内的可燃气体爆炸火焰传播速度研究；
填充不同阻隔防爆材料的可燃气体爆炸火焰传播速度研究；
不同强度磁场下填充不同阻隔防爆材料的可燃气体爆炸火焰传播速度研究。
8.根据权利要求7所述的预混可燃气体爆炸火焰传播速度研究的方法，其特征在于，利用实验加数值模拟，拟合不同条件不同浓度不同种类预混可燃气体爆炸火焰传播速度与管道长度的关系式，根据关系式预测在难以实际测定情况下，对预混可燃气体爆炸火焰传播速度进行数值模拟计算。
9.根据权利要求8所述的预混可燃气体爆炸火焰传播速度研究的方法，其特征在于，预混可燃气体爆炸火焰传播速度与管道长度的关系式如下：
v＝a*lb+c；
其中v代表爆炸火焰传播速度，单位是m/s，l代表管道长度，单位是m，a、b、c是系数；
系数a、b、c通过下式进行所述关系式的拟合：
Figure FDA0002956230440000031
R2代表拟合优度，在[0,1]之间，值越大代表关系式拟合的越好，
Figure FDA0002956230440000032
是回归方程的平均值，
Figure FDA0002956230440000033
代表回归方程的总方差，
Figure FDA0002956230440000034
是所有y值的平均值，
Figure FDA0002956230440000035
是y的总方差，用1减去y对回归方程的方差(未解释离差)与y的总方差的比值，y减去
Figure FDA0002956230440000036
也就是残差，是拟合方程中不能解释的部分，用1减去不能解释的部分，那么剩下的就是解释的部分，也就是说自变量解释了因变量变动的百分比的多少，那么R2的值肯定是越大越好，意味着该模型把y的变动解释得好；
不同浓度不同气体拟合出来的爆炸火焰传播速度和管道长度的关系式的系数a、b、c是不同的，取决于可燃气体所含的能量大小，能量大的可燃气体拟合出来的关系式的系数a就越大，爆炸火焰传播速度随着管道长度越长，产生的击波更大，爆炸火焰传播速度增加越快，爆炸火焰传播速度v与管道长度l是幂函数方程关系，管道长度l是自变量，爆炸火焰传播速度v是幂也就是因变量，b是自变量管道长度l的指数，c是修正方程后为了减小误差的补充数，让调整R2更大，更接近于1，拟合的方程越贴近实验数值。
10.根据权利要求8所述的预混可燃气体爆炸火焰传播速度研究的方法，其特征在于：
对于9.5％甲烷-空气预混可燃气体在圆形管道内爆炸火焰速度传播的研究，v＝8.106*l4.528+2.235；
对于5％丙烷-空气预混可燃气体在圆形管道内爆炸火焰速度传播的研究，v＝17.45*l4.211+1.747。
CN202110226048.4A 2021-03-01 2021-03-01 预混可燃气体爆炸火焰传播速度研究装置和方法 Pending CN112816653A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110226048.4A CN112816653A (zh) 2021-03-01 2021-03-01 预混可燃气体爆炸火焰传播速度研究装置和方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110226048.4A CN112816653A (zh) 2021-03-01 2021-03-01 预混可燃气体爆炸火焰传播速度研究装置和方法

Publications (1)

Publication Number Publication Date
CN112816653A true CN112816653A (zh) 2021-05-18

Family

ID=75862648

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110226048.4A Pending CN112816653A (zh) 2021-03-01 2021-03-01 预混可燃气体爆炸火焰传播速度研究装置和方法

Country Status (1)

Country Link
CN (1) CN112816653A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113325034A (zh) * 2021-07-08 2021-08-31 安徽理工大学 一种煤矿瓦斯煤尘爆炸的试验系统及试验方法
CN116298179A (zh) * 2023-03-29 2023-06-23 北京航天试验技术研究所 测试含能材料气流冲击感度的装置及方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113325034A (zh) * 2021-07-08 2021-08-31 安徽理工大学 一种煤矿瓦斯煤尘爆炸的试验系统及试验方法
CN116298179A (zh) * 2023-03-29 2023-06-23 北京航天试验技术研究所 测试含能材料气流冲击感度的装置及方法

Similar Documents

Publication Publication Date Title
CN112816653A (zh) 预混可燃气体爆炸火焰传播速度研究装置和方法
Peters et al. Structure and similarity of nitric oxide production in turbulent diffusion flames
CN108535446B (zh) 管道瓦斯爆炸引起沉积煤尘二次爆炸的实验装置与方法
CN111255723B (zh) 一种防爆风机防爆性能的试验系统和方法
CN110319982A (zh) 基于机器学习的埋地式燃气管线泄漏判断方法
Pera An experimental database dedicated to the study and modelling of cyclic variability in spark-ignition engines with LES
CN103033606A (zh) 一种测定半敞开空间可燃气体最小点火能的方法
Beerer et al. An experimental ignition delay study of alkane mixtures in turbulent flows at elevated pressures and intermediate temperatures
CN214473306U (zh) 预混可燃气体爆炸火焰传播速度研究装置
Li et al. Prediction of equivalence ratio in pulse combustor from ion current amplitude spectrum
CN117233210B (zh) 一种聚乙烯粉尘爆炸危险性预测方法
Jing et al. Investigation on the characteristics of single-phase gas explosion and gas-coal dust coupling explosion in bifurcated tubes
Kypraiou et al. Proper orthogonal decomposition analysis of a turbulent swirling self-excited premixed flame
Cooper et al. Effect of deflagration-to-detonation transition on pulse detonation engine impulse
Versailles et al. Measurements of the reactivity of premixed, stagnation, methane-air flames at gas turbine relevant pressures
Chan Measurement of concentrations of transient gases using a conventional NDIR analyser
CN111006870A (zh) 一种定容燃烧传热装置
Yatsufusa et al. Development of measurement technic on propagating flame by densely installed ion-probes
CN115048853A (zh) 一种基于bp神经网络的甲烷的jwl状态方程参数标定方法
Wang et al. Study on noise-vibration coupling characteristics of premixed methane–air flame propagation in a tube with an acoustic absorption material
Licong et al. Numerical simulation of shock wave structure in gas explosion
CN116562193B (zh) 旋转爆震发动机燃烧效率分析方法及系统
Boeck et al. A physics-based model for explosions in vented vessel-pipe systems
Farabi et al. Simulation of Non-premixed Turbulent Natural Gas Combustion in Delft-Jet-in-Hot-Coflow Burner
Moinuddin Ansari Enhancing Performance and Reducing Emissions in Natural Gas Aspirated Engines through Machine Learning Algorithm

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