CN114611297B - 一种基于淹水频率的河湖湿地植物修复区域优选方法 - Google Patents

一种基于淹水频率的河湖湿地植物修复区域优选方法 Download PDF

Info

Publication number
CN114611297B
CN114611297B CN202210245128.9A CN202210245128A CN114611297B CN 114611297 B CN114611297 B CN 114611297B CN 202210245128 A CN202210245128 A CN 202210245128A CN 114611297 B CN114611297 B CN 114611297B
Authority
CN
China
Prior art keywords
river
lake
wetland
grid
partition
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
CN202210245128.9A
Other languages
English (en)
Other versions
CN114611297A (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.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
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 China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN202210245128.9A priority Critical patent/CN114611297B/zh
Publication of CN114611297A publication Critical patent/CN114611297A/zh
Application granted granted Critical
Publication of CN114611297B publication Critical patent/CN114611297B/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/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • 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/10Numerical modelling
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A40/00Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
    • Y02A40/10Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in agriculture
    • Y02A40/22Improving land use; Improving water use or availability; Controlling erosion

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Fluid Mechanics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Purification Treatments By Anaerobic Or Anaerobic And Aerobic Bacteria Or Animals (AREA)
  • Cultivation Of Plants (AREA)

Abstract

本发明公开了一种基于淹水频率的河湖湿地植物修复区域优选方法,包括以下步骤:S1、根据河湖支流流量等数据,建立河湖水动力数学模型;S2、采用河湖水动力数学模型,计算各个网格的淹水频率;S3、根据各个网格的淹水频率,对河湖水域和河湖湿地植物修复区域分别进行分区,得到各个网格的河湖水文分区和湿地植物修复分区类型;S4、根据各个网格的河湖水文分区和湿地植物修复分区类型,绘制基于淹水频率的河湖湿地植物修复分区图;S5、根据河湖湿地植物修复分区图,在不同修复分区种植相应的河湖湿地植物;本发明解决了现有缺乏一种对河湖湿地植物修复区域进行定量化优选的方法的问题。

Description

一种基于淹水频率的河湖湿地植物修复区域优选方法
技术领域
本发明涉及河湖湿地植物修复技术领域,具体涉及一种基于淹水频率的河湖湿地植物修复区域优选方法。
背景技术
湿地植物是河湖水生态***的重要组成部分,湿地植物修复重建是河湖水生态***保护的关键措施。湿地植物生长繁殖受到水文情势的显著影响,尤其是对于鄱阳湖、洞庭湖等水文情势剧烈波动的大型湖泊,其湿地植物往往表现为沿水文情势梯度的条带状分布特征。不同类型的湿地植物有其适宜的水文情势条件,基于水文情势条件优选出适宜不同类型湿地植物种植和生长的修复区域,是河湖湿地植物修复重建的关键研究内容。
淹水频率是指一定时间内某区域被水淹没的时间占总时间的比例。淹水频率作为重要的河湖水文情势指标,不仅可以表征干湿交替变化的程度,更是影响河湖湿地植物生长繁殖和群落结构演替的重要指标,例如苦草所需的淹水频率明显高于狗牙根,因此淹水频率是影响河湖湿地植物修复区域选取的关键因素。目前,由于实测水文站点有限,河湖水深实测数据的时空精度不足,难以实现整个河湖淹水频率的定量化分析,导致湿地植物修复区域的选取不合理,严重影响了湿地植物的成活率。
发明内容
针对现有技术中的上述不足,本发明提供的一种基于淹水频率的河湖湿地植物修复区域优选方法解决了现有缺乏一种对河湖湿地植物修复区域进行定量化优选的方法的问题。
为了达到上述发明目的,本发明采用的技术方案为:一种基于淹水频率的河湖湿地植物修复区域优选方法,包括以下步骤:
S1、根据河湖支流流量等数据,建立河湖水动力数学模型;
S2、采用河湖水动力数学模型,计算各个网格的淹水频率;
S3、根据各个网格的淹水频率,对河湖水域和河湖湿地植物修复区域分别进行分区,得到各个网格的河湖水文分区和湿地植物修复分区类型;
S4、根据各个网格的河湖水文分区和湿地植物修复分区类型,绘制基于淹水频率的河湖湿地植物修复分区图;
S5、根据河湖湿地植物修复分区图,在不同修复分区种植相应的河湖湿地植物。
进一步地,所述步骤S1包括以下分步骤:
S11、建立河湖水动力数学模型的控制方程;
S12、设定控制方程的边界条件和初始条件,初步构建河湖水动力数学模型;
S13、根据控制方程、边界条件和初始条件,调整河湖糙率,模拟得到水位模拟值;
S14、根据水位模拟值与水位实测值的拟合情况,对初步构建的河湖水动力数学模型进行验证和修正,得到构建完成的河湖水动力数学模型。
进一步地,所述步骤S11中河湖水动力数学模型的控制方程为:
其中,H为总水深,H=h+ζ,h为底床高程,ζ为自由水面高程,u为流速在x方向上的分量,v为流速在y方向上的分量,w为流速在z方向上的分量,Qu和Qv为水平扩散项,AV为垂向紊流粘性系数,mx和my分别为x和y方向坐标变换因子,m=mxmy,m为坐标变换因子之积,Z为转换后的高程,p为静水压力,f为科氏力系数,g为重力加速度,为对时间t求偏导,/>为对x方向求偏导,/>为对y方向求偏导,/>为对z方向求偏导。
进一步地,所述步骤S12中边界条件包括:河湖支流流量数据、河湖气象数据和河湖风场数据;
所述初始条件包括:河湖初始水位和河湖初始流速。
进一步地,所述步骤S2包括以下分步骤:
S21、设定模拟时段,采用河湖水动力数学模型模拟各个网格在模拟时段的各个时间的水深模拟值;
S22、根据各个网格在模拟时段的各个时间的水深模拟值,计算各个网格的淹水频率。
进一步地,所述步骤S22中各个网格的淹水频率的计算公式为:
其中,SFi为网格i的淹水频率,T为模拟时段,ti为网格i在模拟时段内被淹没的时间之和,单个网格在某个时间的水深模拟值大于0,则该网格在该时间被淹没。
进一步地,所述步骤S3包括以下分步骤:
S31、根据各个网格的淹水频率,对河湖水域进行水文分区,得到河湖水文分区;
S32、根据河湖水文分区和湿地植物生长对淹水的需求,对河湖湿地植物修复区域进行分区,得到河湖湿地植物修复分区;
S33、根据河湖水文分区、河湖湿地植物修复分区和各个网格的淹水频率,得到各个网格的河湖水文分区和湿地植物修复分区类型。
进一步地,所述步骤S31中河湖水文分区包括5类区域,分别为:长期淹没区、高频淹没区、干湿平衡区、高频出露区和长期出露区;
所述步骤S32中河湖湿地植物修复分区包括5类区域,分别为:耐长期水淹植物区、耐高频水淹植物区、干湿平衡植物区、耐高频出露植物区和耐长期出露植物区。
综上,本发明的有益效果为:
(1)、本发明通过河湖水动力数学模型构建、淹水频率计算、河湖水文分区划分,定量化优选出了不同类型湿地植物适宜的修复区域,共计耐长期水淹植物区、耐高频水淹植物区、干湿平衡植物区、耐高频出露植物区和耐长期出露植物区5个分区。通过该方法可有效提高湿地植物成活率,提升湿地植物修复效果,为河湖水生态***保护与修复提供定量化的技术支撑手段。
(2)、本发明通过构建河湖水动力数学模型,模拟得到各个网格的水深数据,据此计算淹水频率,最后根据淹水频率将河湖分为长期淹没区、高频淹没区、干湿平衡区、高频出露区和长期出露区,实现了针对河湖水文特征的定量化分区,可为河湖水文水资源精细化管理提供技术支撑。
附图说明
图1为一种基于淹水频率的河湖湿地植物修复区域优选方法的流程图;
图2为鄱阳湖湿地植物修复分区图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1所示,一种基于淹水频率的河湖湿地植物修复区域优选方法,包括以下步骤:
S1、根据河湖支流流量数据,建立河湖水动力数学模型;
步骤S1包括以下分步骤:
S11、建立河湖水动力数学模型的控制方程;
步骤S11中河湖水动力数学模型的控制方程为:
其中,H为总水深,H=h+ζ,h为底床高程,ζ为自由水面高程,u为流速在x方向上的分量,v为流速在y方向上的分量,w为流速在z方向上的分量,Qu和Qv为水平扩散项,AV为垂向紊流粘性系数,mx和my分别为x和y方向坐标变换因子,m=mxmy,m为坐标变换因子之积,Z为转换后的高程,p为静水压力,f为科氏力系数,g为重力加速度,为对时间t求偏导,/>为对x方向求偏导,/>为对y方向求偏导,/>为对z方向求偏导。
S12、设定控制方程的边界条件和初始条件,初步构建河湖水动力数学模型;
骤S12中边界条件包括:河湖支流流量数据、河湖气象数据和河湖风场数据;
初始条件包括:河湖初始水位和河湖初始流速。
S13、根据控制方程、边界条件和初始条件,调整河湖糙率,模拟得到水位模拟值;
S14、根据水位模拟值与水位实测值的拟合情况,对初步构建的河湖水动力数学模型进行验证和修正,得到构建完成的河湖水动力数学模型。
S2、采用河湖水动力数学模型,计算各个网格的淹水频率;
步骤S2包括以下分步骤:
S21、设定模拟时段,采用河湖水动力数学模型模拟各个网格在模拟时段的各个时间的水深模拟值;
S22、根据各个网格在模拟时段的各个时间的水深模拟值,计算各个网格的淹水频率。
步骤S22中各个网格的淹水频率的计算公式为:
其中,SFi为网格i的淹水频率,T为模拟时段,ti为网格i在模拟时段内被淹没的时间之和,单个网格在某个时间的水深模拟值大于0,则该网格在该时间被淹没。
S3、根据各个网格的淹水频率,对河湖水域和河湖湿地植物修复区域分别进行分区,得到各个网格的河湖水文分区和湿地植物修复分区类型;
步骤S3包括以下分步骤:
S31、根据各个网格的淹水频率,对河湖水域进行水文分区,得到河湖水文分区;
步骤S31中河湖水文分区包括5类区域,分别为:长期淹没区、高频淹没区、干湿平衡区、高频出露区和长期出露区。
S32、根据河湖水文分区和湿地植物生长对淹水的需求,对河湖湿地植物修复区域进行分区,得到河湖湿地植物修复分区;
所述步骤S32中河湖湿地植物修复分区包括5类区域,分别为:耐长期水淹植物区、耐高频水淹植物区、干湿平衡植物区、耐高频出露植物区和耐长期出露植物区。
不同修复分区适宜种植不同类型的湿地植物,例如耐长期水淹植物区适宜种植可耐受长期水淹的苦草等植物。
S33、根据河湖水文分区、河湖湿地植物修复分区和各个网格的淹水频率,得到各个网格的河湖水文分区和湿地植物修复分区类型。
河湖水文分区、河湖湿地植物修复分区和各个网格的淹水频率3者间的关系如表1:
表1
S4、根据各个网格的河湖水文分区和湿地植物修复分区类型,绘制基于淹水频率的河湖湿地植物修复分区图;
绘制时可采用ArcGIS等软件。
S5、根据河湖湿地植物修复分区图,在不同修复分区种植相应的河湖湿地植物。
实验:
鄱阳湖作为我国最大的淡水河湖,流域面积为16.22万km2,湖区面积达到50km2。鄱阳湖是长江中下游大型通江河湖,也是国际重要湿地,其水资源、水环境和水生态地位非常重要。由于湖区地形、江湖关系等影响,鄱阳湖水文特征呈现明显的空间异质性,主河道常年被水淹没,洲滩则处于干湿交替状态。
本实验计算了2003-2018年多年平均边界条件下的鄱阳湖全年淹水频率,绘制了鄱阳湖湿地植物修复分区图,见图2。从图可知,鄱阳湖存在耐长期水淹植物区、耐高频水淹植物区、干湿平衡植物区、耐高频出露植物区、耐长期出露植物区等修复分区,分别占整个湖区面积的47%、12%、19%、13%、9%。鄱阳湖湖区以耐长期水淹植物区为主,所占面积约为整个鄱阳湖的一半;其次为干湿平衡植物区,耐长期出露植物区所占面积最小。

Claims (2)

1.一种基于淹水频率的河湖湿地植物修复区域优选方法,其特征在于,包括以下步骤:
S1、根据河湖支流流量数据,建立河湖水动力数学模型:
S11、建立河湖水动力数学模型的控制方程,
S12、设定控制方程的边界条件和初始条件,初步构建河湖水动力数学模型,所述边界条件包括:河湖支流流量数据、河湖气象数据和河湖风场数据;所述初始条件包括:河湖初始水位和河湖初始流速,
S13、根据控制方程、边界条件和初始条件,调整河湖糙率,模拟得到水位模拟值,
S14、根据水位模拟值与水位实测值的拟合情况,对初步构建的河湖水动力数学模型进行验证和修正,得到构建完成的河湖水动力数学模型;
S2、采用河湖水动力数学模型,计算各个网格的淹水频率:
S21、设定模拟时段,采用河湖水动力数学模型模拟各个网格在模拟时段的各个时间的水深模拟值,
S22、根据各个网格在模拟时段的各个时间的水深模拟值,计算各个网格的淹水频率;
S3、根据各个网格的淹水频率,对河湖水域和河湖湿地植物修复区域分别进行分区,得到各个网格的河湖水文分区和湿地植物修复分区类型:
S31、根据各个网格的淹水频率,对河湖水域进行水文分区,得到河湖水文分区,
S32、根据河湖水文分区和湿地植物生长对淹水的需求,对河湖湿地植物修复区域进行分区,得到河湖湿地植物修复分区,
S33、根据河湖水文分区、河湖湿地植物修复分区和各个网格的淹水频率,得到各个网格的河湖水文分区和湿地植物修复分区类型;
S4、根据各个网格的河湖水文分区和湿地植物修复分区类型,绘制基于淹水频率的河湖湿地植物修复分区图;
S5、根据河湖湿地植物修复分区图,在不同修复分区种植相应的河湖湿地植物;
所述步骤S11中河湖水动力数学模型的控制方程为:
其中,H为总水深,H=h+ζ,h为底床高程,ζ为自由水面高程,u为流速在x方向上的分量,v为流速在y方向上的分量,w为流速在z方向上的分量,Qu和Qv为水平扩散项,AV为垂向紊流粘性系数,mx和my分别为x和y方向坐标变换因子,m=mxmy,m为坐标变换因子之积,Z为转换后的高程,p为静水压力,f为科氏力系数,g为重力加速度,为对时间t求偏导,/>为对x方向求偏导,/>为对y方向求偏导,/>为对z方向求偏导;
所述步骤S22中各个网格的淹水频率的计算公式为:
其中,SFi为网格i的淹水频率,T为模拟时段,ti为网格i在模拟时段内被淹没的时间之和,单个网格在某个时间的水深模拟值大于0,则该网格在该时间被淹没。
2.根据权利要求1所述的基于淹水频率的河湖湿地植物修复区域优选方法,其特征在于,所述步骤S31中河湖水文分区包括5类区域,分别为:长期淹没区、高频淹没区、干湿平衡区、高频出露区和长期出露区;
所述长期淹没区的淹水频率范围为80%≤SF i≤100%,高频淹没区的淹水频率范围为60%≤SF i<80%,干湿平衡区的淹水频率范围为40%≤SF i<60%,高频出露区的淹水频率范围为20%<SF i<40%,长期出露区的淹水频率范围为0%≤SF i≤20%;
所述步骤S32中河湖湿地植物修复分区包括5类区域,分别为:耐长期水淹植物区、耐高频水淹植物区、干湿平衡植物区、耐高频出露植物区和耐长期出露植物区。
CN202210245128.9A 2022-03-14 2022-03-14 一种基于淹水频率的河湖湿地植物修复区域优选方法 Active CN114611297B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210245128.9A CN114611297B (zh) 2022-03-14 2022-03-14 一种基于淹水频率的河湖湿地植物修复区域优选方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210245128.9A CN114611297B (zh) 2022-03-14 2022-03-14 一种基于淹水频率的河湖湿地植物修复区域优选方法

Publications (2)

Publication Number Publication Date
CN114611297A CN114611297A (zh) 2022-06-10
CN114611297B true CN114611297B (zh) 2023-11-24

Family

ID=81863636

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210245128.9A Active CN114611297B (zh) 2022-03-14 2022-03-14 一种基于淹水频率的河湖湿地植物修复区域优选方法

Country Status (1)

Country Link
CN (1) CN114611297B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108108911A (zh) * 2018-01-09 2018-06-01 中国科学院南京地理与湖泊研究所 基于水生态***健康的大型湖泊分区水质目标制定方法
CN111684970A (zh) * 2020-05-27 2020-09-22 中都工程设计有限公司 堤岸生态护坡的构建方法
CN112307420A (zh) * 2020-10-29 2021-02-02 大连海洋大学 一种河口湿地植被生境修复方法
CN113439617A (zh) * 2021-07-02 2021-09-28 广西壮族自治区自然资源生态修复中心 一种生态护岸消落带植物群种的设计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108108911A (zh) * 2018-01-09 2018-06-01 中国科学院南京地理与湖泊研究所 基于水生态***健康的大型湖泊分区水质目标制定方法
CN111684970A (zh) * 2020-05-27 2020-09-22 中都工程设计有限公司 堤岸生态护坡的构建方法
CN112307420A (zh) * 2020-10-29 2021-02-02 大连海洋大学 一种河口湿地植被生境修复方法
CN113439617A (zh) * 2021-07-02 2021-09-28 广西壮族自治区自然资源生态修复中心 一种生态护岸消落带植物群种的设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
谷娟等.鄱阳湖水体淹没频率变化及其湿地植被的响应.生态学报.2018,第38卷(第21期),第7718-7726页. *

Also Published As

Publication number Publication date
CN114611297A (zh) 2022-06-10

Similar Documents

Publication Publication Date Title
CN104933268B (zh) 一种基于一维非恒定流数值模型的洪水分析方法
CN111428401A (zh) 一种堰塞湖溃决过程模拟方法
US20210181375A1 (en) Numerical method for simulating a karez well in association with a groundwater model
CN116502775B (zh) 一种水文序列增强及预测方法
Zhang et al. The dynamic capacity calculation method and the flood control ability of the Three Gorges Reservoir
CN114580283A (zh) 山溪性强潮河口分汊段中长期动力地貌演变数值模拟方法
CN114611297B (zh) 一种基于淹水频率的河湖湿地植物修复区域优选方法
Kuang et al. Impact of large-scale reclamation on hydrodynamics and flushing in Victoria Harbour, Hong Kong
Zeng et al. Evaluating the hydrological effects of the Three Gorges Reservoir based on a large-scale coupled hydrological-hydrodynamic-dam operation model
Bellos et al. 2D flood modelling: The case of Tous dam break
Rojas et al. Macquarie River floodplain flow modeling: implications for ecogeomorphology
Bailly et al. Vegetation patch effects on flow resistance at channel scale
Huang et al. Prediction of design typhoon wind speeds and profiles using refined typhoon wind field model
Tang et al. Two-dimensional water environment numerical simulation research based on EFDC in Mudan River, Northeast China
CN117332709B (zh) 一种建筑风压模拟等效数值方法及模拟装置
CN101000337B (zh) 水土保持半比尺模型试验的实现方法
CN117010727B (zh) 一种基于水平年的径流序列一致性改正方法
CN117852850B (zh) 基于虚拟河道替代闸孔泄流计算的蓄滞洪区洪水调度方法及***
CN114459571B (zh) 一种基于初态水位曲线簇河道初始水位确定方法
CN118095562A (zh) 一种基于水文水动力模型的湖泊汛末蓄水策略智能优化方法
Xin-yu et al. Computer simulation of flood scheduling in large scale flood control systems
Lei Sr et al. Numerical simulation study on water environment of Puzhehei Lake
CN110717247B (zh) 一种考虑梯田影响的坡面产汇流过程模拟方法
Gu Study on the Application of the Drainage Pipe Network and River Channel Coupling Model in Urban Flood Control and Drainage
CN118212089A (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