CN113885086B - 一种井下直流赤道向偶极动源异常自显式超前探测方法 - Google Patents
一种井下直流赤道向偶极动源异常自显式超前探测方法 Download PDFInfo
- Publication number
- CN113885086B CN113885086B CN202110896207.1A CN202110896207A CN113885086B CN 113885086 B CN113885086 B CN 113885086B CN 202110896207 A CN202110896207 A CN 202110896207A CN 113885086 B CN113885086 B CN 113885086B
- Authority
- CN
- China
- Prior art keywords
- distance
- electrode
- receiving
- electrodes
- tunneling
- 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
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
- G01V3/26—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging operating with magnetic or electric fields produced or modified either by the surrounding earth formation or by the detecting device
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- 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
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Abstract
本发明公开一种井下直流赤道向偶极动源异常自显式超前探测方法,该方法将接收电极M或将接收电极M和N固定布置在掘进面上,发射电极A和B布置在巷道底板且发射电极AB连线的垂直平分线为巷道底板中线,发射电极A和B从A1B1处向掘进面方向移动,形成依次位于A1B1,A2B2,...,AiBi...,AnBn位置点的偶极移动源;发射电极A和B在每个位置点发射电流一次,接收电极M或MN对应测量电压一次,进行偶极动源超前探测。本发明提高了井下直流超前探测的准确性,是一种利用稳定电流场对称性自显异常的快速探测方法。
Description
技术领域
本发明属于电和电磁法勘探领域,涉及一种井下直流超前探测技术,具体涉及一种井下直流赤道向偶极动源异常自显式超前探测方法。
背景技术
井下直流电法超前探测,是我国煤矿井下探测掘进工作面前方异常体的主要方法之一。现行的直流超前探测方法,通常采用称为三极法的单极发射—偶极接收装置,将发射电极固定布置在掘进工作面上或附近、接收电极向远离掘进面方向移动。此后,又有了建立在三极法基础上改进的多单极—偶极组合法。由于直流电场近距作用的性质,在发射源激励下异常体的响应,总是通过场从异常体周围向外传递,强度随着离开异常体距离的增大而减弱。将发射电极固定在掘进面附近、接收电极向掘进工作面相反方向移动的方法,获得的异常不仅比较微弱,而且掘进面后方的异常体还会对超前探测形成干扰,造成误判。
发明内容
针对现有技术中的缺陷和不足,本发明提供了一种井下直流赤道向偶极动源异常自显式超前探测方法,以将接收电极固定布置在掘进工作面上尽可能地靠近前方异常体,通过发射电极向掘进工作面方向移动实现超前探测。该方法以电极与巷道间的几何对称关系,形成对称的电流场。在探测过程中,当接收电极测量到超过噪声电平的信号时,便可将掘进面前方的异常体自动地显示出来。
为达到上述目的,本发明采取如下的技术方案:
一种井下直流赤道向偶极动源异常自显式超前探测方法,该方法将接收电极M或将接收电极M和N固定布置在掘进面上且接收电极位于掘进面底边的垂直平分线上,发射电极A和B布置在巷道底板且发射电极A和B连线的垂直平分线为巷道底板中线,发射电极A和B从A1和B1位置处向掘进面方向移动,形成依次位于A1和B1位置点,A2和B2位置点,...,Ai和Bi位置点...,An和Bn位置点的偶极移动源;发射电极A和B在每个位置点发射电流一次,接收电极M或接收电极M和N对应测量电压一次,进行偶极动源超前探测;
当单极接收电极M固定布置在掘进面时,超前探测最大探测距离估算公式为:
当偶极将接收电极M和N固定布置在掘进面时,超前探测最大探测距离估算公式为:
上式中,O′1O是发射电极距AB中点O′1到掘进面底边中点O的距离,O″O是接收电极距MN中点O″到O点的距离,IABmax为最大发射电流,为接收电极距MN观测到的噪声电平,ρ为掘进面前方岩层的电阻率。
本发明还包括如下技术特征:
可选的,所述A1和B1位置点,A2和B2位置点,...,Ai和Bi位置点...,An和Bn位置点之间的间隔越小,探测的分辨率越高;发射电极A和B从A1和B1位置点处开始,发射电极距AB中点O′沿巷道底板中线,以分辨率要求的间隔向掘进面方向移动,电流场透入掘进工作面前方的深度由浅到深。
Di≈O′1O-O′iO,(i=1,2,...,n) (2a)
上式中的O′1O由式(1a)确定,O′iO是发射电极距AB中点O′i到O的距离;
Di≈O′1O-O′iO,(i=1,2,...,n) (2c)
其中的O′1O由式(1b)确定,O′iO是发射电极距AB中点O′i到O的距离;
式中的是Ai到M的距离、/>是Bi到M的距离、/>是Ai到N的距离、/>是Bi到N的距离,其中MO是M到O的距离,NO是N到O的距离,/>是发射电极A和B位于Ai和Bi位置点对应的发射电流,/>是与/>对应的观测电压。/>
可选的,探测结束后,将实测数据中的误差记录作为异常解释的判别标准;
当单极接收电极M固定布置在掘进面时,满足下式(3a)则将实测数据解释为异常:
当接收电极M和N固定布置在掘进面时,满足下式(3b)则将实测数据解释为异常:
式中,Mean±S.D.为平均数±标准差,AiM、BiM、AiN和BiN采用公式(2)中的定义。
可选的,当掘进面前方最大探测距离范围内无异常体时,由于电极布置对称性导致电流场的对称性,接收电极M或接收电极M和N观测到的是噪声电平;当掘进面前方最大探测距离范围内存在异常体时,电流场失去对称性,当接收电极观测到3~5倍噪声电平的电压信号,即可判定掘进面前方存在异常体。
可选的,当掘进面前方最大探测距离范围内存在异常体,但电流场尚未接触到该异常体,接收电极M或接收电极M和N观测到的是噪声电平;随着发射电极A和B向掘进面方向移动,电流场与异常体接触并被扰动,当扰动传递到掘进面上且接收电极观测到3~5倍噪声电平的电压信号,即可判定异常体的存在。
本发明与现有技术相比,有益的技术效果是:
(1)根据直流电场近距作用的原理,将接收电极固定布置在掘进工作面上,最大限度地接近掘进面前方异常体在发射源激励下的响应信号;发射电极向掘进面方向移动、电流场透入掘进面前方的深度由浅到深的观测方式,信噪比随探测距离的增大而增大。
(2)电极布置的几何对称关系,使得只有当掘进面前方存在异常体时,接收电极才会观测到超出噪声的信号。根据信号大于3~5倍噪声电平的标准,可以在施工现场快速判断出前方异常,提高了超前探测的实时性。
(3)给出的最大探测距离估算公式,建立了电极布置的几何空间、发射电流、岩层电阻率以及环境噪声电平等各项要素之间的关系,为井下直流赤道向偶极动源异常自显式超前探测施工参数的确定提供了定量依据。
(4)给出的与发射电极每移动一次所对应的测深点和视电阻率计算公式,是针对本发明设计的算法,具有专门性;根据实测数据误差记录作为异常响应的判别标准,为数据解释的可靠性提供了更多的保证。
(5)本发明的赤道向偶极装置发射,可在极短的巷道内进行施工,所需发射设备轻便,并且没有轴向偶极发射那样的自身探测距离损失;在接收装置方面,有单极和偶极两种接收方式,分别固定布置在掘进工作面的垂直平分线上,不仅能够根据电流场的对称性变化在施工现场快速判断出掘进面前方是否存在地质异常体,及时为井下掘进施工提供可靠资料;而且能适应不同的噪声电平,可以根据不同工况选用。
附图说明
图1为本发明施工布置示意图,其中1a是单极接收,1b是偶极接收。
图2为本发明超前探测测深点示意图,其中2a是单极接收,2b是偶极接收。
图3为实验结果图。
图中:1–掘进面;2–巷道;3–异常体;4–大地。
以下结合说明书附图和具体实施方式对本发明做具体说明。
具体实施方式
为了进一步提高井下直流超前探测的准确性,本发明的赤道向偶极动源异常自显式超前探测方法,接收电极M或接收电极M和N固定布置在掘进面的垂直平分线上,发射电极A和B布置在巷道底板中线的垂直平分线上,从距离掘进面A1B1处沿巷道底板中线、以分辨率要求的间隔,向掘进面方向移动形成一系列的动源AiBi。当掘进面前方超前探测距离内不存在异常体时,电极M或MN接收到的是噪声电平或/>当掘进面前方有异常体存在时,电流场会失去原有的对称性,异常体响应信号将被加载到接收电极M或接收电极M和N上,由此可在施工现场快速判断出掘进面前方的异常体。本发明给出了赤道向偶极动源两种接收方式的最大探测距离估算公式,与AiBi对应的各测深点以及视电阻率计算公式,最大探测距离估算公式中噪声电平的获取方法,探测过程中信噪比的测定方法,利用实测数据中误差记录判别异常的公式等。
遵从上述技术方案,以下给出本发明的具体实施例,需要说明的是本发明并不局限于以下具体实施例,凡在本申请技术方案基础上做的等同变换均落入本发明的保护范围。下面结合实施例对本发明做进一步详细说明。
实施例1:
本实施例提供一种井下直流赤道向偶极动源异常自显式超前探测方法,分别如图1a和图1b所示,该方法将接收电极M或接收电极M和N固定布置在掘进面上且接收电极位于掘进面底边的垂直平分线上,尽可能地接近掘进工作面前方地质异常体、避免巷道顶底板、侧帮附近异常体的影响;发射电极A和B布置在巷道底板且发射电极A和B连线的垂直平分线为巷道底板中线,发射电极A和B从A1B1处沿巷道底板中线以分辨率要求的间隔向掘进面方向移动,信噪比随探测距离的增加而增加,形成依次位于A1和B1位置点,A2和B2位置点,...,Ai和Bi位置点...,An和Bn位置点的偶极移动源;发射电极A和B在每个位置点发射电流一次,接收电极M或接收电极M和N对应测量电压一次,进行偶极动源超前探测;当忽略巷道的影响时,发射电极初始位置与超前探测最大探测距离位置相对掘进面对称,从而发射电极初始位置与掘进面底边中点O间的距离为最大探测距离;
当单极接收电极M固定布置在掘进面时,超前探测最大探测距离估算公式为:
当偶极接收电极MN固定布置在掘进面时,超前探测最大探测距离估算公式为:
上式中,O′1O是发射电极距AB中点O′1到掘进面底边中点O的距离,O″O是接收电极MN中点O″到O点的距离,IABmax为最大发射电流,为将接收电极M和N观测到的噪声电平,ρ为掘进面前方岩层的电阻率。
A1和B1位置点,A2和B2位置点,...,Ai和Bi位置点...,An和Bn位置点之间的间隔越小,探测的分辨率越高;发射电极A和B从A1和B1位置点处开始,发射电极距AB中点O′沿巷道底板中线,以分辨率要求的间隔向掘进面方向移动,电流场透入掘进工作面前方的深度由浅到深。
在上述探测过程中,对于每一个发射电极位于AiBi位置,通过发射机未发送和发送电流时接收机的空采和实采,测定噪声电平或/>和接收电压/>或/>根据信噪比决定是否进行重复观测、对掘进面前方是否存在异常体在现场做出初步判断。当井下作业时间允许时,不要省略此步骤。
由公式(2)计算发射电极位于Ai所对应的测深点和视电阻率,判断掘进面前方异常体是低阻体还是高阻体。
Di≈O′1O-O′iO,(i=1,2,...,n) (2a)
上式中的O′1O由式(1a)确定,O′iO是发射电极距AB中点O′i到O的距离;
Di≈O′1O-O′iO,(i=1,2,...,n) (2c)
其中的O′1O由式(1b)确定,O′iO是发射电极距AB中点O′i到O的距离;
探测结束后,将实测数据中的误差记录作为异常解释的判别标准;
当单极接收电极M固定布置在掘进面时,满足下式(3a)则将实测数据解释为异常:
当偶极接收电极M和N固定布置在掘进面时,满足下式(3b)则将实测数据解释为异常:
式中,Mean±S.D.为平均数±标准差,AiM、BiM、AiN和BiN采用公式(2)中的定义。
当掘进面前方最大探测距离范围内无异常体时,由于电极布置对称性导致电流场的对称性,接收电极M或接收电极M和N观测到的是噪声电平;当掘进面前方最大探测距离范围内存在异常体时,电流场失去对称性,当接收电极观测到3~5倍噪声电平的电压信号,即可判定掘进面前方存在异常体。
当掘进面前方最大探测距离范围内存在异常体,但电流场尚未接触到该异常体,接收电极M或接收电极M和N观测到的是噪声电平;随着发射电极AiBi向掘进面方向移动,电流场与异常体接触并被扰动,当扰动传递到掘进面上且接收电极观测到3~5倍噪声电平的电压信号,即可判定异常体的存在。
实验验证:
为了验证公式的正确性、有效性,利用公式进行了数值模拟计算。假定在巷道底板下20m处有1个半径5m的高阻体,其电阻率为1000Ω.m,围岩的电阻率为100Ω.m,采用赤道向偶极发射、偶极MN接收装置,利用公式(2d)进行计算。图3是相应的计算结果,图3中横坐标为距离,单位为m,纵坐标为视电阻率,单位为Ω.m。由该图可见,经公式计算出来的异常最大值的位置与实际模型位置对应,说明公式是有效的、可用的。
Claims (4)
1.一种井下直流赤道向偶极动源异常自显式超前探测方法,其特征在于,该方法将接收电极M或将接收电极M和N固定布置在掘进面上且接收电极位于掘进面底边的垂直平分线上,发射电极A和B布置在巷道底板且发射电极A和B连线的垂直平分线为巷道底板中线,发射电极A和B从A1和B1位置处向掘进面方向移动,形成依次位于A1和B1位置点,A2和B2位置点,...,Ai和Bi位置点...,An和Bn位置点的偶极移动源;发射电极A和B在每个位置点发射电流一次,接收电极M或接收电极M和N对应测量电压一次,进行偶极动源超前探测;
当单极接收电极M固定布置在掘进面时,超前探测最大探测距离估算公式为:
当偶极将接收电极M和N固定布置在掘进面时,超前探测最大探测距离估算公式为:
上式中,O1′O是发射电极距AB中点O1′到掘进面底边中点O的距离,O″O是接收电极距MN中点O″到O点的距离,IABmax为最大发射电流,为接收电极距MN观测到的噪声电平,ρ为掘进面前方岩层的电阻率;
Di≈O′1O-O′iO,(i=1,2,...,n) (2a)
上式中的O′1O由式(1a)确定,O′iO是发射电极距AB中点O′i到O的距离;
Di≈O′1O-O′iO,(i=1,2,…,n) (2c)
其中的O′1O由式(1b)确定,O′iO是发射电极距AB中点O′i到O的距离;
式中的是Ai到M的距离、/>是Bi到M的距离、/>是Ai到N的距离、/>是Bi到N的距离,其中MO是M到O的距离,NO是N到O的距离,/>是发射电极A和B位于Ai和Bi位置点对应的发射电流,/>是与/>对应的观测电压;
探测结束后,将实测数据中的误差记录作为异常解释的判别标准;
当单极接收电极M固定布置在掘进面时,满足下式(3a)则将实测数据解释为异常:
当接收电极M和N固定布置在掘进面时,满足下式(3b)则将实测数据解释为异常:
式中,Mean±S.D.为平均数±标准差,AiM、BiM、AiN和BiN采用公式(2)中的定义。
2.如权利要求1所述的井下直流赤道向偶极动源异常自显式超前探测方法,其特征在于,所述A1和B1位置点,A2和B2位置点,…,Ai和Bi位置点…,An和Bn位置点之间的间隔越小,探测的分辨率越高;发射电极A和B从A1和B1位置点处开始,发射电极距AB中点O′沿巷道底板中线,以分辨率要求的间隔向掘进面方向移动,电流场透入掘进工作面前方的深度由浅到深。
3.如权利要求1所述的井下直流赤道向偶极动源异常自显式超前探测方法,其特征在于,当掘进面前方最大探测距离范围内无异常体时,由于电极布置对称性导致电流场的对称性,接收电极M或接收电极M和N观测到的是噪声电平;当掘进面前方最大探测距离范围内存在异常体时,电流场失去对称性,当接收电极观测到3~5倍噪声电平的电压信号,即可判定掘进面前方存在异常体。
4.如权利要求1所述的井下直流赤道向偶极动源异常自显式超前探测方法,其特征在于,当掘进面前方最大探测距离范围内存在异常体,但电流场尚未接触到该异常体,接收电极M或接收电极M和N观测到的是噪声电平;随着发射电极A和B向掘进面方向移动,电流场与异常体接触并被扰动,当扰动传递到掘进面上且接收电极观测到3~5倍噪声电平的电压信号,即可判定异常体的存在。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110896207.1A CN113885086B (zh) | 2021-08-05 | 2021-08-05 | 一种井下直流赤道向偶极动源异常自显式超前探测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110896207.1A CN113885086B (zh) | 2021-08-05 | 2021-08-05 | 一种井下直流赤道向偶极动源异常自显式超前探测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113885086A CN113885086A (zh) | 2022-01-04 |
CN113885086B true CN113885086B (zh) | 2023-06-13 |
Family
ID=79010921
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110896207.1A Active CN113885086B (zh) | 2021-08-05 | 2021-08-05 | 一种井下直流赤道向偶极动源异常自显式超前探测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113885086B (zh) |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2020100984A4 (en) * | 2020-06-11 | 2020-08-13 | Xi’an Northwest Nonferrous Geophysical & Geochemical Exploration Co., Ltd | A Method and Apparatus for Ground-tunnel Wide Field Electromagnetic Surveying |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE19842975B4 (de) * | 1998-09-19 | 2004-01-29 | Kaus, Arnim, Dr. | Geoelektrisches Vorauserkundungsverfahren für maschinelle und konventionelle Tunnel- und Hohlraumvortriebe und Messvorrichtung dazu |
US6496008B1 (en) * | 2000-08-17 | 2002-12-17 | Digital Control Incorporated | Flux plane locating in an underground drilling system |
CN101603423B (zh) * | 2009-07-09 | 2013-06-05 | 煤炭科学研究总院西安研究院 | 一种在煤矿巷道内顺层超前探测含水构造的直流电法方法 |
CN101603419A (zh) * | 2009-07-09 | 2009-12-16 | 煤炭科学研究总院西安研究院 | 一种采煤工作面顺煤层透视的矿井直流电法探测方法 |
CN103278855B (zh) * | 2013-05-13 | 2015-06-24 | 江苏大学 | 一种消除巷道和地形对直流勘探视电阻率影响的方法 |
CN103278857B (zh) * | 2013-05-13 | 2015-11-04 | 江苏大学 | 一种井下直流超前探测的电极布置方法 |
CN106772644B (zh) * | 2016-12-14 | 2018-06-26 | 中国矿业大学 | 矿井瞬变电磁三分量探测方法 |
CN111708088B (zh) * | 2020-06-28 | 2022-04-15 | 中国矿业大学 | 一种基于磁梯度张量的瞬变电磁实时动态超前探测方法及系统 |
-
2021
- 2021-08-05 CN CN202110896207.1A patent/CN113885086B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2020100984A4 (en) * | 2020-06-11 | 2020-08-13 | Xi’an Northwest Nonferrous Geophysical & Geochemical Exploration Co., Ltd | A Method and Apparatus for Ground-tunnel Wide Field Electromagnetic Surveying |
Also Published As
Publication number | Publication date |
---|---|
CN113885086A (zh) | 2022-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772644B (zh) | 矿井瞬变电磁三分量探测方法 | |
CN103995296B (zh) | 瞬变电磁法地孔探测方法与装置 | |
CN109613616B (zh) | 一种浅层岩溶隧道地质综合立体探测方法 | |
WO2020078003A1 (zh) | 一种时间域瞬变电磁波测井边界远探测方法 | |
CN108957563B (zh) | 一种隧道施工超前地质探测系统以及探测方法 | |
CN109143378B (zh) | 一种在煤矿巷道内顺层超前探测含水构造的二次时差方法 | |
US20120130641A1 (en) | Marine Source To Borehole Electromagnetic Mapping Of Sub-Bottom Electrical Resistivity | |
CN109343130B (zh) | 一种侧向激发的回线源地井瞬变电磁探测方法及系统 | |
CN107861159A (zh) | 双电偶源地‑井瞬变电磁探测方法 | |
EP3749988B1 (en) | Mapping fracture length using downhole ground penetrating radar | |
CN103174413A (zh) | 一种钻具以及井下随钻探测储层界面和厚度的方法 | |
CN108828678B (zh) | 一种隧道施工超前地质探测系统 | |
CN108828676A (zh) | 一种地面-矿井巷道瞬变电磁三分量探测方法 | |
CN107884834A (zh) | 同源多场瞬变电磁探测方法 | |
CN107065023B (zh) | 一种基于井上下空间立体瞬变电磁数据采集的探测系统 | |
CN113885086B (zh) | 一种井下直流赤道向偶极动源异常自显式超前探测方法 | |
CN113885083B (zh) | 一种井下直流轴向偶极动源异常自显式超前探测方法 | |
CN113900151B (zh) | 一种井下直流单极动源异常自显式超前探测方法 | |
CN109738964B (zh) | 地震波和电磁波联合反演的隧道预报装置、掘进机及方法 | |
CN113885085B (zh) | 一种井下直流轴向偶极动源超前探测方法 | |
CN113885084B (zh) | 一种井下直流单极动源超前探测方法 | |
CN104991281B (zh) | 一种煤层埋藏深度的探测方法和装置 | |
US20210072420A1 (en) | Low frequency complex resistivity measurement in a formation | |
CN204422776U (zh) | 一种识别并提取煤矿巷道侧帮近似平行小断层信号的装置 | |
CN110805433B (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 |