CN113406888A - 基于ladrc的汽轮发电机控制系统及其静态补偿系数确定方法 - Google Patents

基于ladrc的汽轮发电机控制系统及其静态补偿系数确定方法 Download PDF

Info

Publication number
CN113406888A
CN113406888A CN202110945603.9A CN202110945603A CN113406888A CN 113406888 A CN113406888 A CN 113406888A CN 202110945603 A CN202110945603 A CN 202110945603A CN 113406888 A CN113406888 A CN 113406888A
Authority
CN
China
Prior art keywords
input signal
steam turbine
path
output
control system
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.)
Pending
Application number
CN202110945603.9A
Other languages
English (en)
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.)
Henan Institute of Technology
Original Assignee
Henan Institute of Technology
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 Henan Institute of Technology filed Critical Henan Institute of Technology
Priority to CN202110945603.9A priority Critical patent/CN113406888A/zh
Publication of CN113406888A publication Critical patent/CN113406888A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
    • G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
    • G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
    • G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance

Landscapes

  • Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Evolutionary Computation (AREA)
  • Medical Informatics (AREA)
  • Software Systems (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Control Of Turbines (AREA)

Abstract

本发明提供一种基于LADRC的汽轮发电机控制系统及其静态补偿系数确定方法。该控制系统具有调速器死区或汽轮机发电速率约束，包括调速器、汽轮机、电力系统、线性扩张状态观测器，所述控制系统还包括：误差补偿模块；所述误差补偿模块，针对调速器死区，用于将调速器的理论输出与实际输出之间的误差反馈至线性扩张状态观测器，利用线性扩张状态观测器进行估计，使LADRC达到消除调速器死区的目的；针对汽轮机发电速率约束，将汽轮机的理论输出与实际输出的误差作为外部扰动，加入LESO进行估计，从而使LADRC消除GRC的影响。本发明提出的控制系统具有较好的控制性能，采用描述函数法获得补偿系数的取值范围是可行且有效。

Description

基于LADRC的汽轮发电机控制系统及其静态补偿系数确定 方法
技术领域
本发明涉及汽轮发电机技术领域，尤其涉及一种基于LADRC的汽轮发电机控制系统及其静态补偿系数确定方法。
背景技术
随着现代电力系统规模和复杂性的不断增加，由系统振荡引起的广域停电的风险也在增加。因此，已有很多先进的控制方法被用于解决负荷频率控制问题，如惯性控制和下降速率控制策略、自适应模糊增益调度、滑模控制（sliding mode control，SMC）方法和2-自由度（two-degree-of-freedom，TDF）内模控制（internal model control，IMC）的PID（proportional-integral-derivative）整定方法（参考文献Tan W. Unified Tuning ofPID load frequency controller for power systems via IMC[J]. IEEE Transactionson Power Systems, 2010, V25(1):341-350）。但是，这些文献忽略了系统存在非线性特性的情况。而在汽轮发电机控制系统中，存在多种典型非线性特性，例如：调速器死区和发电速率约束(generation rate constraint，GRC)，它们会恶化系统的动态性能，甚至不稳定。目前也已有文献针对这些问题进行了研究。如模糊C均值聚类技术(fuzzy C-meansclustering technique，FCM)、细菌觅食优化算法(bacteria foraging optimizationalgorithm，BFOA)、万有引力搜索算法(gravitational search algorithm，GSA)和模式搜索法(pattern search，PS)和多变量模型预测控制(model based predictive control，MPC)方法（参考文献Mojtaba S, Mohammad R, Ali M. Robust multivariablepredictive based load frequency control considering generation rateconstraint[J]. Electrical Power and Energy Systems, 2013, 46: 405–413）。上述这些文献针对汽轮发电机控制系统存在的非线性问题，提出了多种控制策略给予解决，并进行参数优化整定。不过，其所设计的控制器结构相对复杂，计算量大，工业应用困难。
发明内容
针对目前的控制器结构相对复杂、计算量大的问题，本发明提供一种基于LADRC的汽轮发电机控制系统及其静态补偿系数确定方法。
第一方面，本发明提供一种基于 LADRC的汽轮发电机控制系统，所述控制系统具有调速器死区，包括调速器、汽轮机、电力系统和线性扩张状态观测器，所述控制系统还包括：误差补偿模块；
所述误差补偿模块，用于将调速器的理论输出与实际输出之间的误差反馈至线性扩张状态观测器，利用线性扩张状态观测器进行估计，使LADRC达到消除调速器死区的目的。
进一步地，所述控制系统还包括：状态反馈增益模块；所述误差补偿模块包括上、下两个补偿分支；其中：
控制系统的输出为频率变化Δf；将控制系统的输出作为线性扩张状态观测器的第一路输入信号；
线性扩张状态观测器的输出与原始输入信号u作减法运算后作为状态反馈增益模块的输入；
状态反馈增益模块的输出分两路：一路与控制系统的输出Δf经1/R倍增益后的结果Δf/R作减法运算后输出中间量u’；另一路作为误差补偿模块的下补偿分支的第一路输入信号；R为机组下降特性；
将所述中间量u’直接作为误差补偿模块的上补偿分支的第一路输入信号，将调速器的实际输出直接作为误差补偿模块的上补偿分支的第二路输入信号；
误差补偿模块的上补偿分支的第一路输入信号和第二路输入信号作减法运算后再经过k倍增益后作为误差补偿模块的下补偿分支的第二路输入信号；k为可手动调节的静态补偿系数；
误差补偿模块的下补偿分支的第一路输入信号和第二路输入信号作减法运算后作为线性扩张状态观测器的第二路输入信号。
进一步地，所述中间量u’经过调速器死区后还作为计算调速器的输入信号的第一参考量；调速器的实际输出还分别直接作为计算调速器的输入信号的第二参考量以及直接作为计算计算汽轮机的输入信号的第一参考量；其中，调速器的输入信号的第一参考量与第二参考量作减法运算后作为调速器的输入信号；
汽轮机的输出分为两路，一路直接作为计算汽轮机的输入信号的第二参考量，另一路直接作为计算电力系统的输入信号的第一参考量；其中，汽轮机的输入信号的第一参考量与第二参考量作减法运算后作为汽轮机的输入信号；
电力系统的输出分为两路，一路直接作为计算电力系统的输入信号的第二参考量，另一路经过K p 倍增益后即为控制系统的输出Δf；将负荷扰动ΔP d 作为计算电力系统的输入信号的第三参考量；其中，电力系统的输入信号的第一参考量依次与第二参考量和第三参考量作减法运算后作为电力系统的输入信号；K p 为发电机增益。
第二方面，本发明还提供一种如上述的基于 LADRC的汽轮发电机控制系统的静态补偿系数确定方法，包括：
步骤1：获得非线性元件的描述函数N(X)，并绘出负倒描述函数-1/N(X)的曲线，具体包括：具有死区非线性特性的描述函数N(X)如公式（1），其中，死区非线性的死区宽度和线性输出特性的斜率分别为0.1和1：
Figure 628472DEST_PATH_IMAGE001
(1)
且有：
Figure 109132DEST_PATH_IMAGE002
其中，X表示非线性的振幅，K表示线性输出特性的斜率，Δ表示死区非线性的死区宽度；
步骤2：求传递函数G(jq)，具体为：定义
Figure 76607DEST_PATH_IMAGE003
其中，G L 是LADRC的传递函数；jq表示频率特性，g 1表示线性环节的总传递函数，g 2表示调速器的传递函数，g 3表示汽轮机的传递函数，g 4表示电力系统的传递函数，g 5表示机组下降特性的倒数，g 6表示静态补偿系数k，T G 表示调速器时间常数，T T 表示汽轮机时间常数， T P 表示发电机时间常数，s表示微分算子；
然后通过信号流图和梅逊公式可以获得如公式（2）所示的传递函数G(jq)：
Figure 668125DEST_PATH_IMAGE004
(2)
步骤3：绘出传递函数G(jq)的奈奎斯特曲线，然后应用奈奎斯特稳定判据进行分析，具体为：将G(jq)的奈奎斯特曲线和-1/N(X)的曲线绘制在同一个坐标系中，G(jq)的奈奎斯特曲线与负实轴的交点为
Figure 362412DEST_PATH_IMAGE005
，若G(jq)的奈奎斯特曲线不包围-1/N(X)的曲线，即
Figure 482814DEST_PATH_IMAGE006
，则该系统稳定；若G(jq)的奈奎斯特曲线包围-1/N(X)的曲线，即
Figure 528131DEST_PATH_IMAGE007
，则该系统不稳定。
第三方面，本发明提供一种基于LADRC的汽轮发电机控制系统，所述控制系统具有汽轮机发电速率约束，包括调速器、汽轮机、电力系统和线性扩张状态观测器，所述控制系统还包括：误差补偿模块；
所述误差补偿模块，用于将汽轮机的理论输出与实际输出之间的误差作为外部扰动反馈至线性扩张状态观测器，利用线性扩张状态观测器进行估计，使LADRC消除汽轮机发电速率约束的影响。
进一步地，所述控制系统还包括：状态反馈增益模块；所述误差补偿模块包括上、下两个补偿分支；其中：
控制系统的输出为频率变化Δf；将控制系统的输出作为线性扩张状态观测器的第一路输入信号；
线性扩张状态观测器的输出与原始输入信号u作减法运算后作为状态反馈增益模块的输入；
状态反馈增益模块的输出分两路：一路与控制系统的输出Δf经1/R 倍增益后的结果Δf/R作减法运算后输出中间量u’；另一路作为误差补偿模块的下补偿分支的第一路输入信号；R为机组下降特性；
将所述中间量u’依次经过
Figure 595444DEST_PATH_IMAGE008
和
Figure 144237DEST_PATH_IMAGE009
倍增益后作为误差补偿模块的上补偿分支的第一路输入信号，将汽轮机的实际输出直接作为误差补偿模块的上补偿分支的第二路输入信号；其中，T G 表示调速器时间常数，T T 表示汽轮机时间常数；
误差补偿模块的上补偿分支的第一路输入信号和第二路输入信号作减法运算后再经过k 1倍增益后作为误差补偿模块的下补偿分支的第二路输入信号；k 1为可手动调节的静态补偿系数；
误差补偿模块的下补偿分支的第一路输入信号和第二路输入信号作减法运算后作为线性扩张状态观测器的第二路输入信号。
进一步地，所述中间量u’还直接作为计算调速器的输入信号的第一参考量；
调速器的输出分为两路，一路直接作为计算调速器的输入信号的第二参考量，另一路直接作为计算汽轮机的输入信号的第一参考量；其中，调速器的输入信号的第一参考量与第二参考量作减法运算后作为调速器的输入信号；
汽轮机的实际输出还分别直接作为计算汽轮机的输入信号的第二参考量以及直接作为计算电力系统的输入信号的第一参考量；其中，汽轮机的输入信号的第一参考量与第二参考量作减法运算后作为汽轮机的输入信号；
电力系统的输出分为两路，一路直接作为计算电力系统的输入信号的第二参考量，另一路经过K p 倍增益后即为控制系统的输出Δf；将负荷扰动ΔP d 作为计算电力系统的输入信号的第三参考量；其中，电力系统的输入信号的第一参考量依次与第二参考量和第三参考量作减法运算后作为电力系统的输入信号；K p 为发电机增益。
第四方面，本发明还提供一种如上述的基于LADRC的汽轮发电机控制系统的静态补偿系统确定方法，包括：
步骤1：获得非线性元件的描述函数N 1(X)，并绘出负倒描述函数-1/N 1(X)的曲线，具体包括；具有饱和非线性特性的描述函数N 1(X)如公式（3），其中，饱和非线性的线性区宽度和线性区特性的斜率分别为0.0017和1：
Figure 701120DEST_PATH_IMAGE010
(3)
且有：
Figure 968154DEST_PATH_IMAGE011
其中，X表示非线性的振幅，K表示线性区特性的斜率，a表示饱和非线性的线性区宽度；
步骤2：求传递函数G 1(jq)，具体为：定义
Figure 901474DEST_PATH_IMAGE012
其中，
Figure 304774DEST_PATH_IMAGE013
是LADRC的传递函数；jq表示频率特性，
Figure 32559DEST_PATH_IMAGE014
表示线性环节的总传递函数，
Figure 786888DEST_PATH_IMAGE015
表示调速器的传递函数，
Figure 196004DEST_PATH_IMAGE016
和
Figure 719389DEST_PATH_IMAGE017
均表示汽轮机的传递函数，
Figure 618075DEST_PATH_IMAGE018
表示电力系统的传递函 数，
Figure 594121DEST_PATH_IMAGE019
表示电力系统的传递函数的倒数，
Figure 134824DEST_PATH_IMAGE020
表示调速器的传递函数与汽轮机的传递函数 的乘积，
Figure 512716DEST_PATH_IMAGE021
表示机组下降特性的倒数，
Figure 847882DEST_PATH_IMAGE022
表示静态补偿系数k 1，
Figure 248907DEST_PATH_IMAGE023
表示常数值1，T G 表示调 速器时间常数，T T 表示汽轮机时间常数， T P 表示发电机时间常数，s表示微分算子；
然后通过信号流图和梅逊公式可以获得如公式（4）所示的传递函数G 1(jq)：
Figure 327722DEST_PATH_IMAGE024
(4)
步骤3：绘出传递函数G 1(jq)的奈奎斯特曲线,然后应用奈奎斯特稳定判据进行分析，具体为：将G 1(jq)的奈奎斯特曲线和-1/N 1(X)的曲线绘制在同一个坐标系中，G 1(jq)的奈奎斯特曲线与负实轴的交点为
Figure 294541DEST_PATH_IMAGE025
，若G 1(jq)的奈奎斯特曲线不包围-1/N 1(X)的曲线，即
Figure 800609DEST_PATH_IMAGE026
，则该系统稳定；若G 1(jq)的奈奎斯特曲线包围-1/N 1(X)的曲线，即
Figure 751247DEST_PATH_IMAGE027
，则该系统不稳定。
本发明的有益效果：
相对于其他控制算法来说，线性自抗扰控制算法（linear active disturbancerejection control, LADRC）是一种独立于被控对象模型的普适性控制结构，结构更简单且适用于各种非线性系统，特别地，线性自抗扰控制算法只需整定2个参数，因此很容易为控制工程师所理解。所以，针对存在调速器死区或汽轮机GRC的汽轮发电机控制系统，本发明采用了线性自抗扰控制算法。
然而，线性自抗扰控制算法作为一种一般控制结构，当系统存在非线性特性时，需要将LADRC转化为高阶控制器，但这样会破坏其固有的控制结构，可能会使其失去主动抗扰的特性。因此，针对存在调速器死区或汽轮机GRC的汽轮发电机控制系统，本发明提出了采用基于线性自抗扰控制算法的误差补偿策略的新的汽轮发电机控制系统，使线性自抗扰控制算法能够快速地消除调速器死区或汽轮机GRC的影响。另外，该补偿策略有一个可手动调节的静态补偿系数k或k 1，但是，随着该补偿系数的变大，系统的控制性能也会受到影响，故需要确认该补偿系数的取值范围。所以，对应地还分别提出了采用描述函数法进行稳定验证和获取该补偿策略中静态补偿系数的取值范围。仿真结果显示，本发明提出的误差补偿策略对改善系统的控制性能有良好的表现，且采用描述函数法获取补偿系数的取值范围，为该系数的选取提供参考也是可行的。
附图说明
图1为本发明实施例提供的一种基于 LADRC的汽轮发电机控制系统的结构示意图之一；
图2为本发明实施例提供的在采用误差补偿策略的LADRC的电力系统阶跃响应示意图之一；
图3为本发明实施例提供的G(jq)的奈奎斯特曲线和-1/N(X)的曲线示意图；
图4为本发明实施例提供的在采用误差补偿策略的LADRC的电力系统阶跃响应之二；
图5为本发明实施例提供的在采用误差补偿策略的LADRC的电力系统阶跃响应之三；
图6为本发明实施例提供的一种基于 LADRC的汽轮发电机控制系统的结构示意图之二；
图7为本发明实施例提供的在采用误差补偿策略的LADRC的电力系统阶跃响应示意图之四；
图8为本发明实施例提供的G 1(jq)的奈奎斯特曲线和-1/N 1(X)的曲线示意图；
图9为本发明实施例提供的在采用误差补偿策略的LADRC的电力系统阶跃响应示意图之五。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例1
在汽轮发电机控制系统中，调速器死区非线性是普遍存在的，因而会影响系统的控制性能。本发明实施例针对具有调速器死区非线性的汽轮发电机控制系统系统，将LADRC用到系统中。接着，在LADRC的基础上，提出了一种误差补偿策略，该补偿策略如图1所示。其思想是将调速器理论输出与实际输出的误差反馈回LESO，利用LESO进行估计，使LADRC达到消除死区的目的，从而使系统的控制性能快速地恢复和改善。
如图1所示，本发明实施例提供一种基于 LADRC的汽轮发电机控制系统，所述控制系统具有调速器死区，包括调速器、汽轮机、电力系统和线性扩张状态观测器，所述控制系统还包括：误差补偿模块；
所述误差补偿模块，用于将调速器的理论输出与实际输出之间的误差反馈至线性扩张状态观测器，利用线性扩张状态观测器进行估计，使LADRC达到消除调速器死区的目的，从而使控制系统的控制性能快速地恢复和改善。
作为一种可实施方式，该控制系统还包括状态反馈增益模块；所述误差补偿模块包括上、下两个补偿分支；其中：
控制系统的输出为频率变化Δf；将控制系统的输出作为线性扩张状态观测器的第一路输入信号；
线性扩张状态观测器的输出与原始输入信号u作减法运算后作为状态反馈增益模块的输入；
状态反馈增益模块的输出分两路：一路与控制系统的输出Δf经1/R倍增益后的结果Δf/R作减法运算后输出中间量u’；另一路作为误差补偿模块的下补偿分支的第一路输入信号；为机组下降特性；
将所述中间量u’直接作为误差补偿模块的上补偿分支的第一路输入信号，将调速器的实际输出直接作为误差补偿模块的上补偿分支的第二路输入信号；
误差补偿模块的上补偿分支的第一路输入信号和第二路输入信号作减法运算后再经过k倍增益后作为误差补偿模块的下补偿分支的第二路输入信号；k为可手动调节的静态补偿系数；
误差补偿模块的下补偿分支的第一路输入信号和第二路输入信号作减法运算后作为线性扩张状态观测器的第二路输入信号。
作为一种可实施方式，所述中间量u’经过调速器死区后还作为计算调速器的输入信号的第一参考量；调速器的实际输出还分别直接作为计算调速器的输入信号的第二参考量以及直接作为计算计算汽轮机的输入信号的第一参考量；其中，调速器的输入信号的第一参考量与第二参考量作减法运算后作为调速器的输入信号；
汽轮机的输出分为两路，一路直接作为计算汽轮机的输入信号的第二参考量，另一路直接作为计算电力系统的输入信号的第一参考量；其中，汽轮机的输入信号的第一参考量与第二参考量作减法运算后作为汽轮机的输入信号；
电力系统的输出分为两路，一路直接作为计算电力系统的输入信号的第二参考量，另一路经过K p 倍增益后即为控制系统的输出Δf；将负荷扰动ΔP d 作为计算电力系统的输入信号的第三参考量；其中，电力系统的输入信号的第一参考量依次与第二参考量和第三参考量作减法运算后作为电力系统的输入信号；K p 为发电机增益。
实施例2
上述实施例中的基于 LADRC的汽轮发电机控制系统中，存在一个可手动调节的静态补偿系数k，但是，随着该补偿系数的变大，系统的控制性能也会受到影响，故需要确认该补偿系数的取值范围。因此，对应上述实施例中的基于 LADRC的汽轮发电机控制系统，本发明实施例还提供一种基于 LADRC的汽轮发电机控制系统的静态补偿系数确定方法，本发明称之为描述函数法，该方法包括以下步骤：
S101：获得非线性元件的描述函数N(X)，并绘出负倒描述函数-1/N(X)的曲线，具体包括：具有死区非线性特性的描述函数N(X)如公式（1），其中，死区非线性的死区宽度和线性输出特性的斜率分别为0.1和1：
Figure 633752DEST_PATH_IMAGE001
(1)
且有：
Figure 452148DEST_PATH_IMAGE002
其中，X表示非线性的振幅，K表示线性输出特性的斜率，Δ表示死区非线性的死区宽度。
S102：求传递函数G(jq)，具体为：对于图1所示的单区域汽轮发电机控制系统，定义
Figure 394696DEST_PATH_IMAGE003
其中，G L 是LADRC的传递函数；jq表示频率特性，g 1表示线性环节的总传递函数，g 2表示调速器的传递函数，g 3表示汽轮机的传递函数，g 4表示电力系统的传递函数，g 5表示机组下降特性的倒数，g 6表示静态补偿系数k，T G 表示调速器时间常数，T T 表示汽轮机时间常数， T P 表示发电机时间常数，s表示微分算子。
然后通过信号流图和梅逊公式可以获得如公式（2）所示的传递函数G(jq)：
Figure 504735DEST_PATH_IMAGE004
(2)
S103：绘出传递函数G(jq)的奈奎斯特曲线，然后应用奈奎斯特稳定判据进行分析，具体为：可通过MATLAB将G(jq)的奈奎斯特曲线和-1/N(X)的曲线绘制在同一个坐标系中，G(jq)的奈奎斯特曲线与负实轴的交点为
Figure 925352DEST_PATH_IMAGE005
，若G(jq)的奈奎斯特曲线不包围-1/N(X)的曲线，即
Figure 866763DEST_PATH_IMAGE006
，则该系统稳定；若G(jq)的奈奎斯特曲线包围-1/N(X)的曲线，即
Figure 714633DEST_PATH_IMAGE007
，则该系统不稳定。
为了证实本发明提出的控制系统及其补偿系数确定方法的可行性与准确性，本发明考虑具有如下模型参数的单区域电力系统：
Figure 374285DEST_PATH_IMAGE028
采用LADRC进行控制，其参数选取为：
Figure 270697DEST_PATH_IMAGE029
其中，w c 是控制器带宽，w o 是观测器带宽，补偿系数取为k=0.7，为了显示LADRC在该补偿策略下的控制效果，在t=1秒时加入阶跃信号ΔP d =0.01，且死区取值为0.1，系统的响应曲线如图2所示。其中，黑色虚线表示当系统存在调速器死区时的情况，黑色实线表示系统不存在调速器死区的情况，黑色点线表示系统采用了误差补偿策略后的情况。可以看到，在不采取误差补偿策略时，系统由于死区非线性的存在，其控制性能受到恶化，而在采用误差补偿策略后，系统的控制性能得到了较好的改善。因此，该补偿策略在消除死区非线性方面是切实有效。
为了确定补偿系数的取值范围，采用描述函数法进行稳定性分析，要使系统稳定，则需满足
Figure 597773DEST_PATH_IMAGE006
，进而可得的取值范围0<k<0.91。而当k=0.91时，G(jq)的奈奎斯特曲线包围-1/N(X)的曲线，则系统不稳定，其结果如图3所示。接着，我们通过手动调节补偿系数k，来确认其的取值范围，仿真结果如图4、5所示。显而易见地，当k逐渐变大时，补偿效果逐渐提高，但是增加到0.91时，系统变得不稳定。仿真结果证实，采用描述函数法对具有调速器死区的汽轮发电机控制系统进行稳定性分析，获取静态补偿系数k的取值范围是可行的。
实施例3
在汽轮发电机控制系统中，汽轮机GRC会影响系统的控制性能，出现不稳定。因此，本发明实施例在LADRC的基础上，提出了一种误差补偿策略，其补偿策略如图6所示。其思想是将汽轮机的理论输出与实际输出的误差作为外部扰动，加入LESO进行估计，从而使LADRC消除GRC的影响，使系统的控制性能可以快速地恢复和改善。
如图6所示，本发明实施例提供一种基于LADRC的汽轮发电机控制系统，所述控制系统具有汽轮机发电速率约束，包括调速器、汽轮机、电力系统和线性扩张状态观测器，所述控制系统还包括：误差补偿模块；
所述误差补偿模块，用于将汽轮机的理论输出与实际输出之间的误差作为外部扰动反馈至线性扩张状态观测器，利用线性扩张状态观测器进行估计，使LADRC消除汽轮机发电速率约束的影响。
作为一种可实施方式，该控制系统还包括：状态反馈增益模块；所述误差补偿模块包括上、下两个补偿分支；其中：
控制系统的输出为频率变化Δf；将控制系统的输出作为线性扩张状态观测器的第一路输入信号；
线性扩张状态观测器的输出与原始输入信号u作减法运算后作为状态反馈增益模块的输入；
状态反馈增益模块的输出分两路：一路与控制系统的输出Δf经1/R 倍增益后的结果Δf/R作减法运算后输出中间量u’；另一路作为误差补偿模块的下补偿分支的第一路输入信号；R为机组下降特性；
将所述中间量u’依次经过
Figure 882124DEST_PATH_IMAGE008
和
Figure 763492DEST_PATH_IMAGE009
倍增益后作为误差补偿模块的上补偿分支的第一路输入信号，将汽轮机的实际输出直接作为误差补偿模块的上补偿分支的第二路输入信号；其中，T G 表示调速器时间常数，T T 表示汽轮机时间常数；
误差补偿模块的上补偿分支的第一路输入信号和第二路输入信号作减法运算后再经过k 1倍增益后作为误差补偿模块的下补偿分支的第二路输入信号；k 1为可手动调节的静态补偿系数；
误差补偿模块的下补偿分支的第一路输入信号和第二路输入信号作减法运算后作为线性扩张状态观测器的第二路输入信号。
作为一种可实施方式，所述中间量u’还直接作为计算调速器的输入信号的第一参考量；
调速器的输出分为两路，一路直接作为计算调速器的输入信号的第二参考量，另一路直接作为计算汽轮机的输入信号的第一参考量；其中，调速器的输入信号的第一参考量与第二参考量作减法运算后作为调速器的输入信号；
汽轮机的实际输出还分别直接作为计算汽轮机的输入信号的第二参考量以及直接作为计算电力系统的输入信号的第一参考量；其中，汽轮机的输入信号的第一参考量与第二参考量作减法运算后作为汽轮机的输入信号；
电力系统的输出分为两路，一路直接作为计算电力系统的输入信号的第二参考量，另一路经过K p 倍增益后即为控制系统的输出Δf；将负荷扰动ΔP d 作为计算电力系统的输入信号的第三参考量；其中，电力系统的输入信号的第一参考量依次与第二参考量和第三参考量作减法运算后作为电力系统的输入信号；K p 为发电机增益。
实施例4
上述实施例中的基于 LADRC的汽轮发电机控制系统中，存在一个可手动调节的静态补偿系数k 1，但是，随着该补偿系数的变大，系统的控制性能也会受到影响，故需要确认该补偿系数的取值范围。因此，对应上述实施例中的基于 LADRC的汽轮发电机控制系统，本发明实施例还提供一种基于LADRC的汽轮发电机控制系统的静态补偿系统确定方法，包括以下步骤：
S201：获得非线性元件的描述函数N 1(X)，并绘出负倒描述函数-1/N 1(X)的曲线，具体包括；具有饱和非线性特性的描述函数N 1(X)如公式（3），其中，饱和非线性的线性区宽度和线性区特性的斜率分别为0.0017和1：
Figure 525912DEST_PATH_IMAGE010
(3)
且有：
Figure 441915DEST_PATH_IMAGE011
其中，X表示非线性的振幅，K表示线性区特性的斜率，a表示饱和非线性的线性区宽度。
S202：求传递函数G 1(jq)，具体为：对于图6所示单区域汽轮发电机控制系统，定义
Figure 897167DEST_PATH_IMAGE012
其中，
Figure 265832DEST_PATH_IMAGE013
是LADRC的传递函数；jq表示频率特性，
Figure 831942DEST_PATH_IMAGE014
表示线性环节的总传递函数，
Figure 540135DEST_PATH_IMAGE015
表示调速器的传递函数，
Figure 166288DEST_PATH_IMAGE016
和
Figure 22249DEST_PATH_IMAGE017
均表示汽轮机的传递函数，
Figure 126471DEST_PATH_IMAGE018
表示电力系统的传递函 数，
Figure 17067DEST_PATH_IMAGE019
表示电力系统的传递函数的倒数，
Figure 814122DEST_PATH_IMAGE020
表示调速器的传递函数与汽轮机的传递函数 的乘积，
Figure 891799DEST_PATH_IMAGE021
表示机组下降特性的倒数，
Figure 65291DEST_PATH_IMAGE022
表示静态补偿系数k 1，
Figure 762725DEST_PATH_IMAGE023
表示常数值1，T G 表示调 速器时间常数，T T 表示汽轮机时间常数， T P 表示发电机时间常数，s表示微分算子。
然后通过信号流图和梅逊公式可以获得如公式（4）所示的传递函数G 1(jq)：
Figure 730681DEST_PATH_IMAGE024
(4)
S203：绘出传递函数G 1(jq)的奈奎斯特曲线,然后应用奈奎斯特稳定判据进行分析，具体为：可通过MATLAB将G 1(jq)的奈奎斯特曲线和-1/N 1(X)的曲线绘制在同一个坐标系中，G 1(jq)的奈奎斯特曲线与负实轴的交点为
Figure 826813DEST_PATH_IMAGE025
，若G 1(jq)的奈奎斯特曲线不包围-1/N 1(X)的曲线，即
Figure 7259DEST_PATH_IMAGE026
，则该系统稳定；若G 1(jq)的奈奎斯特曲线包围-1/N 1(X)的曲线，即
Figure 138026DEST_PATH_IMAGE027
，则该系统不稳定。
为了证实本发明提供的控制系统及其补偿系数确定方法的可行性与准确性，本发明考虑具有如下模型参数的单区域电力系统:
Figure 745725DEST_PATH_IMAGE030
LADRC参数选取如下：
Figure 329153DEST_PATH_IMAGE031
其中，w c 是控制器带宽，w o 是观测器带宽，补偿系数取为k 1=1.14，为了显示LADRC在该补偿策略下的控制效果，当GRC=0.0017MW/s时，在t=1秒时加入阶跃信号ΔP d =0.01，系统的响应曲线如图7所示。黑色点线表示系统采用了误差补偿策略后的情况，黑色实线表示系统没有采用误差补偿策略的情况。可以看到，在不采取误差补偿策略时，系统由于受到汽轮机GRC的影响，其控制性能受到恶化，变得不稳定，而在采用误差补偿策略后，系统的控制性能得到了很好的改善和恢复。因此，该补偿策略是切实有效的。
为确定补偿系数的取值范围，先采用描述函数法进行稳定性分析，要使系统稳定，则需满足
Figure 250972DEST_PATH_IMAGE026
，进而可得k 1的取值范围0< k 1<1.24。而当k 1=1.24时，G 1(jq)的奈奎斯特曲线包围-1/N 1(X)的曲线，则系统不稳定，其结果如图8所示。我们通过手动调节补偿系数k 1，来确认其的取值范围，仿真结果如图9所示。可以看到，当k 1逐渐变大到1.24时，系统变得不稳定。仿真结果证实，采用描述函数法对具有汽轮机GRC的汽轮发电机控制系统进行稳定性分析，从而获取静态补偿系数的取值范围是可行的。
本发明针对存在调速器死区或汽轮机GRC的汽轮发电机控制系统，分别提出了基于线性自抗扰控制算法的误差补偿策略，使线性自抗扰控制算法能够快速地消除调速器死区或汽轮机GRC的影响。同时，也分别对应地提出采用描述函数法进行验证和获取该补偿策略中静态补偿系数的取值范围。仿真结果显示，本发明提出的误差补偿策略对改善系统的控制性能有良好的表现，且采用描述函数法获取补偿系数的取值范围，为该系数的选取提供参考也是可行的。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (8)

1.基于LADRC的汽轮发电机控制系统，所述控制系统具有调速器死区，包括调速器、汽轮机、电力系统和线性扩张状态观测器，其特征在于，所述控制系统还包括：误差补偿模块；
所述误差补偿模块，用于将调速器的理论输出与实际输出之间的误差反馈至线性扩张状态观测器，利用线性扩张状态观测器进行估计，使LADRC达到消除调速器死区的目的。
2.根据权利要求1所述的基于LADRC的汽轮发电机控制系统，其特征在于，还包括：状态反馈增益模块；所述误差补偿模块包括上、下两个补偿分支；其中：
控制系统的输出为频率变化Δf；将控制系统的输出作为线性扩张状态观测器的第一路输入信号；
线性扩张状态观测器的输出与原始输入信号u作减法运算后作为状态反馈增益模块的输入；
状态反馈增益模块的输出分两路：一路与控制系统的输出Δf经1/R倍增益后的结果Δf/R作减法运算后输出中间量u’；另一路作为误差补偿模块的下补偿分支的第一路输入信号；R为机组下降特性；
将所述中间量u’直接作为误差补偿模块的上补偿分支的第一路输入信号，将调速器的实际输出直接作为误差补偿模块的上补偿分支的第二路输入信号；
误差补偿模块的上补偿分支的第一路输入信号和第二路输入信号作减法运算后再经过k倍增益后作为误差补偿模块的下补偿分支的第二路输入信号；k为可手动调节的静态补偿系数；
误差补偿模块的下补偿分支的第一路输入信号和第二路输入信号作减法运算后作为线性扩张状态观测器的第二路输入信号。
3.根据权利要求2所述的基于LADRC的汽轮发电机控制系统，其特征在于，
所述中间量u’经过调速器死区后还作为计算调速器的输入信号的第一参考量；调速器的实际输出还分别直接作为计算调速器的输入信号的第二参考量以及直接作为计算汽轮机的输入信号的第一参考量；其中，调速器的输入信号的第一参考量与第二参考量作减法运算后作为调速器的输入信号；
汽轮机的输出分为两路，一路直接作为计算汽轮机的输入信号的第二参考量，另一路直接作为计算电力系统的输入信号的第一参考量；其中，汽轮机的输入信号的第一参考量与第二参考量作减法运算后作为汽轮机的输入信号；
电力系统的输出分为两路，一路直接作为计算电力系统的输入信号的第二参考量，另一路经过K p 倍增益后即为控制系统的输出Δf；将负荷扰动ΔP d 作为计算电力系统的输入信号的第三参考量；其中，电力系统的输入信号的第一参考量依次与第二参考量和第三参考量作减法运算后作为电力系统的输入信号；K p 为发电机增益。
4.如权利要求1至3任一所述的基于 LADRC的汽轮发电机控制系统的静态补偿系数确定方法，其特征在于，包括：
步骤1：获得非线性元件的描述函数N(X)，并绘出负倒描述函数-1/N(X)的曲线，具体包括：具有死区非线性特性的描述函数N(X)如公式（1），其中，死区非线性的死区宽度和线性输出特性的斜率分别为0.1和1：
Figure 364923DEST_PATH_IMAGE001
(1)
且有：
Figure 346786DEST_PATH_IMAGE002
其中，X表示非线性的振幅，K表示线性输出特性的斜率，Δ表示死区非线性的死区宽度；
步骤2：求传递函数G(jq)，具体为：定义
Figure 417510DEST_PATH_IMAGE003
其中，G L 是LADRC的传递函数；jq表示频率特性，g 1表示线性环节的总传递函数，g 2表示调速器的传递函数，g 3表示汽轮机的传递函数，g 4表示电力系统的传递函数，g 5表示机组下降特性的倒数，g 6表示静态补偿系数k，T G 表示调速器时间常数，T T 表示汽轮机时间常数，T P 表示发电机时间常数，s表示微分算子；
然后通过信号流图和梅逊公式可以获得如公式（2）所示的传递函数G(jq)：
Figure 939758DEST_PATH_IMAGE004
(2)
步骤3：绘出传递函数G(jq)的奈奎斯特曲线，然后应用奈奎斯特稳定判据进行分析，具体为：将G(jq)的奈奎斯特曲线和-1/N(X)的曲线绘制在同一个坐标系中，G(jq)的奈奎斯特曲线与负实轴的交点为
Figure 779538DEST_PATH_IMAGE005
，若G(jq)的奈奎斯特曲线不包围-1/N(X)的曲线，即
Figure 994619DEST_PATH_IMAGE006
，则该系统稳定；若G(jq)的奈奎斯特曲线包围-1/N(X)的曲线，即
Figure 287060DEST_PATH_IMAGE007
，则该系统不稳定。
5.基于LADRC的汽轮发电机控制系统，所述控制系统具有汽轮机发电速率约束，包括调速器、汽轮机、电力系统和线性扩张状态观测器，其特征在于，所述控制系统还包括：误差补偿模块；
所述误差补偿模块，用于将汽轮机的理论输出与实际输出之间的误差作为外部扰动反馈至线性扩张状态观测器，利用线性扩张状态观测器进行估计，使LADRC消除汽轮机发电速率约束的影响。
6.根据权利要求5所述的基于LADRC的汽轮发电机控制系统，其特征在于，还包括：状态反馈增益模块；所述误差补偿模块包括上、下两个补偿分支；其中：
控制系统的输出为频率变化Δf；将控制系统的输出作为线性扩张状态观测器的第一路输入信号；
线性扩张状态观测器的输出与原始输入信号u作减法运算后作为状态反馈增益模块的输入；
状态反馈增益模块的输出分两路：一路与控制系统的输出Δf经1/R 倍增益后的结果Δf/R作减法运算后输出中间量u’；另一路作为误差补偿模块的下补偿分支的第一路输入信号；R为机组下降特性；
将所述中间量u’依次经过
Figure 612999DEST_PATH_IMAGE008
和
Figure 307285DEST_PATH_IMAGE009
倍增益后作为误差补偿模块的上补偿分支的第一路输入信号，将汽轮机的实际输出直接作为误差补偿模块的上补偿分支的第二路输入信号；其中，T G 表示调速器时间常数，T T 表示汽轮机时间常数；
误差补偿模块的上补偿分支的第一路输入信号和第二路输入信号作减法运算后再经过k 1倍增益后作为误差补偿模块的下补偿分支的第二路输入信号；k 1为可手动调节的静态补偿系数；
误差补偿模块的下补偿分支的第一路输入信号和第二路输入信号作减法运算后作为线性扩张状态观测器的第二路输入信号。
7.根据权利要求6所述的基于LADRC的汽轮发电机控制系统，其特征在于，
所述中间量u’还直接作为计算调速器的输入信号的第一参考量；
调速器的输出分为两路，一路直接作为计算调速器的输入信号的第二参考量，另一路直接作为计算汽轮机的输入信号的第一参考量；其中，调速器的输入信号的第一参考量与第二参考量作减法运算后作为调速器的输入信号；
汽轮机的实际输出还分别直接作为计算汽轮机的输入信号的第二参考量以及直接作为计算电力系统的输入信号的第一参考量；其中，汽轮机的输入信号的第一参考量与第二参考量作减法运算后作为汽轮机的输入信号；
电力系统的输出分为两路，一路直接作为计算电力系统的输入信号的第二参考量，另一路经过K p 倍增益后即为控制系统的输出Δf；将负荷扰动ΔP d 作为计算电力系统的输入信号的第三参考量；其中，电力系统的输入信号的第一参考量依次与第二参考量和第三参考量作减法运算后作为电力系统的输入信号；K p 为发电机增益。
8.如权利要求5至7任一所述的基于LADRC的汽轮发电机控制系统的静态补偿系统确定方法，其特征在于，包括：
步骤1：获得非线性元件的描述函数N 1(X)，并绘出负倒描述函数-1/N 1(X)的曲线，具体包括；具有饱和非线性特性的描述函数N 1(X)如公式（3），其中，饱和非线性的线性区宽度和线性区特性的斜率分别为0.0017和1：
Figure 630951DEST_PATH_IMAGE010
(3)
且有：
Figure 676267DEST_PATH_IMAGE011
其中，X表示非线性的振幅，K表示线性区特性的斜率，a表示饱和非线性的线性区宽度；
步骤2：求传递函数G 1(jq)，具体为：定义
Figure 71476DEST_PATH_IMAGE012
其中，是LADRC的传递函数；jq表示频率特性，表示线性环节的总传递函数，
Figure 178607DEST_PATH_IMAGE015
表 示调速器的传递函数，
Figure 784031DEST_PATH_IMAGE016
和
Figure 452910DEST_PATH_IMAGE017
均表示汽轮机的传递函数，
Figure 180695DEST_PATH_IMAGE018
表示电力系统的传递函数，
Figure 669445DEST_PATH_IMAGE019
表示电力系统的传递函数的倒数，
Figure 406457DEST_PATH_IMAGE020
表示调速器的传递函数与汽轮机的传递函数的乘 积，
Figure 929842DEST_PATH_IMAGE021
表示机组下降特性的倒数，
Figure 560019DEST_PATH_IMAGE022
表示静态补偿系数k 1，
Figure 801644DEST_PATH_IMAGE023
表示常数值1，T G 表示调速器 时间常数，T T 表示汽轮机时间常数， T P 表示发电机时间常数，s表示微分算子；
然后通过信号流图和梅逊公式可以获得如公式（4）所示的传递函数G 1(jq)：
Figure 14451DEST_PATH_IMAGE024
(4)
步骤3：绘出传递函数G 1(jq)的奈奎斯特曲线,然后应用奈奎斯特稳定判据进行分析，具体为：将G 1(jq)的奈奎斯特曲线和-1/N 1(X)的曲线绘制在同一个坐标系中，G 1(jq)的奈奎斯特曲线与负实轴的交点为
Figure 392343DEST_PATH_IMAGE025
，若G 1(jq)的奈奎斯特曲线不包围-1/N 1(X)的曲线，即
Figure 993088DEST_PATH_IMAGE026
，则该系统稳定；若G 1(jq)的奈奎斯特曲线包围-1/N 1(X)的曲线，即
Figure 925272DEST_PATH_IMAGE027
，则该系统不稳定。
CN202110945603.9A 2021-08-18 2021-08-18 基于ladrc的汽轮发电机控制系统及其静态补偿系数确定方法 Pending CN113406888A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110945603.9A CN113406888A (zh) 2021-08-18 2021-08-18 基于ladrc的汽轮发电机控制系统及其静态补偿系数确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110945603.9A CN113406888A (zh) 2021-08-18 2021-08-18 基于ladrc的汽轮发电机控制系统及其静态补偿系数确定方法

Publications (1)

Publication Number Publication Date
CN113406888A true CN113406888A (zh) 2021-09-17

Family

ID=77688795

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110945603.9A Pending CN113406888A (zh) 2021-08-18 2021-08-18 基于ladrc的汽轮发电机控制系统及其静态补偿系数确定方法

Country Status (1)

Country Link
CN (1) CN113406888A (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030199997A1 (en) * 2002-04-18 2003-10-23 Zhiqiang Gao Scaling and parameterizing a controller
CN105759616A (zh) * 2016-04-15 2016-07-13 浙江工业大学 考虑死区特性的伺服系统有限时间控制方法
CN108828943A (zh) * 2018-06-20 2018-11-16 南通航运职业技术学院 一种带扰动补偿和有限时间收敛的自抗扰控制方法
CN111654047A (zh) * 2020-06-12 2020-09-11 昆明理工大学 一种基于自抗扰的抽水蓄能与电化学储能联合参与电网负荷频率控制方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030199997A1 (en) * 2002-04-18 2003-10-23 Zhiqiang Gao Scaling and parameterizing a controller
CN105759616A (zh) * 2016-04-15 2016-07-13 浙江工业大学 考虑死区特性的伺服系统有限时间控制方法
CN108828943A (zh) * 2018-06-20 2018-11-16 南通航运职业技术学院 一种带扰动补偿和有限时间收敛的自抗扰控制方法
CN111654047A (zh) * 2020-06-12 2020-09-11 昆明理工大学 一种基于自抗扰的抽水蓄能与电化学储能联合参与电网负荷频率控制方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
常帅兵等: "具有调速器死区的LFC系统的广义自抗扰控制", 《系统仿真学报》 *
常帅兵等: "电力系统供电负荷稳定性优化控制研究", 《计算机仿真》 *
王划一等: "《自动控制原理 第3版》", 30 June 2017, 国防工业出版社 *

Similar Documents

Publication Publication Date Title
US8120302B2 (en) Servo control apparatus
Li et al. A novel frequency support control method for PV plants using tracking LQR
Rinaldi et al. Third order sliding mode observer-based approach for distributed optimal load frequency control
Aravindh et al. Design of observer-based non-fragile load frequency control for power systems with electric vehicles
Bernard et al. Decentralized load frequency control in an interconnected power system using Coefficient Diagram Method
CN110943613B (zh) 一种含干扰和电流约束的直流升压变换器控制方法
CN108964127B (zh) 一种双馈风力发电系统故障穿越的控制方法
Liu Predictive control of high-order fully actuated nonlinear systems with time-varying delays
Roy et al. Robust adaptive backstepping excitation controller design for simple power system models with external disturbances
Yousefian et al. Energy function inspired value priority based global wide-area control of power grid
Xiahou et al. Control of DFIG’s rotor-side converter with decoupling of current loops using observer-based fractional-order sliding-mode regulators
Chojaa et al. Optimization of DFIG wind turbine power quality through adaptive fuzzy control
CN110209122B (zh) 一种多轴运动平台的控制方法、装置、介质及设备
JP6067289B2 (ja) 電力系統の縮約モデル作成装置、作成方法及び作成プログラム
Jaramillo-Lopez et al. Adaptive control for a class of uncertain nonlinear systems: application to photovoltaic control systems
Amokrane et al. State observation of unknown nonlinear SISO systems based on virtual input estimation
CN113406888A (zh) 基于ladrc的汽轮发电机控制系统及其静态补偿系数确定方法
Verij Kazemi et al. Adaptive frequency control support of a DFIG based on second‐order derivative controller using data‐driven method
CN109004658B (zh) 一种互联电力系统负荷频率控制方法及装置
CN114362132B (zh) 基于分数阶滑模的互联电力系统负荷频率控制器设计方法
Aggarwal et al. Frequency control enhancement of an islanded microgrid using disturbance-observer-based approach
CN107248742B (zh) 一种无源鲁棒控制方法
CN114188955B (zh) 基于自适应学习的互联微电网复合频率控制方法及装置
CN113078646A (zh) 一种计及逆变电源控制器饱和特性的故障暂态电流解析方法
Kanchanaharuthai et al. A composite nonlinear controller for higher-order models of synchronous generators under external disturbances

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210917

RJ01 Rejection of invention patent application after publication