CN110907936B - 一种水下三维地形匹配定位导航声呐及导航方法 - Google Patents

一种水下三维地形匹配定位导航声呐及导航方法 Download PDF

Info

Publication number
CN110907936B
CN110907936B CN201911153149.2A CN201911153149A CN110907936B CN 110907936 B CN110907936 B CN 110907936B CN 201911153149 A CN201911153149 A CN 201911153149A CN 110907936 B CN110907936 B CN 110907936B
Authority
CN
China
Prior art keywords
channel
navigation
sonar
underwater
array
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.)
Expired - Fee Related
Application number
CN201911153149.2A
Other languages
English (en)
Other versions
CN110907936A (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.)
Harbin Engineering University
Original Assignee
Harbin Engineering University
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 Harbin Engineering University filed Critical Harbin Engineering University
Priority to CN201911153149.2A priority Critical patent/CN110907936B/zh
Publication of CN110907936A publication Critical patent/CN110907936A/zh
Application granted granted Critical
Publication of CN110907936B publication Critical patent/CN110907936B/zh
Expired - Fee Related 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/10Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration
    • G01C21/12Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning
    • G01C21/16Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation
    • G01C21/165Navigation; Navigational instruments not provided for in groups G01C1/00 - G01C19/00 by using measurements of speed or acceleration executed aboard the object being navigated; Dead reckoning by integrating acceleration or speed, i.e. inertial navigation combined with non-inertial navigation instruments
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/539Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section

Landscapes

  • Engineering & Computer Science (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种水下三维地形匹配定位导航声呐,包括控制系统、数据采集处理系统、多通道接收机、多通道稀疏平行接收线阵、多通道信号源、多通道发射机和多通道发射基阵;控制系统的输出端双向连接数据采集处理系统的输入端,数据采集处理系统的输出端分别连接多通道接收机和多通道信号源的输入端,解决了传统声呐地形探测效率低、信息量少、地形匹配的效率与准确率低、无法修正惯导和多普勒计程仪的累计误差,无法提高地形匹配定位导航算法的鲁棒性的问题。

Description

一种水下三维地形匹配定位导航声呐及导航方法
技术领域
本发明涉及海洋工程技术设备领域,特别是一种水下三维地形匹配定位导航声呐及导航方法。
背景技术
AUV等水下无人平台是未来的应用热点,自主定位导航能力是保证其到达指定地点顺利完成任务的前提。为了保持隐蔽性,其往往不能接收到水面及地上空间成熟应用的卫星定位系统信号。目前主流的水下定位导航方法是将惯导和多普勒计程仪信息相结合,但存在着累计误差;水下地形、磁场及重力场等地球物理属性一般不随时间和气候的变化而变化,也难以伪装和隐藏,因此近年来利用这些地球物理属性的匹配定位导航方法得到了人们的重视,成为了解决水下长航时精确定位导航的有效辅助手段之一。
水下地形匹配定位导航技术是近年来的关注热点,将地形匹配定位导航技术与其他导航方式进行组合可以获得更好的导航性能,其技术原理是由惯导和多普勒计程仪提供地理位置信息作为参考,利用自身配载的地形测量传感器实时获取水下地形,将实时测量地形与已存储的水下参考数字地形图进行匹配,以此来修正惯导和多普勒计程仪的累计误差,获得高精度的定位导航性能。现有的测量水下地形的传感器主要包括单波束测深声呐、多普勒测速声呐和多波束测深声呐等,获取的深度数据分别为单点、四点以及线地形等,水底地形信息量较少;采用的地形匹配定位导航算法也往往存在鲁棒性较差、收敛较慢等不足。这些问题限制了地形匹配定位导航功能在AUV上的工程应用。
发明内容
为解决现有技术中存在的问题,本发明提供了一种水下三维地形匹配定位导航声呐及其导航方法,解决了传统声呐地形探测效率低、信息量少、地形匹配的效率与准确率低、无法修正惯导和多普勒计程仪的累计误差,无法提高地形匹配定位导航算法的鲁棒性的问题。
本发明采用的技术方案是,一种水下三维地形匹配定位导航声呐,包括控制系统、数据采集处理系统、多通道接收机、多通道稀疏平行接收线阵、多通道信号源、多通道发射机和多通道发射基阵;控制系统的输出端双向连接数据采集处理系统的输入端,数据采集处理系统的输出端分别连接多通道接收机和多通道信号源的输入端,多通道信号源的输出端连接多通道发射机的输入端,多通道发射机的输出端连接多通道发射基阵的输入端,多通道接收机的输出端连接多通道稀疏平行接收线阵的输入端。
优选地,多通道稀疏平行接收线阵包括第一均匀间距接收线阵、第二均匀间距接收线阵和第三均匀间距接收线阵;第一均匀间距接收线阵、第二均匀间距接收线阵和第三均匀间距接收线阵由等间距的多个接收阵元组成,其间距由水平观测扇面大小决定;两组相邻的平行阵列间距不等,两个间距的比值服从互质关系。
优选地,一种水下三维地形匹配定位导航声呐的导航方法,包括以下步骤
S1:水下三维地形匹配定位导航声呐数据处理流程;
S2:地形匹配处理流程;
S3:多源信息融合导航流程。
优选地,S1包括以下步骤:
S11:每条接收阵在水平方向上进行波束形成,对垂直于航行方向的发射扇面内的回波进行空间细分接收,得到回波的水平角θ;
S12:对相邻接收阵输出的同号波束回波进行相位差求解处理,利用每个样本对应的两组相位差之间的定量关系估计回波的俯仰角ψ;
S13:根据声波发射和回波到达之间的时间差估计声呐距目标之间的距离R;
S14:根据计算得到的回波的水平角θ、回波的俯仰角ψ和声呐距目标之间的距离R,从而得到目标每个回波样本对应的距离,对于水底的连续回波即能够得到连续的(R、θ、ψ)序列,从而获得连续的三维水下地形;
S15:依据步骤S11至S14同样的方法得到水中障碍目标的三维位置信息,为AUV提供避碰依据。
优选地,S2采用了批相关和非线性滤波相结合的算法,当声呐载体当前定位误差较大时,批相关用以迅速捕获当前地理区域,为非线性滤波连续输出地理位置信息跟踪导航缩小搜索空间,提高搜索效率,进而结合粒子滤波等非线性滤波方法跟踪修正位置、速度等导航信息。
本发明水下三维地形匹配定位导航声呐及其导航方法的有益效果如下:
1.利用多通道二维稀疏接收基阵方案并采用了针对性的信号处理方案,实现了水下三维“面”地形的测量,相比于传统的基于均匀间距的二维接收基阵方案,在复杂度有限增加条件下,将地形测量能力从“线”地形扩展为“面”地形;相比于传统的基于均匀间距的二维接收基阵方案,大大降低了系统复杂度。
2.在惯导和多普勒计程仪的基础上,融合地形匹配定位导航技术,利用批相关算法通过和数据库中存储的先验地形相匹配,即可以对声呐载体当前所在的地理位置进行辅助修正,再结合利用粒子滤波等非线性滤波算法跟踪修正位置、速度等导航信息,以实现声呐载体平台的辅助导航,能够有效辅助修正惯导和多普勒计程仪的累计误差,提高了导航系统的性能。
附图说明
图1为本发明水下三维地形匹配定位导航声呐及其导航方法的结构原理框图。
图2为本发明水下三维地形匹配定位导航声呐及其导航方法的硬件结构原理框图。
图3为本发明水下三维地形匹配定位导航声呐及其导航方法的多通道三条稀疏平行接收线阵组成方案。
图4为本发明水下三维地形匹配定位导航声呐及其导航方法的软件部分之基于多通道二维稀疏接收基阵方案的俯仰角解算流程图。
图5为本发明水下三维地形匹配定位导航声呐及其导航方法的软件部分之基于批相关和粒子滤波相结合的面地形匹配定位导航算法流程图。
图6为本发明水下三维地形匹配定位导航声呐及其导航方法的利用极大似然准则下的批相关进行面地形匹配定位的效果图。
图7为本发明水下三维地形匹配定位导航声呐及其导航方法的利用粒子滤波进行匹配导航的效果图。
图8为本发明水下三维地形匹配定位导航声呐及其导航方法的利用极大似然准则下的批相关和粒子滤波结合进行面地形匹配导航的效果图。
附图标记:11-水下三维地形匹配定位导航声呐数据处理流程、12-地形匹配处理流程、13-多源信息融合导航流程、14-水下三维地形匹配定位导航声呐采集数据、15-实测数据异常值检测与剔除、16-水下地形处理、17-局部水下“面”地形图重构、18-环境全源数据(电子海图,多波束数据,GIS数据库)、19-建立水下参考地形图、20-基于极大似然准则的批相关和粒子滤波双重模式的地形匹配导航、21-惯性导航、22-多普勒测速仪、23-航位推算、24-卡尔曼滤波、25-导航信息输出、31-控制系统、32-数据采集处理系统、33-多通道接收机、34-多通道稀疏平行接收线阵、35-多通道信号源、36-多通道发射机、37-多通道发射基阵、38.接收阵元、39-第一均匀间距接收线阵、40-第二均匀间距接收线阵、41-第三均匀间距接收线阵。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
如图1所示,一种水下三维地形匹配定位导航声呐的总系统结构,包括:
S1:水下三维地形匹配定位导航声呐数据处理流程;
S2:地形匹配处理流程;
S3:多源信息融合导航流程。
S1包括:14-水下三维地形匹配定位导航声呐采集数据,15-实测数据异常值检测与剔除,16-水下地形处理,17-局部水下“面”地形图重构;
S2包括:18-环境全源数据(电子海图,多波束数据,GIS数据库),19-建立水下参考地形图,20-基于极大似然准则的批相关和粒子滤波双重模式的地形匹配导航;
S3包括:21-惯性导航,22-多普勒测速仪,23-航位推算,24-卡尔曼滤波,25-导航信息输出。
如图2所示,一种水下三维地形匹配定位导航声呐,,包括控制系统31、数据采集处理系统32、多通道接收机33、多通道稀疏平行接收线阵34、多通道信号源35、多通道发射机36和多通道发射基阵37;控制系统31的输出端双向连接数据采集处理系统32的输入端,数据采集处理系统32的输出端分别连接多通道接收机33和多通道信号源35的输入端,多通道信号源35的输出端连接多通道发射机36的输入端,多通道发射机36的输出端连接多通道发射基阵37的输入端,多通道接收机33的输出端连接多通道稀疏平行接收线阵34的输入端。
一种三维地形的匹配定位导航声呐,硬件系统由多通道信号源、多通道发射机、多通道接收机、数据采集处理系统、控制系统、多通道发射基阵和多通道稀疏平行接收线阵组成,软件主要包括“面”地形估计和“面”地形匹配辅助定位导航功能。
本发明的工作原理是控制系统接收指令控制整个声呐系统开始工作,多通道发射基阵发射覆盖垂直于航行方向的一定角度扇面的声波,沿航行方向的角度扇面可由多通道信号源和发射机接收到的控制系统传输过来的工作参数相控定向发射覆盖。多通道稀疏平行接收线阵接收水底散射回波,经过多通道接收机的滤波、放大等调理后,传输至数据采集系统进行采集和信号处理,实时解算得到海底“面”地形数据。当声呐载体初始位置误差较大时,采用基于批相关的“面”地形匹配定位导航模式(捕获模式),融合惯导和多普勒计程仪得到的声呐载体当前所在的地理位置,缩小搜索空间后再采用基于粒子滤波的“面”地形匹配导航模式(跟踪模式);而当声呐载体初始位置小初始误差时,融合地理位置信息后可直接进入跟踪模式,最终在滤波收敛达到稳态后实时输出导航位置。
本发明的目的是这样实现的:
1.多通道稀疏平行接收线阵由三条平行的线性阵列组成,每条线性阵列由等间距的多通道阵元组成,间距由水平观测扇面大小决定;两组相邻的平行阵列间距不等,两个间距的比值服从互质关系,如图3所示。
2.利用每条接收线阵数据进行相移波束形成,可以得到水平角θ的估计。如图4所示,利用空间同一位置回波在互质阵上的唯一对应关系可以进行声波信号俯仰角估计。根据1的基阵结构,存在以下关系式:
Figure BDA0002284113430000071
Figure BDA0002284113430000072
ψ是目标俯仰角;
Figure BDA0002284113430000073
是间距最小(d1)的两个平行线阵相同水平角方向波束输出的相位差;
Figure BDA0002284113430000074
是间距最大(d2)的两个平行线阵相同水平角方向波束输出的相位差;λ是声波波长。对于根据实际数据求出的
Figure BDA0002284113430000075
Figure BDA0002284113430000076
可以搜索求出一组最佳的整数组合(N1和N2),代入上式即可求出俯仰角ψ。再根据探测信号发射和接收到的(θ,ψ)方向上的回波信号的时间差即可求出声呐与目标之间的距离R,从而实现了水中目标的三维定位。对于水底散射的连续回波即可依次求解出水平角θ方向上连续的水底“线”地形,不同水平角θ上的“线”地形即组成了“面”地形。
3.利用2中测得的“面”地形采取极大似然等准则的批相关匹配算法,结合惯导和多普勒计程仪提供的导航信息,即可大致缩小声呐载体当前所在的地理区域,在此基础上进一步结合粒子滤波等非线性滤波方法,跟踪修正位置、速度等导航信息,如图5所示。
如图6至图8所示,水下三维地形匹配定位导航声呐一般安装于AUV的艏部,探测视角可指向前下方30°左右。声呐工作频率为150kHz,波长为1cm,3条平行的接收线阵间距d1=1.5cm,d2=2.5cm,二者之比为3:5,满足互质关系要求。每条接收线阵包括48个阵元,间距为半波长0.5cm;发射弧形基阵通道数为5,水平覆盖角度扇面60°,垂直覆盖角度扇面20°,且俯仰角可通过控制5个发射通道所发射波形间的时延和相位关系进行调节。一次探测即可以获得60°×20°角度扇面内的“面”地形。将此“面”地形数据与事先获得的地形数据进行最大似然匹配运算得到比较准确的地理位置,并在此基础上采用粒子滤波,实现声呐载体连续的地形匹配定位,从而完成辅助导航。

Claims (1)

1.一种水下三维地形匹配定位导航声呐的导航方法,其特征在于,包括以下步骤:
S1:水下三维地形匹配定位导航声呐数据处理流程;
S2:地形匹配处理流程;
S3:多源信息融合导航流程;
S1包括以下步骤:
S11:每条接收阵在水平方向上进行波束形成,对垂直于航行方向的发射扇面内的回波进行空间细分接收,得到回波的水平角θ;
S12:对相邻接收阵输出的同号波束回波进行相位差求解处理,利用每个样本对应的两组相位差之间的定量关系估计回波的俯仰角ψ;
S13:根据声波发射和回波到达之间的时间差估计声呐距目标之间的距离R;
S14:根据计算得到的回波的水平角θ、回波的俯仰角ψ和声呐距目标之间的距离R,从而得到目标每个回波样本对应的距离,对于水底的连续回波即能够得到连续的(R、θ、ψ)序列,从而获得连续的三维水下地形;
S15:依据步骤S11至S14同样的方法得到水中障碍目标的三维位置信息,为AUV提供避碰依据;
所述S2采用了批相关和非线性滤波相结合的算法,当声呐载体当前定位误差较大时,批相关用以迅速捕获当前地理区域,为非线性滤波连续输出地理位置信息跟踪导航缩小搜索空间,提高搜索效率,进而结合粒子滤波非线性滤波方法跟踪修正位置、速度导航信息;
其中,所述水下三维地形匹配定位导航声呐,包括控制系统(31)、 数据采集处理系统(32)、多通道接收机(33)、多通道稀疏平行接收线阵(34)、多通道信号源(35)、多通道发射机(36)和多通道发射基阵(37);所述控制系统(31)的输出端双向连接数据采集处理系统(32)的输入端,所述数据采集处理系统(32)的输出端分别连接多通道接收机(33)和多通道信号源(35)的输入端,所述多通道信号源(35)的输出端连接多通道发射机(36)的输入端,所述多通道发射机(36)的输出端连接多通道发射基阵(37)的输入端,所述多通道接收机(33)的输出端连接多通道稀疏平行接收线阵(34)的输入端;
所述多通道稀疏平行接收线阵(34)包括第一均匀间距接收线阵(39)、第二均匀间距接收线阵(40)和第三均匀间距接收线阵(41);所述第一均匀间距接收线阵(39)、第二均匀间距接收线阵(40)和第三均匀间距接收线阵(41)由等间距的多个接收阵元(38)组成,其间距由水平观测扇面大小决定;两组相邻的平行阵列间距不等,两个间距的比值服从互质关系。
CN201911153149.2A 2019-11-22 2019-11-22 一种水下三维地形匹配定位导航声呐及导航方法 Expired - Fee Related CN110907936B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911153149.2A CN110907936B (zh) 2019-11-22 2019-11-22 一种水下三维地形匹配定位导航声呐及导航方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911153149.2A CN110907936B (zh) 2019-11-22 2019-11-22 一种水下三维地形匹配定位导航声呐及导航方法

Publications (2)

Publication Number Publication Date
CN110907936A CN110907936A (zh) 2020-03-24
CN110907936B true CN110907936B (zh) 2021-11-16

Family

ID=69818773

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911153149.2A Expired - Fee Related CN110907936B (zh) 2019-11-22 2019-11-22 一种水下三维地形匹配定位导航声呐及导航方法

Country Status (1)

Country Link
CN (1) CN110907936B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113218386B (zh) * 2021-07-08 2021-10-08 深之蓝海洋科技股份有限公司 机器人在液体下建筑内高精度导航的方法和装置
CN114578438B (zh) * 2022-05-05 2022-07-12 成都理工大学 一种自适应水域电磁探测系统
CN116794663B (zh) * 2023-08-22 2023-11-14 北京海卓同创科技有限公司 一种基于多通道发射阵换能器的声呐及其控制方法
CN118226448A (zh) * 2024-03-25 2024-06-21 海底鹰深海科技股份有限公司 声纳及其装配方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1110709C (zh) * 2001-09-13 2003-06-04 中国科学院声学研究所 用于测量海底微地貌的高分辨测深侧扫声纳系统和测量方法
CN100430694C (zh) * 2006-12-31 2008-11-05 哈尔滨工程大学 多波束宽覆盖海底地形地貌探测装置
CN101907711A (zh) * 2010-07-16 2010-12-08 浙江大学 一种自主水下机器人的同时检测与定位方法
CN103047983B (zh) * 2013-01-14 2016-03-09 哈尔滨工程大学 水下机器人的面地形匹配导航方法
CN107314768B (zh) * 2017-07-06 2020-06-09 上海海洋大学 水下地形匹配辅助惯性导航定位方法及其定位系统
EP3432029A1 (en) * 2017-07-21 2019-01-23 Ebeson AB Underwater mapping method and system
CN110208812A (zh) * 2019-05-21 2019-09-06 哈尔滨工程大学 半潜无人航行器海底三维地形探测装置及方法

Also Published As

Publication number Publication date
CN110907936A (zh) 2020-03-24

Similar Documents

Publication Publication Date Title
CN110907936B (zh) 一种水下三维地形匹配定位导航声呐及导航方法
AU2014233495B2 (en) Systems and methods for navigating autonomous underwater vehicles
CN101900558B (zh) 集成声纳微导航的自主式水下机器人组合导航方法
CN111273298B (zh) 基于波浪滑翔器组网技术的水下声学目标定位与跟踪方法
CN110703203A (zh) 基于多声学波浪滑翔机的水下脉冲声定位系统
CN109870694B (zh) 基于多无人艇平台的高精度长基线定位系统
CN109782289B (zh) 一种基于基线几何结构约束的水下航行器定位方法
CN107664758B (zh) 基于长基线或超短基线组网的深海导航定位系统及方法
CN1779485A (zh) 载人潜水器的组合导航定位方法
CN115856898B (zh) 一种适应于全海深的面阵三维成像声呐点位归算方法
CN105115492A (zh) 基于声学多普勒计程仪的水下地形匹配导航系统
Châtillon et al. SAMI: A low-frequency prototype for mapping and imaging of the seabed by means of synthetic aperture
US5812494A (en) Wide-angle, forward-looking bathymetric mapping
CN114200401A (zh) 一种基于网格划分的水下机器人自定位系统及自定位方法
CN111735455A (zh) 基于改进的高斯距离迭代算法对接回收组合导航方法
CN110907937B (zh) 一种基于“t”型阵的掩埋物合成孔径三维成像方法
CN113156444B (zh) 一种基于运动补偿的多波束声呐高精度成像方法
CN113108778B (zh) 一种具备多条带模式的深水多波束测深方法及系统
Wang et al. Kalman filter based acoustic positioning of deep seafloor datum point with two-step systematic error estimation
CN111381212A (zh) 基于子阵划分的虚拟超短基线定位方法
CN113176539B (zh) 一种水声信号噪声多级抑制与稳健定位系统及定位方法
CN102183766B (zh) 一种合成孔径地层层析方法
CN116358544A (zh) 基于声学特征匹配定位修正惯性导航误差的方法和系统
RU2555479C2 (ru) Способ высокоточного координирования подводного комплекса в условиях подледного плавания
CN113009417B (zh) 利用声场干涉特性的海底声学阵列阵形估计方法

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20211116