巫 亮
(1.中煤科工集團(tuán)重慶研究院有限公司,重慶 400037;2.瓦斯災(zāi)害監(jiān)控與應(yīng)急技術(shù)國家重點(diǎn)實(shí)驗(yàn)室,重慶 400039)
近年來,大氣顆粒物的污染愈發(fā)嚴(yán)重,對于大氣顆粒物濃度的準(zhǔn)確檢測成為了迫在眉睫的任務(wù)[1-5]。β射線吸收法是我國環(huán)保部門推薦的自動(dòng)PM10、PM2.5粉塵濃度檢測的兩種方法之一[6]。其檢測原理是利用濾膜上的顆粒物對β射線的衰減來檢測顆粒物的質(zhì)量濃度,具有測量準(zhǔn)確性高、維護(hù)量小和成本低等優(yōu)點(diǎn)[7-10]。
目前,國外的賽默飛世特、國內(nèi)的河北天河、武漢天虹等廠商都具有相對成熟的β射線自動(dòng)監(jiān)測儀器。一些科研院校也對其工作原理、結(jié)構(gòu)形式及控制方法進(jìn)行了較為深入的研究[11-13]。但是,國內(nèi)外對于β射線吸收法的粉塵濃度計(jì)算公式的研究卻少了涉及。目前采用的計(jì)算方法為確定濾膜的質(zhì)量吸收系數(shù)μm之后,分別測量濾膜抽塵前后β射線穿過濾膜時(shí)的強(qiáng)度,根據(jù)β射線吸收定律計(jì)算出濾膜的重量變化,從而計(jì)算出大氣中的粉塵濃度。但是,根據(jù)實(shí)驗(yàn)分析,濾膜在抽塵過程中受負(fù)壓影響,其質(zhì)量吸收系數(shù)μm會(huì)產(chǎn)生變化。而現(xiàn)有的粉塵濃度計(jì)算公式并沒有引入這一影響因素,從而會(huì)影響測試結(jié)果的準(zhǔn)確度。
綜上所述,為了消除β射線吸收法粉塵濃度計(jì)算公式中濾膜質(zhì)量吸收系數(shù)μm變化引起的誤差,提出了質(zhì)量吸收系數(shù)的影響因子ε的計(jì)算方法。并將其引入到現(xiàn)有的粉塵濃度計(jì)算公式中,推導(dǎo)出了考慮濾膜吸收系數(shù)μm變化的基于β射線吸收法的粉塵濃度的計(jì)算新公式。經(jīng)過實(shí)驗(yàn)驗(yàn)證,該公式有效的提高了β射線粉塵濃度檢測設(shè)備的檢測精度,特別是在低濃度情況下的檢測精度。
β射線是從同位元素C14放射性衰變中釋放出來的高速電子流,其電力本領(lǐng)較小、穿透能力較強(qiáng)。當(dāng)β射線穿過物質(zhì)時(shí),其衰減服從指數(shù)分布規(guī)律,其衰減公式如下[14]:
I=I0e-μmρL
(1)
式中:I0與I分別為β射線穿過介質(zhì)前后的射線強(qiáng)度,單位為Bq;L為β射線投射物質(zhì)的厚度,單位為cm;ρ為β射線穿過物質(zhì)的密度,單位為g/cm3;μm為被投射物質(zhì)對β射線的吸收系數(shù),單位為cm2/g。
由式(1)可知,當(dāng)β射線穿透物質(zhì)的質(zhì)量吸收系數(shù)μm只與其密度及厚度有關(guān)。令dm=ρL=m/πr2,dm為被投射物質(zhì)的質(zhì)量面積,得到式(2):
(2)
式中:r為被檢測物質(zhì)的檢測半徑,單位為cm,m為被檢測面積內(nèi)的物質(zhì)質(zhì)量,單位為g。
反解式(2)得到被檢測質(zhì)量m的計(jì)算公式
(3)
在運(yùn)用β射線吸收法檢測吸收粉塵的質(zhì)量時(shí),首先檢測濾紙抽塵前的β射線吸收強(qiáng)度I1,計(jì)算出此時(shí)的濾膜重量m1,式(4)。然后對濾帶進(jìn)行抽塵,抽塵結(jié)束之后重新檢測其β射線吸收強(qiáng)度I2,并計(jì)算出此時(shí)的濾膜重量m2,式(5)。
(4)
(5)
式(5)和式(4)相減得到濾膜抽塵前后的增重重量:
(6)
最后通過測量抽氣的體積V,計(jì)算出大氣中顆粒物的濃度:
(7)
β射線測量大氣顆粒物濃度的準(zhǔn)確性的關(guān)鍵在于質(zhì)量吸收系數(shù)的確定。決定質(zhì)量吸收系數(shù)大小的因素主要有以下4個(gè)因素[12]:β射線穿透的物質(zhì)的物理組成結(jié)構(gòu)、射線源能量大小、射線源到探測器之間的距離以及濾膜的質(zhì)量密度。由式(7)可知,在現(xiàn)有的粉塵濃度計(jì)算公式中,計(jì)算濾膜抽塵前后的濾膜重量時(shí),質(zhì)量吸收系數(shù)μm為一定值。而事實(shí)上,在濾膜抽塵時(shí),濾膜承受1 500 Pa左右的負(fù)壓,此時(shí)濾膜的密度ρ以及厚度L會(huì)產(chǎn)生變化,由式(1)可知,在放射源強(qiáng)度以及測試距離不變的情況下,其質(zhì)量吸收系數(shù)也會(huì)產(chǎn)生變化。為證明質(zhì)量吸收系數(shù)在負(fù)壓作用下的變化設(shè)計(jì)了如下試驗(yàn):使用β射線自動(dòng)監(jiān)測設(shè)備中使用較為廣泛的玻璃纖維卷式濾膜,分別測試濾膜在抽塵前后的β射線強(qiáng)度,抽塵過程中所抽氣體為零氣,抽塵時(shí)間為30 s,測試濾膜采用玻璃纖維紙帶,連續(xù)測試17組。試驗(yàn)結(jié)果如圖1所示。
圖1 抽氣壓力引起的β射線強(qiáng)度變化曲線
由以上試驗(yàn)結(jié)果可知,在抽取零氣且抽氣時(shí)間短的情況下,濾膜并未采集到粉塵,此時(shí)空氣中粉塵濃度應(yīng)該為0,濾膜質(zhì)量前后也應(yīng)該相等。但是,如果以現(xiàn)有式(7)進(jìn)行計(jì)算,空氣中的粉塵濃度顯然不為0,同理,當(dāng)β射線檢測儀設(shè)備測試的大氣中粉塵濃度較低時(shí),其計(jì)算結(jié)果也會(huì)明顯偏大。由以上分析可見,現(xiàn)有的計(jì)算公式與實(shí)際實(shí)驗(yàn)結(jié)果存在著明顯的矛盾。
m=m′
(8)
(9)
(10)
由以上公式可以求出抽氣時(shí)負(fù)壓對濾膜質(zhì)量吸收系數(shù)產(chǎn)生的影響因子ε。對于不同的β射線自動(dòng)監(jiān)測機(jī)構(gòu),由于系統(tǒng)負(fù)壓以及濾膜的材質(zhì)的不同,其影響因子ε的不同,因此需要重復(fù)以上試驗(yàn)并利用式(8)~式(10)計(jì)算相應(yīng)的影響因子。
由以上分析可知,濾膜在抽塵前后質(zhì)量吸收系數(shù)會(huì)有明顯的改變,而不同的β射線檢測設(shè)備由于設(shè)計(jì)流量、系統(tǒng)阻力以及使用紙帶的品牌不同,其影響因子ε也會(huì)有所不同。因此,每臺β射線粉塵濃度檢測設(shè)備在建立粉塵濃度計(jì)算方法之前都需要提前對影響因子ε進(jìn)行計(jì)算,從而更為準(zhǔn)確的消除其產(chǎn)生的誤差。
在求得濾膜質(zhì)量吸收系數(shù)影響因子之后,將其引入到抽塵前濾膜質(zhì)量的計(jì)算公式中:
(11)
抽塵后的濾膜質(zhì)量計(jì)算公式仍然為式(5),而濾膜增重公式修改為式(12)。
(12)
相應(yīng)的,大氣中粉塵濃度計(jì)算公式修改為式(13)。
(13)
至此,新的粉塵濃度計(jì)算公式推導(dǎo)完成,該公式考慮到了濾膜質(zhì)量吸收系數(shù)μm變化對粉塵濃度計(jì)算的影響。在實(shí)際使用中對比現(xiàn)有的計(jì)算公式需要對β射線未穿透濾膜時(shí)的放射強(qiáng)度I0進(jìn)行測試。而I0在放射源強(qiáng)度不變,且接收器與放射源距離不變的情況下是穩(wěn)定的。因此,該公式在實(shí)際運(yùn)用中具有可行性。
將β射線粉塵濃度檢測系統(tǒng)與標(biāo)準(zhǔn)稱重法采樣器進(jìn)行同步采樣,期間記錄β射線的放射強(qiáng)度,分別使用新粉塵濃度計(jì)算公式和現(xiàn)有的粉塵濃度計(jì)算公式計(jì)算大氣中的粉塵濃度,并對稱重采樣器采集的濾膜進(jìn)行后期分析處理。實(shí)驗(yàn)地點(diǎn)為重慶市沙坪壩某小區(qū),采樣地點(diǎn)距離地面3 m,每次采樣時(shí)間為1 h,對大氣進(jìn)行間斷的對比采樣測試。實(shí)驗(yàn)歷時(shí)2個(gè)月,期間天氣變化較大,得到了32組試驗(yàn)數(shù)據(jù),如圖2所示,結(jié)果中粉塵濃度跨度較大,具有一定的代表性。
圖2 粉塵濃度對比實(shí)驗(yàn)結(jié)果
由實(shí)驗(yàn)結(jié)果可知,兩種公式計(jì)算出的粉塵濃度與標(biāo)準(zhǔn)稱重法測試結(jié)果都比較接近,但是現(xiàn)有計(jì)算公式整體計(jì)算結(jié)果比新公式計(jì)算結(jié)果偏高,特別是在粉塵濃度低于20 μg/m3時(shí),現(xiàn)有計(jì)算公式的粉塵濃度相對于標(biāo)準(zhǔn)稱重法的計(jì)算結(jié)果偏高34%~71%,而新公式計(jì)算的粉塵濃度在低粉塵濃度時(shí)與標(biāo)準(zhǔn)稱重法的偏差在-5%~12%之間。新公式與現(xiàn)有公式計(jì)算的粉塵濃度與標(biāo)準(zhǔn)稱重法測量結(jié)果的相關(guān)系數(shù)分別為0.986和0.975,平均測試誤差由8%減小到了3%。
①從β射線衰減原理入手,分析了影響濾膜質(zhì)量吸收系數(shù)變化的原因,并提出了濾膜質(zhì)量吸收系數(shù)μm的影響因子ε的算法。
②在現(xiàn)有β射線粉塵濃度計(jì)算公式的基礎(chǔ)上引入了質(zhì)量吸收系數(shù)的影響因子,推導(dǎo)出了新的β射線粉塵濃度計(jì)算公式。
③與標(biāo)準(zhǔn)稱重方法的參比方法對比實(shí)驗(yàn)表明,新公式在計(jì)算低濃度大氣粉塵濃度時(shí)比現(xiàn)有公式具有更高的精度,并且將β射線法與標(biāo)準(zhǔn)稱重法測試結(jié)果的相關(guān)系數(shù)由0.975提高到了0.986,平均誤差由8%減小到了3%。。
[1] 狄文靜. 基于β射線法的PM2.5連續(xù)在線監(jiān)測的相關(guān)問題研究[D]. 天津:天津大學(xué),2013.
[2] Mousa Vian S M,Najafi A F. Numerical Simulation of Gas-Liquid-Solid Nows in a Hydro Cyclone Separator[J]. Archive of Applied Mechanics,2009,79(5):395-405.
[3] Matsuda K,Fujimura Y,Hayashi K. Deposition Velocity of PM2.5 Sulfate in the Summer Above a Deciduous Forest in Central Japan[J]. Atmospheric Environment,2010,44(36):4582-4587.
[4] 譚韋君,丁萬山. 基于紅外傳感器和ARM的大氣有害氣體濃度監(jiān)測系統(tǒng)[J]. 傳感技術(shù)學(xué)報(bào),2011,24(3):321-324.
[5] 張加宏,韋圓圓,顧芳,等. 高精度光散射氣溶膠質(zhì)量濃度測量系統(tǒng)的信號處理研究[J]. 傳感技術(shù)學(xué)報(bào),2016,29(4):536-544.
[6] HJ 618—2011. 《環(huán)境空氣PM10和PM2.5的測定重量法》. 國家環(huán)境保護(hù)標(biāo)準(zhǔn).
[7] Michiels H,Deutsch F,De Nocker L,et al. Human Health Impacts of PM2.5 and NOXTransport Air Pollution in Belgium[M]. NATO Science for Peace and Security Series C:Environmental Security,2012:565-570.
[8] 吳佳濤,陳文亮,張?jiān)龈?等. 基于β射線技術(shù)測量PM2.5的準(zhǔn)確性研究[J]. 安全與環(huán)境學(xué)報(bào). 2016,16(2):338-342.
[9] 左鍵,趙政. 基于β射線吸收法的PM2.5連續(xù)自動(dòng)監(jiān)測的關(guān)鍵技術(shù)研究[J]. 儀表技術(shù)與傳感器,2017(1):139-146.
[10] Vladutescu D V,Wu Yonghua,Gross B M,et al. Remote Sensing Instruments Used for Measurement and Model Validation of Optical Parameters of Atmospheric Aerosols[J]. IEEE Transaction on Instrumentation and Measurement,2012,61(6):1733-1746.
[11] 張玉鈞,劉文清,鄭朝暉,等. β射線自動(dòng)監(jiān)測大氣中PM10含量的設(shè)計(jì)方案[J]. 光子學(xué)報(bào),2001,30(3):369-374.
[12] 梁艷,張?jiān)龈?陳文亮,等. 基于β射線法的新型PM2.5自動(dòng)監(jiān)測系統(tǒng)研究[J]. 傳感技術(shù)學(xué)報(bào),2014,27(10):1418-1422.
[13] 練德幸. 基于PLC的β射線吸收法粉塵濃度監(jiān)測技術(shù)研究[D]. 湖南:南華大學(xué).
[14] 趙鑫,潘晉孝,劉賓,等. 基于β射線吸收法的PM2.5測量技術(shù)的研究[J]. 電子技術(shù)應(yīng)用,2013,39(9):74-80.