论文部分内容阅读
摘要：为了快捷地评价梨树沟地区的沟系土壤化探成果，分别利用具有自主知识产权的两种不同分辨率的国产高分辨影像（GF1、GF2）进行1∶50000和1∶10000的遥感地质解译工作，依据高分影像的解译成果并结合实地地质调查工作，查明了土壤化探异常的成因，缩短了勘查时间，节约勘查经费。
关键词：GF1、GF2高分辨影像；地质特征；研究；地面查证
引言
四川省青川县梨树沟为一预查探矿权，区内植被发育，基岩出露较少，严重影響了地质填图工作的质量。为了尽快明确勘查区内的找矿前景，本次预查工作以1∶25000沟系土壤测量为先导，结合已经取得的区域地质调查资料，对高分一号（GF1）和高分二号（GF2）数据分别进行1∶50000及1∶10000两种不同比例尺的地质解译工作[1]，指导地表槽探工程对1∶25000沟系土壤测量圈定的矿致异常进行解剖，大大缩短了找矿时间。
1. 区域地质及金矿成矿特征
1.1 区域地质及矿产特征
梨树沟金矿区的位于松潘—甘孜造山带（Ⅰ）之丹巴弧形逆冲叠置岩片（Ⅱ）之上[3]，其北部紧邻秦岭—大别山造山带、南东部与扬子地块相接（图1）。根据已经取得区域地质调查资料[2]，研究区地层可分为元古界地层和古生界地层两部分，古生界为寒武系和志留系地层，中间缺失奥陶系地层，除志留系顶部地层外，其余均遭受了低级区域变质活动；元古界地层为泥质—碎屑建造地层夹火山岩建造，在元古界地层上部见变质程度较深的（石英）片岩类地层产出；寒武系地层原岩为一套复理石建造地层；志留系地层在研究区由老到新分别可分为三个建造：浅海相泥质—碎屑岩建造夹碳酸盐建造，该地层为研究区赋矿地层、以火山岩—泥质（碎屑）建造组成的多旋回组合夹碳酸盐建造、滨海—浅海相碎屑（泥质）建造夹碳酸盐建造。
区内主要构造痕迹均呈北东向，为沿不同地层之间产生的大型层间滑脱带以及沿相同地层中不同岩组之间发生层间滑脱带[3][4]，层间滑脱构造控制了本地区金矿的产出状态。
与研究区处于同一大地构造单元的主要金矿床为四川省广元市太阳坪金矿和四川省平武县金洞沟金矿，均处于丹巴俯冲滑脱构造岩片之上（图1、表1）。上述所有金矿（点）床的赋矿围岩均为志留系黄坪组（Shn）
2 遥感构造解译
2.1 遥感数据处理
本次遥感解译软件采用了ENVI5.1。先对遥感数据进行预处理及解译工作，然后导入ARCGIS10.3进行解译图件的编制。本次采用的遥感数据为2016年1月24日生成的高分一号（GF1）数据，ID号为L1A0001364652；2015年12月26日生成高分二号（GF2），ID号为L1A0001278172。GF1数据范围较梨树沟探矿权范围略大，GF2数据范围与沟系土壤测量矿致异常范围一致。解译区无云覆盖，影像清晰，植被发育，无大面积基岩露头，地质解译程度属中等。在进行遥感数据处理之前，使用ENVI5.1分别对三种遥感数据进行了Gram-Schmidt Pan Sharpening融合处理，使GF1图像分辨率达到2m，GF2图像分辨率达到0.8m，满足地质解译精度的要求。
根据1∶50000地形图和在该地区进行水系测量工作时采集的16个坐标数据，对通过融合处理的两种遥感数据图像进行了正射校正。根据相关性越差，冗余性越低，图像信息越丰富的原则，对GF1和GF2图像进行相关性分析（表2、表3、），由于OLI图像的第一波段为海岸线波段，故将其剔除出计算成果，根据计算数据，GF1图像选取4、3、1（R、G、B）三个波段，GF2选取4、3、1（R、G、B）三个波段进行解译工作。
2.2 地层单元的遥感特征
2.3 线—环构造解译
本次使用GF1影像进行探矿权范围内的构造解译工作，解译精度为1∶50000，用GF2影像进行土壤异常区的1∶10000构造解译工作。
2.3.1 GF1影像的地层解译
根据1∶20万区域地质调查资料，区内主要地层为寒武系邱家河组下段（q3）、志留系滑天坡组下段+中段（Sh1+2）、志留系毛塔子组（Sm），由于岩性的差异，其在GF1影像上的特征各不相同（见表5），从影像特征并结合1∶20万区调地质资料分析，解译区的地层均为断层接触。从图像上来看，砂岩、片岩和凝灰质千枚岩类地层由于抗风化能力强容易形成比较陡峭的正地形，而由泥质类岩石变质而成的千枚岩抗风化能力较弱，容易形成比较低矮浑圆的地形。
2.3.2 GF1影像的构造解译
将GF1影像资按1∶25000的比例尺后行解译，经解译后区内主要线性构造方向可分为北东向、近东西向和北西向三组组，解译出的线性构造有25条，其中北东向八条、近东西向七条、北西向十条，正环形构造一处（图1）。北西向线性构造多为后期产出，将近东西向线性构造错切。区内规模较大的线性构造为北东向的F1、F2和近东西向的F12，其它线性构造规模均较小。从线性构造特征来看其倾向多与地层倾向一致，倾角较陡。解译区内规模较大的线性构造或与沟系土壤测量有关的重要线性构造的特征见表6。
2.3.3 GF2影像的构造解译
本次GF2影像解译比例尺为1∶5000，从GF2号解译成果来看，在解译内解译的线性构造格架与该区在GF1中解译基本相似（图2），由于解译比例尺变大，局部小规模的线性构造可以被清晰判读，而且GF2影像解译结果可以对前期GF1影像的解译的线性构造的位置和展布情况进行厘定，如图2中F16。
3. GF2影像的构造与沟系土壤测量的关系
从图2中可以看出，As、Au组合异常其长轴方向与GF2影像解译区内的线性构造走向一直，而且As、Au元素的高（峰）值点多与解译的线性构造吻合或就在其附近，显示了高度的一致性。说明区内的沟系土壤化探异常与线性构造关系密切。
4. 解译结果的查证
本次先对F7和F17相交处的As、Au异常高值区进行地面地质查证，查证发现F7线性构造为一典型的顺层剪切带，中见石英脉沿层间接触面产出，顺层剪切带中的原岩为炭质板岩的软性岩层被剪切作用改造为构造岩，在该地段布设岩石原生晕剖面一条，经分析在剪切带中的构造岩和石英脉的As含量分别218×10-6、201×10-6，Au含量分别为9×10-9、11×10-9，而其两侧的围岩的As含量一般小于4×10-6，Au含量小于5×10-9，说明构造活动为引起上述化探异常的重要因素。
5. 结论
（1）GF1和GF2影像资料为自主知识产权的高分辨率遥感影像，获取方便快捷，价格合理。在本次地质勘查工作中提取的勘查区构造、岩性与沟系土壤化探异常浓集中心高度耦合，对少数线性构造的解译与地面查证结果一致。
（2）GF1、GF2高分辨率影像的遥感解译结合区域地质调查成果用于矿区的地质调查中，能够很好地指导地面地质填图工作，缩短勘查时间，并能比较精确地分析化探异常与构造的关系，指导化探异常的解剖工作，是一种非常经济可行的工作手段。
参考文献：
[1] 遥感技术在境外地质调查中的应用—以津巴布韦大岩墙为例（付长亮等人，国土资源遥感 2015年12月）.
[2] 1∶20万广元幅区域地质调查报告（四川省地质局1966年）.
[3] 松潘甘孜造山带的造山过程——许志琴、侯立伟等（1991年）.
[4] 四川省广元市太阳坪金矿区补充详查地质报告——中国人民武装警察部队黄金第五支队曹东宏、杨登美 2007.3.
[5] 四川省广元市太阳坪金矿区补充详查地质报告——中国人民武装警察部队黄金第二总队曹东宏，岳连雄，李小江等 2012.7.
[6] 邓莉莉， 高丽英. 川西南至滇东北部磷矿区域地质特征及成矿地质条件[J]. 西部资源， 2016（5）.
[7] 付垒， 赵鹏彬， 于占超. ICA在甘肃黄草滩地区铁染蚀变信息提取中的应用研究[J]. 西部资源， 2016（5）.
关键词：GF1、GF2高分辨影像；地质特征；研究；地面查证
引言
四川省青川县梨树沟为一预查探矿权，区内植被发育，基岩出露较少，严重影響了地质填图工作的质量。为了尽快明确勘查区内的找矿前景，本次预查工作以1∶25000沟系土壤测量为先导，结合已经取得的区域地质调查资料，对高分一号（GF1）和高分二号（GF2）数据分别进行1∶50000及1∶10000两种不同比例尺的地质解译工作[1]，指导地表槽探工程对1∶25000沟系土壤测量圈定的矿致异常进行解剖，大大缩短了找矿时间。
1. 区域地质及金矿成矿特征
1.1 区域地质及矿产特征
梨树沟金矿区的位于松潘—甘孜造山带（Ⅰ）之丹巴弧形逆冲叠置岩片（Ⅱ）之上[3]，其北部紧邻秦岭—大别山造山带、南东部与扬子地块相接（图1）。根据已经取得区域地质调查资料[2]，研究区地层可分为元古界地层和古生界地层两部分，古生界为寒武系和志留系地层，中间缺失奥陶系地层，除志留系顶部地层外，其余均遭受了低级区域变质活动；元古界地层为泥质—碎屑建造地层夹火山岩建造，在元古界地层上部见变质程度较深的（石英）片岩类地层产出；寒武系地层原岩为一套复理石建造地层；志留系地层在研究区由老到新分别可分为三个建造：浅海相泥质—碎屑岩建造夹碳酸盐建造，该地层为研究区赋矿地层、以火山岩—泥质（碎屑）建造组成的多旋回组合夹碳酸盐建造、滨海—浅海相碎屑（泥质）建造夹碳酸盐建造。
区内主要构造痕迹均呈北东向，为沿不同地层之间产生的大型层间滑脱带以及沿相同地层中不同岩组之间发生层间滑脱带[3][4]，层间滑脱构造控制了本地区金矿的产出状态。
与研究区处于同一大地构造单元的主要金矿床为四川省广元市太阳坪金矿和四川省平武县金洞沟金矿，均处于丹巴俯冲滑脱构造岩片之上（图1、表1）。上述所有金矿（点）床的赋矿围岩均为志留系黄坪组（Shn）
2 遥感构造解译
2.1 遥感数据处理
本次遥感解译软件采用了ENVI5.1。先对遥感数据进行预处理及解译工作，然后导入ARCGIS10.3进行解译图件的编制。本次采用的遥感数据为2016年1月24日生成的高分一号（GF1）数据，ID号为L1A0001364652；2015年12月26日生成高分二号（GF2），ID号为L1A0001278172。GF1数据范围较梨树沟探矿权范围略大，GF2数据范围与沟系土壤测量矿致异常范围一致。解译区无云覆盖，影像清晰，植被发育，无大面积基岩露头，地质解译程度属中等。在进行遥感数据处理之前，使用ENVI5.1分别对三种遥感数据进行了Gram-Schmidt Pan Sharpening融合处理，使GF1图像分辨率达到2m，GF2图像分辨率达到0.8m，满足地质解译精度的要求。
根据1∶50000地形图和在该地区进行水系测量工作时采集的16个坐标数据，对通过融合处理的两种遥感数据图像进行了正射校正。根据相关性越差，冗余性越低，图像信息越丰富的原则，对GF1和GF2图像进行相关性分析（表2、表3、），由于OLI图像的第一波段为海岸线波段，故将其剔除出计算成果，根据计算数据，GF1图像选取4、3、1（R、G、B）三个波段，GF2选取4、3、1（R、G、B）三个波段进行解译工作。
2.2 地层单元的遥感特征
2.3 线—环构造解译
本次使用GF1影像进行探矿权范围内的构造解译工作，解译精度为1∶50000，用GF2影像进行土壤异常区的1∶10000构造解译工作。
2.3.1 GF1影像的地层解译
根据1∶20万区域地质调查资料，区内主要地层为寒武系邱家河组下段（q3）、志留系滑天坡组下段+中段（Sh1+2）、志留系毛塔子组（Sm），由于岩性的差异，其在GF1影像上的特征各不相同（见表5），从影像特征并结合1∶20万区调地质资料分析，解译区的地层均为断层接触。从图像上来看，砂岩、片岩和凝灰质千枚岩类地层由于抗风化能力强容易形成比较陡峭的正地形，而由泥质类岩石变质而成的千枚岩抗风化能力较弱，容易形成比较低矮浑圆的地形。
2.3.2 GF1影像的构造解译
将GF1影像资按1∶25000的比例尺后行解译，经解译后区内主要线性构造方向可分为北东向、近东西向和北西向三组组，解译出的线性构造有25条，其中北东向八条、近东西向七条、北西向十条，正环形构造一处（图1）。北西向线性构造多为后期产出，将近东西向线性构造错切。区内规模较大的线性构造为北东向的F1、F2和近东西向的F12，其它线性构造规模均较小。从线性构造特征来看其倾向多与地层倾向一致，倾角较陡。解译区内规模较大的线性构造或与沟系土壤测量有关的重要线性构造的特征见表6。
2.3.3 GF2影像的构造解译
本次GF2影像解译比例尺为1∶5000，从GF2号解译成果来看，在解译内解译的线性构造格架与该区在GF1中解译基本相似（图2），由于解译比例尺变大，局部小规模的线性构造可以被清晰判读，而且GF2影像解译结果可以对前期GF1影像的解译的线性构造的位置和展布情况进行厘定，如图2中F16。
3. GF2影像的构造与沟系土壤测量的关系
从图2中可以看出，As、Au组合异常其长轴方向与GF2影像解译区内的线性构造走向一直，而且As、Au元素的高（峰）值点多与解译的线性构造吻合或就在其附近，显示了高度的一致性。说明区内的沟系土壤化探异常与线性构造关系密切。
4. 解译结果的查证
本次先对F7和F17相交处的As、Au异常高值区进行地面地质查证，查证发现F7线性构造为一典型的顺层剪切带，中见石英脉沿层间接触面产出，顺层剪切带中的原岩为炭质板岩的软性岩层被剪切作用改造为构造岩，在该地段布设岩石原生晕剖面一条，经分析在剪切带中的构造岩和石英脉的As含量分别218×10-6、201×10-6，Au含量分别为9×10-9、11×10-9，而其两侧的围岩的As含量一般小于4×10-6，Au含量小于5×10-9，说明构造活动为引起上述化探异常的重要因素。
5. 结论
（1）GF1和GF2影像资料为自主知识产权的高分辨率遥感影像，获取方便快捷，价格合理。在本次地质勘查工作中提取的勘查区构造、岩性与沟系土壤化探异常浓集中心高度耦合，对少数线性构造的解译与地面查证结果一致。
（2）GF1、GF2高分辨率影像的遥感解译结合区域地质调查成果用于矿区的地质调查中，能够很好地指导地面地质填图工作，缩短勘查时间，并能比较精确地分析化探异常与构造的关系，指导化探异常的解剖工作，是一种非常经济可行的工作手段。
参考文献：
[1] 遥感技术在境外地质调查中的应用—以津巴布韦大岩墙为例（付长亮等人，国土资源遥感 2015年12月）.
[2] 1∶20万广元幅区域地质调查报告（四川省地质局1966年）.
[3] 松潘甘孜造山带的造山过程——许志琴、侯立伟等（1991年）.
[4] 四川省广元市太阳坪金矿区补充详查地质报告——中国人民武装警察部队黄金第五支队曹东宏、杨登美 2007.3.
[5] 四川省广元市太阳坪金矿区补充详查地质报告——中国人民武装警察部队黄金第二总队曹东宏，岳连雄，李小江等 2012.7.
[6] 邓莉莉， 高丽英. 川西南至滇东北部磷矿区域地质特征及成矿地质条件[J]. 西部资源， 2016（5）.
[7] 付垒， 赵鹏彬， 于占超. ICA在甘肃黄草滩地区铁染蚀变信息提取中的应用研究[J]. 西部资源， 2016（5）.