[发明专利]一种直接生成共转换点道集的方法

摘要

本发明提供一种直接生成共转换点道集的方法，假设转换波地震数据已预处理好，已获得准确的动校正参数和计算转换点位置的速度等参数，按炮检距分组，分时窗叠加获得时间和空间变化的共转换点道集。本发明实现简便、适应大数据量处理。

主权项

1.一种直接生成共转换点道集的方法，具体步骤包括：1）采集原始转换波地震数据，经预处理和速度分析，得到转换波动校正和共转换点位置计算需要的参数，并对预处理后转换波地震数据进行动校正，且通过切除消除动校拉伸的影响，步骤1）所述的预处理包括加载观测系统、静校正、波场分离、叠前去噪、振幅恢复、预置共转换点CCP号；2）步骤1）处理后的转换波地震数据按炮检距大小进行分组，每一组对应一个炮检距范围；3）输入一组炮检距范围内的数据；4）按道头中的共转换点CCP号读取转换波用来计算转换点位置的参数；5）对地震道进行分时窗处理，步骤5）所述的分时窗处理是指将地震道分成多个时窗，时窗长度不变，时窗重叠，重叠时窗的重叠长度为时窗内样点数的百分比perc，perc取值为0-100；6）从第一个时窗开始，计算时窗中点的共转换点CCP号，并由这个共转换点号重新获得该时窗用来计算转换点位置的参数，之后重新计算时窗中点相对准确的共转换点CCP号，直接把时窗内的样点值与窗函数相乘叠在共转换点地震道上，即转换波共反射面元的宽度上，相应共转换点地震道每个样点的覆盖次数加上窗函数的数值，步骤6）所述的计算转换点位置的参数是预先速度分析得到的，是转换波均方根速度纵波均方根速度各向异性参数χ_eff、纵横波垂直速度比γ₀和纵横波有效速度比r_eff，步骤6）所述的共转换点CCP号，是指按工区面元划分的，与常规处理中纵波共中心点CMP号一致，步骤6）所述的共转换点CCP号是由下述方法计算的：x_C＝x_M+D式中x_M、x_C为共中心点CMP和共转换点CCP点的水平坐标，D为转换点偏离中心点的水平距离，D=xp-x2]]>式中x_p为共转换点与炮点之间的水平距离，已知共中心点CMP号计算相应的共转换点CCP号jxC=INTn(xCΔxm)=INTn(xM+DΔxm)=jxM+INTn(DΔxm)]]>其中：INT_n表示取最接近的整数，取整隐含做面元化处理,Δx_m为工区面元大小，所述的共转换点与炮点之间的水平距离，采用单层各向同性介质迭代解、解析解、多层介质指定旅行时各向同性近似解、多层VTI介质指定旅行时近似解、倾斜界面迭代解、解析解计算得到，步骤6）所述的直接把时窗内的样点值与窗函数相乘叠在共转换点地震道上，即转换波共反射面元的宽度上，是把这个时窗内的样点叠加在准确的共转换点CCP号上以及相邻的I个共转换点CCP号上，I=0,1,2,…,INT(2·m/Δx+0.5)，m为转换波共反射面元的宽度，INT表示取整，步骤6）所述的转换波共反射面元的宽度采用以下关系式确定：m=1Δm=0[1/(1-Δm)+0.5]Δm≥0.5[1/Δm+0.5]Δm<0.5]]>Δm=(1+r)·Δsr·Δx-INT[(1+r)·Δsr·Δx]]]>其中：Δx为道间距，Δs为炮间距，r为纵横波速度比，INT表示取整；7）重复步骤6）直到该道的每个时窗完成；8）重复步骤4)至7)直到炮检距分组内所有道输入完毕，输出该组的共转换点地震道，步骤8）所述的输出该组的共转换点地震道为共转换点地震道以及共转换点地震道的覆盖次数，并在道头中置上该组炮检距的大小；9）重复步骤3)至8)直到所有的炮检距分组输入完毕，就得到了时间和空间变化的共转换点叠前道集。