CN115391745B - 一种基于概率匹配平均法的降水预报订正方法及系统 - Google Patents

一种基于概率匹配平均法的降水预报订正方法及系统 Download PDF

Info

Publication number
CN115391745B
CN115391745B CN202211322464.5A CN202211322464A CN115391745B CN 115391745 B CN115391745 B CN 115391745B CN 202211322464 A CN202211322464 A CN 202211322464A CN 115391745 B CN115391745 B CN 115391745B
Authority
CN
China
Prior art keywords
sub
precipitation
historical
monitoring points
forecast
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
CN202211322464.5A
Other languages
English (en)
Other versions
CN115391745A (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.)
Guoneng Daduhe Big Data Service Co ltd
Original Assignee
Guoneng Daduhe Big Data Service Co ltd
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 Guoneng Daduhe Big Data Service Co ltd filed Critical Guoneng Daduhe Big Data Service Co ltd
Priority to CN202211322464.5A priority Critical patent/CN115391745B/zh
Publication of CN115391745A publication Critical patent/CN115391745A/zh
Application granted granted Critical
Publication of CN115391745B publication Critical patent/CN115391745B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/10—Devices for predicting weather conditions
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01W—METEOROLOGY
    • G01W1/00—Meteorology
    • G01W1/14—Rainfall or precipitation gauges
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Biodiversity & Conservation Biology (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Atmospheric Sciences (AREA)
  • Ecology (AREA)
  • Pure & Applied Mathematics (AREA)
  • Environmental Sciences (AREA)
  • Mathematical Physics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Databases & Information Systems (AREA)
  • Hydrology & Water Resources (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本申请实施例公开了一种基于概率匹配平均法的降水预报订正方法及系统，属于降水预报技术领域，包括：获取目标区域的多个监测点的相关信息；基于相关信息，对多个监测点进行聚类，将目标区域分为多个子区域；获取目标区域在未来目标时间段的初始降水预报；基于概率匹配平均法，对初始降水预报进行订正，获取订正后降水预报；对于每个子区域，将订正后降水预报分为多个子时间段的单元数据；通过降水订正模型，基于子区域包括的监测点的相关信息，按照时间先后顺序对多个子时间段的单元数据逐次进行二次修订，获取子区域的目标降水预报，具有提高降水预测的准确度的优点。

Description

一种基于概率匹配平均法的降水预报订正方法及系统
技术领域
本发明主要涉及降水预报技术领域，具体地说，涉及一种基于概率匹配平均法的降水预报订正方法及系统。
背景技术
在近几十年来国内外已经对气候变化(包括降水)开展了广泛的研究，其已经成为全球性的问题，并成为国内外最为关注的研究热点。全球气候模式(GCMs)，是研究气候变化机理和预测未来气候变化趋势的重要工具，但由于GCM空间分辨率较为粗糙，为了获得更加精细化的区域气候变化情况，需要对其进行降尺度处理，从而更加精确的预测区域降水变化情况。
目前主要有3种方式，分别为(1)提高GCM模型输出结果精度；(2)将GCM与高分辨率的区域气候模式(RCM)进行嵌套，即动力降尺度方法；(3)建立GCM输出大尺度气象变量与区域气象因子间的统计关系，即统计降尺度方法。
提高GCM模型精度将导致模型计算量呈指数增长，模型计算及运行需要通过超级计算机来进行相关操作。同时动力降尺度方法计算量大，模拟受收入的边界条件影响，会继承GCMs的误差和不足，同样受局地特征的异质性强迫。此外，动力降尺度模型在不同地区应用过程中，还存在着模拟结果与实测值间表现不稳定、数据化结果不能很好反映实际情况等问题。
因此，需要一种基于概率匹配平均法的降水预报订正方法及系统，用于提高降水预测的准确度。
发明内容
本说明书实施例之一提供一种基于概率匹配平均法的降水预报订正方法，所述方法包括获取目标区域的多个监测点的相关信息；基于所述相关信息，对所述多个监测点进行聚类，将所述目标区域分为多个子区域；获取所述目标区域在未来目标时间段的初始降水预报；基于概率匹配平均法，对所述初始降水预报进行订正，获取订正后降水预报；对于每个所述子区域，将所述订正后降水预报分为多个子时间段的单元数据；通过降水订正模型，基于所述子区域包括的监测点的相关信息，按照时间先后顺序对所述多个子时间段的单元数据逐次进行二次修订，获取所述子区域的目标降水预报。
在一些实施例中，所述基于概率匹配平均法，对所述初始降水预报进行订正，获取订正后降水预报，包括：获取所述目标区域的历史观测降水预报；基于所述目标区域的历史观测降水预报确定多个降水值对应的订正阈值；通过插值法基于所述多个降水值对应的订正阈值，确定降水值与订正阈值的对应曲线；基于所述降水值与订正阈值的对应曲线，对所述初始降水预报进行订正，获取所述订正后降水预报。
在一些实施例中，所述基于所述相关信息，对所述多个监测点进行聚类，将所述目标区域分为多个子区域，包括：对于任意两个监测点，基于所述两个监测点的相关信息，确定所述两个监测点之间的相似度；基于所述相似度，对所述多个监测点进行聚类，将所述目标区域分为多个子区域。
在一些实施例中，所述通过降水订正模型，基于所述子区域包括的监测点的相关信息，按照时间先后顺序对所述多个子时间段的单元数据逐次进行二次修订，获取所述子区域的目标降水预报，包括：对于每个所述子时间段，通过所述降水订正模型基于所述子时间段的前一子时间段的二次修订后的单元数据及所述子区域包括的监测点的相关信息，对所述子时间段的单元数据进行二次修订。
在一些实施例中，所述监测点的相关信息至少包括所述监测点的地形信息、历史降水序列及历史动力热物理因子矩阵，其中，所述历史降水序列包括所述监测点在多个历史时间点的降水量，所述历史动力热物理因子矩阵由多个历史动力热物理因子序列组成，其中，一个所述历史动力热物理因子序列包括所述监测点的一种动力热物理因子在所述多个历史时间点的值；所述通过所述降水订正模型基于所述子时间段的前一子时间段的二次修订后的单元数据及所述子区域包括的监测点的相关信息，对所述子时间段的单元数据进行二次修订，包括：基于所述子区域包括的监测点的地形信息、历史降水序列及历史动力热物理因子矩阵和所述子区域在所述子时间段的前一子时间段的动力热物理因子矩阵，预测所述子区域在所述子时间段的动力热物理因子矩阵；通过所述降水订正模型基于所述子时间段的前一子时间段的二次修订后的单元数据及预测的所述子区域在所述子时间段的动力热物理因子矩阵，对所述子时间段的单元数据进行二次修订。
本说明书实施例之一提供一种基于概率匹配平均法的降水预报订正系统，所述系统包括：信息获取模块，用于获取目标区域的多个监测点的相关信息；区域分割模块，用于基于所述相关信息，对所述多个监测点进行聚类，将所述目标区域分为多个子区域；数据获取模块，用于获取所述目标区域在未来目标时间段的初始降水预报；降水订正模块，用于基于概率匹配平均法，对所述初始降水预报进行订正，获取订正后降水预报，还用于对于每个所述子区域，将所述订正后降水预报分为多个子时间段的单元数据，通过降水订正模型，基于所述子区域包括的监测点的相关信息，按照时间先后顺序对所述多个子时间段的单元数据逐次进行二次修订，获取所述子区域的目标降水预报。
在一些实施例中，所述降水订正模块还用于：获取所述目标区域的历史观测降水预报；基于所述目标区域的历史观测降水预报确定多个降水值对应的订正阈值；通过插值法基于所述多个降水值对应的订正阈值，确定降水值与订正阈值的对应曲线；基于所述降水值与订正阈值的对应曲线，对所述初始降水预报进行订正，获取订正后降水预报。
在一些实施例中，所述区域分割模块还用于：对于任意两个监测点，基于所述两个监测点的相关信息，确定所述两个监测点之间的相似度；基于所述相似度，对所述多个监测点进行聚类，将所述目标区域分为多个子区域。
在一些实施例中，所述降水订正模块还用于：对于每个所述子时间段，通过所述降水订正模型基于所述子时间段的前一子时间段的二次修订后的单元数据及所述子区域包括的监测点的相关信息，对所述子时间段的单元数据进行二次修订。
在一些实施例中，所述监测点的相关信息至少包括所述监测点的地形信息、历史降水序列及历史动力热物理因子矩阵，其中，所述历史降水序列包括所述监测点在多个历史时间点的降水量，所述历史动力热物理因子矩阵由多个历史动力热物理因子序列组成，其中，一个所述历史动力热物理因子序列包括所述监测点的一种动力热物理因子在所述多个历史时间点的值；所述降水订正模块还用于：基于所述子区域包括的监测点的地形信息、历史降水序列及历史动力热物理因子矩阵和所述子区域在所述子时间段的前一子时间段的动力热物理因子矩阵，预测所述子区域在所述子时间段的动力热物理因子矩阵；通过所述降水订正模型基于所述子时间段的前一子时间段的二次修订后的单元数据及预测的所述子区域在所述子时间段的动力热物理因子矩阵，对所述子时间段的单元数据进行二次修订。
本说明书提供的一种基于概率匹配平均法的降水预报订正方法及系统，至少具有以下有益效果：
1、基于概率匹配平均法，对初始降水预报进行订正，获取订正后降水预报，相对于单一模式能够提供降水预报技巧，消减系统性的误差；既能保留降水集合平均的空间分布，也能保留集合成员的概率分布，好于集合平均方法；过滤掉了不可预报的小尺度信息，提高降水预报订正数据精准度；
2、通过降水订正模型，基于子区域包括的监测点的相关信息，按照时间先后顺序对多个子时间段的单元数据逐次进行二次修订，获取子区域的目标降水预报，进一步提高目标降水预报的精准度。
附图说明
本申请将以示例性实施例的方式进一步说明，这些示例性实施例将通过附图进行详细描述。这些实施例并非限制性的，在这些实施例中，相同的编号表示相同的结构，其中：
图1是根据本申请一些实施例所示的一种基于概率匹配平均法的降水预报订正系统的应用场景示意图；
图2是根据本申请一些实施例所示的一种基于概率匹配平均法的降水预报订正系统的模块示意图；
图3是根据本申请一些实施例所示的一种基于概率匹配平均法的降水预报订正方法的示例性流程图。
图中，100、应用场景；110、处理设备；120、网络；130、用户终端；140、存储设备。
具体实施方式
为了更清楚地说明本申请的实施例的技术方案，下面将对实施例描述中所需要使用的附图作简单的介绍。显而易见地，下面描述中的附图仅仅是本申请的一些示例或实施例，对于本领域的普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图将本申请应用于其他类似情景。应当理解，给出这些示例性的实施例仅仅是为了使相关领域的技术人员能够更好地理解进而实现本发明，而并非以任何方式限制本发明的范围。除非从语言环境中显而易见或另做说明，图中相同标号代表相同结构或操作。
应当理解，本文使用的“系统”、“装置”、“单元”和/或“模块”是用于区分不同级别的不同组件、元件、部件、部分或装配的一种方法。然而，如果其他词语可实现相同的目的，则可通过其他表达来替换所述词语。
如本申请和权利要求书中所示，除非上下文明确提示例外情形，“一”、“一个”、“一种”和/或“该”等词并非特指单数，也可包括复数。一般说来，术语“包括”与“包含”仅提示包括已明确标识的步骤和元素，而这些步骤和元素不构成一个排它性的罗列，方法或者设备也可能包含其它的步骤或元素。
虽然本申请对根据本申请的实施例的系统中的某些模块或单元做出了各种引用，然而，任何数量的不同模块或单元可以被使用并运行在客户端和/或服务器上。所述模块仅是说明性的，并且所述系统和方法的不同方面可以使用不同模块。
本申请中使用了流程图用来说明根据本申请的实施例的系统所执行的操作。应当理解的是，前面或后面操作不一定按照顺序来精确地执行。相反，可以按照倒序或同时处理各个步骤。同时，也可以将其他操作添加到这些过程中，或从这些过程移除某一步或数步操作。
图1是根据本申请一些实施例所示的一种基于概率匹配平均法的降水预报订正系统的应用场景100示意图。
如图1所示，应用场景100可以包括处理设备110、网络120、用户终端130和存储设备140。
在一些实施例中，处理设备110可以用于处理与降水预报订正相关的信息和/或数据。例如，处理设备110可以获取目标区域的多个监测点的相关信息，基于相关信息，对多个监测点进行聚类，将目标区域分为多个子区域，获取目标区域在未来目标时间段的初始降水预报，基于概率匹配平均法，对初始降水预报进行订正，获取订正后降水预报，对于每个子区域，将订正后降水预报分为多个子时间段的单元数据，通过降水订正模型，基于子区域包括的监测点的相关信息，按照时间先后顺序对多个子时间段的单元数据逐次进行二次修订，获取子区域的目标降水预报。
在一些实施例中，处理设备110可以是区域的或者远程的。例如，处理设备110可以通过网络120访问存储于用户终端130和存储设备140中的信息和/或资料。在一些实施例中，处理设备110可以直接与用户终端130和存储设备140连接以访问存储于其中的信息和/或资料。在一些实施例中，处理设备110可以在云平台上执行。例如，该云平台可以包括私有云、公共云、混合云、社区云、分散式云、内部云等中的一种或其任意组合。在一些实施例中，处理设备110可以包含处理器，处理器可以包含一个或多个子处理器（例如，单芯处理设备或多核多芯处理设备）。仅仅作为范例，处理器可包含中央处理器（CPU）、专用集成电路（ASIC）、专用指令处理器（ASIP）、图形处理器（GPU）、物理处理器（PPU）、数字信号处理器（DSP）、现场可编程门阵列（FPGA）、可编辑逻辑电路（PLD）、控制器、微控制器单元、精简指令集电脑（RISC）、微处理器等或以上任意组合。
网络120可促进应用场景100中数据和/或信息的交换。在一些实施例中，应用场景100中的一个或多个组件（例如，处理设备110、用户终端130和存储设备140）可以通过网络120发送数据和/或信息给应用场景100中的其他组件。例如，处理设备110可以通过网络120从存储设备140获取目标排水系统的历史运行数据。在一些实施例中，网络120可以是任意类型的有线或无线网络。例如，网络120可以包括缆线网络、有线网络、光纤网络、电信网络、内部网络、网际网络、区域网络、广域网络、无线区域网络、都会区域网络、公共电话交换网络、蓝牙网络、ZigBee网络、近场通讯网络等或以上任意组合。
用户终端130可以获取应用场景100中的信息或数据，用户（可以是用户终端130的使用者。在一些实施例中，用户终端130可以通过网络120与应用场景100中的一个或多个组件（例如，处理设备110或存储设备140）进行数据和/或信息的交换。例如，用户终端130可以通过网络120从处理设备110获取子区域的目标降水预报。在一些实施例中，用户终端130可以包括移动装置、平板电脑、笔记本电脑等中的一种或其任意组合。在一些实施例中，移动装置可以包括可穿戴装置、智能行动装置、虚拟实境装置、增强实境装置等或其任意组合。
在一些实施例中，存储设备140可以与网络120连接以实现与应用场景100的一个或多个组件（例如，处理设备110、用户终端130等）通讯。应用场景100的一个或多个组件可以通过网络120访问存储于存储设备140中的资料或指令。在一些实施例中，存储设备140可以直接与应用场景100中的一个或多个组件（如，处理设备110、用户终端130）连接或通讯。
应该注意的是，上述描述仅出于说明性目的而提供，并不旨在限制本申请的范围。对于本领域普通技术人员而言，在本申请内容的指导下，可做出多种变化和修改。可以以各种方式组合本申请描述的示例性的实施例的特征、结构、方法和其他特征，以获得另外的和/或替代的示例性的实施例。例如，存储设备140可以是包括云计算平台的数据存储设备，例如公共云、私有云、社区和混合云等。然而，这些变化与修改不会背离本申请的范围。
图2是根据本申请一些实施例所示的一种基于概率匹配平均法的降水预报订正系统的模块示意图。如图2所示，基于概率匹配平均法的降水预报订正系统可以包括信息获取模块、区域分割模块、数据获取模块及降水订正模块。
信息获取模块可以用于获取目标区域的多个监测点的相关信息。在一些实施例中，监测点的相关信息至少包括监测点的地形信息、历史降水序列及历史动力热物理因子矩阵，其中，历史降水序列包括监测点在多个历史时间点的降水量，历史动力热物理因子矩阵由多个历史动力热物理因子序列组成，其中，一个历史动力热物理因子序列包括监测点的一种动力热物理因子在多个历史时间点的值。
区域分割模块可以用于基于相关信息，对多个监测点进行聚类，将目标区域分为多个子区域。在一些实施例中，区域分割模块还可以用于：对于任意两个监测点，基于两个监测点的相关信息，确定两个监测点之间的相似度；基于相似度，对多个监测点进行聚类，将目标区域分为多个子区域。
数据获取模块可以用于获取目标区域在未来目标时间段的初始降水预报。
降水订正模块可以用于基于概率匹配平均法，对初始降水预报进行订正，获取订正后降水预报，还用于对于每个子区域，将订正后降水预报分为多个子时间段的单元数据，通过降水订正模型，基于子区域包括的监测点的相关信息，按照时间先后顺序对多个子时间段的单元数据逐次进行二次修订，获取子区域的目标降水预报。在一些实施例中，降水订正模块还可以用于：获取目标区域的历史观测降水预报，基于目标区域的历史观测降水预报确定多个降水值对应的订正阈值，通过插值法基于多个降水值对应的订正阈值，确定降水值与订正阈值的对应曲线，基于降水值与订正阈值的对应曲线，对初始降水预报进行订正，获取订正后降水预报。在一些实施例中，降水订正模块还可以用于：对于每个子时间段，通过降水订正模型基于子时间段的前一子时间段的二次修订后的单元数据及子区域包括的监测点的相关信息，对子时间段的单元数据进行二次修订。在一些实施例中，降水订正模块还可以用于：基于子区域包括的监测点的地形信息、历史降水序列及历史动力热物理因子矩阵和子区域在子时间段的前一子时间段的动力热物理因子矩阵，预测子区域在子时间段的动力热物理因子矩阵；通过降水订正模型基于子时间段的前一子时间段的二次修订后的单元数据及预测的子区域在子时间段的动力热物理因子矩阵，对子时间段的单元数据进行二次修订。
关于信息获取模块、区域分割模块、数据获取模块及降水订正模块的更多描述可以参见图3及其相关描述，此处不再赘述。
图3是根据本申请一些实施例所示的一种基于概率匹配平均法的降水预报订正方法的示例性流程图。在一些实施例中，基于概率匹配平均法的降水预报订正方法可以由基于概率匹配平均法的降水预报订正系统执行。如图3所示，基于概率匹配平均法的降水预报订正方法可以包括如下步骤。
步骤310，获取目标区域的多个监测点的相关信息。在一些实施例中，步骤310可以由信息获取模块执行。
相关信息至少包括监测点的位置信息、地形信息、历史降水序列及历史动力热物理因子矩阵，其中，历史降水序列包括监测点在多个历史时间点的降水量，历史动力热物理因子矩阵由多个历史动力热物理因子序列组成，其中，一个历史动力热物理因子序列包括监测点的一种动力热物理因子在多个历史时间点的值。
地形信息可以至少包括高程、坡度及坡向。
多个历史动力热物理因子序列至少包括广义位温历史序列、广义位涡历史序列、力管涡度历史序列、二阶位涡历史序列、二阶湿位涡历史序列、湿热力平流参数历史序列、螺旋度历史序列、热力垂直通量散度历史序列、Q矢量历史序列、位势变形历史序列、对流涡度矢量历史序列、波作用密度历史序列、广义湿位涡历史序列、热力螺旋度历史序列、Q矢量散度历史序列、位势切变形变历史序列、位势伸展形变历史序列、热力散度垂直通量历史序列及涡度散度梯度历史序列。
一种动力热物理因子的历史序列可以包括监测点的一种动力热物理因子在多个历史时间点的值，例如，广义位温历史序列可以包括监测点的广义位温在多个历史时间点的值。又例如，二阶位涡历史序列可以包括监测点的二阶位涡在多个历史时间点的值。
在一些实施例中，信息获取模块可以从处理设备110、用户终端130、存储设备140和/或外部数据源获取监测点的相关信息。
步骤320，基于相关信息，对多个监测点进行聚类，将目标区域分为多个子区域。在一些实施例中，步骤320可以由区域分割模块执行。
在一些实施例中，区域分割模块基于相关信息，对多个监测点进行聚类，将目标区域分为多个子区域，可以包括：
对于任意两个监测点，基于两个监测点的相关信息，确定两个监测点之间的相似度，即基于两个监测点的位置信息、地形信息、历史降水序列及历史动力热物理因子矩阵，确定两个监测点之间的相似度；
基于相似度，对多个监测点进行聚类，将目标区域分为多个子区域。
在一些实施例中，区域分割模块可以通过相似度确定模型确定两个监测点的相关信息，相似度确定模型的输入为两个监测点的位置信息、地形信息、历史降水序列及历史动力热物理因子矩阵，相似度确定模型的输出为两个监测点之间的相似度，其中，相似度确定模型可以为卷积神经网络（CNN）、深度神经网络（DNN）、循环神经网络（RNN）、多层神经网络（MLP）、对抗神经网络（GAN）等一种或其任意组合。例如，相似度确定模型可以为卷积神经网络和深度神经网络组合形成的模型。
在一些实施例中，区域分割模块可以通过K-means算法基于相似度，对多个监测点进行聚类，将目标区域分为多个子区域，对多个监测点进行聚类，将目标区域分为多个子区域。可以理解的，聚类后生成的每个簇中包括的所有的监测点组成一个子区域。
步骤330，获取目标区域在未来目标时间段的初始降水预报。在一些实施例中，步骤330可以由数据获取模块执行。
未来目标时间段可以为未来的某个时间段。
在一些实施例中，预测数据获取模块可以从用户终端130、存储设备140或外部数据源（例如，多卫星降水联合反演IMERG卫星降水系统、GRAPES全球预报系统、中国气象局等）获取目标区域在未来目标时间段的初始降水预报。
步骤340，基于概率匹配平均法，对初始降水预报进行订正，获取订正后降水预报。在一些实施例中，步骤340可以由降水订正模块执行。
在一些实施例中，降水订正模块基于概率匹配平均法，对初始降水预报进行订正，获取订正后降水预报，可以包括：
获取目标区域的历史观测降水预报；
基于目标区域的历史观测降水预报确定多个降水值对应的订正阈值；
通过插值法基于多个降水值对应的订正阈值，确定降水值与订正阈值的对应曲线；
基于降水值与订正阈值的对应曲线，对初始降水预报进行订正，获取订正后降水预报。
假设X1，X2，…，Xi，…，Xn为一系列离散降水值，例如，0.1，1，5，10，15，25，35，50，70，90mm，相应的观测降水累积概率值为P0(X1），P0(X2），…，P0(Xi），…，P0(Xn）。与P0(X1），P0(X2），…，P0(Xi），…，P0(Xn）概率值对应的模式降水预报值记为X
Figure DEST_PATH_IMAGE002
1，X
Figure 596600DEST_PATH_IMAGE002
2，…，X
Figure 326789DEST_PATH_IMAGE002
i，…，X
Figure 434423DEST_PATH_IMAGE002
n，相应的模式预报降水累积概率值记为Pf（X
Figure 406796DEST_PATH_IMAGE002
1），Pf（X
Figure 690010DEST_PATH_IMAGE002
2），…，Pf（X
Figure 574789DEST_PATH_IMAGE002
i），…，Pf（X
Figure 920451DEST_PATH_IMAGE002
n）。对某一特定降水值Xi，假设观测累积概率分布与模式累积概率分布一致，则令Pf（X
Figure DEST_PATH_IMAGE004
i）
Figure DEST_PATH_IMAGE006
P0(Xi），已知Xi、P0(Xi）计算可得X
Figure 227673DEST_PATH_IMAGE002
i的值，以此类推，可得0.1，1，5，10，15，25，35，50，70，90mm各自对应的订正阈值，再采取插值法即可获得任意降水量的订正阈值，从而确定降水值与订正阈值的对应曲线。
步骤350，对于每个子区域，将订正后降水预报分为多个子时间段的单元数据，通过降水订正模型，基于子区域包括的监测点的相关信息，按照时间先后顺序对多个子时间段的单元数据逐次进行二次修订，获取子区域的目标降水预报。在一些实施例中，步骤350可以由降水订正模块执行。
在一些实施例中，降水订正模块通过降水订正模型，基于子区域包括的监测点的相关信息，按照时间先后顺序对多个子时间段的单元数据逐次进行二次修订，获取子区域的目标降水预报，可以包括：
对于每个子时间段，
通过降水订正模型基于子时间段的前一子时间段的二次修订后的单元数据及子区域包括的监测点的相关信息，对子时间段的单元数据进行二次修订。
在一些实施例中，降水订正模块通过降水订正模型基于子时间段的前一子时间段的二次修订后的单元数据及子区域包括的监测点的相关信息，对子时间段的单元数据进行二次修订，可以包括：
基于子区域包括的监测点的地形信息、历史降水序列及历史动力热物理因子矩阵和子区域在子时间段的前一子时间段的动力热物理因子矩阵，预测子区域在子时间段的动力热物理因子矩阵；
通过降水订正模型基于子时间段的前一子时间段的二次修订后的单元数据及预测的子区域在子时间段的动力热物理因子矩阵，对子时间段的单元数据进行二次修订。
在一些实施例中，降水订正模块可以通过任意方式基于子区域包括的监测点的地形信息、历史降水序列及历史动力热物理因子矩阵和子区域在子时间段的前一子时间段的动力热物理因子矩阵，预测子区域在子时间段的动力热物理因子矩阵。例如，降水订正模块可以通过数据预测模型基于子区域包括的监测点的地形信息、历史降水序列及历史动力热物理因子矩阵和子区域在子时间段的前一子时间段的动力热物理因子矩阵，预测子区域在子时间段的动力热物理因子矩阵，数据预测模型的输入可以包括监测点的地形信息、历史降水序列及历史动力热物理因子矩阵和子区域在子时间段的前一子时间段的动力热物理因子矩阵，数据预测模型的输出为子区域在子时间段的动力热物理因子矩阵，其中，数据预测模型可以为卷积神经网络（CNN）、深度神经网络（DNN）、循环神经网络（RNN）、多层神经网络（MLP）、对抗神经网络（GAN）等一种或其任意组合。例如，数据预测模型可以为卷积神经网络和深度神经网络组合形成的模型。
在一些实施例中，降水订正模块可以通过降水订正模型基于子时间段的前一子时间段的二次修订后的单元数据及预测的子区域在子时间段的动力热物理因子矩阵，对子时间段的单元数据进行二次修订，其中，降水订正模型的输入包括子时间段的前一子时间段的二次修订后的单元数据及预测的子区域在子时间段的动力热物理因子矩阵，降水订正模型的输出包括子区域的目标降水预报，其中，降水订正模型可以为卷积神经网络（CNN）、深度神经网络（DNN）、循环神经网络（RNN）、多层神经网络（MLP）、对抗神经网络（GAN）等一种或其任意组合。例如，降水订正模型可以为卷积神经网络和深度神经网络组合形成的模型。
在本申请的另一些实施例中，提供了一种基于概率匹配平均法的降水预报订正装置，包括至少一个处理设备以及至少一个存储设备；至少一个存储设备用于存储计算机指令，至少一个处理设备用于执行计算机指令中的至少部分指令以实现如上的一种基于概率匹配平均法的降水预报订正方法。
在本申请的又一些实施例中，提供了一种计算机可读存储介质，存储介质存储计算机指令，当计算机指令被处理设备执行时实现如上的基于概率匹配平均法的降水预报订正方法。
上文已对基本概念做了描述，显然，对于本领域技术人员来说，上述详细披露仅仅作为示例，而并不构成对本申请的限定。虽然此处并没有明确说明，本领域技术人员可能会对本申请进行各种修改、改进和修正。该类修改、改进和修正在本申请中被建议，所以该类修改、改进、修正仍属于本申请示范实施例的精神和范围。
同时，本申请使用了特定词语来描述本申请的实施例。如“一个实施例”、“一实施例”、和/或“一些实施例”意指与本申请至少一个实施例相关的某一特征、结构或特点。因此，应强调并注意的是，本说明书中在不同位置两次或多次提及的“一实施例”或“一个实施例”或“一个替代性实施例”并不一定是指同一实施例。此外，本申请的一个或多个实施例中的某些特征、结构或特点可以进行适当的组合。
此外，本领域技术人员可以理解，本申请的各方面可以通过若干具有可专利性的种类或情况进行说明和描述，包括任何新的和有用的工序、机器、产品或物质的组合，或对他们的任何新的和有用的改进。相应地，本申请的各个方面可以完全由硬件执行、可以完全由软件（包括固件、常驻软件、微码等）执行、也可以由硬件和软件组合执行。以上硬件或软件均可被称为“数据块”、“模块”、“引擎”、“单元”、“组件”或“系统”。此外，本申请的各方面可能表现为位于一个或多个计算机可读介质中的计算机产品，该产品包括计算机可读程序编码。
计算机存储介质可能包含一个内含有计算机程序编码的传播数据信号，例如在基带上或作为载波的一部分。该传播信号可能有多种表现形式，包括电磁形式、光形式等，或合适的组合形式。计算机存储介质可以是除计算机可读存储介质之外的任何计算机可读介质，该介质可以通过连接至一个指令执行系统、装置或设备以实现通讯、传播或传输供使用的程序。位于计算机存储介质上的程序编码可以通过任何合适的介质进行传播，包括无线电、电缆、光纤电缆、RF、或类似介质，或任何上述介质的组合。
本申请各部分操作所需的计算机程序编码可以用任意一种或多种程序语言编写，包括面向对象编程语言如Java、Scala、Smalltalk、Eiffel、JADE、Emerald、C++、C#、VB.NET、Python等，常规程序化编程语言如C语言、Visual Basic、Fortran 2003、Perl、COBOL 2002、PHP、ABAP，动态编程语言如Python、Ruby和Groovy，或其他编程语言等。该程序编码可以完全在用户计算机上运行、或作为独立的软件包在用户计算机上运行、或部分在用户计算机上运行部分在远程计算机运行、或完全在远程计算机或服务器上运行。在后种情况下，远程计算机可以通过任何网络形式与用户计算机连接，比如局域网（LAN）或广域网（WAN），或连接至外部计算机（例如通过因特网），或在云计算环境中，或作为服务使用如软件即服务（SaaS）。
此外，除非权利要求中明确说明，本申请处理元素和序列的顺序、数字字母的使用、或其他名称的使用，并非用于限定本申请流程和方法的顺序。尽管上述披露中通过各种示例讨论了一些目前认为有用的发明实施例，但应当理解的是，该类细节仅起到说明的目的，附加的权利要求并不仅限于披露的实施例，相反，权利要求旨在覆盖所有符合本申请实施例实质和范围的修正和等价组合。例如，虽然以上所描述的系统组件可以通过硬件设备实现，但是也可以只通过软件的解决方案得以实现，如在现有的服务器或移动设备上安装所描述的系统。
同理，应当注意的是，为了简化本申请披露的表述，从而帮助对一个或多个发明实施例的理解，前文对本申请实施例的描述中，有时会将多种特征归并至一个实施例、附图或对其的描述中。但是，这种披露方法并不意味着本申请对象所需要的特征比权利要求中提及的特征多。实际上，实施例的特征要少于上述披露的单个实施例的全部特征。
最后，应当理解的是，本申请中所述实施例仅用以说明本申请实施例的原则。其他的变形也可能属于本申请的范围。因此，作为示例而非限制，本申请实施例的替代配置可视为与本申请的教导一致。相应地，本申请的实施例不仅限于本申请明确介绍和描述的实施例。

Claims (10)

1.一种基于概率匹配平均法的降水预报订正方法，其特征在于，包括：
获取目标区域的多个监测点的相关信息；
基于所述相关信息，对所述多个监测点进行聚类，将所述目标区域分为多个子区域；
获取所述目标区域在未来目标时间段的初始降水预报；
基于概率匹配平均法，对所述初始降水预报进行订正，获取订正后降水预报；
对于每个所述子区域，
将所述订正后降水预报分为多个子时间段的单元数据；
通过降水订正模型，基于所述子区域包括的监测点的相关信息，按照时间先后顺序对所述多个子时间段的单元数据逐次进行二次修订，获取所述子区域的目标降水预报；
其中，所述监测点的相关信息至少包括监测点的位置信息、地形信息、历史降水序列及历史动力热物理因子矩阵，所述历史动力热物理因子矩阵由多个历史动力热物理因子序列组成，一个历史动力热物理因子序列包括监测点的一种动力热物理因子在多个历史时间点的值，所述多个历史动力热物理因子序列至少包括广义位温历史序列、广义位涡历史序列、力管涡度历史序列、二阶位涡历史序列、二阶湿位涡历史序列、湿热力平流参数历史序列、螺旋度历史序列、热力垂直通量散度历史序列、Q矢量历史序列、位势变形历史序列、对流涡度矢量历史序列、波作用密度历史序列、广义湿位涡历史序列、热力螺旋度历史序列、Q矢量散度历史序列、位势切变形变历史序列、位势伸展形变历史序列、热力散度垂直通量历史序列及涡度散度梯度历史序列；
所述基于所述相关信息，对所述多个监测点进行聚类，将所述目标区域分为多个子区域，包括：
对于任意两个监测点，基于两个监测点的相关信息，确定两个监测点之间的相似度，即基于两个监测点的位置信息、地形信息、历史降水序列及历史动力热物理因子矩阵，确定两个监测点之间的相似度；
基于相似度，对多个监测点进行聚类，将目标区域分为多个子区域。
2.根据权利要求1所述的一种基于概率匹配平均法的降水预报订正方法，其特征在于，所述基于概率匹配平均法，对所述初始降水预报进行订正，获取订正后降水预报，包括：
获取所述目标区域的历史观测降水预报；
基于所述目标区域的历史观测降水预报确定多个降水值对应的订正阈值；
通过插值法基于所述多个降水值对应的订正阈值，确定降水值与订正阈值的对应曲线；
基于所述降水值与订正阈值的对应曲线，对所述初始降水预报进行订正，获取所述订正后降水预报。
3.根据权利要求1所述的一种基于概率匹配平均法的降水预报订正方法，其特征在于，所述基于所述相关信息，对所述多个监测点进行聚类，将所述目标区域分为多个子区域，包括：
对于任意两个监测点，基于所述两个监测点的相关信息，确定所述两个监测点之间的相似度；
基于所述相似度，对所述多个监测点进行聚类，将所述目标区域分为多个子区域。
4.根据权利要求1-3任意一项所述的一种基于概率匹配平均法的降水预报订正方法，其特征在于，所述通过降水订正模型，基于所述子区域包括的监测点的相关信息，按照时间先后顺序对所述多个子时间段的单元数据逐次进行二次修订，获取所述子区域的目标降水预报，包括：
对于每个所述子时间段，
通过所述降水订正模型基于所述子时间段的前一子时间段的二次修订后的单元数据及所述子区域包括的监测点的相关信息，对所述子时间段的单元数据进行二次修订。
5.根据权利要求4所述的一种基于概率匹配平均法的降水预报订正方法，其特征在于，所述监测点的相关信息至少包括所述监测点的地形信息、历史降水序列及历史动力热物理因子矩阵，其中，所述历史降水序列包括所述监测点在多个历史时间点的降水量，所述历史动力热物理因子矩阵由多个历史动力热物理因子序列组成，其中，一个所述历史动力热物理因子序列包括所述监测点的一种动力热物理因子在所述多个历史时间点的值；
所述通过所述降水订正模型基于所述子时间段的前一子时间段的二次修订后的单元数据及所述子区域包括的监测点的相关信息，对所述子时间段的单元数据进行二次修订，包括：
基于所述子区域包括的监测点的地形信息、历史降水序列及历史动力热物理因子矩阵和所述子区域在所述子时间段的前一子时间段的动力热物理因子矩阵，预测所述子区域在所述子时间段的动力热物理因子矩阵；
通过所述降水订正模型基于所述子时间段的前一子时间段的二次修订后的单元数据及预测的所述子区域在所述子时间段的动力热物理因子矩阵，对所述子时间段的单元数据进行二次修订。
6.一种基于概率匹配平均法的降水预报订正系统，其特征在于，包括：
信息获取模块，用于获取目标区域的多个监测点的相关信息；
区域分割模块，用于基于所述相关信息，对所述多个监测点进行聚类，将所述目标区域分为多个子区域；
数据获取模块，用于获取所述目标区域在未来目标时间段的初始降水预报；
降水订正模块，用于基于概率匹配平均法，对所述初始降水预报进行订正，获取订正后降水预报，还用于对于每个所述子区域，将所述订正后降水预报分为多个子时间段的单元数据，通过降水订正模型，基于所述子区域包括的监测点的相关信息，按照时间先后顺序对所述多个子时间段的单元数据逐次进行二次修订，获取所述子区域的目标降水预报；
其中，所述监测点的相关信息至少包括监测点的位置信息、地形信息、历史降水序列及历史动力热物理因子矩阵，所述历史动力热物理因子矩阵由多个历史动力热物理因子序列组成，一个历史动力热物理因子序列包括监测点的一种动力热物理因子在多个历史时间点的值，所述多个历史动力热物理因子序列至少包括广义位温历史序列、广义位涡历史序列、力管涡度历史序列、二阶位涡历史序列、二阶湿位涡历史序列、湿热力平流参数历史序列、螺旋度历史序列、热力垂直通量散度历史序列、Q矢量历史序列、位势变形历史序列、对流涡度矢量历史序列、波作用密度历史序列、广义湿位涡历史序列、热力螺旋度历史序列、Q矢量散度历史序列、位势切变形变历史序列、位势伸展形变历史序列、热力散度垂直通量历史序列及涡度散度梯度历史序列；
所述基于所述相关信息，对所述多个监测点进行聚类，将所述目标区域分为多个子区域，包括：
对于任意两个监测点，基于两个监测点的相关信息，确定两个监测点之间的相似度，即基于两个监测点的位置信息、地形信息、历史降水序列及历史动力热物理因子矩阵，确定两个监测点之间的相似度；
基于相似度，对多个监测点进行聚类，将目标区域分为多个子区域。
7.根据权利要求6所述的一种基于概率匹配平均法的降水预报订正系统，其特征在于，所述降水订正模块还用于：
获取所述目标区域的历史观测降水预报；
基于所述目标区域的历史观测降水预报确定多个降水值对应的订正阈值；
通过插值法基于所述多个降水值对应的订正阈值，确定降水值与订正阈值的对应曲线；
基于所述降水值与订正阈值的对应曲线，对所述初始降水预报进行订正，获取订正后降水预报。
8.根据权利要求6所述的一种基于概率匹配平均法的降水预报订正系统，其特征在于，所述区域分割模块还用于：
对于任意两个监测点，基于所述两个监测点的相关信息，确定所述两个监测点之间的相似度；
基于所述相似度，对所述多个监测点进行聚类，将所述目标区域分为多个子区域。
9.根据权利要求6-8任意一项所述的一种基于概率匹配平均法的降水预报订正系统，其特征在于，所述降水订正模块还用于：
对于每个所述子时间段，
通过所述降水订正模型基于所述子时间段的前一子时间段的二次修订后的单元数据及所述子区域包括的监测点的相关信息，对所述子时间段的单元数据进行二次修订。
10.根据权利要求9所述的一种基于概率匹配平均法的降水预报订正系统，其特征在于，所述监测点的相关信息至少包括所述监测点的地形信息、历史降水序列及历史动力热物理因子矩阵，其中，所述历史降水序列包括所述监测点在多个历史时间点的降水量，所述历史动力热物理因子矩阵由多个历史动力热物理因子序列组成，其中，一个所述历史动力热物理因子序列包括所述监测点的一种动力热物理因子在所述多个历史时间点的值；
所述降水订正模块还用于：
基于所述子区域包括的监测点的地形信息、历史降水序列及历史动力热物理因子矩阵和所述子区域在所述子时间段的前一子时间段的动力热物理因子矩阵，预测所述子区域在所述子时间段的动力热物理因子矩阵；
通过所述降水订正模型基于所述子时间段的前一子时间段的二次修订后的单元数据及预测的所述子区域在所述子时间段的动力热物理因子矩阵，对所述子时间段的单元数据进行二次修订。
CN202211322464.5A 2022-10-27 2022-10-27 一种基于概率匹配平均法的降水预报订正方法及系统 Active CN115391745B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211322464.5A CN115391745B (zh) 2022-10-27 2022-10-27 一种基于概率匹配平均法的降水预报订正方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211322464.5A CN115391745B (zh) 2022-10-27 2022-10-27 一种基于概率匹配平均法的降水预报订正方法及系统

Publications (2)

Publication Number Publication Date
CN115391745A CN115391745A (zh) 2022-11-25
CN115391745B true CN115391745B (zh) 2023-01-20

Family

ID=84128006

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211322464.5A Active CN115391745B (zh) 2022-10-27 2022-10-27 一种基于概率匹配平均法的降水预报订正方法及系统

Country Status (1)

Country Link
CN (1) CN115391745B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116228046B (zh) * 2023-05-09 2023-07-18 成都信息工程大学 一种基于卫星遥感和地理数据的山区空间降水估算方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107632329A (zh) * 2017-09-04 2018-01-26 贵州东方世纪科技股份有限公司 一种虚拟雨量站雨量计算方法
WO2020223951A1 (zh) * 2019-05-09 2020-11-12 中国科学院深圳先进技术研究院 一种多模式集成降水预报方法及装置
WO2022007380A1 (zh) * 2020-07-06 2022-01-13 浙江大学 面向"互联网+"环境机械产品个性化设计模式匹配方法
WO2022163864A1 (ja) * 2021-02-01 2022-08-04 サントリーホールディングス株式会社 情報処理装置、水資源管理方法、情報処理方法、及び記録媒体

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109636032B (zh) * 2018-12-11 2021-08-24 中国科学院深圳先进技术研究院 基于多模式集成的降水预报方法、系统、终端及存储介质
CN112241604A (zh) * 2020-11-02 2021-01-19 国网湖南省电力有限公司 概率匹配串并联耦合多模型的电网暴雨灾害预报订正方法
CN113033957B (zh) * 2021-02-26 2023-10-27 兰州中心气象台(兰州干旱生态环境监测预测中心) 一种多模式降水预报及实时动态检验评估系统
CN114648181B (zh) * 2022-05-24 2022-11-18 国能大渡河大数据服务有限公司 一种基于机器学习的降水预报修正方法及系统

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107632329A (zh) * 2017-09-04 2018-01-26 贵州东方世纪科技股份有限公司 一种虚拟雨量站雨量计算方法
WO2020223951A1 (zh) * 2019-05-09 2020-11-12 中国科学院深圳先进技术研究院 一种多模式集成降水预报方法及装置
WO2022007380A1 (zh) * 2020-07-06 2022-01-13 浙江大学 面向"互联网+"环境机械产品个性化设计模式匹配方法
WO2022163864A1 (ja) * 2021-02-01 2022-08-04 サントリーホールディングス株式会社 情報処理装置、水資源管理方法、情報処理方法、及び記録媒体

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
基于ECMWF模式的定量降水客观订正方法;郑婧等;《高原气象》;20200624(第04期);830-839 *
基于频率匹配法的中国降水多模式预报订正研究;智协飞等;《大气科学学报》;20191128(第06期);814-823 *

Also Published As

Publication number Publication date
CN115391745A (zh) 2022-11-25

Similar Documents

Publication Publication Date Title
US10439594B2 (en) Actually-measured marine environment data assimilation method based on sequence recursive filtering three-dimensional variation
CN112418482B (zh) 一种基于时间序列聚类的云计算能耗预测方法
JP7021732B2 (ja) 時系列予測装置、時系列予測方法及びプログラム
CN108074015B (zh) 一种风电功率超短期预测方法及系统
CN114648181B (zh) 一种基于机器学习的降水预报修正方法及系统
CN115391745B (zh) 一种基于概率匹配平均法的降水预报订正方法及系统
JP2016161314A (ja) 同化データ決定装置、気象予測装置、同化データ決定方法およびプログラム
CN105184398A (zh) 一种电力最大负荷小样本预测方法
CN111105050B (zh) 风机维护计划的生成方法、装置、设备及存储介质
CN109614384A (zh) Hadoop框架下的电力系统短期负荷预测方法
CN116663395A (zh) 基于参数优化的支持向量机回归的降雨量等值面生成方法
CN115438590A (zh) 一种基于bp神经网络的降水预报订正方法及系统
CN113609763B (zh) 基于不确定性的卫星组件布局温度场预测方法
CN110968929A (zh) 风电场风速的预测方法、装置及电子设备
Zjavka Multi-site post-processing of numerical forecasts using a polynomial network substitution for the general differential equation based on operational calculus
CN111323847A (zh) 用于为模拟集成算法确定权重比的方法和设备
CN114118508A (zh) 一种基于时空卷积网络的od市场航空客流量预测方法
CN117634704A (zh) 一种用于旅游景点的选址优化方法、系统、设备及介质
CN110275895B (zh) 一种缺失交通数据的填充设备、装置及方法
CN115047544A (zh) 一种覆冰预警方法及系统
CN117390364B (zh) 一种用于运营状态下桥梁响应的时空域建模拟合方法
CN117851736B (zh) 一种基于模糊自适应寻优融合的气象要素插值方法
CN115600749B (zh) 地下水位预测方法、装置及电子设备
CN107330567A (zh) 基于大数据技术的配网短期负载预测方法
US11263338B2 (en) Data security maintenance method for data analysis application

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