CN113673139A - 一种双波长粒子谱反演方法 - Google Patents

一种双波长粒子谱反演方法 Download PDF

Info

Publication number
CN113673139A
CN113673139A CN202110936121.7A CN202110936121A CN113673139A CN 113673139 A CN113673139 A CN 113673139A CN 202110936121 A CN202110936121 A CN 202110936121A CN 113673139 A CN113673139 A CN 113673139A
Authority
CN
China
Prior art keywords
wavelength
data
extinction coefficient
aerosol
ratio
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.)
Pending
Application number
CN202110936121.7A
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.)
Hefei Zhongke Guangbo Quantum Technology Co ltd
Original Assignee
Hefei Zhongke Guangbo Quantum Technology 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 Hefei Zhongke Guangbo Quantum Technology Co ltd filed Critical Hefei Zhongke Guangbo Quantum Technology Co ltd
Priority to CN202110936121.7A priority Critical patent/CN113673139A/zh
Publication of CN113673139A publication Critical patent/CN113673139A/zh
Pending legal-status Critical Current

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
    • G06F30/25—Design optimisation, verification or simulation using particle-based methods
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00—Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88—Lidar systems specially adapted for specific applications
    • G01S17/95—Lidar systems specially adapted for specific applications for meteorological use
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/10—Numerical modelling
    • 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)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Remote Sensing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Radar, Positioning & Navigation (AREA)
  • General Health & Medical Sciences (AREA)
  • Electromagnetism (AREA)
  • Biophysics (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Artificial Intelligence (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Biomedical Technology (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)
  • Optical Radar Systems And Details Thereof (AREA)

Abstract

本发明涉及粒子谱反演，具体涉及一种双波长粒子谱反演方法，读取双波长激光雷达的原始数据，对波长数据进行预处理后，检验波长数据是否有效，若波长数据无效，则舍弃无效波长数据，并返回S1，否则对波长数据进行消光系数反演，计算消光系数比值，并在数据库中查找对应波长内消光系数比值的气溶胶有效直径，若波长内消光系数比值在数据库中不可查，则舍弃并返回S1，否则输出对应的气溶胶有效直径，通过气溶胶模型进行修正，修正后输出粒子谱，完成反演；本发明提供的技术方案能够有效克服现有技术所存在的硬件设备结构复杂、计算量大的缺陷。

Description

一种双波长粒子谱反演方法
技术领域
本发明涉及粒子谱反演，具体涉及一种双波长粒子谱反演方法。
背景技术
气溶胶粒子的谱分布，是指对应于不同半径气溶胶粒子浓度的分布。研究气溶胶粒子谱可以更好地揭示颗粒物浓度、雾霾成因、粒径谱分布等，对于大气污染治理和气候成因变化有重要的科研作用。目前，主要采集气溶胶粒径谱仪进行探测，但是这种设备只能探测单点的值，若要探测区域范围内的气溶胶粒子谱分布，就需要使用多台设备进行探测。激光雷达作为一种主动遥感设备，可以对空间范围内的气溶胶进行探测，但是气溶胶粒子谱的反演是一大难题。
现有技术中，有通过三波长(355nm、532nm、1064nm)来反演粒子谱分布，硬件设备的结构复杂，成本较高，并且通过GCV正则化计算反演，计算复杂，计算量大。
发明内容
(一)解决的技术问题
针对现有技术所存在的上述缺点，本发明提供了一种双波长粒子谱反演方法，能够有效克服现有技术所存在的硬件设备结构复杂、计算量大的缺陷。
(二)技术方案
为实现以上目的，本发明通过以下技术方案予以实现：
一种双波长粒子谱反演方法，包括以下步骤：
S1、读取双波长激光雷达的原始数据，对波长数据进行预处理后，检验波长数据是否有效；
S2、若波长数据无效，则舍弃无效波长数据，并返回S1，否则对波长数据进行消光系数反演；
S3、计算消光系数比值，并在数据库中查找对应波长内消光系数比值的气溶胶有效直径；
S4、若波长内消光系数比值在数据库中不可查，则舍弃并返回S1，否则输出对应的气溶胶有效直径；
S5、通过气溶胶模型进行修正，修正后输出粒子谱，完成反演。
优选地，S1中所述对波长数据进行预处理，包括：
将截取原始数据中最后一段数据的平均值作为背景值，利用数据序列减去背景值对波长数据进行扣除背景；
通过手动修正原始数据中近地面段数据，对波长数据进行几何因子修正。
优选地，S2中所述对波长数据进行消光系数反演，包括：
求解米散射激光雷达方程，通过假设激光雷达比，用Fernald方法计算得到消光系数。
优选地，所述米散射激光雷达方程如下：
Figure BDA0003212915570000021
式中，P(r)为激光雷达接收的高度r处的大气后向散射回波功率，Pt为激光雷达发射的激光功率，k为激光雷达系统常数，βm(r)为空气分子的后向散射系数，βa(r)为大气气溶胶的后向散射系数，αm(r’)为空气分子的消光系数，αa(r’)为大气气溶胶的消光系数；
大气气溶胶消光系数通过下式计算得到：
Figure BDA0003212915570000022
式中，Sa是大气气溶胶的消光后向散射比，或称激光雷达比，Sm是空气分子的消光后向散射比。
优选地，S3中所述计算消光系数比值，包括：
对计算出两个波长对应的消光系数相除，得到一组比值序列。
优选地，所述数据库中消光系数比值与气溶胶有效直径之间的对应关系通过模拟计算得到。
优选地，所述气溶胶模型为通过历史探测数据统计得出的不同天气、不同气溶胶类型下粒子有效直径的真值范围。
(三)有益效果
与现有技术相比，本发明所提供的一种双波长粒子谱反演方法，通过双波长进行消光系数反演，不仅降低了硬件设备结构的复杂程度，同时也减少了计算量；通过模拟计算出不同波长消光系数比值与气溶胶有效直径之间的对应关系，建立数据库，从而能够准确查找出气溶胶有效直径；通过历史探测数据建立不同气溶胶模型，对反演结果进行修正，提高了反演的精度和准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍。显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的流程示意图；
图2为337nm、1060nm双波长下消光系数比值与气溶胶有效直径之间的关系曲线。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述。显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
一种双波长粒子谱反演方法，如图1所示，S1、读取双波长激光雷达的原始数据，对波长数据进行预处理后，检验波长数据是否有效。
其中，对波长数据进行预处理，包括：
将截取原始数据中最后一段数据的平均值作为背景值，利用数据序列减去背景值对波长数据进行扣除背景；
通过手动修正原始数据中近地面段数据，对波长数据进行几何因子修正。
当设备异常或者探测时极端天气会导致波长数据无效，此时应该舍弃无效波长数据。
S2、若波长数据无效，则舍弃无效波长数据，并返回S1，否则对波长数据进行消光系数反演。
其中，对波长数据进行消光系数反演，包括：
求解米散射激光雷达方程，通过假设激光雷达比，用Fernald方法计算得到消光系数。
米散射激光雷达方程如下：
Figure BDA0003212915570000041
式中，P(r)为激光雷达接收的高度r处的大气后向散射回波功率，Pt为激光雷达发射的激光功率，k为激光雷达系统常数，βm(r)为空气分子的后向散射系数，βa(r)为大气气溶胶的后向散射系数，αm(r’)为空气分子的消光系数，αa(r’)为大气气溶胶的消光系数。
大气气溶胶消光系数通过下式计算得到：
Figure BDA0003212915570000042
式中，Sa是大气气溶胶的消光后向散射比，或称激光雷达比，Sm是空气分子的消光后向散射比。
S3、计算消光系数比值，并在数据库中查找对应波长内消光系数比值的气溶胶有效直径。
S4、若波长内消光系数比值在数据库中不可查，则舍弃并返回S1，否则输出对应的气溶胶有效直径。
其中，计算消光系数比值，包括：
对计算出两个波长对应的消光系数相除，得到一组比值序列。
数据库中消光系数比值与气溶胶有效直径之间的对应关系通过模拟计算得到。
S5、通过气溶胶模型进行修正，修正后输出粒子谱，完成反演。
气溶胶模型为通过历史探测数据统计得出的不同天气、不同气溶胶类型下粒子有效直径的真值范围。
下面以335nm、1064nm双波长和532nm、1064nm双波长对本申请技术方案的反演过程进行描述。
1)双波长激光雷达发射波长为355nm、1064nm，激光雷达原始数据包含355nm波长通道探测数据、1064nm波长通道探测数据。首先，读取双波长激光雷达的原始数据，对波长355nm和波长1064nm的数据进行扣除背景、几何因子修正，然后对修正数据进行检验是否有效，舍弃无效数据，对有效数据进行消光系数反演，消光系数反演使用fernald方法，计算出355nm波长和1064nm波长对应的消光系数后相除，得到一组比值序列，在数据库中查找355nm/1064nm消光系数比值的气溶胶有效直径；若该波长和消光系数比值在数据库内不可查，则舍弃，若可查，则输出查找值并通过气溶胶模型进行修正，修正后输出粒子谱值，反演结束。
2)双波长激光雷达发射波长为532nm、1064nm，激光雷达原始数据包含532nm波长通道探测数据、1064nm波长通道探测数据。首先，读取双波长激光雷达的原始数据，对波长532nm和波长1064nm的数据进行扣除背景、几何因子修正；然后对修正数据进行检验是否有效，舍弃无效数据，对有效数据进行消光系数反演，消光系数反演使用fernald方法，计算出532nm波长和1064nm波长对应的消光系数后进行相除，得到一组比值序列，在数据库中查找532nm/1064nm消光系数比值的气溶胶有效直径；若该波长和消光系数比值在数据库内不可查，则舍弃，若可查，则输出查找值并通过气溶胶模型进行修正，修正后输出粒子谱值，反演结束。
如图2所示，通过模拟计算，分别得到337nm波长、1060nm波长的消光系数比值、后向散射系数比值与粒子有效直径之间的关系，550nm波长、1060nm波长的消光系数比值、后向散射系数比值与粒子有效直径之间的关系(注：图中各曲线与y轴相交从上到下依次表示337nm波长、1060nm波长的消光系数比值与粒子有效直径之间的关系，550nm波长、1060nm波长的消光系数比值与粒子有效直径之间的关系，337nm波长、1060nm波长的后向散射系数比值与粒子有效直径之间的关系，550nm波长、1060nm波长的后向散射系数比值与粒子有效直径之间的关系)。
在0.1-1um的粒子有效直径范围内，337nm波长、1060nm波长的消光系数比值与粒子有效直径呈线性关系，因此使用355nm、1064nm双波长可反演粒子有效直径。本申请技术方案中，为了提高反演的准确性，双波长选择时应选择消光系数比值与气溶胶粒子有效直径呈线性关系的波长组合。
同时，还需要更多的气溶胶模型，建立数据库来提高反演的准确性，并通过水平探测和气溶胶粒径谱仪的探测结果进行对比验证。
以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不会使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (7)

1.一种双波长粒子谱反演方法，其特征在于：包括以下步骤：
S1、读取双波长激光雷达的原始数据，对波长数据进行预处理后，检验波长数据是否有效；
S2、若波长数据无效，则舍弃无效波长数据，并返回S1，否则对波长数据进行消光系数反演；
S3、计算消光系数比值，并在数据库中查找对应波长内消光系数比值的气溶胶有效直径；
S4、若波长内消光系数比值在数据库中不可查，则舍弃并返回S1，否则输出对应的气溶胶有效直径；
S5、通过气溶胶模型进行修正，修正后输出粒子谱，完成反演。
2.根据权利要求1所述的双波长粒子谱反演方法，其特征在于：S1中所述对波长数据进行预处理，包括：
将截取原始数据中最后一段数据的平均值作为背景值，利用数据序列减去背景值对波长数据进行扣除背景；
通过手动修正原始数据中近地面段数据，对波长数据进行几何因子修正。
3.根据权利要求2所述的双波长粒子谱反演方法，其特征在于：S2中所述对波长数据进行消光系数反演，包括：
求解米散射激光雷达方程，通过假设激光雷达比，用Fernald方法计算得到消光系数。
4.根据权利要求3所述的双波长粒子谱反演方法，其特征在于：所述米散射激光雷达方程如下：
Figure FDA0003212915560000011
式中，P(r)为激光雷达接收的高度r处的大气后向散射回波功率，Pt为激光雷达发射的激光功率，k为激光雷达系统常数，βm(r)为空气分子的后向散射系数，βa(r)为大气气溶胶的后向散射系数，αm(r’)为空气分子的消光系数，αa(r’)为大气气溶胶的消光系数；
大气气溶胶消光系数通过下式计算得到：
Figure FDA0003212915560000021
式中，Sa是大气气溶胶的消光后向散射比，或称激光雷达比，Sm是空气分子的消光后向散射比。
5.根据权利要求4所述的双波长粒子谱反演方法，其特征在于：S3中所述计算消光系数比值，包括：
对计算出两个波长对应的消光系数相除，得到一组比值序列。
6.根据权利要求5所述的双波长粒子谱反演方法，其特征在于：所述数据库中消光系数比值与气溶胶有效直径之间的对应关系通过模拟计算得到。
7.根据权利要求6所述的双波长粒子谱反演方法，其特征在于：所述气溶胶模型为通过历史探测数据统计得出的不同天气、不同气溶胶类型下粒子有效直径的真值范围。
CN202110936121.7A 2021-08-16 2021-08-16 一种双波长粒子谱反演方法 Pending CN113673139A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110936121.7A CN113673139A (zh) 2021-08-16 2021-08-16 一种双波长粒子谱反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110936121.7A CN113673139A (zh) 2021-08-16 2021-08-16 一种双波长粒子谱反演方法

Publications (1)

Publication Number Publication Date
CN113673139A true CN113673139A (zh) 2021-11-19

Family

ID=78542951

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110936121.7A Pending CN113673139A (zh) 2021-08-16 2021-08-16 一种双波长粒子谱反演方法

Country Status (1)

Country Link
CN (1) CN113673139A (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115113232A (zh) * 2022-05-31 2022-09-27 合肥中科光博量子科技有限公司 一种基于伽马谱分布的气溶胶部分微物理参数反演方法
CN115482643A (zh) * 2022-08-24 2022-12-16 清华大学 火灾烟雾探测器及其探测方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103234877A (zh) * 2013-04-10 2013-08-07 中国科学院合肥物质科学研究院 一种大气颗粒物粒径谱时空分布激光雷达数据反演方法
CN204086139U (zh) * 2014-09-24 2015-01-07 南京信息工程大学 用于区分雾和霾的双波长激光探测系统
CN108535149A (zh) * 2018-04-17 2018-09-14 泰州职业技术学院 一种多波长激光雷达气溶胶粒子谱分布反演方法及系统
CN109596594A (zh) * 2018-11-27 2019-04-09 南京信息工程大学 基于拉曼-米散射激光雷达的气溶胶消光系数反演方法
CN110441777A (zh) * 2019-07-11 2019-11-12 中山大学 一种基于激光雷达的气溶胶垂直廓线的反演方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103234877A (zh) * 2013-04-10 2013-08-07 中国科学院合肥物质科学研究院 一种大气颗粒物粒径谱时空分布激光雷达数据反演方法
CN204086139U (zh) * 2014-09-24 2015-01-07 南京信息工程大学 用于区分雾和霾的双波长激光探测系统
CN108535149A (zh) * 2018-04-17 2018-09-14 泰州职业技术学院 一种多波长激光雷达气溶胶粒子谱分布反演方法及系统
CN109596594A (zh) * 2018-11-27 2019-04-09 南京信息工程大学 基于拉曼-米散射激光雷达的气溶胶消光系数反演方法
CN110441777A (zh) * 2019-07-11 2019-11-12 中山大学 一种基于激光雷达的气溶胶垂直廓线的反演方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115113232A (zh) * 2022-05-31 2022-09-27 合肥中科光博量子科技有限公司 一种基于伽马谱分布的气溶胶部分微物理参数反演方法
CN115482643A (zh) * 2022-08-24 2022-12-16 清华大学 火灾烟雾探测器及其探测方法
CN115482643B (zh) * 2022-08-24 2024-02-02 清华大学 火灾烟雾探测器及其探测方法

Similar Documents

Publication Publication Date Title
CN113673139A (zh) 一种双波长粒子谱反演方法
CN109117549B (zh) 一种污染物扩散预测方法及系统
CN105628649B (zh) 气体现场监测特征吸收峰提取方法
CN107643252B (zh) Wms检测瓶内氧气浓度的实时扣背景非线性校正方法
Ptashnik et al. Water vapor continuum absorption in near-IR atmospheric windows
CN113933258A (zh) 基于走航监测的VOCs污染物溯源方法、终端及系统
CN112684471A (zh) 一种基于激光雷达的气溶胶微物理特性反演方法
CN112766127B (zh) 基于互补集合模态分解和sg滤波的雷云点电荷定位方法
CN112730307B (zh) 气体浓度非线性测量方法、装置、计算机设备及存储介质
CN107132545B (zh) 气溶胶中的颗粒物对遥感成像的质量影响的分析方法
Zhou et al. XCO 2 satellite retrieval experiments in short-wave infrared spectrum and ground-based validation
Zheng et al. McPrA-A new gas profile inversion algorithm for MAX-DOAS and apply to 50 m vertical resolution
KR102414410B1 (ko) 직교부분 공간투영 기법을 활용한 독성 화학 가스 탐지 장치, 탐지 방법 및 이를 실행시키기 위한 프로그램을 기록한 컴퓨터 판독 가능한 기록 매체
CN116067338A (zh) 基于激光点云数据的地表沉降监测方法及系统
CN108445508B (zh) 一种轻度雾霾时大气消光系数的波段转换方法及系统
Sipeng et al. Correction of the fernald method using real-time average lidar ratios with Mie–Rayleigh–Raman lidar
CN115032655B (zh) 基于可见光和紫外光的双波长激光雷达反演臭氧浓度算法
Uchiyama Method to retrieve single scattering properties of aerosols using multi-wavelength scattering and absorption coefficient data measured by integrating nephelometer and absorption photometer
Pliutau et al. Simulation studies for comparative evaluation of alternative spectral regions for the sensing of CO2 and O2 suitable for the ASCENDS Mission
CN117590360B (zh) 一种多波长激光雷达颜色比标定方法
CN116380811B (zh) 一种基于氧气吸收带的云检测方法
CN116577324B (zh) 一种对工业区污染源扫描成像的超光谱遥感系统
CN116503390B (zh) 一种基于计算机视觉的五金零件缺陷检测方法
CN118032718A (zh) 基于太阳散射光的大气臭氧整层分布廓线反演方法
US20240295659A1 (en) Method, device, apparatus, and storage medium for identifying type of atmospheric aerosol

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