CN104062501A - 双变压器变电站谐波状态估计方法 - Google Patents
双变压器变电站谐波状态估计方法 Download PDFInfo
- Publication number
- CN104062501A CN104062501A CN201410326249.1A CN201410326249A CN104062501A CN 104062501 A CN104062501 A CN 104062501A CN 201410326249 A CN201410326249 A CN 201410326249A CN 104062501 A CN104062501 A CN 104062501A
- Authority
- CN
- China
- Prior art keywords
- node
- phasor
- measurement
- voltage
- state
- 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
Abstract
本发明涉及一种双变压器变电站谐波状态估计方法,该方法将节点电压相量和注入电流相量中各元素分别按量测量和状态量分类,并按先量测量后状态量的顺序排列,然后将节点导纳矩阵相应分块,变换得到量测方程。对于两台三绕组变压器并行运行网络拓扑结构,已知电流量测点编号和电压量测点编号,将节点编号组合写出导纳矩阵的各分块矩阵,再代入量测方程中计算状态相量。该方法估计准确性高,易于实现且实现成本低。
Description
技术领域
本发明涉及电力系统谐波监测技术领域,特别涉及一种双变压器变电站谐波状态估计方法。
背景技术
在现有技术中,变电站谐波状态估计是利用同步相角测量单元(PMU)作为技术基础展开研究。广域测量系统中的同步相角测量单元(PMU)利用GPS系统的高精度授时信号,能实现对电力系统各个节点动态数据的同步采集。理论上,相角测量单元可以为大规模地区电网进行谐波状态估计提供技术支持。但是在实际电网中,只有某些局部电网安装有同步相量量测单元,而大部分地区电网已经安装的监测装置并没有GPS对时功能,而相角测量单元造价昂贵,如果重新安装同步相角测量单元,会大大增加谐波状态估计的成本。从应用范围看,现有技术方案适用于大区域电网,范围广,在谐波情况下,电网状态很容易受到参数分散性的影响。但是在现有技术方案的谐波状态估计中并没有考虑到参数分散性的因素。
发明内容
本发明的目的在于提供一种双变压器变电站谐波状态估计方法,该方法估计准确性高,易于实现且实现成本低。
为实现上述目的,本发明的技术方案是:一种双变压器变电站谐波状态估计方法,按如下方法进行变电站谐波状态估计:
对于两台三绕组变压器并行运行的情况,当有两侧绕组并列,一侧绕组分列运行时,将两台三绕组变压器的参数都归算到同一电压等级侧,并进行星-角变换,得到变压器等效电路:
节点1、节点2、节点3、节点4的电压相量分别为V 1、V 2、V 3、V 4,注入电流相量分别为I 1、I 2、I 3、I 4,节点1经等效阻抗Z 12与节点2连接,节点1经等效阻抗Z 13与节点3连接,节点1经等效阻抗Z 14与节点4连接,节点2经等效阻抗Z 23与节点3连接,节点2经等效阻抗Z 24与节点4连接;
选择节点1、3进行电压、电流全量测;
列出节点电压方程为:
按量测量和状态量节点编号对节点导纳矩阵Y分块,即
令,,,;
式中,Y kk 为节点导纳矩阵中对应谐波电流量测相量和谐波电压量测相量的子矩阵,Y kn 为节点导纳矩阵中对应谐波电流量测相量和谐波电压状态相量的子矩阵,Y nk 为节点导纳矩阵中对应谐波电流状态相量和谐波电压量测相量的子矩阵,Y nn 为节点导纳矩阵中对应谐波电流状态相量和谐波电压状态相量的子矩阵;
,满足系统可观条件;量测方程为:
;
对于两台三绕组变压器并行运行的情况,当有一侧绕组并列,两侧绕组分列运行时,将两台三绕组变压器的参数都归算到同一电压等级侧,并进行星-角变换,得到变压器等效电路:
节点1、节点2、节点3、节点4、节点5的电压相量分别为V 1、V 2、V 3、V 4、V 5,注入电流相量分别为I 1、I 2、I 3、I 4、I 5,节点1分别经等效阻抗Z 12、Z 13、Z 14、Z 15与节点2、节点3、节点4、节点5连接,节点2经等效阻抗Z 23与节点3连接,节点4经等效阻抗Z 45与节点5连接;
选择节点1、2进行电压、电流相量全量测,选择节点4进行电压相量量测,得到量测矩阵,状态矩阵,其中,V 1、V 2、V 4、I 1、I 2为量测装置测得的量测相量,V 3、V 5、I 3、I 4、I 5为待求的状态相量;
列出节点电压方程为:
按量测量和状态量节点编号对节点导纳矩阵Y分块,即
令,,,;
得到
子矩阵Y kn 和Y nk 都可逆,满足系统可观条件;量测方程为:
。
本发明的有益效果是:
1、采用变电站内已有的集成监测装置,运用多路信道同时监测多个量测点,较易实现同步,可行性更强,也大大降低了成本。
2、在局部电网进行状态估计,范围小,状态估计结果受参数分散性影响小,提高了状态估计的准确性。
附图说明
图1是本发明方法的实现流程图。
图2是本发明实施例中两台三绕组变压器并行运行,两侧绕组并列一侧绕组分列的电路图。
图3是图2的等效电路图。
图4是图3的等效电路图。
图5是本发明实施例中两台三绕组变压器并行运行,一侧绕组并列两侧绕组分列的电路图。
图6是图5的等效电路图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步的详细说明。
谐波状态估计理论:
谐波状态估计以节点电压作为状态量,节点电压、支路电流和节点注入电流作为量测量,用注入电流源表示非线性负荷,各次谐波之间实现了解耦。
谐波状态估计数学模型可表示为
Z = HX + ε (1)
式中,Z为m×1维量测相量;X为待定的n×1维状态相量;H为联系量测值与状态变量关系的m×n维量测矩阵,它与系统的拓扑结构、基本导纳矩阵和系统的量测配置方案都有关系;ε为由于随机噪声、变换器误差和通信问题而导致的m×1维量测误差相量;
对于h次谐波,以节点电压V N 作为状态量
1)如果量测量为节点注入电流I N ,则量测矩阵为关联节点电压V N 作和节点注入电流I N 的节点导纳矩阵Y NN ,量测方程为
I N (h) = Y NN (h)V N (h) (2)
2)如果量测量为节点电压量测V N ,则量测方程为
V N (h) = I* V N (h) (3)
I为单位矩阵;
3)如果量测量为支路电流IL,则量测矩阵为关联节点电压V N 和支路电流IL的支路导纳矩阵Y LN ,量测方程为
I L (h) = Y LN (h)V N (h) (4)
当量测矩阵H满秩时,系统完全可观。由于变电站网络结构简单,节点数目少,可以直接通过量测矩阵H满秩的条件来确定量测点数量和位置。
采用最小二乘法求最佳近似解,即有:
X = (H T H)-1 H T Z (5)
变电站谐波状态估计:
在变电站等效网络中,节点电压方程为:
I = YV (6)
式中,I为节点注入电流相量,V为节点电压相量,Y为节点导纳矩阵;
导纳矩阵Y中元素为:
(7)
对上式中电压相量V和注入电流相量I中各元素按量测量和状态量进行分类,并按先量测量后状态量的顺序排列,表示如下:
(8)
式中,I k 为三相谐波电流量测相量矩阵,I n 为三相谐波电流状态相量矩阵,V k 为三相谐波电压量测相量矩阵,V n 为三相谐波电压状态相量矩阵,Y kk 为节点导纳矩阵中对应谐波电流量测相量和谐波电压量测相量的子矩阵,Y kn 为节点导纳矩阵中对应谐波电流量测相量和谐波电压状态相量的子矩阵,Y nk 为节点导纳矩阵中对应谐波电流状态相量和谐波电压量测相量的子矩阵,Y nn 为节点导纳矩阵中对应谐波电流状态相量和谐波电压状态相量的子矩阵;
根据谐波状态估计理论,方程(5)可写为
(9)
当子矩阵Y niku 可逆时,得到量测方程:
(10)
量测矩阵为:
(11)
量测矩阵的行列式为:
(12)
当子矩阵Y kn 和Y nk 都可逆时,量测矩阵满秩,即系统状态完全可观测。
本发明双变压器变电站谐波状态估计方法,实现流程如图1所示。按如下方法进行变电站谐波状态估计:
对于两台三绕组变压器并行运行的情况,当有两侧绕组并列,一侧绕组分列运行时,电路图如图2所示,节点编号是任意的。将两台三绕组变压器的参数都归算到同一电压等级侧,等效电路图如图3所示。将图3中星形连接阻抗进行星-角变换,用三角形接法连接等效阻抗,得到如图4所示的变压器等效电路:
节点1、节点2、节点3、节点4的电压相量分别为V 1、V 2、V 3、V 4,注入电流相量分别为I 1、I 2、I 3、I 4,节点1经等效阻抗Z 12与节点2连接,节点1经等效阻抗Z 13与节点3连接,节点1经等效阻抗Z 14与节点4连接,节点2经等效阻抗Z 23与节点3连接,节点2经等效阻抗Z 24与节点4连接。
对于四端网络,可列出4个节点电压方程,需要4个量测量使网络可观。即对选两个节点进行电压电流量测。
下面对量测节点和未知节点编号进行排列组合,找到量测点的位置。由等效电路图4可以看出,节点1和节点2,节点3和节点4分别对称,因此分以下两种情况进行分析。电压量测节点和电流量测节点相同,根据节点编号直接写出
1)量测节点:1、3;未知节点:2、4
(13)
(14)
2)量测节点:1、2;未知节点:3、4
(15)
(16)
在工程中,两台并联的变压器一般阻抗电压参数相近,即Y 13≈Y 14,Y 23≈Y 24,因此,和的行列式即使不为0,也很小,会影响估计结果的精度。
根据以上分析,在两台变压器两侧并列,一侧分列情况下,在2个并联节点(节点1、节点2)和两个分列节点(节点3、节点4)各选一个进行电压、电流全量测:本实施例选择节点1、3进行电压、电流全量测;
列出节点电压方程为:
(17)
按量测量和状态量节点编号对节点导纳矩阵Y分块,即
(18)
按式(5)定义令
(19)
(20)
(21)
(22)
由式(17)(18)可知,,满足系统可观条件;量测方程为:
(23)
对于两台三绕组变压器并行运行的情况,当有一侧绕组并列,两侧绕组分列运行时,电路图如图5所示,同样,节点编号为任意的,与高中低压侧无关。将两台三绕组变压器的参数都归算到同一电压等级侧,并用三角形接法连接等效阻抗,得到如图6所示的变压器等效电路:
节点1、节点2、节点3、节点4、节点5的电压相量分别为V 1、V 2、V 3、V 4、V 5,注入电流相量分别为I 1、I 2、I 3、I 4、I 5,节点1分别经等效阻抗Z 12、Z 13、Z 14、Z 15与节点2、节点3、节点4、节点5连接,节点2经等效阻抗Z 23与节点3连接,节点4经等效阻抗Z 45与节点5连接;
该网络有5个节点,可列出5个节点电压方程,需要两个节点进行电压电流全量测,一个节点进行电压量测使网络可观。下面对量测节点和未知节点编号进行排列组合,找到量测点的位置。
由前面分析知,单台变压器网络可观需要一个节点电压电流量测,一个节点电压量测。因此,若使图6中网络可观,需要变压器T 1(节点1、2、3)和变压器T 2(节点1、4、5)各有一个进行电压电流全量测。因此有两种量测点配置方式:
1)并列母线(图中节点1)全量测,两台变压器的分列端(图中节点2、3和节点4、5)再各选一点分别进行全量测和电压量测
2)两台变压器的分列端(图中节点2、3和节点4、5)中各选一点全量测,剩余节点中任选一点进行电压量测。
在本实施例中,选择节点1、2进行电压、电流相量全量测,选择节点4进行电压相量量测,得到量测矩阵,状态矩阵,其中,V 1、V 2、V 4、I 1、I 2为量测装置测得的量测相量,V 3、V 5、I 3、I 4、I 5为待求的状态相量;
由节点电压方程I=YV,节点号按照先量测量后状态量的顺序排列
即
(24)
由导纳矩阵性质可知,Y ij = Y ji ,所以矩阵中元素全部用Y ij (i ≤ j)表示。
按量测量和状态量节点编号对节点导纳矩阵Y分块,即
(25)
按式(5)定义令 (26)
(27)
(28)
(29)
由式(4)中导纳矩阵元素的性质和式(24)(25)可得:
(30)
(31)
子矩阵Y kn 和Y nk 都可逆,满足系统可观条件;量测方程为:
。
为避免产生环流,两台变压器三端都并列的情况很少见。当三侧绕组都并列运行时,等效电路图与单台变压器等效电路图完全相同,等效阻抗参数通过并联计算得到。谐波状态估计算法与单台变压器相同。
以上是本发明的较佳实施例,凡依本发明技术方案所作的改变,所产生的功能作用未超出本发明技术方案的范围时,均属于本发明的保护范围。
Claims (1)
1.一种双变压器变电站谐波状态估计方法,其特征在于,按如下方法进行变电站谐波状态估计:
对于两台三绕组变压器并行运行的情况,当有两侧绕组并列,一侧绕组分列运行时,将两台三绕组变压器的参数都归算到同一电压等级侧,并进行星-角变换,得到变压器等效电路:
节点1、节点2、节点3、节点4的电压相量分别为V 1、V 2、V 3、V 4,注入电流相量分别为I 1、I 2、I 3、I 4,节点1经等效阻抗Z 12与节点2连接,节点1经等效阻抗Z 13与节点3连接,节点1经等效阻抗Z 14与节点4连接,节点2经等效阻抗Z 23与节点3连接,节点2经等效阻抗Z 24与节点4连接;
选择节点1、3进行电压、电流全量测;
列出节点电压方程为:
按量测量和状态量节点编号对节点导纳矩阵Y分块,即
令,,,;
式中,Y kk 为节点导纳矩阵中对应谐波电流量测相量和谐波电压量测相量的子矩阵,Y kn 为节点导纳矩阵中对应谐波电流量测相量和谐波电压状态相量的子矩阵,Y nk 为节点导纳矩阵中对应谐波电流状态相量和谐波电压量测相量的子矩阵,Y nn 为节点导纳矩阵中对应谐波电流状态相量和谐波电压状态相量的子矩阵;
,满足系统可观条件;量测方程为:
;
对于两台三绕组变压器并行运行的情况,当有一侧绕组并列,两侧绕组分列运行时,将两台三绕组变压器的参数都归算到同一电压等级侧,并进行星-角变换,得到变压器等效电路:
节点1、节点2、节点3、节点4、节点5的电压相量分别为V 1、V 2、V 3、V 4、V 5,注入电流相量分别为I 1、I 2、I 3、I 4、I 5,节点1分别经等效阻抗Z 12、Z 13、Z 14、Z 15与节点2、节点3、节点4、节点5连接,节点2经等效阻抗Z 23与节点3连接,节点4经等效阻抗Z 45与节点5连接;
选择节点1、2进行电压、电流相量全量测,选择节点4进行电压相量量测,得到量测矩阵,状态矩阵,其中,V 1、V 2、V 4、I 1、I 2为量测装置测得的量测相量,V 3、V 5、I 3、I 4、I 5为待求的状态相量;
列出节点电压方程为:
按量测量和状态量节点编号对节点导纳矩阵Y分块,即
令,,,;
得到
子矩阵Y kn 和Y nk 都可逆,满足系统可观条件;量测方程为:
。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410326249.1A CN104062501B (zh) | 2014-07-10 | 2014-07-10 | 双变压器变电站谐波状态估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410326249.1A CN104062501B (zh) | 2014-07-10 | 2014-07-10 | 双变压器变电站谐波状态估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104062501A true CN104062501A (zh) | 2014-09-24 |
CN104062501B CN104062501B (zh) | 2017-01-25 |
Family
ID=51550305
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410326249.1A Active CN104062501B (zh) | 2014-07-10 | 2014-07-10 | 双变压器变电站谐波状态估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104062501B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110824385A (zh) * | 2019-11-22 | 2020-02-21 | 国网福建省电力有限公司 | 基于低频激励响应的变压器系统电流相量测试系统及方法 |
CN111579849A (zh) * | 2020-04-10 | 2020-08-25 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种谐波电流分布获取方法及装置 |
CN112018760A (zh) * | 2020-08-27 | 2020-12-01 | 四川大学 | 一种基于非同步数据的谐波状态估计方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000004539A (ja) * | 1998-06-12 | 2000-01-07 | Toyo Electric Mfg Co Ltd | 電圧検出形アクティブフィルタの高調波電圧検出装置 |
CN102624000A (zh) * | 2012-04-12 | 2012-08-01 | 河海大学 | 一种基于自动微分的电力系统谐波状态估计方法 |
CN103018534A (zh) * | 2012-12-20 | 2013-04-03 | 广东电网公司电力科学研究院 | 确定谐波电压的方法及系统 |
CN103310121A (zh) * | 2013-07-10 | 2013-09-18 | 湖南大学 | 一种大型光伏电站与配网谐波交互影响分析模型建模方法 |
CN103323688A (zh) * | 2012-03-20 | 2013-09-25 | 华北电力科学研究院有限责任公司 | 一种电力系统的谐波状态估计方法及设备 |
-
2014
- 2014-07-10 CN CN201410326249.1A patent/CN104062501B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2000004539A (ja) * | 1998-06-12 | 2000-01-07 | Toyo Electric Mfg Co Ltd | 電圧検出形アクティブフィルタの高調波電圧検出装置 |
CN103323688A (zh) * | 2012-03-20 | 2013-09-25 | 华北电力科学研究院有限责任公司 | 一种电力系统的谐波状态估计方法及设备 |
CN102624000A (zh) * | 2012-04-12 | 2012-08-01 | 河海大学 | 一种基于自动微分的电力系统谐波状态估计方法 |
CN103018534A (zh) * | 2012-12-20 | 2013-04-03 | 广东电网公司电力科学研究院 | 确定谐波电压的方法及系统 |
CN103310121A (zh) * | 2013-07-10 | 2013-09-18 | 湖南大学 | 一种大型光伏电站与配网谐波交互影响分析模型建模方法 |
Non-Patent Citations (2)
Title |
---|
刘小江 等: "基于网络拓扑的谐波状态估计量测点最优配置", 《电力系统保护与控制》 * |
谭志勇 等: "考虑谐波相序特性的谐波状态估计算法的研究", 《电力系统保护与控制》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110824385A (zh) * | 2019-11-22 | 2020-02-21 | 国网福建省电力有限公司 | 基于低频激励响应的变压器系统电流相量测试系统及方法 |
CN110824385B (zh) * | 2019-11-22 | 2022-03-08 | 国网福建省电力有限公司 | 基于低频激励响应的变压器系统电流相量测试系统及方法 |
CN111579849A (zh) * | 2020-04-10 | 2020-08-25 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种谐波电流分布获取方法及装置 |
CN111579849B (zh) * | 2020-04-10 | 2022-09-16 | 中国南方电网有限责任公司超高压输电公司检修试验中心 | 一种谐波电流分布获取方法及装置 |
CN112018760A (zh) * | 2020-08-27 | 2020-12-01 | 四川大学 | 一种基于非同步数据的谐波状态估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104062501B (zh) | 2017-01-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Della Giustina et al. | Electrical distribution system state estimation: measurement issues and challenges | |
Wang et al. | An alternative method for power system dynamic state estimation based on unscented transform | |
Du et al. | On-line estimation of transmission line parameters, temperature and sag using PMU measurements | |
Zhao et al. | On-line PMU-based transmission line parameter identification | |
CN102135571B (zh) | 超高压/特高压多回路输电线路零序阻抗抗干扰测量方法 | |
Wang et al. | Transmission lines positive sequence parameters estimation and instrument transformers calibration based on PMU measurement error model | |
US20140244189A1 (en) | System For Achieving Real-Time Monitoring and State Estimation in Power Distribution Networks | |
CN106226653A (zh) | 多级电网的电压暂降的传递规律评估系统以及评估方法 | |
CN110927519B (zh) | 基于μPMU量测值的主动配电网故障定位方法 | |
CN202794435U (zh) | 无线高精度同步相位基准装置 | |
Chen et al. | State estimation for distribution systems using micro-synchrophasors | |
Zamora-Cárdenas et al. | State estimation of flexible AC transmission systems considering synchronized phasor measurements | |
CN104749488A (zh) | 一种基于连续数据窗的直流线路时域故障测距方法 | |
Li et al. | Fault location algorithm in distribution networks considering distributed capacitive current | |
CN104062501A (zh) | 双变压器变电站谐波状态估计方法 | |
Meliopoulos et al. | Distributed state estimator advances and demonstration | |
Ahmad Khan et al. | PTP‐based time synchronisation of smart meter data for state estimation in power distribution networks | |
Yang et al. | PMU placement for optimal three-phase state estimation performance | |
CN103972889B (zh) | 一种配电线路阻抗在线辨识方法 | |
CN104078975B (zh) | 单变压器变电站谐波状态估计方法 | |
CN109858061A (zh) | 面向电压功率灵敏度估计的配电网络等值与化简方法 | |
Goh et al. | Dynamic estimation of power system stability in different Kalman filter implementations | |
Kahunzire et al. | Improving distribution network state estimation by means of phasor measurement units | |
Gomathi et al. | Optimal location of PMUs for complete observability of power system network | |
CN104090161A (zh) | 基于集成监测装置的变电站谐波状态估计方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |