CN107728207B - 地震勘探单炮记录声波识别方法和单炮质量检验方法 - Google Patents

地震勘探单炮记录声波识别方法和单炮质量检验方法 Download PDF

Info

Publication number
CN107728207B
CN107728207B CN201710910574.6A CN201710910574A CN107728207B CN 107728207 B CN107728207 B CN 107728207B CN 201710910574 A CN201710910574 A CN 201710910574A CN 107728207 B CN107728207 B CN 107728207B
Authority
CN
China
Prior art keywords
sound wave
road
recognition methods
window
formula
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
Application number
CN201710910574.6A
Other languages
English (en)
Other versions
CN107728207A (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.)
China National Petroleum Corp
BGP Inc
Original Assignee
BGP Inc
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 BGP Inc filed Critical BGP Inc
Priority to CN201710910574.6A priority Critical patent/CN107728207B/zh
Publication of CN107728207A publication Critical patent/CN107728207A/zh
Application granted granted Critical
Publication of CN107728207B publication Critical patent/CN107728207B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种地震勘探单炮记录声波识别方法和单炮质量检验方法。所述声波识别方法包括步骤：确定声波出现的范围；在时窗内，从起始时间沿时间方向依次取三个相邻小时窗，将其作为整体进行滑动，并分别计算三个小时窗内的能量和中间小时窗的频率，筛选出振幅异常的时间位置；在范围A内，将m道作为拟合单元，判断每个单元内振幅异常的时间位置在横向上是否具有一定的连续性，将具有横向连续特征的道进行线性拟合，筛选出拟合后得到直线斜率k的绝对值大于0.6的拟合单元，则该拟合单元内的道为声波道。所述单炮质量检验方法包括采用如上所述的波识别方法来进行单炮质量检验。本发明的声波识别方法能有效区分单、双支声波，效率高。本发明的单炮质量检验方法能够有助于提高采集到的资料的质量。

Description

地震勘探单炮记录声波识别方法和单炮质量检验方法
技术领域
本发明涉及石油地震勘探采集领域，特别地，涉及一种地震勘探单炮记录声波识别方法和单炮质量检验方法。
背景技术
在地震资料采集中，声波为干扰波，声波的强弱直接影响着单炮地震资料品质。目前在地震资料采集阶段还没有对单炮地震资料中的声波进行识别的技术，进而评价单炮地震资料的采集质量。而在地震资料处理阶段，也只是对地震资料进行声波压制处理。对声波的处理的有很多方法，比较常规的方法有分频去噪和线性压噪，也有利用声波具有非常稳定的传播速度的特点，采用非线性变换手段压制声波干扰。这些方法在处理资料的同时，对有效波的影响是不可避免的。对单炮地震记录中声波干扰的识别一直也是一个难点，传统声波处理技术效率不高，只能初略计算出声波干扰强弱，而且对于没有声波干扰的单炮记录，也会错误地计算得到声波能量。
发明内容
针对现有技术中存在的不足，本发明的目的之一在于解决上述现有技术中存在的一个或多个问题。例如，本发明的目的之一在于提供一种声波识别方法，来快速有效地识别地震单炮记录中的声波干扰。
为了实现上述目的，本发明一方面提供了一种地震勘探单炮记录声波识别方法，其特征在于，所述方法包括步骤：根据式1确定声波出现的起始时间t0，并将t0沿时间方向延伸1.5s～3s，得到包含了干扰道的范围A内的每道的第一时窗，其中，式1为t0＝d/v，d为偏移距，v为声波速度；在每个所述的第一时窗内，从起始时间t0沿时间方向依次取150ms～300ms大小的三个相邻的小时窗TW1、TW2和TW3，将三个小时窗作为一个整体沿时间方向进行滑动，直至TW3的底部到达第一时窗的底部，在滑动过程中分别计算三个小时窗内的能量和TW2的频率，筛选出满足式2的TW2能量最大的时间位置，即振幅异常的时间位置，式2为：TW2E/TW1E>1&&TW2E/TW3E>1.1&&TW2F>40，其中，TW1E、TW2E和TW3E分别为三个小时窗TW1、TW2和TW3内的能量，TW2F为TW2的频率；在含有N道的范围A内，将起始于范围A第一边界的m道作为起始拟合单元，起始拟合单元向范围A的第二边界横向移动，每次移动一道，得一个拟合单元，最终得到N-m+1个拟合单元，判断每个单元内振幅异常的时间位置在横向上是否具有一定的连续性，将具有横向连续特征的道进行线性拟合，筛选出拟合后得到直线斜率k的绝对值大于0.6的拟合单元，则该拟合单元内的道为声波道；其中，第一边界和第二边界沿时间方向延伸。
本发明另一方面提供了一种地震勘探单炮质量检验方法，所述方法包括采用如上所述的地震勘探单炮记录声波识别方法来进行单炮质量检验。
与现有技术相比，本发明的地震勘探单炮记录声波识别方法能有效区分单、双支声波，能够计算声波能量强弱，效率高，识别准确率目前可以达到90％以上。本发明的地震勘探单炮质量检验方法有助于提高采集到的资料的质量。
附图说明
通过下面结合附图进行的描述，本发明的上述和其他目的和特点将会变得更加清楚，其中：
图1示出了地震勘探单炮记录声波识别方法的流程示意图。
图2示出了包含了干扰道的范围A内每道的第一时窗的示意图。
图3示出了本发明识别出的双支声波。
图4示出了本发明识别出的单支声波。
具体实施方式
在下文中，将结合附图和示例性实施例详细地描述本发明的地震勘探单炮记录声波识别方法和地震勘探单炮质量检验方法。本发明出现的第一和第二不表示先后顺序，仅用于相互区别。
本发明一方面提供了一种地震勘探单炮记录声波识别方法，所述方法的创新点包括：1)充分利用声波在记录剖面上呈现的振幅异常、频率高、振幅异常处横向连续且近偏排列呈线性关系三个特点，根据每个特点逐一筛选出声波，快速高效；2)发明不仅识别出声波，而且计算出了声波影响道数、能量以及声波类型。
本发明声波识别方法的技术方案包括：充分利用声波在地震单炮记录剖面的三个明显特征，即在声波出现的时间范围内，必然出现振幅异常，频率增高，并且振幅异常的相邻道横向连续且呈线性关系。图1示出了地震勘探单炮记录声波识别方法的流程示意图，从图中可以看出本发明的声波识别方法的具体步骤。本发明实施的技术方案是首先计算得到每一地震道声波出现的时间范围，然后查找单炮记录中每一道中振幅异常、频率高的位置。然后将每一道(又称地震道)的异常位置时间使用最小二乘法线性拟合，拟合得到异常点连线的斜率，判断斜率是否满足声波干扰的斜率特点。对满足声波斜率的异常点，则认为是声波干扰引起。
根据本发明示例性实施例的地震勘探单炮记录声波识别方法可包括以下步骤：
根据式1确定声波出现的起始时间t0，并将t0沿时间方向延伸1.5s～3s，优选的，为2s，得到包含了干扰道(又称声波道或声波干扰道)的范围A内的每道的第一时窗，即声波出现的范围，声波出现的时间范围包括纵向时间范围和横向地震道范围；其中，式1为t0＝d/v，d为偏移距，v为声波速度；其中，声波速度v为360m/s，实际地震资料采集过程中，声波速度可能受不同环境的影响有所变化，为确保第一时窗能够包含声波出现的位置，将v取为360m/s，以此计算出的时间作为时窗的起始位置。该步骤中，由于声波在每道中的持续时间一般不会超过500ms，将t0沿时间方向延伸1.5s～3s(优选为2.0s)作为第一时窗的终止位置，是为了确保第一时窗包含声波，但延伸范围不能太大，太大对于查找声波导致的异常点没有意义。图2示出了包含干扰道的范围A内每道的第一时窗的示意图，图中的A(即范围A)是呈倒V型的区域，范围A里面的道包含了干扰道。
在每个所述的第一时窗内，从起始时间t0沿时间方向依次取150ms～300ms大小的三个相邻的小时窗TW1、TW2和TW3，将三个小时窗作为一个整体沿时间方向进行滑动，如图2所示，直至TW3的底部到达第一时窗的底部，在滑动过程中分别计算三个小时窗内的能量和TW2的频率，筛选出满足式2的TW2能量(即小时窗TW2内的能量)最大的时间位置，即振幅异常的时间位置(又称作异常振幅位置、振幅异常位置，振幅异常点等)，式2为：TW2E/TW1E>1&&TW2E/TW3E>1.1&&TW2F>40，即同时满足TW2E/TW1E>1、TW2E/TW3E>1.1和TW2F>40Hz三个式子，其中，TW1E、TW2E和TW3E分别为三个小时窗TW1、TW2和TW3的能量，TW2F为小时窗TW2的频率。声波干扰引起振幅异常，就是声波的地震道样点的幅值明显比声波到来前后的样点幅值都大，通过三个小时窗内样点能量的比较，就是为了选出声波异常时间位置。在该步骤中，地震道中声波持续时间一般为200ms左右，从t0沿时间方向依次取150ms～300ms大小的三个小时窗，这样找出的振幅异常位置可能有差异，但不会太大；优选的，可从起始时间t0沿时间方向依次取200ms大小的三个相邻的小时窗TW1、TW2和TW3，选择200ms作为小时窗来找声波引起的振幅异常比较合适。
在含有N道的范围A内，将起始于范围A第一边界的m道作为起始拟合单元，起始拟合单元向范围A的第二边界横向移动，每次移动一道，得一个拟合单元，最终得到N-m+1个拟合单元，判断每个单元内振幅异常的时间位置在横向上是否具有一定的连续性，将具有横向连续特征的道进行线性拟合，筛选出拟合斜率k(即拟合后得到直线斜率k)的绝对值大于0.6的拟合单元，则该拟合单元内的道为声波道；其中，第一边界和第二边界沿时间方向延伸，如图2所示，第一边界和第二边界为包含干扰道的范围A左右两边沿时间方向延伸的边界。本步骤可以剔除不是由声波引起的振幅异常点。其中，N、m都为正整数，N大于m，m不小于5。当每一地震道中第一时窗范围内三个时窗滑动完毕后，若不存在振幅异常点，异常振幅位置时间t＝0；异常振幅横向连续，就是说连续的地震道中异常振幅位置时间t都不为0，就是说这些地震道都出现了振幅异常点。声波是一种线性干扰波，当连续的异常振幅值拟合后斜率大于0.6后，就是声波引起。斜率大于0.6的单元内所有道都为干扰道。
在本实施例中，三个小时窗内的能量可根据式3计算，式3为：
其中，σ为能量，a为样点幅值，n为时窗样点数。
地震单炮记录的是离散化的波形数据，记录参数一般是采样率1ms，记录长度6s，每一道的波形就离散化为6001的采样点来记录，每个采样点记录波形的振幅值。样点幅值就是采样点的振幅值。假如地震单炮记录中有100道，那么就记录有100*6001个样点值。
在本实施例中，TW2的频率计算方法可包括快速傅里叶变换(简称FFT，即FastFourier Transform)。
在本实施例中，不同单炮记录，N的大小取决于地震资料采集参数；m可为8～12，优选的，为10，即在出现声波的道范围内，将连续的10道作为1个拟合单元。当m为10时，判断每个单元内振幅异常的时间位置在横向上是否具有一定的连续性，将具有横向连续特征的道进行线性拟合的步骤包括：若单元内10道中至少有6道的振幅的异常时间位置满足递增或递减，则用最小二乘法进行线性拟合得到直线斜率k。单元内10道中需要至少6道来满足条件，这是经验性结果，当10道中有6道都存在振幅异常点，那么可以认为异常振幅具有横向连续的特征，可能是由声波干扰引起，需要通过拟合计算斜率进一步确认。如果少于6道，可能，那么横向连续的条件就不太吻合了，不需要再进行拟合计算斜率。
在本实施例中，直线斜率k的计算公式可为：
式中，xi为道号，为道号均值，ti为振幅异常的时间位置，为时间均值，p为拟合道数(即拟合单元里具有横向连续特征的道数)。
在本实施例中，本发明还可包括步骤：计算声波的总道数和声波能量的大小。
在本实施例中，本发明还可包括步骤：在确定声波出现的起始时间t0之前，计算偏移距d。
在本实施例中，本发明还可包括判断声波类型的步骤；其中，可人眼直接看地震记录剖面就可看出声波类型。
根据本发明的另一个示例性实施例，本发明的地震勘探单炮记录声波识别方法可包括以下步骤：
第一步：根据偏移距和声波速度确定声波出现范围。
利用每一道的偏移距，在声波速度340m/s的基础上，给定声波最大速度360m/s，用公式t0＝d/v(式中t0为时间，d为偏移距，v为声波速度)得到声波出现的起始时间，然后向下延伸2s，获得声波出现的时窗范围，如图2所示出的时窗范围。
第二步：滑动小时窗，确定出满足声波频率的振幅异常位置。
对图2中虚线范围中每一道，从时窗顶部开始，依次取200ms大小的TW1～TW3三个时窗。分别计算每个小时窗内的频率和能量。三个小时窗能量分别记为TW1E、TW2E、TW3E，能量计算公式为：
式中σ为能量，a为样点幅值，n为时窗样点数。记中间小时窗频率为TW2F，频率计算公式为快速傅里叶变换FFT。三个时窗同时以一个样点的间距向下滑动，依次计算小时窗之间的能量比值大小以及TW2F，直至滑动到大时窗底部，在这滑动过程中筛选出满足条件：TW2E/TW1E>1&&TW2E/TW3E>1.1&&TW2F>40的TW2出现的时间位置，从中选择TW2E最大的位置。
第三步：线性拟合异常振幅位置，筛选出横向具有一定连续性且斜率满足条件的道，即为声波道。
在声波出现的道范围内，从左至右，以10道作为一个拟合单位。若这10道中至少有6道的振幅异常点满足递增或递减，则用最小二乘法进行线性拟合得到直线斜率k。这样做的目的是为剔除不是由声波引起的振幅异常点。计算斜率公式为：
式中p为拟合道数，xi为道号，为道号均值，ti为异常振幅位置时间，为时间均值。
每滑动一道，拟合一次，判断拟合直线斜率是否大于0.6，满足条件则认为这些道为声波干扰道记录下满足条件的道号，最后计算影响的总道数和声波能量大小。
图3、4为采用本技术识别出来的双支声波与单支声波，从图中的效果可以看出，此算法能有效区分单、双支声波。表1示出了图3和图4识别出声波的结果。图3示出的单炮地震记录中，图中两个方框内为声波，表1结果显示声波影响地震道52道，声波干扰强弱40.59，类型为双支声波；图4示出的单炮地震记录中，图中方框内为声波，方框内未识别结果，表1结果显示声波影响地震道17道，声波干扰强弱9.64，类型为单支声波。从实际应用的效果可以看出，此算法能有效区分单、双支声波，计算声波能量强弱，运行效率高，识别准确率目前可以达到90％。
表1
线号 点号 死道数 声波影响地震道数 声波干扰强弱 声波类型
图3 1 1414.5 8 52道 40.59 双支声波
图4 1 820.5 2 17道 9.64 单支声波
本发明另一方面提供了一种地震勘探单炮质量检验方法，所述方法包括采用如上所述的地震勘探单炮记录声波识别方法来进行单炮质量检验。将本发明的地震勘探单炮质量监控方法应用到地震资料采集的质量检验(例如质量监控及资料评价)，能够有利于提高采集到的资料的质量。
本发明的地震勘探单炮记录声波识别方法的应用不局限于地震勘探单炮质量检验中，也可将本发明的地震勘探单炮记录声波识别方法集成在单炮质量监控软件中，为全面准确评价采集地震资料质量奠定坚实基础，使单炮监控软件的质控效果更好。
综上所述，本发明的地震勘探单炮记录声波识别方法的优点包括：能够快速有效地识别地震单炮记录中的声波干扰，识别效果非常好，能有效区分单、双支声波，能够计算声波能量强弱，运行效率高，识别准确率目前可以达到90％以上。本发明的地震勘探单炮质量监控方法能够指导采集施工，提高采集到的资料的质量。
尽管上面已经通过结合示例性实施例描述了本发明，但是本领域技术人员应该清楚，在不脱离权利要求所限定的精神和范围的情况下，可对本发明的示例性实施例进行各种修改和改变。

Claims (7)

1.一种地震勘探单炮记录声波识别方法，其特征在于，所述声波识别方法包括步骤：
根据式1确定声波出现的起始时间t0，并将t0沿时间方向延伸1.5s～3s，得到包含干扰道的范围A内每道的第一时窗，其中，式1为t0＝d/v，d为偏移距，v为声波速度；
在每个所述的第一时窗内，从起始时间t0沿时间方向依次取150ms～300ms大小的三个相邻的小时窗TW1、TW2和TW3，将三个小时窗作为一个整体沿时间方向进行滑动，直至TW3的底部到达第一时窗的底部，在滑动过程中分别计算三个小时窗内的能量和TW2的频率，筛选出满足式2的TW2，从中选择TW2E最大的时间位置，即振幅异常的时间位置，式2为：TW2E/TW1E>1且TW2E/TW3E>1.1且TW2F>40，其中，TW1E、TW2E和TW3E分别为三个小时窗TW1、TW2和TW3内的能量，TW2F为TW2的频率；
在范围A内，将起始于范围A第一边界的m道作为起始拟合单元，起始拟合单元向范围A的第二边界横向移动，每次移动一道，得一个拟合单元，最终得到N-m+1个拟合单元，判断每个单元内振幅异常的时间位置在横向上是否具有一定的连续性，将具有横向连续特征的道进行线性拟合，筛选出拟合后得到直线斜率k的绝对值大于0.6的拟合单元，则该拟合单元内的道为声波道；其中，第一边界和第二边界沿时间方向延伸，所述范围A含有N道，
其中，所述三个小时窗内的能量根据式3计算，式3为：
其中，σ为能量，a为样点幅值，n为时窗样点数。
2.根据权利要求1所述的地震勘探单炮记录声波识别方法，其特征在于，所述TW2的频率计算方法包括快速傅里叶变换。
3.根据权利要求1所述的地震勘探单炮记录声波识别方法，其特征在于，所述直线斜率k的计算公式为：
式中，xi为道号，为道号均值，ti为振幅异常的时间位置，为时间均值，p为拟合道数。
4.根据权利要求1所述的地震勘探单炮记录声波识别方法，其特征在于，所述m为10。
5.根据权利要求4所述的地震勘探单炮记录声波识别方法，其特征在于，所述判断每个单元内振幅异常的时间位置在横向上是否具有一定的连续性，将具有横向连续特征的道进行线性拟合的步骤包括：
若单元内10道中至少有6道的振幅的异常时间位置满足递增或递减，则用最小二乘法进行线性拟合得到直线斜率k。
6.根据权利要求5所述的地震勘探单炮记录声波识别方法，其特征在于，所述直线斜率k的计算公式为：
式中，xi为道号，为道号均值，ti为振幅异常的时间位置，为时间均值，p为拟合道数。
7.一种地震勘探单炮质量检验方法，其特征在于，所述方法包括采用如权利要求1～6中任意一项所述的地震勘探单炮记录声波识别方法来进行单炮质量检验。
CN201710910574.6A 2017-09-29 2017-09-29 地震勘探单炮记录声波识别方法和单炮质量检验方法 Active CN107728207B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710910574.6A CN107728207B (zh) 2017-09-29 2017-09-29 地震勘探单炮记录声波识别方法和单炮质量检验方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710910574.6A CN107728207B (zh) 2017-09-29 2017-09-29 地震勘探单炮记录声波识别方法和单炮质量检验方法

Publications (2)

Publication Number Publication Date
CN107728207A CN107728207A (zh) 2018-02-23
CN107728207B true CN107728207B (zh) 2019-03-26

Family

ID=61209213

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710910574.6A Active CN107728207B (zh) 2017-09-29 2017-09-29 地震勘探单炮记录声波识别方法和单炮质量检验方法

Country Status (1)

Country Link
CN (1) CN107728207B (zh)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011005353A1 (en) * 2009-07-06 2011-01-13 Exxonmobil Upstream Research Company Method for seismic interpretation using seismic texture attributes
CN102004264A (zh) * 2010-10-18 2011-04-06 中国石油化工股份有限公司 一种地震采集资料质量定量分析与评价方法
CN103742131A (zh) * 2014-01-20 2014-04-23 电子科技大学 随钻声波井下信号采集与处理***的时差实时提取方法
EP2781936A2 (en) * 2013-03-22 2014-09-24 CGG Services SA System and method for interpolating seismic data
CN104345341A (zh) * 2013-08-08 2015-02-11 中国石油化工股份有限公司 一种基于区域约束的分频段能量地震面波处理方法
CN106094024A (zh) * 2016-06-21 2016-11-09 中国石油天然气集团公司 一种可控震源滑动扫描数据采集处理方法及装置
WO2016205608A1 (en) * 2015-06-17 2016-12-22 Conocophillips Company Seismic azimuthal gradient estimation
CN106646599A (zh) * 2016-12-28 2017-05-10 中国石油化工股份有限公司 针对地表响应因素产生谐波的自动识别与衰减方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011005353A1 (en) * 2009-07-06 2011-01-13 Exxonmobil Upstream Research Company Method for seismic interpretation using seismic texture attributes
CN102004264A (zh) * 2010-10-18 2011-04-06 中国石油化工股份有限公司 一种地震采集资料质量定量分析与评价方法
EP2781936A2 (en) * 2013-03-22 2014-09-24 CGG Services SA System and method for interpolating seismic data
CN104345341A (zh) * 2013-08-08 2015-02-11 中国石油化工股份有限公司 一种基于区域约束的分频段能量地震面波处理方法
CN103742131A (zh) * 2014-01-20 2014-04-23 电子科技大学 随钻声波井下信号采集与处理***的时差实时提取方法
WO2016205608A1 (en) * 2015-06-17 2016-12-22 Conocophillips Company Seismic azimuthal gradient estimation
CN106094024A (zh) * 2016-06-21 2016-11-09 中国石油天然气集团公司 一种可控震源滑动扫描数据采集处理方法及装置
CN106646599A (zh) * 2016-12-28 2017-05-10 中国石油化工股份有限公司 针对地表响应因素产生谐波的自动识别与衰减方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"基于时窗种子提取的海洋地震资料侧面噪声压制方法";王征 等;《石油物探》;20150531;第54卷(第3期);第274-281页

Also Published As

Publication number Publication date
CN107728207A (zh) 2018-02-23

Similar Documents

Publication Publication Date Title
CN104158611B (zh) 基于频谱分析的无线信号干扰检测***及方法
CN103048593B (zh) 一种气体绝缘开关设备绝缘缺陷种类的识别方法
CN104198898B (zh) 基于脉冲序列分析的局部放电发展过程诊断方法
Chan et al. Time-frequency sparsity map on automatic partial discharge sources separation for power transformer condition assessment
CN109034127B (zh) 一种频谱异常检测方法、装置和电子设备
CN108814590B (zh) 一种心电qrs波群的检测方法及其心电分析方法
CN102792173B (zh) 用于检测部分放电的仪器和方法
CN108693448B (zh) 一种应用于电力设备局部放电模式识别***
CN103424183B (zh) 机械振动信号检测异常干扰消除方法
CN106385288B (zh) 一种雷达脉冲跨信道融合检测装置及其检测方法
CN105445699B (zh) 一种非视距误差消除的测距方法及***
CN108344922B (zh) 一种基于相类区分和s变换的输电线路直击雷故障识别方法
CN106845334A (zh) 一种基于数学形态学的新型噪声提取方法
CN106382981A (zh) 一种单站次声波信号识别提取方法
CN111089856A (zh) 一种拉曼光谱弱信号提取的后处理方法
CN103308829B (zh) 一种gis单次局放信号提取与触发时刻调整方法
CN109425894A (zh) 一种地震异常道检测方法及装置
CN103344989A (zh) 可控震源地震记录中脉冲噪声干扰的分析方法
CN106772572B (zh) 一种微地震监测初至的拾取方法
CN109596354B (zh) 基于自适应共振频带识别的带通滤波方法
CN109061591A (zh) 一种基于序贯聚类的时频线谱检测方法
CN111025100A (zh) 变压器特高频局部放电信号模式识别方法与装置
CN107749174B (zh) 抗地铁干扰的停车检测方法、装置、设备及存储介质
CN107728207B (zh) 地震勘探单炮记录声波识别方法和单炮质量检验方法
CN104570118B (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
CB02 Change of applicant information
CB02 Change of applicant information

Address after: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan

Applicant after: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

Address before: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan

Applicant before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20180330

Address after: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei

Applicant after: BGP INC., CHINA NATIONAL PETROLEUM Corp.

Address before: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan

Applicant before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20201112

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Patentee after: BGP Inc., China National Petroleum Corp.

Address before: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei

Patentee before: BGP Inc., China National Petroleum Corp.