Intelligent geological mapping technology and its applications based on the multivariate knowledge and mapping model
-
摘要:
中国1∶25万—1∶5万区域地质调查积累了海量地质数据,如何充分利用这些基础地质数据适时更新编制小比例尺(1∶50万—1∶500万)地质图,及时、客观地反映最新地质调查成果,给地质编图工作者提出了严峻挑战。充分利用地质、制图、数据库等多元知识构建地质编图模型,支持计算机对岩石、地层、地质构造等地质实体及其复杂关系的智能识别,有效取代地质编图人员对地质要素综合、属性继承与更新、制图输出等的处理决策,通过图元与属性联动,最终形成图库一体的地质图空间数据库。以编图知识提取、模型构建为重点,对智能地质编图技术进行了探索,并通过青海省东部1∶50万地质图智能编制实践进行分析验证。结果表明,编图知识和构建的模型符合智能地质编图要求,智能地质编图技术能够大大提高从中大比例尺到小比例尺地质图缩编的效率和质量。
Abstract:In China, 1:250 000-1:50 000 regional geological survey has accumulated a large amount of database and formed a rich geological map with corresponding scale.How to make full use of these data and maps to compile small scale geological maps, and systematically reflect the latest geological survey results, brings a severe challenge to geological cartographers.Based on the big data, the geological map models are constructed through the multi-knowledge of geology, cartography and database, which can realize the intelligent recognition of geological phenomena such as stratum, magmatic rock, geological structure and their complex relations by computer, which can effectively replace the decision-making of geological cartographers on the synthesis and plotting of complicated geological map elements.Finally, the geological map spatial database is formed through the linkage of map elements and attributes.Focusing on knowledge extraction and map compilation model construction, this paper explores the technology of intelligent geological map compilation, and analyzes and verifies it through the 1∶500 000 intelligent geological map compilation practice in the eastern part of Qinghai Province.The mapping knowledge and model meet the requirements of intelligent geological mapping.The technology of intelligent geological mapping can greatly improve the efficiency and quality of geological map reduction from medium-large scale to small scale.
-
大兴安岭地区发育广泛花岗岩和火山岩,其岩石成因及构造背景一直是国内外地质学者研究的热点。最新的测年结果显示,这些花岗岩和火山岩主要形成于中生代,少数形成于新元古代和古生代,并对大兴安岭地区岩浆作用的年代学格架和地球化学属性有了一定的认识[1-3],但是由于大兴安岭地区植被覆盖严重,导致部分地区晚中生代岩浆岩的年代学和地球化学数据较缺乏,阻碍了对该区岩浆作用和区域构造演化的深刻认识。
近年来,地质工作者按岩性、岩石组合特征及区域对比,将大兴安岭中部地区晚中生代火山岩划分为南、北两部分,北部主要包括七一牧场组、吉祥峰组、上库力组、伊列克得组,南部主要包括塔木兰沟组、满克头鄂博组、玛尼吐组、白音高老组、梅勒图组,南、北两部分火山岩岩石组合可进行对比,但是对于同一期火山岩不同岩石获得的年龄数据往往具有较大的变化范围,有的甚至相差几十百万年[4-8],表明这种划分方案的合理性仍值得商榷,或许直接用不同岩浆期次来表示不同时代的火山岩比较合理,即用晚侏罗世中基性火山岩、晚侏罗世酸性火山岩、早白垩世中性火山岩、早白垩世酸性火山岩等分别代表不同时代的火山岩。
本次选取大兴安岭中部柴河—蘑菇气地区出露的早白垩世中性火山岩,在野外调查和岩相学研究的基础上,通过岩石学、年代学、岩石地球化学等的系统研究,确定大兴安岭中部柴河—蘑菇气地区早白垩世中性火山岩的岩石成因特点,探讨其反映的区域构造背景。
1. 地质背景及样品描述
大兴安岭位于中亚造山带东段,大地构造位置处于西伯利亚板块和华北板块之间(图 1)[9],其古生代经历了多个微陆块碰撞拼合时期[10-14],构造演化与古亚洲洋的闭合密切相关,中、新生代受太平洋构造域和蒙古-鄂霍茨克构造域的叠加和改造。
图 1 研究区大地构造位置(据参考文献[9]修改)F1—牡丹江断裂;F2—敦化-密山断裂;F3—伊通-佳木斯断裂;F4—西拉木伦-长春-延吉断裂;F5—贺根山-扎兰屯-黑河断裂;F6—塔源-喜桂图断裂;F7—得尔布干断裂Figure 1. Tectonic map of Chaihe-Moguqi area in central Da Hinggan Mountains柴河—蘑菇气地区位于大兴安岭中部,晚中生代侏罗纪—白垩纪以陆相火山岩地层和陆相碎屑沉积岩地层为主,主要包括中侏罗世万宝组、晚侏罗世塔木兰沟组、木瑞组、满克头鄂博组和早白垩世玛尼吐组、白音高老组、梅勒图组、大磨拐河组。其中早白垩世玛尼吐组以中性火山岩及火山碎屑岩为主,主要分布在天池、蛤蟆沟、柴河、新立屯、蘑菇气、中和、山泉等地区(图 2),主要岩性有灰绿色-灰紫色安山岩(图 3)、英安岩、安山质晶屑凝灰熔岩、凝灰角砾熔岩等。
2. 分析方法
锆石用重液和磁选相结合从粉碎的岩石样品中分选出来,再在双目镜下提纯。锆石被嵌于树脂样靶中并抛光,之后渡上50nm厚的纯金。阴极发光(CL)图像在中国地质科学院北京离子探针中心制成,主要是查明锆石内部结构,以便准确选点。锆石U-Pb定年在天津地质矿产研究所完成。测试仪器为NEPTUNE多接收器电感耦合等离子体质谱仪(MC-ICP-MS)和氟化氩准分子激光剥蚀系统(NEW WAVE 193nm FX)联机[15]。激光剥蚀斑束直径为35μm,激光剥蚀样品的深度为20~40μm。实验中采用氦作为剥蚀物质的载气。锆石年龄计算采用源自澳大利亚Macquarie大学的标准锆石GJ-1作为外标,元素含量采用人工合成玻璃NIST SRM610作为外标[16]。具体实验操作流程见参考文献[17]。样品的同位素比值和元素含量处理采用Glitter(ver4.0)程序计算,普通Pb采用Anderson的3D坐标法校正[18],年龄加权平均值及谐和图采用Isoplot(v. 3.00)[19]绘制。单个测试点的同位素比值和同位素年龄误差为1σ,206Pb/238U年龄加权平均值置信度为95%。
样品的主量和微量元素分析在天津地质矿产研究所和自然资源部哈尔滨矿产资源监督测试中心完成。挑选新鲜未风化样品洗净、烘干后,用玛瑙研钵研磨至200目的粉末,用于主量和微量元素分析。主量元素采用X射线荧光光谱分析(XRF),微量元素的分析采用等离子质谱法(ICP-MS)完成,详细的主量、微量元素分析实验原理和分析步骤见参考文献[20]。
在LA-ICP-MS锆石U-Pb定年的基础上,进行了锆石原位Hf同位素测定工作。锆石原位Hf同位素分析在天津地质矿产研究所完成,测试仪器为NEPTUNE多接收器电感耦合等离子体质谱仪和氟化氩准分子激光剥蚀系统联机,激光束斑直径为50μm。详细的分析方法见参考文献[15]。分析过程中标准锆石GJ-1的176Hf/177Hf值为0.282000 ± 0.000040。用实验测得的176Hf/177Hf和176Lu/177Hf值计算初始176Hf/177Hf值。176Lu的半衰期采用1.865× 10-11/a[21]。计算εHf(t)值时采用的现在的(176Hf/177Hf)DM和(176Lu/177Hf)DM值分别为0.28325和0.0384[22]。Hf模式年龄的计算见参考文献[22-24]。
3. 分析结果
3.1 锆石U-Pb定年
样品Z11-84岩性为安山质角砾凝灰熔岩,共24个锆石分析点,分析结果见表 1。锆石多为短柱状和粒状,少数为长柱状,粒度为50~200μm,长宽比为1:1~3:1。CL图像显示(图 4-a),这些锆石内部具有典型的环带结构,Th/U值为0.65~1.10,指示这些锆石属于岩浆锆石。在锆石U-Pb谐和图(图 4-b)中,样品Z11-84的锆石U-Pb年龄测定数据点全部位于谐和线上及其附近,206Pb/238U年龄加权平均值为139±1Ma(MSWD=0.59),时代为早白垩世,代表了安山质角砾凝灰熔岩的形成年龄。
表 1 研究区早白垩世中性火山岩测年结果Table 1. Data of zircons for the Early Cretaceous intermediate volcanic rocks in study area样品号 岩性 年龄/Ma 采样位置 方法 数据来源 Z11-84 安山质角砾凝灰熔岩 139±1 碾子山北 LA-ICP-MS A P010-5-2 英安岩 141±5 苏河屯 LA-ICP-MS B P9B2-2 安山岩 137±2 青年林场北山 SHRIMP C P11B11-3 角闪玄武安山岩 131±3 河中林场 SHRIMP C TW1 安山岩 133±7 哈布气林场 Rb-Sr D DH2011RZ27 安山岩 141±2 新立屯 LA-ICP-MS E 注:A—本文数据;B—吉林省区域地质矿产调查所提供数据;C—吉林大学刘正宏教授提供数据;D—辽宁省地质矿产调查院提供数据;E—中国地质调查局沈阳地质调查中心提供数据;Z11-84采样点坐标为北纬47°39′51″、东经122°51′34″ 结合笔者收集的年龄数据(表 2)可知,研究区早白垩世中性火山岩的形成年龄介于140~130Ma之间。
表 2 早白垩世中性火山岩测年数据Table 2. Age statistics for the Early Cretaceous intermediate volcanic rocks分析号 Th/10-6 U/10-6 Th/U 同位素比值 年龄/Ma 206Pb/238U 1σ 207Pb/235U 1σ 207Pb/206Pb 1σ 206Pb/238U 1σ 207Pb/235U 1σ 207Pb/206Pb 1σ Z11-84.1 26 33 0.80 0.021 0.0005 0.1767 0.0176 0.0608 0.0053 134 3 165 16 875 447 Z11-84.2 46 42 1.10 0.0222 0.0004 0.1950 0.0161 0.0639 0.0047 141 2 181 15 1120 343 Z11-84.4 27 39 0.69 0.0213 0.0006 0.1506 0.0182 0.0513 0.0033 136 4 142 17 146 598 Z11-84.5 28 33 0.86 0.0217 0.0005 0.1719 0.0189 0.0574 0.0049 138 3 161 18 852 497 Z11-84.6 20 26 0.78 0.0211 0.0006 0.1571 0.0197 0.0539 0.0069 135 4 148 19 181 625 Z11-84.7 41 42 0.98 0.0216 0.0004 0.2040 0.0114 0.0685 0.0048 138 3 189 11 1662 207 Z11-84.8 40 47 0.86 0.0219 0.0003 0.1742 0.0126 0.0576 0.0045 140 2 163 12 514 316 Z11-84.9 33 43 0.78 0.0218 0.0003 0.1635 0.0125 0.0544 0.0042 139 2 154 12 755 322 Z11-84.10 23 31 0.75 0.0219 0.0004 0.1714 0.0103 0.0567 0.0037 140 3 161 10 1997 220 Z11-84.11 27 39 0.70 0.0218 0.0004 0.1474 0.0158 0.0490 0.0049 139 3 140 15 147 587 Z11-84.12 27 32 0.83 0.0216 0.0005 0.2013 0.0115 0.0677 0.0058 138 3 186 11 2164 196 Z11-84.13 40 50 0.79 0.0220 0.0003 0.1660 0.0099 0.0547 0.0030 140 2 156 9 1222 233 Z11-84.14 27 36 0.75 0.0222 0.0004 0.1502 0.0176 0.0491 0.0039 142 2 142 17 152 565 Z11-84.15 25 33 0.76 0.0217 0.0005 0.1671 0.0096 0.0558 0.0055 138 3 157 9 1787 208 Z11-84.16 43 55 0.79 0.0215 0.0005 0.1815 0.0101 0.0611 0.0039 137 3 169 9 2076 194 Z11-84.17 29 41 0.72 0.0219 0.0003 0.1750 0.0125 0.0579 0.0045 140 2 164 12 1447 281 Z11-84.18 48 54 0.87 0.0220 0.0003 0.1899 0.0096 0.0625 0.0042 140 2 177 9 1329 199 Z11-84.19 24 33 0.72 0.0219 0.0005 0.2673 0.0231 0.0887 0.0097 139 3 241 21 1398 365 Z11-84.20 30 39 0.76 0.0215 0.0004 0.1525 0.0147 0.0514 0.0063 137 2 144 14 152 887 Z11-84.21 28 36 0.77 0.0201 0.0008 0.1849 0.0132 0.0668 0.0063 128 5 172 12 1677 271 Z11-84.22 29 40 0.71 0.0220 0.0007 0.1902 0.0174 0.0627 0.0074 140 4 177 16 697 377 Z11-84.23 27 37 0.74 0.0220 0.0004 0.1838 0.0157 0.0607 0.0049 140 3 171 15 1502 339 Z11-84.24 35 53 0.65 0.0220 0.0003 0.1636 0.0087 0.0539 0.0036 140 2 154 8 1034 216 Z11-84.25 42 55 0.77 0.0217 0.0002 0.1806 0.0080 0.0602 0.0025 139 2 169 7 1203 177 3.2 主量和微量元素
研究区早白垩世中性火山岩的主量和微量元素分析结果见表 3。早白垩世中性火山岩SiO2含量在55.1%~69.1%之间,TiO2含量介于0.39%~1.68%之间,Al2O3含量较高,在14.8%~18%之间,TFeO含量介于3.62% ~7.83%之间,MgO含量低,在0.45% ~ 3.74%之间,CaO含量较高,介于1.16%~6.57%之间,Na2O含量在1.73%~4.72%之间,K2O含量在1.66%~ 4.84%之间,整体为一套富钾钠、富铝,贫镁的岩石。在TAS图解(图 5)中,主要落在安山岩、粗面岩(标准矿物石英含量小于20%)内,属亚碱性系列,含铝指数A/CNK为0.80~1.13,属准铝质到弱过铝质岩石;在岩浆系列判别图解(图 6)中,主要落在高钾钙碱性系列。
表 3 早白垩世中性火山岩主量、微量和稀土元素分析结果Table 3. Major, trace and rare earth elements analyses of the Early Cretaceous intermediate volcanic rocks样品号 B1068-1 B2005-2 B2010-1 B2022-1 B2034-1 B2082-2-1 B3007-1 B4063-1 B5049-3 HfP38-2-1 SiO2 58.40 55.40 62.30 68.40 60.90 63.80 69.10 67.00 55.10 60.42 T1O2 0.98 1.68 0.44 0.39 0.76 0.66 0.55 0.58 1.09 0.91 Al2O3 16.90 15.70 15.30 15.70 16.35 16.60 14.80 15.10 16.90 18.00 Fe2O3 3.62 6.65 2.36 1.78 2.29 3.41 2.47 1.99 3.46 2.90 FeO 3.52 1.85 1.20 1.02 2.95 1.12 0.62 2.02 4.70 1.52 MnO 0.10 0.11 0.07 0.04 0.10 0.07 0.06 0.08 0.13 0.11 MgO 2.67 3.15 1.67 0.45 1.96 1.40 0.45 1.38 3.74 1.74 CaO 5.28 6.57 4.06 2.08 5.06 2.74 1.16 2.29 5.78 3.81 Na2O 4.27 3.25 1.73 4.56 3.05 4.50 3.99 3.75 3.37 4.72 K2O 1.66 2.12 3.01 3.80 3.01 3.58 4.84 4.43 2.74 3.36 P2O5 0.27 0.91 0.12 0.08 0.24 0.23 0.12 0.17 0.37 0.37 烧失量 1.78 2.14 7.36 1.09 2.67 1.47 1.29 0.64 1.82 1.94 总计 99.45 99.53 99.62 99.40 99.34 99.58 99.45 99.43 99.20 99.80 σ 2.28 2.33 1.16 2.75 2.05 3.14 2.99 2.79 3.09 3.75 A/CNK 0.92 0.80 1.13 1.02 0.93 1.02 1.06 1.00 0.89 0.98 Na2O+K2O 5.93 5.37 4.74 8.36 6.06 8.08 8.83 8.18 6.11 8.08 Cr 30.20 70.60 17.80 15.90 20.60 16.50 14.20 24.50 30.00 3.21 Rb 24.8 32.4 100 94.9 64.1 86.1 115 76.4 39.4 101.4 Sr 563 807 1480 388 668 571 491 400 671 664.4 Ba 659 953 1190 1020 1020 979 1180 782 664 1065 Th 4.14 3.48 6.02 6.75 4.84 7.32 6.58 8.26 3.6 9.69 U 1.04 0.98 1.39 1.49 1.19 1.95 1.91 2.06 1.18 2.49 Nb 6.76 17.2 6.81 7.34 6.34 9.18 9.25 9.42 6.78 7.69 Ta 0.5 1.1 0.4 0.4 0.41 0.54 0.47 0.64 0.36 0.53 Zr 224 405 235 290 190 247 271 210 158 194.8 Hf 5.89 9.71 6.38 7.91 5.04 6.6 7.24 6.04 4.38 14.91 La 24.2 60.4 26.2 26.3 24.1 32.6 30.6 27.8 19.3 32.62 Ce 50 131 51.7 53.1 48.3 67.6 60.5 55.9 40 62.37 Pr 6.99 16.8 6.6 6.1 5.96 8.1 8 7.21 5.98 7.19 Nd 29.4 68.3 24.5 23.3 23 31.1 32.4 27.3 25.7 29.78 Sm 6.01 12.2 4.44 4.11 4.22 5.51 6.09 4.98 5.4 5.75 Eu 1.59 2.77 1.22 1.21 1.32 1.48 1.56 1.1 1.52 1.96 Gd 4.84 9.27 3.63 3.32 3.62 4.43 4.79 4.08 4.43 4.15 Tb 0.71 1.44 0.56 0.55 0.56 0.67 0.8 0.69 0.78 0.73 Dy 4.15 6.74 3.04 3.03 2.9 3.72 4.76 3.79 4.13 3.67 Ho 0.82 1.19 0.61 0.61 0.56 0.72 0.93 0.72 0.78 0.73 Er 2.31 3.24 1.74 1.78 1.64 2.01 2.66 2 2.28 1.92 Tm 0.35 0.44 0.27 0.28 0.25 0.32 0.4 0.31 0.33 0.3 Yb 2.28 2.76 1.89 1.97 1.62 2.15 2.83 2.06 2.13 2.2 Lu 0.35 0.42 0.31 0.31 0.26 0.34 0.44 0.32 0.32 0.27 Y 22.2 31 16.8 16.6 15.4 20.2 24.9 18.8 21.6 20.3 ∑REE 156.20 347.97 143.51 142.57 133.71 180.95 181.66 157.06 134.68 173.93 LR/HR 7.48 11.43 9.52 9.63 9.37 10.19 7.90 8.90 6.45 9.99 δEu 0.87 0.77 0.90 0.97 1.01 0.89 0.85 0.72 0.92 1.17 注:主量元素含量单位为%,微量和稀土元素含量为10-6 早白垩世中性火山岩微量元素原始地幔标准化蛛网图(图 7-a)显示,样品主要富集大离子亲石元素(如Rb、K、LREE)和Ba、Sr元素,亏损高场强元素(如Ta、Nb、P、Ti),是壳源岩浆或岩浆被地壳物质混染的典型特征。
稀土元素总量(ΣREE)在133.71 × 10-6~ 347.97×10-6之间,平均值为175.22×10-6,轻、重稀土元素比值LREE/HREE在6.45~11.43之间,平均值为9.09,轻、重稀土元素分馏明显,分馏系数(La/Yb)N介于6.50~15.70之间,平均值为9.89。球粒陨石标准化蛛网图(图 7-b)中,均呈右倾配分曲线模式,轻稀土元素富集、重稀土元素亏损,无明显Eu异常,δEu值介于0.72~1.17之间,平均值为0.91。
3.3 锆石原位Hf同位素
在LA-ICP-MS锆石U-Pb定年的基础上,对早白垩世中性火山岩(Z11-84)进行了锆石原位Hf同位素分析,分析结果列于表 4。
表 4 早白垩世中性火山岩锆石原位Hf同位素分析结果Table 4. In situ zircon Hf isotopic compositions of the Early Cretaceous intermediate volcanic rocks分析点 年龄/Ma 176Yb/l77Hf 2σ 176Lu/177Hf 2σ 176Hf/177Hf 2σ 176Hf/177Hf (corr) 2σ εHf(0) εHf(t) 2σ TDM1(Hf)/Ma TDM2(Hf)/Ma fLu/Hf Z11-84.1 134 0.0300 0.0000 0.0008 0.0000 0.282977 0.000019 0.282981 0.000019 7.4 10.26 0.7 382 482 -0.98 Z11-84.2 141 0.0430 0.0001 0.0011 0.0000 0.282859 0.000025 0.282863 0.000025 3.2 6.21 0.9 553 714 -0.97 Z11-84.4 136 0.0446 0.0001 0.0012 0.0000 0.282941 0.000023 0.282944 0.000023 6.1 8.96 0.8 439 556 -0.96 Z11-84.5 138 0.0347 0.0003 0.0010 0.0000 0.282971 0.000024 0.282975 0.000024 7.2 10.12 0.9 393 493 -0.97 Z11-84.6 135 0.0359 0.0001 0.0011 0.0000 0.282975 0.000019 0.282979 0.000019 7.3 10.18 0.7 388 487 -0.97 Z11-84.7 138 0.0458 0.0001 0.0015 0.0000 0.282979 0.000025 0.282983 0.000025 7.5 10.35 0.9 386 479 -0.96 Z11-84.8 140 0.0249 0.0001 0.0007 0.0000 0.282949 0.000023 0.282953 0.000023 6.4 9.40 0.8 421 534 -0.98 Z11-84.9 139 0.0204 0.0000 0.0006 0.0000 0.282965 0.000020 0.282968 0.000020 6.9 9.93 0.7 398 504 -0.98 Z11-84.10 140 0.0224 0.0000 0.0007 0.0000 0.282954 0.000020 0.282958 0.000020 6.6 9.57 0.7 414 525 -0.98 Z11-84.11 139 0.0447 0.0002 0.0012 0.0000 0.282904 0.000023 0.282907 0.000023 4.8 7.73 0.8 492 628 -0.96 Z11-84.12 138 0.0301 0.0001 0.0008 0.0000 0.282932 0.000025 0.282936 0.000025 5.8 8.73 0.9 446 570 -0.98 Z11-84.14 142 0.0297 0.0002 0.0008 0.0000 0.282883 0.000025 0.282886 0.000025 4.0 7.08 0.9 515 666 -0.98 Z11-84.15 138 0.0318 0.0001 0.0008 0.0000 0.282917 0.000032 0.282921 0.000032 5.3 8.23 1.1 467 599 -0.97 Z11-84.16 137 0.0301 0.0001 0.0008 0.0000 0.282994 0.000025 0.282997 0.000025 8.0 10.91 0.9 359 448 -0.98 Z11-84.17 140 0.0411 0.0001 0.0011 0.0000 0.282963 0.000023 0.282967 0.000023 6.9 9.84 0.8 406 510 -0.97 Z11-84.18 140 0.0367 0.0001 0.0010 0.0000 0.282894 0.000030 0.282897 0.000030 4.4 7.42 1.1 503 646 -0.97 Z11-84.19 139 0.0351 0.0001 0.0011 0.0000 0.282991 0.000028 0.282995 0.000028 7.9 10.84 1.0 365 453 -0.97 Z11-84.20 137 0.0341 0.0000 0.0010 0.0000 0.282976 0.000024 0.282980 0.000024 7.4 10.28 0.8 385 483 -0.97 Z11-84.22 140 0.0354 0.0001 0.0010 0.0000 0.282966 0.000028 0.282970 0.000028 7.0 10.00 1.0 399 501 -0.97 碾子山北部早白垩世安山质角砾凝灰岩(Z11-84)中19颗锆石的初始176Hf/177Hf值为0.282863~ 0.282997,εHf(t)值介于+6.21~+10.91之间,其加权平均值为+8.96±0.71,单阶段模式年龄(TDM1)和二阶段模式年龄(TDM2)变化范围分别为359~553Ma和448~714Ma。
由上可见,研究区早白垩世中性火山岩的锆石εHf(t)值均为正值,且具有年轻的二阶段模式年龄(TDM2),二阶段Hf模式年龄集中在新元古代和古生代。在t-εHf(t)图解(图 8)上[27],数据点全部落在球粒陨石演化线和亏损地幔演化线之间,同时落入兴蒙造山带东段范围内,说明早白垩世中性火山岩为地幔来源或来自新增生的年轻地壳物质,与中亚造山带内显生宙火成岩的Hf同位素组成相似[28-29]。
4. 讨论
4.1 早白垩世中性火山岩的岩石成因
研究区早白垩世中性火山岩化学成分显示主要为安山岩、粗面岩等,主量元素总体呈现富钾钠、富铝、贫镁的特征,为准铝质-弱过铝质、高钾钙碱性岩石。研究表明,粗面质火山岩的成因主要包括玄武质岩浆的分异作用[30-32]、玄武质岩浆与流纹质岩浆的混合作用[33-34]和镁铁质地壳熔融[4, 35] 3种观点。研究区早白垩世中性火山岩中铁族元素Cr含量为3.21×10-6~30.20×10-6,远低于地幔橄榄岩源区部分熔融形成的原始玄武质岩浆(Cr=500×10-6~ 600×10-6),且在研究区缺乏同时期的基性岩,因此玄武质岩浆的分异作用很难解释本区大范围分布的中性火山岩。在La-La/Sm图解(图 9)中,早白垩世中性火山岩的La/Sm值与La呈正相关性,表明其为部分熔融形成。
早白垩世中性火山岩富集大离子亲石元素(如Rb、K、LREE)和Ba、Sr元素,亏损高场强元素(如Ta、Nb、P、Ti);富集轻稀土元素(LREE)、亏损重稀土元素(HREE);Ce/Pb值主要介于2.87~4.10之间,平均值为3.36,与地壳平均值3.3相近[36];La/Nb值为2.85~4.24(平均值为3.52),明显高于原始地幔(0.94),接近陆壳的比值(2.2)。这些特征指示,早白垩世中性火山岩应该来源于下地壳岩石的部分熔融。研究区早白垩世安山质角砾凝灰岩(Z11-84)的锆石εHf(t)值均为正值(+6.21~+10.91),且具有年轻的二阶段模式年龄(TDM2=448~714Ma),二阶段Hf模式年龄集中在新元古代和古生代,与大兴安岭中部晚中生代安山岩的Sr-Nd同位素特征(具有低的Sr初始比值和正的εNd(t)值)一致[37],表明早白垩世中性火山岩的岩浆源区主要为地幔或新元古代和古生代期间从地幔新增生的地壳物质。Nb/Ta值集中在13.52~19.68(平均值为16.47),接近球粒陨石和原始地幔的平均值(17.5)[38-39],而高于地壳的Nb/Ta平均值(12~13)[40];Zr/Hf值集中在34.77~ 41.71(平均值为37.40),高于地壳的Zr/Hf平均值(33)[41],与地幔的平均值(37)接近[42],也反映了岩浆源区为从亏损地幔新增生的下地壳物质,所以保留了地幔的部分微量元素的特点。
研究区早白垩世中性火山岩的SiO2含量为55.1%~69.1%,Al2O3为14.8%~18%,MgO为0.45%~ 3.74%,Y为15.4×10-6~31×10-6,Yb为1.62×10-6~ 2.83×10-6,Sr为388×10-6~1480×10-6,无负Eu异常,这些地球化学特征与埃达克岩接近[43-44]。在Yb-Sr/Yb图解(图 10)[45]中,早白垩世中性火山岩部分样品也落入埃达克岩的区域内,结合其K2O的含量为1.66%~4.84%,平均值为3.26%,应为C型埃达克岩,而C型埃达克岩通常认为是玄武质岩浆底侵加厚的地壳底部引起下地壳部分熔融形成的[46]。早白垩世中性火山岩的Nb含量为6.34×10-6~17.20×10-6(平均值为8.68×10-6),Y含量为15.4×10-6~31×10-6(平均值为20.78×10-6),Zr/Y值为7.31~17.47(平均值为11.81),Rb/Sr值为0.04~0.24(平均值为0.13),Ba/Rb值主要集中在10.24~16.85(平均值为12.22),Ba/Sr值为0.80~2.63(平均值为1.60),这些微量元素特征均反映了研究区早白垩世中性火山岩与安第斯地区安山岩和大陆岛弧安山岩相似(表 5),可能与早先俯冲事件带来的流体影响有关[37]。综上所述,研究区早白垩世中性火山岩岩浆源区应为新元古代和古生代期间新增生的下地壳玄武质岩石的部分熔融,并可能受到早先俯冲事件带来的流体的影响。
表 5 早白垩世中性火山岩微量元素特征与不同构造环境安山岩对比[37]Table 5. Trace element characteristics of the Early Cretaceous intermediate volcanic rocks and a comparison with andesites安山岩类 Nb Y Zr/Y Rb/Sr Ba/Rb Ba/Sr 研究区中性火山岩平均值 8.68 20.78 11.81 0.13 12.22 1.60 安第斯安山岩平均值 10.0 15.0 14.6 0.11 10.0 1.05 大陆岛弧安山岩平均值 9.4 22.0 5.4 0.09 8.2 1.16 大洋岛弧低钾安山岩平均值 0.8 25.0 2.2 0.04 18.5 0.61 大洋岛弧其他安山岩平均值 5.0 20.0 4.7 0.08 13.0 0.73 注:Nb和Y的单位为10-6 4.2 早白垩世中性火山岩形成的构造背景
从侏罗纪开始,东北地区的构造格架与其前期历史相比发生了深刻的变化,由东西向构造转变为北东向和北北东向构造,进入环太平洋构造体系和蒙古-鄂霍茨克构造体系的演化阶段[47-48]。对于东北地区大面积晚中生代岩浆岩形成时的构造背景探讨也一直存在争议,主要观点有以下4种:①与地幔柱构造有关[49-51];②与蒙古鄂霍茨克洋的闭合和造山后垮塌有关[52-53];③与古太平洋板块的俯冲作用有关[7, 54];④与蒙古-鄂霍茨克洋闭合和古太平洋俯冲的共同作用有关[5, 55-56]。
研究发现,大兴安岭及其邻区并不存在环状火山岩带,同时中生代火山岩的形成时代范围较大,因此与喷发周期较短的地幔柱构造有关的作用模式很难解释中生代火山岩的时空分布问题[52, 57]。吉黑东部早—中侏罗世(173~190Ma)钙碱性火山岩组合的发现[58-59]和小兴安岭—张广才岭地区同时期双峰式火成岩组合的存在[55, 60-61],表明早侏罗世古太平洋板块向亚洲大陆斜向俯冲,吉黑东部早中侏罗世岩浆作用应该与本次俯冲事件有关。同时,对于吉黑东部火山岩的年代学研究显示,该区尚未发现164~140Ma的火山岩年龄,而145~130Ma的岩浆活动主要分布于松辽盆地以西的大兴安岭及满洲里地区[1, 4, 62-66],表明古太平洋板块的俯冲在中侏罗世晚期—早白垩世早期存在一个间歇期[47],暗示大兴安岭地区中侏罗世—早白垩世早期的岩浆作用与古太平洋板块的俯冲作用无关。早期由于蒙古-鄂霍茨克洋构造演化的研究程度较低,并且通常认为仅存在西伯利亚南缘的北向俯冲[67],而东北地区北缘则处于被动大陆边缘环境,因此很难将大兴安岭地区晚中生代的岩浆活动与蒙古-鄂霍茨克洋的演化联系到一起。近年来,随着蒙古-鄂霍茨克洋缝合带研究程度的提高,尤其是蒙古-鄂霍茨克洋存在向南俯冲的确定[68-70],越来越多的学者认为,大兴安岭地区晚中生代的岩浆活动应与蒙古-鄂霍茨克洋的闭合和造山后垮塌有关[53, 57]。
在(Y+Nb)-Rb构造环境判别图解(图 11)中[71],早白垩世中性火山岩落入后碰撞构造环境。早白垩世早期中性火山岩中部分样品具有C型埃达克岩的特点,表明研究区经历了陆壳加厚的过程。综上所述,大兴安岭中部柴河—蘑菇气地区早白垩世中性火山岩的形成与蒙古-鄂霍茨克洋闭合陆壳加厚之后的岩石圈伸展环境有关。
5. 结论
通过对大兴安岭中部柴河—蘑菇气地区中性火山岩的岩相学、锆石U-Pb年代学、岩石地球化学和锆石原位Hf同位素研究,可以得出如下结论:大兴安岭中部柴河—蘑菇气地区早白垩世中性火山岩的形成年龄介于140~130Ma之间,岩浆源区应为新元古代和古生代期间新增生的下地壳玄武质岩石的部分熔融,并可能受到早先的俯冲事件带来的流体的影响,形成于蒙古-鄂霍茨克洋闭合陆壳加厚之后的岩石圈伸展环境。
致谢: 感谢审稿专家及吉林大学地球探测科学与技术学院王永志教授对本文提出的修改意见。 -
表 1 地质体面要素数据模型
Table 1 Data model of geological polygon
第0列 第1列 第2列 第3列 序号 编图原始资料地质体 成果图地质体 备注 数据项名称 数据项编码 数据项名称 数据项编码 1 地质体标识号 Feature_Id 地质体标识号 Feature_Id 重置Id 2 地质体类型代码 Feature_Type 地质体代号 Feature_Type 依模型转换 3 地质体名称 Geobody_Name 地质体名称 Geobody_Name 依模型转换 4 地质体时代 Geobody_Era 地质体时代 Geobody_Era 属性继承 5 地质体下限年龄值 Geobody_Age1 地质体下限年龄值 Geobody_Age1 属性继承 6 地质体上限年龄值 Geobody_Age2 地质体上限年龄值 Geobody_Age2 属性继承 7 子类型标识 Subtype 子类型标识 Subtype 属性继承 8 地质体主要岩石名称 Rock_Name 属性提取 9 地质体岩石特征 Rock_Character 属性提取 10 地层分区 Stratum_Division 属性提取 11 地质体重要性标识 Sign 重新赋值 12 地质体更新标识 updata 重新赋值 表 2 地质图图示模型分类
Table 2 Classification models by graphic representation of geological map
主要类别 描述 地图定义 图层组成列表、地图窗口位置等 图层定义 图层命名、图层间压盖关系等图层属性设置 页面定义 定义图页面内容、结构,如图名、比例尺、图例、图签、镶图、页面尺寸等地图要素 图形要素 定义地图编图单位,每类图元内容的属性、符号、尺寸、纹理、颜色、方位、形状等内容 符号化规则 定义图内要素的符号、包括符号设定的参数、放置位置、对齐方式、注记要求等 符号库 预定义的各种可用于符号化要素的符号集(基于制图系统平台制作) 颜色方案 图中使用的所有颜色方案、配色要求 表 3 地质体面要素图示模型
Table 3 Graphical model of geological polygon
数据项名称 数据项编码 数据项类型 数据长度 填充模式 F_type Integer 2 图案编号 F_P_number Integer 4 图案高 F_P_height Float 5.2 图案宽 F_P_width Float 5.2 图案颜色 F_P_color Integer 4 图案角度 F_P_angle Float 5.2 图案笔宽 F_P_penwidth Float 5.2 填充颜色 F_color Integer 4 地质体代码 Feature_Type Character 30 表 4 地质图关系模型示例
Table 4 Example of relation model for geological map
Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ 地质年代 成果地质图编图单位 原始地质图编图单位 地质体类型 重要地质体 代 纪 世 地质体代号 地质体名称 地质体代号 地质体名称 子类型标识 重要性标识 新
↓
↓
↓
老Qhal 全新世冲积 Qh$a$l 全新统冲积 0 0 … … … … … … 角度不整合 K1H 河口群 K@1H$3 河口群三组 0 0 K@1H$2 河口群二组 0 0 K@1H$1 河口群一组 0 0 角度不整合 δμJ 侏罗纪闪长玢岩 Jδμ 灰绿色闪长玢岩 1 0 … … … … … … 注:$表示其后第一个字母或数字为上标;@表示其后第一个字母或数字为下标;以下表中同 表 5 地质体接触关系与判别规则
Table 5 Contact relation of geological bodies and the distinguishing rules
序号 类型 判别条件 1 整合接触 1.SubtypeA=0或3或5;SubtypeB=0或3或5
2.A和B两地层单位相邻且其间无“不整合”相隔2 平行不整合接触 1.SubtypeA=0或3或5;SubtypeB=0或3或5
2.A和B两地层单位间有“平行不整合”注释3 角度不整合接触 1.SubtypeA=0或3或5;SubtypeB=0或3或5
2.A和B两地层单位间有“平行不整合”注释4 侵入接触 1.SubtypeA=1或2或3或6
2.A位于B上部(A晚于B形成时代)5 沉积接触 1.SubtypeA=1或2或3或6
2.A位于B下部(A早于B形成时代)6 断层接触 1.A与B地质体间有断层通过;或A与B公共边与断层线重叠的地质界线
2.SubtypeA=47 构造接触 1.A位于B上部(A晚于B形成时代)。“混杂岩”、“岩组”(特殊地质体4)等晚于与其接触地质体 8 水体接触 1.SubtypeB=9 9 其他接触 1.SubtypeB=7或8 10 第四纪松散堆
积物覆盖界线1.SubtypeA=0
2.Feature_TypeA like “q%”(第四纪地质体)注:(1)A—地质体,B—与A相邻地质体;SubtypeA—地质体A的地质体类型;SubtypeB—地质体B的地质体类型;Feature_TypeA——地质体A的地质体代号;(2)地质体类型-Subtype取值意义:0—沉积(火山)岩岩石地层单位;1—侵入岩岩石年代单位;2—侵入岩谱系单位;3—变质岩地(岩)层单位;4—特殊地质体;5—非正式地层单位;6—脉岩;7—戈壁沙漠;8—冰川与终年积雪;9—水体 表 6 临夏地区西秦岭地层分区1∶25万和1∶50万地质图关系模型(部分)
Table 6 1:250 000 and 1:500 000 geological map relation for stratum subarea of West Qinling in Linxia area
Ⅰ Ⅱ Ⅲ Ⅳ Ⅴ 地质年代 1:50万地质图编图单位 1:25万地质图编图单位 地质体类型 重要地质体 代 纪 世 地质体代号 地质体名称 地质体代号 地质体名称 子类型标识 重要性标识 新生代 第四纪 全新世 Qhal 全新世冲积 Qh$a$l 全新统冲积 0 0 Qhalp 全新世冲洪积 Qh$p$a$l 全新统冲洪积 0 0 Qh$p$a$l(T@Ⅰ) 洪冲积(一级阶地) 0 0 Qh$p$a$l(T@Ⅱ) 洪冲积(二级阶地) 0 0 晚更新世 Qp3m 晚更新世马兰组 Qp@3m$2 马兰组上段 0 0 Qp@3m$1 马兰组下段 0 0 Qp3alp 晚更新世冲洪积 Qp@3$p$a$l 上更新统冲洪积 0 0 Qp$3$p$a$l(T@Ⅲ) 洪冲积(三级阶地) 0 0 Qp$3$p$a$l(T@Ⅳ) 洪冲积(四级阶地) 0 0 中更新世 Qp2alp 中更新世冲洪积 Qp@2$p$a$l 中更新统洪积 0 0 Qp$3$p$a$l(T@Ⅳ) 洪冲积(四级阶地) 0 0 早更新世 Qp1alp 早更新世冲洪积 Qp@1$p$a$l 下更新统冲洪积 6 1 Qp1d 早更新世东山组 Qp@1d 东山组 0 0 Q 第四纪成因未分 Q 第四系未分 0 0 q 石英脉 q 石英脉 6 1 γ 花岗岩脉 γ 花岗岩脉 6 0 Σ 超基性岩 Σ 超基性岩 6 1 新近纪 NG 新近纪甘肃群 N@2j 积石组 0 0 N@2h 何王家组 0 0 N@1m 毛沟组 0 0 E@3t 他拉组 0 0 古近纪 E1-2x 古—始新世西柳沟组 Ex 西柳沟组 0 0 EX 古近纪西宁群 Ey 野狐城组 0 0 角度不整合 0 中生代 白垩纪 早白垩世 K1H 河口群 K@1H$3 河口群三组 0 0 K@1H$2 河口群二组 0 0 K@1H$1 河口群一组 0 0 角度不整合 侏罗纪 δμJ 侏罗纪闪长玢岩 Jδμ 灰绿色闪长玢岩 1 0 γδJ 侏罗纪花岗闪长岩 Jγδ 灰白色花岗闪长岩 1 0 ηγJ 侏罗纪二长花岗岩 Jηγ 浅肉红色二长花岗岩 1 0 δοJ 侏罗纪石英闪长岩 Jδο 灰白色石英闪长岩 1 0 … … … … … … … … -
乔秀夫, 丁孝忠. 区域地质调查与编图综合研究——写在新编《中国地质图集》出版之际[J]. 地质通报, 2003, 22(10): 769-774. doi: 10.3969/j.issn.1671-2552.2003.10.004 王家耀. 空间数据自动综合研究进展及趋势分析[J]. 测绘科学技术学报, 2008, 25(1): 1-8. https://www.cnki.com.cn/Article/CJFDTOTAL-JFJC200801002.htm 武芳, 巩现勇, 杜佳威. 地图制图综合回顾与前望[J]. 测绘学报, 2017, 46(10): 1645-1664. doi: 10.11947/j.AGCS.2017.20170287 李仰春, 王永志, 陈圆圆, 等. 智绘地质——新一代智能化地质编图模式及应用[J]. 地质通报, 2020, 39(6): 861-870. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=20200606&flag=1 何文娜, 朱长青, 李仰春, 等. 智能地质体综合技术在地质图缩编中的应用[J]. 地质通报, 2020, 39(6): 871-879. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=20200607&flag=1 何文娜, 朱长青, 王永志, 等. 智绘地质——一种基于MapGIS的地质编图智能化系统[J]. 地球物理学进展, 2019, 34(5): 2030-2036. https://www.cnki.com.cn/Article/CJFDTOTAL-DQWJ201905041.htm Tomlin C D. GIS and Cartographic modeling[M]. California: Esri Press, 2012.
BGS. DiGMapGB digital geological maps of Great Britain Datasets[EB/OL]. [2018-10-24](2019-08-01). http://www.bgs.ac.uk/products/digitalmaps/DiGMapGBMaps.html, 2019.
OneGeology. Portal[EB/OL]. [2018-10-24](2019-08-01). http://www.onegeology.org/portal/home.html, 2019.
Smirnoff A, Huot-Vezina G, Paradis S J, et al. Generalizing geological maps with the GeoScaler software: The case study approach[J]. Computers & Geosciences, 2012, (40): 66-86. http://www.sciencedirect.com/science/article/pii/S0098300411002421
李廷栋. 当前国际地质编图的一些动态[J]. 地质论评, 2013, 2: 208, 216. https://www.cnki.com.cn/Article/CJFDTOTAL-DZLP201302003.htm 叶天竺, 黄崇轲, 邓志奇. 1: 250万中华人民共和国数字地质图空间数据库[J]. 中国地质, 2017, 44(S1): 19-24. https://kns.cnki.net/KCMS/detail/detail.aspx?dbname=CJFDPREP&filename=DIZI2017S1004 庞健峰, 丁孝忠, 韩坤英, 等. 1: 100万中华人民共和国数字地质图空间数据库[J]. 中国地质, 2017, 44(S1): 8-18. https://kns.cnki.net/KCMS/detail/detail.aspx?dbname=CJFDPREP&filename=DIZI2017S1003 中国地质调查局基础调查部. 中国小比例尺地质图概况[J]. 中国地质, 2018, 45(2): 423. https://www.cnki.com.cn/Article/CJFDTOTAL-DIZI201802024.htm 任纪舜, 牛宝贵, 王军, 等. 1: 500万国际亚洲地质图[J]. 地球学报, 2013, 34(1): 24-30. https://www.cnki.com.cn/Article/CJFDTOTAL-DQXB201301004.htm 于萍萍, 陈建平, 柴福山, 等. 基于地质大数据理念的模型驱动矿产资源定量预测[J], 地质通报, 2015, 34(7): 1333-1343. doi: 10.3969/j.issn.1671-2552.2015.07.011 Charlie F. Modeling active database-driven cartography within GIS databases[C]//The 21st International Cartographic Conference of the International Cartographic Association. Durban, South Africa, 2003.
Buckly A, Frye C, Buttenfield B. An Information Model for Maps: Towards Cartographic Production from GIS Databases[EB/OL]. [2018-10-24](2019-08-01). http://icaci.org/files/documents/ICC-proceedings/ICC2005/htm/pdf/oral/TEMA1/SESSION1/AILEEN%20BUCKLEY.pdf, 2019.
陈刚, 姚仲友, 王天刚, 等. 大洋洲地区系列地质图件编制[J]. 地质通报, 2014, 33(2/3): 159-163. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=2014020302&flag=1 王家耀, 钱海忠. 制图综合知识及其应用[J]. 武汉大学学报(信息科学版), 2006, 5: 382-386, 439. https://www.cnki.com.cn/Article/CJFDTOTAL-WHCH200605002.htm 王家耀, 孙群, 王光霞, 等. 地图学原理与方法[M]. 北京: 科学出版社, 2014. van Gasselt S, Nass A. Planetary mapping——The data model's perspective and GIS framework[J]. Planetary and Space Science, 2010, (59): 12. http://www.sciencedirect.com/science/article/pii/S0032063310002849
陈军, 王东华, 商瑶玲, 等. 国家1: 50 000数据库更新工程总体设计研究与技术创新[J]. 测绘学报, 2010, 39(1): 7-10. https://www.cnki.com.cn/Article/CJFDTOTAL-CHXB201001005.htm 王杨刚, 郝丽荣, 黄辉, 等. 基于空间数据和专家知识驱动的地质编图技术研究与应用[J]. 地质通报, 2019, 38(12): 2067-2076. http://dzhtb.cgs.cn/gbc/ch/reader/view_abstract.aspx?file_no=20191215&flag=1 赵鹏大. 地质大数据特点及其合理开发利用[J], 地学前缘, 2019, 26(4): 1-5. https://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201904002.htm 黄崇轲, 李浩川. 中华人民共和国1: 50万数字地质图和空间元数据[J]. 中国地质, 2001, 28(2): 43-48. doi: 10.3969/j.issn.1000-3657.2001.02.008 王家耀, 孙力楠, 成毅. 创新思维改变地图学[J]. 地理空间信息, 2011, 9(2): 1-5. doi: 10.3969/j.issn.1672-4623.2011.02.001 罗建民, 张旗. 大数据开创地学研究新途径: 查明相关关系, 增强研究可行性[J]地学前缘, 2019, 26(4): 6-12. https://www.cnki.com.cn/Article/CJFDTOTAL-DXQY201904003.htm 王家耀, 邹建华. 地图制图数据处理的模型方法[M]. 北京: 解放军出版社, 1992. 王家耀. 地图学与地理信息工程研究[M]. 北京: 科学出版社, 2005. 王光霞. 数字环境下制图综合概念和方法的拓展[J]. 测绘学院学报, 2005, 9: 207-211. doi: 10.3969/j.issn.1673-6338.2005.03.017 邬宽廉, 左群超, 邬可筠. 建立地质图数据库的地质年代代码问题[J]. 地质通报, 2008, 27(7): 1103-1106. doi: 10.3969/j.issn.1671-2552.2008.07.025 DD2006-06. 中国地质调查局地质调查技术标准[S]. 中国地质调查局, 2006. 黄辉, 路彦明, 张大可, 等. 基础地质图更新编制技术研究与数据库建设报告. 廊坊: (原)武警黄金地质研究所, 2015-2017. 黄辉, 路彦明, 张文华, 等. 基础地质数据更新技术流程方案. 廊坊: (原)武警黄金地质研究所, 2015-2017. 张克信, 朱云海, 樊光明. 1: 25万临夏市幅(I48C001001)区域地质调查报告(修测). 武汉: 中国地质大学(武汉)地质调查研究院, 2006.