CN113872190A - 一种提高小信号系统计算效率的降维方法 - Google Patents

Abstract

本发明涉及一种提高小信号系统计算效率的降维方法，包括以下步骤：1)建立永磁直驱风机经柔性直流输电并网系统的详细模型；2)通过特征值分析法对详细模型进行求解，得到状态矩阵的特征值；3)将特征值按照对系统的影响程度大小进行划分，按照时间尺度的方法剔除快速衰减的变量并保留缓慢衰减的变量；4)建立永磁直驱风机经柔性直流输电并网系统的降维模型；5)将详细模型的动态响应曲线与降维模型的动态响应曲线进行对比，验证降维模型的可靠性、一致性与相似性。与现有技术相比，本发明具有显著降低模型维度、提升计算效率、准确性高、实用性广等优点。

Description

一种提高小信号系统计算效率的降维方法

技术领域

本发明涉及电力系统稳定性分析领域，尤其是涉及一种提高小信号系统计算效率的降维方法。

背景技术

面对风电并网的大力发展，大量的风电机组并网的场景，往往采用风电场模型聚合的方法，在电网故障的情况下，将动态特性相近的风电机组划分为一个机群或者聚合为一台风电机组，但是当对大量风机采用详细的数学模型来进行稳定性分析时，由于单台风电机组的构成较为复杂，将会大大增加模型的复杂程度，计算时间加长，占用内存庞大。

目前，国内外学者针对电力系统模型数据降维方法开展了广泛研究。

基于电力系统时间尺度特性的降维方法在传统交流电力系统中发展较为成熟，但在交直流电力系统中应用尚不多见。有学者早在20世纪80年代初就注意到电力系统的多时间尺度特性，并研究了多时间尺度分解的迭代算法。有学者在电力网络为RL电路的假设下推得了互联电力系统的双时间尺度模型，并运用奇异摄动理论研究了同步发电机模型的降维方法，但未考虑负荷的时间尺度特性及负荷模型的降维。有学者首次建立了三时间尺度的电力系统奇异摄动模型，给出了忽略快动态和固定慢动态的降维条件，并对纯交流系统进行了降维研究，但未计算所提出的降维条件成立的奇异摄动参数范围，对于时间尺度较多、且所提取的奇异摄动参数无法准确描述状态变量变化快慢程度的系统而言，容易造成降维前后系统稳定性的不一致。

发明内容

本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种提高小信号系统计算效率的降维方法。

本发明的目的可以通过以下技术方案来实现：

一种提高小信号系统计算效率的降维方法，包括以下步骤：

1)建立永磁直驱风机经柔性直流输电并网系统的详细模型；

2)通过特征值分析法对详细模型进行求解，得到状态矩阵的特征值；

3)将特征值按照对系统的影响程度大小进行划分，按照时间尺度的方法剔除快速衰减的变量并保留缓慢衰减的变量；

4)建立永磁直驱风机经柔性直流输电并网系统的降维模型；

5)将详细模型的动态响应曲线与降维模型的动态响应曲线进行对比，验证降维模型的可靠性、一致性与相似性。

所述的永磁直驱风机经柔性直流输电并网系统为一小信号系统。

所述的步骤1)中，根据永磁直驱风机经柔性直流输电并网系统的拓扑结构构建详细模型，该拓扑结构依次包括直驱风电场、永磁同步发电机、换流器系统、VSC-HVDC以及受端电网，所述的换流器系统由全控整流器、直流稳压电容与全控逆变器三部分组成。

根据永磁直驱风机经柔性直流输电并网系统拓扑结构的动态模型，确定状态变量X为：

其中，X_PMSG为D-PMSG模型的状态变量，X_plls为风机网侧锁相环状态变量，X_PMSG-Ctrl为D-PMSG控制器参数状态变量，X_plls1为VSC-HVDC送端锁相环状态变量，X_HVDC为VSC-HVDC模型的状态变量，X_plls2为VSC-HVDC受端锁相环状态变量，X_HVDC-Ctrl为VSC-HVDC模型控制器的状态变量，Δ表示微分量，ω₁为风力机转子电气转速，ω₂为发电机转子电气转速，δ₁、δ₂分别为风力机于发电机转子相对于额定电气转速下同步旋转参考轴的电气角位移，i_d为PMSG的d轴电流，i_q为PMSG的q轴电流，u_dc为直流电容端电压，i_sd为网侧控制器d轴电流，i_sq为网侧控制器q轴电流，θ_plls为风机网侧锁相环的输出相位，Z_plls为风机网侧锁相环的角度瞬时变化值的负数，w₁、w₂、w₃别为机侧控制器状态变量，x₁、x₂、x₃别为网侧控制器状态变量，θ_plls1为VSC-HVDC送端锁相环的输出相位，Z_plls1为VSC-HVDC送端锁相环的角度瞬时变化值的负数，u_dc1为送端控制器直流电容端电压，i_dc为直流电容电流，u_dc2为受端控制器直流端电容电压，i_s2d为受端输出d轴电流，i_s2q为受端输出q轴电流，θ_plls2,为受端锁相环的输出相位，Z_plls2为受端锁相环的角度瞬时变化值的负数，y₁、y₂、y₃、y₄分别为送端控制器状态变量，z₁、z₂、z₃、z₄别为受端控制器状态变量。

根据永磁直驱风机经柔性直流输电并网系统拓扑结构的动态模型，确定非状态变量Y为：

其中，Y₁为PMSG代数变量，Y₂为交流环节代数变量，Y₃为受端电网代数变量，T₁为作用于轴系风力机质量块的机械转矩，T₂为作用于轴系发电机质量块的电磁转矩的负数，u_d、u_q分别为PMSG机端的d轴和q轴电压，u_cd、u_cq为网侧交流坐标系下的PMSG端d轴与q轴电压，u_s为电网电压幅值，u_s1为送端电网电压幅值，u_c1d、u_c1q分别为送端交流坐标系下的PMSG端d轴与q轴电压，u_c2d、u_c2q分别为受端电网交流坐标系下d轴与q轴电压，u_s2d、u_s2q分别为受端交流坐标系下d轴与q轴电压。

所述的永磁直驱风机经柔性直流输电并网系统的详细模型的表达式为

其中，p为微分算子，A₁为状态变量系数矩阵，B₁为代数变量系数矩阵，C₁为代数变量转化矩阵。

所述的步骤2)中，状态矩阵A的表达式为：

A＝A₁+B₁C₁。

所述的步骤3)具体为：

获取状态矩阵各特征值对应的阻尼比，并设定阈值，阻尼比超过阈值的状态变量作为快速衰减的变量被剔除，反之，阻尼比未超过阈值的状态变量作为缓慢衰减的变量被保留，以此完成降维。

所述的阻尼比的计算式为：

其中，α_i、ω_i分别为第i个特征值λ_i对应的实部和虚部，且λ_i＝α_i±jω_i(i＝1,2,…,n)，当ξ_i≥10％时表明系统阻尼较强，而ξ_i≤5％时表明系统阻尼较弱，当ξ_i≤0时表明系统阻尼为负，出现增幅振荡。

降维模型中主导状态变量的表达式为：

X_r＝[Δi_q,Δi_sq,Δw₃,Δx₃,Δi_s2d,Δi_s2q,Δy₃,Δy₄,Δz₂]

pX_r＝A_r(s)X_r

其中，p为微分算子，i_q为PMSG的q轴电流，i_sq为网侧控制器q轴电流，X_r为主导状态变量，w₃为机侧控制器状态变量，x₃为网侧控制器状态变量，i_s2d为受端输出d轴电流，i_s2q为受端输出q轴电流，y₃、y₄为送端控制器状态变量，z₂为受端控制器状态变量，s为拉普拉斯算子，A_r(s)为运算形式的降维系数矩阵。

与现有技术相比，本发明具有以下优点：

本发明从状态变量衰减速度方面来对高维模型的状态变量进行筛选，不但使高维模型的维度显著降低，提升了计算效率，而且降维模型的准确性仍得到了保证，在工程中面对大量风电场并网采用详细数学模型建模时，具有广泛的实用性。

附图说明

图1为本发明直驱风机经柔性直流输电并网系统的拓扑结构图。

图2为D-PMSG经VSC-HVDC并网系统的拓扑结构图。

图3为机侧控制器模型。

图4为网侧控制器模型。

图5为送端控制器模型。

图6为受端控制器模型。

图7为d₁-q₁旋转坐标系与d₂-q₂旋转坐标系的关系图。

图8为D-PMSG与VSC-HVDC之间的简化交流输电线路模型。

图9为本发明主导状态变量的动态响应曲线，其中，图(9a)为状态变量Δi_q的动态响应曲线，图(9b)为状态变量Δi_sq的动态响应曲线，图(9c)为状态变量Δw₃的动态响应曲线，图(9d)为状态变量Δx₃的动态响应曲线，图(9e)为状态变量Δi_s2d的动态响应曲线，图(9f)为状态变量Δi_s2q的动态响应曲线，图(9g)为状态变量Δy₃的动态响应曲线，图(9h)为状态变量Δy₄的动态响应曲线，图(9i)为状态变量Δz₂的动态响应曲线。

具体实施方式

本发明提供一种提高小信号系统计算效率的降维方法，该方法主要包括以下步骤：

S1、建立永磁直驱风机经柔性直流输电并网的详细模型；

S2、通过特征值分析法对详细模型进行分析求解；

S3、将特征值按照对系统的影响程度大小进行划分，按照时间尺度的方法略去快速衰减的变量并保留缓慢衰减的变量；

S4、建立详细系统的降维模型；

S5、将详细模型的动态响应曲线与降维模型的动态响应曲线进行对比，验证降维模型的可靠性、一致性与相似性。

其主要原理介绍如下：

本发明通过建立系统的详细小信号模型，并对其进行特征值分析，筛选出所需要观测的特征向量，再通过参与因子的大小排序，略去快速衰减的量以筛选出需要保留的特征向量以达到降维的目的。

在小信号系统分析过程中，通常可以通过状态变量的变化而将其分为两类：快动态量与慢动态量，因此可以采用两个或者多个尺度来进行渐近展开求解，可以将快速衰减的量剔除，以降低模型的复杂程度。

假设一个系统线性化之后可由下式描述：

式中，Δx为状态变量，A为状态空间矩阵，B为输入矩阵，u为输入变量，将系统所有特征向量列出，略去其中快速衰减的模式，只保留主要特征值，即缓慢衰减的模式，则上式可写成：

式中：Δx_r为主导特征值(衰减缓慢的变量)；Δx₂为非主导特征值，消去变量Δx₂得到降维系统的状态方程：

式中：s为拉普拉斯算子；A_r(s)为运算形式的降维系数矩阵。

综上可得，通过时间尺度的筛选，在一定程度上减少了模型的维度，降低了减少了系统的计算量

实施例

在本发明中，目标在于通过时间尺度，达到降低详细系统模型维度而又不失去详细模型主要特征，以便计算更加简便快捷。该方法在建立直驱风机经柔性直流输电方式并网的详细状态空间模型的基础上，采用特征值分析法对模型所产生的振荡模式进行分析，得出系统产生的所有振荡模式。基于该详细模型的分析结果，通过对比各振荡模式的阻尼比，选择出对系统影响较大的振荡模式，并基于时间尺度略去快速衰减的模式，可建立降维系统的状态方程，降低了模型的维度，详细系统模型拓扑图如图1所示。

图1具体为D-PMSG经VSC-HVDC并网系统的拓扑结构，本例以1台D-PMSG代表直驱风电场，风力机直接与永磁同步发电机(Permanent Magnet Synchronous Generator，PMSG)相连，通过控制器系统、VSC-HVDC连接至受端电网，其中换流系统由全控整流器、直流稳压电容与全控逆变器三部分组成，除了机侧控制器，其他控制器均采用基于电网电压的空间矢量控制策略。

1、D-PMSG的动态模型

由于D-PMSG的转速较低、轴系刚性相对不足等特点，使得其轴系的自然扭振频率较低，因此，相比于单质量块模型，如图2所示，双质量块模型更能精确地描述D-PMSG轴系的机械特性。

双质量模块的轴系线性化模型为：

其中，ΔT1＝ΔTm，ΔT2＝-ΔTe，p为微分算子，T_J1、T_J2分别为风力机和发电机转子的惯性时间常数；ω₁、ω₂分别为风力机和发电机转子的电气转速，且ω₁＝n_pω_w，ω₂＝n_pω_g，ω_w、ω_g分别为风力机和发电机转子的机械转速，np为PMSG的极对数，在标幺制中，机械转速数值上等于电气转速；ω_b为基准电气转速，取ω_b＝ω_N，ω_N为风力机的额定电气转速；K₁₂为轴系的刚度系数；D₁₁，D₂₂分别为风力机与发电机的自阻尼系数，D₁₂为风力机与发电机之间的互阻尼系数；δ₁、δ₂分别为风力机与发电机转子相对于额定电气转速下同步旋转参考轴的电气角位移，Tm为作用于轴系风力机质量块的机械转矩，Te为作用于轴系发电机质量块的电磁转矩。

2、永磁发电机模型

在建立PMSG的模型时，采用发电机惯例，则经过等量派克变换，在正序d-q旋转坐标系下PMSG的数学模型如下。

式中，u_d和u_q分别为PMSG机端的d轴和q轴电压；i_d和i_q分别为PMSG的d轴和q轴电流；R_s为定子电阻；ω为发电机转子的电气角速度，且ω＝ω₂，ψ_d与ψ_q分别为d轴与q轴PMSG各绕组的磁链，ψ_f为永磁体磁链，L_d和L_q分别是PMSG的d轴和q轴电感，ω₂₀为发电机转子的电气转速初值，i_d0为PMSG的d轴电流初值，i_q0为PMSG的q轴电流，ψ_f0为永磁体磁链。

3、机侧控制器模型

机侧控制器的模型即为PMSG的电压方程，机侧控制器的控制器即MSC通过控制D-PMSG输出有功功率以实现最佳风能跟踪控制，如图3所示，采用i_{d_ref}＝0控制策略，则可以得到MSC直流环节的动态模型为：

其中，w₁、w₂、w₃均为机侧控制器状态变量，k_i1、k_i2均为机侧控制器积分系数，k_p1、k_p2、k_p3均为机侧控制器比例系数。

4、网侧控制器模型

对于网侧控制器的控制器，即GSC，常采用基于电网电压定向的矢量控制策略实现直流电压控制和并网无功控制，如图4所示，则GSC的线性化动态模型为：

其中，x₁、x₂、x₃均为网侧控制器状态变量，h_p1、h_p2、h_p3均为网侧控制器比例系数，h_i1、h_i2、h_i3均为网侧控制器比例系数，i_sd与i_sq分别为网侧交流坐标系下的d轴与q轴分量，i_sd0与i_sq0分别为网侧交流坐标系下的d轴与q轴分量初值，θ_plls为受端锁相环的输出相位，Q_{s_ref}为无功初值，

为网侧交流坐标系下的d轴电压，

为网侧交流坐标系下的d轴电压初值，ω_s0为网侧控制器交流侧的电网角频率的，u_s为电网电压幅值，d₁与d₂为系数，u_sd0为电网d轴电压初值，L为网侧换流器并网电感。

5、直流环节模型

直流环节的动态模型为：

C_dcu_dcpu_dc＝u_di_d+u_qi_q-(u_cdi_sd+u_cqi_sq)

其中，C_dc是直流电容；u_dc是直流电容端电压。

对其进行线性化，则可得直流环节的线性化动态模型：

其中，u_d0为PMSG端d轴电压，u_q0为PMSG端q轴电压，p₀为有功初值，u_dc0为直流电容端电压初值，u_cd、u_cq为网侧交流坐标系下的PMSG端d轴与q轴电压，u_cd0、u_cq0为网侧交流坐标系下的PMSG端d轴与q轴电压初值，i_d0为PMSG端d轴电流初值，i_q0为PMSG端q轴电流初值。

6、送端控制器模型

VSC-HVDC的送端控制器采用基于电网电压定向的空间矢量控制的定有功控制和定无功控制，其控制模型如图5所示，对于送端控制器的控制器，即送端控制器(SEC)，忽略Δω_s1，则其线性化动态模型可表示为：

其中，y₁、y₂、y₃、y₄均为送端控制器的状态变量，u_s1d0为，i_s1d、i_s1q分别为送端控制器交流侧d轴与q轴电流分量，i_s1d0、i_s1q0分别为送端控制器交流侧d轴与q轴电流分量初值，u_s1为送端电网电压幅值，b_p1、b_p2、b_p3、b_p4均为送端控制器比例系数，b_i1、b_i2、b_i3、b_i4均为送端控制器积分系数，θ_plls1为风机网侧锁相环的输出相位，u_c1d、u_c1q分别为送端交流坐标系下的PMSG端d轴与q轴电压，ω_s10为送端环流其交流测电网角频率初值，L₁为送端控制器逆变侧等效电感。

7、受端控制器模型

受端控制器采用基于电网电压定向的空间矢量控制实现定直流电压与定无功功率控制策略，用来保证VSC-HVDC的稳定运行，其控制模型如图6所示，对于受端控制器的控制器，即受端控制器(REC)，忽略Δω_s2，则其线性化动态模型可表示为：

其中，z₁、z₂、z₃、z₄分别为受端控制器状态变量，u_c2d、u_c2q分别为受端电网交流坐标系下d轴与q轴电压，i_s2d、i_s2q分别为受端交流坐标系下的d轴与q轴电流，i_s2d0、i_s2q0分别为受端交流坐标系下的d轴与q轴电流初值，c_i1、c_i2、c_i3、c_i4分别为受端控制器积分系数，c_p1、c_p2、c_p3、c_p4分别为受端控制器比例系数，u_s2d、u_s2q分别为受端交流坐标系下d轴与q轴电压，u_s2d0、u_s2q0分别为受端交流坐标系下d轴与q轴电压初值，θ_plls2为受端锁相环的输出相位，θ_plls20为受端锁相环的输出相位初值，h₁、h₂、e₁、e₂、e₃、e₄、e₅、e₆、d₅、d₆均为系数。

8、D-PMSG与VSC-HVDC之间的接口动态模型

假设有两个正序同步旋转d-q坐标系，分别为d₁-q₁坐标系与d₂-q₂坐标系，d₁轴与d₂轴之间的夹角为θ，两个旋转坐标系之间的变化如图7所示。

忽略交流输电线路的电容效应，则D-PMSG与VSC-HVDC之间的交流输电线路模型可以简化为如图8所示的形式。在交流输电线路中，用电阻R_l表示交流输电线路的等值电阻，电感L_l表示交流输电线路的等值电感。

将基于d₁-q₁旋转坐标系的VSC-HVDC送端变流器及其控制器模型中的物理量变换到d₂-q₂旋转坐标系，则进一步得到接口动态方程，即交流输电线路的动态方程：

其中，u'_s1d、u'_s1q分别为经过坐标变换后的送端控制器交流侧d轴与q轴电压分量，

为电网基准角频率的倒数，R_l为交流输电线路电阻，L_l为交流输电线路电感，L₁为送端控制器逆变侧等效电感。

9、全系统的线性化动态模型

线性化动态模型

D-PMSG经VSC-HVDC并网系统的状态变量为：

根据上述所建的动态模型，可得：

状态矩阵A中含有各个控制器参数，因此，通过调节控制器参数，可以影响状态矩阵A的元素，进而影响状态矩阵A的特征值。

pX＝AX

A＝A₁+B₁C₁

本发明以直驱风机经柔性直流输电方式并网系统为例进行仿真分析，验证本发明提出的降维方法的准确性。采用6台2MW、0.69kV的永磁直驱风电机组，并对其进行等值。风速设定为12m/s，桨距角为0°，风机的直流电压参考值设为1.0(pu)，网侧控制器的无功功率参考值设为0，机侧控制器的参考电流设为0。对所建的系统模型在交流输电距离为100km的情况下进行特征值分析。

特征值分析结果如表1所示。

表1特征值计算结果

其中，λ₂₃,λ₂₄代表低频振荡LFO-1，λ₂₅,λ₂₆,λ₂₇,λ₂₈,λ₂₉,λ₃₀分别代表次同步振荡SSO-1、SSO-2、SSO-3，λ₃₁,λ₃₂为超次同步振荡模式SupSO,λ₃₃,λ₃₄为高频振荡模式HFO，其他模式阻尼较高，属于快速衰减模式，可以略去。

LFO模式的阻尼比为3.73％，小于其稳定所需的阻尼比阈值5％；在次/超同步频段内，SSO-1与SSO-2模式的阻尼比小于0；SSO-3、SupSO与HFO模式的阻尼比均较小，但考虑到振荡频率越高，稳定所要求的阻尼比则越低，而SSO-3、SupSO与HFO阻尼比分别达到了1.91％、0.53％与3.07％。

选取观测的振荡模式是SSO-1、SSO-2、SSO-3与SupSO模式。

表2各振荡模式对应的强相关量和子系统

通过特征值与参与因子的计算，分析结果表明，模型的主导状态变量为Δiq、Δisq、Δw3、Δx3、Δis2d、Δis2q、Δy3、Δy4、Δz2，因此系统模型可降为9阶模型，其计算量相较于原矩阵34阶已大大减少，因此，降维模型中主导状态变量的表达式为：

X_r＝[Δi_q,Δi_sq,Δw₃,Δx₃,Δi_s2d,Δi_s2q,Δy₃,Δy₄,Δz₂]

pX_r＝A_r(s)X_r

系统降维前后主导状态变量的动态响应曲线如图9所示。从图9可见，状态变量降维前后的动态响应曲线虽具有误差，但误差较小，却具有很好的相似性、一致性。

综上，本发明提出了一种提高小信号系统计算效率的降维方法，其目标在于通过时间尺度，达到降低详细系统模型维度而又不失去详细模型主要特征，以便计算更加简便快捷。该方法在建立直驱风机经柔性直流输电方式并网的详细状态空间模型的基础上，采用特征值分析法对模型所产生的振荡模式进行分析，得出系统产生的所有振荡模式。基于该详细模型的分析结果，通过对比各振荡模式的阻尼比，选择出对系统影响较大的振荡模式，并基于时间尺度略去快速衰减的模式，可建立降维系统的状态方程，降低了模型的维度。最后，通过对详细模型以及降维模型的动态响应曲线的对比，验证了降维模型的相似性与一致性。

Claims (10)

1.一种提高小信号系统计算效率的降维方法，其特征在于，包括以下步骤：

1)建立永磁直驱风机经柔性直流输电并网系统的详细模型；

2)通过特征值分析法对详细模型进行求解，得到状态矩阵的特征值；

3)将特征值按照对系统的影响程度大小进行划分，按照时间尺度的方法剔除快速衰减的变量并保留缓慢衰减的变量；

4)建立永磁直驱风机经柔性直流输电并网系统的降维模型；

5)将详细模型的动态响应曲线与降维模型的动态响应曲线进行对比，验证降维模型的可靠性、一致性与相似性。

2.根据权利要求1所述的一种提高小信号系统计算效率的降维方法，其特征在于，所述的永磁直驱风机经柔性直流输电并网系统为一小信号系统。

3.根据权利要求1所述的一种提高小信号系统计算效率的降维方法，其特征在于，所述的步骤1)中，根据永磁直驱风机经柔性直流输电并网系统的拓扑结构构建详细模型，该拓扑结构依次包括直驱风电场、永磁同步发电机、换流器系统、VSC-HVDC以及受端电网，所述的换流器系统由全控整流器、直流稳压电容与全控逆变器三部分组成。

4.根据权利要求1所述的一种提高小信号系统计算效率的降维方法，其特征在于，根据永磁直驱风机经柔性直流输电并网系统拓扑结构的动态模型，确定状态变量X为：

其中，X_PMSG为D-PMSG模型的状态变量，X_plls为风机网侧锁相环状态变量，X_PMSG-Ctrl为D-PMSG控制器参数状态变量，X_plls1为VSC-HVDC送端锁相环状态变量，X_HVDC为VSC-HVDC模型的状态变量，X_plls2为VSC-HVDC受端锁相环状态变量，X_HVDC-Ctrl为VSC-HVDC模型控制器的状态变量，Δ表示微分量，ω₁为风力机转子电气转速，ω₂为发电机转子电气转速，δ₁、δ₂分别为风力机于发电机转子相对于额定电气转速下同步旋转参考轴的电气角位移，i_d为PMSG的d轴电流，i_q为PMSG的q轴电流，u_dc为直流电容端电压，i_sd为网侧控制器d轴电流，i_sq为网侧控制器q轴电流，θ_plls为风机网侧锁相环的输出相位，Z_plls为风机网侧锁相环的角度瞬时变化值的负数，w₁、w₂、w₃别为机侧控制器状态变量，x₁、x₂、x₃别为网侧控制器状态变量，θ_plls1为VSC-HVDC送端锁相环的输出相位，Z_plls1为VSC-HVDC送端锁相环的角度瞬时变化值的负数，u_dc1为送端控制器直流电容端电压，i_dc为直流电容电流，u_dc2为受端控制器直流端电容电压，i_s2d为受端输出d轴电流，i_s2q为受端输出q轴电流，θ_plls2,为受端锁相环的输出相位，Z_plls2为受端锁相环的角度瞬时变化值的负数，y₁、y₂、y₃、y₄分别为送端控制器状态变量，z₁、z₂、z₃、z₄别为受端控制器状态变量。

5.根据权利要求4所述的一种提高小信号系统计算效率的降维方法，其特征在于，根据永磁直驱风机经柔性直流输电并网系统拓扑结构的动态模型，确定非状态变量Y为：

其中，Y₁为PMSG代数变量，Y₂为交流环节代数变量，Y₃为受端电网代数变量，T₁为作用于轴系风力机质量块的机械转矩，T₂为作用于轴系发电机质量块的电磁转矩的负数，u_d、u_q分别为PMSG机端的d轴和q轴电压，u_cd、u_cq为网侧交流坐标系下的PMSG端d轴与q轴电压，u_s为电网电压幅值，u_s1为送端电网电压幅值，u_c1d、u_c1q分别为送端交流坐标系下的PMSG端d轴与q轴电压，u_c2d、u_c2q分别为受端电网交流坐标系下d轴与q轴电压，u_s2d、u_s2q分别为受端交流坐标系下d轴与q轴电压。

6.根据权利要求5所述的一种提高小信号系统计算效率的降维方法，其特征在于，所述的永磁直驱风机经柔性直流输电并网系统的详细模型的表达式为

其中，p为微分算子，A₁为状态变量系数矩阵，B₁为代数变量系数矩阵，C₁为代数变量转化矩阵。

7.根据权利要求1所述的一种提高小信号系统计算效率的降维方法，其特征在于，所述的步骤2)中，状态矩阵A的表达式为：

A＝A₁+B₁C₁。

8.根据权利要求1所述的一种提高小信号系统计算效率的降维方法，其特征在于，所述的步骤3)具体为：

获取状态矩阵各特征值对应的阻尼比，并设定阈值，阻尼比超过阈值的状态变量作为快速衰减的变量被剔除，反之，阻尼比未超过阈值的状态变量作为缓慢衰减的变量被保留，以此完成降维。

9.根据权利要求8所述的一种提高小信号系统计算效率的降维方法，其特征在于，所述的阻尼比的计算式为：

其中，α_i、ω_i分别为第i个特征值λ_i对应的实部和虚部，且λ_i＝α_i±jω_i(i＝1,2,…,n)，当ξ_i≥10％时表明系统阻尼较强，而ξ_i≤5％时表明系统阻尼较弱，当ξ_i≤0时表明系统阻尼为负，出现增幅振荡。

10.根据权利要求1所述的一种提高小信号系统计算效率的降维方法，其特征在于，所述的步骤4)中，降维模型中主导状态变量的表达式为：

X_r＝[Δi_q,Δi_sq,Δw₃,Δx₃,Δi_s2d,Δi_s2q,Δy₃,Δy₄,Δz₂]

pX_r＝A_r(s)X_r

其中，p为微分算子，i_q为PMSG的q轴电流，i_sq为网侧控制器q轴电流，X_r为主导状态变量，w₃为机侧控制器状态变量，x₃为网侧控制器状态变量，i_s2d为受端输出d轴电流，i_s2q为受端输出q轴电流，y₃、y₄为送端控制器状态变量，z₂为受端控制器状态变量，s为拉普拉斯算子，A_r(s)为运算形式的降维系数矩阵。

Images