[发明专利]一种初至波迭代拾取的方法

说明书

技术领域

本发明涉及石油地球物理勘探资料处理技术，是一种初至波迭代拾取的方法。

背景技术

石油地球物理勘探资料处理中的初至拾取，早期通过手工拾取各道的初至时间，这种方法工作量往往非常庞大，道数甚至超过一万道，这样全部用手工拾取就成为一项费时、费工、十分繁杂的工作。随着地震勘探的不断深入和数据处理精度的不断提高，准确快速地拾取地震波初至是非常重要的。在折射波静校正、VSP解释、浅层折射波勘探和层析成像勘探中，常用初至波的旅行时来反演表层速度和厚度等介质参数，特别是在一些表层较为复杂的地区，如黄土塬、沙漠等，初至拾取是各种初至波静校正的基础。早期的初至拾取方法已不能满足需求。

近年来，公开了不同的初至波自动拾取方法。早期Peraldi和Clement公开了拐点校正法，将地震道与参考道进行互相关，归一化的互相关函数极大值对应的时移即初至时间；1980年Hatherrly公开了线性最小平方预测技术和拐点校正相结合的方法，提出首先识别第一个峰值和拐点，然后估算二者的统计差值；1983年Gelchinsky和Shtivelman公开了一种相关法，假设各道的脉冲形状不发生变化，然而实际情况与假设不符；1987年黄成之等采用统计方法将地震初至波记录分成信号和噪音两个部分，并使这两部分统计特征之间的差别为很大；1990年陈遵德等采用模糊检测方法确定初至波的隶属函数进行初至波检测；1991年Lin Zhang采用相似度准则，建立地震波同相轴的目标函数，用最短路径算法求取目标函数极小，达到自动拾取地震剖面同相轴的目的，但此方法建立目标函数的工作量较大；1994年庄东海等采用将地震记录初至拾取看做一个模式识别过程，充分利用地震记录的瞬时(局部)特征和整体特征，用神经网络方法进行地震记录初至拾取；1995年江玉乐等提出同极性能量比值法，即改进的能量比值法；1996年Fabio Boschetti等提出了一种基于分形维的初至检测算法，该方法是基于在地震道随着信号的出现其分形维就发生变化的特征来确定地震道初至；新近发展起来的BP神经网络法利用多参数特征进行模式识别，充分利用地震记录的瞬时(局部)特征和整体特征来确定地震道初至。

2006年潘树林等利用单道边界检测公式和样条插值进行拾取初至，即：

R(r)=|NM(N-M)|---(1)]]>

M=Σp=rr-nxp]]>N=Σp=rr+nxp]]>

式中：R(r)为第r个样点的边界特征，R(r)越大，代表r点的边界特征越强，最大边界特征值所对应的点就是能量边界。该能量比值公式对初至波到达前的干扰波能量变化过于敏感，对于较差的资料初至拾取误差较大。