多极化电磁测深法在胶东大尹格庄后疃矿区深部找矿中的应用

Application of Multi-Polarization Electromagnetic Sounding Method in the Deep Prospecting of the Houtuan Mining Area of Dayingezhuang，Jiaodong

摘要

胶东大尹格庄金矿1 500 m深度以浅矿体基本查明，开展1 500~3 000 m深度的找矿探索既是矿山进一步开发的现实需要，也是新一轮找矿突破战略行动的必然要求。通过在大尹格庄后疃矿区外围实施长度为4 000 m的大深度多极化电磁测深法勘探剖面，采集1~1×10⁵ Hz天然场电磁信号，反演计算得到3 000 m以浅电阻率和磁化率剖面。结合已有地质找矿模型，根据新太古代TTG和胶东岩群变质岩低阻高磁化率和晚侏罗世花岗岩体高阻低磁化率的电磁特性，推断低阻高磁化向高阻低磁化过渡带为控矿的招平断裂带，推测断裂带倾角变缓的扩容空间为有利赋矿部位。经80ZK2101钻孔验证，在2 943.18~3 108.41 m区间发现黄铁绢英岩化蚀变带，其中3 100.06~3 101.06 m深度为金矿体，钻孔验证结果与多极化电磁测深推断成果基本吻合。

Abstract

The shallow orebody located at a depth of 1 500 meters within the Dayingezhuang gold mine in Jiaodong has been extensively delineated.Consequently，exploring the depth range of 1 500 to 3 000 meters is not only essential for providing supplementary resources to support further mine development but also possesses considerable demonstrative value for achieving new breakthroughs in ore exploration.To this end，a 4-kilometer-long multi-polarization electromagnetic sounding profile was conducted at the periphery of the Houtuan area of the Dayingezhuang gold mine.During this study，natural field electromagnetic signals ranging from 1 Hz to 100 kHz were collected.Inversion calculations produced resistivity and magnetization profiles extending to a depth of 3 000 meters.Based on the existing geological prospecting model and the electromagnetic characteristics of Neoproterozoic TTG and Jiaodong Group metamorphic rocks，which exhibit low resistivity and high magnetization，as well as Late Jurassic granite bodies，characterized by high resistivity and low magnetization，the transition zone from low resistance and high magnetization to high resistance and low magnetization is identified as the ore-controlling Zhaoping fault zone，and the expanding space where the dip angle of the fault zone from low-resistivity and high-magnetization to high-resistivity and low-magnetization is interpreted as the ore-controlling Zhaoping fault zone.Additionally，the region where the dip angle of the fault zone becomes more horizontal is inferred to be a favorable site for mineralization.This interpretation is supported by data from drill hole 80ZK2101，which identified a pyrite-silicified alteration zone between depths of 2 943.18 m and 3 108.41 m，with a gold ore body located at depths of 3 100.06 m to 3 101.06 m.The drilling results corroborate the interpretations derived from multi-polarization electromagnetic sounding data.

Graphical abstract

关键词

多极化电磁测深法 / 深部找矿 / 准静态磁化率测深 / 大尹格庄后疃矿区 / 胶东 / 招平断裂带

Key words

multi-polarize electromagnetics sounding method / deep prospecting / quasi-static magnetic explo-ration / Dayingezhuang Houtuan mining area / Jiaodong / Zhaoping fault belt

大尹格庄金矿是分布在胶东地区招远—平度断裂带的特大型金矿床之一。近10年来，研究人员在该矿区开展了多种方法的地球物理探测工作（高帅等，2017；刘向东，2022），深部找矿取得了重要进展。大地电磁测深（MT）法因具有探测深度大、电性差异敏感等优点，在复杂地形和覆盖层条件下能够探测从浅表部至深部的电性结构信息，被广泛应用于深部找矿。研究表明，大地电磁测深及其衍生方法在大尹格庄矿区深部探测中取得了显著成效。研究人员通过大地电磁测深方法实现了-1 000 m以下深部矿体的定位，为区域深部矿体探测提供了技术支撑（邓湘伟等，2011）；在此基础上，结合构造地球化学测量揭示了1 150~1 450 m深度范围内隐伏矿体分布特征，明确了矿体主要赋存在控矿断裂带的构造蚀变带内，并沿断裂带倾向呈条带状展布（杨斌等，2012）。招金矿业股份有限公司采用音频大地电磁测深法（AMT），推断2 000 m深度内的断裂和地层界限，进一步证实了隐伏矿体的赋存规律（李清玉等，2019）。此外，中国冶金地质总局山东正元地质勘查院通过开展可控源音频大地电磁（CSAMT）测深，证实了招平断裂带中段向深部稳定延伸，并提出梯度带是间接找矿标志（刘永昌等，2020）；针对招平断裂带深部特征的研究进一步表明，该断裂带在2 500 m深度呈舒缓波状向下延伸，局部有倾角变化，具有明显的浅部陡、深部缓的变化趋势，且断裂产状变化部位与隐伏矿体赋存具有高度耦合性（刘振等，2023）。基于前人研究，正确认识3 000 m深度内的含矿情况是进一步开展深部找矿的现实需要（宋明春等，2015，2022a；刘向东等，2022）。然而，传统大地电磁测深方法存在一定的局限性，例如：MT和AMT方法依赖天然场源，信号微弱且抗干扰能力不足；CSAMT方法受近场效应和体积效应的影响，对埋深较大或水平展布矿体的识别能力有限（Wang et al.，2010）。此外，通过卡尼亚视电阻率单一参数反演存在非唯一性问题，难以针对复杂深部地下空间矿体进行准确预测（Pak et al.，2017）。

为了获取更多的深部矿体信息，席振铢等（2019）提出了多极化电磁测深法（MPEM）。该方法基于传统大地电磁测深理论，采用共中心点、不同极距的天然电场测量方式，结合中心点垂直磁场，计算不同极距的卡尼亚电阻率，并通过权重重新计算综合电阻率。同时，在中心点和基站同步观测正交水平磁场，结合磁场信号获得不同频率的准静态磁化率，可实现磁化率测深。MPEM方法在实际应用中展现出显著优势。在某金属矿区，该方法能够精确反演磁铁矿体及围岩的电性结构，突破了传统电磁测深方法在深部磁性异常减弱条件下的应用局限，为复杂地质条件下的矿体探测提供了可靠支持（韦洪兰等，2023）。在云南某铅锌矿的探测中，MPEM方法成功识别出3条主要断裂和极低阻异常区，证实了其抗干扰能力较传统方法明显增强（何全华等，2023）。相比传统的CSAMT、MT和AMT等方法，MPEM方法有效降低了地形影响、静态效应和电磁干扰，可同时获得电阻率和磁化率2种物性参数，有利于提高预测精度，降低勘探风险（周胜等，2023）。

后疃矿区位于大尹格庄金矿Ⅱ号矿带外围，处于深部矿体向东延伸及侧伏部位，具备良好的成矿有利条件和找矿潜力（申玉科等，2022）。然而，由于该矿区矿体埋藏深度大，区内次级构造复杂且分布规律不清，传统电磁测深方法在探测过程中面临诸多局限（高学坎等，2024）。为此，本研究采用多极化电磁测深法（MPEM）对该矿区3 000 m深度范围内金矿控矿断裂的电性和磁性异常特征进行探测，旨在揭示金矿控矿断裂的空间分布规律，为深部隐伏金矿体的空间结构解析和靶区预测提供科学依据。

1 矿区地质与地球物理特征

大尹格庄金矿位于胶东半岛西北部，华北克拉通东缘胶北隆起中部，招平断裂带中段，是胶东地区重要的超大型蚀变岩型金矿（图1）。胶北隆起主要由前寒武纪变质基底岩和中生代花岗岩组成，其中前寒武纪变质基底主要岩系有新太古代TTG质片麻岩套、胶东岩群以及少量基性—超基性侵入岩（宋明春等，2018；赵磊等，2023）。胶东岩群呈包体状产于TTG片麻岩套中，由黑云变粒岩、斜长角闪岩和角闪黑云变粒岩夹磁铁石英岩组成（Xie et al.，2014；Yang et al.，2016；Liu et al.，2022）。区内中生代构造—岩浆活动强烈，广泛分布中生代花岗岩，主要有晚侏罗世玲珑花岗岩和栾家河花岗岩，早白垩世早期郭家岭花岗岩和早白垩世晚期伟德山花岗岩（毛先成等，2023；刘向东等，2024）。

大尹格庄矿区内出露的招平断裂带上盘主要为新太古代TTG质片麻岩和胶东岩群变质岩（Mao et al.，2019），下盘主要为晚侏罗世玲珑花岗岩，含大量中生代中基性—酸性脉岩，主要成分为闪长玢岩脉和煌斑岩脉。矿区内断裂发育，主干断裂为NNE向招平断裂，次级断裂包括NW向南周家断裂、NE向大尹格庄断裂和SE向南沟断裂（申玉科等，2022；刘向东等，2024）。大尹格庄金矿床围岩蚀变广泛发育，蚀变带规模大。招平断裂下盘蚀变类型主要为黄铁矿化、绢云母化、硅化、钾长石化和碳酸盐化，这些蚀变与金矿化关系密切，是寻找金矿的重要标志，此外还发育有绿泥石化和高岭土化；上盘常表现为硅化、绿泥石化、钾长石化和碳酸盐化。招平断裂为主要容矿构造，金矿体大多产于主裂面下盘0~90 m范围内的构造蚀变带中，主要赋矿围岩为黄铁绢英岩化碎裂岩和黄铁绢英岩化花岗质碎裂岩（刘向东等，2024）。后疃矿区位于大尹格庄金矿床矿体群外围，是本次开展多极化电磁法测深的工作区。

在大尹格庄矿区采集各类岩（矿）石标本，对其进行电阻率和磁化率测定，得到各类岩（矿）石标本物性参数（表1）。胶东群变质岩位于断裂带上盘，主要由黑云变粒岩和斜长角闪岩组成，二者在铁磁性矿物含量、结构和构造上的差异导致磁化率变化范围较大。斜长角闪岩呈强磁性，黑云变粒岩则表现为弱磁性。

岩石破碎和蚀变过程中，电性和磁性显著变化，破碎与蚀变作用导致电阻率显著降低，强硅化改造作用使部分岩石的电阻率升高。含矿断裂带主要由黄铁绢英岩化碎裂岩、富多金属硫化物矿石、黄铁绢英岩和断层泥组成，整体呈低阻、低磁特征。含矿断裂带的低阻、低磁特征与上盘胶东群变质岩和下盘花岗岩形成明显对比，异常特征显著。这些电磁特性为岩性界面划分、断裂构造识别及矿体定位提供了重要依据。

2 数据采集与处理

数据采集工作使用MPMT-18多极化大地电磁仪完成，主要设备包括MPMT-18主机和ICMS-H磁场传感器，可实现1~1×10⁵ Hz频带范围内的电磁信号采集，适用于地下0~3 000 m深度电阻率和磁化率差异介质的探测。多极化电磁测深装置按图2方式布置，测线走向为110°，测点位置3组电极的布置以及磁棒的H_x 方向与测线方向保持一致，磁棒的H_y 方向与测线方向垂直；基站处磁棒的H_xb 方向与测线方向平行，磁棒的H_yb 方向与测线方向垂直。在数据采集前，上午时段进行了一次平行试验，用于验证野外数据采集时仪器的多个电道和磁道接收的信号是否具有一致性。平行试验的测试结果如图3所示，即在1×10⁵ Hz高频附近，多个电道和磁道接收的电阻率和相位略有差异，其他频段的测试结果趋于一致，且相关度均大于0.5，证明仪器多通道一致性良好。

在数据采集过程中，选用RTK测量仪对测线和测点进行定位，点位误差小于±50 cm。为了抑制静态效应，将电极按照100 m、50 m和25 m 3个电极距沿测线方向平行布置，探测剖面测线长度为4 km。磁棒方向垂直测线方向，水平误差小于±2°，磁棒与主机的距离超过5 m。基站设在距剖面1.1 km外的相对低磁区，如图1所示。观测频率范围为1~1×10⁵ Hz，通过分频多道采样技术实现不同波段电磁信号的同步采集，采样频率f_s为277.780 kHz、4 000 Hz和250 Hz，采样时间为10，300，1 200 s。

完成数据采集后，利用MPMT-18系统中的专用软件处理多极化电磁测深时间序列数据，并采用Bostick反演方法，通过约束系数和反演系数控制反演结果，获得多极化电磁测深拟剖面图（图4）。

3 资料解释与钻孔验证

基于多极化电磁测深拟剖面图和矿区岩（矿）石的电性和磁性特征（表1），得到二维反演剖面图和推断地质剖面图（图5）。图5（a）显示，右侧-2 800 m标高以浅区域视电阻率呈中—低阻特征，磁化率表现为高低交替变化的特征；左侧-900~-3 000 m标高区域视电阻率呈高阻特征，磁化率整体较低。结合物性分析，推测左侧高阻、低磁区域对应晚侏罗世玲珑花岗岩，右侧中—低阻、高低磁交替变化区域为新太古代TTG和胶东岩群变质岩系。

图5（a）中剖面左侧标高-500 m处为高阻异常，推测为变质岩经硅化作用引起的电阻率显著提升。剖面右侧-2 200 m标高处呈现为串珠状低阻异常，推断下盘晚侏罗世玲珑花岗岩体多期活动，造成上盘新太古代TTG岩系受到多期断层活动，导致接触带尤其是上盘老的基底岩体电阻率显著降低（陈大磊等，2022；贺春艳等，2022）；右侧磁化率等值线呈垂向分布，为基性—超基性岩脉侵入且岩脉走向与磁场方向接近垂直的反映。

沿招平断裂分布的断层泥、含多金属硫化物富矿石以及黄铁绢英岩化碎裂岩整体呈现为低阻低磁化特征，与上盘新太古代TTG和胶东岩群变质岩的低阻、高低磁性交替变化特征以及下盘晚侏罗世玲珑花岗岩的高阻、低磁特征差异明显。在-800~-2 900 m标高范围内，电阻率等值线存在一串珠状展布、向南东方向倾斜的中—低阻凹陷和高—低阻渐变组合异常带，伴磁化率高低值交替变化。据此推断该区域为招平断裂带，其向深部稳定延伸，在1800~2600测点附近倾角变陡，在2600测点以远倾角趋缓。

根据以往对胶东地区断裂控矿规律的研究，断裂倾角变缓的扩容空间是成矿有利部位，由浅部至深部矿体构成阶梯式成矿模式（宋明春等，2012，2022b）。本次多极化电磁测深结果显示，招平断裂带在-2 800~-3 000 m标高产状变缓，预测该区域找矿前景良好。

基于多极化电磁测深成果，设计施工了80ZK2101钻孔，孔深为3 120.35 m，开孔标高为181.95 m。钻孔于2 943.18~3 108.41 m区间发现黄铁绢英岩化蚀变带，其中3 100.06~3 101.06 m深度探获厚度为1.00 m的蚀变岩型金矿体，金品位为3.44×10^-6。该结果与根据多极化电磁测深推断的招平断裂带位置［图5（b）］（标高为-2 800~-3 000 m，对应深度为-2 9181.95~3 181.95 m）以及成矿有利部位基本一致。

4 结论

（1）在大尹格庄矿区完成1条多极化电磁测深剖面，采集1~1×10⁵ kHz频率范围的天然电磁场信号，有效反映了垂向3 000 m深度的岩（矿）石电磁特征和地质结构。

（2）采用多极化大地电磁测深法同时提取电阻率和磁化率信息，根据低阻高磁化异常带圈定了招平断裂带位置，将-2 800~-3 000 m标高异常带转折变化指示的断裂倾角变缓位置圈定为有利赋矿部位。

（3）通过80ZK2101钻孔验证，招平断裂带及其控制的矿化蚀变带和金矿体的延伸深度可达3 000 m，研究结果为下一步深部找矿突破提供了依据。

参考文献

[1] Chen Dalei， Wang Runsheng， He Chunyan，et al，2022.Application of integrated geophysical exploration in deep spatial structures：A case study of Jiaodong gold ore concentration area［J］.Geophysical and Geochemical Exploration，46（1）：70-77.

[2] Deng Xiangwei， Dai Xueling， Peng Shenglin，et al，2011.Orientation breakthrough of -1000 meters deep orebodies in Dayingezhuang Houcang segment of Zhaoping fault belt［J］.Gold，32（12）：11-17.

[3] Gao Shuai， Zeng Qingdong， Yu Changming，et al，2017.Application of remote sensing and comprehensive geophysical methods in spatial localization of hidden mentallogenic intrusions in south Zhaoyuan，Shandong Province［J］.Gold Science and Technology，25（5）：1-10.

[4] Gao Xuekan， Xi Zhenzhu， Xiang Yinhe，et al，2024.Deep ore body positioning and prediction in Houtuan mining area of Dayingezhuang gold deposit in Jiaodong［J］.Gold，45（6）：60-64.

[5] He Chunyan， Yao Zheng， Guo Guoqiang，et al，2022.Deep structural features of the Jiaobei terrane of Jiaodong and the North Sulu Orogenic Belt：The inspiration from deep exploration of the geophysics［J］.Progress in Geophysics，37（4）：1392-1404.

[6] He Quanhua， Yan Changning， Huang Xuelin，et al，2023.Application Research of audio magnetotelluric sounding （AMT） and multipolarized magnetotelluric sounding （MPMT） in exploration of a lead zinc mine in Yunnan［J］.World Nonferrous Metals，38（15）：109-111.

[7] Li Qingyu， Yang Bin， Duan Lei，et al，2019.Prospecting and prediction of magnetotelluric method for Dayingezhuang gold mine in Jiaodong［C］//Abstracts of the 9th National Symposium on Metallogenic Theory and Prospecting Methods.Guiyang：Mineral Deposit Geochemistry Committee of Chinese Society of Mineralogy.

[8] Liu Xiangdong，2022.Deep Prospecting Prediction of Typical Gold Deposits in the Zhaoping Gold Belt，Jiaodong Peninsula，China［D］.Beijing：China University of Geosciences（Beijing）.

[9] Liu X D， Ding Z J， Song M C，et al，2022.Geology and minerali-zation of the Dayingezhuang supergiant gold deposit（180 t） in the Jiaodong Peninsula，China：A review［J］.China Geology，5（4）：696-721.

[10] Liu Xiangdong， Yang Zhenliang， Xu Shaohui，et al，2024.Geology and mineralization of the Dayingezhuang supergiant gold deposit in the Jiaodong Peninsula，China［J］.Geological Review，70（Supp.1）：49-51.

[11] Liu Xiangdong， Zhou Mingling， Xu Shaohui，et al，2022.Prospecting prediction and verification at a depth of 3000 m in the Shuiwangzhuang gold deposit，northwestern Jiaodong Peninsula，eastern China［J］.Geological Bulletin of China，41（6）：946-957.

[12] Liu Yongchang， Sun Jing，2020.Application of comprehensive geophysical and geochemical exploration in deep gold exploration in the middle of Zhaoping fault zone，Shandong Province［J］.Mineral Exploration，11（3）：530-539.

[13] Liu Zhen， Du Liming， Mei Zhenhua，et al，2023.Geophysical evidence for the lateral relief pattern of the Zhaoping Fault and its significance［J］.Journal of Guilin University of Te-chnology，43（3）：397-405.

[14] Mao X C， Ren J， Liu Z K，et al，2019.Three-dimensional prospectivity modeling of the Jiaojia-type gold deposit，Jiaodong Peninsula，Eastern China：A case study of the Dayingezhuang deposit［J］.Journal of Geochemical Exploration，203：27-44.

[15] Mao Xiancheng， Wang Chuntan， Liu Zhankun，et al，2023.Three-dimensional metallogenic prediction with integration of structural reconstruction at the Dayingezhuang gold deposit，northwestern Jiaodong Peninsula［J］.Northwestern Geology，56（5）：72-84.

[16] Pak Y C， Li T L， Kim G S，2017.2D data-space cross-gradient joint inversion of MT，gravity and magnetic data［J］.Journal of Applied Geophysics，143：212-222.

[17] Shen Yuke， Guo Tao， Zhicheng Lü，et al，2022.Analysis on the ore-controlling structures and prospecting in the Dayingezhuang gold deposit，Zhaoyuan，Shandong Province［J］.Geology in China，49（1）：215-225.

[18] Song Mingchun，2015.The main achievements and key theory and methods of deep-seated prospecting in the Jiaodong gold concentration area，Shandong Province［J］.Geological Bulletin of China，34（9）：1758-1771.

[19] Song Mingchun， Ding Zhengjiang， Liu Xiangdong，et al，2022b.Structural controls on the Jiaodong type gold deposits and metallogenic model［J］.Acta Geologica Sinica，96（5）：1774-1802.

[20] Song Mingchun， Song Yingxin， Ding Zhengjiang，et al，2018.Jiaodong gold deposits：Essential characteristics and major controversy［J］.Gold Science and Technology，26（4）：406-422.

[21] Song Mingchun， Yang Liqiang， Fan Hongrui，et al，2022a.Current progress of metallogenic research and deep prospecting of gold deposits in the Jiaodong Peniusula during 10 years for Exploration Breakthrough Strategic Action［J］.Geological Bulletin of China，41（6）：903-935.

[22] Song Mingchun， Yi Pihou， Xu Junxiang，et al，2012.A step metallogenetic model for gold deposits in the northwestern Shandong Peninsula，China［J］.Science China（Earth Science），42（7）：992-1000.

[23] Wang R， Masashi E， Di Q Y，et al，2010.The analysis of CSAMT responses of dyke embedded below conductive overburden［J］.Chinese Journal of Geophysics，53（3）：677-684.

[24] Wei Honglan， Long Xia， Xi Zhenzhu，et al，2023.Multi-polarized magnetotellurics system and its application to the exploration in metal mines［J］.Gold，44（1）：68-74.

[25] Xi Zhenzhu， Long Xia， Wang Liang，et al，2019.A multi-polarisation geomagnetic bathymetry method：201910859457［P］.2019-11-22.

[26] Xie S， Xie H， Wang S，et al，2014.Ca.2.9 Ga granitoid magmatism in eastern Shandong，North China Craton：Zircon dating，Hf-in-zircon isotopic analysis and whole-rock geochemistry［J］.Precambrian Research，255：538-562.

[27] Yang Bin， Gao Xing， Peng Xinglin，et al，2012.Oriental prognosis of deep orebodies in Dayingezhuang-Houcang mining area of Zhaoping fault zone［J］.The Chinese Journal of Nonferrous Metals，22（3）：872-879.

[28] Yang L Q， Deng J， Wang Z L，et al，2016.Relationships between gold and pyrite at the Xincheng gold deposit，Jiaodong Peninsula，China：Implications for gold source and deposition in a brittle epizonal environment［J］.Economic Geology，111（1）：105-126.

[29] Zhao Lei， Zhang Rucheng， Sun Weiqing，et al，2023.A review on the research progress of the Early Precambrian supracrustal rocks in the Jiaobei Massif［J］.Acta Petrologica Sinica，39（8）：2211-2237.

[30] Zhou Sheng， Xi Zhenzhu， Long Xia，et al，2023.Principle and application of multipolarization magnetotelluric method［J］.Xinjiang Geology，41（Supp.1）：95.

[31] 陈大磊，王润生，贺春艳，等，2022.综合地球物理探测在深部空间结构中的应用——以胶东金矿集区为例［J］.物探与化探，46（1）：70-77.

[32] 邓湘伟，戴雪灵，彭省临，等，2011.招平断裂带大尹格庄—后仓段-1000 m深部矿体定位突破［J］.黄金，32（12）：11-17.

[33] 高帅，曾庆栋，于昌明，等，2017.遥感及综合物探方法用于山东招远南部隐伏成矿侵入体的空间定位［J］.黄金科学技术，25（5）：1-10.

[34] 高学坎，席振铢，向胤合，等，2024.胶东大尹格庄金矿床后疃矿区深部矿体定位预测［J］.黄金，45（6）：60-64.

[35] 贺春艳，姚铮，郭国强，等，2022.胶东胶北地体与北苏鲁造山带深部结构特征——来自地球物理深部探测的启示［J］.地球物理学进展，37（4）：1392-1404.

[36] 何全华，颜长宁，黄雪林，等，2023.音频大地电磁测深（AMT）与多极化大地电磁测深（MPMT）在云南某铅锌矿勘探中的应用研究［J］.世界有色金属，38（15）：109-111.

[37] 李清玉，杨斌，段磊，等，2019.胶东大尹格庄金矿大地电磁法找矿预测［C］//第九届全国成矿理论与找矿方法学术讨论会论文摘要集.贵阳：中国矿物岩石地球化学学会矿床地球化学专业委员会.

[38] 刘向东，2022.胶东招平金矿带典型金矿床深部成矿预测［D］.北京：中国地质大学（北京）.

[39] 刘向东，杨真亮，徐韶辉，等，2024.胶东大尹格庄超大型金矿床地质特征与成矿作用［J］.地质论评，70（增1）：49-51.

[40] 刘向东，周明岭，徐韶辉，等，2022.胶西北水旺庄金矿床3000 m深度找矿预测［J］.地质通报，41（6）：946-957.

[41] 刘永昌，孙靖，2020.综合物化探在招平断裂带中段金矿深部找矿的应用［J］.矿产勘查，11（3）：530-539.

[42] 刘振，杜利明，梅贞华，等，2023.招平断裂侧伏规律的地球物理证据及其意义［J］.桂林理工大学学报，43（3）：397-405.

[43] 毛先成，王春锬，刘占坤，等，2023.顾及构造改造的胶西北大尹格庄金矿床三维成矿预测［J］.西北地质，56（5）：72-84.

[44] 申玉科，郭涛，吕志成，等，2022.山东招远大尹格庄金矿控矿构造解析与找矿预测［J］.中国地质，49（1）：215-225.

[45] 宋明春，2015.胶东金矿深部找矿主要成果和关键理论技术进展［J］.地质通报，34（9）：1758-1771.

[46] 宋明春，丁正江，刘向东，等，2022b.胶东型金矿床控矿构造及成矿模式［J］.地质学报，96（5）：1774-1802.

[47] 宋明春，宋英昕，丁正江，等，2018.胶东金矿床：基本特征和主要争议［J］.黄金科学技术，26（4）：406-422.

[48] 宋明春，杨立强，范宏瑞，等，2022a.找矿突破战略行动十年胶东金矿成矿理论与深部勘查进展［J］.地质通报，41（6）：903-935.

[49] 宋明春，伊丕厚，徐军祥，等，2012.胶西北金矿阶梯式成矿模式［J］.中国科学（地球科学），42（7）：992-1000.

[50] 韦洪兰，龙霞，席振铢，等，2023.多极化大地电磁系统及其在金属矿勘探中的应用［J］.黄金，44（1）：68-74.

[51] 席振铢，龙霞，王亮，等，2019.一种多极化大地电磁测深方法：201910859457［P］.2019-11-22.

[52] 杨斌，高星，彭省临，等，2012.招平断裂带大尹格庄—后仓段深部矿体定位预测［J］.中国有色金属学报，22（3）：872-879.

[53] 赵磊，张儒诚，孙伟清，等，2023.华北克拉通东部胶北地体早前寒武纪表壳岩系：研究进展与综述［J］.岩石学报，39（8）：2211-2237.

[54] 周胜，席振铢，龙霞，等，2023.多极化大地电磁法原理及应用［J］.新疆地质，41（增1）：95.

基金资助

国家重点研发计划项目“碳酸盐岩容矿的锰铝矿探测预测技术与找矿模型”(2022YFC2903404)

国家重点研发计划项目“低缓磁异常弱信息提取与富铁矿体定位技术示范”(2022YFC2903704)