CN113783208B - 一种双馈型机组风电场虚拟惯量控制方法和系统 - Google Patents

一种双馈型机组风电场虚拟惯量控制方法和系统 Download PDF

Info

Publication number
CN113783208B
CN113783208B CN202111323003.5A CN202111323003A CN113783208B CN 113783208 B CN113783208 B CN 113783208B CN 202111323003 A CN202111323003 A CN 202111323003A CN 113783208 B CN113783208 B CN 113783208B
Authority
CN
China
Prior art keywords
unit
inertia
power
coefficient
wind
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.)
Active
Application number
CN202111323003.5A
Other languages
English (en)
Other versions
CN113783208A (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.)
China Electric Power Research Institute Co Ltd CEPRI
Original Assignee
China Electric Power Research Institute Co Ltd CEPRI
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 China Electric Power Research Institute Co Ltd CEPRI filed Critical China Electric Power Research Institute Co Ltd CEPRI
Priority to CN202111323003.5A priority Critical patent/CN113783208B/zh
Publication of CN113783208A publication Critical patent/CN113783208A/zh
Application granted granted Critical
Publication of CN113783208B publication Critical patent/CN113783208B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/241—The oscillation concerning frequency
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects

Landscapes

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

Abstract

本发明提供了一种双馈型机组风电场虚拟惯量控制方法和系统，包括：基于电网系统的实际频率计算双馈型机组风电场的一次调频附加功率；基于各机组的轮毂中心风速和风轮转速分配一次调频附加功率，得到各机组的惯性额外功率和下垂额外功率；根据各机组的惯性额外功率、下垂额外功率和发电机转速，控制各机组在匹配惯量支撑能力前提下协同参与调频；本发明基于各机组的轮毂中心风速和风轮转速分配一次调频附加功率，使得各机组能够提供不同的转子动能惯量，避免了过度惯量支撑导致一次调频失败和缓解了增加的机械疲劳损伤。

Description

一种双馈型机组风电场虚拟惯量控制方法和系统
技术领域
本发明属于新能源风力发电场电网频率调控技术领域，具体涉及一种双馈型机组风电场虚拟惯量控制方法和系统。
背景技术
提高可再生能源在发电中的比例是解决环境污染和化石能源短缺问题的重要手段和实现碳达峰、碳中和的关键途径。风力发电作为可再生能源的重要形式之一近年来虽然日渐重视并且得到了迅速发展，但由于风电机组输出功率的波动性使得电网源荷之间动态不平衡，导致电网频率偏差增大，严重危及了电力系统的安全稳定运行。
为解决电网中不断升高的风电穿透率问题，提高电力系统的调频能力，需对风电变换器的优化控制，促使风电机组参与频率支撑与惯量响应。目前，针对双馈风电机组的频率支撑控制策略包括虚拟惯量控制、功率备用、附加储能等。但功率备用降低了机组的发电效率，附加储能又增加了发电成本，而良好的虚拟惯量控制技术兼具效率与成本成为调频的主流技术之一，各国也因此制定了相关的风电惯量及调频能力规范。
当前虚拟惯量控制策略更多的关注点在单台机组的调控能力，比如如何防止功率二次跌落，增加惯量支撑能力，以满足调频规范的要求，并未考虑风电场各机组调频的协同作用，最大化风电场各风电机组转子动能所能提供的不同的虚拟惯量支撑能力，并且最小化降低传动链机械疲劳损伤，这也是本发明主要解决的问题。
在一种风电机组一次调频与虚拟惯量协调控制方法及装置中：
根据电网频率确定风电场中各风电机组虚拟惯量控制的附加功率：
Figure 28013DEST_PATH_IMAGE001
根据电网频率确定风电场一次调频控制的附加功率和风电场中各风电机组虚拟惯量控制的附加功率；根据风电场一次调频控制的附加功率确定风电场中各风电机组的功率参考值和桨距角控制值；利用风电场中各风电机组虚拟惯量控制的附加功率对风电场中各风电机组的功率参考值进行修正，获得风电场中各风电机组的功率控制值；按照风电场中各风电机组的桨距角控制值和风电机组的功率控制值调节风电场中各风电机组的桨距角和功率；该发明对风电机组采用一次调频控制和虚拟惯量控制相结合的控制策略，可以改善大规模风电并网系统的频率稳定性。
该发明着重电网一次调频的基础上通过惯量控制附加功率来改善大规模风电并网系统的频率稳定性，没有根据现场实际风况区别利用机组惯量控制附加功率的能力；且该发明未考虑惯量控制的风电场内机组的机械损伤。
发明内容
为克服上述现有技术的不足，本发明提出一种双馈型机组风电场虚拟惯量控制方法，包括：
基于电网系统的实际频率计算双馈型机组风电场的一次调频附加功率；
基于各机组的轮毂中心风速和风轮转速分配所述一次调频附加功率，得到各机组的惯性额外功率和下垂额外功率；
根据各机组的惯性额外功率、下垂额外功率和发电机转速，控制各机组在匹配惯量支撑能力前提下协同参与调频。
优选的，所述根据各机组的惯性额外功率、下垂额外功率和发电机转速，控制各机组在匹配惯量支撑能力前提下协同参与调频，包括：
根据各机组的惯性额外功率、下垂额外功率和发电机转速，分别调整各机组的发电机励磁转矩；
根据各机组的风轮转速、发电机转速和发电机励磁转矩，分别为各机组的传动链施加带通滤波阻尼；
基于各机组的发电机励磁转矩和带通滤波阻尼，控制各机组在匹配惯量支撑能力前提下协同参与调频。
优选的，所述发电机励磁转矩的计算式如下：
Figure 266359DEST_PATH_IMAGE002
式中，
Figure 741203DEST_PATH_IMAGE003
为机组i的发电机励磁转矩，
Figure 165991DEST_PATH_IMAGE004
为机组i的下垂额外功率，
Figure 409891DEST_PATH_IMAGE005
为机组i 的惯性额外功率，
Figure 661880DEST_PATH_IMAGE006
为机组i当前初始功率，
Figure 624020DEST_PATH_IMAGE007
为机组i的发电机转速。
优选的，所述带通滤波阻尼的计算式如下：
Figure 986868DEST_PATH_IMAGE008
式中，
Figure 321160DEST_PATH_IMAGE009
为机组i的带通滤波阻尼，
Figure 744051DEST_PATH_IMAGE010
为机组i的风轮转速，
Figure 662329DEST_PATH_IMAGE007
为机组i的发电机 转速，
Figure 828868DEST_PATH_IMAGE011
为机组i的增益系数；
所述机组i的增益系数
Figure 516201DEST_PATH_IMAGE011
的值根据机组i的发电机励磁转矩动态调整。
优选的，所述根据各机组的风轮转速和发电机转速，分别为各机组的传动链施加带通滤波阻尼之后，且基于各机组的发电机励磁转矩和带通滤波阻尼，控制各机组在匹配惯量支撑能力前提下协同参与调频之前，还包括：
分别为各机组的传动链设计动态相角补偿，补偿控制存在的滞后性。
优选的，所述一次调频附加功率的计算式如下：
Figure 608528DEST_PATH_IMAGE012
式中，
Figure 14102DEST_PATH_IMAGE013
为一次调频附加功率，K为风电场一次调频比例系数，
Figure 718753DEST_PATH_IMAGE014
为电网系统的 参考频率，
Figure 526172DEST_PATH_IMAGE015
为电网系统的实际频率，
Figure 290866DEST_PATH_IMAGE016
为风电场频率监测死区。
优选的，所述基于各机组的轮毂中心风速和风轮转速分配所述一次调频附加功率，得到各机组的惯性额外功率和下垂额外功率，包括：
基于各机组的轮毂中心风速计算各机组的惯性系数的比例系数以及下垂系数的比例系数；
基于各机组的惯性系数的比例系数、下垂系数的比例系数以及风轮转速，将所述一次调频附加功率分配至参与频率响应控制的各机组，得到各机组的惯性额外功率和下垂额外功率。
优选的，所述各机组的惯性系数的比例系数的计算式如下：
Figure 183735DEST_PATH_IMAGE017
式中，
Figure 193542DEST_PATH_IMAGE018
为机组i的惯性系数的比例系数，v i 为机组i的轮毂处实时风速，
Figure 121046DEST_PATH_IMAGE019
为机 组的切入风速，
Figure 56641DEST_PATH_IMAGE020
为机组的切出风速。
优选的，所述各机组的下垂系数的比例系数的计算式如下：
Figure 702386DEST_PATH_IMAGE021
式中，
Figure 716217DEST_PATH_IMAGE022
为机组i的下垂系数的比例系数，a为预设常数，v i 为机组i的轮毂处实时 风速，
Figure 324659DEST_PATH_IMAGE019
为机组的切入风速，
Figure 227893DEST_PATH_IMAGE020
为机组的切出风速。
优选的，所述基于各机组的惯性系数的比例系数、下垂系数的比例系数以及风轮转速，将所述一次调频附加功率分配至参与频率响应控制的各机组，得到各机组的惯性额外功率和下垂额外功率的计算式如下：
Figure 95355DEST_PATH_IMAGE023
Figure 446964DEST_PATH_IMAGE024
Figure 349061DEST_PATH_IMAGE025
式中，
Figure 626459DEST_PATH_IMAGE013
为一次调频附加功率，
Figure 981216DEST_PATH_IMAGE026
为参与频率响应控制机组的数量，
Figure 258221DEST_PATH_IMAGE005
为机组i的惯性额外功率，
Figure 545983DEST_PATH_IMAGE004
为机组i的下垂额外功率，
Figure 761325DEST_PATH_IMAGE018
为机组i的惯性系数的比例系数，
Figure 603380DEST_PATH_IMAGE027
为 机组惯性系数，
Figure 326485DEST_PATH_IMAGE022
为机组i的下垂系数的比例系数，
Figure 203174DEST_PATH_IMAGE028
为机组下垂系数，
Figure 320909DEST_PATH_IMAGE029
为机组i的风轮 转速，
Figure 650260DEST_PATH_IMAGE030
为惯量响应机组的风轮转速下限，
Figure 645897DEST_PATH_IMAGE031
为惯量响应机组风轮转速上限，
Figure 111514DEST_PATH_IMAGE032
为机 组频率与电网系统的实际频率的差。
基于同一发明构思，本申请还提供了一种双馈型机组风电场虚拟惯量控制系统包括：一次调频附加功率模块、功率分配模块和控制模块；
所述一次调频附加功率模块，用于基于电网系统的实际频率计算双馈型机组风电场的一次调频附加功率；
所述功率分配模块，用于基于各机组的轮毂中心风速和风轮转速分配所述一次调频附加功率，得到各机组的惯性额外功率和下垂额外功率；
所述控制模块，用于根据各机组的惯性额外功率、下垂额外功率和发电机转速，控制各机组在匹配惯量支撑能力前提下协同参与调频。
优选的，所述控制模块具体用于：
根据各机组的惯性额外功率、下垂额外功率和发电机转速，分别调整各机组的发电机励磁转矩；
根据各机组的风轮转速、发电机转速和发电机励磁转矩，分别为各机组的传动链施加带通滤波阻尼；
基于各机组的发电机励磁转矩和带通滤波阻尼，控制各机组在匹配惯量支撑能力前提下协同参与调频。
优选的，所述发电机励磁转矩的计算式如下：
Figure 636036DEST_PATH_IMAGE002
式中，
Figure 688568DEST_PATH_IMAGE003
为机组i的发电机励磁转矩，
Figure 487897DEST_PATH_IMAGE004
为机组i的下垂额外功率，
Figure 808020DEST_PATH_IMAGE005
为机组i 的惯性额外功率，
Figure 34602DEST_PATH_IMAGE006
为机组i当前初始功率，
Figure 338544DEST_PATH_IMAGE007
为机组i的发电机转速。
优选的，所述带通滤波阻尼的计算式如下：
Figure 440099DEST_PATH_IMAGE008
式中，
Figure 614728DEST_PATH_IMAGE009
为机组i的带通滤波阻尼，
Figure 12211DEST_PATH_IMAGE010
为机组i的风轮转速，
Figure 272291DEST_PATH_IMAGE007
为机组i的发电机 转速，
Figure 413423DEST_PATH_IMAGE011
为机组i的增益系数；
所述机组i的增益系数
Figure 944024DEST_PATH_IMAGE011
的值根据机组i的发电机励磁转矩动态调整。
优选的，所述一次调频附加功率的计算式如下：
Figure 512408DEST_PATH_IMAGE012
式中，
Figure 259784DEST_PATH_IMAGE013
为一次调频附加功率，K为风电场一次调频比例系数，
Figure 939027DEST_PATH_IMAGE014
为电网系统的 参考频率，
Figure 822670DEST_PATH_IMAGE015
为电网系统的实际频率，
Figure 561956DEST_PATH_IMAGE016
为风电场频率监测死区。
优选的，所述功率分配模块具体用于：
基于各机组的轮毂中心风速计算各机组的惯性系数的比例系数以及下垂系数的比例系数；
基于各机组的惯性系数的比例系数、下垂系数的比例系数以及风轮转速，将所述一次调频附加功率分配至参与频率响应控制的各机组，得到各机组的惯性额外功率和下垂额外功率。
优选的，所述各机组的惯性系数的比例系数的计算式如下：
Figure 295163DEST_PATH_IMAGE017
式中，
Figure 778097DEST_PATH_IMAGE018
为机组i的惯性系数的比例系数，v i 为机组i的轮毂处实时风速，
Figure 47404DEST_PATH_IMAGE019
为机 组的切入风速，
Figure 223171DEST_PATH_IMAGE020
为机组的切出风速。
优选的，所述各机组的下垂系数的比例系数的计算式如下：
Figure 679560DEST_PATH_IMAGE021
式中，
Figure 202070DEST_PATH_IMAGE022
为机组i的下垂系数的比例系数，a为预设常数，v i 为机组i的轮毂处实时 风速，
Figure 325884DEST_PATH_IMAGE019
为机组的切入风速，
Figure 406973DEST_PATH_IMAGE020
为机组的切出风速。
优选的，所述基于各机组的惯性系数的比例系数、下垂系数的比例系数以及风轮转速，将所述一次调频附加功率分配至参与频率响应控制的各机组，得到各机组的惯性额外功率和下垂额外功率的计算式如下：
Figure 616237DEST_PATH_IMAGE023
Figure 440974DEST_PATH_IMAGE024
Figure 941266DEST_PATH_IMAGE025
式中，
Figure 193256DEST_PATH_IMAGE013
为一次调频附加功率，
Figure 624238DEST_PATH_IMAGE026
为参与频率响应控制机组的数量，
Figure 721507DEST_PATH_IMAGE005
为机组i的惯性额外功率，
Figure 819913DEST_PATH_IMAGE004
为机组i的下垂额外功率，
Figure 744269DEST_PATH_IMAGE018
为机组i的惯性系数的比例系数，
Figure 662546DEST_PATH_IMAGE027
为 机组惯性系数，
Figure 829085DEST_PATH_IMAGE022
为机组i的下垂系数的比例系数，
Figure 781998DEST_PATH_IMAGE028
为机组下垂系数，
Figure 375790DEST_PATH_IMAGE029
为机组i的风轮 转速，
Figure 781364DEST_PATH_IMAGE030
为惯量响应机组的风轮转速下限，
Figure 984550DEST_PATH_IMAGE031
为惯量响应机组风轮转速上限，
Figure 791969DEST_PATH_IMAGE032
为机 组频率与电网系统的实际频率的差。
本发明还提供一种计算机设备，包括：一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行时，实现如前所述的方法。
本发明还提供一种计算机可读存储介质，其上存有计算机程序，所述计算机程序被执行时，实现如前所述的方法。
与最接近的现有技术相比，本发明具有的有益效果如下：
本发明提供了一种双馈型机组风电场虚拟惯量控制方法和系统，包括：基于电网系统的实际频率计算双馈型机组风电场的一次调频附加功率；基于各机组的轮毂中心风速和风轮转速分配一次调频附加功率，得到各机组的惯性额外功率和下垂额外功率；根据各机组的惯性额外功率、下垂额外功率和发电机转速，控制各机组在匹配惯量支撑能力前提下协同参与调频；本发明基于各机组的轮毂中心风速和风轮转速分配一次调频附加功率，使得各机组能够提供不同的转子动能惯量，避免了过度惯量支撑导致一次调频失败和缓解了增加的机械疲劳损伤。
本发明进一步在虚拟惯量控制的同时，对传动链施加附加阻尼控制，减轻了机组的机械疲劳损失。
附图说明
图1为本发明提供的一种双馈型机组风电场虚拟惯量控制方法流程示意图；
图2为本发明提供的一种双馈型机组风电场虚拟惯量控制系统结构示意图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步的详细说明。
本发明提供了一种双馈型机组风电场虚拟惯量控制方法和系统，根据不同风速下机组的转子动能惯量支撑能力并不相同，可通过动态分配不同机组的额外功率，避免风电场内不同风轮风速的风电机组提供相同的惯量支撑，导致转速较低、风轮风速较小的机组，过度提供惯量支撑所需能量，使得该风电机组功率支撑达不到要求或转速恢复失败，功率二次跌落，最终惯量调频失败；另一方面，缓解了过度支撑对传动链的机械冲击，减轻了传动链的机械疲劳损伤；相应的转速较高、风轮风速较大的机组虚拟惯量支撑后还有一部分能量富余，造成一定程度的惯量浪费。
实施例1：
本发明提供的一种双馈型机组风电场虚拟惯量控制方法流程示意图如图1所示，包括：
步骤1：基于电网系统的实际频率计算双馈型机组风电场的一次调频附加功率；
步骤2：基于各机组的轮毂中心风速和风轮转速分配一次调频附加功率，得到各机组的惯性额外功率和下垂额外功率；
步骤3：根据各机组的惯性额外功率、下垂额外功率和发电机转速，控制各机组在匹配惯量支撑能力前提下协同参与调频。
本发明提供了一种双馈型机组风电场虚拟惯量控制方法，根据风电场各机组不同惯量支撑能力的动态分配额外功率以解决电网惯性降低、调频能力不足和缓解传动链机械疲劳损伤等问题，具体控制方法如下：
步骤1：确定风电场参与一次调频所需的附加功率
Figure 556662DEST_PATH_IMAGE013
:
风电场主控监测电网系统的实际频率
Figure 449532DEST_PATH_IMAGE015
，并根据电网系统的参考频率
Figure 957874DEST_PATH_IMAGE014
与风 电场频率监测死区
Figure 121264DEST_PATH_IMAGE016
下达一次调频附加功率
Figure 56859DEST_PATH_IMAGE013
指令：
Figure 437025DEST_PATH_IMAGE033
（1）
式（1）中K为风电场一次调频比例系数。
步骤2具体包括：
2-1：确定风电场内参与频率响应控制的机组数量及各机组调频变系数比例:
风电机组主控系统将测量的风轮转速反馈给风电场控制系统，只有当风轮转速在
Figure 749057DEST_PATH_IMAGE030
与
Figure 265489DEST_PATH_IMAGE031
之间的风电机组才会收到风电场主控下达的参与频率响应控制指令。其中，
Figure 870521DEST_PATH_IMAGE030
为惯量响应机组的风轮转速下限，
Figure 737983DEST_PATH_IMAGE031
为惯量响应机组风轮转速上限，设参与频率响 应的机组数量为
Figure 588127DEST_PATH_IMAGE026
。
根据采集的轮毂中心风速确定机组i惯性系数的比例系数
Figure 959065DEST_PATH_IMAGE018
：
Figure 502042DEST_PATH_IMAGE034
（2）
和机组i下垂系数的比例系数
Figure 591221DEST_PATH_IMAGE022
：
Figure 480942DEST_PATH_IMAGE035
（3）
式（2）（3）中，
Figure 971966DEST_PATH_IMAGE019
，
Figure 685844DEST_PATH_IMAGE020
分别为机组的切入、切出风速，
Figure 262319DEST_PATH_IMAGE036
为机组i轮毂处实时风速。
2-2：动态分配参与频率响应控制机组的附加功率:
Figure 454266DEST_PATH_IMAGE037
（4）
Figure 829490DEST_PATH_IMAGE038
（5）
Figure 448690DEST_PATH_IMAGE039
（6）
式（4）（5）（6）中
Figure 512461DEST_PATH_IMAGE026
为参与频率响应控制机组的数量，
Figure 508099DEST_PATH_IMAGE027
为机组惯性系数，
Figure 973715DEST_PATH_IMAGE028
为机 组下垂系数，
Figure 265282DEST_PATH_IMAGE005
为机组i惯性额外功率，
Figure 816349DEST_PATH_IMAGE004
为机组i下垂额外功率，
Figure 350098DEST_PATH_IMAGE029
为参与惯量响应的 机组i的风轮转速，
Figure 670221DEST_PATH_IMAGE032
为机组频率与电网系统的实际频率的差。
步骤3具体为：依据分配的支撑功率风电场各机组协同参与调频。
风电场的主控系统与双馈机组的主控系统之间交换每个风电机组的运行状态，包括机组风轮风速，转速，功率，频率偏差及频率偏差变化率等。风电场主控分配各个机组不同的调频附加功率。风电机组主控接到风电场指令后，迅速通过调节发电机励磁转矩，使其功率增加或减少附加功率，共同调节系统频率。单台机组参与调频的方法如下：
1）根据分配而来的附加功率，调整给定机组i发电机励磁转矩
Figure 631224DEST_PATH_IMAGE003
；
Figure 669587DEST_PATH_IMAGE040
（7）
式（7）中，
Figure 239984DEST_PATH_IMAGE041
为机组i当前初始功率，
Figure 680192DEST_PATH_IMAGE042
为机组i发电机转速。
2）为传动链设计带通滤波阻尼，减缓扭振；
由于电磁转矩脉动将给传动链强烈的扭振冲击，为减缓传动链扭振，在给定的发电机电磁转矩的基础增加适当的与传动链扭转速度相反，且成正比的转矩波动，即：
Figure 812096DEST_PATH_IMAGE043
（8）
式（8）中，
Figure 72176DEST_PATH_IMAGE044
为机组i转矩波动即带通滤波阻尼，
Figure 947728DEST_PATH_IMAGE011
为机组i的增益系数，
Figure 242444DEST_PATH_IMAGE010
为机 组i风轮转速，
Figure 312293DEST_PATH_IMAGE042
为机组i发电机转速。
令阻尼转矩的最大值为调整后的给定发电机励磁转矩
Figure 59669DEST_PATH_IMAGE045
的百分之五，因此动态 分配附加功率下各机组的增益系数
Figure 738912DEST_PATH_IMAGE011
动态调整。
由于转矩扰动在传动链共振频率附近的幅值增益最高，因此传动链附加阻尼转矩通过二阶带通滤波器处理，滤波器的频率为的传动链固有振动频率。
3）为传动链设计动态相角补偿，补偿控制存在的滞后性。
由于在实际的控制系统中，控制回路经历信号采集、逻辑运算、数据传输、指令执行等过程，移相环节的拉普拉斯变换有效地补偿相角偏差。
实施例2：
基于同一发明构思，本发明还提供了一种双馈型机组风电场虚拟惯量控制系统。该系统结构如图2所示，包括：一次调频附加功率模块、功率分配模块和控制模块；
其中，一次调频附加功率模块，用于基于电网系统的实际频率计算双馈型机组风电场的一次调频附加功率；
所述功率分配模块，用于基于各机组的轮毂中心风速和风轮转速分配所述一次调频附加功率，得到各机组的惯性额外功率和下垂额外功率；
所述控制模块，用于根据各机组的惯性额外功率、下垂额外功率和发电机转速，控制各机组在匹配惯量支撑能力前提下协同参与调频。
其中，控制模块具体用于：
根据各机组的惯性额外功率、下垂额外功率和发电机转速，分别调整各机组的发电机励磁转矩；
根据各机组的风轮转速、发电机转速和发电机励磁转矩，分别为各机组的传动链施加带通滤波阻尼；
基于各机组的发电机励磁转矩和带通滤波阻尼，控制各机组在匹配惯量支撑能力前提下协同参与调频。
其中，发电机励磁转矩的计算式如下：
Figure 888134DEST_PATH_IMAGE002
式中，
Figure 627420DEST_PATH_IMAGE003
为机组i的发电机励磁转矩，
Figure 95048DEST_PATH_IMAGE004
为机组i的下垂额外功率，
Figure 312402DEST_PATH_IMAGE005
为机组i 的惯性额外功率，
Figure 581710DEST_PATH_IMAGE006
为机组i当前初始功率，
Figure 491897DEST_PATH_IMAGE007
为机组i的发电机转速。
其中，带通滤波阻尼的计算式如下：
Figure 213865DEST_PATH_IMAGE008
式中，
Figure 234911DEST_PATH_IMAGE009
为机组i的带通滤波阻尼，
Figure 594610DEST_PATH_IMAGE010
为机组i的风轮转速，
Figure 675699DEST_PATH_IMAGE007
为机组i的发电机 转速，
Figure 884963DEST_PATH_IMAGE011
为机组i的增益系数；
所述机组i的增益系数
Figure 709700DEST_PATH_IMAGE011
的值根据机组i的发电机励磁转矩动态调整。
其中，根据各机组的风轮转速和发电机转速，分别为各机组的传动链施加带通滤波阻尼之后，且基于各机组的发电机励磁转矩和带通滤波阻尼，控制各机组在匹配惯量支撑能力前提下协同参与调频之前，还包括：
分别为各机组的传动链设计动态相角补偿，补偿控制存在的滞后性。
其中，一次调频附加功率的计算式如下：
Figure 688020DEST_PATH_IMAGE012
式中，
Figure 172966DEST_PATH_IMAGE013
为一次调频附加功率，K为风电场一次调频比例系数，
Figure 869526DEST_PATH_IMAGE014
为电网系统的 参考频率，
Figure 966795DEST_PATH_IMAGE015
为电网系统的实际频率，
Figure 65201DEST_PATH_IMAGE016
为风电场频率监测死区。
其中，功率分配模块具体用于：
基于各机组的轮毂中心风速计算各机组的惯性系数的比例系数以及下垂系数的比例系数；
基于各机组的惯性系数的比例系数、下垂系数的比例系数以及风轮转速，将所述一次调频附加功率分配至参与频率响应控制的各机组，得到各机组的惯性额外功率和下垂额外功率。
其中，各机组的惯性系数的比例系数的计算式如下：
Figure 222513DEST_PATH_IMAGE017
式中，
Figure 406370DEST_PATH_IMAGE018
为机组i的惯性系数的比例系数，v i 为机组i的轮毂处实时风速，
Figure 74374DEST_PATH_IMAGE019
为机 组的切入风速，
Figure 761707DEST_PATH_IMAGE020
为机组的切出风速。
其中，各机组的下垂系数的比例系数的计算式如下：
Figure 355500DEST_PATH_IMAGE021
式中，
Figure 761073DEST_PATH_IMAGE022
为机组i的下垂系数的比例系数，a为预设常数，v i 为机组i的轮毂处实时 风速，
Figure 465724DEST_PATH_IMAGE019
为机组的切入风速，
Figure 771678DEST_PATH_IMAGE020
为机组的切出风速。
其中，基于各机组的惯性系数的比例系数、下垂系数的比例系数以及风轮转速，将所述一次调频附加功率分配至参与频率响应控制的各机组，得到各机组的惯性额外功率和下垂额外功率的计算式如下：
Figure 536372DEST_PATH_IMAGE023
Figure 429241DEST_PATH_IMAGE024
Figure 937583DEST_PATH_IMAGE025
式中，
Figure 599509DEST_PATH_IMAGE013
为一次调频附加功率，
Figure 269524DEST_PATH_IMAGE026
为参与频率响应控制机组的数量，
Figure 151155DEST_PATH_IMAGE005
为机组i的惯性额外功率，
Figure 197608DEST_PATH_IMAGE004
为机组i的下垂额外功率，
Figure 979620DEST_PATH_IMAGE018
为机组i的惯性系数的比例系数，
Figure 86116DEST_PATH_IMAGE027
为 机组惯性系数，
Figure 484736DEST_PATH_IMAGE022
为机组i的下垂系数的比例系数，
Figure 856853DEST_PATH_IMAGE028
为机组下垂系数，
Figure 493371DEST_PATH_IMAGE029
为机组i的风轮 转速，
Figure 770769DEST_PATH_IMAGE030
为惯量响应机组的风轮转速下限，
Figure 859947DEST_PATH_IMAGE031
为惯量响应机组风轮转速上限，
Figure 513783DEST_PATH_IMAGE032
为机 组频率与电网系统的实际频率的差。
实施例3：
本发明还提供一种计算机设备，该设备包括：一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行时，实现如实施例1中所述的方法，或发明内容中所述的方法。
实施例4：
本发明还提供一种计算机可读存储介质，其上存有计算机程序，所述计算机程序被执行时，实现如实施例1或发明内容中所述的方法。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质（包括但不限于磁盘存储器、CD-ROM、光学存储器等）上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备（系统）、和计算机程序产品的流程图和／或方框图来描述的。应理解可由计算机程序指令实现流程图和／或方框图中的每一流程和／或方框、以及流程图和／或方框图中的流程和／或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和／或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用于说明本发明的技术方案而非对其保护范围的限制，尽管参照上述实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解:本领域技术人员阅读本发明后依然可对申请的具体实施方式进行种种变更、修改或者等同替换，但这些变更、修改或者等同替换，均在申请待批的权利要求保护范围之内。

Claims (7)

1.一种双馈型机组风电场虚拟惯量控制方法，其特征在于，包括：
基于电网系统的实际频率计算双馈型机组风电场的一次调频附加功率；
基于各机组的轮毂中心风速和风轮转速分配所述一次调频附加功率，得到各机组的惯性额外功率和下垂额外功率；
根据各机组的惯性额外功率、下垂额外功率和发电机转速，控制各机组在匹配惯量支撑能力前提下协同参与调频；
所述一次调频附加功率的计算式如下：
Figure 125070DEST_PATH_IMAGE002
式中，
Figure DEST_PATH_IMAGE003
为一次调频附加功率，K为风电场一次调频比例系数，
Figure 207427DEST_PATH_IMAGE004
为电网系统的参考频率，
Figure DEST_PATH_IMAGE005
为电网系统的实际频率，
Figure 432128DEST_PATH_IMAGE006
为风电场频率监测死区；
所述基于各机组的轮毂中心风速和风轮转速分配所述一次调频附加功率，得到各机组的惯性额外功率和下垂额外功率，包括：
基于各机组的轮毂中心风速计算各机组的惯性系数的比例系数以及下垂系数的比例系数；
基于各机组的惯性系数的比例系数、下垂系数的比例系数以及风轮转速，将所述一次调频附加功率分配至参与频率响应控制的各机组，得到各机组的惯性额外功率和下垂额外功率；
所述各机组的惯性系数的比例系数的计算式如下：
Figure 419676DEST_PATH_IMAGE008
式中，
Figure DEST_PATH_IMAGE009
为机组i的惯性系数的比例系数，v i 为机组i的轮毂处实时风速，
Figure 253771DEST_PATH_IMAGE010
为机组的切入风速，
Figure DEST_PATH_IMAGE011
为机组的切出风速；
所述各机组的下垂系数的比例系数的计算式如下：
Figure DEST_PATH_IMAGE013
式中，
Figure 485907DEST_PATH_IMAGE014
为机组i的下垂系数的比例系数，a为预设常数；
所述基于各机组的惯性系数的比例系数、下垂系数的比例系数以及风轮转速，将所述一次调频附加功率分配至参与频率响应控制的各机组，得到各机组的惯性额外功率和下垂额外功率的计算式如下：
Figure 583307DEST_PATH_IMAGE016
Figure 58151DEST_PATH_IMAGE018
Figure 400664DEST_PATH_IMAGE020
式中，
Figure DEST_PATH_IMAGE021
为参与频率响应控制机组的数量，
Figure 192033DEST_PATH_IMAGE022
为机组i的惯性额外功率，
Figure DEST_PATH_IMAGE023
为机组i的下垂额外功率，
Figure 975181DEST_PATH_IMAGE024
为机组惯性系数，
Figure DEST_PATH_IMAGE025
为机组下垂系数，
Figure 983326DEST_PATH_IMAGE026
为机组i的风轮转速，
Figure DEST_PATH_IMAGE027
为惯量响应机组的风轮转速下限，
Figure 893645DEST_PATH_IMAGE028
为惯量响应机组风轮转速上限，
Figure DEST_PATH_IMAGE029
为机组频率与电网系统的实际频率的差。
2.如权利要求1所述的方法，其特征在于，所述根据各机组的惯性额外功率、下垂额外功率和发电机转速，控制各机组在匹配惯量支撑能力前提下协同参与调频，包括：
根据各机组的惯性额外功率、下垂额外功率和发电机转速，分别调整各机组的发电机励磁转矩；
根据各机组的风轮转速、发电机转速和发电机励磁转矩，分别为各机组的传动链施加带通滤波阻尼；
基于各机组的发电机励磁转矩和带通滤波阻尼，控制各机组在匹配惯量支撑能力前提下协同参与调频；
所述发电机励磁转矩的计算式如下：
Figure DEST_PATH_IMAGE031
式中，
Figure 91584DEST_PATH_IMAGE032
为机组i的发电机励磁转矩，
Figure DEST_PATH_IMAGE033
为机组i当前初始功率，
Figure 796366DEST_PATH_IMAGE034
为机组i的发电机转速；
所述带通滤波阻尼的计算式如下：
Figure 495069DEST_PATH_IMAGE036
式中，
Figure DEST_PATH_IMAGE037
为机组i的带通滤波阻尼，
Figure 723925DEST_PATH_IMAGE038
为机组i的风轮转速；
所述机组i的增益系数
Figure DEST_PATH_IMAGE039
的值根据机组i的发电机励磁转矩动态调整。
3.如权利要求2所述的方法，其特征在于，所述根据各机组的风轮转速和发电机转速，分别为各机组的传动链施加带通滤波阻尼之后，且基于各机组的发电机励磁转矩和带通滤波阻尼，控制各机组在匹配惯量支撑能力前提下协同参与调频之前，还包括：
分别为各机组的传动链设计动态相角补偿，补偿控制存在的滞后性。
4.一种双馈型机组风电场虚拟惯量控制系统，其特征在于，包括：一次调频附加功率模块、功率分配模块和控制模块；
所述一次调频附加功率模块，用于基于电网系统的实际频率计算双馈型机组风电场的一次调频附加功率；
所述功率分配模块，用于基于各机组的轮毂中心风速和风轮转速分配所述一次调频附加功率，得到各机组的惯性额外功率和下垂额外功率；
所述控制模块，用于根据各机组的惯性额外功率、下垂额外功率和发电机转速，控制各机组在匹配惯量支撑能力前提下协同参与调频；
所述一次调频附加功率的计算式如下：
Figure 693150DEST_PATH_IMAGE002
式中，
Figure 286942DEST_PATH_IMAGE003
为一次调频附加功率，K为风电场一次调频比例系数，
Figure 210292DEST_PATH_IMAGE004
为电网系统的参考频率，
Figure 180522DEST_PATH_IMAGE005
为电网系统的实际频率，
Figure 738673DEST_PATH_IMAGE006
为风电场频率监测死区；
所述功率分配模块具体用于：
基于各机组的轮毂中心风速计算各机组的惯性系数的比例系数以及下垂系数的比例系数；
基于各机组的惯性系数的比例系数、下垂系数的比例系数以及风轮转速，将所述一次调频附加功率分配至参与频率响应控制的各机组，得到各机组的惯性额外功率和下垂额外功率；
所述各机组的惯性系数的比例系数的计算式如下：
Figure 768946DEST_PATH_IMAGE008
式中，
Figure 911084DEST_PATH_IMAGE009
为机组i的惯性系数的比例系数，v i 为机组i的轮毂处实时风速，
Figure 419425DEST_PATH_IMAGE010
为机组的切入风速，
Figure 346930DEST_PATH_IMAGE011
为机组的切出风速；
所述各机组的下垂系数的比例系数的计算式如下：
Figure 33257DEST_PATH_IMAGE040
式中，
Figure 413423DEST_PATH_IMAGE014
为机组i的下垂系数的比例系数，a为预设常数；
所述基于各机组的惯性系数的比例系数、下垂系数的比例系数以及风轮转速，将所述一次调频附加功率分配至参与频率响应控制的各机组，得到各机组的惯性额外功率和下垂额外功率的计算式如下：
Figure 977653DEST_PATH_IMAGE016
Figure DEST_PATH_IMAGE041
Figure 572714DEST_PATH_IMAGE020
式中，
Figure 679210DEST_PATH_IMAGE021
为参与频率响应控制机组的数量，
Figure 546672DEST_PATH_IMAGE022
为机组i的惯性额外功率，
Figure 911663DEST_PATH_IMAGE023
为机组i的下垂额外功率，
Figure 548181DEST_PATH_IMAGE024
为机组惯性系数，
Figure 841890DEST_PATH_IMAGE025
为机组下垂系数，
Figure 931068DEST_PATH_IMAGE026
为机组i的风轮转速，
Figure 584904DEST_PATH_IMAGE027
为惯量响应机组的风轮转速下限，
Figure 593704DEST_PATH_IMAGE028
为惯量响应机组风轮转速上限，
Figure 42003DEST_PATH_IMAGE029
为机组频率与电网系统的实际频率的差。
5.如权利要求4所述的系统，其特征在于，所述控制模块具体用于：
根据各机组的惯性额外功率、下垂额外功率和发电机转速，分别调整各机组的发电机励磁转矩；
根据各机组的风轮转速、发电机转速和发电机励磁转矩，分别为各机组的传动链施加带通滤波阻尼；
基于各机组的发电机励磁转矩和带通滤波阻尼，控制各机组在匹配惯量支撑能力前提下协同参与调频；
所述发电机励磁转矩的计算式如下：
Figure 900369DEST_PATH_IMAGE031
式中，
Figure 92316DEST_PATH_IMAGE032
为机组i的发电机励磁转矩，
Figure 687114DEST_PATH_IMAGE023
为机组i的下垂额外功率，
Figure 837473DEST_PATH_IMAGE022
为机组i的惯性额外功率，
Figure 901244DEST_PATH_IMAGE033
为机组i当前初始功率，
Figure 913193DEST_PATH_IMAGE034
为机组i的发电机转速；
所述带通滤波阻尼的计算式如下：
Figure 113230DEST_PATH_IMAGE036
式中，
Figure 432827DEST_PATH_IMAGE037
为机组i的带通滤波阻尼，
Figure 983894DEST_PATH_IMAGE038
为机组i的风轮转速，
Figure 517644DEST_PATH_IMAGE034
为机组i的发电机转速，
Figure 854078DEST_PATH_IMAGE039
为机组i的增益系数；
所述机组i的增益系数
Figure 80660DEST_PATH_IMAGE039
的值根据机组i的发电机励磁转矩动态调整。
6.一种计算机设备，其特征在于，包括：一个或多个处理器；
存储器，用于存储一个或多个程序；
当所述一个或多个程序被所述一个或多个处理器执行时，实现如权利要求1至3中任一项所述的方法。
7.一种计算机可读存储介质，其特征在于，其上存有计算机程序，所述计算机程序被执行时，实现如权利要求1至3中任一项所述的方法。
CN202111323003.5A 2021-11-10 2021-11-10 一种双馈型机组风电场虚拟惯量控制方法和系统 Active CN113783208B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111323003.5A CN113783208B (zh) 2021-11-10 2021-11-10 一种双馈型机组风电场虚拟惯量控制方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111323003.5A CN113783208B (zh) 2021-11-10 2021-11-10 一种双馈型机组风电场虚拟惯量控制方法和系统

Publications (2)

Publication Number Publication Date
CN113783208A CN113783208A (zh) 2021-12-10
CN113783208B true CN113783208B (zh) 2022-03-08

Family

ID=78956976

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111323003.5A Active CN113783208B (zh) 2021-11-10 2021-11-10 一种双馈型机组风电场虚拟惯量控制方法和系统

Country Status (1)

Country Link
CN (1) CN113783208B (zh)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109936150B (zh) * 2017-12-15 2024-02-02 国家电网公司 一种虚拟惯量控制的优化控制方法、装置及其控制器
PT3644497T (pt) * 2018-10-24 2022-05-26 Ge Energy Power Conversion Technology Ltd Procedimento e dispositivo para comando de inércia virtual, destinados a centrais elétricas com máquina assíncrona de dupla alimentação
CN112502900B (zh) * 2020-10-16 2021-10-15 浙江工业大学 基于非线性阻尼控制的风电齿轮箱瞬变载荷主动抑制方法
CN113036779B (zh) * 2021-03-03 2022-06-14 山东大学 一种风电参与电网调频的时变调频参数设定方法及系统
CN113224774B (zh) * 2021-05-17 2022-07-26 山东大学 一种电网频率扰动过程风储协调控制方法及系统

Also Published As

Publication number Publication date
CN113783208A (zh) 2021-12-10

Similar Documents

Publication Publication Date Title
CN104917201B (zh) 模拟惯性与超速相结合的双馈风机有功频率控制器及方法
Zhang et al. A frequency regulation strategy for wind power based on limited over-speed de-loading curve partitioning
CN105134485B (zh) 一种双馈风力机组惯性调频主动转速保护控制系统与方法
CN109936152B (zh) 高渗透率风电并网后的电网调频控制方法、风力发电机
CN109861251B (zh) 一种用于微网暂稳态频率优化的双馈风机综合控制方法
CN107910896A (zh) 基于虚拟惯性和桨距角控制风电机组参与电网调频的方法
CN107895955B (zh) 一种风电补偿水轮机水锤效应的协同控制方法
Zhang et al. Comparison of inertia control methods for DFIG-based wind turbines
CN108242814A (zh) 变速变频风电机组的一次调频方法和设备
CN107681689A (zh) 双馈风机在微电网中的频率控制参数选取方法
CN106532739A (zh) 风电机组分频段参与电力系统一次调频方法
CN107069799B (zh) 双馈式风力发电机组的控制方法和系统
CN113054672A (zh) 一种考虑多环节协同控制的风电调频方法及系统
CN112994043A (zh) 自同步双馈风电机组惯量与一次调频的控制方法及系统
CN112332442A (zh) 基于双馈风力发电场虚拟惯量控制策略优化
Jeon et al. PI control loop–based frequency smoothing of a doubly-fed induction generator
CN113783208B (zh) 一种双馈型机组风电场虚拟惯量控制方法和系统
CN113224774A (zh) 一种电网频率扰动过程风储协调控制方法及系统
CN111211582B (zh) 一种基于风电场中双馈风机调频能力的电网综合调频方法
CN105048511A (zh) 一种含可控惯性风力发电机的发电系统惯性综合控制方法
CN112271760A (zh) 适用于直驱风电场交流并网的调频控制方法
CN117060484A (zh) 一种基于风储联合调频系统的改进自适应控制方法
CN117117901A (zh) 一种海上风电柔直系统的频率控制方法
CN110460098A (zh) 基于虚拟质量块的风力机双质块轴系稳定控制方法
CN106952180B (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