CN110700884B - 一种基于解析法比拟预测综放工作面涌水量的方法 - Google Patents

一种基于解析法比拟预测综放工作面涌水量的方法 Download PDF

Info

Publication number
CN110700884B
CN110700884B CN201910950397.3A CN201910950397A CN110700884B CN 110700884 B CN110700884 B CN 110700884B CN 201910950397 A CN201910950397 A CN 201910950397A CN 110700884 B CN110700884 B CN 110700884B
Authority
CN
China
Prior art keywords
water
working surface
water inflow
working face
face
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
Application number
CN201910950397.3A
Other languages
English (en)
Other versions
CN110700884A (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 Coal Energy Research Institute Co Ltd
Original Assignee
China Coal Energy Research Institute 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 Coal Energy Research Institute Co Ltd filed Critical China Coal Energy Research Institute Co Ltd
Priority to CN201910950397.3A priority Critical patent/CN110700884B/zh
Publication of CN110700884A publication Critical patent/CN110700884A/zh
Application granted granted Critical
Publication of CN110700884B publication Critical patent/CN110700884B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21FSAFETY DEVICES, TRANSPORT, FILLING-UP, RESCUE, VENTILATION, OR DRAINING IN OR OF MINES OR TUNNELS
    • E21F17/00Methods or devices for use in mines or tunnels, not covered elsewhere
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells

Landscapes

  • Engineering & Computer Science (AREA)
  • Mining & Mineral Resources (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于解析法比拟预测综放工作面涌水量的方法,包括步骤:步骤一:获取预测工作面和比拟工作面的导水裂隙带高度;步骤二:确定预测工作面和比拟工作面的相关参数;步骤三:基于解析法确定预测工作面涌水量Q1和比拟工作面涌水量Q0之间的比例系数;步骤四:预测综放工作面的涌水量。本发明基于解析法利用综采工作面涌水量对综放工作面涌水量进行预测,考虑涌水量与充水含水层厚度之间的非线性关系,预测结果更准确,为工作面及矿井的排水系统设防提供合理的依据,节约了矿井基建投资费用,保障了矿井的安全高效生产。

Description

一种基于解析法比拟预测综放工作面涌水量的方法
技术领域
本发明属于煤炭开采水害防治技术领域,尤其涉及一种基于解析法比拟预测综放工作面涌水量的方法。
背景技术
随着我国煤炭资源战略重心的逐渐西移,西部侏罗系煤层的开采已经由浅部延伸到深部,顶板涌水量的精准预测预报面临更为严峻的考验。目前开采深部侏罗系煤层的矿井大多为新建矿井,且深部侏罗系煤层尚未大规模开采,矿区水文地质勘探程度普遍较低,所获取的水文地质参数十分有限。传统的涌水量预计主要有解析法、数值法、水文地质条件比拟法以及统计学方法等,但是在缺乏可靠的水文地质参数的情况下使得传统的预计方法往往与实际相差较大。预测值较大时造成排水系统的闲置;预测值较小时容易发生顶板水害事故。因此,深部侏罗系煤层综放开采的涌水量准确预计成为水害防治的关键。
发明内容
为了解决上述问题,本发明提供了一种基于解析法比拟预测综放工作面涌水量的方法,用于对深部侏罗系煤层综放开采的涌水量准确预计。
本发明通过以下技术方案来实现的:
一种基于解析法比拟预测综放工作面涌水量的方法,包括下述步骤:
步骤一:获取预测工作面和比拟工作面的导水裂隙带高度;
步骤二:获取预测工作面和比拟工作面的相关参数,该参数包括预测工作面的充水含水层厚度M1;比拟工作面的充水含水层厚度M0、涌水量Q0
步骤三:基于解析法确定预测工作面涌水量Q1和比拟工作面涌水量Q0之间的比例系数;
步骤四:预测综放工作面的涌水量。
本发明进一步的改进在于,步骤一中,预测工作面和比拟工作面的导水裂隙带高度获取过程如下:
根据《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规范》(2017)、井下或地面实测以及数值模拟等综合确定预测工作面和比拟工作面的导水裂隙带高度值。
本发明进一步的改进在于,步骤二中,预测工作面和比拟工作面的相关参数获取过程如下:
根据勘探期间钻孔的揭露资料,统计导水裂隙带高度范围内的充水含水层厚度,充水含水层厚度为细粒砂岩、中粒砂岩以及粗粒砂岩厚度的总和;对预测工作面和比拟工作面内部及工作面外围500m范围内的钻孔进行统计,并计算含水层厚度的平均值,将含水层厚度的平均值作为充水含水层的厚度,分别得到预测工作面的充水含水层厚度M1和比拟工作面的充水含水层厚度M0;根据比拟工作面回采期间不同回采进尺时所观测的采空区涌水量得到比拟工作面涌水量Q0
本发明进一步的改进在于,步骤三中,基于解析法中的集水廊道法确定预测工作面涌水量Q1和比拟工作面涌水量Q0之间的比例系数。
本发明进一步的改进在于,集水廊道法计算涌水量Q的公式如下:
Figure BDA0002225599880000021
Figure BDA0002225599880000022
其中,L为进水边界长度;K为渗透系数;H为水头高度;M为导水裂隙带范围内的含水层厚度;R为影响半径;r0为引用半径;S为工作面上方充水含水层水头高度下降量。
本发明进一步的改进在于,步骤三的具体过程如下:
预测工作面涌水量Q1和比拟工作面涌水量Q0的比值为比例系数λ,即
Figure BDA0002225599880000031
其中,λ为比例系数;Q1为预测工作面涌水量;Q0为比拟工作面涌水量;L1为预测工作面进水边界长度;L0为比拟工作面进水边界长度;K1为预测工作面渗透系数;K0为比拟工作面渗透系数;H1为预测工作面水头高度;H0为比拟工作面水头高度;M1为预测工作面含水层厚度;M0为比拟工作面含水层厚度;R1为预测工作面影响半径;R0为比拟工作面影响半径。
本发明进一步的改进在于,步骤四中,综放工作面涌水量的计算公式如下:
Q1=λQ0
本发明进一步的改进在于,所述预测方法适用于开采深部侏罗系煤层的矿井。
与现有技术相比,本发明至少具有如下有益的技术效果:
传统的比拟法预测工作面涌水量时,主要考虑工作面的开采面积F和含水层的水位降深S,然而传统的比拟法是基于地质与水文地质条件相似、开采方式基本相同的矿井,但开采方式不同时,上覆岩层的移动规律差别较大,导致充水含水层的厚度及充水条件都发生较大变化,使得传统的预测方法在开采方式不同时不具适用性。本发明提供的基于解析法比拟预测综放工作面涌水量的方法先分别确定预测工作面和比拟工作面的导水裂隙带高度,然后再基于解析法中的集水廊道法确定包含预测工作面和比拟工作面的水文地质信息,通过确定出的水文地质信息求取预测工作面和比拟工作面之间的比例系数λ,最后求得预测工作面的涌水量。通过上述的预测方法相比于传统的比拟法与实际更为接近,为矿井的排水系统设防和防治水措施提供更为精准的指导,保障了矿井安全高效回采。
附图说明
图1是本发明基于解析法比拟预测综放工作面涌水量的方法的流程图;
图2是西北某矿103工作面涌水量预测值与实际值对比图。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明。
如图1所示,本发明提供的基于解析法比拟预测综放工作面涌水量的方法,包括下述步骤:
步骤一确定预测工作面和比拟工作面的导水裂隙带高度。具体过程为:根据《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规范》(2017)、井下或地面实测以及数值模拟等综合确定预测工作面和比拟工作面的导水裂隙带高度值。
步骤二确定预测工作面的充水含水层厚度M1以及比拟工作面的充水含水层厚度M0、涌水量Q0。具体过程为:根据勘探期间钻孔的揭露资料,统计导水裂隙带高度范围内的充水含水层厚度,充水含水层厚度为细粒砂岩、中粒砂岩以及粗粒砂岩厚度的总和。对预测工作面和比拟工作面内部及工作面外围500m范围内的钻孔进行统计,并计算含水层厚度的平均值,将含水层厚度的平均值作为充水含水层的厚度,分别得到预测工作面的充水含水层厚度M1和比拟工作面的充水含水层厚度M0;根据比拟工作面回采期间不同回采进尺时所观测的采空区涌水量得到比拟工作面的涌水量Q0
步骤三确定预测工作面涌水量Q1和比拟工作面涌水量Q0之间的比例系数。具体过程为:首先,根据集水廊道法分别确定预测工作面涌水量Q1和比拟工作面涌水量Q0,集水廊道法计算公式如下:
Figure BDA0002225599880000041
Figure BDA0002225599880000042
其中,L为进水边界长度;K为渗透系数;H为水头高度;M为导水裂隙带范围内的含水层厚度;R为影响半径;r0为引用半径;S为工作面上方充水含水层水头高度下降量。
然后将预测工作面涌水量Q1和比拟工作面涌水量Q0的比值作为比例系数λ,即
Figure BDA0002225599880000051
其中,λ为比例系数;Q1为预测工作面涌水量;Q0为比拟工作面涌水量;L1为预测工作面进水边界长度;L0为比拟工作面进水边界长度;K1为预测工作面渗透系数;K0为比拟工作面渗透系数;H1为预测工作面水头高度;H0为比拟工作面水头高度;M1为预测工作面含水层厚度;M0为比拟工作面含水层厚度;R1为预测工作面影响半径;R0为比拟工作面影响半径。
步骤四预测综放工作面的涌水量Q1。具体过程为:采用如下方法
Figure BDA0002225599880000052
为了解决开采深部侏罗系煤层的矿井地质及水文地质条件相似而开采方式不同,顶板涌水量预计困难的问题。本发明提供的基于解析法比拟预测综放工作面涌水量的方法先分别确定预测工作面和比拟工作面的导水裂隙带高度,然后再基于解析法中的集水廊道法确定包含预测工作面和比拟工作面的水文地质信息,通过确定出的水文地质信息求取预测工作面和比拟工作面之间的比例系数λ,最后求得预测工作面的涌水量。通过上述的预测方法相比于传统的比拟法与实际更为接近,为矿井的排水系统设防和防治水措施提供更为精准的指导,保障了矿井安全高效回采。
实施例:
本实施例以开采深部侏罗系煤层某矿的103综放工作面涌水量预测为例,步骤如下:
步骤一:某矿的103工作面宽度300m,走向长度1100m,煤层厚度平均10m左右,开采方式为综放开采;比拟(201)工作面宽度300m,煤层平均厚度5.92m,开采方式为综采开采。
根据现场实测和数值模拟分析,比拟工作面采后裂采比为23.7;由多元线性回归和数值模拟分析得到预测工作面采后裂采比为25.4。
步骤二:根据勘探期间钻孔的实际揭露情况,分别统计预测工作面和比拟工作面导水裂隙带范围内的充水含水层厚度,得到预测工作面的充水含水层厚度M1=81m、比拟工作面的充水含水层厚度M0=52m;比拟工作面回采期间涌水量Q0随进尺的变化情况如下表:
Figure BDA0002225599880000061
步骤三:将预测工作面和比拟工作面相关参数代入集水廊道法公式:
Figure BDA0002225599880000062
Figure BDA0002225599880000063
其中,L为进水边界长度;K为渗透系数;H为水头高度;M为导水裂隙带范围内的含水层厚度;R为影响半径;r0为引用半径;S为工作面上方充水含水层水头高度下降量。
由于预测工作面和比拟工作面地质及水文地质条件相似,因此在回采进尺相同时公式中进水边界长度L、渗透系数K以及影响半径R相同,因此比例系数
Figure BDA0002225599880000064
预测工作面和比拟工作面充水含水层水头高度为550m,因此
Figure BDA0002225599880000065
步骤四:103工作面涌水量Q1=λQ0=1.51Q0,根据比拟工作面回采期间的涌水量Q0预测出103工作面涌水量Q1如下表:
Figure BDA0002225599880000071
103工作面回采结束后,根据矿井实际涌水量数据,绘制出涌水量预测值和实际涌水量的对比情况见说明书附图2。预测值与实际值吻合度较好,证明本方法准确率较高,运用本发明提供的预计方法,能够满足深部侏罗系煤层综放开采的涌水量预计工作。

Claims (3)

1.一种基于解析法比拟预测综放工作面涌水量的方法,其特征在于,包括下述步骤:
步骤一:获取预测工作面和比拟工作面的导水裂隙带高度;
步骤二:获取预测工作面和比拟工作面的相关参数,该参数包括预测工作面的充水含水层厚度M1;比拟工作面的充水含水层厚度M0、涌水量Q0观测;预测工作面和比拟工作面的相关参数获取过程如下:
根据勘探期间钻孔的揭露资料,统计导水裂隙带高度范围内的充水含水层厚度,充水含水层厚度为细粒砂岩、中粒砂岩以及粗粒砂岩厚度的总和;对预测工作面和比拟工作面内部及工作面外围500m范围内的钻孔进行统计,并计算含水层厚度的平均值,将含水层厚度的平均值作为充水含水层的厚度,分别得到预测工作面的充水含水层厚度M1和比拟工作面的充水含水层厚度M0;根据比拟工作面回采期间不同回采进尺时所观测的采空区涌水量得到比拟工作面涌水量Q0观测
步骤三:基于解析法中的集水廊道法确定预测工作面涌水量Q1解析和比拟工作面涌水量Q0解析之间的比例系数;集水廊道法计算涌水量Q的公式如下:
Figure FDA0002683066540000011
Figure FDA0002683066540000012
其中,L为进水边界长度;K为渗透系数;H为水头高度;M为导水裂隙带范围内的充水含水层厚度;R为影响半径;r0为引用半径;S为工作面上方充水含水层水头高度下降量;
具体过程如下:
预测工作面涌水量Q1解析和比拟工作面涌水量Q0解析的比值为比例系数λ,即
Figure FDA0002683066540000021
其中,λ为比例系数;Q1解析为预测工作面涌水量;Q0解析为比拟工作面涌水量;L1为预测工作面进水边界长度;L0为比拟工作面进水边界长度;K1为预测工作面渗透系数;K0为比拟工作面渗透系数;H1为预测工作面水头高度;H0为比拟工作面水头高度;M1为预测工作面充水含水层厚度;M0为比拟工作面充水含水层厚度;R1为预测工作面影响半径;R0为比拟工作面影响半径;
步骤四:预测综放工作面的涌水量,计算公式如下:
Q1预测=λQ0观测
2.根据权利要求1所述的一种基于解析法比拟预测综放工作面涌水量的方法,其特征在于,步骤一中,预测工作面和比拟工作面的导水裂隙带高度获取过程如下:
根据《建筑物、水体、铁路及主要井巷煤柱留设与压煤开采规范》(2017)、井下或地面实测以及数值模拟综合确定预测工作面和比拟工作面的导水裂隙带高度值。
3.根据权利要求1所述的一种基于解析法比拟预测综放工作面涌水量的方法,其特征在于,所述方法适用于开采深部侏罗系煤层的矿井。
CN201910950397.3A 2019-10-08 2019-10-08 一种基于解析法比拟预测综放工作面涌水量的方法 Active CN110700884B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910950397.3A CN110700884B (zh) 2019-10-08 2019-10-08 一种基于解析法比拟预测综放工作面涌水量的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910950397.3A CN110700884B (zh) 2019-10-08 2019-10-08 一种基于解析法比拟预测综放工作面涌水量的方法

Publications (2)

Publication Number Publication Date
CN110700884A CN110700884A (zh) 2020-01-17
CN110700884B true CN110700884B (zh) 2020-11-06

Family

ID=69198030

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910950397.3A Active CN110700884B (zh) 2019-10-08 2019-10-08 一种基于解析法比拟预测综放工作面涌水量的方法

Country Status (1)

Country Link
CN (1) CN110700884B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113901546A (zh) * 2021-09-26 2022-01-07 中铁第一勘察设计院集团有限公司 隧道涌水量计算方法及装置
CN114264680B (zh) * 2021-11-15 2023-06-13 中煤科工集团西安研究院有限公司 一种基于比拟法的矿井水中氟化物浓度预测方法

Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104200099A (zh) * 2014-09-01 2014-12-10 山东科技大学 基于水文地质台账的矿井涌水量计算方法
CN104200291A (zh) * 2014-09-28 2014-12-10 山东科技大学 一种基于小波变换和arma-svm的涌水量预测方法
CN105160079A (zh) * 2015-08-11 2015-12-16 西南石油大学 一种岩溶隧道岩溶安全厚度计算方法
CN106437844A (zh) * 2016-10-26 2017-02-22 重庆大学 一种超前预报隧道涌水位置的方法
CN107622164A (zh) * 2017-09-25 2018-01-23 中国煤炭地质总局勘查研究总院 一种矿井涌水量预测方法及装置
CN107748928A (zh) * 2017-10-12 2018-03-02 中国煤炭地质总局水文地质局 一种基于空间分析的矿井动涌水量计算方法及系统
CN207598293U (zh) * 2017-11-24 2018-07-10 兖州煤业股份有限公司 矿井工作面涌水监测、预测装置
CN108491986A (zh) * 2018-06-22 2018-09-04 中交公局第工程有限公司 一种隧道涌水量预测方法
CN108627436A (zh) * 2017-03-20 2018-10-09 中国海洋石油集团有限公司 一种基于施工动态监测数据的地下水封油库涌水量的预测方法
CN108955649A (zh) * 2018-05-21 2018-12-07 中国矿业大学 一种煤矿工作面顶板过程涌水量的预测方法
CN109165449A (zh) * 2018-08-28 2019-01-08 中国煤炭地质总局勘查研究总院 矿井涌水边界过水断面流量的计算方法及装置
CN109166633A (zh) * 2018-08-28 2019-01-08 中国煤炭地质总局勘查研究总院 矿井涌水量计算方法及装置
CN109190854A (zh) * 2018-10-29 2019-01-11 六盘水师范学院 一种底板充水多因素比拟预测工作面涌水量的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150001161A1 (en) * 2013-07-01 2015-01-01 Rockwater Resource, LLC Liquid treatment station including plural mobile units and methods for operation thereof

Patent Citations (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104200099A (zh) * 2014-09-01 2014-12-10 山东科技大学 基于水文地质台账的矿井涌水量计算方法
CN104200291A (zh) * 2014-09-28 2014-12-10 山东科技大学 一种基于小波变换和arma-svm的涌水量预测方法
CN105160079A (zh) * 2015-08-11 2015-12-16 西南石油大学 一种岩溶隧道岩溶安全厚度计算方法
CN106437844A (zh) * 2016-10-26 2017-02-22 重庆大学 一种超前预报隧道涌水位置的方法
CN108627436A (zh) * 2017-03-20 2018-10-09 中国海洋石油集团有限公司 一种基于施工动态监测数据的地下水封油库涌水量的预测方法
CN107622164A (zh) * 2017-09-25 2018-01-23 中国煤炭地质总局勘查研究总院 一种矿井涌水量预测方法及装置
CN107748928A (zh) * 2017-10-12 2018-03-02 中国煤炭地质总局水文地质局 一种基于空间分析的矿井动涌水量计算方法及系统
CN207598293U (zh) * 2017-11-24 2018-07-10 兖州煤业股份有限公司 矿井工作面涌水监测、预测装置
CN108955649A (zh) * 2018-05-21 2018-12-07 中国矿业大学 一种煤矿工作面顶板过程涌水量的预测方法
CN108491986A (zh) * 2018-06-22 2018-09-04 中交公局第工程有限公司 一种隧道涌水量预测方法
CN109165449A (zh) * 2018-08-28 2019-01-08 中国煤炭地质总局勘查研究总院 矿井涌水边界过水断面流量的计算方法及装置
CN109166633A (zh) * 2018-08-28 2019-01-08 中国煤炭地质总局勘查研究总院 矿井涌水量计算方法及装置
CN109190854A (zh) * 2018-10-29 2019-01-11 六盘水师范学院 一种底板充水多因素比拟预测工作面涌水量的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
鹤煤九矿涌水水源分析及涌水量预测;杜小冲;《能源与环保》;20181231;第69-72+76页 *

Also Published As

Publication number Publication date
CN110700884A (zh) 2020-01-17

Similar Documents

Publication Publication Date Title
CN111191849B (zh) 一种西部矿区深埋工作面涌水量预测方法
CN102041995B (zh) 复杂油藏水淹状况监测系统
CN107728228B (zh) 水体下采煤覆岩导水裂隙带发育高度预计方法
Cooley et al. Nonlinear‐regression groundwater flow modeling of a deep regional aquifer system
CN104863577B (zh) 利用地震纵波传播时间预测地层孔隙压力的方法
CN102747991A (zh) 一种确定合采井单层产量的方法
CN110359895B (zh) 一种非均质巨厚砂岩水平井分段压裂的探放水方法
CN112508330B (zh) 西部矿区开采扰动下矿井水来源判别方法
CN110700884B (zh) 一种基于解析法比拟预测综放工作面涌水量的方法
CN112132454B (zh) 一种煤层顶板或底板含水层富水性综合评价方法
CN105354365A (zh) 间接充水含水层突水危险性综合评价方法及系统
CN107895216B (zh) 一种快速圈定油气掩盖区古层间氧化前锋线位置的方法
CN114046178A (zh) 一种用于煤矿井下防治水预警系统、方法及应用
CN111415038A (zh) 一种多工作面采空区涌水量随采精细化预测方法
CN113375753B (zh) 一种煤矿综采工作面采动对地下水影响的监测与分析方法
Light et al. Mine-water flow between contiguous flooded underground coal mines with hydraulically compromised barriers
CN110851991B (zh) 一种地下水流数值模拟方法
CN111119992B (zh) 一种煤层顶板疏放水钻孔参数确定方法
Winters et al. Ground water flow parameterization of an Appalachian coal mine complex
CN115861551A (zh) 一种地浸采铀的数字建井方法
CN115907187A (zh) 一种大采高综放导水裂隙带发育高度预测方法
Feng et al. Research on Water-conducting Fractured Zone height under the condition of large mining height in Yushen mining area, China
Martinez et al. Use of numerical groundwater modelling for mine dewatering assessment
CN106503854A (zh) 综放开采煤层上覆岩体跨裂高度预测方法
CN107369103A (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