CN115685355A - 用于随钻方位电磁波电阻率仪器的幅相提取方法及装置 - Google Patents
用于随钻方位电磁波电阻率仪器的幅相提取方法及装置 Download PDFInfo
- Publication number
- CN115685355A CN115685355A CN202211355755.4A CN202211355755A CN115685355A CN 115685355 A CN115685355 A CN 115685355A CN 202211355755 A CN202211355755 A CN 202211355755A CN 115685355 A CN115685355 A CN 115685355A
- Authority
- CN
- China
- Prior art keywords
- signal
- phase
- amplitude
- sampling
- analog
- 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.)
- Granted
Links
- 238000005553 drilling Methods 0.000 title claims abstract description 28
- 238000000605 extraction Methods 0.000 title claims abstract description 16
- 238000005070 sampling Methods 0.000 claims abstract description 49
- 238000012545 processing Methods 0.000 claims abstract description 22
- 238000000034 method Methods 0.000 claims abstract description 17
- 238000006243 chemical reaction Methods 0.000 claims abstract description 15
- 230000003750 conditioning effect Effects 0.000 claims abstract description 15
- 238000001914 filtration Methods 0.000 claims abstract description 9
- 230000003321 amplification Effects 0.000 claims abstract description 6
- 239000000284 extract Substances 0.000 claims abstract description 6
- 238000003199 nucleic acid amplification method Methods 0.000 claims abstract description 6
- 230000000737 periodic effect Effects 0.000 claims abstract description 4
- 230000009467 reduction Effects 0.000 claims abstract description 3
- 238000004364 calculation method Methods 0.000 claims description 10
- 230000008569 process Effects 0.000 claims description 7
- 230000009471 action Effects 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 5
- 238000011896 sensitive detection Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000009434 installation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- 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
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种用于随钻方位电磁波电阻率仪器的幅相提取方法及装置,所述方法包括如下步骤:步骤1、将接收天线感应的微弱电信号利用信号调理电路进行程控增益放大和低频滤波处理,得到满足采样要求的待采样模拟信号;步骤2、根据奈奎斯特采样定理,利用模数转换模块将所述待采样模拟信号进行离散化处理,得到多个周期的数字信号;步骤3、控制器对离散后的连续多个周期信号进行叠加降噪处理,将每个周期内相位相同的采样点进行叠加,叠加后取平均得到一个周期的数字信号,用此信号与方波参考信号进行幅度和相位的提取;所述的方波参考信号由每周期采样点数N决定,再根据接收信号的同相分量和正交分量,进而计算出被测信号幅度和相位信息。
Description
技术领域
本发明涉及油气勘探测量电路技术领域,尤其是一种用于随钻方位电磁波电阻率仪器的幅相提取方法及装置。
背景技术
随着常规油气勘探开发进入中后期,越来越多的复杂油气资源需要采用大斜度井或水平井开发以提高产量和效益,其中随钻测井技术是大斜度井、水平井油气勘探开发不可缺少的重要手段,主要用于地质导向和地层评价,成为当今石油勘探开发中不可替代的技术之一。该技术是在钻井时获取岩性、饱和度、孔隙度等地层评价信息,能提供钻井工程参数以便于有效的钻进储层,并能在恶劣井眼条件下开展测井工作,因而在诸多方面具有优越性。随钻方位电磁波电阻率仪器作为地质导向和地层评价的重要仪器之一,其利用地层电磁学特性测量电阻率信息。所测量的地层电阻率是开展地层含油、气、水等定性评价的重要依据,同时也是定量评价储层含油气饱和度的重要参数之一。该仪器具有高效、准确等特点,能有效对地层未知的探测区域进行评估,从而保证井眼轨迹在期望的储层内,提高储层钻遇率。
传统的随钻电磁波电阻率仪器由于发射天线与接收天线处于同一方向,不具备方位探测能力,无法探测地层边界和走向。而随钻方位电磁波电阻率仪器通过多个发射天线向地层发射出不同频率的电磁波信号,由于不同电阻率地层对电磁波的吸收作用不同,经地层传播、衰减后的信号会携带地层参数信息,通过改变接收天线的倾角和安装位置,对接收到的电磁波信号进行处理和计算,得到不同接收天线的幅度和相位信息,进而反演出地层电阻率、到地层边界的距离和方位,用于实时地质导向和地层评价。
高效提取接收天线感应的电磁波信号的幅度和相位是获取地层电阻率的关键问题。目前商用仪器通常将接收到的高频模拟信号进行程控增益放大、低通滤波、模数转换等处理,得到数字信号,再提取出数字信号的幅度和相位。其中,幅度和相位的提取主要采用数字相敏检波算法,将采集的数字信号与预先存储在仪器内部的与被测信号同频的标准正余弦信号进行乘法和加法运算,得到被测信号的同相分量和正交分量,最终计算出幅度和相位。该方法利用接收信号中的有用信号与参考信号之间的相关性来检测出幅度和相位,可以抑制信号中的噪声,但此种计算方法相对复杂,占用较多的井下处理器资源,降低了处理的实时性、提高了功耗,给仪器的应用带来不便。
在对接收信号进行程控增益放大、低通滤波和模数转换的过程中,由于FPGA(Field Programmable Gate Array,现场可编程门阵列)具有更好的并行处理能力,其可以实现对程控增益放大电路、低通滤波电路和模数转换电路等外围硬件电路的精确控制,同时也可以对采集的数字信号进行处理和计算,因此被优先选用作为随钻方位电磁波电阻率仪器的井下控制器。但采用上述数字相敏检波算法,利用标准正余弦信号作为提取幅度和相位的参考信号时,需要将其预先存储在FPGA内部的ROM(Read-Only Memory)中,在计算过程中,需要调用FPGA内部的乘法器,导致内部电路过于复杂,占用更多的处理器资源,增大电路功耗。
发明内容
本发明的目的在于提供一种随钻方位电磁波电阻率仪器的幅相提取方法及装置,简化传统利用正余弦信号作为参考信号的数字相敏检波计算方法,降低井下处理器资源使用率。从而降低电路功耗,有利于在井下仪器中实现。
本发明的技术方案为:一种用于随钻方位电磁波电阻率仪器的幅相提取方法,包括如下步骤:
步骤1、将接收天线感应的微弱电压信号利用信号调理电路进行程控增益放大和低频滤波处理,得到满足采样要求的待采样模拟信号;
步骤2、根据奈奎斯特采样定理,利用模数转换模块将所述待采样模拟信号进行离散化处理,得到多个周期的数字信号;
步骤3、控制器对离散后的连续多个周期信号进行叠加降噪处理,将每个周期内相位相同的采样点进行叠加,叠加后取平均得到一个周期的数字信号,用此信号与方波参考信号进行幅度和相位的提取;所述的方波参考信号由每周期采样点数决定,再根据接收信号的同相分量和正交分量,进而计算出被测信号的幅度和相位信息。
有益效果:
现有的随钻方位电磁波电阻率仪器利用被测信号与同频标准正余弦信号进行数字相敏检波运算,提取出幅度和相位信息,该方法需预先存储好离散的正余弦信号以及需要调用FPGA内部的乘法器资源,增加了算法和电路的复杂性。本技术提出的方波参考信号避免使用ROM和乘法器,只利用简单的加减法运算,即可提取幅度和相位信息,简化了计算方式,提高了可靠性,更适用于井下高温的恶劣环境。提高幅度和相位提取计算的实时性,进而提高仪器实时地质导向能力。
附图说明
图1为幅度和相位的处理流程和方法;
图2为随钻方位电磁波电阻率仪器信号处理流程。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整的描述,显然,所描述的实施例仅为本发明的一部分实施例,而不是全部的实施例,基于本发明中的实施例,本领域的普通技术人员在不付出创造性劳动的前提下所获得的所有其他实施例,都属于本发明的保护范围。
根据本发明实施例,一种用于随钻方位电磁波电阻率仪器的幅相提取装置,包括:接收天线、信号调理模块、模数转换模块和控制器。
如图1中,所述接收天线用于将接收天线感应到的微弱电磁波信号转换为电压信号,接收天线的输出端与信号调理模块相连,供后续电路处理。
所述信号调理模块,与控制器和模数转换模块相连,共包括两部分:程控增益放大电路和低通滤波电路。程控增益放大电路接收控制器的增益设置参数,实现根据微弱信号的幅度灵活调整增益的放大功能;低通滤波电路根据接收信号的频率,在控制器的作用下选择适当的滤波频率,滤除高频噪声。
所述模数转换模块,与信号调理模块和控制器相连,用于将满足采样要求的模拟信号进行离散化处理,得到数字信号,供后续电路处理和计算。
所述控制器选用FPGA,分别与信号调理模块和模数转换模块相连。一方面实现对信号调理模块中增益的控制以及滤波频率的选择;另一方面实现对模数转换模块中ADC(Analog-to-Digital Converter)芯片采样频率的配置,将模拟信号转换为数字信号,并对采样后的数字信号进行处理和计算,利用与每周期采样点数相同的方波信号实现被测信号幅度和相位的提取。
所述幅度和相位的提取方法,是将离散后多个周期的数字信号进行叠加,将相位相同的采样点进行叠加求和取平均,最终将多个周期的信号处理成一个周期的信号,并将这一个周期的信号与点数相同的离散的方波信号进行运算,得到被测信号的幅度和相位。
具体地,求取被测信号幅度和相位的方法包括以下步骤:
a.根据奈奎斯特采样定理,对频率为f0的被测信号进行离散处理。
具体地,控制器FPGA配置ADC芯片的采样频率为fs,与被测信号之间满足关系:fs=N×f0,其中,N为一个周期内的采样点数,且满足,N>2,f0为被测信号频率,离散后的数字信号可以用x[k]表示,k取0,1,2...m-1,m为采样周期数。
c.对叠加后的信号y[n]与同频率采样点数相同的离散的方波信号r[n]和c[n]进行运算。方波信号如下所示:
具体地,利用y[n]、r[n]和c[n]进行如下运算:
根据本发明的一个具体实施例,以随钻方位电磁波电阻率仪器上接收天线实际接收到的电磁波信号为例,进一步阐述本发明提出的技术方案。
按所述技术方案中的步骤进行幅度和相位的提取,具体操作如下:
步骤a:如附图2所示,随钻方位电磁波电阻率仪器首先利用接收天线将感应到的电磁波信号转换为电压信号,此时电压信号极其微弱,仅在nV级别。
步骤b:对微弱的电信号进行放大、滤波处理。由于随钻方位电磁波电阻率仪器接收天线感应到的信号极其微弱且有多个信号频率,通常在400kHz或2MHz附近,因此步骤b需要在控制器FPGA的作用下结合信号幅度和频率实现不同的增益控制和滤波频率的选择,从而达到ADC芯片的采样要求。
步骤c:根据奈奎斯特采样定理,采样频率应大于被测信号最高频率的两倍,才能无失真地恢复出被测信号,而且为了保证采集精度,实际选用的采样频率通常远高于被测信号最高频率的两倍,此处将采样频率设置为待测信号最高频率的8倍,即二者之间满足关系:fs=8×f0,每周期采样8个点。
假设被测的模拟信号可以表示为:
离散后的多个周期可以表示为:
其中,k取0,1,2...m-1,m为采样周期数。
再根据N=8,确定方波信号:
则以方波提取幅度和相位的方法可以表示为:
进而得到幅度和相位的计算公式:
以上就是随钻方位电磁波电阻率仪器通过采用与每周期采样数相同的方波信号为参考提取接收天线上信号幅度和相位的计算方法。一方面,该方法不用预先将标准的正余弦信号存储在控制器FPGA内的ROM中,且在计算的过程中,只需简单的加减法运算,即可计算出幅度和相位信息,避免使用乘法器资源,降低了井下处理器资源利用率,简化电路结构,提升了仪器的可靠性,使之更加适应井下高温的测试环境。另一方面,可以实现边采集边计算的方式,节约了井下测量时间,这可以使测井仪器及时、高效地判断地层信息,提高了数据处理的实时性,保证仪器在储层内钻进,提高钻遇率。
尽管上面对本发明说明性的具体实施方式进行了描述,以便于本技术领域的技术人员理解本发明,且应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
Claims (6)
1.一种用于随钻方位电磁波电阻率仪器的幅相提取方法,其特征在于,包括如下步骤:
步骤1、将接收天线感应的微弱电压信号利用信号调理电路进行程控增益放大和低频滤波处理,得到满足采样要求的待采样模拟信号;
步骤2、根据奈奎斯特采样定理,利用模数转换模块将所述待采样模拟信号进行离散化处理,得到多个周期的数字信号;
步骤3、控制器对离散后的连续多个周期信号进行叠加降噪处理,将每个周期内相位相同的采样点进行叠加,叠加后取平均得到一个周期的数字信号,用此信号与方波参考信号进行幅度和相位的提取;所述的方波参考信号由每周期采样点数决定,再根据接收信号的同相分量和正交分量,进而计算出被测信号的幅度和相位信息。
6.一种用于随钻方位电磁波电阻率仪器的幅相提取装置,其特征在于,所述装置包括:
接收天线,用于将接收天线感应到的微弱电磁波信号转换为电压信号,接收天线的输出端与信号调理模块相连,供后续信号调理模块电路处理;
信号调理模块,与控制器和模数转换模块相连,共包括两部分:程控增益放大电路和低通滤波电路;程控增益放大电路在控制器的作用下,接收参数信息,实现对微弱信号指定增益的放大功能;低通滤波电路根据接收信号的频率,在控制器的作用下选择适当的滤波频率,滤除高频噪声;
模数转换模块,与信号调理模块和控制器相连,用于将满足采样要求的模拟信号进行离散化处理,得到数字信号,供后续电路处理和计算。
控制器,选用FPGA,分别与信号调理模块和模数转换模块相连,实现对信号调理模块中增益的控制以及滤波频率的选择;以及对模数转换模块中ADC芯片采样频率进行配置,将模拟信号转换为数字信号,并对离散后的数字信号进行处理和计算,利用与每周期采样点数相同的方波信号实现待测信号幅度和相位的提取;
所述幅度和相位的提取,是将离散后多个周期的数字信号进行叠加,将相位相同的采样点进行叠加求和取平均,最终将多个周期的信号处理成一个周期的信号,并将这一个周期的信号与点数相同的离散的方波信号进行运算,得到被测信号的幅度和相位。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211355755.4A CN115685355B (zh) | 2022-11-01 | 2022-11-01 | 用于随钻方位电磁波电阻率仪器的幅相提取方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211355755.4A CN115685355B (zh) | 2022-11-01 | 2022-11-01 | 用于随钻方位电磁波电阻率仪器的幅相提取方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115685355A true CN115685355A (zh) | 2023-02-03 |
CN115685355B CN115685355B (zh) | 2024-02-27 |
Family
ID=85047982
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211355755.4A Active CN115685355B (zh) | 2022-11-01 | 2022-11-01 | 用于随钻方位电磁波电阻率仪器的幅相提取方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115685355B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116856920A (zh) * | 2023-07-06 | 2023-10-10 | 中国科学院地质与地球物理研究所 | 一种随钻方位电磁波电阻率仪器使用方法及仪器 |
CN117270063A (zh) * | 2023-09-13 | 2023-12-22 | 中国科学院地质与地球物理研究所 | 一种信号频率获取方法以及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102704921A (zh) * | 2012-05-28 | 2012-10-03 | 中国石油天然气集团公司 | 一种随钻电磁波电阻率的测量装置和测量方法 |
CN105676008A (zh) * | 2016-01-15 | 2016-06-15 | 中国电力科学研究院 | 一种数字式电场传感器 |
CN106374918A (zh) * | 2016-08-29 | 2017-02-01 | 北京航空航天大学 | 一种基于dsp和fpga的高速数字锁相放大器 |
CN107152276A (zh) * | 2016-03-02 | 2017-09-12 | 中国石油化工股份有限公司 | 用于随钻电磁波电阻率测量的信号处理方法与系统 |
CN112160746A (zh) * | 2020-09-27 | 2021-01-01 | 电子科技大学 | 一种超深电阻率测井的时域测量装置 |
WO2022068635A1 (zh) * | 2020-09-30 | 2022-04-07 | 华为技术有限公司 | 一种信号处理电路以及相关装置 |
-
2022
- 2022-11-01 CN CN202211355755.4A patent/CN115685355B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102704921A (zh) * | 2012-05-28 | 2012-10-03 | 中国石油天然气集团公司 | 一种随钻电磁波电阻率的测量装置和测量方法 |
CN105676008A (zh) * | 2016-01-15 | 2016-06-15 | 中国电力科学研究院 | 一种数字式电场传感器 |
CN107152276A (zh) * | 2016-03-02 | 2017-09-12 | 中国石油化工股份有限公司 | 用于随钻电磁波电阻率测量的信号处理方法与系统 |
CN106374918A (zh) * | 2016-08-29 | 2017-02-01 | 北京航空航天大学 | 一种基于dsp和fpga的高速数字锁相放大器 |
CN112160746A (zh) * | 2020-09-27 | 2021-01-01 | 电子科技大学 | 一种超深电阻率测井的时域测量装置 |
WO2022068635A1 (zh) * | 2020-09-30 | 2022-04-07 | 华为技术有限公司 | 一种信号处理电路以及相关装置 |
Non-Patent Citations (7)
Title |
---|
IQBAL SHAHID: "Power modulation of vortex beams using phase/amplitude adjustable transmissive coding metasurfaces", JOURNAL OF PHYSICS D-APPLIED PHYSICS, vol. 54, no. 3, 21 January 2021 (2021-01-21) * |
和丽真;杨国华;王?;郭庆明;: "基于随钻电磁波仪器的欠采样DPSD设计与实现", 电子测量技术, no. 012 * |
和丽真;杨国华;王?;郭庆明;: "基于随钻电磁波仪器的欠采样DPSD设计与实现", 电子测量技术, no. 012, 31 December 2020 (2020-12-31) * |
和丽真;杨国华;王?;郭庆明;: "基于随钻电磁波仪器的欠采样DPSD设计与实现", 电子测量技术, no. 12, 23 June 2020 (2020-06-23) * |
张金波;湛向;刘二伟;曹爱华;: "基于DDS/SOPC的多路可调谐波信号发生器", 电力自动化设备, no. 11 * |
牛明: "基于脉冲成型的OFDM系统峰均比抑制", 自动化技术与应用, no. 003 * |
王中兴;荣亮亮;林君;尚新磊;段清明;蒋川东;: "基于4倍频采样的数字正交FID信号检测技术", 数据采集与处理, no. 05 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116856920A (zh) * | 2023-07-06 | 2023-10-10 | 中国科学院地质与地球物理研究所 | 一种随钻方位电磁波电阻率仪器使用方法及仪器 |
CN116856920B (zh) * | 2023-07-06 | 2024-04-02 | 中国科学院地质与地球物理研究所 | 一种随钻方位电磁波电阻率仪器使用方法及仪器 |
CN117270063A (zh) * | 2023-09-13 | 2023-12-22 | 中国科学院地质与地球物理研究所 | 一种信号频率获取方法以及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN115685355B (zh) | 2024-02-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115685355B (zh) | 用于随钻方位电磁波电阻率仪器的幅相提取方法及装置 | |
US6777940B2 (en) | Apparatus and method for resistivity well logging | |
CA2706861C (en) | Determination of borehole azimuth and the azimuthal dependence of borehole parameters | |
US7260477B2 (en) | Estimation of borehole geometry parameters and lateral tool displacements | |
US20060173627A1 (en) | Enhanced measurement of azimuthal dependence of subterranean parameters with filters and/or discretely sampled data | |
CN104343440B (zh) | 泥浆压力脉冲信号的检测方法及系统 | |
CN109143374B (zh) | 一种井周散射体成像方法及系统 | |
WO2010129780A2 (en) | Directional resistivity imaging using harmonic representations | |
WO2004070422A2 (en) | Methods of processing magnetotelluric signals | |
US10901109B2 (en) | Detection and avoidance of downhole noise and interference sources | |
US6414492B1 (en) | Method and apparatus for passive detection of geophysical discontinuities in the earth | |
US20240151872A1 (en) | Signal acquisition method and device of azimuthal electromagnetic wave resistivity instrument while drilling | |
CN109493637A (zh) | 一种车位状态检测方法 | |
RU2628660C2 (ru) | Измерения дальности с использованием модулированных сигналов | |
US7855560B2 (en) | Analog front-end coherent in-phase demodulation data acquisition system for resistivity image tools | |
CN110133727A (zh) | 超深油气储层高精度电磁频谱探测方法 | |
US7863902B2 (en) | Coherent in-phase demodulation for resistivity image tools with spread spectrum continuous phase modulation stimulus | |
CN106254007B (zh) | 一种应用于lte230系统的频谱感知方法和装置 | |
Liu et al. | Research on weak signal detection for downhole acoustic telemetry system | |
CN111611686A (zh) | 一种通信信号时频域的检测方法 | |
CN105008950A (zh) | 用于测向和方向不确定度确定的测量装置和测量方法 | |
WO2007114923A2 (en) | Method for detection of signal source using estimation of noise statistics | |
CN108071385B (zh) | 一种随钻电磁波电阻率幅值比和相位差的测量方法和电路 | |
CN109061629B (zh) | 一种基于地质雷达技术的泥石流堆积物厚度探测方法 | |
CN112177601B (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 |