申请/专利权人：吉林大学

申请日：2023-06-04

公开（公告）日：2024-03-08

公开（公告）号：CN117665941A

主分类号：G01V1/36

分类号：G01V1/36;G01V1/28;G01V1/30;G01H17/00

优先权：

专利状态码：在审-实质审查的生效

法律状态：2024.03.26#实质审查的生效;2024.03.08#公开

摘要：本发明涉及一种基于能量突变函数包络的地震尖峰噪声簇识别与压制方法。根据地震数据中的尖峰噪声相比于目标地震信号，其持续时间更短、总能量更低，且在能量突变函数中幅值差异较大的特性，设定了半高窗阈值和能量阈值，对地震数据中的尖峰噪声进行识别，并将其压制。经实验，本方法能够在保护地震信号细节的情况下的同时有效识别并压制尖峰噪声及尖峰噪声簇，特别适用于由脚步声或风浪引起的地震数据去噪。

主权项：1.一种基于能量突变函数包络的地震尖峰噪声簇识别与压制方法，其特征在于，包括如下步骤：a.读取一组包含尖峰噪声的地震数据，记为si，i∈[1，L]，L为地震数据的采样点总数，采样率记为Fs；b.在si中采用人工拾取或其他已有自动拾取方法拾取所有地震事件，每个地震事件对应一个时间序列，将其中持续时间最长的目标地震事件记为Smax，找出Smax的起始点、终止点、和最大值对应的采样点序号分别记为Tbegin，Tmax和Tend，定义Tf＝Tmax-Tbegin，Tb＝Tend-Tmax；c.定义si的能量突变函数为 d.通过样条插值方法求取CFi的峰值包络线Bi，其中插值间距为0.1Fs个采样点；e.令j＝1；f.寻找Bkj中的最大值点，若有多个最大值点取第一个点，将最大值及其对应的采样点序号分别记录为Uj和Mj，定义Pj＝Bkj，kj∈[Mj-Tf，Mj+Tb]，记Djkj＝|Bkj-0.5Uj|，寻找Djkj的最小值及其对应的采样点序号分别记录为Aj1和幅值Cj，若有多个最小值点取第一个点，再寻找Dji的最大值Cj′，若有多个最大值点取第一个点，令Cj＝10Cj′，然后在Dji中寻找满足Mj-Aj1×Mj-Aj2＜0条件的最小值点的采样点序号Aj2，若有多个最小值点取第一个点，定义Pj的半高窗Wj＝[Aj1，Aj2]，半高窗长Zj＝|Aj1-Aj2|，求取Wj内的能量总和 将Bkj置零，再令j＝j+1；g.不断重复步骤f，直到Bkj的最大值为0，令J＝j-1，定义E＝{E1，E2…EJ}，Z＝{Z1，Z2…ZJ}，取能量阈值 和半高窗阈值 其中σ表示标准差函数；h.令p＝1，j＝1；i.对Pj判断其半高窗内能量总和与半高窗长是否满足 若不满足条件5则认为是地震信号，否则认为Pj中含有尖峰噪声，将Pj内最大值对应的采样点序号Mj记为Mp′，半高窗[Aj1，Aj2]记为[Ap1′，Ap2′]，然后令p＝p+1；j.令j＝j+1；k.不断循环步骤i至步骤j，直到j＝J+1，令P＝p-1；l.令p＝1，s′i＝si；m.令 然后令p＝p+1；n.不断循环步骤m，直到p＝P+1，此时s′i即为去除尖峰噪声后的地震数据。

全文数据：

权利要求：

百度查询： 吉林大学 基于能量突变函数包络的地震尖峰噪声簇识别与压制方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD