【技术实现步骤摘要】
基于FX预测滤波的基追踪算法的低信噪比地震数据优化处理方法
[0001]本专利技术涉及地震数据处理
，具体涉及一种基于FX预测滤波的基追踪算法的低信噪比地震数据优化处理方法。
技术介绍
[0002]地震反射系数反演作为一种从采集到的地震资料反演反射系数序列的反褶积方法，是地震处理中一个非常重要的流程。由于地震道的带限特性，反褶积问题通常是病态的。针对能够得到的无限多的反射系数形式，需要施加适当的约束从而获得合适的结果。在此前的研究中，柯西准则和L1范数已被广泛应用于稀疏反射系数反演。前人提出了匹配追踪算法(MP)，匹配追踪算法将信号扩展为一个冗余字典中的小波(或原子)的线性叠加。通过寻找最符合条件的小波并计算它与原始地震道之间的残差，反复迭代直至残差的能量低于预设的阈值，从而得到反射系数序列的最优解。基于多通道的匹配追踪算法(MCMP)则进一步提高了空间的连续性，并且获得了高分辨率的时频谱。在地震反射系数反演的方法中，基追踪算法是一个相对较新的算法。将该方法用于求解欠定线性方程组时，获得的解是唯一并且稀疏的。这些特性说...
【技术保护点】
【技术特征摘要】
1.一种基于FX预测滤波的基追踪算法的低信噪比地震数据优化处理方法，其特征在于,包括以下步骤：1)使用基追踪算法计算反射系数r'；2)对上述步骤1)求得的r'应用FX预测滤波得到滤波后的反射系数r；3)将得到的滤波后的反射系数r与已知的子波做褶积，从而得到合成地震记录d；4)重新计算系数向量m的初始值；5)重复基追踪算法流程，根据上述步骤4)中得到的最新初始值m，迭代最小化目标函数直到获得满意的结果为止。2.如权利要求1所述的基于FX预测滤波的基追踪算法的低信噪比地震数据优化处理方法，其特征在于：所述步骤2)中由于FX预测滤波是在频域中进行，而基追踪算法是在时间域中实现,因此在实际操作中，必须首先利用傅里叶变换将反射信号转换到频域，待预测滤波完成后，再变换回时间域。3.如权利要求1所述的基于FX预测滤波的基追踪算法的低信噪比地震数据优化处理方法，其特征在于：所述步骤2)中的FX预测滤波是基于地震信号图像由频域中一系列线性相干的反射叠加组成这一假设的一种滤波方法；在线性假设下，地震信号是严格的正弦形式，并在频域有可预测性；对某一固定频率f，某一道地震信号可以通过其相邻的地震道预测得到：D
i
(f)＝P1D'
i
‑1(f)+P2D'
i
‑2(f)+
…
+P
L
D'
i
‑
L
(f),
ꢀꢀꢀꢀꢀꢀꢀꢀ
(10)上式10)中P
i
表示长度为L的预测滤波器，D
i
为频域中经过平滑后的某道地震信号，D'
i
则表示平滑之前的该道地震信号；由于输入的地震信号和反演得到的反射系数有相似的空间连续性，根据其可预测性，反射系数可由公式(11)表示：R
i
(f)＝P1R'
i
‑1(f)+P2R'
i
‑2(f)+
…
+P
L
R'
i
‑
L
(f),
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(11)上式11)中，R
i
和R'分别表示频域中平滑之后及之前的反射系数。4.如权利要求1所述的基于FX预测滤波的基追踪算法的低信噪比地震数据优化处理方法，其特征在于，所述基追踪算法的具体流程为：假设地层为一系列具有固定速度和密度的平层，地震道可以被表示为反射系数和子波的褶积结果：d(t)＝w(t)*r(t)+n(t)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(1)；上式(1)中,d(t)为观测到的地震道，w(t)表示地震子波，r(t)为待求的反射系数，n(t)为地震道上存在的噪音；将地层的反射被看作是反射对，偶极子被用于分解地层的反射对，反射对可看做是偶对(r
e
)和奇对(r
o
)的叠加；假设采样率为Δt，两个反射子之间的时间厚度为nΔt，则偶对(r
e
)和奇对(r
o
)分别被表示为：r
e
(t)＝δ(t)+δ(t+nΔt)考虑到整个地震道...
【专利技术属性】
技术研发人员：张帅，肖睿，王若，潘卫红，庞守吉，陆程，祝有海，
申请(专利权)人：中国地质调查局油气资源调查中心，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。