唐曉川 管少斌 胡明考 李峰林 杜曉立 歐陽游 劉珊珊

摘要：介紹了γ測井儀的測量依據(jù)、測量方法、測量標準裝置各項指標和量值傳遞的數(shù)學(xué)模型，討論了γ測井儀測量結(jié)果的4個影響因素及不確定度來源，給出了相應(yīng)的分析評定方法；通過FD-3019型γ測井儀的應(yīng)用實例，評定了FD-3019型γ測井儀測量結(jié)果的擴展不確定度不超過10%，其最大項分量為模型體源標準量值的不確定度，該項標準不確定度為4.0%，統(tǒng)計并分析了11年來γ測井儀鈾換算系數(shù)的不確定度，2007年新檢定規(guī)程發(fā)布之前最大不超過11%，之后最大不超過10%。

關(guān)鍵詞：γ測井儀；不確定度；測井模型標準裝置

中圖分類號 ： TL77文獻標識碼： A

γ測井儀在鈾礦勘查以及鈾礦地浸采冶領(lǐng)域中發(fā)揮著重要的作用，其主要功能是利用檢定得出的鈾換算系數(shù)將所測計數(shù)率換算為鈾含量，目前國內(nèi)常用的γ測井儀主要有FD-3019、HD-4002等[1]。作為γ測井規(guī)范中鈾礦儲量計算的唯一指定儀器，γ測井儀的測量精度問題則直接關(guān)系到后期儲量估算、采冶生產(chǎn)等各個方面，而現(xiàn)階段衡量測量精度的關(guān)鍵指標是測量不確定度，因此評定γ測井儀鈾換算系數(shù)的測量不確定度也愈顯重要。

1 測量依據(jù)、方法及數(shù)學(xué)模型

1.1 測量依據(jù)

本文依據(jù)標準為EJ/J 1-2007《γ測井儀檢定規(guī)程》，測量所用計量標準裝置為測井模型標準裝置和固體鐳標準源。測井模型標準裝置位于石家莊市核工業(yè)航測遙感中心，為國家核行業(yè)最高標準，證書號為[2001]國防計標字第1462號，[2005]國量標核證字第004號，其技術(shù)指標見表1；鐳標準源為英制固體鐳源，其技術(shù)指標見表2[2]。

表1 γ測井模型標準裝置技術(shù)指標

序號 類型 主元素含量 空氣比釋動能率

范圍 標準不確定度 范圍/nGyh-1 標準不確定度

1 本底模型 — — — —

2 鈾模型 eU：≥50×10-6 ≤6% ≥120 ≤6%

3 釷模型 eTh：≥50×10-6 ≥60

4 鉀模型 K：≥3% ≥30

5 混合模型 eU：≥50×10-6

eTh：≥50×10-6 ≥85

表2 固體鐳標準源指標

鐳源號 主要技術(shù)指標 數(shù)量

活度 不確定度 %（k=3）

6-1038 0.95 mg 3.1% 1

1.2 被測對象

FD-3019型測井儀（以下簡稱儀器）。

1.3 環(huán)境條件

溫度：5℃～35℃；濕度：<80%。

1.4 測量方法

將探測器置于鈾模型體源井孔內(nèi)，使探測器靈敏體積的有效幾何中心與井孔活性區(qū)的有效幾何中心重合后進行測量，得到儀器鈾含量靈敏度，再用固體鐳標準源掛在不同距離產(chǎn)生的空氣比釋動能率，測量得到儀器空氣比釋動能率靈敏度[3-4]，按公式（1）計算得到儀器鈾換算系數(shù)。

1.5 數(shù)學(xué)模型

（1）

式中：

——鈾換算系數(shù)，（nGy·h-1）/（0.01%eU），eU為當量鈾含量；

——鈾含量靈敏度，（0.01%eUs）-1；

——空氣比釋動能率靈敏度，s-1/（nGyh-1）。

2 不確定度評定

2.1 不確定度來源分析

由上述數(shù)學(xué)模型和測井測量原理可以看出影響鈾換算系數(shù)的不確定度來源有[5]：

1）測井模型標準裝置量值不確定度；

2）固體鐳標準源量值不確定度；

3）儀器檢定過程中的測量不確定度（包括放射性統(tǒng)計漲落、儀器的角響應(yīng)、儀器能量響應(yīng)非線性、環(huán)境輻射貢獻、固體鐳標準源放置的位置）。

2.2 測井模型標準裝置量值的不確定度（ucul）

用B類標準不確定度評定。

該標準裝置為模型體源的一級標準，其K、U、Th、Ra的含量（標準的量值）是與國家一級標準物質(zhì)進行比較測量，通過化學(xué)分析方法所得含量經(jīng)U-Ra平衡修正、氡泄漏修正、濕度修正后的含量。其復(fù)現(xiàn)量值的不確定度由標準裝置證書給出，ucul=4.0%。

2.3 固體鐳標準源量值的不確定度（ ）

用B類標準不確定度評定。

該標準源為6-1038，由國防科技工業(yè)電離輻射一級計量站檢定，其復(fù)現(xiàn)量值的擴展不確定度為3.1%，包含因子k=3，取固體鐳標準源量值的不確定度 =1.1%。

2.4 儀器檢定過程的測量不確定度

儀器在檢定過程中的不確定度主要來源于放射性統(tǒng)計漲落，儀器的角響應(yīng)、環(huán)境輻射貢獻的變化、固體鐳標準源的量值、固體鐳標準源放置的位置等，且大多數(shù)測量過程中有許多變動情況，如天與天的變動性，月與月的變動性，儀器與儀器的變動性等，說明該隨機過程誤差不僅是儀器或測量標準的特性，而且是整個測量過程的特征。

2.4.1 儀器檢定過程的A類不確定度（umA）

儀器檢定過程的A類不確定度來源只有放射性統(tǒng)計漲落引起的不確定度。

放射性統(tǒng)計漲落引起的不確定度大小與測得的計數(shù)相關(guān)，因而可通過控制測量計數(shù)來控制此項不確定度。實際測量次數(shù)為10次，原始數(shù)據(jù)見表3。

從表3中的統(tǒng)計結(jié)果可以看出，除本底模型體源（F-0-I）和鉀模型體源（KF-6-I）外，各模型體源測量結(jié)果的不確定度最大為0.5%，由于本底模型體源和鉀模型體源受綜合環(huán)境輻射影響較大，計算得到的不確定度相對較大，該項不確定度已在環(huán)境輻射變化中予以考慮，在此不做重復(fù)計算。因此，放射性統(tǒng)計漲落引起的不確定度為0.5%。

表3 儀器測量原始數(shù)據(jù)記錄表

模型號 儀器讀數(shù)平均值

/s-1 相對標準偏差

/% 相對不確定度

/%

F-0-I 11.9 5.9 1.9

KF-6-I 39.5 2.3 0.7

UF-0.03-I 800.4 0.7 0.2

UF-0.2-I 5537.4 0.7 0.2

UF-0.5-I 14574.1 0.8 0.3

UF-1.0-I 28055.9 0.7 0.2

ThF-0.05-I 519.0 1.1 0.4

UThF-0.01-0.03-I 575.4 1.5 0.5

2.4.2 儀器檢定測量過程的B類不確定度（umB）

（1）儀器的角響應(yīng)引入的不確定度ucu2

對于FD-3019型測井儀，由于探測器四周管壁以及屏蔽材料的不均勻，導(dǎo)致儀器測量結(jié)果發(fā)生變化，按照檢定規(guī)程的技術(shù)要求儀器角響應(yīng)不大于1%，則在儀器性能正常時，此項因素所引起的不確定度不超過1.0%。

（2）環(huán)境輻射變化引入的不確定度ucu3

該系列模型體源設(shè)置的場地受天氣因素影響很小，但由于檢定時排風程度不可能一致，環(huán)境氡的影響會有所不同，綜合考慮，環(huán)境因素帶來的不確定度仍定為ucu3=0.5%。

（3）固體鐳標準源放置位置的準確性都屬于隨機性的，檢定時對多個源位置測量，用線性擬合求得有關(guān)參數(shù)，因而該因素的影響將顯著減小，引入的不確定度 一般不會超過1.0%。

（4）儀器能量響應(yīng)非線性引起的不確定度ucu4

對于在本系列模型體源特別是高含量模型體源進行檢定的儀器，應(yīng)考慮儀器死時間、脈沖疊加效應(yīng)等因素導(dǎo)致的儀器響應(yīng)非線性的影響， FD-3019型儀器具有死時間修正，此項因素所引起的不確定度不超過1.0%。

上述各分量彼此獨立，因而儀器的檢定過程B類不確定度umB為：

2.5 合成標準不確定度

2.5.1 主要標準不確定度匯總表

根據(jù)以上分析，主要標準不確定度匯總表見表4。

表4 主要標準不確定度匯總表

不確定度來源 標準不確定度 評定方法

模型體源標準量值的不確定度

4.0% B

固體鐳標準源量值的不確定度

1.1% B

儀器檢定過程的A類不確定度

0.5% A

儀器檢定過程的B類不確定度

1.8% B

2.5.2 合成標準不確定度

上述各不確定度來源之間彼此相互獨立。合成標準不確定度 為：

2.5.3 擴展不確定度

主要標準不確定度均為正態(tài)分布，因此P=0.95時，可取包含因子k=2，則：

取測井儀鈾換算系數(shù)的擴展不確定度Urel為10%，k=2。

2.6 統(tǒng)計比較

筆者通過統(tǒng)計11年來我站所檢定的，來自全國2000多臺/次的各型γ測井儀的檢定數(shù)據(jù)發(fā)現(xiàn)，鈾換算系數(shù)的擴展不確定度均≤11%，其原因是2007年之前依照舊規(guī)程JJG（核工）023-1991《測井儀檢定規(guī)程》檢定的，其向上取整后的擴展不確定度為11%，2007年后采用新規(guī)程EJ/J 1-2007《γ測井儀檢定規(guī)程》的，則為10%。

3 結(jié)論

本文分析了γ測井儀鈾換算系數(shù)的不確定度來源，以FD-3019型γ測井儀為例使用 A類和B類不確定度評定方法對各影響因素的不確定度進行了評定，得出其擴展不確定度最大為10%，并通過統(tǒng)計歷年檢定數(shù)據(jù)證實，依據(jù)EJ/J 1-2007《γ測井儀檢定規(guī)程》確定的γ測井儀不確定度最大為10%。

參考文獻：

[1] EJ/T 611-2005，測井規(guī)范[S].北京：國防科學(xué)技術(shù)工業(yè)委員會，2005.

[2] EJ/J 1—2007，測井儀檢定規(guī)程[S]. 北京：國防科學(xué)技術(shù)工業(yè)委員會，2007.

[3] 田德祥.電離輻射量及其單位[M].北京：原子能出版社，1993.

[4] JJF 1035-2006，電離輻射計量名詞及定義[S].北京：國家質(zhì)量監(jiān)督檢驗檢疫總局，2006.

[5] JJF 1059.1-2012，測量不確定度評定與表示[S]. 北京：國家質(zhì)量監(jiān)督檢驗檢疫總局，2012