CN113779489A - 一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别方法 - Google Patents

一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别方法 Download PDF

Info

Publication number
CN113779489A
CN113779489A CN202110951695.1A CN202110951695A CN113779489A CN 113779489 A CN113779489 A CN 113779489A CN 202110951695 A CN202110951695 A CN 202110951695A CN 113779489 A CN113779489 A CN 113779489A
Authority
CN
China
Prior art keywords
formula
slurry
entropy
calculating
excavation surface
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
CN202110951695.1A
Other languages
English (en)
Other versions
CN113779489B (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.)
Tongji University
Original Assignee
Tongji 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 Tongji University filed Critical Tongji University
Priority to CN202110951695.1A priority Critical patent/CN113779489B/zh
Publication of CN113779489A publication Critical patent/CN113779489A/zh
Application granted granted Critical
Publication of CN113779489B publication Critical patent/CN113779489B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21D—SHAFTS; TUNNELS; GALLERIES; LARGE UNDERGROUND CHAMBERS
    • E21D9/00—Tunnels or galleries, with or without linings; Methods or apparatus for making thereof; Layout of tunnels or galleries
    • E21D9/06—Making by using a driving shield, i.e. advanced by pushing means bearing against the already placed lining
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0635—Risk analysis of enterprise or organisation activities
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/08—Construction

Landscapes

  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Theoretical Computer Science (AREA)
  • Strategic Management (AREA)
  • Economics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Data Mining & Analysis (AREA)
  • Mining & Mineral Resources (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Pure & Applied Mathematics (AREA)
  • Operations Research (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Marketing (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Game Theory and Decision Science (AREA)
  • Educational Administration (AREA)
  • Development Economics (AREA)
  • Quality & Reliability (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Environmental & Geological Engineering (AREA)
  • Databases & Information Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geology (AREA)
  • Excavating Of Shafts Or Tunnels (AREA)

Abstract

本发明提供一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别方法。首先，通过前期地勘以及连续监测获得泥水盾构推进过程中的重要参数，包括泥水压力Pw，水深hw，上覆各土层参数：土层厚度
Figure DDA0003218158110000016
粘聚力
Figure DDA0003218158110000014
有效粘聚力
Figure DDA0003218158110000015
内摩擦角
Figure DDA0003218158110000011
有效内摩擦角
Figure DDA0003218158110000012
进泥浆密度ρin，排泥浆密度ρout，刀盘转速ω，刀盘扭矩TF等，将所得参数进行归一化处理，并计算各因素信息熵。其次，依据所得信息熵计算相应权重。最后，通过
Figure DDA0003218158110000013
计算系统内熵变情况，并以此判断开挖面状态，针对性地给出掘进策略。

Description

一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别 方法
技术领域
本发明涉及一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别方法，旨在利用盾构施工参数实时判断前方开挖面的稳定性，并针对性做出掘进策略调整。
背景技术
近年来，随着城市、城际交通运输的发展，隧道工程的建设愈发火热，其建设环境逐渐由浅层地下空间向深层地下空间演变。深部地下空间常常面临着深覆土、多地层、高水压的复杂环境，其建设难度远远大于浅层空间。在软土富水地层中，通常采用泥水盾构法施工。由于泥水舱内情况较难悉知，不易观察，工程技术人员对开挖面的前方土体稳定性的判别是一大难点，对潜在风险和事故不易做到准确预测和预防。
传统方法大多通过泥水压力的变化情况判断前方开挖面状态，该类方法忽略了其他可能产生影响的因素，事实上，刀盘转速、扭矩、进浆密度、排浆密度等因素都对泥水舱及开挖面系统的稳定性有重要影响。
发明内容
基于背景技术，本发明提出一种新的泥水盾构开挖面稳定状态判别方法，将可能影响稳定性的各参数，通过信息熵的方式有机地结合在一起，给出更加全面、可靠的判别方法。
本发明的目的在于：将影响开挖面的各项复杂因素有机地整合，跨越开挖面状态不可视的障碍，从系统稳定性的角度，给出开挖面状态的判别方法，为实际泥水盾构施工中的安全提供参考和保障。
为了实现发明目的，采用技术方案如下：
本申请是一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别方法，主要适用于泥水盾构场景，面向盾构机操作人员和施工管理人员。通过对掘进过程重要参数的采集提取及处理，分析出泥水循环系统内的熵变情况(即dS)，并以此为基础判断开挖面的稳定状态，从而有利于在系统状态发生改变的前期尽早干预，降低开挖面失稳的风险。
有益效果：
1、能够解决泥水盾构环境因不可视而导致的开挖面状态难以探察的问题。
2、充分考虑泥水盾构掘进过程中可能影响开挖面稳定性的重要因素，并给出简明清晰的判断准则，有助于相关操作人员和管理人员在系统变化前期对掘进策略进行调整和干预，提高作业安全系数。
附图说明
图1为本发明信息采集处理反馈系统示意图
图2为熵变组成示意图(公式一)
图1中，1——传感器(用于测量开挖面各相关参数)；2——信号接收器(用于接收传感器所测数据)；3——泥水舱；4——工作舱；5——地层环境；6——计算机处理及指令反馈模块。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别方法，其特征在于，具体包括以下步骤：
S1、开挖面的稳定状态由其熵变dS情况决定，如图2所示，熵变可按公式一进行确定：
公式一：
Figure BDA0003218158090000031
其中，S为系统总熵值，从进出浆环节引进的熵流为
Figure BDA0003218158090000032
(假定短时间内进出浆流速差异引发物质熵变很小)，系统内部熵变为diS，T为泥水循环系统内平均温度，Q为系统从外界引进的能量，δQ为系统从外界引进能量的变分。
S2、根据公式二计算δQ：
公式二：
Figure BDA0003218158090000033
其中，Qin为进浆带来的能量，Qout为排浆带走的能量，可由公式三、公式四确定；
公式三：
Figure BDA0003218158090000034
公式四：
Figure BDA0003218158090000035
公式三、公式四中，
Figure BDA0003218158090000041
为进浆管浆液平均温度，
Figure BDA0003218158090000042
为排浆管浆液平均温度，qin为进浆流量，qout为排浆流量，k为经验调整系数，可由试验测定。
S3、计算系统内部熵变为diS：
在泥水盾构掘进过程中，按固定采样频率f收集以下数据：
泥水压力Pw，水深hw，上覆各土层参数：土层厚度
Figure BDA0003218158090000043
粘聚力
Figure BDA0003218158090000044
有效粘聚力
Figure BDA0003218158090000045
内摩擦角
Figure BDA0003218158090000046
有效内摩擦角
Figure BDA0003218158090000047
进泥浆密度ρin，排泥浆密度ρout，刀盘转速ω，刀盘扭矩TF等。并将其整理为描述泥水系统平衡状态的状态参量：{si}＝s1,s2,...,sm。将采集到的数据(包括但不限于上述采集数据种类)按公式五进行参数归一化处理；
公式五：
Figure BDA0003218158090000048
按公式六计算各因素信息熵；
公式六：
Figure BDA0003218158090000049
其中，Si为第i个影响因素的信息熵；
Figure BDA00032181580900000410
为第i个因素第j个数据的归一化值；si为第i个因素所有采样点数据之和，即
Figure BDA00032181580900000411
按公式七计算各影响因素的权重；
公式七：
Figure BDA00032181580900000412
其中wi为第i个因素的熵权值。
按公式八计算开挖面体系总熵值，并取微分；
公式八：
Figure BDA0003218158090000051
S4、根据上述计算结果，进行以下判断并向盾构操作系统返回指令。
当dS＝0时，泥水系统相对平衡稳定，较为安全，可持续推进。
当dS＞0时，泥水系统循环重建，会逐渐转化为另一种相对平衡的状态。
当dS＜＜0并保持一定时间时，泥水系统远离平衡态，开挖面状态易出现突变。
实际上，dS随时间处于起伏波动的状态，应待其在一段时间内保持上述某种状态时，再进行干预，将得到更好的效果。
实施例1：
基于上述理论，给出实施例如下：
某泥水盾构隧道在掘进过程中，通过传感器1实时对关键施工参数进行监测，(如图1所示)。关键施工参数包括：土层厚度
Figure BDA0003218158090000052
粘聚力
Figure BDA0003218158090000053
有效粘聚力
Figure BDA0003218158090000054
内摩擦角
Figure BDA0003218158090000055
有效内摩擦角
Figure BDA0003218158090000056
进泥浆密度ρin，排泥浆密度ρout，刀盘转速ω，刀盘扭矩TF等。
将测量值通过信号接收器2进行接收，并将其整理为描述泥水系统平衡状态的状态参量：{si}＝s1,s2,...,sm。将采集到的数据(包括但不限于上述采集数据种类)按公式五进行参数归一化处理；
公式五：
Figure BDA0003218158090000061
按公式六计算各因素信息熵；
公式六：
Figure BDA0003218158090000062
其中，Si为第i个影响因素的信息熵；
Figure BDA0003218158090000063
为第i个因素第j个数据的归一化值；si为第i个因素所有采样点数据之和，即
Figure BDA0003218158090000064
按公式七计算各影响因素的权重；
公式七：
Figure BDA0003218158090000065
其中wi为第i个因素的熵权值。
按公式八计算开挖面体系总熵值，并取微分；
公式八：
Figure BDA0003218158090000066
根据公式二计算δQ：
公式二：
Figure BDA0003218158090000067
其中，Qin为进浆带来的能量，Qout为排浆带走的能量，可由公式三、公式四确定；
公式三：
Figure BDA0003218158090000068
公式四：
Figure BDA0003218158090000069
公式三、公式四中，
Figure BDA00032181580900000610
为进浆管浆液平均温度，
Figure BDA00032181580900000611
为排浆管浆液平均温度，qin为进浆流量，qout为排浆流量，k为经验调整系数，可由试验测定。
根据公式一计算开挖面的熵变情况：
公式一：
Figure BDA0003218158090000071
在工作舱4中计算机处理及指令反馈模块6通过上述步骤进行处理，绘制熵值随掘进时间的变化情况(如图2所示)。
从图2可以观测到，dS＞0，即泥水盾构系统的总熵值处于稳定上升的状态，从总熵值的曲线可以看出，其发展趋势渐趋平缓，稳定增加。从进出浆环节引进的熵流
Figure BDA0003218158090000072
占比较小，说明泥浆循环控制系统工作稳定，系统内部熵变diS，也随时间推进而缓慢增长，且其趋势渐趋平缓，这表示泥水平衡在新条件下将重新构建，开挖面整体处于较稳定状态，可以继续推进。这一参考意见将反馈给盾构机相关操作人员。

Claims (1)

1.一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别方法，其特征在于，具体包括以下步骤：
S1、开挖面的稳定状态由其熵变dS情况决定，熵变可按公式一进行确定：
公式一：
Figure FDA0003218158080000011
其中，S为系统总熵值，从进出浆环节引进的熵流为
Figure FDA0003218158080000012
系统内部熵变为diS，T为泥水循环系统内平均温度，Q为系统从外界引进的能量，δQ为系统从外界引进能量的变分；
S2、根据公式二计算δQ：
公式二：
Figure FDA0003218158080000013
其中，Qin为进浆带来的能量，Qout为排浆带走的能量，可由公式三、公式四确定；
公式三：
Figure FDA0003218158080000014
公式四：
Figure FDA0003218158080000015
公式三、公式四中，
Figure FDA0003218158080000016
为进浆管浆液平均温度，
Figure FDA0003218158080000017
为排浆管浆液平均温度，qin为进浆流量，qout为排浆流量，k为经验调整系数，可由试验测定；
S3、计算系统内部熵变为diS：
在泥水盾构掘进过程中，按固定采样频率f收集以下数据：
泥水压力Pw，水深hw，上覆各土层参数：土层厚度
Figure FDA00032181580800000110
粘聚力
Figure FDA00032181580800000111
有效粘聚力
Figure FDA00032181580800000112
内摩擦角
Figure FDA0003218158080000018
有效内摩擦角
Figure FDA0003218158080000019
进泥浆密度ρin，排泥浆密度ρout，刀盘转速ω，刀盘扭矩TF等；并将其整理为描述泥水系统平衡状态的状态参量：{si}＝s1,s2,...,sm；将采集到的数据按公式五进行参数归一化处理；
公式五：
Figure FDA0003218158080000021
按公式六计算各因素信息熵；
公式六：
Figure FDA0003218158080000022
其中，Si为第i个影响因素的信息熵；
Figure FDA0003218158080000023
为第i个因素第j个数据的归一化值；si为第i个因素所有采样点数据之和，即
Figure FDA0003218158080000024
按公式七计算各影响因素的权重；
公式七：
Figure FDA0003218158080000025
其中wi为第i个因素的熵权值；
按公式八计算开挖面体系总熵值，并取微分；
公式八：
Figure FDA0003218158080000026
S4、根据上述计算结果，进行以下判断并向盾构操作系统返回指令；
当dS＝0时，泥水系统相对平衡稳定，较为安全，可持续推进；
当dS＞0时，泥水系统循环重建，会逐渐转化为另一种相对平衡的状态；
当dS＜＜0并保持一定时间时，泥水系统远离平衡态，开挖面状态易出现突变；
dS随时间处于起伏波动的状态，应待其在一段时间内保持上述某种状态时，再进行干预，将得到更好的效果。
CN202110951695.1A 2021-08-18 2021-08-18 一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别方法 Active CN113779489B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110951695.1A CN113779489B (zh) 2021-08-18 2021-08-18 一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110951695.1A CN113779489B (zh) 2021-08-18 2021-08-18 一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别方法

Publications (2)

Publication Number Publication Date
CN113779489A true CN113779489A (zh) 2021-12-10
CN113779489B CN113779489B (zh) 2024-02-02

Family

ID=78838303

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110951695.1A Active CN113779489B (zh) 2021-08-18 2021-08-18 一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别方法

Country Status (1)

Country Link
CN (1) CN113779489B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117108294A (zh) * 2023-10-16 2023-11-24 山东济矿鲁能煤电股份有限公司阳城煤矿 一种盾构机故障智能监测系统

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010222805A (ja) * 2009-03-23 2010-10-07 Railway Technical Res Inst 湧水に対する掘削面の安定性評価方法
CN102519790A (zh) * 2011-12-07 2012-06-27 上海大学 确定泥水盾构隧道开挖面支护压力的方法
CN102704937A (zh) * 2011-03-28 2012-10-03 同济大学 一种盾构隧道施工沉降模式的识别处理方法
CN104963691A (zh) * 2015-06-03 2015-10-07 华中科技大学 一种复杂地层条件下土压盾构开挖面稳定性预测控制方法
WO2019140794A1 (zh) * 2018-01-16 2019-07-25 天津大学 一种基于量纲分析的盾构掘进总推力计算方法
US20190244312A1 (en) * 2016-10-14 2019-08-08 Sws Engineering S.P.A. Procedure and system for the calculation of the level of risk in the proximity of the excavation front of an underground work
CN110378574A (zh) * 2019-07-01 2019-10-25 山东大学 水下隧道盾构施工开挖面稳定性评价方法、系统及设备
WO2021073497A1 (zh) * 2019-10-18 2021-04-22 中铁隧道局集团有限公司 一种基于扭推比的破碎地层tbm卡机风险预警方法
CN112837182A (zh) * 2020-12-29 2021-05-25 长春工程学院 一种盾构掘进安全稳定性评价方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2010222805A (ja) * 2009-03-23 2010-10-07 Railway Technical Res Inst 湧水に対する掘削面の安定性評価方法
CN102704937A (zh) * 2011-03-28 2012-10-03 同济大学 一种盾构隧道施工沉降模式的识别处理方法
CN102519790A (zh) * 2011-12-07 2012-06-27 上海大学 确定泥水盾构隧道开挖面支护压力的方法
CN104963691A (zh) * 2015-06-03 2015-10-07 华中科技大学 一种复杂地层条件下土压盾构开挖面稳定性预测控制方法
US20190244312A1 (en) * 2016-10-14 2019-08-08 Sws Engineering S.P.A. Procedure and system for the calculation of the level of risk in the proximity of the excavation front of an underground work
WO2019140794A1 (zh) * 2018-01-16 2019-07-25 天津大学 一种基于量纲分析的盾构掘进总推力计算方法
CN110378574A (zh) * 2019-07-01 2019-10-25 山东大学 水下隧道盾构施工开挖面稳定性评价方法、系统及设备
WO2021073497A1 (zh) * 2019-10-18 2021-04-22 中铁隧道局集团有限公司 一种基于扭推比的破碎地层tbm卡机风险预警方法
CN112837182A (zh) * 2020-12-29 2021-05-25 长春工程学院 一种盾构掘进安全稳定性评价方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杨忠年;刘泉维;刘林胜;张艳;: "泥水盾构开挖面稳定性影响因素敏感性分析", 现代隧道技术, no. 02 *
黄程俣;王炳龙;: "大直径泥水盾构穿江底深槽段开挖面稳定性研究", 华东交通大学学报, no. 01 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117108294A (zh) * 2023-10-16 2023-11-24 山东济矿鲁能煤电股份有限公司阳城煤矿 一种盾构机故障智能监测系统
CN117108294B (zh) * 2023-10-16 2024-01-12 山东济矿鲁能煤电股份有限公司阳城煤矿 一种盾构机故障智能监测系统

Also Published As

Publication number Publication date
CN113779489B (zh) 2024-02-02

Similar Documents

Publication Publication Date Title
CN108301825B (zh) 一种高温高压储层孔隙压力分析方法
CN104695937B (zh) 钻井综合提速优化专家系统
CN110617074B (zh) 一种盾构施工中地面沉降量与掘进参数的关联关系方法
CN106593400A (zh) 一种超低渗致密油藏体积压裂定向井布井方法
CN110671095B (zh) 一种地质地层压力智能随钻软测量方法
CN113779489A (zh) 一种基于熵变原理的泥水平衡盾构系统开挖面稳定性判别方法
CN114723138A (zh) 一种基于epdm数据模型的油气藏动态分析系统和方法
CN110984974B (zh) 基于水侵速度确定有水气藏合理采气速度的方法
CN102536223B (zh) 煤层裂隙孔隙度快速计算方法
CN101936157B (zh) 一种利用测井资料检测高压盐水层孔隙压力的方法
CN113946899A (zh) 一种基于实时地质信息的泥水盾构掘进参数预测方法
CN113137226B (zh) 便携式岩土体力学参数钻探测试系统及设备
CN113653496B (zh) 一种泥水盾构隧道掘进全过程稳定方法
CN112819195A (zh) 一种隧道超前钻探地质精细化预报方法及应用
CN108825204A (zh) 基于随钻工程参数的钻时校正方法
CN113719295B (zh) 一种泥水盾构法隧道掘进全过程稳定性智能控制系统
CN115660235B (zh) 一种一井多用煤层气井生产全过程产量预测方法
CN109917489A (zh) 一种地下承压水位确定的新方法
CN110017118A (zh) 一种防漏堵漏系统
CN116894393B (zh) 一种多参量信息融合的顶板含水层富水性判别方法
Liu et al. Performance assessment and structural design of the atmospheric cutterhead of slurry shield machine
CN116992199A (zh) 一种油页岩地层异常压力预测方法
Sobri et al. Assessment of Tube Well Pumping Test Performance on Different Geological Formation
CN106018754A (zh) 一种重构土壤的精细压实模型的构建方法
Mabile et al. An expert system helps in formation recognition

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