CN103278616B - 一种土壤腐蚀性快速评价的多因子方法 - Google Patents

Abstract

本发明公开了一种土壤腐蚀性快速评价的多因子方法，建立评价标准包括对区域内不同土壤样品的信息采集，运用因子分析法从众多土壤腐蚀性因素中筛选出影响土壤腐蚀性的关键因素，运用模糊聚类分析法对土壤腐蚀性等级进行判定。本发明提出了一种将因子分析法和模糊聚类分析相结合的多因子土壤腐蚀性评价方法，该法能够大大减少在评价土壤腐蚀性时所需的工作量以及评价所需时间，同时该方法的评价结果更接近实际，更加具有科学性，为深入研究土壤腐蚀因素间的关系以及对土壤腐蚀性进行准确、快速评价提供新的方法。

Description

一种土壤腐蚀性快速评价的多因子方法

技术领域

本发明属于腐蚀技术领域，具体地说，本发明涉及土壤腐蚀性评价的多因子方法。

背景技术

土壤环境中变电站接地网金属材料受到腐蚀给变电站的建设带来了重大损失，若接地网遭受严重腐蚀,运行中满足不了热稳定要求,容易发生短路造成地网烧断、地电位升高、高压窜入二次回路等事故。美国曾经发生的由于电网腐蚀造成大面积停电事故,已引起世界各国政府和电力企业对电网防腐保护的高度重视。准确评价土壤腐蚀性能够为变电站的建设和维护带来正确的指导，因此评价土壤腐蚀性是一个具有重要现实意义和理论意义的研究课题。但是土壤本身是个由气、液、固三相物质构成的多介质复杂体系，金属在土壤中的腐蚀除了与金属本身的材料有关外，还受到比如土壤的类型、土壤电阻率、含水量、土壤的pH值、可溶性盐含量、土壤的氧化还原电位等土壤理化性能的影响，这些影响因素或单独或相互作用共同影响金属材料在土壤中的腐蚀，评价土壤的腐蚀性是十分困难的，也是研究人员一直在探索的问题。

传统的埋片失重法周期长，短的需要1~2年，长的需要3~5年，费时费力、重现性差，这些特点使它在土壤腐蚀性评价方面不便广泛展开应用，研究较为成熟的是单项指标评价方法，对土壤腐蚀性从土壤电阻率、氧化还原电位、含水量、含盐量、pH值等方面分别进行了评价，但是没有一个土壤因素可单独决定土壤的腐蚀性，土壤环境的复杂性决定了必须探索出更为准确的土壤腐蚀性评价方法。德国的DIN50929标准和美国的ANSIA21.5评价法等多因素评价方法越来越受到研究人员的重视，这些评价方法均先对土壤有关因素分析得出评价，并给出评价指数，然后将这些评价指数累计起来，再给出腐蚀性评价等级，这种打分制的评价方法能够将土壤信息较多的考虑进来，相对于单项指标评价法具有更高的准确性，但此类评价方法的内容针对性过于宽泛，某些参数测量起来较为困难，不利于在实际中广泛应用，同时在该标准中评价标准的测量结果与所评分数间的对应关系存在明显跃迁，郭安祥在陕西省30个不同地区变电站典型土壤的腐蚀性研究中也提出了这些标准的不足之处。近年来，随着数值分析技术的不断发展，越来越多的数值计算方法应用于土壤腐蚀的研究，促进土壤腐蚀性评价方法的不断改进，同时也促进了土壤腐蚀加速试验的发展。

纵观土壤腐蚀性评价的研究现状，对土壤腐蚀性评价缺乏有效地方法，首先从土壤腐蚀性的评价指标来讲，对于土壤腐蚀性评价应该从什么影响因素着手还没有具体的研究，不是对单一的某个指标进行评判就是将全部指标囊括进去，具有一定的盲目性，结果造成两个方面的问题，一是评价不够准确但是方便，另外一个就是评价准确但是不宜开展实施，所以需要提高土壤腐蚀性评价的准确性又要对影响土壤腐蚀性的指标进行筛选，取得关键性因素进行评价，并通过腐蚀加速试验对土壤腐蚀性的评价结果进行验证。

发明内容

本发明的目的在于提供土壤腐蚀性快速评价的多因子方法，以深入研究土壤腐蚀因素间的关系，对土壤腐蚀性进行准确、快速评价，为变电站建设及日常维护提供具有参考价值的建议。

为实现上述目的，本发明采取了以下技术方案：

（1）收集研究区域内不同典型地点土壤样品进行调研分析，统计影响土壤腐蚀性评价准则的各个因素，包括每个土壤样本的各个理化性能数值，以此获取多因子评价方法的信息；

（2）运用因子分析方法从众多的影响因素中确定影响土壤腐蚀性的关键因素，为为多因子评价筛选指标；

（3）通过土壤腐蚀加速试验确定各个土壤样本的腐蚀性等级，为多因子评价方法提供评价基准；

（4）运用模糊聚类分析方法对土壤样本进行聚类分析得到土壤样本的聚类谱系图，根据确定的土壤腐蚀性等级对聚类谱系图进行量化，根据金属在土壤中的腐蚀程度将不同土壤的腐蚀性分为轻微腐蚀、较轻腐蚀、中等腐蚀、较重腐蚀、严重腐蚀五个等级；

（5）引入未知腐蚀性等级的土壤样品的数据信息通过聚类分析确定其土壤腐蚀性等级，从而对评价标准进行评估。

本发明步骤（1）中多因子评价方法的信息来源中所述的土壤信息具体包括：电阻率、氧化还原电位、pH、Na⁺、含水量、含盐量、Cl^-、SO₄ ^2-、CO₃ ^2-、HCO₃ ^-共十个变量

本发明步骤（2）中采用因子分析方法为多因子评价筛选出关键因素，具体步骤如下：

第一步：原始数据的标准化处理

收集到的数据构成一个m×n的矩阵

所选指标都是定量指标，各个指标的含义与计算方法均不同，造成各个指标的量纲也不同，为增加评价结果的科学性和准确性需要对评价指标进行标准化处理，

标准化处理公式：

式中x_ij为原始数据，z_ij为标准化数据

第二步：数据分析前的KMO和Bartlett值检验

进行因子分析需要分析变量间具有较好的相关性，由KMO和Bartlett值来判断能否进行因子分析，其中Bartlett检验用于确认原始变量是否取自于多元正态分布的整体，若数据符合Bartlett检验，则进行KMO检验，KMO检验用于分析原始变量间的偏相关性和简单相关性的相对大小，若KMO过小，则不适合做因子分析。

第三步：对标准化数据求相关矩阵R

得到标准化数据后再进行相应的变换求出相关矩阵，变换公式如下：

i=j时，r_ij=1

i≠j时，r_ij=r_ji

第四步：求相关矩阵的特征值和特征向量

用雅克比算法求得特征值和对应的特征向量

（R-λE）Χ=0

设R的特征跟为λ₁，λ₂，…λ_n，并假定λ₁≥λ₂≥…λ_n≥0，称λ_i为所对应的指标为第i主成分，记trR=λ₁+λ₂+…+λ_n,则λ_i/trR是第i主成分的贡献。记T_i=λ₁+λ₂+…+λ_i,称T_i/trR是前i个主成分的累积贡献。若确定主成分分析的因子累积贡献率为大于85%，则可以对前几个主成分进行提取。

第五步：计算公共因子初始负载及旋转负载

公共因子初始负载的计算方法如下：

其中e_i为特征跟λ_i所对应的特征向量。

由于主成分法确定因子负载比较简单，得到的特殊因子之间并不相互独立，也就是得到的因子负载并不完全符合因子模型的前提条件，系数矩阵B不唯一成为因子负载进行旋转的理论依据，旋转使得因子负载阵结构简化，利于对公共因子进行解释。对初始负载旋转采取最大方差正交旋转法进行旋转。由旋转得到的矩阵可以根据得分判定出影响土壤腐蚀性的关键因素。

本发明步骤（4）采用模糊聚类分析方法开展多因子评价，具体步骤如下：

第一步：选取关键因素后所组成的矩阵，把建立好的标准化矩阵运用求取模糊相似矩阵的传递闭包转化为模糊等价矩阵R。做合成运算M²=M·M，M⁴=M²·M²，…，存在一个自然数k，使得M^2k=M^k·M^k=M^k=（r_ij）_n×n，这时M^k便是一个模糊等价关系,其中2k-1≤n≤2k，k-1<log₂n≤k。

第二步：在建好的模糊等价矩阵的基础上，根据λ值进行聚类，λ截取的规则如下：

参数λ值是设定的聚类模糊度，其取值范围为0-1，得到的矩阵称为λ截距矩阵，该矩阵只有0和1的数值，根据不同的λ值依次求出λ截距矩阵，根据截距矩阵即可进行分类，所有不同的聚类结果形成聚类的谱系图。

本发明步骤（5）中由腐蚀性评价方法判断土壤腐蚀性等级，其具体过程为：以已知腐蚀性的样本为基础得到的谱系图类似于一种评判标准，在此基础上引入新的未知土样，并且为了保证所得到谱系图标准的不变，引进的新样本数目不宜过多，可以每次只引进一个新样，经过与已知腐蚀性等级的样本经过共同聚类分析即可得出聚类结果，从而对土壤的腐蚀性等级进行判定。

与现有技术相比，本发明的土壤腐蚀性评价的多因子方法具有以下的优异性能：

1、本发明通过采用因子分析法对影响土壤腐蚀性的因素进行筛选，提取出关键因素，为土壤腐蚀性评价提供了有针对性且可靠的指标，既节省了传统评价所需的巨大工作量，也避免了进行评价时带有的盲目性；

2、本发明采用了模糊聚类分析方法，对于土壤这种复杂体系，该方法能通过模糊等价关系变换，定量确定各样本的亲疏关系，从而对样本进行科学分类，相比单因子评价在评价中准确性更高；

3、本发明采取土壤腐蚀加速试验箱在进行腐蚀加速试验时更加接近土壤本身的腐蚀状况，并且相对于传统腐蚀试验，节省了大量时间；

4、本发明的土壤腐蚀性评价的多因子方法，因子分析法和模糊聚类分析均可通过计算机软件完成，对于未知土壤的评价也只需土壤关键因素的几项指标，操作方便、能够进行快速评价。

附图说明

图1是本发明总体操作流程图。

图2是本发明中因子分析法的流程图。

具体实施方式

本发明的目的在于提供土壤腐蚀性快速评价的多因子方法，以深入研究土壤腐蚀因素间的关系，对土壤腐蚀性进行准确、快速评价，为变电站建设及日常维护提供具有参考价值的建议。具体步骤如下：

（1）收集研究区域内不同典型地点土壤样品进行调研分析，统计影响土壤腐蚀性评价准则的各个因素，包括每个土壤样本的各个理化性能数值，以此获取多因子评价方法的信息；

（2）运用因子分析方法从众多的影响因素中确定影响土壤腐蚀性的关键因素，为为多因子评价筛选指标；

（3）通过土壤腐蚀加速试验确定各个土壤样本的腐蚀性等级，为多因子评价方法提供评价基准；

（4）运用模糊聚类分析方法对土壤样本进行聚类分析得到土壤样本的聚类谱系图，根据确定的土壤腐蚀性等级对聚类谱系图进行量化，根据金属在土壤中的腐蚀程度将不同土壤的腐蚀性分为轻微腐蚀、较轻腐蚀、中等腐蚀、较重腐蚀、严重腐蚀五个等级；

（5）引入未知腐蚀性等级的土壤样品的数据信息通过聚类分析确定其土壤腐蚀性等级，从而对评价标准进行评估。

本发明步骤（2）中采用因子分析方法为多因子评价筛选出关键因素，具体步骤如下：

第一步：原始数据的标准化处理

收集到的数据构成一个m×n的矩阵

所选指标都是定量指标，各个指标的含义与计算方法均不同，造成各个指标的量纲也不同，为增加评价结果的科学性和准确性需要对评价指标进行标准化处理，

标准化处理公式：

式中x_ij为原始数据，z_ij为标准化数据

第二步：数据分析前的KMO和Bartlett值检验

进行因子分析需要分析变量间具有较好的相关性，由KMO和Bartlett值来判断能否进行因子分析，其中Bartlett检验用于确认原始变量是否取自于多元正态分布的整体，若数据符合Bartlett检验，则进行KMO检验，KMO检验用于分析原始变量间的偏相关性和简单相关性的相对大小，若KMO过小，则不适合做因子分析。

第三步：对标准化数据求相关矩阵R

得到标准化数据后再进行相应的变换求出相关矩阵，变换公式如下：

i=j时，r_ij=1

i≠j时，r_ij=r_ji

第四步：求相关矩阵的特征值和特征向量

用雅克比算法求得特征值和对应的特征向量

（R-λE）Χ=0

设R的特征跟为λ₁，λ₂，…λ_n，并假定λ₁≥λ₂≥…λ_n≥0，称λ_i为所对应的指标为第i主成分，记trR=λ₁+λ₂+…+λ_n,则λ_i/trR是第i主成分的贡献。记T_i=λ₁+λ₂+…+λ_i,称T_i/trR是前i个主成分的累积贡献。若确定主成分分析的因子累积贡献率为大于85%，则可以对前几个主成分进行提取。

第五步：计算公共因子初始负载及旋转负载

公共因子初始负载的计算方法如下：

其中e_i为特征跟λ_i所对应的特征向量。

由于主成分法确定因子负载比较简单，得到的特殊因子之间并不相互独立，也就是得到的因子负载并不完全符合因子模型的前提条件，系数矩阵B不唯一成为因子负载进行旋转的理论依据，旋转使得因子负载阵结构简化，利于对公共因子进行解释。对初始负载旋转采取最大方差正交旋转法进行旋转。由旋转得到的矩阵可以根据得分判定出影响土壤腐蚀性的关键因素。

本发明步骤（4）采用模糊聚类分析方法开展多因子评价，具体步骤如下：

第一步：选取关键因素后所组成的矩阵，把建立好的标准化矩阵运用求取模糊相似矩阵的传递闭包转化为模糊等价矩阵R。做合成运算M²=M·M，M⁴=M²·M²，…，存在一个自然数k，使得M^2k=M^k·M^k=M^k=（r_ij）_n×n，这时M^k便是一个模糊等价关系,其中2k-1≤n≤2k，k-1<log₂n≤k。

第二步：在建好的模糊等价矩阵的基础上，根据λ值进行聚类，λ截取的规则如下：

参数λ值是设定的聚类模糊度，其取值范围为0-1，得到的矩阵称为λ截距矩阵，该矩阵只有0和1的数值，根据不同的λ值依次求出λ截距矩阵，根据截距矩阵即可进行分类，所有不同的聚类结果形成聚类的谱系图。

本发明步骤（5）中由腐蚀性评价方法判断土壤腐蚀性等级，其具体过程为：以已知腐蚀性的样本为基础得到的谱系图类似于一种评判标准，在此基础上引入新的未知土样，并且为了保证所得到谱系图标准的不变，引进的新样本数目不宜过多，可以每次只引进一个新样，经过与已知腐蚀性等级的样本经过共同聚类分析即可得出聚类结果，从而对土壤的腐蚀性等级进行判定。

以上是针对本发明的可行实施例的具体说明，但该实施例并非仅用以限制本发明的专利范围，凡未脱离本发明的等效实施或变更，均应包含于本发明的专利范围中。

Claims (1)

1.一种土壤腐蚀性快速评价的多因子方法，其特征是在方法上建立土壤腐蚀性快速评价准则，包括如下步骤：

（1）统计影响土壤腐蚀性评价准则的各个因素，包括每个土壤样本的各个理化性能数值，获取多因子评价方法的信息来源；

（2）运用因子分析方法从众多的影响因素中确定影响土壤腐蚀性的关键因素，确定多因子评价的指标；

（3）运用模糊聚类分析方法对土壤样本所筛选的关键因素进行聚类分析得到土壤样本的聚类谱系图，根据确定的土壤腐蚀性等级对聚类谱系图进行标准化，将不同土壤的腐蚀性分为轻微腐蚀、较轻腐蚀、中等腐蚀、较重腐蚀、严重腐蚀五个等级；

（4）根据未知腐蚀性等级的土壤样品的数据信息通过聚类分析确定其处于聚类谱系图中的位置从而得到其土壤腐蚀性等级；

所述的步骤（1）中影响土壤腐蚀性评价准则的各个因素具体包括：土壤电阻率、土壤氧化还原电位、pH、pNa、含水量、含盐量、Cl^-、SO₄ ^2-、CO₃ ^2- _、HCO₃ ^-；

所述的步骤（2），因子分析法具体步骤如下：

第一步：原始数据的标准化处理：

从m个土壤样本检测得到n项指标的数据构成一个m×n的矩阵

Х_m×n=

标准化处理公式：

式中x_ij为原始数据，z_ij为标准化数据

第二步：数据分析前的KMO和Bartlett值检验；

第三步：对标准化数据求相关矩阵R：

r_ij

i=j时，r_ij=1

i≠j时，r_ij=r_ji

第四步：求相关矩阵的特征值和特征向量：

（R-λE）Χ=0

第五步：计算公共因子初始负载及旋转负载：

公共因子初始负载的计算方法如下：

B=

对初始负载旋转采取最大方差正交旋转法进行旋转；

所述的步骤（3），模糊聚类分析具体分为以下步骤：

第一步：选取关键因素后所组成的矩阵，进行标准化处理，将建立好的标准化矩阵转化为模糊等价矩阵R，转化的方法是求取模糊相似矩阵的传递闭包；

第二步：在建好的模糊等价矩阵的基础上，根据λ值进行聚类，λ截取的规则如下：

根据不同的λ值依次求出λ截距矩阵，所有不同的聚类结果形成聚类的谱系图；

所述的步骤（4），在引入未知土壤样本进行聚类分析的时候，每次只引进一个新样本。