CN111506867B - 一种基于样条插值的变压器每相电量数据生成方法 - Google Patents

一种基于样条插值的变压器每相电量数据生成方法 Download PDF

Info

Publication number
CN111506867B
CN111506867B CN202010481458.9A CN202010481458A CN111506867B CN 111506867 B CN111506867 B CN 111506867B CN 202010481458 A CN202010481458 A CN 202010481458A CN 111506867 B CN111506867 B CN 111506867B
Authority
CN
China
Prior art keywords
phase
interpolation
function
transformer
power factor
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010481458.9A
Other languages
English (en)
Other versions
CN111506867A (zh
Inventor
彭青宁
邢伟
白杨
刘旭涛
李媛
马进
胡飞虎
吴小曼
王峰
张丽云
杨青
裴林
赵静
李一非
杨伟
陈宇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xixian New District Power Supply Company State Grid Shaanxi Electric Power Co
Xian Jiaotong University
Original Assignee
Xixian New District Power Supply Company State Grid Shaanxi Electric Power Co
Xian Jiaotong University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xixian New District Power Supply Company State Grid Shaanxi Electric Power Co, Xian Jiaotong University filed Critical Xixian New District Power Supply Company State Grid Shaanxi Electric Power Co
Priority to CN202010481458.9A priority Critical patent/CN111506867B/zh
Publication of CN111506867A publication Critical patent/CN111506867A/zh
Application granted granted Critical
Publication of CN111506867B publication Critical patent/CN111506867B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • Operations Research (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Complex Calculations (AREA)

Abstract

低压用户相位识别等应用需要用到变压器二次侧每相电量数据，然而现有系统并不存储该数据。直接使用变压器二次侧每相电压、电流和功率因数值相乘并求和可以得到变压器每相电量数据，但误差较大。为了解决现有方法获得的变压器每相电量数据误差较大的问题，本发明提出了一种基于样条插值的变压器每相电量数据生成方法，分别对采集的电压、电流和功率因数时间序列值进行三次样条插值，从而分别获得电压、电流和功率因数的插值函数；再将电压、电流、功率因数的插值函数相乘并按时间进行积分，就可得到更精确的变压器每相电量数据。

Description

一种基于样条插值的变压器每相电量数据生成方法
技术领域
本发明属于配用电技术领域，特别涉及一种基于样条插值的变压器每相电量数据生成方法。
背景技术
配电网的安全、可靠运行需要各类数据的支持。但在实际中，由于数据采集设备功能不足或系统存储容量有限等原因，配电变压器的部分运行数据并未被采集并存储，尤其是数量众多的低压配电变压器。由于部分运行数据未被采集，导致一些研究和应用受阻。例如，中国专利申请号201810498584.8公开的一种基于频谱分析的电力用户相位识别方法，在其权利要求1中，要求采集变压器关口表电表每一相位的电量数据。然而在实际中，变压器关口表电表并不采集和存储每一相位的电量数据。这也就意味着，必须通过其它方法获得变压器每相的电量数据。
由于变压器关口表电表每15分钟会采集二次侧每相的电压、电流和功率因数值，可以基于这些数据拟合生成每相的电量数据。例如，低压居民户每天只存储一笔电量数据，在识别用户挂载相位时，需要使用的变压器二次侧每相每天电量数据可以按照如下公式进行拟合：
Figure BDA0002517556680000011
其中，96为每天考核点的数量(每15分钟1次)；U、I、
Figure BDA0002517556680000012
分别为96个考核点的电压、电流、功率因数值；T为周期，值为0.25小时(15分钟)。然而，U、I、
Figure BDA0002517556680000013
并不是15分钟期间的平均有效值，而是瞬时有效值(即采集时间点的有效值)。由于采集数据的间隔时间较长，直接采用上述公式拟合生成的每相每天电量数据误差较大，严重影响了低压用户相位识别应用的准确率。
发明内容
为了克服上述现有技术的缺点，提高低压用户相位识别的准确率，本发明的目的在于提供一种基于样条插值的变压器每相电量数据生成方法，基于采集的变压器每相电压、电流和功率因数值，能够更精确地合成变压器每相电量数据。
为了实现上述目的，本发明采用的技术方案是：
一种基于样条插值的变压器每相电量数据生成方法，逐一对变压器的每个相位完成如下操作：
(1)采集该相位n个周期的电压值、电流值和功率因数值；
(2)分别对采集的n个周期的电压值、电流值和功率因数值进行样条插值，得到相应的插值函数；
(3)将电压值的插值函数、电流值的插值函数和功率因数值的插值函数相乘并对该n个周期的时间进行积分计算，获得该相位n个周期的电量数据。
与现有技术相比，本发明对电压、电流、功率因数进行三次样条插值，较为精确地填补了采样间隔时间内的缺失数据，得到了比原始采样数据信息更为丰富的电压、电流、功率因数的插值函数。通过对电压、电流、功率因数的插值函数相乘并按时间进行积分，由此可得到更精确的变压器每相电量数据，从而可以提高基于变压器每相电量数据的低压用户相位识别方法的准确率。
附图说明
图1是本发明流程示意图。
具体实施方式
下面结合附图和实施例详细说明本发明的实施方式。
低压用户相位识别等应用需要用到变压器二次侧每相电量数据，然而现有电表并不采集和存储该数据。直接使用电表采集的变压器二次侧每相电压、电流和功率因数值相乘并求和的方法可以得到变压器每相电量数据，但由于电压、电流和功率因数的采集间隔时间较长，该方法得到的结果误差较大。为了得到更精确的变压器每相电量数据，本发明提供一种基于样条插值的变压器每相电量数据生成方法。
如图1所示，一种基于样条插值的变压器每相电量数据生成方法，包括如下步骤：
(一)采集变压器所有相位n个周期的电压值、电流值和功率因数值；
(二)选择未被计算过的相位A。
(三)使用三次样条插值法计算选择相位电压在这n个周期内的插值函数，其具体步骤如下：
a.将n个周期的采样时刻构建为集合{xi}，其中，i＝0,1,2,……,n-1，x0＜x1＜…＜xn-1，n为采样周期数；xi时刻采集的采样值为yi，与{xi}对应的采样值的集合为{yi}。
b.用n-1段三次多项式在采样数据点之间构建一个三次样条，则插值函数u(x)在x∈[x0,xn-1]中为分段函数，如下公式1-1所示：
Figure BDA0002517556680000031
式中，u1(x)，…，un-1(x)即ui(x)，ui(x)为满足下列条件的三次多项式：
Ⅰ，插值条件：
Figure BDA0002517556680000032
Ⅱ，样条互相连接：
ui(xi)＝ui+1(xi)i＝1,2,…,n-2 (1-3)
其中，ui+1(xi)指将xi带入公式(1-1)右边的函数ui+1(x)中得到值，公式(1-3)表示公式(1-1)中的函数在定义域交集中的函数值相等，也即将i+1带入公式(1-2)第二个公式中时，公式(1-2)的第一个公式和第二个公式相等。
Ⅲ，在xi处，直到二阶导数均连续：
Figure BDA0002517556680000041
ui'(xi)是ui(xi)的一阶导数，ui”(xi)是ui(xi)的二阶导数；
选取边界条件为非扭结边界条件。
Figure BDA0002517556680000042
u1”'(x0)是函数u1(x)在x＝x0时的三阶导数，u2”'(x1)是函数u2(x)在x＝x1时的三阶导数，un-1”'(xn-1)是函数un-1(x)在x＝xn-1时的三阶导数，un-2”'(xn-2)是函数un-2(x)在x＝xn-2时的三阶导数。
c.根据(1-2)、(1-3)、(1-4)和边界条件(1-5)，共得到(4n-4)个方程，求解该方程组可以计算得到分段函数u(x)中共(4n-4)个参数，也即可以得到电压的插值函数u(x)，再令u(x)＝V(t)即可求得选择相位的电压的插值函数V(t)。
(四)与(三)同理，可以求解选择相位电流值的插值函数I(t)。
(五)与(三)同理，可以求解选择相位功率因数值的插值函数P(t)。
(六)计算该积分式
Figure BDA0002517556680000043
其中，T为n个周期的时间长度，t为时间，V(t)为电压值的插值函数，I(t)为电流值的插值函数，P(t)为功率因数值的插值函数。得到选择相位在n个周期的时间长度的用电量。
(七)分别选择其余相位，进行(三)、(四)、(五)和(六)，得到其余相位在在n个周期的时间长度的用电量。
(八)至此，已经求得该变压器三个相位在在n个周期的时间长度的每相用电量。
在本发明的一个实例中，以计算相位A的电压为例，n＝5，采集的电压如下表：
采集时刻 x0＝0 x1＝1 x2＝2 x3＝3 x4＝4
采集电压 y0＝219.0 y1＝221 y2＝220.5 y3＝223 y4＝220.4
此时插值函数为下式所示，共计16个未知数。
Figure BDA0002517556680000051
Ⅰ，插值条件，共计8个独立方程。
Figure BDA0002517556680000052
Ⅱ，样条互相连接，共计4个方程，但方程与I中方程相同。
Figure BDA0002517556680000053
Ⅲ，在xi处，直到二阶导数均连续：共计6个独立方程。
Figure BDA0002517556680000054
选取边界条件为非扭结边界条件。共计2个独立方程。
Figure BDA0002517556680000055
因此，总共得到16个独立方程，可以通过线性代数知识解得u(x)中16个未知数。令求得的解如下式所示
Figure BDA0002517556680000056
也即可以解得相位A的电压的插值函数为：
Figure BDA0002517556680000061
同理对相位A的电流和功率因数进行相同处理得到I(t)和P(t)
Figure BDA0002517556680000062
Figure BDA0002517556680000063
此时计算相位A的用电量时，
Figure BDA0002517556680000064
综上，本发明提出了一种基于样条插值的变压器每相电量数据生成方法，分别对采集的电压、电流和功率因数时间序列值进行三次样条插值，从而分别获得电压、电流和功率因数的插值函数；再将电压、电流、功率因数的插值函数相乘并按时间进行积分，就可得到更精确的变压器每相电量数据。

Claims (2)

1.一种基于样条插值的变压器每相电量数据生成方法，其特征在于，逐一对变压器的每个相位完成如下操作：
(1)采集该相位n个周期的电压值、电流值和功率因数值；
(2)分别对采集的n个周期的电压值、电流值和功率因数值进行样条插值，得到相应的插值函数；
(3)将电压值的插值函数、电流值的插值函数和功率因数值的插值函数相乘并对该n个周期的时间进行积分计算，获得该相位n个周期的电量数据；
其中，所述样条插值为三次样条插值，逐一对电压值、电流值和功率因数值完成如下操作：
a，将n个周期的采样时刻构建为集合{xi}，其中，i＝0,1,2,……,n-1，x0＜x1＜…＜xn-1，n为采样周期数；xi时刻采集的采样值为yi，yi为xi时刻对应的电压值或电流值或功率因数值；
b，构建如公式(1)所示的分段函数：
Figure FDA0004159922250000011
式中，u1(x)，…，un-1(x)即ui(x)，ui(x)为满足下列条件的三次多项式：
Ⅰ，插值条件：
Figure FDA0004159922250000012
Ⅱ，样条互相连接：
ui（xi）＝ui+1（xi）i＝1,2,…,n-2 (3)
Ⅲ，在xi处，直到二阶导数均连续：
Figure FDA0004159922250000021
ui'(xi)是ui(xi)的一阶导数，ui”(xi)是ui(xi)的二阶导数；
c，设置边界条件，根据公式(2)、(3)、(4)和边界条件计算得到插值函数u(x)；
所述边界条件为如公式(5)所示的非扭结边界条件：
Figure FDA0004159922250000022
u1”'(x0)是函数u1(x)在x＝x0时的三阶导数，u2”'(x1)是函数u2(x)在x＝x1时的三阶导数，un-1”'(xn-1)是函数un-1(x)在x＝xn-1时的三阶导数，un-2”'(xn-2)是函数un-2(x)在x＝xn-2时的三阶导数。
2.根据权利要求1所述基于样条插值的变压器每相电量数据生成方法，其特征在于，积分计算过程如下：
Figure FDA0004159922250000023
其中，T为n个周期的时间长度，t为时间，V(t)为电压值的插值函数，I(t)为电流值的插值函数，P(t)为功率因数值的插值函数。
CN202010481458.9A 2020-05-31 2020-05-31 一种基于样条插值的变压器每相电量数据生成方法 Active CN111506867B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010481458.9A CN111506867B (zh) 2020-05-31 2020-05-31 一种基于样条插值的变压器每相电量数据生成方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010481458.9A CN111506867B (zh) 2020-05-31 2020-05-31 一种基于样条插值的变压器每相电量数据生成方法

Publications (2)

Publication Number Publication Date
CN111506867A CN111506867A (zh) 2020-08-07
CN111506867B true CN111506867B (zh) 2023-05-09

Family

ID=71864442

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010481458.9A Active CN111506867B (zh) 2020-05-31 2020-05-31 一种基于样条插值的变压器每相电量数据生成方法

Country Status (1)

Country Link
CN (1) CN111506867B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114859284B (zh) * 2022-06-13 2025-01-28 云南电网有限责任公司电力科学研究院 电流互感器的阶跃响应测试方法和系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006254594A (ja) * 2005-03-10 2006-09-21 Tokyo Electric Power Co Inc:The 保護継電装置
CN103116698A (zh) * 2013-01-23 2013-05-22 合肥工业大学 一种基于三次样条的gm(1,1)模型预测方法
CN104700151A (zh) * 2014-05-26 2015-06-10 国网辽宁省电力有限公司 一种基于三次样条插值曲线拟合的风电功率考核方法
CN105610154A (zh) * 2016-01-18 2016-05-25 云南电网有限责任公司电力科学研究院 一种在线更换电能表电量补偿处理方法及装置
CN108132399A (zh) * 2017-12-15 2018-06-08 国家电网公司 一种提高数字化变电站电能质量分析精度的简化插值算法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10296030B2 (en) * 2015-10-07 2019-05-21 University Of Utah Research Foundation Systems and methods for power system management

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2006254594A (ja) * 2005-03-10 2006-09-21 Tokyo Electric Power Co Inc:The 保護継電装置
CN103116698A (zh) * 2013-01-23 2013-05-22 合肥工业大学 一种基于三次样条的gm(1,1)模型预测方法
CN104700151A (zh) * 2014-05-26 2015-06-10 国网辽宁省电力有限公司 一种基于三次样条插值曲线拟合的风电功率考核方法
CN105610154A (zh) * 2016-01-18 2016-05-25 云南电网有限责任公司电力科学研究院 一种在线更换电能表电量补偿处理方法及装置
CN108132399A (zh) * 2017-12-15 2018-06-08 国家电网公司 一种提高数字化变电站电能质量分析精度的简化插值算法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
宋晓林 ; 任伟 ; 张振铎 ; 阳星 ; 曾翔君 ; .数字化变电站合并单元插值误差对于电能计量的影响.电测与仪表.2017,(第11期),全文. *
张小开 ; 侯永亮 ; .失压情况下电能表的电能计算的新方法.山西电力.2011,(第02期),全文. *
朱家明 ; .两种三次样条插值建立方法及MATLAB运行时间比较.中国培训.2016,(第10期),全文. *
邓重一 ; 廖海生 ; .霍尔传感器特性曲线的样条函数拟合方法.微细加工技术.2008,(第02期),全文. *

Also Published As

Publication number Publication date
CN111506867A (zh) 2020-08-07

Similar Documents

Publication Publication Date Title
CN109799474B (zh) 一种基于数据相关性的多表位表箱电能表相位识别方法及系统
CN102033185A (zh) 一种基于视在功率估算的电能表错接线判断方法
CN104215882B (zh) 一种基于有源单端口网络电阻极性的电压暂降源定位方法
CN107069718A (zh) 基于三相不平衡对线路损耗影响的低压配网线损计算方法
CN108490361A (zh) 一种基于云端反馈的荷电状态SoC计算方法
CN111506867B (zh) 一种基于样条插值的变压器每相电量数据生成方法
CN104268798B (zh) 一种分布式电源及电动汽车对配电网影响的评估方法
CN110850333B (zh) 一种低压配电系统单相接地故障相别识别方法
CN111965479A (zh) 一种基于数据精细化量测比对的线损异常定位方法
CN108595376A (zh) 一种区分定频空调与冲击钻的非侵入负荷辨识方法
CN112180165B (zh) 基于智能表和终端运行数据低压供电线路阻抗测算方法
CN110865328B (zh) 基于ami的智能电表相位识别、拓扑辨识及阻抗估计方法
CN105048445B (zh) 计及多类型分布式电源的主动配电网三相状态估计方法
CN107621562A (zh) 一种用于分析窃电的单相表分流准确判断方法
CN115575884B (zh) 充电站内充电桩充电量对齐及计量误差计算方法
CN107767060B (zh) 配网线路理论线损计算系统及方法
CN114461982B (zh) 一种输电线路保护特性识别及电压暂降持续时间估计方法
CN112505613B (zh) 一种高供高计全失压计量故障损失电量估算方法及系统
CN115712811A (zh) 基于低压拓扑配电网下的阻抗计算模型
CN116579624A (zh) 一种基于改进等值电阻法的含光伏低压线损计算方法
CN112258007B (zh) 一种基于配变容量的配电网同期线损成本量化方法
CN112448389B (zh) 基于小扰动数据辨识负荷模型中电动机比例参数的方法与系统
CN105139091B (zh) 一种基于时间序列法的电容器电容值及其变化趋势预测方法
CN103543318B (zh) 一种周期性负荷电缆导体温度峰值时刻的确定方法
CN113985266A (zh) 多智能电表集中化负荷识别方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant