CN102594215B - 一种光伏电站模型参数辨识方法 - Google Patents

Abstract

本发明公开了一种光伏电站模型参数辨识方法，该方法为：将光伏电站模型划分为外环控制和内环控制两个辨识环节，外环控制参数辨识时其输入数据来自于光伏电站的实际运行数据；内环控制参数辨识时其输入数据来自于对光伏电站逆变器的现场测试试验数据。内环控制参数辨识采用最小二乘法，外环控制参数辨识采用组合优化辨识方法。通过对光伏电站的模型参数辨识，可以得到光伏组件实际最大功率运行点数据、逆变器关键控制参数，从而可使建立的光伏电站仿真模型具有较高的仿真精度，更接近光伏电站实际运行状态。

Description

一种光伏电站模型参数辨识方法

技术领域

本发明属新能源应用领域，涉及一种光伏电站模型参数的辨识方法。

背景技术

光伏发电具有波动性和间歇性，大规模光伏电站并网运行会影响电力系统的安全稳定经济运行。对光伏电站模型参数进行辨识有助于电网规划部门准确分析大型光伏电站接入对电网的影响、有利于统筹安排常规电源和光伏电站的协调配合，适时及时地调整调度计划，合理安排电网运行方式。一方面有效地减轻光伏接入对电网的不利影响，提高电力系统运行的安全性和稳定性，另一方面通过合理布置光伏电站的接入，以充分利用太阳能资源，获得更大的经济效益和社会效益。

为了得到具有较高可靠信度的光伏电站等效模型，需要对模型中的参数包括光伏组件最大功率运行点参数、逆变器最大功率跟踪控制参数、逆变器内环跟踪控制参数等进行辨识。而目前还未见运用光伏电实际运行数据和对逆变器现场测试数据进行模型参数辨识的报道。

发明内容

本光伏电站模型参数辨识方法，将光伏电站模型划分为外环控制和内环控制两个辨识环节，外环控制参数辨识时其输入数据来自于光伏电站的实际运行数据，主要辨识光伏组件最大功率运行点参数、最大功率点跟踪(MPPT)控制参数；内环控制参数辨识时其输入数据来自于对光伏电站逆变器的现场测试试验数据，用于辨识逆变器动态跟踪控制参数。本发明特征在于，该方法包括如下步骤：

A、将光伏电站模型划分为外环控制和内环控制两个辨识环节，外环控制参数辨识，辨识光伏组件最大功率运行点参数、最大功率点跟踪(MPPT)控制参数；内环控制参数辨识，用于辨识逆变器动态跟踪控制参数；

B、对光伏电站输入源数据进行预处理：外环控制参数辨识时其输入数据来自于光伏电站的实际运行数据，应对运行数据进行筛选,去除其中的错误数据，然后将合格的数据添加到辨识数据库；内环控制参数辨识时其输入数据来自于对逆变器的现场施加阶跃信号得到响应曲线数据；

C、内环控制参数辨识采用最小二乘法，其辨识流程：

(1)建立内环开环传递函数：

$W_{\mathrm{ci}}\left(s\right)=\frac{K_{\mathrm{iP}}{K}_{\mathrm{PWM}}}{T_{i}L}\frac{T_{i}s+1}{s^2(1.5T_{s}s+1)}$

其中，K_iP为比例环节系数、K_PWM为PMW控制环节的增益，T_i、T_s积分环节系数；

(2)对现场施加阶跃信号得到的响应曲线数据进行相关转换处理，包括数据的有效性分析、电气量程与工程量程的转换、标幺化处理共三个步骤；

(3)将转换后的响应数据和阶跃输入信号为输入，采用最小二乘法计算方法，对内环传递函数中的参数K_iP、K_PWM及T_i进行辨识；

其中，K_iP为比例环节系数，K_PWM为PWM控制的比例环节系数；T_i积分环节系数；

D、外环控制参数辨识采用组合优化方法，其辨识流程：

(1)建立外环开环传递函数，需要辨识的参数为K_v、T_v、T_gv：

$W_{\mathrm{cv}}\left(s\right)=\frac{0.75K_v(T_{v}s+1)}{C{T}_v{s}^2(T_{\mathrm{gv}}s+1)}$

其中，K_V为比例环节系数，T_V、T_gV积分环节系数；

(2)外环辨识的采用光伏电站的实际运行数据，在进行数据的有效性分析、电气量程与工程量程的转换、标幺化处理后实际输入数据；

(3)将预处理后的光伏电站运行数据为输入，组合优化计算方法；

首先进行初步参数优化计算：建立目标函数方程如下式；通过迭代计算得到辨识参数K_V、T_V、T_gV的初步优化结果；

$F=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}\left[(M_i-S_i){\mathrm{\ω}}_i\right]$

然后进行精细化的参数优化计算：通过建立如下的目标函数方程式，并进行迭代计算，得到辨识参数K_V、T_V、T_gV的最终优化结果；

$F=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\left[(M_i-S_i){\mathrm{\ω}}_i\right]}^2$

其中，M_i是测量数据；S_i是仿真数据；ω_i是权值；

E、采用不同运行条件下的光伏电站运行数据和测试试验数据，重复A—D步骤，以实现对得到的模型内外环参数K_V、T_V、T_gV、K_iP、K_PWM及T_i进行校核，验证模型参数的泛化能力。

本发明的有益效果是，能够运用光伏电实际运行数据和对逆变器现场测试数据进行模型参数辨识，通过对光伏电站的模型参数辨识，还能够得到光伏组件实际最大功率运行点数据、逆变器关键控制参数，从而可使建立的光伏电站仿真模型具有较高的仿真精度，更接近光伏电站实际运行状态。

附图说明 图1为本发明光伏电站参数辨识流程图；

图2为外环参数辨识框图；

图3为光伏电站内环控制结构图；

图4为仿真软件中的光伏电站参数辨识模型图；

图5为简化的电流内环结构图。

具体实施方式：

一种光伏电站模型参数辨识方法，本发明特征在于，辨识方法步骤包括：

A、将光伏电站模型划分为外环控制和内环控制两个辨识环节，外环控制参数辨识，辨识光伏组件最大功率运行点参数、最大功率点跟踪(MPPT)控制参数；内环控制参数辨识，用于辨识逆变器动态跟踪控制参数；

B、对光伏电站输入源数据进行预处理：外环控制参数辨识时其输入数据来自于光伏电站的实际运行数据，应对运行数据进行筛选,去除其中的错误数据，然后将合格的数据添加到辨识数据库；内环控制参数辨识时其输入数据来自于对逆变器的现场施加阶跃信号得到响应曲线数据；

C、内环控制参数辨识采用最小二乘法，其辨识流程：

(1)建立内环开环传递函数：

$W_{\mathrm{ci}}\left(s\right)=\frac{K_{\mathrm{iP}}{K}_{\mathrm{PWM}}}{T_{i}L}\frac{T_{i}s+1}{s^2(1.5T_{s}s+1)}$

其中，K_iP为比例环节系数、K_PWM为PMW控制环节的增益，T_i、T_s积分环节系数；

(2)对现场施加阶跃信号得到的响应曲线数据进行相关转换处理，包括数据的有效性分析、电气量程与工程量程的转换、标幺化处理共三个步骤；

(3)将转换后的响应数据和阶跃输入信号为输入，采用最小二乘法计算方法，对内环传递函数中的参数K_iP、K_PWM及T_i进行辨识；

其中，K_iP为比例环节系数，K_PWM为PWM控制的比例环节系数；T_i积分环节系数；

D、外环控制参数辨识采用组合优化方法，其辨识流程：

(1)建立外环开环传递函数，需要辨识的参数为K_v、T_v、T_gv：

$W_{\mathrm{cv}}\left(s\right)=\frac{0.75K_v(T_{v}s+1)}{C{T}_v{s}^2(T_{\mathrm{gv}}s+1)}$

其中，K_V为比例环节系数，T_V、T_gV积分环节系数；

(2)外环辨识的采用光伏电站的实际运行数据，在进行数据的有效性分析、电气量程与工程量程的转换、标幺化处理后实际输入数据；

(3)将预处理后的光伏电站运行数据为输入，组合优化计算方法；

首先进行初步参数优化计算：建立目标函数方程如下式；通过迭代计算得到辨识参数K_V、T_V、T_gV的初步优化结果；

$F=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}\left[(M_i-S_i){\mathrm{\ω}}_i\right]$

然后进行精细化的参数优化计算：通过建立如下的目标函数方程式，并进行迭代计算，得到辨识参数K_V、T_V、T_gV的最终优化结果；

$F=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\left[(M_i-S_i){\mathrm{\ω}}_i\right]}^2$

其中，M_i是测量数据；S_i是仿真数据；ω_i是权值；

E、采用不同运行条件下的光伏电站运行数据和测试试验数据，重复A—D步骤，以实现对得到的模型内外环参数K_V、T_V、T_gV、K_iP、K_PWM及T_i进行校核，验证模型参数的泛化能力。

本发明光伏电站模型参数的辨识包括现场动态数据获取、数据预处理、光伏电站外环模型构建、光伏电站内环模型构建、采用优化算法进行参数辨识计算等步骤。

而在光伏电站模型参数辨识过程中，现场动态试验获得的数据，并不能直接用于模型参数辨识，需要结合模型结构和采用的辨识方法，进行相应的预处理，才能满足辨识需要。数据预处理在模型辨识中非常重要，也是模型辨识的基础工作。数据预处理情况比较复杂，处理的方法也比较多，往往必须针对具体的研究对象，根据一些经验性规则采用相应的处理方法。

本发明方法实现光伏电站参数辨识的主要实施技术操作要点包括：

1)数据的有效性分析

现场试验完成以后，需对采样数据进行有效性分析。有效性分析主要包括以下几个方面：测点响应时序是否正确；

测点测量值与实际正常值之间是否有较大差距；

测点是否安装正确；

如果试验数据经过分析，不能满足系统辨识的需要，且无法补救，应考虑再次进行试验，以达到系统辨识的要求。

2)电气量程与工程量程的转换

通过数据录波装置得到的现场数据，是用电压数值表示的电气量，需要根据电气量程与工程量程的对应关系，转换为与现场数据对应的工程量。

3)标幺化处理

在电力系统仿真分析中使用的是标幺值，因此需要将辨识数据进行标幺化处理，标幺化处理中的基准值如下：

标幺值等于实际值除以基准值；

频率的基准值为系统额定频率；

功率基准值为额定有功功率；

光照强度的基准值为1000W/m²。

4)初始工况的选择

在数据预处理中，需要选定用于系统参数辨识和模型校核的试验数据，对初始工况一般选取系统稳定运行，没有调整和干扰时的稳定工况作为初始工况，此时系统各参数间存在一定的稳态对应关系，因此这类数据在进行参数辨识和模型校核计算时更容易进行工作点的初始化(operation point)。

5)建立光伏电站参数辨识流程

如图1，为建立的参数辨识的流程图。在光伏电站参数辨识过程中通过不断的修改光伏电站或者逆变器参数，使目标函数达到最小，从而使仿真结果与实际数据的偏差达到最小，从而得到参数的最优结果。

6)光伏电站模型外环参数辨识

如图2所示，优化计算环节和是参数识别的关键部分，用于计算实测数据和仿真数据之间误差的目标函数，其中的参数p取值为1或2；优化计算执行器则收集所有的目标函数并最小化这些目标函数。

图2中，将预处理后的光伏电站运行数据为输入，组合优化计算方法。

首先进行初步参数优化计算：建立目标函数方程如下式。通过迭代计算得到辨识参数K_V、T_V、T_gV的初步优化结果。

$F=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}\left[(M_i-S_i){\mathrm{\ω}}_i\right]$

然后进行精细化的参数优化计算：通过建立如下的目标函数方程式，并进行迭代计算，得到辨识参数K_V、T_V、T_gV的最终优化结果。

$F=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\left[(M_i-S_i){\mathrm{\ω}}_i\right]}^2$

其中，M_i是测量数据；S_i是仿真数据；ω_i是权值；

7)光伏电站模型内环参数辨识

如图3所示，光伏电站内环控制因信号采样的延迟和PWM控制的小惯性特性，iq电流内环结构。其中，K_iP、K_il为比例积分环节系数、K_PWM为PMW控制环节的增益，T_s时间常数，R为阻抗值、L为感抗值。

图3中，Ts为电流内环电流采样周期(即亦为PWM开关周期)为简化分析，暂不考虑eq的扰动，且将PI调节器传递函数写成零极点形式，即

$\begin{array}{cc}{K}_{\mathrm{iP}}+\frac{K_{\mathrm{iI}}}{s}=K_{\mathrm{iP}}\frac{T_{i}s+1}{T_{i}s}& K_{\mathrm{iI}}=\frac{K_{\mathrm{iP}}}{T_i}\end{array}$

按II型系统设计，忽略VSR交流侧电阻，得简化的电流内环结构，如下所示。

得到开环传函：

$W_{\mathrm{ci}}\left(s\right)=\frac{K_{\mathrm{iP}}{K}_{\mathrm{PWM}}}{T_{i}L}\frac{T_{i}s+1}{s^2(1.5T_{s}s+1)}$

而内环控制参数辨识则建立的内环传递函数，采用最小二乘法，其辨识流程主要包括：

(1)对现场施加阶跃信号得到的响应曲线数据进行预处理，包括数据的有效性分析、电气量程与工程量程的转换、标幺化处理共三个步骤；

(2)将预处理后的响应数据和阶跃输入信号为输入，采用最小二乘法计算方法，对内环传递函数中的参数K_iP、K_PWM及T_i进行辨识；

8)参数辨识模型的搭建

如图4所示，构建了仿真软件中的光伏电站参数辨识模型。辨识主要考虑光伏逆变器的有功功率和无功功率特性，因此Measuerment file需要提供实际的逆变器有功和无功功率给comparator模块，再通过功率测量环节将仿真中光伏逆变器有功和无功功率特性提供给comparator模块，根据实际的数据和仿真结果计算目标函数值。

如果需要对逆变器的出口电压特性进行辨识，可以在该模型的基础之上多加一个电压测量环节，将逆变器出口电压提供给comparator模块。再通过Measurement file模块将实际的电压数据提供给comparator模块，即可。

同时，为了使仿真结果跟随实际的日照和温度变化，还需要将实际的日照和温度通过Measuerment file模块输入到光伏模块，让光伏模型在实际的日照和温度条件下运行。

Claims (1)

1.一种光伏电站模型参数辨识方法，其特征在于，辨识方法步骤包括：

A、将光伏电站模型划分为外环控制和内环控制两个辨识环节，外环控制参数辨识，辨识光伏组件最大功率运行点参数、最大功率点跟踪MPPT控制参数；内环控制参数辨识，用于辨识逆变器动态跟踪控制参数；

B、对光伏电站输入源数据进行预处理：外环控制参数辨识时其输入数据来自于光伏电站的实际运行数据，应对运行数据进行筛选,去除其中的错误数据，然后将合格的数据添加到辨识数据库；内环控制参数辨识时其输入数据来自于对逆变器现场施加阶跃信号得到响应曲线数据；

C、内环控制参数辨识采用最小二乘法，其辨识流程：

(1)建立内环开环传递函数：

$W_{\mathrm{ci}}\left(s\right)=\frac{K_{\mathrm{iP}}{K}_{\mathrm{PWM}}}{T_{i}L}\frac{T_{i}s+1}{s^2(1.5T_{s}s+1)}$

其中，K_iP为比例环节系数、K_PWM为PMW控制环节的增益，T_i、T_s积分环节系数；

(2)对现场施加阶跃信号得到的响应曲线数据进行相关转换处理，包括数据的有效性分析、电气量程与工程量程的转换、标幺化处理共三个步骤；

(3)将转换后的响应数据和阶跃输入信号为输入，采用最小二乘法计算方法，对内环传递函数中的参数K_iP、K_PWM及T_i进行辨识；

其中，K_iP为比例环节系数，K_PWM为PWM控制的比例环节系数；T_i积分环节系数；

D、外环控制参数辨识采用组合优化方法，其辨识流程：

(1)建立外环开环传递函数，需要辨识的参数为K_v、T_v、T_gv：

$W_{\mathrm{cv}}\left(\text{s}\right)=\frac{0.75K_v(T_{v}s+1)}{{\mathrm{CT}}_v{s}^2(T_{\mathrm{gv}}s+1)}$

其中，K_V为比例环节系数，T_V、T_gV积分环节系数；

(2)外环辨识的采用光伏电站的实际运行数据，在进行数据的有效性分析、电气量程与工程量程的转换、标幺化处理后实际输入数据；

(3)将预处理后的光伏电站运行数据为输入，组合优化计算方法；

首先进行初步参数优化计算：建立目标函数方程如下式；通过迭代计算得到辨识参数K_V、T_V、T_gV的初步优化结果；

$F=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}\left[(M_i-S_i){\mathrm{\ω}}_i\right]$

然后进行精细化的参数优化计算：通过建立如下的目标函数方程式，并进行迭代计算，得到辨识参数K_V、T_V、T_gV的最终优化结果；

$F=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\left[(M_i-S_i){\mathrm{\ω}}_i\right]}^2$

其中，M_i是测量数据；S_i是仿真数据；ω_i是权值；

E、采用不同运行条件下的光伏电站运行数据和测试试验数据，重复A—D步骤，以实现对得到的模型内外环参数K_V、T_V、T_gV、K_iP、K_PWM及T_i进行校核，验证模型参数的泛化能力。