CN113536211A - 一种用于east托卡马克磁信号快速提取的数值方法 - Google Patents

一种用于east托卡马克磁信号快速提取的数值方法 Download PDF

Info

Publication number
CN113536211A
CN113536211A CN202110840131.0A CN202110840131A CN113536211A CN 113536211 A CN113536211 A CN 113536211A CN 202110840131 A CN202110840131 A CN 202110840131A CN 113536211 A CN113536211 A CN 113536211A
Authority
CN
China
Prior art keywords
magnetic
response signal
magnetic response
east
calculating
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
Application number
CN202110840131.0A
Other languages
English (en)
Other versions
CN113536211B (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.)
Dalian University of Technology
Original Assignee
Dalian University of Technology
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 Dalian University of Technology filed Critical Dalian University of Technology
Priority to CN202110840131.0A priority Critical patent/CN113536211B/zh
Publication of CN113536211A publication Critical patent/CN113536211A/zh
Application granted granted Critical
Publication of CN113536211B publication Critical patent/CN113536211B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/12—Simultaneous equations, e.g. systems of linear equations
    • G—PHYSICS
    • G21—NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21C—NUCLEAR REACTORS
    • G21C17/00—Monitoring; Testing ; Maintaining
    • G21C17/10—Structural combination of fuel element, control rod, reactor core, or moderator structure with sensitive instruments, e.g. for measuring radioactivity, strain
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/10—Nuclear fusion reactors
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00—Energy generation of nuclear origin
    • Y02E30/30—Nuclear fission reactors

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Operations Research (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Plasma & Fusion (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Investigating Or Analyzing Materials By The Use Of Magnetic Means (AREA)
  • Measuring Magnetic Variables (AREA)

Abstract

本发明提供了一种用于EAST托卡马克磁信号快速提取的数值方法。首先根据EAST托卡马克磁探针阵列在空间中的排布情况计算相应的权重系数并储存在控制系统中等待调用。在EAST实验放电过程中，通过磁探针阵列实时测量各处的磁响应信号，将相应位置的磁信号数值与权重系数相乘并通过计算公式提取特定环向分量的复数信号。最后将提取到的信号储存在系统中等待后续分析或传输给其他控制算法进行处理。本发明实现了快速提取EAST放电过程中磁响应信号的特定环向分量，该数值方法计算效率高，达到实时获取磁响应信号特定分量的目的，弥补了传统拟合方法因效率低导致只能通过实验后期处理信号的不足。

Description

一种用于EAST托卡马克磁信号快速提取的数值方法
技术领域
本发明涉及磁约束可控核聚变领域托卡马克装置放电实验数据的快速处理，特别涉及一种用于EAST托卡马克磁响应信号环向分量快速提取的数值方法。
背景技术
随着社会的进步，人类对能源的需求与日俱增，能源短缺问题逐渐凸显。目前人类所使用的能源主要为化石能源，化石能源储量有限且燃烧释放有害气体，因此亟需寻找化石能源的替代品。目前普遍认为核聚变能是解决能源问题的终极方案，其中磁约束法中的托卡马克装置是最有前景的实现方法。我国在此领域处于国际前列，其中中科院等离子体物理研究所的东方超环(EAST)托卡马克是我国重要装置之一。
为了监测以及实时控制EAST托卡马克放电过程中的等离子体状态，需要各种诊断设备来测量等离子体参数，其中磁探针阵列就是用来测量空间中磁响应信号的设备。通过对测量到的磁响应信号进行分析可以判断等离子体状态并采取相应的实时控制手段对等离子体进行控制来达到稳态运行的目的。某些实时控制算法需要磁响应信号的特定环向分量作为输入参数，但是传统的拟合手段效率较低，往往只能在实验后进行拟合分析，无法达到实时提取磁响应信号的特定环向分量的水平，这阻碍了许多实时控制算法的发展。本发明提出的方法具有较高的计算效率，足够实现实时提取磁响应信号特定分量的目的，同时保证了提取信号的准确性，弥补了传统拟合方法只能在实验完成后进行信号处理的不足，为发展各种以磁响应信号特定环向分量为输入的实时控制算法提供了基础，是一种高效且稳定可靠的数值方法。
发明内容
发明目的：实现实时提取EAST放电过程中磁响应信号的特定环向分量的目的，弥补传统拟合方法由于拟合效率较低而只能在实验后期对信号进行处理的不足，为发展多种以磁响应信号特定环向分量为输入参数的实时控制算法提供基础。
本发明的技术方案：
一种用于EAST托卡马克磁信号快速提取的数值方法，该数值方法计算效率高，可以达到实时获取磁响应信号特定分量的目的，同时可以保证提取信号的准确性，是一种高效且稳定可靠的数值方法，具体包括以下步骤：
步骤1：根据EAST托卡马克装置的磁探针阵列空间分布计算相应的权重系数W1、W2…Wn(n为环向磁探针数，通常为偶数)；
步骤1.1：将n个磁探针分为两两一组
Figure BDA0003178567600000021
Figure BDA0003178567600000022
共
Figure BDA0003178567600000023
组，并分别计算每组的权重系数；
步骤1.2：首先计算第一组磁探针所对应的权重系数W1和
Figure BDA0003178567600000024
Figure BDA0003178567600000025
上测量到的磁响应信号S1和
Figure BDA0003178567600000026
满足如下方程组：
acosφ1-bsinφ1＝S1
acosφn/2+1-bsinφn/2+1＝Sn/2+1
其中a、b为待定系数(用于计算所需提取的磁响应信号)，S为不同磁探针测量到的磁响应信号，φ为磁探针所在空间位置所对应的环向角大小。最终提取的磁响应信号环向分量为
Figure BDA0003178567600000027
其中
Figure BDA0003178567600000028
i为虚数符号；
求解如上方程组可得待定系数如下：
Figure BDA0003178567600000029
Figure BDA00031785676000000210
因此相应的权重系数W定义如下：
Figure BDA00031785676000000211
Figure BDA0003178567600000031
步骤1.3：根据步骤1.2中采用的方法计算剩余组磁探针所对应的权重系数；
步骤2：将步骤1中得到的权重系数W1、W2…Wn存储到控制系统中等待调用；
步骤3：在EAST托卡马克装置实验放电过程中，通过磁探针阵列实时测量t时刻空间中的磁响应信号S1、S2…Sn；
步骤4：调用步骤2中存储的权重系数并根据如下公式计算该时刻的磁响应信号环向分量Ct：
Figure BDA0003178567600000032
Figure BDA0003178567600000033
Figure BDA0003178567600000034
Figure BDA0003178567600000035
Figure BDA0003178567600000036
步骤5：将该磁响应信号环向分量Ct存储在系统中留待实验后分析或实时发送到需要该信号的实时控制算法作后续处理。
步骤6：重复进行步骤3-5即可实现实时提取磁响应信号特定环向分量的目的。
本发明的有益效果：本发明实现了实时提取EAST托卡马克装置放电过程中磁响应信号的特定环向分量的目的，弥补了传统拟合方法效率低的不足，计算效率高且同时保证信号准确性，是一种稳定高效的数值方法。
附图说明
图1本发明所适用的EAST托卡马克装置的横截面示意图。
图2本发明所适用的EAST托卡马克装置磁探针阵列环向空间分布示意图(俯视图)。
图3本发明快速提取的磁响应信号与传统拟合得到的磁响应信号的对比。
图4本发明实现实时提取EAST实验过程中磁响应信号特定环向分量的主要流程图。
具体实施方式
下面结合附图和实施例对本发明进行详细说明。
EAST托卡马克装置具有如附图1所示截面位形，其中芯部位置为实验放电过程中包含等离子体的区域。等离子体外围装配多组磁探针阵列用于测量实验过程中产生的磁响应信号。每组磁探针阵列包含多个磁探针，磁探针沿着环向均匀分布，如附图2所示(以其中一组磁探针阵列为例)。环向均匀分布8个磁探针(Sen1、Sen2…Sen8)，相邻磁探针环向相距
Figure BDA0003178567600000041
磁探针测量的信号通常无法直接使用，需要对特定环向分量进行提取。本实施例通过计算各磁探针的权重系数来提高提取效率，可以实现快速提取特定环向分量的目的。首先将8个磁探针分成4组(Sen1-Sen5、Sen2-Sen6、Sen3-Sen7、Sen4-Sen8)，根据下面具体实施步骤1中的方法计算相应的权重系数可得：
Figure BDA0003178567600000042
Figure BDA0003178567600000043
Figure BDA0003178567600000044
Figure BDA0003178567600000045
Figure BDA0003178567600000046
Figure BDA0003178567600000047
Figure BDA0003178567600000051
Figure BDA0003178567600000052
其中φ1-φ8分别为0、
Figure BDA0003178567600000053
π、
Figure BDA0003178567600000054
在实验之前将权重系数存储到控制系统中等待调用。实验放电开始后，通过磁探针阵列实时测量t时刻空间中的磁响应信号S1、S2…S8，并采用上面计算得到的权重系数通过具体实施步骤4中的方法计算该时刻的磁响应信号环向分量Ct，然后将提取信号存储或发送到其他实时控制算法。重复测量、提取和存储发送等步骤多次直到实验放电结束即可达到实时获取磁响应信号特定环向分量的目的。与传统拟合方法(单次计算时间为秒量级)相比，该方法(单次计算时间为毫秒量级)极大提高了提取效率，并且保证了提取结果的准确性，提取结果与传统拟合结果十分接近，如附图3所示。
具体实施步骤如下：
步骤1：根据EAST托卡马克装置的磁探针阵列空间分布计算相应的权重系数W1、W2…Wn(n为环向磁探针数，通常为偶数)；
步骤1.1：将n个磁探针分为两两一组
Figure BDA0003178567600000055
Figure BDA0003178567600000056
共
Figure BDA0003178567600000057
组，并分别计算每组的权重系数；
步骤1.2：首先计算第一组磁探针所对应的权重系数W1和
Figure BDA0003178567600000058
Figure BDA0003178567600000059
上测量到的磁响应信号S1和
Figure BDA00031785676000000510
满足如下方程组：
acosφ1-bsinφ1＝S1
acosφn/2+1-bsinφn/2+1＝Sn/2+1
其中a、b为待定系数(用于计算所需提取的磁响应信号)，S为不同磁探针测量到的磁响应信号，φ为磁探针所在空间位置所对应的环向角大小。最终提取的磁响应信号环向分量为
Figure BDA00031785676000000511
其中
Figure BDA00031785676000000512
i为虚数符号；
求解如上方程组可得待定系数如下：
Figure BDA0003178567600000061
Figure BDA0003178567600000062
因此相应的权重系数W定义如下：
Figure BDA0003178567600000063
Figure BDA0003178567600000064
步骤1.3：根据步骤1.2中采用的方法计算剩余组磁探针所对应的权重系数；
步骤2：将步骤1中得到的权重系数W1、W2…Wn存储到控制系统中等待调用；
步骤3：在EAST托卡马克装置实验放电过程中，通过磁探针阵列实时测量t时刻空间中的磁响应信号S1、S2…Sn；
步骤4：调用步骤2中存储的权重系数并根据如下公式计算该时刻的磁响应信号环向分量Ct：
Figure BDA0003178567600000065
Figure BDA0003178567600000066
Figure BDA0003178567600000067
Figure BDA0003178567600000068
Figure BDA0003178567600000069
步骤5：将该磁响应信号环向分量Ct存储在系统中留待实验后分析或实时发送到需要该信号的实时控制算法作后续处理。
步骤6：重复进行步骤3-5即可实现实时提取磁响应信号特定环向分量的目的。
以上内容是结合优选技术方案对本发明所做的进一步详细说明，不能认定发明的具体实施仅限于这些说明。对本发明所属技术领域的普通技术人员来说，在不脱离本发明的构思的前提下，还可以做出简单的推演及替换，都应当视为本发明的保护范围。

Claims (1)

1.一种用于EAST托卡马克磁信号快速提取的数值方法，其特征在于,该数值方法实现了快速提取EAST托卡马克实验放电过程中磁探针阵列测量到的磁信号的特定分量，具体包括以下步骤：
步骤1：根据EAST托卡马克装置的磁探针阵列空间分布计算相应的权重系数W1、W2…Wn，n为环向磁探针数，为偶数；
步骤1.1：将n个磁探针分为两两一组
Figure FDA0003178567590000011
共
Figure FDA0003178567590000012
组，并分别计算每组的权重系数；
步骤1.2：首先计算第一组磁探针所对应的权重系数W1和
Figure FDA0003178567590000013
Figure FDA0003178567590000014
上测量到的磁响应信号S1和
Figure FDA0003178567590000015
满足如下方程组：
acosφ1-bsinφ1＝S1
acosφn/2+1-bsinφn/2+1＝Sn/2+1
其中，a、b为待定系数，用于计算所需提取的磁响应信号；S为不同磁探针测量到的磁响应信号，φ为磁探针所在空间位置所对应的环向角大小；最终提取的磁响应信号环向分量为
Figure FDA0003178567590000016
其中
Figure FDA0003178567590000017
i为虚数符号；
求解如上方程组得待定系数如下：
Figure FDA0003178567590000018
Figure FDA0003178567590000019
因此相应的权重系数W定义如下：
Figure FDA00031785675900000110
Figure FDA00031785675900000111
步骤1.3：根据步骤1.2中采用的方法计算剩余组磁探针所对应的权重系数；
步骤2：将步骤1中得到的权重系数W1、W2…Wn存储到控制系统中，等待调用；
步骤3：在EAST托卡马克装置实验放电过程中，通过磁探针阵列实时测量t时刻空间中的磁响应信号S1、S2…Sn；
步骤4：调用步骤2中存储的权重系数并根据如下公式计算该时刻的磁响应信号环向分量Ct：
Figure FDA0003178567590000021
Figure FDA0003178567590000022
Figure FDA0003178567590000023
Figure FDA0003178567590000024
Figure FDA0003178567590000025
步骤5：将该磁响应信号环向分量Ct存储在系统中留待实验后分析或实时发送到需要该信号的实时控制算法作后续处理；
步骤6：重复进行步骤3-5，即实现实时提取磁响应信号特定环向分量。
CN202110840131.0A 2021-07-24 2021-07-24 一种用于east托卡马克磁信号快速提取的数值方法 Active CN113536211B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110840131.0A CN113536211B (zh) 2021-07-24 2021-07-24 一种用于east托卡马克磁信号快速提取的数值方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110840131.0A CN113536211B (zh) 2021-07-24 2021-07-24 一种用于east托卡马克磁信号快速提取的数值方法

Publications (2)

Publication Number Publication Date
CN113536211A true CN113536211A (zh) 2021-10-22
CN113536211B CN113536211B (zh) 2024-02-06

Family

ID=78088878

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110840131.0A Active CN113536211B (zh) 2021-07-24 2021-07-24 一种用于east托卡马克磁信号快速提取的数值方法

Country Status (1)

Country Link
CN (1) CN113536211B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116679863A (zh) * 2023-07-28 2023-09-01 中国科学院合肥物质科学研究院 一种适用于east辐射偏滤器充气系统的数据处理系统

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5804965A (en) * 1995-09-25 1998-09-08 The United States Of America As Represented By The United States Department Of Energy Method and apparatus for steady-state magnetic measurement of poloidal magnetic field near a tokamak plasma
RU2191410C2 (ru) * 2000-10-16 2002-10-20 Российский научный центр "Курчатовский институт" Система регистрации рентгеновского излучения в установках типа токамак
CN106650951A (zh) * 2016-12-29 2017-05-10 中国科学技术大学 基于时域全局相似度的多通道测量数据自动筛选方法
CN106991271A (zh) * 2017-03-07 2017-07-28 中国科学院合肥物质科学研究院 一种适用于east偏滤器探针诊断数据处理的软件系统
CN107896412A (zh) * 2017-09-27 2018-04-10 中国科学院合肥物质科学研究院 一种用于测量射频波特性的高频磁探针诊断系统
CN110232205A (zh) * 2019-04-28 2019-09-13 大连理工大学 用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法
CN112930015A (zh) * 2021-03-17 2021-06-08 核工业西南物理研究院 一种等离子体电磁特性诊断用电磁复合探针
CN112992385A (zh) * 2021-05-06 2021-06-18 西南交通大学 一种准环对称仿星器磁场位形的设计方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5804965A (en) * 1995-09-25 1998-09-08 The United States Of America As Represented By The United States Department Of Energy Method and apparatus for steady-state magnetic measurement of poloidal magnetic field near a tokamak plasma
RU2191410C2 (ru) * 2000-10-16 2002-10-20 Российский научный центр "Курчатовский институт" Система регистрации рентгеновского излучения в установках типа токамак
CN106650951A (zh) * 2016-12-29 2017-05-10 中国科学技术大学 基于时域全局相似度的多通道测量数据自动筛选方法
CN106991271A (zh) * 2017-03-07 2017-07-28 中国科学院合肥物质科学研究院 一种适用于east偏滤器探针诊断数据处理的软件系统
CN107896412A (zh) * 2017-09-27 2018-04-10 中国科学院合肥物质科学研究院 一种用于测量射频波特性的高频磁探针诊断系统
CN110232205A (zh) * 2019-04-28 2019-09-13 大连理工大学 用于托卡马克中共振磁扰动控制新经典撕裂模的模拟方法
CN112930015A (zh) * 2021-03-17 2021-06-08 核工业西南物理研究院 一种等离子体电磁特性诊断用电磁复合探针
CN112992385A (zh) * 2021-05-06 2021-06-18 西南交通大学 一种准环对称仿星器磁场位形的设计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
李, 罗家融, 王华忠, 黄勤超: "EAST上极向场线圈位置测量方法的研究与设计", 计算机测量与控制, no. 08 *
李凤娟;朱翔;杨建华;高翔;: "应用高速CCD诊断研究HT-7装置边界等离子体旋转", 核聚变与等离子体物理, no. 03 *
王正汹等: "多螺旋新经典撕裂模大规模模拟研究", 《第十八届全国等离子体科学技术会议摘要集》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116679863A (zh) * 2023-07-28 2023-09-01 中国科学院合肥物质科学研究院 一种适用于east辐射偏滤器充气系统的数据处理系统
CN116679863B (zh) * 2023-07-28 2023-10-24 中国科学院合肥物质科学研究院 一种适用于east辐射偏滤器充气系统的数据处理系统

Also Published As

Publication number Publication date
CN113536211B (zh) 2024-02-06

Similar Documents

Publication Publication Date Title
Fujimoto et al. Mapping neutron star data to the equation of state using the deep neural network
CN107273582B (zh) 一种用于快中子反应堆中子输运燃耗耦合分析的计算方法
CN113536211A (zh) 一种用于east托卡马克磁信号快速提取的数值方法
FR2868866A1 (fr) Procede et dispositif de surveillance du coeur d'un reacteur nucleaire
CN110060788B (zh) 一种通用的热离子空间核反应堆电源瞬态热电特性分析方法
EP3074982B1 (en) Method and system for generating a nuclear reactor core loading distribution
CN109800446B (zh) 一种锂离子电池放电过程电压不一致性估计方法及装置
CN116070505B (zh) 一种基于模型降阶与堆外计数的堆芯功率在线监测方法
CN110739091A (zh) 一种核电厂次临界度直接测量系统及硬件平台和测量方法
CN106500991B (zh) 基于自适应多尺度AVG-Hat变换的轴承故障信号特征提取方法
CN106405622B (zh) 一种快速核截面多普勒展宽方法
JP2016156740A (ja) 共鳴計算プログラム及び解析装置
KR102105037B1 (ko) 원자로 크러드 분석 장치 및 방법
CN116754213A (zh) 一种基于强噪声背景下的电动闸阀故障诊断方法、装置及设备
CN115017741B (zh) 一种燃料电池气体扩散层的重构方法、装置和电子设备
CN108008262B (zh) 一种局部放电快速检测方法
CN106897481B (zh) 利用空间多点功率预测特征的风功率集成预测方法及装置
CN110866346B (zh) 一种用于获取干式空心电抗器固有振动特性的方法及系统
CN114237043A (zh) 基于深度学习的燃气轮机设备传递函数闭环辨识方法
FR3073625B1 (fr) Procede et dispositif electronique de simulation d'une reponse ultrasonore d'une piece metallique, programme d'ordinateur et procede et systeme electronique de controle associes
CN106897790B (zh) 一种选择最优预报格点的风功率预测方法及装置
Ferreira et al. Deep learning for plasma tomography in nuclear fusion
CN117093834B (zh) 一种三相电压系统异常监测方法
CN111400869B (zh) 一种堆芯中子通量时空演变预测方法、装置、介质及设备
Ba et al. DRSN Fan Rolling Bearing Fault Diagnosis Method Based on DWT with Improvement

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