采用音頻大地電磁法間接探測深埋富集鐵礦床

席振銖，朱偉國，張道軍，張良六，馮萬杰，鄧志剛

(1. 中南大學 地球科學與信息物理學院，長沙 410083；2. 中南大學 有色金屬成礦預測教育部重點實驗室，長沙 410083；3. 長沙五維地科勘察技術有限責任公司，長沙 410205；4. 重慶市地質勘查局 川東南地質大隊，重慶 400038)

采用音頻大地電磁法間接探測深埋富集鐵礦床

席振銖1,2，朱偉國1,2，張道軍3，張良六4，馮萬杰1,2，鄧志剛3

(1. 中南大學 地球科學與信息物理學院，長沙 410083；2. 中南大學 有色金屬成礦預測教育部重點實驗室，長沙 410083；3. 長沙五維地科勘察技術有限責任公司，長沙 410205；4. 重慶市地質勘查局 川東南地質大隊，重慶 400038)

由于深埋富集鐵礦礦體的埋深超過500 m，直接運用音頻大地電磁法(AMT法)圈定深埋富集鐵礦礦體異常的縱向分辨率難以實現，為此，根據礦床成因、成礦環境以及礦體賦存空間特征等地質因素，建立地質模型，在地質模型的基礎上構建地電模型；在地電模型的基礎上，應用AMT法達到間接找礦的目的。綦江式鐵礦和寧鄉式鐵礦的勘查實踐結果表明，應用AMT法能夠確定找礦標志、含礦層次和賦礦構造等，從而間接實現探測深埋富集鐵礦床。

綦江式鐵礦；寧鄉式鐵礦；AMT法；地電模型

隨著國民經濟的快速發展，我國對鐵資源的需求日趨增加，但由于國內0～500 m鐵礦資源不足，需加大對深部鐵礦的勘查力度，以緩解我國對鐵礦資源需求的壓力。

相對于幾何電阻率測深法，音頻大地電磁法(Audio-frequency magnetotelluric method, AMT法)[1-2]具有裝置輕便、工作效率高、不受高阻層屏蔽和對低阻層有較高分辨率等優點，是探測深部鐵礦資源的有效方法[3-4]。

在礦產勘查中，AMT法主要用來劃分地層，探測褶皺斷裂構造以及直接探測金屬礦床異常[5-7]，即使面對埋深較大的目標，只要具備了相應的地球物理條件，音頻大地電磁測深法也能夠取得較好的效果。但是由于方法本身限制，隨著趨膚深度的增加，音頻大地電磁測深法的分辨率能力也將相應降低[8-9]。因此，對于埋深較大的富集鐵礦，直接運用 AMT法來圈定鐵礦異常的難度較大，尤其是在低阻覆蓋層下進行探測，降低了 AMT的探測能力，使其無法直接探測出異常的存在[10]。

針對深埋隱伏富集鐵礦床，運用 AMT法難以直接圈定鐵礦異常的情況，本文作者提出根據礦床成因、成礦環境以及礦體賦存空間特征等地質因素，建立地質模型，在地質模型的基礎上構建電性模型；在地電模型的基礎上，應用AMT法達到間接找礦的目的。

1 AMT方法原理

AMT法是通過觀測雷電、太陽風等天然音頻電磁波在大地產生的感應電磁場來研究地下電性結構的一種地球物理方法[11-12]，采用的天然場源的頻率變化范圍為100 kHz～1 Hz，通過測量地表水平電場分量Ex，與正交的水平磁場分量Hy之比，可確定地下介質的電阻率值(ρa)。其計算公式為

式中：f是頻率，Hz；，ρa是視電阻率，?·m；E是電場強度，mV/km；H是磁場強度，nT。在電磁理論中，把電磁場在大地中傳播時，其振幅衰減到初始值 1/e時的深度，定義為趨膚深度(δ)，其公式為

由此可知，趨膚深度將隨電阻率(ρ)和頻率(f )的變化而變化。測量是在和地下研究深度相對應的頻帶上進行的。一般來說，頻率較高的數據反映淺部的電性特征，頻率較低的數據反映較深的地層特征。因此，在一個寬頻帶上觀測電場和磁場信息，并由此計算出視電阻率和相位。

2 綦江式鐵礦找礦應用

2.1 構建綦江式鐵礦找礦地電模型

根據地洼學說理論[13]，調查評價區區域大地構造處于揚子陸塊四川中生代前陸盆地華鎣山碳酸鹽臺地的南延部位。按照綦江式鐵礦成礦演化地質模式(見圖1)，綦江式鐵礦屬內陸湖沼相化學沉積型鐵礦床，也可稱為大陸湖沼型碎屑巖-泥質巖-有機巖層狀赤鐵礦-菱鐵礦礦床[14]。赤鐵礦和菱鐵礦、碎屑巖以及煤層、炭質頁巖等有機巖構成了綦江式鐵礦含礦層位的低阻特征，其上覆地層主要是侏羅系的石英砂巖，下伏地層為三疊系的灰巖，都是呈相對高阻特征。特別是礦體的富集部位與煤層、炭質頁巖等有機巖厚度和品質正相關，在礦體的富集部位，煤層、炭質頁巖等有機巖厚度大、品質好，從而大大降低含礦層位的電阻率，增強了低阻異常強度，擴大了異常范圍；同時，由于礦體早期富集在洼地和后期受褶皺構造，在空間上形成局部地質體，從而引起局部低阻異常。

圖1 綦江式鐵礦成礦演化地質模式：1—由沖積中粗粒巖屑石英砂巖組成的侵蝕基地；2—河水流入洼地沉積由侵蝕工風化碎屑組成的邊灘石英砂石；3—在蛇曲河邊灘沉積之上發育的由煤、碳質頁巖組成的泥炭湖；4—河水帶入砂、泥質、鐵質進入湖泊；5—主要成鐵時期、地表和(或)地下水流攜帶粉砂質進入湖泊，經同生-成巖期后的變化形成菱鐵礦、赤鐵礦和燧石；6—粘土巖、泥質砂巖、含鐵砂巖組成鐵層頂板或與鐵層相變或缺失Fig. 1 Ore-forming geological evolution mode of Qijiang-type iron deposit

綜合上述分析，針對綦江式鐵礦 AMT法間接找礦技術線路如下：1)結合地質資料推斷出侏羅系下統珍珠沖組和三疊系上統須家河組為巨厚低阻特征的找礦標志層位；2)通過標志層位確定含礦層位和褶皺的地下空間展布；3)在低阻標志層位中，確定閉合的相對低阻異常，低阻異常就是礦體最富集的部位，從而實現間接找礦的目的。

2.2 AMT間接探測綦江式鐵礦應用

礦區地處云貴高原與重慶市南部過渡地帶，區內地形大致呈北高南低的趨勢，山脈延伸大體呈南北向，與地質構造線的方向基本一致。野外采用EH-4電導率成像儀開展 AMT測線1條，測線方向 60°，測點58個，點距50 m，電極距50 m，采用連續電磁陣列剖面的工作方式(CEMAP)，測線總長2 600 m (2 600～5 200 m)。

圖2所示為綦江縣某鐵礦AMT電阻率斷面圖，除地表風化低阻層外，地下可分為3個電性變化區，由淺到深，分別定為層①、層②、層③。斷面圖表層下伏出現的高阻層，推斷為侏羅系中統和下統組成的高阻層①，主要由含有石英砂巖高阻巖性組成；1)下伏出現的低阻層，推斷為地層侏羅系下統和三疊系上統組成的低阻層②，主要由赤鐵礦和菱鐵礦、碎屑巖以及煤層、炭質頁巖等有機巖組成；2)下伏出現的高阻層，即斷面圖底部的高阻層，推斷為地層三疊系組成的高阻層③，主要為灰巖。從斷面成果圖又可以看出，高阻層①、低阻層②和高阻層③呈現拱形特征，推斷是東溪背斜引起，背斜軸部在里程3 150 m附近。從斷面圖不難發現在3 100～4 600 m里程段存在閉合低阻異常，推測為煤層、炭質頁巖等有機巖，厚度大和品質好，為赤鐵礦和菱鐵礦富集部位。所以，3 100～4 600 m里程段是找礦最為有利的地段。ZK4101和ZK4103的鉆孔資料揭示礦體厚度在1～5 m之間，品位約在20%左右，品位較低，而ZK4105孔礦體厚度在8～12 m之間，品位約在40%左右，品位較高，可見，從ZK4101鉆孔、ZK4103鉆孔再到ZK4105鉆孔礦體逐漸富集。低阻異常雖然不能直接反映礦體，但是可以反映礦體富集早期的洼地區，間接實現找礦的目的。

3 寧鄉式鐵礦找礦應用

3.1 構建寧鄉式鐵礦找地電模型

調查區大地構造位置位于上揚子陸塊南部碳酸鹽臺地武隆凹褶束[13]。按照寧鄉式鐵礦成礦演化地質模式(見圖3[15])，該類鐵礦為沉積巖型礦床，成礦時代為晚泥盆系黃家磴期，控礦的古地理為鄂西海盆西北部海灣濱海環境，濱海相鄰濱亞相，主礦層賦存于上泥盆統黃家磴組上部鐵質碎屑巖建造，其分布嚴格受含礦地層(上泥盆統)控制[16-17]。泥盆系上統黃家磴組為主要的含赤鐵礦和菱鐵礦礦層位，其巖性主要為泥質白云巖、泥質灰巖、灰巖，砂巖、頁巖、粉砂巖等，其電性為低阻特征；其下伏地層志留系中統羅惹坪組巖性主要為頁巖、砂質頁巖，夾砂巖，整體表現為低阻特征；而其上覆地層為石炭系中統黃龍組，其巖性主要為灰巖夾白云質灰巖及硅質灰巖，表現為高阻特征。

圖2 綦江縣某鐵礦AMT法電阻率斷面深度—里程圖[10]：1—電性分界線；2—物探推測成礦帶；3—斷層破碎帶；4—侏羅紀石英砂巖、粉砂巖帶；5—侏羅紀泥巖、頁巖帶；6—三疊系石英砂巖、粉砂巖帶；7—石英砂巖；8—泥巖Fig. 2 Depth—distance plots of cross section of resistivity from AMT in Qijiang county[10]: 1—Electrical line; 2—Geophysical speculating ore belt; 3—Fault zone; 4—Jurassic quartz sandstone and siltstone; 5—Jurassic mudstone and shale; 6—Triassic quartz sandstone, siltstone; 7—Quartz sandstone; 8—Mudstone

綜合上述分析，針對寧鄉式鐵礦 AMT法間接找礦的技術線路如下：石炭系和泥盆系具有明顯的電性差異，通過界定石炭系和泥盆系的層位空間的變化形態，就能夠確定深部含礦層位褶皺形態、賦存地下空間位置，達到指導深部找礦的目的。

3.2 AMT間接探測寧鄉式鐵礦應用

調查評價區位于渝東、鄂西接壤的巫山山脈，屬中山深切割溶蝕地貌。區內碳酸鹽巖廣布，巖溶地貌發育，地勢劇烈起伏。野外開展AMT測線1條，測線方向134°，平行于地質勘探線，測點251個，點距20 m，電極距20 m，采用連續電磁陣列剖面的工作方式(CEMAP)，測線總長 5 000 m (0～5 000 m)。

圖3 寧鄉式鐵礦成礦演化地質模式[15]Fig. 3 Ore-forming geological evolution mode of Ningxiang-type iron deposit[15]

圖4 巫山縣某鐵礦AMT法電阻率斷面深度—里程圖：1—物探推測地層界線；2—斷層破碎帶；3—三疊系灰巖帶；4—二疊系泥質灰巖、頁巖帶；5—泥盆系炭質頁巖、粘土巖、砂巖、赤鐵礦帶Fig. 4 Depth—distance plots of cross section of resistivity from AMT in Wushan county: 1—Geophysical speculating stratigraphic boundaries; 2—Fault zone; 3—Triassic limestone; 4—Dyasargillaceous limestone and shale; 5—Devonian black cat, claystone,sandstone, haematite

從巫山縣某鐵礦AMT法電阻率斷面圖(見圖4)可以看出：1)根據斷面圖、地層信息結果，以及等值線電性變化，除地表風化低阻層外，地下可分為3個電性變化區，由淺到深，分別定為 L1、L2、L3。斷面圖最右端大約4 000 m到5 000 m里程出現的高阻層，推斷為三疊系下統嘉陵江組和大冶組組成的高阻層L1，主要由灰巖、白云質灰巖高阻巖性組成；①左邊大約里程2 600 m到4 000 m出現的高低阻互層，推斷為二疊系中統、上統以及石炭系中統黃龍組組成的高低阻互層L2，主要含灰色中厚層夾薄層含泥質灰巖及黑色頁巖和煤線以及淺棕色、灰色含鈣質頁巖，形成高低阻互層；②左邊深部，約0～2 000 m里程的低阻層，推斷為由泥盆系上統黃家橙組、中統云臺觀以及志留系中統羅惹坪組組成的低阻層L3，主要含灰黑色炭質頁巖或煤線、灰白、淡紅色粘土巖、砂巖(含鐵質類砂巖)、赤鐵礦(褐鐵礦)等低阻巖層組成。2)根據斷面圖視電阻率變化，可以看出高阻層L2、高低阻互層L2和低阻層L3在小號里程呈現拱形特征，推斷是賀家坪背斜引起，背斜軸部在里程1 850 m附近；同樣高阻層L1、高低阻互層L2之間的分界線在測線中間里程呈向下方向，推斷為石磙槽向斜西翼邊緣形狀。3)根據鉆孔資料，推測巫山層鐵礦主要位于高低阻互層L2底板，低阻層L3頂板，即L2與L3交界層位處，并推出此層位隨著賀家坪背斜西翼的變化方向和石磙槽向斜西翼的延伸方向向深部拓展，且賀家坪背斜東翼向小里程方向向下延伸，其核部鐵礦層標高最大；鐵礦層隨L2與L3交界層位處向大里程趨于深部，其埋深較大。

地面 ZK5101、ZK5102、ZK5104、ZK5106 和ZK5108鉆孔揭示的石炭系和泥盆系的地質界線對應AMT電阻率斷面圖反映的高低組界面，即高阻電性層L2與低阻電性層L3的接觸面；但是，ZK5110鉆孔在1 050 m的孔深位置還是二疊系的灰巖，說明賀家坪背斜東西兩翼并不對稱，西緩東陡，這一地質特征與AMT反映的結果吻合。可見，AMT可以有效確定深部褶皺構造的展布以及深埋礦體的賦存狀態。

4 結論

1) 針對綦江式鐵礦的成因和地質背景，應用AMT法首先推斷出侏羅系下統珍珠沖組和三疊系上統須家河組為巨厚低阻特征的找礦標志層位；然后通過標志層位確定含礦層位和褶皺的地下空間展布；最后在低阻標志層位中，根據有機巖與鐵礦關系確定閉合的相對低阻異常就是礦體最富集的部位，從而實現間接找礦的目的。

2) 針對寧鄉式鐵礦的地質成因和地質背景，石炭系和泥盆系具有明顯的電性差異，應用 AMT法界定石炭系和泥盆系的層位空間的變化形態，確定深部含礦層位褶皺形態、賦存地下空間位置，達到指導深部找礦的目的。

3) 在綦江式鐵礦和寧鄉式鐵礦的勘查實踐表明，在地質模型的基礎上建立電性模型，根據地電模型，應用AMT法能夠間接實現探測深埋富集鐵礦床。

REFERENCES

[1] 李金銘. 地電場與電法勘探[M]. 北京: 地質出版社, 2005:427-440.LI Jin-ming. Geoelectric field and electrical exploration [M].Beijing: Geological Publishing House, 2005: 427-440.

[2] 湯井田, 何繼善. 可控源音頻大地電磁法及其應用[M]. 長沙:中南大學出版社, 2005: 1-6.TANG Jing-tian, HE Ji-shan. Methods and applications of CSAMT [M]. Changsha: Central South University Press, 2005:1-6.

[3] 陳樂壽, 劉國棟. 大地電磁測深研究[M]. 北京: 地震出版社,1984: 1-4.CHEN Le-shou, LIU Guo-dong. Research of magnetotelluric sounding method [M]. Beijing: Seismological Publishing House,1984: 1-4.

[4] 樸化榮. 電磁測深法原理[M]. 北京: 地質出版社, 1990:14-22.PIAO Hua-rong. Electromagnetic sounding principle [M].Beijing: Geological Publishing House, 1990: 14-22.

[5] 孫 燕, 劉建明, 曾慶棟, 張作倫, 張 松, 汪在聰, 褚少雄,葉 杰, 于昌明. 綜合地球物理方法在某金多金屬礦區找礦中的應用[J]. 地球物理進展, 2010. 25(6): 2096-2101.SUN Yan, LIU Jian-ming, ZENG Qing-dong, ZHANG Zuo-lun,ZHANG Song, WANG Zai-cong, CHU Shao-xiong, YE Jie, YU Chang-ming. Application of comprehensive geophysical methods in ore prospecting in one Au polymetallic mineralizing area [J]. Progress in Geophysics, 2010, 25(6): 2096-2101.

[6] 朱光喜. 音頻大地電磁在鐵路隧道工程勘察中的應用[J]. 工程地球物理學報, 2009, 6(3): 294-298.ZHU Guang-xi. Application of audio magnetotelluric in geological engineering exploration of railway tunnel [J]. Chinese Journal of Engineering Geophysics, 2009, 6(3): 294-298.

[7] 王玉和, 王厚臣, 程久龍, 姜國慶. EH4電磁方法在斷層定位及導水性探測中的應用[J]. 山東科技大學: 自然科學版, 2010,29(4): 35-39.WANG Yu-he, WANG Hou-chen, CHENG Jiu-long, JIANG Guo-qing. Application of EH4 electromagnetic method in locating of faults and detection of its water transmissibility [J].Journal of Shandong University of Science and Technology:Natural Science, 2010, 29(4): 35-39.

[8] NABIGHIAN M N. Electromagnetic methods in applied geophysics (Volume 1): Theory [M]. Tulsa: Society of Exploration Geophysics, 1988: 5-10.

[9] 周 軍, 王成龍, 曾 珍. 音頻大地電磁測深在攀西銅礦勘探中的應用[J]. 內蒙古石油化工, 2008, 12(1): 104-106.ZHOU Jun, WANG Cheng-long, ZENG Zhen. Application of audio magnetotelluric sounding in western of Panzhihua for copper exploration [J]. Inner Mongolia Petrochemical Industry,2008, 12(1): 104-106.

[10] 席振銖, 馮萬杰, 李瑞雪, 陳興朋. 低阻覆蓋層對高頻大地電磁測深的影響[J]. 地質與勘探, 2011, 47(4): 673-678.XI Zhen-zhu, FENG Wan-jie, LI Rui-xue, CHEN Xing-peng.Effect of a low-resistivity cover on high-frequency magnetotelluric sounding [J]. Geology and Exploration, 2011,47(4): 673-678.

[11] 考夫曼A A, 凱勒G V. 頻率域與時間域電磁測深[M]. 王建謀, 譯. 北京: 地質出版社, 1987: 112-144.KAUFMAN A A, KELLER G V. Frequency and transient soundings [M]. WANG Jian-mou, transl. Beijing: Geological Publishing House, 1987: 112-144.

[12] 石應俊, 劉國棟, 吳廣耀, 王家映, 王華俊, 王光鍔, 吳乃爵.大地電磁測深法教程[M]. 北京: 地震出版社, 1985: 7-15.SHI Ying-jun, LIU Guo-dong, WU Guang-yao, WANG Jia-ying,WANG Hua-jun, WANG Guang-e, WU Nai-jue. Magnetotelluric sounding method tutorial [M]. Beijing: Seismological Publishing House, 1985: 7-15.

[13] 陳國達. 地臺活化說及其找礦意義[M]. 北京: 地質出版社,1960.CHEN Guo-da. Activation of platform theory and its prospecting significance [M]. Beijing: Geological Publishing House, 1960.

[14] 陳 龍. 重慶綦江縣篆塘角鐵礦地質特征及成因[J]. 四川地質學報, 2009, 29(2): 149-161.CHEN Long. On geological characteristics and genesis of the Zhuantangjiao Fe deposit in Qijiang county [J]. Acta Geologica Sichuan, 2009, 29(2): 149-161.

[15] 徐 喬, 左全周, 劉雄剛, 吳發明, 黃維富, 陳 芳. 重慶市巫山縣桃花-鄧家鐵礦調查評價總體設計及2010年度工作方案[R]. 重慶: 重慶市地質調查院, 2010.XU Qiao, ZUO Quan-zhou, LIU Xiong-gang, WU Fa-ming,HUANG Wei-fu, CHEN Fang. Survey and evaluation of the overall design scheme of the year 2010 of Taohua-Dengjia iron ore in Wushan county Chongqing city [R]. Chongqing:Chongqing Bureau of Geological Exploration, 2010.

[16] 趙一鳴, 畢承思. 寧鄉式沉積鐵礦床的時空分布和演化[J].礦床地質, 2000, 19(4): 350-362.ZHAO Yi-ming, BI Cheng-si. Time-space distribution and evolution of the Ningxiang type sedimentary iron deposits [J].Mineral Deposits, 2000, 19(4): 350-362.

[17] 肖 軍, 劉嚴松, 孫傳敏, 何政偉, 周慶磊. 巫山縣桃花赤鐵礦地質特征及成因探討[J]. 礦物巖石, 2009, 29(3): 69-73.XIAO Jun, LIU Yan-song, SUN Chuan-min, HE Zheng-wei,ZHOU Qing-lei. Study on the geological characteristics and genesis of Taohua hematite deposit in the Wushan county [J].Journal of Mineralogy and Petrology, 2009, 29(3): 69-73.

Indirectly exploit buried deposits of rich iron by audio-frequency magnetotelluric method

XI Zhen-zhu1,2, ZHU Wei-guo1,2, ZHANG Dao-jun3, ZHANG Liang-liu4, FENG Wan-jie1,2, DENG Zhi-gang3
(1. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China;2. Key Laboratory of Metallogenic Prediction of Nonferrous Metals, Ministry of Education,Central South University, Changsha 410083, China;3. Changsha 5D Geo-survey and Technical Co., Ltd., Changsha 410205, China;4. Geological Team of Southeastern Sichuan, Chongqing Bureau of Geological Exploration and Exploitation of Mineral Resources, Chongqing 400038, China)

When the enrichment iron ore is buried more than 500 m, it is difficult to achieve the longitudinal resolution of iron ore abnormalities by using audio-frequency magnetotelluric (AMT) method directly. For this reason, based on the geological factors of the ore genesis, metallogenic environment, ore controlling space, and so on, the geological model was made and then the electrical model was set up, which was used to achieve the aim of indirect ore prospection by AMT method. The prospecting tests results in Qijiang-type iron deposit and Ningxiang-type iron deposit prove that the AMT method can deterimine prespecting mark, ore-bearing strata, ore-bearing structure and so on. So, indirectly detect deep-buried rich iron ore.

Qijiang-type iron deposit; Ningxiang-type iron deposit; AMT method; geoelectric model

P631.3

A

1004-0609(2012)03-0928-06

2011-12-01；

2012-01-04

席振銖，教授，博士；電話：13873150690; E-mail: xizhenzhu@163.com

(編輯 龍懷中)