火炸药学报    2018, Vol. 41 Issue (4): 399-403   DOI: 10.14077/j.issn.1007-7812.2018.04.014
0

引用本文  

王长健, 赵凤起, 孙美, 杨燕京, 古城辉. 基于透过率场分布的推进剂羽流烟雾面积测试方法[J]. 火炸药学报, 2018, 41(4): 399-403. DOI: 10.14077/j.issn.1007-7812.2018.04.014
WANG Chang-jian, ZHAO Feng-qi, SUN Mei, YANG Yan-jing, GU Cheng-hui. Smoke Area Measurement Method of Propellant Plumes Based on Transmittance Field Distributing[J]. Chinese Journal of Explosives & Propellants, 2018, 41(4): 399-403. DOI: 10.14077/j.issn.1007-7812.2018.04.014

基金项目

国家自然科学基金(No.201503163)

作者简介

王长健(1983-), 男, 高级工程师, 从事推进剂羽流特征信号评估工作。E-mail:18991898347@163.com

文章历史

收稿日期:2017-08-23
修回日期:2018-04-24
基于透过率场分布的推进剂羽流烟雾面积测试方法
王长健, 赵凤起, 孙美, 杨燕京, 古城辉     
西安近代化学研究所燃烧与爆炸技术重点实验室, 陕西 西安 710065
摘要: 为表征推进剂装药羽流烟雾的空间分布特征,提出了一种基于烟雾透过率分布的烟雾面积测试方法,即将火药装药燃烧流场烟雾的光学特性(透过率)与烟雾的物理特性(面积参数)同步进行测量,进而分析基于透过率的羽流烟雾分布区域的空间及时域特性。首先利用图像传感器对光学背景板前方的羽流烟雾图像进行采集,然后计算测试区域的烟雾透过率场分布,最后以特定透过率区间作为烟雾区域的判别条件,提取烟雾区域并计算烟雾区域面积;对Al-CMDB推进剂装药羽流烟雾在工作后0.6~2.5 s内的烟雾分布面积特性进行了分析。结果表明,透过率区间在[0,30%]的烟雾区域在测试区域中部的烟雾面积明显大于左、右两个区域;在发动机开始工作后0.60~0.76 s内,透过率在[30%,70%]区间的烟雾区域面积明显大于其他两个透过率区间的面积;在0.76 s之后,透过率在[70%,100%]区间的烟雾面积逐步增大;在1.6 s后不存在透过率小于70%的烟雾区域。该方法可有效测试特定透过率区间内的羽流烟雾空间分布以及其面积随时间变化的特性,也可用于发射药身管武器膛口烟雾面积特性分析。
关键词: 透过率场分布     推进剂装药     羽流烟雾     烟雾透过率     低特征信号推进剂    
Smoke Area Measurement Method of Propellant Plumes Based on Transmittance Field Distributing
WANG Chang-jian, ZHAO Feng-qi, SUN Mei, YANG Yan-jing, GU Cheng-hui     
Science and Technology on Combustion and Explosion Laboratory, Xi'an Modern Chemistry Research Institute, Xi'an 710065, China
Abstract: To characterize the spatial distribution characteristics of plume smoke of propellant charge, a smoke area measurement method based on smoke transmittance distribution was proposed. The opitcal property (transmittance) and space characteristic (area) were measured simultaneously and analyzed. First, the plume smoke images in the front of optical background plate were collected by image sensor. Then the transmittance field distribution in the test area was calculated. Finally, the specific transmittance interval is used as the discriminative condition of the smoke region and the area of smoke region was extracted and the area of smoke region was calculated. The smoke distribution area characteristics of plume smoke of Al-CMDB propellant charge during 0.6-2.5 s after working was analyzed. The results show that the smoke area of transmittance in the[0, 30%] interval in the middle of the test area is bigger than that of left region and right region. In 0.60-0.76 s after the start of the engine work, the smoke area of transmittance in the[30%, 70%] interval is obviously bigger than that of other two transmittance intervals; after 0.76 s, the smoke area of transmittance in the[70%, 100%] interval is gradually increased. After 1.6 s, there is no smoke area with the transmittance of less than 70%. This method can effectively test the spatial distribution of plume smoke in the range of specific transmittance and the characteristics of its area with time change. It can also be used to analyze the smoke area characteristics of the muzzle of the gun propellant barrel weapon.
Keywords: transmittance field distributing     propellant charge     plume smoke     smoke transmittance     low signature propellant    
引言

低特征信号推进剂以及洁净发射药是火药发展的重要方向之一。火药装药在工作过程中产生的流场烟雾,会显著降低武器系统的制导和隐身性能[1-4]。一方面,火药燃烧产生的烟雾会暴露武器系统的位置及运动轨迹;另一方面,流场烟雾对于穿过其中的制导信号有衰减、散射、吸收等多种作用,导致制导信号强度降低,严重时甚至失效。此外,身管武器中发射药燃烧产生的烟雾还会对操作人员或视频制导系统的二次瞄准造成干扰。

火药装药燃烧流场烟雾对武器系统制导信号的干扰程度,不仅与烟雾浓度有关,还和烟雾面积、厚度以及制导信号穿过烟雾的路径有关。目前,国内对火药装药流场烟雾的测试方法多是对流场烟雾透过率的测量,很少涉及到对于流场烟雾场形态的分布、烟雾区域大小等物理参数的研究[5-8]。在火灾预防、发烟剂烟雾特性测试等领域中,通常依据烟雾的形态特征以及灰度阈值对视频或图像中的烟雾区域进行提取并计算烟雾区域面积[9-11]。在这些方法所得到的结果中,烟雾区域内的烟雾透过率范围及分布是不明确的,也无法确定区域内烟雾场对制导信号的影响程度。

本研究提出一种基于透过率场分布的烟雾面积测试方法,将火药装药燃烧流场烟雾的光学特性(即透过率)与烟雾的物理特性(面积参数)同步进行测量,进而分析基于透过率的羽流烟雾分布区域的空间及时域特性。

1 实验 1.1 基本原理

为完成基于透过率的烟雾面积测试,测试装置主要由光学背景板、图像传感器、工控机组成[8]。其中光学背景板提供灰白相间条纹的光学背景,图像传感器对光学背景前方的流场烟雾图像进行采集[12-13]。在测试过程中,测试系统与测试对象(推进剂装药或发射药装药)的现场布置俯视图如图 1所示。光学背景板位于烟雾区域的一侧,图像传感器位于烟雾区域另一侧。

图 1 烟雾面积测试系统原理(俯视图) Figure 1 Principle of smoke area measure system(top view)

烟雾透过率的计算基于朗伯贝尔定律。本方法利用图像传感器对火药装药燃烧流场烟雾进行图像采集。假定羽流烟雾自身为非光源,无光信号辐射的情况下,图像中烟雾区域某点的成像灰度是该区域烟雾透射的光信号和烟雾反射环境光信号这两者在图像传感器的叠加,可表示为:

$ {{I}_{0}}={{I}_{\text{a0}}}\cdot \tau +{{I}_{\text{b0}}} $ (1)

式中:I0为图像传感器接收的信号光强;Ia0为光学背景的信号光强;τ为烟雾的光透过率;Ib0为烟雾反射、折射其他光源的信号光强。

为消除烟雾反射折射其他光源引入的光噪声,背景板设计为灰白相间条纹[9],对于不同灰度的光学背景,传感器接收到的光信号强度可表示为:

$ \left\{ \begin{align} &{{I}_{0}}={{I}_{\text{a0}}}{{\tau }_{0}}+{{I}_{\text{b0}}}~ \\ &{{I}_{1}}={{I}_{\text{a1}}}{{\tau }_{1}}+{{I}_{\text{b1}}} \\ \end{align} \right. $ (2)

一般情况下,烟雾场分布是不均匀的,且随时间发生变化。在特定时刻,当所表征的烟雾场区域足够小时,可以近似认为相邻两个最小单元内,烟雾透过率相同,并且其反射、折射其他光源信号的强度相同,即τ0=τ1,且Ib0=Ib1。在传感器的线性响应区间内,图像像素灰度与该点光信号强度成正比,因此该点的光透过率可表示为:

$ \tau =~\frac{{{I}_{0}}-{{I}_{1}}}{{{I}_{\text{a0}}}-{{I}_{\text{a1}}}} $ (3)

通过公式(3)可获得烟雾图像中烟雾场内任意一点(x, y)处的烟雾透过率。在实际应用中,还可以对最小处理单元的烟雾透过率变化建立一些更为复杂的模型。在此基础上,进行图像中烟雾透过率分布场的几何特征研究,可以获得现实烟雾的透过率分布场对应的空间特性。假设图像烟雾场内任意两点的坐标分别为(x1, y1)、(x2, y2),则现实烟雾场对应两点在成像平面投影长度的函数关系可表示为:

$ L=k\sqrt{{{\left( {{x}_{1}}-{{x}_{2}} \right)}^{2}}+{{\left( {{y}_{1}}-{{y}_{2}} \right)}^{2}}} $ (4)

式中:L为物理空间内两点的距离;k为图像比例尺。

将烟雾图像所对应的透过率矩阵,以U(τaττb)为阈值进行二值化处理(其中τa, τb为常数),即可实现对烟雾图像中特定透过率烟雾区域的识别。该过程可表示为:

$ B\left( i, j \right)=\left\{ \begin{align} &0, \tau \left( i, j \right)>{{\tau }_{\text{b}}}\cup \tau \left( i, j \right)<{{\tau }_{\text{a}}}~ \\ &1, ~\ \ \ \ \ \ \ \ \ {{\tau }_{\text{a}}}\le T\left( i, j \right)\le {{\tau }_{\text{b}}} \\ \end{align} \right. $ (5)

式中:B(i, j)为图像二值化操作后的数值;τ(i, j)为图像中(i, j)像素点对应的烟雾透过率。

通过统计二值图像中数值为1的像素点数量,利用图像的比例尺即可计算烟雾图像中烟雾区域的面积。该过程可表示为:

$ S={{k}^{2}}\cdot N $ (6)

式中:S为图像中特定透过率阈值的烟雾区域面积;k为比例尺;N为满足阈值条件的像素点数量。

1.2 实验设计

测试样品为Al-CMDB推进剂(含铝粉的复合改性双基推进剂),单孔管状药,装药方式为自由装填。发动机为Φ50 mm标准试验缩比发动机,工作压力设计为(7±0.5) MPa,工作时间为0.72 s。

所采用的光学背景板尺寸为2.4 m×0.8 m,条纹宽度为0.2 m,光学背景板左侧距发动机喷口水平距离为2 m,发动机轴线位于背景板高度0.4 m处。图像传感器为加可见光滤波片的日本Photron FASTCAM高速摄像机(响应光谱为0.4~0.7 μm,帧速60 f/s)。外场环境,环境温度(30±10) C°,风力小于3级,少云。

2 结果与讨论 2.1 特定时刻烟雾空间分布的特征分析

通过对测试获得图像序列中的单帧图像进行处理,可获得该时刻的烟雾透过率空间分布数据。图 2(a)为Al-CMDB推进剂装药开始工作第0.62 s时的烟雾可见光透过率分布伪彩图,图 2(b)为该时刻烟雾可见光透过率的直方分布图,图 2(c)为透过率在[0, 30%]区间内的烟雾区域分布图,图 2(d)为透过率在[40%, 70%]区间内的烟雾区域分布图。

图 2 烟雾场透过率分布图 Figure 2 Smoke field transmittance distribution map

图 2(a)可知,在该时刻,羽流烟雾分布并不均匀,浓度高的烟雾主要分布在测试区域的上半部分,该区域透过率偏低;由图 2(c)可知,透过率小于30%的区域基本集中在测试区域的右上角。

图 2(b)可知,在该时刻,测试区域的烟雾透过率大部分分布在[40%, 80%]区间内;通过对图 2(d)分析可知,透过率区间在[40%, 70%]的区域主要分布在测试区域的中部和下部。

2.2 烟雾空间分布的时域特征

发动机开始工作后的测试区域内基于透过率区间、水平区域以及垂直区域上划分的烟雾面积-时间曲线如图 3所示。其中对透过率在[0, 30%]、[30%, 70%]、[70%, 100%]区间的烟雾面积分别统计,结果如图 3(a)所示;将测试区域在垂直方向上等分为上部、中部、下部3块区域,统计透过率在[0, 30%]区间内的烟雾区域面积,结果如图 3(b)所示;将测试区域在水平方向上等分为左、中、右3块区域,统计透过率在[0, 30%]区间内的烟雾区域,结果如图 3(c)所示。

图 3 烟雾面积—时间曲线 Figure 3 Smoke area—time curves

结合试验烟雾图像进行分析可以发现,在发动机工作过程中羽流火焰辐射强度大,背景板及烟雾本身均反射明显的火焰光信号,部分区域在图像传感器中成像灰度饱和,对测试结果中烟雾透过率准确度的影响较大。如需克服羽流火焰对测试结果的影响, 一方面需要增大图像传感器对光强的响应范围;另一方面也可调整背景板的位置, 使其远离羽流火焰, 减少火焰辐射对测试装置的影响。在图 3数据中,t=0.6 s之前的数据可信度较差,不纳入讨论范围。

图 3(a)可以看出,在0.6~0.76 s时,透过率在[30%,70%]内的烟雾区域面积大于其他两个区间内的烟雾面积;当t > 0.76s时,透过率在[70%,100%]区间的烟雾面积大于其他两个区间内的烟雾面积;当t=1.6 s后,透过率在70%以上的烟雾区域基本占据整个测试区域。相当于在发动机工作结束0.88 s后,测试区域内的羽流烟雾基本完全消散。

图 3(b)可以看出,在发动机工作后的0.6~0.93 s时间内,透过率在[0, 30%]区间的烟雾面积整体随时间增加而逐步减小,其中在测试区域上部的烟雾面积始终大于其他两个区域,透过率在[0, 30%]区间烟雾面积在测试区域中部及下部基本相同,且在0.72 s之后基本都降为0。这是由于发动机工作结束后羽流流速减弱,羽流烟雾受热空气影响,有明显的上漂现象,这一点在烟雾图像中也有明显体现。

图 3(c)可以看出,在发动机工作后的0.6~0.93 s时间内,透过率在[0, 30%]区间的烟雾面积整体随时间增加而逐步减小,其中测试区域中部的烟雾面积相对最大,烟雾区域左端的烟雾面积相对最小。这说明羽流烟雾在向前运动的过程中有一个扩散并且稀释的过程,在中部区域(距喷口2.8~3.6 m范围),烟雾扩散起主要作用,透过率在[0, 30%]区间的烟雾面积有一定程度增大;在右部区域(距喷口3.6~4.4 m范围),烟雾稀释起主要作用,烟雾面积较中部区域烟雾面积有一定程度的减小。

3 结论

(1) Al-CMDB推进剂装药羽流烟雾分布并不均匀。从水平方向来看,透过率区间在[0,30%]的烟雾区域在测试区域中部的分布面积明显大于左、右两个区域;在垂直方向来看,透过率区间在[0,30%]的烟雾区域在上部的分布面积明显大于中、下两个区域。

(2) Al-CMDB推进剂装药羽流烟雾在时域分布上具有明显特征。在0.6~0.76 s内,透过率在[30%,70%]区间的烟雾区域面积明显大于其他两个透过率区间的面积;在0.76 s之后,透过率在[70%,100%]区间的烟雾面积逐步增大;在1.6 s后不存在透过率小于70%的烟雾区域。

(3) 该方法能有效测定特定透过率区间的推进剂装药羽流烟雾面积及其空间分布特征,但是在靠近羽流火焰区域,测试结果的精度会受到较大影响。该方法也可用于发射药身管武器膛口烟雾面积特性分析。


参考文献
[1]
李上文, 赵凤起, 徐司雨. 低特征信号固体推进剂技术[M]. 北京: 国防工业出版社, 2013.
[2]
莫红军, 王宁飞. 固体火箭发动机羽烟特征信号的分类及测试评估概述[J]. 固体火箭技术, 2004, 27(3): 229–232.
MO Hong-jun, WANG Ning-fei. Classification and assessment of solid propellant rocket exhaust signature[J]. Journal of Solid Rocket Technology, 2004, 27(3): 229–232. DOI:10.3969/j.issn.1006-2793.2004.03.016
[3]
姚禄玖, 高钧麟, 肖凯涛, 等. 烟幕理论与测试技术[M]. 北京: 国防工业出版社, 2004.
[4]
王宏, 刘桂生, 孙美, 等. 固体火箭发动机羽流特征信号检测技术研究[J]. 固体火箭技术, 2001, 24(2): 64–66.
WANG Hong, LIU Gui-sheng, SUN Mei, et al. Study on the measurement of characteristic signals of plume in rocket motor[J]. Journal of Solid Rocket Technology, 2001, 24(2): 64–66. DOI:10.3969/j.issn.1006-2793.2001.02.016
[5]
Marbach G, Loepfe M, Brupbacher T. An image processing technique for fire detection in video images[J]. Fire Safety Journal, 2006, 41(4): 285–289. DOI:10.1016/j.firesaf.2006.02.001
[6]
乔杨, 徐熙平, 卢常丽, 等. 机器视觉距离目标尺寸自动标定测量系统研究[J]. 兵工学报, 2012, 33(6): 759–763.
QIAO Yang, XU Xi-ping, LU Chang-li, et al. Study on automatic long-distance dimension measurement system based on machine vision[J]. Acta Armamentarii, 2012, 33(6): 759–763.
[7]
郑日恒, BrayK N C, RoggB, 等. CCD互相关粒子图像测速及其在对抗流场测量中的作用[J]. 推进技术, 2000, 21(6): 83–88.
ZHENG Ri-heng, Bray K N C, Rogg B, et al. CCD Cross-correlation particle image velocimetry and its application to counterflow flowfield measurements[J]. Journal of Propulsiton Teachnology, 2000, 21(6): 83–88. DOI:10.3321/j.issn:1001-4055.2000.06.022
[8]
牛群峰, 陈生生, 王莉. 基于虚拟仪器的火焰尺寸测量系统设计[J]. 计测技术, 2016, 36(1): 25–29.
NIU Qun-feng, CHEN Sheng-sheng, WANG Li. Design of flame size measurement system based on virtual instrument[J]. Metrology & Measurement Technology, 2016, 36(1): 25–29.
[9]
李宋, 吴文全. 应用CCD对浓度场的全场瞬时测量[J]. 物理实验, 2004, 24(8): 15–18.
LI Song, WU Wen-quan. Instantaneous measurement of full concentration fields based on CCD[J]. Imentation, 2004, 24(8): 15–18. DOI:10.3969/j.issn.1005-4642.2004.08.004
[10]
章惠, 张露. 低浓度烟雾检测方法研究与仿真[J]. 计算机仿真, 2012, 29(5): 254–257.
ZHANG Hui, ZHANG Lu. Research and simulation of low concentration smoke detection method[J]. Computer Simulation, 2012, 29(5): 254–257. DOI:10.3969/j.issn.1006-9348.2012.05.062
[11]
虞燕凤. 基于视频图像的火灾烟雾检测研究[D]. 大连: 大连海事大学, 2013. http://cdmd.cnki.com.cn/Article/CDMD-10151-1013204643.htm
[12]
王长健, 孙美, 魏志刚, 等. 基于视频图像分析的推进剂羽流烟雾透过率测量技术[J]. 火炸药学报, 2016, 39(5): 79–83.
WANG Chang-jian, SUN Mei, WEI Zhi-gang, et al. Research on measurement technology of plume-smoke transmittance of solid propellant based on video image analysis[J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao), 2016, 39(5): 79–83.
[13]
王长健, 张克旺, 孙美, 等. 基于图像传感器的推进剂发烟量测试方法[J]. 火炸药学报, 2017, 40(1): 97–100.
WANG Chang-jian, ZHANG Ke-wang, SUN Mei, et al. Propellant smoke density test method based on image sensor[J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao), 2017, 40(1): 97–100.