遥感技术作为矿产勘查的一种技术手段由于具有经济、高效等优势，在低植被覆盖区已被广泛应用并取得了较好的成效（张玉君等，2009; 耿新霞等，2010; Liu Lei et al.，2013，2016; 任广利等，2013; 尹芳等，2014; 刘德长等，2017; 王钦军等，2017）。已有研究涉及到众多类型矿床遥感找矿方法及应用，如岩浆分异型铁矿（耿新霞等，2010）、沉积变质型铁矿（王俊峰，2013）、火山岩型铁矿（陈劲松等，2012）、镁铁-超镁铁岩型铜-镍硫化物矿（杨建民等，2007; Liu Lei et al.，2014; Rogge et al.，2014）、斑岩型铜钼矿（Bishop et al.，2011; Pour et al.，2011; 姚佛军等，2012; Liu Lei et al.，2013，2017; Alimohammadi et al.，2015; 王钦军等，2017）、与石英脉有关铜金矿（Liu Lei et al.，2016）、造山带型金矿（Pour et al.，2013; 孙雨等，2015）、变质热液型金矿（张微等，2010）、沉积变质型铜多金属矿（黄智才等，2014）等，在取得良好找矿效果的同时为后续相关矿床找矿提供了方法借鉴。

胡兹达尔-拉斯贝拉铅锌成矿带为特提斯铅锌巨型成矿带的组成部分（图1），是巴基斯坦最为重要的铅锌成矿带，该成矿带的矿床类型主要为喷流-沉积型（Sedimentary exhalative，SEDEX）和密西西比河谷型（Mississippi valley-type，MVT），代表性矿床包括杜达（Duddar）和贡尕（Gunga）等（Ahsan et al.，1997; 张洪瑞等，2010; Song Yucai et al.，2019; 张辉善，2021）。目前针对胡兹达尔-拉斯贝拉铅锌成矿带及相应矿床的研究较少，制约了对该成矿带的理解及其潜力分析。随着“一带一路”倡议的实施，境外地质工作不断深入，大量中国企业在海外投资矿产勘查、开发，特别是胡兹达尔-拉斯贝拉成矿带储量最大的杜达矿床即为中国企业在开发。开展该成矿带遥感地质特征及找矿标志研究，可以为该成矿带后续找矿、中国企业投资提供依据。

因此，本研究以胡兹达尔-拉斯贝拉成矿带为研究区，通过合理选择遥感数据源和图像处理方法，对典型矿床样品光谱特征、影像特征分析，结合关键控矿要素研究，建立与沉积岩有关铅锌矿遥感找矿标志并开展成矿预测，不仅对该区找矿具有指导意义，形成的方法也可以用于特提斯铅锌巨型成矿带找矿勘查工作。

1 研究区地质概况

全球三大成矿域（环太平洋、特提斯、古亚洲）之一的特提斯成矿域，位于欧亚大陆与冈瓦纳大陆的交接部位，被认为是欧亚大陆南部一条东西展布的巨型构造带，也是地球上最年轻的造山褶皱带（张洪瑞等，2010）。特提斯成矿域赋存有一系列斑岩型Cu-Mo-Au矿床、与岩浆热液有关的W-Sn矿床、岩浆型铬铁矿、VMS型Cu-Pb-Zn矿床、浅成低温热液型Au-Hg-Sb矿床及与沉积岩有关的Pb-Zn矿床（张洪瑞等，2010; 吕鹏瑞等，2016，2020; Song Yucai et al.，2019; 张辉善，2021）。

其中，巴基斯坦境内的胡兹达尔-拉斯贝拉铅锌成矿带地处印度板块西北缘向阿拉伯板块俯冲拼合地带（图1）。区内广泛出露侏罗纪—白垩纪沉积地层，侏罗纪地层自下而上包括斯平沃组（Spingwar Formation，Js）、劳瑞莱组（Loralai Formation，Jl）和安吉拉组（Anjira Formation，Ja）。其中，斯平沃组主要为浅海相硅质碎屑岩（砂岩页岩互层夹少量灰岩），劳瑞莱组为厚层灰岩夹页岩，安吉拉组为深海相灰岩与泥岩互层（图2）。

研究区内已知的SEDEX型杜达铅锌矿和贡尕铅锌矿均赋存于中下侏罗统安吉拉组碳酸盐岩和碎屑沉积岩中（图2），矿体呈层状、似层状和网脉状产出，矿石矿物以闪锌矿、方铅矿、黄铁矿、白铁矿为主，脉石矿物主要包括方解石、石英、重晶石、绿泥石、磷灰石、长石、云母和有机碳，围岩蚀变主要有黄铁矿化、硅化、碳酸盐化、重晶石化、绿泥石化和钠长石化。其中黄铁矿化是矿区发育最广泛的蚀变，在围岩和铅锌矿中呈细粒或脉状分布，形成黄铁矿化围岩（Ahsan et al.，1997; Song Yucai et al.，2019; 张辉善，2021）。其中，杜达Pb-Zn矿床是胡兹达尔-拉斯贝拉成矿带储量规模最大的矿床（铅锌资源量1431万t，平均品位Zn 8.6%，Pb3.2%）（张辉善，2021）。野外调查表明层状矿体与安吉拉组地层关系密切，明显受沉积层位控制，矿体和围岩发生同步褶曲。在下部存在网脉矿体，部分与围岩整合接触，部分则为穿层矿体，受同沉积断裂构造控制。贡尕矿床位于胡兹达尔-拉斯贝拉地区北部，铅锌资源量686万t，平均品位Zn 11.4%、Pb2.1%（Ahsan et al.，1997）。含矿层位为安吉拉组下段，岩性为黑色炭质页岩、粉砂岩、泥岩夹灰岩。区内已知的苏迈MVT型铅锌矿床矿体主要产于中下侏罗统劳瑞莱组中，矿体沿灰岩层面或裂隙产出（Song Yucai et al.，2019）。

1 —第四系; 2—上新统:砂岩夹砾岩、页岩; 3—中新统:泥岩; 4—渐新统:砂岩、灰岩; 5—始新统:灰岩; 6—古新统:灰岩、砾岩、砂岩; 7—上白垩统:砂岩; 8—下白垩统:灰岩; 9—侏罗系:灰岩夹页岩; 10—上三叠统-下侏罗统:灰岩与页岩互层; 11—蛇绿岩; 12—地层界线; 13—断裂; 14—铅锌矿床; 15—铅锌矿化点

1 —Quaternary; 2—Pliocene: sandstone with conglomerate and shale; 3—Miocene: mudstone; 4—Oligocene: sandstone and limestone; 5—Eocene: limestone; 6—Paleocene: limestone, conglomerate and sandstone; 7—Upper Cretaceous: sandstone; 8—Lower Cretaceous: limestone; 9—Jurassic: limestone with shale; 10—Upper Triassic to lower Jurassic: limestone interlayered with shale; 11—ophiolites; 12—stratigraphic boundary; 13—fault; 14—Pb-Zn deposit; 15—Pb-Zn mineralized point

受印度板块向西北方向斜切碰撞影响，区内构造极其发育，广泛发育近N—S向、NE—SW向断裂及轴向NNW—SSE向和近N—S向的褶皱带。区内已知铅锌矿矿体以层状、似层状以及顺层脉状为主产出，受构造控制，与向斜地层和次级断裂关系密切。

2 数据与方法

2.1 数据

2.1.1 遥感影像

本研究采用多源遥感数据，即Landsat ETM数据、ASTER数据、资源一号02D高光谱数据。其中ETM数据覆盖整个成矿带，主要开展1∶25万遥感地质解译及信息提取; ASTER数据覆盖侏罗纪含矿地层，重点解译与铅锌矿化有关构造、岩性界线并提取矿化蚀变带; 对中部地区采用资源一号02D高光谱数据进行含矿地层自动提取。

（1）ETM影像。ETM在可见光/近红外（Visible to near infrared，VNIR）和短波红外（Shortwave infrared，SWIR）区间具有6个多光谱波段（空间分辨率30 m）和1个全色波段（空间分辨率15 m）。覆盖胡兹达尔-拉斯贝拉地区需2景ETM数据，轨道号:153/041、153/042，成像时间为2000年9月23日。

（2）ASTER影像。ASTER为美国、日本1999年12月联合发射的Terra卫星所携带的传感器，在VNIR（空间分辨率15 m）和SWIR（空间分辨率30 m）区间共有9个波段，光谱分辨率较高。特别是在SWIR区间有6个波段（1.60~2.43 μm），可依据岩石、矿物光谱特征选择各种假彩色合成及比值组合开展解译。本研究共处理、分析了覆盖胡兹达尔-拉斯贝拉地区重点区域的8景ASTER数据，数据基本覆盖了侏罗纪含矿地层。

（3）资源一号02D高光谱影像。资源一号02D为我国2019年9月12日在太原卫星发射中心成功发射的高光谱传感器，其在VNIR-SWIR区间共有166个波段，空间分辨率30 m，VNIR、SWIR区间光谱分辨率分别为10 nm、20 nm。获取了胡兹达尔-拉斯贝拉地区中部1景资源一号02D高光谱影像，成像时间为2019年12月14日，数据无云、雪，质量较好。

2.1.2 样品测试与分析

为了准确分析区内矿床的遥感特征并建立找矿标志，本研究选择胡兹达尔-拉斯贝拉铅锌成矿带内储量最大、成矿特征较典型的杜达铅锌矿和贡尕铅锌矿作为重点调查对象，开展野外矿化带取样、光谱测试、岩矿鉴定，据此确定矿化带的反射光谱特性。

通过对杜达矿区进行实测剖面调查、取样分析，表明矿区内出露地层主要为白垩系塞姆班组（Ks）、侏罗系安吉拉组（Ja）、劳瑞莱组（Jl）和斯平沃组（Js）（图3）。其中安吉拉组为主要含矿层位，主要岩性为灰岩、炭质泥岩、页岩以及粉砂岩。

采用SR-3500地物光谱仪（350~2500 nm，Spectral Evolution Inc.，USA）对杜达矿床的样品进行光谱测试，使用白色BaSO₄面板进行反射校准参考。本次共测试杜达矿区地表围岩、矿化带、矿体及坑道内样品30件，获得样品光谱（新鲜面、风化面）54条，代表性光谱见图4。安吉拉组各岩性单元中含矿层位的灰岩光谱均具有2200 nm和2350 nm吸收特征组合，可能是由于热液在流体循环过程中携带成矿物质的同时使灰岩发生钠长石化等蚀变，长石部分风化形成了黏土矿物所致。而未含矿层位斯平沃组和劳瑞莱组样品则不具有2200 nm的吸收特征。因此2200 nm和2350 nm吸收特征组合可以作为判断存在铅锌矿化有关蚀变的诊断性特征。

2.2 实验方法

2.2.1 影像增强方法

对ETM、ASTER和资源一号02D高光谱影像分别进行几何校正，利用FLAASH（Fast line-of-sight atmospheric analysis of spectral hypercubes，FLAASH）模型对3种影像进行大气校正（Berk et al.，1998）。

根据SEDEX型、MVT型矿床的蚀变特征及杜达矿区矿化带样品实测反射光谱，含矿层位样品光谱具有2200 nm和2350 nm组合吸收特征，对应于ETM数据的第7波段，ASTER数据的第6、8波段。

1 —白垩系塞姆班组; 2—侏罗系安吉拉组; 3—劳瑞莱组; 4—斯平沃组

1 —Cretaceous Sember Formation; 2—Jurassic Anjira Formation; 3—Loralai Formation; 4—Spingwar Formation

选择假彩色合成、主成分分析等方法处理ETM数据，获得321波段（RGB）自然彩色、741波段（RGB）假彩色合成、主成分分析等影像，开展全区1∶25万构造、岩性、矿化带解译; 采用Crosta法处理ETM数据提取蚀变信息，即利用ETM1457波段组合进行主成分分析，根据特征向量矩阵选取相应主成分提取羟基和碳酸根矿物信息，所选取主成分的特征向量中第5波段和第7波段载荷较大且符号相反，利用ETM1345波段组合进行主成分分析提取铁氧化物信息，所选取主成分的特征向量中第1波段和第3波段载荷较大且符号相反（Crosta et al.，1989）。

利用假彩色合成、相对吸收深度比值（Relative absorption band depth，RBD）等方法处理ASTER数据，获得468波段（RGB）假彩色合成、RBD865（RGB）比值假彩色合成影像，进行重点区矿化信息增强（Rowan et al.，2003，2005，2006）; 采用1234、1346和1348波段主成分分析，从3个主成分分析结果中分别选取第1波段和第2波段、第4波段和第6波段、第4波段和第8波段载荷较大且符号相反的主成分，提取铁氧化物、Al-OH和CO₃^2-异常（张玉君等，2009）。

对于资源一号02D影像，以含矿地层像元光谱为参考，利用光谱角制图（Spectral angle mapper，SAM）法进行该影像上重点含矿地层空间分布制图（Kruse et al.，1993）。

2.2.2 遥感解译及遥感找矿标志

基于各种图像处理方法增强的ETM和ASTER数据，以1∶25万地质图为基础，建立各岩性单元、构造的解译标志，解译研究区岩性、构造。根据典型矿床矿化带样品反射光谱特性、影像特征及其与岩性、构造的关系，建立沉积岩有关铅锌矿床遥感找矿标志，并圈定研究区内可能的含矿地层。

3 结果

3.1 沉积岩有关铅锌矿遥感找矿标志

3.1.1 遥感地质解译结果

基于多种图像处理方法处理ETM数据，与1∶25万地质图对比后建立了各岩性单元的解译标志，开展全区范围1∶25万岩性、构造解译。在ETM主成分123（RGB）假彩色合成影像上（图5），各岩性单元界线清晰，不同地层色调、纹理差异明显。其中侏罗纪安吉拉组、劳瑞莱组和斯平沃组呈暗绿色，地层层理清晰，呈条带状，总体以NNW—SSE向展布，岩层走向上存在密集转折，表明岩层褶皱较发育。蛇绿岩位于研究区中部，浅蓝—深蓝色，近N—S向不规则状延伸。

在ETM主成分123影像上，断裂构造特征明显，规模较大断裂多表现为两侧岩性色调、走向、纹理等特征的差异，部分规模较小断层可见错断地层和褶皱构造的影像特征。根据构造解译结果，区内大断裂以近N—S向、NNW—SSE向为主，次级断裂、节理为NE—SW向、近E—W向。断裂构造在区内分布不均，其中侏罗纪地层内发育较多N—S向、NNW—SSE向顺层断裂，与矿化关系密切，而在白垩纪地层内分布有大量NE—SW向、近E—W向的小断裂和节理（图2、图5）。

3.1.2 遥感找矿标志

（1）含矿层位影像特征。根据含矿层位光谱特征，选择3种影像突出含矿层位（图6）。在杜达矿区周边ETM主成分145（RGB）假彩色合成影像上，容矿地层安吉拉组呈紫色、紫红色（图6a）; 在ASTER468（RGB）假彩色合成影像上，容矿地层呈粉红色，走向整体为NNW—SSE向，因此呈NNW—SSE向延伸的条块状、透镜状影像特征，并伴有与走向垂直的NEE—SWW向密集、平行小冲沟，该岩层与东侧岩层接触界线为一条与岩层走向平行的逆冲断层（图6b）; 在ASTER RBD865比值（RGB）假彩色合成影像上，容矿地层呈青色（图6c）。杜达矿床含矿地层所在层位向NNW向延伸约10 km，在此层位有多个已知矿化点，证实了该含矿层位的稳定性及其影像特征的准确性。

1 —铅锌矿床; 2—铅锌矿化点; 3—地层界线; 4—断裂; 5—水体; 6—侏罗系

1 —Pb-Zn deposit; 2—Pb-Zn mineralized point; 3—stratigraphic boundary; 4—fault; 5—water; 6—Jurassic

在ASTER468假彩色合成影像上（图7a），贡尕矿化带呈浅粉色，条带状N—S向延伸，ASTER比值影像上该带以青色显示（图7b），该带在N—S方向延伸约2 km，宽约200 m，向南宽度逐渐变窄。

（2）蚀变特征。根据ETM1345、1457波段组合主成分分析特征向量矩阵，对两个组合主成分分析PC4分别取反后获得铁氧化物信息和羟基-碳酸根矿物信息影像，按μ（中值）+2σ（标准差）阈值分割后叠加于假彩色合成图像（ETM741）之上。ETM蚀变信息提取结果表明，含矿层位具有明显的含羟基-碳酸根蚀变异常，该异常与解译的含矿层位吻合较好（图8a）。

利用ASTER 1234、1346和1348波段主成分分析，对三个组合主成分分析PC4分别取反、阈值分割后提取铁氧化物、Al-OH和CO₃^2-异常。结果表明，容矿地层具有普遍强Al-OH异常，部分层位碳酸盐化（CO₃^2-），而褐铁矿化（Fe³⁺）相对较弱，仅在矿区南部局部地区有所表现（图8b）。

（3）遥感找矿标志。通过对杜达铅锌矿床、贡尕铅锌矿床及该带部分矿点的控矿要素、反射光谱、影像特征的综合分析，确定了胡兹达尔-拉斯贝拉地区与沉积岩有关铅锌矿床遥感找矿标志（表1）。该区的主要找矿标志为:地层为侏罗系安吉拉组、劳瑞莱组和斯平沃组碳酸盐岩和碎屑沉积岩; 影像上NNW—SSE向条块状、透镜状展布; ETM主成分145影像上呈紫色、紫红色，ASTER468影像上呈粉—粉红色，ASTER RBD865比值影像上呈青色; 具有强Al-OH和零星CO₃^2-蚀变异常，局部弱铁染; 褶皱带内、局部发育NNW—SSE向断裂。

3.2 找矿远景区圈定

利用所建立的找矿标志，基于ETM和ASTER影像处理结果，根据控矿要素、影像特征和成矿特征，圈定主要预测含矿层位分布区30处（图9），结果与已知矿化点吻合较好。圈定的含矿层位分布区主要集中在中部和北部，均位于侏罗纪地层内。其中中部区域内的含矿层位明显呈NNW—SSE向展布，与地层走向一致，区域内已知的矿化点主要位于安吉拉组和劳瑞莱组地层内部，部分与斯平沃组地层关系密切，矿化明显顺层展布，呈层状、透镜状; 同时，顺层方向的断裂构造也与所预测的含矿层位关系密切，表明流体的运移和矿体的就位可能受次级断裂控制。

（a）—ETM主成分145影像;（b）—ASTER468影像;（c）—ASTER RBD865影像; 1—侏罗系界线; 2—断裂; 3—铅锌矿床; 4—铅锌矿化点

(a) —Principle component 1, 4 and 5 of ETM in RGB; (b) —ASTER 4, 6 and 8 in RGB; (c) —ASTER RBD 8, 6 and 5 in RGB; 1—Jurassic boundary; 2—fault; 3—Pb-Zn deposit; 4—Pb-Zn mineralized point

（a）—ASTER468影像;（b）—ASTER RBD865影像; 1—地层界线; 2—断裂; 3—侏罗系; 4—预测含矿层位; 5—铅锌矿床; Q—第四系; E2—始新统; K—白垩系; J—侏罗系

(a) —ASTER 4, 6 and 8 in RGB; (b) —ASTER RBD 8, 6 and 5 in RGB; 1—stratigraphic boundary; 2—fault; 3—Jurassic; 4—ore prospecting targets; 5—Pb-Zn deposit; Q—Quaternary; E2—Eocene; K—Cretaceous; J—Jurassic

北部区域内含矿层位呈近N—S向和近E—W向展布，其中贡尕矿床的矿化地层呈近N—S向延伸，而其北西侧矿化带则表现为近E—W向延伸特征，表明在北部地区受构造影响岩层的产状发生了明显的变化，且北部地区的断裂构造的走向也变为近E—W向延伸，对预测的含矿地层具有控制作用。

（a）—ETM影像提取结果;（b）—ASTER提取的Al-OH和碳酸根结果

(a) —Results extracted from ETM image; (b) —results of Al-OH and carbonate extracted from ASTER image

4 主要找矿远景区影像特征分析

4.1 杜达矿床外围远景区

根据图6、图8中杜达矿床矿化地层影像特征及蚀变特征分析可知，安吉拉组含矿地层沿NNW向稳定延伸近12 km，宽约600 m，中部宽两侧窄，该带中部—南部已发现了多处矿化点（由于前期获取的资料比例尺较小，矿化点位置可能存在一定偏差）。而在该层位北部约4 km长的安吉拉组地层尚未发现矿化点（图8b中绿色椭圆范围），表明该区具有较大的找矿潜力。

4.2 贡尕矿床外围远景区

在贡尕矿床东部和南部，有2个预测含矿地层分布区，内部存在影像特征与杜达、贡尕矿床类似的条带或团块，在图7a中呈粉红色，而在图7b比值假彩色合成影像中呈青色，空间分布不规则，以团块状、带状、环状分布，且规模较大。而在这两个预测区内尚未发现矿化，值得开展野外检查工作。

4.3 Barindri Dhora地区矿化带

杜达矿床北北西方向约30 km处的Barindri Dhora地区，有一处预测含矿层位集中区。在ASTER468假彩色合成影像上（图10a），该区预测含矿层位均呈明显的粉红色条带沿NNW—SSE向延伸，ASTER比值影像上这些条带均以青色显示（图10b），且存在多条平行条带，这些条带为斯平沃组（Js），岩性主要为砂岩页岩互层，夹少量灰岩，该地层的水系为NEE—SWW向平行状水系。该区的断裂构造以NNW—SSE向顺层断裂为主，与预测的含矿层位关系密切，部分断裂构成了含矿层位与围岩的界线。在图10中西北角预测区内已发现的多处矿化点均与遥感预测的含矿地层条带吻合较好，而南部和东部2个预测区内有3条与北西侧预测区内影像特征相近的地层，但尚未发现矿化，具有较好的找矿前景。

1 —铅锌矿床; 2—铅锌矿化点; 3—地层界线; 4—断裂; 5—水体; 6—侏罗系; 7—预测含矿地层

1 —Pb-Zn deposit; 2—Pb-Zn mineralized point; 3—stratigraphic boundary; 4—fault; 5—water; 6—Jurassic; 7—prospecting mineralized strata

以表1找矿标志为基础，本研究利用研究区中部地区2019年12月14日成像的资源一号02D高光谱影像探索国产高光谱影像自动提取含矿地层的可行性。根据已知矿化点的位置，从资源一号02D高光谱影像上提取了含矿地层与灰岩的像元光谱（图11），矿化地层像元光谱吸收位置位于2200 nm和2350 nm附近，表明为含Al-OH矿物占主导所致，与图4中安吉拉组样品实测光谱的吸收特征一致; 而灰岩像元光谱在2335 nm存在强吸收，与研究区野外调查获取的斯平沃组和劳瑞莱组中灰岩样品光谱特征类似（图4）。以图11中含矿地层像元光谱为参考，利用光谱角制图（SAM）法处理资源一号02D高光谱影像获得了该景影像上含矿地层的空间分布（图12），结果中含矿地层呈NNW—SSE向带状延伸，与图10中两种影像分析结果吻合较好，图像中部、西部所圈定的含矿地层预测区均被有效提取。

同时，选择Barindri Dhora地区为典型区，对比三种影像提取的含矿地层结果（图13）。其中图13（a）为ETM1457波段主成分分析提取的含羟基-碳酸根结果，图13（b）为ASTER1346波段主成分分析提取的Al-OH结果，而图13（c）为资源一号02D光谱角制图提取结果。三个结果中北西侧重点区内的主要含矿地层均被有效提取，但东侧、南侧2个重点区内含矿地层提取结果则差异较大，ETM、ASTER结果中南侧重点区一处异常无法提取（图13中南侧白色箭头），东侧重点区内异常仅少量像元呈断续状出现（图13中北东侧白色箭头），资源一号02D所提取的含矿地层则更为完整、连续，与图10吻合更好（图13中白色箭头位置）。虽然三种数据所采用的提取方法不同可能会对结果产生部分影响，但由于ASTER结果优于ETM结果，资源一号02D结果优于ASTER的结果，表明提取结果的差异主要是三种影像光谱分辨率不同占主导，证明了资源一号02D高光谱影像在提取本区含铅锌矿地层方面具有更大的优势。如果后续能够获取覆盖整个胡兹达尔-拉斯贝拉地区的资源一号02D影像，根据表1找矿标志中含矿地层光谱存在2200 nm和2350 nm附近吸收特征以及含矿层位等约束条件（表1），可以在该区开展含矿地层自动提取及找矿预测工作。

（a）—ASTER 468影像;（b）—ASTER RBD865影像; 1—地层界线; 2—断裂; 3—预测含矿地层; 4—铅锌矿化点; K—白垩系; J—侏罗系

(a) —ASTER 4, 6 and 8 in RGB; (b) —ASTER RBD8, 6 and 5 in RGB; 1—stratigraphic boundary; 2—fault; 3—ore prospecting targets; 4—Pb-Zn mineralized point; K—Cretaceous; J—Jurassic

1 —矿化地层; 2—灰岩

1 —Mineralized strata; 2—limestone

5 结论

胡兹达尔-拉斯贝拉铅锌成矿带作为特提斯成矿域铅锌成矿带的重要组成部分，该铅锌成矿带及典型矿床的研究程度相对较低，制约了对该成矿带的认识及找矿勘查研究。本研究根据喷流-沉积型（SEDEX）和密西西比河谷型（MVT）型铅锌矿床的成矿地质背景，通过合理选择遥感数据源和图像处理方法，对该区典型矿床及样品光谱特征、影像特征综合分析，结合关键控矿要素研究，建立了胡兹达尔-拉斯贝拉地区与沉积岩有关铅锌矿遥感找矿标志并开展成矿预测，为该成矿带找矿勘查提供可行性方案。

1 —地层界线; 2—断裂; 3—侏罗系; 4—预测含矿地层; 5—提取的矿化地层; 6—铅锌矿化点; Q—第四系; K₂—上白垩统; K₁—下白垩统; J—侏罗系; Σ—蛇绿岩

1 —Stratigraphic boundary; 2—fault; 3—Jurassic; 4—ore prospecting targets; 5—extracted mineralized stratum; 6—Pb-Zn mineralized point; Q—Quaternary; K₂—Upper Cretaceous; K₁—Lower Cretaceous; J—Jurassic; Σ—ophiolites

（1）在基于ETM数据完成了胡兹达尔-拉斯贝拉成矿带1∶25万岩性、构造解译与信息提取工作的基础上，通过对杜达铅锌矿床及该带部分矿点的控矿要素、反射光谱、影像特征的综合分析，建立了胡兹达尔-拉斯贝拉地区与沉积岩有关铅锌矿床遥感找矿标志。该区的主要找矿标志为:地层为侏罗系安吉拉组、劳瑞莱组和斯平沃组碳酸盐岩和碎屑沉积岩; 影像上NNW—SSE向条块状、透镜状展布; ETM主成分145影像上呈紫色、紫红色等，ASTER468影像上呈粉红色，ASTER RBD865比值影像上呈青色; 具有带状Al-OH蚀变异常，局部CO₃^2-和弱铁染; 褶皱带内、局部发育NNW—SSE向断裂。

（2）利用所建立的找矿标志，基于ETM和ASTER影像处理结果，根据控矿要素、影像特征和找矿标志，圈定主要预测含矿地层分布区30处，结果与已知矿床、矿点吻合较好。圈定的含矿地层分布区主要集中在研究区中部和北部，均位于侏罗纪地层内且矿化带与地层走向一致，矿化主要位于安吉拉组和劳瑞莱组地层内部，部分与斯平沃组地层关系密切，矿化明显顺层展布，呈层状、透镜状。

（3）根据含矿地层的像元光谱吸收特征，利用研究区中部地区资源一号02D高光谱影像提取了含矿地层的空间分布，结果与根据找矿标志分析结果吻合较好，资源一号02D高光谱数据提取结果经对比明显优于ETM和ASTER数据提取结果，证实了在该区利用国产高光谱影像开展与沉积岩有关铅锌矿含矿地层自动提取与找矿的可行性。

（a）—ETM结果;（b）—ASTER结果;（c）—资源一号02D结果; 1—地层界线; 2—断裂; 3—预测含矿地层; 4—遥感提取结果; 5—铅锌矿化点; K—白垩系; J—侏罗系

(a) —ETM result; (b) —ASTER result; (c) —ZY-1-02D result; 1—stratigraphic boundary; 2—fault; 3—ore prospecting targets; 4—extracted mineralized strata by remote sensing; 5—Pb-Zn mineralized point; K—Cretaceous strata; J—Jurassic strata

参考文献