王宇健 陳志祥

（1.山東省煙臺生態(tài)環(huán)境監(jiān)測中心，山東 煙臺 264001；2.山東省核與輻射安全監(jiān)測中心，山東 濟南 250000）

0 引言

環(huán)境γ 劑量率監(jiān)測在輻射環(huán)境監(jiān)測領(lǐng)域有重要意義，可用于積累天然放射性本底水平數(shù)據(jù)，為異常情況和輻射應(yīng)急決策提供參考；監(jiān)測輻射源放射性水平、了解γ輻射變化趨勢，評估輻射源的影響。瞬時γ 劑量率監(jiān)測、累積γ 劑量率監(jiān)測是環(huán)境γ 劑量率監(jiān)測常用的監(jiān)測方式，兩種方法各有特點，互為補充。

目前，對于輻射環(huán)境γ 劑量率測量結(jié)果準(zhǔn)確性的保證多側(cè)重于測量過程的質(zhì)量控制方面，包括規(guī)范儀器操作使用方法、按時檢定、繪制質(zhì)控圖、儀器比對等方法以保證儀器測量穩(wěn)定性、準(zhǔn)確性，以上方法可保證測量結(jié)果的可信度。對于測得的數(shù)據(jù)本身是否在合理水平缺少必要的甄別手段，常通過簡單的數(shù)據(jù)對比，利用技術(shù)人員的經(jīng)驗判定，這種基于簡單對比和經(jīng)驗判斷的評價方式，說服力相對有限。也有人采用3σ 判定法則進行異常值判定，該方法使用較為簡單，但檢驗效果較差。所以，找到一種基于統(tǒng)計學(xué)的，通過定量計算對γ 劑量率（瞬時、TLD）測量結(jié)果進行判定的合理檢驗方法，有助于提升數(shù)據(jù)質(zhì)量和輻射環(huán)境水平描述的合理性、說服性。

1 統(tǒng)計分布假設(shè)與檢驗方法的選擇

1.1 統(tǒng)計分布假設(shè)的選擇

環(huán)境統(tǒng)計學(xué)作為統(tǒng)計學(xué)在環(huán)境領(lǐng)域的應(yīng)用，要求參與統(tǒng)計的測量數(shù)據(jù)能夠代表總體環(huán)境水平，可以通過合理的假設(shè)找到一種適合的統(tǒng)計學(xué)分布并加以檢驗。按照以上原則，該文對試驗條件設(shè)置如下：1）參考《環(huán)境空氣質(zhì)量監(jiān)測點位布設(shè)技術(shù)規(guī)范（試行）》（HJ 664-2013）和TOPSIS方法，在半徑約10km 范圍內(nèi)各布設(shè)6 個γ 劑量率瞬時和TLD 測量點位，且實驗時間段內(nèi)點位固定不變，以保證區(qū)域覆蓋性和代表性。2）測量全過程、各要素嚴(yán)格執(zhí)行相關(guān)標(biāo)準(zhǔn)及質(zhì)控要求，確保數(shù)據(jù)真實、可靠，在理論上保證測量數(shù)據(jù)—檢驗結(jié)果—輻射環(huán)境水平三者的一致性。3）有一定的測量頻次，且測量時間相對固定，在一定時間段內(nèi)可積累足夠數(shù)量并相對連續(xù)的數(shù)據(jù)以供統(tǒng)計分析，保證統(tǒng)計分布假設(shè)與檢驗的準(zhǔn)確性；四是γ 劑量率瞬時和TLD 測量同點位、同頻次進行，且數(shù)據(jù)經(jīng)過計算轉(zhuǎn)換均以nGy/h 為單位，作為相對獨立，但理論上又能反映同一對象及其變化趨勢的兩種樣本數(shù)據(jù)，其樣本間比較結(jié)果具有一定的說明性和驗證性；五是本底情況下，環(huán)境γ 劑量率瞬時與TLD 測量結(jié)果一般相對穩(wěn)定，出現(xiàn)離群值或整體水平較大波動的原因較易辨識。

同時，作為環(huán)境級測量數(shù)據(jù)，γ 劑量率測量受多種因素影響，在測量過程中也應(yīng)注意和盡量避免以下問題對測量結(jié)果及統(tǒng)計分布的影響：1）TLD 測量過程中重復(fù)測量引入的A 類不確定度以及刻度因子引入的B 類不確定度，以及瞬時測量時，遇到施工、修路等本底環(huán)境臨時或永久改變等；2）有研究表明，TLD 測量結(jié)果在一定程度上受環(huán)境溫、濕度影響，具體影響情況為，由于低溫峰和主計量峰部分重疊，熱釋光測量儀在計數(shù)前的“退火”階段，難以完全排除低溫峰對測量計數(shù)的影響，從而難以排除布片過程中環(huán)境溫度對最終測量結(jié)果的影響，另外長時間連續(xù)降水導(dǎo)致的環(huán)境濕度過高對劑量片吸收情況也有影響；3）一般情況下，瞬時γ 劑量率的測量結(jié)果減去了宇宙射線響應(yīng)值，而TLD 難以準(zhǔn)確測量宇宙射線響應(yīng)值而未予去除，因此兩種測量結(jié)果存在一定差值（正常情況下瞬時γ 劑量率小于TLD 測量數(shù)據(jù)），且差值難以準(zhǔn)確計量，一定程度上影響二者比對；4）TLD 測量數(shù)據(jù)反映一個時段γ 劑量率的平均值，瞬時測量數(shù)據(jù)反映測量當(dāng)時的γ 劑量率水平，二者趨勢一致，但單點位、單次測量值對比結(jié)果僅供參考。

綜上所述，在做好測量過程質(zhì)量控制的基礎(chǔ)上，可暫時忽略上述不利因素對統(tǒng)計分布假設(shè)與檢驗的影響，而認(rèn)為輻射環(huán)境γ 劑量率瞬時和TLD 測量數(shù)據(jù)具備準(zhǔn)確、穩(wěn)定、連續(xù)、受環(huán)境本身非輻射水平因素影響較少，且可以積累足夠數(shù)量的數(shù)據(jù)。所以，其數(shù)據(jù)作為樣本，具有一定的代表性，可視為連續(xù)隨機變量，從而對其進行正態(tài)分布假設(shè)，并對該假設(shè)進行檢驗。

1.2 檢驗方法的選擇

正態(tài)分布是統(tǒng)計分析中最重要的分布，檢驗一組數(shù)據(jù)是否服從正態(tài)分布通常有3 種方法：計算偏度系數(shù)和峰度系數(shù)檢驗；夏皮洛-威爾克（S-W）正態(tài)分布檢驗（樣本容量較小時）；P-P 圖及Q-Q 圖圖形檢驗。該文選擇前兩種定量分析方法進行數(shù)據(jù)分析，兩種方法的檢驗結(jié)果在表達正態(tài)分布不同特性方面各有側(cè)重，為更全面分析試驗數(shù)據(jù)的特性，該文同時選用2 種方法進行檢驗。

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

該文結(jié)合試驗數(shù)據(jù)（瞬時、TLD 各6 組共址測量數(shù)據(jù)，實驗時間2017 年三季度～2021 年二季度，測量頻次一季度/次，不考慮缺失值情況下單組樣本數(shù)量16 個），重點研究統(tǒng)計檢驗方法在γ 劑量率瞬時和TLD 測量結(jié)果分析中的應(yīng)用。監(jiān)測數(shù)據(jù)中出現(xiàn)個別缺失值，是環(huán)境監(jiān)測及環(huán)境統(tǒng)計學(xué)中經(jīng)常遇到的情況。尤其在TLD 測量中，劑量片經(jīng)常會因為其布設(shè)樹木被修剪砍伐而導(dǎo)致數(shù)據(jù)缺失。

在某些環(huán)境水平分析應(yīng)用中，為更好體現(xiàn)整體環(huán)境的連續(xù)水平，且考慮測量數(shù)據(jù)具有序列特征，常對原始測量數(shù)據(jù)中一組按照時間排列樣本中的缺失值進行統(tǒng)計學(xué)補全：對于單個缺失值，可使用臨近點均值插補，臨近值個數(shù)一般選1～3；若測量數(shù)據(jù)中出現(xiàn)連續(xù)2 個缺失值，可進行線性插值插補；若單組樣本首、末位序數(shù)樣本出現(xiàn)缺失值，可進行臨近點線性趨勢插補。

需要注意的是，統(tǒng)計學(xué)補全的數(shù)據(jù)相當(dāng)于人為對該組數(shù)據(jù)的分布特性加入了一定程度的干擾，會讓該組數(shù)據(jù)更加符合某種統(tǒng)計學(xué)分布，其檢驗結(jié)果更傾向于統(tǒng)計學(xué)預(yù)期。為更直接地說明正態(tài)分布假設(shè)及其檢驗方法在該試驗中的應(yīng)用，所以不對缺失值進行補全處理。原始數(shù)據(jù)見表1。

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

使用以下2 種方法對表1 中γ 劑量率瞬時與TLD 測量數(shù)據(jù)的正態(tài)分布假設(shè)進行檢驗。

計算偏度和峰度系數(shù)定量檢驗。偏度是統(tǒng)計數(shù)據(jù)分布偏斜方向和程度的度量，定義為樣本的三階標(biāo)準(zhǔn)化距，樣本偏度系數(shù)如公式（1）所示。

峰度是研究數(shù)據(jù)分布陡峭或平滑的統(tǒng)計量，定義為四階中心距除以方差的平方減3，樣本峰度系數(shù)如公式（2）所示。

以Z 評分判定其是否服從正態(tài)分布，Z=偏度系數(shù)or 峰度系數(shù)/偏度系數(shù)標(biāo)準(zhǔn)誤，在α=0.05 的檢驗水平下，如果Z在±1.96 之間，則可以說明數(shù)據(jù)服從正態(tài)分布。

因公式中各參量可直接讀取或通過簡單計算獲取，不再贅述計算過程，計算結(jié)果見表2。

夏皮洛-威爾克（S-W）正態(tài)分布檢驗。S-W 屬于非參數(shù)檢驗，用來檢驗數(shù)據(jù)是否符合正態(tài)分布，樣本容量較?。颖救萘?-50）時適用。具體檢驗步驟如公式（3）所示。

該文中樣本數(shù)據(jù)假設(shè)檢驗顯著性水平α 取0.05，選取點位1 瞬時γ 劑量率為例進行假設(shè)檢驗。

作出假設(shè)H：點位1 瞬時γ 劑量率數(shù)據(jù)服從正態(tài)分布，備擇假設(shè)為H。

a：將歷次測量瞬時γ 劑量率數(shù)據(jù)從小到大排列成：72.4，73.1，74.4，74.5，77.5，79.5，79.9，80.7，82.7，83.1，85.6，86.0，86.9，93.1，93.6，97.5。

b：查閱GB4882 中a值表a～a分別為0.5056，0.3290，0.2521，0.1939，0.1447，0.1005，0.0593，0.0196，同時計算樣本均值，帶入公式，求得W=0.946。

c：對α=0.05，查表得n=16 時，Z=0.887。

d：由于0.946 ＞0.887，所以不拒絕正態(tài)性原假設(shè)H0。

所有點位瞬時及累積γ 劑量率S-W 檢驗結(jié)果見表2。

正態(tài)分布假設(shè)的檢驗結(jié)果。從上述兩種正態(tài)分布定量檢驗方法對γ 劑量率瞬時和TLD 測量各6 組樣本數(shù)據(jù)的檢驗情況來看，兩種方法的檢驗結(jié)果一致度較高，即：瞬時測量全部6 組樣本數(shù)據(jù)及TLD 測量點位1、2、3、6 四組樣本數(shù)據(jù)服從正態(tài)分布假設(shè)，樣本數(shù)據(jù)其偏（峰）度系數(shù)/偏（峰）度系數(shù)標(biāo)準(zhǔn)誤的絕對值均小于1.96 的同時，其夏皮洛-威爾克（S-W）檢驗結(jié)果也服從正態(tài)分布；TLD 測量點位4、5兩組樣本數(shù)據(jù)不服從正態(tài)分布假設(shè)，其中點位4、5 偏度系數(shù)/偏度系數(shù)標(biāo)準(zhǔn)誤的絕對值均大于1.96，同時點位5 夏皮洛-威爾克（S-W）檢驗W 值小于臨界值。

另外不難看出，雖然TLD 測量點位1、2、3、6 四組數(shù)據(jù)樣本總體服從正態(tài)分布，但也存在一定的“離散”傾向，主要體現(xiàn)在其偏度系數(shù)/偏度系數(shù)標(biāo)準(zhǔn)誤的絕對值比較接近1.96、夏皮洛-威爾克（S-W）檢驗結(jié)果比較接近臨界值。對于不符合正態(tài)分布和有一定程度“離散”傾向的樣本數(shù)據(jù)組，均應(yīng)懷疑其存在離群值。

2.2 數(shù)據(jù)分析

對TLD 數(shù)據(jù)出現(xiàn)的不符合正態(tài)分布檢驗的結(jié)果進行追溯和分析發(fā)現(xiàn)，TLD 測量過程中TLD 測量設(shè)備進行過多次校準(zhǔn)，以2018 年1 季度為界，在此之前，各點位TLD 數(shù)據(jù)整體高于后續(xù)組內(nèi)其他數(shù)據(jù)。經(jīng)分析TLD 測量設(shè)備經(jīng)校準(zhǔn)后刻度因子較高，導(dǎo)致了測量數(shù)據(jù)整體正態(tài)性變差，這也是TLD 各點位出現(xiàn)檢驗結(jié)果不符合正態(tài)分布或接近臨界值的根本原因，2018 年1 季度重新檢定并使用新的刻度因子后，測量數(shù)據(jù)整體穩(wěn)定。

表1 γ 劑量率瞬時與TLD 測量數(shù)據(jù)（單位：nGy/h）

表2 正態(tài)分布假設(shè)檢驗結(jié)果

雖然γ 劑量率瞬時與TLD 測量有諸多不同之處，外界環(huán)境影響以及儀器設(shè)備本身穩(wěn)定性差異必然導(dǎo)致最終測量結(jié)果的不同，但二者在同點位較長時間段內(nèi)測量數(shù)據(jù)的變化趨勢總體保持一致。

3 結(jié)論

該文對試驗區(qū)域內(nèi)（同點位、同頻次）環(huán)境γ 輻射空氣吸收劑量率瞬時和TLD 測量結(jié)果進行了正態(tài)分布假設(shè)與兩種不同方法檢驗，檢驗結(jié)果表明，該試驗數(shù)據(jù)總體質(zhì)量良好，其中瞬時測量數(shù)據(jù)的穩(wěn)定性優(yōu)于TLD 測量數(shù)據(jù)，試驗區(qū)域內(nèi)γ 劑量率輻射水平發(fā)展趨勢在合理區(qū)間內(nèi)正常波動。另外，可通過繪制瞬時與累積測量數(shù)據(jù)趨勢圖的方式，進一步驗證二者測量結(jié)果的合理性。

該文選用的相關(guān)假設(shè)、檢驗及相應(yīng)的數(shù)據(jù)追溯方法，可以作為對環(huán)境γ 劑量率瞬時、TLD 測量數(shù)據(jù)質(zhì)量定量評價，以及環(huán)境γ 劑量率輻射水平描述的有效手段。通過對數(shù)據(jù)本身的統(tǒng)計假設(shè)檢驗，也對儀器設(shè)備的穩(wěn)定性，運行狀況提供了參考資料。

兩種不同的假設(shè)檢驗方法不僅適用于環(huán)境γ 劑量率數(shù)據(jù)分析判定，也可在輻射環(huán)境監(jiān)測的其他項目得以應(yīng)用，特別是正常情況下監(jiān)測結(jié)果穩(wěn)定在探測下限（MDC）以上的環(huán)境樣品，例如核電基地周邊環(huán)境介質(zhì)中總放射性分析、環(huán)境介質(zhì)中氚、碳-14 的分析等。隨著核電廠放射性流出物的排放，采用兩種不同的假設(shè)檢驗方法，可有效甄別一組數(shù)據(jù)是否存在異常，判別核電廠對周邊輻射環(huán)境帶來的影響，為輻射環(huán)境相關(guān)決策提供資料。