CN107256312B - 一种基于潮流环境下海湾纳潮变化量计算方法 - Google Patents
一种基于潮流环境下海湾纳潮变化量计算方法 Download PDFInfo
- Publication number
- CN107256312B CN107256312B CN201710443352.8A CN201710443352A CN107256312B CN 107256312 B CN107256312 B CN 107256312B CN 201710443352 A CN201710443352 A CN 201710443352A CN 107256312 B CN107256312 B CN 107256312B
- Authority
- CN
- China
- Prior art keywords
- bay
- river
- reclamation
- tidal
- water
- 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.)
- Expired - Fee Related
Links
Images
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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F22/00—Methods or apparatus for measuring volume of fluids or fluent solid material, not otherwise provided for
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Fluid Mechanics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Other Liquid Machine Or Engine Such As Wave Power Use (AREA)
Abstract
本发明公开了一种基于潮流环境下海湾纳潮变化量计算方法,涉及海湾围垦工程。本发明其包括海湾围垦工程实际占用的海水体积的变化量;再以海湾围垦工程实施后造成的海湾水体潮波变形所导致的海水体积变化量来修正海湾纳潮变化量的误差。本发明给出基于数值模拟计算方法的新方法,为数值模拟计算方法与传统计算方法的改进方法。它真实有效的揭示了围垦工程实施造成的纳潮变化量。在实践中并证明了相关理论,为整个海湾进行围垦工程的设计与研究奠定了基础。
Description
技术领域
本发明涉及海湾围垦工程,具体涉及一种基于潮流环境下海湾纳潮变化量计算方法。
背景技术
近年来,随着我国改革开放的不断推进,经济建设的热潮一浪高过一浪,工程建设行业也得到了前所未有的大发展。沿海地区,尤其是沿海港口城市,土地资源日超紧缺,海涂围垦工程蓬勃发展。然而在河口地区,尤其是海湾内进行围垦会造成海湾纳潮量的改变。
纳潮量是海湾可以接纳的潮水体积,是一个潮周期中的最大水交换量,为海湾高潮水量与低潮水量之差。其大小直接影响着海湾水深的维持以及海湾与外海的水体交换和污染物的迁移扩散,制约着海湾的自净能力和环境容量,对维护海湾良好的生态环境至关重要。因此,研究围垦工程实施后,海湾内纳潮量的变化至关重要。
现有纳潮量的计算方法主要有两种:
(1)传统计算方法。即纳潮量计算可通过海湾高潮时的最大纳水体积与港湾低潮时的最小纳水体积之差求出,公式为:
W=0.5×(S1+S2)(h1-h2)
式中:W为纳潮量;h1、h2分别为高、低潮位;S1、S2分别为高、低潮潮位对应的水域面积。
(2)ADCP测量方法。该方法通过在半封闭海湾的湾口设计一个封闭断面直接走航测量得到进出该断面的水量(纳潮量)。
(3)数值模拟方法。通过建立并验证工程海域潮流数学模型,计算出封闭湾口断面的潮位及流速流向,得到通过湾口断面的流量变化序列,再将这些流量值积分可得到海湾的纳潮量。
对于传统计算方法,该方法将高低潮位之间水体默认为坡度均匀一致的斜面,存在一定的局限性(如图1所示)。同时,在计算围垦工程实施后海湾纳潮变化量时,无法考虑围垦工程造成的潮波变形情况,且当围垦工程顶高程低于高潮位或围垦工程位置水深高于低潮位时则无法计算围垦工程对纳潮量变化的影响。
对于ADCP测量方法。该方法需要进行现场实测工作,除费用高之外,因海湾繁忙的船行以及较差的海况经常难于实施,即便实施测得的结果也难以做到同步而失真,需要做一定的数据修正工作。同时,ADCP测量方法只能测量现状条件下海湾纳潮量,并不能预测围垦工程实施后造成的海湾纳潮量变化。
对于数值模拟计算方法,其计算误差直接取决于模型研究人员的水平。
发明内容
本发明所要解决的技术问题是提供一种基于潮流环境下海湾纳潮变化量的计算方法,旨在给出基于数值模拟计算方法的新方法,为数值模拟计算方法与传统计算方法的改进方法。它真实有效的揭示了围垦工程实施造成的纳潮变化量。
为解决上述技术问题,本发明所采取的技术方案是:一种基于潮流环境下海湾纳潮变化量计算方法,其特征在于:其包括海湾围垦工程实际占用的海水体积的变化量;再以海湾围垦工程实施后造成的海湾水体潮波变形所导致的海水体积变化量来修正海湾纳潮变化量的误差。
进一步的技术方案在于,对海湾围垦工程实施后造成的海湾水体潮波变形所导致的水体体积变化量,是基于建立海湾潮流运动数学模型进行数值模拟,计算出封闭湾口断面以内水体高低潮位变化,进而得到水体体积变化量。
进一步的技术方案在于,所述海湾潮流运动数学模型为:
控制方程如下:
其中:h=η+d,η和d分别表示水面高度和静水深,x和y分别表示横轴和纵轴坐标,t为时间,g为重力加速度,和分别为沿x和y方向的深度平均流速,f为柯氏力系数,ρ为流体密度,ρ0为参考密度,S为点源流量,us与vs为点源流速,Tij为应力项,包括粘性应力、紊流应力和对流等,根据水深平均的流速梯度计算;
Manning系数可以根据底部糙率计算;
进一步的技术方案在于,所述海湾纳潮变化量为:
ΔW=S工(h工1-h工2)+S河1×Δh河1+S河2×Δh河2
式中:ΔW为海湾纳潮变化量;S工为围填海工程面积;h工1为围填海工程附近高潮位(当高潮位低于围垦工程堤顶高程)或围垦工程堤顶高程(当高潮位高于围垦工程堤顶高程);h工2为围填海工程附近低潮位(当低潮位高于围垦工程附近平均水深)或围垦工程附近水深(当低潮位低于围垦工程附近平均水深);S河1、S河2分别为高、低潮潮位对应的河口水域面积(除去围垦工程面积);Δh河1、Δh河2分别数值模拟计算出的高、低潮潮位变化值(海湾平均);Δh河1增加为正值,减小为负值;Δh河2增加为负值,减小为正值。
采用上述技术方案所产生的有益效果在于:本发明给出基于数值模拟计算方法的新方法,为数值模拟计算方法与传统计算方法的改进方法。它真实有效的揭示了围垦工程实施造成的纳潮变化量。在实践中并证明了相关理论,为整个海湾进行围垦工程的设计与研究奠定了基础。
附图说明
下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1是传统计算方法断面示意图;
图2是围垦工程前示意图(图中灰色为海域,白色为陆域或岛屿);
图3是围垦工程后示意图(图中灰色为海域,白色为陆域或岛屿);
图4是围垦工程占用海湾纳潮水体体积示意图(图中灰色水体);
图5是实测潮位及流速流向测站示意图;
图6是网格示意图;
图7是潮位验证图;
图8是流速流向验证图;
图9是表1的特征点布置图。
具体实施方式
下面结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
本发明提供了一种基于潮流环境下海湾纳潮变化量计算方法,其特征在于:其包括海湾围垦工程实际占用的海水体积的变化量;再以海湾围垦工程实施后造成的海湾水体潮波变形所导致的海水体积变化量来修正海湾纳潮变化量的误差。
优选的,对海湾围垦工程实施后造成的海湾水体潮波变形所导致的水体体积变化量,是基于建立海湾潮流运动数学模型进行数值模拟,计算出封闭湾口断面以内水体高低潮位变化,进而得到水体体积变化量。
优选的,所述海湾潮流运动数学模型为:
控制方程如下:
其中:h=η+d,η和d分别表示水面高度和静水深,x和y分别表示横轴和纵轴坐标,t为时间,g为重力加速度,和分别为沿x和y方向的深度平均流速,f为柯氏力系数,ρ为流体密度,ρ0为参考密度,S为点源流量,us与vs为点源流速,Tij为应力项,包括粘性应力、紊流应力和对流等,根据水深平均的流速梯度计算;
Manning系数可以根据底部糙率计算;
优选的,所述海湾纳潮变化量为:
ΔW=S工(h工1-h工2)+S河1×Δh河1+S河2×Δh河2
式中:ΔW为海湾纳潮变化量;S工为围填海工程面积;h工1为围填海工程附近高潮位(当高潮位低于围垦工程堤顶高程)或围垦工程堤顶高程(当高潮位高于围垦工程堤顶高程);h工2为围填海工程附近低潮位(当低潮位高于围垦工程附近平均水深)或围垦工程附近水深(当低潮位低于围垦工程附近平均水深);S河1、S河2分别为高、低潮潮位对应的河口水域面积(除去围垦工程面积);Δh河1、Δh河2分别数值模拟计算出的高、低潮潮位变化值(海湾平均);Δh河1增加为正值,减小为负值;Δh河2增加为负值,减小为正值。
实施例:
一、纳潮变化量分解
纳潮量为海湾高潮时最大纳水体积与海湾低潮时最小纳水体积之差。以某海湾内围垦工程为例(图2和图3),整个海湾面积可分为两部分,即:
S=S工+S河
其纳潮量也可分为两部分:
W=W工+W河
则纳潮变化量也可分为两部分:
ΔW=ΔW工+ΔW河
二、围垦工程占用海湾水体体积
根据上述公式可见,整个海湾进行围垦工程后,海湾纳潮变化量包括两部分,其中一部分为围垦工程实际占用的海水体积,一部分为围垦工程实施后造成的海湾水体潮波变形所导致的水体体积变化量。其中,ΔW工可采用实测资料进行计算,如图4所示。
ΔW工=S工(h工1-h工2)
其中:S工为围填海工程面积(已知);h工1为围填海工程附近高潮位(当高潮位低于围垦工程堤顶高程时),或围垦工程堤顶高程(当高潮位高于围垦工程堤顶高程时);h工2为围填海工程附近低潮位(当低潮位高于围垦工程附近平均水深时)或围垦工程附近水深(当低潮位低于围垦工程附近平均水深时)。
其中,围垦工程附近水深可根据实测水深得到。高、低潮位可根据工程附近潮位站实测潮位资料得到,对应到本工程即为工程附近T4潮位站实测潮位资料(如图5所示)。
本工程围填海面积约为2.95km2,实测大潮两次涨潮过程潮差分别为5.20m和5.66m(工程附近T4站)。则围垦工程所占的潮棱体约0.32×108m3。
三、潮波变形造成纳潮变化量
ΔW河为围垦工程实施后海湾水域潮波变形导致的纳潮变化量,可采用数值模拟方法进行计算。可采用自主开发程序或比较成熟的MIKE21、FVCOM等软件计算。
3.1潮流运动数学模型
控制方程如下:
其中:h=η+d,η和d分别表示水面高度和静水深,x和y分别表示横轴和纵轴坐标,t为时间,g为重力加速度,和分别为沿x和y方向的深度平均流速,f为柯氏力系数,ρ为流体密度,ρ0为参考密度,S为点源流量,us与vs为点源流速,Tij为应力项,包括粘性应力、紊流应力和对流等,根据水深平均的流速梯度计算。
Manning系数可以根据底部糙率计算。
3.2计算域选取及网格剖分
以上述某海区为例进行计算:计算区域西边界取至上游水闸,东西方向长约97km,北方向长约62km。模型模拟区域总面积达6000km2。
如图6所示,采用三角形网格进行剖分,距离工程区较远的海区网格较稀疏,对工程海域所在的海湾进行网格加密。计算域采用的三角形网格边长最大为2000m,最小网格边长为10m,总共有13669个网格节点,25028个三角形单元。
3.3边界条件及概化方法
外海开边界采用潮位过程控制,由中国海域潮汐模型提供。上游主要径流根据水文站实测流量给出。模型中可考虑漫滩、潜堤等动边界的处理。
3.4模型验证
为验证模型的合理性,采用实测水文全潮资料对模型进行率定和验证(测站参见图5)。部分潮位及流速流向验证曲线参见图7和图8,精度满足规程要求。
3.5潮位变化统计
在工程所在海湾内均匀地选取尽量多的特征点(图9)。表1给出了数值模拟计算出来的围填海工程实施前后特征点位置高低潮位变化情况
表1特征点位置高低潮位变化值(m)
3.6潮波变形造成纳潮变化量
潮波变形造成的纳潮变化量为高潮水体变化量与低潮水体变化量之和:
ΔW河=S河1×Δh河1+S河2×Δh河2
S河1、S河2分别为高、低潮潮位对应的河口水域面积(除去围垦工程),其可根据海湾实测水深及工程附近实测高低潮位统计得到。
Δh河1、Δh河2分别为数值模拟计算出的高、低潮潮位变化值(海湾平均)。
其中:Δh河1增加为正值,减小为负值;Δh河2增加为负值,减小为正值。
本文实例计算可得,由于潮波变形造成的纳潮量变化约为0.043×108m3。
四、海湾纳潮变化量计算方法总结
针对围垦工程实施后海湾纳潮变化量,作者提出了一种新的计算方法,该方法基于潮流数值模拟,公式如下,为上述两部分之和:
ΔW=ΔW工+ΔW河
ΔW工=S工(h工1-h工2)
ΔW河=S河1×Δh河1+S河2×Δh河2
ΔW=S工(h工1-h工2)+S河1×Δh河1+S河2×Δh河2
式中:ΔW为由于围填海工程造成的海湾纳潮变化量;S工为围填海工程面积;h工1为围填海工程附近高潮位(当高潮位低于围垦工程堤顶高程)或围垦工程堤顶高程(当高潮位高于围垦工程堤顶高程);h工2为围填海工程附近低潮位(当低潮位高于围垦工程附近平均水深)或围垦工程附近水深(当低潮位低于围垦工程附近平均水深);S河1、S河2分别为高、低潮潮位对应的河口水域面积(除去围垦工程)。Δh河1、Δh河2分别数值模拟计算出的高、低潮潮位变化值(海湾平均)。Δh河1增加为正值,减小为负值;Δh河2增加为负值,减小为正值。
综上所述,本为实例中围填海工程实施后纳潮变化量为0.363×108m3。
Claims (2)
1.一种基于潮流环境下海湾纳潮变化量计算方法,其特征在于:其包括海湾围垦工程实际占用的海水体积的变化量;再以海湾围垦工程实施后造成的海湾水体潮波变形所导致的海水体积变化量来修正海湾纳潮变化量的误差;
对海湾围垦工程实施后造成的海湾水体潮波变形所导致的水体体积变化量,是基于建立海湾潮流运动数学模型进行数值模拟,计算出封闭湾口断面以内水体高低潮位变化,进而得到水体体积变化量;
所述海湾纳潮变化量为:
ΔW=S工(h工1-h工2)+S河1×Δh河1+S河2×Δh河2
式中:ΔW为海湾纳潮变化量;S工为围填海工程面积;h工1为当高潮位低于围垦工程堤顶高程时的围填海工程附近高潮位或当高潮位高于围垦工程堤顶高程时的围垦工程堤顶高程;h工2为当低潮位高于围垦工程附近平均水深时的围填海工程附近低潮位或当低潮位低于围垦工程附近平均水深时的围垦工程附近水深;S河1、S河2分别为除去围垦工程面积的高、低潮潮位对应的河口水域面积;Δh河1、Δh河2分别数值模拟计算出的海湾平均高、低潮潮位变化值;Δh河1增加为正值,减小为负值;Δh河2增加为负值,减小为正值。
2.根据权利要求1所述的一种基于潮流环境下海湾纳潮变化量计算方法,其特征在于:所述海湾潮流运动数学模型为:
控制方程如下:
其中:h=η+d,η和d分别表示水面高度和静水深,x和y分别表示横轴和纵轴坐标,t为时间,g为重力加速度,和分别为沿x和y方向的深度平均流速,f为柯氏力系数,ρ为流体密度,ρ0为参考密度,S为点源流量,us与vs为点源流速,Tij为应力项,包括粘性应力、紊流应力和对流,根据水深平均的流速梯度计算;
Manning系数根据底部糙率计算;风
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710443352.8A CN107256312B (zh) | 2017-06-13 | 2017-06-13 | 一种基于潮流环境下海湾纳潮变化量计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710443352.8A CN107256312B (zh) | 2017-06-13 | 2017-06-13 | 一种基于潮流环境下海湾纳潮变化量计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107256312A CN107256312A (zh) | 2017-10-17 |
CN107256312B true CN107256312B (zh) | 2020-09-25 |
Family
ID=60024043
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710443352.8A Expired - Fee Related CN107256312B (zh) | 2017-06-13 | 2017-06-13 | 一种基于潮流环境下海湾纳潮变化量计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107256312B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108595831B (zh) * | 2018-04-22 | 2021-05-07 | 天津大学 | 一种计算实时潮致混合的动态模拟方法 |
CN111695250B (zh) * | 2020-06-04 | 2022-12-13 | 哈尔滨工程大学 | 一种内潮波特征提取方法 |
CN115100537B (zh) * | 2022-06-27 | 2023-03-28 | 中国科学院地理科学与资源研究所 | 一种基于遥感影像的潮汐能资源评估方法 |
CN116796517B (zh) * | 2023-05-30 | 2024-02-02 | 山东省海洋资源与环境研究院(山东省海洋环境监测中心、山东省水产品质量检验中心) | 一种湾口断面纳潮量计算方法、系统、设备及介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103774605A (zh) * | 2014-01-10 | 2014-05-07 | 河海大学 | 一种提高环抱式港池水体交换能力的设计方法 |
CN106168991A (zh) * | 2016-06-24 | 2016-11-30 | 珠江水利委员会珠江水利科学研究院 | 一种基于水动力数值模拟的感潮河网潮位预报方法 |
CN106709262A (zh) * | 2017-01-11 | 2017-05-24 | 河海大学 | 基于动态概化的潮流能水轮机阵列水动力影响数值模拟方法 |
-
2017
- 2017-06-13 CN CN201710443352.8A patent/CN107256312B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103774605A (zh) * | 2014-01-10 | 2014-05-07 | 河海大学 | 一种提高环抱式港池水体交换能力的设计方法 |
CN106168991A (zh) * | 2016-06-24 | 2016-11-30 | 珠江水利委员会珠江水利科学研究院 | 一种基于水动力数值模拟的感潮河网潮位预报方法 |
CN106709262A (zh) * | 2017-01-11 | 2017-05-24 | 河海大学 | 基于动态概化的潮流能水轮机阵列水动力影响数值模拟方法 |
Non-Patent Citations (5)
Title |
---|
Combined effects of channel dredging,land reclamation and long-range jetties upon the long-term evolution of channel-shoal system in Qinzhou bay,SW China;Yu-Hai Wang;《Ocean Engineering》;20141115;全文 * |
半封闭海湾潮间带部分围垦后纳潮量计算的商榷——以胶州湾为例;杨世伦 等;《海洋科学》;20030831;全文 * |
张坤.钦州湾茅尾海纳潮量和水交换能力的数值模拟研究.《中国优秀硕士学位论文全文数据库 基础科学辑》.2015,A010-44. * |
钦州湾茅尾海纳潮量和水交换能力的数值模拟研究;张坤;《中国优秀硕士学位论文全文数据库 基础科学辑》;20150215;第1-3章 * |
风暴潮模拟中潮位对风拖曳力系数的影响研究;孔俊 等;《海洋预报》;20080229;正文第74-75页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107256312A (zh) | 2017-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107256312B (zh) | 一种基于潮流环境下海湾纳潮变化量计算方法 | |
CN108256137B (zh) | 一种强潮河口湾人工岛作业区港池航道回淤模拟方法 | |
CN109271672B (zh) | 一种河-湖-泵站相互影响作用下的河道水面线计算方法 | |
CN109582996A (zh) | 一种小尺度岸滩剖面与大尺度岸线变化的耦合模拟方法 | |
CN110135069B (zh) | 输水隧洞输水时的泥沙特征获取方法、装置、计算机设备 | |
CN115840975B (zh) | 一种风暴潮增水漫堤预警方法、系统、装置和存储介质 | |
CN108021780A (zh) | 一种基于无规则非结构网格模型的山洪动态仿真方法 | |
CN110955998A (zh) | 一种基于gis的大范围泥石流数值模拟及数值处理方法 | |
CN106759063A (zh) | 一种感潮河段造床流量的计算方法 | |
CN108625337B (zh) | 一种确定潮流界以下沙质河床河段整治水位的方法 | |
Wei et al. | Effects of flood on thermal structure of a stratified reservoir | |
CN110287504B (zh) | 一种模拟近岸海水交换规律的耦合模型的建立方法 | |
Jifu et al. | Influences of the waterway project in the Yangtze River Estuary on sediment transport | |
CN116703072A (zh) | 河库连通河网地区的水量-水质综合调度方法 | |
He et al. | Numerical Simulation of Hydrodynamic Conditions on Closure Gap of Immersed Tunnel in Hong Kong-Zhuhai-Macao Bridge | |
Sharbaty | Two Dimensional simulations of seasonal flow patterns in the Gorgan Bay | |
Wang et al. | Research and application of hydrodynamics modeling of channel in reservoir area-Case of Feilaixia station to Qingyuan station section | |
Huang et al. | Cumulative effect of the perennial reclamation projects on water exchange in Sanmen Bay | |
Zang et al. | Application of ROMS for simulating evolution and migration of tidal sand waves | |
CN117313425B (zh) | 一种年均含沙量的计算方法 | |
Mu et al. | Effects of the large-scale reclamation project on hydro-dynamic environment in the estuary | |
Pan et al. | Evolution of water exchanging period in Lvdao lagoon induced by different water depth | |
Miao et al. | Numerical Simulation of Hydrodynamic Fields Based on Unsteady Model in Tidal Estuary | |
CN118350319B (zh) | 河口最大浑浊带变动趋势预测方法、设备、介质及产品 | |
CN116894401A (zh) | 一种由数值模拟分析获得水库水体咸淡水界面变化规律的方法和应用 |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200925 |