CN104977847B - 一种面向常减压优化的稳态工况判别方法 - Google Patents

一种面向常减压优化的稳态工况判别方法 Download PDF

Info

Publication number
CN104977847B
CN104977847B CN201510379120.1A CN201510379120A CN104977847B CN 104977847 B CN104977847 B CN 104977847B CN 201510379120 A CN201510379120 A CN 201510379120A CN 104977847 B CN104977847 B CN 104977847B
Authority
CN
China
Prior art keywords
time
steady
steady state
state
real
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
CN201510379120.1A
Other languages
English (en)
Other versions
CN104977847A (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.)
NANJING RICHISLAND INFORMATION ENGINEERING Co Ltd
Original Assignee
NANJING RICHISLAND INFORMATION ENGINEERING 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 NANJING RICHISLAND INFORMATION ENGINEERING Co Ltd filed Critical NANJING RICHISLAND INFORMATION ENGINEERING Co Ltd
Priority to CN201510379120.1A priority Critical patent/CN104977847B/zh
Publication of CN104977847A publication Critical patent/CN104977847A/zh
Application granted granted Critical
Publication of CN104977847B publication Critical patent/CN104977847B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Vaporization, Distillation, Condensation, Sublimation, And Cold Traps (AREA)
  • Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
  • Feedback Control In General (AREA)

Abstract

本发明公开了一种面向常减压优化的稳态工况判别方法,通过选择最能表征常减压装置加工工况的参数进行稳态判断,稳态判断根据所判断的数据类型的不同,分为实时稳态判断和历史稳态判断两类:在实时任务中选择实时稳态判断,旨在对实时数据进行平稳段的鉴别,以判断是否可以进行装置的实时优化;在稳态操作知识库的建立时选择历史稳态判断,旨在对历史数据进行稳态查询和分析。判断依据指定时间范围内的数据均值和拟合直线的斜率,避免常规统计方法在常减压这类工况大范围变化情形下的误判问题。本发明尤其在加工原油切换造成生产调整较为频繁,装置瞬时波动较大时,能有效实现切油后的稳态工况及时判断,为装置实时优化奠定基础。

Description

一种面向常减压优化的稳态工况判别方法
技术领域
本发明涉及一种生产过程工况稳态判别方法,尤其涉及在常减压装置进行实时优化前,通过操作参数对常减压装置运行状态进行判断的方法,具体是一种面向常减压优化的稳态工况判别方法。
背景技术
常减压蒸馏装置是炼油加工的第一道工序,通过蒸馏可按产品生产方案将原油分割成相应的直馏汽油、煤油及各种润滑油馏分等半成品,半成品经适当精制和调配等,即可成为合格产品。常减压装置是炼油厂最大的耗能装置之一,约占炼油总耗能的25%~30%,因此通过装置实时优化,可以提高高价值产品收率,降低装置的能耗,对炼厂的经济效益提高具有极其重要的意义。
目前化工领域一般采用统计学方法和手段,选用均值、标准差、变异系数(包括极差变异系数、标准差变异系数)等统计量,以及直方图、控制图等统计技术,对生产运行进行统计稳态分析。一方面,由于抽样的局限性,样本信息不可能是总体信息的完整反应,故即使样本指数呈正态分布,上下控制界限为中心线±3倍标准差时,也会出现判断出错。另一方面,常减压装置是一个典型的多工况加工过程,在不同原油来料情况下,装置工况变化很大。类似常减压装置这样各参数波动范围大、常点较多的情况,若采用上述常规统计方法则会导致筛选过严,部分的稳态将被判定为非稳态,影响装置的实时优化,进而影响企业效益提升。
发明内容
本发明针对背景技术中存在的问题提出一种面向常减压优化的稳态工况判别方法,所述稳态工况判别方法通过选择最能表征常减压装置加工工况的参数进行稳态判断,稳态判断根据所判断的数据类型的不同,分为实时稳态判断和历史稳态判断两类:在实时任务中选择实时稳态判断,旨在对实时数据进行平稳段的鉴别,以判断是否可以进行装置的实时优化;在稳态操作知识库的建立时选择历史稳态判断,旨在对历史数据进行稳态查询和分析。
其中,实时稳态判断的步骤为:
1)以当前时间Ta作为基准,向历史方向前推时间ta,确定数据采集的起始时间Tb,即Tb=Ta-ta,以tb为时间间隔,采集该段时间范围内用以进行稳态判断的装置参数数据;
2)以tc为间隔对采集到的参数进行平滑,即对tc内每tb的取值求算术平均值;
3)依据稳态判断规则判断Tb至Ta时间范围内装置是否处于稳态,稳态判断规则:若为非稳态,则等待tc后,返回步骤1)继续判断;
若为稳态,则输出此稳态开始及结束时间,此次稳态判断结束,允许进行实时优化。
其中,历史稳态判断的步骤为:
1)确定进行稳态判断的历史数据的起始时间T1和结束时间T2
2)计算ΔT=T2-T1,判断ΔT是否大于等于ta
若ΔT<ta,则此次稳态判断结束;
若ΔT≥ta,则以T2作为首个稳态工况判断的结束时间,确定首个稳态工况判断的起始时间Tx=T2-ta,此时循环迭代标识符K=0;
3)采集Tx至T2时间范围内的稳态参数,并对每tc内的所有采集数据点取算术平均值,依据稳态判断规则,判断Tx至T2时间范围内,常减压装置生产工况是否为稳态:
若为非稳态,则转步骤5);
若为稳态,则将数据的开始时间Tx向历史方向前推tc,得到新的数据开始时间Tx=Tx-tc,K=K+1;
4)判断Tx是否在T1和T2时间范围内:
若Tx∈[T1,T2],则返回步骤3),以寻找最长稳态时间段;
否则转步骤5);
5)判断K是否为0:
若K=0,则将稳态判断的结束时间T2向历史方向前推tb,定义新的T2,即T2=T2-tb,返回步骤2);
若K≠0,则输出最长稳态时间段,计算下一次稳态判断的结束时间T2=T2-tb+tc,返回步骤2),对新工况进行稳态判断。
作为一种优选的实施方式,所述ta=3h,tb=5min,tc=30min。
作为一种优选的实施方式,所述稳态判断规则,具体是依据下式:
式中xi(1≤i≤n)为实时或历史采集的某一段时间范围内每个时间点ti的数值,为xi(1≤i≤n)的平均值,γ为该稳态参数对应的允许最大偏差,γ由工艺人员根据生产状况设定;若参与稳态判断的该段时间范围内的所有参数数值均在对应的允许偏差范围内,或者若该段时间范围内每个参数最多仅有一个点超出对应的允许偏差范围γ,但小于1.2γ,则对参数各点(ti,xi)进行线性拟合,若拟合出的直线斜率k∈[-0.1,0.1,]则认为该段时间范围内,常减压装置生产工况处于稳态。
作为一种优选的实施方式,针对参数各点(ti,xi)拟合出的直线为yi=k*ti+b,k为拟合直线的斜率,拟合的目标函数为
优选的,选取的参数为常减压装置的原油加工流量、常压塔顶温度、实际换热终温、常一中回流量、常二线抽出量5个参数。
有益效果:
本发明提出一种面向常减压优化的稳态工况判别方法,通过选择最能表征常减压装置加工工况的参数(优选是原油加工流量、常压塔顶温度、实际换热终温、常一中回流量、常二线抽出量5个操作参数)进行稳态分析,实现对常减压装置运行工况状态的准确判断,尤其在加工原油切换造成生产调整较为频繁,装置瞬时波动较大时,能有效实现切油后的稳态工况的及时判断,为装置实时优化奠定基础。除此之外,通过对常减压装置历史数据的稳态判断,建立稳态操作知识库,为基于相似工况查询的油种切换指导提供数据基础。
附图说明
图1为实时稳态判断流程图;
图2为历史稳态判断流程图。
具体实施案例
下面对本发明的实施例作详细说明,本实施例在以本发明技术方案为前提下进行实施,给出了详细的实施方式和具体的操作过程,但本发明的保护范围不限于下述的实施例。
以某企业常减压装置为例,该企业常减压装置具有炼化企业典型工艺,包括初馏塔、常压塔和减压塔。选取原油加工流量、常压塔顶温度、实际换热终温、常一中回流量、常二线抽出量这5个稳态参数的多个平稳段数据进行分析。根据该企业的具体生产情况,由工艺人员决定操作参数所允许的最大偏差γ分别为:原油加工流量30t/h、常压塔顶温度3.0℃、实际换热终温4.0℃、常一中回流量6.0t/h、常二线抽出量1.0t/h。
1.实时稳态判断
1)以当前时间Ta(2015-05-19 16:00)作为基准,向历史方向前推3h,确定数据采集的起始时间Tb(2015-05-19 13:00),以5min为时间间隔,采集该段时间范围内常减压装置的原油加工流量、常压塔顶温度、实际换热终温、常一中回流量、常二线抽出量5个参数的装置数据,如表1所示:
表1 实时稳态判断参数数据(平滑前)
2)以30min为间隔对采集到的稳态参数数据进行平滑,平滑后的各参数数据如表2所示:
表2 实时稳态参数数据(平滑后)
时间范围 原油加工流量 常压塔顶温度 实际换热终温 常一中回流量 常二线抽出量
2015/5/19 13:00:00~13:30:00 738.19 138.02 294.17 115.16 64.98
2015/5/19 13:30:00~14:00:00 741.91 138.20 292.90 115.14 64.96
2015/5/19 14:00:00~14:30:00 736.25 138.36 293.01 115.05 64.75
2015/5/19 14:30:00~15:00:00 734.31 138.47 294.51 114.97 64.69
2015/5/19 15:00:00~15:30:00 742.95 138.24 296.14 114.95 64.87
2015/5/19 15:30:00~16:00:00 742.72 137.78 295.66 115.06 64.02
3)判断在该段时间范围内(2015-05-19 13:00至2015-05-19 16:00)的原油加工流量、常压塔顶温度、实际换热终温、常一中回流量、常二线抽出量是否满足稳态要求,通过下列3个步骤实现。
(1)判断各参数的取值是否在允许偏差范围内,各参数的偏差如表3所示:
表3 各稳态参数数据(平滑后)与其均值的差值
依据表3中的各参数实际偏差,对比各参数允许的最大偏差可以看出,各参数的值均在对应的允许偏差范围内。
(2)在满足步骤(1)的基础上,对各参数数值(平滑前)进行线性拟合,拟合得到的各直线斜率如表4所示。
表4 各参数线性拟合直线斜率
拟合直线 原油加工流量偏差 常压塔顶温度偏差 实际换热终温偏差 常一中回流量偏差 常二线抽出量偏差
k 0.088 -0.004 -0.082 -0.006 -0.020
由于各参数的拟合曲线斜率k均满足k∈[-0.1,0.1],则认为其满足稳态判断步骤(2)的指标要求。
可见,依据本发明稳态判断条件,该段时间范围内(2015-05-19 13:00至2015-05-19 16:00)装置处于稳态。此时输出稳态起始时间和结束时间,此次稳态判断结束,后续可以进行装置的实时优化。
2.历史数据稳态判断
1)确定进行稳态判断的历史数据的起始时间T1(2015-05-19 12:30)和结束时间T2(2015-05-19 16:00);
2)计算ΔT,ΔT=T2-T1=3.5h,ΔT>3h;
3)以T2(2015-05-19 16:00)作为首个稳态工况判断的结束时间,计算得到首个稳态工况判断的起始时间Tx=T2-3h,即2015-05-19 13:00,此时K=0;
4)以5min为时间间隔,采集Tx(2015-05-19 13:00)至T2(2015-05-19 16:00)时间范围内的原油加工流量、常压塔顶温度、实际换热终温、常一中回流量和常二线抽出量,如表5所示。
表5 历史稳态参数数据(平滑前)
5)对每30min内所有采集数据点取平均值,并依据稳态判断规则判断,得到在Tx至T2时间范围内,常减压装置生产工况为稳态(详见实时数据稳态判断);
6)将数据的开始时间Tx(2015-05-19 13:00)向历史方向前推30min,计算得到新的数据开始时间Tx(2015-05-19 12:30),此时K=K+1=1;
7)判断新计算得到的Tx∈[T1,T2],采集Tx至T2时间范围内装置稳态参数,并以30min为时间间隔对数据采集点进行平滑,计算结果如表6所示:
表6 稳态参数数据(平滑后)
时间范围 原油加工流量 常压塔顶温度 实际换热终温 常一中回流量 常二线抽出量
2015/5/19 12:30:00~13:00:00 738.27 137.73 293.23 115.24 65.05
2015/5/19 13:00:00~13:30:00 738.19 138.02 294.17 115.16 64.98
2015/5/19 13:30:00~14:00:00 741.91 138.20 292.90 115.14 64.96
2015/5/19 14:00:00~14:30:00 736.25 138.36 293.01 115.05 64.75
2015/5/19 14:30:00~15:00:00 734.31 138.47 294.51 114.97 64.69
2015/5/19 15:00:00~15:30:00 742.95 138.24 296.14 114.95 64.87
2015/5/19 15:30:00~16:00:00 742.72 137.78 295.66 115.06 64.02
8)计算各参数与均值的偏差,如表7所示:
表7 各稳态参数的采集值与均值的偏差
由表7可见,各参数在各时间点的采集值均在对应的允许偏差范围内。
9)对各稳态参数在各时间点的采集值(平滑前)进行线性拟合,得到的各直线斜率如表8所示。
表8 稳态参数拟合直线斜率
拟合直线 原油加工流量偏差 常压塔顶温度偏差 实际换热终温偏差 常一中回流量偏差 常二线抽出量偏差
k 0.089 0.005 0.070 -0.006 -0.019
由于各参数指标的拟合曲线斜率均在[-0.1,0.1]范围内,满足稳态判断要求。即在时间范围(2015-05-19 12:30至2015-05-19 16:00)内装置处于稳态。
10)将数据的开始时间Tx(2015-05-19 12:30)向历史方向前推30min,计算得到新的数据开始时间Tx(2015-05-19 12:00),此时K=K+1=2;
11)判断判断K=2≠0,则输出最长稳态时间段(2015-05-19 12:30~16:00,并计算下一次稳态判断的历史数据结束时间T2=T2+25min(2015-05-19 12:25);
12)计算ΔT=T2-T1,ΔT<3h,此次稳态判断结束。
尽管本发明已经参照附图和优选实例进行了说明,但是,对于本领域的技术人员来说,本发明可以有各种更改和变化。本发明的各种更改、变化,和等同物有所附的权利要求书的内容涵盖。
本发明未涉及技术均与现有技术相同,或可采用现有技术实现。

Claims (5)

1.一种面向常减压优化的稳态工况判别方法,其特征在于所述稳态工况判别方法通过选择最能表征常减压装置加工工况的参数进行稳态判断,稳态判断根据所判断的数据类型的不同,分为实时稳态判断和历史稳态判断两类:在实时任务中选择实时稳态判断,旨在对实时数据进行平稳段的鉴别,以判断是否可以进行装置的实时优化;在稳态操作知识库的建立时选择历史稳态判断,旨在对历史数据进行稳态查询和分析;
实时稳态判断的步骤为:
1)以当前时间Ta作为基准,向历史方向前推时间ta,确定数据采集的起始时间Tb,即Tb=Ta-ta,以tb为时间间隔,采集该段时间范围内用以进行稳态判断的装置参数数据;
2)以tc为间隔对采集到的参数进行平滑,即对tc内每tb的取值求算术平均值;
3)依据稳态判断规则判断Tb至Ta时间范围内装置是否处于稳态:
若为非稳态,则等待tc后,返回步骤1)继续判断;
若为稳态,则输出此稳态开始及结束时间,此次稳态判断结束,允许进行实时优化;
历史稳态判断的步骤为:
1)确定进行稳态判断的历史数据的起始时间T1和结束时间T2
2)计算ΔT=T2-T1,判断ΔT是否大于等于ta
若ΔT<ta,则此次稳态判断结束;
若ΔT≥ta,则以T2作为首个稳态工况判断的结束时间,确定首个稳态工况判断的起始时间Tx=T2-ta,此时循环迭代标识符K=0;
3)采集Tx至T2时间范围内的稳态参数,并对每tc内的所有采集数据点取算术平均值,依据稳态判断规则,判断Tx至T2时间范围内,常减压装置生产工况是否为稳态:
若为非稳态,则转步骤5);
若为稳态,则将数据的开始时间Tx向历史方向前推tc,得到新的数据开始时间Tx=Tx-tc,K=K+1;
4)判断Tx是否在T1和T2时间范围内:
若Tx∈[T1,T2],则返回步骤3),以寻找最长稳态时间段;
否则转步骤5);
5)判断K是否为0:
若K=0,则将稳态判断的结束时间T2向历史方向前推tb,定义新的T2,即T2=T2-tb,返回步骤2);
若K≠0,则输出最长稳态时间段,计算下一次稳态判断的结束时间T2=T2-tb+tc
返回步骤2),对新工况进行稳态判断。
2.根据权利要求1所述的一种面向常减压优化的稳态工况判别方法,其特征在于所述ta=3h,tb=5min,tc=30min。
3.根据权利要求1所述的一种面向常减压优化的稳态工况判别方法,其特征在于所述稳态判断规则,具体是依据下式:
<mrow> <mo>-</mo> <mi>&amp;gamma;</mi> <mo>&amp;le;</mo> <msub> <mi>x</mi> <mi>i</mi> </msub> <mo>-</mo> <mover> <mi>x</mi> <mo>&amp;OverBar;</mo> </mover> <mo>&amp;le;</mo> <mi>&amp;gamma;</mi> </mrow>
式中xi(1≤i≤n)为实时或历史采集的某一段时间范围内每个时间点ti的数值,为xi(1≤i≤n)的平均值,γ为该稳态参数对应的允许最大偏差,γ由工艺人员根据生产状况设定;若参与稳态判断的该段时间范围内的所有参数数值均在对应的允许偏差范围内,或者若该段时间范围内每个参数最多仅有一个点超出对应的允许偏差范围γ,但小于1.2γ,则对参数各点(ti,xi)进行线性拟合,若拟合出的直线斜率k∈[-0.1,0.1],则认为该段时间范围内,常减压装置生产工况处于稳态。
4.根据权利要求3所述的一种面向常减压优化的稳态工况判别方法,其特征在于针对参数各点(ti,xi)拟合出的直线为yi=k*ti+b,k为拟合直线的斜率,拟合的目标函数为
5.根据权利要求1所述的一种面向常减压优化的稳态工况判别方法,其特征在于选取的参数为常减压装置的原油加工流量、常压塔顶温度、实际换热终温、常一中回流量、常二线抽出量5个参数。
CN201510379120.1A 2015-07-01 2015-07-01 一种面向常减压优化的稳态工况判别方法 Active CN104977847B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510379120.1A CN104977847B (zh) 2015-07-01 2015-07-01 一种面向常减压优化的稳态工况判别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510379120.1A CN104977847B (zh) 2015-07-01 2015-07-01 一种面向常减压优化的稳态工况判别方法

Publications (2)

Publication Number Publication Date
CN104977847A CN104977847A (zh) 2015-10-14
CN104977847B true CN104977847B (zh) 2018-05-22

Family

ID=54274456

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510379120.1A Active CN104977847B (zh) 2015-07-01 2015-07-01 一种面向常减压优化的稳态工况判别方法

Country Status (1)

Country Link
CN (1) CN104977847B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105389648A (zh) * 2015-10-21 2016-03-09 南京富岛信息工程有限公司 一种常减压装置稳态工况的判别方法
CN110261693A (zh) * 2019-04-16 2019-09-20 南京华盾电力信息安全测评有限公司 一种基于稳态采样技术的在线热力试验方法及***
CN111026031B (zh) * 2019-12-13 2023-01-31 红云红河烟草(集团)有限责任公司 一种卷烟制丝过程数据的稳态识别方法
CN111338310B (zh) * 2020-03-30 2023-02-07 南京富岛信息工程有限公司 一种工业过程稳态工况识别与分类方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102684201A (zh) * 2012-05-30 2012-09-19 华南理工大学 一种基于电压越限概率的含风电场电网无功优化方法
CN103838216A (zh) * 2014-03-07 2014-06-04 华北电力大学(保定) 基于数据驱动案例匹配的电站锅炉燃烧优化方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102684201A (zh) * 2012-05-30 2012-09-19 华南理工大学 一种基于电压越限概率的含风电场电网无功优化方法
CN103838216A (zh) * 2014-03-07 2014-06-04 华北电力大学(保定) 基于数据驱动案例匹配的电站锅炉燃烧优化方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
原油蒸馏装置在线优化软件的工程实施;侯凯锋等;《石油化工》;20001220;第29卷(第12期);第934-938页 *
基于实时运行数据的稳态在线优化***研究;余黎黎;《信息科技辑》;20021215;第1页倒数4-5行,第5页倒数1-3行,第6页,第7页第1-8行,第52页,第53页第1-7行,附图4-4 *
常减压蒸馏装置在线数据校正与优化控制;周传光等;《青岛化工学院学报》;20000630;第21卷(第2期);第134-138页 *
稳态过程在线数据校正技术的工业实施;李博等;《石油化工》;20001020;第29卷(第10期);第768-771页 *
过程运行数据的稳态检测方法综述;刘吉臻等;《仪器仪表学报》;20130831;第34卷(第8期);第1739-1748页 *

Also Published As

Publication number Publication date
CN104977847A (zh) 2015-10-14

Similar Documents

Publication Publication Date Title
CN104977847B (zh) 一种面向常减压优化的稳态工况判别方法
CN101863088B (zh) 一种橡胶混炼过程中门尼粘度的预报方法
CN105676817B (zh) 不同大小样本均值-标准偏差控制图的统计过程控制方法
CN106781486B (zh) 基于浮动车数据的交通状态评价方法
AU2017218149B2 (en) Methods and systems for detection of abnormal karyotypes
CN108171142B (zh) 一种确定复杂工业过程中关键变量因果关系的方法
CN104990893B (zh) 一种基于相似判别的汽油辛烷值检测方法
CN111143413A (zh) 基于数据流概念漂移的异常检测方法
CN104777785A (zh) 一种基于指令域分析的数控加工工艺参数动态优化方法
CN110298578B (zh) 一种混流生产线动态瓶颈预测方法
CN114742481B (zh) 一种基于成分的特钢性能评估方法及***
CN105242709A (zh) 一种模具温度自适应控制方法
CN105868164A (zh) 一种基于有监督的线性动态***模型的软测量建模方法
CN116433085A (zh) 一种轧制过程控制***性能评估方法
CN102457878B (zh) 一种在关键业务预测中剔除异常数据的方法和装置
CN115169453A (zh) 基于密度聚类协同深度残差网络的热连轧宽度预测方法
CN108537249B (zh) 一种密度峰值聚类的工业过程数据聚类方法
CN111625934B (zh) 一种基于d-s证据理论的退火加热过程的多模态识别方法
CN117587403A (zh) 基于化学转化的铝合金模板表面处理方法及***
CN108364095B (zh) 基于数据挖掘的炼钢生产过程钢水质量诊断方法
CN109033590B (zh) 一种基于模锻模具磨损的工艺参数优化方法
CN115202190B (zh) 一种基于工艺指标调节潜质的工况划分和识别方法
CN111145835A (zh) 一种多重熔解曲线的软件分析方法
CN115907348A (zh) 一种多晶硅产品质量检测和控制的方法、装置及设备
CN113822587B (zh) 一种基于总线电流数据的工厂产能评估方法

Legal Events

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