地質內涵法在化探礦化信息提取中的應用探討
——以山東五井-悅莊一帶為例

武斌，張海瑞，聶佩孝，徐友松，張洪濱，史國萍

(1.山東省第四地質礦產勘查院，山東 濰坊 261021；2.山東省海岸帶地質研究中心，山東 濰坊 261021；3.山東省國土測繪院，山東 濟南 250000)

0 引言

勘查地球化學找礦方法作為礦產地質調查中最主要的找礦手段，長期以來為區域地質找礦提供主要的找礦依據，對礦產資源勘查發揮了重要作用[1]。傳統方法由于化探數據處理之后，利用化探數據的高低或結構特征來圈定異常，然后對化探異常進行評價，所以異常評價存在明顯的主觀性，異常評價過程主要受評價者專業基礎和對工作區地質條件掌握程度等人為因素影響，很可能導致對高背景下的大范圍的外帶異常和可能具有找礦前景的隱伏低緩異常的遺漏[2]。深入挖掘勘查地球化學數據隱藏的成礦信息，強調成礦作用與地球化學信息的關聯性，如何定性、定量地對異常進行客觀的綜合評價，消除異常評價的多解性問題，是化探數據處理的前沿課題[3]。傳統的化探成果圖包括地球化學圖、地球化學元素異常圖、綜合異常圖、組合異常圖、異常元素剖析圖。在前人化探數據處理的基礎上，本文運用一種新的化探成果表達方式，從化探數據代表的地質意義出發，強調變量之間的相互關系[4]，利用成礦類型異常圖和成礦強度異常圖來直觀的表達和評價區內的主要礦產的成礦潛力。

1 研究區地質背景

研究區位于華北板塊魯西隆起區魯中隆起沂山-臨朐斷隆之沂山凸起與魯山-鄒平斷隆之博山凸起、魯山凸起的相毗鄰地帶[5]，中心部位為嵩山復式雜巖體(圖1)，外圍出露新太古代及中生代中酸性侵入體和寒武紀、奧陶紀地層，斷裂構造發育，巖漿活動頻繁。研究區內各時代地層發育較齊全，包括古生代長清群、九龍群、馬家溝群、月門溝群，中生代淄博群、萊陽群、青山群，新生代官莊群及新生代第四系。

1—第四系；2—白堊系；3—侏羅系；4—二疊系；5—石炭系；6—奧陶系；7—寒武系；8—蒼山序列侵入巖；9—沂南序列侵入巖；10—傲俫山序列侵入巖；11—推斷斷裂；12—實測斷裂；13—不整合接觸界線；14—嵩山褶皺穹隆圖1 研究區地質簡圖

構造較發育，斷裂構造以脆性斷裂為主，褶皺次之。脆性斷裂按展布方向劃分，主要有NW向、NE向、近EW向、近SN向4組，構成研究區“米”字型斷裂構造格架，均為燕山運動以來的產物(1)山東省地質調查院，山東省1∶5萬五井幅、悅莊幅區域地質調查報告，1992年。；區內規模較大的褶皺構造是由中生代巖體侵位于寒武紀地層形成的穹窿構造--嵩山穹窿，位于研究區中部五井斷裂以東的嵩山周圍，穹窿核部為嵩山雜巖體，四周的寒武紀地層層理產狀向四周傾斜，為中生代巖體就位時向周邊推擠圍巖所形成。

2 樣品采集及分析測試

2.1 樣品采集

研究區地球化學景觀分區為低山丘陵區，沿斷裂、裂隙等發育的溝谷切割較深，風化剝蝕強烈，多為干涸季節性河流，是粗粒碎屑主要供給源，化探樣品均為水系沉積物樣品，采樣點密度為4.32點/km2，共采集樣品3583件，其中重復樣品79件。采樣介質為代表匯水域基巖成分的細沙、粉砂等物質成分，采樣粒度采用-60目粒級；重復樣采集由不同同組人員或質量檢查人員在同點不同時采集，用來衡量采樣誤差。按《地球化學普查規范(1∶50000)》(DZ/T0011—2015)的要求，本次重復采樣數占全區總采樣數的2.3%，五井幅內38件，悅莊幅內41件，重復樣點都均勻分布在研究區內。具體樣品采集情況見表1。

表1 樣品采集情況一覽表

2.2 樣品制備

樣品加工前在小于60℃恒溫箱內進行充分烘干，并在實驗室進行細碎加工，為防止樣品污染，采用無污染的磨樣機進行加工，樣品的加工場所進行嚴格隔離；樣品細碎加工的粒度為-0.074mm(200目)。符合粒度要求的樣品重量均不少于加工前樣品重量的90%，為保證加工粒度要求，質量檢查人員每天在已加工好的樣品中隨機抽出一定的樣品用不銹鋼篩過篩，達到90%的樣品通過-0.074mm(200目)。

2.3 分析方法

根據研究區的區域地球化學特征，結合成礦地質條件及近年化探普查元素分析項目的執行情況，本次研究確定的分析元素為Au，Ag，Cu，Pb，Zn，As，Sb，Hg，W，Sn，Bi，Mo，Ni，Co，Cr，V，Ti，Mn 18種元素，元素的測試分析方法如下：Au元素為化學光譜法，Ag元素為光譜深孔電極法，As，Sb，Hg元素為原子熒光法，Bi，Cu，Pb，Zn，Mo，W，Mn，Cr，Co，Ni，Sn，Ti，V元素為電感耦合等離子體質譜法。

樣品測試通過重復樣、二級標樣等進行控制，經監控測試元素的報出率均為92.36%大于90%，二級標樣精密度合格率為98.48%大于98%，金合格率為93.54%大于90%，內檢分析、異常點檢查和密碼抽查結果的合格率為95.74%大于90%，樣品分析結果均合格可靠。

3 確定異常下限

區內各化探元素數據基本符合正態分布[6]，在此基礎上確定了研究區的異常下限(表2)，異常下限為平均值加兩倍的標準離差求得。

表2 各元素化探特征值一覽表

續表2 各元素化探特征值一覽表

本次研究以地質時代和巖性特征為基礎，劃分地質子區[7]，統計各子區元素特征值(表3)。全區劃分為4個地質子區，分別是寒武-奧陶系、第四系，新太古代基底花崗巖類，中生代侵入巖。各子區所包含的地質體巖性為:第四系主要巖性為黏土、泥質砂巖、砂巖等沉積物；寒武系包括長清群、九龍群，主要巖性為泥巖、砂巖、灰巖、白云巖等，奧陶系包括九龍群和馬家溝群，主要巖性為灰巖、白云巖等；新太古代基底花崗巖類地質子區包括新太古代傲徠山序列花崗巖、萬山莊序列變質巖；中生代侵入巖地質子區包括沂南序列和蒼山序列，其中沂南序列巖性為閃長巖、二長閃長玢巖、石英閃長玢巖，蒼山序列巖性為二長斑巖、石英閃長玢巖。

4 討論

4.1 成礦元素組合分類

成礦元素之間的組合關系，揭示區內經歷的主要成礦地質作用，可以較大程度的反映不同成礦類型的存在，亦可反映地質背景、構造環境、成礦規律等地質特征，因此科學合理的對成礦元素組合分類進行分析，并提取有用的礦化指示信息對成礦預測起著至關重要的作用[8]。本次研究運用SPSS數據處理軟件，結合區內實際地質概況和不同類型礦床的元素組合分類特征，采用R型聚類分析法(圖2)和主成分因子分析法，對大量的化探數據原始變量進行相關性分析[9]，利用各類成礦元素的組合特征來表達客觀的地質成礦類型(表3、表4)。

表4 因子解釋總方差一覽表

圖2 R型聚類分析譜系圖(平方歐式距離組間平均連接)

表3 因子載荷矩陣正交旋轉載荷一覽表

研究區內已有鐵寨、李季等比較成熟的礦產地，其金銅多金屬礦床類型以斑巖-矽卡巖接觸交代型為主[10]，故本次首先選取斑巖-矽卡巖型礦床類型的異常值進行分析研究，區內的矽卡巖型接觸交代礦床與中生代燕山期侵入雜巖關系密切，故本次主要針對斑巖-矽卡巖型和淺成低溫熱型多金屬礦開展成礦強度評價工作。首先對區內的化探數據進行元素相關性分析，運用的方法為R型聚類分析[11-13]和因子分析。其中R型聚類分析利用平方歐氏距離的標準化得分進行各元素組間平均連接，求得各類元素間的譜系圖(圖2)，然后按標準化得分小于16，得到四組元素組合類簇(元素組合大于1種)，分別為Ⅰ類簇為Ni，Cr，V，Ti，Co；Ⅱ類簇為Au，Cu；Ⅲ類簇為Pb，Zn，Ag，Mo；Ⅳ類簇為Sb，As。

因子分析利用主成分分析法對因子載荷矩陣進行正交旋轉，變量選擇取值在0.4以上的可以將其看作能夠具有明顯成因聯系的元素組合。通過計算，可以得出累計方差貢獻率為57.84%，并在旋轉矩陣中根據主要載荷元素得出因子結構式，F1成分對應元素組合為Cr，Ni，Co，V，F2成分對應元素組合為Pb，Sb，Zn，As，F3成分對應元素組合為Au，Cu，Ag，Mo，W，F4成分對應元素組合為Bi，Mn，Ti。

斑巖-矽卡巖礦床系統中各類元素與成礦巖體的親疏關系是不同的，受其地球化學性質影響，普遍存在元素分帶。W，Cu元素相對來說較親巖體[14]，一般來說在成礦巖體近端和成礦巖體中富集成礦；Au元素既可以在巖體中和巖體近端富集，又可在巖體末端富集；而Hg，Sb，As，Pb，Zn等元素則一般在巖體的遠端富集[15]；本次分析測試的元素中Cr，Ni，Co，V，Ti，Mn，Bi元素均為基性—超基性類巖體富集元素[16-17](圖3)。

圖3 斑巖-矽卡巖系統成礦元素分帶特征圖[17]

通過R型聚類分析、因子分析和各元素與地質成礦條件的對應關系[18]，本次研究按照不同礦床類型所對應的元素組合類型將18種分析測試元素分為6類，其中大于兩種元素個數的組合分為4類，具體如下：第一類，Cr，Ni，Co，V，該類為基性元素組合；第二類，Pb，Sb，Zn，As，為斑巖-矽卡巖礦床低溫遠端富集成礦元素；第三類，Au，Cu，Ag，Mo，W，該類元素為親巖體的高溫近端元素；第四類，Bi，Mn，Ti為基性元素特征組合，其組合類型根據聚類分析和因子分析結果劃分。

4.2 礦床類型異常

礦床類型異常值其本質就是利用與不同礦床或是不同地質體相關的成礦指示元素[19]和環境指示元素的富集程度進行定量評價，從而客觀的反映出不同類型礦床的成礦異常強度，其基礎是指示元素的組合分類特征及其分類的置信度。

礦床類型異常值圈定代表礦床類型成礦的異常圖主要用到以下表達式：

式中：Ad—礦床類型異常值，n—元素個數，Ci—某元素礦集系數，xi—某個采樣點中第i種元素含量，Ti—某子區內第i種元素的異常下限，m—某地質子區內采樣點個數。

n的選取原則可分以下3種情況：①該類型礦床相關的指示元素(包括成礦指示元素和環境指示元素)；②主成礦元素和只與該類型礦床相關的指示元素；③所有與該類型礦床相關元素；

斑巖-矽卡巖金銅多金屬礦床類型異常的圈定，n選取Au，Cu，Ag，Mo，W五種組合元素，最終求得該類礦床異常圖如圖4A所示，該類金銅多金屬成礦的高大全異常主要沿NW向的九山斷裂斷續分布，區內九山斷裂最北端的異常具有兩處濃集中心，分別位于鐵寨和李季，而上述兩者也正是鐵寨金銅礦床所在位置，這也揭示了該類礦床類型異常圈定模型具有較高的置信度；區內九山斷裂中段異常位于NE，NW，NWW向斷裂交會處，該異常亦有兩處濃集中心，其中一處位于寺頭村附近，另一處濃集中心位于呂匣附近，呂匣附近的濃集中心的強度值為研究區內的最高值，顯示出其具有較好的斑巖-矽卡巖成礦潛力，但在此處目前為止未有發現的金銅多金屬礦(化)點與之相對應，故可以作為后期找礦勘查的重點；區內九山斷裂南段異常位于下洋河附近，該處異常規模雖然相對較小，但其亦具有兩級濃度分帶，顯示出該處具有一定找礦潛力。

斑巖-矽卡巖鉛鋅多金屬礦床類型異常的圈定，n選取Pb，Zn，Sb，As四種組合元素，最終求得該類礦床異常圖(圖4B)，該類鉛鋅多金屬成礦的異常分為兩類，一類沿北西向九山斷裂的兩側分布，且與九山斷裂的垂向距離具有等距性特點，另一類沿區內南部的前寒武紀花崗巖類基底與古生代蓋層的不整合面呈近東西向的波狀分布，淺成低溫熱液礦化元素沿成礦結構面(不整合面)的富集沉淀是該類異常的主要成因，區內該類礦床成礦的主要尋找方向為各類深、大斷裂兩側及成礦結構面[20]，但因區域上缺少鉛鋅多金屬低溫的成熟礦床(2)山東省地質調查院，山東省1∶25萬淄博幅區域地質調查報告，2003年。，筆者認為在區內鉛鋅多金屬礦床的找礦潛力具有一定局限性。

Cr，Ni，Co，V組合元素的成礦異常作為名義上的一種“礦床類型”，其在區內的主要地質內涵為該類具有基性—超基性特征組合元素的分布強度，從其成礦類型異常圖上來看(圖4C)，其主要沿區內NE向、NW向、NWW向的斷裂分布，這類元素基性—超基性的特性揭示了區內的這幾組斷裂具有深大斷裂的可能，另外根據研究區域上礦產特征及區內的成礦類型異常強度結合后期的礦產檢查工作，初步推斷區內這類元素組合的成礦潛力不大。

Bi，Mn，Ti組合元素成礦異常值的地質含義與上述“Cr，Ni，Co，V”相同(圖4D)，在區內并未顯示有高大全異常，這類組合元素的成礦外帶異常，主要有特定地質體所引起，通過后期的研究工作發現，這類組合主要為區內寒武紀炒米店組的含鐵質結核引起[21]，從圈定的異常圖上亦可看出與寒武系地層空間分布對應性較好，異常強度較低但穩定，均屬于低緩的大范圍異常，該類異常對于成礦的指示意義不大。

1—第四系；2—白堊系；3—侏羅系；4—二疊系；5—石炭系；6—奧陶系；7—寒武系；8—蒼山序列侵入巖；9—沂南序列侵入巖；10—傲俫山序列侵入巖；11—實測斷裂；12—不整合接觸界線；13—嵩山褶皺穹??；14—金銅礦(化)體；15—砂金礦點；A組合類型為Au-Cu-W-Ag-Mo；B組合類型為Pb-Zn-Sb-As；C組合類型為Co-Ni-Cr-V；D組合類型為Bi-Mn-Ti圖4 礦床類型異常圖

4.3 成礦強度異常

對內生礦床而言，要形成一個礦集區、大型或超大型礦床必須要有多種地質作用疊加，單一的地質作用很難形成大型或超大型的礦床。而多種地質作用的疊加改造必然造成強烈的成礦物質交換和各種復雜的元素在空間上的聚集，從而導致組合元素的復雜[22]，利用組合元素的復雜程度來代表地質成礦作用的疊加復雜程度[23]，從而客觀的反映在不同空間位置的成礦強度(圖5)。

1—第四系；2—白堊系；3—侏羅系；4—二疊系；5—石炭系；6—奧陶系；7—寒武系；8—蒼山序列侵入巖；9—沂南序列侵入巖；10—傲俫山序列侵入巖；11—推斷斷裂；12—實測斷裂；13—不整合接觸界線；14—嵩山褶皺穹?。?5—金銅礦(化)體；16—砂金礦點；17—金銅成礦遠景區；18—金銅多金屬找礦靶區；19—成礦強度異常極值圖5 成礦強度異常圖

用成礦強度異常值來圈定代表成礦強度的異常圖主要用到以下表達式：

式中：Ad—成礦強度異常值，n—某子區內所有分析元素個數，Ci—某元素礦集系數(第i種元素含量值與該子區內第i種元素異常下限的比值)，xi—某個采樣點中第i種元素含量，Ti—某子區內第i種元素的異常下限，m—某地質子區內采樣點個數。n—選取原則可分以下4種情況：①代表某子區內所有已分析的元素個數；②代表某子區內成礦元素個數；③代表某子區內成礦元素+成礦指示元素的個數；④代表某子區內成礦元素+指示元素+環境指示元素的個數。

從成礦強度異常圖來看(圖5)，成礦強度的異常帶集中區主要分布在3處。其中最北部一處成礦強度異常帶位于鐵寨-李季一帶的NW向九山斷裂沂南序列侵入巖分布地帶，該處成礦強度異常帶具有3處濃集中心，極值分別為11，11和10，代表該處至少具有10種以上成礦元素富集，反映出該處經歷復雜的地質成礦作用疊加過程，具有較好的成礦潛力；另一處成礦強度異常帶分布在九山斷裂的呂匣-寺頭一帶，該處為NE向斷裂、NWW向斷裂與NW向九山斷裂3組斷裂交會部位，成礦強度異常具有兩處濃集中心，極值分別為8和7，推測該成礦強度異常帶地質成礦疊加作用稍弱于鐵寨-李季一帶，亦具有一定成礦潛力；第三處成礦強度異常帶位于區內東南部下洋河-寶畔臺-張家莊一帶，成礦強度異常呈近東西向展布，長軸延伸方向基本與前寒武花崗巖基底與古生代蓋層的不整合界線平行，下洋河濃集中心極值為9，寶畔臺濃集中心極值為7，張家莊濃集中心極值為5，反映出該處成礦強度異常帶從下洋河-寶畔臺-張家莊地質成礦疊加作用逐步減弱，但該處異常帶具有兩級濃度分帶，亦顯示出具有一定的成礦潛力。

從上述3處成礦強度異常帶的空間分布和規模可以看出均受斷裂構造活動和中生代燕山晚期巖漿活動等重要地質成礦作用疊加影響[24]，且每個異常帶內越靠近九山斷裂其濃集中心的極值越大，說明九山斷裂為區內主要的控礦斷裂。

5 結論

(1)通過賦予數據以客觀的地質內涵的方法，來圈定的成礦類型異常圖、成礦強度異常圖等帶有特定地質意義的異常圖能客觀的凸顯成礦規律性，研究區內已發現的鐵寨、李季等成熟的金銅多金屬礦床與本次化探礦化信息提取后成果表達的異常極值在空間位置上套合較好，很好的驗證了本次研究的可信度，對于其他類型的礦床亦具有實用性。

(2)礦床類型和成礦強度的評價的核心是對成礦信息的分析與把握，關鍵是如何賦予各種成礦信息以客觀的地質內涵，基礎在于詳實的野外地質調查與綜合研究的有機結合，由理論認識創新帶動找礦方法的創新最終實現找礦突破。

(3)矽卡巖型金銅多金屬成礦強度評價成果在研究區內已有多處礦床進行驗證，但是從礦床類型異常圖和成礦強度異常圖上來看，呂匣一帶斷裂交會部位斑巖-矽卡巖型礦床類型異常和成礦強度異常均有顯示，指示該處斑巖-矽卡巖型金銅多金礦的成礦潛力較大，可作為研究區內下一步找礦工作的重點。

(4)根據本次對成礦類型異常圖和成礦強度異常圖的研究，結合研究區的地質特征、成礦背景、成礦環境在研究區內圈定出金銅多金屬成礦預測區1處，預測區內根據已有礦床與成礦強度異常的相對應關系以及異常評價的置信度圈定金銅多金屬找礦靶區3處。