一种适用于配电网在线分析的改进直流潮流算法
【技术领域】
[0001] 本发明涉及电力系统潮流分析领域，具体为适用于配电网在线分析的改进直流潮 流算法
【背景技术】
[0002] 由于可在生能源的间歇性与随机性，其大量接入将对配电网造成广泛的影响，主 要表现在：改变配电网的电压水平、提高配电网的短路容量、继电保护策略的复杂度加大、 影响网络的供电可靠性以及加剧电能质量的恶化等。应对这种情况，2008年国际大电网会 议（CIGRE)配电与分布式发电专委会（C6)的C6. 11项目组提出：未来的配电网应能够使 用灵活的网络拓扑，对分布式能源（distributed energy resource, DER)进行主动管理与 控制，以实现可再生能源的高效友好接入。
[0003] 配网潮流计算是对配电网运行状态分析及分布式能源主动管理的基础，配电网 中分布式电源优化调度，在线协调控制等都需要用到配电网潮流数据。直流潮流算法由 于模型简单、方程组完全线性化、不存在收敛性问题，在电力系统静态安全分析（static Security Analysis)，实时安全调度（real-time security-constrained dispatch)等令页 域已取得了广泛的应用。但是，其模型忽略节点电压幅值偏差，仅能求取节点电压相角的近 似值，且计算精度受网络参数及运行状态影响大等劣势，限制了其进一步的应用。尤其是对 于配电网，线路具有较大的R/X值，采用直流潮流模型会产生的难以忽略的误差，限制了其 在配电网中的应用。
【发明内容】
[0004] 为克服现有标准直流潮流算法计算精度低，不能直接用于配电网潮流分析等问 题，本发明针对辐射状配电网提出一种可应用于在线潮流分析的改进直流潮流算法。其优 点如下：潮流方程为线性方程组，无需迭代求解，也不存在收敛性问题；对线路参数不作任 何假设，适用于任意参数的辐射状网络；不仅可以计算节点电压相角，也可以计算节点电压 幅值，可满足配电网在线分析的实时性与精度要求。
[0005] 本发明采用的技术方案是：
[0006] 1.获取辐射状配电网网络参数及节点功率注入信息，建立辐射状配电网节点关系 模型，获得辐射状配电网各支路有功、无功潮流与节点有功、无功注入的关系，并建立直流 潮流前推算法，具体如下：
[0007] 对于任意辐射状配电网，建立其节点关系：假设根节点为松弛节点，对于每一条支 路连接的两个节点，距离根节点远的节点为距离根节点近的节点的子节点。以图5所示9 节点配电系统为例，假设根节点1为松弛节点，节点2为节点1的子节点，节点3、7为节点 2的子节点，末端节点（8、9、6)没有子节点。
[0008] 对于任意辐射状配电网，忽略对地支路及线路有功、无功损耗，支路（i，j)的有功 及无功潮流可由下式计算：
[0009]
【主权项】
1. 一种适用于配电网在线分析的改进直流潮流算法，其特征在于： 1) 通过建立辐射状配电网节点关系模型，获得辐射状配电网各支路有功、无功潮流与 节点有功、无功注入的关系，并建立支路潮流前推算法获得支路潮流方程，求解各支路有 功、无功潮流； 2) 通过对支路潮流方程的线性化，获得辐射状配电网支路有功、无功潮流与支路连接 的两个节点电压相角的关系，并建立节点电压相角回代算法，求解支路节点电压相角； 3) 通过将支路潮流方程中节点电压幅值写为基值与偏移量之和的形式，并将已求得的 支路节点电压相角带入支路潮流方程，忽略电压幅值二阶偏移量值，得到以节点电压幅值 偏移量为未知数的线性方程组，求解此线性方程组即获得节点电压幅值的偏移量，即求得 所有节点的电压幅值。
2. 根据权利要求1所述的一种适用于配电网在线分析的改进直流潮流算法，其特征在 于步骤1)具体如下： 对于任意辐射状配电网，建立其节点关系：假设根节点为松弛节点，对于每一条支路连 接的两个节点，距离根节点远的节点为距离根节点近的节点的子节点；忽略对地支路及线 路有功、无功损耗，支路（i，j)的有功及无功潮流由下式计算：
式中Pij、Qij分别为从节点i流向节点j的有功、无功潮流；节点j为节点i的子节点；pjj、Qjj分别为节点j的有功、无功注入；kGj表示节点k为节点j的子节点。
3. 根据权利要求1所述的一种适用于配电网在线分析的改进直流潮流算法，其特征在 于，步骤2)具体如下： 对于支路（i，j)，节点电压幅值与相角与支路有功、无功潮流关系如下：
式中PU、QU分别为支路（i，j)上从节点i流向节点j的有功、无功潮流；V 分别为 节点i、j的电压幅值；RiPXu分别为支路（i，j)的串联电阻与电抗；0u为节点i与节点j的电压相角差； 对于配电网，通常有如下假设成立： 1) 正常运行的电网，各节点电压通常在额定电压附近，近似处理为Vj= 1 ; 2) 支路两端相角差很小，可近似认为：0 0,则sin0 0cos0 1; 则可得：
通过支路潮流前推算法，可计算整个配电网所有支路的有功潮流及无功潮流，在此基 础上可得，节点电压相角的计算公式： 0』=0「(XuPu-RuQj (4) 式中9i、9j分别为节点i与节点j的节点电压相角，节点j为节点i的子节点； 以配电网根节点为松弛节点，电压相角设为0,则可以从根节点向末端节点依次回代各 节点电压相角。
4.根据权利要求1所述的一种适用于配电网在线分析的改进直流潮流算法，其特征在 于步骤3)具体如下： 对一个N节点的配电网，其潮流方程为： S=diag[V] ? [Y]、[V]* (5) 式中，S为节点注入列向量；Y为节点导纳矩阵；V为节点电压相量列向量；diag[ ?]表 示以[?]中元素为对角线元素组成的对角阵；[?r中元素为[?]中元素的共轭； 将节点i的电压相量写为如下形式： V^iV^+AV^e'6- (6) 对于所有节点C=1;对于松弛节点与PV节点，AVi=Vsrt_l，Vsrt为其设定的电压幅 值；对于PQ节点，~ 1，故AV产0,为待求量； 那么，式（6)可展开如下： S=diag[Vb+AV] ?Y+ ? [Vb+AV] =diag[Vb] ?Y+ ? [Vb]+diag[AV] ?Y+ ? [Vb] (7) +diag[Vb] ?Y+? [AV]+diag[AV] ?Y+? [AV] 式中，Vb=[C,C,...,G]T，AV= [AVi，...，AVN]T; 0为N阶相角差方阵，其中的元素 0 ij= 0厂0 j; 0i、0」分别为节点i与节点j的节点电压相角，节点j为节点i的子节点； 设Y+为一个N阶方阵，其元素％ ，且节点导纳矩阵对应的元素，％为其共轭 值； 由假设忽略式⑵中的二阶扰动量（A\xAVj)，则式（5)化简为： A S = diag [ A V] ? Y+ ? [Vb] +diag [Vb] ? Y+ ? [ A V] (8) 式中，AS=S_diag[Vb] ?Y+ ? [Vb]; 又因为： diag[AV] ?Y+ ? [Vb] =diag[Y+ ?Vb] ? [AV] (9) 将式（9)带入式（8)并消去PV节点及松弛节点电压，可得： [AS]PQ=As.[AV]PQ (10) 其中，对于松弛节点与PV节点AV=vsrt-i，为已知量，[?]PQS[ ?]消去PV节点及 松弛节点数据后所得，4为（diag[Y+ ?Vb]+diag[Vb] ?Y+)消去PV节点及松弛节点数据所 得； 令：
其中，[AP]PA[AQ]PQ中元素分别为[AS] 元素的实部与虚部，AP与A#元素分 别为^中元素的实部与虚部；将式（11)带入式（10)得：
贝IJ由式（12)可以得出求解PQ节点电压幅值以下四种计算格式； (1) 复功率格式，以下简称S格式 由式（12)的复功率方程可得： [AV]pQ=As_1 ? [AS]pQ (13) 其中[?r1表示矩阵[?]的逆矩阵； 因节点电压幅值为实数，按（13)计算所得电压修正值为虚部很小的复数，忽略所得修 正量的虚部，仅将所得结果的实部作为节点电压幅值的修正量； (2) 有功功率格式，以下简称P格式 由式（12)中有功方程可得： [AV]pq= Ap^1 ? [AP]pq (14) (3) 无功功率格式，以下简称Q格式 由式（12)中有功方程可得： [AVjp^A^1* [AQ]pq (15) (4) 最小二乘格式，以下简称LSQ格式 由式（12)中的有功及无功方程可组成一个超定方程组，求得该超定方程组的最小二 乘解如下：
其中[*]T代表矩阵[?]的转置； 按以上四种格式计算得到PQ节点电压修正量后，即算得整个配网所有的节点的电压 幅值及相角信息。
【专利摘要】本发明公开一种适用于配电网在线分析的改进直流潮流算法，属于电力系统潮流分析领域，该方法能够应用于配电网在线分析，实时优化控制等领域。首先，通过从配电网末端至根节点的节点注入功率前推，求解配电网各个支路有功、无功潮流；其次，通过配电网根节点至配电末端的回代，求解配电网各个节点的电压相角；最后，忽略配电网节点电压幅值偏差二阶项，获得节点电压幅值与节点注入功率线性方程组，并求解该线性方程组，获得节点电压幅值。发明在保留标准直流潮流法优点的同时，对线路参数不作任何假设，适用于任意参数的辐射状电网，不仅可以计算节点电压幅值，也可计算节点电压相角，且求解精度与数值稳定性较标准直流潮流算法有大幅提高。
【IPC分类】G06F19-00
【公开号】CN104573384
【申请号】CN201510037070
【发明人】吴在军, 王洋, 窦晓波, 胡敏强, 黄冬冬, 周欣, 胡靖宜
【申请人】东南大学
【公开日】2015年4月29日
【申请日】2015年1月23日