基于一种改进小波阈值的微震信号降噪方法(MATLAB)

2024-05-10 15:36

本文主要是介绍基于一种改进小波阈值的微震信号降噪方法(MATLAB),希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

微震是指岩体由于在人为扰动或自然原因下受力变形,发生破裂过程中能量积聚而释放的弹性波或应力波。微震信号具有信噪比低、不稳定性、瞬时性和多样性等特点。因此,在任何损坏之前都会出现微小的裂缝,这种微小的裂缝是由岩层中应力和应变的变化引起的。应力和应变发生变化后,会以弹性能量释放的形式产生弹性波,并能被传感器接收。因此,利用微震监测技术在不同空间方位上布设的传感器接收到岩体产生的弹性波信息,对波形事件信息进行分析,处理后确定微震事件的发生事件、空间位置、数量、严重程度以及能量释放等信息,从而推断岩石从微观应变到宏观破裂的发展趋势进而对岩爆灾害进行动态预警。

由于实际微震监测工程环境的复杂性,当微震事件发生时,往往带有背景噪声,甚至噪声淹没了微震信息。对于带有噪声的微震数据往往会舍弃,因为直接使用这部分数据会对后续的相位初至拾取、震源定位、微震信号特征提取、震源机理解释等工作产生很大地影响,这也就降低了数据的利用率。在波形分类过程中发现这部分数据占有较大比重,采集到的纯微震信号也可能含有噪音,为了更大程度的利用监测到的数据,所以对含噪微震信号进行有针对性的降噪是非常必要的。同时对信噪比比较低的未知信号波形是否含有微震信息也可以通过降噪去判断。

鉴于此,提出一种基于一种改进小波阈值的微震信号降噪方法,该方法采用交叉验证方法确定阈值,不依赖于噪声的先验知识,部分代码如下:

function [psih] = wfilth(type, N, a, opt)
% Outputs the FFT of the wavelet of family 'type' with parameters
% in 'opt', of length N at scale a: (psi(-t/a))^.
%
% [Inputs]
%   type: wavelet type 
%   N: number of samples to calculate
%   a: wavelet scale parameter 
%   opt: wavelet options 
%   opt.dt: delta t 
%
% [Outputs]
%   psih: wavelet sampling in frequency domain 
%---------------------------------------------------------------------------------opt = struct(); k = 0:(N-1);xi = zeros(1, N);xi(1:N/2+1) = 2*pi/N*[0:N/2];xi(N/2+2:end) = 2*pi/N*[-N/2+1:-1];psihfn = wfiltfn(type, opt);psih = psihfn(a*xi);% Normalizing psih = psih * sqrt(a) / sqrt(2*pi);% Center around zero in the time domainpsih = psih .* (-1).^k;
%完整代码可由知乎学术咨询获得:https://www.zhihu.com/consult/people/792359672131756032

出图如下:

完整代码可由知乎学术咨询获得:

工学博士,担任《Mechanical System and Signal Processing》《中国电机工程学报》等期刊审稿专家,擅长领域:现代信号处理,机器学习/深度学习,时间序列分析/预测,电气设备(电机,变压器等)表面缺陷检测、电气设备(电机,变压器等)异常检测、电气设备(电机,变压器等)故障诊断与健康管理PHM、电气设备(电机,变压器等)剩余使用寿命预测等。

这篇关于基于一种改进小波阈值的微震信号降噪方法(MATLAB)的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/976850

相关文章

Nginx设置连接超时并进行测试的方法步骤

《Nginx设置连接超时并进行测试的方法步骤》在高并发场景下,如果客户端与服务器的连接长时间未响应,会占用大量的系统资源,影响其他正常请求的处理效率,为了解决这个问题,可以通过设置Nginx的连接... 目录设置连接超时目的操作步骤测试连接超时测试方法:总结:设置连接超时目的设置客户端与服务器之间的连接

Java判断多个时间段是否重合的方法小结

《Java判断多个时间段是否重合的方法小结》这篇文章主要为大家详细介绍了Java中判断多个时间段是否重合的方法,文中的示例代码讲解详细,感兴趣的小伙伴可以跟随小编一起学习一下... 目录判断多个时间段是否有间隔判断时间段集合是否与某时间段重合判断多个时间段是否有间隔实体类内容public class D

Python使用国内镜像加速pip安装的方法讲解

《Python使用国内镜像加速pip安装的方法讲解》在Python开发中,pip是一个非常重要的工具,用于安装和管理Python的第三方库,然而,在国内使用pip安装依赖时,往往会因为网络问题而导致速... 目录一、pip 工具简介1. 什么是 pip?2. 什么是 -i 参数?二、国内镜像源的选择三、如何

IDEA编译报错“java: 常量字符串过长”的原因及解决方法

《IDEA编译报错“java:常量字符串过长”的原因及解决方法》今天在开发过程中,由于尝试将一个文件的Base64字符串设置为常量,结果导致IDEA编译的时候出现了如下报错java:常量字符串过长,... 目录一、问题描述二、问题原因2.1 理论角度2.2 源码角度三、解决方案解决方案①:StringBui

Linux使用nload监控网络流量的方法

《Linux使用nload监控网络流量的方法》Linux中的nload命令是一个用于实时监控网络流量的工具,它提供了传入和传出流量的可视化表示,帮助用户一目了然地了解网络活动,本文给大家介绍了Linu... 目录简介安装示例用法基础用法指定网络接口限制显示特定流量类型指定刷新率设置流量速率的显示单位监控多个

Java覆盖第三方jar包中的某一个类的实现方法

《Java覆盖第三方jar包中的某一个类的实现方法》在我们日常的开发中,经常需要使用第三方的jar包,有时候我们会发现第三方的jar包中的某一个类有问题,或者我们需要定制化修改其中的逻辑,那么应该如何... 目录一、需求描述二、示例描述三、操作步骤四、验证结果五、实现原理一、需求描述需求描述如下:需要在

JavaScript中的reduce方法执行过程、使用场景及进阶用法

《JavaScript中的reduce方法执行过程、使用场景及进阶用法》:本文主要介绍JavaScript中的reduce方法执行过程、使用场景及进阶用法的相关资料,reduce是JavaScri... 目录1. 什么是reduce2. reduce语法2.1 语法2.2 参数说明3. reduce执行过程

C#中读取XML文件的四种常用方法

《C#中读取XML文件的四种常用方法》Xml是Internet环境中跨平台的,依赖于内容的技术,是当前处理结构化文档信息的有力工具,下面我们就来看看C#中读取XML文件的方法都有哪些吧... 目录XML简介格式C#读取XML文件方法使用XmlDocument使用XmlTextReader/XmlTextWr

C++初始化数组的几种常见方法(简单易懂)

《C++初始化数组的几种常见方法(简单易懂)》本文介绍了C++中数组的初始化方法,包括一维数组和二维数组的初始化,以及用new动态初始化数组,在C++11及以上版本中,还提供了使用std::array... 目录1、初始化一维数组1.1、使用列表初始化(推荐方式)1.2、初始化部分列表1.3、使用std::

oracle DBMS_SQL.PARSE的使用方法和示例

《oracleDBMS_SQL.PARSE的使用方法和示例》DBMS_SQL是Oracle数据库中的一个强大包,用于动态构建和执行SQL语句,DBMS_SQL.PARSE过程解析SQL语句或PL/S... 目录语法示例注意事项DBMS_SQL 是 oracle 数据库中的一个强大包,它允许动态地构建和执行