一种地层层序划分方法及装置

摘要

本发明提供了一种地层层序划分方法及装置，方法包括：根据地震剖面的地震反射特征识别二级层序界面、三级层序界面；将所述二级层序界面、三级层序界面分别定位到地震剖面中的一测井的对应深度；在二级层序界面的对应深度点处对测井的测井曲线进行一级重构，生成一级重构曲线；在三级层序界面的对应深度点处对一级重构曲线进行二级重构，生成二级重构曲线；按预设尺度对二级重构曲线进行小波分析，生成地层层序划分结果。本发明对长期层序的准确识别保证了层序曲线中长期低频层序和短期高频层序划分的结果准确可靠，依据地震对测井曲线经过两次重构得到二级重构曲线，进行不同尺度的小波变换分析，实现了层序划分过程中真正的井震结合。

权利要求

1.一种地层层序划分方法，其特征在于，所述的方法包括：

根据地震剖面的地震反射特征识别二级层序界面、三级层序界面；

将所述二级层序界面、三级层序界面分别定位到所述地震剖面中的一测井的对应深度；

在所述二级层序界面的对应深度点处对所述测井的测井曲线进行一级重构，生成一级重构曲线；

在所述三级层序界面的对应深度点处对所述一级重构曲线进行二级重构，生成二级重构曲线；

按预设尺度对所述二级重构曲线进行小波分析，生成地层层序划分结果；

其中，所述的将所述二级层序界面、三级层序界面分别定位到所述地震剖面中的一测井的对应深度包括：根据所述地震剖面中的测井的合成记录的时深关系分别将所述二级层序界面时间点t1，t2，…，tn对应到所述测井的测井曲线上的深度点h1，h2，…，hn；

根据所述合成记录的时深关系将所述三级层序界面的时间点tn1，tn2，…，tnm对应到测井曲线上的深度点hn1，hn2，…，hnm；

其中，tnm为第n个二级层序内第m个三级层序界面的时间点，n表示二级层序的次序，m表示三级层序界面的次序，hnm为第n个二级层序内第m个三级层序界面的深度点，n表示二级层序的次序，m表示三级层序界面的次序，m、n∈N+，N+为正整数；

其中，所述的在所述二级层序界面的对应深度点处对所述测井的测井曲线进行一级重构，生成的一级重构曲线为：SB(h)＝SA(h)+SAmax/8+(SAmax/8)*(-1)n+1，hn

2.如权利要求1所述的地层层序划分方法，其特征在于，所述的在所述三级层序界面的对应深度点处对所述一级重构曲线进行二级重构，生成的二级重构曲线为：SC(h)＝SB(h)-SAmax(n)/16+(SAmax(n)/16)*(-1)m，hnm

3.如权利要求2所述的地层层序划分方法，其特征在于，所述的按预设尺度对所述二级重构曲线进行小波分析，生成地层层序划分结果包括：按预设尺度对所述二级重构曲线进行小波分析，生成小波能量谱系图和小波系数曲线图；

根据所述小波能量谱系图和小波系数曲线图的间断点对所述二级重构曲线进行层序划分。

4.如权利要求3所述的地层层序划分方法，其特征在于，所述按预设尺度对所述二级重构曲线进行小波分析，生成地层层序划分结果进一步包括：对所述二级重构曲线进行6000尺度的小波分析，生成6000尺度的小波能量谱系图，根据该小波能量谱系图中的间断点划分二级层序；

对所述二级重构曲线进行1400尺度的小波分析，生成1400尺度的小波能量谱系图，根据该小波能量谱系图中的间断点划分三级层序；

对所述二级重构曲线进行700尺度的小波分析，生成700尺度的小波能量谱系图，根据该小波能量谱系图中的间断点划分四级层序；

对所述二级重构曲线进行150尺度的小波分析，生成小波系数曲线图，根据小波系数曲线图中的间断点划分五级层序。

5.一种地层层序划分装置，其特征在于，所述的装置包括：

识别模块，用于根据地震剖面的地震反射特征识别二级层序界面、三级层序界面；

定位模块，用于将所述二级层序界面、三级层序界面分别定位到所述地震剖面中的一测井的对应深度；

一级重构模块，用于在所述二级层序界面的对应深度点处对所述测井的测井曲线进行一级重构，生成一级重构曲线；

二级重构模块，在所述三级层序界面的对应深度点处对所述一级重构曲线进行二级重构，生成二级重构曲线；

地层划分模块，用于按预设尺度对所述二级重构曲线进行小波分析，生成地层层序划分结果；

其中，所述的定位模块包括：

二级层序定位单元，根据所述地震剖面中的测井的合成记录的时深关系分别将所述二级层序界面时间点t1，t2，…，tn对应到所述测井的测井曲线上的深度点h1，h2，…，hn；

三级层序定位单元，根据所述合成记录的时深关系将所述三级层序界面的时间点tn1，tn2，…，tnm对应到测井曲线上的深度点hn1，hn2，…，hnm；

其中，tnm为第n个二级层序内第m个三级层序界面的时间点，n表示二级层序的次序，m表示三级层序界面的次序，hnm为第n个二级层序内第m个三级层序界面的深度点，n表示二级层序的次序，m表示三级层序界面的次序，m、n∈N+，N+为正整数；

其中，所述的一级重构曲线为：

SB(h)＝SA(h)+SAmax/8+(SAmax/8)*(-1)n+1，hn

6.如权利要求5所述的地层层序划分装置，其特征在于，所述的二级重构曲线为：SC(h)＝SB(h)-SAmax(n)/16+(SAmax(n)/16)*(-1)m，hnm

7.如权利要求6所述的地层层序划分装置，其特征在于，所述的地层划分模块包括：小波分析单元，用于按预设尺度对所述二级重构曲线进行小波分析，生成小波能量谱系图和小波系数曲线图；

划分单元，用于根据所述小波能量谱系图和小波系数曲线图的间断点对所述二级重构曲线进行层序划分。

8.如权利要求7所述的地层层序划分装置，其特征在于，所述小波分析单元进行小波分析包括：对所述二级重构曲线进行6000尺度的小波分析，生成6000尺度的小波能量谱系图，根据该小波能量谱系图中的间断点划分二级层序；

对所述二级重构曲线进行1400尺度的小波分析，生成1400尺度的小波能量谱系图，根据该小波能量谱系图中的间断点划分三级层序；

对所述二级重构曲线进行700尺度的小波分析，生成700尺度的小波能量谱系图，根据该小波能量谱系图中的间断点划分四级层序；

对所述二级重构曲线进行150尺度的小波分析，生成小波系数曲线图，根据小波系数曲线图中的间断点划分五级层序。

说明书

一种地层层序划分方法及装置

技术领域

[0001] 本发明涉及勘探技术，具体的讲是一种地层层序划分方法及装置。

背景技术

[0002] 层序地层学在油气田勘探开发阶段中已得到广泛的应用，如何对层序进行准确划分是层序地层学研究中最基础、最重要的一个环节。准确的层序划分是进行沉积旋回分析，沉积环境研究的基础，对地层的充填模式及沉积相的研究具有重要意义。层序是一套相对整一的、成因上有联系的、其顶和底面以不整合面或者与这些不整合面可以对比的整合面为界的地层，是反映自然界周期性海平面升降的结果。层序反映着沉积旋回的变换，长周期的沉积旋回在地震剖面上有更明显的响应，而短周期的沉积旋回则在测井曲线上有更敏感的特征。因此，要想准确对层序进行划分，必须结合地震和测井资料进行综合分析。

[0003] 随着计算机技术的发展，层序划分的方法也越来越多，现有技术中以小波变换、频谱属性趋势分析和砂泥比曲线分析为主。

[0004] 小波变换是一种时频分析方法。利用小波变换划分层序，由多个不同周期(尺度)沉积旋回叠加的测井曲线，通过小波变换，被分解成各自周期独立的沉积旋回，以尺度的形式展示出来。研究表明，由多个不同周期沉积旋回叠加的测井曲线，通过考察小波时频能量图局部能量团的变化和多种伸缩尺度的周期性震荡特征，可分析地层的旋回性并与各级层序界面建立对应关系。例如，同一套地层在小波能谱图上表现为稳定的能量团，能量团间的断点以及小波曲线的变化点对应地层中的沉积间歇面或剥蚀面，表征沉积环境的改变，因此可作为层序界面，这就是小波时频分析划分地层的依据。

[0005] 频谱属性趋势分析的主要研究方法是以最大熵方法为基础对测井资料进行最大熵频谱属性分析，其实质就是根据研究目的不同把包含有多种频谱成分的测井曲线利用滑动窗口技术分解成单一频谱成分的资料，从而得到频谱变化属性(PEFA)曲线，利用PEFA技术分析测井资料，可以快速识别地层间断即层序界面。对PEFA曲线进行积分运算得到更有价值的INPEFA曲线，应用其拐点识别层序界面和洪泛面。

[0006] 运用砂泥比曲线划分层序，根据测井曲线采用计算机自动划分高分辨率层序时，首先通过自然伽马测井曲线计算出泥砂比、砂泥比曲线，然后综合利用这些曲线求出地层短期基准面变化曲线，最后根据短期基准面变化曲线，综合其他信息划分出不同级别的层序。

[0007] 现有技术的层序划分方法都是对测井资料进行分析，而在测井曲线上只能对短周期高频层序进行很好的识别，而长周期之间例如二级和三级层序，两者在测井曲线上的特征差异不大，采用现有技术难以通过测井曲线区分是二级还是三级层序。

发明内容

[0008] 为解决现有技术中无法准确对长周期低频层序进行准确识别的问题，本发明提出一种利用不同尺度识别不同级别层序的方法来进行层序划分，本发明实施例提供了一种地层层序划分方法，包括：

[0009] 根据地震剖面的地震反射特征识别二级层序界面、三级层序界面；

[0010] 将所述二级层序界面、三级层序界面分别定位到所述地震剖面中的一测井的对应深度；

[0011] 在所述二级层序界面的对应深度点处对所述测井的测井曲线进行一级重构，生成一级重构曲线；

[0012] 在所述三级层序界面的对应深度点处对所述一级重构曲线进行二级重构，生成二级重构曲线；

[0013] 按预设尺度对所述二级重构曲线进行小波分析，生成地层层序划分结果。

[0014] 本发明实施例中，将二级层序界面、三级层序界面分别定位到所述地震剖面中的一测井的对应深度包括：

[0015] 根据所述地震剖面中的测井的合成记录的时深关系分别将所述二级层序界面时间点t1，t2，…，tn对应到所述测井的测井曲线上的深度点h1，h2，…，hn；

[0016] 根据所述合成记录的时深关系将所述三级层序界面的时间点tn1，tn2，…，tnm对应到测井曲线上的深度点hn1，hn2，…，hnm；

[0017] 其中，tnm为第n个二级层序内第m个三级层序界面的时间点，n表示二级层序的次序，m表示三级层序界面的次序，hnm为第n个二级层序内第m个三级层序界面的深度点，n表示二级层序的次序，m表示三级层序界面的次序，m、n∈N+，N+为正整数。

[0018] 本发明实施例中，所述的在所述二级层序界面的对应深度点处对所述测井的测井曲线进行一级重构，生成的一级重构曲线为：

[0019] SB(h)＝SA(h)+SAmax/8+(SAmax/8)*(-1)n+1，hn

[0020] 其中，SB(h)为重构后的一级重构曲线B在各深度点处的测井值，SA(h)为测井曲线A在各深度点处的原测井值，SAmax为测井曲线的全局最大值，即原始测井曲线A整个深度段的最大测井值，hn为二级层序界面对应的不同深度点，n表示二级层序的次序，N+为正整数。

[0021] 本发明实施例中，所述的在所述三级层序界面的对应深度点处对所述一级重构曲线进行二级重构，生成的二级重构曲线为：

[0022] SC(h)＝SB(h)-SAmax(n)/16+(SAmax(n)/16)*(-1)m，hnm

[0024] 本发明实施例中，所述的按预设尺度对所述二级重构曲线进行小波分析，生成地层层序划分结果包括：

[0025] 按预设尺度对所述二级重构曲线进行小波分析，生成小波能量谱系图和小波系数曲线图；

[0026] 根据所述小波能量谱系图和小波系数曲线图的间断点对所述二级重构曲线进行层序划分。

[0027] 本发明实施例中，所述按预设尺度对所述二级重构曲线进行小波分析，生成地层层序划分结果进一步包括：

[0028] 对所述二级重构曲线进行6000尺度的小波分析，生成6000尺度的小波能量谱系图，根据该小波能量谱系图中的间断点划分二级层序；

[0029] 对所述二级重构曲线进行1400尺度的小波分析，生成1400尺度的小波能量谱系图，根据该小波能量谱系图中的间断点划分三级层序；

[0030] 对所述二级重构曲线进行700尺度的小波分析，生成700尺度的小波能量谱系图，根据该小波能量谱系图中的间断点划分四级层序；

[0031] 对所述二级重构曲线进行150尺度的小波分析，生成小波系数曲线图，根据小波系数曲线图中的间断点划分五级层序。

[0032] 同时，本发明还提供一种地层层序划分装置，包括：

[0033] 识别模块，用于根据地震剖面的地震反射特征识别二级层序界面、三级层序界面；

[0034] 定位模块，用于将所述二级层序界面、三级层序界面分别定位到所述地震剖面中的一测井的对应深度；

[0035] 一级重构模块，用于在所述二级层序界面的对应深度点处对所述测井的测井曲线进行一级重构，生成一级重构曲线；

[0036] 二级重构模块，在所述三级层序界面的对应深度点处对所述一级重构曲线进行二级重构，生成二级重构曲线；

[0037] 地层划分模块，用于按预设尺度对所述二级重构曲线进行小波分析，生成地层层序划分结果。

[0038] 本发明实施例中，定位模块包括：

[0039] 二级层序定位单元，根据所述地震剖面中的测井的合成记录的时深关系分别将所述二级层序界面时间点t1，t2，…，tn对应到所述测井的测井曲线上的深度点h1，h2，…，hn；

[0040] 三级层序定位单元，根据所述合成记录的时深关系将所述三级层序界面的时间点tn1，tn2，…，tnm对应到测井曲线上的深度点hn1，hn2，…，hnm；

[0041] 其中，tnm为第n个二级层序内第m个三级层序界面的时间点，n表示二级层序的次序，m表示三级层序界面的次序，hnm为第n个二级层序内第m个三级层序界面的深度点，n表示二级层序的次序，m表示三级层序界面的次序，m、n∈N+，N+为正整数。

[0042] 本发明实施例中，所述的地层划分模块包括：

[0043] 小波分析单元，用于按预设尺度对所述二级重构曲线进行小波分析，生成小波能量谱系图和小波系数曲线图；

[0044] 划分单元，用于根据所述小波能量谱系图和小波系数曲线图的间断点对所述二级重构曲线进行层序划分。

[0045] 本发明实施例中，小波分析单元进行小波分析包括：

[0046] 对所述二级重构曲线进行6000尺度的小波分析，生成6000尺度的小波能量谱系图，根据该小波能量谱系图中的间断点划分二级层序；

[0047] 对所述二级重构曲线进行1400尺度的小波分析，生成1400尺度的小波能量谱系图，根据该小波能量谱系图中的间断点划分三级层序；

[0048] 对所述二级重构曲线进行700尺度的小波分析，生成700尺度的小波能量谱系图，根据该小波能量谱系图中的间断点划分四级层序；

[0049] 对所述二级重构曲线进行150尺度的小波分析，生成小波系数曲线图，根据小波系数曲线图中的间断点划分五级层序。

[0050] 本发明对长期低频层序利用地震识别，再通过合成记录定位到测井上；这样通过地震对长期层序的准确识别保证了层序曲线中长期低频层序划分的结果准确可靠。并且，本发明依据地震对测井曲线经过两次重构得到二级重构曲线，最后直接对二级重构曲线进行不同尺度的小波变换分析，得到不同级别层序的划分方案，其实现了层序划分过程中真正的井震结合。

[0051] 为让本发明的上述和其他目的、特征和优点能更明显易懂，下文特举较佳实施例，并配合所附图式，作详细说明如下。

附图说明

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

[0053] 图1为本发明公开的一种地层层序划分方法的流程图；

[0054] 图2为本发明实施例中井的地震剖面图；

[0055] 图3为本发明实施例中二级重构曲线及小波变换能量谱系图；

[0056] 图4本发明实施例中二级重构曲线及小波变换能量谱系图；

[0057] 图5本发明实施例中二级重构曲线及小波变换能量谱系图；

[0058] 图6为二级重构曲线及其对应的150尺度小波系数曲线图；

[0059] 图7为本发明公开的一种地层层序划分装置的框图；

[0060] 图8为本发明实施例中的框图；

[0061] 图9为本发明实施例中的框图。

具体实施方式

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

[0063] 如图1所示，本发明提供一种地层层序划分方法，包括：

[0064] 步骤S101，根据地震剖面的地震反射特征识别二级层序界面、三级层序界面；

[0065] 步骤S102，将所述二级层序界面、三级层序界面分别定位到所述地震剖面中的一测井的对应深度；

[0066] 步骤S103，在所述二级层序界面的对应深度点处对所述测井的测井曲线进行一级重构，生成一级重构曲线；

[0067] 步骤S104，在所述三级层序界面的对应深度点处对所述一级重构曲线进行二级重构，生成二级重构曲线；

[0068] 步骤S105，按预设尺度对所述二级重构曲线进行小波分析，生成地层层序划分结果。

[0069] 本发明通过重构测井曲线，将地震剖面上识别的二、三级层序界面信息附加到测井上，使得其在后续小波变换中可以被准确识别，解决了原始测井曲线本身难以准确识别长期层序的难题。对重构后的测井曲线进行小波分析，通过选取不同尺度的不同小波进行小波变换，依据得到的小波时频能量图和小波系数曲线图来分析各个级别的层序旋回。

[0070] 由于1级层序尺度很大，很容易就在地震上识别；而6级层序尺度很小，是基于岩心进行识别，二者几乎都不需要以测井曲线为依据，故本发明不考虑这两种层序的识别，本发明主要针对现有技术基于测井曲线进行2-5级层序划分时，无法准确对长周期低频层序(二级、三级)进行准确识别的问题，提出一种利用不同尺度识别不同级别层序的方法来进行层序划分。借助地震对长周期旋回的准确识别能力，将其识别结果转换到测井曲线上，再通过对测井曲线进行分尺度的小波变换分析，逐级准确识别出不同级次的层序旋回。本发明具体实施时包括如下步骤：

[0071] 步骤1：依据地震剖面区域上明显的上超下削的地震反射特征，识别出n个二级层序界面，划分为n+1个二级层序；同时根据合成记录的时深关系将地震上二级层序界面的时间点t1，t2，…，tn对应到测井曲线上的深度点h1，h2，…，hn(其中n∈N+，N+为正整数)。

[0072] 步骤2：在上述各二级层序界面对应深度点h1，h2，…，hn处对原始测井曲线A进行一级重构，得到曲线B，成为一级重构曲线。曲线B的重构公式为：

[0073] SB(h)＝SA(h)+SAmax/8+(SAmax/8)*(-1)n+1，hn

[0075] 步骤3：依据地震剖面局部明显的上超下截的地震反射特征，识别出三级层序界面；同时根据合成记录的时深关系将地震上三级层序界面的时间点tn1，tn2，…，tnm(tnm为第n个二级层序内第m个三级层序界面的时间点，n表示二级层序的次序，m表示三级层序界面的次序)对应到测井曲线上的深度点hn1，hn2，…，hnm(hnm为第n个二级层序内第m个三级层序界面的深度点，n表示二级层序的次序，m表示三级层序界面的次序)(其中m、n∈N+，N+为正整数)。

[0076] 步骤4：在上述三级层序界面对应深度点hn1，hn2，…，hnm(其中m、n∈N+，N+为正整数)处对曲线B进行二级重构，得到曲线C，称为二级重构曲线。曲线C的重构公式为：

[0077] SC(h)＝SB(h)-SAmax(n)/16+(SAmax(n)/16)*(-1)m，hnm

[0079] 步骤5：对重构的曲线C进行6000尺度的小波分析，得到小波能量谱系图，将小波能量谱系图上的间断点作为二级层序的分界点，进而对重构曲线C中包含的二级层序进行划分。

[0080] 步骤6：对重构的曲线C进行1400尺度的小波分析，得到小波能量谱系图，将小波能量谱系图上的间断点作为三级层序分界点，进而对重构曲线C中包含的三级层序进行划分。

[0081] 步骤7：对重构的曲线C进行700尺度的小波分析，得到小波能量谱系图，将小波能量谱系图上的间断点作为四级层序分界点，进而对重构曲线C中包含的四级层序进行划分。

[0082] 步骤8：对重构的曲线C进行150尺度的小波分析，得到小波系数曲线图，将小波系数曲线图上的间断点作为五级层序的分界点，进而对重构曲线C中包含的五级层序进行划分。

[0083] 上述本发明实施例中，进行小波分析中尺度的选择为经验值，本发明的实现并不以此为限。

[0084] 本发明针对长期低频层序，依据地震剖面上的明显反射特征对其识别，通过合成记录将其转换到测井曲线上对应位置，再结合测井曲线本身良好的识别短周期高频层序的能力，对层序进行划分。具有以下两个显著优点：

[0085] 一、对长期低频层序利用地震识别，再通过合成记录定位到测井上；这样通过地震对长期层序的准确识别保证了层序曲线中长期低频层序划分的结果准确可靠。

[0086] 二、依据地震对测井曲线经过两次重构得到二级重构曲线，最后直接对二级重构曲线进行不同尺度的小波变换分析，得到不同级别层序的划分方案，其实现了层序划分过程中真正的井震结合。

[0087] 下面结合具体的实例对本发明技术方案做进一步详细说明。

[0088] 在辽河坳陷滩海地区选取海南洼陷新生界进行研究为例，海南洼陷新生界包括古近系和新近系，其中古近系发育东一段(d1)、东二段(d2)、东三段(d3)、沙一段(s1)、沙二段(s2)以及沙三段(s3)地层，沙四段(s4)地层不发育；新近系发育馆陶组(Ng)和明化镇组(Nm)地层。进行地层层序划分包括：

[0089] (1)如图2所示，在地震剖面上依据区域和局部的上超下削特征，可在古近系识别出1个二级层序界面(Ng组底)和3个三级层序界面(Ng_S亚段底、d3段底和s2段底)。

[0090] (2)选取地震剖面上的井well1作为实施例，Well1井的测量深度范围为190.00-3600.00m。对well1井制作合成记录，将上述地震上识别出的二、三级层序界面定位到井well1上的对应深度，其中1个二级层序界面Ng组底对应的深度为1167.00m，3个三级层序界面Ng_S亚段底、E3d3段底和E3s2段底对应的深度分别为721.00m、2230.50m和2971.00m，本发明实施例中图2是地震剖面，其左侧值并不是深度值，是时间值；Well1井的单位是深度单位。其中，a：明显削截，b1、b2：区域不整合，c1、c2:局部不整合，d：强反射，e:古近系或中生界底。

[0091] (3)选取井Well1的AC曲线(本实施例中该井的自然伽马测井GR、电阻率测井RT等曲线不完整)作为原始测井曲线A，找出其全局最大值为SAmax＝189us/m，依据公式(1)及二级层序界面深度点计算得到重构的一级重构曲线B。

[0092] (4)分别在以1167.00m为界的两个二级层序段内找出其局部最大值为SAmax(1)＝55us/m，SAmax(2)＝50us/m。依据公式(2)及三级层序界面深度点计算得到重构的二级重构曲线C，如图3中(a)所示，well1井AC曲线二级重构曲线。

[0093] (5)本发明实施例中使用Matlab软件中的小波工具箱对二级重构曲线C进行6000尺度的morl小波(在Matlab小波变换工具箱中，经过对db小波、sym小波、haar小波以及morl小波进行试验分析，发现选择morl小波进行小波变换的分析结果最好)变换，得到小波能量谱系图，如图3中(b)，曲线C6000尺度的小波变换能量谱系图。不同层序的沉积环境不同，反映在小波能量谱系图上即为不同层序形成的能量团之间存在明显间断点，且同一层序内部能量线由下至上收敛。因此，从图3中(b)中可以看出，在标注1处1167.00m处对应着Ng组底界的位置，可以很明显地将能量分为两部分，分别代表两个不同的二级层序，即将well1井新生界划分为两个二级层序，其与海南洼陷本身的层序地层特征相吻合。而针对原始曲线A同样进行6000尺度的morl小波变换，如图3中(c)所示，曲线A6000尺度的小波变换能量谱系图，可以看出，在原始曲线A的6000尺度的能量谱系图上，无法在将整体能量团在1167.00m分为两部分，即难以对二级层序进行准确划分。

[0094] (6)对二级重构曲线C进行1400尺度(经验值)的morl小波分析，得到小波能量谱系图，图4中(a)为well1井AC曲线二级重构曲线，图4中(b)为曲线C1400尺度小波变换能量谱系图，可知，3个三级层序的分界点可以清楚识别，即可以对三级层序进行准确划分；结合二级层序分界点，可将well1井新生界划分为5个三级层序，标注1为721.00m，标注2为2230.50m，标注3为2971.00m。

[0095] (7)在上述二、三级层序准确划分的基础上，开始进行较短周期层序的识别与划分。对二级重构曲线C进行700尺度(经验值)的morl小波分析得到小波能量谱系图，如图5所示。图5中(a)为well1井AC曲线二级重构曲线，图5中(b)为曲线C700尺度小波变换能量谱系图，根据图5中(b)的能量分布，可以识别出4个四级层序分界点；结合前述二、三级层序分界点，可以将well1井新生界划分为9个四级层序，标注1表示深度为487.25m，标注2表示1676.50m，标注3表示2683.75m，标注4表示3450.25m。

[0096] (8)对于短期层序的识别，由于其频率高，而能量谱系图的分辨率较低，因此采用小波系数曲线图进行分析。小波系数曲线的突变处可对应各级层序界面、湖泛面以及相转换面，反映沉积环境的变化。对二级重构曲线C进行150尺度(经验值)的morl小波分析，得到小波系数曲线图，如图所示，曲线a为well1井AC曲线二级重构曲线，b为曲线C150尺度小波系数曲线，根据图6中的小波曲线变化特征，寻找出10个突变点，结合前述二、三、四级层序分界点，可以将well1井新生界分为19个五级层序。

[0097] 综合以上8步的结果，可以准确地对海南洼陷well1井新生界进行层序划分。

[0098] 本发明通过在地震上依据区域不整合和局部不整合分别对二、三级层序进行识别，再选择剖面上所含的井，制作合成记录，将地震上识别的二三级层序分界点定位到测井曲线上对应的深度点。在对应深度点处通过经验公式对测井曲线进行重构，重构的目的是为了将二、三级层序界面特征放大，以此在准确划分长期层序的基础上控制内部短期层序划分的准确性。通过曲线重构得到可以直接用来进行小波变化划分层序的重构曲线，针对不同级别的层序划分，对重构曲线选用不同尺度的小波变换，得到对应的小波能量谱系图或小波系数曲线图，通过能谱图和系数曲线图中的间断点来识别各级层序分界点，最终可对研究区发育的层序旋回进行准确的划分。

[0099] 同时，本发明还公开一种地层层序划分装置，如图7所示，该装置包括：

[0100] 识别模块701，用于根据地震剖面的地震反射特征识别二级层序界面、三级层序界面；

[0101] 定位模块702，用于将所述二级层序界面、三级层序界面分别定位到所述地震剖面中的一测井的对应深度；

[0102] 一级重构模块703，用于在所述二级层序界面的对应深度点处对所述测井的测井曲线进行一级重构，生成一级重构曲线；

[0103] 二级重构模块704，在所述三级层序界面的对应深度点处对所述一级重构曲线进行二级重构，生成二级重构曲线；

[0104] 地层划分模块705，用于按预设尺度对所述二级重构曲线进行小波分析，生成地层层序划分结果。

[0105] 如图8所示，本发明实施例中的定位模块702包括：

[0106] 二级层序定位单元7021，根据所述地震剖面中的测井的合成记录的时深关系分别将所述二级层序界面时间点t1，t2，…，tn对应到所述测井的测井曲线上的深度点h1，h2，…，hn；

[0107] 三级层序定位单元7022，根据所述合成记录的时深关系将所述三级层序界面的时间点tn1，tn2，…，tnm对应到测井曲线上的深度点hn1，hn2，…，hnm；

[0108] 其中，tnm为第n个二级层序内第m个三级层序界面的时间点，n表示二级层序的次序，m表示三级层序界面的次序，hnm为第n个二级层序内第m个三级层序界面的深度点，n表示二级层序的次序，m表示三级层序界面的次序，m、n∈N+，N+为正整数。

[0109] 如图9所示，地层划分模块705包括：

[0110] 小波分析单元7051，用于按预设尺度对所述二级重构曲线进行小波分析，生成小波能量谱系图和小波系数曲线图；

[0111] 划分单元7052，用于根据所述小波能量谱系图和小波系数曲线图的间断点对所述二级重构曲线进行层序划分。

[0112] 本发明公开的装置解决问题的原理与上述方法的实现相似，不再赘述。

[0113] 本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

[0114] 本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

[0115] 这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

[0116] 这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

[0117] 本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。