CN114136320B - 基于特征互补特性的海洋多物理场参数特征定位融合方法 - Google Patents

基于特征互补特性的海洋多物理场参数特征定位融合方法 Download PDF

Info

Publication number
CN114136320B
CN114136320B CN202111375208.8A CN202111375208A CN114136320B CN 114136320 B CN114136320 B CN 114136320B CN 202111375208 A CN202111375208 A CN 202111375208A CN 114136320 B CN114136320 B CN 114136320B
Authority
CN
China
Prior art keywords
parameter
characteristic
parameters
field
physical
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
CN202111375208.8A
Other languages
English (en)
Other versions
CN114136320A (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.)
707th Research Institute of CSIC
Original Assignee
707th Research Institute of CSIC
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 707th Research Institute of CSIC filed Critical 707th Research Institute of CSIC
Priority to CN202111375208.8A priority Critical patent/CN114136320B/zh
Publication of CN114136320A publication Critical patent/CN114136320A/zh
Application granted granted Critical
Publication of CN114136320B publication Critical patent/CN114136320B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01CMEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C21/00Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00
    • G01C21/20Instruments for performing navigational calculations
    • G01C21/203Specially adapted for sailing ships
    • 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
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Automation & Control Theory (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种基于特征互补特性的海洋多物理场参数特征定位融合方法,步骤为:1)在n个海洋物理场参数中,根据参数特征误差带选择出所有具备特征定位融合条件的参数;2):m个参数进行特征定位融合,获取海洋多物理场参数匹配位置lo;3):若步骤1)确定的m个具备特征定位融合条件参数中有k个参数的特征误差带具有明显带状分布形式,k>1,且这k个参数特征误差带的交集与R相等,则采用自由方向组合方法获取海洋多物理场参数匹配位置lo。本发明充分利用海洋重力场、海洋地磁场、海底水深等多物理场参数特征之间存在的互补特性实现海洋多物理场参数特征定位融合,与单一参数匹配定位相比具有更好的定位精度、可靠性和适用性。

Description

基于特征互补特性的海洋多物理场参数特征定位融合方法
技术领域
本发明属于深远海平台水下自主导航技术领域,特别涉及一种基于特征互补特性的海洋海洋多物理场参数特征定位融合方法。
背景技术
海洋物理场包括海洋重力场、海底地形、海洋地磁场,海洋物理场与地理位置具有强相关性,是天然的导航参考资源。海洋物理场导航是水下导航体系建设的重要增量手段,是未来水下导航技术发展的一个重要方向。
地球重力场分布于地球表面及其临近空间,海底地形起伏,地质构造和矿产分布导致的海底密度分布不均匀均可引起重力场异常。地球和周围空间之间存在磁场,海洋中广泛存在的岩浆岩、火山岩、变质岩的分布及其几何展布会引起海洋地磁场分布异常。海底地形是指海水覆盖之下的固体地球表面形态,大洋中地形复杂,水深变化多端,包括绵长的海岭及横向断裂带、高耸的海山、起伏的海丘、深邃的海沟等等。
综上海洋重力场、海洋地磁场和海底地形均具备特征定位的条件,但受限于测量手段发展水平,在可用性和精度方面有所不足。由于海洋重力场、海洋地磁场和海底地形的成因有明显区别,其特征形态以及定位能力也必然有所不同。本文提出了一种利用区域海洋重力场、地磁场和海底地形特征互补特性的海洋多物理场参数特征定位融合方法。
发明内容
本发明的目的在于克服单一海洋多物理场参数特征匹配定位技术的不足,充分利用海洋重力场、海洋地磁场、海底水深等多物理场参数特征之间存在的互补特性,提出一种海洋多物理场参数特征定位融合方法,与单一参数匹配定位相比具有更好的定位精度、可靠性和适用性。
本发明解决其技术问题是通过以下技术方案实现的:
一种基于特征互补特性的海洋多物理场参数特征定位融合方法,其特征在于:其包括如下步骤:
1)在n个海洋物理场参数中,根据参数特征误差带选择出所有具备特征定位融合条件的参数:
①载体沿某个方向航行进入特征区域,首先计算n个参数在该特征区匹配定位的特征误差带:计算区域中心处航迹背景图值与周围各点处航迹背景图值差值标准差,差值标准差小于k倍该参数测量误差的位置点组成的集合即为该参数的特征误差带,k取1.5-2,计算这n个参数在该特征区域的特征误差带,记为Ri(i=1,…,n),R为n个参数特征误差带的交集,即R=R1∩…∩Rn
②对于其中参数X,RX为除参数X外其它n-1个参数的特征误差带的交集,若R是RX的真子集,则参数X具备特征定位融合的条件,选择出所有具备特征定位融合条件的参数,假设共m个,记为Xi,i=1,2,…,m,假设R中有h个位置点,记为rj,j=1,2,…,h;
2):m个参数进行特征定位融合,获取海洋多物理场参数匹配位置lo:
对于参数Xi,载体沿某个方向航行进入特征区域,通过该参数测量传感器获取航迹上测量信息Xi,计算Xi与R中所有位置点rj,j=1,2,…,h处航迹背景图信息Yj差值的标准差,记为δij;参数Xi的测量误差水平记为σi;对于位置点rj,计算valj=(δ1j12j2+…+δmjm)/m,R中使valj最小的位置点即为海洋多物理场参数匹配位置lo;
3):若步骤1)确定的m个具备特征定位融合条件参数中有k个参数的特征误差带具有明显带状分布形式,k>1,且这k个参数特征误差带的交集与R相等,则采用自由方向组合方法获取海洋多物理场参数匹配位置lo,k个参数记为Zi,i=1,…,k。
而且,所述海洋物理场参数包括海洋重力场、地磁场的各种参数以及海底水深参数等。
而且,所述步骤3)采用自由方向组合方法确定匹配位置lo的步骤为:
①确定参数Zi,i=1,2,…,k的定位误差最小方向gi和自由方向ti,参数Zi的自由方向ti为该参数特征误差带带状延伸的方向,定位误差最小方向gi为自由方向ti的垂向。
②载体沿某个方向航行进入特征区域,根据参数Zi的航迹测量信息和背景图信息获取其特征匹配位置wi,i=1,…,k;
③根据Zi参数特征匹配位置wi,在平面坐标系中标出对应点Wi,i=1,…,k,过点Wi沿其自由方向ti作直线,共作k条直线,形成d个交点;连接所有可能的交点形成凸多边形,所有交点或是凸多边形顶点,或在凸多边形内部,凸多边形的内切圆心O位置即为匹配位置lo,或当k条直线相交于1点时,该点位置即为匹配位置lo。
本发明的优点和有益效果为:
本发明的基于特征互补特性的海洋多物理场参数特征定位融合方法,充分利用海洋多物理场参数特征之间存在的互补特性,有效将其定位优势融合,获取比单一参数特征定位精度更高的匹配位置,提升了海洋物理场特征匹配定位的精度和效率。
附图说明
图1为本发明的k=3时自由方向组合法搜索匹配位置示意图;
图2a为磁总场异常背景图三维图;
图2b为水深背景图三维图;
图2c为重力异常背景图三维图;
图3a为磁总场异常背景图等值线与待定位区域示意图;
图3b为水深背景图等值线与待定位区域示意图;
图4为磁总场异常、水深、重力异常特征误差带对比图;
图5a为仿真的水深测量信息一;
图5b为仿真的磁总场异常测量信息一;
图6a为仿真的水深测量信息二;
图6b为仿真的磁总场异常测量信息二;
图7a为仿真的水深测量信息三;
图7b为仿真的磁总场异常测量信息三;
图8为第一次磁总场异常、水深特征匹配位置和多物理场参数定位融合位置图;
图9为第二次磁总场异常、水深特征匹配位置和多物理场参数定位融合位置图;
图10为第三次磁总场异常、水深特征匹配位置和多物理场参数定位融合位置图;
图11a为实施例二区域磁总场异常背景图三维图;
图11b为实施例二区域水深背景图三维图;
图11c为实施例二区域重力异常背景图三维图;
图12a为实施例二磁总场异常背景图等值线与待定位区域;
图12b为实施例二水深背景图等值线与待定位区域;
图13为实施例二磁总场异常、水深、重力异常特征误差带对比图;
图14a为实施例二仿真的水深测量信息一;
图14b为实施例二仿真的磁总场异常测量信息一;
图15a为实施例二仿真的水深测量信息二;
图15b为实施例二仿真的磁总场异常测量信息二;
图16a为实施例二仿真的水深测量信息三;
图16b为实施例二仿真的磁总场异常测量信息三;
图17为第一次磁总场异常、水深特征匹配位置和自由组合方法位置图;
图18为第二次磁总场异常、水深特征匹配位置和自由组合方法位置图;
图19为第三次磁总场异常、水深特征匹配位置和自由组合方法位置图。
具体实施方式
下面通过具体实施例对本发明作进一步详述,以下实施例只是描述性的,不是限定性的,不能以此限定本发明的保护范围。
实施例1
一种基于特征互补特性的海洋多物理场参数特征定位融合方法,包括如下步骤:
1)在n个海洋物理场参数中,根据参数特征误差带选择出所有具备特征定位融合条件的参数,海洋物理场参数包括海洋重力场、地磁场的各种参数以及海底水深参数等。本实施例中采用在地磁总场异常、水深、重力异常三个海洋物理场参数中,根据参数特征误差带选择出所有具备特征定位融合条件的参数,该三个海洋物理场参数的背景图如图2a、图2b、图2c所示;
①载体沿某个方向航行进入特征区域F,如图3a,图3b所示,首先计算3个参数在该特征区匹配定位的特征误差带,参数在该特征区域的特征误差带计算方法是:计算区域中心处航迹背景图值与周围各点处航迹背景图值差值标准差,差值标准差小于k倍该参数测量误差的位置点组成的集合即为该参数的特征误差带,k取1.5-2,本实施例取1.5,地磁总场异常、水深、重力异常特征误差带计算结果见图4。计算这3个参数在该特征区域的特征误差带,记为Ri,i=1,2,3,R为3个参数特征误差带的交集,即R=R1∩…∩Rn
②对于其中参数X,RX为除参数X外其它n-1个参数的特征误差带的交集,若R是RX的真子集,则参数X具备特征定位融合的条件,选择出所有具备特征定位融合条件的参数,假设共m个,记为Xi,i=1,2,…,m,假设R中有h个位置点,记为rj,j=1,2,…,h。本实施例中,对于磁总场异常和水深参数,R是除该参数以外的其它2个参数的特征误差带交集的真子集,磁总场异常和水深参数具备特征定位融合的条件记为X1,X2。R中有h个位置点,记为rj,j=1,2,…,h。
2):m个参数进行特征定位融合,获取海洋多物理场参数匹配位置lo:
对于参数Xi,载体沿某个方向航行进入特征区域,通过该参数测量传感器获取航迹上测量信息Xi,计算Xi与R中所有位置点rj,j=1,2,…,h处航迹背景图信息Yj差值的标准差,记为δij;参数Xi的测量误差水平记为σi;对于位置点rj,计算valj=(δ1j12j2+…+δmjm)/m,R中使valj最小的位置点即为海洋多物理场参数匹配位置lo;
仿真了3种误差情况的磁总场异常和水深参数测量信息,如图5a、图5b、图6a、图6b、图7a、图7b所示。获得了3种测量信息对应的海洋多物理场参数匹配位置lo,以及为了对比,同时获得了磁总场异常特征匹配位置Pm和水深特征匹配位置Pd,如图8、图9、图10所示。磁总场异常特征匹配位置、水深特征匹配位置和多物理场参数定位融合位置误差如表1所示,多物理场参数定位融合位置误差小于磁总场异常特征匹配位置和水深特征匹配位置误差,验证了多物理场参数定位融合方法的有效性。
表1 3次磁总场异常、水深特征匹配位置和多物理场参数定位融合位置误差
Figure BDA0003363688940000051
实施例2
步骤1)、2)同于实施例1,本实施例三个海洋物理场参数的背景图如图11a、图11b、图11c所示。
步骤3)为:若步骤1)确定的m个具备特征定位融合条件参数中有k个参数的特征误差带具有明显带状分布形式,k>1,且这k个参数特征误差带的交集与R相等,则采用自由方向组合方法获取海洋多物理场参数匹配位置lo,k个参数记为Zi,i=1,…,k。
本实施例中步骤1)确定的2个具备特征定位融合条件参数中有2个参数的特征误差带具有明显带状分布形式,且这2个参数特征误差带的交集与R相等,则可以采用一种更为简单直观的方法即自由方向组合方法获取海洋多物理场参数匹配位置lo。2个参数记为Zi,i=1,2,其中Z1为水深参数,Z2为磁总场异常参数。
采用自由方向组合方法确定匹配位置lo的步骤为:
①确定参数Z1、Z2的定位误差最小方向gi和自由方向ti,参数Zi的自由方向ti为该参数特征误差带带状延伸的方向,定位误差最小方向gi为自由方向ti的垂向,如图13所示。
②载体沿某个方向航行进入特征区域,如图12a、图12b所示,根据参数Z1、Z2的航迹测量信息和背景图信息获取其特征匹配位置w1,w2,仿真了3种误差情况的磁总场异常和水深参数测量信息,如图14a、图14b、图15a、图15b、图16a、图16b所示,分别获得了每种情况对应的磁总场异常和水深参数特征匹配位置。
③根据Zi参数特征匹配位置wi,在平面坐标系中标出对应点Wi,i=1,…,k,过点Wi沿其自由方向ti作直线,共作k条直线,形成d个交点;连接所有可能的交点形成凸多边形,所有交点或是凸多边形顶点,或在凸多边形内部,凸多边形的内切圆心O位置即为匹配位置lo,或当k条直线相交于1点时,该点位置即为匹配位置lo。
③根据Zi参数特征匹配位置w1,w2,在平面坐标系中标出对应点W1,W2。过点W1沿其自由方向t1作直线,过点W2沿其自由方向t2作直线,共作2条直线,交于一点,即为匹配位置lo。3种测量误差情况对应的自由方向组合方法如图17、18、19所示。磁总场异常特征匹配位置、水深特征匹配位置和自由方向组合方法位置误差如表2所示,自由方向组合方法位置误差小于磁总场异常特征匹配位置和水深特征匹配位置误差,验证了自由方向组合方法的有效性。
表2 3次仿真水深特征、磁总场异常匹配位置和自由方向组合位置误差
Figure BDA0003363688940000071
尽管为说明目的公开的本发明的实施例和附图,但是本领域的技术人员可以理解,在不脱离本发明及所附权利要求的精神和范围内,各种替换、变化和修改都是可能的,因此本发明的范围不局限于实施例和附图所公开的内容。

Claims (2)

1.一种基于特征互补特性的海洋多物理场参数特征定位融合方法,其特征在于:其包括如下步骤:
1)在n个海洋物理场参数中,根据参数特征误差带选择出所有具备特征定位融合条件的参数:①载体沿某个方向航行进入特征区域,首先计算n个参数在该特征区匹配定位的特征误差带:计算区域中心处航迹背景图值与周围各点处航迹背景图值差值标准差,差值标准差小于k倍该参数测量误差的位置点组成的集合即为该参数的特征误差带,k取1.5-2,计算这n个参数在该特征区域的特征误差带,记为Ri(i=1,…,n),R为n个参数特征误差带的交集,即R=R1∩…∩Rn;②对于其中参数X,RX为除参数X外其它n-1个参数的特征误差带的交集,若R是RX的真子集,则参数X具备特征定位融合的条件,选择出所有具备特征定位融合条件的参数,假设共m个,记为Xi,i=1,2,…,m,假设R中有h个位置点,记为rj,j=1,2,…,h;2):m个参数进行特征定位融合,获取海洋多物理场参数匹配位置lo:对于参数Xi,载体沿某个方向航行进入特征区域,通过该参数测量传感器获取航迹上测量信息的参数Xi,计算Xi与R中所有位置点rj,j=1,2,…,h处航迹背景图信息Yj差值的标准差,记为δij;参数Xi的测量误差水平记为σi;对于位置点rj,计算valj=(δ1j12j2+…+δmjm)/m,R中使valj最小的位置点即为海洋多物理场参数匹配位置lo;3):若步骤1)确定的m个具备特征定位融合条件参数中有k个参数的特征误差带具有明显带状分布形式,k>1,且这k个参数特征误差带的交集与R相等,则采用自由方向组合方法获取海洋多物理场参数匹配位置lo,k个参数记为Zi
Figure QLYQS_1
所述步骤3)采用自由方向组合方法确定匹配位置lo的步骤为:①确定参数Zi,i=1,2,…,k的定位误差最小方向gi和自由方向ti,参数Zi的自由方向/>
Figure QLYQS_2
为该参数特征误差带带状延伸的方向,定位误差最小方向gi为自由方向ti的垂向;②载体沿某个方向航行进入特征区域,根据参数Zi的航迹测量信息和背景图信息获取其特征匹配位置wi,i=1,…k;③根据Zi参数特征匹配位置wi,在平面坐标系中标出对应点Wi,i=1,…k,过点Wi沿其自由方向ti作直线,共作k条直线,形成d个交点;连接所有可能的交点形成凸多边形,所有交点或是凸多边形顶点,或在凸多边形内部,凸多边形的内切圆心O位置即为匹配位置lo,或当k条直线相交于1点时,该点位置即为匹配位置lo。
2.根据权利要求1所述基于特征互补特性的海洋多物理场参数特征定位融合方法,其特征在于:所述海洋物理场参数包括海洋重力场、地磁场的各种参数以及海底水深参数。
CN202111375208.8A 2021-11-19 2021-11-19 基于特征互补特性的海洋多物理场参数特征定位融合方法 Active CN114136320B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111375208.8A CN114136320B (zh) 2021-11-19 2021-11-19 基于特征互补特性的海洋多物理场参数特征定位融合方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111375208.8A CN114136320B (zh) 2021-11-19 2021-11-19 基于特征互补特性的海洋多物理场参数特征定位融合方法

Publications (2)

Publication Number Publication Date
CN114136320A CN114136320A (zh) 2022-03-04
CN114136320B true CN114136320B (zh) 2023-06-20

Family

ID=80390277

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111375208.8A Active CN114136320B (zh) 2021-11-19 2021-11-19 基于特征互补特性的海洋多物理场参数特征定位融合方法

Country Status (1)

Country Link
CN (1) CN114136320B (zh)

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4954837A (en) * 1989-07-20 1990-09-04 Harris Corporation Terrain aided passive range estimation
CN103344954B (zh) * 2013-07-08 2014-07-02 国家海洋局第二海洋研究所 一种基于多来源水深数据融合的海底地形地貌构建方法
CN105043388B (zh) * 2015-06-29 2018-01-05 中国船舶重工集团公司第七0七研究所 基于惯性/重力匹配组合导航的向量搜索迭代匹配方法
CN108225310B (zh) * 2017-12-22 2020-08-25 中国船舶重工集团公司第七0七研究所 一种重力辅助导航航迹规划方法
CN108896040B (zh) * 2018-03-30 2021-03-26 中国空间技术研究院 天空海一体化水下潜器惯性/重力组合导航方法和系统
CN109141426B (zh) * 2018-08-10 2020-11-10 中国空间技术研究院 一种水下重力匹配导航适配区的方法
CN109059964B (zh) * 2018-09-19 2021-07-23 中国船舶重工集团公司第七0七研究所 一种基于“重力峰”的惯性导航与重力测量双校准方法
CN112229404B (zh) * 2020-08-31 2024-02-09 中国空间技术研究院 基于海底地形三维优化原理提高海洋重力场插值精度方法
CN112762924B (zh) * 2020-12-23 2023-07-14 北京机电工程研究所 基于重力梯度-地形异源数据匹配的导航定位方法
CN112729292A (zh) * 2021-01-08 2021-04-30 中国船舶重工集团公司第七0七研究所 基于层次分析方法的多重统计参数的重力适配区选取方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
重力匹配定位技术要素影响分析与误差建模研究;李晓平;《重力测量技术——2014年惯性技术发展动态发展方向研讨会文集》;全文 *

Also Published As

Publication number Publication date
CN114136320A (zh) 2022-03-04

Similar Documents

Publication Publication Date Title
CN106886024B (zh) 深海多波束声线精确跟踪方法
Masoud et al. Tectonic architecture through Landsat-7 ETM+/SRTM DEM-derived lineaments and relationship to the hydrogeologic setting in Siwa region, NW Egypt
CN106444776B (zh) 一种无人艇自主性能评估方法
CN110260044A (zh) 一种海底管道定位方法
CN107367735A (zh) 一种在役油气管道河流穿越段电磁‑声波检测方法
CN114136320B (zh) 基于特征互补特性的海洋多物理场参数特征定位融合方法
CN112378376B (zh) 一种基于传感阵列和测斜仪的海底变形联合监测方法
Tomczak Modern methods of underwater positioning applied in subsea mining
Smith A geographical primer to maritime boundary‐making
Peñaherrera-Palma et al. Rapid seafloor mapping of the northern Galapagos Islands, Darwin and Wolf
Murphy et al. A GIS-based application of drainage basin analysis and geomorphometry in the submarine environment: the Gollum Canyon System, Northeast Atlantic
US20170254915A1 (en) Earth surveying with two or more marine vehicles for improved drilling applications
CN107255810A (zh) 基于单信标测距定位双精度加权融合的航向误差补偿方法
CN111768068A (zh) 一种海上搜救船舶响应能力评估方法
CN114136321B (zh) 基于同源特征的海洋多物理场信息融合定位方法
Shin et al. A study on airborne LiDAR calibration and operation techniques for bathymetric survey
Song et al. Underwater slope measurement using a tilted multibeam sonar head
Wang et al. The Application Of Acoustic Detection Technology In The Investigation Of Submarine Pipelines
CN112444890B (zh) 一种深部铀矿二维氡气差量探测方法
WO2021192082A1 (ja) 位置誤差予測装置、モデル作成装置、位置誤差予測方法、モデル作成方法及びプログラム
Werdaya et al. The Geomorphology of Bintuni Bay, Eastern Indonesia
Roh et al. DETERMINATION OF VERTICAL DATUM LEVEL FOR TIDAL BENCH MARK USING GNSS BUOY OBSERVATIONS
Li et al. A refined localization method for underwater tetherless sensor networks
Ånonsen et al. Water Referenced Doppler Velocity Aiding for AUV Navigation
Ekpa et al. Determination of a Section of Woji Riverbed Depths for Safe Navigation

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