CN109979172A - 一种基于新安江模型的动态山洪临界雨量预报方法 - Google Patents

一种基于新安江模型的动态山洪临界雨量预报方法 Download PDF

Info

Publication number
CN109979172A
CN109979172A CN201910283595.9A CN201910283595A CN109979172A CN 109979172 A CN109979172 A CN 109979172A CN 201910283595 A CN201910283595 A CN 201910283595A CN 109979172 A CN109979172 A CN 109979172A
Authority
CN
China
Prior art keywords
critical
rainfall
critical rainfall
mountain torrents
dynamic
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.)
Withdrawn
Application number
CN201910283595.9A
Other languages
English (en)
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.)
Hebei Meteorological Disaster Prevention Center
Nanjing University of Information Science and Technology
Original Assignee
Hebei Meteorological Disaster Prevention Center
Nanjing University of Information Science and Technology
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 Hebei Meteorological Disaster Prevention Center, Nanjing University of Information Science and Technology filed Critical Hebei Meteorological Disaster Prevention Center
Priority to CN201910283595.9A priority Critical patent/CN109979172A/zh
Publication of CN109979172A publication Critical patent/CN109979172A/zh
Withdrawn legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/14—Rainfall or precipitation gauges
    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B21/00—Alarms responsive to a single specified undesired or abnormal condition and not otherwise provided for
    • G08B21/02—Alarms for ensuring the safety of persons
    • G08B21/10—Alarms for ensuring the safety of persons responsive to calamitous events, e.g. tornados or earthquakes
    • G—PHYSICS
    • G08—SIGNALLING
    • G08B—SIGNALLING OR CALLING SYSTEMS; ORDER TELEGRAPHS; ALARM SYSTEMS
    • G08B31/00—Predictive alarm systems characterised by extrapolation or other computation using updated historic data
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Emergency Management (AREA)
  • Atmospheric Sciences (AREA)
  • Hydrology & Water Resources (AREA)
  • Computing Systems (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Ecology (AREA)
  • Environmental Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于新安江模型的动态山洪临界雨量预报方法，包括以下步骤：S1.采集发生山洪的小流域水文信息；S2.构建三水源新安江模型；S3.确定临界流量；S4.实时预报动态临界雨量。本发明通过三水源新安江模型计算考虑了流域前期水文条件的影响，得到实时动态的山洪临界雨量预报方法。采用二分法迭代计算临界雨量值是收敛的，优化了计算速度。本发明方法可以推广至各气象水文部门，取代目前普遍采用的临界雨量固定值预警指标。该方法有效提高了山洪临界雨量预警指标与自然条件的吻合程度，为进一步提高山洪灾害预警精度提供了有力的技术支撑。

Description

一种基于新安江模型的动态山洪临界雨量预报方法
技术领域
本发明属于山洪灾害防治研究技术领域，具体涉及一种基于新安江模型的 动态山洪临界雨量预报方法。
背景技术
我国山洪灾害发生频次高、损失程度大，并且山洪灾害已成为洪涝灾害中致人死亡的主 要灾种。山洪一般指山丘区中小流域由降雨引起的突发性、暴涨暴落的洪水。国内气象与水 文等各部门正在研发有效的山洪监测预警预报系统，力求使灾害程度达到最小。
目前，各气象部门在山洪灾害预警中多采用临界雨量作为预警指标。在目 前的实际应用中，所采用的山洪临界雨量是一个固定值。由于有前期雨量、流 域湿润程度、河道水位等因素影响，山洪灾害发生所需要的临界雨量不是一个 固定值，而是根据实际情况在随时变化的。采用临界雨量为固定值的预警模式 不能充分考虑到实际情况的影响，不能保证有效的预警精度，极易发生空报、 漏报、误报的现象。
发明内容
本发明所要解决的技术问题是针对上述现有技术的不足，提供一种基于新 安江模型的动态山洪临界雨量预报方法，充分考虑流域前期的降水、蒸散发、 土壤含水量、河道流量等影响，得到实时动态的临界雨量值作为山洪预警指 标。
为实现上述技术目的，本发明采取的技术方案为：
一种基于新安江模型的动态山洪临界雨量预报方法，包括以下步骤：
S1.采集发生山洪的小流域水文信息；
S2.构建三水源新安江模型；
S3.确定临界流量；
S4.实时预报动态临界雨量。
为优化上述技术方案，采取的具体措施还包括：
上述的步骤S1所述水文信息包括降水量、水面蒸发量和流域出口断面流量。
上述的步骤S2所述构建三水源新安江模型包括以下步骤：
S2.1.采用三层模型计算蒸散发，采用蓄满产流模型计算产流；
S2.2.用自由水蓄水库结构将总产流划分为地面径流、壤中流和地下径流；
S2.3.采用线性水库计算坡地汇流，采用滞后演算法和马斯京根分段连续演 算法计算河道汇流。
上述的步骤S3具体为：根据历史流量资料，采用水文频率分析的方法，以 一定年份为重现期的流量值作为临界流量值。
上述的步骤S4具体包括以下步骤：
S4.1.计算模型的预热期；
S4.2.延长假设临界雨量值及雨型分布；
S4.3.二分法迭代得到最终的临界雨量值。
上述的步骤S4.3具体为：
将预热期连同步骤S4.2延长的雨量、蒸发过程，输入新安江模型进行连续 计算，输出得到流量过程结果；
判断延长期的计算流量值的最大值与临界流量值的关系，采用二分法迭代得 到新的临界雨量值，返回步骤S4.2，直到精度满足要求，得到最终的临界雨量 值。
本发明具有以下有益效果：
本发明通过三水源新安江模型计算考虑了流域前期水文条件的影响，得到实 时动态的山洪临界雨量预报方法。
采用二分法迭代计算临界雨量值是收敛的，优化了计算速度。
本发明方法可以推广至各气象水文部门，取代目前普遍采用的临界雨量固定 值预警指标。该方法有效提高了山洪临界雨量预警指标与自然条件的吻合程度， 为进一步提高山洪灾害预警精度提供了有力的技术支撑。
附图说明
图1是本发明实施例的三水源新安江模型流程图；
图2是本发明实施例基于新安江模型的二分法迭代计算临界雨量流程图；
图3是本发明实施例模型预热期输出结果；
图4是本发明实施例各预警级别及时段的模型计算延长期流量过程；
图5是本发明实施例各预警级别及时段的临界雨量值。
具体实施方式
以下结合附图对本发明的实施例作进一步详细描述。
本发明的一种基于新安江模型的动态山洪临界雨量预报方法，包括以下步骤：
S1.采集发生山洪的小流域水文信息；
实施例中，首先需要掌握流域的气候、历史洪水、植被、地貌、地质结构等 特征。需要搜集的流域时间序列水文信息资料包括：降水量，水面蒸发，流域出 口断面流量。
S2.构建三水源新安江模型；
实施例中，根据S1中搜集的历史水文资料，构建三水源新安江模型，对模 型参数进行率定。
构建三水源新安江模型包括以下步骤：
S2.1.采用三层模型计算蒸散发，采用蓄满产流模型计算产流；
S2.2.用自由水蓄水库结构将总产流划分为地面径流、壤中流和地下径流；
S2.3.采用线性水库计算坡地汇流，采用滞后演算法和马斯京根分段连续演 算法计算河道汇流。
三水源新安江模型计算流程如图1所示。
模型的输入是降雨量和观测水面蒸发量，输出为实际蒸散发量和流域出口断 面的流量过程。流程图中框内的变量为模型计算中间变量，在线上的变量为模型 参数。
S3.确定临界流量；
实施例中，根据历史山洪流量信息，确定发生山洪的临界流量。临界流量有 时也称为漫滩流量，它是指发生山洪的流量最小值或称流量阈值。如果临界流量 值不容易确定，也可以根据历史流量资料，采用水文频率分析的方法，以一定年 份为重现期的流量值作为临界流量值。
S4.实时预报动态临界雨量。
实施例中，有了以上步骤中的已经率定好的新安江模型、临界流量，结合实 时观测的降雨、蒸发资料，便可以进行实时动态的山洪临界雨量预报，具体步骤 如下：
S4.1.计算模型的预热期
如果要对当前时刻t2的山洪临界雨量进行预报，那么需要从t2时刻之前的 t1时刻开始运行新安江模型。那么t1至t2时刻便称为模型的预热期。输入t1至 t2时刻的观测降雨、水面蒸发，通过新安江模型便可以计算得到流量过程，以及 t2时刻的初始状态。初始状态反映了流域前期的饱和程度。在新安江模型中，对 后面流量计算有影响的初始状态变量包括张力水蓄水容量W、自由水蓄水容量 S、产流面积比FR、前期各部分流量过程QS、QI、QG等。
为了考虑前期初始状态的影响，只要将预热期连同后期的资料一起输入新安 江模型进行连续计算，便可以方便快捷地使延长期计算结果考虑了前期(预热期) 水文条件的影响。
S4.2.假设临界雨量值的延长及雨型分布
在预热期t1至t2时刻的观测降雨资料后，延长假设一个临界雨量值以及在 各时段的雨量分布。一般在雨期的蒸散发可以忽略不计，所以水面蒸发输入的延 长都可以认为是零。
S4.3.二分法迭代得到最终的临界雨量值
将预热期连同S4.2延长的雨量、蒸发过程，输入新安江模型进行连续计算， 输出得到流量过程结果。判断延长期的计算流量值的最大值与临界流量值的关系， 采用二分法迭代得到新的临界雨量值，返回步骤S4.2，直到精度满足要求，得到 最终的临界雨量值。其流程如图2所示。
下面以河北省内丘县柳林镇的山洪沟为例说明具体实施方式。
在步骤S1中，搜集流域内历史的面平均雨量资料、水面蒸发资料、流量过 程资料。采用ArcGIS软件对数字高程分析得到流域面积大小。
该流域位于半干旱地区，年均降水量较小，历史洪水过程有陡涨抖落现象。
在步骤S2中，根据所搜集的历史资料建立流域的三水源新安江模型，模型 计算时段为1h，并率定模型参数。三水源新安江模型计算流程如图1所示，可 分为蒸散发计算、产流计算、分水源计算和汇流计算四个层次结构，分别编写相 应的计算模块进行调用。
在步骤S3中，这里分为4级(Lv1-Lv4)预警，所以这里从小到大确定4个 各级的临界流量值。
有了以上建模准备和临界流量值的确定，便可以进行实时动态山洪临界雨量 预报计算，步骤如下：
在S4.1步骤中，设以2016年7月18日至8月16日为预热期，分别计算8 月16日各预警级别(Lv1-Lv4)下1h，3h，6h，12h，24h为时段的山洪临界雨 量。
在步骤S4.2中，在预热期资料的基础上对雨量资料进行延长。这里的雨型 先按平均化处理，也就是将临界雨量平均分配到每个计算时段。如有需要也可以 按其他权重雨型进行分配。
在步骤S4.3中，设置初始的临界雨量下限值与上限值，设置迭代精度可以 为1mm。按照图2流程进行二分法迭代计算，最终可以得到各预警级别及时段 的临界雨量预报值。这样的临界雨量值是一个动态输出值，考虑了流域前期水文 条件的影响。
按照步骤S4.1中设置的时间，图3为程序输出的预热期图形结果。图3上 图为降雨过程，中图为计算流量与实测流量过程对比，下图为土壤张力水含量百 分比变化过程。图4为程序输出的临界雨量条件下延长期的预报流量过程图，其 中横坐标为时段数，纵坐标为流量值(m3/s),可以发现迭代过程是收敛的，最终 结果到达临界流量值。图5为程序输出的各预警级别及时段下的临界雨量值，可 以很直观的进行动态山洪临界雨量预警预报。
以上仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实 施例，凡属于本发明思路下的技术方案均属于本发明的保护范围。应当指出， 对于本技术领域的普通技术人员来说，在不脱离本发明原理前提下的若干改进 和润饰，应视为本发明的保护范围。

Claims (6)

1.一种基于新安江模型的动态山洪临界雨量预报方法，其特征在于：包括以下步骤：
S1.采集发生山洪的小流域水文信息；
S2.构建三水源新安江模型；
S3.确定临界流量；
S4.实时预报动态临界雨量。
2.根据权利要求1所述的一种基于新安江模型的动态山洪临界雨量预报方法，其特征在于：步骤S1所述水文信息包括降水量、水面蒸发量和流域出口断面流量。
3.根据权利要求1所述的一种基于新安江模型的动态山洪临界雨量预报方法，其特征在于：步骤S2所述构建三水源新安江模型包括以下步骤：
S2.1.采用三层模型计算蒸散发，采用蓄满产流模型计算产流；
S2.2.用自由水蓄水库结构将总产流划分为地面径流、壤中流和地下径流；
S2.3.采用线性水库计算坡地汇流，采用滞后演算法和马斯京根分段连续演算法计算河道汇流。
4.根据权利要求1所述的一种基于新安江模型的动态山洪临界雨量预报方法，其特征在于：步骤S3具体为：根据历史流量资料，采用水文频率分析的方法，以一定年份为重现期的流量值作为临界流量值。
5.根据权利要求1所述的一种基于新安江模型的动态山洪临界雨量预报方法，其特征在于：步骤S4具体包括以下步骤：
S4.1.计算模型的预热期；
S4.2.延长假设临界雨量值及雨型分布；
S4.3.二分法迭代得到最终的临界雨量值。
6.根据权利要求1所述的一种基于新安江模型的动态山洪临界雨量预报方法，其特征在于：步骤S4.3具体为：
将预热期连同步骤S4.2延长的雨量、蒸发过程，输入新安江模型进行连续计算，输出得到流量过程结果；
判断延长期的计算流量值的最大值与临界流量值的关系，采用二分法迭代得到新的临界雨量值，返回步骤S4.2，直到精度满足要求，得到最终的临界雨量值。
CN201910283595.9A 2019-04-09 2019-04-09 一种基于新安江模型的动态山洪临界雨量预报方法 Withdrawn CN109979172A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910283595.9A CN109979172A (zh) 2019-04-09 2019-04-09 一种基于新安江模型的动态山洪临界雨量预报方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910283595.9A CN109979172A (zh) 2019-04-09 2019-04-09 一种基于新安江模型的动态山洪临界雨量预报方法

Publications (1)

Publication Number Publication Date
CN109979172A true CN109979172A (zh) 2019-07-05

Family

ID=67083901

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910283595.9A Withdrawn CN109979172A (zh) 2019-04-09 2019-04-09 一种基于新安江模型的动态山洪临界雨量预报方法

Country Status (1)

Country Link
CN (1) CN109979172A (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110459036A (zh) * 2019-09-09 2019-11-15 四川省水利科学研究院 一种基于深度学习的山洪预警方法
CN111047213A (zh) * 2019-12-25 2020-04-21 长江水利委员会水文局 基于多要素联合判别的中长期水资源量临界预警指标方法
CN112396297A (zh) * 2020-11-03 2021-02-23 华中科技大学 一种洪水过程遭遇时间和量级发生规律的解析方法和系统
CN116523292A (zh) * 2023-02-03 2023-08-01 中国科学院精密测量科学与技术创新研究院 基于致污临界雨量的农业面源污染风险动态识别预警方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007025962A (ja) * 2005-07-14 2007-02-01 Nec Corp 情報提供システムおよび防災情報提供システムならびにその情報提供方法
CN101604356A (zh) * 2009-06-24 2009-12-16 南京大学 一种不确定性中长期水文预报模型的建立方法
CN102542732A (zh) * 2010-12-08 2012-07-04 中国电信股份有限公司 自然灾害监测预警方法与系统
CN102704439A (zh) * 2012-01-11 2012-10-03 戴会超 一种水库分期设计洪水的推求方法
JP2012203747A (ja) * 2011-03-25 2012-10-22 Yasuaki Iwai 表示装置、描画用プログラム、防災システム
CN103473892A (zh) * 2013-08-23 2013-12-25 成都理工大学 沟道起动型泥石流发生预警方法及应用
CN104408900A (zh) * 2014-11-10 2015-03-11 柳州师范高等专科学校 基于动态优化的神经网络洪水预警装置及其方法
CN105760666A (zh) * 2016-02-04 2016-07-13 南京信息工程大学 一种基于dem的汇水区临界雨量计算方法
CN105894741A (zh) * 2016-05-04 2016-08-24 南京信息工程大学 一种多资源融合的洪水灾害监测预警装置及方法
CN107145678A (zh) * 2017-05-22 2017-09-08 中国水利水电科学研究院 一种河道平面二维水流模型糙率的率定方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2007025962A (ja) * 2005-07-14 2007-02-01 Nec Corp 情報提供システムおよび防災情報提供システムならびにその情報提供方法
CN101604356A (zh) * 2009-06-24 2009-12-16 南京大学 一种不确定性中长期水文预报模型的建立方法
CN102542732A (zh) * 2010-12-08 2012-07-04 中国电信股份有限公司 自然灾害监测预警方法与系统
JP2012203747A (ja) * 2011-03-25 2012-10-22 Yasuaki Iwai 表示装置、描画用プログラム、防災システム
CN102704439A (zh) * 2012-01-11 2012-10-03 戴会超 一种水库分期设计洪水的推求方法
CN103473892A (zh) * 2013-08-23 2013-12-25 成都理工大学 沟道起动型泥石流发生预警方法及应用
CN104408900A (zh) * 2014-11-10 2015-03-11 柳州师范高等专科学校 基于动态优化的神经网络洪水预警装置及其方法
CN105760666A (zh) * 2016-02-04 2016-07-13 南京信息工程大学 一种基于dem的汇水区临界雨量计算方法
CN105894741A (zh) * 2016-05-04 2016-08-24 南京信息工程大学 一种多资源融合的洪水灾害监测预警装置及方法
CN107145678A (zh) * 2017-05-22 2017-09-08 中国水利水电科学研究院 一种河道平面二维水流模型糙率的率定方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
李青: ""基于洪峰模数的山洪灾害雨量预警指标研究"", 《地球信息科学学报》 *
王江婷: ""基于HEC_HMS模型的北方典型小流域山洪预警模拟与研究"", 《中国优秀硕士学位论文全文数据库工程科技I辑》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110459036A (zh) * 2019-09-09 2019-11-15 四川省水利科学研究院 一种基于深度学习的山洪预警方法
CN111047213A (zh) * 2019-12-25 2020-04-21 长江水利委员会水文局 基于多要素联合判别的中长期水资源量临界预警指标方法
CN112396297A (zh) * 2020-11-03 2021-02-23 华中科技大学 一种洪水过程遭遇时间和量级发生规律的解析方法和系统
CN112396297B (zh) * 2020-11-03 2021-06-29 华中科技大学 一种洪水过程遭遇时间和量级发生规律的解析方法和系统
CN116523292A (zh) * 2023-02-03 2023-08-01 中国科学院精密测量科学与技术创新研究院 基于致污临界雨量的农业面源污染风险动态识别预警方法

Similar Documents

Publication Publication Date Title
CN109979172A (zh) 一种基于新安江模型的动态山洪临界雨量预报方法
CN106884405B (zh) 一种无资料地区溃堤型山洪灾害分析评价方法
CN108597189B (zh) 基于临界雨量的分布式中小流域地质灾害及洪水预警方法
CN111582755B (zh) 一种基于多维度集合信息山洪灾害综合风险动态评估方法
CN106202790B (zh) 一种分布式河北模型构建方法及其应用
CN105912770A (zh) 水文实时预报系统
CN108345735A (zh) 一种流域水文模型参数校准方法
CN110274656B (zh) 一种城市内河水位预报预警方法
CN106706033A (zh) 一种海绵城市性能监视系统和方法
CN110009002A (zh) 基于多维降雨特征空间最优决策的山洪快速预警预报方法
Jiang et al. Geomorphic evolution of the Yellow River Delta: Quantification of basin-scale natural and anthropogenic impacts
CN110459036B (zh) 一种基于深度学习的山洪预警方法
CN106529739A (zh) 基于网格化的洪水预测方法
CN115130396A (zh) 一种河道型水库库区的分布式水文模型建模方法
CN108416468A (zh) 一种洪水超前预警预报方法
CN109118718A (zh) 泥石流发生降雨i-d曲线阈值构建方法、流域泥石流预警方法
CN106295859A (zh) 一种无径流资料地区汇流计算方法
CN105160121A (zh) 一种有限元控制的分布式水文模型的建模方法
CN114139258A (zh) 城市内涝建模评估方法、系统及计算机可读存储介质
CN108269199A (zh) 一种面向对象的小水库群时空分布式出流计算方法
CN112528563A (zh) 一种基于svm算法的城市内涝预警方法
CN103810401B (zh) 一种人类活动影响分离的二维径流还原方法
CN113919125A (zh) 一种基于区域产汇流耦合模型系统的防洪预报调度方法
CN106446359A (zh) 基于流式大数据处理方式的降雨径流预测计算方法
Haibo et al. Application of Synthetic Unit Hydrograph on HEC-HMS Model for flood forecasting

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
WW01 Invention patent application withdrawn after publication

Application publication date: 20190705

WW01 Invention patent application withdrawn after publication