王久良，楊超，許曼，楊愛雪，劉德武，趙靖雯，高超

(1.河北省地質(zhì)礦產(chǎn)勘查開發(fā)局第五地質(zhì)大隊， 河北 唐山 063000； 2.河北省海洋地質(zhì)環(huán)境調(diào)查中心， 河北 唐山 063000)

0 引言

長城金礦位于冀東冷口斷裂帶東北端，中心地理坐標(biāo)為東經(jīng)118°50′06″，北緯40°09′19″，是冀東地區(qū)產(chǎn)于碳酸巖地層的代表性金礦床(孫愛群等，2002)。該區(qū)域巖石圈斷裂、煌斑巖等幔源巖漿巖發(fā)育，很可能存在Au富集的上地幔源區(qū)，與膠東西部金礦集區(qū)成礦背景相似，顯示出巨大的找礦潛力(張俊怡，2018)。自1987年發(fā)現(xiàn)以來，該礦床歷經(jīng)多次勘查，探明金屬量約12 t，平均品位5.71×10-6，但目前探明礦體基本賦存在-300 m標(biāo)高以上(河北省地質(zhì)礦產(chǎn)局第五地質(zhì)大隊，2018①)；近些年隨著勘探程度的增加，礦體連續(xù)性及品位均較差，勘查效果不理想，一味地投入鉆探驗證，勘查風(fēng)險以及勘查成本均較大，因而對長城金礦-300 m標(biāo)高以下部位進行預(yù)測以及精準(zhǔn)定位，已成為亟需解決的問題(畢炳坤等，2019)。本文通過總結(jié)分析長城金礦控礦因素，利用可控源音頻大地電磁測深(CSAMT)方法查明了礦區(qū)電阻率異常特征，發(fā)現(xiàn)礦區(qū)礦化蝕變帶與圍巖電阻率差異明顯，根據(jù)綜合研究建立了該礦床找礦模型，在-500～-600 m標(biāo)高間預(yù)測了第二礦化富集帶，為深部勘查提供了重要依據(jù)。

1 成礦地質(zhì)背景

長城金礦位于冷口斷裂帶北東端，該斷裂帶是冀東地區(qū)重要的多金屬成礦帶之一，大地構(gòu)造位置屬華北克拉通北緣，地處華北板塊與西伯利亞板塊之間，馬蘭峪復(fù)背斜與山海關(guān)隆起銜接部位，呈NW向分布于遷安冷口—青龍清河沿一帶(李紀(jì)良，1991)。該斷裂帶北東側(cè)地層主要為太古宙結(jié)晶基底及中元古界沉積巖蓋層，南西側(cè)主要為太古宙深變質(zhì)巖；另有一些燕山期中酸性巖體沿該斷裂帶北側(cè)侵入。區(qū)域上金礦床大多賦存于冷口斷裂帶北東側(cè)的沉積巖地層中，其中NW向構(gòu)造是主控礦構(gòu)造，NE與NW向斷層交匯處往往形成具有工業(yè)意義的礦體，斷裂帶沿線主要有長城金礦、小井峪金礦、清河沿金礦等多處金礦床(董建樂，2002)。

2 礦區(qū)地質(zhì)特征

圖1 長城金礦構(gòu)造位置圖(a)及地質(zhì)簡圖(b)(據(jù)王久良等，2021修改)1—第四系；2—侏羅系后城組；3—薊縣系楊莊組；4—長城系高于莊組三至四段；5—構(gòu)造碎裂巖；6—煌斑巖；7—逆斷層及編號；8—正斷層及編號；9—推測斷層及編號；10—地層產(chǎn)狀；11—礦體位置及編號；12—村莊；13—研究區(qū)

長城金礦區(qū)出露地層為長城系高于莊組，薊縣系楊莊組、霧迷山組和侏羅系后城組地層以及新生界第四系(圖1)。礦區(qū)內(nèi)受冷口斷裂帶影響，構(gòu)造發(fā)育，主要為NW、NE向斷裂。同時，區(qū)內(nèi)存在一背斜褶皺，軸向NE，翼部及軸部被斷裂構(gòu)造錯動，而金礦脈就賦存于上述構(gòu)造交匯處。肖營子巖體位于礦區(qū)北東側(cè)3 km處，而礦區(qū)未見巖漿巖出露(樊秉鴻等，1998)。

礦區(qū)共發(fā)現(xiàn)40余條金礦體，礦床規(guī)模已達(dá)中型，礦區(qū)內(nèi)大部分礦帶受斷裂或蝕變巖帶控制，呈似層狀、透鏡狀和脈狀在礦化帶中展布，其中以Ⅰ礦帶Ⅰ-9號礦體規(guī)模最大，礦體長度360 m，平均厚度3.02 m，平均品位3.33×10-6，走向110°～140°，傾向SW，傾角55°～72°。目前所探明的礦體基本賦存于-300 m標(biāo)高以上，巖石類型主要為自碎角礫巖、蝕變巖以及鐵白云石化白云巖。淺部礦體圍巖蝕變硅化、絹云母化、泥化、碳酸鹽化及鐵白云石化均有發(fā)育，其中以鐵白云石化范圍相對較廣，主要沿隱爆角礫巖的多組裂隙呈迷霧狀、網(wǎng)脈狀分布。深部蝕變以鐵白云石化、硅化最為發(fā)育，其次是絹云母化、泥化、碳酸鹽化。鐵白云石化較強處常伴隨較強硅化；礦石為灰黑色致密塊狀蝕變白云巖時，金含量往往最高，品位可達(dá)20×10-6以上，褐鐵礦化、硅化、鐵白云石化為區(qū)內(nèi)的找礦標(biāo)志(王可南，1988)。

3 巖石電性特征

本次工作從平硐及鉆孔巖芯中采取標(biāo)本，巖性較為單一，圍巖為白云巖，礦體和蝕變帶巖性有角礫巖、硅質(zhì)白云巖、鐵白云石化白云巖，對所采集的物性標(biāo)本測量了電阻率和極化率(表1)?？梢姡簵钋f組純白云巖為相對較高的電阻率以及極化率；不純白云巖(硅質(zhì)白云巖、鐵白云石化白云巖)、角礫巖都擁有較低的電阻率以及極化率。綜合已知礦體產(chǎn)出位置及礦石類型認(rèn)為，礦區(qū)礦體電性特征表現(xiàn)為中低電阻率，而極化率與圍巖差別不大。綜上所述，本區(qū)內(nèi)礦體和礦體蝕變帶電阻率明顯低于圍巖，從物性上來說CSAMT勘查具備物理前提(魏文博等，2006)。

4 CSAMT工作

4.1 CSAMT原理

CSAMT是一種電磁勘探方法，相對于MT、AMT等無源電磁勘探方法，CSAMT采用人工發(fā)射源具有信噪比高、抗干擾能力強的優(yōu)點(柳建新等，2008)。在非均勻介質(zhì)條件下，代表一定范圍地質(zhì)體的綜合反映的電阻率以視電阻率的形式來反映(王建新等，2014)。視電阻率的計算方法有多種類型，本次采用常用的卡尼亞電阻率(ρω)，計算公式為(任宏等，2017；王振亮等，2019；劉誠，2020)：

(1)

式(1)中，ω為角頻率；μ為磁導(dǎo)率；Ex、Hy為正交的電場和磁場。由于介質(zhì)對電磁波有吸收作用，電場衰減到入射時1/e電磁波所傳播的距離稱為趨膚深度δ：

(2)

式(2)中，ρ為介質(zhì)電阻率，f為電磁波頻率。根據(jù)經(jīng)驗，探測深度D與工作頻率和大地電阻率有如下關(guān)系：

(3)

通過上述推導(dǎo)可知，電磁波頻率越大，探測深度越淺，CSAMT就是利用改變電磁波這一原理達(dá)到測深目的。

4.2 CSAMT工作方法

本次可控源音頻大地電磁測深(CSAMT)工作主要布置在礦體南西側(cè)，測線共9條，方向為0°和90°共2個方向，0°方向為西側(cè)和東側(cè)1300線—1900線，90°方向為南側(cè)和北側(cè)8180線～8780線(圖2)。本次可控源音頻大地電磁測深(CSAMT)法野外工作技術(shù)參數(shù)為：供電電極距(AB)=1200 m，接收電極距(MN)=40 m，可測扇區(qū)的夾角(θ)≤30°，工作頻率范圍f=9600～48 Hz，收發(fā)距(r)≥6500 m，點距40 m。每日采集數(shù)據(jù)傳輸至電腦，對照野外記錄確認(rèn)無誤后，利用V8配套的SSMT2000軟件，經(jīng)過傅里葉轉(zhuǎn)換實現(xiàn)時間域轉(zhuǎn)換為頻率域、相位轉(zhuǎn)換、剔除跳點，利用MT2D soft軟件對數(shù)據(jù)進行平滑、靜態(tài)校正、空間濾波處理，選擇合適的層間隔進行Bostick反演得到電阻率數(shù)據(jù)(劉巖波等，2015)。加上測地高程成果，利用sufer軟件繪制單剖面圖，再由已知地質(zhì)情況與剖面成果對照，反復(fù)修正不合理處，最終獲得CSAMT工作成果(孫燕等，2010)。

表1 巖礦石物性特征統(tǒng)計表

圖2 研究區(qū)工程布置圖1—人工堆積；2—第四系；3—侏羅系后城組；4—薊縣系楊莊組未分；5—長城系高于莊組未分；6—構(gòu)造碎裂巖；7—礦化蝕變帶；8—金礦體；9—斷層及編號；10—實測地質(zhì)界線；11—地層產(chǎn)狀；12—見礦鉆孔及編號；13—未見礦鉆孔及編號；14—物探測線及編號；15—測點及編號

4.3 電阻率異常特征

4.3.1 礦區(qū)電阻率總體特征

由0°方位立體切面圖(圖3)可見，該區(qū)電阻率由南向北依次可分為藍(lán)色低阻區(qū)(ρ<100 Ω·m)、淺藍(lán)色中阻區(qū)(100≤ρ<500 Ω·m)及黃—深紅色中高阻區(qū)(ρ≥500 Ω·m)，在礦區(qū)南部從西測線1300到東端1900線，均顯示低阻異常，西向東其異常稍有減弱，低阻異常近地表等值線南西傾，越往深部逐漸變陡近直立且部分段轉(zhuǎn)為北東傾向。測線從南向北8180～8780線，低阻異常范圍逐漸減少，電阻率逐漸升高，至8780線已達(dá)中等電阻率值。根據(jù)已有資料分析，低阻區(qū)應(yīng)為區(qū)域性斷裂即冷口斷裂帶的反映，由南向北即遠(yuǎn)離該斷裂，電阻率逐漸升高。淺藍(lán)色中阻區(qū)應(yīng)為其次級構(gòu)造發(fā)育，巖石破碎，并存在多個扭曲分支現(xiàn)象，及多個異常中心，推測該區(qū)域斷裂應(yīng)存在多期次活動，使地層出現(xiàn)隆起、滑脫、傾向扭轉(zhuǎn)等現(xiàn)象，并將地層分成若干小塊。

4.3.2 典型剖面電阻率特征

由1800和1900線CSAMT反演電阻率剖面圖(圖4)可以看出，電阻率等值線的總體分布形態(tài)向南傾，局部受構(gòu)造活動的影響向北傾，形成倒轉(zhuǎn)背斜。圖中等值線值<100 Ω·m為低阻區(qū)，應(yīng)為冷口斷裂帶的反映?？拷淇跀嗔褞橄鄬Φ妥鑵^(qū)，等值線值介于100～500 Ω·m，推測應(yīng)為構(gòu)造破碎帶或地層層間裂隙發(fā)育所致。圖中的高阻異常體，根據(jù)已知鉆孔巖芯觀察為楊莊組地層的反映。通過區(qū)內(nèi)已有鉆孔資料研究發(fā)現(xiàn)，礦體在電阻率等值線圖上的分布呈現(xiàn)一定的規(guī)律，基本賦存于-300 m標(biāo)高以上，且均分布于等值線值為100～500 Ω·m之間的相對低阻帶中，位于低阻帶各等值線過渡帶上，礦體傾向SW，與等值線傾向基本一致。礦體賦存于低阻區(qū)段，與靠近冷口斷裂帶、構(gòu)造破碎、地層層間裂隙發(fā)育有關(guān)，該區(qū)段可以為成礦流體的運移和存儲提供良好的通道與空間。

圖3 0°方位立體切面圖

圖4 CSAMT反演電阻率剖面圖a—1800線；b—1900線

4.3.3 垂向電阻率異常特征

圖5 垂向電阻率等值線圖

由電阻率等值線垂向圖可知(圖5)，-300 m標(biāo)高以上鉆孔控制見礦位置為低阻帶，電阻率等值線值均在80～500 Ω·m之間。從-300 m標(biāo)高開始，高阻異常區(qū)域開始擴大，電阻率等值線值開始增加，-300～-400 m標(biāo)高之間為次高阻異常帶，相應(yīng)位置等值線值多分布于500～1500 Ω·m之間，研究已有鉆孔資料可知，多數(shù)鉆孔在該標(biāo)高區(qū)間并未見礦，只有個別見礦。通過巖芯觀察，高阻區(qū)段主要為楊莊組致密白云巖，由此可知，高阻異常帶(>500 Ω·m)不易成礦，主要是由于高阻異常帶為白云巖，巖性較致密，層間裂隙不發(fā)育，不利于礦液的運移和存儲。從電阻率剖面圖可以看出，在1800線8620～8820點號-200 m標(biāo)高處，1900線8340～8540點號-200 m標(biāo)高處存在相對低阻帶，與已知礦體賦存的低阻區(qū)類似，顯示較好的找礦前景。電阻率等值線垂向圖顯示，-500～-600 m標(biāo)高間電阻率等值線值又出現(xiàn)降低的現(xiàn)象，等值線值介于80～500 Ω·m之間，為相對低阻帶，呈交替出現(xiàn)，推斷應(yīng)為第二個礦化富集區(qū)，找礦遠(yuǎn)景巨大。為今后在該礦區(qū)投入深部鉆探工程，對該區(qū)段進行深部探礦提供了有力依據(jù)。

5 找礦模型建立及其找礦意義

長城金礦地處冷口斷裂帶北緣，本文通過前期勘查資料及本文研究成果對其成礦特征總結(jié)如下：

(1)通過前期勘查發(fā)現(xiàn)，礦體主要賦存在薊縣系楊莊組粉紅色白云巖地層中(韋龍明等，2001)。

(2)從物探電阻率異常特征看，已知礦體主要賦存在電阻率在100～500 Ω·m范圍內(nèi)。

(3)從構(gòu)造方面來看，礦體主要賦存部位主要為受冷口斷裂帶影響的次級構(gòu)造中(王成輝等，2012)。

長城金礦多年來綜合研究，其成礦類型尚未達(dá)成統(tǒng)一認(rèn)識，但多數(shù)學(xué)者認(rèn)為其成因為卡林型金礦(王寶德等，2003)，本文參照卡林型金礦低溫?zé)嵋呵秩胩攸c及礦區(qū)電阻率特異常特征，首次將長城金礦成礦規(guī)律與礦區(qū)電阻率異常相結(jié)合，建立了礦區(qū)找礦模型(圖6)。

6 結(jié)論

(1)通過綜合分析研究區(qū)地質(zhì)資料并結(jié)合該區(qū)地球物理特征，結(jié)合礦區(qū)電阻率異常特征建立了礦區(qū)找礦模型，并在-500～-600 m標(biāo)高間發(fā)現(xiàn)了第二礦化富集帶，為深部找礦指明了方向，預(yù)測該礦床深部進一步開發(fā)將產(chǎn)生重大的經(jīng)濟效益。

(2)地球物理一直是輔助深部找礦的有力手段，然而深部物性勘查數(shù)據(jù)難以直接解釋，以淺部已知條件為約束，總結(jié)地球物理找礦模型，再向深部推演，用找礦模型擬合電阻率勘查成果得出找礦方向，是較為可行的找礦方法，該方法在長城金礦的應(yīng)用，取得了較好的找礦效果，發(fā)揮了地球物理手段其高效、經(jīng)濟的特點，為今后深部找礦提供了有力的手段支撐。

圖6 電阻率特征找礦模型圖1—冷口斷裂帶；2—完整薊縣系白云巖；3—微裂隙白云巖；4—容礦構(gòu)造；5—已知礦體；6—推斷礦體；7—電阻率范圍/Ω·m

注 釋

① 河北省地質(zhì)礦產(chǎn)局第五地質(zhì)大隊. 2018. 冀東高板河-清河沿一帶中上元古界地層成礦規(guī)律研究[R].1-45.