CN109191543A - 一种交流采样同断面数据生成方法 - Google Patents
一种交流采样同断面数据生成方法 Download PDFInfo
- Publication number
- CN109191543A CN109191543A CN201811093241.XA CN201811093241A CN109191543A CN 109191543 A CN109191543 A CN 109191543A CN 201811093241 A CN201811093241 A CN 201811093241A CN 109191543 A CN109191543 A CN 109191543A
- Authority
- CN
- China
- Prior art keywords
- data
- same
- acquisition device
- frequency
- power grid
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 32
- 238000005070 sampling Methods 0.000 title claims abstract description 27
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 6
- 230000005611 electricity Effects 0.000 claims description 28
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000005516 engineering process Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/203—Drawing of straight lines or curves
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/24—Classification techniques
Abstract
本发明提供一种交流采样同断面数据生成方法,涉及电力系统自动化技术领域。该方法首先使用不同的采集装置采集电网的实时电量数据,同一采集装置采集的实时电量数据按照不同的频率进行分类;再将同一采集装置采集的同一频率下的数据运用三参数正弦波曲线拟合算法进行拟合,得到不同频率下采集数据的分段拟合函数;然后确定一个时间基准点,并求出在时间轴上不同频率下同一时间点的时间,并将求出的时间带入所求得的分段拟合函数中,求同断面时间点的对应数据值;最后将求出的所有同断面的数据在数据展示平台形成同断面数据。本发明提供的交流采样同断面数据生成方法,生成的同断面数据误差小,精度高,能为高级应用软件分析提供更加准确的数据。
Description
技术领域
本发明涉及电力系统自动化技术领域,尤其涉及一种交流采样同断面数据生成方法。
背景技术
目前,随着现代科技水平的不断发展,电子科技行业例如雷达、通信、医疗、仪器等都达到了一个新的高度,对数据的精度要求越来越高,电网内部也越来越智能化,电网中存在很多高级应用软件,准确的数据会使高级应用软件做出合理的决策,但现在的电力测量的实数据为非同断面数据,这种数据误差较大,影响了各类高级应用软件的分析效果,当电网出现故障后,不利于做出正确的决策,不利于电网安全稳定运行,因此,需要提供电网同断面数据,为应用软件分析提供准确的数据。
最小二乘法在解决采集数据处理和误差方面得到广泛应用,并且可提高数据处理的效率和精确度,已成为数据处理的重要的、比较可靠的技术手段。根据两个变量x和y的n组实验数据,(xi,yi),i=1,2,…,n来求得这两个变量的函数关系的近似式(经验公式),这个过程就是曲线拟合。在选择了曲线的数学公式模型之后,使用最小二乘法选择曲线方程中的系数,一般通过对偏差的平方和对每个系数求偏导数,使偏导数为零来建立方程组,经过一系列推导,得到系数公式,然后求得系数。正弦波形的等间距采样序列获得其三参数拟合正弦曲线函数,是一种基本信号处理方法,在许多场合下获得了应用,如评价数据采集系统的有效位数、采集速率、交流增益、通道间延迟、触发特性等,在调制信号的数字化解调和失真度测量中也有应用。
针对三参数正弦波曲线拟合算法,已有众多学者作过大量研究工作,对于这类算法,其相对精度、绝对精度、效率、收敛性以及运行时间等基本可以满足工程应用的要求。
发明内容
针对现有技术的缺陷,本发明提供一种交流采样同断面数据生成方法,实现采集的电网数据的同断面数据的生成,为高级应用软件分析提供完整的数据。
一种交流采样同断面数据生成方法,包括以下步骤:
步骤1、使用电网中不同的采集装置采集电网的实时电量数据,并对同一采集装置采集的电量数据按不同的频率进行拟合;
所述对同一采集装置采集的电量数据按不同的频率进行拟合具体方法为:
步骤1.1、使用电网中的采集装置采集不同频率(f1,f2,f3,……fm)下的电网实时电量数据,并把所采集的数据按照不同的频率进行分类;得到同一采集装置采集的实时电量数据为:[(t1,1,y1,1)、(t1,2,y1,2)、……(t1,n,y1,n)]、[(t2,1,y2,1)、(t2,2,y2,2)、……(t2,n,y2,n)]、……[(tr,1,yr,1)……(tr,i,yr,i)……(tr,n,yr,n)]……[(tm-1,1,ym-1,1)、(tm-1,2,ym-1,2)、……(tm-1,n,ym-1,n)]、[(tm,1,ym,1)、(tm,2,ymm,2)、……(tm,n,ym,n)];
其中,tr,i为频率fr时,采集的第i个电量数据的时间;yr,i为频率fr时,采集的第i个电量数据,i=1、2、...、n,n为电网中同一采集装置采集的电网实时电量数据,r=1、2、...、m,m为不同的频率总数;
步骤1.2、将电网中同一采集装置采集的同一频率下的电网实时电量数据[(tr,1,yr,1)……(tr,i,yr,i)……(tr,n,yr,n)]运用三参数正弦波曲线拟合算法进行拟合,具体方法为:
步骤1.2.1、求同一采集装置在频率为fr时的采样周期;
已知采集装置的采样速率为vr,采样间隔为Δtr,采样个数为n,则每个采样周期的角频率为
步骤1.2.2、求同一采集装置在同一采样频率下采集的电量数据的拟合曲线;
同一采集装置在同一采样频率下采集的数据序列为(tr,1,yr,1)……(tr,i,yr,i)……(tr,n,yr,n),则采集的数据的正弦离散形式为:
yr,i=Brsin(ωrtr,i)+Crcos(ωrtr,i) (2)
其中,Br和Cr为均系数;
选取Br和Cr,使如下公式所示的残差平方和最小:
此时的系数Br=Br0和Cr0为同一频率下采集的电量数据的最小二乘拟合值;
步骤1.2.3、所述选取Br和Cr,使残差平方和最小的具体方法为:
构造如下矩阵:
则将残差平方和公式化简为矩阵形式,如下公式所示:
当最小时,得到未知θr的最小二乘解如下公式所示:
则同一频率下采集的电量数据的拟合函数如下公式所示:
根据拟合函数求出三角函数的最佳拟合曲线里的幅值Ar和相角αr,如下公式所示:
进而得到同一采集装置的第r个周期采集的电量数据的函数关系,如下公式所示:
g(tr)=Arsin(ωrtr+αr)
步骤1.3、重复步骤1.2,得到同一采集装置不同频率下采集数据的分段拟合函数,g(t1),g(t2),…,g(tm);
步骤2、确定电网所有待求时刻的同断面数据;
步骤2.1、确定一个时间基准点,并求出在时间轴上不同频率下同一时间点的时间;
设以频率f1取点的数据周期内t1为基准时间,则以频率fr周期内取点的对应时间点如下公式所示:
步骤2.2、把步骤2.1求出的时间带入由步骤1.3所求得的分段拟合函数的表达式中,求同断面时间点的对应数据值,如下公式所示:
g(tr)=Arsin(ωrtr+αr)
步骤2.3、重复步骤2.2,算出r=1,2,3,…,m时的g(tr),进而求出所有同断面的数据;
步骤3:将求出的所有同断面的数据在数据展示平台形成同断面数据。
由上述技术方案可知,本发明的有益效果在于:本发明提供的一种交流采样同断面数据生成方法,突破传统的的电力测量数据分析方法,利用已采集数据点,通过算法拟合,拟合最优的数学曲线,运用已求的数学曲线,确定在不同频率下同一时刻的数据,实现电网数据的同断面采集,为高级应用软件分析提供完整的数据。本发明方法生成的同断面数据误差小,精度高,能为高级应用软件分析提供更加准确的数据。
附图说明
图1为本发明实施例提供的一种交流采样同断面数据生成方法的流程图;
图2为本发明实施例提供的电网中同一采集装置采集的不同频率下的电网实时电量数据的示意图;
图3为本发明实施例提供的电网中同一采集装置在不同频率下采集的数据的拟合曲线示意图;
图4为本发明实施例提供的电网中同一采集装置产生的同断面数据示意图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
本实施例以某电网为例,使用本发明的交流采样同断面数据生成方法,将该电网的实时测量数据生成同断面数据。
一种交流采样同断面数据生成方法,如图1所示,包括以下步骤:
步骤1、使用电网中不同的采集装置采集电网的实时电量数据,并对同一采集装置采集的电量数据按不同的频率进行拟合;
所述对同一采集装置采集的电量数据按不同的频率进行拟合具体方法为:
步骤1.1、使用电网中的采集装置采集不同频率(f1,f2,f3,……fm)下的电网实时电量数据,并把所采集的数据按照不同的频率进行分类,如图2所示;得到同一采集装置采集的实时电量数据为:[(t1,1,y1,1)、(t1,2,y1,2)、……(t1,n,y1,n)]、[(t2,1,y2,1)、(t2,2,y2,2)、……(t2,n,y2,n)]、……[(tr,1,yr,1)……(tr,i,yr,i)……(tr,n,yr,n)]……[(tm-1,1,ym-1,1)、(tm-1,2,ym-1,2)、……(tm-1,n,ym-1,n)]、[(tm,1,ym,1)、(tm,2,ym,2)、……(tm,n,ym,n)];
其中,tr,i为频率fr时,采集的第i个电量数据的时间;yr,i为频率fr时,采集的第i个电量数据,i=1、2、...、n,n为电网中同一采集装置采集的电网实时电量数据,r=1、2、...、m,m为不同的频率总数;
步骤1.2、将电网中同一采集装置采集的同一频率下的电网实时电量数据[(tr,1,yr,1)……(tr,i,yr,i)……(tr,n,yr,n)]运用三参数正弦波曲线拟合算法进行拟合,具体方法为:
步骤1.2.1、求同一采集装置在频率为fr时的采样周期;
已知采集装置的采样速率为vr,采样间隔为Δtr,采样个数为n,则每个采样周期的角频率为
步骤1.2.2、求同一采集装置在同一采样频率下采集的电量数据的拟合曲线;
同一采集装置在同一采样频率下采集的数据序列为(tr,1,yr,1)……(tr,i,yr,i)……(tr,n,yr,n),则采集的数据的正弦离散形式为:
yr,i=Brsin(ωrtr,i)+Crcos(ωrtr,i) (2)
其中,Br和Cr为均系数;
选取Br和Cr,使如下公式所示的残差平方和最小:
此时的系数Br=Br0和Cr0为同一频率下采集的电量数据的最小二乘拟合值;
步骤1.2.3、所述选取Br和Cr,使残差平方和最小的具体方法为:
构造如下矩阵:
则将残差平方和公式化简为矩阵形式,如下公式所示:
当最小时,得到未知θr的最小二乘解如下公式所示:
则同一频率下采集的电量数据的拟合函数如下公式所示:
根据拟合函数求出三角函数的最佳拟合曲线里的幅值Ar和相角αr,如下公式所示:
进而得到同一采集装置的第r个周期采集的电量数据的函数关系,如下公式所示:
g(tr)=Arsin(ωrtr+αr)
步骤1.3、重复步骤1.2,得到同一采集装置不同频率下采集数据的分段拟合函数,g(t1),g(t2),…,g(tm),如图3所示;
步骤2、确定电网所有待求时刻的同断面数据;
步骤2.1、确定一个时间基准点,并求出在时间轴上不同频率下同一时间点的时间;
设以频率f1取点的数据周期内t1为基准时间,则以频率fr周期内取点的对应时间点如下公式所示:
步骤2.2、把步骤2.1求出的时间带入由步骤1.3所求得的分段拟合函数的表达式中,求同断面时间点的对应数据值,如下公式所示:
g(tr)=Arsin(ωrtr+αr)
步骤2.3、重复步骤2.2,算出r=1,2,3,…,m时的g(tr),进而求出所有同断面的数据;
步骤3:将求出的所有同断面的数据在数据展示平台形成如图4所示的同断面数据。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明权利要求所限定的范围。
Claims (4)
1.一种交流采样同断面数据生成方法,其特征在于:包括以下步骤:
步骤1、使用电网中不同的采集装置采集电网的实时电量数据,并对同一采集装置采集的电量数据按不同的频率进行拟合;
步骤2、求出在时间轴上不同频率下同一时间点的时间,确定电网所有待求时刻的同断面数据;
步骤3、将求出的所有同断面的数据在数据展示平台形成同断面数据。
2.根据权利要求1所述的一种交流采样同断面数据生成方法,其特征在于:步骤1所述对同一采集装置采集的电量数据按不同的频率进行拟合具体方法为:
步骤1.1、使用电网中的采集装置采集不同频率(f1,f2,f3,……fm)下的电网实时电量数据,并把所采集的数据按照不同的频率进行分类;得到同一采集装置采集的实时电量数据为:[(t1,1,y1,1)、(t1,2,y1,2)、……(t1,n,y1,n)]、[(t2,1,y2,1)、(t2,2,y2,2)、……(t2,n,y2,n)]、……[(tr,1,yr,1)……(tr,i,yr,i)……(tr,n,yr,n)]……[(tm-1,1,ym-1,1)、(tm-1,2,ym-1,2)、……(tm-1,n,ym-1,n)]、[(tm,1,ymm,1)、(tm,2,ymm,2)、……(tm,n,ym,n)];
其中,tr,i为频率fr时,采集的第i个电量数据的时间;yr,i为频率fr时,采集的第i个电量数据,i=1、2、...、n,n为电网中同一采集装置采集的电网实时电量数据,r=1、2、...、m,m为不同的频率总数;
步骤1.2、将电网中同一采集装置采集的同一频率下的电网实时电量数据[(tr,1,yr,1)……(tr,i,yr,i)……(tr,n,yr,n)]运用三参数正弦波曲线拟合算法进行拟合;
步骤1.3、重复步骤1.2,得到同一采集装置不同频率下采集数据的分段拟合函数,g(t1),g(t2),…,g(tm)。
3.根据权利要求2所述的一种交流采样同断面数据生成方法,其特征在于:所述骤步1.2的具体方法为:
步骤1.2.1、求同一采集装置在频率为fr时的采样周期;
已知采集装置的采样速率为vr,采样间隔为Δtr,采样个数为n,则每个采样周期的角频率为
步骤1.2.2、求同一采集装置在同一采样频率下采集的电量数据的拟合曲线;
同一采集装置在同一采样频率下采集的数据序列为(tr,1,yr,1)……(tr,i,yr,i)……(tr,n,yr,n),则采集的数据的正弦离散形式为:
yr,i=Brsin(ωr tr,i)+Crcos(ωr tr,i) (2)
其中,Br和Cr为均系数;
选取Br和Cr,使如下公式所示的残差平方和最小:
此时的系数Br=Br0和Cr0为同一频率下采集的电量数据的最小二乘拟合值;
步骤1.2.3、所述选取Br和Cr,使残差平方和最小的具体方法为:
构造如下矩阵:
则将残差平方和公式化简为矩阵形式,如下公式所示:
当最小时,得到未知θr的最小二乘解如下公式所示:
则同一频率下采集的电量数据的拟合函数如下公式所示:
根据拟合函数求出三角函数的最佳拟合曲线里的幅值Ar和相角αr,如下公式所示:
进而得到同一采集装置的第r个周期采集的电量数据的函数关系,如下公式所示:
g(tr)=Arsin(ωrtr+αr)。
4.根据权利要求3所述的一种交流采样同断面数据生成方法,其特征在于:所述骤步2的具体方法为:
步骤2.1、确定一个时间基准点,并求出在时间轴上不同频率下同一时间点的时间;
设以频率f1取点的数据周期内t1为基准时间,则以频率fr周期内取点的对应时间点如下公式所示:
步骤2.2、把步骤2.1求出的时间带入由步骤1.3所求得的分段拟合函数的表达式中,求同断面时间点的对应数据值,如下公式所示:
g(tr)=Arsin(ωrtr+αr)
步骤2.3、重复步骤2.2,算出r=1,2,3,…,m时的g(tr),进而求出所有同断面的数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811093241.XA CN109191543B (zh) | 2018-09-19 | 2018-09-19 | 一种交流采样同断面数据生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811093241.XA CN109191543B (zh) | 2018-09-19 | 2018-09-19 | 一种交流采样同断面数据生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109191543A true CN109191543A (zh) | 2019-01-11 |
CN109191543B CN109191543B (zh) | 2022-12-09 |
Family
ID=64908305
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811093241.XA Active CN109191543B (zh) | 2018-09-19 | 2018-09-19 | 一种交流采样同断面数据生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109191543B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111144705A (zh) * | 2019-12-05 | 2020-05-12 | 国网辽宁省电力有限公司锦州供电公司 | 一种基于带时标信息采集的全网同断面数据处理方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103837795A (zh) * | 2014-02-18 | 2014-06-04 | 国网山东省电力公司 | 基于广域故障录波信息的调度端电网故障诊断方法 |
WO2015196735A1 (zh) * | 2014-06-23 | 2015-12-30 | 华南理工大学 | 基于啮合频率和频谱校正技术的风电齿轮箱阶次跟踪方法 |
CN106204684A (zh) * | 2016-07-13 | 2016-12-07 | 国家海洋信息中心 | 一种基于任务流的海洋专题图件自动化编制方法 |
CN108258722A (zh) * | 2018-01-09 | 2018-07-06 | 国网辽宁省电力有限公司 | 一种提升电网新能源接纳能力的功频技术 |
-
2018
- 2018-09-19 CN CN201811093241.XA patent/CN109191543B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103837795A (zh) * | 2014-02-18 | 2014-06-04 | 国网山东省电力公司 | 基于广域故障录波信息的调度端电网故障诊断方法 |
WO2015196735A1 (zh) * | 2014-06-23 | 2015-12-30 | 华南理工大学 | 基于啮合频率和频谱校正技术的风电齿轮箱阶次跟踪方法 |
CN106204684A (zh) * | 2016-07-13 | 2016-12-07 | 国家海洋信息中心 | 一种基于任务流的海洋专题图件自动化编制方法 |
CN108258722A (zh) * | 2018-01-09 | 2018-07-06 | 国网辽宁省电力有限公司 | 一种提升电网新能源接纳能力的功频技术 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111144705A (zh) * | 2019-12-05 | 2020-05-12 | 国网辽宁省电力有限公司锦州供电公司 | 一种基于带时标信息采集的全网同断面数据处理方法 |
CN111144705B (zh) * | 2019-12-05 | 2023-04-18 | 国网辽宁省电力有限公司锦州供电公司 | 一种基于带时标信息采集的全网同断面数据处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109191543B (zh) | 2022-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106610450B (zh) | 一种单相电能计量芯片 | |
CN102664588A (zh) | 用于频率调制和幅度调制信号的数字化解调方法及装置 | |
CN101915874A (zh) | 一种基于傅立叶变换的谐波检测方法 | |
CN103941088A (zh) | 一种基于三相信号的电力系统频率快速测量方法 | |
CN104697436A (zh) | 一种基于傅里叶级数的圆感应同步器误差模型分析方法 | |
CN108845194B (zh) | 基于幅值波动量峰度检测原理的谐波责任量化方法 | |
CN109284933B (zh) | 一种基于数理统计的电子式互感器状态评估系统及方法 | |
CN104849545A (zh) | 一种电力系统频率测量方法及测量装置 | |
CN101718582A (zh) | 风力发电机组音调测试方法 | |
CN105699738A (zh) | 一种基于pwm的交流信号有效值测量方法 | |
CN109061345B (zh) | 适用于电力系统的有效值测量方法与系统 | |
CN104181391A (zh) | 数字功率计谐波检测的方法 | |
CN106771582A (zh) | 高频脉冲信号的测试方法及测试仪 | |
CN109191543A (zh) | 一种交流采样同断面数据生成方法 | |
CN106990286B (zh) | 基于欧拉方法的四象限谐波电能计量装置及方法 | |
CN104198811B (zh) | 低频信号频率测量方法和装置 | |
CN104407197B (zh) | 一种基于三角函数迭代的信号相量测量的方法 | |
CN112415457A (zh) | 一种电力综合测试仪交流量采集和校准实现方法及系统 | |
CN113177325A (zh) | 复杂地形下校正标准k-ε模型可调参数方法、装置及存储介质 | |
CN103605323B (zh) | 化工生产的离散控制方法及装置 | |
US7205916B1 (en) | Digital method and apparatus for resolving shaft position | |
CN116381333A (zh) | 基于acmd算法的电动汽车充电桩电能计量方法 | |
CN109858170B (zh) | 一种强非线性Modelica系统模型的频率特性估算系统及方法 | |
CN107085133A (zh) | 用于计算单相有效值的方法及装置 | |
CN109521269A (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 |