CN113447074B - 基于分布式光纤技术的高寒区井壁变形状态辨识方法 - Google Patents
基于分布式光纤技术的高寒区井壁变形状态辨识方法 Download PDFInfo
- Publication number
- CN113447074B CN113447074B CN202110642814.5A CN202110642814A CN113447074B CN 113447074 B CN113447074 B CN 113447074B CN 202110642814 A CN202110642814 A CN 202110642814A CN 113447074 B CN113447074 B CN 113447074B
- Authority
- CN
- China
- Prior art keywords
- temperature
- wall
- sensing
- optical fiber
- well wall
- 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
- 239000013307 optical fiber Substances 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000005516 engineering process Methods 0.000 title claims abstract description 11
- 239000002689 soil Substances 0.000 claims abstract description 37
- 239000004567 concrete Substances 0.000 claims abstract description 27
- 230000008014 freezing Effects 0.000 claims abstract description 12
- 238000007710 freezing Methods 0.000 claims abstract description 12
- 238000011161 development Methods 0.000 claims abstract description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 9
- 239000002002 slurry Substances 0.000 claims abstract description 4
- 239000000835 fiber Substances 0.000 claims description 57
- 239000011435 rock Substances 0.000 claims description 32
- 238000010438 heat treatment Methods 0.000 claims description 12
- 238000012544 monitoring process Methods 0.000 claims description 10
- 229920000049 Carbon (fiber) Polymers 0.000 claims description 9
- 239000004917 carbon fiber Substances 0.000 claims description 9
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 6
- 239000000463 material Substances 0.000 claims description 6
- 230000003287 optical effect Effects 0.000 claims description 6
- 208000010392 Bone Fractures Diseases 0.000 claims description 4
- 206010017076 Fracture Diseases 0.000 claims description 4
- 239000003673 groundwater Substances 0.000 claims description 3
- 239000003085 diluting agent Substances 0.000 claims description 2
- 239000003822 epoxy resin Substances 0.000 claims description 2
- 229920000647 polyepoxide Polymers 0.000 claims description 2
- 239000011230 binding agent Substances 0.000 claims 1
- 230000005540 biological transmission Effects 0.000 claims 1
- 238000005070 sampling Methods 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 5
- 238000004062 sedimentation Methods 0.000 abstract 1
- 239000003245 coal Substances 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 2
- 238000000691 measurement method Methods 0.000 description 2
- 239000004033 plastic Substances 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000000945 filler Substances 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004441 surface measurement Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01D—MEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
- G01D21/00—Measuring or testing not otherwise provided for
- G01D21/02—Measuring two or more variables by means not covered by a single other subclass
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21F—SAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
- E21F17/00—Methods or devices for use in mines or tunnels, not covered elsewhere
- E21F17/18—Special adaptations of signalling or alarm devices
Landscapes
- Engineering & Computer Science (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geology (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
一种基于分布式光纤技术的高寒区井壁变形状态辨识方法,利用混凝土井壁沿线布设的应变传感光纤:计算井壁与岩土体的垂直沉降量,确定岩土体含水层变化引起的井壁所受附加应力,根据井壁内外温度差异,建立井壁内、外壁温度与切向应力之间关系,确定井壁受力特征,结合井壁与岩土体温度的分布式测量,确定冻结深度;建立传感光纤所受剪应力与应变梯度之间关系,确定井壁裂隙发育位置;当井壁受含水层水压、冻土压力及注浆压力作用下发生破裂时,地下水或注浆浆液沿裂隙流动并形成井壁温度的降低区,通过对温度降低区的辨识,确定裂隙发育位置,该方法具有分布式监测、辨识简便、精度高的特点。
Description
技术领域
本发明是一种煤矿井壁变形状态辨识方法,特别是基于分布式光纤技术的高寒区井壁变形状态辨识方法。
背景技术
井筒是煤矿生产的咽喉部位,井壁破裂和局部失稳影响煤炭企业的正常生产和井下煤矿工人的生命安全。有关井壁变形破坏状态辨识方面,主要采用几何测量法、传感器监测法及分布式光纤监测法。几何测量法,通过测量混凝土井壁变形量,确定井壁的稳定性,但该方法只限于表面量测、监测精度低。传感器监测法,主要采用电感式、振弦式及光纤光栅式应变和温度传感光纤,监测井壁变形破坏特征,该方法属于点式测量,存在漏检、高寒区中传感器适应性低的问题。分布式光纤监测法,还未形成井壁内部混凝土受力状况、裂隙分布和冻土层深度确定方法。
发明内容
针对上述现有技术存在的问题,本发明提供一种模型简单、易于操作的基于分布式光纤感测技术的高寒区立井井壁变形破坏状态辨识方法,通过井壁深度范围内应变和温度参量的连续量测,获得井壁冻结深度、混凝土受力状态和裂隙发育分布,及时掌握含水层水压、冻结压力及注浆压力等作用下井壁变形破坏状态,确保立井井壁的安全稳定。
本发明采用的技术方案是:一种基于分布式光纤技术的高寒区井壁变形状态辨识方法,具体方法如下:
步骤1)辨识井壁受力特征
1.1)根据井壁和周围岩土体中传感光纤所测的应变分布,分别计算井壁和周围岩土体的沉降量,确定周围岩土体对井壁的附加应力。
有关井壁和周围岩土体的沉降计算公式为
式中,S为基于分布式实测应变的井壁或井壁周围岩土体的沉降量;εm为实测应变;n为采用点数量;li和li+1分别为井壁或井壁周围岩土体中任意段埋深范围。
有关周围岩土体对井壁的附加应力为
式中,q为附加应力;E为井壁混凝土弹性模量;r为井壁的外径;εi和εi+1分别为井壁埋深分别为hi和hi+1的实测应变。
1.2)根据温度传感光纤测量井壁内外温度差异,建立井壁内、外壁温度与切向应力之间关系,确定井壁受力特征
1.3)利用混凝土井壁沿线布设的应变传感光纤,建立应变传感光纤所受剪应力与应变梯度之间关系,确定井壁裂隙发育位置,即
裂隙上、下应变传感光纤剪应力受力方向互为相向,且[σ]=E·ε大于井壁混凝土的抗拉强度σb时,井壁发生裂隙;
步骤2)当高寒区含水层冻胀压力、水压力及注浆压力等作用下井壁发生变形破坏时,利用井壁和周围岩土体沿线布设的应变传感光纤,计算岩土体沉降量对井壁的附加应力,确定井壁变形破坏位置。
步骤3)利用井壁沿线布设的温度传感光纤,当井壁受含水层水压、冻土压力及注浆压力作用下发生破裂时,地下水或注浆浆液沿裂隙流动并形成井壁温度的降低区,通过对温度降低区的辨识,确定裂隙发育位置。
进一步的,温度降低区的确定方法包括:1)基于井壁沿线布设的温度传感光纤,测量井壁温度分布,当测得温度低于零度时,确定井壁冻结深度;或
2)通过含有碳纤维加热丝的温度传感光纤,计算井壁沿线温度梯度,通过温度梯度的特征点,确定井壁冻结深度,温度梯度的计算公式为
TG为温度梯度;k为比例系数;T为温度传感光纤所测岩土体和混凝土温度;T0为通过碳纤维加热丝控制的温度传感光纤基础温度。
进一步的,分别在立井的井壁与岩土体边界、立井外壁与立井内壁之间及立井内壁表面周向布设多根垂向的传感光纤,传感光纤沿立井法向铺设;同时沿立井内壁布设环向的传感光纤,传感光纤沿立井深度设置多个,所有传感光纤以首尾相连方式熔接后,与预留在井口的光缆进行熔接并铺设至监测室,通过光缆与分布式光纤解调仪连接后,测量井壁沿线应变和温度参量。
更进一步的,所采用的传感光纤为一体式应变、温度传感光纤,应变传感光纤为紧套传感光纤,温度传感光纤为松套传感光纤,将温度传感光纤测量数据作为温度补偿,校准应变传感光纤量值。
更进一步的,所述传感光纤与护套之间添加碳纤维加热丝,加热丝连通电源之后可以调节传感光纤的温度值。
更进一步的,所述传感光纤可以设置为单独的应变光纤和温度光纤。
更进一步的,沿立井井壁开凹槽,凹槽的宽度为3mm、深度为5mm,传感光纤嵌装于凹槽内,待传感光纤布设后,采用环氧树脂与稀释剂组合的粘结剂进行充填。
与现有技术相比,本发明采用一种新颖的高寒区立井井壁应变和温度分布式光纤监测方法,通过井壁沿线设置的应变和温度传感光纤,探测出沿着光纤不同位置的应变和温度变化,精准确定井壁冻结深度、井壁内裂隙发育分布及井壁受力状态,实现立井井壁变形的实时判识,及时对井壁裂隙和混凝土强度较弱的位置进行注浆维护,避免灾害发生。该方法具有模型简单、全断面状态辨识及成本低廉等特点,可用于各种施工工艺下立井井壁的裂隙分布、冻结深度及混凝土壁强度的状态辨识,在本方法领域内具有广泛的操作性和实用性。
附图说明
图1是本发明的立井井壁分布式光纤监测方法布置示意图;
图2是是本发明图1的俯视示意图。
图中:1、立井内壁,2、立井外壁,3、光缆,4、一体式应变、温度传感光纤,5、岩土体。
具体实施方式
下面结合附图对本发明的实施例作进一步的描述:
实施步骤:如图1和图2所示,沿立井沿壁和环向开挖宽度为3mm、深度为5mm的凹槽,沿井壁竖向分别在井壁与岩土体5边界、立井外壁2与立井内壁1之间及立井内壁1表面布设一体式应变、温度传感光纤4,应变传感光纤为紧套传感光纤,温度传感光纤为松套传感光纤;同时沿立井内壁1环向位置布设类似的环向传感光纤;待所有传感光纤布设后,采用环氧树脂与稀释剂组合的粘结剂进行充填,确保光纤与井壁及岩土体5紧密结合,所使用的含有碳纤维加热丝的温度传感光纤结构为,在传感光纤与护套之间添加碳纤维加热丝,加热丝连通电源之后可以调节传感光纤的温度值,用于测定井壁的温度值与变形量,将温度传感光纤测量数据作为温度补偿,校准应变传感光纤量值;将所有传感光纤首尾相连方式熔接后,与预留在井口的光缆进行熔接并铺设至监测室,通过光缆与分布式光纤解调仪连接后,测量井壁应变和温度参量。
如图1所示,一种基于分布式光纤技术的高寒区井壁变形状态辨识方法,包括冻土带深度、井壁受力状态及裂隙分布确定方法,其具体方法如下:
步骤1)有关混凝土井壁受力特征辨识方法
1.1)根据井壁和周围岩土体中传感光纤所测的应变分布,分别计算井壁和周围岩土体的沉降量,确定周围岩土体对井壁的附加应力。
有关井壁和周围岩土体的沉降计算公式为
式中,S为基于分布式实测应变的井壁或井壁周围岩土体的沉降量;εm为实测应变;n为采用点数量;li和li+1分别为井壁或井壁周围岩土体中任意段埋深范围。
有关周围岩土体对井壁的附加应力为
式中,q为附加应力;E为井壁混凝土弹性模量;r为井壁的外径;εi和εi+1分别为井壁埋深分别为hi和hi+1的实测应变。
1.2)根据温度传感光纤测量井壁内外温度差异,建立井壁内外壁温度与切向应力之间关系,确定井壁受力特征,
1.3)利用混凝土井壁沿线布设的垂向传感光纤,建立传感光纤所受剪应力与应变梯度之间关系,确定井壁裂隙发育位置,即
裂隙上下传感光纤剪应力方向为互为相向,且[σ]=E·ε大于井壁混凝土的抗拉强度σb时,井壁发生裂隙。
裂隙发育位置的确定精度为0.05m,裂隙大小为0-3cm,针对混凝土井壁的受力状态,将井壁预警级别分为高、中、低,分别对应井壁混凝土受力的弹塑性状态、塑性状态和破裂状态。
步骤2)利用井壁沿线布设的温度传感光纤,当井壁受含水层水压、冻土压力及注浆压力作用下发生破裂时,地下水或注浆浆液沿裂隙流动并形成井壁温度的降低区,通过对温度降低区的辨识,确定裂隙发育位置。
有关冻结深度确定方法为:(1)基于井壁沿线布设的温度传感光纤,测量井壁温度分布,当分布式光纤所测的温度低于零度时,确定井壁冻结深度;
(2)通过含有碳纤维加热丝的温度传感光纤,计算井壁沿线温度梯度,通过温度梯度的特征点,确定井壁冻结深度。
Claims (5)
1.一种基于分布式光纤技术的高寒区井壁变形状态辨识方法,其特征在于,具体方法如下:
步骤1)辨识井壁受力特征
1.1)根据井壁和周围岩土体中传感光纤所测的应变分布,分别计算井壁和周围岩土体的沉降量,确定周围岩土体对井壁的附加应力
有关井壁和周围岩土体的沉降计算公式为:
式中,S为基于分布式实测应变的井壁或井壁周围岩土体的沉降量;εm为实测应变;n为采样点数量;li和li+1分别为井壁或井壁周围岩土体中任意段埋深范围;
有关周围岩土体对井壁的附加应力为:
式中,q为附加应力;E为井壁混凝土弹性模量;r为井壁的外径;εi和εi+1分别为井壁埋深分别为hi和hi+1的实测应变;
1.2)根据温度传感光纤测量井壁内外温度差异,建立井壁内、外壁温度与切向应力之间关系,确定井壁受力特征,计算公式为:
1.3)利用混凝土井壁沿线布设的应变传感光纤,建立应变传感光纤所受剪应力与应变梯度之间关系,确定井壁裂隙发育位置,即
裂隙上、下应变传感光纤剪应力受力方向互为相向,且[σ]=E·ε大于井壁混凝土的抗拉强度σb时,井壁发生裂隙;
步骤2)当高寒区含水层冻胀压力、水压力及注浆压力等作用下井壁发生变形破坏时,利用井壁和周围岩土体沿线布设的应变传感光纤,计算岩土体沉降量对井壁的附加应力,确定井壁变形破坏位置;
步骤3)利用井壁沿线布设的温度传感光纤,当井壁受含水层水压、冻土压力及注浆压力作用下发生破裂时,地下水或注浆浆液沿裂隙流动并形成井壁温度的降低区,通过对温度降低区的辨识,确定裂隙发育位置。
3.根据权利要求1或2所述的一种基于分布式光纤技术的高寒区井壁变形状态辨识方法,其特征在于,分别在立井周围岩土体、井壁与岩土体边界、立井外壁与立井内壁之间及立井内壁表面布设多根垂向传感光纤;同时沿立井内壁表面环向布设传感光纤,传感光纤沿立井深度设置多层;所有传感光纤以首尾相连方式熔接后,与预留在井口的光缆进行熔接并铺设至监测室,通过光缆与分布式光纤解调仪连接后,测量井壁沿线应变和温度参量;所述的传感光纤为一体式应变、温度传感光纤,将温度传感光纤测量数据作为温度补偿,校准应变传感光纤量值。
4.根据权利要求3所述的一种基于分布式光纤技术的高寒区井壁变形状态辨识方法,其特征在于,所述应变传感光纤为紧套传感光纤;温度传感光纤分为两种,一种为松套传感光纤;另一种为在松套传感光纤的传感光纤与护套之间添加碳纤维加热丝,加热丝连通电源之后能够调节传感光纤的温度值。
5.根据权利要求3所述的一种基于分布式光纤技术的高寒区井壁变形状态辨识方法,其特征在于,沿立井井壁开凹槽,凹槽的宽度为3mm、深度为5mm,传感光纤嵌装于凹槽内,待传感光纤布设后,采用环氧树脂与稀释剂组合的粘结剂进行充填。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110642814.5A CN113447074B (zh) | 2021-06-09 | 2021-06-09 | 基于分布式光纤技术的高寒区井壁变形状态辨识方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110642814.5A CN113447074B (zh) | 2021-06-09 | 2021-06-09 | 基于分布式光纤技术的高寒区井壁变形状态辨识方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113447074A CN113447074A (zh) | 2021-09-28 |
CN113447074B true CN113447074B (zh) | 2022-05-17 |
Family
ID=77811081
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110642814.5A Active CN113447074B (zh) | 2021-06-09 | 2021-06-09 | 基于分布式光纤技术的高寒区井壁变形状态辨识方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113447074B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114199679B (zh) * | 2021-12-09 | 2024-05-28 | 南京大学 | 基于光纤拉拔的冻土多物性参数的分布式原位测试方法 |
CN114111712B (zh) * | 2022-01-26 | 2022-04-26 | 天津风霖物联网科技有限公司 | 一种用于监测建筑物沉降的系统及方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103821507A (zh) * | 2014-03-18 | 2014-05-28 | 中国矿业大学 | 立井井壁变形分布式光纤检测方法 |
CN108489435A (zh) * | 2018-03-22 | 2018-09-04 | 安徽理工大学 | 对井壁进行监测的钢弦式传感器系统和变形预警的方法 |
CN108627186A (zh) * | 2018-03-22 | 2018-10-09 | 安徽理工大学 | 对井壁进行监测的光纤传感器系统与变形预警的方法 |
CN109356567A (zh) * | 2018-05-04 | 2019-02-19 | 中国石油集团海洋工程有限公司 | 深水浅部地层井壁稳定性预测方法 |
CN109779623A (zh) * | 2019-02-25 | 2019-05-21 | 中国标准化研究院 | 一种矿井筒监测方法 |
CN111382504A (zh) * | 2020-02-28 | 2020-07-07 | 中国矿业大学 | 一种煤层采动覆岩沉降状态辨识方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2412041A1 (en) * | 2000-06-29 | 2002-07-25 | Paulo S. Tubel | Method and system for monitoring smart structures utilizing distributed optical sensors |
-
2021
- 2021-06-09 CN CN202110642814.5A patent/CN113447074B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103821507A (zh) * | 2014-03-18 | 2014-05-28 | 中国矿业大学 | 立井井壁变形分布式光纤检测方法 |
CN108489435A (zh) * | 2018-03-22 | 2018-09-04 | 安徽理工大学 | 对井壁进行监测的钢弦式传感器系统和变形预警的方法 |
CN108627186A (zh) * | 2018-03-22 | 2018-10-09 | 安徽理工大学 | 对井壁进行监测的光纤传感器系统与变形预警的方法 |
CN109356567A (zh) * | 2018-05-04 | 2019-02-19 | 中国石油集团海洋工程有限公司 | 深水浅部地层井壁稳定性预测方法 |
CN109779623A (zh) * | 2019-02-25 | 2019-05-21 | 中国标准化研究院 | 一种矿井筒监测方法 |
CN111382504A (zh) * | 2020-02-28 | 2020-07-07 | 中国矿业大学 | 一种煤层采动覆岩沉降状态辨识方法 |
Non-Patent Citations (2)
Title |
---|
内层钢板高强钢筋混凝土复合井壁在冻结井筒应用研究;姚直书等;《采矿与安全工程学报》;20180715(第04期);第663-669页 * |
基于光纤应变测试技术的井筒壁后注浆井壁变形监测;郭建伟等;《煤矿安全》;20150320(第03期);第153-156页 * |
Also Published As
Publication number | Publication date |
---|---|
CN113447074A (zh) | 2021-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104315988B (zh) | 采动覆岩变形分布式光纤检测方法 | |
CN102878893B (zh) | 滑坡深部位移监测系统及方法 | |
CN111382504A (zh) | 一种煤层采动覆岩沉降状态辨识方法 | |
CN113447074B (zh) | 基于分布式光纤技术的高寒区井壁变形状态辨识方法 | |
CN104596405B (zh) | 地下雨污管道变形接触式实时监测装置与方法 | |
CN112504336A (zh) | 一种滑坡区管道变形监测系统 | |
CN103438820A (zh) | 一种钻孔剖面岩土体分层变形光纤测量方法 | |
CN111914373B (zh) | 长距离岩石顶管摩阻力计算方法及管岩接触状态检测方法 | |
CN103471747A (zh) | 一种非金属抗浮锚杆杆体应力测试装置 | |
CN204730824U (zh) | 一种分布式沉降测量装置 | |
CN110987039A (zh) | 一种工作面回采后煤层顶底板破坏深度测试方法 | |
CN104404986B (zh) | 一种岩质边坡潜在滑动面的定量识别方法 | |
CN102943493A (zh) | 测量预制桩内力及变形的方法 | |
CN107642111A (zh) | 一种端部敞口型phc管桩桩端阻力测试方法 | |
CN107631752A (zh) | 一种中线式尾矿库监测的装置和方法 | |
CN108225265A (zh) | 一种软土路基变形远程综合监测系统及安装方法 | |
CN104569351B (zh) | 一种用于相似模拟的光纤传感器铺设装置及铺设方法 | |
CN107941137B (zh) | 任意倾角钻孔变形量测方法 | |
CN205785182U (zh) | 一种基于光纤光栅传感技术的软土地基沉降传感器 | |
CN207228185U (zh) | 一种端部敞口型phc管桩桩端阻力测试装置 | |
CN104316029A (zh) | 一种地质沉降监测装置及监测方法 | |
CN107326942A (zh) | Phc管桩结构荷载检测装置及安装方法 | |
CN204154307U (zh) | 一种地质沉降监测装置 | |
CN106643649A (zh) | 一种土体深层沉降和孔隙水压力测量装置及测量方法 | |
CN109186445A (zh) | 无线监测炭质岩边坡表面变形的测试设备及其应用方法 |
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 |