CN111753401B - 基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法 - Google Patents

基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法 Download PDF

Info

Publication number
CN111753401B
CN111753401B CN202010468866.0A CN202010468866A CN111753401B CN 111753401 B CN111753401 B CN 111753401B CN 202010468866 A CN202010468866 A CN 202010468866A CN 111753401 B CN111753401 B CN 111753401B
Authority
CN
China
Prior art keywords
slope
blasting
vibration
rock
cantilever beam
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
CN202010468866.0A
Other languages
English (en)
Other versions
CN111753401A (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.)
Wuhan University of Technology WUT
Original Assignee
Wuhan University of Technology WUT
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 Wuhan University of Technology WUT filed Critical Wuhan University of Technology WUT
Priority to CN202010468866.0A priority Critical patent/CN111753401B/zh
Publication of CN111753401A publication Critical patent/CN111753401A/zh
Application granted granted Critical
Publication of CN111753401B publication Critical patent/CN111753401B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01HMEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
    • G01H17/00Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves, not provided for in the preceding groups
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A10/00TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
    • Y02A10/23Dune restoration or creation; Cliff stabilisation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,包括以下步骤:将岩质高边坡简化为悬臂梁力学模型,建立振动速度预测公式进行现场爆破试验,确定系数k、β1和β2,得到岩质高边坡爆破振动预测公式的具体形式。岩质高边坡爆破开挖过程中,岩质高边坡的表层围岩底部受到约束作用,具有典型的悬臂梁特性,因此将岩质高边坡简化为悬臂梁力学模型,据此建立振动速度预测公式,将影响边坡振动速度的边坡底部约束考虑在内,提高了岩质高边坡地形的爆破坡面振动速度预测的精确度,可以广泛应用于岩质高边坡的爆破开挖工程中,用以指导爆破施工,此发明用于爆破振动技术领域。

Description

基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法
技术领域
本发明涉及爆破振动技术领域,特别是涉及一种基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法。
背景技术
目前,在露天矿山开采、大型水电开挖、岩质深基坑开挖等岩质高边坡工程作业中,仍然以岩体的爆破开挖施工方法为主。其中,岩质高边坡爆破开挖诱发边坡的振动速度是边坡围岩损伤控制的关键指标,可指导实际的爆破工程。
现有常用的爆破振动预测公式为萨道夫斯基预测公式,在平整地形条件下预测地面的爆破振动质点速度具有较高的精度,其具体形式为:
式中K为考虑地质影响及地形影响的系数;α为与地质条件有关的爆破地震波衰减系数;Q为最大单响药量;R为测点与爆源的距离。
由于岩质高边坡地形的爆破振动的传播受到了底部约束的影响,传统的萨道夫斯基预测公式对于岩质高边坡地形的爆破振动速度预测的准确性不高,随着大型水电开挖、岩质深基坑开挖等岩质高边坡工程越来越多,传统萨道夫斯基预测公式已经不能满足实际工程对于爆破振速精确性和安全性的要求。
发明内容
本发明的目的在于至少解决现有技术中存在的技术问题之一,提供一种基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,能够提高岩质高边坡地形的爆破坡面振动速度预测的精确度。
根据本发明的实施例,提供一种基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,包括以下步骤:
S1.将岩质高边坡简化为悬臂梁力学模型,建立振动速度预测公式
式中k为考虑地质影响、边坡底部约束及地形影响的系数,β1为与地质条件有关的爆破地震波衰减系数,β2为边坡底部约束影响因子,Q为最大单响药量,R为测点与爆源的距离,h为测点与爆源的高差,H为岩质高边坡总高度;
S2.进行现场爆破试验,确定系数k、β1和β2,得到岩质高边坡爆破振动预测公式的具体形式。
根据本发明实施例所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,所述步骤S2包括以下步骤:
S210.在边坡坡面进行钻孔装药,记录最大单响药量Q,沿边坡坡面布置爆破振动测试仪,记录爆破振动过程不同高程差位置的震动波形,获得爆破振动速度峰值V,同时测量并记录测点与爆源距离R、测点与爆源高差h、岩质高边坡总高度H;
S220.根据所述步骤S210获得的数据开展二元线性回归分析,确定系数k、β1、β2,以得到岩质高边坡爆破振动预测公式具体形式。
根据本发明实施例所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,所述步骤S220包括以下步骤:
S221.根据所述步骤S1的振动速度预测公式和所述步骤S210的数据计算ln V、和/>
S222.对ln V、和/>进行二元线性回归分析,得到ln V关于/>和/>的线性方程;
S223.根据所述线性方程的斜率和截距,求得系数k、β1、β2,将所述k、β1、β2代入振动速度预测公式,得到岩质高边坡爆破坡面的振动速度预测公式的具体形式。
根据本发明实施例所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,所述步骤S1包括以下步骤:
S101.建立悬臂梁的受迫振动方程;
S102.结合量纲分析方法得到振动速度预测公式。
根据本发明实施例所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,所述步骤S101中悬臂梁的受迫振动方程为
式中f(t)为等效爆破荷载,q(H)sin(pt)为等效梁受围岩的作用荷载,只考虑第一阶振型,则得到梁上任意一点的最大速度峰值与整个梁上最大速度峰值的比值为sin(πh/2H),作为表征边坡各点速度与整个边坡最大速度峰值比值特性的无量纲量。
根据本发明实施例所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,所述步骤S102中,具体步骤为:爆破振动受到测点与爆源的高差h、岩质高边坡总高度H、测点爆心距R、地表岩体质点振动位移μ、地表岩体质点振动加速度a、岩体的固有振动频率f、最大单响药量Q、岩体密度ρ、振动波传播速度c、以及爆轰时间t的影响,由量纲分析得到
V=Φ(Q,R,H,h,f,c,a,t,μ,ρ),
式中取质量Q,速度c,长度R为基本量纲,则有
不同无量纲量πi的乘积和乘方仍为无量纲量,得到
岩质高边坡的振动速度表示为
有益效果:岩质高边坡爆破开挖过程中,岩质高边坡的表层围岩底部受到约束作用,具有典型的悬臂梁特性,因此将岩质高边坡简化为悬臂梁力学模型,据此建立振动速度预测公式,并且通过现场爆破试验记录和测量相关参数,代入振动速度预测公式中获得系数k、β1和β2的值,得到岩质高边坡爆破振动预测公式的具体形式,本基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,将影响边坡振动速度的边坡底部约束考虑在内,提高了岩质高边坡地形的爆破坡面振动速度预测的精确度,可以广泛应用于岩质高边坡的爆破开挖工程中,用以指导爆破施工。
附图说明
下面结合附图对本发明作进一步说明:
图1是本发明实施例现场爆破试验的结构示意图;
图2是本发明实施例爆破振动速度V与高程差h的关系图。
具体实施方式
本部分将详细描述本发明的具体实施例,本发明之较佳实施例在附图中示出,附图的作用在于用图形补充说明书文字部分的描述,使人能够直观地、形象地理解本发明的每个技术特征和整体技术方案,但其不能理解为对本发明保护范围的限制。
在本发明的描述中,需要理解的是,涉及到方位描述,例如上、下、前、后、左、右等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。
在本发明的描述中,若干的含义是一个或者多个,多个的含义是两个以上,大于、小于、超过等理解为不包括本数,以上、以下、以内等理解为包括本数。如果有描述到第一、第二只是用于区分技术特征为目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量或者隐含指明所指示的技术特征的先后关系。
本发明的描述中,除非另有明确的限定,设置、安装、连接等词语应做广义理解,所属技术领域技术人员可以结合技术方案的具体内容合理确定上述词语在本发明中的具体含义。
本发明实施例提供一种基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,包括以下步骤:
S1.将岩质高边坡简化为悬臂梁力学模型,建立振动速度预测公式
式中k为考虑地质影响、边坡底部约束及地形影响的系数,β1为与地质条件有关的爆破地震波衰减系数,β2为边坡底部约束影响因子,Q为最大单响药量,R为测点与爆源的距离,h为测点与爆源的高差,H为岩质高边坡总高度。其中,步骤S1包括步骤S101和步骤S102。
S101.通过简化的悬臂梁力学模型建立悬臂梁的受迫振动方程:
式中f(t)为等效爆破荷载,q(H)sin(pt)为等效梁受围岩的作用荷载,通常对于岩质高边坡这类刚度较大的岩土体结构,只考虑第一阶振型,则求解得到梁上任意一点的最大速度峰值与整个梁上最大速度峰值的比值为sin(πh/2H),以此作为表征边坡各点速度与整个边坡最大速度峰值比值特性的无量纲量。
S102.采用量纲分析方法得到振动速度预测公式:
结合上述动力学分析结果可知岩质高边坡的爆破振动响应跟测点与爆源的高差h、岩质高边坡总高度H、以及两者比值的正弦函数相关。
此外,爆破振动还主要受到地形地貌、地质条件、测点爆心距R、地表岩体质点振动位移μ、地表岩体质点振动加速度a、岩体的固有振动频率f、最大单响药量Q、岩体密度ρ、振动波传播速度c、以及爆轰时间t等因素影响,由量纲分析得到:
V=Φ(Q,R,H,h,f,c,a,t,μ,ρ),
式中取质量Q,速度c,长度R为基本量纲(根据π定理取独立变量为(Q,R,c),则有7个π分量,以πi代表无量纲量),则有:
不同无量纲量πi的乘积和乘方仍为无量纲量,得到:
岩质高边坡的振动速度表示为:
S2.进行现场爆破试验,确定未知系数k、β1和β2,得到岩质高边坡爆破振动预测公式的具体形式,具体地,步骤S2包括步骤S210和步骤S220。
S210.参照图1,在边坡坡面进行钻孔装药,记录最大单响药量Q,沿边坡坡面布置爆破振动测试仪,记录爆破振动过程不同高程差位置的震动波形,获得爆破振动速度峰值V,同时测量并记录测点与爆源距离R、测点与爆源高差h、岩质高边坡总高度H,获得至少3组测试数据;
S220.根据步骤S210获得的数据开展二元线性回归分析,确定系数k、β1、β2,以得到岩质高边坡爆破振动预测公式具体形式。其中,步骤S220还包括以下步骤:
S221.根据步骤S1的振动速度预测公式和步骤S210的数据计算ln V、
S222.对ln V、和/>进行二元线性回归分析,得到ln V关于/>和/>的线性方程;
S223.根据线性方程的斜率和截距,采用自然常数e为底的指数函数,求得系数k、β1、β2的值,将k、β1、β2的值代入振动速度预测公式,得到岩质高边坡爆破坡面的振动速度预测公式的具体形式。
本实施例的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,主要针对传统的萨道夫斯基预测公式在岩质高边坡爆破坡面振动速度预测上的局限性,岩质高边坡爆破开挖过程中,岩质高边坡的表层围岩底部受到约束作用,具有典型的悬臂梁特性,因此将岩质高边坡简化为悬臂梁力学模型,据此建立振动速度预测公式,并且通过现场爆破试验记录和测量相关参数,代入振动速度预测公式中获得系数k、β1和β2的值,得到岩质高边坡爆破振动预测公式的具体形式,以此修正传统的萨道夫斯基预测公式。
在本实施例中,k=106.10,β1=1.55,β2=0.18,从而得到岩质高边坡爆破振动预测公式具体形式将现场实测值、萨道夫斯基预测值和使用本实施例的爆破振动速度预测公式得到的预测值(即修正公式预测值)绘制于同一表格和折线图中,参照表一和图2。
表一
高程差/m 实测值m/s 萨道夫斯基公式预测值m/s 修正公式预测值m/s
1 4.13 4.28 3.91
2 4.01 3.96 4.63
3 4.65 3.83 4.71
4 4.95 3.50 4.61
5 4.59 3.23 4.41
6 4.29 2.82 3.97
7 3.56 2.58 3.67
8 3.23 2.37 3.40
由表一和图2可知,使用本发明实施例的公式的预测精度明显高于目前广泛应用的萨道夫斯基公式,能够实现对岩质高边坡爆破坡面振动速度的精确预测。
本实施例的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,将影响边坡振动速度的边坡底部约束考虑在内,提高了岩质高边坡地形的爆破坡面振动速度预测的精确度,可以广泛应用于岩质高边坡的爆破开挖工程中,用以指导爆破施工,例如根据相关规范(爆破安全规程GB6722-2017)中对结构物的安全振动速度的规定确定V,反推爆破开挖区的最大单响药量Q的限制,具体关系式为:
上面结合附图对本发明实施例作了详细说明,但是本发明不限于上述实施例,在所述技术领域普通技术人员所具备的知识范围内,还可以在不脱离本发明宗旨的前提下作出各种变化。

Claims (6)

1.一种基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,其特征在于,包括以下步骤:
S1.将岩质高边坡简化为悬臂梁力学模型,建立振动速度预测公式
式中k为考虑地质影响、边坡底部约束及地形影响的系数,β1为与地质条件有关的爆破地震波衰减系数,β2为边坡底部约束影响因子,Q为最大单响药量,R为测点与爆源的距离,h为测点与爆源的高差,H为岩质高边坡总高度;
S2.进行现场爆破试验,确定系数k、β1和β2,得到岩质高边坡爆破振动预测公式的具体形式。
2.根据权利要求1所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,其特征在于,所述步骤S2包括以下步骤:
S210.在边坡坡面进行钻孔装药,记录最大单响药量Q,沿边坡坡面布置爆破振动测试仪,记录爆破振动过程不同高程差位置的震动波形,获得爆破振动速度峰值V,同时测量并记录测点与爆源距离R、测点与爆源高差h、岩质高边坡总高度H;
S220.根据所述步骤S210获得的数据开展二元线性回归分析,确定系数k、β1、β2,以得到岩质高边坡爆破振动预测公式具体形式。
3.根据权利要求2所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,其特征在于,所述步骤S220包括以下步骤:
S221.根据所述步骤S1的振动速度预测公式和所述步骤S210的数据计算ln V、和/>
S222.对ln V、和/>进行二元线性回归分析,得到ln V关于/>的线性方程;
S223.根据所述线性方程的斜率和截距,求得系数k、β1、β2,将所述k、β1、β2代入振动速度预测公式,得到岩质高边坡爆破坡面的振动速度预测公式的具体形式。
4.根据权利要求1所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,其特征在于,所述步骤S1包括以下步骤:
S101.建立悬臂梁的受迫振动方程;
S102.结合量纲分析方法得到振动速度预测公式。
5.根据权利要求4所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,其特征在于:所述步骤S101中悬臂梁的受迫振动方程为
式中f(t)为等效爆破荷载,q(H)sin(pt)为等效梁受围岩的作用荷载,只考虑第一阶振型,则得到梁上任意一点的最大速度峰值与整个梁上最大速度峰值的比值为sin(πh/2H),作为表征边坡各点速度与整个边坡最大速度峰值比值特性的无量纲量。
6.根据权利要求5所述的基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法,其特征在于,所述步骤S102中,具体步骤为:爆破振动受到测点与爆源的高差h、岩质高边坡总高度H、测点爆心距R、地表岩体质点振动位移μ、地表岩体质点振动加速度a、岩体的固有振动频率f、最大单响药量Q、岩体密度ρ、振动波传播速度c、以及爆轰时间t的影响,由量纲分析得到
V=Φ(Q,R,H,h,f,c,a,t,μ,ρ),
式中取质量Q,速度c,长度R为基本量纲,则有
不同无量纲量πi的乘积和乘方仍为无量纲量,得到
岩质高边坡的振动速度表示为
CN202010468866.0A 2020-05-28 2020-05-28 基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法 Active CN111753401B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010468866.0A CN111753401B (zh) 2020-05-28 2020-05-28 基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010468866.0A CN111753401B (zh) 2020-05-28 2020-05-28 基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法

Publications (2)

Publication Number Publication Date
CN111753401A CN111753401A (zh) 2020-10-09
CN111753401B true CN111753401B (zh) 2024-04-09

Family

ID=72673590

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010468866.0A Active CN111753401B (zh) 2020-05-28 2020-05-28 基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法

Country Status (1)

Country Link
CN (1) CN111753401B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112611805B (zh) * 2020-12-11 2024-04-05 武汉理工大学 一种基于衰减系数的评价围岩松动圈范围的方法
CN112632759B (zh) * 2020-12-14 2022-09-23 武汉科技大学 一种毫秒延时起爆网路理论单响药量的概率计算方法
CN114723863A (zh) * 2022-01-28 2022-07-08 北京理工大学 露天台阶爆破的振动速度云图预测方法、系统及电子设备

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104390537A (zh) * 2014-09-23 2015-03-04 武汉大学 一种基于爆破振动测试的边坡预裂爆破开挖损伤控制方法
CN104654948A (zh) * 2015-01-16 2015-05-27 武汉科技大学 一种边坡与近接隧道协同爆破施工的爆破控制方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104390537A (zh) * 2014-09-23 2015-03-04 武汉大学 一种基于爆破振动测试的边坡预裂爆破开挖损伤控制方法
CN104654948A (zh) * 2015-01-16 2015-05-27 武汉科技大学 一种边坡与近接隧道协同爆破施工的爆破控制方法

Also Published As

Publication number Publication date
CN111753401A (zh) 2020-10-09

Similar Documents

Publication Publication Date Title
CN111753401B (zh) 基于悬臂梁模型的岩质高边坡爆破坡面振动速度预测方法
CN103150421B (zh) 利用位移监测同时确定抗滑桩桩位和临界深度的方法
CN108538026A (zh) 一种边坡稳定性的多指标综合评价与预警方法
CN103150871A (zh) 利用地下水位与位移实时监测的滑坡预测方法
CN109146142A (zh) 一种回采巷道矿压显现数据的预测方法
CN1991301A (zh) 地下岩移的监测方法
CN108665146A (zh) 大型露天矿山边坡稳定性评价精度确定方法
Xu et al. Blasting vibration control using an improved artificial neural network in the ashele copper mine
CN111784025B (zh) 基于简支梁模型的竖井边墙爆破振动速度预测系统和方法
CN115660487A (zh) 一种隧道光面爆破综合质量评价方法
Wang et al. A TBM cutter life prediction method based on rock mass classification
CN107941137A (zh) 任意倾角钻孔变形量测方法
CN111323192B (zh) 一种基于强夯加固的深埋管道振动衰减效应的测试方法
CN100349014C (zh) 重力勘探数据处理变密度地形校正方法
CN104018482B (zh) 地基极限承载力的测量方法
CN112431236B (zh) 一种新型支护桩危险点的判定方法
CN113804853A (zh) 煤矿瓦斯含量检测方法及系统
CN113031091A (zh) 一种判断高密度直立矿体深部延伸的坑道重力方法
CN105488729A (zh) 一种浮选型多金属尾矿库二次开发利用的资源量调查方法
Wu et al. A new mathematical model for predicting the surface vibration velocity on the step topography
CN110485993A (zh) 一种地震资料套变点预测方法
CN207317727U (zh) 一种任意倾角钻孔变形量测仪
CN108614947A (zh) 一种岩体风化、卸荷的数值判别方法
CN115112771A (zh) 一种测定围岩爆破开挖损伤区的装置及方法
CN113031090B (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