Ar稀释C2H2+2.5O2预混气高频爆轰的端面结构*

赵焕娟,伯玉兰,张英华,严屹然

(北京科技大学资源与工程学院,北京 100083)

乙炔是工业上切割板材最常用的气体之一。乙炔极易燃烧、分解、爆炸，而且破坏力极大，爆炸事故屡屡发生，而乙炔十分不稳定，对它的研究比较少。有学者为了获得清晰的胞格，在C2H2-O2中加入Ar[1-2]，研究Ar对预混气的胞格结构的影响，但大多是定性分析[3]，且限于表征壁面胞格[4]。在我国，学者们在气相胞格爆轰传播现象方面进行了大量研究[5-6]，且大多采用烟膜法进行[7]，在胞格结构表征爆轰、胞格形成机理及其量化规律方面取得了较大进展[8-12]。

Shepherd等[13-14]、Lee等[15-16]都尝试过使用数字图像处理技术进行定量的不稳定度分析，通过数字图像处理技术直接扫描烟膜来获得端面结构。由于非均匀烟灰沉积物产生的误差，导致“灰度不均匀”，所以数字化处理的灰度图像误差较大，因此，必须对灰度进行处理。Lee等[15-16]、赵焕娟等[17-18]求助于手绘记录扫描，得到可以进行数字处理的图片。在进行烟膜手绘记录时，需要加入少量人为主观性判断，以消除少量的杂散线和实验中的误操作导致的线条。

环绕在管内壁的烟膜上仅记录了爆轰管壁处的三波点经过的路径，当前对爆轰结构的研究大多集中在管壁上螺旋爆轰结构上[19]，对其内部发展规律研究较少。Voitsekhovskii等[20]在螺旋爆轰面单头、双头和四头螺旋爆轰的横波运动方面做了大量研究。对于更高频的螺旋爆轰，其横波运动穿过爆轰面的表征方法非常困难，需要实验结果的支持，如何清晰地记录爆轰面的结构是一个关键问题[21]。

本文中，使用在管道末段设置烟熏玻璃的方法，记录C2H2+2.5O2+70%Ar和C2H2+2.5O2+85%Ar预混气高频螺旋爆轰端面结构，通过自编程序对端面记录图像进行数字化处理，分析表征高频螺旋爆轰端面胞格结构。

1 实 验

1.1 方法

实验装置主要包括4个部分，分别是爆轰管道系统、数据采集系统、点火系统和充配气系统，如图1所示。

爆轰管道系统由驱动段和测试段组成。驱动段为长1 000 mm、内径63.5 mm、壁厚10 mm的钢管，测试段为聚碳酸酯透明管。在驱动段左侧法兰端面处安装点火塞，并在靠近点火塞位置布置Shchelkin螺旋增爆器，用来加速爆轰的形成过程，从而在驱动段内形成稳态爆轰。

将测试主体管两端密封，检查测试主体管的密封性，将测试段的一端通过法兰密封，在测试段的另一端安装开槽法兰，将烟熏玻璃(厚4～6 mm)置于法兰末端，保证放置过程中烟熏玻璃无旋转。

制备化学当量比为C2H2+2.5O2+70%Ar和C2H2+2.5O2+85%Ar的预混气，所有的爆炸性气体都已预先在高压瓶由分压的方法制备，通过控制面板将不同气体分别安全输入到充配气系统的预混气高压瓶中，充罐后预混气放置24 h以上，以确保罐内不同气体均匀混合。

将测试段抽真空，前后分别充入引爆气C2H2-O2及实验气体C2H2+2.5O2+70%Ar和C2H2+2.5O2+85%Ar点火，爆轰结束后将管内抽真空，平稳缓慢放入大气，以保证不在烟熏玻璃上形成冲刷痕迹。打开测试段末端法兰，取出烟熏玻璃后进行扫描。

图1 直径63.5 mm爆轰管道的结构简图Fig.1 Schematic diagram of detonation tube with diameter of 63.5 mm

1.2 结果

在初始压力高于5 kPa时，获得自持爆轰传播。为了能够容易区分螺旋头的个数，对于70%Ar的预混气，通过压力传感器把初始压力分别定为5.32、5.74、7.45、10.57 kPa，得到的端面烟膜图如图2所示。

图2 C2H2+2.5O2+70%Ar预混气的烟膜记录Fig.2 Smoked foils records of premixed C2H2+2.5O2+70%Ar

对于含85%Ar预混气，通过压力传感器把初始压力分别定为3.65、14、15.57、16.23 kPa，得到的端面烟膜图如图3所示。

图3 C2H2+2.5O2+85%Ar预混气的烟膜记录Fig.3 Smoked foils records of premixed C2H2+2.5O2+85%Ar

2 结果分析

2.1 端面烟膜数字化处理

为了获得清晰的端面结果，取C2H2+2.5O2+70%Ar预混气初始压力分别为5.32、5.74、7.45、10.57 kPa的端面烟膜图(见图2)，使用图像处理软件Photoshop，通过调整色差和对比度，利用钢笔工具进行描点，然后描边轨迹，填充椭圆工具选择出来的圆，画出三波轨迹，如图4所示。阴影为管道内径范围，直径为63.5 mm，下同。

取C2H2+2.5O2+85%Ar预混气初始压力分别为3.65、14.00、15.57、16.23 kPa的端面烟膜图(见图3)，处理的三波轨迹如图5所示。

由图4～5可知，随着压力的增大，端面结构都变得更加复杂，得到的封闭图形的尺寸也都在变小。胞格尺寸代表波中能量的大小，爆轰的传播有波峰面燃烧的热量和爆轰波的能量，因为爆轰初始压力的增大，分子初始获得的能量多，更多活化分子参与反应，提高分子间反应速率，可以形成更多胞格，表现在胞格尺寸就是变小。

但是，可以发现两种预混气的胞格规则程度不同，需要定量分析。这是因为随着Ar浓度越高，对C2H2稀释越强，爆轰得到稳定的多头爆轰需要的压力更大，所以85%Ar的预混气形成单头爆轰的压力更高。由于单头爆轰、双头爆轰及四头爆轰所给出的胞格数量太少，不具有参考性，所以主要研究多头爆轰随气体浓度及初始压力的关系。

图4 C2H2+2.5O2+70%Ar预混气的端面示意图Fig.4 End-on result of premixed C2H2+2.5O2+70%Ar

图5 C2H2+2.5O2+85%Ar预混气的端面示意图Fig.5 End-on result of premixed C2H2+2.5O2+85%Ar

为了说明多头爆轰端面结构与气体浓度和气体初始压力的关系，将端面结构中每个封闭的图形进行圆的拟合，得到对应封闭图形的圆心及半径。常见的圆曲线的拟合方法主要有平均值法、加权平均法及最小二乘法。平均值法与加权平均法虽然计算过程简单、难度较低，但是二者都只能作到在一定程度上的近似，在离散点分布均匀或者精度要求不高的情况下可以使用；而最小二乘法理论严密、精度最高，拟合效果最佳，是目前圆曲线拟合的主要方法，本文中采用最小二乘法进行圆的拟合。通过MATLAB软件，编制图像识别程序，得到每个圆心的极坐标。相邻圆的两个圆心的距离为圆心距，对圆的半径和圆心距进行数字化处理，得到半径和圆心距的方差。通过半径和圆心距的方差，可以判断出胞格大小的离散程度和胞格分布的离散程度，从而判断出爆轰端面结构的规则程度。

根据爆轰端面实验结果，给出了几种标准图形，如图6所示。根据绘制的端面结果示意图，用MATLAB程序进行图像识别并数字化处理，得到相应的胞格半径方差、相邻胞格圆心距方差及胞格数，见表1。

由表1可知，随着胞格数的增大，胞格半径方差减小，相邻胞格圆心距方差减小，端面胞格尺寸更均匀且分布更规则。程序识别得到的结果和图片结果相符，说明该程序可靠。

图6 预混气端面处理结果示意图Fig.6 Schematic diagram of end-on results of premixed gas

压力状况V(r)/mmV(l)/mmn(a)4.032 65.944 034(b)3.140 55.091 739(c)1.107 02.098 077(d)0.255 10.561 7119

2.2 C2H2+2.5O2+70%Ar端面结构处理

2.2.1胞格半径

以胞格半径尺寸为横坐标，以胞格数量为纵坐标，得到相应的直方图如图7所示。比较胞格数的变化及胞格半径的分布情况，可以判断出出胞格大小的集中程度。

图7 C2H2+2.5O2+70%Ar预混气爆轰端面结构半径等间距直方图Fig.7 Histograms with same radius distance

从直方图可以看出，随着压力的增大，半径大小逐渐变小且更加集中，从纵坐标可以看出胞格的个数逐渐增多，说明随着压力的增大，端面结构胞格大小趋于稳定，最后稳定在1.5～3.0 mm之间。

图8 C2H2+2.5O2+70%Ar预混气爆轰 端面结构圆心距等间距直方图Fig.8 Histograms with same center distance

方差作为随机变量偏离程度的参考量，可以用来量化胞格半径的偏离程度。所以使用统计学方差公式来计算间距数据，可得到半径方差，见表2。

通过方差的小大可以看出来，随着压力的增大，方差逐渐变小，且趋于稳定。说明压力越大，半径的离散度越小，且逐渐趋于稳定，与直方图得出的结论一致。

2.2.2相邻胞格中心距离

以相邻胞格中心距离为横坐标，以胞格个数为纵坐标，可得到相应的直方图。比较纵坐标的变化及圆心距的分布情况，可给出胞格分布的均匀程度。其中，在初始压力为10.57 kPa下得到的直方图，如图8所示。

HIGHLIGHT：Espionage and fraud in cyber is not an armed conflict,says System Expert’s Jonathan Gossels.

可以看出，随着压力的增大，圆心距逐渐趋于集中，且个数逐渐变多，说明随着压力的增大，胞格分布趋于稳定。得到的不同初始压力下圆心距方差，见表2。

比较方差的小大，可以发现，方差随初始压力的增大逐渐变小。说明压力越大圆心距的离散度越小，胞格的分布越均匀，与直方图得出的结论一致。

表2 C2H2+2.5O2+70%Ar在不同压力下 的半径方差和圆心距方差Table 2 Radius variance, center distance variance of C2H2 +2.5O2+70% Ar under different pressures

表3 C2H2+2.5O2+85%Ar在不同压力下 的半径方差和圆心距方差Table 3 Radius variance, center distance variance of C2H2 +2.5O2+85% Ar under different pressures

2.3 C2H2+2.5O2+85%Ar端面结构处理

2.3.1胞格半径

同样，以胞格半径为横坐标，以相应胞格个数为纵坐标，得到相应的直方图。比较纵坐标的变化及半径的分布情况，给出胞格大小的集中程度。其中，在初始压力为16.23 kPa下得到的直方图，如图9所示。

由直方图可以看出，随着压力的增大，半径越来越小，且圆的个数越来越多。说明随着压力的增大，胞格大小更加稳定，胞格半径集中在2～5 mm。同样，可以得到半径的方差，见表3。

通过方差的小大变化，可以看出，随着压力的增大，方差逐渐变小。说明压力越大，半径的离散度越小，胞格大小更加稳定，与直方图得出的结论一致。

2.3.2相邻胞格中心距离

同样，相邻胞格中心的距离为横坐标，以相应的胞格个数为纵坐标，得到相应的直方图。比较纵坐标的变化及相邻胞格中心距离的分布情况，给出胞格分布的均匀程度。其中，在初始压力为16.23 kPa下得到的直方图，如图10所示。

从直方图可以发现，随着压力的增大，圆心距逐渐趋于集中，且个数逐渐变多。说明随着压力的增大，胞格分布趋于稳定。同样，可以得到圆心距的方差，见表3。

图9 C2H2+2.5O2+85%Ar预混气爆轰端面结构 半径等间距直方图Fig.9 Histograms with same radius distance

图10 C2H2+2.5O2+85%Ar预混气爆轰端面结构 圆心距等间距直方图Fig.10 Histograms with same center distance

通过方差的小大可以看出来，随着压力的增大，方差逐渐变小。说明压力越大圆心距的离散度越小，胞格的分布越均匀，与直方图得出的结论一致。

2.4 不同Ar浓度的对比分析

不同预混气体的活化能受到初始压力的影响较大，导致不同初始压力下形成的端面胞格数量差别太多，样本数量不够精确，所以压力较小时得到的胞格数量样本不够。因此，尝试分析胞格数近似的不同预混气的端面结构的稳定性。将C2H2+2.5O2+70%Ar及C2H2+2.5O2+85%Ar预混气在不同初始压力下形成的端面结果，通过编程处理得到的胞格数及半径方差与圆心距方差进行对比分析，如图11～12所示。

由图可知，胞格数相近的C2H2+2.5O2预混气，Ar浓度越高，半径方差及圆心距方差越小。说明胞格的稳定性随着Ar浓度的升高而增强。

图11 不同浓度预混气胞格半径方差随胞格数的变化Fig.11 Variation of radius variance with number of cell for different Ar concentrations

图12 不同浓度预混气相邻胞格圆心距方差随胞格数的变化Fig.12 Variation of center distance variance with number of cell for different Ar concentrations

3 结 论

(1)在内径63.5 mm圆形管道进行爆轰实验，得到不同Ar浓度的乙炔预混气在不同初始压力下形成的端面结构。

(2)使用数字化图像技术，处理实验得到了端面烟膜记录，获得不同初始条件下的C2H2+2.5O2+70%Ar和C2H2+2.5O2+85%Ar预混气的爆轰端面结构。其中70%Ar预混气得到的胞格半径稳定在1.5～3.0 mm，半径方差在0.331 2～3.169 8 mm，圆心距方差在1.218 9～9.983 9 mm；85%Ar预混气得到的胞格半径稳定在2.0～5.0 mm，半径方差在1.354 0～11.779 9 mm，圆心距方差在4.797 7～26.149 0 mm。并且，随着压力的增大，胞格更加稳定且更加均匀。

(3)在胞格数近似的情况下，85%Ar预混气得到的半径方差及圆心距方差比70%Ar的方差小，Ar浓度越高的预混气得到的胞格结构越稳定。