基于響應(yīng)面法的木粉塵最大爆炸壓力試驗(yàn)研究
——以桑木粉塵為例

萬(wàn)杭煒，趙江平

(西安建筑科技大學(xué) 材料與礦資學(xué)院，陜西 西安 710055)

0　引言

粉塵爆炸是指可燃性粉塵與助燃性氣體在密閉或半密閉空間內(nèi)觸及點(diǎn)火源時(shí)發(fā)生的伴有光、熱、聲的快速化學(xué)反應(yīng)過(guò)程[1]。粉塵爆炸會(huì)造成嚴(yán)重的財(cái)產(chǎn)損失或人員傷亡[2]。因此對(duì)粉塵爆炸研究已經(jīng)刻不容緩。最大爆炸壓力是指爆炸過(guò)程中產(chǎn)生的最大壓力，它是衡量爆炸強(qiáng)度的指標(biāo)之一[1]，它是作為爆炸防護(hù)措施的參考依據(jù)[3]。木粉塵是纖維類粉塵中發(fā)生事故最多的一種可燃性粉塵[4]。

針對(duì)木粉塵，國(guó)內(nèi)研究重點(diǎn)主要在職業(yè)健康方面。國(guó)內(nèi)外學(xué)者對(duì)木粉塵爆炸特性也進(jìn)行了相關(guān)研究。例如，陳國(guó)華等[5]研究了紅木粉的最小點(diǎn)火能隨粉塵濃度和粉塵粒徑的變化關(guān)系，表明最小點(diǎn)火能隨著粉塵濃度增加先減小后增大，隨著粒徑的增大而升高；任瑞娥等[6]研究了木粉的最低著火溫度與粉塵粒徑和分散壓力的關(guān)系，結(jié)果表明，最低著火溫度隨著粉塵粒徑的增大而升高，隨著分散壓力的增大而升高；Amyotte等[7]研究了木纖維、聚乙烯樣品和兩者混合時(shí)的爆炸特性，研究表明，粉塵粒徑越小，發(fā)生爆炸可能性越大，產(chǎn)生的爆炸后果越嚴(yán)重。當(dāng)木纖維和聚乙烯混合時(shí)，產(chǎn)生的爆炸威力更大；Clara[8]等對(duì)不同種類的煤粉和生物質(zhì)粉塵在1 m3標(biāo)準(zhǔn)爆炸測(cè)試裝置進(jìn)行了爆炸壓力、爆炸下限等爆炸特性參數(shù)研究。

可知，對(duì)木粉塵爆炸特性研究還比較匱乏。并且學(xué)者一般只是對(duì)爆炸特性參數(shù)的影響因素進(jìn)行研究，而沒(méi)有對(duì)其影響因素進(jìn)行排序和研究發(fā)生粉塵最大爆炸壓力的最適宜條件。基于此，筆者以桑木粉塵為試驗(yàn)對(duì)象，對(duì)其最大爆炸壓力進(jìn)行研究，以期在預(yù)防木粉塵爆炸事故方面為相關(guān)企業(yè)提供一定的參考。

1　試驗(yàn)材料與方法

1.1　試驗(yàn)材料及其表征

選取某家具廠加工木材時(shí)所產(chǎn)生的桑木鋸末為試驗(yàn)原料。木材主要由C，H，O和1% N元素構(gòu)成，化學(xué)主要組成有纖維素、半纖維素和木素，它是一種天然的高分子化合物[9]。木粉的工業(yè)分析如下表1所示。

表1　木粉工業(yè)分析Table 1　Proximate analysis of the wood dust

1.2　試驗(yàn)條件

單因素試驗(yàn)在室內(nèi)溫度為10 ℃，空氣濕度為83%的條件下進(jìn)行；響應(yīng)面試驗(yàn)在室內(nèi)溫度為10 ℃，空氣濕度為71%的條件下進(jìn)行。

1.3　試驗(yàn)預(yù)處理

1.3.1篩分

利用振動(dòng)篩把粉塵過(guò)80，100，120，140，160，200，250，350目篩進(jìn)行篩分，得到其質(zhì)量分?jǐn)?shù)分別為3.52%，30.63%，18.66%，15.49%，15.14%，7.04%，7.39%，2.11%。根據(jù)其質(zhì)量分?jǐn)?shù)可畫出在一定范圍內(nèi)的木粉粒徑分布圖，如圖1所示。

圖1　木粉粒徑分布Fig.1　Diameter distribution of the wood dust

1.3.2干燥

把篩分好的試驗(yàn)原料統(tǒng)一放到烘箱內(nèi)在50 ℃條件下干燥2 h。

1.4　試驗(yàn)設(shè)備

儀器設(shè)備包括東北大學(xué)安全工程研究中心研制的1.2 L Hartmann管、空氣壓縮機(jī)、電子天平、吸塵器、標(biāo)準(zhǔn)篩、振動(dòng)篩，101-2AB型電熱鼓風(fēng)干燥箱，Gems捷邁壓力傳感器(量程為0～5 MPa， 響應(yīng)時(shí)間為微秒級(jí) ，采集頻率為一個(gè)數(shù)據(jù)點(diǎn)1 ms，通過(guò)采集軟件自行設(shè)置)。Hartmann管式粉塵爆炸儀如圖2所示。

圖2　1.2 L Hartmann管式粉塵爆炸裝置Fig.2　Dust explosion device of 1. 2 L Hartmann tube

1.5　單因素試驗(yàn)

1.5.1粒徑

分別稱取100，160，200，250，350目5種不同粒徑粉塵0.9 g，點(diǎn)火延遲時(shí)間設(shè)置為10 s，噴粉壓力為0.6 MPa，比較不同粒徑對(duì)最大爆炸壓力影響。

1.5.2質(zhì)量濃度

分別稱取粒徑為100目的粉塵0.3，0.6，0.9，1.2，1.5 g，其相應(yīng)質(zhì)量濃度分別為250，500，750，1000，1 250 g/m3。點(diǎn)火延遲時(shí)間設(shè)置為10s，噴粉壓力為0.6 MPa，比較不同質(zhì)量濃度對(duì)最大爆炸壓力影響。

1.5.3點(diǎn)火延遲時(shí)間

稱取粒徑大小為120目的粉塵0.9 g，點(diǎn)火延遲時(shí)間分別設(shè)置為6, 8,10,12,14 s，噴粉壓力為0.6 MPa，比較不同點(diǎn)火延遲時(shí)間對(duì)最大爆炸壓力影響。

1.6　響應(yīng)面試驗(yàn)[10-11]

在單因素試驗(yàn)基礎(chǔ)上運(yùn)用Design-Expert 8.0.6軟件中的Box-Behnken Design對(duì)試驗(yàn)進(jìn)行設(shè)計(jì)。

文獻(xiàn)[12]指出，粉塵爆炸跟粉體性質(zhì)、粉塵云性質(zhì)和外界條件有關(guān)。粉體性質(zhì)有粉塵粒度、形狀、比表面積、燃燒熱等。粉塵粒度大小可以通過(guò)篩網(wǎng)篩分控制。粒徑越小，比表面積越大，而其形狀、燃燒熱等由于是本身固有屬性，故選取粉塵粒徑作為變量因素；粉塵云性質(zhì)有粉塵濃度、濕度、含氧量、湍流度等。粉塵云濃度可以通過(guò)粉塵質(zhì)量的多少來(lái)改變。湍流度會(huì)直接或間接導(dǎo)致粉塵云濃度變化。濕度通過(guò)烘箱干燥處理，一方面是為了控制單一變量，另一方面是由于粉塵水分的存在，會(huì)降低粉塵爆炸性，這對(duì)防爆設(shè)備的強(qiáng)度設(shè)計(jì)上會(huì)產(chǎn)生一定的干擾。含氧量由于在實(shí)際工業(yè)過(guò)程中較充足，因此不作考量范圍。所以綜合考慮，選取經(jīng)過(guò)干燥處理的不同質(zhì)量粉塵作為變量因素；外界條件有初始溫度、初始?jí)毫ΑⅫc(diǎn)火源溫度，容器尺寸及形狀等。點(diǎn)火源溫度通過(guò)點(diǎn)火延遲時(shí)間來(lái)改變。初始溫度、壓力不作處理，是因?yàn)榇蠖鄶?shù)粉塵爆炸事故是在常溫常壓下由于外部點(diǎn)火源引起的。粉塵最大爆炸壓力在忽略容器熱損失時(shí)，與容器尺寸和形狀無(wú)關(guān)。因此選擇點(diǎn)火延遲時(shí)間作為變量因素。

綜上，為了研究木粉的最大爆炸壓力并結(jié)合實(shí)際情況，選取粉塵質(zhì)量濃度A，粒徑B，點(diǎn)火延遲時(shí)間C三因素為響應(yīng)變量，最大爆炸壓力Pm(R1)為響應(yīng)值，設(shè)計(jì)三因素三水平響應(yīng)面試驗(yàn)，各設(shè)計(jì)因素水平值與編碼值如表2所示。

1.7　數(shù)據(jù)處理

試驗(yàn)數(shù)據(jù)由Arduino1.8.2軟件和Serial Port Utility軟件共同獲得。單因素試驗(yàn)所得數(shù)據(jù)用MATLAB R2015a軟件畫圖，響應(yīng)面試驗(yàn)所得數(shù)據(jù)用Design-Expert 8.0.6軟件分析。

表2　Box-Behnken Design因素水平及編碼Table 2　Levels of factors and codes in Box-Behnken

2　試驗(yàn)結(jié)果與分析

2.1　單因素試驗(yàn)

2.1.1粒徑對(duì)最大爆炸壓力影響

圖3所示的是在粉塵質(zhì)量濃度為750 g/m3下測(cè)得的最大爆炸壓力隨粒徑變化的曲線。從圖中可以看出，隨著目數(shù)不斷增大，最大爆炸壓力逐漸增大。當(dāng)目數(shù)大于200時(shí)，增大趨勢(shì)已不明顯。這是因?yàn)殡S著目數(shù)增大，粉塵粒徑越小，粉塵比表面積增大，其與氧接觸面積增大，加速了粉塵與氧氣的反應(yīng)，燃燒速率增大，產(chǎn)生的熱量越多，爆炸壓力越大，當(dāng)粉塵比表面積增大到一定程度時(shí)，由于氧氣不足，導(dǎo)致最大爆炸壓力增大不顯著[13]。因此，選擇中心點(diǎn)粒徑為200目。

圖3　最大爆炸壓力隨粒徑變化情況Fig.3　The maximum explosion pressure varies with particle size

2.1.2濃度對(duì)最大爆炸壓力影響

圖4所示的是在粉塵粒徑為100目下測(cè)得的最大爆炸壓力隨質(zhì)量濃度變化的曲線。

圖4　最大爆炸壓力隨質(zhì)量濃度變化情況Fig.4　The maximum explosion pressure varies with mass concentration

從圖可知，隨著粉塵濃度的增加，最大爆炸壓力先增大后減小。初始粉塵質(zhì)量濃度增加時(shí)，最大爆炸壓力增大，當(dāng)質(zhì)量濃度超過(guò)750 g/m3,最大爆炸壓力減小。這是因?yàn)殡S著粉塵濃度增大，分子間有效碰撞次數(shù)越多，化學(xué)反應(yīng)速率加快，產(chǎn)生的能量越多，當(dāng)質(zhì)量濃度達(dá)到一定濃度時(shí)，裝置內(nèi)所含氧濃度不足占主導(dǎo)地位，未發(fā)生反應(yīng)的粉塵越來(lái)越多，并且未反應(yīng)粉塵會(huì)吸收參加反應(yīng)粉塵所產(chǎn)生的燃燒熱，使體系最大爆炸壓力降低[14]。因此，選擇中心點(diǎn)質(zhì)量濃度為750 g/m3。

2.1.3點(diǎn)火延遲時(shí)間對(duì)最大爆炸壓力影響

電熱絲加熱作為噴粉引爆的點(diǎn)火源。點(diǎn)火延遲時(shí)間指從電熱絲加熱到開始噴粉的這段時(shí)間。

從圖5可看出，隨著點(diǎn)火延遲時(shí)間的增大，最大爆炸壓力逐漸上升。根據(jù)焦耳定律公式：Q=I2Rt,當(dāng)I與R不變時(shí)，Q與t成正比關(guān)系，故電阻絲能量隨著時(shí)間的增加而增大。增大點(diǎn)火能量能提高物質(zhì)燃燒速率，從而反應(yīng)燃燒熱增加，使最大爆炸壓力增大[15]。當(dāng)點(diǎn)火延遲時(shí)間超過(guò)10 s時(shí)，最大爆炸壓力增加減慢。因此，選擇中心點(diǎn)點(diǎn)火延遲時(shí)間為10 s。

圖5　最大爆炸壓力隨點(diǎn)火延遲時(shí)間變化情況Fig.5　The maximum explosion pressure varies with the ignition delay time

2.2　響應(yīng)面試驗(yàn)

Box-Behnken Design方案和試驗(yàn)結(jié)果如表3。

表3　響應(yīng)面試驗(yàn)設(shè)計(jì)方案與結(jié)果Table 3　Experimental design and results ofresponse surface methodology

續(xù)表3

利用Design-Expert軟件分析，得到推薦模型為二階模型。對(duì)所得模型進(jìn)行方差分析和顯著性檢驗(yàn)。如表4所示。

表4　方差分析Table 4　Analysis of variance table

注：SS-平方和，DF-自由度，MS-均方，P(Pr>F)-無(wú)顯著影響概率；差異顯著(P<0.05)，差異高度顯著(P<0.01)，差異極顯著(P<0.001)。

擬合回歸方程為R1=0.44+0.03A+0.028B+0.064C-0.008 375AB+0.002 375AC+0.001 125BC-0.033C2-0.007 091B2-0.02C2。

2.2.1顯著性檢驗(yàn)

從表4知，模型P(Pr>F)<0.000 1，表明模型極顯著，可以利用該模型進(jìn)行優(yōu)化。失擬項(xiàng)P(Pr>F)>0.05表明不顯著，擬合性較好。質(zhì)量濃度A(P=0.000 3)，粒徑B(P=0.000 5)，點(diǎn)火延遲時(shí)間C(P<0.000 1)，A2(P=0.000 1)都達(dá)到極顯著水平；C2(P=0.004 5)高度顯著；AB，AC，BC都不顯著。綜上表明，試驗(yàn)因子對(duì)響應(yīng)值不是線性關(guān)系。

用Design-Expert軟件對(duì)回歸方程進(jìn)行誤差統(tǒng)計(jì)分析，如表5所示。

表5　誤差統(tǒng)計(jì)分析表Table 5　Error statistics table

從文獻(xiàn)[16]知，表5中數(shù)據(jù)達(dá)到檢驗(yàn)值，說(shuō)明回歸方程模型可行性好。

結(jié)合F值大小和回歸方程系數(shù)，可得因子對(duì)響應(yīng)值影響程度順序?yàn)椋狐c(diǎn)火延遲時(shí)間>粉塵質(zhì)量濃度>粉塵粒徑。

2.2.2響應(yīng)面分析[17]

根據(jù)回歸方程模型作出各因子交互作用的響應(yīng)面圖，如圖6～8所示。

圖6　點(diǎn)火延遲時(shí)間和粒徑對(duì)最大爆炸壓力的交互影響Fig.6　Effects of ignition delay time and particle size on maximum explosion pressure

圖7　質(zhì)量濃度和粒徑對(duì)最大爆炸壓力的交互影響Fig.7　Effects of mass concentration and particle size on maximum explosion pressure

圖8　點(diǎn)火延遲時(shí)間和質(zhì)量濃度對(duì)最大爆炸壓力的交互影響Fig.8　Effects of ignition delay time and mass concentration on maximum explosion pressure

從圖6-8中可以直觀地看出，相應(yīng)的響應(yīng)面均為凸型曲面，說(shuō)明最大爆炸壓力存在極值，且取得極值即最大值的點(diǎn)的水平投影在所選因素水平范圍內(nèi)。

圖6表明，點(diǎn)火延遲時(shí)間改變所引起最大爆炸壓力的變化較粒徑坡度大，這表明點(diǎn)火延遲時(shí)間的影響大于粒徑，與2.2.1所得結(jié)果相符。同時(shí)整體坡度變化不大，表明粒徑與點(diǎn)火延遲時(shí)間的交互作用不顯著，即BC不顯著，這在表4可體現(xiàn)。

圖7表明，質(zhì)量濃度所引起的最大爆炸壓力的改變較大，說(shuō)明質(zhì)量濃度的影響大于粒徑，與2.2.1所得結(jié)果相符。

圖8表明，點(diǎn)火延遲時(shí)間改變所引起最大爆炸壓力的變化較質(zhì)量濃度坡度大，這表明點(diǎn)火延遲時(shí)間的影響大于質(zhì)量濃度，與2.2.1所得結(jié)果相符。

2.2.3最危險(xiǎn)爆炸強(qiáng)度

最危險(xiǎn)爆炸強(qiáng)度是指在不同水平下測(cè)得的最大爆炸壓力的最大值[1]。Design-Expert軟件為響應(yīng)面分析提供了最佳優(yōu)化方案。由軟件分析得到最佳試驗(yàn)方案為：質(zhì)量濃度840.24g.m-3，粒徑260目，點(diǎn)火延遲時(shí)間12s，最危險(xiǎn)強(qiáng)度預(yù)測(cè)值為0.511 775MPa。

在此試驗(yàn)條件下進(jìn)行驗(yàn)證，得到最大爆炸壓力為0.508MPa，與預(yù)測(cè)值非常接近，說(shuō)明Design-Expert軟件提供的最佳優(yōu)化方案的可行性和有效性。

3　結(jié)論

1)通過(guò)單因素試驗(yàn)得出：桑木粉塵最大爆炸壓力隨著粒徑減小而逐漸增大；隨著粉塵濃度的增大，最大爆炸壓力先增大后減小，且粉塵濃度在750 g/m3取得最大值；隨著點(diǎn)火延遲時(shí)間的增大，最大爆炸壓力逐漸增大。

2)通過(guò)Design-Expert軟件對(duì)Box-Behnken Design試驗(yàn)方案結(jié)果進(jìn)行回歸模型分析并結(jié)合響應(yīng)曲面，得到在粒徑中心點(diǎn)為200目，質(zhì)量濃度中心點(diǎn)為750 g/m3，點(diǎn)火延遲時(shí)間中心點(diǎn)為10s下影響最大爆炸壓力大小順序?yàn)椋狐c(diǎn)火延遲時(shí)間>粉塵質(zhì)量濃度>粒徑。

3)Design-Expert軟件分析取得最危險(xiǎn)爆炸強(qiáng)度試驗(yàn)條件為：質(zhì)量濃度840.24 g/m3，粒徑260目，點(diǎn)火延遲時(shí)間12s，并且最大爆炸壓力預(yù)測(cè)值為0.511 775 MPa，經(jīng)驗(yàn)證可行，這對(duì)預(yù)防粉塵爆炸的防爆設(shè)備強(qiáng)度本質(zhì)安全設(shè)計(jì)提供了一定的參考依據(jù)。

[1]畢明樹，楊國(guó)剛. 氣體和粉塵爆炸防治工程學(xué)[M]. 北京:化學(xué)工業(yè)出版社，2012： 3-90.

[2]多英全，劉垚楠，胡馨升. 2009～2013年我國(guó)粉塵爆炸事故統(tǒng)計(jì)分析研究[J]. 中國(guó)安全生產(chǎn)科學(xué)技術(shù)，2015(2):186-190.

DUO Yingquan，LIU Yaonan，HU Xingsheng.Statistical analysis on dust explosion accidents occurring in China during 2009-2013[J].Journal of Safety Science and Technology，2015(2):186-190.

[3]覃欣欣，李潤(rùn)之，司榮軍. 酚醛樹脂粉塵爆炸危險(xiǎn)性實(shí)驗(yàn)研究[J].工業(yè)安全與環(huán)保，2015(12):12-14.

QIN Xinxin，LI Runzhi，SI Rongjun. Experimental research on hazard of phenolic resin dust explosion [J].Industrial Safety and Environmental Protection，2015(12):12-14.

[4]胡雙啟. 燃燒與爆炸[M]. 北京: 北京理工大學(xué)出版社，2015: 74-79.

[5]陳國(guó)華，江湖一佳，王新華. 紅木粉爆炸特性實(shí)驗(yàn)研究[J].消防科學(xué)與技術(shù)，2016(9):1200-1203.

CHEN Guohua，JIANG Huyijia，WANG Xinhua.Study on characteristics of mahogany wood dust explosion[J].Fire Science and Technology, 2016(9):1200-1203.

[6]任瑞娥，譚迎新. 木粉最低著火溫度的實(shí)驗(yàn)研究[J]. 中國(guó)粉體技術(shù)，2014(5):45-47.

REN Ruie，TAN Yingxin.Experimental study on minimum ignition temperature of wood dust cloud[J].China Powder Science and Techology，2014(5):45-47.

[7]P.R. Amyotte，C.T. Cloney， F.I. Khan，R.C. Ripley，Dust explosion risk moderation for flocculent dusts，J. Loss Prev. Proc. Ind. 25 (2012) 862-869.

[8]Paul R. Amyotte，Chris T.Cloney，F(xiàn)raisal I. Khan et al. Dust explosion risk moderation for flocculent dusts[J].Journal of Loss Prevention in the Process Industries，2012，25(5)：862-869.

[9]葛明裕,等. 木材加工化學(xué)[M]. 哈爾濱：東北林業(yè)大學(xué)出版社, 1995:14-16.

[10]Box G，Wilson K B.On the experimental attainment of optimum conditions[J].Journal of the R oyal Statistical Society Series B:Statistical Methodology， 1951(1): 1-45.

[11]Mason R L ，Gunst R F， Hess J L.Statistical design and analysis ofexperiments with applications to engineering and science[M].The United States of America: John Wiley and Sons Publication，2003.

[12]鄔長(zhǎng)城，薛偉，賈愛忠，等. 燃燒爆炸理論基礎(chǔ)與應(yīng)用. 北京：化學(xué)工業(yè)出版社[M]. 2016:74-78.

[13]Dufaud O，Traoré M，Perrin L，et al. Experimental investigation and modelling of aluminum dusts explosions in the 20 L sphere[J].Journal of Loss Prevention in the Process Industries，2010， 23(2): 226-236.

[14]譚靜怡，袁旌杰，杜兵，等. 含濕煤粉爆炸行為的實(shí)驗(yàn)研究[J]. 消防科學(xué)與技術(shù)，2014(3):247-250.

TAN Jingyi，YUAN Jingjie，DU Bing，et al. Experimental investigation on explosion behaviors of pre-humidification coal dust[J].Fire Science and Technology，2014(3):247-250.

[15]蒯念生，黃衛(wèi)星，袁旌杰，等. 點(diǎn)火能量對(duì)粉塵爆炸行為的影響[J]. 爆炸與沖擊，2012(4):432-438.

KUAI Niansheng，HUANG Weixing，YUAN Jingjie， et al. Influence of ignition energy on dust explosion behavior[J].Explosion and Shock Waves，2012(4):432-438.

[16]李莉，張賽，何強(qiáng)，等. 響應(yīng)面法在試驗(yàn)設(shè)計(jì)與優(yōu)化中的應(yīng)用[J]. 實(shí)驗(yàn)室研究與探索，2015(8):41-45.

LI Li，ZHANG Sai，HE Qiang，et al.Application of response surface methodology in experiment design and optimization[j].Research and Exploration in Laboratory，2015(8):41-45.

[17]姚笛，馬萍，王穎，等. 響應(yīng)面法優(yōu)化玉米芯中木聚糖的提取工藝[J]. 食品科學(xué)，2011(8):111-115.

YAO Di，MA Ping，WANG Ying，et al. Optimization of extraction process for xylan from corncob by response surface methodology[J].Food Science，2011(8):111-115.