CN110032822B - 乏池失去部分冷却后温度及温升速率计算的分析方法 - Google Patents
乏池失去部分冷却后温度及温升速率计算的分析方法 Download PDFInfo
- Publication number
- CN110032822B CN110032822B CN201910322489.7A CN201910322489A CN110032822B CN 110032822 B CN110032822 B CN 110032822B CN 201910322489 A CN201910322489 A CN 201910322489A CN 110032822 B CN110032822 B CN 110032822B
- Authority
- CN
- China
- Prior art keywords
- rri
- ptr
- temperature
- heat
- spent
- 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
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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- 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)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Monitoring And Testing Of Nuclear Reactors (AREA)
- Laminated Bodies (AREA)
Abstract
本发明公开了一种乏池失去部分冷却后温度及温升速率计算的分析方法。该乏池失去部分冷却后温度及温升速率计算的分析方法包括以下步骤:步骤一,创建PTR/RRI换热器的工作原理模型;步骤二,初态参数的计算;步骤三,终态参数的计算。该乏池失去部分冷却后温度及温升速率计算的分析方法通过对大修期间停运一列PTR/RRI热交换器后的乏池温度变化趋势进行计算,根据计算结果合理安排大修期间相关工作,防止由于乏池温度过高导致超出运行技术规范的情况。
Description
技术领域
本发明是关于乏燃料水池失冷(失去部分冷却或失去全部冷却)后的乏燃料水池温度和温升速率进行预测领域,特别是关于一种乏池失去部分冷却后温度及温升速率计算的分析方法。
背景技术
反应堆停堆后,由于燃料组件中裂变产物的衰变,乏燃料组件仍然有一定的热量散出,称之为衰变热。当乏燃料组件卸载至乏燃料水池后,衰变热的积累会导致水池中温度升高。为了避免温度升高超出运行技术规范的要求,需要投运两列PTR/RRI热交换器以带出衰变热,维持乏池温度处于要求的范围内。
乏池的温度与热交换器RRI侧的进口温度、乏燃料组件衰变热、RRI侧水流量、PTR侧水流量、热交换器的效率相关。
乏燃料组件卸载至乏池后的衰变热数据由技术部燃料管理科进行计算,并通过热平衡关系式计算并给出大修期间堆芯和乏燃料池衰变热报告,其中包含RRI热交换器和PTR泵运行关系图(见图1),通过该图并根据机组实际的停堆天数和RRI侧进口温度查询当前机组至少需要投运两列PTR/RRI热交换器或者一列PTR/RRI热交换器。同时该报告中还包括乏池完全失冷后的温度变化曲线及乏池沸腾时间。
在换料大修期间一般会投运两列PTR/RRI热交换器,但是在大修期间可能会由于某些工作需要停运一列,此时所有的衰变热需要由剩余的一列热交换器带出,这种情况下会导致乏池温度的升高,在当前的电厂管理和技术规程中均没有对于停运一列PTR/RRI热交换器后的乏池温度及温升速率变化的计算方法进行规定。
另外,对于大修期间衰变热报告中计算通常采用的均为换热器的名义参数(换热器效率、PTR侧流量、RRI侧流量等),在机组的运行过程中,这些参数可能会发生改变,导致最终计算的结果与实际也会产生一定的偏差。
公开于该背景技术部分的信息仅仅旨在增加对本发明的总体背景的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域一般技术人员所公知的现有技术。
发明内容
本发明的目的在于提供一种步骤简单合理的乏池失去部分冷却后温度及温升速率计算的分析方法,该乏池失去部分冷却后温度及温升速率计算的分析方法通过对大修期间停运一列PTR/RRI热交换器后的乏池温度变化趋势进行计算,根据计算结果合理安排大修期间相关工作,防止由于乏池温度过高导致超出运行技术规范的情况。
为实现上述目的,本发明提供了一种乏池失去部分冷却后温度及温升速率计算的分析方法,包括以下步骤:
步骤一,创建PTR/RRI换热器的工作原理模型;
步骤二,初态参数的计算;
步骤三,终态参数的计算。
在一优选的实施方式中,步骤一中PTR/RRI换热器工作原理为:
通过传热公式根据RRI和PTR侧水流量、比热容计算在当前工况下传递的热量:
RRI侧带走的热量WRRI=(TRRI.out.e-TRRI.in.e)×C×QRRI/3.6/106MW (1)
PTR侧带走的热量WPTR=(TPTR.in.e-TPTR.out.e)×C×QRRI/3.6/106MW (2)总的换热量为We=(WPTR+WRRI)/2 (3)
其中,C为水的比热容,RRI侧进口温度为TRRI.in.e,出口温度为TRRI.out.e,进出口温差为ΔTRRI.e,RRI侧流量为QRRI,PTR侧进口温度为TPTR.in.e,出口温度为TPTR.out.e,进出口温差为ΔTPTR.e,PTR侧流量为QPTR,PTR侧平均温度与RRI侧平均温度偏差为ΔTPTR/RRI.e。
在一优选的实施方式中,步骤二中:
当预期停运一列换热器时乏池的衰变热为W0时:
对于RRI侧
We=k*ΔTRRI.e (4)
W0=k*ΔTRRI.0 (5)
得到RRI侧进出口温差ΔTRRI.0=ΔTRRI.e×W0/We (6)
同理,PTR侧进出口温差ΔTPTR.0=ΔTPTR.e×W0/We (7)
PTR侧平均温度与RRI侧平均温度差ΔTPTR/RRI.0=ΔTPTR/RRI.e×W0/We (8)
由9、10、11方程式得到换热器传递热量与PTR侧进口温度之间的关系
当两列热交换器运行时,每一列热交换器带出的热量为W0/2,此时乏池温度为:
在一优选的实施方式中,步骤三中:
当其中一列热交换器停运后,稳定工况下所有的衰变热将最终由剩余的一列全部带出,那么乏燃料水池最终的温度将会变化为:
对应乏池温度最终升高:
假设在停运一列热交换器后乏池的温升速率与时间t的关系为:
根据温升公式W=C×m×ΔT,通过积分运算最终得到
k1=ΔT·k2
其中m为乏池水装量,单位kg;C为水的比热容,单位J/kg/℃
对温升速率进行积分运算得到乏池的温度与时间t的关系为
其中TPTR.in.1为停运一列热交换器前的乏池温度。
在一优选的实施方式中,步骤一中根据RRI和PTR侧水流量、比热容计算在当前工况下传递的热量的传热公式为对数平均法。
与现有技术相比,根据本发明的乏池失去部分冷却后温度及温升速率计算的分析方法具有如下有益效果:该乏池失去部分冷却后温度及温升速率计算的分析方法通过对大修期间停运一列PTR/RRI热交换器后的乏池温度变化趋势进行计算,根据计算结果合理安排大修期间相关工作,防止由于乏池温度过高导致超出运行技术规范的情况。
附图说明
图1是乏燃料水池冷却——RRI热交换器和PTR泵运行关系图。
图2是根据本发明第一实施方式的乏池失去部分冷却后温度及温升速率计算的分析方法的PTR/RRI换热器的工作模型图。
图3是根据本发明第一实施方式的乏池失去部分冷却后温度及温升速率计算的分析方法的乏池温度和温升速率计算流程图。
图4是根据本发明第一实施方式的乏池失去部分冷却后温度及温升速率计算的分析方法的乏池温度和温升速率变化示意图。
具体实施方式
下面结合附图,对本发明的具体实施方式进行详细描述,但应当理解本发明的保护范围并不受具体实施方式的限制。
除非另有其它明确表示,否则在整个说明书和权利要求书中,术语“包括”或其变换如“包含”或“包括有”等等将被理解为包括所陈述的元件或组成部分,而并未排除其它元件或其它组成部分。
如图3和图4所示,根据本发明优选实施例的乏池失去部分冷却后温度及温升速率计算的分析方法包括以下具体步骤:
步骤一,创建PTR/RRI换热器的工作原理模型;
如图2所示,其中,RRI侧进口温度为TRRI.in.e,出口温度为TRRI.out.e,进出口温差为ΔTRRI.e,RRI侧流量为QRRI,PTR侧进口温度为TPTR.in.e,出口温度为TPTR.out.e,进出口温差为ΔTPTR.e,PTR侧流量为QPTR,PTR侧平均温度与RRI侧平均温度偏差为ΔTPTR/RRI.e。
通过传热公式(对数平均法)根据RRI和PTR侧水流量、比热容可以计算在当前工况下传递的热量:
RRI侧带走的热量WRRI=(TRRI.out.e-TRRI.in.e)×C×QRRI/3.6/106MW (1)
PTR侧带走的热量WPTR=(TPTR.in.e-TPTR.out.e)×C×QRRI/3.6/106MW (2)
总的换热量为We=(WPTR+WRRI)/2 (3)
其中C为水的比热容。
对于现场实际运行工况下,由于温度变化较小,可以认为换热器的效率不发生变化,所以换热器传递的热量与换热器进出口温差可以简单看做线性关系:W=k*ΔT,其中k为换热器效率。
RRI侧进出口温度、温差,PTR侧进出口温度、温差,换热量等参数可以采用换热器的名义参数进行计算,也可以采用实际稳定运行工况下的实际运行参数进行计算。
由于乏燃料组件的衰变热随时间呈近乎指数衰减变化,所以需要考虑预期停运一列换热器时的乏燃料组件衰变热,可以通过专用软件计算得到。
步骤二,初态参数的计算;
当预期停运一列换热器时乏池的衰变热为W0时:
对于RRI侧
We=k*ΔTRRI.e (4)
W0=k*ΔTRRI.0 (5)
可以得到RRI侧进出口温差ΔTRRI.0=ΔTRRI.e×W0/We (6)
同理,PTR侧进出口温差ΔTPTR.0=ΔTPTR.e×W0/We (7)
PTR侧平均温度与RRI侧平均温度差ΔTPTR/RRI.0=ΔTPTR/RRI.e×W0/We (8)
由9、10、11方程式可以得到换热器传递热量与PTR侧进口温度(乏池温度)之间的关系
当两列热交换器(认为两列热交换器的参数完全一样)运行时,每一列热交换器带出的热量为W0/2,此时乏池温度为:
步骤三,终态参数的计算;
当其中一列热交换器停运后,稳定工况下所有的衰变热将最终由剩余的一列全部带出,那么乏燃料水池最终的温度将会变化为:
对应乏池温度最终升高:
当停运一列热交换器,此时由于热量的积累,会导致乏池温度升高,由于乏池温度的升高,由PTR***通过换热器传递至RRI的热量将会增加,虽然此时换热器仍不足以带出全部的热量,但相比初始状态温升速率有所降低,最终乏池温度会稳定在某一温度,该温度相比初始状态升高了ΔT。
由于乏池温升速率和乏池温度并不是线性关系,可以假设在停运一列热交换器后乏池的温升速率与时间t的关系为:
根据温升公式W=C×m×ΔT,通过积分运算最终得到
k1=ΔT·k2
其中m为乏池水装量,单位kg;C为水的比热容,单位J/kg/℃
对温升速率进行积分运算可以得到乏池的温度与时间t的关系为
其中TPTR.in.1为停运一列热交换器前的乏池温度。
综上,该乏池失去部分冷却后温度及温升速率计算的分析方法通过对大修期间停运一列PTR/RRI热交换器后的乏池温度变化趋势进行计算,根据计算结果合理安排大修期间相关工作,防止由于乏池温度过高导致超出运行技术规范的情况。本计算方法不仅仅针对核电厂的PTR/RRI换热器,对于各个行业的双列换热器停运其中一列(停运一列换热器后导致的温度变化对换热器本身的换热系数的影响可忽略不计时)后的最终温升速率和温度均适用。
前述对本发明的具体示例性实施方案的描述是为了说明和例证的目的。这些描述并非想将本发明限定为所公开的精确形式,并且很显然,根据上述教导,可以进行很多改变和变化。对示例性实施例进行选择和描述的目的在于解释本发明的特定原理及其实际应用,从而使得本领域的技术人员能够实现并利用本发明的各种不同的示例性实施方案以及各种不同的选择和改变。本发明的范围意在由权利要求书及其等同形式所限定。
Claims (3)
1.一种乏池失去部分冷却后温度及温升速率计算的分析方法,其特征在于,包括以下步骤:
步骤一,创建PTR/RRI换热器的工作原理模型;
步骤二,初态参数的计算;
步骤三,终态参数的计算;
其中,所述步骤一中PTR/RRI换热器工作原理为:
通过传热公式根据RRI和PTR侧水流量、比热容计算在当前工况下传递的热量:
RRI侧带走的热量WRRI=(TRRI.out.e-TRRI.in.e)×C×QRRI/3.6/106MW (1)
PTR侧带走的热量WPTR=(TPTR.in.e-TPTR.out.e)×C×QRRI/3.6/106MW (2)
总的换热量为We=(WPTR+WRRI)/2 (3)
其中,C为水的比热容,RRI侧进口温度为TRRI.in.e,出口温度为TRRI.out.e,进出口温差为ΔTRRI.e,RRI侧流量为QRRI,PTR侧进口温度为TPTR.in.e,出口温度为TPTR.out.e,进出口温差为ΔTPTR.e,PTR侧流量为QPTR,PTR侧平均温度与RRI侧平均温度偏差为ΔTPTR/RRI.e;
其中,所述步骤二中:
当预期停运一列换热器时乏池的衰变热为W0时:
对于RRI侧
We=k*ΔTRRI.e (4)
W0=k*ΔTRRI.0 (5)
得到RRI侧进出口温差ΔTRRI.0=ΔTRRI.e×W0/We (6)
同理,PTR侧进出口温差ΔTPTR.0=ΔTPTR.e×W0/We (7)
PTR侧平均温度与RRI侧平均温度差ΔTPTR/RRI.0=ΔTPTR/RRI.e×W0/We (8)
由9、10、11方程式得到换热器传递热量与PTR侧进口温度之间的关系
当两列热交换器运行时,每一列热交换器带出的热量为W0/2,此时乏池温度为:
2.如权利要求1所述的乏池失去部分冷却后温度及温升速率计算的分析方法,其特征在于,所述步骤三中:
当其中一列热交换器停运后,稳定工况下所有的衰变热将最终由剩余的一列全部带出,那么乏燃料水池最终的温度将会变化为:
对应乏池温度最终升高:
假设在停运一列热交换器后乏池的温升速率与时间t的关系为:
根据温升公式W=C×m×ΔT,通过积分运算最终得到
k1=ΔT·k2
其中m为乏池水装量,单位kg;C为水的比热容,单位J/kg/℃
对温升速率进行积分运算得到乏池的温度与时间t的关系为
其中TPTR.in.1为停运一列热交换器前的乏池温度。
3.如权利要求1所述的乏池失去部分冷却后温度及温升速率计算的分析方法,其特征在于,所述步骤一中根据RRI和PTR侧水流量、比热容计算在当前工况下传递的热量的传热公式为对数平均法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910322489.7A CN110032822B (zh) | 2019-04-22 | 2019-04-22 | 乏池失去部分冷却后温度及温升速率计算的分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910322489.7A CN110032822B (zh) | 2019-04-22 | 2019-04-22 | 乏池失去部分冷却后温度及温升速率计算的分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110032822A CN110032822A (zh) | 2019-07-19 |
CN110032822B true CN110032822B (zh) | 2023-09-01 |
Family
ID=67239459
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910322489.7A Active CN110032822B (zh) | 2019-04-22 | 2019-04-22 | 乏池失去部分冷却后温度及温升速率计算的分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110032822B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013035918A1 (ko) * | 2011-09-08 | 2013-03-14 | 한전원자력연료 주식회사 | 사용후핵연료 저장조 실시간 감시시스템 및 그 방법 |
JP2013092485A (ja) * | 2011-10-27 | 2013-05-16 | Japan Atomic Energy Agency | 高温ガス炉の冷却材喪失事故時の運転方法 |
CN103514312A (zh) * | 2012-06-28 | 2014-01-15 | 中国核电工程有限公司 | 一种压缩空气***丧失及后果的分析方法 |
CN103824603A (zh) * | 2014-03-10 | 2014-05-28 | 中国人民解放军陆军军官学院 | 一种驱动高温液态金属循环流动冷却堆内部件的方法 |
CN104102782A (zh) * | 2014-07-18 | 2014-10-15 | 上海电力学院 | 压水堆核电站反应堆温度rbf神经网络控制方法 |
CN206907495U (zh) * | 2017-07-21 | 2018-01-19 | 国家电投集团科学技术研究院有限公司 | 核反应堆 |
JP2018136226A (ja) * | 2017-02-22 | 2018-08-30 | 東京電力ホールディングス株式会社 | 原子炉リスク管理装置、原子炉リスク管理方法、及び原子炉リスク管理プログラム |
CN109140563A (zh) * | 2018-08-21 | 2019-01-04 | 成都市新明节能科技有限公司 | 一种基于云计算平台的电锅炉精确供暖***及方法 |
JP2019023847A (ja) * | 2017-07-24 | 2019-02-14 | Koa株式会社 | 基板温度演算方法、基板温度演算装置、基板温度算出プログラム |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9406409B2 (en) * | 2013-03-06 | 2016-08-02 | Nuscale Power, Llc | Managing nuclear reactor spent fuel rods |
CN107863168B (zh) * | 2017-11-15 | 2020-04-07 | 苏州热工研究院有限公司 | 一种核电机组功率运行失去最终热阱的后撤时间分析方法 |
CN108108529B (zh) * | 2017-12-01 | 2021-07-06 | 东方电气集团东方汽轮机有限公司 | 一种简便测量铸造界面换热系数的反算方法 |
-
2019
- 2019-04-22 CN CN201910322489.7A patent/CN110032822B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2013035918A1 (ko) * | 2011-09-08 | 2013-03-14 | 한전원자력연료 주식회사 | 사용후핵연료 저장조 실시간 감시시스템 및 그 방법 |
JP2013092485A (ja) * | 2011-10-27 | 2013-05-16 | Japan Atomic Energy Agency | 高温ガス炉の冷却材喪失事故時の運転方法 |
CN103514312A (zh) * | 2012-06-28 | 2014-01-15 | 中国核电工程有限公司 | 一种压缩空气***丧失及后果的分析方法 |
CN103824603A (zh) * | 2014-03-10 | 2014-05-28 | 中国人民解放军陆军军官学院 | 一种驱动高温液态金属循环流动冷却堆内部件的方法 |
CN104102782A (zh) * | 2014-07-18 | 2014-10-15 | 上海电力学院 | 压水堆核电站反应堆温度rbf神经网络控制方法 |
JP2018136226A (ja) * | 2017-02-22 | 2018-08-30 | 東京電力ホールディングス株式会社 | 原子炉リスク管理装置、原子炉リスク管理方法、及び原子炉リスク管理プログラム |
CN206907495U (zh) * | 2017-07-21 | 2018-01-19 | 国家电投集团科学技术研究院有限公司 | 核反应堆 |
JP2019023847A (ja) * | 2017-07-24 | 2019-02-14 | Koa株式会社 | 基板温度演算方法、基板温度演算装置、基板温度算出プログラム |
CN109140563A (zh) * | 2018-08-21 | 2019-01-04 | 成都市新明节能科技有限公司 | 一种基于云计算平台的电锅炉精确供暖***及方法 |
Non-Patent Citations (1)
Title |
---|
《乏燃料池丧失冷却后燃料操作区域》;王国栋等;《核动力工程》;20161231;第105-106页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110032822A (zh) | 2019-07-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6826868B2 (ja) | 原子力プラントの評価システム及び方法 | |
EP1775732B1 (en) | A method of estimating dryout properties in a nuclear light water reactor | |
CN110032822B (zh) | 乏池失去部分冷却后温度及温升速率计算的分析方法 | |
EP1378916A2 (en) | Structured approach for risk-informing deterministic safety analyses | |
Sierchuła | Analysis of passive residual heat removal system in AP1000 nuclear power plant | |
Kumar et al. | Performance modeling of furnace draft air cycle in a thermal power plant | |
EP2071581A2 (en) | Method and apparatus for determination of safety limit minimum critical power ratio for a nuclear fuel core | |
Skalozubov et al. | Analysis of nuclear safety in diversification of Westinghouse fuel assemblies at WWER-1000 | |
Tewari et al. | Performance evaluation and availability analysis of steam generating system in a thermal power plant | |
JP3203921B2 (ja) | 熱交換器冷却系監視制御装置 | |
JP5122071B2 (ja) | 電力供給計画作成システム及び電力供給計画の作成方法 | |
Ma et al. | System Study: Isolation Condenser 1998-2022 | |
EP2048670B1 (en) | Methods and devices relating to a nuclear light water reactor of the boiling water kind | |
Deswandri et al. | Risk identification of integral pressurized water reactor (IPWR) cooling system using a combination HAZOP, FMEA, and FTA methods | |
CN115095994B (zh) | 一种热泵机组的运行控制方法、热泵机组及存储介质 | |
Kichline et al. | Use of expert judgment to support human reliability analysis of implementing FLEX equipment | |
Schroeder | System Study: Isolation Condenser | |
Juarsa et al. | Backward phenomenon on natural circulation flow based on power differences in FASSIP-01 loop | |
Warinner et al. | SASSYS-1 computer code verification with EBR-II test data | |
Alatrash et al. | Analyses for primary coolant pump coastdown phenomena for Jordan Research and Training Reactor | |
Griveau et al. | RELAP-5 Loss of Forced Cooling (LOFC) Transient Response Modeling for the PB-AHTR | |
Скалозубов et al. | Alternative method of strength analysis for a heat-power equipment under cyclic loads | |
Chang et al. | Sensitivity Study on LOHS Analysis for a demonstration SFR | |
CN118158984A (zh) | 一种数据机房液冷***的控制方法、介质及*** | |
Davis et al. | Modeling the GFR with RELAP5-3D |
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 |