剩余静校正互相关处理

最新推荐文章于 2024-01-02 17:33:34 发布

帅雪人

最新推荐文章于 2024-01-02 17:33:34 发布

阅读量641

点赞数

文章标签： python 开发语言

本文链接：https://blog.csdn.net/Alva_xtt/article/details/123566193

版权

本文介绍了使用互相关算法进行剩余静校正的过程，包括计算相对静校正量和分离绝对静校正量。通过校正前后的图像对比，分析了存在的问题并提出了后续研究思路。

摘要由CSDN通过智能技术生成

文章目录

前言
一、互相关函数与相对静校正量
二、分离相对静校正量
三、校正前后对比
四、原因分析以及后续思路
- 原因分析
- 后续思路

前言

关于剩余静校正已经有了较深的体会，主要完整实现了互相关算法，包括相对静校正量的计算以及每个点的绝对静校正量的分离。

一、互相关函数与相对静校正量

这是使用互相关进行剩余精校正的第一个步骤，主要利用互相关函数计算出每个地震道与模型道的互相关函数值，根据计算出的曲线，寻找二者之间的相对静校正量。

def cross_debug(para_trace_set):
	"""
    Calculate statics based on real seismic data.
    :param para_trace_set: real seismic data
    :return: relative statics for all traces
    """
    # Obtain model trace
    temp_model_trace = deter_model_trace(para_trace_set)
    # Delete useless information
    para_trace_set = np.delete(para_trace_set, [-1, -2, -3], axis=1)

    # Cross-correlation function value
    temp_cross_value = []
    # Cross-correlation process
    # gamma_t = sum(y_i * x_{i + t})
    for i in range(START, END):
        temp_cross_value_per = np.zeros(para_trace_set