CN105319586A - 一种速度谱三维滤波优化速度分析方法 - Google Patents
一种速度谱三维滤波优化速度分析方法 Download PDFInfo
- Publication number
- CN105319586A CN105319586A CN201410316146.7A CN201410316146A CN105319586A CN 105319586 A CN105319586 A CN 105319586A CN 201410316146 A CN201410316146 A CN 201410316146A CN 105319586 A CN105319586 A CN 105319586A
- Authority
- CN
- China
- Prior art keywords
- velocity
- dimensional
- spectrum
- velocity spectrum
- grid
- 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.)
- Granted
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种针对速度谱三维滤波的优化速度分析方法。包括：对选取速度分析线上的每个点都产生速度谱；对速度谱数据进行道头修改，将二维速度谱构建成三维叠加数据体；对构建的三维叠后数据体加三维处理网格，网格面元中心Y坐标范围与速度分析线的面元中心Y坐标范围一致，网格X坐标增量与速度分析线道间距一致，且横测线方向面元个数与生成速度谱时定义的速度扇区扫描个数一致；对该网格化的三维数据体进行三维滤波处理；通过道头恢复，将构建的三维叠加数据体恢复为二维速度谱。本发明通过道头修改将二维速度谱构建成三维叠后数据进行三维滤波处理，更直接而有针对性地提高了速度谱的信噪比，从而提高了速度分析精度。
Description
技术领域
本发明涉及地震资料处理技术领域，特别是涉及到一种针对速度谱三维滤波的优化速度分析方法。
背景技术
地震成像对速度的分析精度要求越来越高，速度谱信噪比的高低直接影响着速度谱能量团聚焦程度，从而制约着速度拾取精度。现存的速度优化方法均是通过对地震道集优化来实现，缺少直接针对速度谱的滤波方法。
发明内容
本发明的目的是为了填补现有技术的空白，提供一种速度谱三维滤波优化速度分析方法，实现直接针对速度谱进行三维滤波，更直观而有针对性的提高速度谱的信噪比，从而达到提高速度分析精度。
本发明将速度分析线上逐点排列的速度谱看成是三维地震叠加数据体，将空间连续分布的速度谱能量团视作此三维叠后数据体中的有效地震信号，通过针对此速度谱构成的数据体的三维滤波处理来改善速度谱的信噪比，增强速度谱能量团的聚焦性，以此来提高速度分析精度。
本发明的目的可通过如下技术措施来实现：
在步骤1中，对选取速度分析线上的每个点都产生速度谱。
在步骤2中，对速度谱数据进行道头修改，将速度分析线方向视为纵测线方向，将速度谱中速度轴方向改为横测线方向，连同时间方向，将二维速度谱构建成三维叠加数据体。
在步骤3中，对构建的三维叠后数据体加三维处理网格，网格面元中心Y坐标范围与速度分析线的面元中心Y坐标范围一致，网格X坐标增量与速度分析线道间距一致，且横测线方向面元个数与生成速度谱时定义的速度扇区扫描个数一致。
在步骤4中，对该网格化的三维数据体进行三维滤波处理。
在步骤5中，通过道头恢复，将构建的三维叠加数据体恢复为二维速度谱。
本发明中的速度谱三维滤波方法通过道头修改将二维速度谱构建成三维叠后数据进行三维滤波处理，更直接而有针对性地提高了速度谱的信噪比，从而提高了速度分析精度。
附图说明
图1为本发明一种实施例的流程图；
图2为采用本发明在该工区实施例中构建的三维数据体示意图；
图3为本发明在该工区实施例中对三维叠加数据体使用的三维处理网格；
图4为在该工区采用常规速度谱法与采用本发明方法产生的速度分析线上某个点的速度谱对比图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂，下文特以胜利油田某工区为例，并配合所附图式，作详细说明如下。
图1为本发明在该工区实施例中速度谱三维滤波方法的流程图。在步骤101，由常规速度谱法生成该工区一条速度分析线上每个点的速度谱。该速度分析线上共有501个点，道间距25米，第一个和最后一个点的（X，Y）坐标分别为（8450,36475），（8450,48975）。在该步骤中，定义的速度扇区扫描个数为121，流程进入到步骤102。
在步骤102，修改速度谱数据道头。步骤101中定义的速度扇区扫描个数为121，即每个点的速度谱是121道，将道号与道间距的乘积赋给每个速度谱道的X坐标。至此，速度谱可以被视为一个三维叠后数据体，速度分析线方向视为纵测线方向，速度谱中速度轴方向改为了横测线方向，第一个和最后一个速度谱点的坐标分别变为（25,36475），（3025,48975），完成了三维数据体构建（纵测线，横测线，时间）。流程进入到步骤103。
在步骤103，对构建的三维叠后数据体加三维处理网格，网格面元中心Y坐标范围与速度分析线的面元中心Y坐标范围一致，网格X坐标增量与速度分析线道间距一致，且横测线方向面元个数与生成速度谱时定义的速度扇区扫描个数一致，按照上述要求，该工区三维处理网格四点坐标为：A（25,36475）、B（25,48975）、C（3025,36475）、D（3025,48975）。流程进入到步骤104。
在步骤104，对速度谱进行三维滤波处理。流程进入到步骤105。
在步骤105，恢复原道头：恢复原速度分析线X坐标，即所有的速度谱道X坐标为8450。滤波后的三维数据体恢复为二维速度谱。至此，该速度谱就可以用于速度分析。流程结束。
图2为采用本发明在该工区实施例中构建的三维叠后数据体示意图。其中，速度分析线方向看作纵测线方向，速度谱中速度轴方向改为横测线方向，与时间方向一起构建成三维叠后数据体。
图3为本发明在该工区实施例中对三维叠加数据体使用的三维处理网格；
图4为在该工区采用常规速度谱法与采用本发明方法产生的速度分析线上某个点的速度谱对比图。其中，左图为常规速度谱法生成的速度谱，右图为该发明方法生成的速度谱。从图4中可以看出，本发明中的速度谱三维滤波方法提高了速度谱的信噪比，增强了复杂地质构造处速度谱的聚焦性，能够提高速度拾取的精度，进而改善复杂地质构造的叠加和偏移成像效果。
Claims (2)
1.一种速度谱三维滤波优化速度分析方法，其特征是：
（1）对选取速度分析线上的每个点都产生速度谱；
（2）对速度谱数据进行道头修改，将速度分析线方向视为纵测线方向，将速度谱中速度轴方向改为横测线方向，连同时间方向，将二维速度谱构建成三维叠加数据体；
（3）对构建的三维叠后数据体加三维处理网格，网格面元中心Y坐标范围与速度分析线的面元中心Y坐标范围一致，网格X坐标增量与速度分析线道间距一致，且横测线方向面元个数与生成速度谱时定义的速度扇区扫描个数一致；
（4）对该网格化的三维数据体进行三维滤波处理；
（5）通过道头恢复，将构建的三维叠加数据体恢复为二维速度谱。
2.根据权利要求1所述的速度谱三维滤波优化速度分析方法，其特征是：在（1）中，选取速度分析线由常规速度谱法生成该工区一条速度分析线。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410316146.7A CN105319586B (zh) | 2014-07-04 | 2014-07-04 | 一种速度谱三维滤波优化速度分析方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410316146.7A CN105319586B (zh) | 2014-07-04 | 2014-07-04 | 一种速度谱三维滤波优化速度分析方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105319586A true CN105319586A (zh) | 2016-02-10 |
| CN105319586B CN105319586B (zh) | 2018-02-02 |
Family
ID=55247389
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410316146.7A Active CN105319586B (zh) | 2014-07-04 | 2014-07-04 | 一种速度谱三维滤波优化速度分析方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105319586B (zh) |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040167650A1 (en) * | 2003-02-21 | 2004-08-26 | Fujitsu Limited | Partial reprojection method and device in three-dimensional CAD system and computer program |
| US20050128198A1 (en) * | 2003-11-17 | 2005-06-16 | Samsung Electronics Co., Ltd. | Method and apparatus for generating three-dimensional finite element mesh |
| US20070100492A1 (en) * | 2005-10-21 | 2007-05-03 | Mamoru Idaka | Three dimensional processing data setting system, method for setting three-dimensional processing data, computer program for setting three-dimensional processing data, medium with three-dimensional processing data stored therein that is readable by computer and laser processing equipment operated by the three-dimensional data |
| CN101852864A (zh) * | 2009-04-03 | 2010-10-06 | 中国石油集团东方地球物理勘探有限责任公司 | 一种利用地表一致性统计频谱分析技术处理海量地震数据的方法 |
| CN103217713A (zh) * | 2012-01-19 | 2013-07-24 | 中国石油化工股份有限公司 | 油气勘探地震速度分析数据优化方法 |
-
2014
- 2014-07-04 CN CN201410316146.7A patent/CN105319586B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20040167650A1 (en) * | 2003-02-21 | 2004-08-26 | Fujitsu Limited | Partial reprojection method and device in three-dimensional CAD system and computer program |
| US20050128198A1 (en) * | 2003-11-17 | 2005-06-16 | Samsung Electronics Co., Ltd. | Method and apparatus for generating three-dimensional finite element mesh |
| US20070100492A1 (en) * | 2005-10-21 | 2007-05-03 | Mamoru Idaka | Three dimensional processing data setting system, method for setting three-dimensional processing data, computer program for setting three-dimensional processing data, medium with three-dimensional processing data stored therein that is readable by computer and laser processing equipment operated by the three-dimensional data |
| CN101852864A (zh) * | 2009-04-03 | 2010-10-06 | 中国石油集团东方地球物理勘探有限责任公司 | 一种利用地表一致性统计频谱分析技术处理海量地震数据的方法 |
| CN103217713A (zh) * | 2012-01-19 | 2013-07-24 | 中国石油化工股份有限公司 | 油气勘探地震速度分析数据优化方法 |
Non-Patent Citations (2)
| Title |
|---|
| 包吉山等: "二维倾斜叠加速度分析", 《石油地球物理勘探》 * |
| 郭树祥: "一种优化数据量以提高地震速度分析精度的方法", 《石油物探》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105319586B (zh) | 2018-02-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN102930509B (zh) | 一种机载激光点云数据的智能化滤波方法 | |
| CN103454686B (zh) | 基于地层切片的小尺度沉积相进行储层预测的方法及系统 | |
| CN104459768B (zh) | 一种基于可视化的三维空间目标地质体追踪方法 | |
| CN103116192A (zh) | 一种缝洞型碳酸盐岩油藏储集体建模方法 | |
| CN107490808B (zh) | 一种高可靠性地震勘探观测系统的建立方法 | |
| CN106526677A (zh) | 一种海上自适应压制鬼波的宽频逆时偏移成像方法 | |
| CN105093319B (zh) | 基于三维地震数据的地面微地震静校正方法 | |
| CN104483705A (zh) | 一种三维剩余静校正方法 | |
| CN103592698A (zh) | 一种基于地震属性的观测系统评价方法及装置 | |
| CN105549077B (zh) | 基于多级多尺度网格相似性系数计算的微震震源定位方法 | |
| CN106154315A (zh) | 一种对称均匀的高密度地震采集方法 | |
| CN103217710A (zh) | 获取共偏移距共方位角域道集的方法及系统 | |
| CN102353986B (zh) | 含逆断层的复杂地质构造的全三维联动地震解释方法 | |
| CN104267432B (zh) | 一种基于规则化的转换波共转换点道集高精度抽取方法 | |
| CN106908836B (zh) | 采集脚印压制方法及系统 | |
| CN106896408B (zh) | 一种角度域叠前时间偏移方法 | |
| CN103576194B (zh) | 一种深度域复杂地表地震构造成像方法 | |
| CN104155690B (zh) | 基于椭球展开的三维地震数据叠加速度求取方法 | |
| CN105738949B (zh) | 一种用于时移地震的九面元一致性并行处理方法 | |
| CN105319586A (zh) | 一种速度谱三维滤波优化速度分析方法 | |
| CN105487106A (zh) | 一种基于高斯射线束目的层能量照明的补炮方法 | |
| CN104765066B (zh) | 地震三维波速扫描聚焦成像方法 | |
| CN105425290B (zh) | 一种叠前时间偏移的方法及装置 | |
| CN109270576A (zh) | 多通道三维地震物理模拟数据采集方法及装置 | |
| CN106338762B (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 |