CN110518590A - 考虑负荷电压静特性的配电网线性潮流计算方法 - Google Patents
考虑负荷电压静特性的配电网线性潮流计算方法 Download PDFInfo
- Publication number
- CN110518590A CN110518590A CN201910719366.7A CN201910719366A CN110518590A CN 110518590 A CN110518590 A CN 110518590A CN 201910719366 A CN201910719366 A CN 201910719366A CN 110518590 A CN110518590 A CN 110518590A
- Authority
- CN
- China
- Prior art keywords
- node
- power
- voltage
- load
- distribution network
- 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.)
- Granted
Links
- 238000004364 calculation method Methods 0.000 title claims abstract description 43
- 230000003068 static effect Effects 0.000 title claims abstract description 29
- 238000000034 method Methods 0.000 claims abstract description 36
- 230000009466 transformation Effects 0.000 claims abstract description 7
- 239000011159 matrix material Substances 0.000 claims description 48
- 239000013256 coordination polymer Substances 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 9
- 238000002347 injection Methods 0.000 claims description 6
- 239000007924 injection Substances 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 5
- 238000003672 processing method Methods 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 2
- 238000011160 research Methods 0.000 claims description 2
- 230000005611 electricity Effects 0.000 claims 1
- 238000004445 quantitative analysis Methods 0.000 abstract description 2
- 238000012360 testing method Methods 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000002940 Newton-Raphson method Methods 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000000844 transformation Methods 0.000 description 2
- 239000007983 Tris buffer Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005206 flow analysis Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000012804 iterative process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004088 simulation 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/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
-
- 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)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
考虑负荷电压静特性的配电网线性潮流计算方法,包括:结合配电网的典型特征,定量分析适用于配电网的简化条件;根据电压的变化将负荷模型划分为恒阻抗分量、恒电流分量和恒功率分量,建立一种考虑电压静特性的ZIP负荷模型,并根据配电网的简化条件,对其线性化处理;将潮流模型存在的非线性项,通过最大和最小电压加权因子进行加权,得到一个线性的广义函数,并结合配电网的简化条件,将潮流模型简化成一种全线性化潮流模型;对全线性化潮流模型进行初等变换,实现对电压幅值和相位角的解耦。根据配电网各节点的已知量,求出网络中各个节点的电压和相角。本发明一种考虑负荷电压静特性的配电网线性潮流计算方法,与现有的线性潮流相比,该方法具有较高的精度和通用性。
Description
技术领域
本发明涉及配电网重构技术领域,具体涉及一种考虑负荷电压静特性的配电网线性潮流计算方法。
背景技术
潮流计算是现代配电系统分析最为基本的部分,是配电系统监测、控制和决策最重要的工具。目前,常用的配电网潮流计算方法主要包括前推回代法、改进牛顿法、隐式Zbus高斯法和回路阻抗法等。但由于这些方法的潮流方程都是非线性的,涉及到迭代过程,因此计算的效率不高,难以满足配电网实时运行和快速分析的需求。
为了提高潮流计算方法的效率,已有许多文献展开了相关的研究,目前典型的处理方法是根据配电系统树状网络的结构,以回路分析法作为潮流计算的理论,利用配电网的典型特征,对回路电压方程进行线性化处理,使其简化为线性代数方程组,最终通过一次公式代入便可求出电压分布结果。或是对基本潮流方程中的三角函数项进行多项式拟合,再利用系统的运行特点将电压幅值和相角解耦,然后对非线性项进行不同程度的泰勒级数展开,最终得到一种全线性的潮流方程。
上述方法虽然有效的减少了潮流计算的时间,但只考虑了恒功率负荷,且相比于迭代法精度较低。然而,在实际潮流计算时,更加通用的是采用具有电压静特性的ZIP负荷模型,要求保持高精度的电压幅值。
针对以上存在的问题,采用简单曲线拟合技术,将电压相关的负荷模型分解为恒阻抗分量和恒电流分量,建立了基于拟合系数的ZI负荷模型,可将负荷模型应用于线性潮流计算中;或是结合ZIP负荷模型的特点,将恒阻抗分量和恒电流分量对应的负荷折算到标称电压下,采用三个特定的矩阵对负荷模型进行统一处理,在进行潮流计算时考虑了负荷模型对电压的依赖。
以上改进方法虽然考虑了负荷的电压静特性,但增加了计算的负担,难以解决环网问题,且在处理病态系统时,计算精度较差。
发明内容
针对现有的线性潮流方法仅考虑恒功率负荷和计算精度低的缺点,本发明提供一种考虑负荷电压静特性的配电网线性潮流计算方法。仿真结果表明,与现有的线性潮流相比,该方法具有较高的精度和通用性。
本发明采取的技术方案为:
在考虑电压静特性的ZIP负荷模型的基础上,利用配电网的典型特征,对潮流方程和负荷模型的非线性项进行线性化处理,以降低潮流计算的复杂程度。同时,提出了一种通过最大和最小电压加权因子的广义公式来求解潮流方程,实现对潮流模型的高精度逼近。并对简化后的潮流方程进行初等变换,实现了电压幅值和相角的二次解耦。
考虑负荷电压静特性的配电网线性潮流计算方法,包括以下步骤:
步骤1:结合配电网的典型特征,定量分析适用于配电网的简化条件;
步骤2:根据电压的变化将负荷模型划分为恒阻抗分量、恒电流分量和恒功率分量,建立一种考虑电压静特性的ZIP负荷模型,并根据配电网的简化条件,对其线性化处理;
步骤3:将潮流模型存在的非线性项1/Vi,通过最大和最小电压加权因子进行加权,得到一个线性的广义函数,并结合配电网的简化条件,将潮流模型简化成一种全线性化潮流模型;
步骤4:对全线性化潮流模型进行初等变换,实现对电压幅值和相位角的解耦。
步骤5:根据配电网各节点的已知量,求出网络中各个节点的电压和相角。
所述步骤1中,配电网的典型特征包括:①、节点电压幅值趋近于1.0p.u.;②、线路两端的相角非常小,使得配电网中所有节点的电压相角与平衡节点的电压相角相差不大;③、R/X比值比较大,一般接近或者大于1。定量分析适用于配电网的简化条件,可以降低潮流方程的复杂程度,提高潮流计算的效率。
所述步骤1中,配电网的简化条件,如公式(1)所示:
所述步骤2中,负荷模型是指当负荷的实际电压偏离额定电压时,负荷的特性也将发生变化。因此,根据电压的变化可将负荷模型划分为恒阻抗分量、恒电流分量和恒功率分量。其中,恒阻抗分量和恒电流分量分别指负荷模型的电压二次项和电压一次项,恒功率分量是指负荷模型的常数项。
所述步骤2中,根据电压的变化,建立一种考虑电压静特性的ZIP负荷模型,如公式(5)所示:
其中,P(V),Q(V)分别表示节点负载有功功率和无功功率。V,VN分别为节点实际电压和额定电压;PN,QN分别表示额定电压下的有功功率和无功功率。CZ,CI,CP分别表示节点有功恒阻抗、恒电流、恒功率负荷的比例系数;C'Z,C'I,C'P分别表示节点无功恒阻抗、恒电流、恒功率负荷的比例系数;各参数满足的约束条件为CZ+CI+CP=1,C'Z+C'I+C'P=1。在实际潮流计算时,采用更加通用的ZIP负荷模型可以更加精确的反映系统的电压分布。
所述步骤2中,ZIP负荷模型线性化处理,以有功功率为例,如公式(2)所示:
其中,P(V)表示节点负载有功功率;PN分别表示额定电压下的有功功率;CZ,CI,CP分别表示节点有功恒阻抗、恒电流、恒功率负荷的比例系数,V为节点实际电压,△V为实际电压与额定电压的差值,通常取值为0~0.1p.u.,其二阶及高阶项的取值相对较小,潮流计算中可忽略不计。PP,PI分别为ZIP负荷模型线性化后有功功率的恒功率负荷和恒电流负荷新的比例系数。对ZIP负荷模型进行线性化处理,可以降低模型的复杂程度,减小潮流计算的负担,提高计算效率。
所述步骤3中,根据有功功率线性化处理的方法,无功功率线性表达式如公式(3)所示:
Q(V)≈QP+QIV (3)
其中,Q(V)表示节点负载有功功率;QP,QI分别为ZIP负荷模型线性化后无功功率的恒功率负荷和恒电流负荷新的比例系数。
所述步骤3中,潮流方程表达式如公式(6)所示:
所述步骤3中,线性的广义函数K如公式(4)所示:
其中,i为节点的编号;Vi为节点i的电压幅值;Vmin,Vmax分别为电压最小加权因子和电压最大加权因子,要求Vmin≤Vi≤Vmax。
所述步骤3中,全线性潮流方程如公式(7)所示:
其中,Vmin,Vmax分别为电压最小加权因子和电压最大加权因子;S,W,R分别为配电网中的平衡节点、PQ节点、PV节点的集合;PPS,PPW,QPW分别平衡节点有功功率、PQ节点有功功率、PQ节点无功功率的恒功率系数;PIS,PIW,QIW分别为平衡节点有功功率、PQ节点有功功率、PQ节点无功功率的恒电流系数;VS、δS分别为平衡节点电压幅值和相角;VW、δW分别为PQ节点电压幅值和相角;VR、δR分别为PV节点电压幅值和相角;diag(PPR)、diag(PPW)、diag(QPW)分别为PV节点有功功率的恒功率负荷系数对角矩阵、PQ节点有功功率的恒功率负荷系数对角矩阵、PV节点无功功率的恒功率负荷系数对角矩阵;GRS和BRS分别为PV节点与平衡节点的互导纳矩阵的实部和虚部;GWS和BWS分别为PQ节点与平衡节点的互导纳矩阵的实部和虚部;GWR和BWR分别为PQ节点与PV节点的互导纳矩阵的实部和虚部;GRR和BRR为PV节点自导纳矩阵的实部和虚部;BRW和GRW分别为PQ节点与PV节点的互导纳矩阵的实部和虚部;GWW和BWW分别为PQ节点自导纳矩阵的实部和虚部。
所述步骤4中,为了实现电压幅值和相角的第二次解耦,可将公式(7)进行两次初等变换如式(8)所示:
其中,为已知项,可以通过节点导纳矩阵和已知节点功率求出;为已知项,可通过以平衡节点、PQ节点、PV节点所组成的导纳矩阵求出;为系统PQ节点和PV节点的相角集合,为PQ节点的电压集合,两者均为未知量。
所述步骤5中,计算网络中各个节点的电压和相角模型,如公式(9)所示:
JX=B (9)
其中,B为(n+m-2)×1阶矩阵,其元素可由节点注入功率、导纳矩阵、平衡节点的电压幅值和PV节点的电压幅值求出,J为(n+m-2)×(n+m-2)阶方阵,其元素可由节点注入功率和导纳矩阵求得。X为系统未知电压和相角的集合。该模型的个数总和为m+n-2刚好等于未知量的个数,因此,可求出网络中各个节点的电压和相角。
本发明一种考虑负荷电压静特性的配电网线性潮流计算方法,技术效果如下:
1)、通用性强:针对现有配电网线性潮流采用恒功率负荷模型,不能准确反映系统潮流电压的分布,本发明在进行潮流。计算时考虑了负荷模型对电压的依赖,建立了更加通用的ZIP负荷模型。
2)、计算效率高:本发明定量分析适用于配电网的典型特征,对ZIP负荷模型和潮流方程的非线性项进行线性化处理,最终得到一种全线性的潮流计算模型。降低了潮流模型的复杂程度,提高了潮流计算的效率,可用于配电系统的快速潮流分析和实时调度。
3)、精度高:本发明提出一种通过最大和最小电压加权因子的广义公式来求解潮流方程,能够在不需要迭代的情况下实现对潮流模型的高精度逼近。
4)、抗干扰能力强:本发明对高R/X比值、过负载系统和弱环网都具有较强的适应性。高R/X比值、弱环网系统、过负载等多种场景下,本发明仍具有较高的精度,抗干扰能力强。
附图说明
图1为IEEE33节点系统结构示意图。
图2为恒功率负荷模型电压幅值分布图。
图3为恒功率负荷模型电压相角分布图。
图4为不同比例系数的ZIP负荷模型电压幅值分布图。
图5为ZIP负荷模型电压幅值和相角误差分布图。
具体实施方式
步骤1:结合配电网的典型特征,定量分析适用于配电网的简化条件;
步骤2:根据电压的变化将负荷模型划分为恒阻抗分量、恒电流分量和恒功率分量,建立一种考虑电压静特性的ZIP负荷模型,并根据配电网的简化条件,对其线性化处理;
步骤3:将潮流模型存在的非线性项1/Vi,通过最大和最小电压加权因子进行加权,得到一个线性的广义函数,并结合配电网的简化条件,将潮流模型简化成一种全线性化潮流模型;
步骤4:对全线性化潮流模型进行初等变换,实现对电压幅值和相位角的解耦。
步骤5:根据配电网各节点的已知量,求出网络中各个节点的电压和相角。
下面结合附图,对优选实例进行详细说明:
2、适用于配电网简化条件如下:
1)由于两个相邻节点之间的距离非常短,所以沿分支的电压降很小。因此,节点i到节点j相角差δij的值可以足够接近于零,则可以对潮流方程的非线性项sinδij和Vjδij进行简化,如公式(1)所示:
其中,i和j表示节点的编号;δi,δj分别为节点i,j的相角;δij为节点i,j的相角差,Vj为节点j的电压幅值。
2)基于电压静特性的ZIP负荷模型存在电压的平方项Vi 2。则定义Vi=1-△Vi,通过忽略二次项△Vi 2得到其近似线性化,如公式(2)所示:
其中,P(V)表示节点负载有功功率,PN表示额定电压下的有功功率;CZ,CI,CP分别表示节点有功恒阻抗、恒电流、恒功率负荷的比例系数,V为节点实际电压,△V为实际电压与额定电压的差值,通常取值为0~0.1p.u.,其二阶及高阶项的取值相对较小,潮流计算中可忽略不计。PP,PI分别为ZIP负荷模型线性化后恒功率负荷和恒电流负荷新的比例系数。对ZIP负荷模型进行线性化处理,可以降低模型的复杂程度,减小潮流计算的负担,提高计算效率。
3)根据有功功率线性化处理的方法,无功功率线性表达式如公式(3)所示:
Q(V)≈QP+QIV (3)
其中,Q(V)表示节点负载有功功率;QP,QI分别为ZIP负荷模型线性化后无功功率的恒功率
4)令双曲函数K=1/Vi,并将其通过最大和最小电压加权因子进行加权,得到一个线性的广义函数如公式(4)所示:
其中,Vi为节点i的电压幅值;Vmin,Vmax分别为电压最小加权因子和电压最大加权因子,要求Vmin≤Vi≤Vmax。
负荷和恒电流负荷新的比例系数。
2、所述建立一种考虑电压静特性的ZIP负荷模型,如公式(5)所示:
其中,P(V),Q(V)分别表示节点负载有功功率和无功功率。V,VN分别为节点实际电压和额定电压;PN,QN分别表示额定电压下的有功功率和无功功率。CZ(C'Z),CI(C'I),CP(C'P)分别表示节点有功(无功)恒阻抗、恒电流、恒功率负荷的比例系数。各参数满足的约束条件为CZ+CI+CP=1,C'Z+C'I+C'P=1。在实际潮流计算时,采用更加通用的ZIP负荷模型可以更加精确的反映系统的电压分布。
3、所述潮流方程表达式如公式(6)所示:
式中:i和j表示节点的编号;Vi和Vj分别为节点i和j的电压幅值;Pi和Qi分别表示配电网中节点i的总注入有功功率和无功功率;PLi和QLi分别表示节点i负载有功功率和无功功率。PGi和QGi表示DG的有功功率和无功功率;Gij和Bij分别表示导纳矩阵的实部和虚部;δij为节点i,j的相角差。
4、简化后的全线性潮流方程如公式(7)所示:
其中,Vmin,Vmax分别为电压最小加权因子和电压最大加权因子;S,W,R分别为配电网中的平衡节点、PQ节点、PV节点的集合;PPS,PPW,QPW分别平衡节点有功功率、PQ节点有功功率、PQ节点无功功率的恒功率系数;PIS,PIW,QIW分别为平衡节点有功功率、PQ节点有功功率、PQ节点无功功率的恒电流系数;VS、δS分别为平衡节点电压幅值和相角;VW、δW分别为PQ节点电压幅值和相角;VR、δR分别为PV节点电压幅值和相角;diag(PPR)、diag(PPW)、diag(QPW)分别为PV节点有功功率的恒功率负荷系数对角矩阵、PQ节点有功功率的恒功率负荷系数对角矩阵、PV节点无功功率的恒功率负荷系数对角矩阵;GRS和BRS分别为PV节点与平衡节点的互导纳矩阵的实部和虚部;GWS和BWS分别为PQ节点与平衡节点的互导纳矩阵的实部和虚部;GWR和BWR分别为PQ节点与PV节点的互导纳矩阵的实部和虚部;GRR和BRR为PV节点自导纳矩阵的实部和虚部;BRW和GRW分别为PQ节点与PV节点的互导纳矩阵的实部和虚部;GWW和BWW分别为PQ节点自导纳矩阵的实部和虚部。
5、进行两次初等变换如式(8)所示:
其中,为已知项,可以通过节点导纳矩阵和已知节点功率求出;为已知项,可通过以平衡节点、PQ节点、PV节点所组成的导纳矩阵求出;为系统PQ节点和PV节点的相角集合,为PQ节点的电压集合,两者均为未知量。
6、计算网络中各个节点的电压和相角模型,如公式(9)所示:
JX=B (9)
其中,B为(n+m-2)×1阶矩阵,其元素可由节点注入功率、导纳矩阵、平衡节点的电压幅值和PV节点的电压幅值求出,J为(n+m-2)×(n+m-2)阶方阵,其元素可由节点注入功率和导纳矩阵求得。X为系统未知电压和相角的集合。该模型的个数总和为m+n-2刚好等于未知量的个数,因此,可求出网络中各个节点的电压和相角。
7、以牛顿拉夫逊法的计算结果作为基准值,引入电压幅值的最大误差和平均误差作为指标,计算方法如公式(10),(11)所示:
△Vmax=max(|Vi LPF-Vi AC|) (10)
其中,n表示节点的个数;i表示节点编号;△Vmax,△Vmean分别为电压幅值的最大误差和平均误差;Vi LPF表示线性潮流的电压幅值,Vi AC表示NR法的电压幅值;max表示最大值。下面将通过实施例,进一步说明本发明技术效果:
以IEEE33节点系统作为案例,对配电系统简化条件和本发明有效性进行分析。IEEE33节点系统如图1所示,该系统基准电压和基准功率分别为12.66kV,10MVA。假设首端为平衡节点,电压幅值为1.0p.u.,电压相角为0,最小电压加权因子Vmin=0.9p.u.,最大电压加权因子Vmax=1.0p.u.。
图2和图3分别为2种潮流计算方法在IEEE33节点系统的电压幅值和电压相角计算结果分布图。可以发现,本发明与牛顿拉夫逊法的计算结果高度一致,证明了本发明简化条件的合理性。为了测试不同比例系数的ZIP负荷模型对电压幅值的影响,以IEEE33节点系统作为算例,本发明采用的ZIP负荷模型的比例系数分别为:A恒功率负荷(Z:I:P=0:0:1);B恒阻抗负荷(Z:I:P=1:0:0);C恒电流负荷(Z:I:P=0:1:0);D组合一(Z:I:P=0.1:0.2:0.7);E组合二(Z:I:P=0.2:0.7:0.1);F组合三(Z:I:P=0.7:0.1:0.2)。从而计算得到的各种组合模型电压幅值如图4所示。可以发现随着ZIP负荷模型中恒阻抗系数的增加,将会逐渐提升系统的电压水平,且对系统末端节点的幅度提升最为明显,并且各种分量比值的变化,都将影响着系统电压幅值的变化。因此,在线性潮流计算过程中,对负荷电压静特性的考虑也是影响潮流分布的重要因素之一。图5为ZIP负荷模型电压幅值和相角误差分布图,从图5中不难看出,不同比例系数的ZIP负荷模型会对电压幅值和相角精度造成不同的影响。其中恒功率负荷模型可以使得电压精度达到最高,而随着恒阻抗系数和恒电流系数的增加,电压的误差也会变得越来越大。但不管比例系数如何变化,本发明方法的电压幅值最大误差仍可以保持在10-3~10-5,相角最大误差仍可以保持在10-2。由此可见,本发明所提方法在处理含ZIP负荷模型的配网系统时,仍具有足够的计算精度。
为了分析不同负载比重对本发明算法精度的影响,通过增加系统中所有线路的负载,来模拟均匀过载的病态系统运行情况。本发明中系统负载比重的调整范围为1.5~3.0。同时,在IEEE33节点系统的基础上闭合8-21、9-15两个联络开关来组成环网进行测试,测试结果如表1所示。可以发现,本发明在高负载情况下仍然可以保持较高的精度。
表1线性潮流算法在不同负载比重下的误差值
在考虑ZIP负荷模型时,两种潮流计算方法在不同的配电系统中的运行时间如表2所示。可以发现,在考虑ZIP负荷模型的情况下,相比较于BFS,本发明将测试系统的运行时间降低了80%左右,计算效率得到了显著提高。
表2潮流算法的运行时间
将本发明在IEEE69、IEEE85、IEEE141节点系统进行测试,表3为各系统电压幅值最大误差和相角最大误差。可以发现,本发明线性法在大型的配电网系统中电压幅值误差的数量级始终可以保持在10-4,相角误差的数量级可以保持在10-2,具有较高的精度。
表3线性潮流算法在不同测试系统中的误差值
Claims (9)
1.考虑负荷电压静特性的配电网线性潮流计算方法,其特征在于包括以下步骤:
步骤1:结合配电网的典型特征,定量分析适用于配电网的简化条件;
步骤2:根据电压的变化将负荷模型划分为恒阻抗分量、恒电流分量和恒功率分量,建立一种考虑电压静特性的ZIP负荷模型,并根据配电网的简化条件,对其线性化处理;
步骤3:将潮流模型存在的非线性项1/Vi,通过最大和最小电压加权因子进行加权,得到一个线性的广义函数,并结合配电网的简化条件,将潮流模型简化成一种全线性化潮流模型;
步骤4:对全线性化潮流模型进行初等变换,实现对电压幅值和相位角的解耦;
步骤5:根据配电网各节点的已知量,求出网络中各个节点的电压和相角。
2.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流计算方法,其特征在于:所述步骤1中,配电网的典型特征包括:①、节点电压幅值趋近于1.0p.u.;②、线路两端的相角非常小,使得配电网中所有节点的电压相角与平衡节点的电压相角相差不大;③、R/X比值比较大,接近或者大于1。
3.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流计算方法,其特征在于:所述步骤1中,配电网的简化条件,如公式(1)所示:
4.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流计算方法,其特征在于:所述步骤2中,根据电压的变化,建立一种考虑电压静特性的ZIP负荷模型,如公式(5)所示:
其中,P(V),Q(V)分别表示节点负载有功功率和无功功率;V,VN分别为节点实际电压和额定电压;PN,QN分别表示额定电压下的有功功率和无功功率;CZ,CI,CP分别表示节点有功恒阻抗、恒电流、恒功率负荷的比例系数;C′Z,C′I,C′P分别表示节点无功恒阻抗、恒电流、恒功率负荷的比例系数;各参数满足的约束条件为CZ+CI+CP=1,C′Z+C′I+C′P=1。
5.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流计算方法,其特征在于:所述步骤2中,ZIP负荷模型线性化处理,以有功功率为例,如公式(2)所示:
其中,P(V)表示节点负载有功功率;PN分别表示额定电压下的有功功率;CZ,CI,CP分别表示节点有功恒阻抗、恒电流、恒功率负荷的比例系数,V为节点实际电压,△V为实际电压与额定电压的差值,通常取值为0~0.1p.u.,其二阶及高阶项的取值相对较小,潮流计算中可忽略不计;PP,PI分别为ZIP负荷模型线性化后恒功率负荷和恒电流负荷新的比例系数;同理,无功功率线性化处理的方法与有功功率一致。
6.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流算法研究,其特征在于:所述步骤3中,线性的广义函数K如公式(4)所示:
其中,i为节点的编号;Vi为节点i的电压幅值;Vmin,Vmax分别为电压最小加权因子和电压最大加权因子,要求Vmin≤Vi≤Vmax。
7.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流计算方法,其特征在于:所述步骤3中,全线性潮流方程如公式(7)所示:
其中,Vmin,Vmax分别为电压最小加权因子和电压最大加权因子;S,W,R分别为配电网中的平衡节点、PQ节点、PV节点的集合;PPS,PPW,QPW分别平衡节点有功功率、PQ节点有功功率、PQ节点无功功率的恒功率系数;PIS,PIW,QIW分别为平衡节点有功功率、PQ节点有功功率、PQ节点无功功率的恒电流系数;VS、δS分别为平衡节点电压幅值和相角;VW、δW分别为PQ节点电压幅值和相角;VR、δR分别为PV节点电压幅值和相角;diag(PPR)、diag(PPW)、diag(QPW)分别为PV节点有功功率的恒功率负荷系数对角矩阵、PQ节点有功功率的恒功率负荷系数对角矩阵、PV节点无功功率的恒功率负荷系数对角矩阵;GRS和BRS分别为PV节点与平衡节点的互导纳矩阵的实部和虚部;GWS和BWS分别为PQ节点与平衡节点的互导纳矩阵的实部和虚部;GWR和BWR分别为PQ节点与PV节点的互导纳矩阵的实部和虚部;GRR和BRR为PV节点自导纳矩阵的实部和虚部;BRW和GRW分别为PQ节点与PV节点的互导纳矩阵的实部和虚部;GWW和BWW分别为PQ节点自导纳矩阵的实部和虚部。
8.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流计算方法,其特征在于:所述步骤4中,对电压幅值和相位角的解耦模型,如公式(8)所示:
其中,为已知项,可以通过节点导纳矩阵和已知节点功率求出;为已知项,可通过以平衡节点、PQ节点、PV节点所组成的导纳矩阵求出;为系统PQ节点和PV节点的相角集合,为PQ节点的电压集合,两者均为未知量。
9.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流计算方法,其特征在于:所述步骤5中,计算网络中各个节点的电压和相角模型,如公式(9)所示:
JX=B (9)
其中,B为(n+m-2)×1阶矩阵,其元素可由节点注入功率、导纳矩阵、平衡节点的电压幅值和PV节点的电压幅值求出,J为(n+m-2)×(n+m-2)阶方阵,其元素可由节点注入功率和导纳矩阵求得;X为系统未知电压和相角的集合;该模型的个数总和为m+n-2刚好等于未知量的个数,因此,可求出网络中各个节点的电压和相角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910719366.7A CN110518590B (zh) | 2019-08-05 | 2019-08-05 | 考虑负荷电压静特性的配电网线性潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910719366.7A CN110518590B (zh) | 2019-08-05 | 2019-08-05 | 考虑负荷电压静特性的配电网线性潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110518590A true CN110518590A (zh) | 2019-11-29 |
CN110518590B CN110518590B (zh) | 2022-12-20 |
Family
ID=68625268
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910719366.7A Active CN110518590B (zh) | 2019-08-05 | 2019-08-05 | 考虑负荷电压静特性的配电网线性潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110518590B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111541252A (zh) * | 2020-06-02 | 2020-08-14 | 国网重庆市电力公司 | 基于历史数据和误差分析的配电网线性化潮流通式最优独立变量选取方法 |
CN111934320A (zh) * | 2020-08-20 | 2020-11-13 | 三峡大学 | 基于线性潮流方程和改进烟花算法的主动配电网快速重构方法 |
CN112290604A (zh) * | 2020-10-15 | 2021-01-29 | 珠海博威电气股份有限公司 | 计及负荷特性的配电网协调优化方法和装置、存储介质 |
CN113224780A (zh) * | 2021-03-23 | 2021-08-06 | 重庆大学 | 一种双极直流配电网潮流的线性化计算方法 |
CN114188948A (zh) * | 2021-12-14 | 2022-03-15 | 国网福建省电力有限公司 | 一种考虑不对称三角形恒功率负荷的潮流计算方法 |
CN114784885A (zh) * | 2022-03-23 | 2022-07-22 | 华北水利水电大学 | 基于Zbus矩阵的梯度法配电网分布式电源选址定容方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018049737A1 (zh) * | 2016-09-18 | 2018-03-22 | 国电南瑞科技股份有限公司 | 一种基于分区负荷控制的安全校正计算方法 |
CN108718091A (zh) * | 2018-07-09 | 2018-10-30 | 国网福建省电力有限公司 | 一种应用于主动配电网的三相极坐标系线性潮流计算方法 |
CN109995075A (zh) * | 2018-09-29 | 2019-07-09 | 国网上海市电力公司 | 一种含分布式电源的主动配电网动态重构方法 |
-
2019
- 2019-08-05 CN CN201910719366.7A patent/CN110518590B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2018049737A1 (zh) * | 2016-09-18 | 2018-03-22 | 国电南瑞科技股份有限公司 | 一种基于分区负荷控制的安全校正计算方法 |
CN108718091A (zh) * | 2018-07-09 | 2018-10-30 | 国网福建省电力有限公司 | 一种应用于主动配电网的三相极坐标系线性潮流计算方法 |
CN109995075A (zh) * | 2018-09-29 | 2019-07-09 | 国网上海市电力公司 | 一种含分布式电源的主动配电网动态重构方法 |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111541252A (zh) * | 2020-06-02 | 2020-08-14 | 国网重庆市电力公司 | 基于历史数据和误差分析的配电网线性化潮流通式最优独立变量选取方法 |
CN111541252B (zh) * | 2020-06-02 | 2023-12-12 | 国网重庆市电力公司 | 基于历史数据和误差分析的配电网线性化潮流通式最优独立变量选取方法 |
CN111934320A (zh) * | 2020-08-20 | 2020-11-13 | 三峡大学 | 基于线性潮流方程和改进烟花算法的主动配电网快速重构方法 |
CN111934320B (zh) * | 2020-08-20 | 2023-08-01 | 三峡大学 | 基于线性潮流方程和改进烟花算法的主动配电网快速重构方法 |
CN112290604A (zh) * | 2020-10-15 | 2021-01-29 | 珠海博威电气股份有限公司 | 计及负荷特性的配电网协调优化方法和装置、存储介质 |
CN113224780A (zh) * | 2021-03-23 | 2021-08-06 | 重庆大学 | 一种双极直流配电网潮流的线性化计算方法 |
CN113224780B (zh) * | 2021-03-23 | 2023-03-14 | 重庆大学 | 一种双极直流配电网潮流的线性化计算方法 |
CN114188948A (zh) * | 2021-12-14 | 2022-03-15 | 国网福建省电力有限公司 | 一种考虑不对称三角形恒功率负荷的潮流计算方法 |
CN114188948B (zh) * | 2021-12-14 | 2023-06-23 | 国网福建省电力有限公司 | 一种考虑不对称三角形恒功率负荷的潮流计算方法 |
CN114784885A (zh) * | 2022-03-23 | 2022-07-22 | 华北水利水电大学 | 基于Zbus矩阵的梯度法配电网分布式电源选址定容方法 |
CN114784885B (zh) * | 2022-03-23 | 2024-06-18 | 华北水利水电大学 | 基于Zbus矩阵的梯度法配电网分布式电源选址定容方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110518590B (zh) | 2022-12-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110518590B (zh) | 考虑负荷电压静特性的配电网线性潮流计算方法 | |
CN101635457B (zh) | 一种基于状态估计残差参数灵敏度的电网参数估计方法 | |
Khazraj et al. | A performance comparison between extended Kalman Filter and unscented Kalman Filter in power system dynamic state estimation | |
Wang et al. | Multi-stage holomorphic embedding method for calculating the power-voltage curve | |
Ma et al. | A hybrid dynamic equivalent using ANN-based boundary matching technique | |
CN107565553A (zh) | 一种基于ukf的配电网抗差动态状态估计方法 | |
CN108565852B (zh) | 一种三阶段递进的大电网电压稳定评估的故障筛选与排序方法 | |
CN108054757B (zh) | 一种内嵌无功和电压的n-1闭环安全校核方法 | |
CN110308366B (zh) | 一种基于正交匹配追踪算法的谐波源定位方法 | |
CN110829444B (zh) | 计及静态频率和电压特性的电网受端系统紧急切负荷方法 | |
CN104794531A (zh) | 基于网损等值负荷模型的改进直流最优潮流方法 | |
CN109617080B (zh) | 基于改进的雅可比矩阵的直角坐标牛顿法潮流计算方法 | |
CN106130025A (zh) | 一种基于回路分析法计算三相配电网潮流的线性方法 | |
CN109446643B (zh) | 一种基于实测数据建立家用电器负荷谐波模型的方法 | |
CN109698505B (zh) | 大电网静态电压稳定在线防控的调控量化映射计算方法 | |
CN110518627A (zh) | 计及风电不确定性的电力系统静态电压稳定概率评估方法 | |
CN106532712A (zh) | 含小阻抗支路电网的补偿法直角坐标牛顿法潮流计算方法 | |
Andreoli et al. | Topology validation via simultaneous state & topology estimation with phasor data processing capability | |
Jo et al. | Kalman-filter-based multilevel analysis to estimate electric load composition | |
CN106953331A (zh) | 一种考虑不确定性和三相不平衡的电力系统谐波潮流算法 | |
CN116599083A (zh) | 直流微电网高频振荡稳定性风险评估方法 | |
CN110781638B (zh) | 基于节点端口伏安特性的最优戴维南等值参数计算方法 | |
Amouzad Mahdiraji et al. | A new method for simplification and reduction of state estimation’s computational complexity in stability analysis of power systems | |
Al Ameri et al. | Fast estimation method for selection of optimal distributed generation size using Kalman filter and graph theory | |
CN113794198A (zh) | 抑制宽频振荡的方法、装置、终端及存储介质 |
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 |