一种用于次同步振荡仿真的直驱风机等值建模方法技术

本发明专利技术提出了一种用于次同步振荡仿真的直驱风机等值建模方法。本方法综合考虑了直驱永磁同步风电机组的控制器结构，根据实际风电场构造及风机参数建立仿真模型，在公共连接节点处施加一个次同步频率的测试电流，获得电压谐波各分量的幅值，根据谐波各分量幅值和控制器参数之间的关系，获得等值风电场的模型参数，实现对风电场的等值建模。

【技术实现步骤摘要】
一种用于次同步振荡仿真的直驱风机等值建模方法
本专利技术属于电力系统领域，特别地涉及一种用于次同步振荡仿真的直驱风机等值建模方法。
技术介绍
风电具有清洁、来源广泛等特点，近几年来发展迅猛。风能汇聚地区，地理上多处于偏远地区，多位于远离负荷中心的电网末端，因此我国风电场一般都采用集中并网再远距离输送的方式传输电能，由此引发的次同步振荡问题严重威胁着电力系统的安全稳定运行。电力系统的次同步振荡属于系统的振荡失稳，是一种由特殊的机电耦合引发的电力系统各部分间持续功率交互现象。传统火电机组引发的次同步振荡，指电力系统和发电机组以低于系统同步频率的某个或多个振荡频率交换显著的能量，振荡频率一般在8～49Hz。直驱风机作为现役的主要机型之一，其与传统火电机组不同，采用全功率式的逆变器作为并网接口，具有响应速度快，系统整体惯性小等特点。而直驱风机由于机械部分被交-直-交换流器隔离，其主导模态不参与系统振荡，而现有研究结果表明，网侧换流器控制的各个参数、锁相环、电网强度等都会影响次同步振荡的特性。同时直驱风机并网系统在发生次同步振荡时会有大量的分布有规律的谐波产生。实际系统中的风电场多由几十台甚至上百风电机组组成，同一地区通常有若干风电场接入。由于风机数目庞大，在分析系统稳定性时，构建系统特征方程会面临特征矩阵维数过高，统计参数所需时间过长等问题。因此，必须对风电场整体进行等值，即将同一风电场内的同型号风电机组等值成一台或几台风机。而现有的等值方法在选取参数时均是基于暂态值或经验值，缺少系统且有理论值支持的等值方法。因此，本专利技术提出一种用于次同步振荡仿真的直驱风机等值建模方法，通过分析直驱风机对谐波的动态响应过程，来获多机聚合后系统的整体参数，达到对风电场进行等值的目的。
技术实现思路
本专利技术提出一种用于次同步振荡仿真的直驱风机等值建模方法，包括以下步骤：步骤1：建立实际风电场的数字仿真平台；步骤2：在风电场出口公共连接点处施加谐波电流，测得谐波电压幅值；步骤3：根据谐波幅值和控制器参数之间的函数关系，获得等值后的参数。进一步，步骤3包括如下子步骤：1)建立系统小扰动分析模型，其包括直驱风机的锁相环、网侧换流器控制的小信号模型；2)根据形成的小信号模型，通过假定扰动输入的方式，逐步计算得出直驱风机dq坐标系下的电压、电流及电压指令值，最后得到风电场出口处各谐波分量的幅值、相位和控制器参数间的函数关系：Anj＝Anj(Kpp,Kip,Kpi,Kii,Kpdc,Kidc,L,C,δ,ωs)γnj＝γnj(Kpp,Kip,Kpi,Kii,Kpdc,Kidc,L,C,δ,ωs)(n＝1,2,3...；j＝1,2)其中Anj，γnj为各频率分量的幅值和相位，Kpp,Kii为锁相环比例、积分系数，Kpdc,Kidc为电压外环比例、积分系数，Kpi,Kii为电流内环比例、积分系数，C,L,Lg为直流侧电容、网侧换流器连接电感、电网等效电感，δ为正弦波扰动幅值，ωs为正弦波扰动频率，x表示各谐波分量幅值、相位和参数之间的函数；3)控制器参数和谐波分量幅值间函数关系的函数形式可以表示成：Kpp＝Kpp(Anj,Kip,Kpi,Kii,Kpdc,Kidc...L,C,δ,ωs)...C＝C(Kpp,Kip,Kpi,Kii,Kpdc,Kidc...L,Anj,δ,ωs)(n＝1,2,3...；j＝1,2)该式为子步骤2)中的函数的反函数，利用matlab中solve函数求解模型参数。附图说明：图1为针对次同步振荡问题的直驱风机等值原理图；图2为直驱风机并网原理图；图3为网侧换流器dq坐标下控制系统框图；图4为锁相环控制系统框图；图5为本专利技术实施例中等值前风电场公共连接点电压频谱；图6为本专利技术实施例中等值后风电场公共连接点电压频谱。具体实施方式：本专利技术提出一种用于次同步振荡仿真的直驱风机等值建模方法。首先需要建立风电场的实时仿真平台，然后在公共连接点处施加谐波电流，测得谐波电压幅值，根据谐波幅值和控制器参数之间的函数关系获得等值后的参数。图1为针对次同步振荡问题的直驱风机等值原理图。包括建立直驱风机控制系统小信号模型、输出谐波的幅值与控制器参数间函数关系、反函数计算得到等值后的控制器参数、得到等值后的单机模型；在风电场出口公共连接点处叠加测试信号电流谐波，提取谐波幅值。建立直驱风机控制系统小信号模型的如下：直驱风机对谐波响应过程：直驱风机的并网结构及网侧换流器控制结构如图2所示，网侧换流器dq坐标下控制系统如图3所示，锁相环控制系统结构如图4所示。其中αβ表示换流器控制系统的坐标系，为网侧换流器端口电压值，uα,uβ为换流器输出电压值，uαref,uβref为电压指令值，iα,iβ为流入换流器电流值，iαref,iβref为电流指令值，udc,udc_ref为直流电压、直流电压指令值，Kpdc,Kidc为电压外环比例、积分系数，Kpi,Kii为电流内环比例、积分系数，C,L,Lg为直流侧电容、网侧换流器连接电感、电网等效电感。(1)锁相环小信号方程系统同步旋转角度为θ，锁相环输出相位为θc：其中θc＝-uq/Ut0(2)Ut0为网侧电压幅值，直驱风机网侧换流器采取电压定向的控制策略，在图4控制方式下可以认为Ut0＝Ud0＝1，则式(2)线性化后有：Δθc＝-Δuq/Ut0+Δθ(3)将(1)线性化后代入(3)中可以得到：公式(4)即为锁相环对次同步频率谐波的响应过程。(2)控制系统的小信号方程由图3得到直驱风机控制器同步坐标系下的小信号方程为：两相同步旋转坐标下的电气量由三相静止坐标下的电气量经过派克变换得到，锁相环在次同步谐波影响下会对这一过程产生影响，考虑谐波分量干扰的情况下控制系统方程为：其中输出谐波的幅值与控制器参数间函数关系。控制器及锁相环对次同步频率分量的响应过程。假设从风机出口处的公共连接点节点传入网侧一个扰动电流，为幅值δ，频率ωs的三相对称正弦波，在后文的推导过程中发现扰动电流的初始相位不影响控制系统的动作过程，不失一般性的可以假设其初始相位为零，A相电流可以表示为：Δia＝δcos(ωst)(8)由于网络电阻较小，相比电抗通常可以忽略，(8)引起的A相电压波动为：将(8)(9)变换到dq坐标系下可以得到：将(11)代入(4)中得到：将(11)(12)代入(6)中得到电压指值：其中Mij为各分量幅值，γij为各分量的初始相位(i,j＝1,2,3)。Δudref，Δuqref与Δud，Δuq叠加作为新的扰动。这些新的扰动分量继续作用于换流器的控制系统。按照式(8)到式(13)的过程，多次迭代后得到三相电压的最终表达式：其中Anj为各分量幅值，γnj为各分量的初始相位(n＝1,2,3...，i＝1,2)。从上述推导可知Anj、γnj为有关直驱风机控制器参数的函数，将其表示成函数形式为：反函数计算得到等值后的控制器参数、得到等值后的单机模型。(1)直驱风机出口处电压谐波的提取如图2所示，在直驱风机出口处注入一个小幅值的电流测试信号Δia＝δcos(ωst)，此时风电场出口处电压具有与式(14)相同的形式：其中Ani可以由傅里叶变换获得，显然各次谐波幅值的获取较相位容易。(2)等值后参数的获取式(15)的反函数形式为：求解本文档来自技高网...

【技术保护点】
1.一种用于次同步振荡仿真的直驱风机等值建模方法，包括以下步骤：步骤1：建立实际风电场的数字仿真平台；步骤2：在风电场出口公共连接点处施加谐波电流，测得谐波电压幅值；步骤3：根据谐波幅值和控制器参数之间的函数关系，获得等值后的参数。

【技术特征摘要】
1.一种用于次同步振荡仿真的直驱风机等值建模方法，包括以下步骤：步骤1：建立实际风电场的数字仿真平台；步骤2：在风电场出口公共连接点处施加谐波电流，测得谐波电压幅值；步骤3：根据谐波幅值和控制器参数之间的函数关系，获得等值后的参数。2.权利要求1所述的一种用于次同步振荡仿真的直驱风机等值建模方法，其特征在于，所述步骤3包括如下子步骤：1)建立系统小扰动分析模型，其包括直驱风机的锁相环、网侧换流器控制的小信号模型；2)根据形成的小信号模型，通过假定扰动输入的方式，逐步计算得出直驱风机dq坐标系下的电压、电流及电压指令值，最后得到风电场出口处各谐波分量的幅值、相位和控制器参数间的函数关系：Anj＝Anj(Kpp,Kip,Kpi,Kii,Kpdc,Kidc,L,C,δ,ωs)γnj＝γnj(Kpp,Kip,K...

【专利技术属性】
技术研发人员：徐衍会，曹宇平，
申请(专利权)人：华北电力大学，
类型：发明
国别省市：北京,11