蓬莱黑岚沟金矿区岩浆作用对成矿制约：岩石地球化学、年代学证据

刘强　李大鹏　张岩　尉鹏飞　蔡娜　张超　谢伟

摘要：黑岚沟金矿区位于胶北地区东部，目前探明金储量占蓬莱地区的70 %以上。矿区内与成矿相关的侵入岩富硅、富碱、贫镁，TFeO（1.22 %～12.36 %，平均值3.38 %）、Rb质量分数小于270×10-6，稀土元素球粒陨石标准化配分曲线具明显的右倾分配模式，具有弱的正Eu异常，w（Rb）/w（Nb）值为7.8～36.0，为A型花岗岩，具有低镁埃达克岩的特征。LA-ICP-MS锆石U-Pb法测得侵入岩结晶年龄分别为（172.8±0.3）Ma、（172.8±0.4）Ma和（174.0±0.4）Ma，均属中侏罗世早期岩浆活动产物。黑岚沟金矿区主要受下地壳和岩石圈地幔控制，成矿物质来源与地壳的部分熔融、岩石圈地幔和软流圈地幔上涌相关，属于壳幔混合来源。

关键词：地球化学；锆石U-Pb法；Pb同位素；黑岚沟；蓬莱；成矿物质来源

中图分类号：TD11 P618.51文献标志码：A开放科学（资源服务）标识码（OSID）：

文章编号：1001-1277（2023）09-0093-10doi：10.11792/hj20230915

黑岚沟金矿区位于胶北地区东部、栖霞—蓬莱金矿带北段，由多个石英脉型金矿床组成，主要包括黑岚沟矿段、初格庄矿段、齐家沟矿段、强家沟矿段和侯格庄矿段等，矿区目前探明金储量占蓬莱地区的70 %以上。前人对黑岚沟金矿区的地球物理特征、岩浆活动的构造背景、成矿物质来源、同位素年龄、成矿规律等［1-7］方面进行了初步研究。目前，该矿区内花岗岩类地球化学演化特征及其指示大地构造背景转化的系统研究成果缺失，制约了对岩浆活动的构造驱动机制和大规模金成矿机理的深入理解。本研究在细致野外地质调查工作基础上，采集黑岚沟金矿区内与成矿关系密切的侵入岩样品及典型矿石，开展岩相学、岩石地球化学和锆石U-Pb年代学等方面的研究，厘定侵入岩的成岩（矿）时代，揭示岩体起源及其形成环境，探讨岩石成因与演化及其对成矿的制约，该研究对于进一步认识胶西北地区燕山晚期岩浆演化与金矿成因具有一定的科学意义。

1 区域地质背景及矿床地质

黑岚沟金矿床位于华北地台郯庐深大断裂带东侧，胶北隆起区的栖霞复背斜北翼东段，其成矿地质条件受控于所属成矿带的地层、构造、岩浆岩。矿区内地层出露较简单，构造运动较频繁，岩浆活动强烈（见图1）。

矿区内矿体多位于岩体的构造破碎蚀变带内，石英脉多以充填为主，蚀变岩以热液交代为主。矿体总体呈脉状、网脉状、透镜状，产状受构造破碎蚀变带控制。矿体在构造破碎蚀变带中的赋存状态具有一定规律性，在构造膨胀部位围岩蚀变越强烈，则金矿化越明显。在构造膨胀部位存在上、下2个结构面（断层面），夹持在结构面之间的蚀变体，是由几种碎裂和蚀变程度不同的蚀变岩石所组成。靠近上、下结构面，岩石破碎、裂隙密集，多穿插脉状或网脉状石英脉，形成夹石英脉的黄铁绢英岩，其构成了矿体的主体。

2 样品及分析方法

2.1 样品采集及特征

样品采集主要在黑岚沟矿段、齐家沟矿段和强家沟—初格庄矿段的井下及岩芯中选择侵入岩样品及矿体围岩。共采集了侵入岩样品10件，矿石样品3件。

黑岚沟金矿区与成矿作用相关的侵入岩主要为燕山早期的侏罗世花岗岩（见图2），在地表分布于矿区东部燕子夼—磁山一带，井下大面积分布为中粗粒二长花岗岩。岩石呈灰白色、灰色，中粗粒结构，块状构造。主要矿物成分及占比为：钾长石30 %～35 %，斜长石30 %～35 %，石英20 %～25 %，黑云母5 %～10 %。

矿区矿石成分简单，金（见图3）赋存于石英脉中，金属矿物主要有自然金、银金礦、黄铁矿、方铅矿、闪锌矿，以及少量黄铜矿、磁黄铁矿；脉石矿物主要有石英、长石、绢云母、绿泥石、方解石等。自然金常与金属硫化物共生，由于部分颗粒较大，形成的明金肉眼可见，但分布不均，往往在矿石金品位很富的地方可见。

2.2 分析方法

全岩主量元素和微量元素分析测试工作在山东省地质科学研究院完成。主量元素测试采用X射线荧光光谱法（XRF），使用仪器为Axiosm AX型X射线荧光光谱仪。微量元素和稀土元素测试分析采用电感耦合等离子体质谱法（ICP-MS），使用的仪器为Element XR等离子体质谱仪。

锆石分选工作在廊坊市峰泽源岩矿检测技术有限公司完成。样品靶的制作、透射光、反射光和阴极发光由北京锆年领航科技有限公司完成，锆石U-Pb年龄测试分析在山东省地质科学研究院的激光剥蚀-多接收电感耦合等离子体质谱仪（LA-ICP-MS）上完成。锆石U-Pb测年数据处理采用ICP-MS DataCal程序，U-Pb年龄谐和图采用Isoplot软件绘制。

3 结果及分析

3.1 岩石地球化学分析

对黑岚沟金矿区侵入岩样品进行了岩石地球化学全分析，结果见表1。

3.1.1 主量元素特征

黑岚沟金矿区侵入岩石按SiO2质量分数可以分为2类，一类样品SiO2质量分数高，为62.08 %～73.48 %，平均值69.51 %；另一类SiO2质量分数较低，仅有51.80 %。TFeO质量分数较低，为1.22 %～12.36 %，平均值3.38 %；Al2O3质量分数13.45 %～18.40 %，平均值15.28 %；Na2O质量分数0.15 %～5.67 %，平均值3.39 %；烧失量较低。将样品进行碱－硅（TAS）图解投图，结果见图4。由图4可知：样品大多投在花岗岩区域内。从岩石系列上看，样品落在了高钾钙碱性—钙碱性系列区域内（见图5），属于偏铝质—过铝质岩石（见图6）。

3.1.2 稀土元素和微量元素特征

稀土元素（REE）质量分数为85.63×10-6～241.78×10-6，平均值157.81×10-6，在稀土元素球粒陨石标准化配分模式图（见图7-a））中具有相似的稀土配分特征，配分曲线显示为明显右倾到相对平坦，w（LREE）/w（HREE）=9.83～22.96，w（La）N/w（Yb）N =14.38～47.19，富集轻稀土元素，相对亏损重稀土元素，轻、重稀土元素分异现象较明显。矿石稀土元素大部分表现出弱的负Eu异常，δEu=0.81～1.15，少部分具有弱的正Eu异常，δEu=1.04～1.15。花岗质围岩样品表现出较弱的负Ce异常，δCe=0.87～0.92，表明该岩体遭受了轻微的蚀变及风化作用。

在微量元素原始地幔标准化蛛网图（见图7-b））中，所有矿石样品均表现出明显的富集大离子亲石元素Rb、K和高场强元素Nd、Ta、Zr、Hf，具有高的Sr/Y值（w（Sr/Y）=11～115.59，平均值50.13）。

3.2 锆石U-Pb年代学

本次研究选择了3件样品进行锆石U-Pb同位素测试分析，锆石阴极发光（CL）图像（见图8）显示，锆石以长柱状、板柱状晶体为主，个别为短柱状，粒径多在100～150 μm，长宽比多为3∶1～3∶2，少量为3∶1～5∶1。锆石晶体轮廓清晰，晶面多数光滑，内部结构清楚，振荡环带、核幔结构较发育，部分锆石颗粒发育裂隙，表现出典型的岩浆锆石特征。

样品测试数据见表2。齐家沟矿段花岗岩样品QJG-6（见图9-a））14个测点的206Pb/238U-207Pb/235U年龄和谐，集中位于谐和线附近，其谐和年龄值为（172.8±0.3）Ma，MSWD=0.9；206Pb/238U加权平均年龄为（171.8±2.5）Ma，MSWD=0.64。齐家沟矿段花岗岩样品QJG-12（见图9-b））14个测点的206Pb/238U-207Pb/235U年龄和谐，集中于谐和线附近，其谐和年龄值为（172.8±0.4）Ma，MSWD=1.7；206Pb/238U加权平均年龄为（172.87±0.47）Ma，MSWD=0.96。强家沟矿段花岗岩样品D114-1（见图9-c））11个测点的206Pb/238U-207Pb/235U年龄和谐，集中于谐和线附近，其谐和年龄值为（174.0±0.4）Ma，MSWD=0.6；206Pb/238U加权平均年龄为（173.94±0.35）Ma，MSWD=0.6。黑岚沟金矿区侵入体的岩浆结晶年龄均属燕山早期中侏罗早期岩浆活动的产物。

4 讨 论

4.1 成岩时代

黑岚沟金矿区存在3期侵入岩体，形成时代分别为古元古代、早侏罗世、早白垩世，对应岩体分别为大柳行序列燕子夼单元二长花岗岩、玲珑序列郭家店单元二长花岗岩和郭家岭序列罗家单元二长花岗岩。本次研究采集的侵入岩样品为井下与矿化相关的侏罗纪岩体，未发生变质，且锆石阴极发光图像岩浆振荡环带清晰，所测样品的数据点均位于U-Pb年龄谐和曲线上，研究结果将齐家沟矿段花岗岩结晶年龄限定为（172.8±0.3）Ma，将强家沟矿段二长花岗岩结晶年龄限定为（174.0±0.4）Ma，表明这2处岩体是中侏罗世岩浆活动的产物，属燕山早期侵入岩。

从鲁西与胶东岩浆侵位时间来看，XU等［14］认为铜石岩体和铜井岩体的锆石U-Pb年龄分别为（177±4）Ma和126～132 Ma，ZHANG等［15］获得了铜石岩体的Rb-Sr等时线年龄为（185.4±2.0）Ma，胡华斌等［16］报道了铜石闪长玢岩的锆石U-Pb年龄为（175.7±3.8）Ma，LAN等［17］报道了铜石岩体的石英二长岩和正长斑岩的年龄为180.1～184.7 Ma，郭谱［18］获得铜井岩体U-Pb年龄为175～178 Ma。这些年龄代表了铜石杂岩体的侵位时间，而本次在黑岚沟金矿区发现的172～174 Ma的岩浆侵位活动与鲁西地区铜石杂岩体的侵位时间基本一致，揭示着胶东地区也存在着与鲁西同时代岩浆活动，从胶东和鲁西地区的成矿物质来源是否相同这一角度思考，值得进一步研究。已有研究認为，玲珑序列岩体的年龄在150～160 Ma，本次获得的年龄（172～174 Ma）要早于玲珑序列岩体的活动范围，这一发现是对胶东地区一直以来没有发现160～170 Ma岩浆活动的补充，具有很好的科学研究意义。

4.2 构造背景

在晚三叠世—晚侏罗世，板块俯冲后碰撞的巨大挤压作用导致华北克拉通南缘地壳增厚，早侏罗世末期，郯庐断裂开始产生并发生左行平移运动，这时候玲珑花岗岩强力侵位，而郯庐断裂带发生了巨大的左行平移，将胶北地块从蚌埠—五河的东部推移到现在的位置。黑岚沟金矿区侵入岩亦随着郯庐断裂带的活动而侵位形成，在其w（K2O）-w（SiO2）图解（见图10-a）中样品基本落入大陆花岗斑岩区域内，在侵入岩w（Rb）-w（Y+Nb）判别图解（见图10-b）上显示其形成于火山弧—同碰撞构造环境；在A/NK-A/CNK图上显示为大陆弧与后造山花岗岩，以上特征显示黑岚沟金矿区侵入岩的岩浆源区主要为重熔型花岗质岩浆，成岩环境处于挤压构造背景。到了早白垩世，郯庐断裂带由左行压扭向右行拉张转化，导致胶北地块郭家岭岩体侵位。同时，郯庐断裂带发生了大规模的拉张，导致沂沭裂谷、胶莱盆地、青山群火山岩、金矿床等形成和隆起区大量脉岩侵入。

从侵入岩的岩石地球化学分析结果来看，主量元素显示了高硅、高碱、富钾、贫钙、贫镁及低钛的特征，w（Fe）/w（Fe+Mg）值较高，w（K2O）/w（Na2O）值和w（K2O）较高，微量元素Rb、Nd、Zr质量分数较高，Ba、Sr、Ti质量分数较低，Eu异常不明显，符合A型花岗岩特征。本次研究的花岗岩具有较高的w（TFeO）（1.22 %～12.36 %，平均值3.38 %），且Rb质量分数（92.10×10-6～269.00×10-6）小于270×10-6，这与高分异的I型花岗岩（低TFeO，<1.00 %）有区别。花岗岩具有较低的P2O5质量分数（平均值0.11 %）和较高的Na2O质量分数（平均值为3.39 %），这与S型花岗岩区别明显。在10 000w（Ga）/w（Al）-w（TFeO）/w（MgO）图解（见图11）中，样品大部分也均落入A型花岗岩区域内，因此可以判定侵入岩属于A型花岗岩。

4.3 岩浆作用对成矿制约

黑岚沟金矿区侵入围岩的主量元素、微量元素分析表明，其均具有高SiO2、高Al2O3、低MgO，富集大离子亲石元素（K、Rb、Ba）和轻稀土元素，亏损高场强元素（P、Ti）和重稀土元素，以及高Sr、低Y和Yb的特征，A/NK-A/CNK图解中样品落入偏铝质—过铝质区域内，在w（K2O）-w（SiO2）图解上样品落入高钾钙碱性—钙碱性区域内。在稀土元素球粒陨石标准化配分模式图上具有较为明显的右倾分配模式，轻重稀土元素分异现象较明显。这些特征显示出低镁埃达克岩的特点（见图12）。花岗岩较高的w（La）N/w（Yb）N值和亏损重稀土元素的特征指示了该岩体的源区残留相中含有石榴石和/或角闪石，同样该岩体稀土元素表现出弱的正Eu异常，这区别于典型的板内花岗岩的“V”型配分模式。花岗岩中w（Rb）/w（Nb）值为7.8～36.0，高于地壳的平均值（4.50）。以上特征表明，黑岚沟金矿区侵入岩主要为壳源物质熔融的产物。

对于胶西北地区金矿或者扩大到胶东地区金矿的成矿物质来源，翟建平等［23-24］认为，成矿流体来源于大气降水，成矿物质来源于围岩，包括花岗岩、变质岩；张连昌等［25］认为，成矿流体来源于与中基性脉岩有关的幔源流体，成矿物质主要来源于地幔；吕古贤等［26-27］认为，晚太古代—元古代的海底基性火山岩及沉积岩系是胶东金矿的初始矿源岩，胶东前寒武纪变质岩系是金矿的中间矿源岩，而主要在中生代期间形成的同源交代重熔花岗岩是直接矿源岩；侯明兰等［3］对蓬莱地区金矿的黄铁矿和方铅矿单矿物进行了系统的Pb同位素研究，认为黑岚沟金矿区的成矿流体与围岩（如新太古代胶东岩群地层）的水岩交换反应作用更强，更多上地壳物质进入了成矿流体，幔源贡献可能稍少，壳源组分混入稍多。DENG等［28］认为，流体和成矿金属并非主要来源于已经变质的地壳岩石或胶东半岛大量的花岗质岩浆作用，而可能来源于下地壳，即早前俯冲事件中被交代和富集的地幔岩石圈。

基于岩石学、地球化学和U-Pb年代学的综合研究，认为黑岚沟金矿区成矿主要受下地壳和岩石圈地幔控制。西太平洋板块向东亚大陆边缘俯冲导致胶东地区岩石圈增厚，加厚下地壳的进一步熔融和侵位形成岩体。由太平洋板块俯冲作用形成的弧后伸展构造环境，地幔岩浆上升造成地壳物质部分熔融产生了酸性岩浆房［28-29］。因此，胶北地区金矿的成矿物质来源与地壳部分熔融、岩石圈地幔和软流圈地幔上涌密切相关，属于壳幔混合来源。

5 结 论

1）黑岚沟金矿区主要赋矿侵入岩显示高钾钙碱性—钙碱性系列花岗岩特征，属偏铝质—过铝质A型花岗岩，具低镁埃达克岩特征，是壳源物质熔融的产物，其岩浆源区主要為重熔型花岗质岩浆，成岩环境处于挤压构造背景。

2）黑岚沟金矿区典型中酸性侵入岩的结晶年龄分别为（172.8±0.3）Ma、（172.8±0.4）Ma和（174.0±0.4）Ma，均属中侏罗世早期岩浆活动的产物，比前人认为的玲珑序列花岗岩的形成年龄早约10 Ma。

3）黑岚沟金矿区成矿主要受下地壳和岩石圈地幔控制，成矿物质来源与地壳部分熔融、岩石圈地幔和软流圈地幔上涌密切相关，属于壳幔混合来源。

致谢：在野外调查期间得到了山东省核工业二七三地质大队等地质同行的大力协助。中国地质大学（北京）李小伟教授对本文提供了很好的建议；本文撰写的过程中引用了大量文献，由于篇幅所限，未能一一列出，在此深表感谢！

［参 考 文 献］

［1］孙永深.黑岚沟矿区金矿成矿规律探讨［J］.黄金，1999，20（6）：13-17.

［2］毛景文，李厚民，王义天，等.地幔流体参与胶东金矿成矿作用的氢氧碳硫同位素证据［J］.地质学报，2005，79（6）：839-857.

［3］侯明兰，蒋少涌，姜耀辉，等.胶东蓬莱金成矿区的S-Pb同位素地球化学和Rb-Sr同位素年代学研究［J］.岩石学报，2006，22（10）：2 525-2 533.

［4］侯明兰，蒋少涌，沈昆，等.胶东蓬莱金矿区流体包裹体和氢氧同位素地球化学研究［J］.岩石学报，2007，23（9）：2 241-2 256.

［5］张宗庆.山东蓬莱黑岚沟金矿田断裂控矿特征浅析［J］.科技创新导报，2008（24）：97-98.

［6］马顺溪，白宜娜，孙永联，等.胶东蓬莱大柳行金矿田燕山金矿床流体包裹体特征和氢-氧同位素研究［J］.地质学报，2020，94（11）：3 391-3 403.

［7］吴晋超，申俊峰，申玉科，等.胶东黑岚沟金矿田黄铁矿热电性研究及深部成矿预测［J］.西北地质，2021，54（2）：111-125.

［8］MIDDLEMOST E A K.Naming materials in the magma/igneous rock system［J］.Earth-Science Reviews，1994，37（3/4）：215-224.

［9］MAITRE R W，BATEMAN P，DUDEK A，et al.A Classification of igneous rocks and glossary of terms：Recommendation of the international union of the geological subcommission on the systematic of igneous rocks［M］.Oxford：Blackwell Scientific Publications，1989.

［10］RICKWOOD P C.Boundary lines within petrologic diagrams which use oxides of major and minor elements［J］.Lithos，1989，22（4）：247-263.

［11］MANIAR P D，PICCOLI P M.Tectonic discrimination of granitoids［J］.Geological Society of American Bulletin，1989，101（5）：635-643.

［12］TAYLOR S R，MCLENNAN S M.The continental crust：Its composition and evolution［M］.Oxford：Blackwell Scientific Publications，1985.

［13］MCDONOUGH W F，SUN S S.The Composition of the Earth［J］.Chemical Geology，1995，120：223-253.

［14］XU Y G，HUANG X L，MA J L，et al.Crust-mantle interaction during the tectono-thermal reactivation of the North China Craton：Constraints from SHRIMP zircon U-Pb chronology and geochemistry of Mesozoic plutons from western Shandong［J］.Contribution to Mineralogy and Petrology，2004，147（6）：750-767.

［15］ZHANG H F，SUN M，ZHOU X H，et al.Geochemical constraints on the origin of Mesozoic alkaline intrusive complexes from the North China Craton and tectonic implications［J］.Lithos，2004，81（1）：297-317.

［16］胡華斌，毛景文，牛树银，等.鲁西平邑地区浅成热液金矿床成矿作用［M］.北京：地震出版社，2007.

［17］LAN T G，FAN H R，SANTOSH M，et al.Crust-mantle interaction beneath the Luxi Block，eastern North China Craton：Evidence from coexisting mantle and crust-derived enclaves in a quartz monzonite pluton［J］.Lithos，2013，177：1-16.

［18］郭谱.鲁西中生代金成矿的地球动力学背景研究［D］.北京：中国地质大学（北京），2014.

［19］COLEMAN R G，PETERMAN Z E.Oceanic plagiogranite［J］.Journal of Geophysical Research Atmospheres，1975，80（8）：1 099-1 108.

［20］PEARCE J A，HARRISS B W，TINDLE A G.Trace element discrimination diagrams for the tectonic interpretation of granitic rocks［J］.Journal of Petroleum Geology，1984，25（4）：956-983.

［21］JOSEPH B W，KENNETH L，CURRIE B W C.A-type Granites：Geochemical characteristics，discrimination and petrogenesis［J］.Contributions to Mineralogy and Petrology，1987，95（4）：407-419.

［22］MARTIN H.Adakitic Magmas：modern analogues of Archaean gra-nitoids［J］.Lithos，1999，46（3）：411-429.

［23］翟建平，徐光平，胡凯.栖霞金矿矿物、流体和同位素特征及意义［J］.矿床地质，1998，17（4）：20-26.

［24］卢焕章，GUHA J，方根保.山东玲珑金矿的成矿流体特征［J］.地球化学，1999，28（5）：421-437.

［25］张连昌，沈远超，李厚民，等.胶东地区金矿床流体包裹体的He、Ar同位素组成及成矿流体来源示踪［J］.岩石学报，2002，18（4）：559-565.

［26］吕古贤，韦昌山，郭涛，等.胶东矿集区金矿成矿地质事件研究初探［J］.黄金地质，2004（2）：1-7.

［27］李秀章，王勇军，李衣鑫，等.胶东蓬莱黑岚沟金矿床黄铁矿微区地球化学特征及对成矿流体的启示［J］.地质通报，2022，41（6）：1 023-1 038.

［28］DENG J，YANG L Q，GROVES D I，et al.An integrated mineral system model for the gold deposits of the giant Jiaodong province，eastern China［J］.Earth-Science Reviews，2020，208（2）：103274.

［29］TAN J，WEI J H，GUO L L，et al.LA-ICP-MS zircon U-Pb dating and phenocryst EPMA of dikes，Guocheng，Jiaodong Peninsula：Implications for North China Craton lithosphere evolution［J］.Science in China（Series D），2008，51：1 483-1 500.

Constraints of magmatism on mineralization in Heilangou gold deposit，Penglai：Evidence from lithogeochemistry and chronology

Liu Qiang1，2，3，Li Dapeng1，2，3 ，Zhang Yan1，2，3，Wei Pengfei1，2，3，Cai Na1，2，3，Zhang Chao  1，2，3，Xie Wei1，2，3

（1.Shandong Institute of Geological Sciences;2.MNR Key Laboratory of Gold Mineralization Processes and Resource Utilization;3.Shandong Provincial Key Laboratory of Metallogenic Geological Process and Resource Utilization）

Abstract：Heilangou gold deposit is located in the east of Jiaobei area，and the proven gold reserves currently account for more than 70 % of that of Penglai area.The intrusive rocks related to mineralization in this area are rich in silicon，and alkali while lacking magnesium，with TFeO（1.22 %-12.36 %，with an average of 3.38 %），and the Rb mass fraction is less than 270×10-6.The REE chondrite-normalized distribution curve has an obvious right-leaning distribution model and weak positive Eu anomaly，and the w（Rb）/w（Nb） value is 7.8-36.0，which is A-type gra-nite，with the characteristic of low-Mg adakites.The crystal ages of intrusive rocks measured by the LA-ICP-MS zircon U-Pb method are respectively（172.8±0.3）Ma，（172.8±0.4）Ma，and（174.0±0.4）Ma，which are all related to magmatic activity in the early Middle Jurassic.Heilangou gold area is mainly controlled by the lower crust and lithosp-heric mantle，and the source of ore-forming material is related to the partial melting of the crust and upwelling of lithospheric mantle and asthenosphere mantle，which belongs to the mixed material source of crust and mantle.

Keywords：geochemistry;zircon U-Pb method;Pb isotope;Heilangou;Penglai;source of ore-forming material

收稿日期：2023-04-30；  修回日期：2023-05-30

基金項目：国家自然科学基金项目（42172094）；山东省自然科学基金项目（ZR2020MD032，ZR2020QD028）；山东省重点研发计划项目（2022CXPT047）；山东省地勘基金项目（鲁勘字〔2021〕7号，鲁勘字〔2022〕9号）

作者简介：刘 强（1983—），男，高级工程师，博士，从事地质矿产研究工作；E-mail：lqiang715@163.com

*通信作者：张 岩（1985—），男，高级工程师，博士，从事地质矿产、矿床学研究工作；E-mail：sddkyzy@163.com