CN113821862B - 一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法 - Google Patents

一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法 Download PDF

Info

Publication number
CN113821862B
CN113821862B CN202111256225.XA CN202111256225A CN113821862B CN 113821862 B CN113821862 B CN 113821862B CN 202111256225 A CN202111256225 A CN 202111256225A CN 113821862 B CN113821862 B CN 113821862B
Authority
CN
China
Prior art keywords
river
sediment
average
slope angle
blocking building
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
CN202111256225.XA
Other languages
English (en)
Other versions
CN113821862A (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.)
Hangzhou Jinghang Canal Second Channel Construction Investment Co ltd
Hohai University HHU
Nanjing Hydraulic Research Institute of National Energy Administration Ministry of Transport Ministry of Water Resources
Original Assignee
Hangzhou Jinghang Canal Second Channel Construction Investment Co ltd
Hohai University HHU
Nanjing Hydraulic Research Institute of National Energy Administration Ministry of Transport Ministry of Water Resources
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 Hangzhou Jinghang Canal Second Channel Construction Investment Co ltd, Hohai University HHU, Nanjing Hydraulic Research Institute of National Energy Administration Ministry of Transport Ministry of Water Resources filed Critical Hangzhou Jinghang Canal Second Channel Construction Investment Co ltd
Priority to CN202111256225.XA priority Critical patent/CN113821862B/zh
Publication of CN113821862A publication Critical patent/CN113821862A/zh
Application granted granted Critical
Publication of CN113821862B publication Critical patent/CN113821862B/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
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/08Probabilistic or stochastic CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/08Fluids
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Civil Engineering (AREA)
  • Computational Mathematics (AREA)
  • Structural Engineering (AREA)
  • Architecture (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Fluid Mechanics (AREA)
  • Mathematical Physics (AREA)
  • Revetment (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法,主要应用于水利工程、航道工程以及河道治理等领域,该方法认河床的泥沙粒径与水流的大小密切相关,且河道中水流的临底流速符合高斯分布。根据水流的分布概率,推导得到河床泥沙沙粒径的概率分布。在此基础上,统计河道的特征流量、特征水位、和河道特性参数、河床阻力系数等,并通过数值积分得到拦河建筑物坝前泥沙淤积体的平均粒径。通过室内试验,建立泥沙粒径与泥沙水坡角的关系式,采用该关系式得到拦河建筑物坝前泥沙淤积体泥沙的水下坡角,该角度即为拦河建筑物坝前泥沙淤积体的平均坡角。

Description

一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法
技术领域
本发明是一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法,主要应用于水利工程、航道工程以及河道治理等领域。
背景技术
泥沙淤积问题是水利工程中需要解决的一项重大课题。河道中修船闸、大坝等建拦河建筑物后,改变的河道的流量过程,降低了流速,容易在上游发生淤积。对于拦河建筑物前的泥沙淤积形态,往往受到拦河建筑的运用方式决定,如:闸门的开启方式、排沙孔的控制运用规程等。拦河建筑物前的泥沙淤积体形态的预测对于建筑物功能的发挥以及后期的调度有重要的作用。
目前,大量的研究主要集中于水库坝前的泥沙淤积形态,常称为冲刷漏斗。例如,有相关学者调查过黄河上游刘家峡、八盘峡、青铜峡以及盐锅峡等水库坝前冲刷漏斗的形态,但就其坡度的确定尚无很有效的方法,这对于水库功能的充分发挥较为不利。对于部分通航河段,拦河建筑物前的泥沙淤积还会影响通航效率。
拦河建筑物前泥沙淤积体的平均坡角的影响因素主要包括:床沙的粒径、建筑物前的水流条件。如果能预估拦河建筑物修建后其上游的床沙粒径和代表性的水流条件,则可以计算得到拦河建筑物前泥沙淤积体最大的平均坡角,为工程设计、管理提供依据。因此,本发明分别提出拦河建筑物前泥沙淤积体粒径的确定方法,在此基础上计算泥沙淤积体平均坡角,该方法简单实用,可为相关设计提供依据。
发明内容
发明目的:本发明的目的在于提供一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法,在实际情况中,需要收集河道的水文特征参数,在此基础上预测坝前泥沙颗粒的平均粒径,并针对性开展室内试验,建立泥沙粒径与水下坡角的关系,综合确定坝前泥沙淤积体的平均坡角。
技术方案:为了实现上述目的,本发明的一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法,包括如下步骤:
步骤1.1,统计拦河建筑物上游河道的多年平均日流量和拦河建筑物前的平均水位,将其作为拦河建筑物的特征下泄流量Q和特征水位Z;
步骤1.2,统计拦河建筑物上游附近河道的平均宽度B和河床平均高程Z0
步骤1.3,根据流速的概率分布,预测拦河建筑物前淤积泥沙的平均粒径D0
步骤1.4,选取不同粒径的泥沙,开展均匀沙条件下的水下坡角实验,拟合得到泥沙粒径与泥沙水下坡角的关系;
步骤1.5,根据步骤1.4中的水下坡角与泥沙粒径的关系式,计算得到该粒径的水下坡角,即得到拦河建筑物坝前泥沙淤积体的平均坡角。
进一步地,所述的预测拦河建筑物前淤积泥沙的平均粒径D0包括如下步骤:
步骤2.1,拦河建筑物坝前泥沙淤积体颗粒运动过程中粒径与水流紊动密切相关,符合如下概率分布
式中,D为拦河建筑物坝前泥沙淤积体的泥沙粒径,为泥沙颗粒附近瞬时的平均流速,ρs为泥沙密度,ρ为水体密度,g为重力加速度,u*为摩阻流速,/>R为水力半径,J为能坡。
步骤2.2,根据泥沙粒径的概率分布确定平均粒径为:
式中,D0为拦河建筑物前淤积泥沙的平均粒径;φ(D)为淤积泥沙的概率分布函数。
积分后可得到泥沙的平均粒径表达式:
式中R为水力半径,u*为摩阻流速,J为能坡,ρs为泥沙密度,ρ为水体密度,g为重力加速度。
进一步地,所述的拦河建筑物的特征下泄流量Q采用下式计算:
则变为:式中n为河道糙率,A为特征下泄流量Q条件下坝前过水面积A=B(Z-Z0),水力半径R=Z-Z0
进一步地,根据所述的特征下泄流量Q计算得到拦河建筑物前泥沙淤积体的平均粒径:
进一步地,所述的拦河建筑物坝前泥沙淤积体的平均坡角的计算过程如下:
在实验室中拟合出均匀沙粒径与水下坡角的关系,再根据该关系式确定泥沙堆积体平均坡角,均匀沙粒径与水下坡角的关系为:
y=34+2.1lnD0
式中,y为泥沙水下坡角,D0的床面泥沙平均粒径;拦河建筑物前泥沙淤积体平均坡角的计算方法为:
有益效果:与现有技术相比,本发明的方法认河床的泥沙粒径与水流的大小密切相关,且河道中水流的临底流速符合高斯分布。根据水流的分布概率,推导得到河床泥沙沙粒径的概率分布。在此基础上,统计河道的特征流量、特征水位、和河道特性参数、河床阻力系数等,并通过数值积分得到拦河建筑物坝前泥沙淤积体的平均粒径。通过室内试验,建立泥沙粒径与泥沙水下坡角的关系式,采用该关系式得到拦河建筑物坝前泥沙淤积体泥沙的水下坡角,该角度即为拦河建筑物坝前泥沙淤积体的平均坡角。
附图说明
图1是本发明一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法流程图;
图2是本发明泥沙粒径与泥沙水下坡角的拟合关系曲线。
具体实施方式
下面结合附图和具体实施例对本发明作更进一步的说明。
如图1所示,一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法,包括如下步骤:
步骤1.1,统计拦河建筑物上游河道的多年平均日流量和拦河建筑物前的平均水位,将其作为拦河建筑物的特征下泄流量Q(单位:m3/s)和特征水位Z(单位:m);
步骤1.2,统计拦河建筑物上游附近河道的平均宽度B和河床平均高程Z0
步骤1.3,根据流速的概率分布,预测拦河建筑物前淤积泥沙的平均粒径D0(单位:m);
步骤1.4,选取不同粒径的泥沙,开展均匀沙条件下的水下坡角实验,拟合得到泥沙粒径与水下坡角的关系;
步骤1.5,根据步骤1.4中的水下坡角与泥沙粒径的关系式,计算得到该粒径的水下坡角,即得到拦河建筑物坝前泥沙淤积体的平均坡角。
确定拦河建筑物坝前泥沙淤积体的粒径的步骤为:
步骤2.1,拦河建筑物坝前泥沙淤积体颗粒运动过程中粒径与水流紊动密切相关,符合如下概率分布
式中,D为拦河建筑物坝前泥沙淤积体的泥沙粒径,为泥沙颗粒附近瞬时的平均流速,ρs为泥沙密度,ρ为水体密度,g为重力加速度,u*为摩阻流速,/>R为水力半径,J为能坡。
步骤2.2,根据泥沙粒径的概率分布确定平均粒径为:
式中,D0为拦河建筑物前淤积泥沙的平均粒径;φ(D)为淤积泥沙的概率分布函数。
积分后可得到泥沙的平均粒径表达式:
式中R为水力半径,u*为摩阻流速,J为能坡,ρs为泥沙密度,ρ为水体密度,g为重力加速度。
步骤2.3:特征流量下拦河建筑物前的特征下泄流量Q采用下式计算:
则可变为:式中n为河道糙率,可根据河道的特性确定,无实测资料条件下可采用/>来进行计算,A为特征下泄流量Q条件下坝前过水面积A=B(Z-Z0),水力半径R可近似为R=Z-Z0
步骤2.4,将步骤2.3中的计算结果代入步骤2.2即可计算得到拦河建筑物前泥沙淤积体的平均粒径。
在实验室中拟合出均匀沙粒径与水下坡角的关系,再根据该关系式确定泥沙堆积体平均坡角。均匀沙粒径与水下坡角的关系为:
y=34+2.1lnD0
式中,y为泥沙水下坡角,D0的床面泥沙平均粒径。拦河建筑物前泥沙淤积体平均坡角的计算方法为:
如图2所示,为泥沙颗粒的水下坡角的拟合公式,只要能预测拦河建筑物前泥沙淤积体的平均粒径,即可得到泥沙淤积体的坡度,能为拦河建筑物的运用方式和设计提供参考。天然河道拦河建筑物前泥沙淤积体的粒径与水流条件密切相关,具有不确定性,难以预测。本发明从水流的特征出发,认为河床泥沙的粒径与水流紊动的概率分布一致,通过水流的紊动概率分布得到了河床泥沙粒径的概率分布,并通过积分预测了平衡条件下河建筑物前泥沙淤积体的平均粒径,通过该方法得到了拦河建筑物坝前泥沙淤积体的平均坡角,该方法简单实用,可以进行快速预测。
实施例
一种拦河建筑物坝前泥沙淤积体平均坡角方法,包括如下步骤:
步骤1,统计拦河建筑物上游河道的多年平均日流量350m3/s和拦河建筑物前的平均水位3m,将其作为拦河建筑物的特征下泄流量和特征水位;
步骤2,统计拦河建筑物上游附近河道的平均宽度40m和河床平均高程-4m;
步骤3,根据河道两侧植物覆盖度和阻水特性,确定河道的阻力系数n=0.02;
步骤4,根据上述步骤中的参数取值,预测拦河建筑物坝前泥沙的平均粒径:
步骤5,选取不同粒径的泥沙,开展均匀沙条件下的水下坡角实验,拟合得到泥沙粒径与水下坡角的关系,本例得到的水下坡角关系为:
y=34+2.1lnD
步骤6,根据步骤5中的水下坡角与泥沙粒径的关系式,计算得到该粒径的水下坡角为34.6度,即河建筑物坝前泥沙淤积体平均坡角为34.6度。

Claims (4)

1.一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法,其特征在于,包括如下步骤:
步骤1.1,统计拦河建筑物上游河道的多年平均日流量和拦河建筑物前的平均水位,将其作为拦河建筑物的特征下泄流量Q和特征水位Z;
步骤1.2,统计拦河建筑物上游附近河道的平均宽度B和河床平均高程Z0
步骤1.3,根据流速的概率分布,预测拦河建筑物前淤积泥沙的平均粒径D0
步骤1.4,选取不同粒径的泥沙,开展均匀沙条件下的水下坡角实验,拟合得到泥沙粒径与泥沙水下坡角的关系;
步骤1.5,根据步骤1.4中的水下坡角与泥沙粒径的关系式,计算得到该粒径的水下坡角,即得到拦河建筑物坝前泥沙淤积体的平均坡角;
所述的拦河建筑物坝前泥沙淤积体的平均坡角的计算过程如下:
拟合出均匀沙粒径与水下坡角的关系,再根据该关系式确定泥沙堆积体平均坡角,均匀沙粒径与水下坡角的关系为:
y=34+2.1lnD0
式中,y为泥沙水下坡角,D0的床面泥沙平均粒径;拦河建筑物坝前泥沙淤积体的平均坡角的计算方法为:
其中,ρs为泥沙密度,ρ为水体密度,g为重力加速度,n为河道糙率,Q为特征下泄流量,B为河道的平均宽度,Z为特征水位,Z0为河床平均高程。
2.根据权利要求1所述的一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法,其特征在于,所述的预测拦河建筑物前淤积泥沙的平均粒径D0包括如下步骤:
步骤2.1,拦河建筑物坝前泥沙淤积体颗粒运动过程中粒径与水流紊动密切相关,符合如下概率分布
式中,D为拦河建筑物坝前泥沙淤积体的泥沙粒径,为泥沙颗粒附近瞬时的平均流速,ρs为泥沙密度,ρ为水体密度,g为重力加速度,u*为摩阻流速,/>R为水力半径,J为能坡;
步骤2.2,根据泥沙粒径的概率分布确定平均粒径为:
式中,D0为拦河建筑物前淤积泥沙的平均粒径;φ(D)为淤积泥沙的概率分布函数;积分后可得到泥沙的平均粒径表达式:
式中R为水力半径,u*为摩阻流速,J为能坡,ρs为泥沙密度,ρ为水体密度,g为重力加速度。
3.根据权利要求2所述的一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法,其特征在于,所述的拦河建筑物的特征下泄流量Q采用下式计算:
则变为:式中n为河道糙率,A为特征下泄流量Q条件下坝前过水面积A=B(Z-Z0),水力半径R=Z-Z0
4.根据权利要求3所述的一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法,其特征在于,根据所述的特征下泄流量Q计算得到拦河建筑物前泥沙淤积体的平均粒径:
CN202111256225.XA 2021-10-27 2021-10-27 一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法 Active CN113821862B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111256225.XA CN113821862B (zh) 2021-10-27 2021-10-27 一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111256225.XA CN113821862B (zh) 2021-10-27 2021-10-27 一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法

Publications (2)

Publication Number Publication Date
CN113821862A CN113821862A (zh) 2021-12-21
CN113821862B true CN113821862B (zh) 2023-08-08

Family

ID=78918911

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111256225.XA Active CN113821862B (zh) 2021-10-27 2021-10-27 一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法

Country Status (1)

Country Link
CN (1) CN113821862B (zh)

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102359862A (zh) * 2011-08-12 2012-02-22 河海大学 粉沙质和淤泥质海岸泥沙运动数值模拟方法
CN106960105A (zh) * 2017-03-31 2017-07-18 河海大学 一种层理淤积河床上桥墩局部冲刷深度计算方法
CN110046469A (zh) * 2019-05-13 2019-07-23 水利部交通运输部国家能源局南京水利科学研究院 多约束条件下水电站坝前河床冲淤变形的计算方法
CN111008357A (zh) * 2019-12-19 2020-04-14 中国水利水电第七工程局有限公司 一种污水管道中泥沙平均淤堵高度的确定方法
CN111350162A (zh) * 2020-03-03 2020-06-30 河海大学 一种高滩深槽水库坝前恢复库容的挖槽设计方法
CN113159994A (zh) * 2021-04-25 2021-07-23 黄河勘测规划设计研究院有限公司 多泥沙河流的船闸淤积处理方法及装置
CN113297753A (zh) * 2021-06-24 2021-08-24 四川大学 山区河流透过性坝沿程水深预测方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11859488B2 (en) * 2018-11-29 2024-01-02 Bp Exploration Operating Company Limited DAS data processing to identify fluid inflow locations and fluid type
CN111507540B (zh) * 2020-04-28 2022-04-26 长江水利委员会长江科学院 一种资料贫乏地区长系列泥沙数据综合重构方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102359862A (zh) * 2011-08-12 2012-02-22 河海大学 粉沙质和淤泥质海岸泥沙运动数值模拟方法
CN106960105A (zh) * 2017-03-31 2017-07-18 河海大学 一种层理淤积河床上桥墩局部冲刷深度计算方法
CN110046469A (zh) * 2019-05-13 2019-07-23 水利部交通运输部国家能源局南京水利科学研究院 多约束条件下水电站坝前河床冲淤变形的计算方法
CN111008357A (zh) * 2019-12-19 2020-04-14 中国水利水电第七工程局有限公司 一种污水管道中泥沙平均淤堵高度的确定方法
CN111350162A (zh) * 2020-03-03 2020-06-30 河海大学 一种高滩深槽水库坝前恢复库容的挖槽设计方法
CN113159994A (zh) * 2021-04-25 2021-07-23 黄河勘测规划设计研究院有限公司 多泥沙河流的船闸淤积处理方法及装置
CN113297753A (zh) * 2021-06-24 2021-08-24 四川大学 山区河流透过性坝沿程水深预测方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
排沙洞进口泥沙冲淤三维数学模型分析计算;王新年;;水利水电技术(第05期);全文 *

Also Published As

Publication number Publication date
CN113821862A (zh) 2021-12-21

Similar Documents

Publication Publication Date Title
Liu et al. Influences of slope gradient on soil erosion
CN113821862B (zh) 一种拦河建筑物坝前泥沙淤积体平均坡角的确定方法
Mahdian Khalili et al. Time evolution effect on the scour characteristics downstream of the sluice gate with the submerged hydraulic jump in a laboratory model
Ghomeshi Reservoir sedimentation modelling
CN108827871B (zh) 一种管式泥沙侵蚀试验装置中泥沙表面切应力确定方法
Chavan et al. Alluvial channel hydrodynamics around tandem piers with downward seepage
Beg et al. Performance analysis of rectangular SIT (sediment invert trap) for stormwater drainage system
Kheder et al. Sediment transport investigations using three dimensional numerical modeling in a large Canal: Marala Ravi Link Canal (Pakistan)
De Villiers et al. Modelling of long-term sedimentation at Welbedacht Reservoir, South Africa
Bakri et al. Analysis of sediment distribution at the intake structure
Omran et al. A model for removing sediments from open channels
Azmeri et al. Evaluation of the Settling Basin of Keumala Weir, Krueng Baro River, Aceh, Indonesia
Yang et al. A graded sediment transport and bed evolution model in estuarine basins and its application to the Yellow River Delta
Ji et al. Application of SOBEK model in the Yellow River estuary
CN109934506B (zh) 一种基于定量化判断航道整治工程类型的选取方法
Zhang et al. Numerical Simulation Study of Multiphase Subsea Pipeline Sand Deposition Law
CN116876600B (zh) 带出流口调节系统的耙吸挖泥船泥舱消能箱及其施工工艺
Laksitaningtyas et al. Preliminary Experiment on Reservoir Sediment Flushing
Altaee et al. Time Variations of Sediment Floc Size and Density by using Settling Column Data
Dorfmann et al. Simulation of hydrodynamic and sediment transport processes–Two Austrian case studies
He et al. Coupling Analysis of Connectivity between Coarse Sediment Source Areas and Deposition Area in the Yellow River, China.
Havrevoll et al. Efficiency of Pressurized Rock Traps for Unlined Hydropower Tunnels. Energies 2021, 14, 4344
Bazarov et al. Experimental evaluation of the throughput of the supply channel of pumping stations
Bouisfi et al. Geometry effect of irrigation storage basin on particles removal efficiency: A computational fluid dynamics study
Xing et al. Field Experimental Study on Concentration Distribution Characteristics of Medium-Coarse Sand in Long-Distance Pipelines

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