基于短波红外技术的西藏珠勒地区蚀变矿物填图.pdf

第45卷第4期2023年8 月四川冶金Sichuan MetallurgyVol.45No.4Aug.,2023文章编号：10 0 1510 8(2 0 2 3)0 40 0 2 50 8基于短波红外技术的西藏珠勒地区蚀变矿物填图余杰,刘学,王琳,王洪飞，司永强（成都理工大学地球科学学院，四川成都6 10 0 59）摘要：蚀变矿物组合对于热液型矿床勘查工作具有重要的指导意义，短波红外技术对于蚀变矿物的识别鉴定具有效率更高的优势。本文使用短波红外技术，基于GF-5号AHSI高光谱影像对明矾石和高岭石进行蚀变填图，结合Landsat8蚀变信息提取结果，通过ASDFeildSpec野外光谱辐射仪对岩石样品进行光谱测试，通过SWIR光谱分析、野外实地路线考察、岩相学显微镜下鉴定矿物等方法验证了短波红外技术对GF-5号高光谱数据蚀变填图的有效性。将蚀变与地层、岩性和构造结合分析，最终得出，珠勒地区明矾石化和高岭石化受断裂构造影响较大，明矾石类型主要为钾明矾石，光谱吸收特征在2 2 0 9nm附近，钾明矾石和高岭石伴生出现可作为珠勒地区高硫型浅成低温热液矿床的找矿标志，为下一步地质找矿提供依据。关键词：西藏；珠勒地区；蚀变；短波红外；高光谱；二次主成分分析中图分类号：P627Mapping of altered minerals in Zhule Area,Tibet basedYU Jie,LIU Xue,WANG Ling,WANG Hongfei,SI Yongqiang(Chengdu University of Technology,Chengdu 610059,Sichun,China)Abstract:Altered mineral assemblages have important guiding significance for the exploration of hydrothermal de-posits.The short-wave infrared technology has the advantage of higher efficiency for the identification of alteredminerals.In this paper,combined with the extraction results of Landsat8 alteration information,short-waveinfrared technology was used to map alteration of alunite and kaolinite based on GF-5 AHSI hyperspectral image.Spectral tests were conducted on rock samples by ASD FeildSpec field spectral radiometer.The effectiveness ofshort-wave infrared technology for alteration mapping of GF-5 hyperspectral data was verified by SWIR spectral a-nalysis,field route investigation and mineral identification under petrographic microscope.By combining the altera-tion with stratigraphic,lithology and structure analysis,it is concluded that the alunite mineralization and kaolinitemineralization in Zhule area are greatly affected by fault structure.The main type of alunite is potassium alum,andthe spectral absorption characteristics are around 22o9 nm.The associated occurrence of potassium alunite and kao-linite can be used as the prospecting indicator of high sulfide epithermal deposit in Zhule area.It provides the basisfor further geological prospecting.Key words:Tibet;Zhule area;alteration;short wave infrared;hyperspectral;secondary principal component analy-sis珠勒矿区地处西藏阿里地区革吉县，大地构造位置处于青藏高原冈底斯一腾冲陆块与喜马拉雅陆块的衔接部位，属于冈底斯成矿带的西段。目前，带内已发现斑岩型、矽卡岩型、浅成低温热液型等多种矿床1-3,开展了水系沉积物地球化学、流体包裹体同基金项目：本文由中国地质调查局项目 战略性矿产重点远景区矿产地质调查（编号：DD20221684)资助。作者简介：余杰（1998 一），女，硕士研究生，地质资源与地质工程（资源与环境遥感）专业，E-mail：12 6 310 7 57 2 q q.c o m。文献标识码：Aon short wave infrared technology位素4-6 、重磁7.8 等一系列工作,对带内不同类型矿床的成矿物质来源、成矿流体和成矿年龄进行了研究。珠勒矿区受复杂地貌条件限制，传统地质工作存在一定滞后性，遥感技术为该区域的地质勘查工作带来了新的契机。胡文泰9 发现珠勒矿区存在明显的26羟基、铁染和碳酸盐蚀变；但基于卫星遥感蚀变填图的精度如何？地面实际矿物光谱特征及矿物类型是否与航天遥感解译信息一致？这些问题不得而知。作者基于实地地质勘查工作，以星一地一体化高光谱遥感技术为手段，开展珠勒矿区的光谱勘查，以期解决区内精细化蚀变填图的问题，确定找矿靶区。Hunt10对不同岩石在可见光一近红外谱段的光谱特征进行了研究，发现光谱特征变化由矿物中离子的电子跃迁、基团振动、电荷振动等引起。络阴离子由于弯曲振动，导致阴离子在2 0 0 0 2 50 0 nm存在诊断特征；阳离子，如 Fe2+、Fe 3+、M n+等，由于发生电子跃迁，在40 0 12 0 0 nm产生诊断性吸收特征11)。短波红外光谱(SWIR,13002 50 0 n m)具有成本低、识别速度快等优点，对火山、热液成矿系统中的中低温矿物具有特殊的敏感性，针对斑岩成矿系统矿床的研究已呈现出显著的优势12-16 。针对高硫化浅成低温热液矿床的研究,Zhouet al.17总结了明矾石的光谱变化规律，唐敏慧等18 、郭娜等19 确定了西藏铁格隆南矿床的高硫化矿体蚀变分带特征。综上所述，作者基于前人研究成果，以GF-5号蚀变填图为基础，结合Landsat8OLI多光谱提取的羟基、铁染蚀变信息结果，在充分野外调研和地面光谱测量后，精准验证GF-5蚀变填图结果，开展地面精细化蚀变填图，完成矿区找矿潜力分析。1研究区地质背景珠勒位于西藏阿里地区革吉县东南约6 0 km处，属革吉镇管辖。研究区属特提斯造山系的东部与阿尔卑斯一喜马拉雅巨型构造带东段，冈底斯一腾冲陆块与喜马拉雅陆块的衔接部位，紧邻雅鲁藏布江结合带西段。研究区出露地层较为简单，主要分布地层新近系、古近系、白垩系，其次为第四系（见图1）。第四系（Qh pal）岩性主要为冲洪积砂、沙砾石、粉砂质粘土、砾石砂层、粉砂、粉土层；新近系喜山阶（z x m 元Ni）主要为浅肉红色二长花岗斑岩；上白垩系燕山阶存在两种岩性，浅灰色中粒二长花岗岩、灰色一浅灰色中粒一中细粒花岗闪长岩；新近系中新统有布噶寺组（Ni6）,岩性为粗面岩一玄武岩，嘎扎村组一段(Nigl）岩性为粗安质晶屑熔结凝灰岩、粗安岩，二段(N1g）岩性为粗安质角砾凝灰岩、安山岩；古近系始新统林子宗群有帕那组（E2p）、典中组四川冶金（Ei-2 d），帕那组岩性主要为火山角砾岩、集块岩，典中组岩性主要为安山岩；上白垩统的江巴组(K2jb）岩性主要为流纹岩一流纹质凝灰岩，上则弄群(K,Z)主要为安山岩、玄武岩。典中组与上覆帕那组呈整合接触，与上覆布寺嘎组呈角度不整合，布噶寺组呈角度不整合覆盖于典中组之上。2kmQhNN,bN.gN,bQhN.g我拉N.gN.gEdQ1N.g7131一第四系；2 一上白垩统竞柱山组；3一新近系二长花岗斑岩；4一上白垩统二长花岗岩5一上白垩统花岗闪长岩；6 一新近系中新统粗面岩一玄武岩；7 一新近系中新统粗安质角砾凝灰岩、安山岩；8 一古近系始新统火山角砾一集块岩；9一古近系始新统安山岩；10 一上白垩统流纹岩一流纹质凝灰岩；11一下白垩统安山岩一玄武岩；12 一性质不明断层；13一推测断层；14一整合地质界线；15一不整合地质界线；16 一地名图1西藏珠勒地区地质简图2数据来源及预处理2.1楼数据来源本文使用Landsat8 OLI多光谱数据以及 GF-5号AHSI高光谱数据。OLI数据为L1级数据，产品号：LC81440382015269LGN01；G F-5号数据为L1级产品，产品号为GF5_AHSI_E81.40_N31.81_20190613_005833_L10000047709。两种数据含云量均小于5%，满足蚀变弱信息的提取。Landsat8卫星的陆地成像仪（Operational LandImager，O LI)传感器的多光谱影像，影像波段参数见表1。OLI多光谱影像共11个波段，其中，可将光一近红外一短波红外范围的波段为bandlband7,光谱分辨率为30 m，常用于提取铁染类、羟基类蚀变矿物2 0 1，全波段（pan)具有较高的空间分第45卷QhK,ZN.gN.gQEPOhEMdN.gE.PNgEd230K4E.P8Eid91415EPEPN,bK.jbN.bK.jbNb5N.g6K.jb10167oKK,Z11KEd12第4期辨率，空间分辨率为15m，热红外波段（TIR）为band10b a n d 11,空间分辨率较低，为10 0 m。本次主要利用可见光到短波红外光谱范围的波段进行铁染、羟基等蚀变异常的提取，影像数据。GF-5号可见短波红外高光谱相机（theAdvanced Hyperspectral Imager，A H SI)数据是兼顾宽幅宽、宽谱段的高光谱遥感影像，能较好地应用于蚀变填图2 1.2 2 。该数据包含近红外(VNIR)波段和短波红外(SWIR)两个光谱范围的波段，光谱范围为40 0 2 50 0 nm，空间分辨率皆为30 m，幅宽达6 0km。光谱分辨率方面，短波红外的光谱分辨率略低于近红外波段，其中，近红外波段的光谱分辨率为5nm，短波红外的光谱分辨率为10 nm。本次使用的数据为短波红外光谱范围的波段，数据光谱分布范围为10 0 42 513nm,共18 0 个波段。表 1 Landsat8 OLI 波段信息波段波段名称BandlCrostaBand2BlueBand3GreenBand4RedBand5NIRBand6SWIR1Band7SWIR2Band8PanBand9CirrusBand1oTIRSBand11TIRS2.2预处理2.2.1Landsat8 OLI多光谱预处理本次工作区所使用的Landsat8 OLI数据采用了辐射校正、FLAASH大气校正、归一化、图像裁剪等预处理。FLAASH大气校正可一定程度上减少大气对于影像数据带来的影响，但该模型校正后得到的影像像元值为反射率的10 0 0 0 倍，因此后续进行了归一化处理得到0 1.0 范围的反射率值，为后续蚀变信息的提取打下基础。2.2.2GF-5 AHSI高光谱预处理GF-5 AHSI的L1级的SWIR波段的数据进行辐射定标、波段剔除、FLAASH大气校正、正射校正、图像裁剪等预处理（见图2）。高光谱数据预处理主要在ENVI5.3软件上完成。辐射定标是利用GF-5AHSI数据的增益(Gain)和偏移值（Offset)，将图像的灰度值转换为Sichuan Metallurgy绝对辐射率值。波段剔除主要剔除水汽吸收波段和低信噪比波段，以此降低数据处理量，并减少噪音等无效信息对有用信息的干扰。GF-5 AHSI高光谱数据的SWIR光谱范围波段受水汽吸收影响严重，且信噪比低波段对蚀变信息的提取作用不大，以上两种波段将被去除。本次提取的蚀变矿物（明矾石、高岭石)主要吸收特征表现在SWIR的130 0 2 50 0nm波谱范围内，去除波段如表2 所示，包含受水汽吸收影响的波段SWIR4350 以及SWIR96112，共2 5个波段，分别对应光谱范围1359 1418 nm、18051940nm；低信噪比波段主要有17 个波段，对应波谱范围为1351443nm、19 48 19 56 n m、19 9 92 0 15n m、2 4452 513n m，分布不集中。表2 珠勒地区GF-5号短波红外波段去除波段类型波段96112水汽吸收4350波长范围/m空间分辨率/m0.430.45300.450.51300.530.59300.640.67300.850.88301.571.65302.112.29300.500.68151.361.383010.611.1910011.512.5110027.波长范围/nm波段数量180519402513591418421350525314351443低信噪比11311411912117218014151剩余波段54951151181221713蚀变信息提取3.1多光谱羟基、铁染蚀变信息提取主成分分析法（Crosta方法）提取羟基、铁染等蚀变信息已经较为成熟2 3。研究区影像中包含了第四系、水体、雪等干扰因素，常规处理方法为“掩膜十主成分分析”，但当掩膜面积较大时，蚀变信息也可能被误除，而二次主成分分析可以有效降低干扰信息影响，增强蚀变弱信息2 4。因此，本文对研究区采用二次主成分分析法一使用Crosta方法、波段比值法提取铁染、羟基蚀变信息的基础上进行主成分分析。代表性铁染蚀变矿物有褐铁矿、黄钾铁矾；羟基类代表性蚀变矿物有AI-OH类的白云母、伊利石、高岭石、明矾石、蒙脱石以及Mg-OH类的绿泥石、绿帘石等等2 5。利用UGSG波谱库中上述几种典型蚀变矿物的光谱曲线，基于OLI进行波谱重采样，结果如图2 所示。1948195619992015244525131004134214261451179719651990194024371741142450280.80.4-0.2-0.010.5地开石伊利石白云母高龄石蒙脱石1.0(d)0.80.40.2-0.0+0.40.60.81.01.21.41.61.82.02.2wavelength(m)a-AI-OH类蚀变矿物的USGS标准光谱曲线；b-Mg-OH类蚀变矿物的USGS标准光谱曲线；c-铁染类蚀变矿物的USGS标准光谱曲线；d-基于OLI重采样的A1-OH蚀变矿物光谱曲线；e-基于OLI重采样的Mg-OH蚀变矿物光谱曲线；f-基于OLI重采样的铁染类蚀变矿物光谱曲线从图2 a可以看出，Al-OH的吸收特征位于2.10m附近，Mg-OH则表现为2.3m附近（图2b)，铁染（主要为Fe3+)则表现为0.4m和0.9 m附近(图 2 c)。从图2 d、e 可以看出，羟基类矿物在0.430.45m(OLI1)、2.112.2 9 m(OLI7)具有吸收特征,在 0.530.59 m(OLI3)、1.57 1.6 5 m(OLI6)具有反射特征，且OLI7为强吸收；OLI6为强反射；而在图 2 f 中,铁染在 0.8 50.8 8 m(OLI5)与 2.112.29m(OLI7)具有强吸收特征，在0.6 40.6 7 m(OLI4)与1.57 1.6 5m(OLI6)具有反射特征,且OLI6 相对于OLI5表现为强反射。由此可见，羟基和铁染都在 OLI6表现为强反射,在OLI7都表现为强吸收，存在信息重叠。因此，使用PCA（O L I 1、OLI3、O L I 6、O L I 7)提取羟基时混合有铁染信息；使用PCA(OLI4、O L I5、O L I6、O L I7)提取铁染信息则混合有羟基信息。而使用波段比值法OLI6/OLI7（下面表示为Vi)可增强羟基蚀变信息；OLI6/OLI5（下面表示为V）可增强铁染蚀变信息。表3主成分特征向量矩阵PCA(OLI1、O L I3、O L I6、O L I7)EigenvectorBand1PC1-0.26907770.429928PC2-0.65971440.556285PC30.050048-0.229008-0.513916PC4-0.699907四川冶金地开石一伊利石高龄石一一白云母一一蒙脱石1.0(a)第45卷一绿帘石一一绿泥石1.0(b)0.8-0.20.01.01.5wavelength(m)图2 典型蚀变矿物USGS波谱库光谱曲线及基于OLI重采样光谱曲线Band3Band6-0.6815710.4722480.673249-0.219849黄钾铁矾一一褐铁矿1.0(c)0.80.6-0.40.20.02.02.50.51.0(e)0.8-0.2-0.00.40.60.81.01.21.41.61.82.02.2Band7-0.527470.1797390.8251920.0923721.0wavelength(m)绿帘石绿泥石wavelength(m)表4主成分特征向量矩阵PCA(OLI4、O L I 5、O L I 6、O L I 7)EigenvectorBand4PC1-0.43712-0.49783-0.59397PC20.61271-0.442460.556234PC30.211941PC40.623373根据上述蚀变信息的反射吸收特征，结合表3特征值可知,羟基蚀变信息应该在OLI1与OLI7同号、OLI3和OLI6同号，且OLI6和OLI7异号的主分量中，即PC3（下面表示为V4），且表现为负异常高特征值。同理，铁染蚀变信息则应该在OLI4与OLI6同号、OLI5与OLI7同号，且OLI6与OLI7异号的主分量中，即PC4（下面表示为V3）,亦表现为负异常高特征值（见表4）。而波段比值法得到的羟基、铁染蚀变皆表现为正异常且高特征值。因此，进行二次主成分分析PCA（O L I 6/O L I 7，O L I 6/OLI5,PC4(OLI1,OLI3,OLI6,OLI7),PC3(OLI4,OLI5,OLI6,OLI7)）。得到特征向量矩阵如表5所示,PC3与Vi和V2呈负相关,与Vs和V呈正相关。因此,PC3主要包含铁染蚀变和羟基蚀变的组合信息。选则二次主成分分析得到的PC3同波段比值法、一次主成分分析得到的主分量进行RGB(PC3，Vi，Vs)假彩色合成，红色表示“铁染十羟基”混合蚀变信息，在此基础上，圈定了蚀变靶区（见图3)3.2高光谱蚀变矿物填图使用最小噪音变换MNF、纯净像元指数PPI、N1.52.02.50.51.00.80.60.40.20.0-0.40.60.81.01.21.41.61.82.02.2wavelength(m)Band5-0.30881-0.49807-0.678990.299551.0wavelength(m)黄钾铁矾揭铁矿Band61.52.02.5Band7-0.456390.3455680.782078-0.24627第4期维可视化、光谱分析等提取端元光谱，再使用光谱角填图法SAM进行明矾石和高岭石的蚀变填图（见图4）。波段剔除主要剔除水汽吸收波段和低信噪比波段，以此降低数据处理量，并减少噪音等无效信息对有用信息的干扰。光谱分析使用的波谱库为实测光谱，结合SAM填图法2 6 ，可以使填图结果更准确。表5二次主成分分析特征向量矩阵Sichuan Metallurgy布，并在断层交界附近集中分布。高岭石分布地层较广，主要集中分布在典中组（E1-2d）和布噶寺组（Ni 6）,零星分布在帕那组（E2）、则弄群（KZ）；岩性方面主要分布于安山岩、粗面岩、玄武岩等中基性火山岩中，构造方面，在F断层和次断层的交界处分布集中，并垂直于Fi断层的走向，呈线状分布(图 5b)。8105（a)298110105(b)8110EigenvectorPC1PC2PC3PC481008100预处理1端元提取1蚀变填图图4GF-5AHSIL1级SWIR数据处理流程图对珠勒地区使用GF-5号AHSI进行明矾石和高岭石的蚀变填图结果如图5a。可见，明矾石化主要分布在研究区西南部，高岭石主要沿NW-SE向分布，北部较为聚集，南部零星分布。对比多光谱提取的“羟基十铁染”蚀变靶区，填图结果较为吻合。结合以上蚀变填图和地层岩性、构造分析：明矾石主要分布在新近系中新统，噶扎村组（Nig）分布较为集中,而在古近系始新统典中组（Ei-2d）零星分布;岩性方面，明矾石主要分布在安山岩中；构造方面，明矾石化蚀变主要沿垂直于NEE向断裂分V10.8197110.5702940.05195-0.082820.011781-0.013250.6249881058110NO”羟基+铁染“蚀变81058110图3 1RGB(PC3,V1,V3)假彩色合成图SWIRGF-5MNFV20.56568-0.82035-0.019060.77712812km8115辐射定标V30.071493大气校正FLAA.SH1PPIN维可视化SAMV40.054484-0.03780.6217040.78044正射校正光谱分析81*10Qh1N7a-GF-5号AHSI高光谱蚀变矿物填图；b-蚀变矿物遥感地质综合图16-明矾石；17-高岭石；18-“羟基十铁染”混合蚀变靶区。其他符号含义同图1图5珠勒地区明矾石、高岭石蚀变矿物分布4野外验证4.1知短波红外光谱测试本次采集并测试岩石样品6 1件，测量仪器为美国ASDFeildSpec便携式光谱扫描仪。通过光谱地质解译裁剪图像810523K4N.b5N689K10K.Z11121516专家(The Spectral Geologist，T SG)进行数据处理与解译，测得研究区蚀变矿物主要有明矾石、铁白云石、黑云母、方解石、绿泥石、水铝石、白云石、绿帘石、角闪石、黄钾铁矾、蒙脱石、绢云母（白云母、钠云母、多硅白云母）、伊利石、菱铁矿、电气石、黔帘石等。野外采集得到的岩石矿物光谱通常是多种岩石矿物的综合结果，为得到较为纯净的单一矿物光谱，需进行光谱解混。通常情况下，大尺度的光谱混合通常被认为是一种线性组合，小尺度的则被认为是非线性组合2 7。对于岩矿信息提取这一类大尺度的光谱混合，本次采用线性光谱解混的方式对野外采集的实测光谱进行光谱解混。最小二乘法是线性光谱解混中一种常用的方法，精度较高2 8 ,采用此方法进行光谱解混，得到图6 两种蚀变矿物的光谱曲线。解混得到的光谱曲线与USGS岩石矿物光谱库中的光谱曲线波形吻合，从光谱角度证明了蚀变矿物识别结果可靠。811017O1830明矾石是一种在pH极低、中低温环境下生成的蚀变矿物2 9 。明矾石的化学成分主要为KAl3(SO4）2（O H)6，因其含有的A1-OH基团，而在2 10 0nm 附近存在诊断性光谱特征11,而由于 K 被 Na和不同程度替代时，吸收特征波长位置发生飘移30 。本次研究区所采集的岩石样品中，含有明矾石的样品一共有5件，吸收特征波长在2 2 0 9 nm附(a)0.910.80.70.60.50.40.30.20.120004.2野外路线调查根据路线地质调查，研究区古近系典中组火山岩分布区为蚀变中心，其地表发育强烈的硅化、绢云母化、泥化、钾化、绿泥石化、黄铁绢英岩化、褐铁矿化。蚀变可分为上下两部分，上部主要为红褐色强硅化褐铁矿化蚀变岩，局部赤铁矿化，新鲜面内多孔洞，泥化类型以明矾石为主，另发育叶腊石、高岭石、石膏等；下部为灰白色、淡黄褐色硅化泥化蚀变岩，四川冶金近(图6 a)，由此可见，珠勒地区的明矾石类型主要为钾明矾石。高岭石也是在酸性低温环境中生成的，与明矾石常一起出现的蚀变矿物。高岭石由于存在Al-OH,在2 2 0 5nm出现强吸收（图6 b），这与明矾石区分易发生混淆。但其在2 315nm、2 355n m、2 38 5nm三处出现连续吸收特征，可区别于明矾石。TK0711G1(b)TK0709GTK0708B1TK0705B1TK08-B1-明矾石USGS标准光谱(2214)(2216)(2209(2209）(2209)（2 2 0 8)21002200Wavelength(nm)图6 珠勒地区地表岩石样品中钾明矾石(a)和高岭石(b)光谱曲线4.3岩相学镜下鉴定将野外采集的地表岩石样品进行加工，打磨成光片、薄片，在显微镜下，根据岩相学方面的知识，在TK0708采样点找到了明矾石（见图7）。第45卷0.61DO311-B1-高岭石USGS标准光谱0.5-0.4-0.30.20.123002400(2315)(2335)(2355)25002000风化面多孔洞，新鲜面长石泥化、绢云母化，局部发育黄钾铁矾。围岩蚀变具有明显的分带现象，从中心向外表现为强硅化、明矾石化、绢云母化、泥化、绿泥石化。21002200Wavelength(nm)230024002500(a)(b）硅报？硅化带火山角砾岩流武岩泥化带(a)野外蚀变带分布；（b)明矾石岩石标本；（c)明矾石镜下照片（TK0705)图7 野外照片及镜下鉴定照片岭石伴生出现，可成为珠勒地区高硫型浅成低温热5结论液矿床的指示标志。(1)短波红外技术能方便快捷地厘定蚀变矿物(2)对比多光谱OLI提取的“羟基十铁染”混合类型,SWIR光谱分析得出珠勒地区主要存在钾明蚀变结果，结合光谱测试分析、野外勘查、镜下鉴定矾石，吸收特征位于2 2 0 9 nm附近。钾明矾石和高等方法，认为GF-5号高光谱影像能满足蚀变矿物精第4期细填图，为下一步的研究起到指导作用。(3）本文基于ASD实测光谱，采用GF-5号AHSI高光谱地SWIR波段影像，提取了明矾石和高岭石的蚀变信息，结合地层、岩性、构造分析得出，明矾石在新近系嘎扎村组的安山岩中，在断层交接处附近集中分布，高岭石则在古近系典中组的安山岩和新近系布噶寺组粗面岩一玄武岩中集中分布，且都在断层的交界处。认为断层构造对此次蚀变有着很大影响。基于短波红外技术提取的蚀变集中分布之处，推测该处具有较大的成矿潜力。参考文献：1黄瀚霄，李光明,刘洪，等。冈底斯成矿带西段首次发现低硫化型浅成低温热液型矿床一罗布真金银多金属矿床J.中国地质，2 0 18，45（3）：6 2 8 6 2 9.2黄黄瀚霄，张林奎，刘洪，等.西藏冈底斯成矿带西段矿床类型、成矿作用和找矿方向J.地球科学，2 0 19，44(6):18761887.3刘洪，黄瀚霄，张林奎，等.西藏冈底斯成矿带西段鲁尔玛晚三叠世斑岩型铜(金)矿点的发现及意义J.沉积与特提斯地质，2 0 2 1，41(4)：59 9 一6 11.4刘洪，夏祥标，黄瀚霄，等.西藏冈底斯成矿带西段学修玛尔幅水系沉积物地球化学统计分析与找矿前景J.桂林理工大学学报，2 0 19，39（4）：8 47 8 55.5文刘洪，李光明，黄瀚霄，等.冈底斯成矿带西段鲁尔玛斑岩型铜（金)矿床的成矿物质来源研究J.矿床地质，2019,38(3):631643.6ZZHUANG Li,LANG Xinghai,RICKLEMAN Daniel,et al.Age and genesis of the Pusangguo skarn Cu-dom-inated polymetallic deposit，,G a n g d e s e m e t a l l o g e n i cbelt,TibetJ.Journal of Asian Earth Sciences,2019,169210227.7胡斌.基于重磁数据的冈底斯成矿带深部结构分析与成矿预测D.中国地质大学（北京)，2 0 19.8YUE Sheng,JIN Sheng,COMEAU Matthew-J,et al.Con-trols on the metallogenesis of the Lhasa-Mozugongkadistrict,Gangdese Belt,Tibetan Plateau:Constraints onmelt distribution and viscosity from the 3-D electrical struc-ture of the lithosphere J.Ore Geology Reviews,2022,145104881.9故胡文泰.基于Aster数据的遥感蚀变信息提取一西藏珠勒地区为例J.科学技术创新，2 0 2 2，（5）：4649.10 GRAHAM-R Hunt.Nearinfrared（1.3-2.4）mspectra of alteration mineralsPotential for use in re-mote sensingJ.Geophysics，19 7 9,44(12):19 7 4Sichuan Metallurgy11甘甫平，王润生，马噶乃.基于特征谱带的高光谱遥感矿物谱系识别J.地学前缘，2 0 0 3，（2）：445一454.12 NA Guo,THOMAS Cudahy,TANG Juxing,et al.Mapping white mica alteration associated with theJiama porphyry-skarn Cu deposit,central Tibet usingfield SWIR spectrometry.2017.13 NA Guo,GUO Wenbo,SHI Weixin,et al.Character-ization of Illite Clays associated with the Sinongduo lowsulfidation epithermal deposit,Central Tibet using fieldSWIR spectrometry.2020:103228.14陈华勇，张世涛，初高彬，等.鄂东南矿集区典型矽卡岩一斑岩矿床蚀变矿物短波红外（SWIR)光谱研究与勘查应用J.岩石学报，2 0 19，35（12）：36 2 9 36 43.15郭娜，郭科，张婷婷，等.基于短波红外勘查技术的西藏甲玛铜多金属矿热液蚀变矿物分布模型研究J.地球学报，2 0 12，33(4)：6 416 53.16汪重午，郭娜，郭科，等.基于短波红外技术的斑岩一矽卡岩型矿床中绿泥石蚀变分布特征研究：以西藏甲玛铜多金属矿为例J.地质与勘探，2 0 14，50（6）：1137-1146.17 YAN Zhou,FAN Feipeng,XING Guangfu,et al.Characteristics and genesis of the Fanshan lithocap,Zhejiang Province:Exploration implications from thelargest alunite deposit of China.2022:105038.18 唐敏惠，周涛发，范裕，等.高硫成矿系统中明矾石的稳定同位素特征J.合肥工业大学学报（自然科学版），2011,34(9)：140 2 140 7.19郭娜，黄一入，郑龙，等.高硫一低硫化浅成低温热液矿床的短波红外矿物分布特征及找矿模型一一以西藏铁格隆南（荣那矿段）、斯弄多矿床为例J.地球学报，2017,38(5):767778.20汪子义，张廷斌，易桂花，等.西藏冲江一厅宫斑岩型铜矿区遥感蚀变信息提取与找矿预测.地质与勘探，2016，52(6)：112 9 1137.21刘银年，孙德新，胡晓宁，等.高分五号可见短波红外高光谱相机设计与研制J.遥感学报，2 0 2 0，2 4（4）：333344.22 BEI Ye,TIAN Shufang,CHENG Qiuming,et al.Ap-plication of Lithological Mapping Based on AdvancedHyperspectral Imager（A H SI）I m a g e r y O n b o a r dGaofen-5（G F-5)Sa t e lli t e J.Re m o t e Se n s i n g,2 0 2 0,12(23).2 3毛晓长，刘文灿，杜建国，等.ETM+和ASTER数据在遥感矿化蚀变信息提取应用中的比较一一以安徽铜陵凤凰山矿田为例J.现代地质，2 0 0 5，（2）：30 9 一314.24吴志春，叶发旺，郭福生，等.主成分分析技术在遥感蚀31-1986.32变信息提取中的应用研究综述J.地球信息科学学报，2 0 18，2 0(11)：16 4416 56.25唐超,周可法,张楠楠，等.基于Landsat-8OLI和AS-TER数据集成和融合的矿化蚀变信息提取：以包古图斑岩型铜矿为例J.地质科技情报，2 0 18，37（6）：211-217.26 F-A Kruse,Lefkoff A-B,Boardman J-W,et al.Thespectral image processing system(SIPS)interactivevisualization and analysis of imaging spectrometer dataJJ.Remote Sensing of Environment,1993,44(2):四川冶金145-163.2 7 张兵.高光谱图像处理与信息提取前沿J.遥感学报，2016,20(5):10621090.2 8 张帅洋，华文深，应家驹，等.高光谱线性解混研究进展J.激光杂志，2 0 2 1，42（0 3）：17 2 1.29 E-Stoffregen Roger.Genesis of acid-sulfate alterationand Au-Cu-Ag mineralization at Summitville,ColoradoJ.Economic Geology,1987,82(6).30郭娜.甲玛斑岩一矽卡岩型铜矿床蚀变矿物组合研究D.成都理工大学，2 0 12.第45卷(上接第 17 页)5结论依据SAE-ChinaJ3201一2 0 13汽车覆盖件板材抗凹性能试验方法，完成DC53D十Z抗凹性能测试，初始刚度和失稳凹陷载荷分别为8 6.3N/mm和87.6N。根据试验工况建立了抗凹试验仿真模型并进行计算，初始刚度仿真结果8 6.3N/mm，仿真误差0.7%，失稳凹陷载荷仿真结果8 9.5N仿真误差2.2%，仿真精度高。基于该模型，完成DC53D十Z、HC180BD十Z、H C 30 0/50 0 D PD 十Z三种汽车板以及0.55mm、0.6 m m 0.6 5m m、0.7 m m 四个厚度的抗凹性能仿真分析，获得了不同厚度汽车板的初始刚度和失稳凹陷载荷数据。通过数据分析：初始刚度与厚度有关，呈线性关系，方程为K;=Ct一97.9，其中C=285.4；失稳凹陷载荷与厚度和屈服强度均有关，失稳凹陷载荷随着厚度和屈服强度的增大而增大，失稳凹陷载荷F与厚度t、屈服强度s呈幂函数关系，方程为F，=C t，其中a=2.14,b=0.134,C=116.7。参考文献：1逐若东，程超，陈自凯.DP490双相钢车门外板应用研究J.锻压技术，2 0 2 0，45（5）：110 115.2冀鹏.HC300/500DP在