CN112525774B - 一种基于涡街流量计频谱的流速、密度和粘度的测量方法 - Google Patents
一种基于涡街流量计频谱的流速、密度和粘度的测量方法 Download PDFInfo
- Publication number
- CN112525774B CN112525774B CN201910882747.7A CN201910882747A CN112525774B CN 112525774 B CN112525774 B CN 112525774B CN 201910882747 A CN201910882747 A CN 201910882747A CN 112525774 B CN112525774 B CN 112525774B
- Authority
- CN
- China
- Prior art keywords
- viscosity
- reynolds number
- density
- frequency spectrum
- vortex shedding
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N11/00—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties
- G01N11/02—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties by measuring flow of the material
- G01N11/04—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties by measuring flow of the material through a restricted passage, e.g. tube, aperture
- G01N11/08—Investigating flow properties of materials, e.g. viscosity, plasticity; Analysing materials by determining flow properties by measuring flow of the material through a restricted passage, e.g. tube, aperture by measuring pressure required to produce a known flow
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N9/00—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
- G01N9/26—Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity by measuring pressure differences
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P5/00—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft
- G01P5/14—Measuring speed of fluids, e.g. of air stream; Measuring speed of bodies relative to fluids, e.g. of ship, of aircraft by measuring differences of pressure in the fluid
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Aviation & Aerospace Engineering (AREA)
- Engineering & Computer Science (AREA)
- Measuring Volume Flow (AREA)
Abstract
本发明公开了一种基于涡街流量计频谱的流速、密度和粘度的测量方法,包含:步骤1、积累训练数据:通过实验或模拟测试涡街流量计的压力波动;步骤2、预处理数据与特征选择;步骤3、通过机器学习算法建立雷诺数和压力信号频谱的关系,然后通过雷诺数给出流速的非线性校正以及压力信号幅度与密度之间的关系;步骤4、预测:使用测试数据进行运行,以预测雷诺数的值;步骤5、根据雷诺数推导出粘度结果。本发明通过涡街流量计测量的压力波动的频谱以估计雷诺数,使用与时间相关的压力信号的形状来实现同时测量流速、密度和粘度;可以通过直接监测流过管道的流体混合物的粘度来提高测量精度,故障检测和预防性维护。
Description
技术领域
本发明涉及流量计测量领域,特别涉及一种基于涡街流量计频谱的流速、密度和粘度的测量方法。
背景技术
目前,往往需要通过流量计来测量工业环境中管道的流速,流量计包括差压、科里奥利、磁力、涡流(涡街流量计)、超声波和传热技术。其中,科里奥利、超声波、差压和涡流是最常见的。涡街流量计主要用于工业管道介质流体的流量测量,如气体、液体、蒸气等多种介质。其特点是压力损失小,量程范围大,精度高,在测量工况体积流量时几乎不受流体密度、压力、温度、粘度等参数的影响。无可动机械零件,因此可靠性高,维护量小。与科里奥利流量计和超声波流量计相比,涡街流量计更便宜、维护要求更低。
目前有些研究(Power-spectrum-analysis-based signal processing systemof vortex flowmeters–K.Xu,Y.Huang,X.Lv,IEEE trans.Instr.Meas.(2006)vol55no3pp1006-1011),仅仅关注一些信号分析,但是不未研究信号谐波与雷诺数的关系。还有研究(On the formation and suppression of vortex‘shedding’at low Reynolds numbers–P.Strykowski and K.Sreenivasan–J.Fluid.Mech.(1990)vol 218pp71-107)指出涡旋脱落频谱的实验和数值观察,它们的测量(热线风速测量)清楚地表明功率谱与雷诺数的函数关系。
还有研究(Secondary frequencies in the wake of a circular cylinderwith vortex shedding–S.Abarbanel,W.S.Don,D.Gottlieb,D.Rudy,C.Townsend–J.Fluid.Mech.(1991)vol 225pp 557-574)比较了几种数值方法,以阐明分谐波的出现条件和涡旋脱落频率的谐波,得出结论,虽然次谐波峰值的存在和幅度取决于外部因素,但谐波的存在仅取决于非流线体的几何形状和雷诺数。
这些论文涉及基本的流体力学,而不是专注于涡街流量计,也都没有涉及阻流体上的压力波动,更没有报道使用压力/速度波动的谱来找到雷诺数的方法。
另,专利申请(EP2072970B1:Method for determining the viscosity of afluid using a vortex flow meter)中记载,通过测量形成第一涡流所需的时间,使用涡旋流量计测量流体粘度,但是该方法受到严格限制,因为它仅允许在泵启动时测量流体的粘度,而不是在工业设置中的连续操作期间。
专利申请(DE102013019872A1:Method and apparatus for determining theviscosity of a flowing fluid in a flow channel)中,使用涡旋脱落来找到流体的粘度。当流体高速流动时,它会形成强烈的漩涡,从而改变系统的传热特性。通过热板和温度计网络监控这种传热,这种仪器需要了解流速和密度以确定粘度,或者需要了解粘度和密度以确定速度。
专利申请(DE102011050716B4 Method and apparatus for on-linemeasurement of viscosity of a fluid)使用涡旋流量计基于流速的知识确定流体粘度,以及还有专利申请(DE69811044T2 A method for measuring the viscosity of afluid)使用已知尺寸的通道上的压降和由涡旋流量计测量的流速来测量流体粘度,它需要比所提出的发明更多的传感器和更复杂的几何形状。此外它不能轻易地安装在现有的装置上。
涡街流量计的一个重要限制是它们测量流速但不提供流体密度和粘度的信息,但是密度对于确定流量的质量流速至关重要;而且粘度也很重要,因为有三个原因:(a)对于测量精度:涡街流量计略微非线性,这种非线性取决于雷诺数Re=ρDv/μ,其中ρ为流体密度,μ为粘度,D为钝体(也称非流线体)直径,v为流速。因此,知道粘度(和密度)会产生雷诺数,然后,可以校正测量值,这提高了速度测量精度。(b)对于故障检测:当使用制冷剂,蒸汽或可能发生相变的流体时,测量粘度可以在发生严重故障(堵塞,腐蚀等等)之前检测到意外的相变(凝胶,冷凝)。(c)用于诊断:即使使用已知成分的流体,监测粘度也可以在检测到超压时加速故障诊断;实际上,超压可能源于管堵塞或粘度增加(由于相或组成变化),当知道粘度时可以在不打开管的情况下立即确定原因。最常见的是,通过假设给定的流体成分并使用列表数据校正温度偏差来估计粘度,但是这种方法需要制表数据且这些数据可能并不总是可用,并且在组成方面不稳健。
基于上述原因,研发一种基于涡街流量计频谱的流速、密度和粘度的测量方法实为必要。
发明内容
本发明的目的在于提供一种基于涡街流量计频谱的流速、密度和粘度的测量方法,通过涡街流量计测量的压力波动的频谱以估计雷诺数,使用与时间相关的压力信号的形状来实现同时测量流速、密度和粘度;本发明可以通过直接监测流过管道的流体混合物的粘度来提高测量精度,故障检测和预防性维护。
为了达到上述目的,本发明通过以下技术方案实现:
一种基于涡街流量计频谱的流速、密度和粘度的测量方法,包含:步骤1、积累训练数据:通过实验或模拟测试涡街流量计的压力波动;步骤2、预处理数据与特征选择;所述特征选择是指将该压力信号频谱相关的量作为特征向量;步骤3、训练:通过机器学习算法建立雷诺数和压力信号频谱的关系,然后通过雷诺数给出流速的非线性校正以及压力信号幅度与密度之间的关系;步骤4、预测:使用测试数据进行运行,以预测雷诺数的值;步骤5、根据步骤4中的雷诺数推导出粘度结果。
优选地,所步骤S5中包含:
St=f·D/v=g(Re) (1);
式中,St为斯特罗哈尔数,当雷诺数Re为一定范围值时,St基本为一常数;f为压力波动频率;D为非流线体的体径;v为流速;当已知雷诺数Re,通过公式(1)算g的结果,再从结果g算流速v;
另外:P=ρ·v2·h(Re) (2)
其中,P为压力波动大小,ρ为流体的密度,当已知雷诺数Re,计算得出h,基于h和v,可以计算流体的密度ρ;
再根据:Re=ρDv/μ (3);
其中,μ为流体的粘度;当已知Re、ρ、D和v,可以算出流体粘度μ。
优选地,所述机器学习算法包含线性回归、最近邻、支持向量机回归、感知器方法以及它们之间的任意组合。
与现有技术相比,本发明的有益效果为:
本发明与现有技术中的仅使用压力幅度和波动频率的方法不同,本发明通过涡街流量计测量的压力波动的频谱以估计雷诺数,使用与时间相关的压力信号的形状来实现同时测量流速、密度和粘度;可以通过直接监测流过管道的流体混合物的粘度来提高测量精度,故障检测和预防性维护。
附图说明
图1为本发明的涡街流量计原理示意图;
图2为本发明的不同雷诺数下的模拟压力波动轨迹示意图;
图3为本发明的不同雷诺数下的压力信号频谱示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。
如图1所示,本发明公开了一种基于涡街流量计频谱的流速、密度和粘度的测量方法,其通过软件升级在典型的涡街流量计上实现。涡街流量计包括阻流体和连接到其上的压力传感器。压力传感器可以连接或不连接到阻流体。压力传感器用于检测由于阻流体的涡流脱落引起的流入相关的压力波动。
如图1-图3结合所示,本发明的基于涡街流量计频谱的流速、密度和粘度的测量方法,该方法包含以下步骤:
步骤1、积累训练数据:通过实验或模拟测试涡街流量计的压力波动;
步骤2、预处理数据与特征选择:
所述预处理数据的过程包含:
(i)对与时间相关的数据(例如压力波动)执行fft(快速傅里叶变换);
(ii)将压力波动频率进行归一化;例如图3所示的不同波动曲线的基波的无量纲频率都等于1,即将因雷诺数所引起的压力波动的基波的频率进行归一化;
(iii)取fft振幅的平方根(目标是放大小峰);
(iv)缩放:缩放0(最小振幅)和1(最大振幅)。
其中,该步骤(iii)还可以通过取fft振幅的对数进行进行预处理,然后在0(最小雷诺数)和1(最大雷诺数)之间进行缩放。
例如,压力信号频谱如图3所示,横坐标是压力波动频率,纵坐标是压力波动傅里叶变化的平方根。
fft不是获取信号频率内容信息的唯一方法,还可以采用其他任意函数进行预处理,例如小波算法、统计测量、递归神经网络和广义变换等,本发明在此不做限制。
所述特征选择是指将该压力信号频谱相关的量作为特征向量。
进一步地,将特征向量X链接到与雷诺数相关的量y,并通过参数优化计算y。其中,可以随便选X和y(只要它们与压力信号、雷诺数有关)。
本实施例,当前实现中,X=sqrt(F),F就是压力波动的傅里叶变化,y=log10(Re)。
步骤3、训练:通过机器学习算法建立雷诺数和压力信号频谱的关系,然后通过雷诺数给出流速的非线性校正以及压力信号幅度与密度之间的关系。其中,所述机器学习算包含例如线性回归(例如岭回归等)、最近邻、支持向量机回归、感知器和其他方法,以及它们之间的任意组合(堆叠或增强)。
步骤4、预测:上述模型(即步骤3中的流速的非线性校正以及压力信号幅度与密度之间的关系)使用新的测试数据进行运行,以预测雷诺数的值。其中,测试数据与训练数据类似地,必须进行预处理,并且必须对输出进行后处理(例如反向缩放和取幂),以恢复雷诺数。
步骤5、根据步骤4中的雷诺数推导出粘度结果,具体如下:
St=f·D/v=g(Re) (1);
式中,St为斯特罗哈尔数,当雷诺数Re为一定范围值时(例如500-150000时),St基本为一常数;f为压力波动频率;D为非流线体的体径;v为流速;因此,如果知道雷诺数Re,就可以算g(即g是通过公式1校准计算的结果),从结果g可以算流速v。
另外:P=ρ·v2·h(Re) (2)
其中,P为压力波动大小,ρ为流体的密度,因此,如果知道雷诺数Re,就可以算h(即h是通过校准计算的结果),从h和v可以算流体的密度ρ。
再根据:Re=ρDv/μ (3);
其中,μ为流体的粘度,因此,如果知道Re、ρ、D和v,就可以算流体粘度μ。
其中,在目前的实施中,岭回归显示出最佳表现,然后雷诺数给出了流速的非线性校正以及压力信号幅度和密度之间的关系。
如图2所示显示了涡街流量计在不同雷诺数下记录的模拟压力波动轨迹。其中,横坐标t为时间,纵坐标p为不同雷诺数下的压力信号。
根据p=(p-pmin)/(pmax-pmin),将压力信号按比例调整到0和1之间。在低雷诺数(例如Re=100)时,压力类似于正弦曲线;且在较高的雷诺数(例如Re=31623)下变得较不规则。
执行这些压力时间序列的傅立叶变换(F=傅立叶(p)),可以看出,当雷诺数增加时,尖峰(谐波)的数量也会增加:然后雷诺数给出了流速的非线性校正以及压力信号幅度和密度之间的关系;最终,粘度是从雷诺数的表达推导出来的(参考上述步骤S4)。
如表1所示为本发明的从模拟给出的预测与实际雷诺数的表格,该表格显示本发明的算法从Re=300到Re=11,000得到的结果是准确的。
表1本发明的从模拟给出的预测与实际雷诺数的表格
综上所述,本发明的基于涡街流量计频谱的流速、密度和粘度的测量方法,使用快速傅里叶变换算法提取信号频谱,并使用该频谱作为特征向量,同时使用神经网络、线性回归、支持向量等机器学习算法预测雷诺数,然后基于雷诺数给出了获取流速的非线性校正以及压力信号幅度与密度之间的关系,最终推导出粘度。
尽管本发明的内容已经通过上述优选实施例作了详细介绍,但应当认识到上述的描述不应被认为是对本发明的限制。在本领域技术人员阅读了上述内容后,对于本发明的多种修改和替代都将是显而易见的。因此,本发明的保护范围应由所附的权利要求来限定。
Claims (2)
1.一种基于涡街流量计频谱的流速、密度和粘度的测量方法,其特征在于,包含:
步骤1、积累训练数据:通过实验或模拟测试涡街流量计的压力波动;
步骤2、预处理数据,获取压力信号频谱内容,并进行特征选择;其中,所述特征选择是指将该压力信号频谱相关的量作为特征向量;
步骤3、训练:通过机器学习算法建立雷诺数和压力信号频谱的关系,然后通过雷诺数给出流速的非线性校正以及压力信号幅度与密度之间的关系;
步骤4、预测:基于流速的非线性校正以及压力信号幅度与密度之间的关系,使用测试数据进行运行,以预测雷诺数的值;
步骤5、根据步骤4中的雷诺数推导出粘度结果,包含:
St=f·D/v=g(Re) (1);
式中,St为斯特罗哈尔数,当雷诺数Re为一定范围值时,St基本为一常数;f为压力波动频率;D为非流线体的体径;v为流速;当已知雷诺数Re,通过公式(1)算g的结果,再从结果g算流速v;
另外:P=ρ·v2·h(Re) (2)
其中,P为压力波动大小,ρ为流体的密度,当已知雷诺数Re,计算得出h,基于h和v,可以计算流体的密度ρ;
再根据:Re=ρDv/μ (3);
其中,μ为流体的粘度;当已知Re、ρ、D和v,可以算出流体粘度μ。
2.如权利要求1所述的基于涡街流量计频谱的流速、密度和粘度的测量方法,其特征在于,
所述机器学习算法包含线性回归、最近邻、支持向量机回归、感知器方法以及它们之间的任意组合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910882747.7A CN112525774B (zh) | 2019-09-18 | 2019-09-18 | 一种基于涡街流量计频谱的流速、密度和粘度的测量方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910882747.7A CN112525774B (zh) | 2019-09-18 | 2019-09-18 | 一种基于涡街流量计频谱的流速、密度和粘度的测量方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112525774A CN112525774A (zh) | 2021-03-19 |
CN112525774B true CN112525774B (zh) | 2021-10-26 |
Family
ID=74975121
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910882747.7A Active CN112525774B (zh) | 2019-09-18 | 2019-09-18 | 一种基于涡街流量计频谱的流速、密度和粘度的测量方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112525774B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116106173B (zh) * | 2023-04-07 | 2023-06-27 | 天津索克汽车试验有限公司 | 一种新能源汽车用机油粘度测量方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5880377A (en) * | 1996-10-15 | 1999-03-09 | Lsi Logic Corporation | Method for low velocity measurement of fluid flow |
WO2004048906A2 (en) * | 2002-11-22 | 2004-06-10 | Cidra Corporation | Method for calibrating a flow meter having an array of sensors |
CN1595072A (zh) * | 2004-06-29 | 2005-03-16 | 浙江大学 | 涡街质量流量测量方法 |
CN105043511A (zh) * | 2015-07-24 | 2015-11-11 | 合肥科迈捷智能传感技术有限公司 | 基于雷诺数的涡街流量计系数修正方法 |
CN108351239A (zh) * | 2015-07-17 | 2018-07-31 | 恩德斯+豪斯流量技术股份有限公司 | 基于漩涡流量测量原理的流量测量装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AR052402A1 (es) * | 2005-07-18 | 2007-03-21 | Vignolo Gustavo Gabriel | Metodo de medida de velocidad en turbo-alimentadores de motores de combustion interna mediante analisis de vibraciones, sofware y equipo para realizarlo |
TW200825382A (en) * | 2006-12-15 | 2008-06-16 | Univ Nat Taiwan | Eddy current precision flow meter |
-
2019
- 2019-09-18 CN CN201910882747.7A patent/CN112525774B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5880377A (en) * | 1996-10-15 | 1999-03-09 | Lsi Logic Corporation | Method for low velocity measurement of fluid flow |
WO2004048906A2 (en) * | 2002-11-22 | 2004-06-10 | Cidra Corporation | Method for calibrating a flow meter having an array of sensors |
CN1595072A (zh) * | 2004-06-29 | 2005-03-16 | 浙江大学 | 涡街质量流量测量方法 |
CN108351239A (zh) * | 2015-07-17 | 2018-07-31 | 恩德斯+豪斯流量技术股份有限公司 | 基于漩涡流量测量原理的流量测量装置 |
CN105043511A (zh) * | 2015-07-24 | 2015-11-11 | 合肥科迈捷智能传感技术有限公司 | 基于雷诺数的涡街流量计系数修正方法 |
Non-Patent Citations (5)
Title |
---|
A facile pressure drop measurement system and its applications to gas-liquid microflow;Riaud, A et al;《MICROFLUIDICS AND NANOFLUIDICS》;20131130;715-724 * |
Analysis of Reynolds number scaling for viscous vortex reconnection;Qionglin Ni et al;《Physics of Fluids》;20121010;105102 /1-12 * |
Cyclones and attractive streaming generated by acoustical vortices;Antoine Riaud;《PHYSICAL REVIEW E》;20140714;013008/1-8 * |
双向内外管差压流量计的数值仿真与实验研究;吉鹏;《 中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》;20160815;C030-106 * |
基于涡街特性的流动分析与参数检测;孙志强;《 中国优秀博硕士学位论文全文数据库 (博士) 信息科技辑》;20070615;I140-63 * |
Also Published As
Publication number | Publication date |
---|---|
CN112525774A (zh) | 2021-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US20180128660A1 (en) | Flow rate determination method and apparatus | |
CN110197049B (zh) | 一种基于瞬变反问题的非金属管道泄漏定位方法 | |
CN110296327B (zh) | 一种基于瞬变流频率响应分析的管道泄漏检测方法 | |
JP2006500557A (ja) | 2相の流体の流れの検出および測定 | |
US20220163136A1 (en) | System and method for monitoring flow rate of regulating valve based on acoustic sensor | |
CN101208589A (zh) | 用于过程控制系统的诊断设备 | |
CN112525774B (zh) | 一种基于涡街流量计频谱的流速、密度和粘度的测量方法 | |
Ravula et al. | Experimental validation of leak and water-ingression detection in low-pressure gas pipeline using pressure and flow measurements | |
US11460818B2 (en) | Evaluation apparatus, evaluation system, and evaluation method | |
CN110186521B (zh) | 基于小波脊特征提取的涡街湿气过读补偿与流量测量方法 | |
CN109783972B (zh) | 基于流固耦合分析计算的止回阀内泄漏流量的监测方法 | |
CN113639987A (zh) | 一种基于瞬变流的阀门状态探测系统和方法 | |
Moazzeni et al. | Measurement and calibration of thermal cross-correlation-based flowmeter operating in harsh environment | |
Miau et al. | On measurement uncertainty of a vortex flowmeter | |
Brahma et al. | Data-based estimation and simulation of compressible pulsating flow with reverse-flow through an orifice | |
CN104808643B (zh) | 一种基于改进的倒双谱分析的控制回路非线性检测方法 | |
CN111521633B (zh) | 基于热扩散原理的气液多相流流型识别装置及方法 | |
Andrjesdottir et al. | Detection of fouling in a heat exchanger by application of a lock-in technique | |
CN114222918A (zh) | 从声信号识别液体流变特性 | |
CN103759894B (zh) | 一种换热器漏流的检测方法及系统 | |
EP3879236B1 (en) | Method for estimating a flow rate using a vortex flow meter | |
TWI633405B (zh) | 即時預檢填料塔液泛現象的方法 | |
Anabtawi et al. | Detection of blade contamination in turbine flowmeters using neural networks | |
Lindsay | Revision History | |
Haushofer et al. | Automated detection of loading and flooding points in packed columns |
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 |