CN109101776A - 基于挡墙侧移监测数据的基坑反分析方法 - Google Patents

基于挡墙侧移监测数据的基坑反分析方法 Download PDF

Info

Publication number
CN109101776A
CN109101776A CN201811169118.1A CN201811169118A CN109101776A CN 109101776 A CN109101776 A CN 109101776A CN 201811169118 A CN201811169118 A CN 201811169118A CN 109101776 A CN109101776 A CN 109101776A
Authority
CN
China
Prior art keywords
soil body
sidesway
standardization
barricade
soil
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.)
Granted
Application number
CN201811169118.1A
Other languages
English (en)
Other versions
CN109101776B (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.)
Institute Of Industrial Technology Chongqing University
Original Assignee
Institute Of Industrial Technology Chongqing University
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 Institute Of Industrial Technology Chongqing University filed Critical Institute Of Industrial Technology Chongqing University
Priority to CN201811169118.1A priority Critical patent/CN109101776B/zh
Publication of CN109101776A publication Critical patent/CN109101776A/zh
Application granted granted Critical
Publication of CN109101776B publication Critical patent/CN109101776B/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/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Structural Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Civil Engineering (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
  • Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)

Abstract

发明提供基于挡墙侧移监测数据的基坑反分析方法。该方法通过开挖初步阶段挡墙侧移监测值与挡墙侧移理论计算值进行比较,逐步反求校核土体参数,把上一开挖阶段求得的新土体参数用于下一步开挖的挡墙侧移预测,如此反复,其预测精度逐步提高,保证开挖到既定深度的挡墙侧移预测值和监测值误差很小。同时该方法将可靠度设计的概念融入到了挡墙侧移的预测中,考虑了土体变异性,更加符合工程实际。本发明计算简便,工程可操作性强,可供基坑有关设计施工人员使用,有效地预测基坑变形,动态调整支护结构设计。本发明可编写成简单的程序软件,适合推广使用。

Description

基于挡墙侧移监测数据的基坑反分析方法
技术领域
本发明涉及土木工程建设领域,特别涉及一种基坑反分析方法。
背景技术
随着城市建设发展,城市人口密度大,而可利用的土地资源越来越有限,从而地下空间的开发利用是具有非常广大的前景的。因此在城市建设中出现了越来越多的深基坑工程。城市中高层建筑深层地下车库、地铁车站、地下商场等深基坑工程呈现高密度、高难度趋势,周围施工环境越来越复杂,对勘察、设计、位移监测与信息处理都要求更高更苛刻。
在基坑工程的建设过程中,基坑的有关监测数据能否有效及时处理,直接影响到基坑后续开挖的支护结构设计的调整,从而关系到基坑工程能否安全施工。现有技术中,有关基坑监测数据进行土体参数反分析的方法多集中于数值分析方法,工程可操作性较弱,计算复杂。
因此亟须提供一种计算简便的基坑反分析方法,供基坑有关设计施工人员使用。
发明内容
本发明的目的是提供基于挡墙侧移监测数据的基坑反分析方法,以解决现有技术中存在的问题。
为实现本发明目的而采用的技术方案是这样的,基于挡墙侧移监测数据的基坑反分析方法,包括以下步骤:
1)在基坑挡墙外侧土体的多个部位进行勘察取样。对样本的土体参数进行统计分析。所述土体参数包括标准化土体剪切强度和标准化土体刚度比得到样本的标准化土体剪切强度比均值标准化土体剪切强度比变异系数COV1、样本的标准化土体刚度比均值和标准化土体刚度比变异系数COV2。其中,Cu为土体不排水剪切强度,σ′v为土体竖向有效应力,E50为土体割线刚度。
2)采用样本的标准化土体剪切强度比均值作为总体标准化土体剪切强度比均值的初始矩估计。采用样本的标准化土体刚度比均值作为总体标准化土体刚度比均值μx2的初始矩估计。得到标准化土体剪切强度比标准差和标准化土体刚度比的标准差
3)对样本土体参数的数据统计值进行正态性检验。当土体参数服从正态分布时,得到标准化土体剪切强度比标准化土体刚度比当土体参数服从对数正态分布时,得到标准化土体剪切强度比x1=exp(λ11x′1),标准化土体刚度比x2=exp(λ22x′2)。其中,x′为当量正态化函数自变量,
4)得到基坑挡墙侧移理论模型如式(1)所示。
式中,B为基坑开挖宽度,m。T为软土层厚度,m。γ为土体容重,kN/m。hi为第i步基坑开挖深度,m。为自然对数形式的系统刚度。EI为挡墙刚度。γw为水的重度kN/m。μw为水位校正因子,取0.8。
5)得到第i步开挖基坑挡墙侧移理论计算值和监测值的误差函数为其中,为第i步开挖基坑挡墙侧移监测值。
令目标函数
6)令误差函数求解x′1和x′2满足目标函数为最小值d(x)min。其中,x′sat1和x′sat2为满足条件的最优解。第i步开挖反演后的标准化土体剪切强度比为xsati和标准化土体刚度比为xsati
7)在基坑开挖过程之中,根据前步工况的量测信息反演确定基坑挡墙侧移理论模型参数将参数值代入基坑挡墙侧移理论模型,对基坑下一步开挖阶段挡墙将要出现的侧移进行预报。根据预报结果,将预报值与实测值进行对比,改正预报参数,采取相应的工程措施实施下一步开挖,直到施工结束。
进一步,步骤3)中,通过极大似然估计法确定样本土体参数的分布类型。
进一步,步骤6)中,当土体参数为正态分布时,反演后的标准化土体剪切强度比标准化土体刚度比当土体参数为对数正态分布时,反演后的标准化土体剪切强度比xsati=exp(λ11x′sati),标准化土体刚度比xsati=exp(λ22x′sati)。
进一步,步骤7)中,使用基于第i步开挖基坑挡墙侧移监测值反演得到的土体参数xsati替换标准化土体剪切强度比xsati替换标准化土体刚度比
本发明的技术效果是毋庸置疑的:
A.在进行基坑土体参数反分析的过程中,将可靠度设计的概念融入到了挡墙侧移的预测中,考虑了土体变异性,更符合工程实际情况;
B.在进行挡墙侧移理论计算值的计算中,采用挡墙侧移预测模型,计算简便;
C.在进行基坑土体参数反分析的过程中,把上一开挖阶段求得的新土体参数用于下一步开挖挡墙侧移的预测,其预测精度不断提高,保证开挖到既定深度的挡墙侧移预测值和监测值误差很小。
附图说明
图1方法流程图;
图2为基坑横截面示意图;
图3为实施例1中土体参数反演变化折线图;
图4为实施例1中基坑挡墙侧移预测值变化图;
图5为实施例2中土体参数反演变化折线图;
图6为实施例2中基坑挡墙侧移预测值变化图。
具体实施方式
下面结合实施例对本发明作进一步说明,但不应该理解为本发明上述主题范围仅限于下述实施例。在不脱离本发明上述技术思想的情况下,根据本领域普通技术知识和惯用手段,做出各种替换和变更,均应包括在本发明的保护范围内。
实施例1:
自然土体是在漫长的地质年代中,经受各种复杂的地质作用而形成,致使土体的性质差异很大。就是在同一地区,同一土体内,也会存在地质特性的强烈变化性,这就决定了土体参数本身的变异性。
现有基坑分析方法通常将土体考虑成各向同性的均质体,难以对工程的实际情况作出客观的评价。本实施例引入以概率为基础的可靠度分析方法,结合实际工程,将土体的标准化土体剪切强度比和标准化土体刚度比作为变量考虑对监测数据反分析,从而有效地预测基坑变形,动态调整支护结构设计。
参见图2,本实施例挡墙贯入硬土层深度为3~5m。本实施例公开基于挡墙侧移监测数据的基坑反分析方法,包括以下步骤:
1)在基坑挡墙外侧土体的多个部位进行勘察取样。对样本的土体参数进行统计分析。参见表1,所述样本的土体参数包括标准化土体剪切强度和标准化土体刚度比统计指标包括标准化土体剪切强度和标准化土体刚度比统计项目包括样本容量、分布区间、均值和变异系数。得到样本的标准化土体剪切强度比均值标准化土体剪切强度比变异系数COV1、样本的标准化土体刚度比均值和标准化土体刚度比变异系数COV2。其中,Cu为土体不排水剪切强度,σ′v为土体竖向有效应力,E50为土体割线刚度。
表1
2)采用样本的标准化土体剪切强度比均值作为总体标准化土体剪切强度比均值的初始矩估计。采用样本的标准化土体刚度比均值作为总体标准化土体刚度比均值的初始矩估计。得到标准化土体剪切强度比标准差和标准化土体刚度比的标准差变异系数COV1=COV2=0.15。
3)通过极大似然估计法判断样本土体参数总体服从正态分布。则标准化土体剪切强度比标准化土体刚度比
4)将基坑开挖宽度B=21m,软土层厚度T=22m,土体容重γ=17.3kN/m;第一步基坑开挖深度h1=2m,基坑支护系统刚度标准化土体剪切强度比标准化土体刚度比水位校正因子μw=0.8,代入基坑挡墙侧移理论计算公式,则可以得到基坑挡墙侧移理论计算值:
5)第一步开挖基坑挡墙侧移监测值为则第一步开挖基坑挡墙侧移理论计算值和监测值的误差函数为令目标函数在本实施例中,基坑监测数据如表2所示。
表2
6)令误差函数g(x)=0,求解x′1和x′2满足目标函数为最小值d(x)min,此步可以由计算机快速解出。满足上述条件的最优解为x′1=1.27,x′2=0.57。则经过第一步开挖基坑挡墙侧移监测值反演后的标准化土体剪切强度比和标准化土体刚度比
7)将基于第一步开挖基坑挡墙侧移监测值反演得到的标准化土体剪切强度比x1=0.298和标准化土体刚度比x2=162.8替换初始的标准化土体剪切强度比和标准化土体刚度比
8)重复步骤,直到基坑开挖结束,根据基坑每步开挖的挡墙侧移监测值不断反演土体参数,表3即为基于每一步开挖基坑挡墙侧移监测值反演得到的标准化土体剪切强度比和标准化土体刚度比,图3为标准化土体剪切强度比和标准化土体刚度比反演变化折线图,从而得到基于每一步开挖基坑挡墙侧移监测值反演得到的基坑挡墙侧移预测值如图4,从图4可以看出随着反演次数的不断增加,基坑挡墙侧移预测值精度越来越高。
表3
实施例2:
本实施例挡墙贯入硬土层深度为3~5m。本实施例公开基于挡墙侧移监测数据的基坑反分析方法,包括以下步骤:
1)在基坑挡墙外侧土体的多个部位进行勘察取样。对样本的土体参数进行统计分析。参见表4,所述土体参数包括标准化土体剪切强度和标准化土体刚度比统计指标包括标准化土体剪切强度和标准化土体刚度比统计项目包括样本容量、分布区间、均值和变异系数。得到样本的标准化土体剪切强度比均值标准化土体剪切强度比变异系数COV1、样本的标准化土体刚度比均值和标准化土体刚度比变异系数COV2
表4
2)采用样本的标准化土体剪切强度比均值作为总体标准化土体剪切强度比均值的初始矩估计。采用样本的标准化土体刚度比均值作为总体标准化土体刚度比均值的初始矩估计。得到标准化土体剪切强度比标准差和标准化土体刚度比的标准差变异系数COV1=COV2=0.15。
3)本实施例样本土体参数服从对数正态分布。则标准化土体剪切强度比x1=exp(λ11x′1),标准化土体刚度比x2=exp(λ22x′2),
其中 可以得到标准化土体剪切强度比x1=exp(λ11x′1)=exp(-1.525+0.149x′1),标准化土体刚度比x2=exp(λ22x′2)=exp(5.000+0.149x′2)。
4)将基坑开挖宽度B=12.2m,软土层厚度T=22m,土体容重γ=19.1kN/m;第一步基坑开挖深度h1=4.7m,基坑支护系统刚度标准化土体剪切强度比x1=exp(λ11x′1),标准化土体刚度比x2=exp(λ22x′2),水位校正因子μw=0.8,代入基坑挡墙侧移理论计算公式,则可以得到基坑挡墙侧移理论计算值:
5)第一步开挖基坑挡墙侧移监测值为则第一步开挖基坑挡墙侧移理论计算值和监测值的误差函数为 令目标函数在本实施例中,基坑监测数据如表5所示。
表5
Stage 1 Stage 2 Stage 3 Stage 4
开挖深度/h(m) 4.7 7.9 11.3 12.2
挡墙侧移/δ<sub>h,M</sub>(mm) 143.9 156.0 163.0 162.1
6)令误差函数g(x)=0,求解x′1和x′2满足目标函数为最小值d(x)min,此步可以由计算机快速解出。满足上述条件的最优解为x′1=-1.16,x′2=3.41。则经过第一步开挖基坑挡墙侧移监测值反演后的标准化土体剪切强度比x1=exp(λ11x′1)=0.131和标准化土体刚度比x2=exp(λ22x′2)=124.7。
7)将基于第一步开挖基坑挡墙侧移监测值反演得到的标准化土体剪切强度比x1=0.131和标准化土体刚度比x2=124.7替换初始的标准化土体剪切强度比和标准化土体刚度比
8)重复步骤,直到基坑开挖结束,根据基坑每步开挖的挡墙侧移监测值不断反演土体参数,表6即为基于每一步开挖基坑挡墙侧移监测值反演得到的标准化土体剪切强度比和标准化土体刚度比,图5为标准化土体剪切强度比和标准化土体刚度比反演变化折线图,从而得到基于每一步开挖基坑挡墙侧移监测值反演得到的基坑挡墙侧移预测值如图6,从图6可以看出随着反演次数的不断增加,基坑挡墙侧移预测值精度越来越高。
表6

Claims (4)

1.基于挡墙侧移监测数据的基坑反分析方法,其特征在于,包括以下步骤:
1)在基坑挡墙外侧土体的多个部位进行勘察取样;对样本的土体参数进行统计分析;所述土体参数包括标准化土体剪切强度和标准化土体刚度比得到样本的标准化土体剪切强度比均值标准化土体剪切强度比变异系数COV1、样本的标准化土体刚度比均值和标准化土体刚度比变异系数COV2;其中,Cu为土体不排水剪切强度,σ′v为土体竖向有效应力,E50为土体割线刚度;
2)采用样本的标准化土体剪切强度比均值作为总体标准化土体剪切强度比均值的初始矩估计;采用样本的标准化土体刚度比均值作为总体标准化土体刚度比均值的初始矩估计;得到标准化土体剪切强度比标准差和标准化土体刚度比的标准差
3)对样本土体参数的数据统计值进行正态性检验;当土体参数服从正态分布时,得到标准化土体剪切强度比标准化土体刚度比当土体参数服从对数正态分布时,得到标准化土体剪切强度比x1=exp(λ11x′1),标准化土体刚度比x2=exp(λ22x′2);其中,x′为当量正态化函数自变量,
4)得到基坑挡墙侧移理论模型如式(1)所示;
式中,B为基坑开挖宽度,m;T为软土层厚度,m;γ为土体容重,kN/m;hi为第i步基坑开挖深度,m;为自然对数形式的系统刚度;EI为挡墙刚度;γw为水的重度kN/m;μw为水位校正因子,取0.8;
5)得到第i步开挖基坑挡墙侧移理论计算值和监测值的误差函数为其中,为第i步开挖基坑挡墙侧移监测值;
令目标函数
6)令误差函数求解x1′和x′2满足目标函数为最小值d(x)min;其中,x′sat1和x′sat2为满足条件的最优解;第i步开挖反演后的标准化土体剪切强度比为xsati和标准化土体刚度比为xsati
7)在基坑开挖过程之中,根据前步工况的量测信息反演确定基坑挡墙侧移理论模型参数将参数值代入基坑挡墙侧移理论模型,对基坑下一步开挖阶段挡墙将要出现的侧移进行预报;根据预报结果,将预报值与实测值进行对比,改正预报参数,采取相应的工程措施实施下一步开挖,直到施工结束。
2.根据权利要求1所述的基于挡墙侧移监测数据的基坑反分析方法,其特征在于:步骤3)中,通过极大似然估计法确定样本土体参数的分布类型。
3.根据权利要求1所述的基于挡墙侧移监测数据的基坑反分析方法,其特征在于:步骤6)中,当土体参数为正态分布时,反演后的标准化土体剪切强度比标准化土体刚度比当土体参数为对数正态分布时,反演后的标准化土体剪切强度比xsati=exp(λ11x′sati),标准化土体刚度比xsati=exp(λ22x′sati)。
4.根据权利要求1所述的基于挡墙侧移监测数据的基坑反分析方法,其特征在于:步骤7)中,使用基于第i步开挖基坑挡墙侧移监测值反演得到的土体参数xsati替换标准化土体剪切强度比xsati替换标准化土体刚度比
CN201811169118.1A 2018-10-08 2018-10-08 基于挡墙侧移监测数据的基坑反分析方法 Active CN109101776B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811169118.1A CN109101776B (zh) 2018-10-08 2018-10-08 基于挡墙侧移监测数据的基坑反分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811169118.1A CN109101776B (zh) 2018-10-08 2018-10-08 基于挡墙侧移监测数据的基坑反分析方法

Publications (2)

Publication Number Publication Date
CN109101776A true CN109101776A (zh) 2018-12-28
CN109101776B CN109101776B (zh) 2022-12-09

Family

ID=64868358

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811169118.1A Active CN109101776B (zh) 2018-10-08 2018-10-08 基于挡墙侧移监测数据的基坑反分析方法

Country Status (1)

Country Link
CN (1) CN109101776B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110397092A (zh) * 2019-05-14 2019-11-01 洛阳市轨道交通有限责任公司 砂卵石基坑位移预测方法
CN111705808A (zh) * 2020-06-30 2020-09-25 河海大学 一种适用于悬挂式深基坑工程的土体参数动态反演分析方法
CN116522742A (zh) * 2023-07-05 2023-08-01 清华大学 基于深度学习的深基坑多层土力学参数反分析方法及组件
CN117150880A (zh) * 2023-07-19 2023-12-01 华东交通大学 一种基坑支护结构的变形预测以及逆向设计方法及系统
CN117150880B (zh) * 2023-07-19 2024-06-11 华东交通大学 一种基坑支护结构的变形预测以及逆向设计方法及系统

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104077451A (zh) * 2014-07-03 2014-10-01 中铁上海工程局集团有限公司 一种用于深厚软土地铁基坑土体参数反演分析的方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104077451A (zh) * 2014-07-03 2014-10-01 中铁上海工程局集团有限公司 一种用于深厚软土地铁基坑土体参数反演分析的方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
CHAO-FENG ZENG ET AL.: "Responses of retaining wall and surrounding ground to pre-excavation dewatering in an alternated multi-aquifer-aquitard system", 《JOURNAL OF HYDROLOGY》 *
WENGANG ZHANG ET AL.: "Influence of groundwater drawdown on excavation responsese - A case history in Bukit Timah granitic residual soils", 《JOURNAL OF ROCK MECHANICS AND GEOTECHNICAL ENGINEERING》 *
杜利峰 等: "辽宁营口沿海地区某深基坑工程数值模拟的分析", 《高原地震》 *
潘洪科 等: "基坑工程稳定性分析与预测的两种方法", 《矿业研究与开发》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110397092A (zh) * 2019-05-14 2019-11-01 洛阳市轨道交通有限责任公司 砂卵石基坑位移预测方法
CN111705808A (zh) * 2020-06-30 2020-09-25 河海大学 一种适用于悬挂式深基坑工程的土体参数动态反演分析方法
CN116522742A (zh) * 2023-07-05 2023-08-01 清华大学 基于深度学习的深基坑多层土力学参数反分析方法及组件
CN117150880A (zh) * 2023-07-19 2023-12-01 华东交通大学 一种基坑支护结构的变形预测以及逆向设计方法及系统
CN117150880B (zh) * 2023-07-19 2024-06-11 华东交通大学 一种基坑支护结构的变形预测以及逆向设计方法及系统

Also Published As

Publication number Publication date
CN109101776B (zh) 2022-12-09

Similar Documents

Publication Publication Date Title
CN111441758B (zh) 页岩油气甜点区的预测方法及装置
US10324229B2 (en) System and method of pore type classification for petrophysical rock typing
KR101618713B1 (ko) 지반정보에 기초한 3차원 공간 모델링 방법
CN109101776A (zh) 基于挡墙侧移监测数据的基坑反分析方法
CN105804730B (zh) 使用历史井数据的资源识别的方法和系统
US10620339B2 (en) Static earth model calibration methods and systems using tortuosity evaluations
CN109375283B (zh) 一种砂岩储层3d渗透率演化史的分析方法
CN109241627A (zh) 概率分级的动态支护方法及自动设计支护方案的装置
US20150285950A1 (en) Systems and Methods for Selecting Facies Model Realizations
Bashari et al. Estimation of deformation modulus of rock masses by using fuzzy clustering-based modeling
Gong et al. A modified rock mass classification system for TBM tunnels and tunneling based on the HC method of China
CN112182694A (zh) 一种基于bim系统的灌浆工程全过程动态分析方法
CN115017791A (zh) 隧道围岩级别识别方法和装置
CN115586086A (zh) 一种基于大数据的井壁失稳分析方法
Paraskevopoulou et al. Assessing the failure potential of tunnels and the impacts on cost overruns and project delays
Wang et al. Data-driven analysis of soil consolidation with prefabricated vertical drains considering stratigraphic variation
CN116862432A (zh) 一种隧道浅埋段施工方案动态调整方法及系统
Yang et al. Quantitative geomorphological constraints on the landform evolution of the current active boundary of the northeastern Tibetan Plateau
CN110414085A (zh) 已开发断块油藏原始油水界面确定方法
CN114066271A (zh) 一种隧道突水灾害监测管理系统
CN111815769B (zh) 逆冲推覆构造带构造的建模方法、计算设备及存储介质
Kovacevic et al. The use of neural networks to develop CPT correlations for soils in northern Croatia
CN116559938A (zh) 一种油气藏裂缝模型建立方法及电子设备
Huang et al. Prediction of design typhoon wind speeds and profiles using refined typhoon wind field model
Ding et al. Quantitative estimation of TBM disc cutter wear from in-situ parameters by optimization algorithm improved back-propagation neural network: A case study of a metro tunnel in Guangzhou, China

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