CN105527452B - 基于同位素标记的堵水截源后底板承压水量的测定方法 - Google Patents

基于同位素标记的堵水截源后底板承压水量的测定方法 Download PDF

Info

Publication number
CN105527452B
CN105527452B CN201510969435.1A CN201510969435A CN105527452B CN 105527452 B CN105527452 B CN 105527452B CN 201510969435 A CN201510969435 A CN 201510969435A CN 105527452 B CN105527452 B CN 105527452B
Authority
CN
China
Prior art keywords
water
mrow
bottom plate
curtain
hole
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
CN201510969435.1A
Other languages
English (en)
Other versions
CN105527452A (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.)
Shandong University of Science and Technology
Original Assignee
Shandong University of Science and Technology
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 Shandong University of Science and Technology filed Critical Shandong University of Science and Technology
Priority to CN201510969435.1A priority Critical patent/CN105527452B/zh
Publication of CN105527452A publication Critical patent/CN105527452A/zh
Application granted granted Critical
Publication of CN105527452B publication Critical patent/CN105527452B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • 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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
    • G01F1/704—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow using marked regions or existing inhomogeneities within the fluid stream, e.g. statistically occurring variations in a fluid parameter
    • 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/001—Full-field flow measurement, e.g. determining flow velocity and direction in a whole region at the same time, flow visualisation

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Fluid Mechanics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于同位素标记的堵水截源后底板承压水量的测定方法，该方法首先通过对采煤工作面底板水文地质条件进行精细探测；然后在井下施工在煤层中布置帷幕截流巷，在帷幕截流巷处布置帷幕线，并依次对注浆孔进行帷幕注浆；利用工作面上下顺槽，每隔一定的距离向所述煤层中打向下的倾斜钻孔作为投放孔，相邻的投放孔之间打观测孔；最后通过公式计算得出底板承压水量和次要水源对底板承压水的补给速度。本发明对受承压水威胁的矿井采用堵水截源后人工排水时选用排水设备及所需费用具有重要指导意义，并且具有安全性高、可靠性好，施工方法简单、成本低等特点，具有较好的技术推广价值。

Description

基于同位素标记的堵水截源后底板承压水量的测定方法
技术领域
本发明涉及煤矿井下承压水量的测定方法，具体涉及一种基于同位素标记的堵水截源后底板承压水量的测定方法。
背景技术
随着煤矿开采深度的不断增加，开采条件、地质条件越来越复杂，工作面回采过程中，往往会遇到承压水的问题，这给安全生产带来很大困难，很容易出现底板突水事故。
目前，本领域的研究者们通常是采用深降强排和带压开采的方式来解决上述问题，其中的深降强排过程中一般会采用堵水截源后将底板承压水排出，这样减少了承压水的排水量，然而在选择排水设备时，一般矿井选择大功率的排水设备，没有考虑到承压水总量和补给速度，势必会造成对设备的浪费，而且大功率的设备耗电也多。
发明内容
本发明的目的在于提供一种基于同位素标记的堵水截源后底板承压水量的测定方法，该方法可以测出承压水的总量和次要水源对承压水的补给速度，具有一定的推广价值。
其技术解决方案包括：
一种堵水截源后底板承压水量的测定方法，依次包括以下步骤：
a对采煤工作面底板水文地质条件进行精细探测；
b井下施工在煤层中布置帷幕截流巷，在所述帷幕截流巷处布置帷幕线，根据所述帷幕线确定注浆孔的位置，并依次对所述注浆孔进行帷幕注浆，其中，所述注浆孔依次穿过位于所述煤层下方的隔水层和承压水含水层，并伸入至所述承压水含水层下面的岩层；
c利用工作面上下顺槽，每隔一定的距离向所述煤层中打向下的倾斜钻孔作为投放孔，所述投放孔穿过所述隔水层，并伸入至所述承压水含水层，向每个投放孔中投放重氧水，并且每个投放孔中的重氧水密度相同且体积相同，记下每个投放孔中的重氧水密度δ；
d在相邻的投放孔之间打向下的倾斜钻孔作为观测孔，所述观测孔穿过所述隔水层，并伸入至所述承压水含水层，所述投放孔在所述承压水含水层中的深度与观测孔在承压水含水层中的深度相同；
e待重氧水与承压水含水层中的承压水混合均匀后，进行测量，测量次数从第一次、第二次…第n次，其中每隔1小时测量一次，每次测量的同位素密度作为δi，通过式(1)计算得出底板承压水量M，
式(1)中：m为投放同位素总量，m3；
δ为最初投放的同位素密度，g/m3；
δi为第i次测量样本中同位素密度，g/m3；
f运用重氧水质量守恒原则，通过式(2)计算出次要水源对底板承压水的补给速度V，
式(2)中：Mi为第i次测量所得承压水的总量，m3；
n为测量的次数，且n≥10。
作为本发明的一个优选方案，步骤a中，精细探测包括对底板承压水主要的涌水通道位置、补给水源的流速、流量及流向进行测定。
作为本发明的另一个优选方案，步骤b中，所述帷幕线与地下水流方向垂直。
优选的，所述投放孔在所述承压水含水层中的深度为1～3m。
本发明所带来的有益技术效果是：
与现有技术相比，本发明通过向投放孔中放置同位素，并且与观测孔配合，通过同位素重氧水的质量守恒原则，借助公式计算得到底板承压水量和次要水源对承压水的补给速度，该方法在帷幕注浆和投放孔、观测孔实施过程中操作简单，所需人员少，一般2～4人即可完成，节约人员成本，并且整个操作过程均在巷道中进行，可以保证工作人员的安全，投放孔和观测孔长度较长，在实验完毕后封堵方便；
该方法借助公式计算得到底板承压水量和次要水源对承压水的补给速度，对受承压水威胁的矿井采用堵水截源后人工排水时选用排水设备及所需费用具有重要指导意义，根据此方法可以精确选取水泵的流量，避免盲目选取大排量的水泵，造成设备浪费，从而节约设备成本和降低水泵的能源消耗；并且在帷幕注浆和投放孔和观测孔实施过程中操作简单，所需人员少，一般2～4人即可完成，节约人员成本；
本发明相对于现有技术具有安全性与可靠性好，施工方法简单、效率高、成本低等有益效果。
附图说明
下面结合附图对本发明做进一步说明：
图1为本发明测定方法的流程图；
图2为帷幕注浆及钻孔布置平面图；
图3为帷幕注浆及钻孔布置剖面图；
图中，1-帷幕截流巷，2-承压水含水层，3-帷幕线，4-观测孔，5-投放孔，6-隔水层，7-岩层。
具体实施方式
下面结合具体实施例对本发明做详细说明。
如图1所示，本发明，基于同位素标记的堵水截源后底板承压水量的测定方法，依次包括以下步骤：
步骤1、收集地质资料，对采煤工作面水文地质条件做精细探测，具体精细探测的内容包括有：是否符合帷幕注浆的要求，查明底板承压水主要的涌水通道位置，并测定补给水源的流速、流量和流向；
步骤2、当上述精细探测达标后，帷幕注浆以实现堵水截源，具体的如图2所示，帷幕注浆切断地下水通道，包括以下子步骤：
井下施工并在煤层开挖帷幕截流巷1，布置帷幕线3，帷幕的位置应该符合以下条件：
(1.1)帷幕线必须布置在主要进水口地段并与地下水流方向直交；
(1.2)要求主要过水断面封堵严密，不产生绕流；
(1.3)做到帷幕形成后不被采矿所破坏；
为确保帷幕注浆达到堵水的效果，在注浆孔的布置上应该做到以下几点：
注浆孔布置的两端应坐落在不透水体上；
(2.1)注浆孔的布置密度应当合理选取，保证注浆孔浆液扩散后可以相互连接；
(2.2)注浆孔的钻进深度以穿透承压水含水层2为原则，即注浆孔一定要伸入承压水含水层2下面的岩层7，注浆孔在岩层7内的最宜深度为3～5m；
(2.3)注浆孔的布置以垂直水流方向为宜；
注浆钻孔完成后立即进行灌浆，形成完整的帷幕；
步骤3、利用工作面上下顺槽，每隔一定的距离向煤层中打向下的倾斜钻孔作为投放孔5，投放孔5有14个，如图2所示，分别编号为T1、T2、T3、T4、T5、T6、T7、T8、T9、T10、T11、T12、T13、T14，这14个投放孔5均穿过隔水层6，并伸入至承压水含水层2，向每个投放孔中投放重氧水，并且每个投放孔中的重氧水密度相同且体积相同，记下每个投放孔中的重氧水密度δ，并且根据投放孔中重氧水的密度和体积可计算得到投放的同位素总量；
步骤4、在相邻的投放孔之间打向下的倾斜钻孔作为观测孔4，分别编号G1、G2、G3、G4…G14，所有的观测孔4均穿过隔水层6，并伸入至承压水含水层2，投放孔5在承压水含水层2中的深度与观测孔在承压水含水层2中的深度相同；
步骤5、提取样本进行计算，计算出堵水截源后底板承压水量，具体如下：
提取承压水样本进行测量，测量分为两步：
(1)检验同位素与承压水是否混合均匀，
首先用质谱仪测出所有观测孔承压水中重氧水的密度，如果密度相差太大，则证明同位素没有混合均匀，然后相隔12小时候再次测量，直到样本中重氧水密度相近为止；
(2)提取样本进行计算，计算出堵水截源后底板承压水量，
证明同位素混合均匀后，第一次全面测量，以后每隔1小时测量一次，重复进行n次(n≥10)；
待重氧水与承压水含水层中的承压水混合均匀后，进行测量，测量次数从第一次、第二次…第n次，其中每隔1小时测量一次，每次测量的同位素密度作为δi，通过式(1)计算得出底板承压水量M，
式中：m为投放同位素总量，m3；
δ为最初投放的同位素密度，g/m3；
δi为第i次测量样本中同位素密度，g/m3；
运用重氧水质量守恒原则，通过式(2)计算出次要水源对底板承压水的补给速度V，
式(2)中：Mi为第i次测量所得承压水的总量，m3；
n为测量的次数，且n≥10。
实施例1：
以某煤矿某一工作面为例，在治理工作面底板承压水过程中，对该煤矿工作面底板堵水截源后底板承压水量的测定方法进行举例说明:
施工方法，包括以下步骤：
步骤一、收集地质资料，对采煤工作面水文地质条件做精细探测，
具体精细探测的内容包括有：
1、勘测此工作面地质是否符合帷幕注浆的要求，勘测结果表明工作面底板承压水的主要水源较为集中；
2、利用承压水示踪实验查明底板承压水主要的涌水通道位置，并测定补给水源的流速、流量和流向。
步骤二、当上述精细探测达标后，帷幕注浆以实现堵水截源，帷幕注浆切断地下水通道，包括以下子步骤：
(1)井下施工并在煤层开挖帷幕截流巷，然后布置帷幕线，帷幕线布置在主要进水口地段并与地下水流方向直交，帷幕线距离工作面停采线5～10m，保证帷幕形成后不被采矿过程所破坏；
(2)为确保帷幕注浆达到堵水的效果，注浆孔布置的两端坐落在上、下顺槽的底板岩层中，注浆孔的布置间距为1.5m，保证注浆孔浆液扩散后相互连接；注浆孔伸入承压水含水层下面的岩层5m；
注浆钻孔完成后立即进行灌浆，形成完整的帷幕；
步骤三、利用工作面上下顺槽，每隔25m向煤层中打向下的倾斜钻孔作为投放孔，投放孔共设14个，分别编号为T1、T2、T3、T4、T5、T6、T7、T8、T9、T10、T11、T12、T13、T14，向每个投放孔中投放密度为21.6％重氧水2500mm；
步骤四、在相邻的投放孔之间打向下的倾斜钻孔作为观测孔，共设14个，分别编号G1、G2、G3、G4…G14，投放孔在承压水含水层中的深度与观测孔在承压水含水层中的深度相同；
步骤五、提取样本进行计算，计算出堵水截源后底板承压水量，具体如下：
提取承压水样本进行测量，测量分为两步：
(1)检验同位素与承压水是否混合均匀，
首先用质谱仪测出所有观测孔承压水中重氧水的密度，第一次测量时，样本密度相差太大，方差为4.51，证明同位素没有混合均匀，然后相隔12小时候进行第二次测量，样本方差为0.19，证明同位素混合均匀。
(2)提取样本进行计算，计算出堵水截源后底板承压水量，
待重氧水与承压水含水层中的承压水混合均匀后，进行测量，测量次数从第1次、第2次…第15次，其中每隔1小时测量一次，每次测量的同位素密度作为δi，通过式(1)计算得出底板承压水量M，
式中：m为投放同位素总量，m3；
δ为最初投放的同位素密度，g/m3；
δi为第i次测量样本中同位素密度，g/m3；
运用重氧水质量守恒原则，通过式(2)计算出次要水源对底板承压水的补给速度V，
式(2)中：Mi为第i次测量所得承压水的总量，m3；
n为测量的次数，且n＝15。
测量结果为：此工作面承压水总量为421.7万t，次要水源对底板承压水的补给速度为26m3/h。据此方法选取水泵，此工作面第一年即为矿务集团节约设备成本105万元，电费15万元。
需要说明的是，在本说明书的教导下本领域技术人员所做出的任何等同方式，或明显变型方式均应在本发明的保护范围内。

Claims (4)

1.一种堵水截源后底板承压水量的测定方法，其特征在于，依次包括以下步骤：
a对采煤工作面底板水文地质条件进行精细探测；
b井下施工在煤层中布置帷幕截流巷，在所述帷幕截流巷处布置帷幕线，根据所述帷幕线确定注浆孔的位置，并依次对所述注浆孔进行帷幕注浆，其中，所述注浆孔依次穿过位于所述煤层下方的隔水层和承压水含水层，并伸入至所述承压水含水层下面的岩层；
c利用工作面上下顺槽，每隔一定的距离向所述煤层中打向下的倾斜钻孔作为投放孔，所述投放孔穿过所述隔水层，并伸入至所述承压水含水层，向每个投放孔中投放重氧水，并且每个投放孔中的重氧水密度相同且体积相同，记下每个投放孔中的重氧水密度δ；
d在相邻的投放孔之间打向下的倾斜钻孔作为观测孔，所述观测孔穿过所述隔水层，并伸入至所述承压水含水层，所述投放孔在所述承压水含水层中的深度与观测孔在承压水含水层中的深度相同；
e待重氧水与承压水含水层中的承压水混合均匀后，进行测量，测量次数从第一次、第二次…第n次，其中每隔1小时测量一次，每次测量的同位素密度作为δi，通过式(1)计算得出底板承压水量M，
<mrow> <mi>M</mi> <mo>=</mo> <mfrac> <mrow> <mi>m</mi> <mo>&amp;times;</mo> <mi>&amp;delta;</mi> </mrow> <msub> <mi>&amp;delta;</mi> <mi>i</mi> </msub> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
式(1)中：m为投放同位素总量，m3；
δ为最初投放的同位素密度，g/m3；
δi为第i次测量样本中同位素密度，g/m3；
f运用重氧水质量守恒原则，通过式(2)计算出次要水源对底板承压水的补给速度V，
<mrow> <mi>V</mi> <mo>=</mo> <mfrac> <mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>i</mi> <mo>=</mo> <mn>1</mn> </mrow> <mi>n</mi> </munderover> <mrow> <mo>(</mo> <msub> <mi>M</mi> <mrow> <mi>i</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>-</mo> <msub> <mi>M</mi> <mi>i</mi> </msub> <mo>)</mo> </mrow> </mrow> <mi>n</mi> </mfrac> <mo>-</mo> <mo>-</mo> <mo>-</mo> <mrow> <mo>(</mo> <mn>2</mn> <mo>)</mo> </mrow> <mo>;</mo> </mrow>
式(2)中：Mi为第i次测量所得承压水的总量，m3；
n为测量的次数，且n≥10。
2.根据权利要求1所述的堵水截源后底板承压水量的测定方法，其特征在于：步骤a中，精细探测包括对底板承压水主要的涌水通道位置、补给水源的流速、流量及流向进行测定。
3.根据权利要求1所述的堵水截源后底板承压水量的测定方法，其特征在于：步骤b中，所述帷幕线与地下水流方向垂直。
4.根据权利要求1所述的堵水截源后底板承压水量的测定方法，其特征在于：所述投放孔在所述承压水含水层中的深度为1～3m。
CN201510969435.1A 2015-12-21 2015-12-21 基于同位素标记的堵水截源后底板承压水量的测定方法 Active CN105527452B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510969435.1A CN105527452B (zh) 2015-12-21 2015-12-21 基于同位素标记的堵水截源后底板承压水量的测定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510969435.1A CN105527452B (zh) 2015-12-21 2015-12-21 基于同位素标记的堵水截源后底板承压水量的测定方法

Publications (2)

Publication Number Publication Date
CN105527452A CN105527452A (zh) 2016-04-27
CN105527452B true CN105527452B (zh) 2018-05-22

Family

ID=55769803

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510969435.1A Active CN105527452B (zh) 2015-12-21 2015-12-21 基于同位素标记的堵水截源后底板承压水量的测定方法

Country Status (1)

Country Link
CN (1) CN105527452B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7069812B2 (ja) * 2018-02-22 2022-05-18 株式会社島津製作所 分析用液監視装置
CN109059591A (zh) * 2018-07-19 2018-12-21 上海四埃美微科技有限公司 一种同步加热管及其制备方法和应用
CN111141928B (zh) * 2020-01-21 2021-09-21 中国矿业大学 一种示踪剂投放装置及其使用方法
CN111781121B (zh) * 2020-07-03 2023-01-31 四川大学 一种地铁地表沉陷预警系统
CN111929746B (zh) * 2020-07-28 2022-04-01 山东科技大学 一种基于示踪技术的探查陷落柱发育程度的检测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4963019A (en) * 1989-07-28 1990-10-16 Illinois State University Groundwater azimuth detection
CN101782591A (zh) * 2010-02-01 2010-07-21 河海大学 一种以温度为示踪剂的地下水流速流向探测方法及装置
CN104696010A (zh) * 2015-03-20 2015-06-10 青岛理工大学 一种矿山突水综合治理方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4963019A (en) * 1989-07-28 1990-10-16 Illinois State University Groundwater azimuth detection
CN101782591A (zh) * 2010-02-01 2010-07-21 河海大学 一种以温度为示踪剂的地下水流速流向探测方法及装置
CN104696010A (zh) * 2015-03-20 2015-06-10 青岛理工大学 一种矿山突水综合治理方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
"浅埋煤层开采矿井涌水量计算方法研究";刘洋;《西安科技大学学报》;20120930;第32卷(第5期);第565-568页 *
"矿井涌水量观测的几种简单方法";朱雷;《煤炭技术》;20070131;第26卷(第1期);第26-24页 *
"矿井涌水量计算方法评述";崔杰;《水力采煤与管道运输》;20091231(第4期);第1-4页 *
"矿井涌水量预测方法的发展";陈酩知 等;《工程地球物理学报》;20090228;第6卷(第1期);第68-71页 *

Also Published As

Publication number Publication date
CN105527452A (zh) 2016-04-27

Similar Documents

Publication Publication Date Title
CN105527452B (zh) 基于同位素标记的堵水截源后底板承压水量的测定方法
Sinclair et al. In situ leaching of copper: Challenges and future prospects
CN107083988B (zh) 基于千米定向钻孔的采动裂隙带瓦斯抽采技术方法
CN104481587B (zh) 大采深、大跨度综放采煤工作面顶板砂岩裂隙水探防方法
CN106194244B (zh) 低透气煤层液相co2相变致裂增透网格式瓦斯抽采方法
CN108331573B (zh) 采煤区的废弃油井治理方法
Sun et al. Multi-field coupling of water inrush channel formation in a deep mine with a buried fault
CN109767136A (zh) 煤层顶板含水层涌突水危险性评价与预测四双工作法
CN102799955B (zh) 突水系数小于0.06MPa/m区底板突水评价三图法
CN110645040A (zh) 一种基于定向钻机的突出煤层“五定”瓦斯治理方法
Tipping et al. Evidence for hydraulic heterogeneity and anisotropy in the mostly carbonate Prairie du Chien Group, southeastern Minnesota, USA
CN108757043A (zh) 一种针对采掘工作面防治水方法
CN109583047B (zh) 一种隧道突水突泥地表塌陷形成时间的预测方法
CN201526309U (zh) 地下水文观测孔
CN108561175B (zh) 煤岩巷关联的掘进巷道利用贯通钻孔进行瓦斯排放的方法
McCulloch et al. Determining cleat orientation of deeper coalbeds from overlying coals
CN113294203A (zh) 一种富水煤层岩巷安全掘进揭煤防治水保护方法
Koleini et al. Hypogenic karstification and conduit system controlling by tectonic pattern in foundation rocks of the Salman Farsi Dam in South-Western Iran
Davis et al. FACTORS INVOLVED IN EVALUATING GROUND WATER IMPACTS OF DEEP COAL MINE DRAINAGE 1
Naithani et al. Geological and geotechnical investigations of Loharinag-Pala hydroelectric project, Garhwal Himalaya, Uttarakhand
D'Andrea In situ leaching research in a copper deposit at the Emerald Isle mine
CN108196312B (zh) 一种水文地质中含水层与隔水层的识别方法
CN112901252A (zh) 深井高瓦斯压力揭煤区的地面钻井掏煤辅助消突方法
Aston et al. The effect of test cavity geology on the in situ permeability of coal measures strata associated with longwall mining
Jeffrey Hydraulic Fracture Geometry from Mineback Mapping

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant