煙塵濃度測(cè)量方法綜述

李 昆，鐘 磊，張洪泉

(1.哈爾濱工程大學(xué)，黑龍江哈爾濱150001;2.中國電子科技集團(tuán)公司第四十九研究所，黑龍江哈爾濱150001)

0 引言

進(jìn)入21世紀(jì)，人類生產(chǎn)生活的各個(gè)方面都會(huì)產(chǎn)生煙塵，鋼鐵和有色金屬冶煉、火力發(fā)電、水泥和石油化工企業(yè)的生產(chǎn)過程，煤礦、采石場(chǎng)等礦山生產(chǎn)中產(chǎn)生的煙塵，車輛和飛機(jī)排放的廢氣，以及垃圾焚燒處理、家庭爐灶、供暖鍋爐排出的煙氣等，都是煙塵污染的主要來源，其中，以燃料燃燒排出的數(shù)量最大，主要成分是未燃燒的碳粒(C)，還含有少量 SiO2，Al2O3，F(xiàn)e2O3，CaO 等。煙塵污染大氣，對(duì)人們身體健康有很大的危害，會(huì)引起心臟病患者死亡率的增加。另一方面，隨著低碳環(huán)保行動(dòng)的深入，降低生活環(huán)境中C的含量成為環(huán)境保護(hù)的重中之重。因此，對(duì)排放源煙塵濃度的測(cè)量就成為環(huán)境監(jiān)測(cè)的一個(gè)重要方面［1］。

1 測(cè)量方法

目前，根據(jù)測(cè)量機(jī)理的不同分為兩類分析方法:取樣法和非取樣法［2］。

1.1 取樣法

取樣法是從待測(cè)區(qū)域中取部分具有代表性的含煙塵氣樣，并將顆粒從樣品中分離出來，再送入隨后的分析測(cè)量系統(tǒng)來測(cè)量煙塵質(zhì)量濃度的方法［3］。

1)濾膜稱重法

濾膜稱重法的基本原理是以規(guī)定的流量采樣，將空氣中的煙塵顆粒沉集于高性能濾膜上，稱濾膜采樣前、后的質(zhì)量，由質(zhì)量差求得沉集的煙塵顆粒質(zhì)量，再根據(jù)采樣空氣體積，計(jì)算出煙塵顆粒的質(zhì)量濃度。

由于受濾膜性能影響，大多測(cè)量采用PM10和PM2.5 2個(gè)標(biāo)準(zhǔn)的煙塵顆粒物。如，武漢分析儀器廠FC—2B型粉塵采樣器，采用濾膜稱重原理，采集作業(yè)場(chǎng)所空氣的煙塵，適用于冶金、礦山、陶瓷等有粉塵危害的作業(yè)場(chǎng)所定點(diǎn)采樣，流量達(dá)到25L/min，精度為2.5級(jí)［4］，時(shí)間精度小于2s。

該方法原理簡(jiǎn)單，測(cè)定數(shù)據(jù)可靠，測(cè)量不受顆粒物物理性質(zhì)的影響。但操作煩瑣費(fèi)時(shí)(一般3～24 h)、噪聲大。

2)β射線吸收法

β射線吸收法測(cè)量裝置由β射線源、濾膜支架及探測(cè)器等組成。當(dāng)含塵樣氣通過濾膜時(shí)，顆粒被過濾在濾膜上，經(jīng)過一段時(shí)間后，轉(zhuǎn)動(dòng)軸帶動(dòng)濾膜移動(dòng)并使被濾顆粒進(jìn)入測(cè)量區(qū)域，測(cè)量區(qū)域上部發(fā)出的β射線透過顆粒介質(zhì)后衰減并被接收，根據(jù)β射線的衰減程度即可確定被濾塵樣的質(zhì)量，進(jìn)而求得被測(cè)粉塵的質(zhì)量濃度。

β射線吸收法是在稱重法基礎(chǔ)上發(fā)展而來的，該方法主要用于煤礦粉塵與工業(yè)燃燒煙塵(主要含C和S)的測(cè)量，以及用于氣溶膠質(zhì)量濃度的監(jiān)測(cè)。2009年，懷化鐵路疾病預(yù)防控制中心的向純海對(duì)β射線吸收法測(cè)量作業(yè)場(chǎng)空氣粉塵(炭黑粉塵)濃度做研究［5］，結(jié)果表明:β射線吸收法反演出的粉塵質(zhì)量與濾膜增重Δm基本一致，濾膜增重與β射線吸收度呈正相關(guān)，且隨采樣時(shí)間t和濾膜增重Δm的增加，β射線吸收度的變異系數(shù)和炭黑粉塵濃度的變異系數(shù)總體趨勢(shì)是減小，當(dāng)采樣時(shí)間 t≥7 min，濾膜增重 Δm≥0.3 mg時(shí)，2種變異系數(shù)均在10%以下。檢測(cè)下限為增重Δm為0.0002 mg。典型代表儀器武漢天虹儀器集團(tuán)的TH—25型大氣顆粒物濃度監(jiān)測(cè)儀，其利用低能β射線的輻射吸收原理，對(duì)空氣總懸浮顆粒物(TSP)、飄塵(PM10)、呼吸性顆粒物(PM5和PM2.5)進(jìn)行濃度監(jiān)測(cè)分析，采樣流量為16.7 L/min，準(zhǔn)確度3%，測(cè)量濃度為0～5.0 mg/m3，重復(fù)性為2%。

該方法測(cè)量的動(dòng)態(tài)范圍寬，準(zhǔn)確度及靈敏度高，且測(cè)量結(jié)果只與粒子的質(zhì)量有關(guān)。但該方法存在安全隱患，同時(shí)，系統(tǒng)要求增加各種屏蔽措施，結(jié)構(gòu)設(shè)備復(fù)雜且昂貴。

3)壓電晶體差頻法

壓電晶體法采用石英諧振器作為敏感元件。其工作原理是使空氣以恒定流量通過切割器，進(jìn)入由高壓放電針和微量石英諧振器組成的靜電采樣器，在高壓電暈放電的作用下，氣流中的顆粒物全部沉降于測(cè)量諧振器的電極表面上，因電極上增加了顆粒物的質(zhì)量，其振蕩頻率發(fā)生變化，根據(jù)頻率變化可測(cè)定煙塵顆粒物的質(zhì)量濃度。

與其他測(cè)煙塵濃度的方法相比，壓電晶體差頻法具有靈敏度很高，石英壓電晶體電極的質(zhì)量靈敏度理論上為180 Hz/μg。假定所測(cè)空氣煙塵濃度為 150 μg/m3，以1 L/min采樣流速采樣2 min，所采煙塵量為0.3 mg，儀器的理論響應(yīng)值為54 Hz，就可準(zhǔn)確測(cè)定。其采樣流量低、采樣時(shí)間短是其他測(cè)塵法無法比擬的;檢測(cè)范圍寬，由于輸出的低頻信號(hào)達(dá)10～105Hz，能滿足大氣煙塵不同濃度的測(cè)定，如果輸出104Hz，被測(cè)煙塵濃度理論上等于28 mg/m3，這樣高的煙塵濃度一般已超出環(huán)境煙塵濃度的范圍。

典型儀器是日本加野麥克斯的3511粉塵計(jì)，該儀器的核心是壓電晶體。將浮游粉塵收集到壓電晶體上，通過石英頻率的變化測(cè)試粉塵質(zhì)量，測(cè)量濃度范圍為0.02～10 mg/m3，靈敏度為 0.005 μg/Hz。

但由于壓電晶體每做完一次測(cè)試后需要重新清潔后才能進(jìn)行下次測(cè)試，所以，這種測(cè)試方法不能進(jìn)行長時(shí)間在線檢測(cè)。

4)微量天平振蕩法

測(cè)量原理是基于錐形元件振蕩微量天平原理，核心部件為錐形元件振蕩器。錐形元件振蕩器在其自然頻率下振蕩，振蕩頻率由振蕩器件的物理特性、參加振蕩的濾膜質(zhì)量和沉積在濾膜上的顆粒物質(zhì)量決定。儀器通過采樣泵和流量計(jì)，使環(huán)境空氣以一恒定的流量通過采樣濾膜，顆粒物則沉積在濾膜上。測(cè)量出一定間隔時(shí)間前、后的2個(gè)振蕩頻率，就能計(jì)算出在這一段時(shí)間里收集在濾膜上顆粒物的質(zhì)量，再除以流過濾膜的空氣的總體積，得到這段時(shí)間內(nèi)空氣中顆粒物的平均濃度。

以本方法為測(cè)量機(jī)理的代表儀器是美國RP公司的TEOM系列RP—1400 a監(jiān)測(cè)儀，該儀器直接和實(shí)時(shí)測(cè)量室內(nèi)(外)環(huán)境中小于10μm的煙塵質(zhì)量濃度。其技術(shù)特點(diǎn)是高精度，小時(shí)平均質(zhì)量濃度為±1.5μg/m3，24 h平均質(zhì)量濃度為 ±0.5μg/m3，高分辨率為 0.01μg/m3。

微量天平振蕩法適用范圍很廣［6］，現(xiàn)代主要用于空間環(huán)境表面污染(分子污染和顆粒物污染)的監(jiān)測(cè)，又因其高靈敏度、高分辨率及實(shí)時(shí)在線監(jiān)測(cè)、輸出數(shù)字化等優(yōu)點(diǎn)在電化學(xué)和生物領(lǐng)域備受關(guān)注。

1.2 非取樣法

非取樣法是利用煙塵顆粒物與射線、光等作用后所產(chǎn)生的衰減、散射等現(xiàn)象來間接測(cè)量煙塵濃度的方法。非取樣法主要有:黑度法、濁度法、光散射法。

1)黑度法［7］

此方法又叫林格曼黑度法。它是基于監(jiān)測(cè)人員用有不同黑色面積的玻璃片對(duì)排放煙塵的黑度進(jìn)行目測(cè)，然后與林格曼黑度(共分六級(jí))對(duì)比后，確定被測(cè)煙塵的黑度，再按林格曼黑度級(jí)與煙塵濃度對(duì)照表得到煙塵排放濃度。這種方法使用簡(jiǎn)單、方便，操作人員很容易掌握使用，但顯然這種方法不夠科學(xué)，也不夠可靠，無法獲得煙塵的絕對(duì)濃度。

以黑度法為測(cè)量機(jī)理的儀器，如青島高通分析儀器有限公司的LGM—A1型林格曼煙塵黑度計(jì)，其測(cè)量精度為0.5%。

黑度法只是粗略了解煙塵的黑度等級(jí)而不需要獲得其絕對(duì)濃度，主要用于煙塵黑度監(jiān)測(cè)，應(yīng)用于鍋爐、工業(yè)爐窯、火電廠及煉焦?fàn)t等場(chǎng)所［8］。

2)濁度法

該方法將光源與探測(cè)器分別安裝在煙道兩側(cè)，光遇到煙塵顆粒后由于吸收、散射等作用使光強(qiáng)衰減，探測(cè)器接收的是顆粒的透射光。根據(jù)郎伯—比爾定律，透射光強(qiáng)與顆粒的大小和濃度相關(guān)，這就為煙塵顆粒物濃度測(cè)量提供了尺度，通過計(jì)算介質(zhì)的濁度，得到煙塵的質(zhì)量濃度。

該方法原理簡(jiǎn)單、技術(shù)成熟，廣泛用于工業(yè)煙囪、煤礦瓦斯監(jiān)測(cè)，但用于濃度測(cè)量時(shí)必須預(yù)先知道被測(cè)對(duì)象的粒徑分布或者平均粒徑，具有一定的局限性，即在濃度極低時(shí)，光強(qiáng)變化不大，濃度極高時(shí)，光強(qiáng)衰減過大，從而信噪比大大降低，因此，在這種特殊情況下，效果較差;當(dāng)煙塵組分發(fā)生變化時(shí)，測(cè)量結(jié)果也會(huì)出現(xiàn)偏差;由于光源、探測(cè)器及反射鏡等需要分立安裝，因此，需要嚴(yán)格對(duì)準(zhǔn);反射鏡等光學(xué)鏡片附著煙塵后，也會(huì)影響測(cè)量結(jié)果。

3)光散射法［9］

光散射法基于光散射原理，當(dāng)光束入射到顆粒(不管是固體顆粒、液滴或者氣泡)上時(shí)，將向空間四周散射，光的各個(gè)散射參數(shù)與煙塵顆粒的濃度密切相關(guān)。將探測(cè)器安裝在某一散射角處，獲得散射光強(qiáng)數(shù)據(jù)后，基于散射理論對(duì)煙塵濃度進(jìn)行反演。

代表儀器有LD—3F防爆袖珍型電腦激光粉塵儀，可以直接讀取質(zhì)量濃度，適用于室內(nèi)外環(huán)境中可吸入顆粒物(PM10)濃度的檢測(cè)，測(cè)定范圍為0.01～100 mg/m3或者0.001～10 mg/m3;靈敏度為0.01或者0.001mg/m3;測(cè)量精度為±10%。

光散射法之所以獲得廣泛應(yīng)用是因?yàn)橄啾绕渌麥y(cè)量方法具有如下顯著優(yōu)點(diǎn):適用性廣，除了測(cè)量固體顆粒(粉末)外，還可以測(cè)量液體顆粒(液滴)、氣體顆粒(氣泡)，而不用知道顆粒的化學(xué)組成;粒徑測(cè)量范圍寬，從幾個(gè)納米(10-3μm)到約103μm，甚至更大;測(cè)量準(zhǔn)確、精度高、重復(fù)性能好，對(duì)單分散系高分子聚合物標(biāo)準(zhǔn)粒子的測(cè)量誤差和重復(fù)性偏差可以限制在1%～2%之內(nèi)。

2 發(fā)展概況

2.1 國外研究情況

1957年，Van de Hulst H C的專著“Light Scattering by Small Particles”出版［10］，成為粒子散射理論的經(jīng)典著作。

20世紀(jì)70年代起，激光雷達(dá)發(fā)展迅速，激光雷達(dá)遙感理論日趨成熟。1983年，F(xiàn)itzgerald基于Mie散射理論分析了低對(duì)流層氣溶膠粒子的相對(duì)濕度對(duì)0.694，10.6μm兩個(gè)波段后向散射系數(shù)的影響［11］。1994年，美國國家航空航天局利用“發(fā)現(xiàn)號(hào)”航天飛機(jī)利用激光雷達(dá)進(jìn)行了大氣探測(cè)實(shí)驗(yàn)并獲得成功，使激光雷達(dá)的探測(cè)范圍大大拓寬［12］。

進(jìn)入21世紀(jì)，煙塵濃度在線監(jiān)測(cè)和大氣污染防治成為焦點(diǎn)，光后向散射法研究的內(nèi)容主要集中在粒子后向散射特性的理論計(jì)算和實(shí)時(shí)測(cè)量的實(shí)現(xiàn)2個(gè)方面。2001年，Boss Emmanuel等人在實(shí)驗(yàn)測(cè)量的基礎(chǔ)上利用體散射函數(shù)對(duì)雨滴的后向散射系數(shù)進(jìn)行了反演，為理論計(jì)算提供了有力的支持。2002年，Tishkovets V P對(duì)薄層粒子系的后向散射理論進(jìn)行了詳細(xì)的推導(dǎo)和計(jì)算，為粒子系的后向散射理論的發(fā)展奠定了基礎(chǔ)。2004年，Griaznov V等人基于Mie散射模型，對(duì)水滴球形粒子含偏心包含物的后向散射特性進(jìn)行了分析，為雷達(dá)亮暗帶現(xiàn)象出現(xiàn)的原因提供了理論依據(jù)［13］。2006年，Michael Mishchenko I對(duì)這一階段有關(guān)后向散射增強(qiáng)等一系列粒子散射問題進(jìn)行了梳理和總結(jié)，并指出了未來的研究方向［14］。2008年，Budak Vladimir P發(fā)展了一種新的邊界處理方法，對(duì)高濃度各向異性的三維粒子系的散射特性進(jìn)行了計(jì)算，進(jìn)一步拓展了散射理論［15］。2009年，Chaikovskaya Ludmila I利用CCD雷達(dá)獲得的實(shí)驗(yàn)數(shù)據(jù)反演多級(jí)散射粒子系的后向散射特性，直接獲得了粒子散射系的極化圖和 Mueller矩陣圖，信息更加豐富［16］。2010年，Sommersten Endre R利用離散坐標(biāo)法和Monte Carlo法對(duì)兩平行層不同粒子系的散射場(chǎng)進(jìn)行了計(jì)算，結(jié)果表明:該方法計(jì)算結(jié)果更加精確，計(jì)算時(shí)間也大大縮短［17］。

2.2 國內(nèi)研究情況

國內(nèi)起步較晚，最早是1979年，上海市激光技術(shù)研究所的孫渝生運(yùn)用后向散射激光多普勒測(cè)速原理，對(duì)管內(nèi)微粒懸浮液自由下落的速度進(jìn)行了測(cè)量［18］。1997年，復(fù)旦大學(xué)的黃興忠和金亞秋，計(jì)算了群聚球形粒子和群聚非球形粒子的極化散射場(chǎng)，得出粒子之間的相干使后向散射明顯增強(qiáng)的結(jié)果。同年，國家海洋局的夏達(dá)英等利用后向散射熒光計(jì)對(duì)試驗(yàn)海區(qū)內(nèi)人工投放的粉煤灰懸浮顆粒的濃度進(jìn)行了測(cè)量，結(jié)果表明:后向散射法優(yōu)于前向法和衰減法［19］。

近年來，我國在顆粒物后向散射理論研究與實(shí)驗(yàn)方面已經(jīng)取得顯著的成果。2001年，上海理工大學(xué)蔡小舒等人對(duì)高濃度顆粒系的后向散射特性進(jìn)行了計(jì)算，指出當(dāng)顆粒濃度增大到一定程度后，決定后向散射強(qiáng)度的是顆粒的粒徑，而與濃度基本無關(guān)［20］。同年，濟(jì)南大學(xué)的王玉茹等人設(shè)計(jì)并搭建了光散射法測(cè)定煙塵濃度的實(shí)驗(yàn)裝置，并對(duì)典型煙塵樣品進(jìn)行了測(cè)量，結(jié)果表明:相對(duì)于光吸收法而言，光散射法具有更高的靈敏度［21］。2003年，北京理工大學(xué)的李麗等人采用輻射傳輸理論建立了水下距離選通脈沖激光成像系統(tǒng)的后向散射光模型，根據(jù)典型系統(tǒng)參數(shù)計(jì)算并繪制了系統(tǒng)的重疊系數(shù)和后向散射光變化曲線［22］。同年，天津大學(xué)精密測(cè)試技術(shù)及儀器國家重點(diǎn)實(shí)驗(yàn)室的蘇翼雄等人通過研究線偏振光照射在生物組織表面時(shí)其散射光的偏振態(tài)特性，提出了一種非接觸方法進(jìn)行生物組織后向散射光信號(hào)拾取的方法，為提高組織成分光譜檢測(cè)的信噪比提供了實(shí)驗(yàn)基礎(chǔ)。2004年，南京理工大學(xué)的谷亞林等人推導(dǎo)并計(jì)算了有限長導(dǎo)體圓柱的后向散射極化矩陣，分析了目標(biāo)的特征極化參數(shù)和入射角之間的關(guān)系［23］。2006年，華中科技大學(xué)的鄧勇等人基于Mie散射理論，建立了雙層散射介質(zhì)的單次后向散射光譜的理論模型，結(jié)果表明:利用后向散射光對(duì)粒徑大小、分布和復(fù)折射率的高靈敏度，可對(duì)實(shí)現(xiàn)癌細(xì)胞的早期識(shí)別［24］。2009年，西安電子科技大學(xué)的吳振森等人利用Monte-Carlo方法對(duì)電磁波在煙塵中的傳輸特性進(jìn)行數(shù)值計(jì)算，給出煙塵中反射率和透射率隨入射角和煙塵厚度變化的數(shù)值結(jié)果［25］。2010年，中科院上海光機(jī)所的王向偉等人基于單個(gè)氣泡的米氏散射理論，利用多層介質(zhì)蒙特卡—羅模擬方法，對(duì)I，IA，IB類海水下的尾流氣泡激光后向散射特性進(jìn)行了模擬，證明了可利用尾流氣泡群的激光后向散射作為魚雷的引信［26］。

3 主要發(fā)展方向

1)計(jì)算煙塵粒子散射特性的光學(xué)常數(shù):目前光學(xué)常數(shù)的測(cè)量主要采用KBr樣片透射法，該方法操作復(fù)雜、被測(cè)粒徑受限，且只能通過紅外波段的結(jié)果外推至可見光波段。

2)計(jì)算復(fù)雜形貌煙塵粒子的后向散射特性:目前構(gòu)建的煙塵粒子模型僅為特殊的非球模型，如橢球、圓柱等。而實(shí)際煙塵粒子形貌各異、表面十分復(fù)雜，因此，按傳統(tǒng)方法難以準(zhǔn)確獲得實(shí)際煙塵粒子的后向散射特性。

3)計(jì)算煙塵粒子系的后向散射特性:目前的計(jì)算是基于Mie散射理論進(jìn)行疊加或復(fù)散射計(jì)算，這僅適用于球形顆粒系，而難以滿足煙塵粒子排放現(xiàn)狀。

4 結(jié)束語

煙塵濃度的在線測(cè)量方法應(yīng)滿足以下要求:

1)采樣速度要足夠快，能進(jìn)行長期、實(shí)時(shí)地監(jiān)測(cè)，能夠滿足生產(chǎn)過程中對(duì)數(shù)據(jù)量的要求。

2)數(shù)據(jù)處理必須實(shí)時(shí)準(zhǔn)確，及時(shí)反映出煙塵顆粒物特性的變化。

3)測(cè)量系統(tǒng)結(jié)構(gòu)簡(jiǎn)單、可靠，能夠在惡劣條件下長期運(yùn)行，便于維護(hù)。

4)測(cè)量系統(tǒng)還應(yīng)具有較好的經(jīng)濟(jì)性，價(jià)格合理。

［1］王 瑋，劉紅杰.TSP—PM10—PM2.5—2型中流量大氣顆粒物采集系統(tǒng)的開發(fā)和應(yīng)用［J］.中國環(huán)境監(jiān)測(cè).2003，19(1):35-37.

［2］王乃寧.顆粒粒徑的光學(xué)測(cè)量技術(shù)及應(yīng)用［M］.北京:原子能出版社，2000:25-30.

［3］GB 3095—1995，空氣環(huán)境質(zhì)量標(biāo)準(zhǔn)［S］.

［4］GB/T5817—2009，粉塵作業(yè)場(chǎng)所危害程度分級(jí)［S］.

［5］向純海.β射線吸收法測(cè)量鐵路作業(yè)場(chǎng)所空氣常見粉塵濃度［J］.鐵道勞動(dòng)安全衛(wèi)生與環(huán)保，2009，36(6):316-319.

［6］李 萍.15MHz溫控石英晶體微量天平的設(shè)計(jì)［D］.蘭州:蘭州大學(xué)，2009:6-13.

［7］包信宗，葛廣英，張衛(wèi)華，等.基于BP的煙氣林格曼黑度級(jí)測(cè)量研究［J］.聊城大學(xué)學(xué)報(bào):自然科學(xué)版，2009，22(2):22-25.

［8］楊大惠，尹常慶，楊傳俊.煙氣黑度監(jiān)測(cè)中應(yīng)注意的問題［J］.環(huán)境監(jiān)測(cè)管理與技術(shù)，2005，17(3):38-39.

［9］張 敏.光散射法測(cè)量超細(xì)顆粒粒度的研究［D］.天津:天津大學(xué).2005:1-12.

［10］Van de Hulst H C.Light scattering by small particles［M］.John Wiley and Sons Inc，1957:1-10.

［11］Fitzgerald.Effect of relative humidity on the aerosol backscattering coefficient at 0.694 and 10.6μm wavelengths［J］.Applied Optics，1984，23(3):411-417.

［12］Kandlikar M，Ramachandran G.Inverse methods for analyzing aerosol spectrometer measurements:A critical review［J］.Journal of Aerosol Science，1999，30(4):413-437.

［13］Griaznov V.Numerical simulation of light backscattering by spheres with off-center inclusion:Application to the lidar case［J］.Applied Optics，2004，43(29):5512-5522.

［14］Michael Mishchenko I.Maxwell’s equations，radiative transfer，and coherent backscattering:A general perspective［J］.Journal of Quantitative Spectroscopy and Radiative Transfer，2006，101(3):540-555.

［15］Budak Vladimir P.On the solution of a vectorial radiative transfer equation in a arbitrary three-dimensional turbid medium with anisotropic scattering［J］.Journal of Quantitative Spectroscopy and Radiative Transfer，2008，109(2):220-234.

［16］Chaikovskaya Ludmila I.New technique to simulate angular patterns of polarized signals frommultiply scattering media in MFOV and CCD lidars［J］.Journal of Quantitative Spectroscopy and Radiative Transfer，2009，110(8):506-525.

［17］Sommersten Endre R.Discrete ordinate and Monte-Carlo simulations for polarized radiative transfer in a coupled system consisting of two media with different refractive indices［J］.Journal of Quantitative Spectroscopy and Radiative Transfer，2010，111(4):616-633.

［18］孫渝生.后向散射激光多普勒測(cè)速裝置［J］.激光，1979(5):34-38.

［19］黃興忠，金亞秋.群聚球形粒子極化散射的 T矩陣數(shù)值模擬［J］.計(jì)算物理，1997，14(1):376-381.

［20］俞大勇，蔡小舒.高濃度超細(xì)顆粒系后向散射研究［J］.中國粉體技術(shù)，2001，7(4):1-3.

［21］王玉茹，馮德振，張智杰.光散射法煙塵濃度的實(shí)驗(yàn)研究［J］.濟(jì)南大學(xué)學(xué)報(bào):自然科學(xué)版，2001，15(9):253-255.

［22］李 麗.水下距離選通成像系統(tǒng)后向散射光的計(jì)算［J］.北京理工大學(xué)學(xué)報(bào)，2003，23(4):488-491.

［23］谷亞林.有限長導(dǎo)體圓柱的后向散射極化分析［J］.電波科學(xué)學(xué)報(bào)，2004，19(3):269-273.

［24］鄧 勇.雙層散射介質(zhì)的單次后向散射光譜分析［J］.光學(xué)學(xué)報(bào)，2006，26(4):595-599.

［25］吳振森.煙塵中電磁波傳輸特性的Monte-Carlo模擬［J］.物理學(xué)報(bào)，2009，58(4):2397-2404.

［26］王向偉.艦船尾流激光后向散射特性研究［J］.光學(xué)學(xué)報(bào)，2010，30(1):14-18.