[0001] 本发明涉及混合储能技术领域，具体涉及一种混合储能VSG控制系统稳定性分析方法及系统。

[0002] 作为我国能源变革关键技术支撑之一的大规模储能技术，因为其可以为电网提供调峰、调频、应急响应等多种辅助服务，近年来受到了业内的广泛关注。分布式混合储能系统由蓄电池与超级电容组成，通过换流器与电网相连，具备功率双向流动能力，具有功率调节迅速、应用模式多样等优点。分布式混合储能系统因为其低惯性、低阻尼特性会对网侧的频率稳定性造成影响。为了使储能换流器控制系统具备如同步电机拥有的转动惯量和阻尼，目前广泛采用虚拟同步机(Virtual Synchronous Generation,VSG)控制技术，使其像同步发电机一样参与到电网频率和电压的调节中。然而目前混合储能VSG控制系统参与调节效果不理想，依然对网侧的频率稳定性造成影响。

[0044] 下面将结合附图对本发明的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0045] 在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。

[0046] 在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，还可以是两个元件内部的连通，可以是无线连接，也可以是有线连接。对于本领域的普通技术人员而言，可以具体情况理解上述术语在本发明中的具体含义。

[0047] 此外，下面所描述的本发明不同实施方式中所涉及的技术特征只要彼此之间未构成冲突就可以相互结合。

[0048] 本发明实施例提供一种混合储能VSG控制系统稳定性分析方法，如图1所示，包括如下步骤：

[0049] 步骤S1：建立混合储能虚拟同步机转子机械方程和无功调节方程。

[0050] 在一具体实施例中，如图2所示，分布式电源点通过DC/DC换流器与直流母线相连，因为其有功出力PDG波动性，需要配置由蓄电池和超级电容组成的混合储能设备，通过调节蓄电池有功功率PB、超级电容有功功率PC来平抑PDG波动。当交流大电网发生频率振荡事故，功率振荡信号通过AC/DC换流器传递至直流侧，需要混合储能具备频率振荡抑制能力。

[0051] 如图3所示，Udc为直流母线电压；Uabc、iabc为混合储能DC/AC换流器交流侧三相电压、电流；eabc为交流电网三相电压；Rf、Lf、Lg、Rg构成滤波电路。VSG控制流程中，ω为虚拟角频率；δ为虚拟功角；Pref为有功功率参考值；Pe为虚拟电磁功率，VSG无定子绕组损耗，Pe又称VSG控制输出功率；Pm为虚拟机械功率；u为VSG中虚拟内电势；Qe为无功功率实际输出值；Qref为无功功率参考值。转子机械方程为：

[0052]

[0053] 式(1)中，H为虚拟转动惯量；Tm、Te分别为VSG机械转矩、电磁转矩，D为虚拟阻尼系数；ω0为额定角频率；在VSG控制中，虚拟转动惯量H使得
混合储能换流器在功率和频率调节过程中具有了惯性，阻尼系数D使得混合储能DC/AC换流器具有抑制交流电网功率振荡的能力。DC/AC换流器控制同样具有励磁调节惯性，无功调节方程为：

[0054]

[0055] 式(2)中，U0为额定电压有效值；Δu为虚拟内电势与额定电压之间偏差；kq为无功调节系数。

[0056] 在本发明实施例中，通过采用VSG控制策略，可以为混合储能换流器增加类似于同比发电机的惯性和阻尼效果。

[0057] 步骤S2：将混合储能虚拟同步机转子机械方程和无功调节方程中的状态变量进行线性化及拉普拉斯变化，并用小扰动信号进行表示。

[0058] 在一具体实施例中，为了验证所提自适应虚拟阻尼VSG控制系统稳定性，本实施例对混合储能VSG控制进行工频小信号分析，将式(1)VSG转子机械方程及式(2)VSG无功调节方程中的状态变量进行线性化及拉普拉斯变化，用小扰动信号进行表示可以得到：

[0059]

[0060] 式(3)中：s为拉普拉斯算子； 为虚拟角频率、虚拟功角、虚拟机械功率、VSG输出功率对应的在直流工作点小扰动量；Se为VSG输出视在功率，因为δ很小，则有sinδ≈δ，可以得到：

[0061] 步骤S3：根据混合储能虚拟同步机转子机械方程和无功调节方程的小扰动信号模型计算得到混合储能VSG控制系统输出有功功率。

[0062] 在一具体实施例中，通过如下步骤计算得到混合储能VSG控制系统输出有功功率：

[0063] 步骤S31：根据混合储能虚拟同步机转子机械方程和无功调节方程的小扰动信号模型计算得到传递函数方程。

[0064] 步骤S32：根据传递函数方程计算得到混合储能VSG控制系统输出有功功率。

[0065] 传递函数用G(s)表示：

[0066]

[0067] 式(4)中， Pe(s)、Pm(s)为VSG输出功率、虚拟机械功率在s域的工频小信号模型。当混合储能发生频率振荡时，VSG控制系统输出有功表示为：

[0068]

[0069] 式(5)中：ΔPm(s)为虚拟机械功率变化量在s域的工频小信号模型。

[0070] 步骤S4：将混合储能VSG控制系统输出有功功率最大值对虚拟阻尼系数进行求导，进行混合储能VSG控制系统稳定性分析。

[0071] 在一具体实施例中，通过如下步骤计算将混合储能VSG控制系统输出有功功率最大值对虚拟阻尼系数进行求导，进行混合储能VSG控制系统稳定性分析：

[0072] 步骤S41：根据混合储能VSG控制系统输出有功功率计算得到VSG控制系统输出有功功率时域内响应函数。

[0073] 步骤S42：根据VSG控制系统输出有功功率时域内响应函数得到有功功率最大值。

[0074] 步骤S43：将混合储能VSG控制系统输出有功功率最大值对虚拟阻尼系数进行求导，分析虚拟阻尼系数对输出功率的影响。

[0075] 具体地，根据式(5)，可以得到Pe时域内响应函数为：

[0076]

[0077] 式(6)中： 根据式(6)，Pe最大值为：

[0078]

[0079] 为了分析虚拟阻尼系数D对输出功率Pe的影响，将式(7)中Pe.max对D进行求导，可得：

[0080]

[0081] 根据式(8)可知：当混合储能VSG控制系统输出功率增加时，虚拟机械功率也相应增加，即ΔPm＞0，Pe.max对D导数＜0；当混合储能VSG控制系统输出功率减少时，虚拟机械功率也相应减少，即ΔPm＜0，Pe.max对D导数＞0。通过自适应调节D可以限制功率正向、反向超调，进一步的抑制混合储能频率振荡。

[0082] 本发明提供的混合储能VSG控制系统稳定性分析方法，包括：建立混合储能虚拟同步机转子机械方程和无功调节方程；将混合储能虚拟同步机转子机械方程和无功调节方程中的状态变量进行线性化及拉普拉斯变化，并用小扰动信号进行表示；根据混合储能虚拟同步机转子机械方程和无功调节方程的小扰动信号模型计算得到混合储能VSG控制系统输出有功功率；将混合储能VSG控制系统输出有功功率最大值对虚拟阻尼系数进行求导，进行混合储能VSG控制系统稳定性分析。通过对VSG控制系统线性化及拉普拉斯变化，分析混合储能VSG控制系统稳定性，并得到虚拟阻尼与储能频率振荡的关系。从而根据虚拟阻尼与储能频率振荡的关系制抑制混合储能频率振荡，到达混合储能VSG控制系统参与调节的理想效果。

[0083] 本发明实施例提供一种混合储能VSG控制系统稳定性分析系统，如图4所示，包括：

[0084] 模型建立模块1，用于建立混合储能虚拟同步机转子机械方程和无功调节方程。详细内容参见上述实施例中步骤S1的相关描述，在此不再赘述。

[0085] 第一计算模块2，用于将混合储能虚拟同步机转子机械方程和无功调节方程中的状态变量进行线性化及拉普拉斯变化，并用小扰动信号进行表示。详细内容参见上述实施例中步骤S2的相关描述，在此不再赘述。

[0086] 第二计算模块3，用于根据混合储能虚拟同步机转子机械方程和无功调节方程的小扰动信号模型计算得到混合储能VSG控制系统输出有功功率。详细内容参见上述实施例中步骤S3的相关描述，在此不再赘述。

[0087] 第三计算模块4，用于将混合储能VSG控制系统输出有功功率最大值对虚拟阻尼系数进行求导，进行混合储能VSG控制系统稳定性分析。详细内容参见上述实施例中步骤S4的相关描述，在此不再赘述。

[0088] 本发明实施例还提供一种计算机设备，如图5所示，该设备可以包括处理器61和存储器62，其中处理器61和存储器62可以通过总线或者其他方式连接，图5以通过总线连接为例。

[0089] 处理器61可以为中央处理器(Central Processing Unit，CPU)。处理器61还可以为其他通用处理器、数字信号处理器(Digital Signal Processor，DSP)、专用集成电路(Application Specific Integrated Circuit，ASIC)、现场可编程门阵列(Field‑Programmable Gate Array，FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件等芯片，或者上述各类芯片的组合。

[0090] 存储器62作为一种非暂态计算机可读存储介质，可用于存储非暂态软件程序、非暂态计算机可执行程序以及模块，如本发明实施例中的对应的程序指令/模块。处理器61通过运行存储在存储器62中的非暂态软件程序、指令以及模块，从而执行处理器的各种功能应用以及数据处理，即实现上述方法实施例中的混合储能VSG控制系统稳定性分析方法。

[0091] 存储器62可以包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需要的应用程序；存储数据区可存储处理器61所创建的数据等。此外，存储器62可以包括高速随机存取存储器，还可以包括非暂态存储器，例如至少一个磁盘存储器件、闪存器件、或其他非暂态固态存储器件。在一些实施例中，存储器62可选包括相对于处理器61远程设置的存储器，这些远程存储器可以通过网络连接至处理器61。上述网络的实例包括但不限于互联网、企业内部网、企业内网、移动通信网及其组合。

[0092] 一个或者多个模块存储在存储器62中，当被处理器61执行时，执行本发明实施提供的混合储能VSG控制系统稳定性分析方法。

[0093] 上述计算机设备具体细节可以对应参阅图1‑图3所示的实施例中对应的相关描述和效果进行理解，此处不再赘述。

[0094] 本领域技术人员可以理解，实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成的程序可存储于一计算机可读取存储介质中，该程序在执行时，可包括如上述各方法的实施例的流程。其中，存储介质可为磁碟、光盘、只读存储记忆体(Read‑Only Memory，ROM)、随机存储记忆体(Random Access Memory，RAM)、快闪存储器(Flash Memory)、硬盘(Hard Disk Drive，缩写：HDD)或固态硬盘(Solid‑State Drive，SSD)等；存储介质还可以包括上述种类的存储器的组合。

[0095] 显然，上述实施例仅仅是为清楚地说明所作的举例，而并非对实施方式的限定。对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。而由此所引申出的显而易见的变化或变动仍处于本发明创造的保护范围之中。