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

摘要

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

主权项

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