文章编号：1004-0609(2012)03-0880-09

铜陵天马山及其外围地区矿床空间信息找矿模型

疏志明^{1, 2}，杨  斌^{1, 2}，王雄军^{1, 2}，叶  珂^{1, 2}，梁恩云^{1, 2}，刘庚寅^{1, 2}

(1. 中南大学 有色金属成矿预测教育部重点实验室，长沙 410083；

2. 中南大学 地球科学与信息物理学院，长沙 410083)

摘  要：基于GIS技术，对铜陵天马山及其外围约100 km²区域进行数字矿床空间信息找矿预测模型的研究。通过工作区信息统计单元的划分、预测区地质信息(包括构造、地层、岩浆岩、地表矿化以及遥感混合蚀变信息等)变量的确定以及编码和赋值，采用特征分析法确定空间网格单元成矿异常有利度模型，根据预测单元计算结果，结合地质分析，圈定找矿靶区14处。

关键词：找矿预测；信息统计；特征分析；天马山及其外围地区；铜陵

中图分类号：P624　　     文献标志码：A

Ore prediction model for spatial information of deposits in Tianmashan and its periphery area, Tongling

SHU Zhi-ming^{1, 2}, YANG Bin^{1, 2}, WANG Xiong-jun^{1, 2}, YE Ke^{1, 2}, LIANG En-yun^{1, 2}, LIU Geng-yin^{1, 2}

(1. Key Laboratory of Metallogenic Prediction of Nonferrous Metals, Ministry of Education,

Central South University, Changsha 410083, China;

2. School of Geosciences and Info-Physics, Central South University, Changsha 410083, China)

Abstract: Based on geographic information system (GIS), the model of space information related with ore prediction in Tianmashan and its periphery area, Tongling of 100 km² was studied. According to the division of information statistics units, the determination of geological information variables which consist of structure, strata, magmatic rocks, mineralization and remote sensing alteration information etc., along with coding and assignment, the method of signature analysis was adopted to ensure the mineralized anomalies favorability model of the space information grid units. Combined with geological analysis, 14 prospecting targets were delineated.

Key words: ore prediction; information statistics; signature analysis; Tianmashan and its periphery area; Tongling

矿产资源定量预测是在对矿床与地质条件之间的关系进行分析的基础上实施的，借用地质单元将矿床与地质信息联系起来建立预测模型，得以对已知和潜在矿床及其相应的资源量作出预测及评价^[1]。20世纪80年代以来，赵鹏大和陈永清^[2]认为“在传统的矿床统计分析方法中，须用样本的观测结果来描述总体特征和确定远景区”。因此，定量预测首先应保证抽样  的随机性，还要保证样品的代表性。为此，通常选择一定大小的网格将整个研究区划分为面积相等、形状相同的“单元”。由于“单元”犹如地质取样中的样  品，用作统一预测和取值范围的基本单位，同时也是进行成矿远景计算、比较、评价的基本单位，因此，单元的大小和形状对预测效果有很大的影响^[3-7]。本 文作者根据空间信息单元化定位预测方法，采用特征分析方法开展铜陵天马山及其外围地区成矿条件及找矿预测研究，对地质、地球物理、地球化学、遥感等多源信息进行有机综合，建立矿床空间信息模型，为找矿靶区优选提供依据。

1  矿区地质概况

本文研究范围包括天马山矿区及其外围约100 km²区域(见图1)。区内主要出露志留系、泥盆系、石炭系、二叠系、三叠系及第四系地层，其中石炭系、二叠系、三叠系为主要赋矿层位。

区内褶皱和断裂构造发育，主要褶皱有铜官山背斜和金口岭向斜等。铜官山背斜斜贯测区中部，长约18 km，背斜轴线总体方向为42°。背斜向北东倾伏，倾伏角约10°。其核部为志留系地层，向两翼依次出露志留系至三叠系地层。背斜倾伏端的北西翼地层倾向北西，倾角30°~50°。南东翼地层浅部倾向南东，倾角约25°，至深部发生倒转，倾向北西，倾角在60°以上。金口岭向斜位于铜官山背斜之北西，长约11 km。褶皱线的方向从南西端的35°向北东逐渐转成45°。向斜核部为三叠系中统石灰岩。区内发育北东向和北西向断裂，主要断裂有白鹤—松树山走向断层、笔西走向断层、宝山断层、老庙基山断层和松树山—尾砂坝断层等。

图1  天马山及其外围地区信息统计单元划分图

Fig. 1  Subdivision map of information statistics units for Tianmashan and its periphery area: 1—Quaternary system; 2—Chishan group; 3—Longtoushan group; 4—Fenshuiling group; 5—Nanlinghu group; 6—Tashan group; 7—Xiaoliangting group; 8—Dalong group; 9—Longtan group; 10—Maokou group; 11—Qixia group; 12—Middle and upper Carboniferous; 13—Wutong group; 14—Zhaishan group; 15—Fentoushan group; 16—Fracture; 17—Conjectural fracture; 18—Fold axis; 19—Small-sized flexure; 20—E-W rise; 21—Diorite; 22—Quartz diorite; 23—Diorite porphyrite; 24—Gossan; 25—Ore occurrence; 26—Scope of blending remote sensing alteration; 27—Grids of statistical forecast unit

区内发育燕山期侵入岩，主要岩体有天鹅抱蛋山岩体、铜官山岩体和金口岭岩体等。天鹅抱蛋山岩体侵入于铜官山背斜倾伏端东侧，出露于马山、青石山、天鹅抱蛋山之间，岩性为闪长岩体，出露形态呈不规则圆形，面积0.7~0.8 km²。铜官山岩体分布在老庙基山、宝山、罗村以及笔山之间，为一石英闪长岩体，呈岩株状产出，地表形态为不规则的椭圆形，出露面积约1.5 km²。金口岭岩体分布在金口岭及其以北地区，为一石英闪长岩体，呈岩株状产出，地表形态呈北东向带状展布，在其南部出现膨大，出露面积约    5 km²。

区内分布有天马山、黄狮涝硫金矿床和铜官山、金口岭铜矿床等。其中天马山矿床受层位控制明显，矿带长1 400 m，主矿体赋存于天鹅抱蛋山岩体外侧中石炭统黄龙组下段白云岩中及其与上、下岩层之间的层间构造带中，其次赋存于黄龙组上段至栖霞组底部的大理岩中，矿体的形态主要呈似层状，其次为透镜状、筒状、脉状等。黄狮涝金矿位于天马山矿床西南部，处于铜官山背斜的南东翼，矿床主要赋存于五通组砂页岩与黄龙组上段大理岩(局部为下段白云岩)之间的纵断裂带中，其次为黄龙组上下岩性段之间的界面上，矿体主要呈似层状和透镜状产出。铜官山铜矿床处于铜官山背斜的北东倾没端，并环绕铜官山岩体分布，包括松树山、老庙基山、小铜官山、老山、宝山、罗村和笔山等矿段，矿区发育矽卡岩化蚀变，主要铜储量分布在铜官山岩体的东南侧接触带。

2  信息统计单元的划分

目前，在国内外的成矿预测中应用较广的是规则网格单元划分法，它能在统一观察和定量的前提下，把众多的地质变量所包含的矿产资源信息量最大限度地反映出来，这有利于矿与非矿地质特征的判断，并且给矿产预测的计算机网格化带来了方便，尤其是在GIS支持下，网格单元的划分及单元中信息的提取非常便利^[8-9]。

本次研究统计单元的划分主要考虑对矿化的显示，同时又考虑了统计计算、地质信息变量的选取和空间分析等要素。根据矿区的实际情况和统计计算的处理能力，采用规则网格法在1:10 000的铜陵天马山矿区地质图上按200 m×200 m的网格将研究区划分为2 580个信息统计单元(见图1)。

3  预测区地质信息变量的确定及  编码

3.1  地质信息变量的确定

天马山及其外围地区矿床的形成受该地区燕山期岩浆岩、石炭系等地层、褶皱断裂构造及多期成矿事件等多种地质因素复合的结果，这些因素与成矿的关联性及其表现形式，是厘定成矿地质信息变量的基础。

3.1.1  构造信息变量

天马山及其外围地区矿床空间分布明显受褶皱、断裂构造及岩体接触带构造的制约。

天鹅抱蛋山岩体、铜官山岩体及已探明的主要硫金矿床和铜矿床均集中分布在铜官山背斜的北东倾伏端附近，显示了该背斜的空间扩容及其控岩控矿作用。该背斜两翼和西南倾伏端附近有多处矿化点、铁帽和蚀变分布，而次级褶曲或隆起叠加部位是构造应力释放和空间扩容的有利部位，也是找矿的有利线索。

控矿断裂构造以层间断裂为代表，这种断裂主要出现在不同层位界面或不同岩性界面，如C₂h¹/D₃w、C₂h²/C₂h¹、C₃c/C₂h²、P₁q/C₃c层位界面，由于存在着显著的岩性和物理化学性质差异，在构造应力作用下往往容易形成层间滑动和层间虚脱，成为矿区主要的容矿构造并制约了矿体的形态产状。加之石炭系黄龙组白云岩中沉积有胶状黄铁矿层，岩浆热液沿此层间滑脱构造迁移时，与黄铁矿层发生叠加改造，往往形成复合成因的层状矿体。

岩体接触带构造的控矿主要表现为岩体与碳酸盐岩接触带部位发育矽卡岩化并产生空间扩容，形成典型接触带矽卡岩型铜矿体，如铜官山矿床。在天马山矿区，在天鹅抱蛋山岩体的超覆部位及接触带的缓倾部位往往见有硫金矿体。

因此，地表构造行迹可作为判断深部隐伏岩体形态变化及成矿有利程度的重要依据，尤其断裂密集区、断裂交汇部、构造转折部及褶皱、断裂构造的复合部是找矿的有利部位。

3.1.2  地层信息变量

石炭系黄龙组不仅是天马山硫金矿、黄狮涝金矿和铜官山铜矿层状矿体的主要赋矿层位，该层位在整个铜陵地区都是一个非常重要的赋矿层位，在区域上完全可以与冬瓜山铜矿、新桥硫铁矿等对比，表现出明显的层控特点，这在很大程度上与该层位中普遍存在的“同生沉积含矿层”有关，其典型标志之一是早世代生成的胶黄铁矿。该胶黄铁矿层不仅扮演了矿源层的角色，对铜、金等金属元素而言，还起到了沉淀剂的作用。

3.1.3  燕山期侵入岩信息变量

燕山期侵入岩是主要控矿因素和找矿标志之一。岩体与矿体的空间关系、岩体及接触带蚀变特征及岩体的含矿性均显示天马山、黄狮涝、铜官山、金口岭等矿床的形成与燕山期岩浆岩有着密切的成因关系。

在天马山矿区，天鹅抱蛋山岩体与硫金矿床空间关系较密切。岩体接触带附近见有受接触带构造控制的金硫矿体；接触带附近发育矽卡岩化蚀变并具有一定的分带性。据测定，矿区石英闪长岩中金含量为   95×10^-9[10]，明显高于铜官山矿区其它岩体。

黄狮涝金矿距天马山和铜官山岩体较远，但矿区内也见有小规模的石英闪长岩和闪长玢岩岩体分布，岩体与矿体空间关系密切，石炭系和泥盆系地层普遍遭受热变质作用。据测定，石英闪长岩和闪长玢岩中金元素含量较高，分别达到78.6×10^-9和46.0×10^-9 (华东地质调查局812队，1990)。

铜官山和金口岭矿床则以矽卡岩型矿化为主。其中铜官山矿床围绕铜官山岩体接触带分布并发育典型的围岩蚀变分带。据测定，铜官山石英闪长岩岩体中铜元素的含量达126.3×10^-6，金口岭岩体中铜含量达550×10^-6[11]。

3.1.4  地表矿化信息变量

在天马山、铜官山等已知矿床外围尚有多处矿点或铁帽分布，这些矿化异常信息是成矿热液活动的重要线索，尤其在有利的岩浆岩、构造和地层复合部位出现的地表矿化信息，可作为深部找矿预测的重要  依据。

3.1.5  遥感混合蚀变信息变量

燕山期侵入岩，尤其是矿床(点)附近的燕山期侵入岩中普遍发育钾化、内矽卡岩化、黄铁绢英岩化等蚀变，而外接触带围岩中则发育矽卡岩化、硅化、黄铁矿化蚀变等，由于断裂系统的贯通，蚀变的范围往往超出岩体或矿体所在位置，在地表有所显示。

基于ETM数据源，本文作者对天马山及其外围地区遥感混合蚀变信息进行了提取。由于该地区岩石出露程度较高，加之蚀变岩石与广泛出露的碳酸盐岩地层及未蚀变花岗岩在光谱特征上反差明显，对遥感蚀变信息的提取十分有利，所提取的混合蚀变分布与现有矿床分布区具有较高的吻合性，可作为隐伏矿体成矿预测的标志之一。

综合以上分析，天马山及其外围地区的成矿系统是燕山期花岗岩、地质构造、有利地层及多成矿事件耦合等多因素复合作用的产物。控制和影响天马山矿区成矿的地质信息变量非常的复杂，但在提取与找矿有关的信息时，有些信息是定性而不是定量的，因而无法进行统计计算。在统计变量选择的过程中，既要考虑选择的变量便于计算机进行空间分析与计算，又要注重选择有利于成矿预测的代表性变量，因此，为了建立空间定位预测模型，使参与叠加分析的所有图层都包含在研究的预测范围内，特选择以下共11个变量进行空间分析：石炭系、遥感蚀变、铁帽、二叠系、三叠系、岩体接触带、泥盆系、北东向断裂、北西向断裂、构造交汇处、褶皱。

3.2  地质信息变量编码设计

为了有效地组织和管理上述地质信息变量，需要依据变量实体之间不同的特征、相似的特征以及不同变量实体的组合特征来对地质变量进行编码。对地质信息变量的编码设计是在分类的基础上进行的。

本文作者选择的11个变量基本上可以分为5类：地层、构造、岩体、蚀变、矿点。在编码过程中要注意对整个系统的数据进行系统设计、统筹安排，使系统数据编码具有较强的系统性。综合考虑以上原则，结合地学空间数据的特点，参考有关国家标准，本文编码体系见表1~6。

4  地质信息变量赋值

地质变量在GIS中是以层的形式存储在数据库中的，因此，不同的地质变量其属性也不同，但基本上可以分为3类：点文件、线文件和区文件。地质变量取值的实质是统计各网格单元内：是否有点分布在网格内；是否有线通过网格；是否有某个层位的区文件覆盖网格。处理这样的变量在以往的研究中通用的取值方法是二态赋值法，即如果该地质变量在某一划分网格中存在，则其值为1，否则为0。在这里，本文作者采用MAPGIS软件中的空间分析模块对各个地质变量进行叠加分析取值。对不同的地质变量其空间分析的方法不同，具体如下：

1) 点变量取值

判断某个网格单元内是否有点分布，如矿点等。在MAPGIS空间分析模块中空间分析菜单下用区空间分析的区对点相交分析就可得到含矿单元的区文件。这个区文件中就包含了所有取值为1的预测矿点。

表1  地层信息编码表

Table 1  Stratigraphic codes

表2  断裂信息编码表

Table 2  Fault codes

表3  岩体信息编码表

Table 3  Granite body codes

表4  蚀变信息编码表

Table 4  Alteration codes

表5  褶皱信息编码表

Table 5  Fold codes

表6  矿点信息编码表

Table 6  Mineral occurrence codes

2) 线变量取值

判断某个网格单元内是否有线通过，如断裂等。在MAPGIS空间分析模块中空间分析菜单下，用区空间分析的区对线相交分析就可得到有断裂通过单元的区文件。

3) 面变量取值

判断某个网格单元内是否有面通过如地层。首先用空间分析模块的条件检索功能，根据地层代号(若无此属性字段，可在编辑模块中根据地层颜色参数统改层号、改当前层、存当前层等功能)将地层分布图分解成几个区文件，每个文件只包含一个地层单位，有几种地层单位(或岩体)就分为几个区文件。判断某个单元内是否有某一地层出露，可用空间分析模块中的检索菜单下的区域内检索功能，在对话框中选择区域条件文件为地层区文件，被检索文件为网格单元区文件，就可生成有某一地层通过单元的区文件。有几个地层区文件就做几次区域内检索并生成相应数量的区   文件。

利用已划分好的网格进行地质变量取值，这样就形成m(2 580)×n(11)数据矩阵，m表示网格数，n表示变量数。作为参考，这里只列出10个已知矿点的变量取值表(表7)。

5  成矿有利度法的数学描述及其  确定

5.1  数学描述

成矿有利度法是希腊和德国地质学家和数学地质学家合作推出的，该方法在1986年意大利国际数学地质讨论会上受到了各国数学地质工作者的好评。其数学表达式为

                             (1)

式中：f为成矿有利度；w_i为第i个找矿标志的权系数；c_i为第i个找矿标志；p(c_i)为第i个找矿标志出现的概率；N为参加估计的找矿标志个数。

从式(1)可以看出，在成矿有利度法的数学表达式中，各找矿标志的权系数的确定是建模的关键。

变量权系数w可根据下列矩阵方程求得

                               (2)

这里λ是(CC^T)的最大特征值，C是m×n矩阵，代表n个地质变量在m个网格单元上的取值，C^T是C的转置矩阵。

地质变量以二态赋值方式赋值，即预测单元内出现为1，否则为0，数值型变量则以实际数值归一化后赋值。地质变量型找矿预测标志出现的概率以统计方法估计，数值型找矿预测标志的概率以归一化数值替代。根据矩阵表8，应用10个已知矿点组成的数据矩阵，采用MATHCAD数学软件就可以计算出权系数w，代入成矿有利度公式，据此就可以确定找矿预测标志的权系数。然后将各找矿预测标志的权系数经正规化变换，使其和为1，由此可建立铜陵天马山找矿预测数学模型。其主要过程包括：

1) 建立地质变量距阵C并求得转置矩阵C^T(见表9)；

2) 根据地质变量距阵C和转置矩阵C^T，加入中间变量B；B=(CC^T)，求得(CC^T)矩阵(见表9)；

3) 调用eigenval( )函数求得特征矩阵，再调用max(eigenval( ), 0)函数求得最大特征值λ；

4) 最后求得对应最大特征值λ的特征向量eigenvec(B, λ)，即权系数向量W。

(5) 根据权系数向量W，结合成矿有利度式(1)求得找矿预测数学模型表达式如下：

f = 0.089 131 w₁+0.125 361 w₂+0.095 864 w₃+

0.127 605 w₄+0.008 015 w₅+0.058 352 w₆+

0.075 665 w₇+0.116 063 w₈+0.102 918 w₉+

0.108 689 w₁₀+0.092 337 w₁₁                             (3)

这样就可以计算出每个网格单元中的成矿有利度。

5.2  信息统计单元成矿有利度的确定

本研究在铜陵天马山研究区内共划分了2 580个网格信息单元，将信息单元的成矿有利度值按0.1的值域划分为9个信息数据组，并分别统计每组数据的频率(见表10)。据此绘制成矿有利度频率分布图(见图2)，从图2所示频率分布点可确定预测单元的找矿信息临界值为0.6。在全区的2 580个单元中，有203个单元的成矿有利度≥0.6，其中有35个单元为已知有矿单元。

表7  地质变量取值表

Table 7  Values of geological variables

表8  地质变量距阵

Table 8  Matrix of geological variables

表9  地质变量(CC^T)距阵

Table 9  Matrix of geological variables (CC^T)

表10  成矿有利度分级表

Table10  Classification of metallogenic favourability

6  矿床空间信息成矿预测模型的  实现

预测成果输出有两种形式：数据输出和图形输出。数据输出是利用预测模型在数学软件中计算后已将预测结果写入到属性数据表中，通过查找数据表即可得到。再利用生成的信息单元数据成图(见图3)。按成矿有利度0.5、0.6和0.7为异常分界点，对预测单元进行了分级，预测单元可分为3级，即Ⅰ级、Ⅱ级和Ⅲ级，其中Ⅰ级预测单元(大于0.7)为成矿条件最有利，找矿标志明显，找矿潜力大；Ⅱ级预测单元(介于0.6和0.7之间)为成矿条件比较有利，找矿标志较明显，找矿潜力较大；Ⅲ级预测单元(介于0.5和0.6之间)为成矿条件较一般，但仍有成矿可能。

图2  成矿有利度频率分布图

Fig. 2  Metallogenic favourability frequency

根据预测单元计算结果及空间关系，结合地质分析，圈定找矿靶区14处，并按找矿潜力的大小分为A、B、C三级，绘制了天马山及其外围地区找矿靶区预测图(见图3)。

7  结论

1) 根据矿区的实际情况和统计计算的处理能力，采用规则网格法在1:10 000的铜陵天马山矿区地质图上按200 m×200 m的网格将研究区划分为2 580个信息统计单元。

2) 为了建立空间定位预测模型，使参与叠加分析的所有图层都包含在研究的预测范围内，特选择以下11个变量进行空间分析：石炭系、遥感蚀变、铁帽、二叠系、三叠系、岩体接触带、泥盆系、北东向断裂、北西向断裂、构造交汇处、褶皱。

3) 根据所绘制成矿有利度频率分布图，由频率 分布点可确定预测单元的找矿信息临界值为0.6。在 全区的2 580个单元中，有203个单元的成矿有利度≥0.6，其中有35个单元为已知有矿单元。

图3  天马山及其外围地区找矿靶区预测图

Fig. 3  Prospecting targets of Tianmashan and its periphery area: 1—Prediction unit of grade Ⅲ; 2—Prediction unit of grade Ⅱ;   3—Prediction unit of grade Ⅰ; 4—Prospecting target and serial numbers

4) 根据预测单元计算结果及空间关系，结合地质分析，圈定找矿靶区14处，并按找矿潜力的大小分为A、B、C三级，绘制了天马山及其外围地区找矿靶区预测图(见图3)。

REFERENCES

[1] 朱裕生, 肖克炎, 丁鹏飞. 成矿预测方法[M]. 北京: 地质出版社, 1997: 1-23.
ZHU Yu-sheng, XIAO Ke-yan, DING Peng-fei. Metallogcnic prognosis methods [M]. Beijing: Geological Publishing House, 1997: 1-23.

[2] 赵鹏大, 陈永清. 基于地质异常单元金矿找矿有利地段圈定与评价[J]. 地球科学—中国地质大学学报, 1999, 24(5): 443-448.
ZHAO Peng-da, CHEN Yong-qing. Geological anomaly unit-based delineation and assessment of preferable gold ore-finding area [J]. Earth Science—Journal of China University of Geosciences, 1999, 24(5): 443-448.

[3] 王於天. 成矿预测单元的基本概念及其划分方法[J]. 地质论评, 1990, 36(6): 24-29.
WANG Yu-tian. Basic concept of geostatistical cells of minerogenic prognosis and its division method [J]. Geological Review, 1990, 36(6): 24-29.

[4] 李新中, 赵鹏大, 肖克炎, 胡光道. 矿床统计预测单元划分的方法与程序[J]. 矿床地质, 1998, 17(4): 369-375.
LI Xin-zhong, ZHAO Peng-da, XIAO Ke-yan, HU Guang-dao. The method and program for the unit partition in the statistical prediction of mineral deposits [J]. Mineral Deposits, 1998, 17(4): 369-375.

[5] 池顺都, 赵鹏大. 应用GIS圈定找矿可行地段和有利地段—以云南元江地区大红山群铜矿床预测为例[J]. 地球科学—中国地质大学学报, 1998, 23(2): 125-128.
CHI Shun-du, ZHAO Peng-da. Delineating permissive ore-finding area and preferable ore-finding area by GIS: An example from the prediction of copper deposits in Yuanjiang area, Yunnan Province [J]. Earth Science—Journal of China University of Geosciences, 1998, 23(2): 125-128.

[6] 陈石羡. 地理信息系统在金属矿产预测中的应用[J]. 地质找矿论丛, 1998, 13(1): 74-83.
CHEN Shi-xian. Application of GIS to prediction of metal mineral resources [J]. Contributions to Geology and Mineral Resources Research, 1998, 13(1): 74-83.

[7] 陈永清, 夏庆霖. 应用地质异常单元圈定矿产资源体潜在地段—以鲁西铜石金矿田为例[J]. 地球科学—中国地质大学学报, 1999, 24(5): 459-463.
CHEN Yong-qing, XIA Qing-lin. Application of geological anomaly unit method to delineation of potential mineral resources areas: An example from Tongshi gold field, Western Shangdong Province [J]. Earth Science—Journal of China University of Geosciences, 1999, 24(5): 459-463.

[8] 谢贵明, 范继璋. 吉林省珲春东部地区金矿综合信息找矿模型及找矿靶区预测[J]. 黄金科学技术, 2000, 8(5): 20-27.
XIE Gui-ming, FAN Ji-zhang. The synthetic information model of prospecting gold deposits and prediction of prospecting gold targets in east Hunchun district of Jilin province [J]. Gold Science and Technology, 2000, 8(5): 20-27.

[9] 曹 瑜, 胡光道, 杨志峰, 沈珍瑶. 基于GIS有利成矿信息的综合[J]. 武汉大学学报: 信息科学版, 2003, 28(2): 167-176.
CAO Yu, HU Guang-dao, YANG Zhi-feng, SHENG Zhen-yao. Synthesis of beneficial ore-forming information based on GIS [J]. Geomatics and Information Science of Wuhan University, 2003, 28(2): 167-176.

[10] 俞沧海, 袁小明. 铜陵天马山硫金矿床地质特征及成因探讨[J]. 矿产与地质, 2002, 16(2): 74-77.
YU Cang-hai, YUAN Xiao-ming. Geological features and study on genesis of Tianmashan sufur-gold deposit [J]. Mineral Resources and Geology, 2002, 16(2): 74-77.

[11] 袁小明. 铜官山矿田铜金成矿模式探讨[J]. 地球学报, 2002, 23(6): 541-546.
YUAN Xiao-ming. A Tentative discussion on the copper and gold metallogenic model of the Tongguanshan ore field [J]. Acth Geoscientia Sinica, 2002, 23(6): 541-546.

(编辑 何学锋)

基金项目：国家科技支撑计划资助项目(2006BAB01B07)；国家重点基础研究发展计划资助项目(2007CB416608)

收稿日期：2011-12-01；修订日期：2012-02-11

通信作者：杨  斌，高级工程师，博士；电话：0731-88836469; E-mail: 903755562@qq.com