山東某地區(qū)水中鍶-90 放射性水平測量數(shù)據(jù)質(zhì)量的統(tǒng)計學檢驗方法研究

左坤靈 匡開宇 徐永時

（山東省核與輻射安全監(jiān)測中心，山東海陽 265100）

鍶-90 是一種純β 人工放射性核素，半衰期為28.8a，主要來源于核爆炸落下灰、核事故的釋放以及核燃料循環(huán)后段設施運行的釋放[1]。鍶-90 屬于高毒性核素，化學性質(zhì)與鈣相似，是典型的親骨核素。鍶易經(jīng)胃腸道、呼吸道、皮膚和傷口進入人體，并大部分蓄積于骨骼中，引起骨骼組織和造血器官的明顯變化，遠期效應可以致癌。因此，開展水中鍶-90 放射性水平測量及測量數(shù)據(jù)質(zhì)量的統(tǒng)計學檢驗方法研究對于及時發(fā)現(xiàn)和控制水體放射性污染，維護居民用水安全具有重要意義。

1. 樣品與方法

1.1 樣品采集

樣品類型為飲用水，每季度采集一次。采樣器具為聚乙烯塑料桶，使用前先用待采水樣清洗3 次后，再將水樣按采樣量采集于清洗后的樣品容器中，如發(fā)現(xiàn)樣品中有顆粒物或沉淀，應盡快分離。水樣采集后，用濃硝酸酸化到pH=1 ～2。水樣保存期一般不得超過2 個月。

1.2 監(jiān)測方法

樣品均采用《水和生物樣品水中鍶-90 的放射化學分析方法》（HJ815-2016）中快速法進行分析。每個樣品取10L 水樣，核素經(jīng)沉淀濃集后，用酸溶解，調(diào)節(jié)pH 后制備成前處理液，通過涂有HDEHP 的聚三氟氯乙烯（簡稱keL-F）色層柱定量吸附釔，使釔與鍶、銫等低價離子分離。再以1.5mol/L 硝酸淋洗色層柱，清除釔以外的其他被吸附的鈰、钷等稀土離子，并以6mol/L 硝酸解吸釔，以草酸釔沉淀進行β 計數(shù)和稱重[2]。

1.3 儀器與質(zhì)控

測量儀器為美國ORTEX 公司的MPC9604 低本底α/β測量儀。儀器每兩年檢定一次，且都在檢定有效期內(nèi)。

2. 統(tǒng)計分布假設檢驗過程與結(jié)果分析

環(huán)境統(tǒng)計學作為統(tǒng)計學在環(huán)境領(lǐng)域的應用，要求參與統(tǒng)計的測量數(shù)據(jù)能夠代表總體環(huán)境水平，可以通過合理的假設找到一種適合的統(tǒng)計學分布并加以檢驗[3]。

水中鍶-90 測量結(jié)果受多種因素影響，其不確定度主要包括以下4 個方面：儀器探測效率的不確定度；樣品化學回收率的不確定度；樣品測量的不確定度；取樣的不確定度[4]。其中儀器探測效率和樣品測量的不確定度占比最大，約80%以上。儀器探測效率的不確定度因釔-90 探測效率步驟多、存在的不確定度因素多、不確定度分量大，一般變化不大；樣品測量的不確定度與樣品計數(shù)有關(guān)，在同樣測量條件下放射性水平越低，不確定度越大；樣品化學回收率的不確定度，主要與分析方法有關(guān)；取樣的不確定度，與取樣量和取樣方式有關(guān)。

實驗樣品化學回收率均在80%左右，取樣量和取樣方式一致。在做好上述測量過程質(zhì)量控制的基礎(chǔ)上，可暫時忽略不確定度對統(tǒng)計分布假設與檢驗的影響，認為水中鍶-90 測量數(shù)據(jù)具備準確、穩(wěn)定、連續(xù)、受環(huán)境本身非輻射水平因素影響較少，且可以積累足夠數(shù)量的數(shù)據(jù)。所以，其數(shù)據(jù)作為樣本，具有一定的代表性，可視為連續(xù)隨機變量，從而對其進行正態(tài)分布假設，并對該假設進行檢驗。

2.1 數(shù)據(jù)預處理

監(jiān)測數(shù)據(jù)中出現(xiàn)個別缺失值，在環(huán)境監(jiān)測及環(huán)境統(tǒng)計學中經(jīng)常遇到[5]。為更好地說明正態(tài)分布假設及其檢驗方法在本次實驗中的應用，對實驗原始數(shù)據(jù)中一組按照時間排列樣本中的單個缺失值，進行臨近點均值插補（臨近序數(shù)首選2,不滿足2 個臨近值則臨近序數(shù)選1）。由于2018年一季度飲用水1 點位和飲用水3 點位未采集到樣品，進行臨近點均值插補。經(jīng)預處理后的數(shù)據(jù)如表1 所示。

表1 飲用水中鍶-90放射性水平測量結(jié)果（單位：mBq/L）

2.2 正態(tài)分布檢驗方法與結(jié)果

檢驗一組數(shù)據(jù)是否服從正態(tài)分布通常有3 種方法[6]。

（1）偏度系數(shù)和峰度系數(shù)檢驗：該方法為定量檢驗，偏度系數(shù)/偏度系數(shù)標準誤的絕對值小于1.96，則偏度符合正態(tài)分布；峰度同理。

（2）夏皮洛-威爾克（S-W）正態(tài)分布檢驗（單組樣本量小于2000：該方法為定量檢驗，平均值置信區(qū)間取95%，顯著性大于0.05 則認為符合正態(tài)分布。

（3）累積概率（P-P 圖）與分位數(shù)（Q-Q 圖）檢驗：該方法為定性檢驗，P-P 圖是根據(jù)測量數(shù)據(jù)的實際累積概率與其假設理論分部累積概率所繪制的散點圖，當滿足時P-P 圖近似一條直線；Q-Q 圖描述的是分位數(shù)，其他與P-P 圖同理。

從上述檢驗情況來看：飲用水中鍶-90 測量數(shù)據(jù)全部服從正態(tài)分布假設。3 種方法的檢驗結(jié)果均服從正態(tài)分布。即樣本數(shù)據(jù)偏（峰）度系數(shù)/偏（峰）度系數(shù)標準誤的絕對值均小于1.96，夏皮洛-威爾克（S-W）正態(tài)分布檢驗（置信度95%）的顯著性也均大于0.05，其P-P 圖和Q-Q圖也均近似一條直線。詳見表2、圖1 和圖2。

圖2 Q-Q圖正態(tài)分布檢驗

3. 結(jié)語

本文通過對實驗區(qū)域內(nèi)飲用水測量結(jié)果進行正態(tài)分布假設及3 種方法檢驗，結(jié)果顯示飲用水服從正態(tài)分布，因此認為相關(guān)實驗數(shù)據(jù)總體質(zhì)量良好，飲用水中鍶-90 測量數(shù)據(jù)在合理區(qū)間內(nèi)正常波動。同時，本文選用的相關(guān)假設、檢驗可以作為對環(huán)境水中鍶-90 測量數(shù)據(jù)質(zhì)量定量評價，以及水中鍶-90 輻射水平描述的有效手段。