CN113109875A - 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法 - Google Patents

一种全波形速度场约束下的盐下碳酸盐岩储层反演方法 Download PDF

Info

Publication number
CN113109875A
CN113109875A CN202110483109.5A CN202110483109A CN113109875A CN 113109875 A CN113109875 A CN 113109875A CN 202110483109 A CN202110483109 A CN 202110483109A CN 113109875 A CN113109875 A CN 113109875A
Authority
CN
China
Prior art keywords
inversion
low
wave impedance
full waveform
longitudinal wave
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
CN202110483109.5A
Other languages
English (en)
Other versions
CN113109875B (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.)
Beijing Research Center of CNOOC China Ltd
CNOOC China Ltd
Original Assignee
Beijing Research Center of CNOOC China Ltd
CNOOC China 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 Beijing Research Center of CNOOC China Ltd, CNOOC China Ltd filed Critical Beijing Research Center of CNOOC China Ltd
Priority to CN202110483109.5A priority Critical patent/CN113109875B/zh
Publication of CN113109875A publication Critical patent/CN113109875A/zh
Application granted granted Critical
Publication of CN113109875B publication Critical patent/CN113109875B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • G01V1/50Analysing data

Landscapes

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

Abstract

本发明涉及一种全波形速度场约束下的盐下碳酸盐岩储层反演方法,其包括以下步骤:1)目标工区岩石物理分析及井震标定,2)全波形速度场约束下低频模型的构建,3)基于部分叠加地震数据进行叠前同时反演,4)反演结果检验及分析。本发明在求取低频模型时,采用全波形速度场约束下的多井插值法进行低频模型构建,准确补充地震所缺低频成分;同时基于岩石物理分析,确定IP‑IS参数(纵波阻抗减横波阻抗)为反演敏感弹性参数;通过叠前反演,得到IP‑IS弹性参数体,进一步提升盐下碳酸盐岩储层预测精度。

Description

一种全波形速度场约束下的盐下碳酸盐岩储层反演方法
技术领域
本发明涉及地质领域,具体的说,涉及一种全波形速度场约束下的盐下碳酸盐岩储层反演方法。
背景技术
目前,反演已广泛应用于勘探阶段的储层定性或定量预测和开发阶段的储量计算、井网部署、油藏动态监测等方面,已成为地下储层刻画过程中必不可少的方法之一。反演过程中,通过低频模型补偿地震资料中缺失的低频信息是开展岩性识别、储层定量解释的重要步骤。一般情况下,低频模型是以地震解释层位和沉积规律为约束,将测井信息在整个数据体范围内进行内插和外推得到的,一定程度上补偿了地震资料中缺失的低频信息。但是,在井间储层横向变化剧烈、井点位置地震反射能量差异较大的区域内,补充的低频信息反而会给反演结果定量解释带来误差。
碳酸盐岩储层在全球油气分布中占有重要地位,碳酸盐岩储集层构成的油气田储量大、单井产量高,容易形成大型油气田。和碎屑岩储层相比,碳酸盐岩储层更为复杂,空间分布随机性强,突变性大,储层纵向巨厚,多期发育,横向沉积相带变化快,非均质性强,储层与围岩阻抗差异不明显。同时,盐丘和围岩间存在巨大的速度差异,造成盐下地震波场十分复杂,地震资料的分辨率和信噪比低,导致盐下碳酸盐岩储层储层识别和预测难度大。
发明内容
针对上述问题,本发明的目的是提供一种全波形速度场约束下的盐下碳酸盐岩储层反演方法,基于全波形反演得到的速度场约束建立低频模型,补充地震所缺低频成分;同时基于岩石物理分析,确定IP-IS参数(纵波阻抗-减横波阻抗参数)为反演敏感弹性参数;通过叠前反演,得到IP-IS弹性参数体,实现盐下碳酸盐岩储层精细预测。
为实现上述目的,本发明提供一种全波形速度场约束下的盐下碳酸盐岩储层反演方法,对于目标工区基于全波形反演和岩石物理分析的结果进行叠前反演得到得到IP-IS弹性参数体(即纵波阻抗-减横波阻抗参数体),实现盐下碳酸盐岩储层精细预测。
所述的全波形速度场约束下的盐下碳酸盐岩储层反演方法,包括以下步骤:
1)目标工区岩石物理分析及井震标定;
2)全波形速度场约束下低频模型的构建,其中,全波形速度场是针对目标工区的特有的低频模型建立方法;
3)基于部分叠加地震数据对目标工区进行叠前反演;
4)将步骤3)得到的反演结果进一步处理,得到IP-IS弹性参数体,作为反演的最终结果,对储层进行精细预测。
所述步骤1)中还包括,在测井资料一致性处理及横波预测的基础上,进行岩石物理分析,明确能够区分储层及非储层的反演敏感弹性参数;根据测井数据、部分叠加地震数据进行井震标定,提取多井子波。
所述步骤2)中采用全波形速度场约束下的多井插值法进行低频模型构建,具体过程包括:
21)首先,速度场应为全波形反演得到的速度场,并对其进行去野值、平滑、插值处理;
22)提取目的层段测井纵波阻抗曲线,对其低通滤波,保留10hz以下有效成分,分别统计目的层段全波形速度以及低通滤波后测井纵波阻抗的均值和方差,利用正态分布参数估计公式,将全波形速度场转化为测井纵波阻抗数据体;
23)应用简单克里格法消除井点位置的纵波阻抗预测误差,得到纵波阻抗低频模型;
24)分别对目的层段测井纵波阻抗与密度,纵波阻抗与横波阻抗进行交会分析,并进行线性回归,利用线性回归公式将步骤23)得到的纵波阻抗低频模型转换为横波阻抗低频模型和密度低频模型,做为叠前地震反演所需的低频模型。
所述步骤22)中,正态分布参数估计公式如下:
result=input*(a/b)+(c-d)*(a/b) (1)
式中,a为低通滤波后测井纵波阻抗的标准差,b为目的层段全波形速度的标准差,c为低通滤波后测井纵波阻抗的均值,d为目的层段全波形速度的均值。
所述步骤3)中基于部分叠加地震数据进行叠前同时反演具体包括:
31)基于部分叠加地震数据体,采用Aki-Richard近似公式进行叠前同时反演计算,得到缺乏低频成分的弹性参数体,包括纵波阻抗、横波阻抗和密度;
32)采用频率域合并的方式,将步骤2)得到的低频模型频带范围内的值,合并到步骤31)得到的反演结果中,补充原始地震缺失的低频成分。
其中,述补充后的新得到的反演结果包括纵波阻抗、横波阻抗和密度。
所述步骤31)中,Aki-Richard近似公式如下:
Figure BDA0003049172810000021
式中,Rpp(θ)为纵波的反射系数,vp和vs分别为上层介质的纵、横波速度,Δvp和Δvs分别为上、下层介质的纵、横波速度差,θ为地震波的入射角。
所述步骤4)中将步骤32)得到的反演结果进一步处理具体为,将步骤3)得到的补充低频成分后的纵波阻抗减去横波阻抗,得到IP-IS弹性参数体,作为反演的最终结果。
本发明由于采取以上技术方案,其具有以下优点:本发明在求取低频模型时,采用全波形速度场约束下的多井插值法进行低频模型构建,准确补充地震所缺低频成分;同时基于岩石物理分析,确定IP-IS参数(纵波阻抗减横波阻抗)为反演敏感弹性参数;通过叠前反演,得到IP-IS弹性参数体,进一步提升盐下碳酸盐岩储层预测精度。
附图说明
图1是本发明叠前反演流程示意图;
图2是盐下碳酸盐岩地震剖面图(局部盐岩厚度超过2400m,地震信号能量衰减严重);
图3(a)是全波形速度剖面图,图3(b)是原始地震剖面图,图3(c)是全波形速度储层段切片属性图;
图4(a)是目标工区岩石物理交会分析图,图4(b)是IP-IS参数直方图;
图5(a)是井震标定图,图5(b)是提取的多井子波示意图;
图6(a)是测井纵波阻抗与横波阻抗交会分析图,图6(b)是测井纵波阻抗与密度交汇分析图;
图7利用不同方法建立的低频模型平面属性图,(a)采用多井反距离加权插值方法建立的低频模型平面属性图,(b)采用全波形速度场约束的多井插值方法建立的低频模型平面属性图;
图8为全波形速度场约束下的盐下碳酸盐岩储层反演盲井连井剖面。
具体实施方式
以下结合附图对本发明的原理和特征进行描述,所举实例只用于解释本发明,并非用于限定本发明的范围。
本发明提供一种全波形速度场约束下的盐下碳酸盐岩储层反演方法,对于目标工区基于全波形反演和岩石物理分析的结果进行叠前反演得到得到IP-IS弹性参数体(即纵波阻抗-减横波阻抗参数体),实现盐下碳酸盐岩储层精细预测。包括以下步骤:
1)目标工区岩石物理分析及井震标定:
在测井资料一致性处理及横波预测的基础上,进行岩石物理分析,明确能够区分储层及非储层的反演敏感弹性参数;
根据测井数据、部分叠加地震数据进行井震标定,提取多井子波;
2)全波形速度场约束下低频模型的构建:
采用全波形速度场约束下的多井插值法进行低频模型构建,具体过程包括:
①首先,速度场应为全波形反演得到的速度场,并对其进行去野值、平滑、插值处理;
②提取目的层段测井纵波阻抗曲线,对其低通滤波,保留10hz以下有效成分,分别统计目的层段全波形速度以及低通滤波后测井纵波阻抗的均值和方差,利用正态分布参数估计公式,将全波形速度场转化为测井纵波阻抗数据体。
③应用简单克里格法消除井点位置的纵波阻抗预测误差,得到纵波阻抗低频模型。
④分别对目的层段测井纵波阻抗与密度,纵波阻抗与横波阻抗进行交会分析,并进行线性回归,利用线性回归公式将步骤③得到的纵波阻抗低频模型转换为横波阻抗低频模型和密度低频模型,做为叠前地震反演所需的低频模型;
3)基于部分叠加地震数据进行叠前同时反演:
基于部分叠加地震数据体,采用Aki-Richard近似公式进行叠前同时反演计算,得到缺乏低频成分的弹性参数体,包括纵波阻抗、横波阻抗和密度;
同时,采用频率域合并的方式,将步骤2)得到的低频模型频带范围内的值,合并到步骤3)得到的反演结果中,补充原始地震缺失的低频成分;
4)反演结果检验及分析:
将步骤3)得到的补充低频成分后的纵波阻抗减去横波阻抗,得到IP-IS弹性参数体,作为反演的最终结果,对储层进行精细预测。
优选地,在所述步骤2)中,正态分布参数估计公式如下:
result=input*(a/b)+(c-d)*(a/b) (1)
式中,a为低通滤波后测井纵波阻抗的标准差,b为目的层段全波形速度的标准差,c为低通滤波后测井纵波阻抗的均值,d为目的层段全波形速度的均值。
优选地,在所述步骤3)中,Aki-Richard近似公式如下:
Figure BDA0003049172810000041
式中,Rpp(θ)为纵波的反射系数,vp和vs分别为上层介质的纵、横波速度,Δvp和Δvs分别为上、下层介质的纵、横波速度差,θ为地震波的入射角。
如图1所示,本发明提供一种全波形速度场约束下的盐下碳酸盐岩储层反演方法。
图2是盐下碳酸盐岩地震剖面图,从图中可以看出,目的储层上方存在巨厚盐岩层,局部厚度超过2400m,盐丘和围岩间存在巨大的速度差异,造成盐下地震信号能量衰减严重,地震波场复杂,资料的分辨率和信噪比大幅降低。
图3(a)是全波形速度剖面图,图3(b)是对应原始地震剖面图,图3(c)是全波行速度储层段切片属性图,从图中可以看出全波形速度模型很好地反映了碳酸盐岩储层分布特征。相对传统地震叠加速度频带集中在0-2hz,全波形速度频带可达0-5hz,能更加充分补充原始地震缺失的低频成分,更准确刻画厚储层的边界及空间展布。
具体实施步骤:
1)目标工区岩石物理分析及井震标定:
在测井资料一致性处理及横波预测的基础上,进行岩石物理分析,明确能够区分储层及非储层的反演敏感弹性参数,从交汇图4(a)及直方图4(b)中可以看出,IP-IS参数是能够区分储层及非储层的反演敏感弹性参数;根据测井数据、部分叠加地震数据进行井震标定(图5(a)),提取稳定的多井子波(图5(b));
2)全波形速度场约束下低频模型的构建:
采用全波形速度场约束下的多井插值法进行低频模型构建,具体过程包括:
①首先,速度场应为全波形反演得到的速度场,并对其进行去野值、平滑、插值处理;
②提取目的层段测井纵波阻抗曲线,对其低通滤波,保留10hz以下有效成分,分别统计目的层段全波形速度以及低通滤波后测井纵波阻抗的均值和方差,利用正态分布参数估计公式,将全波形速度场转化为测井纵波阻抗数据体。
③应用简单克里格法消除井点位置的纵波阻抗预测误差,得到纵波阻抗低频模型。
④分别对目的层段测井纵波阻抗与密度(图6(a)),纵波阻抗与横波阻抗进行交会分析(图6(b)),并进行线性回归,利用线性回归公式将步骤③得到的纵波阻抗低频模型转换为横波阻抗低频模型和密度低频模型,做为叠前地震反演所需的低频模型。图7(b)是采用全波形速度场约束的多井插值方法建立的低频模型平面属性图,相比图7(a)所示采用常规反距离加权插值方法建立的低频模型平面属性,采用全波形速度场约束低频模型更能反映碳酸盐岩储层展布特征;
3)基于部分叠加地震数据进行叠前同时反演:
基于部分叠加地震数据体,采用Aki-Richard近似公式进行叠前同时反演计算,得到缺乏低频成分的弹性参数体,包括纵波阻抗、横波阻抗和密度;
同时,采用频率域合并的方式,将步骤2)得到的低频模型频带范围内的值,合并到步骤3)得到的反演结果中,补充原始地震缺失的低频成分;
4)反演结果检验及分析:
将步骤3)得到的补充低频成分后的纵波阻抗减去横波阻抗,得到IP-IS弹性参数体,作为反演的最终结果,对储层进行精细预测。反演结果如图8所示,从图中可以看出,反演结果与井吻合度高,井间预测性较强。
优选地,在所述步骤2)中,正态分布参数估计公式如下:
result=input*(a/b)+(c-d)*(a/b) (1)
式中,a为低通滤波后测井纵波阻抗的标准差,b为目的层段全波形速度的标准差,c为低通滤波后测井纵波阻抗的均值,d为目的层段全波形速度的均值。
优选地,在所述步骤3)中,Aki-Richard近似公式如下:
Figure BDA0003049172810000061
式中,Rpp(θ)为纵波的反射系数,vp和vs分别为上层介质的纵、横波速度,Δvp和Δvs分别为上、下层介质的纵、横波速度差,θ为地震波的入射角。
本发明仅以上述实施例进行说明,各部件的结构、设置位置及其连接都是可以有所变化的。在本发明技术方案的基础上,凡根据本发明原理对个别部件进行的改进或等同变换,均不应排除在本发明的保护范围之外。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (8)

1.一种全波形速度场约束下的盐下碳酸盐岩储层反演方法,其特征在于,对于目标工区基于全波形反演和岩石物理分析的结果进行叠前反演得到得到IP-IS弹性参数体,实现盐下碳酸盐岩储层精细预测。
2.根据权利要求1所述的全波形速度场约束下的盐下碳酸盐岩储层反演方法,其特征在于,包括以下步骤:
1)目标工区岩石物理分析及井震标定;
2)全波形速度场约束下的目标工区低频模型的构建;
3)基于部分叠加地震数据进行叠前同时反演;
4)将步骤3)得到的反演结果进一步处理,得到IP-IS弹性参数体,作为反演的最终结果,对储层进行精细预测。
3.根据权利要求2所述的全波形速度场约束下的盐下碳酸盐岩储层反演方法,其特征在于,所述步骤1)中还包括,在测井资料一致性处理及横波预测的基础上,进行岩石物理分析,明确能够区分储层及非储层的反演敏感弹性参数;根据测井数据、部分叠加地震数据进行井震标定,提取多井子波。
4.根据权利要求2所述的全波形速度场约束下的盐下碳酸盐岩储层反演方法,其特征在于,所述步骤2)中采用全波形速度场约束下的多井插值法进行低频模型构建,具体过程包括:
21)首先,速度场应为全波形反演得到的速度场,并对其进行去野值、平滑、插值处理;
22)提取目的层段测井纵波阻抗曲线,对其低通滤波,保留10hz以下有效成分,分别统计目的层段全波形速度以及低通滤波后测井纵波阻抗的均值和方差,利用正态分布参数估计公式,将全波形速度场转化为测井纵波阻抗数据体;
23)应用简单克里格法消除井点位置的纵波阻抗预测误差,得到纵波阻抗低频模型;
24)分别对目的层段测井纵波阻抗与密度,纵波阻抗与横波阻抗进行交会分析,并进行线性回归,利用线性回归公式将步骤23)得到的纵波阻抗低频模型转换为横波阻抗低频模型和密度低频模型,做为叠前地震反演所需的低频模型。
5.根据权利要求4所述的全波形速度场约束下的盐下碳酸盐岩储层反演方法,其特征在于,所述步骤22)中,正态分布参数估计公式如下:
result=input*(a/b)+(c-d)*(a/b) (1)
式中,a为低通滤波后测井纵波阻抗的标准差,b为目的层段全波形速度的标准差,c为低通滤波后测井纵波阻抗的均值,d为目的层段全波形速度的均值。
6.根据权利要求2所述的全波形速度场约束下的盐下碳酸盐岩储层反演方法,其特征在于,所述步骤3)中基于部分叠加地震数据进行叠前同时反演具体包括:
31)基于部分叠加地震数据体,采用Aki-Richard近似公式进行叠前同时反演计算,得到缺乏低频成分的弹性参数体,包括纵波阻抗、横波阻抗和密度;
32)采用频率域合并的方式,将步骤2)得到的低频模型频带范围内的值,合并到步骤31)得到的反演结果中,补充原始地震缺失的低频成分。
7.根据权利要求6所述的全波形速度场约束下的盐下碳酸盐岩储层反演方法,其特征在于,在所述步骤31)中,Aki-Richard近似公式如下:
Figure FDA0003049172800000021
式中,Rpp(θ)为纵波的反射系数,vp和vs分别为上层介质的纵、横波速度,Δvp和Δvs分别为上、下层介质的纵、横波速度差,θ为地震波的入射角。
8.根据权利要求6或7所述的全波形速度场约束下的盐下碳酸盐岩储层反演方法,其特征在于,所述步骤4)中将步骤32)得到的反演结果进一步处理具体为,将步骤3)得到的补充低频成分后的纵波阻抗减去横波阻抗,得到IP-IS弹性参数体,作为反演的最终结果。
CN202110483109.5A 2021-04-30 2021-04-30 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法 Active CN113109875B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110483109.5A CN113109875B (zh) 2021-04-30 2021-04-30 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110483109.5A CN113109875B (zh) 2021-04-30 2021-04-30 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法

Publications (2)

Publication Number Publication Date
CN113109875A true CN113109875A (zh) 2021-07-13
CN113109875B CN113109875B (zh) 2022-12-20

Family

ID=76720805

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110483109.5A Active CN113109875B (zh) 2021-04-30 2021-04-30 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法

Country Status (1)

Country Link
CN (1) CN113109875B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117236172A (zh) * 2023-09-15 2023-12-15 中海石油(中国)有限公司深圳分公司 一种针对火成岩的深度学习低频模型构建方法及系统
WO2024103475A1 (zh) * 2022-11-14 2024-05-23 中国石油天然气股份有限公司 一种碳酸盐岩油气藏定量表征方法与系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104570128A (zh) * 2013-10-16 2015-04-29 中国石油化工股份有限公司 一种基于岩石脆性因子的非常规油气储层预测方法
CN110609327A (zh) * 2018-06-14 2019-12-24 中国石油化工股份有限公司 基于叠前地震属性的碳酸盐岩储集相预测方法及装置
WO2020146863A1 (en) * 2019-01-13 2020-07-16 Schlumberger Technology Corporation Seismic image data interpretation system
CN111596354A (zh) * 2020-05-11 2020-08-28 同济大学 一种机器学习框架下考虑空间约束的地震储层预测方法
CN111897011A (zh) * 2020-07-27 2020-11-06 中国石油大学(华东) 储层孔隙特征确定方法、装置及设备
CN112114380A (zh) * 2020-08-11 2020-12-22 中国石油化工股份有限公司 一种用于砂砾岩油藏储层的预测方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104570128A (zh) * 2013-10-16 2015-04-29 中国石油化工股份有限公司 一种基于岩石脆性因子的非常规油气储层预测方法
CN110609327A (zh) * 2018-06-14 2019-12-24 中国石油化工股份有限公司 基于叠前地震属性的碳酸盐岩储集相预测方法及装置
WO2020146863A1 (en) * 2019-01-13 2020-07-16 Schlumberger Technology Corporation Seismic image data interpretation system
CN111596354A (zh) * 2020-05-11 2020-08-28 同济大学 一种机器学习框架下考虑空间约束的地震储层预测方法
CN111897011A (zh) * 2020-07-27 2020-11-06 中国石油大学(华东) 储层孔隙特征确定方法、装置及设备
CN112114380A (zh) * 2020-08-11 2020-12-22 中国石油化工股份有限公司 一种用于砂砾岩油藏储层的预测方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
张晶玉 等: "多属性储层参数反演方法在M油田中的应用", 《油气藏评价与开发》 *
董洪超 等: "基于精细地层格架的古近系储层相控反演技术研究", 《物探化探计算技术》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2024103475A1 (zh) * 2022-11-14 2024-05-23 中国石油天然气股份有限公司 一种碳酸盐岩油气藏定量表征方法与系统
CN117236172A (zh) * 2023-09-15 2023-12-15 中海石油(中国)有限公司深圳分公司 一种针对火成岩的深度学习低频模型构建方法及系统

Also Published As

Publication number Publication date
CN113109875B (zh) 2022-12-20

Similar Documents

Publication Publication Date Title
CA3076280C (en) Method and apparatus for determining thin interlayers
CN109425896B (zh) 白云岩油气储层分布预测方法及装置
CN111596364B (zh) 基于高精度层序地层格架的地震沉积微相组合分析方法
CN104570067B (zh) 一种地球物理勘探中相控地震反演方法
CN110231652B (zh) 一种基于密度的含噪声应用空间聚类的地震相提取方法
CN109521474B (zh) 一种三维双控下的叠前地质统计学反演方法
CN114994759A (zh) 基于gan网络的碳封存箱智能识别方法和系统
CN104769458A (zh) 一种基于柯西分布的叠后波阻抗反演方法
CN113109875B (zh) 一种全波形速度场约束下的盐下碳酸盐岩储层反演方法
CN111506861B (zh) 一种目的层有利区域裂缝强度计算方法
CN114994758B (zh) 碳酸盐岩断控储层的波阻抗提取与结构表征方法和系统
CN111856566B (zh) 湖相滩坝砂体中薄储层预测方法及装置
CN111983683B (zh) 一种少井条件下盐下湖相灰岩储层预测方法和系统
CN108508489B (zh) 一种基于波形微变化匹配的地震反演方法
CN112114380A (zh) 一种用于砂砾岩油藏储层的预测方法
CN110727027B (zh) 多期次河道砂的精细刻画方法
CN115857047A (zh) 一种地震储层综合预测方法
CN113848593A (zh) 一种定量预测含煤地层中岩浆岩侵蚀区的方法
CN110795513B (zh) 河流相源储异位型致密油气甜点区分布的预测方法
CN113806674A (zh) 古河道纵向尺度的量化方法、装置、电子设备及存储介质
CN112147676A (zh) 一种煤层及夹矸厚度预测方法
CN112748467B (zh) 利用地震资料确定Russell流体因子的方法及系统
CN110850504B (zh) 一种基于铀曲线拟阻抗约束的页岩密度参数叠前反演方法
CN113640873B (zh) 砂岩型铀矿叠前地震预测方法及装置
CN111352154B (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