CN108107482A - 一种少井条件下烃源岩多参数预测方法及系统 - Google Patents

一种少井条件下烃源岩多参数预测方法及系统 Download PDF

Info

Publication number
CN108107482A
CN108107482A CN201711315830.3A CN201711315830A CN108107482A CN 108107482 A CN108107482 A CN 108107482A CN 201711315830 A CN201711315830 A CN 201711315830A CN 108107482 A CN108107482 A CN 108107482A
Authority
CN
China
Prior art keywords
organic carbon
data
source rock
hydrocarbon source
target zone
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
CN201711315830.3A
Other languages
English (en)
Other versions
CN108107482B (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 Normal University
Original Assignee
Beijing Normal University
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 Normal University filed Critical Beijing Normal University
Priority to CN201711315830.3A priority Critical patent/CN108107482B/zh
Publication of CN108107482A publication Critical patent/CN108107482A/zh
Application granted granted Critical
Publication of CN108107482B publication Critical patent/CN108107482B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Abstract

本发明公开一种少井条件下烃源岩多参数预测方法及系统，所述方法包括：采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置；在目的层的位置上，从时间域的地震体数据中提取地震属性信息；根据地震属性信息进行属性融合，并结合所述测井数据确定烃源岩位置；对测井数据中的岩性数据进行统计，确定各深度段岩性数据所对应的沉积相模式；利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置；按照有机碳体数据所在位置进行分割，获得多个不同沉积相模式下的有机碳体数据；对多个不同沉积相模式下的有机碳体数据进行组合，获得有机碳的预测值，提高了预测有机碳的精确度。

Description

一种少井条件下烃源岩多参数预测方法及系统
技术领域
本发明涉及烃源岩预测技术领域，特别是涉及一种少井条件下烃源岩多参数预测方法及系统。
背景技术
在盆地勘探初期,预测烃源岩对烃源岩位置的识别、对烃源岩厚度的预测以及对烃源岩质量的评价具有至关重要的作用。在多井地区，地球化学分析方法多应用于烃源岩的预测与评价，但是在少井地区，尤其是海上勘探区，井少直接制约着烃源岩的预测。故研究学者采用井震结合实现对烃源岩进行预测，先根据烃源岩“低频-强振幅”的特征在地震剖面上肉眼确定烃源岩位置，在烃源岩位置内通过地球物理反演与测井结合的方法得到有机碳预测值，但是现有预测过程中，由于现有根据频率高低划分烃源发育位置时，只是单从地震剖面上人为观察，是不定量的科学研究，同时地球物理反演过程中没有考虑沉积相模式的不同，因此上述预测有机碳存在预测精度低的问题。另外，现有在烃源岩位置上通过计算砂地比得到烃源岩厚度时，同一个地区不同的沉积相模式下均采用同一个预测模型，因此预测出的烃源岩的厚度也存在预测精确度低的问题。
基于上述问题，如何提高烃源岩参数预测的精度成为本领域亟需解决的问题。
发明内容
本发明的目的是提供一种少井条件下烃源岩多参数预测方法及系统，以提高少井条件下烃源岩多参数预测的精确度。
为实现上述目的，本发明提供一种少井条件下烃源岩多参数预测方法，所述方法包括：
获取时间域的地震体数据和测井数据；
采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置；
在所述目的层的位置上，从时间域的地震体数据中提取地震属性信息；
根据地震属性信息进行属性融合，并结合所述测井数据确定烃源岩位置；
对所述测井数据中的岩性数据进行统计，确定各深度段岩性数据所对应的沉积相模式；
在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置；
按照所述有机碳体数据所在位置进行分割，获得多个不同沉积相模式下的有机碳体数据；
对多个不同沉积相模式下的有机碳体数据进行组合，获得有机碳的预测值。
可选的，所述采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置，具体包括：
从所述时间域的地震体数据中提取多个相邻的单道地震数据；
对各相邻的单道地震数据进行傅里叶变换，获得多个频率域的单道地震数据；
根据各所述频率域的单道地震数据确定多个低频信息所在位置；
根据各所述低频信息所在位置确定烃源岩所在的目的层的位置。
可选的，所述在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置步骤之前还包括：
根据所述测井数据中的声波测井曲线和深侧向电阻率测井曲线确定有机碳参数预测曲线；
从有机碳参数预测曲线中提取有机碳预测值；
判断有机碳预测值与有机碳预测值对应的有机碳实测值之差是否小于第一设定值，如果小于或等于第一设定值，则选取有机碳参数预测曲线，如果大于第一设定值，则重新确定有机碳参数预测曲线。
可选的，所述方法还包括：
利用所统计的岩性数据中砂岩和泥岩所在深度，确定不同沉积相模式下的砂岩声波时差压实模型；
根据所述声波时差压实模型确定各沉积相模式下泥岩速度和砂岩速度；
获取速度谱；
根据所述速度谱确定目的层的速度；
根据所述目的层的速度、所述泥岩速度和所述砂岩速度确定各沉积相模式的泥岩百分比；
根据所述烃源岩所在目的层的顶底面深度确定各沉积相模式下的目的层的地层厚度；
根据所述泥岩百分比和所述目的层的地层厚度确定各沉积相模式的烃源岩初始厚度；
对各所述沉积相模式的烃源岩初始厚度进行组合，获得目的层的烃源岩厚度。
可选的，对各相邻的单道地震数据进行傅里叶变换，获得多个频率域的单道地震数据的具体公式为：
其中，x(t)为所提取的单道地震数据，t代表时间，f代表频率，w(τ-t,f)为高斯窗口，τ为控制高斯窗口在时间轴即t轴位置的参数，X(t,f)为变换后频率域的单道地震数据。
可选的，所述根据所述测井数据中的声波测井曲线和深侧向电阻率测井曲线确定有机碳参数预测曲线的具体公式为：
其中，R为深侧向电阻率测井曲线值，R基线为深侧向电阻率曲线基线值，△t为声波时差测井曲线值，△t基线为声波时差测井曲线基线值，LOM为热变指数,TOCbackground为背景TOC含量,W(TOC)为有机碳参数预测曲线。
可选的，所述根据所述目的层的速度、所述泥岩速度和所述砂岩速度确定各沉积相模式的泥岩百分比的具体公式为：
其中，Vm某一深度点的泥岩速度，Vs为某一深度点砂岩速度，Vint为目的层层速度，Pm为泥岩百分比。
本发明还提供一种少井条件下烃源岩多参数预测系统，所述系统包括：
第一获取模块，用于获取时间域的地震体数据和测井数据；
目的层的位置确定模块，用于采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置；
地震属性信息提取模块，用于在所述目的层的位置上，从时间域的地震体数据中提取地震属性信息；
烃源岩位置确定模块，用于根据所述地震属性信息进行属性融合，并结合所述测井数据确定烃源岩位置；
沉积相模式确定模块，用于对所述测井数据中的岩性数据进行统计，确定各深度段岩性数据所对应的沉积相模式；
有机碳体数据所在位置确定模块，用于在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置；
有机碳体数据确定模块，用于按照所述有机碳体数据所在位置进行分割，获得多个不同沉积相模式下的有机碳体数据；
有机碳的预测值确定模块，用于对多个不同沉积相模式下的有机碳体数据进行组合，获得有机碳的预测值。
可选的，所述目的层的位置确定模块，具体包括：
提取单元，用于从所述时间域的地震体数据中提取多个相邻的单道地震数据；
频率域的单道地震数据确定单元，用于对各相邻的单道地震数据进行傅里叶变换，获得多个频率域的单道地震数据；
低频信息所在位置确定单元，用于根据各所述频率域的单道地震数据确定多个低频信息所在位置；
目的层位置确定单元，用于根据各所述低频信息所在位置确定烃源岩所在的目的层的位置。
可选的，所述系统还包括：
砂岩声波时差压实模型确定模块，用于利用所统计的岩性数据中砂岩和泥岩所在深度，确定不同沉积相模式下的砂岩声波时差压实模型；
泥岩和砂岩速度确定模块，用于根据所述声波时差压实模型确定各沉积相模式下泥岩速度和砂岩速度；
第二获取模块，用于获取速度谱；
目的层的速度确定模块，用于根据所述速度谱确定目的层的速度；
泥岩百分比确定模块，用于根据所述目的层的速度、所述泥岩速度和所述砂岩速度确定各沉积相模式的泥岩百分比；
目的层的地层厚度确定模块，用于根据所述烃源岩所在目的层的顶底面深度确定各沉积相模式下的目的层的地层厚度；
烃源岩初始厚度确定模块，用于根据所述泥岩百分比和所述目的层的地层厚度确定各沉积相模式的烃源岩初始厚度；
烃源岩厚度确定模块，用于对各所述沉积相模式的烃源岩初始厚度进行组合，获得目的层的烃源岩厚度。
根据本发明提供的具体实施例，本发明公开了以下技术效果：
本发明先采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置；然后在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置；最后对多个不同沉积相模式下的有机碳体数据进行组合，获得有机碳的预测值；提高了预测有机碳的精确度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例少井条件下烃源岩多参数预测方法流程图；
图2为本发明实施例少井条件下烃源岩多参数预测系统结构图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明的目的是提供一种少井条件下烃源岩多参数预测方法及系统，以提高少井条件下烃源岩多参数预测的精确度。
为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。
图1为本发明实施例少井条件下烃源岩多参数预测方法流程图，如图1所示，本发明提供一种少井条件下烃源岩多参数预测方法，所述方法包括：
步骤100：获取时间域的地震体数据和测井数据；所述测井数据包括深度信息、声波测井曲线、深侧向电阻率测井曲线和岩性数据。
步骤200：采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置；步骤200具体包括：
步骤2001：从所述时间域的地震体数据中提取多个相邻的单道地震数据。
步骤2002：对各相邻的单道地震数据进行傅里叶变换，获得多个频率域的单道地震数据；具体公式为：
其中，x(t)为所提取的单道地震数据，t代表时间，f代表频率，w(τ-t,f)为高斯窗口，τ为控制高斯窗口在时间轴即t轴位置的参数，X(t,f)为变换后频率域的单道地震数据。
步骤2003：根据各所述频率域的单道地震数据确定多个低频信息所在位置。
步骤2004：根据各所述低频信息所在位置确定烃源岩所在的目的层的位置。
步骤300：在所述目的层的位置上，从时间域的地震体数据中提取地震属性信息。
步骤500：根据所述地震属性信息进行属性融合，并结合所述测井数据确定烃源岩位置。
步骤500：对所述测井数据中的岩性数据进行统计，确定各深度段岩性数据所对应的沉积相模式。
步骤600：在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置。
在步骤600之前还包括：
步骤6001：根据所述测井数据中的声波测井曲线和深侧向电阻率测井曲线确定有机碳参数预测曲线；具体公式为：
其中，R为深侧向电阻率测井曲线值，R基线为深侧向电阻率曲线基线值，△t为声波时差测井曲线值，△t基线为声波时差测井曲线基线值，LOM为热变指数,TOCbackground为背景TOC含量，W(TOC)为有机碳参数预测曲线。
步骤6002：从有机碳参数预测曲线中提取有机碳预测值。
步骤6003：判断有机碳预测值与有机碳预测值对应的有机碳实测值之差是否小于第一设定值，如果小于或等于第一设定值，则选取有机碳参数预测曲线，如果大于第一设定值，则重新确定有机碳参数预测曲线。
步骤700：按照所述有机碳体数据所在位置进行分割，获得多个不同沉积相模式下的有机碳体数据。
步骤800：对多个不同沉积相模式下的有机碳体数据进行组合，获得有机碳的预测值。
步骤900：利用所统计的岩性数据中砂岩和泥岩所在深度，确定不同沉积相模式下的砂岩声波时差压实模型。
步骤1000：根据所述声波时差压实模型确定各沉积相模式下泥岩速度和砂岩速度。
步骤1100：获取速度谱。
步骤1200：根据所述速度谱确定目的层的速度。
步骤1300：根据所述目的层的速度、所述泥岩速度和所述砂岩速度确定各沉积相模式的泥岩百分比；具体公式为：
其中，Vm某一深度点的泥岩速度，Vs为某一深度点砂岩速度，Vint为目的层层速度，Pm为泥岩百分比。
步骤1400：根据所述烃源岩所在目的层的顶底面深度确定各沉积相模式下的目的层的地层厚度。
步骤1500：根据所述泥岩百分比和所述目的层的地层厚度确定各沉积相模式的烃源岩初始厚度。
步骤1600：对各所述沉积相模式的烃源岩初始厚度进行组合，获得目的层的烃源岩厚度。
图2为本发明实施例少井条件下烃源岩多参数预测系统结构图，如图2所示，本发明还提供一种少井条件下烃源岩多参数预测系统，所述系统包括：
第一获取模块1，用于获取时间域的地震体数据和测井数据。
目的层的位置确定模块2，用于采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置；所述目的层的位置确定模块，具体包括：
提取单元，用于从所述时间域的地震体数据中提取多个相邻的单道地震数据。
频率域的单道地震数据确定单元，用于对各相邻的单道地震数据进行傅里叶变换，获得多个频率域的单道地震数据；具体公式为：
其中，x(t)为所提取的单道地震数据，t代表时间，f代表频率，w(τ-t,f)为高斯窗口，τ为控制高斯窗口在时间轴即t轴位置的参数，X(t,f)为变换后频率域的单道地震数据。
低频信息所在位置确定单元，用于根据各所述频率域的单道地震数据确定多个低频信息所在位置。
目的层位置确定单元，用于根据各所述低频信息所在位置确定烃源岩所在的目的层的位置。
地震属性信息提取模块3，用于在所述目的层的位置上，从时间域的地震体数据中提取地震属性信息。
烃源岩位置确定模块4，用于根据所述地震属性信息进行属性融合，并结合所述测井数据确定烃源岩位置。
沉积相模式确定模块5，用于对所述测井数据中的岩性数据进行统计，确定各深度段岩性数据所对应的沉积相模式。
有机碳体数据所在位置确定模块6，用于在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置。
有机碳体数据所在位置确定模块还包括：
有机碳参数预测曲线确定单元，用于根据所述测井数据中的声波测井曲线和深侧向电阻率测井曲线确定有机碳参数预测曲线；具体公式为：
其中，R为深侧向电阻率测井曲线值，R基线为深侧向电阻率曲线基线值，△t为声波时差测井曲线值，△t基线为声波时差测井曲线基线值，LOM为热变指数,TOCbackground为背景TOC含量,W(TOC)为有机碳参数预测曲线。
有机碳预测值提取单元，用于从有机碳参数预测曲线中提取有机碳预测值。
判断单元，用于判断有机碳预测值与有机碳预测值对应的有机碳实测值之差是否小于第一设定值，如果小于或等于第一设定值，则选取有机碳参数预测曲线，如果大于第一设定值，则重新确定有机碳参数预测曲线。
有机碳体数据所在位置确定单元，用于在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置。
有机碳体数据确定模块7，用于按照所述有机碳体数据所在位置进行分割，获得多个不同沉积相模式下的有机碳体数据。
有机碳的预测值确定模块8，用于对多个不同沉积相模式下的有机碳体数据进行组合，获得有机碳的预测值。
砂岩声波时差压实模型确定模块9，用于利用所统计的岩性数据中砂岩和泥岩所在深度，确定不同沉积相模式下的砂岩声波时差压实模型。
泥岩和砂岩速度确定模块10，用于根据所述声波时差压实模型确定各沉积相模式下泥岩速度和砂岩速度。
第二获取模块11，用于获取速度谱。
目的层的速度确定模块12，用于根据所述速度谱确定目的层的速度。
泥岩百分比确定模块13，用于根据所述目的层的速度、所述泥岩速度和所述砂岩速度确定各沉积相模式的泥岩百分比；具体公式为：
其中，Vm某一深度点的泥岩速度，Vs为某一深度点砂岩速度，Vint为目的层层速度，Pm为泥岩百分比。
目的层的地层厚度确定模块14，用于根据所述烃源岩所在目的层的顶底面深度确定各沉积相模式下的目的层的地层厚度。
烃源岩初始厚度确定模块15，用于根据所述泥岩百分比和所述目的层的地层厚度确定各沉积相模式的烃源岩初始厚度。
烃源岩厚度确定模块16，用于对各所述沉积相模式的烃源岩初始厚度进行组合，获得目的层的烃源岩厚度。
本发明先采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置；然后在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置；最后对多个不同沉积相模式下的有机碳体数据进行组合，获得有机碳的预测值，进而提高了预测有机碳的精确度。
本发明首先利用所统计的岩性数据中砂岩和泥岩所在深度，确定不同沉积相模式下的砂岩声波时差压实模型；然后根据所述声波时差压实模型确定各沉积相模式下泥岩速度和砂岩速度；根据所述目的层的速度、所述泥岩速度和所述砂岩速度确定各沉积相模式的泥岩百分比；根据所述泥岩百分比和所述目的层的地层厚度确定各沉积相模式的烃源岩初始厚度；最后对各所述沉积相模式的烃源岩初始厚度进行组合，获得目的层的烃源岩厚度，进而提高了预测烃源岩厚度的精确度。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。

Claims (10)

1.一种少井条件下烃源岩多参数预测方法，其特征在于，所述方法包括：
获取时间域的地震体数据和测井数据；
采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置；
在所述目的层的位置上，从时间域的地震体数据中提取地震属性信息；
根据所述地震属性信息进行属性融合，并结合所述测井数据确定烃源岩位置；
对所述测井数据中的岩性数据进行统计，确定各深度段岩性数据所对应的沉积相模式；
在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置；
按照所述有机碳体数据所在位置进行分割，获得多个不同沉积相模式下的有机碳体数据；
对多个不同沉积相模式下的有机碳体数据进行组合，获得有机碳的预测值。
2.根据权利要求1所述的方法，其特征在于，所述采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置，具体包括：
从所述时间域的地震体数据中提取多个相邻的单道地震数据；
对各相邻的单道地震数据进行傅里叶变换，获得多个频率域的单道地震数据；
根据各所述频率域的单道地震数据确定多个低频信息所在位置；
根据各所述低频信息所在位置确定烃源岩所在的目的层的位置。
3.根据权利要求1所述的方法，其特征在于，所述在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置步骤之前还包括：
根据所述测井数据中的声波测井曲线和深侧向电阻率测井曲线确定有机碳参数预测曲线；
从有机碳参数预测曲线中提取有机碳预测值；
判断有机碳预测值与有机碳预测值对应的有机碳实测值之差是否小于第一设定值，如果小于或等于第一设定值，则选取有机碳参数预测曲线，如果大于第一设定值，则重新确定有机碳参数预测曲线。
4.根据权利要求1所述的方法，其特征在于，所述方法还包括：
利用所统计的岩性数据中砂岩和泥岩所在深度，确定不同沉积相模式下的砂岩声波时差压实模型；
根据所述声波时差压实模型确定各沉积相模式下泥岩速度和砂岩速度；
获取速度谱；
根据所述速度谱确定目的层的速度；
根据所述目的层的速度、所述泥岩速度和所述砂岩速度确定各沉积相模式的泥岩百分比；
根据所述烃源岩所在目的层的顶底面深度确定各沉积相模式下的目的层的地层厚度；
根据所述泥岩百分比和所述目的层的地层厚度确定各沉积相模式的烃源岩初始厚度；
对各所述沉积相模式的烃源岩初始厚度进行组合，获得目的层的烃源岩厚度。
5.根据权利要求2所述的方法，其特征在于，对各相邻的单道地震数据进行傅里叶变换，获得多个频率域的单道地震数据的具体公式为：
<mrow> <mi>X</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>,</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>=</mo> <munderover> <mo>&amp;Integral;</mo> <mrow> <mo>-</mo> <mi>&amp;infin;</mi> </mrow> <mrow> <mo>+</mo> <mi>&amp;infin;</mi> </mrow> </munderover> <mi>x</mi> <mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <mi>w</mi> <mrow> <mo>(</mo> <mi>&amp;tau;</mi> <mo>-</mo> <mi>t</mi> <mo>,</mo> <mi>f</mi> <mo>)</mo> </mrow> <mo>&amp;times;</mo> <msup> <mi>e</mi> <mrow> <mo>(</mo> <mo>-</mo> <mi>j</mi> <mn>2</mn> <mi>&amp;pi;</mi> <mi>f</mi> <mi>t</mi> <mo>)</mo> </mrow> </msup> <mi>d</mi> <mi>t</mi> </mrow>
其中，x(t)为所提取的单道地震数据，t代表时间，f代表频率，w(τ-t,f)为高斯窗口，τ为控制高斯窗口在时间轴即t轴位置的参数，X(t,f)为变换后频率域的单道地震数据。
6.根据权利要求3所述的方法，其特征在于，所述根据所述测井数据中的声波测井曲线和深侧向电阻率测井曲线确定有机碳参数预测曲线的具体公式为：
其中，R为深侧向电阻率测井曲线值，R基线为深侧向电阻率曲线基线值，△t为声波时差测井曲线值，△t基线为声波时差测井曲线基线值，LOM为热变指数,TOCbackground为背景TOC含量,W(TOC)为有机碳参数预测曲线。
7.根据权利要求4所述的方法，其特征在于，所述根据所述目的层的速度、所述泥岩速度和所述砂岩速度确定各沉积相模式的泥岩百分比的具体公式为：
<mrow> <msub> <mi>P</mi> <mi>m</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>V</mi> <mi>int</mi> </msub> <mo>&amp;CenterDot;</mo> <msub> <mi>V</mi> <mi>s</mi> </msub> <mo>&amp;CenterDot;</mo> <msub> <mi>V</mi> <mi>m</mi> </msub> <mo>-</mo> <msub> <mi>V</mi> <mi>m</mi> </msub> </mrow> <mrow> <msub> <mi>V</mi> <mi>s</mi> </msub> <mo>-</mo> <msub> <mi>V</mi> <mi>m</mi> </msub> </mrow> </mfrac> </mrow>
其中，Vm某一深度点的泥岩速度，Vs为某一深度点砂岩速度，Vint为目的层层速度，Pm为泥岩百分比。
8.一种少井条件下烃源岩多参数预测系统，其特征在于，所述系统包括：
第一获取模块，用于获取时间域的地震体数据和测井数据；
目的层的位置确定模块，用于采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置；
地震属性信息提取模块，用于在所述目的层的位置上，从时间域的地震体数据中提取地震属性信息；
烃源岩位置确定模块，用于根据所述地震属性信息进行属性融合，并结合所述测井数据确定烃源岩位置；
沉积相模式确定模块，用于对所述测井数据中的岩性数据进行统计，确定各深度段岩性数据所对应的沉积相模式；
有机碳体数据所在位置确定模块，用于在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置；
有机碳体数据确定模块，用于按照所述有机碳体数据所在位置进行分割，获得多个不同沉积相模式下的有机碳体数据；
有机碳的预测值确定模块，用于对多个不同沉积相模式下的有机碳体数据进行组合，获得有机碳的预测值。
9.根据权利要求8所述的系统，其特征在于，所述目的层的位置确定模块，具体包括：
提取单元，用于从所述时间域的地震体数据中提取多个相邻的单道地震数据；
频率域的单道地震数据确定单元，用于对各相邻的单道地震数据进行傅里叶变换，获得多个频率域的单道地震数据；
低频信息所在位置确定单元，用于根据各所述频率域的单道地震数据确定多个低频信息所在位置；
目的层位置确定单元，用于根据各所述低频信息所在位置确定烃源岩所在的目的层的位置。
10.根据权利要求8所述的系统，其特征在于，所述系统还包括：
砂岩声波时差压实模型确定模块，用于利用所统计的岩性数据中砂岩和泥岩所在深度，确定不同沉积相模式下的砂岩声波时差压实模型；
泥岩和砂岩速度确定模块，用于根据所述声波时差压实模型确定各沉积相模式下泥岩速度和砂岩速度；
第二获取模块，用于获取速度谱；
目的层的速度确定模块，用于根据所述速度谱确定目的层的速度；
泥岩百分比确定模块，用于根据所述目的层的速度、所述泥岩速度和所述砂岩速度确定各沉积相模式的泥岩百分比；
目的层的地层厚度确定模块，用于根据所述烃源岩所在目的层的顶底面深度确定各沉积相模式下的目的层的地层厚度；
烃源岩初始厚度确定模块，用于根据所述泥岩百分比和所述目的层的地层厚度确定各沉积相模式的烃源岩初始厚度；
烃源岩厚度确定模块，用于对各所述沉积相模式的烃源岩初始厚度进行组合，获得目的层的烃源岩厚度。
CN201711315830.3A 2017-12-12 2017-12-12 一种少井条件下烃源岩多参数预测方法及系统 Active CN108107482B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711315830.3A CN108107482B (zh) 2017-12-12 2017-12-12 一种少井条件下烃源岩多参数预测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711315830.3A CN108107482B (zh) 2017-12-12 2017-12-12 一种少井条件下烃源岩多参数预测方法及系统

Publications (2)

Publication Number Publication Date
CN108107482A true CN108107482A (zh) 2018-06-01
CN108107482B CN108107482B (zh) 2019-05-07

Family

ID=62209624

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711315830.3A Active CN108107482B (zh) 2017-12-12 2017-12-12 一种少井条件下烃源岩多参数预测方法及系统

Country Status (1)

Country Link
CN (1) CN108107482B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112394392A (zh) * 2019-08-13 2021-02-23 中国石油天然气股份有限公司 对烃源岩分布情况进行评价的方法和装置

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105319588A (zh) * 2014-07-25 2016-02-10 中国石油化工股份有限公司 一种利用地震数据进行泥页岩含气量预测方法
CN106291748A (zh) * 2015-06-26 2017-01-04 中国石油化工股份有限公司 一种获取烃源岩有机碳含量的方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105319588A (zh) * 2014-07-25 2016-02-10 中国石油化工股份有限公司 一种利用地震数据进行泥页岩含气量预测方法
CN106291748A (zh) * 2015-06-26 2017-01-04 中国石油化工股份有限公司 一种获取烃源岩有机碳含量的方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
刘军，等: "烃源岩TOC地球物理定量预测新技术及在珠江口盆地的应用", 《成都理工大学学报（自然科学版）》 *
张益明，等: "井震联合预测白云深水区海相烃源岩TOC", 《地球物理学进展》 *
李松峰，等: "稀井区烃源岩有机碳的地球物理预测方法-珠江口盆地恩平凹陷恩平组烃源岩勘探实例", 《石油地球物理勘探》 *
李林，等: "多属性预测在涸西南凹陷烃源岩研究中的应用", 《地球物理学进展》 *
殷杰，等: "利用测井和地震信息识别和预测优质烃源岩__以渤海湾盆地饶阳凹陷沙一段为例", 《天然气地球科学》 *
郭双，等: "致密油藏烃源岩有机碳含量预测方法研究", 《中国锰业》 *
陈宇航，等: "利用地震信息定量预测烃源岩TOC质量分数", 《中南大学学报(自然科学版)》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112394392A (zh) * 2019-08-13 2021-02-23 中国石油天然气股份有限公司 对烃源岩分布情况进行评价的方法和装置
CN112394392B (zh) * 2019-08-13 2023-09-26 中国石油天然气股份有限公司 对烃源岩分布情况进行评价的方法和装置

Also Published As

Publication number Publication date
CN108107482B (zh) 2019-05-07

Similar Documents

Publication Publication Date Title
CN106054248B (zh) 一种基于大面积致密储层地震岩石物理反演方法
CN104502969B (zh) 河道砂岩性油藏的识别方法
CN104698497B (zh) 一种砂泥岩薄互层层位的精细标定方法
CN103233727B (zh) 一种反演地层横波速度径向剖面的方法
CN109425896A (zh) 白云岩油气储层分布预测方法及装置
Verma et al. Estimation of total organic carbon and brittleness volume
CN104267429A (zh) 确定地层压力的方法及装置
CN105527653B (zh) 一种基于地质信息的虚拟井构建方法
CN102109611B (zh) 一种利用低频地震属性预测油藏优质储层的方法
CN103149589A (zh) 火成岩油气勘探方法及装置
CN105468886A (zh) 基于岩石物理学参数计算地层压力的方法
CN106405642B (zh) 基于脱压实声波速度的地震反演储层预测方法
CN105445800A (zh) 一种厚层砂体顶部分异岩性油藏的识别方法
CN104570067A (zh) 一种地球物理勘探中相控地震反演方法
Naseer et al. Porosity prediction of lower cretaceous unconventional resource play, south Indus Basin, Pakistan, using the seismic spectral decomposition technique
CN105116449B (zh) 一种弱反射储层的识别方法
CN104632202A (zh) 确定干粘土三孔隙度测井参数值的方法及装置
Luo et al. An ensemble-based framework for proactive geosteering
Nardone et al. Shallow velocity structure and site effects at Mt. Vesuvius, Italy, from HVSR and array measurements of ambient vibrations
CN108107482B (zh) 一种少井条件下烃源岩多参数预测方法及系统
CN110910267B (zh) 海相碳酸盐岩隔夹层的确定方法及装置
Walls et al. North Sea reservoir characterization using rock physics, seismic attributes, and neural networks; a case history
CN111983683A (zh) 一种少井条件下盐下湖相灰岩储层预测方法和系统
CN114114409A (zh) 海域天然气水合物的岩石物理建模方法、电子设备及介质
Singleton Geophysical data processing, rock property inversion, and geomechanical model building in a Midland Basin development project, Midland/Ector counties, Texas

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