CN109821646A - 一种磨煤机出口煤粉流量软测量方法 - Google Patents

一种磨煤机出口煤粉流量软测量方法 Download PDF

Info

Publication number
CN109821646A
CN109821646A CN201910030816.1A CN201910030816A CN109821646A CN 109821646 A CN109821646 A CN 109821646A CN 201910030816 A CN201910030816 A CN 201910030816A CN 109821646 A CN109821646 A CN 109821646A
Authority
CN
China
Prior art keywords
coal
flow
air
coal pulverizer
pulverizer
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
CN201910030816.1A
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.)
BEIJING GUODIAN ZHISHEN CONTROL TECHNOLOGY Co Ltd
North China Electric Power University
Original Assignee
BEIJING GUODIAN ZHISHEN CONTROL TECHNOLOGY Co Ltd
North China Electric Power 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 BEIJING GUODIAN ZHISHEN CONTROL TECHNOLOGY Co Ltd, North China Electric Power University filed Critical BEIJING GUODIAN ZHISHEN CONTROL TECHNOLOGY Co Ltd
Priority to CN201910030816.1A priority Critical patent/CN109821646A/zh
Publication of CN109821646A publication Critical patent/CN109821646A/zh
Pending legal-status Critical Current

Links

Landscapes

  • Disintegrating Or Milling (AREA)
  • Crushing And Grinding (AREA)

Abstract

本发明公开了一种磨煤机出口煤粉流量软测量方法,所述测量方法包括以下步骤:步骤1:基于一次风的质量平衡和能量平衡方程,建立磨煤机入口一次风温度和流量模型;步骤2:基于磨煤机内部煤的质量平衡方程,建立磨煤机内原煤量和煤粉量模型;步骤3:基于磨煤机入口一次风温度、一次风流量和磨煤机内煤粉量,构建磨出口煤粉流量软测量方法;本发明实现了磨煤机出口煤粉流量的在线估计,提高了制粉系统出力控制精度和调整精细化配风。

Description

一种磨煤机出口煤粉流量软测量方法
技术领域
本发明涉及煤粉流量软测量技术领域,特别是涉及磨煤机出口煤粉流量软测量方法。
背景技术
我国“富煤、少气、贫油”的特征,导致了以煤炭为主的一次能源布局。长期以来,我国煤炭消费量在一次能源消费总量中所占比例约为70%,煤炭消费总量的50%以上来自于以火力发电为主的电力行业。截止2016年末,全国总发电装机容量到达16.5亿千瓦,其中,火电10.5亿千瓦,占总装机容量的64.0%。火电作为我国电力供应的重要保障,对实现节能优先的目标具有决定性作用,也是集中治理燃煤污染、解决环境问题的关键。制粉系统作为火电机组煤炭的输入端,其控制品质直接关系到锅炉燃烧的安全性、稳定性、经济性、环保性和快速性,尤其在新能源电力系统的大环境下,面对电网的频繁调峰调频和煤炭清洁高效利用的双重需求,制粉系统的控制性能至关重要。
目前,在火电厂制粉系统的出力一般通过控制各给煤机转速进行间接控制,控制方法本身较为粗犷且存在一定的滞后和惯性,在这种控制方式下,燃烧系统的风煤比不是总风量与入炉煤量的比值,而是总风量与总给煤量的比值,导致锅炉燃烧过程无法实现精准配风。近年来,随着一系列煤粉流量的检测方法的出现,使得解决这一问题成为可能,但这些方法或对检测元件产生较强的冲刷磨损,或对安装位置要求十分苛刻,均不适用实际工程应用。
因此希望有一种磨煤机出口煤粉流量软测量方法能有效解决现有技术中存在的问题。
发明内容
本发明公开了一种磨煤机出口煤粉流量软测量方法,所述测量方法包括以下步骤:
步骤1:基于一次风的质量平衡和能量平衡方程,建立磨煤机入口一次风温度和流量模型;
步骤2:基于磨煤机内部煤的质量平衡方程,建立磨煤机内原煤量和煤粉量模型;
步骤3:基于磨煤机入口一次风温度和流量、磨煤机内原煤量和煤粉量,构建磨煤机出口煤粉流量软测量方法。
优选地,所述步骤1的磨煤机入口一次风温度和流量模型的微分方程形式为:
式中,Wair为磨煤机入口一次风流量;T1为冷热风门至一次风流量惯性常数;为冷风门最大风流量;为热风门最大风流量;uL为冷风门开度;uH为热风门开度;θin为磨煤机入口一次风温度;T2为冷热风门至一次风温度惯性常数;CL为冷风比热容;θL为冷风温度,即空预器入口风温;CH为热风比热容;θH为热风温度,即空预器出口一次风温;Cin为磨煤机入口一次风比热容。
优选地,所述冷风比热容CL、热风比热容CH和磨煤机入口一次风比热容Cin分别为冷风温度、热风温度、磨煤机入口温度的函数为:
其中,CL的计算过程中分别取温度为0℃和25℃,湿度为d=10g/kg对应的比热容,即1.011kJ/kg·℃和1.012kJ/kg·℃;CH的计算过程中温度分别为200℃和400℃,干空气对应的比热容,即1.011,kJ/kg·℃和1.028,kJ/kg·℃;Cin的计算过程中温度为200℃和300℃,干空气对应的比热容,即1.011kJ/kg·℃和1.019kJ/kg·℃。
优选地,所述步骤2的磨煤机内原煤量和煤粉量模型的微分方程形式为:
I=K1Mc+K2Mpf+K3
式中,Mc为磨煤机内原煤质量;Wc为给煤量;Wpf为磨煤机内煤粉质量;Wpf为磨煤机出口煤粉流量;Kconv为单位时间内原煤向煤粉的转换系数。
优选地,所述磨煤机出口煤粉流量Wpf为一次风携带出磨煤机的煤粉质量流量,其正比于磨煤机内部的煤粉质量和一次风差压,一次风差压与磨煤机入口一次风温度和流量的函数关系为:
Wpf=KpfΔPpaMpf
其中,ΔPpa为一次风差压;Kpf为煤粉流量系数。
优选地,所述磨煤机内原煤质量Mc、单位时间内原煤向煤粉的转换系数Kconv、煤粉流量系数Kpf和磨煤机内煤粉质量Mpf通过磨电流公式进行辨识,磨电流公式为:
I=K1Mpf+K2Mc+K3
其中,I为磨电流,K1、K2、K3为模型待辨识参数。
优选地,所述步骤3的磨煤机出口煤粉流量软测量为基于状态量的代数方程,取磨煤机内原煤量和煤粉量模型中的中间变量Wpf作为磨煤机出口煤粉流量估计信号。
本发明公开的磨煤机出口煤粉流量软测量方法具有以下有益效果:
1.测量精确高。磨煤机是制粉系统中相对独立的环节,基于磨煤机的质量和能量平衡方程能够建立高精度的磨煤机出口煤粉流量估计模型,基于该估计模型提出的磨煤机出口煤粉流量软测量方法具有较高测量精度。
2.便于工程应用。建模过程采用机理分析和数据驱动相结合的建模方法,根据实际磨煤机对象的历史运行数据进行模型辨识,辨识得到的磨煤机模型能够反应实际对象的真实运行状态,有一定的工程应用和推广价值。
附图说明
图1是MPS型中速磨煤机结构示意图。
图2是模型辨识流程图。
图3是模型输出和实际输出对比验证示意图。
图4是磨出口煤粉流量软测量对比验证示意图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
由于磨煤机是制粉系统中相对独立的个体,其存在精确的质量和能量平衡关系,本发明专利正是基于这种质能平衡关系建立磨出口煤粉流量估计模型,实现磨煤机出口煤粉流量的在线估计,为提高制粉系统出力控制精度和调整精细化配风奠定基础。
本发明适用于火电厂MPS(Mill Parter Ship)型中速磨煤机,该磨煤机是由德国Babcock公司设计制造的一种辊盘式中速磨煤机,其电耗低、出力平稳且受碾磨件磨损影响小。
如图1所示,在MPS中速磨煤机内原煤由落煤管进入两个碾磨部件磨辊和磨盘之间,在压紧力的作用下受到挤压和碾磨而被粉碎成煤粉。一次风由热风和冷风经过各自的调节挡板混合至一定温度,通过风环进入中速磨煤机对煤粉进行干燥,同时将一定细度的煤粉带入碾磨区上部的粗粉分离器中进行分离,合格的细粉被一次风带出分离器出口的煤粉管道送到锅炉中燃烧,不合格的煤粉落入落煤管中继续碾磨。
磨煤机入口一次风温度和流量模型,所述的磨入口一次风流量和温度模型是基于一次风的质量平衡和能量平衡建立的,具体微分方程形式如下:
式中,Wair为磨入口一次风流量,kg/s;T1为冷热风门至一次风流量惯性常数,s;冷风门最大风流量,kg/s;热风门最大风流量,kg/s;uL冷风门开度,%;uH热风门开度,%;θin磨入口一次风温度,℃;T2冷热风门至一次风温度惯性常数,s;CL冷风比热容,kJ/kg·℃;θL冷风温度,即空预器入口风温,℃;CH热风比热容,kJ/kg·℃;θH热风温度,即空预器出口一次风温,℃;Cin磨入口一次风比热容,kJ/kg·℃。
所述的冷风比热容CL、热风比热容CH、磨入口一次风比热容Cin分别冷风温度、热风温度、磨入口温度的函数,可查表按下式计算:
其中,CL的计算过程中分别取温度为0℃;和25℃,湿度为d=10g/kg对应的比热容,即1.011kJ/kg·℃和1.012kJ/kg·℃;CH的计算过程中温度为200℃和400℃,干空气对应的比热容,即1.011,kJ/kg·℃和1.028,kJ/kg·℃;Cin的计算过程中温度为200℃和300℃,干空气对应的比热容,即1.011kJ/kg·℃和1.019kJ/kg·℃。
上述模型中,Wair,uL,uH,θin,θL,θH为可测量参数,CL,CH,Cin为中间过程参数,T1和T2为待辨识参数,可通过模型输出的一次风流量和一次风温度与实际输出数据对比辨识得到。
磨内原煤量和煤粉量模型,所述的磨内原煤量和煤粉量模型是基于磨内煤的质量平衡建立的,具体微分方程形式如下:
I=K1Mc+K2Mpf+K3
式中,Mc为磨煤机内原煤质量,kg;Wc为给煤量,kg/s;Mpf为磨内煤粉质量,kg;Wpf为磨出口煤粉流量,kg/s;Kconv为单位时间内原煤向煤粉的转换系数。其中,磨出口煤粉流量Wpf等于一次风携带出磨煤机的煤粉质量流量,其正比于磨煤机内部的煤粉质量和一次风差压,而一次风差压又与磨入口一次风温度和流量存在一定的函数关系,具体计算公式如下:
Wpf=KpfΔPpaMpf
其中,ΔPpa为一次风差压,mbar;Kpf为煤粉流量系数。
由于上述状态方程中的两个状态量(磨内原煤量和煤粉量)均不可测,导致模型辨识过程无法进行。为此,本发明专利构建了一个能够表征磨内存煤量的可测状态——磨电流来完成模型参数的辨识,具体公式如下:
I=K1Mpf+K2Mc+K3
其中,I为磨电流,A,K1、K2、K3为模型待辨识参数。。
上述模型中,Nc,I为可测量参数;Wc,Wpf,ΔPpa为中间过程参数;Mc,Kconv,Mpf,Kpf,K1,K2,K3为待辨识参数。
磨出口煤粉流量软测量方法,取磨内原煤量和煤粉量模型中的中间变量Wpf作为磨出口煤粉流量估计信号,由其表达式可见,其是一次风流量、一次风温度、磨内煤粉量三个状态量的代数计算形式。至此,基于模型状态估计提出一种磨出口煤粉流量软测量方法。
模型辨识及验证
本发明专利提出的磨出口煤粉流量软测量模型可采用多种模型辨识方法进行辨识和验证,本发明中采用实数编码的遗传算法对模型中的待定参数进行辨识,以实际输出和模型输出误差的函数作为个体适应度函数,种群大小为50,终止迭代次数为50,交叉概率为0.9,变异概率为0.3,具体辨识流程如图2所示。
其中,模型适应度函数为:
其中,N为测量数据点的个数;W1和W2为权重系数;Wair为模型输出的一次风流量,kg/s;为实际一次风流量,kg/s;Tin为模型输出的磨入口温度,℃;为实际磨入口温度,℃;I为实际磨电流,℃;为实际磨电流,A。
所述辨识获得的模型参数列表如下表所示:
表1磨出口煤粉流量软测量模型参数列表
为了验证模型在多工况下的准确性和有效性,连续5天的实际运行数据被用来对比验证,采样期间磨煤机负荷率为59.88%-97.59%。如图3所示,模型输出与实际输出基本吻合,所建模型能够准确描述磨煤机在中高负荷段的动态特性。
如图4所为了进一步验证磨出口煤粉流量估计信号的准确性,将模型软测量获得磨出口煤粉流量与实际给煤量进行对比验证,对比结果示磨出口煤粉流量和给煤量基本吻合,且磨出口煤粉流量相对于给煤量迟延约20秒,验证了磨出口煤粉流量软测量方法的准确性。
最后需要指出的是:以上实施例仅用以说明本发明的技术方案,而非对其限制。尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (7)

1.一种磨煤机出口煤粉流量软测量方法,其特征在于,所述测量方法包括以下步骤:
步骤1:基于一次风的质量平衡和能量平衡方程,建立磨煤机入口一次风温度和流量模型;
步骤2:基于磨煤机内部煤的质量平衡方程,建立磨煤机内原煤量和煤粉量模型;
步骤3:基于磨煤机入口一次风温度和流量、磨煤机内原煤量和煤粉量,构建磨煤机出口煤粉流量软测量方法。
2.根据权利要求1所述的磨煤机出口煤粉流量软测量方法,其特征在于:所述步骤1的磨煤机入口一次风温度和流量模型的微分方程形式为:
式中,Wair为磨煤机入口一次风流量;T1为冷热风门至一次风流量惯性常数;为冷风门最大风流量;为热风门最大风流量;uL为冷风门开度;uH为热风门开度;θin为磨煤机入口一次风温度;T2为冷热风门至一次风温度惯性常数;CL为冷风比热容;θL为冷风温度,即空预器入口风温;CH为热风比热容;θH为热风温度,即空预器出口一次风温;Cin为磨煤机入口一次风比热容。
3.根据权利要求2所述的磨煤机出口煤粉流量软测量方法,其特征在于:所述冷风比热容CL、热风比热容CH和磨煤机入口一次风比热容Cin分别为冷风温度、热风温度、磨煤机入口温度的函数为:
其中,CL的计算过程中分别取温度为0℃和25℃,湿度为d=10g/kg对应的比热容,即1.011kJ/kg·℃和1.012kJ/kg·℃;CH的计算过程中温度分别为200℃和400℃,干空气对应的比热容,即1.011kJ/kg·℃和1.028kJ/kg·℃;Cin的计算过程中温度为200℃和300℃,干空气对应的比热容,即1.011kJ/kg·℃和1.019kJ/kg·℃。
4.根据权利要求2所述的磨煤机出口煤粉流量软测量方法,其特征在于:所述步骤2的磨煤机内原煤量和煤粉量模型的微分方程形式为:
I=K1Mc+K2Mpf+K3
式中,Mc为磨煤机内原煤质量;Wc为给煤量;Mpf为磨煤机内煤粉质量;Wpf为磨煤机出口煤粉流量;Kconv为单位时间内原煤向煤粉的转换系数。
5.根据权利要求4所述的磨煤机出口煤粉流量软测量方法,其特征在于:所述磨煤机出口煤粉流量Wpf为一次风携带出磨煤机的煤粉质量流量,其正比于磨煤机内部的煤粉质量和一次风差压,一次风差压与磨煤机入口一次风温度和流量的函数关系为:
Wpf=KpfΔPpaMpf
其中,ΔPpa为一次风差压;Kpf为煤粉流量系数。
6.根据权利要求5所述的磨煤机出口煤粉流量软测量方法,其特征在于:所述磨煤机内原煤质量Mc、单位时间内原煤向煤粉的转换系数Kconv、煤粉流量系数Kpf和磨煤机内煤粉质量Mpf通过磨电流公式进行辨识,磨电流公式为:
I=K1Mpf+K2Mc+K3
其中,I为磨电流,K1、K2、K3为模型待辨识参数。
7.根据权利要求5所述的磨煤机出口煤粉流量软测量方法,其特征在于:所述步骤3的磨煤机出口煤粉流量软测量为基于状态量的代数方程,取磨煤机内原煤量和煤粉量模型中的中间变量Wpf作为磨煤机出口煤粉流量估计信号。
CN201910030816.1A 2019-01-14 2019-01-14 一种磨煤机出口煤粉流量软测量方法 Pending CN109821646A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910030816.1A CN109821646A (zh) 2019-01-14 2019-01-14 一种磨煤机出口煤粉流量软测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910030816.1A CN109821646A (zh) 2019-01-14 2019-01-14 一种磨煤机出口煤粉流量软测量方法

Publications (1)

Publication Number Publication Date
CN109821646A true CN109821646A (zh) 2019-05-31

Family

ID=66861566

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910030816.1A Pending CN109821646A (zh) 2019-01-14 2019-01-14 一种磨煤机出口煤粉流量软测量方法

Country Status (1)

Country Link
CN (1) CN109821646A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112547294A (zh) * 2020-11-23 2021-03-26 西安热工研究院有限公司 一种热态条件下中速磨煤机入口风量的获取方法
CN113420467A (zh) * 2021-08-23 2021-09-21 长沙理工大学 一种融合机理与数据驱动的风粉浓度测量方法
CN116078535A (zh) * 2023-01-09 2023-05-09 西安热工研究院有限公司 一种中速磨煤机及其工作方法

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4177951A (en) * 1978-06-28 1979-12-11 Combustion Engineering Inc. Pulverizer air flow and temperature control
US4540129A (en) * 1982-11-12 1985-09-10 The Babcock & Wilcox Company Pulverizer control system
JPH06137918A (ja) * 1992-10-28 1994-05-20 Chugoku Electric Power Co Inc:The 微粉炭流量の計測方法
CN105388765A (zh) * 2015-12-24 2016-03-09 东南大学 一种中速磨煤机的多变量推断预测控制方法
CN105571653A (zh) * 2015-12-17 2016-05-11 广东电网有限责任公司电力科学研究院 磨煤机风量测量方法和系统
CN107066657A (zh) * 2016-12-23 2017-08-18 北京源深节能技术有限责任公司 磨煤机出口温度优化节能效果评价方法及装置
CN108508747A (zh) * 2018-01-12 2018-09-07 华润电力(贺州)有限公司 配置双进双出磨煤机制粉系统的直流炉机组模型建模方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4177951A (en) * 1978-06-28 1979-12-11 Combustion Engineering Inc. Pulverizer air flow and temperature control
US4540129A (en) * 1982-11-12 1985-09-10 The Babcock & Wilcox Company Pulverizer control system
JPH06137918A (ja) * 1992-10-28 1994-05-20 Chugoku Electric Power Co Inc:The 微粉炭流量の計測方法
CN105571653A (zh) * 2015-12-17 2016-05-11 广东电网有限责任公司电力科学研究院 磨煤机风量测量方法和系统
CN105388765A (zh) * 2015-12-24 2016-03-09 东南大学 一种中速磨煤机的多变量推断预测控制方法
CN107066657A (zh) * 2016-12-23 2017-08-18 北京源深节能技术有限责任公司 磨煤机出口温度优化节能效果评价方法及装置
CN108508747A (zh) * 2018-01-12 2018-09-07 华润电力(贺州)有限公司 配置双进双出磨煤机制粉系统的直流炉机组模型建模方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
曾德良,高珊,胡勇: "《MPS型中速磨煤机建模与仿真》", 《动力工程学报》 *
高耀岿,胡勇,曾德良,刘吉臻: "《基于磨出口温度定值优化的制粉系统多变量预测控制》", 《动力工程学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112547294A (zh) * 2020-11-23 2021-03-26 西安热工研究院有限公司 一种热态条件下中速磨煤机入口风量的获取方法
CN113420467A (zh) * 2021-08-23 2021-09-21 长沙理工大学 一种融合机理与数据驱动的风粉浓度测量方法
CN116078535A (zh) * 2023-01-09 2023-05-09 西安热工研究院有限公司 一种中速磨煤机及其工作方法

Similar Documents

Publication Publication Date Title
CN109821646A (zh) 一种磨煤机出口煤粉流量软测量方法
CN101697179A (zh) 基于正反热平衡关系的电站锅炉燃煤热值趋势测算方法
CN108672025B (zh) 一种钢球磨煤制粉系统的多工况最优解耦控制方法
CN113033873B (zh) 基于入炉煤质测量的脱硫系统入口硫氧化物含量预测方法
CN109442463B (zh) 一种磨煤机煤粉细度调节方法
Gao et al. Modeling of a medium speed coal mill
CN102600959B (zh) 一种双进双出磨煤机制粉系统的协调控制系统及方法
Niemczyk et al. Derivation and validation of a coal mill model for control
Ataş et al. Numerical and experimental analysis of pulverized coal mill classifier performance in the Soma B Power Plant
CN102645523A (zh) 基于制粉系统热平衡的入炉煤收到基水分在线辨识方法
CN109459102B (zh) 一种中速磨煤机进口一次风量测量方法和系统
CN102639937B (zh) 用于监测和控制电厂的系统和关联的方法
CN102235893B (zh) 一次风量测量方法及装置
Pronobis et al. Optimisation of coal fineness in pulverised-fuel boilers
Gao et al. Optimization control of a pulverizing system on the basis of the estimation of the outlet coal powder flow of a coal mill
CN103816987A (zh) 一种双进双出磨煤机出粉量的计算方法
CN103807852B (zh) 煤粉锅炉的二次风配风装置及二次风配风方法
CN113182057B (zh) 一种中速磨煤机出力及单耗试验结果修正计算方法
CN105279573B (zh) 一种火电厂煤耗率经济性分析方法
Zeng et al. Modeling and control of ball mill system considering coal moisture
CN218076900U (zh) 一种低碳排放制粉系统
CN209109379U (zh) 基于煤粉细度在线调整的磨煤机出口温度控制系统
CN103699782B (zh) 一种中速磨制粉系统给煤量软测量方法
CN204084340U (zh) 煤粉锅炉的二次风配风装置
CN115475455A (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20190531