[发明专利]一种基于纵横波匹配的多分量资料联合速度分析方法

摘要

本发明提供了一种基于纵横波匹配的多分量资料联合速度分析方法，属于地球物理勘探数据处理领域。所述方法利用纵横波场统一的反射波时距方程同步在纵波和转换波波场进行多参数速度分析，利用纵波和转换波道集的同步对比，实现纵横波相互间速度趋势的同步对照；利用标志层位，同步估算相对准确的纵横波垂直速度比，并通过剖面压缩，利用深度一致原则，循环速度迭代，求取准确的纵横波速度场。利用本发明，可以准确地对转换波地震数据进行处理，获得了高分辨率、高保真、高信噪比以及能比较清楚地反映出构造和岩性变化特征的地震剖面。

主权项

一种基于纵横波匹配的多分量资料联合速度分析方法，其特征在于：该方法利用纵横波场统一的反射波时距方程同步在纵波和转换横波波场进行多参数速度分析，利用纵波和转换横波道集的同步对比，实现纵横波相互间速度趋势的同步对照；利用标志层位，同步估算相对准确的纵横波垂直速度比，并通过剖面压缩，利用纵横波同一反射界面深度一致原则，循环速度迭代，求取准确的纵横波速度场；所述方法具体包括以下步骤：(1)纵、横波场动态共反射点道集提取步骤一：提取纵波CDP道集数据与转换横波ACP道集数据；步骤二：估计工区内纵横波速度比的最大值γmax与处理过程中所需要考虑的最小深度hmin，计算该深度与速度比下各地震道对应的转换点位置xp1与yp1；步骤三：估计工区内纵横波速度比的最小值γmin并据此计算各地震道对应的渐近转换点位置xp2与yp2，渐近转换点位置采用式(1)计算；xp2=xs+(xr-xs)11+γminyp2=ys+(ys-ys)11+γmin---(1)]]>其中xs、ys、xr、yr分别为炮点与检波点的纵横坐标；经过上述步骤二和步骤三后各地震道所对应的转换点范围已确定，即当扫描速度比γ满足条件γmin≤γ≤γmax，且深度h满足条件h≥hmin时，各地震道各采样点对应的转换点位置(xp,yp)必满足条件：或步骤四：对于给定的ACP道集上的点，根据各道的转换点范围抽取所有转换点位置可能落在该点的地震道，形成超道集；(2)纵横波联合速度建模步骤一：建立纵横波场统一的反射波时距方程，公式如下：tref2=tref02+x2Vref2-2ηx4Vref2[tref02Vref2+μx2]]]>上式中，tref纵横波反射波场反射波到达时间，tref0纵横波反射波场双程反射时间，Vref纵横波反射波场叠加速度，x为激发点到接收点之间的水平距离；μ、η各向异性系数；对于P‑P反射波情况下，两个各向异性系数有线性关系：μ＝1+2η；对于P‑S反射波情况下，两个各向异性系数可分别表示如下：η=(γ0γeff-1)2+8χeff(1+γ0)8γ0(1+γeff)2]]>μ=2η(1+γ0)[(γ0-1)γeff2+2χeff](γ0-1)γeff(γ0γeff-1)+2(1+γ0)χeff]]>其中γ0是垂向速度比，γeff是等效速度比，χeff是转换横波各向异性系数；(3)利用统一的时距方程实现纵横波联合速度分析步骤一：应用统一的时距方程，利用动态抽取的CDP道集与ACP道集，求取纵横波反射波场叠加速度；步骤二：利用纵波和转换横波道集的同步对比，实现纵横波相互间速度趋势的同步对照，并进行动校正；步骤三：利用标志层位的纵横波速度场，同步估算相对较准确的初始纵横波垂直速度比；步骤四：利用求取的纵横波垂直速度比，对纵波动校道集进行拉伸或者对转换横波动校道集进行压缩，并对拉伸或者压缩后的纵波和转换横波动校道集进行标志层位的对比，利用深度一致原则，当标志层位处于同一位置时说明速度比正确，若不处于同一位置，说明速度比不准确，也就是速度不准确，需要重新返回，从步骤一开始进行循环迭代速度分析，所述步骤一是：应用统一的时距方程，利用动态抽取的CDP道集与ACP道集，求取纵横波反射波场叠加速度；步骤五：通过拉平的动校道集、纵横波时间域的压缩对比，通过多次迭代分析，得到准确的纵横波速度比和速度场，为叠加以及偏移提供了可靠的保障；步骤六：对上述拉平的动校道集进行精细的时空变动校切除及叠加，获得多分量地震资料的叠加剖面；步骤七：准确成像：以步骤五得到的叠加速度场作为初始输入，进行偏移速度分析，得到准确的偏移速度场，通过叠前时间偏移得到多分量资料的叠前时间偏移剖面。