黃 勇，解立峰，張紅偉，魯長(zhǎng)波，安高軍，熊春華，陳 群

（1. 南京理工大學(xué)化工學(xué)院，江蘇 南京 210094；2. 常州大學(xué)環(huán)境與安全工程學(xué)院，江蘇 常州 213164；3. 中央軍委后勤保障部油料研究所，北京 102300；4. 常州大學(xué)石油化工學(xué)院，江蘇 常州 213164）

柴油被廣泛用作為大型車(chē)輛、艦船和武器裝備的燃料，一旦車(chē)輛發(fā)生交通事故或者裝甲裝備遭受炮火襲擊，都可能會(huì)引爆柴油造成人員傷亡[1]。因此，美國(guó)陸軍坦克研發(fā)工程中心和俄羅斯國(guó)防部第25國(guó)立油料化學(xué)裝備研究所先后開(kāi)展了防火柴油技術(shù)的研究[2-3]。防火柴油是含有微乳化成分的含水柴油，利用自身的防火性能使得泄漏燃料不易燃燒形成地面池火。近年，魯長(zhǎng)波等在美、俄防火柴油中添加高分子聚合物防霧劑，運(yùn)用防霧劑對(duì)燃料拋撒霧化的抑制作用開(kāi)發(fā)出一種具有防火抑爆功能的新型微乳化柴油[4]；他又通過(guò)理化性能分析、儲(chǔ)存安定性評(píng)定、實(shí)驗(yàn)室爆炸性能評(píng)定和外場(chǎng)爆炸性能評(píng)定，對(duì)新型微乳化柴油和普通軍用柴油進(jìn)行了綜合性能對(duì)比分析，認(rèn)為新型微乳化柴油能有效解決柴油油箱在受到外力打擊或意外撞擊時(shí)發(fā)生燃燒、爆炸所帶來(lái)的安全問(wèn)題[5]。

開(kāi)展液體燃料拋撒和云霧爆炸實(shí)驗(yàn)有助于掌握燃料的拋撒霧化效果及其云霧爆炸威力，它常用于燃料空氣炸藥（FAE）裝置結(jié)構(gòu)參數(shù)、爆炸分散特征和爆轟特性研究。學(xué)者們分別從FAE爆炸拋撒過(guò)程力學(xué)模型及其數(shù)學(xué)描述[6-7]、FAE裝置結(jié)構(gòu)參數(shù)和拋撒方式對(duì)云霧拋撒的影響[8-9]、FAE云霧爆轟特性及其威力[10]等方面進(jìn)行了較為系統(tǒng)的研究。但是，從云霧爆炸火球表面溫度、高溫持續(xù)時(shí)間、截面積等特征參數(shù)角度來(lái)研究燃料云霧的爆轟威力，尤其運(yùn)用相關(guān)評(píng)價(jià)方法對(duì)云霧爆炸火球特征參數(shù)進(jìn)行定量計(jì)算并評(píng)估燃料云霧爆轟性能的研究鮮有報(bào)道。

本文中，開(kāi)展-10#柴油、普通微乳化柴油（防火柴油）和新型微乳化柴油（防火抑爆柴油）的拋撒和云霧爆炸實(shí)驗(yàn)，并運(yùn)用灰色關(guān)聯(lián)分析法對(duì)其爆炸火球的特征參數(shù)進(jìn)行定量計(jì)算，從爆炸火球溫度場(chǎng)角度評(píng)估新型微乳化柴油的云霧爆炸威力及其抑爆性能，最后通過(guò)激波及其高速氣流拋撒實(shí)驗(yàn)研究新型微乳化柴油的抑爆機(jī)理。

1 拋撒和云霧爆炸實(shí)驗(yàn)介紹

1.1 柴油樣品

表1給出了本實(shí)驗(yàn)所用5種柴油樣品的組分構(gòu)成及性質(zhì)。樣品1為-10#柴油；樣品2和樣品3分別為含水質(zhì)量分?jǐn)?shù)為5%和15%的普通微乳化柴油；樣品4和5為在樣品2中添加了質(zhì)量分?jǐn)?shù)為0.2%和0.4%的高分子聚合物防霧劑構(gòu)成的新型微乳化柴油。

表1 柴油樣品的組成和理化性能Table 1 Component and physical & chemical properties of diesel fuel samples

1.2 實(shí)驗(yàn)裝置和測(cè)試儀器

1.2.1 液體燃料拋撒裝置

如圖1所示，液體燃料拋撒裝置采用薄壁圓筒結(jié)構(gòu)，容積約1 L；圓筒上、下端板為PVC塑料板，側(cè)壁為聚酯塑料外殼，上、下端板的厚度和強(qiáng)度遠(yuǎn)大于側(cè)壁，以限制燃料的軸向飛散，促使形成的燃料云霧成餅狀，盡可能地?cái)U(kuò)大云霧的覆蓋面積；中心PVC管中段裝填雷管和炸藥，兩端為阻火泡沫。

圖1 液體燃料拋撒裝置示意圖Fig. 1 Schematic of liquid fuel dispersal device

1.2.2 點(diǎn)火炸藥

中心裝藥即中心管內(nèi)裝填的雷管和炸藥，本實(shí)驗(yàn)的中心裝藥為1支8號(hào)電雷管和50 g海薩爾藥柱。

二次起爆藥柱固定在拋撒裝置附近，用于在柴油爆炸拋撒成云霧后引燃柴油云霧，由1支8號(hào)電雷管和100 g海薩爾藥柱構(gòu)成。

二次起爆藥柱應(yīng)布置在燃料較富集的區(qū)域，文獻(xiàn)[11]中總結(jié)得到了二次起爆藥柱布置點(diǎn)的經(jīng)驗(yàn)公式：

式中：L為距離中心裝藥起爆爆心的二次起爆藥柱位置，(m)；Dc為云霧最大直徑，(m)；Dh為云霧最大空洞直徑，(m)；m0為燃料裝填量，(kg)。

1 L柴油樣品的的質(zhì)量約為0.8 kg，因此，由式(1)～(3)計(jì)算可知，將二次起爆藥柱設(shè)置在距爆心1.5 m處是合適的。

1.2.3 測(cè)試系統(tǒng)

（1）高速照相機(jī)：Fastcam nltima APX型，日本Photron公司，最大拍攝頻率為12 kHz，每幀最大像素?cái)?shù)為 1 504×1 128。

（2）紅外熱成像儀：Mikronscan 7200V型，美國(guó)Mikron公司，溫度響應(yīng)時(shí)間為2 μs。

1.3 實(shí)驗(yàn)場(chǎng)布置

將兩根鐵架豎立在實(shí)驗(yàn)場(chǎng)，頂端用鐵絲相連。在鐵絲中部懸掛燃料拋撒裝置至1.3 m高（由文獻(xiàn)[12]可知燃料拋撒裝置炸高0.8～1.6 m時(shí)可以保證有足夠的空間分散和爆轟，并且具有較好的爆轟狀態(tài)）。在燃料拋撒裝置后側(cè)固定木桿，木桿頂部綁定二次起爆藥柱。高速照相機(jī)、紅外熱成像儀布置在燃料拋撒裝置前方，實(shí)驗(yàn)場(chǎng)布置如圖2所示。

圖2 實(shí)驗(yàn)場(chǎng)布置示意圖Fig. 2 Schematic of the layout of experimental field

1.4 實(shí)驗(yàn)方法

首先中心裝藥起爆拋撒柴油樣品，再運(yùn)用同步控制儀控制二次起爆藥柱在中心裝藥起爆42 ms后起爆，引燃柴油云霧團(tuán)，同時(shí)用高速照相機(jī)、紅外熱成像儀記錄引爆過(guò)程并采集火球的表面溫度。高速照相機(jī)的拍攝頻率為1 000 Hz；紅外熱成像儀的采樣周期為32 ms。

2 實(shí)驗(yàn)結(jié)果與分析

2.1 拋撒與云霧爆炸過(guò)程

運(yùn)用圖像處理軟件對(duì)高速照相機(jī)拍攝記錄的拋撒過(guò)程分幅照片進(jìn)行判讀、處理，量取拋撒最前沿射流、云霧在圖片中的徑向位移，并根據(jù)圖片尺寸與實(shí)際尺寸的比例尺，獲得實(shí)際徑向擴(kuò)展距離，即云霧的徑向擴(kuò)展半徑Rc。圖3是柴油樣品爆炸拋撒的云霧徑向擴(kuò)展半徑Rc隨時(shí)間的變化曲線(xiàn)。

圖3 柴油云霧徑向擴(kuò)展半徑隨時(shí)間的變化Fig. 3 Variation of radial extension radius of diesel fuel cloud with time

從圖3中可以看出，前2, 3 ms云霧徑向擴(kuò)展半徑急劇增大，此時(shí)拋撒液體環(huán)內(nèi)部爆炸產(chǎn)物壓力遠(yuǎn)大于氣動(dòng)阻力，爆炸作用力起主導(dǎo)作用。該階段燃料的黏性等理化性能對(duì)其拋撒作用幾乎沒(méi)有影響。4～8 ms時(shí)，5種柴油的云霧徑向擴(kuò)展半徑趨于平緩，該階段可認(rèn)為是氣動(dòng)阻力起主導(dǎo)作用，此時(shí)液體環(huán)正逐步瓦解。20 ms后，4、5號(hào)柴油的云霧徑向擴(kuò)展半徑相對(duì)較小，這是因?yàn)樗鼈冚^高的黏性抑制了造成液滴破碎的R-T和K-H不穩(wěn)定性作用，它們?nèi)砸暂^大液滴形式運(yùn)動(dòng)、擴(kuò)展，氣動(dòng)阻力較大，拋撒速度較小，分散效果差。

圖4 柴油云霧爆炸過(guò)程Fig. 4 Diesel fuel cloud explosion process

由于拋撒出的柴油液滴與氣體之間有較大的相對(duì)運(yùn)動(dòng)，在液滴表面形成了逆氣流方向的剪切力，液滴被逐層剝離形成細(xì)霧滴群。二次藥柱起爆后，細(xì)霧滴群在爆炸產(chǎn)物高溫作用下迅速汽化、燃燒發(fā)生熱爆燃。圖4是5種柴油云霧在二次起爆藥柱作用下被引燃形成爆炸火球的分幅圖片。從圖4中可以看出：二次起爆藥柱起爆后1 ms（即中心裝藥起爆后第43 ms，下同）時(shí)，柴油云霧被二次起爆藥柱引燃形成火球；200 ms時(shí)，4種微乳化柴油爆炸火球邊緣的火焰已熄滅，1號(hào)柴油云霧仍在劇烈燃燒；400 ms時(shí)，3、4、5號(hào)柴油火球基本已熄滅消散，而1號(hào)柴油火球只有徑向兩端部分火焰熄滅。

2.2 爆炸火球特征參數(shù)

爆炸火球特征參數(shù)主要包括表面溫度、高溫持續(xù)時(shí)間和火球尺寸等。運(yùn)用紅外熱成像儀自帶的MikroSpec軟件對(duì)火球熱成像圖進(jìn)行分析、處理，可以得到具體的爆炸火球特征參數(shù)，如表2所示，其中：Tm為表面最高溫度，T0為火球表面最高溫度時(shí)的表面平均溫度，Δt為各溫度區(qū)間的持續(xù)時(shí)間，S為火球最大截面積；對(duì)于粗略的估算，柴油等碳?xì)淙剂闲纬傻摹昂癜l(fā)光火焰”可以視為黑體[13]，由Stefan-Boltzmann定律可知黑體輻射度J*為

式中：Stefan-Boltzmann常數(shù)σ0=5.67×10-8W/(m2·K4)；T1為黑體絕對(duì)溫度，即火球表面最高平均溫度。

表2 柴油云霧爆炸火球特征參數(shù)Table 2 Characteristic parameters of diesel fuel cloud explosion fireball

由表2可知，前4種柴油云霧爆炸火球的表面最高溫度均超過(guò)了1 273.15 K，其中3號(hào)柴油火球表面最高溫度最大，而火球表面達(dá)到最高溫度時(shí)的表面平均溫度T0與火球表面最高平均溫度T1則由1號(hào)至5號(hào)柴油遞減，均與實(shí)驗(yàn)測(cè)量得到的火球表面平均溫度Ta隨時(shí)間的變化情況（如圖5所示）較為一致，即柴油云霧爆炸火球表面平均溫度隨著柴油樣品中水、防霧劑含量的增加而減小，可見(jiàn)火球表面最高溫度Tm在起爆初期受到炸藥爆炸產(chǎn)物的影響比較明顯，它只體現(xiàn)了爆心附近區(qū)域火球的表面溫度，不能代表火球表面整體情況。同時(shí)也說(shuō)明，采用T0、T1表示火球火焰燃燒劇烈程度是合適的。

圖5 柴油云霧火球表面平均溫度隨時(shí)間的變化Fig. 5 Variation of the surface average temperature of diesel fuel cloud fireball with time

3 抑爆性能評(píng)估

灰色關(guān)聯(lián)分析法是分析系統(tǒng)中各元素之間關(guān)聯(lián)程度或相似程度的方法，無(wú)須確定評(píng)估指標(biāo)的權(quán)重，有利于增強(qiáng)評(píng)估的客觀(guān)性[14]，因此運(yùn)用灰色關(guān)聯(lián)分析法對(duì)爆炸火球特征參數(shù)進(jìn)行定量計(jì)算，可以評(píng)估新型微乳化柴油火球的爆炸威力及其抑爆性能。

3.1 確定指標(biāo)數(shù)據(jù)列

選取4種評(píng)估柴油云霧爆炸火球特性的指標(biāo)：火球表面最高溫度時(shí)的平均溫度T0、高溫（高于1 273.15K）持續(xù)時(shí)間Δt、火球最大截面積S、火球的輻射度J0。其評(píng)價(jià)指標(biāo)值如表2所示，指標(biāo)數(shù)據(jù)列為：

3.2 指標(biāo)數(shù)據(jù)規(guī)范化處理

對(duì)指標(biāo)數(shù)據(jù)列進(jìn)行規(guī)范化處理，以便把原始數(shù)據(jù)化為無(wú)量綱、同級(jí)、正向可加的數(shù)據(jù)。

由于表2中各指標(biāo)數(shù)據(jù)均為極小型數(shù)據(jù)，即評(píng)估指標(biāo)值越小，火球爆炸威力越小，關(guān)聯(lián)度越大，抑爆性能越好，因此采用下式對(duì)指標(biāo)數(shù)據(jù)進(jìn)行標(biāo)準(zhǔn)化

式中：dij為原始數(shù)據(jù)，為指標(biāo)數(shù)據(jù)列中第i個(gè)評(píng)估對(duì)象的第j項(xiàng)指標(biāo)值，μij為標(biāo)準(zhǔn)數(shù)型，Mj和mj分別為第j項(xiàng)指標(biāo)的最大值和最小值。根據(jù)式(5)得到表3。

表3 規(guī)范化后的指標(biāo)值Table 3 Normalized index values

3.3 構(gòu)造最優(yōu)參考數(shù)據(jù)列

根據(jù)以上指標(biāo)數(shù)列構(gòu)造出最優(yōu)參考數(shù)據(jù)列。最優(yōu)參考數(shù)據(jù)列由各項(xiàng)指標(biāo)在以上數(shù)據(jù)列中最優(yōu)的值所構(gòu)成，即 x0=（1，1，1，1）。

3.4 計(jì)算最優(yōu)關(guān)聯(lián)度，并排序評(píng)估

計(jì)算最優(yōu)關(guān)聯(lián)度τi，即xi與最優(yōu)參考數(shù)據(jù)列x0的關(guān)聯(lián)度。關(guān)聯(lián)分析中被比較數(shù)列常記為xi，一般表示為：

對(duì)于一個(gè)參考數(shù)據(jù)列x0，比較數(shù)列為xi，可用下述關(guān)系表示比較點(diǎn)與參考點(diǎn)在各點(diǎn)的差：

式中：δik(k)是第k個(gè)時(shí)刻xi與x0的相對(duì)差值，稱(chēng)為xi對(duì)x0在k時(shí)刻的關(guān)聯(lián)系數(shù)。在實(shí)際使用時(shí)，應(yīng)根據(jù)序列間的關(guān)聯(lián)程度選擇分辨系數(shù)，一般取δ=0.5。因此，由式(7)計(jì)算可得：δ11(1)=δ12(2)=δ13(3)=δ14(4)=0.333 3。

關(guān)聯(lián)系數(shù)只表示各點(diǎn)數(shù)據(jù)間的關(guān)聯(lián)程度，由于關(guān)聯(lián)系數(shù)的數(shù)很多，信息過(guò)于分散，不便于比較，為此有必要將各點(diǎn)的關(guān)聯(lián)系數(shù)集中為一個(gè)值，求平均值便是作為這種信息集中處理的一種方法。最優(yōu)關(guān)聯(lián)度的一般表達(dá)式為

根據(jù)式（8）計(jì)算得到：τ1=0.333 3，τ2=0.567 3，τ3=0.725 0，τ4=0.730 2，τ5=1

由此可得，在被評(píng)估的5種柴油中，1號(hào)柴油的爆炸威力最大，3、4號(hào)柴油的爆炸威力比較接近，而5號(hào)柴油的爆炸威力最小。可見(jiàn)，新型微乳化柴油具有良好的抑爆性能，隨著防霧劑含量的增加，其抑爆性能也越好，而且從3、4號(hào)柴油評(píng)估結(jié)果可以確定，柴油的含水量在15%以下時(shí)，多增加10%的水與加入質(zhì)量分?jǐn)?shù)為0.002的防霧劑的抑爆效果相當(dāng)。

4 抑爆機(jī)理分析

運(yùn)用液體燃料拋撒和成像系統(tǒng)可以記錄柴油柱在激波及其高速氣流作用下的拋撒霧化情況。液體燃料拋撒和成像系統(tǒng)由水平激波管、壓力傳感器、信號(hào)調(diào)理器、時(shí)序控制器、光源控制器、高壓閃光電源、凹面反射鏡和單反照相機(jī)（BULB模式）等組成，如圖6所示。

圖6 液體燃料拋撒和成像系統(tǒng)俯視示意圖Fig. 6 Vertical schematic of liquid fuel dispersal and imaging system

水平激波管管口外部15 mm處自上而下布置有油杯和上、下兩根出油管，上、下出油管之間留有3 mm的空隙以形成油柱。出油管為不銹鋼管，上管的內(nèi)徑為3 mm，下管的內(nèi)徑為5 mm，如圖7所示。

圖8是油柱在1.75馬赫激波作用下的霧化情況，圖9是運(yùn)用Spraytec型實(shí)時(shí)噴霧激光粒度分析儀測(cè)定得到的油樣拋撒霧滴特征平均粒徑，其中DV(10)、DV(50)、DV(90)為累積體積分?jǐn)?shù)為10%、50%、90%時(shí)對(duì)應(yīng)的霧滴粒徑，即小于此粒徑的霧滴體積分別占全部霧滴總體積的10%、50%、90%。D[3,2]為表面積動(dòng)量平均徑，即索態(tài)爾（Stauter）平均直徑。D[4,3]為體積或質(zhì)量動(dòng)量平均徑。

圖7 水平激波管示意圖Fig. 7 Schematic of a horizontal shock tube

圖8中，3號(hào)云霧團(tuán)的顏色比1、2號(hào)稍淺，說(shuō)明霧滴平均粒徑較大，密度較小。4、5號(hào)柴油霧化效果差，呈液絲、液塊狀向前運(yùn)動(dòng)，只有少部分的細(xì)小霧滴從液絲上脫落出來(lái)，而由圖9可見(jiàn)4、5號(hào)柴油的5種特征平均粒徑均顯著高于前3種柴油霧滴粒徑，如4、5號(hào)柴油霧滴索態(tài)爾平均直徑D[3,2]比1號(hào)柴油霧滴D[3,2]大77.7及110.6 μm?？傮w上，圖8反映的拋撒霧化現(xiàn)象與圖9中5種柴油霧滴平均粒徑大小排列是吻合的。

圖8 激波及高速氣流作用下柴油柱的霧化過(guò)程Fig. 8 Atomization process of diesel fuel column by shock wave and high speed airflow

由此可見(jiàn)，液體燃料的運(yùn)動(dòng)黏度和黏性對(duì)燃料的拋撒霧化效果有著至關(guān)重要的影響。尤其是，在燃料中添加高分子聚合物防霧劑后增強(qiáng)了其“霧化抑制效應(yīng)”，防霧劑的高分子鏈變形和松弛造成燃料的黏彈性增大[15-16]，如圖10所示，其中ω是樣品小幅振動(dòng)的角頻率，3號(hào)柴油的黏性模量始終稍大于彈性模量，表明黏性起主要作用，而5號(hào)柴油的彈性模量值顯著大于自身的黏性模量以及3號(hào)柴油的彈性模量，表明它具有較高的彈性，因此新型微乳化柴油在高速拋撒過(guò)程中分子鏈的拉伸將產(chǎn)生抵制外部變形的拉伸應(yīng)力[17-19]，使液滴在高速氣流剪切作用下也不易破碎，液滴分散效果差，抑爆性能好。

圖9 霧滴特征平均粒徑Fig. 9 Characteristic mean particle size of droplet

圖10 柴油樣品的小振幅振蕩流變曲線(xiàn)Fig. 10 Small amplitude oscillatory rheological curves of diesel fuel samples

5 結(jié) 論

（1）拋撒和云霧爆炸實(shí)驗(yàn)中，中心藥柱起爆拋撒20 ms后新型微乳化柴油的拋撒云霧徑向擴(kuò)展半徑顯著小于其他柴油，如在24 ms時(shí)，4、5號(hào)新型微乳化柴油拋撒云霧徑向擴(kuò)展半徑比1號(hào)-10#柴油分別小15.43和15.95 cm。新型微乳化柴油云霧爆炸火球特征參數(shù)也較小，4、5號(hào)柴油的火球的表面最高平均溫度比1號(hào)柴油分別低296.90和336.90 K；火球表面高溫（高于1 273.15 K）持續(xù)時(shí)間方面，4、5號(hào)柴油為140、0 ms，分別比1號(hào)柴油少94和234 ms；4、5號(hào)柴油的火球最大截面積也分別只有1號(hào)柴油的60.10%、53.53%。

（2）運(yùn)用灰色關(guān)聯(lián)分析法對(duì)爆炸火球的表面最高溫度時(shí)的平均溫度、高溫持續(xù)時(shí)間、最大截面積、輻射度等特征參數(shù)進(jìn)行定量計(jì)算，評(píng)估了5種柴油火球的爆炸威力。結(jié)果表明，隨著柴油中的含水量和防霧劑含量的增加，微乳化柴油的爆炸威力越小，其抑爆性能越好。另外，3、4號(hào)柴油的抑爆性能接近，說(shuō)明柴油的含水量在15%以下時(shí)，多增加10%的水與添加0.2%防霧劑的抑爆效果相當(dāng)。

（3）運(yùn)用激波及其高速氣流拋撒霧化柴油樣品，再應(yīng)用激光粒度分析儀測(cè)定拋撒霧化形成的霧滴平均粒徑，4、5號(hào)柴油的霧滴特征平均粒徑顯著大于其他柴油，如4、5號(hào)柴油霧滴索態(tài)爾平均直徑D[3,2]比1號(hào)柴油霧滴D[3,2]大77.7和110.6 μm。又根據(jù)其流變性能確定，新型微乳化柴油中添加防霧劑使其液滴黏彈性增大，在高速氣流剪切作用不易破碎、霧化，液滴分散效果差，抑爆性能好。