一种泥页岩有机碳含量的预测方法和应用与流程

1.本发明涉一种泥页岩有机碳含量的预测方法和应用。


背景技术：

2.有机质丰度是烃源岩评价的核心参数，一般以有机碳含量表示(toc)。高有机质丰度的烃源岩是油气富集的基本物质条件。近年来，以自生自储为特征的页岩油气的勘探开发迅速发展，页岩油气的富集与烃类原位滞留与泥页岩有机质丰度特征息息相关，因此，获取泥页岩中有机质丰度参数是页岩油气勘探开发的基础工作。另一方面，中国陆相断陷湖盆发育规模相对局限，受古气候和物源供给等因素影响较大，沉积地层的岩性和岩相非均质性较强，由此有机碳含量在不同的岩性和岩相样品上表现差异较大，通过地球化学方法获取有机碳含量最为可靠，但相对成本也最高，对于非均质待测样品往往需要测定多个样品，检测成本更高，而且地球化学方法的实验周期也最长，此外还需要粉碎样品因此岩样遭到破坏，通过其他相对价廉、快速且密集的测试参数预测有机碳含量是现下比较普遍的做法。
3.cn105822301a公开了一种有机碳含量的测井预测数学模型建立方法，包括以下步骤：1)获取有机碳数据和测井曲线数据；2)作单井综合柱状图；3)对有机碳数据进行单井深度归位，调整有机碳取样点的深度位置；4)计算每条测井曲线和toc的相关系数；5)如果toc与gr相关系数最高且相关性最好，则进行下一步，否则有机碳含量的测井预测数学模型建立失败；6)作toc和gr的散点图，分别建立toc和gr的线性关系式、二次函数关系式和指数关系式；7)选取toc和gr相关性最好的关系式，该关系式就是预测目的区域多口井下一阶段toc的数学模型。该方法基于测井数据，由于测井数据采样间隔为25厘米，不能满足更低级别的诸如几厘米的纹层样品预测。但是该方法基于测井数据，由于测井数据采样间隔为25厘米，不能满足更低级别的诸如几厘米的纹层样品预测。
4.cn108510211a公开了一种烃源岩有机质丰度评价方法，包括以下步骤：结合分析化验资料，确定因素集；建立评判集，通过升半岭形分布隶属函数式，求取不同参数的隶属度，构建模糊关系矩阵；通过权重矩阵和模糊矩阵进行耦合，根据最大隶属法，评判结果。本发明的优点在于：加强了各个参数之间的联系，也避免了人为判断的影响，提高判断的准确性。但是该方法主要用于烃源岩评价打分方面，解决一般在评价生油岩有机质丰度的时候各个地化参数间出现相互矛盾的情况，并未实现烃源岩有机质丰度预测的目的。
5.cn109507733a公开了一种预测烃源岩有机质丰度的方法及装置。预测方法包括：获取目标区域的测井曲线；根据所述测井曲线计算所述目标区域的总有机碳toc值；根据所述目标区域的toc值，结合测井曲线合成地震记录，进行井震标定；根据所述井震标定的结果，提取井旁道地震属性特征，建立toc与所述地震属性特征的关系式；根据所述关系式建立目标区域toc的空间分布。通过本发明预测烃源岩有机质丰度的方法，能准确识别和评价烃源岩储层的高孔渗区分布和预测优质储层烃源岩的空间位置和物性参数分布。但是该方法基于测井数据，由于测井数据采样间隔为25厘米，不能满足更低级别的诸如几厘米的纹
层样品预测。
6.有机质丰度在自然伽马、密度、声波时差、电阻率、中子孔隙度、地震波几何形态、振幅、波形、频率、能量、相位等地球物理参数上有一定的反映。由此发展的一系列利用测井信息(如
△
log r法、地层体积密度法、自然伽马法、多元线性回归法等)及地震属性信息(如地震多属性反演法)预测有机质丰度信息。基于测井数据的有机质丰度预测模型，由于测井数据采样间隔固定为25厘米，不能满足更低级别的诸如几厘米的陆相页岩纹层样品的预测。地震属性只能应用于平面等问题。


技术实现要素：

7.针对现有技术存在的地球化学方法获取有机碳含量最为可靠，但相对成本也最高、实验周期长，此外还需要粉碎样品因此岩样遭到破坏，然而通过其他相对价廉、快速且密集的测试参数预测有机碳含量又存在不能满足更低级别的诸如几厘米的陆相页岩纹层样品的预测、只能应用于平面等问题，本发明提供一种新的泥页岩有机碳含量的预测方法，本发明的预测方法能够方便、快速和价廉地预测泥页岩的有机碳含量，无需破坏待测样品，且能满足低级别的诸如几厘米的陆相页岩纹层样品的预测。
8.本发明第一方面提供了一种泥页岩有机碳含量的预测方法，包括：
9.步骤a、通过回归方法，对泥页岩样品的有机碳含量和元素相对含量数据进行数据分析，选取数学模型；
10.步骤b、获取待预测的泥页岩的元素相对含量数据，并代入所述数学模型。
11.根据本发明所述的预测方法的一些实施方式，所述回归方法为多重逐步回归方法。
12.根据本发明所述的预测方法的一些实施方式，可以但不限于利用spss数据分析工具进行多重逐步回归方法计算。
13.根据本发明所述的预测方法的一些实施方式，步骤a的泥页岩样品为单井和/或相邻井的同一地质年代泥页岩样品。
14.根据本发明所述的预测方法的一些实施方式，步骤a的泥页岩样品可以选自泥岩和页岩，优选选自泥质白云岩、白云质泥岩、泥质灰岩、钙芒硝泥岩、灰质泥岩、砂质泥岩和页岩中的一种或多种。
15.根据本发明所述的预测方法的一些实施方式，步骤a的泥页岩样品的数量为3-40个，优选为10-35个，更优选为20-30个。
16.根据本发明所述的预测方法的一些实施方式，所述元素包括常量元素和/或微量元素，优选地，所述元素选自钛、钒、铬、锰、铁、钴、镍、铜、锌、镓、锗、锆、铌、钼、钌、铑、钯、银、铟、锡、锑、铪、钽、钨、铼、铂、金、铅、铋、镁、铝、硅、磷、氯、钾、钙、砷、铷、锶、钇、锆、镉、钡和硫元素中的一种或多种。在本发明中，测定元素后可以通过该元素的单质或氧化态形式呈现。例如但不限于mgo、al2o3、sio2、p2o5、s、cl、k2o、cao、tio2、v、cr、mno、fe2o3、co、ni、cu、zn、as、se、rb、sr、y、zr、nb、mo、ag、cd、sn、sb、ba等。
17.根据本发明所述的预测方法的一些实施方式，待预测的泥页岩为步骤a的泥页岩样品的单井和/或相邻井的同一地质年代的泥页岩。
18.根据本发明所述的预测方法的一些实施方式，所述泥页岩样品的有机碳含量获得
方法可以但不限于通过岩石热解仪和/或碳硫仪获得。
19.根据本发明所述的预测方法的一些实施方式，所述元素相对含量的获得方法以方便、快速和价格低廉为目的，例如但不限于通过手持x-射线荧光光谱仪元素扫描获得。
20.根据本发明所述的预测方法的一些实施方式，在步骤a和步骤b之间，该方法还包括判定并修正步骤a选取的相关性数学模型。例如获得数学模型后，判定该数学模型的是否符合地质认知，如果符合，则继续进行步骤b，如果不符合，则应修正为符合地质认知的相关系数最高的模型。
21.本发明第二方面提供了上述的预测方法在油气勘探与开发中的应用。
22.采用本发明的预测方法，能够方便、快速和价廉地对泥页岩有机碳含量进行预测，在优选情况下，本发明预测方法通过对少量代表性样品进行有机碳含量(toc)测试和手持荧光元素扫描(xrf)，获得样品的有机碳含量和元素相对含量信息，利用spss逐步回归方法，建立有机碳含量和元素相对含量的相关性数学模型，对同层段其它样品进行xrf扫描，利用已建立的相关性数学模型计算样品的有机碳含量，进而实现了单井或相邻井的同一地质年代其它泥页岩样品的有机碳含量预测。
附图说明
23.图1为本发明实施例1提供的预测toc与实测toc对比图。
具体实施方式
24.为使本发明更加容易理解，下面将结合实施例来详细说明本发明，这些实施例仅起说明性作用，并不局限于本发明的应用范围。
25.【实施例1】
26.步骤1、选取泥页岩样品。
27.选取江汉盆地潜江凹陷某a井潜三段10米岩心，主要发育纹层/层状泥质白云岩、白云质泥岩、泥质灰岩、钙芒硝泥岩的岩性。
28.步骤2、泥页岩样品及其测试结果。
29.根据岩心观察，人工选取其中25块典型样品作为建模数据来源，同时进行岩石热解分析(本例使用仪器型号为hawktm岩石热解仪，载气为氦气，检测器类型为fid检测器和红外检测器)和手持x-射线荧光光谱仪扫描分析(本例使用仪器型号为布鲁克s1 titantm手持式x荧光光谱仪，采集模式为双段地球化学模式，采集时间为60秒)，获得并整理测试toc和元素相对含量数值(见表1)。
30.步骤3、整理步骤2数据一一对应，输入spss，建立toc和测试元素相关性模型，设置toc为因变量，测试元素作为自变量(mgo、al2o3、sio2、p2o5、s、cl、k2o、cao、tio2、v、cr、mno、fe2o3、co、ni、cu、zn、as、se、rb、sr、y、zr、nb、mo、ag、cd、sn、sb、ba，30种元素)，计算方法为多重逐步回归方法，模型汇总见表2，模型系数见表3，
31.由此可见，自动计算最优化模型为：
32.toc(％)＝2.816 0.593
×
wt
fe2o3-5.754
×
wt
tio2-0.099
×
wt
mgo
ꢀꢀꢀ
式1
33.分析认为，泥页岩中fe2o3主要与黄铁矿、菱铁矿有关，是还原环境指示元素，其丰度与有机质含量成正比，是有机质富集的标志(崔景伟等.2013)，由此与toc成正相关性；
tio2为典型陆源组分，常与泥页岩有机质含量成反比(孙省利等.2003)；mgo与大气水及热水活动有关，对有机质保存具有破坏性作用(曹剑等，2012)，因此与toc含量成反比，由此，该模型符合地质认识。
34.步骤4、利用x-射线荧光光谱仪对toc未知样品进行元素扫描及有机碳含量预测。
35.利用手持x-射线荧光光谱仪(扫描设备与参数与步骤2相同)对同层段其它10块岩心样品进行元素扫描，获取待测样品的元素含量百分比值，整理后代入式1，计算预测toc值，检测结果及计算结果见表4。
36.为检验预测值结果，同时进行了岩石热解分析获得toc实测值，结果显示如图1所示，计算的toc预测值与实测toc值相关性较好。
37.38.39.[0040][0041]
以上所述的仅是本发明的优选实例。应当指出对于本领域的普通技术人员来说，在本发明所提供的技术启示下，作为本领域的公知常识，还可以做出其它等同变型和改进，
也应视为本发明的保护范围。