[发明专利]基于自适应神经模糊推理系统的配电网故障分类方法

权利要求书

1.一种基于自适应神经模糊推理系统的配电网故障分类方法，其特征在于，包括以下步骤：

步骤1：基于仿真软件，搭建典型配电网结构，仿真出各种类型的配电网故障，建立包括各种配电网故障类型的多个训练样本的训练样本集；构造一个递阶自适应神经模糊推理系统；在使用该推理系统对故障进行分类前，使用混合学习算法对构建的递阶模糊推理系统进行训练，确定系统中的前提参数和结论参数,所述故障类型包括:单相接地短路故障的A相接地短路故障、B相接地短路故障以及C相接地短路故障，两相接地短路故障中的AB两相接地短路故障、AC两相接地短路故障以及BC两相接地短路故障，三相短路故障，两相不接地短路故障中AB两相不接地短路故障、AC两相不接地短路故障以及BC两相不接地短路故障；定义M分别代表10种配电网故障,其中,M＝1,2,…,10；具体包括以下子步骤:

步骤1.1：基于PSCAD-EMTDC仿真，搭建典型配电网结构，在其上对配电网的10种故障类型进行仿真，分别仿真出各类故障的三相电流，以此建立包括各配电网故障类型的多个训练样本的训练样本集；

步骤1.2：对训练样本集中所有训练样本进行特征提取，得到每个训练样本的特征向量，本步骤利用统计参数来构造特征量，对训练样本构造特征量ρ_a,b，ρ_a,c和ρ_b,c为：

其中，

s_max＝max(s_p),p＝a,b,c 式三

其中，分别表示a、b、c三相电流的标幺值；表示为第t个样本的p相电流；r为样本数；E(x)为求变量x的数学期望；s_p(p＝a,b,c)为p相电流的标准差；ρ_a,b、ρ_a,c及ρ_b,c分别表示a相与b相电流的相关系数、a相与c相电流的相关系数以及b相与c相电流的相关系数；

则构造的特征向量表示为[x₁,x₂,x₃,x₄,x₅,x₆]，其中，x₄＝ρ_a,b，x₅＝ρ_a,c，x₆＝ρ_b,c；

步骤1.3：将步骤1.2中构造的训练样本数据的特征向量基于matlab模糊系统工具箱中的ANFIS函数建立包括全部输入变量的TS单级模糊模型；基于所建立的TS单级模糊模型分别计算出每种故障的各特征量相对于输出变量的灵敏度，根据灵敏度分别对所有训练样本的各故障的特征量进行重要性排序；建立的基于TS单级模糊模型的ANFIS是一个典型的多输入单输出的系统，其模糊推理规则为：

式中为规则Rⁱ对应的第j个输入变量的模糊集合，为规则Rⁱ对应的第j个输入变量的结论参数，f_i(X)为系统根据规则Rⁱ所得到的输出，i＝1,2,…,n；j＝1,2,…,m，m表示为输入变量数，n为模糊规则数；

由模糊规则公式，定义系统有m个输入变量x₁,x₂,…,x_m，建立TS单级模糊模型，单输出量由n条模糊If-Then规则组成的集合来表示；所述ANFIS的网络结构共分为5层：隶属度函数生成层、规则推理层、模糊化层、去模糊化层和输出层，由各层得到系统的多输入单输出模型为如下的函数表达形式：

其中，f_i表示为第i个规则对应的系统输出，为第i个规则对应的第j个输入变量的隶属度，分别为第i个规则下第j个输入变量的隶属度函数的中心和宽度，m表示为输入变量数，n为模糊规则数；

所述对所有训练样本的各故障的特征量进行重要性排序具体方法是：标号为M(M＝1,2,…,10)的配电网故障，其训练样本数据表示为[x₁(t),x₂(t),x₃(t),x₄(t),x₅(t),x₆(t)；y(t)],t＝1,2,…,r，其中r为训练样本数，x₁(t),x₂(t),x₃(t),x₄(t),x₅(t),x₆(t)表示为由训练样本根据步骤1.2提取的特征量，y(t)表示为模糊系统的希望输出量，对于标号为M的配电网故障其希望输出量为y(t)＝M，执行以下步骤：

步骤1.3.1，将步骤1.2中构造的训练样本数据的特征向量基于matlab模糊系统工具箱中的ANFIS函数建立包括全部输入变量的TS单级模糊模型；设该单级模糊模型由n条T-S模糊规则进行描述，根据步骤1.3中的式七可得第i条模糊规则的形式为：

Rⁱ：如果x₁为且…，且x₆为则

式中为第i个规则对应的第j个输入变量的模糊集合，为第i个规则对应的第j个输入变量的结论参数；

根据步骤1.3中的式十四可知由n条T-S模糊规则组成的模糊模型的输出为

其中，f_i表示为第i个规则对应的系统输出，为第i个规则对应的第j个输入变量的隶属度，n为模糊规则数，m表示为输入变量数，由步骤1.2可得m＝6；

步骤1.3.2，基于所建立的TS单级模糊模型计算出每一个输入变量相对于输出变量的灵敏度；即第j个输入变量x_j相对于输出变量的灵敏度为：

其中t＝1,2,…,r；为第t个样本的模糊模型输出，x_j(t)为第t个样本的第j个输入变量，f_i(t)表示为第t个样本的第i个规则对应的系统输出，分别为第i个规则对应的第j个输入变量的中心和宽度，为第i个规则对应的第j个输入变量的结论参数，r为训练样本数，n为模糊规则数；

步骤1.3.3，将输入变量x_j根据灵敏度S_j由大到小依次排序，得到X(t)＝[x¹(t)x²(t)x³(t)x⁴(t)x⁵(t)x⁶(t)],t＝1,2,…,r，r为训练样本数；

步骤1.4：构建递阶模糊推理系统，采用模糊C均值聚类方法确定递阶模糊推理系统的规则和初始前提参数；采用递阶模糊推理系统，其是一种双输入二层递阶模糊系统，每个子系统都只有两个输入与一个输出，输出与另一个输入变量共同输入到下一个子系统中，直到最后一个输入变量，根据步骤1.3.3得到的关于输入变量相对于输出变量的灵敏度排序结果对输入变量进行配置，得到新的训练样本数据，由模糊C均值的聚类方法计算出每一种故障类型的聚类中心，执行以下步骤：

步骤1.4.1，初始化；给定聚类类别数10，设定误差值ε和循环次数M′，分别任取各类中的某个样本的特征向量作为初始聚类中心，设置初始循环次数q＝1；

步骤1.4.2，如果q＞M′或者直到||V^(q)-V^(q-1)||＜ε转向步骤1.4.6，否则转向步骤1.4.3，其中V^(q)表示为第q次循环的聚类中心，其中，||·||表示矩阵范数；

步骤1.4.3，更新隶属度矩阵U^(q)：

其中，d_i,j为第j个训练样本与第i个聚类中心的欧式距离，p′为模糊指数，一般p′＝2，u_i,j为第j个训练样本对第i个聚类中心的隶属度值；

步骤1.4.4，更新聚类中心矩阵V^(q+1)：

其中，V^(q+1)表示为第q+1次循环的聚类中心，p′为模糊指数，一般p′＝2，u_i,j为第j个训练样本对第i个聚类中心的隶属度值，r为训练样本数；

步骤1.4.5，令q＝q+1，返回步骤1.4.2；

步骤1.4.6，输出(U，V)；

从而得到10条最优模糊规则和初始模型前提参数和

最优模糊规则为：

R¹：如果x₁为且…，且x_m为则式中为模糊集合；

R²：如果x₁为且…，且x_m为则式中为模糊集合；

R¹⁰：如果x₁为且…，且x_m为则式中为模糊集合；

初始模型前提参数和为：

其中，是子系统第j个输入变量对应的第i个隶属度函数的中心；是子系统的第j个输入变量对应的第i个隶属度函数的宽度；V_ij是第i类聚类中心的第j个元素；θ_i是第i个聚类所包含全部训练样本组成的集合；k_i是第i类集合中的元素的个数；

步骤1.5：使用混合学习算法对构建的递阶模糊推理系统进行训练，确定各子系统最终的前提参数和结论参数；

训练递阶模糊推理系统各级子系统规则前提参数和结论参数，从第一级开始逐级进行参数优化；采用径向基神经网络模型和算法实现前提参数的训练，后件参数用最小二乘法穿插在其中进行，这时是对所有的子系统逐级进行的；

步骤2，当步骤1中构建的递阶模糊推理系统训练完毕后，即开始对测试样本进行分类，具体包括以下子步骤：

步骤2.1:采集故障三相电流，得到测试样本；

在仿真软件上，对典型配电网上设置某一种故障类型，采集其故障三相电流，得到测试样本；

步骤2.2：对测试样本进行特征提取，得到测试样本的特征向量；

采用步骤1中的步骤1.2中的统计参数的方法对测试样本进行处理，得到测试样本的特征向量[x₁,x₂,x₃,x₄,x₅,x₆]；

步骤2.3：按照步骤1中的步骤1.3中的步骤1.3.3中的训练样本的特征量的重要性排序结果，对测试样本的特征量进行排序，得到排序后的特征向量[x¹,x²,x³,x⁴,x⁵,x⁶]；

步骤2.4：将排序后的测试样本的特征向量作为步骤1中的步骤1.5训练后的递阶模糊推理系统的输入量，得到分类结果。

2.根据权利要求1所述一种基于自适应神经模糊推理系统的配电网故障分类方法，其特征在于，所述步骤1.5中，对每一个训练样本有6个特征量，因此需要建立5级子系统；对第j(j＝1,2,…,5)级子系统的训练为例作相应说明，相应的训练数据为[x^j(t),x^j+1(t)；g(t)],j＝1；t＝1,2,…,r或[x^j+1(t),y^j-1(t)；g(t)],j＝2,3,4,5；t＝1,2,…,r，其中r为训练样本数；设置每一级子系统的希望输出为系统最终的希望输出量，更新上一次的该级子系统模型，同时用最小二乘法更新该子系统的结论参数，采用梯度下降的误差反传算法训练前提参数，具体步骤如下：

步骤1.5.1，对步骤1.3中的公式十四作等价变化，将结论参数分离出来：

式中，为第i个规则对应的第j个输入变量的隶属度，x_e为第e个训练样本，n为模糊规则数，m为输入变量数，x₀＝1；进而公式二十简化为：

y＝φ(X)D

其中，

φ(X)=φ01(X)φ11(X)...φm1(X)...φ0n(X)φ1n(X)...φmn(X);]]>

D=b01(X)b11(X)...bm1(X)...b0n(X)b1n(X)...bmn(X)T;]]>

步骤1.5.2，设置误差值ε₁和ε₂和循环次数M″，设置初始循环次数p″＝1；

步骤1.5.3，如果p″＞M″或者且其中表示为第p″次循环下第j个输入变量对应的第i个隶属度函数的中心，表示为第p″次循环下第j个输入变量对应的第i个隶属度函数的宽度，转向步骤1.5.7，否则转向步骤1.5.4；

步骤1.5.4，将r个样本点[x^j+1(t),y^j-1(t)],t＝1,2,…,r输入子系统，系统输出为：

其中φ为1×n(m+1)的矩阵,为由步骤1.5.1的式二十得到的矩阵；

误差指标函数为J(D)＝1/2||Y-φD||²，其中Y为子系统的希望输出矩阵；根据最小二乘法原理，使J(D)达到最小，则可得子系统的结论参数D＝[φ^Tφ]^-1φ^TY；

步骤1.5.5，参数优化的第二步，固定结论参数D，采用梯度下降的误差反传算法训练前提参数；取误差函数为E(c,σ)＝1/2||Y-φD||²，其中Y为子系统的希望输出，D为固定结论参数，φ是关于变量(c,σ)的函数；

则训练的前提参数表示为

步骤1.5.6，令p″＝p″+1，返回步骤1.5.3；

步骤1.5.7，子系统训练结束，此时的前提参数和结论参数为子系统最终的前提参数和结论参数。