一种异常初至时间修正方法技术

本发明专利技术公开一种异常初至时间修正方法，包括利用预设的能量比值公式计算每道每一个采样点的能量比值，得到每道的能量比值曲线；将每道的能量比值曲线最大值对应的采样点作为相应地震道的初至时间的样点；利用滑动时窗对初至时间线性拟合，确定可靠初至时间；计算可靠初至时间的正负波形面积比值的均值和初至波峰点能量的均值；计算异常初至道在预设的时窗内每个子波的正负波形面积的比值和每个子波的波峰点能量；确定异常初至道的正确初至时间；滑动着对所有异常初至时间进行计算修正。本发明专利技术可以大大提高拾取初至时间的精准度和道数，为后续利用初至时间的属性解决地震数据的静校正问题提供有效的依据。

【技术实现步骤摘要】
一种异常初至时间修正方法
本专利技术涉及石油地球物理勘探资料处理
，具体是一种异常初至时间修正方法。
技术介绍
随着高密度高效采集技术的应用越来越广泛，勘探的目标区逐渐转向复杂地表条件的地区，低信噪比海量数据的初至拾取工作面临巨大的挑战，传统方法需花费大量的人力和时间进行修正异常初至时间，大量的交互影响拾取精度和效率。因此，初至时间自动拾取过程中修正异常初至时间方法研究对低信噪比大数据量的资料显得非常必要。目前初至时间自动拾取的方法有很多，主要包括拐点校正法、相关法、能量比值法、图像边界检测法、神经网络法、分形维数法等，但是在复杂近地表条件下，由于初至时间能量特征、波形特征、相位特征变化大，又常常受噪声干扰，初至时间自动拾取方法的精度仍然无法满足要求，需要对拾取的初至进行修正，目前的修正技术主要是样条内插技术，该技术仅适用于近地表的厚度、速度空间相对稳定地区，而对于近地表构造复杂和初至波信噪比低的区域将失效。因此，现有技术中亟需一种有效的异常初至时间修正方法，从而为后续利用初至时间的属性解决地震数据的静校正问题提供有效的依据。
技术实现思路
本专利技术主要是克服现有技术中的不足之处，提出一种异常初至时间修正方法，该方法提出新的能量比值公式，提高低信噪比地震资料初至时间拾取精度，根据可靠初至时间的能量和波形面积对异常初至进行精准地修正，提高初至时间自动拾取的精准度和道数，为后续利用初至时间的属性解决地震数据的静校正问题提供有效的依据。本专利技术解决上述技术问题所提供的技术方案是：一种异常初至时间修正方法，包括以下步骤：A、根据采集的地震数据计算每道每一个采样点的能量比值，得到所述道的能量比值曲线；B、从所述道的能量比值曲线中得到最大值对应的采样点，并将采样点作为相应地震道初至时间的样点；C、利用滑动时窗对初至时间线性拟合，根据初至时间与拟合线的距离确定可靠初至时间；D、计算可靠初至时间的正负波形面积、两个面积的比值、比值的均值以及初至波峰点能量的均值；E、计算异常初至道在预设的时窗内每个子波的正负波形面积的比值和每个子波的波峰点能量；F、根据上述得到的数据计算异常初至道在预设的时窗内每个子波的综合评价值，确定异常初至时间的正确初至时间；G、滑动着对所有异常初至时间进行计算修正。进一步的技术方案是，所述步骤A中采集的地震数据利用预设的能量比值公式计算每道每一个采样点的能量比值，其预设的能量比值公式如下：式中：Ri(r)代表第i道第r个采样点的能量比值，w为计算时窗的点数，i为地震道道号，1≤i≤m，r为对应的第i道的采样点点号，1≤r≤n，n为一道的采样点数，α为稳定因子，β为大于0的整数。进一步的技术方案是，所述步骤C的具体过程为：a、以第一道初至时间为起点，滑动时窗长度为L，滑动步长为1，逐道滑动着对拾取的初至时间线性拟合，计算滑动时窗内每道初至时间与拟合直线的时间差，若滑动时窗内某道大于给定的阈值δ，则为异常初至，否则为可靠初至；b、以第一道初至时间为起点，确定第k1道为第一个异常初至时间，逐道滑动着向后寻找，直到第k2道为异常初至，且第k2+1道为可靠初至时，将第k1至k2道之间的初至时间标记为异常初至段，对所有拾取的异常初至时间进行标记。进一步的技术方案是，所述步骤D中可靠初至时间的正负波形面积比值的计算公式如下：式中：为第k1-d道可靠初至时间的正负波形面积比值，和分别为第k1-d道初至时间的正波形面积和负波形面积，且d为大于等于1小于等于l的正整数，l为预设的可靠初至时间的道数。进一步的技术方案是，所述述步骤D中比值的均值的计算公式如下：式中：为与第k1-1道相邻的可靠初至时间的正负波形面积比值的均值，为第k1-d道可靠初至时间的正负波形面积比值，d为大于等于1小于等于l的正整数，l为预设的可靠初至时间的道数。进一步的技术方案是，所述所述述步骤D中初至波峰点能量的均值的计算公式如下：式中：为与第k1-1道相邻的l道可靠初至波峰点能量的均值，为第k1-d道可靠初至波峰点的能量，且d为大于等于1小于等于l的正整数，l为预设的可靠初至时间的道数。进一步的技术方案是，所述步骤E的具体过程为：从确定的可靠初至时间相邻的异常初至时间第k1道开始，逐道计算在预设的时窗内每个子波的正负波形面积比值和子波的波峰点能量为式中:和分别为第k1道预设的时窗内某个子波的正波形面积和负波形面积，给定的时窗为以相邻的可靠初至时间为中心，上下分别取T毫秒为时窗的长度,第k1+1道以第k1道修正后的初至为中心给定时窗，逐道滑动着确定第k1+2至k2+1道的时窗。进一步的技术方案是，所述步骤F中异常初至道在预设的时窗内每个子波的综合评价公式如下：式中：S为子波的综合评价值，最小值对应的子波为第k1道的初至时间，和分别为第k1道预设时窗内某个子波的波峰点能量和正负波形面积比值，和分别为与第k1-1道相邻的l道可靠初至波峰点能量的均值和正负波形面积比值的均值，ε1和ε2为加权系数，ε1+ε2＝1，l为预设的可靠初至时间的道数。进一步的技术方案是，所述步骤F中通过预设时窗计算可靠初至的初至时间与实际的初至时间的差的绝对值来确定异常初至道的正确初至时间。进一步的技术方案是，所述确定异常初至道的正确初至时间的具体过程为：根据上述计算预设时窗内每个子波的综合评价值，确定了第k1道的初至时间，将第k1道作为可靠初至时间，依次计算第k1+1至k2+1道的初至时间，计算第k2+1道的初至时间和实际的初至时间的差的绝对值，若所述差的绝对值小于1，则计算的第k1至k2道的初至时间是正确的，否则确定的第k1至k2道的初至时间是错误的；如果上述计算的第k1至k2道的初至时间是错误的，则以第k2+1道为起始位置，依次计算第k2至k1-1道的初至时间，如果计算的第k1-1道的初至时间和实际的初至时间的差的绝对值小于1，则计算的第k2至k1道的初至时间是正确的，否则确定的第k1至k2道的初至时间是错误的，第k1至k2道的初至时间是0。本专利技术具有益效果：本专利技术可以有效地实现对大数据量、低信噪比地震资料的异常初至时间进行计算修正，准确拾取初至时间，与现有技术相比，可以大大提高拾取初至时间的精准度和道数，为后续利用初至时间的属性解决地震数据的静校正问题提供有效的依据。附图说明图1为本专利技术实施例提供的异常初至时间修正方法的流程图；图2为本专利技术实施例提供的炮集数据的示意图；图3为本专利技术实施例提供的单道原始地震记录的示意图；图4为本专利技术实施例提供的能量比值曲线的示意图；图5为本专利技术实施例提供的单炮初至时间拾取的示意图；图6为本专利技术实施例提供的剔除异常初至时间的示意图；图7为本专利技术实施例提供的对异常初至时间修正的示意图。具体实施方式下面结合实施例和附图对本专利技术做更进一步的说明。本专利技术的一种异常初至时间修正方法，包括以下步骤：步骤101、根据采集的地震数据和预设的能量比值公式计算每道每一个采样点的能量比值，得到所述道的能量比值曲线；其预设的能量比值公式如下：式中：Ri(r)代表第i道第r个采样点的能量比值，w为计算时窗的点数，i为地震道道号，1≤i≤m，r为对应的第i道的采样点点号，1≤r≤n，n为一道的采样点数，α为稳定因子，取值0.1，β为大于0的整数，取值1；其中本文档来自技高网...

【技术保护点】
1.一种异常初至时间修正方法，其特征在于，包括以下步骤：A、根据采集的地震数据计算每道每一个采样点的能量比值，得到所述道的能量比值曲线；B、从所述道的能量比值曲线中得到最大值对应的采样点，并将采样点作为相应地震道初至时间的样点；C、利用滑动时窗对初至时间线性拟合，根据初至时间与拟合线的距离确定可靠初至时间；D、计算可靠初至时间的正负波形面积、两个面积的比值、比值的均值以及初至波峰点能量的均值；E、计算异常初至道在预设的时窗内每个子波的正负波形面积的比值和每个子波的波峰点能量；F、根据上述得到的数据计算异常初至道在预设的时窗内每个子波的综合评价值，确定异常初至时间的正确初至时间；G、滑动着对所有异常初至时间进行计算修正。

【技术特征摘要】
1.一种异常初至时间修正方法，其特征在于，包括以下步骤：A、根据采集的地震数据计算每道每一个采样点的能量比值，得到所述道的能量比值曲线；B、从所述道的能量比值曲线中得到最大值对应的采样点，并将采样点作为相应地震道初至时间的样点；C、利用滑动时窗对初至时间线性拟合，根据初至时间与拟合线的距离确定可靠初至时间；D、计算可靠初至时间的正负波形面积、两个面积的比值、比值的均值以及初至波峰点能量的均值；E、计算异常初至道在预设的时窗内每个子波的正负波形面积的比值和每个子波的波峰点能量；F、根据上述得到的数据计算异常初至道在预设的时窗内每个子波的综合评价值，确定异常初至时间的正确初至时间；G、滑动着对所有异常初至时间进行计算修正。2.根据权利要求1所述的一种异常初至时间修正方法，其特征在于，所述步骤A中采集的地震数据利用预设的能量比值公式计算每道每一个采样点的能量比值，其预设的能量比值公式如下：式中：Ri(r)代表第i道第r个采样点的能量比值，w为计算时窗的点数，i为地震道道号，1≤i≤m，r为对应的第i道的采样点点号，1≤r≤n，n为一道的采样点数，α为稳定因子，β为大于0的整数。3.根据权利要求1所述的一种异常初至时间修正方法，其特征在于，所述步骤C的具体过程为：a、以第一道初至时间为起点，滑动时窗长度为L，滑动步长为1，逐道滑动着对拾取的初至时间线性拟合，计算滑动时窗内每道初至时间与拟合直线的时间差，若滑动时窗内某道大于给定的阈值δ，则为异常初至，否则为可靠初至；b、以第一道初至时间为起点，确定第k1道为第一个异常初至时间，逐道滑动着向后寻找，直到第k2道为异常初至，且第k2+1道为可靠初至时，将第k1至k2道之间的初至时间标记为异常初至段，对所有拾取的异常初至时间进行标记。4.根据权利要求1所述的一种异常初至时间修正方法，其特征在于，所述步骤D中可靠初至时间的正负波形面积比值的计算公式如下：式中：为第k1-d道可靠初至时间的正负波形面积比值，和分别为第k1-d道初至时间的正波形面积和负波形面积，且d为大于等于1小于等于l的正整数，l为预设的可靠初至时间的道数。5.根据权利要求1所述的一种异常初至时间修正方法，其特征在于，所述述步骤D中比值的均值的计算公式如下：式中：为与第k1-1道相邻的可靠初至时间的正负波形面积比值的均值，为第k1-d道可靠初至时间的正负波形面积比值，d为大于等于1小于等于l的正整数，l...

【专利技术属性】
技术研发人员：许银坡，尹成，赵虎，
申请(专利权)人：西南石油大学，
类型：发明
国别省市：四川,51