CN107704696A - 并网光伏电站动态等值阻抗及仿真分析策略 - Google Patents
并网光伏电站动态等值阻抗及仿真分析策略 Download PDFInfo
- Publication number
- CN107704696A CN107704696A CN201710956027.1A CN201710956027A CN107704696A CN 107704696 A CN107704696 A CN 107704696A CN 201710956027 A CN201710956027 A CN 201710956027A CN 107704696 A CN107704696 A CN 107704696A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mfrac
- msup
- photovoltaic
- 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.)
- Pending
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 41
- 238000004458 analytical method Methods 0.000 title claims abstract description 14
- 101100499229 Mus musculus Dhrsx gene Proteins 0.000 claims abstract description 17
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 238000012545 processing Methods 0.000 claims abstract description 9
- 238000005286 illumination Methods 0.000 claims description 23
- 238000010586 diagram Methods 0.000 claims description 16
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000013178 mathematical model Methods 0.000 claims description 11
- 238000002474 experimental method Methods 0.000 claims description 5
- 239000003990 capacitor Substances 0.000 claims description 3
- 238000000034 method Methods 0.000 abstract description 6
- 238000010248 power generation Methods 0.000 description 6
- 238000012937 correction Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
-
- 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)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Microelectronics & Electronic Packaging (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Photovoltaic Devices (AREA)
Abstract
本发明涉及继电保护技术领域,是一种建立并网光伏电站动态等值阻抗及仿真分析策略,其特点是,包括的步骤有:根据并网光伏电站的拓扑结构定义动态等值阻抗模型。通过PSCAD/EMTDC建模仿真及MATLAB数据处理,将仿真结果和实际光伏电站实验测量结果进行比较,验证模型及动态等值阻抗的准确性,其误差在工程允许的6%以下,能够满足工程应用精度要求。该建模方法能清晰地反映出光伏发电系统的组成结构,且动态等值阻抗模型适应光伏发电的波动性和间歇性,满足各种环境和运行方式变化,在优化控制和继电保护方面具有很好的实用参考价值。
Description
技术领域
本发明涉及继电保护技术领域,是一种并网光伏电站动态等值阻抗及仿真分析策略。
背景技术
光伏并网发电是太阳能利用的发展趋势,在对光伏并网系统进行仿真研究时,精确的光伏电站数学模型有利于对整个光伏发电系统进行优化设计,提高系统总体效率,而数学模型中光伏电站动态等值阻抗则对光伏控制系统的优化以及继电保护可靠性的提高起到决定性作用。本建模方法能清晰地反映出光伏发电系统的组成结构,在优化控制和继电保护方面具有很好的实用参考价值。
发明内容
本发明所要解决的技术问题是,克服现有技术的不足,提供一种科学合理,判别准确,自愈性强,可靠性高的并网光伏电站动态等值阻抗及仿真分析策略。
解决其技术问题采用的技术方案是,一种并网光伏电站动态等值阻抗及仿真分析策略,其特征是,它包括的步骤有:
1)根据并网光伏电站的拓扑结构定义动态等值阻抗模型;
搭建光伏并网的拓扑结构,包括光伏电源,升压斩波电路,逆变器,滤波器,变压器,其中:Upv为光伏电源输出的电压,Ipv为光伏电源输出的电流;
a)定义光伏电源到斩波电路出口侧的等值阻抗
搭建光伏电源到斩波电路出口侧的等值电路图,其中,Udc为等值的电压源、Zdc动态等值阻抗,工程用太阳电池模型通常要求仅采用供应商提供的几个重要技术参数,如Isc、Uoc、Im、Um、Pm,就能在一定的精度下复现阵列的特性,并能便于计算机分析;
太阳能电池的I-U方程可简化为
在最大功率点时,U=Um,I=Im,可得
由于在常温条件下远大于1可忽略式中的-1项,解出C1
注意到开路状态下,当I=0时,U=Uoc,并把式(3)代入式(2)得解出C2
把式(3)和式(4)代入到式(1),即可得到光伏电池在标准条件下的输出特性表达式;
当外界条件发生变化时,为了得到不同光照强度及温度下光伏电池的输出特性,则需要对其在标准条件下的输出特性表达式进行修正,根据条件的变化对光伏电池输出电压和电流的变化量进行修正,考虑光强和温度影响的工程等值模型,采用式(5)计算,其中,S为任意光照强度,Sref为标准光照强度;T为环境温度,Tref为标准环境温度;Uoc、Isc、Um和Im为厂商提供的技术参数,其中,Uoc光伏电源的开路电压,Isc光伏电源的短路电流,Um最大功率点电压,Im最大功率点电流;α为在参考日照下的电流变化温度系数,0.015Amps/℃;β为在参考日照下的电压变化温度系数,取0.7V/℃;Rs为光伏列阵的串联电阻,与光伏组件内部单体光伏电池的串联并联方式有关;Upv为光伏电源出口侧的电压,Ipv为光伏电源出口侧的电流,C1、C2为简化光伏电池模型后求解的常系数,
当Ipv=0时,Upv=Uoc,此时的U’oc为任意光照和温度条件下光伏组件的开路电压
当Upv=0时,Ipv=ISC,此时的I’sc为任意光照和温度条件下光伏组件的短路电流
其次搭建Boost斩波电路,在较长的一段时间内,斩波电路的电容充电,其电压值逐渐上升到恒定值,即光伏阵列的开路电压Uocboost=UPVoc;
式中,Iscclose为开关闭合时的负载电流;RVQ为三极管通态时电阻,取值0.2Ω;RVD为斩波电路二极管漏电阻,取值0.2Ω;Isc为光伏组件的短路电流
Iscopen=I'sc (9)
式中Iscopen为开关开路时的负载电流。因此光伏电站短路电流的平均值为
式中D—Boost斩波电路的占空比,取值0.7,
式(11)为光伏电站的等值阻抗的数学模型,式(12)为光伏电站的等值电压源的数学模型,其中Ns为光伏组件的串联数,Np为光伏组件的并联数,
b)定义光伏电源到逆变器出口侧等值阻抗
搭建光伏电源到逆变器出口侧等值电路图,其中Ueq、Zeq分别为等值的电压源和动态等值阻抗,把逆变器等值成一个受控电压源和电阻串联,由于在逆变电路中的IGBT内阻很小,所以Zac可以忽略,利用戴维南定理可求得:
Zeq=βZdc Ueq=βUdc (13)
其中,M为调制度;L2为滤波器电抗值;Ls为网侧阻抗;Us为并网电压;η为逆变器效率;Ppv.m为光伏电源最大点功率;Upv.m为光伏电源最大功率点电压,
C)定义光伏电源到变压器出口侧等值阻抗
搭建光伏电源到变压器出口侧等值电路图,其中,Ups为等值的电压源、Zps为动态等值阻抗,ZL1、ZL2、ZC、为滤波器等值阻抗,ZT为变压器等值阻抗,利用戴维南定理可求得:
2)通过PSCAD/EMTDC建模仿真及MATLAB数据处理,将仿真结果和实际光伏电站实验测量结果进行比较,验证模型及动态等值阻抗的准确性,其误差在工程允许的6%以下,能够满足工程应用精度要求,
搭建66kV配电网典型拓扑图,利用PSCAD建模仿真及MATLAB数据处理比较仿真结果和光伏电站实验测量结果,
搭建光伏电源并入66kV配电网PSCAD仿真模型,光伏电站模块内部包括光伏电源模块,升压斩波电路,逆变器模块,滤波器模块和变压器模块,MPPT模块,控制模块等。线路参数:AB长度为26.38km,线路BC长度为21.99km,线路CD的长度为47.7km;Z1=0.0178+j0.314Ω/km,Z0=0.295+j1.04Ω/km;光伏电源部分参数:设置当前光照强度为400W/m2,温度为25℃;参考光照强度1000W/m2,参考环境温度25℃;在参考日照下的电流变化温度系数α取0.015Amps/℃;在参考日照下的电压变化温度系数β取0.7V/℃;参考厂商提供的技术参数为最大功率点功率260W,最大功率点工作电流8.52A,最大功率点工作电压30.51V,短路电流9.09A,开路电压37.65V;光伏组件的串联数Ns为400,光伏组件的并联数Np为400;光伏逆变器效率97.5%;其他参数:滤波器电抗值L2=0.037Ω,阻抗Ls=0.47+j8.28Ω。
利用PSCAD采集的实时数据与各组件出厂参数计算光伏电站等值阻抗Zps和等值电压源Ups,将测量的光伏电站出口处电压电流波形分别与仿真计算的等效输出电压和等效输出电流比较,由于数学模型中存在并网处不断更新的反馈量,仿真计算波形逐渐增大并趋于稳定,最后接近实验测量波形,为了更直观的验证该动态等值阻抗模型的准确性,定义其相对误差为:
其中,δi为电流相对误差,δu为电压相对误差,I为实验测量电流,Ips为仿真电流,U为实验测量电压,Ups为仿真电压。
本发明的一种并网光伏电站动态等值阻抗及仿真分析策略,按照光伏发电系统的各个组成部分的连接顺序,利用各部分的通用等值模型和戴维南定理进行综合建模,根据光照强度和温度、光伏组件的出厂参数、逆变器参数、滤波器参数、以及并网电压等定义动态等值阻抗模型,能清晰地反映出光伏发电系统的组成结构,且动态等值阻抗模型适应光伏发电的波动性和间歇性,满足各种环境和运行方式变化,在优化控制和继电保护方面具有很好的实用参考价值。
附图说明
图1为一种并网光伏电站动态等值阻抗及仿真分析策略流程图;
图2为光伏并网的拓扑结构图;
图3为光伏电源到斩波电路出口侧等值电路图;
图4为Boost斩波电路图;
图5为光伏电源到逆变器出口侧等值电路图;
图6为光伏电源到变压器出口侧等值电路图;
图7为66kV配电网典型拓扑图;
图8为光伏电源并入66kV配电网PSCAD仿真模型图;
图9为仿真波形与实验测量波形对比图;
图10为电流电压相对误差折线图。
具体实施方式
下面将结合附图和实施例对本发明的技术方案进行清楚、完整地描述。
参照图1和图7,本发明的一种并网光伏电站动态等值阻抗及仿真分析策略,以A相电压电流波形为例,包括如下步骤:
S1:根据并网光伏电站的拓扑结构定义动态等值阻抗模型;
步骤根据并网光伏电站的拓扑结构定义动态等值阻抗模型,如图2所示,具体包括:
S11:定义光伏电源到斩波电路出口侧的等值阻抗
搭建光伏电源到斩波电路出口侧的等值电路图,如图3所示,其中Udc、Zdc分别为等值的电压源和动态等值阻抗。工程用太阳电池模型通常要求仅采用供应商提供的几个重要技术参数,如Isc、Uoc、Im、Um、Pm,就能在一定的精度下复现阵列的特性,并能便于计算机分析;
太阳能电池的I-U方程可简化为
在最大功率点时,U=Um,I=Im,可得
由于在常温条件下远大于1可忽略式中的-1项,解出C1
注意到开路状态下,当I=0时,U=Uoc,并把式(3)代入式(2)得解出C2
把式(3)和式(4)代入到式(1),即可得到光伏电池在标准条件下的输出特性表达式。
当外界条件发生变化时,为了得到不同光照强度及温度下光伏电池的输出特性,则需要对其在标准条件下的输出特性表达式进行修正。有两种较常见的修正方法,一种是根据条件发生变化时,根据一般情况下的开路电压U’oc、短路电流I’sc、最大功率点的电压U’m、最大功率点的电流I’m计算新的C’1和C’2,然后带入公式(1)求出修正后的输出特性表达式。虽然这种修正方法可以得到不同条件下的光伏电池的4个基本性能参数,但计算复杂且计算量大,需要不断重复计算以更新C1和C2的值。本发明采用的修正方法是根据条件的变化对光伏电池输出电压和电流的变化量进行修正,采用考虑光强和温度影响的工程等值模型,如式(5)所示。其中S为任意光照强度,Sref为标准光照强度;T为环境温度,Tref为标准环境温度;Uoc,Isc,Um,Im为厂商提供的技术参数,分别为光伏电源的开路电压,短路电流,最大功率点电压和最大功率点电流;α为在参考日照下的电流变化温度系数,0.015Amps/℃;β为在参考日照下的电压变化温度系数,取0.7V/℃;Rs为光伏列阵的串联电阻,与光伏组件内部单体光伏电池的串联并联方式有关;Upv,Ipv分别为光伏电源出口侧的电压电流。C1,C2为简化光伏电池模型后求解的常系数。
当Ipv=0时,Upv=Uoc,此时的U’oc为任意光照和温度条件下光伏组件的开路电压
当Upv=0时,Ipv=ISC,此时的I’sc为任意光照和温度条件下光伏组件的短路电流
其次搭建Boost斩波电路,如图4所示,在较长的一段时间内,斩波电路的电容充电,其电压值逐渐上升到恒定值,即光伏阵列的开路电压Uocboost=UPVoc;
式中Iscclose为开关闭合时的负载电流;RVQ为三极管通态时电阻,取值0.2Ω;RVD为斩波电路二极管漏电阻,取值0.2Ω;Isc为光伏组件的短路电流。
Iscopen=I'sc (9)
式中Iscopen为开关开路时的负载电流。因此光伏电站短路电流的平均值为
式中D—Boost斩波电路的占空比,取值0.7。
式(11)为光伏电站的等值阻抗的数学模型,式(12)为光伏电站的等值电压源的数学模型,其中Ns为光伏组件的串联数,Np为光伏组件的并联数。
S12:定义光伏电源到逆变器出口侧等值阻抗
搭建光伏电源到逆变器出口侧等值电路图,如图5所示,其中Ueq、Zeq分别为等值的电压源和动态等值阻抗。把逆变器等值成一个受控电压源和电阻串联,由于在逆变电路中的IGBT内阻很小,所以Zac可以忽略,利用戴维南定理可求得:
Zeq=βZdc Ueq=βUdc (13)
其中M为调制度;L2为滤波器电抗值;Ls为网侧阻抗;Us为并网电压;η为逆变器效率;Ppv.m为光伏电源最大点功率;Upv.m为光伏电源最大功率点电压。
S13:定义光伏电源到变压器出口侧等值阻抗;
搭建光伏电源到变压器出口侧等值电路图,如图6所示,
其中Ups、Zps分别为等值的电压源和动态等值阻抗。ZL1、ZL2、ZC、为滤波器等值阻抗,ZT为变压器等值阻抗,利用戴维南定理可求得:
其中:
(1)S2:通过PSCAD/EMTDC建模仿真及MATLAB数据处理,将仿真结果和实际光伏电站实验测量结果进行比较,验证模型及动态等值阻抗的准确性,其误差在工程允许的6%以下,能够满足工程应用精度要求;
步骤S2包括通过PSCAD/EMTDC建模仿真及MATLAB数据处理,将仿真结果和实际光伏电站实验测量结果进行比较,验证模型及动态等值阻抗的准确性,其误差在工程允许的6%以下,能够满足工程应用精度要求,具体包括:
参照图7,搭建66kV配电网典型拓扑图,
利用PSCAD建模仿真及MATLAB数据处理比较仿真结果和光伏电站实验测量结果。
如图8所示,搭建光伏电源并入66kV配电网PSCAD仿真模型,光伏电站模块内部包括光伏电源模块,升压斩波电路,逆变器模块,滤波器模块和变压器模块,MPPT模块,控制模块等。线路参数:AB长度为26.38km,线路BC长度为21.99km,线路CD的长度为47.7km;Z1=0.0178+j0.314Ω/km,Z0=0.295+j1.04Ω/km;光伏电源部分参数:设置当前光照强度为400W/m2,温度为25℃;参考光照强度1000W/m2,参考环境温度25℃;在参考日照下的电流变化温度系数α取0.015Amps/℃;在参考日照下的电压变化温度系数β取0.7V/℃;参考厂商提供的技术参数为最大功率点功率260W,最大功率点工作电流8.52A,最大功率点工作电压30.51V,短路电流9.09A,开路电压37.65V;光伏组件的串联数Ns为400,光伏组件的并联数Np为400;光伏逆变器效率97.5%;其他参数:滤波器电抗值L2=0.037Ω,阻抗Ls=0.47+j8.28Ω。
利用PSCAD采集的实时数据与各组件出厂参数计算光伏电站等值阻抗Zps和等值电压源Ups,将测量的光伏电站出口处电压电流波形分别与仿真计算的等效输出电压和等效输出电流比较。由于数学模型中存在并网处不断更新的反馈量,仿真计算波形逐渐增大并趋于稳定,最后接近实验测量波形,为了更直观的验证该动态等值阻抗模型的准确性,定义其相对误差为:
其中δi为电流相对误差,δu为电压相对误差,I和Ips分别为实验测量电流和仿真电流,U和Ups分别为实验测量电压和仿真电压。
根据上述步骤,如图9所示,以A相电压电流波形为例,仿真计算波形逐渐增大并趋于稳定,最后接近实验测量波形。如图10电流电压误差折线图所示,经计算δi和δu误差在0.1秒前大幅度下降,0.3秒后误差逐渐趋于稳定,最后稳定在工程允许的6%以下,能够满足工程应用精度的要求。
本发明的并网光伏电站动态等值阻抗及仿真分析策略的仿真软件程序的编制,依据自动控制技术和计算机处理技术编制,是本领域技术人员所熟悉的技术。
显然,上述实施例仅仅是为清楚地说明所作的举例,而并非对实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引伸出的显而易见的变化或变动仍处于本发明创造的保护范围之中。
Claims (1)
1.一种并网光伏电站动态等值阻抗及仿真分析策略,其特征是,它包括的步骤有:
1)根据并网光伏电站的拓扑结构定义动态等值阻抗模型;
搭建光伏并网的拓扑结构,包括光伏电源,升压斩波电路,逆变器,滤波器,变压器,其中:Upv为光伏电源输出的电压,Ipv为光伏电源输出的电流;
a)定义光伏电源到斩波电路出口侧的等值阻抗
搭建光伏电源到斩波电路出口侧的等值电路图,其中,Udc为等值的电压源、Zdc动态等值阻抗,工程用太阳电池模型通常要求仅采用供应商提供的几个重要技术参数,如Isc、Uoc、Im、Um、Pm,就能在一定的精度下复现阵列的特性,并能便于计算机分析;
太阳能电池的I-U方程可简化为
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<msup>
<mi>e</mi>
<mfrac>
<msub>
<mi>U</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
</mrow>
</msub>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
在最大功率点时,U=Um,I=Im,可得
<mrow>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<mo>=</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<msup>
<mi>e</mi>
<mfrac>
<msub>
<mi>U</mi>
<mi>m</mi>
</msub>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
由于在常温条件下远大于1可忽略式中的-1项,解出C1
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
</mfrac>
<mo>)</mo>
</mrow>
<msup>
<mi>e</mi>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>U</mi>
<mi>m</mi>
</msub>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>)</mo>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
注意到开路状态下,当I=0时,U=Uoc,并把式(3)代入式(2)得解出C2
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mfrac>
<msub>
<mi>U</mi>
<mi>m</mi>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mfrac>
<msub>
<mi>I</mi>
<mi>m</mi>
</msub>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
</mfrac>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
把式(3)和式(4)代入到式(1),即可得到光伏电池在标准条件下的输出特性表达式;
当外界条件发生变化时,为了得到不同光照强度及温度下光伏电池的输出特性,则需要对其在标准条件下的输出特性表达式进行修正,根据条件的变化对光伏电池输出电压和电流的变化量进行修正,考虑光强和温度影响的工程等值模型,采用式(5)计算,其中,S为任意光照强度,Sref为标准光照强度;T为环境温度,Tref为标准环境温度;Uoc、Isc、Um和Im为厂商提供的技术参数,其中,Uoc光伏电源的开路电压,Isc光伏电源的短路电流,Um最大功率点电压,Im最大功率点电流;α为在参考日照下的电流变化温度系数,0.015Amps/℃;β为在参考日照下的电压变化温度系数,取0.7V/℃;Rs为光伏列阵的串联电阻,与光伏组件内部单体光伏电池的串联并联方式有关;Upv为光伏电源出口侧的电压,Ipv为光伏电源出口侧的电流,C1、C2为简化光伏电池模型后求解的常系数,
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>I</mi>
<mo>=</mo>
<mi>&alpha;</mi>
<mfrac>
<mi>S</mi>
<msub>
<mi>S</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
</mfrac>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mfrac>
<mi>S</mi>
<msub>
<mi>S</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
</mfrac>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>&Delta;</mi>
<mi>U</mi>
<mo>=</mo>
<mo>-</mo>
<mi>&beta;</mi>
<mrow>
<mo>(</mo>
<mi>T</mi>
<mo>-</mo>
<msub>
<mi>T</mi>
<mrow>
<mi>r</mi>
<mi>e</mi>
<mi>f</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>R</mi>
<mi>s</mi>
</msub>
<mi>&Delta;</mi>
<mi>I</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<msup>
<mi>e</mi>
<mfrac>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>U</mi>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>I</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
当Ipv=0时,Upv=Uoc,此时的U’oc为任意光照和温度条件下光伏组件的开路电压
<mrow>
<msub>
<msup>
<mi>U</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<mi>&Delta;</mi>
<mi>U</mi>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
<mi>l</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<mfrac>
<mn>1</mn>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>I</mi>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
当Upv=0时,Ipv=ISC,此时的I’sc为任意光照和温度条件下光伏组件的短路电流
<mrow>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mo>(</mo>
<mrow>
<msup>
<mi>e</mi>
<mfrac>
<mrow>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>U</mi>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>I</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
其次搭建Boost斩波电路,在较长的一段时间内,斩波电路的电容充电,其电压值逐渐上升到恒定值,即光伏阵列的开路电压Uocboost=UPVoc;
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
<mi>c</mi>
<mi>l</mi>
<mi>o</mi>
<mi>s</mi>
<mi>e</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>&times;</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>Q</mi>
</mrow>
</msub>
<mrow>
<mn>3</mn>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,Iscclose为开关闭合时的负载电流;RVQ为三极管通态时电阻,取值0.2Ω;RVD为斩波电路二极管漏电阻,取值0.2Ω;Isc为光伏组件的短路电流
Iscopen=I'sc (9)
式中Iscopen为开关开路时的负载电流。因此光伏电站短路电流的平均值为
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>&times;</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>D</mi>
</mrow>
</msub>
<mrow>
<mn>3</mn>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>D</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
式中D—Boost斩波电路的占空比,取值0.7,
式(11)为光伏电站的等值阻抗的数学模型,式(12)为光伏电站的等值电压源的数学模型,其中Ns为光伏组件的串联数,Np为光伏组件的并联数,
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>p</mi>
</msub>
<msub>
<mi>I</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<msub>
<msup>
<mi>U</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>p</mi>
</msub>
<msub>
<msup>
<mi>I</mi>
<mo>&prime;</mo>
</msup>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>D</mi>
<mo>&times;</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>D</mi>
</mrow>
</msub>
<mrow>
<mn>3</mn>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>D</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mi>&Delta;</mi>
<mi>U</mi>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
<mi>ln</mi>
<mrow>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>I</mi>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>e</mi>
<mfrac>
<mrow>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>U</mi>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>I</mi>
</mrow>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mrow>
<mi>D</mi>
<mo>&times;</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>D</mi>
</mrow>
</msub>
<mrow>
<mn>3</mn>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>D</mi>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>d</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>U</mi>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>I</mi>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>12</mn>
<mo>)</mo>
</mrow>
</mrow>
b)定义光伏电源到逆变器出口侧等值阻抗
搭建光伏电源到逆变器出口侧等值电路图,其中Ueq、Zeq分别为等值的电压源和动态等值阻抗,把逆变器等值成一个受控电压源和电阻串联,由于在逆变电路中的IGBT内阻很小,所以Zac可以忽略,利用戴维南定理可求得:
Zeq=βZdc Ueq=βUdc (13)
其中,M为调制度;L2为滤波器电抗值;Ls为网侧阻抗;Us为并网电压;η为逆变器效率;Ppv.m为光伏电源最大点功率;Upv.m为光伏电源最大功率点电压,
<mrow>
<mi>M</mi>
<mo>=</mo>
<mfrac>
<mrow>
<mn>8</mn>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&eta;N</mi>
<mi>p</mi>
</msub>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
<mo>.</mo>
<mi>m</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msqrt>
<mn>3</mn>
</msqrt>
<msub>
<mi>U</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mo>(</mo>
<mrow>
<msub>
<mi>&omega;L</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&omega;L</mi>
<mi>s</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>U</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
<mo>.</mo>
<mi>m</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>14</mn>
<mo>)</mo>
</mrow>
</mrow>
C)定义光伏电源到变压器出口侧等值阻抗
搭建光伏电源到变压器出口侧等值电路图,其中,Ups为等值的电压源、Zps为动态等值阻抗,ZL1、ZL2、ZC、为滤波器等值阻抗,ZT为变压器等值阻抗,利用戴维南定理可求得:
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>p</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mo>(</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mn>2</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mi>T</mi>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>15</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>p</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>U</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mrow>
<mi>L</mi>
<mn>1</mn>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>Z</mi>
<mi>c</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
其中:
<mrow>
<msub>
<mi>Z</mi>
<mrow>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<mn>2</mn>
<msqrt>
<mn>2</mn>
</msqrt>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>U</mi>
<mo>+</mo>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
<mi>l</mi>
<mi>n</mi>
<mo>(</mo>
<mrow>
<mfrac>
<mn>1</mn>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
</mfrac>
<mo>+</mo>
<mfrac>
<mrow>
<mi>&Delta;</mi>
<mi>I</mi>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<mfrac>
<mrow>
<msub>
<mi>&eta;N</mi>
<mi>p</mi>
</msub>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<msub>
<mi>P</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
<mo>.</mo>
<mi>m</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msqrt>
<mn>3</mn>
</msqrt>
<msub>
<mi>U</mi>
<mi>s</mi>
</msub>
</mrow>
</mfrac>
<mo>(</mo>
<mrow>
<msub>
<mi>&omega;L</mi>
<mn>2</mn>
</msub>
<mo>+</mo>
<msub>
<mi>&omega;L</mi>
<mi>s</mi>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msubsup>
<mi>U</mi>
<mi>s</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
<mrow>
<msub>
<mi>N</mi>
<mi>p</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>s</mi>
<mi>c</mi>
</mrow>
</msub>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<msub>
<mi>C</mi>
<mn>1</mn>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msup>
<mi>e</mi>
<mfrac>
<mrow>
<mo>-</mo>
<mi>&Delta;</mi>
<mi>U</mi>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mn>2</mn>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>o</mi>
<mi>c</mi>
</mrow>
</msub>
</mrow>
</mfrac>
</msup>
<mo>-</mo>
<mn>1</mn>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>+</mo>
<mi>&Delta;</mi>
<mi>I</mi>
<mo>)</mo>
</mrow>
<mo>&times;</mo>
<mrow>
<mo>(</mo>
<mi>D</mi>
<mo>&times;</mo>
<mfrac>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>D</mi>
</mrow>
</msub>
<mrow>
<mn>3</mn>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>Q</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>R</mi>
<mrow>
<mi>V</mi>
<mi>D</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>+</mo>
<mo>(</mo>
<mrow>
<mn>1</mn>
<mo>-</mo>
<mi>D</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>N</mi>
<mi>s</mi>
</msub>
<msub>
<mi>U</mi>
<mrow>
<mi>p</mi>
<mi>v</mi>
<mo>.</mo>
<mi>m</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>17</mn>
<mo>)</mo>
</mrow>
</mrow>
2)通过PSCAD/EMTDC建模仿真及MATLAB数据处理,将仿真结果和实际光伏电站实验测量结果进行比较,验证模型及动态等值阻抗的准确性,其误差在工程允许的6%以下,能够满足工程应用精度要求,
搭建66kV配电网典型拓扑图,利用PSCAD建模仿真及MATLAB数据处理比较仿真结果和光伏电站实验测量结果,
搭建光伏电源并入66kV配电网PSCAD仿真模型,光伏电站模块内部包括光伏电源模块,升压斩波电路,逆变器模块,滤波器模块和变压器模块,MPPT模块,控制模块等。线路参数:AB长度为26.38km,线路BC长度为21.99km,线路CD的长度为47.7km;Z1=0.0178+j0.314Ω/km,Z0=0.295+j1.04Ω/km;光伏电源部分参数:设置当前光照强度为400W/m2,温度为25℃;参考光照强度1000W/m2,参考环境温度25℃;在参考日照下的电流变化温度系数α取0.015Amps/℃;在参考日照下的电压变化温度系数β取0.7V/℃;参考厂商提供的技术参数为最大功率点功率260W,最大功率点工作电流8.52A,最大功率点工作电压30.51V,短路电流9.09A,开路电压37.65V;光伏组件的串联数Ns为400,光伏组件的并联数Np为400;光伏逆变器效率97.5%;其他参数:滤波器电抗值L2=0.037Ω,阻抗Ls=0.47+j8.28Ω。
利用PSCAD采集的实时数据与各组件出厂参数计算光伏电站等值阻抗Zps和等值电压源Ups,将测量的光伏电站出口处电压电流波形分别与仿真计算的等效输出电压和等效输出电流比较,由于数学模型中存在并网处不断更新的反馈量,仿真计算波形逐渐增大并趋于稳定,最后接近实验测量波形,为了更直观的验证该动态等值阻抗模型的准确性,定义其相对误差为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&delta;</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mrow>
<mi>I</mi>
<mo>-</mo>
<msub>
<mi>I</mi>
<mrow>
<mi>p</mi>
<mi>s</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>I</mi>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&delta;</mi>
<mi>u</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<mi>U</mi>
<mo>-</mo>
<msub>
<mi>U</mi>
<mrow>
<mi>p</mi>
<mi>s</mi>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mo>/</mo>
<mi>U</mi>
<mo>&times;</mo>
<mn>100</mn>
<mi>%</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>18</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,δi为电流相对误差,δu为电压相对误差,I为实验测量电流,Ips为仿真电流,U为实验测量电压,Ups为仿真电压。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710956027.1A CN107704696A (zh) | 2017-10-15 | 2017-10-15 | 并网光伏电站动态等值阻抗及仿真分析策略 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710956027.1A CN107704696A (zh) | 2017-10-15 | 2017-10-15 | 并网光伏电站动态等值阻抗及仿真分析策略 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107704696A true CN107704696A (zh) | 2018-02-16 |
Family
ID=61183782
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710956027.1A Pending CN107704696A (zh) | 2017-10-15 | 2017-10-15 | 并网光伏电站动态等值阻抗及仿真分析策略 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107704696A (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109842155A (zh) * | 2019-03-05 | 2019-06-04 | 协鑫能源工程有限公司 | 光伏组件串并联损失评估方法、装置以及计算机设备 |
CN109888820A (zh) * | 2019-02-14 | 2019-06-14 | 国网甘肃省电力公司嘉峪关供电公司 | 用于计算新能源机组的短路电流贡献值的方法 |
CN110135118A (zh) * | 2019-06-11 | 2019-08-16 | 华北电力大学(保定) | 一种汽轮发电机阻尼有限元模型的校验方法 |
CN110581567A (zh) * | 2019-08-08 | 2019-12-17 | 国网山东省电力公司济南市历城区供电公司 | 实时跟踪内阻匹配的功率传输方法和系统 |
CN111914421A (zh) * | 2020-08-04 | 2020-11-10 | 合肥工业大学 | 针对不同故障机理的光伏支路统一建模方法 |
CN113241736A (zh) * | 2021-04-12 | 2021-08-10 | 中国电力科学研究院有限公司 | 一种用于新能源送出线路保护适应性分析的方法及系统 |
CN114841067A (zh) * | 2022-05-06 | 2022-08-02 | 南京轻机包装机械有限公司 | 一种基于matlab的全功率风力发电变流器IGBT开路故障仿真方法 |
CN114967823A (zh) * | 2022-06-07 | 2022-08-30 | 河北工程大学 | 基于改进黑寡妇的光伏最大功率点跟踪的控制方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101882896A (zh) * | 2010-07-12 | 2010-11-10 | 山东电力研究院 | 大规模光伏电站动态等值阻抗的建模方法 |
US20130088086A1 (en) * | 2011-10-07 | 2013-04-11 | Sma Solar Technology Ag | Circuitry arrangement for a solar power plant comprising a dc voltage source for an offset voltage |
CN103310121A (zh) * | 2013-07-10 | 2013-09-18 | 湖南大学 | 一种大型光伏电站与配网谐波交互影响分析模型建模方法 |
CN103730910A (zh) * | 2013-12-20 | 2014-04-16 | 南京南瑞集团公司 | 一种大规模光伏电站并网的动态等值方法 |
CN105425111A (zh) * | 2015-12-16 | 2016-03-23 | 国网安徽省电力公司 | 一种利用过渡电阻特征的输电线路故障测距方法 |
-
2017
- 2017-10-15 CN CN201710956027.1A patent/CN107704696A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101882896A (zh) * | 2010-07-12 | 2010-11-10 | 山东电力研究院 | 大规模光伏电站动态等值阻抗的建模方法 |
US20130088086A1 (en) * | 2011-10-07 | 2013-04-11 | Sma Solar Technology Ag | Circuitry arrangement for a solar power plant comprising a dc voltage source for an offset voltage |
CN103310121A (zh) * | 2013-07-10 | 2013-09-18 | 湖南大学 | 一种大型光伏电站与配网谐波交互影响分析模型建模方法 |
CN103730910A (zh) * | 2013-12-20 | 2014-04-16 | 南京南瑞集团公司 | 一种大规模光伏电站并网的动态等值方法 |
CN105425111A (zh) * | 2015-12-16 | 2016-03-23 | 国网安徽省电力公司 | 一种利用过渡电阻特征的输电线路故障测距方法 |
Non-Patent Citations (1)
Title |
---|
侯瑞鹏 等: "继电保护整定计算中光伏发电站阻抗的计算方法研究", 《内蒙古电力技术》 * |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109888820A (zh) * | 2019-02-14 | 2019-06-14 | 国网甘肃省电力公司嘉峪关供电公司 | 用于计算新能源机组的短路电流贡献值的方法 |
CN109888820B (zh) * | 2019-02-14 | 2021-11-02 | 国网甘肃省电力公司嘉峪关供电公司 | 用于计算新能源机组的短路电流贡献值的方法 |
CN109842155A (zh) * | 2019-03-05 | 2019-06-04 | 协鑫能源工程有限公司 | 光伏组件串并联损失评估方法、装置以及计算机设备 |
CN110135118A (zh) * | 2019-06-11 | 2019-08-16 | 华北电力大学(保定) | 一种汽轮发电机阻尼有限元模型的校验方法 |
CN110135118B (zh) * | 2019-06-11 | 2023-04-14 | 华北电力大学(保定) | 一种汽轮发电机阻尼有限元模型的校验方法 |
CN110581567A (zh) * | 2019-08-08 | 2019-12-17 | 国网山东省电力公司济南市历城区供电公司 | 实时跟踪内阻匹配的功率传输方法和系统 |
CN111914421A (zh) * | 2020-08-04 | 2020-11-10 | 合肥工业大学 | 针对不同故障机理的光伏支路统一建模方法 |
CN111914421B (zh) * | 2020-08-04 | 2024-02-27 | 合肥工业大学 | 针对不同故障机理的光伏支路统一建模方法 |
CN113241736A (zh) * | 2021-04-12 | 2021-08-10 | 中国电力科学研究院有限公司 | 一种用于新能源送出线路保护适应性分析的方法及系统 |
CN114841067A (zh) * | 2022-05-06 | 2022-08-02 | 南京轻机包装机械有限公司 | 一种基于matlab的全功率风力发电变流器IGBT开路故障仿真方法 |
CN114967823A (zh) * | 2022-06-07 | 2022-08-30 | 河北工程大学 | 基于改进黑寡妇的光伏最大功率点跟踪的控制方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107704696A (zh) | 并网光伏电站动态等值阻抗及仿真分析策略 | |
Aashoor et al. | A variable step size perturb and observe algorithm for photovoltaic maximum power point tracking | |
CN104778352B (zh) | 基于stft算子的七参数光伏电池输出特性建模方法 | |
CN103715719B (zh) | 一种适用于电力系统全过程动态仿真的光伏模型建模方法 | |
Reisi et al. | Classification and comparison of maximum power point tracking techniques for photovoltaic system: A review | |
Iqbal et al. | Design and simulation of a PV System with battery storage using bidirectional DC-DC converter using Matlab Simulink | |
Yusivar et al. | Buck-converter photovoltaic simulator | |
Zheng et al. | Modeling and analysis of the AC/DC hybrid micro-grid with bidirectional power flow controller | |
Hekss et al. | Nonlinear control strategy of single phase half bridge shunt active power filter interfacing renewable energy source and grid | |
CN107332241A (zh) | 用于评估含分布式新能源配电网可靠性的方法及系统 | |
Plakhtii et al. | Research of operating modes and features of integration of renewable energy sources into the electric power system | |
CN106970677B (zh) | 一种太阳能逆变器快速mppt方法 | |
Kane et al. | Characterization of PV cell-environmental factors consideration | |
Villarejo et al. | Comparison of central vs distributed inverters: application to photovoltaic systems | |
Vika | Modelling of photovoltaic modules with battery energy storage in Simulink/Matlab: With in-situ measurement comparisons | |
Navada et al. | Modelling of a solar photovoltaic power plant for power system studies | |
Ding et al. | Improved global maximum power point tracking method based on voltage interval for PV array under partially shaded conditions | |
CN110456223B (zh) | 一种含分布式电源和电动汽车的配电网短路电流测算方法 | |
Assia et al. | Detailed modeling of two diode photovoltaic module using MATLAB simulik | |
Yang et al. | An improved variable step size MPPT method for photovoltaic array under partial shading conditions | |
Lu et al. | Research and simulation on photovoltaic power system maximum power control | |
Abutaima et al. | Design of Improved Incremental Conductance with Fast Intelligent (FI) Based MPPT Technique for Solar PV System | |
Li et al. | A piecewise function based equivalent aggregation of photovoltaic array considering partial shading situation | |
Kong et al. | Dynamic equivalent impedance of photovoltaic plant and simulation analysis | |
Zhang et al. | Parameter identification and optimization of PI controllers for islanded DC micro-grids |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20180216 |
|
WD01 | Invention patent application deemed withdrawn after publication |