乙炔是工业上切割板材最常用的气体之一。乙炔极易燃烧、分解、爆炸，而且破坏力极大，爆炸事故屡屡发生，而乙炔十分不稳定，对它的研究比较少。有学者为了获得清晰的胞格，在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所示。

另外，建筑施工现场往往区域大，作业点面多，仅靠标准化设施防护及安全管理人员巡查还不够，还应采取视频监控形式，安全管理人员通过各点摄像头全场监控，以此做到全方位无死角盯防。

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

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

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

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

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

图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所示。

RVOT组患者右胸导联(V3R、V4R、V5R)R波振幅均明显小于LVOT组，而LVOT组患者V3R、V4R、V5R的S波振幅均明显小于RVOT组。见表2。

采用第二代带精英策略的非支配排序遗传算法（NSGA-II），对船舶储能系统复合储能装置进行设备选型和优化配置，确定磷酸铁锂电池和超级电容的型号及串并联个数。

图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，下同。

责任护士首先对患者个人情况以及心理状态进行全面评估,了解老人住院期间不适感,冠心病发作时给患者带来多种痛苦及负面情绪影响程度。尽量接纳老人合理要求,增加适当的心理疏导,责任护士每日与患者沟通时间不少于30 min。对于儿女不在身边照看的老人,安排专人照看,以免患者寂寞,产生焦虑、紧张的心理,养成孤僻习惯,封闭自己[3]。交流过程中,给予患者适当的鼓励以及积极的行为暗示,使老年患者保持积极乐观的心态,增强战胜疾病的勇气和信心。

NFC手机支付是一种把手机变成虚拟卡的卡模拟技术，是一种近场高频无线通信技术，允许设备间点对点非接触式数据交互，由非接射频RFID技术演变而来。NFC设备靠近非接读卡器后NFC设备进入工作状态，NFC设备与读卡器根据APDU协议建立交互，读卡器将敏感数据发给NFC设备的安全载体进行处理，NFC设备将响应数据发给读卡器。

由图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软件，编制图像识别程序，得到每个圆心的极坐标。相邻圆的两个圆心的距离为圆心距，对圆的半径和圆心距进行数字化处理，得到半径和圆心距的方差。通过半径和圆心距的方差，可以判断出胞格大小的离散程度和胞格分布的离散程度，从而判断出爆轰端面结构的规则程度。

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

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

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

表1 预混气在不同压力状况下的半径方差、圆心距方差及胞格数 Table 1 Radius variance, center distance variance and number of cells under different pressures

压力状况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所示。比较胞格数的变化及胞格半径的分布情况，可以判断出出胞格大小的集中程度。

通过教师的引导，学生可以在体会中发现，由于代词it和连词and的使用，（b）的表达较之于（a）更加自然和连贯。教师可顺便总结代词和连接词对语句连贯的衔接功能，使学生无意识的语篇知识转化为语篇意识。

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

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

“你的父亲，”她又说，“是一件外套。那个时候，他穿着外套来到我们家，就是睡觉也不脱下。一天夜里，我鼓足勇气伸出手在那件外套上一摸，发现里面什么也没有。直到多年之后我才弄清事情的真相。”[2]

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

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

普林斯顿大学的一个研究小组的专家给出了问题的答案。该研究小组正热衷于宇宙大爆炸理论的研究。早在1946年，美国物理学家伽莫夫提出了大爆炸宇宙模型，预言了作为大爆炸遗迹的电磁辐射背景存在的可能性。1953年，他指出宇宙如此古老以致变得异常寒冷，当年爆炸后残存的辐射温度可能只有5K。

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

2.2.2 相邻胞格中心距离

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

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

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

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

p/kPaV(r)/mmV(l)/mm5.323.169 89.983 95.742.446 17.080 47.452.145 45.511 710.570.674 82.189 612.710.331 21.218 9

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

p/kPaV(r)/mmV(l)/mm8.4011.779 926.149 010.505.065 910.481 214.001.910 56.992 615.571.905 45.861 616.231.354 04.797 7

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

2.3.1 胞格半径

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

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

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

2.3.2 相邻胞格中心距离

2.3.1 吸附时间 从图2可见，AB-8大孔树脂与提取液混合，在吸附2 h时吸附率最高为32.9%，时间延长吸附率没有升高。因此，吸附时间以2 h为宜。

同样，相邻胞格中心的距离为横坐标，以相应的胞格个数为纵坐标，得到相应的直方图。比较纵坐标的变化及相邻胞格中心距离的分布情况，给出胞格分布的均匀程度。其中，在初始压力为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。并且，随着压力的增大，胞格更加稳定且更加均匀。

AWPNM-CLQE算法中，步骤(3)挖掘Lk和Nk，步骤(4)挖掘强正关联规则，步骤(5)挖掘强负关联规则.过程Pretreatment()对文档集进行预处理；函数MiningAWPL1()挖掘L1；过程PruningNotQ()剪除不含查询词项Q的候选2_项集C2，以及小于minPR的正项集和小于minNR的负项集；过程MiningAWP ()挖掘Lk和Nk，并输出Lk和Nk.

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

参考文献：

[1] STREHLOW R A, LIAUGMINAS R, WATSON R H, et al. Transverse wave structure in detonations[J]. Symposium on Combustion, 1967,11(1):683-692.

[2] 喻健良,高远,闫兴清,等.高浓度氩气稀释气体爆轰波临界管径和临界间距关系[J].爆炸与冲击,2015,35(4):603-608.

YU Jianliang, GAO Yuan, YAN Xingqing, et al. Correlation between the critical tube diameter and annular interval for detonation wave in high-concentration argon diluted mixtures[J]. Explosion and Shock Waves, 2015,35(4):603-608.

[3] BIVOL G Y, GOLOVASTOV S V, GOLUB V V. Prechamber initiation of gaseous detonation in a channel[J]. Combustion Science and Technology, 2016,188(7):1165-1179.

[4] 张博,白春华.气相爆轰动力学特征研究进展[J].中国科学:物理学、力学、天文学,2014,44(7):665-681.

ZHANG Bo, BAI Chunhua. Research progress on the dynamic characteristics of gaseous detonation[J]. Scientia Sinica: Physica, Mechanica and Astronomica, 2014,44(7):665-681.

[5] 姜宗林,滕宏辉.气相规则胞格爆轰波起爆与传播统一框架的几个关键基础问题研究[J].中国科学:物理学、力学、天文学,2012,42(4):421-435.

JIANG Zonglin, TENG Honghui. Research on some fundamental problems of the universal framework for regular gaseous detonation initiation and propagation[J]. Scientia Sinica: Physica, Mechanica and Astronomica, 2012,42(4):421-435.

[6] 张博, LEE J H S, 白春华.高浓度氩气稀释对C2H2-2.5O2气体直接起爆临界能量影响的实验研究[J].高压物理学报,2012,26(1):55-62.

ZHANG Bo, LEE J H S, BAI Chunhua. Experiment investigation of the influence of highly argon dilution on the critical initiation energy for direct initiation of C2H2-2.5O2 mixtures[J]. Chinese Journal of High Pressure Physics, 2012,26(1):55-62.

[7] 韩学斌.管道内甲烷-空气预混气体火焰传播特性和结构特征的研究[D].北京:北京化工大学,2011.

[8] TROTSYUK A V, FOMIN P A, VASIL’EV A A. Numerical study of cellular detonation structures of methane mixtures[J]. Journal of Loss Prevention in the Process Industries, 2015,36:394-403.

[9] 徐晓峰,解立峰,彭金华,等.环氧丙烷-空气混合物爆轰波胞格结构的研究[J].爆炸与冲击,2004,24(2):158-162.

XU Xiaofeng, XIE Lifeng, PENG Jinhua, et al. Study on the detonation cellular structure of propylene epoxide-air[J]. Explosion and Shock Waves, 2004,24(2):158-162.

[10] 王昌建,徐胜利.直管内胞格爆轰的基元反应数值研究[J].爆炸与冲击,2005,25(5):405-416.

WANG Changjian, XU Shengli. Numerical study on cellular detonation in a straight tube based on detailed chemical reaction model[J]. Explosion and Shock Waves, 2005,25(5):405-416.

[11] 刘岩,武丹,王健平.低马赫数下斜爆轰波的结构[J].爆炸与冲击,2015,35(2):203-207.

LIU Yan, WU Dan, WANG Jianping. Structure of oblique detonation wave at low inflow Mach number[J]. Explosion and Shock Waves, 2015,35(2):203-207.

[12] TENG H, NG H D, KANG L, et al. Evolution of cellular structures on oblique detonation surfaces[J]. Combustion & Flame, 2014,162(2):470-477.

[13] SHEPHERD J E, MOEN I O, MURRAY S B, et al. Analyses of the cellular structure of detonations[J]. Symposium on Combustion, 1986,21(1):1649-1658.

[14] SHEPHERD J E, TIESZEN S R. Detonation cellular structure and image proces[C]∥The 21st International Symposium on Combustion. Munich, West Germany, 1986.

[15] LEE J J, FROST D L, LEE J H S, et al. Digital signal processing analysis of soot foils[C]∥KUHL A, LEYER J, BORISOV A, et al. Dynamic aspects of detonations. American Institute of Aeronautics and Astronautics, 1993:182-202.

[16] LEE J J, GARINIS D, FROST D L, et al. Two-dimensional autocorrelation function analysis of smoked foil patterns[J]. Shock Waves, 1995,5(3):169-174.

[17] ZHAO H J, LEE J H S, LEE J, et al. Quantitative comparison of cellular patterns of stable and unstable mixtures[J]. Shock Waves, 2016,26(5):621-633.

[18] 赵焕娟, LEE J H S, 张英华,等.边界条件对甲烷预混气爆轰特性的影响[J].爆炸与冲击,2017,37(2):201-207.

ZHAO Huanjuan, LEE J H S, ZHANG Yinghua, et al. Effects of boundary conditions on premixed CH4+2O2 detonation characteristics[J]. Explosion and Shock Waves, 2017,37(2):201-207.

[19] LEE J H S. The detonation phenomenon[M]. New York: Cambridge University Press, 2008.

[20] VOITSEKHOVSKII B V, MITROFANOV V V, TOPCHIYAN M E. Structure of the detonation front in gases (survey)[J]. Combustion, Explosion, and Shock Waves, 1969,5(3):267-273.

[21] 赵焕娟, LEE J H S, 张英华,等.爆轰波三波点擦除烟迹表面积碳机制[J].工程科学学报,2017,39(3):335-341.

ZHAO Huanjuan, LEE J H S, ZHANG Yinghua, et al. Precise mechanism of triple point passage removing soot on soot-coated surface[J]. Chinese Journal of Engineering, 2017,39(3):335-341.