焦倉文 鄧 明 陸士立 袁 明 馬艷芳

小口徑γ能譜測井儀研制

焦倉文1,2鄧 明1陸士立2袁 明2馬艷芳2

1（中國地質(zhì)大學(xué)（北京） 地球物理與信息技術(shù)學(xué)院 北京 100083）2（核工業(yè)北京地質(zhì)研究院 中核集團(tuán)鈾資源勘查與評價技術(shù)重點(diǎn)實(shí)驗(yàn)室 北京 100029）

為實(shí)現(xiàn)在鈾、釷混合型礦床上準(zhǔn)確測定鈾含量之目的，開發(fā)了直徑45 mm的小口徑自然γ能譜測井儀。該γ能譜儀采用閃爍晶體探測器鍺酸鉍(Bi3Ge4O12, BGO)，晶體為?30 mm×100 mm；為減小死時間，從系統(tǒng)角度提出旨在提高數(shù)據(jù)吞吐率的基于雙端口RAM加雙CPU結(jié)構(gòu)的儀器設(shè)計方案。針對自動穩(wěn)譜難題，軟件上采用基于133Ba自穩(wěn)參考源的實(shí)時四能窗穩(wěn)譜算法；硬件上基于溫度粗調(diào)輔以光電倍增管(Photomultiplier Tube, PMT)高壓調(diào)節(jié)的思想，測井儀開機(jī)階段根據(jù)測井環(huán)境溫度設(shè)定PMT高壓初值，測井過程中根據(jù)譜漂移程度實(shí)時微調(diào)高壓，兩者結(jié)合自動穩(wěn)譜精度有大幅度提高，結(jié)合Pb-Sn組合屏蔽低能γ射線解決了高含量鈾礦體強(qiáng)γ射線對測井儀自動穩(wěn)譜系統(tǒng)干擾的難題。經(jīng)放射性標(biāo)準(zhǔn)模型站對譜儀進(jìn)行多模型標(biāo)定、死時間修正、水和鐵套管吸收系數(shù)修正等檢測，以及對能譜儀野外實(shí)際測試結(jié)果對比，結(jié)果表明，該小口徑自然γ能譜測井儀檢查測井異常相對誤差不大于10%。

小口徑，鍺酸鉍閃爍晶體探測器，自動穩(wěn)譜

自然γ能譜測井是在鉆孔內(nèi)對巖石自然γ放射線進(jìn)行能譜測量與分析，通過分別測定地層內(nèi)鈾、釷、鉀含量來研究井剖面地層性質(zhì)的測井方法。目前，國內(nèi)γ能譜測井儀在石油工業(yè)的應(yīng)用主要是通過自然γ能譜測井所測量的U、Th、K含量來識別巖性、研究沉積環(huán)境、生油層、尋找儲集層、確定粘土含量等目的，特點(diǎn)是含量測量范圍小、測井儀直徑大（如中國石油測井公司SNGR5410、山東勝利石油工程技術(shù)公司SL6329、北京環(huán)鼎科技的SGS-1C自然γ能譜測井儀，直徑一般大于90 mm），不適宜鈾礦勘查中小口徑、高放射性含量的測量要求。國外Mount Sopris（美國）、RG（英國）、Geovista（英國）等公司生產(chǎn)的自然γ能譜測井儀直徑大多在50?60 mm，多數(shù)采用鍺酸鉍(Bi3Ge4O12, BGO)、NaI和BaLr3閃爍晶體探測器。我國一些熱液型鈾礦床屬于鈾、釷混合型礦床，對于此類礦床需要采用γ能譜測井方法分別測定礦石的鈾、釷含量，才能準(zhǔn)確計算鈾的儲量[1]，且鉆孔裸眼直徑都為60 mm，因此提出開發(fā)小口徑γ能譜測井儀。

在γ能譜測井中，目前普遍采用能量分辨率較高、價格便宜的NaI(TI)閃爍晶體作為γ射線全譜測量用探測器。鈾礦勘查中由于受鉆探孔直徑的限制，晶體尺寸較小，隨之帶來的問題是道計數(shù)率偏低造成統(tǒng)計漲落誤差大。經(jīng)過試驗(yàn)比較和實(shí)際測試后，本項(xiàng)目選用探測效率較高的BGO晶體，其與同體積的NaI(Tl)閃爍晶體探測器相比，計數(shù)率可提高2?4倍。

1 儀器設(shè)計方案

為實(shí)現(xiàn)連續(xù)γ能譜測井的目的，小口徑γ能譜測井儀N451（核工業(yè)北京地質(zhì)研究院）在硬件設(shè)計上采用雙微處理器(Central Processing Unit, CPU)加雙端口RAM的架構(gòu)。為提高數(shù)據(jù)吞吐率，采取使采樣保持電路(Sample and Hold, SH)、閾值比較電路(Threshold Compare Circuit)、A/D轉(zhuǎn)換電路、地址鎖存電路(Address Latch)以及采樣保持電容瀉放電電路等構(gòu)成一個硬件閉環(huán)系統(tǒng)的方案，其整機(jī)原理框圖如圖1所示。為保證能譜數(shù)據(jù)幀的完整性，兩微處理器之間設(shè)計有同步機(jī)制，就數(shù)據(jù)幀本身而言，就好似一個先進(jìn)先出(First In First Out, FIFO)存儲器，最先被微處理器1寫進(jìn)去的能譜數(shù)據(jù)幀則最先被微處理器2取走[2]。微處理器1承擔(dān)高速γ脈沖采集，微處理器2實(shí)現(xiàn)數(shù)據(jù)讀取、并道處理以及自動穩(wěn)譜等功能，并實(shí)時將譜數(shù)據(jù)通過異步串行口發(fā)送到地面系統(tǒng)。

硬件電路實(shí)現(xiàn)γ脈沖數(shù)據(jù)吞吐率達(dá)333 k·s?1，達(dá)到減小死時間的目的。為提高能譜儀的微分非線性，采用高位數(shù)AD轉(zhuǎn)換器，通過并道方案進(jìn)行平滑處理。

圖1 小口徑N451 γ能譜測井儀原理框圖Fig.1 Functional block diagram of slim borehole logger N451.

2 自動穩(wěn)譜

γ能譜測井儀自動穩(wěn)譜精度決定了核素含量的解算精度。γ能譜測井不同于車載γ譜儀、航空γ譜儀，后者在自動穩(wěn)譜方面主要針對溫度引起的譜漂移問題[3?4]，采用天然40K的1.460 MeV能峰和208Tl的2.614 MeV能峰等作為基準(zhǔn)進(jìn)行穩(wěn)譜，優(yōu)點(diǎn)是沒有自穩(wěn)源引起的本底干擾。而鉆孔中進(jìn)行γ能譜測井，晶體探測器更接近圍巖礦體且放射性強(qiáng)度變化劇烈，要求自動穩(wěn)譜反應(yīng)迅速，對于薄礦層尤其如此。測井儀譜漂內(nèi)在的影響因素主要是由光電倍增管和晶體探測器引起；光電倍增管的增益隨溫度、高壓的波動以及γ脈沖計數(shù)率而變化；晶體探測器的光輸出脈沖幅度隨環(huán)境溫度升高而降低。針對BGO溫度系數(shù)大的缺點(diǎn)[5]，提出相應(yīng)的自動穩(wěn)譜策略。在傳統(tǒng)的γ能譜儀中，多采用附加參考源的穩(wěn)譜方法[6?11]，根據(jù)特定能量的光電峰對應(yīng)道址變化而實(shí)時加以調(diào)整。γ能譜測井儀N451硬件上采用高精度電壓基準(zhǔn)結(jié)合數(shù)字電位器控制光電倍增管(Photomultiplier Tube, PMT)高壓輸出的方案，軟件上輔以133Ba自穩(wěn)源實(shí)現(xiàn)自動穩(wěn)譜功能，如圖1所示。針對兩能窗穩(wěn)譜精度差的情況[8]，采用基于自穩(wěn)參考源的四能窗實(shí)時穩(wěn)譜軟件設(shè)計方案[9]，這里涉及到參考源能窗寬度選擇(表1)、峰形識別(圖2)以及自穩(wěn)參考指標(biāo)調(diào)整閾值的確定等一系列問題。

N451 γ能譜測井儀內(nèi)置的133Ba自穩(wěn)基準(zhǔn)源有較高的分支比，其γ射線能量主要是356 keV，其他較弱的有302 keV、79 keV等，用作穩(wěn)譜基準(zhǔn)的能峰是356 keV峰。結(jié)合133Ba峰形的特點(diǎn)，基于自穩(wěn)源的356 keV能峰的各能窗計數(shù)率的比例關(guān)系如圖2所示。

表1 133Ba 356 keV能峰的能窗選擇Table 1 Energy windows 133Ba 356 keV signature.

圖2 133Ba在356 keV能峰的能窗記數(shù)比例關(guān)系Fig.2 Count ratio of the energy windows centered on 133Ba at 356 keV.

本底的扣除：定義Ncor=[N1?N4]/3，這里選取常數(shù)3的依據(jù)是，D=[N4u+N4d]/2、C=[N3u+N3d]/2、B=[N2u+N2d]/2、A=[N1u+N1d]/2（腳標(biāo)代表相應(yīng)能窗的上道和下道道址，圖2中未標(biāo)出），根據(jù)幾何關(guān)系A(chǔ)D/BC=2.86，為計算方便近似為3，這里關(guān)系N2?Ncor≈N3應(yīng)成立。定義自動穩(wěn)譜參考指標(biāo)CFF′=[(N2?Ncor) ?2N3]/[(N2?Ncor)+2N3]，驗(yàn)證其有效性是通過設(shè)定不同PMT高壓值后觀察CFF′指標(biāo)的變化而獲知的，如表2所示，在不同PMT高壓條件下，依據(jù)各自穩(wěn)能窗計數(shù)分別按上式計算CFF′值，負(fù)值說明譜向上漂移，正值則說明譜已向低端漂移（初始參考HV=1 110 V），說明其對譜漂是敏感的，可以作為自穩(wěn)參考調(diào)整指標(biāo)。

表2 驗(yàn)證自動穩(wěn)譜參考指標(biāo)CFF′的有效性Table 2 Verification of parameter CFF′ for auto-stabilization through HV.

經(jīng)理論計算結(jié)合工程實(shí)際[10]，確定自穩(wěn)參考調(diào)整閥值CFF'為3%并設(shè)定自穩(wěn)調(diào)整周期為4 s，通過實(shí)際測試，結(jié)果表明自穩(wěn)控制精度有了大幅度的提高。試驗(yàn)用溫控箱為廣東東莞貝爾公司的BE-HL-150H3溫濕度控制儀（溫度?20 oC至150 oC，穩(wěn)定度±0.5 oC），無自穩(wěn)與具備自動穩(wěn)譜功能的試驗(yàn)測試結(jié)果對比如圖3所示。圖3(a)和(b)分別為無自穩(wěn)時和有自穩(wěn)時溫度引起的譜漂。

由圖3(a)可知，無自穩(wěn)時小口徑γ能譜測井儀N451溫度從0 oC升到60 oC，就133Ba的356 keV能峰而言，道漂移25道，而208Tl的2614 keV能峰漂移200道。隨著溫度升高，BGO晶體分辨率有所降低，能峰幅度減小；相對應(yīng)地由圖3(b)可知，自穩(wěn)狀態(tài)下在20?45 oC內(nèi)，譜漂校正效果有明顯的改善，對208Tl的2614 keV能峰的譜漂控制在±2道內(nèi)（就512道而言）。

當(dāng)?shù)V體的鈾含量很高時，礦石的γ射線會嚴(yán)重干擾自動穩(wěn)譜系統(tǒng)，使其無法正常工作。N451測井儀在UF-1.0模型（1%鈾含量）中所測的譜如圖4(a)所示，因模型的低能γ射線很強(qiáng)，133Ba峰將被淹沒，自穩(wěn)受到強(qiáng)烈的干擾。為此，采取了兩項(xiàng)措施：一是增強(qiáng)133Ba源；二是在BGO晶體外加一個Pb-Sn組合屏以吸收低能γ射線。在U-1.0模型譜中其結(jié)果如圖4(b)所示，133Ba峰高高突起在U模型的低能γ射線譜之上，壓低了外部射線對儀器自穩(wěn)的干擾。

需說明的是，實(shí)際應(yīng)用時要綜合考慮測井速度、分辨礦層厚度兩方面因素。通常情況下測井深度采樣間隔5 cm，以每4 s作為自穩(wěn)調(diào)整周期，測井速度3 m·min?1，N451數(shù)據(jù)采集速度為1 spectrum·s?1計，則相當(dāng)于1 point·s?1，即4 point·s?1（4個spectrum） =80 cm·s?1，在上述測井速度條件下，N451能夠基本不失真地采集能譜數(shù)據(jù)所能分辨的礦層厚度為20 cm。

圖3 溫度對測井儀N451能譜漂移的影響(a) 無自穩(wěn)時，(b) 有自穩(wěn)時Fig.3 Spectra drift of N451 due to temperature. (a) Without stabilization, (b) With stabilization

圖4 強(qiáng)低能γ射線對N451 133Ba自穩(wěn)能峰的干擾(a) 無Pb-Sn組合屏，(b) 有Pb-Sn組合屏Fig.4 Influence of low energy γ ray on the 133Ba signature ofN451. (a) Without Pb-Sn filter, (b) With Pb-Sn filter

3 標(biāo)定與野外試驗(yàn)

在核工業(yè)放射性勘查計量站的不同含量U、Th、K基準(zhǔn)模型、混合基準(zhǔn)模型、本底基準(zhǔn)模型上(圖5)，對γ能譜測井儀N451(圖6)進(jìn)行數(shù)個模型井和鐳點(diǎn)源的標(biāo)定；還分別進(jìn)行了鐵、水吸收系數(shù)的測定及修正；隨后在江西相山地區(qū)鈾、釷混合礦床進(jìn)行了七口井的γ能譜測井。由于目前沒有可以直接比對的γ能譜測井儀，γ能譜測井解釋結(jié)果與作為核工業(yè)儲量計算標(biāo)準(zhǔn)的FD-3019 γ總量測井儀測井結(jié)果做了對比，工業(yè)礦(鈾含量＞0.05%) 32.6 m，礦化段(鈾含量0.03%?0.05%) 7.5 m，異常段(鈾含量0.01%?0.03%) 31.4 m，總計米百分含量(m%)的相對誤差為?4.23%。FD-3019相對于中子活化分析結(jié)果的相對誤差為33.65%；N451相對于中子活化分析結(jié)果的相對誤差為26.18%。表3和表4分別是鉆孔ZK80-37的FD-3019和N451 γ能譜測井結(jié)果，對比表明γ能譜測井結(jié)果良好，符合《EJ/T 611-2005 γ測井規(guī)范》中“檢查測井異常面積或米百分值允許相對誤差應(yīng)不大于10%”的規(guī)定。

圖5 標(biāo)定γ能譜測井儀的飽和基準(zhǔn)模型Fig.5 Saturated module body sources for γ spectrometry logger calibration.

圖6 小口徑γ能譜測井儀N451Fig.6 Slim borehole γ spectrometry logger N451.

表3 鉆孔ZK80-37中孔深824?831 m段的分層解釋結(jié)果(FD-3019)Table 3 Layer interpretation result from 824?831 m of ZK80-37 (FD-3019).

表4 鉆孔ZK80-37中記錄孔深824?831 m段的分層解釋結(jié)果表(N451)Table 4 Layer interpretation result from 824?831 m of ZK80-37 (N451).

4 結(jié)語

(1) 在小口徑γ能譜測井晶體探測器尺寸受限的情形下，采用基于雙端口RAM加雙CPU結(jié)構(gòu)的硬件設(shè)計方案，實(shí)現(xiàn)了高采集數(shù)據(jù)吞吐率下減小死時間和降低統(tǒng)計漲落引起測量誤差的目標(biāo)。

(2) 由于鉆孔內(nèi)晶體探測器更接近圍巖礦體，有射線強(qiáng)度變化劇烈的特殊性，自動穩(wěn)譜在硬件上采用高精度電壓基準(zhǔn)加高位數(shù)字電位器調(diào)整PMT高壓輸出的方案，軟件上采用四能窗實(shí)時自穩(wěn)算法，實(shí)現(xiàn)了對N451穩(wěn)譜精度的大幅度提高，同時擴(kuò)大能譜測量線性范圍。

(3) 采用Pb-Sn組合屏，屏蔽低能γ射線對自穩(wěn)峰的干擾實(shí)現(xiàn)從本底模型到1%鈾含量模型的自動穩(wěn)譜。調(diào)整自穩(wěn)源強(qiáng)度以期達(dá)到引起死時間最少、譜本底盡可能低而又不影響自穩(wěn)效果的目的，優(yōu)化基于峰形的穩(wěn)譜算法進(jìn)一步提高了自動穩(wěn)譜精度。致謝 本測井儀研制過程中得到了陸士立研究員的悉心指導(dǎo)及袁明高工的幫助，對在相山野外γ能譜測井中江西核工業(yè)261大隊(duì)銀泳兵工程師等給予的支持與合作表示感謝。

1 湯彬, 吳永鵬, 張雄杰, 等. 高精度能譜測井與234mPa特征γ射線的“直接鈾定量”技術(shù)[J]. 核技術(shù), 2012, 35(10): 745?750 TANG Bin, WU Yongpeng, ZHANG Xiongjie, et al. Direct uranium quantitative detection in high precision γ-logging by 1.001 MeV234mPa γ-ray[J]. Nuclear Techniques, 2012, 35(10): 745–750

2 焦倉文, 陸士立, 鄧明, 等. 基于雙端口RAM的高速γ脈沖幅度分析電路設(shè)計[J]. 核電子學(xué)與探測技術(shù)，2012, 32(6): 706–708 JIAO Cangwen, LU Shili, DENG Ming, et al. The design of high speed gamma ray pulses height analysis circuit[J]. Nuclear Electronics & Detection Technology, 2012, 32(6): 706–708

3 羅耀耀. 航空γ能譜系統(tǒng)中自動穩(wěn)譜的研究與實(shí)現(xiàn)[D].成都理工大學(xué), 2010 LUO Yaoyao. Research and realization on the auto-spectral method for airborne gamma-ray spectrometry system[D]. Chengdu University of Technology, 2010

4 吳永鵬, 賴萬昌, 葛良全, 等. 多道γ能譜儀中的特征峰穩(wěn)譜技術(shù)[J]. 物探與化探, 2003, 27(2): 131–134 WU Yongpeng, LAI Wanchang, GE Liangquan, et al. A new spectrum stabilization technique in portable gamma ray spectrometry[J]. Geophysical & Geochemical Exploration, 2003, 27(2): 131–134

5 曹西征. BGO晶體探頭的溫度特性研究[D]. 東北師范大學(xué), 2004 CAO Xizheng. Research on the temperature characteristics of BGO crystal[D]. Northeast Normal University, 2004

6 谷鳴, 王功慶, 楊錦晴, 等. γ能譜穩(wěn)譜獲取的研究[J].核技術(shù), 1996, 19(12): 712–716 GU Ming, WANG Gongqing, YANG Jinqing, et al. A study on acquiring steady gamma ray spectra[J]. Nuclear Techniques, 1996, 19(12): 712–716

7 敖奇, 魏義祥, 屈建石. NaI譜儀數(shù)字穩(wěn)譜方法設(shè)計[J].核技術(shù)與探測技術(shù), 2009, 29(1): 19–22 AO Qi, WEI Yixiang, QU Jianshi. Design of digital stabilizing method for NaI spectroscopy[J]. Nuclear Electronics & Detection Technology, 2009, 29(1): 19–22

8 鞏建忠, 王衛(wèi)華. 一種用于自然γ能譜測井儀器的自動穩(wěn)譜技術(shù)[J]. 石油儀器, 2006, 20(5): 6–7 GONG Jianzhong, WANG Weihua. A new automatic spectra-stabilization method for natural gamma ray spectra-log[J]. Petroleum Instrument, 2006, 20(5): 6–7

9 魯保平, 張惠芳. 四能窗穩(wěn)譜技術(shù)在巖性密度測井儀中的應(yīng)用[J]. 測井技術(shù), 2008, 32(1): 76–79 LU Baoping, ZHANG Huifang. Spectrum stabilization with 4 energy windows in litho-density logging tool[J]. Well Logging Technology, 2008, 32(1): 76–79

10 李勇平, 方建國. 數(shù)控核測井儀中的實(shí)時穩(wěn)譜程序設(shè)計[J]. 核技術(shù), 1996, 19(3): 155–159 LI Yongping, FANG Jianguo. Real time program design of spectrum stabilization for the computerized nuclear logging instruments[J]. Nuclear Techniques, 1996, 19(3): 155–159

11 柏林, 田彥民, 范廣軍. 基于多道的自然γ能譜測井儀的研制[J]. 艦船科學(xué)技術(shù), 2008, 30(6): 291–293, 311 BAI Lin, TIAN Yanmin, FAN Guangjun. Development of a natural logging gamma ray instrument based on multi-channel[J]. Ship Science and Technology, 2008, 30(6): 291–293, 311

CLC TL817

Research and development of slim borehole gamma ray spectrometry logger

JIAO Cangwen1,2DENG Ming1LU Shili2YUAN Ming2MA Yanfang2
1(School of Geophysics and Information Technology, China University of Geosciences, Beijing 100083, China) 2(China National Nuclear Corporation, Key Laboratory of Uranium Resources Exploration and Evaluation Technology, Beijing Research Institute of Uranium Geology, Beijing 100029, China)

Background: γ spectrometry logging is one method for accurately computing uranium content in uranium-thorium mixed type deposit for slim borehole application. Purpose: The aim is to design slim borehole γ spectrometry logger N451 with ?30 mm×100 mm high detection efficiency Bi4Ge3O12(BGO) scintillation detector, whose diameter and length are 45 mm and 1 300 mm, respectively. Methods: By applying the configuration of dual-microprocessor plus dual port RAM to improving the data throughput rate, the aim of reducing dead time was achieved. As to auto-stabilization, by means of 4-window algorithm with133Ba reference source, plus temperature-based coarse set while power on and fine tune Photomultiplier Tube (PMT) high voltage while logging, the auto-stabilization accuracy was considerably improved. Using Pb-Sn mask method, N451 overcomes the interference on auto-stabilization system from the high level uranium deposit low energy γ-ray. Results: After multiple saturated modules and226Ra point source calibration and Fe as well as H2O coefficient correction, field logging was surveyed. Conclusion: The uranium content relative error was less than 10%. The slim borehole γ spectrometry logger is up to the specification.

Slim borehole, Bi3Ge4O12(BGO) crystal, Automatic spectrum stabilization

TL817

10.11889/j.0253-3219.2015.hjs.38.040403

國防預(yù)研基金(No.41603070203)資助

焦倉文，男，1969年出生，1992年畢業(yè)于長春地質(zhì)學(xué)院，高級工程師，現(xiàn)為在職博士生，研究領(lǐng)域?yàn)榈厍蛭锢頊y井

鄧明，E-mail: dengming@cugb.edu.cn

2014-11-13，

2014-12-15