CN106547968B - 吸收塔近壁逃逸现象量化方法 - Google Patents
吸收塔近壁逃逸现象量化方法 Download PDFInfo
- Publication number
- CN106547968B CN106547968B CN201610945568.XA CN201610945568A CN106547968B CN 106547968 B CN106547968 B CN 106547968B CN 201610945568 A CN201610945568 A CN 201610945568A CN 106547968 B CN106547968 B CN 106547968B
- Authority
- CN
- China
- Prior art keywords
- wall
- flue gas
- flow velocity
- absorption tower
- quantifying
- 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
Links
- 238000010521 absorption reaction Methods 0.000 title claims abstract description 32
- 238000000034 method Methods 0.000 title claims abstract description 29
- 239000007921 spray Substances 0.000 claims abstract description 36
- UGFAIRIUMAVXCW-UHFFFAOYSA-N Carbon monoxide Chemical compound [O+]#[C-] UGFAIRIUMAVXCW-UHFFFAOYSA-N 0.000 claims abstract description 32
- 239000003546 flue gas Substances 0.000 claims abstract description 32
- 238000005507 spraying Methods 0.000 claims abstract description 29
- 238000004088 simulation Methods 0.000 claims description 11
- 239000000779 smoke Substances 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 abstract description 2
- 101100356682 Caenorhabditis elegans rho-1 gene Proteins 0.000 abstract 2
- 238000004364 calculation method Methods 0.000 description 6
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 239000012530 fluid Substances 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 230000001419 dependent effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- FKLFBQCQQYDUAM-UHFFFAOYSA-N fenpiclonil Chemical compound ClC1=CC=CC(C=2C(=CNC=2)C#N)=C1Cl FKLFBQCQQYDUAM-UHFFFAOYSA-N 0.000 description 2
- 238000007689 inspection Methods 0.000 description 2
- 230000010349 pulsation Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000002457 bidirectional effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 229910052602 gypsum Inorganic materials 0.000 description 1
- 239000010440 gypsum Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 230000021715 photosynthesis, light harvesting Effects 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 230000001681 protective effect Effects 0.000 description 1
- 230000000541 pulsatile effect Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000005514 two-phase flow Effects 0.000 description 1
Images
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
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
- B01D53/34—Chemical or biological purification of waste gases
- B01D53/46—Removing components of defined structure
- B01D53/48—Sulfur compounds
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D53/00—Separation of gases or vapours; Recovering vapours of volatile solvents from gases; Chemical or biological purification of waste gases, e.g. engine exhaust gases, smoke, fumes, flue gases, aerosols
- B01D53/34—Chemical or biological purification of waste gases
- B01D53/74—General processes for purification of waste gases; Apparatus or devices specially adapted therefor
- B01D53/80—Semi-solid phase processes, i.e. by using slurries
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D2251/00—Reactants
- B01D2251/40—Alkaline earth metal or magnesium compounds
- B01D2251/404—Alkaline earth metal or magnesium compounds of calcium
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D2258/00—Sources of waste gases
- B01D2258/02—Other waste gases
- B01D2258/0283—Flue gases
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B01—PHYSICAL OR CHEMICAL PROCESSES OR APPARATUS IN GENERAL
- B01D—SEPARATION
- B01D2259/00—Type of treatment
- B01D2259/12—Methods and means for introducing reactants
- B01D2259/126—Semi-solid reactants, e.g. slurries
Abstract
本发明涉及电力生产，公开了一种吸收塔近壁逃逸现象量化方法，该方法包括：计算喷淋截面的烟气平均流速v1和喷淋密度ρ1；计算近壁区域的烟气平均流速v2和喷淋密度ρ2；根据喷淋截面的烟气平均流速v1、喷淋截面的喷淋密度ρ1、近壁区域的烟气平均流速v2以及近壁区域的喷淋密度ρ2计算近壁逃逸系数A。本发明可以对近壁区域逃逸现象进行量化计算，从而有效衡量吸收塔近壁逃逸现象的程度，以便在改造中采取相应的技术措施。
Description
技术领域
本发明涉及电力生产，具体地，涉及一种吸收塔近壁逃逸现象量化方法。
背景技术
目前在火电厂脱硫主要采用石灰石-石膏湿法脱硫，在湿法脱硫中，最常见的技术为喷淋空塔，吸收塔喷淋层设计时，多采用喷嘴均匀布置的方法，这种布置，在喷淋截面的主体部分，能达到均匀喷淋，但在靠近吸收塔塔壁的区域，由于喷淋覆盖率较低，区域阻力较小，烟气流速较快，形成局部区域的“烟气走廊”，此部分烟气的脱硫效果较差，从而影响整体的脱硫效率。
早期的喷淋空塔湿法脱硫技术，在设计时大多未考虑烟气的近壁逃逸问题，由于目前环保标准的日益提高，大部分喷淋空塔湿法脱硫需要进行改造，为了达到较高的脱硫效率，必须减小或消除吸收塔近壁逃逸问题，尤其是要达到超低排放的燃煤机组。这就需要有一个量化的指标来直观的衡量吸收塔近壁逃逸的程度，从而在改造过程中采取相应的技术措施。
发明内容
本发明的目的是提供一种吸收塔近壁逃逸现象量化方法，该吸收塔近壁逃逸现象量化方法可以对近壁区域逃逸现象进行量化计算，从而有效衡量吸收塔近壁逃逸现象的程度，以便在改造中采取相应的技术措施。
为了实现上述目的，本发明提供一种吸收塔近壁逃逸现象量化方法，该方法包括：计算喷淋截面的烟气平均流速v1和喷淋密度ρ1；计算近壁区域的烟气平均流速v2和喷淋密度ρ2；根据喷淋截面的烟气平均流速v1、喷淋截面的喷淋密度ρ1、近壁区域的烟气平均流速v2以及近壁区域的喷淋密度ρ2计算近壁逃逸系数A。
通过上述技术方案，采用本发明提供的吸收塔近壁逃逸现象量化方法，计算出喷淋截面的烟气平均流速v1、喷淋截面的喷淋密度ρ1、近壁区域的烟气平均流速v2和近壁区域的喷淋密度ρ2，并根据v1、ρ1、v2和ρ2计算近壁逃逸系数A。本发明可以对近壁区域逃逸现象进行量化计算，从而有效衡量吸收塔近壁逃逸现象的程度，以便在改造中采取相应的技术措施。
本发明的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明，但并不构成对本发明的限制。在附图中：
图1是本发明提供的吸收塔近壁逃逸现象量化方法的流程图；
图2是本发明提供的近壁区域选取示意图。
具体实施方式
以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明，并不用于限制本发明。
图1是本发明提供的吸收塔近壁逃逸现象量化方法的流程图。如图1所示，本发明提供一种吸收塔近壁逃逸现象量化方法，该方法包括：计算喷淋截面的烟气平均流速v1和喷淋密度ρ1；计算近壁区域的烟气平均流速v2和喷淋密度ρ2；根据喷淋截面的烟气平均流速v1、喷淋截面的喷淋密度ρ1、近壁区域的烟气平均流速v2以及近壁区域的喷淋密度ρ2计算近壁逃逸系数A。
喷淋截面设置有多个喷嘴，但是不同吸收塔喷淋截面上设置的喷嘴数量和喷嘴距离可能均不相同，因此可能导致喷淋截面的喷淋密度ρ1和近壁区域的喷淋密度ρ2均不相同。所述喷淋截面的喷淋密度ρ1和所述近壁区域的喷淋密度ρ2根据喷嘴的配置情况计算。所述喷淋截面的烟气平均流速v1和所述近壁区域的烟气平均流速v2和入口烟气量有关。
所述喷淋截面的烟气平均流速v1以及近壁区域的烟气平均流速v2通过数值模拟计算。
数值模拟也叫计算机模拟。依靠电子计算机，结合有限元或有限容积的概念，通过数值计算和图像显示的方法，达到对工程问题和物理问题乃至自然界各类问题研究的目的。
具体地，本发明利用CFD商用软件FLUENT进行数值模拟。
计算流体动力学(Computational Fluid Dynamics简称CFD)是通过计算机数值计算和图像显示，对包含有流体流动和传热等相关物理现象的系统所做的分析。CFD的基本思想可以归纳为：把原来在时间域及空间域上连续的物理量场，如速度场和压力场，用一系列有限离散点上的变量值的集合来代替，通过一定的原则和方式建立起关于这些离散点上变量之间关系的代数方程组，然后求解代数方程组获得场变量的近似值。
FLUENT软件的求解主要包括以下几个步骤：
1.确定几何形状，生成计算网格(用GAMBIT，也可以读入其它指定程序生成的网格)
2.输入网格、检查网格
3.选择求解器
4.选择求解的方程：层流或湍流(或无粘流)，化学组分或化学反应，传热模型等。确定其它需要的模型如：风扇、热交换器、多孔介质等模型。
5.确定流体物性、指定边界条件
6.条件计算控制参数、流场初始化
7.迭代计算
8.检查结果保存结果，后处理等。
本实施例模拟对象的网格采用GAMBIT进行划分。所做计算为吸收塔内烟气与喷淋的气液两相流动，因此计算域选取为吸收塔浆液上方烟气流经通道。整体模型均采用了四、六面体网格，共划分2416293个体网格。
本实施例数值模拟边界条件设置如下：从增压风机流出烟气的总体积流量为1245400Nm3/h，速度为10～12m/s(模拟时取10m/s)，密度为0.925kg/m3，进吸收塔设计温度为110℃。吸收塔内配有四层喷淋层，第四喷淋层在原第一喷淋层下方，有90个喷嘴，每个喷嘴的流量为83.3m3/h。其中靠壁面的29个喷嘴为单向实心锥喷嘴，其余61个为上下双向中空锥喷嘴。
本实施例选用统关模拟模型，并应用可实现的k-ε双方程模型进行模拟计算。k-ε模型是双方程模型中的最为广泛应用的湍流模型，它引入了表示各向同性小尺度涡机械能转化成热能速率的脉动动能耗散率ε。其k-ε双方程主要的特点就是它是通过求偏微分方程来考虑湍流物理量的输运过程，即通过求解偏微分方程确定脉动特征速度与平均场速度梯度的关系，而不是直接将两者联系起来；另外，模型中的特征长度不是由经验确定，而是以耗散尺度作为特征长度，并由求解相应的偏微分方程得到。相对来说可实现的k-ε双方程更为简单，使用方便，目前用于计算有旋的均匀剪切流、平面混合流、平面射流、圆形射流等流动，而且都取得了与试验数据比较一致的结果，适合工程问题的研究。
直角坐标系下，稳态的标准k-ε模型的通用控制方程如下：
在得到v1、ρ1、v2以及ρ2之后，所述近壁逃逸系数A可以通过以下公式得出：
其中，v2大于v1，ρ1大于ρ2。
所述近壁逃逸系数A越大，近壁逃逸现象越严重；所述近壁逃逸系数A越小，近壁逃逸现象越轻微，A为0时无近壁逃逸现象。
以上是最佳量化公式，另外，本发明还提供以下公式的变形也同样可以得到近壁逃逸系数A：
其中，v2大于v1，ρ1大于ρ2，近壁逃逸系数A越大，近壁逃逸现象越严重；近壁逃逸系数A越小，近壁逃逸现象越轻微，A为1时无近壁逃逸现象。
其中，v2大于v1，近壁逃逸系数A越大，近壁逃逸现象越严重；近壁逃逸系数A越小，近壁逃逸现象越轻微，A为1时无近壁逃逸现象。
其中，ρ1大于ρ2，近壁逃逸系数A越大，近壁逃逸现象越严重；近壁逃逸系数A越小，近壁逃逸现象越轻微，A为1时无近壁逃逸现象。
图2是本发明提供的近壁区域选取示意图。如图2所示，吸收塔喷淋截面上，D为吸收塔内径，d为圆环形近壁区域的宽度，r为喷嘴的喷射圆半径。其中圆环的宽度d的选取根据吸收塔喷淋层布置情况确定，本发明优选地，近壁区域圆环的宽度d为最外层喷嘴喷射圆半径的0.5～1倍(即r的0.5～1倍)。
通过上述技术方案，采用本发明提供的吸收塔近壁逃逸现象量化方法，计算出喷淋截面的烟气平均流速v1、喷淋截面的喷淋密度ρ1、近壁区域的烟气平均流速v2和近壁区域的喷淋密度ρ2，并根据v1、ρ1、v2和ρ2计算近壁逃逸系数A。本发明可以对近壁区域逃逸现象进行量化计算，从而有效衡量吸收塔近壁逃逸现象的程度，以便在改造中采取相应的技术措施。
以上结合附图详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种简单变型，这些简单变型均属于本发明的保护范围。
另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合，为了避免不必要的重复，本发明对各种可能的组合方式不再另行说明。
此外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明的思想，其同样应当视为本发明所公开的内容。
Claims (5)
1.一种吸收塔近壁逃逸现象量化方法，其特征在于，该方法包括：
计算喷淋截面的烟气平均流速v1和喷淋密度ρ1；
计算近壁区域的烟气平均流速v2和喷淋密度ρ2；
根据喷淋截面的烟气平均流速v1、喷淋截面的喷淋密度ρ1、近壁区域的烟气平均流速v2以及近壁区域的喷淋密度ρ2通过以下公式计算近壁逃逸系数A：
其中，v2大于v1，ρ1大于ρ2。
2.根据权利要求1所述的吸收塔近壁逃逸现象量化方法，其特征在于，所述喷淋截面的烟气平均流速v1以及近壁区域的烟气平均流速v2通过数值模拟的方式计算。
3.根据权利要求1所述的吸收塔近壁逃逸现象量化方法，其特征在于，所述喷淋截面的喷淋密度ρ1和所述近壁区域的喷淋密度ρ2根据喷嘴的配置情况计算。
4.根据权利要求1所述的吸收塔近壁逃逸现象量化方法，其特征在于，所述近壁逃逸系数A越大，近壁逃逸现象越严重；所述近壁逃逸系数A越小，近壁逃逸现象越轻微。
5.根据权利要求1所述的吸收塔近壁逃逸现象量化方法，其特征在于，所述近壁区域宽度为最外层喷嘴喷射圆半径的0.5～1倍。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610945568.XA CN106547968B (zh) | 2016-10-26 | 2016-10-26 | 吸收塔近壁逃逸现象量化方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610945568.XA CN106547968B (zh) | 2016-10-26 | 2016-10-26 | 吸收塔近壁逃逸现象量化方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106547968A CN106547968A (zh) | 2017-03-29 |
| CN106547968B true CN106547968B (zh) | 2020-05-05 |
Family
ID=58393792
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610945568.XA Active CN106547968B (zh) | 2016-10-26 | 2016-10-26 | 吸收塔近壁逃逸现象量化方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106547968B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109012068A (zh) * | 2018-10-17 | 2018-12-18 | 上海外高桥第二发电有限责任公司 | 一种用于治理白色烟羽的排放优化方法 |
| CN116272363B (zh) * | 2023-05-17 | 2023-08-01 | 浙江浙能迈领环境科技有限公司 | 船舶废气混合式脱硫系统及其方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013006125A (ja) * | 2011-06-22 | 2013-01-10 | Babcock Hitachi Kk | 湿式排煙脱硫装置及び湿式排煙脱硫方法 |
| CN103301748A (zh) * | 2013-06-08 | 2013-09-18 | 莱芜市泰钢工程设计研究有限公司 | 湿式喷淋吸收塔添加有机催化剂烟气脱硫脱汞系统及方法 |
| CN103463961A (zh) * | 2013-09-09 | 2013-12-25 | 武汉钢铁(集团)公司 | 气流均布式湿法烟气脱硫吸收塔 |
| CN204816211U (zh) * | 2015-07-13 | 2015-12-02 | 大唐环境产业集团股份有限公司 | 一种气液逃逸阻止装置 |
| CN204952618U (zh) * | 2015-09-11 | 2016-01-13 | 中国神华能源股份有限公司 | 一种脱硫吸收塔 |
-
2016
- 2016-10-26 CN CN201610945568.XA patent/CN106547968B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2013006125A (ja) * | 2011-06-22 | 2013-01-10 | Babcock Hitachi Kk | 湿式排煙脱硫装置及び湿式排煙脱硫方法 |
| CN103301748A (zh) * | 2013-06-08 | 2013-09-18 | 莱芜市泰钢工程设计研究有限公司 | 湿式喷淋吸收塔添加有机催化剂烟气脱硫脱汞系统及方法 |
| CN103463961A (zh) * | 2013-09-09 | 2013-12-25 | 武汉钢铁(集团)公司 | 气流均布式湿法烟气脱硫吸收塔 |
| CN204816211U (zh) * | 2015-07-13 | 2015-12-02 | 大唐环境产业集团股份有限公司 | 一种气液逃逸阻止装置 |
| CN204952618U (zh) * | 2015-09-11 | 2016-01-13 | 中国神华能源股份有限公司 | 一种脱硫吸收塔 |
Non-Patent Citations (2)
| Title |
|---|
| 单塔双循环湿法烟气脱硫传质过程模拟及结构优化研究;秦明臣;《万方学位论文》;20161020;第65-70页 * |
| 脱硫吸收塔流场优化及研究;魏书洲等;《华北电力技术》;20160316;第19页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106547968A (zh) | 2017-03-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| You et al. | Experimental and numerical investigations of shell-side thermo-hydraulic performances for shell-and-tube heat exchanger with trefoil-hole baffles | |
| Yang et al. | Dimensional characteristics of wind effects on the performance of indirect dry cooling system with vertically arranged heat exchanger bundles | |
| Zhang et al. | Convective heat transfer in cross-corrugated triangular ducts under uniform heat flux boundary conditions | |
| Lu et al. | Experimental study of crosswind effects on the performance of small cylindrical natural draft dry cooling towers | |
| Missirlis et al. | Numerical development of a heat transfer and pressure drop porosity model for a heat exchanger for aero engine applications | |
| Mishra et al. | Numerical investigation of air suction through the louvers of a funnel due to high velocity air jet | |
| Kong et al. | Wind leading to improve cooling performance of natural draft air-cooled condenser | |
| Gao et al. | Experimental research on circumferential inflow air and vortex distribution for wet cooling tower under crosswind conditions | |
| CN106547968B (zh) | 吸收塔近壁逃逸现象量化方法 | |
| Liu et al. | Numerical study on the effect of jet spacing on the Swirl flow and heat transfer in the turbine airfoil leading edge region | |
| CN104516998A (zh) | 基于双面壁面边界条件的超音速火焰喷涂的仿真方法 | |
| Chokphoemphun et al. | Experimental study and neural networks prediction on thermal performance assessment of grooved channel air heater | |
| Wang et al. | Flow dead zone analysis and structure optimization for the trefoil-baffle heat exchanger | |
| Sun et al. | A high-precision method for calculating the pressure drop across wire mesh filters | |
| Wang et al. | Physical and numerical study on unsteady shedding behaviors of ventilated partial cavitating flow around an axisymmetric body | |
| Zheng et al. | Thermal performance and heat transfer reliability analysis in helically corrugated helical tube | |
| Kashani et al. | Heat and mass transfer in natural draft cooling towers | |
| Shao et al. | Study of improvements on flow maldistribution of double tube-passes shell-and-tube heat exchanger with rectangular header | |
| Zhou et al. | Numerical simulation of droplets re-entrainment in baffle demister | |
| Ligrani et al. | Spatially resolved heat transfer and friction factors in a rectangular channel with 45-deg angled crossed-rib turbulators | |
| Reuter et al. | Computational fluid dynamics analysis of cooling tower inlets | |
| Hassani et al. | Experimental and numerical analysis of curved turbulators in different arrangements through a rectangular channel | |
| Ke et al. | The optimum angle of attack of delta winglet vortexgenerators on heat transfer performance of finned flat tube bank with considering nonuniform fin temperature | |
| CN108506622A (zh) | 一种基于弧形导流片的低阻力三通构件 | |
| Abushammala et al. | Maximizing mass transfer using highly curved helical pipes: a CFD investigation |
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 |