沉积旋回自动划分方法、装置、电子设备及介质与流程

1.本发明属于地球物理勘探技术领域，具体涉及一种沉积旋回自动划分方法、装置、电子设备及存储介质。


背景技术：

2.沉积旋回是指按颗粒从大到小、比重从大到小的顺序先后分层沉积而成岩层的规律，在地层剖面上，表现为在从老到新的顺序，岩层依粒度从粗到细，依次为粗砂岩—中砂岩—细砂岩—粉砂岩—泥岩。沉积旋回注重岩性沉积旋回，沉积旋回是指沉积作用和沉积条件按相同的次序不断重复沉积而组成的一个层序。沉积旋回是在地层垂直方向上相似岩性与岩相，例如岩石的颜色/粒度等的规律性周期重复。储层中的沉积旋回一般规模较小。
3.在进行油气勘探的过程中，在对目标底层进行沉积旋回划分时，通常需要先对沉积旋回进行分析。常规的沉积旋回划分方法是依据层序地层学或高分辨率层序地层学为指导，以研究区内钻井资料和地震资料为基础对沉积旋回界面和性质进行判定，但该方法工作量巨大且人为干扰过多，存在多解性。
4.此外，由于沉积物粒度的变化会导致地震剖面的频谱能量发生变化，也可以通过小波变换对地震信号的频谱能量进行分析，也可以对沉积旋回进行分析判定，但地震资料往往缺失高频成分，分辨率较低，所以对小规模的储层沉积旋回进行判定仍存在一定困难。
5.因此，本领域需要提供一种新的沉积旋回自动划分方法，减低工作量并减少认为干扰，提高划分精度。


技术实现要素：

6.针对现有方法划分沉积旋回存在工作量巨大，人为干扰多，存在多解性及分辨率低等问题，本文提出一种沉积旋回自动划分方法，该方法主要是以高分辨率层序地层学为指导，通过优选测井曲线，对优选曲线进行频谱分析并提取差值曲线实现开发阶段密井网条件下多井旋回自动化追踪，批量生成各口井的分层结果，提高工作效率及正确性，减少人工干预和调整。
7.根据本发明的一个方面，提供一种沉积旋回自动划分方法，包括：
8.根据储层类型选取对岩性敏感的测井曲线；
9.对所述测井曲线提取差值曲线；
10.根据所述差值曲线的变化趋势划分沉积旋回。
11.进一步地，该方法还包括对选取的所述测井曲线进行预处理，包括剔除奇异值、校正、滤波和标准化。
12.进一步地，通过对测井曲线值计算自相关函数，获得差值曲线。
13.进一步地，测井曲线值的自相关函数为：
14.rn(r1,r2,r3,...,rn)
15.对该函数进行外推得到函数的第n+1个值(r
n+1
)，
[0016][0017]
通过计算自相关函数的每个值与实测数据之间的差值获得差值曲线xn。
[0018]
进一步地，通过实测测井数据值与计算所得自相关函数值可得每一个深度点测井数据的差值，根据所得的差值曲线xn的走势，对沉积旋回进行划分。
[0019]
进一步低，该方法还包括对所得的沉积旋回进行成图。
[0020]
根据本发明的另一方面，提供一种沉积旋回自动划分装置，包括：
[0021]
选择模块，根据储层类型选取对岩性敏感的测井曲线；
[0022]
差值模块，对所述测井曲线提取差值曲线；
[0023]
划分模块，根据所述差值曲线的变化趋势划分沉积旋回。
[0024]
具体地，所述差值模块通过对测井曲线值计算自相关函数，通过计算自相关函数的每个值与实测数据之间的差值，获得差值曲线；
[0025]
测井曲线值的自相关函数为：
[0026]rn
(r1,r2,r3,...,rn)
[0027]
对该函数进行外推得到函数的第n+1个值(r
n+1
)，
[0028][0029]
通过计算自相关函数的每个值与实测数据之间的差值获得差值曲线xn。
[0030]
根据本发明的另一方面，提供一种电子设备，所述电子设备包括：
[0031]
存储器，存储有可执行指令；
[0032]
处理器，所述处理器运行所述存储器中的所述可执行指令，以实现所述的沉积旋回自动划分方法。
[0033]
根据本发明的另一方面，提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的沉积旋回自动划分方法。
[0034]
本发明对比已有的技术具有以下创新点：
[0035]
通过自相关函数的计算，提高频谱分辨率；通过提取差值曲线，自动划分沉积旋回，大大较少工作量并减少人为干扰，使结果更加准确。
附图说明
[0036]
通过结合附图对本公开示例性实施方式进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施方式中，相同的参考标号通常代表相同部件。
[0037]
图1为根据本发明的沉积旋回自动划分方法流程图。
[0038]
图2为根据本发明实施例的原始测井曲线与计算所得差值曲线。
[0039]
图3为根据本发明实施例的人工划分沉积旋回。
[0040]
图4为根据本发明实施例的自动划分沉积旋回。
[0041]
图5为根据本发明实施例的技术流程图。
具体实施方式
[0042]
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施方
式所限制。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
[0043]
常规的旋回分析以研究区内钻井资料和地震资料为基础对沉积旋回界面和性质进行判定，但该方法工作量巨大且人为干扰过多，存在多解性。此外，通过小波变换对地震信号的频谱能量进行分析，也可以对沉积旋回进行分析判定，但地震资料往往缺失高频成分，分辨率较低。
[0044]
本发明提出一种沉积旋回自动划分方法，该方法主要是以高分辨率层序地层学为指导，通过优选测井曲线，对优选曲线进行频谱分析并提取差值曲线实现开发阶段密井网条件下多井旋回自动化追踪，批量生成各口井的分层结果，提高工作效率及正确性，减少人工干预和调整。
[0045]
如图1所示，本发明提出了一种沉积旋回自动划分方法，包括：
[0046]
根据储层类型选取对岩性敏感的测井曲线；
[0047]
对所述测井曲线提取差值曲线；
[0048]
根据所述差值曲线的变化趋势划分沉积旋回。
[0049]
具体地，本发明创建了一种依据高分辨率层序地层学，对测井曲线进行频谱分析，并对曲线提取差值曲线进行沉积旋回的自动划分，详细内容如下：
[0050]
首先，选取对岩性敏感的测井曲线。
[0051]
优选地，还可以对原始测井曲线进行预处理，主要包括剔除奇异值，校正，滤波和标准化等。
[0052]
接下来，对测井曲线值进行频谱估计。对原始测井曲线值计算自相关函数rn(r1,r2,r3,...,rn)，对该函数进行外推得到函数的第n+1个值(r
n+1
)，其符合条件为：
[0053][0054]
其中
[0055][0056][0057]
令可得求取的梯度矩阵d(r)
[0058][0059]
求取梯度矩阵的解，即
[0060][0061]
其中为梯度矩阵的最后一行，即
[0062][0063][0064]
这样，可以求得r
n+1
，依次求解，通过计算自相关函数的每个值与实测数据之间的差值即可获得差值曲线xn。
[0065]
(3)根据所得差值曲线xn的走势，上粗下细为正旋回，上细下粗为反旋回，对沉积旋回进行划分。
[0066]
(4)对所得沉积旋回进行成图。
[0067]
为便于理解本发明实施例的方案及其效果，以下给出具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
[0068]
实施例1
[0069]
参照图2-图5，描述本发明的一个具体应用实施例。
[0070]
首先，根据储层类型选取敏感曲线；对选取的敏感曲线进行曲线预处理；接下来，对敏感曲线提相应的差值曲线；最后根据差值曲线的变化趋势划分沉积旋回。
[0071]
接下来，本实施例通过对比人工划分沉积旋回和根据差值曲线自动划分沉积旋回的结果对比说明：
[0072]
从川孝603井旋回特征及砂组发育特征看(参照图3)，沙溪庙组沉积期，js2时期沉积水体最浅，砂岩最为发育的，期次为砂组较为发育的下沙溪庙组时期，js1沉积时期水体最深，沉积能量最弱，主要以泥质沉积为主。蓬莱镇组沉积期，jp1、jp2时期沉积期，工区北东和北西物源区处分别处于三角洲前缘，发育水下分流河道微相为主的砂岩，砂岩最为发育的，jp3沉积时期水体最深，主要发育前三角洲—滨浅湖，沉积能量最弱，主要以泥质沉积为主。
[0073]
根据cx603井的gr曲线对中期、长期、短期沉积旋回进行人工划分，具体结果如图3所示，同时对gr曲线进行频谱分析，提取gr曲线的差值曲线，根据差值曲线的趋势形态和转折点(参照图2)进行自动追踪，从而确定旋回界面，实现沉积旋回的自动划分(参照图4)。
[0074]
对比蓬三段，蓬二段和蓬一段的结果表面，中长期旋回自动划分结果与人工划分
结果吻合度更高，短期旋回更精细。
[0075]
实施例2
[0076]
本实施例提供一种沉积旋回自动划分装置，包括：
[0077]
选择模块，根据储层类型选取对岩性敏感的测井曲线；
[0078]
差值模块，对所述测井曲线提取差值曲线；
[0079]
划分模块，根据所述差值曲线的变化趋势划分沉积旋回。
[0080]
选择模块、差值模块、划分模块依次通信连接。选择模块将选取的敏感测井曲线发送至差值模块，差值模块对敏感测井曲线提取差值曲线，并将提取结果发送至划分模块，划分模块根据所述差值曲线的变化趋势划分沉积旋回。
[0081]
实施例3
[0082]
本实施例提供一种电子设备包括，该电子设备包括：存储器，存储有可执行指令；处理器，处理器运行存储器中的可执行指令，以实现上述的沉积旋回自动划分方法。
[0083]
根据本公开实施例的电子设备包括存储器和处理器。
[0084]
该存储器用于存储非暂时性计算机可读指令。具体地，存储器可以包括一个或多个计算机程序产品，该计算机程序产品可以包括各种形式的计算机可读存储介质，例如易失性存储器和/或非易失性存储器。该易失性存储器例如可以包括随机存取存储器(ram)和/或高速缓冲存储器(cache)等。该非易失性存储器例如可以包括只读存储器(rom)、硬盘、闪存等。
[0085]
该处理器可以是中央处理单元(cpu)或者具有数据处理能力和/或指令执行能力的其它形式的处理单元，并且可以控制电子设备中的其它组件以执行期望的功能。在本公开的一个实施例中，该处理器用于运行该存储器中存储的该计算机可读指令。
[0086]
本领域技术人员应能理解，为了解决如何获得良好用户体验效果的技术问题，本实施例中也可以包括诸如通信总线、接口等公知的结构，这些公知的结构也应包含在本公开的保护范围之内。
[0087]
有关本实施例的详细说明可以参考前述各实施例中的相应说明，在此不再赘述。
[0088]
实施例4
[0089]
本实施例提供一种计算机可读存储介质，该计算机可读存储介质存储有计算机程序，该计算机程序被处理器执行时实现所述的沉积旋回自动划分方法。
[0090]
根据本公开实施例的计算机可读存储介质，其上存储有非暂时性计算机可读指令。当该非暂时性计算机可读指令由处理器运行时，执行前述的本公开各实施例方法的全部或部分步骤。
[0091]
上述计算机可读存储介质包括但不限于：光存储介质(例如：cd－rom和dvd)、磁光存储介质(例如：mo)、磁存储介质(例如：磁带或移动硬盘)、具有内置的可重写非易失性存储器的媒体(例如：存储卡)和具有内置rom的媒体(例如：rom盒)。
[0092]
本领域技术人员应理解，上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果，并不意在将本发明的实施例限制于所给出的任何示例。
[0093]
以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。