CN102230973A - 一种三维分步傅立叶粘滞声波深度偏移方法 - Google Patents
一种三维分步傅立叶粘滞声波深度偏移方法 Download PDFInfo
- Publication number
- CN102230973A CN102230973A CN2011100920113A CN201110092011A CN102230973A CN 102230973 A CN102230973 A CN 102230973A CN 2011100920113 A CN2011100920113 A CN 2011100920113A CN 201110092011 A CN201110092011 A CN 201110092011A CN 102230973 A CN102230973 A CN 102230973A
- Authority
- CN
- China
- Prior art keywords
- depth
- stack
- omega
- sound wave
- partiald
- 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.)
- Pending
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种三维分步傅立叶粘滞声波深度偏移方法,处理步骤为:对叠前地震数据进行预处理,建立深度速度模型;使用频谱比法求取品质因子Q值;输入叠前或叠后地震数据体、深度速度模型、品质因子Q值模型,输出粘滞声波波动方程叠前或叠后深度偏移结果;和常规波动方程深度偏移结果进行对比,检验结果正确性;将对比的两个数据体都转换到时间域,进行频谱分析,检验分析粘滞声波波动方程深度偏移方法提高分辨率的能力;与传统的声波偏移技术相比,偏移结果能有效地补偿岩石吸收对地震波能量的影响,获得更真实、更精确地下成像效果,使最终成像的能量真正体现地下介质的岩性和物性变化,为AVO分析和利用地震数据进行油气检测提供可靠资料。
Description
技术领域:
本发明涉及一种三维分步傅立叶粘滞声波深度偏移方法,属于地震勘探数据处理过程中的偏移成像技术,具体说是一种对地震波的吸收和频散进行补偿的深度偏移方法。
背景技术:
现有的商业软件使用的粘弹性补偿是目前工业界仍广泛使用的反Q滤波技术。但地震波的衰减不仅与传播时间有关,与子波频率、速度场分布、射线路径、地下构造展布等都有关系,仅利用反Q滤波不能恢复地震波场的本来特征。随着勘探技术的不断发展,地震成像的目标由完成构造成像任务逐步转移到寻找低幅构造以及岩性异常检测上来。地下介质的粘弹性行为,作为影响地震波传播过程中振幅、频率、相位等变化的主要因素,需要在成像过程中考虑和补偿。
发明内容
本发明的目的是提供一种三维分步傅立叶粘滞声波深度偏移方法,以粘滞性声波波动方程为理论基础,将地震波的吸收和频散考虑到偏移延拓算子中去,在波场延拓同时,补偿地震波的吸收和频散。
本发明的目的是这样实现的:一种三维分步傅立叶粘滞声波深度偏移方法,主要利用三维分步傅立叶粘滞声波深度偏移技术,处理步骤为:
(1)对叠前地震数据进行预处理,建立深度速度模型;
(2)根据实际地质情况,使用频谱比法求取品质因子Q值;
(3)输入叠前或叠后地震数据体、深度速度模型、品质因子Q值模型,输出粘滞声波波动方程叠前或叠后深度偏移结果;
(4)和常规波动方程深度偏移结果进行对比,检验结果正确性;
(5)将对比的两个数据体都转换到时间域,进行频谱分析,检验分析粘滞声波波动方程深度偏移方法提高分辨率的能力。
频率-空间域粘滞声波方程:
其中ω是角频率,ρ是介质密度,M(x,z,ω)是复粘弹性模量。
本发明的优点是:与传统的声波偏移技术相比,偏移结果能有效地补偿岩石吸收对地震波能量的影响,获得更真实、更精确地下成像效果,使最终成像的能量真正体现地下介质的岩性和物性变化,为AVO分析和利用地震数据进行油气检测提供可靠资料。
附图说明
图1为本发明的简单流程图。
图2为本发明实施例中显示的常规波动方程深度偏移测试结果。
图3为本发明实施例中显示的粘滞声波波动方程深度偏移测试结果。
具体实施方式
如图1所示,本发明将地震波的吸收和频散考虑到偏移延拓算子中去,在偏移成像的同时,补偿地震波的吸收和频散,其具体实现步骤包括:
1、对叠前地震数据进行预处理,建立深度速度模型。
2、根据本地区实际地质情况,使用频谱比法求取品质因子Q值。
3、输入叠前或叠后地震数据体、深度速度模型、品质因子Q值模型,输出粘滞声波波动方程叠前或叠后深度偏移结果。
4、和常规波动方程深度偏移结果进行对比,检验结果正确性。
5、将对比的两个数据体都转换到时间域,进行频谱分析,检验分析粘滞声波波动方程深度偏移方法提高分辨率的能力。
频率一空间域粘滞声波方程:
其中ω是角频率,ρ是介质密度,M(x,z,ω)是复粘弹性模量。
本发明采用的具体方法原理如下:
对于地震波在粘弹性介质中传播时的衰减,通常用地层的品质因子Q来描述,Q值是岩石的物理性质之一,定义为地震波传播一个波长的距离后,其能量衰减的程度。非完全弹性介质对能量的损耗也可以通过假设弹性模量为复数来表征,此时Q值可用复模量的实部与虚部之比来定义,即:
M(ω)=MR(ω)+iMI(ω),复速度v(ω)表示为:
ω是角频率,ρ是介质密度,M(x,z,ω)是复粘弹性模量。
根据Kjartansson(1979)提出的关于地震波在粘弹性介质中传播的相关理论,相速度v(ω)、参考角频率ω0及其对应相速度v0、介质常数γ、品质因子Q和复波数KC有如下关系式成立:
对上式做三维傅立叶变换,得到下式:
式中kx,ky分别是x,y对应的波数,进一步得到:
令上式可分裂为下面两个方程,分别代表相移项和时移项。
利用分步傅立叶方法在频率波数和频率空间双域交替求解上式,完成粘滞声波深度偏移。本发明既可以应用于叠后记录的处理,也可以应用于叠前记录的处理。本发明允许Q值在空间域逐点变化,在一些高分辨率勘探地区会发挥作用有用。
如图2、图3所示:在某探区的实际三维数据进行测试,处理范围:INLINE方向有269条线;CROSSLINE方向有289个CDP点;采样间隔:2毫秒;记录长度:7000毫秒;偏移深度6000米;深度采样间隔为10米。上图显示的处理结果都是纯波剖面,没有经过任何处理修饰。可以看出对于浅中层的小断块,小砂体来说,分辨率都有很大改进,时间域频谱分析结果,主频提高15HZ,频带明显展宽。
Claims (2)
1.一种三维分步傅立叶粘滞声波深度偏移方法,主要利用三维分步傅立叶粘滞声波深度偏移技术,其特征在于具体步骤为:
(1)对叠前地震数据进行预处理,建立深度速度模型;
(2)根据实际地质情况,使用频谱比法求取品质因子Q值;
(3)输入叠前或叠后地震数据体、深度速度模型、品质因子Q值模型,输出粘滞声波波动方程叠前或叠后深度偏移结果;
(4)和常规波动方程深度偏移结果进行对比,检验结果正确性;
(5)将对比的两个数据体都转换到时间域,进行频谱分析,检验分析粘滞声波波动方程深度偏移方法提高分辨率的能力。
2.根据权利要求1所述的一种三维分步傅立叶粘滞声波深度偏移方法,其特征在于频率-空间域粘滞声波方程:
其中ω是角频率,ρ是介质密度,M(x,z,ω)是复粘弹性模量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100920113A CN102230973A (zh) | 2011-04-13 | 2011-04-13 | 一种三维分步傅立叶粘滞声波深度偏移方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2011100920113A CN102230973A (zh) | 2011-04-13 | 2011-04-13 | 一种三维分步傅立叶粘滞声波深度偏移方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102230973A true CN102230973A (zh) | 2011-11-02 |
Family
ID=44843559
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2011100920113A Pending CN102230973A (zh) | 2011-04-13 | 2011-04-13 | 一种三维分步傅立叶粘滞声波深度偏移方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102230973A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103777240A (zh) * | 2012-10-25 | 2014-05-07 | 中国石油化工股份有限公司 | 一种可控震源分频扫描频带补偿方法 |
RU2526794C2 (ru) * | 2012-08-29 | 2014-08-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Российский государственный университет нефти и газа имени И.М.Губкина" (РГУ нефти и газа имени И.М.Губкина) | Способ определения упругих свойств горных пород на основе пластовой адаптивной инверсии сейсмических данных |
CN105388518A (zh) * | 2014-09-04 | 2016-03-09 | 中国石油化工股份有限公司 | 一种质心频率与频谱比联合的井中地震品质因子反演方法 |
CN106324663A (zh) * | 2015-06-17 | 2017-01-11 | 中国石油化工股份有限公司 | 一种品质因子的获取方法 |
CN106443786A (zh) * | 2016-11-14 | 2017-02-22 | 中国科学院地质与地球物理研究所 | 基于地面接收的反射地震资料的q值场建模方法 |
CN106662664A (zh) * | 2014-06-17 | 2017-05-10 | 埃克森美孚上游研究公司 | 快速粘声波和粘弹性全波场反演 |
CN107024716A (zh) * | 2016-02-01 | 2017-08-08 | 中国石油化工股份有限公司 | 一种地震波场吸收补偿成像方法及系统 |
-
2011
- 2011-04-13 CN CN2011100920113A patent/CN102230973A/zh active Pending
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2526794C2 (ru) * | 2012-08-29 | 2014-08-27 | Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Российский государственный университет нефти и газа имени И.М.Губкина" (РГУ нефти и газа имени И.М.Губкина) | Способ определения упругих свойств горных пород на основе пластовой адаптивной инверсии сейсмических данных |
CN103777240A (zh) * | 2012-10-25 | 2014-05-07 | 中国石油化工股份有限公司 | 一种可控震源分频扫描频带补偿方法 |
CN103777240B (zh) * | 2012-10-25 | 2016-08-10 | 中国石油化工股份有限公司 | 一种可控震源分频扫描频带补偿方法 |
CN106662664A (zh) * | 2014-06-17 | 2017-05-10 | 埃克森美孚上游研究公司 | 快速粘声波和粘弹性全波场反演 |
CN105388518A (zh) * | 2014-09-04 | 2016-03-09 | 中国石油化工股份有限公司 | 一种质心频率与频谱比联合的井中地震品质因子反演方法 |
CN105388518B (zh) * | 2014-09-04 | 2018-01-05 | 中国石油化工股份有限公司 | 一种质心频率与频谱比联合的井中地震品质因子反演方法 |
CN106324663A (zh) * | 2015-06-17 | 2017-01-11 | 中国石油化工股份有限公司 | 一种品质因子的获取方法 |
CN106324663B (zh) * | 2015-06-17 | 2018-10-02 | 中国石油化工股份有限公司 | 一种品质因子的获取方法 |
CN107024716A (zh) * | 2016-02-01 | 2017-08-08 | 中国石油化工股份有限公司 | 一种地震波场吸收补偿成像方法及系统 |
CN107024716B (zh) * | 2016-02-01 | 2019-04-02 | 中国石油化工股份有限公司 | 一种地震波场吸收补偿成像方法及系统 |
CN106443786A (zh) * | 2016-11-14 | 2017-02-22 | 中国科学院地质与地球物理研究所 | 基于地面接收的反射地震资料的q值场建模方法 |
CN106443786B (zh) * | 2016-11-14 | 2018-04-20 | 中国科学院地质与地球物理研究所 | 基于地面接收的反射地震资料的q值场建模方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102230973A (zh) | 一种三维分步傅立叶粘滞声波深度偏移方法 | |
CN102478667B (zh) | 频散avo反演方法 | |
CN102112894A (zh) | 用地震表面波的波形评估土壤性质 | |
CN104155693A (zh) | 储层流体流度的角道集地震响应数值计算方法 | |
CN102156297B (zh) | 基于砂岩油藏叠后地震数据的流体替换方法 | |
CN104570072A (zh) | 一种粘弹性介质中的球面pp波反射系数建模方法 | |
CN103412327A (zh) | 一种裂缝性储层的粘弹性参数提取方法 | |
CN101598803A (zh) | 一种直接得到转换波叠加剖面的方法 | |
Chen et al. | Numerical simulation of frequency-dependent seismic response and gas reservoir delineation in turbidites: A case study from China | |
CN103487831A (zh) | Avo地震正演计算方法 | |
Lin et al. | Effect of lateral heterogeneity on surface wave testing: Numerical simulations and a countermeasure | |
Luo et al. | Generation of a pseudo-2D shear-wave velocity section by inversion of a series of 1D dispersion curves | |
Zhang et al. | Retrieval of shallow S-wave profiles from seismic reflection surveying and traffic-induced noise | |
CN102998703A (zh) | 基于地表一致性反褶积进行储层预测的方法及设备 | |
CN103364834A (zh) | 一种利用叠前地震频散分析预测储层渗透率的方法 | |
CN104007465B (zh) | 基于弹性反演的纯纵波数据刻画储层发育细节方法 | |
CN104714253A (zh) | 一种基于弥散黏滞性波动方程的avo/ava分析方法 | |
Lin et al. | Assessment of ground improvement with improved columns by surface wave testing | |
Tremblay et al. | Practical considerations for array-based surface-wave testing methods with respect to near-field effects and shear-wave velocity profiles | |
Zhu | Seismic modeling, inversion, and imaging in attenuating media | |
Wang et al. | A new insight into the hydrocarbon indicator of tight sandstone reservoir | |
Nzikou Mamboukou | Forward modelling and inversion of the ultrasonic wave propagation through a homogeneous and porous rock | |
Hailemikael et al. | From ambient vibration data analysis to 1D ground-motion prediction of the Mj 5.9 and the Mj 6.5 Kumamoto earthquakes in the Kumamoto alluvial plain, Japan | |
Joh et al. | Wavenumber unwrapping for the automated analysis of SASW measurements | |
Zhang et al. | Simulated various characteristic waves in acoustic full waveform relating to cement bond on the secondary interface |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20111102 |