CN112467790A - 一种抑制mmc互联变换器虚拟同步机功率震荡的方法 - Google Patents

一种抑制mmc互联变换器虚拟同步机功率震荡的方法 Download PDF

Info

Publication number
CN112467790A
CN112467790A CN202011299319.0A CN202011299319A CN112467790A CN 112467790 A CN112467790 A CN 112467790A CN 202011299319 A CN202011299319 A CN 202011299319A CN 112467790 A CN112467790 A CN 112467790A
Authority
CN
China
Prior art keywords
power
virtual
mmc
grid
synchronous machine
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011299319.0A
Other languages
English (en)
Other versions
CN112467790B (zh
Inventor
冯仰敏
杨沛豪
赵勇
李立勋
张恩享
王燕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian Thermal Power Research Institute Co Ltd
Original Assignee
Xian Thermal Power Research Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian Thermal Power Research Institute Co Ltd filed Critical Xian Thermal Power Research Institute Co Ltd
Priority to CN202011299319.0A priority Critical patent/CN112467790B/zh
Publication of CN112467790A publication Critical patent/CN112467790A/zh
Application granted granted Critical
Publication of CN112467790B publication Critical patent/CN112467790B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381Dispersed generators
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24Arrangements for preventing or reducing oscillations of power in networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/36Arrangements for transfer of electric power between ac networks via a high-tension dc link
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J5/00Circuit arrangements for transfer of electric power between ac networks and dc networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
    • Y02E60/60Arrangements for transfer of electric power between AC networks or generators via a high voltage DC link [HVCD]

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

本发明公开了一种抑制MMC互联变换器虚拟同步机功率震荡的方法,主要包括:以交直流混合微电网MMC互联变换器交流侧角度出发,得到电压方程;建立MMC互联变换器功率输出方程;建立MMC互联变换器功率调节关系式;得到基于功率传输原则的虚拟同步机控制机械方程和虚拟同步机无功励磁控制方程;对MMC互联变换器功率输出化简方程中的有功功率求一阶、二阶导数;得到虚拟同步机参考输出功率表达式;虚拟同步机参考输出功率表达式建立小信号模型并进行拉氏变换,计算其特征根;得到系统自然震荡转速和阻尼比设计自适应虚拟惯量和自适应虚拟阻尼系数,带入基于功率传输原则的虚拟同步机控制机械方程中,起到抑制MMC互联变换器虚拟同步机功率震荡的作用。

Description

一种抑制MMC互联变换器虚拟同步机功率震荡的方法
技术领域
本发明涉及一种抑制MMC互联变换器虚拟同步机功率震荡的方法,该方法一种抑制混合微电网MMC互联变换器虚拟同步机功率震荡的方法,采用自适应虚拟参数代替传统虚拟参数,增加控制系统对的功率震荡的抑制能力。
背景技术
交直流混合微电网中MMC互联变换器接口变换器通常采用下垂控制,常规下垂控制不具备惯性和阻尼环节,无法维持分布式电源高渗透率,虚拟同步机(VirtualSynchronous Generation,VSG)控制通过模仿同步发电机特性能够为系统提供惯性和阻尼支撑,实现一次调频和一次调压。虚拟同步机暂态变换过程中缓冲能量变化量采用虚拟转动惯量J变化表征,各种摩擦阻力采用虚拟阻尼系数D变化量表征。
但常规虚拟转动惯量J为定值,J选值过小,系统响应时间可以减少但无法对功率震荡起到抑制作用;J选值过大,系统可以对功率波动起到抑制作用但增加了响应时间。同理,虚拟阻尼系数D选取也会影响功率震荡抑制性能。
发明内容
本发明的目的在于提供一种抑制MMC互联变换器虚拟同步机功率震荡的方法,具体应用自适应虚拟惯量和自适应阻尼系数来抑制混合微电网MMC互联变换器虚拟同步机功率调节过程中的功率震荡问题,通过小信号模型分析额定虚拟参数取值范围。
本发明采取如下技术方案来实现的:
一种抑制MMC互联变换器虚拟同步机功率震荡的方法,包括以下步骤:
1)以交直流混合微电网MMC互联变换器交流侧角度出发,得到电压方程;
3)根据步骤1)MMC互联变换器交流侧电压方程,建立MMC互联变换器功率输出方程,并进行化简;
3)根据交直流混合微电网MMC互联变换器控制系统中,交流微电网、直流微电网瞬时有功功率变化量相同,建立MMC互联变换器功率调节关系式;
4)将步骤3)MMC互联变换器功率调节关系式中交流微电网有功出力分解为稳态功率变化量、瞬时功率变化量,直流微电网下垂控制模型有功出力调节量分解为稳态功率变化量、动态功率变化量;
5)根据步骤3)MMC互联变换器功率调节关系式和步骤4)交流微电网、直流微电网有功出力分解式得到基于功率传输原则的虚拟同步机控制机械方程和虚拟同步机无功励磁控制方程;
6)根据发生功率震荡时,角频率变化情况,将步骤5)虚拟同步机控制机械方程中的虚拟惯量和虚拟阻尼修改为自适应虚拟惯量和自适应虚拟阻尼系数;
7)对步骤2)MMC互联变换器功率输出化简方程中的有功功率求一阶、二阶导数;
8)将步骤7)中MMC互联变换器功输出有功功率一阶、二阶导数和步骤4)交流微电网有功出力分解表达式带入步骤5)基于功率传输原则的虚拟同步机控制机械方程中,得到虚拟同步机参考输出功率表达式;
9)对步骤8)中虚拟同步机参考输出功率表达式建立小信号模型并进行拉氏变换,计算其特征根;
10)分析步骤9)虚拟同步机参考输出功率小信号模型和其特征根,得到系统自然震荡转速和阻尼比;
11)参考同步发电机振荡频率,得到额定虚拟惯量取值范围,在考虑阻尼系数的基础上,利用最优二阶系统分析法,考虑阻尼系数的基础上,设置步骤10)阻尼比和额定阻尼系数;
12)根据步骤11)额定虚拟惯量取值范围和阻尼比和额定阻尼系数,设计步骤6)自适应虚拟惯量和自适应虚拟阻尼系数,带入步骤5)基于功率传输原则的虚拟同步机控制机械方程中,起到抑制MMC互联变换器虚拟同步机功率震荡的作用。
本发明进一步的改进在于,步骤1)以交直流混合微电网MMC互联变换器交流侧角度出发,得到电压方程,
Figure BDA0002786337740000031
其中:
Figure BDA0002786337740000032
为交流电网电压e0对应的矢量;
Figure BDA0002786337740000033
为MMC互联变换器交流侧电压Uac对应的矢量;
Figure BDA0002786337740000034
为交流电网电流i0对应的矢量。
本发明进一步的改进在于,步骤2)的具体实现方法为:根据步骤1)MMC互联变换器交流侧电压方程,建立MMC互联变换器功率传输方程:
Figure BDA0002786337740000035
其中:Rf、Xf为滤波电路电阻值及感抗;δ为交流电网电压矢量
Figure BDA0002786337740000036
与MMC互联变换器交流侧电压矢量
Figure BDA0002786337740000037
之间的相角差;交流电网电压e0类比于同步电机电枢电动势;MMC互联变换器交流侧电压Uac类比于同步电机端电压;在交直流混合微电网MMC互联变换器控制系统中,通过控制δ正负,实现功率双向流动,当δ>0,
Figure BDA0002786337740000038
超前
Figure BDA0002786337740000039
MMC互联变换器运行于逆变模式,功率由直流微电网向交流微电网传输;当δ<0,
Figure BDA00027863377400000310
滞后
Figure BDA00027863377400000311
MMC互联变换器运行于整流模式,功率由交流微电网向直流微电网传输;当δ=0,
Figure BDA00027863377400000312
Figure BDA00027863377400000313
同相位,交流微电网、直流微电网之间无功率交换;MMC互联变换器,满足Rf<<Xf,对MMC互联变换器功率传输方程进行化简:
Figure BDA00027863377400000314
交流电网电压矢量与MMC互联变换器交流侧电压矢量
Figure BDA00027863377400000315
之间的相角差δ很小,则有:
Figure BDA00027863377400000316
有功功率表示为:
Figure BDA00027863377400000317
本发明进一步的改进在于,步骤3)的具体实现方法为:根据交直流混合微电网MMC互联变换器控制系统中,交流微电网、直流微电网瞬时有功功率变化量相同,建立MMC互联变换器功率调节关系式:Pacref-Pac=Pdc-Pdcref=ΔP;
其中:Pac为交流微电网输出功率实际值;Pacef为交流微电网输出功率实际值;Pdc为直流微电网输出功率实际值;Pdcref为直流微电网输出功率实际值;ΔP为MMC互联变换器功率调节量。
本发明进一步的改进在于,步骤4)的具体实现方法为:将步骤3)MMC互联变换器功率调节关系式中交流微电网有功出力分解为稳态功率变化量、瞬时功率变化量:
Figure BDA0002786337740000041
其中:kω为交流电网下垂调节系数;kω(ω-ω0)是稳态有功变化量;
Figure BDA0002786337740000042
为虚拟转子惯量吸收或者发出的瞬时有功功率,在虚拟同步机控制中为交流频率提供惯性环节;直流微电网下垂控制模型有功出力调节量分解为稳态功率变化量、动态功率变化量:
Figure BDA0002786337740000043
其中:kudc为直流电网下垂调节系数;Udc为直流母线电压实际值;Udc0为直流母线电压初始值;Cdc为直流侧电容值;kudc(Udc-Udc0)是稳态有功变化量,
Figure BDA0002786337740000044
是直流电容的充放电功率,属于动态功率波动。
本发明进一步的改进在于,步骤5)的具体实现方法为:根据步骤3)MMC互联变换器功率调节关系式和步骤4)交流微电网、直流微电网有功出力分解式得到基于功率传输原则的虚拟同步机控制机械方程和虚拟同步机无功励磁控制方程
Figure BDA0002786337740000045
其中:J为虚拟转动惯量;为交流微电网侧角频率当前值;为交流微电网侧角频率初始值;D为虚拟阻尼系数;kudc为直流微电网下垂调节系数;Udc为直流侧母线电压当前值;Udc0为直流侧母线电压初始值;Cdc为直流侧电容值;拟同步机感应内电势由两部分组成:一部分是空载时,虚拟励磁电压对于的空载电动势,一部分是由无功功率偏差产生,虚拟同步机无功励磁控制方程为:E=E0+kq(Qref-Q);
其中:E为虚拟同步机感应内电势有效值;E0为励磁空载电动势有效值;kq为无功电压下垂控制系数;Qref为无功功率参考值;Q为无功功率当前值;综合虚拟同步机有功频率控制计算得到的虚拟转子角频率和相角差δ,得到交直流混合微电网MMC互联变换器等效交流输出端口三相电压调制信号为:
Figure BDA0002786337740000051
本发明进一步的改进在于,步骤6)的具体实现方法为:根据发生功率震荡时,角频率变化情况,在角频率增加阶段,增加J同时适当减少D;在角频率减少阶段,减少J同时适当增加D的自适应控制方法;该方法保证系统响应速度,并且加快功率进入稳定状态将步骤5)虚拟同步机控制机械方程中的虚拟惯量和虚拟阻尼修改为自适应虚拟惯量和自适应虚拟阻尼系数:
Figure BDA0002786337740000052
其中:kj为虚拟惯量调节系数;J0为额定虚拟转动惯量;ω0为额定机械角频率;kd为虚拟阻尼调节系数;D0为额定虚拟阻尼系数;在自适应虚拟同步机控制系统中,虚拟惯量调节系数kj表征虚拟转动惯量跟随频率偏差调节的能力,kd表征虚拟阻尼系数跟随频率偏差的能力。
本发明进一步的改进在于,步骤7)的具体实现方法为:对步骤2)MMC互联变换器功率输出化简方程中的有功功率求一阶、二阶导数:
Figure BDA0002786337740000061
本发明进一步的改进在于,步骤8)的具体实现方法为:将步骤7)中MMC互联变换器功输出有功功率一阶、二阶导数和步骤4)交流微电网有功出力分解表达式带入步骤5)基于功率传输原则的虚拟同步机控制机械方程中,得到虚拟同步机参考输出功率表达式:
Figure BDA0002786337740000062
本发明进一步的改进在于,步骤9)的具体实现方法为:对步骤8)中虚拟同步机参考输出功率表达式建立小信号模型并进行拉氏变换:
Figure BDA0002786337740000063
计算其特征根:
Figure BDA0002786337740000064
步骤10)的具体实现方法为:分析步骤9)虚拟同步机参考输出功率小信号模型和其特征根,得到系统自然震荡转速和阻尼比:
Figure BDA0002786337740000065
步骤11)的具体实现方法为:参考同步发电机振荡频率:0.628rad/s≤ωS≤15.7rad/s,得到额定虚拟惯量取值范围:
Figure BDA0002786337740000066
在考虑阻尼系数的基础上,利用最优二阶系统分析法,考虑阻尼系数的基础上,设置步骤10)阻尼比和额定阻尼系数:
Figure BDA0002786337740000067
步骤12)的具体实现方法为:根据步骤11)额定虚拟惯量取值范围:
Figure BDA0002786337740000071
和阻尼比和额定阻尼系数:
Figure BDA0002786337740000072
设计步骤6)自适应虚拟惯量J和自适应虚拟阻尼系数D,带入步骤5)基于功率传输原则的虚拟同步机控制机械方程中,起到抑制MMC互联变换器虚拟同步机功率震荡的作用。
与现有技术相比,本发明至少具有如下有益的技术效果:
1.本发明提出一种在角频率增加阶段,增加J同时适当减少D;在角频率减少阶段,减少J同时适当增加D的自适应控制方案,该方案可以保证虚拟同步机响应速度,并且加快功率进入稳定状态。可以有效抑制混合微电网MMC互联变换器虚拟同步机功率震荡。
2.本发明采用小信号模型分析得到的额定虚拟惯量、额定阻尼系数取值更为合理。
附图说明
图1为MMC互联变换器主电路拓扑图;
图2为基于功率传输原则的MMC互联变换器虚拟同步机控制框图;
图3为虚拟同步机功率与角频率变化曲线,图3(a)为发生功率震荡时,图3(b)为角频率减少阶段;
图4为自适应虚拟同步机控制框图;
图5为MMC互联变换器自适应虚拟同步机控制模型;
图6为在0.5s时刻交流微电网负荷增加3kW,1s时刻3kW负荷切除,有功功率动态仿真波形;
图7为在0.5s时刻交流微电网负荷增加3kW,1s时刻3kW负荷切除,频率动态仿真波形;
图8为在0.5s时刻交流微电网负荷增加3kW,1s时刻3kW负荷切除,自适应虚拟惯量变化曲线仿真波形;
图9为在0.5s时刻交流微电网负荷增加3kW,1s时刻3kW负荷切除,自适应阻尼系数变化曲线仿真波形
图10为在0.5s时刻直流微电网负荷增加3kW,1s时刻3kW负荷切除,有功功率动态仿真波形;
图11为在0.5s时刻直流微电网负荷增加3kW,1s时刻3kW负荷切除,频率动态仿真波形;
图12为在0.5s时刻直流微电网负荷增加3kW,1s时刻3kW负荷切除,自适应虚拟惯量变化曲线仿真波形;
图13为在0.5s时刻直流微电网负荷增加3kW,1s时刻3kW负荷切除,自适应阻尼系数变化曲线仿真波形。
具体实施方式
下面通过附图,对本发明的技术方案做进一步的详细描述。
如图1所示,以交直流混合微电网MMC互联变换器交流侧角度出发,交流侧电压方程可表示为:
Figure BDA0002786337740000081
式(1)中:
Figure BDA0002786337740000082
为交流电网电压e0对应的矢量;
Figure BDA0002786337740000083
为MMC互联变换器交流侧电压Uac对应的矢量;
Figure BDA0002786337740000084
为交流电网电流i0对应的矢量。MMC互联变换器功率传输方程可表示为:
Figure BDA0002786337740000085
式(2)中:Rf、Xf为滤波电路电阻值及感抗;δ为交流电网电压矢量
Figure BDA0002786337740000091
与MMC互联变换器交流侧电压矢量
Figure BDA0002786337740000092
之间的相角差。MMC互联变换器,通常满足Rf<<Xf,则式(2)可简化为:
Figure BDA0002786337740000093
式(1)、式(3)所描述MMC互联变换器交流侧电压方程、功率传输方程分别类比于同步电机的电压方程、功率方程。交流电网电压e0类比于同步电机电枢电动势;MMC互联变换器交流侧电压Uac类比于同步电机端电压。同步电机根据电枢电动势和端电压相角差的正负,即,可以作为电动机或者发电机。同样,在交直流混合微电网MMC互联变换器控制系统中,通过控制δ正负,实现功率双向流动。当δ>0,
Figure BDA0002786337740000094
超前
Figure BDA0002786337740000095
MMC互联变换器运行于逆变模式,功率由直流微电网向交流微电网传输;当δ<0,
Figure BDA0002786337740000096
滞后
Figure BDA0002786337740000097
MMC互联变换器运行于整流模式,功率由交流微电网向直流微电网传输;当δ=0,
Figure BDA0002786337740000098
Figure BDA0002786337740000099
同相位,交流微电网、直流微电网之间无功率交换。
如图2所示,在交直流混合微电网MMC互联变换器控制系统中,交流微电网、直流微电网瞬时有功功率变化量相同,存在:
Pacref-Pac=Pdc-Pdcref=ΔP (4)
式(4)中,Pac为交流微电网输出功率实际值;Pacef为交流微电网输出功率实际值;Pdc为直流微电网输出功率实际值;Pdcref为直流微电网输出功率实际值;ΔP为MMC互联变换器功率调节量。
交流微电网中有功出力调节量可表示为:
Figure BDA00027863377400000910
式(5)中:kω为交流电网下垂调节系数;kω(ω-ω0)是稳态有功变化量。
Figure BDA00027863377400000911
为虚拟转子惯量吸收或者发出的瞬时有功功率,在虚拟同步机控制中为交流频率提供惯性环节。
直流微电网中有功出力调节量可表示为:
Figure BDA0002786337740000101
式(6)中:kudc为直流电网下垂调节系数;Udc为直流母线电压实际值;Udc0为直流母线电压初始值;Cdc为直流侧电容值。kudc(Udc-Udc0)是稳态有功变化量,
Figure BDA0002786337740000102
是直流电容的充放电功率,属于动态功率波动。
基于功率传输原则的虚拟同步机控制机械方程为:
Figure BDA0002786337740000103
式(7)中:J为虚拟转动惯量;为交流微电网侧角频率当前值;为交流微电网侧角频率初始值;D为虚拟阻尼系数;kudc为直流微电网下垂调节系数;Udc为直流侧母线电压当前值;Udc0为直流侧母线电压初始值;。MMC互联变换器虚拟同步机控制中可以通过直接控制交流频率和直流电压实现交直流混合微电网有功功率调节控制,均衡混合微电网有功负荷。
虚拟同步机感应内电势由两部分组成:一部分是空载时,虚拟励磁电压对于的空载电动势,一部分是由无功功率偏差产生,虚拟同步机无功励磁控制方程为:
E=E0+kq(Qref-Q) (8)
式(8)中:E为虚拟同步机感应内电势有效值;E0为励磁空载电动势有效值;kq为无功电压下垂控制系数;Qref为无功功率参考值;Q为无功功率当前值。
综合虚拟同步机有功频率控制计算得到的虚拟转子角频率和相角差δ,可以得到交直流混合微电网MMC互联变换器等效交流输出端口三相电压调制信号为:
Figure BDA0002786337740000111
如图3所示,当发生功率震荡时,角频率增加阶段ω>ω0,其中a阶段dω/dt<0,c阶段dω/dt>0,角频率增加阶段需要增加虚拟转动惯量J来约束角频率的增加,以防止ω过快增加从而造成更大超调。角频率减少阶段ω<ω0,其中b阶段dω/dt<0,d阶段dω/dt>0,需要减少虚拟转动惯量J使功率尽快恢复至稳定值。
如图4所示,本发明提出一种在角频率增加阶段,增加J同时适当减少D;在角频率减少阶段,减少J同时适当增加D的自适应控制方法。该方法可以保证系统响应速度,并且加快功率进入稳定状态。
自适应虚拟惯量和虚拟阻尼系数可表示为:
Figure BDA0002786337740000112
式(10)中:kj为虚拟惯量调节系数;J0为额定虚拟转动惯量;ω0为额定机械角频率;kd为虚拟阻尼调节系数;D0为额定虚拟阻尼系数。
在自适应虚拟同步机控制系统中,虚拟惯量调节系数kj表征虚拟转动惯量跟随频率偏差调节的能力,kd表征虚拟阻尼系数跟随频率偏差的能力。kj、kd取值过大会引起虚拟转动惯量和虚拟阻尼比超过其最大值,系统响应速度降低,不利于功率调节控制。所以需要对kj、kd进行限制,满足自适应虚拟同步机控制系统响应速度要求。因此在选择kj、kd时,需要充分的综合考虑系统动态性能,以及整体阻尼需求,从而做出灵活选择。
为了得到的额定虚拟参数取值范围,需要建立小信号模型。
根据式(3)中有功功率表达式,且交流电网电压矢量与MMC互联变换器交流侧电压矢量
Figure BDA0002786337740000121
之间的相角差δ很小,则有:
Figure BDA0002786337740000122
式(3)中有功功率又可以表示为:
Figure BDA0002786337740000123
对上式分别求一阶、二阶导数:
Figure BDA0002786337740000124
将式(12)和式(5)带入式(7)中可得:
Figure BDA0002786337740000125
对式(13)建立小信号模型并进行进行拉氏变换可得:
Figure BDA0002786337740000126
特征根为:
Figure BDA0002786337740000127
为了保证控制系统稳定,需要两个特征根都位于复平面的左半部分,又因为虚拟同步机阻尼系数D恒为正,为保证特征根实部为负,需保持虚拟惯量J≥0。在功率震荡过程中,虚拟同步机输出功率响应特性可等效为一个典型的二阶传递函数,根据式(15)可以得到系统自然震荡转速和阻尼比为:
Figure BDA0002786337740000128
参考同步发电机振荡频率:0.628rad/s≤ωS≤15.7rad/s,则额定虚拟惯量为:
Figure BDA0002786337740000131
在交直流混合微电网MMC互联变换器虚拟同步机控制系统中,为获得较快响应速度和较少超调量,在考虑阻尼系数的基础上,利用最优二阶系统分析法,取阻尼比和额定阻尼系数为:
Figure BDA0002786337740000132
如图5所示,为了验证本文所提交直流混合微电网MMC互联变换器自适应虚拟同步机控制的有效性,交直流混合微电网仿真模型采用Matlab/Simulink计算机仿真平台搭建模型,MMC互联变换器自适应虚拟同步机控制原理为:首先通过获取交流侧的母线电压与电流值,接着经过功率计算环节获取MMC互联变换器实际产生有功功率P与实际产生的无功功率Q,然后将功率作为参考信号,送至自适应虚拟同步机控制系统中,自适应虚拟同步机控制输出电压矢量,通过驱动电路控制MMC互联变换器内的IGBT导通、关断。仿真参数如下表:
表4-1仿真参数
Figure BDA0002786337740000133
Figure BDA0002786337740000141
如图6所示,MMC互联变换器工作于逆变模式,0.5s时刻,功率由MMC互联变换器从直流侧向交流侧传输,交流微电网负荷增加3kW,MMC互联变换器流过+3kW有功功率,MMC互联变换器逆变工作状态输出有功功率上升,采用常规虚拟同步机控制得到的响应功率出现大范围向上波动,且峰值较高,功率震荡严重,而采用本文所提自适应虚拟参数虚拟同步机控制得到的响应功率波形波动减少很多,峰值相应减少,可更快稳定至功率目标值,解决了功率震荡问题。1s时刻,3kW负荷切除,采用常规虚拟同步机控制得到的响应功率出现大范围向下波动,功率震荡严重,极易造成设备脱网,而采用本文所提自适应虚拟参数虚拟同步机控制得到的响应输出功率波形更为平滑,输出功率波形几乎无超调,很快稳定至目标功率,功率震荡趋近于0,调节时间进一步缩短,更适用于交直流混合微电网MMC互联变换器功率控制中。
如图7所示,0.5s时刻,当交流微电网负荷增加3kW,交流频率下降响应负荷增加,且频率调节波动明显,采用本文所提自适虚拟同步机控制频率能够稳定至目标值。1s时刻交流微电网3kW负荷切除,能够很快稳定至目标值。
如图8所示,本文所提自适应虚拟惯量J可以随MMC互联变换器功率变换而自适应调节,功率波动偏差增大,J相应增大;功率波动偏差减小,J相应减小。自适应虚拟惯量调节具有连续光滑特性,可以减少交直流混合微电网MMC互联变换器功率调节过程中的功率震荡问题,加快功率稳定至目标值,有效提高了混合微电网的暂态稳定性能。
如图9所示,自适应阻尼系数可以随MMC互联变换器功率变换而自适应调节,功率波动偏差增大,阻尼系数相应增大;功率波动偏差减小,阻尼系数相应减小。自适应阻尼系数调节具有连续光滑特性,可同时满足有功功率和频率超调较小的要求。
如图9所示,MMC互联变换器工作于整理模式,0.5s时刻,直流微电网负荷增加3kW,功率由MMC互联变换器从交流侧向直流侧传输,MMC互联变换器整理工作状态输出有功功率下降,采用常规虚拟同步机控制得到的响应功率出现大范围向下波动,波动范围较大,功率震荡严重,极易造成微电网内部设备脱网。而采用本文所提自适应虚拟参数虚拟同步机控制得到的响应功率波形波动减少很多,波动得到了很好的控制,可更快稳定至功率目标值,解决了功率震荡问题。1s时刻,3kW负荷切除,虚拟同步机响应负荷变化,输出功率上升,采用常规虚拟同步机控制得到的响应功率出现大范围向上波动,峰值较高且功率震荡严重,而采用本文所提自适应虚拟参数虚拟同步机控制得到的响应输出功率波形更为平滑,输出功率波形几乎无超调,很快稳定至目标功率,功率震荡趋近于0,调节时间进一步缩短,更适用于MMC互联变换器功率调节控制系统中。
如图10所示,0.5s时刻,当直流微电网负荷增加5kW,采用传统虚拟同步机控制频率调节波动明显,采用本文所提自适虚拟同步机控制的频率较为缓慢过渡稳定值,避免频率瞬时跌落而产生的换流器逆变失败。1s时刻直流微电网3kW负荷切除,采用常规虚拟同步机控制的频率响应明显较长,1.2s才能稳定,频率调节过程中快速的恢复且稳定至目标值。
如图11所示,MMC互联变换器工作于整流状态时,自适应虚拟惯量J可以随MMC互联变换器功率变换而自适应调节,功率波动偏差增大,J相应增大;功率波动偏差减小,J相应减小。自适应虚拟惯量调节具有连续光滑特性,可以减少交直流混合微电网MMC互联变换器功率调节过程中的功率震荡问题。
如图12所示,MMC互联变换器工作于整流状态时,自适应阻尼系数可以随MMC互联变换器功率变换而自适应调节,功率波动偏差增大,阻尼系数相应增大;功率波动偏差减小,阻尼系数相应减小。自适应阻尼系数调节具有连续光滑特性,可同时满足有功功率和频率超调较小的要求。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何限制,凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化,均仍属于本发明技术方案的保护范围内。图2为基于功率传输原则的MMC互联变换器虚拟同步机控制框图。

Claims (10)

1.一种抑制MMC互联变换器虚拟同步机功率震荡的方法,其特征在于,包括以下步骤:
1)以交直流混合微电网MMC互联变换器交流侧角度出发,得到电压方程;
3)根据步骤1)MMC互联变换器交流侧电压方程,建立MMC互联变换器功率输出方程,并进行化简;
3)根据交直流混合微电网MMC互联变换器控制系统中,交流微电网、直流微电网瞬时有功功率变化量相同,建立MMC互联变换器功率调节关系式;
4)将步骤3)MMC互联变换器功率调节关系式中交流微电网有功出力分解为稳态功率变化量、瞬时功率变化量,直流微电网下垂控制模型有功出力调节量分解为稳态功率变化量、动态功率变化量;
5)根据步骤3)MMC互联变换器功率调节关系式和步骤4)交流微电网、直流微电网有功出力分解式得到基于功率传输原则的虚拟同步机控制机械方程和虚拟同步机无功励磁控制方程;
6)根据发生功率震荡时,角频率变化情况,将步骤5)虚拟同步机控制机械方程中的虚拟惯量和虚拟阻尼修改为自适应虚拟惯量和自适应虚拟阻尼系数;
7)对步骤2)MMC互联变换器功率输出化简方程中的有功功率求一阶、二阶导数;
8)将步骤7)中MMC互联变换器功输出有功功率一阶、二阶导数和步骤4)交流微电网有功出力分解表达式带入步骤5)基于功率传输原则的虚拟同步机控制机械方程中,得到虚拟同步机参考输出功率表达式;
9)对步骤8)中虚拟同步机参考输出功率表达式建立小信号模型并进行拉氏变换,计算其特征根;
10)分析步骤9)虚拟同步机参考输出功率小信号模型和其特征根,得到系统自然震荡转速和阻尼比;
11)参考同步发电机振荡频率,得到额定虚拟惯量取值范围,在考虑阻尼系数的基础上,利用最优二阶系统分析法,考虑阻尼系数的基础上,设置步骤10)阻尼比和额定阻尼系数;
12)根据步骤11)额定虚拟惯量取值范围和阻尼比和额定阻尼系数,设计步骤6)自适应虚拟惯量和自适应虚拟阻尼系数,带入步骤5)基于功率传输原则的虚拟同步机控制机械方程中,起到抑制MMC互联变换器虚拟同步机功率震荡的作用。
2.根据权利要求1所述的一种抑制MMC互联变换器虚拟同步机功率震荡的方法,其特征在于,步骤1)以交直流混合微电网MMC互联变换器交流侧角度出发,得到电压方程,
Figure FDA0002786337730000021
其中:
Figure FDA0002786337730000022
为交流电网电压e0对应的矢量;
Figure FDA0002786337730000023
为MMC互联变换器交流侧电压Uac对应的矢量;
Figure FDA0002786337730000024
为交流电网电流i0对应的矢量。
3.根据权利要求2所述的一种抑制MMC互联变换器虚拟同步机功率震荡的方法,其特征在于,步骤2)的具体实现方法为:根据步骤1)MMC互联变换器交流侧电压方程,建立MMC互联变换器功率传输方程:
Figure FDA0002786337730000025
其中:Rf、Xf为滤波电路电阻值及感抗;δ为交流电网电压矢量
Figure FDA00027863377300000212
与MMC互联变换器交流侧电压矢量
Figure FDA00027863377300000213
之间的相角差;交流电网电压e0类比于同步电机电枢电动势;MMC互联变换器交流侧电压Uac类比于同步电机端电压;在交直流混合微电网MMC互联变换器控制系统中,通过控制δ正负,实现功率双向流动,当δ>0,
Figure FDA0002786337730000026
超前
Figure FDA0002786337730000027
MMC互联变换器运行于逆变模式,功率由直流微电网向交流微电网传输;当δ<0,
Figure FDA0002786337730000028
滞后
Figure FDA0002786337730000029
MMC互联变换器运行于整流模式,功率由交流微电网向直流微电网传输;当δ=0,
Figure FDA00027863377300000210
Figure FDA00027863377300000211
同相位,交流微电网、直流微电网之间无功率交换;MMC互联变换器,满足Rf<<Xf,对MMC互联变换器功率传输方程进行化简:
Figure FDA0002786337730000031
交流电网电压矢量与MMC互联变换器交流侧电压矢量
Figure FDA0002786337730000032
之间的相角差δ很小,则有:
Figure FDA0002786337730000033
有功功率表示为:
Figure FDA0002786337730000034
4.根据权利要求3所述的一种抑制MMC互联变换器虚拟同步机功率震荡的方法,其特征在于,步骤3)的具体实现方法为:根据交直流混合微电网MMC互联变换器控制系统中,交流微电网、直流微电网瞬时有功功率变化量相同,建立MMC互联变换器功率调节关系式:Pacref-Pac=Pdc-Pdcref=ΔP;
其中:Pac为交流微电网输出功率实际值;Pacef为交流微电网输出功率实际值;Pdc为直流微电网输出功率实际值;Pdcref为直流微电网输出功率实际值;ΔP为MMC互联变换器功率调节量。
5.根据权利要求4所述的一种抑制MMC互联变换器虚拟同步机功率震荡的方法,其特征在于,步骤4)的具体实现方法为:将步骤3)MMC互联变换器功率调节关系式中交流微电网有功出力分解为稳态功率变化量、瞬时功率变化量:
Figure FDA0002786337730000035
其中:kω为交流电网下垂调节系数;kω(ω-ω0)是稳态有功变化量;
Figure FDA0002786337730000036
为虚拟转子惯量吸收或者发出的瞬时有功功率,在虚拟同步机控制中为交流频率提供惯性环节;直流微电网下垂控制模型有功出力调节量分解为稳态功率变化量、动态功率变化量:
Figure FDA0002786337730000037
其中:kudc为直流电网下垂调节系数;Udc为直流母线电压实际值;Udc0为直流母线电压初始值;Cdc为直流侧电容值;kudc(Udc-Udc0)是稳态有功变化量,
Figure FDA0002786337730000038
是直流电容的充放电功率,属于动态功率波动。
6.根据权利要求5所述的一种抑制MMC互联变换器虚拟同步机功率震荡的方法,其特征在于,步骤5)的具体实现方法为:根据步骤3)MMC互联变换器功率调节关系式和步骤4)交流微电网、直流微电网有功出力分解式得到基于功率传输原则的虚拟同步机控制机械方程和虚拟同步机无功励磁控制方程
Figure FDA0002786337730000041
其中:J为虚拟转动惯量;为交流微电网侧角频率当前值;为交流微电网侧角频率初始值;D为虚拟阻尼系数;kudc为直流微电网下垂调节系数;Udc为直流侧母线电压当前值;Udc0为直流侧母线电压初始值;Cdc为直流侧电容值;拟同步机感应内电势由两部分组成:一部分是空载时,虚拟励磁电压对于的空载电动势,一部分是由无功功率偏差产生,虚拟同步机无功励磁控制方程为:E=E0+kq(Qref-Q);
其中:E为虚拟同步机感应内电势有效值;E0为励磁空载电动势有效值;kq为无功电压下垂控制系数;Qref为无功功率参考值;Q为无功功率当前值;综合虚拟同步机有功频率控制计算得到的虚拟转子角频率和相角差δ,得到交直流混合微电网MMC互联变换器等效交流输出端口三相电压调制信号为:
Figure FDA0002786337730000042
7.根据权利要求6所述的一种抑制MMC互联变换器虚拟同步机功率震荡的方法,其特征在于,步骤6)的具体实现方法为:根据发生功率震荡时,角频率变化情况,在角频率增加阶段,增加J同时适当减少D;在角频率减少阶段,减少J同时适当增加D的自适应控制方法;该方法保证系统响应速度,并且加快功率进入稳定状态将步骤5)虚拟同步机控制机械方程中的虚拟惯量和虚拟阻尼修改为自适应虚拟惯量和自适应虚拟阻尼系数:
Figure FDA0002786337730000043
其中:kj为虚拟惯量调节系数;J0为额定虚拟转动惯量;ω0为额定机械角频率;kd为虚拟阻尼调节系数;D0为额定虚拟阻尼系数;在自适应虚拟同步机控制系统中,虚拟惯量调节系数kj表征虚拟转动惯量跟随频率偏差调节的能力,kd表征虚拟阻尼系数跟随频率偏差的能力。
8.根据权利要求7所述的一种抑制MMC互联变换器虚拟同步机功率震荡的方法,其特征在于,步骤7)的具体实现方法为:对步骤2)MMC互联变换器功率输出化简方程中的有功功率求一阶、二阶导数:
Figure FDA0002786337730000051
9.根据权利要求8所述的一种抑制MMC互联变换器虚拟同步机功率震荡的方法,其特征在于,步骤8)的具体实现方法为:将步骤7)中MMC互联变换器功输出有功功率一阶、二阶导数和步骤4)交流微电网有功出力分解表达式带入步骤5)基于功率传输原则的虚拟同步机控制机械方程中,得到虚拟同步机参考输出功率表达式:
Figure FDA0002786337730000052
10.根据权利要求9所述的一种抑制MMC互联变换器虚拟同步机功率震荡的方法,其特征在于,步骤9)的具体实现方法为:对步骤8)中虚拟同步机参考输出功率表达式建立小信号模型并进行拉氏变换:
Figure FDA0002786337730000053
计算其特征根:
Figure FDA0002786337730000054
步骤10)的具体实现方法为:分析步骤9)虚拟同步机参考输出功率小信号模型和其特征根,得到系统自然震荡转速和阻尼比:
Figure FDA0002786337730000061
步骤11)的具体实现方法为:参考同步发电机振荡频率:0.628rad/s≤ωS≤15.7rad/s,得到额定虚拟惯量取值范围:
Figure FDA0002786337730000062
在考虑阻尼系数的基础上,利用最优二阶系统分析法,考虑阻尼系数的基础上,设置步骤10)阻尼比和额定阻尼系数:
Figure FDA0002786337730000063
步骤12)的具体实现方法为:根据步骤11)额定虚拟惯量取值范围:
Figure FDA0002786337730000064
和阻尼比和额定阻尼系数:
Figure FDA0002786337730000065
设计步骤6)自适应虚拟惯量J和自适应虚拟阻尼系数D,带入步骤5)基于功率传输原则的虚拟同步机控制机械方程中,起到抑制MMC互联变换器虚拟同步机功率震荡的作用。
CN202011299319.0A 2020-11-18 2020-11-18 一种抑制mmc互联变换器虚拟同步机功率震荡的方法 Active CN112467790B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011299319.0A CN112467790B (zh) 2020-11-18 2020-11-18 一种抑制mmc互联变换器虚拟同步机功率震荡的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011299319.0A CN112467790B (zh) 2020-11-18 2020-11-18 一种抑制mmc互联变换器虚拟同步机功率震荡的方法

Publications (2)

Publication Number Publication Date
CN112467790A true CN112467790A (zh) 2021-03-09
CN112467790B CN112467790B (zh) 2023-02-24

Family

ID=74836701

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011299319.0A Active CN112467790B (zh) 2020-11-18 2020-11-18 一种抑制mmc互联变换器虚拟同步机功率震荡的方法

Country Status (1)

Country Link
CN (1) CN112467790B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113193569A (zh) * 2021-04-22 2021-07-30 南方电网科学研究院有限责任公司 具有短期频率支撑和直流侧振荡抑制功能的电容能量控制方法
CN113507143A (zh) * 2021-07-08 2021-10-15 东北电力大学 基于改进vsg技术的混合微电网ic自适应控制策略
CN113629783A (zh) * 2021-07-20 2021-11-09 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 一种含功率反馈的自适应惯量下垂控制方法
CN115036942A (zh) * 2022-08-10 2022-09-09 西安热工研究院有限公司 一种混合储能vsg频率振荡抑制方法和装置

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106972519A (zh) * 2017-04-27 2017-07-21 湖南大学 直流输电系统直流侧谐振的有源阻尼控制装置及方法
US20170261280A1 (en) * 2016-03-10 2017-09-14 Sapphire Defense Group LLC Enhanced metal-metal-matrix composite weapon barrels and ways of making the same
CN109446653A (zh) * 2018-10-30 2019-03-08 华中科技大学 混合型mmc状态空间、稳态解析和小信号模型的建模方法
CN109861246A (zh) * 2018-12-24 2019-06-07 燕山大学 一种基于vsg的光伏微网动态频率稳定控制方法
CN109861279A (zh) * 2019-01-24 2019-06-07 太原理工大学 一种适用于虚拟同步发电机的转动惯量自适应控制方法
CN110120676A (zh) * 2019-06-18 2019-08-13 国网新疆电力有限公司电力科学研究院 基于模拟同步电机特性的mmc-hvdc换流站功率控制方法及系统

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20170261280A1 (en) * 2016-03-10 2017-09-14 Sapphire Defense Group LLC Enhanced metal-metal-matrix composite weapon barrels and ways of making the same
CN106972519A (zh) * 2017-04-27 2017-07-21 湖南大学 直流输电系统直流侧谐振的有源阻尼控制装置及方法
CN109446653A (zh) * 2018-10-30 2019-03-08 华中科技大学 混合型mmc状态空间、稳态解析和小信号模型的建模方法
CN109861246A (zh) * 2018-12-24 2019-06-07 燕山大学 一种基于vsg的光伏微网动态频率稳定控制方法
CN109861279A (zh) * 2019-01-24 2019-06-07 太原理工大学 一种适用于虚拟同步发电机的转动惯量自适应控制方法
CN110120676A (zh) * 2019-06-18 2019-08-13 国网新疆电力有限公司电力科学研究院 基于模拟同步电机特性的mmc-hvdc换流站功率控制方法及系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
ALIAKBAR JAMSHIDIFAR: "Small-Signal Dynamic DQ Model of Modular Multilevel Converter for System Studies", 《IEEE TRANSACTIONS ON POWER DELIVERY》 *
JIANGBO CHEN: "Harmonic Suppression Scheme of L-type New Energy Grid-connected Inverter Based on Harmonic Injection Method", 《2019 IEEE 8TH INTERNATIONAL CONFERENCE ON ADVANCED POWER SYSTEM AUTOMATION AND PROTECTION (APAP)》 *
张宝锋: "电化学储能在新能源发电侧的应用分析", 《热力发电》 *
赵 峰: "孤网变频系统电量频率跟踪方法", 《电工技术》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113193569A (zh) * 2021-04-22 2021-07-30 南方电网科学研究院有限责任公司 具有短期频率支撑和直流侧振荡抑制功能的电容能量控制方法
CN113507143A (zh) * 2021-07-08 2021-10-15 东北电力大学 基于改进vsg技术的混合微电网ic自适应控制策略
CN113507143B (zh) * 2021-07-08 2022-09-20 东北电力大学 基于改进vsg技术的混合微电网ic自适应控制策略
CN113629783A (zh) * 2021-07-20 2021-11-09 内蒙古电力(集团)有限责任公司内蒙古电力科学研究院分公司 一种含功率反馈的自适应惯量下垂控制方法
CN115036942A (zh) * 2022-08-10 2022-09-09 西安热工研究院有限公司 一种混合储能vsg频率振荡抑制方法和装置

Also Published As

Publication number Publication date
CN112467790B (zh) 2023-02-24

Similar Documents

Publication Publication Date Title
CN112467790B (zh) 一种抑制mmc互联变换器虚拟同步机功率震荡的方法
WO2022077847A1 (zh) 一种用于混合微电网mmc互联变换器的虚拟同步机控制方法
CN112270075B (zh) 一种mmc互联变换器虚拟同步机小信号模型分析方法
CN112398166B (zh) 储能一次调频虚拟同步机参数分析方法
CN112003323A (zh) 一种利用自适应虚拟参数提高风电并网一次调频性能的方法
CN105958552B (zh) 适应不平衡电网及负载条件的虚拟同步发电机控制方法
CN102074967B (zh) 一种具有并网特性的储能型风电场控制方法
CN112467784B (zh) 一种混合微网换流器自适应虚拟同步机控制方法
CN110429655B (zh) 基于同步机三阶模型的储能单元主动支撑控制方法及系统
CN103311949B (zh) 大功率光伏逆变器的低电压穿越控制方法
CN108429289B (zh) 一种基于虚拟同步发电机的控制方法及系统
CN111900762B (zh) 一种自适应vsg微电网逆变器控制方法及系统
CN104638679A (zh) 一种采用自适应调节的虚拟同步发电机频率控制方法
Torres et al. Inverter-based diesel generator emulator for the study of frequency variations in a laboratory-scale autonomous power system
CN105552948B (zh) 一种基于柔性直流输电系统的电网调频方法
CN108199396A (zh) 储能逆变器虚拟励磁闭环控制系统及其设计方法
WO2022198764A1 (zh) 海上风电并网系统的控制方法和系统
CN108695885B (zh) 对称电压暂降下改进的虚拟同步发电机控制方法
CN112467789B (zh) 一种基于功率传输原则的混合微网虚拟同步机控制方法
CN112398167A (zh) 一种提高微网储能一次调频性能的方法
Tan et al. Research on primary frequency regulation of wind turbine based on new nonlinear droop control
CN104808147A (zh) 一种三相短路下双馈发电机仿真方法及其系统
CN111327074A (zh) 一种基于虚拟同步机技术的直流侧混合储能协调控制方法
CN109066735B (zh) 一种不平衡电网电压下的双馈风力发电系统及其控制方法
CN112510696B (zh) 微电网暂态电压稳定性调控系统及调控方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant