针对储层含油气性检测的角度域叠前偏移数据处理方法

摘要

本发明涉及一种针对储层含油气性检测的角度域叠前偏移数据处理方法，包括以下步骤：1)模拟地面地震，数据采集器采集地震波信号，按振幅相对保真处理流程进行预处理；2)输入设备根据储层预测需要选定局部目标成像范围，并对角度域成像数组进行初始化；3)数据采集器将步骤1)预处理后的地震波数据逐道输入到处理器，处理器根据中心点位置与偏移孔径限制，确定该道地震波对应的最大成像范围，并判断该道最大成像范围是否在选定的局部目标成像范围内，若为是，执行步骤4)，若为否，输入下一道等。与现有技术相比，本发明具有根据地震波局部方向特征实现局部角度域叠前时间偏移成像，直接在偏移过程中一步获得高保真的入射角域共成像点道集等优点。

主权项

1.一种针对储层含油气性检测的角度域叠前偏移数据处理方法，其特征在于，包括以下步骤：1）模拟地面地震，数据采集器采集地震波信号，按振幅相对保真处理流程进行预处理；2）输入设备根据储层预测需要选定局部目标成像范围，并对角度域成像数组进行初始化；3）数据采集器将步骤1）预处理后的地震波数据逐道输入到处理器，处理器根据中心点位置与偏移孔径限制，确定该道地震波对应的最大成像范围，并判断该道最大成像范围是否在选定的局部目标成像范围内，若为是，执行步骤4），若为否，输入下一道；4）处理器对该道地震波进行每个成像时间孔径内所有点的角度域成像处理，判断该道地震波所有的成像时间的角度域成像处理是否完成，若为是，执行步骤5），若为否，继续执行步骤4）；5）处理器判断所有的地震波数据是否都理处完成，若为是，执行步骤6），若为否，返回步骤3）；6）得到局部成像空间的角度域共成像点道集和部分角度叠加剖面，基于角度域共成像点道集进行振幅随入射角变化AVA分析或反演，或基于部分角度叠加剖面进行弹性阻抗EI反演，获得反映储层及其含油气性的属性数据与图表，并通过显示器显示；所述的步骤4）中的处理器对该道地震波进行每个成像时间孔径内所有点的角度域成像处理过程如下：（1）从角度域成像空间开始时刻进行偏移脉冲响应计算：根据当前成像时间τ的偏移孔径范围以及中心点坐标，确定当前成像时间成像点横向分布区间；（2）在成像点横向分布区间内，针对每个成像点（x,y,τ），根据非双曲双平方根DSR时差方程计算炮点-成像点单程射线走时t_s以及成像点-接收点单程射线走时t_r，再利用公式：计算入射与散射射线参数p_s与p_r，再根据Snell定律计算射线的起飞角α_s与α_r，二维各向同性介质情况下按γ_2D＝（α_s+α_r）/2计算得到入射角；在三维各向同性介质情况下，通过${\gamma }_{3D}=\frac{1}{2}{\cos }^{-1}(\cos {\alpha }_s\cos {\alpha }_r+\sin {\alpha }_s\sin {\alpha }_r\cos \beta )$计算得到入射角，其中β是入射射线所在垂直平面与散射射线所在垂直平面的夹角；对于二维情况，通过如下公式计算几何扩散校正权系数：$w_{2D}=\frac{\sqrt{\cos {\alpha }_{s_0}\cos {\alpha }_{r_0}}}{v}(\frac{t_r}{t_s}+\frac{t_s}{t_r})$对于三维情况，通过如下公式计算几何扩散校正权系数：$w_{3D}=\frac{\sqrt{\cos {\alpha }_{s_0}\cos {\alpha }_{r_0}}}{v}\sqrt{t_s{t}_r}(\frac{t_r}{t_s}+\frac{t_s}{t_r})(\frac{1}{t_s}+\frac{1}{t_r})$其中与分别对应入射射线与散射射线在地面与垂向的夹角，它们能由各自的射线参数以及地面速度按Snell定律计算得到；3）将输入地震道对应走时的振幅乘以几何扩散校正权系数后叠加到角度域成像空间（x,y,γ,τ），其中x,y,τ为时间域成像坐标，γ为入射角；完成当前成像时间孔径内所有点的角度域成像后，进入下一个成像时间，返回步骤1）。