一种各向异性速度分析和动校正方法

摘要

本发明涉及一种高精度各向异性速度分析和动校正方法；地震数据准备，分选好CMP地震数据道集；采用双曲方程对小炮检距反射波数据进行无拉伸影响的速度分析；通过对速度谱进交互解释，获得反射波双曲速度；根据反射波双程时间和双曲速度，采用非剥层优化方法计算层速度；根据最小二乘优化原理计算出唯一的层速度；计算垂直均方根速度；利用全炮检距数据进行无拉伸影响的各向异性速度分析，得到各向异性参数；计算更高精度垂直均方根速度和各向异性参数；进行各向异性动校正；利用已经计算出的双曲速度和炮检距范围，自动完成垂直均方根速度的计算，进而进行各向异性速度分析、动校正，计算获得的参数不受动校正拉伸影响，具有高精度、高效率特点。

主权项

一种高精度各向异性速度分析和动校正方法，其特征在于：(1)地震数据准备，包括原始地震数据读入、地震数据编辑、观测系统加载和预处理等；分选好CMP地震数据道集；(2)给定小炮检距范围、速度分析窗长、步长，采用双曲方程对小炮检距反射波数据进行无拉伸影响的速度分析；采用无拉伸影响速度分析方法进行双曲速度分析，通过对速度谱进交互解释，获得反射波双曲速度；(3)根据反射波双程时间和双曲速度，采用非剥层优化方法计算层速度；根据双程时间、叠加速度和炮检距范围，重构反射波旅行时间随炮检距的变化，进而计算反射波旅行时间及其斜率随炮检距的变化，利用某层顶界面、底界面反射波时间及其斜率随炮检距的变化，采用下述公式计算该层的层速度 v n = Δx Δt dx dt 其中，x为炮检距，t为反射波旅行时间，dx/dt为反射波旅行时间的斜率，Δx为相同斜率dx/dt对应的该层顶、底反射的炮检距差，Δt为相同斜率dx/dt对应的该层顶、底反射时间差；由于斜率随炮检距的变化而变化，因此不同的斜率值均能计算层速度，根据最小二乘优化原理计算出唯一的层速度。(4)计算垂直均方根速度；(5)利用全炮检距数据进行无拉伸影响的各向异性速度分析，得到各向异性参数；在进行双曲速度分析、各向异性速度分析时，沿着双曲速度、各向异性参数所确定的反射波时间随炮检距变化轨迹，直接从原始地震数据中取出数据，形成矩阵数据；取出的数据的宽度为道数，长度为子波主周期的0.5至1.5倍之间；通过对此矩阵数据进行叠加判别参数是否合适，由于不受动校正拉伸的影响，因此参数估计的精度得到提高；(6)根据计算出的双曲速度、各向异性参数，采用非剥层优化方法计算更高精度垂直均方根速度和各向异性参数；(a)根据计算得到的反射波双程旅行时间、双曲速度、各向异性参数，利用各向异性时距方程，计算出全炮检距反射波旅行时间及其斜率随炮检距的变换；(b)采用非剥层优化层速度计算方法计算层速度；(c)计算垂直均方根速度；(d)根据垂直均方根速度、双程旅行时间、反射波旅行时间随炮检距的变化，利用各向异性时距方程，采用最小二乘优化的方法计算各向异性参数。(7)进行各向异性动校正。