[发明专利]地震信号的数据处理方法、装置、电子设备及存储介质

权利要求书

1.一种地震信号的数据处理方法，其特征在于，包括：

利用奇异谱分析方法对原始地震信号数据进行滤波，得到初始滤波数据；

确定所述原始地震信号数据和所述初始滤波数据之间的差异性和相似性；

基于所述差异性、所述相似性、预先设置的加权强度控制参数和局部相似性阈值计算加权矩阵；所述加权强度控制参数用于预测不稳定扰动的能量衰减，所述局部相似性阈值用于确定信号保留程度；所述加权矩阵用于预测不稳定扰动的位置和能级；

基于所述加权矩阵对所述原始地震信号数据再次进行滤波处理；

所述确定所述原始地震信号数据和所述初始滤波数据之间的差异性和相似性，包括：

调用绝对偏差计算关系式计算所述原始地震信号数据和所述初始滤波数据之间的绝对偏差U，U＝[u_i,j]，所述绝对偏差计算关系式为：u_i,j为所述绝对偏差对应矩阵中(i，j)位置处的元素，d_ij为所述原始地震信号数据对应矩阵中(i，j)位置处的元素，为所述初始滤波数据对应矩阵中(i，j)位置处的元素；

调用局部相似性计算关系式计算所述原始地震信号数据和所述初始滤波数据之间的局部相似性V，V＝[v_i,j]，所述局部相似性计算关系式为：

式中，v_i,j为局部相似性对应矩阵中(i，j)位置处的元素，ω为局部平滑窗口函数，τ为当前位置(i，j)沿i移动τ个单位，κ为当前位置(i，j)沿j移动κ个单位，d_τ,κ为原始地震信号数据对应矩阵中(τ，κ)位置处的元素，为初始滤波数据对应矩阵中(τ，κ)位置处的元素；

所述基于所述差异性、所述相似性、预先设置的加权强度控制参数和局部相似性阈值计算加权矩阵，包括：

调用加权元素计算关系式计算所述加权矩阵，所述加权元素计算关系式为：

式中，w_i,j为所述加权矩阵在位置(i,j)处的元素，u_i,j为用于反映所述差异性的绝对偏差对应矩阵在位置(i,j)处的元素，ε为归一化的绝对偏差，v_i,j为用于反映所述相似性的局部相似性对应矩阵在位置(i,j)处的元素，η为所述局部相似性阈值，P为所述加权强度控制参数。

2.根据权利要求1所述的地震信号的数据处理方法，其特征在于，所述基于所述加权矩阵对所述原始地震信号数据再次进行滤波处理，包括：

计算所述加权矩阵对所述原始地震信号数据的哈达玛积，得到一次地震信号修改数据；

沿时间方向，对所述一次地震信号修改数据进行正向一维傅里叶变换，得到频率切片；

对所述频率切片依次执行Hankel矩阵嵌入操作、降秩操作、平均化操作，得到二次地震信号修改数据；

沿时间方向，对所述二次地震信号修改数据进行反向一维傅里叶变换，得到所述原始地震信号数据的数据处理结果。