CN117457024A - 数据处理方法、装置、计算机可读存储介质及计算机设备 - Google Patents

数据处理方法、装置、计算机可读存储介质及计算机设备 Download PDF

Info

Publication number
CN117457024A
CN117457024A CN202210843496.3A CN202210843496A CN117457024A CN 117457024 A CN117457024 A CN 117457024A CN 202210843496 A CN202210843496 A CN 202210843496A CN 117457024 A CN117457024 A CN 117457024A
Authority
CN
China
Prior art keywords
frequency domain
domain waveforms
waveforms
waveform
target
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.)
Pending
Application number
CN202210843496.3A
Other languages
English (en)
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
China Petroleum Logging Co Ltd
Original Assignee
China National Petroleum Corp
China Petroleum Logging Co Ltd
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 China National Petroleum Corp, China Petroleum Logging Co Ltd filed Critical China National Petroleum Corp
Priority to CN202210843496.3A priority Critical patent/CN117457024A/zh
Publication of CN117457024A publication Critical patent/CN117457024A/zh
Pending legal-status Critical Current

Links

Classifications

    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/48—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
    • G10L25/51—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for comparison or discrimination
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/18—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being spectral information of each sub-band
    • G—PHYSICS
    • G10—MUSICAL INSTRUMENTS; ACOUSTICS
    • G10L—SPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • G10L25/21—Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters the extracted parameters being power information

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computational Linguistics (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)

Abstract

本发明公开了一种数据处理方法、装置、计算机可读存储介质及计算机设备。其中，该方法包括：获取目标井不同深度点对应的多个第一时域波形；分别获取与多个第一时域波形对应的多个第一频域波形；分别获取与多个第一频域波形对应的参照频域波形；分别基于多个第一频域波形和对应的参照频域波形，确定多个第一频域波形中的异常第一频域波形；将异常第一时域波形从多个第一时域波形中清除出去，得到多个第二时域波形；基于多个第二时域波形，确定目标井的目标时域波形。本发明解决了测井过程中因异常数据造成的数据分析结果不准确且重新采集数据所耗费的成本过高的技术问题。

Description

数据处理方法、装置、计算机可读存储介质及计算机设备
技术领域
本发明涉及地球物理测井领域，具体而言，涉及一种数据处理方法、装置、计算机可读存储介质及计算机设备。
背景技术
在相关技术中，在利用声波测井过程中可能会出现数据缺失或者数据错误，而这些异常数据会影响后续对声波数据的分析和处理，而重新采集数据又会增加采集的数据的成本。
因此，在相关技术中，存在测井过程中因异常数据造成的数据分析结果不准确且重新采集数据所耗费的成本过高的技术问题。
针对上述的问题，目前尚未提出有效的解决方案。
发明内容
本发明实施例提供了一种数据处理方法、装置、计算机可读存储介质及计算机设备，以至少解决测井过程中因异常数据造成的数据分析结果不准确且重新采集数据所耗费的成本过高的技术问题。
根据本发明实施例的一个方面，提供了一种数据处理方法，包括：获取目标井不同深度点对应的多个第一时域波形；分别获取与多个第一时域波形对应的多个第一频域波形；分别获取与多个第一频域波形对应的参照频域波形；分别基于多个第一频域波形和对应的参照频域波形，确定多个第一频域波形中的异常第一频域波形；将异常第一时域波形从多个第一时域波形中清除出去，得到多个第二时域波形；基于多个第二时域波形，确定目标井的目标时域波形。
可选地，分别获取与多个第一频域波形对应的参照频域波形，包括：基于多个第一频域波形对应的深度点的顺序，将多个第一频域波形划分为多个频域波形组；分别确定与多个第一频域波形对应的频域波形组；分别确定与多个第一频域波形对应的频域波形组所包括的深度点中，位于中间位置处的目标深度点；分别基于目标深度点对应的第一频域波形，确定与多个第一频域波形对应的参照频域波形。
可选地，分别基于目标深度点对应的第一频域波形，确定与多个第一频域波形对应的参照频域波形，包括：在目标深度点的个数为1的情况下，将与目标深度点对应的第一频域波形确定为与多个第一频域波形对应的参照频域波形；在目标深度点的个数为2的情况下，对与目标深度点对应的两个第一频域波形进行平均处理，得到第一频域平均波形，并将第一频域平均波形确定为与多个第一频域波形对应的参照频域波形。
可选地，分别基于多个第一频域波形和对应的参照频域波形，确定多个第一频域波形中的异常第一频域波形，包括：分别获取多个第一频域波形和对应参照频域波形之间的差异系数；在差异系数大于第一预定阈值的情况下，确定对应的第一频域波形为异常第一频域波形。
可选地，基于多个第二时域波形，确定目标井的目标时域波形，包括：分别获取与多个第二时域波形对应的多个第二频域波形；对多个第二频域波形中信号幅度不大于第二预定阈值的幅度进行置零，得到多个目标频域波形；分别确定与多个目标频域波形对应的多个目标时域波形。
可选地，对第二频域波形中信号幅度不大于第二预定阈值的幅度进行置零，得到目标频域波形，包括：采用更新第二预定阈值的方式，对第二频域波形按照预定次数循环迭代进行置零操作，得到目标频域波形，其中，更新第二预定阈值包括：基于第二频域波形对应的最大幅值和预定参数更新第二预定阈值。
根据本发明实施例的另一方面，还提供了一种数据处理装置，包括：第一获取模块，用于获取目标井不同深度点对应的多个第一时域波形；第二获取模块，用于分别获取与多个第一时域波形对应的多个第一频域波形；第三获取模块，用于分别获取与多个第一频域波形对应的参照频域波形；第一确定模块，用于分别基于多个第一频域波形和对应的参照频域波形，确定多个第一频域波形中的异常第一频域波形；清除模块，用于将异常第一时域波形从多个第一时域波形中清除出去，得到多个第二时域波形；第二确定模块，用于基于多个第二时域波形，确定目标井的目标时域波形。
可选地，第三获取模块包括：划分单元，用于基于多个第一频域波形对应的深度点的顺序，将多个第一频域波形划分为多个频域波形组；第一确定单元，用于分别确定与多个第一频域波形对应的频域波形组；第二确定单元，用于分别确定与多个第一频域波形对应的频域波形组所包括的深度点中，位于中间位置处的目标深度点；第三确定单元，用于分别基于目标深度点对应的第一频域波形，确定与多个第一频域波形对应的参照频域波形。
根据本发明实施例的另一方面，还提供了一种计算机可读存储介质，计算机可读存储介质包括存储的程序，其中，在程序运行时控制计算机可读存储介质所在设备执行上述任意一项的数据处理方法。
根据本发明实施例的另一方面，还提供了一种计算机设备，包括：存储器和处理器，存储器存储有计算机程序；处理器，用于执行存储器中存储的计算机程序，计算机程序运行时使得处理器执行上述任意一项的数据处理方法。
在本发明实施例中，通过将对目标井在不同深度点处采集得到的时域波形转换为频域波形，并将频域波形与参照频域波形进行比对，确定出存在明显异常数据的频域波形数据，将这些存在异常数据的频域波形在时域中对应的波形数据进行清除，达到了清除采集得到的时域波形数据中异常数据的目的，从而实现了从原始的时域波形中确定出异常数据并清除异常数据的技术效果，进而解决了测井过程中因异常数据造成的数据分析结果不准确且重新采集数据所耗费的成本过高技术问题。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1是根据本发明实施例的数据处理方法的流程图；
图2是根据本发明可选实施方式的声波测井坏道自适应识别与重建方法的流程图；
图3是根据本发明可选实施方式的差异系数分布示意图；
图4是根据本发明可选实施方式的坏道识别结果示意图；
图5是根据本发明可选实施方式的坏道深度段置零示意图；
图6是根据本发明可选实施方式的坏道深度段数据与重建结果对比图；
图7是根据本发明实施例的数据处理装置的结构框图。
具体实施方式
为了使本技术领域的人员更好地理解本发明方案，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分的实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都应当属于本发明保护的范围。
需要说明的是，本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象，而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换，以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。此外，术语“包括”和“具有”以及他们的任何变形，意图在于覆盖不排他的包含，例如，包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元，而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
术语说明
声波测井，指利用声波在不同岩石中传播时，速度、幅度及频率的变化等声学特性不相同来研究钻井的地质剖面，判断固井质量的一种测井方法。
差异系数，也称变差系数、离散系数、变异系数，用CV表示。它是一组数据的标准差与其均值的百分比，是测算数据离散程度的相对指标，是一种相对差异量数。
根据本发明实施例，提供了一种数据处理的方法实施例，需要说明的是，在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行，并且，虽然在流程图中示出了逻辑顺序，但是在某些情况下，可以以不同于此处的顺序执行所示出或描述的步骤。
图1是根据本发明实施例的数据处理方法的流程图，如图1所示，该方法包括如下步骤：
步骤S102，获取目标井不同深度点对应的多个第一时域波形；
步骤S104，分别获取与多个第一时域波形对应的多个第一频域波形；
步骤S106，分别获取与多个第一频域波形对应的参照频域波形；
步骤S108，分别基于多个第一频域波形和对应的参照频域波形，确定多个第一频域波形中的异常第一频域波形；
步骤S110，将异常第一时域波形从多个第一时域波形中清除出去，得到多个第二时域波形；
步骤S112，基于多个第二时域波形，确定目标井的目标时域波形。
通过上述步骤，将对目标井在不同深度点处采集得到的时域波形转换为频域波形，并将频域波形与参照频域波形进行比对，确定出存在明显异常数据的频域波形数据，将这些存在异常数据的频域波形在时域中对应的波形数据进行清除，达到了清除采集得到的时域波形数据中异常数据的目的，从而实现了从原始的时域波形中确定出异常数据并清除异常数据的技术效果，进而解决了测井过程中因异常数据造成的数据分析结果不准确且重新采集数据所耗费的成本过高技术问题。
作为一种可选的实施例，分别获取与多个第一频域波形对应的参照频域波形时，可以采用多种方式，例如，可以采用以下方式：基于多个第一频域波形对应的深度点的顺序，将多个第一频域波形划分为多个频域波形组；分别确定与多个第一频域波形对应的频域波形组；分别确定与多个第一频域波形对应的频域波形组所包括的深度点中，位于中间位置处的目标深度点；分别基于目标深度点对应的第一频域波形，确定与多个第一频域波形对应的参照频域波形。在确定与第一频域波形对应的参照频域波形时，考虑到深度点位置距离较大时对应获得的波形数据也会有一定程度上的差异，为了使确定出来的参照频域波形具有足够可靠的参照性，可以采用先按照深度点对频域波形进行分组再确定出每组频域波形对应的参照频域波形的方式，例如，可以将每5个深度点对应的频域波形确定为一个频域波形组，以这种方式，可以在不同深度范围的采用不同的参照频域波形，避免了由于深度点位置差异较大却采用某一固定参照频域波形而造成异常数据判断失误的情况。同时，为了对频域波形组中的每个第一频域波形都确定出一个参照性较好的参照频域波形，还可以从频域波形组所包括的深度点中确定出位于中间位置处的目标深度点，并将该目标深度点对应的频域波形确定为参照频域波形。
作为一种可选的实施例，分别基于目标深度点对应的第一频域波形，确定与多个第一频域波形对应的参照频域波形时，可以采用多种方式，例如，可以采用以下方式：在目标深度点的个数为1的情况下，将与目标深度点对应的第一频域波形确定为与多个第一频域波形对应的参照频域波形；在目标深度点的个数为2的情况下，对与目标深度点对应的两个第一频域波形进行平均处理，得到第一频域平均波形，并将第一频域平均波形确定为与多个第一频域波形对应的参照频域波形。确定目标深度点对应的参照频域波形时，若目标深度点为一个，则可以直接将该深度点对应的频域波形确定为参照频域波形，若目标深度点位两个，则可以将这两个目标深度点对应的频域波形做平均处理，将平均后的频域波形确定为参照频域波形。
作为一种可选的实施例，分别基于多个第一频域波形和对应的参照频域波形，确定多个第一频域波形中的异常第一频域波形时，可以采用多种方式，例如，可以采用以下方式：分别获取多个第一频域波形和对应参照频域波形之间的差异系数；在差异系数大于第一预定阈值的情况下，确定对应的第一频域波形为异常第一频域波形。在确定了与第一频域波形对应的参照频域波形之后，可以通过计算上述两者之间的差异系数的方式来判断第一频域波形与参照频域波形的差异情况，若计算得到的差异系数大于第一预定阈值，则可以认为该第一频域波形与参照频域波形相比差异较大，即该第一频域波形的波形数据是异常的，存在明显的数据缺失或数据错误的情况，进而可以通过上述方法确定出第一频域波形中全部的异常数据。
作为一种可选的实施例，基于多个第二时域波形，确定目标井的目标时域波形时，可以采用多种方式，例如，可以采用以下方式：分别获取与多个第二时域波形对应的多个第二频域波形；对多个第二频域波形中信号幅度不大于第二预定阈值的幅度进行置零，得到多个目标频域波形；分别确定与多个目标频域波形对应的多个目标时域波形。在从第一时域波形中清除了存在异常数据的时域波形之后，在剩余的第二时域波形中仍可能存在受异常数据影响的部分数据，例如，对于存在异常数据的波形，其异常数据可能会降低附近波形的振幅值，使其它波形中呈现出弱振幅，因此，为了消除异常数据带来的影响，可以采用第二预定阈值来对与第二时域波形对应的第二频域波形进行幅值上的判断，将第二频域波形上幅值小于第二预定阈值的部分确定为是受异常数据影响的，并通过将这部分数据置零的方式来消除异常数据在第二时域波形或第二频域波形中的影响。
作为一种可选的实施例，对第二频域波形中信号幅度不大于第二预定阈值的幅度进行置零，得到目标频域波形时，可以采用多种方式，例如，可以采用以下方式：采用更新第二预定阈值的方式，对第二频域波形按照预定次数循环迭代进行置零操作，得到目标频域波形，其中，更新第二预定阈值包括：基于第二频域波形对应的最大幅值和预定参数更新第二预定阈值。为了达到对于目标时域波形的精度要求，可以不断地利用新得到的频域波形对应的幅值和预定参数来更新第二预定阈值，并利用更新后的第二预定阈值对上述新得到的频域波形进行新一轮的操作置零，当第二频域波形经过预定次数的指令操作之后，就可以认为得到的是满足精度要求的目标频域波形，这也就意味着与该目标频域波形对应的就是目标时域波形，即与目标井对应的目标时域波形。
需要说明的是，在上述将第一时域波形转换至频域时可以采用傅里叶变换的方式，而将第一频域波形转换至时域时则可以采用傅里叶逆变换的方式。针对第二时域波形可以利用二维傅里叶变换的方式将其转换至频域，而将目标频域波形转换为目标时域波形时，可以采用二维傅里叶逆变换的方式。
基于上述实施例及可选实施例，本发明提出一种可选实施方式，下面进行说明。
在声波测井过程中，由于测速过快、遇阻遇卡等原因常会导致采集所得到的声波阵列数据中出现数据缺失或数据错误的问题，即出现坏道数据，坏道数据会对后续声波数据处理解释带来不利影响，而重新采集数据也会使成本大幅增加，因此，本发明可选实施方式提出一种声波测井坏道自适应识别与重建的新型技术流程，该方法包括针对性的坏道数据自动识别以及基于凸集投影技术的声波数据自适应重建等。
图2是根据本发明可选实施方式的声波测井坏道自适应识别与重建方法的流程图，如图2所示，本发明可选实施方式具体包括以下步骤：
(1)执行声波测井(如数字声波测井、变密度声波测井、阵列声波测井等)采集，从某一深度开始，取其原始数据X(t)输入分析流程；
(2)本发明可选实施方式通过判断波形数据中明显不一致的波形数据来进行坏道数据识别，由于坏道波形数据的振幅、形状等都会与正常数据有一定差异，所以通过选取参照波形数据，并将声波测井所得的每道波形数据沿深度方向逐一与参照数据计算差异系数，通过差异系数大小来判定波形振幅与参考波形振幅的差距。具体操作为将(1)所得数据X(t)，进行快速傅里叶变换得到频谱数据S(k)，为简化参照波形数据选取流程，本发明可选实施方式将所得傅里叶变换后的数据沿深度方向，每L个深度点设置为一个窗口，取窗口中间数据作为参照数据rf(k)；
(3)定义参照数据rf(k)以后，将窗口L内每一个深度数据与rf(k)做比较，计算差异系数，计算公式如下所示：
其中，Si(k)是波形函数，rfw(k)是选取的参照波形数据，Ci是差异系数。
图3是根据本发明可选实施方式的差异系数分布示意图，如图3所示，通过差异系数Ci大小可以判定波形振幅与参考波形振幅的差距，可以看出在该图中Ci高值主要分布在深度4940m-5060m之间；
(4)分析步骤(3)差异系数Ci的大小分布，根据实际情况选取合适的坏道判定值m，声波测井数据坏道识别中m一般取0.2-4；如果Ci>m，说明该数据与参照数据振幅具有较大的差异，则将该数据判定为坏道数据。
图4是根据本发明可选实施方式的坏道识别结果示意图，如图4所示，当Ci＞m时，该数据与参照数据振幅具有较大的差异，因此将该数据中Ci＞m的部分判定为坏道数据，该图中将坏道判定值m赋值为3.5。
图5是根据本发明可选实施方式的坏道深度段置零示意图，如图5所示，坏道数据基本全被置零，这也进一步证明了本发明可选实施方式中的坏道识别技术能够较为准确的识别出坏道数据；
(5)步骤(2)(3)(4)完成了声波测井数据中的坏道数据识别过程，为了完成声波波列坏道数据的重建过程，本发明可选实施方式先根据坏道识别的结果将原始波形中的坏道数据全部置零处理，以此来消除坏道数据对声波波形数据的不利影响。具体操作为根据(4)的坏道识别结果先将原始声波波形数据X(t)中的坏道数据进行置零处理得到新的波形数据X0(t)，再对置零处理后的波形数据X0(t)进行二维傅里叶变换得到S0(k)；
(6)设置参数λ(i)，如下方公式所示，式中a表示需要重建的数据占总数据的百分比，一般取0.5-0.9；b表示步长；i＝1，2，…，a/b，a/b表示迭代次数；c表示噪声所占百分比，一般取0.01；
λ(i)＝[a，-b，c]
(7)设置阈值参数Threshod(i)如下方公式所示，判断S0(k)和Threshod(i)的大小，将小于等于Threshod的S0(k)置零，以此来消除弱振幅数据影响；
Threshod(i)＝λ(i)×max(|S0(k)|)
(8)对步骤(7)所得数据进行二维傅里叶变换逆变换，重建声波测井波形数据；
(9)重复步骤(6)、(7)、(8)，迭代计算a/b次，直至迭代到理想的精度，完成坏道数据重建。
图6是根据本发明可选实施方式的坏道深度段数据与重建结果对比图，如图6所示，图6以λ(i)＝[0.5，-0.01，0.01]为例，假设需要重建的数据占总数据的50％，噪声占比1％，循环迭代步长为0.01，迭代次数为50。原始数据与重建后的数据对比图可以看出原始数据中的坏道数据部分都得到了较好的重建效果，这也证明了该发明方法中的坏道重建技术能够较为准确的实现坏道数据重建。
根据本发明实施例，还提出了一种数据处理装置，图7是根据本发明实施例的数据处理装置的结构框图，如图7所示，该装置包括：第一获取模块71，第二获取模块72，第三获取模块73，第一确定模块74，清除模块75和第二确定模块76，下面对该装置进行说明。
第一获取模块71，用于获取目标井不同深度点对应的多个第一时域波形；第二获取模块72，连接至上述第一获取模块71，用于分别获取与多个第一时域波形对应的多个第一频域波形；第三获取模块73，连接至上述第二获取模块72，用于分别获取与多个第一频域波形对应的参照频域波形；第一确定模块74，连接至上述第三获取模块73，用于分别基于多个第一频域波形和对应的参照频域波形，确定多个第一频域波形中的异常第一频域波形；清除模块75，连接至上述第一确定模块74，用于将异常第一时域波形从多个第一时域波形中清除出去，得到多个第二时域波形；第二确定模块76，连接至上述清除模块75，用于基于多个第二时域波形，确定目标井的目标时域波形。
作为一种可选的实施例，第三获取模块73包括：划分单元，用于基于多个第一频域波形对应的深度点的顺序，将多个第一频域波形划分为多个频域波形组；第一确定单元，用于分别确定与多个第一频域波形对应的频域波形组；第二确定单元，用于分别确定与多个第一频域波形对应的频域波形组所包括的深度点中，位于中间位置处的目标深度点；第三确定单元，用于分别基于目标深度点对应的第一频域波形，确定与多个第一频域波形对应的参照频域波形。
根据本发明实施例，还提供了一种计算机可读存储介质，计算机可读存储介质包括存储的程序，其中，在程序运行时控制计算机可读存储介质所在设备执行上述任意一项的数据处理方法。
根据本发明实施例，还提供了一种计算机设备，包括：存储器和处理器，存储器存储有计算机程序；处理器，用于执行存储器中存储的计算机程序，计算机程序运行时使得处理器执行上述任意一项的数据处理方法。
上述本发明实施例序号仅仅为了描述，不代表实施例的优劣。
在本发明的上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
在本申请所提供的几个实施例中，应该理解到，所揭露的技术内容，可通过其它的方式实现。其中，以上所描述的装置实施例仅仅是示意性的，例如所述单元的划分，可以为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，单元或模块的间接耦合或通信连接，可以是电性或其它的形式。
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
另外，在本发明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。
所述集成的单元如果以软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。

Claims (10)

1.一种数据处理方法，其特征在于，包括：
获取目标井不同深度点对应的多个第一时域波形；
分别获取与所述多个第一时域波形对应的多个第一频域波形；
分别获取与所述多个第一频域波形对应的参照频域波形；
分别基于所述多个第一频域波形和对应的所述参照频域波形，确定所述多个第一频域波形中的异常第一频域波形；
将所述异常第一时域波形从所述多个第一时域波形中清除出去，得到多个第二时域波形；
基于所述多个第二时域波形，确定所述目标井的目标时域波形。
2.根据权利要求1所述的方法，其特征在于，所述分别获取与所述多个第一频域波形对应的参照频域波形，包括：
基于所述多个第一频域波形对应的深度点的顺序，将所述多个第一频域波形划分为多个频域波形组；
分别确定与所述多个第一频域波形对应的频域波形组；
分别确定与所述多个第一频域波形对应的频域波形组所包括的深度点中，位于中间位置处的目标深度点；
分别基于所述目标深度点对应的第一频域波形，确定与所述多个第一频域波形对应的参照频域波形。
3.根据权利要求2所述的方法，其特征在于，所述分别基于所述目标深度点对应的第一频域波形，确定与所述多个第一频域波形对应的参照频域波形，包括：
在所述目标深度点的个数为1的情况下，将与所述目标深度点对应的第一频域波形确定为与所述多个第一频域波形对应的参照频域波形；
在所述目标深度点的个数为2的情况下，对与所述目标深度点对应的两个第一频域波形进行平均处理，得到第一频域平均波形，并将所述第一频域平均波形确定为与所述多个第一频域波形对应的参照频域波形。
4.根据权利要求1所述的方法，其特征在于，所述分别基于所述多个第一频域波形和对应的所述参照频域波形，确定所述多个第一频域波形中的异常第一频域波形，包括：
分别获取所述多个第一频域波形和对应参照频域波形之间的差异系数；
在差异系数大于第一预定阈值的情况下，确定对应的第一频域波形为所述异常第一频域波形。
5.根据权利要求1至4中任一项所述的方法，其特征在于，所述基于所述多个第二时域波形，确定所述目标井的目标时域波形，包括：
分别获取与所述多个第二时域波形对应的多个第二频域波形；
对所述多个第二频域波形中信号幅度不大于第二预定阈值的幅度进行置零，得到多个目标频域波形；
分别确定与所述多个目标频域波形对应的多个目标时域波形。
6.根据权利要求5所述的方法，其特征在于，所述对所述第二频域波形中信号幅度不大于第二预定阈值的幅度进行置零，得到目标频域波形，包括：
采用更新所述第二预定阈值的方式，对所述第二频域波形按照预定次数循环迭代进行置零操作，得到所述目标频域波形，其中，更新所述第二预定阈值包括：基于所述第二频域波形对应的最大幅值和预定参数更新所述第二预定阈值。
7.一种数据处理装置，其特征在于，包括：
第一获取模块，用于获取目标井不同深度点对应的多个第一时域波形；
第二获取模块，用于分别获取与所述多个第一时域波形对应的多个第一频域波形；
第三获取模块，用于分别获取与所述多个第一频域波形对应的参照频域波形；
第一确定模块，用于分别基于所述多个第一频域波形和对应的所述参照频域波形，确定所述多个第一频域波形中的异常第一频域波形；
清除模块，用于将所述异常第一时域波形从所述多个第一时域波形中清除出去，得到多个第二时域波形；
第二确定模块，用于基于所述多个第二时域波形，确定所述目标井的目标时域波形。
8.根据权利要求7所述的装置，其特征在于，所述第三获取模块包括：
划分单元，用于基于所述多个第一频域波形对应的深度点的顺序，将所述多个第一频域波形划分为多个频域波形组；
第一确定单元，用于分别确定与所述多个第一频域波形对应的频域波形组；
第二确定单元，用于分别确定与所述多个第一频域波形对应的频域波形组所包括的深度点中，位于中间位置处的目标深度点；
第三确定单元，用于分别基于所述目标深度点对应的第一频域波形，确定与所述多个第一频域波形对应的参照频域波形。
9.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质包括存储的程序，其中，在所述程序运行时控制所述计算机可读存储介质所在设备执行权利要求1至6中任意一项所述的数据处理方法。
10.一种计算机设备，其特征在于，包括：存储器和处理器，
所述存储器存储有计算机程序；
所述处理器，用于执行所述存储器中存储的计算机程序，所述计算机程序运行时使得所述处理器执行权利要求1至6中任意一项所述的数据处理方法。
CN202210843496.3A 2022-07-18 2022-07-18 数据处理方法、装置、计算机可读存储介质及计算机设备 Pending CN117457024A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210843496.3A CN117457024A (zh) 2022-07-18 2022-07-18 数据处理方法、装置、计算机可读存储介质及计算机设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210843496.3A CN117457024A (zh) 2022-07-18 2022-07-18 数据处理方法、装置、计算机可读存储介质及计算机设备

Publications (1)

Publication Number Publication Date
CN117457024A true CN117457024A (zh) 2024-01-26

Family

ID=89587885

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210843496.3A Pending CN117457024A (zh) 2022-07-18 2022-07-18 数据处理方法、装置、计算机可读存储介质及计算机设备

Country Status (1)

Country Link
CN (1) CN117457024A (zh)

Similar Documents

Publication Publication Date Title
CN105323018B (zh) 用于在实时示波器上执行联合抖动和振幅噪声分析的方法
US7557572B1 (en) Method of detecting and compensating for faults in an MRI process
CN117131336B (zh) 一种电子设备连接器数据处理方法
CN111912521A (zh) 一种非平稳信号的频率检测方法和存储介质
EP3686617A1 (en) Jitter determination method and measurement instrument
CN106680873B (zh) 一种自动衡量地震数据谐波噪声强弱的振幅谱比值方法
CN112904415B (zh) 基于重构源谱一致性的q层析反演方法、系统及处理设备
CN117457024A (zh) 数据处理方法、装置、计算机可读存储介质及计算机设备
CN109471203B (zh) 提高地震资料分辨率效果评判优选方法
CN113375065B (zh) 管道泄漏监测中趋势信号的消除方法及装置
CN111401108B (zh) 斜率参数提取方法、装置及计算机可读存储介质
CN113726318A (zh) 基于wm的局部放电白噪声自适应抑制方法
CN110470396B (zh) 光谱平滑降噪的处理方法
CN111399039B (zh) 斜率参数提取方法、装置及存储介质
CN111399038B (zh) 斜率参数提取方法、装置及计算机可读存储介质
CN111007369B (zh) 一种超高频电磁波信号到达时间差计算方法和装置
CN111856559B (zh) 基于稀疏贝叶斯学习理论的多道地震谱反演方法及系统
US6873923B1 (en) Systems and methods for performing analysis of a multi-tone signal
CN110007342A (zh) 一种用于低信噪比地震信号的时频域直接拾取初至方法及系统
CN112014884A (zh) 压制近炮点强能量噪声的方法及装置
CN117668496B (zh) 一种土工格栅应变测试方法及系统
EP4000523A1 (en) Method of forming modifying data related to data sequence of data frame including electroencephalogram data, processing method of electroencephalogram data and electroencephalogram apparatus
CN113050167B (zh) 地震数据的噪声压制方法及装置
CN118300673B (zh) 一种光纤监测方法、装置、设备、存储介质及产品
CN112241024B (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