CN111758197B - 获取交流电力网中源荷的网损功率分量的对称方法 - Google Patents
获取交流电力网中源荷的网损功率分量的对称方法 Download PDFInfo
- Publication number
- CN111758197B CN111758197B CN201980002748.6A CN201980002748A CN111758197B CN 111758197 B CN111758197 B CN 111758197B CN 201980002748 A CN201980002748 A CN 201980002748A CN 111758197 B CN111758197 B CN 111758197B
- Authority
- CN
- China
- Prior art keywords
- node
- power
- network
- source load
- alternating current
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 239000011159 matrix material Substances 0.000 claims abstract description 36
- 238000004364 calculation method Methods 0.000 claims abstract description 12
- 238000004590 computer program Methods 0.000 claims description 6
- 238000009825 accumulation Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—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/007—Arrangements for selectively connecting the load or loads to one or several among a plurality of power lines or power sources
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—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/001—Methods to deal with contingencies, e.g. abnormalities, faults or failures
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—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/04—Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/50—Controlling the sharing of the out-of-phase component
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S40/00—Systems for electrical power generation, transmission, distribution or end-user application management characterised by the use of communication or information technologies, or communication or information technology specific aspects supporting them
- Y04S40/20—Information technology specific aspects, e.g. CAD, simulation, modelling, system security
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中,所述根据所述交流电力网中节点源荷功率和支路导纳建立节点源荷功率关于节点平移电压和节点电压相位的线性表达式的方法具体为:
按照如下关系式建立节点源荷功率关于节点平移电压和节点电压相位的线性表达式:
其中,i和k均为交流电力网中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为所述交流电力网中节点的总个数;Pi和Qi分别为接于节点i的源荷有功功率和源荷无功功率,且统称为节点i的源荷功率;所述Pi等于接于节点i的电源有功功率减去负荷有功功率,所述Qo等于接于节点i的电源无功功率减去负荷无功功率;gik和bik分别是连接在节点i和节点k之间的支路ik的电导和电纳,且统称为支路ik的导纳;θi和θk分别为节点i和节点k的电压相位;υi和υk分别为节点i和节点k的平移电压,且都是平移-1.0后的标幺值电压。
步骤S102中,所述根据所述节点源荷功率关于节点平移电压和节点电压相位的线性表达式建立交流电力网稳态的线性对称模型的方法具体为:
按照如下关系式建立交流电力网稳态的线性对称模型:
[P1Q1…PiQi…PnQn]T=(G*,*)[θ1υ1…θiυi…θnυn]T
且(G*,*)先置零、再扫描支路按下式累加构建:G2i-1,2i-1=G2i-1,2i-1-bij,G2i-1,2i=G2i-1,2i+gij,G2i-1,2j-1=G2i-1,2j-1+bij,G2i-1,2j=G2i-1,2j-gij,G2i,2i-1=G2i,2i-1-gij,G2i,2i=G2i,2i-bij,G2i,2j-1=G2i,2j-1+gij,G2i,2j=G2i,2j+bij。
其中,i和j均为交流电力网中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为所述交流电力网中节点的总个数;P1和Q1分别为接于节点1的源荷有功功率和源荷无功功率,且统称为节点1的源荷功率;所述P1等于接于节点1的电源有功功率减去负荷有功功率,所述Q1等于接于节点1的电源无功功率减去负荷无功功率;Pi和Qi分别为接于节点i的源荷有功功率和源荷无功功率,且统称为节点i的源荷功率;所述Pi等于接于节点i的电源有功功率减去负荷有功功率,所述Qi等于接于节点i的电源无功功率减去负荷无功功率;Pn和Qn分别为接于节点n的源荷有功功率和源荷无功功率,且统称为节点n的源荷功率;所述Pn等于接于节点n的电源有功功率减去负荷有功功率,所述Qn等于接于节点n的电源无功功率减去负荷无功功率;gij和bij分别是连接在节点i和节点j之间的支路ij的电导和电纳,且统称为支路ij的导纳;θ1、θi和θn分别为节点1、节点i和节点n的电压相位;v1、υi和υn分别为节点1、节点i和节点n的平移电压,且都是平移-1.0后的标幺值电压;(G*,*)是2n×2n维全节点导纳矩阵;G2i-1,2i-1、G2i-1,2i、G2i-1,2j-1、G2i-1,2j、G2i,2i-1、G2i,2i、G2i,2j-1、G2i,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逆矩阵建立全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式的方法具体为:
按照如下关系式建立全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式:
[θ1v1…θiυi…θnυn]T=(a*,*)[P1Q1…PiQi…PnQn]T
(a*,*)=(G*,*)+
其中,i为交流电力网中节点的编号,且属于连续自然数的集合{1,2,…,n};n为所述交流电力网中节点的总个数;θ1、θi和θn分别为节点1、节点i和节点n的电压相位;υ1、υi和vn分别为节点1、节点i和节点n的平移电压,且都是平移-1.0后的标幺值电压;P1和Q1分别为接于节点1的源荷有功功率和源荷无功功率,且统称为节点1的源荷功率;所述P1等于接于节点1的电源有功功率减去负荷有功功率,所述Q1等于接于节点1的电源无功功率减去负荷无功功率;Pi和Qi分别为接于节点i的源荷有功功率和源荷无功功率,且统称为节点i的源荷功率;所述Pi等于接于节点i的电源有功功率减去负荷有功功率,所述Qi等于接于节点i的电源无功功率减去负荷无功功率;Pn和Qn分别为接于节点n的源荷有功功率和源荷无功功率,且统称为节点n的源荷功率;所述Pn等于接于节点n的电源有功功率减去负荷有功功率,所述Qn等于接于节点n的电源无功功率减去负荷无功功率;(G*,*)是2n×2n维全节点导纳矩阵;上标符号+是求M-P逆矩阵的运算符;(a*,*)是所述全节点导纳矩阵(G*,*)的M-P逆矩阵。
步骤S104中,所述根据所述全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式建立网损功率关于全网节点源荷功率的对称代数表达式的方法具体为:
基于网损功率表达式:PL=∑ik∈Ωgik[(θi-θk)2+(υi-υk)2],按照如下关系式建立网损功率关于全网节点源荷功率的对称代数表达式:
其中,i、k和h均为交流电力网中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为所述交流电力网中节点的总个数;ik是连接在节点i和节点k之间的支路;Ω是交流电力网中所有支路构成的集合;gik和bik分别是连接在节点i和节点k之间的支路ik的电导和电纳,且统称为支路ik的导纳;PL是交流电力网的网损功率;a2i-1,2h-1、a2k-1,2h-1、a2i-1,2h、a2k-1,2h、a2i,2h-1、a2k,2h-1、a2i,2h、a2k,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列的元素;Ph和Qh分别为接于节点h的源荷有功功率和源荷无功功率,且统称为节点h的源荷功率;所述Ph等于接于节点h的电源有功功率减去负荷有功功率,所述Qh等于接于节点h的电源无功功率减去负荷无功功率。
步骤S105中,所述根据所述网损功率关于全网节点源荷功率的对称代数表达式和Shapley值定理建立获取源荷的网损功率分量的对称代数计算式的方法具体为:
按照如下关系式建立源荷的网损功率分量的对称代数计算式:
其中,i、j、k和h均为交流电力网中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为所述交流电力网中节点的总个数;ik是连接在节点i和节点k之间的支路;Ω是交流电力网中所有支路构成的集合;PLj是交流电力网的网损功率中归属接于节点j的源荷的分量,简称源荷的网损功率分量;gik和bik分别是连接在节点i和节点k之间的支路ik的电导和电纳,且统称为支路ik的导纳;a2i-1,2h-1、a2k-1,2h-1、a2i-1,2h、a2k-1,2h、a2i,2h-1、a2k,2h-1、a2i,2h、a2k,2h、a2i-1,2j-1、a2k-1,2j-1、a2i-1,2j、a2k-1,2j、a2i,2j-1、a2k,2j-1、a2i,2j、a2k,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列的元素;Ph和Qh分别为接于节点h的源荷有功功率和源荷无功功率,且统称为节点h的源荷功率;所述Ph等于接于节点h的电源有功功率减去负荷有功功率,所述Qh等于接于节点h的电源无功功率减去负荷无功功率;Pj和Qj分别为接于节点j的源荷有功功率和源荷无功功率,且统称为节点j的源荷功率;所述Pj等于接于节点j的电源有功功率减去负荷有功功率,所述Qk等于接于节点j的电源无功功率减去负荷无功功率。
上述源荷的网损功率分量的对称代数计算式适用交流电力网中全部节点的源荷功率,全部节点的源荷功率都被等同对待,这正是称本发明为获取交流电力网中源荷的网损功率分量的对称方法的缘故。此外,该对称代数计算式是源荷功率的全变量(而非增量)表达式,它因此对源荷功率的大范围变化都准确。这种网损功率分量与源荷功率之间的对称精准关系解决了现有的获取交流电力网调度中网损功率的方法不能准确跟踪源荷功率的动态变化、无法计入无功影响和结果不唯一的问题。
本发明实施例提供的一种计算机可读存储介质,是存储有计算机程序的介质。所述计算机程序可以为源代码程序、对象代码程序、可执行文件或某些中间形式等。所述计算机程序被处理器执行时实现如上实施例所述获取交流电力网中源荷的网损功率分量的对称方法的步骤。所述计算机可读存储介质可以包括能够携带所述计算机程序的任何实体或装置,例如U盘、移动硬盘、光盘、计算机存储器、随机存取存储器等。
以上所述实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围,均应包含在本发明的保护范围之内。
Claims (2)
1.一种获取交流电力网中源荷的网损功率分量的对称方法,其特征在于,包括:
根据交流电力网中节点源荷功率和支路导纳建立节点源荷功率关于节点平移电压和节点电压相位的线性表达式;其中,按照如下关系式建立节点源荷功率关于节点平移电压和节点电压相位的线性表达式:
其中,i和k均为交流电力网中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为所述交流电力网中节点的总个数;Pi和Qi分别为接于节点i的源荷有功功率和源荷无功功率,且统称为节点i的源荷功率;gik和bik分别是连接在节点i和节点k之间的支路ik的电导和电纳,且统称为支路ik的导纳;θi和θk分别为节点i和节点k的电压相位;υi和υk分别为节点i和节点k的平移电压,且都是平移-1.0后的标幺值电压;
根据所述节点源荷功率关于节点平移电压和节点电压相位的线性表达式建立交流电力网稳态的线性对称模型;按照如下关系式建立交流电力网稳态的线性对称模型:
[P1 Q1…Pi Qi…Pn Qn]T=(G*,*)[θ1 v1…θi vi…θn un]T
且(G*,*)先置零、再扫描支路按下式累加构建:G2i-1,2i-1=G2i-1,2i-1-bij,G2i-1,2i=G2i-1,2i+gij,G2i-1,2j-1=G2i-1,2j-1+bij,G2i-1,2j=G2i-1,2j-gij,G2i,2i-1=G2i,2i-1-gij,G2i,2i=G2i,2i-bij,G2i,2j-1=G2i,2j-1+gij,G2i,2j=G2i,2j+bij;
其中,i和j均为交流电力网中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为所述交流电力网中节点的总个数;P1和Q1分别为接于节点1的源荷有功功率和源荷无功功率,且统称为节点1的源荷功率;Pi和Qi分别为接于节点i的源荷有功功率和源荷无功功率,且统称为节点i的源荷功率;Pn和Qn分别为接于节点n的源荷有功功率和源荷无功功率,且统称为节点n的源荷功率;gij和bij分别是连接在节点i和节点j之间的支路ij的电导和电纳,且统称为支路ij的导纳;θ1、θi和θn分别为节点1、节点i和节点n的电压相位;υ1、υi和υn分别为节点1、节点i和节点n的平移电压,且都是平移-1.0后的标幺值电压;(G*,*)是2n×2n维全节点导纳矩阵;G2i-1,2i-1、G2i-1,2i、G2i-1,2j-1、G2i-1,2j、G2i,2i-1、G2i,2i、G2i,2j-1、G2i,2j都是所述全节点导纳矩阵(G*,*)中的元素;
根据所述交流电力网稳态的线性对称模型,利用M-P逆矩阵建立全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式;其中,按照如下关系式建立全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式:
[θ1 v1…θi vi…θn vn]T=(a*,*)[P1 Q1…Pi Qi…Pn Qn]T
(a*,*)=(G*,*)+
其中,i为交流电力网中节点的编号,且属于连续自然数的集合{1,2,…,n};n为所述交流电力网中节点的总个数;θ1、θi和θn分别为节点1、节点i和节点n的电压相位;v1、υi和υn分别为节点1、节点i和节点n的平移电压,且都是平移-1.0后的标幺值电压;P1和Q1分别为接于节点1的源荷有功功率和源荷无功功率,且统称为节点1的源荷功率;Pi和Qi分别为接于节点i的源荷有功功率和源荷无功功率,且统称为节点i的源荷功率;Pn和Qn分别为接于节点n的源荷有功功率和源荷无功功率,且统称为节点n的源荷功率;(G*,*)是2n×2n维全节点导纳矩阵;上标符号+是求M-P逆矩阵的运算符;(a*,*)是所述全节点导纳矩阵(G*,*)的M-P逆矩阵;
根据所述全网节点平移电压和节点电压相位关于全网节点源荷功率的线性对称矩阵表达式建立网损功率关于全网节点源荷功率的对称代数表达式;其中,基于网损功率表达式:PL=∑ik∈Ωgik[(θi-θk)2+(υi-υk)2],按照如下关系式建立网损功率关于全网节点源荷功率的对称代数表达式:
其中,i、k和h均为交流电力网中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为所述交流电力网中节点的总个数;ik是连接在节点i和节点k之间的支路;Ω是交流电力网中所有支路构成的集合;gik和bik分别是连接在节点i和节点k之间的支路ik的电导和电纳,且统称为支路ik的导纳;PL是交流电力网的网损功率;a2i-1,2h-1、a2k-1,2h-1、a2i-1,2h、a2k-1,2h、a2i,2h-1、a2k,2h-1、a2i,2h、a2k,2h都是2n×2n维全节点导纳矩阵的M-P逆矩阵中的元素;Ph和Qh分别为接于节点h的源荷有功功率和源荷无功功率,且统称为节点h的源荷功率;
根据所述网损功率关于全网节点源荷功率的对称代数表达式和Shapley值定理建立获取源荷的网损功率分量的对称代数计算式;其中,按照如下关系式建立源荷的网损功率分量的对称代数计算式:
其中,i、j、k和h均为交流电力网中节点的编号,且都属于连续自然数的集合{1,2,…,n};n为所述交流电力网中节点的总个数;ik是连接在节点i和节点k之间的支路;Ω是交流电力网中所有支路构成的集合;PLj是交流电力网的网损功率中归属接于节点j的源荷的分量,简称源荷的网损功率分量;gik和bik分别是连接在节点i和节点k之间的支路ik的电导和电纳,且统称为支路ik的导纳;a2i-1,2h-1、a2k-1,2h-1、a2i-1,2h、a2k-1,2h、a2i,2h-1、a2k,2h-1、a2i,2h、a2k,2h、a2i-1,2j-1、a2k-1,2j-1、a2i-1,2j、a2k-1,2j、a2i,2j-1、a2k,2j-1、a2i,2j、a2k,2j都是2n×2n维全节点导纳矩阵的M-P逆矩阵中的元素;Ph和Qh分别为接于节点h的源荷有功功率和源荷无功功率,且统称为节点h的源荷功率;Pj和Qj分别为接于节点j的源荷有功功率和源荷无功功率,且统称为节点j的源荷功率。
2.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1所述获取交流电力网中源荷的网损功率分量的对称方法的步骤。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2019/073440 WO2020154847A1 (zh) | 2019-01-28 | 2019-01-28 | 获取交流电力网中源荷的网损功率分量的对称方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111758197A CN111758197A (zh) | 2020-10-09 |
CN111758197B true CN111758197B (zh) | 2024-04-02 |
Family
ID=71839848
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201980002748.6A Active CN111758197B (zh) | 2019-01-28 | 2019-01-28 | 获取交流电力网中源荷的网损功率分量的对称方法 |
Country Status (3)
Country | Link |
---|---|
US (1) | US11165251B2 (zh) |
CN (1) | CN111758197B (zh) |
WO (1) | WO2020154847A1 (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111952967B (zh) * | 2020-08-11 | 2022-07-08 | 广东电网有限责任公司广州供电局 | 一种多微电网系统停电故障恢复方法、系统及设备 |
CN114362160B (zh) * | 2022-01-11 | 2024-01-30 | 浙江华云电力工程设计咨询有限公司 | 一种基于比例均分原则的碳排放及绿电追踪方法、装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103580032A (zh) * | 2012-07-20 | 2014-02-12 | 比亚迪股份有限公司 | 电网补偿系统及其控制方法 |
CN104319772A (zh) * | 2014-10-21 | 2015-01-28 | 深圳大学 | 交直流混联电力网的损耗功率灵敏度的对称获取方法 |
CN104995810A (zh) * | 2014-11-18 | 2015-10-21 | 深圳大学 | 交流电力网中源荷同变的对称功率传输系数的获取方法 |
WO2016077995A1 (zh) * | 2014-11-18 | 2016-05-26 | 深圳大学 | 交流电力网的发电和用电有功实时调度的对称方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6940735B2 (en) * | 2003-11-14 | 2005-09-06 | Ballard Power Systems Corporation | Power converter system |
US7619327B2 (en) * | 2007-07-06 | 2009-11-17 | Hamilton Sundstrand Corporation | Hybrid electromechanical power transfer system |
CN102403724B (zh) * | 2011-11-09 | 2013-10-16 | 深圳大学 | 交直流混联电力网中节点电压灵敏度的对称获取方法 |
CN103956733B (zh) * | 2014-04-25 | 2016-05-18 | 深圳大学 | 电力网中节点到支路的有功功率传输系数的对称获取方法 |
-
2019
- 2019-01-28 CN CN201980002748.6A patent/CN111758197B/zh active Active
- 2019-01-28 US US16/630,253 patent/US11165251B2/en active Active
- 2019-01-28 WO PCT/CN2019/073440 patent/WO2020154847A1/zh active Application Filing
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103580032A (zh) * | 2012-07-20 | 2014-02-12 | 比亚迪股份有限公司 | 电网补偿系统及其控制方法 |
CN104319772A (zh) * | 2014-10-21 | 2015-01-28 | 深圳大学 | 交直流混联电力网的损耗功率灵敏度的对称获取方法 |
CN104995810A (zh) * | 2014-11-18 | 2015-10-21 | 深圳大学 | 交流电力网中源荷同变的对称功率传输系数的获取方法 |
WO2016077995A1 (zh) * | 2014-11-18 | 2016-05-26 | 深圳大学 | 交流电力网的发电和用电有功实时调度的对称方法 |
Also Published As
Publication number | Publication date |
---|---|
WO2020154847A1 (zh) | 2020-08-06 |
US11165251B2 (en) | 2021-11-02 |
CN111758197A (zh) | 2020-10-09 |
US20210226446A1 (en) | 2021-07-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Solving OPF using linear approximations: fundamental analysis and numerical demonstration | |
Yang et al. | LMP revisited: A linear model for the loss-embedded LMP | |
Wang et al. | Power system state estimation via feasible point pursuit: Algorithms and Cramér-Rao bound | |
CN111758197B (zh) | 获取交流电力网中源荷的网损功率分量的对称方法 | |
Groß et al. | On the steady-state behavior of a nonlinear power system model | |
Li et al. | Robust data-driven and fully distributed volt/var control for active distribution networks with multiple virtual power plants | |
Wang et al. | Model order reduction for transient simulation of active distribution networks | |
Liu et al. | A network analysis modeling method of the power electronic converter for hardware-in-the-loop application | |
Wang et al. | A two-stage method for assessment of voltage stability in power system with renewable energy | |
Sandoval et al. | Three-phase second-order analytic probabilistic load flow with voltage-dependent load | |
Bhowmik et al. | Power flow analysis of power system using power perturbation method | |
Zhu et al. | Fractional modelling and simulation for single‐phase PWM rectifier | |
Mazumder et al. | Stability analysis of micropower network | |
Liu et al. | Interval power flow calculation algorithm for multi‐terminal dc distribution networks considering distributed generation output uncertainties | |
Zhang | A projection iterative algorithm for the Signorini problem using the boundary element method | |
Gerstner et al. | A domain decomposition approach for solving dynamic optimal power flow problems in parallel with application to the German transmission grid | |
Murthy et al. | Topology based approach for efficient load flow solution of radial distribution networks | |
CN111758196B (zh) | 获取交流电力网中源荷的节点平移电压分量的对称方法 | |
Jiang et al. | Improving performance of the non-iterative direct state estimation method | |
CN111758195B (zh) | 获取交流电力网中源荷的支路均方电流分量的对称方法 | |
Schneider et al. | Dynamic state estimation utilizing high performance computing methods | |
Archundia-Aranda et al. | Harmonic penetration method for radial distribution networks | |
Sun et al. | Covariance analysis of lav robust dynamic state estimation in power systems | |
CN106021191A (zh) | 一种基于分解滤波的网络等值算法 | |
Ebrahimi et al. | Performance verification of parametric average-value model of line-commutated rectifiers under unbalanced conditions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |