CN112394391B - 雷克子波分数阶导数的地震相定量表征方法 - Google Patents

雷克子波分数阶导数的地震相定量表征方法 Download PDF

Info

Publication number
CN112394391B
CN112394391B CN201911335706.2A CN201911335706A CN112394391B CN 112394391 B CN112394391 B CN 112394391B CN 201911335706 A CN201911335706 A CN 201911335706A CN 112394391 B CN112394391 B CN 112394391B
Authority
CN
China
Prior art keywords
seismic
thin interbed
fractional
reservoir
waveform
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
Application number
CN201911335706.2A
Other languages
English (en)
Other versions
CN112394391A (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 Offshore Oil Corp CNOOC
CNOOC China Ltd Tianjin Branch
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Tianjin Branch
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 China National Offshore Oil Corp CNOOC, CNOOC China Ltd Tianjin Branch filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201911335706.2A priority Critical patent/CN112394391B/zh
Publication of CN112394391A publication Critical patent/CN112394391A/zh
Application granted granted Critical
Publication of CN112394391B publication Critical patent/CN112394391B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/24Recording seismic data
    • G01V1/245Amplitude control for seismic recording
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/282Application of seismic models, synthetic seismograms

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种雷克子波分数阶导数的地震相定量表征方法,包括:一,采用高精度合成记录方法,对薄互层储层精细标定,确定薄互层储层对应的地震响应;二,对比薄互层储层对应的地震响应与替换后地震响应,分析薄互层储层对地震同相轴振幅响应的主控因素;三,将得到的薄互层地震响应作为单波叠合响应;四,利用雷克子波多种分数阶导数波形,表征薄互层储层对应地震反射同相轴的振幅、频率和波形信息,代表不同的薄互层储层结构;五,利用匹配追踪算法,将不同子波分数阶导数波形与地震层析约束下的地震波形进行匹配处理,找到地震波形所对应的分数阶导数,以表征地震道对应的振幅、波形和频率,完成薄互层储层地震相定性半定量描述。

Description

雷克子波分数阶导数的地震相定量表征方法
技术领域
本发明涉及石油勘探地震数据处理与解释技术领域,尤其涉及一种雷克子波分数阶导数的地震相定量表征方法。
背景技术
随着勘探开发进程的不断推进,勘探开发目标逐渐转向包括薄互层储层油藏、岩性油藏等在内的复杂油气藏。其中,新近系馆陶组薄互层储层因埋藏浅、分布范围广,逐渐成为渤海勘探的重点。渤海庙西探区馆陶组储层以薄互层储层为主,储层单层厚度薄,5m以下砂体占比超过80%,泥岩盖层薄且发育不稳定,导致薄互层之间调谐干涉严重,现有地震资料表现为杂乱弱连续反射或空白反射,难以开展单砂体尺度的储层预测工作。研究表明,薄互层储层发育程度与其对应地震波的波形、频率和振幅有密切关系。因此,寻求能够表征薄互层储层发育程度的地震相地球物理表征参数,定性或半定量预测薄互层储层发育程度及平面展布,成为渤海庙西探区勘探过程中的一项重要课题。
发明内容
本发明的主要目的在于克服现有技术存在的上述缺点,而提供一种雷克子波分数阶导数的地震相定量表征方法,其根据地震地层学,利用已钻井资料所揭示的沉积旋回特征标定地震资料所揭示的时频旋回特征,利用匹配追踪算法,将薄互层的地震波形与反射子波的分数阶导数进行匹配,找到地震波对应的分数阶导数,达到薄互层储层地震相定性半定量描述的目的。
本发明的目的是由以下技术方案实现的。
本发明雷克子波分数阶导数的地震相定量表征方法,包括以下步骤:
第一步,采用高精度合成记录方法,对薄互层储层进行精细标定,确定薄互层储层对应的地震响应;通过岩石物理参数替换,将薄互层储层中的砂岩逐一替换为泥岩,分析替换后对应的地震响应特征;
第二步,在第一步的基础上,对比薄互层储层对应的地震响应与替换后的地震响应,分析薄互层储层对地震同相轴振幅响应的主控因素;
第三步,根据单波具有多种频率和波形的性质,将第二步得到的薄互层地震响应作为一个或多个单波的叠合响应;
第四步,利用雷克子波的多种分数阶导数波形,拟合不同薄互层地震响应对应的波形,形成不同地层结构的薄互层储层地震集合,使多种分数阶导数能够表征薄互层储层对应地震反射同相轴的振幅、频率和波形信息,从而代表不同的薄互层储层结构;
第五步,利用匹配追踪算法,将不同子波的分数阶导数波形与地震层析约束下的地震波形进行匹配处理,找到地震波形所对应的分数阶导数,得到地震相分数阶导数分布图,以此来表征地震道所对应的振幅、波形和频率,完成薄互层储层地震相定性半定量描述。
前述的雷克子波分数阶导数的地震相定量表征方法,其中,所述
第一步,通过计算机软件,制作高精度合成地震记录,确定薄互层储层的地震响应特征;时深关系为地震波反射时间与地层海拔深度的对应关系,在得到准确的时深关系的基础上,通过计算机软件,以薄互层储层砂层组为储层研究单位,按照由上而下的顺序以及泥岩趋势线,逐个将砂岩的声波时差曲线和密度曲线替换为泥岩的曲线特征,观察替换后的合成记录变化情况;
第二步,分析薄互层储层砂层组的地震响应特征,得到薄互层储层每一个砂体对合成记录波形的贡献度,从而得到薄互层储层地震响应的主控因素;将地震波波形看做是每一个薄互层储层反射的叠加效果,每个砂体对地震波形均有一定的贡献度;
第三步,根据单波具有多种频率和波形的性质,将第二步得到的薄互层地震响应作为一个或多个单波的叠合响应;
第四步,根据第二步和第三步的认识,计算不同频率的雷克子波的分数阶导数波形,建立不同频率的雷克子波分数阶导数波形和薄互层储层地震波形的关系,从而拟合不同薄互层地震响应对应的波形,形成不同地层结构的薄互层储层地震集合,使多种分数阶导数能够表征薄互层储层对应地震反射同相轴的振幅、频率和波形信息,从而代表不同的薄互层储层结构;
第五步,根据匹配追踪算法,通过计算机程序,利用最小二乘法将不同子波的分数阶导数波形与地震层析约束下的地震波形进行匹配处理,找到与薄互层储层地震响应波形最为相近的雷克子波的分数阶导数,以此来表征地震道所对应的振幅、波形和频率,完成薄互层储层地震相定性半定量描述。
前述的基于雷克子波分数阶导数的地震相定量表征方法,其中:所述第五步的具体步骤包括,
(1)在薄互层储层地震反射特征精细标定的基础上,利用地震同相轴的振幅、波形和频率属性反映薄互层储层的发育情况,通过岩石物理分析及测井曲线岩性替换,对井点处薄互层储层中的砂岩逐一替换成泥岩,分析不同地层结构的薄互层储层的地震响应特征;
(2)参考地层沉积的旋回性,将砂泥岩突变界面看做岩性过渡带,这种过渡带对应的波形是由几个具有零阶和一阶分数阶导数的子波合成的,因而可以假设一个地层分界面阻抗的数学模型是由任意阶分数阶导数的因果信号和非因果信号组成;其中,因果信号和非因果信号的数学表达式如下:
Figure BDA0002330865420000031
Figure BDA0002330865420000032
其中Γ为Gamma函数,其表达式如下所示,
Figure BDA0002330865420000033
参数深度z和分数阶α生成了因果信号
Figure BDA0002330865420000034
非因果信号
Figure BDA0002330865420000036
和均衡信号
Figure BDA0002330865420000035
当α=0时,因果信号变为不连续的突变信号,此时表现为传统的两种介质的波阻抗突变,当α≥1时,因果信号是可微分的;当0<α<1时,因果信号是连续但不可微分的;当-1<α≤0时,因果信号既不连续又不可微分;
一维地震数据模型f可用分数阶样条函数表达如下:
Figure BDA0002330865420000041
其中ai为权系数,zi为分数阶导数为αi的位置。对f求导数,即变为模型的的反射系数:
Figure BDA0002330865420000042
众所周知,合成地震记录是反射系数与子波的褶积,即:
Figure BDA0002330865420000043
其中,
Figure BDA0002330865420000044
数学上,这种表达式称为子波
Figure BDA0002330865420000045
的分数阶导数,函数f的分数阶导数的定义如下:
Figure BDA0002330865420000046
由(1)式和(6)式可知,合成地震记录可以由子波的各分数阶导数加权而成,即:
Figure BDA0002330865420000047
(3)利用匹配追踪算法,将不同子波的分数阶导数波形与地震层析约束下的地震波形进行匹配处理,找到地震波形所对应的分数阶导数,以此来表征地震道所对应的波形和频率,达到薄互层储层半定量地震描述的目的。
本发明雷克子波分数阶导数的地震相定量表征方法的有益效果,根据地震地层学,利用已钻井资料所揭示的沉积旋回特征标定地震资料所揭示的时频旋回特征,并以此确定薄互层的研究尺度。在此基础上,利用匹配追踪算法,将薄互层的地震波形与反射子波的分数阶导数进行匹配,找到地震波对应的分数阶导数,以此表征地震道所对应的振幅、波形和频率,达到薄互层储层地震相定性半定量描述的目的。
具体实施方式
本发明雷克子波分数阶导数的地震相定量表征方法,其包括以下步骤:
第一步,采用高精度合成记录方法,对薄互层储层进行精细标定,确定薄互层储层对应的地震响应;通过岩石物理参数替换,将薄互层储层中的砂岩逐一替换为泥岩,分析替换后对应的地震响应特征;通过计算机软件,制作高精度合成地震记录,确定薄互层储层的地震响应特征;时深关系为地震波反射时间与地层海拔深度的对应关系,在得到准确的时深关系的基础上,通过计算机软件,以薄互层储层砂层组为储层研究单位,按照由上而下的顺序以及泥岩趋势线,逐个将砂岩的声波时差曲线和密度曲线替换为泥岩的曲线特征,观察替换后的合成记录变化情况。
第二步,在第一步的基础上,对比薄互层储层对应的地震响应与替换后的地震响应,分析薄互层储层对地震同相轴振幅响应的主控因素;分析薄互层储层砂层组的地震响应特征,得到薄互层储层每一个砂体对合成记录波形的贡献度,从而得到薄互层储层地震响应的主控因素;将地震波波形看做是每一个薄互层储层反射的叠加效果,每个砂体对地震波形均有一定的贡献度。
第三步,根据单波具有多种频率和波形的性质,将第二步得到的薄互层地震响应作为一个或多个单波的叠合响应;
第四步,利用雷克子波的多种分数阶导数波形,拟合不同薄互层地震响应对应的波形,形成不同地层结构的薄互层储层地震集合,使多种分数阶导数能够表征薄互层储层对应地震反射同相轴的振幅、频率和波形信息,从而代表不同的薄互层储层结构;即根据第二步和第三步的认识,计算不同频率的雷克子波的分数阶导数波形,建立不同频率的雷克子波分数阶导数波形和薄互层储层地震波形的关系,从而拟合不同薄互层地震响应对应的波形,形成不同地层结构的薄互层储层地震集合,使多种分数阶导数能够表征薄互层储层对应地震反射同相轴的振幅、频率和波形信息,从而代表不同的薄互层储层结构;
第五步,利用匹配追踪算法,通过计算机程序,利用最小二乘法将不同子波的分数阶导数波形与地震层析约束下的地震波形进行匹配处理,找到与薄互层储层地震波形所对应的分数阶导数,得到地震相分数阶导数分布图,以此来表征地震道所对应的振幅、波形和频率,完成薄互层储层地震相定性半定量描述;具体步骤包括:
(1)在薄互层储层地震反射特征精细标定的基础上,利用地震同相轴的振幅、波形和频率属性反映薄互层储层的发育情况,通过岩石物理分析及测井曲线岩性替换,对井点处薄互层储层中的砂岩逐一替换成泥岩,分析不同地层结构的薄互层储层的地震响应特征;
(2)参考地层沉积的旋回性,将砂泥岩突变界面看做岩性过渡带,这种过渡带对应的波形是由几个具有零阶和一阶分数阶导数的子波合成的,因而可以假设一个地层分界面阻抗的数学模型是由任意阶分数阶导数的因果信号和非因果信号组成;其中,因果信号和非因果信号的数学表达式如下:
Figure BDA0002330865420000061
Figure BDA0002330865420000062
其中Γ为Gamma函数,其表达式如下所示,
Figure BDA0002330865420000063
参数深度z和分数阶α生成了因果信号
Figure BDA0002330865420000064
非因果信号
Figure BDA0002330865420000066
和均衡信号
Figure BDA0002330865420000065
当α=0时,因果信号变为不连续的突变信号,此时表现为传统的两种介质的波阻抗突变,当α≥1时,因果信号是可微分的;当0<α<1时,因果信号是连续但不可微分的;当-1<α≤0时,因果信号既不连续又不可微分;
一维地震数据模型f可用分数阶样条函数表达如下:
Figure BDA0002330865420000071
其中ai为权系数,zi为分数阶导数为αi的位置。对f求导数,即变为模型的的反射系数:
Figure BDA0002330865420000072
众所周知,合成地震记录是反射系数与子波的褶积,即:
Figure BDA0002330865420000073
其中,
Figure BDA0002330865420000074
数学上,这种表达式称为子波
Figure BDA0002330865420000075
的分数阶导数。函数f的分数阶导数的定义如下:
Figure BDA0002330865420000076
由(1)式和(6)式可知,合成地震记录可以由子波的各分数阶导数加权而成,即:
Figure BDA0002330865420000077
(3)利用匹配追踪算法,将不同子波的分数阶导数波形与地震层析约束下的地震波形进行匹配处理,找到地震波形所对应的分数阶导数,以此来表征地震道所对应的波形和频率,达到薄互层储层半定量地震描述的目的。
本发明利用雷克子波不同的分数阶导数波形表征已钻井所揭示的薄互层储层组合的波形特征,利用匹配追踪算法,将具有特定地质含义的波形进行平面匹配预测,达到薄互层储层地震相定性半定量描述的目的。本发明具有物理含义明确、精度高、运算快捷的特点。
以上所述,仅是本发明的较佳实施例而已,并非对本发明作任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。

Claims (3)

1.一种雷克子波分数阶导数的地震相定量表征方法,其特征在于,包括以下步骤:
第一步,采用高精度合成记录方法,对薄互层储层进行精细标定,确定薄互层储层对应的地震响应;通过岩石物理参数替换,将薄互层储层中的砂岩逐一替换为泥岩,分析替换后对应的地震响应特征;
第二步,在第一步的基础上,对比薄互层储层对应的地震响应与替换后的地震响应,分析薄互层储层地震同相轴振幅响应的主控因素;
第三步,根据单波具有多种频率和波形的性质,将第二步得到的薄互层地震响应作为一个或多个单波的叠合响应;
第四步,利用雷克子波的多种分数阶导数波形,拟合不同薄互层地震响应对应的波形,形成不同地层结构的薄互层储层地震集合,使多种分数阶导数能够表征薄互层储层对应地震反射同相轴的振幅、频率和波形信息,从而代表不同的薄互层储层结构;
第五步,利用匹配追踪算法,将不同子波的分数阶导数波形与地震层析约束下的地震波形进行匹配处理,找到地震波形所对应的分数阶导数,得到地震相分数阶导数分布图,以此来表征地震道所对应的振幅、波形和频率,完成薄互层储层地震相定量描述。
2.根据权利要求1所述的雷克子波分数阶导数的地震相定量表征方法,其特征在于,具体步骤如下:
第一步,通过计算机软件,制作高精度合成地震记录,确定薄互层储层的地震响应特征;时深关系为地震波反射时间与地层海拔深度的对应关系,在得到准确的时深关系的基础上,通过计算机软件,以薄互层储层砂层组为储层研究单位,按照由上而下的顺序以及泥岩趋势线,逐个将砂岩的声波时差曲线和密度曲线替换为泥岩的曲线特征,观察替换后的合成记录变化情况;
第二步,分析薄互层储层砂层组的地震响应特征,得到薄互层储层每一个砂体对合成记录波形的贡献度,从而得到薄互层储层地震响应的主控因素;将地震波波形看做是每一个薄互层储层反射的叠加效果,每个砂体对地震波形均有一定的贡献度;
第三步,根据单波具有多种频率和波形的性质,将薄互层地震响应作为一个或多个单波的叠合响应;
第四步,计算不同频率的雷克子波的分数阶导数波形,建立不同频率的雷克子波分数阶导数波形和薄互层储层地震波形的关系,从而拟合不同薄互层地震响应对应的波形,形成不同地层结构的薄互层储层地震集合,使多种分数阶导数能够表征薄互层储层对应地震反射同相轴的振幅、频率和波形信息,从而代表不同的薄互层储层结构;
第五步,根据匹配追踪算法,通过计算机程序,利用最小二乘法将不同子波的分数阶导数波形与地震层析约束下的地震波形进行匹配处理,找到与薄互层储层地震响应波形最为相近的雷克子波的分数阶导数,以此来表征地震道所对应的振幅、波形和频率,完成薄互层储层地震相定量描述。
3.根据权利要求2所述的雷克子波分数阶导数的地震相定量表征方法,其特征在于:薄互层储层地震相定量描述的具体步骤包括:
(1)在薄互层储层地震反射特征精细标定的基础上,利用地震同相轴的振幅、波形和频率属性反映薄互层储层的发育情况,通过岩石物理分析及测井曲线岩性替换,对井点处薄互层储层中的砂岩逐一替换成泥岩,分析不同地层结构的薄互层储层的地震响应特征;
(2)参考地层沉积的旋回性,将砂泥岩突变界面看做岩性过渡带,这种过渡带对应的波形是由几个具有零阶和一阶分数阶导数的子波合成的,因而可以假设一个地层分界面阻抗的数学模型是由任意阶分数阶导数的因果信号和非因果信号组成;其中,因果信号和非因果信号的数学表达式如下:
Figure FDA0003640274240000021
Figure FDA0003640274240000022
其中Γ为Gamma函数,其表达式如下所示,
Figure FDA0003640274240000023
其中,z表示深度,t表示时间,R表示实数集,深度z和分数阶导数α生成了因果信号
Figure FDA0003640274240000024
非因果信号
Figure FDA0003640274240000025
和均衡信号
Figure FDA0003640274240000026
当α=0时,因果信号变为不连续的突变信号,此时表现为传统的两种介质的波阻抗突变,当α≥1时,因果信号是可微分的;当0<α<1时,因果信号是连续但不可微分的;当-1<α≤0时,因果信号既不连续又不可微分;
一维地震数据模型f(z)可用分数阶样条函数表达如下:
Figure FDA0003640274240000027
其中ai为权系数,zi为分数阶导数为αi的深度,对z求导数,即变为模型的反射系数:
Figure FDA0003640274240000028
其中,
Figure FDA0003640274240000029
ΛC为正反射系数集合,ΛA为负反射系数集合,众所周知,合成地震记录是反射系数与子波的褶积,即:
Figure FDA0003640274240000031
其中*表示褶积运算符号,数学上,这种表达式称为子波φ(z)的分数阶导数,f的分数阶导数Dαf的定义如下:
Figure FDA0003640274240000032
由(1)式和(6)式可知,合成地震记录可以由子波的各分数阶导数加权而成,即:
Figure FDA0003640274240000033
(3)利用匹配追踪算法,将不同子波的分数阶导数波形与地震层析约束下的地震波形进行匹配处理,找到地震波形所对应的分数阶导数,以此来表征地震道所对应的波形和频率,达到薄互层储层定量地震描述的目的。
CN201911335706.2A 2019-12-23 2019-12-23 雷克子波分数阶导数的地震相定量表征方法 Active CN112394391B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911335706.2A CN112394391B (zh) 2019-12-23 2019-12-23 雷克子波分数阶导数的地震相定量表征方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911335706.2A CN112394391B (zh) 2019-12-23 2019-12-23 雷克子波分数阶导数的地震相定量表征方法

Publications (2)

Publication Number Publication Date
CN112394391A CN112394391A (zh) 2021-02-23
CN112394391B true CN112394391B (zh) 2022-08-02

Family

ID=74603754

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911335706.2A Active CN112394391B (zh) 2019-12-23 2019-12-23 雷克子波分数阶导数的地震相定量表征方法

Country Status (1)

Country Link
CN (1) CN112394391B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104062681A (zh) * 2013-03-22 2014-09-24 中国石油化工股份有限公司 一种基于分数阶导数的地震层位追踪预处理方法

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102928875B (zh) * 2012-11-07 2016-05-11 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于分数阶傅里叶域的子波提取方法
EP3076205B1 (en) * 2015-03-31 2023-06-14 CGG Services SAS Method for survey data processing compensating for visco-acoustic effects in tilted transverse isotropy reverse time migration
CN109143340B (zh) * 2018-08-20 2020-03-10 中国海洋石油集团有限公司 一种基于常q模型的粘弹介质地震波模拟方法及系统
CN110471104B (zh) * 2019-08-26 2021-03-16 电子科技大学 基于智能特征学习的叠后地震反射模式识别方法
CN110471105B (zh) * 2019-09-10 2021-05-04 中国海洋石油集团有限公司 一种层位约束下基于地震解析道的相干增强断层识别方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104062681A (zh) * 2013-03-22 2014-09-24 中国石油化工股份有限公司 一种基于分数阶导数的地震层位追踪预处理方法

Also Published As

Publication number Publication date
CN112394391A (zh) 2021-02-23

Similar Documents

Publication Publication Date Title
US8923093B2 (en) Determining the quality of a seismic inversion
EP1746443B1 (en) Method of estimating elastic parameters and rock composition of underground formations using seismic data
Yanhu et al. A method of seismic meme inversion and its application
Azevedo et al. Geostatistical rock physics AVA inversion
EP2420863A2 (en) Determining a Position of a Geological Layer Relative to a Wavelet Response in Seismic Data
US10310117B2 (en) Efficient seismic attribute gather generation with data synthesis and expectation method
CN111142156B (zh) 基于相控的地震强屏蔽时频信息提取及剥离方法
CN113031068B (zh) 一种基于反射系数精确式的基追踪叠前地震反演方法
Zhang et al. Viabilities of seismic ray impedance and elastic impedance for hydrocarbon-sand discrimination
CN111506861B (zh) 一种目的层有利区域裂缝强度计算方法
CN102393532A (zh) 地震信号反演方法
CN113109875B (zh) 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法
CN112394391B (zh) 雷克子波分数阶导数的地震相定量表征方法
CN109143335B (zh) 一种合成地震记录的制作方法、系统、介质及设备
CN109471203B (zh) 提高地震资料分辨率效果评判优选方法
Volkov et al. From qualitative interpretation to quantitative analysis: prediction of properties of geological bodies by using the spectral decomposition attribute–case study of Achimov turbidity system in West-Siberia
Jinshui et al. Model-data-driven seismic inversion method based on small sample data
CN110389381B (zh) 一种基于地震属性的砂体储层预测方法及装置
Trickett In search of the vibroseis first arrival
Li et al. A lithofacies prediction method based on stochastic simulation
CN113311482B (zh) 一种基于叠前谱反演优化的高分辨率中深层储层预测方法
CN113009579B (zh) 地震数据反演方法及装置
Wu et al. The Application of Stochastic Seismic Inversion Method of Geostatistical Inversion in Putaohua X Area
Kelly et al. The effect of improved, low-frequency bandwidth in full-waveform inversion for velocity
CN115343765A (zh) 绕射波成像数据能量的定量刻度方法、介质及设备

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant