王會(huì)敏，盧金祥，周 舟，雷 靂，夏慶霖，楊 朋，龔 銀

（1.中國(guó)地質(zhì)大學(xué)（武漢）資源學(xué)院，湖北武漢 430084；2.湖北省地質(zhì)局第七地質(zhì)大隊(duì)，湖北宜昌 443100）

0 引言

識(shí)別潛在的地球化學(xué)過(guò)程對(duì)于礦產(chǎn)勘探、土地利用規(guī)劃以及自然和人為因素的環(huán)境評(píng)估等一系列 工 作 非 常 重 要（Borovec， 1996； Evansa et al.，1996； Thoms et al.， 2010； Xie et al.， 2010； Sajn et al.， 2011； Liu et al.， 2013； Wang et al.， 2014）。近年來(lái)，地球化學(xué)數(shù)據(jù)的處理和與成礦或環(huán)境污染有關(guān)的地球化學(xué)異常的識(shí)別受到越來(lái)越多的關(guān)注。經(jīng)典統(tǒng)計(jì)學(xué)分析是一種比較流行的基于頻率和高斯定律的統(tǒng)計(jì)方法，其應(yīng)用前提是數(shù)據(jù)服從正態(tài)分布，但地球化學(xué)原始數(shù)據(jù)極少能夠滿(mǎn)足這一前提（Vistelius， 1960； Reimann and Filzmoser， 2000），因而在進(jìn)行經(jīng)典統(tǒng)計(jì)學(xué)分析之前需要將地球化學(xué)數(shù)據(jù)進(jìn)行變換以服從近似正態(tài)分布（Box and Cox，1964）。Cheng et al.（1994）提 出 了 濃 度-面 積（Concentration-Area， C-A）模型，隨后在分析勘查地球化學(xué)數(shù)據(jù)方面得到普及，基于C-A 模型又建立了能譜-面積多重分形模型（Spectrum-Area， S-A）（Cheng et al.， 2000）和局部奇異性分析（Local Singularity Analysis， LSA）（Cheng， 2007），分形/多重分形模型同時(shí)考慮了地球化學(xué)模式的頻率和空間變化，因此對(duì)地球化學(xué)異常的識(shí)別效果比較明顯。目前，在GIS 支持下對(duì)地球化學(xué)數(shù)據(jù)進(jìn)行分形/多重分形建模，特別是LSA 建模是地球化學(xué)找礦的一個(gè)熱點(diǎn)研究方向（ Zuo et al.， 2009； Carranza，2010； Zuo， 2011； Zuo et al.， 2013； Huang and Zhao， 2015； Chen et al.， 2016。Zuo and Wang，2016）。成分?jǐn)?shù)據(jù)是指樣品數(shù)據(jù)變量總和為定值的具有復(fù)雜性質(zhì)的一類(lèi)數(shù)據(jù)，由此產(chǎn)生的閉合效應(yīng)會(huì)使數(shù)據(jù)之間產(chǎn)生偽相關(guān)（Aitchison， 1982）、地球化學(xué)數(shù)據(jù)的單個(gè)成分是按整個(gè)樣品的比例或百分比來(lái)測(cè)量的，因此具有成分?jǐn)?shù)據(jù)的性質(zhì)。傳統(tǒng)的多元統(tǒng)計(jì)方法，如主成分分析（Principal Component Analysis， PCA）、聚類(lèi)分析和因子分析在處理地球化學(xué)數(shù)據(jù)時(shí)難免會(huì)出現(xiàn)元素間偽相關(guān)的弊端，進(jìn)而造成結(jié)果難以解釋?zhuān)ǜ秮嘄垼?2019； 鄭超杰， 2021）。因此，加性對(duì)數(shù)比變換（Additive Log-Ratio Transformation， ALR）（Aitchison， 1982）、中心對(duì) 數(shù) 比 變 換（Centered Log-Ratio Transformation，CLR）（Aitchison， 1982）和等距對(duì)數(shù)比變換（Isometric Log-Ratio Transformation， ILR）（Egozcue et al.，2003）被用于“打開(kāi)數(shù)據(jù)”，將數(shù)據(jù)從Aitchison 空間轉(zhuǎn)換至歐式（Euclidean）空間。在這三種變換方法中，ILR 方法相比ALR 和CLR 可以有效消除成分?jǐn)?shù)據(jù)產(chǎn)生的閉合效應(yīng)，更好地解釋數(shù)據(jù)的組成性質(zhì)（Liu et al.， 2020）。

黃陵背斜作為湖北省重要的金礦集區(qū)，從20世紀(jì)60 年 代 發(fā) 現(xiàn) 以 來(lái)（向 萌 等， 2021； 周 舟 等，2021 ），找礦工作主要依托老鄉(xiāng)報(bào)礦、傳統(tǒng)水系沉積物測(cè)量方法圈定的化探異常等。盡管前人對(duì)本區(qū)的金礦開(kāi)展了富有成效的研究，但受限于已發(fā)現(xiàn)金礦的規(guī)模較小和開(kāi)采深度較淺，找礦工作一直未取得重大突破，其成礦潛力尚有待進(jìn)一步查明。本文利用ILR變換和多重分形濾波分解方法對(duì)研究區(qū)20世紀(jì)80 年代及本世紀(jì)初完成的1∶5 萬(wàn)水系沉積物樣品中的7 種元素進(jìn)行分析，目的主要是：（1）通過(guò)比較ILR 變換和Log 變換的結(jié)果，探討成分?jǐn)?shù)據(jù)在與成礦相關(guān)的地球化學(xué)異常研究中的閉合效應(yīng)的影響；（2）基于PCA 得到的結(jié)果，利用S-A 分形模型識(shí)別與金礦化相關(guān)的化探異常，結(jié)合地質(zhì)因素及以往的勘查工作情況進(jìn)行圈定成礦遠(yuǎn)景區(qū)，為該研究區(qū)的進(jìn)一步找礦工作提供方向。

1 地質(zhì)背景及礦化特征

黃陵背斜區(qū)域上位于揚(yáng)子板塊的北部，揚(yáng)子板塊和華夏板塊共同組成了華南板塊（Guo et al.，1989； Ji et al.， 2014； Han et al.， 2017）。揚(yáng)子板塊北部以秦嶺-大別-蘇魯造山帶為界與華北克拉通相鄰，西北部以龍門(mén)山斷裂與松潘-甘孜地塊相鄰（Chen et al.， 2013）（圖1a）。揚(yáng)子板塊的基底主要為古元古代和新元古代巖石，零星分布太古代巖石（Qiu et al.， 2000）。太古代和古元古代巖石主要出露在黃陵背斜區(qū)域（ Gao et al.， 2011； Guo et al.，2014），少部分?jǐn)鄶嗬m(xù)續(xù)地出露在黃土嶺雜巖、漁洞子雜巖、斗嶺雜巖和后河雜巖以及大紅山組、河口組和陽(yáng)坡組（Han et al.， 2017）。越來(lái)越多的證據(jù)表明太古代和古元古代巖石及其繼承的鋯石顆粒遍布整個(gè)揚(yáng)子克拉通（Zhang and Zheng， 2007）。

圖1 區(qū)域大地構(gòu)造簡(jiǎn)圖（a， 據(jù)Cui et al.， 2022）和黃陵背斜核部地質(zhì)礦產(chǎn)簡(jiǎn)圖（b）Fig.1 Regional tectonic map （a， modified from Cui et al.， 2022） and geological and mineral map of the core of the Huangling anticline （b）

黃陵地塊以背斜構(gòu)造的形式坐落于華南板塊的北緣，宜昌市的西北部（圖1b），地塊核部為太古-古元古代崆嶺雜巖，主要由正片麻巖、角閃巖、蛇紋巖、偏晶巖、石英巖、大理巖和稀有麻粒巖組成（馬大銓等， 1997； Gao et al.， 1999； Guo et al.， 2015）。崆嶺雜巖遭受角閃巖相區(qū)域變質(zhì)作用（江磷生等，2002； 姜繼圣， 1986； 劉玉秀， 1987； 凌文黎等，2001 ），古元古代和太古代之間呈現(xiàn)出平行不整合接觸關(guān)系（熊成云等， 2004）。太古代變質(zhì)巖漿巖厚度大于3600 m，由TTG 片麻巖、花崗片麻巖和局部保留的斜長(zhǎng)角閃巖（在TTG 片麻巖中通常以構(gòu)造透鏡體和邊界的形式出現(xiàn)）。TTG 片麻巖的原巖為T(mén)TG（3.0 ～ 2.9 Ga）（高山和張本仁， 1990； Guo et al.， 2015 ），花崗片麻巖的原巖為高鉀A 型花崗巖（2.7 ～ 2.6 Ga）（Chen et al.，2013），斜長(zhǎng)角閃巖的原巖為中基性火山巖（富公勤等， 1993）。古元古代變質(zhì)表殼巖厚度為2500 m，為典型的孔茲巖系列巖石，包括石英-硅線(xiàn)石-石榴石（狹義孔茲巖類(lèi)）、石墨片巖、變粒巖（2.4 ～ 2.1 Ga）、長(zhǎng)英質(zhì)片麻巖和大理石（姜繼圣， 1986； 李福喜和聶學(xué)武， 1987）。孔茲巖系列的原巖為陸源碎屑沉積巖（李福喜和聶學(xué)武， 1987），形成于長(zhǎng)期克拉通化時(shí)期的淺海陸棚、瀉湖和潮坪相（邱鳳等， 2015），孔茲巖系列中還大量賦存著中小型石墨礦和石榴石礦。新元古代黃陵花崗雜巖體被新元古代蓮沱地層不整合覆蓋，由一套石英閃長(zhǎng)巖-二長(zhǎng)花崗巖-花崗閃長(zhǎng)巖單元組成，年齡大約為850 ～ 819 Ma。研究區(qū)內(nèi)斷裂構(gòu)造發(fā)育，均具有多期性與繼承性活動(dòng)特點(diǎn)，以霧渡河、板倉(cāng)河斷層為代表的NW 向斷裂最為醒目，霧渡河大斷裂是區(qū)內(nèi)輝鉬礦、磁鐵礦、稀有放射性礦產(chǎn)重要的導(dǎo)礦、控礦構(gòu)造。

黃陵背斜現(xiàn)核部金礦分布廣泛，但規(guī)模比較小，目前已發(fā)現(xiàn)70多個(gè)礦（化）點(diǎn)。核北部金礦賦存于變質(zhì)TTG 巖系和表殼巖系中，具有代表性的有水月寺、巴山寺和白竹坪三個(gè)礦集區(qū)；核南部金礦主要賦存于黃陵花崗雜巖中，具有代表性的金礦有馬滑溝、拐子溝和洪武村。這些礦床大都表現(xiàn)為以下特征：（1）礦體厚度較窄，一般為0.2 ～ 1 m，但品位較高，一般為（5 ～ 50）×10-6；（2）礦體主要賦存于主斷裂兩側(cè)的次級(jí)脆韌性復(fù)合剪切帶內(nèi)，主要為NW 向；（3）礦體呈脈狀、透鏡狀、不規(guī)則塊狀、塊狀等產(chǎn)于黃鐵礦和黃銅礦中（熊成云等， 1998）。圍巖蝕變主要有黃鐵（銅）化、硅化、絹云母化和綠泥石化，其中硅化和黃鐵礦化與金礦化關(guān)系最為密切（熊成云等， 1998）。

2 數(shù)據(jù)和方法

2.1 數(shù)據(jù)

20世紀(jì)80 ～ 90年代，湖北省鄂西地質(zhì)大隊(duì)對(duì)研究區(qū)部分圖幅進(jìn)行1∶5 萬(wàn)區(qū)域礦產(chǎn)調(diào)查，水系沉積物樣品布設(shè)以4 點(diǎn)/km2～ 8 點(diǎn)/km2為基本采樣密度。以盡量控制最大受水面積為原則，對(duì)I、II 級(jí)水系進(jìn)行系統(tǒng)采樣，對(duì)III 級(jí)水系適當(dāng)控制，在水系沉積物采樣難以控制的區(qū)域，以土壤樣品加以控制。樣品由鄂西地質(zhì)大隊(duì)試驗(yàn)站及湖北省地質(zhì)局測(cè)試，測(cè)試方法為等離子光譜分析、X 射線(xiàn)熒光方法、極譜法、化學(xué)光譜法等。21 世紀(jì)初，湖北省地調(diào)院對(duì)黃陵部分地區(qū)開(kāi)展1∶5 萬(wàn)礦產(chǎn)遠(yuǎn)景調(diào)查，水系沉積物測(cè)量各項(xiàng)工作按照《1∶5 萬(wàn)地球化學(xué)普查規(guī)范》DT/T0011-91執(zhí)行，樣品由武漢綜合巖礦測(cè)試中心分析測(cè)試，測(cè)試方法為等離子體質(zhì)譜法、發(fā)射光譜法和原子熒光法。綜合兩次礦調(diào)成果，研究區(qū)共包含化探數(shù)據(jù)18857個(gè)，共同的測(cè)試元素包括Au、Ag、Cu、Pb、Zn、Mo、Sb，元素含量單位Au為×10-9，其它為×10-6。對(duì)這七種元素進(jìn)行參數(shù)統(tǒng)計(jì)，其中研究區(qū)內(nèi)元素均值除以全國(guó)范圍內(nèi)水系沉積物均值可以得到富集系數(shù)，是反映元素是否具有富集特征及成礦潛力的重要指標(biāo)，元素的標(biāo)準(zhǔn)差與平均含量的比值為變異系數(shù)。由表1 可知，研究區(qū)內(nèi)除了Sb 元素外，其余元素富集系數(shù)均大于1.00，表現(xiàn)出富集特征，有利于成礦作用的發(fā)生，富集系數(shù)小于1.00 的元素指示貧化的特征。對(duì)數(shù)變換可以在一定程度上降低異常值的影響，使數(shù)據(jù)比原始數(shù)據(jù)分布更加符合正態(tài)分布。

表1 研究區(qū)化探元素參數(shù)統(tǒng)計(jì)表Table 1 Statistical table of geochemical exploration element parameters in the study area

2.2 基于ILR變換的主成分分析方法

主成分分析（PCA）是多元統(tǒng)計(jì)分析中最常用的方法之一，它可以將幾個(gè)相關(guān)變量組合成一個(gè)變量，從而將數(shù)據(jù)集的維數(shù)根據(jù)變量的協(xié)方差或相關(guān)性降低為不相關(guān)的主成分，得到的主成分表示多維變 量 之 間 的 相 互 關(guān) 系（Jolliffe， 2002； 李 娟 等，2021； 竇小雨等， 2022 ）。

Egozcue et al.（2003）提出了等距對(duì)數(shù)比變換（ILR），其核心思想是將單行空間映射到歐式空間，使得兩個(gè)成分向量在單行空間中的距離等于經(jīng)過(guò)等距對(duì)數(shù)比變換后得到的歐式距離。定義數(shù)據(jù)集Sm= [x1，x2，……，xm]， ILR變換公式為：

式中xi> 0，i= 1， 2， …，m。，k為任意常數(shù)，Sm數(shù)據(jù)空間由m個(gè)行向量構(gòu)成。

從地質(zhì)學(xué)角度來(lái)看，ILR 變換得到的結(jié)果不能得到合理的解釋?zhuān)驗(yàn)榈玫降男伦兞繑?shù)量（m-1）比原始變量數(shù)量（m）少一個(gè)，并且與原始變量沒(méi)有直接的關(guān)系。為了更好地解釋基于ILR 變換的主成分分析結(jié)果，將基于ILR 變換后的數(shù)據(jù)進(jìn)行主成分分析之后，再用下面正交基把結(jié)果變換到CLR空間：

V= （v1， …，vm-1）是一個(gè)行向量vi的（m-1）×m的矩陣。

2.3 能譜-面積分形模型

Cheng et al.（1994）提出了一個(gè)基本的C-A 分形模型，該模型是關(guān)于濃度值(ρ)與等值線(xiàn)值ρ所包圍的面積（A）之間的關(guān)系：

式中，A(ρ)為等值線(xiàn)值ρ所圍成的面積；v是定義ρ值對(duì)應(yīng)的含量臨界值；a1和a2是分形維數(shù)，為正值，A(ρ)/ρ可以通過(guò)雙對(duì)數(shù)坐標(biāo)圖利用最小二乘法擬合直線(xiàn)斜率得到； ∝表示比例。

Cheng et al.（1999）將C-A 模型擴(kuò)展為S-A 模型，用于在頻率域表示能譜密度-面積之間的關(guān)系：

式中，S為譜能密度； A( )≥S 表示能譜密度大于值S 的區(qū)域面積；β為各向異性尺度指數(shù)；d為總體集中度的參數(shù)； ∝表示比例。

S-A 模型基于冪律函數(shù)，該函數(shù)由二維線(xiàn)性廣義尺度不變性推導(dǎo)而來(lái)，表示在二維頻率域由波數(shù)組成的集合面積與光譜能量密度高于S( )A(> S) 之間的關(guān)系。該模型基于頻域不同的各向異性尺度特性，將空間模式分解為多個(gè)分量。它已被證明是分解混合地球化學(xué)模式的有力工具（Panahi et al.，2004； Zuo and Xia， 2009； Zuo， 2011）。

3 結(jié)果和討論

3.1 成分?jǐn)?shù)據(jù)分析

Log 變換能夠在某種程度上減少異常值的影響，得到的數(shù)據(jù)相比原始數(shù)據(jù)更加對(duì)稱(chēng)，ILR 變換能夠用于打開(kāi)原始數(shù)據(jù)，因此對(duì)研究區(qū)化探數(shù)據(jù)的七種元素使用R 語(yǔ)言程序（Templ et al.， 2017）先分別進(jìn)行Log 變換和ILR 變換，隨后進(jìn)行主成分分析。圖2 為原始數(shù)據(jù)分別基于Log 變換和ILR 變換的主成分析得到的雙標(biāo)圖，PC1 和PC2 兩個(gè)主成分分別累計(jì)占據(jù)Log 變換數(shù)據(jù)和ILR 變換數(shù)據(jù)總方差的62.199%和56.925%（表2），相對(duì)應(yīng)的特征值都大于1，表明這些成分是相對(duì)重要且有意義的（Kaiser，1958）。

表2 研究區(qū)化探元素主成分分析參數(shù)統(tǒng)計(jì)表Table 2 Statistical table of principal component analysis parameters of geochemical exploration elements in the study area

ILR 變換數(shù)據(jù)的PC1 組成和對(duì)數(shù)變換數(shù)據(jù)的PC1 組成有著明顯區(qū)別（圖2）。Log 變換后的數(shù)據(jù)變量呈現(xiàn)一邊倒的趨勢(shì)，PC1 成分所有變量均為正載荷，不能夠反應(yīng)元素之間的真實(shí)關(guān)系，PC2 變量中Au 占最大的載荷，與Ag 呈負(fù)相關(guān)，與實(shí)際情況不符。由于化探數(shù)據(jù)具有成分?jǐn)?shù)據(jù)的性質(zhì)，使用傳統(tǒng)的Log 變換會(huì)導(dǎo)致偽相關(guān)的現(xiàn)象（劉向沖等，2017；Liu et al.， 2020），ILR 變換通過(guò)構(gòu)建標(biāo)準(zhǔn)正交基以達(dá)到消除閉合效應(yīng)的目的，能夠更直觀地解釋數(shù)據(jù)的組成性質(zhì)（Filzmoser et al.， 2011），其結(jié)果顯示了兩個(gè)不同的組合。第一主成分表現(xiàn)為正的Au、Ag、Cu，且Au 元素占有最大的載荷，指示與研究區(qū)金成礦相關(guān)的元素組合，相比Log 變換的主成分分析，元素組合更加清晰。對(duì)第一主成分PC1（Au-Ag-Cu）得分進(jìn)行反距離加權(quán)插值（Inverse Distance Weight， IDW）插值得到，與Au 單元素異常圖進(jìn)行對(duì)比可以看出（圖3a，b），經(jīng)過(guò)成分?jǐn)?shù)據(jù)變換的PC1 主成分得分的高值區(qū)與礦床（點(diǎn)）有更好的空間吻合關(guān)系，并且增強(qiáng)了可能由深部礦化引起的弱異常。

3.2 能譜-面積分形模型分析

地球化學(xué)景觀是巖石形成過(guò)程、礦化作用和人類(lèi)活動(dòng)相互疊加形成的混合空間模式，不同的地球化學(xué)分區(qū)對(duì)應(yīng)不同的地區(qū)化學(xué)背景，相對(duì)應(yīng)不同的巖石類(lèi)型（Zuo， 2011）。點(diǎn)數(shù)據(jù)到柵格圖的插值計(jì)算是勘查地球化學(xué)和環(huán)境地球化學(xué)最常規(guī)的一個(gè)流程，有幾種插值方法，如反距離加權(quán)（IDW）、克里金和樣條擬合，這些方法都集成在商業(yè)地理信息系統(tǒng)軟件中。在這些方法中，IDW 是目前最受歡迎的移動(dòng)平均插值方法之一，其前提假設(shè)是鄰近觀測(cè)值比遠(yuǎn)處觀測(cè)值對(duì)插值點(diǎn)的貢獻(xiàn)更大（Zuo et al.，2009）。ILR變換后PC1的柵格圖是在Arcgis軟件中利用IDW 完成的，隨后的S-A 方法是基于GeoDAS軟件完成的。傅里葉變換用來(lái)把PC1的空間模式轉(zhuǎn)換到頻率域中，得到功率譜密度和相位兩個(gè)分量（Cheng et al.， 2000）。能譜密度（S）和單元格數(shù)量（≥S）之間的關(guān)系繪制在雙對(duì)數(shù)圖中（圖4），在擬合的雙對(duì)數(shù)曲線(xiàn)中可以看出分段信息，主要分為兩段，每段的R2都大于90%，證明擬合效果較好，左邊的線(xiàn)（y= -0.7552x+ 6.799）代表異常，右邊的線(xiàn)（y= -0.9024x+ 7.3953）代表背景。

圖4 第一主成分（ILR轉(zhuǎn)換）S-A模型雙對(duì)數(shù)圖Fig.4 Log-Log plot of first principal component （ILR transform） S-A model

利用反傅里葉變換函數(shù)將這兩個(gè)頻率分量轉(zhuǎn)換回空間域，得到兩種空間模式：背景圖（圖3c）和異常圖（圖3d）。從圖3c 可以看出，PC1 背景值高的區(qū)域主要與花崗巖有關(guān)，背景值低的主要與地層有關(guān)。此外，從異常圖（圖3d）可以看出，有效抑制了復(fù)雜地質(zhì)背景的影響，更加凸顯了礦致異常，大部分已知金礦床分布在高異常區(qū)。此外，結(jié)合異常表現(xiàn)形式和地質(zhì)特征和以往勘查工作程度，本研究圈定3處異常遠(yuǎn)景區(qū)。位于板倉(cāng)河斷層西北段以南的1號(hào)Au-Ag-Cu組合異常屬于增強(qiáng)了的弱異常，與前人圈定的1個(gè)分散流一級(jí)異常、1個(gè)分散流二級(jí)異常重疊，但新異常的面積更大，也更加明顯；位于馬滑溝金礦西側(cè)的NWW向斷層西端以南的2號(hào)異常（區(qū)內(nèi)僅有1 處面積較小的重砂異常，主要發(fā)育黃涼河組和廟灣混雜巖等老基底）和位于板倉(cāng)河斷層?xùn)|南段以南的3號(hào)異常（區(qū)內(nèi)以前未圈定任何異常，發(fā)育新元古代花崗閃長(zhǎng)巖、斑狀黑云母花崗閃長(zhǎng)巖，以及閃長(zhǎng)玢巖脈，部分地段被震旦系覆蓋）屬于新發(fā)現(xiàn)的異常，其異常強(qiáng)度較低，但濃集模式卻與含礦地段相似，異常源埋深大而造成的地表異常弱，則指示了深部可能存在金礦（化）體。

4 結(jié)論

許多實(shí)例研究表明，地球化學(xué)數(shù)據(jù)一般不服從正態(tài)分布，這給經(jīng)典統(tǒng)計(jì)學(xué)的應(yīng)用帶來(lái)了問(wèn)題，因?yàn)榻y(tǒng)計(jì)的基礎(chǔ)是基于正態(tài)分布假設(shè)的。對(duì)數(shù)變換是將原始數(shù)據(jù)轉(zhuǎn)換為對(duì)數(shù)空間以獲得正態(tài)分布的常用方法之一，它可以在一定程度上減少異常值的影響。分形模型同時(shí)考慮了地球化學(xué)數(shù)據(jù)的頻率分布和空間分布，多重分形將分形拓展到頻率域當(dāng)中，在分解混合地球地球化學(xué)模式的應(yīng)用中占據(jù)重要地位。地球化學(xué)具有成分?jǐn)?shù)據(jù)的性質(zhì)，在分析前應(yīng)先進(jìn)行打開(kāi)；否則可能會(huì)得到有偏差（元素之間偽相關(guān)）的結(jié)果。本文以鄂西黃陵背斜的水系沉積物地球化學(xué)數(shù)據(jù)為例，對(duì)上述討論進(jìn)行了驗(yàn)證，并且為該地區(qū)下一步的找礦勘查工作提供方向。通過(guò)本次研究，得到以下結(jié)論：

（1）地球化學(xué)數(shù)據(jù)應(yīng)該用ILR 轉(zhuǎn)換方法打開(kāi)，ILR 轉(zhuǎn)換后的數(shù)據(jù)相比Log 轉(zhuǎn)換后的數(shù)據(jù)得到的主成分分析結(jié)果元素組合更加清晰，其中第一主成分與研究區(qū)金礦化密切相關(guān)。

（2）地球化學(xué)模式是多重過(guò)程相互疊加的結(jié)果，S-A模型既考慮了頻率分布和空間幾何關(guān)系，能夠從地球化學(xué)模式當(dāng)中分解出礦致異常，壓制復(fù)雜背景的影響。

（3）S-A 方法得到的異常圖壓制了由復(fù)雜背景造成的異常，突出了與金礦相關(guān)的異常。據(jù)異常分布、地質(zhì)特征和以往的勘查工作程度，本文圈定了3處找礦遠(yuǎn)景區(qū)，指示該研究區(qū)下一步金礦的勘查工作。

致謝：感謝湖北省地質(zhì)局第七地質(zhì)大隊(duì)和中國(guó)地質(zhì)大學(xué)（武漢）的同事和老師在野外地質(zhì)調(diào)查和數(shù)據(jù)收集、處理過(guò)程中提供的幫助。

［附中文參考文獻(xiàn)］

竇小雨， 郭小剛， 周宏， 王秦， 茍瑞， 郭東寶.2022.甘肅白銀廠銅多金屬礦田地球化學(xué)場(chǎng)特征及信息提取［J］.地質(zhì)與勘探， 58（6）：1139-1153.

富公勤， 袁海華， 李世麟.1993.黃陵斷隆北部太古界花崗巖-綠巖地體的發(fā)現(xiàn)［J］.礦物巖石， 13（1）： 5-13.

付亞龍.2019.成分?jǐn)?shù)據(jù)處理方法研究-以甘肅省尖山-平口峽地區(qū)1：20萬(wàn)巖屑地球化學(xué)數(shù)據(jù)為例［D］.西安： 長(zhǎng)安大學(xué)： 1-59.

高山， 張本仁.1990.揚(yáng)子地臺(tái)北部太古宙TTG 片麻巖的發(fā)現(xiàn)及其意義［J］.地球科學(xué)（中國(guó)地質(zhì)大學(xué)學(xué)報(bào)）， 15（6）： 675-679.

姜繼圣.1986.黃陵變質(zhì)地區(qū)的同位素地質(zhì)年代及地殼演化［J］.長(zhǎng)春地質(zhì)學(xué)院學(xué)報(bào)， （3）： 1-11.

江磷生， 周忠友， 陳鐵龍.2002.黃陵地區(qū)的幾個(gè)主要基礎(chǔ)地質(zhì)問(wèn)題［J］.湖北地礦， 16（1）： 8-13.

李福喜， 聶學(xué)武.1987.黃陵斷隆北部崆嶺群地質(zhì)時(shí)代及地層劃分［J］.湖北地質(zhì)， 1（1）： 33-46.

李娟， 王勝建， 田玉昆， 周惠， 劉策， 薛宗安.2022.基于多源數(shù)據(jù)統(tǒng)計(jì)分析的頁(yè)巖脆性定量評(píng)價(jià)方法-以鄂西地區(qū)牛蹄塘組頁(yè)巖為例［J/OL］.中國(guó)地質(zhì).https：// kns.cnki.net/kcms/ detail /11.1167.P.20220110.1840.008.html.

凌文黎， 王歆華， 程建萍.2001.揚(yáng)子北緣晉寧期望江山基性巖體的地球化學(xué)特征及其構(gòu)造背景［J］.礦物巖石地球化學(xué)通報(bào)， 20（4）：218-221.

劉向沖，王文磊，裴英茹.2017.西藏多龍礦集區(qū)水系沉積物地球化學(xué)數(shù)據(jù)定量分析與解釋?zhuān)跩］.地質(zhì)力學(xué)學(xué)報(bào)， 23（5）： 695-706.

劉玉秀.1987.宜昌黃陵背斜北部崆嶺群變質(zhì)作用的礦物共生分析［J］.湖北地質(zhì)， 1（1）： 79-84.

馬大銓?zhuān)?李志昌， 肖志發(fā).1997.鄂西崆嶺雜巖的組成、時(shí)代及地質(zhì)演化［J］.地球?qū)W報(bào)， 18（3）： 10-18.

邱鳳， 高建營(yíng)， 裴銀， 張美云.2015.黃陵背斜石墨礦地質(zhì)特征及成礦規(guī)律［J］.資源環(huán)境與工程， 29（3）： 280-285.

向萌， 胡勝華， 聶開(kāi)紅， 盧金祥， 楊朋， 周舟.2021.鄂西黃陵背斜核部金礦地球化學(xué)特征及成因探討［J］.資源環(huán)境與工程， 35（6）：787-793.

熊成云， 韋昌山， 金光富， 李文羨， 向文金.1998.鄂西黃陵背斜核部中段金礦基本特征及成礦規(guī)律［J］.華南地質(zhì)與礦產(chǎn)， （1）：32-40.

熊成云， 韋昌山， 金光富， 譚文清， 李文羨.2004.鄂西黃陵背斜地區(qū)前南華紀(jì)古構(gòu)造格架及主要地質(zhì)事件［J］.地質(zhì)力學(xué)學(xué)報(bào)，10（2）： 97-112.

楊金香.2007.黃陵地區(qū)孔茲巖系巖石學(xué)特征及變質(zhì)作用探討［J］.資源環(huán)境與工程， 21（3）： 226-231.

鄭超杰.2021.基于成分?jǐn)?shù)據(jù)及機(jī)器學(xué)習(xí)在阿舍勒地區(qū)的綜合找礦研究［D］.桂林：桂林理工大學(xué)：1-188

周舟， 牟宗玉， 雷靂， 蔣達(dá)源， 聶開(kāi)紅， 向萌， 楊朋.2021.揚(yáng)子陸塊北緣六沖坪金礦方解石Sm-Nd 等時(shí)線(xiàn)年齡及其地質(zhì)意義［J］.資源環(huán)境與工程， 35（6）： 802-806.