CN111478362B - 一种光伏电站和电网的一体化仿真方法、系统和存储介质 - Google Patents

一种光伏电站和电网的一体化仿真方法、系统和存储介质 Download PDF

Info

Publication number
CN111478362B
CN111478362B CN202010277195.XA CN202010277195A CN111478362B CN 111478362 B CN111478362 B CN 111478362B CN 202010277195 A CN202010277195 A CN 202010277195A CN 111478362 B CN111478362 B CN 111478362B
Authority
CN
China
Prior art keywords
node
grid
photovoltaic power
power
power station
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
CN202010277195.XA
Other languages
English (en)
Other versions
CN111478362A (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.)
Beijing Kedong Electric Power Control System Co Ltd
Training Center of State Grid Hebei Electric Power Co Ltd
Original Assignee
Beijing Kedong Electric Power Control System Co Ltd
Training Center of State Grid Hebei Electric Power Co Ltd
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 Beijing Kedong Electric Power Control System Co Ltd, Training Center of State Grid Hebei Electric Power Co Ltd filed Critical Beijing Kedong Electric Power Control System Co Ltd
Priority to CN202010277195.XA priority Critical patent/CN111478362B/zh
Publication of CN111478362A publication Critical patent/CN111478362A/zh
Application granted granted Critical
Publication of CN111478362B publication Critical patent/CN111478362B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

本发明公开了一种光伏电站和电网的一体化仿真方法、系统和存储介质,方法包括将并网光伏电站的逆变器输出节点、滤波器节点以及并网点作为新增PQ注入节点;根据原来系统PQ节点和PV节点的注入有功功率偏差、新增PQ节点注入的有功功率无功功率、原系统PQ节点的电压幅值修正量以及原系统中除平衡节点外其他节点的相角修正量,构建状态变量值的修正方程并利用PQ分解法求解获得状态变量值的修正量;最终获得各个光伏电站并入电网后内部运行状态数据。本发明将光伏电站的设备融入到电网仿真之中,考虑了电网电压和相角对逆变器、滤波器和变压器的有功和无功的影响,由此可以获得准确的变电站内和电网的功率分布和运行状态。

Description

一种光伏电站和电网的一体化仿真方法、系统和存储介质
技术领域
本发明属于电力系统光伏并网仿真领域,特别涉及一种光伏电站和电网的一体化仿真方法、系统和存储介质。
背景技术
目前的计算机仿真计算分析含光伏并网的电网仿真时,为了简化分析,将并网点作为PQ节点,取近似估计值为并网节点的有功和无功注入,进行全网潮流计算后得出并网点的电压大小和相角,然后回带计算变压器和滤波器的功率损耗和逆变器的输出功率。由于光伏板输出的直流功率在逆变器上存在一定的损耗,且逆变器输出有功功率和无功功率在滤波器和变压器上都存在损耗,损耗的大小与滤波器、变压器的具体结构和并网点电压都有关系,且状态变量并网点电压未知,所以并网点的注入有功功率和无功功率无法准确获得。所以目前在处理光伏并网是将并网点作为PQ节点给出近似估计的并网点的注入有功功率和无功功率进行全网潮流计算的方法存在一定的误差。这种方法实质是把光伏电站和电网分解开来计算,忽略了电网电压和相角对逆变器、滤波器和变压器的有功和无功的影响。
发明内容
本发明旨在针对目前在处理光伏并网是将并网点作为PQ节点给出近似估计的并网点输出有功和无功功率进行全网潮流计算的方法存在一定的误差,这种方法实质是把光伏电站和电网分解开来计算,忽略了电网电压和相角对逆变器、滤波器和变压器的有功和无功的影响,本发明提出一种光伏电站和电网一体化仿真方法。
为实现上述技术目的,本发明采用以下技术方案。
本发明提供一种光伏电站和电网的一体化仿真方法,包括以下步骤:
一种光伏电站和电网的一体化仿真方法,所述方法包括以下步骤:
将并网光伏电站的逆变器输出节点、滤波器节点以及并网点作为新增PQ节点,获取新增PQ节点注入的有功功率偏差、新增PQ节点注入的无功功率偏差、原系统m个PQ节点的电压幅值修正量;
根据原来系统PQ节点和PV节点的注入有功功率偏差、新增PQ节点注入的有功功率偏差、新增PQ节点注入的无功功率偏差、原系统m个PQ节点的电压幅值修正量以及原系统中除平衡节点外其他节点的相角修正量,构建状态变量值的修正方程;
利用PQ分解法求解所构建的状态变量值的修正方程获得状态变量值的修正量;
根据状态变量值的修正量获得光伏电站并入电网后内部运行状态数据。
进一步地,光伏电站并入电网后各新增节点注入有功功率偏差ΔPi的表达式如下:
Figure BDA0002445221380000021
其中ΔPi为光伏电站并入电网后各新增PQ节点注入有功功率偏差即节点i注入有功功率偏差,Ui为节点i的电压幅值,Uj为节点j的电压幅值,Gij为节点i和节点j之间的电导,θij为节点i和节点j之间的相角差,Bij为节点i和节点j之间的电纳,PG,i为节点i设备注入有功功率,PD,i为负荷消耗有功功率。
再进一步地,当节点i为逆变器接入节点时,PG,i=Pinv,i=ηi·Ph,i,PD,i取0,其中,i为光伏电站编号,i=1,2,…,s,s为接入的光伏电站总数,Pinv,i为第i个光伏电站中逆变器的输出有功功率,ηi为第i各光伏电站中逆变器的效率,Nph,i为第i个光伏电站中逆变器的输入直流功率。
进一步地,光伏电站并入电网后各新增节点注入无功功率偏差的计算公式如下:
Figure BDA0002445221380000031
其中ΔQi为各新增PQ节点注入无功功率偏差,Ui为节点i的电压幅值,Uj为节点j的电压幅值,Gij为节点i和节点j之间的电导,θij为节点i和节点j之间的相角差,Bij为节点i和节点j之间的电纳,QG,i为节点i设备发出无功功率;QD,i为负荷消耗无功功率。
进一步地,当节点i为逆变器接入节点时,
Figure BDA0002445221380000032
QD,i取0,其中i为光伏电站编号,i=1,2,…,s,s为接入的光伏电站总数,Qinv,i为第i个光伏电站中逆变器的输出无功功率,
Figure BDA0002445221380000033
为第i个光伏电站的并网点电压,
Figure BDA0002445221380000034
为第i个光伏电站的并网点电流相量。
构建的状态变量值的修正方程如下:
Figure BDA0002445221380000041
Figure BDA0002445221380000042
为原来系统PQ节点和PV节点的注入有功功率偏差,n-1维;
Figure BDA0002445221380000043
为第1个光伏电站中并网新增节点注入有功功率偏差,k1维;
Figure BDA0002445221380000044
为第s个光伏电站并网新增节点注入有功功率偏差,ks维;
Figure BDA0002445221380000045
为原系统m个PQ节点的注入无功功率偏差,m维;
Figure BDA0002445221380000046
为第1个光伏电站中并网新增节点注入无功功率偏差,k1维;
Figure BDA0002445221380000051
为第s个光伏电站并网后新增节点注入无功功率偏差,ks维;
Figure BDA0002445221380000052
为原系统m个PQ节点的电压幅值修正量,m维;
Figure BDA0002445221380000053
为光伏电站1接入电网新增的节点电压幅值修正量,k1维;
Figure BDA0002445221380000054
为光伏电站s接入电网新增的节点电压幅值修正量,ks维;
Figure BDA0002445221380000055
为原系统中除平衡节点外其他节点的相角修正量,n-1维;
Figure BDA0002445221380000056
为光伏电站1接入电网新新增节点电压相角修正量,k1维;
Figure BDA0002445221380000057
为光伏电站s接入电网新增的节点电压相角修正量,ks维。
进一步地,根据状态变量值的修正量获得各个光伏电站并入电网后内部运行状态数据的方法如下:
Figure BDA0002445221380000061
其中向量Um为原系统m个PQ节点的电压幅值,向量U1为光伏电站1接入电网新增的所有节点电压幅值、向量Us为伏电站s接入电网新增的所有节点电压幅值、向量θn-1为原系统中除平衡节点外其他节点的相角,向量θ1为光伏电站1接入电网所有新增节点电压相角,向量θs为光伏电站s接入电网所有新增节点电压相角,
Figure BDA0002445221380000062
Figure BDA0002445221380000063
分别为第l和l+1次修正后的状态变量值,
Figure BDA0002445221380000071
为第l次状态变量值的修正值。
另一方面,本发明提供一种光伏电站和电网的一体化仿真系统,所述系统包括:新电网构建模块,用于将并网光伏电站的逆变器输出节点、滤波器节点以及并网点作为新增PQ节点;
参数获取模块,用于获取新增PQ节点注入的有功功率偏差、新增PQ节点注入的无功功率偏差和原系统m个PQ节点的电压幅值修正量;
模型建立和求解模块,用于根据原来系统PQ节点和PV节点的注入有功功率偏差、新增PQ节点注入的有功功率偏差和新增PQ节点注入的无功功率偏差、原系统m个PQ节点的电压幅值修正量以及原系统中除平衡节点外其他节点的相角修正量构建状态变量值的修正方程,利用PQ分解法求解所构建的修正方程获得状态变量值的修正量;根据确定的状态变量值的修正量获得各个光伏电站并入电网后内部运行状态数据。
本发明提供一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如以上技术方案提供的方法的步骤。
有益技术效果:
本发明通过增加逆变器节点、滤波器节点和变压器节点等节点,扩大电网仿真涉及节点及支路的范围,将光伏电站的设备融入到电网仿真之中,作为电网的一部分,将光伏电站和电网作为一个整体去分析计算,考虑了电网电压和相角对逆变器、滤波器和变压器的有功和无功的影响,由此可以获得准确的变电站内和电网的功率分布和运行状态。
附图说明
图1是本发明具体实施例有s座光伏电站接入电网的示意图;
图2是本发明具体实施例增加光伏电站设备后组成的新电网结构示意图;
图3是本发明具体实施例的结构示意图;
图4是本发明具体实施例的计算结果图。
具体实施方式
为更好地理解本发明,以下给出名词解释:
PQ节点:这类节点的有功功率P和无功功率Q是给定的,节点电压和相位(V,δ)是待求量。
PV节点:这类节点的有功功率P和电压幅值V是给定的,节点的无功功率Q和电压相位δ是待求量,这类节点必须有足够的可调无功容量,用以维持给定的电压幅值,因而又称之为“电压控制节点”。
平衡节点的电压幅值和相位是给定的,而其注入有功功率和无功功率是待求量。
本发明提出一种光伏电站和电网一体化仿真方法,通过增加逆变器节点、滤波器节点和变压器节点等节点,扩大电网仿真涉及节点及支路的范围,将光伏电站的设备融入到电网仿真之中,作为电网的一部分,将光伏电站和电网作为一个整体去分析计算电网潮流分布和运行状态。由于电网中的变压器可以作为线性元件,常用的电网LCL滤波器也为线性元件,所以滤波器和变压器完全可以作为电网的一部分,与原电网组成新的线性网络。逆变器输出节点作为PQ注入节点处理,光伏出力Pph和并网电流
Figure BDA0002445221380000091
作为控制变量,可得逆变器输出功率
Figure BDA0002445221380000092
其他像滤波器和变压器新增的节点可以视为联络节点,注入功率为0。
具体实施例提供的一种光伏电站和电网一体化仿真方法包括以下步骤:
1.确定光伏电站并入的大电网的网络结构、节点类型及参数,假定未接入光伏前原系统中有n个节点,其中m个PQ节点,1个平衡接点,其余为PV节点;所述光伏电站包括光伏板、逆变器、滤波器和变压器,所述光伏板依次连接逆变器、滤波器和变压器后连接到电网;
2.确定光伏电站数量s和他们的并网点1至并网点s;
3.确定各个光伏电站的逆变器和滤波器参数;
4.原系统中电网有n个节点,将光伏电站的逆变器输出节点和滤波器节点也作为电网的一部分节点,若接入第i个光伏电站后在并网侧新增ki个节点,则光伏接入系统后,新增
Figure BDA0002445221380000093
个节点,新电网总共有
Figure BDA0002445221380000094
个节点。
5.接入第i个光伏电站后在并网侧新增ki个节点,光伏接入系统后,新增
Figure BDA0002445221380000095
个节点,光伏接入后新增节点全部作为PQ节点处理。但是由于逆变器输出节点的无功Qinv,i无法确定,所以逆变器节点只能作为P节点处理。
6.此时系统共有
Figure BDA0002445221380000101
个节点的注入有功功率已知,
Figure BDA0002445221380000102
节点注入无功功率已知,可列如下方程:
Figure BDA0002445221380000103
Figure BDA0002445221380000104
个。
Figure BDA0002445221380000105
Figure BDA0002445221380000106
其中,ΔPi为光伏电站并入电网后各新增节点注入有功功率偏差,ΔQi为光伏电站并入电网后各新增节点注入无功功率偏差;PG,i,QG,i为设备发出的有功功率和无功功率;PD,i,QD,i为负荷消耗功率。当节点i为逆变器接入节点时,PG,i=Pinv,i=ηi·Pph,i,PD,i=QD,i=0,其中Pph,i为光伏板传输到逆变器的直流功率,Qinv,i逆变器输出的交流无功功率。
根据第i个光伏电站的光照和温度可以得出光伏板传输到逆变器的直流功率Pph,i,这里的Pph,i为已知量。
求出逆变器输出的交流有功功率:Pinv,i=ηi·Pph,i i=1,2,…,s,其中,Pinv,i为第i个逆变器的输出有功功率,ηi为第i个逆变器的效率,Pph,i为第i个逆变器的输入直流功率。逆变器输出的交流无功功率Qinv,i目前无法求出具体数值,它受并网点电压和电流
Figure BDA0002445221380000107
的影响,
Figure BDA0002445221380000108
其中,Qinv,i为第i个逆变器的输出的无功功率,
Figure BDA0002445221380000109
分别为第i个光伏电站的并网点电压和电流相量。此时的ΔPi和ΔQi共有
Figure BDA00024452213800001010
个方程组,但是新电网有
Figure BDA00024452213800001011
个节点,有
Figure BDA0002445221380000111
个待求变量,由于约束不够此时无法进行求解,必须寻找新的约束条件。
7.在光伏并网时,通过对逆变器的控制,使得并网点注入电流
Figure BDA0002445221380000112
跟踪并网点电压
Figure BDA0002445221380000113
使得
Figure BDA0002445221380000114
的大小和相位与
Figure BDA0002445221380000115
形成固定的关系,可以表示为
Figure BDA0002445221380000116
此时可以根据
Figure BDA0002445221380000117
Figure BDA0002445221380000118
可以求得
Figure BDA0002445221380000119
和它的相邻节点电压之间的关系。
Figure BDA00024452213800001110
其中
Figure BDA00024452213800001111
表示与并网点
Figure BDA00024452213800001112
的相邻的光伏电站侧节点,z为相邻节点之间的阻抗。对于每个光伏并网节点均可列出
Figure BDA00024452213800001113
可以列出s个方程。加上ΔPi和ΔQi
Figure BDA00024452213800001114
个方程组,总共可以列出
Figure BDA00024452213800001115
个方程组,与未知数
Figure BDA00024452213800001116
相等,故可以求解出所有未知量。
8.采用牛顿拉夫逊法求解此方程组,状态变量值的修正方程如下:
Figure BDA00024452213800001117
Figure BDA0002445221380000121
为原来系统PQ节点和PV节点的注入有功功率偏差,n-1维。
Figure BDA0002445221380000122
为第1个光伏电站中并网逆变器、滤波器等新增节点注入有功功率偏差,k1维。
Figure BDA0002445221380000123
为第s个光伏电站中并网逆变器、滤波器等新增节点注入有功功率偏差,ks维。
Figure BDA0002445221380000124
为原系统m个PQ节点的注入无功功率偏差,m维。
Figure BDA0002445221380000125
为第1个光伏电站中并网逆变器、滤波器等新增节点注入无功功率偏差,k1维。
Figure BDA0002445221380000126
为第s个光伏电站并网后新增逆变器、滤波器等新增节点注入无功功率偏差,ks维。
Figure BDA0002445221380000127
为原系统m个PQ节点的电压幅值修正量,m维。
Figure BDA0002445221380000128
为光伏电站1接入电网新增的节点电压幅值修正量,k1维。
Figure BDA0002445221380000131
为光伏电站s接入电网新增的节点电压幅值修正量,ks维。
Figure BDA0002445221380000132
为原系统中除平衡节点外其他节点的相角修正量,n-1维
Figure BDA0002445221380000133
为光伏电站1接入电网新新增节点电压相角修正量,k1维。
Figure BDA0002445221380000134
为光伏电站s接入电网新增的节点电压相角修正量,ks维。
获得各个光伏电站并入电网后内部运行状态数据的方法如下:
Figure BDA0002445221380000135
Figure BDA0002445221380000141
为第l次状态变量值的修正量
Figure BDA0002445221380000142
Figure BDA0002445221380000143
分别为第l和l+1次修正后的状态变量值收敛条件
Figure BDA0002445221380000144
当中的最大值小于ε时,求得的
Figure BDA0002445221380000145
即为最终解,可以据此计算各种关心的依从变量。
上述过程为多个光伏电站同时接入电网的一般化系统分析方法,下面给出一个简单的算例(一个光伏电站接入三节点电网系统)。
实施例如图3所示,一座光伏内部含有滤波器和隔离变压器的光伏电站经节点1接入电网。电网包含三个节点,节点1为PQ母线,负荷功率为P1+jQ1=2+j1;节点2是PV母线,有功注入P2=0.5,节点电压给定值为1.01;节点3是Vθ母线,电压是V3=1.0∠0°。控制变量:逆变器输入有功P6=4.07,并网电流I41=4.16,相位比
Figure BDA0002445221380000146
落后3°。滤波器参数:X1=0.001,X2=0.0015,Y1=0.0005+j0.03,变压器参数:Zm=30+j160,RT=0.01,XT=0.25。
光伏接入电网后原3节点系统扩充到6节点,4,5节点做常规PQ节点出力,他们的注入功率均为0,6节点较为特殊,光伏逆变器输出的交流功率P6=4.07已知,但它的无功注入和电压均为代求变量,此时进行潮流计算时,节点6缺少一个未知量,相应的ΔQi方程缺少一个,无法求解。但是控制变量注入电流
Figure BDA0002445221380000151
大小已知,它和节点1电压相量
Figure BDA0002445221380000152
的相位关系已知,所以可以利用此条件求出节点4和节点1之间的关系
Figure BDA0002445221380000153
此时方程数与未知数个数相等,可以进行求解。其中,
Figure BDA0002445221380000154
表示节点i的电压向量,
Figure BDA0002445221380000155
表示节点i流向节点j的电流,
Figure BDA0002445221380000156
表示节点i对地电流。计算结果如图4所示。实施例:一种光伏电站和电网的一体化仿真系统,所述光伏电站包括光伏板、逆变器、滤波器和变压器,所述光伏板依次连接逆变器、滤波器和变压器后连接到电网;所述系统包括:
新电网构建模块,用于将并网光伏电站的逆变器输出节点、滤波器节点以及并网点作为新增PQ节点;
参数获取模块,用于获取新增PQ节点注入的有功功率偏差、新增PQ节点注入的无功功率偏差和原系统PQ节点的电压幅值修正量;
模型建立和求解模块,用于根据原来系统PQ节点和PV节点的注入有功功率偏差、新增PQ节点注入的有功功率偏差和新增PQ节点注入的无功功率偏差、原系统PQ节点的电压幅值修正量以及原系统中除平衡节点外其他节点的相角修正量构建状态变量值的修正方程,利用PQ分解法求解所构建的修正方程获得状态变量值的修正量;根据确定的状态变量值的修正量获得各个光伏电站并入电网后内部运行状态数据。
需要说明的是,本实施例各个模块的实现方法如上述技术方案所述,在此不再赘述。
本发明通过增加逆变器节点、滤波器节点和变压器节点等节点,扩大电网仿真涉及节点及支路的范围,将光伏电站的设备融入到电网仿真之中,作为电网的一部分,将光伏电站和电网作为一个整体去分析计算,由此可以获得准确的变电站内和电网的功率分布和运行状态。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。

Claims (10)

1.一种光伏电站和电网的一体化仿真方法,其特征在于,所述方法包括以下步骤:
将并网光伏电站的逆变器输出节点、滤波器节点以及并网点作为新增PQ节点,获取新增PQ节点注入的有功功率偏差、新增PQ节点注入的无功功率偏差和原系统PQ节点的电压幅值修正量;
根据原系统PQ节点和PV节点的注入有功功率偏差、新增PQ节点注入的有功功率偏差、新增PQ节点注入的无功功率偏差、原系统PQ节点的电压幅值修正量以及原系统中除平衡节点外其他节点的相角修正量,构建状态变量值的修正方程;
利用PQ分解法求解所构建的状态变量值的修正方程,获得状态变量值的修正量;
根据状态变量值的修正量获得光伏电站并入电网后内部运行状态数据,所述状态变量值表示为:
Figure FDA0003153821040000011
其中向量Um为原系统m个PQ节点的电压幅值,向量U1为光伏电站1接入电网新增的所有节点电压幅值、向量Us为光伏电站s接入电网新增的所有节点电压幅值、向量θn-1为原系统中除平衡节点外其他节点的相角,向量θ1为光伏电站1接入电网所有新增节点电压相角,向量θs为光伏电站s接入电网所有新增节点电压相角。
2.根据权利要求1所述的一种光伏电站和电网的一体化仿真方法,其特征在于,光伏电站并入电网后各新增PQ节点注入有功功率偏差的表达式如下:
Figure FDA0003153821040000021
其中ΔPi为光伏电站并入电网后各新增PQ节点注入有功功率偏差,Ui为节点i的电压幅值,Uj为节点j的电压幅值,Gij为节点i和节点j之间的电导,θij为节点i和节点j之间的相角差,Bij为节点i和节点j之间的电纳,PG,i为节点i设备注入有功功率,PD,i为负荷消耗有功功率。
3.根据权利要求2所述的一种光伏电站和电网的一体化仿真方法,其特征在于,若节点i为逆变器接入节点,PG,i=Pinv,i=ηi·Pph,i,PD,i取0,其中i为光伏电站编号,i=1,2,…,s,s为接入的光伏电站总数,Pinv,i为第i个光伏电站中逆变器的输出有功功率,ηi为第i各光伏电站中逆变器的效率,Pph,i为第i个光伏电站中逆变器的输入直流功率。
4.根据权利要求1所述的一种光伏电站和电网的一体化仿真方法,其特征在于,光伏电站并入电网后各新增PQ节点注入无功功率偏差的计算公式如下:
Figure FDA0003153821040000031
其中ΔQi为光伏电站并入电网后各新增PQ节点注入无功功率偏差,Ui为节点i的电压幅值,Uj为节点j的电压幅值,Gij为节点i和节点j之间的电导,θij为节点i和节点j之间的相角差,Bij为节点i和节点j之间的电纳,QG,i为节点i设备发出无功功率;QD,i为负荷消耗无功功率。
5.根据权利要求4所述的一种光伏电站和电网的一体化仿真方法,其特征在于,若节点i为逆变器接入节点,
Figure FDA0003153821040000032
QD,i取0,其中i为光伏电站编号,i=1,2,…,s,s为接入的光伏电站总数,Qinv,i为第i个光伏电站中逆变器的输出无功功率,
Figure FDA0003153821040000033
为第i个光伏电站的并网点电压,
Figure FDA0003153821040000034
为第i个光伏电站的并网点电流相量。
6.根据权利要求5所述的一种光伏电站和电网的一体化仿真方法,其特征在于,控制并网点注入电流
Figure FDA0003153821040000035
跟踪并网点电压
Figure FDA0003153821040000036
使得
Figure FDA0003153821040000037
的大小和相位与
Figure FDA0003153821040000038
形成固定的关系,表示为
Figure FDA0003153821040000039
Yg,i为并网点电压和并网点电流相量的关系系数。
7.根据权利要求1所述的一种光伏电站和电网的一体化仿真方法,其特征在于,构建的状态变量值的修正方程如下:
Figure FDA0003153821040000041
Figure FDA0003153821040000042
为原来系统PQ节点和PV节点的注入有功功率偏差,n-1维;
Figure FDA0003153821040000043
为第1个光伏电站中并网新增节点注入有功功率偏差,k1维;
Figure FDA0003153821040000044
为第s个光伏电站并网新增节点注入有功功率偏差,ks维;
Figure FDA0003153821040000045
为原系统m个PQ节点的注入无功功率偏差,m维;
Figure FDA0003153821040000046
为第1个光伏电站中并网新增节点注入无功功率偏差,k1维;
Figure FDA0003153821040000051
为第s个光伏电站并网后新增节点注入无功功率偏差,ks维;
Figure FDA0003153821040000052
为原系统m个PQ节点的电压幅值修正量,m维;
Figure FDA0003153821040000053
为光伏电站1接入电网新增的节点电压幅值修正量,k1维;
Figure FDA0003153821040000054
为光伏电站s接入电网新增的节点电压幅值修正量,ks维;
Figure FDA0003153821040000055
为原系统中除平衡节点外其他节点的相角修正量,n-1维;
Figure FDA0003153821040000056
为光伏电站1接入电网新新增节点电压相角修正量,k1维;
Figure FDA0003153821040000057
为光伏电站s接入电网新增的节点电压相角修正量,ks维。
8.根据权利要求7所述的一种光伏电站和电网的一体化仿真方法,其特征在于,根据状态变量值的修正量获得各个光伏电站并入电网后内部运行状态数据的方法如下:
Figure FDA0003153821040000061
Figure FDA0003153821040000062
Figure FDA0003153821040000063
分别为第l和l+1次修正后的状态变量值,
Figure FDA0003153821040000064
为第l次状态变量值的修正量。
9.一种光伏电站和电网的一体化仿真系统,其特征在于,所述系统包括:
参数获取模块:将并网光伏电站的逆变器输出节点、滤波器节点以及并网点作为新增PQ节点,获取新增PQ节点注入的有功功率偏差、新增PQ节点注入的无功功率偏差和原系统PQ节点的电压幅值修正量;
修正方程构建模块:根据原系统PQ节点和PV节点的注入有功功率偏差、新增PQ节点注入的有功功率偏差、新增PQ节点注入的无功功率偏差、原系统PQ节点的电压幅值修正量以及原系统中除平衡节点外其他节点的相角修正量,构建状态变量值的修正方程;
修正量获取模块:利用PQ分解法求解所构建的状态变量值的修正方程,获得状态变量值的修正量;
运行状态数据获取模块:根据状态变量值的修正量获得光伏电站并入电网后内部运行状态数据;所述状态变量值表示为:
Figure FDA0003153821040000071
其中向量Um为原系统m个PQ节点的电压幅值,向量U1为光伏电站1接入电网新增的所有节点电压幅值、向量Us为伏电站s接入电网新增的所有节点电压幅值、向量θn-1为原系统中除平衡节点外其他节点的相角,向量θ1为光伏电站1接入电网所有新增节点电压相角,向量θs为光伏电站s接入电网所有新增节点电压相角。
10.一种计算机可读存储介质,所述计算机可读存储介质存储有计算机程序,其特征在于,所述计算机程序被处理器执行时实现如权利要求1~8任意一项权利要求所述方法的步骤。
CN202010277195.XA 2020-04-10 2020-04-10 一种光伏电站和电网的一体化仿真方法、系统和存储介质 Active CN111478362B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010277195.XA CN111478362B (zh) 2020-04-10 2020-04-10 一种光伏电站和电网的一体化仿真方法、系统和存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010277195.XA CN111478362B (zh) 2020-04-10 2020-04-10 一种光伏电站和电网的一体化仿真方法、系统和存储介质

Publications (2)

Publication Number Publication Date
CN111478362A CN111478362A (zh) 2020-07-31
CN111478362B true CN111478362B (zh) 2021-09-17

Family

ID=71751522

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010277195.XA Active CN111478362B (zh) 2020-04-10 2020-04-10 一种光伏电站和电网的一体化仿真方法、系统和存储介质

Country Status (1)

Country Link
CN (1) CN111478362B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103199524A (zh) * 2013-03-12 2013-07-10 武汉大学 一种适应多种分布式电源接入的潮流计算方法
CN106410835A (zh) * 2016-10-12 2017-02-15 北京科东电力控制系统有限责任公司 一种多端柔性直流输电系统的电网仿真实现方法
CN107171329A (zh) * 2017-07-11 2017-09-15 国网江苏省电力公司宿迁供电公司 一种含新能源并网的配电网潮流计算方法
CN108808681A (zh) * 2017-05-02 2018-11-13 南京理工大学 基于混合注入模型的光伏并网潮流计算方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103199524A (zh) * 2013-03-12 2013-07-10 武汉大学 一种适应多种分布式电源接入的潮流计算方法
CN106410835A (zh) * 2016-10-12 2017-02-15 北京科东电力控制系统有限责任公司 一种多端柔性直流输电系统的电网仿真实现方法
CN108808681A (zh) * 2017-05-02 2018-11-13 南京理工大学 基于混合注入模型的光伏并网潮流计算方法
CN107171329A (zh) * 2017-07-11 2017-09-15 国网江苏省电力公司宿迁供电公司 一种含新能源并网的配电网潮流计算方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
A new method of detecting PV Grid-connected Inverter Islanding based on the frequency variation;Zhen-yue Wang等;《2011 International Conference on Materials for Renewable Energy & Environment》;20110627;全文 *
基于混合注入模型的光伏并网潮流计算研究;王立国等;《现代电力》;20140228;第31卷(第1期);第45-51页 *
潮流计算和机电暂态仿真中风光储联合发电系统的实用等值方法;王皓怀等;《中国电机工程学报》;20120105;第32卷(第1期);第1-8页 *

Also Published As

Publication number Publication date
CN111478362A (zh) 2020-07-31

Similar Documents

Publication Publication Date Title
Bernstein et al. Linear power-flow models in multiphase distribution networks
CN109149620B (zh) 一种自储能多端柔直系统控制方法及系统
Lai et al. Distributed power control for DERs based on networked multiagent systems with communication delays
CN105514971A (zh) 一种适用于各种运行模式微电网的潮流计算方法
CN103077480B (zh) 一种电力系统的安全校核方法
CN104794541B (zh) 一种基于模拟退火和锥优化的配电网运行优化方法
CN104156542B (zh) 一种基于隐式投影的有源配电系统稳定性仿真方法
CN108054757A (zh) 一种内嵌无功和电压的n-1闭环安全校核方法
Abouloifa et al. Output‐feedback nonlinear adaptive control strategy of the single‐phase grid‐connected photovoltaic system
CN112467748A (zh) 三相不平衡主动配电网双时标分布式电压控制方法及系统
Lian et al. Steady‐state solutions of a voltage source converter with dq‐frame controllers by means of the time‐domain method
CN109638871B (zh) 考虑风电接入的大规模交直流混联系统主网划分方法
CN110137948B (zh) 一种双馈风电场次同步振荡风险的评估方法及装置
CN113890039B (zh) 一种多端柔性直流配电网潮流调度优化方法
CN111327048A (zh) 一种含三端snop的配电网鲁棒运行优化方法
CN102780220A (zh) 含有pv恒定型分布式电源的配电网潮流计算方法
CN111478362B (zh) 一种光伏电站和电网的一体化仿真方法、系统和存储介质
CN113206517A (zh) 孤岛微电网频率和电压恢复控制方法、装置、设备和介质
CN112531756A (zh) 用于储能系统电量平衡的分布式控制方法及系统、设备
CN104377721A (zh) 一种电网电压不平衡时vsc-hvdc优化控制方法
CN112234598A (zh) 一种电磁暂态仿真初始化方法
Nicolini et al. PLL parameters tuning guidelines to increase stability margins in multiple three‐phase converters connected to weak grids
Gupta et al. Coordinated load flow solution for coupled transmission-distribution system incorporating load modeling
CN115133540A (zh) 一种配电网无模型的实时电压控制方法
CN103701119A (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