CN106777558B - 一种泵站引渠及前池水体泥沙参考浓度的确定方法 - Google Patents
一种泵站引渠及前池水体泥沙参考浓度的确定方法 Download PDFInfo
- Publication number
- CN106777558B CN106777558B CN201611071246.3A CN201611071246A CN106777558B CN 106777558 B CN106777558 B CN 106777558B CN 201611071246 A CN201611071246 A CN 201611071246A CN 106777558 B CN106777558 B CN 106777558B
- Authority
- CN
- China
- Prior art keywords
- silt
- water body
- formula
- pumping plant
- approach channel
- 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
Links
Classifications
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Separation Of Suspended Particles By Flocculating Agents (AREA)
Abstract
本发明涉及一种泵站引渠及前池水体泥沙参考浓度的确定方法,包括以下步骤:1)计算得出泵站引渠及前池水体的流速随时间变化的关系式;2)计算泵站引渠及前池水体的最大谢尔兹数Θm、摩擦因子f和粗糙高度kN;3)计算泵站引渠及前池水体的泥沙沉降速度和泥沙在输沙层内的沉降时间与流动周期的比值Ψ;4)计算泵站引渠及前池水体的泥沙对流速的相位漂移ψ和泥沙相位的残留α;5)计算泵站引渠及前池水体的泥沙体积分数参考浓度。
Description
技术领域
本发明涉及一种泵站引渠及前池水体泥沙参考浓度的确定方法,属于水利工程泵站领域。
背景技术
在泥沙设计的数值模拟中,很关键的一步是给出泥沙浓度方程的底部边界条件,一般都在距离壁面2倍泥沙粒径高度的位置提出参考浓度C。边界条件的精度决定了数值模拟的结果精度。
精准知晓包含相位差作用和考虑泥沙体积守恒的泥沙参考浓度的大小,是准确模拟非恒定流动条件下水体挟沙率和输沙率、判断水泵机组磨蚀的重要基础。现有工程中广泛使用的参考浓度都是基于恒定流动的理论推导出来的,不包含相位差作用即泥沙运动相对于水动力条件所需要的时间响应,也未能考虑泥沙体积的守恒。这导致出现的问题是:无法体现泥沙运动相对流速的相位漂移和运动的泥沙总量在流动速度减少后的相位残留,给计算带来极大的误差;泥沙浓度随着流动条件的增大也难以达到饱和状态,参考浓度随流动强度太大从而失效。
发明内容
针对上述问题,本发明的目的是提供一种在非恒定流动条件下确定相关相位差参数和泥沙体积守恒定律来获得泵站引渠及前池水体泥沙参考浓度的方法。
为实现上述目的,本发明采用以下技术方案:一种泵站引渠及前池水体泥沙参考浓度的确定方法,包括以下步骤:
1)计算得出泵站引渠及前池水体的流速随时间变化的关系式U(t);
2)计算泵站引渠及前池水体的最大谢尔兹数Θm、摩擦因子f和粗糙高度kN;
3)计算泵站引渠及前池水体的泥沙沉降速度和泥沙在输沙层内的沉降时间与流动周期的比值Ψ;
4)计算泵站引渠及前池水体的泥沙对流速的相位漂移ψ和泥沙相位的残留α;
5)计算泵站引渠及前池水体的泥沙体积分数参考浓度。
所述步骤1)中,这里采用2阶Stokes波动理论得出泵站引渠及前池水体的流速随时间变化的关系式。
所述步骤2)中,采用如下的式(1)~(3)计算泵站引渠及前池泥沙的最大谢尔兹数Θm、摩擦因子f和粗糙高度kN,具体过程如下:
其中U是泵站引渠及前池水体的流速,下标m代表最大值;D是泥沙粒径;T是流动周期;Θ是谢尔兹数,下标m代表最大值;f是摩擦因子;s是泥沙与水密度的比值;所在地的重力加速度g;kN是动床面粗糙高度;
先假定最大谢尔兹数Θm<1,由式(3)得到粗糙高度的值kN,代入式(2),得到摩擦因子的值f;摩擦因子的值代入式(1)得到的最大谢尔兹数的值Θm若大于1则假定不成立;
根据最大谢尔兹数Θm>1,由式(3)得到粗糙高度kN=5ΘmD代入式(2),并联合式(1)得到
从式(4)解得摩擦因子的值f,代入式(1)得最大谢尔兹数的值Θm,代入式(3)得粗糙高度的值kN。
所述步骤3)中,泵站引渠及前池水体的泥沙沉降速度和泥沙在输沙层内的沉降时间与流动周期的比值Ψ,具体计算过程如下:
从如下的式(5)得到泥沙沉降速度的值,代入如下的式(6)得到泥沙在输沙层内的沉降时间与流动周期的比值Ψ;
其中ν是水的运动粘性系数。
所述步骤4)中,泵站引渠及前池水体的泥沙对流速的相位漂移ψ和泥沙相位的残留α,具体过程如下:
采用如下提供的式(7)和式(8),计算泥沙对流速的相位漂移ψ和泥沙相位的残留α
α=exp(-0.2/Ψ)。 (8)
所述步骤5)中,泵站引渠及前池泥沙的泥沙体积分数参考浓度的获取方法如下:
将所述步骤3)和所述步骤4)获得的泥沙对流速的相位漂移ψ、泥沙相位的残留α和泥沙沉降速度w代入下述式(9)得到泥沙体积分数参考浓度随时间变化的关系式;将式(1)得到的谢尔兹数代入式(9),即可得到最终的泥沙参考浓度随时间变化和相应的函数图像,式(9)如下所示:
其中F代表水体的挟沙能力。
本发明由于采取以上技术方案,其具有以下优点:1、通过泥沙体积守恒的考虑,避免参考浓度随流动强度太大而失效。2、通过所确定的相关相位差参数进行设计计算,可以反应非恒定流动过程中泥沙相对水动力条件所需要的时间响应。3、本发明通过所确定的相关相位差参数进行设计计算,体现了泥沙运动相对流速的相位漂移。4、本发明通过所确定的相关相位差参数进行设计计算,体现了具有运动的泥沙总量在流动速度减少后的相位残留。
附图说明
图1本发明工作过程采用的非恒定流动速度过程,U/Um是无量纲的流速,t/T是无量纲的时间;
图2是本发明与恒定流理论确定的参考浓度比较,流动强度Um=0.6m/s,C是泥沙参考浓度,实线代表本发明,点划线代表传统方法;
图3是本发明与恒定流理论确定的参考浓度比较,流动强度Um=1.0m/s,实线代表本发明,点划线代表传统方法;
图4是本发明最大泥沙体积参考浓度与流动强度的比较,Θm是最大谢尔兹数,Cm是最大泥沙体积参考浓度。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明提供的一种泵站引渠及前池水体泥沙参考浓度的确定方法,包括以下步骤:
1)计算得出泵站引渠及前池水体的流速随时间变化的关系式U(t),具体可以利用2阶Stokes波动理论得出泵站引渠及前池水体流速随时间变化的关系式,因该部分属于本领域的公知方法,因此不进行具体说明。
2)计算泵站引渠及前池水体的最大谢尔兹数Θm、摩擦因子f和粗糙高度kN,具体可通过下述的式(1)~(3)进行计算:
式中,U是泵站引渠及前池水体的流速,下标m代表最大值;D是泥沙粒径;T是流动周期;Θ是谢尔兹数,下标m代表最大值;f是摩擦因子;s是泥沙与水密度的比值;g是重力加速度;kN是动床面粗糙高度。
在具体求解时,先假定最大谢尔兹数Θm<1,由式(3)得到粗糙高度的值kN,代入式(2),得到摩擦因子的值f;摩擦因子的值代入式(1)得到的最大谢尔兹数的值Θm若大于1则假定不成立。
根据最大谢尔兹数Θm>1,由式(3)得到粗糙高度kN=5ΘmD代入式(2),并联合式(1)得到
从式(4)解得摩擦因子的值f,代入式(1)得最大谢尔兹数的值Θm,代入式(3)得粗糙高度的值kN。
3)计算泵站引渠及前池水体的泥沙沉降速度和泥沙在输沙层内的沉降时间与流动周期的比值Ψ,具体计算过程如下:
由下述式(5)计算得到泵站前池水体的泥沙沉降速度w,将w的值代入式下述式(6)得泥沙在输沙层内的沉降时间与流动周期的比值Ψ,运用的式(5)和式(6)如下:
式中,ν是水的运动粘性系数。
4)计算泵站引渠及前池水体的泥沙对流速的相位漂移ψ和泥沙相位的残留α,具体过程如下:
利用下述式(7)得到泵站前池的泥沙对流速的相位漂移ψ,由下述式(8)得泥沙相位的残留α的值:
α=exp(-0.2/Ψ)。 (8)
5)计算泵站引渠及前池水体的泥沙体积分数参考浓度C(t)。
将谢尔兹数Θ、泥沙对流速的相位漂移ψ、泥沙相位的残留α和泥沙沉降速度w代入下述式(9)得到泥沙体积分数参考浓度随时间变化的关系式:
其中F代表水体的挟沙能力。
下面结合附图和实例对本发明进行详细的描述。以下实例用于说明本发明,但不用来限制本发明的范围。
1)计算得出泵站引渠及前池水体的流速随时间变化的关系式。
已知泥沙所处气压为一个标准大气压、水温为20℃、最大流动速度Um=0.6m/s,流动周期T=5s,泥沙粒径D=1.0×10-4m,泥沙与水密度的比值s=2.65,所在地的重力加速度g=9.8m/s2,水的运动粘性系数ν=1.0×10-6m2/s,临界谢尔兹数Θcr取0.05。采用2阶Stoke波动理论得出流速与时间的关系式为:
U(t)=0.48cos[2π(t/T-0.214)]+0.12cos[4π(t/T-0.214)] (10)
其函数图像见图1。
2)计算最大谢尔兹数Θm、摩擦因子f和粗糙高度kN。
计算最大谢尔兹数Θm、摩擦因子f、粗糙高度kN,采用本发明提供的式(1)~(3)。
先假定最大谢尔兹数Θm<1,由式(3)得到粗糙高度kN=5.0×10-4m;代入式(2),得到摩擦因子f=8.0×10-3;代入式(1)得到最大谢尔兹数Θm=0.89<1。
3)计算泥沙沉降速度和泥沙在输沙层内的沉降时间与流动周期的比值Ψ。
采用本发明提供的式(6)和式(5),计算泥沙沉降速度w和泥沙在输沙层内的沉降时间与流动周期的比值Ψ。
从式(5)得到泥沙沉降速度w=8.4×10-3m/s,代入式(6)得到泥沙在输沙层内的沉降时间与流动周期的比值Ψ=1.5×10-1。
4)计算泥沙对流速的相位漂移ψ和泥沙相位的残留α。
采用本发明提供的式(7)和式(8),计算泥沙对流速的相位漂移ψ和泥沙相位的残留α。
根据泥沙在输沙层内的沉降时间与流动周期的比值Ψ,从式(7)得泥沙对流速的相位漂移ψ=1.2×10-1s,从式(8)得泥沙相位的残留α=0.27。
5)计算泥沙体积分数参考浓度,采用本发明提供的式(9)和式(1)。
计算泥沙体积分数参考浓度C,采用本发明提供的式(9)和式(1)。
将步骤3)和步骤4)获得的泥沙对流速的相位漂移ψ、泥沙相位的残留α和泥沙沉降速度w代入式(9)得到
将式(1)得到的谢尔兹数Θ代入式(11)即可得到最终的泥沙参考浓度随时间变化的关系式。如图2是此种情况下泥沙参考浓度与恒定流理论的泥沙浓度的函数图像。
下面将最大流动速度Um增加到1.0m/s,重复步骤2)~步骤5)的过程,得到谢尔兹数Θm=2.86>1、摩擦因子f=9.3×10-3、粗糙高度kN=1.43×10-3m、泥沙对流速的相位漂移ψ=5.0×10-1s、泥沙相位的残留α=0.63。图3是此种情况下泥沙参考浓度与恒定流理论的泥沙浓度的函数图像。
从上述具体实施例中得出的说明书附图我们还可以更直观的体现本发明的以下特征。
1.体现了泥沙运动相对流速的相位漂移。
图2本发明的实线C落后于图1的流速t/T=0.024相位。如图1中U/Um=0的t/T=0时刻、最大值U/Um=1的t/T=0.21时刻;由于相位漂移,对应图2中实线C最小值的是t/T=0.024时刻、最大值的是t/T=0.234时刻。由于流动强度增大,图3本发明的实线C落后于图1的流速相位幅度有所增大。而代表恒定流理论的点划线则与图1的流速同相位。
2.具有运动的泥沙总量在流动速度减少后的相位残留。
非恒定流动中,已悬浮泥沙的沉降是需要时间的。相位残留的存在使得泥沙体积浓度一般不会随着瞬时的流速减小而减为0。图2、图3代表本发明的实线最小值分别为0.05和0.19,而代表恒定流理论的点划线最小值为0。从图3到图2,相位差作用减小,本发明的结果(实线)越来越趋向于恒定流动的结果(虚线)。
3.避免了参考浓度随流动强度增大而变无穷大。
如附图4所示,由于整体体积守恒的考虑,泥沙体积参考浓度最大值Cm不会随着流动强度Θm的增加而无限增大,而是趋向于一个饱和状态。
上述各实施例仅用于对本发明的目的、技术方案和有益效果进行了进一步详细说明,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (4)
1.一种泵站引渠及前池水体泥沙参考浓度的确定方法,包括以下步骤:
1)计算得出泵站引渠及前池水体的流速随时间变化的关系式U(t);
2)计算泵站引渠及前池水体的最大谢尔兹数Θm、摩擦因子f和动床面粗糙高度kN;
3)计算泵站引渠及前池水体的泥沙沉降速度w,以及泥沙在输沙层内的沉降时间与流动周期T的比值Ψ;
4)计算泵站引渠及前池水体的泥沙对流速的相位漂移ψ和泥沙相位的残留α,具体计算过程如下:
采用如下提供的式(7)和式(8),计算泥沙对流速的相位漂移ψ和泥沙相位的残留α
α=exp(-0.2/Ψ) (8)
5)计算泵站引渠及前池水体的泥沙体积分数参考浓度C(t),具体计算过程如下:
将所述步骤3)和所述步骤4)获得的泥沙对流速的相位漂移ψ、泥沙相位的残留α和泥沙沉降速度w代入下述式(9)得到泥沙体积分数参考浓度随时间变化的关系式;将式(1)得到的谢尔兹数代入式(9),即可得到最终的泥沙参考浓度随时间变化和相应的函数图像,式(9)如下所示:
其中F代表水体的挟沙能力;w是泥沙沉降速度;t是时间;U是泵站引渠及前池水体流速,下标m代表最大值;Θcr是临界谢尔兹数。
2.如权利要求1所述的一种泵站引渠及前池水体泥沙参考浓度的确定方法,其特征在于:所述步骤1)中,这里采用2阶Stokes波动理论得出泵站引渠及前池水体的流速随时间变化的关系式。
3.如权利要求2所述的一种泵站引渠及前池水体泥沙参考浓度的确定方法,其特征在于:所述步骤2)中,采用如下的式(1)~(3)计算泵站引渠及前池水体的最大谢尔兹数Θm、摩擦因子f和动床面粗糙高度kN,具体过程如下:
其中U是泵站引渠及前池水体的流速,下标m代表最大值;D是泥沙粒径;T是流动周期;Θ是谢尔兹数,下标m代表最大值;f是摩擦因子;s是泥沙密度与水密度的比值;所在地的重力加速度g;kN是动床面粗糙高度;
先假定最大谢尔兹数Θm<1,由式(3)得到动床面粗糙高度kN的值,代入式(2),得到摩擦因子f的值;摩擦因子f的值代入式(1),得到的最大谢尔兹数Θm的值若大于1,则假定不成立;
根据最大谢尔兹数Θm>1,由式(3)得到动床面粗糙高度kN=5ΘmD,代入式(2),并联合式(1)得到
从式(4)解得摩擦因子f的值,代入式(1)得最大谢尔兹数Θm的值,代入式(3)得动床面粗糙高度kN的值。
4.如权利要求3所述的一种泵站引渠及前池水体泥沙参考浓度的确定方法,其特征在于:所述步骤3)中,泵站引渠及前池水体的泥沙沉降速度,以及泥沙在输沙层内的沉降时间与流动周期T的比值Ψ,具体计算过程如下:
从如下的式(5)得到泥沙沉降速度的值,代入如下的式(6)得到泥沙在输沙层内的沉降时间与流动周期的比值Ψ;
其中ν是水的运动粘性系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611071246.3A CN106777558B (zh) | 2016-11-29 | 2016-11-29 | 一种泵站引渠及前池水体泥沙参考浓度的确定方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611071246.3A CN106777558B (zh) | 2016-11-29 | 2016-11-29 | 一种泵站引渠及前池水体泥沙参考浓度的确定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106777558A CN106777558A (zh) | 2017-05-31 |
CN106777558B true CN106777558B (zh) | 2019-08-09 |
Family
ID=58905330
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611071246.3A Active CN106777558B (zh) | 2016-11-29 | 2016-11-29 | 一种泵站引渠及前池水体泥沙参考浓度的确定方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106777558B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111259606B (zh) * | 2018-11-14 | 2021-12-14 | 中国农业大学 | 一种泵站引渠及前池含沙动床面流速分布确定方法和装置 |
CN111291459B (zh) * | 2018-11-21 | 2022-02-08 | 中国农业大学 | 一种泵站引渠及前池泥沙通量的确定方法及*** |
CN111208274B (zh) * | 2018-11-21 | 2022-12-09 | 中国农业大学 | 一种泵站引渠及前池水体泥沙浓度分布确定方法及*** |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7159841B2 (en) * | 2002-11-07 | 2007-01-09 | The United States Of America As Represented By The United States Department Of Energy | Piezoelectric axial flow microvalve |
CN101265914A (zh) * | 2008-05-12 | 2008-09-17 | 杭州小尔机电设计有限公司 | 两相流潜水泵疏浚***中加压站最优数量的确定方法 |
CN102032935A (zh) * | 2010-12-07 | 2011-04-27 | 杭州电子科技大学 | 城市排水汇流管网污水泵站流量的软测量方法 |
US9139791B2 (en) * | 2008-02-13 | 2015-09-22 | Hydrocoal Technologies, Llc | Processing device for improved utilization of fuel solids |
WO2016024882A1 (ru) * | 2014-08-12 | 2016-02-18 | Ооо "Фармасапфир" | СПОСОБ ОБРАБОТКИ ТРУЩИХСЯ ПОВЕРХНОСТЕЙ ДЕТАЛЕЙ ИЗ МОНОКРИСТАЛЛА НА ОСНОВЕ α-Аl2Оз |
CN105613220A (zh) * | 2015-12-30 | 2016-06-01 | 中国农业科学院农业资源与农业区划研究所 | 定量灌溉控制*** |
-
2016
- 2016-11-29 CN CN201611071246.3A patent/CN106777558B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7159841B2 (en) * | 2002-11-07 | 2007-01-09 | The United States Of America As Represented By The United States Department Of Energy | Piezoelectric axial flow microvalve |
US9139791B2 (en) * | 2008-02-13 | 2015-09-22 | Hydrocoal Technologies, Llc | Processing device for improved utilization of fuel solids |
CN101265914A (zh) * | 2008-05-12 | 2008-09-17 | 杭州小尔机电设计有限公司 | 两相流潜水泵疏浚***中加压站最优数量的确定方法 |
CN102032935A (zh) * | 2010-12-07 | 2011-04-27 | 杭州电子科技大学 | 城市排水汇流管网污水泵站流量的软测量方法 |
WO2016024882A1 (ru) * | 2014-08-12 | 2016-02-18 | Ооо "Фармасапфир" | СПОСОБ ОБРАБОТКИ ТРУЩИХСЯ ПОВЕРХНОСТЕЙ ДЕТАЛЕЙ ИЗ МОНОКРИСТАЛЛА НА ОСНОВЕ α-Аl2Оз |
CN105613220A (zh) * | 2015-12-30 | 2016-06-01 | 中国农业科学院农业资源与农业区划研究所 | 定量灌溉控制*** |
Non-Patent Citations (3)
Title |
---|
《Effects of wave asymmetry on sheet flow》;Wilson;《Coastal Engineering》;19951231;第25卷(第3-4期);第191-204页 * |
《泵站前池水沙流的LBM—LES数值模拟研究》;唐学林 等;《现代节水高效农业与生态灌区建设(上)》;20101231;第799-810页 * |
《考虑颗粒动态尺度影响的泥沙扩散系数模型建立及应用》;张自超 等;《农业工程学报》;20161123;第32卷(第22期);第129-137页 * |
Also Published As
Publication number | Publication date |
---|---|
CN106777558A (zh) | 2017-05-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106777558B (zh) | 一种泵站引渠及前池水体泥沙参考浓度的确定方法 | |
CN106777557B (zh) | 一种泵站引渠及前池水体挟沙率的确定方法 | |
CN107341320B (zh) | 一种抽水蓄能电站的尾水调压室设置判别方法 | |
CN101275459B (zh) | 控制井中的流 | |
CN106351645B (zh) | 一种抽油机井动液面连续测量方法及装置 | |
CN110991760A (zh) | 一种高气液比产水气井临界携液流速预测方法 | |
CN104990450B (zh) | 基于dcs的高位集水冷却塔水位维持*** | |
CN106709153B (zh) | 一种泵站前池扬沙率的确定方法 | |
CN102128658B (zh) | 一种完全分离型油、气、水多相流量计 | |
Xu et al. | Flow pattern and anti-silt measures of straight-edge forebay in large pump stations. | |
CN105068582B (zh) | 基于dcs的高位集水冷却塔水位维持控制方法 | |
CN110941868A (zh) | 一种水电站调压室临界稳定断面的计算方法 | |
CN102619791B (zh) | 潮汐式虹吸真空抽气方法和装置 | |
CN104050341B (zh) | 基于液相与气相耦合的调压室水位波动模拟方法及*** | |
CN105547885A (zh) | 气井冲蚀实验***及其冲蚀方法和流束调节装置 | |
CN111046569B (zh) | 一种空化水射流喷嘴结构设计方法 | |
CN101469659B (zh) | 水机气蚀根治方法 | |
CN204729282U (zh) | 隔膜泵试验装置 | |
CN108343401A (zh) | 抽油机调衡节能装置 | |
ElGammal et al. | Numerical investigations on intake tube design of micro kaplan hydro-turbine system | |
CN111967143B (zh) | 一种适于深水气井的临界携液流量预测方法 | |
Kramer et al. | Minimising the air demand of micro-hydro impulse turbines in counter pressure operation | |
Oh et al. | OWC design to increase wave energy absorption efficiency in wave conversion systems | |
Tang et al. | Effect of air injector on the airlift performance in air-water-solid three-phase flow | |
Tan et al. | Numerical study on cylindrical particle solid–liquid two-phase flow in a slurry pump |
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 |