CN113009531A - 一种小尺度高精度的低空对流层大气折射率模型 - Google Patents
一种小尺度高精度的低空对流层大气折射率模型 Download PDFInfo
- Publication number
- CN113009531A CN113009531A CN202110198298.1A CN202110198298A CN113009531A CN 113009531 A CN113009531 A CN 113009531A CN 202110198298 A CN202110198298 A CN 202110198298A CN 113009531 A CN113009531 A CN 113009531A
- Authority
- CN
- China
- Prior art keywords
- atmospheric
- grid
- refractive index
- constraint
- elevation
- 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
Links
Images
Classifications
-
- 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
- G01S19/00—Satellite radio beacon positioning systems; Determining position, velocity or attitude using signals transmitted by such systems
- G01S19/01—Satellite radio beacon positioning systems transmitting time-stamped messages, e.g. GPS [Global Positioning System], GLONASS [Global Orbiting Navigation Satellite System] or GALILEO
- G01S19/13—Receivers
- G01S19/35—Constructional details or hardware or software details of the signal processing chain
- G01S19/37—Hardware or software details of the signal processing chain
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
Landscapes
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Biochemistry (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Health & Medical Sciences (AREA)
- Computer Networks & Wireless Communication (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种小尺度高精度的低空对流层大气折射率模型。它利用研究区域内的GNSS观测数据及配置的气象产品作为输入值,然后将研究区域低对流层划分若干个独立的三维网格;并假设每个独立网格内的大气折射率是均匀不变的,将两个测站间的斜路径对流层延迟值表示为单位网格内的大气折射率与斜路径长度的乘积;通过附加水平约束、垂直约束及先验值约束观测方程,得到最终的层析观测模型。本发明克服了现有技术对流层经验模型已经不能满足精密定位的需求的缺点;具有能实时测量层析网格内的大气折射率,测量精度较高的优点。
Description
技术领域
本发明涉及GNSS高精度定位及GNSS气象学技术领域,更具体地说它是一种小尺度高精度的低空对流层大气折射率模型,特别涉及到GNSS网高精度实时快速解算及极端天气研究。
背景技术
对流层是大气中离我们最近的一层,也是最稠密的一层。虽然对流层对整个大气圈而言只是很薄的一层,但它集中了大气质量的80%以上。大气中的水汽几乎都集中在这里,云、雨、雪等天气现象都发生在这一层,对流层与人类的关系最为密切。对对流程建模分析在GNSS数据处理中是提高定位精度的重要课题,也是地球科学、大地测量学应用和气象应用的课题。
对流层大气传播延迟是影响GNSS测量精度的主要误差之一。由于真实对流层大气的复杂性,对流层延迟的研究变得很困难。随着GNSS定位精度要求的提高,对流层经验模型已经不能满足精密定位的需求。因此,如何建立一个更加准确且符合区域特性的对流层延迟改正模型成为一个当前研究的热点问题。
发明内容
本发明的目的是为了提供一种小尺度高精度的低空对流层大气折射率模型,能实时测量层析网格内的大气折射率,测量精度较高。
为了实现上述目的,本发明的技术方案为:一种小尺度高精度的低空对流层大气折射率模型,其特征在于:利用研究区域内的GNSS观测数据及配置的气象产品作为输入值,然后将研究区域低对流层划分若干个独立的三维网格;并假设每个独立网格内的大气折射率是均匀不变的,将两个测站间的斜路径对流层延迟值表示为单位网格内的大气折射率与斜路径长度的乘积;通过附加水平约束、垂直约束及先验值约束观测方程,得到最终的层析观测模型。
在上述技术方案中,研究区域低对流层网格划分,具体方法为:
垂直网格划分包括均匀网格划分和非均匀网格划分;水平方向上网格划分,保持不同的GNSS测站分布在不同的网格内,水平距离保持5-8km。
在上述技术方案中,两个测站间的斜路径对流层延迟值的具体计算方法如下:
式(1)中:h1和h2分别为两个测站高度,N表示大气折射率,h表示高程,α表示两测站与水平线的夹角;
其中,N通过式(2)计算得到:
式(2)中:P表示大气压强,e表示水汽压,T表示温度,k1,k2,k3是常数;
联合公式(1)和(2)进一步得到STD:
式(3)中:h1和h2分别为两个测站高度,h表示高程,α表示两测站与水平线的夹角,P表示大气压强,e表示水汽压,T表示温度,k1,k2,k3是常数;
大气压强与大气温度,水汽压与大气压强在垂直方向上存在的关系以及大气温度随高程变化分布表示为:
T=T1+β(h-h1) (6)
式(4)、(5)、(6)中:P,T,e分别为大气压强,大气温度和水汽压;P1,T1,e1分别表示为地表处大气压强,大气温度和水汽压;h表示高程;β表示大气温度递减率;h1表示地表处的高程;γ表示大气混合比指数;g1表示重力加速度;Rd表示干大气气体常数;
把公式(4)、(5)和(6)带入公式(3),进一步得到STD:
式(7)中:h1和h2分别为两个测站高度,h表示高程,α表示两测站与水平线的夹角,k1,k2,k3是常数;T为大气温度;P1,T1,e1分别表示为地表处大气压强,大气温度和水汽压;h表示高程;β表示大气温度递减率;h1表示地表处的高程;γ表示大气混合比指数;g0表示重力加速度;Rd表示干大气气体常数;
在上述技术方案中,设研究区域内所划分的经度方向的网格数为nl,纬度方向网格数为nn,高度方向的网格数为nh,每个网格中的独立的大气折射率值为Ni,j,k,其单第q条信号所穿越的网格的长度为其单位为km,则第q条观测信号所满足的观测方程为:
式(8)中:STDq为两个GNSS测站间的斜路径对流层延迟,i,j,k表示网格的三维坐标。
在上述技术方案中,在水平方向,利用高斯距离加权函数来获取水平约束系数B,高斯距离加权函数的具体表达式如下:
式(9)中:下标i,j,k表示三维网格的坐标;nl和nn分别为在东西方向和南北方向上层析网格划分的个数;di,j,k表示未知水汽密度网格和已知水汽密度网格之间的距离;δ为平滑因子,是一个常数。
在上述技术方案中,通过对探空以及掩星资料的研究,大气折射率在空间上呈现出类似高斯指数变化的特性:
式(10)中:N(h)表示高程h处的大气折射率;hz表示大气折射率标高;N1是一个常数;h1为测站高度;
根据公式(10),垂直方向上相邻两层的大气折射率之间的关系表示为:
式(11)中:下标i,j,k表示网格在东西、南北和高程方向的坐标;hk表示第k层析网格的高度;hk+1表示第k+1层析网格的高度;hz表示大气折射率标高;h1为测站高度。
在上述技术方案中,通过附加水平约束、垂直约束及先验值约束观测方程,得到最终的层析观测模型;将层析模型、水平约束及垂直约束用矩阵的形式表示,最终的层析观测模型表示为:
式(12)中:STD表示斜路径延迟;N表示大气折射率,S表示GNSS信号在网格内穿越的距离;Chori和Cvert分别为水平约束和垂直约束系数;ΔSTD表示斜路径延迟误差向量;Δhori和Δvert分别为水平约束误差向量和垂直约束误差向量。
本发明中的小尺度是指水平横跨距离100km以内;本发明中的低空对流层是指距海平面距离3km以下的高度;本发明的大气折射率偏差优于5N,精度高。
本发明具有如下优点:
(1)本发明新模型利用研究区域内的GNSS观测数据及配置的气象产品作为输入值,然后将研究区域划分若干个独立的三维网格,并假设每个独立网格内的大气折射率是均匀不变的,将两个测站间的斜路径对流层延迟值可以表示为单位网格内的大气折射率与斜路径长度的乘积,再通过附加水平约束、垂直约束及先验值约束;本发明所述的新的低对流层模型的优势是:为所研究的区域提供实时高精度的大气折射率三维分布信息,为GNSS基线的对流层估计提供高精度的先验信息值,同时也可以为小尺度GNSS气象学提供初始值;
(2)本发明利用卡尔曼滤波方法可以估算出每个独立网格内的大气折射率。
附图说明
图1为本发明的工艺流程图。
具体实施方式
下面结合附图详细说明本发明的实施情况,但它们并不构成对本发明的限定,仅作举例而已。同时通过说明使本发明的优点更加清楚和容易理解。
本发明基于小区域内的气象资料及GNSS观测数据建立的对流层参考模型,它能较全球模型更能可靠地反映局部大气实际特性。本发明首先利用积分的方法获取两个测站间的斜路径对流层延迟,将研究区域划分成若干个独立的网格,利用层析反演的方法获取独立网格内的大气折射率。
参阅附图可知:一种小尺度高精度的低空对流层大气折射率模型,利用研究区域内的GNSS观测数据及配置的气象产品作为输入值(即取两个测站之间STD之差作为输入值),然后将研究区域低对流层划分若干个独立的三维网格;并假设每个独立网格内的大气折射率是均匀不变的,将两个测站间的斜路径对流层延迟值表示为单位网格内的大气折射率与斜路径长度的乘积;通过附加水平约束、垂直约束及先验值约束观测方程,得到最终的层析观测模型(即低空对流层大气折射率模型)。
本发明所述的低空对流层大气折射率模型利用卡尔曼滤波方法对层析观测模型进行解算,便可以得到层析网格内的大气折射率(即每个独立网格内的大气折射率)。
原始观测的GNSS数据,需要配置轨道和钟差文件,才能进行精密单点定位解算,获得ZTD、测站坐标、模糊度及接收机钟差。把接收机钟差、模糊度及测站坐标带入精密单点定位的公式,可以获得每颗卫星对应的STD。
进一步地,研究区域低对流层网格划分,具体方法为:
垂直方向上,网格底层的高度为地表,网格顶层为区域内GNSS高程最高的高度;
垂直网格划分,可以分为包括均匀网格划分和非均匀网格划分;水平方向上网格划分,尽量保持不同的GNSS测站分布在不同的网格内,水平距离保持5-8km。
进一步地,获得两个测站之间的斜路径对流层延迟STD(Slant troposphericdelay);两测站间的STD利用气象数据和积分模型可以得到,两个测站间的斜路径对流层延迟值的具体计算方法如下:
式(1)中,h1和h2分别为两个测站高度,N表示大气折射率,h表示高程,α表示两测站与水平线的夹角;
其中,N通过式(2)计算得到:
式(2)中,P表示大气压强,e表示水汽压,T表示温度,k1,k2,k3是常数;
联合公式(1)和(2)可进一步得到STD:
式(3)中:h1和h2分别为两个测站高度,h表示高程,α表示两测站与水平线的夹角,P表示大气压强,e表示水汽压,T表示温度,k1,k2,k3是常数;
大气压强与大气温度,水汽压与大气压强在垂直方向上存在的关系以及大气温度随高程变化可以分布表示为:
T=T1+β(h-h1) (6)
其中,P,T,e分别为大气压强,大气温度和水汽压;P1,T1,e1分别表示为地表处大气压强,大气温度和水汽压;h表示高程;β表示大气温度递减率;h1表示地表处的高程;γ表示大气混合比指数;g1表示重力加速度;Rd表示干大气气体常数;
把公式(4)、(5)和(6)带入公式(3),进一步得到STD:
其中,式(7)中:h1和h2分别为两个测站高度,h表示高程,α表示两测站与水平线的夹角,P表示大气压强,e表示水汽压,k1,k2,k3是常数;P,T,e分别为大气压强,大气温度和水汽压;P1,T1,e1分别表示为地表处大气压强,大气温度和水汽压;h表示高程;β表示大气温度递减率;h1表示地表处的高程;γ表示大气混合比指数;g1表示重力加速度;Rd表示干大气气体常数;
进一步地,STD可以表达网格内的大气折射率与穿越网格距离长度乘积的累积和;设研究区域内所划分的经度方向的网格数为nl,纬度方向网格数为nn,高度方向的网格数为nh,每个网格中的独立的大气折射率值为Ni,j,k,其单第q条信号所穿越的网格的长度为其单位为km,则第q条观测信号所满足的观测方程为:
式(8)中:STDq为两个GNSS测站间的斜路径对流层延迟,i,j,k表示网格的三维坐标;其中,STDq表示其中任意一条GNSS信号所获得的STD。
更进一步地,有很多网格内没有信号穿越,故本发明中的层析方程(8)通常是无法直接利用平差的方法进行解算,需要在公式(8)的基础上附加一定的约束条件,才能对公式(8)进行求解;
在水平方向,通常利用高斯距离加权函数来获取水平约束系数B,高斯距离加权函数的具体表达式如下:
式(9)中:下标i,j,k表示三维网格的坐标;nl和nn分别为在东西方向和南北方向上层析网格划分的个数;di,j,k表示未知水汽密度网格和已知水汽密度网格之间的距离;δ为平滑因子,是一个常数。
更进一步地,为了建立垂直方向上信号穿越的网格和未有信号穿越的网格之间的联系,增加层析方程的适定性;本发明通过对探空以及掩星资料的研究,发现了大气折射率在空间上呈现出类似高斯指数变化的特性:
式(10)中:N(h)表示高程h处的大气折射率;hz表示大气折射率标高;N1是一个常数;h1为测站高度;
根据公式(10),垂直方向上相邻两层的大气折射率之间的关系可表示为:
式(11)中:下标i,j,k表示网格在东西、南北和高程方向的坐标;hk表示第k层析网格的高度;hk+1表示第k+1层析网格的高度;hz表示大气折射率标高;h1为测站高度;N表示大气折射率。
更进一步的,本发明利用斜路径对流层延迟反演出不同网格内的大气折射率;
通过附加水平约束、垂直约束及先验值约束观测方程,得到最终的层析观测模型;将层析模型、水平约束及垂直约束用矩阵的形式表示,最终的层析观测模型表示为:
式(12)中:STD表示斜路径延迟;N表示大气折射率;S表示GNSS信号在网格内穿越的距离;Chori和Cvert分别为水平约束和垂直约束系数;ΔSTD表示斜路径延迟误差向量;Δhori和Δvert分别为水平约束误差向量和垂直约束误差向量;
利用卡尔曼滤波方法对层析观测模型进行解算,便可以得到层析网格内的大气折射率(即每个独立网格内的大气折射率)。
验证实例:
现将本发明所述的小尺度高精度的低空对流层大气折射率模型应用于乌东德水电站实时形变监测。与传统的建模方法相比,本发明所述的建模在高程方向上精度提高了70%;另外,使用本发明后,卫星残差由原来的分米级降低到现在的厘米级。
其它未说明的部分均属于现有技术。
Claims (7)
1.一种小尺度高精度的低空对流层大气折射率模型,其特征在于:利用研究区域内的GNSS观测数据及配置的气象产品作为输入值,然后将研究区域低对流层划分若干个独立的三维网格;并假设每个独立网格内的大气折射率是均匀不变的,将两个测站间的斜路径对流层延迟值表示为单位网格内的大气折射率与斜路径长度的乘积;通过附加水平约束、垂直约束及先验值约束观测方程,得到最终的层析观测模型。
2.根据权利要求1所述的小尺度高精度的低空对流层大气折射率模型,其特征在于:研究区域低对流层网格划分,具体方法为:
垂直网格划分包括均匀网格划分和非均匀网格划分;水平方向上网格划分,保持不同的GNSS测站分布在不同的网格内,水平距离保持5-8km。
3.根据权利要求2所述的小尺度高精度的低空对流层大气折射率模型,其特征在于:两个测站间的斜路径对流层延迟值的具体计算方法如下:
式(1)中:h1和h2分别为两个测站高度,N表示大气折射率,h表示高程,α表示两测站与水平线的夹角;
其中,N通过式(2)计算得到:
式(2)中:P表示大气压强,e表示水汽压,T表示温度,k1,k2,k3是常数;
联合公式(1)和(2)进一步得到STD:
式(3)中:h1和h2分别为两个测站高度,h表示高程,α表示两测站与水平线的夹角,P表示大气压强,e表示水汽压,T表示温度,k1,k2,k3是常数;
大气压强与大气温度,水汽压与大气压强在垂直方向上存在的关系以及大气温度随高程变化分布表示为:
T=T1+β(h-h1) (6)
式(4)、(5)、(6)中:P,T,e分别为大气压强,大气温度和水汽压;P1,T1,e1分别表示为地表处大气压强,大气温度和水汽压;h表示高程;β表示大气温度递减率;h1表示地表处的高程;γ表示大气混合比指数;g1表示重力加速度;Rd表示干大气气体常数;
把公式(4)、(5)和(6)带入公式(3),进一步得到STD:
式(7)中:h1和h2分别为两个测站高度,h表示高程,α表示两测站与水平线的夹角,k1,k2,k3是常数;T为大气温度;P1,T1,e1分别表示为地表处大气压强,大气温度和水汽压;h表示高程;β表示大气温度递减率;h1表示地表处的高程;γ表示大气混合比指数;go表示重力加速度;Rd表示干大气气体常数;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110198298.1A CN113009531A (zh) | 2021-02-22 | 2021-02-22 | 一种小尺度高精度的低空对流层大气折射率模型 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110198298.1A CN113009531A (zh) | 2021-02-22 | 2021-02-22 | 一种小尺度高精度的低空对流层大气折射率模型 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113009531A true CN113009531A (zh) | 2021-06-22 |
Family
ID=76406296
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110198298.1A Pending CN113009531A (zh) | 2021-02-22 | 2021-02-22 | 一种小尺度高精度的低空对流层大气折射率模型 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113009531A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113534213A (zh) * | 2021-07-26 | 2021-10-22 | 中国电子科技集团公司第五十四研究所 | 一种公里级区域大气相位不一致性高精度建模与修正方法 |
CN114415208A (zh) * | 2021-12-29 | 2022-04-29 | 自然资源部大地测量数据处理中心 | 一种附加外部数据集信息的地基gnss对流层顶探测方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108008416A (zh) * | 2017-12-04 | 2018-05-08 | 武汉大学 | 一种估算斜路经对流层延迟的积分方法 |
CN108345009A (zh) * | 2018-02-08 | 2018-07-31 | 中国石油大学(华东) | 无先验信息约束的gps三维水汽层析方法和装置 |
CN109543353A (zh) * | 2018-12-29 | 2019-03-29 | 广东电网有限责任公司 | 三维水汽反演方法、装置、设备和计算机可读存储介质 |
-
2021
- 2021-02-22 CN CN202110198298.1A patent/CN113009531A/zh active Pending
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108008416A (zh) * | 2017-12-04 | 2018-05-08 | 武汉大学 | 一种估算斜路经对流层延迟的积分方法 |
CN108345009A (zh) * | 2018-02-08 | 2018-07-31 | 中国石油大学(华东) | 无先验信息约束的gps三维水汽层析方法和装置 |
CN109543353A (zh) * | 2018-12-29 | 2019-03-29 | 广东电网有限责任公司 | 三维水汽反演方法、装置、设备和计算机可读存储介质 |
Non-Patent Citations (2)
Title |
---|
万蓉;郑国光;于胜杰;付志康;: "基于观测约束的地基GPS三维水汽层析技术研究" * |
夏朋飞: "联合空基和地基GNSS观测反演大气水汽方法研究" * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113534213A (zh) * | 2021-07-26 | 2021-10-22 | 中国电子科技集团公司第五十四研究所 | 一种公里级区域大气相位不一致性高精度建模与修正方法 |
CN114415208A (zh) * | 2021-12-29 | 2022-04-29 | 自然资源部大地测量数据处理中心 | 一种附加外部数据集信息的地基gnss对流层顶探测方法 |
CN114415208B (zh) * | 2021-12-29 | 2023-04-28 | 自然资源部大地测量数据处理中心 | 一种附加外部数据集信息的地基gnss对流层顶探测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jenkins et al. | Interactions between ice and ocean observed with phase-sensitive radar near an Antarctic ice-shelf grounding line | |
CN110031877B (zh) | 一种基于grnn模型的区域nwp对流层延迟改正方法 | |
CN105069295B (zh) | 基于卡尔曼滤波的卫星以及地面降水测量值同化方法 | |
CN105842692B (zh) | 一种insar测量中的大气校正方法 | |
CN109917424B (zh) | 多因子约束下的nwp反演对流层延迟的残差改正方法 | |
CN112034490B (zh) | 一种nwp反演对流层延迟改进方法 | |
CN113009531A (zh) | 一种小尺度高精度的低空对流层大气折射率模型 | |
Richter et al. | Ice flow velocities over Vostok Subglacial Lake, East Antarctica, determined by 10 years of GNSS observations | |
Pierdicca et al. | Excess path delays from Sentinel interferometry to improve weather forecasts | |
CN113960635B (zh) | 一种顾及日变化的对流层延迟改正方法 | |
CN113639893B (zh) | 一种基于多气象因子的近地加权平均温度信息获取方法 | |
CN112711022B (zh) | 一种GNSS层析技术辅助的InSAR大气延迟改正方法 | |
CN116858290B (zh) | 基于大型无人机的深远海海面高度观测定标方法和系统 | |
Seko et al. | The meso-γ scale water vapor distribution associated with a thunderstorm calculated from a dense network of GPS receivers | |
Bakış et al. | Analysis and comparison of spatial rainfall distribution applying different interpolation methods in Porsuk river basin, Turkey | |
CN115755103A (zh) | 一种抗差自适应的gnss水汽层析方法 | |
CN103760582A (zh) | 一种遮挡环境下卫星双差观测结构的优化方法 | |
CN116466376A (zh) | 数值预报模式辅助的实时ppp改善方法 | |
CN109164439A (zh) | 一种电波射线上大气折射率计算方法 | |
Zhang et al. | A new four-layer inverse scale height grid model of china for zenith tropospheric delay correction | |
CN107991646A (zh) | 基于云架构的甚低频导航电波传播预测模型精化方法 | |
CN114417565A (zh) | 一种基于地面气象数据的对流层折射率梯度剖面建模方法 | |
Wen et al. | Ice thickness over the southern limit of the Amery Ice Shelf, East Antarctica, and reassessment of the mass balance of the central portion of the Lambert Glacier-Amery Ice Shelf system | |
Jensen et al. | The effect of different tropospheric models on precise point positioning in kinematic mode | |
Farah | Assessment of UNB3M neutral atmosphere model and EGNOS model for near-equatorial-tropospheric delay correction |
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 |