Analysis of Ore-forming Conditions and Prediction of Prospecting Targets in Hongqigou Area, Benxi County, Liaoning Province
-
摘要:
通过对辽宁省本溪县红旗沟地区区域典型矿床成矿相关岩体的成矿属性及构造、围岩等成矿条件对比分析,认为研究区内热液石英脉型金(铜)矿床、矽卡岩型、斑岩型铜金钼多金属矿床、与基性岩有关的铜金钴镍金属矿床的极具找矿潜力。运用有利成矿条件信息的赋值统计方法,对研究区进行金铜钼等找矿靶区预测工作。
Abstract:Based on large quantities of Liaoning province red flag trench typical ore deposit mineralization related to regional mineralization properties of rock mass and structure, metallogenic conditions, such as comparative analysis of surrounding rock, in that study hydrothermal quartz vein type gold (copper) ore deposit, skarn type, porphyry copper gold molybdenum polymetallic deposit, related to the rock copper gold cobalt nickel metal deposit prospecting potential. Gold, copper, molybdenum and other prospecting targets in the study area are predicted by the statistical method of assigning the information of favorable metallogenic conditions.
-
青藏高原的形成与演化经历了多个洋盆的闭合以及陆陆碰撞过程,由此形成了高原上多个近EW向延伸的构造缝合带,将青藏高原划分为多个次级地块(Zhu et al., 2011; Kapp et al., 2019)。其中,班公湖–怒江缝合带(以下简称班–怒带)横亘于青藏高原中部,是中生代班公湖–怒江特提斯洋(以下简称班-怒洋)构造演化的残迹,其EW向延伸达2 000 km以上,构成了拉萨地块与羌塘地块之间的地质界线(图1a)。
图 1 青藏高原构造单元划分图(a)、南羌塘地块中—晚侏罗世侵入岩展布图(b)、卡易错地区地质简图(c)1.第四系;2.上三叠统日干配错组;3.古新统—始新统牛堡组;4.中—晚侏罗世花岗岩;5. 构造岩快;6. 断层;7. 角度不整合;8. 湖泊;9. 采样点;10. 锆石U-Pb年龄(本文);11. 锆石U-Pb年龄(引用);JSSZ.金沙江缝合带;LSSZ.龙木措–双湖缝合带;BNSZ.班公湖–怒江缝合带;IYZSZ.印度–雅鲁藏布缝合带;LT.拉萨地块;SQ.南羌塘地块;NQ.北羌塘地块;年龄数据引自Li等(2014)、Liu等(2014)、Wu等(2016)、Sun等(2020)、Yang等(2021)Figure 1. (a) Tectonic subdivision of the Tibetan Plateau, (b) distribution of Middle-Late Jurassic intrusive rocks within Southern Qiangtang block and (c) simplified geological map of the Kayico area, Tibet中生代以来,受班–怒洋俯冲闭合过程的影响,在班–怒带上及其两侧发育了大规模的火山岩浆活动,这些多期次岩浆岩记录了洋盆开合过程的信息,是反演区域构造–岩浆演化的关键,为揭示班–怒洋俯冲闭合过程提供了重要约束(李永飞等, 2005; Zhu et al., 2011, 2016; Pan et al., 2012; Wu et al., 2016, 2019a, 2019b; 刘海永等, 2019; 吴浩等, 2020)。然而,尽管国内外学者先后在青藏高原中部开展了大量的地质研究工作,但是关于区域上中生代多期次、多样性岩浆活动的成因与深部动力学过程尚存有较大的争议,是众多学者关注的热点(Kapp et al., 2007; Pan et al., 2012; Zhu et al., 2016; Fan et al., 2017)。近年来,越来越多的岩浆岩研究资料表明羌塘地块的南缘出露着大规模的中—晚侏罗世岩浆岩(图1b),并显示复杂的地球化学组成,仅在中酸性侵入岩中就先后识别了普通钙碱性I型花岗岩、高分异型花岗岩、富Na埃达克岩与富K埃达克岩等岩石类型(Li et al., 2014; Wu et al., 2018)。然而,羌塘地块南缘中—晚侏罗世多样性岩石类型之间究竟有何成因联系、形成于何种构造环境、反映了怎样的深部动力学过程尚不明确,亟待进一步研究。不仅如此,尽管关于班-怒洋的闭合时限仍存有争议(Kapp et al., 2007; Wu et al., 2019b; Fan et al., 2021),但是当前的研究普遍认为青藏高原中部中-晚侏罗世岩浆岩形成于洋壳俯冲背景,是班–怒洋洋壳俯冲消减引发的弧型岩浆活动(Li et al., 2014; Wu et al., 2016, 2018)。因此,查明羌塘地块南缘中—晚侏罗世岩浆作用过程,不仅对认识班–怒洋构造演化史具有重要的指示意义,同时对探讨俯冲带弧型岩浆起源与演化、壳幔物质循环与交换过程同样具有重要的约束。
藏北卡易错地区发育着大规模的晚侏罗世酸性侵入岩,为揭示班-怒洋俯冲过程、探讨俯冲带复杂的构造–岩浆活动提供了理想的研究对象(图1b)。本次在详细野外地质调查工作的基础上,对藏北卡易错地区出露的花岗闪长岩体进行了系统的岩石学、地质年代学、地球化学和同位素的研究工作,以此确定花岗闪长岩的形成时代与岩石成因,进一步对比区域岩浆岩研究资料,共同约束区域构造–岩浆过程,为揭示班–怒洋俯冲过程、探讨青藏高原早期形成与演化史提供新的约束。
1. 地质背景与岩石学特征
研究表明青藏高原自中生代以来经历了多个地块的闭合、碰撞过程,由此形成了青藏高原上近EW向延伸的多个构造缝合带,并将青藏高原从北至南划分为可可西里–松潘–甘孜地块、北羌塘地块、南羌塘地块、拉萨地块和喜马拉雅地块等多个次级地块(Zhu et al., 2011; Kapp et al., 2019)(图1a)。其中,班-怒带夹持于南羌塘地块和拉萨地块之间,大量的岩浆岩在南羌塘地块的南缘出露着大规模的中-晚侏罗世岩浆岩,该期岩浆岩以中酸性侵入岩为主,具有持续时间长、地球化学组分多样的特征,为探讨青藏高原早期形成与演化过程提供了理想的研究对象。
本次研究区卡易错地区位于日土县NE方向约45 km,大地构造位置处于班-怒带以北、南羌塘地块的南缘。区内构造–地层格架近NW–SE向展布,其研究区西南部主要以上三叠统日干配错组(T3r)灰岩夹砂岩和古新统—始新统牛堡组(E1-2n)砂、砾岩为主;而研究区东北部主要以酸性侵入岩为主。此外,受构造作用的影响,区内出露着大量的灰岩和砾岩的构造岩块(图1c)。前人已经对区内花岗岩体进行了初步的年代学和地球化学的研究工作,研究认为卡易错岩体形成于中—晚侏罗世(168~160 Ma),其地球化学组成指示岩石类型以高分异型花岗岩为主,是古老的变火成岩地壳部分熔融并经历广泛结晶分异作用的产物(Li et al., 2014; Liu et al., 2014)。
本次研究的花岗闪长岩(E 80°6′25″;N 33°28′20″)呈岩株状侵入于构造岩块之中(图2a),出露规模长约为10 m、宽约为3 m,岩石整体呈灰黑色,块状构造,中粗粒花岗结构(图2b),矿物组成以长石、角闪石、石英为主,粒度在0.5~2 mm之间,副矿物有锆石、磷灰石等(图2c、图2d)。
2. 分析方法
本次锆石U-Pb、全岩主微量地球化学与锆石Lu-Hf同位素测试分析工作均在武汉上谱分析科技有限责任公司完成。锆石U-Pb同位素定年和微量元素含量利用LA-ICP-MS同时分析完成,GeolasPro激光剥蚀系统由COMPexPro 102 ArF 193 nm准分子激光器和MicroLas光学系统组成,ICP-MS型号为Agilent 7700e,详细的仪器参数和分析流程见Zong等(2017)和李艳广等(2023)。分析数据的离线处理(包括对样品和空白信号的选择、仪器灵敏度漂移校正、元素含量及U-Pb同位素比值和年龄计算)采用软件ICPMSDataCal(Liu et al., 2008, 2010)完成。锆石样品的U-Pb年龄谐和图绘制和年龄加权平均计算采用Isoplot/Ex_ver3(Ludwig, 2003)完成。全岩主量元素含量利用日本理学PrimusⅡ X射线荧光光谱仪(XRF)分析完成,微量元素含量利用Agilent 7700e ICP-MS分析完成。原位微区锆石Lu-Hf同位素比值测试利用激光剥蚀多接收杯等离子体质谱(LA-MC-ICP-MS)完成。激光剥蚀系统为Geolas HD(Coherent,德国),MC-ICP-MS为Neptune Plus(Thermo Fisher Scientific,德国)。
3. 分析结果
3.1 锆石U-Pb定年与Lu-Hf同位素
花岗闪长岩锆石整体呈无色透明的长柱状,晶形完好,粒径为200~300 μm,长宽比为2∶1~3∶1(图3)。本次共对20颗锆石进行测试分析工作(表1),20颗锆石测点206Pb/238U年龄均集中在162~154 Ma之间。在谐和图上(图3),所有测点都落在谐和线上或附近区域,获得锆石206Pb/238U年龄加权平均值为(158.4±1.8)Ma(MSWD=0.15),这与前人在卡易错花岗岩体中获得的168~160 Ma的年龄信息基本一致,代表了花岗闪长岩的形成时代。
表 1 卡易错花岗闪长岩LA-ICP-MS锆石U-Pb定年分析结果Table 1. LA-ICP-MS U-Pb dating results for zircons of Kayico granodiorites点号 同位素比值(1σ) 年龄比值(Ma) 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 207Pb/206Pb 1σ 207Pb/235U 1σ 206Pb/238U 1σ 18T12-01 0.04954 0.00219 0.17017 0.00781 0.02491 0.00066 173 59 160 7 159 4 18T12-02 0.04929 0.00266 0.16824 0.00927 0.02475 0.00067 162 78 158 8 158 4 18T12-03 0.04965 0.00191 0.16944 0.00690 0.02475 0.00064 179 50 159 6 158 4 18T12-04 0.04989 0.00233 0.16974 0.00824 0.02467 0.00064 190 66 159 7 157 4 18T12-05 0.04939 0.00438 0.17064 0.01519 0.02505 0.00071 166 147 160 13 159 4 18T12-06 0.05092 0.00237 0.16985 0.00810 0.02419 0.00066 237 61 159 7 154 4 18T12-07 0.04885 0.00211 0.16623 0.00740 0.02467 0.00067 141 56 156 6 157 4 18T12-08 0.04955 0.00218 0.17037 0.00777 0.02493 0.00066 174 59 160 7 159 4 18T12-09 0.04888 0.00299 0.16655 0.01042 0.02471 0.00064 142 95 156 9 157 4 18T12-10 0.05012 0.00238 0.17289 0.00852 0.02501 0.00065 201 67 162 7 159 4 18T12-11 0.04908 0.00453 0.16867 0.01561 0.02492 0.00071 152 154 158 14 159 4 18T12-12 0.04889 0.00190 0.16671 0.00681 0.02472 0.00065 143 50 157 6 157 4 18T12-13 0.04940 0.00280 0.17026 0.00980 0.02499 0.00069 167 83 160 9 159 4 18T12-14 0.04970 0.00173 0.16964 0.00637 0.02475 0.00064 181 43 159 6 158 4 18T12-15 0.04905 0.00219 0.16888 0.00775 0.02496 0.00068 150 59 158 7 159 4 18T12-16 0.04950 0.00233 0.17051 0.00831 0.02498 0.00066 172 66 160 7 159 4 18T12-17 0.05017 0.00176 0.17246 0.00651 0.02492 0.00064 203 44 162 6 159 4 18T12-18 0.04892 0.00463 0.16803 0.01594 0.02490 0.00071 144 159 158 14 159 4 18T12-19 0.05000 0.00226 0.17212 0.00808 0.02496 0.00066 195 62 161 7 159 4 18T12-20 0.04862 0.00207 0.17018 0.00751 0.02538 0.00068 130 56 160 7 162 4 此外,对8颗获得谐和年龄的锆石进行原位Lu-Hf同位素测试(表2),8颗锆石测点的176Yb/177Hf在
0.012865 ~0.037412 之间,176Lu/177Hf在0.000427 ~0.001145 之间,表明所测锆石放射成因Hf的积累很少(吴福元等, 2007),测定的176Hf/177Hf(0.282602 ~0.282669 )值可以用来代替锆石的初始176Hf/177Hf值。8个测点Hf同位素地壳模式年龄TDMC为1090 ~1222 Ma,对应的εHf(t)值变化范围在−2.66~−0.27之间。表 2 卡易错花岗闪长岩锆石Lu-Hf同位素组成Table 2. Lu-Hf isotopes of zircons from the Kayico granodiorites.点号 年龄 (Ma) 176Hf/177Hf 1σ 176Lu/177Hf 1σ 176Yb/177Hf 1σ εHf(0) 1σ εHf(t) 1σ TDM1 TDMC fLu/Hf 18T12-01 159 0.282608 0.000008 0.000562 0.000003 0.017391 0.000117 −5.81 0.59 −2.38 0.60 901 1207 −0.98 18T12-02 158 0.282602 0.000009 0.000996 0.000005 0.031211 0.000271 −6.03 0.61 −2.66 0.62 921 1222 −0.97 18T12-03 158 0.282609 0.000010 0.000842 0.000026 0.026596 0.000831 −5.77 0.61 −2.39 0.62 907 1207 −0.97 18T12-04 154 0.282618 0.000010 0.000867 0.000013 0.027469 0.000376 −5.46 0.62 −2.17 0.63 895 1191 −0.97 18T12-05 159 0.282641 0.000009 0.000427 0.000002 0.012865 0.000067 −4.65 0.61 −1.20 0.62 853 1142 −0.99 18T12-06 157 0.282643 0.000009 0.000608 0.000005 0.019402 0.000131 −4.57 0.60 −1.18 0.61 854 1139 −0.98 18T12-07 159 0.282631 0.000008 0.000532 0.000001 0.016431 0.000077 −4.98 0.59 −1.55 0.60 868 1161 −0.98 18T12-08 159 0.282669 0.000009 0.001145 0.000043 0.037412 0.001346 −3.64 0.60 −0.27 0.61 829 1090 −0.97 3.2 地球化学
本次共采集4件花岗闪长岩样品进行全岩主微量元素地球化学分析工作,分析结果见表3。4件样品的地球化学组成相对均一,其SiO2含量为62.6%~65.2%,Al2O3含量为15.9%~16.6%,TiO2含量为0.68%~0.81%,TFe2O3含量为4.53%~5.37%,MgO含量为2.06%~2.34%,全碱(Na2O+K2O)含量为5.40%~5.56%,Na2O/K2O值为0.94~1.10,Mg#为50~52。在岩石类型判别图解中(图4a),样品均在花岗闪长岩区域;在K2O-SiO2图解中(图4b),样品显示中钾–高钾钙碱性的特征;在A/NK-A/CNK图解中(图3c),样品整体显示弱过铝质的特征(A/CNK=1.02~1.05)。以上主量元素地球化学组成表明样品整体显示弱过铝质中钾-高钾钙碱性花岗闪长岩的特征。
表 3 卡易错花岗闪长岩全岩主量(%)和微量(10−6)元素分析结果Table 3. Whole-rock major (%) and trace (10−6) element contents of Kayico granodiorites元素 T12h1 T12h2 T12h3 T12h4 元素 T12h1 T12h2 T12h3 T12h4 SiO2 63.6 65.2 64.5 62.6 Zr 248 230 244 245 TiO2 0.78 0.68 0.71 0.81 Nb 13.8 13.1 13.7 14.1 Al2O3 16.1 15.9 16.2 16.6 Sn 3.58 5.26 4.55 4.30 TFe2O3 5.37 4.53 4.82 5.12 Cs 17.5 21.1 14.6 20.0 MnO 0.08 0.06 0.07 0.07 Ba 457 383 401 428 MgO 2.34 2.06 2.07 2.33 La 29.8 32.7 32.1 30.2 CaO 4.48 4.44 4.38 4.83 Ce 66.7 68.5 63.0 61.6 Na2O 2.70 2.74 2.83 2.74 Pr 6.97 7.45 7.02 6.90 K2O 2.86 2.78 2.57 2.77 Nd 26.37 28.1 27.3 25.4 P2O5 0.15 0.13 0.14 0.15 Sm 5.26 5.68 5.75 5.24 LOI 0.99 1.11 1.30 1.30 Eu 1.11 1.05 1.14 1.16 SUM 99.5 99.7 99.5 99.3 Gd 5.08 4.65 4.55 4.66 Li 56.7 49.2 63.6 54.6 Tb 0.83 0.83 0.78 0.74 Be 1.97 2.31 2.18 2.07 Dy 4.95 5.01 4.54 4.49 Sc 14.9 12.7 13.8 14.8 Ho 1.00 1.01 0.91 0.92 V 74.0 61.3 67.1 74.1 Er 2.83 2.87 2.72 2.60 Cr 40.0 32.8 37.6 37.7 Tm 0.45 0.43 0.40 0.41 Co 15.4 11.5 12.9 12.8 Yb 2.70 2.78 2.46 2.42 Ni 16.6 14.3 15.2 15.5 Lu 0.41 0.44 0.37 0.37 Cu 104 45.7 44.1 45.7 Hf 6.64 6.24 6.79 6.29 Zn 56.6 43.1 45.7 51.6 Ta 1.02 1.10 1.05 0.98 Ga 19.2 18.8 19.1 19.5 Tl 1.33 1.40 1.13 1.34 Rb 170 185 169 178 Pb 8.40 9.57 7.76 10.2 Sr 188 181 194 203 Th 14.4 16.9 15.5 12.4 Y 28.2 28.0 26.8 24.8 U 2.18 2.88 2.52 2.23 图 4 卡易错花岗岩TAS图解(Middlemost, 1994)(a)、K2O-SiO2图解(Le Maitre et al., 1989; Rickwood, 1989)(b)和A/NK-A/CNK图解(c)(Shand, 1943)Figure 4. (a) TAS classification diagram, (b) K2O vs. SiO2 diagram and (c) A/NK vs. A/CNK diagram of Kayico granitic rocks在球粒陨石标准化稀土元素配分曲线中,样品呈轻稀土元素富集、重稀土元素亏损的右倾模式([La/Yb]N=7.90~9.34),同时具有不同程度的Eu负异常(Eu/Eu*=0.62~0.72)(图5a)。在原始地幔标准化蛛网图中,样品具有明显Nb、Ta等高场强元素以及Ba、Sr、Eu等大离子亲石元素的亏损(图5b)。
图 5 岩石球粒陨石标准化稀土元素模式图(a)和原始地幔标准化微量元素蛛网图(b)(标准化值引自Sun 等1989)Figure 5. (a) Chondrite-normalized rare earth element and (b) primitive-mantle–normalized multi-element patterns4. 讨论
4.1 中-晚侏罗世岩浆活动
近年来,国内外学者已经对青藏高原中部出露的岩浆岩开展了系统的研究工作,报道了大量的年代学和地球化学数据(Li et al., 2014; Liu et al., 2014; Wu et al., 2016; 刘海永等, 2022)。越来越多的研究表明在南羌塘地块南缘的卡易错、材玛、青草山、改则、荣玛、高保约一带存在一期大规模的中—晚侏罗世岩浆活动(约165~150 Ma),这些岩浆岩整体呈带状近平行于班–怒带出露,EW向延伸近上千千米,岩石类型整体以中酸性侵入岩为主,具有持续时间长、分布范围广的特征(图1b)。前人的研究在卡易错花岗岩体中获得了168~160 Ma的锆石U-Pb年龄信息(Liu et al., 2014; Li et al., 2014),本次锆石U-Pb定年工作在卡易错花岗闪长岩中获得了158 Ma的年龄,表明区内花岗岩形成于中—晚侏罗世之交,与南羌塘地块上大规模发育的中酸性侵入岩形成时代相近,应该是区域上同一期构造–岩浆活动的产物。
此外,Wu等(2018)对南羌塘地块上发育的中—晚侏罗世侵入岩进行了系统的地球化学对比研究,研究发现该期侵入岩整体显示复杂的地球化学特征,根据岩石的主、微量元素含量和同位素地球化学特征,可以划分为普通钙碱性I型花岗岩、高分异花岗岩、富Na埃达克岩与富K埃达克岩等不同的岩石类型,反应了区域上花岗质岩石复杂的岩浆源区与成岩过程。而根据现有地球化学资料,卡易错花岗岩可以划分为普通钙碱性I型花岗岩和高分异型花岗岩两类(Li et al., 2014; Liu et al., 2014),然而二者之间有何成因联系尚不明确。
4.2 构造背景
大量的岩浆岩研究资料表明青藏高原中部中生代岩浆作用整体沿班–怒带展布,主要发育在混杂带上及其两侧的南羌塘地块南缘和拉萨地块北缘,显示与与班–怒洋密切的时空联系,系统的研究工作已经初步建立中生代多期次岩浆作用与班–怒洋俯冲闭合过程之间的成因联系(Zhu et al., 2011, 2016; Kapp et al., 2019; Wu et al., 2019a, 2019b)。尽管对于班–怒洋的闭合时限一直存在争论,然而现有的争议均认为班–怒洋闭合与拉萨–羌塘地块陆陆碰撞应不早于白垩纪,而青藏高原中部侏罗纪构造演化主要受班–怒洋俯冲消减作用的影响(Liu et al., 2022)。此外,班–怒带上蛇绿混杂岩系统的年代学和地球化学的研究同样表明班–怒洋洋盆在中—晚侏罗世仍存有一定的规模(范建军等, 2019; 李志军等, 2019; 唐跃等, 2021)。
卡易错花岗闪长岩显著的Nb、Ta等元素亏损(图5),显示与俯冲成因的弧型岩浆岩相似的地球化学特征。近年来,不同学者对南羌塘地块上发育的中—晚侏罗世岩浆岩开展了大量的研究工作,陆续取得了众多的研究进展(Liu et al., 2014; Li et al., 2016; Wu et al., 2016)。其中,中酸性侵入岩整体显示陆缘弧的地球化学特征,而洋壳熔融成因的埃达克质花岗岩的识别则为深部俯冲洋壳的存在提供了最直接的岩石学证据(Li et al., 2016)。此外,区域上还陆续报道了中—晚侏罗世俯冲成因的钙碱性弧型安山岩与OIB型辉绿岩(李小波等, 2015; 董宇超等, 2016; Li et al., 2016),如此复杂多样的岩石组合反映了俯冲带上复杂的源区物质组成与循环过程。不仅如此,最近的研究提出班–怒洋的初始俯冲起始于晚三叠世末—早侏罗世(Qian et al., 2020; Liu et al., 2022),进一步表明班–怒洋中—晚侏罗世应处于洋壳俯冲背景。综上所述,南羌塘地块南缘出露的中—晚侏罗世岩浆岩应该形成于洋壳俯冲背景,是班-怒洋北向俯冲至南羌塘地块之下引发的弧型岩浆活动。
4.3 岩石类型与成因
花岗岩的类型划分与岩石成因长期以来一直以来是众多地质学家关注的热点研究问题,其中Chappell等(1974, 1992)根据花岗岩的岩浆源区物质组成和成岩构造环境的差异将花岗岩划分为I、S、M、A型4类,该分类方式已经被广泛运用于花岗岩成因与演化的研究工作(王亮等, 2022; 孙巍等, 2022)。Li等(2014)和Liu等(2014)对区内花岗岩岩石开展了年代学和地球化学的研究工作,识别了普通钙碱性I型和高分异型两类花岗岩。笔者在新获得的卡易错花岗闪长岩研究资料的基础上,进一步收集整理了前人已报道的卡易错花岗岩体的数据资料,以此准确约束卡易错地区花岗质岩石的岩石类型与成因。
卡易错花岗岩整体具有低的Zr+Ce+Nb+Y含量,然而根据其(Na2O+K2O)/CaO和FeOt/MgO值的不同,区内岩石可以划分为高分异与未分异两组(图6a、图6b )。不同于前人报道的高分异型花岗岩,本次采集的花岗闪长岩样品显示低的(Na2O+K2O)/CaO和FeOt/MgO值,具有未分异花岗岩的特征。在P2O5-SiO2和Th-Rb图解中,花岗闪长岩样品均显示出I型花岗岩的演化趋势(图6c、图6d)。结合花岗闪长岩弱过铝质的特征(图4c),卡易错花岗闪长岩应该属于普通钙碱性I型花岗岩。I型花岗岩一般认为起源于变火成岩下地壳部分熔融或者幔源物质对变沉积岩下地壳的改造(Petford et al., 1996; Chappell et al., 2001; Li et al., 2007),研究区内乃至整个南羌塘地块南缘均未发现大规模幔源岩浆活动成因的基性岩,结合花岗闪长岩相对均一的锆石原位Hf同位素组成(εHf(t)=−2.66~−0.27),卡易错花岗闪长岩难以解释为幔源物质与变沉积岩下地壳熔体混合的产物。同时,花岗闪长岩具有高的CaO/Na2O以及低的Al2O3/TiO2和Rb/Ba、Rb/Sr值,进一步指示着其起源于玄武质火成岩下地壳的部分熔融(图7a、 图7b)。此外,花岗闪长岩中不同程度的Eu、Sr、Ba等元素的负异常一般认为是成岩过程中存在着长石类矿物的结晶分异,而Nb、Ta等元素的亏损则和金红石/榍石的结晶分离相关。综上所述,卡易错花岗闪长岩应该是南羌塘地块之下古老的变火成岩下地壳熔体经历一定结晶分异作用的产物。
图 6 卡易错花岗岩(Na2O/K2O)/CaO-Zr+Ce+Nb+Y图解(a)、FeOt/MgO-Zr+Ce+Nb+Y图解(b)(Whalen et al., 1987)、P2O5-SiO2图解(c)和Th-Rb图解(d)(Li et al., 2007)Figure 6. Geochemical classification diagrams of Kayico granitic rocks: (a)(Na2O/K2O)/CaO vs. Zr+Ce+Nb+Y diagram, (b) FeOt/MgO vs. Zr+Ce+Nb+Y diagram, (c) P2O5 vs. SiO2diagram, (d) Th vs. Rb diagram of Kayico granitic rocks图 7 卡易错花岗岩CaO/Na2O-Al2O3/TiO2图解(a)、Rb/Ba-Rb/Sr图解(b)(Sylvester, 1998)、 Rb/Sr-Sr图解(c)、Ba-Sr图解(d)(Rollinson, 1993)Amp.闪石;Bi.黑云母;Kfs.钾长石;Pl.斜长石Figure 7. (a) CaO/Na2O vs. Al2O3/TiO2 diagram, (b) Rb/Ba vs. Rb/Sr diagram, (c) Rb/Sr vs. Sr diagram, (d) Ba vs. Sr diagram of Kayico granitic rocks卡易错高分异型花岗岩具有与钙碱性I型花岗岩形成时代一致,在空间上密切共生,并相似的锆石Hf同位素特征,表明二者具有相似的岩浆源区物质组成(Li et al., 2014; Liu et al., 2014)。然而,高分异型花岗岩显示强烈的亏损Eu、Sr、Ba等元素(图5),指示着岩浆在侵位过程中经历了强烈的结晶分异作用。在Rb/Sr-Sr和Ba-Sr图解中,卡易错两类花岗岩之间显示出明显的长石类矿物分离结晶趋势(图7c、图7d),进一步表明卡易错不同类型花岗岩地球化学的差异是后期岩浆侵位过程中经历不同演化过程的产物。近年来,晶粥体模型(MUSH)被广泛应用于解释同期共生花岗闪长岩与高分异型花岗岩之间的成因联系(Bachmann et al., 2004; Hildreth, 2004),该模型提出花岗质岩浆在浅层岩浆房中发生矿物结晶时,形成一种晶体与液体共存的晶粥体,其外围的物质结晶形成富含斑晶的花岗质岩石,而中心的残余岩浆则形成高分异的高硅花岗岩(Wu et al., 2017)。晶粥体模型无疑为卡易错地区空间上相伴生、时代上相一致的普通钙碱性I型花岗岩和高分异型花岗岩提供了合理的成因解释。因此,研究认为,在中—晚侏罗世班–怒洋持续北向俯冲过程中,底侵的幔源玄武质岩浆诱发南羌塘地块之下古老的变火成岩下地壳物质发生重熔,形成的熔体上升侵位、并在浅层岩浆房内发生显著低压结晶分异作用,其外围的晶粥体与内部的熔体分别冷凝形成了卡易错地区普通钙碱性I型和高分异型两类花岗质岩石(图8)。
5. 结论
(1)锆石U-Pb年龄表明,卡易错花岗闪长岩形成于158 Ma,与前人在卡易错岩体中获得的168~160 Ma的年龄相一致,是南羌塘地块南缘中—晚侏罗世构造-岩浆活动的产物。
(2)卡易错花岗闪长岩具有钙碱性I型花岗岩的地球化学特征,以及相对富集的锆石Hf同位素组成,是班-怒洋中-晚侏罗世北向俯冲背景下南羌塘地块古老的变火成岩下地壳熔融的产物。
(3)结合现有研究资料,晶粥体模式为卡易错地区钙碱性I型与高分异型两类花岗岩提供了合理的成因解释,二者是壳源熔体在浅层岩浆房内经历结晶分异后不同端元冷凝的产物。
致谢:野外工作得到西藏自治区地质调查院刘海永博士和吉林大学罗安波博士等人的帮助,审稿专家对稿件提出的宝贵意见对论文质量提高至关重要,在此一并致以衷心的感谢。
-
图 1 研究区构造位置图(a)和区域地质图(b)
1. 白垩系小岭组;2.侏罗系;3.二叠系;4.石炭系;5.奥陶系;6.寒武系;7. 南华系桥头组;8. 青白口系南芬组;9. 青白口系钓鱼台组;10. 青白口系;11. 下元古界含铁碎屑岩建造;12. 下元古界含有色贵金属碎屑岩-碳酸盐岩建造;13. 下元古界含硼火山岩建造;14. 新太古界含铜金绿岩建造;15. 中太古界含铁建造;16. 早白垩世花岗斑岩;17. 早白垩世碱长花岗岩;18. 早白垩世似斑状正长花岗岩;19. 早白垩世闪长岩;20. 早白垩世二长花岗岩;21. 中侏罗世花岗闪长岩;22. 晚三叠世粗面岩;23. 古元古代片麻状花岗岩;24. 太古代混合岩;25. 太古代混合花岗岩;26. 实测断层;27. 推测隐伏断层;28. 整合地质界线;29. 不整合地质界线;30. 火山机构;31. 铁矿床;32. 铜矿床;33. 铜钼矿床;34. 铜多金属矿床;35. 金矿(床);36. 铜金矿;37. 乡镇驻地;38. 研究区范围
Figure 1. (a) Structural location map and (b) regional geological map of the study area
图 2 岩体成因类型Na2O-K2O图解(a)(据Collins et al.,1982)和成矿专属性判别Rb/Sr-SiO2图解(b)(据Blevin et al.,1995)
Figure 2. (a) Diagram of Na2O-K2O and (b) diagram of Rb/Sr-SiO2 discrimination of metallogenic specificity
表 1 典型矿床特征及有利成矿条件统计表
Table 1 Statistical table of typical deposit characteristics and favorable metallogenic conditions
典型矿床名称 成矿及
伴生元素成因类型 岩浆岩 控矿构造 容矿构造 赋矿围岩 矿源层 化探异常 航磁异常 黄柏峪铜钼矿床 主Cu、Mo; 矽卡岩型 1.γπJ2 1.EW∧NE断裂;
2.向斜构造1.岩体与围岩接触带;2.EW向断裂蚀变带 1.Qbn泥灰岩、Єj灰岩、Om灰岩;2.矽卡岩 1.Ar2;2.γPt1 发育1∶20万水系沉积CuPbZnAgAs综合异常 无 马鹿沟铜多金属
矿床主Cu、Mo、Au、 W、 Tfe;伴生V、Ga、In、Ge、Np 矽卡岩型 1.ξγK1;2.κγπK1 1.EW∧SN∧NE∧NW向断裂;2.倒转背斜构造 1.岩体与围岩接触带;2.近SN向层间裂隙 1.Єz灰岩、Om;2.矽卡岩 1.ξγK1;2.κγK1 发育1∶20万水系沉积物铜铅锌钼组合异常 无 套袖峪铜铅锌多金属矿床 主Cu、Pb、Zn、Sn、S,伴生银、 镉、铟等 矽卡岩型 1.ηγK1(隐伏);2.δK1 1.EW∧SN∧NW断裂;2.背斜构造 1.NW向及SN向的次级断裂;2.层间裂隙、节理 1.Qbn泥灰岩、Єj角砾状灰岩 1.ηγK1、2.δK1 发育1∶20万水系沉积物铜铅锌钼组合异常 无 白水金
矿床主Au 热液石英脉+破碎蚀变岩型 1.ηγK1;2.δK1 1.EW∧SN∧NW断裂 1.EW向断层破碎带;2.SN向张性节理裂隙 1.Qbd石英岩、角砾岩;2.Pt12大理岩,片岩 1.ηγK1、2.δK1 发育1∶20万水系沉积物AuCuPb异常 无 红旗沟
金矿床主Au,伴生Cu、Pb、Ag 热液石英脉型 1.γδJ2 1.NNE∧NNW断裂 1.SN向节理;2.NW向节理;3.NE向节理 1.(γδJ2) 1.γδJ2 发育1∶20万水系沉积物金铜钼组合异常 有正
异常沙金沟
铜金矿主Cu、Au;伴生Ag、Zn、Bi、Sr 热液石英脉型 1.δK1;2.γPt1;3.δ、χ 1.NW∧EW向2. 组断裂构造交接处 1.NNW向张性裂隙;2.NE向张性裂隙 1.Qbd石英岩;2.Qbn页岩、泥灰岩 1.γPt1;2.δK1 发育1∶20万水系沉积物金铜钼组合异常 有正
异常表 2 岩体全岩地球化学数据(主量元素:%;稀土和微量元素: 10−6)
Table 2 Whole-rock geochemical data of rock mass (Major elements: %; REE and Trace elements: 10−6)
岩体名称代号 样号 岩性 Na2O MgO Al2O3 SiO2 P2O5 K2O CaO TiO2 MnO 柳树底下岩体
γπJ2QP37S1 花岗斑岩 0.48 0.29 10.05 80.5 0.07 6.65 0.11 0.12 0.12 QP37S3 花岗斑岩 0.42 0.61 12.47 75.52 0.11 7.95 0.04 0.28 0.28 马鹿沟ξγK1 Pm105XWG13 正长花岗岩 4.36 0.7 13.88 69.9 0.12 4.18 0.72 0.47 0.091 Pm105XWG14 正长花岗岩 4.78 0.6 14.74 68.76 0.13 4.61 0.66 0.45 0.094 荫士王沟
κγπ K1Pm506XWG2 碱长花岗斑岩 4.86 0.75 15.32 68.44 0.11 4.39 0.49 0.38 0.09 Pm506XWG3 碱长花岗斑岩 3.74 0.26 12.92 75.52 0.01 4.53 0.03 0.10 0.06 北大山ηγK1 QP43S1 二长花岗岩 5.24 0.06 5.24 74.96 0.06 5.47 0.06 0.008 0.008 QP43S3 二长花岗岩 4.64 0.87 4.64 70.46 0.08 4.59 0.08 0.039 0.039 QP39S1 二长花岗岩 4.7 0.98 4.7 71.02 0.15 4.89 0.15 0.036 0.036 红旗沟γδJ2 C-G6 花岗闪长岩 4.76 0.34 16.63 69.41 0.06 2.32 3.00 0.17 0.08 QP34S2 花岗闪长岩 5.61 1.22 17.60 65.04 0.25 2.40 3.21 0.33 0.06 QP34S2-3 花岗闪长岩 5.31 1.03 17.36 66.82 0.23 2.88 3.65 0.30 0.06 沙金沟γPt1 PM511G1 片二长花岗岩 3.85 0.54 12.97 75.8 0.02 4.05 0.17 0.08 0.03 QP30D2 二长花岗岩 3.46 0.37 13.54 73.84 0.07 4.14 1.28 0.11 0.03 白石砬子δK1 PM007G1δK1 角闪石闪长岩 2.84 6.08 17 46.95 0.05 0.69 12.6 1.525 0.124 PM007G2δK1 角闪石闪长岩 3.88 3.00 17.19 52.67 0.5 3.61 6.84 1.726 0.102 柳树底下岩体
γπJ2QP37S1 花岗斑岩 0.34 1.72 1.72 968.8 2.65 5.75 5.382 1.77 145.7 QP37S3 花岗斑岩 0.34 1.87 1.87 890.9 3.79 7.00 7.708 3.78 218.4 马鹿沟ξγK1 Pm105XWG13 正长花岗岩 2.77 1.8 0.9 609.39 5.24 12.52 41.307 3.72 139.05 Pm105XWG14 正长花岗岩 2.65 1.56 0.8 905.96 3.57 12.20 34.280 3.64 148.5 荫士王沟
κγπ K1Pm506XWG2 碱长花岗斑岩 4.07 2.21 0.91 737.61 4.32 8.53 33.950 4.01 143.15 Pm506XWG3 碱长花岗斑岩 2.25 1.01 0.52 169.29 1.86 3.19 58.362 0.17 220.15 北大山ηγK1 QP43S1 二长花岗岩 1.07 0.14 1.07 485.9 1.58 6.21 20.690 2.52 216.2 QP43S3 二长花岗岩 1.07 0.97 1.07 917.8 3.79 12.49 28.130 6.01 250.1 QP39S1 二长花岗岩 0.93 1.9 0.93 1127 4.90 5.50 0.050 5.20 121.9 红旗沟γδJ2 C-G6 花岗闪长岩 2.62 1.58 0.35 1053 1.72 4.38 8.050 3.85 55.3 QP34S2 花岗闪长岩 1.45 1.29 2.06 873.2 4.61 4.92 6.013 2.62 40.1 QP34S2-3 花岗闪长岩 1.28 1.26 0.35 878.3 5.40 7.18 5.716 3.33 29.9 沙金沟γPt1 PM511G1 片二长花岗岩 2.02 1.31 0.37 515.64 1.34 7.13 5.620 4.35 105.44 QP30D2 二长花岗岩 0.73 1.33 1.04 1209.2 15.80 0.60 6.500 1.00 123.9 白石砬子δK1 PM007G1δK1 角闪石闪长岩 10.776 6.76 1.18 1601 18.00 27.06 21.700 32.52 77.6 PM007G2δK1 角闪石闪长岩 9.774 5.25 0.57 109 1.72 7.34 56.740 1.80 199 岩体名称代号 样号 岩性 Sr V Zr Hf Ta Li Th U Rb/Sr 数据来源 柳树底下岩体γπJ2 QP37S1 花岗斑岩 150.7 34.58 105.4 4.04 4.04 4.04 1.15 0.72 1.03 ① QP37S3 花岗斑岩 152.3 54.93 53.6 2.55 2.55 2.55 1.66 0.93 0.70 ① 马鹿沟
ξγK1Pm105XWG13 正长花岗岩 186.1 20.98 382.1 10.01 1.89 20.26 15.83 2.29 0.75 ② Pm105XWG14 正长花岗岩 188.3 21.26 329.1 10.56 1.25 16.52 14.08 2.22 0.79 ② 荫士王沟
κγπ K1Pm506XWG2 碱长花岗斑岩 229.1 23.32 303.5 8.62 1.91 17.27 5.405 2.10 0.62 ② Pm506XWG3 碱长花岗斑岩 23.1 4.37 226.6 9.46 4.11 7.53 7.013 4.52 9.52 ② 北大山
ηγK1QP43S1 二长花岗岩 201.6 10.87 64.5 3.10 3.14 15.90 27.79 3.99 1.07 ① QP43S3 二长花岗岩 446.3 20.55 160.1 6.14 4.33 46.58 22.43 3.94 0.56 ① QP39S1 二长花岗岩 540.6 24.18 89.6 3.11 0.83 0.32 3.82 1.18 0.23 ① 红旗沟
γδJ2C-G6 花岗闪长岩 1015.4 9.03 124 3.25 0.45 19.23 0.78 0.78 0.05 ② QP34S2 花岗闪长岩 1162.0 69.50 119.1 4.53 1.35 9.98 1.74 0.56 0.03 ① QP34S2-3 花岗闪长岩 1212.0 60.88 6.4 0.55 1.31 7.34 0.84 0.47 0.02 ① 沙金沟γPt1 PM511G1 片二长花岗岩 104.5 9.27 81.78 2.22 0.62 4.12 11.91 6.58 1.01 ② QP30D2 二长花岗岩 488.8 13.70 138.7 12.50 8.14 20.00 10.2 2.00 0.25 ① 白石砬子
δK1PM007G1δK1 角闪石闪长岩 988.1 146.67 234 4.35 1.44 15.76 4.17 1.45 0.08 ② PM007G2δK1 角闪石闪长岩 28.2 5.01 257 6.90 4.49 13.66 9.54 2.74 7.06 ② 注:①.数据来源于辽宁地质勘查院(2014);②.数据来源于沈阳地质调查中心(2015)。 表 3 成矿要素权重赋值表
Table 3 Weight assignment table of ore-forming elements
要素名称 要素属性 热液型金铜矿床相关要素权重赋值 矽卡岩型铜金矿床相关要素权重赋值 斑岩型铜钼金矿床相关要素权重赋值 a1 岩体γδJ2 1 1 1 a2 岩体δK1 1 1 1 a3 岩体ξγK1 0 1 0.5 a4 岩体ηγK1 0.5 1 0.5 b1 4组以上断裂交汇处 0.5 1 1 b2 2~3组断裂交汇处 1 0.5 0.5 b3 褶皱、派生张裂、环状放射状裂隙 1 0.5 1 c1 元古代片麻状花岗岩γPt1 0.5 0 0.5 c2 钓鱼台组石英岩 1 0 0.5 c3 南芬组泥灰岩 0.5 0.5 0.5 c4 寒武系奥陶系灰岩 0 1 0.5 c5 岩体γδJ2、岩体δK1 1 0 1 h1 1∶20万水系铜异常 1 1 0.5 h2 1∶20万水系金异常 1 1 0.5 h3 1∶20万水系钼异常 0 0.5 0.5 总权重 10 10 10 表 4 三种类型矿床预测区网格权值统计表
Table 4 Statistical table of grid weights of predicted target areas of three types of ore deposits
石英脉型金铜矿床预测区网格加权值 矽卡岩型铜金床矿预测区网格加权值 斑岩型铜钼金矿床预测区网格加权值 x1 x2 x3 x4 x5 x6 x7 x8 x1 x2 x3 x4 x5 x6 x7 x8 x1 x2 x3 x4 x5 x6 x7 x8 y10 2.5 3.5 4.5 6 5.5 3.5 1 1 y10 1.5 2.5 3.5 4 4 3.5 2 1.5 y10 2 3 3.5 4.5 4.5 3 1 1.5 y9 2.5 4.5 4.5 4.5 5.5 3.5 2 2 y9 2 2.5 2.5 3 4.5 3.5 2 3 y9 2.5 3.5 3.5 3.5 4 3 2 2.5 y8 5 5 3.5 7 6 8 3.5 2 y8 3 3 2.5 5 5 6 5 2 y8 4.5 4.5 3 6 5.5 7 4 2 y7 2 3 5.5 8 5.5 7 2.5 1 y7 1.5 2 3.5 6 5.5 5.5 3.5 1.5 y7 2.5 3.5 5 7 5.5 6.5 3 1.5 y6 4 3 5.5 5.5 7.5 6.5 1.5 1 y6 1.5 1.5 4.5 5.5 7.5 5.5 3 2 y6 3.5 3 5.5 5.5 8 6 2 1 y5 2.5 4.5 3.5 4.5 4.5 3.5 1.5 1 y5 1.5 2.5 2.5 4 5 5 3.5 2 y5 2 4 3.5 5 4.5 4.5 2.5 1 y4 0.5 2.5 3 3 2.5 1.5 1.5 0 y4 1 1.5 3 4 4.5 4 3.5 1 y4 0.5 2 3.5 3.5 3 3.5 3 0.5 y3 1.5 3.5 4 1.5 2 1 1 0 y3 1.5 3 2.5 2 2.5 1.5 2.5 1 y3 1 3 3.5 1.5 1.5 1.5 2 0.5 y2 3 5.5 6 3.5 1.5 0 1 1 y2 2.5 5.5 4.5 3.5 2 1 1.5 1.5 y2 2.5 4.5 4.5 3 1.5 0.5 1.5 1.5 y1 3 3.5 5 6 2 0 0 0 y1 3 2.5 4 5.5 2 1 1 1 y1 2.5 2.5 3.5 5 2 0.5 0.5 0.5 -
毕献武, 胡瑞忠, 叶造军, 等. A型花岗岩类与铜成矿关系研究——以马厂箐铜矿为例[J]. 中国科学(D辑), 1999, 29(6): 490-495 BI Xianwu, HU Ruizhong, YE Zaojun, et al. Study on the relationship between A-type granitoids and copper mineralization: A case study of Machangqing Copper Mine[J]. Science in China (Series D), 1999, 29(6): 490-495.
董立军. 辽宁省本溪县陈英-兰河峪铜多金属矿找矿潜力分析[J]. 有色矿冶, 2017, 33(1): 9-12 doi: 10.3969/j.issn.1007-967X.2017.01.003 DONG Lijun. Analysis on the Ore Prospecting Potential of Copper Multi-metal Resources of Chenying-Lanheyu District in Benxi, Liaoning[J]. Non-Ferrous Mining and Metallurgy, 2017, 33(1): 9-12. doi: 10.3969/j.issn.1007-967X.2017.01.003
贺高品, 叶慧文. 辽东-吉南地区早元古代两种类型变质作用及其构造意义[J]. 岩石学报, 1988, 14(2): 152-162 HE Gaopin, YE Huiwen. Two types of Early Proterozoic metamorphism and its tectonic significance in Eastern Liaoning and Southern Jilin Area[J]. Acta PetrologicaSinica, 1988, 14(2): 152-162.
K-51-XXIX(宽甸)幅矿产图说明书1∶20000[R]. 辽宁省地质局区域地质测量队二分队, 1967. 李强, 张捷, 廖勇, 等. 辽东裂谷的基本特征及含矿建造划分[J]. 有色矿冶, 2007, 23(4): 9-12 doi: 10.3969/j.issn.1007-967X.2007.04.004 LI Qiang, ZHANG Jie, LIAO Yong, et al. Basic features of Liaodong rife and categorization of its Ore-bearing formation [J]. Non-Ferrous Mining and Metallurgy, 2007, 23(4): 9-12. doi: 10.3969/j.issn.1007-967X.2007.04.004
李安石. 辽宁省侵入岩与成矿[M]. 北京: 地质出版社, 1993 LI Anshi. Intrusive rocks and mineralization in Liaoning Province[M]. Beijing: Geological Publishing House, 1993.
路远发. GeoKit: 一个用VBA构建的地球化学工具软件包[J]. 地球化学, 2004, 33(5): 460-464 LU Yuanfa. GeoKit: A geochemical toolkit for Microsoft Excel[J]. Geochemistry, 2004, 33(5): 460-464.
刘绍平. 数学地质方法及应用[M]. 北京: 石油工业出版社, 2011 LIU Shaoping. LIU Shaoping.The methods and applications on mathematical geology[M]. Beijing: Petroleum Industry Press, 2011.
刘志宏. 黑龙江省翠宏山钨钼锌多金属矿床地质特征及成因[D]. 长春: 吉林大学, 2009. LIU Zhihong. Geological characteristics and origin of deposit in Cuihongshan W, Mo, Zn Polymetatic Deposit [D]. Changchun: Jilin University, 2009.
刘泽增. 本溪白水金矿成矿区带地质特征及远景[J]. 有色矿冶, 2012, 28(5): 1-3 LIU Zezeng. Geological Characteristics and Mining Prospect of Baishui Gold Mine Belt in Benxi City[J]. Non-ferrous Mining and Metallurgy, 2012, 28(5): 1-3.
梁帅, 翟富荣, 杨占兴, 等. 辽东二棚甸子-万宝地区铜多金属矿床控矿条件分析及成矿预测[J]. 矿床地质, 2013, 32(2): 415-426 doi: 10.3969/j.issn.0258-7106.2013.02.015 LIANG Shuai, ZHAI Furong, YANG Zhanxing, et al. Ore-controlling conditions and metallogenic prognosis of copper polymetallic deposits in Erpengdianzi-Wanbao area of eastern Liaoning Province[J]. Mineral Deposits, 2013, 32(2): 415-426. doi: 10.3969/j.issn.0258-7106.2013.02.015
辽宁地质勘查院. 辽宁1∶5万桥头镇幅、本溪市幅、南芬幅、铁山街幅区域地质调查报告[R]. 2014. 石绍山, 尤洪喜, 成龙, 等. 辽宁本溪地区关门山岩体地球化学特征及构造环境[J]. 地质与资源, 2014, 34(5): 441-445 SHI Shaoshan, YOU Hongxi, CHENG Long, et al. Geochemistry and tectonic environment of the Guanmenshan intrusive rocks in Benxi area, Liaoning Province [J]. Geology and Eesources, 2014, 34(5): 441-445.
沈阳地质调查中心. 辽宁1∶5万本溪县、草河掌、田师傅、南孤山子幅区域地质调查告[R]. 2015. 孙宁, 孔凡乾, 韦龙明, 等. 广东市麻布岗地区银多金属矿床地质特征及找矿方向[J]. 有色金属(矿山部分), 2015, 67(1): 43-46 SUN Ning, KONG Fanqian, WEI Longming, et al. Geological characteristics and exploration potential of Mabugang Ag-Polymetallic Deposits, Guangdong Province[J]. Non-Ferrous Metals (Mining Section), 2015, 67(1): 43-46.
吴迪. 本溪黄柏峪铜钼矿床成矿地质条件浅析[J]. 有色矿冶, 2005, 21(3): 8-11 doi: 10.3969/j.issn.1007-967X.2005.03.003 WU di. The Ore-forming Condition of the Cu-Mo Deposit in Huangbaiyu Area of Benxi[J]. Non-Ferrous Mining and Metallurgy, 2005, 21(3): 8-11. doi: 10.3969/j.issn.1007-967X.2005.03.003
王世称. 综合信息矿产预测理论与方法体系新进展[J]. 地质通报, 2010, 29(10): 1399-1403 doi: 10.3969/j.issn.1671-2552.2010.10.002 WANG Shicheng. The new development of theory and method of synthetic information mineral resources prognosis[J]. Geological Bulletin of China, 2010, 29(10): 1399-1403. doi: 10.3969/j.issn.1671-2552.2010.10.002
叶天竺, 吕志成, 庞振山, 等. 勘查区找矿预测理论与方法(总论)[M]. 北京: 地质出版社, 2014. YE Tianzhu, LV Zhicheng, PANG Zhenshan, et al. Theory and method of prospecting prediction in exploration area (General) [M]. Beijing: Geological Publishing House, 2014.
杨永胜. 大兴安岭中北段与金铜钼矿有关岩浆岩成矿专属性及红彦地区成矿预测[D]. 武汉: 中国地质大学武汉, 2017 YANG Yongsheng. Metallogenic Specialization of Gold, Copper and Molybdenum Mineralized Magmatic Rocks in the Middle-Northern Great Xing'an Range and Metallogenic Prediction in Hongyan Area[D]. Wuhan: China University of Geosciences , 2017.
杨佳佳, 林楠. 基于证据加权模型的综合信息成矿预测[J]. 地质学报, 2016, 90(10): 2908-2918 YANG Jiajia, LIN Nan. Comprehensive Information Prospecting Based on the Weighted Evidence Model[J]. Acta GeologicaSinica, 2016, 90(10): 2908-2918.
朱国宝. 辽宁东部地区中侏罗世黄柏峪岩体地球化学特征及其构造环境[J]. 资源信息与工程, 2019, 34(5): 20-22 doi: 10.3969/j.issn.2095-5391.2019.05.006 ZHU Guobao. Geochemical characteristics and tectonic environment of the Middle Jurassic Huangbaiyu Pluton in eastern Liaoning[J]. Resource Information and Engineering, 2019, 34(5): 20-22. doi: 10.3969/j.issn.2095-5391.2019.05.006
赵鹏大, 孟宪国. 地质异常与矿产预测[J]. 地球科学, 1993, 18(01): 39-47 doi: 10.3321/j.issn:1000-2383.1993.01.001 ZHAO Pengda, MENG Xianguo. Geological anomaly and mineral prediction[J]. Earth Science, 1993, 18(01): 39-47. doi: 10.3321/j.issn:1000-2383.1993.01.001
Blevin P L, Chappell B W. Chemistry, origin and evolution of mineralized granites in the Lachlan Fold Belt, Australia; the metallogeny of I- and S-type granites[J]. Economic Geology, 1995, 90(6): 1604-1619. doi: 10.2113/gsecongeo.90.6.1604
Collins W J, Beams S D, White A J R, et al. Nature and origin of A type granites with particular reference to Southeastern Australia[J]. Contributions to Mineralogy and Petrology, 1982, 80, 189-200. doi: 10.1007/BF00374895
-
期刊类型引用(9)
1. 解铭威,周慧荻,陈耸,王向荣. 银川市生态系统服务价值评估及多情景模拟. 水土保持研究. 2025(01): 294-304 . 百度学术
2. 洪桦,张渝. 基于土地利用变化的生态系统服务价值演变研究——以福州滨海新城为例. 农业与技术. 2024(04): 107-112 . 百度学术
3. 龚芯磊 ,陈鸿申 ,胡世敏 ,李正荣 ,聂坤 ,梁龙 . 特定模型与GIS结合的生态系统服务功能评价及分区研究——以黔北芙蓉江流域绥正道地区为例. 贵州地质. 2024(03): 328-338 . 百度学术
4. 陈彦珍,张丽华,王亚军,宋学云,曾乐,郝爱华,贺甜田. 不同枸杞品种的区域差异性研究. 宁夏农林科技. 2024(07): 46-52 . 百度学术
5. 张永生,李建国,赵广臣. 农田回归用水引入鸣翠湖后水质改善的数值模拟研究. 水电能源科学. 2024(11): 51-54 . 百度学术
6. 朱青青,刘超,沈艳,马红彬,谭松伟,王国会,李燕,李千飞,李国强. 宁夏罗山草地生态系统服务价值的地形效应. 应用生态学报. 2024(12): 3267-3274 . 百度学术
7. 欧阳渊,刘洪,张景华,唐发伟,张腾蛟,黄勇,黄瀚霄,李富,陈敏华,宋雯洁. 西南山区生态地质调查技术方法研究. 西北地质. 2023(04): 218-242 . 本站查看
8. 江山,石绍山,郭常来,冯雨林,孙家全,孙秀波,周丽. 大凌河流域1998~2019年NDVI时空变化及其对气温和降水的响应. 西北地质. 2023(04): 254-262 . 本站查看
9. 高媛,李谦. 湿地生态系统服务价值评估研究现状及展望. 绿色科技. 2023(18): 62-68 . 百度学术
其他类型引用(1)