基于钻井资料定量划分层序地层的方法及装置
阅读说明：本技术 基于钻井资料定量划分层序地层的方法及装置 (Method and device for quantitatively dividing sequence stratum based on drilling data ) 是由 郭泽清 张璐 刘卫红 郝翠果 于 2019-10-11 设计创作，主要内容包括：本发明提供了一种基于钻井资料定量划分层序地层的方法及装置。所述方法包括：对自然伽马测井曲线进行快速傅立叶变换,得到频率曲线图；确定频率曲线图上的主频,得到波长及各级旋回厚度；根据旋回厚度对自然伽马测井曲线进行数字滤波,得到各级旋回；计算钻井同一深度段内的砂地比,得到砂地比曲线；对自然伽马测井曲线进行连续小波变换,得到旋回响应特征曲线；利用自然伽马测井曲线、砂地比曲线及旋回响应特征曲线,确定单井的各级旋回及旋回界面；根据单井的各级旋回及旋回界面,建立联井对比剖面及层序对比格架。本发明具有设计合理,实用性强,制作成本低等优点,实现定量划分层序地层的级别和层序界面,提高划分结果的可靠性和准确性。(The invention provides a method and a device for quantitatively dividing sequence strata based on well drilling data. The method comprises the following steps: carrying out fast Fourier transform on the natural gamma logging curve to obtain a frequency curve graph; determining the main frequency on the frequency curve graph to obtain the wavelength and the thickness of each stage of convolution; performing digital filtering on the natural gamma well logging curve according to the convolution thickness to obtain each level of convolution; calculating the sand-to-ground ratio in the same depth section of the drilled well to obtain a sand-to-ground ratio curve; carrying out continuous wavelet transformation on the natural gamma well logging curve to obtain a convolution response characteristic curve; determining various levels of gyrus and gyrus interfaces of the single well by utilizing a natural gamma logging curve, a sand-to-ground ratio curve and a gyrus response characteristic curve; and establishing a well-connection comparison section and a sequence comparison framework according to each level of the gyrations and the gyrating interfaces of the single well. The method has the advantages of reasonable design, strong practicability, low manufacturing cost and the like, realizes quantitative classification of the level and the sequence interface of the sequence stratum, and improves the reliability and the accuracy of the classification result.)
技术领域
本发明涉及沉积学研究和油气勘探技术领域，尤指一种基于钻井资料定量划分层序地层的方法及装置。
背景技术
在对北美克拉通显生宙地层研究的基础上，Sloss等(1949)首次以“不整合限定的地层单元”引进了层序的概念。1977年《地震地层学》(AAPG26专辑)的出版正式宣告了层序地层学的诞生。Miall(1995)将层序地层学誉为沉积地质学继沉积—响应模型和将板块构造理论应用于沉积盆地分析之后的第三次革命。现在，层序地层学多个学派并存，造成模式的多样化。但是层序地层学的核心理念是在年代框架中研究地层的叠置样式和变化趋势，所以各学派对层序的定义都强调：1)旋回性，即层序代表岩石记录中的地层旋回，无论时间上是否对称；2)时间框架，即绘制同时期的岩相或沉积体系；3)成因上相联系的地层，即与选定的观察级别相比，体系域中没有重大沉积间断；4)可容纳空间和沉积作用的相互影响。
高分辨率层序地层学是对地层记录中反映基准面变化旋回的时间地层单元二元化分；该理论体系及其技术的关键是识别地层中多几次的基准面旋回和等时对比，通过对基准面旋回的不同层次性分析，实现不同级次的层序地层划分与对比，从而构建起高分辨率层序地层格架。
传统的层序地层界面识别主要依靠定性的相面法：根据野外露头上的不整合面，地震剖面上的削蚀、上超及顶超三种终止形式，以及钻井资料的岩性和电测曲线的突变面等界面来定性的识别层序界面和层序划分。依据这些方法识别层序界面和层序划分，因人而异，主观性太强，在操作上具有随意性，结果存在多解性。近年来，单独应用小波分析等定量化手段识别层序界面取得一定效果，但是缺乏其他定量化手段进行综合分析和判断，具有片面性，同时也未形成一套行之有效的技术思路和软件。
发明内容
为了解决上述问题，本发明实施例提供一种基于钻井资料定量划分层序地层的方法，所述方法包括：
对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；
确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度；
根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回；
计算钻井同一深度段内的砂地比，得到砂地比曲线；
对同一深度内的自然伽马测井曲线进行连续小波变换，得到多条不同频带范围内的旋回响应特征曲线；
利用数字滤波处理后的自然伽马测井曲线、砂地比曲线及小波变换后得到的旋回响应特征曲线，确定单井的各级旋回及旋回界面；
根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架。
可选的，在本发明一实施例中，所述对自然伽马测井曲线进行快速傅立叶变换之前还包括：对自然伽马测井曲线进行预处理，所述预处理包括环境校正及奇异点消除。
可选的，在本发明一实施例中，所述确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度包括：对由多个不同周期的沉积旋回叠加而成的自然伽马测井曲线进行数学变换，自然伽马测井曲线被分解成各自独立的周期旋回，以频率的形式记录成频谱曲线；其中，频谱曲线中频率值低，则沉积周期长，地层旋回厚度大；频谱曲线频率值高，则沉积周期短，地层旋回厚度小。
可选的，在本发明一实施例中，所述根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回包括：根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，将多个不同周期旋回信息的叠合的自然伽马测井曲线分解为不同周期的旋回曲线，其中，旋回曲线的最小值深度点对应层序界面，旋回曲线的最大值深度点对应最大湖泛面。
可选的，在本发明一实施例中，所述计算钻井同一深度段内的砂地比，得到砂地比曲线包括：根据钻井的同一深度段内的岩性剖面，以不同厚度单元为计算单元，计算各计算单元内累计砂岩厚度或地层厚度；将得到的累计砂岩厚度或地层厚度数据放到对应深度段中间深度处，以得到砂地比曲线。
本发明实施例还提供一种基于钻井资料定量划分层序地层的装置，所述装置包括：
频率曲线模块，用于对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；
旋回厚度模块，确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度；
各级旋回模块，根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回；
砂地比曲线模块，用于计算钻井同一深度段内的砂地比，得到砂地比曲线；
旋回响应特征曲线模块，用于对同一深度内的自然伽马测井曲线进行连续小波变换，得到多条不同频带范围内的旋回响应特征曲线；
单井旋回模块，用于利用数字滤波处理后的自然伽马测井曲线、砂地比曲线及小波变换后得到的旋回响应特征曲线，确定单井的各级旋回及旋回界面；
层序对比格架模块，用于根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架。
可选的，在本发明一实施例中，所述装置还包括：预处理模块，用于对自然伽马测井曲线进行预处理，所述预处理包括环境校正及奇异点消除。
可选的，在本发明一实施例中，所述旋回厚度模块包括：旋回厚度单元，用于对由多个不同周期的沉积旋回叠加而成的自然伽马测井曲线进行数学变换，自然伽马测井曲线被分解成各自独立的周期旋回，以频率的形式记录成频谱曲线；其中，频谱曲线中频率值低，则沉积周期长，地层旋回厚度大；频谱曲线频率值高，则沉积周期短，地层旋回厚度小。
可选的，在本发明一实施例中，所述各级旋回模块包括：各级旋回单元，用于根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，将多个不同周期旋回信息的叠合的自然伽马测井曲线分解为不同周期的旋回曲线，其中，旋回曲线的最小值深度点对应层序界面，旋回曲线的最大值深度点对应最大湖泛面。
可选的，在本发明一实施例中，所述砂地比曲线模块包括：厚度计算单元，用于根据钻井的同一深度段内的岩性剖面，以不同厚度单元为计算单元，计算各计算单元内累计砂岩厚度或地层厚度；砂地比曲线单元，用于将得到的累计砂岩厚度或地层厚度数据放到对应深度段中间深度处，以得到砂地比曲线。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；
确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度；
根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回；
计算钻井同一深度段内的砂地比，得到砂地比曲线；
对同一深度内的自然伽马测井曲线进行连续小波变换，得到多条不同频带范围内的旋回响应特征曲线；
利用数字滤波处理后的自然伽马测井曲线、砂地比曲线及小波变换后得到的旋回响应特征曲线，确定单井的各级旋回及旋回界面；
根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架。
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；
确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度；
根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回；
计算钻井同一深度段内的砂地比，得到砂地比曲线；
对同一深度内的自然伽马测井曲线进行连续小波变换，得到多条不同频带范围内的旋回响应特征曲线；
利用数字滤波处理后的自然伽马测井曲线、砂地比曲线及小波变换后得到的旋回响应特征曲线，确定单井的各级旋回及旋回界面；
根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架。
本发明基于油气田测井、录井资料定量划分陆源碎屑岩层序地层的方法，可以在编程和现有软件的基础上快速的、批量的划分单井的高分辨率层序地层，具有结构简单，设计合理，实用性强，制作成本低等优点，在油气勘探实践中，可以实现定量划分层序地层的级别和层序界面，提高划分结果的可靠性和准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例一种基于钻井资料定量划分层序地层的方法的流程图；
图2为本发明一具体实施例中基于钻井资料定量划分层序地层方法的流程图；
图3为本发明实施例一种基于钻井资料定量划分层序地层的装置的结构示意图。
具体实施方式
本发明实施例提供一种基于钻井资料定量划分层序地层的方法及装置。
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示为本发明实施例一种基于钻井资料定量划分层序地层的方法的流程图，图中所示方法包括：
步骤S1，对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；
步骤S2，确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度；
步骤S3，根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回；
步骤S4，计算钻井同一深度段内的砂地比，得到砂地比曲线；
步骤S5，对同一深度内的自然伽马测井曲线进行连续小波变换，得到多条不同频带范围内的旋回响应特征曲线；
步骤S6，利用数字滤波处理后的自然伽马测井曲线、砂地比曲线及小波变换后得到的旋回响应特征曲线，确定单井的各级旋回及旋回界面；
步骤S7，根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架。
在本实施例中，本发明建立工区测井、录井数据库，对自然伽马测井曲线预处理；对自然伽马曲线进行快速傅立叶变换(FFT)得到频率曲线图；求取频率曲线图上的主频，进而得到相应的波长，得出各级旋回厚度；对自然伽马曲线进行数字滤波处理，逐级过滤掉高频成分，在深度域得到各级旋回；计算同一深度段内的砂地比，得到砂地比曲线；选用某种小波对同一深度内的自然伽马曲线进行连续小波变换，获得多条不同频带范围内的旋回响应特征曲线；综合利用滤波后的测井曲线、砂地比曲线和小波变换后得到的旋回曲线，确定单井的各级旋回及旋回界面；建立联井对比剖面，建立层序对比格架。
作为本发明的一个实施例，对自然伽马测井曲线进行快速傅立叶变换之前还包括：对自然伽马测井曲线进行预处理，所述预处理包括环境校正及奇异点消除。
作为本发明的一个实施例，确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度包括：对由多个不同周期的沉积旋回叠加而成的自然伽马测井曲线进行数学变换，自然伽马测井曲线被分解成各自独立的周期旋回，以频率的形式记录成频谱曲线；其中，频谱曲线中频率值低，则沉积周期长，地层旋回厚度大；频谱曲线频率值高，则沉积周期短，地层旋回厚度小。
作为本发明的一个实施例，根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回包括：根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，将多个不同周期旋回信息的叠合的自然伽马测井曲线分解为不同周期的旋回曲线，其中，旋回曲线的最小值深度点对应层序界面，旋回曲线的最大值深度点对应最大湖泛面。
作为本发明的一个实施例，计算钻井同一深度段内的砂地比，得到砂地比曲线包括：根据钻井的同一深度段内的岩性剖面，以不同厚度单元为计算单元，计算各计算单元内累计砂岩厚度或地层厚度；将得到的累计砂岩厚度或地层厚度数据放到对应深度段中间深度处，以得到砂地比曲线。
在本发明一具体实施例中，如图2所示为基于钻井资料定量划分层序地层方法的流程图，本发明方法包括：1、建立工区测井、录井数据库，自然伽马测井曲线预处理；2、利用MATLAB语言编程，对某口井的某个深度段内的自然伽马曲线进行快速傅立叶变换(FFT)来进行频谱分析，得到频率曲线图，将自然伽马曲线从深度域转换为频率域；3、频率曲线图上功率较大(高点处)的频率对应于曲线的主频，进而求主频倒数既可以得到相应的波长，得出各级旋回厚度；4、根据已经确定的旋回厚度，选择适当的滤波参数，对自然伽马曲线进行数字滤波处理，逐级过滤掉高频成分，保留与上述计算的旋回厚度相对应的低频成分，就可以在深度域得到各级旋回；5、将上述同一口井的录井岩性剖面输入Bendlink软件，在执行宏命令窗口，利用CarbonScript脚本语言编制程序计算某一深度段内的砂地比，得到砂地比曲线；6、对上述同一口井的自然伽马曲线，选用某种小波对同一深度内进行连续小波变换，获得图形化的相似系数CWTa,b,和多条不同频带范围内的旋回响应特征曲线；7、综合利用滤波后的测井曲线、砂地比曲线和小波变换后得到的旋回曲线，确定单井的各级旋回及旋回界面；8、利用工区内所选的重点探井建立联井对比剖面，建立层序对比格架。本发明在油气勘探实践中，可以实现定量划分层序地层的级别和层序界面，提高划分结果的可靠性和准确性。方法具体过程包括：
(1)建立工区测井、录井数据库，自然伽马测井曲线预处理。
(2)利用MATLAB语言编程对某口井的某个深度段内的自然伽马曲线进行快速傅立叶变换(FFT)来进行频谱分析，得到频率曲线图，将自然伽马曲线从深度域转换为频率域。
(3)频率曲线图上功率较大(高点处)的频率对应于曲线的主频，进而求主频倒数既可以得到相应的波长，得出各级旋回厚度。
(4)根据已经确定的旋回厚度，选择适当的滤波参数，对自然伽马曲线进行数字滤波处理，逐级过滤掉高频成分，保留与上述计算的旋回厚度相对应的低频成分，就可以在深度域得到各级旋回。
(5)将上述同一口井的录井岩性剖面输入Bendlink软件，在执行宏命令窗口，利用CarbonScript脚本语言编制程序计算某一深度段内的砂地比，得到砂地比曲线。
(6)对上述同一口井的自然伽马曲线，选用某种小波对同一深度内进行连续小波变换，获得图形化的相似系数CWTa,b,和多条不同频带范围内的旋回响应特征曲线。
(7)综合利用滤波后的测井曲线、砂地比曲线和小波变换后得到的旋回曲线，确定单井的各级旋回及旋回界面。
(8)利用工区内所选的重点探井建立联井对比剖面，建立层序对比格架。
其中，在步骤(1)中，建立工区数据库以后，首先对测井曲线进行预处理，包括环境校正和奇异点消除。由于在测井时受到的各种影响因素较多，往往包含一些非地质因素，如井壁滑塌、泥浆入侵等，因此在使用之前要经过环境校正。由于原始测井曲线记录的质量不同，来源不同，可能存在一些强烈干扰真实情况的奇异点，如果不消除这些点，就会干扰对真实情况的判断。
其中，在步骤(2)中，频谱分析的方法较多，其中最常用的是快速傅立叶变换法。傅里叶变换函数如果以连续函数得出，变换函数为x(t)，则傅里叶变换由下式得出：
式(1)中t为时间，f为频率，i为虚数单位。
如果变换函数为N项离散时间序列xm，则其离散傅里叶变换为：
式(2)中Xk为频谱值(k＝0,1,2,…,N-1)。
得到的Xk为复数，没有直观意义，习惯上总是把其转换为能量，构成直观的能量－频率图，对应于频率fk的谐振能量为：
Pk＝(real(Xk))2+(imag(Xk))2 (k＝0,1,2,…,N-1) (3)
其中，在步骤(3)中，由多个不同周期的沉积旋回叠加而成的自然伽马测井曲线，通过数学变换，被分解成各自独立的周期旋回，以频率的形式记录成频谱曲线。频率值低，表明该沉积周期长，表现为地层旋回厚度大；反之频率值高，表明该沉积周期短，地层旋回厚度就小。
其中，在步骤(4)中，在频谱域分析得到米兰柯维奇旋回厚度之后，应该返回到深度域，指出对应于不同旋回在深度剖面上的位置，这样才能用于区域对比。滤波器的时间响应为零相位，保持原始信号的深度特征不变。以上述井的同一深度段的自然伽马测井曲线进行滤波处理，将多个不同周期旋回信息的叠合的测井曲线分解为不同周期的旋回曲线，曲线的最小值深度点对应层序界面，最大值深度点对应最大湖泛面。
其中，在步骤(5)中，以上述井的同一深度段内的岩性剖面，分别以不同厚度单元(例如以10m，20m，30m分别代表短期、中期、长期旋回)为计算单元，计算该单元内累计砂岩厚度/地层厚度，将得到的数据放到该深度段中间深度处，把这些点连接起来即得到砂地比曲线。曲线的最大值深度点对应层序界面，最小值深度点对应最大湖泛面。
其中，在步骤(6)中，在连续小波变换中用得最多的是墨西哥帽小波和Morlet小波。Morlet小波作为小波母函数，它是一个周期函数，并且是带Guass包络调制的、频率为ω0的谐波。测井曲线选用Morlet小波进行连续小波变换，其尺度因子a与周期T之间有如下关系：尺度a越小表示小波被压缩，空间范围上观察窗口小，对应于测井信号的短周期分量。尺度a越大表示小波被伸展，空间范围上观察窗口大，对应于测井信号的长周期分量。测井信号经小波变换处理得到清晰的频率结构，并可探测到各频率段之间的突变区，形成特定的界面响应。通过考察多尺度下明显周期性振荡特征与各级层序界面建立一定的对应关系，作为层序划分的参考。
本发明基于油气田测井、录井资料定量划分陆源碎屑岩层序地层的方法，可以在编程和现有软件的基础上快速的、批量的划分单井的高分辨率层序地层，具有结构简单，设计合理，实用性强，制作成本低等优点，在油气勘探实践中，可以实现定量划分层序地层的级别和层序界面，提高划分结果的可靠性和准确性。
如图3所示为本发明实施例一种基于钻井资料定量划分层序地层的装置的结构示意图，图中所示装置包括：
频率曲线模块10，用于对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；
旋回厚度模块20，确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度；
各级旋回模块30，根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回；
砂地比曲线模块40，用于计算钻井同一深度段内的砂地比，得到砂地比曲线；
旋回响应特征曲线模块50，用于对同一深度内的自然伽马测井曲线进行连续小波变换，得到多条不同频带范围内的旋回响应特征曲线；
单井旋回模块60，用于利用数字滤波处理后的自然伽马测井曲线、砂地比曲线及小波变换后得到的旋回响应特征曲线，确定单井的各级旋回及旋回界面；
层序对比格架模块70，用于根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架。
作为本发明的一个实施例，装置还包括：预处理模块，用于对自然伽马测井曲线进行预处理，所述预处理包括环境校正及奇异点消除。
作为本发明的一个实施例，旋回厚度模块包括：旋回厚度单元，用于对由多个不同周期的沉积旋回叠加而成的自然伽马测井曲线进行数学变换，自然伽马测井曲线被分解成各自独立的周期旋回，以频率的形式记录成频谱曲线；其中，频谱曲线中频率值低，则沉积周期长，地层旋回厚度大；频谱曲线频率值高，则沉积周期短，地层旋回厚度小。
作为本发明的一个实施例，各级旋回模块包括：各级旋回单元，用于根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，将多个不同周期旋回信息的叠合的自然伽马测井曲线分解为不同周期的旋回曲线，其中，旋回曲线的最小值深度点对应层序界面，旋回曲线的最大值深度点对应最大湖泛面。
作为本发明的一个实施例，砂地比曲线模块包括：厚度计算单元，用于根据钻井的同一深度段内的岩性剖面，以不同厚度单元为计算单元，计算各计算单元内累计砂岩厚度或地层厚度；砂地比曲线单元，用于将得到的累计砂岩厚度或地层厚度数据放到对应深度段中间深度处，以得到砂地比曲线。
基于与上述一种基于钻井资料定量划分层序地层的方法相同的申请构思，本发明还提供了上述一种基于钻井资料定量划分层序地层的装置。由于该一种基于钻井资料定量划分层序地层的装置解决问题的原理与一种基于钻井资料定量划分层序地层的方法相似，因此该一种基于钻井资料定量划分层序地层的装置的实施可以参见一种基于钻井资料定量划分层序地层的方法的实施，重复之处不再赘述。
本发明基于油气田测井、录井资料定量划分陆源碎屑岩层序地层的方法，可以在编程和现有软件的基础上快速的、批量的划分单井的高分辨率层序地层，具有结构简单，设计合理，实用性强，制作成本低等优点，在油气勘探实践中，可以实现定量划分层序地层的级别和层序界面，提高划分结果的可靠性和准确性。
本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现以下步骤：
对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；
确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度；
根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回；
计算钻井同一深度段内的砂地比，得到砂地比曲线；
对同一深度内的自然伽马测井曲线进行连续小波变换，得到多条不同频带范围内的旋回响应特征曲线；
利用数字滤波处理后的自然伽马测井曲线、砂地比曲线及小波变换后得到的旋回响应特征曲线，确定单井的各级旋回及旋回界面；
根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架。
本发明实施例还提供一种计算机可读存储介质，其上存储有计算机程序，所述计算机程序被处理器执行时实现以下步骤：
对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；
确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度；
根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回；
计算钻井同一深度段内的砂地比，得到砂地比曲线；
对同一深度内的自然伽马测井曲线进行连续小波变换，得到多条不同频带范围内的旋回响应特征曲线；
利用数字滤波处理后的自然伽马测井曲线、砂地比曲线及小波变换后得到的旋回响应特征曲线，确定单井的各级旋回及旋回界面；
根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架。
基于与上述一种基于钻井资料定量划分层序地层的方法相同的申请构思，本发明还提供了上述一种计算机设备及一种计算机可读存储介质。由于该一种计算机设备及一种计算机可读存储介质解决问题的原理与一种基于钻井资料定量划分层序地层的方法相似，因此该一种计算机设备及一种计算机可读存储介质的实施可以参见一种基于钻井资料定量划分层序地层的方法的实施，重复之处不再赘述。
本发明基于油气田测井、录井资料定量划分陆源碎屑岩层序地层的方法，可以在编程和现有软件的基础上快速的、批量的划分单井的高分辨率层序地层，具有结构简单，设计合理，实用性强，制作成本低等优点，在油气勘探实践中，可以实现定量划分层序地层的级别和层序界面，提高划分结果的可靠性和准确性。
本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
- 上一篇：一种医用注射器针头装配设备
- 下一篇：一种基于能谱测井特征谱峰的铀矿定量刻度系数求法