一种地震数据去噪方法、装置及设备制造方法及图纸

技术编号:26341407 阅读:36 留言:0更新日期:2020-11-13 20:21
本说明书实施例提供一种地震数据去噪方法、装置及设备。所述方法包括:接收待处理地震数据;所述待处理地震数据中包含噪声数据;根据预设模型波场将待处理地震数据转化为待处理地震波场;所述预设模型波场为平面波波场;所述待处理地震波场包括至少一个采样点;所述采样点对应有至少一个局部斜率;利用局部斜率计算待处理地震波场对应于各个采样点的连续性评估值;所述连续性评估值,用于衡量待处理地震波场在采样点处的连续性;若采样点的连续性评估值小于连续性阈值,标记所述采样点为噪声数据点;基于噪声数据点去除待处理地震数据中的噪声数据。上述方法,基于地震数据的连续性,提高了地震数据的信噪比,有利于后续生产开发的进行。

A method, device and equipment for seismic data denoising

【技术实现步骤摘要】
一种地震数据去噪方法、装置及设备
本说明书实施例涉及地震信号处理
,特别涉及一种地震数据去噪方法、装置及设备。
技术介绍
在地质勘探领域中,往往通过在目标工区进行人工爆破产生地震波,并利用检波器等地震信号勘探仪器接收所述地震波,并基于所接受到的信号数据确定所述目标工区的地质构造,以便于后续生产开发工作的进行。但是,在实际应用中,由于仪器精度、电磁干扰以及外界环境的影响,使得所获取到的地震数据中往往包含有一定的噪声数据。噪声数据是与目标工区的地质构造不相干的无意义数据,在进行地下构造解释、地球物理反演等操作时会对实验结果的精确性造成一定的影响。当所获取的地震数据的信噪比较低时,噪声会对有效信号进行掩盖,使得无法对所述地震数据进行利用,从而影响后续勘探及生产工作的进行。因此,在利用所接收到的地震数据之前,需要去除所述地震数据中的噪声以获取其中的有效信号。
技术实现思路
本说明书实施例的目的是提供一种地震数据去噪方法、装置及设备,以解决如何去除地震数据中的噪声的问题。为了解决上述技术问题,本说明书实施例提出一本文档来自技高网...

【技术保护点】
1.一种地震数据去噪方法,其特征在于,包括:/n获取待处理地震数据;所述待处理地震数据中包含噪声数据;/n根据预设模型波场将所述待处理地震数据转化为待处理地震波场;所述预设模型波场为平面波波场;所述待处理地震波场包括至少一个采样点;所述采样点对应有至少一个局部斜率;/n利用所述局部斜率计算待处理地震波场对应于各个采样点的连续性评估值;所述连续性评估值,用于衡量待处理地震波场在所述采样点处的连续性;/n若采样点的连续性评估值小于连续性阈值,标记所述采样点为噪声数据点;/n基于所述噪声数据点去除所述待处理地震数据中的噪声数据。/n

【技术特征摘要】
1.一种地震数据去噪方法,其特征在于,包括:
获取待处理地震数据;所述待处理地震数据中包含噪声数据;
根据预设模型波场将所述待处理地震数据转化为待处理地震波场;所述预设模型波场为平面波波场;所述待处理地震波场包括至少一个采样点;所述采样点对应有至少一个局部斜率;
利用所述局部斜率计算待处理地震波场对应于各个采样点的连续性评估值;所述连续性评估值,用于衡量待处理地震波场在所述采样点处的连续性;
若采样点的连续性评估值小于连续性阈值,标记所述采样点为噪声数据点;
基于所述噪声数据点去除所述待处理地震数据中的噪声数据。


2.如权利要求1所述的方法,其特征在于,所述根据预设模型波场将所述待处理地震数据转化为待处理地震波场,包括:
利用公式计算预测波场,式中,为待处理地震波场,A为第一平滑窗口函数,x为待处理地震波场的位置,κ为第一平滑窗口函数的位移,P′为候选待处理地震波场,f(t-σx)表示对应于时刻t、局部斜率σ和位置x的预设模型波场。


3.如权利要求1所述的方法,其特征在于,所述利用所述局部斜率计算待处理地震波场对应于各个采样点的连续性评估值,包括:
计算预测波场在所述采样点处分别对应于各个局部斜率的斜率分布概率;所述斜率分布概率用于衡量预测波场在采样点处沿各个局部斜率分布的概率;
利用所述斜率分布概率计算对应于所述采样点的连续性评估值。


4.如权利要求3所述的方法,其特征在于,所述根据所述局部斜率计算预测波场对应于所述采样点的斜率分布概率,包括:
利用公式计算斜率分布概率,式中,Q为斜率分布概率,B为第二平滑窗口函数,t为采样点对应的时刻,τ为第二平滑窗口函数的位移,Q′为候选斜率分布概率,为预测波场。


5.如权利要求3所述的方法,其特征在于,所述利用所述斜率分布概率计算对应于所述采样点的连续性评估值,包括:
利用公式Ce=-∑σ(Q)/∫Qlogb(Q/∑σ(Q))dσ计算连续性评估值,式中,∑σ(Q)为斜率分布概率沿局部斜率方向的求和结果,Ce为连续性评估值,σ为局部斜率,Q为斜率分布概率,b为预设底数。


6.如权利要求1所述的方法,其特征在于,所述若采样点的连续性评估值小于连续性阈值,标记所述采样点为噪声数据点之前,还...

【专利技术属性】
技术研发人员:黄炜霖王润秋高菲尚国军周延鑫
申请(专利权)人:中国石油大学北京
类型:发明
国别省市:北京;11

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1