层次分析法在冀北地区多金属矿找矿预测中的应用

范素英

(河北省遥感中心，石家庄 050021)

层次分析法在冀北地区多金属矿找矿预测中的应用

范素英

(河北省遥感中心，石家庄 050021)

冀北地区大地构造单元包括华北陆块、天山-兴蒙造山系及中国东部造山-裂谷系，区内太古代变质岩广泛分布，海西期及燕山期岩浆侵入、喷发活动频繁，围岩蚀变发育，是寻找多金属矿产的有利区域。在分析遥感地质解译要素、蚀变遥感异常与矿化关系的基础上，以中分遥感解译的线、环及蚀变遥感异常、侵入岩、赋矿地层及矿产分布信息为判断因子，运用层次分析法建立找矿预测模型，圈出找矿靶区29个，其中3个找矿靶区内已发现5处大型-特大型矿床，6个找矿靶区内发现12处中型矿床。结果表明，在冀北地区运用层次分析法进行多金属矿找矿预测及找矿靶区圈定，具有较好的实用性。

层次分析法； 找矿预测； 冀北

0 引言

国内外学者先后利用层次分析法[1]在不同地区开展了找矿预测研究。肖巧艳等[2]在应用层次分析法集成多元信息的基础上, 应用Kriging +Natural Breaks和C-A分形法圈定了找矿靶区; 周廷全等[3]以烃生成、储集、圈闭及运移等条件为判断因子，对楚雄盆地有利油气勘探区带进行了评价; 刘超等[4]建立了综合信息方法的层次结构模型。王永军等[5-6]把层次分析法运用于张家口地区金矿、铀矿成矿预测。

冀北地区矿产受华北陆块北缘东段成矿带和太行山成矿带控制，目前该区已发现蔡家营银铅锌矿、东坪金矿、撒岱沟门钼矿、寿王坟铜钼矿、北岔沟门铅锌矿、峪耳崖金矿等大型矿床。本文以中分遥感地质解译信息、蚀变遥感异常信息及地质、矿产资料为判断因子，在冀北地区应用层次分析方法进行多金属矿找矿预测并圈定找矿靶区，以期科学评估冀北矿产资源潜力，为进一步探索该区找矿潜力及科学部署矿产资源勘查工作提供依据。

1 研究区概况及数据源

1.1 研究区概况

冀北地区位于河北省北部，包括承德、张家口2个地级市，地理坐标跨E115°00′～117°50′，N40°35′～42°00′。冀北Ⅰ级大地构造单元为天山-兴蒙造山系、华北陆块区和中国东部造山-裂谷系。具活动带性质的天山-兴蒙造山系位于康保―围场深断裂以北，质变主旋回在海西期； 康保―围场深断裂以南为相对稳定的华北陆块区，古元古代末形成结晶基底，上覆平缓的沉积盖层，可进一步分为阴山-冀北地块和晋冀地块； 中国东部造山-裂谷系发生在中-新生代，以强烈的断裂构造、岩浆活动叠加在早期构造行迹之上(图1)。

图1 构造纲要图

天山-兴蒙造山系仅出露二叠纪沉积岩-火山岩。华北陆块区地层出露较全，主要有前寒武纪变质岩、中新元古代―古生代沉积岩、中生代火山岩―碎屑岩等三套岩系。前震旦纪变质岩主要为片麻岩、变粒岩、斜长角闪岩夹磁铁石英岩为主的变质岩建造，变质程度麻粒岩―角闪岩相，次为古元古代片岩、砂岩、变质火山岩、大理岩等浅变质岩建造，变质程度绿片岩相； 中元古代―新元古代―古生代沉积岩为盖层组合，自下而上划分为长城纪滨海相富铁碎屑岩―富镁碳酸盐岩建造、蓟县纪滨海相―浅海相碳酸盐岩建造、青白口纪滨海相碳酸盐岩―碎屑岩建造、寒武―奥陶纪浅海相碎屑岩―碳酸盐岩建造、石炭纪海陆交互相含煤建造和二叠纪陆相碎屑岩建造，缺失晚奥陶纪―早石炭纪沉积。中国东部造山系主要表现为断裂构造岩浆活动和火山岩-沉积岩断陷盆地2个方面。中生代岩浆侵入活动受断裂控制，形成NE向构造岩浆带和EW向构造岩浆带； 中生代火山-沉积岩基本以尚义、赤城、隆化断裂为界，断裂以北为冀北陆相火山沉积断陷带，断裂以南为燕山陆相火山沉积凹陷带。

冀北地区矿藏丰富，矿产受华北陆块北缘东段成矿带和太行山成矿带控制。华北陆块北缘东段成矿带是纬向成矿带，岩浆活动强烈，主要矿产有铁、铜、铅锌、金、银、钼矿等。太行山成矿带走向NNE，是一个中生代岩浆作用的成矿带。本区位于太行山成矿带的北段，主要矿产有铀、钼、铅锌、银、金矿等，矿床成因与中生代侵入岩或火山岩有关。

冀北地区矿产主要成矿方式为热液型矿床、岩浆型矿床、沉积-变质型铁矿和沉积型矿床。热液型矿床和岩浆型矿床围岩蚀变较强，主要蚀变有硅化、绢云母化、蛇纹石化、碳酸盐化、绿泥石化、青磐岩化以及萤石化、叶腊石、褐铁矿化、黄钾铁矾化等； 沉积-变质铁矿、沉积矿床没有围岩蚀变。

1.2 数据源及其预处理

本研究所用遥感数据为Landsat7的ETM+数据，景内无云、雪覆盖，影像层次感强，纹理清楚，满足遥感信息提取的要求。

遥感数据预处理包括大气校正、几何纠正、假彩色合成和数据融合等手段，其目的是提高图像的可解性，确保成果的精确度。大气校正是利用PCI软件Flaash模块进行，用以消除或减少大气、光照影响，还原地物的真实的地表反射率、辐射率和地表温度等参数。几何纠正控制点主要来源于1∶5万地形图，采用多项式法拟合，校正结果与地形资料吻合情况较好，满足遥感制图精度要求。假彩色合成波段组合采用ETM B7，B4，B2组合，合成影像具有较好的反差和饱和度，能够体现岩性和构造信息。数据融合采用PAN方法，改善了图像分辨率，有较好的目视效果，有利于目视解译。

2 遥感信息提取

2.1 遥感地质要素分析与成矿模型构建

遥感地质要素是指ETM图像上显示的与目标地质体光谱特征相关的线、环、带影像和遥感异常信息[7-8]。本研究在分析已有区域地质矿产资料和已知典型矿床特征的基础上，通过人机交互的图像解译和数字处理，提取与矿带、矿田、矿化蚀变或矿体赋存相关的遥感地质要素(图2)和蚀变异常信息，同时结合对地质、矿产和物化探资料的综合研究，建立遥感找矿模型，进行找矿预测，圈定了成矿有利区带或找矿远景。

1)线要素是指与成矿、控矿、导矿和容矿相关的线性构造信息，如断裂构造、脆―韧性剪切强变形带等。冀北地区构造格局以EW和NNE断裂为主，NW向次之。EW向断裂是新元古代一古生代构造运动的产物，切割深度大，由北而南依次为康保―围场断裂、丰宁―隆化断裂、大庙断裂、尚义―赤城-平泉断裂、下板城断裂和兴隆―青龙断裂。EW向断裂早期伴有基性-超基性岩及中酸性岩浆活动，矿床以岩浆型铁、磷、钒、钛、铂、铬、镍矿为主； 中生代时期活动强烈，控制了侏罗纪-白垩纪岩浆活动及与之有成因联系的矽卡岩型和斑岩型铁、铜、金、铅锌、银及钼等金属矿产。NNE向断裂主要受滨太平洋构造域影响，断裂性质以压扭性居多，规模较大，同时强烈的造山运动诱发大规模中、酸性岩浆活动，沿断裂有大量中生代花岗岩侵入和火山喷发，形成了铁、铜、铅锌、钼矿等多金属矿产。NW向断裂不如前两者发育，较重要的有张家口断裂、松枝口断裂等。

2)环要素是指与岩浆侵入、火山机构或构造等相关的圆形-椭圆形影像信息。按环形构造的成因，一般分为侵入岩类环形构造、火山机构类环形构造、构造类环形构造和成因不明环形构造4种类型。研究区内以具有岩浆成因的侵入岩类环形构造、火山机构环形构造为主，与区域成矿作用关系密切。其中侵入岩类环形构造包括花岗岩、闪长岩、基性岩和隐伏岩体环形构造等类型，与之有关的矿床类型则有岩浆型、斑岩型、矽卡岩型、热液型铁、铜、金、铅锌、银、钼、钒、钛、铂、铬、镍及磷矿等金属、非金属矿床，如岩浆型矾山磷铁矿、斑岩型―矽卡岩型寿王坟铜钼矿、岩浆热液型牛圈银矿和蔡家营铅锌银矿等。火山机构类环形构造包括火山口、火山通道环形构造和浅成、超浅成次火山岩环形构造等类型，分布于冀北火山岩盆地内，成矿类型为火山热液型铅锌、银、锰及萤石矿，如小扣花营铅锌银矿、相广锰矿及郝家楼萤石矿等。

1.前寒武变质岩； 2.长城纪高于庄组； 3.蓟县纪雾迷山组； 4.中生代火山岩； 5.基性-超基性岩；6.元古宙花岗岩； 7.晚古生代花岗岩； 8.中生代花岗岩； 9.大型断裂； 10.中型断裂；11. 侵入岩环形构造； 12. 火山机构环形构造； 13. 成因不明环形构造

3)带要素是指与赋矿岩层、矿源层相关的面状地层单元信息。冀北地区赋矿地层主要有前寒武纪变质岩、长城纪高于庄组白云岩、蓟县纪雾迷山组白云岩、中生代火山岩等。前寒武纪变质岩是沉积变质铁矿、金矿赋矿围岩； 中生代火山岩是萤石、珍珠岩、沸石等非金属矿和铅锌、银、锰等多金属矿赋矿围岩； 矽卡岩类矿床产于花岗岩与长城纪高于庄组白云岩、蓟县纪雾迷山组白云岩等碳酸盐类岩石接触带部位，以铁、铜、铅锌、银及钼矿为主。

2.2 蚀变信息提取

矿化蚀变信息提取的方法较多，如主成分分析法、比值法或光谱角法等。本研究采用主成分分析法提取了羟基、铁染遥感异常信息。以B1，B4，B5，B7这4个波段的第4主分量提取羟基异常； 以B1，B3，B4，B5这4个波段的第4主分量提取铁染异常。处理中对数据做了必要的去干扰处理，较大程度上减少了环境对异常提取的影响。

羟基异常与高岭土化、绢云母化、绿泥石化、碳酸盐岩化等含羟基或碳酸根基团的蚀变矿物有关； 铁染异常与铁氧化物的蚀变矿物有关。异常受EW向构造和NNE向构造控制，形成了具相同走向的2组异常带。这2组异常带在区域上构成斜十字形形态，其中心在赤城附近(图3)。其中EW向异常带以铁染组合为主，部分地段为羟基组合； NNE向异常带分布于蔚县、宣化、赤城、丰宁上黄旗一线，主要受乌龙沟-上黄旗断裂控制，异常类型以羟基为主，局部为铁染组合。

图3 研究区羟基和铁染组合遥感异常分布图

3 找矿远景预测与分析

3.1 找矿远景预测模型

层次分析法原理简单，有较严格的数学依据，已被广泛应用于复杂系统的决策、评价、分析和预测等工程。本研究采用层次分析法进行找矿远景预测，选取了与成矿关系密切的遥感信息或地学信息做为预测因子，将其按影响程度分级，构造判断矩阵，求解预测因子权重，构建网格，计算每个网格单元的预测指数。根据预测模型进行评价分析，预测成矿有利地段，圈定找矿靶区。

3.1.1 建立层次结构模型

遵循科学性、可度量性、可比较性及可操作性的原则性，筛选预测因子，建立预测因子指标体系，形成层次结构模型(图4)。模型中以找矿远景预测(A)为总目标层； 以矿(化)点(B1)、侵入岩(B2)、断裂构造(B3)、环形构造(B4)、赋矿地层(B5)、遥感异常(B6)等预测因子为准则层，其中B1为地学信息预测因子，B2—B6为遥感信息预测因子； 以各预测因子的指标为指标层。

图4 找矿靶区评价指标体系层次结构图

3.1.2 确定预测因子分级

按照各因子对成矿作用的贡献进行分级，并赋以分值。

3.1.2.1 矿(化)点(B1)

0处=0,1处=1,2处=3,多处=5。

3.1.2.2 侵入岩(B2)

0处=0,1处=4。

3.1.2.3 断裂(B3)

0条=0,1条=2,2条=4,多条=6。

3.1.2.4 环形构造(B4)

花岗岩类环形构造0个=0，1个=1，多个=3； 火山机构环形构造0个=0，1个=1，多个=3； 构造类、性质不明环形构造0个=0，多于1个=1。

3.1.2.5 赋矿地层(B5)

0处=0，变质岩=碳酸盐=火山岩=4。

3.1.2.6 异常(B6)

0处=0，羟基异常=铁染异常=4。

3.1.3 构建判断矩阵

判断矩阵是表示本层所有因子对上一层某一个因子相对重要性的比较，即从矩阵的第二层开始运用9标度对从属于上一层中每个要素的同层各要素间进行两两比较，如模型中的要素i相对于要素j对上层要素的重要程度，1表示i与j同等重要，3表示i比j略重要，5表示i比j重要，7表示i比j重要很多，9表示i比j极其重要，可以用Bi/Bj表示该重要程度，评定预测因子相对重要性，两两比较后可以得到矩阵，如表1所示。

表1 判断矩阵

3.1.4 判断矩阵特征向量和最大特征根计算

层次分析法确定权重系数的问题，可归结为判断矩阵的特征向量和最大特征根计算问题。

3.1.4.1 特征向量的计算

运用方根法计算特征向量W的分量Wi

1)计算判断矩阵B每一行元素的乘积Mi，即

(1)

式中：n为判断矩阵阶数；bij为判断矩阵第i行第j列元素。

2)计算Mi的n次方根，即

(2)

(3)

式中：Wi为判断矩阵特征向量，即为近似的权向量。则W=[W1，W2，…，Wn]T即为所求特征向量。

3.1.4.2 最大特征根λmax的计算

计算与特征向量W对应的最大特征根λmax，即

(4)

式中： (BW)i为向量BW的第i个分量；n为判断矩阵的阶数。

3.1.5 一致性检验

判断矩阵是否具有满意的一致性，直接影响到由判断矩阵得到的权向量是否真实地反映各比较因素之间的客观权重。因此，需要对判断矩阵进行一致性检验。

1)计算一致性指标CI，即

CI=(λmax-n)/(n-1)，

(5)

式中：λmax为判断矩阵最大特征值；n为判断矩阵的阶数。

2)计算一致性比率CR。一致性指标CI与平均一致性指标RI之比，即

CR=CI/RI。

(6)

3)平均随机一致性指标RI是多次重复进行随机判断矩阵特征值的计算后取算术平均值得到的，对于1—6阶矩阵，查表，取值如表2所示。

表2 平均随机一致性指标RI

对于二阶以上的判断矩阵，当满足判断矩阵的一致性比率CR<0.10时，即可认为判断矩阵具有令人满意的一致性，由此得出的权重值是合理的； 否则应对判断矩阵作适当修正，再重新计算矩阵的权向量，直至一致性检验通过为止。

3.1.6 建立找矿远景预测模型

引入成矿预测指数E进行找矿远景预测评价，成矿预测指数E定义为单位网格面积内各种预测因子与其对应权值乘积的总和，可用以下数学模型描述，即

(7)

式中：E为成矿预测指数；Ri第i个预测因子；Wi第i个预测因子权重值。

3.2 找矿远景预测

3.2.1 最大特征根及所对应的特征向量的计算

以遥感异常、赋矿地层、环形构造、断层、侵入岩、矿(化)点等遥感信息和地学信息为预测指标，列为横行和纵列，运用9标度法，对每个因子进行两两比较后构造判断矩阵(表3)。

表3 靶区评价因子判断矩阵①

①λmax=6.281 1；CI=0.056 22；RI=1.24；CR=0.045<0.1 满足一致性检验。

求解判断矩阵的最大特征根λmax及所对应的特征向量(表4)，特征向量就是预测因子的权值。

表4 预测因子权值

3.2.2 网格的构建

对研究区进行预测评价、圈定找矿远景需要构建合适的网格。研究区面积约76 300 km2，工作比例尺为1∶50万，将网格单元设置为4 km×4 km。

3.2.3 找矿远景预测

运用GIS空间分析功能，通过对各个地质环境因子现状图的空间叠加、计算，按数学模型计算出成矿预测指数E。经统计，E分布在0.2～10.5之间，再利用GIS空间分析模块，求得综合指数E频率拟合曲线，根据曲线突变点(E分别为0.8，2.8，4.2)，确定冀北地区多金属矿找矿预测阈值为4.2。当E≥4.2时定为成矿有利区；E<4.2为次此要或无成矿区。

3.3 找矿远景预测分析

成矿有利区成矿条件较好，岩浆岩、赋矿地层、断裂构造、环形构造、异常、矿点等遥感或地学成矿要素多于3项。一般成矿区成矿条件一般，遥感或地学成矿要素不会多于3项。用包络线将成矿有利区圈出来，即为遥感找矿靶区。冀北地区共圈出29个遥感找矿靶区(图5)。

根据资料，研究区已知大型-特大型矿床5处，中型矿床12处，这些矿田或矿产地都能在相应的找矿靶区找到对应的位置，与圈定的遥感找矿靶区高度吻合。蔡家营靶区内有蔡家营大型银铅锌矿床、撒岱沟门靶区内有撒岱沟门大型钼矿床、小营盘靶区内有东坪金矿和下双台金矿及小营盘金矿大型床； 峪耳崖靶区内有峪耳崖中型金矿床、寿王坟靶区内有寿王坟中型铜钼矿床、牌楼靶区内有北岔沟门中型铅锌矿床、姑子沟矿靶区内有温家沟中型金矿床、小张家口靶区内有黄土梁金矿和后沟金矿中型矿床、小营盘靶区内有于沟岩金矿和水泉沟金矿及南冷沟金矿中型矿床。以与成矿关系密切的遥感要素叠加地学成矿要素为判断因子，采用层次分析法圈定的遥感找矿靶区效果显著，定性、定位可靠，作为找矿远景预测具有较高的实用性。

1．棋盘山找矿靶区； 2．牌楼找矿靶区； 3．牛圈子找矿靶区； 4．撒岱沟门找矿靶区； 5．北头营找矿靶区； 6．八达营找矿靶区；7．姑子沟找矿靶区； 8．胡麻营找矿靶区； 9．大庙找矿靶区； 10．杨树岭找矿靶区； 11．小寺沟找矿靶区； 12．山神庙找矿靶区； 13．挬罗台找矿靶区； 14．寿王坟找矿靶区； 15．蘑菇峪找矿靶区； 16．峪耳崖找矿靶区；17．挂兰峪找矿靶区； 18．炭头山找矿靶区；19．蔡家营找矿靶区； 20．美人沟找矿靶区； 21．青羊沟找矿靶区； 22．小张家口找矿靶区； 23．小营盘找矿靶区； 24．庞家堡找矿靶区； 25．八宝山找矿靶区；26．杏园找矿靶区； 27．口前找矿靶区； 28．矾山找矿靶区； 29．坊口找矿靶区

4 结论

本研究以Landsat7 ETM+数据为信息源，对冀北进行遥感地质解译，建立解译标志，提取线、环、带等遥感要素和羟基、铁染异常，获取了与成矿作用密切相关的控矿构造、岩浆信息、赋矿地层以及羟基、铁染信息。通过分析冀北地区铁、锰、铅锌、银、铜、钼、金矿多金属矿成矿、控矿地质条件，总结遥感要素与铁、锰、铅锌、银、铜、钼、金矿多金属矿化关系，尝试建立以与成矿关系密切的遥感线、环、带要素和遥感异常为遥感判断因子，结合地质、矿产地等地学判断因子，运用层次分析法进行找矿远景预测。

层次分析法原理简单，有较严格的数学依据，广泛应用于复杂系统的决策、评价、分析、预测等工程。本研究采用层次分析法进行找矿远景预测，是将线、环、赋矿岩层、遥感蚀变等遥感信息和地质矿(化)点等地学信息作为判断因子，根据专家意见和分析者的判断，将每个层次判断因子的重要性进行两两比较，而后，利用数学方法计算反映每一层次判断因子的相对重要性次序的权值，根据预测模型进行评价分析，预测成矿有利地段，在冀北地区共圈出29个遥感找矿靶区。通过将圈定的遥感找矿靶区与已知矿产地质资料对比，一些已知的矿田或矿产地与圈定的遥感找矿靶区高度吻合。这表明，采用层次分析法圈定的找矿靶区方法简便，效果显著，在区域性找矿远景预测中具有实用意义和推广价值。

[1] 徐建华.现代地理学中的数学方法[M].北京:高等教育出版社,1996,224-230. XU J H.Mathematical Methods in Contemporary Geography[M].Beijing:Higher Education Press,1996,224-230.

[2] 肖巧艳,王功文,张寿庭,等.豫西南杜关-云阳钼多金属成矿预测研究[J].现代地质,2011,25(1):94-100. Xiao Q Y,Wang G W,Zhang S T,et al.Research on molybdenum polymetallic metallogenic prognosis in Duguan-Yunyang district,southwest Henan Province[J].Geoscience,2011,25(1):94-100.

[3] 周廷全,胡建武.基于层次分析法的楚雄盆地油气地质异常评价[J].高校地质学报,2008,14(2):237-242. Zhou T Q,Hu J W.Evaluation of oil-gas geological anomalies based on the hierarchical analysis for Chuxiong basin, Yunnan Province[J].Geological Journal of China Universities,2008,14(2):237-242.

[4] 刘超,王于天,陈爱菊.层次分析法在综合信息矿产资源预测中的应用[J].长春地质学院学报,1994,24(2):222-228. Liu C,Wang Y T,Cheng A J.Application of AHP in the predication of mineral resources by comprehensive information[J].Journal of Changchun University of Earth Sciences,1994,24(2):222-228.

[5] 王永军,李名松,全旭东,等.基于GIS的层次分析法在张家口北部地区金矿成矿预测中的应用[J].地质科技情报,2007,26(4):7-10. Wang Y J,Li M S,Quan X D,et al.Application of GIS-based analytic hierarchy process for minerogenetic prediction in northern Zhangjiakou region[J].Geological Science and Technology Information,2007,26(4):7-10.

[6] 卢辉雄,王永军,汪冰,等.基于GIS层次分析法在沽源地区铀成矿预测中的应用[J].地球科学进展,2014,29(8):968-973. Lu H X,Wang Y J,Wang B,et al.Application of gis-based analytic hierarchy process for uranium minerogenetic prediction in Guyuan Region[J].Advances in Earth Science,2014,29(8):968-973.

[7] 于学正,曾朝铭,燕云鹏,等.遥感资料应用技术要求[M].北京:地质出版社,2010. Yu X Z,Zeng C M,Yan Y P,et al.Application Technique Requirement of Remote Sensing Data[M].Beijing:Geological Publishing House,2010.

[8] 范素英,郑国庆.河北省矿产资源潜力评价遥感资料应用研究[M].武汉:中国地质大学出版社,2015. Fan S Y,Zheng G Q.The Study of Mineral Resources Potential Evaluation Based on Remote Sensing Data in Hebei Province[M].Wuhan:China University of Geosciences Publishing House,2015.

(责任编辑： 邢宇)

Application of analytic hierarchy process method to ore-prospecting prognosis in northern Hebei

Fan Suying

(CenterofHebeiRemoteSensing,Shijiazhuang050021，China)

North Hebei geotectonic unit includes North China craton, Tianshan - Xingmeng orogenic system and China's eastern orogenic mountains-rift system. Archean metamorphic rocks are widely distributed, the magmatic intrusion-eruption activities of Hercynian and Yanshanian period are frequent and, what is more, wallrock alterations are widely spread in this area; therefore, polymetallic deposits are likely to be found in this area. According to the relationship between the remote sensing geological interpretation factors, alteration remote sensing anomalies and the mineralization, the author selected the linear and ringed structures in medium space resolution remote sensing images and alteration remote sensing anomalies, intrusive rocks, ore formation and mineral distribution information as judgment factors, and established the prospecting prognositic models by AHP (Analysis Hierarchical Process). As a result, 29 prospecting target areas were delineated, 5 superlarge ore deposits were found in 3 prospecting target areas, and 12 medium-sized ore deposits were found in 6 prospecting target areas. The results indicate that polymetallic ore prospecting prediction and delineation of prospecting targets can achieve good effect by AHP in northern Hebei area.

analytic hierarchy process; ore-prospecting prognosis； northern Hebei

10.6046/gtzyyg.2017.02.18

范素英. 层次分析法在冀北地区多金属矿找矿预测中的应用[J].国土资源遥感,2017,29(2)：125-131.(Fan S Y.Application of analytic hierarchy process method to ore-prospecting prognosis in northern Hebei[J].Remote Sensing for Land and Resources,2017,29(2):125-131.)

2015-12-10；

2015-12-20

中国地质调查局地质调查项目“河北省矿产资源潜力评价”(编号： 1212011121001)资助。

范素英(1966-)，女，教授级高工，主要从事遥感测绘和遥感地质研究。Email： fansuying2012@126.com。

TP 79

A

1001-070X(2017)02-0125-07