中子伽馬能譜測井資料處理研究

阮凡　楊?！」鶜g歡　戴家才

摘 要：地層元素測井能夠通過中子與地層作用反饋回來的伽馬射線信息得到地層元素的質(zhì)量分?jǐn)?shù)，文章采用蒙特卡羅的方法先模擬制作了地層中Fe、Ca、Si、Al、C、O六種元素的標(biāo)準(zhǔn)非彈譜，然后通過制作已知成分的混合譜求出了地層的歸一化因子F和相對(duì)靈敏度因子Srj，并用這兩個(gè)參數(shù)解譜出來了不同礦物含量地層的元素含量。

關(guān)鍵詞：標(biāo)準(zhǔn)非彈譜；解譜；蒙特卡羅；地層元素

1 中子能譜解譜

1.1 最小二乘法解譜理論

中子與地層中的元素作用主要是非彈性散射反應(yīng)和俘獲反應(yīng)，并釋放出相應(yīng)的伽馬射線，測井中的伽馬射線用儀器探測器分多道記錄，探測器接收到的伽馬射線是地層中每一種元素與中子作用釋放出來的的伽馬射線的線性疊加。例如混合地層由m中元素組成，接收器分n道接收混合伽馬射線，則理論上第i到接收到的伽馬計(jì)數(shù)率Ci為：

2 蒙特卡羅模擬

用MCNP-5建立地層及儀器模型，根據(jù)標(biāo)準(zhǔn)刻度井群建立了標(biāo)準(zhǔn)井模型，井模型的具體參數(shù)如下：高度200cm，直徑212cm，模擬無限大地層，井眼直徑20cm，采用BGO探測器，尺寸為4.5×100mm，探測器貼近井壁，中子源為14Mev各向同性的點(diǎn)源，中子源和探測器中間隔有45×100mm的鉛屏蔽器，防止中子直接打在探測器上，地層填充為分別為Fe、Ca、Si、Al、C、O六種單質(zhì)元素及由不同含量的Fe2O3、CaCO3、Al2O3、SiO2四種氧化物組成的混合地層，用MCNP-5進(jìn)行中子能量0.1Mev截?cái)?，模擬得到元素的標(biāo)準(zhǔn)非彈譜和混合地層的混合譜。能譜接收區(qū)為0.6-8Mev。得到Ca、Si、Fe、Al、O、C的標(biāo)準(zhǔn)非彈能譜如下圖（1-6）所示：

由于能譜的高低能區(qū)的干擾較多出現(xiàn)譜型的畸變，所以在擬合是要進(jìn)行剔除，選取有效能區(qū)為0.8-7 Mev?；旌献V共進(jìn)行了4組成分的模擬

分析試驗(yàn)結(jié)果可以得知計(jì)算誤差的大小和元素特征道所在位置有很大的關(guān)系一般處于能量道兩端或者元素特征道重疊較多的計(jì)算誤差會(huì)比較大，同時(shí)模型的建立與計(jì)數(shù)結(jié)果的處理與最終的結(jié)果也有很大的關(guān)系。地層礦物和地層元素之間有著穩(wěn)定的數(shù)學(xué)關(guān)系，在確定了礦物成分及其化學(xué)式的基礎(chǔ)上可以把地層元素含量轉(zhuǎn)化為地層礦物含量，從而達(dá)到了分析地層組分的目的。

3 結(jié)束語

文章用蒙特卡羅模擬中子伽馬能譜解譜，初步證實(shí)了最小二乘法解譜的可靠性，但是由于模擬探測器計(jì)數(shù)時(shí)高低能量段的計(jì)數(shù)誤差較大造成了解譜誤差的增大，有必要進(jìn)行模型的改進(jìn)，由于文章的模型為理想的井眼模型，沒有地層水、套管、水泥環(huán)、儀器外殼等眾多因素的影響，實(shí)際地層元素測井中井中的干擾較多素較多，將模型應(yīng)用于實(shí)際解譜還有待進(jìn)一步研究和方法的校正。

參考文獻(xiàn)

[1]龐巨豐.地層元素中子俘獲伽馬能譜測井解釋理論和方法[J].測井技術(shù)，1998，22（2）：116-119.

[2]龐巨豐.地層元素測井中中子-伽馬能譜解析理論及方法[J].同位素，2006，19（2）：70-74.

[3].D_D何雄英.中子發(fā)生器n_密度測井的蒙特卡羅模擬研究[J].原子核物理評(píng)論，2013，6（2）：151-154.

[4]王新光.脈沖中子全譜飽和度測井?dāng)?shù)據(jù)處理方法研究[D].中國石油大學(xué)，2010.

[5]趙國志.蒙特卡羅方法在核測井中的研究與應(yīng)用[D].黑龍江大學(xué)，2010.endprint