一种用于广义流量仪表运行误差远程校准的虚拟支路修正拓扑模型的制作方法

本发明属于流量表计量领域，尤其是一种用于广义流量仪表运行误差远程校准的虚拟支路修正拓扑模型。

背景技术：

运行中的电能计量装置准确性始终是居民用户和国网公司最关注的热点问题。目前，国网公司采用统计分析方法或者添加在线检测设备等手段，来实现智能电能表运行误差的监督与评价。虽然统计分析法改变了过去居民用单相电能表由定期轮换为抽检，但仍是对整体批次的控制；通过添加在线监测设备必将为企业带来采购、维护成本的增加，该方法维护成本高且有局限性。

在实际运行中实际的电能表集群存在损耗，目前，国网公司采用的修正方法是引入定值的线损，依据能量守恒原理处理线损及损耗，其计算结果较为确定，但是实际运行中的耗损为变化值，现有方法存在偏差较大。

本专利针对国网公司4.8亿具在线运行的智能电能表运用全数据样本思维方式，挖掘电力用户用电信息采集系统海量电能数据，无需借助外部标准仪器，开展智能电能表远程校验技术研究，使得供用电双方都能实时了解电能计量装置的准确性，满足客户期望性需求，提高企业形象。

经过检索，未发现相近技术的已公开专利文献。

技术实现要素：

本发明的目的在于克服现有技术的不足之处，提供一种用于广义流量仪表运行误差远程校准的虚拟支路修正拓扑模型。

本发明解决其技术问题是采取以下技术方案实现的：

一种用于广义流量仪表运行误差远程校准的虚拟支路修正拓扑模型，其特征在于：

⑴广义流量仪表集群安装形成树形拓扑结构；

⑵引入虚拟支路，修正广义流量仪表集群拓扑模型；

⑶根据修正后的广义流量仪表树形拓扑结构，以及流量守恒基础模型，获得本集群的流量守恒算法模型

α为修正经验。

而且，所述广义流量仪表包括电能表、水表、燃气表以及热量表。

而且，所述虚拟支路代替线路耗损的耗损值之和，该虚拟支路包括虚拟广义流量仪表和虚拟负载。

而且，所述流量守恒基础模型为：yε＝-η，其中，

yi为仪表m1～mn-1第i次测量的结果向量，

ε＝(ε1ε2…εn-1)^t，η＝(y1,0ε0y2,0ε0…ym,0ε0)^t；

输入任一仪表的相对误差δ0，则有可能通过所有仪表的多次测量结果推算出其他各仪表的相对误差δj,j＝1,2,…,n-1。

而且，所述虚拟支路的流出为集群总表m0的读数增量减去各个分表mj(j＝1,…,n)的和，并将这个值乘以一个由经验值判定的损耗系数；加入虚拟支路后，方程相应变化为α为修正经验。

本发明的优点和积极效果是：

1、本发明利用流量守恒原理和大数据分析技术，在不添加标准设备、不改变表计结构、不改变电能表集群拓扑的情况下，实现海量在线运行电能表误差的远程检测。

2、本发明建立基于虚拟支路的修正拓扑模型，线损为变化值，关注计算结果及变化，，实现在线运行电能表误差测量可靠性100％，更精益化地指导电能表更换工作。

3、本技术不仅可以实现居民家庭电能表运行误差的全面远程检测与诊断，还可以推广到全部广义流量计量，如天然气、水、热等领域，应用前景广阔。

附图说明

图1为本发明误差远程校验流程；

图2为树形拓扑下广义流量仪表集群示意图；

图3为引入虚拟支路的电能表集群示意图。

具体实施方式

下面结合附图并通过具体实施例对本发明作进一步详述，以下实施例只是描述性的，不是限定性的，不能以此限定本发明的保护范围。

为了便于清楚描述本方法，首先给出如下定义：

流：流入或流出某一封闭曲面s的标量(图2箭头所示)。流入为正，流出为负。对流做如下约定：流不会堆积且任何时刻都遵守守恒定律。即穿过任意封闭曲面s的流代数和等于零。流以从0开始的整数标记。

广义流量：简称流量，流对时间的积分。由流守恒可以推出流量守恒。即任意时间段内，穿过封闭曲面s的流量增量代数和等于零。

广义流量仪表：简称流量仪表或流量计(图2中m0～m4)，记录对应流的瞬时流量。流量计以mj标记，j为对应流的编号。

集群：穿过某一封闭曲面s的所有流所对应的流量仪表的集合。

抄表：获取某一时刻集群内所有仪表的读数。逻辑上，抄表动作应返回各仪表相同时刻的读数。但实际上，若各仪表读数返回的时间间隔足够短，则可认为是相同时刻。

测量：一次测量是对集群间隔一段时间的两次抄表，其结果是前后两次读数的增量。

一种用于广义流量仪表运行误差远程校准的虚拟支路修正拓扑模型，本方法依靠日冻结电量，通过算法计算电能表误差，实现对电能表误差的实时远程监控。利用电能表采集系统抄表数据对电能表进行远程校准。

本方法将电能视为一种满足守恒定律的广义流量，从而将电能表与水表、气表等常规流量仪表抽象为一类广义流量仪表。此类仪表具有的共同特征是安装后形成的集群具有树形拓扑。因此在流量守恒约定下，同一时段内总表的实际流量增量等于各分表的实际流量增量之和。由于实际流量增量可以读数增量与相对误差表示，因此可以得到一个包含所有仪表读数增量与相对误差的方程。如果已知集群中任意一个仪表的相对误差而将其他仪表的相对误差作为未知量，并且注意到通过增加测量时段可以使方程达到或超过未知数的数量，那么可以通过求解方程组确定其他个仪表的相对误差。该自主式算法通过集群内部仪表相互比对计算误差，无需借助外部标准仪器。

当树形拓扑流管存在泄漏时，必须对方程组进行修正。在电能计量中，泄漏表现为线路漏电、电阻损耗以及电表功耗。在本方法中，损耗被视为虚拟支路的负载，并引入虚拟电能表。利用智能电能表的电参数测量功能，确定损耗修正项。给出修正后的算法流程以及仿真结果。

本方法不仅可以应用在智能电表，可以应用智能水表、智能燃气表以及智能热量表等广义流量仪表的远程校准当中。

1.理想情况下的流量守恒：

相对误差δ：

其中x为通过某仪表的流量在一个时间段内的实际增量，y为同一时间段内该仪表读数的增量。

某封闭曲面s定义流量仪表集群as＝{mj|流j穿过封闭曲面s}，集群中流量仪表的个数ns＝n。设流量计mj第i次测量结果为yij，与之对应的实际值为xij。其中i＝1,2,…,m，m为测量次数。根据流量守恒约定有

若δj表示流量计mj的相对误差，则由定义节相对误差的定义可以得到将其代入(2)式，有成立。不失一般性，假定δ0已知，则有

令

代入(3)式有

(5)式可写成矩阵形式

yε＝-η(6)

其中

yi为仪表m1～mn-1第i次测量的结果向量。

ε＝(ε1ε2…εn-1)^t，η＝(y1,0ε0y2,0ε0…ym,0ε0)^t。

(6)式表明，若已知集群内某一仪表的相对误差δ0，则可以通过所有仪表的多次测量结果推算出其他各仪表的相对误差δj,j＝1,2,…,n-1。

2.添加虚拟支路计算损耗：

实际的电能表集群存在损耗(电能表损耗、漏电损耗、线路电阻损耗等)，我们引入虚拟支路的概念，代替问题中以线路损耗为主题的损耗值之和。该支路包括虚拟电能表和虚拟负载。集群总损耗可等效为虚拟负载的能耗。假想的封闭曲面s定义了一个电能表集群，电能表以符号mj(j＝0,1,…,n)表示。mn为虚拟电能表。约定流入s的能量为正，流出为负。设在第i个测量时段ti内流过mj的电能为xi,j(i＝1,2,…,m，为标识每个测量时段的序号；m为测量次数)，则根据能量守恒定律有下式成立。

通常认为第i个测量时段内，mj的读数增量yi,j与xi,j具有如下关系。

yi,j＝(1+δj)xi,j(10)

(10)式中，δj为mj的(平均)相对误差。若令

则(11)式可化为

不失一般性，假定已知m0的相对误差δ0，从而可知ε0。规定虚拟电能表mn的相对误差δn＝0，从而有εn＝1及xi,n＝yi,n成立。因此(2-22)式化为

3.虚拟支路中损耗的表达

(13)式中左边第三项xi,n表示测量时段ti内虚拟负载的能耗(集群总的损耗)，(17)式和(19)式方程组都为能量守恒方程组。它们表明在能量守恒的前提下，若已知集群内某一仪表的相对误差δ0，则有可能通过多次测量结果yi,j推算出其他所有电能表的εj，进而得到相对误差δj(j＝1,2,…,n-1)。在实际计算中，虚拟支路的流出为集群总表m0的读数增量减去各个分表mj(j＝1,…,n)的和，并将这个值乘以一个由经验值判定的损耗系数。该损耗系数作为电能表集群的能量损耗，依靠下文中所述循环迭代方法求得。对此，加入虚拟支路后，方程相应变化为α为修正经验，是一定范围内的随机值。

数据筛选：

该方法首先需要从用电信息采集系统中获取电能表相关抄表数据，该部分不再赘述。

该步骤通过填补遗漏数据、消除异常数据、平滑噪声数据，以及纠正不一致数据，去掉数据中的噪音、填充空值、丢失值和处理不一致数据。首先我们需要获取整个树形结构系统中所有表计的测量结果，在一起测量具有不确定性的情况下，我们可以尽量多的获取数据，以保证在合理的计算误差下结果的准确性和完备性要求。数据量最小应不小于电表数量加一。数据筛选有如下步骤：1.实时性检查。2.完整性检查和缺失值处理。3.可视化方法发现异常数据4.正交性检验

实时性检查：我们检查同一组表计抄表值是否在同一时间点抄读的。若不满足此要求，数据无法使用。

数据集成：数据集成是在原有抄表数据的基础上，将电能表抄表数据添加到该表计历史抄表数据中。

完整性检查和缺失值处理：首先我们需要对数据进行完整性检查。完整性检查的主要工作为缺失值处理。对于缺失值处理，传统的方法有直接删除，或利用平均值、中值、分位数、众数、随机值等替代。这样做效果一般，因为等于人对数据进行了修改。我们采用的方法通过历史误差和其他误差量做预测模型来算出缺失值，计算缺失值后添加到算法中进行二次验证，判断缺失值的情况是否可以满足缺失要求。具体的来说，若一组数据中缺失值较多，缺失量过大，我们会直接将该数据舍弃掉，因为如果引入会产生较大噪声，对结果造成不良影响。按照经验值判断，如果缺失值多于10％，我们更倾向于等待更多抄表数据。若能满足要求，我们首先将缺失采集数据电能表读数添加到虚拟支路的用电量中，使用原模型计算电能表误差情况。若数据添加至虚拟支路更新后计算结果保持在合理范围之内，我们再根据该台区平均线损经验占比，从虚拟支路的电量中扣除线损对数据进行填充。填充后，利用填充的表格重新计算误差情况，若满足两次其他表计变化不大且缺失表计与历史计算值偏差不大的条件，我们视为其结果填充正确。不失一般性，我们将缺失值电能表作为表1，在测量次数m不小于n-1的情况下，将(14)式写成矩阵形式

yε＝-η'(15)

其中yi为第i次测量结果向量；ε＝(ε2ε3…εn-1)^t；η'＝η-y·1＝(y1,0ε0-α(y1,0x1,n)-y11y2,0ε0-α(y2,0x2,n)-y12…yn-1,0ε0-α(y1,0xn-1,n)-y1n)^t。由上文所述，同理方程(15)可解。若ε＝(ε2ε3…εn-1)^t满足所述要求，那么方程合理，根据虚拟支路值恢复该结果。

数据分析(借助可视化工具)发现异常数据：将电能表每日抄表数据作为一个向量，那么多次抄读可以组成一个向量组。在多维数据分析技术中，通过将各个向量可视化表示，从中找出异常变动的特殊向量。使用crystalanalysis的可视化工具方法，可以发现特殊离群点值，对此点值做特殊分析，进一步增强数据监控，优化核心数据。另外，利用聚类算法将该向量聚类分析，作为在可视化表示发现异常数据基础上的一种特殊发现方法。

正交性检验：误差分析算法要求方程组各方程间具有弱相关性，否则方程组呈病态导致计算结果误差超出可接受范围。在满足足够的要求时方程才可计算。

计算方法：基于测量不确定度的电能表误差迭代计算方法

由于本方法测量的是电能表一段时间内的平均误差，测量的结果只是被测量的估计值，测量过程中的不确定因素会导致测量的不确定性。对于已知的系统效应进行修正后的测量结果仍然只是被测量的估计值，考虑不确定度评价方法，保证使在不完全理想的复现或样本代表性不够等的情况下，测量仍然能反映真实情况。

迭代法是一种不断用变量的旧值递推新值的过程，迭代算法是解决问题的一种基本方法。其基本思想是从变量的原值推出它的一个新值。在我们的方法中，迭代被应用于对虚拟电路修正经验α的估计。

在此首先给出假设：在一个台区内误差超差表的数量小于总表数量的5％，即为我们认为计量器具的可靠性在95％。

一般而言，计量器具的可靠性会随着时间的增加而降低，因此对于该比率在不同台区我们执行浮动的检定可靠性目标范围。根据台区内表计平均使用时间长短情况，我们对该值进行修正。

我们的方法简述如下：在给定的初始经验值α＝0的情况下，计算电能表集群误差情况，若该计算结果满足可信度要求，即计算的超差的电能表数量小于规定范围内，则该算法结束。否则向上调整α的值区间，返回算法求解。

联合不确定度评价方法：

数学模型:δ＝1/(y^-1η)-1

由于y,η满足独立性，则有：

那么我们有联合扩展不确定度u＝kuc取k＝2，有：

若求解的电能表误差符合实际预期，随后对被测量进行独立重复观测，通过所得到的一系列测得值，用统计分析方法获得实验标准偏差s(x)，当用算术平均值作为被测量估计值时，被测量估计值的不确定度按计算：

利用贝塞尔公式法，在重复性条件下对同一被测量独立重复观测n次，得到n个测得值xi＝(i＝1,2,...,n)，被测量x的最佳估计是n个独立测得值的算数平均值按如下公式计算：

单个测得值xk的实验方差s²(x)：

单个测得值xk的实验标准差s(x)：

被测量估计值的不确定度

尽管为说明目的公开了本发明的实施例和附图，但是本领域的技术人员可以理解：在不脱离本发明及所附权利要求的精神和范围内，各种替换、变化和修改都是可能的，因此，本发明的范围不局限于实施例和附图所公开的内容。