張選同
(新疆大學地質與礦業工程學院,新疆烏魯木齊830046)
地面伽馬能譜測量理論基礎的綜述
張選同*
(新疆大學地質與礦業工程學院,新疆烏魯木齊830046)
伽馬能譜測量是地球物理勘探的方法之一。從伽馬射線的發現,伽馬射線與物質的3種相互作用形式,放射性元素在巖石或土壤中的分布規律,以及能譜測量影響的因素等方面進行基礎論述,從而總結伽馬能譜測量在地質勘探中的應用原理。
伽馬能譜;伽馬射線;相互作用;放射性元素分布規律;影響因素
近年來,隨著國民經濟的發展和國防事業的增強,對礦產資源的需求不斷增長。然而礦產資源的稀缺性、礦床區域地質、經濟條件的落后以及交通不便等因素造成了礦產資源勘探的重大難題。伽馬能譜測量是利用各種巖石和土壤的放射性元素含量差異來尋找礦產或查明地下地質構造的一種勘探方法。伽馬能譜勘探具有高效、輕便、預測精度高等優點。因此,在地質勘查中獲得廣泛應用。最初是在放射性地質勘查領域使用,后來在非放射性礦產和地質填圖中使用,近年來開始用于環境監測和保護。
伽馬射線是1900年由法國化學家和物理學家P.V.維拉德發現的,到1903年被物理學家盧瑟福命名[1]。它具有很強的穿透力,其可以穿透幾十厘米厚的巖石和數百米的空氣[2]。伽馬射線主要由以下過程產生[3]:①α或β衰變;②核反應;③由基態的直接激發;④來自較高激發態的γ躍遷。
盡管自然界有很多天然放射性同位素,其中絕大多數是一次核衰變后形成穩定核素,如87Rb(銣)、40K(鉀)等;有些核素按一定順序相繼衰變構成放射性系,天然放射性系列主要是鈾系、釷系和錒鈾系[4-5]。天然放射性核衰變主要有α衰變和β衰變[6]。在α或β衰變之后,形成的子核往往處于激發態,當原子核從激發態向較低能級躍遷時,同時釋放出伽馬射線,這種現象稱為伽馬衰變[7]。伽馬測量的對象主要是衰變的伽馬射線。
伽馬射線與物質的相互作用是伽馬能譜測量的物理基礎。伽馬射線是一種波長極短的電磁波,與物質相互作用和帶電粒子不同,是一種不帶電荷的光子。其與物質相互作用除了具有穿透外,還有吸收和散射。伽馬射線與物質的相互作用有3種形式:光電效應、康普頓效應和電子對效應[4-6,8-9]。
2.1光電效應
伽馬光子與原子殼層軌道上的內層電子碰撞時,將全部的能量傳遞給原子中某個束縛的電子,使電子從原子中發射出來,發射的電子稱為光電子,伽馬光子消失,該過程稱為光電效應。產生光電效應后的原子處于激發態,當外層電子充填內層電子空缺時,多余的能量以輻射特征X射線或俄歇電子等形式放出,原子回到正常狀態。
光電子的飛行方向,一般與伽馬光子的運動方向和能量(hν)有關。根據動量和能量守恒原理,在光子的入射方向和180°方向上是沒有光電子發射。光電效應是吸收低能量過程,當伽馬光子能量較低的時候,光電子垂直于伽馬光子入射的方向飛出概率較大;隨著入射伽馬光子能量不斷的增加,射出的光電子逐漸靠近伽馬光子入射方向。
2.2康普頓效應
伽馬光子與物質相互作用時除了產生光電效應之外,還有康普頓效應。康普頓效應產生兩種散射方式:其一是彈性散射,另一種是彈性碰撞。彈性碰撞是伽馬光子將部分能量傳給電子,電子與入射光子成φ角度射出,稱為反沖電子,被減弱的伽馬光子也改變了方向,與入射光子成θ角度散射出去,該光子稱為散射光子,這種作用稱為康普頓效應。康普頓效應主要吸收中等能量的伽馬光子。天然放射性核素放出的伽馬光子與巖石作用的主要形式是康普頓效應。
2.3電子對效應
隨著入射光子的能量增高,光電效應和康普頓效應逐漸減弱。當入射光子的能量大于1.02 MeV時,伽馬光子和物質產生電子對效應。電子對效應是指入射的伽馬光子在核庫侖場作用下能量完全被吸收,電子從負能級躍遷到正能級,同時產生一對粒子(正、負電子),且隨著伽馬光子能量的增加逐漸成為主要的作用方式。
從圖1可以看出:隨著入射伽馬光子能量的變化,3種效應所占比例是不同的[8]。一般來說①對于低能量的伽馬光子和原子序數高的吸收物質,主要形式是光電效應;②對于中等能量的伽馬光子和原子序數低的吸收物質,主要形式是康普頓效應;③對于高能量的伽馬光子和原子序數高的吸收物質,主要形式是電子對效應。

圖1 伽馬射線與物質的相互作用(引自IAEA,2003)
伽馬能譜測量是利用各種巖石和土壤中放射性元素含量差異所引起的放射性異常來尋找有用的礦產或查明地下地質構造的一種地球物理方法。根據巖石和土壤中的放射性元素分布特征,探討伽馬能譜測量的物理基礎。
3.1放射性元素在巖石中分布規律
不同類型的巖石,其放射性元素的含量是不相同的[10]。即使是不同時代的同一種巖性的巖石或不同地區的同一種巖性的巖石,其放射性元素含量也是不相同的[11],雖然這些放射性元素在巖石中含量差異很大,但是其具有以下分布規律[4,7,11-13]:
(1)在巖漿巖中,鈾、釷和鉀的含量往往隨著SiO2的含量增加而增長[14],如圖2所示。例如長英質巖石比超基性巖和鐵鎂質巖的放射性元素含量都要高[15]。在火成巖侵入的最后階段形成的巖石,由于溫度的驟降,不再具有這種地球化學行為,但是鉀含量很高,而鈾和釷含量卻相對的減少[15]。鈾和釷含量從基性巖到酸性巖呈增長的趨勢,釷隨著鉀的含量增加也呈顯著性的增長[16]。鈾、釷、鉀含量在巖漿巖中差異很大,在酸性巖漿巖中最高,中性巖稍高,基性、超基性巖中最低。

圖2 火成巖中K、U、Th平均含量變化隨著酸性(Si含量)相應增加(引自Dickson,Scott,1997)
(2)在花崗巖侵入體中,不同期次、不同巖相及不同脈中放射性元素含量都有差異。在巨大的花崗巖中,自外向內,自上向下,放射性元素的含量也有明顯差異;花崗巖中穿插的偉晶巖脈有較高的放射性,但細晶的放射性較低。由于放射性衰變原因,花崗巖體年代越新,其放射性越高,在晚期花崗巖中鈾含量較高(表1)。

表1 我國不同時期花崗巖中鈾含量
(3)沉積巖的放射性元素的含量比巖漿巖和變質巖都低。但是,不同的沉積巖中的鈾、釷、鉀含量也有很大的差異,其中頁巖和粘土巖中的含量較高,與酸性巖接近;砂巖中放射性元素含量變化大,主要取決于礦物顆粒組分;在石灰巖、石膏、鹽巖中的鈾、釷、鉀含量都比較低。在石灰巖中釷含量低,但是石灰巖中釷的含量會隨著雜質成份的增多而增加[17]。當沉積巖中富含磷、鐵和鋁的氫氧化物或有機物時,也會有相對較高的放射性。
(4)變質巖中的鈾、釷、鉀含量與變質前原巖的含量和變質程度有關。變質程度較高的變質巖放射性元素含量相對低,鈾、釷比值也隨著變質程度的增加而有規律的減小[18]。鈾、釷、鉀含量在炭質頁巖、白云母化—黑云母化片麻巖中相對含量都比較高,但是在大理巖、石英巖中含量偏低。這是由于在變質過程中,鈾、釷、鉀等發生了遷移,所以變質巖的放射性元素含量常比原巖的低。
3.2放射性元素在土壤中分布規律
土壤是巖石風化作用的產物,主要由礦物質和有機質組成,放射性元素的含量取決于形成土壤的基巖中放射性元素含量[4,7]。在表生風化條件下,鈾和釷因氧化作用而發生分離,鈾的活動性大,可被淋蝕搬運,在還原環境中富集[19],形成放射性異常;釷以機械風化、遷移為主,通常富集在土壤和風化巖石的殘積物中[20]。
綜上所述,根據巖石、土壤中的放射性元素的分布特征,為伽馬能譜測量提供了地質測量的依據,通過伽馬能譜測量能夠反應放射性元素分布規律。
地面伽馬能譜測量受深度的限制,所以僅能反映放射性元素在地表或近地表的富集的特征。因此,在野外伽馬能譜測量過程中,測得的能譜數據受幾何條件變化、覆蓋層、濕度、基質和溫度等條件的影響[21]。
4.1地形的影響
地形的幾何條件可能影響放射性元素含量的測量[8]。地面伽馬能譜儀是在平面幾何條件為2π的混凝土模型上標定的,而在實際測量工作中,地形起伏,測點通常會在凸起部位或是凹陷部位,這樣測得的值會不同,這是由于地面對探測器所張的立體角不同所致[11]。因此,在露頭很小、山脊、山谷和陡峭的斜坡上等不平的地表伽馬能譜測量時,可能對測量值有不同程度的影響[8]。為了克服地形的影響,應盡量把伽馬能譜儀放在比較平坦的位置進行測量或是對地形進行修正。
對臺階、沖溝、山包等不平坦的地形作相應修正,或將伽馬能譜探測器提高,可有效地減少地形對測量結果的影響[22],也可以用二維地形改正的方法對地形進行修正[23]。地形起伏程度、坡度大小和山坡的長短等因素可能對航空伽馬能譜測量結果產生影響,為此采用了比值校正法對地形校正[24]。為了克服地形的影響,應盡量把伽馬能譜儀放在比較平坦的位置上測量。
4.2覆蓋層的影響
覆蓋層為碎屑殘積、土壤層及植被等可能都會對地面伽馬能譜測量產生影響。根據初級光子注量率隨土壤厚度變化,地面伽馬能譜測量土壤的有效深度在100cm以內,土壤的密度也會對地面伽馬能譜測量產生影響[25]。如果地表風化碎屑、土壤層是原地殘積的而不是搬運的,那么地面伽馬能譜測量可認為是母巖物質[26],能夠反映測區的地質情況。所有類型巖石的風化作用都會導致鉀的丟失,但是對于長英質巖石,同樣會發生鈾和釷的丟失[15]。稠密的植被也會對伽馬輻射有明顯的衰減作用。覆蓋層厚度的影響,一般是隨著覆蓋層厚度的增加伽馬強度減弱[11]。
4.3濕度的影響
伽馬測量受土壤的濕度變化影響,當濕度較高時伽馬射線通量明顯衰減,伽馬測量值降低[21,27]。而孔隙度對能譜測量影響較小,隨著含水量增加能譜強度呈線性降低。用蒙特卡羅程序MCNP4B模擬計算了土壤樣品具有不同含水量時,獲得伽馬探測效率受含水量影響;高能伽馬射線探測受土壤含水量的影響小,而低能的受含水量影響大。在高溫干燥時比低溫濕潤時的伽馬測量值要高,主要是由于地表濕度增大,加強了水分對伽馬的吸收作用,使得實際測量值較低。因此,雨后不要立即進行能譜測量,否則影響到測量結果。
4.4溫度的影響
在其他條件不變時,溫度的變化能使光電倍增管增益發生變化,從而使儀器計數率也發生變化。隨著溫度升高,鈾、釷、鉀和總道參量的讀數都有增高,在20℃~40℃出現一段平緩段,而總道較短[28]。因此,在地面伽馬能譜測量的過程中,應當考慮溫度影響而采用能譜自穩措施,從而提高測量精度。
天然放射性核素中,鈾(238U)衰變成214Bi的過程中釋放出1.76 MeV的能量;同樣釷(232Th)衰變成208Th的過程釋放出2.61 MeV的能量,分別用來計算U和Th的含量[14];而40K衰變釋放的1.46 MeV的能量用來計算K的含量。不同放射性元素釋放出不同特征能量的伽馬射線[5,29],鈾、釷和鉀的平衡衰變系都有自己特征的線譜,傳統的伽馬能譜數據采集和處理方法是檢測3個或4個相對較寬的能量窗。鉀能量窗檢測的是40K衰變釋放出能量為1.46 MeV伽馬射線,鈾系和釷系伽馬射線能量窗都是其衰變子體核素產生的,這些能量窗是普遍認為最適合鈾、釷和鉀的測量,反映出總的放射性含量[8]。
因此,每一種伽馬輻射體都有自己特定能量的伽馬射線。在野外勘探過程中,伽馬能譜儀可以直接測量巖石或土壤中特定能量的伽馬射線,從而就能夠確定具有該譜線特征的放射性元素的存在。再將測得的伽馬強度與標準樣品的伽馬強度進行對比和計算,就可以確定該元素在巖石、礦石或土壤中的含量,這就是伽馬能譜測量的基本原理。
[1]Gerward L.Paul Villard and His Discovery of Gamma Rays[J]. Physics in Perspective,1999,1(4):367-383.
[2]Billings S D.Geophysics of Soil Mapping Using Airborne Gamma-Ray Spectrometry[D].The University of Sydney NSW,1998.
[3]田東風,龔健,伍鈞,等.核材料γ特征譜的探測和分析技術[M].北京:國防工業出版社,2005.
[4]程業勛,王南萍,侯勝利.核輻射場與放射性勘查[M].北京:地質出版社,2005.
[5]魯挑建,姜啟明.放射性地球物理勘查[M].哈爾濱工程大學出版社,2009.
[6]章曄,華榮洲,石柏慎.放射性方法勘查[M].北京:原子能出版社,1990.
[7]周蓉生,П.A.瓦岡諾夫.核方法原理及應用[M].北京:原子能出版社,1994.
[8] International Atomic Energy Agency.Guidelines for Radioelement Mapping Using Gamma Ray Spectrometry Data[M].IAEA,2003.
[9]陳伯顯,張智.核輻射物理及探測學[M].哈爾濱工程大學出版社,2011.
[10]Galbraith JH,Saunders DF.Rock Classification by Characteristics of Aerial Gamma-ray Measurements[J].Journal of Geochemical Exploration,1983,18(1):49-73.
[11]石玉春,吳燕玉,李秀季.放射性物探[M].原子能出版社,1986.
[12]傅良魁.應用地球物理教程:電法放射性地熱[M].地質出版社,1991.
[13]賈文懿,唐紅,方方,等.放射性勘查工作手冊[M].北京:地質出版社,1993.
[14]Wilford J,Minty B.The Use of Airborne Gamma-ray Imagery for Mapping Soils and Understanding Landscape Processes[J]. Developments in Soil Science,2006,31:207-610.
[15]Dickson B,Scott K.Interpretation of Aerial Gamma-ray Surveys-adding the Geochemical Factors[J].AGSO 10urnal of Australian Geology&Geophysics,1997,17(2):187-200.
[16]王劍鋒.鈾地球化學教程[M].北京:原子能出版社,1986.
[17]Galbraith JH,Saunders DF.Rock Classification by Characteristics of Aerial Gamma-ray Measurements[J].Journal of Geochemical Exploration,1983,18(1):49-73.
[18]劉峰,馬玉孝.變質巖填圖單位與放射性元素分布之間的關系[J].礦物巖石,1998(3):20-25.
[19]劉英俊,曹勵明,李兆膦.元素地球化學[M].科學出版社,1984.
[20]王繼倫,李善芳,齊文秀,等.中國金礦物探,化探方法技術的研究與應用[M].地質出版社,1997.
[21]于銘強,周鐳庭,吳慧山,等.放射性測量新技術(第一集)[M].地質出版社,1983.
[22]成智威,蘇川英,馬懷成,等.地形對車載就地HPGeγ譜儀測量結果的影響[J].輻射防護,2005(1):40-44,55.
[23]萬建華,熊盛青,范正國.航空伽馬能譜測量方法技術現狀與展望[J].物探與化探,2012(3):386-391.
[24]熊盛青.航空伽瑪能譜異常解釋中地形影響及校正方法研究[J].地學前緣,1998(2):172.
[25]汪洋.現場伽馬能譜測量影響因素及修正方法研究[D].中國地質大學(北京),2008.
[26] Beamish D.Environmental Radioactivity in the UK:the Airborne Geophysical View of Dose Rate Estimates[J].Journal of Environmental Radioactivity,2014,138(0):249-263.
[27]韓耀昭,劉騰耀.高精度放射性測量及其應用[J].物探與化探.1992(5):345-351.
[28]盧焱,戴麗君,王平.地面γ能譜測量影響因素初探[J].長春地質學院學報,1994(1):92-95.
[29]李舟波,孟令順,梅忠武.資源綜合地球物理勘查[M].地質出版社,2004.
P623
A
1004-5716(2016)08-0129-04
2015-08-12
2015-08-14
張選同(1982-),男(漢族),新疆阿克蘇人,新疆大學地質與礦業工程學院在讀碩士研究生,研究方向:礦床學,礦物學,巖石學。