青藏高原Soumi-NPP和MODIS积雪范围产品的对比分析

Soumi-NPP（Soumi Polar-orbiting Partnership）卫星作为接替服役超期的Terra、Aqua卫星,其积雪范围产品在青藏高原的精度尚未被评价。以Soumi-NPP积雪范围产品为研究对象,利用气象台站点数据并结合更高分辨率的Landsat-8 OLI数据,评价该产品的精度,并与MODIS（Moderate Resolution Imaging Spectroradiometer）积雪范围产品进行对比分析。结果表明：使用气象台站进行数据验证时,NPP、MOD与MYD三种积雪范围产品的总精度均较高,但三者积雪漏分误差都较大,其中MYD的漏分误差最大,为64.2%;当雪深小于5 cm时,三种积雪范围产品的积雪分类精度都较低,雪深大于等于5 cm时,NPP积雪范围产品的积雪分类精度最高,为82.3%,MOD与MYD的精度分别为77.1%和69.4%;利用Landsat-8 OLI数据验证时,Soumi-NPP积雪范围产品的Kappa系数最高,其均值为0.707,为高度一致性。而MOD10A1与MYD10A1的Kappa系数较低,分别为0.476与0.557,为中等一致性;Soumi-NPP积雪范围产品的Kappa系数大多在0.6以上,精度比较稳定,而MODIS积雪范围产品的Kappa系数波动较大,精度稳定性较差。Soumi-NPP积雪范围产品相较于MODIS积雪范围产品,其精度有了较大的提升,为准确监测青藏高原积雪范围提供了一个更优的选择。     

NOAA IMS雪冰产品在青藏高原积雪监测中的适用性分析

在系统评估青藏高原积雪观测典型气象站历史定位坐标精度基础上,利用站点雪深资料对NOAA IMS 4 km和1 km分辨率雪冰产品在青藏高原的精度和适用性进行了验证和评估,定量分析了IMS 4 km到1 km空间分辨率提高和气象站历史定位与GPS定位坐标之间的差异对青藏高原IMS积雪监测精度的影响。结果表明：青藏高原个别气象站历史坐标与当前GPS接收机定位之间存在较大的差异,如安多气象站经度偏小0.6°,纬度偏大0.08°。IMS 4 km雪冰产品在青藏高原的总精度介于76.4%~83.2%,平均为80.1%,积雪分类精度介于35.8%~60.7%,平均为47.2%,平均误判率为17.1%,平均漏判率为45.5%,总体上呈现地面观测的积雪日数越多、平均雪深越大,其总体监测精度越低,而积雪分类精度越高的特点。IMS分辨率从4 km到1 km总体精度平均提高了2.9%,积雪分类精度平均提高了0.9%,主要是由于个别站点的精度提升较大引起的,对高原多数台站积雪监测精度的改进和提升很小。除个别台站外,目前气象站历史坐标和GPS定位坐标之间的差异,对IMS 4 km积雪监测精度验证结果没有影响。然而,今后随着卫星遥感技术的发展,更高时空分辨率的遥感积雪产品将用于积雪监测和研究,精确的地面观测站坐标信息是对这些遥感数据开展精度验证与实际应用的前提。     

青藏高原不同土地覆盖类型下积雪面积判别算法优化北大核心CSCD

MODIS V006版本数据仅提供了归一化积雪指数(NDSI),而用户往往关心的是直观的积雪分类,包括积雪范围或积雪覆盖率。美国国家冰雪数据中心推荐全球积雪范围最佳的NDSI阈值为0.4,但是青藏高原地形复杂多样,积雪斑块化特征明显,单一阈值并不能精确地判识不同下垫面上的积雪。青藏高原被称为地球的第三极,是中国三大稳定积雪区之一,蕴藏了大量的淡水资源。随着全球气候变暖,青藏高原地区积雪融化时间提前,冰川融水增加,影响河流水量,造成洪涝灾害,进而影响人类正常生产生活,因此通过确定不同下垫面阈值,改善传统阈值的积雪高估低估现象,提高积雪识别精度,进而更准确地探究青藏高原积雪状况,显得尤为迫切。本文以青藏高原为研究对象,首先生成MODIS逐日无云NDSI序列并进行验证;其次对应站点雪深数据与NDSI序列,证实在下垫面为林地和非林地的区域,去云NDSI序列与站点雪深均有良好的对应关系,确定不同下垫面最优阈值范围;最后在最优阈值范围内通过混淆矩阵确定最优阈值。计算得出,林地NDSI=0.03时,总体精度最高为94.02%,在该NDSI之下,高估误差OE和低估误差UE分别为1.21%和4.60%;非林地NDSI=0.26时,总体精度OA最高为94.27%,在该NDSI之下,高估误差OE和低估误差UE分别为0.51%和5.03%。因此选取优化后林地阈值为NDSI=0.03,非林地阈值为NDSI=0.26。为避免地面常规观测资料尺度上的局限性,本文采用高精度的Landsat 8 OLI卫星数据识别结果,作为“真值”对优化后阈值的判别结果进行“像元—像元”级别的验证。在定量验证中,优化后NDSI阈值对MOD10A1 V006积雪判别结果的总体精度OA为84.21%,高估误差OE为5.33%,低估误差UE为10.46%;传统阈值对MOD10A1 V006积雪判别结果的总体精度OA为82.86%,高估误差OE为1.48%,低估误差UE为15.66%。可以看出在定量验证中,优化后阈值的积雪判别精度更高。同时在定性验证中,积雪大面积集中的区域,新的阈值与传统阈值提取效果均相对较好;积雪相对分散破碎的区域,优化后阈值能提取出大量积雪,传统阈值则不能。这表明考虑不同土地覆盖类型下的NDSI阈值优化可以有效地提高青藏高原积雪判别精度,为NDSI在积雪识别中的应用提供有力的支撑,有助于更准确地了解该地区积雪分布状况。     

亚洲中部干旱区积雪时空变异遥感分析

亚洲中部干旱区的大尺度遥感积雪信息研究,可在跨界河流水资源分配利用方面提供数据支持,对国家重大战略的生态安全保障有重要作用。采用数据融合方法,将MOD10A2和MYD10A2数据进行融合去云处理,结合气象站点积雪数据评估去云后的积雪识别精度;提取积雪覆盖率(SCP)与积雪日数(SCD)信息,分析SCP与SCD年际、年内变化差异;结合数字高程模型,分析不同高程带下SCP的时空变化规律。结果表明:(1)MOD10A2与MYD10A2融合去云处理,可有效去除云的干扰,准确提取亚洲中部干旱区积雪变化信息。(2)年内SCP最大值范围为55.7%~77.4%,最小值范围为1.6%~2.9%,融雪期SCP下降速率具有明显地域差异,总体SCP呈缓慢增加趋势。(3)总体SCD呈略微下降趋势,32.2%的区域呈下降趋势,30.9%的区域呈增加趋势,36.9%的区域保持稳定不变。(4)海拔1 000 m以下,SCP年内随季节变化呈U型,年际变化显著;1 000~4 000 m区域,SCP年内均随季节的变化呈现出V型,年际变化呈现出稳定性波动;6 000 m以上为永久性积雪,季节、时空变化差异性均不明显。     

基于MODIS的青藏高原季节性积雪去云方法可行性比较研究

青藏高原地处中纬度地区,季节性积雪分布破碎,地面观测站点稀少,中分辨率成像光谱仪（MODIS）可为该地区提供每日积雪监测数据,然而云是光学遥感最大的影响因素,为研究MODIS每日积雪产品去云方法在青藏高原的适用性,根据原理将去云方法归纳总结为五大类,并对每种方法的“潜在假设”开展分析讨论。结果显示：基于时间连续性的方法适用性强,去云效果明显,上下午积雪连续的平均概率为72.5%,而2~5d的连续积雪的概率为5.6%~43%不等,可靠性差;临近像元法可去除零散分布的云,平均正确率达到95.5%,但去除云量较少;基于高程的去云算法在山区适用性好,而在高原腹地由于坡度较小而错判概率较大;采用被动微波遥感数据进行去云则依赖于微波对云的识别率,往往误差较大;采用数学方法拟合积雪边界在积雪破碎、降雪融雪较快的青藏高原地区,物理意义较弱。通过分析研究表明,青藏高原地区MODIS日积雪产品的去云,需综合多种算法的区域适用性,充分考虑青藏高原地形及积雪本身的特征,逐步完善每日积雪去云工作。     

基于MODIS的青藏高原逐日无云积雪产品算法

青藏高原积雪对高亚洲地区水和能量循环起着重要的反馈和调节作用,其变化影响着融雪性河流流量,对下游水资源和经济活动具有重要影响。中分辨率成像光谱仪（MODIS）具有较高的时空分辨率,被广泛应用于积雪遥感动态监测,然而光学遥感积雪受云层影响严重,且青藏高原地区水汽分布不均,局地对流活跃,积雪的赋存时间变化快,这给高原地区逐日积雪监测及其气候学制图带来挑战。在考虑青藏高原地形和积雪分布特征情况下,结合现有的云覆盖下积雪判别算法,采用8个不同方法的组合,逐步实现MODIS逐日无云积雪算法。选取2009年10月1日-2011年4月30日两个积雪季为研究期,并采用145个地面台站观测雪深数据对去云算法各步骤过程开展精度验证,结果表明：当积雪深度>3 cm时,逐日无云积雪产品总分类精度达到96.6%,积雪分类精度达83%,积雪判对概率（召回率）达到89.0%,算法可实现青藏高原地区逐日无云积雪动态监测和积雪覆盖气候学数据重建,对高亚洲地区的水、生态和灾害等全球环境变化影响研究具有重要的意义。     

东北地区MODIS和AMSR-E积雪产品验证及对比

通过2002-2008年6个积雪季节的Terra-Aqua/MODIS积雪产品(MOD10A2、MOD10C2)和Aqua/AMSR-E雪水当量产品,分析了东北地区积雪覆盖面积的变化特征,以研究区气象站点观测的积雪数据为真实值来验证两种产品积雪信息的精度,探讨了云覆盖、土地利用类型和雪深对积雪覆盖精度的影响.结果表明:云的存在对微波数据积雪识别的影响较小,在积雪量较多的12月至次年的2月随云量百分比的变化,MOD10A2积雪覆盖面积比例大体出现负变化.因此,在有云情况下AMSR-E数据反演积雪精度最好.对比草地、耕地、林地和居民地4种土地覆盖类型对监测积雪覆盖精度的影响,发现林地对其影响最大,在林区3种积雪产品的积雪识别精度分别为55.8%、81.2%、85.4%;雪深对AMSR-E积雪产品识别精度影响较小,总体精度为97.8%;积雪深度对MOD10A2积雪产品识别精度影响较大,总体精度为57.3%.MOD10A2、MOD10C2和AMSR-E 3种积雪产品的总体反演精度分别为69.3%、76.6%、76.3%.有必要开发适用于东北地区的积雪覆盖算法,提高估算精度,为能量平衡估算、气候模型、农业生产、土壤墒情监测服务.     

基于HJ-1B卫星数据的积雪面积制图算法研究

积雪是影响气候变化的重要因子, 采用更高时空分辨率的环境减灾卫星遥感数据进行积雪制图算法的研究, 对推进我国自主遥感卫星在积雪监测领域的应用具有重要意义. 采用环境减灾HJ-1B卫星数据, 以青海省果洛藏族自治州达日县为研究区, 应用归一化差值积雪指数(NDSI)法建立了基于HJ-1B卫星数据的积雪面积制图算法, 并比较MODIS与HJ-1B积雪图精度. 结果表明: 研究区HJ-1B积雪制图合理的NDSI阈值为0.37, 总分类精度达到97.97%; 与"真值"影像比较, HJ-1B积雪图K_hat系数为0.911, 高于MODIS的0.817. 说明该研究建立的基于HJ-1B积雪制图算法精度可靠, 适合对研究区积雪进行实时动态监测. HJ-1B更高的空间分辨率对提高研究区积雪覆盖面积监测精度具有重要的使用价值, 但是地形因素是影响HJ-1B数据积雪分类精度的一个重要原因, 随着坡度的增加, 分类误差也随之增大, 尤其是多测误差增加比较显著.     

基于AVHRR影像的北半球积雪识别算法

针对2000年前北半球较高时空分辨率和高精度的历史积雪范围数据缺失问题,利用NOAA-AVHRR地表反射率数据,以Landsat-5 TM生成的积雪范围影像作为参考真值,优化基于多指标的多级决策树积雪识别算法的阈值,并结合云雪混淆区分技术,生成了北半球AVHRR 1981—1999年L1级逐日积雪范围数据集。此外,针对AVHRR在高纬度地区数据完全缺失和低纬度地区数据部分缺失问题,利用微波雪深数据集进行填充,生成了北半球L2级逐日积雪范围数据集。最后,利用北半球1981—1999年间2 546个气象台站记录的雪深数据和939景Landsat-5 TM参考积雪范围影像作为验证数据,对AVHRR积雪范围数据集进行了精度验证。结果表明：L1级和L2级数据集的总体精度分别为81.8%和82.2%,用户精度分别为83.7%和83.8%,生产者精度分别为81.7%和84.2%,说明算法精度较高,错分误差和漏分误差均比较均衡。进一步利用Landsat-5 TM参考积雪范围影像对L2级数据集进行面上精度评估,发现L2级数据集的总体精度为90.3%,用户精度为90.2%,生产者精度为99.1%,L2级数据集精度较高。生成的北半球历史数据集可为全球积雪变化研究提供有效数据补充。     

MODIS逐日积雪产品去云算法研究

由于积雪和云的反射特性, 使用光学遥感监测积雪受到天气的严重干扰, 对研究区云量的分析表明, 无论是MOD10A1还是MYD10A1, 云都是影响该产品对研究区积雪进行实时监测的最大影响因素. 综合不同去云方法, 利用MODIS逐日积雪产品和被动微波数据AMSR-E雪水当量产品, 生成了MODIS逐日无云积雪图像, 并利用研究区85个地面气象观测台站提供的雪深数据对合成的单日无云积雪产品进行验证. 结果表明: 当积雪深度>3 cm时, 新产品的积雪分类精度达到91.7%, 该产品对实时监测青藏高原积雪动态变化具有重要的使用价值.     

祁连山区MODIS积雪反照率产品的精度验证及云下积雪反照率估算研究

积雪反照率在全球气候和能量收支平衡模型中起着重要的作用. 利用祁连山地区大冬树垭口站点反照率实测数据对由TM/ETM+得到的反照率数据进行标定, 然后将TM/ETM+反照率数据通过升尺度对MODIS逐日积雪反照率(SAD)产品在晴空条件下的精度进行了验证. 同时, 发展了一个基于MODIS SAD与AMSR-E SWE数据融合并结合Noah积雪反照率参数化方案估算MODIS SAD数据云下积雪反照率的算法, 通过统计分析纠正了云对积雪反照率的影响, 对云下积雪反照率进行了验证分析. 结果表明:MODIS SAD产品在祁连山地区的精度要低于大面积积雪覆盖的平坦地区(如格陵兰岛), 其平均绝对误差及均方根误差分别为0.0548和0.0727; 云下积雪反照率估算方法可以有效地获取云覆盖下积雪像元的反照率值, 纠正后的无云MODIS SAD数据与地面观测值有较好的一致性, 其平均绝对误差为0.078.     

南岭西段金鸡岭复式花岗岩基地质及岩浆动力

对南岭西段金鸡岭复式花岩岩基进行的Rb-Sr同位素定年研究，确定Ⅰ阶段黑云母二长花岗岩的等时线年龄为169．5MaⅡ、Ⅲ阶段黑云母花岗岩和二云母花岗岩为150．7Ma，其ISr值分别为0．7163和0．7206，表明该复式岩基属燕山早期。金鸡岭复式花岗岩基岩石属钙碱性系列，从Ⅰ阶段到Ⅱ、Ⅲ阶段表现出有规律的演化特征：岩石由准铝质（Ⅰ），演化为过铝质（Ⅱ）和强过铝质（Ⅲ）；钾长石有序度由0．33（Ⅰ），增高到0．69（Ⅱ）和0．80（Ⅲ）；斜长石油An37（Ⅰ）降低到An26（Ⅱ）和An7(Ⅲ)；黑云母由铁质黑云母（Ⅰ）向铁叶云母（Ⅱ）和铁白云母（Ⅲ）方向演化；氧化物（SiO2,Al2O3,TFeO,MgO,TiO2,CaO)-DI图解上呈良好的线性演化关系；成岩温度逐阶段降低，由745℃（Ⅰ）降低到673℃（Ⅱ）至505℃（Ⅲ）。采用地质地球化学方法估算出金鸡岭岩基的侵位深度约为6．3km，成岩压力为180MPa，具中深成相特征，属S型花岗岩并形成于华南板块内部的碰撞构造环境。     

藏北改则--班戈地区下白垩统多巴组的高频旋回沉积特征与盆地性质

西藏北部改则——班戈地区的下白垩统多巴组中所发育的Ⅴ级旋回表现出3种垂向沉积序列:潮坪相碎屑岩沉积序列(Ⅰ)、浅海相碳酸盐岩沉积序列(Ⅲ)和碳酸盐岩、陆源碎屑岩混合沉积序列(Ⅱ)。其中序列Ⅰ-1和Ⅲ是正常的沉积水体向上变浅的沉积序列,而序列Ⅱ具有明显的陆源碎屑先进积后退积的特征,可进一步分成3个亚类:Ⅱ-1、Ⅱ-2和Ⅱ-3。横向上从南部物源区向北呈Ⅰ→Ⅱ-1→Ⅱ-2→Ⅱ-3→Ⅲ的变化过程。结合盆地构造背景和沉积特征,本文提出多巴组高频旋回的形成与弧背盆地南部冲断带的冲断作用相对应。冲断作用期,物源区粗碎屑物质大量进积,大于盆地基底沉降速率,沉积水体逐渐变浅;冲断间歇期,物源区碎屑物质的供应迅速减少,基底沉降速率逐渐大于碎屑物质的进积速率,沉积水体又有所加深。     

东秦岭秋树湾铜钼矿流体包裹体和稳定同位素特征及其地质意义

秋树湾铜钼矿是东秦岭钼矿带上典型的受斑岩体控制的矽卡岩-斑岩角砾岩筒复合型矿床,矿体赋存于成矿母岩花岗岩及矽卡岩和角砾岩筒中。根据矿物共生组合、矿石组构、围岩蚀变及脉体的穿插关系,可划分为早(Ⅰ)、中(Ⅱ)、晚(Ⅲ)3个矿化期,再将Ⅰ期细分为干矽卡岩-钾长石化-石英阶段(Ⅰ₁)、爆破角砾岩阶段(Ⅰ₂)、湿矽卡岩阶段(Ⅰ₃)、磁铁矿阶段(Ⅰ₄);Ⅱ期分为斑岩型铜(钼)矿阶段(Ⅱ_b)和石英硫化物阶段(Ⅱ_s);Ⅲ期为方解石、重晶石、石英阶段(Ⅲ)。流体包裹体可划分为S型含子矿物多相包裹体、L型纯液相包裹体、C型含CO₂三相包裹体、W型气液两相包裹体、G型纯气相包裹体5种类型。按时间先后顺序,成矿流体的温度、盐度、氧化还原环境具有规律性的演化特征。均一温度范围: Ⅰ期为222~406℃,Ⅱ期为152~315℃,Ⅲ期为119~189℃;盐度w(NaCl_eq): Ⅰ期介于4.2%~36.5%,Ⅱ期为3.3%~34.8%,Ⅲ期为4.2%~11.9%。激光拉曼光谱及群体包裹体成分分析结果表明,第Ⅰ期流体以H₂O、CO₂、CH₄、H₂S为主,表现为还原环境;第Ⅱ期流体以H₂O、CO₂、N₂、O₂、SO₄^2-、Cl^-、F^-为主,为氧化环境,暗示流体源于岩浆。流体包裹体岩相学及包裹体测温表明,流体由早期的高温、高盐度、含CO₂的H₂O-NaCl-CO₂体系的岩浆流体在成矿Ⅰ期发生沸腾作用和相分离,伴随着流体沸腾、CO₂逸失、温度下降、大气水的加入、盐度下降等过程,导致大量金属硫化物沉淀。在成矿Ⅱ、Ⅲ期成矿体系趋于开放,流体存在大气降水混入,逐渐演化为晚期的低盐度、中低温度、贫CO₂的流体体系。H、O、S同位素结果表明有地幔流体参与成矿作用。     

云南羊拉铜矿床成矿流体特征及成矿过程

云南羊拉铜矿床位于金沙江构造带中部,是中-晚三叠世金沙江洋盆向西俯冲闭合碰撞造山过程中形成的一个大型铜矿床。矿体多呈层状、似层状顺层产出,但明显受层间破碎带和滑脱带控制。从流体包裹体研究入手,讨论了该矿床成矿流体的特征、演化以及流体不混溶(沸腾)作用与成矿的关系。流体包裹体研究表明,干夕卡岩阶段(Ⅰ)、湿夕卡岩磁铁矿阶段(Ⅱ)、石英硫化物阶段(Ⅲ)以及方解石硫化物阶段(Ⅳ)中发育多种类型的包裹体,主要为气液水两相包裹体和含子矿物多相包裹体,纯液相水包裹体次之,少见纯气相有机质包裹体。其中,含子矿物多相包裹体发育于Ⅰ阶段石榴石、Ⅱ阶段绿帘石,尤其是Ⅲ阶段石英中。Ⅰ、Ⅱ阶段成矿流体具有高温、高盐度特征,均一温度分别为413~593 ℃和336~498 ℃,盐度分别为19.1%~49.7% NaCleq和15.7%~53.3% NaCleq;Ⅲ阶段成矿流体均一温度为148~398 ℃,并具有低盐度(2.1%~9.6% NaCleq)与高盐度(35.5%~65.3% NaCleq)共存的特征;Ⅳ阶段成矿流体具有低温(132~179 ℃)、低盐度(3.4%~10.4% NaCleq)特征。根据流体包裹体的微观特征并结合矿区的宏观地质特征,认为流体不混溶(沸腾)是导致本矿区金属沉淀成矿的主要机制。     

Accuracy verification of MODIS daily cloud-free snow cover products over the Tibetan Plateau北大核心CSCD

Snow cover is of great hydrological,ecological,and climatic significance in the Tibetan Plateau. MODIS snow products are widely used at present but are seriously affected by clouds. Scholars at home and abroad have developed a variety of cloud removal products for raw MODIS daily snow products,but the accuracy of these products in the Tibetan Plateau has not been evaluated comprehensively. Therefore,this paper uses Landsat-8 data with high resolution as the reference value to conduct systematic verification of three datasets of cloud-free snow products released on a daily basis. The results show that compared with the two sets of products （M*D10A1GL06 and MODIS_Dysno_Cloudfree）,which are produced based on raw MODIS daily snow cover product realized by NSIDC（National Snow and Ice Data Center）,the MODIS CGF SCE product produced based on MODIS surface reflectance data,has a great advantage in snow identification accuracy. The MODIS CGF SCE product optimized the NDSI threshold for different land cover types. Although the accuracy of snow identification was significantly improved,the problem of large snow identification error in forest areas was still not effectively resolved,and there was a high underestimate error. © 2022 Science Press (China).     

低渗透储层孔隙结构影响因素及其定量评价——以准噶尔盆地金龙2地区二叠系上乌尔禾组二段为例

利用薄片、扫描电镜、物性及压汞等资料,对金龙2地区三叠系上乌尔禾组二段(P_3w_2)低渗(含砾)砂岩储层孔隙结构及其影响因素进行了研究。结果表明:①乌二段属于低孔低渗储层,根据压汞曲线参数特征,将其孔隙结构划分为三类。②成岩作用控制了储层孔隙结构,定量计算结果表明,乌二段视压实率平均为53.54%、视胶结率平均为47.53%、视溶蚀率平均为13.84%、视微孔隙率平均为60.67%,表现为中等压实、中—强胶结、弱溶蚀、微孔发育等特征;压实作用、胶结作用、微孔隙发育主要控制储层孔隙结构的形成。③引入成岩综合指数来定量表征各种成岩作用的综合强度,其与储层孔隙结构参数(排驱压力、分选系数)以及储层品质因子(RQI)具有较好的统计相关关系,Ⅰ类储层成岩综合指数大于8%,Ⅱ类储层成岩综合指数为2%～8%,Ⅲ类储层成岩综合指数小于2%。因此,可以利用成岩综合指数定量评价储层孔隙结构。     

大兴安岭北段宝兴沟金矿床成矿流体特征及矿床成因

宝兴沟金矿床是大兴安岭北部上黑龙江成矿带内大型金矿床之一,矿体主要产于下侏罗统二十二站组砂岩与早白垩世石英闪长岩、闪长玢岩内外接触带内,其热液成矿作用可划分为：黄铁矿±毒砂-石英（Ⅰ）、多金属硫化物-石英（Ⅱ）及少硫化物-碳酸盐（Ⅲ）3个阶段。流体包裹体岩相学研究表明：Ⅰ阶段矿石主要发育气液两相（LV）、少量含CO₂三相包裹体（H_CO2）及富气相包裹体（FV）;Ⅱ阶段矿石中主要发育LV及少量H_CO2包裹体;Ⅲ阶段矿石中只发育LV包裹体。测温结果显示：Ⅰ、Ⅱ阶段包裹体总体均一温度峰值集中于225.00~300.00℃,盐度（w（NaCl））为2.00%~10.00%;Ⅲ阶段均一温度峰值集中于175.00~225.00℃,盐度为4.00%~8.00%;成矿流体为简单的含CO₂中低温、低盐度的NaCl-H₂O热液体系,总体具有从成矿早期到晚期均一温度、盐度逐渐降低的特征。氢、氧同位素分析结果显示,Ⅰ、Ⅱ阶段成矿流体δD_SMOW为-131.00‰~-108.00‰、δ¹⁸O_SMOW为1.00‰~4.00‰,Ⅲ阶段δD_SMOW为-108.00‰、δ¹⁸O_SMOW为-1.89‰,表明早期以岩浆水为主,晚期逐渐演化为与大气降水混合热液。矿石中黄铁矿（毒砂）δ³⁴S_V-CDT为1.50‰~4.20‰,显示其物质来源以深源岩浆为主。综合分析认为,区内金成矿作用与早白垩世（石英）闪长岩、闪长玢岩侵入活动有直接关系,矿床属中低温岩浆热液成因类型。     

小秦岭东桐峪金矿床的流体包裹体研究

东桐峪金矿床位于小秦岭金矿田的中西部,其含金石英脉受韧性剪切构造带的控制。该矿床的构造-成矿过程可划分为4个阶段:Ⅰ黄铁矿-乳白色石英脉阶段;Ⅱ灰白色石英-黄铁矿阶段;Ⅲ石英-多金属硫化物阶段;Ⅳ石英-碳酸盐阶段。相对于小秦岭地区的其他金矿床,东桐峪金矿床的流体包裹体研究资料相对缺乏。文章表明,该矿床内的流体包裹体类型主要为CO₂-H₂O包裹体和水溶液包裹体,见少量纯液相CO₂包裹体。显微测温表明,Ⅰ阶段的构造-成矿流体以中温、富CO₂等挥发分为特征,包裹体均一温度为221~392℃,盐度w(NaCl_eq)为5.5%~7.9%,密度为0.84~0.93 g/cm³;Ⅱ阶段和Ⅲ阶段以CO₂-H₂O±CH₄流体为主,包裹体均一温度为205~350℃(Ⅱ阶段)和224~271℃(Ⅲ阶段),盐度w(NaCl_eq)集中于5.1%~7.1%,密度为0.83~0.96 g/cm³;Ⅳ阶段的流体演化为中-低温、低盐度的盐水溶液体系,包裹体均一温度为175~185℃。文章对该矿床各成矿阶段的压力进行了估算,Ⅰ、Ⅱ、Ⅲ阶段的流体最小捕获压力分别为123~160 MPa、160~170 MPa、170 MPa左右。