[发明专利]一种射线参数域纵波与转换波匹配方法

权利要求书

1.一种射线参数域纵波与转换波匹配方法，其特征在于，所述的方法包括：

对由地震采集得到并经过偏移成像的纵波与转换波CIP道集进行射线追踪，从偏移距域转换到射线参数域；

在所述射线参数域的纵波与转换波CIP道集上标定和拾取对应同相轴的双程旅行时t和τ，并根据t和τ计算出初始的纵横波速度比模型γ(t)；

对每一所述纵波与转换波CIP道集叠加形成纵波和转换波的地震道，在所述纵横波速度比模型γ(t)的基础上加入预定范围的扰动量形成一系列新的压缩比函数，利用所述新的压缩比函数将转换波地震道压缩到纵波时间，并将压缩后的地震道振幅进行重构；

由压缩及重构的所述转换波及纵波地震道得到能量谱，根据所述能量谱的极值点位置修正所述的纵横波速度比模型γ(t),然后根据修正后的所述的纵横波速度比模型将转换波数据匹配到纵波时间得到匹配数据；

利用Gabor变换根据所述匹配数据产生转换波的时频谱，对每个Gabor变换时窗内的频谱一一进行修正，并利用反Gabor变换生成在时间t的转换波地震数据。

2.根据权利要求1所述的射线参数域纵波与转换波匹配方法，其特征在于，对于地震反射能量较弱、相关谱结果无法清晰指示的目的层区域，在得到所述匹配数据的基础上对局部速度比函数的梯度进行扫描。

3.根据权利要求2所述的射线参数域纵波与转换波匹配方法，其特征在于，进行扫描的局部速度比函数的梯度范围是-1.5到1.5，速度比扰动量从-0.2到0.2。

4.根据权利要求2或3所述的射线参数域纵波与转换波匹配方法，其特征在于，所述的射线追踪为弯曲射线追踪，纵波及转换波的双程旅行时T_pp及T_PS分别为：

Tpp=Σi=0k(up(i-1)+up(i))dli]]>

TPS=12[Σi=0k(up(i)+up(i+1))dli+Σi=k2k(us(i)+us(i+1))dli]]]>

其中dl_i是在水平坐标x_i和x_i-1之间的射线路径，u_p(i)和u_s(i)分别是第i层介质的纵波和横波慢度：其中u_p(i)＝u_p(2k-i)u_s(i)＝u_s(2k-i)，i=1,Λ,2k-1。