本发明涉及新能源发电技术中的风电场等值建模技术领域，具体涉及一种基于相似度同调的双馈风电机动态等值方法。
背景技术：
随着现代电力工业的迅速发展，大规模交直流互联电力系统出现，多机(上千台机)、大电网(几千条线路、几千条母线)特征越来越突出，这使得电力系统规划设计和运行方式的计算变得异常的繁琐复杂。一方面，详细完成这些计算需要耗费大量的时间和空间，另一方面，受到仿真规模和硬件设备的限制，上述计算难以在实际中实现。为此需要对原始系统进行某种合理的简化，从而减小系统规模。实际中，对一个大电力系统的动态性能研究一般只对某一个区域最感兴趣，该区域称之为研究系统，而与此区域较远的区域，研究中只要计及其对研究区域的影响，其内部不必详细描写，往往可做降阶及简化，以节省研究的人力和物力，该拟作简化的区域称为外部系统。把不需要详细分析的外部系统进行等值化简十分必要。
风电作为最具发展潜力和规模化开发的新能源，从20世纪70年代起就受到了世界各国的广泛关注。随着风电的迅猛发展，风电场在规模和数目上都出现了爆发性的增长，因此，大规模风电场并网对电力系统的影响也与日俱增。在研究大规模风电接入对电力系统动态特性的影响时，需要一个能准确表征风电场整体动态特性的风电场模型。风电场通常由几十台甚至上百台风电机组组成，其中每一台双馈风机的模型均达到十几阶，采用详细模型对风电场进行建模仿真效率将十分低下，因此有必要对大规模风电场进行简化聚合建模。
目前，大电网中通常将网络等值到较高的电压等级上，而把处于较低电压等级的风电场等值为一个有功功率为负的负荷，这样的处理忽略掉了风电场的动态特性。而随着风电技术的迅猛发展，单机容量已经达到兆瓦级，风力发电容量在电网总发电容量中所占比例越来越大，其动态特性对电网运行特性的影响也日益显著，当进行动态等值时，风电场集中的地区的特性已经不能再简单的忽略掉，因此，急需深入研究风电场的等值化简方法。
分群和参数聚合是风电场动态等值的重要内容。目前，如何对风电场同调特性进行判别的相关文献较少，相关文献提出按照特征量同调分群的方法，但是该类方法存在特征量的选取是否完备的问题，该方法是否存在普适性是阻碍其应用的瓶颈。基于受扰轨线的相似度分群是对多机系统分群常用的方法，常以发电机的功角作为观测量，而与同步电机不同的是，双馈风机中由于电力电子设备的存在，使得其对外的动态特性和传统同步发电机相差较大，因此在对风电机进行同调判别时，对于采用哪些变量的受扰轨线进行比较还没有形成共识。此外，由于双馈风机采用解耦控制，因此对外不体现惯量，其对系统的影响仅仅以发出功率的多少来表征，目前鲜有对于风电机参数聚合的方法研究的相关文献。因此，急需一种适用于风电场的同调判别以及参数聚合的动态等值方法，以满足含有大规模风电场的电力系统的动态等值需求。
技术实现要素：
本发明的目的是针对上述现有技术的不足，提供了一种基于相似度同调的双馈风电机动态等值方法，以实现风电场等值建模过程中双馈风电机同调机群的精确划分以及动态参数的合理聚合，能够大幅降低风电机组的规模并保留风电机组的动态特性。
本发明的目的可以通过如下技术方案实现：
一种基于相似度同调的双馈风电机动态等值方法，所述方法包括以下步骤：
步骤s10、利用prony算法将双馈风电机的输出功率拟合为振幅、相位、频率和阻尼的指数函数的线性组合；
步骤s20、根据步骤s10得到的线性组合，利用相似理论确定双馈风电机间的相似振荡模式；
步骤s30、用相似振荡模式来定义双馈风电机间同调相似度量化指标，判别风电场中的同调机群；
步骤s40、以双馈风电机输出的有功功率为权重得到等值机的等效暂态内电势，并将同群的双馈风电机等效暂态内电势母线分别通过复变比移相变压器连接到等值机的等效暂态内电势母线；
步骤s50、计算各等值机的等效参数，通过等效双馈风电机阻抗将等值机等效暂态内电势母线连接到等值机。
进一步地，步骤s10得到的线性组合公式如下：
其中，
进一步地，步骤s20的具体过程为：双馈风电机输出功率第i个振荡模式的能量ei由该振荡模式下各采样点模值的平方和得到：
ei＝aiexp(jθi)·exp(n·δt·(αi+j2πfi))(2)
其中，n＝0,1,...,j-1，j表示prony算法分解的阶数；
设定一个能量阈值v，能量阈值v的数值应高于99％，具体数值能够根据同调精度要求设定，将所有振荡模式按照能量从高到低排序，并依次累加直到所得能量不小于总能量
利用上述方法消除等值前后振荡模式中的冗余项后，能够获得包含绝大多数信息的等值前l阶振荡模式，l远小于prony算法分解的阶数j，简化了相似振荡模式的寻找过程；
按照特征信息中频率、阻尼、能量三者欧氏距离最小的原则，以其中一台双馈风电机输出功率曲线l阶振荡模式为基准，在其他双馈风电机输出功率曲线m阶振荡模式中寻找相似振荡模式，对于i＝1,2,…,l，k＝1,2,…,m，计算：
其中，fai表示参考双馈风电机输出功率曲线第i个振荡模式的频率，fbk表示其他双馈风电机输出功率曲线第k个振荡模式的频率，αai表示fai的衰减因子，αbk表示fbk的衰减因子，eai表示参考双馈风电机输出功率曲线第i个振荡模式的能量，ebk表示其他双馈风电机输出功率曲线第k个振荡模式的能量，对于振荡模式i，计算获得最小的dik，则参考双馈风电机输出功率曲线振荡模式i与其他双馈风电机输出功率曲线振荡模式k为一对相似振荡模式，在振荡模式中去除振荡模式k，用于下一对相似振荡模式的寻找过程，最终得到l对相似振荡模式。
进一步地，步骤s30的具体过程为：根据步骤s20得到各双馈风电机的l阶相似振荡模式后，第i个相似振荡模式之间的频率相似度q(fi)和阻尼相似度q(αi)分别为：
其中，fai表示参考双馈风电机输出功率曲线第i个振荡模式的频率，fbi表示其他双馈风电机输出功率曲线第i个振荡模式的频率，αai表示fai的衰减因子，αbi表示fbi的衰减因子，考虑频率和阻尼的相似度权重，因此第i个相似振荡模式之间的相似度为：
其中a、b分别表示频率和阻尼的相似度权重；
振荡模式中各个分模式的能量与总能量的比值作为该分模式能量在相似度计算中占全部振荡模式能量的权重：
考虑能量权重后所有相似振荡模式的相似度为：
其中ei表示双馈风电机输出功率第i个振荡模式的能量，计算出双馈风电机间同调相似度量化指标后，确定风电场中的同调机群。
进一步地，步骤s40的具体过程为：首先，将所述各双馈风电机的定子、转子电压方程表示如下：
其中，u、i、ψ分别表示双馈风电机的电压、电流和磁链；r、l分别表示双馈风电机的电阻和电感；下标中的s、r分别表示双馈风电机的定子、转子；下标中的d、q分别表示直轴、交轴分量；lm表示双馈风电机的互感；p为微分算子；s表示转差率；
式中，
其中，xss＞＞rs，所以双馈风电机的定子压降能够忽略不计，定义双馈风电机等效暂态内电势如下：
e'＝us-jxssis(12)
其中，us表示双馈风电机的机端电压，is表示双馈风电机的定子电流，j代表相位超前90度的符号；
由于双馈风电机采用的是解耦控制，因此对外不体现惯量，双馈风电机对系统的影响仅仅通过其注入的功率来表征，采用以输出功率为权重的聚合方法，设同调机群中双馈风电机的个数为n，等值机的等效暂态内电势
其中，pi表示同调机群中第i个双馈风电机的有功功率；等值机的等效暂态内电势由同调机群中所有的双馈风电机决定，根据双馈风电机的外特性，每个双馈风电机所占权重取决于它们的输出功率；
其中，所述复变比移相变压器的变比为：
进一步地，步骤s50中，各等值机的等效参数为：
其中，n表示同调机群中双馈风电机的个数，p表示双馈风电机的有功功率，q表示双馈风电机的无功功率，xs表示双馈风电机的定子电抗，xr表示双馈风电机的转子电抗，hg表示双馈风电机中感性发电机的转子惯性时间常数，ht表示双馈风电机中风力机的转子惯性时间常数，ks表示双馈风电机中风力机的轴系刚度系数，下标中的eq分别代表等值机参数。
本发明与现有技术相比，具有如下优点和有益效果：
1、本发明的一种基于相似度同调的双馈风电机动态等值方法，能够有效地对大规模双馈风电机进行同调机群划分，并对同调机组进行合理参数聚合。该方法能够在保留风电机组动态特性的前提下显著降低风电机组的规模，从而简化外部区域，以便突出重点关注的待研究区域。该方法首先利用prony算法将风电场输出功率拟合为振幅、相位、频率和阻尼的指数函数的线性组合，接着利用相似理论确定相似振荡模式，用相似振荡模式来定义风电机间同调相似度量化指标，判别风电场中的同调机组；然后以双馈风电机输出有功功率为权重得到等值机的等效暂态内电势，并将同群的双馈风电机暂态内电势母线分别通过复变比移相变压器连接到等效暂态内电势母线；最后计算各等值机组的等效参数，通过等效双馈风电机阻抗将等效暂态内电势连接到等值机。这样，满足了含大规模风电机组的大电网高精度等值化简需求，该聚合方法更能反映风电场的真实动态特性，也为科学、正确地分析大规模风电场并网对大电网动态特性的影响奠定了基础。
2、本发明的基于相似度同调的双馈风电机动态等值方法对风电场进行动态等值，能够大幅降低风电机组的规模并保留风电机组的功率特性，更能反映等值后风电场的真实特性，从而满足了含大规模风电机组的大电网高精度等值化简需求。
附图说明
图1为本发明实施例基于相似度同调的双馈风电机动态等值方法的流程图。
图2为本发明实施例构造等值机的步骤图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述，但本发明的实施方式不限于此。
实施例：
本实施例提供了一种基于相似度同调的双馈风电机动态等值方法，所述方法的流程图如图1所示，包括以下步骤：
步骤s10、利用prony算法将双馈风电机的输出功率特性曲线拟合为振幅、相位、频率和阻尼的指数函数的线性组合；线性组合公式如下：
其中，
步骤s20、根据步骤s10得到的线性组合，利用相似理论确定双馈风电机间的相似振荡模式；具体过程为：双馈风电机输出功率第i个振荡模式的能量ei由该振荡模式下各采样点模值的平方和得到：
ei＝aiexp(jθi)·exp(n·δt·(αi+j2πfi))(2)
其中，n＝0,1,...,j-1，j表示prony算法分解的阶数；
设定一个能量阈值v，能量阈值v的数值应高于99％，具体数值能够根据同调精度要求设定，此处将将能量阈值v的数值设为99.9％，所有振荡模式按照能量从高到低排序，并依次累加直到所得能量不小于总能量
利用上述方法消除等值前后振荡模式中的冗余项后，能够获得包含绝大多数信息的等值前l阶振荡模式，l远小于prony算法分解的阶数j，简化了相似振荡模式的寻找过程；
按照特征信息中频率、阻尼、能量三者欧氏距离最小的原则，以其中一台双馈风电机输出功率曲线l阶振荡模式为基准，在其他双馈风电机输出功率曲线m阶振荡模式中寻找相似振荡模式，对于i＝1,2,…,l，k＝1,2,…,m，计算：
其中，fai表示参考双馈风电机输出功率曲线第i个振荡模式的频率，fbk表示其他双馈风电机输出功率曲线第k个振荡模式的频率，αai表示fai的衰减因子，αbk表示fbk的衰减因子，eai表示参考双馈风电机输出功率曲线第i个振荡模式的能量，ebk表示其他双馈风电机输出功率曲线第k个振荡模式的能量，对于振荡模式i，计算获得最小的dik，则参考双馈风电机输出功率曲线振荡模式i与其他双馈风电机输出功率曲线振荡模式k为一对相似振荡模式，在振荡模式中去除振荡模式k，用于下一对相似振荡模式的寻找过程，最终得到l对相似振荡模式。
步骤s30、用相似振荡模式来定义双馈风电机间同调相似度量化指标，判别风电场中的同调机群；具体过程为：根据步骤s20得到各双馈风电机的l阶相似振荡模式后，第i个相似振荡模式之间的频率相似度q(fi)和阻尼相似度q(αi)分别为：
其中，fai表示参考双馈风电机输出功率曲线第i个振荡模式的频率，fbi表示其他双馈风电机输出功率曲线第i个振荡模式的频率，αai表示fai的衰减因子，αbi表示fbi的衰减因子，考虑频率和阻尼的相似度权重，因此第i个相似振荡模式之间的相似度为：
其中a、b分别表示频率和阻尼的相似度权重，此处认为频率和阻尼的相似度权重相等，得到第i个相似振荡模式之间的相似度为：
q(λi)＝0.5q(fi)+0.5q(αi)
振荡模式中各个分模式的能量与总能量的比值作为该分模式能量在相似度计算中占全部振荡模式能量的权重：
考虑能量权重后所有相似振荡模式的相似度为：
其中ei表示双馈风电机输出功率第i个振荡模式的能量，计算出双馈风电机间同调相似度量化指标后，确定风电场中的同调机群。
步骤s40、以双馈风电机输出的有功功率为权重得到等值机的等效暂态内电势，并将同群的双馈风电机等效暂态内电势母线分别通过复变比移相变压器连接到等值机的等效暂态内电势母线；具体过程为：首先，将所述各双馈风电机的定子、转子电压方程表示如下：
其中，u、i、ψ分别表示双馈风电机的电压、电流和磁链；r、l分别表示双馈风电机的电阻和电感；下标中的s、r分别表示双馈风电机的定子、转子；下标中的d、q分别表示直轴、交轴分量；lm表示双馈风电机的互感；p为微分算子；s表示转差率；
式中，
其中，xss＞＞rs，所以双馈风电机的定子压降能够忽略不计，定义双馈风电机等效暂态内电势如下：
e'＝us-jxssis(12)
其中，us表示双馈风电机的机端电压，is表示双馈风电机的定子电流，j代表相位超前90度的符号；
由于双馈风电机采用的是解耦控制，因此对外不体现惯量，双馈风电机对系统的影响仅仅通过其注入的功率来表征，采用以输出功率为权重的聚合方法，设同调机群中双馈风电机的个数为n，等值机的构造流程如图2所示，等值机的等效暂态内电势
其中，pi表示同调机群中第i个双馈风电机的有功功率；等值机的等效暂态内电势由同调机群中所有的双馈风电机决定，根据双馈风电机的外特性，每个双馈风电机所占权重取决于它们的输出功率；
其中，所述复变比移相变压器的变比为：
步骤s50、计算各等值机的等效参数，通过等效双馈风电机阻抗将等值机等效暂态内电势母线连接到等值机，其中各等值机的等效参数为：
其中，n表示同调机群中双馈风电机的个数，p表示双馈风电机的有功功率，q表示双馈风电机的无功功率，xs表示双馈风电机的定子电抗，xr表示双馈风电机的转子电抗，hg表示双馈风电机中感性发电机的转子惯性时间常数，ht表示双馈风电机中风力机的转子惯性时间常数，ks表示双馈风电机中风力机的轴系刚度系数，下标中的eq分别代表等值机参数。
以上所述，仅为本发明专利较佳的实施例，但本发明专利的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明专利所公开的范围内，根据本发明专利的技术方案及其发明专利构思加以等同替换或改变，都属于本发明专利的保护范围。