CN112649887B - 基于钻井资料定量划分层序地层的方法及装置 - Google Patents
基于钻井资料定量划分层序地层的方法及装置 Download PDFInfo
- Publication number
- CN112649887B CN112649887B CN201910962403.7A CN201910962403A CN112649887B CN 112649887 B CN112649887 B CN 112649887B CN 201910962403 A CN201910962403 A CN 201910962403A CN 112649887 B CN112649887 B CN 112649887B
- Authority
- CN
- China
- Prior art keywords
- gyratory
- curve
- thickness
- natural gamma
- frequency
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000005553 drilling Methods 0.000 title claims abstract description 48
- 238000000034 method Methods 0.000 title claims abstract description 43
- 238000001914 filtration Methods 0.000 claims abstract description 50
- 230000009466 transformation Effects 0.000 claims abstract description 42
- 230000004044 response Effects 0.000 claims abstract description 35
- 230000008021 deposition Effects 0.000 claims description 30
- 238000012545 processing Methods 0.000 claims description 28
- 238000001228 spectrum Methods 0.000 claims description 25
- 230000000737 periodic effect Effects 0.000 claims description 18
- 238000004590 computer program Methods 0.000 claims description 17
- 230000015572 biosynthetic process Effects 0.000 claims description 13
- 238000004364 calculation method Methods 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 8
- 230000007613 environmental effect Effects 0.000 claims description 8
- 230000008030 elimination Effects 0.000 claims description 7
- 238000003379 elimination reaction Methods 0.000 claims description 7
- 238000000638 solvent extraction Methods 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 5
- 238000004519 manufacturing process Methods 0.000 abstract description 6
- 238000013461 design Methods 0.000 abstract description 5
- 238000005755 formation reaction Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 238000004458 analytical method Methods 0.000 description 5
- 239000011435 rock Substances 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 150000002500 ions Chemical class 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000003595 spectral effect Effects 0.000 description 2
- 238000010183 spectrum analysis Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000003628 erosive effect Effects 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V5/00—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity
- G01V5/04—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging
- G01V5/06—Prospecting or detecting by the use of ionising radiation, e.g. of natural or induced radioactivity specially adapted for well-logging for detecting naturally radioactive minerals
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- High Energy & Nuclear Physics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种基于钻井资料定量划分层序地层的方法及装置。所述方法包括：对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；确定频率曲线图上的主频，得到波长及各级旋回厚度；根据旋回厚度对自然伽马测井曲线进行数字滤波，得到各级旋回；计算钻井同一深度段内的砂地比，得到砂地比曲线；对自然伽马测井曲线进行连续小波变换，得到旋回响应特征曲线；利用自然伽马测井曲线、砂地比曲线及旋回响应特征曲线，确定单井的各级旋回及旋回界面；根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架。本发明具有设计合理，实用性强，制作成本低等优点，实现定量划分层序地层的级别和层序界面，提高划分结果的可靠性和准确性。
Description
技术领域
本发明涉及沉积学研究和油气勘探技术领域，尤指一种基于钻井资料定量划分层序地层的方法及装置。
背景技术
在对北美克拉通显生宙地层研究的基础上，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、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
Claims (6)
1.一种基于钻井资料定量划分层序地层的方法，其特征在于，所述方法包括：
对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；
确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度；
根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回；
计算钻井同一深度段内的砂地比，得到砂地比曲线；
对同一深度内的自然伽马测井曲线进行连续小波变换，得到多条不同频带范围内的旋回响应特征曲线；
利用数字滤波处理后的自然伽马测井曲线、砂地比曲线及小波变换后得到的旋回响应特征曲线，确定单井的各级旋回及旋回界面；
根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架；
其中，所述计算钻井同一深度段内的砂地比，得到砂地比曲线包括：
根据钻井的同一深度段内的岩性剖面，以不同厚度单元为计算单元，计算各计算单元内累计砂岩厚度或地层厚度；
将得到的累计砂岩厚度或地层厚度数据放到对应深度段中间深度处，以得到砂地比曲线；
其中，所述确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度包括：对由多个不同周期的沉积旋回叠加而成的自然伽马测井曲线进行数学变换，自然伽马测井曲线被分解成各自独立的周期旋回，以频率的形式记录成频谱曲线；其中，频谱曲线中频率值低，则沉积周期长，地层旋回厚度大；频谱曲线频率值高，则沉积周期短，地层旋回厚度小；
其中，所述根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回包括：根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，将多个不同周期旋回信息的叠合的自然伽马测井曲线分解为不同周期的旋回曲线，其中，旋回曲线的最小值深度点对应层序界面，旋回曲线的最大值深度点对应最大湖泛面。
2.根据权利要求1所述的方法，其特征在于，所述对自然伽马测井曲线进行快速傅立叶变换之前还包括：对自然伽马测井曲线进行预处理，所述预处理包括环境校正及奇异点消除。
3.一种基于钻井资料定量划分层序地层的装置，其特征在于，所述装置包括：
频率曲线模块，用于对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；
旋回厚度模块，确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度；
各级旋回模块，根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回；
砂地比曲线模块，用于计算钻井同一深度段内的砂地比，得到砂地比曲线；
旋回响应特征曲线模块，用于对同一深度内的自然伽马测井曲线进行连续小波变换，得到多条不同频带范围内的旋回响应特征曲线；
单井旋回模块，用于利用数字滤波处理后的自然伽马测井曲线、砂地比曲线及小波变换后得到的旋回响应特征曲线，确定单井的各级旋回及旋回界面；
层序对比格架模块，用于根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架；
其中，所述砂地比曲线模块包括：
厚度计算单元，用于根据钻井的同一深度段内的岩性剖面，以不同厚度单元为计算单元，计算各计算单元内累计砂岩厚度或地层厚度；
砂地比曲线单元，用于将得到的累计砂岩厚度或地层厚度数据放到对应深度段中间深度处，以得到砂地比曲线；
其中，所述旋回厚度模块包括：旋回厚度单元，用于对由多个不同周期的沉积旋回叠加而成的自然伽马测井曲线进行数学变换，自然伽马测井曲线被分解成各自独立的周期旋回，以频率的形式记录成频谱曲线；其中，频谱曲线中频率值低，则沉积周期长，地层旋回厚度大；频谱曲线频率值高，则沉积周期短，地层旋回厚度小；
其中，所述各级旋回模块包括：各级旋回单元，用于根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，将多个不同周期旋回信息的叠合的自然伽马测井曲线分解为不同周期的旋回曲线，其中，旋回曲线的最小值深度点对应层序界面，旋回曲线的最大值深度点对应最大湖泛面。
4.根据权利要求3所述的装置，其特征在于，所述装置还包括：预处理模块，用于对自然伽马测井曲线进行预处理，所述预处理包括环境校正及奇异点消除。
5.一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述计算机程序时实现以下步骤：
对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；
确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度；
根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回；
计算钻井同一深度段内的砂地比，得到砂地比曲线；
对同一深度内的自然伽马测井曲线进行连续小波变换，得到多条不同频带范围内的旋回响应特征曲线；
利用数字滤波处理后的自然伽马测井曲线、砂地比曲线及小波变换后得到的旋回响应特征曲线，确定单井的各级旋回及旋回界面；
根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架；
其中，所述计算钻井同一深度段内的砂地比，得到砂地比曲线包括：
根据钻井的同一深度段内的岩性剖面，以不同厚度单元为计算单元，计算各计算单元内累计砂岩厚度或地层厚度；
将得到的累计砂岩厚度或地层厚度数据放到对应深度段中间深度处，以得到砂地比曲线；
其中，所述确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度包括：对由多个不同周期的沉积旋回叠加而成的自然伽马测井曲线进行数学变换，自然伽马测井曲线被分解成各自独立的周期旋回，以频率的形式记录成频谱曲线；其中，频谱曲线中频率值低，则沉积周期长，地层旋回厚度大；频谱曲线频率值高，则沉积周期短，地层旋回厚度小；
其中，所述根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回包括：根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，将多个不同周期旋回信息的叠合的自然伽马测井曲线分解为不同周期的旋回曲线，其中，旋回曲线的最小值深度点对应层序界面，旋回曲线的最大值深度点对应最大湖泛面。
6.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现以下步骤：
对自然伽马测井曲线进行快速傅立叶变换，得到频率曲线图；
确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度；
根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回；
计算钻井同一深度段内的砂地比，得到砂地比曲线；
对同一深度内的自然伽马测井曲线进行连续小波变换，得到多条不同频带范围内的旋回响应特征曲线；
利用数字滤波处理后的自然伽马测井曲线、砂地比曲线及小波变换后得到的旋回响应特征曲线，确定单井的各级旋回及旋回界面；
根据单井的各级旋回及旋回界面，建立联井对比剖面及层序对比格架；
其中，所述计算钻井同一深度段内的砂地比，得到砂地比曲线包括：
根据钻井的同一深度段内的岩性剖面，以不同厚度单元为计算单元，计算各计算单元内累计砂岩厚度或地层厚度；
将得到的累计砂岩厚度或地层厚度数据放到对应深度段中间深度处，以得到砂地比曲线；
其中，所述确定所述频率曲线图上的主频，得到相应的波长及各级旋回厚度包括：对由多个不同周期的沉积旋回叠加而成的自然伽马测井曲线进行数学变换，自然伽马测井曲线被分解成各自独立的周期旋回，以频率的形式记录成频谱曲线；其中，频谱曲线中频率值低，则沉积周期长，地层旋回厚度大；频谱曲线频率值高，则沉积周期短，地层旋回厚度小；
其中，所述根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，在深度域得到各级旋回包括：根据所述旋回厚度对所述自然伽马测井曲线进行数字滤波处理，将多个不同周期旋回信息的叠合的自然伽马测井曲线分解为不同周期的旋回曲线，其中，旋回曲线的最小值深度点对应层序界面，旋回曲线的最大值深度点对应最大湖泛面。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910962403.7A CN112649887B (zh) | 2019-10-11 | 2019-10-11 | 基于钻井资料定量划分层序地层的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910962403.7A CN112649887B (zh) | 2019-10-11 | 2019-10-11 | 基于钻井资料定量划分层序地层的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112649887A CN112649887A (zh) | 2021-04-13 |
| CN112649887B true CN112649887B (zh) | 2024-03-01 |
Family
ID=75342684
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910962403.7A Active CN112649887B (zh) | 2019-10-11 | 2019-10-11 | 基于钻井资料定量划分层序地层的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112649887B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114002750B (zh) * | 2021-09-27 | 2023-03-24 | 中国石油大学(北京) | 页岩层序识别方法、装置、电子设备及存储介质 |
| CN116859485A (zh) * | 2023-06-27 | 2023-10-10 | 成都理工大学 | 一种高频层序识别与对比方法、系统、设备及终端 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105182444A (zh) * | 2015-09-09 | 2015-12-23 | 中国海洋石油总公司 | 一种高分辨率层序地层格架约束的地质统计学反演方法 |
| CN106707344A (zh) * | 2016-12-12 | 2017-05-24 | 中国石油天然气股份有限公司 | 一种地层层序划分方法及装置 |
| CN108254786A (zh) * | 2018-01-02 | 2018-07-06 | 中国石油天然气集团公司 | 一种储层层序识别方法及装置 |
-
2019
- 2019-10-11 CN CN201910962403.7A patent/CN112649887B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105182444A (zh) * | 2015-09-09 | 2015-12-23 | 中国海洋石油总公司 | 一种高分辨率层序地层格架约束的地质统计学反演方法 |
| CN106707344A (zh) * | 2016-12-12 | 2017-05-24 | 中国石油天然气股份有限公司 | 一种地层层序划分方法及装置 |
| CN108254786A (zh) * | 2018-01-02 | 2018-07-06 | 中国石油天然气集团公司 | 一种储层层序识别方法及装置 |
Non-Patent Citations (4)
| Title |
|---|
| 利用米兰科维奇旋回划分柴达木盆地第四系层序地层;郭少斌;陈成龙;;地质科技情报(第04期);第27-30页 * |
| 小波变换在饶阳凹陷层序地层格架建立中的应用;张勇 等;中国煤炭地质;第22卷(第6期);摘要，第50-54页，图1-图2 * |
| 川西白马庙气田上侏罗统蓬莱镇组高分辨率层序地层对比;淡卫东;张昌民;尹太举;程绪彬;李连民;;沉积学报(第05期);第708-715页 * |
| 柴达木盆地东部三湖地区四系米兰柯维奇旋回分析;王永军 等;地球物理学进展;第22卷(第2期);第545-546、549页，图1-图3 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112649887A (zh) | 2021-04-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Lu et al. | Seismic spectral decomposition using deconvolutive short-time Fourier transform spectrogram | |
| Liu et al. | Seismic data interpolation beyond aliasing using regularized nonstationary autoregression | |
| CN109709607B (zh) | 一种预测薄层砂岩储层厚度方法及装置 | |
| CN103758511B (zh) | 一种井下逆时偏移成像识别隐蔽储层的方法及装置 | |
| CN112649887B (zh) | 基于钻井资料定量划分层序地层的方法及装置 | |
| CN105954802A (zh) | 一种岩性数据体的转换方法及装置 | |
| US20230058742A1 (en) | Method and storage medium for quantitative reconstruction of paleowater depth based on milankovitch cycles | |
| CN107894613A (zh) | 弹性波矢量成像方法、装置、存储介质及设备 | |
| CN103852788A (zh) | 一种基于复地震道分解和重构的地震相位和频率校正方法 | |
| CN102721979A (zh) | 一种基于地震资料的薄层自动解释及厚度预测方法和装置 | |
| CN103969685B (zh) | 一种薄互层地震信号的处理方法 | |
| RU2490677C2 (ru) | Способ комплексной обработки геофизических данных и технологическая система "литоскан" для его осуществления | |
| Fan et al. | Multisource least-squares reverse-time migration with structure-oriented filtering | |
| CN112213782B (zh) | 分相位地震数据的处理方法、装置和服务器 | |
| CN117251674A (zh) | 一种地震相分类方法、装置、电子设备及介质 | |
| Xie et al. | An improved ant-tracking workflow based on divided-frequency data for fracture detection | |
| Huang et al. | Frequency–Space-Dependent Smoothing Regularized Nonstationary Predictive Filtering | |
| CN107884830A (zh) | 沉积地层旋回划分方法及装置 | |
| CN108318936B (zh) | 一种地层划分处理方法和装置 | |
| CN104345347B (zh) | 一种用于致密含气砂岩储层预测的测井曲线恢复方法 | |
| CN114428298A (zh) | 断溶体分带性的识别方法、装置、电子设备及存储介质 | |
| CN107871028B (zh) | 一种层序地层的级别划分方法及装置 | |
| Kopylova et al. | Statistical analysis of precision water level data from observations in a seismoactive region: Case study of the YuZ-5 well, Kamchatka | |
| Famiani et al. | Detection of hidden faults within the Fucino basin from single-station ambient noise measurements: The case study of the Trasacco fault system | |
| CN111562629A (zh) | 一种基于等效孔隙截面指数的饱和度确定方法及装置 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |