WO2020154849A1 - 获取交流电力网中源荷的支路圴方电流分量的对称方法 - Google Patents

获取交流电力网中源荷的支路圴方电流分量的对称方法 Download PDF

Info

Publication number
WO2020154849A1
WO2020154849A1 PCT/CN2019/073442 CN2019073442W WO2020154849A1 WO 2020154849 A1 WO2020154849 A1 WO 2020154849A1 CN 2019073442 W CN2019073442 W CN 2019073442W WO 2020154849 A1 WO2020154849 A1 WO 2020154849A1
Authority
WO
WIPO (PCT)
Prior art keywords
node
power
source
branch
load
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.)
Ceased
Application number
PCT/CN2019/073442
Other languages
English (en)
French (fr)
Inventor
江辉
彭建春
王贵斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shenzhen University
Original Assignee
Shenzhen University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Shenzhen University filed Critical Shenzhen University
Priority to US16/630,304 priority Critical patent/US20210384728A1/en
Priority to PCT/CN2019/073442 priority patent/WO2020154849A1/zh
Priority to CN201980002743.3A priority patent/CN111758195B/zh
Publication of WO2020154849A1 publication Critical patent/WO2020154849A1/zh
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for AC mains or AC distribution networks
    • H02J3/28—Arrangements for balancing of the load in networks by storage of energy
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for AC mains or AC distribution networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for AC mains or AC distribution networks
    • H02J3/003—Load forecast, e.g. methods or systems for forecasting future load demand
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2103/00—Details of circuit arrangements for mains or AC distribution networks
    • H02J2103/30—Simulating, planning, modelling, reliability check or computer assisted design [CAD] of electric power networks
    • H02J2103/35—Grid-level management of power transmission or distribution systems, e.g. load flow analysis or active network management
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2105/00—Networks for supplying or distributing electric power characterised by their spatial reach or by the load
    • H02J2105/50—Networks for supplying or distributing electric power characterised by their spatial reach or by the load for selectively controlling the operation of the loads
    • H02J2105/52—Networks for supplying or distributing electric power characterised by their spatial reach or by the load for selectively controlling the operation of the loads for limitation of the power consumption in the networks or in one section of the networks, e.g. load shedding or peak shaving
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—ELECTRIC POWER NETWORKS; CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2105/00—Networks for supplying or distributing electric power characterised by their spatial reach or by the load
    • H02J2105/61—Load identification

Definitions

  • the present invention relates to the field of electric power engineering, and in particular to a symmetric method and a computer-readable storage medium for obtaining branch current components of (electric) source (negative) loads in an AC power network.
  • the branch mean square current (the square of the effective value) and the source load power is the key to efficiently realize the safety correction of the AC power grid and ensure its safe operation.
  • the branch mean square current component of the source load is a new tool for deep-level realization of AC power network safety correction, which is in urgent need of research and development.
  • the existing AC power network security correction methods either first construct a security correction optimization model, and then obtain a correction scheme through optimization; or first obtain the sensitivity of the branch current to the source load power, and then implement the approximate linear relationship expressed by the sensitivity. Because the former is an optimized model, it is not only unable to guarantee a reliable safety correction scheme, and the calculation amount of the optimization solution is always large. The latter is inaccurate due to the local linear characteristic of the sensitivity of the branch current to the source load power, and then Lead to repeated corrections.
  • the embodiments of the present invention provide a symmetric method and computer-readable storage medium for obtaining the branch current components of the source load in the AC power network, and aim to solve the problem of unreliability and low efficiency of the existing AC power network security correction method.
  • the first aspect of the embodiments of the present invention provides a symmetric method for obtaining the branch current components of the source load in an AC power network, including:
  • the M-P inverse matrix is used to establish a linear symmetric matrix expression of the node translation voltage and the node voltage phase of the entire network with respect to the source load power of the nodes of the entire network;
  • a symmetric algebraic formula for obtaining the branch mean square current component of the source load is established according to the symmetric algebraic expression and the Shapley value theorem of the branch mean square current with respect to the source load power of the nodes in the entire network.
  • the second aspect of the embodiments of the present invention provides a computer-readable storage medium, the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, the above-mentioned branch circuit of the source load in the AC power network is obtained. Steps of the symmetric method of square current components.
  • the above symmetric method for obtaining the branch current component of the source load in the AC power network is implemented according to the symmetric algebraic expression of the branch mean square current with respect to the source load power of the whole network node and the Shapley value theorem to establish the branch to obtain the source load.
  • the symmetrical algebraic calculation formula of the mean square current component realizes the acquisition of the branch mean square current component of the source load in the AC power network.
  • the symmetrical algebraic calculation formula for the branch mean square current component of the source load is applicable to the mean square current of all branches in the AC power network, the source load power of all nodes, and the source load power of all nodes are treated equally.
  • the branch mean square current component of is symmetrical to all source loads; on the other hand, the symmetrical algebraic calculation formula of the branch mean square current component of the source load is the full variable (not incremental) expression of the source load power Therefore, the large-scale variation of the source load power is accurate, and the amount of calculation is reduced.
  • This symmetrical and precise relationship between the branch mean square current component and the source-load power solves the problems of time-consuming, unreliable, inaccurate, and low-efficiency existing AC power network security correction methods.
  • FIG. 1 is an implementation flowchart of a method for obtaining a symmetrical current component of a branch circuit of a source load in an AC power network according to an embodiment of the present invention
  • Fig. 2 is a schematic structural diagram of a general model of an AC power network provided by an embodiment of the present invention.
  • an embodiment of the present invention provides a symmetric method for obtaining branch current components of source loads in an AC power network, including:
  • Step S101 establishing a linear expression of the node source load power with respect to the node translation voltage and the node voltage phase according to the known node source load power and branch admittance in the AC power network;
  • Step S102 establishing a steady-state linear symmetric model of the AC power network according to the linear expression of the node source-load power with respect to the node translation voltage and the node voltage phase;
  • Step S103 according to the steady-state linear symmetric model of the AC power network, use the M-P inverse matrix to establish a linear symmetric matrix expression of the node translation voltage and the node voltage phase of the entire network with respect to the node source load power of the entire network;
  • Step S104 establishing a symmetric algebraic expression of the branch mean square current with respect to the source load power of the whole network node according to the linear symmetric matrix expression of the translational voltage of the whole network node and the node voltage phase with respect to the source load power of the whole network node;
  • Step S105 Establish a symmetric algebraic formula for obtaining the branch mean square current component of the source load according to the symmetric algebraic expression of the branch mean square current with respect to the source load power of the whole network node and the Shapley value theorem.
  • the mean square current of all branches and the source load power of all nodes in the AC power network are calculated according to the above symmetric algebraic calculation formula, and the branch mean square current component of all node source loads can be obtained, thereby realizing the branch of the source load in the AC power network The acquisition of the mean square current component.
  • This symmetrical and accurate relationship between the branch mean square current component and the source-load power solves the problems of time-consuming, unreliable, inaccurate, and low-efficiency methods for the existing AC power network security correction symmetrical method.
  • step S101 the method of establishing a linear expression of the node source load power with respect to the node translation voltage and the node voltage phase according to the node source load power and branch admittance in the AC power network is specifically:
  • i and k are the number of nodes in the network AC power, and both belong to the set of consecutive natural numbers ⁇ 1,2, ..., n ⁇ ; n is the total number of said AC power network nodes; P i and Q i are source is connected to node i and the active charge-charge reactive power source, and charge the power source referred to as node i; P i is equal to the power supply active power to the node i minus load active power, said Q i It is equal to the reactive power of the power supply connected to the node i minus the reactive power of the load; g ik and b ik are the conductance and susceptance of the branch ik connected between the node i and the node k, and are collectively called the branch ik admittance; [theta] i and [theta] k are the phase voltage of the node i and the node k; V i and V k are offset in the voltage node i and node k
  • step S102 the method for establishing a steady-state linear symmetric model of the AC power network according to the linear expression of the node source-load power with respect to the node translation voltage and the node voltage phase is specifically:
  • i and j are the numbers of nodes in the AC power network, and both belong to the set of continuous natural numbers ⁇ 1,2,...,n ⁇ ; n is the total number of nodes in the AC power network; P 1 and Q 1 respectively Is the source load active power and source load reactive power connected to node 1, and collectively referred to as the source load power of node 1.
  • the P 1 is equal to the active power of the power source connected to node 1 minus the active power of the load, the Q 1 power supply to the node is equal to the reactive power by subtracting a load reactive power;
  • P i and Q i are connected to a source node i and the active charge-charge reactive power source, and charge the power source referred to as node i;
  • P i is equal to the active power between the power of node i by subtracting the load active power, Q i is equal to the power supply to the node i by subtracting the reactive power load reactive power;
  • P n and Q n are connected to a node
  • the source charge active power and source charge reactive power of n are collectively referred to as the source charge power of node n;
  • the P n is equal to the active power of the power supply connected to node n minus the active power of the load, and the Q n is equal to the active power connected to the node
  • the above-mentioned steady-state model of the power grid is linear, and all the node source load powers are included in the model and are treated equally. This is why it is called the linear symmetric model.
  • step S103 the method of using the M-P inverse matrix to establish the linear symmetric matrix expression of the node translation voltage and the node voltage phase of the entire network with respect to the source load power of the entire network node according to the linear symmetric model of the AC power network steady state is specifically:
  • i is the number of the node in the AC power network, and belongs to the set of continuous natural numbers ⁇ 1,2,...,n ⁇ ; n is the total number of nodes in the AC power network; ⁇ 1 , ⁇ i and ⁇ n are respectively The voltage phases of node 1, node i and node n; v 1 , v i and v n are the translational voltages of node 1, node i and node n respectively, and they are all standard unit voltages after translation -1.0; P 1 and Q 1 is the source load active power and source load reactive power connected to node 1, and collectively referred to as the source load power of node 1.
  • the P 1 is equal to the active power of the power supply connected to node 1 minus the active power of the load, so said Q is equal to 1 to the node 1 is the reactive power by subtracting power load reactive power;
  • P i and Q i are connected to a source node i and a source charge active reactive power source charge, and referred to as node i charge power;
  • P i is equal to the active power between the power of node i by subtracting the load active power, Q i is equal to the reactive power supply to the node i minus the load reactive power;
  • P n and Q n are The source load active power and source load reactive power connected to node n are collectively referred to as the source load power of node n;
  • the P n is equal to the active power of the power supply connected to node n minus the active power of the load, and the Q n is equal to The reactive power of the power supply connected to the node n minus the reactive power of the load;
  • step S104 the symmetric algebraic expression of the branch mean square current with respect to the source load power of the whole network node is established according to the linear symmetric matrix expression of the translational voltage of the whole network node and the node voltage phase with respect to the source load power of the whole network node
  • the specific method is:
  • i, k, and h are the numbers of nodes in the AC power network, and they all belong to the set of continuous natural numbers ⁇ 1,2,...,n ⁇ ; n is the total number of nodes in the AC power network; g ik and b Ik are the conductance and susceptance of the branch ik connected between node i and node k, and are collectively referred to as the admittance of branch ik; Is the square current of the branch ik connected between node i and node k; a 2i-1,2h-1 , a 2k-1,2h-1 , a 2i-1,2h , a 2k-1,2h , A 2i, 2h-1 , a 2k, 2h-1 , a 2i, 2h , a 2k, 2h are respectively the MP inverse matrix of the 2n ⁇ 2n-dimensional full-node admittance matrix in the 2i-1 row 2h-1 Column, row 2k-1,
  • step S105 the method for establishing a symmetric algebraic formula for obtaining the branch mean square current component of the source load according to the symmetric algebraic expression of the branch mean square current with respect to the source load power of the nodes in the entire network and the Shapley value theorem is specifically as follows: :
  • i, j, k, and h are the numbers of nodes in the AC power network, and they all belong to the set of continuous natural numbers ⁇ 1,2,...,n ⁇ ; n is the total number of nodes in the AC power network; Is the component of the source load connected to node j in the square current of the branch ik connected between node i and node k, referred to as the branch mean square current component of the source load; g ik and b ik are respectively connected to The conductance and susceptance of the branch ik between node i and node k are collectively referred to as the admittance of branch ik; a 2i-1,2h-1 , a 2k-1,2h-1 , a 2i-1,2h , a 2k-1,2h-1 , a 2k,2h-1 , a 2i,2h , a 2k,2h , a 2i-1,2h , a 2i-1,2h
  • the above-mentioned symmetrical algebraic calculation formula for the branch mean square current component of the source load applies to the mean square current of all branches and the source load power of all nodes in the AC power network.
  • the source load power of all nodes is treated equally.
  • the invention is a symmetric method to obtain the branch mean square current component of the source load in the AC power network.
  • the symmetrical algebraic formula is a full variable (not incremental) expression of the source load power, so it is accurate for a wide range of source load power changes. This symmetrical and precise relationship between the branch mean square current component and the source-load power solves the problems of time-consuming, unreliable, inaccurate, and low-efficiency existing AC power network security correction methods.
  • a computer-readable storage medium provided by an embodiment of the present invention is a medium storing a computer program.
  • the computer program may be a source code program, an object code program, an executable file, or some intermediate form.
  • the steps of the symmetric method for obtaining the branch current component of the source load in the AC power network as described in the above embodiment are realized.
  • the computer-readable storage medium may include any entity or device capable of carrying the computer program, such as a U disk, a mobile hard disk, an optical disk, a computer memory, a random access memory, and the like.

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

一种获取交流电力网中源荷的支路圴方电流分量的对称方法，首先根据交流电力网中节点源荷功率和支路导纳建立节点源荷功率关于节点平移电压和节点电压相位的线性表达式；再据此建立交流电力网稳态的线性对称模型；接着利用M-P逆矩阵建立全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式；然后根据该线性对称矩阵表达式建立支路均方电流关于全网节点源荷功率的对称代数表达式；最后利用Shapley值定理建立获取源荷的支路均方电流分量的对称代数计算式，从而实现交流电力网中源荷的支路均方电流分量的获取。这种源荷的支路均方电流分量为实现交流电力网安全校正、保障其安全运行提供了一种新的高效准确工具。

Description

获取交流电力网中源荷的支路圴方电流分量的对称方法 技术领域
本发明涉及电力工程领域，尤其涉及一种获取交流电力网中(电)源(负)荷的支路圴方电流分量的对称方法和计算机可读存储介质。
背景技术
在交流电网中，支路均方电流(有效值的平方)与源荷功率的深层次简洁精准关系，是高效实现交流电力网安全校正、保障其安全运行的关键。源荷的支路均方电流分量是一种深层次的实现交流电力网安全校正的新工具、亟待研发。
现有的交流电力网安全校正方法，要么先构建安全校正优化模型，再通过优化得到校正方案实现；要么先获取支路电流对源荷功率的灵敏度，再利用灵敏度表达的近似线性关系实现。前者因为是优化模型而不仅无法保障可靠得到安全校正方案、且优化求解的计算量总是很大，后者因为支路电流对源荷功率的灵敏度的局部线性特征而使安全校正不准确，继而导致反复校正。
因此，现有的交流电力网安全校正方法要么费时且不可靠、要么不准确和低效能。
发明内容
本发明实施例提供一种获取交流电力网中源荷的支路圴方电流分量的对称方法和计算机可读存储介质，旨在解决现有的交流电力网安全校正方法不可靠且效率低的问题。
本发明实施例第一方面提供了一种获取交流电力网中源荷的支路圴方电流分量的对称方法，包括：
根据交流电力网中节点源荷功率和支路导纳建立节点源荷功率关于节点平 移电压和节点电压相位的线性表达式；
根据所述节点源荷功率关于节点平移电压和节点电压相位的线性表达式建立交流电力网稳态的线性对称模型；
根据所述交流电力网稳态的线性对称模型，利用M-P逆矩阵建立全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式；
根据所述全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式建立支路均方电流关于全网节点源荷功率的对称代数表达式；
根据所述支路均方电流关于全网节点源荷功率的对称代数表达式和Shapley值定理建立获取源荷的支路均方电流分量的对称代数计算式。
本发明实施例第二方面提供了一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述获取交流电力网中源荷的支路圴方电流分量的对称方法的步骤。
上述获取交流电力网中源荷的支路圴方电流分量的对称方法在实施过程中根据支路均方电流关于全网节点源荷功率的对称代数表达式和Shapley值定理建立获取源荷的支路均方电流分量的对称代数计算式，实现交流电力网中源荷的支路均方电流分量的获取。一方面，由于源荷的支路均方电流分量的对称代数计算式适用交流电力网中全部支路的均方电流和全部节点的源荷功率、全部节点的源荷功率都被等同对待，源荷的支路均方电流分量因此对所有源荷都是对称的；另一方面，由于源荷的支路均方电流分量的对称代数计算式是源荷功率的全变量(而非增量)表达式，因此对源荷功率的大范围变化都准确，且减少了计算量。这种支路均方电流分量与源荷功率之间的对称精准关系解决了现有的交流电力网安全校正方法费时且不可靠、不准确和低效能的问题。
附图说明
为了更清楚地说明本发明实施例技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根 据这些附图获得其他的附图。
图1是本发明实施例提供的一种获取交流电力网中源荷的支路圴方电流分量的对称方法的实现流程图；
图2是本发明实施例提供的交流电力网通用模型的结构示意图。
具体实施方式
以下描述中，为了说明而不是为了限定，提出了诸如特定系统结构、技术之类的具体细节，以便透彻理解本发明实施例。然而，本领域的技术人员应当清楚，在没有这些具体细节的其它实施例中也可以实现本发明。在其它情况中，省略对众所周知的系统、装置、电路以及方法的详细说明，以免不必要的细节妨碍本发明的描述。
为了说明本发明所述的技术方案，下面通过具体实施例来进行说明。
请参见图1和图2，本发明实施例提供的一种获取交流电力网中源荷的支路圴方电流分量的对称方法，包括：
步骤S101,根据已知的交流电力网中节点源荷功率和支路导纳建立节点源荷功率关于节点平移电压和节点电压相位的线性表达式；
步骤S102,根据所述节点源荷功率关于节点平移电压和节点电压相位的线性表达式建立交流电力网稳态的线性对称模型；
步骤S103,根据所述交流电力网稳态的线性对称模型，利用M-P逆矩阵建立全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式；
步骤S104,根据所述全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式建立支路均方电流关于全网节点源荷功率的对称代数表达式；
步骤S105,根据所述支路均方电流关于全网节点源荷功率的对称代数表达式和Shapley值定理建立获取源荷的支路均方电流分量的对称代数计算式。
对交流电力网中全部支路均方电流和全部节点源荷功率都按照上述对称代 数计算式计算，即可得到全部节点源荷的支路均方电流分量，从而实现交流电力网中源荷的支路均方电流分量的获取。这种支路均方电流分量与源荷功率之间的对称精准关系解决了现有的交流电力网安全校正对称方法费时且不可靠、不准确和低效能的问题。
步骤S101中，所述根据所述交流电力网中节点源荷功率和支路导纳建立节点源荷功率关于节点平移电压和节点电压相位的线性表达式的方法具体为：
按照如下关系式建立节点源荷功率关于节点平移电压和节点电压相位的线性表达式：
Figure PCTCN2019073442-appb-000001
Figure PCTCN2019073442-appb-000002
其中，i和k均为交流电力网中节点的编号，且都属于连续自然数的集合{1,2,…,n}；n为所述交流电力网中节点的总个数；P i和Q i分别为接于节点i的源荷有功功率和源荷无功功率，且统称为节点i的源荷功率；所述P i等于接于节点i的电源有功功率减去负荷有功功率，所述Q i等于接于节点i的电源无功功率减去负荷无功功率；g ik和b ik分别是连接在节点i和节点k之间的支路ik的电导和电纳，且统称为支路ik的导纳；θ i和θ k分别为节点i和节点k的电压相位；v i和v k分别为节点i和节点k的平移电压，且都是平移-1.0后的标幺值电压。
步骤S102中，所述根据所述节点源荷功率关于节点平移电压和节点电压相位的线性表达式建立交流电力网稳态的线性对称模型的方法具体为：
按照如下关系式建立交流电力网稳态的线性对称模型：
[P 1Q 1…P iQ i…P nQ n] T＝(G *,*)[θ 1v 1…θ iv i…θ nv n] T
且(G *,*)先置零、再扫描支路按下式累加构建：G 2i-1,2i-1＝G 2i-1,2i-1-b ij，G 2i-1,2i＝G 2i-1,2i+g ij，G 2i-1,2j-1＝G 2i-1,2j-1+b ij，G 2i-1,2j＝G 2i-1,2j-g ij，G 2i,2i-1＝G 2i,2i-1-g ij，G 2i,2i＝G 2i,2i-b ij，G 2i,2j-1＝G 2i,2j-1+g ij， G 2i,2j＝G 2i,2j+b ij。
其中，i和j均为交流电力网中节点的编号，且都属于连续自然数的集合{1,2,…,n}；n为所述交流电力网中节点的总个数；P 1和Q 1分别为接于节点1的源荷有功功率和源荷无功功率，且统称为节点1的源荷功率；所述P 1等于接于节点1的电源有功功率减去负荷有功功率，所述Q 1等于接于节点1的电源无功功率减去负荷无功功率；P i和Q i分别为接于节点i的源荷有功功率和源荷无功功率，且统称为节点i的源荷功率；所述P i等于接于节点i的电源有功功率减去负荷有功功率，所述Q i等于接于节点i的电源无功功率减去负荷无功功率；P n和Q n分别为接于节点n的源荷有功功率和源荷无功功率，且统称为节点n的源荷功率；所述P n等于接于节点n的电源有功功率减去负荷有功功率，所述Q n等于接于节点n的电源无功功率减去负荷无功功率；g ij和b ij分别是连接在节点i和节点j之间的支路ij的电导和电纳，且统称为支路ij的导纳；θ 1、θ i和θ n分别为节点1、节点i和节点n的电压相位；v 1、v i和v n分别为节点1、节点i和节点n的平移电压，且都是平移-1.0后的标幺值电压；(G *,*)是2n×2n维全节点导纳矩阵；G 2i-1,2i-1、G 2i-1,2i、G 2i-1,2j-1、G 2i-1,2j、G 2i,2i-1、G 2i,2i、G 2i,2j-1、G 2i,2j都分别是所述全节点导纳矩阵(G *,*)中第2i-1行第2i-1列、第2i-1行第2i列、第2i-1行第2j-1列、第2i-1行第2j列、第2i行第2i-1列、第2i行第2i列、第2i行第2j-1列、第2i行第2j列的元素。
上述电力网稳态模型是线性的，且全部节点源荷功率都被列入该模型中、都被等同对待，这正是称之为线性对称模型的缘故。
步骤S103中,所述根据所述交流电力网稳态的线性对称模型，利用M-P逆矩阵建立全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式的方法具体为：
按照如下关系式建立全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式：
[θ 1v 1…θ iv i…θ nv n] T＝(a *,*)[P 1Q 1…P iQ i…P nQ n] T
(a *,*)＝(G *,*) +
其中，i为交流电力网中节点的编号，且属于连续自然数的集合{1,2,…,n}；n为所述交流电力网中节点的总个数；θ 1、θ i和θ n分别为节点1、节点i和节点n的电压相位；v 1、v i和v n分别为节点1、节点i和节点n的平移电压，且都是平移-1.0后的标幺值电压；P 1和Q 1分别为接于节点1的源荷有功功率和源荷无功功率，且统称为节点1的源荷功率；所述P 1等于接于节点1的电源有功功率减去负荷有功功率，所述Q 1等于接于节点1的电源无功功率减去负荷无功功率；P i和Q i分别为接于节点i的源荷有功功率和源荷无功功率，且统称为节点i的源荷功率；所述P i等于接于节点i的电源有功功率减去负荷有功功率，所述Q i等于接于节点i的电源无功功率减去负荷无功功率；P n和Q n分别为接于节点n的源荷有功功率和源荷无功功率，且统称为节点n的源荷功率；所述P n等于接于节点n的电源有功功率减去负荷有功功率，所述Q n等于接于节点n的电源无功功率减去负荷无功功率；(G *,*)是2n×2n维全节点导纳矩阵；上标符号+是求M-P逆矩阵的运算符；(a *,*)是所述全节点导纳矩阵(G *,*)的M-P逆矩阵。
步骤S104中，所述根据所述全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式建立支路均方电流关于全网节点源荷功率的对称代数表达式的方法具体为：
基于支路均方电流表达式常识：
Figure PCTCN2019073442-appb-000003
按照如下关系式建立支路均方电流关于全网节点源荷功率的对称代数表达式：
Figure PCTCN2019073442-appb-000004
其中，i、k和h均为交流电力网中节点的编号，且都属于连续自然数的集合{1,2,…,n}；n为所述交流电力网中节点的总个数；g ik和b ik分别是连接在节点i和节点k之间的支路ik的电导和电纳，且统称为支路ik的导纳；
Figure PCTCN2019073442-appb-000005
是连接在节点i和 节点k之间的支路ik的圴方电流；a 2i-1,2h-1、a 2k-1,2h-1、a 2i-1,2h、a 2k-1,2h、a 2i,2h-1、a 2k,2h-1、a 2i,2h、a 2k,2h都分别是2n×2n维全节点导纳矩阵的M-P逆矩阵中第2i-1行第2h-1列、第2k-1行第2h-1列、第2i-1行第2h列、第2k-1行第2h列、第2i行第2h-1列、第2k行第2h-1列、第2i行第2h列、第2k行第2h列的元素；P h和Q h分别为接于节点h的源荷有功功率和源荷无功功率，且统称为节点h的源荷功率；所述P h等于接于节点h的电源有功功率减去负荷有功功率，所述Q h等于接于节点h的电源无功功率减去负荷无功功率。
步骤S105中，所述根据所述支路均方电流关于全网节点源荷功率的对称代数表达式和Shapley值定理建立获取源荷的支路均方电流分量的对称代数计算式的方法具体为：
按照如下关系式建立获取源荷的支路均方电流分量的对称代数计算式：
Figure PCTCN2019073442-appb-000006
其中，i、j、k和h均为交流电力网中节点的编号，且都属于连续自然数的集合{1,2,…,n}；n为所述交流电力网中节点的总个数；
Figure PCTCN2019073442-appb-000007
是连接在节点i和节点k之间的支路ik的圴方电流中归属接于节点j的源荷的分量，简称源荷的支路均方电流分量；g ik和b ik分别是连接在节点i和节点k之间的支路ik的电导和电纳，且统称为支路ik的导纳；a 2i-1,2h-1、a 2k-1,2h-1、a 2i-1,2h、a 2k-1,2h、a 2i,2h-1、a 2k,2h-1、a 2i,2h、a 2k,2h、a 2i-1,2j-1、a 2k-1,2j-1、a 2i-1,2j、a 2k-1,2j、a 2i,2j-1、a 2k,2j-1、a 2i,2j、a 2k,2j都分别是2n×2n维全节点导纳矩阵的M-P逆矩阵中第2i-1行第2h-1列、第2k-1行第2h-1列、第2i-1行第2h列、第2k-1行第2h列、第2i行第2h-1列、第2k行第2h-1列、第2i行第2h列、第2k行第2h列、第2i-1行第2j-1列、第2k-1行第2j-1列、第2i-1行第2j列、第2k-1行第2j列、第2i行第2j-1列、 第2k行第2j-1列、第2i行第2j列、第2k行第2j列的元素；P h和Q h分别为接于节点h的源荷有功功率和源荷无功功率，且统称为节点h的源荷功率；所述P h等于接于节点h的电源有功功率减去负荷有功功率，所述Q h等于接于节点h的电源无功功率减去负荷无功功率；P j和Q j分别为接于节点j的源荷有功功率和源荷无功功率，且统称为节点j的源荷功率；所述P j等于接于节点j的电源有功功率减去负荷有功功率，所述Q j等于接于节点j的电源无功功率减去负荷无功功率。
上述源荷的支路均方电流分量的对称代数计算式适用交流电力网中全部支路的均方电流和全部节点的源荷功率，全部节点的源荷功率都被等同对待，这正是称本发明为获取交流电力网中源荷的支路均方电流分量的对称方法的缘故。此外，该对称代数计算式是源荷功率的全变量(而非增量)表达式，它因此对源荷功率的大范围变化都准确。这种支路均方电流分量与源荷功率之间的对称精准关系解决了现有的交流电力网安全校正方法费时且不可靠、不准确和低效能的问题。
本发明实施例提供的一种计算机可读存储介质，是存储有计算机程序的介质。所述计算机程序可以为源代码程序、对象代码程序、可执行文件或某些中间形式等。所述计算机程序被处理器执行时实现如上实施例所述获取交流电力网中源荷的支路圴方电流分量的对称方法的步骤。所述计算机可读存储介质可以包括能够携带所述计算机程序的任何实体或装置，例如U盘、移动硬盘、光盘、计算机存储器、随机存取存储器等。
以上所述实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围，均应包含在本发明的保护范围之内。

Claims (7)

  1. 一种获取交流电力网中源荷的支路圴方电流分量的对称方法，其特征在于，包括：
    根据交流电力网中节点源荷功率和支路导纳建立节点源荷功率关于节点平移电压和节点电压相位的线性表达式；
    根据所述节点源荷功率关于节点平移电压和节点电压相位的线性表达式建立交流电力网稳态的线性对称模型；
    根据所述交流电力网稳态的线性对称模型，利用M-P逆矩阵建立全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式；
    根据所述全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式建立支路均方电流关于全网节点源荷功率的对称代数表达式；
    根据所述支路均方电流关于全网节点源荷功率的对称代数表达式和Shapley值定理建立获取源荷的支路均方电流分量的对称代数计算式。
  2. 根据权利要求1所述的获取交流电力网中源荷的支路圴方电流分量的对称方法，其特征在于，所述根据所述交流电力网中节点源荷功率和支路导纳建立节点源荷功率关于节点平移电压和节点电压相位的线性表达式的方法具体为：
    按照如下关系式建立节点源荷功率关于节点平移电压和节点电压相位的线性表达式：
    Figure PCTCN2019073442-appb-100001
    Figure PCTCN2019073442-appb-100002
    其中，i和k均为交流电力网中节点的编号，且都属于连续自然数的集合{1,2,…,n}；n为所述交流电力网中节点的总个数；P i和Q i分别为接于节点i的源荷有功功率和源荷无功功率，且统称为节点i的源荷功率；g ik和b ik分别是连接在节点i和节点k之间的支路ik的电导和电纳，且统称为支路ik的导纳；θ i和θ k分别为节点i和节点k的电压相位；υ i和υ k分别为节点i和节点k的平移电压，且都是 平移-1.0后的标幺值电压。
  3. 根据权利要求1所述的获取交流电力网中源荷的支路圴方电流分量的对称方法，其特征在于，所述根据所述节点源荷功率关于节点平移电压和节点电压相位的线性表达式建立交流电力网稳态的线性对称模型的方法具体为：
    按照如下关系式建立交流电力网稳态的线性对称模型：
    [P 1Q 1…P iQ i…P nQ n] T＝(G *,*)[θ 1υ 1…θ iυ i…θ nυ n] T
    且(G *,*)先置零、再扫描支路按下式累加构建：G 2i-1,2i-1＝G 2i-1,2i-1-b ij，G 2i-1,2i＝G 2i-1,2i+g ij，G 2i-1,2j-1＝G 2i-1,2j-1+b ij，G 2i-1,2j＝G 2i-1,2j-g ij，G 2i,2i-1＝G 2i,2i-1-g ij，G 2i,2i＝G 2i,2i-b ij，G 2i,2j-1＝G 2i,2j-1+g ij，G 2i,2j＝G 2i,2j+b ij；
    其中，i和j均为交流电力网中节点的编号，且都属于连续自然数的集合{1,2,…,n}；n为所述交流电力网中节点的总个数；P 1和Q 1分别为接于节点1的源荷有功功率和源荷无功功率，且统称为节点1的源荷功率；P i和Q i分别为接于节点i的源荷有功功率和源荷无功功率，且统称为节点i的源荷功率；P n和Q n分别为接于节点n的源荷有功功率和源荷无功功率，且统称为节点n的源荷功率；g ij和b ij分别是连接在节点i和节点j之间的支路ij的电导和电纳，且统称为支路ij的导纳；θ 1、θ i和θ n分别为节点1、节点i和节点n的电压相位；υ 1、υ i和υ n分别为节点1、节点i和节点n的平移电压，且都是平移-1.0后的标幺值电压；(G *,*)是2n×2n维全节点导纳矩阵；G 2i-1,2i-1、G 2i-1,2i、G 2i-1,2j-1、G 2i-1,2j、G 2i,2i-1、G 2i,2i、G 2i,2j-1、G 2i,2j都是所述全节点导纳矩阵(G *,*)中的元素。
  4. 根据权利要求1所述的获取交流电力网中源荷的支路圴方电流分量的对称方法，其特征在于，所述根据所述交流电力网稳态的线性对称模型，利用M-P逆矩阵建立全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式的方法具体为：
    按照如下关系式建立全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式：
    [θ 1υ 1…θ iυ i…θ nυ n] T＝(a *,*)[P 1Q 1…P iQ i…P nQ n] T
    (a *,*)＝(G *,*) +
    其中，i为交流电力网中节点的编号，且属于连续自然数的集合{1,2,…,n}；n为所述交流电力网中节点的总个数；θ 1、θ i和θ n分别为节点1、节点i和节点n的电压相位；υ 1、υ i和υ n分别为节点1、节点i和节点n的平移电压，且都是平移-1.0后的标幺值电压；P 1和Q 1分别为接于节点1的源荷有功功率和源荷无功功率，且统称为节点1的源荷功率；P i和Q i分别为接于节点i的源荷有功功率和源荷无功功率，且统称为节点i的源荷功率；P n和Q n分别为接于节点n的源荷有功功率和源荷无功功率，且统称为节点n的源荷功率；(G *,*)是2n×2n维全节点导纳矩阵；上标符号+是求M-P逆矩阵的运算符；(a *,*)是所述全节点导纳矩阵(G *,*)的M-P逆矩阵。
  5. 根据权利要求1所述的获取交流电力网中源荷的支路圴方电流分量的对称方法，其特征在于，所述根据所述全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式建立支路均方电流关于全网节点源荷功率的对称代数表达式的方法具体为：
    基于支路均方电流表达式常识：
    Figure PCTCN2019073442-appb-100003
    按照如下关系式建立支路均方电流关于全网节点源荷功率的对称代数表达式：
    Figure PCTCN2019073442-appb-100004
    其中，i、k和h均为交流电力网中节点的编号，且都属于连续自然数的集合{1,2,…,n}；n为所述交流电力网中节点的总个数；g ik和b ik分别是连接在节点i和节点k之间的支路ik的电导和电纳，且统称为支路ik的导纳；
    Figure PCTCN2019073442-appb-100005
    是连接在节点i和节点k之间的支路ik的圴方电流；a 2i-1，2h-1、a 2k-1,2h-1、a 2i-1,2h、a 2k-1,2h、a 2i,2h-1、a 2k,2h-1、a 2i,2h、a 2k,2h都是2n×2n维全节点导纳矩阵的M-P逆矩阵中的元素；P h和Q h分别为接于节点h的源荷有功功率和源荷无功功率，且统称为节点h的源 荷功率。
  6. 根据权利要求1所述的获取交流电力网中源荷的支路圴方电流分量的对称方法，其特征在于，所述根据所述支路均方电流关于全网节点源荷功率的对称代数表达式和Shapley值定理建立获取源荷的支路均方电流分量的对称代数计算式的方法具体为：
    按照如下关系式建立获取源荷的支路均方电流分量的对称代数计算式：
    Figure PCTCN2019073442-appb-100006
    其中，i、j、k和h均为交流电力网中节点的编号，且都属于连续自然数的集合{1,2,…,n}；n为所述交流电力网中节点的总个数；
    Figure PCTCN2019073442-appb-100007
    是连接在节点i和节点k之间的支路ik的圴方电流中归属接于节点j的源荷的分量，简称源荷的支路均方电流分量；g ik和b ik分别是连接在节点i和节点k之间的支路ik的电导和电纳，且统称为支路ik的导纳；a 2i-1,2h-1、a 2k-1,2h-1、a 2i-1,2h、a 2k-1,2h、a 2i,2h-1、a 2k,2h-1、a 2i,2h、a 2k,2h、a 2i-1，2j-1、a 2k-1,2j-1、a 2i-1,2j、a 2k-1,2j、a 2i,2j-1、a 2k,2j-1、a 2i,2j、a 2k,2j都是2n×2n维全节点导纳矩阵的M-P逆矩阵中的元素；P h和Q h分别为接于节点h的源荷有功功率和源荷无功功率，且统称为节点h的源荷功率；P j和Q j分别为接于节点j的源荷有功功率和源荷无功功率，且统称为节点j的源荷功率。
  7. 一种计算机可读存储介质，所述计算机可读存储介质存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现如权利要求1至6任一项所述获取交流电力网中源荷的支路圴方电流分量的对称方法的步骤。
PCT/CN2019/073442 2019-01-28 2019-01-28 获取交流电力网中源荷的支路圴方电流分量的对称方法 Ceased WO2020154849A1 (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
US16/630,304 US20210384728A1 (en) 2019-01-28 2019-01-28 Symmetric method for obtaining branch-mean-square-current components induced by sources and loads at individual buses in ac power networks
PCT/CN2019/073442 WO2020154849A1 (zh) 2019-01-28 2019-01-28 获取交流电力网中源荷的支路圴方电流分量的对称方法
CN201980002743.3A CN111758195B (zh) 2019-01-28 2019-01-28 获取交流电力网中源荷的支路均方电流分量的对称方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2019/073442 WO2020154849A1 (zh) 2019-01-28 2019-01-28 获取交流电力网中源荷的支路圴方电流分量的对称方法

Publications (1)

Publication Number Publication Date
WO2020154849A1 true WO2020154849A1 (zh) 2020-08-06

Family

ID=71839905

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2019/073442 Ceased WO2020154849A1 (zh) 2019-01-28 2019-01-28 获取交流电力网中源荷的支路圴方电流分量的对称方法

Country Status (3)

Country Link
US (1) US20210384728A1 (zh)
CN (1) CN111758195B (zh)
WO (1) WO2020154849A1 (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6904372B2 (en) * 2001-12-20 2005-06-07 Abb Technology Ag Method of the identification of weak and/or strong nodes of an electric power system
CN103795063A (zh) * 2014-03-01 2014-05-14 华北电力大学 基于源荷协同系数的线路过负荷紧急控制系统及其方法
CN104362619A (zh) * 2014-10-21 2015-02-18 深圳大学 交流电力网中支路均方电流的策动分量的获取方法
CN104995810A (zh) * 2014-11-18 2015-10-21 深圳大学 交流电力网中源荷同变的对称功率传输系数的获取方法

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2446943C2 (de) * 1974-10-02 1982-06-16 Brown, Boveri & Cie Ag, 6800 Mannheim Schaltungsanordnung zur schnellen Erfassung und Kompensation der Verbraucherblindleistungsaufnahme in ein- oder mehrphasigen Wechselstromnetzen
WO2018209477A1 (zh) * 2017-05-15 2018-11-22 深圳大学 获取直流电力网功率传输系数的均衡电导补偿型偏心方法
CN109478781B (zh) * 2017-05-15 2021-08-31 深圳大学 获取直流电力网功率传输系数的均衡电导补偿型对称方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6904372B2 (en) * 2001-12-20 2005-06-07 Abb Technology Ag Method of the identification of weak and/or strong nodes of an electric power system
CN103795063A (zh) * 2014-03-01 2014-05-14 华北电力大学 基于源荷协同系数的线路过负荷紧急控制系统及其方法
CN104362619A (zh) * 2014-10-21 2015-02-18 深圳大学 交流电力网中支路均方电流的策动分量的获取方法
CN104995810A (zh) * 2014-11-18 2015-10-21 深圳大学 交流电力网中源荷同变的对称功率传输系数的获取方法

Also Published As

Publication number Publication date
CN111758195B (zh) 2023-09-15
CN111758195A (zh) 2020-10-09
US20210384728A1 (en) 2021-12-09

Similar Documents

Publication Publication Date Title
Zamzam et al. Beyond relaxation and Newton–Raphson: Solving AC OPF for multi-phase systems with renewables
Jain et al. Single‐phase single‐stage multifunctional grid interfaced solar photo‐voltaic system under abnormal grid conditions
Farivar et al. Local voltage control in distribution systems: An incremental control algorithm
Zaheeruddin et al. Intelligent fractional-order-based centralized frequency controller for microgrid
Baradar et al. AC power flow representation in conic format
Bunin et al. Experimental real-time optimization of a solid oxide fuel cell stack via constraint adaptation
CN105745809B (zh) 获取多端直流电力网非线性有功潮流的对称方法
CN107210605B (zh) 获取多端直流电力网线性有功潮流的对称方法
CN113452028B (zh) 低压配电网概率潮流计算方法、系统、终端和存储介质
Weng et al. Fixed-point theorem-based voltage stability margin estimation techniques for distribution systems with renewables
Kamel et al. Modeling and analysis of voltage and power control devices in current injections load flow method
CN114977189B (zh) 一种基于DistFlow线性潮流的配电网潮流计算方法及系统
CN108347057A (zh) 一种含lcc-mmc混合直流电网的交直流电力系统潮流计算方法
CN104680249B (zh) 一种基于改进离散粒子群算法的pmu配置方法
Kumar et al. A non-iterative technique for determination of solar cell parameters from the light generated IV characteristic
CN102780220A (zh) 含有pv恒定型分布式电源的配电网潮流计算方法
WO2020154847A1 (zh) 获取交流电力网中源荷的网损功率分量的对称方法
CN106877309A (zh) 主从控制模式下直流微电网的稳定性判定方法
WO2020154850A1 (zh) 获取交流电力网中源荷的节点平移电压分量的对称方法
CN121150165A (zh) 一种考虑暂态过电压约束的lcc-hvdc运行范围优化方法及系统
WO2020154849A1 (zh) 获取交流电力网中源荷的支路圴方电流分量的对称方法
CN105068590A (zh) 光伏发电功率控制方法和系统
Ebeed et al. Constraints violation handling of GUPFC in Newton-Raphson power flow
CN119474648A (zh) 一种电力系统潮流计算方法及系统
CN117194458A (zh) 一种电力碳强度的更新方法、装置、设备及存储介质

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 19913439

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 19913439

Country of ref document: EP

Kind code of ref document: A1

32PN Ep: public notification in the ep bulletin as address of the adressee cannot be established

Free format text: NOTING OF LOSS OF RIGHTS PURSUANT TO RULE 112(1) EPC (EPO FORM 1205A DATED 25/01/2022)

122 Ep: pct application non-entry in european phase

Ref document number: 19913439

Country of ref document: EP

Kind code of ref document: A1