多次波压制技术与应用(1-2).ppt

基于波动方程的多次波压制法通过波动方程模拟实际波场或反演地震数据来预测多次波，然后把多次波从地震数据中减去; 目前主要有三种方法，波场外推法、自由表面多次波压制法和反散射级数法， 波动外推法是通过波场外推来模拟波在水层中的传播，是模型驱动的， 而SRME和反散射级数法是数据驱动的。 * * * * * * * Weglein将多次波压制技术分为两类，第一类是基于信号分析的滤波方法，第二类是基于波动方程的预测减去法； 滤波方法是主要利用一次波和多次波之间在波传播运动学方面的差异，通过各种变换手段，将一次波和多次波分离开，进而滤除多次波。 * 预测反褶积的期望输出是输入序列向前时移的形式，即给定输入x(t),需要预测它在某一将来(t+m）的值，这里的m是预测步长。 * 在理想的无噪音条件下，预测反褶积对输出的分辨率可以通过调节预测步长来控制。单位预测步长意味着最高的分辨率，而较大的预测步长意味着较低的分辨率。但实际上，这会受到资料信噪比的影响。 短算子产生小振幅尖脉冲并有相对高频的尾巴，长算子使谱进一步白化，进一步靠拢尖脉冲相应谱。 * 在压制多次波时，可以用输入道的自相关来确定压制多次波的合适预测步长，预测步长要选大于自相关的第一部分，它代表地震子波， 算子长度要选择能包括自相关中第一个孤立能量团。 * 在t-x域多次波周期性随偏移距增大逐渐变差。但是将它变到τ-p域，在每一个P道上的多次波都具有较好的周期性，这样我们可以沿着P轴改变预测步长以达到衰减多次波的目的。 * F-K 域去多次波方法是通过二维 Fourier 变换将地震资料从时间—空间（x,t）域，变换到频率—波数（f,k）域， 一次波和多次波能量被分在两个不同的区域里，对f-k域中与多次波能量相当的区域冲零即可压制多次波。 * 最小二乘方法是计算倾斜叠加、抛物线Radon变换、双曲Radon 变换、多项式Radon 变换的一般方法，可以消除褶积截断效应，提高变换的分辨率(聚焦度)，但分辨率和精度仍达不到要求，常规Radon 变换成像存在“蝴蝶”状弥散现象，为进一步提高稳定性、精度和分辨率，应该采用高精度Radon变换。 * 高精度Radon变换求取的目标函数为： 对M求导可以得到 * 这是一个模型数据，一次波与多次波在远炮检距处的时差只有20ms， * 多次波去除前有效信号几乎被多次波淹没，多次波去除后才露出它的庐山真面目。 * 在利用高精度Radon变换进行多次波压制时,Radon变换的曲率参数的确定和Radon变换域中压制区域的确定对于多次波的压制效果影响很大，不合适的参数将会严重影响多次波的压制效果。而常规的批量模块无法直观地确定这些参数，也无法即时让处理员看到处理效果，从而影响了资料处理的效果和效率。因此在资料处理过程中通过交互的方式来确定这些参数就非常重要。 我们自主开发的高精度Radon变换交互分析软件。该交互分析软件功能丰富，使用方便，1：自动化程度高，通过在输入数据上交互定义抛物线，就可以在Radon变换参数卡上自动得到相应的变换参数）。 * 2：通过在输入数据上定义抛物线，就可以在Radon谱上自动得到该抛物线所对应的曲率的位置 * 3：对于输入数据，可以直接计算并显示Radon谱，通过与原始数据的对比分析，就可以方便地分析出数据中一次波和多次波各自的范围； * 4：根据分析的结果，可以在计算出的Radon谱数据上交互地定义多次波切除范围，然后直接进行交互多次波压制，并直接查看压制后的效果 * 5：既可以对多次波压制前后的数据进行拉帘对比显示，也可以将输入数据、要压制的多次波数据和压制多次波后的数据同时并排显示，可以直观地看出处理效果。 * 2.5 聚束滤波压制多次波 * 采用Lagrange乘子法可以获得滤波器为： 有效信号由下式估计： 2.5 聚束滤波压制多次波 * 技术特点 　　聚束滤波方法是一种包括信号和相关噪音的模型拟合的处理方法，由于变换域里采样和截断问题转化为模型拟合问题，没有因数字化导致的假象。 　　聚束滤波模型能够包括振幅和相位随偏移距的变化，及远偏移距处因切除直达波和折射波造成的记录道损失。这样的设计特别有利于避免畸变。 2.5 聚束滤波压制多次波 * 　　聚束滤波本身包括时间延迟、振幅和相位随偏移距变化的参数模型，但是这些模型主要是由聚束滤波的滑动工作时窗中的主要同相轴控制的。 　　聚束滤波器通过自动修正参数模型来适应数据并采用强制性回归方法克服数据中信号的干涉。因此特别适合需要进行时间延迟、振幅和相位随偏移距变化精细分析的处理。 2.5 聚束滤波压制多次波 * 原始CMP道集 Radon后（Ω） 聚束滤波压制后 2.5 聚束滤波压制多次波 * 原始速度谱 Radon后（Ω） 聚束滤波压制后 2.5