CN103487829A - 两参数转换波各向异性叠前时间偏移方法 - Google Patents

两参数转换波各向异性叠前时间偏移方法 Download PDF

Info

Publication number
CN103487829A
CN103487829A CN201310406741.5A CN201310406741A CN103487829A CN 103487829 A CN103487829 A CN 103487829A CN 201310406741 A CN201310406741 A CN 201310406741A CN 103487829 A CN103487829 A CN 103487829A
Authority
CN
China
Prior art keywords
transformed wave
time migration
parameter
wave
anisotropy
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
Application number
CN201310406741.5A
Other languages
English (en)
Other versions
CN103487829B (zh
Inventor
罗红明
何光明
李忠
陈爱萍
金德刚
刘奇琳
王栋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd filed Critical Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
Priority to CN201310406741.5A priority Critical patent/CN103487829B/zh
Publication of CN103487829A publication Critical patent/CN103487829A/zh
Application granted granted Critical
Publication of CN103487829B publication Critical patent/CN103487829B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开一种两参数转换波各向异性叠前时间偏移方法,包括以下步骤:(a)获取转换波叠前时间偏移道集;(b)通过扫描转换波偏移速度谱来获取转换波叠前时间偏移速度Vps;(c)通过扫描转换波等效各向异性参数谱来获取转换波等效各向异性参数Kappa;(d)在成像点位置进行两参数转换波各向异性叠前时间偏移处理。

Description

两参数转换波各向异性叠前时间偏移方法
技术领域
本发明涉及石油天然气地震勘探领域,更具体地讲,涉及一种适用于二维和三维转换波叠前时间偏移处理的高效率的两参数转换波各向异性叠前时间偏移方法。
背景技术
随着石油天然气勘探开发的不断深入,地震勘探技术已从常规的三维地震勘探向三维三分量勘探方向发展,而转换波各向异性精确偏移成像已经成为转换波成果评价的关键技术和难题。
具体地讲,转换波各向异性叠前偏移成像方法研究是石油地球物理勘探中长期存在并且目前尚未完全解决的前沿性研究课题,其中,各向异性参数的构建方法直接影响着最终的成像质量。因此,为了获取转换波精细偏移成像的效果,目前主要使用四个参数的转换波各向异性叠前时间偏移成像方法,并且取得了很多理论和实际应用成果。目前关于转换波各向异性叠前时间偏移的研究大多存在以下两个方面的问题。
第一,转换波各向异性叠前时间偏移成像参数繁多、获取繁琐的问题。
如上所述,转换波各向异性叠前时间偏移成像目前主要使用四个参数的方法,需要先后获取转换波速度Vps、等效速度比参数Gamma0、有效速度比参数Gamma_effect以及各向异性参数Eta,每次更新都需要对四个参数重新计算和更新,操作极其繁琐。
第二,各参数对成像分辨率影响不一致。
在转换波叠前时间偏移中,上述四个参数对偏移成像效果影响作用不一致。影响作用由大到小依次为转换波速度Vps、等效速度比参数Gamma0、有效速度比参数Gamma_effect和各向异性参数Eta。因此,转换波速度Vps最为关键,而各向异性参数Eta常常影响有限。
发明内容
因此,本发明的目的在于提供一种快速高效的两参数转换波各向异性叠前时间偏移方法。
根据本发明的一方面,提供一种两参数转换波各向异性叠前时间偏移方法,包括以下步骤:(a)获取转换波叠前时间偏移道集;(b)通过扫描转换波偏移速度谱来获取转换波叠前时间偏移速度Vps;(c)通过扫描转换波等效各向异性参数谱来获取转换波等效各向异性参数Kappa;(d)在成像点位置进行两参数转换波各向异性叠前时间偏移处理。
优选地,针对每个成像点执行步骤(a)至步骤(d)。
优选地,在步骤(d)中,根据以下等式进行两参数转换波各向异性叠前时间偏移处理:
t ps 2 = t ps 0 2 + ( 2 h ) 2 V ps 2 - 2 k eff ( 2 h ) 4 V ps 2 [ t ps 0 2 + ( 0.1 + 2.7 k eff ) · ( 2 h ) 2 ] ,
其中,tps为成像点位置转换波旅行时,tps0为成像点位置转换波自激自收旅行时间,keff为转换波等效各向异性参数,Vps为转换波叠前时间偏移速度,2h为转换波炮点到检波点的距离。
与传统技术相比,根据本发明的两参数转换波各向异性叠前时间偏移方法在考虑的相关影响信息因素丰富的基础上,简化了输入参数,使得偏移成像效率更高。
附图说明
通过下面结合附图对实施例进行的描述,本发明的这些和/或其他方面和优点将会变得清楚和更易于理解,其中:
图1是示出根据本发明的示例性实施例的两参数转换波各向异性叠前时间偏移方法的流程图。
具体实施方式
现在对本发明实施例进行详细的描述,其示例表示在附图中,其中,相同的标号始终表示相同部件。下面通过参照附图对实施例进行描述以解释本发明。
图1是示出根据本发明的示例性实施例的两参数转换波各向异性叠前时间偏移方法的流程图。
参照图1,在步骤S101中,获取转换波叠前时间偏移道集。这里,可是有各种现有方法直接获得转换波叠前时间偏移道集。
接下来,在步骤S102中,通过扫描转换波偏移速度谱来获取转换波叠前时间偏移速度Vps。可通过扫描转换波偏移速度谱,然后在谱上进行自动或人工拾取成像点位置的速度值,来获取转换波叠前时间偏移速度Vps。
然后,在步骤S103中,通过扫描转换波等效各向异性参数谱来获取转换波等效各向异性参数Kappa。如上所述,可通过扫描转换波等效各向异性参数谱,然后在谱上进行自动或人工拾取成像点位置的等效各向异性参数,来获取转换波等效各向异性参数Kappa。根据本发明的示例性实施例,可同时执行步骤S102和S103,或者先执行步骤S102然后执行步骤S103,或者先执行步骤S103然后执行步骤S102。
最后,在步骤S104中,在成像点位置进行两参数转换波各向异性叠前时间偏移处理。具体地讲,可根据以下等式(1)进行两参数转换波各向异性叠前时间偏移处理:
t ps 2 = t ps 0 2 + ( 2 h ) 2 V ps 2 - 2 k eff ( 2 h ) 4 V ps 2 [ t ps 0 2 + ( 0.1 + 2.7 k eff ) · ( 2 h ) 2 ] - - - ( 1 ) ,
其中,tps为成像点位置转换波旅行时,tps0为成像点位置转换波自激自收旅行时间,keff为转换波等效各向异性参数,Vps为转换波叠前时间偏移速度,2h为转换波炮点到检波点的距离。这里,tps0和2h可根据现有技术容易地确定,Vps和keff分别为在步骤S102和S103中获取的参数。
根据本发明的示例性实施例,可针对各个成像点执行步骤S101至S104,从而在所有成像点位置完成两参数转换波各向异性叠前时间偏移处理。
如上所述的两参数转换波各向异性叠前时间偏移方法特别适用于转换波叠前时间偏移成像处理,具有计算步骤简便、计算效率高和处理效果好等特点,在目前西部地区的转换波震资料处理中具有广阔的应用前景。
将理解:可按照硬件、软件或硬件和软件的组合的形式来实现如上所述的本发明的示例性实施例。可在非暂时性计算机可读存储介质中存储任何所述软件。非暂时性计算机可读存储介质存储一个或多个程序(软件模块),所述一个或多个程序包括指令,其中,当由在电子装置中的一个或多个处理器执行所述指令时,所述指令促使电子装置执行本发明的方法。任何所述软件可按照易失性或非易失性存储器(诸如,不论是否是可擦除的或可重写的存储装置,比如只读存储器(ROM))的形式被存储,或者按照存储器(诸如,随机存取存储器(RAM)、存储芯片、装置或集成电路)的形式被存储,或者在光学或磁可读介质(诸如,CD、DVD、磁盘或磁带等)上被存储。将理解:存储装置和存储介质是适于存储包括指令的程序的机器可读存储器的实施例,其中,当所述指令被执行时,所述指令实施本发明的实施例。因此,本发明的实施例提供一种包括代码的程序和存储所述程序的机器可读存储器,其中,所述代码用于实施如上所述的根据本发明的示例性实施例的两参数转换波各向异性叠前时间偏移方法。进一步地说,所述程序可经由任何介质(诸如通过有线或无线连接携带的通信信号)被电子地传输,并且本发明的实施例适当地包含所述程序。
虽然已经显示和描述了一些实施例,但是本领域技术人员应该理解,在不脱离本发明的原理和精神的情况下,可以对这些实施例进行修改,本发明的范围由权利要求及其等同物限定。

Claims (3)

1.一种两参数转换波各向异性叠前时间偏移方法,包括以下步骤:
(a)获取转换波叠前时间偏移道集;
(b)通过扫描转换波偏移速度谱来获取转换波叠前时间偏移速度Vps;
(c)通过扫描转换波等效各向异性参数谱来获取转换波等效各向异性参数Kappa;
(d)在成像点位置进行两参数转换波各向异性叠前时间偏移处理。
2.根据权利要求1所述的两参数转换波各向异性叠前时间偏移方法,其中,针对每个成像点执行步骤(a)至步骤(d)。
3.根据权利要求1所述的两参数转换波各向异性叠前时间偏移方法,其中,在步骤(d)中,根据以下等式进行两参数转换波各向异性叠前时间偏移处理:
t ps 2 = t ps 0 2 + ( 2 h ) 2 V ps 2 - 2 k eff ( 2 h ) 4 V ps 2 [ t ps 0 2 + ( 0.1 + 2.7 k eff ) · ( 2 h ) 2 ] ,
其中,tps为成像点位置转换波旅行时,tps0为成像点位置转换波自激自收旅行时间,keff为转换波等效各向异性参数,Vps为转换波叠前时间偏移速度,2h为转换波炮点到检波点的距离。
CN201310406741.5A 2013-09-09 2013-09-09 两参数转换波各向异性叠前时间偏移方法 Active CN103487829B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310406741.5A CN103487829B (zh) 2013-09-09 2013-09-09 两参数转换波各向异性叠前时间偏移方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310406741.5A CN103487829B (zh) 2013-09-09 2013-09-09 两参数转换波各向异性叠前时间偏移方法

Publications (2)

Publication Number Publication Date
CN103487829A true CN103487829A (zh) 2014-01-01
CN103487829B CN103487829B (zh) 2016-05-04

Family

ID=49828198

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310406741.5A Active CN103487829B (zh) 2013-09-09 2013-09-09 两参数转换波各向异性叠前时间偏移方法

Country Status (1)

Country Link
CN (1) CN103487829B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110208854A (zh) * 2019-06-28 2019-09-06 中国地质调查局沈阳地质调查中心 一种vti介质中等效各向异性参数的获取方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6128580A (en) * 1998-04-17 2000-10-03 Bp Amoco Corporation Converted-wave processing in many-layered anisotropic media
GB2402214A (en) * 2003-05-29 2004-12-01 Pgs Americas Inc Seismic velocity update for anisotropic depth migration
CN102540251A (zh) * 2010-12-16 2012-07-04 中国石油天然气集团公司 二维的横向各向异性介质hti叠前深度偏移建模方法及装置
CN102721977A (zh) * 2012-05-31 2012-10-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 高精度转换波各向异性叠加速度分析方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6128580A (en) * 1998-04-17 2000-10-03 Bp Amoco Corporation Converted-wave processing in many-layered anisotropic media
GB2402214A (en) * 2003-05-29 2004-12-01 Pgs Americas Inc Seismic velocity update for anisotropic depth migration
CN102540251A (zh) * 2010-12-16 2012-07-04 中国石油天然气集团公司 二维的横向各向异性介质hti叠前深度偏移建模方法及装置
CN102721977A (zh) * 2012-05-31 2012-10-10 中国石油集团川庆钻探工程有限公司地球物理勘探公司 高精度转换波各向异性叠加速度分析方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
马昭军等: "《叠前时间偏移在三维转换波资料处理中的应用》", 《石油物探》 *
黄中玉等: "《层状各向异性介质转换波克希霍夫叠前时间偏移》", 《地球物理学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110208854A (zh) * 2019-06-28 2019-09-06 中国地质调查局沈阳地质调查中心 一种vti介质中等效各向异性参数的获取方法
CN110208854B (zh) * 2019-06-28 2020-08-14 中国地质调查局沈阳地质调查中心 一种vti介质中等效各向异性参数的获取方法

Also Published As

Publication number Publication date
CN103487829B (zh) 2016-05-04

Similar Documents

Publication Publication Date Title
Maggi et al. An automated time-window selection algorithm for seismic tomography
CN102221709B (zh) 基于地层参数信息的速度分析与动校正方法
CN102721977B (zh) 高精度转换波各向异性叠加速度分析方法
CN103439738B (zh) 地震勘探单炮记录异常道识别方法
CN104730574B (zh) 构建近地表结构模型的方法
CN103105624B (zh) 基于数据库技术的纵横波时差定位方法
CN104267429A (zh) 确定地层压力的方法及装置
CN104199103B (zh) 转换波静校正方法及装置
CN105629303A (zh) 基于岩石物理的叠前裂缝定量预测方法及系统
CN108828668B (zh) 一种叠前时间偏移数据处理方法及装置
Liu et al. Stochastic inversion of facies and reservoir properties based on multi-point geostatistics
CN104459768A (zh) 一种基于可视化的三维空间目标地质体追踪方法
CN104656139A (zh) 一种转换波的静校正方法及设备
CN102520446A (zh) 高精度三维层序自动识别系统
CN104502996A (zh) 一种密度曲线校正方法及系统
CN106772593A (zh) 绕射波的成像方法及装置
CN106405644A (zh) 裂缝确定方法和装置
CN104133244B (zh) 一种震源横波勘探的近地表静校正处理方法
CN105445793B (zh) 一种确定坏道数据的方法及装置
CN104267432A (zh) 一种基于规则化的转换波共转换点道集高精度抽取方法
CN103513279B (zh) 一种基于地震波波动方程的照明分析计算方法及计算装置
CN102914797A (zh) 一种获得地层各向异性系数的方法及装置
CN103487829B (zh) 两参数转换波各向异性叠前时间偏移方法
CN105445788B (zh) 一种基于模型和全局寻优的速度谱自动解释方法
CN103576201A (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
C14 Grant of patent or utility model
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20180223

Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee after: BGP INC., CHINA NATIONAL PETROLEUM Corp.

Address before: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan

Patentee before: CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20200917

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Patentee after: BGP Inc., China National Petroleum Corp.

Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189

Patentee before: BGP Inc., China National Petroleum Corp.

TR01 Transfer of patent right