CN110460059A - 一种电力系统潮流计算方法 - Google Patents
一种电力系统潮流计算方法 Download PDFInfo
- Publication number
- CN110460059A CN110460059A CN201910684342.2A CN201910684342A CN110460059A CN 110460059 A CN110460059 A CN 110460059A CN 201910684342 A CN201910684342 A CN 201910684342A CN 110460059 A CN110460059 A CN 110460059A
- Authority
- CN
- China
- Prior art keywords
- load
- power
- node
- value
- voltage
- 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
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
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
本发明提供一种电力系统潮流计算方法,包括如下步骤:根据潮流计算模型中各负荷节点变电站最近一次有载调压分接开关动作引起的电气量变化值确定ZP负荷模型参数,在形成节点导纳矩阵时,将各负荷节点ZP负荷模型恒阻抗分量的导纳值加入节点导纳矩阵对应对角元,在设定各节点电压初值后,采用ZP负荷模型中恒功率分量功率值计算负荷节点恒功率分量功率的不平衡量和电压不平衡量,进行潮流计算。
Description
技术领域
本发明涉及电力系统分析计算技术领域,具体涉及一种电力系统潮流计算方法。
背景技术
电力系统潮流计算是研究电力系统稳态运行情况的一种基本电气计算。它的任务是根据给定的运行条件和网络结构确定整个系统的运行状态,如各母线上的电压(幅值及相角)、网络中的功率分布以及功率损耗等。电力系统潮流计算的结果是电力系统稳定计算和故障分析的基础,通过电力系统潮流计算已比较运行方式或规划供电方案的可行性、可靠性和经济性。同时,为了实时监控电力系统的运行状态,也需要进行大量而快速的潮流计算。
基于牛顿-拉夫逊的潮流计算是电力系统分析计算的一种常用方法,该方法通常将负荷节点设置为PQ节点,即负荷的有功功率和无功功率为恒定值。而在实际系统下,负荷的有功功率和无功功率随负荷节点电压值的变化而变化,因此采用恒定功率负荷模型进行潮流计算得到的结果并不符合实际。
传统的牛顿-拉夫逊在处理负荷节点时主要的不足在于:在计算功率的不平衡量时,将负荷功率与电压初值计算功率做差值,没有考虑到负荷的功率-电压特性。
发明内容
针对传统的牛顿-拉夫逊在处理负荷节点时在计算功率的不平衡量时,将负荷功率与电压初值计算功率做差值,没有考虑到负荷的功率-电压特性,分析结果不准确的问题,本发明提供一种电力系统潮流计算方法。
本发明的技术方案是:
一种电力系统潮流计算方法,包括如下步骤:根据潮流计算模型中各负荷节点变电站最近一次有载调压分接开关动作引起的电气量变化值确定ZP负荷模型参数,在形成节点导纳矩阵时,将各负荷节点ZP负荷模型恒阻抗分量的导纳值加入节点导纳矩阵对应对角元,在设定各节点电压初值后,采用ZP负荷模型中恒功率分量功率值计算负荷节点恒功率分量功率的不平衡量和电压不平衡量,进行潮流计算。
优选地,所述的根据潮流计算模型中各负荷节点变电站最近一次有载调压分接开关动作引起的电气量变化值确定ZP负荷模型参数的步骤包括:
根据有载调压分接开关变化引起变电站低压侧的特征量计算ZP负荷模型中恒阻抗分量消耗的有功、无功功率占负荷功率的比值bp、bq:
其中,式中UA、UB为主变最近一次有载调压分接开关动作前后主变低压侧电压值,ΔP1为有载调压分接开关动作后低压侧有功功率瞬时变化值,ΔP2为功率瞬时变化后恢复到功率稳定的变化量,ΔQ1为有载调压分接开关变化后低压侧无功功率瞬时变化值,ΔQ2为功率瞬时变化后恢复到功率稳定的变化量,Pi、Qi为i节点负荷有功、无功功率。
优选地,根据ZP负荷模型参数计算节点导纳矩阵的具体步骤包括:
设定N个节点的电力网络,其网络方程可表示为:其中,为节点电压列向量,为节点注入电流列向量,Y为N阶节点导纳矩阵,其非对角元Yij可表示为:其中,zij为节点i,j的支路阻抗,且有i≠j;
非负荷节点的对角元Yii等于各个节点与该节点所连导纳的总和:当j=0时,zij为节点对地导纳;
负荷节点i的对角元Yii等于各个节点与该节点所连导纳的总和再加上ZP负荷模型中恒阻抗分量的导纳值:
式中,i≠j,当j=0时,zij为节点对地导纳,YLi可由下式计算得出:
式中,Ui为负荷节点电压值,Pi、Qi为负荷功率,bp、bq该节点ZP负荷模型中恒阻抗分量消耗的有功、无功功率占负荷功率的比值。
优选地,设定各节点电压初始值的步骤包括:
设定各节点电压初值实部虚部fi (0),通常设定为节点电压额定值,fi (0)为0。
优选地,根据设定的电压初始值或电压修正值计算负荷节点功率不平衡量和电压不平衡量的步骤包括:
将各节点电压初值或上次修正计算值代入下式ei、fi,求负荷节点恒功率分量功率的不平衡量ΔPi、ΔQi和电压不平衡量ΔUi:
式中,(1-bp)Pi、(1-bq)Qi分别为去除ZP负荷模型中恒阻抗分量负荷的负荷有功、无功功率,ei、fi分别为节点电压实部、虚部,Gij为节点导纳矩阵第i行第j列元素实部,Bij为节点导纳矩阵第i行第j列元素虚部。
优选地,计算各节点电压修正量及电压修正值的步骤具体包括:
将各节点电压初值代入,计算修正方程式系数矩阵元素,系数矩阵元素表达式如下:
当j=i
当j≠i
Rij=Sij=0
计算各节点电压的修正量:
计算各节点电压修正值公式如下:
优选地,判断是否满足收敛条件具体为:
若不满足上式的收敛条件,则以各节点电压修正值作为初值,返回步骤:根据设定的电压初始值或电压修正值计算负荷节点功率不平衡量和电压不平衡量进行下一轮迭代,ε为收敛精度,可根据实际要求设定。
从以上技术方案可以看出,本发明具有以下优点:实际工作过程中,由于负荷存在功率-电压特性,当系统电压变化时负荷功率随之变化,本方法解决了传统的方法当电压变化时,负荷功率保持不变的问题,本发明计算潮流结果负荷节点功率随电压降低而降低,较为符合实际情况,为后续的电力系统的分析提供可靠的参考。
此外,本发明设计原理可靠,结构简单,具有非常广泛的应用前景。
由此可见,本发明与现有技术相比,具有突出的实质性特点和显著地进步,其实施的有益效果也是显而易见的。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,对于本领域普通技术人员而言,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的一种电力系统稳定性分析方法的示意性流程图。
图2是本发明实施例的算例模型图。
具体实施方式
为了使本技术领域的人员更好地理解本发明中的技术方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
如图1所示,一种电力系统潮流计算方法,包括如下步骤:
S1:确定ZP负荷模型参数;
本步骤中,利用各负荷节点变电站与当前时间最近一次有载调压分接开关变化确定ZP负荷模型参数,步骤具体包括:
根据有载调压分接开关变化引起变电站低压侧的特征量计算ZP负荷模型中恒阻抗分量消耗的有功、无功功率占负荷功率的比值bp、bq:
其中,式中UA、UB为主变最近一次有载调压分接开关动作前后主变低压侧电压值,ΔP1为有载调压分接开关动作后低压侧有功功率瞬时变化值,ΔP2为功率瞬时变化后恢复到功率稳定的变化量,ΔQ1为有载调压分接开关变化后低压侧无功功率瞬时变化值,ΔQ2为功率瞬时变化后恢复到功率稳定的变化量,Pi、Qi为i节点负荷有功、无功功率。
S2:根据ZP负荷模型参数计算节点导纳矩阵;
具体步骤包括:
设定N个节点的电力网络,其网络方程可表示为:其中,为节点电压列向量,为节点注入电流列向量,Y为N阶节点导纳矩阵,其非对角元Yij可表示为:其中,zij为节点i,j的支路阻抗,且有i≠j;
非负荷节点的对角元Yii等于各个节点与该节点所连导纳的总和:当j=0时,zij为节点对地导纳;
负荷节点i的对角元Yii等于各个节点与该节点所连导纳的总和再加上ZP负荷模型中恒阻抗分量的导纳值:
式中,i≠j,当j=0时,zij为节点对地导纳,YLi可由下式计算得出:
式中,Ui为负荷节点电压值,Pi、Qi为负荷功率,bp、bq该节点ZP负荷模型中恒阻抗分量消耗的有功、无功功率占负荷功率的比值。
S3:设定各节点电压初始值;设定各节点电压初值实部虚部fi (0),通常设定为节点电压额定值,fi (0)为0。
S4:根据设定的电压初始值或电压修正值计算负荷节点功率不平衡量和电压不平衡量;
将各节点电压初值或上次修正计算值代入下式ei、fi,求负荷节点恒功率分量功率的不平衡量ΔPi、ΔQi和电压不平衡量ΔUi:
式中,(1-bp)Pi、(1-bq)Qi分别为去除ZP负荷模型中恒阻抗分量负荷的负荷有功、无功功率,ei、fi分别为节点电压实部、虚部,Gij为节点导纳矩阵第i行第j列元素实部,Bij为节点导纳矩阵第i行第j列元素虚部。
S5:将各节点电压初值代入,求计算修正方程式系数矩阵元素;
将各节点电压初值代入,计算修正方程式系数矩阵元素,系数矩阵元素表达式如下:
当j=i
当j≠i
Rij=Sij=0
S6:计算各节点电压修正量;
计算各节点电压的修正量:
S7:计算电压修正值修正各节点电压初值;计算各节点电压修正值公式如下:
S8:判断是否满足收敛条件若是,执行步骤S9;
若不满足上式的收敛条件,则各节点电压修正值作为初值,返回步骤S4:根据电压修正值计算负荷节点恒功率分量功率不平衡量和电压不平衡量进行下一轮迭代,ε为收敛精度,可根据实际要求设定。
S9:输出计算结果。
实施例1
本发明实施例如图2所示选用IEEE9节点系统来进行测试,对比分析现有技术和本发明潮流计算结果。
表1给出了各负荷节点ZP负荷模型参数;表2给出了BUS1发电机电压为100%额定值时,现有技术潮流计算功率、电压值;表3给出了BUS1发电机电压为90%额定值时,现有技术和本发明潮流计算功率、电压值比较。
表1各负荷节点ZP负荷模型参数
表2 BUS1发电机电压为额定值100%时潮流计算电压值比较
表3 BUS1发电机电压为额定值90%时潮流计算电压值比较
由表2和表3可看出,降低BUS1发电机电压后,系统各节点电压均下降。由于负荷存在功率-电压特性,当电压降低时,负荷功率响应降低。比较表2、表3现有技术计算潮流结果,负荷节点BUS5、BUS6、BUS8有功、无功功率并没有随电压降低而降低,本发明计算潮流结果负荷节点功率随电压降低而降低,较为符合实际情况。
尽管通过参考附图并结合优选实施例的方式对本发明进行了详细描述,但本发明并不限于此。在不脱离本发明的精神和实质的前提下,本领域普通技术人员可以对本发明的实施例进行各种等效的修改或替换,而这些修改或替换都应在本发明的涵盖范围内/任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (4)
1.一种电力系统潮流计算方法,其特征在于:根据潮流计算模型中各负荷节点变电站最近一次有载调压分接开关动作引起的电气量变化值确定ZP负荷模型参数,在形成节点导纳矩阵时,将各负荷节点ZP负荷模型恒阻抗分量的导纳值加入节点导纳矩阵对应对角元,在设定各节点电压初值后,采用ZP负荷模型中恒功率分量功率值计算负荷节点恒功率分量功率的不平衡量和电压不平衡量,进行潮流计算。
2.如权利要求1所述的一种电力系统潮流计算方法,其特征在于:对于潮流计算中的负荷节点i,其对应ZP负荷模型中恒阻抗分量消耗的有功、无功功率占负荷功率的比值bp、bq可根据下式计算:
式(1)、式(2)中UA、UB为主变最近一次有载调压分接开关动作前后主变低压侧电压值,ΔP1为有载调压分接开关动作后低压侧有功功率瞬时变化值,ΔP2为功率瞬时变化后恢复到功率稳定的变化量,ΔQ1为有载调压分接开关变化后低压侧无功功率瞬时变化值,ΔQ2为功率瞬时变化后恢复到功率稳定的变化量,Pi、Qi为i节点负荷有功、无功功率。
3.如权利要求1所述的一种电力系统潮流计算方法,其特征在于:在形成节点导纳矩阵时,对于负荷节点i,其对角元Yii等于各个节点与该节点所连导纳的总和加上ZP负荷模型中恒阻抗分量的导纳值。
4.如权利要求2所述的一种电力系统潮流计算方法,其特征在于:在计算功率不平衡量时,采用ZP负荷模型中恒功率分量功率值计算负荷节点恒功率分量功率的不平衡量
式中,(1-bp)Pi、(1-bq)Qi分别为去除ZP负荷模型中恒阻抗分量负荷的负荷有功、无功功率,ei、fi分别为节点电压实部、虚部,Gij为节点导纳矩阵第i行第j列元素实部,Bij为节点导纳矩阵第i行第j列元素虚部。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910684342.2A CN110460059B (zh) | 2019-07-26 | 2019-07-26 | 一种电力系统潮流计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910684342.2A CN110460059B (zh) | 2019-07-26 | 2019-07-26 | 一种电力系统潮流计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110460059A true CN110460059A (zh) | 2019-11-15 |
CN110460059B CN110460059B (zh) | 2023-02-03 |
Family
ID=68483566
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910684342.2A Active CN110460059B (zh) | 2019-07-26 | 2019-07-26 | 一种电力系统潮流计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110460059B (zh) |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030040846A1 (en) * | 2001-05-21 | 2003-02-27 | Christian Rehtanz | Stability prediction for an electric power network |
WO2003090328A1 (en) * | 2002-04-22 | 2003-10-30 | The Tokyo Electric Power Company, Incorporated | Method and system for on-line dynamical screening of electric power system |
US20060229767A1 (en) * | 2005-04-08 | 2006-10-12 | Chang Gung University | Method for calculating power flow solution of a power transmission network that includes unified power flow controllers |
US20130289905A1 (en) * | 2012-04-23 | 2013-10-31 | Institute Of Nuclear Energy Research Atomic Energy Council, Executive Yuan | Distribution power flow analysis system and method |
CN103632046A (zh) * | 2013-11-20 | 2014-03-12 | 国网河南省电力公司南阳供电公司 | 一种电网潮流计算方法 |
CN104052063A (zh) * | 2014-07-03 | 2014-09-17 | 国家电网公司 | 一种基于计及静态负荷特性的节点电压稳定指标的电网电压稳定监测方法 |
CN106066918A (zh) * | 2016-06-06 | 2016-11-02 | 国网山东省电力公司经济技术研究院 | 基于同伦的含分布式电源和非线性负荷的短路电流计算方法 |
CN106295160A (zh) * | 2016-08-04 | 2017-01-04 | 山东大学 | 交直流互联电网戴维南等值参数在线计算方法 |
-
2019
- 2019-07-26 CN CN201910684342.2A patent/CN110460059B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20030040846A1 (en) * | 2001-05-21 | 2003-02-27 | Christian Rehtanz | Stability prediction for an electric power network |
WO2003090328A1 (en) * | 2002-04-22 | 2003-10-30 | The Tokyo Electric Power Company, Incorporated | Method and system for on-line dynamical screening of electric power system |
US20060229767A1 (en) * | 2005-04-08 | 2006-10-12 | Chang Gung University | Method for calculating power flow solution of a power transmission network that includes unified power flow controllers |
US20130289905A1 (en) * | 2012-04-23 | 2013-10-31 | Institute Of Nuclear Energy Research Atomic Energy Council, Executive Yuan | Distribution power flow analysis system and method |
CN103632046A (zh) * | 2013-11-20 | 2014-03-12 | 国网河南省电力公司南阳供电公司 | 一种电网潮流计算方法 |
CN104052063A (zh) * | 2014-07-03 | 2014-09-17 | 国家电网公司 | 一种基于计及静态负荷特性的节点电压稳定指标的电网电压稳定监测方法 |
CN106066918A (zh) * | 2016-06-06 | 2016-11-02 | 国网山东省电力公司经济技术研究院 | 基于同伦的含分布式电源和非线性负荷的短路电流计算方法 |
CN106295160A (zh) * | 2016-08-04 | 2017-01-04 | 山东大学 | 交直流互联电网戴维南等值参数在线计算方法 |
Non-Patent Citations (3)
Title |
---|
张泽中: "考虑负荷模型的潮流计算在低电压治理中的应用研究", 《中国优秀硕士学位论文全文数据库工程科技II辑》 * |
王淳等: "负荷电压静特性模型对潮流计算的影响分析", 《南昌大学学报(工科版)》 * |
邱丽萍等: "综合负荷模型对大区互联电网稳定特性的影响", 《电网技术》 * |
Also Published As
Publication number | Publication date |
---|---|
CN110460059B (zh) | 2023-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109861303A (zh) | 一种含母线电压补偿和负荷功率动态分配的直流微电网协调控制方法 | |
KR101132107B1 (ko) | 분산전원이 연계된 전력계통의 전압/무효전력 제어 시스템 및 이를 위한 방법 | |
US20060282239A1 (en) | Method of setting-up steady state model of VSC-based multi-terminal HVDC transmission system | |
CN104778367B (zh) | 基于单一状态断面的广域戴维南等值参数在线计算方法 | |
CN107492888B (zh) | 一种直流受端交流系统的电压薄弱区域评估方法 | |
CN106066918B (zh) | 一种含分布式电源和非线性负荷的短路电流计算方法 | |
CN109802381A (zh) | 一种基于模糊控制的直流微电网多源动态协调控制方法 | |
CN104701858B (zh) | 计及分区动态无功储备的无功电压控制方法 | |
CN110535180B (zh) | 基于灵活资源与传统调压设备协调控制配电网电压的方法 | |
CN106374498A (zh) | 一种考虑二次电压频率控制的微电网潮流计算方法 | |
CN108493952A (zh) | 一种基于模糊自适应补偿的交流微电网无功均分控制方法 | |
CN101527455A (zh) | 基于潮流模块交替迭代的互联电网分布式潮流计算方法 | |
CN109586301A (zh) | 基于同步相量量测等值的智能配电网就地电压控制方法 | |
CN110165667A (zh) | 计及静态电压安全约束的输配协同无功优化方法及系统 | |
CN106253305A (zh) | 一种针对分布式电源的全分散式孤岛无功控制方法 | |
CN107994562A (zh) | 一种计及负载特性的柴储型车载电源系统的稳定性设计方法 | |
CN108988401B (zh) | 多通道交直流混联断面输电能力及其约束条件的评估方法 | |
CN107302224B (zh) | 一种基于内点法的多端直流输电系统换流站控制方法 | |
Zhao et al. | Reactive power optimization considering dynamic reactive power reserves | |
CN110460059A (zh) | 一种电力系统潮流计算方法 | |
CN109830964B (zh) | 一种微电网动态功率平衡方法及系统 | |
CN106356867A (zh) | 一种节点无功电压灵敏度的确定方法 | |
CN109802448A (zh) | 一种可再生能源最大消纳容量分析计算方法 | |
Wu et al. | Research on improvement of line loss algorithm based on three-phase unbalance degree | |
Ciric et al. | Integration of the dispersed generators in the distribution management system |
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 |