基于Python 編程的航空伽馬全譜分析方法的應用研究

馮春圓，王會波，周子陽，蔡文軍，楊春雷，朱林

（1.核工業航測遙感中心 河北省航空探測與遙感技術重點實驗室 中核鈾資源地球物理勘查技術中心（重點實驗室），河北 石家莊 050002；2.中核地礦科技集團有限公司，北京100013；3.華業（北京）科技有限公司，北京100013；4.中國冶金地質總局一局，河北 三河 065201）

傳統的航空伽馬能譜數據處理方法主要為國際原子能機構（1991）［1］推薦的“標準三窗法”，簡稱“三窗法”（TWM），即所謂的窗數據處理法，這種方法在地表勘查和淺部鈾礦床勘探過程中曾發揮了重要作用，而且仍在地質礦產勘查等許多領域得到較為廣泛的應用。該方法也是眾多數據處理軟件常用的數據處理方法，如oasis montaj、AGRS GeoProbe 等軟件，均提供標準三窗法數據處理模塊。隨著航空伽馬能譜測量向放射性本底調查、環境測量、核事故應急航空伽馬能譜測量、地下水及溫泉調查以及大氣氡評價和修正等方面的應用［2-9］。數據處理方法也從單一的窗數據處理走向譜數據處理發展。

航空伽馬全能譜分析方法，簡稱“全譜法”（FSA），是基于全能譜數據進行的一種數據處理方法。就是利用航空伽馬能譜測量中所有能譜數據，而不是常規的窗數據，直接進行飛機和儀器本底、宇宙輻射本底、航測高度等的各項修正，而將大氣氡視為變量，與地面輻射源的鉀、鈾、釷元素含量以及其他核污染源（人工核素）的濃度或活度一起進行計算的分析方法［10］。該方法能有效解決窗數據處理中大氣氡修正引入較大不確定度問題，前人已在原理和應用方面作了大量的研究［11-14］。……