CN116299696A - 一种烃源岩泥质含量、toc、孔隙度、含水饱和度同时定量预测方法 - Google Patents

Abstract

本发明公开了一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测方法，所述步骤包括如下：步骤一：根据实际烃源岩的岩石物理特征构建相应的精确烃源岩地震岩石物理模型，可以得到烃源岩弹性参数与物性参数之间的表达式；步骤二：对精确烃源岩地震岩石物理模型进行泰勒一阶近似展开，得到线性烃源岩地震岩石物理模型；步骤三：将线性烃源岩地震岩石物理模型代入Gray近似中，得到由烃源岩泥质含量、TOC、孔隙度、含水饱和度表征的纵波反射系数；步骤四：利用叠前地震AVO反演得到的烃源岩泥质含量、TOC、孔隙度、含水饱和度预测结果，从而实现从地震数据中同时反演预测烃源岩的泥质含量、TOC、孔隙度、含水饱和度。

Description

一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测 方法

技术领域

本发明涉及一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测方法技术领域，具体为一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测方法。

背景技术

粘土含量、TOC、孔隙度和含水饱和度是烃源岩的关键物性参数，粘土含量是指烃源岩中黏土矿物的体积含量，根据黏土含量可以确定岩石的岩性，粘土含量较高的岩石更有可能是烃源岩，TOC(total organic carbon)表示有机碳含量，是指岩石中所有有机质含有的碳元素的总和占岩石总重量的百分比，与有机质的含量之间有一定的比例关系，可以用于评价岩石的生烃能力，即判定岩石是否为烃源岩，孔隙度和含水饱和度是表征烃源岩孔隙空间及含流体情况的关键物性参数，目前，现有的烃源岩预测技术缺少能够同时预测烃源岩四个关键物性参数的有效方法，与本发明最相近似的实现方案是利用实际测井数据对TOC和多个弹性参数，如纵波速度、声波阻抗、杨氏模量、泊松比等开展交会分析，找到与TOC最相关的弹性参数定义为TOC敏感弹性参数，而后通过地震反演得到敏感弹性参数剖面，再通过敏感参数与TOC之间的拟合数学关系，将敏感弹性参数剖面转换为TOC剖面，在现有技术中不能利用地震数据同时反演预测烃源岩泥质含量、TOC、孔隙度、含水饱和度。

发明内容

本发明技术方案从岩石物理出发，将精确烃源岩地震岩石物理模型线性化，得到烃源岩体积模量、剪切模量、密度同泥质含量、TOC、孔隙度、含水饱和度的线性表达关系式，进而将该关系式代入Gray近似中，得到由烃源岩泥质含量、TOC、孔隙度、含水饱和度表示的纵波反射系数，依据该反射系数构建地震AVO反演的目标函数，从而实现从地震数据中同时反演预测烃源岩的泥质含量、TOC、孔隙度、含水饱和度，解决了背景技术中所提出的问题。

为实现以上目的，本发明通过以下技术方案予以实现：一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测方法，所述步骤包括如下：

步骤一：根据实际烃源岩的岩石物理特征构建相应的精确烃源岩地震岩石物理模型，可以得到烃源岩弹性参数与物性参数之间的表达式；

如图2-5所示，为精确烃源岩地震岩石物理模型和线性烃源岩地震岩石物理模型表征的烃源岩体积模量、剪切模量、密度与物性参数泥质含量、TOC、孔隙度、含水饱和度之间的对应关系，图中虚线为依据表达式二线性烃源岩地震岩石物理模型建立的关系，可以发现，二者近乎吻合，说明可以用线性烃源岩地震岩石物理模型代替精确烃源岩地震岩石物理模型，这样做的目的是为了方便将烃源岩泥质含量、TOC、孔隙度、含水饱和度代入纵波反射系数中；

步骤二：对精确烃源岩地震岩石物理模型进行泰勒一阶近似展开，得到线性烃源岩地震岩石物理模型；

步骤三：将线性烃源岩地震岩石物理模型代入Gray近似中，得到由烃源岩泥质含量、TOC、孔隙度、含水饱和度表征的纵波反射系数，Gray近似式是由体积模量、剪切模量、密度表达的纵波反射系数；

步骤四：利用叠前地震AVO反演从四个不同角度叠加的地震数据中同时反演泥质含量、TOC、孔隙度、含水饱和度，如图6所示为四个角度的地震数据，如图7所示为依据表达式五利用叠前地震AVO反演得到的烃源岩泥质含量、TOC、孔隙度、含水饱和度预测结果利用叠前地震AVO反演得到的烃源岩泥质含量、TOC、孔隙度、含水饱和度预测结果。

优选的，所述步骤一中烃源岩弹性参数与物性参数之间的表达式一为：

其中，K,μ,ρ分别为烃源岩的体积模量、剪切模量、密度；V_c,T,

S_w分别为烃源岩的泥质含量、TOC、孔隙度、含水饱和度。

优选的，所述步骤二中对精确烃源岩地震岩石物理模型进行泰勒一阶近似展开得到表达式二：

其中，V_c0,T₀,

S_w0为泰勒一阶近似展开的参考点，/>

优选的，所述步骤三中Gray近似式是由体积模量、剪切模量、密度表达的纵波反射系数，表达式三为：

其中，θ是入射角，

分别表示体积模量反射系数、剪切模量反射系数、密度反射系数，k＝(V_s/V_p)²是烃源岩横波速度与纵波速度比值的平方，常常根据工区特征，为其赋予一个经验值。

优选的，所述

代入Gray近似中，即可得到由泥质含量、TOC、孔隙度、含水饱和度表示的纵波反射系数表达式四为：

其中，ΔV_c,ΔT,

ΔS_w分别为泥质含量差分、TOC差分、孔隙度差分和含水饱和度度差分，代表反射界面上下层泥质含量、TOC、孔隙度、含水饱和度的差异。

优选的，所述表达式四，在贝叶斯理论框架下，构建了同时预测泥质含量、TOC、孔隙度、含水饱和度的叠前地震AVO反演目标函数表达式五：

其中，

S代表地震数据，G代表耦合式(4)中权重系数的地震子波矩阵，D代表差分矩阵，/>

代表地震数据的噪音方差，/>

代表泥质含量、TOC、孔隙度、含水饱和度的方差。

本发明提供了一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测方法。具备以下有益效果：

本发明技术方案从岩石物理出发，将精确烃源岩地震岩石物理模型线性化，得到烃源岩体积模量、剪切模量、密度同泥质含量、TOC、孔隙度、含水饱和度的线性表达关系式，进而将该关系式代入Gray近似中，得到由烃源岩泥质含量、TOC、孔隙度、含水饱和度表示的纵波反射系数，依据该反射系数构建地震AVO反演的目标函数，从而实现从地震数据中同时反演预测烃源岩的泥质含量、TOC、孔隙度、含水饱和度。

附图说明

图1为本发明流程示意图；

图2为本发明泥质含量变化示意图；

图3为本发明TOC变化示意图；

图4为本发明孔隙度变化示意图；

图5为本发明含水饱和度示意图；

图6为本发明四个角度地震数据，(a)代表0°-8°叠加的地震数据，(b)代表8°-16°叠加的地震数据，(c)代表16°-24°叠加的地震数据，(d)代表24°-32°叠加的地震数据；

图7为本发明烃源岩泥质含量、TOC、孔隙度、含水饱和度预测结果剖面，(a)代表泥质含量反演剖面，(b)代表TOC反演剖面，(c)代表孔隙度反演剖面，(d)代表含水饱和度反演剖面。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

请参阅图1-7，本发明实施例提供一种技术方案：一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测方法，所述步骤包括如下：

步骤一：根据实际烃源岩的岩石物理特征构建相应的精确烃源岩地震岩石物理模型，可以得到烃源岩弹性参数与物性参数之间的表达式；

步骤二：对精确烃源岩地震岩石物理模型进行泰勒一阶近似展开，得到线性烃源岩地震岩石物理模型；

步骤三：将线性烃源岩地震岩石物理模型代入Gray近似中，得到由烃源岩泥质含量、TOC、孔隙度、含水饱和度表征的纵波反射系数，Gray近似式是由体积模量、剪切模量、密度表达的纵波反射系数；

步骤四：利用叠前地震AVO反演从四个不同角度叠加的地震数据中同时反演泥质含量、TOC、孔隙度、含水饱和度，利用叠前地震AVO反演得到的烃源岩泥质含量、TOC、孔隙度、含水饱和度预测结果。

优选的，所述步骤一中烃源岩弹性参数与物性参数之间的表达式一为：

其中，K,μ,ρ分别为烃源岩的体积模量、剪切模量、密度；V_c,T,

S_w分别为烃源岩的泥质含量、TOC、孔隙度、含水饱和度。

优选的，所述步骤二中对精确烃源岩地震岩石物理模型进行泰勒一阶近似展开得到表达式二：

其中，V_c0,T₀,

S_w0为泰勒一阶近似展开的参考点，/>

优选的，所述步骤三中Gray近似式是由体积模量、剪切模量、密度表达的纵波反射系数，表达式三为：

其中，θ是入射角，

分别表示体积模量反射系数、剪切模量反射系数、密度反射系数，k＝(V_s/V_p)²是烃源岩横波速度与纵波速度比值的平方，常常根据工区特征，为其赋予一个经验值。

优选的，所述

代入Gray近似中，即可得到由泥质含量、TOC、孔隙度、含水饱和度表示的纵波反射系数表达式四为：

其中，ΔV_c,ΔT,

ΔS_w分别为泥质含量差分、TOC差分、孔隙度差分和含水饱和度度差分，代表反射界面上下层泥质含量、TOC、孔隙度、含水饱和度的差异。

优选的，所述表达式四，在贝叶斯理论框架下，构建了同时预测泥质含量、TOC、孔隙度、含水饱和度的叠前地震AVO反演目标函数表达式五：

其中，

S代表地震数据，G代表耦合式(4)中权重系数的地震子波矩阵，D代表差分矩阵，/>

代表地震数据的噪音方差，/>

代表泥质含量、TOC、孔隙度、含水饱和度的方差。

通过将精确烃源岩地震岩石物理模型线性化，得到烃源岩体积模量、剪切模量、密度同泥质含量、TOC、孔隙度、含水饱和度的线性表达关系式，进而将该关系式代入Gray近似中，得到由烃源岩泥质含量、TOC、孔隙度、含水饱和度表示的纵波反射系数，依据该反射系数构建地震AVO反演的目标函数，从而实现从地震数据中同时反演预测烃源岩的泥质含量、TOC、孔隙度、含水饱和度。

以上显示和描述了本发明的基本原理和主要特征和本发明的优点,对于本领域技术人员而言，显然本发明不限于上述示范性实施例的细节，而且在不背离本发明的精神或基本特征的情况下，能够以其他的具体形式实现本发明。因此，无论从哪一点来看，均应将实施例看作是示范性的，而且是非限制性的，本发明的范围由所附权利要求而不是上述说明限定，因此旨在将落在权利要求的等同要件的含义和范围内的所有变化囊括在本发明内。不应将权利要求中的任何附图标记视为限制所涉及的权利要求。

此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。

Claims (6)

1.一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测方法，其特征在于：所述步骤包括如下：

步骤一：根据实际烃源岩的岩石物理特征构建相应的精确烃源岩地震岩石物理模型，可以得到烃源岩弹性参数与物性参数之间的表达式；

步骤二：对精确烃源岩地震岩石物理模型进行泰勒一阶近似展开，得到线性烃源岩地震岩石物理模型；

步骤三：将线性烃源岩地震岩石物理模型代入Gray近似中，得到由烃源岩泥质含量、TOC、孔隙度、含水饱和度表征的纵波反射系数，Gray近似式是由体积模量、剪切模量、密度表达的纵波反射系数；

步骤四：利用叠前地震AVO反演从四个不同角度叠加的地震数据中同时反演泥质含量、TOC、孔隙度、含水饱和度，利用叠前地震AVO反演得到的烃源岩泥质含量、TOC、孔隙度、含水饱和度预测结果。

2.根据权利要求1所述一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测方法，其特征在于：所述步骤一中烃源岩弹性参数与物性参数之间的表达式一为：

其中，K,μ,ρ分别为烃源岩的体积模量、剪切模量、密度；V_c,T,

S_w分别为烃源岩的泥质含量、TOC、孔隙度、含水饱和度。

3.根据权利要求1所述一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测方法，其特征在于：所述步骤二中对精确烃源岩地震岩石物理模型进行泰勒一阶近似展开得到表达式二：

其中，V_c0,T₀,

S_w0为泰勒一阶近似展开的参考点，/>

4.根据权利要求1所述一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测方法，其特征在于：所述步骤三中Gray近似式是由体积模量、剪切模量、密度表达的纵波反射系数，表达式三为：

其中，θ是入射角，

分别表示体积模量反射系数、剪切模量反射系数、密度反射系数，k＝(V_s/V_p)²是烃源岩横波速度与纵波速度比值的平方，常常根据工区特征，为其赋予一个经验值。

5.根据权利要求4所述一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测方法，其特征在于：所述

代入Gray近似中，即可得到由泥质含量、TOC、孔隙度、含水饱和度表示的纵波反射系数表达式四为：

其中，ΔV_c,ΔT,

ΔS_w分别为泥质含量差分、TOC差分、孔隙度差分和含水饱和度度差分，代表反射界面上下层泥质含量、TOC、孔隙度、含水饱和度的差异。

6.根据权利要求5所述一种烃源岩泥质含量、TOC、孔隙度、含水饱和度同时定量预测方法，其特征在于：所述表达式四，在贝叶斯理论框架下，构建了同时预测泥质含量、TOC、孔隙度、含水饱和度的叠前地震AVO反演目标函数表达式五：

其中，

S代表地震数据，G代表耦合式(4)中权重系数的地震子波矩阵，D代表差分矩阵，/>

代表地震数据的噪音方差，/>

代表泥质含量、TOC、孔隙度、含水饱和度的方差。

Images