一种地震波初至拾取的方法
【专利摘要】本发明公开了一种地震波初至拾取的方法，根据地震剖面选取直达波，作为参考地震子波；参考地震子波与地震数据进行互信息量的计算；在整个时间轴上搜索互信息量极大值，此极大值所对应的时间即为所求的此道地震波初至。该方法包括根据地震剖面选取直达波，作为参考地震子波；地震子波与地震数据进行互信息量的计算；搜索得到互信息量的极大值，此极大值所对应的时间即为所求的地震波初至。本发明的地震波初至拾取的方法所得结果精度高，在低信噪比地区也能取得很好的结果。
【专利说明】一种地震波初至拾取的方法
【技术领域】
[0001] 本发明属于石油勘探领域，尤其涉及一种地震波初至拾取的方法。
【背景技术】
[0002] 石油勘探中地震资料的静校正处理是准确获取地下构造信息的重要处理环节，而 初至拾取则是静校正处理中一项重要的支撑技术。近年来，石油地震勘探工作逐步转向近 地表地震地质条件复杂的地区。在山区，地形复杂，沟梁纵横交错，高程强烈变化；表层岩性 变化剧烈，低速带岩石构造复杂，地震波速度纵向与横向变化快；在黄土覆盖区，巨厚黄土 层影响原始地震资料的采集质量。同时，干燥松散黄土层与质密黄土层地震波速度差别较 大，激发与接收条件难以一致。上述特征给在山区、黄土覆盖区中地震数据的初至拾取带来 极大的困难。
[0003] 现有的初至拾取方法大体可以分为以下几种：
[0004] 1.基于地震记录的瞬时特征属性的方法，主要包括瞬时振幅、瞬时频率、瞬时相位 等特征属性，这类方法主要有极值法、能量比值法，这类方法是对地震数据进行直接计算， 使用起来简单、方便，并且不容易出错，但是这类方法也存在其劣势就是对噪声比较敏感， 在处理信噪比高的资料时效果较好，但对于信噪比低的数据，则难以提取初至波。
[0005] 2.基于地震记录中各地震道整体特征的方法，主要是相关法、约束初至波提取法 等，这类方法在一定程度上可以有效的抑制噪声，但比较容易受地震道相关性等因素影响， 但从地表情况非常复杂的地质环境中采集来的地震资料，用这类方法处理效果比较差，且 受初至波之后的续至波影响大，与提取的结果与子波的选取关系密切。
[0006] 3.基于人工神经网络初至波提取方法，这类方法对于抑制噪声起到很大作用，且 能够将噪声干扰下的有效信号提取出来，但这种方法前期需要进行大量的处理分析，系统 自学习速度慢，并且算法实现起来比较复杂、执行效率很低。
[0007] 4.基于分形维数的方法，这类方法对于处理复杂的地震资料效果极好，并且抗噪 能力强，提取的初至波精度较高，但这类方法的精度与运算速度与多维数的选取关系很大， 计算量非常大。
[0008] 5.基于图形图像学的方法，主要是边缘检测法、边界追踪法，这类方法对于初至波 附近的噪声非常敏感，当初至波附近的噪声非常强时，算法效果就会变差。
【发明内容】
[0009] 本发明的目的在于提供一种地震波初至拾取的方法，旨在解决现有的初至拾取方 法存在的结果精度低，在低信噪比地区不能取得很好的结果的问题。
[0010] 本发明是这样实现的，一种地震波初至拾取的方法包括：
[0011] 步骤一、根据地震剖面选取直达波，作为参考地震子波；
[0012] 步骤二、地震子波与地震数据进行互信息量的计算；
[0013] 步骤三、搜索得到互信息量的极大值，此极大值所对应的时间即为所求的地震波 初至。
[0014] 进一步，参考地震子波与地震数据进行互信息量计算的具体过程如下：
[0015] 地震剖面中待处理的每一道地震记录作为一个时间序列，简称为处理序列，参考 地震子波作为另外一个时间序列，简称为参考序列；
[0016] 步骤一、定义参考序列为X (η)，处理序列为y (η)，n = 0,1，......L-1，把序列x的 值域平均等分为M个区间，M任意确定，统计落入每个小区间内的样本点数，计算序列x (η) 的概率分布P(x);
[0017] 步骤二、按照步骤一的方法计算序列y (η)的概率分布ρ (y);
[0018] 步骤三、将序列X(η)和序列y (η)张成二维空间，这样序列X(η)和y (η)在此张 成的二维空间等分为M2个小区间，统计落入这些小区间的样本数，利用这些样本点数，计算 p(x，y);
[0019] 步骤四、根据公式
【权利要求】
1. 一种地震波初至拾取的方法，其特征在于，该地震波初至拾取方法包括： 步骤一、根据地震剖面选取直达波，作为参考地震子波； 步骤二、参考地震子波与地震数据进行互信息量的计算； 步骤三、在整个时间轴上搜索互信息量极大值，此极大值所对应的时间即为所求的此 道地震波初至。
2. 如权利要求1所述的地震波初至拾取的方法，其特征在于，参考地震子波与地震数 据进行互信息量计算的具体过程如下： 地震剖面中待处理的每一道地震记录作为一个时间序列，简称为处理序列，参考地震 子波作为另外一个时间序列，简称为参考序列； 步骤一、定义参考序列为X(η)，处理序列为y(η)，n= 0,1，......L_l，把序列x的值域 平均等分为M个区间，M任意确定，统计落入每个小区间内的样本点数，计算序列x(η)的概 率分布P(X); 步骤二、按照步骤一的方法计算序列y(η)的概率分布P(y); 步骤三、将序列X(η)和序列y(η)张成二维空间，这样序列X(η)和y(η)在此张成的二 维空间等分为M2个小区间，统计落入小区间的样本数，利用样本点数，计算ρ(X，y); 步骤四、根据公式//<^ = ^^^ = _1>^〇1(^^〇计算信源熵11(幻和!^)，利用 公式所奸）=五[八砂)]=-ΣΣ夕(戈八）1〇gAu.)计算联合熵H(XY); I J 步骤五、根据所求的H(X)，H(Y)和H(XY)，利用公式I(X;Y) =H(X)+H(Y)-H(XY)计算 互信息量。
【文档编号】G06F19/00GK104459789SQ201410735186
【公开日】2015年3月25日 申请日期:2014年12月8日 优先权日:2014年12月8日
【发明者】翟明岳 申请人:翟明岳