多元统计分析在塔里木瓦吉里塔格钒钛磁铁矿床多元素相关性及成矿预测研究中的应用

2016-12-20赵亚莉刘浩滢厉子龙邹思远孙浩伟励音骐杨树锋陈汉林

浙江大学学报（理学版） 2016年6期

关键词：辉石塔格磁铁矿

赵亚莉，刘浩滢,，厉子龙，邹思远，孙浩伟，励音骐，杨树锋，陈汉林，成军

(1. 渤海大学数理学院，辽宁锦州121000； 2. 浙江大学地球科学学院，浙江杭州310027； 3. 河北地质大学资源学院，河北石家庄 050031； 4. 新疆八一钢铁集团有限责任公司，新疆乌鲁木齐830022)

表1 巴楚瓦吉里塔格辉长岩中主量元素Pearson相关系数

表2 巴楚瓦吉里塔格辉石岩中主量元素Pearson相关系数

表3 巴楚瓦吉里塔格橄辉岩中主量元素Pearson相关系数

表4 巴楚瓦吉里塔格橄榄石中主量元素Pearson相关系数

表5 单斜辉石中主量元素Pearson相关系数

表6 磁铁矿中主量元素Pearson相关系数

表7 钛铁矿中主量元素Pearson相关系数

表8 钛磁铁矿中主量元素Pearson相关系数

表9 钻孔中Fe-Ti-V数据相关统计值

表10 Fe-Ti-V氧化物Pearson相关系数

表11 lnFe2O3T、lnTiO2和lnV2O5回归分析结果

多元统计分析在塔里木瓦吉里塔格钒钛磁铁矿床多元素相关性及成矿预测研究中的应用

赵亚莉1，刘浩滢1,2，厉子龙2*，邹思远2,3，孙浩伟2，励音骐2，杨树锋2，陈汉林2，成军4

塔里木巴楚县瓦吉里塔格地区出露有与塔里木早二叠世大火成岩省的形成密切相关的大型钒钛磁铁矿床. 基于多元统计分析对该矿区露头和大量钻孔岩心样品中Fe-Ti-V元素和部分主量元素数据开展了元素地球化学相关性研究. 采用SPSS软件对大量数据进行相关分析和回归分析，计算各元素间的相关系数，Fe、Ti、V 3种元素的相关性分析结果表明: 矿床中Fe、Ti和V之间均呈正相关, 成矿元素与造岩元素之间存在负相关关系；建立伴生钒元素与铁和钛元素之间关系的回归方程V=-4.984+0.360 Fe+0.984 Ti，由Fe和Ti元素质量分数可以估算伴生组分V的质量分数.利用Matlab软件的可视化功能，绘制了Fe-Ti-V元素质量分数的三维立体图，建立空间模型，发现Fe2O3T、TiO2和V2O5品位较高的位置相对一致，均主要分布于研究区北部中间一带，富集在辉石岩中. 该空间模型的建立有助于识别矿体和非矿体的分布范围和富集趋势，并进行成矿预测. 此项研究对钒元素含量的精确估算和矿藏勘探有一定的指导意义.

相关分析；回归分析；成矿预测；钒钛磁铁矿；塔里木大火成岩省；SPSS；Matlab

0 引言

大火成岩省是指在一段相对短暂的地质时期内，由于持续或脉动式的岩浆作用在板内或者与之具有相似地球化学亲和性的构造环境中所形成的大规模岩浆岩建造，岩性主要为铁镁质喷出岩或侵入岩[1].在大火成岩省的形成过程中，往往伴随着地球深部物质和能量的强烈作用和交换，成矿物质的聚集导致元素的富集和矿床的形成，因此大火成岩省的形成往往与一个大的成矿系统相对应.其形成的矿床类型主要有：Ti-Fe氧化物型矿床、Cu-Ni硫化物型矿床和热液型矿床[2-3].大火成岩省的岩浆活动与钒钛磁铁矿床有密切的成因关系，如峨眉山大火成岩省与攀枝花和红格钒钛磁铁矿床的成因有联系[4].

目前，在塔里木大火成岩省中唯一正在开采的大型钒钛磁铁矿床位于新疆巴楚县瓦吉里塔格地区[5]，已有大量学者和地质工作者对该区域进行了地质调查和专题研究.发现瓦吉里塔格是一个集铁、钒、钛、稀土、磷、金刚石为一体的矿床，产出有大型钒钛磁铁矿床、碳酸岩型稀土矿、碳酸岩型磷灰石矿及含金刚石的金伯利岩[6-7].与基性-超基性杂岩体密切相关的航磁异常以及岩体中辉石岩、橄辉岩及铁矿石的比重较大，瓦吉里塔格地区的这些地球物理特征都为该地区的矿藏勘查工作提供了重要的基础资料.此外，通过对瓦吉里塔格钒钛磁铁矿床中基性-超基性岩石地球化学的研究，证明其原始岩浆为演化的富Fe-Ti的玄武质岩浆，来自于富集的地幔柱源区[8].

但是，对瓦吉里塔格钒钛磁铁矿床多种元素之间的相关性和内在联系及铁-钛-钒在空间上的富集规律和分布特征尚不清楚，也未见针对该矿区的多元统计分析的相关工作，阻碍了深入探讨成矿机理和矿藏勘探研究.

多元统计分析是研究变量之间相互依赖关系的统计规律性的方法，主要包括相关分析、回归分析、因子分析、聚类分析和趋势面分析等.其中，相关分析是通过建立合理的指标来衡量变量之间关系密切程度的统计方法[9].而回归分析是处理变量间关系并且进行系统推断的一种统计方法和技术，通过分析大量样本数据，确定变量间的数学关系、建立回归方程，对所确定的数学关系式的可信程度进行检验，估计回归效果，讨论每个自变量对回归的贡献大小，选择最优回归方程，并对因变量进行预测[10].

相关分析和回归分析可以很大程度地从数据中还原地质信息，在开展元素的地球化学研究时，通过判断金属元素复杂数据的相关关系来研究矿区的地质特征[11].前人通过相关分析挖掘变量中的有用信息，利用相关系数实现矿产资源的预测[12]，应用多元回归分析和逐步回归分析，通过建立回归方程，验证并解释矿床成因、岩浆来源、矿产分布规律和金属质量分数预测等[13-15]，同时在生物成矿研究上也有所应用[16].

因此，本文通过对塔里木巴楚县瓦吉里塔格钒钛磁铁矿床含矿岩体的主量元素以及铁、钛、钒3种元素氧化物之间进行相关分析和回归分析，建立回归方程，从而精细刻画元素之间的相关性与岩浆作用过程的联系，预测伴生组分钛的质量分数，并采用Matlab建立三维立体模型以预测矿体分布规律.

1 区域地质背景与矿区地质概况

塔里木早二叠世大火成岩省是塔里木地质演化历史中岩浆活动最为强烈、影响范围最广的一次地质热事件，大规模溢流玄武岩的覆盖面积超过25×104km2[17]，主要分布于塔里木盆地的中西部地区(见图1)，其周缘被天山造山带、昆仑造山带和阿尔金造山带围绕[18-19].大量的年代学和地球化学研究表明，塔里木大火成岩省的岩浆演化主要划分为早期喷发的大陆溢流玄武岩(290～285 Ma)和晚期侵位的基性-超基性岩及碱性正长岩类侵入岩体和岩脉(284～274 Ma)2个阶段.且各种岩石类型发育的时间序列如下：库普库兹曼组玄武岩→开派兹雷克玄武岩→隐爆角砾岩岩筒→瓦吉里塔格地区基性-超基性层状杂岩体→辉绿岩脉、超基性岩脉→正长岩[20].

瓦吉里塔格地区位于新疆喀什地区巴楚县境内，处于塔里木盆地中央隆起带的西北边缘.区内出露一套由基性-超基性层状侵入岩体(其中主要包括橄辉岩、辉石岩和辉长岩)、基性岩墙群、金伯利质隐爆角砾岩和碱性岩体组成的火成杂岩体(即瓦吉里塔格杂岩体)，并产出大型钒钛磁铁矿床[18,21-24].该杂岩体长约5 km，宽1.5～3 km，地表出露面积约12 km2(见图2)，在地形上呈南北轴稍长、东西端收敛的鸭梨形状[25].其中，基性-超基性层状岩体内部单元主要岩石类型自下而上为橄辉岩、辉石岩和辉长岩，各单元岩石呈层状或似层状产出[26].

本文研究的钒钛磁铁矿床位于瓦吉里塔格地区南部(见图2)，据航磁资料，瓦吉里塔格杂岩体是引起该地区强负磁异常的主要磁性体，也是前人找矿的主要目标物[25].矿区内主要矿化类型为钒钛磁铁矿化，多分布于基性-超基性岩体中，呈中等至弱浸染状，当磁铁矿富集到一定程度时，即形成钒钛磁铁矿体.根据相关地质资料的介绍，截至目前已控制全铁品位大于15%的铁矿石量超过4亿t，同时伴生TiO2金属量2 000万t，平均品位6.99%，伴生V2O5金属量15 000 t，平均品位0.15%.其成因类型属于岩浆晚期分异型，矿石类型为辉长岩-辉石岩-橄辉岩型，矿体主要赋存于辉石岩相带，与我国四川著名的红格钒钛磁铁矿床类似.

图1 塔里木盆地大地构造位置图(a)及塔里木区域地质简图(b)(据文献[22]修改)Fig.1 Geotectonic distribution of Tarim Basin(a) and regional geological sketched map(modified after Ref.[22])

图2 瓦吉里塔格矿区地质简图(a)及钻孔所在位置示意图(b)(据文献[19]修改)Fig.2 Geological sketched map of the Wajilitag area(a) and location diagram of drills (b) (modified after Ref.[19])

2 多元统计分析方法

相关分析、回归分析均通过SPSS软件(SPSS 19.0)实现.SPSS软件(Statistical Product and Service Solutions)是集数据录入、资料编辑、数据管理、统计分析和图形绘制等为一体的统计功能完善的软件,可以进行相关分析、因子分析、聚类分析、回归分析等复杂的多因素统计分析[27].

Matlab除了具有符号计算以及数值计算的功能外，还具有强大的数据可视化能力[28].利用Matlab(2014a)三维图形显示矿体特征，能更加直观地反映矿体的变化趋势、矿化深度、品位及形态，对找矿具有指导意义[29].

将采样点的数据以矩阵的形式输入，首先将零散数据规则化，即通过插值的方法，将数据转换成规则的网格数据；然后调用surf绘图函数绘图[28].

3 样品主量元素和钻孔中铁、钛和钒元素的多元统计分析

先对塔里木巴楚县瓦吉里塔格钒钛磁铁矿基性-超基性层状侵入岩体中主要岩类(包括橄辉岩、辉石岩、辉长岩)中的9种氧化物[30]进行相关分析和回归分析. 以SiO2、Al2O3、MnO、MgO、CaO、Na2O和K2O为变量，分别对Fe2O3T和TiO2进行逐步回归分析；随后对橄榄石、单斜辉石、磁铁矿、钛铁矿以及伴生的钛磁铁矿等矿物成分[30]之间进行逐步回归分析；同时，选取瓦吉里塔格南部钒钛磁铁矿区4条勘探线剖面16个钻孔共计2 603个数据(主要据新疆地矿局第十一地质大队瓦吉里塔格地质详查报告，2009)，对其中的Fe2O3T、TiO2、V2O53种氧化物质量分数进行多元线性回归分析，得到一系列经验回归方程，通过方程对伴生微量元素钒的质量分数进行估测.并根据钻孔坐标和Fe2O3T、TiO2、V2O5的含量，利用Matlab绘制矿体的空间展布特征并进行预测.

3.1 瓦吉里塔格层状杂岩体不同岩类的主量元素相关分析

分别对瓦吉里塔格含钒钛磁铁矿层状侵入岩体中辉长岩、辉石岩和橄辉岩的主量元素数据进行相关性分析(见表1～3)知，辉长岩中Fe2O3T与SiO2、Al2O3、MgO、Na2O相关性显著，且与MgO呈正相关，而与SiO2、Al2O3、Na2O呈负相关(见表1)，表明随着岩浆结晶分异的作用，不断有磁铁矿析出.辉石岩的相关性分析表明，Fe2O3T与SiO2、TiO2与SiO2和MgO均呈显著负相关(见表2)，说明随着SiO2含量的增加，不断有磁铁矿和钛铁矿形成.橄辉岩中，TiO2与SiO2呈显著负相关(见表3)，表示结晶过程中有钛铁矿不断析出.

表1 巴楚瓦吉里塔格辉长岩中主量元素Pearson相关系数

Table 1 Pearson correlation index for major elements of gabbros in the Wajilitag ore deposit of Tarim Basin

注 **指在0.01 水平上显著相关； *指在0.05 水平上显著相关.

表2 巴楚瓦吉里塔格辉石岩中主量元素Pearson相关系数

Table 2 Pearson correlation index for major elements of pyroxenites in the Wajilitag ore deposit of Tarim Basin

注 **指在0.01 水平上显著相关；*指在0.05 水平上显著相关.

表3 巴楚瓦吉里塔格橄辉岩中主量元素Pearson相关系数

Table 3 Pearson correlation index for major elements of olivine pyroxenites in the Wajilitag ore deposit of Tarim Basin

注 **指在0.01 水平上显著相关；*指在0.05 水平上显著相关.

此外，对橄榄石、单斜辉石、钛铁矿、磁铁矿、钛磁铁矿等单矿物中氧化物数据进行相关性分析(见表4～8).橄榄石中SiO2、Fe2O3T、MnO、MgO 4种氧化物之间均呈显著相关性(见表4)，其中SiO2与MgO相关性较好(0.641)，表明在橄榄石形成过程中结晶分异占主导；单斜辉石中Fe2O3T与除MnO以外的其他氧化物均有较好的相关性，TiO2与SiO2、Al2O3、Fe2O3T、MgO和Na2O相关性显著(见表5)；在磁铁矿、钛铁矿和钛磁铁矿中(见表6～8)，Fe2O3T与TiO2均存在明显的负相关(相关系数依次为-0.957，-0.818和-0.996)，造成这一现象的原因可能与Fe2O3T与TiO2在矿物形成过程中呈现类质同象有关.

表4 巴楚瓦吉里塔格橄榄石中主量元素Pearson相关系数

Table 4 Pearson correlation index for major elements of olivine in the Wajilitag ore deposit of Tarim Basin

注 **指在0.01水平上显著相关；*指在0.05水平上显著相关.

表5 单斜辉石中主量元素Pearson相关系数

Table 5 Pearson correlation index for major elements of clinopyroxene

注 **指在0.01 水平上显著相关；*指在0.05水平上显著相关.

表6 磁铁矿中主量元素Pearson相关系数

Table 6 Pearson correlation index for major elements of magnetite

注 **指在0.01水平上显著相关；*指在0.05水平上显著相关.

表7 钛铁矿中主量元素Pearson相关系数

Table 7 Pearson correlation index for major elements of ilmenite

注 **指在0.01水平上显著相关.

表8 钛磁铁矿中主量元素Pearson相关系数

Table 8 Pearson correlation index for major elements of coexisting magnetite-ilmenite

注 **指在0.01水平上显著相关；*指在0.05水平上显著相关.

3.2 瓦吉里塔格矿区钻孔中Fe2O3T、TiO2和V2O5质量分数的相关性分析和回归分析

本文详细统计了瓦吉里塔格钒钛磁铁矿床4个勘探线剖面16个钻孔2 603组(N)样品数据每个变量的平均值、最大值、最小值、方差等基本信息(见表9).对其中的Fe2O3T、TiO2、V2O53个变量进行相关分析(见表10)以及多元线性回归分析(见表11)，且对钻孔中辉石岩、辉长岩、橄辉岩等岩石中3种元素进行逐步回归分析(见表11).

表9 钻孔中Fe-Ti-V数据相关统计值

Table 9 Statistical values of borehole data of Fe-Ti-V

在所有钻孔样品数据中Fe2O3T的质量分数为4.08%～48.72%，平均值是15.94%，方差8.11，数据变化比较稳定；对TiO2质量分数的统计得出平均质量分数是6.40%，最小值为1.39%，最大值为19.40%，方差2.48，变化很小，数据比较集中；微量元素V2O5的质量分数多在0.12%附近，最小不低于0.03%，最大值小于0.94%(见表9).

表10 Fe-Ti-V氧化物Pearson相关系数

Table 10 Pearson correlation index of Fe-Ti-V oxide

注 **指在0.01水平上显著相关.

从可决系数AdjustedR2上对比2个回归方程，不难发现，矿体所在位置的Fe2O3T、TiO2和V2O5数据所建立的回归方程拟合效果更好，这是由于矿体所在位置3种元素的质量分数相对稳定，且正相关性显著.由以上多元线性回归方程，可以由Fe2O3T以及TiO2的质量分数计算得V2O5的估计值，对比实验分析测试值与经验方程计算的估计值，发现绝大多数数据的误差在0.03%以内，但是仍有部分数据误差较大(见图3)，为得到更准确的经验方程还需要进一步的研究.总体来看，所得回归方程有一定的应用价值，将所得经验方程应用于野外现场X射线荧光光谱测量中，通过在野外露头直接测得的铁和钛质量分数由回归方程估算钒的质量分数，可以大大提高勘查工作的效率.

表11 lnFe2O3T、lnTiO2和lnV2O5回归分析结果

Table 11 Regression analytic results of Fe-Ti-V oxides

注表中Fe、Ti、V分别为Fe2O3T、TiO2、V2O5的对数值.

图3 V2O5观测值与预测值对比Fig.3 Comparison of the observed values and predicted values of V2O5 contents

4 应用Matlab数据可视化预测矿体潜在富集区域

在Matlab软件中，利用surf等绘图函数绘制三维曲面图.首先对矿区元素质量分数数据和钻孔坐标数据进行整理并进行网格化；然后输入函数H=surf(X，Y，Z，C), 此处X，Y为数据点的空间坐标，Z为高程值，C为铁(或钛或钒)质量分数；最后得到三维曲面图，且通过色标的变化体现不同位置金属元素的富集程度.

本文对研究区内体积约为3.5×107m3的含矿岩体分层，以约30 m的深度为一层，针对每一层数据绘制三维曲面，将11张曲面图整合到一起，体现出整个岩体的Fe2O3T、TiO2和V2O5质量分数展布特征(见图4～6)，从展布图上可以直观地看到Fe2O3T、TiO2、V2O5质量分数富集区域以及潜在富集区的空间分布规律，结合矿区地质情况，进一步预测矿体走势.

图4 瓦吉里塔格矿区Fe2O3T质量分数空间分布Fig.4 Spatial distribution of Fe2O3T values in the Wajilitag ore deposit

图5 瓦吉里塔格矿区TiO2质量分数空间分布Fig.5 Spatial distribution of TiO2 values in the Wajilitag ore deposit

图6 瓦吉里塔格矿区V2O5质量分数空间分布Fig.6 Spatial distribution of V2O5 values in the Wajilitag ore deposit

在Fe2O3T质量分数空间分布立体图中(见图4)，图标尺颜色的变化表示Fe2O3T质量分数的变化趋势，颜色越红表明Fe2O3T的品位越高.不难发现，研究区的北部中间一带Fe2O3T品位较高，钻孔zk24、zk21、zk32和zk31附近，海拔1 000 m到接近地表位置，Fe2O3T品位平均在20%左右，主要分布于辉石岩岩相中.相比之下，研究区的东北角和西南角，尤其是接近地表部分Fe2O3T品位较低.据此，下一步可以考虑向研究区的西北方向继续布钻(见图4)，证实是否为潜在的铁富集区.在TiO2质量分数空间分布立体图中，图标尺颜色的变化表示TiO2品位的变化趋势，在研究区北部中间一带(zk24和zk32附近)海拔1 050～1 150 m处及海拔1 000 m以下的位置钛质量分数几乎在8%以上且比较富集，此外研究区东北部zk24钻孔中930～1 030 m层位上钛元素质量分数可达到9%～10%，已达工业品位，具有开采价值，可以考虑在此开采甚至向更深的位置施工(见图5).在V2O5质量分数空间分布立体图中(见图6)，V2O5质量分数整体偏低，从颜色变化来看，在研究区中间区域，钻孔zk12、zk21、zk22、zk32和zk31附近钒相对富集，且北部边缘1 050 m以上位置元素质量分数在0.16%以上，可以推测这是伴生的钒的富集区(见图6).

对比图4～6，铁、钛和钒3种元素在北部中间一带质量分数普遍高，富集区域相同，与铁、钛、钒3种元素呈显著正相关的计算结果一致，导致这种现象的原因可能是铁、钛和钒3种元素同为过渡元素可呈类质同象，具有相似的迁移、沉淀的物理化学特性及富集机理.同时，在已知的勘探区内，钒钛磁铁矿矿化区主要集中在北部区域，赋存于辉石岩中.综合研究区的地质特征，该含矿的基性-超基性层状岩体产状近水平，向北展布达2～3 km，结合航磁异常的资料，该研究区向北的层状岩体出露区域仍处于强负磁异常的区域.因此，推测钒钛磁铁矿体在深部具有明显的向西北富集延伸的趋势，下一步值得往北(即瓦吉里塔格中区)布钻和勘查，以验证该区是否为钒钛磁铁矿富集的远景区之一.

5 结论

通过数理统计分析(包括相关分析和回归分析)方法对塔里木大火成岩省中巴楚瓦吉里塔格钒钛磁铁矿床中全岩主量成分、单矿物成分和V-Ti-Fe元素间进行较深入的分析，并开展V-Ti-Fe元素富集特征的三维立体建模，得到以下3点结论：

(1)辉长岩中Fe2O3T与MgO显著正相关，主要由橄榄石结晶作用引起.辉石岩中TiO2和MgO呈负相关，由于它们之间具有不同的富集机制，此成矿元素与造岩元素之间存在负相关关系.通过对钛铁矿、磁铁矿、钛磁铁矿中Fe2O3T和TiO2的分析发现，Fe与Ti成线性相关，且相关系数均为负值，造成这种关系的原因可能是Fe与Ti的类质同象置换.瓦吉里塔格钒钛磁铁矿床中Fe、Ti、V均成正相关，表明它们主要受控于钛磁铁矿和钛铁矿，说明钒、钛和铁在矿石中的富集机理相似.

(2)对综合矿区大量钻孔数据中Fe2O3T、TiO2和V2O5质量分数进行多元线性回归分析，通过对回归方程的复相关系数检验、F统计量检验、T统计量检验，证明所得经验方程的拟合效果良好；应用逐步回归分析方法把Fe、Ti元素代入回归方程得到伴生微量元素V的回归估算值，用回归估算值代替实测值以提高工作效率，同时对野外现场矿藏勘探有一定的参考意义.

(3)通过Matlab三维立体建模，可以快速直观地分析矿石的品位和分布.在瓦吉里塔格钒钛磁铁矿中，Fe2O3T、TiO2和V2O5质量分数相对较高的区域主要在辉石岩，在元素质量分数空间分布立体图中通过颜色的变化体现研究区岩体中金属元素的富集程度、矿体走向和向北延伸趋势，对矿体勘查具有重要作用.

长安大学汤中立院士、新疆地质矿产勘查开发局第二地质大队冯晶荣总工、新疆地矿局第十一地质大队和新疆八一钢铁股份有限公司瓦吉里塔格矿山项目管理部提供了宝贵的矿区地质图、勘探线剖面图和相关矿石样品及钻孔中金属元素分析数据等，在此一并致谢！

[1] BRYAN S E, ERNST R E. Revised definition of Large Igneous Provinces (LIPs) [J]. Earth-Science Reviews,2008,86(1/2/3/4):175-202.

[2] 肖龙，PIRAJNO F，何琦.试论大火成岩省与成矿作用[J].高校地质学报，2007，13(2):148-160. XIAO Long, PIRAJNO F, HE Qi. A preliminary discussion on large igneous provinces and associated ore deposits [J]. Geological Journal of China Universities,2007,13(2):148-160.

[3] PIRAJNO F, ERNST R E, BORISENKO A S, et al. Intraplate magmatism in Central Asia and China and associated metallogeny [J]. Ore Geology Reviews,2009,35:114-136.

[4] ZHOU M F, ROBINSON P T, LESHER C M, et al. Geochemistry, petrogenesis, and metallogenesis of the Panzhihua gabbroic layered intrusion and associated Fe-Ti-V oxide deposits, Sichuan Province, SW China [J]. Journal of Petrology,2005,46(11):2253-2280.

[5] 高玉山.瓦吉里塔格钒钛磁铁矿地质特征及找矿建议[J].新疆钢铁，2007(2):8-9. GAO Yushan. Geological characteristics of Wajilitag vanadic titanomagnetite deposit and its prospecting recommendations[J]. Xinjiang Iron and Steel,2007(2):8-9.

[6] 苏黎.新疆巴楚瓦吉里塔格金伯利岩矿物中岩浆包裹体岩浆研究[J].西北地质科学，1991(2):33-46. SU Li. A research on magmatic inclusions in minerals from kimberlite in Bachu county, Xinjiang, China[J]. Northwest Geoscience,1991(2):33-46.

[7] 宋海涛.新疆巴楚县瓦吉尔塔格基性-超基性岩成矿背景探讨[J].新疆钢铁，2009(1)：54-55. SONG Haitao. Discussion the metallogenic geological background of Wgrtag basic-ultrabasic rocks at Bachu area in Xinjiang [J]. Xinjiang Iron and Steel,2009(1):54-55.

[8] CAO J, WANG C Y, XING C M, et al. Origin of the early Permian Wajilitag igneous complex and associated Fe-Ti oxide mineralization in the Tarim Large Igneous Province, NW China[J]. Journal of Asian Earth Sciences,2014,84:51-68.

[9] 魏立力,马江洪，颜荣芳.概率统计引论[M].北京：科学出版社，2012. WEI Lili, MA Jianghong, YAN Rongfang. An Introduction to Probability and Statistics[M]. Beijing: Science Press,2012.

[10] 何晓群，刘文卿.应用回归分析[M].北京：中国人民大学出版社,2002. HE Xiaoqun, LIU Wenqing. Applied Regression Analysis[M]. Beijing: China Renmin University Press,2002.

[11] CLOSS L G, NICHOL I. The role of factor and regression analysis in the interpretation of geochemical reconnaissance data [J]. Canadian Journal or Earth Sciences,1975(12):1316-1330.

[12] 谈小生.相关系数及其在矿产资源预测中的应用[J].长春地质学院学报，1991，21(4):469-474. TAN Xiaosheng. The relative coefficients and its application to mineral resource prediction [J]. Journal of Changchun University of Earth Science,1991,21(4):496-474.

[13] 付勇,汪立今，柴凤梅，等.多元线性回归和逐步回归分析在白石泉Cu-Ni硫化物矿床研究中的应用[J].地学前缘，2009，16(1):373-380. FU Yong, WANG Lijin, CHAI Fengmei, et al. The application of multivariate linear regression combined with stepwise regression analysis for Cu-Ni sulfide deposit of Baishiquan [J]. Earth Science Frontiers,2009,16(1):373-380.

[14] 刘永刚，张渊逊，吴志之，等.一元线性回归分析方法在玉龙铜矿床的应用及验证[J].中国矿山工程，2012，41(6):4-7. LIU Yonggang, ZHANG Yuanxun, WU Zhizhi, et al. Application and verification of single linear regression analysis method in Yulong copper deposit [J]. China Mine Engineering,2012,41(6):4-7.

[15] 江林.多元统计分析在核勘查数据处理中的应用研究[D].成都：成都理工大学，2011. JIANG Lin. Multivariate Statistical Analysis Applied and Studied in Nuclear Survey Data Processing [D]. Chengdu: Chengdu University of Technology，2011.

[16] 朱弟成，林丽，朱利东.多元回归分析在生物成矿研究中的应用[J].成都理工学院学报，1999，26(3):308-312. ZHU Dicheng, LIN Li, ZHU Lidong. The application of multiple regression analysis to the study of biomineralization [J]. Journal of Chengdu University of Technology,1999,26(3):308-312.

[17] 杨树锋，陈汉林，冀登武，等.塔里木盆地早-中二叠世岩浆作用过程及地球动力学意义[J].高校地质学报，2005，11(4)：504-511. YANG Shufeng, CHEN Hanlin, JI Dengwu, et al. Geological process of early to middle Permian magmatism in Tarim Basin and its geodynamic significance[J]. Geological Journal of China Universities,2005,11(4):504-511.

[18] 余星.塔里木早二叠世大火成岩省的岩浆演化与深部地质作用[D].杭州：浙江大学，2009. YU Xing. Magma Evolution and Deep Geological Process of Early Permian Tarim Large Igneous Province [D]. Hanghzou: Zhejiang University,2009.

[19] 励音骐.塔里木早二叠世大火成岩省岩浆动力学及含矿性研究[D].杭州：浙江大学，2013. LI Yinqi. Study on Magma Dynamics and Ore Potential of the Early Permian Tarim Large Igneous Province [D]. Hangzhou: Zhejiang University,2013.

[20] LI Z L, CHEN H L, SONG B, et al. Temporal evolution of the Permian Large Igneous Province in Tarim Basin in Northwestern China [J]. Journal of Asian Earth Sciences,2011,42:917-927.

[21] ZHANG C L, LI X H, LI Z X, et al. A Permian layered intrusive complex in the Western Tarim Block, Northwestern China: Product of a ca. 275 Ma mantle plume?[J]. Journal of Geology,2008,116(3):269-287.

[22] 邹思远，厉子龙，任钟元，等.塔里木柯坪地区二叠系沉积岩碎屑锆石U-Pb定年和Hf同位素特征及其对塔里木块体地质演化的限定[J].岩石学报，2013，29(10)：3369-3388. ZOU Siyuan, LI Zilong, REN Zhongyuan, et al. U-Pb dating and Hf isotopic compositions of detrital zircons from Permian sedimentary rocks in Keping area of Tarim Basin, Xinjiang, China: Constraints on geological evolution of Tarim Block [J]. Acta Petrologica Sinica,2013,29(10):3369-3388.

[23] ZHANG C L, XU Y G, LI Z X, et al. Diverse Permian magmatism in the Tarim Block, NW China: Genetically linked to the Permian Tarim mantle plume?[J]. Lithos,2010,119(3/4):537-552.

[24] ZHOU M F, ZHAO J H, JIANG C Y, et al. OIB-like, heterogeneous mantle sources of Permian basaltic magmatism in the western Tarim Basin, NW China: Implications for a possible Permian Large Igneous Province [J]. Lithos,2009,113(3/4):583-594.

[25] 李昌年，路凤香，陈美华.巴楚瓦吉里塔格火成杂岩体岩石学研究[J].新疆地质，2001，19(1)：38-42. LI Changnian, LU Fengxiang, CHEN Meihua. Research on petrology of the Wajilitag complex body in North edge in the Talimu Basin [J]. Xinjiang Geology,2001,19(1):38-42.

[26] 徐义刚，王焰，卫荀，等.与地幔柱有关的成矿作用及其主控因素[J].岩石学报，2013，29(10):3307-3322. XU Yigang, WANG Yan, WEI Xun, et al. Mantle plume-related mineralization and their principal controlling factors [J]. Acta Petrologica Sinica,2013,29(10):3307-3322.

[27] 卢纹岱.SPSS统计分析[M].北京：电子工业出版社，2011. LU Wendai. SPSS Statistic Analysis[M]. Beijing: Publishing House of Electronics Industry,2011.

[28] 苏金明，傅荣华，胡卸文.基于Matlab的地学可视化实现[J].地质灾害与环境保护，2009，20(1):83-86. SU Jinming, FU Ronghua, HU Xiewen. The implementation of geo-visualization based on Matlab [J]. Journal of Geological Hazards and Environment Preservation,2009,20(1):83-86.

[29] 齐建全，荣辉.Matlab三维可视化编程在砂岩型铀矿勘查中的应用[J].铀矿地质，2013，29(6):352-356. QI Jianquan, RONG Hui. The application of Matlab based 3D visualization in the exploration of sandstone-type Uranium deposit [J]. Uranium Geology,2013,29(6):352-356.

[30] LI Y Q, LI Z L, CHEN H L, et al. Mineral characteristics and metallogenesis of the Wajilitag layered mafic-ultramafic intrusion and associated Fe-Ti-V oxide deposit [J]. Journal of Asian Earth Sciences,2012,49:161-174.

ZHAO Yali1, LIU Haoying1,2, LI Zilong2, ZOU Siyuan2,3, SUN Haowei2, LI Yinqi2, YANG Shufeng2, CHEN Hanlin2, CHENG Jun4

(1.DepartmentofMathematicsandPhysics,BohaiUniversity,Jinzhou121000,LiaoningProvince,China; 2.SchoolofEarthSciences,ZhejiangUniversity,Hangzhou310027,China; 3.TheCollegeofResources,HebeiGEOUniversity,Shijiazhuang050031,China; 4.XinjiangBayiSteelCompany,Urumqi830022,China)

Multivariate statistical analyses for Fe-Ti-V oxide ore deposit in Wajilitag area of Tarim Basin and its application for metallogenic prediction. Journal of Zhejiang University(Science Edition), 2016,43(6):716-725

The giant Wajilitag Fe-Ti-V oxide ore deposit outcropping in the Bachu area of Tarim Basin, is closely linked with the formations of the Early Permian Tarim Large Igneous Province. The correlations of major elements and metal elements (Fe, Ti and V) of the rock and ore samples from both outcrops and drill holes were studied by multivariate statistical analyses. The relationship between Fe-Ti and major elements was studied by the method of correlation analysis and regression analysis using the SPSS software. Correlation analysis of the drill hole data shows that Fe, Ti and V contents have positive correlations with each other. The equation of V=-4.984+0.360 Fe+0.984 Ti among Fe, Ti and V elements is also established to predict vanadium content when Fe and Ti contents are known. The Fe-Ti-V three-dimensional cubic charts are drawn to establish the spatial distribution mode and predict the potential ore body using the visualization of Matlab. Our study shows that the contents of Fe2O3T、TiO2and V2O5are mostly enriched in pyroxenites in the middle area of the northern part of the Wajilitag ore deposit. This study addresses an effective role in delineating and tracing the distribution trend and the enriched potential of Fe, Ti and V ore bodies.

correlation analysis; regression analysis; metallogenic prediction; Fe-Ti-V oxide ore deposit; Tarim Large Igneous Province; SPSS; Matlab