本发明涉及变电站小电流接地故障判定领域，具体来说是一种基于突变逻辑阵的配网小电流接地故障判定方法。
背景技术：
供电可靠性是一项衡量配电网智能化程度的重要指标，国网“十三五”行动计划也将这一指标纳入了工作报告。
馈线自动化作为提高供电可靠性的重要手段在当前配电网自动化主站系统中已广泛应用，并在日常故障处理和恢复中发挥着重要作用。但是，目前馈线自动化功能大多只能处理馈线短路故障，而对小电流接地故障尚无准确高效的处理手段。小电流接地故障占整个馈线故障的80％以上，提升馈线自动化的小电流接地故障处理能力，是提高供电可靠性的重要手段
目前小电流接地故障定位的方法有两种，一种是无源法，一种是有源法。所谓无源法，是指通过提取故障发生时刻前后的暂态特征量，并经过比对分析从而确定故障位置，如零序电流法、功率方向法、比幅比相法、负序电流法等。该种方法对终端设备采集精度和频率要求很高，需单独安装配置高精度采样装置和抗饱和互感器。所谓有源法，是指通过在接地变压器中性点或者馈线出线处注入特殊信号，通过终端采集到相应信号而综合判定故障位置，如s注入法、加信传递函数法、端口故障诊断法等。该种方法需加装特殊的信号发生器，并需跟特殊的终端配合以检测相应信号。
由上可知，无论是有源法还是无源法皆需要安装特殊的终端设备以配合特殊的采样要求。而当前配电网中已安装大量的采样终端设备如ftu、dtu、故障指示器等，这些设备大都不具备特殊的采样及判定要求。如若安装特殊的采样判定设备，则大部分需要停电安装，进一步降低了供电可靠性，且特殊采样判定设备大部分在初期试制阶段，虽对部分典型线路具有较好的定位效果，但是普遍推广效果欠佳，且价格昂贵，维护繁琐。故此，尽最大可能的利用现有设备和采样数据解决小电流接地故障定位问题，尽量不安装特殊设备，避免重复建设浪费资源，以便更好的建设节约型社会具有重要的现实意义。
技术实现要素：
本发明的目的在于解决变电站小电流接地故障的定位问题，结合配电网特点，给出一种基于突变逻辑阵的配网小电流接地故障判定方法。
用于实现上述目的的技术方案如下：
一种基于突变逻辑阵的配网小电流接地故障判定方法，该方法包括如下步骤：
步骤一、根据配电网设备父子关系生成最小故障判定区域；
步骤二、根据馈线设备三相实时电流数据生成馈线实时电流阵；
步骤三、根据馈线设备三相历史电流数据生成馈线历史电流阵；
步骤四、根据馈线实时电流阵和馈线历史电流阵生成馈线突变电流阵，并根据馈线突变电流阵和馈线对地电容电流阈值生成馈线突变电流逻辑阵；
步骤五、将馈线突变电流逻辑阵非零元素提取构成突变电流逻辑表，根据设备父子关系由父到子的顺序对突变电流逻辑表进行重新排序，得到突变电流逻辑顺序表，并根据突变电流逻辑顺序表进行故障判定。
规定从变电站出线开关开始，沿潮流方向，到各配电变压器为止，某一开关和其所有的子节点设备(开关和配变)构成一个最小故障判定区域。最小故障判定区域是可以进行故障定位的最小范围。具体描述为：
sma＝[ma1,ma2,…,mai,…,man]mai＝[ai1,bi1,bi2,…,bij,…bih]
其中，sma为馈线最小故障判定区域集合，mai为第i个最小故障判定区域，i＝1,2,3,…,n，n为馈线所包含的最小故障判定区域数；ai1为第i个最小故障判定区域的父节点设备，bij为第i个最小故障判定区域中，以ai1为父节点的第j个子节点设备，j＝1,2,3,…,h，h为以ai1为父节点的所有子节点设备数。
定义馈线实时电流阵rcm描述一条配电馈线所有开关和配变的a、b、c三相实时运行状态，具体描述如下：
其中，rak为馈线所包含的第k个开关或配变的a相实时电流值，rbk为馈线所包含的第k个开关或配变的b相实时电流值，rck为馈线所包含的第k个开关或配变的c相实时电流值，k＝1,2,…,m，m为馈线所包含的开关和配变总数。
定义馈线历史电流阵hcm描述一条配电馈线所有开关和配变的a、b、c三相历史运行状态，具体描述如下：
其中，hak、hbk、hck分别为馈线所包含的第k个开关或配变的a、b、c相历史电流值，该历史时刻可根据终端类型设定。
当接地故障发生时，可以用故障发生后的馈线实时电流阵与故障发生前的馈线历史电流阵计算故障引起的稳态电流变化，具体描述为：
acm＝rcm-hcm
其中，acm为馈线突变电流阵，其也是一个3×m矩阵。为计算方便，用馈线突变逻辑阵alm表述馈线电流的突变情况。
其中，lak为馈线第k个设备的a相突变电流逻辑值，lbk为馈线第k个设备的b相突变电流逻辑值，lck为馈线第k个设备的c相突变电流逻辑值；aak、abk、ack分别为馈线突变电流阵acm中第k个设备的a、b、c相电流变化值。α为阈值系数，取值范围为(0,1)。mv为馈线对地电容电流阈值。
若馈线突变逻辑阵alm是零阵，则表示该馈线没有发生接地故障；若alm是非零阵，则表示该馈线发生了接地故障，且非零元素所在行即表示接地故障发生的相别。
若馈线突变逻辑阵alm中的第f(f＝1,2,3)行元素非零，则将所有非零元素提取，构成突变电流逻辑表acl，具体描述为：
acl＝[cl1,cl2,…,cle,…,clg]
其中，cle为alm第f行中第e个逻辑值为1的元素代表的设备编号，e＝1,2,3,…,g，g为alm第f行中非零元素总数。
根据设备父子关系由父到子的顺序对acl进行重新排序，得到突变电流逻辑顺序表acsl，具体描述如下：
acsl＝[cs1,cs2,…,csw,…,csg]
其中，csw为设备编号，w＝1,2,3,…,g，csw-1是设备csw的父节点设备编号，csw+1是设备csw的子节点设备编号。
综合上述分析，接地故障发生在突变电流逻辑顺序表acsl最后一个元素csg为父节点的最小故障判定区间内。
本发明有益技术效果：该方法针对永久性接地故障具有较好的适应性，可依靠当前电网中的传统采集终端实现接地故障的分析定位而不需要安装特殊的采集终端，具有一定的推广价值。
附图说明
图1为典型配电馈线图；
具体实施方式
下面对本发明进行详细的说明：
(1)首先，根据配电网设备父子关系生成最小故障判定区域；
所述的最小故障判定区域是可以进行故障定位的最小范围，规定如下：从变电站出线开关开始沿潮流方向，到各配电变压器为止，某一开关和其所有的子节点设备构成一个最小故障判定区域，所述的的子节点设备包括开关和配电变压器；具体描述为：
sma＝[ma1,ma2,…,mai,…,man]mai＝[ai1,bi1,bi2,…,bij,…bih]
其中，sma为馈线最小故障判定区域集合，mai为第i个最小故障判定区域，i＝1,2,3,…,n，n为馈线所包含的最小故障判定区域数；ai1为第i个最小故障判定区域的父节点设备，bij为第i个最小故障判定区域中，以ai1为父节点的第j个子节点设备，j＝1,2,3,…,h，h为以ai1为父节点的所有子节点设备数。
(2)其次，根据馈线设备三相实时电流数据生成馈线实时电流阵；
定义馈线实时电流阵rcm描述一条配电馈线所有开关和配电变压器的a、b、c三相实时运行状态，具体描述如下：
其中，rak、rbk、rck分别为馈线所包含的第k个开关或配电变压器的a相、b相、c相实时电流值，k＝1,2,…,m，m为馈线所包含的开关和配电变压器总数。
(3)然后，根据馈线设备三相历史电流数据生成馈线历史电流阵；
定义馈线历史电流阵hcm描述一条配电馈线所有开关和配电变压器的a、b、c三相历史运行状态，具体描述如下：
其中，hak、hbk、hck分别为馈线所包含的第k个开关或配电变压器的a相、b相、c相历史电流值，该历史时刻可根据终端类型设定。
(4)再次，根据馈线实时电流阵和馈线历史电流阵生成馈线突变电流阵，并根据馈线突变电流阵和电容电流阈值生成馈线突变电流逻辑阵；
定义馈线突变电流阵acm描述接地故障发生时，用故障发生后的馈线实时电流阵与故障发生前的馈线历史电流阵计算故障引起的稳态电流变化，具体描述为：
acm＝rcm-hcm
其中，acm为馈线突变电流阵，其也是一个3×m矩阵，为计算方便，用馈线突变逻辑阵alm表述馈线电流的突变情况；
其中，lak、lbk、lck为馈线所包含的第k个开关或配电变压器的a相、b相、c相突变电流逻辑；aak、abk、ack分别为馈线突变电流阵acm中馈线所包含的第k个开关或配电变压器的a相、b相、c相电流变化值；α为阈值系数，取值范围为(0,1)；mv为馈线对地电容电流阈值。
(5)最后，将馈线突变电流逻辑阵非零元素提取构成突变电流逻辑表，根据设备父子关系由父到子的顺序对突变电流逻辑表进行重新排序，得到突变电流逻辑顺序表，并根据突变电流逻辑顺序表进行故障判定。
若馈线突变逻辑阵alm是零阵，则表示该馈线没有发生接地故障；若馈线突变逻辑阵alm是非零阵，则表示该馈线发生了接地故障，且非零元素所在行即表示接地故障发生的相别；
若馈线突变逻辑阵alm中的第f(f＝1,2,3)行元素非零，则将所有非零元素提取，构成突变电流逻辑表acl，具体描述为：
acl＝[cl1,cl2,…,cle,…,clg]
其中，cle为alm第f行中第e个逻辑值为1的元素代表的设备编号，e＝1,2,3,…,g，g为alm第f行中非零元素总数；
根据设备父子关系由父到子的顺序对acl进行重新排序，得到突变电流逻辑顺序表acsl，具体描述如下：
acsl＝[cs1,cs2,…,csw,…,csg]
其中，csw为设备编号，w＝1,2,3,…,g，csw-1是设备csw的父节点设备编号，csw+1是设备csw的子节点设备编号；
综合上述分析，接地故障发生在突变电流逻辑顺序表acsl最后一个元素csg为父节点的最小故障判定区间内。
如图1所示配电网典型馈线，#1为变电站母线，s1为变电站出线开关，a～j为安装了ftu、dtu、故障指示器等采集装置的分段开关，dt1～dt5为配电变压器。假设在分段开关d、h、i之间发生了a相永久性接地故障，故事发生后各设备a、b、c三相的实时量测数据ra、rb、rc，故障发生前5分钟的历史量测数据ha、hb、hc如表1所示。
表1故障前后设备量测值(a)
根据(2)～(5)生成馈线突变电流阵acm为：
取馈线对地电容电流阈值mv为10，阈值系数α为0.6，则馈线突变逻辑阵alm为：
因只有馈线突变逻辑阵alm的第一行元素中具有非零值，故接地故障发生在该馈线的a相，根据设备父子关系由父到子的顺序对突变电流逻辑表alm进行重新排序，得到突变电流逻辑顺序表可得突变电流逻辑顺序表acsl为：
acsl＝[s1,a,b,d]
则可知，故障发生在以开关d为父节点的最小故障判定区域内，根据最小故障判定区域，可知发生接地故障的区域为开关d、h、i围成的区域。