CN112600201B - 基于多维全纯嵌入法的高维静态电压稳定边界计算方法 - Google Patents
基于多维全纯嵌入法的高维静态电压稳定边界计算方法 Download PDFInfo
- Publication number
- CN112600201B CN112600201B CN202011437687.7A CN202011437687A CN112600201B CN 112600201 B CN112600201 B CN 112600201B CN 202011437687 A CN202011437687 A CN 202011437687A CN 112600201 B CN112600201 B CN 112600201B
- Authority
- CN
- China
- Prior art keywords
- node
- power
- voltage
- multidimensional
- multivariable
- 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 58
- 230000003068 static effect Effects 0.000 title claims abstract description 51
- 238000004364 calculation method Methods 0.000 title claims abstract description 10
- 238000013507 mapping Methods 0.000 claims abstract description 32
- 230000014509 gene expression Effects 0.000 claims abstract description 18
- 239000000243 solution Substances 0.000 claims description 12
- 238000002347 injection Methods 0.000 claims description 7
- 239000007924 injection Substances 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 6
- 238000002940 Newton-Raphson method Methods 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 150000001875 compounds Chemical class 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- H02J3/06—Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2101/00—Indexing scheme relating to the type of digital function generated
- G06F2101/08—Powers or roots
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Economics (AREA)
- Power Engineering (AREA)
- Public Health (AREA)
- Computing Systems (AREA)
- Operations Research (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Human Resources & Organizations (AREA)
- Marketing (AREA)
- Primary Health Care (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明公开了一种基于多维全纯嵌入法的高维静态电压稳定边界计算方法,属于电力系统分析领域,首先运用含多个物理映射因子的多维全纯嵌入法求解电力系统潮流方程,从而得到节点电压的多变量幂级数形式的解析表达式;然后基于多维柯西‑阿达马定理,推导出高维静态电压稳定边界的解析表达式,从而计算出高维静态电压稳定边界。本发明运用含多个物理映射因子的多维全纯嵌入法并结合多维柯西‑阿达马定理,以解析形式快速准确地计算出高维静态电压稳定边界。
Description
技术领域
本发明属于电力系统分析技术领域,具体涉及一种基于多维全纯嵌入法的高维静态电压稳定边界计算方法。
背景技术
如何快速准确地计算电力系统高维静态电压稳定边界对于评估电力系统电压稳定性具有重要意义。
传统的连续潮流法在各个方向上不断地采用迭代型预测-修正策略来求解高维静态电压稳定边界,但此类方法在面对大型电力系统时会产生巨大的计算负担,难以实现在线应用。后来,提出了如切曲面法,轨道法等技术,但此类技术只用到了高维静态电压稳定边界上有限个点的信息,导致其准确性难以保证。此外,基于加勒金法的参数多项式方法可全局地逼近高维静态电压稳定边界,具有一定的探索价值。但总的来说,上述所有方法在求解电力系统潮流方程时大多采用传统迭代型潮流求解算法,如牛顿-拉夫逊法等,而牛顿-拉夫逊法在一定程度上存在发散问题或产生不可运行解的情况,这会间接影响上述方法结果的可靠性。
作为一种非迭代型潮流求解算法,全纯嵌入法在嵌入形式适当时可以保证求得电力系统的可运行电压解。因此,往往结合全纯嵌入法和帕德近似解析延拓技术来求解某一方向上的静态电压稳定边界,但此类方法无法直接计算出高维静态电压稳定边界且在面临大型电力系统时往往费时。总的来说,现有技术方法难以快速准确地求解电力系统的高维静态电压稳定边界。
发明内容
本发明的目的在于,提供一种基于多维全纯嵌入法的高维静态电压稳定边界计算方法,快速准确地求得电力系统的高维静态电压稳定边界。
本发明提供一种基于多维全纯嵌入法的高维静态电压稳定边界计算方法,其特征在于,包括以下步骤:
S1、基于多维全纯嵌入法在电力系统潮流方程中嵌入多个物理映射因子;
S2、求解节点电压的多变量幂级数形式的解析表达式;
S3、基于多维柯西-阿达马定理,计算出高维静态电压稳定边界。
进一步地,步骤S1具体包括:
S11、写出复数形式的电力系统潮流方程;
S12、基于多维全纯嵌入法,在潮流方程中嵌入多个物理映射因子。
进一步地,步骤S11中复数形式的电力系统潮流方程为:
式中,Yik表示节点导纳矩阵中第i行第k列的元素,表示Yik的共轭,Pi和Si分别表示节点i的有功功率注入和复功率注入,Vi表示节点i的电压,表示Vi的共轭,Vk表示节点i的相邻节点k的电压,表示Vk的共轭,Vsw表示平衡节点的电压,表示PV节点i电压的特定幅值,N表示电力系统的节点总数。
进一步地,步骤S12具体为:
基于多维全纯嵌入法,在潮流方程中嵌入多个物理映射因子s1,s2,…,s2D,得到含多个物理映射因子的潮流方程,含多个物理映射因子的潮流方程如下:
式中,2D表示2D个不同的维度,s1,s2,…,s2D表示2D个全纯物理映射因子,Qi表示节点i的无功功率注入量,Vi(s1,s2,…,s2D)和Qi(s1,s2,…,s2D)分别表示多变量幂级数形式的节点i处的电压和无功功率全纯函数:
式中,Vi[n1,n2,…,n2D]和Qi[n1,n2,…,n2D]分别表示节点i的电压多变量幂级数和无功功率多变量幂级数的系数。
进一步地,步骤S2具体包括:
S21、计算含多个物理映射因子的潮流方程的初始解;
S22、建立节点电压多变量幂级数系数之间的递归关系;
S23、求解节点电压多变量幂级数形式的解析表达式。
进一步地,步骤S21具体为:
将s1=0,s2=0,…,s2D=0代入含多个物理映射因子的潮流方程中:
计算出节点电压多变量幂级数系数的初始解,即第0阶系数Vi[0,0,…,0]。
进一步地,步骤S22具体为:
使含多个物理映射因子的潮流方程两侧得同阶s1s2…s2D的系数一一对应相等,建立节点电压多变量幂级数的当前阶系数和之前阶系数的递归关系。
进一步地,步骤S23具体为:
根据节点电压多变量幂级数的当前阶系数和之前阶系数的递归关系,代入节点电压多变量幂级数系数的初始解,通过递归关系求得节点电压多变量幂级数的所有系数,进而得到节点电压多变量幂级数形式的解析表达式:
进一步地,步骤S3具体包括为:
S31、运用多维柯西-阿达马定理处理节点电压多变量幂级数表达式;
S32、计算出电力系统高维静态电压稳定边界。
进一步地,步骤S3具体为:
由多维柯西-阿达马定理可知,多变量幂级数形式的全纯函数的收敛域{sr1,sr2,…,sr2D}可由下式确定:
式中,Vi[n1,n2…,n2D]表示节点电压多变量幂级数的系数,sr1,sr2,…,sr2D表示节点电压多变量幂级数在各个方向上的收敛半径;
基于多维柯西-阿达马定理,推导出高维静态电压稳定边界与节点电压多变量幂级数的前后项比例之间的解析关系:
本发明的有益效果是:本发明的基于多维全纯嵌入法的高维静态电压稳定边界计算方法采用含多个物理映射因子的多维全纯嵌入法求解潮流方程,可得节点电压的多变量幂级数形式的解析表达式;根据多维柯西-阿达马定理推导出的高维静态电压稳定边界与节点电压多变量幂级数的前后项比例之间的解析关系,能够快速准确地计算电力系统的高维静态电压稳定边界。
附图说明
图1是本发明实施例提供的一种基于多维全纯嵌入法的高维静态电压稳定边界计算方法的流程示意图;
图2是本发明实施例提供的一种表示含m个PQ(负荷)节点,p个PV(发电机)节点和1个平衡节点的电力系统的示意图;
图3是本发明实施例提供的一种IEEE 39节点标准电力系统的网络结构示意图;
图4是本发明实施例提供的一种IEEE 39节点标准电力系统的区域划分示意图;
图5是本发明实施例提供的通过多维全纯嵌入法求得的IEEE 39节点标准电力系统的节点7电压示意图;
图6是本发明实施例提供的一种IEEE 39节点标准电力系统的高维静态电压稳定边界示意图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。此外,下面所描述的本发明各个实施方式中所涉及到的技术特征只要彼此之间未构成冲突就可以相互组合。
本发明的目的是解决现有方法所存在的难以快速准确计算电力系统高维静态电压稳定边界的技术问题:通过在电力系统潮流方程中嵌入多个物理映射因子,求得节点电压的多变量幂级数形式的解析表达式,之后根据多维柯西-阿达马定理推导出高维静态电压稳定边界与节点电压多变量幂级数的前后项比例之间的解析关系,从而快速准确地计算电力系统高维静态电压稳定边界。
图1是本发明实施例提供的一种基于多维全纯嵌入法的高维静态电压稳定边界计算方法的流程示意图,包括以下步骤:
S1、基于多维全纯嵌入法在电力系统潮流方程中嵌入多个物理映射因子;
在本发明实施例中,步骤S1可以通过以下方式实现:
图2表示含m个负荷PQ节点,p个发电机PV节点和1个平衡节点的电力系统,其潮流方程为:
其中,Yik是节点导纳矩阵中第i行第k列的元素,表示Yik的共轭,Pi和Si表示节点i的有功功率注入和复功率注入,Vi是节点i的电压,表示Vi的共轭,Vk为节点i的相邻节点k的电压,表示Vk的共轭,Vsw为平衡节点的电压,代表PV(发电机)节点i电压的特定幅值,N为电力系统的节点总数。
接着,在公式(1)中嵌入多个物理映射因子s1,s2,…,s2D,得到含多个物理映射因子的电力系统潮流方程:
其中,2D表示2D个不同的维度,s1,s2,…,s2D表示2D个全纯物理映射因子,可分别用来独立调整电力系统不同节点的负荷或发电功率,Yik是节点导纳矩阵中第i行第k列的元素,Pi和Si表示节点i的有功功率注入和复功率注入,Qi表示节点i的无功功率注入量,Vi是节点i的电压,表示Vi的共轭,Vk为节点i的相邻节点k的电压,表示Vk的共轭,Vsw为平衡节点的电压,代表PV(发电机)节点i电压的特定幅值,N为电力系统的节点总数,Vi(s1,s2,…,s2D)和Qi(s1,s2,…,s2D)分别为多变量幂级数形式的节点i处的电压和无功功率全纯函数,即:
其中,s1,s2,…,s2D表示2D个全纯物理映射因子,可分别用来独立调整电力系统不同节点的负荷或发电功率,Vi[n1,n2,…,n2D]和Qi[n1,n2,…,n2D]分别为节点i的电压多变量幂级数和无功功率多变量幂级数的系数。
经过步骤S1,可得到含多个物理映射因子的电力系统潮流方程。
S2、求解节点电压的多变量幂级数形式的解析表达式;
在本发明实施例中,步骤S2可以通过以下方式实现:
将s1=0,s2=0,…,s2D=0代入含多个物理映射因子的潮流方程(2)中:
根据公式(4)计算出节点电压多变量幂级数系数的初始解,即第0阶系数Vi[0,0,…,0]。
然后,将含多个物理映射因子的潮流方程两侧同阶s1s2…s2D的系数一一对应相等,可建立节点电压多变量幂级数的当前阶系数(未知量)和之前阶系数(已知量)的递归关系。根据节点电压多变量幂级数的当前阶系数和之前阶系数的递归关系,代入节点电压多变量幂级数系数的初始解(即第0阶系数),即可通过递归关系求得节点电压多变量幂级数的所有系数,从而得到节点电压多变量幂级数形式的解析表达式,即:
经过步骤S2,可得到电力系统节点电压多变量幂级数形式的解析表达式。
S3、基于多维柯西-阿达马定理,快速准确地计算出高维静态电压稳定边界;
在本发明实施例中,步骤S3可以通过以下方式实现:
由多维柯西-阿达马定理可知,多变量幂级数形式的全纯函数的收敛域{sr1,sr2,…,sr2D}可由下式确定:
其中,Vi[n1,n2…,n2D]为节点电压多变量幂级数的系数,sr1,sr2,…,sr2D为节点电压多变量幂级数在各个方向上的收敛半径。
接下来推导高维静态电压稳定边界与节点电压多变量幂级数的前后项比例之间的解析关系。
首先,定义节点电压多变量幂级数前后项比例的上界为L:
接着,给定一个任意小的正数ξ>0,则存在n使得下式满足:
那么类似公式(8),从第m阶项到第n阶项满足公式(9):
其中,m<n。将公式(9)相乘可得:
整理可得:
由于ξ是一个任意小的正数,于是对公式(11)两侧取极限并结合公式(7)可得:
类似公式(7)至公式(12)的推导过程,对于节点电压多变量幂级数前后项比例的下限,可得:
其中,inf表示当n趋于无穷大时的下界。
结合公式(12)和公式(13)可得:
因此,如果节点电压多变量幂级数前后项比例的极限存在,即公式(15)满足:
那么由公式(14)和公式(15)可得:
结合公式(6),可得:
经过步骤S3,可以快速准确地计算出电力系统高维静态电压稳定边界。
下面通过在图3所示的IEEE 39节点标准电力系统中的实施例,并结合附图,对本发明的技术方案作进一步具体的说明。
首先,将图3所示的IEEE 39节点标准电力系统按照几何位置划分为3个区域,如图4所示。然后,根据上述步骤S1中介绍的多维全纯嵌入法,在电力系统潮流方程中嵌入3个物理映射因子s1,s2,s3,分别用于调整区域1、区域2和区域3内的电力系统运行状态,即负荷和发电功率。
接着,根据上述步骤S2中介绍的递归关系求解节点电压的多变量幂级数形式的解析表达式其中,s1,s2,s3为3个独立的物理映射因子,分别用于调整区域1、区域2和区域3内的电力系统运行状态,即负荷和发电功率,Vi[n1,n2,n3]表示节点电压多变量幂级数系数。由步骤2所得的节点7的电压如图5所示,其中s3=1。
最后,基于上述步骤S3中引用的多维柯西-阿达马定理,推导出电力系统高维静态电压稳定边界的解析表达式其中,和分别表示节点电压多变量幂级数的第n阶和第(n+1)阶系数,sr1,sr2,sr3为节点电压多变量幂级数在各个方向上的收敛半径,共同构成电力系统高维静态电压稳定边界。由步骤3所得高维静态电压稳定边界如图6所示,其中s3=1,静态电压稳定的参考边界取自图5。由图6可知,随着节点电压多变量幂级数的阶数上升,所计算出的静态电压稳定边界越来越趋近参考边界,且在n=9时就已非常切合参考边界了,即节点电压多变量幂级数的10次项比上9次项。此实施例证明了通过所提出的方法,仅需要节点电压多变量幂级数的前面少数几项,即可快速准确地计算出电力系统的高维静态电压稳定边界。
需要指出,根据实施的需要,可将本申请中描述的各个步骤/部件拆分为更多步骤/部件,也可将两个或多个步骤/部件或者步骤/部件的部分操作组合成新的步骤/部件,以实现本发明的目的。
尽管本发明较多地使用了电力系统、多维全纯嵌入法、潮流方程、多变量幂级数、高维静态电压稳定边界、解析表达式等术语,但并不排除使用其它术语的可能性。使用这些术语仅仅是为了更方便地描述和解释本发明的本质;把它们解释成任何一种附加的限制都是与本发明精神相违背的。
本领域的技术人员容易理解,以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种基于多维全纯嵌入法的高维静态电压稳定边界计算方法,其特征在于,包括以下步骤:
S1、基于多维全纯嵌入法在电力系统潮流方程中嵌入多个物理映射因子;
S2、求得节点电压的多变量幂级数形式的解析表达式;具体包括:
S21、计算含多个物理映射因子的潮流方程的初始解;
S22、建立节点电压多变量幂级数系数之间的递归关系;
S23、求得节点电压多变量幂级数形式的解析表达式;
S3、基于多维柯西-阿达马定理,计算出高维静态电压稳定边界;具体包括:
S31、运用多维柯西-阿达马定理处理节点电压多变量幂级数表达式;
S32、进而计算出电力系统高维静态电压稳定边界。
2.根据权利要求1所述的基于多维全纯嵌入法的高维静态电压稳定边界计算方法,其特征在于,步骤S1具体包括:
S11、写出复数形式的电力系统潮流方程;
S12、基于多维全纯嵌入法,在潮流方程中嵌入多个物理映射因子。
4.根据权利要求3所述的基于多维全纯嵌入法的高维静态电压稳定边界计算方法,其特征在于,步骤S12具体为:
基于多维全纯嵌入法,在潮流方程中嵌入多个物理映射因子s1,s2,…,s2D,得到含多个物理映射因子的潮流方程,含多个物理映射因子的潮流方程如下:
式中,2D表示2D个不同的维度,s1,s2,…,s2D表示2D个全纯物理映射因子,Qi表示节点i的无功功率注入量,Vi(s1,s2,…,s2D)和Qi(s1,s2,…,s2D)分别表示多变量幂级数形式的节点i处的电压和无功功率全纯函数:
式中,Vi[n1,n2,…,n2D]和Qi[n1,n2,…,n2D]分别表示节点i的电压多变量幂级数和无功功率多变量幂级数的系数。
6.根据权利要求5所述的基于多维全纯嵌入法的高维静态电压稳定边界计算方法,其特征在于,步骤S22具体为:
使含多个物理映射因子的潮流方程两侧的同阶s1 s2…s2D的系数一一对应相等,建立节点电压多变量幂级数的当前阶系数和之前阶系数的递归关系。
8.根据权利要求7所述的基于多维全纯嵌入法的高维静态电压稳定边界计算方法,其特征在于,步骤S3具体为:
由多维柯西-阿达马定理可知,多变量幂级数形式的全纯函数的收敛域{sr1,sr2,…,sr2D}可由下式确定:
式中,sup表示当n趋于无穷大时的上界,Vi[n1,n2…,n2D]表示节点电压多变量幂级数的系数,sr1,sr2,…,sr2D表示节点电压多变量幂级数在各个方向上的收敛半径;
基于多维柯西-阿达马定理,推导出高维静态电压稳定边界与节点电压多变量幂级数的前后项比例之间的解析关系:
式中,sr1,sr2,…,sr2D为节点电压多变量幂级数在各个方向上的收敛半径,共同构成电力系统高维静态电压稳定边界。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011437687.7A CN112600201B (zh) | 2020-12-07 | 2020-12-07 | 基于多维全纯嵌入法的高维静态电压稳定边界计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011437687.7A CN112600201B (zh) | 2020-12-07 | 2020-12-07 | 基于多维全纯嵌入法的高维静态电压稳定边界计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112600201A CN112600201A (zh) | 2021-04-02 |
CN112600201B true CN112600201B (zh) | 2022-06-03 |
Family
ID=75191596
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011437687.7A Active CN112600201B (zh) | 2020-12-07 | 2020-12-07 | 基于多维全纯嵌入法的高维静态电压稳定边界计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112600201B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113794204B (zh) * | 2021-08-19 | 2024-03-26 | 西安交通大学 | 电压稳定性分析方法、装置、电子设备及存储介质 |
CN114336635B (zh) * | 2022-01-05 | 2023-07-25 | 国网内蒙古东部电力有限公司通辽供电公司 | 基于常项值和先验节点的全纯嵌入潮流计算方法、装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107133406A (zh) * | 2017-05-10 | 2017-09-05 | 东北电力大学 | 一种电力系统静态电压稳定域边界的快速搜索方法 |
CN109167366A (zh) * | 2018-08-30 | 2019-01-08 | 国网山西省电力公司电力科学研究院 | 一种电力系统静态电压稳定临界点的计算方法 |
WO2019033820A1 (zh) * | 2017-08-17 | 2019-02-21 | 东北电力大学 | 一种电力系统静态电压稳定域边界快速搜索的优化模型 |
CN110957728A (zh) * | 2018-09-27 | 2020-04-03 | 国网甘肃省电力公司电力科学研究院 | 一种三相有源配电网的全纯嵌入潮流方法及装置 |
CN111682545A (zh) * | 2020-05-27 | 2020-09-18 | 杭州电子科技大学 | 基于helm的配电网二阶电压灵敏度计算方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US11271398B2 (en) * | 2018-07-13 | 2022-03-08 | University Of Tennessee Research Foundation | Voltage stability assessment, control and probabilistic power flow based on multi-dimensional holomorphic embedding techniques |
-
2020
- 2020-12-07 CN CN202011437687.7A patent/CN112600201B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107133406A (zh) * | 2017-05-10 | 2017-09-05 | 东北电力大学 | 一种电力系统静态电压稳定域边界的快速搜索方法 |
WO2019033820A1 (zh) * | 2017-08-17 | 2019-02-21 | 东北电力大学 | 一种电力系统静态电压稳定域边界快速搜索的优化模型 |
CN109167366A (zh) * | 2018-08-30 | 2019-01-08 | 国网山西省电力公司电力科学研究院 | 一种电力系统静态电压稳定临界点的计算方法 |
CN110957728A (zh) * | 2018-09-27 | 2020-04-03 | 国网甘肃省电力公司电力科学研究院 | 一种三相有源配电网的全纯嵌入潮流方法及装置 |
CN111682545A (zh) * | 2020-05-27 | 2020-09-18 | 杭州电子科技大学 | 基于helm的配电网二阶电压灵敏度计算方法 |
Non-Patent Citations (2)
Title |
---|
A Multi-Dimensional Holomorphic Embedding Method to Solve AC Power Flows;Chengxi Liu等;《IEEE Access》;20171109;全文 * |
基于HELM网损灵敏度的配电网无功优化研究;赵亦岚等;《浙江电力》;20200922;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN112600201A (zh) | 2021-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN112600201B (zh) | 基于多维全纯嵌入法的高维静态电压稳定边界计算方法 | |
CN112636358B (zh) | 一种基于多变量商差法的电力系统潮流计算方法 | |
CN109167366B (zh) | 一种电力系统静态电压稳定临界点的计算方法 | |
CN112821392B (zh) | 一种基于全纯嵌入法的静态电压稳定边界计算方法 | |
CN116706921B (zh) | 基于hhl算法的量子牛顿-拉夫逊法潮流计算方法和系统 | |
CN106532712B (zh) | 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 | |
CN111049144B (zh) | 基于定变化率矩阵的有源配电网潮流计算pv节点处理方法 | |
CN110148946B (zh) | 一种基于辅助因子两步求解的孤岛微电网潮流计算方法 | |
Lewis et al. | A comparison of some particle-in-cell plasma simulation methods | |
CN108649585B (zh) | 一种电力系统静态电压稳定域边界快速搜索的直接法 | |
CN114566969A (zh) | 一种适合研究目的使用的直角坐标牛顿法潮流计算方法 | |
CN111181166B (zh) | 一种预测校正的不确定性仿射潮流方法 | |
CN111639463B (zh) | 一种基于XGBoost算法的电力系统扰动后频率特征预测方法 | |
CN110649624B (zh) | 一种电力系统潮流并行计算方法 | |
Kumar et al. | A simplified approach for load flow analysis of radial distribution network | |
Ndanusa et al. | Predictor–corrector methods of high order for numerical integration of initial value problems | |
Veseth et al. | Polarizabilities and photoionization cross sections of OH and HF | |
CN112072634B (zh) | 一种基于潮流嵌入技术的潮流计算方法 | |
CN110717145A (zh) | 基于对称稀疏矩阵技术的分段对称反向高斯-约当消元法 | |
Heinrichs | An adaptive least-squares spectral collocation method with triangular elements for the incompressible Navier–Stokes equations | |
CN112383063B (zh) | 一种基于最小二乘法拟合修正的输配协同潮流计算方法 | |
CN112383064B (zh) | 一种基于迭代变量历史值更新的输配协同潮流计算方法 | |
Sambariya et al. | Stable reduction methods of linear dynamic systems in frequency domain | |
CN111030122B (zh) | 首次迭代雅可比矩阵改变的直角坐标牛顿法潮流计算方法 | |
KR20030071938A (ko) | 데이터 마이닝을 위한 최적의 의사 결정 나무 선택 장치및 그 방법 |
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 |