[发明专利]基于频谱特征的转换横波处理方法

说明书

技术领域

本发明总体来说涉及多波地震勘探技术，更具体地讲，涉及一种基于频谱特征的转换横波高分辨率处理方法。

背景技术

转换横波由于其动力学特性，在频率方面与纵波存在较大的差异：转换横波的主频低、频带范围窄(通常情况下主频为16Hz-20Hz，频带范围为10Hz-40Hz左右)，而纵波的主频较高、频带范围较宽(通常情况下主频能达到30Hz-40Hz、频带范围为5Hz-70Hz，甚至更宽)。并且转换横波在传播过程中，能量、频率衰减较大，造成转换横波成像的纵向分辨率较低，且在频率特征上表现为17Hz-18Hz的强能量频率特征，无法与纵波进行较好地匹配，不能满足多波地震资料联合解释的需求。

为了使转换横波和纵波进行较好地联合解释，通常是以纵波为标准，对转换横波进行高分辨处理。目前针对于低频率、低信噪比、低能量的转换横波，常用的高分辨处理方法主要为反褶积方法，但是无论是叠前反褶积方法还是叠后反褶积方法均是采用提取子波、压缩子波以达到提高分辨率的目的。但是对转换横波执行上述反褶积方法之后，容易造成高频噪声能量增强，信噪比较大程度降低，特别是不能有效地压制17Hz-18Hz频率的强能量以及拓宽频带范围。

发明内容

本发明的目的是提出一种基于频谱特征的转换横波处理方法，以解决转换横波的主频低、频带范围窄无法与纵波进行匹配的技术问题。

本发明的一方面提供一种基于频谱特征的转换横波处理方法，所述方法包括：(a)对时间域的转换横波数据进行变换，得到频率域的转换横波数据；(b)对频率域的转换横波数据的频谱范围内的频率点设定对应的权系数；(c)基于设定的权系数对所述频率点对应的频率成分能量进行调整；(d)对调整后的频率域的转换横波数据进行变换，得到时间域的调整后的转换横波数据；(e)判断调整后的转换横波数据的特征与时间域的纵波数据的特征是否匹配；(f)如果调整后的转换横波数据的特征与时间域的纵波数据的特征不匹配，则执行步骤(g)：基于调整后的转换横波数据的特征与时间域的纵波数据的特征之间的误差确定出导致所述误差产生的频率点，并更新所述频率点的权系数，然后返回执行步骤(c)；(h)如果调整后的转换横波数据的特征与时间域的纵波数据的特征匹配，则执行步骤(i)：输出时间域的调整后的转换横波数据。

可选地，步骤(b)可包括：对频率域的转换横波数据的频谱范围内包含的全部频率点分别设定对应的权系数。

可选地，步骤(b)可包括：对从频率域的转换横波数据的频谱范围内选中的多个频率点设定对应的权系数，然后基于所述多个频率点的权系数利用插值法得到频谱范围内除所述多个频率点之外的其他频率点的权系数。

可选地，步骤(b)可包括：(b1)从频率域的转换横波数据的频谱范围内选中多个频率点；(b2)对所述多个频率点设定对应的权系数；(b3)将所述多个频率点作为区间的端点，将所述频谱范围划分为多个区间；(b4)基于各区间的端点对应的权系数，利用插值法得到各个区间内的除区间的端点之外的其他频率点的权系数。

可选地，步骤(b1)可包括：从频率域的转换横波数据的频谱范围内等频率间隔选中多个频率点，或者从频率域的转换横波数据的频谱范围内任意选中多个频率点。

可选地，所述方法可还包括：获取与时间域的转换横波数据处于相同的目的层的时间域的纵波数据，并对时间域的纵波数据进行变换，得到频率域的纵波数据，其中，对频率域的转换横波数据的频谱范围内的一个频率点设定对应的权系数的步骤可包括：(b11)将频率域的转换横波数据在所述一个频率点对应的频率成分能量与频率域的纵波数据在所述一个频率点对应的频率成分能量进行比较；(b22)如果频率域的转换横波数据在所述一个频率点对应的频率成分能量小于频率域的纵波数据在所述一个频率点对应的频率成分能量，则将所述一个频率点的权系数设置为大于1；(b33)如果频率域的转换横波数据在所述一个频率点对应的频率成分能量等于频率域的纵波数据在所述一个频率点对应的频率成分能量，则将所述一个频率点的权系数设置为等于1；(b44)如果频率域的转换横波数据在所述一个频率点对应的频率成分能量大于频率域的纵波数据在所述一个频率点对应的频率成分能量，则将所述一个频率点的权系数设置为小于1。

可选地，步骤(e)可包括：(e1)将调整后的转换横波数据的特征与时间域的纵波数据的特征进行相似度度量；(e2)根据所述相似度的大小判断调整后的转换横波数据的特征与时间域的纵波数据的特征是否匹配。