利用襯度法繪制化探綜合異常圖

李　歡　徐國志　劉宏哲　馮　盼

(河北地質大學研究生院)

?

利用襯度法繪制化探綜合異常圖

李歡徐國志劉宏哲馮盼

(河北地質大學研究生院)

摘要在化探掃面工作的基礎上，現階段我國正進行的整裝勘查區1∶50 000區域地質調查工作均配備了地球化學測量工作，該工作已成為我國找礦工作的重點。在分析常規綜合異常圖繪制方法的基礎上，以內蒙古錫林浩特某地的區域地質礦產調查工作為例，詳細分析了利用襯度法繪制綜合異常圖的具體流程。實踐表明，采用該方法繪制的綜合異常圖不僅圖面清晰，而且有利于指導找礦工程的布設。

關鍵詞襯度法綜合異常圖單元素襯度值總襯度值

地球化學測量可分為系統采樣、測試分析、數據整理等3個部分，其中，數據整理工作具有承上啟下的關鍵性作用，既是對前期野外工作的系統總結，也可指導下一步的找礦工作。本研究在詳細分析現階段我國地質行業化探綜合異常圖繪制方法的基礎上，提出采用襯度法繪制化探綜合異常圖的思路。

1常規綜合異常圖繪制方法

現階段化探工作中所涉及異常圖主要包括單元素異常圖、組合異常圖和綜合異常圖。單元素異常圖的繪制步驟包括原始數據整理、離散數據網格化和平面等值線繪制。組合異常圖是建立在單元素異常圖之上，依據工作區所屬成礦帶元素特征、地質背景、元素遷移規律、元素富集與貧化特征，結合各元素分析數據的相關性特征，最終確定元素組合，將每一種組合中的相關單元素異常進行添加，即可形成組合異常圖。綜合異常圖是指將所有分析元素的異常特征全部顯示于同一幅異常圖中，該方法是基于已完成的組合異常圖，將所有繪制完成的組合異常圖進行疊加，即可形成綜合異常圖。

2襯度法繪制綜合異常圖思路

綜合異常圖是在組合異常圖的基礎上，研究幾組元素的空間分布規律的綜合性圖件[1-4]。該定義在一定程度上與組合異常圖無太大差異，僅增加了元素數量，同時也增加了圖面負擔，也無法完全代表研究區內所有元素的地球化學特征。為在同一幅綜合異常圖中更全面、系統地顯示研究區內所有元素的地球化學特征，需將所有元素的含量值進行疊加。實驗室內對元素進行常規分析時，某些元素在地殼中的豐度值相對較高，分析值可達10-6數量級；某些元素在地殼中的豐度值相對較低，分析值可達10-9數量級。若僅簡單的對所有元素的分析值進行疊加，無疑會減弱某些含量相對較低(10-6數量級)元素的異常強度。為消除數量級不同而造成的影響，在利用襯度法繪制綜合異常圖時，提出了“單元素襯度值(single element contrast value，SECV)”和“總襯度值(Total contrast value，TCV)”的概念。單元素襯度值可定義為某元素分析數據與該元素異常下限的比值，當單元素襯度值≥1以上時，表明該采樣點的某元素分析結果大于異常下限，可形成該元素異常，將該值定義為異常襯度值(Anomaly contrast，AC)；當單元素襯度值<1時，表明該采樣點的某元素異常下限大于分析結果，無法形成元素異常。總襯度值可定義為某采樣點所有分析元素異常襯度值之和與異常元素數量的比值。單元素襯度值為分析結果與異常下限的比值，屬于無量綱參數，故通過取平均值得到的總襯度值也屬于無量綱參數，因而有利于不同含量級別元素之間的疊加。所謂襯度法，即通過對分析數據進行處理，得到每個采樣點的總襯度值，以總襯度值≥1為異常下限，勾繪等值線，繪制綜合異常圖，通過設置不同類型的步長(表1)，繪制綜合異常圖進行比對，綜合考慮成圖效果，同時結合實際成礦條件，認為以2n為步長生成的綜合異常圖既美觀又清晰。

表1　步長設置對比

3應用實例

內蒙古錫林浩特某地區已完成了區域地質礦產調查工作，化探樣品分析時共分析了16種元素，且通過化探方法初步圈定了規模相當的鉛鋅多金屬礦產。通過對該區域化探樣品分析數據的整理和圖件的編制，采用襯度法繪制綜合異常圖，流程見圖1。

圖1　綜合異常圖繪制流程

3.1原始數據整理

將化探采樣點位和樣品分析數據存儲于Excel表中，其中X坐標位于第1列，Y坐標位于第2列，16種元素分析數據位列其后。

3.2確定異常下限

將上式計算結果轉換為非對數形式，即可得到異常下限值。

上述計算過程可能在繪制單元素異常圖時已完成。對于異常下限的確定，建議采用SPSS軟件進行數據處理。

3.3數據分析

(1)計算單元素襯度值，即某元素的分析結果與該元素異常下限的比值，可利用Excel表完成計算。

(2)區分采樣點有無異常，即當單元素襯度值≥1時，在計算總襯度值時采用該數據；當單元素襯度值<1時，在計算總襯度值時舍棄該數據，并將其單元素襯度值取為0。通過區分采樣點有無異常，可得到16種元素的單元素襯度值表。該過程可采用Excel表中IF函數實現。

(3)計算總襯度值，該值為某采樣點所有分析元素異常襯度值之和與異常元素數量的比值，可在單元素襯度值表中進行計算。異常襯度值數量可通過Excel表中的COUNTIF函數實現。異常襯度值之和即為單元素襯度值表中所有數值之和。通過上述運算，在計算總襯度值時僅需進行簡單的函數運算即可完成。在總襯度值計算完畢后，新建一個Excel表，其中，X坐標位于第1列，Y坐標位于第2列，總襯度值Z位于第3列，將該Excel表另存為 “.txt”格式。

3.4離散數據網格化

對分析數據進行網格化處理，使原本并無關聯的相鄰采樣點產生聯系，是繪制等值線圖的核心。本研究利用MapGIS軟件的DTM分析功能完成離散數據的網格化處理。

3.5平面等值線圖繪制

平面等值線圖繪制是對網格化數據文件的圖面顯示，將數學模型轉換為直觀圖形，該過程將形成初步的綜合異常圖。詳細步驟為①MapGIS軟件→空間分析→DTM分析→Grd模型→平面等值線圖繪制，選取已經生成的網格化數據文件(*.Grd)；②參數設置，該操作步驟至關重要，不僅關系到所生成的綜合異常圖的正確性，而且關系到圖面整體的美觀性，可在設置等值線參數對話框中選擇等值線套區、繪制色階、保留邊界線、等值線光滑處理(光滑度為高程度)、軸向標尺。由于本研究采用的是 1∶ 50 000 區域地質礦產調查工作所得的化探分析數據，故在制圖幅面中應進行數據投影變換。在等值層值中，等值線層分段參數起始值設置為1，設置等值層值時，以2n為步長無法設置步長增，故需手動填寫等值層值，具體數值可參考表1。注記參數設置時要求每條等值線都須有標記，設置完成后如圖2所示。上述參數設置完畢后，點擊“確定”，即可生成等值線圖，并保存為相應的點、線、面文件。

圖2　等值線參數設置界面

3.6圖面修飾

將生成的等值線圖與地質圖進行套合，調整圖中字體大小、角度、疊壓情況，使圖面總體美觀。

4結語

結合內蒙古錫林浩特某地區域地質礦產調查工作成果，提出了利用襯度法繪制化探綜合異常圖的基本思路，并對具體的繪制步驟進行了詳細分析。研究表明，利用襯度法繪制化探綜合異常圖不僅操作簡單，而且可減輕圖面負擔，提高圖幅的清晰度，有一定的實用價值。

參考文獻

[1]馬振東.礦產資源地球化學模型建立與定量預測研究——以長江中下游銅多金屬成礦帶為例[R].武漢：中國地質大學(武漢)，2008．

[2]馬振東.礦產資源地球化學模型建立與定量預測研究——以西藏岡底斯銅多金屬成礦帶為例[R].武漢：中國地質大學(武漢)，2009．

[3]馬振東.礦產資源地球化學模型建立與定量預測研究方法指南[R].北京：中國地質調查局發展研究中心，2010．

[4]中華人民共和國地質礦產部.DZ/T 0011—911∶50 000地球化學普查規范[S].北京：中國標準出版社，1992.

(收稿日期2015-11-01)

李歡(1990—)，男，碩士研究生，050031 河北省石家莊市裕華區槐安東路136號。