CN110717138B - 一种拖曳式瞬变电磁数据预处理方法 - Google Patents
一种拖曳式瞬变电磁数据预处理方法 Download PDFInfo
- Publication number
- CN110717138B CN110717138B CN201910861795.8A CN201910861795A CN110717138B CN 110717138 B CN110717138 B CN 110717138B CN 201910861795 A CN201910861795 A CN 201910861795A CN 110717138 B CN110717138 B CN 110717138B
- Authority
- CN
- China
- Prior art keywords
- data
- measuring point
- correlation
- point data
- transient electromagnetic
- 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
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- 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/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种拖曳式瞬变电磁数据预处理方法,包括:获取全时拖曳瞬变电磁数据,从第一个测点开始,依次计算该测点数据与其它测点的相关度值并记录,直至滑动完成所有测点数据的相关度计算;从第一个测点开始,将与该测点相关度值大于设置阈值的其他测点数据进行相关加权叠加,权值为对应的相关度值,将叠加后的数据记为新的测点数据,直至完成所有测点数据的相关加权叠加计算;所有相关叠加后的新测点数据为双极性梯形波发射的数据,将其正极性周期发射关断后对应的数据减去负极性周期发射关断后对应的数据后除以二,得到预处理后的数据。本发明能够在不减少测点数量的同时有效的提高测点数据信噪比。
Description
技术领域
本发明涉及地球物理信号处理与分析技术研究领域,特别是涉及一种拖曳式瞬变电磁数据预处理方法。
背景技术
拖曳式瞬变电磁探测是指由人工或者全地形车牵引的地面瞬变电磁系统,其发射系统能够在拖曳移动过程中连续进行工作,接收系统同步采集全过程全波形信号数据;相比于传统地面瞬变电磁定点布局探测,该拖曳式探测装置能够更高效快速的对地下地质结构进行探测,因此越发受到科研工作者的重视。然而,拖曳式瞬变电磁装置由于连续移动探测的特殊性,采集的数据易受到工频噪声、尖峰噪声、随机噪声、运动噪声的影响,由于在连续运动,其数据处理不能够采用地面瞬变电磁法的定点多次叠加的方法来提高数据信噪比;又因为其移动速度相对较慢,测线距离短,也无法采用航空瞬变电磁法相邻测点数据直接叠加的方法,会丢失横向分辨率。因此,研究适用于拖曳式瞬变电磁装置的数据预处理方法,在保证数据解释横向分辨率的同时提高数据信噪比,对拓展拖曳式瞬变电磁装置的应用范围,实现地下空间的高效探测与高精度解释具有重要意义。
发明内容
本发明所要解决的技术问题在于提供一种拖曳式瞬变电磁数据预处理方法,解决现有的方法不能满足拖曳式瞬变电磁数据的处理,没法保证数据解释横向分辨率的同时提高数据信噪比的问题。
本发明是这样实现的,
一种拖曳式瞬变电磁数据预处理方法,该方法包括以下步骤:
S1、获取瞬变电磁仪器连续拖曳移动发射时采集的全时拖曳瞬变电磁数据,根据发射频率与接收机采样率提取每个发射周期下接收的瞬变电磁数据,每个发射周期采集的数据记为一个测点数据;
S2、设置相关度阈值,从第一个测点开始,依次计算该测点数据与其它测点的相关度值并记录,随后滑至下一个测点数据进行相关度计算,直至滑动完成所有测点数据的相关度计算;
S3、从第一个测点开始,将与该测点相关度值大于设置阈值的其他测点数据进行相关加权叠加,权值为对应的相关度值,将叠加后的数据记为新的测点数据,随后滑至下一个测点数据进行相关加权叠加计算,直至完成所有测点数据的相关加权叠加计算;
S4、所有相关叠加后的新测点数据为双极性梯形波发射的数据,将其正极性周期发射关断后对应的数据减去负极性周期发射关断后对应的数据后除以二,得到的数据记为拖曳式瞬变电磁预处理后的数据;
S5、所有测点完成上述计算后,完成拖曳式瞬变电磁数据预处理。
进一步地,所述步骤S1中的瞬变电磁仪器是由人工或者全地形车牵引的拖曳式瞬变电磁系统,探测线圈放置于移动平台上,采用中心回线配置,发射系统控制发射线圈在拖曳过程中连续发射,接收系统通过接收线圈同步采集全波形信号数据。
进一步地,所述步骤S1中的测点数据,其数据长度Lm由发射频率fv与接收机采样率fs定义,将拖曳式数据中每Lm个数据记为一个测点数据,表示为:
进一步地,所述步骤S2相关度计算公式满足:
进一步地,所述步骤S2中相关度阈值是经验阈值系数,只有当计算的其它相邻测点相关度高于这个阈值时,才进行相关加权叠加,通过降低阈值增加相关叠加数目,提高视电阻率成像纵向分辨率,通过提高阈值增加成像横向异常分辨率,阈值选取的范围为0.5到0.9之间。
进一步地,所述步骤S3中相关加权叠加具体包括:通过计算目标测点与全部测点的相关度值,选取相关度值高于设置阈值的测点,将其作为权系数分别乘以对应的测点数据,叠加到目标测点,最后除以叠加的权系数之和进行数据归一化,得到目标测点叠加值Um满足:
其中U为测点数据,Uk为对应的第k个测点数据序列,N为总测点数目,rk为第k测点数据与第m测点数据的相关度值,rmin为相关度阈值,Um为加权叠加后的第m测点的数据。
进一步地,所述步骤S4中双极性梯形波发射为:拖曳式瞬变电磁发射系统的连续发射波形,双极性梯形波发射周期为工频周期的偶数倍以抑制工频噪声的干扰,正负极性周期发射关断后对应的数据相减满足:
本发明与现有技术相比,有益效果在于:本发明提出拖曳式瞬变电磁探测数据预处理方法,采用滑动相关加权叠加算法,首先将拖曳式瞬变电磁整体数据按发射周期与采样率的乘积转换为多测点数据,随后采用相关加权叠加算法,依次滑动计算各个测点叠加后的值来保留数据的横向异性,最后通过正反双极性周期数据相减来进一步抑制噪声,完成数据预处理。该方法能够在不减少测点数量的同时有效的提高测点数据信噪比,保证了后期数据解释的横向分辨率与准确度,具有较大的实际应用价值。
附图说明
图1为拖曳式瞬变电磁探测装置示意图;
图2为拖曳式瞬变电磁实测测线数据图;
图3为拖曳式瞬变电磁数据预处理方法流程图;
图4为拖曳式瞬变电磁实测数据单个周期测点图;
图5为拖曳式瞬变电磁实测数据测点相关度图;
图6为滑动相关加权叠加算法示意图;
图7为拖曳式瞬变电磁单个周期仅极性相减处理后原始实测数据图;
图8为拖曳式瞬变电磁单个周期直接叠加处理后原始实测数据图;
图9为拖曳式瞬变电磁单个周期预处理后原始实测数据图;
图10为拖曳式瞬变电磁预处理后视电阻率成像图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白,以下结合实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。
基于滑动相关加权叠加算法的拖曳式瞬变电磁探测数据预处理方法,包括以下步骤:
S1、获取瞬变电磁仪器连续拖曳移动发射时采集的全时拖曳瞬变电磁数据,根据发射频率与接收机采样率提取每个发射周期下接收的瞬变电磁数据,每个发射周期采集的数据记为一个测点数据;
S2、设置相关度阈值,从第一个测点开始,依次计算该测点数据与其它测点的相关度值并记录,随后滑至下一个测点数据进行相关度计算,直至滑动完成所有测点数据的相关度计算;
S3、从第一个测点开始,将与该测点相关度值大于设置阈值的其他测点数据进行相关加权叠加,权值为对应的相关度值,将叠加后的数据记为新的测点数据,随后滑至下一个测点数据进行相关加权叠加计算,直至完成所有测点数据的相关加权叠加计算;
S4、所有相关叠加后的新测点数据为双极性梯形波发射的数据,将其正极性周期发射关断后对应的数据减去负极性周期发射关断后对应的数据后除以二,得到的数据记为拖曳式瞬变电磁预处理后的数据;
S5、所有测点完成上述计算后,即完成拖曳式瞬变电磁数据预处理;
所述步骤S1中的瞬变电磁仪器是由人工或者全地形车牵引的拖曳式瞬变电磁系统,探测线圈放置于移动平台上,采用中心回线配置,发射系统控制发射线圈在拖曳过程中连续发射,接收系统通过接收线圈同步采集全波形信号数据;
所述步骤S1中的测点数据,其数据长度Lm由发射频率fv与接收机采样率fs定义,将拖曳式数据中每Lm个数据记为一个测点数据;
所述步骤S2的相关度计算公式满足:
所述步骤S2中的相关度阈值是经验阈值系数,只有当计算的其它相邻测点相关度高于这个阈值时,才进行相关加权叠加,合理的选取阈值有利于在保证信噪比的同时提高后期反演成像解释横向异常分辨率,阈值选取的范围为0.5到0.9之间;
所述步骤S3中的相关加权叠加是通过计算目标测点与全部测点的相关度值,选取相关度值高于设置阈值的测点,将其作为权系数分别乘以对应的测点数据,叠加到目标测点,最后除以叠加的权系数之和进行数据归一化,得到目标测点叠加值Um满足:
其中U为测点数据,Uk为对应的第k个测点数据序列,N为总测点数目,rk为第k测点数据与第m测点数据的相关度值,rmin为相关度阈值,Um为加权叠加后的第m测点的数据;
所述步骤S4中的双极性梯形波发射是拖曳式瞬变电磁发射系统的连续发射波形,双极性梯形波发射周期为工频周期的偶数倍以抑制工频噪声的干扰,正负极性周期发射关断后对应的数据相减满足:
基于滑动相关加权叠加算法的拖曳式瞬变电磁探测数据预处理方法的详细说明:
(1)、本实施例中拖曳式瞬变电磁探测装置采用入如图1所示的结构,该装置采用中心回线方式探测,包括拖曳车5,以及设置在拖曳平台上的发射系统4和接收系统3,发射系统4在拖曳过程中连续发射,接收系统3同步采集全波形信号数据,发射线圈2为4匝,尺寸为2m×2m,发射频率为12.5Hz,接收线圈1为1匝,半径为0.25m。发射波形为双极性梯形波,接收系统采样率为156.25K,在一段道路上进行拖曳式测量后采集到的整体测线数据X(n),拖曳式瞬变电磁测线数据见图2,数据预处理流程如图3所示;
(2)、根据发射周期和接收机采样率得出一个发射周期对应的数据量为12500个,将原始实测数据每12500个记为一个测点数据U,图4为单个发射周期对应的测点数据图;
(3)、按照下面相关度计算公式,从第一个测点开始,依次计算该测点数据与其它测点的相关度值r并记录,随后滑至下一个测点数据进行相关度计算,直至滑动完成所有测点数据的相关度计算,图5为第二十个测点与其相邻的测点(第一至第一百测点)的相关度图;
(4)、设置相关度阈值rmin,只有当计算的其它相邻测点相关度高于这个阈值时,才进行后面的相关加权叠加,合理的选取阈值有利于在保证信噪比的同时提高后期反演成像解释横向异常分辨率,阈值选取的范围一般为0.6到0.9之间,本文将阈值设置为0.7;
(5)、进行相关加权叠加,将与目标m测点相关度计算值高于阈值的测点进行加权叠加得到Um,权系数rk为计算的相关度值,相关加权叠加算法示意图如图6所示,计算公式如下;
(7)、所有测点完成上述计算后,即完成了拖曳式瞬变电磁数据的预处理;
为了验证该数据预处理方法的有效性,随机选取一个测点数据(第二二十测点)进行分析,图7为该测点仅通过双极性波数据相减处理得到的数据图形,图8为采用直接叠加法处理后的数据结果,即每隔100个测点数据直接叠加合成为一个测点数据,图9为采用相关加权叠加算法预处理后的数据结果;
基于滑动相关加权叠加的拖曳式瞬变电磁探测数据预处理方法,能够在不减少测点数量的同时有效的提高测点数据信噪比,保证后期数据解释的横向分辨率与准确度。对比图7、图8与图9可知,拖曳式瞬变电磁数据在不经过叠加处理时晚期数据信噪比较低,然而采用直接叠加的方法,虽然叠加100次后晚期信噪比明显提高,但是原始数据的幅度减小明显,数据严重失真,采用滑动相关加权叠加算法预处理后的数据,在提高信噪比的同时,保持了信号的真实性,为后期数据解释奠定了重要基础;图10为采用滑动相关加权叠加算法预处理后的拖曳式瞬变电磁数据视电阻率解释结果。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种拖曳式瞬变电磁数据预处理方法,其特征在于,该方法包括以下步骤:
S1、获取瞬变电磁仪器连续拖曳移动发射时采集的全时拖曳瞬变电磁数据,根据发射频率与接收机采样率提取每个发射周期下接收的瞬变电磁数据,每个发射周期采集的数据记为一个测点数据;
S2、设置相关度阈值,从第一个测点开始,依次计算该测点数据与相邻第一至第一百的其它测点的相关度值并记录,随后滑至下一个测点数据进行相关度计算,直至滑动完成所有测点数据的相关度计算;所述步骤S2中相关度阈值是经验阈值系数,只有当计算的其它相邻测点相关度高于这个阈值时,才进行相关加权叠加,通过降低阈值增加相关叠加数目,提高视电阻率成像纵向分辨率,通过提高阈值增加成像横向异常分辨率,阈值选取的范围为0.5到0.9之间;
S3、从第一个测点开始,将与该测点相关度值大于设置阈值的其他测点数据进行相关加权叠加,权值为对应的相关度值,将叠加后的数据记为新的测点数据,随后滑至下一个测点数据进行相关加权叠加计算,直至完成所有测点数据的相关加权叠加计算;
S4、所有相关叠加后的新测点数据为双极性梯形波发射的数据,将其正极性周期发射关断后对应的数据减去负极性周期发射关断后对应的数据后除以二,得到的数据记为拖曳式瞬变电磁预处理后的数据;
S5、所有测点完成上述计算后,完成拖曳式瞬变电磁数据预处理。
2.按照权利要求1所述的方法,其特征在于,所述步骤S1中的瞬变电磁仪器是由人工或者全地形车牵引的拖曳式瞬变电磁系统,探测线圈放置于移动平台上,采用中心回线配置,发射系统控制发射线圈在拖曳过程中连续发射,接收系统通过接收线圈同步采集全波形信号数据。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910861795.8A CN110717138B (zh) | 2019-09-12 | 2019-09-12 | 一种拖曳式瞬变电磁数据预处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910861795.8A CN110717138B (zh) | 2019-09-12 | 2019-09-12 | 一种拖曳式瞬变电磁数据预处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110717138A CN110717138A (zh) | 2020-01-21 |
CN110717138B true CN110717138B (zh) | 2023-06-30 |
Family
ID=69210379
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910861795.8A Active CN110717138B (zh) | 2019-09-12 | 2019-09-12 | 一种拖曳式瞬变电磁数据预处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110717138B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111796334B (zh) * | 2020-07-20 | 2021-08-17 | 吉林大学 | 一种高分辨率拖曳式瞬变电磁数据预处理方法 |
CN111983698B (zh) * | 2020-08-20 | 2023-05-12 | 吉林大学 | 基于成分相位反相的拖曳瞬变电磁数据工频噪声抑制方法 |
CN115542408B (zh) * | 2022-12-05 | 2023-03-28 | 成都理工大学 | 一种海洋瞬变电磁数据预处理与快速定区成像方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2388611A1 (en) * | 2010-05-19 | 2011-11-23 | Jeol Resonance Inc. | Production of spectra with improved signal-to-noise ratio |
CN102495431A (zh) * | 2011-11-16 | 2012-06-13 | 中煤科工集团西安研究院 | 瞬变电磁数据静校正的数据处理方法 |
CN105158808A (zh) * | 2015-07-01 | 2015-12-16 | 长安大学 | 一种浅海瞬变电磁海空探测及其解释方法 |
CN106814402A (zh) * | 2016-12-22 | 2017-06-09 | 中石化石油工程技术服务有限公司 | 瞬变电磁信号叠前去噪方法 |
CN109507739A (zh) * | 2018-12-12 | 2019-03-22 | 吉林大学 | 城市地下空间快速高精度拖曳式阵列电磁探测装置及探测方法 |
CN109597132A (zh) * | 2018-12-26 | 2019-04-09 | 中国科学院声学研究所 | 一种浅海磁性源瞬变电磁探测装置及探测方法 |
CN109946744A (zh) * | 2019-04-01 | 2019-06-28 | 吉林大学 | 一种基于可控源补偿的瞬变电磁勘探系统及方法 |
-
2019
- 2019-09-12 CN CN201910861795.8A patent/CN110717138B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2388611A1 (en) * | 2010-05-19 | 2011-11-23 | Jeol Resonance Inc. | Production of spectra with improved signal-to-noise ratio |
CN102495431A (zh) * | 2011-11-16 | 2012-06-13 | 中煤科工集团西安研究院 | 瞬变电磁数据静校正的数据处理方法 |
CN105158808A (zh) * | 2015-07-01 | 2015-12-16 | 长安大学 | 一种浅海瞬变电磁海空探测及其解释方法 |
CN106814402A (zh) * | 2016-12-22 | 2017-06-09 | 中石化石油工程技术服务有限公司 | 瞬变电磁信号叠前去噪方法 |
CN109507739A (zh) * | 2018-12-12 | 2019-03-22 | 吉林大学 | 城市地下空间快速高精度拖曳式阵列电磁探测装置及探测方法 |
CN109597132A (zh) * | 2018-12-26 | 2019-04-09 | 中国科学院声学研究所 | 一种浅海磁性源瞬变电磁探测装置及探测方法 |
CN109946744A (zh) * | 2019-04-01 | 2019-06-28 | 吉林大学 | 一种基于可控源补偿的瞬变电磁勘探系统及方法 |
Non-Patent Citations (6)
Title |
---|
Tingting Lin ; Yang Zhang ; Ling Wan ; Yongxing Qu ; Jun Lin.A Para-Whole Space Model for Underground Magnetic Resonance Sounding Studies.《IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing》.2016,第9卷(第9期),264 - 271. * |
卢超,吴伟,杨琳瑜,彭应秋.超声相关加权合成断面成像检测方法研究.南昌航空工业学院学报(自然科学版).2004,18(03),25-28. * |
戚志鹏 ; 李貅 ; 郭建磊 ; 孙怀凤 ; 姚伟华 ; 吴琼 ; .基于微分电导的航空瞬变电磁合成孔径快速成像方法研究.地球物理学进展.2015,30(04),1903-1911. * |
林君 ; 张洋 ; .地面磁共振探水技术的研究现状与展望.仪器仪表学报.2016,37(12),2657-2670. * |
桑圣华 ; .基于矿井瞬变电磁深度格式网格化新方法的研究.能源技术与管理.2013,38(05),142-143. * |
赵虎 ; 李瑞 ; 王军 ; .浅水域拖曳式瞬变电磁系统在公路勘察中的应用研究.路基工程.2018,(S1),111-114、119. * |
Also Published As
Publication number | Publication date |
---|---|
CN110717138A (zh) | 2020-01-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110717138B (zh) | 一种拖曳式瞬变电磁数据预处理方法 | |
CN104020495B (zh) | 一种基于探地雷达的地下管线参数自识别方法 | |
CN105158808B (zh) | 一种浅海瞬变电磁海空探测及其解释方法 | |
CN105549097B (zh) | 一种瞬变电磁信号工频及其谐波干扰消除方法及装置 | |
CN102520444B (zh) | 一种叠后地震波中绕射波信息提取方法 | |
CN102819043B (zh) | 阵列信号随机噪声自适应模型去噪方法 | |
CN105785455B (zh) | 一种基于b样条插值的二维地面核磁共振反演方法 | |
CN105785439B (zh) | 小尺度非均匀地质体空间分布位置的预测方法和装置 | |
CN104678435A (zh) | 一种提取Rayleigh面波频散曲线的方法 | |
CN101598809A (zh) | 一种自适应消除线性规则噪声以及多次波干扰的方法 | |
CN101650443B (zh) | 视电阻率的反向传播网络计算方法 | |
CN102692647B (zh) | 一种高时间分辨率的地层含油气性预测方法 | |
CN103364832A (zh) | 一种基于自适应最优核时频分布的地震衰减定性估计方法 | |
CN103777247A (zh) | 一种瞬变电磁雷达探测系统及探测方法 | |
CN101923176B (zh) | 一种利用地震数据瞬时频率属性进行油气检测的方法 | |
CN103163554A (zh) | 利用零偏vsp资料估计速度和q值的自适应波形的反演方法 | |
CN105629317A (zh) | 一种基于站间传递函数的大地电磁噪声压制方法 | |
CN102176064A (zh) | 煤矿井下瞬变电磁仪去除电磁干扰信号的方法 | |
CN105044793A (zh) | 一种多道瞬变电磁探测数据的反演方法和装置 | |
CN105445801A (zh) | 一种消除二维地震资料随机噪音的处理方法 | |
CN110554428A (zh) | 一种基于变分模态分解的地震波低频能量变化率提取方法 | |
CN107807393A (zh) | 基于地震干涉法的单台站集初至波增强方法 | |
CN103217709B (zh) | 一种提高地震数据信噪比和分辨率的面波衰减方法 | |
Sun et al. | An efficient preprocessing method to suppress power harmonic noise for urban towed transient electromagnetic measurement data | |
CN105093318A (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 |