CN112901146A - 声波测井探测中坏道的识别方法及装置 - Google Patents
声波测井探测中坏道的识别方法及装置 Download PDFInfo
- Publication number
- CN112901146A CN112901146A CN201911117129.XA CN201911117129A CN112901146A CN 112901146 A CN112901146 A CN 112901146A CN 201911117129 A CN201911117129 A CN 201911117129A CN 112901146 A CN112901146 A CN 112901146A
- Authority
- CN
- China
- Prior art keywords
- time
- waveform
- bad track
- identifying
- logging
- 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
- 238000000034 method Methods 0.000 title claims abstract description 72
- 238000001514 detection method Methods 0.000 title claims abstract description 43
- 238000012545 processing Methods 0.000 claims description 29
- 238000004590 computer program Methods 0.000 claims description 14
- 230000010355 oscillation Effects 0.000 claims description 14
- 238000004458 analytical method Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 208000011580 syndromic disease Diseases 0.000 claims description 3
- 239000000523 sample Substances 0.000 claims 6
- 238000003384 imaging method Methods 0.000 abstract description 20
- 230000006870 function Effects 0.000 description 14
- 238000010586 diagram Methods 0.000 description 13
- 238000004891 communication Methods 0.000 description 12
- 238000012937 correction Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 5
- 239000000872 buffer Substances 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000005404 monopole Effects 0.000 description 3
- 230000010356 wave oscillation Effects 0.000 description 3
- 238000005553 drilling Methods 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000010295 mobile communication Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- 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
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
本发明提供一种声波测井探测中坏道的识别方法及装置，方法包括：获取目的层段的测井数据；对所述测井数据进行解析得到直达波的到时曲线；以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道。本发明能够有效提高测井数据中坏道的识别准确率和快速性，进而提高测井数据的反射波信噪比和反射体成像精度。
Description
技术领域
本发明涉及油气勘探技术领域，具体涉及一种声波测井探测中坏道的识别方法及装置。
背景技术
脆性地层垮塌、泥岩地层膨胀、钻井速度过快等原因均会使得测井坏境变得恶劣，产生仪器磕碰井壁的现象，对声波测井的波形产生严重影响，导致深度点的波形出现紊乱，变为无效数据，通常这些出现紊乱的波形被称为坏道。坏道会严重干扰声波测井资料中有效信息的提取。因此，识别与校正坏道是开展声波测井处理的关键步骤，尤其是对波形信噪比要求非常高的远探测声波测井数据处理。
在地震勘探领域，地震资料采集环境破坏或地震采集仪器故障等原因均会使得采集数据出现坏道。特别对于海上地震勘探，钻井平台噪声、过船干扰等因素也会给地震采集数据引入噪声，使得一些地震道成为坏道。现有技术中，根据坏道形成原因及波形特征，提出了一系列坏道识别与校正方法，如初至起跳判断法、50Hz干扰道识别法、主频异常识别法和极性反转判断法等。但是，这些方法都是针对地震采集过程中形成的坏道提出的识别方法，由于声波测井波形采集环境及波形特征与地震勘探采集资料的巨大差异，使用现有的坏道识别方法识别声波测井数据中的坏道时，会导致识别准确性降低的问题，进而降低了基于测井数据的反射体成像精度。
因此亟需一种适用于声波测井波形的坏道识别方法。
发明内容
针对现有技术中的问题，本发明提供一种声波测井探测中坏道的识别方法及装置，能够提高坏道的识别准确率和快速性，进而提高声波测井数据的反射波信噪比和反射体成像精度。
为解决上述技术问题，本发明提供以下技术方案：
第一方面，本发明提供一种声波测井探测中坏道的识别方法，包括：
获取目的层段的测井数据；
对所述测井数据进行解析得到直达波的到时曲线；
以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；
根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道。
进一步的，在所述获取目的层段的测井数据并对该测井数据进行解析得到直达波的到时曲线之后，还包括：
对所述测井数据进行解析得到直达波的震荡结束时间；
以所述震荡结束时间至整道波形终止时间为第二时间段，确定所述第二时间段内目的层段中各个深度点的波形幅度值；
根据所述第二时间段内目的层段中各个深度点的波形幅度值和预设的第二临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第二临界值，则确定该深度点的波形为坏道。
进一步的，还包括：
根据所述坏道相邻的非坏道对所述坏道进行校正。
其中，所述根据所述坏道相邻的非坏道对所述坏道进行校正，包括：
将所述坏道中的波形替换为相邻的非坏道中的波形。
进一步的，在所述获取目的层段的测井数据之前，还包括：
通过阵列声波测井的方式获取目标区块的测井数据；
其中，所述目标区块包括目的层段。
其中，所述对所述测井数据进行解析得到直达波的到时曲线，包括：
对所述测井数据进行转换处理得到声波的时差曲线；
根据所述时差曲线确定到时曲线。
其中，所述对所述测井数据进行转换处理得到声波的时差曲线，包括：
采用慢度-时间相关的方式对所述测井数据进行转换处理得到声波的时差曲线。
其中，所述第一临界值是根据波形幅度曲线平均值确定的。
第二方面，本发明提供一种声波测井探测中坏道的识别装置，包括：
获取单元，用于获取目的层段的测井数据；
第一解析单元，用于对所述测井数据进行解析得到直达波的到时曲线；
第一处理单元，用于以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；
第一识别单元，用于根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道。
进一步的，还包括：
第二解析单元，用于对所述测井数据进行解析得到直达波的震荡结束时间；
第二处理单元，用于以所述震荡结束时间至整道波形终止时间为第二时间段，确定所述第二时间段内目的层段中各个深度点的波形幅度值；
第二识别单元，用于根据所述第二时间段内目的层段中各个深度点的波形幅度值和预设的第二临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第二临界值，则确定该深度点的波形为坏道。
进一步的，还包括：
校正单元，用于根据所述坏道相邻的非坏道对所述坏道进行校正。
其中，所述校正单元，包括：
校正子单元，用于将所述坏道中的波形替换为相邻的非坏道中的波形。
进一步的，还包括：
采集单元，用于通过阵列声波测井的方式获取目标区块的测井数据；
其中，所述目标区块包括目的层段。
其中，所述第一解析单元，包括：
第一处理子单元，用于对所述测井数据进行转换处理得到声波的时差曲线；
第二处理子单元，用于根据所述时差曲线确定到时曲线。
其中，所述第一处理子单元，包括：
转换模块，用于采用慢度-时间相关的方式对所述测井数据进行转换处理得到声波的时差曲线。
其中，所述第一临界值是根据波形幅度曲线平均值确定的。
第三方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的声波测井探测中坏道的识别方法的步骤。
第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述的声波测井探测中坏道的识别方法的步骤。
由上述技术方案可知，本发明提供一种声波测井探测中坏道的识别方法及装置，通过获取目的层段的测井数据；对所述测井数据进行解析得到直达波的到时曲线；以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道，能够提高测井数据中坏道的识别准确率和快速性，进而提高测井数据的反射波信噪比和反射体成像精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中声波测井探测中坏道的识别方法的第一种流程示意图。
图2为本发明实施例中声波测井探测中坏道的识别方法的第二种流程示意图。
图3为本发明实施例中声波测井探测中坏道的识别方法的第三种流程示意图。
图4为本发明实施例中声波测井探测中坏道的识别方法的第四种流程示意图。
图5是X1井根据1号时间窗波形幅度剔除坏道示意图。
图6是X1井根据2号时间窗波形幅度剔除坏道及校正坏道示意图。
图7是X1井坏道校正前后偏移成像结果对比图。
图8为本发明实施例中的声波测井探测中坏道的识别装置的结构示意图。
图9为本发明实施例中的电子设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明提供一种声波测井探测中坏道的识别方法的实施例，参见图1，所述声波测井探测中坏道的识别方法具体包含有如下内容：
S101：获取目的层段的测井数据；
可以理解的是，可以通过常规测井、电成像测井和阵列声波测井中至少一种方式获取测井数据。
常规测井方式的测井数据主要包括九条常规测井曲线：自然伽马曲线、井径曲线、自然电位曲线、纵波时差曲线、补偿中子曲线、密度曲线、微电极电阻率曲线、深电阻率曲线和浅电阻率曲线。在测井过程中可直接获得这九条曲线。
通过处理电成像测井方式的测井数据，能够得到电阻率在井壁的周向成像图，用于指示井旁裂缝洞穴等反射体的发育程度和方位。
阵列声波测井方式的测井数据包括单极阵列波形和偶极阵列波形资料，通过处理阵列声波测井方式的测井数据，能够获取相应的单极纵波远探测偏移成像结果和偶极横波远探测偏移成像结果。
需要说明的是，目的层段的测井数据为通过上述方式获取的测井数据中的被测井中的目的层段所对应的数据。
本实施例中的测井数据是通过阵列声波测井得到的。
S102：对所述测井数据进行解析得到直达波的到时曲线；
可以理解的是，到时曲线为各个地震波的到达时间组成的曲线。
在本步骤中，对测井数据进行解析得到直达波的到时曲线，具体是对测井数据进行转换处理得到声波的时差曲线，其中，采用慢度-时间相关的方式处理测井数据得到首波的时差曲线，根据得到的时差曲线确定到时曲线。
根据得到的时差曲线确定到时曲线的计算方式如下：
其中，Time为到时曲线中的到达时间，单位为μs；Slowness为时差曲线中的时差，单位为μs/ft；sBg为起始计算深度，单位为ft；sEnd为终止计算深度，单位为ft；终止计算深度sEnd和起始计算深度sBg之差为仪器源距；ds代表深度采样点间隔，单位为ft。
S103：以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；
在本步骤中，以测井数据中的零时刻为起始时间，以到时曲线中的最小值时刻为截止时间，将该时间段设置为第一时间段，确定第一时间段内目的层段中各个深度点的波形幅度值，其中，波形幅度值采用下式进行计算：
其中，Amp1为第一时间段内波形幅度值，单位为mV；TimeMin为到时曲线中的最小值时刻，单位为μs；Amp1Wv(t)为第一时间段内t时刻的波形幅度，单位为mV；dt为波形时间采样间隔，单位为μs；Num1为第一时间段内波形采样点数。
S104：根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道。
在本步骤中，判断第一临界值分别与各个深度点的波形幅度值的大小，确定波形幅度值大于第一临界值的深度点，将该确定的深度点对应的波形为坏道。
其中，第一临界值是根据波形幅度曲线平均值确定的，具体的，第一临界值设定为第一时间段内对应的波形幅度曲线平均值的多倍值，公式如下：
TAmp1＝Coff1·AmpAvg1；
其中，TAmp1为第一时间段内对应的第一临界值，单位为mV；AmpAvg1为第一时间段内对应的波形幅度曲线平均值，单位为mV；Coff1为系数，本实施例中Coff1＝2。
从上述描述可知，本发明实施例提供的声波测井探测中坏道的识别方法，通过获取目的层段的测井数据；对所述测井数据进行解析得到直达波的到时曲线；以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道，能够提高测井数据中坏道的识别准确率和快速性，进而提高测井数据的反射波信噪比和反射体成像精度。
在本发明的一实施例中，参见图2，所述声波测井探测中坏道的识别方法的步骤S101之后还包含有步骤S202、步骤S203和步骤S204，具体包含有如下内容：
S202：对所述测井数据进行解析得到直达波的震荡结束时间；
在本步骤中，对测井数据进行解析，并统计直达波的震荡结束时间。本实施例中的测井数据为阵列声波测井数据，对于阵列声波测井数据中的单极阵列波形，直达波震荡结束时间为7ms；对于阵列声波测井数据中的偶极阵列波形，直达波震荡结束时间为8ms；通常结束时间不直接选择整道波形终止时间，而是选择终止时间之前0.5ms。
S203：以所述震荡结束时间至整道波形终止时间为第二时间段，确定所述第二时间段内目的层段中各个深度点的波形幅度值；
在本步骤中，以测井数据中的震荡结束时间为起始时间，以整道波形终止时间为截止时间，将该时间段设置为第二时间段，确定第二时间段内目的层段中各个深度点的波形幅度值，其中，波形幅度值采用下式进行计算：
其中，Amp2为第二时间段内波形幅度值，单位为mV；TimeBg为第二时间段内的起始时间，单位为μs；TimeEnd为第二时间段内的结束时间，单位为μs；Amp2Wv(t)为第二时间段内t时刻的波形幅度，单位为mV；dt为波形时间采样间隔，单位为μs；Num2为第一时间段内波形采样点数。
S204：根据所述第二时间段内目的层段中各个深度点的波形幅度值和预设的第二临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第二临界值，则确定该深度点的波形为坏道。
在本步骤中，判断第二临界值分别与各个深度点的波形幅度值的大小，确定波形幅度值大于第二临界值的深度点，将该确定的深度点对应的波形为坏道。
其中，第二临界值是根据波形幅度曲线平均值确定的，具体的，第二临界值设定为第二时间段内对应的波形幅度曲线平均值的多倍值，公式如下：
TAmp2＝Coff2·AmpAvg2；
其中，TAmp2为第一时间段内对应的第一临界值，单位为mV；AmpAvg2为第一时间段内对应的波形幅度曲线平均值，单位为mV；Coff2为系数，本实施例中Coff1＝3。
在本发明的一实施例中，参见图3，所述声波测井探测中坏道的识别方法中后还包含有步骤S301，具体包含有如下内容：
S301：根据所述坏道相邻的非坏道对所述坏道进行校正。
在本步骤中，采用坏道相邻的非坏道对坏道进行校正，能够提高测井数据的可用性，还能够提高测井数据的反射波信噪比和反射体成像精度。
在对坏道进行校正时，具体是将坏道中的波形替换为相邻的非坏道中的波形。通常选择坏道上部相邻的非坏道来替换坏道波形。若连续出现多道坏道，则将最上端坏道相邻的非坏道对坏道实施全部替换。
在本发明的一实施例中，参见图4，所述声波测井探测中坏道的识别方法的步骤S101之前还包含有步骤S100，具体包含有如下内容：
S100：通过阵列声波测井的方式获取目标区块的测井数据；
其中，所述目标区块包括目的层段。
可以理解的是，可以通过常规测井、电成像测井和阵列声波测井中至少一种方式获取测井数据。
需要说明的是，目的层段的测井数据为通过上述方式获取的测井数据中的被测井中的目的层段所对应的数据。
在本步骤中，通过阵列声波测井的方式获取目标区块的测井数据，该目标区块包括目的层段，进而实现获取目的层段的测井数据。
为进一步地说明本方案，本发明提供一种声波测井探测中坏道的识别方法的具体实例，具体包含有如下内容：
S1：搜集目标区块相关测井数据，包括常规测井、电成像测井及阵列声波测井等资料。
参见图5中第二道包含阵列声波测井资料中偶极原始波形数据，从图5中可看到一些明显的坏道，比如在深度7930m和7950m附近。
S2：提取目的层段阵列声波测井资料中直达波到时曲线，并统计到时最小值。
S3：以零时刻为起始时间，以直达波到时最小值为截止时间，开1号时间窗，计算窗内波形幅度值。
其中，1号时间窗的起始时间与终止时间与第一时间段的起始时间与终止时间相同。
参见图5中第四道包含1号时间窗对应的波形幅度曲线，坏道所在深度位置的波形幅度出现峰值。
S4：根据1号时间窗对应的波形幅度曲线设定幅度临界值，对比目的层段内各深度点波形幅度值与临界值大小关系，若超出临界值，则判断该深度点波形为坏道。
X1井1号时间窗对应的幅度临界值设定为波形幅度均值的2倍，可以看到波形幅度曲线(图5第四道)中的坏道对应的峰值会明显超出幅度临界值。图5中第五道为1号时间窗对应的坏道统计曲线，曲线中数值1代表非坏道，数值2代表坏道。图5中第三道为根据1号时间窗波形幅度剔除坏道后的波形。
S5：统计直达波震荡结束时间，以该时间为起始时间，整道波形终止时间为结束时间开2号时间窗，计算窗内波形幅度值。
其中，2号时间窗的起始时间与终止时间与第二时间段的起始时间与终止时间相同。
参见图6中第四道为2号时间窗对应的波形幅度曲线，坏道所在深度位置的波形幅度出现峰值。
S6：根据2号时间窗对应的波形幅度曲线设定能量临界值，对比各深度点波形幅度值与临界值大小关系，若超出临界值，则判断该深度点波形为坏道。
X1井2号时间窗对应的幅度临界值设定为波形幅度均值的3倍，可以看到波形幅度曲线(图6第四道)中的坏道对应的峰值会明显超出幅度临界值。图6中第五道为2号时间窗对应的坏道统计曲线。图6中第三道为根据2号时间窗波形幅度剔除坏道后的波形，这是在根据1号时间窗波形幅度剔除坏道的基础上实施的。
S7：将坏道相邻的非坏道中波形替换坏道波形，实现坏道校正。
从上述描述可知，图6中第六道显示了坏道校正后的波形，与图5第一道中的原始波形相比，波形中一些明显的坏道得到了校正，比如深度7930m附近坏道。图7中第二道为没有实施坏道剔除与校正过程，采用常规远探测处理流程得到的声波远探测偏移成像结果，在深度7920m、7930m、7950m以及7980m可观察到明显的由坏道引起的反射体假象。图7中第三道为实施坏道剔除与校正之后，采用上述实施例中的识别方法和校正方法处理后得到的偏移成像结果，可以看到基本消除了坏道引起的反射体假象，有效反射体(圆框内)成像精度明显提升。
本发明实施例提供一种能够实现所述声波测井探测中坏道的识别方法中全部内容的声波测井探测中坏道的识别装置的具体实施方式，参见图8，所述声波测井探测中坏道的识别装置具体包括如下内容：
获取单元20，用于获取目的层段的测井数据；
第一解析单元30，用于对所述测井数据进行解析得到直达波的到时曲线；
第一处理单元40，用于以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；
第一识别单元50，用于根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道。
进一步的，还包括：
第二解析单元60，用于对所述测井数据进行解析得到直达波的震荡结束时间；
第二处理单元70，用于以所述震荡结束时间至整道波形终止时间为第二时间段，确定所述第二时间段内目的层段中各个深度点的波形幅度值；
第二识别单元80，用于根据所述第二时间段内目的层段中各个深度点的波形幅度值和预设的第二临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第二临界值，则确定该深度点的波形为坏道。
进一步的，还包括：
校正单元90，用于根据所述坏道相邻的非坏道对所述坏道进行校正。
其中，所述校正单元，包括：
校正子单元，用于将所述坏道中的波形替换为相邻的非坏道中的波形。
进一步的，还包括：
采集单元10，用于通过阵列声波测井的方式获取目标区块的测井数据；
其中，所述目标区块包括目的层段。
其中，所述第一解析单元，包括：
第一处理子单元，用于对所述测井数据进行转换处理得到声波的时差曲线；
第二处理子单元，用于根据所述时差曲线确定到时曲线。
其中，所述第一处理子单元，包括：
转换模块，用于采用慢度-时间相关的方式对所述测井数据进行转换处理得到声波的时差曲线。
其中，所述第一临界值是根据波形幅度曲线平均值确定的。
本发明提供的声波测井探测中坏道的识别装置的实施例具体可以用于执行上述实施例中的声波测井探测中坏道的识别方法的实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。
从上述描述可知，本发明实施例提供的声波测井探测中坏道的识别装置，通过获取目的层段的测井数据；对所述测井数据进行解析得到直达波的到时曲线；以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道，能够提高测井数据中坏道的识别准确率和快速性，进而提高测井数据的反射波信噪比和反射体成像精度。
本申请提供一种用于实现所述声波测井探测中坏道的识别方法中的全部或部分内容的电子设备的实施例所述电子设备具体包含有如下内容：
处理器(processor)、存储器(memory)、通信接口(Communications Interface)和总线；其中，所述处理器、存储器、通信接口通过所述总线完成相互间的通信；所述通信接口用于实现相关设备之间的信息传输；该电子设备可以是台式计算机、平板电脑及移动终端等，本实施例不限于此。在本实施例中，该电子设备可以参照实施例用于实现所述声波测井探测中坏道的识别方法的实施例及用于实现所述声波测井探测中坏道的识别装置的实施例进行实施，其内容被合并于此，重复之处不再赘述。
图9为本申请实施例的电子设备9600的系统构成的示意框图。如图9所示，该电子设备9600可以包括中央处理器9100和存储器9140；存储器9140耦合到中央处理器9100。值得注意的是，该图9是示例性的；还可以使用其他类型的结构，来补充或代替该结构，以实现电信功能或其他功能。
一实施例中，声波测井探测中坏道的识别功能可以被集成到中央处理器9100中。其中，中央处理器9100可以被配置为进行如下控制：获取目的层段的测井数据；对所述测井数据进行解析得到直达波的到时曲线；以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道。
从上述描述可知，本申请的实施例提供的电子设备，通过获取目的层段的测井数据；对所述测井数据进行解析得到直达波的到时曲线；以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道，能够提高测井数据中坏道的识别准确率和快速性，进而提高测井数据的反射波信噪比和反射体成像精度。
在另一个实施方式中，声波测井探测中坏道的识别装置可以与中央处理器9100分开配置，例如可以将声波测井探测中坏道的识别配置为与中央处理器9100连接的芯片，通过中央处理器的控制来实现声波测井探测中坏道的识别功能。
如图9所示，该电子设备9600还可以包括：通信模块9110、输入单元9120、音频处理器9130、显示器9160、电源9170。值得注意的是，电子设备9600也并不是必须要包括图9中所示的所有部件；此外，电子设备9600还可以包括图9中没有示出的部件，可以参考现有技术。
如图9所示，中央处理器9100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器9100接收输入并控制电子设备9600的各个部件的操作。
其中，存储器9140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器9100可执行该存储器9140存储的该程序，以实现信息存储或处理等。
输入单元9120向中央处理器9100提供输入。该输入单元9120例如为按键或触摸输入装置。电源9170用于向电子设备9600提供电力。显示器9160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。
该存储器9140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器9140还可以是某种其它类型的装置。存储器9140包括缓冲存储器9141(有时被称为缓冲器)。存储器9140可以包括应用/功能存储部9142，该应用/功能存储部9142用于存储应用程序和功能程序或用于通过中央处理器9100执行电子设备9600的操作的流程。
存储器9140还可以包括数据存储部9143，该数据存储部9143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器9140的驱动程序存储部9144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。
通信模块9110即为经由天线9111发送和接收信号的发送机/接收机9110。通信模块(发送机/接收机)9110耦合到中央处理器9100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。
基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块9110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)9110还经由音频处理器9130耦合到扬声器9131和麦克风9132，以经由扬声器9131提供音频输出，并接收来自麦克风9132的音频输入，从而实现通常的电信功能。音频处理器9130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器9130还耦合到中央处理器9100，从而使得可以通过麦克风9132能够在本机上录音，且使得可以通过扬声器9131来播放本机上存储的声音。
本发明的实施例还提供能够实现上述实施例中的声波测井探测中坏道的识别方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的声波测井探测中坏道的识别方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：获取目的层段的测井数据；对所述测井数据进行解析得到直达波的到时曲线；以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道。
从上述描述可知，本发明实施例提供的计算机可读存储介质，通过获取目的层段的测井数据；对所述测井数据进行解析得到直达波的到时曲线；以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道，能够提高测井数据中坏道的识别准确率和快速性，进而提高测井数据的反射波信噪比和反射体成像精度。
虽然本发明提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
本领域技术人员应明白，本说明书的实施例可提供为方法、装置(系统)或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、装置(系统)和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。术语“上”、“下”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；可以是机械连接，也可以是电连接；可以是直接相连，也可以通过中间媒介间接相连，可以是两个元件内部的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。本发明并不局限于任何单一的方面，也不局限于任何单一的实施例，也不局限于这些方面和/或实施例的任意组合和/或置换。而且，可以单独使用本发明的每个方面和/或实施例或者与一个或更多其他方面和/或其实施例结合使用。
最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围，其均应涵盖在本发明的权利要求和说明书的范围当中。
Claims (18)
1.一种声波测井探测中坏道的识别方法，其特征在于，包括：
获取目的层段的测井数据；
对所述测井数据进行解析得到直达波的到时曲线；
以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；
根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道。
2.根据权利要求1所述的声波测井探测中坏道的识别方法，其特征在于，在所述获取目的层段的测井数据并对该测井数据进行解析得到直达波的到时曲线之后，还包括：
对所述测井数据进行解析得到直达波的震荡结束时间；
以所述震荡结束时间至整道波形终止时间为第二时间段，确定所述第二时间段内目的层段中各个深度点的波形幅度值；
根据所述第二时间段内目的层段中各个深度点的波形幅度值和预设的第二临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第二临界值，则确定该深度点的波形为坏道。
3.根据权利要求1或2所述的声波测井探测中坏道的识别方法，其特征在于，还包括：
根据所述坏道相邻的非坏道对所述坏道进行校正。
4.根据权利要求3所述的声波测井探测中坏道的识别方法，其特征在于，所述根据所述坏道相邻的非坏道对所述坏道进行校正，包括：
将所述坏道中的波形替换为相邻的非坏道中的波形。
5.根据权利要求1所述的声波测井探测中坏道的识别方法，其特征在于，在所述获取目的层段的测井数据之前，还包括：
通过阵列声波测井的方式获取目标区块的测井数据；
其中，所述目标区块包括目的层段。
6.根据权利要求1所述的声波测井探测中坏道的识别方法，其特征在于，所述对所述测井数据进行解析得到直达波的到时曲线，包括：
对所述测井数据进行转换处理得到声波的时差曲线；
根据所述时差曲线确定到时曲线。
7.根据权利要求6所述的声波测井探测中坏道的识别方法，其特征在于，所述对所述测井数据进行转换处理得到声波的时差曲线，包括：
采用慢度-时间相关的方式对所述测井数据进行转换处理得到声波的时差曲线。
8.根据权利要求1所述的声波测井探测中坏道的识别方法，其特征在于，所述第一临界值是根据波形幅度曲线平均值确定的。
9.一种声波测井探测中坏道的识别装置，其特征在于，包括：
获取单元，用于获取目的层段的测井数据；
第一解析单元，用于对所述测井数据进行解析得到直达波的到时曲线；
第一处理单元，用于以测井数据中的零时刻至所述到时曲线中的最小值时刻为第一时间段，确定所述第一时间段内目的层段中各个深度点的波形幅度值；
第一识别单元，用于根据所述第一时间段内目的层段中各个深度点的波形幅度值和预设的第一临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第一临界值，则确定该深度点的波形为坏道。
10.根据权利要求9所述的声波测井探测中坏道的识别装置，其特征在于，还包括：
第二解析单元，用于对所述测井数据进行解析得到直达波的震荡结束时间；
第二处理单元，用于以所述震荡结束时间至整道波形终止时间为第二时间段，确定所述第二时间段内目的层段中各个深度点的波形幅度值；
第二识别单元，用于根据所述第二时间段内目的层段中各个深度点的波形幅度值和预设的第二临界值对各个深度点的波形进行坏道识别；其中，若所述深度点的波形幅度值大于第二临界值，则确定该深度点的波形为坏道。
11.根据权利要求9或10所述的声波测井探测中坏道的识别装置，其特征在于，还包括：
校正单元，用于根据所述坏道相邻的非坏道对所述坏道进行校正。
12.根据权利要求11所述的声波测井探测中坏道的识别装置，其特征在于，所述校正单元，包括：
校正子单元，用于将所述坏道中的波形替换为相邻的非坏道中的波形。
13.根据权利要求9所述的声波测井探测中坏道的识别装置，其特征在于，还包括：
采集单元，用于通过阵列声波测井的方式获取目标区块的测井数据；
其中，所述目标区块包括目的层段。
14.根据权利要求9所述的声波测井探测中坏道的识别装置，其特征在于，所述第一解析单元，包括：
第一处理子单元，用于对所述测井数据进行转换处理得到声波的时差曲线；
第二处理子单元，用于根据所述时差曲线确定到时曲线。
15.根据权利要求14所述的声波测井探测中坏道的识别装置，其特征在于，所述第一处理子单元，包括：
转换模块，用于采用慢度-时间相关的方式对所述测井数据进行转换处理得到声波的时差曲线。
16.根据权利要求9所述的声波测井探测中坏道的识别装置，其特征在于，所述第一临界值是根据波形幅度曲线平均值确定的。
17.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至8任一项所述的声波测井探测中坏道的识别方法的步骤。
18.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至8任一项所述的声波测井探测中坏道的识别方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911117129.XA CN112901146B (zh) | 2019-11-15 | 2019-11-15 | 声波测井探测中坏道的识别方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911117129.XA CN112901146B (zh) | 2019-11-15 | 2019-11-15 | 声波测井探测中坏道的识别方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112901146A true CN112901146A (zh) | 2021-06-04 |
| CN112901146B CN112901146B (zh) | 2023-07-25 |
Family
ID=76103037
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911117129.XA Active CN112901146B (zh) | 2019-11-15 | 2019-11-15 | 声波测井探测中坏道的识别方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112901146B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113550747A (zh) * | 2021-08-31 | 2021-10-26 | 中国地质大学（北京） | 一种浊沸石胶结砂砾岩储层测井识别方法及装置 |
| CN115234225A (zh) * | 2022-07-29 | 2022-10-25 | 中国石油大学(华东) | 一种声波远探测数据质量检测方法 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103388473A (zh) * | 2012-05-10 | 2013-11-13 | 中国石油集团长城钻探工程有限公司 | 在声波测井中消除直达波干扰的方法、系统及声波测井仪 |
| CN103424769A (zh) * | 2013-04-26 | 2013-12-04 | 上海申丰地质新技术应用研究所有限公司 | 一种采空区联合多波地震勘探方法 |
| US20170108607A1 (en) * | 2014-07-02 | 2017-04-20 | Halliburton Energy Services Inc. | Peak Analysis of Ultrasonic Waveforms for Cement Bond Logging |
| US20190018161A1 (en) * | 2016-04-01 | 2019-01-17 | Halliburton Energy Services, Inc. | High precision acoustic logging processing for compressional and shear slowness |
-
2019
- 2019-11-15 CN CN201911117129.XA patent/CN112901146B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103388473A (zh) * | 2012-05-10 | 2013-11-13 | 中国石油集团长城钻探工程有限公司 | 在声波测井中消除直达波干扰的方法、系统及声波测井仪 |
| CN103424769A (zh) * | 2013-04-26 | 2013-12-04 | 上海申丰地质新技术应用研究所有限公司 | 一种采空区联合多波地震勘探方法 |
| US20170108607A1 (en) * | 2014-07-02 | 2017-04-20 | Halliburton Energy Services Inc. | Peak Analysis of Ultrasonic Waveforms for Cement Bond Logging |
| US20190018161A1 (en) * | 2016-04-01 | 2019-01-17 | Halliburton Energy Services, Inc. | High precision acoustic logging processing for compressional and shear slowness |
Non-Patent Citations (1)
| Title |
|---|
| 刘娇;朱强;李其斌;王颖;: "基于时频分析的声波时差检测方法研究", 微计算机信息, no. 10, pages 67 - 69 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113550747A (zh) * | 2021-08-31 | 2021-10-26 | 中国地质大学（北京） | 一种浊沸石胶结砂砾岩储层测井识别方法及装置 |
| CN115234225A (zh) * | 2022-07-29 | 2022-10-25 | 中国石油大学(华东) | 一种声波远探测数据质量检测方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112901146B (zh) | 2023-07-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5792311B2 (ja) | ピッチラグを推定すること | |
| CN112901146B (zh) | 声波测井探测中坏道的识别方法及装置 | |
| CN108375785B (zh) | 裂缝带位置校正方法及装置 | |
| CN107605470A (zh) | 一种纵横波径向速度变化成像方法 | |
| KR101483513B1 (ko) | 음원위치추적장치 및 음원위치추적방법 | |
| CN112487613A (zh) | 一种地层波走时的确定方法和装置 | |
| CN111596379B (zh) | 一种地震观测系统异常识别方法和装置 | |
| CN109633744B (zh) | 地震子波的提取方法、装置、设备及存储介质 | |
| KR20180049182A (ko) | 누락 데이터에 대한 사인곡선 보간 | |
| WO2023123971A1 (zh) | 基于vsp的深度域地震剖面层位标定方法及装置 | |
| WO2015124047A1 (zh) | 一种解析波形的方法和装置 | |
| CN112255679B (zh) | 地震资料绕射深度偏移处理方法及装置 | |
| CN112130209B (zh) | 岩溶储层预测方法及装置 | |
| CN112750458B (zh) | 一种触屏声音检测方法和装置 | |
| KR101703104B1 (ko) | 퇴적물 시료의 음파전달속도 및 음파감쇠를 측정하는 음파측정 방법 및 시스템 | |
| CN109655269B (zh) | 一种马达震动检测方法及装置 | |
| CN113672840B (zh) | 斯通利波衰减提取方法及装置 | |
| CN112817042A (zh) | 沉积储层中油气识别方法及装置 | |
| CN112253100B (zh) | 一种固井质量的确定方法和装置 | |
| CN115932969B (zh) | 波场分离方法及装置 | |
| CN118091749B (zh) | 基于动态时间规整的微地震事件初至拾取方法及系统 | |
| CN114721053B (zh) | 一种井中微地震事件的初至拾取方法、装置及存储介质 | |
| CN110632650B (zh) | 全波形反演数据的分时窗提取方法及装置 | |
| CN109143330A (zh) | 一种优化可控震源获取单炮质量的方法及系统 | |
| CN111366974B (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 |