[发明专利]随机函数中值滤波技术方法

说明书

本发明公开了随机函数中值滤波技术方法，包括：对地震记录中要去除的线性干扰进行干扰波拉平；对拉平后的数据利用随机函数进行重新排序，并记录下排序前后地震道的相对关系；对排序后的地震记录进行中值滤波处理；对排序并进行中值滤波处理后的地震道进行反随机排序，恢复其原始的地震道顺序；对反随机排序后的数据进行反干扰波拉平，得到去除线性干扰后的地震记录。本发明方法可实现线性干扰的压制，可以提高分布式光纤传感技术在VSP采集中的应用，提高资料信噪比，满足高精度VSPCDP成像，为井旁精细构造解释提供参考依据。

技术领域

本发明涉及油田开发技术领域，特别是涉及一种随机函数中值滤波技术方法。

背景技术

线性干扰是地震资料采集中普遍存在的一种规律干扰波，影响地震资料叠加成像，降低资料信噪比。它通常是由地表的干扰源发出干扰，由地震检波器信号接收而形成。另外，当地震波经过地表障碍物的时候，形成新的二次源，产生干扰，被地震检波器接收从而形成干扰波。

干扰波在不同的空间域会表现出不同的形式，但干扰波总能在其中的一个域以线性干扰的形式存在，方便预测及压制。线性干扰具有固定的视速度，能量衰减不大，其频率特征也较固定，因此常用的线性干扰去除方法有FK滤波、τ-р域滤波、中值滤波等。

随着光纤技术的发展，分布式光纤声波传感技术被应用到VSP资料采集中，大大提高了采集效率，降低采集成本。但是，利用分布式光纤传感技术进行VSP资料采集时，当井源距较小时，地震波能量传递到铠装光缆，沿着光缆传播，形成能量较强的线性干扰，我们称之为缆波，严重干扰资料信噪比。

并且由于缆波能量强，视速度与部分有效波视速度相当，常规的FK滤波、τ-р域滤波、以及种植滤波压制缆波效果不理想。主要表现在当FK滤波、τ-р域滤波采用的视速度跨度较小时，其缆波能量不能有效压制，当视速度跨度加大时，对有效信号的损伤较大。中值滤波存在类似的问题，其跨度较小时，对有效信号损伤较大，当跨度较大时，对缆波的压制效果不理想。

分布式光纤传感技术采集到的缆波干扰是线性干扰的一种。实现线性干扰的压制，可以提高分布式光纤传感技术在VSP采集中的应用，提高资料信噪比，满足高精度VSPCDP成像，为井旁精细构造解释提供参考依据。

发明内容

本发明的目的在于，针对现有技术的不足，提供随机函数中值滤波技术方法。

本发明采用的技术方案如下。

随机函数中值滤波技术方法，其特征在于，包括下列步骤：

步骤1：对地震记录中要去除的线性干扰进行干扰波拉平。

步骤2：对拉平后的数据利用随机函数进行重新排序，并记录下排序前后地震道的相对关系；

步骤3：对排序后的地震记录进行中值滤波处理；

步骤4：按照步骤2的对应关系，对排序并进行中值滤波处理后的地震道进行反随机排序，恢复其原始的地震道顺序；

步骤5：对反随机排序后的数据进行反干扰波拉平，得到去除线性干扰后的地震记录。

作为优选技术方案，步骤2中的随机函数采用的是线性余同法。

作为优选技术方案，随机函数采用的公式如下：

x_n＝(ax_n-1+b)mod m

r_n＝x_n/m

tr_n＝x_n modΔtr+tr₀；