Zircon U-Pb dating of the Yuangezhuang pluton in Muping of Shandong Province and its constraints on mineralization of Cu-Mo deposits
-
摘要:
牟平院格庄岩体位于胶东东部，为燕山晚期形成的与斑岩型-矽卡岩型铜钼多金属矿床有关的岩体之一。获得院格庄岩体崖西、后野、西上寨3个单元二长花岗岩LA-ICP-MS锆石U-Pb定年结果分别为（118.10±0.66）Ma、（118.52±0.78）Ma、（118.80±0.67）Ma，皆为早白垩世晚期。岩石地球化学研究表明：岩石属于高钾钙碱性、准铝质系列（K2O/Na2O值为1.00~1.29、A/CNK=0.98~1.02）。轻重稀土元素分馏明显（（La/Yb）N=26.91~35.75），δEu=0.84~0.93，具弱负铕异常。高场强元素Nb、Ti、P元素强烈亏损，大离子亲石元素Rb、K、U富集。地球化学特征显示院格庄岩体为I型花岗岩，为壳幔同溶型花岗岩。该岩体对其周缘铜钼矿空间分布和成矿时代有显著的制约作用，具有铜钼多金属矿的成矿专属性。
Abstract:The Yuangezhuang Pluton located in the east of Jiaodong area,is formed in the late Yanshanian and is one of the granites associated with Cu-Mo polymetallic deposits hosted inporphyry and skarn. LA-ICP-MS zircon U-Pb dating of Yaxi unite,Houye unite and Xishangzhai unite from the Yuangezhuang monzonitic granite yields (118.10±0.66) Ma,(118.52±0.78) Ma,and (118.80±0.67) Ma respectively,indicating the Late Epoch of Early Cretaceous. The results of petrochemistry indicate the pluton is belong to high-K calc-alkaline serie (K2O/Na2O=1.00-1.29,A/CNK=0.98-1.02). The results of trace elements are characterized by dramatic fractionation between light REE and heavy (REE(La/Yb)N=26.91~35.75 unapparent Eu anomalies (δEu=0.84-0.93),enriched in LILE (i.e. Rb,K,U) and loss in HSFE (i.e. Nb,Ti,P). The geochemical characteristics indicate the Yuangezhuang pluton is an I type granite,crust-mantle syntectic granite. The pluton has strong constrain on spatial distribution and metallogenic epoch of Cu-Mo deposits.
-
1. 引言
中生代以来，胶东地区经历了多期次构造岩浆活动和成矿作用（杨立强等，2014；丁正江等，2015；宋明春等，2018），李洪奎等（2013）提出将成矿作用分为三期，分别为挤压环境下的金钼矿成矿作用、挤压—伸展环境的金银矿成矿作用和后造山环境的金多金属矿成矿作用。胶东特有的成矿环境造就了丰富的金及多金属矿产资源，其中累计探明金矿资源储量4500 t（宋英昕等，2017）。研究区地处山东省烟台市牟平区西南，华北克拉通东缘（谢士稳等，2015），区内左行走滑深大断裂—桃村断裂将本为一体的院格庄岩体和牙山岩体切割为两个独立单元（Cui et al., 1994）（图 1），近年来随着其内部及周缘尚家庄钼矿、杏山北钼矿、孔辛头铜钼矿等多金属矿床的发现（李杰等，2013；李超等，2016），为该区多金属矿找矿注入了新活力。本文以中国地质调查局天津地质调查中心“胶东成矿区栖霞—牟平地区地质矿产调查”所属子项目“山东1∶5万水道、观水、冯家幅矿产地质调查”为依托，在详细野外调查基础上，借助LA-ICP-MS锆石U-Pb测年及地球化学的系统研究，进一步探讨院格庄岩体的成因机制，分析其与周缘铜钼多金属矿床成因、成矿时间和空间上的联系，为胶东铜钼多金属矿成矿地质背景研究提供基础地质资料，为该区开展多金属矿产勘查提供技术支持。
![]()
图 1 胶东主要金和铜钼矿床分布图（据李超等, 2016）Figure 1. Regional geological map of Jiaodong area showing the distribution of main gold, copper and molybdenum deposits(after Li Chao et al., 2016)2. 区域地质概况
研究区地层出露较为简单，主要为古元古代荆山群一套碳酸盐岩建造组合和中生代莱阳群一套砾岩、砂砾岩组合，两者呈角度不整合接触，局部断层接触。区内断裂构造发育，经历多期叠加和改造，按其展布方向大致可分为NE、NW及近EW向3组，其发育程度差别悬殊。比较完整和突出的是NE断裂构造，以朱吴断裂与桃村断裂为代表。朱吴断裂位于研究区SE侧，全长96 km，在区内出露长约9 km，宽20~100 m，走向20°左右，倾向NW，倾角70°，具波状弯曲特征，经历多次活动，中生代燕山晚期活动最为强烈，发生了左旋错动（张勇等, 2013）；桃村断裂产出于研究区W侧，全长86 km，在区内出露长约4 km，宽20~100 m，走向30~45°，倾向SE，倾角60~80°，经历了压—张—压扭性的构造活动，以晚期左行压扭性影响最大，其左行走滑位移达30 km（黄永华等，2007；张岳桥等，2007），为院格庄岩体的西界。
区内侵入岩分布广泛，主要有早白垩世伟德山序列含斑二长花岗岩（院格庄岩体）、晚侏罗世玲珑序列弱片麻状二长花岗岩（鹊山岩体）及新太古代栖霞序列条带状英云闪长质片麻岩，其中，院格庄岩体被挟持在桃村和朱吴断裂之间，自外向内划分为3个单元：崖西斑状中粒含角闪二长花岗岩、后野巨斑状中粒含角闪二长花岗岩和西上寨含巨斑中粒黑云二长花岗岩（王世进等, 2009），三者呈同心环状产出，半椭圆状展布。其中斑状中粒含角闪二长花岗岩与荆山群野头组大理岩接触带附近常见矽卡岩化蚀变，周边产出有下雨村金矿、孔辛头铜钼矿和杏山北钼矿床（图 2）。
![]()
图 2 院格庄地区区域地质图（据山东省第三地质矿产勘查院, 1990）Figure 2. Regional geological map of the Yuangezhuang area((afte No.3 Exploration Institute of Geology and Mineral Resources, Shandong Province, 1990)3. 院格庄岩体岩相特征
院格庄岩体位于下雨村西北，西邻回里东至王家疃，南接要捷村北达南周山，面积77.52 km2，自外向内可划分为3个单元：崖西单元、后野单元和西上寨单元（图 2）。各单元岩石特征进行分述如下：
崖西单元：位于院格庄岩体边部，侵入于荆山群和玲珑序列九曲单元中，界线清晰。岩性为似斑状中粒含角闪二长花岗岩，浅肉红色—灰色，似斑状结构，不等粒半自形粒状结构，块状构造；斑晶成分为钾长石（15%~20%），大小9~20 mm；基质主要矿物为石英（20%~23%），斜长石（41%~55%）、钾长石（12%~23%）、黑云母（2%~6%）、角闪石（3%~8%），粒径2.1~5.0 mm。该单元大斑晶较少，一般情况下小斑晶较多且密集（图 3a）；在斑晶中，可见其俘获的黑云母微晶；矿物颗粒相对较细，角闪石自形程度高；其内见细粒闪长质包体（图 3b），呈椭圆状产出，一般直径5~20 cm，局部呈群出现。
![]()
图 3 院格庄岩体野外地质特征a—崖西单元钾长石斑晶密集分布；b—崖西单元闪长质包体；c—后野单元钾长石巨斑呈定向组构；d—后野单元钾长石巨斑；e—西上寨单元岩石野外地质特征；f—西上寨单元钾长石斑晶Figure 3. Geological characteristics of Yuangezhuang plutona-K-feldspar phenocrysts are densely distributed Yaxi unite; b-Diorite inclusions in Yaxi unite granite; c-K-feldspar phenocrysts of Houye unite are in a directional arrangement; d-K-feldspar phenocrysts of Houye unite; e-The geological characteristics of Xishangzai unite granite in the field; f-K-feldspar phenocrysts in Xishangzai unite后野单元：与崖西单元未见明显界线，岩性为巨斑状中粒含角闪二长花岗岩，浅肉红色，似斑状结构，不等粒半自形粒状结构，块状构造；钾长石巨斑占10%~20%，粒径粗大，一般2 cm×4 cm~3 cm×6 cm，大者可达4 cm×8 cm或更大，略具定向性；基质主要矿物成分为斜长石（37%~56%）、钾长石（19%~35%）、石英（18%~28%）、黑云母（2%~5%）、角闪石（＜5%），粒径一般2~5 mm。该单元钾长石斑晶增大，局部见其沿长轴方向定向排列（图 3c）；矿物颗粒均匀性差，黑云母微晶增多，为大斑晶所俘获；局部充填长英质细脉（图 3d）；闪长质包体零星产出。
西上寨单元：位于院格庄岩体中心，同后野单元无明显界线。地表出露部分由于形成的释重节理和风化作用，呈现上部风化下部完整的形态，且界线明显（图 3e）。岩性为含巨斑细中粒含黑云二长花岗岩，浅肉红色，似斑状结构，块状构造；钾长石斑晶大小一般为1 cm×2 cm，大者可达3 cm×6 cm，含量在10%~15%，大小悬殊，但分布均匀；基质矿物成分：斜长石（30%~49%）、钾长石（21%~30%）、石英（21%~34%）、黑云母（＜5%）、角闪石（＜5%），粒径为2~5 mm。该单元斑晶大，分布均匀，含量明显减少(图 3f)；闪长质包体极为少见。
4. 样品的采集与分析方法
本次用于锆石U-Pb测年的3件样品，分别位于崖西单元、后野单元、西上寨单元中心部位的新鲜花岗岩，分别编号2017N1、2017N2、2017N3，岩性特征明显。取样地理坐标分别为：37°19′55.01″ N, 121°23′59″E; 37°17′11.53″N, 121°25′29.46″E; 37°18′28.54″N, 121°22′30.72″E(图 2)。
岩石样品单矿物挑选、制靶、CL照相等项目在廊坊市宏信地质勘查技术服务有限公司完成，在双目镜下挑出锆石，每件样品选取代表性锆石200颗进行制靶；LA-ICP-MS锆石U-Pb年龄测定在中国地质调查局天津地质调查中心实验室进行，测试采用的激光剥蚀束斑直径为32 μm，激光剥蚀60 s, 实验中采用He作为剥蚀物质的载气。锆石年龄采用国际标准锆石91500作为外标，元素含量采用NISTSRM610作为外标。29Si作为内标元素锆石中SiO2的质量分数为32.8%，分析方法见文献（Yuan et al., 2004）；普通铅校正采用Anderson推荐的方法（Anderson, 2002）；样品的同位素比值及元素含量计算采用ICP-MS-DATECAL程序（Liu et al., 2008），年龄计算及谐和图的绘制采用Isoplot程序（Ludwing，2003）。
本次地球化学分析样品分别于崖西、后野、西上寨单元内各取2件，其中崖西单元样品编号为YX1、YX2；后野单元样品编号HY1、HY2；西上寨单元样品编号为XSZ1、XSZ2（图 2）。新鲜岩石进行主量、微量和稀土分析，在山东省地质科学院测试中心完成，主量元素使用熔片法X-射线荧光光谱法（XRF）测试，分析准确度和精度优于5%；微量元素及稀土元素利用酸溶法制备样品，使用等离子体质谱法（ICP-MS）测试，分析准确度和精度优于10%；烧失量采用重量法（GR）测试；Fe2O3、FeO采用容量法（VOL）测试，As、Sb采用氢化物-原子荧光光谱法（HG-AFS）测试，分析准确度和精度优于10%。
5. 分析结果
5.1 锆石U-Pb年龄
本次选取的3件样品中锆石粒度在150~200 μm，长短轴比一般为1.5∶1~3∶1，呈浅褐黄色，半透明或透明，自形短柱状及双锥状晶体，晶棱及晶面清楚，晶型良好。2017N1内选取了26个锆石开展LA-ICP-MS U-Pb同位素测年分析，锆石震荡环带清晰，Th/U比值0.50~1.57，为典型的岩浆锆石。206Pb/238U表面年龄加权平均值（118.10±0.66）Ma（MSWD=1.6）（图 4，图 5，表 1）；2017N2内选取了20个锆石开展LA-ICP-MS U-Pb同位素测年分析，锆石震荡环带清晰，Th/U比值0.56~1.35，均大于0.1，锆石震荡环带明显，为岩浆锆石，206Pb/238U表面年龄加权平均值（118.52±0.78）Ma（MSWD=1.8）（图 6，图 7，表 1）；2017N3内选取了29个锆石微区开展LA-ICP-MS U-Pb同位素测年分析，锆石震荡环带清晰，Th/U比值0.26~1.35，均大于0.1，属于岩浆锆石，206Pb/238U表面年龄加权平均值（118.80±0.67）Ma（MSWD=1.9）（图 8，图 9，表 1）。3件样品的锆石206Pb/238U年龄与207Pb/235U年龄在误差允许范围内基本一致，在206Pb/238U-207Pb/235U谐和图（图 5，图 7，图 9）上分析点均落在谐和线及其附近。实验结果显示院格庄岩体锆石的加权平均年龄在118.10~118.80 Ma，谐和年龄在115~122 Ma。
![]()
图 4 崖西单元样品2017N1锆石阴极发光图像Figure 4. Representative CL images and U-Pb ages of zircons from samples 2017N1 collected from Yaxi unite![]()
图 5 崖西单元样品2017N1 LA-ICP-MS锆石U-Pb年龄谐和年龄图Figure 5. LA-ICP-MS zircon U-Pb concordia diagram of samples 2017N1 form Yaxi pluton表 1 院格庄岩体样品的LA-ICP-MS锆石U-Pb同位素分析结果Table 1. LA-ICP-MS zircon U-Pb isotope analytical results of Yuangezhuang pluton
![]()
图 6 后野单元样品2017N2锆石阴极发光图像Figure 6. Representative CL images and U-Pb ages of zirons from samples 2017N2 collected from Houye unite![]()
图 7 后野单元样品2017N2 LA-ICP-MS锆石U-Pb谐和年龄图Figure 7. LA-ICP-MS zircon U-PbConcordia diagram of samples 2017N2 form Houye unite![]()
图 8 西上寨单元花岗岩样品2017N3锆石阴极发光图像及U-Pb年龄Figure 8. Representative CL images and U-Pb ages of zirons from samples 2017N3 collected from Xishangzhai unite![]()
图 9 西上寨单元花岗岩样品2017N3 LA-ICP-MS锆石U-Pb谐和年龄图Figure 9. LA-ICP-MS zircon U-Pb Concordia diagram of samples 2017N3 from Xishangzai unite5.2 主量元素
岩石主量元素分析结果列于表 2，从表中可以看出，SiO2含量为68.8%~72.28%、Al2O3为14.39%~16.03%、MgO为0.85%~1.77%，K2O/Na2O比值为1.00~1.29，在花岗岩TAS图解（图 10a）上，样品数据主要落于花岗岩区域。岩石总体富碱（K2O+Na2O=7.23~7.79）。在SiO2-（K2O+Na2O）图解上，样品均落在亚碱性岩区域，属高钾钙碱性系列（图 10b、c）。铝饱和指数ASI（A/CNK=0.98~1.02，＜1.1）显示为准铝质花岗岩（图 10d），为富铝（14.39%~16.03%），富硅（66.49%~72.28%，>65%）花岗岩。
表 2 院格庄岩体的主量元素(%)、微量元素和稀土元素(10-6)化学成分Table 2. Chemical compositions of major elements (%), trace elements and REE elements (10-6) of the Yuangezhuang pluton
![]()
图 10 岩石地球化学分类图a—花岗岩TAS分类图（Middlemost, 1994）; b—碱性与亚碱性系列判别图（Irvine and Baragar, 1971）; c—岩石系列SiO2-K2O图解（Peccerillo and Taylor, 1976）; d—岩石类型A/CNK-A/NK图解（Maniar and Piccoli, 1989）Figure 10. Geochemical classification diagrams of granitea-TSA classification diagram of granite(Middlemost, 1994);b-Alkaline and alkaline diagiam for granite(Irvine and Baragar, 1971);c-K2O vs. SiO2 diagram for granite (Peccerillo and Taylor, 1976);d-A/NK vs. A/CNK diagram for granite (Maniar and Piccoli, 1989)5.3 稀土和微量元素
3个侵入岩单元样品的稀土及微量元素分析结果列于表 2，岩石稀土元素总量ΣREE介于131.23×10-6~186.78×10-6，LREE含量为122.01×10-6~175.66×10-6，HREE含量为7.82×10-6~11.12×10-6，δEu=0.84~0.93，δCe=0.87~1.05，LREE/HREE介于13.23~18.84，(La/Yb)N为26.91~35.75显示出明显的轻重稀土分流明显，轻稀土富集的特征。在球粒陨石稀土元素配分图解（图 11）中可以看出，岩石呈明显的轻稀土富集、重稀土分布平坦的右倾型模式，负铕异常不明显，暗示了斜长石的分离结晶作用较弱。
![]()
Figure 11. Chondrite-normalized REE patterns(a, chondrite normalized values after Boynton, 1984) and primitive mantle-normalized traceelement spidergrams (b, primitive mantle normalized values after Sun et al., 1989) for the Yuangezhuang pluton在原始地幔标准化微量元素蛛网图中，样品均表现出Rb、K、U等大离子亲石元素相对富集，Nb、Ta、Ti等高场强元素相对亏损的特征，显示弧岩浆岩特征，其中Nb-Ta-Ti的亏损可能是由于含Ti矿物（如钛铁矿、榍石）的分离结晶作用引起的，Eu的亏损则是由于斜长石的分离结晶作用产生的。利用相容性较差的几个强不活动元素比值可以判断岩浆的源区，样品Rb/Nb介于5.26~9.55，平均7.12，La/Nb介于1.79~4.80，平均2.81，接近大陆壳特征（Weaver, 1991），说明岩体起源于上地壳岩石的部分熔融，而暗色包体的存在，暗示了幔源物质的混入。
6. 讨论
6.1 成因类型、构造环境分析
院格庄岩体各单元产出空间相邻、无明显界线。岩石为似斑状结构，斑晶皆为钾长石，从边缘到中心有小-大-大、密-疏-少的变化特征。闪长质包体主要分布于岩体边缘。该岩体应为同一次岩浆侵入作用下，由于岩浆结晶温度的不同而展现的特征差异，由外而内分别与边缘相、过渡相、中心相相对应。三个侵入单元获得的LA-ICP-MS锆石U-Pb加权平均年龄皆为（118±）Ma，属于伟德山期花岗岩（王世进等，2009）。岩石主量元素显示具有高硅高碱高铝饱和指数（A/NCK>1.1）的特征，同时P2O5（0.1%~0.18%）、MgO（0.85%~1.77%）、TiO2（0.29%~0.41%）含量偏低，显示了I型花岗岩的特点。Rb/Sr值绝大多数为0.05~0.59，数值较低，说明该花岗岩是壳幔混熔成因的深成型花岗岩类，这与岩体中普遍发育的指示地幔来源的暗色包体相符（宋明春等，2003;王对兴等，2019）。在成因类型辨别图中，（Zr+Nb+Ce+Y）值171.42×10-6~242.51×10-6，平均207.44×10-6，小于350×10-6，样品落入OGT范围内（图 12a、b）；而10000Ga/Al值2.22~2.72，平均2.43，小于2.6，岩石样品基本落入I & S型花岗岩内（图 12c~f），与I型花岗岩的特征相对应。
![]()
图 12 Ⅰ型花岗岩成因类型判别a—(Zr+Nb+Ce+Y)-FeO/MgO图解; b—(Zr+Nb+Ce+Y)-(K2O+Na2O)/CaO图解; c—10000Ca/Al-FeO/MgO图解; d—10000Ca/Al-(K2O+Na2O)图解; e—10000Ca/Al-Nb图解; f—10000Ca/Al-Zr图解（据Whalen et al., 1987）Figure 12. Diagram showing the genetic type of the I-type granitea-Diagram of (Zr+Nb+Ce+Y) versus FeO/MgO; b-Diagram of (Zr+Nb+Ce+Y) versus (K2O+Na2O)/CaO; c-Diagram of 10000Ca/Al versus FeO/MgO; d-Diagram of 10000Ca/Al versus (K2O+Na2O); e-Diagram of 10000Ca/Al versus Nb; f-Diagram of 10000Ca/Al versus Zr(after Whalen et al., 1987)Pichavant et al.（1992）指出，磷在强过铝质的熔体中，具有高的溶解度，并随着分异程度的增加含量升高；而在偏铝质或弱过铝质的熔体中，具有很低的溶解度，且随分异程度的增加含量降低。在保证演化样品属于同一岩套的前提下，P2O5和SiO2的关系可以作为鉴别I型和S型花岗岩的一种有效手段。从图 13可以看出，院格庄岩体P2O5含量随着SiO2含量的升高，呈明显降低的趋势，显示出I型花岗岩的特征。
院格庄岩体岩石为准铝质、高钾钙碱性系列的I型花岗岩，类似于活动大陆边缘花岗岩的岩石组合特征。在Nb-Y和(Y +Nb) -Rb构造环境判别图解（图 14）上，样品点分别落入火山弧+碰撞花岗岩、火山弧花岗岩区域，亦反映出与火山弧构造环境的亲缘性。
![]()
图 14 (Y+Nb)-Rb及Y-Nb构造环境判别图解（据Pearce et al., 1984）Figure 14. The diagrams of (Y+Nb) vs. Pb and Y vs. Nb(after Pearce et al., 1984)6.2 岩体与成矿作用的关系
伟德山期侵入岩属胶东中生代第三期岩浆活动，与燕山晚期金及多金属成矿期向对应（丁正江等，2013），近年来发现的铜钼多金属矿床（点）均与此期岩浆活动有关。其中孔辛头矽卡岩型铜钼矿就位于院格庄岩体崖西单元与荆山群陡崖组大理岩接触部位形成的矽卡岩矿体（图 15），辉钼矿样品Re-Os加权平均值年龄为（117.82±0.82）Ma，黄铜矿Re-Os等时线年龄为（118.4±3.2）Ma（李超等，2016），孔辛头铜钼矿床的成矿时代为早白垩世，与院格庄岩体成岩时代一致，表明了院格庄岩体为该矿床的成矿母岩（李超等, 2016）。与院格庄岩体同源的牙山岩体，其SHRIMP锆石测年为（117.7±2.9）Ma（邱连贵等，2008），亦形成于早白垩世，为壳幔混熔型花岗岩（韩宗珠等，2011）。尚家庄斑岩型钼矿床产于牙山岩体中，对尚家庄钼矿的辉钼矿进行Re-Os同位素测年，其成矿时代为（116.4±1.6）Ma，与牙山岩体成岩年龄相一致；辉钼矿中Re含量17.74×10-6~23.72×10-6，平均20.21×10-6，指示成矿物质为壳幔混合来源（李杰等, 2013），说明牙山岩体为尚家庄钼矿的成矿母岩。区内早白垩世岩体周缘可形成矽卡岩型矿床，在岩体内可形成斑岩型矿床，具有铜钼矿成矿专属性。
![]()
图 15 牟平孔辛头铜钼矿床18线地质剖面图（据山东省冶金地质勘探公司第三勘探队, 1979修改）Figure 15. Sketch geological section along No.18 exploration line in the Kong Xintou Cu-Mo deposit, Muping(modified from No.3 Geological Party of Metallurgical Geology Exploration Company Shandong Province, 1979)已有研究表明，伟德山期花岗岩是由壳源酸性岩浆与幔源基性岩浆混合形成的（宋明春等, 2010），花岗岩中大量的闪长质包体记录了岩浆混合的痕迹（张华峰等, 2004）。上述典型矿床成矿物质来源均显示了混合源的岩浆特征，由此分析，伟德山期花岗岩与铜、钼矿床成矿关系密切，而院格庄岩体更是该区铜钼多金属矿的成矿母岩，该期岩体的形成，为铜钼多金属矿成矿提供了热流体和成矿物质。
由前述可知，研究区多金属矿床与院格庄岩体为同一时期岩浆活动的产物，成岩时代限定了成矿时代的下限，即多金属矿床成矿时代不早于118 Ma。白垩纪以来，太平洋板块向欧亚板块俯冲，胶东地区处于弧后拉张的构造环境之下，岩石圈地幔减薄，软流圈上涌，促使壳幔相互作用，形成壳幔混合岩浆并脉动式上侵，而岩浆中携带的大量的成矿物质由于重力分异作用在岩体侵位边部与碳酸盐岩接触交代，在构造有利部位富集成矿，形成矽卡岩型矿床；而处于中部的含矿流体，在上升过程中，随着温压条件的降低，含矿流体内原物理化学平衡条件受到破坏，成矿物质在裂隙等构造薄弱地带沉积成矿，形成斑岩型矿床。
7. 结论
(1) 院格庄岩体3个岩相带LA-ICP-MS锆石U-Pb加权平均年龄分别为（118.10±0.66）Ma、（118.52±0.78）Ma、（118.80±0.67）Ma，为同一期岩浆侵入体。
(2) 院格庄岩体具有高钾钙碱性、准铝质系列（K2O/Na2O值为1.00~1.29、A/CNK=0.98~1.02）特征，岩石轻重稀土元素分馏明显，轻稀土富集，具有弱的负铕异常，微量元素具富集大离子亲石元素，亏损高场强元素特征。为I型花岗岩，属于壳幔混熔型。
(3) 院格庄岩体对其周缘铜钼成矿具有专属性，制约着矿床的产出时代和空间分布。
致谢: 审稿专家为论文提供了宝贵意见；天津地质调查中心田杰鹏博士和何江涛博士在试验测试工作方面给予了帮助，在此一并诚挚感谢！ -
![]()
图 1 胶东主要金和铜钼矿床分布图（据李超等, 2016）
Figure 1. Regional geological map of Jiaodong area showing the distribution of main gold, copper and molybdenum deposits(after Li Chao et al., 2016)
![]()
图 2 院格庄地区区域地质图（据山东省第三地质矿产勘查院, 1990）
Figure 2. Regional geological map of the Yuangezhuang area((afte No.3 Exploration Institute of Geology and Mineral Resources, Shandong Province, 1990)
![]()
图 3 院格庄岩体野外地质特征
a—崖西单元钾长石斑晶密集分布；b—崖西单元闪长质包体；c—后野单元钾长石巨斑呈定向组构；d—后野单元钾长石巨斑；e—西上寨单元岩石野外地质特征；f—西上寨单元钾长石斑晶
Figure 3. Geological characteristics of Yuangezhuang pluton
a-K-feldspar phenocrysts are densely distributed Yaxi unite; b-Diorite inclusions in Yaxi unite granite; c-K-feldspar phenocrysts of Houye unite are in a directional arrangement; d-K-feldspar phenocrysts of Houye unite; e-The geological characteristics of Xishangzai unite granite in the field; f-K-feldspar phenocrysts in Xishangzai unite
![]()
图 4 崖西单元样品2017N1锆石阴极发光图像
Figure 4. Representative CL images and U-Pb ages of zircons from samples 2017N1 collected from Yaxi unite
![]()
图 5 崖西单元样品2017N1 LA-ICP-MS锆石U-Pb年龄谐和年龄图
Figure 5. LA-ICP-MS zircon U-Pb concordia diagram of samples 2017N1 form Yaxi pluton
![]()
图 6 后野单元样品2017N2锆石阴极发光图像
Figure 6. Representative CL images and U-Pb ages of zirons from samples 2017N2 collected from Houye unite
![]()
图 7 后野单元样品2017N2 LA-ICP-MS锆石U-Pb谐和年龄图
Figure 7. LA-ICP-MS zircon U-PbConcordia diagram of samples 2017N2 form Houye unite
![]()
图 8 西上寨单元花岗岩样品2017N3锆石阴极发光图像及U-Pb年龄
Figure 8. Representative CL images and U-Pb ages of zirons from samples 2017N3 collected from Xishangzhai unite
![]()
图 9 西上寨单元花岗岩样品2017N3 LA-ICP-MS锆石U-Pb谐和年龄图
Figure 9. LA-ICP-MS zircon U-Pb Concordia diagram of samples 2017N3 from Xishangzai unite
![]()
图 10 岩石地球化学分类图
a—花岗岩TAS分类图（Middlemost, 1994）; b—碱性与亚碱性系列判别图（Irvine and Baragar, 1971）; c—岩石系列SiO2-K2O图解（Peccerillo and Taylor, 1976）; d—岩石类型A/CNK-A/NK图解（Maniar and Piccoli, 1989）
Figure 10. Geochemical classification diagrams of granite
a-TSA classification diagram of granite(Middlemost, 1994);b-Alkaline and alkaline diagiam for granite(Irvine and Baragar, 1971);c-K2O vs. SiO2 diagram for granite (Peccerillo and Taylor, 1976);d-A/NK vs. A/CNK diagram for granite (Maniar and Piccoli, 1989)
![]()
图 11 院格庄岩体球粒陨石标准化稀土元素分布型式图(a, 球粒陨石标准化值据Boynton, 1984)及原始地幔标准化微量元素蛛网图(b, 原始地幔标准化值据Sun et al., 1989)
Figure 11. Chondrite-normalized REE patterns(a, chondrite normalized values after Boynton, 1984) and primitive mantle-normalized traceelement spidergrams (b, primitive mantle normalized values after Sun et al., 1989) for the Yuangezhuang pluton
![]()
图 12 Ⅰ型花岗岩成因类型判别
a—(Zr+Nb+Ce+Y)-FeO/MgO图解; b—(Zr+Nb+Ce+Y)-(K2O+Na2O)/CaO图解; c—10000Ca/Al-FeO/MgO图解; d—10000Ca/Al-(K2O+Na2O)图解; e—10000Ca/Al-Nb图解; f—10000Ca/Al-Zr图解（据Whalen et al., 1987）
Figure 12. Diagram showing the genetic type of the I-type granite
a-Diagram of (Zr+Nb+Ce+Y) versus FeO/MgO; b-Diagram of (Zr+Nb+Ce+Y) versus (K2O+Na2O)/CaO; c-Diagram of 10000Ca/Al versus FeO/MgO; d-Diagram of 10000Ca/Al versus (K2O+Na2O); e-Diagram of 10000Ca/Al versus Nb; f-Diagram of 10000Ca/Al versus Zr(after Whalen et al., 1987)
![]()
图 14 (Y+Nb)-Rb及Y-Nb构造环境判别图解（据Pearce et al., 1984）
Figure 14. The diagrams of (Y+Nb) vs. Pb and Y vs. Nb(after Pearce et al., 1984)
![]()
图 15 牟平孔辛头铜钼矿床18线地质剖面图
（据山东省冶金地质勘探公司第三勘探队, 1979修改）
Figure 15. Sketch geological section along No.18 exploration line in the Kong Xintou Cu-Mo deposit, Muping
(modified from No.3 Geological Party of Metallurgical Geology Exploration Company Shandong Province, 1979)
表 1 院格庄岩体样品的LA-ICP-MS锆石U-Pb同位素分析结果
Table 1 LA-ICP-MS zircon U-Pb isotope analytical results of Yuangezhuang pluton
下载: 导出CSV
表 2 院格庄岩体的主量元素(%)、微量元素和稀土元素(10-6)化学成分
Table 2 Chemical compositions of major elements (%), trace elements and REE elements (10-6) of the Yuangezhuang pluton
下载: 导出CSV
-
Andersen T. 2002. Correction of common lead in U-Pb analyses that do not report 204Pb[J]. Chemical Geology, (1/2): 59-79.
Boynton W V. 1984. Geochemistry of the rare earth elements: Meteorite studies[C]//Henderson P(ed.). Rare Earth Element Geochemistry[J]. Elsevier, 63-114.
Chappell B W, White A J R. 1992. I-and S-type granites in the Lachlan Fold Belt[J]. Geological Society of America Special Papers, 272: 1-26. http://www.tandfonline.com/servlet/linkout?suffix=CIT0011&dbid=16&doi=10.1080%2F00206814.2017.1377121&key=10.1017%2FS0263593300007720
Cui Kerui, Shi Yangshen. 1994. Offset of rock bodies on bodies on both sides of Muping-Qingdao fault zone in eastern Shandong, China[J]. Journal of Nanjing University(Natural Science edition), 30(4): 648-661. http://en.cnki.com.cn/Article_en/CJFDTotal-NJDZ404.013.htm
Ding Zhengjiang, Sun Fengyue, Liu Fulai, Liu Jianhui, Liu Dianhao, Zhang Pijian, Du Shengxian, Li Bing. 2013. U-Pb dating of zircons from the weideshan molybdenum copper polymetallic deposits in Jiaodong peninsula, China, and its geological significance[J]. Acta Petrologica Sinica, 29(2): 607-618(in Chinese with English abstract). http://www.oalib.com/paper/1474764
Ding Zhengjiang, Sun Fengyue, Liu Fulai, Liu Jianhui, Peng Qiming, Ji Pan, Li Bile, Zhang Pijian. 2015. Mesozoicgeodynamic evolution and metallogenic series of major metaldeposits in Jiaodong Peninsula, China[J]. Acta Petrologica Sinica, 31(10): 3045-3080(in Chinese with English abstract). http://www.zhangqiaokeyan.com/academic-journal-cn_acta-petrologica-sinica_thesis/0201252016378.html
Han Zongzhu, Zhang He, Tian Yuan, Deng Shenggui, Zhang Kaidi, Wang Shujie, Gao Fang. 2011. Geochemisty and granodiorites from the Qixia in Shandong Province[J]. Periodical of Ocean University of China, 41(6): 97-103(in Chinese with English abstract).
Huang Yonghua, You Huichuan, Song Yisheng, Guo Yugui, Huan Wenlin. 2007. Recent activity of faults in Jiaodong Peninsula, Shandong[J]. Technology for Earthquake Disaster Prevention, 2(1): 39-49(in Chinese with English abstract). http://www.researchgate.net/publication/288939454_Recent_Activities_of_Faults_in_Jiaodong_Peninsula_Shandong_Province_China
Irvine T N, Baragar W R A. 1971. A guide to the chemical classification of the common volcanic rocks[J]. Canadian Journal of Earth Sciences, 8: 523-548. doi: 10.1139/e71-055
Li Chao, Pei Haoxiang, Wang Denghong, Zhou Limin, Zhao Hong, Li Xinwei, Qu wenjun, Du Andao. 2016. Age and source constraints for Kongxintou copper molybdenum deposit Shandong from Re-Os isotope in molybdenite and chalcopyrite[J]. Acta Geologica Sinica, 90(2): 240-249(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/ http://search.cnki.net/down/default.aspx?filename=DZXE201602004&dbcode=CJFD&year=2016&dflag=pdfdown
Li Hongkui, Shi Wenge, Li Yifan, Li Luyi, Han Daicheng, Cao Lili, Liu Jimei. 2013. Study on gold mineralization ages in Jiaodong Area, Shandong Province[J]. Gold Science and Technology, 21(3): 1-9(in Chinese with English abstract). http://www.en.cnki.com.cn/Article_en/CJFDTOTAL-HJKJ201303002.htm
Li Jie, Song Mingchun, Wang Meiyun, Li Shiyong, Zhou Mingling, Ni Shijun, Zhang Chengjiang, Ding Zhengjiang, Yue Yaopo. 2013. The molybdenite Re-Os age and genetic analysis of the Shangjiazhuang Mo deposit inJiaodong area[J]. Geology in China, 40(5): 1612-1621(in Chinese with English abstract). http://www.researchgate.net/publication/282371811_The_molybdenite_Re-Os_age_and_genetic_analysis_of_the_Shangjiazhuang_Mo_deposit_in_Jiaodong_area
Liu Yongsheng, Hu Zhaochu, Gao Shan, Gunther Drtlef, Xu Juan, Gao Changgui, Chen Haihong. 2008. In situ analysis of major and trace elements of anhydrous minerals by LA-ICP-MS without applying an internal standard[J]. Chenmical Geology, 257(1/2): 34-43. http://www.sciencedirect.com/science/article/pii/S0009254108003501
Ludwig K R. 2003. User's Manual for Isoplot 3.0: A Geochronological Toolkit for Microsoft Excel[M]. Berkeley Geochronology Center.Special Publication, 4: 1-71.
Maniar P D, Piccoli P M. 1989. Tectonic discrimination of granitoids[J]. Geological Society of America Bulletin, 101: 635-643. doi: 10.1130/0016-7606(1989)101<0635:TDOG>2.3.CO;2
Middlemost E A K. 1994. Naming materials in the magma/igneous rock system[J]. Earth Science Review, 37: 215-224. doi: 10.1016/0012-8252(94)90029-9
Pearce J A, Harris B W, Tindle A G. 1984. Trace element discrimination diagrams for the tectonic interpretation of granitic rocks[J]. Journal of Petrology, 25: 956-983. doi: 10.1093/petrology/25.4.956
Peccerillo A, Taylor S R. 1976. Geochemistry of Eocene calc-alkaline volcanic rocks from the Kastamonu area, northern Turkey[J]. Contribution to Mineralogy and Petrology, 58: 63-81. doi: 10.1007/BF00384745
Pichanant M, Montel J M, Richard L R. 1992. Apatite solubility in peraluminous liquids: Experimental data and an extension of the Harrison-Watson model[J]. Geochimica et Cosmochimica Acta, 56(10): 3855-3861. doi: 10.1016/0016-7037(92)90178-L
Qiu Liangui, Ren Fenglou, Cao Zhongxiang, Zhang Yueqiao. 2008. Late Mesozoicmagmatic activities and their constraints on geotectonics of Jiaodong region[J]. Geotectonica et Metallogenia, 32 (1): 117-123(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/ http://search.cnki.net/down/default.aspx?filename=DGYK200801016&dbcode=CJFD&year=2008&dflag=pdfdown
Song Mingchun, Yin Pihou, Cui Shuxue, Xu Junxiang, Zhou Mingling, Jiang Hongli, Huang Tailing, Jiao Xiumei, Wan Guopu, Cao Chunguo. 2013. Thermal upliftingextensionore-forming theory and its prospecting significance in Jiaodong gold deposit[J]. Shandong Land and Resources, 29(7): 1-12(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-SDDI201307005.htm
Song Mingchun, Song Yingxin, Li Jie, Li Shiyong. 2015. Metallogenic Series of gold and nonferrous metal deposits Related to Cretaceous granites in Eastern Shandong Peninsula, China[J]. Geotectonica et Metallogenia, 39(5): 828-843(in Chinese with English abstract). http://www.cqvip.com/QK/90781X/201505/666468467.html
Song Mingchun, Song Yingxin, Ding Zhengjiang, Li Shiyong. 2018. Jiaodong gold deposits: Essential characteristics and major controversy[J]. Gold Science and Technology, 26(4): 406-422(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTotal-HJKJ201804006.htm
Song Yingxin, Song Mingchun, Ding Zhengjiang, Wei Xufeng, Xu Shaohui, Li Jie, Tan Xianfeng, Li Shiyong, Zhang Zhaolu, Jiao Xiumei, Hu Hong, Cao Jia. 2017. Major advances on deep prospecting in Jiaodong gold ore cluster and its metallogenic characteristics[J]. Gold Science and Technology, 25(3): 4-18(in Chinese with English abstract). http://www.cqvip.com/QK/96123A/201703/672697683.html
Sun S S, McDonough W F. 1989. Chemical and isotopic systematics of oceanic basalts: Implications for mantle composition and processes[C]//Saunders A D, Norry M J(eds.). Magmatism in the Ocean Basins. Geological Society of London Special Publication, 42: 313-345.
Wang Duixing, Guan Qi, GaoWanli, Li Chunlin, Zhang Juquan, Zhao Kaihua. 2019. Zircon U-Pb age, geochemistry and genesis of the Early Cretaceous mafic enclaves and granites in Tiantai area, eastern Zhejiang Province[J]. Geology in China, 46(6): 1512-1529(in Chinese with English abstract).
Wang shijin, Wang laiming, Wan Yusheng, Zhang Chengji, Song Zhiyong, Wang Jinguang. 2009. Age and stage division of intrusive rocks in eastern Shandong Province[J]. Shandong Land and Resources, 25(12): 8-20(in Chinese with English abstract).
Weaver B L. 1991. The origin of ocean island basalt end-member composition: Trace element and isotopic constraints[J]. Earth Planet, 104: 381-397. doi: 10.1016/0012-821X(91)90217-6
Whalen J B, Currie K L, Chappell B W. 1987. A-type granites: Geochemical characteristics, discrimination and petrogenesis[J]. Contributions to Mineralogy and Petrology, 95(4): 407-419. doi: 10.1007/BF00402202
Xie Shiwen, Wang Shijin, Xie Hangqiang, Liu Shoujie, Dong Chunyan, Ma Mingzhu, Ren Peng, Liu Dunyi. 2015. Petrogensis of ca. 2.7 Ga TTG rocks in the Jiaodong terrane, North China craton and its geological implications[J]. Acta Petrologica Sinica, 31(10): 2974-2990(in Chinese with English abstract).
Yang Liqiang, Deng Jun, Wang Zhongliang, Zhang Liang, Guo Linnan, Song Mingchun, Zheng Xiaoli. 2014. Mesozoic gold metallogenic system of the Jiaodong gold province, eastern China[J]. Acta Petrologica Sinica, 30(9): 2447-2467(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-YSXB201409001.htm
Yuan Honglin, Gao Shan, Liu Xiaoming, Li Huiming, Detlef Gunther, Wu Fuyuan. 2004. Accurate U-Pb age and Trace element determinations of zircon by laser ablation-inductively coupled Plasma-Mass spectrometry[J]. Geostandards and Geoanalytical Research, 28(3): 353-370. doi: 10.1111/j.1751-908X.2004.tb00755.x
Zhang Huafeng, Zhai Mingguo, He Zhongfu, Peng Peng, Xu Baoliang. 2004. Petrogenesis and implications of the sodium-rich granites from the Kunyushan complex eastern Shandong Province[J]. Acta Petrologica Sinica, 20(3): 369-380 (in Chinese with English abstract). http://d.wanfangdata.com.cn/Periodical/ysxb98200403001
Zhang Yong, Ren Fenglou, Gong Shuyun, Zhang Xianrong, Lin Manman, Liu Shanshan. 2013. Cretaceous stressfield of the Muping-Jimo fault belt and its implication for tectonic evolution[J]. Marine Geology & Quaternary Geology, 33(2): 79-85(in Chinese with English abstract). http://adsabs.harvard.edu/abs/2014MGQG...33...79Z
Zhang Yueqiao, Li Jinliang, Zhang Tian, Yuan Jiayin. 2007. Late Mesozoickinematic history of the Muping-Jimo fault zone in JiaodongPeninsula, Shandong Province, East China[J]. Geological Review, 53(3): 289-300(in Chinese with English abstract). http://en.cnki.com.cn/Article_en/CJFDTOTAL-DZLP200703000.htm
丁正江, 孙丰月, 刘福来, 刘建辉, 刘殿浩, 张丕建, 杜圣贤, 李兵. 2013. 胶东伟德山地区铜钼多金属矿锆石U-Pb法测年及其地质意义[J]. 岩石学报, 29(2): 607-618. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201302019.htm 丁正江, 孙丰月, 刘福来, 刘建辉, 彭齐鸣, 纪攀, 李碧乐, 张丕建. 2015. 胶东中生代动力学演化及主要金属矿床成矿系列[J]. 岩石学报, 31(10): 3045-3080. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201510011.htm 韩宗珠, 张贺, 田元, 邓声贵, 张凯棣, 王淑杰, 高芳. 2011. 山东省栖霞牙山花岗闪长岩岩石地球化学特征及成因探讨[J]. 中国海洋大学学报, 41(6): 97-103. https://www.cnki.com.cn/Article/CJFDTOTAL-QDHY201106018.htm 黄永华, 尤惠川, 宋毅盛, 郭玉贵, 环文林. 2007. 山东胶东半岛地区断裂最新活动性研究[J]. 震灾防御技术, 2(1): 39-49. doi: 10.3969/j.issn.1673-5722.2007.01.006 李超, 裴浩翔, 王登红, 周利敏, 赵鸿, 李欣慰, 屈文俊, 杜安道. 2016. 山东孔辛头铜钼矿成矿时代及物质来源: 来自黄铜矿、辉钼矿Re-Os同位素证据[J]. 地质学报, 90(2): 240-249. doi: 10.3969/j.issn.0001-5717.2016.02.004 李洪奎, 时文革, 李逸凡, 李璐邑, 韩代成, 曹丽丽, 刘继梅. 2013. 山东胶东地区金矿成矿时代研究[J]. 黄金科学技术, 21(3): 1-9. doi: 10.3969/j.issn.1005-2518.2013.03.001 李杰, 宋明春, 王美云, 李世勇, 周明岭, 倪师军, 张成江, 丁正江, 岳跃破. 2013. 胶东尚家庄钼矿床Re-Os同位素年龄及其地质意义[J]. 中国地质, 40(5): 1612-1621. doi: 10.3969/j.issn.1000-3657.2013.05.024 邱连贵, 任凤楼, 曹忠祥, 张岳桥. 2018. 胶东地区晚中生代岩浆活动及对大地构造的制约[J]. 大地构造与成矿学, 32(1): 117-123. https://www.cnki.com.cn/Article/CJFDTOTAL-DGYK200801016.htm 宋明春, 伊丕厚, 崔书学, 徐军祥, 周明岭, 姜洪利, 黄太岭, 焦秀美, 万国普, 曹春国. 2013. 胶东金矿"热隆-伸展"成矿理论及其找矿意义[J]. 山东国土资源, 29(7): 1-12. doi: 10.3969/j.issn.1672-6979.2013.07.001 宋明春, 宋英昕, 李杰, 李世勇. 2015. 胶东与白垩纪花岗岩有关的金及有色金属矿床成矿系列[J]. 大地构造与成矿学, 39(5): 828-843. https://www.cnki.com.cn/Article/CJFDTOTAL-DGYK201505007.htm 宋明春, 宋英昕, 丁正江, 李世勇. 2018. 胶东金矿床: 基本特征和主要争议[J]. 黄金科学技术, 26(4): 406-422. https://www.cnki.com.cn/Article/CJFDTOTAL-HJKJ201804006.htm 宋英昕, 宋明春, 丁正江, 魏绪峰, 徐韶辉, 李杰, 谭现峰, 李世勇, 张照录, 焦秀美, 胡弘, 曹佳. 2017. 胶东金矿集区深部找矿重要进展及成矿特征[J]. 黄金科学技术, 25(3): 4-18. https://www.cnki.com.cn/Article/CJFDTOTAL-HJKJ201703003.htm 王对兴, 管琪, 高万里, 李春麟, 张聚全, 赵凯华. 2019. 浙东天台地区早白垩世花岗岩及暗色包体锆石U-Pb年龄、地球化学及其成因[J]. 中国地质, 46(6): 1512-1529. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI201906018.htm 王世进, 王来明, 万渝生, 张成基, 宋志勇, 王金光. 2009. 鲁东地区侵入岩形成时代和期次划分[J]. 山东国土资源, 25(12): 8-20. doi: 10.3969/j.issn.1672-6979.2009.12.004 谢士稳, 王世进, 颉颃强, 刘守偈, 董春艳, 马铭株, 任鹏, 刘敦一. 2015. 华北克拉通胶东地区~2.7Ga TTG岩石的成因及地质意义[J]. 岩石学报, 31(10): 2974-2990. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201510008.htm 杨立强, 邓军, 王中亮, 张良, 郭林楠, 宋明春, 郑小礼. 2014. 胶东中生代金成矿系统[J]. 岩石学报, 30(9): 2447-2467. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB201409001.htm 张华锋, 翟明国, 何中甫, 彭澎, 许保良. 2004. 胶东昆嵛山杂岩中高锶花岗岩地球化学成因及其意义[J]. 岩石学报, 20(3): 369-380. https://www.cnki.com.cn/Article/CJFDTOTAL-YSXB200403001.htm 张勇, 任凤楼, 龚淑云, 张现荣, 林曼曼, 刘珊珊. 2013. 牟平-即墨断裂带白垩纪构造应力场及转化机制[J]. 海洋地质与第四纪地质, 33(2): 79-85. https://www.cnki.com.cn/Article/CJFDTOTAL-HYDZ201302013.htm 张岳桥, 李金良, 张田, 袁嘉音. 2007. 胶东半岛牟平-即墨断裂带晚中生代运动学转换历史[J]. 地质论评, 53(3): 289-300. doi: 10.3321/j.issn:0371-5736.2007.03.001 -
期刊类型引用(2)
1. 邹占春，刘俊玉，陈建，李景波，丁正江，孙丽莎，李双飞，唐名鹰，张蕾，王欣然. 胶东西南部三合山岩体岩石成因和构造背景：年代学、地球化学及Sr-Nd-Pb-Hf同位素约束. 中国地质. 2024(04): 1387-1410 . 
本站查看
2. 刘欢，刘晓敏，马玉涛，程鑫，张鲁涛，赵凤臣，王清江，胡培强. 烟台围子山省级自然保护区地质遗迹类型与评价. 山东国土资源. 2022(11): 37-44 . 百度学术
其他类型引用(1)
计量
- 文章访问数: 3002
- HTML全文浏览量: 825
- PDF下载量: 3766
- 被引次数: 3
