-
0 引言
-
地球物理勘探可以在深部隐伏矿勘探中发挥重要作用。尤其是地质情况复杂的地区，采用多种地球物理勘探方法相结合，可以起到较好找矿效果 （程华等，2024；赵磊等，2024）。这是因为在矿石品位高低含量不同、规模大小不一、矿床不连续等情况下，往往使得矿石自身电阻率与围岩差异不大，仅仅依靠电阻率测量无法取得较好的找矿效果（谭义东，1991；刘涛等，2011；殷保全，2014；吴新刚等， 2016；吕玉增等，2023）。相比于单一的视电阻率测量方法，可以同时测量视电阻率和激电参数的激电法在大多情况下能够探测到浸染状金属硫化物引起的异常（向晓松和贾继标，2014；张晓东等，2017； 郝海强等，2023；农观海等，2023）。由于激电法受地形影响相对较小，激电异常可靠性比电阻率异常更高（陈颖祥和梁鸿喜，1977；刘晓宇，2018；吕翔， 2022；Revil et al.，2022）。
-
然而，目前大量激电测量为时间域脉冲式供电，深部目标体二次场是否能充分达到饱和而被探测感知，与充电时长密不可分（石昆法和莫焰玉， 1982；陆学村，2012）；加之需要人工多次改变 AB供电极距，以达到不同深度的探测目的（张彦和陈兴峰，2019；郭洪娟等，2022）。这一技术现状使得传统的激电测量方法具有施工难度大、效率低、成本高、深部信号微弱、探测深度有限等不足之处。
-
针对以上缺陷，陈儒军等（2020）在仪器和方法上进行了大量创新，独创性地结合扩频通信信号处理技术和微弱地球物理信号采集技术，研发出扩频激电探测技术与装备，并将这项技术命名为扩频激电（SSIP）。该项探测技术具有大深度、高精度、抗干扰、成本低、效率高等特点，在1000 m以浅的探测范围内可以测得激电参数，能有效识别有色金属矿致异常（Liu et al.，2016，2017）。该技术已在国内西藏、甘肃、云南、贵州、河南和江苏等地多个矿产勘查项目进行了有效应用（刘卫强和陈儒军，2016；陈儒军等，2020；Ahmad et al.，2024）。
-
目前，扩频激电三级测深技术在碳酸型锰矿探测应用较少，笔者以湖南省永州市零陵区东湘桥矿区锰矿勘查为例，引入扩频激电测量技术对锰矿矿床进行了勘探。经后期多个钻孔揭露，有效验证了激电异常，异常特征与矿床地下埋深、形态特征高度吻合性，为后期在相似地质背景下采用该技术方法进行地质找矿起到了较好的示范作用。因此，采用扩频激电法在祁零盆地开展碳酸锰矿勘探研究具有重要的应用价值。
-
1 地质与地球物理特征
-
1.1 矿区地质特征
-
矿区位于祁零盆地荞麦冲—东湘桥复式向斜的中南部（图1），整体构造为一向斜，地层东西两翼基本对称发育，出露地层主要有石炭系（C）、二叠系 （P）及第四系（Q4）。其中，石炭系地层为壶天群 （CPH），分布于矿区东西侧，即向斜左右两翼，主要为白云岩和白云质灰岩。二叠系地层有上统龙潭组（P3l）、中统孤峰组（P2g）、小江边组（P2x）和栖霞组 （P2q）。上统龙潭组（P3l）地层岩性以泥岩和页岩为主，位于矿区中部，亦为向斜核部地层；中统孤峰组 （P2g）—栖霞组（P2q）地层则呈平行整合接触关系，对称分布于向斜两翼。中统孤峰组（P2g）地层为一套碳酸盐岩相，为碳酸锰矿的主要赋存层位；小江边组（P2x）和栖霞组（P2q）地层主要岩性为硅质岩、硅质灰岩和白云质灰岩；第四系（Q）地层呈风化残坡积物多沿山坡及沟谷分布，一般厚0～20 m。
-
图1 研究区地质图
-
矿区构造较简单，总体为一轴向北东的东湘桥向斜，伴随零星规模较小的近东西向断裂构成了区内总体构造框架。东湘桥向斜在区域上位于荞麦冲—东湘桥复式向斜的南段，该向斜构造控制了区内锰矿的展布。
-
1.2 矿体地质特征
-
（1）地层层位对碳酸锰矿的控制
-
有利的赋矿地层存在，是区内碳酸锰矿产出的基本条件之一。区内碳酸锰矿主要赋存于二叠系孤峰组上段上部地层中，岩性为一套富含泥质和藻类生物的含泥灰岩、泥灰岩；向下一套地层为含硅质的较纯灰岩，锰含量较低，基本未见碳酸锰矿层的发育。因此，二叠系孤峰组上段地层含生物屑的不纯碳酸盐岩控制了碳酸锰矿的产出和分布。
-
（2）构造对碳酸锰矿的控制
-
构造对碳酸锰矿的控制，不仅表现在区域成矿上，而且对特定矿床的控制上也起到了重要作用。据该区域地质资料表明，在区域上祁零盆地属于成锰盆地，该盆地的发育和形成受多组深大断裂控制。在盆地内发育一组呈雁列产出的近南北向次级向斜构造，向斜构造控制了区内锰矿床的产出，同时也控制了赋锰岩系和碳酸锰矿层的分布及矿层的产状、形态。由此可见，构造对区内锰矿床具有分级控制的特点。
-
（3）沉积相对区内锰矿床的控制
-
孤峰期沉积了一套富含铁锰质的硅质岩、碳酸盐岩和泥岩，由于相邻古陆（江南古陆和雪峰古陆） 及火山活动区的物质经强烈化学分解形成氧化硅、锰、铁流入海盆，加之海水表层、海底分别繁育大量放射虫、硅质海绵，故形成锰（铁锰）矿。矿石中的菱锰矿多与藻类生物呈球粒状分布，显示出两者的密切关系。区内孤峰期为泻湖潮坪沉积环境。
-
1.3 地球物理电性特征
-
区内由浅至深不同地层主要岩性分别为第四系（Q4）含砾黏土；上二叠统龙潭组（P3l）泥岩、页岩夹砂岩；中二叠统孤峰组上段（P2g2）泥灰岩、黄铁矿、锰矿；下段（P2g1）硅质灰岩、白云质灰岩；中二叠统小江边组（P2x）硅质岩、硅质灰岩；中二叠统栖霞组（P2q）钙质灰岩夹页岩；石炭系壶天群组（CPH）白云质灰岩等。不同地层岩性电性特征如表1所示。
-
据表1 可知，含矿地层由于不同岩石电阻率差异较大，致该地层综合电阻率与上下相邻地层有一定的差异。若仅仅依靠单一视电阻率测量，则很难判断含矿地层在地下深部的埋藏特征，其找矿效果也将大受影响。但含锰灰岩和与之伴生的黄铁矿均表现为极为突出的高极化特征，致赋矿地层综合极化率达14.2%，而其余地层综合极化率均在1.4% 及以下，两者形成10倍的差异。
-
因此，含矿地层具有中阻、高极化的电性特征，为采用激电测量进行深部找矿奠定了良好的物性前提。
-
2 基本原理
-
扩频激电（SSIP）系统采集方式为无线分布式高精度阵列采集，发送频率范围为1/256～8192 Hz，根据调制阶数的不同频率数量范围为 4～256。系统通过发送机发送 M序列伪随机扩频信号，大大提高了频域分辨率及信号抗干扰能力；通过无线分布式阵列高精度接收系统实现了大深度、高效率和高精度扩频激电信号检测。在方法上，利用相对相位谱法对电磁耦合感应进行校正（王书民和雷达，2002； 陈儒军等，2004；李栋等，2019），获取激电相位响应和 Cole-Cole 模型参数；通过扩频信号发射、Robust 和相关技术压制电磁干扰，增强信噪比；通过带地形的多频点多装置组合精密反演，实现探测目标的精密成像，得到丰富的地电结构信息；最后通过多频点多参数反演结果的分析，得到勘探成果与地质体的最佳耦合，实现成矿有利区的高可信度预测。
-
3 技术参数设置与反演
-
扩频激电三极测深装置，数据采用无线分布式阵列接收。施工时，先根据剖面长度和点距计算出测点数量，然后一次性在整条剖面上布置全部测量电极和采集站，再根据设计的供电点逐点供电，每次供电时，所有通道同时测量，当所有设计供电点完成后，整条剖面测量完成。其采集的数据量为常规人工跑极方式的10倍以上，大大提高了反演精度和可靠性，实现地下目标的精细探测，具体工作布设如图2所示。
-
其中，测点距为 10 m，测量极距 MN 为 20 m，最小—最大供电极距 A0为 15～1500 m。发送组合频率为4个的5阶扩频全方波信号，基准频率1/16 Hz，供电电流1. 0～8.5 A。
-
扩频激电数据采集完成后，由 ZondRes2D 软件反演。该软件针对扩频激电观测方式进行了多次改进，可以对任意供电电极和测量电极组合的扩频激电数据进行反演。本次选择的反演方法为Occam 反演，迭代次数为10。根据视极化率的变化范围及测区岩矿石电性参数测量结果，本文将极化率的上限设置为 10%。在剖分方面，第一层厚度为 10 m，增加因子为 10%，总层数为 30 层，最大反演深度为 823 m。反演结果表明电阻率和极化率反演误差都在可接受范围内，可以将电阻率和极化率反演结果用于解释和钻孔定位。
-
4 成果推断解释
-
由于探测目标体极化率特征相比电阻率更加明显，其异常推断解释以极化率反演结果为主，电阻率反演结果为辅。
-
从极化率反演剖面图（图3）中可以看出：在平距 125 m（标高+150 m）、200 m（标高+50 m）、300～450 m 测段（标高-100 m）、600 m（标高 0 m）、725 m （标高+50 m）和800 m（标高+175 m）处，出现了多个断断续续的高极化异常体，异常强度达 10. 0%，而背景围岩平均为 1. 0%，两者相差 10 倍。异常体强烈稳定，整体上呈下凹带状特征。
-
在电阻率剖面图（图4）中的平距 100～900 m，标高-50 m 至地表区域，为连续稳定的低阻区域，-50 m 标高以下，则为高阻区域。其连续稳定的低阻区为向斜核部的上二叠统龙潭组（P3l）地层；盆状高阻基底则指示中二叠统孤峰下段—石炭系壶天群组多个高阻地层；其高—低过渡区域（500～800 Ω·m），则为含矿地层（中二叠统孤峰上段）埋深处。
-
5 钻探揭露情况
-
对T51号勘探剖面布设了ZKT5101～ZKT5106、ZK1701、ZK3401等8个钻孔，孔深在121.5～316.2 m。分别在埋深 222.5 m、284. 0 m、281.2 m、187.7 m、 85. 0 m、59.6 m、281.5 m和117.9 m处的灰岩中，揭露到厚0.9～2.8 m锰矿体，经样品测试分析（表2），其品位 9.3%～15.2% 锰矿矿脉 2 条。所揭露的矿脉较为连续稳定，其埋藏范围和形态与激电异常带基本吻合。
-
图2 扩频激电双边三极测深示意图
-
图3 T51线扩频激电极化率（a）与地质剖面图（b）
-
图4 T51线扩频激电电阻率（a）和地质剖面图（b）
-
6 结论
-
（1）区内锰矿的产出和分布与多种成矿因素有关，与成矿关系密切的主要有地层岩性、岩相和构造等因素。其中地层产状决定了矿床的大致空间分布范围，岩性组分的复杂程度影响了矿石的综合物性特征，而构造则对矿床在地下的形态起到了重塑作用。
-
（2）由于电阻率相对于围岩差异较小，与整个含矿层相比其厚度比例极小，在无法形成一定规模电阻率异常情况下，结合极化率参数，能较好的识别矿层的空间展布。因此，多参数测量是寻找锰矿矿藏的有效途径。
-
（3）经钻探验证，扩频激电法在寻找碳酸岩型锰矿矿床方面是有效的，极化率异常与矿脉揭露深度、形态轮廓十分吻合，电阻率异常较好的反应矿区地层展布信息，同时，扩频激电法工作与传统激电测量方法相比效率快、成本低、激电异常可靠，值得推广使用。
-
参考文献
-
Ahmad J， Chen R J， Ijaz A， Muhammad Y， Shahid A S， Osama A R， Farid U， Shah F， Li R. 2024. Spread Spectrum Induced Polarization （SSIP） Survey for the Qiushuwan Copper-Molybdenum Deposits in Southern Henan Province， China［J］. Minerals， 14（9）： 934-947.
-
Liu W Q， Chen R J， Cai H Z， Luo W B， André R. 2017. Correlation analysis for spread-spectrum induced-polarization signal processing in electromagnetically noisy environments［J］. Geophysics， 82（5）： 243-256.
-
Liu W Q， Chen R J， Cai H Z， Luo W B. 2016. Robust statistical methods for impulse noise suppressing of spread spectrum induced polarization data， with application to a mine site， Gansu Province， China［J］. Journal of Applied Geophysics， 135： 397-407.
-
Revil A， Pierre V， Su Z Y， Chen R J. 2022. Induced polarization as a tool to assess mineral deposits： A review［J］. Minerals， 12（5）： 571-584.
-
陈儒军， 何继善， 白宜诚， 汤井田. 2004. 多频激电相对相位谱研究［J］. 中南大学学报（自然科学版）， （1）： 106-111.
-
陈儒军， 刘春明， 何兰芳， 唐冬春， 柳卓， 杨云见. 2020. 阵列式扩频激电法关键技术研究及应用［D］. 长沙： 中南大学.
-
陈颖祥， 梁鸿喜. 1977. 起伏地形条件下激电异常的定性分析和定量计算［J］. 地质与勘探， （1）： 67-72.
-
程华， 李水平， 白德胜， 曹杰， 孙进， 谢彦军， 荆鹏， 宋永利. 2024. 河南省嵩县九仗沟—东湾金矿区深部地球物理特征与找矿预测［J］. 矿产勘查， 15（4）： 600-611.
-
郭洪娟， 聂冀强， 赵卫东. 2022. 洪山岩体多金属成矿远景区激电测深探测深度浅析［J］. 世界有色金属， （20）： 101-103.
-
郝海强， 刘志远， 杨闪， 王建国， 张明科， 谢静博. 2023. 激发极化法在西非加纳某锰矿勘查中的应用［J］. 矿产勘查， 14 （7）： 1106-1113.
-
李栋， 杨帆， 高鹏举. 2019. 相位激电法电磁耦合效应的校正［J］. 矿产勘查， 10（2）： 316-320.
-
刘涛， 陈卫， 陈伟民， 侯朝勇. 2011. 内蒙古嘎仙镍钴矿区物探找矿方法技术组合及应用［J］. 矿产勘查， 2（6）： 772-779.
-
刘卫强， 陈儒军. 2016. 稳健统计用于扩频激电数据预处理与脉冲噪声压制［J］. 地球物理学进展， 31（3）： 1332-1341.
-
刘晓宇. 2018. 激电测深中的地形改正［D］. 北京： 中国地质大学（北京）.
-
陆学村. 2012. 激电二次场衰减分布规律实验研究［D］. 邯郸： 河北工程大学.
-
吕翔. 2022. 浅谈激电中梯及激电测深物探方法在地形复杂、切割剧烈地区的实用性［J］. 世界有色金属， （15）： 226-228.
-
吕玉增， 程一鸣， 卢光辉， 姚双秋. 2023. 综合物探技术在桂西沉积型铝土矿勘查中的应用［J］. 矿产与地质， 37（5）： 1039-1045.
-
农观海， 黄凯， 邓耀辉， 卢卯， 杜蔺. 2023. 黔西北大坪子铅锌矿区电法组合异常特征及找矿潜力分析［J］. 矿物学报， 43 （6）： 903-911.
-
石昆法， 莫焰玉. 1982. 对激电二次场时间特性影响因素的实验研究［J］. 地质与勘探， （1）： 53-57.
-
谭义东. 1991. 物探找金模式的建立和研究［J］. 地质与勘探， （1）： 55.
-
王书民， 雷达. 2002. 相位激电法（偶极—偶极）单频电磁耦合校正方法［J］. 物探与化探， （1）： 57-59， 63.
-
吴新刚， 陆桂福， 杨亚斌. 2016. 铁锌矿床物探异常特征及勘查实例［J］. 物探化探计算技术， 38（3）： 347-352.
-
向晓松， 贾继标. 2014. 双频激电法在新疆金沙滩金矿勘查中的应用［J］. 矿产勘查， 5（1）： 53-57.
-
殷保全. 2014. 综合物探方法在斑岩铜矿中的应用效果［J］. 新疆有色金属， 37（5）： 20-22.
-
赵磊， 李香资， 谢国兴， 权知心， 岑朝正， 郜朋飞， 张猛. 2024. 内蒙古索家沟银多金属矿床地球物理特征及找矿前景［J］. 矿产勘查， 15（2）： 264-272.
-
张晓东， 方捷， 张定源， 谢恩才， 凌其怀， 孙瑛， 冯海林. 2017. 激发极化法在东溪金矿接替资源勘查中的应用［J］. 物探与化探， 41（3）： 445-451.
-
张彦， 陈兴峰. 2019. 五极纵轴测深和对称四极激电测深在多金属找矿中的对比研究［J］. 世界有色金属， （3）： 61-62.
-
摘要
本文在充分研究勘查区地质背景的基础上，首次采用扩频激电法在祁零盆地开展锰矿资源的探测。通过双边三极测深反演结果分析，发现了一条较连续的相对高极化异常带，异常带强度与围岩相差10倍，突出而稳定；电阻率反演结果则表现为“U”型盆状高阻基底包裹低阻的特征，推测其高极化异常带由地下矿床引起；而高阻盆底包裹低阻则指示地层受向斜构造控制，与已知的向斜构造一致。经钻探勘探，发现2条较为连续的锰矿矿脉，厚0. 9～2. 8 m，品位9. 3%～15. 2%，与黄铁矿成伴生关系，一同赋存于灰岩中。通过SSIP实验成果以及钻探验证情况，得出视极化率异常与矿脉揭露深度、形态轮廓十分吻合，视电阻率异常较好地反应矿区地层展布信息，扩频激电法在碳酸锰矿产资源勘查方面的有效性得到验证。
Abstract
On the basis of fully studying the geological background conditions of the exploration area， this article first uses spread-spectrum induced polarization （SSIP） to explore manganese ore resources in the Qiling Basin. Through the analysis of the inversion results of SSIP sounding using the forward and inversion pole-dipole array， a relatively continuous high polarization anomaly zone was discovered， with an intensity 10 times higher than that of the surrounding rock， protruding and stable； The inversion results of apparent resistivity show the characteristic of U-shaped basin shaped high resistive basement enveloping low resistive zone， suggesting that its high polarization anomaly zone is caused by underground mineral deposit； The inclusion of high resistive basin bottom and low resistive zone indicates that the strata are controlled by syncline structures， consistent with known syncline structures. Through drilling and exploration， two relatively continuous manganese ore veins were discovered， with a thickness of 0. 9 m to 2. 8 m and a grade of 9. 3% to 15. 2%. They are associated with pyrite and coexist in limestone. Through experimental results of SSIP survey and drilling verification， it was found that the polarization anomaly is very consistent with the depth and shape of the ore vein exposure， and the resistivity anomaly reflects the distribution of the mining area strata well. Therefore， the effectiveness of the SSIP method in the exploration of manganese carbonate mineral resources is verified.
