CN110245434A - Miha纯气动操作条件下氢气传输调控模型建模方法 - Google Patents
Miha纯气动操作条件下氢气传输调控模型建模方法 Download PDFInfo
- Publication number
- CN110245434A CN110245434A CN201910528046.3A CN201910528046A CN110245434A CN 110245434 A CN110245434 A CN 110245434A CN 201910528046 A CN201910528046 A CN 201910528046A CN 110245434 A CN110245434 A CN 110245434A
- Authority
- CN
- China
- Prior art keywords
- hydrogen
- residual oil
- reactor
- interface
- miha
- 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.)
- Granted
Links
- 239000001257 hydrogen Substances 0.000 title claims abstract description 87
- 229910052739 hydrogen Inorganic materials 0.000 title claims abstract description 87
- IEYFRQLMJACEQY-UHFFFAOYSA-N 16-iodo-3-methylhexadecanoic acid Chemical compound OC(=O)CC(C)CCCCCCCCCCCCCI IEYFRQLMJACEQY-UHFFFAOYSA-N 0.000 title claims abstract description 21
- 101100264174 Mus musculus Xiap gene Proteins 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 title claims abstract description 19
- 125000004435 hydrogen atom Chemical class [H]* 0.000 title 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 claims abstract description 71
- 230000005540 biological transmission Effects 0.000 claims abstract description 40
- 150000002431 hydrogen Chemical class 0.000 claims abstract description 25
- 238000013461 design Methods 0.000 claims abstract description 12
- 239000007788 liquid Substances 0.000 claims description 14
- 238000005984 hydrogenation reaction Methods 0.000 claims description 13
- 239000012071 phase Substances 0.000 claims description 13
- 238000006243 chemical reaction Methods 0.000 claims description 11
- 239000007791 liquid phase Substances 0.000 claims description 11
- 238000004458 analytical method Methods 0.000 claims description 7
- 239000007789 gas Substances 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 6
- 238000009792 diffusion process Methods 0.000 claims description 3
- 230000004907 flux Effects 0.000 claims description 3
- 239000012528 membrane Substances 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 230000005501 phase interface Effects 0.000 claims 1
- 230000000704 physical effect Effects 0.000 abstract description 5
- 230000003321 amplification Effects 0.000 abstract description 3
- 238000003199 nucleic acid amplification method Methods 0.000 abstract description 3
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Production Of Liquid Hydrocarbon Mixture For Refining Petroleum (AREA)
Abstract
本发明涉及MIHA纯气动操作条件下氢气传输调控模型建模方法，考虑了两种输氢机制对氢气传输量的影响，建立的模型综合反映了反应器结构、体系物性、操作参数等对氢气传输的影响，可用于预测氢气传输量随反应器设计参数的变化，有利于MIHA反应器的设计和放大。
Description
技术领域
本发明属于反应器、建模技术领域，具体涉及MIHA纯气动操作条件下氢气传输调控模型建模方法。
背景技术
氢气泡在渣油体系停留时间内，氢气传输量是决定渣油反应效率的关键参数。氢气在渣油中的传输一般认为由亨利定律和菲克定律共同决定，但渣油加氢反应器微界面体系中氢气的传输过程尚未见报导，因此需建立分析渣油加氢反应器内的氢气传输模型，从而实现对MIHA反应器设计及MIHA的反应体系设计的指导。
发明内容
本发明的目的在于提供MIHA纯气动操作条件下氢气传输调控模型建模方法，以研究反应器结构、体系物性以及操作参数对氢气传输的影响，从而实现对MIHA反应器设计及MIHA的反应体系设计的指导。
MIHA微气泡形成可采用三种方式，即：纯液动、纯气动以及气液联动。纯液动和纯气动操作条件下，体系运行以及微气泡形成所需能量完全由液体机械能或气体静压能提供；气液联动操作条件下，气体静压能和液体机械能同时提供体系运行及微气泡形成所需能量。本发明探讨了纯气动操作条件下氢气传输调控模型建模方法。
为实现上述目的，本发明采用如下技术方案：
一种MIHA纯气动操作条件下氢气传输调控模型建模方法，包括如下步骤：
(1)分析MIHA纯气动操作条件下渣油加氢反应体系中氢气传输的影响参数；
假定在渣油加氢反应体系中，渣油的主要成分为含硫沥青质，其反应速率和转化率由体系液相中氢分子总传输量NT决定；基于压力传输和界面传输两种机制确定氢气在反应器内的质量传递；
NT＝Na+Np(3)
其中Na为界面传输量，mol；Np为压力传输量，mol；
(2)获取基于界面传输的氢气质量传递；
假定体系中氢气气泡大小均一，且通过相界面的传输速率相同，那么，通过气液相界面所传输的氢分子总量Na由体系中稳定存在的气泡总数量n0、及渣油在反应器内的平均停留时间te所决定，如方程(2)所示：
通过对渣油在反应器内的运动过程分析，渣油平均停留时间te表示为：
其中，VL为渣油液相主体体积，m3；QL为新鲜渣油流量，m3/s；
VL由氢气-渣油两相相对运动决定，通过对气液两相运动的分析后得：
H0为反应器高度，m；S0为反应器横截面积，m2；v32为气泡上升速度，m/s；
将式(4)代入式(3)得：
体系中稳定存在的气泡总数量n0由气体体积守恒及气泡运动分析得到，即：
其中QG为氢气流量，m3/s；d32为气泡Sauter平均直径，m；v32为气泡上升速度，m/s；
相界面的传输速率基于菲克定律获取，即：
其中Jz＝0指气泡液膜内侧边界处氢气传质通量，mol/m2；DL为氢气在渣油中的液相扩散系数，m2/s；
CH *为气液界面处氢气饱和浓度，mol/m3；
C* H＝Pm/HA (8)
其中Pm为操作压力，Pa；HA为亨利系数，Pa.m3/mol；
将方程(5)～(8)代入方程(2)后得：
(3)获取基于压力传输的氢气质量传递；
压力输氢量NP及其极限值均由Henry定律决定，与渣油液相主体体积VL、操作压力Pm及渣油平均停留时间te有关，假定NP为极限溶解量，即：
Np＝VL×CH * (10)
当渣油的摩尔质量为M1kg/mol时，CH *与渣油中氢气的摩尔分数x有如下关系：
其中，ρL为渣油密度，kg/m3；
由方程(4)、(9)、(11)得：
(4)获取氢分子总传输量；
由方程(1)(9)(12)得：
本发明的另一目的在于提供上述方法构建的MIHA纯气动操作条件下氢气传输调控模型。
本发明的又一目的在于提供上述方法设计的反应器。
本发明的反应器结构可参见发明人在先申请的专利CN106187660A，本发明中不再赘述。本发明中利用构建的模型反应器结构、体系物性以及操作参数对氢气传输的影响，从而可基于氢气传输效果进行相关的反应器结构参数设计。
本发明的方法针对MIHA建立了纯气动操作条件下氢气传输调控模型，考虑了两种输氢机制对氢气传输量的影响，建立的模型综合反映了反应器结构、体系物性、操作参数等对氢气传输的影响，可用于预测氢气传输量随反应器设计参数的变化，有利于MIHA反应器的设计和放大。
附图说明
图1是纯气动条件下操作压力Pm对单一界面传输量Na的影响；
图2是纯气动条件下气泡大小d32对单一界面传输量Na的影响；
图3是纯气动条件下操作压力Pm对单一压力输氢量NP的影响；
图4是纯气动条件下气泡大小d32对单一压力输氢量NP的影响；
图5是纯气动条件下输氢总量NT随操作压力Pm的变化关系；
图6是纯气动条件下输氢总量NT随气泡大小d32的变化关系。
具体实施方式
下面结合附图说明和具体实施方式对本发明的技术方案作进一步阐述。
实施例1
一种MIHA纯气动操作条件下氢气传输调控模型建模方法，包括如下步骤：
(1)分析MIHA纯气动操作条件下渣油加氢反应体系中氢气传输的影响参数；
假定在渣油加氢反应体系中，渣油的主要成分为含硫沥青质，其反应速率和转化率由体系液相中氢分子总传输量NT决定；基于压力传输和界面传输两种机制确定氢气在反应器内的质量传递；
NT＝Na+Np (5)
其中Na为界面传输量，mol；Np为压力传输量，mol；
(2)获取基于界面传输的氢气质量传递；
假定体系中氢气气泡大小均一，且通过相界面的传输速率相同，那么，通过气液相界面所传输的氢分子总量Na由体系中稳定存在的气泡总数量n0、及渣油在反应器内的平均停留时间te所决定，如方程(2)所示：
通过对渣油在反应器内的运动过程分析，渣油平均停留时间te表示为：
其中，VL为渣油液相主体体积，m3；QL为新鲜渣油流量，m3/s；
VL由氢气-渣油两相相对运动决定，通过对气液两相运动的分析后得：
H0为反应器高度，m；S0为反应器横截面积，m2；v32为气泡上升速度，m/s；
将式(4)代入式(3)得：
体系中稳定存在的气泡总数量n0由气体体积守恒及气泡运动分析得到，即：
其中QG为氢气流量，m3/s；d32为气泡Sauter平均直径，m；
相界面的传输速率基于菲克定律获取，即：
其中Jz＝0指气泡液膜内侧边界处氢气传质通量，mol/m2；DL为氢气在渣油中的液相扩散系数，m2/s；
CH*为气液界面处氢气饱和浓度，mol/m3；
其中Pm为操作压力，Pa；HA为亨利系数，Pa.m3/mol；
将方程(5)～(8)代入方程(2)后得：
(3)获取基于压力传输的氢气质量传递；
压力输氢量NP及其极限值均由Henry定律决定，与渣油液相主体体积VL、操作压力Pm及渣油平均停留时间te有关，假定NP为极限溶解量，即：
Np＝VL×CH * (10)
当渣油的摩尔质量为M1kg/mol时，CH *与渣油中氢气的摩尔分数x有如下关系：
其中，ρL为渣油密度，kg/m3；
由方程(4)、(9)、(11)得：
(4)获取氢分子总传输量；
由方程(1)(9)(12)得：
实施例2
本实施例基于实施例1的建模方法，针对具体的反应器结构和反应体系研究纯气动条件下操作条件、气泡尺度对氢气传输的影响。反应器内气泡尺度受反应器结构、物性参数、操作条件影响，其在纯气动操作条件下的调控模型可参考申请人的在先专利CN109684769A。
计算条件如下：
反应器设计参数：
反应器高度H0＝4.0m；反应器横截面积S0＝0.0013m2；
渣油密度ρL＝800Kg/m3；渣油的摩尔质量M1＝0.416kg/mol；
操作压力Pm＝0-25MPa；
新鲜渣油流量QL＝5.83×10-7m3/s；氢气流量QG＝5.56×10-4m3/s；
结果如图1～6所示。
图1和图2反映了渣油加氢体系中气泡尺度和操作压力对界面输氢量的影响。
由图1可看出，当气泡直径为毫米级时(如1.0mm和3.0mm)，界面输氢量随操作压力的增大而增加的趋势不显著；但当气泡直径为微米级(小于1.0mm)时，该趋势随气泡尺度的减小而愈加显著。图2亦体现了该变化规律。由图1还可知，当操作压力较小(如小于3.0MPa)时，气泡尺度对界面输氢量的影响较为复杂。
图3和图4反映了渣油加氢体系中气泡尺度和操作压力对压力输氢量的影响。
由图3可知，当体系中气泡直径为0.05mm时，操作压力由常压增至25MPa时，压力输氢量增加值约1.0mol，此时增大操作压力对于体系氢传输量的增加几乎可以忽略；但是，当气泡直径增大至3.0mm时，压力输氢量增加约9.0mol，增幅是常压时数十倍。由此可见，对于微米级气泡，升高操作压力对压力输氢量的影响不大。
由4图可知：当操作压力一定时，随着体系中气泡尺度的增大，压力输氢量呈三个不同的变化趋势：当气泡尺度约小于0.10mm或大于0.30mm时，气泡尺度对压力输氢量的影响较小；当气泡尺度介于0.10和0.30mm之间时，压力输氢量随气泡尺度的增大呈快速增加的趋势。当体系中气泡尺度一定时，压力输氢量随操作压力的增大而增大。
图5体现当渣油加氢体系中气泡尺度不同时，操作压力对体系输氢总量的影响。由此图可知：当体系中气泡尺度一定时，体系输氢总量随操作压力的增大而增加；当操作压力一定时，体系输氢总量随气泡尺度的增大而减小，该趋势当操作压力大于5.0MPa时是确定的，但当操作压力小于5.0MPa时，该趋势则较为复杂。图6也体现了这一变化趋势。由图5和图6还可知，当体系中气泡分别为微米级和毫米级时，前者输氢总量远大于后者。
由图1～6的结果可知，当渣油加氢体系中气泡尺度为毫米级时，升高操作压力有利于界面输氢和压力输氢量的增加，从而使体系氢传输总量增加；当渣油加氢体系中气泡尺度为微米级时，界面输氢量决定了体系中氢传输总量，其大小受操作压力的影响较小，也即，若在能形成微界面体系，则体系反应所需氢的量可以在较低操作压力下由界面输氢过程提供。
Claims (3)
1.一种MIHA纯气动操作条件下氢气传输调控模型建模方法，其特征在于，包括如下步骤：
(1)分析MIHA纯气动操作条件下渣油加氢反应体系中氢气传输的影响参数；
假定在渣油加氢反应体系中，渣油的主要成分为含硫沥青质，其反应速率和转化率由体系液相中氢分子总传输量NT决定；基于压力传输和界面传输两种机制确定氢气在反应器内的质量传递；
NT＝Na+Np (1)
其中Na为界面传输量，mol；Np为压力传输量，mol；
(2)获取基于界面传输的氢气质量传递；
假定体系中氢气气泡大小均一，且通过相界面的传输速率相同，那么，通过气液相界面所传输的氢分子总量Na由体系中稳定存在的气泡总数量n0、及渣油在反应器内的平均停留时间te所决定，如方程(2)所示：
通过对渣油在反应器内的运动过程分析，渣油平均停留时间te表示为：
其中，VL为渣油液相主体体积，m3；QL为新鲜渣油流量，m3/s；
VL由氢气-渣油两相相对运动决定，通过对气液两相运动的分析后得：
H0为反应器高度，m；S0为反应器横截面积，m2；v32为气泡上升速度，m/s；
将式(4)代入式(3)得：
体系中稳定存在的气泡总数量n0由气体体积守恒及气泡运动分析得到，即：
其中QG为氢气流量，m3/s；d32为气泡Sauter平均直径，m；
相界面的传输速率基于菲克定律获取，即：
其中Jz＝0指气泡液膜内侧边界处氢气传质通量，mol/m2；DL为氢气在渣油中的液相扩散系数，m2/s；
CH *为气液界面处氢气饱和浓度，mol/m3；
其中Pm为操作压力，Pa；HA为亨利系数，Pa.m3/mol；
将方程(5)～(8)代入方程(2)后得：
(3)获取基于压力传输的氢气质量传递；
压力输氢量NP及其极限值均由Henry定律决定，与渣油液相主体体积VL、操作压力Pm及渣油平均停留时间te有关，假定NP为极限溶解量，即：
Np＝VL×CH * (10)
当渣油的摩尔质量为M1 kg/mol时，CH *与渣油中氢气的摩尔分数x有如下关系：
其中，ρL为渣油密度，kg/m3；
由方程(4)、(9)、(11)得：
(4)获取氢分子总传输量；
由方程(1)(9)(12)得：
2.权利要求1所述方法构建的MIHA纯气动操作条件下氢气传输调控模型。
3.权利要求1所述方法设计的反应器。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910528046.3A CN110245434B (zh) | 2019-06-18 | 2019-06-18 | Miha纯气动操作条件下氢气传输调控模型建模方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910528046.3A CN110245434B (zh) | 2019-06-18 | 2019-06-18 | Miha纯气动操作条件下氢气传输调控模型建模方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110245434A true CN110245434A (zh) | 2019-09-17 |
| CN110245434B CN110245434B (zh) | 2023-06-27 |
Family
ID=67887885
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910528046.3A Active CN110245434B (zh) | 2019-06-18 | 2019-06-18 | Miha纯气动操作条件下氢气传输调控模型建模方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110245434B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113035285A (zh) * | 2021-04-01 | 2021-06-25 | 南京延长反应技术研究院有限公司 | 微气泡尺寸对油品加氢脱硫效果影响的计算方法 |
| CN113075091A (zh) * | 2021-02-24 | 2021-07-06 | 南京延长反应技术研究院有限公司 | 一种确定px氧化制ta体系中气泡大小的方法 |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103745088A (zh) * | 2013-12-12 | 2014-04-23 | 西北大学 | 一种煤焦油加氢氢耗的动力学计算方法 |
| CN105975685A (zh) * | 2016-05-03 | 2016-09-28 | 华东理工大学 | 一种用于渣油延迟焦化过程的建模和优化方法 |
| CA2931753A1 (fr) * | 2015-06-04 | 2016-12-04 | IFP Energies Nouvelles | Procede d'exploitation d'un gisement d'hydrocarbures contenant des composes organo-soufres au moyen d'un modele thermo-cinetique et d'une simulation de reservoir compositionnelle |
| CN106202910A (zh) * | 2016-07-07 | 2016-12-07 | 华东理工大学 | 一种渣油加氢处理装置的收率实时预测方法及其应用 |
| CN107346378A (zh) * | 2017-08-30 | 2017-11-14 | 南京大学 | 微界面强化反应器传质速率构效调控模型建模方法 |
| US20170372196A1 (en) * | 2016-06-22 | 2017-12-28 | Saudi Arabian Oil Company | Systems and methods for rapid prediction of hydrogen-induced cracking (hic) in pipelines, pressure vessels, and piping systems and for taking action in relation thereto |
| CN109684769A (zh) * | 2019-01-09 | 2019-04-26 | 南京大学 | Miha纯气动操作条件下气泡尺度调控模型建模方法 |
| CN109887551A (zh) * | 2019-02-04 | 2019-06-14 | 南京大学 | Miha纯气动操作条件下传质调控模型建模方法 |
-
2019
- 2019-06-18 CN CN201910528046.3A patent/CN110245434B/zh active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103745088A (zh) * | 2013-12-12 | 2014-04-23 | 西北大学 | 一种煤焦油加氢氢耗的动力学计算方法 |
| CA2931753A1 (fr) * | 2015-06-04 | 2016-12-04 | IFP Energies Nouvelles | Procede d'exploitation d'un gisement d'hydrocarbures contenant des composes organo-soufres au moyen d'un modele thermo-cinetique et d'une simulation de reservoir compositionnelle |
| CN105975685A (zh) * | 2016-05-03 | 2016-09-28 | 华东理工大学 | 一种用于渣油延迟焦化过程的建模和优化方法 |
| US20170372196A1 (en) * | 2016-06-22 | 2017-12-28 | Saudi Arabian Oil Company | Systems and methods for rapid prediction of hydrogen-induced cracking (hic) in pipelines, pressure vessels, and piping systems and for taking action in relation thereto |
| CN106202910A (zh) * | 2016-07-07 | 2016-12-07 | 华东理工大学 | 一种渣油加氢处理装置的收率实时预测方法及其应用 |
| CN107346378A (zh) * | 2017-08-30 | 2017-11-14 | 南京大学 | 微界面强化反应器传质速率构效调控模型建模方法 |
| CN109684769A (zh) * | 2019-01-09 | 2019-04-26 | 南京大学 | Miha纯气动操作条件下气泡尺度调控模型建模方法 |
| CN109887551A (zh) * | 2019-02-04 | 2019-06-14 | 南京大学 | Miha纯气动操作条件下传质调控模型建模方法 |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113075091A (zh) * | 2021-02-24 | 2021-07-06 | 南京延长反应技术研究院有限公司 | 一种确定px氧化制ta体系中气泡大小的方法 |
| CN113075091B (zh) * | 2021-02-24 | 2023-01-10 | 南京延长反应技术研究院有限公司 | 一种确定px氧化制ta体系中气泡大小的方法 |
| CN113035285A (zh) * | 2021-04-01 | 2021-06-25 | 南京延长反应技术研究院有限公司 | 微气泡尺寸对油品加氢脱硫效果影响的计算方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110245434B (zh) | 2023-06-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110245434A (zh) | Miha纯气动操作条件下氢气传输调控模型建模方法 | |
| CN206285902U (zh) | 使浆料匀速流出搅拌罐的装置 | |
| CN104929087B (zh) | 一种阻尼阻抗式调压装置 | |
| CN106869918A (zh) | 海上油井产能测试实时调整方法 | |
| CN110385094A (zh) | 一种液体水合物快速连续化制备装置 | |
| CN106933262B (zh) | 一种在线压力信号解析的多相流场量测及气-液传质增效控制方法 | |
| CN106278836A (zh) | 中等浓度甲醛与甲缩醛合成聚甲氧基二甲醚的装置和方法 | |
| CN102222168A (zh) | 一种深水钻井立管参激横向振动分析方法 | |
| CN103896977A (zh) | 一种生产乙烯基烷氧基硅烷的方法及其设备 | |
| CN106370410B (zh) | 一种液膜发生装置及其使用方法 | |
| CN204385208U (zh) | 一种用于生物反应器的摇动平台 | |
| CN201184184Y (zh) | 油井油气水连续分测装置 | |
| CN214681643U (zh) | 换热介质均布的列管式有机硅流化床反应器 | |
| CN201728110U (zh) | 用于生产硼氢化钠的硼氢化反应釜 | |
| CN104357324A (zh) | 一种用于生物反应器的摇动平台 | |
| CN207067098U (zh) | 泡排剂性能的加压评价装置 | |
| CN201212393Y (zh) | 一种称重式油井计量装置 | |
| CN208097739U (zh) | 重氮反应废气水封分离箱 | |
| CN103263884B (zh) | 一种合成反应器和氯胺法合成甲基肼装置和方法 | |
| CN207483756U (zh) | 一种平行反应器的罐压控制装置 | |
| CN207271221U (zh) | 连续进料气液反应装置 | |
| CN201154945Y (zh) | 一种用于甲烷液相氧化反应的鼓泡塔 | |
| CN205635431U (zh) | 一种羟基钴胺素连续制备装置 | |
| CN206778389U (zh) | 一种3‑乙硫基丁醛的生产装置 | |
| CN205903875U (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 |