CN110518590B - 考虑负荷电压静特性的配电网线性潮流计算方法 - Google Patents

考虑负荷电压静特性的配电网线性潮流计算方法 Download PDF

Info

Publication number
CN110518590B
CN110518590B CN201910719366.7A CN201910719366A CN110518590B CN 110518590 B CN110518590 B CN 110518590B CN 201910719366 A CN201910719366 A CN 201910719366A CN 110518590 B CN110518590 B CN 110518590B
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.)
Active
Application number
CN201910719366.7A
Other languages
English (en)
Other versions
CN110518590A (zh
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.)
China Three Gorges University CTGU
Original Assignee
China Three Gorges University CTGU
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 China Three Gorges University CTGU filed Critical China Three Gorges University CTGU
Priority to CN201910719366.7A priority Critical patent/CN110518590B/zh
Publication of CN110518590A publication Critical patent/CN110518590A/zh
Application granted granted Critical
Publication of CN110518590B publication Critical patent/CN110518590B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • 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

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)所示：
Figure BDA0002155043200000021
所述步骤2中，负荷模型是指当负荷的实际电压偏离额定电压时，负荷的特性也将发生变化。因此，根据电压的变化可将负荷模型划分为恒阻抗分量、恒电流分量和恒功率分量。其中，恒阻抗分量和恒电流分量分别指负荷模型的电压二次项和电压一次项，恒功率分量是指负荷模型的常数项。
所述步骤2中，根据电压的变化，建立一种考虑电压静特性的ZIP负荷模型,如公式(5)所示：
Figure BDA0002155043200000031
其中，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)所示：
Figure BDA0002155043200000032
其中，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)所示：
Figure BDA0002155043200000033
所述步骤3中，线性的广义函数K如公式(4)所示：
Figure BDA0002155043200000041
其中，i为节点的编号；Vi为节点i的电压幅值；
Figure BDA0002155043200000042
Vmin，Vmax分别为电压最小加权因子和电压最大加权因子，要求Vmin≤Vi≤Vmax。
所述步骤3中，全线性潮流方程如公式(7)所示：
Figure BDA0002155043200000043
其中，
Figure BDA0002155043200000044
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)所示：
Figure BDA0002155043200000045
其中，
Figure BDA0002155043200000046
为已知项，可以通过节点导纳矩阵和已知节点功率求出；
Figure BDA0002155043200000047
为已知项，可通过以平衡节点、PQ节点、PV节点所组成的导纳矩阵求出；
Figure BDA0002155043200000051
为系统PQ节点和PV节点的相角集合，
Figure BDA0002155043200000052
为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)所示:
Figure BDA0002155043200000061
其中，i和j表示节点的编号；δi，δj分别为节点i,j的相角；δij为节点i,j的相角差，Vj为节点j的电压幅值。
2)基于电压静特性的ZIP负荷模型存在电压的平方项Vi 2。则定义Vi＝1-△Vi，通过忽略二次项△Vi 2得到其近似线性化，如公式(2)所示：
Figure BDA0002155043200000062
其中，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)所示：
Figure BDA0002155043200000071
其中，Vi为节点i的电压幅值；
Figure BDA0002155043200000072
Vmin，Vmax分别为电压最小加权因子和电压最大加权因子，要求Vmin≤Vi≤Vmax。
负荷和恒电流负荷新的比例系数。
2、所述建立一种考虑电压静特性的ZIP负荷模型，如公式(5)所示：
Figure BDA0002155043200000073
其中，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)所示：
Figure BDA0002155043200000074
式中：i和j表示节点的编号；Vi和Vj分别为节点i和j的电压幅值；Pi和Qi分别表示配电网中节点i的总注入有功功率和无功功率；PLi和QLi分别表示节点i负载有功功率和无功功率。PGi和QGi表示DG的有功功率和无功功率；Gij和Bij分别表示导纳矩阵的实部和虚部；δij为节点i,j的相角差。
4、简化后的全线性潮流方程如公式(7)所示：
Figure BDA0002155043200000081
其中，
Figure BDA0002155043200000082
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)所示：
Figure BDA0002155043200000083
其中，
Figure BDA0002155043200000084
为已知项，可以通过节点导纳矩阵和已知节点功率求出；
Figure BDA0002155043200000085
为已知项，可通过以平衡节点、PQ节点、PV节点所组成的导纳矩阵求出；
Figure BDA0002155043200000086
为系统PQ节点和PV节点的相角集合，
Figure BDA0002155043200000087
为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)
Figure BDA0002155043200000091
其中，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线性潮流算法在不同负载比重下的误差值
Figure BDA0002155043200000101
在考虑ZIP负荷模型时，两种潮流计算方法在不同的配电系统中的运行时间如表2所示。可以发现，在考虑ZIP负荷模型的情况下，相比较于BFS，本发明将测试系统的运行时间降低了80％左右，计算效率得到了显著提高。
表2潮流算法的运行时间
Figure BDA0002155043200000102
将本发明在IEEE69、IEEE85、IEEE141节点系统进行测试，表3为各系统电压幅值最大误差和相角最大误差。可以发现，本发明线性法在大型的配电网系统中电压幅值误差的数量级始终可以保持在10-4，相角误差的数量级可以保持在10-2，具有较高的精度。
表3线性潮流算法在不同测试系统中的误差值
Figure BDA0002155043200000103

Claims (9)

1.考虑负荷电压静特性的配电网线性潮流计算方法，其特征在于包括以下步骤：
步骤1：结合配电网的典型特征，定量分析适用于配电网的简化条件；
步骤2：根据电压的变化将负荷模型划分为恒阻抗分量、恒电流分量和恒功率分量，建立一种考虑电压静特性的ZIP负荷模型，并根据配电网的简化条件，对其线性化处理；
步骤3：将潮流模型存在的非线性项1/Vi，通过最大和最小电压加权因子进行加权，得到一个线性的广义函数，并结合配电网的简化条件，将潮流模型简化成一种全线性化潮流模型；
步骤4：对全线性化潮流模型进行初等变换，实现对电压幅值和相位角的解耦；
步骤5：根据配电网各节点的已知量，求出网络中各个节点的电压和相角。
2.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流计算方法，其特征在于：所述步骤1中，配电网的典型特征包括：①、节点电压幅值趋近于1.0p.u.；②、线路两端的相角非常小，使得配电网中所有节点的电压相角与平衡节点的电压相角相差不大；③、R/X比值比较大，接近或者大于1。
3.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流计算方法，其特征在于：所述步骤1中，配电网的简化条件，如公式(1)所示：
Figure FDA0003921179660000011
4.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流计算方法，其特征在于：所述步骤2中，根据电压的变化，建立一种考虑电压静特性的ZIP负荷模型,如公式(5)所示：
Figure FDA0003921179660000012
其中，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)所示：
Figure FDA0003921179660000021
其中，P(V)表示节点负载有功功率；PN分别表示额定电压下的有功功率；CZ,CI,CP分别表示节点有功恒阻抗、恒电流、恒功率负荷的比例系数，V为节点实际电压，△V为实际电压与额定电压的差值，通常取值为0～0.1p.u.，其二阶及高阶项的取值相对较小，潮流计算中可忽略不计；PP，PI分别为ZIP负荷模型线性化后恒功率负荷和恒电流负荷新的比例系数；同理，无功功率线性化处理的方法与有功功率一致。
6.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流计算方法，其特征在于：所述步骤3中，线性的广义函数K如公式(4)所示：
Figure FDA0003921179660000022
其中，i为节点的编号；Vi为节点i的电压幅值；
Figure FDA0003921179660000023
Vmin，Vmax分别为电压最小加权因子和电压最大加权因子，要求Vmin≤Vi≤Vmax。
7.根据权利要求1所述考虑负荷电压静特性的配电网线性潮流计算方法，其特征在于：所述步骤3中，全线性潮流方程如公式(7)所示：
Figure FDA0003921179660000024
其中，
Figure FDA0003921179660000025
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)所示：
Figure FDA0003921179660000031
其中，
Figure FDA0003921179660000032
为已知项，可以通过节点导纳矩阵和已知节点功率求出；
Figure FDA0003921179660000033
为已知项，可通过以平衡节点、PQ节点、PV节点所组成的导纳矩阵求出；
Figure FDA0003921179660000034
为系统PQ节点和PV节点的相角集合，
Figure FDA0003921179660000035
为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刚好等于未知量的个数，因此，可求出网络中各个节点的电压和相角。
CN201910719366.7A 2019-08-05 2019-08-05 考虑负荷电压静特性的配电网线性潮流计算方法 Active CN110518590B (zh)

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 CN110518590A (zh) 2019-11-29
CN110518590B true 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)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111541252B (zh) * 2020-06-02 2023-12-12 国网重庆市电力公司 基于历史数据和误差分析的配电网线性化潮流通式最优独立变量选取方法
CN111934320B (zh) * 2020-08-20 2023-08-01 三峡大学 基于线性潮流方程和改进烟花算法的主动配电网快速重构方法
CN112290604B (zh) * 2020-10-15 2023-06-30 珠海博威电气股份有限公司 计及负荷特性的配电网协调优化方法和装置、存储介质
CN113224780B (zh) * 2021-03-23 2023-03-14 重庆大学 一种双极直流配电网潮流的线性化计算方法
CN114188948B (zh) * 2021-12-14 2023-06-23 国网福建省电力有限公司 一种考虑不对称三角形恒功率负荷的潮流计算方法
CN114784885A (zh) * 2022-03-23 2022-07-22 华北水利水电大学 基于Zbus矩阵的梯度法配电网分布式电源选址定容方法

Citations (3)

* Cited by examiner, † Cited by third party
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 国网上海市电力公司 一种含分布式电源的主动配电网动态重构方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
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 国网上海市电力公司 一种含分布式电源的主动配电网动态重构方法

Also Published As

Publication number Publication date
CN110518590A (zh) 2019-11-29

Similar Documents

Publication Publication Date Title
CN110518590B (zh) 考虑负荷电压静特性的配电网线性潮流计算方法
CN101635457B (zh) 一种基于状态估计残差参数灵敏度的电网参数估计方法
Moghavvemi et al. Technique for assessment of voltage stability in ill-conditioned radial distribution network
CN110957728B (zh) 一种三相有源配电网的全纯嵌入潮流方法及装置
CN107565553A (zh) 一种基于ukf的配电网抗差动态状态估计方法
CN108565852B (zh) 一种三阶段递进的大电网电压稳定评估的故障筛选与排序方法
CN110829444B (zh) 计及静态频率和电压特性的电网受端系统紧急切负荷方法
CN108054757A (zh) 一种内嵌无功和电压的n-1闭环安全校核方法
CN112653134B (zh) 基于helm电压灵敏度的配电网电压稳定判定方法
Syai'in et al. Real-time unbalanced load flow development using direct-ZBR method and modified Lambda iteration for on-line monitoring and control
Schwerdtner et al. Adaptive sampling for structure-preserving model order reduction of port-Hamiltonian systems
CN113283784B (zh) 基于区间优化算法的可靠性逆问题模型及其求解方法
CN107994567B (zh) 一种广义快速分解状态估计方法
Karimi et al. Impact of load modeling in distribution state estimation
CN109888788B (zh) 一种电力系统最优潮流的求解方法
CN109193639B (zh) 一种电力系统抗差估计方法
CN112511056A (zh) 一种基于相量测量的鲁棒发电机动态状态估计方法
Zhu et al. Effect of load modeling on voltage stability
CN110632847A (zh) 一种自耦pd协同控制理论新方法
Savin et al. Covariance based uncertainty analysis with unscented transformation
CN110412353B (zh) 一种基于功率损耗的光伏电站三相变压器阻抗计算方法
Pazderin et al. An energy flow calculation method based on energy meters data
CN107658880B (zh) 基于关联矩阵运算的快速分解法系数矩阵计算方法
Fu et al. Analysis of Subsynchronous Oscillation Caused by Multiple VSCs with Different Dynamics Under Strong Grid Connections
Kataoka A smooth power flow model of electric power system with generator reactive power limits taken into consideration

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