CN110134906A - 一种非对称功率分布下堆芯dnbr计算方法 - Google Patents

一种非对称功率分布下堆芯dnbr计算方法 Download PDF

Info

Publication number
CN110134906A
CN110134906A CN201910168390.6A CN201910168390A CN110134906A CN 110134906 A CN110134906 A CN 110134906A CN 201910168390 A CN201910168390 A CN 201910168390A CN 110134906 A CN110134906 A CN 110134906A
Authority
CN
China
Prior art keywords
analyzed
power
component
reactor core
channel
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
CN201910168390.6A
Other languages
English (en)
Other versions
CN110134906B (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.)
China Nuclear Power Engineering Co Ltd
Original Assignee
China Nuclear Power 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 China Nuclear Power Engineering Co Ltd filed Critical China Nuclear Power Engineering Co Ltd
Priority to CN201910168390.6A priority Critical patent/CN110134906B/zh
Publication of CN110134906A publication Critical patent/CN110134906A/zh
Application granted granted Critical
Publication of CN110134906B publication Critical patent/CN110134906B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • GPHYSICS
    • G21NUCLEAR PHYSICS; NUCLEAR ENGINEERING
    • G21CNUCLEAR REACTORS
    • G21C17/00Monitoring; Testing ; Maintaining
    • G21C17/02Devices or arrangements for monitoring coolant or moderator
    • G21C17/038Boiling detection in moderator or coolant
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E30/00Energy generation of nuclear origin
    • Y02E30/30Nuclear fission reactors

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Plasma & Fusion (AREA)
  • Monitoring And Testing Of Nuclear Reactors (AREA)

Abstract

本发明属于核电站安全技术领域,具体涉及一种非对称功率分布下堆芯DNBR计算方法。包括步骤S101,按照平面直角坐标系对反应堆的堆芯进行划分,筛选待分析组件并调整到堆芯的第一象限内;按照平面直角坐标系对待分析组件中的热组件进行划分,筛选出热组件中的待分析通道并调整到热组件的第四象限内;步骤S102,对待分析组件功率采用中心化的原则进行重新调整与布置,保证全堆芯组件功率的归一化;步骤S103,对热组件中的待分析通道的功率分布进行重新调整与布置,保证热组件的第四象限内的待分析通道的功率的归一化;步骤S104,得到用于DNBR计算的堆芯的第一象限的组件功率分布情况、热组件的第四象限的功率分布情况,调用DNBR计算程序完成DNBR计算。

Description

一种非对称功率分布下堆芯DNBR计算方法
技术领域
本发明属于核电站安全技术领域,具体涉及一种非对称功率分布下堆芯DNBR计算方法。
背景技术
核电厂在进行反应堆堆芯热工水力设计以及事故分析中,堆芯的偏离泡核沸腾比(DNBR)的计算至关重要。DNBR是评判I类工况、II工况以及个别III类、IV类工况是否满足验收准则的重要参数。
在通常的稳态堆芯子通道分析中,堆芯功率分布基本上都是对称的。但在核电厂燃料循环换料操作以及诸如主蒸汽管道断裂(MSLB)等事故中,堆芯的功率分布是非对称的。因此需要研究一种能够用于此非对称功率分布的堆芯DNBR计算方法。
发明内容
针对目前无法在非对称功率分布的情况下进行堆芯DNBR计算方法的弊端,本发明的目的是提供一种堆芯DNBR计算方法,能够在堆芯功率分布不对称的情况下,利用现有的DNBR计算程序进行堆芯DNBR的计算。
为达到以上目的,本发明采用的技术方案是一种非对称功率分布下堆芯DNBR计算方法,用于对反应堆的堆芯进行调整与部属,配合DNBR计算程序完成DNBR计算,所述堆芯包括若干组件,所述组件内包括若干放置有一根燃料棒的通道,包括以下步骤:
步骤S101,按照平面直角坐标系对反应堆的堆芯进行划分,筛选出全堆芯功率最高的1/4所述组件作为待分析组件,调整到所述堆芯的第一象限内;
按照平面直角坐标系对所述待分析组件中的功率最高的热组件进行划分,筛选出所述热组件中包含有功率最高的1/4所述燃料棒的通道作为待分析通道,并调整到所述热组件的第四象限内;
步骤S102,将所述热组件置于所述堆芯的平面直角坐标系的原点的位置,其它所述待分析组件根据功率值从大到小在所述堆芯的第一象限内围绕所述热组件进行环向分布;
步骤S103,对所述热组件中的所述待分析通道的功率分布进行重新调整与布置,保证所述热组件的第四象限内的所述待分析通道的功率的归一化;
步骤S104,在所述步骤S102与所述步骤S103的基础上,得到最终可用于DNBR计算的堆芯的所述堆芯的第一象限内的待分析堆芯的功率分布情况,以及得到最终可用于DNBR计算的所述热组件的第四象限的功率分布情况,调用所述DNBR计算程序,最终完成所述DNBR计算;所述待分析堆芯为所述堆芯的第一象限内的所述待分析组件的一部分,所述待分析堆芯包括所述热组件。
进一步,
所述步骤S102共需进行三次调整,第一次调整是将堆芯的所述堆芯的第一象限内的、包括所述热组件在内的所有所述待分析组件进行归一化,保证所述待分析组件的功率的平均值为1.0;
第二次调整是将所述热组件的功率调整至原全堆芯组件中的最大值,并赋予所述堆芯的第一象限内的其它所述待分析组件相应的调整因子,保证所述堆芯的第一象限内的堆芯的归一化;
第三次调整是将所述热组件与其它所述待分析组件的位置根据所述中心化的原则进行重新布置。
进一步,
所述步骤S103共需进行二次调整,第一次调整是将所述热组件的第四象限内的所述待分析通道的功率进行归一化,保证所述待分析通道的功率的平均值在所述归一化的前后保持一致;第二次调整是将所述最热通道中的燃料棒的功率调整至原最高值,将其它所述待分析通道赋予相应的调整因子,保证所述热组件的第四象限内的所述待分析通道的功率的归一化;所述最热通道是指所述待分析通道中功率最高的一个。
进一步,在所述步骤S104中,根据所述步骤S102与所述步骤S103,可得到用于DNBR计算的所述待分析堆芯在所述堆芯的第一象限的组件功率分布情况和所述热组件的第四象限的功率分布情况,此时直接调用所述DNBR计算程序,可完成所述DNBR计算。
进一步,在所述步骤S102中,对所述待分析组件功率进行重新调整与布置,采用中心化的原则;在分布过程中保证全堆芯的所述组件的功率的归一化;所述归一化是指对全堆芯的各个所述组件的功率进行调整时始终保持全部所述组件的总功率不变。
进一步,在所述步骤S103中,所述归一化是指对所述热组件的第四象限内的所述待分析通道的功率进行调整时始终保持所述热组件的第四象限内的全部所述待分析通道的总功率不变。
本发明的有益效果在于:
1.能够在堆芯的功率分布处于非对称的情况下进行DNBR计算。
2.能够在对堆芯进行燃料循环时计算得到相应的DNBR结果,进一步保证了换料操作的安全进行。
3.充分利用了现有的DNBR计算程序输入文件,无需另外进行输入文件的模型修改。
附图说明
图1是本发明具体实施方式中所提供的一种非对称功率分布下堆芯DNBR计算方法的流程图;
图2是本发明具体实施方式中的步骤S101中,按照平面直角坐标系对反应堆的堆芯进行划分,并将待分析组件2调整到堆芯的第一象限内以后的堆芯布置示意图(也就是全堆芯各组件功率分布示意图);其中,组件1、待分析组件2上的数值代表其功率;
图3是本发明具体实施方式中的步骤S101中,按照平面直角坐标系对热组件3进行划分,并将热组件3中的待分析通道5调整到热组件3的第四象限内以后的通道功率分布示意图;其中,通道4、待分析通道5、最热通道6上的数值代表其功率;
图4是本发明具体实施方式中的步骤S102中的第一次调整后,堆芯的第一象限内的待分析组件2的功率分布示意图;
图5是本发明具体实施方式中的步骤S102中的第二次调整后,堆芯的第一象限内的待分析组件2的功率分布示意图;
图6是本发明具体实施方式中的步骤S102中的第三次调整后,堆芯的第一象限内的待分析组件2的功率分布示意图;
图7是本发明具体实施方式中的步骤S103中的第一次调整后,热组件3的第四象限内的待分析通道5的功率分布示意图;
图8是本发明具体实施方式中的步骤S103中的第二次调整后,热组件3的第四象限内的待分析通道5的功率分布示意图;
图9是本发明具体实施方式中经过步骤102、步骤103后所得到的热组件3的第四象限内的待分析通道5的功率分布示意图;也就是最终用于DNBR计算的1/4热组件的功率分布的示意图;
图10是本发明具体实施方式中经过步骤102、步骤103后所得到的最终用于DNBR计算的待分析堆芯(1/4堆芯)的功率分布的示意图;
图中:1-组件,2-待分析组件,3-热组件,4-通道,5-待分析通道,6-最热通道,7-控制棒(用于堆芯功率的控制)。
具体实施方式
下面结合附图和实施例对本发明作进一步描述。
本发明提供的一种非对称功率分布下堆芯DNBR计算方法(见图1),包括以下步骤:
步骤S101,按照平面直角坐标系对反应堆的堆芯进行划分,反应堆的堆芯包括若干横截面为矩形的组件1,组件1内包括若干垂直于横截面、内部放置有一根燃料棒的通道4;梳理上游专业提供的全堆芯以及热组件的功率分布(此处上游专业指的是物理专业),完成标准分布区域转换,筛选出全堆芯功率最高的1/4组件作为待分析组件2,将待分析组件2调整到堆芯的第一象限内,见图2;
按照平面直角坐标系对待分析组件2中的功率最高的热组件3进行划分(热组件3就是全部的待分析组件2中的功率最高的一个),筛选出热组件3中包含有功率最高的燃料棒所在的通道4作为待分析通道5,将待分析通道5调整到热组件3的第四象限内,见图3;
步骤S102,对待分析组件2功率进行重新调整与布置,采用中心化的原则,将热组件3置于堆芯的平面直角坐标系的原点的位置,其它待分析组件2根据功率值从大到小在堆芯的第一象限内围绕热组件3进行环向分布,在分布过程中保证全堆芯的热组件的功率的归一化;归一化是指对全堆芯的各个热组件的功率进行调整时始终保持全部热组件的总功率不变(例如某些组件的功率调高一些,相应的另一些组件的功率就调低一些,但所有组件的总功率整体上不发生变化);
其中,步骤S102共需进行三次调整,第一次调整是将堆芯的第一象限内的、包括热组件3在内的所有待分析组件2进行归一化,保证待分析组件2的功率的平均值为1.0,见图4;
第二次调整是将热组件3的功率调整至原全堆芯组件中的最大值,并赋予堆芯的第一象限内的其它待分析组件2相应的调整因子,保证堆芯的第一象限内的堆芯的归一化,见图5;
第三次调整是将热组件3与其它待分析组件2的位置根据中心化的原则进行重新布置,见图6。
步骤S103,对热组件3中的待分析通道5的功率分布进行重新调整与布置,保证热组件3的第四象限内的待分析通道5的功率的归一化;归一化是指对热组件3的第四象限内的待分析通道5的功率进行调整时始终保持热组件3的第四象限内的全部待分析通道5的总功率不变;
步骤S103共需进行二次调整,第一次调整是将热组件3的第四象限内的待分析通道5的功率进行归一化,保证待分析通道5的功率的平均值在归一化的前后保持一致,见图7;第二次调整是将最热通道6中的燃料棒的功率调整至原最高值(最热通道6是指待分析通道5中功率最高的一个),将其它待分析通道5赋予相应的调整因子,保证热组件3的第四象限内的待分析通道5的功率的归一化,见图8。
步骤S104,在步骤S102与步骤S103的基础上,得到最终可用于DNBR计算的堆芯的第一象限内的待分析堆芯的功率分布情况(见图10),以及得到最终可用于DNBR计算的热组件3的第四象限的功率分布情况(见图9),最后调用原DNBR计算程序,最终完成DNBR计算。
本发明所述的装置并不限于具体实施方式中所述的实施例,本领域技术人员根据本发明的技术方案得出其他的实施方式,同样属于本发明的技术创新范围。

Claims (6)

1.一种非对称功率分布下堆芯DNBR计算方法,用于对反应堆的堆芯进行调整与部属,配合DNBR计算程序完成DNBR计算,所述堆芯包括若干组件(1),所述组件(1)内包括若干放置有一根燃料棒的通道(4),包括以下步骤:
步骤S101,按照平面直角坐标系对反应堆的堆芯进行划分,筛选出全堆芯功率最高的1/4所述组件(1)作为待分析组件(2),调整到所述堆芯的第一象限内;
按照平面直角坐标系对所述待分析组件(2)中的功率最高的热组件(3)进行划分,筛选出所述热组件(3)中包含有功率最高的1/4所述燃料棒的通道(4)作为待分析通道(5),并调整到所述热组件(3)的第四象限内;
步骤S102,将所述热组件(3)置于所述堆芯的平面直角坐标系的原点的位置,其它所述待分析组件(2)根据功率值从大到小在所述堆芯的第一象限内围绕所述热组件(3)进行环向分布;
步骤S103,对所述热组件(3)中的所述待分析通道(5)的功率分布进行重新调整与布置,保证所述热组件(3)的第四象限内的所述待分析通道(5)的功率的归一化;
步骤S104,在所述步骤S102与所述步骤S103的基础上,得到最终可用于DNBR计算的所述堆芯的第一象限内的待分析堆芯的功率分布情况,以及得到最终可用于DNBR计算的所述热组件(3)的第四象限的功率分布情况,调用所述DNBR计算程序,最终完成所述DNBR计算;所述待分析堆芯为所述堆芯的第一象限内的所述待分析组件(2)的一部分,所述待分析堆芯包括所述热组件(3)。
2.如权利要求1所述的方法,其特征是:
所述步骤S102共需进行三次调整,第一次调整是将所述堆芯的第一象限内的、包括所述热组件(3)在内的所有所述待分析组件(2)进行归一化,保证所述待分析组件(2)的功率的平均值为1.0;
第二次调整是将所述热组件(3)的功率调整至原全堆芯组件中的最大值,并赋予所述堆芯的第一象限内的其它所述待分析组件(2)相应的调整因子,保证所述堆芯的第一象限内的堆芯的归一化;
第三次调整是将所述热组件(3)与其它所述待分析组件(2)的位置根据所述中心化的原则进行重新布置。
3.如权利要求1所述的方法,其特征是:
所述步骤S103共需进行二次调整,第一次调整是将所述热组件(3)的第四象限内的所述待分析通道(5)的功率进行归一化,保证所述待分析通道(5)的功率的平均值在所述归一化的前后保持一致;第二次调整是将所述最热通道(6)中的燃料棒的功率调整至原最高值,将其它所述待分析通道(5)赋予相应的调整因子,保证所述热组件(3)的第四象限内的所述待分析通道(5)的功率的归一化;所述最热通道(6)是指所述待分析通道(5)中功率最高的一个。
4.如权利要求1所述的方法,其特征是:在所述步骤S104中,根据所述步骤S102与所述步骤S103,可得到用于DNBR计算的所述待分析堆芯在所述堆芯的第一象限的功率分布情况和所述热组件(3)的第四象限的功率分布情况,此时直接调用所述DNBR计算程序,完成所述DNBR计算。
5.如权利要求1所述的方法,其特征是:在所述步骤S102中,对所述待分析组件(2)功率进行重新调整与布置,采用中心化的原则;在分布过程中保证全堆芯的所述组件的功率的归一化;所述归一化是指对全堆芯的各个所述组件的功率进行调整时始终保持全部所述组件的总功率不变。
6.如权利要求1所述的方法,其特征是:在所述步骤S103中,所述归一化是指对所述热组件(3)的第四象限内的所述待分析通道(5)的功率进行调整时始终保持所述热组件(3)的第四象限内的全部所述待分析通道(5)的总功率不变。
CN201910168390.6A 2019-03-06 2019-03-06 一种非对称功率分布下堆芯dnbr计算方法 Active CN110134906B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910168390.6A CN110134906B (zh) 2019-03-06 2019-03-06 一种非对称功率分布下堆芯dnbr计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910168390.6A CN110134906B (zh) 2019-03-06 2019-03-06 一种非对称功率分布下堆芯dnbr计算方法

Publications (2)

Publication Number Publication Date
CN110134906A true CN110134906A (zh) 2019-08-16
CN110134906B CN110134906B (zh) 2024-07-12

Family

ID=67568517

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910168390.6A Active CN110134906B (zh) 2019-03-06 2019-03-06 一种非对称功率分布下堆芯dnbr计算方法

Country Status (1)

Country Link
CN (1) CN110134906B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110633454A (zh) * 2019-09-19 2019-12-31 中国核动力研究设计院 一种基于修正法的chf关系式dnbr限值统计学确定方法
CN110727920A (zh) * 2019-09-19 2020-01-24 中国核动力研究设计院 一种基于分组法的chf关系式dnbr限值统计学确定方法
CN111553022A (zh) * 2020-04-30 2020-08-18 中国核动力研究设计院 瞬变外力对动态自反馈条件下堆芯dnbr限值修正系数的整体式求解方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02251791A (ja) * 1989-03-27 1990-10-09 Nuclear Fuel Ind Ltd 原子炉炉心
WO1995011509A1 (fr) * 1993-10-22 1995-04-27 Japan Atomic Energy Research Institute Cycle du combustible nucleaire
JP2006208391A (ja) * 2006-03-08 2006-08-10 Toshiba Corp 燃料集合体および原子炉の炉心
CN103400618A (zh) * 2013-07-01 2013-11-20 中国核电工程有限公司 一种vver核电站反应堆堆芯非对称装载方法
CN107977504A (zh) * 2017-11-28 2018-05-01 广东核电合营有限公司 一种非对称堆芯燃料管理计算方法、装置及终端设备

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH02251791A (ja) * 1989-03-27 1990-10-09 Nuclear Fuel Ind Ltd 原子炉炉心
WO1995011509A1 (fr) * 1993-10-22 1995-04-27 Japan Atomic Energy Research Institute Cycle du combustible nucleaire
JP2006208391A (ja) * 2006-03-08 2006-08-10 Toshiba Corp 燃料集合体および原子炉の炉心
CN103400618A (zh) * 2013-07-01 2013-11-20 中国核电工程有限公司 一种vver核电站反应堆堆芯非对称装载方法
CN107977504A (zh) * 2017-11-28 2018-05-01 广东核电合营有限公司 一种非对称堆芯燃料管理计算方法、装置及终端设备

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110633454A (zh) * 2019-09-19 2019-12-31 中国核动力研究设计院 一种基于修正法的chf关系式dnbr限值统计学确定方法
CN110727920A (zh) * 2019-09-19 2020-01-24 中国核动力研究设计院 一种基于分组法的chf关系式dnbr限值统计学确定方法
CN110727920B (zh) * 2019-09-19 2022-08-19 中国核动力研究设计院 一种基于分组法的chf关系式dnbr限值统计学确定方法
CN110633454B (zh) * 2019-09-19 2022-10-21 中国核动力研究设计院 一种基于修正法的chf关系式dnbr限值统计学确定方法
CN111553022A (zh) * 2020-04-30 2020-08-18 中国核动力研究设计院 瞬变外力对动态自反馈条件下堆芯dnbr限值修正系数的整体式求解方法

Also Published As

Publication number Publication date
CN110134906B (zh) 2024-07-12

Similar Documents

Publication Publication Date Title
CN110134906A (zh) 一种非对称功率分布下堆芯dnbr计算方法
CN103378601A (zh) 一种基于bang-bang控制的一次调频方法及装置
CN107066612A (zh) 一种基于SparkJoin操作的自适应数据倾斜调整方法
CN104358625B (zh) 一种变推力发动机推力闭环控制的方法
CN107243299A (zh) 一种反应装置和一种反应系统以及一种液相加氢反应方法
CN103943161B (zh) 用于重水反应堆的多相堆芯设计和钍基燃料
CN110322975A (zh) 一种核反应堆初始堆芯无源启动方法
CN103106931B (zh) 一种适用于超临界水冷堆双流程堆芯的堆内构件
CN110826205B (zh) 一种适用于压水堆Pin-by-pin计算的能群结构优化方法
Sekimoto et al. Power level control of CANDLE reactor without control rods
CN106294759A (zh) 一种自动获取群体用户作息时间的方法及系统
CN106886458A (zh) 一种集群资源分配方法和装置
CN106779389A (zh) 一种温室气体核算标准管理模块
Song et al. Calculation of Thermal-Fluid Parameters for the Pre-designed VHTR Core Based on CORONA
DE1221373B (de) Einrichtung zum Regeln von Kernreaktoren
CN105677977B (zh) 一种通用锂离子动力电池单体有限元仿真前处理方法
Hykes et al. Coupled CASMO5-VIPRE Analysis of Non-Uniform Void Distributions in BWR Lattices
EP2472523A2 (en) Optimized nuclear fuel bundle supports and methods of creating reactor cores with the same
Kobzar et al. Experimental investigations on the efficiency of heat exchange intensifying spacer grids
CN105943067B (zh) 一种x射线质量对比评估装置及方法
Sun et al. Preliminary analysis of effect of Vanadium Self-Powered Neutron Detectors on AP1000 reactor core physical parameters
ATZ et al. Ducted Assembly Steady State Heat Transfer Software
Wang et al. Coupled Three Dimensional Neutronics/Thermal-Hydraulics Code STTA for SCWR Core Transient Analysis
Bae et al. Thermal-hydraulic characteristics for fuel channels using burnable poison in the CANDU reactor
DE102010053263A1 (de) Vorrichtung zum Mischen von heißen Gasströmen

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