單 健 劉曉翔 肖德濤 劉 衛(wèi) 劉 彥

1（南華大學(xué) 核科學(xué)技術(shù)學(xué)院 氡湖南省重點(diǎn)實(shí)驗(yàn)室 衡陽 421001）2（中國科學(xué)院上海應(yīng)用物理研究所 嘉定園區(qū) 上海 201800）

大氣顆粒物中210Pb的活度-粒徑分布測(cè)量

單 健1劉曉翔1肖德濤1劉 衛(wèi)2劉 彥1

1（南華大學(xué) 核科學(xué)技術(shù)學(xué)院 氡湖南省重點(diǎn)實(shí)驗(yàn)室 衡陽 421001）2（中國科學(xué)院上海應(yīng)用物理研究所 嘉定園區(qū) 上海 201800）

用M235大氣顆粒物分級(jí)采樣器采樣獲得不同范圍粒徑的可吸入大氣顆粒物樣品，然后分別用β分析方法和γ分析方法測(cè)量不同粒徑尺度樣品中的210Pb的活度分布。得到大氣顆粒物樣品中210Pb的體積比活度在不同粒徑上的分布范圍為0.004 4–1.832 1mBq·m?3，最大值均落在0.49–0.95μm粒徑范圍。兩種分析方法所測(cè)得的數(shù)據(jù)具有近似的活度-粒徑分布關(guān)系，其離散性變化規(guī)律也基本一致，都是在粒徑范圍為0.49–1.5μm之間的標(biāo)準(zhǔn)差最大，β計(jì)數(shù)法測(cè)量的數(shù)據(jù)離散性整體比γ計(jì)數(shù)法測(cè)量的數(shù)據(jù)偏大。

210Pb活度，活度-粒徑分布，比活度

大氣顆粒物中210Pb是鈾系衰變鏈中的產(chǎn)物，其主要是由大氣中的222Rn衰變而來。210Pb在大氣環(huán)境中的含量不多，但是半衰期為22.26a，容易累積而導(dǎo)致其濃度升高。并且210Pb的子體210Bi和210Po會(huì)在一定時(shí)間內(nèi)與其達(dá)到放射性平衡，210Bi和210Po的半衰期相對(duì)較短，但是相對(duì)的放射性較強(qiáng)。放射衛(wèi)生界對(duì)大氣中222Rn短壽命子體的研究比較多，而對(duì)于222Rn長(zhǎng)壽命子體的研究很少，尤其是對(duì)于大氣中210Pb的研究更少[1–3]。

210Pb是一種示蹤核素，其可以被單獨(dú)運(yùn)用，或與7Be核素結(jié)合運(yùn)用于大氣環(huán)境污染監(jiān)測(cè)中。在一個(gè)理想的封閉空間中（相對(duì)于核素的封閉系統(tǒng)），可以從210Pb/222Rn、210Bi/210Pb和210Po/210Pb的活度比中得到氣溶膠的滯留時(shí)間。但大量的研究表明這三對(duì)活度比獲得的滯留時(shí)間不一致。從粒徑分級(jí)氣溶膠中這些核素的分布中看出，粒徑較小的氣溶膠滯留時(shí)間的計(jì)算結(jié)果一致，這可能意味著更大粒徑尺度的氣溶膠來自揚(yáng)塵。本研究使用不同的測(cè)量方法對(duì)大氣顆粒物中210Pb的活度粒徑分布情況進(jìn)行測(cè)量，為進(jìn)一步開展利用210Pb等自然產(chǎn)生的核素進(jìn)行大氣環(huán)境示蹤研究提供基礎(chǔ)[4–10]。

1 測(cè)量210Pb的方法及原理

1.1 γ分析方法及原理

通過γ分析法測(cè)定210Pb活度，只需要對(duì)樣品進(jìn)行簡(jiǎn)單的處理就可以放入高純鍺(HPGe) γ譜儀中測(cè)量。根據(jù)210Pb發(fā)射γ粒子在能量為46.5keV處產(chǎn)生的全能峰面積，以及刻度好的效率曲線和210Pb的γ發(fā)射幾率等數(shù)據(jù)可以確定210Pb活度。本實(shí)驗(yàn)測(cè)量使用了GMX50-N型(HPGe) γ譜儀，相對(duì)效率為50%，配套有反康普頓系統(tǒng)。主探測(cè)器由一個(gè)N型HPGe探頭，反符合由一個(gè)?22.86cm的環(huán)狀NaI探測(cè)器和一個(gè)?7.62cm的塞子探測(cè)器組成。210Pb釋放γ射線的能量處在低能區(qū)，實(shí)驗(yàn)測(cè)量中必須盡量降低本底和康普頓坪。因此選擇在低能區(qū)探測(cè)效率較高，并配套反康普頓譜儀的GMX50-N型(HPGe) γ譜儀作為實(shí)驗(yàn)裝置(圖1)。

圖1 GMX50-N型(HPGe) γ譜儀Fig.1 GMX50-N type (HPGe) γ spectrometer.

1.2 β分析方法及原理

由于210Pb自身發(fā)射的β粒子能量很低，最大能量分別是16.5keV和63keV，因此β譜分析方法是

通過測(cè)定其第一代子體210Bi發(fā)射的最大能量為1.16MeV的β射線來推算210Pb活度。以HNO3-HClO4-HCl消解體系破壞濾膜，使用穩(wěn)定核素206Pb為載體，利用強(qiáng)堿性陰離子交換樹脂實(shí)現(xiàn)樣品中210Pb與衰變子體210Bi-210Po分離，通過稱量PbSO4沉淀，確定全程化學(xué)回收率。將樣品置于干燥處放置一個(gè)月后，待210Pb-210Bi達(dá)成放射性平衡，在MDS-4低本底計(jì)數(shù)器上測(cè)量210Bi的β粒子活度，再用公式計(jì)算可間接得到210Pb的活度。另外，樣品中除含有210Pb外，還有214Pb與212Pb等其它Pb的同位素可能會(huì)形成干擾，將處理好的樣品放置一個(gè)月后，214Pb (26.8min)與212Pb (10.6 h)已基本衰變完，就不會(huì)對(duì)結(jié)果產(chǎn)生干擾了[10–11]。

2 大氣顆粒物樣品的采集與測(cè)量

2.1 大氣顆粒物樣品的采集

本實(shí)驗(yàn)選擇大流量的大氣顆粒物分級(jí)采樣器(圖2)，采樣流量大致為1.13m3·min?1，由于撞擊式的采樣方式在采樣的顆粒物達(dá)到一定量時(shí)，收集效率會(huì)明顯下降，所以通常將采樣時(shí)間設(shè)置為30–60h。采樣點(diǎn)周圍地勢(shì)開闊，無高大建筑物和工業(yè)廠房，周圍農(nóng)田、草地和樹木較多，幾乎沒有揚(yáng)塵，加上采樣點(diǎn)離地面有20多米，土壤揚(yáng)塵的影響可以忽略。采樣前注意天氣的變化趨勢(shì)，確保采樣期間內(nèi)沒有下雨。采樣膜必須提前一天放入恒溫恒濕箱體中作稱重前的平衡，然后稱量各級(jí)采樣膜的重量。采樣結(jié)束后又把采樣膜放入箱體中，在相同的溫濕度條件下，同樣平衡24h以上后再稱重。

圖2 M235沖擊式大流量大氣顆粒物分級(jí)采樣器的切割頭(a)和泵(b)Fig.2 Cutting head (a) and the pump (b) of M235 high flow impact atmospheric particulates grade sampler.

2.2 γ分析方法的實(shí)驗(yàn)過程

用蒸餾水清洗實(shí)驗(yàn)容器，把稱重后的采樣濾膜小心地剪碎于燒杯中。在燒杯中加入王水、2mL 高氯酸，把燒杯放到恒溫加熱器上，以 180°C的溫度加熱至冒白煙，重復(fù)以上過程兩次。用0.5mol·L?1的鹽酸溶解殘留物，然后用定量濾紙過濾樣品。把過濾后的溶液樣品轉(zhuǎn)移至專用的γ測(cè)量樣品盒內(nèi)，用蒸餾水定容為18.21g。密封γ測(cè)量樣品盒，完成210Pb測(cè)量樣品的制樣。

對(duì)于弱放射性樣品長(zhǎng)時(shí)間的測(cè)量是獲得足夠計(jì)數(shù)以減少統(tǒng)計(jì)誤差的一般方法，因此計(jì)數(shù)系統(tǒng)的長(zhǎng)時(shí)間穩(wěn)定性很重要，實(shí)驗(yàn)選取的γ譜儀穩(wěn)定性較好，在44.08–48.72keV能量范圍的本底計(jì)數(shù)率低于0.0054情況下，實(shí)驗(yàn)樣品單次測(cè)量時(shí)間一般在6×104s以上，這樣基本可以保證足夠的計(jì)數(shù)，確保統(tǒng)計(jì)誤差較低。210Pb的γ放射性活度測(cè)量與計(jì)算過程如圖3所示。

圖3 210Pb的γ放射性活度測(cè)量與計(jì)算過程Fig.3 Measurement and calculation process of γ radioactivity of210Pb.

本實(shí)驗(yàn)選擇了規(guī)格為?70mm×30mm的聚乙烯圓柱形樣品盒，每次樣品定量為18.21g（樣品的液面高度0.5cm），這樣既能測(cè)量足夠多的樣品以達(dá)到盡可能高的計(jì)數(shù)率，又不至于因?yàn)闃悠诽穸饑?yán)重的自吸收現(xiàn)象。譜儀的本底水平是最關(guān)鍵的制約因素，本實(shí)驗(yàn)使用的測(cè)量系統(tǒng)本底較低，其中本底測(cè)量時(shí)間432000s中有431340.7s的活時(shí)間，能區(qū)為44.08–48.72keV的本底計(jì)數(shù)為2436，平均本底計(jì)數(shù)率為0.0056。針對(duì)18.21g的樣品（樣品的液面高度0.5cm）能量為46.5keV的γ射線，根據(jù)效率刻度軟件計(jì)算得到的探測(cè)效率為0.11563。最后，根據(jù)《GB/T 11743-2013 土壤中放射性核素的γ能譜分析方法》中的公式(式(1))計(jì)算得探測(cè)下限為0.108Bq。

式中，tb是本底譜測(cè)量時(shí)間，s；Nb是本底譜中相對(duì)應(yīng)于某一全能峰的本底計(jì)數(shù)；K是相對(duì)于置信度為0.95的取值，K=1.645。

2.3 β分析的實(shí)驗(yàn)過程

運(yùn)用β分析法測(cè)量210Pb的活度，需要先分析樣品中210Bi的活度。相對(duì)于γ分析法，樣品需經(jīng)過較為復(fù)雜的化學(xué)方法處理，才能對(duì)其進(jìn)行測(cè)量。

將濾膜取出，剪碎后裝入100mL燒杯中，準(zhǔn)確加入2mL Pb2+載體溶液化學(xué)以確定回收率。為充分破壞濾膜并使得到的溶液為鹽酸體系，以HNO3-HClO4-HCl消解體系破壞濾膜，再用漏斗過濾出含有樣品的清液，將濾液控制在80mL左右。將浸泡在去離子水中的201×7型強(qiáng)堿性陰離子交換樹脂裝入玻璃交換柱的儲(chǔ)液池中，玻璃交換柱的上下兩端均用玻璃纖維封堵。用100mL 1mol·L?1HCl溶液以約0.22mL·min?1（每10–15 s一滴）速度過柱后待用。將過濾好的樣品液轉(zhuǎn)入交換柱儲(chǔ)液池內(nèi)，控制過柱速度約為0.28mL·min?1（每15–20s一滴）。用140mL 1.0mol·L?1鹽酸溶液以相同速率繼續(xù)過柱，最后再加入80mL去離子水，用100mL燒杯接洗脫液。

將溶液加熱蒸干，然后加兩滴2mol·L?1H2SO4溶液，再加入1mL飽和Na2SO4溶液，立刻攪拌直至出現(xiàn)大量白色沉淀，然后將溶液靜置。取過濾用濾紙?jiān)?5°C恒溫箱烘干1h，用微量電子天平（精確度為0.00001g）立即稱量。沉淀采用可拆卸式不銹鋼漏斗抽濾，依次用去離子水和無水乙醇洗滌沉淀。完畢后，放入干燥箱45 °C烘干1h后稱量。抽濾前、后的濾紙質(zhì)量差即為回收的PbSO4沉淀質(zhì)量，與之前加入的鉛載體溶液中206Pb質(zhì)量相比較即可確定化學(xué)回收率。樣品放置于干燥處一個(gè)月，達(dá)到平衡后在MDS-4低本底計(jì)數(shù)器上計(jì)210Bi的β粒子活度，再利用公式可間接得到210Pb的活度。

Pb回收率可以按照式(2)進(jìn)行計(jì)算：式中，Y是標(biāo)準(zhǔn)源中Pb的化學(xué)回收率；ms是標(biāo)準(zhǔn)源中PbSO4的凈重，mg；Cc是Pb (PbSO4)在載體溶液中的質(zhì)量濃度，mg·mL?1；Vc是加入載體溶液的體積，mL。利用MDS-4低本底計(jì)數(shù)器（計(jì)數(shù)氣體：10%甲烷，90%氬氣）對(duì)已知活度的標(biāo)準(zhǔn)源（241Am、90Sr和210Pb標(biāo)準(zhǔn)源）測(cè)量720min，可以得到儀器的探測(cè)效率。

3 結(jié)果與討論

3.1210Pb放射性測(cè)量的原始值及分析

本研究用γ計(jì)數(shù)法測(cè)量了上海市嘉定區(qū)冬季的不同粒徑大氣顆粒物的樣品，分別得到了這些樣品中210Pb的活度含量，并計(jì)算了對(duì)應(yīng)的體積比活度，分析了其體積比活度的粒徑分布特征，得到總懸浮顆粒物(Total suspended particulate, TSP)。根據(jù)樣品的粒徑分級(jí)情況及測(cè)量的活度數(shù)據(jù)按粒徑范圍分組處理并求平均值或求和，得到如表1所示的數(shù)據(jù)。

表1 γ計(jì)數(shù)法測(cè)量上海市嘉定區(qū)冬季大氣顆粒物中210Pb的活度粒徑分布數(shù)據(jù)(mBq·m?3)Table 1 Activity-size distribution of210Pb in the atmospheric particulates measured by γ count method in Shanghai in winter.

本研究用β計(jì)數(shù)法測(cè)量了衡陽市蒸湘區(qū)不同粒徑大氣顆粒物的樣品，分別得到了這些樣品中210Pb的活度含量，并計(jì)算了對(duì)應(yīng)的體積比活度，分析了其體積比活度的粒徑分布特征。得到如表2所示的數(shù)據(jù)。其中HY-1樣品的采集時(shí)間為春季，天氣晴朗，日平均氣溫24 °C；HY-2、HY-3、HY-4樣品的采集時(shí)間為秋季，日平均氣溫約為24 °C。

本研究還于夏季在上海市嘉定區(qū)進(jìn)行了采樣，用β計(jì)數(shù)法測(cè)量數(shù)據(jù)如表3所示。采樣時(shí)天氣晴朗，日平均氣溫約為34 °C。

表2 β計(jì)數(shù)法測(cè)量衡陽市蒸湘區(qū)春秋季大氣顆粒物中210Pb的活度粒徑分布數(shù)據(jù)(mBq·m?3)Table 2 Activity-size distribution of210Pb in the atmospheric particulates measured by β count method in Hengyang in spring and fall.

表3 β計(jì)數(shù)法測(cè)量上海市嘉定區(qū)夏季大氣顆粒物中210Pb的活度粒徑分布數(shù)據(jù)(mBq·m?3)Table 3 Activity-size distribution of210Pb in the atmospheric particulates measured by β count method in Shanghai in summer.

3.2 γ計(jì)數(shù)法與β計(jì)數(shù)法測(cè)量結(jié)果的比較

表1?3中的數(shù)據(jù)表明，大氣顆粒物中210Pb體積比活度的最大值均落在0.49–0.95μm粒徑范圍，此粒徑范圍內(nèi)210Pb的體積比活度含量占TSP（指粒徑＜100μm的大氣顆粒物）中總含量均超過了37%，最高的是表3中的數(shù)據(jù)，為61.66%。三組數(shù)據(jù)中PM3.0中的210Pb占TSP的總含量分別為95.18%、96.18%、94.89%。也就是說，210Pb幾乎全部吸附在可吸入的大氣顆粒物上，可以認(rèn)為210Pb是很容易進(jìn)入人體并可能產(chǎn)生影響的。這可以在一定程度上證明，目前國際上的許多專家認(rèn)為，210Pb、210Bi、210Po之類的氣載放射性對(duì)人類健康的影響可能比先前估計(jì)的更大，這一觀點(diǎn)是正確的。

把分時(shí)間段的10批樣品的測(cè)量結(jié)果進(jìn)行比較，做出了如圖4所示的柱狀比較圖，從圖4中可以看出，利用γ譜測(cè)量方法得到的數(shù)據(jù)與低本底β計(jì)數(shù)法所得到的衡陽與上海兩地區(qū)樣品中210Pb活度隨粒徑分布具有幾乎相同的規(guī)律。

HY-2、HY-3、HY-4樣品是用γ計(jì)數(shù)法測(cè)量的同一時(shí)間段同一采樣點(diǎn)的樣品，SH-1、SH-2、SH-3、SH-4樣品是用β計(jì)數(shù)法測(cè)量的同一時(shí)間段同一采樣點(diǎn)的樣品。用實(shí)驗(yàn)標(biāo)準(zhǔn)差分別對(duì)這兩組數(shù)據(jù)進(jìn)行離散性分析得到圖5。從圖5中可以看出，兩種方法所測(cè)得的數(shù)據(jù)離散性比較一致，都是在粒徑范圍為0.49–1.5μm之間的標(biāo)準(zhǔn)差最大，在粒徑小于0.49μm或大于1.5μm的范圍內(nèi)標(biāo)準(zhǔn)差都比較小。

圖4 大氣顆粒物中210Pb體積比活度粒徑分布柱狀圖Fig.4 Activity-size distribution of210Pb in the atmospheric particulates.

圖5 γ計(jì)數(shù)法與β計(jì)數(shù)法的數(shù)據(jù)離散性對(duì)比Fig.5 γ counting discrete compared with the data β counting method.

4 結(jié)語

(1) 實(shí)驗(yàn)結(jié)果表明：上海市嘉定區(qū)冬季的大氣顆粒物中TSP的210Pb體積比活度為1.1827–1.5292mBq·m?3，平均值為1.4537mBq·m?3，夏季大氣顆粒物中TSP的210Pb體積比活度平均值為0.3108mBq·m?3；衡陽市秋季的大氣顆粒物中TSP的210Pb體積比活度為3.1558–4.5579 mBq·m?3，春季大氣顆粒物中TSP的210Pb體積比活度為1.6614mBq·m?3。這主要是受溫度、濕度、降雨、風(fēng)速、風(fēng)向等自然因素的影響。

(2) β計(jì)數(shù)法測(cè)量的數(shù)據(jù)的離散性整體比γ計(jì)數(shù)法測(cè)量的數(shù)據(jù)偏大，這表示β計(jì)數(shù)法的穩(wěn)定性較γ計(jì)數(shù)法差，這主要是因?yàn)棣掠?jì)數(shù)法需要進(jìn)行復(fù)雜的化學(xué)前處理，其數(shù)據(jù)精度受人為操作、實(shí)驗(yàn)條件影響。

1 盧正勇. 氣溶膠科學(xué)引論[M]. 北京: 原子能出版社, 2000: 10 LU Zhengyong. Introduction to aerosol science[M]. Beijing: Atomic Energy Press, 2000: 10

2 Porstendorfer J. Physical parameters and dose factors of the radon and thoron decay products[J]. Radiation Protection Dosimetry, 2001, 94(4): 365–373. DOI: 10.1093/oxfordjournals.rpd.a006512

3 Mohamed A, El-Hussein A. Comparison of outdoor activity size distribution of220Rn and222Rn progeny[J]. Applied Radiation and Isotopes, 2005, 62(6): 955–959. DOI: 10.1016/j.apradiso.2005.01.003

4 Baskaran M. Po-210 and Pb-210 as atmospheric tracers and global atmospheric Pb-210 fallout: a review[J]. Journal of Environmental Radioactivity, 2011, 102(5): 500–513. DOI: 10.1016/j.jenvrad.2010.10.007

5 Yue W S, Li X L, Liu J F, et al. Characterization of PM2.5 in the ambient air of Shanghai city by analyzing individual particles[J]. Science of the Total Environment, 2006, 368(2): 916–925. DOI: 10.1016/j.scitotenv.2006. 03.043

6 Baulig A, Sourdeval M, Meyer M, et al. Biological effects of atmospheric particles on human bronchial epithelial cells: comparison with diesel exhausts particles[J]. Toxicology in Vitro, 2003, 17(5–6): 567–573. DOI: 10.1016/S0887-2333(03)00115-2

7 Neas L M. Fine particulate matter and cardiovascular disease[J]. Fuel Processing Technology, 2000, 65: 57–67

8 Du J Z, Huang D K. A comment on “Measurement of the radioactivity of238U?226Ra?210Pb?228Th?232Th?228Ra?137Cs and40K in tea using gamma-spectrometry” by S. Harb[J]. Journal of Radio Analytical and Nuclear Chemistry, 2007, 274(1): 63-66. DOI: 10.1007/s10967-006-6901-5

9 Li X, Zhu J, Guo P, et al. Preliminary studies on the source of PM10 aerosol particles in the atmosphere of Shanghai city by analyzing single aerosol particles[J]. Nuclear Instruments and Methods in Physics Research, 2003, B210: 412–417

10 單健, 肖德濤, 劉衛(wèi), 等. 大氣顆粒物滯留時(shí)間的確定方法及初步實(shí)驗(yàn)驗(yàn)證[J]. 物理學(xué)報(bào), 2010, 59(8): 5358–5363 SHAN Jian, XIAO Detao, LIU Wei, et al. A method of estimating residence time of atmospheric aerosol and its preliminary experimental verification[J]. Acta Physica Sinica, 2010, 59(8): 5358–5363

11 黃展常, 吳其反, 唐志輝, 等. 環(huán)境水中210Pb核素放射性活度的快速分析方法[J]. 核技術(shù), 2013, 36(7): 070301. DOI: 10.11889/j.0253-3219.2013.hjs.36.070301 HUANG Zhanchang, WU Qifan, TANG Zhihui, et al. A method for rapid determination of210Pb radioactivity in environmental water[J]. Nuclear Techniques, 2013, 36(7): 070301. DOI: 10.11889/j.0253-3219.2013.hjs.36.070301

Measurement of the activity-size distribution of210Pb in atmospheric particulates

SHAN Jian1,2LIU Xiaoxiang1XIAO Detao1LIU Wei2LIU Yan1

1(Radon Provincial Key Laboratory, Institute of Nuclear Science and Technology, University of South China, Hengyang 421001, China) 2(Shanghai Institute of Applied Physics, Chinese Academy of Sciences, Jiading Campus, Shanghai 201800, China)

Background:210Pb in air has yielded insight on the sources and transit times of aerosols. In an ideal enclosed air mass (closed system with respect to these nuclides), the residence times of aerosols can be obtained from the activity ratios of210Pb/222Rn,210Bi/210Pb, and210Po/210Pb. Through measuring the activity-size distribution of210Pb in the atmospheric particulates, it is conducive to the use of210Pb tracer study in atmosphere. Purpose: This study aims to determine the activity-size distribution of210Pb in the atmospheric particulates. Methods: Inhalable particles samples at different particle size levels were collected by M235 impact large flow classification sampler. The activity distribution of210Pb in different particle size scales was determined by detecting β rays of210Bi (the first generation progeny of210Pb), and radioactivity of210Pb was determined by using γ spectrometer, respectively. Results: Experimental results show that the distribution range of the specific activity (mBq·m?3) of210Pb in atmospheric particulates is 0.004 4–1.832 1mBq·m?3and the maximum are within 0.49–0.95μm size distribution. Conclusion: The activity-size distribution of210Pb in the atmospheric particulates is mainly due to the effects of temperature, humidity, rainfall, wind speed, wind direction and other natural factors; the discrete of β counting measured data is larger than the discrete of γ counting measured data.

Activity of210Pb, Activity-particle size distribution, Specific activity

TL99

10.11889/j.0253-3219.2015.hjs.38.110303

國家自然科學(xué)基金(No.11375083)資助

單健，男，1975年出生，2010年于南華大學(xué)獲博士學(xué)位，教授，從事氣載放射性測(cè)量及示蹤技術(shù)的研究

2015-09-22，

2015-10-13

CLC TL99