专利名称：非侵入式用电设备单元电流在线量测方法
技术领域：
本发明涉及一种用电设备单元电流在线量测方法。特别是涉及一种能够在不影响电 力用户正常生产、生活的情况下得到其内部主要用电设备单元电流参数的非侵入式用电 设备单元电流在线量测方法。
背景技术：
电力负荷分解数据对全社会(包括电力公司和电力用户)节能和电力系统规划及运
行都具有十分重要的意义。
用电设备单元电流参数是非侵入式电力负荷分解的数据基础，参数的获取目前主要
采用离线统计的方法。离线统计方法主要适合于实验室研究或具有较好独立测量条件的 电力用户。实验室研究就是有针对性地购买一些用电设备，对其进行大量的量测和统计； 但由于用电设备的种类繁多，实验室无法购买所有种类的用电设备，如果用户的测量条 件允许也可在用户内部对用电设备进行大量的量测和统计。
离线统计方法无疑是一种十分科学、准确获取用电设备单元电流参数的方法，但由 于电力负荷的复杂多样以及量测条件的限制，无法实现对每一类用电设备都进行单独大 量采样统计得到其单元电流。
发明内容
本发明所要解决的技术问题是，提供一种能够在不影响电力用户正常生产、生活的 情况下得到其内部主要用电设备单元电流近似参数的非侵入式用电设备单元电流在线量 测方法。
本发明所采用的技术方案是 一种非侵入式用电设备单元电流在线量测方法，用于 电力系统负荷能耗分解与监测、用电分析、负荷节能、负荷预测和负荷建模，包括如下 步骤
第一步在电力负荷入口处监测电压、电流信息，并判断是否有用电设备投入，如 果没有则继续监测，如果有则转入下一步；
第二步记录用电设备投入前、后电力负荷的稳态电压、电流数据； 第三步对上一步所记录的数据进行谐波分析得到电流表达式； 第四步计算用电设备投入前后电力负荷稳态电流的差值；
第五步对上一步所计算的电流的差值进行单元化，得到用电设备的单元电流参数。 第三步所得到的电流表达式是《(O = c。s(A欲+《,J
式中，a )是某一用电设备的单元电流，下标a表示用电设备；AeAT，即A:是一个 正整数；A^表示单元电流中第A:次谐波分量的角频率；《，。表示单元电流中第A次谐波 分量的初相角；^，。=1，当^:>1时，A,。是一个小于l的非负实数。
第四步所述的用电设备投入前后电力负荷稳态电流的差值是
=Ai,腳.[1. cos(W + ^腳)+…+ a汰腳.cos一 + U +…] —.[1.cos一 + 6>n) +... + aft .cos(Aw + ^汰)h—] =/al. [1 cos(纽+《J +…+ "。* . cos(A:w, +《J +…]
式中，U0、 、 &"、气_和^_分别表示用电设备投入使用后A、
A、 c^和^的新值；/。(0是计算得到的该用电设备稳态工作电流的瞬时值；^是用电 设备稳态工作电流基波分量的幅值；《,是用电设备稳态工作电流基波分量的初相角；aat 表示用电设备稳态工作电流第A:次谐波分量幅值与基波分量幅值的比例系数；《t表示用
电设备稳态工作电流中第A次谐波分量的初相角。
第五步所述的用电设备的单元电流参数是，根据稳态工作电流可得到这一用电设备 的单元电流，如下所示
= 1.cos一 + 6>al) +…+.cos(to + +
式中，ao是这一用电设备的单元电流；^是用电设备稳态工作电流基波分量的初
相角；"。t表示用电设备稳态工作电流第A次谐波分量幅值与基波分量幅值的比例系数； 《A表示用电设备稳态工作电流中第A:次谐波分量的初相角。
本发明的非侵入式用电设备单元电流在线量测方法，能够在不影响电力用户正常生 产、生活的情况下得到其内部主要用电设备单元电流的参数，其更有利于非侵入式电力 负荷分解的开展。
图1是本发明的非侵入式用电设备单元电流在线量测方法的流程图； 图2是空调启动全过程负荷实测电压、电流波形图； 图3是空调启动前负荷实测稳态电压、电流波形图； 图4是空调启动后负荷实测稳态电压、电流波形图。
具体实施例方式
下面结合实施例和附图对本发明的非侵入式用电设备单元电流在线量测方法做出详细说明。
如图1所示，本发明的非侵入式用电设备单元电流在线量测方法，是用于电力系统 负荷能耗分解与监测、用电分析、负荷节能、负荷预测和负荷建模，包括如下步骤
第一步在电力负荷入口处监测电压、电流信息，并判断是否有用电设备投入，如 果没有则继续监测，如果有则转入下一步；
第二步记录用电设备投入前、后电力负荷的稳态电压、电流数据；
第三步对上一步所记录的数据进行谐波分析得到电流表达式；
本发明进行非侵入式用电设备单元电流在线量测时，当供电电压符合国标规定时， 用电设备的单元电流，即其电流基波幅值的标幺值为l时的电流具有一定的统计规律，可 表示成如下形式
《")=Z aM cos(■ +《,。) (1)
式中，aO是某一用电设备的单元电流，下标"表示用电设备；A:eAT，即yt是一个 正整数；A:w表示单元电流中第A:次谐波分量的角频率；《，。表示单元电流中第A次谐波 分量的初相角；《M=1，当&>1时，A，。是一个小于l的非负实数。
用电设备投入使用前电力负荷的稳态电流可表示成如下形式
/,(0 =.[1.cos一 + + —!""汰.cos(A:欣+《》+] (2)
式中，Z,(O表示电力负荷稳态电流的瞬时值；下标/表示电力负荷；^表示电流负荷 稳态电流基波分量的幅值；o是电力负荷稳态电流基波分量的角频率；^是电力负荷稳 态电流基波分量的初相角；"w表示电力负荷稳态电流第)fe次谐波分量幅值与基波分量幅 值的比例系数；^表示电力负荷稳态电流第)t次谐波分量的初相角；y^y表示电力负荷稳
态电流第A次谐波分量的角频率。
第四步计算用电设备投入前后电力负荷稳态电流的差值； '当电力负荷内部有一主要用电设备投入使用后，其新的稳态电流如下-
《腳(,)=.[1.cos(纽+ 腳)+... + or汰，.cos(A:欲+ 6>汰拜)+…](3)
式中，(,腳(O、厶,，、A,卿、 腳和&,腳分别表示用电设备投入使用后i,(0、 /,i、 6>n、 c^和A的新值。
该用电设备的稳态工作电流可近似用其投入使用后的负荷稳态电流减去其投入使用
前的负荷稳态电流得到，如式(4)所示
= [1. eos(w, +《,腳)+…+ a汰，.cos(tor + 6>汰，)+…]
—.[1.cosO/ + (9") +... + ortt cos(A:c;, + (9议)+…] (4) =.[1.cos(M +《》+... + 0 。*.cos(A:纽
5式中，"O是计算得到的该用电设备稳态工作电流的瞬时值；d是用电设备稳态工 作电流基波分量的幅值；《t是用电设备稳态工作电流基波分量的初相角；a&表示用电 设备稳态工作电流第A:次谐波分量幅值与基波分量幅值的比例系数；《4表示用电设备稳
态工作电流中第A次谐波分量的初相角。
第五步对上一步所计算的电流的差值进行单元化，得到用电设备的单元电流参数。 根据其稳态工作电流可得到这一用电设备的单元电流，如下所示
= 1.cos(纽+ <9。!) +…+ ofa4.cos(to + 。 +... (5)
式中，ao是这一用电设备的单元电流；《,是用电设备稳态工作电流基波分量的初
相角；a。4表示用电设备稳态工作电流第A:次谐波分量幅值与基波分量幅值的比例系数； &表示用电设备稳态工作电流中第A:次谐波分量的初相角。
因此，当负荷内部有主要用电设备投入，且其它用电设备没有变化时，只要得到其 投入前、后的负荷稳态电流就可近似的计算出这一用电设备的单元电流。
图2是实测得到的空调启动全过程负荷电压、电流波形，这一过程可分为如图中所 示的三个部分
(1) 空调启动前稳态；
(2) 空调启动暂态；
(3) 空调启动后稳态； 」 图3是图2中第(1)部分的电压、电流实测波形。经谐波分析，其电压表达式为
、(0 = 225.09V^[1.0 cosOO + 0.015 cos(3w -167.33°) + 0.003 cos(5W — 121.13°)
+ 0.0006cos(7W +149.41°) + 0.0004cos(9w +10.86°) + 0扁6cos(1 krf + 2.12°)]
电流表达式式为
z;0) = 6.8V^[1.0cos(W + 3.11°) + 0.18cos(3w — 24.74°) + 0.067 cos(5W — 53.06°) + 0.04cos(7w - 2.82°) + 0.032cos(9w + 2.08°) + 0.013cos(1 l纽+131.25°)] 图4是图2中第(3)部分的电压、电流实测波形，其电压表达式为 v3 0) = 221.06V^[1.0 cos(W) + 0.016 cos(3 W — 174.91°) + 0.0025 cos(5w — 119.93°)
+ 0.0006cos(7W -141.22°) + 0.0005cos(9w +10.79°) + 0.0006cos(1 + 4.55°)]
电流表达式式为
/30) = 17.75V^[1.0cos(W —1.22°) + 0.09cos(3w — 55.73°) + 0.037 cos(5W — 118.13°) + 0.0095 cos(7W + 0.3°) + 0.0095 cos(9w + 33.74°) + 0.0052 cos(l kyf - 85.37°)]
空调稳态电流可用其投入前后的负荷稳态电流差值来得到
=10.99V^[1.0cos(纽—3.91°) + 0.076cos(3纽—105.6°) + 0.057cos(5w -159.31°) + 0.0096cos(7纽-7.84°) + 0.01cos(9纽—48.44°) + 0.0064cos(1 l纽—23.35°)]其中，L(O表示空调稳态电流。 那么，空调的单元电流为
= l.Ocos(纽-3.91°) + 0.076cos(3c^ —105.6°) + 0.057 cos(5M —159.31°) + 0.0096cos(7纽—7.84°) + 0.01cos(9w, — 48.44°) + 0.0064 cos(l l纽—23.35°)
cxo表示空调的单元电流，即其基波标幺值为i的电流。
根据求得的空调单元电流，对图2中的第(3)部分进行负荷分解，结果如表1所示:
表l
类别真实比例分解结果分解偏差pQPQPQ
电阻类19. 2%0%23%0%3. 8%0%
空调类61. 8%231. 7%58. 7%223. 9%-3. 1%-7. 8%
计算机8.2%25. 4%8.2%25 . 60/00%0. 2%
日光灯10. 8%-157. 1%10. 1%-149. 5%-0. 7%7. 6%
从分解结果看，这种非侵入式用电设备单元电流参数获取方法可较好的为电力负荷 能耗分解提供数据基础，促进负荷分解的开展。
权利要求
1.一种非侵入式用电设备单元电流在线量测方法，其特征在于，用于电力系统负荷能耗分解与监测、用电分析、负荷节能、负荷预测和负荷建模，包括如下步骤第一步在电力负荷入口处监测电压、电流信息，并判断是否有用电设备投入，如果没有则继续监测，如果有则转入下一步；第二步记录用电设备投入前、后电力负荷的稳态电压、电流数据；第三步对上一步所记录的数据进行谐波分析得到电流表达式；第四步计算用电设备投入前后电力负荷稳态电流的差值；第五步对上一步所计算的电流的差值进行单元化，得到用电设备的单元电流参数。
2. 根据权利要求l所述的非侵入式用电设备单元电流在线量测方法，其特征在于， 第三步所得到的电流表达式是(o = S "*,。
cos(^+《,。)式中，aO是某一用电设备的单元电流，下标"表示用电设备；A:eiV'，即A:是一个 正整数；A:w表示单元电流中第A:次谐波分量的角频率；《，。表示单元电流中第A次谐波 分量的初相角；气。=1，当A〉1时，A，。是一个小于1的非负实数。
3. 根据权利要求l所述的非侵入式用电设备单元电流在线量测方法，其特征在于， 第四步所述的用电设备投入前后电力负荷稳态电流的差值是=厶,腳.[1 cos(纽+ ^腳)+…+ or汰腳 cos(A:欲+ 6>议，)+…] -/;1 [1 cos(crf + +…+ a汰 cos(A:w, + &) +] =. [1. cos(纽+ cos(A:纽+《A) +…]式中'O /,一 、 A,腳、气腳和^，分别表示用电设备投人使用后/,(0、 vA、 c^和^的新值；/。")是计算得到的该用电设备稳态工作电流的瞬时值；/。,是用电 设备稳态工作电流基波分量的幅值；^是用电设备稳态工作电流基波分量的初相角；《。4 表示用电设备稳态工作电流第A次谐波分量幅值与基波分量幅值的比例系数；表示用电设备稳态工作电流中第A次谐波分量的初相角。
4. 根据权利要求l所述的非侵入式用电设备单元电流在线量测方法，其特征在于， 第五步所述的用电设备的单元电流参数是，根据稳态工作电流可得到这一用电设备的单 元电流，如下所示《(,)=1. cos(erf + 。 + . + ", cos(to/ + 。 + .. 式中，《(O是这一用电设备的单元电流；^是用电设备稳态工作电流基波分量的初相角；"ai表示用电设备稳态工作电流第A次谐波分量幅值与基波分量幅值的比例系数； &表示用电设备稳态工作电流中第A:次谐波分量的初相角。
全文摘要
本发明公开一种非侵入式用电设备单元电流在线量测方法，用于电力系统负荷能耗分解与监测、用电分析、负荷节能、负荷预测和负荷建模，有如下步骤在电力负荷入口处监测电压、电流信息，并判断是否有用电设备投入，如果没有则继续监测，如果有则转入下一步；记录用电设备投入前、后电力负荷的稳态电压、电流数据；对上一步所记录的数据进行谐波分析得到电流表达式；计算用电设备投入前后电力负荷稳态电流的差值；对上一步所计算的电流的差值进行单元化，得到用电设备的单元电流参数。本发明能够在不影响电力用户正常生产、生活的情况下得到其内部主要用电设备单元电流的参数，其更有利于非侵入式电力负荷分解的开展。
文档编号G01R19/00GK101576580SQ20091006913
公开日2009年11月11日 申请日期2009年6月4日 优先权日2009年6月4日
发明者余贻鑫, 葛少云, 黄纯华, 鹏 黎 申请人:天津天大求实电力新技术股份有限公司