一种地震速度谱自动拾取方法
【技术领域】
[0001] 本发明属于油气及煤层气地震勘探与开发领域，具体涉及一种地震速度谱自动拾 取方法。
【背景技术】
[0002] 地震速度谱（以下简称速度谱）与频谱的概念类似，是把地震波的能量相对于波 速的变化关系的曲线称为速度谱。在地震勘探中，速度谱通常指多次覆盖技术中的叠加速 度谱。
[0003] 在地震资料处理及后期的地震数据解释中将非常广泛地利用到速度数据，而速度 数据的来源就是速度谱的拾取。目前速度谱的拾取都是采用人工交互方式完成的，时间长、 效率低。
[0004] 目前现有的速度谱自动拾取手段有基于蒙特卡洛法的速度谱自动拾取等，但是其 技术复杂，效率低，不能全局解释。
【发明内容】
[0005] 本发明的目的在于解决上述现有技术中存在的难题，提供一种地震速度谱自动拾 取方法，将速度谱剖面数据转换成二维网格数据，并利用一定的技术手段搜索出每行峰值 数据，再设置时窗寻找区域最大值，最终利用约束条件完成速度曲线的连接。
[0006] 本发明是通过以下技术方案实现的：
[0007] -种地震速度谱自动拾取方法，包括：
[0008] (1)创建网格数据：将以标准地震数据格式存储的速度谱数据转换成二维网格数 据；
[0009] (2)设置计算时窗并在计算时窗区域内确定最大值位置：
[0010] 定义计算时窗的宽度w，并在计算时窗的区域内搜寻最大值Vmax，并以第一个计 算时窗中的最大值所在行的位置为速度曲线的起始位置；
[0011] 所述计算时窗的区域为wXk ;
[0012] (3)定义速度搜索阈值vh，将计算时窗向下移动一个单元，然后在计算时窗的区 域内搜寻最大值Vmax，如果找到的最大值Vmax大于阈值vh，则记录下该最大值Vmax的位 置，即其对应的二维数组的下标η和k的值，如果找到的最大值Vmax小于阈值vh，则放弃记 录；
[0013] -个单元是指计算时窗的宽度w ;
[0014] ⑷完成速度曲线连接，形成未插值的速度曲线；
[0015] 重复第三步，直到计算时窗依次移动到网格数据的底部，完成所有区域最大值 Vmax所在位置的记录，并将各个位置依次连接，形成未插值的速度曲线；
[0016] (5)对步骤（4)得到的未插值的速度曲线进行速度曲线插值获得最终的速度曲 线。
[0017] 所述步骤（1)是这样实现的：
[0018] 标准地震数据格式存储的速度谱数据中，S(n)表示样点值，η表示数组序列，将标 准地震数据格式存储的速度谱数据转换成网格数据后，网格数据中用S(n，k)表示样点值， η和k分别表示y轴方向与X轴方向的序列。
[0019] 所述步骤（2)中的计算时窗的宽度选择范围为3~20行。
[0020] 所述步骤（3)中阈值vh的取值为第一个计算时窗中的最大值Vmax的二分之一。
[0021] 所述步骤（5)是采用线性插值的方式对未插值的速度曲线进行插值，即对每两个 相邻的位置采用线性插值公式得到这两个位置之间的曲线，然后依次完成所有相邻位置的 计算即得到最终的速度曲线。
[0022] 与现有技术相比，本发明的有益效果是：利用本发明方法进行的速度谱自动拾取， 实现简单、计算速度快，有利于缩减地震数据处理周期，此发明能够为常规油气及煤层气的 勘探、开发服务。
【附图说明】
[0023] 图1速度谱数据的转换
[0024] 图2计算时窗内获取最大值
[0025] 图3计算时窗的移动
[0026] 图4定位阈值的最大值
[0027] 图5速度曲线连接
[0028] 图6本发明方法的步骤框图。
【具体实施方式】
[0029] 下面结合附图对本发明作进一步详细描述：
[0030] 如图6所示，本发明地震速度谱自动拾取方法主要分为以下五个步骤：
[0031] (1)第一步：创建网格数据：
[0032] 速度谱数据一般都是以标准地震数据格式（SEG-Y格式）存储的，SEG-Y格式数据 的特点是以列的方式记录数据，在本发明中为了便于后续步骤的计算，必须先转换成二维 网格数据。其示意如图1所示，SEG-Y数据中S(n)表示样点值，η表示数组序列，网格数据 中S(n，k)表示样点值，η和k分别表示y轴方向与X轴方向的序列。
[0033] (2)第二步：设置计算时窗并在计算时窗区域内确定最大值位置
[0034] 定义计算时窗（w)宽度（一般3~20行为佳（对于地震数据处理来说，一般会把 时窗定义在3~20之间，一般时窗小计算时间就长但精度高，时窗大计算时间就短但精度 低，从程序界面上进行设置即可）），并在计算时窗的区域内（wXk)搜寻最大值Vmax，并以 第一个时窗中的最大值所在行位置为速度曲线的起始位置。示意如图2所示。
[0035] (3)第三步：移动计算时窗，并重复第二步：
[0036] 定义速度搜索阈值（vh)(-般用第一个时窗中的最大值Vmax的二分之一）。然 后将计算时窗向下移动一个单元（此单元即是计算时窗的宽度w)，见图3。并重复第二步。 找到大于阈值vh的最大值Vmax并记录下位置（即记录其对应的二维数组的下标：n和k 值），如果本时窗区域的最大值小于vh，则放弃记录。见图4所示。
[0037] (4)完成速度曲线连接
[0038] 重复第三步（每次重复是用相同的阈值），直到计算时窗依次移动到网格数据底 部，记录下所有区域最大值Vmax所在位置，并依次连接，所连接的曲线即为未插值的速度 曲线，如图5所示。
[0039] (5)对速度曲线插值
[0040] 以上四步完成了速度曲线的自动拾取，但此曲线只记录了每一次计算时窗中最大 值的位置，对于网格数据还有很多行空缺。本发明采用线性插值的方式对曲线进行插值 (即对每两个相邻点采用下面的公式得到这两点之间的曲线，然后依次完成所有相邻点的 计算即得到最终的速度曲线）。得到最终的速度曲线。
[0041] 线性插值法是指使用连接两个已知量的直线来确定在这两个已知量之间的一个 未知量的值的方法。
[0042] 假设已知坐标（x0, y0)与（xl，yl)，要得到[x0, xl]区间内某一位置X在直线上 的值，则有：
[0044] 假设方程两边的值为a，那么这个值就是插值系数一从xO到X的距离与从xO到 xl距离的比值。由于X值已知，所以可以从公式得到a的值
[0048] 这样，在代数上就可以表示成为：
[0049] y = (1- a ) y〇+ a yl
[0050] 或者，
[0051] y = y0+ a (yl-y〇)
[0052] 这样通过a就可以直接得到y。实际上，即使x不在x0到xl之间并且a也不是 介于〇到1之间，这个公式也是成立的。在这种情况下，这种方法叫作线性外插。
[0053] 上述技术方案只是本发明的一种实施方式，对于本领域内的技术人员而言，在本 发明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本 发明上述【具体实施方式】所描述的方法，因此前面描述的方式只是优选的，而并不具有限制 性的意义。
【主权项】
1. 一种地震速度谱自动拾取方法，其特征在于：所述方法包括： (1) 创建网格数据：将以标准地震数据格式存储的速度谱数据转换成二维网格数据； (2) 设置计算时窗并在计算时窗区域内确定最大值位置： 定义计算时窗的宽度w，并在计算时窗的区域内搜寻最大值Vmax，并以第一个计算时 窗中的最大值所在行的位置为速度曲线的起始位置； 所述计算时窗的区域为wXk; (3) 定义速度搜索阈值vh，将计算时窗向下移动一个单元，然后在计算时窗的区域内 搜寻最大值Vmax，如果找到的最大值Vmax大于阈值vh，则记录下该最大值Vmax的位置，即 其对应的二维数组的下标η和k的值，如果找到的最大值Vmax小于阈值vh，则放弃记录； 一个单元是指计算时窗的宽度w ; (4) 完成速度曲线连接，形成未插值的速度曲线； 重复第三步，直到计算时窗依次移动到网格数据的底部，完成所有区域最大值Vmax所 在位置的记录，并将各个位置依次连接，形成未插值的速度曲线； (5) 对步骤（4)得到的未插值的速度曲线进行速度曲线插值获得最终的速度曲线。2. 根据权利要求1所述的地震速度谱自动拾取方法，其特征在于：所述步骤（1)是这 样实现的： 标准地震数据格式存储的速度谱数据中，S (η)表示样点值，η表示数组序列，将标准地 震数据格式存储的速度谱数据转换成网格数据后，网格数据中用S(n，k)表示样点值，η和 k分别表示y轴方向与X轴方向的序列。3. 根据权利要求1所述的地震速度谱自动拾取方法，其特征在于：所述步骤（2)中的 计算时窗的宽度选择范围为3~20行。4. 根据权利要求1所述的地震速度谱自动拾取方法，其特征在于：所述步骤（3)中阈 值vh的取值为第一个计算时窗中的最大值Vmax的二分之一。5. 根据权利要求1所述的地震速度谱自动拾取方法，其特征在于：所述步骤（5)是采 用线性插值的方式对未插值的速度曲线进行插值，即对每两个相邻的位置采用线性插值公 式得到这两个位置之间的曲线，然后依次完成所有相邻位置的计算即得到最终的速度曲 线。
【专利摘要】本发明提供了一种地震速度谱自动拾取方法，属于油气及煤层气地震勘探与开发领域。本方法包括：(1)创建网格数据：将以标准地震数据格式存储的速度谱数据转换成二维网格数据；(2)设置计算时窗并在计算时窗区域内确定最大值位置；(3)定义速度搜索阈值vh，将计算时窗向下移动一个单元，然后在计算时窗的区域内搜寻最大值Vmax，如果找到的最大值Vmax大于阈值vh，则记录下该最大值Vmax的位置，即其对应的二维数组的下标n和k的值，如果找到的最大值Vmax小于阈值vh，则放弃记录；(4)完成速度曲线连接，形成未插值的速度曲线；(5)对步骤(4)得到的未插值的速度曲线进行速度曲线插值获得最终的速度曲线。
【IPC分类】G01V1/30, G01V1/32
【公开号】CN105572733
【申请号】CN201410524637
【发明人】宋志翔, 杨强, 刘百红, 岳承琪
【申请人】中国石油化工股份有限公司, 中国石油化工股份有限公司石油物探技术研究院
【公开日】2016年5月11日
【申请日】2014年10月8日