發布時間：所屬分類：工程師職稱論文瀏覽：1次
摘 要： 摘 要 為了構建深部砂巖型鈾礦化信息提取技術，本文論述了基于放射性地球物理方法提取深部鈾礦化信息的理論基礎，采用氡異常二維柵格襯度、分形濾波、氡與計算氡比值等方法對我國二連盆地及準噶爾盆地實測土壤氡濃度數據進行了處理，同時，應用證據權重法，集成地質、
摘 要 為了構建深部砂巖型鈾礦化信息提取技術，本文論述了基于放射性地球物理方法提取深部鈾礦化信息的理論基礎，采用氡異常二維柵格襯度、分形濾波、氡與計算氡比值等方法對我國二連盆地及準噶爾盆地實測土壤氡濃度數據進行了處理，同時，應用證據權重法，集成地質、土壤氡、車載伽馬全譜、航空伽馬能譜等成礦有利信息圖層，開展了多元信息定量化提取技術應用研究，通過數據處理，表層信息得到了有效壓制，排除了氡異常的漂移或假異常，處理后的異常信息與礦體在地表的垂直投影及控礦構造正相關，為我國北方深部砂巖型鈾礦找礦工作提供了技術支撐.
關鍵詞 土壤氡濃度; 氡襯度; 分形濾波; 多元信息提取; 深部鈾礦化信息
0 引 言
砂巖型鈾礦產于中-新生代沉積盆地且淺部鈾礦已基本查明，尋找深部砂巖型鈾礦正處在全面攻堅階段，各種勘查手段正在向高精度、多元化、精細化發展且不斷創新提高. 我國的鄂爾多斯盆地、伊犁盆地、二連盆地、巴丹吉林-巴音戈壁盆地、塔里木盆地、松遼盆地、準噶爾盆地、四川盆地等沉積盆地具有較大的砂巖型鈾礦找礦潛力，隨著大營鈾礦的發現，引起新一輪深部勘查熱潮，但是，深部勘查技術仍是當前面臨的難題. 從已發現的砂巖型鈾礦床來看，產于這些沉積盆地的砂巖型鈾礦大多是隱伏盲礦，上覆蓋層較厚( 一般 200 m 左 右) ，且大多沉積盆地具有深部( 礦 體 埋 深 大 于 500 m) 找礦潛力. 然而，砂巖型鈾礦較厚的上覆蓋層卻增加了放射性物探找礦難度，因為砂巖鈾礦體的強放射性異常信息在上覆蓋層傳遞過程中經歷了圍巖及蓋層介質的層層阻隔，到達地表時表征深部鈾礦化信息十分微弱，甚至被表層信息掩蓋，深部放射性地球物理找礦標志不明顯. 伽馬能譜測量方法基本被認為“無效”，具有深穿透能力的氡及其子體測量方法，由于受地表淺層信息及氣象等因素的干擾，指示深部鈾礦化信息不明顯，常出現假異常或異常漂移，由于這些原因，放射性地球物理勘查方法也曾經一度被鈾礦勘查領域內地質工作者懷疑甚至摒棄. 進入 21 世紀后，隨著砂巖型鈾礦的不斷發現及放射性地球物理方法在已知礦體上方地表的實驗分析研究及對蓋層地質環境認識，提出了剖面 “雙峰”、平面“圈閉類環狀”的表征深部鈾礦化信息的氡異常特征( 李必紅等，2007) ，這些氡異常特征在隨后發現的砂巖型鈾礦上的不斷試驗中被證實，但出現這種情況的原因還在探索中. 在未知區域的應用中，如何識別和提取指示深部鈾礦信息的放射性地球物理找礦標志是砂巖型鈾礦“攻深找盲”的關鍵問題和難點. 解決問題可從兩方面入手，第一是采集方面，第二是數據處理方面，本文著重介紹數據處理方面.
1 表征深部鈾礦化信息的放射性理論基礎
關于放射性測量表征深部鈾礦化信息的機理研究已經經歷了較長的時間，主要是圍繞鈾衰變子體從礦體中釋放、遷移及富集的過程，涉及粒子反沖( Porstendrfer，1994) 、團簇理論( Mandich，1996; Hunter and Jarrold，1999; ButterveckDempewolf et al. ，2000; 賈文懿等，2000; 樂仁昌等，2002) 、擴散( Semkow et al. ，1990) 、對流( Al-Tamimi and Abumurad， 2001) 、接力傳遞( 吳慧山等，1997) 等方面的理論或觀點，這些理論或觀點大多經過實驗模擬，基本是目前比較公認的.
鈾系 列 分 鈾 組 ( 238 U、234 Th、…) 和 鐳 組 ( 226 Ra、222 Rn、 218 Po、214 Pb、214 Bi、…210 Po、…) ，鈾礦放射性勘查主要關注鐳組核素，其中，氡( 222 Rn) 是鈾238 U 子體中唯一的以氣體存在的核素，是深部找礦重要的示蹤劑. 表層土壤中氡由兩部分構成: 其一是原生氡，由測點位置有限小空間內土壤中正常放射性鈾系核素238 U，…，238 Ra 的衰變產生的; 其二是外來氡，由非測量點位置的放射性鈾系核素238 U，…，238 Ra 衰變產生的，包括深部接力傳遞遷移( 吳慧山等，1997) ，地表運移和大氣氡的滲入. 對于野外現場實測的每個測點，其土壤中氡是空間氡信息疊加的結果，如何在這種疊加的氡信息中提取來自深部遷移的氡是深部找礦的關鍵，筆者認為以數據采集和數據處理為突破口，不斷創新理論和技術手段是一種可行路徑.
2 深部鈾礦化信息提取技術
2. 1 氡異常分形濾波技術
分形濾波技術是目前數據處理技術研究熱點，成秋明等提出并發展了分形技術在地球化學中的應用( Cheng et al. ， 1994; Cheng，2003; 陳永清等，2006; 成秋明等，2009) ，本論文在此基礎上，將分形濾波技術應用于實測土壤氡濃度的數據處理中，通過傅立葉變換將異常與背景分離，選取合適的濾波器去掉背景值，再進行逆變換后的數據壓制了背景值的干擾，使異常分帶性明顯，可有效地提取深部氡濃度信息.
如圖 1，內蒙古某一深部砂巖型鈾礦礦區實測的土壤氡濃度數據經過分形濾波處理結果，原始土壤氡濃度等值圖 ( 圖 1a) 顯示礦體在地表的垂直投影與氡濃度低暈場對應，氡異常暈或氡高暈場處在礦體側面，經過分形濾波處理后的等值圖上，礦體在地表垂直投影范圍內顯示密度較大的串珠狀或線狀異常，表明分形濾波技術處理后的氡濃度異常密度可用于評價深部鈾礦化信息.
2. 2 氡濃度與計算氡濃度比值法
伽馬能譜測量表征淺層信息，伽馬能譜現場實測的鈾當量含量通常是通過測量的氡子體214 Bi( 1. 76 MeV) 特征能量譜線換算獲得的，然而，實測獲得的數據信息中包含深部的信息量小，深部信息提取難度大. 假設土壤層中氡氣沒有外來的，全是土壤層正常放射性核素衰變產生，且全部析出的氡氣都保存在土壤中，則，伽馬能譜獲得的鈾當量含量與土壤氡氣來自土壤正常核素產生的氡氣濃度呈線性關系. 也就是說，伽馬能譜所獲鈾當量含量與土壤氡氣測量的氡濃度背景值呈線性關系. 于是，采用伽馬能譜鈾當量含量數據計算氡濃度( 稱之為計算氡) ，反應淺層氡濃度分布信息特征，這樣，通過實測氡濃度與計算氡濃度的比值來評價氡的深部異常信息.
應用實測氡濃度與計算氡濃度比值來評價氡濃度異常時，一般當實測氡濃度比計算氡濃度高時，可能由深部放射性核素引起，即異常有可能與深部礦化有關.
采用這種方法處理了新疆準噶爾某一地區土壤氡濃度，如圖 2，可以看出，土壤氡濃度與計算氡濃度比值越大的數據分布區域對深部鈾成礦信息越有利，但是，計算氡濃度是假設浮土層為飽和情況，而實際土壤空氣與大氣之間是流動的，往往是處于動態的平衡狀態，所以這種方法還存在缺陷.
2. 3 氡異常二維柵格襯度技術
襯度已被廣泛應用于醫學、光學、電子學、生命科學等領域，用于突出目標信息或顯示信息異常，也是放射性異常評價的一個重要指標，測點氡濃度與測點周邊一定范圍內各氡濃度平均值的比值可突出深部異常信息( LI et al. ，2014) . 本論文是基于 ArcGis 平臺和襯度概念，將實測的散點數據插值成二維柵格數據，選取適當的滑動濾波器，對柵格數據進行滑動平均濾波，獲取數據區域范圍內背景值柵格數據，對二維柵格數據進行小波變換處理，獲取數據范圍內異常柵格數據，通過異常柵格數據和背景柵格數據的比，進行深部鈾礦化信息評價. 筆者對二連盆地實測土壤氡濃度進行處理，達到了比較滿意的效果，圖 3 為二連盆地已知砂巖型鈾礦區土壤氡濃度及氡異常襯度等值圖，原始氡濃度等值圖( 圖 3a) 氡異常暈與已揭露的工業孔的漂移較大，在土壤氡濃度襯度異常等值圖中( 圖 3b) 顯示明顯的兩條串珠狀線性氡襯度異常，該異常該地段的地質構造對應，且主礦體在地表的 垂直投影處在這兩條串珠狀線性異常中間靠近北側異常.
2. 4 多元信息定量提取
對于出露地表的淺層鈾礦，采用伽馬能譜測量或者伽馬總量測量即可確定礦體位置和規模. 對于深部鈾礦體，單一方法往往存在多解性，深部鈾礦化信息在傳遞至地表的過程中受遷移營力、蓋層介質環境、核素衰變、地表信息的干擾等條件的約束，很難垂直傳遞，加大了找礦難度，采用多元信息提取可極大減少多解性帶來的找礦不確定性.
相關知識推薦：地球物理學期刊有哪些
為此，筆者在二連盆地某一礦區聯合開展車載伽馬能譜和瞬時測氡技術研究，車載伽馬能譜測量線距 250 ～ 500 m，測量周期 5 s，車速 20 km/h，土壤氡氣測量采用點線網格為 100 m × 500 m. 將該地區實測土壤氡濃度，車載伽馬能譜鈾當量含量、釷當量含量、鉀含量經過數據處理，形成深部成礦有利的氡信息圖層，鈾當量含量的變異系數、差量鈾、集成鈾圖層、等伽馬能譜信息圖層，結合該地區調研的航空伽馬能譜信息圖層和地層組合熵圖層，組成綜合信息系列圖層. 這里，鈾當量含量變異系數是統計單元內鈾當量含量均方差與鈾當量含量均值的比值. 差量鈾是指統計窗口內各測點鈾當量含量與該窗口內鈾當量含量的平均值的差值. 集成鈾是指鈾富集系數、鈾襯度、全譜分解分形鈾三者歸一化后的加和信息，鈾富集系數為鈾、鈾釷比、鈾鉀比三者的乘積，全譜分形分解鈾是通過奇異值分解伽馬全譜處理后( Hovgaard and Grasty，1997) ，將全譜數據中噪聲和突變等干擾信號除去，應用分形方法提取異常場的方法.
采用證據權重方法進行這些信息圖層綜合處理后，得到后驗概率如圖 4，其中后驗概率越大，成礦概率越大，成礦越有利，從圖 4 可以看出已揭露的工業孔基本處在后驗概率大于 0. 65 的范圍內.
3 結論與探討
3. 1 深部鈾礦化信息傳遞至地表富集，其過程復雜多變，影響因素較多，除了蓋層地質環境外，還與鈾及其子體核素的衰變、遷移以及鈾鐳平衡等因素有關，從放射性地球物理實測數據中提取深部鈾礦化信息關鍵是如何壓制淺層信息.
3. 2 實測土壤氡濃度是土壤本身固有氡和外來氡的信息疊加，可通過數據處理來提取來自深部的氡異常信息. 采用氡異常二維柵格襯度、分形濾波、氡濃度與計算氡濃度比值對實測土壤氡濃度數據進行了處理，表明該方法具有較好的應用效果; 多元信息提取是鈾礦“攻深找盲”的趨勢，具有科學價值和研究意義.
3. 3 放射性地球物理的數據處理只是提取深部鈾礦化信息的解決途徑之一，數據采集也是提取深部鈾礦信息一種途徑，而且，深部鈾礦目標層結構，蓋層地質環境、核素遷移也是深部鈾礦勘查需要解決的問題，可進一步開展鈾子體核素在礦體蓋層空間遷移的三維動態模擬及放射性地球物理勘查、深穿透地球化學勘查、普通地球物理勘查融合式研究. 但是，深穿透地球化學勘查過程繁瑣，周期長，還有待提高.——論文作者：李必紅，吳慧山，趙丹，楊龍泉
