CN114063177B - 一种大地电磁数据去噪方法及系统 - Google Patents

一种大地电磁数据去噪方法及系统 Download PDF

Info

Publication number
CN114063177B
CN114063177B CN202111357256.4A CN202111357256A CN114063177B CN 114063177 B CN114063177 B CN 114063177B CN 202111357256 A CN202111357256 A CN 202111357256A CN 114063177 B CN114063177 B CN 114063177B
Authority
CN
China
Prior art keywords
magnetotelluric
time sequence
field component
determining
values
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
CN202111357256.4A
Other languages
English (en)
Other versions
CN114063177A (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.)
Beijing Research Institute of Uranium Geology
Original Assignee
Beijing Research Institute of Uranium Geology
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 Beijing Research Institute of Uranium Geology filed Critical Beijing Research Institute of Uranium Geology
Priority to CN202111357256.4A priority Critical patent/CN114063177B/zh
Publication of CN114063177A publication Critical patent/CN114063177A/zh
Application granted granted Critical
Publication of CN114063177B publication Critical patent/CN114063177B/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
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
    • 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

Abstract

本发明涉及一种大地电磁数据去噪方法及系统。该方法包括根据大地电磁数据获取大地电磁时间序列；根据大地电磁时间序列确定误差椭球体；根据所述误差椭球体对所述大地电磁时间序列进行离群点剔除，并修复大地电磁时间序列；根据修复后的大地电磁时间序列重构大地电磁数据。本发明能够有效去除大地电磁数据中的噪声并提高数据处理效率。

Description

一种大地电磁数据去噪方法及系统
技术领域
本发明涉及大地电磁数据预处理技术领域，特别是涉及一种大地电磁数据去噪方法及系统。
背景技术
随着社会的发展，人为活动增多，测区附近矿山、采石场、高压线、电气设备、雷达站、无线电台分布密集，大地电磁噪声类型复杂多样，且能量幅值强、相关性好、频率覆盖范围极宽，现有频率域处理的方法对该类强噪声干扰无能为力。大地电磁噪声问题一直困扰着广大大地电磁研究者,如何消除大地电磁信号中的噪声干扰,提高大地电磁测深数据质量,是一大难点。
大地电磁去噪的方法总的来说分为频率域去噪方法和时间域去噪方法,频率域去噪方法中的短时傅里叶变换通过采用滑动窗口截取信号进行傅里叶变换,从而获得任意时刻的频谱,但短时Fourier变换不能在时间域和频率域两个方向同时获得最高的分辨能力。因此,它难以满足对非平稳性的大地电磁信号进行高精度分析的要求。
对时间序列进行去噪也是大地电磁去噪的一个重要方向。Jones(1989)用中值滤波对时间域中持续时间较短的脉冲进行修复。Egbert(1997)用多脉冲响应函数与相应磁道卷积得到的预测数据来填补时间域电场中向下的脉冲。Hattingh(1989)提出了处理连续的宽频谱噪声的相关自适应去噪方法。汤井田和李晋(2010)在总结矿集区强干扰规律的基础上，提出用广义形态滤波器对大地电磁时间域信号中的矿集区强干扰波形进行提取和分离，处理结果表明广义形态滤波器能有效提取噪声波形，在一定程度上改善受噪声污染的视电阻率曲线。然而，广义形态滤波器在提取噪声波形的同时，也提取了时间域曲线中的缓变化，造成低频信息的丢失。针对这一问题，汤井田(2012)提出一种基于数学形态学top-hat变换的大地电磁时间域噪声压制方案，利用top-hat变换对波峰和波谷的检测能力，采用直线型结构元素，对大地电磁时间域信号进行去噪，去噪结果表明，所提方法能够去除噪声波形并保留时间域信号的缓变化，但该方法需要针对不同宽度的波形选取不同长度的结构元素，对含有复杂波形的时间序列难以将所有噪声去除。
针对电磁干扰严重的问题，亟需一种新的检测大地电磁时间序列中电磁干扰的方法或系统。
发明内容
本发明的目的是提供一种大地电磁数据去噪方法及系统，能够自动检测大地电磁时间序列中的离群值，计算方法简便高效，对后期视电阻率和相位曲线的修复效果明显，进而能够有效去除大地电磁数据中的噪声并提高数据处理效率。
为实现上述目的，本发明提供了如下方案：
一种大地电磁数据去噪方法，包括：
根据大地电磁数据获取大地电磁时间序列；所述大地电磁时间序列的大小为n*m；其中，m为场分量的个数；n为采样点数；
根据大地电磁时间序列确定误差椭球体；
根据所述误差椭球体对所述大地电磁时间序列进行离群点剔除，并修复大地电磁时间序列；
根据修复后的大地电磁时间序列重构大地电磁数据。
可选地，所述根据大地电磁时间序列确定误差椭球体，具体包括：
提取所述大地电磁时间序列中每个场分量的中位数值；
根据每个场分量的中位数值确定m维的误差椭球体在对应的场分量的原点以及轴长度；
根据场分量的原点以及轴长度确定m维的误差椭球体。
可选地，所述根据每个场分量的中位数值确定m维的误差椭球体在对应的场分量的原点以及轴长度，具体包括：
利用公式Am＝FmMAD(Rm)或公式Am＝Fmσ(Rm)确定轴长度；
其中，Am为第m个场分量的轴长度，MAD为中位数值绝对偏差，σ为标准偏差，Fm为第m个场分量的误差因子，Rm为第m个场分量的时间序列值。
可选地，所述根据所述误差椭球体对所述大地电磁时间序列进行离群点剔除，并修复大地电磁时间序列，具体包括：
根据大地电磁时间序列中的一组值以及相应的轴长度确定比较比值；
若比较比值大于1，则大地电磁时间序列中的一组值位于误差椭球体之外，该组值剔除；若比较比值小于1，则大地电磁时间序列中的一组值位于误差椭圆之内，该组值保留。
一种大地电磁数据去噪系统，包括：
大地电磁时间序列获取模块，用于根据大地电磁数据获取大地电磁时间序列；所述大地电磁时间序列的大小为n*m；其中，m为场分量的个数；n为采样点数；
误差椭球体确定模块，用于根据大地电磁时间序列确定误差椭球体；
离群点剔除模块，用于根据所述误差椭球体对所述大地电磁时间序列进行离群点剔除，并修复大地电磁时间序列；
大地电磁数据重构模块，用于根据修复后的大地电磁时间序列重构大地电磁数据。
可选地，所述误差椭球体确定模块具体包括：
中位数值提取单元，用于提取所述大地电磁时间序列中每个场分量的中位数值；
原点以及轴长度确定单元，用于根据每个场分量的中位数值确定m维的误差椭球体在对应的场分量的原点以及轴长度；
m维的误差椭球体确定单元，用于根据场分量的原点以及轴长度确定m维的误差椭球体。
可选地，所述原点以及轴长度确定单元具体包括：
轴长度确定子单元，用于利用公式Am＝FmMAD(Rm)或公式Am＝Fmσ(Rm)确定轴长度；
其中，Am为第m个场分量的轴长度，MAD为中位数值绝对偏差，σ为标准偏差，Fm为第m个场分量的误差因子，Rm为第m个场分量的时间序列值。
可选地，所述离群点剔除模块具体包括：
比较比值确定单元，用于根据大地电磁时间序列中的一组值以及相应的轴长度确定比较比值；
若比较比值大于1，则大地电磁时间序列中的一组值位于误差椭球体之外，该组值剔除；若比较比值小于1，则大地电磁时间序列中的一组值位于误差椭圆之内，该组值保留。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明所提供的一种大地电磁数据去噪方法及系统，根据大地电磁时间序列确定误差椭球体，根据所述误差椭球体对所述大地电磁时间序列进行离群点剔除，并修复大地电磁时间序列；即对记录的时间序列的不同统计学特性来定义一个误差椭球，从而处理大地电磁数据噪声干扰问题。自动检测大地电磁时间序列中的离群值，计算方法简便高效，对后期视电阻率和相位曲线的修复效果明显。并且，利用误差椭球体对所述大地电磁时间序列进行离群点剔除加快了对大地电磁时间序列中离群点剔除的速度，可以极大地提高数据处理效率。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明所提供的一种大地电磁数据去噪方法流程示意图；
图2为离群点剔除前的大地电磁时间序列示意图；
图3为大地电磁数据可视化示意图；
图4为处理前后的时间序列对比示意图；
图5为原始数据的视电阻率和相位曲线示意图；
图6为处理后数据的视电阻率和相位曲线示意图；
图7为本发明所提供的一种大地电磁数据去噪系统结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种大地电磁数据去噪方法及系统，能够有效去除大地电磁数据中的噪声并提高数据处理效率。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明所提供的一种大地电磁数据去噪方法流程示意图，如图1所示，本发明所提供的一种大地电磁数据去噪方法，包括：
S101，根据大地电磁数据获取大地电磁时间序列；所述大地电磁时间序列的大小为n*m；其中，m为场分量的个数；n为采样点数；场分量Ex、Ey、Hx、Hy、Hz数据都含不相干噪声；
S102，根据大地电磁时间序列确定误差椭球体；
S102具体包括：
提取所述大地电磁时间序列中每个场分量的中位数值；
根据每个场分量的中位数值确定m维的误差椭球体在对应的场分量的原点以及轴长度；
根据场分量的原点以及轴长度确定m维的误差椭球体。
所述根据每个场分量的中位数值确定m维的误差椭球体在对应的场分量的原点以及轴长度，具体包括：
利用公式Am＝FmMAD(Rm)或公式Am＝Fmσ(Rm)确定轴长度；
其中，Am为第m个场分量的轴长度，MAD为中位数值绝对偏差，σ为标准偏差，Fm为第m个场分量的误差因子，Rm为第m个场分量的时间序列值。
S103，根据所述误差椭球体对所述大地电磁时间序列进行离群点剔除，并修复大地电磁时间序列；
S103具体包括：
根据大地电磁时间序列中的一组值以及相应的轴长度确定比较比值；
若比较比值大于1，则大地电磁时间序列中的一组值位于误差椭球体之外，该组值剔除；若比较比值小于1，则大地电磁时间序列中的一组值位于误差椭圆之内，该组值保留。
根据椭球体的这些特性估计一组值(Sn,1,...,Sn,m)是位于椭圆内部还是外部，即：
其中当Cn＞1位于误差椭球体之外，被归类为异常值，反之，当Cn＜1位于误差椭圆之内，被归类为正常值。
建议使用值Am＝3σ或Am＝4MAD，以达到近95％的置信区间。但是，在很多情况下，需要通过反复尝试并结合处理结果来验证选定的参数是否合理。Fm值的调整是需要用户来完成。
S104，根据修复后的大地电磁时间序列重构大地电磁数据。
作为一个具体的实施例，一个大地电磁数据集包括电道记录(Ex，Ey)和磁道记录(Hx，Hy)，采样率为18.315赫兹，该数据集中每个分量包括6177600个采样数据。
如图2所示，各电、磁道的时间序列均出现了很多饱和信号。如果通过肉眼检查并进行手动编辑，耗时很长并且极易出错。
利用本发明所提供的方法对数据进行处理，数据绘图将对应于一个四维误差椭球，很难可视化。因此，将可视化分为两个图形，一个用于Ex、Ey，另一个用于Hx、Hy(图3)。使用参数Fm＝1.5σ(对于Ex、Ey)和Fm＝2σ(对于Hx、Hy)进行处理时，该数据的97.65％被分类为非离群值。
时间序列处理后结果显示(图4)，大地电磁电、磁道时间序列中的离群值被剔除。
处理前后的时间序列使用同一软件、选取相同的参数进行处理。图5((a)中横坐标:秒，纵坐标:欧姆米；(b)中横坐标:秒，纵坐标:度)显示了原始数据的视电阻率和相位曲线。可以看到，原始数据小于100秒的视电阻率和相位曲线比较平滑，但是低频数据质量较差。图6((a)中横坐标:秒，纵坐标:欧姆米；(b)中横坐标:秒，纵坐标:度)显示了经过剔除离群点数值处理的结果，得到的曲线数据效果更好，低频可以达到1000秒。这说明了算法是实用的、有效的。
图7为本发明所提供的一种大地电磁数据去噪系统结构示意图，如图7所示，本发明所提供的一种大地电磁数据去噪系统，包括：
大地电磁时间序列获取模块701，用于根据大地电磁数据获取大地电磁时间序列；所述大地电磁时间序列的大小为n*m；其中，m为场分量的个数；n为采样点数；
误差椭球体确定模块702，用于根据大地电磁时间序列确定误差椭球体；
离群点剔除模块703，用于根据所述误差椭球体对所述大地电磁时间序列进行离群点剔除，并修复大地电磁时间序列；
大地电磁数据重构模块704，用于根据修复后的大地电磁时间序列重构大地电磁数据。
所述误差椭球体确定模块702具体包括：
中位数值提取单元，用于提取所述大地电磁时间序列中每个场分量的中位数值；
原点以及轴长度确定单元，用于根据每个场分量的中位数值确定m维的误差椭球体在对应的场分量的原点以及轴长度；
m维的误差椭球体确定单元，用于根据场分量的原点以及轴长度确定m维的误差椭球体。
所述原点以及轴长度确定单元具体包括：
轴长度确定子单元，用于利用公式Am＝FmMAD(Rm)或公式Am＝Fmσ(Rm)确定轴长度；
其中，Am为第m个场分量的轴长度，MAD为中位数值绝对偏差，σ为标准偏差，Fm为第m个场分量的误差因子，Rm为第m个场分量的时间序列值。
所述离群点剔除模块703具体包括：
比较比值确定单元，用于根据大地电磁时间序列中的一组值以及相应的轴长度确定比较比值；
若比较比值大于1，则大地电磁时间序列中的一组值位于误差椭球体之外，该组值剔除；若比较比值小于1，则大地电磁时间序列中的一组值位于误差椭圆之内，该组值保留。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (5)

1.一种大地电磁数据去噪方法，其特征在于，包括：
根据大地电磁数据获取大地电磁时间序列；所述大地电磁时间序列的大小为n*m；其中，m为场分量的个数；n为采样点数；
根据大地电磁时间序列确定误差椭球体；
根据所述误差椭球体对所述大地电磁时间序列进行离群点剔除，并修复大地电磁时间序列；
根据修复后的大地电磁时间序列重构大地电磁数据；
所述根据大地电磁时间序列确定误差椭球体，具体包括：
提取所述大地电磁时间序列中每个场分量的中位数值；
根据每个场分量的中位数值确定m维的误差椭球体在对应的场分量的原点以及轴长度；
根据场分量的原点以及轴长度确定m维的误差椭球体；
所述根据每个场分量的中位数值确定m维的误差椭球体在对应的场分量的原点以及轴长度，具体包括：
利用公式Am＝FmMAD(Rm)或公式Am＝Fmσ(Rm)确定轴长度；
其中，Am为第m个场分量的轴长度，MAD为中位数值绝对偏差，σ为标准偏差，Fm为第m个场分量的误差因子，Rm为第m个场分量的时间序列值；
所述根据所述误差椭球体对所述大地电磁时间序列进行离群点剔除，并修复大地电磁时间序列，具体包括：
根据大地电磁时间序列中的一组值以及相应的轴长度确定比较比值；
若比较比值大于1，则大地电磁时间序列中的一组值位于误差椭球体之外，该组值剔除；若比较比值小于1，则大地电磁时间序列中的一组值位于误差椭圆之内，该组值保留。
2.一种大地电磁数据去噪系统，用于实施如权利要求1所述的一种大地电磁数据去噪方法；其特征在于，包括：
大地电磁时间序列获取模块，用于根据大地电磁数据获取大地电磁时间序列；所述大地电磁时间序列的大小为n*m；其中，m为场分量的个数；n为采样点数；
误差椭球体确定模块，用于根据大地电磁时间序列确定误差椭球体；
离群点剔除模块，用于根据所述误差椭球体对所述大地电磁时间序列进行离群点剔除，并修复大地电磁时间序列；
大地电磁数据重构模块，用于根据修复后的大地电磁时间序列重构大地电磁数据。
3.根据权利要求2所述的一种大地电磁数据去噪系统，其特征在于，所述误差椭球体确定模块具体包括：
中位数值提取单元，用于提取所述大地电磁时间序列中每个场分量的中位数值；
原点以及轴长度确定单元，用于根据每个场分量的中位数值确定m维的误差椭球体在对应的场分量的原点以及轴长度；
m维的误差椭球体确定单元，用于根据场分量的原点以及轴长度确定m维的误差椭球体。
4.根据权利要求3所述的一种大地电磁数据去噪系统，其特征在于，所述原点以及轴长度确定单元具体包括：
轴长度确定子单元，用于利用公式Am＝FmMAD(Rm)或公式Am＝Fmσ(Rm)确定轴长度；
其中，Am为第m个场分量的轴长度，MAD为中位数值绝对偏差，σ为标准偏差，Fm为第m个场分量的误差因子，Rm为第m个场分量的时间序列值。
5.根据权利要求2所述的一种大地电磁数据去噪系统，其特征在于，所述离群点剔除模块具体包括：
比较比值确定单元，用于根据大地电磁时间序列中的一组值以及相应的轴长度确定比较比值；
若比较比值大于1，则大地电磁时间序列中的一组值位于误差椭球体之外，该组值剔除；若比较比值小于1，则大地电磁时间序列中的一组值位于误差椭圆之内，该组值保留。
CN202111357256.4A 2021-11-16 2021-11-16 一种大地电磁数据去噪方法及系统 Active CN114063177B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111357256.4A CN114063177B (zh) 2021-11-16 2021-11-16 一种大地电磁数据去噪方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111357256.4A CN114063177B (zh) 2021-11-16 2021-11-16 一种大地电磁数据去噪方法及系统

Publications (2)

Publication Number Publication Date
CN114063177A CN114063177A (zh) 2022-02-18
CN114063177B true CN114063177B (zh) 2023-09-29

Family

ID=80273486

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111357256.4A Active CN114063177B (zh) 2021-11-16 2021-11-16 一种大地电磁数据去噪方法及系统

Country Status (1)

Country Link
CN (1) CN114063177B (zh)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1995003557A1 (en) * 1993-07-21 1995-02-02 Western Atlas International, Inc. Method of determining formation resistivity utilizing combined measurements of inductive and galvanic logging instruments
WO2001052218A1 (en) * 2000-01-11 2001-07-19 Solipsys Corporation Method and system for sharing vehicle telemetry data among a plurality of users over a communications network
JP2002328162A (ja) * 2001-05-01 2002-11-15 Mitsubishi Electric Corp 目標追尾装置
CN101825723A (zh) * 2004-03-17 2010-09-08 维斯特恩格科地震控股有限公司 海上地震测量方法和系统
EP2327994A1 (en) * 2009-11-27 2011-06-01 BAE Systems PLC Improvements in or relating to tracking radio signal sources
AU2015213395A1 (en) * 2009-03-17 2015-09-10 Schlumberger Technology B.V. Relative and absolute error models for subterranean wells
CN107870001A (zh) * 2017-11-27 2018-04-03 东南大学 一种基于椭球拟合的磁力计校正方法
US9945970B1 (en) * 2011-08-29 2018-04-17 Seismic Innovations Method and apparatus for modeling microseismic event location estimate accuracy
CN111158045A (zh) * 2020-01-06 2020-05-15 中国石油化工股份有限公司 油藏改造微地震事件散点聚类分析方法及系统
WO2021002761A1 (en) * 2019-07-02 2021-01-07 Equinor Energy As Improved inversions of geophysical data

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6873265B2 (en) * 2001-09-14 2005-03-29 Quakefinder Llc Satellite and ground system for detection and forecasting of earthquakes

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO1995003557A1 (en) * 1993-07-21 1995-02-02 Western Atlas International, Inc. Method of determining formation resistivity utilizing combined measurements of inductive and galvanic logging instruments
WO2001052218A1 (en) * 2000-01-11 2001-07-19 Solipsys Corporation Method and system for sharing vehicle telemetry data among a plurality of users over a communications network
JP2002328162A (ja) * 2001-05-01 2002-11-15 Mitsubishi Electric Corp 目標追尾装置
CN101825723A (zh) * 2004-03-17 2010-09-08 维斯特恩格科地震控股有限公司 海上地震测量方法和系统
AU2015213395A1 (en) * 2009-03-17 2015-09-10 Schlumberger Technology B.V. Relative and absolute error models for subterranean wells
EP2327994A1 (en) * 2009-11-27 2011-06-01 BAE Systems PLC Improvements in or relating to tracking radio signal sources
US9945970B1 (en) * 2011-08-29 2018-04-17 Seismic Innovations Method and apparatus for modeling microseismic event location estimate accuracy
CN107870001A (zh) * 2017-11-27 2018-04-03 东南大学 一种基于椭球拟合的磁力计校正方法
WO2021002761A1 (en) * 2019-07-02 2021-01-07 Equinor Energy As Improved inversions of geophysical data
CN111158045A (zh) * 2020-01-06 2020-05-15 中国石油化工股份有限公司 油藏改造微地震事件散点聚类分析方法及系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴曲波 ; 张志勇 ; 柯丹 ; 范利飞 ; .点源二维直流电阻率法正演模拟.工程地球物理学报.2014,(01),全文. *
王福昌 ; 万永革 ; 钱小仕 ; 张丽娟 ; 张梅东 ; .由地震分布丛集性给出断层参数的一种新方法.地球物理学报.2013,(02),全文. *

Also Published As

Publication number Publication date
CN114063177A (zh) 2022-02-18

Similar Documents

Publication Publication Date Title
Zhong et al. Partial discharge signal denoising based on singular value decomposition and empirical wavelet transform
Song et al. Second generation wavelet transform for data denoising in PD measurement
CN108458871B (zh) 一种基于改进经验小波变换的齿轮箱故障识别方法
CN114492538B (zh) 一种城市中压配电电缆局部放电信号去噪方法
CN111308285B (zh) 一种窄带干扰降噪方法
CN112014773B (zh) 一种用于检测小电流接地系统电缆早期故障的方法
Abdullah et al. Power quality analysis using bilinear time-frequency distributions
Zhong et al. An improved wavelet spectrum segmentation algorithm based on spectral kurtogram for denoising partial discharge signals
CN113325277A (zh) 一种局部放电处理方法
CN108802175A (zh) 一种基于emd的脉冲涡流检测信号降噪方法
CN114063177B (zh) 一种大地电磁数据去噪方法及系统
CN107610055B (zh) 傅里叶变换光谱仪干涉图噪声检测及抑制方法
CN111239565B (zh) 基于分层式去噪模型的充油套管局部放电脉冲信号处理方法及系统
CN111553308A (zh) 一种电力变压器局部放电信号的重构方法
CN117169664A (zh) 一种局放传感器抗干扰系统及方法
Zhang et al. Application of signal processing techniques to on-line partial discharge detection in cables
CN103078661A (zh) 一种基于迭代门限的扩频系统干扰抑制方法
CN113702821A (zh) 一种提取gis局部放电信号的方法及系统
Li et al. A Supraharmonic measurement method based on matrix pencil with high-order mixed mean cumulant
CN112817056A (zh) 一种大地电磁信号去噪方法及系统
CN104361892A (zh) 一种干扰信号与语音调制信号混叠的窄带干扰识别方法
Zhang et al. Simultaneous denoising and preserving of seismic signals by multiscale time-frequency peak filtering
CN117520916A (zh) 一种电缆局部放电信号去噪方法
CN110569812B (zh) 一种故障信号的包络解调方法及包络解调系统
Hu et al. A new technique for extracting partial discharge signals in on-line monitoring with wavelet analysis

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