风电汇集地区动态无功控制交互影响的分析方法技术

本发明专利技术提供一种风电汇集地区动态无功控制交互影响的分析方法，该分析方法包括：根据待分析的风电汇集地区的网架结构形成节点导纳矩阵及网络方程；根据待分析的风电汇集地区的动态无功补偿装置的控制方式建立动态无功补偿装置的动态方程及输出方程；根据网络方程、动态方程及输出方程形成待分析的风电汇集地区的线性化微分方程及线性化状态矩阵；根据线性化状态矩阵获取在动态无功补偿装置的不同控制方式下的特征根；通过对特征根进行参与向量分析生成待分析的评风电汇集地区态无功控制相互作用的分析结果。

【技术实现步骤摘要】
风电汇集地区动态无功控制交互影响的分析方法
本专利技术是关于风电场无功电压控制技术，具体地，是关于一种风电汇集地区动态无功控制交互影响的分析方法。
技术介绍
风能因其可再生、分布广、无污染等特性，已经成为世界可再生能源发展的重要方向。由于我国风能资源主要集中在“三北”地区，而负荷中心多分布在东部地区，决定了目前我国所采取的风电大规模集中接入，长距离输送供给负荷的方式，而这种方式给电力系统带来很多挑战，其中，风电场接入系统的无功电压问题就是比较突出的一个。为应对此问题，《风电接入电网技术规定》中要求并网风电场应配置无功—电压控制系统。目前，在内蒙、冀北等风电汇集地区已经配备自动电压控制系统，风电场自动电压控制(AVC)子站作为二级电压控制模块，将上级调度的电压指令转换为无功指令，下发给动态无功补偿装置的控制器，从而调节并网点母线电压；但在新疆、甘肃等地区的风电场内，动态无功补偿装置则直接接收调度的电压指令来调节并网点母线电压。为避免上述事故的发生，需对静止无功补偿器(SVC)控制方式的风电汇集地区进行电压稳定性评价，现有的电压稳定性评价方法存在的缺点是：该方法仅针对动态无功补偿装置的感性支路恒无功控制方式在系统有电容扰动时汇集地区电压的稳定性，无法评价SVC采用恒电压控制时的电压稳定性，而目前风电场中SVC控制方式多样，除恒无功控制方式，还有高压侧恒电压控制、低压侧恒电压控制等方式。因此，如何能够更加全面、准确地评价风电汇集地区中系统在动态无功补偿装置的不同控制方式下的电压稳定性是当前需要解决的问题。
技术实现思路
本专利技术实施例的主要目的在于提供一种风电汇集地区动态无功控制交互影响的分析方法，以更加全面、准确地评价风电汇集地区中系统在动态无功补偿装置的不同控制方式下的电压稳定性。为了实现上述目的，本专利技术实施例提供一种风电汇集地区动态无功控制交互影响的分析方法，该方法包括：根据待分析的风电汇集地区的网架结构形成节点导纳矩阵及网络方程；根据待分析的风电汇集地区的动态无功补偿装置的控制方式建立所述动态无功补偿装置的动态方程及输出方程；根据所述网络方程、动态方程及输出方程形成待分析的风电汇集地区的线性化微分方程及线性化状态矩阵；根据所述线性化状态矩阵获取在动态无功补偿装置的不同控制方式下的特征根；通过对所述特征根进行参与向量分析生成待分析的评风电汇集地区态无功控制相互作用的分析结果。在一实施例中，上述的网络方程为：I＝YU，其中，I为系统中各母线的注入电流矩阵、U为系统中各母线的节点电压矩阵，Y为系统的所述节点导纳矩阵。在一实施例中，上述的动态方程为：其中，Bs1i为第一状态变量；T1i为采样延时时间；Ui为从系统中得到的控制点的实际电压；Bs2i第二状态变量；KIi为控制器积分系数；KPi为控制器比例系数；Urefi为SVC的参考电压；BSVci为第三状态变量，即SVC输出的导纳；T2i为SVC触发时间常数。在一实施例中，上述的输出方程为：Ii＝BiUi，其中，Ui为控制点的节点电压、Ii为控制点的注入电流，Bi为动态无功补偿装置在控制点的等效导纳。在一实施例中，上述的根据所述网络方程、动态方程及输出方程形成待分析的风电汇集地区的线性化微分方程及线性化状态矩阵，包括：根据二阶同步发电机的控制模型、所述动态方程及输出方程生成状态方程；对所述网络方程及状态方程分别进行线性化，并将线性化后的状态方程代入线性化后的网络方程，生成所述的线性化微分方程及线性化状态矩阵。在一实施例中，上述的二阶同步发电机的控制模型为：及其中，δ为功率角；ω为发电机转子角频率；ω0为发电机转子额定角频率；TJ为惯性时间常数；Pm为输入功率；Pe为电磁功率；D为阻尼系数；Ud为定子d轴电压；E'd为d轴暂态电势；x'q为q轴暂态电抗；Iq为定子q轴电流；r为定子电阻；Id为定子d轴电流；Uq为定子q轴电压；E'q为q轴暂态电势；x'd为d轴暂态电抗。在一实施例中，上述的状态方程为：其中，为发电机状态变量的微分矩阵；XG为发电机状态变量矩阵；UG为发电机的输出电压矩阵；IG为发电机的输出；为SVC状态变量的微分矩阵；XSVCi为SVC状态变量矩阵；USVCi为SVC的输出电压矩阵；ISVCi为SVC的输出电流矩阵；AG、BG、CG、DG、ASVCi、BSVCi、CSVCi、DSVCi为各矩阵前面的系数矩阵，无确定的物理含义。在一实施例中，上述的线性化微分方程为：其中，△X＝[△XG…△XSVCi…]，A＝diag{AG…ASVCi...}，B＝diag{BG…BSVCi…}，C＝diag{-CG…-CSVCi…}，D＝Y-diag{DG…DSVCi…}。在一实施例中，上述的线性化状态矩阵为：其中，A＝diag{AG…ASVCi...}；B＝diag{BG…BSVCi…}；D＝Y-diag{DG…DSVCi…}；C＝diag{-CG…-CSVCi…}。在一实施例中，上述的根据所述线性化状态矩阵获取在动态无功补偿装置的不同控制方式下的特征根，包括：根据所述线性化状态矩阵分别获取在动态无功补偿装置的恒无功控制方式下的特征根λg(A)＝{λ1,...,λn}及恒电压控制方式下的特征根λv(A)＝{λ1,...,λn}。本专利技术实施例的有益效果在于，本专利技术实施例的风电汇集地区动态无功控制交互影响的分析方法，从动态无功补偿装置的控制方式入手研究风电汇集地区电压振荡的原因，通过小干扰稳定分析法求解系统的特征根，并分析系统在动态无功补偿装置的不同控制方式下特征根的变化，可以分析系统在动态无功补偿装置的不同控制方式下的稳定性，对汇集地区动态无功控制方式的选择起到指导意义。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1为根据本专利技术实施例的风电汇集地区动态无功控制交互影响的分析方法的流程图；图2为典型风电汇集地区拓扑结构图；图3为SVC恒电压控制框图；图4为SVC采用不同无功控制方式时的根轨迹示意图；图5为SVC高压侧恒电压控制下的时域仿真图；图6为SVC低压侧恒电压控制下的时域仿真图；图7为SVC恒无功控制下的时域仿真图。具体实施方式下面将结合本专利技术实施例中的附图，对本专利技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本专利技术一部分实施例，而不是全部的实施例。基于本专利技术中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本专利技术保护的范围。本专利技术实施例提供一种风电汇集地区动态无功控制交互影响的分析方法。以下结合附图对本专利技术进行详细说明。本专利技术实施例提供一种风电汇集地区动态无功控制交互影响的分析方法，如图1所示，该分析方法主要包括以下各步骤：步骤S101：根据待分析的风电汇集地区的网架结构形成节点导纳矩阵及网络方程；步骤S102：根据待分析的风电汇集地区的动态无功补偿装置的控制方式建立动态无功补偿装置的动态方程及输出方程；步骤S103：根据网络方程、动态方程及输出方程形成待分析的风电汇集地区的线性化微分本文档来自技高网...

【技术保护点】
一种风电汇集地区动态无功控制交互影响的分析方法，其特征在于，所述的分析方法包括：根据待分析的风电汇集地区的网架结构形成节点导纳矩阵及网络方程；根据待分析的风电汇集地区的动态无功补偿装置的控制方式建立所述动态无功补偿装置的动态方程及输出方程；根据所述网络方程、动态方程及输出方程形成待分析的风电汇集地区的线性化微分方程及线性化状态矩阵；根据所述线性化状态矩阵获取在动态无功补偿装置的不同控制方式下的特征根；通过对所述特征根进行参与向量分析生成待分析的评风电汇集地区态无功控制相互作用的分析结果。

【技术特征摘要】
1.一种风电汇集地区动态无功控制交互影响的分析方法，其特征在于，所述的分析方法包括：根据待分析的风电汇集地区的网架结构形成节点导纳矩阵及网络方程；根据待分析的风电汇集地区的动态无功补偿装置的控制方式建立所述动态无功补偿装置的动态方程及输出方程；所述动态方程为：其中，Bs1i为第一状态变量，T1i为采样延时时间，Ui为从系统中得到的控制点的实际电压，Bs2i第二状态变量，KIi为控制器积分系数，KPi为控制器比例系数，Urefi为SVC的参考电压，BSVci为第三状态变量，即SVC输出的导纳，T2i为SVC触发时间常数；所述输出方程为：Ii＝BiUi，其中，Ui为控制点的节点电压，Ii为控制点的注入电流，Bi为动态无功补偿装置在控制点的等效导纳；根据所述网络方程、动态方程及输出方程形成待分析的风电汇集地区的线性化微分方程及线性化状态矩阵，其中，包括根据二阶同步发电机的控制模型、所述动态方程及输出方程生成状态方程；所述二阶同步发电机的控制模型为：及其中，δ为功率角，ω为发电机转子角频率，ω0为发电机转子额定角频率，TJ为惯性时间常数，Pm为输入功率，Pe为电磁功率，D为阻尼系数，Ud为定子d轴电压，E'd为d轴暂态电势，x'q为q轴暂态电抗，Iq为定子q轴电流，r为定子电阻，Id为定子d轴电流，Uq为定子q轴电压，E'q为q轴暂态电势，x'd为d轴暂态电抗；所述状态方程为：其中，为发电机状态变...

【专利技术属性】
技术研发人员：谢欢，周红婷，吴涛，曹天植，李善颖，赵峰，梁浩，姚谦，付宏伟，
申请(专利权)人：华北电力科学研究院有限责任公司，国家电网公司，国网冀北电力有限公司电力科学研究院，
类型：发明
国别省市：北京;11