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

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

【技术实现步骤摘要】
一种少井条件下烃源岩多参数预测方法及系统
本专利技术涉及烃源岩预测
，特别是涉及一种少井条件下烃源岩多参数预测方法及系统。
技术介绍
在盆地勘探初期,预测烃源岩对烃源岩位置的识别、对烃源岩厚度的预测以及对烃源岩质量的评价具有至关重要的作用。在多井地区，地球化学分析方法多应用于烃源岩的预测与评价，但是在少井地区，尤其是海上勘探区，井少直接制约着烃源岩的预测。故研究学者采用井震结合实现对烃源岩进行预测，先根据烃源岩“低频-强振幅”的特征在地震剖面上肉眼确定烃源岩位置，在烃源岩位置内通过地球物理反演与测井结合的方法得到有机碳预测值，但是现有预测过程中，由于现有根据频率高低划分烃源发育位置时，只是单从地震剖面上人为观察，是不定量的科学研究，同时地球物理反演过程中没有考虑沉积相模式的不同，因此上述预测有机碳存在预测精度低的问题。另外，现有在烃源岩位置上通过计算砂地比得到烃源岩厚度时，同一个地区不同的沉积相模式下均采用同一个预测模型，因此预测出的烃源岩的厚度也存在预测精确度低的问题。基于上述问题，如何提高烃源岩参数预测的精度成为本领域亟需解决的问题。
技术实现思路
本专利技术的目的是提供一种少井条件下烃源岩多参数预测方法及系统，以提高少井条件下烃源岩多参数预测的精确度。为实现上述目的，本专利技术提供一种少井条件下烃源岩多参数预测方法，所述方法包括：获取时间域的地震体数据和测井数据；采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置；在所述目的层的位置上，从时间域的地震体数据中提取地震属性信息；根据地震属性信息进行属性融合，并结合所述测井数据确定烃源岩位置；对所述测井数据中的岩性数据进行统计，确定各深度段岩性数据所对应的沉积相模式；在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置；按照所述有机碳体数据所在位置进行分割，获得多个不同沉积相模式下的有机碳体数据；对多个不同沉积相模式下的有机碳体数据进行组合，获得有机碳的预测值。可选的，所述采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置，具体包括：从所述时间域的地震体数据中提取多个相邻的单道地震数据；对各相邻的单道地震数据进行傅里叶变换，获得多个频率域的单道地震数据；根据各所述频率域的单道地震数据确定多个低频信息所在位置；根据各所述低频信息所在位置确定烃源岩所在的目的层的位置。可选的，所述在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置步骤之前还包括：根据所述测井数据中的声波测井曲线和深侧向电阻率测井曲线确定有机碳参数预测曲线；从有机碳参数预测曲线中提取有机碳预测值；判断有机碳预测值与有机碳预测值对应的有机碳实测值之差是否小于第一设定值，如果小于或等于第一设定值，则选取有机碳参数预测曲线，如果大于第一设定值，则重新确定有机碳参数预测曲线。可选的，所述方法还包括：利用所统计的岩性数据中砂岩和泥岩所在深度，确定不同沉积相模式下的砂岩声波时差压实模型；根据所述声波时差压实模型确定各沉积相模式下泥岩速度和砂岩速度；获取速度谱；根据所述速度谱确定目的层的速度；根据所述目的层的速度、所述泥岩速度和所述砂岩速度确定各沉积相模式的泥岩百分比；根据所述烃源岩所在目的层的顶底面深度确定各沉积相模式下的目的层的地层厚度；根据所述泥岩百分比和所述目的层的地层厚度确定各沉积相模式的烃源岩初始厚度；对各所述沉积相模式的烃源岩初始厚度进行组合，获得目的层的烃源岩厚度。可选的，对各相邻的单道地震数据进行傅里叶变换，获得多个频率域的单道地震数据的具体公式为：其中，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所述的方法，其特征在于，所述采用傅里叶变换，根据时间域的地震体数据确定烃源岩所在的目的层的位置，具体包括：从所述时间域的地震体数据中提取多个相邻的单道地震数据；对各相邻的单道地震数据进行傅里叶变换，获得多个频率域的单道地震数据；根据各所述频率域的单道地震数据确定多个低频信息所在位置；根据各所述低频信息所在位置确定烃源岩所在的目的层的位置。3.根据权利要求1所述的方法，其特征在于，所述在不同沉积相模式下，利用有机碳参数预测曲线和地质统计学反演技术确定烃源岩位置上不同沉积相模式下的有机碳体数据所在位置步骤之前还包括：根据所述测井数据中的声波测井曲线和深侧向电阻率测井曲线确定有机碳参数预测曲线；从有机碳参数预测曲线中提取有机碳预测值；判断有机碳预测值与有机碳预测值对应的有机碳实测值之差是否小于第一设定值，如果小于或等于第一设定值，则选取有机碳参数预测曲线，如果大于第一设定值，则重新确定有机碳参数预测曲线。4.根据权利要求1所述的方法，其特征在于，所述方法还包括：利用所统计的岩性数据中砂岩和泥岩所在深度，确定不同沉积相模式下的砂岩声波时差压实模型；根据所述声波时差压实模型确定各沉积相模式下泥岩速度和砂岩速度；获取速度谱；根据所述速度谱确定目的层的速度；根据所述目的层的速度、所述泥岩速度和所述砂岩速度确定各沉积相模式的泥岩百分比；根据所述烃源岩所在目的层的顶底面深度确定各沉积相模式下的目的层的地层厚度；根据所述泥岩百分比和所述目的层的地层厚度确定各沉积相模式的烃源岩初始厚度；对各所述沉积相模式的烃源岩初始厚度进行组合，获得目的层的烃源岩厚度。5.根据权利要求2所述的方法，其特征在于，对各相邻的单道地震数据进行傅里叶变换，获得多个频率域的单道地震数据的具体公式为：其中，x(t)为所提取的单道地震数据，t代表时间，f代表频率，w(τ-t,f)为高斯窗口，τ为控制高斯窗口在时间轴即t轴位置的参数，X(t,f)为变换后频率域的单道地震数据。6.根据权利要求3所述的方法，其特征在于，所述根据所述测井数据中的声波测井曲线和深侧向电阻率测井...

【专利技术属性】
技术研发人员：张金亮，李娜，刘金水，申雯龙，许广臣，
申请(专利权)人：北京师范大学，
类型：发明
国别省市：北京,11