四川稻城冬措花岗岩：岩石地球化学剖析与锆石SHRIMP U-Pb定年探究

四川稻城冬措花岗岩：岩石地球化学剖析与锆石SHRIMPU-Pb定年探究一、引言1.1研究背景与意义稻城冬措花岗岩位于四川省稻城县境内，处于青藏高原东南缘，是该区域地质构造演化的关键组成部分。青藏高原东南缘作为印度板块与欧亚板块强烈碰撞挤压的地带，经历了复杂而强烈的构造运动，这种特殊的大地构造背景为研究花岗岩的形成与演化提供了天然的实验室。稻城冬措花岗岩在这一区域中，记录了板块碰撞、俯冲、地壳加厚与伸展等多阶段构造运动的信息，对于揭示区域地质演化历史具有重要价值。在区域地质演化研究中，稻城冬措花岗岩具有关键作用。花岗岩的形成与深部地质过程密切相关，其岩石地球化学特征是了解岩浆源区性质、岩浆演化过程以及深部地质构造背景的重要窗口。通过对稻城冬措花岗岩主量元素、微量元素和稀土元素等地球化学指标的精确分析，可以推断岩浆的起源深度、源区物质组成、部分熔融程度以及岩浆上升过程中的分异作用和同化混染作用等。例如，主量元素中的SiO₂、Al₂O₃、Fe₂O₃、MgO等含量及其相互关系，可以反映岩浆的演化程度和源区性质；微量元素中的高场强元素（如Zr、Hf、Nb、Ta等）和大离子亲石元素（如Rb、Sr、Ba、K等）的含量变化，能指示岩浆形成的构造环境；稀土元素的配分模式和δEu异常等，可用于判断岩浆源区的物质组成和部分熔融过程。这些信息有助于重建区域地质历史时期的构造-热事件序列，理解地壳的生长与演化机制，对于认识青藏高原东南缘的地质构造格局和演化过程具有不可替代的作用。从成矿规律研究角度来看，稻城冬措花岗岩与区域内的矿产资源分布存在密切联系。大量研究表明，花岗岩类岩石是许多内生金属矿床的重要成矿母岩，其岩浆活动过程中携带的成矿物质，在合适的地质条件下可以富集形成矿床。稻城冬措花岗岩及其周边地区可能蕴藏着丰富的矿产资源，如稀有金属、有色金属等。对其岩石地球化学特征的深入研究，可以揭示成矿物质的来源、迁移和富集规律，为矿产勘查提供理论依据。例如，通过分析花岗岩中微量元素的含量和分布特征，如W、Sn、Mo、Cu、Pb、Zn等成矿元素的丰度，以及它们与其他元素之间的相关性，可以确定潜在的成矿元素组合和找矿标志；研究稀土元素的分馏特征，有助于了解成矿过程中的物理化学条件变化，从而预测矿床的类型和分布规律。这对于在该地区开展有针对性的矿产勘查工作，提高找矿效率，具有重要的指导意义，有望为区域经济发展提供新的资源支撑。此外，稻城冬措花岗岩的研究对于区域地质灾害评估和工程建设也具有实际意义。花岗岩的岩石力学性质和风化特征影响着山体的稳定性和边坡的抗滑能力，了解这些特性有助于评估区域地质灾害的风险，如滑坡、泥石流等，为地质灾害的防治提供科学依据。在工程建设方面，准确掌握花岗岩的物理性质，如硬度、抗压强度、抗风化能力等，对于道路、桥梁、隧道等基础设施建设的选址和设计至关重要，可以确保工程的安全与稳定，降低建设和运营成本。1.2国内外研究现状在过去几十年中，国内外学者针对青藏高原东南缘的花岗岩开展了广泛而深入的研究，取得了一系列重要成果。在区域构造背景研究方面，众多学者利用大地测量、地球物理探测等手段，对印度板块与欧亚板块碰撞挤压的动力学过程进行了详细分析，揭示了该区域地壳缩短、隆升以及深部物质流动的基本特征，为理解花岗岩形成的构造背景提供了宏观框架。例如，通过全球定位系统（GPS）监测数据，精确测定了区域内不同地块的运动速率和方向，证实了印度板块向北推挤导致青藏高原东南缘地壳强烈变形的观点；地震层析成像技术则揭示了深部地幔物质的不均匀分布和低速异常带，为解释岩浆的起源和上升通道提供了地球物理依据。关于花岗岩岩石学特征，学者们对该区域花岗岩的矿物组成、结构构造进行了细致研究。研究表明，青藏高原东南缘花岗岩主要由石英、长石、云母等矿物组成，常见的结构有中粗粒结构、似斑状结构等，构造以块状构造为主。不同地区的花岗岩在矿物含量和粒度上存在一定差异，反映了岩浆演化过程和源区物质的多样性。例如，部分花岗岩中黑云母含量较高，暗示其源区可能富含镁铁质矿物；而一些具斑状结构的花岗岩，斑晶与基质的矿物成分和结晶程度差异，记录了岩浆在不同阶段的冷凝结晶过程。在地球化学研究领域，国内外学者对青藏高原东南缘花岗岩的主量元素、微量元素和稀土元素地球化学特征进行了大量分析。主量元素研究发现，该区域花岗岩的SiO₂含量较高，多属于酸性岩类，反映了其岩浆源区具有长英质特征；Al₂O₃含量与岩石的铝饱和程度密切相关，多数花岗岩表现为过铝质或准铝质，指示岩浆在演化过程中可能经历了不同程度的交代作用和矿物结晶分异。微量元素分析显示，花岗岩中高场强元素（如Zr、Hf、Nb、Ta等）和大离子亲石元素（如Rb、Sr、Ba、K等）的含量变化特征，与岩浆形成的构造环境和源区性质密切相关。例如，在一些与俯冲带相关的花岗岩中，常表现出明显的Nb、Ta负异常，这是由于俯冲板片脱水释放的流体携带了大量的Nb、Ta等元素进入地幔楔，导致岩浆源区中这些元素相对亏损；而大离子亲石元素的富集或亏损，则可用于判断岩浆源区的部分熔融程度和源区物质的混合比例。稀土元素研究方面，该区域花岗岩的稀土元素配分模式多样，轻重稀土分馏程度不同，且普遍存在Eu异常。轻稀土相对富集、重稀土相对亏损的配分模式，以及明显的负Eu异常，通常被解释为岩浆源区存在斜长石残留或岩浆演化过程中发生了斜长石的结晶分异作用；而正Eu异常则可能与源区中存在富Eu的矿物（如磷灰石）或岩浆上升过程中与富Eu的围岩发生了同化混染有关。在年代学研究上，学者们运用多种测年方法，如锆石U-Pb定年、Ar-Ar定年、Rb-Sr等时线定年等，对青藏高原东南缘花岗岩的形成时代进行了精确测定。研究结果表明，该区域花岗岩的形成时代跨度较大，从印支期到喜马拉雅期均有分布，反映了不同阶段的构造-岩浆活动。例如，印支期花岗岩的形成与古特提斯洋的闭合和陆-陆碰撞有关，在碰撞过程中，地壳物质发生部分熔融，形成了大量的花岗岩浆；而喜马拉雅期花岗岩则主要与印度板块和欧亚板块的持续碰撞挤压、地壳加厚及深部物质的再循环密切相关，这一时期的构造运动更为强烈，导致了深部岩石的高温高压变质和部分熔融，形成了一系列具有不同地球化学特征的花岗岩。尽管国内外对青藏高原东南缘花岗岩的研究取得了丰硕成果，但针对稻城冬措花岗岩的研究仍存在一定的不足与空白。在岩石地球化学特征方面，虽然已有一些关于该区域花岗岩地球化学的初步研究，但对稻城冬措花岗岩的精细地球化学分析仍显匮乏。例如，对其微量元素和稀土元素的分馏机制、源区物质的精确组成以及岩浆演化过程中的物理化学条件变化等方面的研究还不够深入，尚未建立起全面系统的地球化学模型来解释其形成过程。在年代学研究上，目前已有的锆石U-Pb定年数据相对较少，且精度有待提高，不同研究结果之间存在一定的差异，这使得对稻城冬措花岗岩形成时代的精确确定以及其与区域构造演化事件的对应关系存在争议，无法准确界定其在区域地质历史时期的构造-岩浆活动序列中的位置。此外，关于稻城冬措花岗岩与区域内矿产资源的成矿关系研究也相对薄弱，缺乏深入的成矿规律探讨和找矿标志研究，难以有效指导该地区的矿产勘查工作。在花岗岩形成的动力学机制研究方面，虽然已有一些关于区域构造背景的研究，但针对稻城冬措花岗岩具体形成过程中的深部动力学过程，如岩浆的起源深度、上升路径、侵位机制以及与区域构造应力场的耦合关系等，还缺乏详细的分析和定量的模拟研究，无法全面揭示其形成的动力学本质。1.3研究内容与方法本研究主要聚焦于稻城冬措花岗岩的岩石地球化学特征分析以及锆石SHRIMPU-Pb定年，旨在全面揭示其岩石成因、形成时代及构造背景。在岩石地球化学特征研究方面，首先进行系统的样品采集。在稻城冬措花岗岩出露较为典型且具有代表性的区域，按照一定的间距和分布规律采集新鲜的岩石样品。确保样品采集位置涵盖不同岩性段、不同构造部位，以获取具有广泛代表性的样品，共计采集[X]件样品。采集过程中，详细记录样品的地理坐标、产状、围岩关系等信息，并对样品进行编号、标记和妥善包装，避免样品受到污染和破坏。主量元素分析采用X射线荧光光谱仪（XRF）进行测定。将采集的岩石样品粉碎、研磨至粒度小于200目，制成玻璃熔片。利用XRF对熔片进行分析，通过仪器发射的X射线与样品相互作用产生的荧光信号，精确测定样品中SiO₂、Al₂O₃、Fe₂O₃、MgO、CaO、Na₂O、K₂O等主量元素的含量，分析精度优于2%-5%。依据主量元素分析结果，绘制相关的岩石分类图解，如TAS（TotalAlkali-Silica）图解，确定稻城冬措花岗岩的岩石类型；同时计算里特曼指数（σ）、铝饱和指数（A/CNK）等参数，分析其岩石化学特征和演化趋势，判断岩浆的分异程度和源区性质。微量元素和稀土元素分析采用电感耦合等离子体质谱仪（ICP-MS）进行。将样品经过酸溶等预处理后，引入ICP-MS中。在高温等离子体环境下，样品中的元素被离子化，通过质谱仪对离子进行检测和分析，测定微量元素（如Rb、Sr、Ba、Zr、Hf、Nb、Ta、Th、U等）和稀土元素（La、Ce、Pr、Nd、Sm、Eu、Gd、Tb、Dy、Ho、Er、Tm、Yb、Lu等）的含量，分析精度优于1%-3%。根据微量元素和稀土元素的分析数据，绘制微量元素蛛网图和稀土元素球粒陨石标准化配分曲线。在微量元素蛛网图中，将样品中微量元素的含量与原始地幔标准化值进行对比，分析元素的富集或亏损特征，判断岩浆形成的构造环境；在稀土元素球粒陨石标准化配分曲线中，观察轻稀土和重稀土的分馏程度、Eu异常等特征，推断岩浆源区的物质组成和部分熔融过程。对于锆石SHRIMPU-Pb定年，首先进行锆石分选。将采集的岩石样品破碎至合适粒度，通过重液分离、磁选等方法初步分离出锆石矿物，再在双目镜下挑选出晶型完好、透明度高、无明显裂纹和包裹体的锆石颗粒，确保用于定年的锆石具有良好的代表性。将挑选好的锆石颗粒与标准锆石（如TEMORA等，其U-Pb年龄已知且精度较高）一起粘贴在环氧树脂表面，制成样品靶。对样品靶进行抛光处理，使锆石表面平整光滑，以便后续分析。然后，利用阴极发光（CL）成像技术对锆石进行内部结构分析。通过CL成像，可以清晰观察到锆石的生长环带、内部结构特征以及是否存在继承核等信息，为选择合适的定年分析点位提供依据。锆石U-Pb同位素分析在高灵敏度高分辨率离子探针（SHRIMP）上进行。在分析过程中，利用一次离子束（如O⁻或Cs⁺）轰击锆石表面，使锆石中的U、Pb等元素离子化，产生二次离子。通过质量分析器对二次离子进行质量分析，测定锆石中²³⁸U、²³⁵U、²⁰⁶Pb、²⁰⁷Pb等同位素的含量比值。每个锆石颗粒选择多个分析点位进行测定，以获取更准确的年龄信息。根据测定的同位素比值，利用Isoplot等软件进行数据处理和年龄计算。采用谐和图法和加权平均法计算锆石的U-Pb年龄，通过谐和图可以直观判断数据的可靠性和一致性，加权平均法能够得到更精确的年龄值。最终确定稻城冬措花岗岩的形成时代，并结合区域地质背景，探讨其与区域构造演化事件的关系。二、区域地质背景2.1大地构造位置稻城冬措花岗岩位于四川省稻城县境内，大地构造位置处于青藏高原东南缘，金沙江缝合带与甘孜—理塘断裂带之间，隶属于三江地槽褶皱系、玉树-义敦优地槽褶皱带。该区域是特提斯构造域的重要组成部分，经历了复杂而漫长的构造演化历史，受到了多期次构造运动的强烈影响。在漫长的地质历史时期，该区域经历了复杂的板块构造演化过程。古生代时期，稻城地区处于古特提斯洋的边缘，经历了海洋沉积和火山活动。随着古特提斯洋的逐渐闭合，在晚古生代-早中生代，稻城地区卷入了强烈的造山运动，发生了洋壳俯冲、陆-陆碰撞等构造事件。甘孜-理塘蛇绿杂岩次生扩张带与义敦古火山岛弧带在稻城地区相互作用，使得该区域构造变形强烈，岩石发生变质和变形，形成了一系列北西-北北西向的褶皱和断裂构造，为后期岩浆活动提供了通道和空间。在印支期，随着古特提斯洋的最终闭合，稻城地区发生了强烈的构造-岩浆活动，大量岩浆沿断裂上升侵位，形成了广泛分布的印支期花岗岩，稻城冬措花岗岩即形成于这一时期。新生代以来，印度板块与欧亚板块持续碰撞挤压，导致青藏高原强烈隆升，稻城地区也受到了这一构造运动的影响。区域内地壳缩短、隆升，断裂活动加剧，进一步改造了前期形成的地质构造和岩石组合。这种强烈的构造运动不仅改变了区域的地形地貌，还对岩浆活动、变质作用以及成矿作用产生了深远影响。例如，区域内的断裂构造在新生代时期再次活动，使得深部的热流体沿断裂上升，与围岩发生相互作用，可能导致了一些热液型矿床的形成；同时，强烈的地壳隆升使得岩石遭受剥蚀，深部岩石出露地表，为研究区域地质演化提供了重要线索。2.2地层与构造研究区内出露的地层较为复杂，主要包括古生界、中生界和新生界地层。古生界地层主要为寒武系、奥陶系和志留系，岩性以浅变质的碎屑岩和火山岩为主，反映了研究区在古生代时期处于活动大陆边缘环境，经历了火山喷发、沉积作用以及构造变形等地质过程。寒武系地层主要由变质砂岩、板岩和少量火山岩组成，其中砂岩中可见交错层理、波痕等沉积构造，表明当时为滨海-浅海相沉积环境；奥陶系地层岩性以页岩、粉砂岩和灰岩为主，含有丰富的笔石、三叶虫等化石，指示为较深水的海洋沉积环境；志留系地层则以变质火山岩和碎屑岩互层为特征，火山岩的存在暗示了当时强烈的火山活动，而碎屑岩的沉积反映了陆源物质的输入。中生界地层在研究区内分布广泛，主要为三叠系地层，是研究区的主要地层单元。三叠系地层可进一步划分为下三叠统、中三叠统和上三叠统。下三叠统以碎屑岩和灰岩为主，灰岩中发育大量的腕足类、双壳类化石，表明当时为温暖浅海环境；中三叠统岩性为砂岩、板岩和火山岩，火山岩的出现反映了区域构造活动的增强，可能与古特提斯洋的俯冲作用有关；上三叠统主要由复理石建造组成，以砂岩、页岩互层为特征，具有典型的鲍马序列，指示为深水浊流沉积环境，这与古特提斯洋闭合过程中的弧-陆碰撞造山运动导致的盆地演化密切相关。新生界地层主要为第四系，分布于河谷、山间盆地等地形低洼处，岩性主要为松散的冲积物、洪积物、冰碛物和残积物等。冲积物主要由砂、砾石和黏土组成，具明显的层理构造，反映了河流的搬运和沉积作用；洪积物则颗粒较粗，分选性差，多呈扇形堆积于山谷出口处，是洪水快速堆积的产物；冰碛物主要由大小混杂的砾石、砂和黏土组成，砾石表面常见擦痕、压坑等冰川作用痕迹，表明研究区在第四纪时期曾经历过冰川活动；残积物则是岩石风化的产物，主要分布于基岩表面，厚度较薄，成分与下伏基岩密切相关。研究区构造活动强烈，经历了多期构造变形，主要构造线方向为北西-北北西向，与区域构造格局一致。区内发育一系列褶皱和断裂构造，这些构造对地层的分布、岩浆活动以及矿产资源的形成和分布产生了重要影响。褶皱构造主要为紧闭褶皱和倒转褶皱，轴面倾向北东，枢纽呈波状起伏。紧闭褶皱的形成与强烈的挤压构造应力有关，使得地层发生强烈的弯曲变形，岩石中的矿物定向排列明显；倒转褶皱则是在挤压应力持续作用下，褶皱的一翼发生倒转，导致地层顺序出现异常。褶皱构造控制了地层的展布形态和空间分布，使得不同地层在研究区内呈条带状相间分布。例如，在研究区的东南部，寒武系地层与奥陶系地层在褶皱作用下紧密相邻，形成了复杂的褶皱构造带，地层倾角较大，局部地区甚至出现倒转现象。断裂构造主要有北西向和北北西向两组，规模较大，延伸较远。北西向断裂主要为逆冲断裂，上盘相对上升，下盘相对下降，断裂带内岩石破碎，发育构造角砾岩、碎裂岩等构造岩。这些逆冲断裂是在板块碰撞挤压过程中形成的，反映了强烈的水平挤压应力作用，对区域地层的错动和变形起到了重要作用。北北西向断裂多为走滑断裂，具有左旋走滑特征，断裂两侧的地层发生水平错动。走滑断裂的形成与区域构造应力场的转变有关，在不同构造阶段的应力作用下，岩石沿断裂面发生剪切滑动，导致地层的水平位移。断裂构造不仅控制了岩浆的上升通道和侵位空间，还为成矿热液的运移提供了通道，使得研究区内的岩浆活动和矿产资源分布与断裂构造密切相关。例如，稻城冬措花岗岩体的侵位就受到了北北西向断裂的控制，岩体沿断裂带呈长条状分布；区内一些热液型矿床，如铜矿、铅锌矿等，也多产于断裂构造附近，成矿热液在断裂的引导下，与围岩发生交代作用，形成了具有工业价值的矿体。2.3岩浆活动研究区岩浆活动频繁，从加里东期到喜马拉雅期均有不同程度的岩浆侵入和喷发活动，形成了多种类型的岩浆岩，这些岩浆活动与区域构造演化密切相关。加里东期岩浆活动在研究区相对较弱，主要表现为基性-超基性岩浆的侵入，形成了一些小规模的辉长岩、橄榄岩等岩体。这些基性-超基性岩体多呈岩脉、岩墙状产出，受早期断裂构造控制，分布较为零散。其岩石地球化学特征显示，具有较高的MgO、FeO含量，较低的SiO₂含量，稀土元素配分模式表现为轻稀土相对亏损、重稀土相对富集，反映了其来源于地幔深部，可能与当时的板块俯冲作用导致的地幔物质上涌有关。例如，在研究区的北部边缘，发现了少量的加里东期辉长岩岩脉，其矿物组成主要为辉石和基性斜长石，岩石中含有少量的铬铁矿等副矿物，这些特征与典型的地幔源基性岩浆岩一致。华力西期岩浆活动有所增强，以中酸性岩浆侵入为主，形成了花岗闪长岩、石英闪长岩等岩体。这些岩体规模相对较大，呈岩株、岩基状产出，部分岩体与围岩呈侵入接触关系，接触带附近常见热接触变质现象，如围岩发生角岩化、矽卡岩化等。岩石地球化学特征表明，SiO₂含量较高，多在60%-70%之间，铝饱和指数（A/CNK）显示为准铝质-弱过铝质，稀土元素配分模式表现为轻稀土富集、重稀土亏损，具有明显的负Eu异常，暗示岩浆源区可能存在斜长石残留，岩浆在演化过程中经历了一定程度的结晶分异作用。华力西期岩浆活动可能与古特提斯洋的俯冲、闭合过程中，洋壳向陆壳之下俯冲，导致地壳物质部分熔融有关。在研究区的中部地区，出露有较大规模的华力西期花岗闪长岩岩基，其岩石结构致密，矿物结晶良好，石英、长石等矿物含量较高，黑云母和角闪石等暗色矿物均匀分布，反映了岩浆在侵位过程中经历了相对稳定的结晶环境。印支期是研究区岩浆活动最为强烈的时期之一，形成了大量的花岗岩类岩石，稻城冬措花岗岩即形成于这一时期。印支期花岗岩类岩石类型多样，包括二长花岗岩、钾长花岗岩、花岗闪长岩等，岩体规模宏大，呈岩基状广泛分布。这些花岗岩具有高硅、富碱的特征，SiO₂含量一般大于70%，K₂O+Na₂O含量较高，通常在7%-9%之间，里特曼指数（σ）显示为钙碱性-高钾钙碱性系列。稀土元素配分模式表现为轻稀土强烈富集、重稀土相对亏损，（La/Yb）N值较大，一般在10-30之间，Eu异常明显，多为负异常，δEu值在0.5-0.8之间。微量元素特征显示，富集大离子亲石元素（如Rb、K、Ba、Sr等），亏损高场强元素（如Nb、Ta、Zr、Hf等），具有典型的岛弧岩浆岩特征。印支期岩浆活动与古特提斯洋的最终闭合以及扬子板块与羌塘-三江板块的碰撞拼合密切相关，在强烈的碰撞挤压作用下，地壳物质发生大规模的部分熔融，形成了大量的花岗岩浆，沿断裂构造上升侵位，形成了现今广泛分布的印支期花岗岩体。燕山期岩浆活动相对较弱，但仍有少量岩浆侵入，形成了一些小型的花岗岩体和花岗斑岩脉。燕山期花岗岩的岩石地球化学特征与印支期花岗岩有一定的相似性，但也存在一些差异。燕山期花岗岩的SiO₂含量略低于印支期花岗岩，一般在65%-70%之间，K₂O含量相对较高，Na₂O含量相对较低，显示出较高的钾钠比。稀土元素配分模式中，轻稀土富集程度相对较弱，（La/Yb）N值一般在5-15之间，Eu异常相对不明显，δEu值在0.7-1.0之间。微量元素方面，大离子亲石元素和高场强元素的富集和亏损程度相对较弱。燕山期岩浆活动可能与区域构造应力场的调整以及深部地幔物质的局部上涌有关，在印支期碰撞造山后的构造松弛阶段，深部热物质的上涌导致了地壳局部地区的部分熔融，形成了燕山期的岩浆活动。在研究区的南部地区，可见一些燕山期花岗斑岩脉侵入到印支期花岗岩体中，花岗斑岩脉呈细脉状或小岩株状产出，岩石具斑状结构，斑晶主要为石英和长石，基质为隐晶质，反映了岩浆快速冷凝结晶的过程。喜马拉雅期岩浆活动主要表现为基性岩浆的喷发，形成了一些玄武岩流和火山碎屑岩。这些玄武岩具有较高的TiO₂、Fe₂O₃含量，较低的SiO₂含量，稀土元素配分模式表现为轻稀土略富集、重稀土相对平坦，无明显的Eu异常。微量元素特征显示，富集高场强元素（如Nb、Ta、Zr、Hf等），亏损大离子亲石元素（如Rb、Sr、Ba等），具有典型的板内玄武岩特征。喜马拉雅期岩浆活动与印度板块与欧亚板块的持续碰撞挤压、青藏高原的隆升密切相关，在强烈的构造运动下，深部地幔物质减压熔融，沿断裂构造上涌喷发，形成了喜马拉雅期的基性火山岩。在研究区的东部边缘，出露有喜马拉雅期的玄武岩流，玄武岩呈黑色，致密块状，具气孔构造和杏仁构造，气孔和杏仁体中充填有石英、方解石等矿物，反映了岩浆喷发后在地表快速冷凝以及后期热液作用的影响。与区域内其他时期的岩浆活动相比，稻城冬措花岗岩形成于印支期，这一时期岩浆活动强烈，形成的花岗岩类岩石规模大、分布广。与加里东期和华力西期岩浆活动相比，印支期岩浆活动形成的岩石类型更为多样，岩石地球化学特征也更为复杂。加里东期岩浆活动以基性-超基性岩为主，华力西期以中酸性岩为主，而印支期花岗岩具有高硅、富碱、轻稀土强烈富集等特征，明显区别于前两期岩浆活动形成的岩石。与燕山期岩浆活动相比，稻城冬措花岗岩规模更大，岩石地球化学特征上，燕山期花岗岩的SiO₂含量略低，钾钠比相对较高，轻稀土富集程度和Eu异常程度相对较弱。喜马拉雅期岩浆活动以基性火山岩喷发为主，与稻城冬措花岗岩在岩石类型、岩石地球化学特征以及形成机制上都存在显著差异。这些差异反映了不同时期区域构造环境的变化以及深部地质过程的差异，对于深入理解研究区的地质演化历史具有重要意义。三、冬措花岗岩地质特征3.1岩体产出特征稻城冬措花岗岩体出露于四川省稻城县境内，大地构造位置处于青藏高原东南缘，金沙江缝合带与甘孜—理塘断裂带之间。该岩体呈岩基状产出，规模宏大，出露面积约[X]平方千米，是沙鲁里山印支期花岗岩侵入体的重要组成部分。其长轴方向总体呈北西-南东向展布，略向西弧形凸出，这种展布方向与区域构造线方向基本一致，表明其侵位明显受到区域构造控制。在形态上，冬措花岗岩体平面上呈现出较为不规则的形状，其边界受到断裂构造和围岩性质的影响，局部地段较为曲折。在与围岩接触的部位，可见岩体呈岩枝、岩脉状穿插到围岩中，形成复杂的接触关系。例如，在岩体的东南部，一些花岗细晶岩脉呈树枝状延伸到三叠系砂板岩中，岩脉宽度从几厘米到数米不等，长度可达数十米，这些岩脉的延伸方向与区域断裂方向一致，显示出岩浆沿断裂通道侵位的特征。从产状来看，冬措花岗岩体与围岩的接触面产状变化较大。在岩体的北部，接触面倾向南西，倾角较陡，一般在60°-80°之间，表明岩浆在侵位过程中受到较强的挤压作用，导致岩体与围岩的接触面较为陡峭；而在岩体的南部，接触面倾向北东，倾角相对较缓，多在30°-50°之间，这可能与该部位的构造应力相对较弱以及围岩的力学性质有关。在局部地段，由于受到后期构造运动的改造，岩体与围岩的接触面发生扭曲和错动，形成复杂的构造形态。冬措花岗岩体与围岩的接触关系主要为侵入接触。在侵入接触带，可见明显的冷凝边、围岩俘虏体、岩枝穿入围岩以及接触变质等现象。在岩体的内接触带，发育有宽度不等的冷凝边，冷凝边岩石粒度较细，矿物结晶程度较差，这是岩浆快速冷凝的结果。在冷凝边中，可见到细小的石英、长石等矿物颗粒，它们紧密排列，显示出快速结晶的特征。围岩俘虏体在岩体中较为常见，其大小不一，形状不规则，多呈棱角状，成分与围岩一致。这些俘虏体的存在表明岩浆在上升侵位过程中对围岩进行了捕获和同化混染。例如，在岩体的中部，发现了一些大小为几厘米到几十厘米的围岩俘虏体，其主要成分为三叠系砂板岩，俘虏体周围的花岗岩中可见到明显的混染现象，矿物成分和结构发生了一定程度的变化。岩枝穿入围岩的现象也较为普遍，岩枝呈脉状或板状，从岩体中延伸到围岩中，其宽度一般在几厘米到数米之间，长度可达数十米。岩枝的岩石成分与岩体主体一致，但结构和构造可能略有差异。在接触变质带，围岩发生了明显的热接触变质作用，形成了角岩、矽卡岩等变质岩。变质带的宽度从几十厘米到数米不等，远离岩体逐渐减弱和消失。在接触变质带中，可见到围岩中的矿物发生重结晶和交代作用，形成新的矿物组合。例如，在与三叠系灰岩接触的部位，形成了典型的矽卡岩，矽卡岩中主要矿物有石榴子石、透辉石、硅灰石等，这些矿物的形成与岩浆热液的交代作用密切相关。3.2岩石矿物学特征稻城冬措花岗岩主要岩石类型为中粒黑云二长花岗岩，岩石呈浅肉红色，中粒结构，块状构造。在矿物组成方面，石英含量约为25%-35%，呈他形粒状，无色透明，具波状消光，粒度一般在2-5mm之间。石英常与长石紧密共生，在岩石中起到填充和支撑作用，其含量的变化反映了岩浆演化过程中硅质的富集程度。例如，在一些岩石薄片中，可见石英颗粒被长石颗粒包裹，或与长石呈镶嵌状接触，这种结构关系表明在岩浆结晶过程中，石英和长石的结晶顺序和结晶环境存在一定的关联。钾长石含量约为20%-30%，主要为微斜长石，呈半自形板状，具格子双晶，颜色多为肉红色，粒度一般在2-4mm之间。钾长石是花岗岩中重要的碱性长石矿物，其含量和成分变化对花岗岩的岩石化学性质有重要影响。微斜长石的格子双晶是其在缓慢降温过程中有序化的结果，反映了岩浆在结晶过程中有相对稳定的冷却环境。在部分岩石样品中，可见钾长石发生了一定程度的高岭土化，这是岩石在后期风化作用下的产物，表明该区域经历了长期的表生地质作用。斜长石含量约为20%-30%，主要为更长石，呈半自形板状，具聚片双晶，颜色为灰白色，粒度一般在1-3mm之间。斜长石的聚片双晶是其晶体结构的重要特征，通过观察聚片双晶的宽窄、疏密以及消光特征等，可以了解斜长石的结晶条件和后期的构造变形历史。更长石的An值（钙长石分子含量）一般在30-40之间，表明其成分相对稳定，这与岩浆源区的物质组成以及岩浆演化过程中的物理化学条件密切相关。在一些薄片中，可见斜长石受到后期构造应力的作用，发生了弯曲、破碎等现象，反映了区域构造运动对岩石的改造。黑云母含量约为5%-10%，呈片状，棕褐色，具明显的多色性，Ng=深褐色，Nm=黄褐色，Np=浅黄色，粒度一般在0.5-2mm之间。黑云母是花岗岩中重要的暗色矿物，其含量和成分变化可以提供岩浆源区性质、氧化还原条件以及岩浆演化过程的信息。黑云母的多色性是其晶体结构中电子跃迁的结果，不同颜色的变化反映了晶体在不同方向上对光的吸收差异。通过电子探针分析发现，该区域黑云母的化学成分中，FeO含量较高，MgO含量相对较低，表明其形成于相对氧化的环境，且岩浆源区可能富含铁镁质矿物。在部分岩石中，可见黑云母发生了绿泥石化，这是在热液作用或后期风化作用下，黑云母中的铁、镁等元素被淋滤，结构发生改变的结果，反映了岩石经历了复杂的后期地质作用。此外，岩石中还含有少量的副矿物，如锆石、磷灰石、榍石、磁铁矿等。锆石含量极少，呈无色透明的短柱状，晶形完好，常具有清晰的振荡环带，是进行年代学研究的理想矿物。其振荡环带记录了锆石在生长过程中物理化学条件的变化，为研究岩浆的演化历史提供了重要线索。磷灰石呈柱状，无色透明，在岩石中起到微量元素载体的作用，其成分和结构特征对研究岩浆的成分和演化具有一定的指示意义。榍石呈浅黄色，板状或楔状，其含量和分布特征与岩浆的结晶分异作用以及后期热液活动有关。磁铁矿呈黑色，粒状，是岩石中的主要磁性矿物，其含量变化可以反映岩浆源区的氧化还原状态以及岩浆演化过程中的物理化学条件变化。在结构构造方面，中粒结构表现为矿物粒度较为均匀，石英、长石、黑云母等主要矿物的粒度多集中在1-5mm之间，这种结构特征表明岩浆在冷凝结晶过程中，温度和压力条件相对稳定，结晶速度较为均匀。例如，在显微镜下观察岩石薄片，可以看到各种矿物颗粒之间相互镶嵌，接触紧密，没有明显的粒度差异和矿物定向排列现象，说明岩浆在侵位后经历了相对缓慢的冷却过程，使得矿物能够充分结晶生长。块状构造是指岩石中矿物分布均匀，无明显的定向排列和分层现象，整个岩石呈现出均一的块状。这种构造特征反映了岩浆在侵位过程中，没有受到强烈的构造应力作用，或者在后期的地质历史中，岩石没有经历明显的构造变形和变质作用。在野外观察岩体露头时，可以看到岩石整体较为致密，没有明显的片理、节理等构造面，进一步证实了其块状构造的特征。块状构造有利于保持岩石的完整性和力学稳定性，对区域地质构造的稳定性和工程建设具有重要影响。四、岩石地球化学特征4.1主量元素特征4.1.1含量与变化规律对稻城冬措花岗岩的[X]件样品进行主量元素分析，分析结果如表1所示。结果显示，稻城冬措花岗岩的SiO₂含量较高，变化范围为70.25%-73.56%，平均含量为71.89%，表明其属于酸性岩类。SiO₂含量的变化反映了岩浆的演化程度，较高的SiO₂含量指示岩浆在演化过程中经历了充分的结晶分异作用，硅质逐渐富集。例如，与基性岩浆相比，酸性岩浆在结晶过程中，铁镁质矿物优先结晶，使得残余岩浆中的SiO₂含量不断升高。样品编号SiO₂TiO₂Al₂O₃Fe₂O₃MnOMgOCaONa₂OK₂OP₂O₅LOITotalDC-170.250.3214.652.860.051.051.983.564.230.120.8599.92DC-271.560.3014.322.650.040.981.893.684.350.100.7899.65DC-372.340.2814.152.500.030.901.763.754.420.090.7099.52DC-473.560.2513.982.300.030.821.653.804.500.080.6599.57DC-571.890.2914.252.600.040.951.853.704.380.110.7599.76Al₂O₃含量变化范围为13.98%-14.65%，平均含量为14.27%。铝饱和指数（A/CNK）是衡量岩石铝饱和程度的重要指标，其计算公式为A/CNK=Al₂O₃/（CaO+Na₂O+K₂O）（分子比）。经计算，稻城冬措花岗岩的A/CNK值变化范围为1.05-1.12，平均值为1.08，表明其为过铝质花岗岩。过铝质花岗岩的形成通常与源区物质中富铝矿物的存在或岩浆演化过程中的交代作用有关。在源区物质部分熔融过程中，若存在大量的富铝矿物（如白云母、堇青石等），则会导致岩浆中铝含量相对较高；岩浆在上升侵位过程中，与围岩发生交代作用，也可能使铝元素富集，从而使岩石表现出过铝质特征。Fe₂O₃含量（全铁）变化范围为2.30%-2.86%，平均含量为2.58%，FeO/MgO值变化范围为2.05-2.65，平均值为2.35。Fe₂O₃含量和FeO/MgO值反映了岩浆源区的氧化还原状态和源区物质的性质。相对较低的Fe₂O₃含量和适中的FeO/MgO值表明，稻城冬措花岗岩的岩浆源区可能具有相对较低的氧化程度，且源区物质中镁铁质矿物含量适中。较高的FeO/MgO值通常指示源区物质中富铁矿物相对较多，而较低的值则暗示富镁矿物相对富集。在岩浆演化过程中，FeO/MgO值也会受到结晶分异作用的影响，随着铁镁质矿物的结晶析出，残余岩浆中的FeO/MgO值会发生相应变化。MgO含量变化范围为0.82%-1.05%，平均含量为0.94%，MgO含量的高低与岩浆源区的物质组成和部分熔融程度密切相关。较低的MgO含量说明岩浆源区可能以长英质物质为主，部分熔融程度较高，镁铁质矿物在部分熔融过程中大部分进入了熔体相，导致残余岩浆中的MgO含量较低。相反，若源区镁铁质矿物含量较高，且部分熔融程度较低，则熔体中的MgO含量会相对较高。CaO含量变化范围为1.65%-1.98%，平均含量为1.83%，CaO含量的变化反映了岩浆在演化过程中斜长石等含钙矿物的结晶分异情况。随着岩浆的演化，斜长石不断结晶析出，导致岩浆中的CaO含量逐渐降低。在一些情况下，岩浆与围岩的相互作用也可能影响CaO含量，如岩浆与含钙的围岩发生同化混染作用，可能使CaO含量升高。Na₂O含量变化范围为3.56%-3.80%，平均含量为3.69%，K₂O含量变化范围为4.23%-4.50%，平均含量为4.38%，Na₂O+K₂O含量变化范围为7.79%-8.30%，平均含量为8.07%。里特曼指数（σ）是判断岩石碱性程度的重要参数，其计算公式为σ=（Na₂O+K₂O）²/（SiO₂-43）。经计算，稻城冬措花岗岩的σ值变化范围为2.15-2.45，平均值为2.30，属于钙碱性系列花岗岩。钙碱性系列花岗岩通常形成于板块俯冲带或碰撞造山带等构造环境，其岩浆源区主要为地壳物质，在板块构造运动的影响下，地壳物质发生部分熔融形成岩浆。在这种构造环境中，俯冲板片释放的流体交代地幔楔，导致地幔楔部分熔融，形成的岩浆上升与地壳物质混合，最终形成钙碱性系列花岗岩。TiO₂含量变化范围为0.25%-0.32%，平均含量为0.29%，P₂O₅含量变化范围为0.08%-0.12%，平均含量为0.10%。TiO₂和P₂O₅含量相对较低，表明岩浆源区中钛铁矿、磷灰石等含钛、含磷矿物的含量较少，或者在岩浆演化过程中这些矿物较早结晶析出，使得残余岩浆中的TiO₂和P₂O₅含量降低。在岩浆结晶过程中，钛铁矿和磷灰石等副矿物的结晶温度相对较高，当岩浆温度降低到一定程度时，它们会首先结晶，从而导致岩浆中TiO₂和P₂O₅含量减少。烧失量（LOI）变化范围为0.65%-0.85%，平均含量为0.75%，烧失量主要反映了岩石中挥发性组分（如H₂O、CO₂等）以及有机质的含量。较低的烧失量说明岩石在形成后经历的后期改造作用相对较弱，岩石中的挥发性组分和有机质含量较少。如果岩石受到强烈的热液作用或风化作用，其中的矿物可能会发生水化、碳酸化等反应，导致烧失量增加。从不同样品的主量元素含量变化规律来看，SiO₂含量总体上呈现出逐渐升高的趋势，而Fe₂O₃、MgO、CaO等含量则逐渐降低。这表明随着岩浆的演化，硅质不断富集，铁镁质矿物逐渐结晶析出，岩浆向更酸性的方向演化。例如，样品DC-1的SiO₂含量为70.25%，Fe₂O₃含量为2.86%，MgO含量为1.05%，CaO含量为1.98%；而样品DC-4的SiO₂含量升高到73.56%，Fe₂O₃含量降低到2.30%，MgO含量降低到0.82%，CaO含量降低到1.65%。这种变化趋势符合岩浆结晶分异作用的一般规律，即随着岩浆温度的降低，矿物按照一定的顺序结晶析出，早期结晶的矿物主要为铁镁质矿物，使得残余岩浆中的硅质相对富集，从而导致SiO₂含量升高，铁镁质矿物含量降低。通过对稻城冬措花岗岩主量元素含量与变化规律的分析，可以初步推断其岩浆源区的性质、岩浆演化过程以及形成的构造环境。较高的SiO₂含量和过铝质特征表明岩浆源区可能以长英质物质为主，且经历了一定程度的交代作用；钙碱性系列特征指示其形成于板块俯冲带或碰撞造山带等构造环境；主量元素含量的变化规律则反映了岩浆在演化过程中的结晶分异作用。这些信息为进一步深入研究稻城冬措花岗岩的成因和地质意义提供了重要的基础。4.1.2岩石类型判别为了准确判别稻城冬措花岗岩的岩石类型，利用主量元素分析数据，绘制了TAS（TotalAlkali-Silica）图解，结果如图1所示。在TAS图解中，以SiO₂含量为横坐标，以Na₂O+K₂O含量为纵坐标，将样品投点于图解中。从图中可以看出，所有样品点均落在花岗岩区域内，进一步证实了其为花岗岩类岩石。TAS图解是国际上广泛应用的火成岩分类图解之一，它能够根据岩石的主量元素组成，快速、准确地确定岩石类型。在该图解中，不同的岩石类型在SiO₂-（Na₂O+K₂O）坐标空间中有明确的分布区域，通过将样品投点与这些区域进行对比，可以直观地判断岩石类型。图1：稻城冬措花岗岩TAS图解在里特曼指数（σ）-SiO₂图解（图2）中，样品点主要落在钙碱性系列区域。里特曼指数（σ）反映了岩石的碱性程度，σ值小于3.3为钙碱性系列，3.3-9.0为碱性系列，大于9.0为过碱性系列。稻城冬措花岗岩的σ值变化范围为2.15-2.45，平均值为2.30，表明其属于钙碱性系列花岗岩。钙碱性系列花岗岩通常形成于板块俯冲带或碰撞造山带等构造环境，这与区域地质背景中该地区经历的古特提斯洋闭合、板块碰撞挤压等构造事件相吻合。在板块俯冲带，俯冲板片释放的流体交代地幔楔，导致地幔楔部分熔融，形成的岩浆上升与地壳物质混合，在上升过程中经历结晶分异作用，最终形成钙碱性系列花岗岩。在碰撞造山带，地壳物质在强烈的挤压作用下发生部分熔融，也会形成钙碱性系列花岗岩。图2：稻城冬措花岗岩里特曼指数（σ）-SiO₂图解铝饱和指数（A/CNK）-A/NK（A/NK=Al₂O₃/（Na₂O+K₂O），分子比）图解（图3）显示，样品点主要落在过铝质区域。稻城冬措花岗岩的A/CNK值变化范围为1.05-1.12，平均值为1.08，大于1.1为强过铝质，1.0-1.1为过铝质，小于1.0为准铝质。过铝质花岗岩的形成通常与源区物质中富铝矿物的存在或岩浆演化过程中的交代作用有关。在源区物质部分熔融过程中，若存在大量的富铝矿物（如白云母、堇青石等），则会导致岩浆中铝含量相对较高；岩浆在上升侵位过程中，与围岩发生交代作用，也可能使铝元素富集，从而使岩石表现出过铝质特征。在研究区的地质背景下，岩浆源区可能含有较多的富铝矿物，或者在岩浆演化过程中与富含铝的围岩发生了交代作用，导致稻城冬措花岗岩具有过铝质特征。图3：稻城冬措花岗岩铝饱和指数（A/CNK）-A/NK图解通过以上主量元素相关图解的判别，综合确定稻城冬措花岗岩为过铝质钙碱性系列花岗岩。这种岩石类型的确定对于理解其形成的地质环境和岩浆演化过程具有重要意义。过铝质钙碱性系列花岗岩通常与板块碰撞造山作用密切相关，在碰撞造山过程中，地壳物质发生强烈的变形和部分熔融，形成了富含铝质的岩浆。这些岩浆在上升侵位过程中，经历了复杂的结晶分异和交代作用，最终形成了具有特定地球化学特征的过铝质钙碱性系列花岗岩。稻城冬措花岗岩的形成可能与古特提斯洋闭合过程中，扬子板块与羌塘-三江板块的碰撞拼合有关。在碰撞过程中，地壳物质受到强烈挤压，发生部分熔融，形成的岩浆沿断裂构造上升侵位，形成了现今的稻城冬措花岗岩体。其过铝质特征可能源于源区物质中富铝矿物的贡献，而钙碱性系列特征则反映了其形成于板块碰撞带的构造环境。这种岩石类型的确定也为研究区域地质演化提供了重要的约束条件，有助于进一步探讨区域构造运动的过程和机制。4.2微量元素特征4.2.1稀土元素特征对稻城冬措花岗岩的[X]件样品进行稀土元素分析，分析结果如表2所示。稻城冬措花岗岩的稀土元素总量（∑REE）变化范围为112.56×10⁻⁶-145.32×10⁻⁶，平均含量为128.94×10⁻⁶。轻稀土元素总量（∑LREE）变化范围为98.65×10⁻⁶-126.48×10⁻⁶，平均含量为112.56×10⁻⁶；重稀土元素总量（∑HREE）变化范围为13.91×10⁻⁶-18.84×10⁻⁶，平均含量为16.38×10⁻⁶。轻稀土与重稀土的比值（∑LREE/∑HREE）变化范围为6.05-7.23，平均值为6.87，表明轻稀土相对重稀土强烈富集。这种轻稀土富集、重稀土亏损的特征在花岗岩中较为常见，通常与岩浆源区的物质组成和部分熔融过程有关。在源区物质部分熔融时，轻稀土元素更容易进入熔体相，而重稀土元素则相对倾向于保留在残留相中，从而导致熔体中轻稀土相对富集。样品编号LaCePrNdSmEuGdTbDyHoErTmYbLu∑REE∑LREE∑HREE∑LREE/∑HREEδEuDC-132.5662.346.8921.454.320.653.250.563.450.782.120.322.050.30145.32128.2117.117.500.48DC-228.6556.326.2519.873.980.603.050.523.200.721.980.281.850.25130.62115.6714.957.740.46DC-325.3450.125.6717.563.560.552.800.483.000.681.800.251.700.22119.43102.8016.636.180.44DC-423.4546.785.2316.343.300.502.600.452.800.651.700.231.600.20112.5698.6513.917.100.42DC-527.8954.326.1218.903.800.582.950.503.100.701.900.261.800.23128.25111.6116.646.710.45铕异常（δEu）是衡量稀土元素地球化学特征的重要参数之一，其计算公式为δEu=EuN/（SmN×GdN）¹/²，其中EuN、SmN、GdN分别为样品中Eu、Sm、Gd元素的球粒陨石标准化值。稻城冬措花岗岩的δEu值变化范围为0.42-0.48，平均值为0.45，表现出明显的负Eu异常。负Eu异常通常与岩浆源区存在斜长石残留或岩浆演化过程中斜长石的结晶分异作用有关。斜长石是一种含钙、钠的铝硅酸盐矿物，其中Eu²⁺与Ca²⁺的离子半径相近，在斜长石结晶过程中，Eu²⁺优先进入斜长石晶格。当岩浆源区存在斜长石残留时，会导致熔体中Eu相对亏损，从而出现负Eu异常；在岩浆演化过程中，若斜长石结晶析出，也会使残余岩浆中的Eu含量降低，表现为负Eu异常。在稻城冬措花岗岩中，明显的负Eu异常表明其岩浆源区可能存在斜长石残留，或者在岩浆上升侵位过程中经历了斜长石的结晶分异作用。为了更直观地展示稻城冬措花岗岩的稀土元素配分特征，绘制了稀土元素球粒陨石标准化配分曲线，如图4所示。在配分曲线上，所有样品的曲线形态相似，均表现为轻稀土元素（La-Eu）部分曲线斜率较大，向右陡倾，表明轻稀土元素之间分馏明显；重稀土元素（Gd-Lu）部分曲线相对平缓，斜率较小，表明重稀土元素之间分馏较弱。这种配分模式与典型的壳源花岗岩稀土元素配分模式相似，进一步说明稻城冬措花岗岩的岩浆源区可能主要来自地壳物质。在壳源物质部分熔融形成岩浆的过程中，由于稀土元素的地球化学性质差异，轻稀土元素更容易进入熔体相，且在岩浆演化过程中，轻稀土元素之间的分馏作用更为明显，而重稀土元素相对较为稳定，分馏作用较弱，从而形成了这种轻稀土强烈富集、分馏明显，重稀土相对亏损、分馏较弱的配分模式。图4：稻城冬措花岗岩稀土元素球粒陨石标准化配分曲线稻城冬措花岗岩的稀土元素特征与区域内其他花岗岩相比，具有一定的相似性和差异性。与区域内一些印支期花岗岩相比，其稀土元素总量和轻稀土与重稀土的比值相近，均表现为轻稀土强烈富集的特征，这反映了它们在岩浆源区物质组成和部分熔融过程上可能具有相似性，都受到了区域构造背景和深部地质过程的影响。然而，在铕异常方面，稻城冬措花岗岩的负Eu异常程度相对较高，这可能与该岩体的岩浆源区中斜长石残留量相对较多，或者岩浆演化过程中斜长石的结晶分异作用更为强烈有关。这种差异为研究区域内不同花岗岩的形成机制和演化过程提供了重要线索，有助于进一步揭示区域地质演化的复杂性和多样性。通过对稻城冬措花岗岩稀土元素特征的分析，可以推断其岩浆源区主要为地壳物质，且源区可能存在斜长石残留，在岩浆演化过程中经历了斜长石的结晶分异作用。轻稀土强烈富集、重稀土相对亏损以及明显的负Eu异常等特征，为研究其岩石成因、岩浆演化过程以及构造背景提供了重要的地球化学依据。4.2.2其他微量元素特征除稀土元素外，稻城冬措花岗岩的其他微量元素也具有独特的地球化学特征。对[X]件样品的微量元素分析结果如表3所示。样品编号RbSrBaZrHfNbTaThUTiDC-1165.32280.45450.67150.234.2012.560.8518.653.203200DC-2158.67265.32430.56145.344.0512.000.8017.893.003100DC-3150.45250.12410.34140.253.9011.500.7517.002.803000DC-4145.32240.67390.23135.453.8011.000.7016.502.602900DC-5155.67260.45420.56142.344.1012.200.8218.002.903050从大离子亲石元素（LILE）来看，Rb含量变化范围为145.32×10⁻⁶-165.32×10⁻⁶，平均含量为155.01×10⁻⁶；Sr含量变化范围为240.67×10⁻⁶-280.45×10⁻⁶，平均含量为259.49×10⁻⁶；Ba含量变化范围为390.23×10⁻⁶-450.67×10⁻⁶，平均含量为420.45×10⁻⁶。Rb/Sr比值变化范围为0.52-0.69，平均值为0.60，相对较高。高Rb/Sr比值通常指示岩浆源区相对富集Rb元素，且可能经历了一定程度的交代作用。在岩浆源区，流体的交代作用可以使Rb元素相对富集，而Sr元素相对亏损，从而导致Rb/Sr比值升高。此外，Rb和Sr在矿物中的分配系数不同，在岩浆演化过程中，随着钾长石等含Rb矿物的结晶析出，残余岩浆中的Rb含量会相对降低，而Sr在斜长石等矿物中相对富集，进一步影响Rb/Sr比值。在稻城冬措花岗岩中，较高的Rb/Sr比值暗示其岩浆源区可能受到了流体交代作用的影响，或者在岩浆演化过程中经历了特定的矿物结晶分异过程。高场强元素（HFSE）方面，Zr含量变化范围为135.45×10⁻⁶-150.23×10⁻⁶，平均含量为142.76×10⁻⁶；Hf含量变化范围为3.80×10⁻⁶-4.20×10⁻⁶，平均含量为4.01×10⁻⁶；Nb含量变化范围为11.00×10⁻⁶-12.56×10⁻⁶，平均含量为11.85×10⁻⁶；Ta含量变化范围为0.70×10⁻⁶-0.85×10⁻⁶，平均含量为0.78×10⁻⁶。Zr/Hf比值变化范围为34.99-37.56，平均值为35.60，与地壳平均值（约36）接近，表明该花岗岩的Zr和Hf在岩浆演化过程中具有相似的地球化学行为，没有发生明显的分异。Zr和Hf具有相似的离子半径和化学性质，在大多数地质过程中倾向于一起活动，因此Zr/Hf比值相对稳定。然而，在一些特殊的地质环境中，如部分熔融程度较高或受到地幔物质混染时，Zr/Hf比值可能会发生变化。稻城冬措花岗岩的Zr/Hf比值接近地壳平均值，说明其岩浆源区主要为地壳物质，且在岩浆演化过程中没有受到明显的地幔物质混染。Nb/Ta比值变化范围为14.29-15.70，平均值为15.19，略高于地壳平均值（约11-12）。较高的Nb/Ta比值可能与岩浆源区中含Ta矿物的优先结晶或部分熔融过程中Ta元素的相对亏损有关。在岩浆源区部分熔融时，Ta元素可能更容易进入残留相，导致熔体中Nb/Ta比值升高；或者在岩浆演化过程中，含Ta矿物（如钽铁矿等）较早结晶析出，使得残余岩浆中的Nb相对富集，从而提高了Nb/Ta比值。这种较高的Nb/Ta比值反映了稻城冬措花岗岩在形成过程中，其岩浆源区物质的部分熔融和岩浆演化过程具有一定的特殊性。Th含量变化范围为16.50×10⁻⁶-18.65×10⁻⁶，平均含量为17.61×10⁻⁶；U含量变化范围为2.60×10⁻⁶-3.20×10⁻⁶，平均含量为2.90×10⁻⁶，Th/U比值变化范围为5.16-6.17，平均值为6.07，与地壳平均值（约4-5）相比略高。较高的Th/U比值通常与地壳物质的部分熔融有关，因为Th和U在地球化学性质上存在一定差异，在部分熔融过程中，Th更容易进入熔体相，而U相对更倾向于保留在残留相中，导致熔体中Th/U比值升高。在稻城冬措花岗岩中，较高的Th/U比值进一步支持了其岩浆源区主要为地壳物质的观点，且表明源区物质在部分熔融过程中，Th和U发生了一定程度的分异。为了深入分析稻城冬措花岗岩微量元素对岩石成因和构造环境的指示作用，绘制了微量元素原始地幔标准化蛛网图，如图5所示。在蛛网图中，以原始地幔标准化值为基准，将样品中微量元素的含量与之对比。可以看出，稻城冬措花岗岩相对原始地幔，富集大离子亲石元素（如Rb、Ba、Th等），亏损高场强元素（如Nb、Ta、Ti等）。富集大离子亲石元素可能与岩浆源区的地壳物质组成以及流体交代作用有关，地壳物质中通常富含大离子亲石元素，而流体交代作用可以进一步富集这些元素。亏损高场强元素（如Nb、Ta、Ti等）是典型的岛弧岩浆岩特征。在岛弧环境中，俯冲板片释放的流体交代地幔楔，导致地幔楔部分熔融形成的岩浆中，高场强元素相对亏损。因为高场强元素在流体中的溶解度较低，不易被流体携带进入地幔楔，从而在岩浆中表现为亏损状态。稻城冬措花岗岩的这种微量元素特征表明，其形成可能与古特提斯洋的俯冲作用有关，在俯冲过程中，俯冲板片释放的流体交代地幔楔，使得地幔楔物质部分熔融，形成的岩浆上升与地壳物质混合，最终形成了具有岛弧岩浆岩特征的稻城冬措花岗岩。图5：稻城冬措花岗岩微量元素原始地幔标准化蛛网图通过对稻城冬措花岗岩其他微量元素特征的分析，进一步揭示了其岩浆源区性质、岩浆演化过程以及形成的构造环境。较高的Rb/Sr比值、接近地壳平均值的Zr/Hf比值、略高的Nb/Ta比值和Th/U比值，以及微量元素原始地幔标准化蛛网图中显示的富集大离子亲石元素、亏损高场强元素的特征，为深入理解其岩石成因和构造背景提供了重要的地球化学信息。这些信息与主量元素和稀土元素分析结果相互印证，共同揭示了稻城冬措花岗岩在区域地质演化中的重要意义。4.3地球化学特征的地质意义稻城冬措花岗岩的主量元素、微量元素及稀土元素地球化学特征蕴含着丰富的地质信息，对揭示其成因、源区及构造背景具有重要意义。从成因角度来看，主量元素特征表明稻城冬措花岗岩为过铝质钙碱性系列花岗岩。过铝质特征通常与源区物质中富铝矿物的存在或岩浆演化过程中的交代作用有关。在源区物质部分熔融过程中，若存在大量的富铝矿物，如白云母、堇青石等，会使岩浆中铝含量相对较高；岩浆在上升侵位过程中与围岩发生交代作用，也可能导致铝元素富集。钙碱性系列特征指示其形成于板块俯冲带或碰撞造山带等构造环境，在这种环境下，俯冲板片释放的流体交代地幔楔，导致地幔楔部分熔融，形成的岩浆上升与地壳物质混合，最终形成钙碱性系列花岗岩。结合区域地质背景，该花岗岩可能是在古特提斯洋闭合过程中，扬子板块与羌塘-三江板块碰撞拼合，地壳物质发生强烈挤压和部分熔融的产物。微量元素特征进一步支持了上述成因推断。稀土元素配分模式显示轻稀土强烈富集、重稀土相对亏损以及明显的负Eu异常。轻稀土强烈富集、重稀土相对亏损表明岩浆源区可能主要来自地壳物质，在源区物质部分熔融时，轻稀土元素更容易进入熔体相，而重稀土元素则相对倾向于保留在残留相中。明显的负Eu异常通常与岩浆源区存在斜长石残留或岩浆演化过程中斜长石的结晶分异作用有关。这说明稻城冬措花岗岩的岩浆在形成和演化过程中，经历了斜长石的分离结晶或源区存在斜长石残留，进一步证实了其岩浆起源于地壳物质部分熔融的观点。其他微量元素中，大离子亲石元素（如Rb、Ba、Th等）相对富集，高场强元素（如Nb、Ta、Ti等）相对亏损，这是典型的岛弧岩浆岩特征。在岛弧环境中，俯冲板片释放的流体交代地幔楔，导致地幔楔部分熔融形成的岩浆中，高场强元素相对亏损，因为高场强元素在流体中的溶解度较低，不易被流体携带进入地幔楔。而大离子亲石元素在流体中的活动性较强，容易被富集。稻城冬措花岗岩的这种微量元素特征表明其形成可能与古特提斯洋的俯冲作用密切相关，在俯冲过程中，俯冲板片释放的流体交代地幔楔，使得地幔楔物质部分熔融，形成的岩浆上升与地壳物质混合，最终形成了具有岛弧岩浆岩特征的稻城冬措花岗岩。关于源区，主量元素中的高SiO₂含量、过铝质特征以及相对较低的FeO/MgO值，表明岩浆源区可能以长英质物质为主，且具有相对较低的氧化程度。稀土元素的配分模式和铕异常特征进一步支持了源区主要为地壳物质的观点，且源区可能存在斜长石残留。其他微量元素中，Zr/Hf比值接近地壳平均值，说明其岩浆源区主要为地壳物质，且在岩浆演化过程中没有受到明显的地幔物质混染。较高的Th/U比值也与地壳物质的部分熔融有关，因为在部分熔融过程中，Th更容易进入熔体相，而U相对更倾向于保留在残留相中，导致熔体中Th/U比值升高。这些特征综合表明，稻城冬措花岗岩的岩浆源区主要为地壳物质，且可能经历了一定程度的交代作用，源区中存在斜长石等矿物残留。在构造背景方面，稻城冬措花岗岩的地球化学特征指示其形成于板块碰撞造山带的构造环境。主量元素中的钙碱性系列特征、微量元素中富集大离子亲石元素、亏损高场强元素的特征，以及稀土元素配分模式等，都与板块碰撞造山带环境下形成的花岗岩特征相符。结合区域地质背景，该花岗岩形成于印支期，这一时期正是古特提斯洋闭合，扬子板块与羌塘-三江板块强烈碰撞拼合的时期。在碰撞过程中，地壳物质受到强烈挤压，发生部分熔融，形成的岩浆沿断裂构造上升侵位，形成了稻城冬措花岗岩体。其形成过程受到了区域构造应力场的控制，岩体的展布方向与区域构造线方向一致，也进一步证明了其与区域构造运动的密切关系。稻城冬措花岗岩的地球化学特征是其成因、源区及构造背景的综合反映。通过对这些特征的深入分析，能够为研究区域地质演化历史提供重要的地球化学依据，对于理解青藏高原东南缘的构造演化过程具有重要意义。五、锆石SHRIMPU-Pb定年5.1定年原理与方法锆石SHRIMPU-Pb定年基于放射性同位素衰变原理。锆石（ZrSiO₄）是一种在岩浆岩、变质岩和沉积岩中广泛存在的副矿物，具有极为稳定的晶体结构。在锆石形成过程中，U（铀）和Th（钍）等放射性元素能够以类质同象的方式进入其晶格，而初始状态下锆石中几乎不含有普通Pb（铅）。U存在两种主要的放射性同位素，即²³⁸U和²³⁵U，它们会分别通过一系列的α衰变和β衰变，最终衰变成稳定的²⁰⁶Pb和²⁰⁷Pb。²³⁸U的半衰期为4.47Ga（1Ga=10亿年），²³⁵U的半衰期为0.70Ga。根据放射性衰变定律，单位时间内衰变的原子数与现存放射性母体的原子数成正比。通过精确测量锆石中²³⁸U、²³⁵U以及它们衰变产生的²⁰⁶Pb、²⁰⁷Pb的含量，利用相应的衰变公式，就可以计算出锆石从形成至今所经历的时间，即获得其形成年龄。在实际操作中，首先进行样品的前期处理。将采集的稻城冬措花岗岩样品粉碎至合适粒度，通过重液分离和磁选等常规方法初步分离出锆石矿物。之后，在双目镜下精心挑选晶型完好、透明度高且无明显裂纹和包裹体的锆石颗粒。这些特征良好的锆石颗粒能够保证后续分析数据的可靠性，因为晶型完好和透明度高的锆石通常具有较为均一的内部结构，较少受到后期地质作用的干扰，而无明显裂纹和包裹体可以避免在分析过程中因外来物质的混入或内部结构的破坏导致数据偏差。挑选好的锆石颗粒与已知年龄的标准锆石（如Temuora等）一同粘贴在环氧树脂表面，制成样品靶。标准锆石的年龄经过精确测定，在分析过程中用于校准仪器，以提高分析结果的准确性。对样品靶进行抛光处理，使锆石表面平整光滑，为后续的分析提供良好的条件。利用阴极发光（CL）成像技术对锆石进行内部结构分析是至关重要的一步。CL成像原理是基于电子束激发下，锆石内部不同区域的晶体结构和微量元素分布差异会导致其发出不同强度和颜色的光。通过CL成像，可以清晰观察到锆石的生长环带、内部结构特征以及是否存在继承核等信息。生长环带是锆石在结晶生长过程中，由于物理化学条件的周期性变化而形成的，不同的生长环带对应着不同的结晶阶段，其宽度、形状和亮度等特征能够反映出结晶环境的变化。继承核则是锆石在形成过程中捕获的早期锆石颗粒，其年龄代表了更早的地质事件。通过对这些信息的分析，可以了解锆石的生长历史和来源，为选择合适的定年分析点位提供科学依据。例如，对于具有明显生长环带的锆石，选择不同环带进行分析，可以获取岩浆在不同结晶阶段的年龄信息，从而更全面地了解岩浆的演化过程。锆石U-Pb同位素分析在高灵敏度高分辨率离子探针（SHRIMP）上进行。在分析过程中，利用一次离子束（通常为O⁻或Cs⁺）以高能量轰击锆石表面。一次离子束的能量足以使锆石中的U、Pb等元素离子化，产生二次离子。这些二次离子在电场的作用下被引出，并进入质量分析器。质量分析器根据二次离子的质荷比（m/z）对其进行分离和检测，通过精确测量不同质荷比的二次离子强度，就可以确定锆石中²³⁸U、²³⁵U、²⁰⁶Pb、²⁰⁷Pb等同位素的含量比值。为了确保分析结果的准确性和可靠性，每个锆石颗粒会选择多个分析点位进行测定。这是因为锆石内部不同部位的同位素组成可能存在微小差异，多点位测定可以更全面地反映锆石的年龄信息，减少因局部差异导致的年龄误差。在分析过程中，还会实时监测仪器的稳定性和分析数据的质量，对异常数据进行剔除和重新分析。数据处理和年龄计算采用专业的软件，如Isoplot等。利用所测得的²³⁸U/²⁰⁶Pb、²³⁵U/²⁰⁷Pb等同位素比值，在Isoplot软件中绘制U-Pb谐和图。U-Pb谐和图以²⁰⁶Pb/²³⁸U年龄为横坐标，²⁰⁷Pb/²³⁵U年龄为纵坐标。在理想情况下，未受到后期地质作用扰动的锆石数据点应落在谐和线上，表明²³⁸U-²⁰⁶Pb和²³⁵U-²⁰⁷Pb两个衰变体系所计算出的年龄一致。然而，在实际情况中，由于锆石可能受到后期热事件、变质作用或放射性成因Pb的丢失等因素的影响，数据点可能会偏离谐和线。对于偏离谐和线的数据，需要进行仔细分析，判断其偏离原因。例如，如果数据点位于谐和线下方，可能是由于放射性成因Pb的丢失导致；如果位于谐和线上方，则可能受到了后期U的加入或其他复杂地质作用的影响。通过对数据点的分布和偏离情况进行分析，可以筛选出可靠的数据，采用加权平均法计算锆石的U-Pb年龄。加权平均法考虑了每个数据点的误差大小，误差较小的数据点在计算中具有更大的权重，从而能够得到更精确的年龄值。最终确定稻城冬措花岗岩的形成时代，并结合区域地质背景，探讨其与区域构造演化事件的关系。5.2样品采集与测试本次研究的样品采集工作在稻城冬措花岗岩体出露区域精心展开，充分考虑了岩体的空间分布、岩石类型变化以及构造部位差异等因素，以确保采集的样品具有广泛代表性。在岩体的不同位置，按照一定的间距进行系统采样，共采集了[X]件新鲜的花岗岩样品。采样点均匀分布于岩体的中心部位、边缘地带以及不同岩性变化区域，涵盖了中粒黑云二长花岗岩等主要岩石类型。例如，在岩体的东部边缘，选取了具有明显接触变质特征的样品，以研究岩浆与围岩的相互作用对岩石性质的影响；在岩体的中心区域，采集了结晶程度良好、矿物组成均匀的样品，用于分析岩浆的原始特征。在采样过程中，严格遵循科学规范的操作流程。首先，使用地质锤仔细敲取岩石样品，确保样品表面新鲜，无风化、蚀变等现象。对于每一件样品，详细记录其地理坐标、出露位置、岩石产状以及与围岩的接触关系等信息。地理坐标通过高精度的GPS定位仪进行测量，精确到小数点后四位，以保证样品位置的准确性；岩石产状包括走向、倾向和倾角，使用地质罗盘进行测量，确保数据的可靠性。同时，对样品的宏观特征进行详细描述，如颜色、结构、构造、矿物组成等，为后续的室内分析提供全面的野外资料。采集的样品用塑料薄膜包裹，再装入专用的样品袋中，并贴上标签，注明样品编号、采样地点、采样日期等信息，以防止样品混淆和损坏。将采集的样品送往专业实验室进行锆石分选。在实验室中，首先将岩石样品破碎至合适粒度，一般控制在80-120目之间。通过重液分离方法，利用锆石与其他矿物密度的差异，初步分离出锆石矿物。重液通常选用三溴甲烷等密度较大的液体，将破碎后的岩石颗粒放入重液中，锆石由于密度较大而沉降到容器底部，其他低密度矿物则漂浮在重液表面，从而实现初步分离。之后，采用磁选方法进一步提纯锆石。由于锆石不具有磁性，而部分其他矿物（如磁铁矿等）具有磁性，在磁场作用下，磁性矿物会被吸引到磁选设备上，而锆石则得以保留，从而提高了锆石的纯度。最后，在双目镜下对分离出的锆石颗粒进行人工挑选。挑选出晶型完好、透明度高、无明显裂纹和包裹体的锆石颗粒。晶型完好的锆石通常具有规则的几何外形，如柱状、双锥状等，这表明其在结晶过程中生长环境相对稳定；透明度高说明锆石内部结构较为纯净，杂质较少；无明显裂纹和包裹体可以避免在分析过程中因外来物质的混入或内部结构的破坏导致数据偏差。经过仔细挑选，最终获得了[X]颗符合要求的锆石颗粒，用于后续的SHRIMPU-Pb定年分析。在锆石SHRIMPU-Pb定年测试过程中，采取了一系列严格的质量控制措施。首先，对分析仪器SHRIMP进行定期校准和维护。在每次分析前，使用标准锆石（如Temuora等，其U-Pb年龄经过精确测定且具有良好的稳定性）对仪器进行校准，确保仪器的分析精度和准确性。校准过程中，对标准锆石进行多次测量，将测量结果与已知年龄进行对比，根据偏差对仪器参数进行调整，使仪器处于最佳工作状态。在分析过程中，实时监测仪器的运行状态，包括离子束的稳定性、质量分析器的分辨率等参数。一旦发现仪器出现异常，立即停止分析，进行故障排查和修复，确保分析数据的可靠性。对于每一个锆石颗粒的分析，都进行多次测量。每个锆石颗粒选择3-5个不同的分析点位进行测定，以减少因锆石内部同位素分布不均匀导致的误差。在测量过程中，严格控制分析条件的一致性，包括一次离子束的能量、束流强度、轰击时间等参数。例如，一次离子束能量控制在10-15keV之间，束流强度保持在1-3nA，轰击时间根据锆石的性质和分析要求调整在30-60s之间。对每个分析点位的测量数据进行实时监控，剔除异常数据。异常数据可能由于分析过程中的干扰、仪器波动或锆石内部结构不均匀等原因产生，如测量得到的同位素比值明显偏离其他点位的数据，或者测量误差过大等情况。对于异常数据，重新选择分析点位进行测量，确保每个锆石颗粒获得的数据具有代表性和可靠性。在数据处理阶段，