沂南賀家溝地區化探異常特征及找礦方向研究

武斌，孫天柱*，張海瑞，劉繼勇，王秀芬

(1.山東省第四地質礦產勘查院，山東 濰坊 261021；2.山東省地礦工程勘察院，山東 濟南 250000)

0 引言

沂南賀家溝地區HS1化探綜合異常是根據1∶5萬半程幅水系沉積物測量工作圈定，研究區地球化學景觀分區屬于低山丘陵區[1]，地勢總體起伏不大。賀家溝地區距離沂南金礦區直線距離僅28km，屬于其外圍地區[2]。在研究區內開展了1∶5萬水系沉積物測量工作，圈定了主攻礦種Au，Cu等多金屬異常，分析了樣品測量數據，總結了該區的化探異常特征，明確了下一步的找礦方向。

1 地質背景

山東省沂南賀家溝地區位于鄌郚-葛溝斷裂西側，魯中隆起的東南邊緣(圖1)[3]，區域內斷裂構造發育，致使該地區地層、巖石等發生破碎，形成脆弱地帶，給后期巖漿活動和成礦熱液活動提供了有利條件[4]。研究區內出露地層相對簡單(圖2)，寒武紀李官組、朱砂洞組等，在該區東部呈NW—SE向條帶狀展布，其中佟家莊組地層與下伏新太古代蔣峪單元巖漿巖呈角度不整合接觸，與上覆李官組呈平行不整合接觸[5]。發育震旦紀土門群佟家莊組[6]，東北部為第四紀沉積物覆蓋，且向東第四系厚度逐漸增厚。

區內巖漿巖以新太古代巖漿巖為主，包括王家溝單元石英閃長巖，蔣峪單元二長花崗巖；中生代燕山晚期巖漿巖多分布在該區北部，主要為大有單元中細粒閃長巖、大朝陽單元石英閃長玢巖、于山單元二長花崗斑巖。

2 地球化學特征

2.1 元素異常下限

本次水系沉積物樣品各元素數據基本符合正態分布[7]，故可以全區確定異常下限(表1)[8]。異常下限由平均值加兩倍的標準離差求得[9]。

1—二級單元界線；2—三級單元界線；3—四級單元界線；4—五級單元界線；5—斷層及推測斷層；6—不整合界線；7—凹陷區；8—隆起區；9—研究區；10—昌邑-大店斷裂；11—安丘-莒縣斷裂；12—鄌郚-葛溝斷裂圖1 研究區大地構造位置示意圖(據山東省大地構造劃分圖，2014)

2.2 化探異常特征

由賀家溝地區水系沉積物測量綜合異常圖可以看出(圖2)，該綜合異常中Au，Ag，Cu等主要元素的異常空間分布與中生代燕山晚期侵入巖出露范圍基本一致。Au，Cu，Mo，Sb等元素異常套合性較好，關聯性較強。該綜合異常以Au為主元素，特征組合元素為Cu，Ag，As，Mo，Sb，W等；另有Mn，Bi，Sn，Hg元素異常。主元素Au在該綜合異常中異常面積最大為2.35km2,三級濃度分帶較為明顯，Au異常與北部圖幅邊界整體呈倒“凹”型，向北未封閉延伸出界，在Au元素異常中具有4個濃集中心，其中規模最大的濃集中心峰值為10.1×10-9，平均值為3.6×10-9，襯度2.13；因為Au元素的地球化學屬性較穩定，其活動性相對較低[10]，不易形成次生富集現象，這為異常檢查工作提供了思路。Ag異常由3處組成，整體與Au異常交叉嵌套，最北部異常亦未封閉，面積0.61km2，最高具有二級濃度分帶，峰值為213.3×10-9，平均值為152.0×10-9，襯度1.44。Cu異常亦由3處組成，異常形狀與Ag異常相似，且與Au，Ag異常套合性均較好，面積0.48km2，最高具有二級濃度分帶，峰值為156.0×10-6，平均值為8.9×10-6，襯度2.08。其他元素異常特征參數值見表2。

2.3 組合分類特征

元素之間的組合關系可以有效反映異常的存在，揭示區內經歷的主要地質作用[11]、地質背景、構造環境、成礦規律等地質特征，因此科學合理的分析和提取元素組合分類，對地球化學異常的解釋、成礦預測等至關重要[12]。用SPSS數據處理軟件，結合區內實際地質概況，采用因子分析法對大量的原始變量進行降維分析，從而更好地提取區內找礦信息。

表1 主要元素地球化學參數

對水系沉積物測量樣品數據進行相關性分析[13]，得出相關系數矩陣(表3)[14]。由表3可知，Au與As相關系數為0.397，呈一般相關；Ag與Au，Sn，Cu，Bi，Sb，As，Mo，W，Mn等元素相關系數集中在0.116～0.213，呈弱相關；Hg與Au，Sb元素呈弱相關；Sn與Ag，Cu，Bi，Sb，As，Mo，W相關系數為0.143～0.191，呈弱相關；Cu與Sb相關系數為0.36，呈一般相關；Bi與As，W，Mn元素相關系數為0.301～0.488，呈一般相關，另與Au，Ag，Sn，Cu，Mo也有一定相關性；Sb與Cu相關系數為0.360，呈一般相關，另與Au，Ag，Hg，Sn，Bi，As，W等呈相關性；As與Au，Bi，W，Mn相關系數為0.397～0.488，呈一般相關，另與Ag，Sn，Sb，Mo等元素也具相關性。從表3可以看出，Au，As，W，Bi，Mn元素組相關性均較為顯著。

1—第四系沉積物；2—寒武紀朱砂洞組；3—寒武紀李官組；4—震旦紀佟家莊組；5—花崗斑巖；6—石英閃長玢巖；7—閃長巖；8—二長花崗巖；9—石英閃長巖；10—角度不整合界線；11—平行不整合界線；12—斷裂；13—產狀；14—Au異常；15—一級濃度分帶；16—二級濃度分帶；17—三級濃度分帶；18—Ag異常；19—As異常；20—Bi異常；21—Cu異常；22—Hg異常；23—Mn異常；24—Mo異常；25—Sb異常；26—Sn異常；27—W異常圖2 賀家溝地區綜合異常特征簡圖(據劉繼勇等，2019)

表2 賀家溝一帶綜合異常特征參數

表3 賀家溝地區化探元素相關系數矩陣

對區內水系樣品化學分析數據利用主成分分析法作出旋轉因子三維空間載荷分布圖(圖3)[15]，并計算得到主成分特征值、貢獻率和累計貢獻率，結果如表4所示[16]，F1，F2，F3主因子累計貢獻率達51.52%，說明F1，F2，F3因子能夠代表大部分水系樣品數據信息。

圖3 旋轉因子三維空間載荷分布圖解

從圖3、表4中可以看出，所有水系樣品數據受組件1、組件2、組件3主成分因子影響，各元素在其三維空間載荷分布圖上明顯分為3個集團。

表4 主成分因子解釋總方差及元素組合特征

集團1包括Au，W，Bi，Mn，As元素，這一元素組合表現出親硫的地球化學親合性，指示層控、低溫熱液富集成礦特征；集團2包括Ag，Sn，Cu，Sb，Mo元素，這一元素組合表現出親硫、親酸性巖的地球化學親合性，指示氣化—熱液階段的富集成礦特征；集團3只有Hg元素，且在主成分因子組件1、組件2、組件3上載荷值均為負值，說明Hg元素成礦指示性較差，這與其揮發性、絡合性的地球化學特征有關。由此可以看出Au，W，Bi，Mn，As這一組合元素可以作為該地區Au的地球化學找礦指示元素。

2.4 元素分布特征

為了探討元素在不同地質體的分布特征及可能與地質成礦有關的地球化學過程，將全區不同地質單元匯水域內水系沉積物各元素的統計特征值列表(表5)[17]，并制作每個地質單元中各元素均值相對全區各元素均值比值的對數曲線圖(圖 4)[18]，可以看出元素的富集、貧化與地層、巖漿巖具明顯的相關性[19-20]。

表5 全區及各地質單元內元素含量特征值參數

圖4 各地質分區中元素含量變化標準曲線圖

單就地表巖石類型代表的找礦時限來看(暫不考慮各時期地質構造情況下)，寒武紀地層區元素分布具有一定波動特征，多數元素出現小波峰，自寒武系以后基本呈下行或者平穩態勢；第四紀地層元素分布較為平穩，波動不大；新太古代嶧山序列和傲徠山序列侵入巖幅值變化亦不大；燕山晚期的沂南序列侵入巖中大多數元素出現曲線圖中的最高波峰。由上述分析可以初步認為該區寒武紀地層與中生代燕山晚期侵入巖出露區的主攻礦種元素富集特征較為明顯。

3 結論

(1)賀家溝地區Au與As，W，Bi、Mn異常相關性明顯，而W，Bi具有親氧、親酸的親合性，推測與中酸性侵入巖具有很大關系，其中As作為探途元素意義顯著，這可能與區內構造中的褐鐵礦化有關，因表生作用下氫氧化鐵對As有一定吸附作用。異常查證工作階段在該異常區內發現了多處Au礦化，礦化體多成透鏡狀沿中生代侵入巖與其圍巖接觸帶不均勻分布，其礦化品位亦不均勻，刻槽樣Au元素的品位為(1.0～10.5)×10-6。

(2)區內主攻礦種Au的找礦方向應主要圍繞Au元素異常的濃集中心及其周圍的斷裂構造帶和中生代侵入巖及其圍巖蝕變帶展開，現階段已經取得較好成效。

(3)本文通過研究區內的化探異常的特征，從地球化學勘查角度分析區內成礦有利部位，明確了區內找礦方向，提高找礦效率，為下一步圈定找礦靶區提供了可靠的依據。