CN110416999A - 一种考虑风电参与调频的风电场聚合频率响应模型构建方法 - Google Patents

一种考虑风电参与调频的风电场聚合频率响应模型构建方法 Download PDF

Info

Publication number
CN110416999A
CN110416999A CN201910592347.2A CN201910592347A CN110416999A CN 110416999 A CN110416999 A CN 110416999A CN 201910592347 A CN201910592347 A CN 201910592347A CN 110416999 A CN110416999 A CN 110416999A
Authority
CN
China
Prior art keywords
wind
frequency
frequency response
power plant
parameter
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
CN201910592347.2A
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.)
Southeast University
Original Assignee
Southeast University
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 Southeast University filed Critical Southeast University
Priority to CN201910592347.2A priority Critical patent/CN110416999A/zh
Publication of CN110416999A publication Critical patent/CN110416999A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply
    • 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
    • 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/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/48Controlling the sharing of the in-phase component
    • 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
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Business, Economics & Management (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Wind Motors (AREA)

Abstract

本发明公开了一种考虑风电参与调频的风电场聚合频率响应模型构建方法,该方法通过单个风电机组一次调频控制策略和电力系统小信号分析原理,建立单个风电机组的频率响应模型,随后通过参数等值方法建立风电场机群聚合频率响应模型。该聚合模型无需建立复杂的电磁暂态模型,便能够对风电场具有不同参数的单个风机的频率响应特性进行等效的描述,在保证精度的基础上对原有方法计算效率进行了提高。

Description

一种考虑风电参与调频的风电场聚合频率响应模型构建方法
技术领域
本发明涉及电力系统安全稳定控制领域,具体涉及一种考虑风电参与调频的风电场聚合频率响应模型。
背景技术
当前电网中风电接入比例的不断上升,使风电机组在参与电网一次调频,提供一次调频资源方面具有更重要的意义,因此需要对风电机组的频率响应特性进行建模,以便分析风电机组在系统出现频率扰动状态下的动态响应过程。现有方法中对风电机组的惯量和下垂控制响应进行了等效建模,但对于频率控制回路与惯性参数间的耦合特性考虑不足,亦缺少对频率响应过程的传递函数描述。同时,风电场中不同风电机组运行状态及控制参数的不同,也对新的等效和聚合建模方法提出了需求。
因此,建立一个更加精确的风电场频率响应模型,对分析含风电的电力系统频率响应特征,显得很有必要。
发明内容
本发明提供一种考虑风电参与调频的风电场聚合频率响应模型构建方法,解决了现有技术中对风电场风机机群的一次调频过程频率响应特性描述方法的不足。
为达到上述目的,本发明采用的方法是:一种考虑风电参与调频的风电场聚合频率响应模型构建方法,包括如下步骤:
1)针对单个变速风电机组,建立采用虚拟惯量控制和基于桨距角降载控制
的协调频率控制模型,通过小信号分析建立单个风电机组的频率响应模型;
2)针对多个风电机组组成的风电场,通过加权平均等值建立多风电机组聚
合的频率响应模型。
作为本发明的一种改进,所述的虚拟惯量控制策略通过惯性和下垂控制对系统的频率变化做出响应,调整风电机组的电磁功率输出,基于虚拟惯量的控制方法的传递函数表示为
基于桨距角的频率控制策略通过调整桨距角来改变风电机组功率输出,进而对系统频率变化做出响应,基于桨距角频率控制下的桨距角偏差和频率偏差关系式为
Δβ=kbΔf (2)
风电机组的机械转矩可表示为:
根据小信号分析原理,求得的风电机组机械转矩变化量可表示为:
式(3)代入式(4)可得:
进一步,求得的风电机组电磁转矩变化量可表示为:
风电机组的功率摇摆方程为:
2HtsΔω(s)=ΔTm(s)-ΔTe(s) (7)
通过公式(5)(6)和(7),得到频率响应模型传递函数表示:
其中Ht为风机惯性时间常数,ω为转子转速,Pe为风机的有功出力,Pe1为虚拟惯量控制的功率变化量,Δf为频率变化量,kv和Rv分别为虚拟惯量参数和下垂控制参数,kb为缩放因子,ρ为空气密度,r为叶片长度,v为风速,Cp为风机的性能系数,λ为叶尖速比,β为桨距角,Pbase为风机的额定功率,kp为缩放因子,
作为本发明的一种改进,所述步骤2)中风电场风机机群聚合频率响应模型,其建模过程为:
2-1)将多个风机进行参数等值与聚合,用一个等值参数的响应特性代替多个风机的响应特性,可用式(10)表达:
其中ai,bi,ci和qi为风机机群中第i个风机频率响应模型的参数,kwi表述第i个风机额定功率占风电场总额定功率的比例,N为风电场机群中风机的总数;a,b,c和q为风电场聚合频率响应模型的参数。
将式(10)分解为:
其中
将式(10)写为如下的形式:
将上述各式的参数k,d,p,q表示为一个集合X={k,d,p,q};
2-2)确定等值参数X;
定义加权平均系数λi,并将风电场聚合参数X和各单机参数Xi间的关系表示为
比较式(13)和式(14),可知对于X中的参数k,d,p,令λi=kwi,易得
因此,根据加权平均思想可得
因此等值参数X={k,d,p,q}可以计算求得。
有益效果:
本发明提供了针对风电场风机的单机频率响应模型和风电场机群的聚合频率响应模型,考虑风电机组参与调频的风电场频率响应特性进行了公式化描述,对理解风电机组在功率扰动过程中的频率响应过程、电力系统的频率稳定在线预测具有指导作用。
附图说明
图1为本发明的系统流程图;
图2为本发明内容中风电场中风电机组综合虚拟惯量控制和基于桨距角降载控制的组合频率控制策略模型;
图3为本发明内容中的对风电场多机频率响应模型等值过程的示意图;
图4为本发明实施例实施效果比较图。
具体实施方式
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将结合附图对本发明的具体实施方式作进一步地说明。
如图1所述为本发明的系统流程图本发明实施例提供的一种考虑风电参与调频的风电场聚合频率响应模型构建方法,包括如下步骤:
1)针对单个变速风电机组,建立采用虚拟惯量控制和基于桨距角降载控制
的协调频率控制模型,通过小信号分析建立单个风电机组的频率响应模型;
2)针对多个风电机组组成的风电场,通过加权平均等值建立多风电机组聚
合的频率响应模型。
对于步骤1)所述采用虚拟惯量控制和基于桨距控制的降载控制的控制策略模型如附图1所示。其中Ht为风机惯性时间常数,ω为转子转速,Tm和Te分别为机械转矩和电磁转矩,Pe为风机的有功出力,Δf为频率变化量,Δβ为桨距角变化量,kv和Rv分别为虚拟惯量参数和下垂控制参数,kb为缩放因子。如图2所示,虚拟惯量控制回路通过惯性和下垂控制对系统的频率变化做出响应,调整风电机组的电磁功率输出,基于桨距的控制通过调整桨距角对系统频率变化做出响应,经过风机桨距角调整风电机功率输出。
步骤1)中所述的虚拟惯量控制方式的传递函数为
Pe1为虚拟惯量控制的功率变化量。
步骤1)所述的基于桨距角频率控制下的桨距角偏差和频率偏差关系式为
Δβ=kbΔf (2)
步骤1)所述的风电机组的机械转矩可表示为:
步骤1)所述的风电机组机械转矩变化量可表示为:
步骤1)式(3)代入式(4)可得:
步骤1)所述的风电机组电磁转矩变化量可表示为:
步骤1)所述的风电机组功率摇摆方程为:
2HtsΔω(s)=ΔTm(s)-ΔTe(s) (7)
步骤1)所述单个风电机组频率响应模型融合了虚拟惯量控制和桨距控制方式,整体描述了风机的电磁功率变化量与频率变化量的对应关系。通过公式(5)(6)和(7),本频率响应模型的传递函数如下式:
其中ρ为空气密度,r为叶片长度,v为风速,Cp为风机的性能系数,λ为叶尖速比,β为桨距角,Pbase为风机的额定功率,kp为缩放因子,
步骤2)所述风电场风机机群聚合频率响应模型,其建模过程为:
(2-1)将多个风机进行参数等值与聚合,如附图3所示,用一个等值参数
的响应特性代替多个风机的响应特性,可用式(10)表达:
其中ai,bi,ci和qi为风机机群中第i个风机频率响应模型的参数,kwi表述第i个风机额定功率占风电场总额定功率的比例,N为风电场机群中风机的总数;a,b,c和q为风电场聚合频率响应模型的参数。
为简化分析,将式(11)分解为:
其中
因此式(10)也可写为如下的形式:
将上述各式的参数k,d,p,q表示为一个集合X={k,d,p,q}。
(2-1)确定等值参数X。
为使等值后的单机频率响应规律能尽可能地接近原始多机群聚合后的频率响应规律,需要选择最优的X取值。考虑到单个风电机组的功率占风电场总功率的比重对于这一风电机组的特性对整体频率响应特性有重要影响,定义加权平均系数λi,并将风电场聚合参数X和各单机参数Xi间的关系表示为
比较式(13)和式(14),可知对于X中的参数k,d,p,令λi=kwi,易得
因此,根据加权平均思想可得
因此等值参数X={k,d,p,q}可以计算求得。
将本发明建立的风电场聚合频率响应模型的动态频率响应特性分别与聚合之前的多机频率响应模型以及实际详细模型的动态频率响应特性进行比较,即可对本发明提出的聚合频率响应模型进行有效性验证。如图4所示,由结果可以看出,本发明提出的聚合频率响应模型与聚合之前的多机频率响应模型以及实际详细模型在频率变化率、频率最低点以及频率稳态值等频率特性指标具有较高的吻合度,表明本发明提出的聚合频率响应模型具有较高的精度。

Claims (3)

1.一种考虑风电参与调频的风电场聚合频率响应模型构建方法,其特征在于,包括如下步骤:
1)针对单个变速风电机组,建立采用虚拟惯量控制和基于桨距角降载控制的协调频率控制模型,通过小信号分析建立单个风电机组的频率响应模型;
2)针对多个风电机组组成的风电场,通过加权平均等值建立多风电机组聚合的频率响应模型。
2.根据权利要求1所述的一种考虑风电参与调频的风电场聚合频率响应模型构建方法,其特征在于:所述的虚拟惯量控制回路通过惯性和下垂控制对系统的频率变化做出响应,调整风电机组的电磁功率输出,基于虚拟惯量的控制方法的传递函数表示为
基于桨距角的频率控制策略通过调整桨距角来改变风电机组功率输出,进而对系统频率变化做出响应,基于桨距角控制方法的表达式为
Δβ=kbΔf (2)
风电机组的机械转矩可表示为:
根据小信号分析原理,求得的风电机组机械转矩变化量可表示为:
式(3)代入式(4)可得:
进一步,求得的风电机组电磁转矩变化量可表示为:
风电机组的功率摇摆方程为:
2HtsΔω(s)=ΔTm(s)-ΔTe(s) (7)
通过公式(5)(6)和(7),得到频率响应模型传递函数表示:
式中
其中Ht为风机惯性时间常数,ω为转子转速,Pe为风机的有功出力,Pe1为虚拟惯量控制的功率变化量,Δf为频率变化量,kv和Rv分别为虚拟惯量参数和下垂控制参数,kb为缩放因子,ρ为空气密度,r为叶片长度,v为风速,Cp为风机的性能系数,λ为叶尖速比,β为桨距角,Pbase为风机的额定功率,kp为缩放因子,
3.根据权利要求1所述的一种考虑风电参与调频的风电场聚合频率响应模型构建方法,其特征在于:所述步骤2)中风电场风机机群聚合频率响应模型,其建模过程为:
2-1)将多个风机进行参数等值与聚合,用一个等值参数的响应特性代替多个风机的响应特性,可用式(10)表达:
其中ai,bi,ci和qi为风机机群中第i个风机频率响应模型的参数,kwi表述第i个风机额定功率占风电场总额定功率的比例,N为风电场机群中风机的总数;a,b,c和q为风电场聚合频率响应模型的参数。
将式(10)分解为:
其中
将式(10)写为如下的形式:
将上述各式的参数k,d,p,q表示为一个集合X={k,d,p,q};
2-2)确定等值参数X;
定义加权平均系数λi,并将风电场聚合参数X和各单机参数Xi间的关系表示为
比较式(13)和式(14),可知对于X中的参数k,d,p,令λi=kwi,易得,
因此,根据加权平均思想可得
从而计算得到等值参数X={k,d,p,q}。
CN201910592347.2A 2019-07-03 2019-07-03 一种考虑风电参与调频的风电场聚合频率响应模型构建方法 Pending CN110416999A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910592347.2A CN110416999A (zh) 2019-07-03 2019-07-03 一种考虑风电参与调频的风电场聚合频率响应模型构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910592347.2A CN110416999A (zh) 2019-07-03 2019-07-03 一种考虑风电参与调频的风电场聚合频率响应模型构建方法

Publications (1)

Publication Number Publication Date
CN110416999A true CN110416999A (zh) 2019-11-05

Family

ID=68358826

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910592347.2A Pending CN110416999A (zh) 2019-07-03 2019-07-03 一种考虑风电参与调频的风电场聚合频率响应模型构建方法

Country Status (1)

Country Link
CN (1) CN110416999A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111371104A (zh) * 2020-03-20 2020-07-03 国网内蒙古东部电力有限公司 一种基于风储联合发电系统的电网频率稳定控制方法
CN111541238A (zh) * 2020-04-10 2020-08-14 东方电气风电有限公司 一种风力发电机组调频验证模型建模方法
CN115833102A (zh) * 2022-12-08 2023-03-21 南方电网数字电网研究院有限公司 基于模型预测控制的风电场频率快速响应控制方法

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103840476A (zh) * 2014-03-04 2014-06-04 国家电网公司 一种基于单机等值模型的高周切机整定方法
US20150260159A1 (en) * 2014-03-11 2015-09-17 Gamesa Innovation & Technology, S.L. Wind turbine inertia control system
CN105656034A (zh) * 2016-02-29 2016-06-08 合肥工业大学 孤岛微电网频率控制的仿真建模方法
CN106849088A (zh) * 2017-02-17 2017-06-13 三峡大学 一种基于浆距控制的风电有功/频率耦合电力系统频率特性计算方法
CN107240918A (zh) * 2017-07-27 2017-10-10 江苏大学 一种风电接入的电力系统等值简化方法
CN107482652A (zh) * 2017-08-22 2017-12-15 山东大学 一种电力系统低频减载虚拟实验实现方法和系统
CN107679769A (zh) * 2017-10-25 2018-02-09 东南大学 含风电的电力系统频率响应模型建立方法和频率特性指标计算方法
US20190072070A1 (en) * 2016-03-01 2019-03-07 Wind Power Technology Center Of Gansu Electric Power Company Control system for enhancing frequency support of wind turbine generator
CN109560573A (zh) * 2018-12-10 2019-04-02 国网青海省电力公司 一种变速风电机组频率控制器参数的优化方法及装置

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103840476A (zh) * 2014-03-04 2014-06-04 国家电网公司 一种基于单机等值模型的高周切机整定方法
US20150260159A1 (en) * 2014-03-11 2015-09-17 Gamesa Innovation & Technology, S.L. Wind turbine inertia control system
CN105656034A (zh) * 2016-02-29 2016-06-08 合肥工业大学 孤岛微电网频率控制的仿真建模方法
US20190072070A1 (en) * 2016-03-01 2019-03-07 Wind Power Technology Center Of Gansu Electric Power Company Control system for enhancing frequency support of wind turbine generator
CN106849088A (zh) * 2017-02-17 2017-06-13 三峡大学 一种基于浆距控制的风电有功/频率耦合电力系统频率特性计算方法
CN107240918A (zh) * 2017-07-27 2017-10-10 江苏大学 一种风电接入的电力系统等值简化方法
CN107482652A (zh) * 2017-08-22 2017-12-15 山东大学 一种电力系统低频减载虚拟实验实现方法和系统
CN107679769A (zh) * 2017-10-25 2018-02-09 东南大学 含风电的电力系统频率响应模型建立方法和频率特性指标计算方法
CN109560573A (zh) * 2018-12-10 2019-04-02 国网青海省电力公司 一种变速风电机组频率控制器参数的优化方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
冯家铨: "基于桨距控制的电力系统频率响应特性", 《电力科学与工程》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111371104A (zh) * 2020-03-20 2020-07-03 国网内蒙古东部电力有限公司 一种基于风储联合发电系统的电网频率稳定控制方法
CN111371104B (zh) * 2020-03-20 2023-03-24 国网内蒙古东部电力有限公司 一种基于风储联合发电系统的电网频率稳定控制方法
CN111541238A (zh) * 2020-04-10 2020-08-14 东方电气风电有限公司 一种风力发电机组调频验证模型建模方法
CN111541238B (zh) * 2020-04-10 2022-05-24 东方电气风电股份有限公司 一种风力发电机组调频验证模型建模方法
CN115833102A (zh) * 2022-12-08 2023-03-21 南方电网数字电网研究院有限公司 基于模型预测控制的风电场频率快速响应控制方法
CN115833102B (zh) * 2022-12-08 2023-08-25 南方电网数字电网研究院有限公司 基于模型预测控制的风电场频率快速响应控制方法

Similar Documents

Publication Publication Date Title
CN106910142B (zh) 一种含风电有功-频率耦合作用的电力系统频率特性计算方法
CN103138256B (zh) 一种新能源电力消纳全景分析系统及方法
CN102545268B (zh) 一种在风电受限状态下的大电网有功功率实时控制方法
CN110416999A (zh) 一种考虑风电参与调频的风电场聚合频率响应模型构建方法
CN107800146A (zh) 兼顾一次调频和超低频振荡抑制的调速器参数优化方法
CN109560573B (zh) 一种变速风电机组频率控制器参数的优化方法及装置
CN108964032A (zh) 一种建立含风电的系统频率响应模型的方法
CN107482649A (zh) 一种基于分频控制的两域互联系统负荷频率控制方法
Pulgar Painemal Wind farm model for power system stability analysis
CN107679769A (zh) 含风电的电力系统频率响应模型建立方法和频率特性指标计算方法
CN107453410A (zh) 负荷扰动的双馈风机参与风柴微网调频控制方法
CN109936152A (zh) 高渗透率风电并网后的电网调频控制方法、风力发电机
CN106786759B (zh) 与同步发电机组相协调的风电机组一次调频方法
CN106712058A (zh) 双馈风机风电场参与电力系统一次调频的协调控制方法
CN103595046A (zh) 以柴油机侧为主的风柴混合电力系统负荷频率控制方法
CN110417032A (zh) 一种双馈风机参与系统调频的多目标优化控制方法
CN107346889B (zh) 考虑一二次调频及最小频率偏差的负荷削减优化模型构建方法
CN112072677A (zh) 一种基于分数阶pid的抽水蓄能与电化学储能联合参与电网负荷频率控制方法
Malmir et al. Controlling megawatt class WECS by ANFIS network trained with modified genetic algorithm
OĞUZ et al. Adaptive neuro-fuzzy inference system to improve the power quality of variable-speed wind power generation system
CN114759618A (zh) 基于动态能量确定双馈风电并网系统稳定性的方法及系统
CN113241778A (zh) 一种基于多区域互联电网的agc控制方法
CN103606939A (zh) 基于滑模控制的风柴混合电力系统负荷频率控制方法
CN103595045B (zh) 风机参与调频的风柴混合电力系统负荷频率协调控制方法
CN110968958B (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20191105

RJ01 Rejection of invention patent application after publication