CN109470582B - 露天矿山岩体结构面等效抗剪强度参数分级确定方法 - Google Patents

露天矿山岩体结构面等效抗剪强度参数分级确定方法 Download PDF

Info

Publication number
CN109470582B
CN109470582B CN201811450877.5A CN201811450877A CN109470582B CN 109470582 B CN109470582 B CN 109470582B CN 201811450877 A CN201811450877 A CN 201811450877A CN 109470582 B CN109470582 B CN 109470582B
Authority
CN
China
Prior art keywords
side slope
slope
step side
normal stress
combined
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
CN201811450877.5A
Other languages
English (en)
Other versions
CN109470582A (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.)
Zhejiang University ZJU
Original Assignee
Zhejiang University ZJU
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 Zhejiang University ZJU filed Critical Zhejiang University ZJU
Priority to CN201811450877.5A priority Critical patent/CN109470582B/zh
Publication of CN109470582A publication Critical patent/CN109470582A/zh
Application granted granted Critical
Publication of CN109470582B publication Critical patent/CN109470582B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/24—Investigating strength properties of solid materials by application of mechanical stress by applying steady shearing forces
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014—Type of force applied
    • G01N2203/0025—Shearing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

一种露天矿山岩体结构面等效抗剪强度参数分级确定方法，将基于滑面法向应力分布规律，针对不同等级的露天矿山边坡，在边坡稳定性分级分析的基础上，依据不同级别潜在滑移面上法向应力分布规律，合理确定结构面粘聚力与内摩擦角。本发明显著提高了露天矿山边坡潜在滑移面抗剪强度参数精细化取值的可靠性。

Description

露天矿山岩体结构面等效抗剪强度参数分级确定方法
技术领域
本发明涉及一种露天矿山岩体结构面等效抗剪强度参数分级确定方法，特别的是本发明在边坡稳定性分级分析的基础上，依据不同级别潜在滑移面上法向应力分布规律，提供了一种合理确定结构面粘聚力与内摩擦角的技术方案，显著提高了露天矿山边坡潜在滑移面精细化取值的可靠性，本发明属于工程技术领域。
背景技术
岩质边坡稳定性问题是岩石力学与工程地质学的一个重要研究课题，也是工程实践中迫切要求解决的问题。随着矿业的发展和露天开采深度的加大，露天矿山边坡的稳定性已成为直接关系到矿山安全生产与发展的重大问题。矿山边坡的稳定性受岩体结构面控制，岩体结构面的力学特性决定了矿山边坡发生潜在破坏的可能性，其中结构面抗剪强度研究是岩体稳定性分析的重要内容。许多学者认为结构面直剪试验是确定岩体结构面抗剪强度参数最直接有效的方法。特别是对于有一定重要性的岩土工程，在资金、技术条件允许的范围内，采用直剪试验方法确定结构面抗剪强度至关重要。目前，人们大量采用Mohr-Coulomb线性准则描述岩体的强度，但对于天然粗糙节理面，剪切强度的增长率随正应力的逐渐增加、凸起部分逐渐被剪切和扩张角的减小而减小，从而导致其剪切强度与正应力呈非线性关系，因此，Mohr-Coulomb准则对岩体强度的描述有一定局限性。Barton在仔细研究这些结构面的直剪特性和试验结果的基础上，提出了用于估算岩体结构面抗剪强度的非线性经验公式：
Figure BDA0001886636410000021
式中：
Figure BDA0001886636410000022
为基本摩擦角，JRC为结构面粗糙度系数，JCS为结构面壁岩强度，σn为法向应力。
Hoek和Brown认为，JRC-JCS模型非常适用于岩质边坡工程中的结构面抗剪强度估算，该节理剪切强度准则具有较好的工程应用价值，并在工程实践中广泛使用。
目前，等效Mohr-Coulomb的抗剪强度估算方法(切线等效法)是一种常用的通过JRC-JCS模型获取结构面抗剪强度参数的方法。由于不同法向荷载条件下，结构面的抗剪强度参数随之发生变化。因此，只有在对作用在潜在滑移面上的真实法向受荷状态准确分析的前提下，才能准确确定潜在滑移面上的抗剪强度参数。
截止目前，一些专家、学者已经陆续开展了滑面应力分布形式方面的研究：
Bell(1968)将滑面正应力分布假设为含两个参数的函数。
Yang(2001)将滑面正应力分布用二次函数来逼近。
Zhu(2002)用三次拉格朗日多项式来逼近滑动面的正应力分布。
朱大勇(2004)将滑面正应力通过初始函数和修正函数来联合表示。
杨明成(2004)假定条块底正应力函数为具有两个待定参数的函数。
郑宏(2007)将滑面正应力划分为滑体体积力的贡献分量和条间作用力的贡献分量两个部分。
卢坤林(2012)讨论了滑面应力的构成，并论证了两个正应力分量的贡献程度。
上述有关滑面应力分布形式的研究丰富了人们对滑动面上受荷条件的认识，其主要目的是为了从滑面正应力的角度提出结构面安全系数的评价方法。
然而，截止目前，人们在使用等效Mohr-Coulomb的抗剪强度估算方法(切线等效法)获取结构面抗剪强度参数时，并没有对结构面上真实作用的法向应力条件进行分析，获得的结构面抗剪强度参数可能与真实受力条件下结构面抗剪强度参数存在较大的差异。
此外，露天矿山边坡稳定性评价与公路、铁路、建筑、水利等工程边坡相比，具有鲜明的特色和复杂性。大型露天矿山按构成要素与规模大小可划分为总体边坡、组合台阶边坡、台阶边坡3个层次，需分别评价它们的整体稳定性和局部稳定性。若潜在滑移面是同一组结构面，人们往往采用相同的结构面抗剪强度参数来评价不同规模大小的边坡稳定性。然而，不同规模边坡所对应的滑动面(或潜在滑移面)上受到的法向应力分布特征也完全不同，即作用在结构面上的真实应力条件不同。依据非线性JRC-JCS模型，不同法向应力条件下获取的结构面抗剪强度参数也将完全不同。只有当抗剪强度参数满足Mohr-Coulomb准则时，人们用相同的抗剪强度参数来评价不同规模大小的边坡稳定性才是准确的。然而，当结构面剪切强度与正应力呈非线性关系时，常规方法将不再适用，否则边坡稳定性评价结果将产生较大误差。
发明内容
为了克服JRC-JCS模型获取结构面抗剪强度参数的不足，本发明将基于滑面法向应力分布规律，针对不同等级的露天矿山边坡提出一种露天矿山岩体结构面等效抗剪强度参数分级确定方法，在边坡稳定性分级分析的基础上，依据不同级别潜在滑移面上法向应力分布规律，合理确定结构面粘聚力与内摩擦角，显著提高了露天矿山边坡潜在滑移面精细化取值的可靠性。
本发明解决其技术问题所采用的技术方案是：
一种露天矿山岩体结构面等效抗剪强度参数分级确定方法，包括以下步骤：
(1)对露天矿山边坡进行现场调查，依据边坡的构成要素与规模大小，将露天矿山边坡划分为三个层次：总体边坡、组合台阶边坡、台阶边坡；在对边坡岩体结构特征进行精细化描述的基础上，按照位置匹配性原则和规模匹配性原则，采用赤平投影原理分层次分析露天矿山边坡稳定性、判断露天矿山边坡破坏模式、确定露天矿山边坡破坏模型；
(2)采用矿用三维激光扫描测量系统采集矿山边坡坡面几何形态数据，确定总体边坡、组合台阶边坡、台阶边坡的典型剖面位置，并分别获得总体边坡、组合台阶边坡、台阶边坡的坡面形态函数G1(x)、G2i(x)、G3j(x)，其中i、j分别是组合台阶边坡和台阶边坡的个数。依据现场勘探与潜在滑移面粗糙度测量结果，分别建立总体边坡、组合台阶边坡、台阶边坡的潜在滑移面形态函数S1(x)、S2i(x)、S3j(x)。依据已知的滑体平均密度ρ，计算总体边坡、组合台阶边坡、台阶边坡单位宽度的重量w1(x)、w2i(x)、w3j(x)：
Figure BDA0001886636410000051
(3)根据工程类比，结合露天矿山生产爆破现状和矿岩实际条件确定振动加速度a；对可能作用在总体边坡、组合台阶边坡、台阶边坡的荷载进行调查，确定作用在总体边坡、组合台阶边坡、台阶边坡坡面上的水平荷载函数Qx1(x)、Qx2i(x)、Qx3j(x)以及竖向荷载函数Qy1(x)、Qy2i(x)、Qy3j(x)；
(4)通过如下公式分别计算总体边坡、组合台阶边坡、台阶边坡滑面的正应力σ1、σ2i、σ3j：
Figure BDA0001886636410000052
式中：λ1、λ2、λ3是应力修正系数，取值区间为1.0～1.5；
(5)根据获得的总体边坡、组合台阶边坡、台阶边坡滑面正应力分布函数，分别计算总体边坡、组合台阶边坡、台阶边坡对应滑面上不同位置处的法向应力值；对不同位置处的法向应力值进行统计分析，计算平均法向应力值σ1 *、σ2i *、σ3j *；
(6)根据总体边坡、组合台阶边坡、台阶边坡作用在对应滑移面上的法向应力值的平均值σ1 *、σ2i *、σ3j *，现场精细化测定得到的潜在滑移面粗糙度系数JRCn、壁岩强度JCSn以及岩体结构面基本摩擦角
Figure BDA0001886636410000053
分别计算滑移面上的等效内摩擦角
Figure BDA0001886636410000054
Figure BDA0001886636410000061
Figure BDA0001886636410000062
Figure BDA0001886636410000063
(7)将步骤(6)求解得到的等效内摩擦角
Figure BDA0001886636410000067
带入如下Mohr-Coulomb公式计算总体边坡、组合台阶边坡、台阶边坡作用在对应滑移面的粘聚力c1、c2i、c3j。
Figure BDA0001886636410000065
本发明的有益效果主要表现在：在边坡稳定性分级分析的基础上，依据不同级别潜在滑移面上法向应力分布规律，合理确定结构面粘聚力与内摩擦角，显著提高了露天矿山边坡潜在滑移面精细化取值的可靠性。
具体实施方式
下面对本发明作进一步描述。
一种露天矿山岩体结构面等效抗剪强度参数分级确定方法，包括以下步骤：
(1)对露天矿山边坡进行现场调查，依据边坡的构成要素与规模大小，将露天矿山边坡划分为三个层次：总体边坡、组合台阶边坡、台阶边坡；在对边坡岩体结构特征进行精细化描述的基础上，按照位置匹配性原则和规模匹配性原则，采用赤平投影原理分层次分析露天矿山边坡稳定性、判断露天矿山边坡破坏模式、确定露天矿山边坡破坏模型；
(2)采用矿用三维激光扫描测量系统采集矿山边坡坡面几何形态数据，确定总体边坡、组合台阶边坡、台阶边坡的典型剖面位置，并分别获得总体边坡、组合台阶边坡、台阶边坡的坡面形态函数G1(x)、G2i(x)、G3j(x)，其中i、j分别是组合台阶边坡和台阶边坡的个数。依据现场勘探与潜在滑移面粗糙度测量结果，分别建立总体边坡、组合台阶边坡、台阶边坡的潜在滑移面形态函数S1(x)、S2i(x)、S3j(x)。依据已知的滑体平均密度ρ，计算总体边坡、组合台阶边坡、台阶边坡单位宽度的重量w1(x)、w2i(x)、w3j(x)：
Figure BDA0001886636410000071
(3)根据工程类比，结合露天矿山生产爆破现状和矿岩实际条件确定振动加速度a；对可能作用在总体边坡、组合台阶边坡、台阶边坡的荷载进行调查，确定作用在总体边坡、组合台阶边坡、台阶边坡坡面上的水平荷载函数Qx1(x)、Qx2i(x)、Qx3j(x)以及竖向荷载函数Qy1(x)、Qy2i(x)、Qy3j(x)；
(4)滑面正应力的构成中由滑体体积力和坡面外力对滑面正应力的贡献占主导地位，而总体边坡、组合台阶边坡、台阶边坡条块间作用力可忽略不计；通过如下公式分别计算总体边坡、组合台阶边坡、台阶边坡滑面的正应力σ1、σ2i、σ3j：
Figure BDA0001886636410000081
式中：λ1、λ2、λ3是应力修正系数，取值区间为1.0～1.5；
(5)根据获得的总体边坡、组合台阶边坡、台阶边坡滑面正应力分布函数，分别计算总体边坡、组合台阶边坡、台阶边坡对应滑面上不同位置处的法向应力值(不少于30个位置)；对不同位置处的法向应力值进行统计分析，计算平均法向应力值σ1 *、σ2i *、σ3j *；
(6)根据总体边坡、组合台阶边坡、台阶边坡作用在对应滑移面上的法向应力值的平均值σ1 *、σ2i *、σ3j *，现场精细化测定得到的潜在滑移面粗糙度系数JRCn、壁岩强度JCSn以及岩体结构面基本摩擦角
Figure BDA0001886636410000082
分别计算滑移面上的等效内摩擦角
Figure BDA0001886636410000083
Figure BDA0001886636410000084
Figure BDA0001886636410000085
Figure BDA0001886636410000086
(7)将步骤(6)求解得到的等效内摩擦角
Figure BDA0001886636410000087
带入如下Mohr-Coulomb公式计算总体边坡、组合台阶边坡、台阶边坡作用在对应滑移面的粘聚力c1、c2i、c3j。
Figure BDA0001886636410000091

Claims (1)

1.一种露天矿山岩体结构面等效抗剪强度参数分级确定方法，其特征在于，所述方法包括以下步骤：
(1)对露天矿山边坡进行现场调查，依据边坡的构成要素与规模大小，将露天矿山边坡划分为三个层次：总体边坡、组合台阶边坡、台阶边坡；在对边坡岩体结构特征进行精细化描述的基础上，按照位置匹配性原则和规模匹配性原则，采用赤平投影原理分层次分析露天矿山边坡稳定性、判断露天矿山边坡破坏模式、确定露天矿山边坡破坏模型；
(2)采用矿用三维激光扫描测量系统采集矿山边坡坡面几何形态数据，确定总体边坡、组合台阶边坡、台阶边坡的典型剖面位置，并分别获得总体边坡、组合台阶边坡、台阶边坡的坡面形态函数G1(x)、G2i(x)、G3j(x)，其中i、j分别是组合台阶边坡和台阶边坡的个数，依据现场勘探与潜在滑移面粗糙度测量结果，分别建立总体边坡、组合台阶边坡、台阶边坡的潜在滑移面形态函数S1(x)、S2i(x)、S3j(x)，依据已知的滑体平均密度ρ，计算总体边坡、组合台阶边坡、台阶边坡单位宽度的重量w1(x)、w2i(x)、w3j(x)：
Figure FDA0002604992970000011
(3)根据工程类比，结合露天矿山生产爆破现状和矿岩实际条件确定振动加速度a；对可能作用在总体边坡、组合台阶边坡、台阶边坡的荷载进行调查，确定作用在总体边坡、组合台阶边坡、台阶边坡坡面上的水平荷载函数Qx1(x)、Qx2i(x)、Qx3j(x)以及竖向荷载函数Qy1(x)、Qy2i(x)、Qy3j(x)；
(4)通过如下公式分别计算总体边坡、组合台阶边坡、台阶边坡滑面的法向应力σ1、σ2i、σ3j：
Figure FDA0002604992970000021
式中：λ1、λ2、λ3是应力修正系数，取值区间为1.0～1.5；
(5)根据获得的总体边坡、组合台阶边坡、台阶边坡滑面法向应力分布函数，分别计算总体边坡、组合台阶边坡、台阶边坡对应滑面上不同位置处的法向应力值；对不同位置处的法向应力值进行统计分析，计算平均法向应力值σ1 *、σ2i *、σ3j *；
(6)根据总体边坡、组合台阶边坡、台阶边坡作用在对应滑移面上的法向应力值的平均值σ1 *、σ2i *、σ3j *，现场精细化测定得到的潜在滑移面粗糙度系数JRCn、壁岩强度JCSn以及岩体结构面基本摩擦角
Figure FDA0002604992970000025
分别计算滑移面上的等效内摩擦角
Figure FDA0002604992970000026
Figure FDA0002604992970000022
Figure FDA0002604992970000023
Figure FDA0002604992970000024
(7)将步骤(6)求解得到的等效内摩擦角
Figure FDA0002604992970000032
带入如下Mohr-Coulomb公式计算总体边坡、组合台阶边坡、台阶边坡作用在对应滑移面的粘聚力c1、c2i、c3j；
Figure FDA0002604992970000031
CN201811450877.5A 2018-11-30 2018-11-30 露天矿山岩体结构面等效抗剪强度参数分级确定方法 Active CN109470582B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811450877.5A CN109470582B (zh) 2018-11-30 2018-11-30 露天矿山岩体结构面等效抗剪强度参数分级确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811450877.5A CN109470582B (zh) 2018-11-30 2018-11-30 露天矿山岩体结构面等效抗剪强度参数分级确定方法

Publications (2)

Publication Number Publication Date
CN109470582A CN109470582A (zh) 2019-03-15
CN109470582B true CN109470582B (zh) 2020-10-27

Family

ID=65674681

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811450877.5A Active CN109470582B (zh) 2018-11-30 2018-11-30 露天矿山岩体结构面等效抗剪强度参数分级确定方法

Country Status (1)

Country Link
CN (1) CN109470582B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111157368A (zh) * 2020-02-17 2020-05-15 中国电建集团西北勘测设计研究院有限公司 一种开挖支护条件下岩土抗剪强度参数精细化获取方法
CN113298939A (zh) * 2021-04-08 2021-08-24 湖南科技大学 露天矿山边坡岩体结构面抗剪强度分级确定系统及方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103135128A (zh) * 2013-01-31 2013-06-05 重庆大学 一种地震荷载作用下三维边坡稳定性预测方法
CN107330224A (zh) * 2017-07-24 2017-11-07 中国地质大学（武汉） 一种分条间作用力倾角非假定的边坡稳定性分析条分法
CN108445188A (zh) * 2018-04-16 2018-08-24 浙江大学 基于中智区间函数的岩体结构面粗糙度系数尺寸效应下边坡稳定性表达方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7451666B2 (en) * 2006-01-13 2008-11-18 Jr Johanson, Inc. Apparatus and test procedure for measuring the cohesive, adhesive, and frictional properties of bulk granular solids

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103135128A (zh) * 2013-01-31 2013-06-05 重庆大学 一种地震荷载作用下三维边坡稳定性预测方法
CN107330224A (zh) * 2017-07-24 2017-11-07 中国地质大学（武汉） 一种分条间作用力倾角非假定的边坡稳定性分析条分法
CN108445188A (zh) * 2018-04-16 2018-08-24 浙江大学 基于中智区间函数的岩体结构面粗糙度系数尺寸效应下边坡稳定性表达方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于滑面正应力分布的爆炸震动荷载作用下边坡稳定性分析;刘华丽 等;《岩土力学》;20130831;第34卷;全文 *
大型露天矿山边坡稳定性等精度评价方法;杜时贵;《岩石力学与工程学报》;20180630;第37卷(第6期);第4节至第6节 *

Also Published As

Publication number Publication date
CN109470582A (zh) 2019-03-15

Similar Documents

Publication Publication Date Title
CN109470581B (zh) 露天矿山边坡岩体结构面抗剪强度分级确定方法
CN107917865B (zh) 一种致密砂岩储层多参数渗透率预测方法
KR101642951B1 (ko) Ｇｉｓ 기반 실시간 지진피해 예측 방법
CN108072902B (zh) 一种碳酸盐岩岩石物理模型的建模方法
CN109470582B (zh) 露天矿山岩体结构面等效抗剪强度参数分级确定方法
CN109543326B (zh) 露天矿山边坡岩体结构面法向应力分级确定方法
Bargawa et al. Iron ore resource modeling and estimation using geostatistics
Ceballos et al. Relationship between RMRb and GSI based on in situ data
CN108665146B (zh) 大型露天矿山边坡稳定性评价精度确定方法
CN110940736A (zh) 岩石边坡爆破开挖损伤区岩体力学参数确定方法及系统
Sousa Evaluation of joints in granitic outcrops for dimension stone exploitation
Abe et al. Roughness of fracture surfaces in numerical models and laboratory experiments
Singh et al. Bedform effect on the reorganization of surface and subsurface grain size distribution in gravel bedded channels
CN108776854B9 (zh) 大型露天矿山边坡稳定性等精度评价方法
Ghazi et al. Empirical relationships of shear wave velocity, SPT-N value and vertical effective stress for different soils in Mashhad, Iran
Kulatilake et al. Spectral-based peak-shear-strength criterion for rock joints
Chern et al. Behavior of regular triangular joints under cyclic shearing
Moore Field procedures for the headcut erodibility index
Geier et al. Determination of consolidation properties, selection of computational methods, and estimation of potential error in mine tailings settlement calculations
Feng et al. A new method of block shape classification
lakovlev et al. Determination of joint mechanical parameters for stability analysis in low stress open pit mines
Gasc-Barbier et al. Morphological and mechanical analysis of natural marble joints submitted to shear tests
Yim et al. Critical strain concept-based simple method for pre-evaluation of tunnel face safety using RMR
Rutten et al. Experimental study on rock cutting with a pickpoint
Abe et al. Numerical simulation of a series of flume tests with dry and wet sands by using depth averaged material point method

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