MapGIS與Excel相結(jié)合在新疆恰爾隆一帶化探數(shù)據(jù)處理中的應(yīng)用

陳仕炎

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

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

本文以新疆恰爾隆一帶1∶5萬水系沉積物地球化學(xué)測量為例，從背景值和異常下限的確定、不同地球化學(xué)圖件繪制等方面簡述了化探數(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ū)中南部），河流水量比較充沛，為開展水系沉積物測量的有利地段。

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

2 成圖準(zhǔn)備

2.1 數(shù)據(jù)準(zhǔn)備

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

2.2 地球化學(xué)背景值及異常下限確定

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

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

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

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

3 生成地球化學(xué)異常分布圖

（1）進(jìn)入 MapGIS 6.7主菜單，點(diǎn)擊空間分析模塊中“DTM分析”子模塊。

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

（3）在平面等值線圖繪制命令中導(dǎo)入＊.Grd文件，進(jìn)行等值線圖的繪制。

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

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

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

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

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

4 繪制綜合異常圖

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

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

5 結(jié)語

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

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

［2］張釗，韋龍明，陳三明，等.MAPGIS在地質(zhì)填圖及化探數(shù)據(jù)處理中的應(yīng)用［J］.物探化探計(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］紀(jì)宏金，林瑞慶，周永昶.關(guān)于若干化探數(shù)據(jù)處理方法的討論［J］.地質(zhì)與勘探，2001，37（4）：56－59.

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

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

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

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