[发明专利]一种角度域叠前时间偏移方法

摘要

本发明实施例公开了一种角度域叠前时间偏移方法，属于地震勘探中反射地震资料处理技术领域，该方法首先通过偏移速度场和地震数据求取每个成像点处的拟成像地层倾角和偏移幅值，并在此基础上建立二维数组，得到拟成像倾角道集，然后借助该二维数组在每个成像点处求取拟成像地层倾角的两个边界，来确定精确的倾角成像区，在该精确的倾角成像区内通过二维数组对偏移幅值进行叠加，生成地震成像剖面，从而避免了成像过程中掺入对成像没有贡献的噪音信号，提高了成像剖面的信噪比，从而使成像更为准确。

主权项

1.一种角度域叠前时间偏移方法，其特征在于，该方法包括：步骤1：获取叠前地震数据，对所述地震数据进行预处理，得到选定成像线的偏移速度场；步骤2：根据所述偏移速度场和所述地震数据，在所述选定成像线对应的每个成像点处求取拟成像地层倾角和偏移幅值；步骤3：以时间深度和所述拟成像地层倾角为维度定义二维数组，形成时间深度维和拟成像倾角维，所述二维数组在所述时间深度维上包括若干时间采样点，在所述拟成像地层倾角维上包括若干倾角采样点，在对应每个所述时间采样点和所述倾角采样点存放所述偏移幅值，得到所述选定成像线的存放有所述偏移幅值的所述二维数组，所述步骤3包括：定义所述选定成像线上等间距的若干离散点为CDP，在所述选定成像线上的每个CDP处，定义所述二维数组A[m][n]存放所述偏移幅值，所述二维数组的第一维为所述拟成像地层倾角维，所述拟成像地层倾角维具有m个所述倾角采样点，m定义为奇数，定义角度间隔为Δφ，则m可表示为m＝π/Δφ+1，定义正整数P，将表达式m/2+1的数值取整后赋予P，在所述拟成像地层倾角维中，m＝P位置存放0度的所述拟成像地层倾角；所述二维数组的第二维为所述时间深度维，所述时间深度维具有n个所述时间采样点数，定义时间深度为T，定义时间采样间隔为Δt，则n可表示为n＝T/Δt；在每个所述CDP处，将所述成像点的所述偏移幅值根据所述成像点的所述时间深度值和所述拟成像地层倾角对应的所述时间采样点和所述倾角采样点存放到所述二维数组中，得到所述选定成像线的存放有所述偏移幅值的所述二维数组；步骤4：根据存放有所述偏移幅值的所述二维数组，在每个所述时间采样点处求取所述倾角采样点的两个边界，确定所述二维数组在所述拟成像地层倾角维上的叠加区域，所述步骤4包括：对所述二维数组A[m][n]上每个所述时间采样点，令i为所述倾角采样点的索引值，令bi是对应所述倾角采样点的所述偏移幅值；a)求取振幅加和关于x的函数其中x∈[1,P]，x为正整数，f₀为f(x)绝对值的平均值，对f(x)求取x的一次导数f'(x)，当f'(x)的绝对值不大于0.001*f₀时，记录下x的数值，该数值为所述时间采样点处所述倾角采样点的所述两个边界中的一个边界；b)求取振幅加和关于x的函数其中x∈[P,m]，x为正整数，f₀为f(x)绝对值的平均值，对f(x)求取x的一次导数f'(x)，当f'(x)的绝对值不大于0.001*f₀时，记录下x的数值，该数值为所述时间采样点处所述倾角采样点的所述两个边界中的另一个边界；c)对其他所述CDP位置处的所述拟成像倾角道集重复步骤a)与b)，可得到每个所述CDP位置处的所述二维数组的所述倾角采样点的两个边界，确定所述二维数组在所述拟成像地层倾角维上的叠加区域；步骤5：在每个所述时间采样点处，对所述叠加区域内的所述二维数组存放的所述偏移幅值分别进行叠加，得到所述选定成像线的成像数据体；步骤6：根据所述成像数据体显示所述选定成像线的地震成像剖面。