[发明专利]基于层次分析法的成矿预测方法

权利要求书

1.一种基于层次分析法的成矿预测方法，其特征在于，所述方法包括以下步骤：

基于遥感影像数据进行地质特征解译，得到地质特征解译结果图；

基于遥感影像数据进行蚀变异常区圈定，得到蚀变异常区圈定结果图；

基于已有的地质资料，通过地球化学水系沉积物测量，进行蚀变带内的地球化学异常区圈定，得到地球化学异常区圈定结果图；

将地质特征解译结果、蚀变异常区圈定结果和地球化学异常区圈定结果作为成矿影响因子，利用层次分析法计算各成矿影响因子的权重；

将地质特征解译结果图、蚀变异常区圈定结果图和地球化学异常区圈定结果图叠加，根据各成矿影响因子的权重，计算叠加图中各像素点的综合权重值，并根据所述综合权重值，圈定成矿预测区。

2.根据权利要求1所述方法，其特征在于，在进行地质特征解译和蚀变异常区圈定前，还包括预处理所述遥感影像数据，具体包括：

地质特征解译前，融合不同时段、不同空间分辨率的遥感影像数据；对融合后的遥感影像进行拉伸处理，增加图像的可视性；

遥感蚀变异常区圈定前，统一遥感影像数据中各波段的空间分辨率，并将各波段合并；对统一空间分辨率后的各波段进行大气校正；对大气校正后的各波段进行干扰地物剔除。

3.根据权利要求2所述方法，其特征在于，所述基于遥感影像数据进行地质特征解译，具体包括：在遥感影像上分别建立断裂、褶皱、环形构造和岩体接触带的解译标识，结合已有的地质图和地形图资料，通过目视解译对解译标识进行归纳整理，得到地质特征解译结果图。

4.根据权利要求3所述方法，其特征在于，所述基于遥感影像数据进行蚀变异常区圈定，具体为：选用ASTER遥感影像数据作为所述遥感影像数据，通过主成分分析法进行蚀变异常区圈定，具体包括：

基于蚀变带内的主要蚀变矿物，从ASTER 1至ASTER 9波段中选取的四个波段及四个波段分别对应的各像素点，建立协方差矩阵；

求出所述协方差矩阵的特征值和特征向量；利用所述特征向量作为行，求出变换矩阵；

利用所述变换矩阵，再对每个像素点进行变换，求得图像的四个主分量PC1、PC2、PC3和PC4；所述每个主分量均为四维向量；

根据蚀变矿物的波谱特征，从四个主分量中选定与所述蚀变矿物相对应的主分量；

统计所述选定的主分量的像素灰度平均值和标准离差；

根据所述像素灰度的平均值和标准离差，确定蚀变矿物的一级、二级和三级蚀变异常区的阈值，根据所述阈值分别圈定蚀变矿物的一级、二级和三级蚀变异常区，得到蚀变异常区圈定结果图。

5.根据权利要求4所述方法，其特征在于，在基于已有的地质资料，通过地球化学水系沉积物测量，进行蚀变带内的地球化学异常区圈定完成之后，还包括步骤：基于所述地球化学异常区圈定的结果，通过比较各地球化学异常区的异常信息显著程度，将多个地球化学异常区分为一级、二级和三级地球化学异常区。

6.根据权利要求5所述方法，其特征在于，所述将地质特征解译结果、蚀变信息提取结果和地球化学异常圈定结果作为成矿影响因子，利用层次分析法计算各成矿影响因子的权重，具体包括：

建立成矿影响因子的层次结构模型；

根据所述层次结构模型建立判断矩阵P；

根据所述判断矩阵P计算各成矿影响因子的权重。

7.根据权利要求6所述方法，其特征在于，所述建立成矿影响因子的层次结构模型，具体包括：

将地质特征解译结果、蚀变信息提取结果和地球化学异常圈定结果作为一级成矿影响因子；

将断裂、褶皱、环形构造、岩体接触带、一级蚀变异常区、二级蚀变异常区、三级蚀变异常区、一级地球化学异常区、二级地球化学异常区和三级地球化学异常区作为二级成矿影响因子；

将断裂、褶皱、环形构造和岩体接触带分别对应的0～0.3km缓冲区、0.3～0.6km缓冲区和0.6～1.0km缓冲区作为三级成矿影响因子。