[发明专利]一种初至波迭代拾取的方法

摘要

本发明是石油物探资料处理的初至波迭代拾取的方法，计算每一道每个采样点处的能量比值，得到能量比值曲线，确定初至预测值，以此为起点找到其中最大振幅值，计算相邻道的初至时间差的绝对值和的平均值，取一道和前后两道初至时间，排列计算可信度系数，对可信度系数为0的初至利用前后相邻4个可靠初至进行拉格朗日多项式拟合插值，得到初至预测值。本发明对于初至波形变化较大、各种波相互干扰，同时外界干扰等的影响的复杂地区，可以获得准确的初至时间，解决了初至拾取方法的难题。

主权项

1.一种初至波迭代拾取的方法，特点是通过以下技术方案实现：1)激发并采集记录单炮地震数据；2)利用以下公式计算每一道每个采样点处的能量比值：R(r)=(NM+A2W)4×(M-N)2r2---(2)]]>M=Σp=1rxp2]]>N=Σp=rr+Wxp2]]>A=1nΣr=1n|xr|]]>式中：每一道的第一个采样点到当前采样点的长度为长时窗，r为长时窗的结束样点，n为一道的采样点数，N为当前点之后W个采样点的能量之和，W为短时窗的计算点数，A为所用道的所有点振幅绝对值和的平均值，M为所用道的长时窗内所有采样点的能量之和；3)取任一道，利用步骤2)中公式(2)计算该道每个采样点处的能量比值，得到能量比值曲线，在能量比值曲线上确定门槛值对应的时间点，在对应的原始地震道上该时间点为初至预测值；4)以初至预测值为起点，长度缺省30个采样点为第一时窗，找到其中最大振幅值，修正波峰点的能量参数与最大振幅的乘积，作为振幅门槛值，以初至预测值为起点，长度是第一时窗的1/2为第二时窗，第二时窗内第一个振幅大于振幅门槛值的位置，就是所要找的初至波峰点，如果在第二个时窗内找不到满足条件的振幅大于振幅门槛值的位置，则该道初至不存在；5)计算相邻道的初至时间差的绝对值和的平均值式中：t₁，t₂，...，t_m是第一道到最后一道的初至时间，m为一炮的总道数；以五道为一组，任一道初至时间t_i，i为道号，i大于等于1小于等于m，满足下面条件中一个，该道初至的可信度系数为1，如果均不满足，该道初至的可信度系数为0：第一个条件：j＝i-4，i-3，i-2，i-1；第二个条件：j＝i-3，i-2，i-1，i+1；第三个条件：j＝i-2，i-1，i+1，i+2；第四个条件：j＝i-1，i+1，i+2，i+3；第五个条件：j＝i+1，i+2，i+3，i+4；t_j为t_i相邻道的初至时间，j为道号，j大于等于1小于等于m，门槛值k为大于0小于等于2，缺省情况下为1。6)每炮除了前两道和后两道，任取一道初至时间t_i，对该道和相邻的前后四道初至时间t_i-2，t_i-1，t_i，t_i+1，t_i+2排序，如果该道排序后的下标与排序前下标的差的绝对值等于2，该道初至的可信度系数为0，否则该道初至的可信度系数为1；7)通过步骤5)和步骤6)，每一炮除了前两道和最后两道外每道初至都得到两个可信度系数，对于任一道，如果得到的两个可信度系数相同，该道初至可信度系数不变，如果两次得到的可信度系数不同，则该道初至的可信度系数为0，对于每炮的前两道和后面两道取步骤5)得到可信度系数，每一道可信度系数等于1的初至为可靠初至；8)对可信度系数为0的初至利用前后相邻4个可靠初至进行拉格朗日多项式拟合插值，得到初至预测值；对于边缘的可信度系数为0的初至，则利用二次曲线拟合插值，得到边缘的初至预测值；初至预测值向上找最近的波峰点，向下找第一个波峰点，如果向上最近波峰点的振幅值小于向下第一个波峰点的振幅值和步骤4)中修正波峰点的能量参数的乘积，初至预测值下方的第一个波峰点为初至波峰点，否则初至预测值向上最近的波峰点为初至波峰点。