CN105045962B - 边采边复耕地区动态施工标高确定方法 - Google Patents
边采边复耕地区动态施工标高确定方法 Download PDFInfo
- Publication number
- CN105045962B CN105045962B CN201510357534.4A CN201510357534A CN105045962B CN 105045962 B CN105045962 B CN 105045962B CN 201510357534 A CN201510357534 A CN 201510357534A CN 105045962 B CN105045962 B CN 105045962B
- Authority
- CN
- China
- Prior art keywords
- stage
- mrow
- absolute altitude
- area
- earth
- 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
Landscapes
- Position Fixing By Use Of Radio Waves (AREA)
- Steroid Compounds (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及边采边复耕地区动态施工标高确定方法，属于采矿技术、土地利用和土地复垦技术领域。该方法包括：利用概率积分法分阶段动态沉陷预计，同时考虑地表原始高程，模拟各阶段地表动态沉陷情况，然后根据当地的自然社会经济情况，确定边采边复耕地区复垦标高，结合边采边复空间布局中耕地区的范围以及当地潜水位埋深，确定耕地动态施工位置范围，最后根据复垦标高和地表后续下沉，确定边采边复耕地区动态施工标高。本发明综合考虑了边采边复工程施工时，需复垦耕地区的范围、地表沉陷情况、后续下沉以及复垦标高，从而确定了边采边复耕地区动态施工标高，简单易懂，便于指导工程施工，为边采边复技术的实地实施提供了有力保障。
Description
技术领域
本发明属于采矿技术、土地利用和土地复垦技术领域，特别涉及中高潜水位平原矿区边采边复工程施工过程中耕地区动态施工标高的确定方法。
背景技术
据不完全统计，新中国成立以来，我国已开采出煤炭约445亿吨，对社会经济发展做出了巨大的贡献。然而煤炭的开采也带来了诸多问题，我国96％的煤炭产量来自于地下开采，并多采用走向长壁全部垮落法管理顶板，土地沉陷系数大，采煤沉陷土地占煤炭开采破坏土地的91％。据保守估计，我国由于煤炭资源开采累计沉陷面积超过100万公顷，且以每年3.0万～4.7万公顷的速度增加。因此，采煤沉陷土地复垦已成为中国最具特色的难题。
自20世纪80年代，中国的土地复垦在利用不同方法修复采煤沉陷土地方面，取得了显著的成效。现在通常使用的方法有四种：直接修整法(土地平整)、疏排法、挖深垫浅法和利用煤矸石、粉煤灰、湖泥等充填材料进行充填复垦。但这些方法与环境领域的“先污染、后治理”的情况一样，大都是“先破坏、后复垦”，致力于治理稳沉的土地，也就是说煤炭开采结束，土地沉陷达到最终状态后，再采取措施进行复垦。这时环境受到严重破坏，对土地的大部分损伤也已经形成，尤其是在中高潜水位平原矿区，土地稳沉后，将有大量肥沃的土地沉入水中而无法挽救。因此，为了最大限度地保护土地资源，缩短复垦时间，实现及时复垦，井工矿山边开采边复垦技术——边采边复技术应运而生。
目前对于边采边复技术的复垦时机、复垦空间布局等先后有了不断深入的探索和研究，并取得了丰硕的成果。中国矿业大学(北京)已申请的专利——用于采煤沉陷地边采边复的基于土方平衡的基塘布局方法，在土方平衡的基础上确定了最佳的耕地区与挖深区范围，为复垦空间布局的确定提供了一定的科学与理论依据。但该方法仅确定了最终的复垦空间布局，并不能很好地指导整个动态施工过程，因此在复垦空间布局(即最终耕地区和挖深区范围)确定的前提下，如何确定某一时刻需要复垦的耕地范围，以及耕地的施工标高，从而更好的指导工程施工，对于边采边复技术思想的实地实施至关重要。
发明内容
本发明的目的是为解决上述难题，提出边采边复耕地区动态施工标高的确定方法，是在采用边采边复技术确定了复垦布局(耕地区与水域)与最终复垦标高的基础上，通过分析复垦时间节点地面沉陷情况、后续下沉、最终复垦标高，确定各复垦时间节点的动态复垦施工标高，从而科学指导边采边复技术的实施。
本发明提出的边采边复耕地区动态施工标高的确定方法，该方法主要适用于中高潜水位平原矿区，包括利用概率积分法分阶段动态沉陷预计，同时考虑地表原始高程，模拟各阶段地表动态沉陷情况，然后确定边采边复耕地区复垦标高，结合边采边复空间布局中耕地区的范围以及当地潜水位埋深，确定耕地动态施工位置范围，最后根据复垦标高和地表后续下沉，确定边采边复耕地区动态施工标高。
该方法具体包括以下步骤：
1)分阶段地表采煤沉陷动态模拟：根据地下煤炭赋存条件、开采计划，以月或年为单位将开采划分为多个阶段，利用概率积分法分阶段进行动态沉陷预计，并与煤炭开采前地表的原始高程进行叠加分析，模拟各阶段地表采煤沉陷情况。具体如下：设地表任意点P的坐标为(x,y)，则各阶段i(i＝1,2,3,……n)点P的高程值，如表达式(1)所示：
式(1)中：Hi(x,y)为阶段i地表任意点P的高程，单位m；
H0(x,y)为地表任意点P的原始高程，单位m；
Wi(x,y)为阶段i地表任意点P的下沉值，单位mm。
2)确定边采边复耕地区复垦标高：边采边复耕地区的复垦标高由农作物生长的临界深度、地下潜水位、外河水位及5年一遇洪水位决定，如表达式(2)所示：
HR＝hL+max{HD,HW,HH} (2)
式(2)中：HR为充填耕地区复垦标高，单位m；hL为农作物生长临界深度，单位m；
max{}为取括号中的最大者；HD为地下潜水位高程，单位m；
HW为外河水位高程，单位m；HH为5年一遇洪水位高程，单位m。
3)确定各阶段需复垦耕地区的位置范围：根据边采边复空间布局中的耕地区范围DG，步骤1)中各阶段地表采煤沉陷情况，以及当地潜水位埋深，确定边采边复耕地区范围内各阶段的积水区域Di，则各阶段需复垦耕地区的范围Fi如表达式(3)所示：
4)确定各阶段需复垦耕地区的施工标高：根据步骤2)中确定的边采边复耕地区复垦标高，同时考虑地表的后续下沉，确定各阶段需复垦耕地区的施工标高，如表达式(4)所示：
HSi(xQ,yQ)＝HR+(Hi(xQ,yQ)-Hn(xQ,yQ)) (4)
式(4)中：HSi(xQ,yQ)为阶段i需复垦耕地区任意点Q的施工标高，单位m；
Hi(xQ,yQ)为阶段i需复垦耕地区任意点Q的高程，单位m；
Hn(xQ,yQ)为阶段n需复垦耕地区任意点Q的高程，单位m。
本发明主要具有以下技术优点：
本发明通过分阶段的沉陷预计和地表原始高程，清晰的模拟地表动态采煤沉陷情况，在此基础上，综合考虑当地的自然社会经济情况，确定边采边复耕地区复垦标高，结合边采边复空间布局中的耕地区范围、当地潜水位埋深和耕地区复垦标高，综合确定边采边复耕地区动态施工范围及其施工标高。
本发明用于指导边采边复工程施工，方法简单易懂，易于实地操作，是边采边复技术思想实地实施的有力保障。
附图说明
图1为边采边复耕地区动态施工标高确定方法流程图。
图2为点P、Q位置及各阶段需复垦耕地区位置范围图。
图3为阶段2需复垦耕地区施工标高剖面图。
具体实施方式：
本发明提出的边采边复耕地区动态施工标高的确定方法，结合附图及实施例详细说明如下：
本发明的方法主要适用于中高潜水位平原矿区，包括利用概率积分法分阶段动态沉陷预计，同时考虑地表原始高程，模拟各阶段地表动态沉陷情况，然后根据当地的自然社会经济情况，确定边采边复耕地区复垦标高，结合边采边复空间布局中耕地区的范围以及当地潜水位埋深，确定耕地动态施工位置范围，最后根据复垦标高和地表后续下沉，确定边采边复耕地区动态施工标高。
该方法具体包括以下步骤：
1)分阶段地表采煤沉陷动态模拟：根据地下煤炭赋存条件、开采计划，以月或年为单位将开采划分为多个阶段，利用概率积分法分阶段进行动态沉陷预计，并与煤炭开采前地表的原始高程进行叠加分析，模拟各阶段地表采煤沉陷情况。具体如下：设地表任意点P的坐标为(x,y)，则各阶段i(i＝1,2,3,……n)点P的高程值，如表达式(1)所示：
式(1)中：Hi(x,y)为阶段i地表任意点P的高程，单位m；
H0(x,y)为地表任意点P的原始高程，单位m；
Wi(x,y)为阶段i地表任意点P的下沉值，单位mm。
2)确定边采边复耕地区复垦标高：边采边复耕地区的复垦标高由农作物生长的临界深度、地下潜水位、外河水位及5年一遇洪水位决定，如表达式(2)所示：
HR＝hL+max{HD,HW,HH} (2)
式(2)中：HR为充填耕地区复垦标高，单位m；hL为农作物生长临界深度，单位m；
max{}为取括号中的最大者；HD为地下潜水位高程，单位m；
HW为外河水位高程，单位m；HH为5年一遇洪水位高程，单位m。
3)确定各阶段需复垦耕地区的位置范围：根据边采边复空间布局中的耕地区范围DG，步骤1)中各阶段地表采煤沉陷情况，以及当地潜水位埋深，确定边采边复耕地区范围内各阶段的积水区域Di，则各阶段需复垦耕地区的范围Fi如表达式(3)所示：
4)确定各阶段需复垦耕地区的施工标高：根据步骤2)中确定的边采边复耕地区复垦标高，同时考虑地表的后续下沉，确定各阶段需复垦耕地区的施工标高，如表达式(4)所示：
HSi(xQ,yQ)＝HR+(Hi(xQ,yQ)-Hn(xQ,yQ)) (4)
式(4)中：HSi(xQ,yQ)为阶段i需复垦耕地区任意点Q的施工标高，单位m；
Hi(xQ,yQ)为阶段i需复垦耕地区任意点Q的高程，单位m；
Hn(xQ,yQ)为阶段n需复垦耕地区任意点Q的高程，单位m。
实施例：
本实施例为东部某一高潜水位平原矿区，区内地势平坦，地表自然高程在+43.2m～+44.6m之间，平均+43.9m，地形坡度为2‰，地下潜水位埋深为1.5m左右，煤层平均厚度5.0m，埋藏深度800m，工作面长250m，工作面推进长度1050m，工作面推进速度5m/天。
本实施例基于采煤沉陷地边采边复技术思想，在计算机软件环境下模拟确定充填区施工标高，其方法流程如图1所示，包括以下步骤：
1)分阶段地表采煤沉陷动态模拟：根据地下煤炭赋存条件、开采计划，以月或年为单位将开采划分为多个阶段，本实施例中，工作面开采共210天，因此以月为单位划分为7个阶段，每个阶段开采150m。利用概率积分法分阶段进行动态沉陷预计，并与煤炭开采前地表的原始高程进行叠加分析，模拟各阶段地表采煤沉陷情况。本实施例中，利用概率积分法对7个阶段分别进行沉陷预计，地表任意点P的坐标为(810,-237)，位置如图2中所示，Fi表示阶段i需复垦地区的位置范围，i＝2,3,4,5,6,7点P在阶段4预计的下沉值为W4(810,-237)＝500mm，点P的原始高程为H0(810,-237)＝43.7m，因此点P在阶段4煤炭开采后的高程H4(810,-237)＝43.7-500/1000m＝43.2m。。
2)确定边采边复耕地区复垦标高：边采边复耕地区的复垦标高由农作物生长的临界深度、地下潜水位、外河水位及5年一遇洪水位决定，本实施例中，农作物生长临界深度hL＝1.0m，地下潜水位高程HD＝42.4m，外河水位高程HW＝41.8m，5年一遇洪水位高程HH＝42.1m，因此充填耕地区复垦标高HR＝1.0+max{42.4,41.8,42.1}m＝43.4m。
3)确定各阶段需复垦耕地区的位置范围：根据边采边复空间布局中的耕地区范围DG，步骤1)中各阶段地表采煤沉陷情况，以及当地潜水位埋深，确定边采边复耕地区范围内各阶段的积水区域Di，进而确定各阶段需复垦耕地区的范围Fi。本实施例中，地下潜水位高程HD＝42.4m，结合步骤1)中对7个阶段地表采煤沉陷情况的模拟，得出阶段3地表开始出现积水，因此F1＝D2-D1＝0，其他各阶段需复垦耕地区的范围如图2所示。
4)确定各阶段需复垦耕地区的施工标高：根据步骤2)中确定的边采边复耕地区复垦标高，同时考虑地表的后续下沉，确定各阶段需复垦耕地区的施工标高。本实施例中，阶段2需复垦耕地区任意点Q原始高程H0(223，-125)＝44.1m，阶段2预计的下沉值为W2(223，-125)＝700mm，阶段7预计的下沉值为W7(223，-125)＝1800mm，因此H2(223，-125)＝43.4m，H7(223，-125)＝42.3m。则阶段2需复垦耕地区任意点Q的施工标高HS2(223，-125)＝43.4+(43.4-42.3)m＝44.5m。边采边复阶段2需复垦耕地区施工标高剖面如图3所示，图中A-B表示原始地面，C表示当地潜水位，D表示可挖掘取土线，E表示煤层，曲线1表示阶段2下沉盆地范围，曲线2表示最终下沉盆地范围，带点的区域表示阶段2需复垦地区。
Claims (1)
1.一种边采边复耕地区动态施工标高确定方法，其特征在于，该方法包括：利用概率积分法分阶段动态沉陷预计，同时考虑地表原始高程，模拟各阶段地表动态沉陷情况，然后根据当地的自然社会经济情况，确定边采边复耕地区复垦标高，结合边采边复空间布局中耕地区的范围以及当地潜水位埋深，确定耕地动态施工位置范围，最后根据复垦标高和地表后续下沉，确定边采边复耕地区动态施工标高；
具体包括以下步骤：
1)分阶段地表采煤沉陷动态模拟：根据地下煤炭赋存条件、开采计划，以月或年为单位将开采划分为多个阶段，利用概率积分法分阶段进行动态沉陷预计，并与煤炭开采前地表的原始高程进行叠加分析，模拟各阶段地表采煤沉陷情况；具体如下：设地表任意点P的坐标为(x,y)，则各阶段i点P的高程值，如表达式(1)所示：
<mrow>
<msub>
<mi>H</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>H</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>.</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mfrac>
<mrow>
<msub>
<mi>W</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>x</mi>
<mo>,</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
</mrow>
<mn>1000</mn>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式(1)中：Hi(x,y)为阶段i地表任意点P的高程，单位m；
H0(x,y)为地表任意点P的原始高程，单位m；
Wi(x,y)为阶段i地表任意点P的下沉值，单位mm；
i＝1,2,3,……n；
2)确定边采边复耕地区复垦标高：边采边复耕地区的复垦标高由农作物生长的临界深度、地下潜水位、外河水位及5年一遇洪水位决定，如表达式(2)所示：
HR＝hL+max{HD,HW,HH} (2)
式(2)中：HR为耕地区复垦标高，单位m；hL为农作物生长临界深度，单位m；
max{}为取括号中的最大者；HD为地下潜水位高程，单位m；
HW为外河水位高程，单位m；HH为5年一遇洪水位高程，单位m；
3)确定各阶段需复垦耕地区的位置范围：根据边采边复空间布局中的耕地区范围DG，步骤1)中各阶段地表采煤沉陷情况，以及当地潜水位埋深，确定边采边复耕地区范围内各阶段的积水区域Di，则各阶段需复垦耕地区的范围Fi如表达式(3)所示：
<mrow>
<msub>
<mi>F</mi>
<mi>i</mi>
</msub>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>i</mi>
<mo>+</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>i</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>,</mo>
<mn>3</mn>
<mo>,</mo>
<mo>...</mo>
<mi>n</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mi>G</mi>
</msub>
<mo>-</mo>
<msub>
<mi>D</mi>
<mi>n</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mi>n</mi>
</mrow>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
4)确定各阶段需复垦耕地区的施工标高：根据步骤2)中确定的边采边复耕地区复垦标高，同时考虑地表的后续下沉，确定各阶段需复垦耕地区的施工标高，如表达式(4)所示：
HSi(xQ,yQ)＝HR+(Hi(xQ,yQ)-Hn(xQ,yQ)) (4)
式(4)中：HSi(xQ,yQ)为阶段i需复垦耕地区任意点Q的施工标高，单位m；
Hi(xQ,yQ)为阶段i需复垦耕地区任意点Q的高程，单位m；
Hn(xQ,yQ)为阶段n需复垦耕地区任意点Q的高程，单位m。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510357534.4A CN105045962B (zh) | 2015-06-25 | 2015-06-25 | 边采边复耕地区动态施工标高确定方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510357534.4A CN105045962B (zh) | 2015-06-25 | 2015-06-25 | 边采边复耕地区动态施工标高确定方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105045962A CN105045962A (zh) | 2015-11-11 |
| CN105045962B true CN105045962B (zh) | 2018-04-06 |
Family
ID=54452505
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510357534.4A Active CN105045962B (zh) | 2015-06-25 | 2015-06-25 | 边采边复耕地区动态施工标高确定方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105045962B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108663014B (zh) * | 2018-03-30 | 2019-07-30 | 中国水利水电科学研究院 | 一种垄沟耕作坡耕地田块相对高程数据计算方法 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102889083A (zh) * | 2012-09-18 | 2013-01-23 | 中国矿业大学(北京) | 基于gis栅格单元的采煤沉陷地表土剥离时空确定方法 |
| CN103235853A (zh) * | 2013-04-21 | 2013-08-07 | 中国矿业大学(北京) | 用于采煤沉陷地边采边复的基于土方平衡的基塘布局方法 |
-
2015
- 2015-06-25 CN CN201510357534.4A patent/CN105045962B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102889083A (zh) * | 2012-09-18 | 2013-01-23 | 中国矿业大学(北京) | 基于gis栅格单元的采煤沉陷地表土剥离时空确定方法 |
| CN103235853A (zh) * | 2013-04-21 | 2013-08-07 | 中国矿业大学(北京) | 用于采煤沉陷地边采边复的基于土方平衡的基塘布局方法 |
Non-Patent Citations (2)
| Title |
|---|
| 《Development and management of land reclamation in China》;Wei Wang等;《Ocean & Coastal Management》;20141231;第102卷;第415-425页 * |
| 《试论井工煤矿边开采边复垦技术》;胡振琪等;《煤炭学报》;20130228;第38卷(第2期);第301-307页的2.1，2.3，2.3.1，2.3.3，3.1 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105045962A (zh) | 2015-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Hu et al. | Optimization of concurrent mining and reclamation plans for single coal seam: a case study in northern Anhui, China | |
| CN101422094B (zh) | 平原矿区厚煤层开采塌陷地动态预复垦方法 | |
| WO2016045440A1 (zh) | 一种采矿方法 | |
| CN103643956B (zh) | 一种采煤塌陷区煤矸石充填复垦表层覆土厚度设计方法 | |
| CN107882561A (zh) | 一种高潜水位采煤沉陷区超前治理方法 | |
| CN107989612B (zh) | 一种高潜水位煤矿区动态预复垦新方法 | |
| CN103235853B (zh) | 用于采煤沉陷地边采边复的基于土方平衡的基塘布局方法 | |
| CN106437719A (zh) | 煤矿沉陷盆地的修复方法 | |
| Careddu et al. | Promoting ecological sustainable planning for natural stone quarrying. The case of the Orosei Marble Producing Area in Eastern Sardinia | |
| Li et al. | Innovation for sustainable mining: Integrated planning of underground coal mining and mine reclamation | |
| CN109372511A (zh) | 一种利用黄河泥沙动态充填采煤塌陷地的方法 | |
| CN104730224A (zh) | 沟谷地形下煤炭开采对地表径流影响的实验装置及方法 | |
| CN103541376B (zh) | 采煤沉陷区地基在重复开采条件下的基础变形预测方法 | |
| Elsheikh | Mitigation of groundwater level deterioration of the Nubian Sandstone aquifer in Farafra Oasis, Western Desert, Egypt | |
| CN113077088A (zh) | 城中采煤沉陷区水土资源空间格局重构方法 | |
| CN102877846B (zh) | 基于单元法的边采边复表土剥离时机确定方法 | |
| CN104564069A (zh) | 一种基于方格网法的地面动态沉陷预测与复垦方法 | |
| CN103742151B (zh) | 应用于边采边复技术的土壤重构方法 | |
| CN105045962B (zh) | 边采边复耕地区动态施工标高确定方法 | |
| CN102918950B (zh) | 一种采煤沉陷地复垦施工敏感区的确定方法 | |
| CN111400659A (zh) | 采-选-充一体化开采技术耕地保护率计算方法 | |
| CN105808935B (zh) | 多煤层开采下地面积水时间的确定方法 | |
| CN103226738B (zh) | 一种基于预期复垦耕地率和复垦成本的复垦时机优选方法 | |
| Hu et al. | Introduction to concurrent mining and reclamation for coal mines in China | |
| Zhang et al. | Study on the geological environmental disturbance and reclamation technology of underground coalmines: a case study, Xinjiang, China |
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 |