对应分析发在地浸铀矿山水质评价中地应用.docVIP

对应分析法在地浸铀矿山水质评价中的应用 利广杰核工业北京化工冶金研究院，北京101149） 究重点[1]Method)又称关联分析法，是由法国人Benzencri于1969年提出的。对应分析的关键是使型和型分析中变量点和样品点的载荷能反映在相同的公因子轴上，因此可以用相同的因子轴去同时表示变量和样品。由于型因子分析和型因子分析都能反映一个整体的不同侧面，它们之间一定存在内在联系。对应分析就是通过一个过渡矩阵将两者有机的结合起来。对应分析法最大特点是能把众多的样品和众多的变量同时作到同一张图解上，将样品的大类及其属性在图上直观而又明了地表示出来，具有直观性。另外，它还可以从因子载荷图上对样品进行直观的分类，而且能够指示分类的主要参数(主因子)以及分类的依据，是一种直观、简单、方便的多元统计方法[2]。其方法的具体步骤为： 1.1 由有个样品个变量的原始资料阵出发，计算标度化的概率矩阵： 其中： 1.2 计算过渡矩阵； 1.3 进行因子分析 1.3.1 进行型因子分析 计算协方差距阵的特征值，在按照其累积百分比，取其前个特征值，。型因子载荷矩阵。 变量 因子 在两因子轴平面上作变量点图。 1.3.2 型因子分析 由于型因子分析和因子分析型的非零特征值相同，型与型的具有相同的特征值，故用型分析中所得的特征值计算其对应与矩阵相应的单位特征向量，，从而得到型因子载荷矩阵。 样品 因子 1 2 在与型相应的因子平面上作样品点图[3]。 2 应用实例 以新疆某地浸铀矿山地下水质为例，自2009年7月至2010年6月共有12个样品，对每个样品分析8个变量（即8种化学成分）。8个变量依次为U、Ca2+ 、Mg2+、Fe3+、∑Fe、SO42-、pH、Eh。这些监测数据构成了一个具有12行8列的原始数据矩阵，见表1。 依据进行型因子分析，计算矩阵的特征值。因此将表1作标准化处理，然后求出8个变量两两之间的相关系数，形成一个8×8阶的方阵，即相关矩阵，见表2。 由matlap[]软件算出各主因子对所有变量的方差贡献，全部方差贡献的总和：，前三个较大的值分别为，。它们在方差总贡献中所占的百分比即累积百分比表3。 由于前两个因子的方差贡献占方差总贡献的89.7％，因此前两个因子已能较好地反映出各变量的变化情况。 表1 新疆某地浸矿山化学成分随时间变化的数据2009.7～2010.6 变量 日期 U/mg/L Ca2+ /mg/L Mg2+/mg/L Fe3+ ∑Fe SO42-/g/L pH Eh 09-7-29 5.67 617.3 128.96 10.7 135.2 9.69 2.35 378 09-8-29 5.535 590.48 128.96 10.2 124.8 9.99 2.44 301 09-9-29 5.67 603.9 120.9 12.5 147.68 10.43 2.15 285 09-10-29 5.67 630.74 96.72 13.1 149.76 10.83 2.48 277 09-11-29 6.3 657.58 88.66 12.8 147.68 11.33 2.42 278 09-12-29 6.3 630.74 128.96 13.3 156 11.77 2.37 266 10-1-29 5.99 590.48 120.9 12.9 151.84 12.16 2.11 261 10-3-1 6.3 563.64 137.02 13.8 166.4 12.42 2.11 256 10-3-29 5.67 523.4 137 14.1 174.72 12.86 1.91 291 10-4-29 6.3 536.8 128.96 15.5 172.6 14.21 1.89 284 10-5-29 6.3 590.48 120.9 14.3 170.56 15.57 1.83 275 10-6-29 5.67 603.9 120.9 16.2 179.92 15.35 1.85 283 表2 8种水化学成分的相关矩阵 2.1435E-05 8.281E-06 3.0709E-05 3.4512E-05 03.951E-05 -2.27E-06 -0.000127 8.281E-06 0-0.0011334 -0.0001569 -0.0006752 -0.0002241 0.0001101 -0.000636 3.0709E-05 -0.0011334 00.00