本发明属于感知方法，具体涉及一种基于分布式实时数据仓库的台区停电及故障实时感知方法。
背景技术：
：电网在生产生活中扮演着越来越重要的角色。随着科技的发展，我国的电网由原来的单发电机组电网逐渐向分布式电网发展。分布式电网具有多个独立的发电设备，每个分布式网络内也分布着不同的用电设备。这种电网的缺陷在于：上级电网往往只收集下级电网上报的相关数据，而不关注下级电网内部运行的数据。如果发生事故，例如停电，跳闸，输电线路损毁等，那么剧烈的电网数据变动可以体现在下级电网上报的数据中。如果只是下级电网的输电质量产生了波动，例如电压波动、三相电压不平衡、谐波所占比例过高、频率偏差等，而这些输电质量问题尚未构成电网事故，那么由于子电网的“内部消化”作用，则上级电网可能完全无法了解下级电网是否产生了故障。而且由于电网互相是连通的，一个子电网产生的故障很可能会波及周围电网，因此故障范围大小也是需要感知的内容。现有技术尚无法解决上述问题，因此需要一种针对分布式电网电能质量进行感知的方法。技术实现要素：本发明针对现有技术的缺陷，提供一种基于分布式实时数据仓库的台区停电及故障实时感知方法。本发明是这样实现的：一种基于分布式实时数据仓库的台区停电及故障实时感知方法，包括下述步骤：步骤一：构建数据矩阵；步骤二：采样或接受上报数据；步骤三：核算标志位形成完整数据矩阵；步骤四：事故感知；步骤五：故障感知。如上所述的一种基于分布式实时数据仓库的台区停电及故障实时感知方法，其中，所述的步骤一包括，构建下述结构的数据矩阵z＝[aij]上述数据矩阵的行有无限多，列共有10列；每列数据的物理意义分别为，第一列ai1表示采样或者上报数据的次数；第二列ai2表示采样或上报的节点编号；第三列ai3表示谐波所占比例标志位；第四列ai4表示电压波动比例标志位；第五列ai5表示三相电压不平衡标志位；第六列ai6表示频率偏差值标志位；第七列ai7表示谐波所占比例采样或上报数据值；第八列ai8表示电压波动比例采样或上报数据值；第九列ai9表示三相电压不平衡采样或上报数据值；第十列ai10表示频率偏差值采样或上报数据值。如上所述的一种基于分布式实时数据仓库的台区停电及故障实时感知方法，其中，所述的步骤二包括，当下层数据为基层电网数据时，进行采样；当下层数据为非基层电网数据时，接受下级电网的上报数据，无论采集数据还是上报的数据，均符合下述要求，采样或上报的数据包括下述四种，即分别针对谐波所占比例、电压波动比例、三相电压不平衡、频率偏差值进行采样，采样频率根据需要进行，以电网频率的10-20倍频率为采样频率，上报频率与采样频率相同，本申请的采样或上报内容都不是直接取得所需数据，而是采集相关基础数据然后根据基础数据计算得到的，具体为：针对谐波所占比例：对电网电能采样，计算谐波所占比例，该比例用百分比表示，该所占比例为采样数据，电压波动比例：采集电网电压，将电网电压与标准值进行相减并取绝对值，然后将差值与标准值相除，得到百分比数值，得到的百分比数值为本申请采样的内容，三相电压不平衡：用电压负序分量与正序分量的方均根，并百分比表示，该百分比值为本申请的采样内容，频率偏差：采集电网频率，用该频率与标准频率的差作为频率偏差的采样内容，由于本申请只需要使用其数值部分，因此所述的采样内容为频率偏差的绝对值。如上所述的一种基于分布式实时数据仓库的台区停电及故障实时感知方法，其中，所述的步骤三包括，设定阈值，所述的阈值，按照下述标准设定：所述的谐波所占比例的阈值为5％，所述的电压波动比例的阈值为5％，所述的三相电压不平衡的阈值2％，所述的频率偏差值的阈值为0.2hz，根据步骤二中实际采样或上报数据的值，以上述阈值为标准，若实际采样或上报数据的值小于等于阈值，则置相应标志位为1，否则置相应标志位为0。如上所述的一种基于分布式实时数据仓库的台区停电及故障实时感知方法，其中，所述的步骤四包括，步骤4.1：形成事故感知矩阵取数据矩阵中的前6列数据形成事故感知矩阵，在事故感知矩阵中，第一列数据ai1的数值相同，即取出同次采样或上报数据的不同节点的所有标志位数据，形成事故感知矩阵，步骤4.2：感知事故范围取事故矩阵后四列，如果同列中有任意两行的数据均为0，那么该两行对应节点的特定事故相关性增加1，并判定有事故，否则相关性保持不变，步骤4.3：输出事故范围将步骤4.2感知到的结果输出，如果步骤4.2感知到的结果为无事故，那么执行步骤五；如果步骤4.2感知到的结果为有事故，则输出事故类型和事故范围。如上所述的一种基于分布式实时数据仓库的台区停电及故障实时感知方法，其中，所述的步骤五包括，步骤5.1：形成故障感知矩阵取数据矩阵中的前2列数据和后4列数据形成故障感知矩阵，在事故障知矩阵中，第一列数据ai1的数值相同，即取出同次采样或上报数据的不同节点的所有采样或上报数据，形成故障感知矩阵，步骤5.2：感知事故范围由于故障相关性的计算是对称的，即a节点与b节点的故障相关性，与b节点与a节点的故障相关性是相同的，因此本步骤后续计算中只计算a→b的故障相关性，不重复计算b→a的故障相关性，步骤5.2.1：谐波所占比例故障相关性对于第三列数据，若第a行第3列的元素ca3的值，和第b行第3列的元素cb3的值均为0时，判定无故障，否则判定有故障，并用下述公式计算其中，a、b表示感知故障矩阵的行数，且a≠b，表示感知故障矩阵中第a行第3列的元素ca3的值的平方，表示感知故障矩阵中第b行第3列的元素cb3的值的平方；||表示取绝对值，若的值小于等于0.25，则判定第a节点和第b节点的谐波所占比例故障相关性存在故障相关性，步骤5.2.2：电压波动比例故障相关性对于第四列数据，若第a行第4列的元素ca4的值，和第b行第4列的元素cb4的值均为0时，判定无故障；否则判定有故障；当判定结果为有故障时，如果ca4和cb4中有任意一个的值为0而另一个不为0，那么令当ca4和cb4均不为0时，并用下述公式计算其中ca4表示第a行第4列的元素的值；cb4第b行第4列的元素的值，若的值小于等于0.2，则判定第a节点和第b节点的电压波动比例存在故障相关性，步骤5.2.3：三相电压不平衡故障相关性对于第五列数据，若第a行第5列的元素ca5的值，和第b行第5列的元素cb5的值均为0时，判定无故障，否则判定有故障，并用下述公式计算其中ca5表示第a行第5列的元素的值；cb5第b行第5列的元素的值，若的值小于等于0.25，则判定第a节点和第b节点的三相电压不平衡存在故障相关性，步骤5.2.4：频率偏差值故障相关性对于第六列数据，若第a行第6列的元素ca6的值，和第b行第6列的元素cb6的值均为0时，判定无故障，否则判定有故障，并用下述公式计算其中ca6表示第a行第6列的元素的值；cb6第b行第6列的元素的值，||表示取绝对值，若的值大于等于0.8，则判定第a节点和第b节点的三频率偏差值存在故障相关性，步骤5.3：故障范围输出若步骤5.2判定的结果为无故障，则输出无故障；若步骤5.2判定的结果为有故障，则输出故障类型及故障范围。本发明的显著效果是：本发明针对特定电网参数进行收集或上报，并在收集或上报的时候就进行事故判断，节约了大量运算工作量。在事故和故障两个层面，分别感知影响范围，不仅能够了解分布式电网事故或故障类型，同时也能感知影响范围。具体实施方式一种基于分布式实时数据仓库的台区停电及故障实时感知方法，包括下述内容：步骤一：构建数据矩阵构建下述结构的数据矩阵z＝[aij]上述数据矩阵的行有无限多，列共有10列；每列数据的物理意义分别为，第一列ai1表示采样或者上报数据的次数；第二列ai2表示采样或上报的节点编号；第三列ai3表示谐波所占比例标志位；第四列ai4表示电压波动比例标志位；第五列ai5表示三相电压不平衡标志位；第六列ai6表示频率偏差值标志位；第七列ai7表示谐波所占比例采样或上报数据值；第八列ai8表示电压波动比例采样或上报数据值；第九列ai9表示三相电压不平衡采样或上报数据值；第十列ai10表示频率偏差值采样或上报数据值。上述数据矩阵之所以有无限多行，是因为随着时间的推移，子节点的采样或者上报数据会不断的填写到该数据矩阵中，因此该数据矩阵理论上有无限多行。在实际操作中，由于存储空间的限制，一般采用24小时备份制度，即24小时内的数据直接存储在运算计算机中，超过24小时的数据备份到服务器的备份硬盘中。3个月以上的数据根据实际需要，可以考虑压缩后归档存储，也可以考虑直接丢弃该部分数据。步骤二：采样或接受上报数据当下层数据为基层电网数据时，进行采样；当下层数据为非基层电网数据时，接受下级电网的上报数据。无论采集数据还是上报的数据，均符合下述要求。采样或上报的数据包括下述四种，即分别针对谐波所占比例、电压波动比例、三相电压不平衡、频率偏差值进行采样。采样频率根据需要进行，一般以电网频率的10-20倍频率为采样频率。虽然采样频率可以根据需要调整，但是只要是同一次评价活动，该采样频率一旦确定就不再更改。上报频率与采样频率相同。本申请的采样或上报内容都不是直接取得所需数据，而是采集相关基础数据然后根据基础数据计算得到的，具体为针对谐波所占比例：对电网电能采样，计算谐波所占比例，该比例用百分比表示，该所占比例为采样数据。电压波动比例：采集电网电压，将电网电压与标准值进行相减并取绝对值，然后将差值与标准值相除，得到百分比数值，得到的百分比数值为本申请采样的内容。三相电压不平衡：用电压负序分量与正序分量的方均根，并百分比表示，该百分比值为本申请的采样内容。频率偏差：采集电网频率，用该频率与标准频率的差作为频率偏差的采样内容。由于本申请只需要使用其数值部分，因此所述的采样内容为频率偏差的绝对值。步骤三：核算标志位形成完整数据矩阵设定阈值，所述的阈值，按照下述标准设定：所述的谐波所占比例的阈值为5％，所述的电压波动比例的阈值为5％，所述的三相电压不平衡的阈值2％，所述的频率偏差值的阈值为0.2hz。根据步骤二中实际采样或上报数据的值，以上述阈值为标准，若实际采样或上报数据的值小于等于阈值，则置相应标志位为1，否则置相应标志位为0。例如某节点在某时刻上报数据如下表所示谐波所占比例电压波动比例三相电压不平衡频率偏差值3％6％2.5％0.5hz那么根据上述阈值判断标准，则相应标志位的值分别为假设该节点为分布式系统的第5个子节点，上报数据为第384次上报数据，那么形成的该条数据矩阵为384510003％6％2.5％0.5hz步骤四：事故感知步骤4.1：形成事故感知矩阵取数据矩阵中的前6列数据形成事故感知矩阵。在事故感知矩阵中，第一列数据ai1的数值相同。即取出同次采样或上报数据的不同节点的所有标志位数据，形成事故感知矩阵。例如数据矩阵的一部分如下表所示384111113％2％1％0.1hz384210113％6％1％0.1hz384300116％6％1％0.1hz384411113％2％1％0.1hz384510003％6％2.5％0.5hz385111111％1％1％0.1hz385211111％1％1％0.1hz385311111％1％1％0.1hz385411111％2％1％0.1hz385511111％1％1％0.1hz386111111％1％1％0.1hz386211011％1％3％0.1hz386311011％1％3％0.1hz386410111％6％1％0.1hz取384次上报数据形成事故感知矩阵，则得到的384次事故感知矩阵如下所示3841111138421011384300113844111138451000步骤4.2：感知事故范围取事故矩阵后四列，如果同列中有任意两行的数据均为0，那么该两行对应节点的特定事故相关性增加1，并判定有事故，否则相关性保持不变。仍以384次上报数据形成事故感知矩阵为例，该事故感知矩阵为3841111138421011384300113844111138451000第三列只有1个数据0，则所有节点的相关性均保持不变；第四列有3个数据0，则节点2、节点3和节点5，三个节点中任意两个的电压波动比例事故具备相关性；第五列、第六列与第三列类似，也只有1个数据0，所有节点的相关性均保持不变。如果事故感知矩阵中后四列的数据均为1，则判定无事故。步骤4.3：输出事故范围将步骤4.2感知到的结果输出。如果步骤4.2感知到的结果为无事故，那么执行步骤五；如果步骤4.2感知到的结果为有事故，则输出事故类型和事故范围。仍以384次上报数据形成事故感知矩阵为例，判定结果为有事故，事故类型为电压波动比例事故，事故范围为节点2、节点3和节点5。步骤五：故障感知步骤5.1：形成故障感知矩阵取数据矩阵中的前2列数据和后4列数据形成故障感知矩阵。在事故障知矩阵中，第一列数据ai1的数值相同。即取出同次采样或上报数据的不同节点的所有采样或上报数据，形成故障感知矩阵。例如数据矩阵的一部分如下表所示取385次上报数据形成故障感知矩阵，则得到的385次故障感知矩阵如下所示38511％1％1％0.1hz38521％1％1％0.1hz38531％1％1％0.1hz38541％2％1％0.1hz38551％1％1％0.1hz步骤5.2：感知事故范围由于故障相关性的计算是对称的，即a节点与b节点的故障相关性，与b节点与a节点的故障相关性是相同的，因此本步骤后续计算中只计算a→b的故障相关性，不重复计算b→a的故障相关性。步骤5.2.1：谐波所占比例故障相关性对于第三列数据，若第a行第3列的元素ca3的值，和第b行第3列的元素cb3的值均为0时，判定无故障，否则判定有故障，并用下述公式计算其中，a、b表示感知故障矩阵的行数，且a≠b，表示感知故障矩阵中第a行第3列的元素ca3的值的平方，表示感知故障矩阵中第b行第3列的元素cb3的值的平方；||表示取绝对值。若的值小于等于0.25，则判定第a节点和第b节点的谐波所占比例故障相关性存在故障相关性。步骤5.2.2：电压波动比例故障相关性对于第四列数据，若第a行第4列的元素ca4的值，和第b行第4列的元素cb4的值均为0时，判定无故障；否则判定有故障；当判定结果为有故障时，如果ca4和cb4中有任意一个的值为0而另一个不为0，那么令当ca4和cb4均不为0时，并用下述公式计算其中ca4表示第a行第4列的元素的值；cb4第b行第4列的元素的值。若的值小于等于0.2，则判定第a节点和第b节点的电压波动比例存在故障相关性。步骤5.2.3：三相电压不平衡故障相关性对于第五列数据，若第a行第5列的元素ca5的值，和第b行第5列的元素cb5的值均为0时，判定无故障，否则判定有故障，并用下述公式计算其中ca5表示第a行第5列的元素的值；cb5第b行第5列的元素的值，若的值小于等于0.25，则判定第a节点和第b节点的三相电压不平衡存在故障相关性。步骤5.2.4：频率偏差值故障相关性对于第六列数据，若第a行第6列的元素ca6的值，和第b行第6列的元素cb6的值均为0时，判定无故障，否则判定有故障，并用下述公式计算其中ca6表示第a行第6列的元素的值；cb6第b行第6列的元素的值，||表示取绝对值。若的值大于等于0.8，则判定第a节点和第b节点的三频率偏差值存在故障相关性。步骤5.3：故障范围输出若步骤5.2判定的结果为无故障，则输出无故障；若步骤5.2判定的结果为有故障，则输出故障类型及故障范围。例如前述例子中，存在故障，故障分别为谐波所占比例故障，故障范围为节点1、节点2、节点3、节点4、节点5；电压波动比例故障，故障范围为节点1、节点2、节点3、节点4、节点5；三相电压不平衡故障，故障范围为节点1、节点2、节点3、节点4、节点5；频率偏差值故障，故障范围为节点1、节点2、节点3、节点4、节点5。当前第1页1 2 3
技术特征：1.一种基于分布式实时数据仓库的台区停电及故障实时感知方法，其特征在于，包括下述步骤：
步骤一：构建数据矩阵；
步骤二：采样或接受上报数据；
步骤三：核算标志位形成完整数据矩阵；
步骤四：事故感知；
步骤五：故障感知。
2.如权利要求1所述的一种基于分布式实时数据仓库的台区停电及故障实时感知方法，其特征在于：所述的步骤一包括，
构建下述结构的数据矩阵
z＝[aij]
上述数据矩阵的行有无限多，列共有10列；每列数据的物理意义分别为，
第一列ai1表示采样或者上报数据的次数；
第二列ai2表示采样或上报的节点编号；
第三列ai3表示谐波所占比例标志位；
第四列ai4表示电压波动比例标志位；
第五列ai5表示三相电压不平衡标志位；
第六列ai6表示频率偏差值标志位；
第七列ai7表示谐波所占比例采样或上报数据值；
第八列ai8表示电压波动比例采样或上报数据值；
第九列ai9表示三相电压不平衡采样或上报数据值；
第十列ai10表示频率偏差值采样或上报数据值。
3.如权利要求2所述的一种基于分布式实时数据仓库的台区停电及故障实时感知方法，其特征在于：所述的步骤二包括，
当下层数据为基层电网数据时，进行采样；当下层数据为非基层电网数据时，接受下级电网的上报数据，无论采集数据还是上报的数据，均符合下述要求，
采样或上报的数据包括下述四种，即分别针对谐波所占比例、电压波动比例、三相电压不平衡、频率偏差值进行采样，
采样频率根据需要进行，以电网频率的10-20倍频率为采样频率，
上报频率与采样频率相同，
本申请的采样或上报内容都不是直接取得所需数据，而是采集相关基础数据然后根据基础数据计算得到的，具体为：
针对谐波所占比例：对电网电能采样，计算谐波所占比例，该比例用百分比表示，该所占比例为采样数据，
电压波动比例：采集电网电压，将电网电压与标准值进行相减并取绝对值，然后将差值与标准值相除，得到百分比数值，得到的百分比数值为本申请采样的内容，
三相电压不平衡：用电压负序分量与正序分量的方均根，并百分比表示，该百分比值为本申请的采样内容，
频率偏差：采集电网频率，用该频率与标准频率的差作为频率偏差的采样内容，由于本申请只需要使用其数值部分，因此所述的采样内容为频率偏差的绝对值。
4.如权利要求3所述的一种基于分布式实时数据仓库的台区停电及故障实时感知方法，其特征在于：所述的步骤三包括，
设定阈值，所述的阈值，按照下述标准设定：
所述的谐波所占比例的阈值为5％，
所述的电压波动比例的阈值为5％，
所述的三相电压不平衡的阈值2％，
所述的频率偏差值的阈值为0.2hz，
根据步骤二中实际采样或上报数据的值，以上述阈值为标准，若实际采样或上报数据的值小于等于阈值，则置相应标志位为1，否则置相应标志位为0。
5.如权利要求4所述的一种基于分布式实时数据仓库的台区停电及故障实时感知方法，其特征在于：所述的步骤四包括，
步骤4.1：形成事故感知矩阵
取数据矩阵中的前6列数据形成事故感知矩阵，在事故感知矩阵中，第一列数据ai1的数值相同，即取出同次采样或上报数据的不同节点的所有标志位数据，形成事故感知矩阵，
步骤4.2：感知事故范围
取事故矩阵后四列，如果同列中有任意两行的数据均为0，那么该两行对应节点的特定事故相关性增加1，并判定有事故，否则相关性保持不变，
步骤4.3：输出事故范围
将步骤4.2感知到的结果输出，
如果步骤4.2感知到的结果为无事故，那么执行步骤五；
如果步骤4.2感知到的结果为有事故，则输出事故类型和事故范围。
6.如权利要求5所述的一种基于分布式实时数据仓库的台区停电及故障实时感知方法，其特征在于：所述的步骤五包括，
步骤5.1：形成故障感知矩阵
取数据矩阵中的前2列数据和后4列数据形成故障感知矩阵，在事故障知矩阵中，第一列数据ai1的数值相同，即取出同次采样或上报数据的不同节点的所有采样或上报数据，形成故障感知矩阵，
步骤5.2：感知事故范围
由于故障相关性的计算是对称的，即a节点与b节点的故障相关性，与b节点与a节点的故障相关性是相同的，因此本步骤后续计算中只计算a→b的故障相关性，不重复计算b→a的故障相关性，
步骤5.2.1：谐波所占比例故障相关性
对于第三列数据，若第a行第3列的元素ca3的值，和第b行第3列的元素cb3的值均为0时，判定无故障，否则判定有故障，并用下述公式计算
其中，a、b表示感知故障矩阵的行数，且a≠b，表示感知故障矩阵中第a行第3列的元素ca3的值的平方，表示感知故障矩阵中第b行第3列的元素cb3的值的平方；||表示取绝对值，
若的值小于等于0.25，则判定第a节点和第b节点的谐波所占比例故障相关性存在故障相关性，
步骤5.2.2：电压波动比例故障相关性
对于第四列数据，若第a行第4列的元素ca4的值，和第b行第4列的元素cb4的值均为0时，判定无故障；否则判定有故障；当判定结果为有故障时，如果ca4和cb4中有任意一个的值为0而另一个不为0，那么令当ca4和cb4均不为0时，并用下述公式计算
其中ca4表示第a行第4列的元素的值；cb4第b行第4列的元素的值，
若的值小于等于0.2，则判定第a节点和第b节点的电压波动比例存在故障相关性，
步骤5.2.3：三相电压不平衡故障相关性
对于第五列数据，若第a行第5列的元素ca5的值，和第b行第5列的元素cb5的值均为0时，判定无故障，否则判定有故障，并用下述公式计算
其中ca5表示第a行第5列的元素的值；cb5第b行第5列的元素的值，
若的值小于等于0.25，则判定第a节点和第b节点的三相电压不平衡存在故障相关性，
步骤5.2.4：频率偏差值故障相关性
对于第六列数据，若第a行第6列的元素ca6的值，和第b行第6列的元素cb6的值均为0时，判定无故障，否则判定有故障，并用下述公式计算
其中ca6表示第a行第6列的元素的值；cb6第b行第6列的元素的值，||表示取绝对值，
若的值大于等于0.8，则判定第a节点和第b节点的三频率偏差值存在故障相关性，
步骤5.3：故障范围输出
若步骤5.2判定的结果为无故障，则输出无故障；
若步骤5.2判定的结果为有故障，则输出故障类型及故障范围。
技术总结本发明属于感知方法，具体涉及一种基于分布式实时数据仓库的台区停电及故障实时感知方法。一种基于分布式实时数据仓库的台区停电及故障实时感知方法，包括下述步骤：步骤一：构建数据矩阵；步骤二：采样或接受上报数据；步骤三：核算标志位形成完整数据矩阵；步骤四：事故感知；步骤五：故障感知。本发明的显著效果是：本发明针对特定电网参数进行收集或上报，并在收集或上报的时候就进行事故判断，节约了大量运算工作量。在事故和故障两个层面，分别感知影响范围，不仅能够了解分布式电网事故或故障类型，同时也能感知影响范围。
技术研发人员：孙伟;郑高峰;张淑娟;汪玉;卞真旭;邵珺伟;唐亮;梁晓伟;姜海辉;杨睿;张良;陶晓峰
受保护的技术使用者：国网安徽省电力有限公司电力科学研究院;国网安徽省电力有限公司;北明软件有限公司
技术研发日：2020.11.17
技术公布日：2021.03.12
