CN105572733B - 一种地震速度谱自动拾取方法 - Google Patents
一种地震速度谱自动拾取方法 Download PDFInfo
- Publication number
- CN105572733B CN105572733B CN201410524637.0A CN201410524637A CN105572733B CN 105572733 B CN105572733 B CN 105572733B CN 201410524637 A CN201410524637 A CN 201410524637A CN 105572733 B CN105572733 B CN 105572733B
- Authority
- CN
- China
- Prior art keywords
- rate curve
- computation window
- maximum value
- interpolation
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种地震速度谱自动拾取方法,属于油气及煤层气地震勘探与开发领域。本方法包括:(1)创建网格数据:将以标准地震数据格式存储的速度谱数据转换成二维网格数据;(2)设置计算时窗并在计算时窗区域内确定最大值位置;(3)定义速度搜索阈值vh,将计算时窗向下移动一个单元,然后在计算时窗的区域内搜寻最大值Vmax,如果找到的最大值Vmax大于阈值vh,则记录下该最大值Vmax的位置,即其对应的二维数组的下标n和k的值,如果找到的最大值Vmax小于阈值vh,则放弃记录;(4)完成速度曲线连接,形成未插值的速度曲线;(5)对步骤(4)得到的未插值的速度曲线进行速度曲线插值获得最终的速度曲线。
Description
技术领域
本发明属于油气及煤层气地震勘探与开发领域,具体涉及一种地震速度谱自动拾取方法。
背景技术
地震速度谱(以下简称速度谱)与频谱的概念类似,是把地震波的能量相对于波速的变化关系的曲线称为速度谱。在地震勘探中,速度谱通常指多次覆盖技术中的叠加速度谱。
在地震资料处理及后期的地震数据解释中将非常广泛地利用到速度数据,而速度数据的来源就是速度谱的拾取。目前速度谱的拾取都是采用人工交互方式完成的,时间长、效率低。
目前现有的速度谱自动拾取手段有基于蒙特卡洛法的速度谱自动拾取等,但是其技术复杂,效率低,不能全局解释。
发明内容
本发明的目的在于解决上述现有技术中存在的难题,提供一种地震速度谱自动拾取方法,将速度谱剖面数据转换成二维网格数据,并利用一定的技术手段搜索出每行峰值数据,再设置时窗寻找区域最大值,最终利用约束条件完成速度曲线的连接。
本发明是通过以下技术方案实现的:
一种地震速度谱自动拾取方法,包括:
(1)创建网格数据:将以标准地震数据格式存储的速度谱数据转换成二维网格数据;
(2)设置计算时窗并在计算时窗区域内确定最大值位置:
定义计算时窗的宽度w,并在计算时窗的区域内搜寻最大值Vmax,并以第一个计算时窗中的最大值所在行的位置为速度曲线的起始位置;
所述计算时窗的区域为w×k;
(3)定义速度搜索阈值vh,将计算时窗向下移动一个单元,然后在计算时窗的区域内搜寻最大值Vmax,如果找到的最大值Vmax大于阈值vh,则记录下该最大值Vmax的位置,即其对应的二维数组的下标n和k的值,如果找到的最大值Vmax小于阈值vh,则放弃记录;
一个单元是指计算时窗的宽度w;
(4)完成速度曲线连接,形成未插值的速度曲线;
重复第三步,直到计算时窗依次移动到网格数据的底部,完成所有区域最大值Vmax所在位置的记录,并将各个位置依次连接,形成未插值的速度曲线;
(5)对步骤(4)得到的未插值的速度曲线进行速度曲线插值获得最终的速度曲线。
所述步骤(1)是这样实现的:
标准地震数据格式存储的速度谱数据中,S(n)表示样点值,n表示数组序列,将标准地震数据格式存储的速度谱数据转换成网格数据后,网格数据中用S(n,k)表示样点值,n和k分别表示y轴方向与x轴方向的序列。
所述步骤(2)中的计算时窗的宽度选择范围为3~20行。
所述步骤(3)中阈值vh的取值为第一个计算时窗中的最大值Vmax的二分之一。
所述步骤(5)是采用线性插值的方式对未插值的速度曲线进行插值,即对每两个相邻的位置采用线性插值公式得到这两个位置之间的曲线,然后依次完成所有相邻位置的计算即得到最终的速度曲线。
与现有技术相比,本发明的有益效果是:利用本发明方法进行的速度谱自动拾取,实现简单、计算速度快,有利于缩减地震数据处理周期,此发明能够为常规油气及煤层气的勘探、开发服务。
附图说明
图1速度谱数据的转换
图2计算时窗内获取最大值
图3计算时窗的移动
图4定位阈值的最大值
图5速度曲线连接
图6本发明方法的步骤框图。
具体实施方式
下面结合附图对本发明作进一步详细描述:
如图6所示,本发明地震速度谱自动拾取方法主要分为以下五个步骤:
(1)第一步:创建网格数据:
速度谱数据一般都是以标准地震数据格式(SEG-Y格式)存储的,SEG-Y格式数据的特点是以列的方式记录数据,在本发明中为了便于后续步骤的计算,必须先转换成二维网格数据。其示意如图1所示,SEG-Y数据中S(n)表示样点值,n表示数组序列,网格数据中S(n,k)表示样点值,n和k分别表示y轴方向与x轴方向的序列。
(2)第二步:设置计算时窗并在计算时窗区域内确定最大值位置
定义计算时窗(w)宽度(一般3~20行为佳(对于地震数据处理来说,一般会把时窗定义在3~20之间,一般时窗小计算时间就长但精度高,时窗大计算时间就短但精度低,从程序界面上进行设置即可)),并在计算时窗的区域内(w×k)搜寻最大值Vmax,并以第一个时窗中的最大值所在行位置为速度曲线的起始位置。示意如图2所示。
(3)第三步:移动计算时窗,并重复第二步:
定义速度搜索阈值(vh)(一般用第一个时窗中的最大值Vmax的二分之一)。然后将计算时窗向下移动一个单元(此单元即是计算时窗的宽度w),见图3。并重复第二步。找到大于阈值vh的最大值Vmax并记录下位置(即记录其对应的二维数组的下标:n和k值),如果本时窗区域的最大值小于vh,则放弃记录。见图4所示。
(4)完成速度曲线连接
重复第三步(每次重复是用相同的阈值),直到计算时窗依次移动到网格数据底部,记录下所有区域最大值Vmax所在位置,并依次连接,所连接的曲线即为未插值的速度曲线,如图5所示。
(5)对速度曲线插值
以上四步完成了速度曲线的自动拾取,但此曲线只记录了每一次计算时窗中最大值的位置,对于网格数据还有很多行空缺。本发明采用线性插值的方式对曲线进行插值(即对每两个相邻点采用下面的公式得到这两点之间的曲线,然后依次完成所有相邻点的计算即得到最终的速度曲线)。得到最终的速度曲线。
线性插值法是指使用连接两个已知量的直线来确定在这两个已知量之间的一个未知量的值的方法。
假设已知坐标(x0,y0)与(x1,y1),要得到[x0,x1]区间内某一位置x在直线上的值,则有:
假设方程两边的值为α,那么这个值就是插值系数一从x0到x的距离与从x0到x1距离的比值。由于x值已知,所以可以从公式得到α的值
同样,
这样,在代数上就可以表示成为:
y=(1-α)y0+αy1
或者,
y=y0+α(y1-y0)
这样通过α就可以直接得到y。实际上,即使x不在x0到x1之间并且α也不是介于0到1之间,这个公式也是成立的。在这种情况下,这种方法叫作线性外插。
上述技术方案只是本发明的一种实施方式,对于本领域内的技术人员而言,在本发明公开了应用方法和原理的基础上,很容易做出各种类型的改进或变形,而不仅限于本发明上述具体实施方式所描述的方法,因此前面描述的方式只是优选的,而并不具有限制性的意义。
Claims (5)
1.一种地震速度谱自动拾取方法,其特征在于:所述方法包括:
(1)创建网格数据:将以标准地震数据格式存储的速度谱数据转换成二维网格数据;
(2)设置计算时窗并在计算时窗区域内确定最大值位置:
定义计算时窗的宽度w,并在计算时窗的区域内搜寻最大值Vmax,并以第一个计算时窗中的最大值所在行的位置为速度曲线的起始位置;
所述计算时窗的区域为w×k;
(3)定义速度搜索阈值vh,将计算时窗向下移动一个单元,然后在计算时窗的区域内搜寻最大值Vmax,如果找到的最大值Vmax大于阈值vh,则记录下该最大值Vmax的位置,即其对应的二维数组的下标n和k的值,如果找到的最大值Vmax小于阈值vh,则放弃记录;
一个单元是指计算时窗的宽度w;
(4)完成速度曲线连接,形成未插值的速度曲线;
重复第三步,直到计算时窗依次移动到网格数据的底部,完成所有区域最大值Vmax所在位置的记录,并将各个位置依次连接,形成未插值的速度曲线;
(5)对步骤(4)得到的未插值的速度曲线进行速度曲线插值获得最终的速度曲线。
2.根据权利要求1所述的地震速度谱自动拾取方法,其特征在于:所述步骤(1)是这样实现的:
标准地震数据格式存储的速度谱数据中,S(n)表示样点值,n表示数组序列,将标准地震数据格式存储的速度谱数据转换成网格数据后,网格数据中用S(n,k)表示样点值,n和k分别表示y轴方向与x轴方向的序列。
3.根据权利要求1所述的地震速度谱自动拾取方法,其特征在于:所述步骤(2)中的计算时窗的宽度选择范围为3~20行。
4.根据权利要求1所述的地震速度谱自动拾取方法,其特征在于:所述步骤(3)中阈值vh的取值为第一个计算时窗中的最大值Vmax的二分之一。
5.根据权利要求1所述的地震速度谱自动拾取方法,其特征在于:所述步骤(5)是采用线性插值的方式对未插值的速度曲线进行插值,即对每两个相邻的位置采用线性插值公式得到这两个位置之间的曲线,然后依次完成所有相邻位置的计算即得到最终的速度曲线。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410524637.0A CN105572733B (zh) | 2014-10-08 | 2014-10-08 | 一种地震速度谱自动拾取方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410524637.0A CN105572733B (zh) | 2014-10-08 | 2014-10-08 | 一种地震速度谱自动拾取方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105572733A CN105572733A (zh) | 2016-05-11 |
CN105572733B true CN105572733B (zh) | 2018-11-13 |
Family
ID=55883074
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410524637.0A Active CN105572733B (zh) | 2014-10-08 | 2014-10-08 | 一种地震速度谱自动拾取方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105572733B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US10996358B2 (en) * | 2017-08-18 | 2021-05-04 | Saudi Arabian Oil Company | Image-guided velocity interpolation using a mask cube |
CN107957594A (zh) * | 2017-11-15 | 2018-04-24 | 中国石油集团东方地球物理勘探有限责任公司 | 地震数据的椭圆校正方法、动校正方法及速度谱计算方法 |
CN112540404B (zh) * | 2019-09-20 | 2024-04-12 | 中国石油化工股份有限公司 | 一种基于深度学习的自动速度分析方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073064A (zh) * | 2009-11-25 | 2011-05-25 | 中国石油天然气集团公司 | 一种利用相位信息提高速度谱分辨率的方法 |
CN103984021A (zh) * | 2014-05-22 | 2014-08-13 | 中国海洋大学 | 一种基于系列速度误差剖面的偏移速度分析方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
AU2012388240B2 (en) * | 2012-08-20 | 2016-05-19 | Landmark Graphics Corporation | Methods and systems of incorporating pseudo-surface pick locations in seismic velocity models |
-
2014
- 2014-10-08 CN CN201410524637.0A patent/CN105572733B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102073064A (zh) * | 2009-11-25 | 2011-05-25 | 中国石油天然气集团公司 | 一种利用相位信息提高速度谱分辨率的方法 |
CN103984021A (zh) * | 2014-05-22 | 2014-08-13 | 中国海洋大学 | 一种基于系列速度误差剖面的偏移速度分析方法 |
Non-Patent Citations (3)
Title |
---|
"偏移速度分析与速度反演方法评述";刘百红 等;《CT理论与应用研究》;20131031;第22卷(第4期);第727-740页 * |
"用人工神经网络自动拾取速度谱";董林平 等;《石油地球物理勘探》;19960630;第31卷;第98-103页 * |
"蒙特卡洛法速度谱的自动拾取";张正和 等;《大庆石油地质与开发》;20020630;第21卷(第3期);第79-80页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105572733A (zh) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106353792B (zh) | 一种适用于水力压裂微震震源定位的方法 | |
CN103399346B (zh) | 一种井震联合初始波阻抗建模方法 | |
CN103439738B (zh) | 地震勘探单炮记录异常道识别方法 | |
CN102073067B (zh) | 一种提高地震数据初至波自动拾取效率的方法 | |
CN105572733B (zh) | 一种地震速度谱自动拾取方法 | |
CN112180433B (zh) | 地震初至波拾取方法及装置 | |
CN102901985A (zh) | 一种适用于起伏地表的深度域层速度修正方法 | |
CN102222366A (zh) | 复杂空间曲面的拟合方法 | |
CN112150582A (zh) | 一种面向多模态数据的地质剖面图近似表达方法 | |
CN109100791B (zh) | 基于纵横向空间约束的速度反演方法 | |
CN103268632A (zh) | 一种机载激光雷达扫描生成地形信息的方法 | |
CN105589096B (zh) | 一种基于d-s证据理论的沉积相带划分方法 | |
CN109490961B (zh) | 起伏地表无射线追踪回折波层析成像方法 | |
CN109188516B (zh) | Radon域能量扫描叠加的微地震事件定位方法 | |
CN103513288B (zh) | 一种二维网格数据的补偿方向滤波方法 | |
CN107664771B (zh) | 一种基于相似性系数的微地震全波形定位方法 | |
CN103777244A (zh) | 一种地震裂缝属性体的定量分析方法 | |
CN105572732A (zh) | 一种属性变化率逐次晋级的裂缝发育带检测方法 | |
CN104280766B (zh) | 一种利用局部数据同相轴斜率的直接偏移方法 | |
CN104570100B (zh) | 多子波克希霍夫地震数据偏移方法 | |
CN106646608B (zh) | 一种用于观测系统参数论证的层位模型建立方法 | |
CN110263763B (zh) | 基于岩性中点正弦插值模型分析旋回地层的方法 | |
CN103984021B (zh) | 一种基于系列速度误差剖面的偏移速度分析方法 | |
CN111123356B (zh) | 基于初至信息的异常道智能判识方法 | |
CN113325467A (zh) | 一种基于槽波频散特征的微地震震源定位方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |