一种泥页岩成因分析方法及装置

摘要

本发明实施例提供一种泥页岩成因分析方法及装置。方法包括：获取待勘探区域中多口勘探井的测井数据，勘探井中包括至少一口岩心井和多口非岩心井；根据测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架；分别对每口勘探井进行标定，根据标定结果将岩心井的天文旋回地层格架和各非岩心井的天文旋回地层格架进行连接获得连接关系；获取岩心井对应的泥页岩成因，根据岩心井的泥页岩成因和连接关系获得所述非岩心井的泥页岩成因。装置用于执行上述方法。本发明实施例通过岩心井来分析非岩心井的泥页岩成因，无需对非岩心井进行取心分析，从而降低了分析成本。

权利要求

1.一种泥页岩成因分析方法，其特征在于，包括：获取待勘探区域中多口勘探井的测井数据，所述多口勘探井中包括至少一口岩心井和多口非岩心井；

根据所述测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架；

分别对每口勘探井进行标定，根据标定结果将所述岩心井的天文旋回地层格架和各非岩心井的天文旋回地层格架进行连接，获得连接关系；

获取岩心井对应的泥页岩成因，根据所述岩心井的泥页岩成因和所述连接关系获得所述非岩心井的泥页岩成因。

2.根据权利要求1所述的方法，其特征在于，所述方法还包括：获取所述岩心井的岩心资料，对所述岩心资料进行分析获得所述岩心井的泥页岩成因。

3.根据权利要求1所述的方法，其特征在于，所述根据所述测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架，包括：对所述测井数据进行预处理，所述预处理包括插值、去趋势化和异常值剔除；

对预处理后的测井数据进行频谱分析和显著性检验，获得不同时间尺度上所述测井数据中各种信号的信号值；

利用小波分析从所述测井数据中的多种信号中获取目标信号；

对所述目标信号进行轨道调谐分析，并根据分析结果建立所述勘探井的天文旋回地层格架。

4.根据权利要求1所述的方法，其特征在于，所述分别对每口勘探井进行标定，包括：利用预设测年方法对各勘探井进行测年，根据测年结果确定各勘探井的锚点，通过所述锚点对各勘探井进行标定；

所述预设测年方法为石英ESR测年、火山岩同位素测年或地磁极性测年。

5.根据权利要求1所述的方法，其特征在于，所述分别对每口勘探井进行标定，包括：根据所述勘探井的标志层、岩性特征、沉积旋回和电性特征对所述勘探井进行标定。

6.一种泥页岩成因分析装置，其特征在于，包括：获取模块，用于获取待勘探区域中多口勘探井的测井数据，所述多口勘探井中包括至少一口岩心井和多口非岩心井；

格架建立模块，用于根据所述测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架；

标定模块，用于分别对每口勘探井进行标定，根据标定结果将所述岩心井的天文旋回地层格架和各非岩心井的天文旋回地层格架进行连接，获得连接关系；

第一成因分析模块，用于获取岩心井对应的泥页岩成因，根据所述岩心井的泥页岩成因和所述连接关系获得所述非岩心井的泥页岩成因。

7.根据权利要求6所述的装置，其特征在于，所述装置，还包括：第二成因分析模块，用于获取所述岩心井的岩心资料，对所述岩心资料进行分析获得所述岩心井的泥页岩成因。

8.根据权利要求6所述的装置，其特征在于，所述格架建立模块，具体用于：对所述测井数据进行预处理，所述预处理包括插值、去趋势化和异常值剔除；

对预处理后的测井数据进行频谱分析和显著性检验，获得不同时间尺度上所述测井数据中各种信号的信号值；

利用小波分析从所述测井数据中的多种信号中获取目标信号；

对所述目标信号进行轨道调谐分析，并根据分析结果建立所述勘探井的天文旋回地层格架。

9.根据权利要求6-8任一项所述的装置，其特征在于，所述标定模块，具体用于：利用预设测年方法对各勘探井进行测年，根据测年结果确定各勘探井的锚点，通过所述锚点对各勘探井进行标定；

所述预设测年方法为石英ESR测年、火山岩同位素测年或地磁极性测年。

10.根据权利要求6-8任一项所述的装置，其特征在于，所述标定模块，具体用于：根据所述勘探井的标志层、岩性特征、沉积旋回和电性特征对所述勘探井进行标定。

说明书

一种泥页岩成因分析方法及装置

技术领域

[0001] 本发明涉及页岩油气勘探技术领域，具体而言，涉及一种泥页岩成因分析方法及装置。

背景技术

[0002] 在传统的砂岩或碳酸盐岩研究中，通常根据测井相来在无岩心的井中进行成因判断。即利用各类型沉积体在测井上的沉积旋回特征，建立各类型沉积旋回的特有测井相模板，从而利用测井相模板实现无岩心井的成因判定。和砂岩及碳酸盐岩垂向上频繁的岩性变化不同，泥页岩的岩性在垂向上较为稳定，这就导致泥页岩在垂向上的测井响应差异不明显，使目前还没有成熟的各成因类型泥页岩的测井相模板。因此，目前油气勘探开发中的泥页岩的成因分析依赖于岩心分析和测试，但该方法导致泥页岩成因分析只能在岩心井中开展，非岩心井的泥页岩成因没法分析。

[0003] 现有技术中，可以通过对泥页岩进行取心从而获得到取心资料，根据取心资料来分析泥页岩的成因。泥页岩取心成本昂贵，导致泥页岩取心资料稀缺。在钻井过程中，如需取心，需要先暂停钻井，然后用专门的设备取心。这将迅速增加钻井的周期和成本。据粗略估计，每取1米岩心的成本为2-5万元。对于油气勘探和有利层段评价而言，需要在明确厚达数十米甚至数百米岩心成因的基础上进行有利层段评价。因此，如仅依赖于岩心，将极大地增加油气勘探成本。

发明内容

[0004] 有鉴于此，本发明实施例的目的在于提供一种泥页岩成因分析方法及装置，以解决上述技术问题。

[0005] 第一方面，本发明实施例提供了一种泥页岩成因分析方法，包括：

[0006] 获取待勘探区域中多口勘探井的测井数据，所述多口勘探井中包括至少一口岩心井和多口非岩心井；

[0007] 根据所述测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架；

[0008] 分别对每口勘探井进行标定，根据标定结果将所述岩心井的天文旋回地层格架和各非岩心井的天文旋回地层格架进行连接，获得连接关系；

[0009] 获取岩心井对应的泥页岩成因，根据所述岩心井的泥页岩成因和所述连接关系获得所述非岩心井的泥页岩成因。

[0010] 进一步地，所述方法还包括：

[0011] 获取所述岩心井的岩心资料，对所述岩心资料进行分析获得所述岩心井的泥页岩成因。

[0012] 进一步地，所述根据所述测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架，包括：

[0013] 对所述测井数据进行预处理，所述预处理包括插值、去趋势化和异常值剔除；

[0014] 对预处理后的测井数据进行频谱分析和显著性检验，获得不同时间尺度上所述测井数据中各种信号的信号值；

[0015] 利用小波分析从所述测井数据中的多种信号中获取目标信号；

[0016] 对所述目标信号进行轨道调谐分析，并根据分析结果建立所述勘探井的天文旋回地层格架。

[0017] 进一步地，所述分别对每口勘探井进行标定，包括：

[0018] 利用预设测年方法对各勘探井进行测年，根据测年结果确定各勘探井的锚点，通过所述锚点对各勘探井进行标定；

[0019] 所述预设测年方法包括石英ESR测年、火山岩同位素测年或地磁极性测年。

[0020] 进一步地，所述分别对每口勘探井进行标定，包括：

[0021] 根据所述勘探井的标志层、岩性特征、沉积旋回和电性特征对所述勘探井进行标定。

[0022] 第二方面，本发明实施例提供了一种泥页岩成因分析装置，包括：

[0023] 获取模块，用于获取待勘探区域中多口勘探井的测井数据，所述多口勘探井中包括至少一口岩心井和多口非岩心井；

[0024] 格架建立模块，用于根据所述测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架；

[0025] 标定模块，用于分别对每口勘探井进行标定，根据标定结果将所述岩心井的天文旋回地层格架和各非岩心井的天文旋回地层格架进行连接，获得连接关系；

[0026] 第一成因分析模块，用于获取岩心井对应的泥页岩成因，根据所述岩心井的泥页岩成因和所述连接关系获得所述非岩心井的泥页岩成因。

[0027] 进一步地，所述装置，还包括：

[0028] 第二成因分析模块，用于获取所述岩心井的岩心资料，对所述岩心资料进行分析获得所述岩心井的泥页岩成因。

[0029] 进一步地，所述格架建立模块，具体用于：

[0030] 对所述测井数据进行预处理，所述预处理包括插值、去趋势化和异常值剔除；

[0031] 对预处理后的测井数据进行频谱分析和显著性检验，获得不同时间尺度上所述测井数据中各种信号的信号值；

[0032] 利用小波分析从所述测井数据中的多种信号中获取目标信号；

[0033] 对所述目标信号进行轨道调谐分析，并根据分析结果建立所述勘探井的天文旋回地层格架。

[0034] 进一步地，所述标定模块，具体用于：

[0035] 利用预设测年方法对各勘探井进行测年，根据测年结果确定各勘探井的锚点，通过所述锚点对各勘探井进行标定；

[0036] 所述预设测年方法包括石英ESR测年、火山岩同位素测年或地磁极性测年。

[0037] 进一步地，所述标定模块，具体用于：

[0038] 根据所述勘探井的标志层、岩性特征、沉积旋回和电性特征对所述勘探井进行标定。

[0039] 第三方面，本发明实施例提供一种电子设备，包括：处理器、存储器和总线，其中，[0040] 所述处理器和所述存储器通过所述总线完成相互间的通信；

[0041] 所述存储器存储有可被所述处理器执行的程序指令，所述处理器调用所述程序指令能够执行第一方面的方法步骤。

[0042] 第四方面，本发明实施例提供一种非暂态计算机可读存储介质，包括：

[0043] 所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行第一方面的方法步骤。

[0044] 本发明实施例通过测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架，然后将岩心井和非岩心井的天文旋回地层格架进行连接，并根据岩心井的泥页岩成因获得非岩心井的泥页岩成因，无需对非岩心井进行取心分析，从而降低了分析成本。

[0045] 本发明的其他特征和优点将在随后的说明书阐述，并且，部分地从说明书中变得显而易见，或者通过实施本发明实施例了解。本发明的目的和其他优点可通过在所写的说明书、权利要求书、以及附图中所特别指出的结构来实现和获得。

附图说明

[0046] 为了更清楚地说明本发明实施例的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，应当理解，以下附图仅示出了本发明的某些实施例，因此不应被看作是对范围的限定，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他相关的附图。

[0047] 图1为本发明实施例提供的一种泥页岩成因分析方法流程示意图；

[0048] 图2为本发明实施例提供的东营凹陷B井频谱分析图；

[0049] 图3为本发明实施例提供的东营凹陷B井浮动天文旋回格架示意图；

[0050] 图4为本发明实施例提供的A井和B井的天文旋回地层格架连接示意图；

[0051] 图5为本发明实施例提供的一种泥页岩成因分析装置结构示意图；

[0052] 图6为本发明实施例提供的电子设备的结构框图。

具体实施方式

[0053] 下面将结合本发明实施例中附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。通常在此处附图中描述和示出的本发明实施例的组件可以以各种不同的配置来布置和设计。因此，以下对在附图中提供的本发明的实施例的详细描述并非旨在限制要求保护的本发明的范围，而是仅仅表示本发明的选定实施例。基于本发明的实施例，本领域技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。

[0054] 应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。同时，在本发明的描述中，术语“第一”、“第二”等仅用于区分描述，而不能理解为指示或暗示相对重要性。

[0055] 图1为本发明实施例提供的一种泥页岩成因分析方法流程示意图，如图1所示，该方法包括：

[0056] 步骤101：获取待勘探区域中多口勘探井的测井数据，所述多口勘探井中包括至少一口岩心井和多口非岩心井。

[0057] 在具体的实施过程中，由于天文旋回参数的作用范围大，在一个盆地及其周缘地区日晒量基本相同，因此对盆地深水区各区域的泥页岩的影响作用相同。基于此，在天文旋回作用下，同一时间内的深水区相距较远(数十公里到数百公里)的不同井沉积的泥页岩具有相同的成因机制。因此，可以获取待勘探区域中多口勘探井的测井数据，对每口勘探井进行GR测井分析，保证，GR曲线的测试密度不大于0.125米，能够满足旋回地层学分析的精度要求。应当说明的是，在多口勘探井中，应当有至少一口岩心井，剩下的勘探井为非岩心井，对于岩心井来说，还可以获取到该岩心井对应的岩心资料。

[0058] 步骤102：根据所述测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架。

[0059] 在具体的实施过程中，分别利用测井数据对相应勘探井进行米兰科维奇旋回分析。米兰科维奇旋回有三个参数，分别是：1)偏心率(周期约为0.1、0.4、1.3、2百万年)，表征地球绕太阳运动的轨迹。2)斜率(周期为39、41千年)，即地球自转轴的倾角(变化范围22o02'-24o30')，表现为黄道面与赤道面的夹角。3)岁差(周期为19、23千年)，指地球自转轴绕黄道轴的旋转运动。不同级别天文旋回通过以下机理控制沉积作用。偏心率值增大或减小控制地球与太阳之间的距离，偏心率值较高阶段比偏心率值较低阶段日晒量高20％左右。斜率的大小影响到地球上不同纬度和不同季节的气候差异程度大小。比如，如果斜率为零，太阳将一直直射赤道，地球上不存在季节变化；斜率增加，季节性增强。在受偏心率调制的情况下，一个岁差周期内的北半球(或南半球)近日点的季节经历一个“夏季-冬季-夏季”的周期变化，从而影响太阳能量的分布。以上参数一起作用于盆地，共同控制泥页岩沉积作用。因此，根据对每口勘探井的分析结果可以建立天文旋回地层格架。

[0060] 步骤103：分别对每口勘探井进行标定，根据标定结果将所述岩心井的天文旋回地层格架和各非岩心井的天文旋回地层格架进行连接，获得连接关系。

[0061] 在具体的实施过程中，分别对每口勘探井进行标定，应当说明的是，可以根据测年来确定勘探井的锚点，然后通过锚点对勘探井进行标定，也可以通过选择全区稳定分布的标志层进行标定，根据标定结果将所述岩心井的天文旋回地层格架和各非岩心井的天文旋回地层格架进行连接，获得连接关系。

[0062] 步骤104：获取岩心井对应的泥页岩成因，根据所述岩心井的泥页岩成因和所述连接关系获得所述非岩心井的泥页岩成因。

[0063] 在具体的实施过程中，通过岩心井的岩心资料可以获取该岩心井对应的泥页岩成因，然后分析岩心井的取心层段中各成因机制泥页岩与天文旋回格架的耦合规律，在明确泥页岩岩心沉积机制及其在天文旋回格架下的演化规律基础上，即可在横跨岩心井与非岩心井天文旋回格架的约束下，分析非岩心井的泥页岩沉积机制。可以理解的是，由于天文旋回参数的作用范围大，在一个盆地及其周缘地区日晒量基本相同，因此对盆地深水区各区域的泥页岩的影响作用相同。基于此，在天文旋回作用下，同一时间内的深水区相距较远(数十公里到数百公里)的不同井沉积的泥页岩具有相同的成因机制。因此，可以通过分析岩心井中泥页岩的成因，进而在同一层段内推广到邻近的无岩心的井中。最终实现在无岩心的井中，通过泥页岩的成因机制来优选有利储集层段。因此，通过岩心井的泥页岩成因以及岩心井与非岩心井之间的天文旋回地层格架的连接关系，可以对应获得到非岩心井的泥页岩成因。

[0064] 本发明实施例通过测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架，然后将岩心井和非岩心井的天文旋回地层格架进行连接，并根据岩心井的泥页岩成因获得非岩心井的泥页岩成因，无需对非岩心井进行取心分析，从而降低了分析成本。

[0065] 在上述实施例的基础上，所述方法还包括：

[0066] 获取所述岩心井的岩心资料，对所述岩心资料进行分析获得所述岩心井的泥页岩成因。

[0067] 在具体的实施过程中，通过岩心井中的取心层的岩心资料，可以对该岩心井的泥页岩成因进行分析，其具体分析方法与现有传统方法一致，本发明实施例对此不做具体限定。

[0068] 本发明实施例根据岩心资料分析岩心井的泥页岩成因，并对岩心井和非岩心井建立天文旋回地层格架的连接关系，从而根据岩心井的泥页岩成因以及连接关系获得非岩心井的泥页岩成因，无需对非岩心井进行取心操作便可获取到泥页岩成因，降低了分析成本。

[0069] 在上述实施例的基础上，所述根据所述测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架，包括：

[0070] 对所述测井数据进行预处理，所述预处理包括插值、去趋势化和异常值剔除；

[0071] 对预处理后的测井数据进行频谱分析和显著性检验，获得不同时间尺度上所述测井数据中各种信号的信号值；

[0072] 利用小波分析从所述测井数据中的多种信号中获取目标信号；

[0073] 对所述目标信号进行轨道调谐分析，并根据分析结果建立所述勘探井的天文旋回地层格架。

[0074] 在具体的实施过程中，建立各岩心井和非岩心井的天文旋回格架，细分为以下4个步骤，分别是：

[0075] 1)数据选择和预处理：含油气盆地测井数据采样间隔一般为0.125m，能够满足旋回地层学分析的精度要求，故这里采用GR曲线进行米兰科维奇旋回分析。在开展高频旋回研究之前，利用Forward测井解释平台的环境校正模块对GR数据进行环境影响校正，消除其他非轨道因素对结果可能造成的影响。同时对GR数据进行插值和去趋势化处理，并去掉异常极值，尽量使轨道因素影响的结果显著。

[0076] 2)频谱分析和显著性检验：将一个时间序列分解为不同时间尺度上的贡献值，更好的理解时间序列中各种信号变化性产生的物理过程。在频谱分析后，应用自噪声模型或ARI噪声模型对结果进行显著性检验。本发明实施例选择0.05的显著性水平，也就是高于95％的置信度。通过将位于置信度之上的计算结果与理论周期进行比较，确定实际周期。

[0077] 3)小波分析：采用Torrence(1998)提供的连续小波变换工具箱，利用小波分析来分析时间序列中目标信号变化，通过将时间序列分解到时频空间，确定目标信号的主要状态以及这些状态随时间的变化情况。

[0078] 4)轨道调谐：首先以地球轨道要素的理论计算变化曲线为目标曲线，其次对时间控制点之间的古气候数据插值以构建一个初始时间尺度，所用数据由独立的地质年代研究所得，再次反复应用数字滤波从古气候记录中提取相应的轨道频率组分，将目标信号的曲线与目标曲线进行对比分析，如不吻合，则通过添加新的时间控制点使滤波曲线与相应的理论轨道曲线协调一致，从而保持气候记录与地球轨道要素的变化在周期和相位上一致。

[0079] 基于以上分析可以建立勘探井的天文旋回地层格架。

[0080] 在上述实施例的基础上，所述分别对每口勘探井进行标定，包括：

[0081] 利用预设测年方法对各勘探井进行测年，根据测年结果确定各勘探井的锚点，通过所述锚点对各勘探井进行标定；

[0082] 所述预设测年方法包括石英ESR测年、火山岩同位素测年或地磁极性测年。

[0083] 在具体的实施过程中，为了能够获得岩心井的天文旋回地层格架与非岩心井之间的天文旋回地层格架之间的关系，在分别为勘探井建立了天文旋回地层格架之后，还需要对每个勘探井进行标定，从而能够获取到岩心井与非岩心井之间的对应关系。其中可以通过预设测年方法对各个勘探井进行测年，根据测年结果确定各勘探井的锚点，通过锚点对各个勘探井进行标定。具体的预设测年方法可以是：

[0084] 1)选择经风化、破碎和再沉积以后才开始形成氧空位和E'心的石英，进行ESR石英测年。

[0085] 2)火山岩或火山灰同位素测年：根据不同的研究对象，采用不同的同位素测年方法，锆石已成为同位素年代学研究中最常用和最有效的方法之一，而很老的岩石年代学可以用K-Ar法测定。

[0086] 3)地磁极性年龄：利用最新的Lourens04年表，测定岩石年龄。

[0087] 本发明实施例通过锚点来对各个勘探井进行标定，从而能够较为准确的获得岩心井的天文旋回地层格架和非岩心井的天文旋回地层格架之间的连接关系。

[0088] 在上述实施例的基础上，所述分别对每口勘探井进行标定，包括：

[0089] 根据所述勘探井的标志层、岩性特征、沉积旋回和电性特征对所述勘探井进行标定。

[0090] 在具体的实施过程中，为了能够获得岩心井的天文旋回地层格架与非岩心井之间的天文旋回地层格架之间的关系，在分别为勘探井建立了天文旋回地层格架之后，还需要对每个勘探井进行标定，从而能够获取到岩心井与非岩心井之间的对应关系。其中可以选择各勘探井的标志层，标志层岩性稳定、厚度均匀且易于识别。以标志层为基础，再结合地层的岩性特征，沉积旋回和电性特征等，详细进行地层的对比和划分。这里主要采用测井曲线解释法，选择全井段：①在全油区分布广泛且稳定；②储层的岩性和物性特征明显；③在全油区范围内，储层的测井特征明显且无明显变化或变化有规律。

[0091] 本发明实施例通过标志仓来对各个勘探井进行标定，从而能够较为准确的获得岩心井的天文旋回地层格架和非岩心井的天文旋回地层格架之间的连接关系。

[0092] 本发明实施例通过以中国东营凹陷始新统的A井(非取心井)和B井(取心井)为例，对A井的泥页岩成因进行分析：

[0093] 首先，开展这两口井的钻井和测井工作，并获取B井在深水泥页岩层段的岩心资料。

[0094] 其次，分别对这两口井进行米兰科维奇旋回分析，具体包括：选取GR数据进行预处理，然后进行频谱分析，在95％的置信度范围进行显著性检验，其次进行小波分析，图2为本发明实施例提供的东营凹陷B井频谱分析图，如图2所示。结合小波能量谱与频谱分析可知，小波图上的4个能量高的频带与频谱分析的主要周期有较好的对应，如下表所示：

[0095]

[0096] 在完成数据选择和预处理、频谱分析和显著性检验、小波分析、轨道调谐分析后，分别建立A井和B井的浮动旋回地层格架，图3为本发明实施例提供的东营凹陷B井浮动天文旋回格架示意图，如图3所示。B井泥页岩取心段可细分为5个125千年的偏心率旋回(自下而上依次编号1-5)和27个23千年的岁差旋回(自下而上依次编号1-27)。

[0097] 然后，在A井和B井中的锚点年龄44.6Ma的标定下，即“连接”了A井和B井的天文旋回格架。分别对A井的偏心率旋回和岁差旋回编号，保证两口井同一时期沉积旋回的编号相同，图4为本发明实施例提供的A井和B井的天文旋回地层格架连接示意图，如图4所示。

[0098] 最后，分析B井泥页岩的取心特征，垂向上依次划分出多期静水泥页岩和重力流沉积泥页岩。分析B井天文旋回与泥页岩沉积机制的耦合关系，发现125千年偏心率旋回的最大值附近的三个23千年岁差旋回对应于重力流泥页岩，而其它时期的泥页岩为静水沉积，从而明确了各岁差旋回的泥页岩沉积成因。因此，A井与B井相对应的岁差旋回也具有同样的沉积成因，从而判定出A井的泥页岩沉积成因。

[0099] 图5为本发明实施例提供的一种泥页岩成因分析装置结构示意图，如图5所示，该装置包括：获取模块501、格架建立模块502、标定模块503和第一成因分析模块504，其中，[0100] 获取模块501用于获取待勘探区域中多口勘探井的测井数据，所述多口勘探井中包括至少一口岩心井和多口非岩心井；格架建立模块502用于根据所述测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架；标定模块503用于分别对每口勘探井进行标定，根据标定结果将所述岩心井的天文旋回地层格架和各非岩心井的天文旋回地层格架进行连接，获得连接关系；第一成因分析模块504用于获取岩心井对应的泥页岩成因，根据所述岩心井的泥页岩成因和所述连接关系获得所述非岩心井的泥页岩成因。

[0101] 在上述实施例的基础上，所述装置，还包括：

[0102] 第二成因分析模块，用于获取所述岩心井的岩心资料，对所述岩心资料进行分析获得所述岩心井的泥页岩成因。

[0103] 在上述实施例的基础上，所述格架建立模块，具体用于：

[0104] 对所述测井数据进行预处理，所述预处理包括插值、去趋势化和异常值剔除；

[0105] 对预处理后的测井数据进行频谱分析和显著性检验，获得不同时间尺度上所述测井数据中各种信号的信号值；

[0106] 利用小波分析从所述测井数据中的多种信号中获取目标信号；

[0107] 对所述目标信号进行轨道调谐分析，并根据分析结果建立所述勘探井的天文旋回地层格架。

[0108] 在上述实施例的基础上，所述标定模块，具体用于：

[0109] 利用预设测年方法对各勘探井进行测年，根据测年结果确定各勘探井的锚点，通过所述锚点对各勘探井进行标定；

[0110] 所述预设测年方法包括石英ESR测年、火山岩同位素测年或地磁极性测年。

[0111] 在上述实施例的基础上，所述标定模块，具体用于：

[0112] 根据所述勘探井的标志层、岩性特征、沉积旋回和电性特征对所述勘探井进行标定。

[0113] 所属领域的技术人员可以清楚地了解到，为描述的方便和简洁，上述描述的装置的具体工作过程，可以参考前述方法中的对应过程，在此不再过多赘述。

[0114] 综上所述，本发明实施例通过测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架，然后将岩心井和非岩心井的天文旋回地层格架进行连接，并根据岩心井的泥页岩成因获得非岩心井的泥页岩成因，无需对非岩心井进行取心分析，从而降低了分析成本。

[0115] 本实施例公开一种计算机程序产品，所述计算机程序产品包括存储在非暂态计算机可读存储介质上的计算机程序，所述计算机程序包括程序指令，当所述程序指令被计算机执行时，计算机能够执行上述各方法实施例所提供的方法，例如包括：获取待勘探区域中多口勘探井的测井数据，所述多口勘探井中包括至少一口岩心井和多口非岩心井；根据所述测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架；分别对每口勘探井进行标定，根据标定结果将所述岩心井的天文旋回地层格架和各非岩心井的天文旋回地层格架进行连接，获得连接关系；获取岩心井对应的泥页岩成因，根据所述岩心井的泥页岩成因和所述连接关系获得所述非岩心井的泥页岩成因。

[0116] 本实施例提供一种非暂态计算机可读存储介质，所述非暂态计算机可读存储介质存储计算机指令，所述计算机指令使所述计算机执行上述各方法实施例所提供的方法，例如包括：获取待勘探区域中多口勘探井的测井数据，所述多口勘探井中包括至少一口岩心井和多口非岩心井；根据所述测井数据对相应的勘探井进行米兰科维奇旋回分析，并根据分析结果对各勘探井建立天文旋回地层格架；分别对每口勘探井进行标定，根据标定结果将所述岩心井的天文旋回地层格架和各非岩心井的天文旋回地层格架进行连接，获得连接关系；获取岩心井对应的泥页岩成因，根据所述岩心井的泥页岩成因和所述连接关系获得所述非岩心井的泥页岩成因。

[0117] 请参照图6，图6为本发明实施例提供的电子设备的结构框图。电子设备可以包括分析装置601、存储器602、存储控制器603、处理器604、外设接口605、输入输出单元606、音频单元607、显示单元608。

[0118] 所述存储器602、存储控制器603、处理器604、外设接口605、输入输出单元606、音频单元607、显示单元608各元件相互之间直接或间接地电性连接，以实现数据的传输或交互。例如，这些元件相互之间可通过一条或多条通讯总线或信号线实现电性连接。所述分析装置601包括至少一个可以软件或固件(firmware)的形式存储于所述存储器602中或固化在分析装置601的操作系统(operating system，OS)中的软件功能模块。所述处理器604用于执行存储器602中存储的可执行模块，例如分析装置601包括的软件功能模块或计算机程序。

[0119] 其中，存储器602可以是，但不限于，随机存取存储器(Random Access Memory，RAM)，只读存储器(Read Only Memory，ROM)，可编程只读存储器(Programmable Read-Only Memory，PROM)，可擦除只读存储器(Erasable Programmable Read-Only Memory，EPROM)，电可擦除只读存储器(Electric Erasable Programmable Read-Only Memory，EEPROM)等。其中，存储器602用于存储程序，所述处理器604在接收到执行指令后，执行所述程序，前述本发明实施例任一实施例揭示的流过程定义的服务器所执行的方法可以应用于处理器604中，或者由处理器604实现。

[0120] 处理器604可以是一种集成电路芯片，具有信号的处理能力。上述的处理器604可以是通用处理器，包括中央处理器(Central Processing Unit，简称CPU)、网络处理器(Network Processor，简称NP)等；还可以是数字信号处理器(DSP)、专用集成电路(ASIC)、现成可编程门阵列(FPGA)或者其他可编程逻辑器件、分立门或者晶体管逻辑器件、分立硬件组件。可以实现或者执行本发明实施例中的公开的各方法、步骤及逻辑框图。通用处理器可以是微处理器或者该处理器604也可以是任何常规的处理器等。

[0121] 所述外设接口605将各种输入/输出装置耦合至处理器604以及存储器602。在一些实施例中，外设接口605，处理器604以及存储控制器603可以在单个芯片中实现。在其他一些实例中，他们可以分别由独立的芯片实现。

[0122] 输入输出单元606用于提供给用户输入数据实现用户与所述服务器(或本地终端)的交互。所述输入输出单元606可以是，但不限于，鼠标和键盘等。

[0123] 音频单元607向用户提供音频接口，其可包括一个或多个麦克风、一个或者多个扬声器以及音频电路。

[0124] 显示单元608在所述电子设备与用户之间提供一个交互界面(例如用户操作界面)或用于显示图像数据给用户参考。在本实施例中，所述显示单元608可以是液晶显示器或触控显示器。若为触控显示器，其可为支持单点和多点触控操作的电容式触控屏或电阻式触控屏等。支持单点和多点触控操作是指触控显示器能感应到来自该触控显示器上一个或多个位置处同时产生的触控操作，并将该感应到的触控操作交由处理器604进行计算和处理。

[0125] 所述外设接口605将各种输入/输入装置耦合至处理器604以及存储器602。在一些实施例中，外设接口605，处理器604以及存储控制器603可以在单个芯片中实现。在其他一些实例中，他们可以分别由独立的芯片实现。

[0126] 输入输出单元606用于提供给用户输入数据实现用户与处理终端的交互。所述输入输出单元606可以是，但不限于，鼠标和键盘等。

[0127] 可以理解，图6所示的结构仅为示意，所述电子设备还可包括比图6中所示更多或者更少的组件，或者具有与图6所示不同的配置。图6中所示的各组件可以采用硬件、软件或其组合实现。

[0128] 在本申请所提供的几个实施例中，应该理解到，所揭露的装置和方法，也可以通过其它的方式实现。以上所描述的装置实施例仅仅是示意性的，例如，附图中的流程图和框图显示了根据本发明的多个实施例的装置、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或代码的一部分，所述模块、程序段或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意，在有些作为替换的实现方式中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。

[0129] 另外，在本发明各个实施例中的各功能模块可以集成在一起形成一个独立的部分，也可以是各个模块单独存在，也可以两个或两个以上模块集成形成一个独立的部分。

[0130] 所述功能如果以软件功能模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤。而前述的存储介质包括：U盘、移动硬盘、只读存储器(ROM，Read-Only Memory)、随机存取存储器(RAM，Random Access Memory)、磁碟或者光盘等各种可以存储程序代码的介质。

[0131] 以上所述仅为本发明的优选实施例而已，并不用于限制本发明，对于本领域的技术人员来说，本发明可以有各种更改和变化。凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步定义和解释。

[0132] 以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应所述以权利要求的保护范围为准。

[0133] 需要说明的是，在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、物品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、物品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括所述要素的过程、方法、物品或者设备中还存在另外的相同要素。