用于判别灯影组二段和灯影组四段的装置及方法与流程

本发明涉及石油勘探技术领域，特别涉及一种用于判别灯影组二段和灯影组四段的装置及方法。

背景技术：

四川盆地是重要的油气勘探地区之一，四川盆地中灯影组分布广泛。以四川盆地震旦系灯影组为例，根据生物、岩性、电性特征，将该灯影组沿纵向划分为四段：灯影组一段，厚度为60-100m，贫藻特征；灯影组二段，厚度为400-500m，富藻特征；灯影组三段，厚度为10-60m，贫藻特征；灯影组四段，厚度为200-300m，富藻特征。可见，灯影组二段和灯影组四段的生物、岩性、电性特征相似，而钻井进入灯影组的200m以内便完钻，因此导致部分钻井完钻的层位存在争议，制约了油气勘探的研究工作。基于上述，有必要提供一种用于判别灯影组二段和灯影组四段的方法。

相关技术中，通过生物、岩性、电性特征等个人经验来判别灯影组二段和灯影组四段，但是由于受沉积古地理影响，使不同部位的震旦系灯影组的岩性、厚度和生物组合特征不完全一致。或者受桐湾运动的影响，使灯影组四段遭受不通过程度的淋滤和剥蚀，造成地层厚度差异较大。因此，不能准确地判别灯影组二段和灯影组四段。

技术实现要素：

本发明实施例提供了一种用于判别灯影组二段和灯影组四段的装置及方法，可解决上述技术问题。技术方案如下：

一方面，本发明实施例提供了一种用于判别灯影组二段和灯影组四段的方法，所述方法包括：

获取标准灯影组二段和标准灯影组四段的岩心样本库；

根据所述岩心样本库，分别获取标准灯影组二段和标准灯影组四段的全岩元素分布特征；

根据所述岩心样本库，分别获取标准灯影组二段和标准灯影组四段的碳同位素分布频率特征；

根据所述岩心样本库，分别获取标准灯影组二段和标准灯影组四段的碳同位素和氧同位素的关系特征；

获取目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征；

若所述目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征均符合标准灯影组二段或标准灯影组四段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征，确定所述目标岩心样本属于灯影组二段或灯影组四段。

在一种可能的设计中，所述根据所述岩心样本库，分别获取标准灯影组二段和标准灯影组四段的全岩元素分布特征，包括：

根据所述岩心样本库，分别获取标准灯影组二段所对应的每种岩心样本的全元素及质量百分比；

根据所述标准灯影组二段所对应的每种岩心样本的全元素及质量百分比，获取标准灯影组二段所对应的每种岩心样本的全元素及质量百分比之间的第一关系示意图；

根据所述岩心样本库，分别获取标准灯影组四段所对应的每种岩心样本的全元素及质量百分比；

根据所述标准灯影组四段所对应的每种岩心样本的全元素及质量百分比，获取标准灯影组四段所对应的每种岩心样本的全元素及质量百分比之间的第二关系示意图。

在一种可能的设计中，所述根据所述岩心样本库，分别获取标准灯影组二段和标准灯影组四段的碳同位素分布频率特征，包括：

根据所述岩心样本库，获取标准灯影组二段所对应的每种岩心样本的碳同位素的值及频率；

根据所述标准灯影组二段所对应的每种岩心样本的碳同位素的值及频率，获取所述标准灯影组二段所对应的多种岩心样本的碳同位素的值及频率的第三关系示意图；

根据所述岩心样本库，获取标准灯影组四段所对应的每种岩心样本的碳同位素的值及频率；

根据所述标准灯影组四段所对应的每种岩心样本的碳同位素的值及频率，获取所述标准灯影组四段所对应的多种岩心样本的碳同位素的值及频率的第四关系示意图。

在一种可能的设计中，所述根据所述岩心样本库，分别获取标准灯影组二段和标准灯影组四段的碳同位素和氧同位素的关系特征，包括：

根据所述岩心样本库，获取标准灯影组二段所对应的每个岩心样本的碳同位素的值和氧同位素的值；

根据所述岩心样本库，获取标准灯影组四段所对应的每个岩心样本的碳同位素的值和氧同位素的值；

以氧同位素的值为横坐标，碳同位素的值为纵坐标，获取标准灯影组二段所对应的每个岩心样本和标准灯影组四段所对应的每个岩心样本的散点图。

在一种可能的设计中，所述方法还包括：

若所述目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征中至少一种不符合标准灯影组二段或标准灯影组四段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征，确定所述目标岩心样本不属于灯影组二段或灯影组四段。

另一方面，本发明实施例提供了一种用于判别灯影组二段和灯影组四段的装置，所述装置包括：

第一获取模块，用于获取标准灯影组二段和标准灯影组四段的岩心样本库；

第二获取模块，用于根据所述岩心样本库，分别获取标准灯影组二段和标准灯影组四段的全岩元素分布特征；

第三获取模块，用于根据所述岩心样本库，分别获取标准灯影组二段和标准灯影组四段的碳同位素分布频率特征；

第四获取模块，用于根据所述岩心样本库，分别获取标准灯影组二段和标准灯影组四段的碳同位素和氧同位素的关系特征；

第五获取模块，用于获取目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征；

第一确定模块，用于若所述目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征均符合标准灯影组二段或标准灯影组四段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征，确定所述目标岩心样本属于灯影组二段或灯影组四段。

在一种可能的设计中，所述第二获取模块，包括：

第一获取单元，用于根据所述岩心样本库，分别获取标准灯影组二段所对应的每种岩心样本的全元素及质量百分比；

第二获取单元，用于根据所述标准灯影组二段所对应的每种岩心样本的全元素及质量百分比，获取标准灯影组二段所对应的每种岩心样本的全元素及质量百分比之间的第一关系示意图；

第三获取单元，用于根据所述岩心样本库，分别获取标准灯影组四段所对应的每种岩心样本的全元素及质量百分比；

第四获取单元，用于根据所述标准灯影组四段所对应的每种岩心样本的全元素及质量百分比，获取标准灯影组四段所对应的每种岩心样本的全元素及质量百分比之间的第二关系示意图。

在一种可能的设计中，所述第三获取模块，包括：

第五获取单元，用于根据所述岩心样本库，获取标准灯影组二段所对应的每种岩心样本的碳同位素的值及频率；

第六获取单元，用于根据所述标准灯影组二段所对应的每种岩心样本的碳同位素的值及频率，获取所述标准灯影组二段所对应的每种岩心样本的碳同位素的值及频率的第三关系示意图；

第七获取单元，用于根据所述岩心样本库，获取标准灯影组四段所对应的每种岩心样本的碳同位素的值及频率；

第八获取单元，用于根据所述标准灯影组四段所对应的每种岩心样本的碳同位素的值及频率，获取所述标准灯影组四段所对应的每种岩心样本的碳同位素的值及频率的第四关系示意图。

在一种可能的设计中，所述第四获取模块，包括：

第九获取单元，用于根据所述岩心样本库，获取标准灯影组二段所对应的每个岩心样本的碳同位素的值和氧同位素的值；

第十获取单元，用于根据所述岩心样本库，获取标准灯影组四段所对应的每个岩心样本的碳同位素的值和氧同位素的值；

第十一获取单元，用于以氧同位素的值为横坐标，碳同位素的值为纵坐标，获取标准灯影组二段所对应的每个岩心样本和标准灯影组四段所对应的每个岩心样本的散点图。

在一种可能的设计中，所述装置还包括：

第二确定模块，用于若所述目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征中至少一种不符合标准灯影组二段或标准灯影组四段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征，确定所述目标岩心样本不属于灯影组二段或灯影组四段。

本发明实施例提供的技术方案带来的有益效果至少包括：

本发明实施例提供的用于判别灯影组二段和灯影组四段的方法，通过获取目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征，并与标准灯影组二段和标准灯影组四段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征进行比对，能够准确地、快速地、经济地确认目标岩心样本属于灯影组二段或灯影组四段，利于油气勘探工作的研究。

附图说明

为了更清楚地说明本发明实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是本发明实施例提供的用于判别灯影组二段和灯影组四段的方法流程图；

图2是本发明实施例提供的用于判别灯影组二段和灯影组四段的装置结构示意图；

图3是实施例1提供的hs1井灯影组二段岩心样本的全岩元素与质量百分比之间的直方图；

图4是实施例1提供的z6井灯影组二段岩心样本的全岩元素与质量百分比之间的直方图；

图5是实施例1提供的gs7井灯影组四段岩心样本的全岩元素与质量百分比之间的直方图；

图6是实施例1提供的mx9井灯影组四段岩心样本的全岩元素与质量百分比之间的直方图；

图7是实施例1提供的bhs1井目标岩心样本的全岩元素与质量百分比之间的直方图；

图8是实施例1提供的灯影组二段和灯影组四段分别对应的岩心样本、bhs1井目标岩心样本的碳同位素的值与碳同位素的频数之间的直方图；

图9是实施例1提供的灯影组二段和灯影组四段分别对应的岩心样本、bhs1井目标岩心样本的碳同位素和氧同位素之间的散点图。

具体实施方式

除非另有定义，本发明实施例所用的所有技术术语均具有与本领域技术人员通常理解的相同的含义。为使本发明的目的、技术方案和优点更加清楚，下面将结合附图对本发明实施方式作进一步地详细描述。

一方面，本发明实施例提供了一种用于判别灯影组二段和灯影组四段的方法，如附图1所示，该方法包括：

步骤11、获取标准灯影组二段和标准灯影组四段的岩心样本库。

可以向终端中直接输入标准灯影组二段和标准灯影组四段的岩心样本库数据，终端也可以通过采集设备来获取标准灯影组二段和标准灯影组四段的岩心样本库数据。

标准灯影组二段和标准灯影组四段的岩心样本可以直接采集无争议钻井岩心或野外露头的岩心，在取样时应尽量避免采集孔洞发育的样本，选择致密的泥晶白云岩样本。

步骤12、根据岩心样本库，分别获取标准灯影组二段和标准灯影组四段的全岩元素分布特征。

步骤12包括但不限于以下子步骤：

步骤121、根据岩心样本库，分别获取标准灯影组二段所对应的每种岩心样本的全元素及质量百分比。

需要说明的是，在岩心样本库中，标准灯影组二段对应有多种岩心样本，每种岩心样本对应有多个岩心样本。获取的每种岩心样本的全元素对应的质量百分比可以为多个岩心样本的平均数据。

通过《gb/t14506.30-2010硅酸盐岩石化学分析方法第30部分:44个元素量测定》标准提供的方法可以检测每种岩心样本的全元素及其质量百分比。

其中，每个岩心样本的全元素包括但不限于：钙、镁、硅、铁、铝、钠、钾、磷、锰、钛。

在灯影组二段和灯影组四段中，镁的质量百分比与钙的质量百分比之比存在差异，因此，也可以基于镁的质量百分比与钙的质量百分比之比这一特征来进行分析。

步骤122、根据标准灯影组二段所对应的每种岩心样本的全元素及质量百分比，获取标准灯影组二段所对应的每种岩心样本的全元素及质量百分比之间的第一关系示意图。

示例地，第一关系示意图为以元素为横坐标，以质量百分比为纵坐标的直方图，以直观、清晰地观察标准灯影组二段所对应的每种岩心样本的全元素的分布情况。

步骤123、根据岩心样本库，分别获取标准灯影组四段所对应的每种岩心样本的全元素及质量百分比。

步骤123与步骤121的方法相同，此处不再赘述。

步骤124、根据标准灯影组四段所对应的每种岩心样本的全元素及质量百分比，获取标准灯影组四段所对应的每种岩心样本的全元素及质量百分比之间的第二关系示意图。

示例地，第二关系示意图为以元素为横坐标，以质量百分比为纵坐标的直方图，以直观、清晰地观察标准灯影组四段所对应的每种岩心样本的全元素的分布情况。

步骤13、根据岩心样本库，分别获取标准灯影组二段和标准灯影组四段的碳同位素分布频率特征。

灯影组属于富藻地层，生物对环境具有良好的敏感性，同时，生物对碳同位素具有“分馏”作用，因此不同时代地层的生物兴衰程度存在差异，这使得不同段的灯影组的碳同位素分布频率不同。基于此，可以通过碳同位素分布频率特征来区别灯影组二段和灯影组四段。

步骤13包括但不限于以下子步骤：

步骤131、根据岩心样本库，获取标准灯影组二段所对应的每种岩心样本的碳同位素的值及频率。

根据《有机物和碳酸盐岩碳、氧同位素分析方法》标准提供的方法检测标准灯影组二段所对应的每个岩心样本的碳同位素的值，并统计每种岩心样本的碳同位素的值出现的频率。通过获取每种岩心样本中每个岩心样本的碳同位素的值，并取平均值，得到每种岩心样本的碳同位素的值。

步骤132、根据标准灯影组二段所对应的每种岩心样本的碳同位素的值及频率，获取标准灯影组二段所对应的多种岩心样本的碳同位素的值及频率的第三关系示意图。

示例地，第三关系示意图为以每种岩心样本所对应的碳同位素的值为横坐标，以频率为纵坐标的直方图，以直观、清晰地观察标准灯影组二段所对应的多种岩心样本的碳同位素的值与频率的关系。

步骤133、根据岩心样本库，获取标准灯影组四段所对应的每种岩心样本的碳同位素的值及频率。

同理，根据《syt5238-2008有机物和碳酸盐岩碳、氧同位素分析方法》标准提供的方法检测标准灯影组四段所对应的每个岩心样本的碳同位素的值，并统计每种岩心样本的碳同位素的值出现的频率。通过获取每种岩心样本中每个岩心样本的碳同位素的值，并取平均值，得到每种岩心样本的碳同位素的值。

步骤134、根据标准灯影组四段所对应的每种岩心样本的碳同位素的值及频率，获取标准灯影组四段所对应的多种岩心样本的碳同位素的分布频率的质量百分比及频率的第四关系示意图。

示例地，第四关系示意图为以每种岩心样本所对应的碳同位素的值为横坐标，以频率为纵坐标的直方图，以直观、清晰地观察标准灯影组四段所对应的多种岩心样本的碳同位素的值与频率的关系。

步骤14、根据岩心样本库，分别获取标准灯影组二段和标准灯影组四段的碳同位素和氧同位素的关系特征。

步骤14包括但不限于以下子步骤：

步骤141、根据岩心样本库，获取标准灯影组二段所对应的每个岩心样本的碳同位素的值和氧同位素的值。

根据《syt5238-2008有机物和碳酸盐岩碳、氧同位素分析方法》标准提供的方法获取标准灯影组二段所对应的每个岩心样本的碳同位素的值和氧同位素的值。并去掉异常数据点，碳同位素(δ¹³cpdb)小于-12‰，氧同位素(δ¹⁸opdb)小于-10‰或大于10‰的数值。

其中，δ¹³cpdb是将样品的c13稳定同位素值进行了国际标准pdb的处理，δ¹⁸opdb是将样品的o18稳定同位素值进行了国际标准pdb的处理。

步骤142、根据岩心样本库，获取标准灯影组四段所对应的每个岩心样本的碳同位素的值和氧同位素的值。

参见与步骤141相同的方法，获取标准灯影组四段所对应的每个岩心样本的碳同位素的值和氧同位素的值，并去掉异常数据点。

步骤143、以氧同位素的值为横坐标，碳同位素的值为纵坐标，获取标准灯影组二段所对应的每个岩心样本和标准灯影组四段所对应的每个岩心样本的散点图。

在散点图中，标准灯影组二段和标准灯影组四段所对应的每个岩心样本的散点形状可以不同，以能够容易地区分出灯影组二段和灯影组四段两个区域所对应的散点。

步骤15、获取目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征。

基于步骤12-步骤14同样的方法，获取目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征。

即获取目标岩心样本的元素与质量百分比之间的直方图、碳同位素种类与频率之间的直方图、碳同位素的值与氧同位素的值之间的散点图。

步骤16、若目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征均符合标准灯影组二段或标准灯影组四段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征，确定目标岩心样本属于灯影组二段或灯影组四段。

若目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征均与标准灯影组二段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征相匹配，则确定目标岩心样本属于灯影组二段。

若目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征均与标准灯影组四段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征相匹配，则确定目标岩心样本属于灯影组四段。

需要说明的是，“相匹配”指的是接近、类似、或在对应范围内。

对于步骤11和步骤15的顺序不作具体限定。基于步骤11-步骤16，能够实现对灯影组二段和灯影组四段的判别。

可选地，本发明实施例提供的方法还包括：

若目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征中至少一种不符合标准灯影组二段或标准灯影组四段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征，确定目标岩心样本不属于灯影组二段或灯影组四段。

此时，可能是岩心样本库不符合要求，需要重新筛选岩心样本库。

本发明实施例提供的用于判别灯影组二段和灯影组四段的方法，通过获取目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征，并与标准灯影组二段和标准灯影组四段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征进行比对，能够准确地、快速地、经济地确认目标岩心样本属于灯影组二段或灯影组四段，利于油气勘探工作的研究。

另一方面，本发明实施例提供了一种用于判别灯影组二段和灯影组四段的装置，如附图2所示，该装置包括：

第一获取模块21，用于获取标准灯影组二段和标准灯影组四段的岩心样本库。

第二获取模块22，用于根据岩心样本库，分别获取标准灯影组二段和标准灯影组四段的全岩元素分布特征。

第三获取模块23，用于根据岩心样本库，分别获取标准灯影组二段和标准灯影组四段的碳同位素分布频率特征。

第四获取模块24，用于根据岩心样本库，分别获取标准灯影组二段和标准灯影组四段的碳同位素和氧同位素的关系特征。

第五获取模块25，用于获取目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征。

第一确定模块26，用于若目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征均符合标准灯影组二段或标准灯影组四段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征，确定目标岩心样本属于灯影组二段或灯影组四段。

可选地，第二获取模块22，包括：

第一获取单元，用于根据岩心样本库，分别获取标准灯影组二段所对应的每种岩心样本的全元素及质量百分比；

第二获取单元，用于根据标准灯影组二段所对应的每种岩心样本的全元素及质量百分比，获取标准灯影组二段所对应的每种岩心样本的全元素及质量百分比之间的第一关系示意图；

第三获取单元，用于根据岩心样本库，分别获取标准灯影组四段所对应的每种岩心样本的全元素及质量百分比；

第四获取单元，用于根据标准灯影组四段所对应的每种岩心样本的全元素及质量百分比，获取标准灯影组四段所对应的每种岩心样本的全元素及质量百分比之间的第二关系示意图。

可选地，第三获取模块23，包括：

第五获取单元，用于根据岩心样本库，获取标准灯影组二段所对应的每种岩心样本的碳同位素的值及频率；

第六获取单元，用于根据标准灯影组二段所对应的每种岩心样本的碳同位素的值及频率，获取标准灯影组二段所对应的每种岩心样本的碳同位素的值及频率的第三关系示意图；

第七获取单元，用于根据岩心样本库，获取标准灯影组四段所对应的每种岩心样本的碳同位素的值及频率；

第八获取单元，用于根据标准灯影组四段所对应的每种岩心样本的碳同位素的值及频率，获取标准灯影组四段所对应的每种岩心样本的碳同位素的值及频率的第四关系示意图。

可选地，第四获取模块24，包括：

第九获取单元，用于根据岩心样本库，获取标准灯影组二段所对应的每个岩心样本的碳同位素的值和氧同位素的值；

第十获取单元，用于根据岩心样本库，获取标准灯影组四段所对应的每个岩心样本的碳同位素的值和氧同位素的值；

第十一获取单元，用于以氧同位素的值为横坐标，碳同位素的值为纵坐标，获取标准灯影组二段所对应的每个岩心样本和标准灯影组四段所对应的每个岩心样本的散点图。

可选地，本发明实施例提供的装置还包括：

第二确定模块，用于若目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征中至少一种不符合标准灯影组二段或标准灯影组四段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征，确定目标岩心样本不属于灯影组二段或灯影组四段。

本发明实施例提供的用于判别灯影组二段和灯影组四段的装置，通过第五获取模块获取目标岩心样本的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征，并与第二获取模块、第三获取模块、第四获取模块获取的标准灯影组二段和标准灯影组四段的全岩元素分布特征、碳同位素分布频率特征、碳同位素和氧同位素的关系特征进行比对，能够准确地、快速地、经济地确认目标岩心样本属于灯影组二段或灯影组四段，利于油气勘探工作的研究。

以下将通过一些实施例进一步地描述本发明。

实施例1

本实施例采用本发明实施例提供的方法对bhs1井底所处的层位进行判别，方法如下：

对无争议的灯影组二段、灯影组四段进行采样分析，避开孔洞发育的样品，选择颜色均一、致密的样品。灯影组二段对应的岩心样本包括：hs1井灯影组二段岩心样本30个、z6井灯影组二段岩心样本64个、mx8井灯影组二段岩心样本101个、w117井灯影组二段岩心样本174个。灯影组四段对应的岩心样本包括：gs7井灯影组四段岩心样本131个、mx9井灯影组四段岩心样本122个、mx8井灯影组四段岩心样本81个、w117井灯影组四段岩心样本11个。目标岩心样本包括102个bhs1井岩心样本。

分别获取hs1井灯影组二段岩心样本、z6井灯影组二段岩心样本、gs7井灯影组四段岩心样本、mx9井灯影组四段岩心样本、以及bhs1井目标岩心样本的全岩元素的质量百分比的平均值，参见表1。基于五种岩心样本中各元素的质量百分比的平均值，作出全元素与质量百分比之间的直方图，分别参见图3、图4、图5、图6、图7。

表1

由表1可知，灯影组四段中mg/ca的比值小于0.7，灯影组二段中mg/ca的比值大于0.7。灯影组四段中硅元素的质量百分比大于10％，灯影组二段中硅元素的质量百分比小于10％。灯影组四段中铝元素的质量百分比大于1％，灯影组二段中铝元素的质量百分比小于1％。而bhs1井的岩心样本中mg/ca的比值、硅元素的质量百分比、铝元素的质量百分比均与灯影组四段的特征相匹配。并且，由图3至图6可以看出，灯影组四段所对应的岩心样本富含钙、镁、硅，灯影组二段所对应的岩心样本富含钙、镁，而图7中bhs1井的岩心样本富含钙、镁、硅。基于上述，可以确定bhs1井的岩心样本在全岩元素分布特征方面接近灯影组四段。

根据岩心样本库，分别获取标准灯影组二段和标准灯影组四段的碳同位素分布频率特征。

分别获取hs1井灯影组二段岩心样本、z6井灯影组二段岩心样本、gs7井灯影组四段岩心样本、mx9井灯影组四段岩心样本、以及bhs1井目标岩心样本的碳同位素的平均值，具体参数详见表2。并获取碳同位素的频率，建立碳同位素的值与碳同位素的频率之间的直方图，参见附图8。

表2

由表2可知，灯影组四段的碳同位素的值均小于灯影组二段的碳同位素的值，bhs1井目标岩心样本的碳同位素的值较小，更接近灯影组四段。由附图8可知，bhs1井目标岩心样本的碳同位素的值的频率更接近灯影组四段。因此，bhs1井目标岩心样本在碳同位素分布频率特征方面与灯影组四段的碳同位素分布频率特征相匹配。

分别获取mx8井、w117井灯影组二段、灯影组四段、gs7灯影组四段每个岩心样本的碳同位素和氧同位素数据，以及bhs1井每个目标岩心样本中碳同位素和氧同位素数据。以碳同位素值作为纵坐标、氧同位素值作为横坐标生成散点图，根据已知样本建立灯影组二段与灯影组四段的分布确定判别图版，参见附图9。其中，灯影组二段可以简写为灯二段(d2)，灯影组四段可以简写为灯四段(d4)。

由附图9可以得出，灯影组二段与灯影组四段的散点特征差异明显，但仍存在小部分重叠区，而bhs1井目标岩心样本所对应的散点大部分落在了灯影组四段，少数落在重叠区，因此可以确定bhs1井目标岩心样本更接近灯影组四段。

综上，确定bhs1井目标岩心样本属于灯影组四段。

上述所有可选技术方案，可以采用任意结合形成本公开的可选实施例，在此不再一一赘述。

以上所述仅为本发明的说明性实施例，并不用以限制本发明的保护范围，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。