2. 重庆大学输配电装备及系统安全与新技术国家重点实验室 重庆 沙坪坝区 400044;
3. 重庆工程职业技术学院电气工程学院 重庆 沙坪坝区 400037;
4. 国家电网长沙供电公司 长沙 410000
2. State Key Laboratory of Power Transmission Equipment & System Security and New Technology, Chongqing University Shapingba Chongqing 400044;
3. Electrical Engineering College, Chongqing Vocational Institute of Engineering Shapingba Chongqing 400037;
4. State Grid Changsha Electric Power Supply Company Changsha 410000
互联电网中对内网进行研究分析时,有必要考虑与之紧密联系的外网对内网的影响[1]。由于电力市场的竞争关系以及电网分层分区的运行管理特点,内网一般只能获取外网的典型运行方式数据,如大方式和小方式运行数据,而难以获取外网的实时同步运行数据[2, 3],给全网一体化分析带来困难。对此,非拓扑静态等值方法在外网实时同步运行数据未知的情况下计算外网等值参数,进行一体化分析以保留外网对内外的影响,如挂等值机法[4]和戴维南等值法[5, 6, 7],基于Ward等值电路或REI等值电路的方法[3, 8, 9, 10]。上述方法可分为单端口和多端口两类,本文研究对象为前者。这类方法,特别是戴维南等值法,已广泛应用于内网稳态分析[5, 6, 7]中。在这些应用中,不精确或者不恰当的非拓扑静态等值方法在电力系统分析中会产生较大的误差甚至错误的计算结果[11, 12],一个精确有效的非拓扑静态等值方法应尽可能地考虑如下两个因素。
1) 外网元件特性的全面性。
外网对内网电压和功率的支撑作用主要通过线路、变压器、发电机、负荷以及无功补偿装置等元件实现,而这些元件有着截然不同的特性。因此,需要将这些元件的不同特性通过等值参数保留在外网等值电路中,而现有单端口等值电路均没有全面考虑外网元件的特性[5, 6, 7, 11, 12],影响了其等值的有效性和精确性。
2) 外网等值参数的物理约束。
外网等值参数的合理性和准确性是单端口非拓扑静态等值方法得到有效应用的关键[7, 10]。为解决上述问题,文献[10]提出考虑了电阻、电抗的一般物理约束,电阻和电抗均大于0,电抗大于电阻;文献[7]提出电阻、电抗的上下限约束,但未给出其上下限取值方法。以上方法还存在约束范围宽泛,上下限取值困难的问题,所以仍有等值参数不合理和不准确的情况。因此,在单端口非拓扑静态等值方法中需进一步研究更为有效的等值参数物理约束。
针对上述问题本文做出如下研究:首先通过合理的等值参数将外网元件的特性保留在外网等值电路中,提出基于外网元件特性的单端口等值电路;基于该等值电路、边界节点处的多时段实时状态信息和外网的典型运行方式数据,建立了等值参数的量测方程和物理约束,从而提出带约束的静态等值优化方法;最后,利用新英格兰39节点系统分析本文等值参数的合理性和准确性。
1 考虑外网元件特性的单端口等值电路设计在单端口互联电网中,外网通过单个边界节点与内网紧密相联。由于外网通常包含了发电机、线路、变压器、负荷以及无功补偿装置等元件,外网对内网的电压和功率支撑作用均通过这些元件实现,而这些元件有着截然不同的特性。因此,需要将这些元件的不同特性尽可能全面地保留在单端口等值电路中,以充分反映外网对内网的影响。为此,本文提出考虑外网元件特性的单端口外网等值电路,如图 1所示。
在图 1中,${\dot U^t}$和${\dot I^t}$分别代表第t个时段边界节点处的实时状态信息:边界节点电压相量和注入内网的等效电流相量。外网等值电路包含了如下4类等值元件:等值电压源、等值阻抗支路、等值注入电流、等值对地支路,它们是考虑外网元件特性全面性的具体体现,其作用如下:
1) 等值电压源是对外网发电机的等值,其所在母线电压用等值参数$\dot E$表示。等值电压源在潮流计算中作为PV节点,被用于模拟外网发电机对边界节点及内网的电压和功率支撑作用。
2) 等值阻抗支路,其等值参数为Z,是对外网线路/变压器的等值。在实际电力系统中,线路/变压器的主要作用是功率的输送和分配,其特性可通过阻抗参数Z保留在外网等值电路中。
3) 等值注入电流,用等值参数$\dot I$表示,是对外网负荷在边界节点b处的等值。负荷是外网中的重要元件,将其特性保留在外网等值电路中,有助于提高静态等值精度。
4) 等值对地支路电纳,其等值参数为B,表征原外部网络中各种对地支路在边界节点b处的等值,包括线路/变压器模型中的对地支路、并联电容器和并联电抗器等。
综上,在外网等值电路中有必要将外网中各主要元件的特性通过以上等值元件保留在等值电路中,以全面反映外网对内网的电压和功率支撑。
2 考虑等值参数物理约束的静态等值优化方法本文基于等值电路、边界处多时段实时状态信息和外网典型运行方式数据,分别建立待求等值参数的量测方程和物理约束,并进一步提出带约束的静态等值优化模型及其求解流程。
2.1 等值参数量测方程基于边界处实时状态信息较易获得边界节点电压和内网等效注入电流,如图 1中所示${\dot U^t}$、${\dot I^t}$。在边界节点b处应用KCL,有:
$\frac{{\dot E - {{\dot U}^t}}}{Z} = {\dot I^t} + \dot I + {\text{j}}B{\dot U^t}$ | (1) |
式中,j是虚数的符号。
将式(1)左右两边同乘Z,移项并将各复数按实部和虚部展开,等式两端的实部和虚部分别相等,可以得到如下量测方程:
$\begin{gathered} f_1^t(x) = {E_{\operatorname{Re} }} - {U_{\operatorname{Re} }} - R(I_{\operatorname{Re} }^t + {I_{\operatorname{Re} }} - BU_{\operatorname{Im} }^t) + \hfill \\ \;\;\;\;\;\;\;\;\;\;\;\;\;\;\;X(I_{\operatorname{Im} }^t + {I_{\operatorname{Im} }} + BU_{\operatorname{Re} }^t) = 0 \hfill \\ \end{gathered} $ | (2) |
$\begin{gathered} f_2^t({\mathbf{x}}) = {E_{\operatorname{Im} }} - {U_{\operatorname{Im} }} - R(I_{\operatorname{Im} }^t + {I_{\operatorname{Im} }} + BU_{\operatorname{Re} }^t) - \hfill \\ \;\;\;\;\;\;\;\;\;\;\;\;\;X(I_{\operatorname{Re} }^t + {I_{\operatorname{Re} }} - BU_{\operatorname{Im} }^t) = 0 \hfill \\ \end{gathered} $ | (3) |
式中,x=[R,X,ERe,EIm,IRe,IIm,B]T,下标“Re”和“Im”分别代表了各变量的实部和虚部。
2.2 等值参数物理约束本节利用外网典型运行方式数据,借助优秀的拓扑法[13]计算外网等值参数的典型值,并提出等值参数物理约束的建立方法。
在静态等值的拓扑法领域,获得符合实际物理意义的外网等值参数,可通过外网的详细信息,如拓扑结构、电力元件参数及运行状态等数据计算得到。而在非拓扑法领域,虽然内网无法获得外网的实时同步数据,但可获得外网典型运行方式下的拓扑结构、电力元件参数及运行状态。因此,可基于外网典型运行数据,结合本文等值电路的特点,采用现有拓扑法计算得到物理意义明确的外网等值参数的典型值,且这些典型值能够较为准确地反映外网元件在典型运行方式下的运行状态,进而利用该等值参数典型值建立合理、有效的等值参数上下限约束。
现利用外网的大方式和小方式典型运行数据,采用文献[13]中方法计算本文外网等值电路中的等值阻抗Z、等值负荷电流幅值I和等值对地电纳B的典型值,即获得大、小方式下的两组外网等值参数向量为:
${{\mathbf{x'}}_{bom}} = [{R_{bom}},{X_{bom}},{I_{bom}},{B_{bom}}]$ | (4) |
${{\mathbf{x'}}_{som}} = [{R_{som}},{X_{som}},{I_{som}},{B_{som}}]$ | (5) |
式中,下标bom代表大方式,som代表小方式。
根据继电保护的知识,大方式下外网的负荷水平较高,发电机、线路、变压器等设备一般大量投运,此时外网的负荷节点与发电机节点之间的等值阻抗最小,外网负荷电流幅值最大,对地电纳最大。小方式下外网的负荷水平较低,外网中部分设备退出运行,此时外网的负荷节点与发电机节点之间的等值阻抗最大,外网负荷电流幅值最小,对地电纳最小。根据上述关系可建立外网等值参数R、X、IRe、IIm、B的物理约束为:
$\left\{ {\begin{array}{*{20}{c}} {{R_{bom}} \leqslant R \leqslant {R_{som}}} \\ {{X_{bom}} \leqslant X \leqslant {X_{som}}} \\ {{I_{som}} \leqslant \sqrt {{I_{\operatorname{Re} }}^2 + {I_{\operatorname{Im} }}^2} \leqslant {I_{bom}}} \\ {{B_{som}} \leqslant B \leqslant {B_{bom}}} \end{array}} \right.$ | (6) |
由于等值参数的物理约束来源于等值参数的典型值,但考虑到外网的当前运行方式较外网的典型运行方式有较大变化的情况,当前运行方式下的外网等值参数可能不在物理约束的范围内。此时可在外网等值阻抗、等值负荷电流幅值的典型值上扩大或缩小一定比例,以保证所设置的物理约束能够包含当前外网等值参数的实际解。在式(6)中各不等式的两边乘以柔性系数向量α和β,其中各个量可分别在0.8~1.0和1.0~1.2中取值,柔性系数的具体取值可在等值参数计算中进行调整。另外,根据发电机机端电压的调压要求对等值电压源电压幅值E的取值限定在[Emin,Emax]。
结合上述分析及等值参数的量测方程和物理约束,即式(2)、式(3)和式(6),可建立等值参数物理约束的单端口非拓扑静态等值优化模型为:
$\left\{ \begin{gathered} J = \min \left( {\sum\limits_{t = 1}^m {\sum\limits_{i = 1}^2 {{{(f_i^t({\mathbf{x}}))}^2}} } } \right) \hfill \\ {\text{s}}{\text{.t}}{\text{. }}\;{\alpha _1}{R_{bom}} \leqslant R \leqslant {\beta _1}{R_{som}} \hfill \\ \;\;\;\;\;{\alpha _2}{X_{bom}} \leqslant X \leqslant {\beta _2}{X_{som}} \hfill \\ \;\;\;\;\;{\alpha _3}{I_{som}} \leqslant \sqrt {{I_{\operatorname{Re} }}^2 + {I_{\operatorname{Im} }}^2} \leqslant {\beta _3}{I_{bom}} \hfill \\ \;\;\;\;\;{\alpha _4}{B_{som}} \leqslant B \leqslant {\beta _4}{B_{bom}} \hfill \\ \;\;\;\;\;{E_{min}} \leqslant \sqrt {{E_{\operatorname{Re} }}^2 + {E_{\operatorname{Im} }}^2} \leqslant {E_{max}} \hfill \\ \end{gathered} \right.$ | (7) |
式中,αj∈α,βj∈β,j=1,2,3,4。
待求变量x包含3个等值元件参数(R,X,B)和4个等值元件状态(ERe,EIm,IRe,IIm)。由式(2)和式(3)可知,独立量测方程数始终小于待求变量数。本文采用多时段的实时状态信息求解该模型,假设在多个时段间隔内外网等值元件参数近似不变的要求,而等值元件状态可能变化,但在等值精度允许范围内,仍然可将其作为常数进行估计。综合考虑等值参数的可估性和等值精度,取m=4,并运用内点法对上述模型求解。等值参数求解的具体流程如图 2所示,其中变量初值可根据外网典型参数进行选取,如取上下限的平均值。
为验证本文方法的有效性,选择本文方法、无约束方法、戴维南等值法[5]、挂等值机法[4]进行对比分析。无约束方法是在本文方法的基础上去掉所有约束,其他环节保持不变;戴维南等值法中考虑了电阻、电抗的一般物理约束。现利用新英格兰39节点系统对上述方法进行等值参数的全面性和合理性分析以及静态安全分析。
以新英格兰39节点系统中的节点16为边界,节点19、20、21、22、23、24、33、34、35、36为外网,其余节点为内网形成单端口的互联电网算例。假定全网负荷以及发电机有功出力按照0.15%的步长递增,然后分别进行全网潮流计算获得多个时段的状态信息。并通过增减发电机出力、增减负荷、投切线路和无功补偿装置等措施形成该系统的大方式和小方式典型运行数据。
3.2 评价指标定义相对误差er和安全误差es这两个能够反映外网等值精度的指标分别为:
${e_r} = \left| {\frac{{x - {x^{{\text{eq}}}}}}{x}} \right| \times 100\% $ | (8) |
${e_{\text{s}}} = \left| {\frac{{x - {x^{{\text{eq}}}}}}{{{S_{{\text{base}}}}}}} \right| \times 100\% $ | (9) |
式中,x和xeq分别代表真值和估计值,真值是在仿真条件下,内外均采用详细模型的全网潮流计算值,而估计值是在内网采用详细模型,外网采用等值模型下的全网潮流计算值;Sbase表示系统基准功率,本文算例中Sbase=305 MVA。
3.3 数据分析1) 外网元件特性全面性和等值参数合理性分析。
现利用本文方法以及上述方法对本文算例的外网进行等值参数计算,如表 1所示。从表中可以看出本文方法通过4类等值元件对外网元件进行等值且等值参数最为合理;无约束方法计算的等值电阻小于0,且等值负荷为发出有功,不具备合理性。此外,戴维南等值法和挂等值机法缺乏对外网元件特性的全面反映。
2) 静态安全分析。
在本文算例中进行线路或发电机开断,分析等值前后内网潮流的相对误差和安全误差。在开断内网支路15-16的情况下,内网潮流的相对误差及其安全误差如图 4a~图 4e所示。可以看出,本文方法的误差明显小于其他两种方法,如图 4c中挂等值机法为273.23%,戴维南等值法为121.99%,本文方法仅为2.59%;图 4e中挂等值机法为25.52%,戴维南等值法为11.39%,而本文方法仅为0.24%。
表 2列出了在内网中不同支路或发电机开断的情况下各方法的评价指标最大值。由表 2可见本文方法的相对误差和安全误差最小,其等值精度均要高于其他两种方法。根据3.2节对相对误差的定义可知本文方法的估计值非常接近真值,证明了本文方法的正确性,而其他两种方法的相对误差较大,偏离了真值。如当开断支路16-17时,各方法支路无功潮流的相对误差最大值分别为:挂等值机法289.21%,戴维南等值法106.91%,本文方法仅2.08%。
综上分析,本文方法能较为全面地反映外网主要元件的特性,所计算的等值参数更加合理有效,能够在内网运行方式发生较大变化后为其提供适当的电压和功率支撑,等值精度明显提高。
4 结论本文方法能够较为全面地反映外网主要元件的特性,进而准确地模拟实际外网对内网的电压和功率支撑作用,提高了静态安全分析的精度。本文方法仅适用于单端口的非拓扑静态等值,因此本文后续的研究工作是根据本文的研究思路对多端口的非拓扑静态等值方法展开深入研究,以提高非拓扑静态等值方法的工程实用性。
[1] | 余娟, 张勉, 朱柳, 等. 考虑灵敏度一致性的外网静态等 值新理论研究[J]. 中国电机工程学报, 2013, 33(10): 64-70. YU Juan, ZHANG Mian, ZHU Liu, et al. New theory on external network static equivalent based on sensitivity consistency[J]. Proceedings of the CSEE, 2013, 33(10): 64-70. |
[2] | POON K K Y, EMAMI R, BOSE A, et al. External data exchange issues for state estimation in power systems[J]. IEEE Transactions on Power Systems, 2012, 27(2): 849-856. |
[3] | 余娟, 朱柳, 颜伟, 等. 基于内网实测信息的多端口外网 静态等值方法[J]. 电力自动化设备, 2013, 33(6): 32-37. YU Juan, ZHU Liu, YAN Wei, et al. Multi-port external network static equivalent on the basis of internal network measured information[J]. Electric Power Automation Equipment, 2013, 33(6): 32-37. |
[4] | 张海波, 张伯明, 王俏文, 等. 不同外网等值模型对EMS 应用效果影响的试验研究[J]. 电网技术, 2006, 30(3): 1-6. ZHANG Hai-bo, ZHANG Bo-ming, WANG Qiao-wen, et al. A test investigationon effeet of different externai network equivalent models on EMS applications[J]. Power System Technology, 2006, 30(3): 1-6. |
[5] | VU K, BEGOVIC M M, NOVOSEL D, et al. Use of local measurements to estimate voltage-stability margin[J]. IEEE Transactions on Power Systems, 1999, 14(3): 1029-1035. |
[6] | 刘明松, 张伯明, 姚良忠, 等. 基于PMU和改进戴维南等 值模型的电压稳定在线监视[J]. 电力系统自动化, 2009, 33(10): 6-10. LIU Ming-song, ZHANG Bo-ming, YAO Liang-zhong, et al. Online voltage stability monitoring based on PMU and improved Thevenin equivalent model[J]. Automation of Electric Power Systems, 2009, 33(10): 6-10. |
[7] | 刘道伟, 韩学山, 任玲玉, 等. 基于轨迹灵敏度的戴维南 等效参数迭代优化辨识[J]. 中国电机工程学报, 2010, 30(S1): 37-42. LIU Dao-wei, HAN Xue-shan, REN Ling-yu, et al. Identification of Thevenin equivalent parameters using iterative optimization approach based on the trajectory sensitivity[J]. Proceedings of the CSEE, 2010, 30(S1): 37-42. |
[8] | 颜伟, 李世明, 陈俊, 等. 基于内网实测信息的两端口外 网静态等值参数估计方法[J]. 中国电机工程学报, 2011, 31(13): 101-106. YAN Wei, LI Shi-ming, CHEN Jun, et al. Parameter estimation for two-port external network static equivalent model on the basis of internal network measured information[J]. Proceedings of the CSEE, 2011, 31(13): 101-106. |
[9] | 李智, 杨洪耕. REI等值和变距线性回归在地区电网电压 无功优化控制中的应用[J]. 电力系统保护与控制, 2012, 40(1): 18-22. LI Zhi, YANG Hong-geng. Application of REI equivalent and varying-intercept linear regression in voltage and reactive power optimization control in regional power grid[J]. Power System Protection and Control, 2012, 40(1): 18-22. |
[10] | 郭文鑫, 李世明, 朱柳, 等. 考虑量测误差计及等值参 数物理约束的外网静态等值方法[J]. 电网技术, 2012, 36(10): 216-221. GUO Wen-xin, LI Shi-ming, ZHU Liu, et al. A method for static equivalence of external network considering measurement error and physical constraints of equivalent parameters[J]. Power System Technology, 2012, 36(10): 216-221. |
[11] | 李卫星, 牟晓明, 李志民. 电力系统戴维南等值参数的 解析与思考[J]. 中国电机工程学报, 2012, 32(S1): 28-34. LI Wei-xing, MOU Xiao-ming, LI Zhi-min. Analysis and thoughts on the thevenin equivalent parameters for power systems[J]. Proceedings of the CSEE, 2012, 32(S1): 28-34. |
[12] | WANG Yun-fei, PORDANJANI I R, LI Wei-xing, et al. Voltage stability monitoring based on the concept of coupled single-port circuit[J]. IEEE Transactions on Power Systems, 2011, 26(4): 2154-2163. |
[13] | YU J, ZHANG M, LI W Y. Static equivalent method based on component particularity representation and sensitivity consistency[J]. IEEE Transactions on Power Systems, 2014, 29(5): 2400-2408. |