一种新的石油地震勘探数据去噪方法
【技术领域】
[0001] 本发明属于地震数据勘探技术领域，涉及一种新的石油地震勘探数据去噪方法。
【背景技术】
[0002] 石油勘探地震数据去噪，是石油地震勘探数据处理中非常重要的一个环节。去噪 方法一般是根据不同的变换域中信号和噪声的差异进行。常用的变换域有f_k域，f-x域， τ-P域和频域等。如果噪声和信号在所述变换域中差异不明显，或者无差异，则去噪无法进 行。例如，经过常用的低通滤波器（可以选用Butterworth滤波器实现）后，噪声和信号在 频域将无法再进一步分离，但数据中仍然存在噪声，只不过此时噪声和信号（数据）具有相 同的振幅谱。
【发明内容】
[0003] 本发明的目的在于提供一种新的石油地震勘探数据去噪方法，解决了现有技术常 用的低通滤波器噪声和信号在频域将无法再进一步分离，数据中仍然存在噪声的问题。
[0004] 本发明所采用的技术方案是按照以下步骤进行：
[0005] 步骤1 :首先采集得到的石油地震勘探数据，经过本领域常规的预处理；
[0006] 步骤2 :然后利用快速傅立叶变换分析地震数据的频谱，确定有用信号的频谱范 围，进而明确低通滤波器的截止频率，利用此截止频率对地震数据进行低通滤波，滤除地震 数据的高频噪声；
[0007] 步骤3 :对于滤除高频噪声之后的地震数据，送入分数傅立叶变换（FrFT)域进行 变换，在进行FrFT变换时，首先按照一定的时间间隔扫描将其分为不同的阶数，在分别计 算不同阶数所对应的FrFT变换；
[0008] 步骤4 :根据扫描所得到的FrFT域变换，选择出信号最为集中的变换数据和其所 对应的阶数Btl;
[0009] 步骤5 :将步骤3中进行FrFT域变换后的数据，再次进行反变换，反变换的阶数 为-?，经此变换之后的数据即为所求数据。
[0010] 进一步，所述步骤2中的低通滤波方式采用巴特沃斯滤波器实现。
[0011] 进一步，所述步骤3中FrFT变换的阶数范围为[0,4]，时间间隔小于0.5秒。
[0012] 本发明的有益效果是石油勘探数据中的噪声的进一步滤除。
【附图说明】
[0013] 图1是本发明方法流程示意图；
[0014] 图2是Ricker子波的FrFT变换图。
【具体实施方式】
[0015] 下面结合附图和【具体实施方式】对本发明进行详细说明。
[0016] 如图1所示，本发明方法按照以下步骤进行：
[0017] 步骤1 :首先在野外采集得到的石油地震勘探数据，经过本领域常规的预处理（预 处理可采用数据解编、道编辑等处理方法）；
[0018] 步骤2 :利用FFT (快速傅立叶变换）分析地震数据的频谱，确定有用信号的频谱 范围，进而明确低通滤波器的截止频率，利用此截止频率对地震数据进行低通滤波（可以 利用巴特沃斯滤波器实现，也可以选择其他的实现形式），滤除地震数据的高频噪声；
[0019] 步骤3 :对于滤除高频噪声之后的地震数据，送入分数傅立叶变换（FrFT)域进行 变换，FrFT变换的阶数范围为[0,4]。在进行FrFT变换时，首先进行变换扫描，即在[0,4] 范围内按照一定的时间间隔（时间间隔优选小于0.5秒）分别计算不同阶数所对应的FrFT 变换，其计算公式为：
【主权项】
1. 一种新的石油地震勘探数据去噪方法，其特征在于，该新的石油地震勘探数据去噪 方法包括以下步骤： 步骤一：首先采集得到的石油地震勘探数据，预处理； 步骤二：然后利用快速傅立叶变换分析地震数据的频谱，确定有用信号的频谱范围，进 而明确低通滤波器的截止频率，利用此截止频率对地震数据进行低通滤波，滤除地震数据 的尚频噪声； 步骤三：对于滤除高频噪声之后的地震数据，送入分数傅立叶变换FrFT域进行变换， 在进行FrFT变换时，首先按照时间间隔扫描将其分为不同的阶数，在分别计算不同阶数所 对应的FrFT变换； 步骤四：根据扫描所得到的FrFT域变换，选择出信号最为集中的变换数据和其所对应 的阶数aQ; 步骤五：将步骤三中进行FrFT域变换后的数据，再次进行反变换，反变换的阶数 为-%，经此变换之后的数据即为所求数据。
2. 如权利要求1所述新的石油地震勘探数据去噪方法，其特征在于：所述步骤二中的 低通滤波方式采用巴特沃斯滤波器实现。
3. 如权利要求1所述新的石油地震勘探数据去噪方法，其特征在于：步骤三中FrFT变 换的阶数范围为[〇,4]，时间间隔小于0.5秒。
【专利摘要】本发明公开了一种新的石油地震勘探数据去噪方法，采集得到的石油地震勘探数据，利用快速傅立叶变换分析地震数据的频谱，确定有用信号的频谱范围，进而明确低通滤波器的截止频率，利用此截止频率对地震数据进行低通滤波，滤除地震数据的高频噪声；对于滤除高频噪声之后的地震数据，送入FrFT域进行变换，在进行FrFT变换时，根据扫描所得到的FrFT域变换，选择出信号最为集中的变换数据和其所对应的阶数a0；将FrFT域变换后的数据，再次进行反变换，反变换的阶数为-a0，经此变换之后的数据即为所求数据。发明的有益效果是石油勘探数据中的噪声的进一步滤除。
【IPC分类】G01V1-36
【公开号】CN104597502
【申请号】CN201410735187
【发明人】翟明岳
【申请人】翟明岳
【公开日】2015年5月6日
【申请日】2014年12月8日