地球化學(xué)綜合編圖

  （一）編圖試點(diǎn)研究簡(jiǎn)要史料。
 
  1973年3月，湖南地質(zhì)局物探大隊編制了《地球物理（化學(xué)）探礦圖件編制規范》。
 
  1974年4月，冶金部地質(zhì)司在廣西桂林舉行冶金地質(zhì)系統區域地球化學(xué)編圖及新方法、新技術(shù)交流會(huì )。
 
  1978年4月，物探所任天祥在《物探與化探》上發(fā)表了《關(guān)于1：20萬(wàn)區域化探資料整理利用方法的討論》，提出編制1：20萬(wàn)和1: 50萬(wàn)地區域化學(xué)圖的方法。
 
  1980年10月，謝學(xué)錦等發(fā)表了《皖浙贛邊區區域化探的資料整理與異常評價(jià)》，編制了皖浙贛邊區Cu、Pb、Zn等12個(gè)元素的地球化學(xué)圖。
 
  同年，桂林冶金地質(zhì)研究所與廣東地質(zhì)勘探公司940隊開(kāi)展了“廣東大陸區域化探編圖及找礦研究”。
 
  1985年10月，地礦部物化探所、江西物探大隊和北京計算中心合作完成的《贛西北地球化學(xué)試驗圖集》由地質(zhì)出版社出版，這是我國公開(kāi)出版的第一份地球化學(xué)圖集。
 
  1987年8月，地礦部物探所研制的IBMPC/XT應用開(kāi)發(fā)及通用成圖系統，移植和開(kāi)發(fā)了20多個(gè)常用化探數據處理程序，并研制了一套微機成圖系統。
 
  1988～1992年期間，北京計算中心承擔地質(zhì)行業(yè)科學(xué)技術(shù)發(fā)展基金資助的科研項目（劉如英、金諄浩、鄔寬廉、鄔可筠等），開(kāi)發(fā)研制了“計算機地球化學(xué)圖編圖系統”（小比例尺）(GCMCS)，在我國率先比較系統地研究小比例尺計算機地球化學(xué)編圖系統。
 
  1990年11月，地礦部北京計算中心（劉如英、金諄浩、鄔可筠）與遼寧物探隊（白希杰等）合作，首次用1: 20萬(wàn)區域化探數據為遼寧編制全省1：50萬(wàn)地球化學(xué)圖，編制了37種元素的等值線(xiàn)圖、色塊圖及異常圖，以及8種元素的綜合異常圖，可以全面反映遼寧省有關(guān)元素的分布特征和規律，是一份全省性的基礎地球化學(xué)資料。
 
  1993年，受地礦部勘技司委托，北京計算中心劉如英、金諄浩、趙玉祥等在中國地質(zhì)大學(xué)（北京）舉辦培訓班，推廣應用“計算機地球化學(xué)編圖系統（小比例尺）”。
 
  1993年12月，成都地質(zhì)學(xué)院完成的地面放射性測量工作站通過(guò)專(zhuān)家鑒定，以微機或筆記本機為平臺，實(shí)現了地面放射性測量數據的采集、傳輸、處理、成圖成像和解釋全過(guò)程的計算機化，該系統為國內首創(chuàng )。
 
  20世紀90年代后，大部分?。▍^、市）按?。▍^、市）域編制了?。▍^、市）地球化學(xué)圖，如浙江、新疆、湖北、貴州等（楊永忠，1997），同時(shí)按成礦區帶和大區域編制了區帶地球化學(xué)圖，如東昆侖、北祁連、內蒙古北山、華南等區域性地球化學(xué)圖。主要采用色塊、等值線(xiàn)和數據加密插值的成圖方法。
 
  （二）編制彩色地球化學(xué)圖件中一些問(wèn)題的處理。
 
  插值：由于輸出的彩色圖件為全面色，每個(gè)像元點(diǎn)都應具有數據及其對應的顏色，因此輸出彩圖之前必須求出每個(gè)像元點(diǎn)的插值。處理中常用的插值算法包括距離平方倒數、雙線(xiàn)性方法、雙三次樣條等。單一的插值算法往往存在缺點(diǎn)。距離平方倒數方法插值不連續，并且，當插值點(diǎn)與數據點(diǎn)的距離太近時(shí)，可能出現誤差；而雙三次樣條則在插值中出現高于數據中的最大值和低于數據中的最小值、甚至出現負值的不合理情況。為了避免這些問(wèn)題使用聯(lián)合算法插值，成功地解決了這一問(wèn)題。
 
  色階劃分：對已插值數據進(jìn)行分級，并對各數量級指定顯示的顏色。分級方法包括以標準差的倍數作為相鄰級別的差值、用累積頻率（或分位數）分級、使用等比級數分級等。
 
  色彩的調配：抖動(dòng)算法。由于彩色墨水只有4種顏色(青C、品紅M、黃Y和黑K)，由CMYK的不同組合至多可以獲得8種顏色（包括白色）。為了顯示更多的顏色，必須用加大象元、犧牲分辨率的方法。例如使用4×4的抖動(dòng)矩陣，每個(gè)象元可以獲得大約5 000種顏色，但此時(shí)對于分辨率為300 dpi的Novajet繪圖儀，分辨率降到300/4=75 dpi。
 
  (三）全國范圍的地球化學(xué)編圖。
 
  進(jìn)入20世紀90年代，全國區域化探掃面計劃宜測地區的任務(wù)基本完成，積累了豐富而系統的高質(zhì)量數據，許多?。▍^、市）據之編制了?。▍^、市）域的地球化學(xué)圖和圖冊。于是，編制全國范圍地球化學(xué)圖的任務(wù)就提上日程。
 
  1993年，中國勘查技術(shù)院立設立“中國小比例尺地球化學(xué)圖編制”項目，由地科院物化探研究所組織實(shí)施，在謝學(xué)錦和任天祥主持下，向運川負責編圖方法技術(shù)研究。這次編圖以1：20萬(wàn)和1：50萬(wàn)區域化探掃面數據為基礎，以1：5萬(wàn)和1: 2.5萬(wàn)圖幅為統計單元計算元素含量算術(shù)平均值，作為編圖的數據源。采用改進(jìn)的距離倒數加權法對數據進(jìn)行加密（插值），將數據劃分為16個(gè)色階，生成彩色像元，通過(guò)彩色噴墨繪圖儀繪制圖件。分別編制了39種元素的全國1：500萬(wàn)地球化學(xué)圖（1:2.5萬(wàn)圖幅數據集）和1：
 
  1 000萬(wàn)全國地球化學(xué)圖集（1:5萬(wàn)圖幅數據集）。編圖成果于1996年提交，這是我國首次編制的全國性39種元素的地球化學(xué)圖件。
 
  2001年，中國地質(zhì)調查局下達新的全國性編圖任務(wù)，即“全國區域地球化學(xué)系列圖編制”項目，由牟緒贊、奚小環(huán)負責，中國地質(zhì)調查局發(fā)展研究中心組織實(shí)施，謝學(xué)錦、任天祥等為顧問(wèn)，向運川為編圖組長(cháng)。這次編圖的特點(diǎn)是，直接以化探掃面的原始數據為編圖依據，首次匯集了全國28個(gè)?。▍^、市）的1: 20萬(wàn)和1: 50萬(wàn)區域化探資料，涉及1：20萬(wàn)圖幅1299個(gè)，1：50萬(wàn)圖幅18個(gè)，數據點(diǎn)約142萬(wàn)個(gè)，包括39種元素及氧化物的測試數據，覆蓋面積達650×104 km2。最終，分別編制出了39種元素與氧化物的1：500萬(wàn)全國地球化學(xué)圖和1：1 200萬(wàn)地球化學(xué)圖集。并且，在編圖同時(shí)建立全國區域地球化學(xué)數據庫。項目歷時(shí)4年，2005年6月提交編圖成果。
 
  由于采用原始數據編圖，數據量大，地球化學(xué)數據處理又有其特殊性，出現了一系列技術(shù)問(wèn)題需要解決。編圖工作采用中國地質(zhì)調查局研制開(kāi)發(fā)的“多元地學(xué)空間數據管理與分析系統GeoExpl”作為編圖平臺，使數據處理、編圖和輸出等過(guò)程一體化操作。第一，要分別對各元素進(jìn)行系統誤差的處理，以便得到很好的圖面效應。為此，制定了4條系統誤差校正的原則：①分析誤差源，所展示的數據誤差與周邊數據值具有明顯的臺階狀。
 
  ②數據誤差在空間上具有區域性特點(diǎn)，省、大區域、圖幅或分析批次。③在數據值的分布上，掩蓋了地球化學(xué)特征和地質(zhì)特征展布的延續性和規律性。④在數據校正方法上，盡可能地選擇線(xiàn)性校正，通過(guò)簡(jiǎn)單的計算可以復原數據。第二，數據網(wǎng)格化處理。大區域編圖難以利用原始點(diǎn)位直接構圖，需要對原始數據進(jìn)行網(wǎng)格化處理，通過(guò)對不同網(wǎng)格化參數的選擇、對比和研究，確定了效果較好的網(wǎng)格化處理方法，即：網(wǎng)格距4 kmX4 km的網(wǎng)度；數距搜索方式采用搜索半徑為12 km的圓域；數據計算方法采用改進(jìn)的距離倒數加權法，即指數距離冪倒數加權法。第三，地球化學(xué)專(zhuān)題圖的生成。應用GeoExpl提供的專(zhuān)題圖生成模塊，制作各元素等量線(xiàn)、等值區圖，以及相關(guān)的特征圖等。采用累計頻率與對數分級相結合，同時(shí)參考數據的特征值如平均值、中位數等，將數據分為19級，繼而劃分級次和色階，分級統計計算特征值，利用系統提供的專(zhuān)題圖創(chuàng )建模塊，生成地球化學(xué)等值線(xiàn)和等值區，分別以線(xiàn)圖層和區圖層的形式，作為系統顯示窗的圖形顯示數據。用這套方法制成的圖件，應用了GIS技術(shù)和現代數據處理技術(shù)，是近十幾年來(lái)計算機技術(shù)在地球化學(xué)數據管理、處理、成圖方法研究的結晶。