陳仕炎

（長沙有色冶金設(shè)計研究院有限公司，長沙410011）

MapGIS地理信息系統(tǒng)是武漢中地信息工程有限公司研制的，集數(shù)字制圖、數(shù)據(jù)庫管理以及空間分析為一體的空間地理信息系統(tǒng)。它以其功能強大而廣泛應(yīng)用于城市規(guī)劃、測繪、地質(zhì)勘查、資源管理等領(lǐng)域，并已成為我國各領(lǐng)域數(shù)字化建設(shè)的首選軟件之一。在地質(zhì)找礦方面，該系統(tǒng)在提升地質(zhì)勘查技術(shù)手段上發(fā)揮了重要作用，大大提高了工作效率，而且保證了空間數(shù)據(jù)的精確性。

本文以新疆恰爾隆一帶1∶5萬水系沉積物地球化學測量為例，從背景值和異常下限的確定、不同地球化學圖件繪制等方面簡述了化探數(shù)據(jù)處理的方法和步驟，通過MapGIS軟件與Excel電子表格結(jié)合使用在化探數(shù)據(jù)處理過程中取得了良好的效果。

1 工作概況

工作區(qū)位于新疆西昆侖山西段北坡，塔里木盆地南緣。地勢南高北低，海拔一般為3 000～5 500m，最大相對高差在1 500～2 000m左右。區(qū)內(nèi)水系較發(fā)育，均屬內(nèi)陸河的源頭水系，主干水系為塔什庫爾干河和恰爾隆河。支流水系比較發(fā)育，水系組合格局多呈現(xiàn)為樹枝狀（調(diào)查區(qū)北部）及網(wǎng)格狀（調(diào)查區(qū)中南部），河流水量比較充沛，為開展水系沉積物測量的有利地段。

本次化探取樣共計2 706件，化探樣品測試由新疆地礦局中心實驗室測試完成，分析項目包括：Au、Ag、As、Sb、Cu、Pb、Zn、W、Sn、Mo、Bi、Ni、Co、U、Th等15個元素。

2 成圖準備

2.1 數(shù)據(jù)準備

成圖前逐點采集樣點的平面坐標。數(shù)據(jù)排列順序為X坐標、Y坐標、元素含量（見表1），圖面坐標與實際坐標要進行轉(zhuǎn)換。分別建立各元素含量數(shù)據(jù)的電子表格文件，在Excel電子表格中另存文件，文件類型為＊.txt。

2.2 地球化學背景值及異常下限確定

考慮到方法的實用性、有效性和易操作性，通過幾種方法在工作區(qū)的試驗對比，迭代法確定的背景值及異常下限較低，更有利于突出弱異常。因此，本工作區(qū)背景值和異常下限的確定選用迭代法。

應(yīng)用迭代法在Excel中計算異常下限的過程：（1）計算全區(qū)各元素原始數(shù)據(jù)的均值（X1）和標準偏差（Sd1）；（2）按X1＋nSd1的條件剔除一批高值后獲得一個新數(shù)據(jù)集，再計算此數(shù)據(jù)集的均值（X2）和標準偏差（Sd2）；（3）重復第二步，直至無特高值點存在，求出最終數(shù)據(jù)集的均值（X）和標準偏差（Sd），則X做為背景值C0，X＋nSd（n根據(jù)情況選1.5或2，3）做為異常下限Ca。

表1 測區(qū)內(nèi)部分水系沉積物樣品原始數(shù)據(jù)

表2 工作區(qū)元素地球化學特征值及參數(shù)表

3 生成地球化學異常分布圖

（1）進入 MapGIS 6.7主菜單，點擊空間分析模塊中“DTM分析”子模塊。

（2）在“Grd模型”下拉菜單中單擊“離散數(shù)據(jù)網(wǎng)格化”命令，根據(jù)需要修改網(wǎng)格參數(shù)調(diào)整網(wǎng)格疏密程度，選擇網(wǎng)格化方法，導入＊.txt文件，生成網(wǎng)格數(shù)據(jù)（＊.Grd文件）并保存（圖1）。

（3）在平面等值線圖繪制命令中導入＊.Grd文件，進行等值線圖的繪制。

圖1 離散數(shù)據(jù)網(wǎng)格化顯示視窗

（4）打開生成的網(wǎng)格數(shù)據(jù)文件，在設(shè)置等值線參數(shù)對話框中（圖2），選擇等值線套區(qū)、保留邊界線、繪制色階和等值線光滑處理等，并分別對等值層的起始值、終止值、步長值和線參數(shù)等項進行設(shè)置，點擊“確定”按鈕即可顯示元素地球化學異常分布圖，結(jié)果如圖3所示。

（5）保存所生成的點文件、線文件和區(qū)文件。

圖2 設(shè)置等值線參數(shù)視窗

圖3 測區(qū)部分地段金元素地球化學異常分布圖

4 繪制綜合異常圖

依照上述生成元素地球化學異常分布圖的步驟，并且根據(jù)工作區(qū)地球化學特征，確定元素組合，把元素組合中各元素的單元素異常添加在一起，然后在MapGIS編輯子系統(tǒng)創(chuàng)建工程，并在工程管理器窗口依次添加各元素的地球化學異常分布圖等相關(guān)文件；最后點擊窗口中“1∶1”按鈕，在窗口位置便可看到各元素地球化學等值線圖疊加在一起，形成地球化學綜合異常分布圖。測區(qū)部分地段綜合異常圖見圖4。

圖4 測區(qū)部分地段化探綜合異常圖

5 結(jié)語

Excel電子表格具有強大的數(shù)理統(tǒng)計功能，可以快速統(tǒng)計計算出元素地球化學特征參數(shù)及異常下限，而MapGIS軟件在化探數(shù)據(jù)處理中具有獨到的功用，如利用它自動生成各元素的地球化學等值線圖，直接編制綜合異常圖等，具有操作簡單、圖件精確、快速高效的特點。

［1］袁義生，劉應(yīng)忠，羅明學，等.應(yīng)用 MAPGIS制作地球化學圖單元素異常圖及綜合異常圖［J］.貴州地質(zhì)，2007，24（4）：156.

［2］張釗，韋龍明，陳三明，等.MAPGIS在地質(zhì)填圖及化探數(shù)據(jù)處理中的應(yīng)用［J］.物探化探計算技術(shù)，2010，32（2）：221－224.

［3］趙榮軍.不同方法在欒川北部化探數(shù)據(jù)處理中的應(yīng)用［J］.地質(zhì)與勘探，2006，42（3）：67－71.

［4］吳錫生.化探數(shù)據(jù)處理方法［M］.北京：地質(zhì)出版社，1992：25－26.

［5］李明，崔永文，李俊，等.MapGIS與Excel相結(jié)合在遼北化探數(shù)據(jù)處理中的應(yīng)用［J］.科研探索與知識創(chuàng)新，2009，12（下）：83－84.

［6］紀宏金，林瑞慶，周永昶.關(guān)于若干化探數(shù)據(jù)處理方法的討論［J］.地質(zhì)與勘探，2001，37（4）：56－59.

［7］宋運紅，賈大成，李健，等.河北承德三溝地區(qū)化探數(shù)據(jù)處理與異常特征［J］.吉林大學學報：地球科學版，2008，38（增）：167－169.

［8］戴慧敏，宮傳東，鮑慶中，等.區(qū)域化探數(shù)據(jù)處理中幾種異常下限確定方法的對比——以內(nèi)蒙古查巴奇地區(qū)水系沉積物為例［J］.物探與化探，2010，34（6）：782－786.

［9］春乃芽.如何利用Excel處理化探數(shù)據(jù)［J］.物探化探計算技術(shù)，2006，28（3）：272－276.

［10］楊世瑩.Excel數(shù)據(jù)統(tǒng)計與分析范例［M］.北京：中國青年電子出版社，2005：18－27.