戴慧敏, 趙 君, 楊忠芳, 宮傳東, 鄭春穎, 孫中任

1)中國(guó)地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院, 北京 100083; 2)沈陽(yáng)地質(zhì)礦產(chǎn)研究所, 遼寧沈陽(yáng) 110034; 3)遼陽(yáng)市鏵子鎮(zhèn)政府, 遼寧遼陽(yáng) 111300

基于地球化學(xué)背景的多圖幅系統(tǒng)誤差校正
——以區(qū)域地球化學(xué)調(diào)查數(shù)據(jù)Au元素為例

戴慧敏1, 2), 趙 君2), 楊忠芳1), 宮傳東3), 鄭春穎2), 孫中任2)

1)中國(guó)地質(zhì)大學(xué)(北京)地球科學(xué)與資源學(xué)院, 北京 100083; 2)沈陽(yáng)地質(zhì)礦產(chǎn)研究所, 遼寧沈陽(yáng) 110034; 3)遼陽(yáng)市鏵子鎮(zhèn)政府, 遼寧遼陽(yáng) 111300

針對(duì)區(qū)域地球化學(xué)圖件中出現(xiàn)的含量等值線(xiàn)環(huán)繞子區(qū)邊界現(xiàn)象, 以 Au的多個(gè)圖幅拼接為例, 使用了分幅標(biāo)準(zhǔn)化法、襯度返回法及定和化法 3種方法對(duì)地質(zhì)大調(diào)查以來(lái)新方法獲得的某成礦帶區(qū)域地球化學(xué)數(shù)據(jù)進(jìn)行系統(tǒng)誤差校正, 校正后編制地球化學(xué)等量線(xiàn)圖顯示, 所采用的3種校正方法對(duì)于研究區(qū)數(shù)據(jù)存在的系統(tǒng)誤差都不能達(dá)到較好的調(diào)平效果。針對(duì)研究區(qū)數(shù)據(jù)特點(diǎn), 嘗試在分幅標(biāo)準(zhǔn)化法和襯度返回法基礎(chǔ)上提出基于地球化學(xué)背景的誤差校正方法, 進(jìn)一步對(duì) 4種方法校正后數(shù)據(jù)參數(shù)特征及編制地球化學(xué)等量線(xiàn)圖進(jìn)行對(duì)比?；诘厍蚧瘜W(xué)背景的誤差校正方法對(duì)研究區(qū)達(dá)到了較好的數(shù)據(jù)校正效果, 并對(duì)校正前后數(shù)據(jù)進(jìn)行相關(guān)對(duì)比和檢驗(yàn), 滿(mǎn)足化探數(shù)據(jù)中存在的系統(tǒng)誤差屬于線(xiàn)性誤差這一基本假設(shè), 在實(shí)際應(yīng)用中, 可供地球化學(xué)數(shù)據(jù)系統(tǒng)誤差校正靈活地使用, 以取得更好的應(yīng)用效果。

地球化學(xué)圖; 系統(tǒng)誤差; 校正; 地球化學(xué)背景

以成礦區(qū)帶或以省區(qū)為單位大中比例尺區(qū)域性地球化學(xué)勘查工作的相繼開(kāi)展, 在成礦區(qū)帶獲得了海量的區(qū)域地球化學(xué)數(shù)據(jù), 對(duì)這些數(shù)據(jù)資料采用先進(jìn)的計(jì)算技術(shù)進(jìn)行開(kāi)發(fā)利用與分析研究(李隨民等, 2005), 為找尋各種金屬礦產(chǎn), 特別為找尋足以影響一個(gè)國(guó)家國(guó)計(jì)民生的特大型礦床提供全新、有效的思路與方法(當(dāng)然它也為地質(zhì)及成礦理論研究提供新思路與新線(xiàn)索)。成礦區(qū)帶的1:20萬(wàn)地球化學(xué)勘查數(shù)據(jù)都是分圖幅開(kāi)展, 不同圖幅由不同單位(或同一單位在不同時(shí)期)、不同時(shí)期完成, 這些數(shù)據(jù)也都是符合分析要求才被驗(yàn)收的, 但不同單位、不同期次完成的海量區(qū)域地球化學(xué)數(shù)據(jù)仍不可避免地存在誤差, 這種誤差被稱(chēng)為系統(tǒng)誤差, 在繪制地球化學(xué)圖時(shí)系統(tǒng)誤差的存在往往表現(xiàn)為明顯的含量等值線(xiàn)環(huán)繞子區(qū)邊界的現(xiàn)象, 由于系統(tǒng)誤差存在的系統(tǒng)性和隨機(jī)性, 在化探分析中很難被控制或者排除, 所以經(jīng)常對(duì)分析所得數(shù)據(jù)進(jìn)行處理來(lái)減少或消除系統(tǒng)誤差的影響。針對(duì)不同圖幅間的系統(tǒng)誤差先后提出了襯度法(應(yīng)祥熙, 2006)、分區(qū)標(biāo)準(zhǔn)化(紀(jì)宏金等, 1993;焦保權(quán)等, 2009)、分幅標(biāo)準(zhǔn)化(陳永清等, 1995)、平差法(陳建國(guó)等, 1997)、C型變換法(紀(jì)宏金等, 2001;陳明等, 1999)、定和化法(時(shí)艷香等, 2005)、歸一化法(劉大文, 2004)以及由分區(qū)標(biāo)準(zhǔn)化法發(fā)展而來(lái)的基于極限邊界原理和最小二乘法(紀(jì)宏金等, 2005;代永剛等, 2005)等眾多方法, 紀(jì)宏金等(2001)針對(duì)系統(tǒng)誤差校正的幾種方法的效果進(jìn)行了對(duì)比, 認(rèn)為各個(gè)方法均存在優(yōu)缺點(diǎn), 周永恒等(2011)也在上述諸多方法基礎(chǔ)上提出基于地質(zhì)背景分析的區(qū)域化探數(shù)據(jù)系統(tǒng)誤差校正, 筆者認(rèn)為該方法對(duì)于研究區(qū)面積較大、地質(zhì)背景復(fù)雜區(qū)域具有一定的局限性。紀(jì)宏金等(2001)在化探數(shù)據(jù)系統(tǒng)誤差校正、解決相鄰圖幅拼接問(wèn)題假想認(rèn)為: 某元素在不同圖幅的含量有均值不等, 均方差不等現(xiàn)象, 則校正后應(yīng)滿(mǎn)足不同圖幅均值應(yīng)相等、均方差也應(yīng)相等。基于上述原理, 本文使用了分幅(分區(qū))標(biāo)準(zhǔn)化法、襯度法、定和化法等不同方法對(duì)系統(tǒng)誤差進(jìn)行校正, 對(duì)比其效果,在此基礎(chǔ)上嘗試提出新的誤差校正方法。

研究區(qū)位于大興安嶺中段(圖1), 發(fā)育大規(guī)模中酸性侵入巖和火山巖, 局部出露中、新元古界變質(zhì)基底和古生代海相-海陸交互相火山-沉積地層, 大面積分布中生代陸相火山巖和沉積巖, 尤其是侏羅—白堊系的火山碎屑凝灰?guī)r、熔結(jié)凝灰?guī)r、礫巖、砂巖、粉砂巖、玄武巖、安山玄武巖, 玄武質(zhì)凝灰?guī)r等, 是本區(qū)最主要的地層。在研究區(qū)西部主要為侏羅系上統(tǒng)滿(mǎn)克頭鄂博組、瑪尼吐組及白音高老組,巖系為陸相火山巖系, 受北東向斷陷盆地及火山機(jī)構(gòu)控制; 東部為白堊系玄武質(zhì)凝灰?guī)r, 玄武巖、安山玄武巖等。古生代花崗巖類(lèi)(花崗閃長(zhǎng)巖、花崗巖及花崗斑巖)較發(fā)育, 主要為石炭紀(jì)—二疊紀(jì)花崗巖,中生代花崗巖以小面積分布于研究區(qū)西南部。研究區(qū)構(gòu)造主要受蒙古—鄂霍茨克板塊俯沖作用的影響,發(fā)育得爾布干等 NE─NNE向和 NW向斷裂構(gòu)造,巖漿活動(dòng)強(qiáng)烈。地球化學(xué)特征具有明顯的3個(gè)分區(qū)特征, 西部火山巖-次火山巖區(qū)的金銀鉛鋅地球化學(xué)區(qū), 中部侵入巖區(qū)的鎢鉍鉬地球化學(xué)區(qū), 東部北東向的玄武巖、安山玄武巖區(qū)的銅鈷鎳銀金多元素地球化學(xué)區(qū)(金浚等, 2011)。

1 地球化學(xué)數(shù)據(jù)的不同校正方法效果對(duì)比

研究區(qū)數(shù)據(jù)為1: 20萬(wàn)成礦區(qū)帶上不同單位、不同期次完成的區(qū)域地球化學(xué)調(diào)查數(shù)據(jù), 編制地球化學(xué)圖, 不同圖幅拼接處出現(xiàn)了明顯的含量等值線(xiàn)環(huán)繞子區(qū)邊界的現(xiàn)象(圖2), 從圖2看出, 系統(tǒng)誤差存在2個(gè)區(qū)塊, 其中B區(qū)塊由B1、B2、B3和B4四個(gè)圖幅構(gòu)成, C區(qū)塊由C1和C2兩個(gè)圖幅組成, 對(duì)各個(gè)區(qū)塊進(jìn)行地球化學(xué)參數(shù)統(tǒng)計(jì)顯示(表1), 不同區(qū)塊地球化學(xué)各個(gè)背景均值差異明顯, 針對(duì)上述現(xiàn)象分別采用分幅標(biāo)準(zhǔn)化法、襯度返回法、定和化法3種方法進(jìn)行系統(tǒng)誤差校正。

1.1 地球化學(xué)數(shù)據(jù)的分幅標(biāo)準(zhǔn)化法誤差校正

對(duì)于系統(tǒng)誤差, 假定測(cè)試值與真值(或期望值)之間具有如下關(guān)系:

其中, xij是某元素在第i個(gè)分析室對(duì)第j個(gè)樣品的測(cè)試值, xj0為第j個(gè)樣品的真值, ki為第i個(gè)分析室批量測(cè)試某元素時(shí)的系統(tǒng)偏倚系數(shù)。

因此, 系統(tǒng)偏倚系數(shù)可定義為測(cè)試值與真值之比:

設(shè)Xi, Si分別是利用第 i個(gè)分析室測(cè)試的某元素?cái)?shù)據(jù)計(jì)算的該元素的樣本均值與均方差, 則有標(biāo)準(zhǔn)化變換:

其中, xij為某元素在第i個(gè)分析室第j個(gè)樣品的測(cè)試值, tij為xij的標(biāo)準(zhǔn)化數(shù)值, m為分析室數(shù), n為樣品數(shù)。

將(1)式代人(3)式, 經(jīng)換算整理可得:

其中, X0、S0分別為真值的樣本均值與均方差。

上述則是運(yùn)用分幅標(biāo)準(zhǔn)化數(shù)據(jù)編制區(qū)域地球化學(xué)圖的統(tǒng)計(jì)學(xué)原理。

運(yùn)用該方法對(duì)研究區(qū)地球化學(xué)數(shù)據(jù)修正結(jié)果見(jiàn)表2及圖3, B區(qū)塊的四個(gè)圖幅相對(duì)于原始數(shù)據(jù)已有明顯改善, 校正后滿(mǎn)足不同圖幅間均值和均方差相等, 但是圖幅拼接處仍然具有明顯的等量線(xiàn)環(huán)繞子區(qū)邊界現(xiàn)象, C區(qū)塊兩個(gè)圖幅由原來(lái)的數(shù)據(jù)相對(duì)于全區(qū)的低背景, 修正后成為相對(duì)于全區(qū)的較高背景,圖幅拼接處也仍然存在等量線(xiàn)環(huán)繞子區(qū)邊界現(xiàn)象,認(rèn)為研究區(qū)不適合使用該誤差校正方法。

圖1 研究區(qū)地質(zhì)略圖(據(jù)陳衍景等, 2012; 內(nèi)蒙古自治區(qū)地質(zhì)礦產(chǎn)局, 1991修改)Fig. 1 Simplified geological map of the study area(modified after CHEN et al., 2012; Bureau of Geology and Mineral Resources of Inner Mongolia, 1991)

圖2 研究區(qū)誤差校正前Au地球化學(xué)圖Fig. 2 Primary geochemical map of gold element in the study area

圖3 分幅標(biāo)準(zhǔn)化法修正的Au地球化學(xué)圖Fig. 3 Corrected geochemical map of gold element after sheet standardization

1.2 地球化學(xué)數(shù)據(jù)的定和化法誤差校正

設(shè)對(duì)n個(gè)樣品分析了m個(gè)元素, 以矩陣形式將其樣品實(shí)測(cè)值記為X=(xij)n×m, 其中xij為第i個(gè)樣品第j個(gè)變量的實(shí)測(cè)值(i=1, 2, …, n; j=1, 2, …, m)。令對(duì)角陣:

如果樣品分析準(zhǔn)確, 則應(yīng)有:

這就是地球化學(xué)數(shù)據(jù)的定和特點(diǎn)。

在實(shí)際工作中, 由于受系統(tǒng)誤差影響, 式(3)未必成立。此時(shí), 若將矩陣X作如下變換:

則矩陣Y中n個(gè)樣品的各變量總和滿(mǎn)足定和特點(diǎn), 這就是地球化學(xué)數(shù)據(jù)的定和化。

當(dāng)然, 在化探實(shí)踐中, 不可能得到地殼中全部元素的分析值, 因而這種定和化處理方法也會(huì)存在誤差。但實(shí)驗(yàn)表明, 只要保證分析元素(特別是常量元素)足夠多, 就能近似滿(mǎn)足式(5)。

其中S為各個(gè)分析元素含量總和, kj為第j個(gè)元素的實(shí)測(cè)值, 則該方法可對(duì)消除系統(tǒng)誤差發(fā)揮積極作用。

表1 不同地球化學(xué)背景界限含量值及各背景均值/(ng/g)Table 1 Limit content and background average values in different geochemical backgrounds/(ng/g)

表2 不同誤差校正方法校正后均值X/(ng/g)及均方差STable 2 Means X/(ng/g) and standard deviations S in data after correction by different systematic error correction methods

在區(qū)內(nèi)27624個(gè)網(wǎng)格化樣品中, 包含了SiO2、Al2O3、K2O、Na2O、CaO、MgO、TFe、P、Ba、Ag、As、Au、Cu、Zn、B、Be、Bi、Cd、Co、Cr、F、Hg、La、Li、Mn、Mo、Nb、Ni、Pb、Sb、Sn、Sr、Th、Ti、U、V、W、Y、Zr等39個(gè)元素和組分, 根據(jù)上述原理, 對(duì)工作區(qū)數(shù)據(jù)首先計(jì)算S值, S值分布情況見(jiàn)圖 4, 進(jìn)而使用式(4)對(duì)原始數(shù)據(jù)做定和化變換, 經(jīng)過(guò)定和化處理后的 Au含量等量線(xiàn)圖(圖 5), 從地球化學(xué)含量圖可以看出等量線(xiàn)環(huán)繞子區(qū)邊界現(xiàn)象依然明顯存在, 定和化方法無(wú)法將元素空間含量分布校正到更合理狀態(tài)。

1.3 襯度返回法誤差校正

圖4 研究區(qū)樣品S值分布圖Fig. 4 Histogram of the total content of 39 elements in the study area

圖5 定和化法修正的Au元素地球化學(xué)圖Fig. 5 Corrected geochemical map of gold element after addition

將整幅圖數(shù)據(jù)按 2個(gè)臺(tái)階水平歸類(lèi), 每一個(gè)臺(tái)階求其平均值, 并求出每個(gè)點(diǎn)含量的襯度值。然后,將每點(diǎn)襯度乘上全區(qū)的平均值, 得出校正含量值。該套數(shù)據(jù)即為平抑系統(tǒng)誤差后的修正值。將“襯度返回法”調(diào)平后(圖 6)與調(diào)平前(圖 2)兩張地球化學(xué)圖進(jìn)行對(duì)比來(lái)看, 該方法誤差校正結(jié)果良好, 圖幅間明顯的等量線(xiàn)環(huán)繞子區(qū)邊界現(xiàn)象基本不存在, 尤其以 C1及 C2區(qū)塊校正后地球化學(xué)區(qū)域分帶連續(xù), 而B(niǎo)1、B2、B3和B4四個(gè)圖幅組成的區(qū)塊與其他圖幅系統(tǒng)誤差雖已明顯改善, 地球化學(xué)區(qū)域分布特征具有一定的分帶連續(xù)性, 基本不存在圖幅拼接的等量線(xiàn)環(huán)繞子區(qū)邊界現(xiàn)象, 但被校正區(qū)塊地球化學(xué)場(chǎng)的高低背景變化不明顯, 從地球化學(xué)區(qū)域分布分帶特征看, 該區(qū)塊的南部與東北部均為高背景, 區(qū)域地質(zhì)特征顯示該區(qū)塊東南部(B2、B3、B4)高背景應(yīng)該為研究區(qū)北東向高背景的一部分, 應(yīng)屬于北東向玄武巖帶對(duì)應(yīng)的地球化學(xué)場(chǎng), 說(shuō)明將該區(qū)塊簡(jiǎn)單地進(jìn)行襯度返回法校正等量線(xiàn)環(huán)繞子區(qū)邊界現(xiàn)象消失, 地球化學(xué)圖顯示具有明顯、連續(xù)的高、低地球化學(xué)分帶特征, 但是與鄰區(qū)相同地質(zhì)背景下的地球化學(xué)場(chǎng)不對(duì)應(yīng)。

2 基于地球化學(xué)背景的誤差校正

通過(guò)上述各種誤差校正方法對(duì)數(shù)據(jù)質(zhì)量的改善結(jié)果對(duì)比, 三種方法對(duì)本區(qū)的數(shù)據(jù)誤差校正均不能達(dá)到理想狀態(tài), 這可能與研究區(qū)數(shù)據(jù)存在誤差區(qū)塊的數(shù)據(jù)固有的分布特征有關(guān), B區(qū)塊數(shù)據(jù)在研究區(qū)地球化學(xué)圖上及數(shù)據(jù)分布狀態(tài)可以看出該區(qū)塊數(shù)據(jù)分布極為集中, 呈明顯的非正態(tài)分布(圖 7), 將B區(qū)塊數(shù)據(jù)的四幅數(shù)據(jù)編制地球化學(xué)圖(圖8), 從圖可以看出B4幅與其它3幅在地球化學(xué)背景上又存在明顯差異, 系統(tǒng)誤差校正時(shí)將四幅作為一個(gè)整體不合理, 因此對(duì)各個(gè)區(qū)塊分別進(jìn)行校正, 將研究區(qū)存在系統(tǒng)誤差的數(shù)據(jù)剔除編制 A區(qū)塊地球化學(xué)圖(圖9)可以看出, 地球化學(xué)圖則不存在等量線(xiàn)環(huán)繞子區(qū)邊界現(xiàn)象, 并且地球化學(xué)高、低背景分布連續(xù), 具有明顯的地球化學(xué)分帶規(guī)律, 因此可認(rèn)為A區(qū)數(shù)據(jù)為真值, 將其作為標(biāo)準(zhǔn)區(qū)塊, 將B區(qū)塊和C區(qū)塊向A區(qū)塊校正, 校正方法如下。

圖6 襯度返回法修正的Au地球化學(xué)圖Fig. 6 Corrected geochemical map of gold element after contrast deduction

圖7 B區(qū)塊校正前Au元素?cái)?shù)據(jù)分布直方圖Fig. 7 Primary histogram of gold element in block B

圖8 B區(qū)塊四幅Au地球化學(xué)圖Fig. 8 Primary geochemical map of gold element in block B

圖9 A區(qū)塊Au地球化學(xué)圖Fig. 9 Primary geochemical map of gold element in block A

按地球化學(xué)數(shù)據(jù)分布累頻的25%和75%含量值,分別以小于等于25%含量值為低背景、25%～75%含量值為背景、大于等于 75%含量為高背景, 將地球化學(xué)數(shù)據(jù)分為三個(gè)背景, 被校正的六個(gè)區(qū)塊及A區(qū)塊各背景界限含量值及各個(gè)區(qū)塊對(duì)應(yīng)各個(gè)地球化學(xué)背景的平均值見(jiàn)表 1, 從各個(gè)背景含量范圍可以看出, B區(qū)塊的各背景含量明顯高于A區(qū)塊相應(yīng)背景含量值, C區(qū)塊則明顯低于A區(qū)塊相應(yīng)含量值, 不同區(qū)塊在背景、低背景、高背景的均值又存在明顯差異。根據(jù)上述參數(shù)特征求各個(gè)地球化學(xué)背景的襯值,再乘以A區(qū)塊相應(yīng)的地球化學(xué)背景的平均值得到校正含量值, 校正后數(shù)據(jù)編制地球化學(xué)圖(圖 10), 從校正后的地球化學(xué)圖看, 存在系統(tǒng)誤差的 2個(gè)大區(qū)塊(B、C區(qū)塊)內(nèi)包括的 6個(gè)小區(qū)塊(B1、B2、B3、B4及C1、C2)內(nèi)部及與A區(qū)塊之間等量線(xiàn)環(huán)繞子區(qū)邊界現(xiàn)象消失, 地球化學(xué)圖顯示具有明顯、連續(xù)的高、低地球化學(xué)分帶特征, 對(duì)校正后數(shù)據(jù)重新計(jì)算均值與標(biāo)準(zhǔn)區(qū)域數(shù)據(jù)均值對(duì)比(表2)可以看出, 定和化方法修正后各個(gè)區(qū)塊數(shù)據(jù)均值基本沒(méi)有變化, 分幅標(biāo)準(zhǔn)化法修正后各個(gè)區(qū)塊數(shù)據(jù)均值相等為 0.9,均方差除 B2區(qū)塊略高外, 其它區(qū)塊基本相等, 襯度返回法校正后各個(gè)區(qū)塊數(shù)據(jù)均值差異較小, 但與標(biāo)準(zhǔn)區(qū)塊數(shù)據(jù)均值 0.85相比差別較大, 在襯度返回法基礎(chǔ)上改進(jìn)的地球化學(xué)背景法校正后各個(gè)區(qū)塊數(shù)據(jù)均值差異小, 并與標(biāo)準(zhǔn)區(qū)塊數(shù)據(jù)均值基本相當(dāng), 地球化學(xué)背景法修正后數(shù)據(jù)與原始數(shù)據(jù)分布散點(diǎn)圖(圖11)檢驗(yàn)顯示, 校正后數(shù)據(jù)與原始數(shù)據(jù)具有明顯的線(xiàn)性特征, 結(jié)合本區(qū)地質(zhì)背景特征顯示, B區(qū)塊校正后與其東北部和西南部地球化學(xué)高背景區(qū)形成北東向地球化學(xué)高背景帶, 該區(qū)恰為本研究區(qū)北東向分布的玄武巖區(qū)及極為發(fā)育的北東向斷裂構(gòu)造帶, 系統(tǒng)誤差校正后的地球化學(xué)背景與地質(zhì)背景的吻合, 也證明了使用該方法比較有效。

圖10 基于地球化學(xué)背景法校正后Au地球化學(xué)圖Fig. 10 Corrected geochemical map of gold element based on geochemical background

圖11 修正后Au數(shù)據(jù)與原始Au數(shù)據(jù)分布散點(diǎn)圖Fig. 11 Scatter plots of corrected data and primary data

3 結(jié)論

不同的地球化學(xué)數(shù)據(jù)誤差校正方法前提不同,具有各自的優(yōu)缺點(diǎn), 修正數(shù)據(jù)的結(jié)果不同。

(1)使用3種不同的誤差校正方法在本區(qū)修正后數(shù)據(jù)參數(shù)特征及編制地球化學(xué)圖顯示, 任何一種方法不能同時(shí)滿(mǎn)足校正后圖幅數(shù)據(jù)均值、均方差與相鄰標(biāo)準(zhǔn)區(qū)塊均值、均方差接近及編制地球化學(xué)圖基本不存在等量線(xiàn)環(huán)繞子區(qū)邊界現(xiàn)象的兩個(gè)基本條件,表明某一種校正方法可能在某一區(qū)域或者對(duì)某種分布狀態(tài)的數(shù)據(jù)是適用的、有效的。

(2)地球化學(xué)背景系統(tǒng)誤差校正法是以分幅標(biāo)準(zhǔn)化法和襯度返回法為理論基礎(chǔ), 將各個(gè)地球化學(xué)背景為單元計(jì)算修正系數(shù)進(jìn)行調(diào)平, 保證了被修正數(shù)據(jù)修正前后原有的地球化學(xué)高、低背景區(qū)域分布不發(fā)生變化, 修正后與標(biāo)準(zhǔn)數(shù)據(jù)具有同步地球化學(xué)背景及分布特征。校正結(jié)果顯示既滿(mǎn)足校正后數(shù)據(jù)均值與標(biāo)準(zhǔn)區(qū)塊數(shù)據(jù)均值基本相當(dāng), 編制地球化學(xué)圖基本不存在等量線(xiàn)環(huán)繞子區(qū)邊界現(xiàn)象, 地球化學(xué)高、低背景帶分布連續(xù), 并與原始數(shù)據(jù)表現(xiàn)為明顯的線(xiàn)性相關(guān), 與地質(zhì)背景反映的地球化學(xué)特點(diǎn)吻合,被校正區(qū)域以外的其它區(qū)域原有的地球化學(xué)場(chǎng)未發(fā)生改變。

(3)上述各種誤差校正方法對(duì)比認(rèn)為, 任何誤差校正方法應(yīng)滿(mǎn)足交界處兩邊具有相近或相同的均值和均方差, 同時(shí), 還應(yīng)滿(mǎn)足校正后的地球化學(xué)圖面上不存在等量線(xiàn)環(huán)繞子區(qū)邊界現(xiàn)象, 標(biāo)準(zhǔn)區(qū)域應(yīng)保持原有地球化學(xué)場(chǎng)不發(fā)生改變, 校正后地球化學(xué)特征應(yīng)與地質(zhì)背景所反映的地球化學(xué)場(chǎng)一致。

陳建國(guó), 郭小蘭, 劉春華. 1997. 區(qū)域地球化學(xué)圖件拼接中的圖幅平差法[J]. 地球科學(xué)─中國(guó)地質(zhì)大學(xué)學(xué)報(bào), 22(6): 619-622.

陳明, 李金春. 1999. 化探背景與異常識(shí)別的問(wèn)題與對(duì)策[J]. 地質(zhì)與勘探, 35(2): 25-29.

陳衍景, 張成, 李諾, 楊永飛, 鄧軻. 2012. 中國(guó)東北鉬礦床地質(zhì)[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 42(5): 1223-1268.

陳永清, 紀(jì)宏金. 1995.標(biāo)準(zhǔn)化地球化學(xué)圖的編制方法及應(yīng)用效果[J]. 長(zhǎng)春地質(zhì)學(xué)院學(xué)報(bào), 25(2): 216-221.

代永剛, 郝立波, 陸繼龍. 2005. 區(qū)域化探中多圖幅系統(tǒng)誤差的校正[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 35(1): 128-130.

紀(jì)宏金, 連長(zhǎng)云, 杜慶豐. 1993. 地球化學(xué)數(shù)據(jù)的標(biāo)準(zhǔn)化與襯度變換[J]. 物探化探計(jì)算技術(shù), 15(1): 19-25.

紀(jì)宏金, 林瑞慶, 周永昶. 2001. 關(guān)于若干化探數(shù)據(jù)處理方法的討論[J]. 地質(zhì)與勘探, 37(4): 55-59.

紀(jì)宏金, 時(shí)艷香, 代永剛, 郝立波, 陸繼龍. 2005. 基于邊界極限原理的多圖幅系統(tǒng)誤差校正方法[J]. 物探化探計(jì)算技術(shù), 27(2): 154-155

焦保權(quán), 白榮杰, 孫淑梅, 潘志恒, 李世平. 2009. 地球化學(xué)分區(qū)標(biāo)準(zhǔn)化方法在區(qū)域化探信息提取中的應(yīng)用[J]. 物探與化探, 33(2): 165-169.

金浚, 王京彬, 黃建軍, 陳偉民, 李占龍, 張守林, 黃儉和, 石明. 2011. 內(nèi)蒙古大興安嶺中北段成礦環(huán)境、找礦方向及勘查技術(shù)方法研究[M]. 北京: 地質(zhì)出版社: 21-24.

李隨民 姚書(shū)振. 2005. 基于 MAPGIS的分形方法確定化探異常[J]. 地球?qū)W報(bào), 26(2): 187-190.

劉大文. 2004. 區(qū)域地球化學(xué)數(shù)據(jù)的歸一化處理及應(yīng)用[J]. 物探與化探, 28(3): 273-279.

內(nèi)蒙古自治區(qū)地質(zhì)礦產(chǎn)局. 1991. 內(nèi)蒙古自治區(qū)區(qū)域地質(zhì)志[M].北京: 地質(zhì)出版社.

時(shí)艷香, 紀(jì)宏金, 陸繼龍. 2005. 地球化學(xué)數(shù)據(jù)的定和化及其在系統(tǒng)誤差校正中的應(yīng)用[J]. 物探化探計(jì)算技術(shù), 27(1): 48-50.

應(yīng)祥熙. 2006. “標(biāo)準(zhǔn)正態(tài)轉(zhuǎn)換法”用于化探數(shù)據(jù)特征的“反推”以及系統(tǒng)誤差的調(diào)平[J]. 地質(zhì)與勘探, 42(5): 84-89.

周永恒, 李守義, 張璟, 徐山. 2011. 基于地質(zhì)背景分析的區(qū)域化探數(shù)據(jù)系統(tǒng)誤差校正[J]. 吉林大學(xué)學(xué)報(bào)(地球科學(xué)版), 41(3): 753-758.

References:

Bureau of Geology and Mineral Resources of Inner Mongolia Autonomous Region. 1991. Regional Geology of Inner Mongolia Autonomous Region[M]. Beijing: Geological Publingshing House(in Chinese).

CHEN Jian-guo, GUO Xiao-lan LIU Chun-hua.1997. Map adjustment method for merging geochemical map[J]. Earth Science-Journal of China University of Geosciences, 22(6): 619-622(in Chinese with English abstract).

CHEN Ming, LI Jin-chun. 1999. Problem and countermeasure of recognition of geochemical background and anomaly[J]. Geology and Prospecting, 35(2): 25-29(in Chinese with English abstract).

CHEN Yan-jing, ZHANG Cheng, LI Nuo, YANG Yong-fei, DENG Ke. 2012. Geology of the Mo Deposits in Northeast China[J]. Journal of Jilin University(Earth Science Edition), 42(5): 1223-1268(in Chinese with English abstract).

CHEN Yong-qing, JI Hong-jin. 1995. The Compiling method of standardized geochemical map and its application effects[J]. Journal of Changchun University of Earth Sciences, 25(2): 216-221(in Chinese with English abstract).

DAI Yong-gang, HAO Li-bo, LU Ji-long. 2005. The Systemic Error Correction of Maps in Regional Geochemical Prospecting[J]. Journal of Jilin University(Earth Science Edition), 35(1): 128-130(in Chinese with English abstract).

JI Hong-jin, LIAN Chang-yun, DU Qing-feng. 1993. Standardizatiaon and contrast transformation for chemical data[J]. Computing Techniques for Geophysical and Geochemical Exploration, 15(1): 19-25(in Chinese with English abstract).

JI Hong-jin, LIN Rui-qing, ZHOU Yong-chang. 2001. A Disscussion about some data processing methods in geochemical exploration[J]. Geology and Prospecting, 37(4): 55-59(in Chinese with English abstract).

JI Hong-jin, SHI Yan-xiang, DAI Yong-gang, HAO Li-bo, LU Ji-long. 2005. A multi-map systemic error correction method based on border limiting principle[J]. Computing Techniques for Geophysical and Geochemical Exploration, 27(2): 154-155(in Chinese with English abstract).

JIAO Bao-quan, BAI Rong-jie, SUN Shu-mei, PAN Zhi-heng, LI Shi-ping. 2009. The application of geochemical zoning standardized method to the extraction of reginal geochemical information[J]. Geophysical & Geochemical Exploration, 33(2): 165-169(in Chinese with English abstract).

JIN Jun, WANG Jing-bin, HUANG Jain-jun, CHEN Wei-min, LI Zhan-long, ZHANG Shou-lin, HUANG Jian-he, SHI Ming. 2011. The metallogenic environment, prospecting direction and the prospecting technology and means study in the middle and north area of Da Hinggan Mts Inner Mongolia[M]. Beijing: Geological Publishing House: 21-24(in Chinese).

LI Sui-min, YAO Shu-zhen. 2005. The Determination of Thresholds by the Fractal Method Based on MAPGIS[J]. Acta Geoscientica Sinica, 26(2): 187-190(in Chinese with English abstract).

LIU Da-wen. 2004. The normalization of regional geochemical data and its application[J]. Geophysical & Geochemical Exploration, 28(3): 273-279(in Chinese with English abstract).

SHI Yan-xiang, JI Hong-jin, LU Ji-long. 2005. Application of geochemical data and its application systematic error correction[J]. Computing Techniques for Geophysical and Geochemical Exploration, 27(1): 48-50(in Chinese with English abstract).

YING Xiang-xi. 2006. Using “conversion method of standard normal distribution” to “deduction” of geochemical prospecting data and revise of system error[J]. Geology and Prospecting, 42(5): 84-89(in Chinese with English abstract).

ZHOU Yong-heng, LI Shou-yi, ZHANG Jing, XU Shan. 2011. The Systemic Error Correction of Regional Geochemical Mapping Data Based on Geological Background Analysis[J]. Journal of Jilin University(Earth Science Edition), 41(3): 753-758(in Chinese with English abstract).

Multi-map Systematic Error Correction Method Based on Geochemical Background: A Case Study on Gold of Regional Geochemical Survey

DAI Hui-min1, 2), ZHAO Jun2), YANG Zhong-fang1), GONG Chuan-dong3), ZHENG Chun-ying2), SUN Zhong-ren2)
1) School of Earth Sciences and Resources, China University of Geosciences(Beijing), Beijing 100083; 2) Shenyang Institute of Geology and Mineral Resources, Shenyang, Liaoning 110034; 3) Huazi Town Government of Liaoyang, Liaoyang, Liaoning 111300

To eliminate the phenomenon that content contour often circles the subregion in geochemical maps, the authors took the gold geochemical map as an example and used the division standardized method, the contrast deduction and the addition method to revise the systematic error of the regional geochemical data of the Greater Khingan Mountains metallogenic belt obtained by the new method since the geological survey, with a comparative study of the application effects. The result shows that the three methods couldn’t be applied to revising the systematic error between the analytical data of several different sections. According to the data characteristics of the study area, the geochemical background systematic error correction method is put forward based on the data division standardized method and the contrast deduction, and the four methods were compared with each other for their application effects through the data parameters, characteristics and geochemical contour maps. This means can effectively remove the systematic errors in maps-merging, and the revised data can meet the basic assumption that the systematic error in geochemical data belongs to linear error. In practical application, this method can be flexibly used for systematic error correction of chemical data so as to get better effect.

geochemical map; systematic error; correction; geochemical background

P596; O241.1

A

10.3975/cagsb.2014.05.16

本文由中國(guó)地質(zhì)調(diào)查局成礦帶區(qū)域地球化學(xué)調(diào)查示范項(xiàng)目“區(qū)域化探方法技術(shù)研究與成果集成”(編號(hào): 基[2010]礦評(píng)03-03-36)資助。

2014-04-22; 改回日期: 2014-06-27。責(zé)任編輯: 張改俠。

戴慧敏, 女, 1979年生。博士研究生, 高級(jí)工程師。主要從事應(yīng)用地球化學(xué)研究工作。通訊地址: 110034, 沈陽(yáng)市皇姑區(qū)黃河北大街280號(hào)沈陽(yáng)地質(zhì)調(diào)查中心。E-mail: daihuimin78@126.com。