CN110309581B - 一种水下潜标位置综合校准测量点快速优化布局方法 - Google Patents
一种水下潜标位置综合校准测量点快速优化布局方法 Download PDFInfo
- Publication number
- CN110309581B CN110309581B CN201910567230.9A CN201910567230A CN110309581B CN 110309581 B CN110309581 B CN 110309581B CN 201910567230 A CN201910567230 A CN 201910567230A CN 110309581 B CN110309581 B CN 110309581B
- Authority
- CN
- China
- Prior art keywords
- calibration
- submerged buoy
- signal
- measuring points
- points
- 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
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
- G01S5/00—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations
- G01S5/18—Position-fixing by co-ordinating two or more direction or position line determinations; Position-fixing by co-ordinating two or more distance determinations using ultrasonic, sonic, or infrasonic waves
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Development Economics (AREA)
- Radar, Positioning & Navigation (AREA)
- Game Theory and Decision Science (AREA)
- Geometry (AREA)
- Remote Sensing (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开了一种水下潜标位置综合校准测量点快速优化布局方法,所述布局方法包括以下步骤:步骤一:建立潜标位置校准模型;步骤二:确定不同声信号传播距离的时延测量误差关系式;步骤三:构建潜标位置单点校准精度数学模型;步骤四:选取潜标预估位置区域,在此区域内选取N个点利用步骤一至步骤三所述方法构造区域校准精度目标函数,当目标函数取得最小值时,此时的(xi,yi)即为测量点位置的最优值;步骤五:采用人工蜂群算法求解步骤四中的目标函数F,获得测量点优化布局结果。本发明的方法具有更符合实际情况、优化速度更快、测量更加准确等优点。
Description
技术领域
本发明属于水下潜标布局领域,提出了一种水下潜标位置综合校准测量点快速优化布局方法。
背景技术
开发和利用海洋资源必须首先探知海洋。为实现声学信息的采集及海洋探测等目的,潜标系统应运而生并且被广泛应用于实践当中。潜标系统是一种获取海洋环境信息的有效手段,它可以在恶劣的海洋环境条件下长期、连续、隐蔽地在水下工作,对过往船只、水中运动目标和各种海况下的环境噪声信息进行收集。同时潜标系统独立自主工作能力较好,自动化程度高。近年来,潜标系统已经被广泛应用于国防军事、海洋科学研究、水下工程前期调查、海洋开发等各个领域。潜标作为整个潜标系统的支撑点,在工作前必须首先对潜标位置进行校准,其校准精度直接影响整个系统的工作性能。因此,潜标位置校准问题一直是各国学者研究的重点问题。
目前海底潜标位置校准方法主要包括投放位置标示法、垂线相交法、超短基线定位校准法以及绝对校准法,其中绝对校准法是一种基于到达时间TOA(time of arrival)的校准方法。它主要采用测量船搭载船载声呐对待测潜标进行测定,测量船配备有卫星定位系统,然后通过测量信号由测量船到待测潜标之间的时延信息,采用几何交汇的方法解算得到待测潜标的绝对位置。传统的绝对校准过程中,一般选取成矩形的四个测量点上的相关信息利用球交汇模型进行解算。通常情况下,待测潜标深度由潜标自带的测深系统测出,当待测潜标深度已知时,球交汇模型退化为圆交汇模型。测量过程中,潜标入水后,由于水下水流环境十分复杂,导致潜标的工作位置距离投放位置有着一定偏差。此时,测量点布局将直接影响到潜标位置的校准精度。传统的测量点布局手段一般采用固定时延测量误差和蒙特卡洛方法进行,然而时延测量误差与声信号传播距离直接相关,采用固定误差分析不符合实际情况,同时蒙特卡洛方法耗时过长,不适合实时分析校准精度。此外,由于潜标偏离预定投放位置一定距离,此时若仍然采用预定投放点单点为中心点进行测量点布局与实际情况并不相符。
发明内容
本发明的目的是提出一种水下潜标位置综合校准测量点快速优化布局方法,以解决现有的水下潜标布局方法潜标位置精度差的问题,该方法能够以潜标预估区域为中心区域快速且合理地优化潜标位置校准的测量点布局,从而有效地提高潜标位置的校准精度。
本发明通过以下技术方案实现:一种水下潜标位置综合校准测量点快速优化布局方法,所述布局方法包括以下步骤:
步骤一:建立潜标位置校准模型;
步骤二:确定不同声信号传播距离的时延测量误差关系式;
步骤三:构建潜标位置单点校准精度数学模型;
步骤四:选取潜标预估位置区域,在此区域内选取N个点利用步骤一至步骤三所述方法构造区域校准精度目标函数,即:
其中,K为N个点中的一个点,HDOPk为第K点的校准精度,当目标函数取得最小值时,此时的(xi,yi)即为测量点位置的最优值;
步骤五:采用人工蜂群算法求解步骤四中的目标函数F,获得测量点优化布局结果。
进一步的,在步骤一中,具体的,选取四个测量点组成矩形,其校准方程为:
其中,(x,y)为待测阵元位置坐标,i为测量点序号,(xi,yi)为测量点位置坐标,ti为声信号由测量船到待测潜标所需的时延,h为为测量船声头与待测潜标之间深度差,c为海中声速。
进一步的,在步骤二中包括以下步骤:
步骤二一:基于Bellhop声场模型确定声速梯度以及声传播损失TL与声信号传播距离的关系,将确立的TL代入信噪比计算公式确定信噪比SNR:
SNR=SL-NL-TL
其中,SL为声源级,NL为海洋环境噪声级;
步骤二二:根据克拉美罗下界,获得时延测量误差dt的均方误差下界表示式:
其中,T1为信号持续时间,fH为信号的频率上限,fL为信号的频率下限,SNR为信噪比,T1、fH、fL参数由潜标位置校准系统确定,信噪比由上述Bellhop声场模型确定。
进一步的,在步骤三中,具体的:根据校准模型写出各项参数误差与校准精度dx、dy之间的关系:
其中,E()表示取期望,
用HDOP表述校准精度,则有:
本发明的有益效果在于:本发明提出了一种水下潜标位置综合校准测量点快速优化布局方法,该方法相对于传统方法的优势主要在于:
1)通过Bellhop声场模型确定不同声信号传播距离的时延测量误差关系式,使得时延测量误差更加准确,结果更加可信。
2)通过构建潜标位置校准精度数学模型进行校准精度分析,使得运算速度较蒙特卡洛方法得到了极大地提升,能够对测量点布局进行快速优化。
3)选择潜标预估位置区域构造区域校准精度目标函数而不是以投放点进行单点分析,与事实吻合,所获结果更贴近实际情况。
附图说明
图1为本发明的一种水下潜标位置综合校准测量点快速优化布局方法的方法流程图;
图2为声速梯度及传播损失与距离关系,其中图2(a)为声速梯度图;图2(b)为传播损失与距离关系图;
图3为水下潜标位置综合校准测量点快速优化布局结果示意图。
具体实施方式
下面将结合本发明实施例中的附图对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
参照图1所示,本发明通过以下技术方案实现:一种水下潜标位置综合校准测量点快速优化布局方法,所述布局方法包括以下步骤:
步骤一:建立潜标位置校准模型;
步骤二:确定不同声信号传播距离的时延测量误差关系式;
步骤三:构建潜标位置单点校准精度数学模型;
步骤四:选取潜标预估位置区域,在此区域内选取N个点利用步骤一至步骤三所述方法构造区域校准精度目标函数,即:
其中,K为N个点中的一个点,HDOPk为第K点的校准精度,当目标函数取得最小值时,此时的(xi,yi)即为测量点位置的最优值;
步骤五:采用人工蜂群算法求解步骤四中的目标函数F,获得测量点优化布局结果。
参照图1所示,在本部分优选实施例中,在步骤一中,具体的,选取四个测量点组成矩形,其校准方程为:
其中,(x,y)为待测阵元位置坐标,i为测量点序号,(xi,yi)为测量点位置坐标,ti为声信号由测量船到待测潜标所需的时延,h为为测量船声头与待测潜标之间深度差,c为海中声速。
参照图1所示,在本部分优选实施例中,在步骤二中包括以下步骤:
步骤二一:基于Bellhop声场模型确定声速梯度以及声传播损失TL与声信号传播距离的关系,将确立的TL代入信噪比计算公式确定信噪比SNR:
SNR=SL-NL-TL
其中,SL为声源级,NL为海洋环境噪声级;
步骤二二:根据克拉美罗下界,获得时延测量误差dt的均方误差下界表示式:
其中,T1为信号持续时间,fH为信号的频率上限,fL为信号的频率下限,SNR为信噪比,T1、fH、fL等参数由潜标位置校准系统确定,信噪比由上述Bellhop声场模型确定。
参照图2所示,在本部分优选实施例中,在步骤三中,具体的:由于步骤一所述校准模型误差来源于各项参数测量误差,即:声速测量误差dc,时延测量误差dti,深度差测量误差dh,测量点位置GPS定位误差dxi、dyi。
根据校准模型写出各项参数误差与校准精度dx、dy之间的关系:
其中,E()表示取期望,
用HDOP表述校准精度,则有:
下面提供一个具体的实施算例:
采用仿真数据对本发明所设计的水下潜标位置综合校准测量点快速优化布局方法进行验证,并对过程结果进行说明。
首先给出各参数如下:测量船声头入水6m,测量船运动速度为8m/s,声信号频率上限fL为10kHz,频率下限fH为11kHz。声速测量误差dc为0.5m/s,时延测量误差dti由Bellhop声学模型确定。发射声源级SL为195dB,噪声级NL为三级海况时噪声谱级63dB。深度差测量误差dh为0.2m,测量点位置GPS定位误差为2m。潜标预估位置区域为以坐标原点(0,0)为圆心,半径为25m的圆形区域。潜标深度为200m,选取区域内5000个点进行解算。首先给出Bellhop声学模型确定的声速梯度及声传播损失与距离关系图,如图2所示。最后即可获得测量点优化布局结果如图3所示。解算所用时间约为285s,区域校准精度约为2.49m。可以得到四个最佳测量点位置分别为(449.1,449.3)、(449.1,-449.3)、(-449.1,-449.3)、(-449.1,449.3)。
仿真数据处理结果说明本发明所设计的方法可以对测量点布局进行快速优化,从而有效提高潜标位置校准精度。
Claims (4)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910567230.9A CN110309581B (zh) | 2019-06-27 | 2019-06-27 | 一种水下潜标位置综合校准测量点快速优化布局方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910567230.9A CN110309581B (zh) | 2019-06-27 | 2019-06-27 | 一种水下潜标位置综合校准测量点快速优化布局方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110309581A CN110309581A (zh) | 2019-10-08 |
CN110309581B true CN110309581B (zh) | 2022-11-01 |
Family
ID=68076311
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910567230.9A Active CN110309581B (zh) | 2019-06-27 | 2019-06-27 | 一种水下潜标位置综合校准测量点快速优化布局方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110309581B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111046537B (zh) * | 2019-11-26 | 2022-12-23 | 中国船舶重工集团有限公司第七一0研究所 | 一种潜标系统布放时水下分离系留索强度仿真方法 |
CN115508779B (zh) * | 2022-10-31 | 2023-09-29 | 浙江大学 | 一种高速落点目标的定位方法及装置 |
CN116106875B (zh) * | 2023-04-07 | 2023-07-14 | 北京神州普惠科技股份有限公司 | 岸基阵坐标联合校准方法、系统、电子设备及存储介质 |
Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103529441A (zh) * | 2012-07-02 | 2014-01-22 | 中国科学院声学研究所 | 一种被动合成孔径目标信号检测和分辨方法及系统 |
CN105703866A (zh) * | 2015-12-29 | 2016-06-22 | 南京世海声学科技有限公司 | 一种水下线列阵的使用网络恢复时钟的同步采集方法 |
CN106546954A (zh) * | 2016-11-09 | 2017-03-29 | 哈尔滨工程大学 | 一种深海超短基线定位精度检验方法 |
CN106959429A (zh) * | 2017-04-25 | 2017-07-18 | 中国人民解放军海军工程大学 | 基于多正交信号的水下区域导航标校方法及系统 |
EP3193187A1 (en) * | 2016-01-13 | 2017-07-19 | ST SPORTSERVICE GmbH | Method for calibrating a local positioning system based on time-difference-of-arrival measurements |
CN107064876A (zh) * | 2017-06-20 | 2017-08-18 | 哈尔滨工业大学 | 基于无线传感网络和超声波的室内定位系统及定位方法 |
WO2017178778A1 (en) * | 2016-04-11 | 2017-10-19 | Statsports Group Limited | System for enhanced indoor and outdoor position measurement |
CN107272004A (zh) * | 2017-07-24 | 2017-10-20 | 哈尔滨工程大学 | 一种基于水下信标位置修正的单信标测距定位方法 |
CN107479031A (zh) * | 2017-07-24 | 2017-12-15 | 哈尔滨工程大学 | 一种隐蔽的水下被动定位方法 |
CN107678032A (zh) * | 2017-07-21 | 2018-02-09 | 哈尔滨工程大学 | 一种基于虚拟收发信标的单信标测距定位方法 |
CN107990891A (zh) * | 2016-10-26 | 2018-05-04 | 中国科学院沈阳自动化研究所 | 基于长基线和信标在线标定的水下机器人组合导航方法 |
CN108650629A (zh) * | 2018-04-26 | 2018-10-12 | 兰州理工大学温州泵阀工程研究院 | 一种基于无线通信基站的室内三维定位算法 |
CN109241576A (zh) * | 2018-08-14 | 2019-01-18 | 西安电子科技大学 | 基于蚁狮算法的稀疏天线方向图综合方法 |
CN109357677A (zh) * | 2018-10-22 | 2019-02-19 | 哈尔滨工程大学 | 一种应用于水下单信标导航的航路规划方法 |
CN109490927A (zh) * | 2018-12-26 | 2019-03-19 | 天津水运工程勘察设计院 | 一种水下整平架定位系统及其定位方法 |
CN109738902A (zh) * | 2019-03-06 | 2019-05-10 | 哈尔滨工程大学 | 一种基于同步信标模式的水下高速目标高精度自主声学导航方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB201200831D0 (en) * | 2012-01-18 | 2012-02-29 | Sensewhere Ltd | Improved positioning system |
US9568584B2 (en) * | 2012-03-15 | 2017-02-14 | Soongsil University Research Consortium Techno-Park | Wireless positioning server using clock offset calibration positioning method using the same |
-
2019
- 2019-06-27 CN CN201910567230.9A patent/CN110309581B/zh active Active
Patent Citations (16)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103529441A (zh) * | 2012-07-02 | 2014-01-22 | 中国科学院声学研究所 | 一种被动合成孔径目标信号检测和分辨方法及系统 |
CN105703866A (zh) * | 2015-12-29 | 2016-06-22 | 南京世海声学科技有限公司 | 一种水下线列阵的使用网络恢复时钟的同步采集方法 |
EP3193187A1 (en) * | 2016-01-13 | 2017-07-19 | ST SPORTSERVICE GmbH | Method for calibrating a local positioning system based on time-difference-of-arrival measurements |
WO2017178778A1 (en) * | 2016-04-11 | 2017-10-19 | Statsports Group Limited | System for enhanced indoor and outdoor position measurement |
CN107990891A (zh) * | 2016-10-26 | 2018-05-04 | 中国科学院沈阳自动化研究所 | 基于长基线和信标在线标定的水下机器人组合导航方法 |
CN106546954A (zh) * | 2016-11-09 | 2017-03-29 | 哈尔滨工程大学 | 一种深海超短基线定位精度检验方法 |
CN106959429A (zh) * | 2017-04-25 | 2017-07-18 | 中国人民解放军海军工程大学 | 基于多正交信号的水下区域导航标校方法及系统 |
CN107064876A (zh) * | 2017-06-20 | 2017-08-18 | 哈尔滨工业大学 | 基于无线传感网络和超声波的室内定位系统及定位方法 |
CN107678032A (zh) * | 2017-07-21 | 2018-02-09 | 哈尔滨工程大学 | 一种基于虚拟收发信标的单信标测距定位方法 |
CN107479031A (zh) * | 2017-07-24 | 2017-12-15 | 哈尔滨工程大学 | 一种隐蔽的水下被动定位方法 |
CN107272004A (zh) * | 2017-07-24 | 2017-10-20 | 哈尔滨工程大学 | 一种基于水下信标位置修正的单信标测距定位方法 |
CN108650629A (zh) * | 2018-04-26 | 2018-10-12 | 兰州理工大学温州泵阀工程研究院 | 一种基于无线通信基站的室内三维定位算法 |
CN109241576A (zh) * | 2018-08-14 | 2019-01-18 | 西安电子科技大学 | 基于蚁狮算法的稀疏天线方向图综合方法 |
CN109357677A (zh) * | 2018-10-22 | 2019-02-19 | 哈尔滨工程大学 | 一种应用于水下单信标导航的航路规划方法 |
CN109490927A (zh) * | 2018-12-26 | 2019-03-19 | 天津水运工程勘察设计院 | 一种水下整平架定位系统及其定位方法 |
CN109738902A (zh) * | 2019-03-06 | 2019-05-10 | 哈尔滨工程大学 | 一种基于同步信标模式的水下高速目标高精度自主声学导航方法 |
Non-Patent Citations (7)
Title |
---|
Juhyun Pyo.On-Site Calibration-Based Estimation Method of Forward Seabed Elevation Using Forward Scan Sonar.《IEEE Sensors Journal》.2019,8832-8844. * |
Kenneth G. Foote.Organizing error sources in a standard on active-sonar calibration by standard target.《OCEANS 2018 MTS/IEEE Charleston》.2019,1-5. * |
周天立.深海声学潜标声基元位置校准系统显控软件设计与实现.《中国优秀硕士学位论文全文数据库 信息科技辑》.2018,I138-437. * |
张涛.基于惯导及声学浮标辅助的水下航行器导航定位系统.《中国惯性技术学报》.2017,741-745. * |
曹俊.水下载体的声学定位系统测距误差分析.《南京大学学报(自然科学)》.2016,816-824. * |
梁国龙.基于稀疏信号重构的阵元位置误差校正方法.《声学学报》.2017,677-684. * |
郭锦标.多正交信号水下导航系统位置校准的CRLB研究.《舰船电子工程》.2019,143-147. * |
Also Published As
Publication number | Publication date |
---|---|
CN110309581A (zh) | 2019-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11733381B2 (en) | Sound velocity profile inversion method based on inverted multi-beam echo sounder | |
CN106595834B (zh) | 一种获得深海大深度声场水平纵向相关性的方法 | |
CN110309581B (zh) | 一种水下潜标位置综合校准测量点快速优化布局方法 | |
CN107664758B (zh) | 基于长基线或超短基线组网的深海导航定位系统及方法 | |
CN109738902B (zh) | 一种基于同步信标模式的水下高速目标高精度自主声学导航方法 | |
CN111896962B (zh) | 一种海底应答器定位方法、系统、存储介质及应用 | |
CN107167224B (zh) | 一种船舶辐射噪声的测量方法 | |
CN110260044B (zh) | 一种海底管道定位方法 | |
CN113671443B (zh) | 基于掠射角声线修正的水声传感器网络深海目标定位方法 | |
CN110133627B (zh) | 水声定位导航系统阵元位置校准测量点间距优化方法 | |
CN112540348A (zh) | 一种基于空间尺度上的声线修正算法在长基线水声定位系统上的应用 | |
CN111220146B (zh) | 一种基于高斯过程回归学习的水下地形匹配定位方法 | |
Xin et al. | A TOA/AOA underwater acoustic positioning system based on the equivalent sound speed | |
CN110132281B (zh) | 一种基于询问应答模式的水下高速目标高精度自主声学导航方法 | |
CN109991670B (zh) | 一种水下目标定位方法 | |
CN108387872B (zh) | 基于最大偏移量法的超短基线定位优化方法 | |
CN108761470B (zh) | 一种基于拖缆形态方程解析的目标定位方法 | |
CN117574627A (zh) | 一种多波束测深系统测线长度计算模型的构建方法 | |
CN116840786A (zh) | 一种面向长基线水下定位系统的声线修正方法 | |
CN110531319B (zh) | 一种基于测量点优化布局的水声应答器位置高精度校准方法 | |
CN116973899B (zh) | 一种基于凸聚类的冰下多目标定位方法和装置 | |
CN117146830B (zh) | 一种自适应多信标航位推算和长基线的紧组合导航方法 | |
CN109490868B (zh) | 一种基于分布式垂直线列阵的海上目标运动分析方法 | |
Huang et al. | An improvement of long baseline system using particle swarm optimization to optimize effective sound speed | |
CN115629389A (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 |