CN113162072B - 基于光伏电站运行特性和光照强度的变下垂系数控制方法 - Google Patents

基于光伏电站运行特性和光照强度的变下垂系数控制方法 Download PDF

Info

Publication number
CN113162072B
CN113162072B CN202110449595.9A CN202110449595A CN113162072B CN 113162072 B CN113162072 B CN 113162072B CN 202110449595 A CN202110449595 A CN 202110449595A CN 113162072 B CN113162072 B CN 113162072B
Authority
CN
China
Prior art keywords
photovoltaic power
power station
active
reactive
droop
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
CN202110449595.9A
Other languages
English (en)
Other versions
CN113162072A (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.)
State Grid Corp of China SGCC
Wuhan University WHU
Global Energy Interconnection Research Institute
Original Assignee
State Grid Corp of China SGCC
Wuhan University WHU
Global Energy Interconnection Research Institute
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 State Grid Corp of China SGCC, Wuhan University WHU, Global Energy Interconnection Research Institute filed Critical State Grid Corp of China SGCC
Priority to CN202110449595.9A priority Critical patent/CN113162072B/zh
Publication of CN113162072A publication Critical patent/CN113162072A/zh
Application granted granted Critical
Publication of CN113162072B publication Critical patent/CN113162072B/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/24Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/241The oscillation concerning frequency
    • 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/12Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load
    • H02J3/16Circuit arrangements for ac mains or ac distribution networks for adjusting voltage in ac networks by changing a characteristic of the network load by adjustment of reactive power
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/30Reactive power compensation

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Electrical Variables (AREA)

Abstract

本发明提供一种基于光伏电站运行特性和光照强度的变下垂系数控制方法。该方法在光伏发电有功‑频率和无功‑电压下垂控制的基础上,结合太阳辐照度和光伏电站运行特性,对有功‑频率和无功‑电压下垂控制设置随太阳辐照度和光伏电站运行特性变化的下垂系数,从而让光伏发电更加充分的参与到电力系统调频和调压中。本发明通过增大有功‑频率特性曲线和无功‑电压特性曲线在运行点附近的下垂系数,减小频率和电压偏差,提高电力系统的稳定性;通过增大有功或无功调节裕量较大的光伏电站的下垂系数,使裕量较大的光伏电站在调频和调压过程中承担更多任务,各光伏电站充分发挥调频和调压能力,提高光伏电站间协调控制的精度。

Description

基于光伏电站运行特性和光照强度的变下垂系数控制方法
技术领域
本申请涉及光伏发电参与调频和调压领域,具体地涉及一种基于光伏电站运行特性和光照强度的变下垂系数控制方法。
背景技术
光伏发电因储量丰富、成本效益高、清洁环保的优点近年来得到迅速发展,并网容量不断增加。但光伏发电系统作为静止元件,缺乏转动惯量,且一般处于最大功率跟踪模式,不参与系统调频;且光伏发电通过逆变器接入系统,不具备传统同步机的电压支撑能力,不参与电压调节。伴随大规模光伏电力的接入,常规发电机组在电力系统中的占比降低,电力系统惯性减小,一定程度上削弱了电力系统的调频能力及应对功率短缺与频率波动的能力。同时,由于光伏本身不参与并网点电压调节,需要在并网点设置无功补偿器,增加了运行和维护成本,同时没有充分利用光伏逆变器的无功容量。为减小电力系统安全运行的风险,要求光伏机组也要具备调频和调压能力。
现有技术中,较多使用的是给光伏发电增加有功-频率和无功-电压下垂特性,即以系统频率偏差和并网点电压偏差为输入,光伏逆变器自行调节其有功和无功输出,最终减小系统频率偏差和并网点电压偏差,属于有差调节。由于光伏发电通过逆变器并网,其参与调速和调压都是通过控制系统完成,不存在实质的调速和调压装置。因此光伏发电具有可变下垂系数的优点,而现有技术中的定下垂系数控制方法并没有充分发挥这个优势,在系统频率偏差或并网点电压偏差较小时,光伏逆变器增发或减发的有功、无功量较小,不能充分发挥光伏发电的调频和调压潜力,系统频率和并网点电压也均会存在较大偏差,有失稳风险。此外,由于各光伏电池所能输出的最大功率受光照条件的影响较大,受到天气特性和光伏运行特性等差异性影响,各光伏场站的有功和无功裕量也各有不同。常规采用的平均分配策略不能充分发挥不同发电能力光伏参与频率调节和电压调节的潜力。
因此,需要提供一种新型的基于光伏电站运行特性和光照强度的变下垂系数控制方法。
发明内容
为了克服现有技术的不足,本发明的目的是提出一种基于光伏电站运行特性和光照强度的变下垂系数控制方法,能在考虑光伏电站运行特性和天气特性的情况下充分发挥各光伏电站参与频率调节和电压调节的潜力,并提高系统的稳定性和良好暂态响应。
为实现上述目的,本发明所采用的解决方案为:
一种基于光伏电站运行特性和光照强度的变下垂系数控制方法,其包括以下步骤:
步骤1:确定有功-频率控制基本模型和无功-电压控制基本模型,具体包括以下步骤:
步骤1-1:建立有功-频率下垂曲线的基本函数模型;
Figure BDA0003038229980000021
式中:P为光伏电站的有功出力;f为系统频率;fref为频率参考值;P0为光伏电站在系统频率f等于频率参考值fref时的有功出力;e为自然对数;α1为与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数调节因子;α2为与光伏电站的有功调节速率有关的有功-频率下垂曲线的下垂系数调节因子;
步骤1-2:建立无功-电压下垂曲线的基本函数模型;
Figure BDA0003038229980000022
式中:Q为光伏电站的无功出力;U为并网点电压;Uref为电压参考值;Q0为光伏电站在并网点电压U等于电压参考值Uref时的无功出力;β1为与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数调节因子;β2为与光伏电站的无功调节速率有关的无功-电压下垂曲线的下垂系数调节因子;
步骤1-3:将所述步骤1-1的有功-频率下垂曲线的基本函数模型对系统频率f求导数得到光伏电站的有功-频率下垂曲线的下垂系数Kp
Figure BDA0003038229980000023
步骤1-4:将所述步骤1-2的无功-电压下垂曲线的基本函数模型对并网点电压U求导数得到光伏电站的无功-电压下垂曲线的下垂系数Kq
Figure BDA0003038229980000024
步骤2:采集光伏电站的运行特性和太阳辐照度,根据所述光伏电站的运行特性和太阳辐照度,计算光伏电站的有功调节裕量ΔPmargin、光伏电站的有功调节速率ΔPspeed、光伏电站的无功调节裕量ΔQmargin和光伏电站的无功调节速率ΔQspeed
步骤3:根据所述步骤2中获得的光伏电站的有功调节裕量ΔPmargin获得所述步骤1-1中与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数调节因子α1;根据所述步骤2中获得的光伏电站的有功调节速率ΔPspeed获得所述步骤1-1中与光伏电站的有功调节速率有关的有功-频率下垂曲线的下垂系数调节因子α2;根据所述与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数调节因子α1、与光伏电站的有功调节速率有关的有功-频率下垂曲线的下垂系数调节因子α2、系统频率f以及频率参考值fref,得到步骤1-3所述的光伏电站的有功-频率下垂曲线的下垂系数Kp
步骤4:根据所述步骤2中获得的光伏电站的无功调节裕量ΔQmargin获得所述步骤1-2中与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数调节因子β1;根据所述步骤2中获得的光伏电站的无功调节速率ΔQspeed获得所述步骤1-2中与光伏电站的无功调节速率有关的无功-电压下垂曲线的下垂系数调节因子β2;根据所述与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数调节因子β1、与光伏电站的无功调节速率有关的无功-电压下垂曲线的下垂系数调节因子β2、并网点电压U以及电压参考值Uref,得到步骤1-4所述的光伏电站的无功-电压下垂曲线的下垂系数Kq
步骤5:检测系统频率f,根据AGC控制指令、频率参考值fref以及所述步骤3得到的有功-频率下垂曲线的下垂系数Kp获得光伏电站的有功参考值Pref;检测并网点电压U,根据AVC控制指令、电压参考值Uref以及所述步骤4得到的无功-电压下垂曲线的下垂系数Kq获得光伏电站的无功参考值Qref,输出获得的光伏电站的有功参考值Pref和无功参考值Qref给光伏并网逆变器,完成调频和调压任务。
可优选的是,所述步骤3中获得的与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数调节因子为:
Figure BDA0003038229980000031
式中:α1-j为第j个光伏电站与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数的调节因子;
Figure BDA0003038229980000032
为第j个光伏电站的有功调节裕量。
可优选的是,所述步骤3中获得的与光伏电站的有功调节速率有关的有功-频率下垂曲线的下垂系数调节因子为:
Figure BDA0003038229980000041
式中:α2-j为第j个光伏电站与光伏电站的有功调节速率有关的有功-频率下垂曲线的下垂系数的调节因子;ηp为有功调节速率系数;
Figure BDA0003038229980000042
为第j个光伏电站的无功调节速率。
所述有功调节速率系数ηp为传统下垂控制下的最大下垂系数除以标准测试环境下的光伏电站的有功调节速率。
可优选的是,所述步骤4中获得的与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数调节因子为:
Figure BDA0003038229980000043
式中:β1-j为第j个光伏电站与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数的调节因子;
Figure BDA0003038229980000044
为第j个光伏电站的无功调节裕量。
可优选的是,所述步骤4中获得的与光伏电站的无功调节速率有关的无功-电压下垂曲线的下垂系数调节因子为:
Figure BDA0003038229980000045
式中:β2-j为第j个光伏电站与光伏电站的无功调节速率有关的无功-电压下垂曲线的下垂系数的调节因子;ηq为无功调节速率系数,
Figure BDA0003038229980000046
为第j个光伏电站的无功调节速率。
所述无功调节速率系数ηq为传统下垂控制下的最大下垂系数除以标准测试环境下的光伏电站的无功调节速率。
可优选的是,所述步骤5中的AGC控制指令设置为系统频率f等于频率参考值fref时的光伏电站的有功出力P0,所述步骤5中获得的光伏电站的有功参考值Pref为:
Pref=Kp(fref-f)+P0
所述步骤5中的AVC控制指令设置为并网点电压U等于电压参考值Uref时的光伏电站的无功出力Q0,所述步骤5中获得的光伏电站的无功参考值Qref为:
Qref=Kq(Uref-U)+Q0
进一步,所述步骤2中光伏电站的有功调节裕量的具体步骤为:
步骤2-1:将所述步骤2中获取的太阳辐照度根据光伏电池的输出特性,得到光伏电站的最大有功出力:
Figure BDA0003038229980000051
式中:
Figure BDA0003038229980000052
为第j个光伏电站的最大有功出力;n为光伏电站个数;
步骤2-2:根据所述步骤2中获取的光伏电站的运行特性和步骤2-1中得到的各光伏电站的最大有功出力,得到各光伏电站的有功调节裕量:
Figure BDA0003038229980000053
式中:
Figure BDA0003038229980000054
为第j个光伏电站的有功调节裕量;-P0j为第j个光伏电站在系统频率等于频率参考值时的有功出力;
所述步骤2中的光伏电站的有功调节速率为:
Figure BDA0003038229980000055
式中:
Figure BDA0003038229980000056
为第j个光伏电站的有功调节速率;Sj为第j个光伏电站的太阳辐照度;λ为光伏功率调节速率-太阳辐照度特性常数;εp-j为第j个光伏电站的与有功相关的光伏电池特性;
所述步骤2中的光伏电站的无功调节裕量为:
Figure BDA0003038229980000057
式中:
Figure BDA0003038229980000058
为第j个光伏电站的无功调节裕量;
Figure BDA0003038229980000059
为第j个光伏电站的逆变器总容量;Q0j为第j个光伏电站的在并网点电压等于电压参考值时的无功出力;-P0j为第j个光伏电站在系统频率等于频率参考值时的有功出力。
所述步骤2中的光伏电站的无功调节速率为:
Figure BDA00030382299800000510
式中:
Figure BDA00030382299800000511
为第j个光伏电站的无功调节速率;εq-j为第j个光伏电站的与无功相关的光伏电池特性;
Figure BDA00030382299800000512
为第j个光伏电站的有功调节速率。
与现有技术相比,本发明的有益效果在于:
通过设置与光伏电站运行特性、太阳辐照度、系统频率和并网点电压相关的变下垂曲线的下垂系数,一方面在进行有功-频率和无功-电压控制时根据系统频率和并网点电压的偏差大小,调节下垂曲线的下垂系数,增大有功-频率特性曲线和无功-电压特性曲线在运行点附近的下垂曲线的下垂系数,在系统频率和并网点电压偏差较小时充分使用光伏发电的有功和无功裕量,进一步减小频率和电压偏差,提高了系统的动态稳定性;另一方面,综合考虑了各光伏电站的太阳辐照度和光伏电站运行特性的影响,根据一次调频和调压过程中各光伏电站在当前太阳辐照度下的有功调节裕量、有功调节速率、无功调节裕量和无功调节速率对其出力进行合理分配,使裕量较大的光伏电站承担更多的调频和调压任务,各光伏电站充分发挥调频和调压能力,提高控制精度。
附图说明
图1为本发明实施例的基于光伏电站运行特性和光照强度的变下垂系数控制方法的控制框图;
图2为本发明实施例的有功-频率和无功-电压变下垂曲线的下垂系数控制流程图;
图3为本发明实施例中理想的有功-频率和无功-电压下垂特性曲线;
图4为发明本实施例中光伏电站有功-频率特性对比图;
图5为本发明实施例中光伏电站无功-电压特性对比图;
图6为本发明实施例中系统频率响应对比图;
图7为本发明实施例中光伏电站并网点电压响应对比图。
具体实施方式
以下,参照附图对本发明的实施方式进行说明。
本发明实施例提供了一种基于光伏电站运行特性和光照强度的变下垂系数控制方法,具体步骤包括:
步骤1:确定有功-频率控制基本模型和无功-电压控制基本模型,为充分发挥光伏发电参与系统调频和调压的潜力,在系统频率偏差或光伏并网点偏差较小时,尽可能多使用光伏发电的有功和无功裕量,进一步减小频率和电压偏差,即要求在正常运行点附近的下垂系数较大;与此同时,为充分发挥调节裕量较大的光伏发电快速调节潜力,在参与系统调频和调压的过程中,调节裕量较大的光伏场站应承担更多的调节任务,即要求调节裕量较大的场站下垂系数较大;如图1所示,本发明的控制方法框图中包含有功-频率控制基本模型A和无功-电压控制基本模型B,具体包括以下步骤:
步骤1-1:建立P-f(有功-频率)下垂曲线的基本函数模型;
Figure BDA0003038229980000071
式中:P为光伏电站的有功出力;f为系统频率;fref为频率参考值;P0为光伏电站在系统频率f等于频率参考值fref时的有功出力;e为自然对数;α1为与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数调节因子;α2为与光伏电站的有功调节速率有关的有功-频率下垂曲线的下垂系数调节因子;
步骤1-2:建立Q-U(无功-电压)下垂曲线的基本函数模型;
Figure BDA0003038229980000072
式中:Q为光伏电站的无功出力;U为并网点电压;Uref为电压参考值;Q0为光伏电站在并网点电压U等于电压参考值Uref时的无功出力;β1为与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数调节因子;β2为与光伏电站的无功调节速率有关的无功-电压下垂曲线的下垂系数调节因子;
步骤1-3:将步骤1-1的式(1)对系统频率f求导数得到光伏电站的有功-频率下垂曲线的下垂系数Kp
Figure BDA0003038229980000073
式中:Kp12,f)为光伏电站的有功-频率下垂曲线的下垂系数;
步骤1-4:将步骤1-2的式(2)对并网点电压U求导数得到光伏电站的无功-电压下垂曲线的下垂系数Kq
Figure BDA0003038229980000074
式中:Kq12,U)为光伏电站的无功-电压下垂曲线的下垂系数;
步骤2:采集光伏电站的运行特性和太阳辐照度,根据采集的光伏电站的运行特性和太阳辐照度,计算光伏电站的有功调节裕量ΔPmargin、光伏电站的有功调节速率ΔPspeed、光伏电站的无功调节裕量ΔQmargin和光伏电站的无功调节速率ΔQspeed,光伏电站的运行特性包含:光伏功率调节速率-太阳辐照度特性常数、光伏电站的与有功相关的光伏电池特性、光伏电站的与无功相关的光伏电池特性、光伏电站的逆变器总容量;
步骤2-1:将采集的太阳辐照度根据光伏电池的输出特性,得到光伏电站的最大有功出力:
Figure BDA0003038229980000081
步骤2-2:根据步骤2-1中的式(5),得到各光伏电站的最大有功出力:
Figure BDA0003038229980000082
式中:
Figure BDA0003038229980000083
为第j个光伏电站的最大有功出力;n为光伏电站个数;
步骤2-3:根据采集的光伏电站的运行特性和步骤2-2中的式(6)得到各光伏电站的有功调节裕量:
Figure BDA0003038229980000084
式中:
Figure BDA0003038229980000085
为第j个光伏电站的有功调节裕量;-P0j为第j个光伏电站在系统频率等于频率参考值时的有功出力;
进一步,计算步骤2中的光伏电站的有功调节速率:
Figure BDA0003038229980000086
式中:
Figure BDA0003038229980000087
为第j个光伏电站的有功调节速率;Sj为第j个光伏电站的太阳辐照度;λ为光伏功率调节速率-太阳辐照度特性常数;εp-j为第j个光伏电站的光伏电池特性;
进一步,计算步骤2中的光伏电站的无功调节裕量:
Figure BDA0003038229980000088
式中:
Figure BDA0003038229980000089
为第j个光伏电站的无功调节裕量;
Figure BDA00030382299800000810
为第j个光伏电站的逆变器总容量;Q0j为第j个光伏电站的在并网点电压等于电压参考值时的无功出力;-P0j为第j个光伏电站在系统频率等于频率参考值时的有功出力。
更进一步,计算步骤2中的光伏电站的无功调节速率为:
Figure BDA00030382299800000811
式中:
Figure BDA00030382299800000812
为第j个光伏电站的无功调节速率;εq-j为第j个光伏电站的与无功相关的光伏电池特性;
Figure BDA00030382299800000813
为第j个光伏电站的有功调节速率。
步骤3:根据步骤2中获得的光伏电站的有功调节裕量计算与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数的调节因子;根据步骤1中式(1)可知,当系统频率取最小值时,频率偏差最大,此时光伏电站的有功出力为
Figure BDA0003038229980000091
Figure BDA0003038229980000092
带入式(1)并结合式(7)可得与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数调节因子为:
Figure BDA0003038229980000093
根据传统下垂控制下的最大下垂系数除以标准测试环境下的光伏电站的有功调节速率获得有功调节速度系数ηp
进而可得到当前天气状况下的有功-频率下垂曲线的下垂系数Kp的最大值Kp_max
Figure BDA0003038229980000094
当系统频率等于频率参考值时,式(3)可以写作
Figure BDA0003038229980000095
此时结合式(11)-(12),可以得到与光伏电站的有功调节速率有关的有功-频率下垂曲线的下垂系数调节因子为:
Figure BDA0003038229980000096
步骤4:根据步骤2中获得的光伏电站的无功调节裕量计算与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数的调节因子;根据步骤1中式(2)可知,当并网点电压取最小值时,电压偏差最大,此时光伏电站的无功出力为
Figure BDA0003038229980000097
Figure BDA0003038229980000098
带入式(2)并结合式(9)可得与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数调节因子为:
Figure BDA0003038229980000099
根据传统下垂控制下的最大下垂系数除以标准测试环境下的光伏电站的无功调节速率获得无功调节速度系数ηq
进而可得到当前天气状况下的无功-电压下垂曲线的下垂系数Kq的最大值Kq_max
Figure BDA00030382299800000910
当并网点电压等于电压参考值时,式(4)写作
Figure BDA00030382299800000911
此时结合式(14)-(15),可得到与光伏电站的无功调节速率有关的无功-电压下垂曲线的下垂系数调节因子:
Figure BDA00030382299800000912
将式(11)、(13)、(14)、(16)确定的下垂曲线的下垂系数调节因子,分别带入式(3)、(4)得到有功-频率下垂曲线的下垂系数Kp和无功-电压下垂曲线的下垂系数Kq
进一步,确定有功-频率下垂控制依据公式(17)和无功-电压下垂控制依据公式(18)
ΔP=Kp(fref-f) (17)
ΔQ=Kq(Uref-U)。 (18)
步骤5:检测系统频率f,根据AGC控制指令以及步骤4确定的有功-频率下垂控制依据公式(17)获得光伏电站的有功参考值Pref,AGC控制指令设置为系统频率f等于频率参考值fref时的光伏电站的有功出力P0;检测并网点电压U,根据AVC控制指令以及步骤4确定的无功-电压下垂控制依据公式(18)获得光伏电站的无功参考值Qref,AVC控制指令设置为并网点电压U等于电压参考值Uref时的光伏电站的无功出力Q0,输出光伏电站的有功参考值和无功参考值给光伏并网逆变器,完成调频和调压任务。
Pref=Kp(fref-f)+P0 (19)
Qref=Kq(Uref-U)+Q0。 (20)
图3为本发明理想的有功-频率和无功-电压下垂曲线。有功-频率和无功-电压的下垂曲线的下垂系数在运行点附近较大,充分使用光伏电站的有功和无功裕量,减小系统频率和并网点电压偏差。且裕量较大的光伏电站的有功-频率和无功-电压的下垂系数较大,承担更多的调频和调压任务,充分发挥光伏电站的调频和调压能力。
提供一具体实施例,说明采用本发明的方法得到各光伏场站的功率曲线,实现协调控制效果的过程。
步骤1:本实施例的光伏并网系统中共有1个火电厂和2个光伏电站,按照图2所示的控制流程,首先建立以上有功-频率、无功-电压下垂曲线的基本函数模型。
步骤2:获取两个光伏电站的太阳辐照度和运行特性:C、D两个光伏电站处太阳辐照度均为1000W/m2,温度均为25℃。因此可得到2个光伏电站的最大输出有功功率为Pmax=[300 300](MW),当前有功出力值为P=[200 140](MW),即频率波动前的有功出力P0,无功出力值为Q=[60 40](MVar),即电压波动前的无功出力Q0,两光伏电站经变压器并网的并网点电压参考值均为330kV,频率参考值为50Hz。
步骤3:设两光伏电站逆变器参数相同,两光伏电站的有功调节速率与太阳辐照度成正比:ΔPspeed=[60 60](MW/s)。计算两光伏电站在当前天气特性下的有功调节裕量为ΔPmargin=[100 160](MW)。根据以上计算方式,可得α1=[200 320],α2=[4 2.5],则两光伏电站的有功-频率下垂曲线的基本函数为:
Figure BDA0003038229980000111
Figure BDA0003038229980000112
两光伏电站有功-频率特性对比图如图4所示。在传统下垂控制下其下垂系数均为常数值,而在本发明的变下垂系数控制中,有功-频率特性曲线在运行点附近的下垂系数较大,充分利用了光伏发电的有功裕量。而光伏电站C的有功调节裕量较小,光伏电站D的有功调节裕量较大,因此光伏电站C的下垂系数小于光伏电站D,充分发挥了光伏电站D的调频能力。
步骤4:设两光伏电站逆变器参数相同,两光伏电站的无功调节速率与有功调节速率成正比:ΔQspeed=[80 80](MVar/s)。计算两光伏电站在当前太阳辐照度下的无功调节裕量为ΔQmargin=[163.6 225.3](MVar)。根据以上计算方式,可得β1=[327.2 450.6],β2=[0.4 0.25],则两光伏电站的无功-电压下垂曲线的基本函数为:
Figure BDA0003038229980000113
Figure BDA0003038229980000114
两光伏电站无功-电压特性对比图如图5所示,在传统下垂控制下其下垂系数均为常数值,而在本发明的变下垂系数控制中,无功-电压特性曲线在运行点附近的下垂系数较大,充分利用了光伏发电的无功裕量。而光伏电站C的无功调节裕量较小,光伏电站D的无功调节裕量较大,因此光伏电站C的下垂曲线的下垂系数小于光伏电站D,充分发挥了光伏电站D的调压能力。
步骤5:检测系统频率和并网点电压,结合AGC和AVC的控制指令及下垂系数,输出光伏并网逆变器的有功无功参考值。系统频率响应对比图和光伏电站并网点电压响应对比图分别如图6和图7所示,由图6系统频率响应图可知,当采用传统下垂控制时,减去负荷后系统频率发生波动,上升了50.04Hz;当采用本发明的变下垂系数控制时,在相同扰动下,系统频率仅增加到50.02Hz。与传统下垂控制相比,本发明的变下垂系数控制降低了系统频率偏差,提高了系统的稳定性。由图7光伏电站并网点电压响应图可知,当采用传统下垂控制时,增加负荷后并网点电压下降到323kV左右;当采用本发明的变下垂系数控制时,在相同扰动下,并网点电压保持在330kV左右。与传统下垂控制相比,本发明的变下垂系数控制降低了并网点电压偏差,提高了系统的稳定性。
与传统下垂控制相比,本发明提出的基于光伏电站运行特性和光照强度的变下垂系数控制方法进一步减小了频率和电压偏差,提高了系统的动态稳定性,同时充分发挥光伏电站调频和调压能力,使裕量较大的光伏场站承担更多的调频和调压任务,提高了控制精度。
以上所述的实施例仅是对本发明的优选实施方式进行描述,并非对本发明的范围进行限定,在不脱离本发明设计精神的前提下,本领域普通技术人员对本发明的技术方案做出的各种变形和改进,均应落入本发明权利要求书确定的保护范围内。

Claims (10)

1.一种基于光伏电站运行特性和光照强度的变下垂系数控制方法,其特征在于,其包括以下步骤:
步骤1:确定有功-频率控制基本模型和无功-电压控制基本模型,具体包括以下步骤:
步骤1-1:建立有功-频率下垂曲线的基本函数模型;
Figure FDA0003597235540000011
式中:P为光伏电站的有功出力;f为系统频率;fref为频率参考值;P0为光伏电站在系统频率f等于频率参考值fref时的有功出力;e为自然对数;α1为与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数调节因子;α2为与光伏电站的有功调节速率有关的有功-频率下垂曲线的下垂系数调节因子;
步骤1-2:建立无功-电压下垂曲线的基本函数模型;
Figure FDA0003597235540000012
式中:Q为光伏电站的无功出力;U为并网点电压;Uref为电压参考值;Q0为光伏电站在并网点电压U等于电压参考值Uref时的无功出力;β1为与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数调节因子;β2为与光伏电站的无功调节速率有关的无功-电压下垂曲线的下垂系数调节因子;
步骤1-3:将所述步骤1-1的有功-频率下垂曲线的基本函数模型对系统频率f求导数得到光伏电站的有功-频率下垂曲线的下垂系数Kp
Figure FDA0003597235540000013
步骤1-4:将所述步骤1-2的无功-电压下垂曲线的基本函数模型对并网点电压U求导数得到光伏电站的无功-电压下垂曲线的下垂系数Kq
Figure FDA0003597235540000014
步骤2:采集光伏电站的运行特性和太阳辐照度,根据所述光伏电站的运行特性和太阳辐照度,计算光伏电站的有功调节裕量ΔPmargin、光伏电站的有功调节速率ΔPspeed、光伏电站的无功调节裕量ΔQmargin和光伏电站的无功调节速率ΔQspeed
步骤3:根据所述步骤2中获得的光伏电站的有功调节裕量ΔPmargin获得所述步骤1-1中与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数调节因子α1;根据所述步骤2中获得的光伏电站的有功调节速率ΔPspeed获得所述步骤1-1中与光伏电站的有功调节速率有关的有功-频率下垂曲线的下垂系数调节因子α2;根据所述与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数调节因子α1、与光伏电站的有功调节速率有关的有功-频率下垂曲线的下垂系数调节因子α2、系统频率f以及频率参考值fref,得到步骤1-3所述的光伏电站的有功-频率下垂曲线的下垂系数Kp
步骤4:根据所述步骤2中获得的光伏电站的无功调节裕量ΔQmargin获得所述步骤1-2中与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数调节因子β1;根据所述步骤2中获得的光伏电站的无功调节速率ΔQspeed获得所述步骤1-2中与光伏电站的无功调节速率有关的无功-电压下垂曲线的下垂系数调节因子β2;根据所述与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数调节因子β1与光伏电站的无功调节速率有关的无功-电压下垂曲线的下垂系数调节因子β2、并网点电压U以及电压参考值Uref,得到步骤1-4所述的光伏电站的无功-电压下垂曲线的下垂系数Kq
步骤5:检测系统频率f,根据AGC控制指令、频率参考值fref以及所述步骤3得到的有功-频率下垂曲线的下垂系数Kp获得光伏电站的有功参考值Pref;检测并网点电压U,根据AVC控制指令、电压参考值Uref以及所述步骤4得到的无功-电压下垂曲线的下垂系数Kq获得光伏电站的无功参考值Qref,输出获得的光伏电站的有功参考值Pref和无功参考值Qref给光伏并网逆变器,完成调频和调压任务。
2.根据权利要求1所述的基于光伏电站运行特性和光照强度的变下垂系数控制方法,其特征在于,所述步骤3中获得的与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数调节因子为:
Figure FDA0003597235540000021
式中:α1-j为第j个光伏电站与光伏电站的有功调节裕量有关的有功-频率下垂曲线的下垂系数调节因子;
Figure FDA0003597235540000022
为第j个光伏电站的有功调节裕量。
3.根据权利要求2所述的基于光伏电站运行特性和光照强度的变下垂系数控制方法,其特征在于,所述步骤3中获得的与光伏电站的有功调节速率有关的有功-频率下垂曲线的下垂系数调节因子为:
Figure FDA0003597235540000031
式中:α2-j为第j个光伏电站与光伏电站的有功调节速率有关的有功-频率下垂曲线的下垂系数调节因子;ηp为有功调节速率系数;
Figure FDA0003597235540000032
为第j个光伏电站的有功调节速率;
所述有功调节速率系数ηp为传统下垂控制下的最大下垂系数除以标准测试环境下的光伏电站的有功调节速率。
4.根据权利要求1所述的基于光伏电站运行特性和光照强度的变下垂系数控制方法,其特征在于,所述步骤4中获得的与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数调节因子为:
Figure FDA0003597235540000033
式中:β1-j为第j个光伏电站与光伏电站的无功调节裕量有关的无功-电压下垂曲线的下垂系数调节因子;
Figure FDA0003597235540000034
为第j个光伏电站的无功调节裕量。
5.根据权利要求4所述的基于光伏电站运行特性和光照强度的变下垂系数控制方法,其特征在于,所述步骤4中获得的与光伏电站的无功调节速率有关的无功-电压下垂曲线的下垂系数调节因子为:
Figure FDA0003597235540000035
式中:β2-j为第j个光伏电站与光伏电站的无功调节速率有关的无功-电压下垂曲线的下垂系数调节因子;ηq为无功调节速率系数;
Figure FDA0003597235540000036
为第j个光伏电站的无功调节速率;
所述无功调节速率系数ηq为传统下垂控制下的最大下垂系数除以标准测试环境下的光伏电站的无功调节速率。
6.根据权利要求1所述的基于光伏电站运行特性和光照强度的变下垂系数控制方法,其特征在于,所述步骤5中的AGC控制指令设置为系统频率f等于频率参考值fref时的光伏电站的有功出力P0,所述步骤5中获得的光伏电站的有功参考值Pref为:
Pref=Kp(fref-f)+P0
所述步骤5中的AVC控制指令设置为并网点电压U等于电压参考值Uref时的光伏电站的无功出力Q0,所述步骤5中获得的光伏电站的无功参考值Qref为:
Qref=Kq(Uref-U)+Q0
7.根据权利要求1所述的基于光伏电站运行特性和光照强度的变下垂系数控制方法,其特征在于,获得所述步骤2中光伏电站的有功调节裕量ΔPmargin的具体步骤为:
步骤2-1:将所述步骤2中获取的太阳辐照度根据光伏电池的输出特性,得到各光伏电站的最大有功出力:
Figure FDA0003597235540000041
式中:
Figure FDA0003597235540000042
为第j个光伏电站的最大有功出力;n为光伏电站个数;
步骤2-2:根据所述步骤2中获取的光伏电站的运行特性和步骤2-1中得到的各光伏电站的最大有功出力,得到各光伏电站的有功调节裕量:
Figure FDA0003597235540000043
式中:
Figure FDA0003597235540000044
为第j个光伏电站的有功调节裕量;P0j为第j个光伏电站在系统频率等于频率参考值时的有功出力。
8.根据权利要求1所述的基于光伏电站运行特性和光照强度的变下垂系数控制方法,其特征在于,所述步骤2中的光伏电站的有功调节速率为:
Figure FDA0003597235540000045
式中:
Figure FDA0003597235540000046
为第j个光伏电站的有功调节速率;Sj为第j个光伏电站的太阳辐照度;λ为光伏功率调节速率-太阳辐照度特性常数;εp-j为第j个光伏电站的与有功相关的光伏电池特性。
9.根据权利要求1所述的基于光伏电站运行特性和光照强度的变下垂系数控制方法,其特征在于,所述步骤2中的光伏电站的无功调节裕量为:
Figure FDA0003597235540000047
式中:
Figure FDA0003597235540000048
为第j个光伏电站的无功调节裕量;
Figure FDA0003597235540000049
为第j个光伏电站的逆变器总容量;Q0j为第j个光伏电站的在并网点电压等于电压参考值时的无功出力;P0j为第j个光伏电站在系统频率等于频率参考值时的有功出力。
10.根据权利要求1所述的基于光伏电站运行特性和光照强度的变下垂系数控制方法,其特征在于,所述步骤2中的光伏电站的无功调节速率为:
Figure FDA0003597235540000051
式中:
Figure FDA0003597235540000052
为第j个光伏电站的无功调节速率;εq-j为第j个光伏电站的与无功相关的光伏电池特性;
Figure FDA0003597235540000053
为第j个光伏电站的有功调节速率。
CN202110449595.9A 2021-04-25 2021-04-25 基于光伏电站运行特性和光照强度的变下垂系数控制方法 Active CN113162072B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110449595.9A CN113162072B (zh) 2021-04-25 2021-04-25 基于光伏电站运行特性和光照强度的变下垂系数控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110449595.9A CN113162072B (zh) 2021-04-25 2021-04-25 基于光伏电站运行特性和光照强度的变下垂系数控制方法

Publications (2)

Publication Number Publication Date
CN113162072A CN113162072A (zh) 2021-07-23
CN113162072B true CN113162072B (zh) 2022-08-09

Family

ID=76870396

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110449595.9A Active CN113162072B (zh) 2021-04-25 2021-04-25 基于光伏电站运行特性和光照强度的变下垂系数控制方法

Country Status (1)

Country Link
CN (1) CN113162072B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104578144A (zh) * 2014-12-15 2015-04-29 浙江大学 微电网中光伏逆变器反下垂控制方法
CN104810822A (zh) * 2015-05-22 2015-07-29 中国电力科学研究院 一种微网双向dc/dc变化下垂系数的控制方法
CN105162167A (zh) * 2015-09-30 2015-12-16 国网山东省电力公司电力科学研究院 一种基于自适应下垂控制的风光储微网调频方法
CN110783958A (zh) * 2019-11-13 2020-02-11 安徽工程大学 一种基于动态下垂系数的vsg控制方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102161812B1 (ko) * 2019-06-18 2020-10-05 연세대학교 산학협력단 유효 전력 제어를 이용한 태양광 발전 시스템의 mppt 제어를 위한 장치 및 방법

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104578144A (zh) * 2014-12-15 2015-04-29 浙江大学 微电网中光伏逆变器反下垂控制方法
CN104810822A (zh) * 2015-05-22 2015-07-29 中国电力科学研究院 一种微网双向dc/dc变化下垂系数的控制方法
CN105162167A (zh) * 2015-09-30 2015-12-16 国网山东省电力公司电力科学研究院 一种基于自适应下垂控制的风光储微网调频方法
CN110783958A (zh) * 2019-11-13 2020-02-11 安徽工程大学 一种基于动态下垂系数的vsg控制方法

Also Published As

Publication number Publication date
CN113162072A (zh) 2021-07-23

Similar Documents

Publication Publication Date Title
CN107732956A (zh) 变功率跟踪轨迹的两级式光伏并网系统低电压穿越方法
Zhong et al. DC-side synchronous active power control of two-stage photovoltaic generation for frequency support in Islanded microgrids
CN108683212B (zh) 一种基于功率解耦的混合储能型虚拟同步发电机控制方法
CN107658899A (zh) 一种风光柴微网频率协调控制方法
CN101860044A (zh) 风电场无功电压的协调控制方法
CN112332421A (zh) 基于自适应下垂控制的光伏电站参与电网电压调节方法
CN108123494A (zh) 基于最优转速功率追踪的双馈风机参与电网调频控制方法
CN110071526A (zh) 一种pmsg风机参与电网频率调节的自适应下垂控制方法
Zheng et al. Modeling and analysis of the AC/DC hybrid micro-grid with bidirectional power flow controller
CN108491025A (zh) 一种离线系数拟合的光伏发电系统主动减载方法
CN111446725A (zh) 一种用于微电网的混合储能调频控制方法
CN104850166A (zh) 一种采用滑模层极值搜索控制的光伏发电最大功率追踪方法
CN113162072B (zh) 基于光伏电站运行特性和光照强度的变下垂系数控制方法
Zhu et al. $\mu $-Synthesis Robust Control of Variable Speed Wind Turbine Generators for Participating in Microgrid Frequency Regulation
CN110112723B (zh) 一种直流微电网离网状态下的下垂控制方法
CN112865131A (zh) 一种参与电网频率调节的光伏逆变器运行控制方法和系统
CN113193570B (zh) 考虑光强和运行特性的光伏一次调频功率控制方法及装置
Lou et al. Analysis of primary frequency response based on overspeed and pitch control reserve and coordinated control strategy
CN115882524A (zh) 一种提升频率响应能力的风电机组控制参数整定方法
CN115940257A (zh) 一种光储发电虚拟惯量控制方法
He et al. Active LVRT strategy based on improved pv engineering model and power feedforward
CN110994669A (zh) 一种针对光伏电站集中式逆变器的控制方法及系统
CN112311020A (zh) 一种风电场无功优化调度方法和系统
CN110943459A (zh) 基于电压响应的多无功补偿装置并列运行控制方法及系统
Abo-Alela Design and Performance of Hybrid Wind-Solar Energy Generation System for Efficiency Improvement

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