一种少井地区的有效烃源岩丰度的预测方法与流程

本发明涉及地质层预测技术领域，尤其涉及一种少井地区的有效烃源岩丰度的预测方法。

背景技术：

有效烃源岩是指既有油气的生成又能够排出油气的岩层。近年来，一些油田的勘探失利，使人们意识到不能把研究重点局限在烃源岩的分布范围上，应该进一步深入了解烃源岩中有效部分即有机质丰度高、生烃强度大、排烃效率高、排烃厚度大的“有效的”烃源岩层的分布范围。尤其是在测井较少勘探程度较低的地区，有效烃源岩丰度与分布预测的难度更大。而且长期以来，对于有效烃源岩丰度下限值的判定，业内也存在较大争议，下限值划定太低会导致勘探资金的浪费，下限值划定太高，必然导致一些有利区被否定。

在多井地区，利用地化分析资料结合层序地层学、地震相、沉积相分析，可以对有效烃源岩丰度与分布进行预测。在少井地区，测井资料稀少，利用地化分析资料结合层序地层学、地震相、沉积相分析，对有效烃源岩丰度与分布进行预测只能圈定烃源岩发育的主要范围，无法计算得到精确烃源岩分布的厚度和toc的等值线图。预测的精度不够，只能进行定性-半定量预测，无法做到较为精确的定量预测。

技术实现要素：

有鉴于此，本发明的目的在于提供一种少井地区的有效烃源岩丰度的预测方法。本发明结合测井、地震和地化分析数据准确得到了少井地区有效烃源岩的丰度(有效烃源岩厚度、toc含量和ro)平面分布图，解决了少井地区有效烃源岩丰度预测结果不准，且预测不全面的问题。

为了实现上述发明目的，本发明提供以下技术方案：

本发明提供了一种少井地区的有效烃源岩丰度的预测方法，所述有效烃源岩丰度包括有效烃源岩的厚度、toc含量和镜质体反射率，包括以下步骤：

(1)采用油源对比分析烃源岩样本与原油样本，确定研究区内主力烃源岩发育地层；

(2)有效烃源岩厚度及分布的获取：

利用录井岩屑资料，得到所述主力烃源岩发育地层的单井泥岩百分比；

采用时深转化将测井数据从深度域转化为时间域，得到单井地震属性，以单井泥岩百分比为因变量，以单井地震属性为自变量，建立地震属性与泥岩百分比之间的线性预测模型；

将所述线性预测模型应用于主力烃源岩发育地层，得到主力烃源岩发育地层泥岩百分比的预测结果；

利用地震解释技术，得到主力烃源岩发育地层的总厚度，所述主力烃源岩发育地层的总厚度与所述主力烃源岩发育地层泥岩百分比预测结果相乘得到研究区有效烃源岩的厚度及分布；

(3)有效烃源岩toc含量分布的获取：

采用岩石热解地球化学分析确定主力烃源岩发育地层toc下限值；将单井toc含量与测井曲线值进行相关分析，选取有效测井曲线，建立单井toc含量与有效测井曲线之间的多元回归预测模型；

将所述多元预测模型应用于主力烃源岩发育地层，得到主力烃源岩发育地层toc分布曲线；

利用测井曲线计算得到地球物理参数曲线，对所述主力烃源岩发育地层的toc分布曲线与地球物理参数曲线进行相关分析，基于最小二乘原理，建立一元拟合方程；

利用地球物理反演技术，反演得到各地球物理参数数据体，并依据所述一元拟合方程将所述地球物理参数数据体转化成为主力烃源岩发育地层toc数据体；剔除所述主力烃源岩发育地层toc数据体中toc值<主力烃源岩发育地层toc下限值的数据，得到研究区有效烃源岩toc含量分布图；

(4)有效烃源岩ro分布的获取：

根据研究区镜质组反射率ro的平面等值线图，圈定2>ro>0.5的区域，得到研究区有效烃源岩ro分布图。

优选地，所述步骤(1)中油源对比分析的方法包括以下步骤：利用色谱质谱联用仪分别对所述烃源岩样本和原油样本进行检测，得到烃源岩样本与原油样本中生物标志物的m/z217和m/z191质量色谱图，比较烃源岩样本与原油样本质量色谱图的相似程度，找到与原油样本质量色谱图相似的烃源岩样本，即主力烃源岩发育地层。

优选地，所述生物标志物包括异构烷烃、甾族和多环萜类。

优选地，所述主力烃源岩发育地层的单井泥岩百分比的获取方法，包括以下步骤：利用录井岩屑资料得到录井岩性柱状图，在所述录井岩性柱状图上读取主力烃源岩发育地层的顶底界面深度，计算主力烃源岩发育地层顶底界面深度差，即为主力烃源岩发育地层的总厚度；

在所述录井岩性柱状图上读取主力烃源岩发育地层泥岩层的顶底界面深度，计算主力烃源岩发育地层泥岩层的顶底界面深度差，即为各个泥岩层厚度；将所有泥岩层厚度求和，即为主力烃源岩发育地层泥岩层总厚度；

所述主力烃源岩发育地层的单井泥岩百分比为所述主力烃源岩发育地层泥岩层总厚度与所述主力烃源岩发育地层总厚度的比值。

优选地，所述主力烃源岩发育地层toc下限值的确定方法包括：采集烃源岩井壁取芯或岩屑样本，利用lecocs-344红外碳硫测定仪测量得到所述烃源岩井壁取芯或岩屑样本的toc含量记为w(toc)；

利用岩石热解分析法，获得所述烃源岩井壁取芯或岩屑样本的热解s1值，s1值为300℃时检测的单位质量烃源岩井壁取芯或岩屑样本的烃含量；分析s1/w(toc)-w(toc)散点交会图，其中散点数据趋势线峰值处对应的toc值，即为主力烃源岩发育地层toc下限值。

优选地，所述步骤(3)中的地球物理参数包括：纵波阻抗、密度、横纵波速度比、泊松比、纵波速度和横波速度。

本发明提供了一种少井地区的有效烃源岩丰度的预测方法，本发明利用油源对比方法确定主力烃源岩发育地层，利用测井曲线计算主力烃源岩发育地层toc分布曲线。进而分析地球物理参数与toc参数之间的相关关系得到地球物理参数-toc参数的线性拟合模型，根据线性拟合模型结合地球物理反演技术得到toc分布数据体，从而得到有效烃源岩toc含量的平面分布，该方法能够定量预测有效烃源岩toc含量，能够直观得到有效烃源岩厚度与toc含量的平面分布。根据单井泥岩百分含量与地震属性值建立泥岩百分含量平面分布的预测模型，利用地震解释技术得到的烃源岩发育地层厚度，与泥岩百分含量平面分布相乘得到烃源岩厚度分布。本发明的预测有助于圈定盆地内生烃中心，具有较高的准确度。本发明解决了少井地区有效烃源岩丰度预测不全面及不准确的缺陷。

同时，本发明还提供了一种准确的有效烃源岩丰度下限值的判定方法，降低成本的同时也避免了资源的浪费。

附图说明

图1为油源对比分析结果；

图2为主力烃源岩发育地层泥岩百分比的预测结果；

图3为主力烃源岩发育地层有效烃源岩总厚度分布图；

图4为主力烃源岩发育地层toc分布曲线；

图5为有效烃源岩层toc含量分布等值线图；

图6为有效烃源岩层的ro分布图；

图7为有效烃源岩层的丰度预测图。

具体实施方式

本发明提供了一种少井地区的有效烃源岩丰度和厚度的预测方法，所述有效烃源岩丰度包括有效烃源岩的厚度、toc含量和镜质体反射率，包括以下步骤：

(1)采用油源对比分析烃源岩样本与原油样本，确定研究区内主力烃源岩发育地层；

(2)有效烃源岩厚度及分布的获取：

利用录井岩屑资料，得到所述主力烃源岩发育地层的单井泥岩百分比；

采用时深转化将测井数据从深度域转化为时间域，得到单井地震属性，以单井泥岩百分比为因变量，以单井地震属性为自变量，建立地震属性与泥岩百分比之间的线性预测模型；

将所述线性预测模型应用于主力烃源岩发育地层，得到主力烃源岩发育地层泥岩百分比的预测结果；

利用地震解释技术，得到主力烃源岩发育地层的总厚度，所述主力烃源岩发育地层的总厚度与所述主力烃源岩发育地层泥岩百分比预测结果相乘得到研究区有效烃源岩的厚度及分布；

(3)有效烃源岩toc含量分布的获取：

采用岩石热解地球化学分析确定主力烃源岩发育地层toc下限值；将单井toc含量与测井曲线值进行相关分析，选取有效测井曲线，建立单井toc含量与有效测井曲线之间的多元回归预测模型；

将所述多元预测模型应用于主力烃源岩发育地层，得到主力烃源岩发育地层toc分布曲线；

利用测井曲线计算得到地球物理参数曲线，对所述主力烃源岩发育地层的toc分布曲线与地球物理参数曲线进行相关分析，基于最小二乘原理，建立一元拟合方程；

利用地球物理反演技术，反演得到各地球物理参数数据体，并依据所述一元拟合方程将所述地球物理参数数据体转化成为主力烃源岩发育地层toc数据体；剔除所述主力烃源岩发育地层toc数据体中toc值<主力烃源岩发育地层toc下限值的数据，得到研究区有效烃源岩toc含量分布图；

(4)有效烃源岩ro分布的获取：

根据研究区镜质组反射率ro的平面等值线图，圈定2>ro>0.5的区域，得到研究区有效烃源岩ro分布图。

本发明提供的预测方法首先采用油源对比分析烃源岩样本与原油样本，确定研究区内主力烃源岩发育地层。

在本发明中，所述步骤(1)中油源对比分析的方法优选包括以下步骤：利用色谱质谱联用仪分别对所述烃源岩样本和原油样本进行检测，得到烃源岩样本与原油样本中生物标志物的m/z217和m/z191质量色谱图，比较烃源岩样本与原油样本质量色谱图的相似程度，找到与原油样本质量色谱图相似的烃源岩样本，即主力烃源岩发育地层。在本发明中，所述生物标志物优选包括异构烷烃、甾族和多环萜类。在本发明中，所述色谱的分析参数优选包括：进样方式为样品直接或用正己烷溶解后分流或不分流注入；汽化室温度为300～350℃；柱始温为70～100℃；柱终温为300～350℃；程序升温速率不超过4℃/min；载气为氦气。在本发明中，所述质谱的分析参数优选包括：离子化方式为电子轰击；电子能量为70ev；离子源温度为150～250)℃；分辨率为大于500或全质量范围为1个质量单位；扫描方式为全扫描或多离子检测；扫描质量范围为50～600uo；扫描速率为全扫描每秒不低于300u；多离子检测每个循环的时间不超过2s。在本发明中，所述色谱质谱联用仪中色谱和质谱接口温度优选为250～300℃。

在本发明中，所述相似程度的比较方法优选包括当所述烃源岩样本与原油样本特征峰的峰值数目大致相似，同一个特征峰对应的时间相差不大，整个质量色谱图的变化趋势相当，则确定该样本为主力烃源发育地层。

确定研究区内的主力烃源岩发育地层后，本发明预测有效烃源岩厚度及分布。在本发明中，所述有效烃源岩厚度的获取优选采用主力烃源岩发育地层的总厚度与主力烃源岩发育地层泥岩百分比预测结果相乘得到研究区有效烃源岩的厚度及分布。

在本发明中，所述主力烃源岩发育地层的总厚度优选通过以下步骤得到：利用petrel软件中的manualinterpretation功能，人工在地震剖面上识别画出烃源岩发育地层顶底界面；利用petrel软件中的calculate功能，输入公式thickness＝surface1(烃源岩顶界面)-surface2(烃源岩底界面)，计算得到烃源岩发育地层总厚度。

在本发明中，所述主力烃源岩发育地层泥岩百分比预测结果优选采用包括以下步骤的方法得到：利用录井岩屑资料，得到所述主力烃源岩发育地层的单井泥岩百分比；采用时深转化将测井数据从深度域转化为时间域，得到单井地震属性；以单井泥岩百分比为因变量，以单井地震属性为自变量，建立地震属性与泥岩百分比之间的线性预测模型；将所述线性预测模型应用于主力烃源岩发育地层，得到主力烃源岩发育地层泥岩百分比的预测结果。

在本发明中，所述主力烃源岩发育地层单井泥岩百分比的获取方法，优选包括以下步骤：利用录井岩屑资料得到录井岩性柱状图，在所述录井岩性柱状图上读取主力烃源岩发育地层的顶底界面深度，计算主力烃源岩发育地层顶底界面深度差，即为主力烃源岩发育地层的总厚度；在所述录井岩性柱状图上读取主力烃源岩发育地层泥岩层的顶底界面深度，计算主力烃源岩发育地层泥岩层的顶底界面深度差，即为各个泥岩层厚度；将所有泥岩层厚度求和，即为主力烃源岩发育地层泥岩层总厚度；所述主力烃源岩发育地层的单井泥岩百分比为所述主力烃源岩发育地层泥岩层总厚度与所述主力烃源岩发育地层总厚度的比值。本发明对所述录井岩屑资料的获取没有特殊的限定，采用本领域技术人员熟知的获取方式得到即可。在本发明中，所述录井岩性柱状图优选根据录井岩屑资料，利用resform软件，根据录井岩屑资料中的岩性与深度的对应关系，得到录井岩性柱状图。

在本发明中，所述单井地震属性优选采用包括以下步骤的方法得到：采用petrel软件中的geologicalinterpretation系统中的wellcorrelation模块进行时深转换，利用petrel软件中的seismicinterpretation系统中的multitraceattributes模块进行地震属性的提取，确定单井的地震属性。在本发明中，所述multitraceattributes模块的参数设置优选包括将windowlength设置为11；将numberofiteration设置为50。

在本发明中，所述线性预测模型为式i所示：

式i中，pm为泥岩百分含量，li为地震属性值，i为地震属性li与泥岩百分含量pm之间的相关系数，p为常数，n为测井个数。在本发明中，所述i和p可以利用spss(统计产品与服务解决方案)软件，对多个测井计算得到地震属性值和泥岩百分含量的值进行多元回归分析得到。

确定主力烃源岩发育地层后，获取主力烃源岩发育地层有效烃源岩toc含量分布。在本发明中，所述效烃源岩toc含量分布的获取优选包括以下步骤：采用岩石热解地球化学分析确定主力烃源岩发育地层toc下限值；将主力烃源岩发育地层toc值与测井曲线值进行相关分析，选取有效测井曲线，建立主力烃源岩发育地层toc含量与有效测井曲线之间的多元回归预测模型；将所述多元预测模型应用于主力烃源岩发育地层，得到主力烃源岩发育地层toc分布曲线；利用测井曲线计算得到地球物理参数曲线，对所述主力烃源岩发育地层的toc分布曲线与地球物理参数曲线进行相关分析，基于最小二乘原理，建立一元拟合方程；利用地球物理反演技术，反演得到各地球物理参数数据体，并依据所述一元拟合方程将所述地球物理参数数据体转化成为主力烃源岩发育地层toc数据体；剔除所述主力烃源岩发育地层toc数据体中toc值<主力烃源岩发育地层toc下限值的数据，得到研究区有效烃源岩toc含量分布图。

在本发明中，所述主力烃源岩发育地层toc下限值的确定方法包括：采集烃源岩井壁取芯或岩屑样本，利用lecocs-344红外碳硫测定仪测量得到所述烃源岩井壁取芯或岩屑样本的toc含量记为w(toc)，利用岩石热解分析法，获得所述烃源岩井壁取芯或岩屑样本的热解s1值，s1值为300℃时检测的单位质量烃源岩井壁取芯或岩屑样本中的烃含量；分析s1/w(toc)-w(toc)散点交会图，其中散点数据趋势线峰值处对应的toc值，即为主力烃源岩发育地层toc下限值。

本发明将单井toc含量与测井曲线值进行相关分析，选取有效测井曲线，建立所述单井toc含量与所述有效测井曲线之间的多元回归预测模型。在本发明中，所述单井toc含量的获取方法优选包括利用lecocs-344红外碳硫测定仪测量单井toc含量。在本发明中，对所述相关分析的方法没有特殊的限定，采用本领域技术人员熟知的分析方法即可，具体的，如以测井曲线作为a变量，以单井toc含量作为b变量，利用spss软件中相关分析模块，计算得到a、b变量之间的皮尔斯相关系数。在本发明中，所述ssps软件中，优选将相关参数选择为pearson，显著性检验方式选择为双尾检验(t)。

在本发明中，所述测井曲线优选包括自然伽马曲线、中子孔隙度曲线、密度曲线、声波时差曲线、电阻率曲线，自然电位曲线。

在本发明中，所述多元回归预测模型优选利用spss(statisticalproductandservicesolutions)软件进行多元线性回归分析获得。

在本发明中，所述多元回归预测模型优选为式ii所示：

式ii中w(toc)为烃源岩样本toc含量，bi、d为常数，mi是为有效测井曲线，个数为m；bi、d利用spss(statisticalproductandservicesolutions)软件中多元线性回归分析模块求取。

得到多元回归预测模型后，本发明将所述多元回归预测模型应用于主力烃源岩发育地层，得到主力烃源岩发育地层toc分布曲线。

在本发明中，当测井曲线缺少自然伽马曲线、中子孔隙度曲线、密度曲线和自然电位曲线时，并且利用多元回归模型建立预测模型计算得到的toc分布曲线与实际测量的烃源岩toc含量相关性不佳时，优选采用超声波曲线和地层真电阻率曲线重叠法计算得到主力烃源岩发育地层toc分布曲线。在本发明中，所述相关性的计算方法优选为利用spss软件中的相关分析模块，把测量得到的toc含量看成a变量，把计算得到的toc含量看做是b变量。计算测量得到的toc含量(a变量)与计算得到的toc含量(b变量)之间的皮尔斯相关系数，相关分析软件为spss(statisticalproductandservicesolutions)软件，相关分析模块参数选择：相关参数类型选择为pearson，即皮尔森相关系数。显著性检验方式选择为：双尾检验(t)。在本发明中，当所述相关系数小于0.5时，可以认定为相关性不佳。在本发明中，所述超声波曲线和地层真电阻率曲线重叠法为本领域技术人员熟知的计算方法。

在本发明中，当地层真电阻率曲线缺失的情况下，优选采用深侧向电阻率曲线代替，主要的计算公式如式iii所示：

式iii中，r为地层真电阻率，r基线为地层真电阻率基线值，δt为声波时差曲线，δt基线为声波时差曲线基线值，lom热变指数，tocbackground为toc含量背景值。

得到主力烃源岩发育地层toc分布曲线后，本发明利用测井曲线计算得到地球物理参数曲线，对所述主力烃源岩发育地层的toc分布曲线与地球物理参数曲线进行相关分析，基于最小二乘原理，建立一元拟合方程。

在本发明中，优选采用geoview软件中的elog模块中的crossplot功能进行地球物理曲线的计算。在本发明中，所述geoview软件中elog模块中crossplot功能的具体参数设置优选包括：打开crossplot功能，选取主力烃源岩发育地层的toc分布曲线为a变量，地球物理参数(密度、纵波阻抗、泊松比、横纵波速度比等)为b变量，得到地球物理参数。在本发明中，优选采用regression功能中的leastsquares功能进行所述地球物理曲线与主力烃源岩发育地层的toc分布曲线之间的相关性分析。在本发明中，所述regression功能中的leastsquares功能的具体参数设置优选包括numberofiteration为50；typeofnorm(1.0-2.0)为2.0；normalizedstandarderror为0。

在本发明中，所述一元拟合方程优选式iv所示：

式iv中，d为常数，d为地球物理参数。在本发明中，所述常数d、bi的获取方法优选为利用spss(statisticalproductandservicesolutions)通过一元回归分析得到。

本发明利用地球物理反演技术，反演得到各种地球物理参数数据体，并依据所述一元拟合方程将所述地球物理参数数据体转化成为toc数据体；然后剔除所述toc数据体中toc<主力烃源岩发育地层toc下限值的数据，得到有效烃源岩的toc含量。

在本发明中，所述地球反演技术针对的数据优选为利用叠后地震数据进行纵波阻抗反演，得到叠后纵波阻抗数据体。在本发明中，所述地球物理参数数据体转化为toc数据体的方法优选为在jason软件中利用functionmod模块将toc数据体通过地球物理参数与toc参数之间的一元拟合方程转成成为toc参数数据体。

确定主力烃源岩发育地层后，根据研究区镜质组反射率ro的平面等值线图，圈定2>ro>0.5的区域，得到研究区有效烃源岩ro分布图。

下面结合实施例对本发明提供少井地区的有效烃源岩丰度的预测方法进行详细的说明，但是不能把它们理解为对本发明保护范围的限定。

实施例1

以东海盆地丽水-椒江凹陷为例：

(1)采用油源对比分析烃源岩样本与原油样本中，确定研究区内主力烃源岩发育地层；油源对比分析结果如图1所示。从图1可以看出：原油与烃源岩样本的生物标志物分布特征很相似，表明具有良好的亲缘关系。

(2)有效烃源岩厚度及分布的获取：

利用录井岩屑资料，得到所述主力烃源岩发育地层的单井泥岩百分比分别为0.7、0.65、0.68、0.75、0.58；

采用时深转化将测井数据从深度域转化为时间域，得到单井地震属性为均方根频率为：155.3hz、167.2hz、170.9hz、149.7hz、158.6hz；均方根振幅为：1011.3、770.9、690.9、1126.1、1140.1；以单井泥岩百分比为因变量，以单井地震属性为自变量，建立地震属性与泥岩百分比之间的线性预测模型为r＝-0.00045x-0.00924y+1.89；其中，r为泥岩百分含量；x为均方根振幅；y为均方根频率；

将所述线性预测模型应用于主力烃源岩发育地层，得到主力烃源岩发育地层泥岩百分比的预测结果见图2；从图2可以看出：主力烃源岩发育地层泥岩百分比较高，表示在烃源岩发育地层中，有较厚的泥岩，即有效烃源岩厚度较大。

利用地震解释技术，得到主力烃源岩发育地层的总厚度分布，将所述主力烃源岩发育地层的总厚度与所述主力烃源岩发育地层泥岩百分比预测结果相乘得到主力烃源岩发育地层有效烃源岩的厚度分布图，见图3，从图3可以看出：丽水-椒江凹陷西次凹有效烃源岩分布厚度最大，烃源岩厚度大于400m，烃源岩条件好，预测结果符合地质规律。

(3)有效烃源岩toc含量分布的获取：

采用岩石热解地球化学分析确定主力烃源岩发育地层toc下限值为0.5％；将单井toc含量与测井曲线值进行相关分析，选取有效测井曲线(中子孔隙度曲线、自然伽马曲线、密度曲线、声波时差曲线、地层真电阻率曲线)，建立单井toc含量与有效测井曲线之间的多元回归预测模型为：w(toc)＝3.532*cn-0.014*gr-0.025*rt+2.336*den-0.037*ac-0.787，其中cn为中子孔隙度测井曲线、gr为自然伽马测井曲线、rt地层真电阻率测井曲线、ac为声波时差测井曲线、den为密度测井曲线；

将所述多元预测模型应用于主力烃源岩发育地层，得到主力烃源岩发育地层toc分布曲线，见图4。从图4可以看出：利用多元回归分析方法计算得到的单井toc分布曲线与实际测量得到的toc数据拟合较好。表明预测结果具有可靠性。

利用测井曲线计算得到地球物理参数曲线，对所述主力烃源岩发育地层的toc分布曲线与地球物理参数曲线进行相关分析，基于最小二乘原理，建立一元拟合方程：w(toc)＝5.179*10-5*p-impedance+0.483，其中p-impedance为纵波阻抗(地球物理参数)；

利用地球物理反演技术，反演得到各地球物理参数数据体，并依据所述一元拟合方程将所述地球物理参数数据体转化成为主力烃源岩发育地层toc数据体；剔除所述主力烃源岩发育地层toc数据体中toc值<主力烃源岩发育地层toc下限值的数据，得到有效烃源岩层toc含量分布图，参见图5。从图5可以看出：在丽水-椒江凹陷西次凹北部靠近灵峰古潜山附近toc含量分布为1.5～3％，具有很高的toc含量，分布较好的烃源岩。

(4)有效烃源岩ro分布的获取：

根据研究区镜质组反射率ro的平面等值线图，圈定2>ro>0.5的区域，得到有效烃源岩层ro分布图，结果见6，从图6可以看出：在丽水-椒江凹陷西次凹不仅具有较厚的烃源岩，并且烃源岩的toc含量较高，且具有较高的成熟度(ro>0.5)，所以丽水-椒江凹陷西次凹具有很好的烃源岩条件，丽水-椒江凹陷东次凹烃源岩条件较西次凹差。

将得到的有效烃源岩厚度分布、toc含量分布及ro参数拟合到一张图中，得到有效烃源层的丰度预测图，参见图7。从图7可以看出：颜色的深浅代表烃源岩厚度的大小，等值线的分布代表toc含量的分布、其中整个数据分布的范围是0.5<ro<2的范围。

从实施例可以看出，本发明的有效烃源岩丰度结果准确，解决了少井地区有效烃源岩丰度预测不全面及不准确的缺陷。

以上所述仅是本发明的优选实施方式，应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。